黃子睿， 陳忠范， 黃東升， 何敏娟

(1 東南大學(xué)土木工程學(xué)院， 南京 211189；2 南京林業(yè)大學(xué)生物質(zhì)材料國(guó)家地方聯(lián)合工程研究中心， 南京 210037；3 阿伯特大學(xué)土木與環(huán)境工程系， 埃德蒙頓 T6G1H9； 4 同濟(jì)大學(xué)土木工程學(xué)院， 上海 200092)

0 引言

工程竹是原竹經(jīng)工業(yè)化制造而成的竹基復(fù)合材料。根據(jù)膠合單元不同，工程竹通常有重組竹、集成竹、竹膠板等。在工程竹制造過(guò)程中，對(duì)原竹材進(jìn)行了篩選，僅截取竹莖中材質(zhì)比較均勻的部分，剔除其中的缺陷，并將不同部位進(jìn)行分散搭配，經(jīng)高溫干燥處理后，再進(jìn)行膠合。因此，工程竹不僅顯著減少了原竹幾何尺寸和力學(xué)性能的變異性，且強(qiáng)度較高，收縮和翹曲小，適用于制作建筑梁、柱和樓板等基本構(gòu)件，是一種理想的建筑結(jié)構(gòu)材料[1]。

近十年來(lái)，我國(guó)學(xué)者針對(duì)工程竹結(jié)構(gòu)的基礎(chǔ)和應(yīng)用研究開(kāi)展了許多工作[2-5]，發(fā)展了輕型工程竹結(jié)構(gòu)、工程竹框架結(jié)構(gòu)、鋼-竹組合結(jié)構(gòu)等多種結(jié)構(gòu)形式，取得了一系列成果。工程竹作為一種可再生、綠色環(huán)保材料，也已逐漸被人們接受，我國(guó)各地已陸續(xù)建造了許多現(xiàn)代竹結(jié)構(gòu)建筑[6]。由于目前還沒(méi)有工程竹結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)，現(xiàn)代竹結(jié)構(gòu)建筑設(shè)計(jì)與建造主要參考現(xiàn)行《木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)》(GB 50005—2017)[7](簡(jiǎn)稱(chēng)木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn))。然而，研究證實(shí)，工程竹材料力學(xué)性質(zhì)與木材力學(xué)性質(zhì)并不完全相同。通常，木材的順紋受壓比例極限與其順紋受壓強(qiáng)度差異并不太大，有的文獻(xiàn)甚至以理想彈塑性關(guān)系描述木材的順紋受壓應(yīng)力-應(yīng)變關(guān)系[8]。因此，采用線彈性理論計(jì)算木構(gòu)件的極限承載力與其實(shí)際承載力比較接近，且偏于安全，在工程上可以接受。然而，竹材的受壓應(yīng)力-應(yīng)變曲線則表現(xiàn)出顯著的非彈性特點(diǎn)，其受壓比例極限約為其極限受壓強(qiáng)度的50%[1]。工程竹的這種非彈性特點(diǎn)，對(duì)受彎構(gòu)件或壓彎構(gòu)件的力學(xué)行為影響顯著[9-14]。準(zhǔn)確預(yù)估工程竹結(jié)構(gòu)的極限承載力，應(yīng)該考慮材料的非彈性性質(zhì)。若采用木結(jié)構(gòu)設(shè)計(jì)公式[15-17]，將會(huì)顯著低估工程竹構(gòu)件的極限承載力。

近年來(lái)，筆者所在課題組以重組竹和集成竹為研究對(duì)象，系統(tǒng)開(kāi)展了材料基本力學(xué)性能、受彎構(gòu)件和壓彎構(gòu)件的極限承載力非彈性分析方法研究，建立了此類(lèi)構(gòu)件極限承載力與變形計(jì)算方法[1,5,9-13]。本文系統(tǒng)總結(jié)了這些成果，給出了相應(yīng)的承載力與變形計(jì)算主要公式，以期為工程竹結(jié)構(gòu)極限狀態(tài)設(shè)計(jì)提供依據(jù)。

1 非彈性應(yīng)力-應(yīng)變關(guān)系

重組竹一般選用毛竹為原料，首先將竹莖切成長(zhǎng)2m左右的竹條，在80℃左右環(huán)境中烘干至含水率11%～16%，再將竹條碾成竹絲束，平行組坯，浸漬酚醛樹(shù)脂膠，在一定壓力和溫度下膠合而成。集成竹是將長(zhǎng)約3m、寬30mm、厚20mm的竹條在80℃環(huán)境中烘干至含水率6%～11%，再將酚醛樹(shù)脂均勻涂抹于竹片表面，同方向組坯膠合而成。

文獻(xiàn)[1]和文獻(xiàn)[11]的研究表明，不論是重組竹還是集成竹，其順紋應(yīng)力-應(yīng)變關(guān)系均表現(xiàn)出類(lèi)似的特性，即順紋受拉呈完全彈性，而順紋受壓在彈性極限內(nèi)呈線性，超過(guò)彈性極限后呈明顯非線性特性。為了建立工程竹受彎構(gòu)件和壓彎構(gòu)件的非線性分析模型，文獻(xiàn)[1]采用如圖1所示二次曲線模擬工程竹的順紋本構(gòu)關(guān)系：

圖1 工程竹順紋應(yīng)力-應(yīng)變關(guān)系

(1)

其中：

式中：λ1，λ2，λ3均為材料常數(shù)；fcu和εcu分別為極限抗壓強(qiáng)度和極限壓應(yīng)變；fce和εce分別為比例極限壓應(yīng)力和比例極限壓應(yīng)變；ftu和εtu分別為極限抗拉強(qiáng)度和極限拉應(yīng)變；E為順紋楊氏模量。

2 受彎構(gòu)件極限承載力

木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)給出的受彎構(gòu)件承載力計(jì)算公式為M/Wn≤fm，其中，M，Wn，fm分別為受彎構(gòu)件控制截面的彎矩、截面凈抵抗矩以及材料抗彎強(qiáng)度設(shè)計(jì)值?？梢?jiàn)，木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)中極限承載力的計(jì)算采用了線彈性力學(xué)模型。由于工程竹受壓本構(gòu)關(guān)系非線性特性明顯，故利用線彈性模型計(jì)算工程竹極限承載力會(huì)導(dǎo)致較大誤差。試驗(yàn)[9]表明，工程竹受彎構(gòu)件破壞時(shí)，破壞截面受壓區(qū)靠近構(gòu)件上表面一定范圍內(nèi)的材料處于受壓屈曲狀態(tài)，使構(gòu)件產(chǎn)生非彈性響應(yīng)；而截面的其他區(qū)域，均處于彈性工作狀態(tài)，如圖2所示。因此，工程竹受彎構(gòu)件的極限彎矩Mu可由下式得到：

(2)

式中：f(y)為受壓區(qū)截面應(yīng)力分布；yc，yp，yt分別為彈性受壓區(qū)(ECZ)、非彈性受壓區(qū)(PCZ)和彈性受拉區(qū)(ETZ)的高度(圖2)，可通過(guò)截面平衡條件和幾何條件求得；kp為構(gòu)件曲率。

圖2 極限狀態(tài)下受彎構(gòu)件截面的應(yīng)力與應(yīng)變分布

(3)

計(jì)算結(jié)果[9-11]證明，采用式(2)計(jì)算重組竹與集成竹受彎構(gòu)件的承載力具有較高的精度。然而，式(3)是一個(gè)與材料極限應(yīng)變相關(guān)的參數(shù)，不便于工程應(yīng)用。為此，文獻(xiàn)[13]采用梯形分布代替受壓屈曲區(qū)的非線性應(yīng)力分布，得到了工程竹受彎構(gòu)件的極限彎矩計(jì)算公式：

(4)

式中b，h分別為梁截面寬度和高度。

以課題組研發(fā)的空心竹樓蓋單向板為例(圖3)，經(jīng)計(jì)算表明，式(4)的計(jì)算結(jié)果與試驗(yàn)結(jié)果誤差不大，且偏于安全，可以滿足工程計(jì)算精度要求。

圖3 空心竹樓蓋單向板試驗(yàn)與計(jì)算值比較[11]

3 壓彎構(gòu)件

3.1 單偏心受壓構(gòu)件

各國(guó)主要規(guī)范，如NDS-2015[14]，CSA086-14[15]和Eurocode 5[16]，給出了以軸向壓力-彎矩方程表達(dá)的構(gòu)件失效準(zhǔn)則，以校核壓彎構(gòu)件是否滿足極限狀態(tài)下的承載力要求。盡管有些規(guī)范或標(biāo)準(zhǔn)給出的失效準(zhǔn)則考慮了構(gòu)件的非線性特征，但極限應(yīng)力仍是基于線彈性理論給出的結(jié)果。如木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)給出的單向偏心受壓柱極限承載力驗(yàn)算公式：

(5)

式中：N和M0分別為偏心受壓構(gòu)件承受的軸向壓力和彎矩；e0為軸向壓力的初始偏心距；An和Wn分別截面凈面積和凈抵抗矩；fc和fm分別為材料抗壓強(qiáng)度與抗彎強(qiáng)度設(shè)計(jì)值。

與木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)類(lèi)似的還有歐洲規(guī)范Eurocode 5[16]，該規(guī)范給出的單向偏心受壓柱承載力驗(yàn)算公式與式(5)類(lèi)似，只不過(guò)采用應(yīng)力為基本變量。顯然，式(5)依據(jù)線彈性理論，將軸向壓力和彎矩產(chǎn)生的荷載效應(yīng)進(jìn)行線性疊加，未考慮材料的非彈性響應(yīng)。有些設(shè)計(jì)標(biāo)準(zhǔn)考慮了壓彎構(gòu)件的非線性響應(yīng)，將壓彎構(gòu)件的破壞準(zhǔn)則表達(dá)為軸向壓力與彎矩的二次曲線。如加拿大木結(jié)構(gòu)設(shè)計(jì)手冊(cè)[18]給出的在偏心壓力作用下中長(zhǎng)柱的極限承載力計(jì)算公式：

(6)

式中：Pf為軸向荷載；Pr為軸向抗力極限值；Mr為抵抗彎矩極限值；PE為Euler臨界荷載。

事實(shí)上，受壓構(gòu)件的承載力與其長(zhǎng)細(xì)比密切相關(guān)。對(duì)于細(xì)長(zhǎng)柱，整體失穩(wěn)是其主要破壞模式，早在1744年，Euler就提出了細(xì)長(zhǎng)柱的極限承載力計(jì)算公式[19]。當(dāng)柱的長(zhǎng)細(xì)比介于短柱和細(xì)長(zhǎng)柱之間(即中長(zhǎng)柱)時(shí)，其極限承載力取決于材料抗壓強(qiáng)度和長(zhǎng)細(xì)比，柱破壞時(shí)，非彈性特征顯著。

為了計(jì)算承受軸向壓力中長(zhǎng)柱的極限承載力，Engesser于1889年提出了切線模量理論[20]，該理論以切線模量Et，即應(yīng)力-應(yīng)變曲線某點(diǎn)的斜率，代替彈性模量E，得到的柱承載力計(jì)算公式和Euler公式具有相同的形式。此外，Von Karman于1919年提出割線模量理論[21]，該理論考慮材料卸載彈性應(yīng)力-應(yīng)變關(guān)系，得到的承載力公式也與Euler公式具有相同的形式，只不過(guò)以割線模量Er代替Euler公式中的彈性模量E。試驗(yàn)結(jié)果[10]表明，Engesser和von Karman公式能較好預(yù)估軸心受壓工程竹中長(zhǎng)柱的極限承載力，但不能用于計(jì)算偏心受壓柱的極限承載力。

文獻(xiàn)[21]在Euler和Timoshenko梁理論基礎(chǔ)上，建立了彎曲構(gòu)件和壓彎構(gòu)件統(tǒng)一模型，即梁-柱模型，研究了桿件偏心受壓或在軸向壓力與端部彎矩的聯(lián)合作用下，構(gòu)件極限承載力計(jì)算問(wèn)題，給出的壓彎構(gòu)件控制方程為彎矩關(guān)于曲率的二階微分方程，材料的非線性則通過(guò)破壞面的彎矩-曲率-軸力關(guān)系表示。對(duì)端部同時(shí)承受軸力和彎矩聯(lián)合作用的壓彎構(gòu)件，當(dāng)軸力恒定時(shí)，采用梁-柱模型很容易得到微分方程的解。對(duì)于偏心受壓柱，由于計(jì)算截面彎矩是軸力的函數(shù)，難以得到方程解析解，需采用數(shù)值方法對(duì)一系列荷載增量進(jìn)行逐步計(jì)算，以逼近極限承載力真實(shí)值。采用這種方法進(jìn)行結(jié)構(gòu)設(shè)計(jì)計(jì)算非常繁瑣。文獻(xiàn)[22]提出一種基于應(yīng)變的非彈性分析方法，該方法假設(shè)破壞截面上的應(yīng)變已知，根據(jù)應(yīng)力-應(yīng)變關(guān)系得出破壞截面的應(yīng)力分布，從而計(jì)算破壞截面的撓度和內(nèi)力。文獻(xiàn)[10]根據(jù)中長(zhǎng)柱偏心受壓試驗(yàn)結(jié)果，考慮工程竹拉、壓應(yīng)力-應(yīng)變關(guān)系的差異和極限狀態(tài)下受壓應(yīng)力-應(yīng)變關(guān)系非彈性特點(diǎn)，將破壞截面分為彈性區(qū)域和受壓屈曲區(qū)域(圖4)，得到中長(zhǎng)壓彎構(gòu)件彎曲破壞時(shí)極限承載力計(jì)算公式：

圖4 單偏壓柱及其破壞截面應(yīng)力分布

(7)

為了計(jì)算柱破壞截面的撓度，文獻(xiàn)[10]還提出了一種虛擬塑性鉸法。該方法假定破壞面的撓度由彈性撓度和非線性變形撓度兩部分組成，彈性撓度是指荷載達(dá)到比例極限時(shí)產(chǎn)生的撓度，由Euler方程求得；非線性變形撓度是指外荷載由比例極限增加至極限荷載過(guò)程中產(chǎn)生的撓度，在此過(guò)程中，破壞截面附近一定范圍內(nèi)由于非彈性壓曲，沿柱縱向產(chǎn)生一定長(zhǎng)度的塑性鉸。通過(guò)塑性鉸發(fā)展過(guò)程中的幾何關(guān)系分析，可得偏心受壓柱破壞時(shí)破壞截面(1/2柱高)處的非線性變形撓度δp：

(8)

式中：Lp為塑性鉸長(zhǎng)度；ytu為截面破壞時(shí)受拉區(qū)高度。

可見(jiàn)，非線性變形撓度取決于塑性鉸的長(zhǎng)度。假設(shè)塑性鉸長(zhǎng)度等于破壞截面高度，即Lp=h，則非線性撓度為：

(9)

3.2 雙向偏心受壓構(gòu)件

雙向偏心受壓構(gòu)件的非彈性分析模型研究較多[23-27]，大多數(shù)研究都假設(shè)材料具有理想彈塑性應(yīng)力-應(yīng)變關(guān)系以減少計(jì)算難度。文獻(xiàn)[17]采用了與單向偏壓構(gòu)件類(lèi)似的分析方法，建立了工程竹中長(zhǎng)柱在雙向偏心受壓作用下的極限承載力分析模型。該模型假定柱破壞時(shí)，最外層纖維拉斷，受壓區(qū)纖維屈曲。根據(jù)平截面假定，將柱的破壞截面分為彈性受拉區(qū)、彈性受壓區(qū)和非彈性受壓區(qū)，由幾何條件得到破壞截面應(yīng)變分布方程，從而進(jìn)一步得到極限承載力計(jì)算公式。由平截面假定可知，應(yīng)變分布方程可表示為：

ε(x,y)=Px+Qy+R

(10)

式中P，Q，R為應(yīng)變分布平面的外法線向量。

根據(jù)文獻(xiàn)[17]的研究結(jié)果，工程竹雙向偏心受壓構(gòu)件破壞時(shí)，破壞截面上的非彈性壓曲區(qū)有兩種可能的幾何形狀，即三角形或梯形。以非彈性受壓區(qū)為三角形為例，圖5給出了破壞截面的應(yīng)變和應(yīng)力分布。

圖5 非彈性受壓區(qū)為三角形時(shí)的截面應(yīng)變和應(yīng)力分布

由平衡條件可得上述兩種情況下雙向偏心受壓構(gòu)件的極限承載力計(jì)算公式：

非彈性受壓區(qū)為三角形時(shí)：

(11)

非彈性受壓區(qū)為梯形時(shí)：

3(b1+b2)2-(b1-b2)2]Q-3(2b-b1-

b2)R}-αa(2b-b1-b2)f

(12)

圖6 雙向偏心受壓柱截面上的彈性和非彈性區(qū)域

由于求解高階方程的困難，通過(guò)平衡條件和幾何條件確定彈性受壓區(qū)和塑性受壓區(qū)的分界線方程在數(shù)學(xué)上很難實(shí)現(xiàn)，故難以給出應(yīng)變平面的外法線向量顯式解析解。試驗(yàn)和分析證實(shí)，雙向偏心受壓柱截面的破壞始于塑性受壓區(qū)最外邊緣纖維的屈服，彈性和非彈性區(qū)域分界線為直線且與中性軸平行，它始于塑性受壓區(qū)的最遠(yuǎn)端，在截面發(fā)生損傷與破壞過(guò)程中，逐步平移至截面內(nèi)部。因此，式(11)和式(12)提供了一個(gè)計(jì)算雙向偏心受壓構(gòu)件非線性響應(yīng)的漸近法，即通過(guò)應(yīng)變分布平面方程，得到彈性和非彈性區(qū)域的分界線方程，逐漸將分界線從受壓區(qū)最外點(diǎn)向截面內(nèi)平移，計(jì)算相應(yīng)的外荷載和受拉區(qū)最外端的應(yīng)力，直至受拉區(qū)達(dá)到抗拉極限，停止計(jì)算，此時(shí)得到的外荷載值即為雙向偏心受壓柱的極限承載力。計(jì)算表明，上述方法可以得到滿意的結(jié)果。計(jì)算與試驗(yàn)結(jié)果對(duì)比如表1所示。

雙向偏心受壓柱極限承載力試驗(yàn)值與計(jì)算值對(duì)比[17] 表1

4 結(jié)語(yǔ)

工程竹順紋受拉和受壓應(yīng)力-應(yīng)變差異顯著，受壓應(yīng)力-應(yīng)變關(guān)系呈明顯非線性特征，因此，工程竹極限承載力與變形計(jì)算應(yīng)考慮其非彈性響應(yīng)。本文系統(tǒng)總結(jié)了筆者所在課題組關(guān)于工程竹受彎構(gòu)件和偏心受壓構(gòu)件極限承載力與變形的計(jì)算方法，列出了相關(guān)計(jì)算公式。其中，受彎構(gòu)件和單向偏心受壓構(gòu)件可通過(guò)本文公式直接計(jì)算相應(yīng)的極限承載力，雙向偏心受壓構(gòu)件則需要通過(guò)逐步漸進(jìn)法進(jìn)行計(jì)算。就作為實(shí)用計(jì)算公式而言，式(2)，(3)，(7)，(11)，(12)還需要進(jìn)一步簡(jiǎn)化，以下工作還有待進(jìn)一步開(kāi)展：

(1)受壓屈曲區(qū)應(yīng)力不均勻系數(shù)α的取值問(wèn)題。研究[9,11,17]表明，不論彎曲問(wèn)題還是偏心受壓?jiǎn)栴}，α值的計(jì)算公式完全相同，通過(guò)計(jì)算分析，給定α某一確定的近似值，則式(2)，(3)，(7)，(11)，(12)將大為簡(jiǎn)化，從而為得到工程竹極限承載力實(shí)用計(jì)算公式創(chuàng)造條件。

(2)各國(guó)木結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)(規(guī)范)只給出了偏心受壓構(gòu)件的承載力校核公式，沒(méi)有給出極限承載力計(jì)算公式。本文所列公式均可計(jì)算偏心受壓構(gòu)件的極限承載力。因此，通過(guò)進(jìn)一步分析、簡(jiǎn)化，給出偏心受壓構(gòu)件承載力實(shí)用計(jì)算公式，對(duì)編制工程竹結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)具有理論意義。

(3)本文討論的偏心受壓構(gòu)件承載力計(jì)算公式均在受拉區(qū)纖維應(yīng)力達(dá)到抗拉強(qiáng)度的假定下導(dǎo)出，即假設(shè)偏壓構(gòu)件的破壞模式為大偏心受壓破壞。當(dāng)偏心距較小時(shí)，受拉區(qū)纖維在構(gòu)件破壞時(shí)達(dá)不到抗拉極限強(qiáng)度。因此，對(duì)于小偏心受壓構(gòu)件的非彈性分析，還需進(jìn)一步研究。