有限元模型就是有单元和节点组成的,正确分析完之后总能得到单元解和节点解,至于单元解和节点解的关系,可以参考这篇文章
ANSYS中的节点解与单元解是怎么回事?下次别说你还不懂
下面我们通过一个悬臂梁的例子来看看单元解和节点解的关系,以及如何将单元解和节点解保存为txt文件。
悬臂梁尺寸为10*10*100mm,右端面固定,左端面受到竖直向下的作用力
选择solid185单元,材料弹性模型E=2e11,泊松比为0.3,六面体网格,划分网格并添加边界条件,结果如下:
计算结果如下:
节点位移和节点等效应力
单元等效应力云图
保存节点解和单元解,这里以节点应力和单元应力为例,其他解可以参考下面这篇博客APDL获取节点和单元的结果
获取节点解使用*GET命令流,获取节点应力结果命令流如下:
序号 结果 APDL语句 说明
1 X方向应力 *GET,SX,NODE,i,S,X *GET,变量名,类型,编号,应力,X方向
2 Y方向应力 *GET,SY,NODE,i,S,Y *GET,变量名,类型,编号,应力,Y方向
3 Z方向应力 *GET,SZ,NODE,i,S,Z *GET,变量名,类型,编号,应力,Z方向
4 XY剪应力 *GET,SXY,NODE,i,S,XY *GET,变量名,类型,编号,应力,XY方向
5 YZ剪应力 *GET,SYZ,NODE,i,S,YZ *GET,变量名,类型,编号,应力,YZ方向
6 XZ剪应力 *GET,SXZ,NODE,i,S,XZ *GET,变量名,类型,编号,应力,XZ方向
7 主应力 *GET,S1,NODE,i,S,1 *GET,变量名,类型,编号,应力,主应力
8 第二主应力 *GET,S2,NODE,i,S,2 *GET,变量名,类型,编号,应力,第二主应力
9 第三主应力 *GET,S3,NODE,i,S,3 *GET,变量名,类型,编号,应力,第三主应力
10 von Mises应力 *GET,SVON,NODE,i,S,EQV *GET,变量名,类型,编号,应力,Mises等效应力
具体命令流如下:
/POST1 !后处理
SET,1 !从结果文件中读出所指定的数据集
*GET,MaxNodeNum,NODE,,NUM,MAX !获取节点号
*GET,MinNodeNum,NODE,,NUM,MIN !获取小节点号
*GET,NodeNum,NODE,,COUNT !获取节点数
*DIM,Stress,ARRAY,NodeNum,2 !定义2维数组NodeNum*2
*DO,i,MinNodeNum,MaxNodeNum
Stress(i-MinNodeNum+1,1) = i
*GET,Stress(i-MinNodeNum+1,2),NODE,i,S,EQV
*ENDDO
*CFOPEN,myseqv100,TXT !打开一个文件
*VWRITE,Stress(1,1),Stress(1,2)
(F20.0,,F20.5) !格式输出20表示数据总字符宽度,5表示小数部分宽度
*CFCLOS !关闭一个文件
上面代码因为使用了*VWRITE命令,需要将其保存为txt文件,用导入的方式运行代码。运行完之后在ANSYS的工作目录下生成myseqv100.txt文件,文件结果如下:
列是节点号,第二列是节点von Mises等效应力,同时我们在ANSYS菜单Parameters->Scalar Parameters中可以看到节点数