张系斌结构力学简明教程讲义教案课件3杆系静力分析24.pptxVIP

3.4 杆系结构单元分析子程序3.4.3 平面杆系结构设计3.4.4 空间杆系结构设计 Select Case (Plane) Case (1) 平面桁架元素赋值 Case (2) 平面梁柱元素赋值 Case (3) · · · · · · Case Default 出错信息 End Select Select Case (Space) Case (1) 空间桁架元素赋值 Case (2) 空间梁柱元素赋值 Case (3) 交叉梁元素赋值 Case Default 出错信息 End Select3.4 杆系结构单元分析子程序3.4.5 有关单元等效结点荷载设计和进一步的考虑1) 单元等效结点荷载设计同仿单元刚度。2) 从各类单元刚度元素的计算，可看到要用到长度、 单元弹性特性、单元截面特性等数据。因此，要 确定存放它们的数据结构。要将它们作为出口。3) 为计算单元等效结点荷载元素，首先要建立各种 荷载情况等效荷载表达式，它们可由积分或载常 数表得到。然后要解决荷载信息的存放结构，也 要将它们作为出口量。4) 单元刚度矩阵、等效结点荷载矩阵都应先清零。4.1 杆系结构整体分析 首先就全刚结点平面刚架进行讨论，然后推广。4.1.1 总的思路 在单元特性搞清后，将单元拼装回去。在结点处位移自动协调基础上，如果全部结点平衡，则求得的结点位移将是实际结构的解。因此，整体分析就是设法建立结点平衡方程。4.1.2 坐标转换 组成结构的杆件可以各个方向，单元分析对局部坐标，因此，必须将物理量转为统一坐标——整体坐标。1) 力的转换关系4.1 杆系结构整体分析2) 位移转换关系3) 转换矩阵 转换矩阵是正交矩阵。 如果记 称为整体单元刚度矩阵4.1 杆系结构整体分析4) 杆端力转换5) 杆端位移转换本节以后的讨论认为都是对整体坐标的6) 刚度方程的转换 则 这就是整体坐标下的单元刚度方程。②①③12222 图中蓝色的表示结点荷载（已知），红色的表示杆端力（未知的）， 、 分别1、2单元杆端力子矩阵。对1、4结点“荷载”含有未知反力。11 由示意图可见，结构结点的平衡方程为4.1 杆系结构整体分析4.1.3 结点平衡方程的建立1) 一简单例子（如图） 图中有两套编号，红的是单元杆端编号，黑的是结构整体编号。1-1) 结点示意1-2) 结点平衡②①1 杆系结构整体分析 从例图可见，其全部结点平衡方程为若记4.1 杆系结构整体分析若引入矩阵记号，则结点平衡方程可改写作式中[I]、[0]分别为单位和零矩阵。 这一结论虽然是由一个例子得到的，但是显然对一切结构都是成立的。问题在于不同结构，[A]矩阵是不同的。②①③1222114.1 杆系结构整体分析4.1.4 杆端位移用结点位移来表示仍以简单例子来说明若记 由结点、杆端位移的协调条件，可得[? ]、[? ]的对应关系为 式中 [A]T是前面力关系[A]的转置，因此[A]T称为位移转换矩阵。②①③122211 这就是整体刚度方程，它的物理实质是结点平衡。[K]称作结构刚度矩阵（或整体刚度矩阵），[R]称作综合等效结点荷载矩阵，它由两部分组成。4.1 杆系结构整体分析4.1.5 整体刚度方程——结点平衡若记 引入位移转换关系，则②①③12221整体结点码1单元个数4.1 杆系结构整体分析4.1.6 整体刚度矩阵的建立 若将[A]按单元分成图示三个子矩阵 由此可见，整体刚度矩阵可由各单元整体刚度矩阵装配累加得到。为说明如何装配，先将单元刚度矩阵进行分割 则由矩阵乘法可证明，[A]i[k]i[A]iT的结果是，将刚度矩阵子矩阵按整体结点码 r 、s 送整体刚度矩阵相应位置。这一装配规则称为“对号入座”。 则③④①②4.1 杆系结构整体分析1) 任意结构情况 上面结论是通过具体例子（全刚结点平面刚架）得到的，由虚位移原理或势能原理进行整体分析（见讲义），可得任意结构其结论同此例。2) 结点位移编号 对图示三铰刚架，当仅用一种单元（梁柱自由是单元）时结点位移编号如图所示。 如果按结点顺序，对结点非零位移进行依次编号，这一序号称作结点位移码。为便于计算机处理并减少结构刚度矩阵的阶次，将零位移的号码变为零。3) 单元定位向量 按单元局部结点码顺序，将结点位移码排成的向量，称作单元的定位向量。③③④④①①②②③④①②4.1 杆系结构整体分析 对图示刚架各单元的定位向量为如何获得带铰的单元刚度矩阵和等效荷载矩阵①(0,0,1,3,4,5)①(0,0,1,2,3)② (0,0,,6,7,8)② (0,0,2,10,11,12)③(1,2,3,4,5)③(3,4,5,6,7,8)④(4,5,6,7,8)④(6,7,9,10,