第三讲有限元一般过程.pptVIP

第三讲 最小势能原理与有限元分析的一般过程 §3-1 最小势能原理 一、结构势能的定义 ——结构势能 ——内力势能，即结构弹性应变能 ——外力势能，数值上为在假定位移下外力的虚功的负值。 结构弹性应变能可以表示为： v —— 结构或弹性体的体积 * 外力虚功可以表示为： —— 结构外力已知表面。 V —— 结构或弹性体的体积 ——结构内任意点的位移向量。 —— 结构体力载荷向量。 —— 结构面力载荷向量。 * 二、最小势能原理的描述 在满足协调条件及给定位移边界条件的一切可能位移中，满足平衡条件的位移使结构势能取极小值。 最小势能原理的变分表达式为： * §3-2 有限元分析的一般过程 一、结构的离散化 将结构或弹性体人为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。 这一步要解决以下几个方面的问题: 1、选择一个适当的参考系，既要考虑到工程设计习惯，又要照顾到建立模型的方便。 2、根据结构的特点，选择不同类型的单元。对复合结构可能同时用到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。 3、根据计算分析的精度、周期及费用等方面的要求，合理确定单元的尺寸和阶次。 4、根据工程需要，确定分析类型和计算工况。要考虑参数区间及确定最危险工况等问题。 5、根据结构的实际支撑情况及受载状态，确定各工况的边界约束和有效计算载荷。 * 在有限元法中通常选择多项式函数作为单元位移插值函数，并利用节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。 位移插值函数需要满足相容（协调）条件，采用多项式形式的位移插值函数，这一条件始终可以满足。 但近年来有人提出了一些新的位移插值函数，如：三角函数、样条函数及双曲函数等，此时需要检查是否满足相容条件。 二、选择位移插值函数 1、位移插值函数的要求 * 形函数的性质： 1）相关节点处的值为 1，不相关节点处的值为 0。 2）形函数之和恒等于 1。 2、位移插值函数的收敛性（完备性）要求： 1） 位移插值函数必须包含常应变状态。 2）位移插值函数必须包含刚体位移。 3、复杂单元形函数的构造 对于高阶复杂单元，利用节点处的位移连续性条件求解形函数，实际上是不可行的。因此在实际应用中更多的情况下是利用形函数的性质来构造形函数。 以阶梯轴的形函数为例 两个形函数分别为 在 节点有： 在 节点有： 在任何点有： 这里我们称 为 的相关节点， 为 的相关节点，其它点均为不相关节点。 * 使用最小势能原理，需要计算结构势能，由弹性应变能和外力虚功两部分构成。结构已经被离散，弹性应变能可以由单元弹性应变能叠加得到，外力虚功中的体力、面力都是分布在单元上的，也可以采用叠加计算。 1、计算单元弹性应变能 —— 单元体积。 由几何关系 代入前式有： 令： 称单元刚度矩阵，简称单刚。 这样单元弹性应变能可以表示为： 三、 单元分析 目的：计算单元弹性应变能和外力虚功。 * 2、计算单元外力功 1) 体力虚功 令： 称单元等效体力载荷向量。 单元体力虚功可以表示为： 2) 表面力虚功 ——单元上外力已知的表面 ,注意！这里只考虑结构的边界表面。 令： 称单元等效面力载荷向量。 单元表面力虚功可以表示为： * 从前面推导可以看出： 单元弹性应变能可计算的部分只有单元刚度矩阵，单元外力虚功可计算的部分只有单元等效体力载荷向量和等效面力载荷向量。在实际分析时并不需要进行上述推导，只需要将假定的位移插值函数代入本节推导得出的单元刚度矩阵、等效体力载荷向量和等效面力载荷向量的计算公式即可。 所以我们说有限元分析的第三步是计算单元刚度矩阵、等效体力载荷向量和等效面力载荷向量。 几点说明： 1）单元刚度矩阵具有正定性、奇异性和对称性三各重要特性。所谓正定性指所有对角线元素都是正数，其物理意义是位移方向与载荷方向一致；奇异性是说单元刚度矩阵不满秩是奇异矩阵，其物理意义是单元含有刚体位移；对称性是说单元刚度矩阵是对称矩阵，程序设计时可以充分利用。 2）按照本节公式计算的单元等效体力载荷向量和等效面力载荷向量称为一致载荷向量。实际分析时有时也采用静力学原理计算单元等效体力载荷向量和等效面力载荷向量，实际应用表明在大多数情况下，这样做可以简化计算，同时又基本上不影响分析结果。 * 四、整体分析 目的：计算整个结构的势能，代入最小势能原理：