关于ansys命令流的网格划分
浏览：93890 回答：1
之前划分的是一个有缺陷（通孔）的管道，差不多如下图所示，如果我只是想建立一个小凹槽的缺陷，该如何修改我的命令流？
[img]file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/258910249/QQ/WinTemp/RichOle/4ZYCNLT$V6[~5AYQ1R]_4SG.jpg[/img]
/prep7
/title,001
nl_crack=1 !裂纹轴向长度占几个轴向单元长度，即宽
nc_crack=3 !裂纹周向长度占几个周向单元，即长
nr_crack=1/3 !裂纹径向长度占管道壁厚的比例，即高
nsin=10 !正弦激励周期数
!
!
le=10 !单元长度,轴向
l1=0.5e3 !缺陷在管道轴向的 开始位置
l3=0.5e3 !从缺陷开始左端完整管道长度
lcrack=nl_crack*le !缺陷长度
l2=l1+lcrack !缺陷在管道轴向的 结束位置
l=l1+l3 !管道长度
ri=104.5 !内径
ro=109.5 !外径
h=5 !壁厚
nv1=24 !周向划分的份数，也即一个圆周上节点的个数。
letheta=360/nv1 !管道周向划分，单元角度数
nr=3 !径向划分份数
theta=nc_crack*letheta !缺陷周向角度
thick=nr_crack*h !裂纹深度,穿孔
h1=ro-thick
!
et,1,45
mp,ex,1,2.1e5 !弹性模量
mp,nuxy,1,0.28 !泊松比
mp,dense,1,7.932e-8 !密度
!!*****************************建模开始**************************************
CYL4,0,0,ri,0,Ro,180,l1 !建激励端上半个模型
WPOFFS,0,0,l2
cyl4,0,0,ri,0,Ro,180,l3-lcrack !建远端上半个模型
vsymm,y,all !镜像复制
cyl4,0,0,ri,0,h1,360,-lcrack
cyl4,0,0,h1,theta,ro,360,-lcrack
wpcsys,,0 !转换工作平面到0坐标系
vglue,all !对体进行粘接（GLUE）
nummrg,all !合并相同的项
numcmp,all
/PNUM,VOLUM,1 !/PNUM, Label, KEY 打开控制实体编号显示
!!*****************************建模结束**************************************
!!*****************************网格划分开始**********************************
!
lsel,s,loc,x,ri !根据坐标选择一组线子集，轴向划分，单元长le
lsel,a,loc,x,ro !同时根据坐标选择一组线子集
lsel,a,loc,x,-ri
lsel,a,loc,x,-ro
lesize,all,le,, !对所选择的线段设置网格单元
csys,1
lsel,s,loc,x,ri
lsel,a,loc,x,ro
lsel,r,loc,y,theta !选择裂纹轴向轮廓
lesize,all,le,,
lsel,s,loc,x,ri !周向划分
lsel,a,loc,x,ro
lsel,u,loc,y,0
lsel,u,loc,y,180
lsel,u,loc,y,theta
lesize,all,,letheta,, !以15度为单元长度 来划分
csys,0 !径向划分
lsel,s,loc,x,ri,ro
lsel,u,loc,x,ri
lsel,u,loc,x,ro
csys,1
lsel,u,loc,y,1,theta !选择介于1度到theta度以外的角度
lesize,all,,,nr !径向每条线段的分段数为nr
csys,0
lsel,s,loc,x,-ro,-ri
lsel,u,loc,x,-ri
lsel,u,loc,x,-ro
lesize,all,,,nr
csys,1 !裂纹处的上半部分管道径向划分，x为轴向，y为周向，z为径向
lsel,s,loc,z,l1
lsel,a,loc,z,l2
lsel,r,loc,y,theta-letheta/5,theta+letheta/5 !?
lesize,all,,,nr
vsel,all !选择所有的体
vatt,1,1,1 !给所选择的体设置网格划分单元属性为材料号、实常数设置号、单元类型号为1
MSHAPE,0,3D !指定生成六面体单元
allsel,all !选择所有实体
vsweep,all
csys,1
NROTAT,all !转换节点坐标系到柱坐标系
allsel,all
save,modeled,db,,
[img]file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/258910249/QQ/WinTemp/RichOle/4ZYCNLT$V6[~5AYQ1R]_4SG.jpg[/img]
/prep7
/title,001
nl_crack=1 !裂纹轴向长度占几个轴向单元长度，即宽
nc_crack=3 !裂纹周向长度占几个周向单元，即长
nr_crack=1/3 !裂纹径向长度占管道壁厚的比例，即高
nsin=10 !正弦激励周期数
!
!
le=10 !单元长度,轴向
l1=0.5e3 !缺陷在管道轴向的 开始位置
l3=0.5e3 !从缺陷开始左端完整管道长度
lcrack=nl_crack*le !缺陷长度
l2=l1+lcrack !缺陷在管道轴向的 结束位置
l=l1+l3 !管道长度
ri=104.5 !内径
ro=109.5 !外径
h=5 !壁厚
nv1=24 !周向划分的份数，也即一个圆周上节点的个数。
letheta=360/nv1 !管道周向划分，单元角度数
nr=3 !径向划分份数
theta=nc_crack*letheta !缺陷周向角度
thick=nr_crack*h !裂纹深度,穿孔
h1=ro-thick
!
et,1,45
mp,ex,1,2.1e5 !弹性模量
mp,nuxy,1,0.28 !泊松比
mp,dense,1,7.932e-8 !密度
!!*****************************建模开始**************************************
CYL4,0,0,ri,0,Ro,180,l1 !建激励端上半个模型
WPOFFS,0,0,l2
cyl4,0,0,ri,0,Ro,180,l3-lcrack !建远端上半个模型
vsymm,y,all !镜像复制
cyl4,0,0,ri,0,h1,360,-lcrack
cyl4,0,0,h1,theta,ro,360,-lcrack
wpcsys,,0 !转换工作平面到0坐标系
vglue,all !对体进行粘接（GLUE）
nummrg,all !合并相同的项
numcmp,all
/PNUM,VOLUM,1 !/PNUM, Label, KEY 打开控制实体编号显示
!!*****************************建模结束**************************************
!!*****************************网格划分开始**********************************
!
lsel,s,loc,x,ri !根据坐标选择一组线子集，轴向划分，单元长le
lsel,a,loc,x,ro !同时根据坐标选择一组线子集
lsel,a,loc,x,-ri
lsel,a,loc,x,-ro
lesize,all,le,, !对所选择的线段设置网格单元
csys,1
lsel,s,loc,x,ri
lsel,a,loc,x,ro
lsel,r,loc,y,theta !选择裂纹轴向轮廓
lesize,all,le,,
lsel,s,loc,x,ri !周向划分
lsel,a,loc,x,ro
lsel,u,loc,y,0
lsel,u,loc,y,180
lsel,u,loc,y,theta
lesize,all,,letheta,, !以15度为单元长度 来划分
csys,0 !径向划分
lsel,s,loc,x,ri,ro
lsel,u,loc,x,ri
lsel,u,loc,x,ro
csys,1
lsel,u,loc,y,1,theta !选择介于1度到theta度以外的角度
lesize,all,,,nr !径向每条线段的分段数为nr
csys,0
lsel,s,loc,x,-ro,-ri
lsel,u,loc,x,-ri
lsel,u,loc,x,-ro
lesize,all,,,nr
csys,1 !裂纹处的上半部分管道径向划分，x为轴向，y为周向，z为径向
lsel,s,loc,z,l1
lsel,a,loc,z,l2
lsel,r,loc,y,theta-letheta/5,theta+letheta/5 !?
lesize,all,,,nr
vsel,all !选择所有的体
vatt,1,1,1 !给所选择的体设置网格划分单元属性为材料号、实常数设置号、单元类型号为1
MSHAPE,0,3D !指定生成六面体单元
allsel,all !选择所有实体
vsweep,all
csys,1
NROTAT,all !转换节点坐标系到柱坐标系
allsel,all
save,modeled,db,,
ansys_learner 
