有限元方法(有FORTRAN程序)飞箭系列.pdf

FEPG 参考手册 基于 FEPG 的 有 限 元 方 法 北京飞箭软件有限公司 2003 年 3 月 前 言 本教材的目的是为了帮助 FEPG 系统的用户更好地理解和掌握该 系统，从有限元方法原理，从有限元软件的程序结构和数据结构，从 有限元的输入数据形式，以及软件设计方法等方面对系统进行具体阐 述。 众所周知，有限元方法是结构力学家发明的，其基本原理是基于 变分法与分片插值多项式两项技术，但是有限元技术的推广与研究的 深入使得今天的有限元方法已突破了这两项技术的范围。例如，为了 把有限元方法应用于非结构力学领域，如流体力学领域，人们不得不 采用虚位移原理即“弱形式”取代变分原理，为了求解无界区域问题 或者奇点问题，采用非多项式的形函数（基函数）将具有更好的精度， 因此形函数不得不突破分片插值多项式的限制。 FEPG 系统是基于虚位移原理（即弱形式）而不是变分原理，要 求用户书写弱形式的微分方程表达式，因此本教材的第一章通过对不 同领域的微分方程表达式推导其弱形式，让用户对微分方程弱形式的 推导有一般性的了解，但是一个好的弱形式可能是专家经过多年的研 究后取得的成果，并不是轻而易举就可得到的。 教材的第二章给出各种常用的分片插值多项式形函数表达式及 其推导的一般性方法，但本系统允许用户使用非多项式形函数。第三 章给出本系统有限元程序所需的输入数据的一般形式。第四章论述了 FEPG 系统的有限元计算程序结构，并给出全部 FORTRAN 源程序及其 详细的说明。 目 录 第一章 偏微分方程的“弱”形式 ——虚位移原理 §1.1 偏微分方程的弱解形式……………………………………………1 1.1.1 问题的提出……………………………………………………………1 1.1.2 偏微分方程弱解的积分形式虚位移原理…………………………2 §1.2 稳态问题弱解的积分形式…………………………………………2 1.2.1 分部积分公式…………………………………………………………3 1.2.2二维稳态热传导问题的“弱”形式……………………………………4 1.2.3三维线弹性小变形静态问题的“弱”形式……………………………6 1.2.4 三维稳态渗流问题的“弱”形式……………………………………11 1.2.5 二维粘性不可压缩流体稳态 Navier_Stokes 方程的“弱”形式…14 1.2.6 三维静电场问题的“弱”形式………………………………………16 1.2.7 三维柱坐标静电场问题的“弱”形式………………………………18 §1.3 瞬态问题的“弱”形式……………………………………………19 1.3.1 三维线弹性小变形动态问题的“弱”形式…………………………20 1.3.2 三维瞬态热传导问题的“弱”形式…………………………………23 1.3.3 二维粘性不可压缩流体瞬态 Navier_Stokes 方程的“弱”形式…25 §1.4 求解解梯度的最小二乘法………………………………………27 1.4.1 已知位移求应力………………………………………………………28 1.4.2 已知温度求热流密度…………………………………………………29 1.4.3 已知电势求电场强度…………………………………………………30 第二章 分片多项式的形函数 §2.1 插值函数与单元类型………………………………………………31 2.1.1 一维 Lagrange 单元…………………………………………………31 2.1.2 二维单元……………………………………………………………33 2.1.3 三维单元…………………………………