铁木辛柯梁理论
理论与物理
欧拉-伯努利梁与铁木辛柯梁的区别
老师，铁木辛柯梁和欧拉-伯努利梁有什么区别？
最大的区别是考虑了剪切变形。欧拉-伯努利梁假设“截面始终垂直于中性轴”，而铁木辛柯梁去掉了这个假设。允许截面倾斜于中性轴。
数学上：
- 欧拉-伯努利梁: 转角 = 挠度的微分 → $\theta = dw/dx$
- 铁木辛柯梁: 转角 ≠ 挠度的微分 → $\theta \neq dw/dx$，其差值即为剪切变形
剪切应变 $\gamma$：
$dw/dx$ 是梁轴的斜率，$\theta$ 是截面的转角。这个差值就是剪切变形吧。
没错。欧拉-伯努利梁强制 $\gamma = 0$（无剪切变形）。铁木辛柯梁通过解除这个约束，可以描述更一般的梁行为。
控制方程
请告诉我铁木辛柯梁的微分方程。
是两个联立微分方程：
这里 $A_s = \kappa A$ 是有效剪切面积，$\kappa$ 是剪切修正系数。
剪切修正系数 $\kappa$ 是什么？
梁理论假设截面剪应力均匀分布，但实际的剪应力分布是抛物线状的。$\kappa$ 就是用来修正这个差异的系数。
| 截面形状 | $\kappa$ |
|---|---|
| 矩形截面 | 5/6 ≈ 0.833 |
| 圆形截面 | 6/7 ≈ 0.857 |
| 薄壁圆管 | 1/2 = 0.5 |
| 工字形截面（腹板剪切） | $A_w / A$（腹板面积/总面积） |
$\kappa$ 小于1…也就是说有效剪切面积小于全截面面积呢。
是的。工字钢的剪切主要由腹板承担，所以有效剪切面积接近腹板面积。翼缘主要贡献弯曲刚度，对剪切刚度贡献很小。
挠度的分解
铁木辛柯梁的挠度可以分解为弯曲和剪切部分吗？
可以。总挠度：
简支梁中央受集中载荷的情况：
第一项是欧拉-伯努利梁的挠度，第二项是剪切附加项呢。
剪切挠度比率：
这里 $r = \sqrt{I/A}$ 是截面的回转半径。$r/L$ 越大（梁越粗短），剪切变形的贡献越大。
对于钢材（$E/G \approx 2.6$）的矩形截面（$\kappa = 5/6$），当 $L/h = 10$ 时，剪切挠度占总挠度的百分之几？
计算后约为3%。$L/h = 5$ 时约为12%。$L/h = 3$ 时约为32%。$L/h < 5$ 时剪切变形就不可忽略这个经验法则就是从这个计算来的。
有限元法中的铁木辛柯梁单元
有限元法中的铁木辛柯梁单元有什么特别之处？
铁木辛柯梁单元将 $w$ 和 $\theta$ 作为独立变量处理。欧拉-伯努利梁单元有 $\theta = dw/dx$ 的约束，而铁木辛柯梁单元没有这个约束。
但需要注意。常规的2节点单元（线性插值）会发生剪切自锁。
剪切自锁？
这是与欧拉-伯努利梁的剪切自锁相反的现象。在铁木辛柯梁单元中试图表现纯弯曲变形时，会产生寄生的剪切应变，导致单元自锁。结果是挠度被低估。
对策：
- 减缩积分（1点高斯积分）—— 避免剪切自锁
- 假定应变法 —— 剪切应变的独立近似
- 高阶单元（3节点以上）—— 自锁现象自然消除
欧拉-伯努利梁单元有剪切自锁，铁木辛柯梁单元有剪切自锁…会发生相反的问题呢。
是的。这是有限元单元设计的永恒主题。欧拉-伯努利梁忽略剪切来避免自锁，铁木辛柯梁包含剪切并与自锁作斗争。实用的解决方案是减缩积分。
总结
我来整理一下铁木辛柯梁理论。
要点：
- 考虑剪切变形 — $\gamma = dw/dx - \theta \neq 0$
- 剪切修正系数 $\kappa$ — 取决于截面形状。矩形为5/6
- $L/h < 5$ 时显著 — 粗短梁、夹层板、组合梁
- 注意剪切自锁 — 用减缩积分避免
- 欧拉-伯努利梁的上位替代 — 铁木辛柯梁在 $L/h \to \infty$ 时收敛于欧拉-伯努利梁
是不是说，拿不准的时候就用铁木辛柯梁？
没错。铁木辛柯梁包含了欧拉-伯努利梁。对于 $L/h$ 大的细长梁，会得到与欧拉-伯努利梁相同的结果，所以一直使用铁木辛柯梁也没问题。因此 Abaqus 和 Ansys 的默认梁单元就是铁木辛柯梁。
铁木辛柯梁理论的诞生
斯蒂芬·铁木辛柯在1921年的论文《关于棱柱杆横向振动微分方程的剪切修正因子》中，将剪切变形引入了梁理论。他在俄国革命后流亡美国，并以在斯坦福大学革新工程教育而闻名。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重，越难启动，一旦启动也越难停止。建筑物在地震中摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的性质就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓的紧固力…都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在三维空间坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——故意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物会在地震后一直摇晃下去。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为 mm 时，载荷、弹性模量也需统一为 MPa/N 系单位 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm 制下为 tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm 制下用 N，m 制下也用 N 统一 |
数值解法与实现
铁木辛柯梁单元的公式化
请告诉我铁木辛柯梁单元的刚度矩阵。
对于2节点的铁木辛柯梁单元，仅考虑弯曲自由度（$w_1, \theta_1, w_2, \theta_2$）的情况：
这里 $\Phi = 12EI/(GA_s L^2)$ 是剪切变形参数。
令 $\Phi = 0$（无剪切变形）就与欧拉-伯努利梁的刚度矩阵一致了呢！
确认得很完美。$\Phi$ 可以看作是相对于欧拉-伯努利梁的修正量。$\Phi$ 越大（梁越粗短），剪切变形的影响越大，与欧拉-伯努利梁的差异也越大。
剪切自锁的详细说明
请再详细说明一下剪切自锁的机理。
考虑用2节点的线性插值单元来表现纯弯曲（$M$ = 常数）的情况。如果弯矩是常数，那么 $\theta$ 也应该是常数（$d\theta/dx = M/(EI)$），但 $w$ 的线性插值意味着 $dw/dx$ 也是常数。于是 $\gamma = dw/dx - \theta = \text{常数} \neq 0$。
明明是纯弯曲却产生了剪切应变…这就是寄生剪切呢。
是的。这个剪切应变会储存多余的能量，导致变形被低估（响应过于刚硬）。
减缩积分（1点高斯积分）能解决的原因：$\gamma$
相关主题
なった
詳しく
報告