《铁路轨道有限元分析.doc》由会员分享，可在线阅读，更多相关《铁路轨道有限元分析.doc（17页珍藏版）》请在金锄头文库上搜索。
1、首先，是用AutoCAD或ANSYS将轨道的几何模型建好，转成IGES格式，保存成文件，这一步工作比较繁琐但没有技术麻烦，为了方便大家，现将其上传。见附件。 217052-RAIL60.IGES (0.0k) 接着建模型导入ANSYS： FileImportIGES 定义单元1：Mesh200-用于对平面的单元划分，注意要在Options中设置单元为4节点平面单元（QUAD 4-NODE） 定义单元2：Solid45-用于实体的单元类型。 单元的形成是二维拖拉成三维！ 定义轨道材料常数：弹性模量210GPa，泊松比0.3。 进行平面单元划分之前，先进行合并重节点及压缩节点操作。操作完毕后，平面
2、的线如图所示。 轨道的受力主要在轨头上，故而上面的单元密度应该控制细一点。另外，由于轨道结构比较长（本例采用5跨来分析），二维单元不能太细，否则总的单元数将是极其惊人的。 尝试过用自由网格划分的办法，可是拖拉后，三维单元效果不理想。不是畸变单元太多就是总的单元数太大。最终决定用映射方法进行网格划分。 总体的思路就是将轨道分成三部分：轨头，腹肋和底座。先将现有的面积删除然后重新定义两根线，再由线重新定义面积。同时将L1和L15一分为二！如图所示。 对线进行网格划分控制，采用设置等分分数的方法。设置数据如下： 10 1 1 4 4 4 4 4 10 10 4 2 2 4 4 2 2 16 然后将设
3、置线的关联。 对轨头部分：上面三条线关联起来，下面三条线也关联起来，大概就是上下等分数为12，左右等分数为4。 对腹肋部分：上下等分数为4，左右等分数为10。 对支座部分：上面五条线关联起来，上下等分数为16，左右等分数为2。 设置好后，用映射方法划分网格，如图所示。 接下来是将二维网格拖拉成三维。再做这一步之前，要切记做一件非常非常重要的事情：就是将刚才对线的关联取消！否则拖拉操作将不能成功。将二维网格保存成rail_2Dmesh.db。 由于下面的工作具有很强的规律性，故而采用命令流方法。二维生成三维模型，并对轨枕处施加固定约束的命令流如下： RESUME,rail_2Dmesh,db /
4、prep7 NUMMRG,ALL, , , ,LOW NUMCMP,ALL EXTOPT,ESIZE,30,-0.3 EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,1 REAL,_Z4 ESYS,0 VOFFST,1,568, VOFFST,2,568, VOFFST,3,568, NUMMRG,ALL, , , ,LOW NUMCMP,ALL EXTOPT,ESIZE,20, EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,1 REAL,_Z4 ESYS,0 VOFFST,1,-568, VOFFST,2,-568, VOFFST,3
5、,-568, NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,1136 VGEN,2,5, , ,1136 VGEN,2,6, , ,1136 NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,-568 VGEN,2,5, , ,-568 VGEN,2,6, , ,-568 NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,1704 VGEN,2,5, , ,1704 VGEN,2,6, , ,1704 NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGE
6、N, ,ALL, , ,-284, , , , ,1 NUMMRG,ALL, , , ,LOW NUMCMP,ALL NSEL,S,LOC,X,199,369 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,767,937 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,1335,1420 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,-369,-199 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL N
7、SEL,S,LOC,X,-937,-767 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,-1420,-1335 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL SAVE,rail_3Dmesh_withoutLoad,db FINISH 三维轨道单元模型如图所示： 到此，建模工作已完成。本人认为最重要的部分是在后处理中，本分析的后处理涉及到结果数据的获取及数据文件的输出（不是List命令那样哦，是按照自己想要的格式写到文件中，供其他软件进行进一步数据读写及数据处理）。 时间过得真快，肚子都饿了。暂时先到这里吧
8、！看看有没有人感兴趣，而且也想看看斑竹的意见。免得做一堆无用功，浪费大家的时间。 想要继续的朋友，顶一顶吧。 RESUME,rail_2Dmesh,db ! 重新开始，导入二维网格划分后的模型 /prep7 ! 进入前处理 NUMMRG,ALL, , , ,LOW ! 合并重节点，保留低编号 NUMCMP,ALL ! 压缩编号 EXTOPT,ESIZE,30,-0.3 ! 拖拉设置，在拖拉方向产生30个单元， ! 最小的单元是最大单元的0.3， ! 负号表示中间的单元小，两头的单元大 EXTOPT,ACLEAR,0 ! 拖拉后不删除原来的面 EXTOPT,ATTR,0,0,0 MAT,1 !
9、拖拉后，体的材料特性 REAL,_Z4 ! 拖拉后，体的实常数（默认） ESYS,0 ! 单元坐标系（默认） VOFFST,1,568, ! 将面1向正法线方向拖拉568个单位 VOFFST,2,568, ! 将面2向正法线方向拖拉568个单位 VOFFST,3,568, ! 将面3向正法线方向拖拉568个单位 NUMMRG,ALL, , , ,LOW NUMCMP,ALL EXTOPT,ESIZE,20, ! 拖拉设置，在拖拉方向产生20个单元， ! 各个单元大小相等 EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,1 REAL,_Z4 ESYS,0 VOFFST
10、,1,-568, ! 将面1向负法线方向拖拉568个单位 VOFFST,2,-568, ! 将面2向负法线方向拖拉568个单位 VOFFST,3,-568, ! 将面3向负法线方向拖拉568个单位 NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,1136 ! 体复制，将4号体复制一份到沿正法线方向 ! 1136单位处（以现在位置为起点的相对距离） VGEN,2,5, , ,1136 ! 体复制，将5号体复制一份到沿正法线方向 ! 1136单位处（以现在位置为起点的相对距离） VGEN,2,6, , ,1136 ! 体复制，将6号体复制一份到沿正法线方向
11、 ! 1136单位处（以现在位置为起点的相对距离） NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,-568 ! 体复制，将4号体复制一份到沿负法线方向 ! 568单位处（以现在位置为起点的相对距离） VGEN,2,5, , ,-568 ! 体复制，将5号体复制一份到沿负法线方向 ! 568单位处（以现在位置为起点的相对距离） VGEN,2,6, , ,-568 ! 体复制，将6号体复制一份到沿负法线方向 ! 568单位处（以现在位置为起点的相对距离） NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN,2,4, , ,1704
12、VGEN,2,5, , ,1704 VGEN,2,6, , ,1704 NUMMRG,ALL, , , ,LOW NUMCMP,ALL VGEN, ,ALL, , ,-284, , , , ,1 ! 体平移，将所有体沿负法线方向平移 ! 284个单位（以现在位置为起点的相对距离） NUMMRG,ALL, , , ,LOW NUMCMP,ALL NSEL,S,LOC,X,199,369 ! 轨枕处的约束 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,767,937 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,1335,1420 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,-369,-199 NSEL,R,LOC,Y,-79 D,ALL,ALL ALLSEL,ALL NSEL,S,LOC,X,-937,-767 NSEL,R,LOC,Y,-79 D,ALL,ALL AL
