数值模拟第二讲杆单元要点.pptVIP

按公式计算杆应力： 得： 2.2.3 3-D空间中杆单元 1-D空间杆单元 3- D空间杆单元 换变标坐 基本思想 局 部 总 体 每节点1个自由度 每节点3个自由度 节点上向量的坐标变换和单元刚度矩阵的坐标变换原理与2-D情况相同，变换矩阵由杆的3个方向数组成。 3-D总体坐标系下杆单元每节点有3个平动自由度（位移分量），单元有6个自由度。单元刚度矩阵为6×6方阵。 杆单元目 标：通过杆单元特性方程的建立，初步掌握有限元法单元分析的过程和原理，了解杆系结构分析的原理。 2.1、一维等截面杆单元及其刚度矩阵 L— 杆长 A— 截面积 E— 弹性模量 考虑一个2节点一维等截面杆单元： 单元上的力学量和基本关系如下： ——杆单元位移 ——杆单元应变 ——杆单元应力 2.1.1 直接法导出单元特性(方法一) 杆单元伸长量： 应变—位移关系： 应力—应变关系： 下面通过二种方法研究杆单元的单元特性。 应 变： 应 力： 杆内力： 则杆的轴向刚度： 轴向拉压变形模式下，该杆单元的行为与弹簧单元相同，因此杆单元的刚度矩阵为： 比照弹簧元的刚度方程，写出杆单元的刚度方程为： 2.1.2 公式法导出单元特性(方法二) 步骤1) ： 假设单元上近似位移函数——位移模式 假设单元上的位移为简单多项式函数。有限元中用插值法通过节点位移（待定参数）定义单元位移函数。 对杆单元，引入局部坐标： 则单元近似位移函数（线性位移模式）： 定义2个节点的插值函数（形函数）： 矩阵形式： 由于该杆单元只有2个未知位移分量，因此单元上假设的简单位移函数采用一次多项式。故对单元的节点位移进行线性插值。 该式是有限元法中最重要的关系式之一，通过该式把单元上的近似位移分布函数用节点位移来表示，为进行单元层次上的分析打下了基础。 单元应变： ——单元应变矩阵 单元应力： 用虚位移原理导出单元刚度方程 虚位移原理（虚功原理） 弹性体受力平衡时，若发生虚位移，则外力虚功等于弹性体内的虚应变能。 步骤2) ： 单元应变与单元应力 步骤3) ： 弹性体的虚位移： 假想在弹性体上发生的，满足位移许可条件（内部连续，边界协调）的微小、任意位移场。 可以理解为某个位移场的微小扰动（变分）。 虚位移的特征： 1）假想的，与真实位移无关； 2）几何上是许可的：连续、协调； 3）微小、任意大小。 对于杆单元，定义虚位移如下： 节点虚位移： 单元虚位移分布： 节点力（外力）虚功： 单元虚应变能： 则单元虚应变： 对杆单元应用虚位移原理，可得到： 考虑到 的任意性，立刻得到： 这就是刚度矩阵的一般形式，可推广到其他类型的单元。 ——杆单元刚度矩阵 对于上面的杆单元： 与前面直接法得到的公式相同！ 2.1.3 关于杆单元的讨论 1）在单元局部坐标系下，每个节点一个未知位移分量和 一个自由度，单元共有2个自由度。 2）单元刚度矩阵元素的物理意义 刚度方程中令： 则： 单元刚度矩阵的第i（i=1,2)列元素表示当维持单元的第i个自由度位移为１，其它自由度位移为０时，施加在单元上的节点力分量。 （也可以用此方法直接导出杆单元的刚度矩阵元素） ３）单元刚度矩阵对称、奇异、主对角元素恒正。 单元刚度方程 2.1.4 举例 例1 求图示２段杆中的应力。 解：结构分为２个杆单元，单元之间在节点２铰接。 ２个杆单元的刚度矩阵分别为： 参考前面弹簧系统的方法，装配２杆系统的有限元方程（平衡方程）如下： 引入边界位移约束和载荷： 方程化为： 上述方程组中删除第１，３个方程，得到： 解得： 即位移解为： 单元1应力： 单元2应力： 提示： 1）本例中单元应力的计算采用了材料力学中的方法，与采用有限元单元应力公式 的结果相同。 2）对锥形杆，单元截面积可用平均值。 3）求应力之前需要求出节点位移 ——有限元位移法。 例2：求杆两端的支反力 已知： 解： 先检查杆右端与墙壁是否接触。计算右端的自由 伸长： 所以，右端间隙将闭合，即节点3与刚性墙壁接触。 参照前面的讨论，可直接写出2单元系统平衡方程： 载荷与边界条件： 系统平衡方程为： 分离出第2个方程： 即： 得到： 节点位移列式：