<?xml version="1.0" encoding="ISO-8859-1"?> <OpenGeoSysProject> <mesh axially_symmetric="true">SimpleAxisymmetricCreepWithAnalyticSolutionMFront.vtu</mesh> <geometry>boundary.gml</geometry> <processes> <process> <name>ThermoMechanics</name> <type>THERMO_MECHANICS</type> <integration_order>2</integration_order> <constitutive_relation> <type>MFront</type> <behaviour>PowerLawLinearCreep</behaviour> <material_properties> <material_property name="YoungModulus" parameter="E"/> <material_property name="PoissonRatio" parameter="nu"/> <material_property name="PowerLawFactor" parameter="A"/> <material_property name="PowerLawEnergy" parameter="Q"/> <material_property name="PowerLawExponent" parameter="n"/> <material_property name="ReferenceStress" parameter="sigma_f"/> <material_property name="LinearLawFactor" parameter="A2"/> <material_property name="LinearLawEnergy" parameter="Q"/> <material_property name="SaltGrainSize" parameter="Dgrain"/> </material_properties> </constitutive_relation> <process_variables> <displacement>displacement</displacement> <temperature>temperature</temperature> </process_variables> <secondary_variables> <secondary_variable internal_name="sigma" output_name="sigma"/> <secondary_variable internal_name="epsilon" output_name="epsilon"/> </secondary_variables> <specific_body_force>0 0</specific_body_force> </process> </processes> <media> <medium> <phases> <phase> <type>Solid</type> <properties> <property> <name>density</name> <type>Constant</type> <value>0.</value> </property> <property> <name>specific_heat_capacity</name> <type>Constant</type> <value>100</value> </property> <property> <name>thermal_conductivity</name> <type>Constant</type> <value>1</value> </property> <property> <name>thermal_expansivity</name> <type>Constant</type> <value>4e-5</value> </property> </properties> </phase> </phases> </medium> </media> <time_loop> <processes> <process ref="ThermoMechanics"> <nonlinear_solver>basic_newton</nonlinear_solver> <convergence_criterion> <type>PerComponentDeltaX</type> <norm_type>NORM2</norm_type> <abstols>5e-12 1e-16 1e-16</abstols> </convergence_criterion> <time_discretization> <type>BackwardEuler</type> </time_discretization> <time_stepping> <type>FixedTimeStepping</type> <t_initial>0</t_initial> <t_end>100</t_end> <timesteps> <pair> <repeat>100</repeat> <delta_t>0.1</delta_t> </pair> </timesteps> </time_stepping> </process> </processes> <output> <type>VTK</type> <prefix>SimpleAxisymmetricCreepWithAnalyticalSolutionMFront</prefix> <timesteps> <pair> <repeat>10</repeat> <each_steps>100</each_steps> </pair> </timesteps> <variables> <variable>displacement</variable> <variable>temperature</variable> <variable>sigma</variable> <variable>epsilon</variable> </variables> <suffix>_ts_{:timestep}_t_{:time}</suffix> </output> </time_loop> <parameters> <parameter> <name>A</name> <type>Constant</type> <value>0.18</value> </parameter> <parameter> <name>n</name> <type>Constant</type> <value>5.0</value> </parameter> <parameter> <name>sigma_f</name> <type>Constant</type> <value>1</value> </parameter> <parameter> <name>Q</name> <type>Constant</type> <value>54000</value> </parameter> <parameter> <name>A2</name> <type>Constant</type> <value>0.0</value> </parameter> <parameter> <name>Dgrain</name> <type>Constant</type> <value>5e-2</value> </parameter> <parameter> <name>E</name> <type>Constant</type> <value>25000</value> </parameter> <parameter> <name>nu</name> <type>Constant</type> <value>0.27</value> </parameter> <parameter> <name>temperature_ic</name> <type>Constant</type> <value>373.15</value> </parameter> <parameter> <name>temperature_bc</name> <type>Constant</type> <value>373.15</value> </parameter> <parameter> <name>displacement0</name> <type>Constant</type> <values>0 0</values> </parameter> <parameter> <name>dirichlet0</name> <type>Constant</type> <value>0</value> </parameter> <parameter> <name>m_pressure</name> <type>Constant</type> <value>5.0</value> </parameter> </parameters> <process_variables> <process_variable> <name>displacement</name> <components>2</components> <order>1</order> <initial_condition>displacement0</initial_condition> <boundary_conditions> <boundary_condition> <geometrical_set>cpreepBGRA</geometrical_set> <geometry>bottom</geometry> <type>Neumann</type> <component>1</component> <parameter>m_pressure</parameter> </boundary_condition> <boundary_condition> <geometrical_set>cpreepBGRA</geometrical_set> <geometry>top</geometry> <type>Dirichlet</type> <component>1</component> <parameter>dirichlet0</parameter> </boundary_condition> </boundary_conditions> </process_variable> <process_variable> <name>temperature</name> <components>1</components> <order>1</order> <initial_condition>temperature_ic</initial_condition> <boundary_conditions> <boundary_condition> <geometrical_set>cpreepBGRA</geometrical_set> <geometry>top</geometry> <type>Dirichlet</type> <parameter>temperature_bc</parameter> </boundary_condition> <boundary_condition> <geometrical_set>cpreepBGRA</geometrical_set> <geometry>bottom</geometry> <type>Dirichlet</type> <parameter>temperature_bc</parameter> </boundary_condition> </boundary_conditions> </process_variable> </process_variables> <nonlinear_solvers> <nonlinear_solver> <name>basic_newton</name> <type>Newton</type> <max_iter>50</max_iter> <linear_solver>general_linear_solver</linear_solver> </nonlinear_solver> </nonlinear_solvers> <linear_solvers> <linear_solver> <name>general_linear_solver</name> <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis> <eigen> <solver_type>BiCGSTAB</solver_type> <precon_type>DIAGONAL</precon_type> <max_iteration_step>10000</max_iteration_step> <error_tolerance>1e-17</error_tolerance> </eigen> </linear_solver> </linear_solvers> </OpenGeoSysProject>