材料成形有限元法_有限元程序设计.doc

华 中 科 技 大 学 教 案 章 节 第一章 第一节 周次 第三周 课次 第一次 教 案 内 容 有限元程序设计的基本思想： 一般采用FORTRAN77语言，FORTRAN90语言也开始逐渐被采用。 程序流程标准化（通用性、可读性）。 程序模块化（可移植性、可读性、可调性）。 程序变量名统一化（可读性、可调性）。 有限元程序的一般过程： 最小势能原理（位移法的数学基础）： （1） 其中：(为构形的应变；(为构形的应力；u为构形的虚位移；P为构形上所受的面力；G为构形的体力；(P为构形的势能。 对于线弹性问题而言，构形就是线弹性体。某构形处于平衡状态时，其势能应为最小，即 （2） （3） 单元离散过程： 几何关系（单元内任意点的位移u与应变(之间的关系） （4） 其中L为单元几何关系算子。 本构关系（单元内任意点的应变(与应力(之间的关系） （5） 其中为弹性矩阵。 插值关系（单元内任意点的位移u与单元结点位移ue之间的关系） （6） 其中N为单元插值形函数矩阵 对于某一有限元单元来说，由式（3）~式（6）可得： 华 中 科 技 大 学 教 案 章 节 第一章 第一节 周次 第三周 课次 第一次 教 案 内 容 （7） 其中B=LN。 由于的任意性，得 （8） 最终得到单元刚度方程 （9） 其中：单元刚度矩阵 单元外载荷向量 单元集成过程 （10） （11） 刚度方程约束过程 给出所要求解问题的边界条件。 刚度方程求解过程 由式（11）解出所有结点的位移U。 有限元回代过程 求单元的应变 （12） 求单元的应力 （13） 华 中 科 技 大 学 教 案 章 节 第一章 第二节 周次 第四周 课次 第二次 教 案 内 容 单元几何关系 三维情况，单元的应变与位移的关系为： （1） 即 （2） 二维情况单元的应变与位移的关系可由三维情况退化而来， 平面应力和平面应变情况： （3） 华 中 科 技 大 学 教 案 章 节 第一章 第二节 周次 第四周 课次 第二次 教 案 内 容 轴对称情况 （4） 单元本构关系 三维情况单元的应力应变关系为 （5） 其中 （6） 二维情况可由三维情况退化而来， 平面应力（） （7） 平面应变（） （8） 华 中 科 技 大 学 教 案 章 节 第一章 第二节 周次 第四周 课次 第二次 教 案 内 容 轴对称情况（） （9） 华 中 科 技 大 学 教 案 章 节 第一章 第三节 周次 第五周 课次 第三次 教 案 内 容 有限元法中常采用的单元类型及其形状函数 线性单元（一次单元）： 一维线单元 （1） 其中，l为单元长度；x1和x2为单元结点坐标；x为单元内任意点的坐标。 二维面单元 三角形单元 （2） 其中，A为单元面积 （3） （4） yi、zi（i =1,2,3）为单元结点坐标。 四边形单元 （5） 四边形四个结点的自然坐标为 （r1，s1）=（-1，-1）；（r2，s2）=（1，-1） （r3，s3）=（1，1）；（r4，s4）=（-1，1） 华 中 科 技 大 学 教 案 章 节 第一章 第三节 周次 第五周 课次 第三次 教