有大佬能帮忙看看这个焊接空心球的ANSYS命令流做的对吗？
浏览：405
finish
/clear
/prep7
! 参数定义
c=1 ! 套管厚度系数
cc=1.4 ! 套管直径系数
ccc=0.3 ! 套管高度系数（从球顶点开始算起）
T=127 ! 钢管直径
ts=10 ! 球壁厚
D=350 ! 球径
tt=ts+10 ! 钢管壁厚
tr=c*ts ! 套管厚（10）
T1=cc*T ! 套管直径（177.8）
Hr=ccc*D ! 套管高（105）
Ht=1.5*D ! 钢管高
! 材料属性定义
mp,dens,1,7.85e-6 ! kg/mm3
mp,ex,1,201433 ! N/mm2
mp,prxy,1,0.3
TB,miso,1,1,4,0
TBTEMP,0
TBPT,,360/201433,360
TBPT,,0.02,360
TBPT,,0.14,503.7
TBPT,,0.16,503.7 ! 限制屈服应力
! 几何建模
sphere,D/2-ts,D/2
wprota,,270
cylind,,T1/2,,D/2+Hr ! 连接板
cylind,T/2-tt,T/2,,Ht ! 钢管
cylind,T1/2-tr,T1/2,,D/2+Hr ! 套管
vptn,all
wpcsys,-1
wprota,,90
vsbw,all
wpcsys,-1
csys
vsel,s,loc,y,-D/2,0
vdele,all,,,1
allsel
! 进一步分割体
vsbw,all
wprota,,,45
vsbw,all
wpcsys,-1
csys,5 ! 柱坐标系
vsel,u,loc,y,0,45
vdele,all,,,1
allsel
csys,2
vsel,s,loc,x,0,D/2-3*ts
vdele,all,,,1
allsel
csys
wpoff,,,-tr/2
vsbw,all
wpcsys,-1
! 选择并删除不需要的体
csys,2
vsel,u,loc,x,0,D/2
csys,5
vsel,r,loc,x,0,T/2-tt
vdele,all,,,1
allsel
csys,2
vsel,u,loc,x,0,D/2
csys,5
vsel,u,loc,x,0,T/2
csys
vsel,u,loc,z,-tr/2,0
vdele,19,,,1
allsel
! 合并和压缩编号
nummrg,all
numcmp,all
! 分配材料属性
csys
vsel,s,loc,y,D/2,Ht
csys,5
vsel,r,loc,x,0,T/2
vatt,1
vsel,inve
vatt,1
allsel
vplot,all
csys
wpcsys,-1
! 网格划分设置
vsel,s,loc,z,-tr/2,0
csys,2
vsel,r,loc,x,0,D/2
csys,5
vsel,r,loc,x,T/2-tt,T/2
csys
/pnum,area,1
/pnum,line,1
aslv,s,1
asel,r,loc,z,0
lesize,74,,,4
! 划分面网格
asel,s,,,55
et,1,200,7
esize,6
mshape,0,2d ! 四边形单元
mshkey,1
mshmid,0
amesh,all
allsel
! 划分体网格
vsel,s,loc,z,-tr/2,0
csys,2
vsel,r,loc,x,0,D/2
csys,5
vsel,r,loc,x,T/2-tt,T/2
csys
!vsel,s,,,4
et,2,solid45
mshape,0,3d ! 六面体单元
mshkey,1
mshmid,0
vsweep,all
vsel,inve
vsweep,all
allsel
! 对称操作
vsel,s,loc,z,-tr/2,0
csys,2
vsel,u,loc,x,0,D/2
csys,5
vsel,r,loc,x,T/2,T1/2-tr
vsel,inve
csys
wpcsys,-1
wprota,,,45
csys,4
vsymm,z,all
allsel
nummrg,all
wprota,,,315
allsel
vplot,all
vsymm,z,all
vsymm,y,all
vsymm,x,all
nummrg,all
! 施加约束
csys
wpcsys,-1
asel,s,loc,y,-525
da,all,all
allsel
finish
! 求解设置
/config,nproc,2 ! 推荐使用2核
/solu
!eqslv,amg
eqslv,sparse ! 对于17万单元，稀疏求解器更稳定
antype,0
nropt,auto ! N-R法
nlgeom,on ! 考虑几何非线性
cnvtol,f,100,0.01
!cnvtol,u,1
nsubst,100
outres,all,all
asel,s,loc,y,Ht
da,all,uy,-0.05*D ! 顶部施加向下的位移荷载
allsel
vplot,all
acel,,9.8 ! 重力加速度
allsel
eresx,no
solve
finish
! 后处理
/post26
nsel,s,loc,y,Ht
*get,ncount,node,0,count ! 获得该面上的节点总数
*SET,nNow,0 ! 定义变量nNow
nsol,6,ndnext(nNow),u,y,ydisplacement ! 节点竖向位移（mm）
add,6,6,,,,,,-1 ! 对位移取绝对值
*do,i,1,ncount
*SET,nNow,ndnext(nNow)
rforce,8,nNow,f,y,yforce
add,7,7,8,,,,,,-0.001
*enddo
prod,7,7,,,,,,1
xvar,6
/axlab,x,displacement(mm)
/axlab,y,force(kN)
plvar,7 !绘制变量
prvar,6,7 !
! 显示等效塑性应变
!set,last
!plnsol,eppl,eqv
