由有限元网格模型得到几何模型

从四面体单元再产生几何实体的示例命令流（如果要转换六面体网格，可以通过修改下面的命令流实现。）：
1. 转换未变形的网格：
esel,s,mat,,1 !假设单元材料号为1，先选择单元
*get,enum1,elem,0,count
*get,el,elem,0,num,min
ksel,all
*get,nkp1,kp,0,num,max
ksel,none
asel,none
vsel,none
/prep7
*do,i,1,enum1
*do,j,1,3
nd=nelem(el,j)
knum1=nkp1+4*(i-1)+j
k,knum1,nx(nd),ny(nd),nz(nd) !
*enddo
nd=nelem(el,j+2)
k,knum1+1,nx(nd),ny(nd),nz(nd)
v,knum1-2,knum1-1,knum1,knum1+1
el=elnext(el)
*enddo
numm,kp
cm,aal1,area
asel,r,ext
cm,aext,area
vdel,all
va,all
cmsel,s,aal1
cmsel,u,aext
adel,all,,,1

2. 如果需要转换变形后的网格：
esel,s,mat,,1 !假设单元材料号为1，先选择单元
*get,enum1,elem,0,count
*get,el,elem,0,num,min
ksel,all
*get,nkp1,kp,0,num,max
ksel,none 给13亿中国人们更多听觉关怀！！
asel,none
vsel,none
/post1
set,last !输出最后一步变形后的网格
/prep7
SHPP,OFF
UPGEOM,1,LAST,LAST,'filname','rst',' ' !假设工作文件名为filname
*do,i,1,enum1
*do,j,1,3
nd=nelem(el,j)
knum1=nkp1+4*(i-1)+j
k,knum1,nx(nd),ny(nd),nz(nd) !
*enddo
nd=nelem(el,j+2)
k,knum1+1,nx(nd),ny(nd),nz(nd)
v,knum1-2,knum1-1,knum1,knum1+1
el=elnext(el)
*enddo
numm,kp
cm,aal1,area
asel,r,ext
cm,aext,area
vdel,all
va,all
cmsel,s,aal1
cmsel,u,aext
adel,all,,,1