空间梁单元切线刚度矩阵精确分析方法.docVIP

空间梁单元切线刚度矩阵精确分析方法 　　摘要：为有效进行空间刚架结构后屈曲分析，提出一种新的空间梁单元切线刚度矩阵的精确分析方法。首先用直接法建立梁单元杆端力与杆端位移增量的关系式，然后根据矩阵微分理论求单元杆端力关于杆端位移的导数，最后在求导结果表达式中令杆端位移增量为0，即可得到梁单元切线刚度矩阵。对六层和二十层空间刚架结构进行了后屈曲分析。结果表明：所得的空间梁单元切线刚度矩阵具有足够精度，可有效用于大型空间刚架结构的后屈曲分析。 　　关键词：空间梁单元；切线刚度矩阵；空间刚架；非线性屈曲；弹塑性屈曲 　　中图分类号：TU393.3 文献标志码：A 　　0 引 言 　　在空间刚架结构的屈曲分析中，单元切线刚度矩阵的预测精度对分析结果的正确性有着关键性影响，特别是对于结构的后屈曲分析。若单元切线刚度矩阵预测精度不够，则很难通过荷载-位移平衡路径的极限点，因此也很难准确地确定结构的极限荷载。建立空间梁单元切线刚度矩阵有多种途径，如势能变分方法[1-2]、Timosheko梁柱理论[3]、弹性刚度矩阵与几何刚度矩阵的叠加方法[4]以及一些其他方法[5-7]。 　　势能变分方法[1]建立切线刚度矩阵常常要忽略一些高阶项，导致切线刚度矩阵的预测精度受到影响。Timosheko梁柱理论直接由平衡方程来导出切线刚度矩阵，平衡方程中的力与位移关系可用超越函数表示，分析方法简单，切线刚度矩阵更精确。 　　Yang等[4]首先采用虚功方法来建立梁单元的几何刚度矩阵，然后由弹性刚度矩阵与几何刚度矩阵的叠加建立切线刚度矩阵。文献[5]中基于Prandtl-Renss增量理论，推导出三维梁单元的弹塑性切线刚度矩阵。文献[6]中基于非线性问题的一般平衡方程和空间梁单元的非线性几何方程，推导应力-应变一般线弹性关系下的空间梁单元显式切线刚度矩阵，该刚度矩阵中包含了由初应力和初应变产生的初应力刚度矩阵。文献[7]中使用复合变量求导方法对内力近似求导得到切线刚度矩阵。上述这些方法所建立的空间梁单元切线刚度矩阵并没有通过大型算例来验证其有效性。 　　为有效进行空间刚架结构后屈曲分析，本文中提出一种新的空间梁单元切线刚度矩阵的分析方法。首先在变形状态下用直接法建立梁单元杆端力与杆端位移增量的关系式，然后根据矩阵微分理论求单元杆端力关于杆端位移增量的导数，最后令杆端位移增量为0，即可得到梁单元切线刚度矩阵。本文中笔者建立切线刚度矩阵时，没有忽略任何高阶项，因此所建立的切线刚度矩阵是精确的。通过六层和二十层空间刚架结构的后屈曲分析，验证了本文空间梁单元切线刚度矩阵的有效性。 　　1 空间梁单元杆端向量变换矩阵 　　空间刚架结构的精确非线性分析，需要利用节点空间大转动分析理论[8-10]。节点空间大转动情况的转动次序不满足交换律，即节点总量转角不能由转角增量线性叠加累积计算。 　　节点空间大转动分析需要建立3种局部坐标系：节点坐标系、单元坐标系和单元端面坐标系。这3种局部坐标系对整体坐标系的旋转矩阵分别称为节点方向矩阵、单元方向矩阵和单元端面方向矩阵。图1中给出了结构变形构形下空间梁单元的单元坐标系、单元端面坐标系和结构坐标系。 　　1.1 节点方向矩阵 　　节点坐标系是一个刚性固定于节点上的笛卡儿坐标系，节点发生转动时，该坐标系方向（节点方向）随之变化。通常令初始构形下的节点坐标系与整体坐标系平行，故初始节点方向矩阵q0为3阶单位矩阵，q0=I3，I3为3阶单位矩阵。节点发生转动后，节点方向矩阵q按照节点空间大转动的坐标旋转变换矩阵进行更新，即 　　式中：θi为转角矩阵θ的元素，i=1，2，3；θx，θy，θz分别为x，y，z方向的转角。 　　1.2 单元方向矩阵 　　单元坐标系以Oxeyeze表示（图1），单元坐标系的xe轴通过单元j端和k端截面形心。在初始构形下，单元坐标系的ye，ze轴分别与单元截面2个主轴方向平行。单元方向矩阵由单元坐标系的xe，ye，ze轴对整体坐标系x，y，z轴的方向余弦向量构成。初始构形下的单元方向矩阵r0表示为 　　式中：r01，r02，r03分别为初始构形下的单元坐标系xe，ye，ze轴对整体坐标系x，y，z轴的方向余弦向量。 　　在变形构形下，单元截面形心轴发生弯曲，单元2个端截面也发生了转动和位移，单元坐标系的方向发生了变化。在变形构形下，单元坐标系的xe轴对整体坐标系x，y，z轴的方向余弦向量r1为 　　式中：u为整体坐标系下单元j端和k端节点线位移集合向量，u=（ujx，ujy，ujz，ukx，uky，ukz）T；x为增量步开始时单元j端和k端节点坐标的集合向量；t1=（-I3，I3）。 　　变形构形下的单元方向矩阵r可表示为 　　式中：r1，r2，r3分别为变形构形下单元坐标