2. 复合材料力学基础与冲击响应

复合材料冲击仿真的可信度，首先取决于力学模型是否与材料真实行为匹配。与各向同性金属不同，碳纤维层合板在面内呈正交各向异性，在厚度方向由基体主导，冲击下同时出现纤维拉伸/压缩、基体开裂、层间分层与面内剪切等多种损伤模式。本节围绕各向异性弹性力学、层合板理论、冲击动力学方程、应变率效应与损伤起始判据五条主线，建立后续建模所需的力学基础。

一、各向异性弹性力学：从广义胡克定律到工程常数

一般各向异性材料的应力—应变关系为 σ_i = C_ijε_j，刚度矩阵 C 最多有 21 个独立常数。单层板（单向带）具有一个弹性对称面，退化为正交各向异性，独立常数减为 9 个。在材料主轴 1-2-3 下（1 为纤维方向，2 为面内垂直纤维方向，3 为厚度方向），柔度矩阵可写成工程常数形式：

ε1 = σ1/E1 - ν21σ2/E2 - ν31σ3/E3
ε2 = -ν12σ1/E1 + σ2/E2 - ν32σ3/E3
ε3 = -ν13σ1/E1 - ν23σ2/E2 + σ3/E3
γ12 = τ12/G12,  γ23 = τ23/G23,  γ31 = τ31/G31

典型 T700/环氧单层板示例参数（教学假设值）：E1=135 GPa，E2=E3=9.5 GPa，ν12=ν13=0.30，ν23=0.45，G12=G13=5.0 GPa，G23=3.3 GPa，ρ=1600 kg/m³。可见 E1/E2≈14，纤维方向刚度远高于横向，这正是冲击下损伤沿纤维与基体呈现明显方向性的根源。仿真中若把层合板简化为各向同性，冲击背面的分层与纤维断裂位置将出现系统性偏差。

二、层合板理论：从单层刚度到整体响应

经典层合板理论（CLT）假设直法线保持且忽略厚度方向正应力，将第 k 层的偏轴刚度矩阵 Q̄^(k) 沿厚度积分，得到面内力—应变关系：

[N]   [A B] [ε⁰]
[ ] = [   ] [  ]
[M]   [B D] [κ ]

其中 A 为拉伸刚度，D 为弯曲刚度，B 为耦合刚度。对称铺层（如 [0/90]_s）时 B=0，拉伸与弯曲解耦；非对称铺层会产生拉弯耦合，冲击下更容易出现翘曲与早期分层。铺层顺序直接决定弯曲刚度分布，例如 [0₂/90₂]_s 与 [90₂/0₂]_s 面内刚度相同，但冲击背面拉应力分布不同，损伤形貌差异明显。因此建模时必须按真实铺层逐层定义厚度、角度与材料方向。

三、冲击动力学基本方程

冲击问题必须计入惯性效应，控制方程为动量守恒：

ρ ∂²u/∂t² = ∇·σ + f

与准静态不同，冲击下应力波在层合板内传播，波速由材料刚度与密度决定。沿纤维方向纵波速 c₁=√(E1/ρ)≈9.2 km/s，横向约 2.4 km/s，厚度方向更低。当冲击持续时间与波在厚度方向往返时间可比时，厚度方向应力波反射会显著影响背面拉应力，进而决定分层与纤维断裂位置。LS-DYNA 显式积分中，时间步长受最小单元尺寸与波速控制：Δt ≤ L_min/c，因此冲击区网格尺寸直接决定计算成本与能否捕捉高频响应。

四、应变率效应

聚合物基体与部分纤维对应变率敏感。冲击下应变率可达 10²~10⁴ s⁻¹，基体屈服强度与模量随率提高，表现为“率强化”。工程上常用 Cowper-Symonds 模型描述：

σ_y = σ₀ [1 + (ε̇/D)^(1/p)]

示例参数（教学假设）：D=8000 s⁻¹，p=4。当 ε̇=1000 s⁻¹ 时，动态屈服约为静态的 1.6 倍。若忽略率效应，仿真会低估冲击初期承载能力，导致损伤起始时间偏早、分层面积偏大。LS-DYNA 中可通过 *MAT_PLASTIC_KINEMATIC 的 C、P 参数或 *MAT_024 的率相关选项引入。

五、损伤起始判据

层合板常用 Hashin 三维判据区分纤维与基体、拉伸与压缩四种模式。以纤维拉伸为例：

F_f^t = (σ1/X_T)² + (τ12/S12)² + (τ13/S13)²  ≥ 1 起始

基体拉伸：F_m^t = (σ2/Y_T)² + (τ12/S12)² + (τ23/S23)² ≥ 1。层间分层则用内聚力模型（CZM），以牵引—分离律描述界面，峰值牵引力对应起始，断裂能 G_Ic、G_IIc 控制软化段。示例参数：X_T=2200 MPa，Y_T=55 MPa，S12=90 MPa，G_Ic=0.3 N/mm，G_IIc=0.8 N/mm。仿真中需将判据与单元特征尺寸匹配，避免网格依赖导致损伤面积随网格细化而虚假增大。

综上，本节建立的各向异性本构、层合板刚度、冲击动力学方程、率相关模型与损伤判据，是后续在 LS-DYNA 中定义材料卡、铺层、接触与失效准则的直接依据。理解这些关系，才能在仿真中正确设置参数并合理解读损伤结果。