Abaqus铁路轨道建模crtsⅠ型轨道模型CRTSⅡ型轨道模型crtsⅢ型轨道模型轨道不平顺模拟轨道不平顺插件车轨耦合车轨地基耦合模型动力响应分析轨道弹簧批量施加。兄弟们搞轨道建模的应该都懂CRTS这玩意儿就像变形金刚不同型号差得不是一星半点。今天咱们就唠唠Abaqus里怎么折腾这堆轨道模型特别是那要命的弹簧批量操作保不准能帮各位省两瓶生发液。先说轨道不平顺这事儿手工建模能累出腱鞘炎。我去年搞了个Python插件核心就这行代码profile [A*math.sin(2*math.pi*x/L phase) for x in np.arange(0,length,step)]这可不是普通的正弦波波长L取20m对应轨检车典型病害幅值A别超过4mm否则计算直接崩。记得用random.uniform给phase加扰动不然所有波峰对齐能让你动力响应曲线变成心电图室同款。Abaqus铁路轨道建模crtsⅠ型轨道模型CRTSⅡ型轨道模型crtsⅢ型轨道模型轨道不平顺模拟轨道不平顺插件车轨耦合车轨地基耦合模型动力响应分析轨道弹簧批量施加。CRTS III型板式轨道建模时CA砂浆层千万别傻乎乎建实体。实测用connector单元模拟剪切连接刚度系数取120MN/m³这样既能保证计算速度又能吃得住35吨轴重。关键代码长这样for i in range(1,slab_num1): slab mdb.models[Model-1].parts[Slab-str(i)] slab.engineeringFeatures.ConnectorSection( nameShear-str(i), translationalTypeAXIAL)车轨耦合这块有个坑轮轨接触的摩擦系数千万别照搬教科书。实测潮湿钢轨μ0.3时计算结果最贴近轨检车数据但记得在轮缘接触区域局部调高到0.45不然曲线段轮轨横向力能差出30%最近发现的骚操作是弹簧批量施加大法。用Abaqus的findAt配合numpy.linspace30秒搞定2000个弹簧points [(x, y, z) for x in np.linspace(0, length, 2000)] springs part.Set(namesprings, facespart.faces.findAt(coordinatespoints))重点是要在findAt前关掉提示窗不然弹窗能卡死程序。地基弹簧刚度别均匀分布离桥台20m范围内刚度得按1.2倍系数递增这样基础沉降模拟才不会翻车。最后说个血泪教训动力分析后处理千万别直接用默认的时程输出。用这段代码提取轨头加速度能省80%存储空间odb.steps[Step-1].setValues(timeIncrement0.002) field odb.steps[Step-1].frames[1].fieldOutputs[U].getSubset( regionrail_head_set)记住采样间隔别小于0.002秒否则300km/h工况下会漏掉关键频段。搞完这些你的轨道模型差不多能在普通工作站上跑起来了——当然咖啡还是得备足毕竟再优化也得算个把小时不是
设置法向刚度120e6
Abaqus铁路轨道建模crtsⅠ型轨道模型CRTSⅡ型轨道模型crtsⅢ型轨道模型轨道不平顺模拟轨道不平顺插件车轨耦合车轨地基耦合模型动力响应分析轨道弹簧批量施加。兄弟们搞轨道建模的应该都懂CRTS这玩意儿就像变形金刚不同型号差得不是一星半点。今天咱们就唠唠Abaqus里怎么折腾这堆轨道模型特别是那要命的弹簧批量操作保不准能帮各位省两瓶生发液。先说轨道不平顺这事儿手工建模能累出腱鞘炎。我去年搞了个Python插件核心就这行代码profile [A*math.sin(2*math.pi*x/L phase) for x in np.arange(0,length,step)]这可不是普通的正弦波波长L取20m对应轨检车典型病害幅值A别超过4mm否则计算直接崩。记得用random.uniform给phase加扰动不然所有波峰对齐能让你动力响应曲线变成心电图室同款。Abaqus铁路轨道建模crtsⅠ型轨道模型CRTSⅡ型轨道模型crtsⅢ型轨道模型轨道不平顺模拟轨道不平顺插件车轨耦合车轨地基耦合模型动力响应分析轨道弹簧批量施加。CRTS III型板式轨道建模时CA砂浆层千万别傻乎乎建实体。实测用connector单元模拟剪切连接刚度系数取120MN/m³这样既能保证计算速度又能吃得住35吨轴重。关键代码长这样for i in range(1,slab_num1): slab mdb.models[Model-1].parts[Slab-str(i)] slab.engineeringFeatures.ConnectorSection( nameShear-str(i), translationalTypeAXIAL)车轨耦合这块有个坑轮轨接触的摩擦系数千万别照搬教科书。实测潮湿钢轨μ0.3时计算结果最贴近轨检车数据但记得在轮缘接触区域局部调高到0.45不然曲线段轮轨横向力能差出30%最近发现的骚操作是弹簧批量施加大法。用Abaqus的findAt配合numpy.linspace30秒搞定2000个弹簧points [(x, y, z) for x in np.linspace(0, length, 2000)] springs part.Set(namesprings, facespart.faces.findAt(coordinatespoints))重点是要在findAt前关掉提示窗不然弹窗能卡死程序。地基弹簧刚度别均匀分布离桥台20m范围内刚度得按1.2倍系数递增这样基础沉降模拟才不会翻车。最后说个血泪教训动力分析后处理千万别直接用默认的时程输出。用这段代码提取轨头加速度能省80%存储空间odb.steps[Step-1].setValues(timeIncrement0.002) field odb.steps[Step-1].frames[1].fieldOutputs[U].getSubset( regionrail_head_set)记住采样间隔别小于0.002秒否则300km/h工况下会漏掉关键频段。搞完这些你的轨道模型差不多能在普通工作站上跑起来了——当然咖啡还是得备足毕竟再优化也得算个把小时不是
