固体力学的定义与研究范畴
固体力学(Solid Mechanics)是研究固体材料在外力、温度、时间等作用下的变形、应力分布以及破坏规律的学科。它以 连续介质假设 为基础,把宏观材料视为连续体,忽略原子尺度的离散性,从而能够用偏微分方程描述材料的力学行为。固体力学的核心目标是:
- 预测:在给定的载荷和约束条件下,材料或结构的位移、应变和应力分布。
- 解释:揭示材料内部微观结构(如晶格、缺陷)与宏观力学响应之间的关系。
- 设计:为工程结构(桥梁、航空器、微电子器件等)提供安全可靠的设计依据。
2. 基本概念与基本假设
| 概念 | 含义 | 常用符号 |
|---|---|---|
| 位移 (Displacement) | 质点从原始位置到新位置的矢量变化 | u |
| 应变 (Strain) | 位移梯度的对称部分,描述形变程度 | ε |
| 应力 (Stress) | 单位面积上作用的内力,描述内部受力状态 | σ |
| 弹性模量 (Elastic Modulus) | 描述材料弹性刚度的系数 | E, G, K |
| 泊松比 (Poisson’s Ratio) | 横向应变与轴向应变的比值 | ν |
基本假设
- 连续介质假设:材料在宏观尺度上连续,无孔隙、裂纹等离散缺陷。
- 小变形假设(线性理论):位移梯度足够小,使得几何非线性可以忽略。
- 材料各向同性或各向异性:材料的力学性质在空间方向上是否相同。
示例:在桥梁梁段的受力分析中,通常采用小变形、各向同性的弹性假设,以简化计算并获得足够精确的结果。
3. 研究范畴
固体力学的研究范围广泛,可划分为以下几个主要方向:
3.1 弹性力学
- 线弹性:遵循胡克定律,σ = C : ε(C 为弹性刚度张量)。
- 各向同性弹性:仅需两个独立常数(如 E 与 ν)即可描述。
- 各向异性弹性:需要 21 个独立常数(如单晶体的弹性常数矩阵)。
案例:金属板在均匀拉伸下的应力分布,可通过平面应力模型求解,得到 σ_x = E/(1-ν²)·ε_x 等式。
3.2 塑性力学
- 屈服准则:如 von Mises、Tresca 判据,用于确定材料何时进入塑性阶段。
- 硬化模型:线性硬化、指数硬化、各向同性/各向异性硬化等。
- 数值实现:常用 返回映射算法(Return Mapping)在有限元软件中实现塑性积分。
示例:在汽车碰撞仿真中,钢板的塑性行为通过 von Mises 屈服准则与线性硬化模型描述,以捕捉局部凹陷和能量吸收。
3.3 粘弹性与粘塑性
- 粘弹性模型:如 Maxwell、Kelvin‑Voigt、标准线性固体模型,用于描述材料的时效响应。
- 粘塑性模型:结合粘弹性与塑性特性,适用于土壤、沥青等材料的长期加载。
案例:道路沥青层在日夜温度循环下的应力松弛,可用标准线性固体模型预测其压实度变化。
3.4 断裂力学
- 能量释放率 (G) 与 应力强度因子 (K):评价裂纹扩展的驱动力。
- 线性弹性断裂力学 (LEFM):适用于脆性材料或小裂纹。
- 弹塑性断裂力学:考虑塑性区对裂纹扩展的阻碍作用。
示例:航空铝合金机翼的疲劳裂纹增长率常用 Paris 法则 da/dN = C·(ΔK)^m 进行预测。
3.5 多场耦合
- 热‑力耦合:温度场影响材料的弹性模量、热膨胀系数,引起热应力。
- 电‑力耦合(压电材料):电场与机械场相互作用。
- 磁‑力耦合(磁致伸缩材料):磁场引起形变。
案例:在微机电系统(MEMS)中,热‑力耦合导致的热膨胀误差需要通过耦合有限元分析进行校正。
4. 基本方程
固体力学的核心方程组包括:
平衡方程(力的平衡)
[
\nabla\cdot \boldsymbol{\sigma} + \mathbf{b} = \rho ,\ddot{\mathbf{u}}
]
其中 b 为体力密度,ρ 为密度,(\ddot{\mathbf{u}}) 为加速度。几何方程(位移‑应变关系)
[
\boldsymbol{\varepsilon} = \frac{1}{2}\left(\nabla\mathbf{u} + (\nabla\mathbf{u})^{\mathrm{T}}\right)
]本构方程(应力‑应变关系)
[
\boldsymbol{\sigma} = \mathsf{C}:\boldsymbol{\varepsilon} \quad (\text{线弹性})
]
对于塑性、粘弹性等情况,本构方程会加入内部变量和历史依赖。边界条件
- 位移边界:(\mathbf{u} = \bar{\mathbf{u}}) 在 (\Gamma_u) 上。
- 力边界:(\boldsymbol{\sigma}\cdot\mathbf{n} = \bar{\mathbf{t}}) 在 (\Gamma_t) 上。
数值实现:上述方程在有限元法(FEM)中通过离散化转化为线性或非线性代数方程组,常见的求解流程包括网格划分、单元刚度矩阵组装、全局方程求解以及后处理。
5. 常用分析方法
- 解析法:适用于几何、载荷、材料均较为简化的情形,如梁、板、圆柱的闭式解。
- 半解析法(如 变分法、积分变换):用于处理边界条件复杂但仍可分离变量的问题。
- 数值法:
- 有限元法(FEM):最为通用,可处理任意几何、非线性材料、复杂边界。
- 边界元法(BEM):适合无限或半无限域的弹性问题。
- 谱方法、网格无关法(Meshless):在高精度或大变形场景中有优势。
6. 实际工程案例
6.1 桥梁主梁的弹性分析
- 模型:采用 Euler‑Bernoulli 梁理论,假设跨中受均布荷载 (q)。
- 关键方程:
[
\frac{\mathrm{d}^2}{\mathrm{d}x^2}\left(EI\frac{\mathrm{d}^2 w}{\mathrm{d}x^2}\right)=q
] - 结果:最大挠度 (w_{\max}= \frac{5qL^4}{384EI}),最大弯矩 (M_{\max}= \frac{qL^2}{8})。
6.2 航空结构的塑性极限分析
- 目标:确定在极限载荷下结构是否出现塑性屈服。
- 方法:使用 极限平衡法(Limit Analysis)结合 有限元,构建塑性势能函数并求解最小上界。
6.3 微电子封装的热‑力耦合
- 情境:芯片功耗产生的热流导致封装材料膨胀,引起焊点应力。
- 步骤:
- 先进行热传导分析得到温度场 (\theta(\mathbf{x}))。
- 将温度场代入本构方程 (\sigma = C:(\varepsilon - \alpha\theta I)),其中 (\alpha) 为热膨胀系数。
- 计算焊点的等效应力,评估是否超过材料的屈服强度。
7. 小结
固体力学通过 平衡方程、几何方程和本构方程 的耦合,系统地描述固体材料在各种外部作用下的响应。其研究范畴涵盖 弹性、塑性、粘弹性、断裂、多场耦合 等多个层面,既有理论解析的深度,也有数值模拟的广度。掌握这些基本概念与方法,是进行结构设计、材料选型以及可靠性评估的前提。
进一步阅读:
- 《弹性力学》—— S. Timoshenko
- 《塑性力学》—— J. Lubliner
- 《断裂力学》—— A. A. Griffith
通过系统学习和实践应用,读者能够在工程实际中准确预测材料行为,提升设计安全性与经济性。