线性有限元方法.pptVIP

第一章 线性有限元方法 线性本构关系（1/3） 线性本构关系（2/3） 线性本构关系（3/3） 几何关系 控制方程（1/6） 控制方程（2/6） 控制方程（3/6） 控制方程（4/6） 控制方程（5/6） 控制方程（6/6） 有限元方程的建立（1/3） 有限元方程的建立（2/3） 有限元方程的建立（3/3） 思考题 谢谢！ 请提问、进行讨论！ * 目的：以固体的弹性静力学为例，介绍线性有限元方程建立过程中的关键环节。 地位：以后章节的预备知识 线性本构关系 几何关系 控制方程 有限元方程的建立 有限单元形状函数的构造 有限单元形状函数插值精度 四边形单元向三角形单元的退化 内容： 线弹性本构关系（Hooke定律） Hooke定律的张量形式： Hooke定律的矩阵形式： 对各向同性材料，Dijkl与[D]的显式很容易导得，例如 [D]的元素个数：原始矩阵为6×6＝36个； 根据对称性变成21个； 对于正交各向异性（复合）材料为9个； 对于横观各向同性（复合）材料为5个； 对于各向同性材料为2个。 拉梅（Lame）常数与工程参数（杨氏模量、泊松比）的关系。 若定义平均正应变和平均正应力（或称八面体正应力）分别为 体积、剪切模量与工程参数（杨氏模量、泊松比）的关系。 再定义应力偏量和应变偏量分别为 本构关系又可写成 在塑性力学中广泛使用。 注：这些本构方程都可化成轴对称、平面应力、平面应变等问题的具体形式。 小变形问题的几何关系为： 应变－位移关系的张量形式： 应变－位移关系的矩阵形式： 微分算子矩阵 （作用于形函数 N） 有限元插值 注：可化成轴对称，平面应力、应变等问题的具体形式。 微分（强）形式的控制方程：（应用力学常用） 1）张量形式： 2）矩阵形式： 分布体积力 与几何关系中相同的微分算子矩阵（6×3） 注：只有在笛卡尔坐标系中才有这样的关系 几何关系 边界条件： 1）自然（面力）边界条件： 2）本质（位移）边界条件： [T]与[B]相似，但用法向余弦替换其中的偏微分。 on on 注：自然边界条件将包含在积分形式的控制方程中。 注：位移边界条件就是通常所说的约束条件，必须予以强制满足。 积分（弱）形式的控制方程：（有限元法常用） 1）虚功（位移）原理： 2）势能原理： 表面分布力 注：上述都是仅以位移为基本变量的原理（弱形式）；变分形式因问题会有所不同，如梁和板问题。 积分（弱）形式的控制方程： 3）基于混合变量的广义变分原理－三变量变分原理 注：原我校教授匡震邦老师在这方面做出了重要贡献，于1983年获得国家自然科学二等奖。 胡－鹫（简称HW）变分原理的广义泛函为 三类变分宗量 可以证明，HW变分原理等价于材料本构关系、几何关系、微分控制方程及面力边界条件。 积分（弱）形式的控制方程： 4）基于混合变量的广义变分原理－两变量变分原理 海－耐（简称HR）变分原理的变分为 两类变分宗量 可以看出，HR变分原理等价于几何关系、微分控制方程及面力边界条件。 注：1）－4）都含有位移做为基本变量，可统称为势能类变分原理。 5）余能原理 特点：以（应）力为基本变量 积分（弱）形式的控制方程： 5）Galerkin加权残值法－数学上更加通用和一般 一般步骤： 对微分控制方程、边界条件等的残差在积分意义下进行加权平均； 识别（确定）权函数； 进行求解。 派生的方法包括：最小二乘法、配点法、子域法、迦辽金法和有限体积法等。 适用的问题包括：热传导、流体力学、空气动力学、电磁场、声场等几乎所有用二阶偏微分方程描述的问题。 本构关系、几何关系及控制方程，无论是采用解析方法，还是采用半解析方法，对求解一个固体力学问题都必须用到。有限元方法与这些方法的最主要区别在于，它一开始就着眼于通过数值途径获得结果，因而，在离散的有限单元上进行插值逼近是有限元方法的主要特征。 单元上的位移场 单元的形状插值函数 单元上的结点位移 单元刚度矩阵 单元等效载荷 单元位移 单元应变 单元应力 例如，利用虚功原理 单元刚度矩阵 单元等效载荷 有限元方法是基于单元的，这是与其它类型数值方法的本质区别，它的优点和缺点都是因之而产生。 有限元方程通常都是在单元层面上建立的，但必须经过组装，在整个结构层面上施加约束、加以求解。 经组装 总体刚度矩阵 总体等效载荷 注：在单元刚度矩阵逐个组装成总体刚度矩阵的过程中，首先要确定总体刚度矩阵的存贮