*A1=T33*(l1+2.0*u1)*l1/u1/(l1+u1)/(l1+u1)/(m*m*m)/2.0;
      *B1=0.0;
      *C1=-(l1+2.0*u1)*T33/u1/(l1+u1)/(m*m)/2.0;
      *D1=0.0;
      *A2=T33*(l1+2.0*u1)*l1/u1/(l1+u1)/(l1+u1)/(m*m*m)/2.0;
      *C2=-(l1+2.0*u1)*T33/u1/(l1+u1)/(m*m)/2.0;