现代设计方法第13章 杆系结构有限元法.pptVIP

4） 因 消除整体刚度矩阵中对应的行和列以及载荷向量中对应元素 解得 于是位移向量 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 5） 计算各杆应力，由式（13-24） 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 13.3 平面梁单元 结点位移向量 结点力向量 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 图13-6 平面梁单元 13.3.1 单元的位移模式 轴向位移取为坐标的一次式，挠度取为坐标的三次多项式： （13-25） 写为矩阵形式： 式中： （a） 为待定常数。 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 将对应于梁单元的结点轴向位移记为： 对应于梁单元的结点挠度和转角记为： 将梁单元的结点坐标代入式（13-25）并注意到 ： 可得 ： （b） (c ) 由式（b）可获得待定常数 （d） 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 （e） 将（d）代入（a），可得用结点位移表示的位移模式 （f） 式中 可改写式（f）为 （13-26） (a) 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 式中： 形函数表示为 位于第1，4行、列 位于2，3，5，6行、列 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 （13-27） 13.2.2 单元的应变和应力 梁单元的应变可以分为拉压应变 和弯曲应变 。如不考虑剪切应变对于梁挠度的影响，则 （13-28） （13-29） 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 单元的应力为： （13-30） 13.2.3 单元刚度矩阵 单元内的虚应变为 单元内应力在虚应变上作的虚功为 （g） 单元外力作的虚功有两部分：梁上分布载荷 在虚位移上作的虚功和单元结点力 在结点虚位移上作的虚功，即 (h) 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 由（g）（h）两式得出 即 （i） 令 （i）式改写为 （13-31） 式中： （j） 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 单元刚度矩阵为 将（13-29）代入上式，积分得 （13-33） 式中 是梁截面对主轴的惯性矩，A是梁的截面积 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 13.3.4 等效结点力的计算 由式（j）知，按静力等效原理将分布载荷移置到单元结点上的等效结点力为 （13-34） 对于分布轴向力、分布横向力、分布弯矩等的等效结点力可以分别计算。这里给出几种分布载荷的计算结果，其它载荷情况可以自行推导。也可以由材料力学手册查得相应的梁端反力的计算式并反号作为等效结点载荷。 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 重庆大学机械工程学院 * （c） （c） 第13章 杆系结构 重庆朝天门大桥 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 13.1 一维杆单元 13.1.1 一维杆单元及其位移模式 一维杆单元如图13.1所示，单元长度为L，横截面积为A，材料弹性模量为E。单元的两个结点记为i、j，单元只承受轴向力。取一维坐标O x为由i 指向 j，结点位移向量记为 结点力向量记为 　　　　　　　　　　　　　　　　　　　　　　 线性位移模式为 图13.1 一维杆单元 现代设计方法——第13章 杆系结构 * 重庆大学机械工程学院 （13-3） 13.1.2 一维杆单元的形函数 利用单元两结点处的位移条件： 解出 并回代入（13-3），按结点位移 整理