索膜结构中滑移索地研究.pdfVIP

第四届全国土木工程研究生学术论坛论文集 31 索膜结构中滑移索的研究 刘 凯1，高维成2，刘宗仁1 (1．哈尔滨工业大学土木工程学院，哈尔滨150090l2．哈尔滨工业大学航天科学与力学系．哈尔滨150001) 摘要：针对索膜施工过程中索在膜套中(或膜表面上)滑移的现象．提出了目前索膜结构分析软件均未考虑 索在膜套中(或在膜表面上)滑移的问题。采用非线性有限元的方法并利用线弹性模型和非线性滑移刚度模型 建立了一种适合于索膜结构的新的滑移索单元。本文首先对现有的滑移索单元计算方法进行了简单介绍，指 出了现在对滑移索数值模型存在的一些问题，其次利用空间两节点直线索单元原长与端部的关系推导了滑 移索的滑移刚度；然后提出了采用修正的线性滑移刚度法迭代求解索在膜套中(或在膜表面上)滑移的算法， 最后通过算例验证了本文所述程序的正确性。 关键词：索膜结构；滑移索I施工 自20世纪90年代开始，张力结构的发展呈现出一些新的特点，把索与膜结合起来形成的柔性张力 结构——索膜结构，获得了快速发展。钢索是膜结构重要结构性构件之一，特别是张拉膜结构[1]，它是 通过张拉钢索使薄膜结构获得刚度。在实际工程中．张拉膜柔性边界节点[1](主要指膜与柔性索连接的 节点)采用最广泛的是膜套节点。当膜套中索的张力值大于索与膜套之间的摩擦力时，索就会在膜套中 滑移，此时必须采用滑移索单元对索膜结构进行分析，但现在所有的索膜结构分析软件均未考虑索在膜 套中或在膜表面上滑移的问题，这会导致设计索力不均等一系列问题。由于非线性有限元方法在结构 分析中被广泛采用，所以本文基于此方法，采用线弹性模型和非线性滑移刚度模型建立了一种适合于索 膜结构的新的滑移索[2]单元。 1 滑移索单元计算方法简介 因滑移索跟离散索相比，不仅要满足节点平衡条件，还要满足相容性条件，所以常规的有限元方法 不便求解。目前，对于滑移索的数值模型主要有三种方法：a．线弹性模型Is-s]。这种模型采用线性理论 求解，较为简单，但由于索结构的变形呈现高度的几何非线性，所以线弹性模型精度不高，使用范围受到 限制。b．非线性模型[6。7]。这种模型考虑了滑移索的非线性，求解精度较高，但计算过程相当复杂，不便 于应用。c．其他模型陋引。 2 索膜结构滑移索数值模型 2．1基本假设 a．索是理想柔性的，且满足虎克定律-b．节点为理想无摩擦的铰接节点；c．索在膜套中或膜表面 上是无摩擦滑移。 2．2索单元的基本方程 由于索一膜结构的曲率一般不是很大，当索单元长度较短时，用直线单元来模拟不会产生很大的误 差，并且计算简单；同时若欲提高精度也可以通过加密单元来实现啪，因此结构分析时采用了两节点直 线索单元。假设索单元内任一点的位移厂一[“口硼]r可用单元的节点位移以一 [“， u。 叫。 “。 可： 7．PJz]r和该点的坐标X。一ky z]丁来表示，则该单元的位移函数为： 收稿日期：2006-03—28 修改稿日期：2006-06-16 基金项目：国家自然科学基金资助项目 作者简介：刘凯(1978～)，男，河北省献县人，博士研究生．主要从事大跨施工方面的研究． 32 第四届全国土木工程研究生学术论坛论文集 厂=[N]乱 FNl0 0 N2 0 0] 0 0 N。0 N：0 式中：EN]一l l为索单元的形函数矩阵，其中N，=垒二!，N：=三。在局 0 0 N1 0 0 N2 I 部坐标系中，z。即为索段长度￡。每一迭代步内节点位移增量为： 允一[△l‘lAvlAwlAu2加2△训2]丁 (1) 2．3 单元原长与端部张力的关系 索单元张力：丁=To+AT—To-t--E地 (2) 式中To为索元的初