有限元平面问题MATLAB程序

计算力学(有限元方法部分) 程序设计

程序说明书

程序1：平面问题的有限元分析

文件名：h01.m

算例文本：h01.txt

输出文本：result1.txt

使用前请先将h01.txt放入默认的文本读取路径（我的要求与m文件在同一文件夹内）！

文本输入顺序：

材料信息（编号、弹性模量、泊松比）

注意：材料编号须按1、2、3、4……的顺序排列

节点信息（编号、x坐标、y坐标）

注意：节点编号须按1、2、3、4……的顺序排列

约束信息（约束节点号、x方向有无约束、y方向有无约束、x方向位移、y 方向位移）有约束处填一正数，无约束处填0。无约束处请勿输入位移。

单元信息（厚度、材料编号、节点编号，若为3节点单元，则第四个编号填0）

注意：单元编号须按1、2、3、4……的顺序排列

集中力（作用节点号、x方向力、y方向力）

线荷载（作用边上的两个节点号、x方向分布力、y方向分布力）

面荷载（作用单元号、x方向分布力、y方向分布力）

程序可调部分：

4-6行中可以指定输出哪些图像（按顺序依次为节点、单元图像，x、y方向位移图像，xx、yy、xy方向应力图像），第7行中可以指定输入的.txt文本名称。

文本输出内容：

结点位移信息（节点号、x方向位移、y方向位移）

单元形心处的应变信息（单元号、x方向正应变、y方向正应变、xy方向工程切应变）

单元形心处的应力信息（单元号、x方向正应力、y方向正应力、xy方向切应力）

本程序附有三角形单元自动加密前处理部分h01auto.m，其算例文本：

h01coarse.txt，输出文本：h01new.txt。它可以适用于本题的要求，在已给定本题所需全部信息的情况下将已有的单元加密为三角形单元。其输出文本可直接作为上面程序的输入文本。

h01.m h01.txt h01auto.m h01coarse.txt

欢迎交流与提问！附上邮箱：***************。祝力学学习顺利！