命令流:
/prep7
et,1,63
mp,ex,1,700000
mp,prxy,1,0.33
k,1,0,0,0
k,2,0,500,0
k,3,100,500,0
k,4,100,0,0
a,1,2,3,4 !area1
k,5,220,500,15
k,6,220,0,15
a,3,4,6,5 !area2
k,,290,500,15
k,,290,0,15
a,5,6,8,7 !area3
k,,390,500,15
k,,390,0,15
a,7,8,10,9 !area4
k,,460,500,15
k,,460,0,15
a,9,10,12,11 !area5
asel,s,,,2
type,1
mat,1
esize,25
amesh,all
*get,a2_min,node,,num,min !计算该面上结点的最小编号
*get,a2_max,node,,num,max !计算该面上结点的最大编号
*set,n_a2,a2_max-a2_min+1 !计算该面上的结点总数
allsel
asel,s,,,4
type,1
mat,1
esize,25
amesh,all
allsel
asel,s,,,4
cm,a_4,area
allsel
cmsel,s,a_4,area
nsla,r,1
*get,a4_min,node,,num,min
*get,a4_max,node,,num,max
*set,n_a4,a4_max-a4_min+1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!! area2厚度随x坐标的变化方程 !!!!!
!!!!!! y=-.25*x+75.0 !!!!!
!!!!!! area4厚度随x坐标的变化方程 !!!!!
!!!!!! y=.4*x-96.0 !!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*dim,h_num,,n_a2+n_a4
*do,i,a2_min,a2_max,1
h_num(i)=-.25*nx(i)+75.0
*enddo
*do,i,a4_min,a4_max,1
h_num(i)=.4*nx(i)-96.0
*enddo
rthick,h_num
allsel
*get,emax,elem,,num,max !获取单元总数
r,emax+1,50
r,emax+2,20
r,emax+3,60
allsel
asel,s,,,1
type,1
mat,1
real,emax+1
amesh,all
allsel
asel,s,,,3
type,1
mat,1
real,emax+2
amesh,all
allsel
asel,s,,,5
type,1
mat,1
real,emax+3
amesh,all
ansys图: |