檢視 Elementelasticconstant(6,6,Nel) 的原始碼

Elementelasticconstant(6,6,Nel) are the elastic constants of each element. The strain tensor is in the form of <br>
<math> \Sigma_{i,j=1}^{6} c_{i,j} </math>  <br>
Nel is the element number. These variables save the elastic constants of each elements.   <br>

For example: For wurtzite lattice <br>

 elementelasticconstant(1:6,1:6,totalelements) =0.0
 do i=1,totalelements
   elementelasticconstant(1,1,i) = c11
   elementelasticconstant(1,2,i) = c12
   elementelasticconstant(1,3,i) = c13
   elementelasticconstant(2,2,i) = c11
   elementelasticconstant(3,3,i) = c33
   elementelasticconstant(2,1,i) = c12
   elementelasticconstant(2,3,i) = c13
   elementelasticconstant(3,1,i) = c13
   elementelasticconstant(3,2,i) = c13
   elementelasticconstant(3,3,i) = c33  
   elementelasticconstant(4,4,i) = c44
   elementelasticconstant(5,5,i) = c44
   elementelasticconstant(6,6,i) = c66
   ...
   ...
   ...
 end do