管道纵向模态频散曲线及波形分析.zip

FINISH /CLEAR,START /config,nres,10000 ! !最多结果子步数 !设置材料参数--以mm为单位建模，所以材料属性要换算单位 !材料属性-钢筋(弹性模量200Gpa，密度7800，泊松比0.30） /PREP7 ET,1,42 et,2,45 MP,EX,1,200E9 MP,DENS,1,7800 MP,PRXY,1,0.3 R1=0.006 !半径 R2=0.0054 L=1.2 L1=0.6 L2=0.1 A_Z=0.003 !缺陷沿z轴方向的长度 HZ=100000 HZN=5 W=3.14159265357 PI=ACOS(-1) !π的计算=3.1415926 !划分单元大小 DA_Z=0.003 !Z方向每1mm划分一个单元 !划分一个小部分 CYL4,0,0,0,0,R2,90 CYL4,0,0,R2,0,R1,15 csys,1 AGEN,6,2,,,,15,0 !面积复制命令,包含本身所复制的次数,为现有的坐标系统下复制到其他位置,,, APTN,ALL !面分割 NUMMRG,ALL NUMCMP,ALL csys,0 ARSYM,X,ALL !复制一组面积 ,对称于轴, ALLSEL ARSYM,Y,ALL APTN,ALL NUMMRG,ALL NUMCMP,ALL ALLSEL CM,AREA1,AREA !定义组元，将几何元素分组形成组元 , 由字母数字组成的组元名 , 组元的类型（volu, area, line, kp, elem, node SAVE lsel,s,LENGTH,,R2 LESIZE,ALL,,,6 !为线指定网格尺寸 ALLSEL LSEL,U,LENGTH,,R2 LESIZE,ALL,,,1 ALLSEL SAVE vext,AREA1,,,0,0,L1 !通过给定偏移量由面生成体 allsel asel,s,loc,z,l1-0.001,l1+0.001 cm,area2,area vext,area2,,,0,0,A_Z allsel asel,s,loc,z,l1+A_Z-0.001,l1+A_Z+0.001 cm,area3,area vext,area3,,,0,0,L2 allsel asel,s,loc,z,l1+A_Z+l2-0.001,l1+A_Z+l2+0.001 cm,area4,area vext,area4,,,0,0,A_Z allsel asel,s,loc,z,L1+L2+2*A_Z-0.001,L1+L2+2*A_Z+0.001 cm,area5,area vext,area5,,,0,0,l-l1-l2-2*A_Z ALLSEL VPTN,ALL NUMMRG,ALL NUMCMP,ALL SAVE !划分单元 ASEL,S,AREA,,AREA1 MSHAPE,0,2D MSHKEY,1 TYPE,1 AMAP,7,1,2,3 AMAP,26,19,2,3 AMAP,28,19,39,3 AMAP,27,1,39,3 ASEL,U,AREA,,7 ASEL,U,AREA,,26,28,1 AMESH,ALL ALLSEL NUMMRG,ALL NUMCMP,ALL LSEL,S,LOC,X,0 LSEL,R,LOC,Y,0 LESIZE,ALL,DA_Z,, ALLSEL MSHAPE,0,3D MSHKEY,1 TYPE,2 VSWEEP,ALL SAVE !删除一些单元形成缺陷 VSEL,S,LOC,Z,L1,L1+A_Z vsel,u,volu,,35 vsel,u,volu,,54,56,1 VCLEAR,all VDELE,all, ALLSEL SAVE ! 加载计算 /solu ANTYPE,TRANS ! 静力分析, 瞬态分析 TRNOPT,FULL ! 指定瞬态分析选项 tintp,,,,1.0 ! 定义短暂集成参数 OUTRES,NSOL,ALL ! 规定写入数据库的求解信息,节点自由度,每一步都写 NSEL,S,LOC,Z,l D,ALL,ALL,0 ALLSEL csys,1 nsel,s,loc,z,0 nsel,r,loc,x,r1 CM,JIAZAIJIEdian,node ALLSEL CSYS,0 SAVE TIMINT,ON AUTOTS,ON ! 设定自动时间步长 ! 5个周期，共20微妙。设为125个载荷步 *DO,I,1,125,1 TIME,(4E-7)*I DELTIM,4E-7,4E-7,4E-7 KBC,1 D,JIAZAIJIEdian,UZ,(0.5*(1-cos(2*W*HZ*(4E-7)*I/HZN))*sin(2*W*HZ*(4E-7)*I))*0.0000005 LSWRITE,I *ENDDO TIME,0.8E-3 DELTIM,4E-7,4E-7,4E-7 ddele,JIAZAIJIEdian,uZ LSWRITE,126 ALLSEL SAVE LSSOLVE,1,126 SAVE FINISH