钢轨裂纹扩展仿真案例：FRANC3D联合（DLOAD+UTRACLOAD）子程序操作汇总 原创

案例介绍：在ABAQUS中利用DLOAD+UTRACLOAD子程序来模拟铁路钢轨承受列车车轮的滚动接触载荷，不需要建立车轮，这两个子程序可以模拟车轮经过钢轨时的垂向载荷+切向载荷，对于只研究钢轨受力的同学，该方法节省不少工作量和计算时间。模拟结果如下图所示，对于子程序的代码，其他博主已经有偿详细介绍了有相关内容( https://www.jishulink.com/post/1975039 、 https://www.jishulink.com/video/c16437、https://www.jishulink.com/video/c245885)，为了不扰乱市场，在此本人不做介绍。当然本案例的代码是本人所写，有需要的话可以在我主页找我联系，同样有偿分享。为了了解大家对此案例的需求量，在此设置本站最低价格作为统计，大家按需下单即可。

本人在此分享一个针对DLOAD+UTRACLOAD对钢轨加载，并利用FRANCD进行钢轨裂纹扩展仿真操作流程的案例，详细步骤如下：

1.首先需要在ABAQUS中设置好钢轨模型，建立好插入裂纹的local模型单元集（set），如下图所示。

2.给钢轨添加底部约束和两端对称约束，轨顶上添加压强+表面载荷，方式都是用户自定义大小都设置成1即可，提交作业时需要用到DLOAD+UTRACLOAD子程序代码。并确保代码可以跑出下图的效果，即移动载荷实现滚入到滚出local模型。

3.确保ABAQUS的模型没问题后，需要创建一个inp文件给FRANCD，可以新建一个Crack1作业，可以不选子程序，直接点“写入输出文件”即可，在ABAQUS的工作路径中就会生成Crack1.inp文件。

4.通过Element Group的方式切割子模型，并保留连接面的网格。

5.在local模型的轨顶正中间，插入一个直径为2mm的圆形表面裂纹。

6.提交静态分析，先计算初始裂纹的应力强度因子，取名为C0，到最后一步时，选择Merge nodes，勾选Write files but NOT run analysis。会在FRANC3D的工作路径生成C0_full.inp，以及WriteDptfile.py.