#====================================================================
# Example: Consolidation of a poroelastic sphere subject to impulse
# pressure load on its outer surface boundary
#--------------------------------------------------------------------------------------------------------------------
# Problem Description:
# Coupled solution using QDC_solid and QDC_pressure
#--------------------------------------------------------------------------------------------------------------------
# number_of_nodal_points = 73
# number_of_spatial_dimensions = 2
# number_of_degrees_of_freedom = 3
# number_of_element_groups = 2
# number_of_elements = 60
# number_of_surface_load_cases = 1
# number_of_time_steps = 20
#====================================================================
DEFINE_PROBLEM \
name = "sphere_stagger" \
title = " ----- Consolidation of a poroelastic sphere ------ "
#--------------------------------------------------------------------------------------------------------------------
NODAL_COORDINATES \
generation_type = cylindrical \
reference_coordinates x0=0 y0=0 z0=0
1 0 0.00 0.00
2 4 0.50 0.00 0.00
3 0 4.00 00.0
4 0 4.00 90.0
5 0 0.50 90.0
7 1 8 8
SURFACE_LOADS \
number_of_loads = 1 \
1 1.00 1.00
1 1 9 17 1
8 1 8
NODAL_BOUNDARY_CONDITIONS \
variable = solid_displacement \
generation_type = type_1
1 0 1 1
2 1 0 1
9 0 0 1
66 1 1 0
73 0 1 0
NODAL_BOUNDARY_CONDITIONS \
variable = pressure \
generation_type = type_1
9 8 1
73 0 1
DEFINE_MATERIAL_MODEL \
name = "dummy" \
number_of_material_sets = 1
Stress_Model \
Material_type = Linear \
Material_name = Linear_elastic
Material_Set_Number = 1 \
Youngs_modulus = 1.0e4 \
Poissons_ratio = 0.25 \
solid_mass_density = 2.0 \
fluid_mass_density = 1.0 \
Porosity = 0.30
Scalar_Diffusion_Model \
Material_type = Linear \
Material_name = Scalar_diffusion
Material_Set_Number = 1 \
Mass_density = 1.0 \
Compressibility = 0.0e-6 \
Permeability \
type = isotropic \
k_11 = 1.60e-3
DEFINE_REGION \
name = "Porous_solid" \
element_type = continuum \
element_shape = four_node_quad \
analysis_type = axisymmetric \
strain_displacement = bbar_mean \
number_of_output_sets = 1 \
material_model = "dummy"
1 1 2 3 11 10 1
7 1 1 8 7 8
57 1 1 2 10 18
58 1 1 18 26 34
59 1 1 34 42 50
60 1 1 50 58 66
Field_Output
58 0 17
Body_Force b_x1=0.00 b_x2=0.00 h=0.00
ELEMENT_GROUP \
name = "Solid_equation" \
element_name = QDC_solid
ELEMENT_GROUP \
name = "Diffusion_equation" \
element_name = QDC_pressure
END
TRANSLATOR \
output_format = "femgv"
MESH_DUMP /
nodal_step = 0 /
nodal_frequency = 5 /
spatial_step = 0 /
spatial_frequency = 5 /
final_time = on
NODAL_HISTORY
1 0 3
35 0 3
73 0 2
TIME_SEQUENCE \
number_of_time_steps = 20 \
time_step = 0.025e+0
INITIALIZE_V0
DEFINE_STAGGER name="Group1" \
include_element_group(s) = "Solid_equation" \
equation_type= parabolic \
alpha=1.0 beta=0.0 \
nonlinear_iterations \
jacobian_matrix = on \
max_number_of_iterations = 0 \
variable= Solid_displacement
DEFINE_STAGGER name="Group2" \
include_element_group(s) = "Diffusion_equation" \
equation_type= parabolic \
alpha=1.0 beta=0.0 \
nonlinear_iterations \
jacobian_matrix = on \
max_number_of_iterations = 0 \
variable= Pressure
STAGGER_CONTROL \
max_number_of_iterations = 3 \
convergence_check = on \
convergence_tol_sol = 1.0e-4 \
convergence_tol_cg = 1.0e-9
c PRINT_map
RUN_SOLVER \
mode = execution
STOP
#-------------------------------------------------------------------------------------------------------------------
Notes . .
Notes . .