ANSYS套箱围堰计算代码.DOC

/PREP7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建模一 !*单位cm,kgf !*建立模型 ET,1,LINK1 ET,2,BEAM3,,,,,,1 !*在单元坐标系中输出力和力矩 KEYOPT,2,10,0 !*和SFBEAM配合,指定以长度为单位,荷载相对于I,J ET,3,BEAM4,,,,,,1 KEYOPT,3,10,0 !*和SFBEAM配合,指定以长度为单位,荷载相对于I,J ET,4,SHELL63,1,,0,,0 !*仅考虑弯曲刚度 MP,EX,1,2.06E6 MP,PRXY,1,0.3 !*POISSONS RATIO IS 0.3 MP,DENS,1,7.850E-3 NOD1_R=1 REL1_R=1 ELE1_R=1 NUMSTR,NODE,NOD1_R NUMSTR,ELEM,ELE1_R R,REL1_R+1,25.48, 396.6,283.82,12.6,10 !*顺桥向模板,竖肋,2[10 R,REL1_R+2,25.48,283.3,396.6,10,12.6 !*横桥向模板,竖肋,2[10 R,REL1_R+3,60.97,7481,363.9,12.4,28 !*模板,套箍,I28b R,REL1_R+4,0.6 !*THICKNESS = 0.4cm !*用实常数显示单元 /ESHAPE,1 !*顺桥向模板 K,1,,,120 K,2,,,160 K,3,,,280 K,4,,,430 K,5,,,530 K,6,,,740 KGEN,2,1,6,1,1288.8 !*定义5个面 A,1,7,8,2 A,2,8,9,3 A,3,9,10,4 A,4,10,11,5 A,5,11,12,6 !*分配面上的单元属性 ASEL,S,AREA,,1,5,1 AATT,1,REL1_R+4,4,0 !*指定单元形状为四边形 MSHAPE,0 !*指定映射网格划分 MSHKEY,1 !*打开单元形状检查 SHPP,ON LESIZE,1,,,36 LESIZE,2,,,4 LESIZE,5,,,4 LESIZE,8,,,6 LESIZE,11,,,4 LESIZE,14,,,6 AMESH,ALL ASEL,ALL !*建立竖肋 LGEN,2,4,16,3,35.8 LGEN,18,17,21,1,71.6 LSEL,S,LINE,,17,106,1 LATT,1,REL1_R+1,3 LMESH,ALL LSEL,ALL !*建立套箍 LSEL,S,LINE,,3,15,3 LATT,1,REL1_R+3,3 LMESH,ALL LSEL,ALL *GET,NELEM1_R,ELEM,0,COUNT !*获取单元总数 *GET,NNODE1_R,NODE,0,COUNT !*获取节点总数 EGEN,2,NNODE1_R+NOD1_R-1,ELE1_R,ELE1_R+NELEM1_R-1,1,,,,,,,1747.2 !*横桥向模板 KGEN,2,1,6,1,,1747.2 !*定义5个面 A,1,121,122,2 A,2,122,123,3 A,3,123,124,4 A,4,124,125,5 A,5,125,126,6 !*分配面上的单元属性 ASEL,S,AREA,,6,10,1 AATT,1,REL1_R+4,4,0 !*指定单元形状为四边形 MSHAPE,0 !*指定映射网格划分 MSHKEY,1 !*打开单元形状检查 SHPP,ON LESIZE,107,,,48 AMESH,ALL ASEL,ALL !*建立竖肋 LGEN,2,4,16,3,,36.4 LGEN,24,118,122,1,,72.8 LSEL,S,LINE,,118,237,1 LATT,1,REL1_R+2,3 LMESH,ALL LSEL,ALL !*建立套箍 LSEL,S,LINE,,109,117,2 LATT,1,REL1_R+3,3 LMESH,ALL LSEL,ALL *GET,NELEM2_R,ELEM,0,COUNT !*获取单元总数 *GET,NNODE2_R,NODE,0,COUNT !*获取节点总数 EGEN,2,NNODE2_R+NOD