1、
又参考了一下 ANSYS 帮助文件中的 APDL 编程部分，以及网上的一些意见，
可以有更方
便地得到节点坐标，单元内的节点列表的方法。在一个文本文件中写几条命令
就可以
了。而且可以指定输入写入的格式，不会出现两个数据首尾相接，在 MATLAB
中无法区
分的问题。
得 到节点坐标和单元内包含的节点的程序如下，可以在 ANSYS 的帮助中找到
每一个命
令的详细解释。将下面的程序拷贝到一个文本文件中 lename.txt，保存，去
掉后
缀名 TXT，再拷贝到工作目录下。在命令提示符下输入*use lename，生产
一个 geomle
.txt 文件，打开可以看到里面的数据
! Get the coordinates of each node
! Get the node list of each element
! By:LiuXiaoqin(liuxqsmile@gmail.com), at nuaa,
*get,nodenum,node,,num,max ! 获得节点的数目
*dim,nodepos,array,nodenum,3 ! nodepos 存放节点的坐标
*do,i,1,nodenum,1
*get,nodepos(i,1),node,i,loc,x !获得节点的 X 坐标
*get,nodepos(i,2),node,i,loc,y
*get,nodepos(i,3),node,i,loc,z
*enddo
!以上可以使用 nx,ny,nz 代替
*get,elemnum,elem,,num,max !得到单元的总数目
*dim,elemlist,array,elemnum,6 !单元包含的节点列表，指定每个单元包含
6 个节点
，根据情况修改
*do,i,1,elemnum,1
*do,ii,1,6,1
*get,elemlist(i,ii),elem,i,node,ii !获得节点编号
*enddo
*enddo
*cfopen,geomle,txt !打开文件，写入数据
*vwrite,0
(F8.0,' Coordinates of each node')
*vwrite,sequ,nodepos(1,1),nodepos(1,2),nodepos(1,3)
(F8.0,3e16.8)
*vwrite,0
(F8.0,' Nodes No. of each element')
*vwrite,sequ,elemlist(1,1),elemlist(1,2),elemlist(1,3),elemlist(1,4),el
评论0