$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 $p_{1}$ $p_{2}$ $p_{3}$ $p_{4}$  .... $p_{N1}$
type_2  Number_of_input_parameters_N2 $p_{1}$ $p_{2}$ $p_{3}$ $p_{4}$  .... $p_{N2}$
type_3  Number_of_input_parameters_N3 $p_{1}$ $p_{2}$ $p_{3}$ $p_{4}$  .... $p_{N3}$
 ....
 .... top total layers
type_m  Number_of_input_parameters_Nm $p_{1}$ $p_{2}$ $p_{3}$ $p_{4}$  .... $p_{Nm}$  

Type:
1: means  latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (T - $p_{1}$ ) * $p_{2}$ 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 - $p_{1}$ ) * $p_{2}$ 
   latticeY (T) = latticeY(Tref) + latticeY(Tref) * (T - $p_{1}$ ) * $p_{3}$ 
   latticeZ (T) = latticeZ(Tref) + latticeZ(Tref) * (T - $p_{1}$ ) * $p_{4}$
   Number_of_input_parameters = 4 
3: If T> $p_{3}$ 
   latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (p_{3} - $p_{1}$ ) * $p_{2}$  + latticeXYZ(Tref) * (T - $p_{3}$ ) * $p_{4}$ 
   If T<=$p_{3}$ 
   latticeXYZ (T) = latticeXYZ(Tref) + latticeXYZ(Tref) * (T - $p_{1}$ ) * $p_{2}$  
   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