第十二讲有限元方法初步.pptVIP

南京航空航天大学飞机设计研究所 结构有限元分析方法简介 结构有限元分析实例 如何计算： 弯矩、剪力和轴向力 弯曲应力、剪切应力和拉伸应力 弯曲应变、剪切应变和拉伸应变 结构变形 有限元分析基本思想 难以建立复杂结构外力与位移、内力、应力、应变之间的函数关系。 因此，把大型复杂的结构分解为许多单元，这些单元是简单结构。 例如，简单的梁。 在这些单元中，外力与位移、应力、应变、应变之间的关系很容易确定。 最后，我们把这些单元组合在一起，就完成了大型复杂结构的力学分析难题。 有限元模型 对于每个单元，需要建立如下物理量的关系 外力 变形（角位移/线位移） 应力/应变 建立外力与位移的关系 五个基本步骤： 构造位移形函数 指定节点边界条件 导出应变/位移的关系 导出应力/位移的关系 应用虚功原理建立求解方程组 实例分析 实例分析（续） 构造位移形函数 实例分析（续） 边界条件 实例分析（续） 把公式(2)代入公式 (1) 实例分析（续） 对位移求导，获得应变/位移关系 实例分析（续） 应用弹性原理，获得应力/位移关系 实例分析（续） 虚功原理 实例分析（续） 由虚功原理，建立方程： 实例分析（续） 单元坐标系转换为世界坐标系 单元坐标系与世界坐标系 需要进行坐标变换 单元坐标系转换为世界坐标系 单元坐标系转换为世界坐标系 位移 单元坐标系转换为世界坐标系 世界坐标系下单元力/位移关系 单元坐标系转换为世界坐标系 世界坐标系下的单元刚度矩阵 总刚度矩阵 单元和节点编码 总刚度矩阵 总刚度矩阵 总刚度矩阵 总刚度矩阵 总刚度矩阵 约束处理 如何约束问题 ? 约束处理 分解 计算单元应力、应变 计算单元坐标系下的位移 计算单元应变 计算单元应力、应变 计算单元应力 计算单元轴向力 FEA工作流程总结 * 南京航空航天大学 飞机设计研究所 约束 节点 单元 外力 f1 f2 u1 u2 x 节点 单元 变形后的形状 x u(x) u(0) = u1 u(L) = u2 变形后 获得满足边界条件的位移多项式 应变 应力 弹性模量 内力的功 外力的功 可以表示为梁的截面面积 A 功 = 应力 x 应变 x 体积 刚度矩阵 X Y Z 世界坐标系 u v 单元坐标系 u v 单元坐标系 FX FZ fu fv q fu1 fv1 fv2 fu2 Fx1 Fz1 Fx2 Fz2 力 q u1 v1 v2 u2 X1 Z1 X2 Z2 q 单元坐标系下 世界坐标系 单元坐标系下节点与节点位移的关系 单元坐标系下单元刚度矩阵 世界坐标系下单元刚度矩阵 世界坐标系下节点与节点位移的关系 单元号 1 2 3 节点号 1 2 3 -45 90 L m P [kN] 1,3 1,1 3,1 3,3 2,3 2,2 3,2 3,3 2,1 1,1 2,1 2,2 节点编号为 a 和 b a b i 对第 i 个单元 第 i 个单元在总刚中的位置 总刚 节点位移向量 节点 力向量 求节点位移 求力 例如， 对单元1 负号代表压应力 对单元1 负号代表受压 对单元1 对单元1