08有限单元法课程设计终稿.doc

有限单元法课程设计 学校：辽宁工程技术大学 学院：力学与工程学院 专业班级：工程力学08—2班 设计组成员： 于永军（0816020224） 孙中学（0816020214） 杨立明（0816020219） 钱琳琳（0816020211） 张琳琳（0816020226） 设计题目：连续体平面问题的有限元分析 【题目】：正方形薄板四周受均匀载荷的作用，该结构在边界上受正向分布压力，P=1kN/m，同时在沿对角线y轴上受一对集中压力，载荷为2kN，板厚t=1，泊松比ν=0，见下图： Y 2kN 1kN/m x o 2kN 二，理论依据与分析：此问题，为弹性力学里的平面应力问题，在板的内部，到处都有【1】：σz=0,τyz=0,τxz=0;σx=f1(x,y),σy=f2(x,y),τxy=f3(x,y)，应力具有这种性质的问题，称为平面应力问题。弹性薄板在工程中应用很广泛，对于一些简单的情况，如等厚、单跨、无大孔口，外形规则（如矩形，圆形等）的薄板，已有一些解答和表格可资利用盘【2】。由于连续平板的连续性，仅需要取其在第一象限的四分之一部分研究计算，然后做出一些辅助线将平板分成若干部分，在为每个部分选择分子单元，采用此模型化为4个全等的三角形单元，利用其对称性，四分之一的边界约束，荷载可等效如图所示。 1kN/m 程序原理及实现：用FORTRAN程序的实现。有节点信息文件NODE.IN和单元信息文件ELEMENT.IN,经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件BASIC.IN生成。该程序的特点如下：问题类型：可用于计算弹性力学平面应力问题和平面应变问题。单元类型：采用常应变三角形单元。位移模式：用线性位移模式。载荷类型：节点载荷，非节点载荷应先转换为等效节点载荷。材料性质：弹性体由单一的均匀材料组成。约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有三个自由度的独立约束。方程求解：针对半带宽刚度方程的GUASS消去法。输出文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT,IN。结果文件：输出一般的结果文件DATA.OUT。 程序的原理如框图： 开始  输入数据（子程序READ_IN） BASIC.IN(基本信息文件) NODE.IN(节点信息文件) ELEMEENT.IN(单元信息文件) ↓ 形成单元刚度矩阵（子程序FORM_KE） ↓ 以半带存储方式形成整体刚度矩阵 （BAND K） ↓ 形成节点载荷向量（子程序FORM_P） ↓ 处理边界条件（子程序DO_BC） ↓ 求解方程获得节点位移（子程序SOLVE） ↓ 计算单元及节点应力（子程序） ↓ 输出文件 DATA.OUT ↓ 结束 说明：（1）主要变量： ID： 问题类型码：，ID=1时的平面应力问题，ID=2时的平面应变问题 N_NODE： 节点个数 N_LOAD： 节点荷载个数 N_DOF： 自由度，N_DOF=N_NODE*2（平面问题） N_ELE： 单元个数 N_BAND： 矩阵半带宽 N_BC： 有约束的节点个数 PE：