LY12鋁合金在沖擊絕熱壓縮下的泊松比

譚華,俞宇穎,戴誠達,譚葉

(中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室,四川綿陽 621900)

LY12鋁合金在沖擊絕熱壓縮下的泊松比

譚華,俞宇穎,戴誠達,譚葉

(中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室,四川綿陽 621900)

對泊松比的意義進行討論,給出了單軸應變加載下線彈性區的泊松比與縱向應力和橫向應力的關系,以及準彈性區的泊松比與聲速的關系。通過實驗測量LY12鋁合金在沖擊壓縮下的縱波和體波聲速,計算在20～131 GPa沖擊壓力范圍內的泊松比。實驗結果表明:LY12鋁合金在發生沖擊熔化前泊松比基本保持為常數,約為0.32左右;發生沖擊熔化后混合相區的泊松比快速增加;當接近完全熔化時泊松比趨于理論極限值0.5.泊松比測量為確定金屬的沖擊熔化壓力區間提供了一種有用的方法。

固體力學;泊松比;聲速;單軸應變;準彈性;沖擊熔化

0 引言

在單軸應變加載下,現有實驗技術能夠對縱向應力和應變進行精確測量,但是依然難以對橫向應力進行直接測量。雖然在低應力加載下可以用錳銅計測量橫向應力,但是橫向計的精確標定相當困難。由于應力計襯墊材料的影響和實驗安裝問題,特別是安裝橫向錳銅計的縫隙方向與沖擊波的傳播方向平行,縫隙中的填充材料在縱向沖擊應力作用下產生的應力擾動會對橫向錳銅計產生附加的干擾[1],影響橫向應力測量,導致橫向應力實測數據發生異常并常常受到各種質疑。本文推導了在單軸應變加載下線彈性區材料橫向應力與泊松比的關系。因此,只要能夠設法測量單軸應變加載下線彈性區的泊松比,就可以利用泊松比計算線彈性變形下的橫向應力,對不同填充材料對橫向錳銅計測量及標定的影響開展研究。

另一方面,基于縱向模量、體積模量和剪切模量與泊松比之間的一般性關系,以及單軸應變加載下的模量與縱波、體波或橫波聲速的關系,可以將泊松比與聲速直接聯系起來。經過近些年的發展,我們能夠比較精確地測量沿著Hugoniot線的聲速,以及從Hugoniot狀態卸載時沿著準彈性卸載路徑的聲速[1-2],使實驗測量和研究沿著沖擊絕熱線的泊松比成為可能。

1 基本模型

1.1 線彈性變形下的泊松比與縱向應力和橫向應力的關系

在力學中,泊松比是用單軸應力彈性變形條件下的橫向應變εy與縱向應變εx之比來定義的[3],

取壓縮應變為正,拉伸應變為負,以保證泊松比具有正值。在泊松比的直觀定義中,桿的邊側(橫向)處于自由應力狀態(如圖1所示),縱向應變和橫向應變均是由縱向應力σx引起的。利用自由邊側條件下縱向彈性壓縮模量即楊氏模量E的定義,可將材料的橫向應變εy表示為

圖1 泊松比的定義Fig.1 A sketch of the definition of Poisson's ratio

在單軸應變實驗中,材料處于三軸應力作用下且橫向應變等于0.顯然不能直接根據(1)式計算單軸應變下的泊松比。在單軸應變線彈性變形條件下,假設施加于橫向(y方向)的應力為σy,則它在x方向產生附加的縱向彈性應變ε′x,

因此三軸應力加載線彈性變形條件下x軸方向的總彈性應變等于上述應變的代數疊加,即

(8)式給出了單軸應變線彈性變形下的縱向應力與橫向應力及泊松比的關系。這樣,只要能夠同時測量泊松比和縱向應力σx,就可以知道橫向應力σy,這為研究單軸應變條件下線彈性材料的橫向應力提供了一種方法。

1.2 單軸應變加載下的泊松比與聲速

眾所周知,平面單軸應變加載下的縱向應力σx與平均應力及縱向偏應力S可表示為

對于各向同性材料,平均應力為(σx+2σy)/3=p,偏應力為2(σx-σy)/3.按照聲速的定義,即小擾動應力波的傳播速度[4]。這種小擾動應力波既可以是等熵的,也可以是非等熵的,即準等熵的。在小擾動應力作用下,最大剪應變γ與縱向應變εx相等[3]:γ=εx.根據(9)式容易得到一維應變加載下的縱波聲速cl與體波聲速cb和橫波聲速ct的普遍性關系,

式中:最大分解剪應力 τ=τmax=(σx-σy)/2=3S/4.根據體積模量、剪切模量與楊氏模量和泊松比的一般性關系[3]:

因此,只要能夠測量精確聲速,就能確定泊松比。顯然,當固體材料發生完全熔化處于液相狀態時,材料失去抗剪能力:τ=0,因此G=0或ct=0.或者處于彈-塑性屈服轉變狀態時,剪應力達到極值:τ=τmax,因此也有G=0或ct=0成立。在這兩種情況下均有ν=1/2,這就是泊松比極限值。

在強沖擊壓縮下,只要材料不處于完全沖擊熔化狀態,從沖擊壓縮狀態卸載時或受到再加載時金屬材料將表現出彈-塑性轉變,這種彈性稱為準彈性。Yu等[2]的實驗結果表明,沿著準彈性卸載路徑,剪切模量和體積模量并不保持為常數,線彈性變形假設在準彈性區不再成立。因此準彈性下的應力-應變關系及變形特性不能用線彈性關系來描述。也就是說(8)式在準彈性區不成立。但(17)式及(18)式是從一般的彈性力學關系推導出來的,因此依然能夠通過聲速測量計算準彈性狀態下的剪切模量和泊松比。

(10)式中的聲速是熱力學聲速或歐拉聲速,從一維應變實驗測量到的聲速是拉氏聲速。拉氏聲速a與歐拉聲速c的關系[4-5]為

式中:ρ0表示初始態密度。(17)式和(18)式用拉氏聲速可表示為

總之,聲速測量為我們提供了一種確定準彈性區的泊松比的重要方法。由于Hugoniot狀態是卸載路徑的起始點,因此(17)式～(21)式也可用來計算沿著沖擊絕熱線的泊松比。

1.3 Hugoniot彈性極限下的泊松比及屈服強度

Hugoniot彈性極限被定義為材料在沖擊加載下發生彈-塑性屈服的起始應力。在理想彈性-塑性模型假設下,并忽略彈性變形到塑性變形的轉變過程,若達到Hugoniot彈性極限σHEL時材料的應變為εHEL,則對應的流體靜水壓p=KεHEL.假定上、下屈服面以流體靜水壓線對稱,則屈服強度Y與σHEL及流體靜水壓強p之間的關系為

在沖擊加載實驗中,通過多臺階樣品實驗測量獲得粒子速度剖面并判定Hugoniot彈性極限對應的粒子速度uHEL,結合拉氏聲速aHEL計算Hugoniot彈性極限為

根據Hugoniot彈性的定義,此時材料的應力狀態恰好位于上屈服面,由屈服強度的定義,結合(8)式得到Hugoniot彈性極限時的屈服強度為

2 LY12鋁合金沿著Huogoniot的泊松比實驗測量

在高達數十至數百吉帕沖擊加載下,處于固相區的金屬材料依然會表現出某種彈性特征,這種彈性稱為“準彈性”。它表明Hugoniot狀態不一定在屈服面上。然而金屬材料在高壓沖擊加載下發生了彈-塑性屈服也是不爭的事實,二者之間似乎存在“矛盾”。

強沖擊加載應力波的前沿可達亞納秒甚至皮秒量級,應變率高達109s-1甚至更高,這種極端應力-應變率過程發生在沖擊波陣面上;但沖擊波加載的終態即Hugoniot狀態是熱力學平衡態,其應變率接近于0.從沖擊波陣面到波后Hugoniot狀態,應變率發生了急劇改變。應變率的急速改變導致剪應力的松弛[6],這可能是導致Hugoniot狀態偏離屈服面的主要原因。由此看來,用橫向應力計測量的Hugoniot狀態下的橫向應力并不代表在沖擊壓縮屈服面上的橫向應力,或者說根據錳銅計測得的橫向應力計算的縱向應力與橫向應力之差,并不能代表沖擊加載下材料的屈服強度。

圖2給出了Yu等[2]測量的國產商用LY12鋁在20～100 GPa壓力范圍內6個Hugoniot壓力點(LY1～LY6)下的拉氏縱波聲速與粒子速度的關系,以及從沖擊壓縮狀態卸載時聲速隨粒子速度的變化。實驗中利用中國工程物理研究院流體物理研究所的二級輕氣炮驅動鋁飛片到高速撞擊鋁樣品產生強沖擊壓縮;用VISAR激光速度干涉儀或DISAR位移干涉儀測量鋁樣品與LiF窗口的界面粒子速度剖面。根據粒子速度剖面計算鋁樣品在沖擊加載下和沿著卸載路徑的聲速[2,4]。這 6個數據點的Hugoniot狀態列于表1中。圖3中對高于100 GPa的沖擊壓力點,用鉭飛片進行非對稱碰撞實驗,聲速的計算方法見文獻[2,4],最高壓力點為131 GPa,已經進入鋁的沖擊熔化固液混合相區。

圖2 實驗測量的LY12鋁沿著Hugoniot的拉氏聲速以及從沖擊壓縮狀態卸載時沿著準彈性-塑形卸載路徑的拉氏聲速隨粒子速度的變化Fig.2 The measured Lagrange sound velocities of LY12 Al alloy materials along the Hugoniot and the change of Lagrange sound velocity with particle velocity along the quasi elastic-plastic release path

由于缺乏LY12鋁的零壓聲速數據,表中列出了2024鋁合金和6061鋁合金的零壓聲速數據供參考。根據聲速數據利用(21)式計算了LY12鋁合金在沖擊壓縮下的泊松比νH,列于表1的最后一列。從20～100 GPa范圍內6發對稱碰撞實驗的結果來看,這一壓力范圍內LY12鋁在沖擊壓縮下的泊松比基本保持常數,其平均值約等于0.32,如圖3的虛線所示。

表1 LY12鋁合金沿著沖擊絕熱線的聲速測量結果與泊松比Tab.1 Measured sound velocities along the Hugoniot and the corresponding Poisson's ratio

圖3 鋁的泊松比隨加載壓力的變化Fig.3 Poisson's ratios of LY12 alloy materials vs.shock pressures

3 討論

Reinhart等測量了6061-T6鋁的聲速并報道了根據聲速計算了泊松比[8]。他們的結果表明,6061-T6鋁合金在43～115 GPa沖擊壓力范圍內的泊松比的變化很小,約為0.33,與本文的結果非常接近。Brown等用光分析法測量了鉭在約150～400 GPa壓力區的聲速[10]。在鉭發生沖擊熔化(約300 GPa)前的固相區,他們假定們泊松比與沖擊壓力的關系計算得到的縱波聲速隨壓力的變化與他們的實驗測量結果符合。考慮到光分析法本身的不確定度,實際也可以認為固相區的泊松比近似為常數。

沿著沖擊絕熱線的體波聲速易于根據實測沖擊絕熱線數據和Gruneisen物態方程計算得到[2,4],而橫波聲速測量較困難。沖擊壓縮下金屬材料在固相區的泊松比基本保持為常數這一事實,為利用泊松比和Hugoniot數據估算沖擊高壓下的縱波和橫波聲速提供了極大方便。根據實測的和泊松比和Hugoniot彈性極限可以直接計算Hugoniot彈性極限時的屈服強度YHEL.文獻[9]給出了2024鋁的實測值σHEL=0.6 GPa及YHEL=0.29 GPa,按照ν=0.33用(28)式計算得到2024鋁在Hugoniot彈性極限時的屈服強度為0.30 GPa,與實驗測量結果十分接近。

鋁在125 GPa沖擊加載下發生沖擊熔化。表1中列出了McQueen利用光分析法測量的125 GPa壓力下的縱波聲速[9],以及本文根據鋁的沖擊絕熱線和Gruneisen物態方程計算的在125 GPa壓力點的體波聲速,據此計算的泊松比也列于表1和圖3.根據Yu等[2]測量的LY12鋁在131 GPa沖擊加載下的粒子速度波剖面的結構,雖然已經進入鋁的沖擊熔化固液混合相區,但從波剖面依然可以觀察到顯著的準彈性-塑性卸載特征,并非像通常認為的那樣一旦發生沖擊熔化材料就進入流體狀態不存在彈-塑性卸載特征。這表明即使發生了沖擊熔化,只要熔化沒有完成,材料依然具有一定的強度。根據實測波剖面計算了LY12鋁的縱波和體波聲速,得到131 GPa下的泊松比約為0.4,說明部分沖擊熔化使LY12鋁的泊松比迅速增加。泊松比的增大意味著橫波聲速下降,導致材料的剪切模量減小,意味著材料抵抗剪切加載能力的下降。泊松比的這一特性為判定金屬的沖擊熔化提供了一種方法。有趣的是,Reinhart等在研究6061-T6鋁合金的強度特性時也發現了類似的結果[8]。他們根據6061-T6鋁在145 GPa的縱波和體波聲速計算得到泊松比為0.47,而6061鋁發生完全沖擊熔化時壓力大約在160 GPa左右,此時根據聲速計算的泊松比達到了0.5.因此,聲速測量和泊松比為我們提供了一種確定從起始沖擊熔化到完全沖擊熔化的壓力區間的有用方法。

4 小結

本文從泊松比的定義出發,推導了在單軸應變加載條件下的泊松比與縱波聲速、橫波聲速及體波聲速之間的關系;將泊松比的應用范圍從經典的單軸應力彈性加載狀態拓展到了單軸應變彈性和和準彈性加載狀態。

根據實驗測量在 20～131 GPa沖擊壓縮下LY12鋁合金的縱波和體波聲速,計算了LY12鋁合金的泊松比,并與國外類似的研究結果進行了比較。發現LY12鋁合金的泊松比在固相區基本保持為常數0.32,與國外同類鋁合金的結果非常接近;在發生沖擊熔化后,在固-液混合相區LY12鋁合金的泊松比隨沖擊壓力迅速增加;發生完全熔化后泊松比趨于理論極限值0.5.

泊松比研究為根據物態方程給出的體波聲速計算沿著沖擊絕熱線的縱波聲速和橫波聲速提供了方法,也提供了一種確定從起始沖擊熔化到完全沖擊熔化的壓力區間有用方法。

References)

[1] Kanel G I,Razorenov S V,Fortov V E.Shock-wave phenomena and the properties of condensed matter[M].NY:Springer-Verlag New York Inc,2004.

[2] Yu Y Y,Tan H,Hu J,et al.Determination of effective shear modulus of shock-compressed LY12 Al from particle velocity profile measurements[J].Journal of Applied Physics,2008,103: 103529.

[3] 楊桂通.彈塑性力學引論[M].北京:清華大學出版社,2000.

YANG Gui-tong.Introduction to elastic-plastic mechanics[M]. Beijing:Tsinghua University Press,2000.(in Chinese)

[4] 譚華.實驗沖擊波物理導引[M].北京:國防工業出版社, 2007.

TAN Hua.Introduction to experimental shock-eave physics[M]. Beijing:National Defense Industry Press,2007.(in Chinese)

[5] 譚華.拉氏坐標與守恒方程[J].高能量密度物理,2012(1): 5-10.

TAN Hua.Lagrangian coordinate and the conservation laws[J]. High Energy Density Phyics,2012(1):5-10.(in Chinese)

[6] Bringa E M,Caro A,Wang Y M,et al.Ultrahigh strength in nanocrystalline materials under shock loading[J].Science.2005, 309(5742):1838-1841.

[7] Marsh S P.Shock L.Hugoniot data[M].Berkley,CA:University of California Press,1980.

[8] Reinhart W D,Asay J R,Chhabildas C,et al.Investigation of 6061-T6 aluminum strength properties to 160 GPa[J].AIP Conference Proceedings,2009,1195:977-980.

[9] McQueen R G,Hopson J W,Fritz J N.Optical technique for determining rarefaction wave velocities in very high pressure[J].Review of Scientific Instruments,1982,53(2):245.

[10] Brown J M,Shaner J W.Rarefaction velocities in shocked Tantalum and the high pressure melting point[M]∥Asay J R,Graham R A,Straub G K.Shock Waves in Condensed Matter-1983.Amsterdam,Holland:North Holland Physics Publishing,1984: 91-94.

Poisson's Ratio of LY12 Aluminum Alloy under Shock Compression

TAN Hua,YU Yu-ying,DAI Cheng-da,TAN Ye
(National Key Laboratory of Shock Wave and Detonation Physics,Institute of Fluid Physics, China Academy of Engineering Physics,Mianyang 621900,Sichuan,China)

The physical meaning of Poisson's ratio is discussed.The relationships among the Poisson's ratio and the longitudinal and transverse stresses under uniaxial strain shock loadings,and the relationship between the Poisson's ratio and the sound velocity in quasi-elastic deformation region are presented. Poisson's ratios of LY12 aluminum alloy materials are calculated based on the sound velocities measured at shock pressures from 20 GPa to 131 GPa.It is found that,while the Poisson's ratios keep almost constant,i.e.,0.32,before shock-melting happens,it increases rapidly in the shock-melting solid-liquid mixed phase region and approachs the theoretical limit of Poisson's ratio,i.e.,0.5,as shock-melting completes,which may provide us a useful way for identifying the shock induced melting zone.

solid mechanics;Poisson's ratio;sound velocity;uniaxial strain;quasi-elasticity;shockinduced melting

O347.5

:A

:1000-1093(2014)08-1218-05

10.3969/j.issn.1000-1093.2014.08.013

2013-08-20

國家自然科學基金項目(10972206)

譚華(1942—),男,研究員,博士生導師。E-mail:zhn_tan@263.net