本页目录

力学 V · 能量法与弹性力学入门

层次：本科核心 + 研究生入口｜ 这一页把前四页统一起来，并搭起通往有限元的桥。前面用的是"针对每种载荷推一个公式"的方法。能量法换一个视角：从系统的应变能出发，一举得到变形、超静定解、以及冲击响应。更重要的是——它是有限元法的理论基础（第二十一页）。

学习层：让一座小桁架同时通过能量、单位载荷和平衡三本账

具体工程情境：检查起吊夹具的 V 形小桁架

把一个起吊夹具的两根斜撑简化为平面 V 形桁架：左右支点相距 \(b=0.80\,\mathrm{m}\)，顶点高 \(h=0.60\,\mathrm{m}\)，两根杆的面积均为 \(A=300\,\mathrm{mm^2}\)，材料为 \(E=210\,\mathrm{GPa}\)。顶点承受向下 \(P=12.0\,\mathrm{kN}\) 和小的水平侧载 \(H=0.80\,\mathrm{kN}\)。工程师要给出顶点位移、杆力和储能，还要证明算出来的刚度确实是一个可逆、稳定的线弹性模型。

必须提交的预测

在几何和材料不变时，\(P\) 加倍，应变能 \(U\) 会变为约 2 倍、4 倍，还是不变？
顶点向下的位移，用 Castigliano 对 \(P\) 求导和用顶点向下的单位载荷法计算，结果应相等还是互为相反数？
对称 V 桁架的 \(2\times2\) 自由节点刚度矩阵，哪一组条件能给出正定证书：只看 \(K_{12}=0\)，还是 \(K_{11}>0\) 且 \(\det K>0\)？

正式公式桥：应变能、单位载荷与平衡复算

取顶点自由位移 \(\mathbf u=[u_x,u_y]^T\)，外载向量为 \(\mathbf f=[H,-P]^T\)。对每一根从支点指向顶点的杆，单位方向为 \(\mathbf n_e\)、长度为 \(L_e\)，则

\[ k_e=\frac{EA}{L_e},\qquad \mathbf K=\sum_e k_e\mathbf n_e\mathbf n_e^T,\qquad \mathbf K\mathbf u=\mathbf f. \]

杆伸长、轴力和应变能分别是

\[ \delta_e=\mathbf n_e^T\mathbf u,\qquad N_e=k_e\delta_e,\qquad U=\sum_e\frac{N_e^2L_e}{2EA}=\frac12\mathbf f^T\mathbf u. \]

把 \(P\) 视为“向下的正载荷”，Castigliano 给出

\[ \frac{\partial U}{\partial P}=\mathbf u^T[0,-1]^T=-u_y. \]

若 \(\bar{\mathbf u}\) 由顶点向下的单位力 \(\bar{\mathbf f}=[0,-1]^T\) 产生，单位载荷法给出同一个向下位移

\[ \delta_P=\sum_e\frac{N_e\bar N_eL_e}{EA}=\mathbf f^T\mathbf K^{-1}\bar{\mathbf f}. \]

实验台还把每根杆的力投影回顶点，检查 \(\mathbf r_f=\mathbf f-\sum_eN_e\mathbf n_e\)；再把反力总和与外载相加，验证整体平衡。对于对称的 \(2\times2\) 矩阵，Sylvester 判据 \(K_{11}>0\) 且 \(\det K>0\) 是正定证书；它比“求出了一个数值位移”更强，因为它保证 \(U=\tfrac12\mathbf u^T\mathbf K\mathbf u>0\)（除零位移外）。

误区与模型边界

\(U\) 对载荷是二次函数，所以载荷加倍时能量加四倍、位移只加两倍；不能把能量账当成力账。
单位载荷法的虚载荷必须沿要找的位移方向；若把向下的 \(P\) 改成向上的广义力，导数符号也会随约定改变。
正定只证明当前支撑和几何下的线弹性自由度没有机构模态；它不证明杆件不屈曲、不屈服，也不覆盖铰接间隙、接触脱开和大变形。
这里的桁架杆只承受轴力，节点是理想铰；真实夹具若有刚性节点、偏心连接或弯曲刚度，必须把梁/框架自由度加入模型。
平衡残差为零是必要条件，不是材料、边界和几何都正确的充分条件；正定证书、量纲检查和独立的单位载荷复算要一起看。

对应 lab 的静态 fallback

无 JavaScript 时仍可核对：默认 \(b=0.80\,\mathrm{m}\)、\(h=0.60\,\mathrm{m}\)、\(A=300\,\mathrm{mm^2}\)、\(E=210\,\mathrm{GPa}\)、\(P=12.0\,\mathrm{kN}\)（向下）、\(H=0.80\,\mathrm{kN}\)（向右）。

账本项	公式/默认读数	单位或判定
斜杆长度 \(L_e\)	\(\sqrt{(b/2)^2+h^2}=0.721110\)	m
\(K_{11},K_{12},K_{22}\)	\(5.3763\times10^7,\ 0,\ 1.2097\times10^8\)	N/m（对称几何下 \(K_{12}=0\)）
顶点位移 \((u_x,u_y)\)	\((0.01488,-0.09920)\)	mm
杆轴力 \((N_L,N_R)\)	\((-7.932,-6.490)\)	kN；压缩为负
应变能 \(U\)	\(0.6012\)	J
Castigliano / 单位载荷	\(-u_y=0.09920\)	mm；对应向下 \(P\) 的位移
顶点平衡残差	\(\\|\mathbf r_f\\|\approx0\)	N，浮点容差内
左 / 右支座反力	\((-4.400,6.600)/(3.600,5.400)\)	kN；由两根杆轴力回算，\(x\) 右、\(y\) 上
整体平衡残差	\(\mathbf f+\mathbf R_L+\mathbf R_R\approx\mathbf0\)	N，浮点容差内
正定证书	\(K_{11}>0,\ \det K>0\)	通过

脚本可用时，移动水平载荷会让两根杆力不再相等；改变高度会同时改变杆长、刚度和正定矩阵的特征值。结果揭示后调参不会重新锁回预测门。

迁移问题

如果右支点改成一根柔性拉杆，顶点多出一个水平约束未知量，你会如何选择冗余力并写出协调条件？请说明怎样用 \(\partial U/\partial X=0\) 和一个独立的反力平衡检查避免把“能量极小”误当成“边界画对了”。

一、应变能

弹性体在均匀、小变形、线弹性本构下变形时储存的能量，单位体积应变能可写为：

\[u = \frac{1}{2}\sigma\varepsilon = \frac{\sigma^2}{2E}\]

各种基本变形的总应变能：

\[U_{\text{拉压}} = \frac{N^2 L}{2EA}, \qquad U_{\text{弯曲}} = \int\frac{M^2}{2EI}dx, \qquad U_{\text{扭转}} = \int\frac{T^2}{2GI_p}dx\]

共同形式：在均匀、线弹性杆件或截面内力模型成立时，\(U = \dfrac{(\text{内力})^2 \times \text{长度}}{2 \times \text{刚度}}\)；变截面、塑性、接触和大变形需要回到积分或非线性本构。

二、卡氏定理：求变形的利器

Castigliano 第二定理：

\[\delta_i = \frac{\partial U}{\partial F_i}\]

应变能对某个力求偏导，等于该力作用点沿该力方向的位移。

为什么好用：

不必解微分方程，只需积分再求导；
特别适合复杂结构（曲杆、刚架、桁架）与组合变形；
想求某处没有外力作用的位移？加一个虚设力 \(F_0\)，求导后令 \(F_0=0\) 即可。这个技巧极其实用。

单位载荷法（莫尔积分）是它的等价形式，土木与机械都常用：

\[\delta = \int \frac{M \bar{M}}{EI}dx\]

\(\bar{M}\) 是单位虚载荷引起的弯矩。

三、超静定问题的统一解法

第一页说超静定需要"平衡 + 变形协调 + 本构"三条腿。能量法把它变成一个标准流程：

力法（柔度法）：

解除多余约束，代之以未知力 \(X_1, X_2, \dots\)，得到静定基本系统；
写出总应变能 \(U(X_i)\)；
用变形协调条件——被解除处的位移为零：

\[\frac{\partial U}{\partial X_i} = 0\]

解出 \(X_i\)。

这是一个漂亮的转化：把"几何协调"变成了"能量对多余未知力取极值"。

最小势能原理是它的一般化：对满足边界条件的可容许位移，稳定平衡态使总势能在该状态附近取局部最小值；一般平衡首先给出驻值，非稳定平衡不满足这个“最小”表述。

\[\Pi = U - W_{\text{外力功}} \to \min\]

这就是有限元法的出发点（第二十一页）：把位移场用有限个节点参数近似，让总势能对这些参数取极小，就得到刚度方程 \(\mathbf{K}\mathbf{u}=\mathbf{F}\)。从这个意义上说，有限元法是能量法的数值化。

四、冲击与动载

在无耗散、线弹性、单自由度等效且接触过程明确的理想化下，能量法处理冲击特别方便：假设冲击过程中动能全部转为应变能。

对重物 \(W\) 从高度 \(h\) 落下冲击弹性体：

\[K_d = 1 + \sqrt{1 + \frac{2h}{\delta_{st}}}\]

\(K_d\) 是动荷系数，\(\delta_{st}\) 是静态挠度。动应力 \(= K_d \times\) 静应力。

两个极其重要的推论：

① 突加载荷（\(h=0\)）时 \(K_d = 2\)——突然放上的载荷产生的应力，是缓慢放置的两倍。这个"2 倍"必须记住。

② 结构越"软"（\(\delta_{st}\) 越大），动荷系数越小。

在上述落体—线性弹簧模型里，增加柔度会增大可储存的变形能、降低理想峰值力；这不是“越软越安全”的普遍结论，真实设计还要满足强度、行程、稳定、碰撞和疲劳约束。

工程上到处是这个原理的应用：

汽车的溃缩吸能区（故意做软）；
弹簧垫圈、缓冲器、橡胶减震块；
登山绳必须有弹性（静力绳用于冲坠会伤人）；
在材料、截面、连接边界和允许行程相近的比较中，长螺栓往往比短螺栓更能容纳冲击变形（同样变形量下可储存的应变能更大）；仍需校核预紧、屈服、疲劳和连接滑移。

"刚"未必安全——这是与"强度越高越好"同类的直觉误区。

五、弹性力学入门

材料力学靠假设（平截面假定、单向应力等）简化问题；弹性力学则从一般的三维方程出发：

三组方程（正是第一页"三条腿"的严格形式）：

\[\underbrace{\sigma_{ij,j} + f_i = 0}_{\text{平衡}}, \qquad \underbrace{\varepsilon_{ij}=\tfrac{1}{2}(u_{i,j}+u_{j,i})}_{\text{几何}}, \qquad \underbrace{\sigma_{ij}=\lambda\varepsilon_{kk}\delta_{ij}+2\mu\varepsilon_{ij}}_{\text{本构}}\]

15 个未知量（6 应力 + 6 应变 + 3 位移）、15 个方程，加上边界条件构成完整的边值问题。

可直接写出的解析解有限（例如厚壁筒、圣维南扭转、孔边应力集中、赫兹接触等经典边界），复杂几何因而常需要数值解——这正是有限元存在的理由。

两个必须知道的经典结果：

① 圣维南原理：在线弹性、小变形、边界条件兼容且远离载荷施加区和几何不连续处，应力的远场分布对载荷的局部施加方式不敏感，主要由合力和合力矩决定。

这个原理的实用价值极大：它使我们能用简化的载荷模型分析远场，但“约一个特征尺寸”只是与几何和载荷分布相关的数量级经验，不是普适常数；尖角、孔、接触和边界层附近仍必须保留局部模型。

② 赫兹接触：在经典、无摩擦、光滑弹性体和局部椭圆/圆形接触近似下，接触区应力分布由压力、曲率和弹性常数决定；某一常用切应力分量的最大值可在表面以下约 \(0.5a\)（\(a\) 为接触半宽）的量级出现。这个位置随摩擦、椭圆率、材料和接触边界改变，不能当成所有轴承的固定深度；它只提供滚动接触次表面疲劳的理想解释。

六、要点

应变能的统一形式：\((\text{内力})^2 L/(2\times\text{刚度})\)；
卡氏定理用求导代替解方程，虚设力技巧可求任意点位移；
能量法把超静定的几何协调转为能量极值；对稳定平衡，最小势能原理是有限元的出发点；
在线性弹簧的突加载荷模型中动荷系数 = 2；柔度增加有时能降低理想峰值，但必须同时检查强度、行程、稳定和接触；
圣维南原理在合适的远场条件下弱化局部载荷细节；经典赫兹接触可解释次表面剪应力，但不替代真实轴承接触分析。

线一结束：从力到应力、到变形、到失稳、到能量。下一线让结构动起来——动力学与振动，机械工程中最常见的现场问题来源。