(钢结构设计)3_4空间杆系有限元法.ppt

空间杆系有限元法也称空间桁架位移法。 空间杆系有限元法是计算精度最高的一种方法，适用于各种类型、各种平面形状、不同边界条件的网架，静力荷载、地震作用、温度应力等工况均可计算。 能考虑网架与下部支承结构的共同工作 。 计算程序见下表。 3.4.1网架计算基本假定 网架的节点为空间铰接节点，杆件只承受轴力； 结构材料为完全弹性，在荷载作用下网架变形很小，符合小变形理论。 3.4.2单元刚度矩阵 一等截面空间桁架杆件ij如图所示，设局部直角坐标系为 ， 轴与ij杆平行。 杆端力向量为： 杆端位移向量为： 杆端力和位移的关系可写为 结构分析中为方便杆端力和位移的叠加，应采用统一坐标系，即结构整体坐标xyz。这样需对局部坐标系下的单元刚度矩阵进行坐标转换。 设杆件ij （即 轴）与整体坐标x，y，z轴夹角的余弦分别为l，m，n。由图25所示的几何关系可以得出 式中lij——ij杆的长度 令 分别表示杆件ij在整体坐标系中的节点力，节点位移和单元刚度矩阵。 在整体坐标系中ij杆节点力和节点位移间的关系力为： 两坐标系之间的转换关系为 式中[T]——坐标转换矩阵 并注意到[T]-1=[T]T，得到整体坐标下ij杆节点力和位移的关系为： 得到杆件ij在整体坐标系中的单刚矩阵 ： 3.4.3结构总刚度矩阵及总刚度方程 建立了杆件单元刚度矩阵之后，即可按照变形协调及节点内外力平衡条件建立结构的总刚度矩阵及相应的总刚度方程。 对公式变换为： {Fi} ，{Fj}——分别为杆件ij在整体坐标系下i，j点的杆端力列阵； {δi}，{δj}——分别为杆件ij在整体坐标系下i，j点的位移列阵； [Kii]，[Kjj]——分别为杆件ij在i端，j端发生单位位移时，在i端，j端产生的内力； [Kij]，[Kji]——分别为杆件ij在j端，i端发生单位位移时，在i端，j端产生的内力。 以图26所示的空间桁架节点 3 为例，说明总刚矩阵及总刚方程的建立。该桁架共有9个单元，5个节点，单元及节点编号如图示。相交于节点3的杆件有⑥⑦⑧⑨。 变形协调条件为连于同一节点上的杆端位移相等 ，即： 内外力平衡条件为汇交于同一节点的杆端内力之和等于该节点上的外荷载，即： 连于节点3的杆端力与各节点位移关系为： 整理得： 上式就是节点3得内外力平衡方程，对网架中得所有节点，逐点列出平衡方程，联立起来便为结构总刚度方程，表达式为： 对于本例，总刚度矩阵中的第7行至第9行的元素表示如下： 总刚矩阵具有下列特点： 矩阵具有对称性，计算时不必将所有元素列出，只列出上三角或下三角即可。 矩阵具有稀疏性，网架结构每一节点所连杆件数量有限，总刚矩阵中除主对角及其附近元素为非零元素外，其余均为零元素。 3.4.4总刚矩阵中边界条件的处理方法 未引入边界条件前，总刚矩阵[K]是奇异的，不能进行求解。引入结构边界条件消除刚体位移后，总刚矩阵为正定矩阵。 3.4.6 杆件内力 引入边界条件后，求解公式，得出各节点的位移值，由公式和公式可得出ij杆端内力为 将公式展开并代入公式整理可得杆件内力表达式为 式中 N——杆件轴力，以拉为正。 3.4.7 空间杆系有限元法计算步骤 (1)根据网架结构、荷载对称性选取计算简图，并对其节点和杆件进行编号，为减小总刚矩阵带宽，节点编号应遵循相邻节点号差最小的原则。 (2)计算杆件单元长度及杆件与整体坐标轴夹角余弦； (3)初选各杆的截面积； (4)建立局部和整体坐标系下的单元刚度矩阵； (5)集合总刚矩阵，为减小矩阵容量，宜采用变带宽一维存贮方式； (6)建立荷载列阵； (7)引入边界条件对总刚度方程进行处理； (8)求解总刚度方程，得出各节点位移值； (9)根据节点位移计算杆件内力； (10)按杆件内力调整杆件截面，并重新计算，迭代次数宜不超过4～5次。 * * 3.4空间杆系有限元法 网架杆件 节点位移 单元刚度矩阵 总刚度矩阵 总刚度方程 节点位移值 杆件内力 单元内力与节点位移间关系 引入边界条件 节点平衡及变形协调条件 基本单元 基本未知量 奥运会场馆 鸟巢 图3.24 ij杆的杆端轴力和位移 局部直角坐标下 图3.25 杆件在整体坐标中 整体坐标 坐标转换 奥运会场所 图3.26 单元及节点编号 位移为零 弹性约束 指定位移 处理方法 * * *