有限元第二讲杆单元预案.ppt

2 杆单元 一、一维等截面杆单元及其刚度矩阵 二、二维空间中的杆单元 （平面桁架单元） * * L— 杆长 A— 截面积 E— 弹性模量 考虑一个2节点一维等截面杆单元： （一）直接法导出单元特性 杆单元伸长量： 应力—应变关系： ——杆单元位移 ——杆单元应变 ——杆单元应力 应变—位移关系： 应变： 应力： 杆内力： 杆的轴向刚度： 轴向拉压变形模式下，该杆单元的行为与弹簧单元相同，因此杆单元的刚度矩阵为： 比照弹簧元的刚度方程，写出杆单元的刚度方程为： （二）公式法导出杆单元特性 1、单元上假设近似位移场——位移模式 单元上位移假设为简单多项式函数： 用插值法把多项式中的待定系数 转化为节点位移从而得到插值形式的假设位移函数——单元位移模式如下： 上式中： 单元位移模式写成矩阵形式： 注意：采用一次多项式是因为单元只有2个轴向位移分量，对应2个多项式系数。 2、单元应变： ——单元应变矩阵 3、单元应力： 4、应用弹性体虚功原理导出单元刚度方程。 虚功原理 弹性体受力平衡时，若发生虚位移，则外力虚功等于弹性体内的虚应变能。 ——平衡条件 对于杆单元，定义虚位移如下： 节点虚位移： 单元虚位移： 节点力（外力）虚功： 单元虚应变能： 则单元虚应变： 对杆单元应用虚位移原理，得： 考虑到 的任意性，立刻得到： 这就是刚度矩阵的一般形式，可推广到其他类型的单元。 ——杆单元刚度矩阵 对于上面的杆单元： 与前面直接法得到的公式相同！ （三）关于杆单元的讨论 1）在单元坐标系下，每个节点一个未知位移分量，单元共有2个自由度。 2）单元刚度矩阵元素的物理意义： 刚度方程中令： 则： 所以，单元刚度矩阵的第i（i=1,2)列元素表示当维持单元的第i个自由度位移为１，其它自由度位移为０时，施加在单元上的节点力分量。（也可以用此方法直接导出杆单元的刚度矩阵元素，试练习） ３）单元刚度矩阵性质：对称、奇异、主对角元素恒正。 单元刚度方程 （四）举例 例1 求图示２段杆中的应力。 解：分２个杆单元，单元之间在节点２连接。各单元的刚度矩阵分别为： 参考前面弹簧系统的方法，装配２杆系统的有限元方程（平衡方程）如下： 引入边界位移约束和载荷： 则系统平衡方程化为： 上述方程组中删除第１，３个方程，得到： 解得： 位移解： 单元1应力： 单元2应力： 提示： 1）本例中单元应力的计算采用了材料力学中的方法,与采用有限元单元应力公式 的结果相同。 2）对锥形杆，单元截面积可用平均值。 3）求应力之前需要求出节点位移——有限元位移法。 例2： 已知： 求：杆两端的支反力 解 （一）2-D空间中杆单元 1-D空间杆单元 2- D空间杆单元 坐标变换 原来1-D空间中的杆坐标系作为局部坐标系 每节点2个dof 每节点一个dof 总体 局部 节点位移向量的坐标变换： 向量的坐标变换矩阵为： 显然是正交阵，即： 单元节点位移向量的变换式如下： 或 单元节点力的变换为： 刚度矩阵的坐标变换 局部坐标系下杆单元的刚度方程为： 扩充到4自由度形式： 写成矩阵符号形式： 利用前面的向量坐标变换式，得： 考虑到变换矩阵的正交性，得：