ANSYS影响面计算与绘制方法介绍

finish
/clear
/prep7
et,1,shell181
mp,ex,1,3.0e4
mp,prxy,1,0.2
mp,dens,1,2500e-12
sectype,1,shell
secdata,100
blc4,,,3000,1000
agen,4,all,,,3000
aglue,all
kwpave,2
wprota,,,90
blc4,,,2000,1000
wpcsys,-1
asel,s,,,2
agen,3,all,,,3000
allsel,all
aglue,all
!=================
esize,300
amesh,all
!==================
nsel,s,loc,x,0
nsel,a,loc,x,3000*4
d,all,ux,0,,,,uz,uy
nsel,s,loc,z,-2000
d,all,all,0
!==================
allsel,all
!获取最大节点数目       
nodemax=ndinqr(0,13)   
!==================
!获取加载面的最大节点数目
nsel,s,loc,z,0
nmax2=ndinqr(0,13)
*get,n0,node,,num,min
*dim,JZnode,,nmax2
JZnode(1)=n0
*do,i,2,nmax2
jZnode(i)=ndnext(JZnode(i-1))
*enddo
finish
!==================
!加载求解
/solu
allsel,all
outres,all,all
*do,i,1,nmax2
fdele,all,all
f,jznode(i),fz,-1
solve
*enddo
finish
!====================
/post1
!节点202的UZ影响面
!节点42的UZ影响面
!单元135的MX影响面
*dim,Uz1,,nmax2
*dim,UZ2,,nmax2
*dim,MXX,,nmax2
*do,i,1,nmax2
set,i
Uz1(i)=uz(202)
UZ2(i)=uz(42)
*get,MXX(i),elem,135,smisc,4
*enddo
/graphics
*do,i,1,nodemax
dnsol,i,u,z,0
*enddo
!绘制202 UZ影响面
*do,i,1,nmax2
dnsol,jznode(i),u,z,uz1(i)
*enddo
plnsol,u,z
!绘制42 UZ影响面
*do,i,1,nmax2
dnsol,jznode(i),u,z,uz2(i)
*enddo
plnsol,u,z
!绘制单元135 MX影响面
*do,i,1,nmax2
dnsol,jznode(i),u,z,mxx(i)
*enddo
plnsol,u,z
/device,vector,1
plnsol,u,z