《安徽建筑工程学院计算结构力学7》精选·课件.pptVIP

程序设计 　　在上一章节解线性代数方程组程序设计的基础上，已知{Δ}，求单元杆端力{F’}和支座反力{R}，只需在上章主程序中增加调用求单元位移{δ}的子程序SUB. YWY(讲义P115或参考教材P203) 　　以及求{F’}和{R}的子程序SUB. QDL (讲义P116或参考教材P204)即可。 这两个子程序中只有一个新数组 FLZ(NJZ,3):支座反力数组 * 首页 上页 返回 下页 * 第七章 结构内力和反力的计算 对结构刚度方程 [K]{Δ}={P} 求解得到{Δ}后，根据单元定位向量即可得到各单元位移{δ}，而进行结构设计还需了解结构的内力及反力情况。 结构的内力可通过单元的杆端力求出。 杆端力:结点对杆端的作用力。当某结点是支座结点时,杆端力就是支座对杆端的作用力,也就是支座反力。 单元的杆端力和支座反力通常由两部分组成： 1.单元两端的位移所引起的杆端力{Fδ} 2.单元荷载所引起的杆端力{FQ} 结构设计要求了解在局部坐标系下的杆端力(绘内力图)和整体坐标系下的支座反力(基础设计) 本章介绍这两种单元杆端力的计算 1.杆端力{Fδ} 求解结构方程组[K]{Δ}={P}得到结构未知量{Δ}。 可直接根据MWe,从{Δ}中取出{δ}e (位移协调原理,有限元离散的准则),然后则可由单元刚度方程,求出杆端力： 7-1 杆端位移所引起的杆端力和支座反力 程序设计实现这一步,可从XSA(NE,7)中取出XS(7),然后由(4-4-17)式形成[k],即可做上式乘法。 2.支座反力 支座结点对杆端的作用力称为支座反力,可通过检查单元定位向量找出支座结点及相应的分量,若MWe相应的分量为零,则该分量对应的杆端力分量就是支座反力相应的某一分量： ○ e 注意 ①这样求出的{δ}与{F}均在整体坐标系 ②{F}与[K]、{δ}成正比,刚度大的就多承担一些力。 7-2 单元荷载所引起的 杆端力和支座反力 以下例来说明(既考虑单元荷载又考虑杆端位移)： 例1:求图示简支梁的杆端力(支座反力)。 (a) 求{F} (b) (c) 显然：(a)=(b)+(c) (b)图中: (c)图中: 关键是(c)中{Fδ}的求法： 式中{F0}为固端反力,各类杆件在各种荷载作用下的固端反力列在附表(5-7)。 结构刚阵经单元定位向量[0 l 0 2]T处理后得: 写出结构方程组[K]{Δ}={P}： 由单刚方程可求出(c)图中{Fδ}: 由(a)=(b)+(c)得到杆端力: 结论：在单元荷载和杆端位移共同作用下的单元杆端力计算公式为: 例2:求杆端力{F}。 (a) (c) (b) 求{F} 首先求出(c)图中的{Fδ}： ∵ MW(e)=[0 0 0 1] ∴ 结构刚度方程为 再由单刚方程[K]{δ}={Fδ}可求出{Fδ} 由(1)式: 内力图:M图 Q图 7-3 在局部坐标系下的杆端力计算 绘制结构的内力图时,需要求出在局部坐标系下的杆端力。前面介绍了由{Δ}中取出{δ}(通过MW),则表明所求的杆端力是在整体坐标下。 在局部系下的杆端力计算有两种方法： 4. 支座反力要在整体坐标系下计算,对于支座结点,还应有: {F}=[T]T{F}’ (2) 再由前面(7-2-1)式介绍，可求支座反力Rj 如有单元荷载,还应叠加： 上式中固端反力总是在局部下 如有单元荷载,还应有： 整体系下的杆端力,对支座结点,可由7-2-1 式直接求出Rj 若存有整体系下的单刚,则采用方法二较为节省(可省去方法一中的第4步变换) 如果未在整体系下的存贮单刚,则采用方法一 本程序用方法二,因程序存有XSA(NE,7)是在整体坐标系下的。程序设计见PSTDY中的子过程 SUB. QDL 讨论