趙文杰 楊紹普 任學紅 文桂林

摘? ?要：采用RPL100電子式蠕變疲勞試驗機，對高鐵齒輪鋼18CrNiMo7-6開展了不同應變幅值下的對稱應變循環(huán)實驗，研究了齒輪鋼在對稱應變循環(huán)下的循環(huán)軟硬化行為;對齒輪鋼18CrNiMo7-6在不同應力工況下開展了非對稱應力循環(huán)下單軸棘輪行為實驗，研究了齒輪鋼棘輪變形的演化規(guī)律. 結果顯示，在實驗研究循環(huán)周次內齒輪鋼18CrNiMo7-6在對稱應變控制循環(huán)下表現(xiàn)出循環(huán)軟化的特征，在非對稱應力控制循環(huán)下呈現(xiàn)出衰減的棘輪應變率和常棘輪應變率兩階段. 通過在Ohno-Abdel-Karim非線性隨動硬化律中引入與累積塑性應變相關的棘輪參數(shù)的演化方程，建立了修正的循環(huán)彈塑性本構模型. 模擬結果顯示，該模型能夠合理預測齒輪鋼18CrNiMo7-6的循環(huán)軟化特征和棘輪行為的演化規(guī)律.

關鍵詞：高鐵齒輪鋼;循環(huán)軟化;棘輪行為;本構模型

Abstract：To study the cyclic softening/hardening behavior of high-speed railway gear steel 18CrNiMo7-6， RPL 100 electronic creep fatigue tester was used to carry out the uniaxial symmetrical strain cycling experiments on the gear steel under different strain amplitudes. Then， the ratchetting behaviour of gear steel 18CrNiMo7-6 under different stress conditions was tested to study the ratchetting evolution of gear steel under asymmetric stress cycling. The results show that the gear steel（18CrNiMo7-6） exhibits the cyclic softening characteristics under strain cycling and exhibits two stages of decaying ratcheting strain rate and constant ratcheting strain rate under asymmetric stress cycling. Then，a modified elasto-plastic cyclic constitutive model was established by introducing the evolution equation of ratchetting parameters related to the cumulative plastic strain into the Ohno-Abdel-Karim nonlinear kinematic hardening rule. The simulated results show that the model can reasonably predict the cyclic softening characteristics and the ratchetting evolution? of gear steel（18CrNiMo7-6）.

Key words：high-speed railway gear steel;cyclic softening;ratchetting behaviour;constitutive models

齒輪傳動系統(tǒng)是高速列車走行部關鍵部件，在服役過程中遭受輪軌隨機激勵、電機諧波振動及齒輪嚙合交變循環(huán)載荷的作用，在極端工況下齒面或齒根局部會發(fā)生塑性變形的累積，進而發(fā)生破壞. 因此對齒輪材料開展循環(huán)變形行為研究是很有必要的.

近年來，國內外很多學者開展了對典型金屬材料的循環(huán)變形實驗及其本構模型的理論研究[1-4]，并且現(xiàn)有的實驗研究主要集中在不銹鋼、鋁合金和軌道鋼等典型工程材料的單軸和多軸棘輪行為研究以及工程材料的循環(huán)軟硬化特性. 但是對于高鐵齒輪鋼 18CrNiMo7-6循環(huán)變形行為的實驗研究，目前為止相關的研究報道還未見. 在本構模型的理論研究方面，基于已有的實驗結果，眾多學者基于Armstrong-Frederick非線性隨動硬化模型[5]，發(fā)展和提出了更加合理的能夠描述材料循環(huán)軟硬化行為和棘輪行為的隨動硬化模型，以此來提高模型的預測能力，比較具有代表性且應用比較廣泛的模型有Chaboche模型[6]、Ohno-Wang模型[7]、Jiang-Sehitoglu模型[8]、Ohno-Abdel-Karim模型[9]、Ohno-Kang模型[10]等. 但是各本構模型都有自己的優(yōu)缺點，目前為止沒有一種模型可以對所有材料的各種實驗條件和加載工況進行準確的模擬，而且關于高鐵齒輪鋼18CrNiMo7-6這種軟化材料的棘輪行為的預測，相關的研究工作還比較少.

因此，在本文中首先對齒輪鋼18CrNiMo7-6展開單拉實驗、對稱應變循環(huán)實驗和非對稱應力循環(huán)實驗研究，然后在循環(huán)彈塑性模型理論框架下，在Ohno-Abdel-Karim非線性隨動硬化律中引入與累積塑性應變相關的棘輪參數(shù)的演化方程，建立了修正的Ohno-Abdel-Karim隨動硬化模型，最后應用該模型模擬了齒輪鋼18CrNiMo7-6的循環(huán)變形行為，驗證了改進的模型的合理性.

1? ?實? ?驗

本實驗所采用的齒輪鋼18CrNiMo7-6是進口材料，從德國西門子齒輪輪齒根部原位取樣切割加工得到，試樣取樣位置及機加工流程如圖1所示. 材料為18CrNiMo7-6，其主要的化學成分（質量分數(shù)，%）為C：0.15～0.21;Si：≤0.40;Mn：0.50～0.90;P：≤0.035;S≤0.035;Cr：1.50～1.80;Ni：1.40～1.70;Mo：0.25～0.35. 將試樣加工成標距段長度為30 mm，直徑為8 mm的標準螺紋試樣. 試樣加工尺寸如圖2所示.

本文的實驗機采用RPL100電子式蠕變疲勞實驗機，如圖3所示. 具體的實驗內容包括單拉實驗、對稱應變循環(huán)實驗和非對稱應力循環(huán)實驗.

對齒輪鋼18CrNiMo7-6進行了單拉實驗，加載的應變率分別為0.002% s-1和0.02% s-1，得到的應力應變關系曲線如圖4所示，可以看出齒輪鋼18CrNiMo7-6在不同加載率下的單拉實驗曲線基本上是重合的，可以認為齒輪鋼18CrNiMo7-6黏性很小，基本上可以忽略. 同時可由實驗計算出齒輪鋼的彈性模量E=210 GPa，屈服強度為610 MPa，強度極限為1 224 MPa. 還可以看出齒輪鋼在達到強度極限之前會發(fā)生明顯的應變強化，之后會進入頸縮階段.

對齒輪鋼18CrNiMo7-6在應變速率0.2% s-1下開展不同應變幅值下的循環(huán)軟硬化行為實驗，應變幅值分別設置為±0.5%和±0.6%，得到的循環(huán)應力應變關系曲線以及響應的應力幅值與循環(huán)周次的關系如圖5所示. 由圖5可知，齒輪鋼18CrNiMo7-6在應變控制循環(huán)實驗中表現(xiàn)出循環(huán)軟化特性，即響應的應力幅值隨循環(huán)周次的增加逐漸減小. 同時外加應變幅值越大，齒輪鋼的循環(huán)軟化速率越快. 此外該材料的循環(huán)軟化速率會隨著循環(huán)周次增加而逐漸降低，但在實驗研究的循環(huán)周次內并未出現(xiàn)軟化行為飽和現(xiàn)象.

對齒輪鋼18CrNiMo7-6在加載應力速率為100 MPa/s下開展了非對稱應力循環(huán)下棘輪行為實驗，主要研究了不同應力工況對齒輪鋼的棘輪變形的影響. 圖6（a）給出了在應力工況 （200±800）MPa下的應力應變滯回環(huán)的演化曲線，可以看出滯回環(huán)寬度較窄，在實驗周次內應力應變循環(huán)滯回環(huán)的寬度基本保持不變，會隨著循環(huán)周次向前演化，但滯回環(huán)向前演化的速率比較緩慢.

圖6（b）給出了齒輪鋼18CrNiMo7-6在平均應力為200 MPa，不同應變幅值下的棘輪應變演化曲線;圖6（c）給出了齒輪鋼18CrNiMo7-6在應力幅值為750 MPa，不同平均應力下的棘輪應變演化曲線. 從圖中可以看出，18CrNiMo7-6在實驗研究周次內，棘輪應變演化曲線可以分為兩個階段，即棘輪應變率的衰減階段和常棘輪應變率階段，但是在實驗周次內未出現(xiàn)加速棘輪應變率階段. 還可以看出齒輪鋼18CrNiMo7-6的棘輪應變會受到應力加載水平的影響，即平均應力和應力幅值越大，棘輪應變越大.

2? ?循環(huán)彈塑性本構模型

2.1? ?主控方程

由于齒輪鋼18CrNiMo7-6的變形行為與加載速率無關，可認為是一種無黏性特性的材料，因此在本文中采用循環(huán)彈塑性本構模型對齒輪鋼的循環(huán)變形行為進行模擬. 在初始各向同性假設和小變形框架下，循環(huán)彈塑性本構模型的主控方程為[9]：

2.2? ?隨動硬化律和各向同性軟化律

本文基于Ohno-Abdel-Karim模型[9]提出了一種修正的非線性隨動硬化模型來描述齒輪鋼18CrNiMo7-6的棘輪行為，在該模型中總背應力被劃分為M個分量，即

2.3? ?材料參數(shù)的確定

背應力分量的個數(shù)M取值越大，理論上模擬的結果會越精確，但是會增大計算量，在本文中取M = 8. ri和ξi的值可以由單拉塑性應變和應力曲線的實驗結果確定，因為齒輪鋼18CrNiMo7-6為循環(huán)軟化材料，在計算時有必要除去循環(huán)軟化對材料單拉曲線的影響，同時在曲線變化比較大的過渡段需要多選一些點，然后由下式計算得出.

與棘輪參數(shù)有關的參數(shù)μ0、e1和e2只會影響應力循環(huán)下棘輪應變的演化規(guī)律，基本不會影響應變循環(huán)實驗曲線的模擬，可由任一種實驗工況下的棘輪應變曲線確定. 最后由實驗確定的齒輪鋼18CrNiMo7-6的材料參數(shù)如表1所示.

3? ?模擬結果

采用彈塑性循環(huán)本構模型和表1中所示參數(shù)，首先對齒輪鋼18CrNiMo7-6在應變率為0.002% s-1下的單拉曲線進行了模擬，結果如圖7所示. 由圖7可以看出，經(jīng)過修正的模型可以很好地模擬材料的單拉曲線. 圖8（a）給出修正的模型在應變幅值為±0.6%時循環(huán)應力應變曲線的模擬結果，圖8（b）給出了應變循環(huán)幅值分別為0.5%和0.6%時齒輪鋼18CrNiMo7-6響應的應力幅值的模擬結果. 由圖8可以看出，修正的模型不僅可以很好地模擬齒輪鋼的應變控制循環(huán)下的應力應變曲線，還可以對齒輪鋼的循環(huán)軟化特征進行很好的預測，即能對齒輪鋼的循環(huán)軟化速率隨著循環(huán)周次增加而逐漸下降的現(xiàn)象進行合理模擬. 同時該模型也能對應變幅值越大，齒輪鋼的循環(huán)軟化越大的現(xiàn)象進行合理的模擬.

由于棘輪參數(shù)μ的演化對單拉曲線和應變控制的循環(huán)實驗的預測基本沒有影響，只會對棘輪應變的演化規(guī)律有影響. 因此，圖9給出了改進的模型與Ohno-Abdel-Karim模型預測的齒輪鋼18CrNiMo7-6棘輪應變演化曲線的對比結果. 其中，圖9（a）給出了平均應力200 MPa，應力幅值為800 MPa時應力應變滯回環(huán)的模擬曲線，可以看出滯回環(huán)曲線隨著循環(huán)周次向前演化，但是模擬的滯回環(huán)曲線的寬度相比實驗略顯寬. 圖9（b）和9（c）給出了改進的模型對齒輪鋼在不同工況下棘輪行為的預測結果. 圖9（d）和9（e）給出了Ohno-Abdel-Karim模型預測的棘輪行為的演化曲線，通過對比可以發(fā)現(xiàn)，原Ohno-Abdel-Karim模型預測的棘輪應變?yōu)槌＜啈兟市袨?，而改進的模型不僅可以對棘輪應變逐漸增大的現(xiàn)象進行合理的模擬，還可以對棘輪應變率由大逐漸變小，最后保持為一個接近于零的常數(shù)的演化規(guī)律進行合理的預測.

4? ?結? ?論

高鐵齒輪鋼18CrNiMo7-6變形行為與加載速率無關，因此是一種無黏性特性的材料，該材料在對稱應變循環(huán)實驗中表現(xiàn)出循環(huán)軟化的特性，并且加載的應變幅值越大，材料的軟化效應越明顯. 在非對稱應力循環(huán)實驗中會出現(xiàn)棘輪行為，在實驗研究周次內表現(xiàn)出棘輪應變率衰減段和常棘輪應變率段兩階段，但在實驗周次內未出現(xiàn)加速棘輪應變率段. 通過在Ohno-Abdel-Karim隨動硬化模型中引入棘輪參數(shù)演化方程，不僅能夠對高鐵齒輪鋼18CrNiMo7-6單拉曲線和應變循環(huán)下材料的軟化特性進行合理的預測，而且能夠有效地預測棘輪應變和棘輪應變率的演化規(guī)律.

參考文獻

[1]? ? ZHU Y L，KANG G Z，YU C. A finite cyclic elasto-plastic constitutive model to improve the description of cyclic stress-strain hysteresis loops[J]. International Journal of Plasticity，2017，95：191—215.

[2]? CHABOCHE J L，NOUAILHAS D. Constitutive modeling of ratchetting effects，part I：experimental facts and properties of the classical models[J]. Journal of Engineering Materials & Technology，1989，111（4）：384—392.

[3]? ? 邱博，闞前華，劉宇杰，等. 高強軌道鋼雙軸壓-扭棘輪行為的實驗和模擬[J]. 工程力學，2015，32（7）：229—235.

QIU B，KAN Q H，LIU Y J，et al. Experiment and simulation on the biaxial compression-torsion ratchetting behavior of high strength rail steel[J]. Engineering Mechanics，2015，32（7）：229—235. （In Chinese）

[4]? ? ZHAO W J，YANG S P，WEN G L，et al. Fractional-order visco-plastic constitutive model for uniaxial ratcheting behaviors[J]. Applied Mathematics and Mechanics，2019，40（1）：49—62.

[5]? ? FREDERICK C O ，ARMSTRONG P J . A mathematical representation of the multiaxial Bauschinger effect[J]. Materials at High Temperatures，2007，24（1）：1—26.

[6]? ? CHABOCHE J L. On some modifications of kinematic hardening to improve the description of ratchetting effects[J]. International Journal of Plasticity，1991，7（7）：661—678.

[7]? ? OHNO N，WANG J D. Kinematic hardening rules with critical state of dynamic recovery，part Ⅱ：application to experiments of ratchetting behavior[J]. International Journal of Plasticity，1993，9（3）：391—403.

[8]? ? JIANG Y Y，SEHITOGLU H. Cyclic ratchetting of 1070 steel under multiaxial stress states[J]. International Journal of Plasticity，1994，10（5）：579—608.

[9]? ?ABDEL-KARIM M ，OHNO N . Kinematic hardening model suitable for ratchetting with steady-state[J]. International Journal of Plasticity，2000，16（3/4）：225—240.

[10]? KANG G，OHNO N，NEBU A. Constitutive modeling of strain range dependent cyclic hardening[J]. International Journal of Plasticity，2003，19（10）：1801—1819.

[11]? WU D L，XUAN F Z，GUO S J，et al. Uniaxial mean stress relaxation of 9%～12% Cr steel at high temperature：experiments and viscoplastic constitutive modeling[J]. International Journal of Plasticity，2016，77：156—173.

[12]? CHABOCHE J L，NOUAILHAS D. Constitutive modeling of ratchetting effects，part Ⅱ：possibilities of some additional kinematic rules[J]. Journal of Engineering Materials & Technology，1989，111（4）：384—416.

[13]? CHABOCHE J L. A review of some plasticity and viscoplasticity constitutive theories[J]. International Journal of Plasticity，2008，24（10）：1642—1693.