回弹分析
回弹的理论基础
回弹是什么
两者不同。都是非线性座屈的一种，但荷载-变位曲线的形状不同。
穿跃（snap-through） — 荷载达到极大值（极限点）后，荷载下降而变位继续增加。荷载-变位曲线呈"山"形。Riks法可追跡。
回弹（snap-back） — 荷载下降，变位也回复。荷载-变位曲线呈"S形"或"环形"。通常的Riks法难以追跡。
变位会回复？这是说结构像弹簧一样反弹？
在准静态意义上是的。例如对浅拱进行加载，在某一点处会急剧反转并"啪"地飞向另一侧。此过程的荷载-变位曲线呈S形。动态上是瞬间的，但准静态上荷载和变位都会经过一个往返路径。
回弹的物理
哪些结构会产生回弹？
典型例子：
| 结构 | 现象 | 特点 |
|---|---|---|
| 浅拱 | 加载下反转 | 最经典的回弹 |
| 浅圆顶 | 外压凹陷 | 贝壳结构的回弹 |
| 双稳定贝壳 | 两个稳定形状间的转换 | 故意设计的回弹（变形结构） |
| 混凝土破坏 | 裂纹扩展时的荷载-变位 | 软化区域的回弹 |
| 分离（分层） | 界面破坏扩展 | 能量释放导致的回弹 |
混凝土破坏也会产生回弹？
是的。在混凝土拉伸试验中，裂纹扩展时荷载和变位都会减少。这就是回弹。要正确评估破坏能，需要准确追跡此路径。
数学分类
从数学上如何区分穿跃和回弹？
可以通过荷载-变位曲线的雅可比矩阵（切线矩阵）的性质来分类：
极限点（穿跃） — 荷载制御的雅可比矩阵 $\partial F / \partial u = 0$。荷载达到极值。变位单调增加。
回弹点 — 变位制御的雅可比矩阵 $\partial u / \partial F = 0$ 加上 $\partial F / \partial u = 0$。荷载和变位都转折。
几何上，荷载-变位曲线的切线变为竖直（$du/d\lambda = 0$）就是回弹点。穿跃处切线只是水平（$d\lambda/du = 0$）。
所以变位制御也追不过回弹，就像荷载制御追不过极限点一样。
完全正确。要通过回弹，需要荷载和变位都不作为制御量，需要另外的制御参数。
回弹的能量论
能否从能量角度解释回弹？
回弹可理解为应变能的释放。结构储存了应变能，某个触发点（微小荷载增加）会导致储存的能量瞬间释放。
荷载-变位曲线"折返"部分的面积对应能量释放量。回弹越剧烈，释放能量越大，动态响应（振动、冲击）越明显。
脆性破坏也会有相同的现象吧。
完全是的。格里菲思理论的能量释放率正是与回弹概念同根的。裂纹扩展时，储存的弹性能瞬间释放导致裂纹不稳定传播。
总结
让我总结一下回弹的理论。
要点：
- 回弹是荷载和变位都回复的不稳定现象 — 与穿跃（仅荷载回复）不同
- 浅拱、圆顶、破坏问题中出现 — 涉及广泛的结构问题
- 荷载制御和变位制御都无法追跡 — 需要特殊数值手法
- 理解为应变能释放 — 蓄积的应变能急剧释放
- 伴随动态响应 — 即使准静态分析也会追跡到动力学转换过程
荷载制御和变位制御都追不了，真是很棘手的问题。
正是因为这样，回弹的数值追跡是非线性力学中最具挑战性的问题之一。
回弹和变位制御的反向运动
回弹是荷载-变位曲线的变位方向也会"回复"的急剧转换点，与穿跃不同。外径20mm钢球压入软钢板的试验中，荷载峰值后变位会反向转换。1973年Crisfield、Willam、Riks（独立地）证明了弧长法可以追跡此路径，但物理上理解为"急剧能量释放"。
回弹的数值计算手法
回弹的数值追跡手法
回弹用通常的Riks法无法追跡吧。该怎么办？
通常的Riks法（球形弧长法）用弧长作为约束条件，但在回弹点处弧长方向不确定，容易追跡失败。有几种替代手法。
Cylindrical Arc-Length法
圆柱弧长法仅在变位空间定义弧长（不含荷载参数）：
相当于从通常Riks法的弧长约束中除去 $\psi^2 \Delta\lambda^2$ 项。
为什么这样就能追跡回弹？
因为回弹点处变位"回复"，荷载-变位空间的弧长就变得不唯一。但只在变位空间定义弧长，就能自然地追跡变位折返的方向。不过当曲线弯曲极端时收敛性会下降。
变位制御的技巧
选择特定自由度的变位作为制御量。不用在回弹点折返的自由度，而要选择单调增加的自由度作为制御量。
但未必总能找到这样单调增加的自由度啊。
拱的回弹中，虽然荷载点的变位会折返，但拱的端部水平变位可能一直单调增加。以这个水平变位作为制御量，就能追跡荷载点变位折返的路径。
Abaqus中可以用子选项监测荷载制御外的自由度，利用该自由度的单调性来追跡。具体地，可以通过调整 CONTROLS 中的 FIELD 参数，或利用 STATIC, RIKS 的节点/自由度指定。
能量制御法
能量释放率作为制御量的方法。特别对破坏力学的回弹（裂纹扩展）有效：
制御裂纹面积的增量 $\Delta A$，调整荷载使每一增量的能量释放率 $G$ 等于临界值 $G_c$。
这是专门针对破坏问题的手法吧。
是的。用CZM（内聚力区模型）进行层间剥离分析时如果出现回弹特别有效。Abaqus中与 *STATIC, STABILIZE 或 VCCT 组合使用。
动态分析的替代
也能用动态方式解回弹吗？
可以。回弹本质上是动力学现象，用显式或隐式动态分析追跡实际转换过程是最自然的方法。
优点：
- 没有收敛问题（显式法时）
- 能得到回弹时的振动、冲击响应
- 不需要特殊的弧长制御
缺点：
- 无法直接获得准静态荷载-变位曲线
- 减衰设定会影响结果
- 计算成本大（特别是隐式动态分析）
想要准静态荷载-变位曲线用弧长法，想要动态响应用动态分析，这样分开对吧。
完全正确。实务中通常两个都要做。先用弧长法把平衡路径掌握清楚，再用动态分析评估回弹后的实际响应。
人工粘性稳定化
Abaqus的 *STATIC, STABILIZE 能通过回弹吗？
某种程度上可以，但要小心。稳定化法在回弹点加人工粘性力防止发散。不过粘性太大会扭曲回弹的路径。
用稳定化法时的检查：
- ALLSD（稳定化散逸能）/ ALLIE（内部能）< 5%
- 与无稳定化Riks法的结果对比验证
- 稳定化系数的参数研究
总结
回弹的数值手法，整理一下。
要点：
- 通常Riks法对回弹无力 — 球形弧长法的局限
- 圆柱弧长法 — 仅在变位空间约束弧长
- 单调变位制御 — 找到不产生回弹的自由度来制控
- 动态分析 — 最自然的方法。追跡实际转换过程
- 稳定化法 — 实用但必须用能量比验证
- 手法组合最有效 — 不要固执于单一手法
回弹追跡的Crisfield弧长法
用FEM追跡回弹需要用Crisfield弧长法来同时制控变位方向。通常的弧长法同时增加荷载和变位，但回弹时变位符号反转，所以需要加入符号检测逻辑的"Modified Riks"法。Abaqus的δu_arc设定在回弹点通过时会自动反转变位增量方向。
回弹的实务应用
回弹分析的实务
实务中进行回弹分析的步骤是什么？
回弹往往难以提前预测。经常是"分析后才发现出现了回弹"的情况。
实务流程
1. 线性座屈分析 — 先用通常的特征值座屈掌握概况
2. 非线性分析（Riks法）试行 — 不收敛时可能是回弹
3. 检查荷载-变位曲线 — 看有没有变位折返
4. 选择手法 — 从圆柱弧长法、变位制控、动态分析、稳定化法中选择
5. 验证结果 — 能量平衡、网格收敛性
Riks法不收敛是回弹的信号吧。
没错。Riks法"往复振荡"（在同一点重复迭代）时就是回弹的征兆。
Riks法（先试行）
```
*STEP, NLGEOM=YES, INC=5000
*STATIC, RIKS
0.005, 2.0, 1e-15, 0.01, ,
```
重点：最小增量极小，最大增量受限制。
稳定化法（Riks法失败时）
```
*STEP, NLGEOM=YES, INC=5000
*STATIC, STABILIZE, ALLSDTOL=0.05
0.005, 1.0, 1e-15, 0.01
```
ALLSDTOL=0.05 将稳定化能的许容比设为5%。
动态分析（需要回弹后响应时）
```
*STEP, NLGEOM=YES
*DYNAMIC, APPLICATION=QUASI-STATIC
0.01, 1.0, 1e-8, 0.1
```
QUASI-STATIC选项自动调整数值减衰。抑制惯性效应同时通过回弹。
破坏问题的回弹
混凝土拉伸试验中的回弹怎么处理？
混凝土引张软化区，如果试验机刚度不足就会出现不稳定破坏（回弹）。FEM中也会出现同样现象。
对策：
- 裂纹张开位移（CMOD）制控 — 以裂纹张开量作为制控量。这是选择单调变位的一个例子
- 内聚力区模型 — 用CZM单元稳定地建模软化
- 粘性正则化 — 材料模型中加微小粘性进行数值稳定化。Abaqus的混凝土损伤塑性模型中内置
粘性正则化会"软化"回弹吗？
会的。粘性太大回弹就会消失，响应变得过于光滑。粘性参数的灵敏度分析是必须的，要确保结果不依赖粘性值。
双稳定结构的设计
双稳定结构（故意用回弹的结构）的设计呢？
双稳定结构中，两个稳定形状间的能量壁垒决定了回弹的"触发荷载"。设计要点：
- 两个稳定形状的能量差 — 一个明显低能，则容易回弹且难以返回
- 回弹的速度和冲击 — 动态分析评估回弹时的加速度、冲击力
- 重复回弹的疲劳 — 重复回弹会导致材料疲劳
实务检查清单
请给出回弹分析的检查清单。
用多种手法比较结果很重要。单一手法不够可靠。
回弹是数值上最麻烦的问题，所以用2种独立手法得到相同结果是最佳实践。
混凝土剪切破坏回弹
混凝土梁的剪切破坏中，荷载急剧下降后变位会反向转换的回弹挙动常见。1970年代Bazant发现了这一现象，与"尺寸效应（Size effect）"密切相关，梁越大回弹越剧烈。现在的混凝土结构物FEM破坏分析（裂纹有限元法）默认选用Crisfield法。
回弹的软件比较
回弹分析的工具
哪些求解器能追跡回弹？
回弹需要高端非线性分析功能。
| 功能 | Abaqus | Nastran | Ansys | LS-DYNA |
|---|---|---|---|---|
| Riks法 | *STATIC, RIKS | SOL 106 BUCKLE | Arc-Length | — |
| 圆柱弧长 | *CONTROLS中调整 | 有限支持 | ARCTRM | — |
| 稳定化法 | *STABILIZE | 有限支持 | STABILIZE | — |
| 动态分析（替代） | Explicit/Implicit | SOL 129/400 | Transient | 标准 |
| CZM（破坏） | *COHESIVE | — | CZM单元 | *COHESIVE |
Abaqus选项最多。
就回弹而言Abaqus优势明显。*STATIC, RIKS实现稳定，还有多个替代手段：*STATIC, STABILIZE 和 *DYNAMIC, QUASI-STATIC。遇到回弹问题应首选Abaqus。
LS-DYNA能自然地用动态分析追跡回弹。准静态荷载-变位曲线难得，但动态响应想要就最优。
选定指南
总结一下？
回弹是"非线性分析力"的考验。
完全同意。能否正确追跡回弹是测量非线性分析求解器实力的一个标尺。
OpenSees混凝土非线性分析
OpenSees是太平洋地震工程研究中心开发的开源结构分析框架，具有专为包含回弹在内的非线性分析设计的丰富材料模型。Concrete01、02、07材料模型能精确再现混凝土压缩、拉伸后座屈特性，被世界地震工程研究者广泛使用。
回弹的先端研究
回弹研究的最前线
回弹分析的先端研究是什么？
有3个活跃方向。
准脆性材料的回弹
混凝土等准脆性材料的回弹是？
混凝土拉伸软化、岩体破坏、陶瓷断裂...都会在荷载-变位曲线上出现回弹。
最大课题是网格依赖性。用FEM处理拉伸软化时，变形局部化于单个单元宽度，网格越细脆性响应越明显。
正则化手法：
- 裂纹带模型（crack band model） — Bazant(1983)提出。用网格大小补正软化参数
- 非局部模型 — 应力或应变空间平均化以防局部化
- 梯度损伤模型 — 包含应变梯度的高阶连续体模型
不正则化的话回弹的荷载-变位曲线就网格依赖。
是的。破坏能 $G_f$ 要保持一定来调整软化参数。位相场法（位相场模型）也是破坏回弹正则化的有力新手法，迅速发展中。
多稳定性（多稳定结构）
能有3个以上稳定形状的结构吗？
多稳定性可用形状记忆合金或预应力复合材料实现。各稳定形状间的转换就是回弹。
应用例：
- 可编程结构 — 外部刺激下任意切换形状
- 力学超材料 — 众多单元格回弹来制御宏观非线性响应
- 能量吸收 — 回弹滞后散逸能量
力学超材料的应用有意思。
加州理工和瑞士ETH的小组领导这个研究。用3D打印制作的多稳定超材料，每个单元格都会回弹，宏观上能实现可编程应力-应变曲线。
计算手法的进化
回弹的数值追跡手法在进化吗？
广义弧长法的改良在继续：
分岐追跡很难吧。
非常难。分岐点处多条平衡路径交差，选哪条路依赖微小扰动。数值上用特征值分析确定分岐方向，沿该方向加微小扰动后用Riks法追跡。
总结
回弹先端研究整理一下。
回弹是"结构不稳定性"的最复杂形态。正确追跡能力标志着非线性结构力学的成熟度。
回弹在超弹性材料的应用
具有回弹特性的超弹性材料（形状记忆合金等）能表现"负刚度效应"，展示极强振动吸收性。负刚度装置与正刚度并联结合时整体刚度≈0，外部振动无法传入结构。这被用于精密工厂、电子显微镜安装台的振动隔离，是回弹逆向应用的典范。
回弹的故障排除
回弹分析的故障
回弹分析遇到的故障和对策是什么？
回弹是非线性分析中最难收敛的问题。
Riks法振荡
Riks法在回弹点附近往复振荡无法收敛。
这是Riks法在回弹点弧长方向不确定的表现。
对策（按顺序试）：
1. 减小最大增量 — 0.01 → 0.005 → 0.001
2. 减小最小增量 — $10^{-12}$ → $10^{-15}$
3. 切换到圆柱弧长法 — 仅在变位空间约束
4. 切换到稳定化法 — *STATIC, STABILIZE
5. 切换到动态分析 — *DYNAMIC, QUASI-STATIC
从1开始逐一试。
通常的Riks法设定调整（1、2）就能解决。不行的话尝试手法切换（3～5）。不要固执于一种手法，灵活切换是回弹分析的诀窍。
稳定化法使回弹消失
用稳定化法后荷载-变位曲线变光滑，回弹消失了。
稳定化系数太大。人工粘性吸收了回弹的能量释放。
对策：
- 降低稳定化系数（FAQTOR：0.001 → 0.0001 → 0.00001分步）
- 检查ALLSD/ALLIE比（目标1%以下）
- 与Riks法结果对比验证
稳定化系数太小又会不收敛吧。
是的。稳定化系数需要找"甜蜜点"。太小不稳定，太大回弹被抑制。参数研究确定合适范围是必须的。
动态分析中回弹后振动过大
动态分析追跡回弹后出现大振动，后续分析不稳定。
回弹时能量释放转成运动能引起振动。对策：
- 增加减衰 — Rayleigh减衰（α、β）或Abaqus的*DYNAMIC, APPLICATION=MODERATE DISSIPATION增加数值减衰
- 分阶段卸载荷载 — 回弹点前先卸载再重载
- 质量缩放 — 增大质量降低固有频率，增大时间步长（注意惯性效应）
破坏分析中网格依赖性过大
混凝土破坏分析中改变网格后回弹形状完全变了。
对局部化缺乏正则化。
检查和对策：
- 用了裂纹带模型吗 — 软化曲线应随网格大小调整
- 特征长度正确吗 — Abaqus自动计算，但特殊单元形状时可能不准
- 破坏能 $G_f$ 守恒吗 — 换网格后 $G_f$ 应保持不变确认
破坏能守恒就能抑制网格依赖。
理想上是。实际上完全独立于网格很难，2～3个网格水平比较结果的差在10～20%以内就认可的实务做法。
总结
回弹分析故障排除整理。
回弹是"最难的非线性问题"，要谨慎对待。
完全同意。能正确追跡回弹就是非线性结构力学高手。慢慢来，一个手法接一个手法地尝试积累经验。
弧长法在回弹点发散
弧长法在回弹点发散往往是增量Δl过大"跨越"的结果。将Δl减为1/10，在回弹点附近（荷载差分符号改变处）增加增量数。Abaqus加入STABILIZATION参数后人工减衰稳定回弹点通过，之后残余路径追跡成为可能。
价值
详细
错误