郭瑞奇，任輝啟，龍志林，吳祥云，姜錫權(quán)

（1. 湘潭大學(xué)土木工程與力學(xué)學(xué)院，湖南 湘潭 411105；2. 軍事科學(xué)院國防工程研究院，河南 洛陽 471023；3. 合肥姜水動態(tài)力學(xué)實驗技術(shù)有限公司，安徽 合肥 230031）

混凝土作為房屋、道路、橋梁等工程結(jié)構(gòu)中應(yīng)用最為廣泛的建筑材料，對其各項力學(xué)性能進(jìn)行充分研究是保證結(jié)構(gòu)設(shè)計安全可靠的基礎(chǔ)?？紤]到爆炸、沖擊等突發(fā)事件，一些重要的建筑設(shè)施在設(shè)計時不僅要考慮其承擔(dān)的常規(guī)載荷作用，同時還要考慮其所承擔(dān)的強(qiáng)動載作用。因此，需要對高應(yīng)變率下混凝土材料的各項動態(tài)力學(xué)性能進(jìn)行充分研究。

從目前對混凝土材料的動態(tài)力學(xué)性能研究狀況來看，大多數(shù)仍以動態(tài)壓縮實驗為主[1-2]。然而，混凝土是一種拉壓強(qiáng)度嚴(yán)重不對稱的材料，在靜載作用下其抗拉強(qiáng)度約比抗壓強(qiáng)度小一個量級。而且混凝土在壓力作用下的破碎和拉伸載荷作用下的斷裂是其兩種主要的破壞形式[3]，所以很有必要對混凝土的動態(tài)拉伸力學(xué)性能進(jìn)行研究。

由于分離式Hopkinson 壓桿(split Hopkinson pressure bar，SHTB)實驗裝置結(jié)構(gòu)簡單，發(fā)展較為成熟，因此對混凝土動態(tài)拉伸力學(xué)性能的研究多以使用SHPB 裝置進(jìn)行劈裂[4-6]和層裂[7-10]實驗的形式進(jìn)行，如圖1 所示。其中用于間接測量混凝土材料拉伸性能的巴西圓盤試件在加載過程中處于一種復(fù)雜的應(yīng)力狀態(tài)，且拉應(yīng)力場并不均勻，不能反映材料一維應(yīng)力狀態(tài)的動態(tài)拉伸強(qiáng)度[11]。層裂是混凝土材料的一種典型的沖擊拉伸破壞形式，且試件在破壞過程中處于一維應(yīng)力狀態(tài)。但對于混凝土材料而言，初始加載波有可能對試件造成一定的壓縮損傷，此時材料所表現(xiàn)出的層裂強(qiáng)度已不再是初始材料的特性，而是受壓縮波損傷弱化后的材料特性[12]。

圖1 使用Hopkinson 壓桿設(shè)備對混凝土材料進(jìn)行動態(tài)拉伸實驗[4, 10]Fig.1 Dynamic tensile experiments of concrete materials with Hopkinson pressure bar apparatus[4, 10]

為了研究混凝土材料在一維應(yīng)力狀態(tài)下的動態(tài)拉伸力學(xué)性能，目前已研制出了不同的大直徑SHTB實驗裝置，例如張凱等[13]研制的 ? 75 mm 桿-桿型臥式SHTB 實驗裝置和Levi-Hevroni 等[14]研制的使用落錘加載的 ? 100 mm 立式SHTB 裝置。這種使用管狀子彈/落錘撞擊法蘭盤產(chǎn)生拉伸波的方式存在著套管的對心同軸問題，而且法蘭盤與入射桿的連接會對拉伸波產(chǎn)生影響，不易控制拉伸波形。針對以上問題，姜錫權(quán)等[15]發(fā)明了一種新型霍普金森拉桿裝置，這種新型拉桿使用實心子彈撞擊空心入射管封頭部位，從而產(chǎn)生拉伸波。本文中通過數(shù)值模擬方法研究該新型拉桿中的封頭、銷釘結(jié)構(gòu)對應(yīng)力波傳播的影響。最后，根據(jù)計算結(jié)果對設(shè)備的局部構(gòu)件進(jìn)行改進(jìn)，并建立圓環(huán)狀的細(xì)觀骨料模型開展混凝土材料的動態(tài)直拉模擬研究。

1 實驗裝置介紹

新型霍普金森拉桿裝置如圖2 所示，裝置主要由吸能裝置、圓柱狀實心子彈、倒錐形封頭、空心入射管、實心入射桿、銷釘以及實心透射桿組成。其中空心入射管內(nèi)徑為100 mm，外徑為145 mm，相對應(yīng)的子彈和實心入射桿直徑為100 mm，倒錐形封頭小端直徑為91 mm。

圖2 使用實心子彈撞擊空心入射管封頭部位來產(chǎn)生拉伸波的大直徑SHTB 設(shè)備Fig.2 Large diameter SHTB apparatus utilizing a solid bullet to strike end cap for generation of tensile stress wave

實驗時將試件通過膠粘的方式與入射桿和透射桿連接，然后壓縮空氣發(fā)射系統(tǒng)驅(qū)動撞擊桿（子彈）撞擊倒錐形封頭帶動空心入射管產(chǎn)生一個拉伸脈沖，拉伸脈沖通過銷釘連接傳播至實心入射桿。試件在該拉伸入射脈沖的作用下高速變形直至被拉斷，與此同時向入射桿回傳反射脈沖并向透射桿傳播透射脈沖。緩沖裝置起到吸能的作用，使空心入射管停止運(yùn)動。實驗裝置如圖3所示。

基于一維應(yīng)力波假定，通過應(yīng)變片所測量到入射桿上的入射應(yīng)變波 εI(t) 、反射應(yīng)變波 εR(t) 和透射桿上的透射應(yīng)變波 εT(t) 可求得試件的平均應(yīng)力 σs(t) 、應(yīng)變率 ε˙s(t) 和平均應(yīng)變 εs(t) ，即三波法公式[16]：

式中： A0為桿的橫截面積，As和ls分別為試件的橫截面積和長度，E0和c0分別為入射桿的彈性模量和波速。

2 大直徑SHTB 裝置的數(shù)值研究

首先驗證使用四面體單元對桿件進(jìn)行網(wǎng)格劃分是否會對應(yīng)力波的傳播造成影響，然后建立大直徑SHTB 裝置的有限元模型并進(jìn)行加載計算，對裝置中封頭和銷釘?shù)挠嬎憬Y(jié)果以及空心入射管中的拉伸波形進(jìn)行分析。

2.1 計算前的模型驗證

由于SHPB 裝置構(gòu)造簡單，入射桿和透射桿多為標(biāo)準(zhǔn)的實心長圓柱體，在數(shù)值模擬中通常使用六面體單元(LS-DYNA 中對應(yīng)的Solid164 單元)對這種規(guī)則的幾何體進(jìn)行映射網(wǎng)格劃分[17]。但對于圖2～3 所示的SHTB 裝置而言，其結(jié)構(gòu)相對復(fù)雜，入射管封頭附近的排氣孔以及銷釘部位分別需要開槽和開孔，因此使用四面體單元(LS-DYNA 中對應(yīng)的Solid168 單元)對這種結(jié)構(gòu)進(jìn)行自由網(wǎng)格劃分。

一般在正式實驗前都會進(jìn)行空打測試來進(jìn)行實驗標(biāo)定，也就是移除試件，將入射桿和透射桿對齊，使用子彈直接撞擊入射桿。如果所產(chǎn)生的入射波能順利傳播至透射桿，沒有反射波或反射波微弱到可以忽略不計的地步，則認(rèn)為入射桿和透射桿已經(jīng)調(diào)平可以開始正式實驗。

在正式進(jìn)行數(shù)值計算前，首先使用四面體單元對實心圓柱桿件進(jìn)行網(wǎng)格劃分，并使用背景網(wǎng)格映射法建立混凝土細(xì)觀骨料模型，如圖4 所示，將各個材料均設(shè)置為同類線彈性模型，與桿件材料保持一致。所有的接觸均設(shè)置為理想固連接觸，即不考慮連接失效問題。然后在桿端施加峰值100 MPa，持續(xù)時間480 μs 的半正弦拉伸入射波進(jìn)行加載。這樣做的目的，一方面是為了檢驗使用四面體單元劃分桿件是否會對應(yīng)力波的傳播造成影響，另一方面是為了以此來檢驗?zāi)Ｐ椭兴O(shè)置的接觸類型，網(wǎng)格劃分質(zhì)量以及所建立的混凝土細(xì)觀骨料模型的有效性。

應(yīng)力波穿過試件前后的應(yīng)力云圖如圖5 所示，可以看到應(yīng)力波穿過試件前后其波陣面一直保持為平截面，試件的網(wǎng)格劃分和接觸設(shè)置對應(yīng)力波的傳播沒有影響。入射桿和透射桿中的波形如圖6 所示，可以看到反射波微弱到可以忽略不計的地步，透射波與入射波相比，其峰值和持續(xù)時間保持不變。因此，可以認(rèn)為本文中的接觸設(shè)置、網(wǎng)格劃分質(zhì)量及有限元模型是合理可靠的。

圖5 入射桿和透射桿的軸向應(yīng)力云圖Fig.5 Axial stress distribution in incident and transmitted bars

圖6 入射桿和透射桿上的應(yīng)力波形Fig.6 Stress waveforms in incident and transmitted bar

2.2 大直徑SHTB 裝置有限元模型和計算分析

在驗證了四面體單元分析應(yīng)力波傳播問題的有效性以后，下面建立與圖3 對應(yīng)的大直徑SHTB 設(shè)備各部位的有限元模型，如圖7 所示。共有4 個構(gòu)件，均設(shè)置為線彈性模型，取彈性模量E=200 GPa，密度 ρ0=7.8×103kg/m3，泊松比v=0.3。對于LS-DYNA 顯示動力學(xué)分析中的接觸碰撞問題，通常在接觸物體之間使用ASTS 自動面面接觸[18-19]。因此對于本文中的大直徑SHTB 設(shè)備，在子彈與封頭、封頭與入射管、入射管與銷釘、入射桿與銷釘之間設(shè)置4 對自動面面接觸。

圖7 大直徑SHTB 設(shè)備各部位的有限元模型Fig.7 Finite element models for several parts of the large-diameter SHTB apparatus

倒錐形封頭在子彈以10 m/s 的速度沖擊下與空心入射管碰撞擠壓，作用過程中會使得封頭的一小部分?jǐn)D出入射管(圖8(a))。但這種情況在彈性變形范圍內(nèi)是可逆的，隨著子彈中的卸載波作用完畢，子彈與封頭碰撞的結(jié)束，封頭與入射管的位移逐漸趨于一致(圖8(c))，并最終會恢復(fù)至碰撞前的狀態(tài)(圖8(b))。

圖8 子彈撞擊倒錐形封頭計算結(jié)果Fig.8 Calculation results for the bullet striking the reverse taper-shaped end cap

為了求出封頭在撞擊面上的質(zhì)點應(yīng)力，從而預(yù)估出入射管中的拉伸應(yīng)力波幅值，將子彈與封頭的碰撞過程看作兩根截面尺寸相同的彈性桿共軸撞擊過程。其密度 ρ0均為7 800 kg/m3，波速c0為5 190 m/s，子彈撞擊速度 v2為10 m/s，封頭初始速度 v1為零，則由一維桿中應(yīng)力波的初等理論[20]：

可求出碰撞后封頭的撞擊面上質(zhì)點的應(yīng)力為202.41 MPa，因此可認(rèn)為撞擊后入射管中拉伸波應(yīng)力幅值約為200 MPa。對于大直徑SHTB設(shè)備而言，幅值約為200 MPa 的拉伸波形能滿足大多數(shù)巖土類脆性材料的動態(tài)拉伸實驗，此時銷釘受力情況如圖9～10 所示。銷釘在空心入射管和實心入射桿的作用下處于一種雙剪切應(yīng)力狀態(tài)，可以看到，在10 m/s 的撞擊速度下，銷釘?shù)淖畲笾鲬?yīng)力不超過130 MPa，不超過鋼材料的屈服強(qiáng)度。銷釘在距離中間部位5 cm 以內(nèi)受力較為均勻，從7.5 cm 處受力開始逐漸下降。

l

圖9 銷釘最大主應(yīng)力云圖Fig.9 Maximum principal stress nephograms of dowels

圖10 銷釘受力情況Fig.10 Stresses born by dowels

2.3 大直徑SHTB 裝置應(yīng)力波形分析

子彈以10 m/s 的速度撞擊倒錐形封頭，在空心入射管和實心入射桿中產(chǎn)生的拉伸波形如圖11 所示，可以看到入射管中的拉伸波形為平均應(yīng)力幅值185 MPa、升時130 μs 且?guī)в懈哳l振蕩的梯形波。考慮到封頭部位與入射管的碰撞擠壓過程中會損失一定的能量，所以與2.2 節(jié)中計算得到的封頭撞擊面質(zhì)點應(yīng)力相比，入射管中實際拉伸波形的應(yīng)力幅值略低于200 MPa。通過銷釘后，實心入射桿中的拉伸波形相對光滑，平均應(yīng)力幅值沒有發(fā)生明顯衰減，升時延長35 μs。圖12 所示的應(yīng)力波通過銷釘前后的應(yīng)力云圖也表明，通過銷釘傳入到入射桿中的拉伸波在桿截面分布較為均勻，波陣面保持為平截面，沒有發(fā)生畸變。

圖11 大直徑SHTB 設(shè)備中的拉伸波形Fig.11 Tensile stress waveforms in the large diameter SHTB apparatus

圖12 應(yīng)力波通過銷釘前后的應(yīng)力云圖Fig.12 Stress nephograms of the wave propagation process around the dowels part

自20 世紀(jì)70 年代以來，有研究者開始在入射桿前端添加一個波阻抗相對較低的波形整形器來改善入射波形[21]，這樣做一方面可以降低大直徑桿中應(yīng)力波形的彌散震蕩，另一方面也有利于實現(xiàn)恒應(yīng)變率加載，及早實現(xiàn)試件中應(yīng)力分布的均勻化。到目前為止，紙片、橡膠、塑料以及一些金屬均作為整形器的材料被應(yīng)用于SHPB 實驗當(dāng)中，本文中采用LS-DYNA 中的Johnson-Cook 模型來模擬紫銅波形整形器，材料參數(shù)參照文獻(xiàn)[22]，研究其對于大直徑SHTB 設(shè)備的整形效果，為以后的實際實驗提供指導(dǎo)。其中整形器直徑為15 mm，厚度為2 mm，其有限元模型如圖13 所示。在整形器與封頭和子彈之間分別增設(shè)兩對自動面面接觸，撞擊速度仍設(shè)置為10 m/s，計算得到的波形如圖14 所示。與圖11 中子彈直接撞擊產(chǎn)生的拉伸波形相比，使用紫銅波形整形器以后在入射管中可以得到升時延長且相對光滑的半正弦拉伸入射波。

圖13 紫銅整形器有限元模型剖視圖Fig.13 Sectional drawing of the finite element model for the red copper pulse shaper

圖14 整形后入射管和入射桿中的拉伸波形Fig.14 Shaped tensile waveforms in incident tube and incident bar

3 大直徑SHTB 設(shè)備改進(jìn)優(yōu)化及混凝土細(xì)觀骨料模型動態(tài)直拉數(shù)值模擬

3.1 大直徑SHTB 裝置的改進(jìn)優(yōu)化

Zhang 等[23]使用SHPB 實驗裝置對混凝土類材料進(jìn)行動態(tài)壓縮測試時，分別制作了實心試件和空心試件進(jìn)行對比，結(jié)果顯示使用空心試件測得的動態(tài)強(qiáng)度增強(qiáng)因子要低于實心試件，表明由慣性效應(yīng)引起的徑向約束會對混凝土類材料的動態(tài)強(qiáng)度測試結(jié)果產(chǎn)生重要影響。Li 等[24]通過進(jìn)一步的數(shù)值模擬研究發(fā)現(xiàn)，采用空心試件可以有效減輕試件橫向慣性效應(yīng)的影響，所測得的動態(tài)強(qiáng)度增強(qiáng)因子也相對真實。

最近的研究結(jié)果[25]表明，在使用大直徑SHTB 裝置對實心和空心混凝土試件進(jìn)行測試時，混凝土材料的動態(tài)拉伸強(qiáng)度也同樣會被慣性效應(yīng)所影響?；谝陨涎芯?，同時考慮到2.3 節(jié)中使用波形整形技術(shù)同樣可以在入射管中得到光滑的半正弦入射波，可以將銷釘連接的實心入射桿所簡化，并且把實心透射桿也改造成空心圓管，直接對圓環(huán)狀的混凝土試件進(jìn)行拉伸實驗。入射管和透射管與試件的連接端設(shè)置有內(nèi)螺紋，可與配套的夾具相連接，這樣也同時兼顧了對膠粘的實心圓柱狀混凝土試件或掛接的啞鈴狀混凝土試件進(jìn)行動態(tài)拉伸實驗。以啞鈴狀試件的掛接形式為例，圖15(a)所示為局部改進(jìn)后的大直徑SHTB 設(shè)備示意圖，圖15(b)和圖15(c)則分別為帶有內(nèi)螺紋的入射管和透射管端面的實物圖。

圖15 局部改進(jìn)的大直徑SHTB 設(shè)備Fig.15 Partly improved large diameter SHTB apparatus

下面，針對改進(jìn)后的SHTB 裝置，進(jìn)行混凝土動態(tài)拉伸數(shù)值模擬研究，并與文獻(xiàn)[25]中的實驗結(jié)果進(jìn)行對比分析。一方面為了驗證采用空心入射管和透射管這種結(jié)構(gòu)形式的有效性，為以后的實際實驗提供參考和幫助，另一方面，也為混凝土細(xì)觀骨料模型的動態(tài)拉伸數(shù)值模擬研究提供思路。

3.2 圓環(huán)狀混凝土細(xì)觀骨料模型的建立及材料模型介紹

混凝土是一種非均勻材料，其內(nèi)部含有大量隨機(jī)分布的粗骨料顆粒，采用數(shù)值模擬方法對混凝土材料進(jìn)行研究時若建立均勻的有限元模型則與實際情況不符。近年來，已有諸多學(xué)者將砂漿和粗骨料分開建模以從細(xì)觀層次上研究混凝土的力學(xué)性能，例如：金瀏等[26]對單邊缺口的混凝土試件及L 形試件在不同加載速率下的動態(tài)拉伸破壞模式進(jìn)行了細(xì)觀數(shù)值研究；吳成等[27]通過對混凝土細(xì)觀骨料模型侵徹的模擬得出，影響靶板抗侵徹能力的主要因素是砂漿種類、粗骨料種類和粗骨料體積分?jǐn)?shù)；鄧勇軍等[28]基于混凝土細(xì)觀骨料模型研究了細(xì)觀因素（骨料隨機(jī)分布、骨料強(qiáng)度、砂漿強(qiáng)度等）對剛性彈體正侵徹偏轉(zhuǎn)角度的影響。

前期的工作中，使用混合實現(xiàn)方法能夠快速有效地建立混凝土的細(xì)觀骨料模型[29]，并能夠適用于SHPB 數(shù)值模擬[17]。但是對于形狀較為復(fù)雜的混凝土試件，其對應(yīng)的邊界條件也較為復(fù)雜，因此為了便于建立圓環(huán)狀的混凝土細(xì)觀骨料模型，并減少單元數(shù)量，提高計算效率，本文中將使用背景網(wǎng)格映射法進(jìn)行混凝土細(xì)觀建模。首先，建立所需要圓環(huán)狀試件的幾何模型，并對試件進(jìn)行均勻的網(wǎng)格劃分，提取相應(yīng)的節(jié)點信息和單元信息予以保存。然后，需要生成骨料和砂漿的幾何信息，在這里，將骨料尺寸設(shè)定在5～20 mm 范圍內(nèi)。定義3 個位于(0,1)之間的獨(dú)立隨機(jī)數(shù) nrdm1、 nrdm2、 nrdm3，在三維笛卡爾坐標(biāo)系中，隨機(jī)球形骨料的球心坐標(biāo)可以表示為：

式中： (xi,yi,zi) 為正在投放的第i 個骨料的球心坐標(biāo)， (xn,yn,zn) 為已經(jīng)生成的前i?1 個骨料的球心坐標(biāo),如果不滿足式(5)要求則需要重新生成坐標(biāo)。

設(shè)一共生成了 m 顆骨料，現(xiàn)在有了每個骨料的球心坐標(biāo) (xn,yn,zn) 和半徑 Rn，通過節(jié)點信息可以求得每個單元的中心坐標(biāo) (xe,ye,ze) ，通過式(6)對每個單元進(jìn)行材料屬性判別。

也就是當(dāng)單元與骨料球心之間的距離小于骨料半徑時，則該單元為骨料單元。當(dāng)進(jìn)行了 m 次判定后，單元與 m 個骨料的距離均大于骨料半徑時，該單元為砂漿單元，最后生成的細(xì)觀模型如圖16(b)～(c)所示，其對應(yīng)的幾何模型如圖16(a)所示。

從Chen 等[30]的混凝土細(xì)觀層裂模擬和Xu 等[31]的混凝土細(xì)觀劈裂模擬來看，使用的都是Malvar等[32-33]提出的CONCRETE_DAMAGE_REL3(Mat_072R3)材料模型，也常常被稱為K&C 模型。該模型使用了3 個失效面來描述混凝土材料的塑性性能，考慮了損傷破壞效應(yīng)，并可以通過定義拉伸動態(tài)強(qiáng)度增強(qiáng)因子曲線和壓縮動態(tài)強(qiáng)度增強(qiáng)因子曲線來定義拉伸、壓縮應(yīng)變率效應(yīng)。本文中主要參照Chen 等[30]的混凝土細(xì)觀層裂模擬工作來分別定義砂漿和骨料的動態(tài)強(qiáng)度因子曲線。

3.3 結(jié)果分析

圖16 使用背景網(wǎng)格映射法建立的混凝土細(xì)觀骨料模型Fig.16 Mesoscale concrete models established by the background grid mapping method

圖17 數(shù)值模擬應(yīng)力波形與實驗結(jié)果[25]對比Fig.17 Comparison between simulated and experimental[25] waveforms

在入射桿端直接添加幅值17 MPa 的半正弦入射波進(jìn)行加載，所得到的原始波形如圖17(a)所示。分別提取入射波、反射波和透射波，并將入射波與反射波相疊加與透射波進(jìn)行對比，以驗證試件的應(yīng)力均勻性，結(jié)果如圖17(b)所示，實驗測量的波形[25]如圖17(c)所示?？梢钥吹?，在有足夠升時的半正弦波作用下，數(shù)值模擬和實驗結(jié)果中試件均能滿足應(yīng)力均勻性要求。在滿足試件應(yīng)力均勻性的情況下，可以從透射波峰值看出試件的動態(tài)拉伸強(qiáng)度。從圖17(b)和圖17(c)的對比中可以看到，數(shù)值模擬的透射波峰值要高于實驗的結(jié)果，這是由于在數(shù)值模擬中試件模型的建立考慮了粗骨料的影響，實驗所采用的為不考慮粗骨料的水泥砂漿試件。粗骨料通常作為混凝土材料中的增強(qiáng)相，其抗拉強(qiáng)度要高于水泥砂漿的強(qiáng)度，從而提高了試件的整體抗拉強(qiáng)度。由于實驗中的砂漿試件強(qiáng)度較低，在拉伸入射波的作用下迅速破壞，試件的波阻抗也隨之降低，這會導(dǎo)致透射波的持續(xù)時間也相對較短。除此之外，從反射波和透射波的形態(tài)來看，數(shù)值模擬結(jié)果與實驗結(jié)果符合良好。

圓環(huán)狀混凝土試件數(shù)值模擬的破壞結(jié)果與實驗的破壞結(jié)果[25]對比如圖18 所示。結(jié)合圖17 中的波形分析可以知道，由于所建立的混凝土細(xì)觀骨料模型動態(tài)抗拉強(qiáng)度較高，所以與實驗結(jié)果相比，數(shù)值模擬中圓環(huán)狀混凝土試件的破壞程度相對較低，表面的裂紋相對零散，并未完全融匯貫通，試件也沒有被完全拉斷。但從破壞的形態(tài)和裂紋位置來看，在試件受力均勻的情況下，裂紋主要分布在試件的中間部位，且主裂紋呈“Z”字樣或曲線樣，與實驗結(jié)果[25]相同。

圖18 圓環(huán)狀混凝土試件數(shù)值模擬破壞結(jié)果與實驗結(jié)果[25]對比Fig.18 Comparison between simulated and experimental[25] phenomena

另外，由于粗骨料拉伸強(qiáng)度大于砂漿基體，在數(shù)值模擬中裂紋都是沿著粗骨料邊緣發(fā)展的。從圖18(b)～(c)可以看到，混凝土試件中大部分的粗骨料都處于彈性狀態(tài)，幾乎未發(fā)生塑性變形，因此若要充分發(fā)揮骨料的抗拉強(qiáng)度，還需要進(jìn)一步增強(qiáng)砂漿基體本身的強(qiáng)度。

從以上分析可知，采用空心圓管式的大直徑SHTB 設(shè)備對圓環(huán)狀混凝土試件進(jìn)行直接拉伸是一種行之有效的實驗方法，所建立的圓環(huán)狀細(xì)觀模型能夠定性地描述混凝土材料的動態(tài)拉伸力學(xué)行為。

4 結(jié)論與展望

運(yùn)用數(shù)值模擬方法對新型分離式霍普金森拉桿裝置進(jìn)行了計算和力學(xué)分析，提出了優(yōu)化改進(jìn)措施，并建立了圓環(huán)狀的細(xì)觀骨料模型開展動態(tài)拉伸模擬，得到以下結(jié)論。

（1）通過實心子彈撞擊倒錐形封頭的方式帶動空心入射管運(yùn)動能夠有效產(chǎn)生拉伸入射波，這種結(jié)構(gòu)形式下子彈以10 m/s 的速度撞擊能產(chǎn)生幅值約185 MPa 的拉伸入射波，足以滿足大多數(shù)混凝土、巖石類材料的動態(tài)拉伸實驗，因此子彈速度的主要適用范圍為0～10 m/s。封頭部位的變形和銷釘?shù)氖芰Ψ秶趶椥苑秶鷥?nèi)，使得裝置滿足多次重復(fù)使用的要求。

（2）通過在封頭部位的撞擊面粘貼波形整形器可以在入射管中產(chǎn)生理想的拉伸入射波，因此，即使采用銷釘連接空心入射管和實心入射桿不會使拉伸入射波產(chǎn)生明顯畸變，考慮到試件的橫向慣性效應(yīng)，同時便于配套夾具的連接，可以將銷釘連接的實心入射桿所改進(jìn)，并且把實心透射桿也改造成空心圓管，并在與試件的連接端面設(shè)置內(nèi)螺紋以兼顧對膠粘的實心圓柱狀混凝土試件或掛接的啞鈴狀混凝土試件進(jìn)行動態(tài)拉伸實驗。

（3）使用背景網(wǎng)格映射法建立了圓環(huán)狀混凝土細(xì)觀骨料模型，通過對其動態(tài)拉伸的模擬結(jié)果與實驗結(jié)果的對比分析驗證了改進(jìn)后的大直徑SHTB 裝置結(jié)構(gòu)的有效性，同時為混凝土細(xì)觀骨料模型的動態(tài)拉伸模擬提供了思路。

在驗證了本文所建立的混凝土細(xì)觀骨料模型用于SHTB 動態(tài)拉伸模擬的有效性之后，可以此為基礎(chǔ)，研究骨料含量、骨料的隨機(jī)分布以及骨料與砂漿的基礎(chǔ)物理參數(shù)等細(xì)觀因素對整個混凝土試件動態(tài)拉伸力學(xué)性能的影響，這將是我們今后進(jìn)一步研究的問題。另外，對于局部改進(jìn)后的大直徑SHTB 設(shè)備，如何針對混凝土、巖石等脆性材料設(shè)計出合理可靠的夾具，以及夾具與桿件以螺紋的形式連接對應(yīng)力波傳播的影響問題都值得進(jìn)行深入研究。