固体力学问题的数学描述基础
固体力学研究固体材料在外部载荷作用下的变形、应力分布及稳定性问题。为了从理论上精确预测材料的力学行为，必须建立一套严谨的数学描述体系。这一体系主要由几何方程、物理方程（本构关系）和平衡方程（或运动方程）三部分组成，辅以边界条件，共同构成了求解固体力学问题的完整框架。
几何描述与变形分析
几何描述主要关注物体在受力前后的形状变化。在连续介质力学中，通常采用拉格朗日描述（物质描述）或欧拉描述（空间描述）。对于大变形问题，区分参考构型（Reference Configuration）和当前构型（Current Configuration）至关重要。
变形梯度张量 $\mathbf{F}$ 是描述局部变形的核心量，定义为位置矢量的偏导数：
$$ \mathbf{F} = \frac{\partial \mathbf{x}}{\partial \mathbf{X}} $$
其中 $\mathbf{X}$ 为参考构型中的坐标，$\mathbf{x}$ 为当前构型中的坐标。基于变形梯度，可以进一步定义应变张量。对于小变形假设，通常使用柯西应变张量（Cauchy Strain Tensor）$\boldsymbol{\varepsilon}$，其分量定义为位移场 $\mathbf{u}$ 的对称梯度：
$$ \varepsilon_{ij} = \frac{1}{2} \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) $$
对于大变形情况，则需使用格林-拉格朗日应变张量或阿尔曼西应变张量等有限应变度量，以准确捕捉几何非线性效应。
物理方程与本构关系
物理方程建立了应力与应变之间的定量联系，即本构关系（Constitutive Relations）。这是体现材料特性的关键环节，不同的材料模型对应不同的本构方程。
线弹性材料：遵循胡克定律。对于各向同性材料，广义胡克定律将应力张量 $\boldsymbol{\sigma}$ 与应变张量 $\boldsymbol{\varepsilon}$ 线性关联：
$$ \sigma_{ij} = \lambda \varepsilon_{kk} \delta_{ij} + 2\mu \varepsilon_{ij} $$
其中 $\lambda$ 和 $\mu$ 为拉梅常数，$\delta_{ij}$ 为克罗内克符号。非线性材料：包括超弹性材料、弹塑性材料等。例如，J2 塑性理论通过引入屈服函数和流动法则，描述了材料在超过屈服极限后的塑性流动行为。
粘弹性材料：应力不仅取决于当前应变，还取决于应变的历史，通常通过积分型或微分型本构方程描述，如 Maxwell 模型或 Kelvin-Voigt 模型。
平衡方程与运动方程
平衡方程基于牛顿第二定律，描述了物体内部力的平衡状态。在静态问题中，体力 $\mathbf{b}$ 和表面力 $\mathbf{t}$ 需满足以下平衡条件：
$$ \nabla \cdot \boldsymbol{\sigma} + \mathbf{b} = \mathbf{0} $$
在动态问题中，需考虑惯性力，方程变为：
$$ \nabla \cdot \boldsymbol{\sigma} + \mathbf{b} = \rho \mathbf{a} $$
其中 $\rho$ 为材料密度，$\mathbf{a}$ 为加速度。该方程确保了物体在任意微元体上的力矩和合力平衡，是求解位移场或应力场的基础控制方程。
边界条件与求解策略
一个完整的固体力学问题必须包含适当的边界条件，以保证解的唯一性。边界条件通常分为三类：
- 位移边界条件（Dirichlet 条件）：在部分边界 $\Gamma_u$ 上，指定位移值 $\mathbf{u} = \bar{\mathbf{u}}$。
- 力边界条件（Neumann 条件）：在部分边界 $\Gamma_t$ 上，指定面力 $\mathbf{t} = \bar{\mathbf{t}}$，且满足 $\mathbf{t} = \boldsymbol{\sigma} \cdot \mathbf{n}$，其中 $\mathbf{n}$ 为边界外法向量。
- 混合边界条件：在复杂结构中，边界可能同时包含位移和力约束，或者存在接触边界条件。
在实际工程应用中，解析解往往难以获得，因此常采用数值方法（如有限元法 FEM）将上述偏微分方程组离散化。通过引入变分原理（如最小势能原理），可以将强形式转化为弱形式，从而构建线性或非线性代数方程组进行求解。
总结
固体力学问题的数学描述是一个多层次的系统。从几何变形到物理响应，再到力学平衡，每一步都依赖于精确的数学定义和假设。正确选择本构模型、准确施加边界条件，是确保计算结果物理意义正确且工程应用可靠的前提。掌握这一基础框架，是深入理解复杂结构响应及开发先进数值算法的关键。