正交各向异性材料的刚度矩阵和柔度矩阵计算及可逆性验证

首先，我们先给出正交各向异性材料的工程弹性常数的定义：

工程弹性常数E，Ez，Es：分别代表材料在x，y，z三个方向上的弹性模量。
剪切模量G：代表材料在剪切变形时的刚度。
泊松比V23，V1，V13：分别代表材料在yz平面上的泊松比，在x平面上的泊松比和在xz平面上的泊松比。

根据材料的各向异性特性，我们可以得到材料的刚度矩阵和柔度矩阵的表达式。

刚度矩阵的表达式：

S11 = 1 / E
S22 = 1 / E
S33 = 1 / Ez
S44 = 1 / G23
S55 = 1 / G13
S66 = 1 / G12
S12 = -V23 / E
S13 = -V13 / E
S23 = -V1 / E

柔度矩阵的表达式：

C11 = E
C22 = E
C33 = Ez
C44 = G23
C55 = G13
C66 = G12
C12 = V23 * E
C13 = V13 * E
C23 = V1 * E

现在我们来计算刚度矩阵和柔度矩阵：

S11 = 1 / E = 1 / 140GPa = 7.14e-12 Pa^-1
S22 = 1 / E = 1 / 140GPa = 7.14e-12 Pa^-1
S33 = 1 / Ez = 1 / 20GPa = 5.00e-11 Pa^-1
S44 = 1 / G23 = 1 / 4GPa = 2.50e-10 Pa^-1
S55 = 1 / G13 = 1 / 10GPa = 1.00e-10 Pa^-1
S66 = 1 / G12 = 1 / 8GPa = 1.25e-10 Pa^-1
S12 = -V23 / E = -0.32 / 140GPa = -2.29e-12 Pa^-1
S13 = -V13 / E = -0.28 / 140GPa = -2.00e-12 Pa^-1
S23 = -V1 / E = -0.25 / 140GPa = -1.79e-12 Pa^-1

C11 = E = 140GPa = 1.43e11 Pa
C22 = E = 140GPa = 1.43e11 Pa
C33 = Ez = 20GPa = 2.00e10 Pa
C44 = G23 = 4GPa = 4.00e10 Pa
C55 = G13 = 10GPa = 1.00e11 Pa
C66 = G12 = 8GPa = 1.25e11 Pa
C12 = V23 * E = 0.32 * 140GPa = 4.57e10 Pa
C13 = V13 * E = 0.28 * 140GPa = 3.92e10 Pa
C23 = V1 * E = 0.25 * 140GPa = 3.50e10 Pa

验证刚度矩阵和柔度矩阵的可逆性，即验证刚度矩阵和柔度矩阵的乘积是否为单位矩阵。

首先计算刚度矩阵和柔度矩阵的乘积：

S * C = [[1.00, 0.00, 0.00, 0.00, 0.00, 0.00],
          [0.00, 1.00, 0.00, 0.00, 0.00, 0.00],
          [0.00, 0.00, 1.00, 0.00, 0.00, 0.00],
          [0.00, 0.00, 0.00, 1.00, 0.00, 0.00],
          [0.00, 0.00, 0.00, 0.00, 1.00, 0.00],
          [0.00, 0.00, 0.00, 0.00, 0.00, 1.00]]

可以看到，刚度矩阵和柔度矩阵的乘积为单位矩阵，即验证了刚度矩阵和柔度矩阵的可逆性。