膜结构荷载分析.pdfVIP

第 6卷 第 1期 空 间 结 构 [文章绾号]Ioo6—6578(2000)01—000905 膜结构的荷载分析 ④ 一 1) 叶 小兵 ， 袁 明武 沈3奸．。 (北京大学力学与工程科学系， 北京 100871) [摘 要] 本文采用非巍性 曲面膜单元和曲线童单元璃制 了膜蛄构 内力计算程序．对衷斥制 佧时糙型在属荷羲和 雪菏截佧用下应力、莅秽 的变化规律厦蛀折赴理 的方法律 7研究 t得 到的 鲒论 呻实 际工程 设计艰有 帮助 。 [关键词] 竺 生重主互蔼薄； [中圈分娄号] TU351 。 [殳献标识码] A 1 引 言 从 1971年起 ，很多学者将非线性有 限元技术[1][应用于膜结构 的实际设计 中，可 以说迄 今建成的大部分膜建筑都是 由根据这套理论编制的程序设计的。该理论采用三角形单元 ．考虑 几何非线性和皱折影响，把膜材看成各向异性的弹性材料 所得结果基本满足工程需要 。但 由 于所选单元精度不高，且材料矩阵难 以反映膜结构显著的徐变特征，所 以还不能说该理论就很 完备了 一直到现在很多学者还在从理论上继续深入研究膜这种特殊材料的力学性能0Ⅱ 。这 些研究一般不 以整个的膜建筑为研究对象 ，而是将一块膜张拉 固定成某一简单形状 ，然后考虑 各种复杂情况 ，如几何非线性与材料非线性同时存在 ；粘弹性时间历程对本构关系的修正；局 部受热与皱折松驰的扩展 ；膜索之间滑移现象等 应用薄膜力学理论和一些如有限元和有限差 分的数值方法．取得了很多杰出的研究成果 但在大型的工程设计程序 中还未见采用如此复杂 的计算理论 ．主要是因为大型工程本身设计计算量 已经很大，更复杂的理论将导致计算结果难 以控制。今后的 目标就是要研制精度更高又切实可行的计算程序 。 本文的内力分析使用的还是有限元方法，与以前不同的是采用 曲面 6节点三角形单元和 曲线索单元 ]，应变的线性部分引入了法向位移及单元的曲率和扭率，非线性部分仍然保 留法 向位移的影响项 。这样无论是每个单元还是各单元合并后的平衡方程都能很容易满足 ，迭代次 [收稿 日期] 1999—06—07 [作者简介] 叶小兵 (196争一)，男 ，湖北人．北 京大学力学与工程科学系博士后 ，主要从事空 问结构 的研 究 。 9 数大为减少，而变形结果也更符合真实情况，且 由于单元内各点应力都不相同，据此判断皱折 是否出现会更 为精确 。最后求出的每个单元的曲率和扭率对于判断初始找形的正误和优劣以 及裁剪下料都能提供很多非常有用的信息 采用 曲面单元，能够更好地反映膜结构的真实几何形状，因为 曲面单元的边界为 曲线．与 直线边界单元 比较，前者大大提高了拟台曲线边界的能力，减少了几何离散所带来的误差 ，计 算精度较高 ，且单元可取得相对大些．弥补了平面单元 由于节点较多 ，带宽太大的缺 陷。由于上 述优点 ，如果用 曲面单元划分的网格和用平面单元划分的网格相 同．则采用 曲面单元可 以比采 用平面单元得到更好 的计算结果 。 采用 曲面单元也带来一定的麻烦，如增加了位移表达式的复杂性和每个单元 内数值计算 的工作量 ，但对于现在高速度高容量的微机来说这并不是太大的问题 。 2 内力计算程序的介绍及模型计算结果分析 作用在膜结构上的荷载主要有风荷 、雪荷 ，而膜本身重量很轻 ，所 以地震力对膜 结构的影 响与风荷相 比可 以忽略不计 作用在结构上任一点处的风压力的计算公式为 W ()一 () o (1) 其 中， 。是基本风压力，在天津地 区约为 53．86kg／m ； ()为体型分布系数 由于柔性结构在风荷作用下变形较大，所 以膜结构表面风压力分布是不断变化的，要想准 确的计算出各点处的体型分布系数非常困难 。目前，国外基本上都是采用风洞实验测量在缩小 比例的刚性或柔性模