膜结构初始形态确定有限元分析.docVIP

膜结构初始形态确定有限元分析 　　摘要：介绍索膜结构的特征和基本体系，寻找初始平衡态，确定初始形态的非线性有限元方法以及原理，建立相应的几何方程、物理方程和平衡方程；应用一个典型的膜结构工程作为算例分析，说明膜结构形态确定分析在ANSYS程序中实现的一般过程，最终确定膜结构的初始形态。 　　 关键词：膜结构；初始形态；非线性有限元；ANSYS 　　 膜建筑是20世纪70年代发展university起来的一种新型建筑体系，是空间结构的重要组成部分。柔性的膜结构的结构组成材料本身没有受压性能，只能通过施加预应力，使膜或者加强索获得必要的张拉刚度，从而形成抵抗外部荷载的结构抗力。基于这点，结构形态确定就是膜结构设计过程很重要一个环节。在给定的边界条件下，所施加的预应力的分布和大小与所形成的结构初始形状是相互联系的，合理的确定这一初始形状和相应的自平衡预应力，就是膜结构的初始形态分析。膜结构的初始形态即结构正常使用时的形态，忽略其自身的微小的自重和自平衡预张力，结构不承受任何外部荷载，此时也成为结构零态。 　　 一、膜结构及其形态确定理论概述 　　 （一）非线性有限元法找形的基本原理 　　 非线性有限元法找形的基本原理先将索膜结构离散成若干单元，然后针对索膜结构的小应变、大位移特点，应用几何非线性理论，建立以节点位移为基本未知量的非线性有限元方程组，最后用迭代计算方法并结合边界条件求解。 　　 （二）平衡形态分析及基本假定 　　 由于初始平衡状态是纯力学平衡问题，与所采用的材料无关，故在计算过程中采用小弹性模量法，以便使结构自由变形，达到平衡。先假定如下： 　　 1．索离散化为空间铰接两结点杆元，膜体离散化为3结点三角形平面应力单元； 　　 2．索膜之间没有相对滑动； 　　 3．荷载等效化为结点力； 　　 4．索膜张拉为小应变，符合HOOKE定律； 　　 5．膜材为正交异性弹性材料，变形前后材料主轴始终保持垂直。 　　 （三）小弹性模量曲面自平衡迭代法分析 　　 采用小弹性模量曲面自平衡迭代法，在实际计算中常常把材料弹性模量取得很小，为实际值的1%~0.1%，采用小弹性模量法使得初始位移向目标位移过渡过程中产生的附加应力很小，可以忽略不计，因此，最终得到的目标曲面便可以保持初始设定的预应力状态。同时，在非线性求解过程中，小弹性模量对加强收敛也有较明显的作用。 　　 采用三角形膜单元为例，非线性有限元法的基本公式： 　　 忽略材料非线性，只考虑几何线性： 　　 {Δσ }=D{△ε}{△ε}=（[BL+BNL]{ΔU} 　　 式中{Δσ}单元应力增量矢量；{△ε}为单元应变增量；D为材料本构矩阵．BL 为线性应变位移转换矩阵；BNL为非线性应变位移转换矩阵；{Δu}为单元节点位移增量矢量。 　　 时步n的应力则表示为：σn=σn-1+Δσn-1 　　 根据虚功原理，推导出非线性位移法找形的UL方程: 　　 ([KL+KNL]){ΔU}={R}-{F} 　　 式中：{R}为t+△t时刻的荷载等效节点力矢量，此项在找形过程中可忽略不计；{F}为t时刻的单元应力节点力矢量；{KL}为线性应变增量刚度矩阵；{KNL}为非线性几何刚度矩阵或初应力矩阵。 　　 二、ANSYS方法的膜结构形态确定 　　 ANSYS的基本思想是属于非线性有限元法。以如下膜结构算例阐述分析步骤： 　　 结构外形为正方形，对角线距离为10m,高度为4m,膜面的初始预张力σ=20N/cm,张拉刚度Et=2550N/cm,剪切刚度Gt=800N/cm,泊松比γ=0.3.结构4角点固定，4边为柔性索边界，边索的初始预拉力均为30KN,EA=3*104KN。 　　 膜结构形态确定的实现步骤： 　　 1．建立几何模型。一般由底向上建模，即先建关键点，然后连线，再成面。 　　 2．选择单元类型，输入实常数。本文采用索单元(Iink10)和薄膜单元(Shell41)，且都只受拉；输入索和膜预应力，索用初应变，膜用降低温度法，索的面积，膜的厚度及温度膨胀系数，索膜的虚拟弹性模量，为真实值的10-3，进行网格划分。 　　 3．设定边界条件、各个支座控制点的提升位移，进行第一次找形，保留第一次找形的几何信息。 　　 4．将各个边界控制点的支座位移重新设为零，输人材料的真实弹性模量，重新设定结构索和膜的真实预应力状态，进行第二次找形。 　　 5．重复（4）若干次，直到结构自平衡迭代求解的响应结果达到设定精度。 　　 6．求解结束，最终形态就是找形求解的结果。 　　 　　 图1：自平衡后迭代形态 　　 　　 　　 　　 图2：第一次找形应力分布 　　 　　 　　 图3：最终自平衡迭代应力分布 　　 表1ans