张拉整体结构非线性有限元找形研究.pdfVIP

张拉整体结构非线性有限元找形的研究 关增伟，付朋辉 大连理工大学土木水利学院，辽宁大连 （116023） E-mail ：fuphui@ 摘 要：张拉整体结构是一种新的结构体系，应用在建筑、生物、航天等各个领域。系统 研究了利用通用非线性有限元程序解决张拉整体找形的问题，提出了几种解决非线性有限 元程序不收敛的方法。经算例证明，这几种方法是非常有效的，有重要的工程实用价值。 关键词：张拉整体；找形；非线性有限元 1．引言 张拉整体的概念是由美国著名建筑师富勒提出的。盖格尔在富勒张拉整体穹顶的基础 上，发明了索穹顶，从而使张拉整体的概念首次应用到了大跨度建筑工程中。由于张拉整 体结构固有的符合自然规律的特点，最大限度的利用了材料和截面的特性，因此可以用尽 量少的钢材建造超大跨度的建筑。张拉整体理论上可以做到非常大的跨度，并且随着结构 跨度的增大,其单位面积自重和造价的增加并不大。所以张拉整体在建筑工程中极具前景， 被预测为本世纪最有魅力的结构体系之一[1][2] 。 在给定边界条件下，所施加的预应力的分布和大小同所形成的结构初始形状是相互关 联的。如何合理的确定张拉整体结构的初始平衡状态（即张拉整体结构的找形）是其设计 中的关键技术问题，迄今对这一问题的研究一直是张拉整体结构研究中的焦点。本文的主 要工作就是采用非线性有限元方法求解张拉整体找形这一问题。 2 ．常用找形方法 张拉整体找形的方法有很多，但是都有其局限性，至今还没有一种成熟的方法。其中 力密度法、动力松弛法在国外应用得比较广泛。 2.1 力密度法 力密度法是七十年代 Linkwitz 和 Schek 提出用于索网及膜结构的找形方法，在给定的 结构拓扑条件及边界条件下，力密度法可使找形问题线性化，从而避免了坐标初始值问题 及其它方法致命的收敛性问题[3] 。用力密度法的找形结果在很大程度上依赖于设计人的经 验，平衡形状取决于设计者给定的力密度值，这就在很大程度上限制了它在张拉整体找形 中的应用。 2.2 动力松弛法 动力松弛法早已经成功的应用于膜结构和索网结构，Motro 和 Belkacem 将此通用的找 形方法应用于张拉整体结构[4] 。当结构的节点数较少时，动力松弛法有很好的收敛性，但 是随着节点数的增加收敛效率变差。同样的，如果结构中索杆长比不恒定，此方法将变得 繁琐，因而制约了非规则结构形式方面的应用。 - 1 - 其他还有几何分析法、非线性方程求解法、力的平衡分析法、能量法、坐标缩减法等 [5]，这些方法也都有其局限性，不能广泛使用。 3 ．非线性有限元方法找形理论 非线性有限元在各工程领域都有广泛的应用，同样它也可以应用于张拉整体结构的找 形分析和荷载分析。并且现在有许多通用大型有限元软件，大多工程师对其也比较熟悉， 本文就是寻求采用通用有限元程序来求解张拉整体找形的方法。 3.1 张拉整体结构非线性有限元法基本原理 张拉整体结构是一种柔性结构。在外荷载作用下，结构的变形具有小应变、大位移的 特点，属于几何非线性问题，通常采用修正的拉格朗日描述（U.L.描述），基于虚位移原理 建立张拉整体结构的非线性有限元方程。以 t 时刻的结构形状作为参考构形，则 t 时刻结 构形状和位移增量由以下方程描述 {} { } x [N ]X e （1） { } { } ∆u [N ]∆d e （2 ） 式中，{} { } x 、 ∆u 为单元内任意一点坐标和位移增量向量，[N ]为单元的形函数矩阵，