李月強(qiáng),衣 娜,席 豐

(1.山東建筑大學(xué)土木工程學(xué)院,山東 濟(jì)南 250101; 2.山東軍之星建筑設(shè)計(jì)有限公司,山東 濟(jì)南 250022)

鋼柱抗爆響應(yīng)分析單自由度模型適用性評(píng)估*

李月強(qiáng)1,衣 娜2,席 豐1

(1.山東建筑大學(xué)土木工程學(xué)院,山東 濟(jì)南 250101; 2.山東軍之星建筑設(shè)計(jì)有限公司,山東 濟(jì)南 250022)

為評(píng)估單自由度(SDOF)模型在結(jié)構(gòu)抗爆設(shè)計(jì)中的適用性，分別采用SDOF模型和通用有限元軟件ANSYS/LS-DYNA對(duì)簡(jiǎn)支鋼柱承受爆炸荷載時(shí)的動(dòng)力響應(yīng)進(jìn)行模擬；對(duì)比二者計(jì)算結(jié)果，并以有限元模擬為準(zhǔn)，分析SDOF模型的適用范圍。研究表明:可按照自由振動(dòng)階段SDOF模型位移結(jié)果的振幅大小，將其位移響應(yīng)劃分為有限變形階段、臨界階段、失穩(wěn)破壞階段，有限變形階段SDOF模型與有限元結(jié)果基本一致；截面高寬比、翼緣寬厚比對(duì)鋼柱動(dòng)力破壞形式有重要影響，高寬比越大、翼緣的寬厚比越小，越容易發(fā)生平面外彎扭失穩(wěn)；在SDOF模型中通過(guò)假定塑性鉸分布長(zhǎng)度計(jì)算塑性階段應(yīng)變及應(yīng)變率，采用隨時(shí)間變化的應(yīng)變率計(jì)算Cowper-Symonds本構(gòu)關(guān)系中的應(yīng)力放大系數(shù)是可行的。

爆炸荷載;鋼柱;等效單自由度模型;有限元;應(yīng)變率

由于爆炸荷載傳播速度快、峰值壓力大、作用時(shí)間短，因此難以精確描述所引起的動(dòng)力響應(yīng)。為了方便結(jié)構(gòu)抗爆工程應(yīng)用，理想的方法是提出一些簡(jiǎn)化的分析模型，等效單自由度(single degree of freedom model, SDOF)模型就是其中之一。GB 50009-2012《建筑結(jié)構(gòu)荷載規(guī)范》中規(guī)定：由炸藥、燃?xì)?、粉塵等引起的爆炸荷載宜按等效靜力荷載采用，在確定該靜力荷載時(shí)按單自由度體系強(qiáng)迫振動(dòng)的方法分析得到構(gòu)件的內(nèi)力[1]。但是在單自由度模型的推導(dǎo)及求解過(guò)程中，引入了一些理論假設(shè)，使該模型只能近似反映結(jié)構(gòu)的實(shí)際行為，因此需要對(duì)它在抗爆工程中的適用性進(jìn)行評(píng)估。

對(duì)SDOF模型有過(guò)許多研究。A.A.Nassr等[2]通過(guò)爆炸實(shí)驗(yàn)，驗(yàn)證了SDOF模型中兩個(gè)重要的理論假設(shè)：動(dòng)力響應(yīng)的第一振型占主導(dǎo)地位；結(jié)構(gòu)構(gòu)件的變形在彈性階段結(jié)束后立即進(jìn)入塑性階段，不考慮彈塑性變形。實(shí)驗(yàn)數(shù)據(jù)和SDOF模型計(jì)算結(jié)果有比較好的一致性，但在SDOF模型中采用單一應(yīng)變率計(jì)及應(yīng)變率效應(yīng)；僅對(duì)實(shí)驗(yàn)工況進(jìn)行了模擬，沒(méi)有討論SDOF模型的適用范圍。

本文中，利用ANSYS/LS-DYNA[3]對(duì)鋼柱在爆炸荷載作用下的動(dòng)力響應(yīng)進(jìn)行精確求解，并以此為依據(jù)，詳細(xì)討論SDOF的適用性；在SDOF模型中，通過(guò)假定塑性鉸的分布長(zhǎng)度，計(jì)算塑性階段的應(yīng)變及應(yīng)變率，更合理地考慮應(yīng)變率效應(yīng)。

1 動(dòng)力響應(yīng)分析

分析圖1(a)所示簡(jiǎn)支鋼柱，先受到軸向載荷N的作用，然后受到橫向三角形脈沖載荷p(t)作用，如圖2所示。分別采用SDOF模型和有限元軟件進(jìn)行模擬，對(duì)比二者的柱中點(diǎn)位移、應(yīng)變、應(yīng)變率結(jié)果，分析SDOF模型在抗爆工程應(yīng)用中的適用性。

1.1SDOF模型分析

眾所周知，對(duì)于彈塑性梁的動(dòng)力響應(yīng)，當(dāng)假定其變形模態(tài)后，可簡(jiǎn)化為單自由度模型進(jìn)行分析。對(duì)于梁柱構(gòu)件，當(dāng)計(jì)入p-δ效應(yīng)后，同樣可歸結(jié)為SDOF模型分析問(wèn)題。

圖1(b)所示體系的運(yùn)動(dòng)方程為[4]：

(1)

(2)

(3)

塑性階段φ(ξ)采用在柱中點(diǎn)形成塑性鉸的雙直線：

φ(ξ)=1-2|ξ|ξ=z/L-1/2

(4)

式(1)～(2)及各系數(shù)，彈性階段可由哈密頓原理推導(dǎo)得出；塑性階段在柱中點(diǎn)形成塑性鉸，直接用平衡法列出運(yùn)動(dòng)微分方程，進(jìn)而得出各系數(shù)。SDOF模型的數(shù)值計(jì)算，采用Visual C++6.0編程[5-7]求解，先將運(yùn)動(dòng)方程離散為增量形式，用線加速度法將加速度、速度表示為位移增量形式，進(jìn)而寫(xiě)出控制方程，通過(guò)迭代可求出集中質(zhì)量的位移-時(shí)間歷程曲線。求解中為了考慮應(yīng)變率效應(yīng)，采用Cowper-Symonds本構(gòu)關(guān)系以及鋼材理想彈塑性假設(shè)描述應(yīng)力應(yīng)變關(guān)系：

(5)

SDOF模型中考慮軸向荷載N的3方面影響：(1)使中點(diǎn)處彎矩增大，通過(guò)等效橫向荷載η體現(xiàn)[4]；(2)改變了橫截面中性軸的位置，增大抗彎剛度EI；(3)根據(jù)壓彎構(gòu)件穩(wěn)定理論[10],即ymax=y0/(1-N/NE)，y0為簡(jiǎn)支梁的最大撓度，ymax為壓彎構(gòu)件最大撓度，軸向荷載使最大位移增大，也可認(rèn)為是減小了彈性彎曲剛度，(2)、(3)通過(guò)彈性剛度系數(shù)K體現(xiàn)。

彈性階段柱中點(diǎn)處截面邊緣軸向應(yīng)變?yōu)椋?/p>

κ(ξ,t)=2εmax(ξ,t)/h=φ″(ξ)y(t)

(6)

εmax=hφ″(ξ)y(t)/2=4.8hy(t)/L2

(7)

式中：κ、ε、h分別為曲率、應(yīng)變、橫截面高度[2]。塑性階段柱中點(diǎn)處截面邊緣軸向應(yīng)變?yōu)椋?/p>

εmax=[4y(t)/L](1/l)(h/2)=2hy(t)/(lL)

(8)

式中：l為塑性鉸的分布長(zhǎng)度,假定為2h[11]。柱中點(diǎn)處截面邊緣等效塑性應(yīng)變率為：

(9)

1.2DYNA計(jì)算模型

利用有限元軟件ANSYS/LS-DYNA進(jìn)行分析，采用全積分實(shí)體單元Solid 164建模，材料選用計(jì)及應(yīng)變率效應(yīng)的隨動(dòng)硬化理想彈塑性模型，材料參數(shù)為：密度7 850 kg/m3,彈性模量206 GPa，泊松比0.3，屈服強(qiáng)度345 MPa，Cowper-Symonds模型參數(shù)D=40,q=5，失效應(yīng)變0.2；采用三角形爆炸荷載，峰值壓力p隨不同工況而定，為基準(zhǔn)壓力p0=800 kPa的倍數(shù)，持續(xù)時(shí)間td=3 ms，軸向荷載N為鋼柱靜態(tài)軸向承載力的25%[4]，加載方式前者為瞬態(tài)類型，后者為同時(shí)考慮動(dòng)力釋放和瞬態(tài)分析類型[3]；為了方便實(shí)現(xiàn)鉸接約束并接近真實(shí)情況，在鋼柱兩端增加厚20 mm的端板，一端約束端板截面高度中心線上所有節(jié)點(diǎn)3個(gè)方向的位移，即ux=uy=uz=0，另一端約束相應(yīng)位置節(jié)點(diǎn)兩個(gè)方向的位移，即uy=uz=0[4]。

1.3DYNA模型與SDOF模型結(jié)果比較

鋼柱選用3種H型截面，分別為HM150×100(HM柱)、HW150×150(HW柱)、HN200×100(HN柱)[12]，計(jì)算長(zhǎng)度3 m，荷載同前，各工況只改變改變峰值壓力p。以下，*表示失穩(wěn)破壞。

柱中點(diǎn)處位移如圖3所示，具體數(shù)據(jù)見(jiàn)表1。通過(guò)比較，將SDOF模型的位移結(jié)果分為3個(gè)階段：有限變形階段、臨界階段、失穩(wěn)破壞階段。有限變形階段，即彈性變形和有限塑性變形，兩種方法所得ymax相差較小且SDOF的小于DYNA的，ymin相差較大且SDOF的大于DYNA的(ymax、ymin指自由振動(dòng)階段一個(gè)周期中的最大、最小位移)。這是由于SDOF模型比DYNA模型的總體剛度大，另外還有MP取值的影響。參照鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范，在抗爆設(shè)計(jì)中，應(yīng)該制定具體的破壞準(zhǔn)則，對(duì)爆炸荷載下構(gòu)件的位移限值或承載能力做出規(guī)定，本文中以鋼柱達(dá)到平面內(nèi)極限承載力時(shí)的爆炸荷載為pu，所對(duì)應(yīng)自由振動(dòng)階段的平衡位移為最終位移yu。DYNA所得HM、HW、HN柱pu分別為4.4p0、5.0p0、5.5p0，yu分別為175、236、170 mm； SDOF對(duì)應(yīng)的yu及振幅分別為164.5、9.5 mm，217、3 mm，158.5、11.5 mm。結(jié)合表1可以看出，兩種方法所得自由振動(dòng)階段的位移振幅隨著載荷強(qiáng)度的增大而減小。振幅大說(shuō)明截面彈性區(qū)大，承載力可繼續(xù)提高；振幅小則截面塑性區(qū)大，構(gòu)件臨近破壞。因此可以用振幅的大小界定各階段的范圍或制定相應(yīng)破壞準(zhǔn)則。有限變形階段與臨界階段分界的振幅限值，對(duì)HM、HW、HN柱，約為h/16、h/50、h/17。

臨界階段，DYNA位移結(jié)果顯示鋼柱已破壞，具體形式為平面外彎扭失穩(wěn)，如圖4所示；而SDOF結(jié)果則顯示鋼柱仍有一定承載能力。這是由于DYNA模型為空間模型，爆炸荷載引起的振動(dòng)使當(dāng)變形達(dá)到一定值時(shí)，鋼柱出現(xiàn)平面外位移繼而屈曲，此時(shí)截面并沒(méi)有完全進(jìn)入塑性；而SDOF模型則是以全截面屈服為承載力極限的，假定截面塑性可以完全發(fā)揮。以DYNA結(jié)果為準(zhǔn)，一般認(rèn)為SDOF模型不能用于確定鋼柱極限承載力pu。如前所述，DYNA所得HM、HW、HN柱的pu分別為4.4p0、5.0p0、5.5p0，而SDOF所得pu分別為4.7p0、5.0p0、6.5p0，可見(jiàn)不同截面柱兩模型所得pu的差大小不同，HN柱差別最大，HM柱次之，HW柱最小，其中包含截面高寬比、翼緣寬厚比的影響。DYNA模型中顯示HW柱受壓區(qū)破壞先于平面外屈曲，如圖4(a)所示，說(shuō)明柱失穩(wěn)之前已接近全截面屈服，主要是因?yàn)樵摻孛娓邔挶容^小、翼緣寬厚比較大，可以有效約束截面扭轉(zhuǎn)；而另外兩個(gè)截面柱的破壞形式如圖4(b)所示，平面外屈曲先于全截面屈服?？梢?jiàn)，爆炸荷載下截面高寬比越大，翼緣寬厚比越小的H型截面柱越容易發(fā)生平面外彎扭失穩(wěn)。

失穩(wěn)破壞階段兩種方法都能計(jì)算出鋼柱的失穩(wěn)破壞，不同的是SDOF模型中認(rèn)為破壞的起因是全截面屈服喪失承載力，而DYNA模型則一般為彎扭失穩(wěn)。結(jié)合臨界階段可以得出這樣的結(jié)論：若SDOF模型計(jì)算出鋼柱失穩(wěn)破壞則DYNA的結(jié)果與之相同，相反不一定成立。

另外，SDOF模型中應(yīng)變通過(guò)簡(jiǎn)單公式由位移得出，這種近似處理將導(dǎo)致兩種模型的應(yīng)變計(jì)算誤差大于位移計(jì)算誤差。在有限變形階段，SDOF模型的應(yīng)變結(jié)果比DYNA的結(jié)果小，這同樣是由于SDOF體系剛度較大。

由式(8)～(9)還可看出,SDOF中塑性鉸長(zhǎng)度對(duì)應(yīng)變及應(yīng)變率的計(jì)算有顯著影響。本文中假定塑性鉸長(zhǎng)度為兩倍的截面高度，是一個(gè)常數(shù)，這也就造成了3種截面柱的應(yīng)變誤差差別較大，尤其是HW柱5.0p0工況，SDOF的結(jié)果偏小，這說(shuō)明此時(shí)HW柱的塑性鉸長(zhǎng)度取值偏大。隨著載荷強(qiáng)度的增大，應(yīng)變誤差也在增大，于是可以認(rèn)為：爆炸載荷強(qiáng)度增大、截面塑性發(fā)展增大，塑性鉸分布長(zhǎng)度減小。若能確定具體工況的塑性鉸長(zhǎng)度，兩種模型的應(yīng)變計(jì)算會(huì)有較好的一致性。柱中點(diǎn)截面邊緣壓應(yīng)變曲線如圖5所示，具體數(shù)據(jù)見(jiàn)表2。圖5中，應(yīng)變?yōu)槭軌阂砭壿S向應(yīng)變：圖5(a)、(c)中，DYNA結(jié)果下降段是因?yàn)槠矫嫱馀まD(zhuǎn)變形使得受壓翼緣受拉，部分抵消了壓應(yīng)變；圖5(b)中，DYNA結(jié)果的陡降段是由于相應(yīng)單元在其應(yīng)變達(dá)到失效應(yīng)變0.2后被刪除引起的。

SDOF模型中應(yīng)變率是由應(yīng)變?cè)隽看_定的，而應(yīng)變率又會(huì)影響到塑性恢復(fù)力和體系的剛度，所以模型中應(yīng)變率的計(jì)算是至關(guān)重要的。柱中點(diǎn)截面應(yīng)變率曲線如圖6所示，具體數(shù)據(jù)見(jiàn)表3?？梢钥闯?，利用式(9)計(jì)算應(yīng)變率，與DYNA的結(jié)果基本一致：趨勢(shì)相符，數(shù)值偏小，且各種工況下兩個(gè)模型的應(yīng)變率誤差不同。這里也有塑性鉸分布長(zhǎng)度的影響，本文中側(cè)重點(diǎn)不在塑性鉸，因此不討論。應(yīng)變率曲線圖中第1個(gè)峰值說(shuō)明考慮應(yīng)變率效應(yīng)后屈服強(qiáng)度有所提高。

pp0εmax/10-4SDOFDYNAΔεmax/10-4εmin/10-4SDOFDYNAΔεmin/10-4HM150×100εmax/10-4SDOFDYNAΔεmax/10-4εmin/10-4SDOFDYNAΔεmin/10-4HW150×150εmax/10-4SDOFDYNAΔεmax/10-4εmin/10-4SDOFDYNAΔεmin/10-4HN200×1000.5913-4-9-5-4812-4-8-4-4812-4-8-4-41.01822-4-20-14-61821-3-17-13-41721-4-18-13-52.06769-2252056369-6221935060-108533.0184181315213715168178-101341340127147-208996-74.0374377-33513501339362-23313326-13239263-24205216-115.0******6551069-4146501060-410406444-38380405-25

表3 柱中點(diǎn)截面最大應(yīng)變率Table 3 Mid-span maximum strain rates

2 結(jié) 論

通過(guò)有限元模擬與SDOF模型結(jié)果的比較，著重討論了爆炸荷載作用下簡(jiǎn)支鋼柱SDOF模型的適用性，主要結(jié)論有如下幾點(diǎn)。

(1)可將SDOF模型計(jì)算的位移結(jié)果分為3個(gè)區(qū)段：有限變形階段、臨界階段、失穩(wěn)破壞階段?？贡O(shè)計(jì)中可根據(jù)自由振動(dòng)階段位移振幅大小劃分各階段或建立破壞準(zhǔn)則。有限變形階段，SDOF模型適于工程應(yīng)用，該有限變形對(duì)于HM150×100柱、HW150×150柱、HN200×100柱，約為1.1h、1.4h、0.8h(h為截面高度)。

(2)截面高寬比、翼緣寬厚比對(duì)鋼柱動(dòng)力破壞形式有重要影響，高寬比越大、翼緣寬厚比越小，越容易發(fā)生平面外彎扭失穩(wěn)。一般情況下，SDOF模型不能用于確定鋼柱平面內(nèi)極限承載力，而當(dāng)截面高寬比較小、翼緣寬厚比較大，構(gòu)件平面外扭轉(zhuǎn)變形能有效約束時(shí)，SDOF模型與有限元模擬所得極限承載力差別較小。

(3)本文中SDOF模型采用隨時(shí)間變化的應(yīng)變率計(jì)算Cowper-Symonds本構(gòu)關(guān)系中的應(yīng)力放大系數(shù)，該應(yīng)變率與有限元模擬結(jié)果差別較小，說(shuō)明通過(guò)假定塑性鉸分布長(zhǎng)度計(jì)算塑性階段的應(yīng)變及應(yīng)變率是可行的。

[1] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.GB 50009-2012 建筑結(jié)構(gòu)荷載規(guī)范[S].北京:中國(guó)建筑工業(yè)出版社,2012.

[2] Nassr A A, Razaqpur A G. Single and multi degree of freedom analysis of steel beams under blast loading[J]. Nuclear Engineering and Design, 2012,242:63-77.

[3] LS-DYNA keyword user’s manual: Version 971[Z]. Livermore, California: Livermore Software Technology Corporation, 2007.

[4] Nassr A A, Razaqpur A G. Strength and stability of steel beam columns under blast load[J] .International Journal of Impact Engineering, 2013,55(5):34-48.

[5] 席豐,張?jiān)?脈沖載荷作用下鋼梁動(dòng)力響應(yīng)及反常行為的應(yīng)變率效應(yīng)[J].爆炸與沖擊,2012,32(1):34-42.

Xi Feng, Zhang Yun. The effects of strain rate on the dynamic response and abnormal Behavior of steel beams under pulse loading[J]. Explosion and Shock Waves, 2012,32(1):34-42.

[6] 劉鋒,席豐.子彈撞擊作用下固支淺圓拱的彈塑性動(dòng)力響應(yīng)[J].爆炸與沖擊,2005,25(4):361-367.

Liu Feng, Xi Feng. Elastic-plastic dynamic response of a fully clamped shallow arch subjected to projectile impact[J]. Explosion and Shock Waves, 2005,25(4):361-367.

[7] 譚浩強(qiáng).C程序設(shè)計(jì)(第四版)[M].北京:清華大學(xué)出版社,2010.

[8] Nassr A A, Razaqpur A G. Dynamic response of steel columns subjected to blast loading[J]. Journal of Structural Engineering, 2014,140(7):165-180.

[9] Nassr A A, Razaqpur A G. Experimental performance of steel beams under blast loading[J]. Journal of Performance of Constructed Facilities, 2012,26(26):600-619.

[10] 陳驥.鋼結(jié)構(gòu)穩(wěn)定理論與設(shè)計(jì)(第六版)[M].北京:科學(xué)出版社,2014.

[11] 余同希,斯壯 W J.塑性結(jié)構(gòu)的動(dòng)力學(xué)模型[M].北京:北京大學(xué)出版社，2002.

[12] 陳紹蕃,顧強(qiáng).鋼結(jié)構(gòu)(第二版)[M].北京:中國(guó)建筑工業(yè)出版社,2007.

Abstract: For the evaluation of the applicability of the single degree of freedom (SDOF) model in the structural antiknock design, the dynamic response of the simply supported steel column under explosion load was simulated using both the SDOF model and the ANSYS/LS-DYNA in this paper. By the comparison of the two calculation results, the scope of application of the SDOF model was analyzed according to the finite element simulation. The results show that the displacements calculated using the SDOF model can be divided into three different phases including the finite deformation, in which the SDOF model agrees well with the DYNA simulation, the critical deformation, and the buckling failure deformation, according to the amplitude size in the free vibration. The ratio of the cross section’s depth to its width and that of the flange’s width to its thickness have significant effect on the dynamic failure forms of the steel column, namely the bigger the ratio of the depth to the width and the smaller the ratio of the width to the thickness, the more prone it is for the buckling to suffer out-of-plane bending and twisting. In the SDOF model, it is feasible to calculate the strain and the strain rate in the plastic deformation phase by assuming the plastic hinge distribution length and the stress-magnified coefficient in the Cowper-Symonds constitutive relation by using the time-dependent strain rate.

Keywords: blast load; steel column; single degree of freedom model; finite element; strain rate

(責(zé)任編輯 丁 峰)

Assessmentonsingledegreeoffreedommodelinsteelcolumnanalysisofanti-detonation

Li Yueqiang1, Yi Na2, Xi Feng1

(1.CivilEngineeringInstitute,ShandongJianzhuUniversity,Jinan250101,Shandong,China; 2.ShandongJunzhixingArchitecturalDesignLimitedLiabilityCompany,Jinan250022,Shandong,China)

O342;TU391國(guó)標(biāo)學(xué)科代碼1301565

A

10.11883/1001-1455(2017)05-0957-07

2016-01-11；

2016-05-28

國(guó)家自然科學(xué)基金項(xiàng)目(11272189)

李月強(qiáng)(1986— )，男，碩士研究生;

席 豐,xifeng@sdjzu.edu.cn。