通过双重子结构化对大规模结构进行非线性模型降阶

为何将大型数字结构缩小很重要

工程师常常模拟大型结构（如工厂、桥梁或飞机框架）在风、地震或机械作用下的振动与摆动。这些数字化试验可能包含数十万运动点，在强大的计算机上运行也需要数小时乃至数天。本文提出了一种方法，将如此庞大的模型缩减为仍能表现出原始行为的较小模型，即便结构存在棘手的强非线性节点和真实、复杂的阻尼形式。

把大结构拆成小块

出发点是观察到大多数结构由重复部件组成：相似的框架、楼层或面板。该方法不是一次性处理整个结构，而是将其划分为子结构。每个子结构单独分析，然后通过它们共享边界处的力将其重新连接。这一理念称为子结构化，长期以来已被用于较为简单的线性系统，其中响应与施加的载荷成正比。本工作的贡献在于提出了一种处理更现实行为的方式：某些节点或连接呈现非线性行为，且能量耗散（阻尼）不再遵循简化的教科书式模型。

用简单模式捕捉复杂运动

为了在不丢失重要物理特性的情况下缩减每个子结构的规模，作者采用了称为非线性固有模态（nonlinear normal modes）的概念。从本质上讲，模态是结构偏好振动的特征方式。对于线性系统，这些模态是直观且规律的模式。当运动幅度变大或节点表现为以三次响应而非简单线性响应的刚性弹簧时，这些模式会弯曲变形。文章遵循一套数学方案，将每个非线性模态表示为在全体可能运动空间中的光滑曲面。子结构中每一点的运动用界面处少量关键位移和速度的多项式表示。这样就能把庞大的变量集转成极为紧凑的描述，同时保留节点的非线性特征。

保持静态平衡与真实阻尼

该方法将每个子结构的响应分解为动态部分（非线性模态所在）和静态部分（处理由界面力引起的缓慢变形）。对于静态部分，方法借鉴了称为双重 Craig–Bampton 方法的既有框架。在该方法中，子结构之间的兼容性通过界面力来强制，而不是直接将边界位移粘合在一起。这带来了更小的矩阵和更灵活的拼接方式。本工作的一项重要改进是将一般形式的阻尼直接保留在方程中，而不是假定阻尼仅与质量或刚度成比例。因此，降阶模型能够真实地模拟配备附加阻尼器或以非均匀方式耗散能量的结构。

在数字化工业建筑上的验证

为证明该方法的实用性，作者将其应用于一个钢结构工业建筑的详细模型。建筑框架包含以扭转弹簧建模的节点，其抗力随旋转的立方增长，属于强非线性形式。建筑在水平方向受到接近其某一固有振动频率的正弦激励。首先，用标准时间步进算法求解完整有限元模型，耗时数百秒并占用数百兆字节内存。随后将建筑划分为重复的框架子结构和剩余部分。对于框架，仅保留四个非线性模态，聚焦于最关键节点的水平运动和扭转。求解该降阶系统得到的位移历时几乎与完整模型重合，同时计算时间约减少三分之二，内存占用大幅下降。

为何较少模态仍能给出可信答案

研究还探讨了精度如何随非线性模态的数量和选择而变化。当仅使用一个模态时，预测运动的误差较大。加入第二个直接涉及具有立方行为节点的模态后，误差显著下降，突显了纳入非线性最强处自由度的重要性。使用三个和四个模态时，误差继续下降到工程设计中被认为非常小的水平，同时模型保持紧凑。第二组仿真在外部加入阻尼器，产生高度非比例的阻尼模式。即便在这一更苛刻的情形下，降阶模型仍能紧跟全模型解，并继续在时间和内存上带来可观节省。

这对未来数字化结构意味着什么

通俗地说，本文展示了如何将笨重的数字化建筑变为一个敏捷的替代体，其在受振时的反应几乎与原结构相同，即便节点行为复杂、非线性且能量损失不规则。通过结合子结构化、非线性振动模态和考虑阻尼的表述，该方法为非常大规模结构的快速而可靠模拟打开了大门。这能帮助工程师运行更多情景分析、优化设计并探索新材料和装置，而不受过高计算成本的制约。

引用: Flores, P.A. Nonlinear model reduction for large-scale structures via dual substructuring. Sci Rep 16, 9286 (2026). https://doi.org/10.1038/s41598-026-38015-7

关键词: 结构动力学, 模型降阶, 非线性振动, 有限元分析, 子结构化