非均匀残余应力有限元模拟步骤

1.根据残余应力的分布，拟合出其分布的多项式F=f(x)
2.建模并划分网格
3.点击写入输入文件，将模型信息写入inp文件：
4.讲inp文件中表示节点坐标信息的部分复制出来，拷贝到0.txt文件中，其第一列表示节点编号，右边三列表示节点对应的坐标
5.讲inp文件中描述单元信息的部分拷贝到1.txt文件中，第一列表示对应的单元编号，后面几列表示与单元想联系的节点的编号
6.用MATLAB对0.txt和1.txt文件分别进行处理，根据节点的坐标，以及拟合出的公式，计算出各个节点处的应力的大小。

7.计算各个单元所受的应力，方法是将与单元联系的各个节点的应力求和再求平均，其结果可以看成该单元处的应力，并将单元的编号以及对应的应力大小保存到文件3.txt中，如下图所示
7,.新建一个excel文件，并将3.txt文件进行导入，结果如下：
8.所建的模型名称为Part-1，其实例名为Part-1-1，因此，得把第一列的单元
编号前加上模型实例名称Part-1-1：，然后利用excel的拖拽功能，将第一列的编号按顺序进行编排，其结果如下：
9.将文件以.CSV格式进行保存，文件名取00.CSV
10,得到文件00.CSV,并将其改名为00.dat：
11.将文件00.dat放入abaqus的工作目录中，并编辑模型关键字，在其中加入语句：initial conditions,type=stess,input=00.dat
12.提交job进行运算，得到如下的结果：。