檢視 $Thermallattice 的原始碼

$Thermallattice is the command to setup the thermal expansion coefficient of material, which is used with the command [[$solveheatstrain]], [[$heatstrainsteps]]. 

 The command is relative complication, type common format is 
 
 $Thermallattice
 type_1  Number_of_input_parameters_N1 <math>p_{1}</math> <math>p_{2}</math> <math>p_{3}</math> <math>p_{4}</math>  .... <math>p_{N1}</math>
 type_2  Number_of_input_parameters_N2 <math>p_{1}</math> <math>p_{2}</math> <math>p_{3}</math> <math>p_{4}</math>  .... <math>p_{N2}</math>
 type_3  Number_of_input_parameters_N3 <math>p_{1}</math> <math>p_{2}</math> <math>p_{3}</math> <math>p_{4}</math>  .... <math>p_{N3}</math>
  ....
  .... top total layers
 type_m  Number_of_input_parameters_Nm <math>p_{1}</math> <math>p_{2}</math> <math>p_{3}</math> <math>p_{4}</math>  .... <math>p_{Nm}</math>  

 Type:
 1: means  latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (T - <math>p_{1}</math> ) * <math>p_{2}</math> and 
     Number_of_input_parameters =2 
 2: means x x y z axis have different thermal expansion coefficient. So 
    latticeX (T) = latticeX(Tref) + latticeX(Tref) * (T - <math>p_{1}</math> ) * <math>p_{2}</math> 
    latticeY (T) = latticeY(Tref) + latticeY(Tref) * (T - <math>p_{1}</math> ) * <math>p_{3}</math> 
    latticeZ (T) = latticeZ(Tref) + latticeZ(Tref) * (T - <math>p_{1}</math> ) * <math>p_{4}</math>
    Number_of_input_parameters = 4 
 3: If T> <math>p_{3}</math> 
    latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (p_{3} - <math>p_{1}</math> ) * <math>p_{2}</math>  + latticeXYZ(Tref) * (T - <math>p_{3}</math> ) * <math>p_{4}</math> 
    If T<=<math>p_{3}</math> 
    latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (T - <math>p_{1}</math> ) * <math>p_{2}</math>  
    Number_of_input_parameters = 4 

For example, If we have 6 regions 
 $Thermallattice
 1 2  300.0 0.01 
 1 2  300.0 0.03
 2 4  300.0 0.01 0.015 0.012  
 3 4  300.0 0.01  400.0 0.1 
 1 2  300.0 0.01 
 1 2  300.0 0.03