Thermomechanical staggered integration scheme

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 Global Parameters [REMOVED]).

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

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)

The mechanical time integration scheme is either quasi static or alpha generalized family. The default _theta parameter value is 1.0 ([0,1]).

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 ElementProperties).
• Add thermal boundary conditions:
  • ElementProperties and ConvInteraction for convection.
  • ElementProperties and FluxInteraction for fixed thermal flux.
  • Prescribed Temperature.

• 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 Global Parameters [REMOVED])