【期刊】索膜结构基于非线性有限元法的ANSYS找形分析.pdfVIP

■ - 。碡 ． 鬻 鬻 _ 索膜结构基于非线性有限元法的ANSYS找形分析 郑 昊 ，古娟妮 ，王 娜。 (1．广东工业大学，广东 广州 510090； 2．广东省城市建设高级技工学校，广东 广州510650； 3．华南理工大学土木与交通学院，广东 广州 510640) 摘【 要】 针对索膜结构极强的几何非线性特点，介绍了非线性有限元的基本思想和索膜结构的找形 原理。考虑到 目前有限元通用程序的广泛使用，着重研究了这一方法在ANSYS上的实现过程。 【关键词】 膜结构；找形；非线性有限元法；ANSYS 【中图分类号】 TU383 【文献标识码】 B 膜结构是一种非传统的结构形式，其造型优美，极富现 {△s)=(【B +口 】){AU} 代气息；结构轻巧，有极强的空间跨越能力；易于建造搬迁， 式中{A 伪 单元应力增量矢量；{△8)为单元应变增 有明显的经济效益。因此膜结构建筑成为近几十年蓬勃发 量；D为材料本构矩阵．B 为线性应变位移转换矩阵；BNL为 展的一种新型的大跨度空间形式。 非线性应变位移转换矩阵；{Au}为单元节点位移增量矢量。 膜结构的设计也迥异于传统的刚性结构。由于膜材的 根据虚功原理，推导出非线性位移法找形的UL方程： 柔性特征，其本身没有抗压刚度和抗弯刚度，需赋予一定的 (【I+ M】){hU}={R)一{F} 预拉力才能形成确定的空问曲面形状和抵抗外荷载的能力 ， 式中：{R沩 t+△t时刻的荷载等效节点力矢量，此项在 成为真正的结构。这种在一定预应力作用下，找出的一个既 找形过程中可忽略不计；{F)为t时刻的单元应力节点力矢 符合建筑美观，又满足边界条件和力学平衡的过程就是找形 量；{ 。)为线性应变增量刚度矩阵；{NL)为非线性几何刚 分析过程 、 度矩阵或初应力矩阵。 1 膜结构的非线性有限元法找形 2 ANSYS找形及算例分析 1．1 非线性有限元法简介 2．1 ANSYS找形分析原理 1970年，E ·Haug和 G ·H ·Powell提出了一种基于 ANSYS是一种基于非线性有限元思想的通用有限元软 Newton—Raphson非线性迭代的索膜结构找形方法。它是针 件，可用于索膜结构的找形分形。其基本分析原理是：先用 对索膜结构具有强烈的几何非线性的特点，在小应变、大位 小弹性模量技术，将 目标节点提升到指定高度，用支座移动 移的情况下，首先将膜结构离散成由节点和三角形单元构成 法进行初步找形，目标点固定，其它点连动，得到结构的近似 的空问结构，设定一初始应力分布，采用拉格朗日法建立非 平衡形状。在此几何位形基础上更新节点坐标，释放预应 线性方程组，结合边界条件迭代求解。由于是非线性求解 ， 力 ，重新设定索膜结构的真实材料参数和预应力，进行 自平 无法避免收敛 问题，往往可以通过分段提升的方法来保证收 衡迭代求解。循环若干次，释放掉不平衡力，直至应力分布 敛性。现有的一些通用有限元软件，如 ANSYS，SAP2000， 均匀度达到要求 。 ADIANA等都有用这种方法来进行膜结构的找形。 由于索膜结构中索单元和膜单元均只能承受拉力，故在 1．2 非线性有限元法找形的基本原理