弹性力学平面问题的有限元法.ppt

弹性力学平面问题的有限元法 一、建立有限元模型 离散化，前处理 二、单元分析 位移假设，几何方程，物理方程， 等效节点力，单元刚度矩阵 三、总刚度矩阵、总载荷列阵、引入支撑 最小势能原理 集成技术（叠加法、对号入座） 四、求解、结果分析，后处理 一、建立有限元模型 离散化：几何、载荷 弹性体 有限元模型 分布载荷 等效节点载荷 注意：在有限元模型中，单元与单元、单元与支 座只在节点相连，载荷作用在节点上。 二、单元分析 以平面三角形常应变单元为例 u v x y o i j m i i j j m m v v u u 注意：一个单元的节点数目（ 3 ）； 一个节点的自由度数目（可能位移， 2 ）。 二、单元分析 问题：当节点位移一定，相应的节点力为多大？ u v x y o i j m i i j j m m v v u u 节点力与 节点位移 之间的比 例常数为 刚度系数 目标：建立单元刚度矩阵。 方法：弹性力学的位移法加虚功原理， 皆用节点位移表示。 u v x y o i j m i i j j m m v v u u v u 1 、位移假设 ： 设各点的位移， 即位移场为 y a x a a v y a x a a u 6 5 4 3 2 1 ? ? ? ? ? ? 6 2 1 , , a a a ? 其中 为待定系数，由节点 坐标和节点位移确定。 将位移写成矩阵形式： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 6 5 4 3 2 1 1 0 0 0 0 0 0 1 a a a a a a y x y x v u } ]{ [ } { a H f s ? 设单元各节点坐标和节点位移为已知，问单元内 任一点的位移？ u v x y o i j m i i j j m m v v u u v u y a x a a v y a x a a u 6 5 4 3 2 1 ? ? ? ? ? ? 在节点上满足： m m m m m m j j j j j j i i i i i i y a x a a v y a x a a u y a x a a v y a x a a u y a x a a v y a x a a u 6 5 4 3 2 1 6 5 4 3 2 1 6 5 4 3 2 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 写成矩阵形式： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 6 5 4 3 2 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 a a a a a a y x y x y x y x y x y x v u v u v u m m m m j j j j i i i i m m j j i i } ]{ [ } { a A b e ? ? u v x y o i j m i i j j m m v v u u v u 解线性代数方程组，得 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 6 5 4 3 2 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 a a a a a a y x y x y x y x y x y x v u v u v u m m m m j j j j i i i i m m j j i i } ]{ [ } { a A b e ? ? } { ] [ } { 1 e b A a ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 6 5 4 3 2 1 1 0 0 0 0 0 0 1 a a a a a a y x y x v u } ]{ [ } { a H f s ? 代入 得 1 6 1 6 6 6 2 1 2 } { ] [ ] [ } { ? ? ? ? ? ? e b s A H f ? 1 6 6 2 1 2 } { ] [ } { ? ? ? ? e f N f ? 1 ] ][ [ )] , ( [ ? ?