天然不含節(jié)竹管在軸向壓縮沖擊載荷下沖擊響應(yīng)與吸能特性的數(shù)值分析

收稿日期：2022-02-25""" 修回日期：2022-07-04

基金項目：國家自然科學(xué)基金資助項目（No.11872218；11572161）；浙江省“近海結(jié)構(gòu)沖擊安全防護與健康監(jiān)測”重點科技創(chuàng)新團隊（No.2013TD21）

通信作者：繆馥星，教授。E-mail：miaofuxing@nbu.edu.cn

引用格式：

周盈俠，繆馥星，徐挺，等.天然不含節(jié)竹管在軸向壓縮沖擊載荷下沖擊響應(yīng)與吸能特性的數(shù)值分析［J］.應(yīng)用力學(xué)學(xué)報，2024，41（3）：629-641.

ZHOU Yingxia，MIAO Fuxing，XU Ting，et al.Numerical analysis of impact response and energy absorption characteristics of natural internode bamboo tubes under axial compression impact load［J］.Chinese journal of applied mechanics，2024，41（3）：629-641.

文章編號：1000-4939（2024）03-0629-13

摘" 要：天然不含節(jié)竹管在工程中應(yīng)用較廣，研究此類天然薄壁管結(jié)構(gòu)的抗沖擊能力具有重要的意義。利用分離式霍普金森壓桿（split Hopkinson pressure bar，SHPB），研究天然竹管的抗沖擊性能并進行有限元分析。在ABAQUS中應(yīng)用三維Hashin和Puck單向纖維復(fù)合材料失效準則，結(jié)合應(yīng)變率效應(yīng)對強度的修正建立竹管本構(gòu)模型，對沖擊載荷下的天然不含節(jié)竹管進行模擬，得到的動態(tài)應(yīng)力-應(yīng)變曲線和破壞模式與SHPB實驗接近，驗證有限元模型的準確性。結(jié)果表明：實驗與模擬的極限抗壓強度誤差最小為2.19%、最大為8.68%；在不同的沖擊壓縮載荷下，竹管出現(xiàn)端面壓潰、沿軸向劈裂等破壞形式；竹管外壁節(jié)點處的位移、速度和加速度響應(yīng)均隨沖擊載荷的增加而增加；隨沖擊載荷升高，入射能、反射能、吸收能和比吸能隨時間均呈上升趨勢，透射能隨時間變化不大；不含節(jié)竹管的峰值載荷遠低于薄壁圓管 GFRP （glass fiber reinforcement plastic），但比吸能較為接近，說明不含節(jié)竹管并不適用于需承受大沖擊能量的高速沖擊防護領(lǐng)域。初步結(jié)果可為竹管用于工程結(jié)構(gòu)防護給予理論參考。

關(guān)鍵詞：不含節(jié)竹管；軸向沖擊；沖擊響應(yīng)；吸能特性；有限元分析

中圖分類號：TU366.1" 文獻標志碼：A

DOI：10.11776/j.issn.1000-4939.2024.03.016

Numerical analysis of impact response and energy absorption characteristics of natural internode bamboo

tubes under axial compression impact load

ZHOU Yingxia，MIAO Fuxing，XU Ting，JIN Yue

（Key Laboratory of Impact and Safety Engineering，Ministry of Education，Ningbo University，315211 Ningbo，China）

Abstract：Natural internode bamboo tubes are commonly used in engineering，and it is of great significance to investigate the impact resistance of such natural thin-walled tube structures.The impact resistance of natural bamboo tubes was investigated，and a finite element analysis was carried out with a split Hopkinson pressure bar （SHPB）.The three-dimensional Hashin and Puck unidirectional fiber composite failure criterion was applied in ABAQUS and combined with the correction of the strain rate effect on strength to establish a bamboo tube constitutive model to simulate the natural internode bamboo tubes under impact loading.The dynamic stress-strain curves and damage modes obtained were close to those obtained from the SHPB experiments，proving the accuracy of the finite element model.The results show that the experimental and simulated ultimate compressive strengths have a minimum error of

2.19 percent and a maximum error of 8.68 percent; under different impact compression loads，the bamboo tubes show end face crushing，axial splitting，and other forms of damage; displacement，velocity，and acceleration responses at the nodes of the outer wall of the bamboo tubes increase with the increase of the impact load.The incident energy，reflected energy，absorbed energy，and specific absorbed energy all increase with time，whereas the transmitted energy does not.The peak load of internode bamboo tubes is significantly lower than that of GFRP thin-walled round tubes，but the specific absorbed energy is similar，indicating that internode bamboo tubes are not suitable for high-speed impact protection fields that require large impact energy.The preliminary findings may provide some theoretical guidance for the use of bamboo tubes in engineering structural protection.

Key words：internode bamboo tube; axial impact; impact response; energy absorption characteristic; finite element analysis

竹子是一種天然的纖維材料，也是一種薄壁結(jié)構(gòu)，具有不同于金屬類材料及其薄壁管的獨特的力學(xué)性能［1-6］，工程應(yīng)用較廣，研究竹管在沖擊載荷下的響應(yīng)及吸能性能，對竹管在建筑工程中的應(yīng)用有重要的理論指導(dǎo)意義。

目前，關(guān)于竹材力學(xué)性能的研究較多，譬如，文獻［2］通過準靜態(tài)實驗，研究不同加載方向下竹材的力學(xué)性能，以及通過分離式霍普金森壓桿（split Hopkinson pressure bar，SHPB）實驗，研究應(yīng)變率對竹材力學(xué)性能的影響，并分析了竹材的破壞機理。文獻［7］通過動態(tài)拉伸實驗，發(fā)現(xiàn)了竹材的彈塑性較弱，易在屈服階段發(fā)生脆性斷裂；且發(fā)現(xiàn)竹材的比吸能與其含水率成正比，跟竹齡成反比。隨含水率升高，比吸能可增加到18.77J/g；隨竹齡增加，比吸能可減小到3.71J/g。文獻［8］通過落錘試驗，發(fā)現(xiàn)在軸向沖擊載荷作用下，竹管的損傷主要為屈曲鼓包和加載端的開裂破壞，其吸能特性隨徑厚比的減小而增加；同時發(fā)現(xiàn)動態(tài)載荷下竹管比吸能較靜態(tài)試驗下的大，隨徑厚比不同，約在13～38J/g。

在數(shù)值模擬方面，文獻［9-10］建立1/4、文獻［11］建立1/2不含節(jié)竹管模型，應(yīng)用脆性斷裂本構(gòu)和對稱邊界條件，在軸向壓縮載荷作用下，竹管出現(xiàn)軸向劈裂破壞，與實驗結(jié)果相比，其最大誤差分別為6%和16%。文獻［12］建立了不含節(jié)竹管模型，應(yīng)用各向異性本構(gòu)，在軸向載荷下，得出最大Von Mises應(yīng)力和剪應(yīng)力均出現(xiàn)在加載端，表明該區(qū)域發(fā)生壓潰變形且易萌生裂紋。文獻［13］通過改變代表性體積元中的纖維含量，推導(dǎo)出了不同纖維百分比下的各向異性參數(shù)，用于竹管軸向壓縮模擬，發(fā)現(xiàn)軸向間隔排列的外脊和內(nèi)隔膜可增加竹管抗屈曲性能。另外，根據(jù)結(jié)構(gòu)仿生原理，進行薄壁管和晶胞管的仿竹設(shè)計時，多是基于金屬材料（鋁合金［14-17］、不銹鋼［18］）的設(shè)計及優(yōu)化的結(jié)果，并未能揭示天然竹管的沖擊響應(yīng)及吸能特性。然而，關(guān)于天然竹管在軸向壓縮沖擊載荷下沖擊響應(yīng)與吸能特性的分析尚未見報道。

本研究運用三維Hashin和Puck單向復(fù)合材料損傷VUMAT材料子程序構(gòu)建竹管材料模型，分析其在沖擊壓縮載荷下的沖擊響應(yīng)和吸能特性。對比不同沖擊載荷下有限元模擬和實驗中試樣的破壞形態(tài)和動態(tài)應(yīng)力-應(yīng)變曲線，驗證有限元模型的有效性。并分析不同徑長比竹管在沖擊壓縮載荷下的破壞模式、管外壁節(jié)點的位移、速度、加速度響應(yīng)和吸能特性。本研究的初步結(jié)果可為竹管用于工程結(jié)構(gòu)防護給予理論參考。

1" 竹管的軸向壓縮沖擊實驗

采用SHPB實驗設(shè)備，其材質(zhì)均為高彈鋼，桿徑為74mm，其子彈、入射桿和透射桿的長分別為：400、2700和1500mm。記竹管的高、外徑和壁厚分別為h、d和δ，并定義徑長比為d/h，厚徑比為δ/d，挑選出各沖擊氣壓下重復(fù)性較好的實驗數(shù)據(jù)進行分析，幾何尺寸如表1所示。分別選取0.2、0.3和0.4MPa 3種氣壓，進行沖擊壓縮實驗。如圖1所示，在選取的氣壓下，子彈撞擊入射桿，入射桿上的電阻應(yīng)變片測得入射波信號εI（t）和反射波信號εR（t），透射桿上半導(dǎo)體應(yīng)變片測得透射波信號εT（t），經(jīng)超動態(tài)應(yīng)變儀和示波器，得到?jīng)_擊加載波形，如圖2所示。實驗中采用銅環(huán)作為波形整形器，用以降低橫向彌散效應(yīng)。天然竹管的波阻抗與鋼的相差較大，用靈敏系數(shù)為110的半導(dǎo)體應(yīng)變片采集透射波信號，可避免因?qū)嶒灅悠放c投射桿波阻抗差距較大而導(dǎo)致的透射波信號微弱的問題。

根據(jù)SHPB實驗的一維應(yīng)力波和均勻化假定，得到試樣的應(yīng)變率ε·S（t）、應(yīng)變εS（t）、應(yīng)力σS（t）和平均應(yīng)變率ε-·分別為［19］

ε·S（t）=C0h［εI（t）－εR（t）－εT（t）］（1）

εS（t）=C0h∫t0［εI（t）－εR（t）－εT（t）］dt（2）

σS（t）=E0A02AS［εI（t）+εR（t）+εT（t）］（3）

ε-·=∫ε0ε·（ε）dεε（4）

式（1）～（4）中：A0、AS分別為壓桿和試樣的橫截面積；C0、E0為壓桿的縱波速和彈性模量。消去式（2）～（3）中的時間參數(shù)，即可求得試樣的動態(tài)應(yīng)力-應(yīng)變曲線。

根據(jù)文獻［20］，SHPB實驗中，一維應(yīng)力波假定需滿足

λ=cI·τIgt;5D（5）

式中：λ是入射波寬度；D是桿徑； cI、τI分別為入射波波速、傳播時間。這里波速約為5063m/s，從圖2中可看出，入射波的傳播時間大約從200μs到550μs，共350μs。即λ約為1.77m，大于5倍桿徑0.37m。經(jīng)檢驗，近似滿足一維應(yīng)力波假定。

對最長試樣1-1進行應(yīng)力平衡的分析，試樣中的波速cs約為4300m/s，一次反射所需時間2h/cs約為34.47μs，入射波上升時間約從200μs到375μs，共175μs，大于應(yīng)力波在試樣破壞前所需反射時間137.88μs，近似滿足應(yīng)力平衡假定。同理，其他稍短試樣也均近似滿足應(yīng)力平衡條件。

2" 不含節(jié)竹管軸向壓縮沖擊的有限元分析模型

根據(jù)SHPB實驗過程及其原理［19］，建立竹管在軸向沖擊載荷作用下的有限元模型，如圖3所示。其中，壓桿的彈性模量、密度和泊松比分別為200GPa、7800kg/m3和0.3。選取樣

品0.4MPa情形下的實驗樣品尺寸作為模擬的幾何尺寸進行模擬。記EX、EY、EZ、Xt、Xc、Yt、Yc分別為X、Y、Z方向的彈性模量、X方向拉伸強度和壓縮強度、Y方向拉伸強度和壓縮強度，GXY、GXZ、GYZ、μXY、μXZ、μYZ分別為XOY、XOZ、YOZ面內(nèi)剪切模量和泊松比，SXY、SXZ分別為XOY、XOZ平面剪切強度，ρ為密度，材料參數(shù)如表2所示［13，21-23］。根據(jù)文獻［24-25］測得的竹纖維平均體積分數(shù)為22.745%，本研究選取文獻［13］中竹纖維體積分數(shù)為20%時的參數(shù)進行模擬分析。鋪層方式為［0］5，選用C3D8R單元，設(shè)置通用接觸，所建立的1/4模型中入射桿、透射桿和d/h為0.49、1.03、5.98，竹管的單元數(shù)分別為4050、2250和19950、9250、1800。實驗測得的應(yīng)力脈沖作為數(shù)值模擬加載條件，如圖4所示。根據(jù)一維應(yīng)力波的傳播理論，不考慮試樣與壓桿接觸面間的摩擦效應(yīng)［19］。應(yīng)用ABAQUS，分析d/h為0.49、1.03和5.98竹管分別在3種載荷下的沖擊響應(yīng)及其吸能特性。由于竹管的材料參數(shù)受竹齡、生長地、含水率等的影響較大，本研究結(jié)合文獻［13，21-23］參數(shù)，并在合理范圍內(nèi)微調(diào)，使模擬與實驗的破壞和動態(tài)應(yīng)力-應(yīng)變曲線接近，初步說明竹管的沖擊壓縮有限元計算模型的可靠性。

由于竹子不同方向的強度和剛度值各不相同；同一方向的抗拉和抗壓強度存在差異；在壓縮載荷下會產(chǎn)生軸向劈裂等，這些力學(xué)行為與單向復(fù)合材料較類似。研究者們已提出許多預(yù)測復(fù)材失效的理論，其中Hashin三維和Puck準則預(yù)測單向橫觀各向同性復(fù)材失效比較準確［26］。采用參考文獻［27］中的三維Hashin和Puck損傷起始準則、單元刪除和剛度退化準則，文獻［28］的應(yīng)變率效應(yīng)，建立竹管損傷本構(gòu)模型。

記Ftf、Fcf、Ftm、Fcm分別為纖維縱向拉伸失效、纖維縱向壓縮失效、基體橫向拉伸失效和基體橫向壓縮失效的表征參數(shù)，則

Ftf（σ）=σXXXt2+σXYSXY2+σXZSXZ2（6）

Fcf（σ）=σXXXc（7）

Ftm（σ）=σXX2Xt2+σYY2YtYc+σXYSXY2+

σYY1Yt+1Yc （8）

Fcm（σ）=σXX2Xt2+σYY2YtYc+σXYSXY2+

σYY1Yt+1Yc （9）

式（6）～（9）中：σXX和σYY為正應(yīng)力分量，σXY和σXZ為剪應(yīng)力分量。當材料的纖維縱向拉伸、纖維縱向壓縮、基體橫向拉伸、基體橫向壓縮的表征參數(shù)不小于1時，即為計算單元開始損傷，小于1時，即為計算單元無損傷。

模型中應(yīng)用的應(yīng)變率效應(yīng)［28］如下所示。

S=S01+Clnε·ε·0（10）

式中：ε·為當前應(yīng)變率；ε·0為初始應(yīng)變率，一般為1s-1［29］；S0為初始強度；C為應(yīng)變率效應(yīng)系數(shù)，取0.1［28］。

在有限元計算中，為評估計算網(wǎng)格的敏感性，以0.4MPa，d/h=0.49的竹管為例，選取不同網(wǎng)格尺寸進行網(wǎng)格敏感性分析，如表3所示。發(fā)現(xiàn)隨竹管全模型單元數(shù)增加，相對誤差減小，算例2與算例1的相對誤差僅相差0.25%，算例2的相對誤差已趨于穩(wěn)定，如圖5所示，4種網(wǎng)格下的應(yīng)力-應(yīng)變曲線基本一致。綜合考慮實驗和模擬的極限抗壓強度的相對誤差、單元數(shù)量和計算時長，最終選取算例2的網(wǎng)格尺寸用于后續(xù)計算。

將模型中的入射應(yīng)力σI與反射應(yīng)力σR相加的疊加應(yīng)力σI+σR與透射應(yīng)力σT比較，可得到在3種沖擊載荷下不同d/h竹管的實驗與模擬的動態(tài)應(yīng)力平衡圖，進行應(yīng)力平衡條件的驗證。譬如，0.4MPa沖擊氣壓下，d/h=0.49情形的動態(tài)應(yīng)力平衡如圖6所示。由圖6可知，實驗與模擬中的動態(tài)應(yīng)力狀態(tài)均近似滿足其應(yīng)力平衡條件。

以沖擊氣壓0.4MPa，d/h=5.98竹管為例，取入射桿中距入射加載端面100、200、300和400mm，坐標分別為J（100，0，37）、K（200，0，37）、L（300，0，37）和M（400，0，37）4個節(jié)點，采集的入射波信號如圖7所示。

由表4知，4個節(jié)點的入射波波形寬度與峰值變化比均小于5%。根據(jù)文獻［30］，可忽略橫向彌散效應(yīng)的影響。

3 "軸向壓縮沖擊載荷下不含節(jié)竹管的沖擊力學(xué)行為分析

3.1" 動態(tài)應(yīng)力-應(yīng)變關(guān)系

竹管在不同應(yīng)變率下的動態(tài)應(yīng)力-應(yīng)變曲線如圖8所示。由圖8知，這與壓縮載荷下，脆性材料的破壞應(yīng)變約在1%或更低［31］的結(jié)論相一致。整體上，本研究計算分析模型的數(shù)值模擬結(jié)果與實驗結(jié)果吻合良好。其中，實驗顯示竹管在應(yīng)變率60s-1，εb為

0.00605，即竹管在很短的時間內(nèi)破壞。不妨記極限抗壓強度所對應(yīng)的應(yīng)變?yōu)棣舃，當應(yīng)變小于εb時，應(yīng)力隨應(yīng)變的增加而上升，應(yīng)力-應(yīng)變曲線表現(xiàn)為凸向上的特征；當應(yīng)變大于εb時，模擬結(jié)果較實驗結(jié)果的切線模量大，這主要是由于模擬計算中，當材料達到失效應(yīng)力后，計算單元損傷失效，不再承受載荷，而實驗中，材料在失效后還能承受一定大小的載荷，且計算模型中的竹管材料參數(shù)與實驗中樣品的材料參數(shù)有差異。

在破壞應(yīng)變εb之前，實驗曲線的斜率高于模擬曲線，這主要是因為選取的參數(shù)與實驗試樣有所差異；實際建模時是根據(jù)文獻［13］取了平均體積分數(shù)為20%時彈性模量整體賦值計算，沒有考慮彈性模量沿壁厚變化。

由表5可見，極限抗壓強度的相對誤差，最大為8.68%，最小為2.19%。這比文獻［9-11，32］應(yīng)用脆性斷裂本構(gòu)，竹管在軸向載荷下的劈裂破壞的模擬結(jié)果與實驗的誤差最小為2%、最大為20%的結(jié)論略好。

3.2" 沖擊破壞模式

如圖9所示，計算模型的預(yù)測結(jié)果與此實驗結(jié)果的一致性較好。

在0.2MPa沖擊氣壓下，竹管受壓均出現(xiàn)輕微的壓潰破壞。

2000 μs" from experiment and simulation

其中，d/h為0.52和1.08情形的竹管沿軸向劈裂破壞，且劈裂線呈平行性；隨著沖擊載荷增加50%，壓潰破壞更明顯，沖擊端出現(xiàn)擠壓外翻式的破壞，d/h為0.52和1.05情形的竹管劈裂條數(shù)增加，d/h為6.69情形的竹管壓潰成圓片狀；當沖擊載荷增加1倍時，d/h為0.49和1.03情形的竹管出現(xiàn)劈裂、面外翻折斷裂和端面壓潰等破壞現(xiàn)象，d/h為5.98情形的竹管被壓縮成薄碎片狀。

4" 軸向壓縮沖擊載荷下不含節(jié)竹管的沖擊響應(yīng)

d/h為0.49、1.03和5.98竹管分別在3種沖擊載荷下的位移、速度和加速度響應(yīng)曲線分別如圖10～12所示。

本節(jié)進行竹管在軸向沖擊壓縮載荷作用下的位移、速度和加速度響應(yīng)分析，初步研究天然竹管的抗沖擊性能。在d/h為0.49、1.03和5.98竹管上每隔h/4取點，依次記為點A、B、C；D、E、F；G、H、I，如圖3所示。

總的來看，當竹管入射端受到軸向沖擊時，產(chǎn)生入射波，并迅速沿竹管的徑向和軸向傳播，在竹管末端與透射桿接觸面發(fā)生反射，在入射波和反射波的共同作用下，竹管迅速變形，吸收沖擊能量，當達到竹管的極限強度時，出現(xiàn)劈裂破壞。

在同一幾何尺寸竹管中，沖擊載荷越大，位移、速度和加速度幅值越大，應(yīng)力波的傳播速度越快，同一節(jié)點處，曲線開始響應(yīng)時間越早。

具體來看，位移響應(yīng)中，在510μs后，位移幅值明顯增加，且沖擊載荷越大，位移響應(yīng)的幅值越大；對比3個不同位置處，靠近入射端面的管外壁節(jié)點處的位移響應(yīng)幅值最大，靠近透射端面的位移響應(yīng)幅值最小，說明沖擊載荷增大時，應(yīng)力波先傳到的地方，竹管發(fā)生壓潰變形，并向外擴展，吸收沖擊能量。速度響應(yīng)中，沖擊載荷越大，同一節(jié)點，同響應(yīng)時間下的速度幅值越大，增速越快。加速度響應(yīng)有明顯的振蕩現(xiàn)象，在響應(yīng)早期，沖擊載荷越大，加速度幅值越大，震蕩越劇烈，后期趨于平緩。

5" 軸向壓縮沖擊載荷下的不含節(jié)竹管的吸能特性

5.1" 能量計算

在SHPB試驗中，在高壓氣體作用下，子彈以一定的速度撞擊入射桿，產(chǎn)生的動能記為入射能WI，當入射桿中的應(yīng)力波傳入試樣時，一部分反射回入射桿，記為反射能WR，一部分傳遞到透射桿，記為透射能WT，還有部分能量由竹管變形所吸收，記為吸收能WS。吸收能與其自身質(zhì)量m的比值為比吸能Em。其中入射能、反射能、透射能、吸收能和比吸能的計算如式（11）～（16）［33］

WI=A0C0E0∫T0σ2I（t）dt（11）

WR=A0C0E0∫T0σ2R（t）dt（12）

WT=A0C0E0∫T0σ2T（t）dt（13）

WS=WI－WR－WT（14）

C0=E0ρ（15）

Em=WSm（16）

式（11）～（16）中：A0、C0和E0分別為壓桿的橫截面積、縱波波速和彈性模量；σI（t）、σR（t）和σT（t）分別為t時刻的入射、反射和透射應(yīng)力。

5.2" 能量耗散特性分析

在沖擊載荷作用下，SHPB系統(tǒng)中的入射能、反射能、透射能、吸收能和比吸能的計算結(jié)果和時程曲線，如表6和圖13所示。經(jīng)計算，d/h為0.49、1.03、5.98的竹管完整模型的質(zhì)量分別為15.86、7.77、1.46g。

由表6知，在相同的沖擊載荷下，入射能相等。在相同d/h情況下，隨沖擊載荷升高，入射能、反射能、吸收能和比吸能均增加，透射能變化不大。文中比吸能最大為6.64J/g，最小為0.60J/g。

由圖13知，隨沖擊載荷升高，入射能、反射能和吸收能均呈先上升趨勢，再逐漸趨于平緩。透射能隨時間的變化不大，近似于直線。相同響應(yīng)時間下，反射能遠大于透射能。原因為中空圓管試樣與壓桿間的接觸面積較小，當入射桿中的應(yīng)力波到達試樣時，大部分能量反射回入射桿，即入射能，只有小部分能量經(jīng)試樣傳播到透射桿，而且其中部分能量被竹管的變形所吸收，所以透射能變化不大，幾乎保持直線狀態(tài)。

如圖14所示，選取本研究SHPB實驗與文獻［34］落錘實驗中，高度較接近的竹管試樣進行對比，發(fā)現(xiàn)在SHPB實驗中竹管的峰值載荷Fmax較落錘實驗大，到達Fmax后曲線下降速度較落錘實驗快，F(xiàn)max所對應(yīng)的位移較落錘實驗小。所以竹管在這兩種加載條件下的載荷-位移曲線有明顯不同，吸能特性也會有所變化。

本研究用SHPB進行天然竹管的動態(tài)力學(xué)性能的研究，側(cè)重點在材料屬性上。而落錘的加載方式是垂直方向，側(cè)重于結(jié)構(gòu)屬性。SHPB實驗可滿足一維應(yīng)力波和應(yīng)力均勻化假定，實驗數(shù)據(jù)是在試樣整體相對均勻受力情況下測得的。

文獻［35］中試樣玻璃鋼纖維增強塑料（glass fiber reinforcement plastic，GFRP）薄壁管和天然竹管均為纖維增強復(fù)合材料，且形狀均為圓管狀。將實驗試樣進行對比，文獻［35］中應(yīng)變率為

160s-1，厚徑比δ/d為0.08時的GFRP薄壁圓管，與本研究中應(yīng)變率為150s-1，δ/d為0.09的不含節(jié)竹管試樣3-2的應(yīng)變率和厚徑比都比較接近，圖15和表7為GFRP薄壁圓管與天然不含節(jié)竹管的載荷-位移曲線和吸能特性對比。

由圖15知，在SHPB加載情形下，GFRP薄壁圓管的Fmax明顯遠大于不含節(jié)竹管，且大于落錘［36］和準靜態(tài)［37］加載時的Fmax。在落錘［36］和準靜態(tài)［37］加載條件下，GFRP薄壁管Fmax所對應(yīng)的位移接近，且均大于SHPB加載下的。由表7可知，在應(yīng)變率和δ/d都比較接近的情況下，GFRP薄壁圓管的比吸能與天然不含節(jié)竹管比較接近，說明天然不含節(jié)竹管具有一定的吸能防護能力；其吸收能和Fmax遠大于天然不含節(jié)竹管，說明天然不含節(jié)竹管并不適用于承受大沖擊能量的高速沖擊防護領(lǐng)域。

6" 結(jié)束語

本研究建立了不含節(jié)竹管受軸向沖擊載荷的有限元分析模型，通過與天然竹管的SHPB實驗結(jié)果對比，驗證了有限元分析模型的有效性。進而模擬研究不含節(jié)竹管的沖擊響應(yīng)與吸能特性。分析得出：天然不含節(jié)竹管在軸向沖擊壓縮載荷下的破壞模式隨竹管徑長比的變化有所不同，展現(xiàn)出了脆性材料沖擊破壞模式的特點，還有竹管會沿著軸向發(fā)生劈裂破壞的現(xiàn)象。建立的有限元分析模型模擬的沖擊破壞過程與實驗結(jié)果一致度較好；模擬的竹管在不同應(yīng)變率下的動態(tài)應(yīng)力-應(yīng)變關(guān)系，其中極限抗壓強度誤差最小為2.19%，最大為8.68%，與實驗結(jié)果吻合較好。模擬給出的竹管比吸能最大為6.64J/g，最小為0.60J/g。

不含節(jié)竹管外壁節(jié)點處的位移、速度和加速度響應(yīng)幅值均隨沖擊載荷的增加而增加，且加速度響應(yīng)出現(xiàn)明顯的振蕩現(xiàn)象。隨沖擊載荷升高，入射能、反射能、吸收能和比吸能均增加，而透射能隨響應(yīng)時間的變化不大。不含節(jié)竹管的峰值載荷遠低于GFRP薄壁圓管，但比吸能較為接近，說明竹管具有一定的防護能力，但并不適用于承受大沖擊能量的高速沖擊防護領(lǐng)域。

參考文獻：

［1］" 王志威，張曉偉，姚利宏，等.竹材宏觀微觀力學(xué)性能試驗與數(shù)值模擬［J］.林業(yè)工程學(xué)報，2022，7（1）：31-37.

WANG Zhiwei，ZHANG Xiaowei，YAO Lihong，et al.Experimental study and numerical simulation on the macro and micro mechanical properties of bamboo［J］.Journal of forestry engineering，2022，7（1）：31-37（in Chinese）.

［2］" HU D Y，SONG B，DANG L W，et al.Effect of strain rate on mechanical properties of the bamboo material under quasi-static and dynamic loading condition［J］.Composite structures，2018，200：635-646.

［3］" 任海青，張東升，潘雁紅.竹材抗壓動態(tài)破壞過程分析［J］.南京林業(yè)大學(xué)學(xué)報（自然科學(xué)版），2007，31（2）：47-50.

REN Haiqing，ZHANG Dongsheng，PAN Yanhong.Dynamic compressive mechanical behavior of bamboo［J］.Journal of Nanjing Forestry University（natural sciences edition），2007，31（2）：47-50（in Chinese）.

［4］" LAKKAD S C，PATEL J M.Mechanical properties of bamboo，a natural composite［J］.Fibre science and technology，1981，14（4）：319-322.

［5］" 何子奇，周緒紅，肖灑，等.毛竹受壓構(gòu)件性能試驗研究［J］.建筑結(jié)構(gòu)學(xué)報，2018，39（增刊2）：233-241.

HE Ziqi，ZHOU Xuhong，XIAO Sa，et al.Experimental investigation on bamboo columns under axial compression［J］.Journal of building structures，2018，39（S2）：233-241（in Chinese）.

［6］" GYANSAH L，AKINWONMI A S，AFFAM M.The fracture behaviour of fresh bamboo under uniaxial compressive loading condition［J］.Research journal of applied sciences，engineering and technology，2010，2（8）：720-726.

［7］" ZOU M，WEI C G，LI J Q，et al.The energy absorption of bamboo under dynamic axial loading［J］.Thin-walled structures，2015，95：255-261.

［8］" ZUHRI M Y M，LIAO Y，WANG Q Y，et al.The energy absorbing properties of bamboo-based structures［J］.Journal of sandwich structures amp; materials，2019，21（3）：1032-1054.

［9］" 趙章榮，傅萬四，劉星雨，等.結(jié)構(gòu)用竹材機械性能測試及數(shù)值模擬研究［J］.木材加工機械，2014，25（1）：11-13.

ZHAO Zhangrong，F(xiàn)U Wansi，LIU Xingyu，et al.Study on bamboo culm used for structure mechanical performance tests and numerical simulation［J］.Wood processing machinery，2014，25（1）：11-13（in Chinese）.

［10］ZHAO Z R，F(xiàn)U W S，HAN W，et al.Study on bamboo culm used for structure axial compression performance numerical simulation［J］.Advanced materials research，2014，842：13-17.

［11］FU W S，ZHAO Z R，HAN W，et al.Research on finite element model for parallel to bamboo culms axial shear［J］.Applied mechanics and materials，2014，477-478：986-989.

［12］AKINLABI E T，ANANE-FENIN K，AKWADA D R.Bamboo：the multipurpose plant［M］.Cham：Springer International Publishing，2017：225-226.

［13］CUI J H，QIN Z，MASIC A，et al.Multiscale structural insights of load bearing bamboo：a computational modeling approach［J］.Journal of the mechanical behavior of biomedical materials，2020，107：103743.

［14］ZOU M，XU S C，WEI C G，et al.A bionic method for the crashworthiness design of thin-walled structures inspired by bamboo［J］.Thin-walled structures，2016，101：222-230.

［15］許述財，鄒猛，魏燦剛，等.仿竹結(jié)構(gòu)薄壁管的軸向耐撞性分析及優(yōu)化［J］.清華大學(xué)學(xué)報（自然科學(xué)版），2014，54（3）：299-304.

XU Shucai，ZOU Meng，WEI Cangang，et al.Axial crash worthiness analysis and optimization of a bionic thin-walled tube based on bamboo structure［J］.Journal of Tsinghua University（science and technology），2014，54（3）：299-304（in Chinese）.

［16］SONG J F，XU S C，WANG H X，et al.Bionic design and multi-objective optimization for variable wall thickness tube inspired bamboo structures［J］.Thin-walled structures，2018，125：76-88.

［17］FU J，LIU Q，LIUFU K M，et al.Design of bionic-bamboo thin-walled structures for energy absorption［J］.Thin-walled structures，2019，135：400-413.

［18］CHEN B C，ZOU M，LIU G M，et al.Experimental study on energy absorption of bionic tubes inspired by bamboo structures under axial crushing［J］.International journal of impact engineering，2018，115：48-57.

［19］王禮立.應(yīng)力波基礎(chǔ)［M］.2版.北京：國防工業(yè)出版社，2005.

［20］ZHANG S R，WANG X H，WANG C，et al.Compressive behavior and constitutive model for roller compacted concrete under impact loading：considering vertical stratification［J］.Construction and building materials，2017，151：428-440.

［21］曾其蘊，李世紅，鮑賢镕.竹節(jié)對竹材力學(xué)強度影響的研究［J］.林業(yè)科學(xué)，1992，28（3）：247-252.

ZENG Qiyun，LI Shihong，BAO Xianrong.Effect of bamboo nodal on mechanical properties of bamboo wood［J］.Scientia silvae Sinicae，1992，28（3）：247-252（in Chinese）.

［22］SAKARAY H，TOGATI N V V K，REDDY I V R.Investigation on properties of bamboo as reinforcing material in concrete［J］.International journal of engineering research and applications，2012，2（1）：77-83.

［23］邵卓平，黃盛霞，吳福社，等.毛竹節(jié)間材與節(jié)部材的構(gòu)造與強度差異研究［J］.竹子研究匯刊，2008，27（2）：48-52.

SHAO Zhuoping，HUANG Shengxia，WU Fushe，et al.A study on the difference of structure and strength between internodes and nodes of moso bamboo［J］.Journal of bamboo research，2008，27（2）：48-52（in Chinese）.

［24］CHEN K，LONG L C.Analysis of the effects of fiber gradient distribution on the mechanical properties of Moso bamboo［J］.Applied mechanics and materials，2014，590：13-18.

［25］LONG L C，WANG Z K，CHEN K.Analysis of the hollow structure with functionally gradient materials of Moso bamboo［J］.Journal of wood science，2015，61（6）：569-577.

［26］SHARMA A，MISHRA R，JAIN S，et al.Deformation behavior of single and multi-layered materials under impact loading［J］.Thin-walled structures，2018，126：193-204.

［27］MIAO H R，WU Z Y，YING Z P，et al.The numerical and experimental investigation on low-velocity impact response of composite panels：effect of fabric architecture［J］.Composite structures，2019，227：111343.

［28］荊臻，田常錄，吳健，等.帶孔復(fù)合材料層板動態(tài)拉伸破壞的應(yīng)變率效應(yīng)［J］.南京航空航天大學(xué)學(xué)報，2019，51（1）：55-62.

JING Zhen，TIAN Changlu，WU Jian，et al.Strain rate effects on dynamic tensile failure of composite laminated plates with holes［J］.Journal of Nanjing University of Aeronautics amp; Astronautics，2019，51（1）：55-62（in Chinese）.

［29］葉輝，朱艷榮，蒲永鋒.纖維增強復(fù)合材料應(yīng)變率效應(yīng)的數(shù)值仿真［J］.吉林大學(xué)學(xué)報（工學(xué)版），2019，49（5）：1622-1629.

YE Hui，ZHU Yanrong，PU Yongfeng.Numerical simulation of strain rate effect of fiber reinforced composites［J］.Journal of Jilin University（engineering and technology edition），2019，49（5）：1622-1629（in Chinese）.

［30］陳博斐，邱欣，陳江瑛.混凝土SHPB束桿裝置沖擊實驗的數(shù)值仿真［J］.應(yīng)用力學(xué)學(xué)報，2017，34（1）：125-130.

CHEN Bofei，QIU Xin，CHEN Jiangying.Numerical simulation of concrete impact test using split Hopkinson pressure bundled bars device［J］.Chinese journal of applied mechanics，2017，34（1）：125-130（in Chinese）.

［31］陳為農(nóng)，宋博.分離式霍普金森（考爾斯基）桿：設(shè)計、試驗和應(yīng)用［M］.北京：國防工業(yè)出版社，2018.

［32］趙章榮，楊光，傅萬四，等.整竹開裂性能數(shù)值模擬研究［J］.木材加工機械，2015，26（6）：4-6.

ZHAO Zhangrong，YANG Guang，F(xiàn)U Wansi，et al.Study on bamboo culm cracking performancenumerical simulation［J］.Wood processing machinery，2015，26（6）：4-6（in Chinese）.

［33］ZOU Y C，XIONG C，YIN J H，et al.Experimental and modeling studies of stress wave propagation and energy dissipation mechanism in layered composite structures［J］.Shock and vibration，2021，2021：9912709.

［34］UMER R，BALAWI S，RAJA P，et al.The energy-absorbing characteristics of polymer foams reinforced with bamboo tubes［J］.Journal of sandwich structures amp; materials，2014，16（1）：108-122.

［35］李善恩，繆馥星，周風(fēng)華，等.玻璃鋼纖維增強塑料薄壁管抗沖擊性能的實驗研究［J］.復(fù)合材料學(xué)報，2018，35（12）：3324-3330.

LI Shanen，MIAO Fuxing，ZHOU Fenghua，et al.Experimental study on impact performance of glass fiber reinforced plastic thin-walled tubes［J］.Acta materiae compositae Sinica，2018，35（12）：3324-3330（in Chinese）.

［36］PALANIVELU S，VAN PAEPEGEM W，DEGRIECK J，et al.Experimental study on the axial crushing behaviour of pultruded composite tubes［J］.Polymer testing，2010，29（2）：224-234.

［37］MAHDI E，SEBAEY T A.An experimental investigation into crushing behavior of radially stiffened GFRP composite tubes［J］.Thin-walled structures，2014，76：8-13.

（編輯" 呂茵）