摘要： 形狀記憶合金在受到強(qiáng)沖擊荷載作用時會發(fā)生相變，而相變對其結(jié)構(gòu)件的動態(tài)力學(xué)響應(yīng)影響顯著。本文采用同時考慮靜水壓力和偏應(yīng)力影響的相變臨界準(zhǔn)則，推導(dǎo)了增量型的相變本構(gòu)模型?；趶V義特征理論，得到了復(fù)雜應(yīng)力狀態(tài)下特征波速的解析表達(dá)式。特征波速不僅和材料本身的力學(xué)參數(shù)（如拉壓不對稱性和混合相的模量）有關(guān)，還和材料所處的應(yīng)力狀態(tài)有關(guān)。對因相變導(dǎo)致體積膨脹的TiNi 合金而言，拉壓不對稱性系數(shù)的增大會提高慢波的波速，而對快波幾乎沒有影響。在相變橢圓的短軸處（α = 90°），慢波的波速最低，并隨混合相無量綱模量的增大而顯著減低，混合相無量綱模量由2 增加至5 時，波速降低幅度為36.2%，而快波的波速達(dá)到最大值c0，和混合相的模量無關(guān)；在相變橢圓的長軸處（α = 180°），慢波的速度達(dá)到最大值，而快波的波速達(dá)到最小值c2。

關(guān)鍵詞： 相變；復(fù)合應(yīng)力波；薄壁管；TiNi 合金；偽彈性

中圖分類號： O347.4 國標(biāo)學(xué)科代碼： 13015 文獻(xiàn)標(biāo)志碼： A

固體材料在受到高速撞擊、爆炸等強(qiáng)沖擊荷載作用時，可能會發(fā)生屈服甚至相變。1956 年，Bancroft 等[1] 采用炸藥爆炸加載的方法發(fā)現(xiàn)了鐵的α 相至ε 相的轉(zhuǎn)變，從而迎來了沖擊相變的研究熱潮。沖擊加載下加載點(diǎn)的信息是通過應(yīng)力波的形式向外傳播的，相變信息也將由應(yīng)力波攜帶并傳播。形狀記憶合金（shape memory alloy， SMA）是一種典型的相變材料，在一定條件下處于偽彈性狀態(tài)或形狀記憶狀態(tài)。由于SMA 兼具功能性和結(jié)構(gòu)強(qiáng)度，在抗沖擊方面具有獨(dú)特的優(yōu)越性，在軍事及工程領(lǐng)域得到越來越多的關(guān)注[2-3]。在強(qiáng)沖擊載荷作用下，相變后材料的力學(xué)性能會發(fā)生根本性變化，相變會強(qiáng)烈改變介質(zhì)中沖擊波的波形，在材料的內(nèi)部會形成三波結(jié)構(gòu)波陣面以及卸載沖擊波，對材料和結(jié)構(gòu)的動態(tài)響應(yīng)和破壞特性產(chǎn)生重大影響。因此，研究相變材料內(nèi)應(yīng)力波的傳播規(guī)律，對探索新的現(xiàn)象和機(jī)制、設(shè)計更加有效的抗沖擊結(jié)構(gòu)具有重要意義[4]。

目前，關(guān)于沖擊相變對應(yīng)力波傳播規(guī)律的影響方面的研究，仍以一維沖擊加載為主。Chen 等[5] 和Bekker 等[6] 針對沖擊載荷作用下半無限長桿中相變對應(yīng)力波的傳播影響開展了研究，但未涉及相關(guān)的卸載問題。王文強(qiáng)等[7] 利用CdS 型相變本構(gòu)模型[8] 研究了加、卸載條件下相邊界的演化規(guī)律。Tang 等[9]提出了一種通過沖擊加載制備梯度材料的方法，適當(dāng)調(diào)整卸載的應(yīng)力邊界條件，可以使處于混合相的材料變?yōu)樽兠芏忍荻炔牧?，這種材料的材料組分連續(xù)變化，在材料內(nèi)部并沒有宏觀上的邊界。Dai 等[10] 通過分析和計算的方法，研究了相變應(yīng)力波的傳播演化規(guī)律，發(fā)現(xiàn)卸載相邊界存在分岔、加速等奇異現(xiàn)象，研究結(jié)果表明這些現(xiàn)象和一階相變的強(qiáng)非線性有關(guān)。此外，Niemczura 等[11] 通過實(shí)驗(yàn)方法研究了沖擊拉伸條件下多晶NiTi 合金相邊界的傳播規(guī)律?；趦?nèi)變量相變本構(gòu)模型，Zhu 等[12] 對矩形載荷作用下桿件內(nèi)相變對應(yīng)力波的傳播進(jìn)行了系統(tǒng)研究。和常規(guī)彈塑性材料不同的是，處于偽彈性狀態(tài)的相變材料在強(qiáng)間斷卸載后會產(chǎn)生卸載沖擊波，經(jīng)過多次相互作用后會在沖擊端產(chǎn)生拉應(yīng)力區(qū)。Liu 等[13] 借助于碲鎘汞紅外探測裝置研究了溫度對應(yīng)力波的影響，結(jié)果表明溫度和相變相互耦合，會強(qiáng)烈影響應(yīng)力波的傳播。

大量的研究表明，不同加載條件下的相變現(xiàn)象會有所不同。Plietsch 等[14] 通過一維拉壓實(shí)驗(yàn)，發(fā)現(xiàn)TiNi 合金在拉伸和壓縮條件下存在強(qiáng)度差異效應(yīng)，即拉-壓不對稱性。Orgéas 等[15] 對各向同性的多晶TiNi 合金開展了一系列的實(shí)驗(yàn)，包括準(zhǔn)靜態(tài)拉伸、壓縮和剪切，指出拉-壓不對稱性是TiNi 合金固有的性質(zhì)。在復(fù)雜應(yīng)力狀態(tài)下，Lexcellent 等[16] 針對薄壁圓管進(jìn)行了拉-扭組合加載實(shí)驗(yàn)，從晶體學(xué)角度出發(fā)給出了應(yīng)力空間中的相變臨界面預(yù)測模型，理論模型和實(shí)驗(yàn)結(jié)果基本吻合。此外，還有學(xué)者對TiNi 薄壁圓管開展了雙軸比例加載及非比例加載下的力學(xué)性能實(shí)驗(yàn)[17-22]，這些實(shí)驗(yàn)多局限于準(zhǔn)靜態(tài)條件下，目的是尋求可以精確反映相變行為的本構(gòu)模型，而在復(fù)合應(yīng)力下相變材料內(nèi)應(yīng)力波傳播方面的研究并不多。最近，Song 等[23-24] 建立了薄壁圓管內(nèi)相變復(fù)合應(yīng)力波傳播的控制方程，發(fā)現(xiàn)在圓管的端部突加純扭載荷至材料處于混合相狀態(tài)，其內(nèi)會先出現(xiàn)彈性壓縮波，展現(xiàn)出相變影響下截然不同的應(yīng)力波圖譜，但材料完全進(jìn)入馬氏體相后應(yīng)力波的傳播規(guī)律，卻沒有涉及。考慮到材料進(jìn)入二相后的硬化效應(yīng)，Wang 等[25-26]采用數(shù)值方法重新研究了這一問題，發(fā)現(xiàn)了波系結(jié)構(gòu)的變化現(xiàn)象和波前狀態(tài)的自調(diào)整過程。馬氏體相區(qū)的彈性波追趕混合相區(qū)的復(fù)合應(yīng)力波，從而在管內(nèi)形成沖擊波，Cui 等[27] 把混合相區(qū)和馬氏體相區(qū)統(tǒng)一起來進(jìn)行考慮，對TiNi 合金管內(nèi)復(fù)合應(yīng)力波的傳播現(xiàn)象進(jìn)行了數(shù)值研究。

管作為一種典型的構(gòu)件，在抗沖擊、吸能領(lǐng)域有著廣泛的應(yīng)用。TiNi 材料在相變過程中不僅形狀會發(fā)生改變，而且體積也會發(fā)生變化[4， 28]。材料進(jìn)入混合相后復(fù)合應(yīng)力波的波速并不是恒定的，而和材料所處的應(yīng)力狀態(tài)有關(guān)，這在以往的文獻(xiàn)中并沒有過多的探討。了解管內(nèi)應(yīng)力波的傳播速度，可以用于判斷結(jié)構(gòu)缺陷的位置、評估結(jié)構(gòu)損傷程度，為結(jié)構(gòu)安全評估和維護(hù)提供有力支撐。本文以處于偽彈性狀態(tài)的TiNi 合金薄壁圓管為研究對象，基于廣義特征理論構(gòu)建相變應(yīng)力波波速的解析表達(dá)式，并分析影響應(yīng)力波波速的因素。

1 控制方程

1.1 基本方程

2.1 拉壓不對稱性對波速的影響

對于TiNi 合金而言，相變過程中由于伴隨著體積發(fā)生變化，在宏觀力學(xué)性能上呈現(xiàn)出拉-壓不對稱性，從而使得相變橢圓在σ?τ 平面上沿著σ 軸移動。相變橢圓關(guān)于σ 軸對稱，橢圓的上半部分如圖3 所示，其中，τM 為應(yīng)力橢圓上的最大剪應(yīng)力；k1 為相變起始時的相變功相關(guān)參數(shù)。

圖4 顯示了拉壓不對稱性對應(yīng)力波波速的影響，圖中：實(shí)線和虛線分別代表無量綱快波和慢波；黑色、藍(lán)色和紫色線分別代表α=0， 0.1和0.2，其中泊松比ν=0.3，Em/E=0.1。從圖中可以看出，α的變化對快波的影響并不大；而慢波的波速卻隨著α的增大而增大。同時，快波的波速是先增大后減小，在切應(yīng)力最大處快波的波速最大，其值為c0；而慢波的波速隨著參數(shù)α的增大先減小后增大，應(yīng)力狀態(tài)在σ軸（α=0o或180°）上慢波速度達(dá)到最大為c2。此外，對于不考慮拉壓不對稱的材料（α=0），關(guān)于橢圓短軸（τ 軸）對稱的應(yīng)力狀態(tài)，即應(yīng)力狀態(tài)（σ，τ）和（?σ，τ），快波和慢速的波速并不改變，如常規(guī)的彈塑性材料；而對于拉壓不對稱材料（α≠0，根據(jù)式（29） 和圖4，即使應(yīng)力狀態(tài)關(guān)于短軸對稱，快波和慢波的波速卻并不相同：以拉壓不對稱系數(shù)α=0.2 的材料為例，當(dāng)α=30°時，無量綱慢波和快波的速度分別為0.285 和0.755，而α=150°無量綱慢波和快波的速度分別為0.322 和0.774。

2.2 混合相模量對波速的影響

圖5 給出了混合相模量對相變應(yīng)力波波速的影響，圖中虛線表示無量綱慢波波速，實(shí)線表示無量綱快波波速。從圖5 上可以看出，混合相的模量對波速的影響程度和應(yīng)力狀態(tài)高度相關(guān)。對于快波，當(dāng)α = 90°時，快波的速度為c0，和混合相的模量無關(guān)；α 增大或者減小時，快波波速顯著降低—當(dāng)α 增大至150°時，混合相無量綱模量由2 增加至5 時，無量綱快波波速由0.851降至0.783，降低幅度為8.0%。對于慢波，隨著混合相模量的增大，慢波的波速顯著減低—當(dāng)α = 90°時，混合相無量綱模量由2 增加至5 時，無量綱慢波波速由0.445 降低至0.284，降低幅度為36.2%。

3 結(jié) 論

針對TiNi 合金相變后存在剪切變形和體積膨脹的特性，利用綜合考慮靜水壓力和偏應(yīng)力的相變臨界準(zhǔn)則，得到了增量型的相變本構(gòu)模型及處于偽彈性狀態(tài)下TiNi 合金薄壁管內(nèi)應(yīng)力波傳播的控制方程。借助于廣義特征理論，得到了應(yīng)力波傳播速度的函數(shù)表達(dá)式。研究結(jié)果表明，相變區(qū)內(nèi)應(yīng)力波波速不僅與材料本身的力學(xué)參數(shù)有關(guān)，還與材料所處的應(yīng)力狀態(tài)有關(guān)。對于相變后體積膨脹的TiNi合金而言，由于材料拉壓不對稱性的存在，應(yīng)力狀態(tài)關(guān)于相變橢圓短軸對稱時，快波和慢波的波速并不具有對稱性。在相變橢圓的短軸處（α = 90°），慢波的波速最低，并隨混合相無量綱模量的增大而顯著減低，混合相無量綱模量由2 增加至5 時，波速降低幅度為36.2%，而快波的波速達(dá)到最大值c0，和混合相的模量無關(guān)；在相變橢圓的長軸處（α = 180°），慢波的速度達(dá)到最大值，而快波的波速達(dá)到最小值c2。

參考文獻(xiàn)：

[1]BANCROFT D， PETERSON E L， MINSHALL S. Polymorphism of iron at high pressure [J]. Journal of Applied Physics，1956， 27（3）： 291–298. DOI： 10.1063/1.1722359.

[2]GUIDA M， SELLITTO A， MARULO F， et al. Analysis of the impact dynamics of shape memory alloy hybrid composites for"advanced applications [J]. Materials， 2019， 12（1）： 153. DOI： 10.3390/ma12010153.

[3]GUPTA A K， VELMURUGAN R， JOSHI M， et al. Studies on shape memory alloy-embedded GFRP composites for improved"post-impact damage strength [J]. International Journal of Crashworthiness， 2019， 24（4）： 363–379. DOI： 10.1080/13588265.2018.1452549.

[4]唐志平. 相變應(yīng)力波 [M]. 北京： 科學(xué)出版社， 2022： 194–247.

TANG Z P. Stress waves with phase transition [M]. Beijing： Science Press， 2022： 194–247.

[5]CHEN Y C， LAGOUDAS D C. Impact induced phase transformation in shape memory alloys [J]. Journal of the Mechanics"and Physics of Solids， 2000， 48（2）： 275–300. DOI： 10.1016/S0022-5096（99）00044-7.

[6]BEKKER A， JIMENEZ-VICTORY J C， POPOV P， et al. Impact induced propagation of phase transformation in a shape"memory alloy rod [J]. International Journal of Plasticity， 2002， 18（11）： 1447–1479. DOI： 10.1016/S0749-6419（02）00025-6.

[7]王文強(qiáng)， 唐志平. 沖擊下宏觀相邊界的傳播 [J]. 爆炸與沖擊， 2000， 20（1）： 25–31. DOI： 10.3321/j.issn：1001-1455.2000.01.005.

WANG W Q， TANG Z P. Propagation of macroscopic phase boundary under shock loading [J]. Explosion and Shock Waves，2000， 20（1）： 25–31. DOI： 10.3321/j.issn：1001-1455.2000.01.005.

[8]TANG Z P， GUPTA Y M. Shock-induced phase transformation in cadmium sulfide dispersed in an elastomer [J]. Journal of"Applied Physics， 1988， 64（4）： 1827–1837. DOI： 10.1063/1.341782.

[9]TANG Z P， DAI X. A preparation method of functionally graded materials with phase transition under shock loading [J].Shock Waves， 2006， 15（6）： 447–452. DOI： 10.1007/s00193-006-0048-8.

[10]DAI X Y， TANG Z P， XU S L， et al. Propagation of macroscopic phase boundaries under impact loading [J]. International"Journal of Impact Engineering， 2004， 30（4）： 385–401. DOI： 10.1016/s0734-743x（03）00090-3.

[11]NIEMCZURA J， RAVI-CHANDAR K. Dynamics of propagating phase boundaries in NiTi [J]. Journal of the Mechanics and"Physics of Solids， 2006， 54（10）： 2136–2161. DOI： 10.1016/j.jmps.2006.04.003.

[12]ZHU P P， DAI H H. Wave propagation in a shape memory alloy bar under an impulsive loading [J]. Journal of Applied"Mechanics， 2016， 83（10）： 104502. DOI： 10.1115/1.4034115.

[13]LIU Y G， SHEN L Y， CHEN Y J， et al. Thermomechanical coupling effect on the phase transition wave propagation in an"SMA TiNi bar subjected to shock loading [J]. International Journal of Mechanical Sciences， 2022， 235： 107710. DOI： 10.1016/j.ijmecsci.2022.107710.

[14]PLIETSCH R， EHRLICH K. Strength differential effect in pseudoelastic NiTi shape memory alloys [J]. Acta Materialia， 1997，45（6）： 2417–2424. DOI： 10.1016/S1359-6454（96）00354-0.

[15]ORGéAS L， FAVIER D. Stress-induced martensitic transformation of a NiTi alloy in isothermal shear， tension and"compression [J]. Acta Materialia， 1998， 46（15）： 5579–5591. DOI： 10.1016/S1359-6454（98）00167-0.

[16]LEXCELLENT C， BLANC P. Phase transformation yield surface determination for some shape memory alloys [J]. Acta"Materialia， 2004， 52（8）： 2317–2324. DOI： 10.1016/j.actamat.2004.01.022.

[17]GRABE C， BRUHNS O T. Path dependence and multiaxial behavior of a polycrystalline NiTi alloy within the pseudoelastic"and pseudoplastic temperature regimes [J]. International Journal of Plasticity， 2009， 25（3）： 513–545. DOI： 10.1016/j.ijplas.2008.03.002.

[18]MEHRABI R， ANDANI M T， KADKHODAEI M， et al. Experimental study of NiTi thin-walled tubes under uniaxial tension，torsion， proportional and non-proportional loadings [J]. Experimental Mechanics， 2015， 55（6）： 1151–1164. DOI：10.1007/s11340-015-0016-2.

[19]WANG X M， ZHOU Q T， LIU H， et al. Experimental study of the biaxial cyclic behavior of thin-wall Tubes of NiTi shape"memory alloys [J]. Metallurgical and Materials Transactions A， 2012， 43（11）： 4123–4128. DOI： 10.1007/s11661-012-1225-2.

[20]FARAJPOUR M R， SHAHIDI A R， FARAJPOUR A. A nonlocal continuum model for the biaxial buckling analysis of"composite nanoplates with shape memory alloy nanowires [J]. Materials Research Express， 2018， 5（3）： 035026. DOI： 10.1088/2053-1591/aab3a9.

[21]SITTNER P， HARA Y， TOKUDA M. Experimental study on the thermoelastic martensitic transformation in shape memory"alloy polycrystal induced by combined external forces [J]. Metallurgical and Materials Transactions A， 1995， 26（11）：2923–2935. DOI： 10.1007/bf02669649.

[22]SUN Q P， LI Z Q. Phase transformation in superelastic NiTi polycrystalline micro-tubes under tension and torsion–from"localization to homogeneous deformation [J]. International Journal of Solids and Structures， 2002， 39（13/14）： 3797–3809.DOI： 10.1016/S0020-7683（02）00182-8.

[23]SONG Q Z， TANG Z P. Combined stress waves with phase transition in thin-walled tubes [J]. Applied Mathematics and"Mechanics， 2014， 35（3）： 285–296. DOI： 10.1007/s10483-014-1791-7.

[24]宋卿爭. 復(fù)合加載下NiTi 合金力學(xué)特性和相變波的研究 [D]. 合肥： 中國科學(xué)技術(shù)大學(xué)， 2014： 95–120.

SONG Q Z. Mechanical properties and phase transition waves of NiTi alloy under combined stresses [D]. Hefei： University of"Science and Technology of China， 2014： 95–120.

[25]WANG B，" TANG Z P. Study on the propagation of coupling shock waves with phase transition under combined tension-torsion impact loading [J]. Science China Physics， Mechanics amp; Astronomy， 2014， 57（10）： 1977–1986. DOI： 10.1007/s11433-014-5468-3.

[26]王波. 相變材料及聚合物中的復(fù)合應(yīng)力波研究 [D]. 合肥： 中國科學(xué)技術(shù)大學(xué)， 2017： 23–41.

WANG B. Research on stress waves of phase transition material and polymer under combined stress [D]. Hefei： University of"Science and Technology of China， 2017： 23–41.

[27]CUI S T， LIANG L L. Influence of phase transformation on stress wave propagation in thin-walled tubes [J]. Waves in"Random and Complex Media. DOI： 10.1080/17455030.2022.2164631.

[28] LAGOUDAS D C. Shape memory alloys： modeling and engineering applications [M]. New York： Springer， 2008.

[29]李永池. 波動力學(xué) [M]. 合肥： 中國科學(xué)技術(shù)大學(xué)出版社， 2015： 210–223.

LI Y C. Wave mechanics [M]. Hefei： China University of Science and Technology Press， 2015： 210–223.

[30]AURICCHIO F， PETRINI L. A three-dimensional model describing stress-temperature induced solid phase transformations：thermomechanical coupling and hybrid composite applications [J]. International Journal for Numerical Methods in"Engineering， 2004， 61（5）： 716–737. DOI： 10.1002/nme.1087.

[31]QIDWAI M A， LAGOUDAS D C. On thermomechanics and transformation surfaces of polycrystalline NiTi shape memory"alloy material [J]. International Journal of Plasticity， 2000， 16（10/11）： 1309–1343. DOI： 10.1016/S0749-6419（00）00012-7.

[32]郭揚(yáng)波， 唐志平， 徐松林. 一種考慮靜水壓力和偏應(yīng)力共同作用的相變臨界準(zhǔn)則 [J]. 固體力學(xué)學(xué)報， 2004， 25（4）： 417–422.DOI： 10.3969/j.issn.0254-7805.2004.04.009.

GUO Y B， TANG Z P， XU S L. A critical criterion for phase transformation considering both hydrostatic pressure and"deviatoric stress effects [J]. Acta Mechanica Solida Sinica， 2004， 25（4）： 417–422. DOI： 10.3969/j.issn.0254-7805.2004.04.009.

（責(zé)任編輯 王小飛）