有限元程序设计课程设计.docVIP

重庆大学本科学生课程设计任务书 课程设计题目 有限元程序设计 学院 航空航天学院 专业 工程力学 年级 已知参数和设计要求： 独立完成有限元程序设计。 独立选择计算算例，并能通过算例判断程序的正确性。 独立完成程序设计报告，报告内容包括理论公式、程序框图、程序本体、计算算例，算例结果分析、结论等。 学生应完成的工作： 复习掌握有限单元法的基本原理。 掌握弹性力学平面问题3节点三角形单元或4节点等参单元有限元方法的计算流程，以及单元刚度矩阵、等效节点载荷、节点应变、节点应力和高斯积分等的计算公式。 用Fortran语言编写弹性力学平面问题3节点三角形单元或4节点等参单元的有限元程序。 在Visual Fortran 程序集成开发环境中完成有限元程序的编辑和调试工作。 利用编写的有限元程序，计算算例，分析计算结果。 撰写课程设计报告。 目前资料收集情况（含指定参考资料）： 王勖成，有限单元法，北京：高等教育出版社，2002。 O.C. Zienkiewicz, R. L. Taylor, Finite Element Method, 5th Edition, McGraw-Hall Book Company Limited, 2000。 张汝清，董明，结构计算程序设计，重庆：重庆大学出版社，1988。 课程设计的工作计划： 第1周星期一上午：教师讲解程序设计方法，程序设计要求和任务安排。 第1周星期一至星期二完成程序框图设计。 第1周星期三至第2周星期四完成程序设计。 第2周星期五完成课程设计报告。 任务下达日期 2014 年 6 月 6 日 完成日期 2014 年 6 月30 日 指导教师 （签名） 学 生 （签名） 前言。 随着科技的发展，理论分析越来越不能满足各行各业的需求，而大量的实验又耗资巨大，所以依靠计算机的分析方法越来越流行。基于有限单元法的abaqus等软件就是通过数值方法来解决工程问题的一种高效的方法。 有限元法的基本思想是“化整为零”，把一个复杂的结构分解成相对简单的“单元”，各单元之间通过结点相互连接。单元内的物理量由单元结点上的物理量按一定的假设内插得到，这样就把一个复杂结构从无限多个自由度简化为有限个单元组成的结构。我们只要分析每个单元的力学特性，然后按照有限元法的规则把这些单元“拼装”成整体，就能够得到整体结构的力学特性。 本文要探讨的事有限单元法的核心，利用fortran编程来编写二维有限单元法的计算过程，并通过abaqus来验证其准确性，然后在用abaqus详细分析一些二维的问题，从而较深对有限单元法的理解。 理论依据。 虚位移原理和虚应原理分别是平衡方程与力的边界条件和几何方程与位移边界条件的等效积分形式。将物理方程引入虚位移原理和虚应力原理可以分别导出最小位能原理和最小余能原理，它们本质上和等效积分的伽辽金“弱”形式相一致，这是建立弹性力学有限原方程一般表达格式的理论基础。 弹性力学平面问题有限元分析步骤为： 单元定义 单元离散化 单元推导 总装 联立方程组求解 平面4、8 节点有限元公式及计算原理 1、平面问题基本公式。 (1)、有限元平衡方程： (2)、单元刚度矩阵： (3)、单元等效结点载荷： 2.平面问题计算步骤。 （1）通过Serendipity四边形单元格式构造插值函数。对于4节点单元，插值函数为： （i=1,2,3,4） 对于8节点单元，插值函数为： （i=1,2,3,4） （2）通过， 平面应力问题： 平面应变问题： 式中E 和v 分别为材料Young 氏模量和Poisson 比。 计算出D和B矩阵。 （3）单元矩阵的变换： 所以面积微元可表示成。 （4）运用高斯积分公式求积分公式近似数值积分解： 其中权函数 （5）单元刚度矩阵与载荷列向量可表示为： ， 将其代入整体刚度矩阵与载荷列向量可得整体刚度矩阵K和载荷列向量P。 （6）总势能为：。 根据最小势能原理可得： 运用上式即可求得求得各节点位移。 轴对称问题。 对于轴对称问题唯一不同于平面问题的是轴对称问题多了一个旋转分量，从而使问题变得更加复杂。下面简单介绍一下轴对称问题的区别于平面问题的特点。 轴对称应变：