张力结构的非线性有限元分析Ξ.pdfVIP

第 14 卷第 3 期 计 算 力 学 学 报 V o l. 14 N o. 3 1997 年 8 月 CH IN ESE JOU RNAL O F COM PU TA TLONAL M ECHAN ICS A ugust 1997 张力结构的非线性有限元分析 　　　　董　明　　　夏绍华　　　　　钱若军　　沈祖炎　　　 ( ) ( ) ( ) 　　　　 同济大学, 上海, 200092 　 河海大学, 南京, 210024 　　　　　　　 同济大学, 上海, 200092 摘　要　张力结构中的索结构分析通常采用非线性二节点直线或二次曲线单元, 但是在大 跨度索结构, 尤其是索穹顶结构的分析中, 常用的单元已不能满足精度的要求。本文提出了 考虑自重作用下有初始垂度的五结点非线性空间曲线元模型, 放弃了一些特殊的假定, 考虑 了应变表达式中高阶量的影响, 推导出了适合于弹性大位移几何非线性分析的L agrange 方 程及其切线刚度矩阵, 编制了相应的计算程序。若干数值例题验证了本文理论的正确性, 数 值解具有极高的精度。本文提出的模型可适合于任意索结构的分析。 关键词　索穹顶结构; 几何非线性; 非线性有限元 分类号　TU 394; O 242 21 1　引　言 张力结构的非线性有限元分析无论在国内外都均已进行了大量的研究, 经过长期的研究 和比较, 在工程中已被广泛采用的预应力索网、索桁架及单索体系等张力结构的有限元分析方 法也已趋于成熟, 在各种不同的数学力学模式中, 较多地被采用的是二节点一维非线性空间直 线元, 文献[1 ]、[5 ]等均给出了相关的有限元基本方程及单元刚度矩阵的表达式, 这种模式是 基于如下的假定: 即不计索元的自重, 索元是直线元。这个假定对索段较短的索网、索桁架等结 构是合适的, 然而, 随着大跨度张力结构的发展, 尤其是索穹顶结构的出现, 对于结构中较长的 索段如索穹顶中的脊索和环索等, 上述的假定显然会引起不能容忍的误差。因此, 有必要充分 考虑索元自重导致的初垂度及初应变和初应力的影响, 进一步研究更为精确的数学力学模式, 才能正确地进行形状分析(FORM F IND IN G) , 预应力态和荷载态分析。 本文提出了一种迄今为止最为精确的五结点非线性空间曲线铰结单元模型, 由于在建模 时放弃了一些特殊的假定, 因此这种模型便更具有一般性。数字算例表明, 本文提出的模式是 精确, 具有极高的精度, 并有较快的收敛速度和较好的稳定性。故可用于索穹顶结构及其它任 意张力结构的分析。 2　基本假定和坐标系的建立 2 1　基本假定 考虑到张力结构的力学特性, 本文采用了如下假定: 国家自然科学基金资助 　收稿日期: 　董明: 男, 1973 年生, 硕士研究生 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 3 期　　 　　　　　　　　董　明等: 张力结构的非线性有限元分析　　　　　　　　　　 269 1 索之间的交叉节点为无摩擦的理想铰接节点; 2 索被认为是完全柔性的, 即既不能承受任何弯矩, 也不能承受任何压力; 3 索在受拉时材料仍处于线弹性阶段, 即应力应变关系仍满足虎克定律; 4 考虑索段自重的影响。 2 2　坐标系的建立 如 同一般, 张力结构是在结构的整体坐标系O X Y Z 中定义, 而对于各个单元可分别建立 单元局部坐标系 o ′x y z 。整体坐标系和单元局部坐标系均为迪卡尔右手直角坐标系, 如图 1 所 示。而沿单元弧向可建立弧向坐标系S 。单元节点的位移向量、荷载向量都是在弧向坐标系S 中