楊 磊，薛曉宏，姜繼果，黃正業(yè)，姜子麒

(中鐵二十局集團第六工程有限公司，陜西 西安 710032)

0 引言

隨著我國經(jīng)濟的發(fā)展，城市化進程不斷加速，土地資源日趨緊缺，超高層、大跨度建筑開始受到青睞，由于鋼管混凝土具有承載力高、塑性和韌性好、施工方便、耐火性能好、經(jīng)濟效益好、輕質(zhì)化等優(yōu)點[1-2]，且矩形鋼管混凝土柱變形能力強、截面慣性矩大、節(jié)點構(gòu)造簡單及使用性強[3]，因此被廣泛應用于現(xiàn)代化建筑結(jié)構(gòu)中。理想狀態(tài)下，鋼管和混凝土應緊密接觸，共同承受外部荷載；但實際施工過程中存在混凝土振搗不密實，鋼管內(nèi)部有加勁肋、內(nèi)隔板及后期混凝土養(yǎng)護不合理等因素的影響，不可避免地會導致核心混凝土內(nèi)部產(chǎn)生缺陷[4-5]。

自1957年 Kloppel等[6]開始研究圓鋼管混凝土軸心受壓柱，并探討該類構(gòu)件承載力的計算方法起，對鋼管混凝土力學性能的研究越來越多[7-11]。目前對帶缺陷鋼管混凝土柱承載能力的研究中，對在缺陷隨機分布下的情況研究較少。葉勇[12]設置上、下貫通的均勻脫空缺陷，通過改變脫空率研究缺陷對軸壓承載能力的影響。張戎令等[13]設置的缺陷位置只考慮中心缺陷和邊側(cè)缺陷。劉開[14]在設置缺陷時，考慮了角部空洞、脫空、中央長方體、中央正方體和中央圓柱體等缺陷類型，并根據(jù)不同缺陷率設置不同的缺陷尺寸。門有冬[15]分別研究了孔洞缺陷、環(huán)形脫空缺陷和球冠形脫空缺陷對鋼管混凝土承載力的影響。李永進等[16]以脫空率和脫空類型為主要參數(shù)研究長期荷載對帶脫空缺陷混凝土受力性能的影響，其中脫空類型有通長均勻脫空缺陷和通長球冠形脫空缺陷。張傳欽等[17]以軸壓比、彎扭比為主要試驗參數(shù)，研究脫空缺陷對壓彎扭構(gòu)件力學性能的影響，其設計的缺陷類型為球冠形缺陷，以脫空處核心混凝土邊緣到鋼管內(nèi)壁的距離與構(gòu)件截面外直徑的比值來定義缺陷率。廖飛宇等[18]設置環(huán)向脫空缺陷研究鋼管混凝土構(gòu)件在壓彎扭復合受力作用下的滯回性能。

綜上所述，現(xiàn)階段對帶缺陷鋼管混凝土柱力學性能的研究中，模擬缺陷時人為干擾成分較多，忽略了缺陷實際存在的隨機性。本文基于ABAQUS/Python的二次開發(fā)，實現(xiàn)了不同缺陷大小與位置的隨機生成及結(jié)果自動提取；得到不同隨機缺陷類型對應的結(jié)構(gòu)承載力，從而總結(jié)缺陷對結(jié)構(gòu)承載力的影響規(guī)律。

1 建立有限元模型

1.1 材料本構(gòu)關系選擇

本文所建模型中鋼材的本構(gòu)模型選用理想型彈塑性模型，采用Von-Mises屈服準則，其應力-應變關系曲線如圖1所示[3]。核心混凝土本構(gòu)模型采用混凝土損傷塑性模型，該模型可以很好地模擬混凝土構(gòu)件的塑性與非線性特性，其假定材料的彈性行為是各向同性和線性的，而單軸拉伸和壓縮損傷塑性部分的力學行為由損傷塑性描述。

圖1 鋼管本構(gòu)關系

設計模型的基本參數(shù)為：方鋼管邊長為400mm，壁厚6mm，長1.2m；上、下受力剛性板邊長為440mm，厚50mm?；炷练侨毕菁蠁卧弯摴艿闹饕嬎銋?shù)為：混凝土密度2 600kg/m3，彈性模量3.7×104MPa，泊松比0.2；鋼管屈服強度330MPa，彈性模量2.1×105MPa，泊松比0.3。缺陷區(qū)域的單元假定為正常彈性模量的1/20，即為1.5×103MPa，端部受力剛性板的彈性模量取鋼材的100倍，不考慮其非線性。依據(jù) GB 50010—2010《混凝土結(jié)構(gòu)設計規(guī)范》中相關公式和參數(shù)及相關文獻計算核心混凝土塑性損傷模型參數(shù)，膨脹角為40°，偏心率為0.1，fb0/fc0為1.16，不變量應力比K為0.667，黏滯系數(shù)為0.000 5。依照模型尺寸確定混凝土各頂點處的坐標值：xbot=-197，xup= 197，ybot=-197，yup= 197，zbot= 0，zup= 1 200(見圖2)。

圖2 鋼管混凝土柱核心區(qū)混凝土模型

1.2 單元選取

鋼管的厚度尺寸相較于整體尺寸很小，其計算單元采用4節(jié)點減縮積分格式的殼單元來模擬，殼單元厚度方向采用9個積分點的Simpson積分方法，核心混凝土和兩端受力板采用8節(jié)點減縮積分格式的三維實體單元。

1.3 網(wǎng)格劃分

考慮到在分析接觸時，剛度更大的鋼管為主表面，核心混凝土為從屬表面，從屬表面的網(wǎng)格密度不能低于主表面[3]。為確定結(jié)果的網(wǎng)格穩(wěn)定性，采用試算法確定網(wǎng)格密度，模型建立后，先采用較粗糙的網(wǎng)格進行模擬，然后不斷精細網(wǎng)格布置，直至最終確定，網(wǎng)格單元數(shù)-承載力曲線如圖3所示。

圖3 網(wǎng)格單元數(shù)-承載力曲線

1.4 界面接觸

鋼管混凝土模型的相互作用采用表面與表面接觸，定義接觸屬性時，法線方向的接觸采用“硬”接觸，切線方向摩擦公式定義為“罰”[19]，當達到極限臨界摩擦值時，界面間開始發(fā)生相對滑移，伴隨剪應力傳遞停止。

1.5 加載方式與邊界條件

在鋼管柱兩端設置大剛度受力板以提高計算的收斂性，設置其與鋼管混凝土為綁定連接。在下部受力板施加三向位移約束。加載方式為在上部剛性受力板上設置參考點，將參考點與上部鋼板耦合，對特征點進行位移加載。

2 隨機有限元二次開發(fā)

構(gòu)建缺陷時，主要考慮缺陷的大小和位置，假定缺陷的形狀為球形，因此主要構(gòu)件參數(shù)有球體半徑大小和圓心坐標?；诖耍肞ython語言創(chuàng)建隨機缺陷的步驟如圖4所示。

圖4 隨機缺陷生成流程

基于此開發(fā)鋼管混凝土有限元模型混凝土缺陷隨機生成軟件，操作界面及生成缺陷如圖5所示。在此操作界面上可選擇設置缺陷的形狀，模擬次數(shù)，缺陷模式，混凝土模型在三維空間中x，y，z方向的上、下限，生成隨機缺陷位置和形狀的模式與參數(shù)。設置完成后，點擊繪圖可在繪圖區(qū)生成隨機缺陷模型，點擊導出，可自動生成缺陷模型腳本文件。

圖5 隨機缺陷生成軟件界面及生成缺陷

3 計算結(jié)果分析

3.1 承載力結(jié)果統(tǒng)計分析

在ABAQUS中建立好鋼管混凝土模型后，在缺陷生成軟件中確定好缺陷的設計參數(shù)，然后自動生成Python腳本，在ABAQUS中運行腳本，進行隨機有限元計算。鑒于最終計算得到的結(jié)果較多且篇幅有限，截取部分隨機缺陷數(shù)據(jù)如表1所示。

表1 部分隨機缺陷數(shù)據(jù)

計算共分2組，對應不同缺陷半徑的均值和方差，第1組均值為60，方差為40；第2組均值為30，方差為10。2組缺陷對應的鋼管混凝土有限元模型承載力頻率直方圖及得到擬合的曲線如圖6所示，承載力數(shù)值統(tǒng)計特征如表2所示。用K-S法對2組數(shù)據(jù)進行正態(tài)性檢驗，確定其分布是否滿足正態(tài)分布，顯著性水平為0.05，檢驗結(jié)果為：①第1組 均值為2 817.85kN，標準差為225.61，P=0.000 181(P表示數(shù)據(jù)的顯著性)；②第2組 均值為 3 008.15kN，標準差為66.04，P=0.000 002；兩組數(shù)據(jù)的顯著性均<0.05，不滿足正態(tài)分布。

圖6 承載力頻率直方圖

表2 承載力統(tǒng)計特征

由圖6可知，曲線最高點均偏向于x軸右邊，位于右半部分曲線比正態(tài)分布曲線更陡，而左半部分曲線較平緩，其尾線比右半部分曲線更長，無限延伸到接近x軸；2組數(shù)據(jù)均滿足平均數(shù)<中位數(shù)<眾數(shù)的關系，考慮2組數(shù)據(jù)滿足偏態(tài)分布。計算得第1組數(shù)據(jù)的峰度為0.92，偏度為-1.07；第2組數(shù)據(jù)的峰度1.52，偏度為-1.33。對于呈現(xiàn)偏態(tài)分布的數(shù)據(jù)，可利用Box-Cox變化減弱其偏態(tài)狀況，更接近于正態(tài)分布[20]，用x表示原始數(shù)據(jù)，y表示變化后得到的數(shù)據(jù)，其計算公式如下：

(1)

分別對這2組數(shù)據(jù)進行Box-Cox變化，通過Minitab軟件先確定參數(shù)λ的估計值如圖7所示。由圖7可知，2組數(shù)據(jù)中λ的最優(yōu)值均>5，最后得到2組數(shù)據(jù)的λ值分別為7.6，24。對變換后的數(shù)據(jù)通過SPSS分析軟件，利用K-S法進行正態(tài)性檢驗，得到P1=0.074>0.05，P2=0.052>0.05，說明2組變換后的數(shù)據(jù)均服從正態(tài)分布。

圖7 Box-Cox變換中的參數(shù)λ

3.2 缺陷率對軸壓承載力的影響

為研究缺陷率大小對結(jié)構(gòu)承載力的影響，應先控制缺陷的分布范圍，本文主要通過缺陷的球心坐標來確定缺陷的位置分布。針對2組參數(shù)下隨機生成的球心坐標劃分的范圍為：x方向-98.5～98.5mm，y方向-98.5～98.5mm，z方向500～700mm；即選取的缺陷分布范圍大致位于核心混凝土中心區(qū)域，具體如表3所示。

表3 中心區(qū)域隨機缺陷

根據(jù)表3中計算出各類缺陷的缺陷率及對應結(jié)構(gòu)的承載力折減系數(shù)，其中折減系數(shù)定義為帶缺陷結(jié)構(gòu)軸向承載力與無缺陷結(jié)構(gòu)軸向承載力的比值，經(jīng)計算不帶缺陷鋼管混凝土模型的承載力為4 781.07kN，得到缺陷率-折減系數(shù)折線如圖8所示。

圖8 缺陷率-折減系數(shù)折線

由圖8可知，隨著缺陷率的增大，折減系數(shù)整體呈下降趨勢，即構(gòu)件軸向承載力隨缺陷率的增大而減小。當缺陷率較小時，即在缺陷率≤0.005%的情況下，折線的斜率最大，平均斜率達到14.6，這表明在這一范圍內(nèi)結(jié)構(gòu)的軸向承載力下降最快；當缺陷率為0.005%～0.01%時，折線的平均斜率為8.2，當缺陷率>0.01%時，折線的平均斜率為5.1，這表明隨著缺陷率的增大，軸向承載力雖下降，但下降趨勢逐漸變緩。即使缺陷率很小的情況下，其軸向承載力也減小近40%，核心區(qū)混凝土內(nèi)部缺陷的存在對結(jié)構(gòu)承載力的影響很大，不容忽視。

3.3 缺陷位置對軸壓承載力的影響

為研究缺陷位置對結(jié)構(gòu)軸壓承載力的影響，先將結(jié)構(gòu)內(nèi)核心混凝土劃分成不同區(qū)域，再選取1個缺陷率范圍，該范圍內(nèi)的隨機缺陷在各區(qū)域內(nèi)均有分布，每個隨機缺陷均對應有結(jié)構(gòu)的軸向承載力，通過計算各區(qū)域軸向承載力的平均值并進行比較，從而找出缺陷位置對結(jié)構(gòu)軸向承載力的影響。區(qū)域劃分以缺陷球心坐標為依據(jù)，共有以下幾部分：①區(qū)域ax方向-100～100mm，y方向-100～100mm，z方向300～900mm，即混凝土的中心區(qū)域；②區(qū)域b 相較于區(qū)域a，z軸數(shù)值不變，x，y軸坐標位于區(qū)域a對應位置以外，即混凝土中間部分靠近兩側(cè)的區(qū)域；③區(qū)域cx方向-100～100mm，y方向-100～100mm，z方向0～300mm，900～1 200mm，即混凝土兩端靠近中間的區(qū)域；④區(qū)域dz方向坐標與區(qū)域c一致，x，y方向坐標位于區(qū)域c對應位置以外，即混凝土兩端靠近兩側(cè)的位置。各區(qū)域?qū)糠志唧w信息如表4～7所示。缺陷不同區(qū)域的結(jié)構(gòu)承載力直方圖如圖9所示。

由圖9可知，當缺陷位于c區(qū)域時，結(jié)構(gòu)承載力普遍較大，其承載力>3 050kN的頻率可達40%；當缺陷位于區(qū)域a時，結(jié)構(gòu)承載力主要集中在2 400～2 800kN，其頻率可達60%；當缺陷位于區(qū)域b時，結(jié)構(gòu)承載力主要集中在2 900～2 950kN,其頻率為40%；當缺陷位于區(qū)域d時，結(jié)構(gòu)承載力主要集中在2 700～2 900kN,其頻率為60%。通過計算，各區(qū)域?qū)S向承載力的平均值如下：區(qū)域a為2 662.86kN，區(qū)域b為2 863.66kN，區(qū)域c為3 001.98kN，區(qū)域d為2 873.54kN，通過比較可得c>d>b>a。當缺陷位于結(jié)構(gòu)兩端靠近中間部位時，缺陷對結(jié)構(gòu)軸向承載力的影響最??；當缺陷位于結(jié)構(gòu)兩端兩側(cè)位置和結(jié)構(gòu)中間兩側(cè)位置時，即當缺陷位于邊部位置時，無論其在結(jié)構(gòu)兩端還是結(jié)構(gòu)中間位置，其對結(jié)構(gòu)軸向承載力的影響基本相同；當缺陷位于結(jié)構(gòu)的中心部位時，結(jié)構(gòu)的軸向承載力平均值最低，對結(jié)構(gòu)軸向承載力的影響最大。

表4 區(qū)域a隨機缺陷

表5 區(qū)域b隨機缺陷

表6 區(qū)域c隨機缺陷

表7 區(qū)域d隨機缺陷

圖9 缺陷不同區(qū)域的結(jié)構(gòu)承載力直方圖

4 結(jié)語

1)通過對大量承載力數(shù)值進行統(tǒng)計分析可發(fā)現(xiàn)，帶缺陷鋼管混凝土柱承載力數(shù)值維持在一定的范圍內(nèi)，且服從偏態(tài)分布。

2)混凝土內(nèi)部缺陷的存在對鋼管混凝土柱軸向承載力的影響很大，在缺陷率很小的情況下，軸向承載力的折減就能達到40%以上；結(jié)構(gòu)承載力隨缺陷率的增大而降低，當缺陷率較小時，結(jié)構(gòu)承載力下降很快，隨著缺陷率的增大，結(jié)構(gòu)承載力下降開始變緩。

3)當缺陷位于核心混凝土中心部位時，對結(jié)構(gòu)的承載力影響最大；位于核心混凝土兩端靠近中間位置的缺陷對結(jié)構(gòu)承載力影響最小，當缺陷位于核心混凝土邊部位置時，對結(jié)構(gòu)承載力的影響位于上述2種情況之間。