1、建立二维可变形实体，赋予简单的属性
2、之后定义相关的材料信息
E=100MPa μ=0
然后是赋予截面属性
3、建立分析步
这里我们为了看到更多的步骤，直接采用固定增量步
（约束完全，一定收敛）
4、施加载荷与约束
然后对顶边施加用户自定义位移
5、创建需要的子程序
打开VS，将帮助文档中关于disp的部分复制粘贴进去
6、然后加上自己定义的东西
这里说明一下，U(1),U(2),U(3)并不是指x,y,z的变形量
而是传递的自由度，以及对应的一阶导，二阶导
所以这里的U(1)指的是y方向的变形量
7、划分网格
采用默认的网格划分方式，种子大小为100
8、提交计算
