有限元分析 第二讲 杆单元讲解.pptVIP

2 杆单元 ? 一、一维等截面杆单元及其刚度矩阵 ? ? ? ? ? ? ? j i u u d 单元节点位移： L — 杆长 A — 截面积 E — 弹性模量 考虑一个 2 节点一维等截面杆单元： ? ? ? ? ? ? ? j i f f f 单元节点力： （一）直接法导出单元特性 ? 杆单元伸长量： i j u u ? ? ? 应力 — 应变关系： dx du ? ? ? ? E ? ) ( ) ( ) ( x x x u u ? ? ? ? ? ? ? —— 杆单元位移 —— 杆单元应变 —— 杆单元应力 应变 — 位移关系： ? 应变： L ? ? ? ? 应力： L E E ? ? ? ? ? ? 杆内力： ? ? ? ? ? ? ? ? k L EA L EA A F ? ? 杆的轴向刚度： L EA k ? ? 轴向拉压变形模式下，该杆单元的行为与弹簧单元相同， 因此杆单元的刚度矩阵为： ? 比照弹簧元的刚度方程，写出杆单元的刚度方程为： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? j i j i j i u u L EA u u k k k k f f 1 1 1 1 L EA k ? （二）公式法导出杆单元特性 1 、单元上假设近似位移场 —— 位移模式 单元上位移假设为简单多项式函数： 0 1 ( ) u x a a x ? ? 用插值法把多项式中的待定系数 转化为节点位移 从而得到插值形式的假设位移函数 —— 单元位移模式如 下： 0 1 , a a ( ) ( ) ( ) i i j j u x N x u N x u ? ? ( ) 1 ( ) i j x N x L x N x L ? ? ? 上式中： 形函数。 中称为形状函数，简称 是插值基函数，有限元 j i N N , 单元位移模式写成矩阵形式： ? ? Nd ? ? ? ? ? ? ? ? j i j i u u N N u ( ) ( ) ( ) i i j j u x N x u N x u ? ? 注意：采用一次多项式是因为单元只有 2 个轴向位 移分量，对应 2 个多项式系数。 ? ? 称为单元形函数矩阵。 , j i N N ? N 式中 2 、 单元应变： Bd d N ? ? ? ? ? ? ? ? ? dx d dx du ? —— 单元应变矩阵 B 3 、 单元应力： Bd E E ? ? ? ? 4 、应用弹性体虚功原理导出单元刚度方程。 ? 虚功原理 弹性体受力平衡时，若发生 虚位移 ，则外力虚功等于 弹性体内的虚应变能。 —— 平衡条件 ? ? Nd ? ? ? ? ? ? ? ? j i j i u u N N u ? ? ? ? L L x N x N dx d j i / 1 / 1 ) ( ) ( ? ? ? B 对于杆单元，定义虚位移如下： ? ? ? ? ? ? ? j i u u ? ? ? d 节点虚位移： 单元虚位移： d N ? ? ? u 节点力（外力）虚功： f d T ? 单元虚应变能： d B B d Bd B d T ? ? ? ? ? ? ? ? ? ? ? ? ? dV E dV E dV V V V T T T T ? ? ? ?? 则单元虚应变： d B ? ? ?? ? ? ) ( u dx d 对杆单元应用虚位移原理，得： d ? d B B d f d T ? ? ? ? ? ? ? ? ? ? dV E V T T ? ? 考虑到 的任意性，立刻得到： kd d B B f T ? ? ? ? ? ? ? ? ? ? ? dV E V ? ? V dV E B B k T 这就是刚度矩阵的一般形式，可推广到其他类型的单元。 —— 杆单元刚度矩阵 对于上面的杆单元： 与前面直接法得到的公式相同！ （三）关于杆单元的讨论 1 ）在单元坐标系下，每个节点一个未知位移分量，单元 共有 2 个自由度。 2 ）单元刚度矩阵元素的物理意义： 刚度方程中令： 则： 所以，单元刚度矩阵的第 i （ i=1,2) 列元素表示当维持单元 的第 i 个自由度位移为１，其它自由度位移为０时，施加 在单元上的节点力分量。（也可以用此方法直接导出杆单 元的刚度矩阵元素，试练习） ３）单元刚度