索托節(jié)點(diǎn)的數(shù)值模擬方法研究

楊維國, 王萌, 王亞, 阮志勇

(北京交通大學(xué) 土木建筑工程學(xué)院，北京 100044)

?

索托節(jié)點(diǎn)的數(shù)值模擬方法研究

楊維國, 王萌, 王亞, 阮志勇

(北京交通大學(xué) 土木建筑工程學(xué)院，北京 100044)

為了更加準(zhǔn)確地模擬索托結(jié)構(gòu)中的索托節(jié)點(diǎn)，本文利用ANSYS軟件進(jìn)行整體桁架的模擬并進(jìn)行相關(guān)的有限元分析，提出一種新的索托節(jié)點(diǎn)數(shù)值模擬方法，即不變剛度短梁法，并與其他3種數(shù)值模擬方法(自由度耦合方法、自由度耦合加變剛度彈簧單元方法、點(diǎn)-點(diǎn)接觸方法)進(jìn)行分析對比。結(jié)果表明：4種方法對桁架結(jié)構(gòu)撓度和構(gòu)件內(nèi)力引起的差別并不明顯，結(jié)果基本吻合，內(nèi)力走勢一致。由于自由度耦合加變剛度彈簧單元方法已得到驗(yàn)證，故說明本文提出的不變剛度短梁法是合理可行的；自由度耦合加變剛度彈簧單元方法和不變剛度短梁法能更好地模擬索托節(jié)點(diǎn)處的摩擦力行為，而且這兩種方法理論明確且計(jì)算簡單方便，是相對較好的索托節(jié)點(diǎn)數(shù)值模擬方法。

索托節(jié)點(diǎn)；數(shù)值模擬；不變剛度短梁法；自由度耦合；變剛度彈簧；點(diǎn)-點(diǎn)接觸；摩擦力

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20160623.0827.006.html

索托結(jié)構(gòu)是一種新型大跨空間結(jié)構(gòu)形式[1-2]，是在總結(jié)斜拉結(jié)構(gòu)和懸索結(jié)構(gòu)優(yōu)缺點(diǎn)的基礎(chǔ)上提出的全新結(jié)構(gòu)形式。索托結(jié)構(gòu)是指索體從結(jié)構(gòu)的底部穿過，將結(jié)構(gòu)托起，索托結(jié)構(gòu)外形類似于斜拉結(jié)構(gòu)，但受力模式卻完全不同。索托結(jié)構(gòu)中索為結(jié)構(gòu)提供豎向彈性支撐，同時(shí)也有效地減小了索對主體結(jié)構(gòu)的有害水平分力，使結(jié)構(gòu)能夠?qū)崿F(xiàn)更大的跨度。

索托結(jié)構(gòu)設(shè)計(jì)中如何對索托節(jié)點(diǎn)進(jìn)行準(zhǔn)確模擬是一個(gè)關(guān)鍵問題。對索托結(jié)構(gòu)進(jìn)行理論計(jì)算需確定索托節(jié)點(diǎn)接觸問題的種類，運(yùn)用合適的求解方法來滿足工程精度及收斂性要求?？紤]到索托節(jié)點(diǎn)計(jì)算的非線性及對計(jì)算工作量的要求，需要在索托結(jié)構(gòu)體系的計(jì)算中對索托節(jié)點(diǎn)采用簡化處理的方法。既不能違背索托結(jié)構(gòu)的基本性能，又要達(dá)到簡化計(jì)算的目的[3]。索托結(jié)構(gòu)中，拉索可以在索托節(jié)點(diǎn)處自由地滑動使兩側(cè)索力趨于平衡，但摩擦力是不可避免的因素，所以索托節(jié)點(diǎn)兩側(cè)索力不相等。根據(jù)索托預(yù)應(yīng)力節(jié)點(diǎn)足尺試驗(yàn)可知[3]，拉索摩擦損失隨張拉力的增大而增大，損失值約為張拉力的1%～5%，如何模擬索托節(jié)點(diǎn)處的摩擦力是準(zhǔn)確模擬索托節(jié)點(diǎn)受力行為最主要的問題。

在現(xiàn)有的3種模擬方法(自由度耦合方法[4]、自由度耦合加變剛度彈簧單元方法[5]、點(diǎn)-點(diǎn)接觸方法[3])基礎(chǔ)上，提出了不變剛度短梁法。以桁架結(jié)構(gòu)為例，對4種方法進(jìn)行桁架結(jié)構(gòu)撓度、構(gòu)件內(nèi)力與索托節(jié)點(diǎn)摩擦力行為的對比和評價(jià)。最后確定模擬索托節(jié)點(diǎn)最準(zhǔn)確的方法。

1　模型參數(shù)與荷載條件

1.1模型參數(shù)

本文以北京蘆城體校曲棍球館索托結(jié)構(gòu)為例進(jìn)行計(jì)算[6-7]，取三榀索托桁架為研究對象，結(jié)構(gòu)整體示意如圖1所示。該索托桁架為體外施力式索托桁架結(jié)構(gòu)[8]，索體從主體鋼結(jié)構(gòu)的底部穿過，拉索兩端錨固在邊柱上，斜向索與水平向的夾角(所謂索夾角)θ為20°。為了保證邊柱不會因?yàn)樗鞯膹埨l(fā)生側(cè)向位移，柱主要以軸力為主，在每根邊柱從柱頂往外側(cè)拉一根索，拉索與豎直方向的夾角為17°，從地面錨固到邊柱為11 m，見圖1。

桁架平面布置和立面布置如圖 2所示，上弦桿

外側(cè)一列從左到右分別為S1到S31，內(nèi)側(cè)一列從左到右分別為S32到S62，下弦桿從左到右分別為X1到X30，桁架上部豎桿從左到右分別為P1到P32，桁架上部斜桿從左到右分別為Q1到Q32。腹桿編號從左到右分別為F1到F124。桁架各桿件及索的具體材料屬性如表1 所示。

圖1　索托桁架結(jié)構(gòu)計(jì)算模型示意圖(單位：m)Fig.1　Calculation model of cable supported truss structure(unit:m)

表1　結(jié)構(gòu)構(gòu)件材料屬性

1.2荷載條件

1.2.1桁架荷載

結(jié)構(gòu)自重由軟件分析程序自動考慮，恒載Gk取1.5 kN/m2(不含自重)，按 《建筑結(jié)構(gòu)荷載規(guī)范》 取不上人屋面活荷載Qk取0.5 kN/m2，計(jì)算桿件內(nèi)力和結(jié)構(gòu)撓度分別采用以下荷載工況：

工況1：1.35Gk+0.98Qk=2.515 kN/m2

工況2：1.0Gk+1.0Qk=2.0 kN/m2

計(jì)算桿件內(nèi)力時(shí)，荷載值取2.515 kN/m2，索托桁架間距為6 m，桁架節(jié)間間距為3 m，換算成節(jié)點(diǎn)荷載等于22.635 kN。由于計(jì)算模型是三榀桁架，桁架之間用鋼管鏈接，桁架結(jié)構(gòu)水平剛度得到提高，這與真實(shí)地模擬實(shí)際工程中結(jié)構(gòu)的水平剛度足夠大相符，不會發(fā)生側(cè)向失穩(wěn)，故只考慮豎向力作用。

1.2.2拉索張力

為使索與主體鋼結(jié)構(gòu)協(xié)調(diào)工作，需要對索施加適當(dāng)?shù)膹埩Α＝Y(jié)構(gòu)在施工和工作階段是兩種不同的狀態(tài)：1)施工狀態(tài)，即初始平衡狀態(tài)，此時(shí)結(jié)構(gòu)受自重、部分施工荷載和預(yù)應(yīng)力的作用；2)工作狀態(tài)，結(jié)構(gòu)在自重、外部荷載和預(yù)應(yīng)力作用下的平衡狀態(tài)。考慮到結(jié)構(gòu)在工作狀態(tài)下的撓度和應(yīng)力比初始狀態(tài)下大，對拉索施加張力時(shí)只考慮結(jié)構(gòu)在工作狀態(tài)下的受力情況，即結(jié)構(gòu)在自重、外部荷載和預(yù)應(yīng)力作用下的受力狀態(tài)，結(jié)構(gòu)內(nèi)力和撓度分別用荷載工況1和工況2計(jì)算。

為了找到合理的張拉值，本文通過有限元分析，在純桁架、拉索零張力與某數(shù)值張力3種情況下對結(jié)構(gòu)的撓度、桿件與索的應(yīng)力比進(jìn)行控制對比。計(jì)算結(jié)果如表2所示，經(jīng)過對比可得，桁架結(jié)構(gòu)在拉索350 kN張拉下，主體結(jié)構(gòu)的撓度以及內(nèi)力得到很好的改善，索與主體鋼結(jié)構(gòu)能良好地協(xié)調(diào)工作。說明桁架各桿件截面，拉索截面以及施加的張力取值得當(dāng)，所以本文通過初始應(yīng)變法來對索施加350 kN的張拉力。

表2　不同情況下結(jié)構(gòu)的變形以及桿件應(yīng)力

2　四種索托節(jié)點(diǎn)數(shù)值模擬方法

本文用ANSYS建立三榀桁架，如圖3所示。為了更真實(shí)地模擬實(shí)際工程中索托桁架結(jié)構(gòu)的受力狀態(tài)，索托結(jié)構(gòu)中索的力學(xué)分析模型要滿足大位移、小應(yīng)變假定，計(jì)算時(shí)考慮大變形和應(yīng)力剛化的影響[9-11]在ANSYS索托桁架模型中，根據(jù)構(gòu)件受力不同，采用不同的單元模擬。桁架上弦桿和下弦桿在荷載作用下，除了產(chǎn)生軸力以外，還會產(chǎn)生彎矩和剪力，故采用Beam188單元；邊柱同樣采用Beam188單元；桁架腹桿主要受軸力，彎矩和剪力很小，故采用Link8單元；拉索由于只能承受拉力，故采用Link 10單元。

圖3　有限元模型立體圖Fig.3　The space diagram of finite element model

2.1自由度耦合法[5]

自由度耦合法是將連續(xù)索與其他構(gòu)件的連接處分開設(shè)置節(jié)點(diǎn)，并將此兩節(jié)點(diǎn)沿兩側(cè)拉索夾角平分線方向上的自由度耦合在一起，通過ANSYS有限元方法求解。這種方法能解決連續(xù)索滑移問題，操作簡單但模擬不了索托節(jié)點(diǎn)處拉索與主體鋼結(jié)構(gòu)的摩擦接觸。

2.2自由度耦合加彈簧的方法[6]

利用自由度耦合和變剛度彈簧單元來處理滑移索摩擦問題，即將拉索與桁架結(jié)構(gòu)接觸的節(jié)點(diǎn)與桁架上對應(yīng)節(jié)點(diǎn)在切線方向建立變剛度彈簧單元，徑向進(jìn)行自由度耦合。在索滑動時(shí)，彈簧的內(nèi)力即為滑道兩端索力之差(即該處的摩擦力)在彈簧方向上的分量。若無摩擦滑動，則彈簧剛度為零。若有摩擦滑移，則切線方向相當(dāng)于布置了一有剛度的彈簧，彈簧剛度由試驗(yàn)得到的摩擦力與滑移距離關(guān)系確定，索托節(jié)點(diǎn)模型如圖4所示。

圖4　索托節(jié)點(diǎn)構(gòu)造有限元模型示意圖Fig.4　The finite element model of cable supported joint

文獻(xiàn)[12]對預(yù)應(yīng)力索-鋼轉(zhuǎn)換節(jié)點(diǎn)進(jìn)行力學(xué)性能分析、足尺模型試驗(yàn)研究和ANSYS模擬分析，總結(jié)了預(yù)應(yīng)力損失、摩擦系數(shù)的變化規(guī)律。預(yù)應(yīng)力損失的試驗(yàn)測試結(jié)果可供此類預(yù)應(yīng)力鋼結(jié)構(gòu)工程參考。

應(yīng)用試驗(yàn)結(jié)果，采用具有非線性功能的單向彈簧單元Combin39來模擬變剛度單元，并對其輸入力-變形曲線，實(shí)現(xiàn)利用自由度耦合和變剛度彈簧單元來處理滑移索摩擦問題。

2.3點(diǎn)-點(diǎn)接觸法[3]

點(diǎn)-點(diǎn)接觸方法是將拉索與桁架接觸的節(jié)點(diǎn)、桁架上對應(yīng)的節(jié)點(diǎn)分開建立，在兩段拉索夾角平分線的方向上兩個(gè)節(jié)點(diǎn)間有一個(gè)較小的距離，然后用Contact52單元將兩點(diǎn)連接起來。Contact52單元的一些主要輸入?yún)?shù)有：法向剛度KN、粘性剛度KS、初始條件START、摩擦系數(shù)MU。

2.4不變剛度短梁法

本文提出一種新的索托節(jié)點(diǎn)數(shù)值模擬方法，即不變剛度短梁法。該方法的基本原理是：短梁與主梁鋼結(jié)構(gòu)的鏈接是剛接，短梁本身有一定的抗彎剛度但很小，在對拉索施加預(yù)拉力的時(shí)候，短梁呈彎曲變形，由于短梁本身具有抗彎能力所以使得兩側(cè)索力不相等。短梁的抗彎能力用來模擬索托節(jié)點(diǎn)處的摩擦力，采用Beam188/189單元來模擬短梁。

2.4.1短梁法理論

如圖5(a)所示，在初始狀態(tài)時(shí)，兩側(cè)索的初始張拉均為T0，短梁AO還沒變形，結(jié)構(gòu)處于平衡狀態(tài)。圖5(b)中，右側(cè)索施加張力時(shí)，由于短梁已有彎曲變形，短梁本身的抗彎能力使兩側(cè)索力不相等。右側(cè)索的索力為T2=T0+F，左側(cè)索的索力為T1，T0

(1)

短梁頂端O的水平位移為左邊索在(T1-T0)作用下的伸長量：

(2)

采用結(jié)構(gòu)力學(xué)的位移法，求得短梁頂端O點(diǎn)所受到的剪力為

(3)

式(3)體現(xiàn)了短梁的抗彎剛度與拉索幾何參數(shù)及抗拉剛度的關(guān)系。假如短梁截面是矩形，而且b=h(寬=高)，那么短梁截面尺寸為

(4)

從式(4)可知，該方法的實(shí)質(zhì)是根據(jù)摩擦損失值FS和拉索在索托節(jié)點(diǎn)處的滑移量Δ來確定短梁的剛度(短梁的長度和截面尺寸)。該方法不用建立復(fù)雜接觸，較為簡便靈活，故本文提出此方法并進(jìn)行研究，驗(yàn)證其可行性。

圖5　拉索受力簡圖Fig.5　Diagram of cable force analysis

2.4.2短梁法應(yīng)用

把上述短梁理論應(yīng)用在本文索托桁架結(jié)構(gòu)中的索托節(jié)點(diǎn)，如圖6所示，索夾角為θ，給拉索的初始張拉力為T0，在結(jié)構(gòu)自重以及外荷載作用下，結(jié)構(gòu)有豎向位移，拉索長度變長，使拉索的預(yù)應(yīng)力變大，由于索托節(jié)點(diǎn)處摩擦力的存在，導(dǎo)致兩側(cè)索力不相等，即T2>T1>T0。因?yàn)榻Y(jié)構(gòu)具有兩個(gè)索托節(jié)點(diǎn)，所以索托節(jié)點(diǎn)處拉索的滑移量就是水平索在(T1-T0)作用下的伸長量的一半，即

(5)

圖6　有限元模型中拉索的受力狀態(tài)Fig.6　Cable load condition in a static finite element model

如圖7所示，短梁是在兩側(cè)索夾角平分線方向上(所謂徑向)建立的。摩擦力使兩側(cè)索力不相等，即T2>T1，預(yù)應(yīng)力摩擦損失在切線方向上的分量，即短梁上A點(diǎn)所受到的剪力為

(6)

為了使結(jié)構(gòu)受力更合理，索的效率更高，一般索夾角θ都小于60°，水平索的伸長量近似等于短梁端部A點(diǎn)在其切線方向上的位移，也是索托節(jié)點(diǎn)處拉索的滑移量，即

(7)

圖7　索托節(jié)點(diǎn)處短梁法示意圖Fig.7　The short beam method in cable-supported joint

確定了摩擦損失值FS和拉索在索托節(jié)點(diǎn)處的滑移量Δ就可以通過式(4)來確定短梁的長度和截面尺寸(假如短梁截面是矩形，b=h)。短梁的長度和截面尺寸要滿足一定的條件：

1)短梁的抗彎剛度足夠小，為了模擬兩側(cè)索力的摩擦損失，一般工程中摩擦損失力約為張拉力的1%～5%；在拉索摩擦損失力和拉索滑移量已知的情況下，可以通過式(8)得到幾組短梁的幾何參數(shù)(l，b，h)。

2)短梁的抗壓剛度應(yīng)足夠大，能為主體鋼結(jié)構(gòu)提供足夠大的支撐力，但短梁本身不被壓縮變形；另外，短梁長度不應(yīng)太長，否則會影響到斜索與水平向的夾角。

本文中的索托桁架結(jié)構(gòu)，拉索的初始張拉力T0=350 kN、θ=20°，在結(jié)構(gòu)自重以及外荷載作用下，拉索預(yù)應(yīng)力在無摩擦情況下達(dá)到約775 kN；假如要模擬摩擦力損失為3%，即T2-T1=775×3%=23.25 kN，通過式(6)算出FS≈23 kN。水平索(即中間的拉索)力增加量為T1-T0=775 kN×97%-350 kN≈400 kN，索托節(jié)點(diǎn)處拉索的滑移量約為Δ=0.020 m，通過式(8)得到不同短梁長度的不同截面尺，如表3所示，ANSYS分析結(jié)果如表4所示。

表3　不同短梁幾何參數(shù)

表4　不同幾何參數(shù)短梁的計(jì)算分析結(jié)果

從表4可知，不同的短梁幾何參數(shù)組，得到的結(jié)構(gòu)撓度、索力以及拉索的滑移量也不同，但差別不大，短梁的抗壓剛度越大，結(jié)果差別越小。因?yàn)槎塘汗潭ㄔ谥髁轰摻Y(jié)構(gòu)的的節(jié)點(diǎn)上，而非剛度較大的物體上，所以結(jié)構(gòu)在自重以及外荷載作用下，桁架本身也有變形，使桁架上對應(yīng)的索托節(jié)點(diǎn)有了位移變化，導(dǎo)致兩側(cè)索力有所偏差。因此，只要短梁的抗壓剛度足夠大，結(jié)構(gòu)撓度、拉索滑移量以及拉索支撐力都差別不大，只需調(diào)整索力摩擦損失值到合理值即可。本文中，先選出一組抗壓剛度足夠大的短梁幾何參數(shù)，用ANSYS計(jì)算并分析與需要模擬的索力摩擦損失是否接近，用此方法循環(huán)直到結(jié)果接近，選出較合適的尺寸。

經(jīng)計(jì)算建議短梁在壓力(拉索給主體結(jié)構(gòu)的支撐力)作用下，軸向壓縮量不宜大于1 mm。本文中，拉索提供的支撐力約(T1+T2)sinθ/2≈265kN，短梁的抗壓剛度應(yīng)大于265 MPa，經(jīng)過計(jì)算和挑選，最終選擇短梁截面尺寸為b=h=0.010 4m，長度l=0.080 m。

3　4種模擬方法的結(jié)論對比

本文通過ANSYS建立了4種索托節(jié)點(diǎn)模擬方法(自由度耦合加變剛度彈簧方法已得到驗(yàn)證)，得到索托桁架結(jié)構(gòu)的桿件內(nèi)力，結(jié)構(gòu)豎向位移，索力摩擦損失以及拉索的滑移量，并對計(jì)算結(jié)果進(jìn)行了分析和對比，最后對每種方法進(jìn)行評估，旨在找出較簡單且能滿足工程精度要求的索托節(jié)點(diǎn)建模方法。

由于索托桁架結(jié)構(gòu)左右對稱，作用荷載也是對稱的，左右對稱桿件受力基本相同，故本文只列出左半部桁架桿件的內(nèi)力及左邊斜向索和水平索的索力。

3.1結(jié)構(gòu)整體變形及索力的對比

經(jīng)ANSYS計(jì)算可得利用4種不同索托節(jié)點(diǎn)模擬方法計(jì)算得到的桁架結(jié)構(gòu)整體最大豎向位移，同時(shí)計(jì)算出索托節(jié)點(diǎn)處拉索的滑移量及索力，如表5所示。

表5　桁架結(jié)構(gòu)整體和索托節(jié)點(diǎn)的分析結(jié)果

根據(jù)表5的數(shù)據(jù)可得如下結(jié)論：

1) 不同索托節(jié)點(diǎn)建模方法計(jì)算得到的結(jié)構(gòu)整體豎向位移和變形差別不大，4種方法均能為桁架提供的支撐力。

2) 自由度耦合方法不能模擬索托節(jié)點(diǎn)處的非線性摩擦，所以兩側(cè)索相等，由于拉索可以無摩擦滑移，索托節(jié)點(diǎn)處的滑移量最大。

3) 自由度耦合加彈簧法、點(diǎn)-點(diǎn)接觸法及短梁法都可以模擬索托節(jié)點(diǎn)處的索力摩擦損失，計(jì)算兩側(cè)索力，索力摩擦損失以及滑移量差別不大，均可進(jìn)行參數(shù)控制，更接近工程中索托節(jié)點(diǎn)的摩擦情況。

3.2索托桁架結(jié)構(gòu)構(gòu)件內(nèi)力的對比

經(jīng)ANSYS計(jì)算可得4種不同索托節(jié)點(diǎn)模擬方法下索托桁架結(jié)構(gòu)構(gòu)件內(nèi)力，桿件類別分別為：上弦桿、下弦桿、腹桿、平直桿、平斜桿，對各種桿件的內(nèi)力值進(jìn)行分析和對比。

3.2.1上弦桿內(nèi)力

根據(jù)圖8對比分析可知：4種方法計(jì)算得到的上弦桿內(nèi)力變化趨勢基本一致，結(jié)果基本吻合，并且4種方法所得上弦桿內(nèi)力差別不明顯，不同方法影響較小。

3.2.2腹桿內(nèi)力

由于腹桿受壓桿和受拉桿相間存在，本文中腹桿內(nèi)力分兩組列出：A組編號是1、2、5、6、…、57、58、61、62；B組編號是：3、4、7、8、…、55、56、59、60 。計(jì)算結(jié)果如圖9所示。

由圖9可得：腹桿的內(nèi)力走勢與文獻(xiàn)[9]里的結(jié)論相符，桁架腹桿受拉與受壓桿件相間分布，腹桿在索托節(jié)點(diǎn)處內(nèi)力降低。4種方法計(jì)算出來的腹桿內(nèi)力值非常相近，說明4種方法對腹桿內(nèi)力影響較小。

圖8　4種索托節(jié)點(diǎn)模擬方法計(jì)算的上弦桿的內(nèi)力值Fig.8　Internal forces of upper chord under calculated by four methods

圖9　4種索托節(jié)點(diǎn)模擬方法計(jì)算的腹桿內(nèi)力值Fig.9　Internal forces of web member calculated by four methods

3.2.3下弦桿

由圖10可得：與上弦桿和腹桿內(nèi)力相比，4種方法對下弦桿影響相對明顯，但內(nèi)力結(jié)果基本吻合，變化趨勢也一樣。因此4種方法對下弦桿內(nèi)力影響較小。

圖10　四種索托節(jié)點(diǎn)模擬方法計(jì)算的下弦桿內(nèi)力值Fig.10　Internal forces of lower chord calculated by four methods

3.2.4桁架上部的平直桿和平斜桿內(nèi)力

由圖11和12可得桁架上部的平直桿都是受拉桿，平斜桿都是受壓桿，這兩種桿件受力較小。這兩種桿件從桁架端部到桁架跨中，桿件內(nèi)力慢慢變大，4種方法計(jì)算出來的內(nèi)力變化趨勢一樣，結(jié)果基本吻合。

3.3評價(jià)4種索托節(jié)點(diǎn)模擬方法

1) 自由度耦合方法在一定程度上簡化了索托節(jié)點(diǎn)問題，并被一些實(shí)際工程設(shè)計(jì)所采用，但是由于沒有考慮索單元與主體結(jié)構(gòu)之間的非線性摩擦力，計(jì)算分析方法略顯粗糙，與實(shí)際情況有一定差別。

2) 自由度耦合和變剛度彈簧單元來處理滑移索摩擦問題的方法，彈簧采用Combin39單元，可以模擬不同的摩擦系數(shù)。該方法理論明確，計(jì)算易收斂，并在2008年奧運(yùn)會羽毛球館的實(shí)際工程算例中得到驗(yàn)證，是目前為止較為合理的方法[13]。

圖11　四種索托節(jié)點(diǎn)模擬方法計(jì)算的平直桿內(nèi)力值Fig.11　Internal forces of straight rod calculated by four methods

圖12　四種索托節(jié)點(diǎn)模擬方法計(jì)算的平斜桿內(nèi)力值Fig.12　Internal forces of inclined rod calculated by four methods

3) 點(diǎn)-點(diǎn)接觸法屬于非線性接觸分析，一般的通用有限元軟件可以提供接觸單元，利用接觸單元處理接觸滑移，但這種方法考慮非線性接觸行為，計(jì)算工作量大，計(jì)算方法較復(fù)雜，適用于滑移索的單個(gè)節(jié)點(diǎn)分析，不適用于帶滑移索的整體結(jié)構(gòu)體系分析，不便于一般工程技術(shù)人員掌握。

4) 短梁方法：通過本文4種方法的對比和分析可知短梁方法是可行的，可以通過改變短梁的長度和截面面積來調(diào)整索托節(jié)點(diǎn)處拉索的摩擦損失力。短梁剛度不變，計(jì)算更簡單，計(jì)算效率較高，并且該方法的理論明確，工程人員更容易理解和掌握。但這種方法需要在拉索的張拉力和摩擦情況已知的情況下才能確定短梁的幾何參數(shù)。

4　結(jié)論

1) 本文通過對4種索托節(jié)點(diǎn)數(shù)值模擬方法(自由度耦合方法、自由度耦合加變剛度彈簧單元方法、點(diǎn)-點(diǎn)接觸方法和不變剛度短梁法)的分析計(jì)算可知：由于索托節(jié)點(diǎn)處的摩擦力小(約為張拉力的1% ～ 5%)，所以4種方法對桁架結(jié)構(gòu)撓度和構(gòu)件內(nèi)力引起的差別不是很明顯，結(jié)果基本吻合，內(nèi)力走勢一樣。其中，自由度耦合加變剛度彈簧單元方法在2008年奧運(yùn)會羽毛球館的實(shí)際工程算例中得到驗(yàn)證，進(jìn)而說明本文新提出的不變剛度短梁法是合理可行的。

2) 通過對4種索托節(jié)點(diǎn)模擬方法的評價(jià)可知：自由度耦合加彈簧方法可以模擬不同的摩擦系數(shù)；短梁方法可以通過改變短梁的長度和截面面積來調(diào)整索托節(jié)點(diǎn)處拉索的摩擦損失力。這兩種方法能更好地模擬索托節(jié)點(diǎn)處的摩擦力，理論明確且計(jì)算簡單方便，是相對較好的索托節(jié)點(diǎn)數(shù)值模擬方法。

[1]胡沛銘. 索托桁架結(jié)構(gòu)延性設(shè)計(jì)研究[D]. 北京: 北京交通大學(xué), 2014.

HU Peiming. Research on the ductility design of cable-supported truss structure[D]. Beijing: Beijing Jiaotong University, 2014.

[2]MALERBA P G, GALLI P, DI DOMIZIO M, et al. A large span roof made of cable stayed arches[C]//Proceedings of ARCH′10-the 6th International Conference on Arch Bridges. Fuzhou, China, 2010.

[3]陳丹. 索托結(jié)構(gòu)關(guān)鍵問題研究[D]. 北京: 北京交通大學(xué), 2010.

CHEN Dan. Study on key problems of cable-supported structure[D]. Beijing: Beijing Jiaotong University, 2010.

[4]張阿鵬. 索托結(jié)構(gòu)靜力和抗震性能的研究[D]. 西安: 西安建筑科技大學(xué), 2007.

ZHANG Apeng. Study on static and aseismic behaviors of cable-held structure[D]. Xi’an: Xi’an University of Architecture and Technology, 2007.

[5]楊維國, 鄒清, 李興坡, 等. 索托結(jié)構(gòu)中拉索預(yù)應(yīng)力水平的優(yōu)化研究[J]. 中國礦業(yè)大學(xué)學(xué)報(bào), 2012, 41(2): 219-224.

YANG Weiguo, ZOU Qing, LI Xingpo, et al. An optimized pre-stress level for the design of cable supported trusses[J]. Journal of China university of mining & technology, 2012, 41(2): 219-224.

[6]席根喜, 徐國彬. 索托結(jié)構(gòu)在蘆城體校曲棍球館的應(yīng)用[J]. 空間結(jié)構(gòu), 2010, 16(3): 60-63.

XI Genxi, XU Guobin. Application of cable-held structure in hockey stadium of Beijing Lucheng Physical Training School[J]. Spatial structures, 2010, 16(3): 60-63.

[7]侯俊. 北京蘆城體校曲棍球訓(xùn)練館索托結(jié)構(gòu)施工技術(shù)[J]. 工程施工技術(shù), 2007(6): 89-93.

HOU Jun. Construction technologies of cable-supported structure in Beijing Lucheng physical training school hockey training hall[J]. Construction technology, 2007(6): 89-93.

[8]谷鵬. 索托桁架結(jié)構(gòu)性能影響因素分析[D]. 北京: 北京交通大學(xué), 2011.

GU Peng. Study on the influencing factors of cable-supported truss structural performance[D]. Beijing: Beijing Jiaotong University, 2011.

[9]王新敏. ANSYS工程結(jié)構(gòu)數(shù)值分析[M]. 北京: 人民交通出版社, 2007: 6-14.

[10]SONG M K, KIM S H, CHOI C K. Enhanced finite element modeling for geometric non-linear analysis of cable-supported structures[J]. Structural engineering and mechanics, 2006, 22(5): 575-597.

[11]LI Fengchen, ZHANG Lina, TIAN Shizhu, et al. Analysis of identified cable force of cable-stayed bridge[C]//Proceedings of the 2nd International Conference on Transportation Engineering. Chengdu, China, 2009: 3020-3025.

[12]張國軍, 葛家琪, 谷鵬, 等. 預(yù)應(yīng)力大懸挑鋼結(jié)構(gòu)索-鋼轉(zhuǎn)換節(jié)點(diǎn)足尺模型試驗(yàn)與設(shè)計(jì)研究[J]. 建筑結(jié)構(gòu), 2010, 40(12): 34-40, 83.

ZHANG Guojun, GE Jiaqi, GU Peng, et al. Full-scale model experiment and research of large span prestressed cable-member transfer connection[J]. Building structure, 2010, 40(12): 34-40, 83.

[13]王樹, 張國軍, 葛家琪, 等. 2008奧運(yùn)會羽毛球館預(yù)應(yīng)力損失對結(jié)構(gòu)體系影響分析[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2007, 28(6): 45-51.

WANG Shu, ZHANG Guojun, GE Jiaqi, et al. Influence of prestress loss on structural behavior of the badminton gymnasium for 2008 Olympic Games[J]. Journal of building structures, 2007, 28(6): 45-51.

本文引用格式：

楊維國, 王萌, 王亞,等. 索托節(jié)點(diǎn)的數(shù)值模擬方法研究[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2016, 37(8): 1050-1056.

YANG Weiguo, WANG Meng, WANG Ya, et al. Study of numerical simulation methods for cable-supported joints[J]. Journal of Harbin Engineering University, 2016, 37(8): 1050-1056.

Study of numerical simulation methods for cable-supported joints

YANG Weiguo, WANG Meng, WANG Ya, RUAN Zhiyong

(School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China)

To more accurately simulate cable-supported joints, in this study, we simulated an overall truss by using ANSYS model, and carried out finite element method (FEM) analyses. As a new method, the constant stiffness short beam method was advanced for simulating cable supported joints. It was also compared with other three methods, the degree-of-freedom coupling method, the degree-of-freedom coupling with variable stiffness spring method, and the joint-joint contact method. The results show that: there are no obvious differences in the deformation and internal forces for the four methods, which are basically same. The trend of the internal forces is consistent. Because the free coupling with variable stiffness spring method has been verified, the proposed constant stiffness short beam method we proved is both reasonable and feasible. The degree-of-freedom coupling with variable stiffness spring method and the constant stiffness short beam method can better simulate the friction forces in cable-supported joints, which represents a better way to model cable-supported nodes, and the theory underlying these two methods is clear and convenient.

cable-supported joints; numerical simulation; constant stiffness short beam method; degree-of-freedom coupling method; variable stiffness spring; joint-joint contact method; friction force

2015-05-09.網(wǎng)絡(luò)出版日期：2016-06-23.

國家自然科學(xué)基金項(xiàng)目(51178041)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2011JBM260)；中國航空規(guī)劃建設(shè)發(fā)展有限公司項(xiàng)目(技13研-51).

楊維國(1973-), 男, 教授, 博士;

王亞(1990-), 女, 博士研究生.

王亞，E-mail:wangya_73@163.com.

10.11990/jheu.201505025

TU393.3

A

1006-7043(2016)08-1050-07