第三章节 杆单元.ppt

2杆单元 2.1引例——弹簧单元 弹簧是宏观上最简单的弹性元件。 1、一个弹簧单元的分析 2.2 杆单元 一、一维等截面杆单元及其刚度矩阵 二、二维空间中的杆单元 （平面桁架单元） 上述方程组中删除第１，３个方程，得到： 解得： 位移解： 单元1应力： 单元2应力： 提示： 1）本例中单元应力的计算采用了材料力学中的方法,与采用有限元单元应力公式 的结果相同。 2）对锥形杆，单元截面积用平均值。 3）求应力之前需要求出节点位移——有限元位移法。 例2： 已知： 求：杆两端的支反力 解 （一）2-D空间中杆单元 1-D空间杆单元 2- D空间杆单元 坐标变换 * 2个节点： 节点位移： 节点力： 弹簧刚度： 已知弹簧力——位移关系： 弹簧力，拉伸为正 — 弹簧伸长 考虑弹簧的特性和平衡关系有： 写成矩阵形式： 矩阵符号形式： ——弹簧单元刚度方程 上式中： 单元节点力向量 单元节点位移向量 弹簧单元的刚度矩阵 刚度方程讨论： 1） 有什么特点？ 2） 中元素代表什么含义？ 3）上面方程可以求解吗？为什么？ 2、弹簧系统 各单元的特性分别为： 单元1： 单元2： 按两种方法装配系统特性： 方法1： 分别考虑节点1，2，3的力平衡条件（节点力与节点外载荷的平衡）： 把单元特性代入，得到： 上面方程写成矩阵形式： 或 （弹簧系统的平衡方程） —— 弹簧系统的结构总刚度矩阵 —— 系统节点载荷列阵 讨论：（1） 有那些特点和性质？ （2）上述方程能求解吗？ —— 系统节点位移列阵 方法2： 将单元刚度方程扩大到整体规模： 将上面的矩阵方程叠加，得到： 代入前面节点平衡条件，得系统节点平衡方程： 3）给定载荷和约束条件下的求解 设边界条件为： 则节点平衡方程为： 该方程展开后分为2个部分： 未知量为2个节点位移： 一个支反力： 解上面方程得： 3、举例：弹簧系统 已知条件： 求：（a） 系统总刚度矩阵 （b） 节点2，3的位移 （c） 节点1、4的反力 （d） 弹簧2中的力 解： （a） 各单元的刚度矩阵为： 应用前面的叠加方法，直接得到弹簧系统的总刚度矩阵： 或 总刚度矩阵特征：对称，奇异、带状、稀疏 由前面的做法，可得到弹簧系统的节点平衡方程： （b）：先施加位移边界条件 将 带入平衡方程后，第2，3方程为： 求解得： （c)：由第1，4个方程求得支反力 (d)：弹簧2内力 （拉力） 4、练习题 对图示弹簧系统，求其总刚度矩阵 解 要点回顾 1、弹簧单元刚度方程的建立 弹簧变形平衡 2、弹簧系统的集成 1）列节点平衡方程法 单元特性 系统节点平衡条件 系统平衡方程 2）单元方程扩大相加法 相加 系统节点平衡条件 单元特性 系统平衡方程 L— 杆长 A— 截面积 E— 弹性模量 考虑一个2节点一维等截面杆单元： （一）直接法导出单元特性 杆单元伸长量： 应变—位移关系： 应力—应变关系： ——杆单元位移 ——杆单元应变 ——杆单元应力 应变： 应力： 杆内力： 则杆的轴向刚度： 轴向拉压变形模式下，该杆单元的行为与弹簧单元相同， 因此杆单元的刚度矩阵为： 比照弹簧元的刚度方程，写出杆单元的刚度方程为： （二）公式法导出杆单元特性 假设单元上近似位移函数——位移模式 一个单元上的位移假设为简单多项式函数。有限元中用插值法通过节点位移（待定参数）定义单元位移函数。 对杆单元，引入局部坐标： 则单元近似位移函数（线性位移模式）： 定义2个节点的插值函数（形函数）： 矩阵形式： 单元应变： ——单元应变矩阵 单元应力： 下面应用弹性体虚位移原理导出单元刚度方程。 虚位移原理（虚功原理） 弹性体受力平衡时，若发生虚位移，则外力虚功等于弹性体内的虚应变能。 对于杆单元，定义虚位移如下： 节点虚位移： 单元虚位移分布： 节点力（外力）虚功： 单元虚应变能： 则单元虚应变： 对杆单元应用虚位移原理，得： 考虑到 的任意性，立刻得到： 这就是刚度矩阵的一般形式，可推广到其他类型的单元。 ——杆单元刚度矩阵 对于上面的杆单元： 与前面直接法得到的公式相同！ （三）关于杆单元的讨论 1）在单元坐标系下，每个节点一个未知位移分量—— 一个自由度，单元共有2个自由度。 2）单元刚度矩阵元素的物理意义： 刚度方程中令： 则： 所以，