裴云亮柳鋒任參邴卿德

(1.山東建筑大學(xué) 土木工程學(xué)院，山東 濟南250101；2.山東建大建筑規(guī)劃設(shè)計研究院，山東 濟南250000)

0 引言

網(wǎng)殼結(jié)構(gòu)屬于缺陷敏感型結(jié)構(gòu)，初始缺陷對其穩(wěn)定極限承載力的影響較復(fù)雜。 初始缺陷可分為整體缺陷和桿件缺陷，其中整體缺陷主要指節(jié)點位置偏差，而桿件缺陷包括桿件初彎曲、內(nèi)力初偏心以及殘余應(yīng)力等。 MORRIS[1]發(fā)現(xiàn)初始幾何缺陷可使結(jié)構(gòu)的臨界荷載下降35%，穩(wěn)定分析中需考慮缺陷的影響；對于考慮節(jié)點偏差和不考慮節(jié)點偏差的結(jié)構(gòu)，桿件初彎曲將明顯降低結(jié)構(gòu)承載力，并且能夠改變結(jié)構(gòu)的塑性發(fā)展程度[2-3]；與僅考慮結(jié)構(gòu)整體初始幾何缺陷的計算結(jié)果相比，桿件初彎曲會進一步降低結(jié)構(gòu)的穩(wěn)定承載力，且當(dāng)桿件初彎曲幅值較大時，桿件初彎曲對結(jié)構(gòu)穩(wěn)定承載力的影響不容忽視[4]；馬騰飛等[5]在考慮節(jié)點缺陷和桿件缺陷的基礎(chǔ)上，通過削減不同位置處的桿件截面，使被削減的桿件在網(wǎng)殼結(jié)構(gòu)整體失穩(wěn)之前屈曲，不同位置處的桿件屈曲會不同程度的降低網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載力；趙陽等[6]和田偉等[7-8]對網(wǎng)殼結(jié)構(gòu)整體穩(wěn)定分析時發(fā)現(xiàn)，結(jié)構(gòu)整體失穩(wěn)之前存在桿件失穩(wěn)的網(wǎng)殼結(jié)構(gòu)，有必要引入桿件缺陷，考慮其對桿件失穩(wěn)的影響，若不考慮桿件失穩(wěn)會過高估計結(jié)構(gòu)的承載能力，局部構(gòu)件的失穩(wěn)退出并不會導(dǎo)致結(jié)構(gòu)整體承載力的突然喪失，結(jié)構(gòu)仍會繼續(xù)承載。

為了不影響結(jié)構(gòu)的安全和正常使用，規(guī)范對上述缺陷限值是有明確規(guī)定的。 在實際的大跨網(wǎng)殼結(jié)構(gòu)中不可避免的存在桿件幾何缺陷，文章以某大跨度網(wǎng)殼結(jié)構(gòu)實際工程為研究對象，基于對結(jié)構(gòu)中彎曲桿件的實測數(shù)據(jù)，研究幾何缺陷對實際網(wǎng)殼結(jié)構(gòu)穩(wěn)定承載性能的影響。

1 實際大跨網(wǎng)殼結(jié)構(gòu)現(xiàn)場檢測和設(shè)計復(fù)核

某網(wǎng)殼結(jié)構(gòu)是一個跨度為151.8 m、矢高為45 m的雙層球面網(wǎng)殼，此網(wǎng)殼共有10 302 根桿件，桿件材料為Q235B 鋼，共有2 647 個節(jié)點，節(jié)點采用的是螺栓球節(jié)點，網(wǎng)殼結(jié)構(gòu)如圖1 所示。 在使用過程中，發(fā)現(xiàn)該網(wǎng)殼結(jié)構(gòu)存在部分彎曲的桿件，如圖2所示。 為了解桿件彎曲對網(wǎng)殼結(jié)構(gòu)性能的影響，對此網(wǎng)殼結(jié)構(gòu)進行了現(xiàn)場檢測和設(shè)計復(fù)核工作。

圖1 網(wǎng)殼結(jié)構(gòu)圖

圖2 桿件彎曲圖

1.1 桿件尺寸與變形檢測

采用三維激光掃描儀、游標(biāo)卡尺等測量工具對網(wǎng)殼進行檢測。 網(wǎng)殼結(jié)構(gòu)共有11 種桿件規(guī)格，分別為Φ75.5×3.75、Φ88.5×4、Φ114×4、Φ140×4、Φ159×6、Φ159×8、Φ159×10、Φ180×10、Φ180×12、Φ219×10、Φ219×12，桿件的幾何尺寸、截面規(guī)格與設(shè)計圖紙相符，尺寸偏差滿足要求[9]。 在對網(wǎng)殼結(jié)構(gòu)變形檢測時，發(fā)現(xiàn)其中有16 根桿件產(chǎn)生撓曲變形，這16 根桿件實際分布位置如圖3 所示，桿件各參數(shù)見表1。

圖3 彎曲桿件位置示意圖

表1 彎曲桿件參數(shù)表

1.2 設(shè)計復(fù)核

采用3d3s 鋼結(jié)構(gòu)設(shè)計軟件設(shè)計復(fù)核了網(wǎng)殼結(jié)構(gòu)，此結(jié)構(gòu)在建模時支座采用固定鉸支座。 對結(jié)構(gòu)進行彈性靜力分析和彈性全過程穩(wěn)定分析，網(wǎng)殼結(jié)構(gòu)荷載取值如下:網(wǎng)殼上、下弦永久荷載分別為0.5 和0.1 kN/m2；活荷載、風(fēng)荷載和雪荷載均為0.5 kN/m2；溫度荷載為±40 ℃；地震荷載為地震設(shè)防烈度7 度(0.15g)，主要荷載組合為:(1) 恒荷載＋活荷載；(2) 1.20 恒荷載＋1.40 活荷載；(3) 1.35 恒荷載＋140×0.70 活荷載；(4) 1.20 恒荷載＋140×0.50活荷載＋1.30 水平地震；(5) 1.20 恒荷載＋1.40 風(fēng)荷載；(6) 1.20 恒荷載＋1.40 活荷載＋1.40×0.60 風(fēng)荷載；(7) 1.20 恒荷載＋1.40×0.70 活荷載＋1.40 風(fēng)荷載；(8) 1.20 恒荷載＋1.40 活荷載＋1.40×0.60 溫度荷載。

對設(shè)計結(jié)構(gòu)在上述8 種荷載組合下分別進行計算，復(fù)核結(jié)果如下:桿件最大強度應(yīng)力比為0.89，最大穩(wěn)定應(yīng)力比為0.98，結(jié)構(gòu)最大撓跨比為1/1 494，整體穩(wěn)定系數(shù)為24，因此網(wǎng)殼結(jié)構(gòu)在設(shè)計荷載作用下，桿件最大應(yīng)力比、結(jié)構(gòu)最大撓跨比以及安全系數(shù)均滿足要求[10]，結(jié)構(gòu)具有良好的設(shè)計承載力，桿件出現(xiàn)彎曲變形不是因為設(shè)計承載力不足引起的。 為探討彎曲的桿件對大跨網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載能力的影響，對網(wǎng)殼結(jié)構(gòu)進行穩(wěn)定極限承載分析。

2 缺陷模態(tài)的確定

目前，在鋼結(jié)構(gòu)整體穩(wěn)定分析中主要考慮結(jié)構(gòu)幾何缺陷的影響，選用缺陷模態(tài)法引入初始幾何缺陷，主要有隨機缺陷模態(tài)法和一致缺陷模態(tài)法兩種方法[11]。 隨機缺陷模態(tài)法能夠較好的反映結(jié)構(gòu)的實際工作性能，但計算量大，不適合網(wǎng)殼結(jié)構(gòu)非線性分析；一致缺陷模態(tài)法概念清晰、計算簡便，認(rèn)為缺陷模態(tài)是結(jié)構(gòu)失穩(wěn)時的位移傾向，有利于了解結(jié)構(gòu)的整體穩(wěn)定性，因此選用一致缺陷模態(tài)法確定結(jié)構(gòu)的幾何缺陷分布。

2.1 桿件缺陷

假定上述16 根彎曲桿件的形狀為正弦半波曲線，由式(1)表示為

式中l(wèi)為桿件長度，mm；δ為彎曲幅值，其值取桿件彎曲實測值，彎曲方向為實測方向。

正弦半波曲線是一種理想化的簡單彎曲形狀，采用多段梁法模擬桿件彎曲[7-8，12]，其余鋼桿件按理想直桿模擬。

在有限元分析中，桿件的劃分段數(shù)會影響到計算結(jié)果的精確性，文獻[4，7，13]是通過考察單桿有限元分析模型來確定劃分段數(shù)，由于桿件與結(jié)構(gòu)之間的剛度相互影響，為了更加符合結(jié)構(gòu)真實受力狀態(tài)，將彎曲桿件劃分單元后引入到整體結(jié)構(gòu)中進行線性屈曲分析，并通過比較鄰近特征值變化的大小來確定桿件合理的劃分段數(shù)。 將分析模型中彎曲的桿件分別劃分為1、2、4、6、8、10、12 段，各模型在第1 階屈曲模態(tài)的特征值如圖4 所示。 隨著桿件劃分單元數(shù)量的增加，在第1 階屈曲模態(tài)中，特征值越來越接近，當(dāng)單元數(shù)量＞8 個時，特征值相差不大，考慮到桿件劃分的單元越多結(jié)構(gòu)計算量越大，因此在分析結(jié)構(gòu)時對桿件采用8 單元的劃分。

圖4 第1 階屈曲模態(tài)特征值圖

2.2 整體缺陷

對網(wǎng)殼全過程分析時，初始幾何缺陷分布可采用無缺陷結(jié)構(gòu)的最低階屈曲模態(tài)[10]，但由于特征值屈曲分析是線性分析，剛度矩陣是建立在結(jié)構(gòu)未受載時的初始結(jié)構(gòu)上，因而求得的最低階屈曲模態(tài)不一定反映結(jié)構(gòu)在整個非線性分析過程中的變形趨勢[14]，尤其考慮到桿件彎曲的影響，這就可能不是結(jié)構(gòu)最不利的初始缺陷分布形式；文獻[15]提到大跨空間結(jié)構(gòu)在失穩(wěn)前可能已經(jīng)產(chǎn)生很強的幾何非線性甚至材料非線性，導(dǎo)致結(jié)構(gòu)在發(fā)生整體失穩(wěn)時的變形模態(tài)已不同于第1 階線性屈曲模態(tài)。 因此將特征值屈曲分析得到的多階屈曲模態(tài)分別作為缺陷分布形式引入模型中進行幾何非線性分析，通過比較極限承載力系數(shù)(結(jié)構(gòu)穩(wěn)定極限承載力與荷載標(biāo)準(zhǔn)值的比值)大小確定最不利的整體缺陷。

2.3 最不利缺陷模態(tài)

將帶有彎曲的桿件引入到模型中，對網(wǎng)殼結(jié)構(gòu)進行線性屈曲分析，由于屈曲模態(tài)是網(wǎng)殼結(jié)構(gòu)失穩(wěn)時的位移傾向，對于跨度較大且形式復(fù)雜的網(wǎng)殼結(jié)構(gòu)來說，高階屈曲模態(tài)的影響不可忽略，選取網(wǎng)殼前30 階屈曲模態(tài)作為初始幾何缺陷分布，缺陷幅值取網(wǎng)殼實測撓度的最大值，通過修正初始網(wǎng)殼模型的節(jié)點坐標(biāo)，形成30 個模型樣本，然后分別對這30 個模型進行幾何非線性分析，并通過比較極限承載力系數(shù)的大小，確定最不利的缺陷模態(tài)。 經(jīng)計算得出第23 階屈曲模態(tài)作為初始幾何缺陷分布時對應(yīng)的極限承載力系數(shù)最小，所以第23 階屈曲模態(tài)是最不利的缺陷模態(tài)，如圖5 所示，網(wǎng)殼結(jié)構(gòu)發(fā)生局部凹陷現(xiàn)象。

圖5 最不利缺陷模態(tài)圖/m

3 帶桿件彎曲的網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載分析

為了研究桿件彎曲對網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載力的影響，分別建立以下3 種網(wǎng)殼分析模型:模型1 不考慮任何初始幾何缺陷；模型2 將結(jié)構(gòu)第1 階屈曲模態(tài)作為整體初始幾何缺陷，缺陷幅值取L/300(L為網(wǎng)殼跨度)；模型3 考慮彎曲桿件，采用最不利缺陷模態(tài)建立模型。

3.1 穩(wěn)定極限承載力

采用有限元軟件ANSYS，建立網(wǎng)殼有限元模型，桿件選用空間梁單元beam188，材料為Q235B鋼，其本構(gòu)關(guān)系采用理想彈塑性模型，分別對3 種模型進行幾何、材料雙重非線性分析，圖6 給出了這3 種模型的極限承載力系數(shù)—豎向位移曲線，極限承載力系數(shù)和相應(yīng)豎向位移見表2。

圖6 極限承載力系數(shù)—豎向位移曲線圖

表2 各模型極限荷載表

模型2 考慮整體缺陷，相比于模型1，其穩(wěn)定極限承載力下降了34.1%；模型3 考慮桿件彎曲缺陷與整體缺陷，相比于模型2，其穩(wěn)定極限承載力又下降了28.4%，可見桿件彎曲缺陷與整體缺陷對網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載力影響較大。

3.2 失穩(wěn)機理

為研究彎曲桿件對網(wǎng)殼結(jié)構(gòu)失穩(wěn)過程的影響，對模型3 在加載過程中發(fā)生失穩(wěn)的桿件進行路徑跟蹤，確定桿件失穩(wěn)的順序及桿件軸力的變化，探討網(wǎng)殼結(jié)構(gòu)的失穩(wěn)機理。 桿件失穩(wěn)順序及桿件失穩(wěn)時的荷載值見表3，軸力負(fù)號表示桿件受壓[16]。

在加載初期，網(wǎng)殼桿件的軸力和節(jié)點位移均隨荷載的增加呈線性增長，當(dāng)極限承載力系數(shù)增加到4.20 時，15 號桿件發(fā)生失穩(wěn)，同時觀察到15 號相鄰兩根桿件極限承載力系數(shù)—軸力曲線的斜率發(fā)生了變化，可見15 號桿件失穩(wěn)后，向其相鄰桿件卸載，結(jié)構(gòu)內(nèi)力重分布，由圖6 可知，此時網(wǎng)殼結(jié)構(gòu)的極限承載力系數(shù)—豎向位移曲線處于上升階段，結(jié)構(gòu)仍能繼續(xù)承載。 15 號失穩(wěn)桿件以及相鄰桿件的極限承載力系數(shù)—軸力關(guān)系曲線和局部放大圖如圖7所示。

表3 失穩(wěn)桿件及失穩(wěn)時荷載與內(nèi)力表

圖7 15 號失穩(wěn)桿件的軸力變化及放大圖

當(dāng)極限承載力系數(shù)增加到4.31 時，7 和9 號兩根桿件同時失穩(wěn)，此時7 和9 號相鄰桿件極限承載力系數(shù)—軸力曲線的斜率發(fā)生了變化，可見兩個桿件失穩(wěn)后，向其相鄰桿件卸載，結(jié)構(gòu)內(nèi)力重分布。 由圖6 可知，此時網(wǎng)殼結(jié)構(gòu)的極限承載力系數(shù)—豎向位移曲線處于上升階段，結(jié)構(gòu)仍能繼續(xù)承載。 7 和9號失穩(wěn)桿件以及相鄰桿件的極限承載力系數(shù)—軸力關(guān)系曲線如圖8 所示。

圖8 7 和9 號失穩(wěn)桿件的軸力變化及放大圖

當(dāng)極限承載力系數(shù)增加到4.42 時，2 號桿件發(fā)生失穩(wěn)，同時觀察到其相鄰兩根桿件極限承載力系數(shù)—軸力曲線的斜率發(fā)生了變化，可見2 號桿件失穩(wěn)后，向其相鄰桿件卸載，結(jié)構(gòu)內(nèi)力重分布。 由圖6可知，此時網(wǎng)殼結(jié)構(gòu)的極限承載力系數(shù)—豎向位移曲線處于上升階段，結(jié)構(gòu)仍能繼續(xù)承載。 2 號失穩(wěn)桿件以及相鄰桿件的極限承載力系數(shù)—軸力關(guān)系曲線和局部放大圖如圖9 所示。

圖9 2 號失穩(wěn)桿件的軸力變化及放大圖

當(dāng)極限承載力系數(shù)增加到4.51 時，11 號桿件發(fā)生失穩(wěn)，同時觀察到其相鄰兩根桿件極限承載力系數(shù)—軸力曲線的斜率發(fā)生了變化，可見11 號桿件失穩(wěn)后，向其相鄰桿件卸載，結(jié)構(gòu)內(nèi)力重分布。 由圖6可知，此時網(wǎng)殼結(jié)構(gòu)的極限承載力系數(shù)—豎向位移曲線處于上升階段，結(jié)構(gòu)仍能繼續(xù)承載。 11 號失穩(wěn)桿件以及相鄰桿件的極限承載力系數(shù)—軸力關(guān)系曲線和局部放大圖如圖10 所示。

圖10 11 號失穩(wěn)桿件的軸力變化及放大圖

當(dāng)極限承載力系數(shù)增加到4.59 時，8 和10 號兩根桿件同時發(fā)生失穩(wěn)，同時觀察到其相鄰桿件極限承載力系數(shù)—軸力曲線的斜率發(fā)生了變化，其軸力不再隨荷載線性增加，可見兩桿件失穩(wěn)后，向其相鄰桿件卸載，結(jié)構(gòu)內(nèi)力重分布。 由圖6 可知，此時網(wǎng)殼結(jié)構(gòu)的極限承載力系數(shù)—豎向位移曲線處于上升階段，結(jié)構(gòu)仍能繼續(xù)承載。 8 和10 號失穩(wěn)桿件以及相鄰桿件的極限承載力系數(shù)—軸力關(guān)系曲線和局部放大圖如圖11 所示。

當(dāng)極限承載力系數(shù)增加到4.73 時，網(wǎng)殼節(jié)點位移較大，結(jié)構(gòu)整體出現(xiàn)偏移，發(fā)生整體失穩(wěn)破壞，各桿件的軸力均隨著荷載的卸載而減小，如圖12所示。

由上述分析可知，網(wǎng)殼結(jié)構(gòu)的最終破壞源于帶有彎曲桿件的率先失穩(wěn)，桿件的失穩(wěn)使結(jié)構(gòu)發(fā)生內(nèi)力重分布，結(jié)構(gòu)不會突然喪失承載力，仍能繼續(xù)承載，隨著荷載的增加，失穩(wěn)桿件數(shù)量的增多，最終導(dǎo)致結(jié)構(gòu)整體失穩(wěn)。

圖11 8 和10 號失穩(wěn)桿件的軸力變化及放大圖

圖12 網(wǎng)殼結(jié)構(gòu)達到極限承載時的變形圖/m

4 結(jié)論

通過對帶有彎曲桿件的網(wǎng)殼結(jié)構(gòu)進行幾何、材料雙重非線性承載力分析，得出以下主要結(jié)論:

(1) 網(wǎng)殼結(jié)構(gòu)的最終破壞源于帶有彎曲桿件的率先失穩(wěn)，桿件失穩(wěn)后向其相鄰桿件卸載，致使結(jié)構(gòu)內(nèi)力發(fā)生重分布，結(jié)構(gòu)不會突然喪失整體承載力而是繼續(xù)承載。

(2) 網(wǎng)殼結(jié)構(gòu)考慮了桿件彎曲缺陷與整體缺陷，較僅考慮整體缺陷時的承載力下降了28.4%，桿件彎曲缺陷對網(wǎng)殼結(jié)構(gòu)穩(wěn)定極限承載力的影響較大。