弹性有限元法及应用.ppt

弹性有限元法及应用 徐戊矫 重庆大学材料科学与工程学院 主要学习内容 有限元法的基础 弹性力学的平面问题 轴对称问题 1. 有限元法的基础 1.1 有限单元法能做什么？ 1.2 弹性力学与其它力学学科之间的关系 1.3 弹性力学的任务、研究方法和基本假设 1.4 弹性力学基本量和基本方程 1.5 微分方程的等效积分形式和加权余量法 1.6 变分原理和里兹方法 1.7 能量原理 1.8 力学分析方法概述 1.1 有限单元法能做什么？ 1.2 弹性力学与各力学学科分支之间的关系 1.3 弹性力学的任务、研究方法和基本假设 固体的“完全弹性”特征是：对应于一定的温度，存在着应力和应变之间一一对应的关系，和时间无关，和加载历史无关。 弹性力学6条基本假设： (1) 连续性假设 (2) 均匀性假设 (3) 各向同性假设 (4) 完全弹性假设 (5) 小变形假设 假定物体在外界因素（例如载荷或温度变化等）作 用下所产生的位移远小于物体原来的尺寸，因而应变分量和转角都远小于1。 在研究物体的平衡时，可不考虑由于变形所引起的物体尺寸和位置的变化；在建立几何方程和物理方程时，可以略去应变转角的二次幂或二次乘积以上的项，使得到的基本方程是线性偏微分方程组。 (6) 无初始应力的假设 假定物体处于自然状态，即在外界因素（载荷和温度变化）作用之前，物体内部是没有应力的。根据这条假设，由弹性力学方法求得的应力仅仅是由于载荷或温度变化的作用产生的。若物体内有初始应力存在，则当物体受外界因素作用时，其内实际存在的应力应等于初始应力加上用弹性力学方法所求得的应力。 1.4 弹性力学的基本量和基本方程 1.5 微分方程的等效积分形式和加权余量法 2 弹性力学的平面问题 2.1 弹性力学平面问题的有限元格式 2.2 有限元解的性质和收敛性 2.3 面积坐标为自然坐标时三角形单元的插值函数及单元矩阵的计算 2.1 弹性力学平面问题的有限元格式 2.2 有限元解的性质和收敛性 2.3 面积坐标为自然坐标时三角形单元的插值函数及单元矩阵的计算 3 弹性力学的轴对称问题 3.1 3节点三角形环状单元的插值函数及应力应变矩阵 3.2 3节点环状单元的单元刚度矩阵 3.3 3节点环状单元的等效节点载荷 ③ 单元应力 在单元应变求得之后，根据物理方程，很方便地得到单元应力 与应变矩阵 相同，应力矩阵 也是常量阵，即3节点三角形单元中各点的应力是相同的。 ④ 单元刚度矩阵 在求得单元应力之后，下面的工作是建立单元应力和单元节点力之间的关系。对单元来讲，单元节点力是节点作用在单元上的外力，则将引起其相应的内力。设此单元各节点有虚位移，则在单元内有相应的虚应变。 单元节点力： 单元内应力： 单元节点虚位移： 单元内虚应变： 因此，外力虚功为： 单元内力虚功为： 根据虚功原理： 单元分析小结 几何方程 物理方程 梳理单元分析的全过程，强调单元分析的目的 单元刚度矩阵是对称矩阵。 单元刚度矩阵中任一元素 的物理意义：当单元的第 j 个节点发生单位位移而其余节点位移为零时，需在单元第 i 个节点位移方向上施加的节点力的大小。 单元刚度矩阵中的每个元素反映了单元刚性的大小，称为刚度系数。 单元刚度矩阵的特性： 1）对称性 各种形式的单元都普遍具有这种对称性质。 2）奇异性 单元刚度矩阵 是奇异的，它不存在逆矩阵。 单元刚度矩阵奇异的物理解释是：单元处于平衡时，节点力 相互不是独立的，它们必须满足平衡方程，因此它们是线性相关的。另外，即使给定平衡的单元节点力，也不能确定单元节点位移，因为单元还可以有任意的刚体位移。 3）主元恒正 主元恒正的物理意义是要在某节点发生节点位移，施加在节点位移方向上的节点力必须与节点位移同向。 ⑤ 单元等效节点载荷列阵 看上式，载荷向量 中的每个元素都是作用在节点上的载荷。若结构受有不直接作用在节点的载荷，则应向节点进行移置，这种移置应根据弹性力学圣维南原理，按照静力等效原则进行。即原载荷和移置到节点的载荷在任何虚位移上所作的虚功都相等。 1）均质等厚单元的自重 2）均布侧压 3）x向均布力 4）x向三角形分布载荷 （3）整体分析： ① 整体刚度矩阵的集成 取出其中的一行来看， 某一个单元中第 j 个节点上的节点力等于该单元三个节点 i, j, m 的节点位移在节点 j 上引起的节点力的迭加。 在整