新论文：有限元携手物理引擎，倒塌模拟新思路 / 四六文摘

论文链接：

https://www.mdpi.com/2076-3417/10/12/4408

图1 某真实复杂图书馆倒塌全过程模拟

推送浓缩版

迈阿密事故分析后，很多读者对我们可视化倒塌模拟的技术比较感兴趣，希望我们介绍一下这个方法，所以我们整理了本文。

合理、真实的建筑倒塌废墟虚拟场景构建对地震应急和救援训练至关重要，但是现有研究中针对结构倒塌后的废墟分布研究较少。

现有的有限元与物理引擎方法存在以下两个方面的问题：（1）有限元方法模拟结构强非线性、超大变形倒塌过程计算成本高，且难以考虑块体碰撞。（2）基于物理引擎的倒塌模拟方法虽然具有结构倒塌的大变形模拟优势，但是在结构小变形阶段计算不准确。而结构的倒塌过程往往对于倒塌起始时刻的初始状态十分敏感，因而仅采用物理引擎方法难以保证结构倒塌模式的真实性。

因此，我们提出了联合有限单元法与物理引擎的建筑倒塌及废墟场景模拟框架——在小变形阶段使用有限元程序模拟，在大变形阶段使用基于物理引擎的方法进行模拟。得到的主要结论如下：

（1）所提出的联合模拟框架有效结合了有限元方法与物理引擎的优点。

（2）联合模拟方法考虑了结构构件的碰撞，并可以模拟结构坍塌后废墟的分布。

（3）利用振动台实验验证了所提出的联合模拟方法比单独使用物理引擎的方法更准确。

（4）利用真实复杂图书馆建筑的案例分析验证了联合模拟方法的适用性。

研究背景

地震中由建筑物倒塌产生的废墟使大量人员遭到埋压，严重威胁人员生命安全。相关统计表明，埋压是地震伤亡的主要因素之一。合理、可靠的建筑倒塌废墟虚拟场景构建，可以增强受训者的临场感觉，提升地震应急人员的快速反应和应急处置能力。

图2 地震后建筑倒塌废墟

左图图源：http://tse4.mm.bing.net/th/id/OIP.yUIqpebwWM2eW0sztQQyJQAAAA?pid=Api&rs=1

右图图源：

http://tse4.mm.bing.net/th/id/OIP.s1MWTej7TRRC2gkQzlUHKQHaFC?pid=Api&rs=1

现有基于数值模拟的倒塌后废墟分布研究中，常用的方法包括：离散元（DEM）、应用元（AEM）和物理引擎（Physics Engine）等（注意：有限元（FEM）因为结构强非线性、超大变形倒塌过程计算成本高，且难以考虑块体碰撞，而较少用于倒塌后废墟分布研究）。

图3 有限元大变形发生单元畸变

其中，物理引擎是一种用于计算复杂物理效应的计算机图形学技术，可模拟真实世界刚体的碰撞、并且部分物理引擎对电脑计算能力要求相对较低。近年来，在结构倒塌过程及废墟分布模拟中逐渐得到应用。

那么，物理引擎是否可以横扫倒塌及废墟模拟了呢？

答案当然是否定的，相关研究及实验表明，基于物理引擎建筑倒塌模拟方法虽然具有结构倒塌的大变形模拟优势，但是在结构小变形阶段计算不准确。而结构的倒塌过程往往对于倒塌起始时刻的初始状态十分敏感，因而仅采用物理引擎方法难以保证结构倒塌模式的真实性。

图 4 物理引擎模拟难以保证模拟的真实性

那么，如何同时保证结构小变形准确性与大变形倒塌模拟的可行性呢？

图5 非正经配图

考虑到有限元方法在建筑小变形阶段的非线性响应模拟精度的优势，我们提出了联合有限单元法与物理引擎的建筑倒塌及废墟场景模拟框架——在小变形阶段使用有限元程序模拟，在大变形阶段使用基于物理引擎的方法进行模拟。

有限元与物理引擎联合模拟框架

图6 联合模拟框架

图中倾斜加粗的字体表示所用到的程序，倾斜字体表示所用到的该程序的功能。

该框架包含四个模块：

模块一：有限元模型建立与转换

使用有限元软件建立结构的有限元模型，将有限元几何模型转化为Blender中的Bullet物理引擎几何模型，并对几何模型进行刚体化。

模块二：小变形模拟与倒塌起始时刻确定

开展有限元非线性时程分析，获得结构在倒塌小变形阶段的结构响应。确定结构的倒塌起始时刻，提取结构有限元模型在倒塌起始时刻及之前各个时刻各点的位移和速度。

模块三：数据迁移与大变形倒塌模拟

根据有限元计算数据映射方法，将倒塌起始时刻结构的位移数据和速度数据映射到Bullet物理引擎模型中。然后采用BCB（Bullet Constraints Builder）预处理工具对几何模型进行预处理（具体方法参考论文），从而得到可用于物理引擎计算的模型。最后，使用Bullet物理引擎进行计算，获得结构的大变形的计算数据。

模块四：数据整合渲染与可视化

大变形倒塌分析完成后，使用Blender的渲染功能分别对结构小变形和大变形结果进行渲染，最终将渲染好的视频拼接得到结构倒塌的全过程。

关键技术

注意：此部分过于细节，长话短说，这里仅简要介绍，具体实现详见论文。亦可跳过此节，从案例分析看起。

为了实现上述联合模拟框架，需要解决以下三个关键问题：

（1）物理引擎与有限元的几何模型及材料映射方法

物理引擎的倒塌模拟方法为刚体弹簧元，其单元与目前主流的有限元软件的单元有较大的差别，因此需要提出一种针对物理引擎与有限元软件的模型转换方法。

（2）倒塌起始时刻确定方法

要实现有限元与物理引擎的联合计算，需要根据有限元模拟结果确定有限元与物理引擎联合计算的切换时刻 (倒塌起始时刻)。

（3）有限元计算数据映射方法

有限元计算完成后，需要将有限元小变形计算数据迁移到物理引擎中。

3.1

几何模型及材料映射方法

（1）几何模型映射

采用Blender中提供的基本网格——立方体(cube)作为基础网格进行物理引擎模型构建，几何模型映射分两步：1）形状变换与位置变换：利用Blender提供的基本网格进行形状变换与位置变换，得到与有限元模型单元相同形状和位置的几何单元。2）单元刚体化：将几何模型转换为受物理引擎控制的刚体模型。具体实现详见论文。

图 7 Blender中实体结构建立方法

（2）BCB中材料属性映射方法

BCB会根据Blender中单元名称对单元进行分组。因此，在Blender中建立几何模型时，需要对单元的名称进行特别命名，以确保有限元中材料相同的单元被BCB划分到同一材料组中。具体命名方法详见论文。

3.2

倒塌起始时刻确定方法

本文的联合模拟方法为充分发挥有限元与物理引擎的优势，在倒塌起始时刻Tc前，采用有限元模型计算；在倒塌起始时间点Tc之后，采用物理引擎模型计算。因此，倒塌起始时刻的确定十分重要。部分研究表明可将竖向位移为1m作为倒塌判据。由于这一倒塌判据物理意义明确，因此本文将其选为倒塌起始时刻Tc的确定依据。

3.3

有限元计算数据映射方法

有限元数据映射包括位移场映射与速度场映射。为了实现从有限元到物理引擎BCB的位移场与速度场映射，需要使用Blender的关键帧插入技术。关键帧可定义物体某种状态的起止点，比如：定义位移起止点的位移关键帧与定义运动模式起止点的运动模式关键帧。Blender中有两种运动模式：按照动画系统运动，即物体的运动受位移关键帧控制；按照物理引擎模拟运动，即物体运动受物理引擎计算控制。

（1）位移场映射方法

位移场映射仅需使用位移关键帧，结构到达倒塌起始时刻Tc之前，其构件运动由位移关键帧控制，按照有限元计算结果运动。Blender仅仅完成有限元计算结果的显示，Bullet物理引擎未参与计算。

具体的位移场映射方法如下图所示。

图 8 位移场映射方法

（2）速度场映射方法

速度场映射的目的保证在倒塌起始时刻Tc，有限元模拟切换到BCB模拟时，不仅保证位移一致（采用上一节位移场映射方法），而且速度也保持一致。有限元计算的速度场数据的映射需要使用位移关键帧与运动模式关键帧。

Blender中的速度映射只能通过相邻位移关键帧进行位移差分得到。如果直接使用有限元倒塌起始时刻的位移场xj及倒塌起始时刻前一时刻的位移场xi进行速度映射，那么得到的是近似的割线速度，而非倒塌起始时刻有限元模拟的真实速度。而物理引擎中下一时刻物体的位置，由上一时刻物体的位置和速度确定，因此采用割线速度会导致误差。

因此，本文提出了一种虚拟位移向量的计算方法，将虚拟位移向量xvirtual替换xi，即可实现倒塌关键帧的真实速度场与位移场的映射。其大致思路如图9所示，具体实现详见论文。