金 瀏,祝華杰,杜修力
(北京工業(yè)大學城市減災與防災防護教育部重點實驗室,北京 100124)
鋼筋混凝柱作為建筑結(jié)構(gòu)的主要受力構(gòu)件,僅僅受到軸力作用的情況很少。在地震等水平荷載作用下,鋼筋混凝土柱可能處于壓-彎-剪-扭復合受力狀態(tài)下[1]。而扭轉(zhuǎn)作用的存在會加劇結(jié)構(gòu)的破壞,特別是一些處在特殊位置的柱,如框架結(jié)構(gòu)的邊柱、角柱和曲線梁橋的柱墩。因此,有必要探究鋼筋混凝土柱在扭轉(zhuǎn)作用及壓-彎-剪-扭復合受力狀態(tài)的力學性能。
實際中,扭轉(zhuǎn)作用通常與剪力、軸力或彎矩作用同時發(fā)生。但鋼筋混凝土構(gòu)件的純扭破壞性能研究是其包括扭轉(zhuǎn)作用的復合受力性能研究的基礎,有必要對其進行研究。在過去的幾十年里,Koutchoukali和Belarbi[2],F(xiàn)ang和Shiau[3],Chalioris[4],Chiu等[5],Bernardo和Lopes[6?7]對鋼筋混凝土結(jié)構(gòu)純扭轉(zhuǎn)的力學性能進行了一些試驗研究。探究了混凝土強度[2?3,5?6]、配箍率[5?6]、截面高寬比[4?5]、縱筋的布置[4]等因素對構(gòu)件純扭轉(zhuǎn)破壞的影響。研究發(fā)現(xiàn),縱筋率的提高對鋼筋混凝土構(gòu)件純扭破壞的延性行為和承載能力無較大改善[4];鋼筋混凝土構(gòu)件的純扭轉(zhuǎn)破壞具有一定的延性,但延性行為僅發(fā)生在縱筋與箍筋配筋強度比較窄的范圍內(nèi)[5?6]。Wang等[1]和Prakash等[8]對鋼筋混凝土方/圓柱進行了純扭作用的低周循環(huán)往復試驗,發(fā)現(xiàn)滯回曲線捏縮現(xiàn)象明顯,峰值扭矩后的承載力退化速度快。總的來說,鋼筋混凝土構(gòu)件的純扭轉(zhuǎn)破壞表現(xiàn)出較強的脆性行為。
實際上,鋼筋混凝土結(jié)構(gòu)(梁、柱、節(jié)點)的尺寸效應現(xiàn)象已被很多學者證實[9?13]。事實上,鋼筋混凝土構(gòu)件的扭轉(zhuǎn)破壞類似于剪切破壞。關于鋼筋混凝土結(jié)構(gòu)的剪切行為,大量學者進行了無/有腹筋梁的剪切試驗[9?10]與數(shù)值模擬研究[11?12],證明了抗剪強度存在尺寸效應,并且修正了名義抗剪強度計算公式,考慮了結(jié)構(gòu)尺寸的影響。目前,對于鋼筋混凝土構(gòu)件扭轉(zhuǎn)破壞的尺寸效應研究還較少。Ba?ant等[14]對最大截面尺寸為300 mm的素混凝土梁進行了純扭轉(zhuǎn)試驗,試驗結(jié)果表明素混凝土梁的扭轉(zhuǎn)破壞表現(xiàn)出較強的脆性行為,其名義扭轉(zhuǎn)強度存在明顯的尺寸效應。Kirane等[15]基于Ba?ant等[14]的試驗,進行了無/有腹筋梁(最大截面尺寸為250 mm)純扭轉(zhuǎn)破壞的數(shù)值模擬研究,發(fā)現(xiàn)其名義扭轉(zhuǎn)強度依然存在尺寸效應??梢园l(fā)現(xiàn),目前對鋼筋混凝土構(gòu)件純扭破壞尺寸效應的試驗研究大多集中在小尺寸素混凝土試件上,其尺寸效應規(guī)律認識不足,縱筋率、配箍率等因素對鋼筋混凝土純扭尺寸效應的影響有待深入研究。
綜上所述,目前對鋼筋混凝土構(gòu)件純扭及尺寸效應研究很少,對其破壞機理及尺寸效應規(guī)律認識不足。本文的目的是通過建立鋼筋混凝土圓/方柱純扭破壞三維細觀數(shù)值模型,探究縱筋率、配箍率、截面形狀對鋼筋混凝土柱抗扭強度尺寸效應的影響規(guī)律。最后,對比并修正了Ba?ant尺寸效應律,建立了全結(jié)構(gòu)尺寸范圍內(nèi)的名義抗扭強度尺寸效應預測公式。
鋼筋混凝土柱的尺寸效應主要來源于:① 混凝土材料的不均勻性;② 鋼筋與混凝土復雜的相互作用。對于第一個方面,需要建立一個能反映混凝土非均質(zhì)性的細觀尺度數(shù)值模型,其中混凝土被認為是由骨料、砂漿基體和界面過渡區(qū)組成的三相材料[16?17]。對于第二個方面,需要考慮鋼筋與混凝土之間復雜的非線性相互作用。
由此,建立了鋼筋混凝土柱三維細觀數(shù)值模型,如圖1所示。根據(jù)Cao等[18]的工作,經(jīng)過必要的簡化后,在數(shù)值模擬模型中增加了加載裝置。具體來說,柱底部是固定的,柱頂部的扭轉(zhuǎn)荷載通過施加轉(zhuǎn)角位移實現(xiàn)。
圖1 鋼筋混凝土柱三維細觀數(shù)值模型Fig. 1 Three-dimensional mesoscopic numerical model of reinforced concrete (RC) columns
混凝土各組分采用八節(jié)點六面體減縮積分單元進行劃分,鋼筋采用梁單元進行離散。采用Fortran程序,將粗骨料顆粒按照“取-放”的方法隨機投放在砂漿基質(zhì)中[19?20]。由于骨料形狀對混凝土的宏觀應力-應變行為影響不大[21],因此在三維數(shù)值模型中將骨料形狀簡化為球形。此外,骨料的分布由Fuller級配曲線確定[21?22],采用二級配,骨料的等效粒徑分別為30 mm和12 mm,骨料的體積分數(shù)為45%。界面過渡區(qū)是包裹在粗骨料周圍的薄層,其實際厚度在20 μm~50 μm[23?24]??紤]到模型計算的效率,該模型中界面過渡區(qū)厚度設為1 mm,與?avija等[25]的工作類似。首先建立起素混凝土柱三維細觀數(shù)值模型,參考文獻[9]工作,將有限元網(wǎng)格在三維結(jié)構(gòu)上進行投影,然后根據(jù)各組分在網(wǎng)格上的相對位置判定單元類型(如骨料、砂漿基質(zhì)、界面單元),并在此基礎上插入鋼筋,進而在鋼筋單元和混凝土單元間建立非線性彈簧單元。
骨料的強度明顯高于砂漿基質(zhì)和界面過渡區(qū)。在靜荷載作用下,一般認為試件失效時骨料不會破壞。因此,在細觀模型中,骨料被設置為彈性體[26]。砂漿基質(zhì)和界面過渡區(qū)采用塑性損傷本構(gòu)模型[27?28]來描述其力學性能。需要說明的是,為了緩解由于應變軟化引起的網(wǎng)格敏感性問題,在拉伸應力-應變曲線的下降部分,采用拉伸應力-位移曲線替代應力-應變曲線[29?30],從而保證了單元破壞時斷裂能的唯一性。
鋼筋通常被認為是均質(zhì)材料,其力學性能可以用理想彈塑性本構(gòu)模型來描述。此外,利用鋼筋與混凝土之間的非線性彈簧連接來定義鋼筋與混凝土之間的粘結(jié)滑移關系。采用GB 50010?2010[31]推薦的粘結(jié)滑移(τ-s)本構(gòu)模型來定義鋼筋與混凝土的相互作用,如圖2所示。表1給出了圖2中各參數(shù)的含義及確定方法。
圖2 粘結(jié)滑移本構(gòu)模型Fig. 2 The bond-slip curve
表1 粘結(jié)滑移本構(gòu)模型中采用的參數(shù)Table 1 Parameters utilized in the bond-slip model
參考Prakash等[8]的試驗,建立了鋼筋混凝土柱的三維細觀數(shù)值模型,驗證所建立模型的合理性和準確性。試驗分別對鋼筋混凝土圓柱和方柱進行了純扭載荷測試。同時,對試件頂部施加大小為7%的混凝土軸向承載力的軸向荷載。圓柱的直徑為600 mm,方柱的邊長為560 mm。圓柱和方柱的長細比為6,且具有相同的縱筋率和配箍率。更多的試驗信息詳見文獻[8]。
首先,需要確定混凝土各細觀組分和鋼筋的力學參數(shù)。其中力學參數(shù)包括抗壓強度、抗拉強度、彈性模量、泊松比、斷裂能、剪脹角和橫/縱向鋼筋屈服強度。表2給出了混凝土細觀組分和鋼筋的力學參數(shù)。其中“*”為試驗實測值(如鋼筋、砂漿基質(zhì)的力學參數(shù)根據(jù)試驗[8]取得),“^”為反復試算值,其他參數(shù)為默認值。界面過渡區(qū)的力學參數(shù)采用反復試算的方法得到[9,11?12]。建立標準圓柱體混凝土(?150 mm × 300 mm)三維細觀數(shù)值模型,通過對砂漿基質(zhì)的力學參數(shù)進行不同程度的折減(70%~85%)反復試算,當各組分的力學參數(shù)取表2所示參數(shù)時,混凝土單軸抗壓模擬值為34.4 MPa,試驗值為34.5 MPa,說明所取參數(shù)的合理性。此外,為了探究網(wǎng)格敏感性對模擬結(jié)果的影響,分別對網(wǎng)格尺寸為2 mm、3 mm和5 mm的鋼筋混凝土方柱進行了模擬,模擬與試驗結(jié)果對比如圖3所示。結(jié)果表明,不同網(wǎng)格尺寸下的扭矩-轉(zhuǎn)角關系差異較小。考慮到數(shù)值計算的效率,后續(xù)模擬選擇3 mm的網(wǎng)格尺寸。
表2 混凝土細觀組分及鋼筋力學參數(shù)Table 2 Mechanical parameters of the meso-components of concrete and reinforcing bars
圖3為試驗柱與模擬柱破壞模式的對比及扭矩-扭轉(zhuǎn)曲線的對比。可以發(fā)現(xiàn),斜裂縫的角度約為40°,與試驗結(jié)果相似;峰值扭矩的模擬值與試驗值的差異非常小。需要說明的是,在試驗中試件受循環(huán)荷載,而在模擬中試件受單調(diào)荷載。由于損傷累積和疲勞響應,試驗中扭矩-轉(zhuǎn)角曲線下降部分的軟化速率比模擬中要快??傮w而言,該模型的模擬結(jié)果與實驗結(jié)果吻合較好,說明了三維細觀數(shù)值方法的準確性和適用性。
圖3 模擬結(jié)果與試驗結(jié)果的對比Fig. 3 Comparison of simulation results and test results
基于上述細觀數(shù)值模擬方法,建立了64根鋼筋混凝土柱在純扭轉(zhuǎn)荷載作用下的數(shù)值模型。一共有32個圓柱和32個方柱,具有不同的結(jié)構(gòu)尺寸(D(B) ×h = 200 mm×600 mm、400 mm×1200 mm、800 mm×2400 mm和1000 mm×3000 mm),不同的縱筋率(ρst= 0%、0.64%、1.28%和1.92%)和不同配箍率(ρsv= 0.283%、0.377%、0.565%和1.131%)。通過改變縱筋數(shù)量來改變縱筋率,通過改變箍筋間距來改變配箍率。綜上所述,探討了以下3個因素對鋼筋混凝土柱純扭轉(zhuǎn)破壞尺寸效應的影響:① 截面形狀;② 縱筋率;③ 配箍率。
模擬的試件可分為素混土柱、無腹筋混凝土柱、有腹筋混凝土柱,下面依次分析其破壞模式和鋼筋屈服情況。圖4展示了素混凝土柱純扭破壞過程。模擬結(jié)果表明,在扭矩的作用下,試件先在混凝土柱表面中點附件被拉裂,并迅速延伸形成螺旋斜裂縫,最終試件被劈裂成兩半。破壞模式與Ba?ant等[14]的試驗一致。
圖4 素混凝土柱破壞過程Fig. 4 Failure process of plain concrete columns
圖5展示了無腹筋柱與素混凝土柱的破壞模式的對比,可以發(fā)現(xiàn),加入縱筋后,試件破壞時的扭轉(zhuǎn)斜裂縫增多,但繼續(xù)提高縱筋率對破壞模式無明顯改善。從圖6的無腹筋混凝土柱破壞過程可以看出,縱筋直到試件破壞時也未屈服。這是因為沒有橫向鋼筋約束的縱筋只能通過銷栓作用抗扭,銷栓作用比較小且不可靠,縱筋的抗扭作用未得到充分發(fā)揮。
圖5 無腹筋混凝土柱破壞模式Fig. 5 Failure mode of the RC columns without stirrups
圖6 無腹筋混凝土柱破壞過程Fig. 6 Failure process of RC columns without stirrups
圖7展示了有腹筋柱純扭破壞過程。與無腹筋柱相比,因為有腹筋的約束,縱筋和腹筋在扭矩達到峰值前都可以屈服,充分發(fā)揮了鋼筋的抗扭作用,抗扭能力較無腹筋構(gòu)件顯著提升。
圖7 有腹筋混凝土柱純扭破壞過程Fig. 7 Failure process of RC columns with stirrups
圖8展示了構(gòu)件在相同結(jié)構(gòu)尺寸下不同配箍率下的破壞模式圖。隨著配箍率的增加,構(gòu)件破壞時的扭轉(zhuǎn)斜裂縫增多。并且,斜裂縫與水平面的夾角變大,這是因為隨著配箍率的提高,其縱筋與箍筋的配筋強度比改變,斜裂縫與水平面的傾角改變。
圖8 有腹筋混凝土柱破壞模式Fig. 8 Failure mode of the RC columns with stirrups
圖9為無腹筋柱的扭矩-轉(zhuǎn)角曲線。對于素混土,扭矩-轉(zhuǎn)角曲線的下降段是幾乎垂直的,為典型的脆性破壞。加入縱筋后,其扭矩在峰值扭矩后并沒有突然下降,在一定程度上改善了素混凝土的脆性。并且小尺寸柱的曲線下降段較大尺寸柱更為平緩,說明縱向鋼筋對小尺寸構(gòu)件的加固效果優(yōu)于大尺寸柱。但是,縱向鋼筋對抗扭能力并沒有太大的影響。圖10為有腹筋柱的扭矩-轉(zhuǎn)角曲線。對于圓柱,在同一結(jié)構(gòu)尺寸下,隨著腹筋率的提高峰值扭矩顯著提升,T-θ曲線的下降段變緩。而對于方柱,在同一結(jié)構(gòu)尺寸下,不同腹筋率的峰值扭矩有所改變,但T-θ曲線的下降段趨勢改變較小。對比圓柱和方柱的T-θ曲線的下降段,方柱的下降段更陡,這是因為方形箍筋對混凝土的約束不均勻,相同腹筋率下,箍筋對混凝土的約束作用圓柱優(yōu)于方柱。
圖9 無腹筋混凝土柱的扭矩-轉(zhuǎn)角曲線Fig. 9 Torque-twist curves of the RC columns without stirrups
圖10 有腹筋混凝土柱的扭矩-轉(zhuǎn)角曲線Fig. 10 Torque-twist curves of the RC columns with stirrups
圖11給出了名義扭轉(zhuǎn)強度與結(jié)構(gòu)尺寸的關系。名義抗轉(zhuǎn)強度的定義為[15]:式中:νu為名義扭轉(zhuǎn)強度;Tmax為峰值扭矩;Wtp為受扭塑形抵抗矩,對于圓形截面,Wtp=2/3×πR3,R為圓形截面的半徑,對于矩形截面,Wtp= 1/6×b2(3h?b),b和h分別為矩形的短邊和長邊。顯然,隨著結(jié)構(gòu)尺寸的增加,鋼筋混凝土柱的名義扭轉(zhuǎn)強度降低。從圖11的趨勢線看出:① 方柱的名義扭轉(zhuǎn)強度下降趨勢更快,說明方柱的尺寸效應強于圓柱;② 無腹筋柱的名義扭轉(zhuǎn)強度下降更快,說明無腹筋柱的尺寸效應強于有腹筋柱。當然,僅從名義扭轉(zhuǎn)強度的下降趨勢中無法準確地得出此規(guī)律,需要一個理論的尺寸效應律公式,下面將詳細說明。
圖11 名義扭轉(zhuǎn)強度與結(jié)構(gòu)尺寸的關系Fig. 11 Relationship between nominal torsional strength and structural size
圖12給出了軟化率的計算公式。軟化率為ν-θ曲線中名義抗扭強度峰值νu到0.85νu斜率的絕對值[9],軟化率可以直觀地反應鋼筋混凝土柱扭轉(zhuǎn)破壞的脆性行為。
圖12 軟化率Fig. 12 Definition of softening rate
圖13為有腹筋混凝土柱軟化率與結(jié)構(gòu)尺寸的關系。可以發(fā)現(xiàn),大尺寸柱的軟化率大于小尺寸柱,即大尺寸柱的扭轉(zhuǎn)破壞更具脆性,在一定程度上導致了名義扭轉(zhuǎn)強度的尺寸效應。對比圓柱和方柱的軟化率,可以發(fā)現(xiàn),方柱的軟化率大于圓柱,表明方柱的脆性行為大于圓柱。這是因為圓形箍筋對混凝土的約束作用優(yōu)于方形箍筋,圓柱在扭轉(zhuǎn)破壞時具有更好的延性。并且,隨著配箍率的提高,軟化率降低,試件的延性得到提升,名義扭轉(zhuǎn)強度的尺寸效應得到削弱。這里需要說明的是,配箍率為1.131%的方柱為超筋柱,破壞模式屬于超筋破壞,導致軟化率大于其他試件。
圖13 有腹筋柱的軟化率Fig. 13 Softening rate k of RC columns with stirrups
目前,國內(nèi)外學者已提出了許多混凝土尺寸效應理論,如統(tǒng)計尺寸效應理論[32]、斷裂力學尺寸效應理論[33]、邊界尺寸效應理論[34]。其中,Ba?ant[33]基于斷裂力學提出的混凝土材料層次的尺寸效應理論被較為廣泛地接受,其尺寸效應理論分為2種:即Type-1型尺寸效應和Type-2型尺寸效應。
Type-1型尺寸效應出現(xiàn)在宏觀裂縫剛剛生成即破壞的混凝土結(jié)構(gòu)中。素混凝土柱的扭轉(zhuǎn)破壞是脆性的,破壞發(fā)生在混凝土表面的裂紋萌生處。因此,基于斷裂能的尺寸效應僅僅是由應變局域化造成的。這種情況下,大尺寸構(gòu)件的名義強度在雙對數(shù)圖中將無限趨近于一條不為零的水平漸近線[15]。Type-1型尺寸效應公式為[15]:
式中:νu∞為無限大尺寸(D→∞)的名義抗扭強度;r和Db為材料參數(shù)。由于扭轉(zhuǎn)數(shù)據(jù)較少,參考Kirane等[15]工作,暫定r = 1.5。通過對素混凝土柱模擬數(shù)據(jù)的擬合,得到圓柱的νu∞= 3.13 MPa,Db= 62.52 mm和方柱的νu∞= 2.91 MPa,Db=93.21 mm。
Type-2型尺寸效應出現(xiàn)在裂縫或缺口顯著大于斷裂過程區(qū)尺寸的混凝土結(jié)構(gòu)中。Type-2型尺寸效應公式為[33]:
式中:ν0和D0為材料參數(shù)。擬合得到圓柱的ν0=4.19 MPa,D0= 1525.15 mm和方柱的ν0= 4.26 MPa,D0= 1205.94 mm。
圖14為模擬數(shù)據(jù)與Type-1型和Type-2型尺寸效應律的對比。從圖14可以看出,在模擬數(shù)據(jù)的尺寸范圍內(nèi),數(shù)據(jù)點與Type-1型和Type-2型尺寸效應曲線都較為吻合。但在全尺寸范圍內(nèi),兩條曲線存在以下兩個問題:① 扭轉(zhuǎn)破壞的尺寸效應理論屬于Type-1型尺寸效應,大尺寸構(gòu)件的扭轉(zhuǎn)強度應趨于不為零的水平漸近線[15],Type-2型尺寸效應曲線不能滿足;② 當結(jié)構(gòu)尺寸無限小(D→0)時,Type-1型尺寸效應曲線中的扭轉(zhuǎn)強度趨于無窮大,這與實際不符。因此,需要對尺寸效應公式進行修正。
圖14 尺寸效應律Fig. 14 Size effect law
基于以上的分析,對Ba?ant尺寸效應律公式進行了修正,修正公式為:
式中:ν0和D0為材料參數(shù);νu∞為無限大尺寸(D→∞)的名義抗扭強度。Kim[35]對Ba?ant尺寸效應公式進行了類似的修正,提出了新的抗剪強度尺寸效應預測公式,定義νu∞= 0.47ft。實際上,νu∞的確定應該采用統(tǒng)計的方法,因為缺少扭轉(zhuǎn)數(shù)據(jù),暫且定義νu∞= 0.5ft。通過對素混凝土柱模擬數(shù)據(jù)的擬合,得到圓柱的ν0= 4.58 MPa,D0=245.21 mm和方柱的ν0= 4.80 MPa,D0= 157.11 mm。修正后的尺寸效應公式如圖14所示,對全尺寸的抗扭強度有了更好的預測。
圖15為模擬數(shù)據(jù)與新尺寸效應律公式的擬合結(jié)果。與Ba?ant尺寸效應律公式類似,雙對數(shù)圖中的水平線代表彈性或彈塑性理論(強度準則),表明材料無尺寸效應;斜率為?1/2的斜線代表線彈性斷裂力學理論,表明材料有強尺寸效應。對比雙對數(shù)圖中的數(shù)據(jù)點,發(fā)現(xiàn)無腹筋柱的數(shù)據(jù)更靠近?1/2斜線,表明無腹筋柱的尺寸效應更強,加入腹筋可有效削弱尺寸效應;不同縱筋率的無腹筋柱的數(shù)據(jù)點分布比較集中,表明加入縱筋對尺寸效應無顯著影響;大腹筋率試件的數(shù)據(jù)點更靠近水平線,說明提高腹筋率可削弱尺寸效應;方柱的數(shù)據(jù)點較圓柱更接近?1/2斜線,表明方柱的尺寸效應強于圓柱。
圖15 名義抗扭強度雙對數(shù)圖Fig. 15 Double logarithmic curve of nominal torsional strength
Ba?ant等[14]和Kirane等[15]進行了RC柱純扭轉(zhuǎn)尺寸效應試驗和模擬工作,試驗和模擬數(shù)據(jù)如表3所示。將試驗和模擬數(shù)據(jù)與新尺寸效應律公式進行對比,如圖16所示??梢园l(fā)現(xiàn),新建立的尺寸效應公式與試驗數(shù)據(jù)吻合良好,驗證了所建立公式的準確性。但由于試驗數(shù)據(jù)數(shù)量較少,所建公式的合理性與準確性還需更多相關試驗作進一步驗證。
表3 純扭試驗數(shù)據(jù)庫Table 3 Database of pure torsion tests
圖16 尺寸效應律公式的驗證Fig. 16 Verification of size effect law formula
本文通過三維細觀數(shù)值模擬方法,探究了鋼筋混凝土柱的純扭破壞機理及尺寸效應。分析了縱筋率、箍筋率、截面形狀對承載能力、變形能力和尺寸效應的影響。主要結(jié)論為:
(1) 鋼筋混凝土柱的名義抗扭強度存在尺寸效應,由于箍筋對方柱的約束作用弱于圓柱,導致方柱的尺寸效應強于圓柱。
(2) 名義抗扭強度的尺寸效應會隨著箍筋率的增加而減弱,且圓柱的削弱作用強于方柱,而縱向配筋對抗扭能力和尺寸效應的影響可以忽略。
(3) 修正后的尺寸效應律能夠較好地描述全尺寸范圍內(nèi)的鋼筋混凝土柱抗扭強度尺寸效應現(xiàn)象。
本研究進行了鋼筋混凝土柱純扭荷載的數(shù)值模擬,是進行鋼筋混凝土柱包括扭轉(zhuǎn)作用復合受力研究的基礎。關于鋼筋混凝土柱壓-扭、彎-扭、壓-彎-剪-扭等組合受力性能的研究將在后續(xù)工作中展開,探究軸壓比、剪跨比、扭彎比等對構(gòu)件力學性能和尺寸效應的影響。