3. 柔性体动力学基础

柔性体仿真的数学根基是连续介质力学。与刚体只用有限个广义坐标描述位姿不同，柔性体在空间中占据连续区域，其质量、位移、应变、应力都是空间点与时间的场量。理解这一节的关键，是把“连续体”看成无穷多质点的集合，并建立描述其变形与受力的统一语言。

连续介质假设与基本场量

连续介质假设认为：物质在宏观尺度上连续分布，不存在空隙，因此可以定义密度 ρ(X,t)、位移 u(X,t)、速度 v(X,t) 等场量，并对它们求导和积分。这里 X 是参考构型（未变形状态）中的物质坐标，x 是当前构型中的空间坐标，映射 x = φ(X,t) 描述物质点的运动。位移场定义为 u = x − X，它把“变形”转化为一个可微的向量场。

在柔性体动力学中，我们通常同时使用两套构型：参考构型用于定义材料属性与应变，当前构型用于计算实际受力与几何非线性。二者之间的桥梁是变形梯度张量 F = ∂x/∂X。F 的雅可比行列式 J = det(F) 表示局部体积变化率，J=1 对应不可压缩材料，J≠1 则体现可压缩性。这一张量是后续所有应变、应力度量的出发点。

小变形与大变形假设

小变形假设要求位移梯度远小于 1，即 ‖∂u/∂X‖ ≪ 1。此时可近似认为参考构型与当前构型重合，应变与位移呈线性关系，平衡方程可建立在未变形构型上，刚度矩阵为常数。这一假设适用于金属结构在弹性极限内的振动、小振幅疲劳分析等场景，计算效率高，但一旦位移达到结构尺寸的百分之几，误差便迅速增大。

大变形假设则保留 F 的完整非线性，应变与位移不再线性，平衡方程必须建立在当前构型上，并区分几何非线性与材料非线性。典型场景包括橡胶密封圈压缩、薄膜充气、软体机器人弯曲、生物软组织受力等。此时若仍用小变形公式，会严重低估或高估刚度，甚至导致仿真发散。选择哪种假设，取决于最大位移与特征长度之比，以及是否出现大转动、大曲率或屈曲。

应变与应力度量

应变描述局部变形的程度。小变形下使用工程应变 ε = ½(∇u + ∇uᵀ)，其分量线性叠加，便于直接组装刚度矩阵。大变形下必须使用有限应变张量，最常用的是 Green-Lagrange 应变 E = ½(FᵀF − I)，它定义在参考构型上，且对刚体转动不产生虚假应变。与之对应的还有 Almansi 应变 e = ½(I − F⁻ᵀF⁻¹)，定义在当前构型上。

应力同样有两套度量。Cauchy 应力 σ 定义在当前构型单位面积上，是真实物理应力，但用于参考构型积分时不便。第一类 Piola-Kirchhoff 应力 P 把当前构型的面力映射到参考构型面积上，P = JσF⁻ᵀ，但它不对称。第二类 Piola-Kirchhoff 应力 S = F⁻¹P 对称且与 Green 应变能量共轭，因此在虚功原理和有限元列式中最常使用。三者之间的转换关系必须与所选构型严格对应，否则能量不守恒。

虚功原理与动力学方程

虚功原理是连接连续介质力学与数值离散的桥梁。它不直接求解强形式的平衡微分方程，而是要求对任意虚位移 δu，内力虚功等于外力虚功：∫ S : δE dV₀ = ∫ ρ₀ b·δu dV₀ + ∫ t·δu dA。这一弱形式降低了对场函数连续性的要求，天然适合有限元离散。

将位移场用形函数插值 u = N q，代入虚功原理并整理，即可得到柔性体的动力学方程 M q̈ + C q̇ + K(q) q = F(t)。小变形时 K 为常数矩阵，可直接模态叠加；大变形时 K 随 q 变化，需在每个时间步更新切线刚度并迭代求解。质量矩阵 M 和阻尼矩阵 C 的构造方式（一致质量或集中质量、Rayleigh 阻尼或比例阻尼）会直接影响仿真稳定性与频率精度。

教学示例：悬臂梁的大变形判断

取一根长度 L = 0.5 m、截面 20 mm × 20 mm 的钢制悬臂梁，弹性模量 E = 210 GPa，密度 ρ = 7850 kg/m³。末端施加横向集中力 F。若 F = 50 N，按小变形理论计算末端挠度约为 0.7 mm，与 L 之比约 0.14%，小变形假设成立，K 可取常数。若 F 增大到 5000 N，小变形公式给出约 70 mm，但实际已进入大变形区间，需用 Green 应变与第二类 PK 应力迭代求解，末端挠度会明显小于线性预测值，且梁轴线出现大转动。判断依据可先估算无量纲挠度 δ/L：小于 1% 可优先用小变形，1%～5% 需谨慎，超过 5% 必须采用大变形列式。

掌握这些基础后，后续的有限元离散、刚柔耦合约束以及时间积分算法才有统一的力学依据。实际建模时，先根据几何与载荷判断变形量级，再选择应变与应力度量，最后用虚功原理组装方程，这条路径能有效避免因假设误用导致的刚度失真与仿真发散。