基于LS-DYNA的EPS墊層棚洞耗能減震分析★

袁 松 黎良仆

(四川省交通運輸廳交通勘察設計研究院，四川 成都 610017)

基于LS-DYNA的EPS墊層棚洞耗能減震分析★

袁 松 黎良仆

(四川省交通運輸廳交通勘察設計研究院，四川 成都 610017)

以某在建項目為背景，基于LS-DYNA建立了數值計算模型，模擬了落石沖擊棚洞結構的全過程，揭示了EPS緩沖材料緩沖耗能作用，結果表明：EPS墊層棚洞較傳統砂墊層棚洞抗落石沖擊性能有一定提高，能廣泛應用于實際工程中。

鋼筋混凝土棚洞，EPS聚苯乙烯墊層，落石，動力有限元

0 引言

我國西部地區地勢險峻，落石災害時有發生，對依山而建的大量公路造成了嚴重威脅。棚洞作為落石防護結構，耐久性高、能夠與環境協調發展，故而在西部山區公路廣泛應用。

JTG/T D70—2010公路隧道設計細則[1]中對棚洞洞頂回填做出了明確規定：當存在落石崩塌危害時，洞頂回填土厚度不宜小于1.5 m，對落石災害嚴重區域應酌情加厚。這樣導致設計時截面尺寸增大，配筋量增加，不利于結構抗震，建設成本也不易控制。通過在棚洞頂部混凝土板與回填土層之間鋪設一層緩沖材料，既減輕棚洞自重，又能吸收落石沖擊能量，保護棚洞結構，故成為該領域研究熱點。Zarnani等[2]首先在擋土墻與土體之間填充EPS(聚苯乙烯)層，孫洪利等[3]在超前支護樁后設置EPS墊層，結果表明采用EPS墊層后最大土壓力都大大減小；何思明等[4]將落石沖擊EPS墊層棚洞看成是正交各向異性復合板沖擊問題，求解了EPS墊層棚洞鋼筋混凝土板的沖擊動力響應理論解。縱觀上述研究，缺少實際工程下的整體性研究。本文在前人基礎上，基于LS-DYNA，結合實際在建項目，模擬落石沖擊棚洞結構全過程，以期為工程實際提出建議。

1 有限元數值模擬實驗

研究對象取四川省交通運輸廳交通勘察設計研究院設計的G245線烏斯河至甘洛縣城段公路改建工程棚洞，按12 m×6 m一跨足尺模型和實際配筋進行建模。實際模型配筋圖如圖1所示。

由于落石沖擊棚洞涉及復雜的功能轉化，研究過程中忽略對結果影響較小的因素，作出如下假設：1)不考慮沖擊過程中落石的破壞和回填土層飛濺的現象；2)不考慮落石斜向沖擊棚洞的情況。

1.1 單元及材料

表1 HJC C30混凝土材料模型材料參數及取值

表2 材料參數

材料密度/kg·m-3彈模/MPa屈服極限/MPa泊松比鋼筋78002×1054000.3砂土2000390.1480.3EPS泡沫203.50.140.1落石30002.5×104—0.2

1.2 邊界條件

砂墊層與泡沫、泡沫與棚洞頂板之間采用面面自動接觸，上部T型梁與下部頂端梁之間采用面面普通接觸，頂端梁下部節點全部約束，固定所有自由度。

1.3 模態分析

沖擊問題中阻尼對結果影響較大。LS-DYNA 中阻尼力通過式(1)計算施加到模型的每個節點[8]。

(1)

其中,Ds為阻尼系數；m為結構質量；v為運動速度；Ds的算法由式(2)給出：

(Ds)critical=2ωmin

(2)

其中,ωmin為該結構的基頻，故需進行模態分析得到結構的基頻。

本文采用整體式模型，將鋼筋分布于整個單元中，采用Solid65單元，鋼筋的貢獻度通過調整單元的材料力學參數得以實現[9]，混凝土選擇線彈性本構關系。由此得到結構基頻為：30.298 Hz。

棚洞結構一階振型見圖3。

根據ωn=f×2π，可由計算結果得到ωmin，再根據式(2)求得阻尼系數。

綜上建立有限元模型如圖4所示，整體模型共有161 855個節點，46 211個Beam單元，115 640個Solid單元。

2 數值實驗結果分析

計算可得棚洞結構受落石沖擊作用下任意時間各點處鋼筋應力值、撓度時程曲線如圖5所示。分別計算4 t落石從40 m高度落下沖擊鋼筋混凝土棚洞，EPS墊層填充厚度分別取0.3 m，0.6 m(即1.5 m砂土+0.3 m EPS;1.5 m砂土+0.6 m EPS)，各工況下落石動能時程曲線繪制如圖6所示。

如圖6所示，EPS不同填充厚度下落石能量衰減的曲線有一定差異，0.6 m EPS的曲線向右上方偏移，兩條曲線均出現了三個明顯的拐點，本文認為這是由EPS材料力學性質決定：彈性段吸收較小能量，進入屈服階段后，能量—時間曲線斜率逐漸變陡，表明進入能量吸收的主要階段，而后斜率逐漸降低，材料進入硬化階段，材料被逐漸壓實，耗能水平有所下降。0.6 m EPS曲線與0.3 m EPS曲線相比，拐點有所滯后，即材料進入硬化階段時，吸收能量比0.3 m EPS多，整個耗能階段的時間更長，故隨著填充緩沖材料厚度的增加，結構耗能更加明顯。

通過計算，棚洞結構不填充EPS墊層時，4 t落石從40 m落高落下，T梁主筋應力值達到401.2 MPa，跨中撓度值為0.042 m，其結果見表3，達到HRB400鋼筋屈服極限以及梁撓度限值δ<12/300=0.04 m，即達到棚洞結構極限承載能力。填充0.6 m EPS后，4 t落石從50 m落下沖擊棚洞，棚洞T型梁撓度值達到0.041 m，由此可見，采用1.5 m砂土+0.6 m EPS墊層棚洞的安全性能是1.5 m砂土墊層棚洞的1.2倍。

表3 數值實驗結果對比

3 結語

本文以實際在建工程為背景，重點研究砂墊層+不同厚度EPS的復合墊層在落石沖擊荷載下棚洞抗落石沖擊性能和耗能減震原理，基于數值模擬實驗得到如下結論：1)棚洞頂部設置砂+EPS的復合墊層形式能夠有效降低落石沖擊下結構受力鋼筋的應變值與主梁撓度值。2)EPS墊層耗能減震能力隨著EPS填充厚度的提高而提高，棚洞結構抗落石沖擊性能得到加強，這與EPS材料本身力學性能有關。3)以規范要求最低砂土墊層厚度為基準，采用填充EPS墊層可以在較大落石沖擊能量下有效降低結構自重及建設成本，是一種安全、經濟的棚洞墊層材料。

[1] JTG/T D70—2010，公路隧道設計細則[S].

[2] Zarnani S,Bathurst R J.Numerical parametric study of a ex-panded polystyrene(EPS)geofoam seismic buffers[J].Canadian Geotechnical Journal,2009,46(3):318-338.

[3] 孫洪利,劉元雪,葉四橋,等.超前支護樁樁后EPS墊層優化設置的數值模擬[J].重慶交通大學學報(自然科學版),2011(4):786-789，879.

[4] 王東坡,何思明,歐陽朝軍,等.滾石沖擊荷載下棚洞鋼筋混凝土板動力響應研究[J].巖土力學,2013(3):881-886.

[5] 汪 敏,石少卿,陽友奎.柔性棚洞在落石沖擊作用下的數值分析[J].工程力學,2014(5):151-157.

[6] 熊益波，陳劍杰，胡永樂，等.混凝土Johnson-Holmquist本構模型關鍵參數研究[J].工程力學，2012(1):121-127.

[7] N Kishi,H Konno,K Ikeda,etc..Prototype impact tests on ultimate impact resistance of PC rock-sheds[J].International Journal of Impact Engineering,2002,27(9):969-985.

[8] LS-DYNA ?KEYWORD USER’S MANUAL[Z].2005.

[9] 賈 濤，劉世忠.基于ANSYS實體模型對小跨徑鋼筋混凝土簡支梁舊橋的模態分析[J].蘭州交通大學學報，2008(3):9-12.

Study on the energy dissipation of shed with EPS cushion by LS-DYNA★

Yuan Song Li Liangpu

(SichuanCommunicationSurveying&DesignInstitute,Chengdu610017,China)

Taking the in-built project as the background, the thesis establishes numerical calculation model on the basis of LS-DYNA, simulates the whole process of rockfall shocking hanger tunnel, and demonstrates the damping and energy dissipation effect of EPS cushioning materials. Results show that: comparing to sandy cushion hanger tunnel, the shock-resisting performance of EPS cushion hanger tunnel has been improved, which can be widely applied in actual engineering.

steel reinforce concrete hanger tunnel, EPS polystyrene cushion, rockfall, dynamic finite element

1009-6825(2017)02-0050-02

2016-11-04★:四川省交通科技項目(項目編號：2012B4-2)

袁 松(1983- )，男，碩士，工程師； 黎良仆(1990- )，男，助理工程師

TU375.4

A