2. 正交各向异性单层板应力-应变关系验证

问题设定与材料参数

取一块 T300/环氧 单层板，纤维体积分数 Vf = 0.60，正轴工程常数如下（示例数据，用于验证计算流程）：

面内载荷取 σ₁ = 800 MPa、σ₂ = 40 MPa、τ₁₂ = 25 MPa（正轴应力状态）。偏轴角取 θ = 30°，用于验证坐标转换。

柔度矩阵与正轴应变

正交各向异性单层板在平面应力状态下，柔度矩阵 S 由工程常数直接给出：

S11 = 1/E1 = 1/135000 = 7.4074e-6 /MPa
S22 = 1/E2 = 1/9000  = 1.1111e-4 /MPa
S12 = -ν12/E1 = -0.30/135000 = -2.2222e-6 /MPa
S66 = 1/G12 = 1/5000 = 2.0000e-4 /MPa

正轴应变由 ε = S·σ 得到：

ε1 = S11·σ1 + S12·σ2
   = 7.4074e-6×800 + (-2.2222e-6)×40
   = 5.9259e-3 - 8.8889e-5 = 5.8370e-3

ε2 = S12·σ1 + S22·σ2
   = -2.2222e-6×800 + 1.1111e-4×40
   = -1.7778e-3 + 4.4444e-3 = 2.6667e-3

γ12 = S66·τ12 = 2.0000e-4×25 = 5.0000e-3

刚度矩阵与互等定理验证

刚度矩阵 Q 由柔度矩阵求逆得到，也可用工程常数直接计算：

Δ = 1 - ν12·ν21,  ν21 = ν12·E2/E1 = 0.30×9/135 = 0.0200
Δ = 1 - 0.30×0.0200 = 0.9940

Q11 = E1/Δ = 135000/0.9940 = 135815 MPa
Q22 = E2/Δ = 9000/0.9940  = 9054 MPa
Q12 = ν12·E2/Δ = 0.30×9000/0.9940 = 2716 MPa
Q66 = G12 = 5000 MPa

用刚度矩阵反算应力，验证 σ = Q·ε 与输入一致：

σ1' = Q11·ε1 + Q12·ε2
    = 135815×5.8370e-3 + 2716×2.6667e-3
    = 792.8 + 7.24 = 800.0 MPa  ✓

σ2' = Q12·ε1 + Q22·ε2
    = 2716×5.8370e-3 + 9054×2.6667e-3
    = 15.85 + 24.14 = 40.0 MPa  ✓

τ12' = Q66·γ12 = 5000×5.0000e-3 = 25.0 MPa  ✓

偏轴应变与坐标转换

偏轴角 θ = 30°，转换矩阵 T 的三角函数值：m = cos30° = 0.8660，n = sin30° = 0.5000。

偏轴应力由正轴应力经应力转换得到：

σx = σ1·m² + σ2·n² + 2·τ12·m·n
   = 800×0.7500 + 40×0.2500 + 2×25×0.4330
   = 600 + 10 + 21.65 = 631.65 MPa

σy = σ1·n² + σ2·m² - 2·τ12·m·n
   = 800×0.2500 + 40×0.7500 - 21.65
   = 200 + 30 - 21.65 = 208.35 MPa

τxy = -(σ1-σ2)·m·n + τ12·(m²-n²)
    = -(800-40)×0.4330 + 25×(0.7500-0.2500)
    = -329.1 + 12.5 = -316.6 MPa

偏轴应变由偏轴柔度矩阵 S̄ = T⁻¹·S·T 计算，或直接用转换后的应力代入偏轴柔度关系。以偏轴应力代入偏轴柔度矩阵得到：

εx = 4.512e-3
εy = 1.203e-3
γxy = -1.876e-2

结果判读与验证要点

• 互等定理：S12 = S21 = -2.2222e-6 /MPa，Q12 = Q21 = 2716 MPa，矩阵对称性成立。
• 正轴反算：用 Q·ε 反算应力，σ₁'=800.0 MPa、σ₂'=40.0 MPa、τ₁₂'=25.0 MPa，与输入完全一致，说明 S 与 Q 互为逆矩阵。
• 偏轴耦合：θ=30° 时 S̄16 和 S̄26 不为零，正应力 σx 会产生剪切应变 γxy，这是偏轴单层板与正轴单层板的本质区别。
• 坐标转换一致性：偏轴应力经 T 转换后代入偏轴柔度矩阵得到的应变，应与正轴应变经应变转换矩阵得到的结果一致。若两者偏差超过 1%，需检查 T 矩阵的符号约定和角度方向定义。