充气膜结构研究进展.docVIP

充气膜结构研究进展 　　摘 要：本文从充气膜结构的结构设计原理入手，综述了其形态分析、荷载分析、剪裁分析等方面的研究现状与发展方向。 　　关键字：充气膜结构;形态分析;荷载分析;剪裁分析 　　 　　充气膜结构是以性能优良的薄膜为材料，通过向薄膜构成的密闭空间内充气，利用空气压力支撑膜面，从而形成具有一定刚度、能够覆盖大跨度空间的结构体系。 　　由于膜材所特有的非线性力学特点以及膜结构整体所表现的柔性、张力与形态的统一性，其结构设计原理显著区别于传统结构，属于大形变条件下应变和应力问题[1]。主要包括四个阶段：方案设计、形态分析、荷载分析、剪裁分析。其中，找形分析是基础，荷载分析是关键，剪裁分析是目标和归宿。有关充气膜结构的主要研究工作也就集中在这三者之上[2-4]。 　　 　　1形态分析 　　 　　又称找形分析、找形，目的是寻找满足边界条件和初应力平衡条件的结构形状。初始平衡态的寻找是形态分析的关键，力密度法、动力松弛法和非线性有限元法是索膜结构初始形态分析的主要方法。其中，非线性有限元法在我国相关领域内应用最为广泛。 　　陆鉴恒等人[5]针对膜结构找形中最小曲面的确定问题，采用动力松弛法，对迭代参数进行分析和简化，使迭代参数的简化只跟时间步长有关。从算例数据可得出，在收敛范围内，迭代次数n随着迭代步时间步Δt的增加大体呈先减少再增加的趋势，最小值在T/4附近。并发现：a.动态阻尼动力松弛法的两个参数是相互联系的，跟每一时间步质点对应的周期有关；b.参数的取值：虚拟质量为任意常数，时间步长与对应时刻的质点周期对应，取值范围为(0，T/π)，建议取T/4左右；c.此方法简化了参数的选择，明确了参数选择的物理意义。简化虽然增加了迭代的次数，但是在可接受的范围内，且误差比较说明提出的方法计算精度高，结果可靠，值得尝试和进一步研究改进。 　　东南大学的周树路等人[6] 则针对力密度法的找形过程进行改进，避开其中“力密度”的概念，直接引入膜面应力和索拉力作为初始条件，以节点不平衡力作为控制误差，避免了传统力密度法需要反复试算力密度取值的弊端，使找形计算过程简洁高效。据此编制找形程序，通过复杂算例验证了该算法的正确性和普适性。 　　鉴于力密度法原理简单但找形结果往往不能满足精度要求；非线性有限元法结果精度高但存在确定初始坐标问题和非线性系统的收敛问题。针对这两种方法的不足，温世峰等人[7]在综合以上两种方法后得到了混合法对膜结构进行找形。从两个实际工程的找形结果来看要比力密度法更精确些，在找形过程上又避免了非线性有限元法建模的繁琐和不收敛问题，而且还加快了平衡迭代的速度，是一种可以有效提高找形精确度和找形速度的简便而切实可行的方法。表1为三种方法的找形结果对比。 　　 　　 　　2荷载分析 　　 　　荷载分析是指在己求得预应力作用下稳定的结构初始形状的基础上，确定膜结构在荷载作用下的应力和位移是否满足要求。膜结构的控制荷载通常是风荷载和雪荷载，在荷载作用下膜结构的变形较大，且随着形状的改变，荷载分布也在改变，因此要精确计算结构的变形和应力。 　　在找形分析所确定的初始状态基础上考虑自然和人为的荷载作用。结构静力分析考虑的主要荷载是风荷载、雪荷载、自重、温度荷载及悬挂荷载。前面所述的找形分析的三类方法都可以用于索膜结构的静力分析，其中非线性有限元法是最常用的方法。非线性有限元的基本思想已在上文中详细阐述，这里不再重述。静力分析非线性方程组的求解主要有增量法和迭代法两种。 　　膜结构作为一种风敏感型结构，因此流体-结构耦合作用现象在其风致动力响应方面表现非常明显。刘振华等人[8]针对膜结构本身的特点，给出了流体-结构耦合作用下流体的基本控制方程的ALE(Arbitrary Lagrangian-Eulerian)描述、结构的基本控制方程和求解耦合方程的直接求解算法。采用Adina软件实现了简单张拉膜结构的流体-结构耦合作用风致动力响应数值模拟，并与已有风洞试验结果进行了比较。对比结果表明，给出的方法可以准确地计算出柔性结构考虑流体-结构耦合作用的风致动力响应，可为相关的科学研究和工程设计提供理论依据。表2、表3分别为不同于应力时，位移风振系数与实验结果比较情况。 　　 　　 　　由于《膜结构技术规程》与工程中现行作法对雪载荷分布形式的考虑存在偏差，田?B等人[9]首先给出根据膜结构期望使用年限重新确定基本雪压的方法，再采用最大熵法可靠度理论计算张拉膜结构在雪荷载满跨堆积和半跨堆积时各节点的最大可能变形(位移)，最后利用ANSYS求出相应工况的膜面应力分布，并与现行确定性方法所求得的膜面应力分布做对比分析。经过研究表明：可以用现行工程设计方法设计膜面承受满跨雪载的工况；但对于雪载半跨分布情况，工程中