$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