考慮縱向滑移的組合梁隨機(jī)有限元可靠度分析

賈布裕, 余曉琳, 顏全勝

(華南理工大學(xué) 土木與交通學(xué)院, 廣東 廣州 510640)

考慮縱向滑移的組合梁隨機(jī)有限元可靠度分析

賈布裕, 余曉琳, 顏全勝

(華南理工大學(xué) 土木與交通學(xué)院, 廣東 廣州 510640)

摘要：針對(duì)鋼-混組合梁構(gòu)造復(fù)雜、材料多樣、不確定性強(qiáng)的特點(diǎn),從隨機(jī)不確定性的角度來(lái)研究鋼-混組合梁，基于Newmark幾何模型,將鋼-混組合梁的混凝土板、鋼梁以及剪力連接件作為整體單元考慮,通過(guò)采用T.L列式增量法,提出并推導(dǎo)了考慮縱向滑移效應(yīng)的非線性組合梁?jiǎn)卧Ｐ?。結(jié)合直接微分法,對(duì)組合梁的結(jié)構(gòu)響應(yīng)梯度進(jìn)行了計(jì)算分析；并由FORM法計(jì)算得到可靠度指標(biāo),對(duì)鋼-混組合梁算例進(jìn)行了隨機(jī)有限元可靠度分析。計(jì)算結(jié)果表明：在實(shí)際工程中,并不需要完全剛性連接來(lái)保證結(jié)構(gòu)的可靠性；在保證安全方面,考慮幾何非線性來(lái)計(jì)算結(jié)構(gòu)的可靠度是有必要的。

關(guān)鍵詞：鋼-混組合梁；縱向滑移；非線性組合梁?jiǎn)卧恢苯游⒎址?；隨機(jī)有限元；可靠度分析；幾何非線性

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

網(wǎng)絡(luò)出版日期：2015-11-04.

余曉琳(1978-), 女, 副教授.

近二三十年來(lái),隨著鋼-混組合結(jié)構(gòu)在工程領(lǐng)域的廣泛應(yīng)用,學(xué)者們開(kāi)展了各種關(guān)于鋼-混組合梁的研究工作,取得了豐富的成果,主要集中在以下幾個(gè)方面：非線性分析[1-2]、縱向滑移效應(yīng)[3-7]、橋面板長(zhǎng)期時(shí)變效應(yīng)[8]、剪力連接件抗剪[9]以及有限元單元模擬[3 ,10]等。但以上的研究都局限于確定性的結(jié)構(gòu)分析,對(duì)鋼-混組合結(jié)構(gòu)的不確定性分析,如結(jié)構(gòu)隨機(jī)性分析則鮮有人涉及。結(jié)構(gòu)隨機(jī)性包括結(jié)構(gòu)材料本身物理性質(zhì)以及幾何特性的隨機(jī)不確定性和外界對(duì)其施加作用的隨機(jī)不確定性。鋼-混組合梁涉及多種材料,單元構(gòu)造復(fù)雜,隨機(jī)不確定性因素較多,如何將鋼-混組合梁確定性結(jié)構(gòu)的計(jì)算分析手段和用于分析隨機(jī)性的可靠度分析方法結(jié)合起來(lái),對(duì)將鋼-混組合梁進(jìn)行隨機(jī)有限元可靠度分析,是目前研究的難點(diǎn)之一。本文以推導(dǎo)的考慮縱向滑移效應(yīng)的非線性組合梁?jiǎn)卧Ｐ蜑榛A(chǔ),采用直接微分法,對(duì)鋼-混組合梁進(jìn)行隨機(jī)有限元可靠度分析,并對(duì)數(shù)值算例進(jìn)行了計(jì)算分析。

1有限元模型

1.1　幾何位移模式

如圖1所示的典型的鋼-混組合梁截面,頂板是混凝土橋面板,下托梁是工字鋼梁,兩者通過(guò)剪力連接件相連。圖中坐標(biāo)x軸為順橋向,z軸為橫橋向,y軸為豎向。計(jì)算采用如下基本假定：

1) 橫截面對(duì)稱(chēng)于yox平面；2) 不考慮組合梁的豎向掀起(即橋面板和鋼主梁的豎向撓度一致)；3) 忽略剪切變形影響；4) 剪力連接件沿梁均勻分布。

由于不考慮組合梁的豎向掀起影響,鋼主梁和橋面板的豎向撓度可用同一個(gè)豎向位移表示。

圖1　幾何模型以及單元位移示意圖Fig.1　Geometric model and element displacements

組合單元t時(shí)刻橋面板和鋼主梁截面上任一點(diǎn)的位移和位移增量可表示為

(1)

組合梁t時(shí)刻縱向滑移位移以及增量的表達(dá)式為

(2)

其中,h=ys-yc,值得注意的是,式(2)適合于描述幾何線性問(wèn)題,但從工程實(shí)踐運(yùn)用角度來(lái)看,橋面板和鋼主梁之間的滑移變形(比如轉(zhuǎn)動(dòng))被看作為是有限的,所以式(2)也適合用于描述組合梁幾何非線性的縱向滑移位移[5]。

1.2　單元模式

將組合梁沿軸向進(jìn)行一維離散,單元模式并沒(méi)有采用常規(guī)的8自由度法(兩端節(jié)點(diǎn)各4自由度),因?yàn)槿舨扇?自由度法,則u的形函數(shù)的插值為一次多項(xiàng)式,而v的形函數(shù)的插值為三次多項(xiàng)式(v′為二次), 就會(huì)發(fā)生“滑移鎖定”問(wèn)題[3]。為避免此類(lèi)問(wèn)題的發(fā)生, 可采用增加內(nèi)置自由度的方法。采取在單元內(nèi)部u位移量增置2自由度的方法,最后得到了10自由度的單元模式(如圖2所示)。t時(shí)刻單元節(jié)點(diǎn)位移向量以及節(jié)點(diǎn)位移增量向量為

組合梁的位移場(chǎng)為{uc,us,v},則t時(shí)刻單元內(nèi)的位移和增量通過(guò)形函數(shù)可表示為

其中,形函數(shù)N為

各個(gè)構(gòu)件格林應(yīng)變?cè)隽靠杀硎緸?/p>

(3)

由虛功原理結(jié)合基于T.L列式增量法得到平衡方程：

(4)

圖2　單元節(jié)點(diǎn)自由度示意圖Fig.2　Element nodal degrees of freedom

2組合梁結(jié)構(gòu)隨機(jī)有限元分析

在有限元分析過(guò)程中,首先由控制方程中得到主要反應(yīng)量,如單元節(jié)點(diǎn)位移反應(yīng)量,在得到主要反應(yīng)量后,單元應(yīng)變或應(yīng)力即可從主要反應(yīng)量通過(guò)顯式的表達(dá)式得到,而主要反應(yīng)量也是基本隨機(jī)變量的函數(shù),所以可將基于隨機(jī)有限元的極限狀態(tài)方程G寫(xiě)為主要反應(yīng)量d和基本隨機(jī)變量θ的函數(shù)：

(5)

(6)

(7)

式中:K為切線剛度矩陣。由式(5)可得極限狀態(tài)方程函數(shù)G對(duì)基本變量θ的梯度計(jì)算公式：

(8)

接著計(jì)算標(biāo)準(zhǔn)正態(tài)空間極限狀態(tài)函數(shù)G(u)的梯度向量：

式中:Jθ,u表示Nataf分布對(duì)獨(dú)立標(biāo)準(zhǔn)正態(tài)分布的Jacobi矩陣,θg則由求得,其中的和可直接通過(guò)功能函數(shù)形式獲得,而可通過(guò)式(7)計(jì)算得到的結(jié)構(gòu)主要反應(yīng)量d對(duì)任一隨機(jī)參數(shù)θ的梯度值集合向量獲得。最后求得uG后,采用改進(jìn)HL-RF迭代法,由FORM法計(jì)算得到一次可靠度指標(biāo)β和設(shè)計(jì)試驗(yàn)點(diǎn)u*。

為了避免“滑移鎖定”的發(fā)生,本文在單元內(nèi)部增置了內(nèi)部節(jié)點(diǎn)自由度,而在有限元求解過(guò)程中,需要通過(guò)靜力凝聚法將內(nèi)部節(jié)點(diǎn)的抵抗內(nèi)力集合到結(jié)構(gòu)抵抗力。靜力凝聚法通常用于有內(nèi)部節(jié)點(diǎn)的單元有限元求解,其目的是為了縮減有限元系統(tǒng)方程求解規(guī)模,根據(jù)靜力凝聚法,可知結(jié)構(gòu)的控制平衡方程為

(9)

進(jìn)行靜力凝聚法后式(9)可等效寫(xiě)為

通過(guò)成功救治的案例，我們回顧本例的治療經(jīng)過(guò)，復(fù)習(xí)相關(guān)文獻(xiàn)資料，我們認(rèn)為，急性化膿性中耳炎A族鏈球菌感染合并顱內(nèi)感染，需要多學(xué)科合作，由于患者多伴有基礎(chǔ)病，基礎(chǔ)病的治療應(yīng)給以足夠重視；一旦出現(xiàn)顱內(nèi)并發(fā)癥，應(yīng)當(dāng)及時(shí)行中耳乳突病灶的清除及術(shù)腔引流；抗感染治療聯(lián)合用藥，根據(jù)藥敏試驗(yàn)，選擇敏感的，能透過(guò)血腦屏障及的抗生素，靜脈給藥，療程足夠。

(10)

(11)

式(11)兩邊對(duì)di求偏導(dǎo) (固定隨機(jī)變量θ)：

(12)

將式(12)對(duì)隨機(jī)變量θ進(jìn)行偏導(dǎo),可得

(13)

同樣式(11)兩邊對(duì)隨機(jī)變量θ求偏導(dǎo)有

(14)

(15)

結(jié)合式(13)～(15),整理可得

(16)

假設(shè)整個(gè)單元處于彈性,而Pint為上文推導(dǎo)的單元內(nèi)力,其表達(dá)式為

(17)

式中：BLc、BLs為格林應(yīng)變?cè)隽康木€性表示,BN為格林應(yīng)變?cè)隽康姆蔷€性表示,Nf為滑移增量,式(17) 左右兩邊對(duì)隨機(jī)變量θ進(jìn)行求偏導(dǎo),最后整理可得

在此暫不考慮節(jié)點(diǎn)坐標(biāo)的隨機(jī)性,只先考慮隨機(jī)變量θ為材料參數(shù),按類(lèi)似前節(jié)的推導(dǎo)方法,可得

式中：F(x,θ)為根據(jù)上文所推導(dǎo)的組合梁?jiǎn)卧?jié)點(diǎn)截面力,即

式中：fc(x,θ)為上層混凝土軸向力,fs(x,θ)為下層鋼主梁軸向力,Mcs(x,θ)為整個(gè)組合截面的彎矩,ω(x,θ)為連接件單位長(zhǎng)度剪力,而以上各項(xiàng)需要通過(guò)截面積分得到,由于d固定條件等價(jià)于ε固定條件,則最后可以得到基于截面層次的梯度表達(dá)式：

上述兩式右邊項(xiàng)的條件求導(dǎo)以及連接件單位長(zhǎng)度剪力ω(x,θ)的條件求導(dǎo)的計(jì)算需要借助于材料層次的梯度分析。材料應(yīng)力對(duì)任一隨機(jī)變量θ進(jìn)行偏導(dǎo),可得

3數(shù)值算例分析

根據(jù)上述方法,采用C++語(yǔ)言以及Matlab程序,結(jié)合自編的組合梁確定性有限元程序編制了適用于組合結(jié)構(gòu)的隨機(jī)有限元分析程序,并對(duì)2個(gè)組合梁算例進(jìn)行了隨機(jī)有限元的可靠度計(jì)算分析。

3.1　簡(jiǎn)支組合梁算例

如圖3所示的簡(jiǎn)支組合梁,在上層構(gòu)件上作用著均布荷載,并在上下2層中心沿縱向方向分別作用一對(duì)大小相等方向相反的集中力。上下2種材料均為線彈性,上層構(gòu)件的彈模為12 GPa,下層構(gòu)件的彈模為8 GPa,連接件的剛度通過(guò)定義一個(gè)無(wú)量綱剛度系數(shù)αL來(lái)表示：

式中：L為結(jié)構(gòu)長(zhǎng)度,k為連接件剛度,E1和E2、A1和A2、I1和I2分別為上下2種構(gòu)件的彈模,面積以及截面慣性矩、h為沿截面高度上下2個(gè)構(gòu)件質(zhì)心之間的距離。選取上下2個(gè)構(gòu)件的面積、慣性矩、彈模和連接件的連接剛度以及作用在上層結(jié)構(gòu)的均布荷載作為隨機(jī)變量,隨機(jī)變量的特征如表1所示。

表1　隨機(jī)變量特征表Table 1　Characteristics of structural random variables

為了分析幾何非線性以及連接件剛度對(duì)可靠度的影響,對(duì)幾何線性、幾何非線性以及5種連接件剛度工況下的結(jié)構(gòu)可靠度進(jìn)行了計(jì)算分析,計(jì)算結(jié)果見(jiàn)表2和圖4。

表2　多種工況下的可靠度結(jié)果.Table 2　Reliability results of different conditions

圖4　考慮幾何線性、非線性及不同連接剛度的可靠度結(jié)果Fig.4　Reliability results of geometric linear, non-linear and different connection stiffness

從表2以及圖4中可知：幾何線性工況下的可靠度明顯比幾何非線性工況下的可靠度高,所以在保證安全方面,考慮結(jié)構(gòu)的幾何非線性來(lái)計(jì)算結(jié)構(gòu)的可靠度是有必要的；無(wú)論是幾何線性還是幾何非線性,結(jié)構(gòu)的可靠度隨著連接件剛度的變大而迅速增大,但當(dāng)連接件的剛度達(dá)到一定數(shù)值時(shí),連接件剛度的增加對(duì)結(jié)構(gòu)可靠度的增加影響甚小,換言之,只要連接件的剛度達(dá)到一定程度,就可以保證達(dá)到和完全剛性連接組合結(jié)構(gòu)基本一致的可靠度。所以在實(shí)際工程,并不需要完全剛性連接來(lái)保證結(jié)構(gòu)的可靠性,相反可以通過(guò)計(jì)算得出一定的連接剛度來(lái)設(shè)計(jì)連接件,這樣可有效減少材料費(fèi)用。

3.2　連續(xù)組合梁算例

如圖5所示,該結(jié)構(gòu)為一根兩跨組合連續(xù)梁,組合梁上部構(gòu)件為一塊混凝土板,下部構(gòu)件為一塊工字鋼,兩者通過(guò)連接件聯(lián)結(jié)在一起,各跨的跨中各作用豎向集中力P。材料為非線性,各構(gòu)件的材料本構(gòu)關(guān)系見(jiàn)圖5?？紤]材料本構(gòu)參數(shù)以及單跨跨中豎向集中荷載P的隨機(jī)性,對(duì)結(jié)構(gòu)進(jìn)行了可靠度分析,隨機(jī)變量特征見(jiàn)表3。單元?jiǎng)澐謹(jǐn)?shù)為14,設(shè)所有單元的fc1、εc1、fc2、εc2、fc4、εc4、fs、εsh、fmax、δh分別完全相關(guān)。功能函數(shù)定為：g=0.012-Δ (Δ為單跨跨中撓度)。

圖5　集中荷載作用下連續(xù)組合梁及材料本構(gòu)關(guān)系Fig.5　The continuous composite beam with concentrated force, and material constitutive relationships of components

表3　隨機(jī)變量特征表Table 3　Characteristics of structural random variables

首先考慮無(wú)軸向壓力的工況,采用所提方法計(jì)算得到的可靠度指標(biāo)為β= 3.890 6,失效概率為Pf=5.00×10-5。為了驗(yàn)證結(jié)果的正確性,采用MC法,抽樣106次進(jìn)行可靠度計(jì)算,MC法所得到的失效概率為Pf=1.00×10-5,可靠度指標(biāo)β=4.264 9,其結(jié)果和所提方法結(jié)果比較吻合。對(duì)有軸向壓力的工況同樣進(jìn)行了可靠度計(jì)算分析,軸向壓力P1仍然分別取為0、500、1 000、1 500 kN,其分別對(duì)應(yīng)的可靠度結(jié)果如圖6所示。正如圖6所示,由于軸向壓力的作用,結(jié)構(gòu)的撓度將增大,其對(duì)應(yīng)的可靠度將隨之減小。

圖6　不同軸向壓力的可靠度結(jié)果Fig.6　Reliability results of different axial forces

4結(jié)論

1)組合梁結(jié)構(gòu)構(gòu)造復(fù)雜,不確定性程度高,因此采用隨機(jī)可靠度思想進(jìn)行分析是非常有必要的；

2)首先建立了確定性的組合梁?jiǎn)卧Ｐ?能考慮縱向滑移效應(yīng)和非線性影響,更加真實(shí)、準(zhǔn)確地反映了組合梁的受力性能；

3)將基于直接微分法的隨機(jī)有限元法直接應(yīng)用于所推導(dǎo)的非線性組合梁?jiǎn)卧?得出了考慮縱向滑移效應(yīng)的組合梁結(jié)構(gòu)響應(yīng)梯度表達(dá)式。在其基礎(chǔ)上編制的組合梁隨機(jī)有限元程序能有效完成對(duì)組合梁斜拉橋的隨機(jī)有限元分析；

4)所提的隨機(jī)有限元法直接以確定性有限元為基礎(chǔ),相比目前較為流行的響應(yīng)面法,在精度和穩(wěn)定性方面具有獨(dú)特優(yōu)勢(shì)(尤其是復(fù)雜或非線性程度較高的結(jié)構(gòu))；

5)算例結(jié)果表明：連接件的剛度對(duì)提升組合梁結(jié)構(gòu)的可靠度具有重要影響,其隨著連接件的剛度變大而增大,但連接件的剛度達(dá)到一定界限值時(shí),可靠度不再增大,因此設(shè)計(jì)時(shí)需注意到此點(diǎn),無(wú)需過(guò)于保守,避免浪費(fèi)；計(jì)算組合梁的可靠度時(shí),應(yīng)考慮幾何非線性的影響,考慮幾何非線性得到的結(jié)果和考慮幾何線性得到的結(jié)果存在明顯差異。幾何非線性降低了組合梁結(jié)構(gòu)的可靠度值,這點(diǎn)在實(shí)際工程中需引起重視；

6)組合梁結(jié)構(gòu)是一個(gè)復(fù)雜的組合體,其失效模式眾多,如還有抗彎極限、豎向抗剪極限等問(wèn)題，今后還需完善這些方面的可靠度研究。

參考文獻(xiàn)：

[1]BATTINI J M, NGUYEN Q H, HJIAJ M. Non-linear finite element analysis of composite beams with interlayer slips[J]. Computers and Structures, 2009, 87(13/14): 904-912.

[2]DALL′ASTA A, ZONA A. Non-linear analysis of composite beams by a displacement approach[J]. Computers and Structures, 2002, 80(27-30): 2217-2228.

[3]DALL′ASTA A, ZONA A. Slip locking in finite elements for composite beams with deformable shear connection[J]. Finite Elements in Analysis and Design, 2004, 40(13/14): 1907-1930.

[4]DALL′ASTA A, ZONA A. Three-field mixed formulation for the non-linear analysis of composite beams with deformable shear connection[J]. Finite Elements in Analysis and Design, 2004, 40(4): 425-448.

[5]GIRHAMMAR U A, PAN D H. Exact static analysis of partially composite beams and beam-columns[J]. International Journal of Mechanical Sciences, 2007, 49(2): 239-255.

[6]SOUSA JR J B M, OLIVEIRA C E M, DA SILVA A R. Displacement-based nonlinear finite element analysis of composite beam-columns with partial interaction[J]. Journal of Constructional Steel Research, 2010, 66(6): 772-779.

[7]RANZI G, DALL′ASTA A, RAGNI L, et al. A geometric nonlinear model for composite beams with partial interaction[J]. Engineering Structures, 2010, 32(5): 1384-1396.

[8]FRAGIACOMO M, AMADIO C, MACORINI L. Finite-element model for collapse and long-term analysis of steel-concrete composite beams[J]. Journal of Structural Engineering, 2004, 130(3): 489-497.

[9]KIM B, WRIGHT H D, CAIRNS R. The behaviour of through-deck welded shear connectors: an experimental and numerical study[J]. Journal of Constructional Steel Research, 2001, 57(12): 1359-1380.

[10]ZONA A, RAGNI L, DALL′ASTA A. Finite element formulation for geometric and material nonlinear analysis of beams prestressed with external slipping tendons[J]. Finite Elements in Analysis and Design, 2008, 44(15): 910-919.

[11]李偉. 大跨度斜拉橋靜動(dòng)力可靠度分析[D]. 廣州: 華南理工大學(xué), 2010: 53-55.LI Wei. Structural static and seismic reliability analysis of long-span cable-stayed bridges[D]. Guangzhou: South China University of Technology, 2010: 53-55.

Stochastic finite element reliability analysis of

composite beams considering longitudinal slip

JIA Buyu, YU Xiaolin, YAN Quansheng

(School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, China)

Abstract：Steel-concrete composite beams are characterized by complex structure, material variety, and strong uncertainty, thus it is necessary to analyze them from the view of randomness. Based on the Newmark geometric model, considering the concrete slab, steel beams, and shear connectors as a whole element,and making use of the T.L. column incremental method, a nonlinear composite beam element model was developed by taking the longitudinal slip effect into consideration. Then, based on the model, and the characteristics of composite beams, structural response gradients of composite beams were obtained by the direct differentiation method (DDM), while the reliability index was achieved using the FORM method. For proof of reliability, an example steel-concrete composite beam was analyzed by the stochastic finite element program. The results show that, for practical engineering, it does not need a completely rigid connection to ensure structural reliability. However for security, it is necessary to consider geometric nonlinearity in the calculation of structure reliability.

Keywords：steel-concrete composite beam; longitudinal slip; nonlinear composite beam element; direct differentiation method; stochastic finite element; reliability analysis; geometric nonlinearity

通信作者：余曉琳,E-mail:xlyul@scut.edu.cn.

作者簡(jiǎn)介：賈布裕(1983-), 男, 博士后;

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目(51208208)；中國(guó)博士后科學(xué)基金資助項(xiàng)目(2013M542174)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金資助項(xiàng)目(2014ZZ0019,,2015ZM114).

收稿日期：2014-04-07.

中圖分類(lèi)號(hào)：U448.27

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1006-7043(2015)12-1554-06

doi:10.11990/jheu.201404024