有没有人正交各向异性材料三维单元刚度矩阵代码？ 100

fact = 1.0-V12*V21-V23*V32-V31*V13-2.0*V21*V32*V13 factor = 1.0 / fact C(1,1) = E11qs * ( 1.0 - V23*V32 ) * factor C(1,2) = E22qs * ( V12 + V32*V13 ) * factor C(1,3) = E33qs * ( V13 + V12*V23 ) * factor C(2,1) = C(1,2) C(2,2) = E22qs * ( 1.0 - V13*V31 ) * factor C(2,3) = E33qs * ( V23 + V21*V13 ) * factor C(3,1) = C(1,3) C(3,2) = C(2,3) C(3,3) = E33qs * ( 1.0 - V12*V21 ) * factor C(4,4) = G12qs C(5,5) = G23qs C(6,6) = G13qs