有限元-计算结构力学-大作业.docx

有限元-计算结构力学-大作业 有限元-计算结构力学-大作业 PAGE 有限元-计算结构力学-大作业 SHANGHAI JIAO TONG UNIVERSITY 平面应力问题解的Matlab实现 姓 名: heiya168 学 号: 帆哥 班 级： 指导老师：  TOC \o 1-3 \h \z \u 目录 TOC \o 1-3 \h \z \u 1 绪论 1 2 平面问题的四节点四边形单元 2 单元的构造 2 等参变换 5 边界条件的处理——置“1”法 8 3 有限元分析流程 10 程序原理和流程 10 使用的函数 11 文件管理 11 数据文件格式 11 4 算例——开方孔的矩形板拉伸分析 13 问题的具体参数与载荷 13 Matlab程序计算 13 ANSYS建模计算 15 误差分析 17 5 总结 18 参考文献 19 附录 20  PAGE 25 绪论 有限元方法(finite element method)，是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将它应用于工程技术中，可成为工程设计和分析的可靠工具。 弹性体在载荷作用下，其基本方程可写成以下的三类方程和两种边界条件。 平衡方程——应力与外载荷的关系；几何方程——应变位移关系；物理方程——应力应变关系；力的边界条件；几何边界条件。应用最小位能原理，并利用上述关系，最终建立由刚度方程，节点位移和等效节点载荷所构成的求解方程。带入边界条件求解方程，就可以得出弹性力学问题的一般性解答。 本次大作业基于有限元方法的基本原理，使用Matlab这一平台，针对平面应力问题，采用四节点四边形单元编写了求解单元节点位移的程序。主要内容包括：1）介绍有限元的基本原理；2）编程基本思路及流程介绍；3）程序原理及说明； 4）具体算例 这四个部分。 平面问题的四节点四边形单元 单元的构造 单元的几何和节点描述 平面4节点矩形单元如图4-6所示，单元的节点位移共有8个自由度。节点的编号为1、2、3、4，各自的位置坐标为(xi, yi), i=1,2,3,4，各个节点的位移(分别沿x方向和y方向)为(ui,vi), i=1,2,3,4。 图 平面4节点矩形单元 若采用无量纲坐标 （） 则单元4个节点的几何位置为 （） 将所有节点上的位移组成一个列阵，记作；同样，将所有节点上的各个力也组成一个列阵，记作，那么 （） 若该单元承受分布外载，可以将其等效到节点上，也可以表示为如式所示的节点力。利用函数插值、几何方程、物理方程以及势能计算公式，可以将单元的所有力学参量用节点位移列阵及相关的插值函数来表示；下面进行具体的推导。 单元位移场的表达 从图2-1可以看出，节点条件共有8个，即x方向4个，y方向4个，因此，x和y方向的位移场可以各有4个待定系数，即取以下多项式作为单元的位移场模式 （） 它们是具有完全一次项的非完全二次项，以上两式中右端的第四项是考虑到x方向和y方向的对称性而取的，除此外xy项还有个重要特点，就是“双线性”，当x或y不变时，沿y或x方向位移函数呈线性变化，这与前面的线性项最为相容，而或项是二次曲线变化的。因此，未选或项。 由节点条件，在x=xi，y=yi处，有 （） 将式代入式中，可以求解出待定系数a0，…，a3和b0，…，b3，然后代回式(4-52)中，经整理后有 （） 其中 （） 如以无量纲坐标系来表达，式可以写成 （） 将式写成矩阵形式，有 （） 其中，为该单元的形状函数矩阵。 单元应变场的表达 由弹性力学平面问题的几何方程(矩阵形式)，有单元应变的表达 （） 其中几何矩阵B(x,y)为 = （） 式(4-59)中的子矩阵为 （） 单元应力场的表达 由弹性力学中平面问题的物理方程，可得到单元的应力表达式