带裂纹板受拉力状态下的边界位移输出

上边界位移获取思路：
首先获取上边界线，再获取上边界线上的节点，再提取节点的位移，输出成txt文件后导出到MATLAB中绘图

!选择37号线上的所有节点
alls,all
lsel,r,,,37
nsll
!获取37号线的节点总数，node_num
*get,node_num,node,0,count

!获取目前节点集中最小节点号，node_min
*get,node_min,node,,num,min

!定义一维数组，维度为节点数目node_num
*dim,A,,node_num        ! 储存37号线上的节点
*dim,Weiyi,,node_num    ! 储存37号线上节点位移
!将最小节点号赋值给一个临时变量
temp_node=node_min

!数组的第一个数为最小节点号
A(1)=temp_node

!循环，利用ndnext函数，逐步获取节点号，存储
*do,i,2,node_num
temp_node=ndnext(temp_node)
A(i)=temp_node
*enddo

!此时数组A存储了该节点集的所有节点号
!后续只需循环遍历该数组即可得到每个节点上的位移
*do,i,1,node_num,1
Weiyi(i) = uy(A(i))
*enddo
FINISH

! 这里需要创建一个宏文件来运行生成位移的txt文件，否则会报错
*create, datawrite, mac
*cfopen, E:\FEM_Files\ANSYS_Files\Crack\Weiyi_Output, txt
*vwrite, Weiyi(1,1)
(f20.15)
*cfclose
*end
Datawrite
FINISH