/filename,ex5-1 !设定工作名
/title,The Transient Analysis of the pipe !设定工作标题
/prep7 !进入前处理模块
*set,dia,139.7 !定义变量参数
*set,thick,7.72
*set,hole_dia,12.7
*set,hole_length,25.4
*set,hole_num,10
*set,hole_ang,18
!
et,1,solid95 !选择单元类型SOLID95
mp,ex,1,2.1e5 !定义弹性模量
mp,prxy,1,0.3 !定义泊松比
mp,dens,1,0.00785 !定义密度
wpoff,0,0,-hole_length !沿-z方向偏移工作平面hole_length
cyl4,0,0,dia/2-thick,,dia/2,,hole_length !生成一个圆柱体
!CYL4,XCENTER,YCENTER,RAD1,THETA1,RAD2,THETA2,DEPTH
wpoff,0,0,hole_length/2 !沿+z方向偏移工作平面hole_length
wpro,,,90 !绕y轴旋转工作平面90度
cyl4,0,0,0,,hole_dia/2,,100 !生成圆柱体
!
csys,1 !激活柱坐标系
*do,i,1,hole_num !do循环开始
vgen,2,2,,,,i*360/hole_num,0 !循环复制生成周向的圆柱体
!VGEN,ITIME,NV1,NV2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE
*enddo !do循环结束
!
vovlap,all !对所有的体进行搭接操作
!将分离的同等级对象，变为数个对象的连续体，重叠边界所围成的区域自成一个对象。
allsel,all !选择全部图元
/replot
*get,max_number,volu,0,num,max !将体图元的最大编号置于参数max_number中
!*GET,Par,Entity,ENTNUM,Item1,IT1NUM,Item2,IT2NUM
vdele,1,max_number-1,1,1 !删除编号为1~max_number-1的体图元
!VDELE,NV1,NV2,NINC,KSWP
numcmp,all !压缩编号
!
vgen,5,1,,,,-hole_ang,-hole_length,0 !复制生成轴向的圆柱体
vglue,all !粘接所有的体
!将不连续的同等级对象变为连续体，仍为多个对象。vadd，all为变成一个对象。
/replot
!mesh the model
lsel,s,loc,z,0 !选择z=0的线
lsel,a,loc,z,-5*hole_length !同时选中z=-5*hole_length的线
lesize,all,,,10 !设定线的分割数10
smrtsize,8 !设定智能划分等级8
allsel !选择全部图元
mshkey,0 !设定自由网格划分方式
mshape,1,3d !设定为四面体网格
vmesh,all !对体进行网格划分
!node component
nsel,s,loc,z,0 !选择z=0的节点
cm,zuoduanmian,node !生成节点组件
nsel,s,loc,z,-5*hole_length !选择z=-5*hole_length的节点
cm,youduanmian,node !生成节点组件
!
!apply the displacement
nsel,s,loc,z,-hole_length !选择z=-hole_length的节点
nsel,a,loc,z,-4*hole_length !同时选中z=-4*hole_length的节点
