ABAQUS车轨耦合动力学模型。 读研期间建立的高精度车轨耦合动力学模型，下图依次为整车模型、轮轨力、钢轨加速度时程曲线及钢轨位移时程曲线。
最近在整理读研期间做的车轨耦合动力学模型，顺手把几个关键实现点整理成文。这个模型通过ABAQUS搭建，包含了整车-轨道-地基的完整系统，主要用来研究轮轨接触振动和结构疲劳问题。（模型图、数据曲线图见文首）
整车模型采用多刚体动力学思路，车体、转向架、轮对之间用弹簧阻尼连接。这里有个小坑：直接导入三维模型会导致自由度爆炸。我的处理办法是先用HyperMesh做几何清理，保留关键连接点拓扑结构，再通过*CONNECTOR SECTION定义铰接关系。
ABAQUS车轨耦合动力学模型。 读研期间建立的高精度车轨耦合动力学模型，下图依次为整车模型、轮轨力、钢轨加速度时程曲线及钢轨位移时程曲线。
轨道子系统是整个模型的核心。钢轨用三维实体单元建模时发现计算量太大，后来改用欧拉梁单元+接触面组合方案。这里有个有意思的细节：用Python脚本自动生成扣件系统的非线性弹簧特性矩阵，直接嵌入到*SPRING卡片的INPUT参数里。
def create_fastener_spring(): stiffness_curve = [ (0.0, 0.0), (0.5, 1.2e6), (1.0, 2.5e6), (2.0, 3.8e6) ] with open('fastener.inp','w') as f: f.write('*SPRING,ELSET=Fastener_Spring\n') for strain, stress in stiffness_curve: f.write(f'{strain}, {stress}\n')轮轨接触模型用了Hertz理论+显式动力学求解。实测发现单纯用ABAQUS自带的接触算法会出现"数值颤振"，特别是在轮缘接触区域。后来在接触属性里加了个动态摩擦系数公式才算稳定：
*SURFACE INTERACTION, NAME=WHEEL_RAIL *FRICTION, DEPVAR=1 0.3, 0.0 ! 静摩擦系数 0.25, 0.5 ! 动摩擦系数随滑移速度变化求解器参数设置是另一个重头戏。动力分析步长需要兼顾计算效率和精度，通过试算发现当Δt=1e-5秒时，钢轨振动频率到2000Hz还能保持稳定。这里有个取巧的设定：把车轮作为解析刚体，用*RIGID BODY约束转动自由度，计算速度直接提升40%。
仿真结果里钢轨加速度曲线（图3）出现的高频毛刺，后来发现是扣件阻尼参数不够导致的。调整阻尼比从0.03增加到0.07后，2000Hz以上的振动分量明显衰减。位移时程曲线（图4）里的低频波动，倒是和现场实测的轨道不平顺数据吻合上了。
这个模型现在被师弟用来做轨道裂纹扩展预测，据说把损伤力学模块加进去之后，单工况计算时间从12小时飙到30小时。看来下次得试试子结构模态综合法了...
