《European Journal of Mechanics - A/Solids》：Cascade of restabilizations along linear fundamental paths of axially soft multi-degree-of-freedom statical systems
编辑推荐：
分析了轴向极软的受压梁状结构从非平凡基本路径出发的分岔。研究人员证明，对于一大类受Euler Elastica梁置于Winkler地基上启发的系统，在某些情况下这些系统表现出多重再稳定（multiple restabilizations），即当载荷增加时它们反
分析了轴向极软的受压梁状结构从非平凡基本路径出发的分岔。研究人员证明，对于一大类受Euler Elastica梁置于Winkler地基上启发的系统，在某些情况下这些系统表现出多重再稳定（multiple restabilizations），即当载荷增加时它们反复失去并重新获得稳定性。当浏览一个以轴向刚度参数参数化的系统族时，再稳定的次数会变化，从零（不存在临界载荷，即系统对任意载荷值均稳定），到单次再稳定（存在唯一的不稳定区间），到两次再稳定（即存在两个不稳定区间），……直至N?1次再稳定，其中N?1为梁状系统的横向自由度数目。这些结果可推广到连续梁，此时N→∞。以下考虑成立。1. 不可伸长梁状结构的临界载荷取决于长度。当梁可伸长时，随着载荷增加，长度随载荷减小，因此存在不止一次发生分岔的情形。一个简单的图，如图1所示，提供了有用的工具来预先确定是否以及发生多少次再稳定。2. 地基的作用对于多重再稳定的发生至关重要。实际上，这与同一载荷下不同波长的多个临界模态的共存有关。置于弹性地基上的简支可伸长梁清楚地说明了同样的支撑驱动的再稳定级联。3. 再稳定的机制可追溯到切向刚度特征值在载荷参数临界值处的合并（coalescence）。它们组织了沿非平凡基本路径的稳定性丧失与恢复的序列。因此，该力学现象仍然依赖于Feodosiev所发现的机制，但由于存在多个自由度，它重复多次，从而产生再稳定的级联。4. 此处提供的简单解释可能消除与多重再稳定相关的神秘性，这一现象自80年代以来一直吸引学者，迄今为止通常被视为反直觉的。另一方面，此处仅检测到一个充分而非必要条件，即不排除存在其他类别结构的其他再稳定机制，能够重新打开这一谜团。本研究应通过未来工作加以实施，例如沿以下方向开展。1. 此处假定基本非平凡路径为线性，以获得其简单表达式。然而，由于再稳定需要非常大的缩短，假定轴向应变的非线性本构律将更为现实。考虑这一点的最简单模型将包括一个三次项，这将使基本路径变为非线性。除了基本路径不再可能获得解析表达式之外，通过仍使用图1中的图解构造可以猜测，硬化型定律将通过使再稳定在低缩短水平下可行和/或跳过其中一些再稳定，来促进再稳定的触发。更复杂的本构律，例如表达非局部行为的本构律，如Challamel等（2025）、Luongo等（2025）和Elishakoff等（2026）所使用的，也可以探索以检测它们对再稳定的影响。2. 本文研究的简支情况有助于以简单方式说明该现象，因为它允许获得闭式解。然而，它从另一个角度受到限制，因为已发现可伸长梁的屈曲模态与不可伸长梁的屈曲模态相同（尽管临界载荷不同）。这被认为是正弦模态的一个特性，其波长随载荷分段恒定。怀疑不同的边界条件，在梁不可伸长时允许临界模态随载荷连续变化，会在梁可伸长时改变其波长。此外，对于特定边界条件，可能会发生有限次数的再稳定。为回答最后这些问题，需要进一步的工作。
研究背景：弹性结构屈曲中预临界非线性效应近年来受到关注，尤其是Feodosiev模型（轴向可变形倒立摆受压）的研究表明，当摆的轴向刚度逐渐减小时，两个临界载荷相互接近、合并然后消失。该现象与预临界应变密切相关，预临界应变使刚度依赖于载荷乘子μ，从而导致非线性特征值问题，这在经典处理中通常被忽略。已有研究表明，在变长度一维结构中存在屈曲载荷消失和再稳定的现象，且最多观察到两次再稳定。研究人员提出疑问：是否存在表现出三次或更多次再稳定的结构，即再稳定级联。为此，研究人员开展了一项关于轴向极软多自由度静力系统线性基本路径上再稳定级联的研究。
研究人员考虑了一类m+n自由度系统，其总势能由弹性势能和外力功组成，其中轴向位移为大变形、弯曲位移为小变形，灵感来自Euler Elastica可伸长梁。研究人员证明了对于一大类受Euler Elastica梁置于Winkler地基上启发的系统，在某些情况下表现出多重再稳定，即当载荷增加时反复失去并重新获得稳定性。再稳定次数从零到N?1次不等，其中N?1为梁状系统的横向自由度数目，这些结果可推广到连续梁（N→∞）。该研究的重要意义在于揭示了再稳定级联的机制，消除了多重再稳定的神秘性，并指出地基（支撑）作用是多重再稳定发生的关键因素。
关键技术方法：研究人员采用总势能驻值原理建立平衡方程，通过引入增量变量对平衡方程线性化，得到非线性特征值问题A(μ)?=0。以支撑可伸长N摆为原型模型，通过有限差分方法将其与连续简支可伸长梁在Winkler地基上的屈曲问题进行类比。引入无量纲刚度参数κ（轴向-弯曲刚度比）和χ（支撑-弯曲刚度比），将特征方程分解为N?1个关于μ的二次多项式，每个对应一个屈曲模态。通过数值求解特征方程实根确定临界载荷，绘制(κ,μ)平面上的稳定性图。样本队列来源为离散链式模型和连续梁模型，无实验队列。
研究结果：
1. 引言与Feodosiev模型解释：通过Feodosiev模型说明了单次再稳定的机制——当轴向刚度较大时存在两个临界载荷，基本路径在Pc1处失稳、在Pc2处恢复稳定；随着轴向刚度减小，两临界载荷合并；更小时系统对任意载荷均稳定。将该机制推广到具有多条屈曲分支的结构，若直线与特征曲线多次相交则发生多重再稳定。
2. 一类轴向柔顺弹性系统：建立了m+n自由度系统的总势能表达式，推导了平衡方程和线性化分岔方程，指出分岔条件为非线性特征值问题A(μ)?=0，其特征值数目可能小于、等于或大于m+n，且特征值可随参数变化而合并。
3. 支撑可伸长N摆：以N根刚性伸缩杆组成的简支链为原型，推导了总势能、平衡方程、基本路径和屈曲特征值问题，得到四阶差分方程及其与连续简支梁在Winkler地基上屈曲问题的对应关系。通过引入矩阵D和D2，将屈曲算子表示为D的二次多项式算子。
4. 稳定性图与再稳定级联：通过对三重摆（N=3）和四重摆（N=4）的数值分析，发现无支撑链（χ=0）仅发生单次再稳定（Feodosiev型）；足够强的横向支撑（χ=200或2500）可使各模态不稳定区间分离，沿特定κ值的基本路径出现S,U,S,U,...交替模式，三重摆发生两次再稳定，四重摆发生三次再稳定。模态按波数从高到低依次激活，再稳定级联随N增大而加深。
5. 结论与未来工作：研究结论为：不可伸长梁状结构的临界载荷取决于长度，可伸长时长度随载荷减小，故存在多次分岔；地基作用是多重再稳定发生的关键，与同一载荷下不同波长多个临界模态的共存有关；再稳定机制追溯为切向刚度特征值在临界载荷参数处的合并；简单的解释可能消除多重再稳定的神秘性，但仅检测到充分非必要条件。未来工作方向包括：采用非线性本构律（含三次项）使基本路径非线性化；研究不同边界条件对再稳定次数的影响；探索其他再稳定机制。
