基于Matlab语言的按平面三角形单元划分的结构有限元程序设计.doc

基于Matlab语言的按平面三角形单元划分的结构有限元程序设计
专业: 建筑与土木工程
班级: 建工研12-2
姓名:
学号:
基于Matlab语言的按平面三角形单元划分结构有限元程序设计
有限单元发及Matlab语言概述
1. 有限单元法
随着现代工业、生产技术的发展,不断要求设计高质量、高水平的大型、复杂和精密的机械及工程结构。为此目的,人们必须预先通过有效的计算手段,确切的预测即将诞生的机械和工程结构,在未来工作时所发生的应力、应变和位移因此,需要寻求一种简单而又精确的数值分析方法。有限单元法正是适应这种要求而产生和发展起来的一种十分有效的数值计算方法。
有限元法把一个复杂的结构分解成相对简单的“单元”,各单元之间通过结点相互连接。单元内的物理量由单元结点上的物理量按一定的假设内插得到,这样就把一个复杂结构从无限多个自由度简化为有限个单元组成的结构。我们只要分析每个单元的力学特性,然后按照有限元法的规则把这些单元“拼装”成整体,就能够得到整体结构的力学特性。
有限单元法基本步骤如下:
(1)结构离散:结构离散就是建立结构的有限元模型,又称为网格划分或单元划分,即将结构离散为由有限个单元组成的有限元模型。在该步骤中,需要根据结构的几何特性、载荷情况等确定单元体内任意一点的位移插值函数。
(2)单元分析:根据弹性力学的几何方程以及物理方程确定单元的刚度矩阵。
(3)整体分析:把各个单元按原来的结构重新连接起来,并在单元刚度矩阵的基础上确定结构的总刚度矩阵,形成如下式所示的整体有限元线性方程:
①
式中,是载荷矩阵, 是整体结构的刚度矩阵, 是节点位移矩阵。
(4)载荷移置:根据静力等效原理,将载荷移置到相应的节点上,形成节点载荷矩阵。
(5)边界条件处理:对式①所示的有限元线性方程进行边界条件处理。
(6)求解线性方程:求解式①所示的有限元线性方程,得到节点的位移。在该步骤中,若有限元模型的节点越多,则线性方程的数量就越多,随之有限元分析的计算量也将越大。
(7)求解单元应力及应变根据求出的节点位移求解单元的应力和应变。
(8)结果处理与显示进入有限元分析的后处理部分,对计算出来的结果进行加工处理,并以各种形式将计算结果显示出。
2. Matlab简介
在用有限元法进行结构分析时,将会遇到大量的数值计算,因而在实用上是一定要借助于计算机和有限元程序,才能完成这些复杂而繁重的数值计算工作。而Matlab是当今国际科学界最具影响力和活力的软件。它起源于矩阵运算,并已经发展成一种高度集成的计算机语言。它提供了强大的科学计算,灵活的程序设计流程,高质量的图形可视化与界面设计,便捷的与其他程序和语言接口的功能。Matlab在各国高校与研究单位起着重大的作用。
“工欲善其事,必先利其器”。如果有一种十分有效的工具能解决在教学与研究中遇到的问题,那么Matlab语言正是这样的一种工具。它可以将使用者从繁琐、无谓的底层编程中解放出来,把有限的宝贵时间更多地花在解决问题中,这样无疑会提高工作效率。目前,Matlab已经成为国际上最流行的科学与工程计算的软件工具,现在的Matlab已经不仅仅是一个“矩阵实验室”了,它已经成为了一种具有广泛应用前景的全新的计算机高级编程语言了,有人称它为“第四代”计算机语言,它在国内外高校和研究部门正扮演着重要的角色。Matlab语言的功能也越来越强大,不断适应新的要求提出新的解决方法。可以预见,在科学运算、自动控制与科学绘图领域Matlab语言将长期保持其独一无二的地位。为此,本例采用Matlab语言编程,以利用其快捷强大的矩阵数值计算功能。
问题描述
一矩形薄板,一边固定,承受150kN集中荷载,结构简图及按平面三角形单元划分的有限元模型图如下所示。
材料参数:弹性模量;泊松比:;薄板厚度。在本例中,所取结构模型及数据主要用于程序设计理论分析,与工程实际无关。
参数输入:
单元个数NELEM=4
节点个数NNODE=6
受约束边界点数NVFIX=2
节点荷载个数NFORCE=1
弹性模量YOUNG=2e8
泊松比POISS=
厚度THICK=
单元节点编码数组LNODS =
节点坐标数组COORD =
节点力数组FORCE =[4 0 -150]
约束信息数组FIXED =
以上数值数据为程序运行前输入的初始数据,存为“”文本格式,同时必须放在Matlab工作目录下,路径不对程序不能自动读取指定初始文件,运行出错。初始数据文本文件输入格式如下图:
Matlab语言程序源代码:
程序中变量说明
NNODE 单元节点数
NPION 总结点数
NELEM 单元数
NVFI