两种索网结构找形分析的方法及其比较.docVIP

两种索网结构找形分析的方法及其比较 　　摘要：针对大型空间索网结构分析中的首要问题——索网结构的找形分析，本文基于文献[1]提出的四节点等参曲线索单元模型理论，编制了相应程序进行了实际结构找形分析计算；并对经典的找形分析方法——动力松弛法提出改进措施，同样编制了相应程序进行了实例计算。两种方法都和大型通用有限元分析软件ANSYS计算结果进行对比，证明了理论的正确性，为把这两种方法推广到索网结构的静力和动力分析做好了准备。 　　关键词：索网结构，找形分析，四节点等参元，动力松弛法 　　中图分类号： TU3 文献标识码： A 文章编号： 　　近年来随着社会经济的发展，索网结构作为一种美观、经济、实用的结构形式得到了广泛的应用。与普通的混凝土结构相比，索网结构没有初始刚度，结构需要通过预张拉实现其初始刚度的建立，从而承受外荷载作用。在边界约束条件和结构几何条件确定的前提下，确定一个和结构预应力条件相适应的几何构形，就构成了空间柔性结构中的初始平衡问题，即初始形态分析问题或称为找形分析问题[2]。 　　利用四节点等参曲线索单元模型理论进行找形分析 　　现今索网结构找形分析方法主要集中在有限元方法上。各种有限元方法的主要区别在于其单元类型的不同，有二节点直杆单元、三节点等参元以及五节点等参元等[3][4][5]。文献[1]中提出了一种新型的单元类型：四节点等参曲线索单元模型，并推导了该模型比较完整的相关理论，采用该单元模型可以取得计算精确度和计算代价两者较好的平衡。 　　1.1 四节点等参曲线索单元模型理论 　　这里给出该理论的要点，详细的理论推导过程可以参见文献[1]： 　　理论推导基于如下假定： 　　⑴ 索始终处于弹性工作阶段，应力应变关系符合虎克定律； 　　⑵ 索是理想柔性的，只能承受拉力，不能承受压力和弯矩； 　　⑶ 大位移小应变假定； 　　依照上述假定进行公式推导，可以得到轴向应变的表达式为：,其中，，。 　　应力表达式为,其中为索单元的轴向应力；为索材料的弹性模量；为索单元的轴向应变；为索单元的初始轴向应力。 　　根据虚功原理可推得全量形式的单元平衡方程，然后将全量形式的单元平衡方程进行微分，得到增量形式的微分方程： 。其中，为单元的切线刚度矩阵，的具体表达式： 　　 　　 　　 　　1.2 四节点等参曲线索单元模型理论找形分析 　　利用四节点等参元模型理论进行索网结构找形分析是一个典型的有限元分析过程，所以具体的计算步骤应遵循传统的有限元分析步骤。这些步骤简而言之就是划分单元、确定加载步然后在每一加载步下循环直到位移值满足精度要求为止。需要强调的是找形分析时结构的初始状态是零状态，求解结果是结构找形完成时的状态，结构外荷载则是边界点的位移值，分步加载就相当于分步移动边界节点到指定位置，这是找形分析和其他索网分析过程例如静力分析的不同点。 　　利用动力松弛法进行找形分析 　　对于柔性结构的形状确定，动力松弛法最大的优点是能够由任意初始非平衡状态得到平衡状态，并且与其他方法如非线性有限元法相比，该方法不需要组装单元刚度矩阵，计算简单，容易取得较好的稳定性和收敛性，因此动力松弛法相对而言具有其自身的优势，适合于索网结构的找形分析。 　　2.1 动力松弛法基本公式推导 　　根据牛顿第二运动定律，时刻节点在方向的运动可以表示为：，用中心差分法表示该式为：，则时刻节点在方向的坐标改变为：，同理可得到节点在方向的速度分量和坐标改变。 　　根据当前结构几何位置可以求得任意索段的内力：其中，，分别为索段初始状态内力、弹性刚度和长度。 　　根据上述求得的索段内力，可以求得新的节点不平衡力如下： 　　 　　其中求和符号表示对汇交于节点的所有连杆求和，与节点相连的连杆的另一端为，连杆内力为，为杆长。表示作用于节点的方向外荷载。 　　根据新的节点不平衡力可求得新的速度分量及新的结构坐标，迭代过程如此继续。 　　在此过程中，按照下式计算结构的动能：。当结构的动能达到极大值时，所有的速度分量被设置为零，运动从更新了的当前几何位置重新开始运动，如此反复直到结构的动能和节点不平衡力足够小为止，此时认为结构已经达到我们期望的平衡位置，找形过程结束。 　　2.2 保证动力松弛法收敛的改进措施 　　在上述动力松弛法公式中，时间增量、节点质量、节点速度和节点坐标等都可以不是真正物理意义上的量值，通常被称作虚量。我们可以从数值计算的角度，并从稳定性、收敛性和精确性三个方面来研究动力松弛法参数的选取。 　　① 关于时间增量 　　通常的做法是为了简化计算，将取为1，代入节点速度公式可得，继续求得新的结构坐标，这样迭代继续进行。 　　但是我们可以把取为更小的值，即将1s分成若干等份，比如5份或者10份