12月28日，埃隆·马斯克的“无聊”公司展示了其新概念的电动汽车隧道，让隧道这个词在全世界火了一把。不过就现场情况来看，“以240km的时速风驰电掣”还只是一张大饼。向地下进军，挖隧道的想法是对的，但是目前解决大城市病的出路依然集中在轨道交通身上。
截至2017年底，全国投入运营隧道达14547座，总长15326km。今天，小编为大家分享一种打通基于达索3DEXPERIENCE平台的隧道地层结构模型数值分析流程，实现从BIM模型到计算模型的无缝衔接。
对于隧道专业而言，应用达索3DEXPERIENCE平台进行地层结构模型数值分
析的最大障碍就是土体塑性参数及本构模型不能添加。所以本文重点就是（敲黑板）实现模型塑性参数及本构模型的添加和将计算过程流程化。
添加塑性参数及本构模型
隧道结构数值分析经常使用地层
结构模型，地层结构模型中的土体材料为弹塑性，需要添加弹塑性参数及相对应的本构模型，这是模型能否进行
弹塑性计算的关键。
达索3DEXPERIENCE平台SIMULIA模块现有功能仅能添加材料弹性参数，由于缺少塑性参数及本构模型，对于隧道地层结构模型的弹塑性分析无法实现。
传统ABAQUS中添加塑性材料参数及本构模型
达索3DEXPERIENCE平台SIMU LIA模块现有材料参数
达索3DEXPERIENCE平台中的SIMULIA模块基于ABAQUS开发，因此二者计算模型文件的数据结构相同，对于一些在3DEXPERIENCE平台中无法实现的功能，可充分参考传统ABAQUS计算模型的数据结构，通过对3DEXPERIENCE平台的计算模型数据结构进行更改来实现。
达索平台SIMULIA模块的计算模型文件为.inp文件，此文件基于Python语言，因此工作重点就是应用Python语句，将.inp文件中相关的塑性参数和本构模型进行更改或添加。
计算模型数据文件内容繁杂，人工编辑难度较大，在掌握文件数据结构的基础上，可通过开发程序实现文件更改。对于模型塑性参数及本构模型的更改与添加，可通过下图流程实现。具体步骤为：依托达索3DEXPERIENCE平台建模模块CATIA创建BIM模型；将BIM模型切换到计算模块SIMULIA，生成计算模型；将计算模型导出为．inp文件，与此同时在.txt文件中编写塑性材料参数及本构模型；编写.py程序文件。执行此文件将.inp文件中的材料参数根据实际情况替换成.txt文件中的塑性材料参数及本构模型；生成新的.inp文件。
计算模型文件材料参数更新流程
以上流程实现的重点是在充分了解.inp文件数据结构的基础上，通过Python编写执行程序.py文件。如图4、图5所示，．py文件执行完成后，新生成的.inp文件中弹性材料参数已经更改为实际塑性材料参数，并添加了M—C塑性本构模型，更改后的数据结构和传统ABAQUS中直接添加塑性参数及本构模型生成的数据结构一致。
虽然在达索3DEXPERIENCE平台SIMULIA模块中没有相关塑性材料参数及本构模型的添加界面，但通过程序更改计算模型文件数据结构可以达到和传统ABAQUS相同的效果，这就为隧道地层结构模型在达索软件中的应用扫清了障碍。
初始及更新后的.inp文件材料参数
定义计算流程
传统数值分析软件塑性材料参数及本构模型在操作界面中完成，过程简单明了，而达索平台因为自身的局限性，需要首先更改生成的计算模型文件，然后再将更改后的文件导回软件进行计算，计算后如结果不能满足要求一般还会涉及参数的调整，如此反复，操作复杂，因此有必要将此过程流程化，减少操作，提高计算效率。
达索3DEXPERIENCE平台中的ProcessComposer模块，可将分散的操作过程进行组合，将逐步实现的功能流程化，减少中间操作。利用这一功能，可将隧道地层结构模型数值分析中计算模型文件的导出、修改、导入等单独执行的步骤进行整合，使之流程化，当进行数值分析时，只需执行此流程即可，操作简单，计算高效。此外，在此过程中还可以将计算过程中需要调整的参数提取出来，便于修改。
以添加地层结构模型地质塑性参数及本构模型为例，通过Processcomposer模块将计算过程中各分步进行整合。具体实现步骤为：下载本地文件，包括初始.inp文件、包含有参数信息的.fxt文件以及执行材料更新过程的.py文件等，.py文件是主程序文件，.inp和.txt文件是程序执行过程中被调用的文件；设置关键参数，可将计算过程中的关键参数提取出来便于后期修改；执行.py文件，将初始.inp文件按.txt文件中的材料参及数本构模型进行更新；生成并上传更新后的.inp文件；提交.inp文件到计算器，执行计算；生成并上传计算结果文件(.odb文件)；打开生成的.odb文件，查看计算结果；当结果不满足要求时，更改关键参数重新执行计算。
计算过程中各分步整合后的流程
依托某铁路项目隧道工程BIM模型及地质BIM模型对以上方法进行验证。分析模型选用隧道地层结构模型，分析过程中的本构模型采用M-C准则，分析内容为隧道结构拱顶沉降，主要实现过程如下。
选取待分析段落BIM模型
在达索3DEXPERIENCE平台中已建有某铁路项目隧道工程BIM模型以及地质BIM模型，在建模模块中打开此模型并截取待分析段落的隧道主体结构及地质体。由于BIM模型体现的是三维真实场景，因此较传统模式下的计算分析其计算精度更高。
需要计算的隧遁BIM模型
转换BIM模型为计算模型
将截取后的隧道主体结构及地质体模型切换到计算模块SIMULIA，按照添加材料参数一划分单元网格一定义分析类型一定义接触面一添加荷载一定义边界条件的顺序生成计算模型，由于计算模块中不能添加地质体塑性参数，故此步仅需添加弹性参数，塑性参数及本构关系根据前文方法在后续步骤中添加。
基于隧道BIM模型生成的计算模型
计算过程中需要的基础文件包括3部分：（1）通过计算模型导出的.inp文件；（2）定义材料参数的.txt文件；（3）执行材料更新的.py文件。其中.py文件是通过Python编写的执行文件，通过调用计算模型文件．inp和材料参数文件.txt，生成更新材料参数的.inp文件。
计算流程
下载本地文件一设置关键参数一执行.py文件一上传更新后的.inp文件一执行计算一上传计算结果文件一查看计算结果。
在达索3DEXPERIENCE平台Processcomposer模块中执行前文中定义的计算流程，并生成计算结果，当结果不能满足要求时，更改关键材料参数重新执行计算。
本案例应用本文提出的塑性参数添加方法及将计算过程流程化的方法，顺利完成了隧道地层结构模型数值计算，计算结果与实际情况相符，验证了此方法的可行性。
与传统计算方法对比分析
基于BIM模型进行数值分析的传统方法如下：将BIM模型从BIM建模软件中以一定的格式导出；将导出后的BIM模型导入传统数值分析软件；在传统数值分析软件中添加弹塑性参数、划分网格以及进行其他边界条件的定义；执行计算。
此时，BIM建模软件为达索软件，数值分析软件采用ABAQUS，BIM模型由达索软件的导出格式为.stp。由于BIM建模软件和数值分析软件的接口不同，当BIM模型导入数值分析软件后出现了严重的数据丢失，后续计算无法执行。
