有限元大作业matlab课程设计例子范文

1.题目概况
矩形板尺寸如下图1，板厚为5mm。

材料弹性模量为
松比μ= 0.27 。

施加约束和载荷并讨论：
图
1 计算简图
1.1基本数据
E = 2⨯105N/mm2，泊
序号载荷约束备注42 向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
1.2分析任务/分析工况
讨论板上开孔、切槽等对于应力分布的影响。

（载荷约束组合不变）。

提示：各种圆孔，椭圆孔随大小、形状、数量，分布位置变化引起的应力分布变化；各种形状，大小的切槽及不同位置引起应力分布的变化等，选择二至三种情况讨论，并思考其与机械零部件的构型的相对应关系。

2.模型建立
2.1单元选择及其分析
由于平板长宽分别为300x100，故可取网格单元大小为1。

如图：
2.2模型建立及网格划分
模型按单元为1 划分后的网格大小如图所示：
2.3载荷处理
向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
3.计算分析
3.1位移分布及其分析
（1）位移分布如图：。

计算力学(有限元方法部分) 程序设计程序说明书程序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欢迎交流与提问！附上邮箱：***************。

祝力学学习顺利！。

【MATLAB 算例】4.8.1(1) 基于4节点四面体单元的空间块体分析(Tetrahedron3D4Node)如图4-22所示的一个块体，在右端面上端点受集中力F 作用。

基于MATLAB 平台，计算各个节点位移、支反力以及单元的应力。

取相关参数为：10110Pa,=0.25E μ=⨯，5=110N F ⨯。

图4-22 一个空间块体的分析解答：对该问题进行有限元分析的过程如下。

（1）结构的离散化与编号将结构离散为5个4节点四面体单元，单元编号及节点编号和坐标如图4-22所示，连接关系见表4-8，节点的坐标见表4-9。

表4-8 单元连接关系单元号 节点号 1 2 3 4 51 42 6 1 43 7 6 7 5 1 6 7 84 1 4 6 7表4-9 节点的坐标节点节点坐标/mxyz 1 2 3 4 5 6 7 80 0 0 0.2 0 0 0 0.8 0 0.2 0.8 0 0 0 0.6 0.2 0 0.6 0 0.8 0.6 0.20.80.6节点位移列阵[]111222888 Tu v w u v w u v w =q (4-190)节点外载列阵34780 0 0 0 TT T T T⎡⎤=⎣⎦F F F F F(4-191)其中34785000 00110N ⎡⎤⎡⎤⎢⎥⎢⎥====⎢⎥⎢⎥⎢⎥⎢⎥-⨯⎣⎦⎣⎦F F F F约束的支反力列阵12560000TTT T T ⎡⎤=⎣⎦R R R R R（4-192其中1256112255661256 x x x x y y y y z z z z R R R R R R R R R R R R ⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥====⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦R R R R总的节点载荷列阵12345678 TT T T T T T T T⎡⎤=+==⎣⎦P F R R R R F F R R F F (4-193)（2）计算各单元的刚度矩阵(以国际标准单位)首先在MA TLAB 环境下，输入弹性模量E 、泊松比NU ，然后针对单元1和单元2，分别5次调用函数Tetrahedron3D4Node_Stiffness ，就可以得到单元的刚度矩阵k1(6×6) ～ k5(6×6)。

有限元大作业——钢架结构分析选题人：日期：2016年6月2日目录：第一章：问题重述 (1)一、题目内容： (1)二、题目要求： (1)第二章：有限元法手工求解 (2)一、平面两单元离散化 (2)二、单元分析 (2)三、单元组装 (5)四、边界条件引入及组装总体方程 (5)五、求解整体刚度方程;计算节点2的位移和转角 (6)六、求节点1、3支撑反力 (6)七、设定数据;求解结果 (7)八、绘制轴力图、弯矩图、剪力图 (8)第三章、matlab编程求解： (9)一、总体流程图绘制： (9)二、输入数据： (9)三、计算单元刚度矩阵： (10)四、建立总体刚度矩阵： (10)五、计算未约束点位移： (10)六、计算支反力： (10)七、输出数据： (10)八、编程： (10)第四章有限元求解 (10)一、预处理 (11)二、模型建立： (12)二、分析计算 (14)三、求解结果 (15)四、绘制图像 (16)第五章结果比较 (19)第六章心得体会 (19)一、王小灿： ...................................................................................................................... 错误!未定义书签。

二、孙明哲： ...................................................................................................................... 错误!未定义书签。

三、张国威 .......................................................................................................................... 错误!未定义书签。