張佳毅，鄒貽權(quán)

(湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院，湖北 武漢 430068)

近年來隨著我國綜合國力的增強(qiáng)和計算機(jī)技術(shù)的進(jìn)步，結(jié)構(gòu)理論分析方法和施工技術(shù)等得到長足的進(jìn)步，由于我國體育事業(yè)的發(fā)展，各類輕質(zhì)、高強(qiáng)、大跨建筑隨之出現(xiàn)。因索結(jié)構(gòu)具有自重輕、節(jié)省材料、施工便捷能夠適應(yīng)多樣化的建筑造型等優(yōu)點引起了廣大設(shè)計人員的興趣，越來越多的應(yīng)用索結(jié)構(gòu)屋蓋的大型體育場館相繼出現(xiàn)[1]。但索結(jié)構(gòu)屬于柔性張拉結(jié)構(gòu)，必須施加一定的預(yù)應(yīng)力才能形成承受外部荷載所需的剛度和形狀[2]。因此確定索結(jié)構(gòu)零狀態(tài)、初始態(tài)、荷載態(tài)的幾何形狀是索結(jié)構(gòu)的結(jié)構(gòu)設(shè)計的過程中需要解決的關(guān)鍵問題。索結(jié)構(gòu)其類型主要有：單層索系、雙層索系、索桁體系等結(jié)構(gòu)體系。其中單層輪輻式索網(wǎng)結(jié)構(gòu)應(yīng)用得最晚，該結(jié)構(gòu)形式在賈比爾.艾哈邁德體育場才得到首次應(yīng)用。通過賈比爾.艾哈邁德體育場曾無法通過驗收的現(xiàn)象可以從側(cè)面體現(xiàn)單層輪輻式索網(wǎng)結(jié)構(gòu)的相關(guān)理論方法與施工技術(shù)仍有所欠缺[3]。在單層輪輻式索網(wǎng)結(jié)構(gòu)的相關(guān)研究中，其更多的是研究改進(jìn)力學(xué)找形方法并利用有限元分析軟件加以實現(xiàn)。

在現(xiàn)有的找形方法用于索桿體系時常使用平衡矩陣?yán)碚揫4]，力密度方法[5]和動力松弛方法[6]進(jìn)行計算時。這些方法的研究主要應(yīng)用于索膜結(jié)構(gòu)和索網(wǎng)結(jié)構(gòu)，針對單層輪輻式索網(wǎng)結(jié)構(gòu)的效率較低。

隨著計算機(jī)技術(shù)的發(fā)展，越來越多的設(shè)計人員開始使用參數(shù)化設(shè)計工具進(jìn)行設(shè)計。相比于傳統(tǒng)有限元分析軟件，RHINO和Grasshopper所提供的參數(shù)化平臺可以大大地提高設(shè)計建模和分析的效率，更易實現(xiàn)單層輪輻式索網(wǎng)結(jié)構(gòu)形狀確定的找形方法[7]。

本文擬對單層輪輻式索網(wǎng)結(jié)構(gòu)(圖1)的找形方法和原理進(jìn)行較為全面的闡述，并在參數(shù)化設(shè)計平臺RHINO和Grasshopper上結(jié)合算例證明方法的適用性。

圖 1 輪輻式單層索網(wǎng)結(jié)構(gòu)軸測圖

1 單層輪輻式索網(wǎng)結(jié)構(gòu)初始形態(tài)確定

對于索網(wǎng)結(jié)構(gòu)的初始形狀的確定問題，其實質(zhì)是在索網(wǎng)的幾何構(gòu)形中形成合適的預(yù)應(yīng)力分布使之成為一個自適應(yīng)平衡體系。單層輪輻式索網(wǎng)結(jié)構(gòu)與索膜結(jié)構(gòu)及正交索網(wǎng)結(jié)構(gòu)不同，通常索網(wǎng)結(jié)構(gòu)的初始形狀確定的問題是在給定的邊界條件下進(jìn)行，尋找符合初始態(tài)平衡條件的幾何形狀和與之相對應(yīng)的預(yù)應(yīng)力分布[8]。

然而，對于單層輪輻式索網(wǎng)結(jié)構(gòu)而言，其初始態(tài)幾何形狀具有一定的獨特性。在結(jié)構(gòu)的平面投影和外圈壓環(huán)的節(jié)點坐標(biāo)確定的條件下，則可以確定唯一的內(nèi)環(huán)索節(jié)點坐標(biāo)，從而確定唯一的單層輪輻式索網(wǎng)結(jié)構(gòu)的空間幾何形狀。

因此單層輪輻式索網(wǎng)結(jié)構(gòu)在已知支座節(jié)點X、Y、Z及內(nèi)環(huán)索節(jié)點X、Y坐標(biāo)的條件下，可合理利用該幾何特性，確定其初始狀態(tài)下的幾何位形。該幾何特性可被稱為“共點三索共平面”原理。

1.1 共點三索共平面確定原理

“共點三索共平面”原理主要利用的是拉索的節(jié)點平衡。對于兩根共節(jié)點的懸索，其索力大小均大于零時，索力可在一維空間平衡，形成一維的索結(jié)構(gòu)(圖2a)公式為：

F1+F2=0

(1)

對于三根共節(jié)點且不共線的懸索，其索力大小均大于零時，因“共節(jié)點的三根索共面”，索力可在二維平面上平衡，形成平面的索結(jié)構(gòu)(圖2b)，公式為：

(2)

對于共節(jié)點且不共面的懸索，其索力均大于零時，索力可在三維空間平衡，形成空間的索結(jié)構(gòu)(圖2c)，公式為：

圖 2 索結(jié)構(gòu)

(3)

1.2 單層輪輻式初始態(tài)幾何確定方法

針對本文的研究對象單層輪輻式索網(wǎng)結(jié)構(gòu)(圖3)，其內(nèi)環(huán)索節(jié)點均連接兩根環(huán)索與一根徑向拉索，由“三索共面原理”可知，三根通過環(huán)索節(jié)點連接的索必定共面。例如，拉索A′B′、BB′、B′C′通過節(jié)點B′連接，當(dāng)拉索BB′、A′B′、B′C′保持張力時，線段BB′、A′B′、B′C′必然在同一平面上，即點A、A′、B′、C′四點共面。

圖 3 輪輻式單層索網(wǎng)結(jié)構(gòu)節(jié)點示意圖

根據(jù)“共點三索共平面”原理，在確定平面投影尺寸的情況下通過外圈壓環(huán)的Z坐標(biāo)依次確定內(nèi)圈環(huán)索節(jié)點的Z坐標(biāo)，通過迭代循環(huán)的方式，最終確定索網(wǎng)的初始幾何位形。具體流程如下：

步驟1 確定初始己知量：

1)己知支座節(jié)點的X、Y、Z坐標(biāo)，即XA、YA、ZA、XB、YB、ZB……其中支座A為最高點，支座F為最低點，指定ZF為最低點，指定ZF=0；

2)給定內(nèi)圈環(huán)索節(jié)點的X、Y坐標(biāo)，即XA′、YA′、XB′、YB′……

步驟2 第一次幾何位形計算:

1)假定第一次迭代調(diào)整因子α1確定ZA′的數(shù)值，從而確定ZB′。

假定第一次迭代調(diào)整因子α1，其取值范圍為(0,1],α1通常取1。定義ZA′=α1ZA，可通過ZA得到ZA′。由于節(jié)點B′和節(jié)點G關(guān)于X軸′對稱(圖4a)，且節(jié)點A、A′、B′及G′共面，從而可以確定節(jié)點B′和節(jié)點G′的初始Z坐標(biāo)。

2)確定節(jié)點C′的初始Z坐標(biāo)ZC′。

通過步驟2-1可知節(jié)點A′、節(jié)點B及B′的Z坐標(biāo)(圖4b)，由于節(jié)點B、A′、B′及C′共面，從而可以確定節(jié)點C′的初始Z坐標(biāo)ZC′。

圖 4 步驟示意圖

3)重復(fù)上述操作，確定節(jié)點D′、節(jié)點E′及F′的Z坐標(biāo)，得到1/4結(jié)構(gòu)的所有節(jié)點坐標(biāo)。

4)驗算節(jié)點F坐標(biāo)誤差

由于1/4結(jié)構(gòu)的所有節(jié)點坐標(biāo)均可通過“共點三索共平面”原理確定，節(jié)點A到節(jié)點E分別在平面B′A′G′,C′B′A′,D′C′B′,E′D′C′及F′E′D′內(nèi)(圖4c)。在此條件下，需驗證3)中節(jié)點F′坐標(biāo)誤差，以確定節(jié)點F′是否在平面E′F′M′內(nèi)。因節(jié)點E′、M′關(guān)于Y軸對稱，可知兩點Z坐標(biāo)相同，從而得到平面E′F′M′。

設(shè)L1為該次支座節(jié)點F至E′F′M′的距離，若點F位于平面E′F′M′上方，L取正值，進(jìn)入步驟3；若點F位于平面E′F′M′下方，L取負(fù)值，增大1)中α1取值，重復(fù)1)，直至L為正。

步驟3 第二次幾何位形計算:

假定第二次迭代調(diào)整因子α2，其取值范圍為(0,1],α2通常取0.05。定義ZA′=α2ZA，其余運(yùn)算步驟與步驟2中相同。

設(shè)L2為該次支座節(jié)點F至E′F′M′的距離，若點F位于平面E′F′M′上方，L2取正值，減小α2的取值，重復(fù)步驟3，直至L2為負(fù)；若點F位于平面E′F′M′下方，進(jìn)入步驟4。

步驟4 迭代確定點A′的Z坐標(biāo):

定義第(i+2)次的調(diào)整因子為α3=(α1+α2)/2.其余運(yùn)算步驟與步驟二中相同。設(shè)L3為該次支座節(jié)點F至E′F′M′的距離。

若L3小于誤差限制值[Um]，證明節(jié)點位于平面E′F′M′內(nèi)，滿足“共點三索共平面”原理，終止運(yùn)算；

若L3大于誤差限制值[Um]且L3大于零，則調(diào)整α1=α3，重復(fù)步驟4，直至L3小于限制值[Um]；

若L3大于誤差限制值[Um]且L3小于零，則調(diào)整α2=α3，重復(fù)步驟4，直至L3小于限制值[Um]；

步驟5 生成完整結(jié)構(gòu):

由上述操作可得1/4結(jié)構(gòu)的所有節(jié)點坐標(biāo)，根據(jù)邊界處的對稱條件生成完整的結(jié)構(gòu)模型。

1.3 四支座簡化算例

本研究通過Grasshopper可視化編程實現(xiàn)迭代流程，針對簡化的四支座簡化算例進(jìn)行迭代結(jié)果的實現(xiàn)。并與理論解析的四支座簡化算例的理論結(jié)果進(jìn)行對比，證明該“單層輪輻式索網(wǎng)結(jié)構(gòu)初始形態(tài)方法”的可行性。

本算例為四支座單層輪輻式索網(wǎng)結(jié)構(gòu)，算例計算條件(圖5)。

圖 5 四支座簡化算例示意圖

1.3.1算例理論解析該算例結(jié)構(gòu)關(guān)于XY軸對稱，可知ZA′=ZC′，ZB′=ZD′。

(4)

將節(jié)點A、A′坐標(biāo)代入式(4)中得直線AA′表達(dá)方程：

(5)

因直線B′D′與AA′共面，且直線B′D′與Y軸平行，將節(jié)點B′坐標(biāo)代入式(5)中，得B′點坐標(biāo)表達(dá)式：

ZB′=2ZA′-2500

(6)

將節(jié)點B、B′坐標(biāo)代入式(4)中得直線BB′表達(dá)方程：

(7)

因直線B′B′與A′C′共面，且直線A′C′與X軸平行，將節(jié)點B′坐標(biāo)代入式(5)中，得A′點坐標(biāo)表達(dá)式，同時聯(lián)立式(6)：

ZA′=1666.67

(8)

最終可得：

(9)

1.3.2初始態(tài)幾何確定方法實現(xiàn)算例本研究在grasshopper平臺中利用可視化編程的方式實現(xiàn)“單層輪輻式索網(wǎng)初始態(tài)幾何確定方法”對四支座算例進(jìn)行找形，在本算例中，取第一次迭代調(diào)整因子α1=1，第二次迭代調(diào)整因子α2=0.05，誤差限制值[Um]=0.05,迭代程序如圖6所示。

圖 6 迭代程序示意圖

經(jīng)過12次迭代，最終得到結(jié)果(圖7)：

圖 7 找形結(jié)果示意圖

(10)

“單層輪輻式索網(wǎng)初始態(tài)幾何確定方法”實現(xiàn)的四支座算例結(jié)果與理論解析計算結(jié)果可知，ZA′的誤差值為0.0079%，ZB′的誤差值為0.0045%，該誤差值結(jié)果證明“單層輪輻式索網(wǎng)初始態(tài)幾何確定方法”的可靠性，單層輪輻式索網(wǎng)結(jié)構(gòu)初始態(tài)幾何可用該方法找形。

2 單層輪輻式索網(wǎng)結(jié)構(gòu)零狀態(tài)形狀確定

單層輪輻式索網(wǎng)結(jié)構(gòu)初始態(tài)幾何位形可由上節(jié)所述“單層輪輻式索網(wǎng)初始態(tài)幾何確定方法”進(jìn)行確定，因此單層輪輻式索網(wǎng)結(jié)構(gòu)零狀態(tài)形狀確定的問題屬于初始態(tài)幾何已知問題。針對該類問題目前較為適用的是逆迭代的方法。

2.1 逆迭代基本原理

逆迭代法最開始是用在張弦梁結(jié)構(gòu)的幾何形狀確定的問題上[9]，然后經(jīng)過發(fā)展推廣應(yīng)用在了索網(wǎng)結(jié)構(gòu)找形分析上。逆迭代是基于相似性原理的分析方法，其基本原理為：在給定預(yù)應(yīng)力的條件下，在一次迭代后將所得的節(jié)點位移反向作用在初始位形上，更新節(jié)點坐標(biāo)，進(jìn)行下一次迭代。隨著迭代次數(shù)的增加，在預(yù)應(yīng)力和自重的作用下結(jié)構(gòu)變形后的位形，該迭代時的幾何位形，可被認(rèn)定為零狀態(tài)幾何[10]。其迭代流程如下：

1)給定索結(jié)構(gòu)的初始態(tài)幾何X,定義結(jié)構(gòu)的零狀態(tài)幾何X0；

2)給定索的預(yù)張力P，初始態(tài)幾何X在索自重和預(yù)張力(通常不平衡)作用下結(jié)構(gòu)產(chǎn)生初始位移d1，結(jié)構(gòu)變形成為X1；

3)將位移d1反向作用于X得到初始迭代幾何XK=X-d1；

4)初始迭代幾何Xk在索自重和預(yù)張力作用下結(jié)構(gòu)產(chǎn)生位移dk，結(jié)構(gòu)變形成為Xk+1；

5)再次將變形dk反向作用于Xk得到迭代幾何Xk+2=Xk+1-dk；

6)判斷迭代幾何Xk+2初始幾何X的位移是否小于誤差極限值[Dm]，若小于誤差極限值[Dm]，迭代幾何Xk+2即為所求零狀態(tài)幾何X0；若大于誤差極限值[Dm]，則令k=k+2,重新進(jìn)入步驟4)。

2.2 零狀態(tài)四支座簡化算例

本研究針對2.3節(jié)中簡化的四支座簡化算例(圖8)進(jìn)行相應(yīng)的零狀態(tài)找形的研究，通過grasshopper可視化編程實現(xiàn)2.1節(jié)所述迭代流程。該算例的彈性模量均為1.6×105N/mm2，拉索強(qiáng)度等級為1670 MPa，徑向索BB′設(shè)置主動索力F=300 kN。

圖 8 零狀態(tài)四支座簡化算例示意圖

2.2.1算例理論解析由“共節(jié)點三索共平面”原理可知，徑向索與兩根環(huán)索在相連節(jié)點處平衡，即：

(11)

由結(jié)構(gòu)關(guān)于XY對稱可知，F(xiàn)A′B′=FA′D′=FB′C′=FC′D′,代入得：

(12)

由FBB′=300 kN,代入式(12)得：

(13)

2.2.2逆迭代法算例實現(xiàn)根據(jù)2.1節(jié)中所述逆迭代流程，基于grasshopper參數(shù)化設(shè)計平臺及基于該平臺開發(fā)的有限元分析插件karamba編制零狀態(tài)找形程序。因通過誤差極限值[Dm]判斷逆迭代找形是否成功，在已知預(yù)應(yīng)力分布的條件下，應(yīng)根據(jù)逆迭代求解零狀態(tài)幾何的預(yù)應(yīng)力與已知的預(yù)應(yīng)力進(jìn)行對比，判斷逆迭代法確定零狀態(tài)是否可行。在該算例中設(shè)置最大迭代次數(shù)為20次，誤差極限值[Dm]=1.0 mm。

表1 索力對比表

該算例迭代20次后，迭代結(jié)束，此時迭代幾何X22與初始幾何X對比。迭代結(jié)果通過對比可知，該算例理論索力與逆迭代索力的索力誤差值最大為2.04%，誤差在可接受的范圍內(nèi)。并且隨著誤差極限值[Dm]的減小和迭代次數(shù)的增加，計算結(jié)果將更為精確。

3 結(jié)束語

針對單層輪輻式索網(wǎng)結(jié)構(gòu)初始狀態(tài)和零狀態(tài)找形方法進(jìn)行了系統(tǒng)的研究，并通過參數(shù)化平臺RHINO和 Grasshopper結(jié)合簡化算例對找形方法進(jìn)行驗證。研究得出以下結(jié)論：

1)針對單層輪輻式索網(wǎng)結(jié)構(gòu)的初始狀態(tài)幾何應(yīng)根據(jù)單層輪輻式索網(wǎng)結(jié)構(gòu)的幾何特性進(jìn)行迭代處理進(jìn)行幾何位形確定；

2)逆迭代法可用于索桿體系的零狀態(tài)找形問題，同樣適用于單層輪輻式索網(wǎng)結(jié)構(gòu)的零狀態(tài)找形問題；

3)參數(shù)化平臺RHINO和Grasshopper結(jié)合基于該平臺的各類插件可實現(xiàn)單層輪輻式索網(wǎng)結(jié)構(gòu)的找形方法。