王江波， 李汶峰， 任文科， 王曉東， 李繼亮， 高光發(fā)

(南京理工大學(xué) 機(jī)械工程學(xué)院，南京 210094)

分離式霍普金森壓桿(split Hopkinson pressure bar，SHPB)[1]自從建立以來(lái)，被廣泛應(yīng)用于混凝土[2]、陶瓷[3]、金屬[4]、巖石[5]、復(fù)合材料[6]、聚合物[7-8]和泡沫材料[9]等的動(dòng)態(tài)力學(xué)響應(yīng)測(cè)試，可獲得材料在102～104s-1應(yīng)變率范圍內(nèi)的應(yīng)力應(yīng)變曲線。為了確保試驗(yàn)結(jié)果的準(zhǔn)確性，試驗(yàn)過(guò)程中需要滿足恒應(yīng)變率加載及試件內(nèi)部應(yīng)力均勻等基本條件。

在大多數(shù)情況下，常規(guī)的梯形入射脈沖無(wú)法滿足恒應(yīng)變率變形的要求，同時(shí)還會(huì)產(chǎn)生一系列高頻分量的應(yīng)力波，每一個(gè)頻率分量的應(yīng)力波在傳播時(shí)由于速度不同，會(huì)導(dǎo)致入射波形發(fā)生彌散效應(yīng)[10]，這種現(xiàn)象不僅會(huì)影響試驗(yàn)時(shí)的動(dòng)態(tài)應(yīng)力平衡，而且有可能會(huì)掩蓋材料本身的力學(xué)特性。因此，為了能得到更準(zhǔn)確的試驗(yàn)數(shù)據(jù)，需要對(duì)入射波的波形進(jìn)行修正，研究表明主要有改變撞擊桿的形狀[11]、使用預(yù)加載桿技術(shù)和脈沖整形器技術(shù)[12]。由于前兩種方法有很多弊端如撞擊桿形狀加工困難、很難實(shí)現(xiàn)恒應(yīng)變率加載等，所以波形整形器就成為應(yīng)用最廣泛的脈沖整形技術(shù)，它是通過(guò)對(duì)入射波形的調(diào)整來(lái)改變反射波形，其原理是在撞擊桿和入射桿端面之間放置一個(gè)片狀材料，只要選取適當(dāng)?shù)牟牧虾统叽?，就可以最大程度減小SHPB試驗(yàn)中波形的彌散效應(yīng)，同時(shí)可以實(shí)現(xiàn)動(dòng)態(tài)應(yīng)力平衡和恒應(yīng)變率加載等要求。

由于波形整形器使用簡(jiǎn)單，能對(duì)入射波形進(jìn)行有效的整形，故很多學(xué)者都對(duì)此進(jìn)行了研究。Frew等[13-14]基于材料體積不可壓縮和均勻變形理論，得到了整形器的分析模型，并利用整形器技術(shù)進(jìn)行了巖石、玻璃陶瓷和鋼的SHPB試驗(yàn)，結(jié)果表明，使用整形器后，試件處于動(dòng)態(tài)應(yīng)力平衡狀態(tài)，試驗(yàn)的大多數(shù)時(shí)間具有恒定的應(yīng)變率，還可以通過(guò)調(diào)整整形器尺寸控制入射波的上升持續(xù)時(shí)間，分析模型與試驗(yàn)結(jié)果一致，同時(shí)還可以預(yù)測(cè)試件應(yīng)力平衡和恒應(yīng)變率持續(xù)時(shí)間。Naghdabadi等[15]進(jìn)行了整形器對(duì)入射波影響規(guī)律的試驗(yàn)研究和數(shù)值仿真，給出了不同材料SHPB試驗(yàn)時(shí)整形器尺寸的選擇準(zhǔn)則。Bagher Shemirani 等[16]使用紫銅整形器進(jìn)行了混凝土材料的SHPB桿試驗(yàn)，并對(duì)試驗(yàn)過(guò)程中動(dòng)態(tài)應(yīng)力平衡、恒應(yīng)變率狀態(tài)進(jìn)行了分析，并且給出撞擊桿長(zhǎng)度和速度對(duì)入射波形的影響規(guī)律。

從以上分析可知，學(xué)者們只是針對(duì)某一口徑SHPB裝置或特定材料進(jìn)行了整形器對(duì)入射波構(gòu)形的影響規(guī)律研究，而關(guān)于使用整形器的入射波構(gòu)形是否滿足相似律的研究未曾出現(xiàn)。研究此問(wèn)題可將不同桿徑SHPB裝置的入射波整形技術(shù)及理論進(jìn)行統(tǒng)一，并且對(duì)新型SHPB裝置的設(shè)計(jì)具有非常重要的科學(xué)意義?；诖耍疚睦昧烤V分析理論研究了使用整形器的SHPB入射波形構(gòu)形是否滿足嚴(yán)格的幾何相似律，并根據(jù)試驗(yàn)和數(shù)值模擬結(jié)果對(duì)此進(jìn)行了驗(yàn)證，然后基于直徑為14.5 mm的SHPB裝置研究了整形器直徑、厚度以及撞擊桿長(zhǎng)度和速度對(duì)入射波形的影響規(guī)律，最后，通過(guò)對(duì)各影響因素的定量分析，得到了對(duì)入射波構(gòu)形的兩個(gè)重要拐點(diǎn)T1和T2的無(wú)量綱表達(dá)式。

1 量綱分析

整形器為紫銅材料，由于在試驗(yàn)過(guò)程中整形器的塑性應(yīng)變遠(yuǎn)大于彈性應(yīng)變，加載過(guò)程中整形器的彈性加載時(shí)間遠(yuǎn)小于塑性加載時(shí)間，因此在量綱分析中不考慮整形器的彈性階段，即可以忽略整形器的楊氏模量和泊松比，則其應(yīng)力應(yīng)變曲線可以近似等效為線性硬化模型。設(shè)整形器的塑性硬化模量為E0，屈服強(qiáng)度為σs，密度為ρs，厚度為h，直徑為d。假設(shè)SHPB裝置桿直徑為Db，撞擊桿長(zhǎng)度為L(zhǎng)b，撞擊速度為v，密度為ρb，楊氏模量為Eb。當(dāng)SHPB桿裝置的長(zhǎng)度遠(yuǎn)大于直徑時(shí)，我們可認(rèn)為其滿足一維應(yīng)力假設(shè)，故可以不考慮桿材料的泊松比。因此，我們可以給出入射桿中測(cè)量點(diǎn)截面上的應(yīng)力函數(shù)表達(dá)式

σ=f(Eb,ρb,Lb,Db;E0,σs,ρs,h,d;v,t)

(1)

基于量綱分析理論，我們選擇質(zhì)量(M)，長(zhǎng)度(L)和時(shí)間(T)作為基本量綱，式(1)中所有參數(shù)的量綱列于表1中。

表1 桿和整形器參數(shù)的量綱

(2)

如果沒有使用整形器時(shí)，撞擊桿的直徑和撞擊速度對(duì)無(wú)量綱應(yīng)力并無(wú)影響，可以忽略式(2)中與整形器參數(shù)無(wú)關(guān)的前兩項(xiàng)，同時(shí)由于整形器的厚度很小，不考慮應(yīng)力波在整形器中的傳播時(shí)間，即忽略整形器的慣性效應(yīng)，則可以不考慮密度的影響，則可以得到

(3)

式中：σ*=ρbCbv/2；t*=Lb/Cb?？s比模型中桿材料和整形器與原模型相同時(shí)，即可得到

(4)

式中：下標(biāo)m為原模型；p為縮比模型。

因此，對(duì)于撞擊桿速度v相同時(shí)，要使使用整形器的SHPB入射波構(gòu)形滿足幾何相似律，則必須有

(5)

2 試驗(yàn)研究和數(shù)值模擬

2.1 試驗(yàn)方法

圖1 SHPB裝置系統(tǒng)示意圖

(6)

式中：Ab為桿的橫截面積；As和Ls為試件的橫截面積和長(zhǎng)度

為了驗(yàn)證SHPB裝置的相似律問(wèn)題，以及后續(xù)研究整形器對(duì)入射波形的影響規(guī)律，試驗(yàn)選用常用直徑為14.5 mm的SHPB裝置，該裝置撞擊桿、入射桿和透射桿的規(guī)格尺寸，如表2所示。

表2 SHPB裝置中桿的尺寸和材料

整形器材料的選取與待測(cè)試件材料的性能相關(guān)，一般選用紫銅、黃銅等材料，在加載過(guò)程中其應(yīng)力應(yīng)變關(guān)系可以近似為線性硬化模型。本文主要針對(duì)金屬、高聚物等一些延性材料SHPB試驗(yàn)時(shí)整形器尺寸的選取進(jìn)行研究。采用合適材料和尺寸的整形器不僅可以減小試驗(yàn)過(guò)程中波形彌散對(duì)試驗(yàn)結(jié)果的影響，而且可以實(shí)現(xiàn)試件內(nèi)部的動(dòng)態(tài)應(yīng)力平衡和試驗(yàn)過(guò)程中的恒應(yīng)變率加載。

2.2 數(shù)值模擬模型

2.2.1 整形器材料力學(xué)參數(shù)

選取T2紫銅作為整形器材料，為了整形器仿真材料模型的準(zhǔn)確性，對(duì)紫銅材料進(jìn)行了準(zhǔn)靜態(tài)壓縮試驗(yàn)。紫銅試件準(zhǔn)靜態(tài)壓縮的應(yīng)力應(yīng)變曲線，如圖2所示。根據(jù)相關(guān)標(biāo)準(zhǔn)可以測(cè)得其屈服強(qiáng)度、楊氏模量、以及塑性切線模量，同時(shí)利用密度計(jì)測(cè)量了其密度，如表3所示。

圖2 紫銅材料應(yīng)力應(yīng)變曲線及各參數(shù)的計(jì)算

表3 紫銅整形器的主要力學(xué)參數(shù)

2.2.2 數(shù)值仿真幾何模型

利用有限元軟件ABAQUS對(duì)使用整形器的SHPB試驗(yàn)進(jìn)行了數(shù)值仿真，模型包括撞擊桿、整形器和入射桿。本文主要研究整形器對(duì)入射波的影響規(guī)律，故在建模過(guò)程中忽略了透射桿。由于整形器尺寸較小，為了保證桿與整形器接觸端網(wǎng)格尺寸相匹配，桿的有限元模型采用過(guò)渡網(wǎng)格劃分方法。因此桿的仿真模型包括細(xì)化網(wǎng)格區(qū)(Ⅰ)、過(guò)渡網(wǎng)格區(qū)(Ⅱ)和粗化網(wǎng)格區(qū)(Ⅲ)。

以直徑為14.5 mm的SHPB裝置為例，其有限元模型如圖3所示。模型中細(xì)化網(wǎng)格區(qū)域長(zhǎng)度為30 mm，過(guò)渡網(wǎng)格區(qū)域長(zhǎng)度為20 mm。同時(shí)，為了排除網(wǎng)格尺寸對(duì)數(shù)值模擬結(jié)果的影響，細(xì)化區(qū)的網(wǎng)格為邊長(zhǎng)為0.1 mm的正方形，過(guò)渡區(qū)網(wǎng)格為邊長(zhǎng)為0.725 mm的正三角形，粗化區(qū)網(wǎng)格為0.725 mm×1.450 mm的長(zhǎng)方形。

圖3 ABAQUS軟件中模型的建立及網(wǎng)格的劃分

在整個(gè)試驗(yàn)過(guò)程桿始終處于彈性狀態(tài)，因此數(shù)值仿真模型中入射桿和撞擊桿采用彈性模型，其主要參數(shù)如表4所示(整形器材料采用雙線性模型，其主要參數(shù)已列于表3)。

表4 仿真模型中桿材料的力學(xué)參數(shù)

根據(jù)文獻(xiàn)[17]對(duì)銅試件SHPB研究中界面摩擦因數(shù)的研究并對(duì)此進(jìn)行了微調(diào)，整形片與撞擊桿接觸的一端未添加潤(rùn)滑劑，數(shù)值仿真中設(shè)置摩擦因數(shù)為0.12；而整形片另一端與入射桿通過(guò)潤(rùn)滑劑粘住，摩擦系數(shù)為0.03，同時(shí)對(duì)無(wú)整形器SHPB試驗(yàn)也進(jìn)行了數(shù)值仿真，模型中各部分的接觸定義為面-面接觸。有無(wú)整形片試驗(yàn)與仿真的對(duì)比圖，如圖4所示。從圖4可知，兩者之間有著很好的重復(fù)性，因此可以說(shuō)明仿真模型和材料參數(shù)是準(zhǔn)確可靠的。

圖4 撞擊桿速度為16 m/s時(shí)有無(wú)整形器的入射波形試驗(yàn)和仿真對(duì)比圖

3 SHPB試驗(yàn)入射波相似律的驗(yàn)證

為了驗(yàn)證量綱分析中SHPB試驗(yàn)的幾何相似性，基于直徑為14.5 mm的SHPB裝置，開展了不同縮比尺寸的SHPB試驗(yàn)和數(shù)值仿真研究，具體方案如表5所示。原模型中(幾何縮比為1)撞擊桿長(zhǎng)度為400 mm，整形器的厚度為0.8 mm，直徑為8 mm，縮比模型與原模型材料相同且撞擊桿速度也相同為16 m/s，幾何縮比分別為1.5和0.5。

表5 SHPB裝置的縮比試驗(yàn)與數(shù)值仿真幾何參數(shù)

圖5 SHPB試驗(yàn)中入射波的幾何相似性試驗(yàn)與數(shù)值仿真結(jié)果

根據(jù)SHPB試驗(yàn)相似律量綱分析的結(jié)果和圖5中不同縮比模型的入射波形，當(dāng)原型中桿和整形器材料與縮比模型完全相同，且撞擊速度相同時(shí)，入射波的無(wú)量綱正應(yīng)力與無(wú)量綱時(shí)間之間的關(guān)系圖即無(wú)量綱入射波形基本相同，說(shuō)明考慮整形器的SHPB試驗(yàn)滿足幾何相似律，同時(shí)驗(yàn)證了量綱分析的準(zhǔn)確性。因此，當(dāng)不同直徑的裝置滿足SHPB試驗(yàn)中的兩個(gè)基本假設(shè)后，后續(xù)研究整形器對(duì)入射波的影響規(guī)律可以適用于任何口徑的SHPB裝置。

4 整形器和撞擊桿參數(shù)對(duì)入射波形的影響

對(duì)于入射波的整形問(wèn)題，一般只考慮加載階段，波形的卸載部分對(duì)試驗(yàn)基本沒有影響，因此可以不考慮入射波形卸載部分。SHPB試驗(yàn)中有無(wú)整形器典型的入射波形，如圖6所示。從圖6可知，我們可以發(fā)現(xiàn)，使用整形器的入射波形可以分成兩個(gè)關(guān)鍵節(jié)點(diǎn)和三個(gè)加載區(qū)域，分別為：Ⅰ.彈性加載區(qū)域，即試驗(yàn)時(shí)整形器處于彈性變形的階段；Ⅱ.塑性加載區(qū)域，即試驗(yàn)時(shí)整形器為塑性變形的階段；Ⅲ.剛性加載區(qū)域，即整形器完全塑性變形后直至波形加載結(jié)束。

圖6 入射波形中的關(guān)鍵拐點(diǎn)和加載區(qū)域

由于脈沖整形技術(shù)中影響入射波構(gòu)形的因素較多以及試驗(yàn)時(shí)所處應(yīng)力場(chǎng)的復(fù)雜狀況，所以很難給出整形后入射波形的具體方程式。因此，本文根據(jù)整形器在其工作時(shí)所處的應(yīng)力應(yīng)變狀態(tài)，對(duì)其尺寸對(duì)入射波構(gòu)形兩個(gè)重要拐點(diǎn)T1和T2的影響規(guī)律進(jìn)行了具體分析。事實(shí)上，整形器對(duì)入射波構(gòu)形的整形就是改變T1和T2拐點(diǎn)的位置以及兩拐點(diǎn)之間曲線的形狀(見圖6中Ⅱ區(qū)域)，然后根據(jù)待測(cè)試件的力學(xué)性能調(diào)整入射波的構(gòu)形使試驗(yàn)在恒應(yīng)變率加載下進(jìn)行。然而，圖6中Ⅱ區(qū)域波形很難給出具體的表達(dá)式，故我們對(duì)其起點(diǎn)和終點(diǎn)進(jìn)行了具體的分析，為SHPB試驗(yàn)時(shí)整形器的選取提供一定的參考作用。

當(dāng)考慮某個(gè)影響因素對(duì)入射波某一點(diǎn)無(wú)量綱應(yīng)力的影響時(shí)，可以忽略掉變量時(shí)間的影響，因此，類似于第1章中的量綱分析，我們得到入射波特定點(diǎn)的無(wú)量綱應(yīng)力表達(dá)式為

(7)

后續(xù)研究中針對(duì)式(7)中各因素對(duì)入射波兩個(gè)拐點(diǎn)T1和T2的影響進(jìn)行了具體的分析。

4.1 整形器直徑對(duì)入射波形的影響

基于直徑為14.5 mm的SHPB裝置，進(jìn)行了一系列不同直徑整形器的SHPB試驗(yàn)和數(shù)值仿真。其中，撞擊桿的速度為16 m/s，撞擊桿長(zhǎng)度為400 mm，整形器的厚度為0.8 mm。SHPB試驗(yàn)中使用不同整形器直徑時(shí)的入射波形圖，如圖7所示。

圖7 不同整形器直徑對(duì)應(yīng)的入射波形圖

從圖7可知，當(dāng)整形器厚度和撞擊桿速度不變時(shí)，入射波的T1拐點(diǎn)隨著整形器直徑的增大而增大，這是由于隨著整形器直徑的增大，整形器壓縮過(guò)程中到達(dá)其屈服點(diǎn)時(shí)的面積也隨之增大，而在加載過(guò)程中，整形器和桿端面的軸向力相等，所以導(dǎo)致入射波的T1拐點(diǎn)無(wú)量綱應(yīng)力逐漸增大，同時(shí)T1拐點(diǎn)對(duì)應(yīng)的無(wú)量綱時(shí)間也變大；但整形器直徑對(duì)T2拐點(diǎn)沒有明顯的影響，即整形器直徑主要影響塑性加載區(qū)域Ⅱ的起始點(diǎn)，對(duì)其終點(diǎn)沒有影響。

4.2 整形器厚度對(duì)入射波形的影響

為了分析整形器厚度對(duì)入射波形的影響，開展了不同整形器厚度的SHPB試驗(yàn)和數(shù)值仿真，整形器的直徑為8 mm，撞擊桿長(zhǎng)度和速度同4.1節(jié)，試驗(yàn)和數(shù)值仿真結(jié)果如圖8所示。

圖8 不同整形器厚度對(duì)應(yīng)的入射波形

當(dāng)整形器厚度從0.5 mm增加到1.2 mm時(shí)，入射波形的T1拐點(diǎn)無(wú)量綱應(yīng)力和無(wú)量綱時(shí)間也隨之減小，這是由于壓縮過(guò)程中整形器到達(dá)屈服點(diǎn)(入射波形的T1拐點(diǎn))時(shí)，整形器的厚度越大，其對(duì)應(yīng)的當(dāng)前橫截面積越小，又因?yàn)樵囼?yàn)過(guò)程中整形器和入射桿軸向力相等，故導(dǎo)致T1拐點(diǎn)無(wú)量綱應(yīng)力降低。而T2拐點(diǎn)的無(wú)量綱時(shí)間隨著整形器厚度的增加逐漸增大，波形的塑性加載上升沿逐漸變緩，其原因是整形器在受壓過(guò)程中為均勻變形，隨著厚度的增加，其達(dá)到整形器完全塑性變形時(shí)所需的時(shí)間更長(zhǎng)，故會(huì)導(dǎo)致入射波形的上升沿時(shí)間變長(zhǎng)，即T2拐點(diǎn)右移。

4.3 撞擊桿長(zhǎng)度對(duì)入射波形的影響

在SHPB試驗(yàn)中，入射波的脈沖持續(xù)時(shí)間與撞擊桿的長(zhǎng)度有關(guān)，為研究撞擊長(zhǎng)度對(duì)入射波形的影響，進(jìn)行了不同撞擊桿長(zhǎng)度的SHPB試驗(yàn)與數(shù)值模擬，其中撞擊桿速度為16m/s，整形器直徑和厚度分別為8 mm和0.8 mm，數(shù)值仿真和試驗(yàn)結(jié)果如圖9所示。

圖9 不同撞擊桿長(zhǎng)度的入射波形

從圖9可知，當(dāng)撞擊桿的能量足夠大時(shí)，即入射波形存在剛性加載區(qū)域Ⅲ，整個(gè)上升沿階段不同撞擊桿長(zhǎng)度的入射波形基本一致，不同撞擊桿長(zhǎng)度下的彈性加載區(qū)域Ⅰ和塑性加載區(qū)域Ⅱ重合，故當(dāng)撞擊桿動(dòng)能足夠時(shí)，撞擊桿長(zhǎng)度對(duì)入射波的兩個(gè)關(guān)鍵拐點(diǎn)T1和T2沒有影響，也就是說(shuō)，撞擊桿長(zhǎng)度對(duì)于整形器的整形可以忽略。而當(dāng)撞擊桿動(dòng)能較小時(shí)，入射波應(yīng)力未達(dá)到理論應(yīng)力峰值就會(huì)出現(xiàn)卸載情況，但對(duì)于整形器尺寸和撞擊桿速度相同而撞擊桿長(zhǎng)度不同對(duì)應(yīng)的入射波，其波形與存在剛性加載區(qū)域的入射波形加載區(qū)域Ⅰ和Ⅱ重合，這與Naghdabadi研究中的結(jié)論一致。

4.4 撞擊桿速度對(duì)入射波形的影響

SHPB試驗(yàn)中撞擊桿的速度對(duì)試件的應(yīng)變率有著直接的影響，一般來(lái)說(shuō)，撞擊桿的速度越大，試驗(yàn)時(shí)的應(yīng)變率就越高，故研究撞擊速度對(duì)入射波整形規(guī)律的影響就顯得很重要，撞擊速度分別為8～20 m/s時(shí)對(duì)應(yīng)的入射波形，如圖10所示。其中整形器的直徑和厚度保持不變，分別為8 mm和0.8 mm，撞擊桿的長(zhǎng)度為400 mm。

圖10 不同撞擊桿速度對(duì)應(yīng)的波形圖

根據(jù)波動(dòng)力學(xué)理論，不同撞擊速度的入射波應(yīng)力峰值為σ*=ρbCbv/2，即隨著撞擊速度的增加，入射波的應(yīng)力峰值也逐漸增大。對(duì)入射波形進(jìn)行無(wú)量綱化后，隨著撞擊速度的增大，入射波形的T1拐點(diǎn)的無(wú)量綱應(yīng)力逐漸減小(見圖10)，這是由于在對(duì)應(yīng)力進(jìn)行無(wú)量綱化時(shí)除以理論應(yīng)力峰值導(dǎo)致的，實(shí)際上入射波T1點(diǎn)應(yīng)力值大致相等；入射波形的T2拐點(diǎn)的無(wú)量綱時(shí)間隨著撞擊速度的增大逐漸減小，主要原因是因?yàn)樽矒羲俣雀邥r(shí)導(dǎo)致整形器的變形速率加快，導(dǎo)致其完全塑性變形的時(shí)間減小，故T2拐點(diǎn)的無(wú)量綱時(shí)間逐漸減小。

4.5 T1拐點(diǎn)的定量分析

對(duì)于入射波形T1拐點(diǎn)的取法，具體如下：根據(jù)圖4中試驗(yàn)和仿真的對(duì)比結(jié)果可知，兩個(gè)波形在T1點(diǎn)附近的交叉點(diǎn)約為入射波形第一個(gè)波動(dòng)部分的中點(diǎn)，故為了統(tǒng)一試驗(yàn)和仿真中T1點(diǎn)的取法，統(tǒng)一采用入射波形第一個(gè)震蕩部分的中點(diǎn)為T1點(diǎn)。根據(jù)4.1節(jié)～4.4節(jié)的分析可知，整形器的直徑和厚度以及撞擊桿的速度都會(huì)對(duì)T1拐點(diǎn)有一定的影響，而撞擊桿的長(zhǎng)度對(duì)其沒有影響，同時(shí)，假設(shè)式(6)中的三個(gè)影響因素σs/σ*,d/Db,h/Lb，是相互解耦的，則式(7)可以變?yōu)?/p>

(8)

因此，為了得到T1拐點(diǎn)的無(wú)量綱應(yīng)力表達(dá)式，根據(jù)圖7、圖8和圖10得到了入射波形的T1拐點(diǎn)的無(wú)量綱應(yīng)力與各無(wú)量綱參數(shù)的對(duì)應(yīng)關(guān)系，如圖11所示。從圖11可知，T1拐點(diǎn)的無(wú)量綱應(yīng)力σs/σ*與無(wú)量綱直徑d/Db呈二次函數(shù)關(guān)系，與速度有關(guān)項(xiàng)σs/σ*和整形器無(wú)量綱厚度呈線性關(guān)系，因此我們對(duì)此分別進(jìn)行了多項(xiàng)式和線性擬合，可以得到T1拐點(diǎn)的無(wú)量綱應(yīng)力表達(dá)式為(0<σ/σ*<1)

圖11 T1拐點(diǎn)無(wú)量綱應(yīng)力與各影響因素的關(guān)系

(9)

對(duì)于坐標(biāo)系中的一點(diǎn)，要想確定其準(zhǔn)確位置，需要縱坐標(biāo)和橫坐標(biāo)來(lái)共同確定，故得到入射波T1無(wú)量綱時(shí)間與式(7)右邊的三種無(wú)量綱參數(shù)的對(duì)應(yīng)關(guān)系，如圖12所示。同樣也對(duì)此進(jìn)行了擬合，得到T1拐點(diǎn)的無(wú)量綱時(shí)間表達(dá)式

圖12 T1拐點(diǎn)無(wú)量綱時(shí)間與各影響因素的關(guān)系

(10)

4.6 T2拐點(diǎn)的定量分析

對(duì)于入射波形T2拐點(diǎn)的確定，由于其理論無(wú)量綱應(yīng)力值為1，故選取T2時(shí)即入射波形的震蕩部分中點(diǎn)為1時(shí)的點(diǎn)為T2點(diǎn)。T2拐點(diǎn)的影響因素主要有整形器的厚度和撞擊桿的速度，撞擊桿的長(zhǎng)度和整形器直徑對(duì)其的影響可以忽略不計(jì)，因此式(7)中的K(d/Db)為常數(shù)，而T2拐點(diǎn)的無(wú)量綱應(yīng)力都為1，故只考慮無(wú)量綱時(shí)間t/t*與影響因素h/Lb和σs/σ*的定量關(guān)系。T2拐點(diǎn)無(wú)量綱時(shí)間與各影響因素的對(duì)應(yīng)關(guān)系圖，如圖13所示。發(fā)現(xiàn)無(wú)量綱時(shí)間與兩個(gè)影響因素都呈線性關(guān)系，同樣認(rèn)為無(wú)量綱參數(shù)h/Lb與σs/σ*是相互解耦的，并對(duì)此過(guò)程進(jìn)行了線性擬合，得到T2拐點(diǎn)的無(wú)量綱時(shí)間表達(dá)式(0

圖13 T2拐點(diǎn)無(wú)量綱時(shí)間與各影響因素的關(guān)系

(11)

綜上所述，入射波T1拐點(diǎn)的無(wú)量綱應(yīng)力和時(shí)間可由下面方程組計(jì)算

(12)

入射波T2拐點(diǎn)的無(wú)量綱應(yīng)力和時(shí)間的表達(dá)式如下

(13)

5 結(jié) 論

基于直徑為14.5 mm的SHPB裝置，利用量綱分析理論對(duì)考慮整形器的SHPB試驗(yàn)的相似律問(wèn)題進(jìn)行了研究，在此基礎(chǔ)上，分析了整形器直徑和厚度、撞擊桿長(zhǎng)度和速度對(duì)入射波形的影響規(guī)律，根據(jù)量綱分析的結(jié)果，得到了入射波形兩個(gè)重要拐點(diǎn)T1和T2的無(wú)量綱應(yīng)力和時(shí)間的定量表達(dá)式，主要有以下結(jié)論：

(1)假設(shè)SHPB裝置的原模型和縮比模型中撞擊桿和整形器材料相同且滿足物理相似，當(dāng)撞擊桿的速度保持不變時(shí)，考慮整形器的SHPB裝置滿足嚴(yán)格的幾何相似律。

(2)使用整形器的SHPB試驗(yàn)可以延長(zhǎng)入射波的上升沿時(shí)間，T1拐點(diǎn)的無(wú)量綱應(yīng)力σ/σ*隨著整形器無(wú)量綱直徑d/Db的增大而變大，但隨著撞擊桿速度v和整形器無(wú)量綱厚度h/Lb的增大而逐漸減??；T2拐點(diǎn)的無(wú)量綱時(shí)間與整形器無(wú)量綱厚度h/Lb呈正比，與速度呈反比；撞擊桿的長(zhǎng)度主要影響入射波的持續(xù)時(shí)間，對(duì)入射波的上升沿沒有影響。

(3)SHPB試驗(yàn)脈沖整形技術(shù)是產(chǎn)生一個(gè)與待測(cè)試件應(yīng)力響應(yīng)(透射波形)相似的入射波形，同時(shí)入射波形的T1拐點(diǎn)與反射波平臺(tái)段的起始點(diǎn)呈正相關(guān)，當(dāng)試驗(yàn)應(yīng)變率增大時(shí)，應(yīng)增加整形器的直徑和撞擊桿的速度；T1拐點(diǎn)的無(wú)量綱應(yīng)力σ/σ*和無(wú)量綱時(shí)間t/t*與影響參數(shù)σs/σ*，d/Db和h/Lb分別呈線性、二次函數(shù)和線性關(guān)系；T2拐點(diǎn)的無(wú)量綱應(yīng)力σ/σ*始終為1，無(wú)量綱時(shí)間t/t*與σs/σ*無(wú)量綱參數(shù)h/Lb與σs/σ*都是呈線性關(guān)系的。