ANSYS命令流（APDL）源代码：简支梁桥梁的建立（beam4、beam188的形变及自振频率）（耦合&弹簧连接梁墩）

finish /clear /CONFIG,NPROC,2 /prep7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!桥梁!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! local,11,0,,,,0,90,0 !桥梁习惯局部坐标系 CSYS,11 et,1,beam188 !定义单元属性 Bridge_E=3.5E10 !桥梁弹性模量 Bridge_DENS=2.5E3 !桥梁密度 MP,Ex,1,Bridge_E ! 材料属性 MP,DENS,1,Bridge_DENS MP,PRXY,1,0.3 SECTYPE,1,BEAM,MESH SECREAD,mysec,SECT,,MESH !读取截面 secoffset,cent k,1,,,-1.1 k,2,32,,-1.1 k,3,,,-10 l,1,2 TYPE,1 !指定单元类型 MAT,1 latt,1,,1,,3,,1 lesize,all,,,64 lmesh,all !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!桥墩!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! local,12,0,,-10,,0,-90,0 !桥梁习惯局部坐标系 CSYS,12 et,1,beam188 !定义单元属性 !混凝土材料属性 pier_L=10 !桥墩长度 pier_NUM=2 !桥墩数量 pier_EL=0.2 !单元长度 pier_NE=pier_L/pier_EL !每座桥墩单元数量 pier_H=1.1 !桥梁中心距箱梁底部的高度，通过ASUM查质心的Y坐标 pier_EX=3.55e10 !弹性模量 pier_DENS=2500 !密度+加二期恒载等效密度 pier_Prxy=0.2 MP,EX,2,pier_Ex MP,PRXY,2,pier_Prxy MP,DENS,2,pier_DENS SECTYPE,2,BEAM,MESH, !截面类型为beam截面，该截面编号为2； SECOFFSET,CENT,,, !截面形心位于中间 SECREAD,dunbox,SECT,,MESH !读取截面 k,4, k,5,0,0,10 l,4,5 k,6,1,0,0 K,7,32 K,8,32,0,10 K,9,33,0,0 L,7,8 TYPE,1 !指定单元类型 MAT,2 latt,2,,1,,6,,2 !!!!!!最后一个1还要改成2 lesize,all,,,50 lmesh,all !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!支座弹簧!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ET,3,COMBIN14 KEYOPT,3,2,1 !定义x方向弹簧 ET,4,COMBIN14 KEYOPT,4,2,2 !定义Y向弹簧 ET,5,COMBIN14 KEYOPT,5,2,3 !定义Z向弹簧 ET,6,COMBIN14 KEYOPT,6,2,4 !定义rotx方向弹簧 !定义弹簧单元实常数 Mac_Spring_EX1=1E16 Mac_Spring_EY1=1E16 Mac_Spring_EZ1=1E16 Mac_Spring_ErotX1=1E16 Mac_Spring_CX1=1E16 Mac_Spring_CY1=1E16 Mac_Spring_CZ1=1E16 Mac_Spring_CrotX1=1E16 R,1,Mac_Spring_EX1,Mac_Spring_CX1 R,2,Mac_Spring_EY1,Mac_Spring_CY1 R,3,Mac_Spring_EZ1,Mac_Spring_CZ1 R,4,Mac_Spring_ErotX1,Mac_Spring_CrotX1 TYPE,3 REAL,1 E,1,131 E,2,232 TYPE,4 REAL,2 E,1,131 E,2,232 TYPE,5 REAL,3 E,1,131 E,2,232 TYPE,6 REAL,4 E,1,131 E,2,232 CSYS,0 NSEL,S,LOC,Y,-10 D,all,all ALLSEL,ALL ACEL,,9.806 !自重应力 /SOLU ANTYPE,0, !定义静态分析类型 AUTOTS,ON OUTRES,ALL,all !输入控制 NLGEOM,ON !大变形效应 PSTRES,ON !预应力效应 EQSLV,PCG,1.0E-5, !求解器选择 TIME,1 !第一步应力 NSUBST,5,10,5 !子步数设置 allsel,all SOLVE /triad,off !不显示整体坐标系 /eshape,1 eplot /post1 PLDISP,u,sum !显示总位移结果 /dscale,all,100 !显示位移放大系数