第四章杆系结构静力分析程序编制原理.docVIP

第四章 杆系结构静力分析程序编制原理 4.1 杆系有限元分析程序流程 根据第二章有限元分析的过程，可以得到图4-1所示的典型杆系有限元分析基本程序流程图。 图4-1 有限元分析基本流程图 考虑到桥梁结构有限元分析一般要按照施工过程进行分析，有限元分析基本流程图中的若干功能模块必须结合施工仿真要求做局部调整。本章的论述中给出的流程图是结合施工分析程序而绘制的，是经过调整后的桥梁有限元施工分析程序流程，如图4-2所示。为了便于解释，对程序模块进行了编号，以表示各级子程序的主要调用关系。编号中每两位表示一个级别，例如：6号子程序调用601号子程序，601号子程序调用60101号子程序。 本章仅讨论图4-1中的有限元分析基本模块，图4-2中与施工有关的程序模块将在第五章中再作论述。 4.2 空间杆件系统单元坐标系内的单元刚度矩阵 空间杆件同时承受轴力、弯矩、剪力和扭矩的作用，对两节点空间梁单元，每个节点有6个自由度，杆端位移和内力向量为： (4-1) (4-2) 在右手坐标系中，它们的正方向遵从右手法则，如图4-3、4-4所示。 图4-3 杆端位移 图4-2 按施工过程进行桥梁结构分析的有限元程序 图4-4 杆端内力 对于线性问题，对杆单元、扭转单元、梁单元三种单元的单元刚度矩阵进行组合，可以得到空间梁单元的单元刚度矩阵： (4-3) 4.3 转换矩阵 4.3.1转换矩阵的概念 一般情况下单元坐标系与结构坐标系是不一致的，而节点平衡方程是在结构坐标系下建立的，这样就需要将单元坐标系内的杆端力、节点位移和单元刚度矩阵转换到结构坐标系内。两者之间的关系通过转换矩阵表示。 设在单元坐标系中的力学量以右上角撇号标出，单元坐标系与结构坐标系中的节点平衡方程分别为： (4-4) (4-5) 任意向量在单元坐标系与结构坐标系中分量之间的关系为： (4-6) 式中：，为单元坐标系坐标轴单位方向向量在结构坐标系中的方向余弦： 对于具有12个自由度的和以及和，有如下变换关系： (4-7) (4-8) 式中：，容易验证具有性质。 综合式(4-4)、(4-5)、(4-7)、(4-8)可以得到结构坐标系下的单元刚度矩阵： (4-9) 式(4-9)表明了单元坐标系与结构坐标系中的单元刚度矩阵之间的关系，称为转换矩阵。 4.3.2转换矩阵的计算 空间梁单元的单元坐标系原点在节点，轴取梁中心线，从节点指向节点为正方向，如图4-5所示。轴的单位方向向量为： (4-10) 式中：——单元的长度。 图4-5 单元坐标系坐标轴方向向量计算示意图 和轴应沿断面的惯性主轴方向。一般指定一点，通过、、三点来确定平面，同时指明轴的正方向指向点所在的一侧。注意点不能在杆件轴线及其延长线上。轴的方位及正方向由右手法则确定。记向量 (4-11) 作向量 (4-12) 显然向量指向轴的正方向，于是轴的单位方向向量为： (4-13) 式中：是向量的模。 最后可以得到轴的单位方向向量为： (4-14) 得到的局部坐标轴单位方向向量、、的各分量即为矩阵中的方向余弦，进而可以得到转换矩阵。 形成转换矩阵的程序流程见图4-6。 图4-6 计算转换矩阵的子程序流程图 4.4 带刚臂单元 使用杆系模式计算桥梁结构时，理论上的结构简化计算图式往往会与实际的结构物不相符合，主要有两方面原因：一是节点上所有杆件的轴线未必会交于同一点；二是杆件进入节点附近时常常和刚性很大的节点块连在一起。处理这样的问题一般将该区域假定成刚性域，用带刚臂单元来近似代替。下图为需要采用带刚臂单元的一个例子。 图4-7 需要采用带刚臂单元的情况 通过对一般的单元刚度矩阵作相应变换，可以在程序中引入带刚臂单元。 图4-8为一个带刚臂单元，左右节点号分别为、，刚臂与直杆的连结点为、。单元坐标系由节点、、确定，另外设立“真实”单元坐标系，由节点、、确定，点在、的连线或延长线之外，并且在“真实”单元的主平面内。在坐标系内，单元两端的刚臂长度为。 图4-8 带刚臂单元 由于刚臂本身没有变形，在坐标系中、点与、点的位移关系为： (4-15) (4-16) 令由坐标系到坐标系的转换矩阵为，则坐标系内直杆的位移为： (4-17) 注意不满足的关系。 以杆件实际内力状态作为力状态，在坐标系内、点发生位移作为第一种虚位移状态，对应的杆端力为；在坐标系内、点发生位移作为第二种虚位移状态，对应的杆端力为。由于刚臂段不发生任何变形，所以由虚功原理可以得到： (4-18) 结合式(4-17)有 (4-19) 令杆件在坐标系中的单元刚度矩阵为，则 (4-20) 将式（4-17）、（4-20）代入式（4-19）有： (4-21