膜結構景觀的計算方法
　　
　　一、膜結構景觀是大跨空間結構的主要形式之一，新穎獨特的建筑造型、優良的物理力學性能、簡潔結構體系，使其必將成為我國二十一世紀空間結構發展的主要方向。膜結構景觀作為一種柔性張力體系，與傳統的剛性結構在設計計算上有諸多不同之處，其設計計算的一般過程是：初始形態設計～靜力性能分析～風振響應分析～裁剪設計計算，其中每一環節都是不可或缺的。
　　
　　經驗告訴我們，一種新興結構的推廣和發展，與其計算理論方法、相應的設計軟件的成熟有相當大的關系。目前國內已經出現了諸如上海八萬人體育場等少數應用膜結構景觀的建筑，但遺憾的是其膜部分的設計、計算皆是國外公司所為，我國自行設計的大型膜結構景觀還很少。這當然是國外控制其關鍵技術，也是國內尚無成熟的計算理論方法及相應的計算設計程序的緣由。正因為如此，本文作者一方面進行膜結構景觀的計算理論方法的研究，一方面編制了相應的計算設計程序，旨在能為膜結構景觀在我國的應用發展盡微薄之力。 　　二、計算理論方法
　　
　?。ㄒ唬┓蔷€性有限元基本方程
　　
　　有限單元法是對工程結構進行數值分析方法，特別是在計算機應用越來越普及的今天。膜結構景觀實際設計中，索及桁架等加強、邊緣構件的應用是必不可少的，因此本文程序中包含了膜單元，索單元，桿單元，梁單元。其中空間膜單元定義為三結點的三角形等參元，考慮節點的ｘｙｚ三個方向的位移，但只計及面內的正應力σx、σy和剪應力Z xy。應用Ｕ•Ｌ法列式，可以得結構有限元基本迭代方程為：其中{R}為外荷載向量；{F}為ｔ時刻單元應力節點等效力向量；[kL]為線性應變增量剛度矩陣；[KNL]為非線性應變增量剛度矩陣，非線性方程組的求解，采用增量形式的 Full Newton-Raphson 法。具體推導見文獻[1]。
　　
　　（二）初始形態設計
　　
　　本文這里使用了“形態”這個詞匯，所謂的“形”就是幾何意義上的形狀，所謂的“態”就是結構的內力分布狀態。一種“形”對應一種“態”；反之亦然，一種“態”必然有一種“形”與它對應。膜結構景觀這類柔性結構有一個顯著的特點，即結構材料本身不具有剛度，由這些材料組成的構造體系還只是機構，只有當對其施加了預張力，它才具有了抵抗外荷載的結構剛度。這里應強調的是，不同的張力分布對應著不同的平衡形狀。這就是本文所說的“形”和“態”。當然，這是一個動態過程，對這個過程的研究就是—初始形態設計。這是膜結構景觀與傳統剛性結構計算的一個顯著區別，也膜結構景觀設計計算中的一個關鍵問題。
　　
　　膜結構景觀的設計中，在找到初始形態之前，并不能準確確定膜結構景觀的初始形狀和與之對應的預張力分布狀態，也就是說這時有兩個未知數：一個是初始形狀，一個是預張力分布狀態。這時我們會給定其中的一個而求解另一個，從而產生兩種思路。思路產生比較自然，是將初始形狀作為已知數，把初始預張力作為外荷載施加到結構上，求解達到平衡時的狀態。此法求解比較直接方便，普通的靜力計算程序就可以適用，并可以得到與設計者給定的曲面形狀相近的結果。但是它的缺點是終得到的預張力分布不再是設計給定的，當膜結構景觀曲面比較復雜不規則時，預張力的分布將會非常不均勻，造成施工安裝的困難和對受力性能的不利。
　　
　　思路亦是本文將采用的方法，是將初始預張力的分布狀態作為已知數，而把與之對應的平衡的形狀作為未知數來求解。此時的目的就是要得到具有給定初始預張力分布的一個平衡的初始形態，為了保證終得到的預張力分布即為初始假定的預張力，我們將舍棄變形協調條件和材料本構關系。也就是設計給定控制點的位置，以及預張力的分布狀態，然后尋求在此條件下的與之相對應的平衡的曲面形狀。此種方法的優點是終得到的初始形態的預張力分布即為初始假定的預張力，這會給膜結構景觀的受力性能和施工安裝帶來的好處。前面的分析可知本文的初始形態設計過程中，將忽略變形協調條件和材料本構關系，方程中與變形協調條件及材料本構關系有關的項均將被舍棄。因此，線性剛度在成形計算過程中將不起作用，結構的剛度僅由非線性剛度構成，如果忽略外荷載作用，則膜結構景觀的初始形態分析的有限元基本方程即為：                                                   （2）求解膜結構景觀初始形態的已知條件是，初始預張力（包括索的初始預拉力）的大小及分布、結構的邊界控制點位置，要求解的是結構上其余節點的坐標。具體的理論推導和計算步驟詳見文獻[1]。
　　
　?。ㄈ╈o力性能分析
　　
　　在推導得到含空間膜單元的非線性有限元基本方程后，膜結構景觀的靜力性能分析就變得簡單了，直接應用式（1）即可，但需要指出的是，由于膜材為非抗壓性材料，為此本文引入了由單元主應力判斷褶皺的方法。設單元主應力為σ1、σ2，且σ1>σ2，則：    σ1>0，σ2 >0，單元正常工作。σ1    <0，σ2 <0，單元退出工作。 σ1    >0，σ2 <0，單元為單向受拉，以轉換的應力和剛度代入方程重新迭代。
　　
　?。ㄋ模╋L振響應分析 　　很好的先將膜結構景觀離散化，根據風的概率統計特性，將風速模擬成時間的函數，即人工仿真生成風速時程。然后利用風的空間相關特性，將在結構各個結點處生成的不相關的風速轉換成空間相關的風速場。應用Morison公式將風速轉變為風壓，合理地確定風荷載，并作用于相應的單元結點上，利用有限元法在時間域內直接求解運動微分方程并求得結構的響應。在推導有限元基本方程的迭代公式時，在每一時間步長中考慮了結構的非線性特性、風與結構相互耦合作用的影響、風壓分布系數（體型系數）等等因素。由響應值中求得所要的統計信息，如結構振動的位移、速度、加速度的均值和均方差，以及相應的功率譜，并從中獲得膜結構景觀的風振響應特性。膜結構景觀質量很輕且剛度較弱，因此非線性風振響應分析將是其設計計算中的又一個關鍵問題、難點問題。本文應用隨機模擬法的思路，發展建立了一種可以適用于大型膜結構景觀非線性風振響應分析的方法。該方法的思路是，對工程結構設計計算來說，風作用的大小一般以風壓來表示，本文采用Morison公式來計算風荷載，它可更合理地反映物體的表面壓力。當不考慮風與結構耦合作用時的風荷載公式為：
　　
　　（3）之和是加速流中靜止的彈性結構所受到的總慣性力，總稱為慣性項。式中 是空氣質量密度，A為面積，H為膜厚度， 為風壓分布系數， 為風速，包含了平均風速和脈動風速?？梢钥闯?，Morison公式比一般公式增加了開頭的兩項，如果忽略前兩項，則即為我們通常所說的風荷載公式。如果考慮風與結構的耦合作用，則上式變為：
　　
　　（4） 分別為結構振動的位移、速度和加速度，而結構的運動方程為：
　　
　　（5）文中質量矩陣M采用集中質量矩陣，阻尼矩陣C采用Rayleigh阻尼，非線性的剛度矩陣K如前文所述，而荷載項P采用式（4）。顯然，這是一個復雜的非線性動力方程組，不僅剛度矩陣是隨位移變化而變化的，而且荷載項里還包含了結構的加速度和速度項。為此本文利用Newrk法和Newton-Raphson迭代法的思想，推導了膜結構景觀在風荷載作用下，考慮風與結構耦合作用的非線性動力增量平衡方程，終整理得：                           （6）其中： 方程中個參數的具體含義及取法可參見文獻[1]。
　　
　?。ㄎ澹┎眉粼O計計算
　　
　　膜結構景觀的曲面經初始形態設計得到，一般為復雜不規則、不可展的空間曲面形式，而且是由有限元的離散節點構成。因此，這里存在一個如何將平面的膜材料拼接成空間曲面的裁剪設計問題。裁剪設計是膜結構景觀設計中的一個關鍵問題。裁剪下料圖的準確與否直接關系到施工安裝后的平整度，也就是膜結構景觀的靜、動力計算的前提—初始形態設計得到的初始狀態與實際施工安裝后是否吻合，進而影響靜動力的計算結果。
　　
　　本文發展建立了一種應用廣義泛函變分取值，得到膜結構景觀曲面上的測地線，然后依據測地線計算裁剪線從而生成裁剪圖的方法。具體方法詳見文獻[1]，其過程一般如下：
　　
　　1、由初始形態設計程序求得膜結構景觀的空間曲面離散點坐標。
　　
　　2、根據建筑和結構上的要求確定測地線控制點的位置，計算出測地線軌跡坐標。這里建筑上的要求指的一是測地線的布置要美觀；二是相鄰兩條測地線間距離要控制在膜材幅寬的范圍內。而結構上的要求是要盡量避免結構在受荷時拉力大的方向與測地線方向垂直。
　　
　　3、依據測地線求出裁剪線。
　　
　　4、生成裁剪圖。當然，終的施工下料圖還要考慮繩邊索套，焊縫寬度及初始預張力、溫度應力、徐變等因素的影響。