Differences

This shows you the differences between two versions of the page.

--- doc:user:integration:scheme:staggeredtm [2013/07/12 15:24] – created joris
+++ doc:user:integration:scheme:staggeredtm [2016/03/30 15:23] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
+====== Thermomechanical staggered integration scheme ======
+===== Description =====
+==== Old Metafor Version <= 2422 ====
+Used if ''MDE_NDYT=2''. The mechanical scheme must be quasi-static  (''MDE_NDYN=0'') or alpha genealized (''MDE_NDYN=2''). The thermal scheme is trapezoidal (''MDE_THTY=1'') or generalized midpoint (''MDE_THTY=0''). (see [[doc:user:integration:general:parameters]]).
+For staggered schemes can be used (depending on ''MDE_NDYT''):
+  * ''MDE_NDYT = 2'' : thermomechanical staged isothermal without reevaluation of stresses at the end of a step
+  * ''MDE_NDYT = 3'' : thermomechanical staged adiabatic without reevaluation of stresses at the end of a step
+  * ''MDE_NDYT = 4'' : thermomechanical staged isothermal with reevaluation of stresses at the end of a step
+  * ''MDE_NDYT = 5'' : thermomechanical staged adiabatic with reevaluation of stresses at the end of a step
+==== New Metafor Version > 2422 ====
+<code>
+tiMech = AlphaGeneralizedTimeIntegration(metafor)
+tiTher = TrapezoidalThermalTimeIntegration(metafor)
+tiTher.setTheta(_theta)
+ti = StaggeredTmTimeIntegration(metafor)
+ti.setIsAdiabatic(False)
+ti.setWithStressReevaluation(False)
+ti.setMechanicalTimeIntegration(tiMech)
+ti.setThermalTimeIntegration(tiTher)
+metafor.setTimeIntegration(ti)
+</code>
+The mechanical time integration scheme is either quasi static or alpha generalized family. The default _theta  parameter value is 1.0 ([0,1]).
+===== Data set =====
+To change a test case running with quasiStaticAlgorithm into a thermomechanical test :
+  * Replace elements by their thermal version : ''Volume2DElement'' becomes ''TmVolume2DElement''.
+  * Replace materials by their thermal version : ''EvpIsoHHypoMaterial'' becomes ''TmEvpIsoHHypoMaterial''.
+  * Add new material dependency (optional) : ''matset[1].depend(P_HARD_1, fonction, Key(TO|RE))''.
+  * Add thermal parameters of material law (see materials).
+  * Add initial temperature (optional) (''metafor.getInitialConditionSet().define(no,type,TO|AB,val)'').
+  * Add temperature to contact matrices (optional) (see [[doc:user:elements:general:start|ElementProperties]]).
+  * Add thermal boundary conditions:
+    * ''ElementProperties'' and ''ConvInteraction'' for [[doc:user:elements:boundaries:convection|convection]].
+    * ''ElementProperties'' and ''FluxInteraction'' for [[doc:user:elements:boundaries:flux|fixed thermal flux]].
+    * [[doc:user:conditions:temperature]].
+===== Additional parameters =====
+  * ''MDE_THTY'' :
+    * = 0 generalized midpoint rule for the integration of transient thermal part (set ''MDR_THET'').
+    * = 1 generalized trapezoid rule for the integration of transient thermal part (set ''MDR_THET'').
+  * ''MDE_NDYT'' : must be equal to 2, 3, 4 or 5 to run the staged scheme.
+(see [[doc:user:integration:general:parameters]])