ANSYS apdl命令流笔记2-------单元排序

前言

通过do命令与if命令的配合使用，来对图中面单元进行排序。

一、划分网格

finish
/clear
/prep7
blc4,,,9,9  !建立一个四边形
et,1,13		!单元类型为plane13
esize,,3	!网格尺寸
amesh,all

如图所示

示例：pandas 是基于NumPy 的一种工具，该工具是为了解决数据分析任务而创建的。

二、单元排序

1.引入库

allsel,all
*get,emax,elem,,num,max		!获取最大单元编号
*get,d,elem,,count			!获取单元数
*dim,ne,,d
*dim,ex,,d
*dim,ey,,d			
*dim,eorder,,d		!建立四个维数为d的数组
mine=0				!定义一个参数
*do,j,1,d
	esel,u,elem,,mine	!每次选择时减去编号为mine的单元
	*get,c,elem,,count	!重新获取新选择区域的单元数
	ii=0			!设置一个参数，用于计算单元个数
	*do,I,1,emax
		*if,esel(i),eq,1,then	!判断单元是否被选择
			ii=ii+1
			ne(ii)=i	!将单元编号存于数组
		*endif
	*enddo	!结束循环，具体可参考上一篇笔记
	*do,I,1,c
		*get,ey(i),elem,ne(i),cent,y !得到第i号单元y坐标
		*get,ex(i),elem,ne(i),cent,x !得到第i号单元x坐标
	*enddo
	miny=1e10
	minx=1e10	!定义两个较大的数
	*do,I,1,c
		*if,ex(i),lt,minx,then 	!判断单元x坐标是否小于minx
			miny=ey(i)
			minx=ex(i)
			mine=ne(i)
		*else
		*if,ex(i),eq,minx,then
		*if,ey(i),lt,miny,then	!判断单元y坐标是否小于miny
			miny=ey(i)
			minx=ex(i)
			mine=ne(i)
		*endif
		*endif
		*endif
	*enddo		!结束循环
	eorder(j)=mine	!将判断后的坐标值最小单元号存于数组
*enddo	

2.读取结果

预期结果

先按照x坐标最小，当x坐标相同是，y坐标最小的原则，结果依次应该是1、4、7、2、5、8、3、6、9

实际结果

输入命令

*status,eorder	!显示排序后的结果

与预期结果一致。

总结

看到这里有的人可能会说，这几个单元的排序就需要这么多的命令，那复杂的实体单元是不是更复杂？其实并不是这样的，这是一整套的排序命令，复杂的实体单元与这差不多，只是多了一个Z方向的有关命令。