顧邦平， 胡 雄， 楊振生， 賴金濤， 潘 龍(.上海海事大學(xué) 物流工程學(xué)院,上海 2006；2.紹興文理學(xué)院 機(jī)械與電氣工程學(xué)院,浙江 紹興 2000；.浙江大學(xué) 浙江省先進(jìn)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室,杭州 0027)

消除構(gòu)件加工制造過程中的殘余應(yīng)力是機(jī)械制造工業(yè)領(lǐng)域中的一項(xiàng)重要研究課題。振動(dòng)時(shí)效技術(shù)因?yàn)榫哂性O(shè)備簡單、處理時(shí)間短、處理效果好等特點(diǎn)而成為備受關(guān)注的殘余應(yīng)力消除技術(shù)。國內(nèi)外專家學(xué)者已經(jīng)對振動(dòng)時(shí)效工藝和理論開展了深入的研究[1-5]，并取得了豐富的理論成果和應(yīng)用實(shí)例。通過國內(nèi)外專家學(xué)者的研究發(fā)現(xiàn)，影響振動(dòng)時(shí)效效果的關(guān)鍵工藝參數(shù)主要包括：激振頻率、激振振幅和激振時(shí)間。目前對于激振頻率的選擇主要采用掃頻法，即對構(gòu)件進(jìn)行掃頻振動(dòng)，尋找構(gòu)件振動(dòng)幅值較大的共振頻率，然后在其共振頻率或者亞共振頻率下進(jìn)行振動(dòng)時(shí)效處理。然而通過掃頻法確定振動(dòng)時(shí)效的激振頻率時(shí)，通常并不考慮材料的殘余應(yīng)力分布狀態(tài)以及振型的分布規(guī)律，極大地限制了振動(dòng)時(shí)效消除和均化殘余應(yīng)力的效果。

針對傳統(tǒng)掃頻法確定振動(dòng)時(shí)效激振頻率存在的不足，沈華龍等[6]在研究振動(dòng)時(shí)效工藝時(shí)，采用有限元軟件對構(gòu)件進(jìn)行數(shù)值模態(tài)分析，獲取構(gòu)件各階位移振型和模態(tài)頻率，然后根據(jù)振型的分布規(guī)律指導(dǎo)振動(dòng)時(shí)效的現(xiàn)場試驗(yàn)，取得了不錯(cuò)的效果。然而構(gòu)件的位移振型相對于應(yīng)變振型對構(gòu)件的局部變化并不敏感，并不能直接去研究構(gòu)件某些關(guān)鍵部位的應(yīng)變分布[7-9]，并且構(gòu)件的高階應(yīng)變振型對構(gòu)件的局部變化相對于低階應(yīng)變振型更加敏感。此外，國內(nèi)外專家學(xué)者也對激振頻率大于1 kHz的高頻振動(dòng)時(shí)效[10-11]和超聲振動(dòng)時(shí)效[12-13]開展了相關(guān)的研究。因此獲取構(gòu)件的應(yīng)變振型能夠更準(zhǔn)確的反映振動(dòng)能量在構(gòu)件表面的分布規(guī)律，這對于指導(dǎo)振動(dòng)時(shí)效的現(xiàn)場試驗(yàn)是有益的。

本文采用ANSYS有限元軟件，對構(gòu)件進(jìn)行數(shù)值模態(tài)分析，通過分析應(yīng)變振型和殘余應(yīng)力的分布規(guī)律，確定了基于應(yīng)變振型的振動(dòng)時(shí)效激振頻率的確定方法和步驟，并通過實(shí)驗(yàn)驗(yàn)證該方法的提出對振動(dòng)時(shí)效效果改善的有效性。

1 數(shù)值模態(tài)分析

1.1 應(yīng)變振型控制方程的推導(dǎo)

文中將以位移振型為基礎(chǔ)，推導(dǎo)應(yīng)變振型與位移振型的關(guān)系，建立應(yīng)變振型的控制方程，為后續(xù)的數(shù)值模態(tài)分析奠定基礎(chǔ)。一個(gè)典型的有限單元是由節(jié)點(diǎn)j，k，m和相應(yīng)的連接直線來定義的，則該單元內(nèi)任意一點(diǎn)的位移ui可以表示為[14]

ui=Niai

(1)

式中：Ni是基于節(jié)點(diǎn)位置事先給定的函數(shù)，稱為形狀函數(shù)矩陣;ai表示第i個(gè)單元所有節(jié)點(diǎn)位移組成的節(jié)點(diǎn)位移向量。

當(dāng)已知單元內(nèi)任意位置的位移，則相應(yīng)的可以確定出單元內(nèi)任意位置的應(yīng)變εi，其表達(dá)式為

εi=Siui

(2)

式中：Si為微分算子矩陣。若系統(tǒng)被劃分為n個(gè)單元，則共有n個(gè)式(2)這樣的方程，可以表示為

(3)

可以簡寫為

ε=Su

(4)

式中：ε為n個(gè)單元的應(yīng)變向量;u為n個(gè)單元的節(jié)點(diǎn)位移向量。式(4)建立的應(yīng)變與位移的關(guān)系是在單元坐標(biāo)系中建立的，考慮到各單元在連接節(jié)點(diǎn)上位移是連續(xù)變化的，則單元節(jié)點(diǎn)位移向量可以通過式(5)的坐標(biāo)變換轉(zhuǎn)換到總體坐標(biāo)系中

u=Tug

(5)

式中：T為坐標(biāo)變換矩陣;ug為在總體坐標(biāo)系中的節(jié)點(diǎn)位移向量。將式(5)代入式(4)中，可以得到總體坐標(biāo)系中的應(yīng)變方程為

(6)

采用有限元模型，在總體坐標(biāo)系中，系統(tǒng)的無阻尼運(yùn)動(dòng)方程為

(7)

(Kg-ω2M)Umax=Fmax

(8)

采用模態(tài)疊加法，節(jié)點(diǎn)位移的響應(yīng)可以表示為

Umax=φYφTFmax

(9)

將式(9)代入到式(6)中，可以得到

(10)

令ψ=STφ，對位移振型進(jìn)行微分運(yùn)算可以得到應(yīng)變振型，而矩陣S為微分算子，因此ψ稱為應(yīng)變模態(tài)振型矩陣，ψr為第r階應(yīng)變模態(tài)振型，與第r階位移模態(tài)振型φr相對應(yīng)。從能量的角度，第r階應(yīng)變模態(tài)振型反應(yīng)的是構(gòu)件第r階固有能量分布狀態(tài)，與固有頻率和位移模態(tài)振型一樣也是反映構(gòu)件固有特性的一個(gè)參數(shù)。

1.2 模型建立及材料設(shè)置

物理模型：數(shù)值模態(tài)分析的試件為后續(xù)實(shí)驗(yàn)中采用的45#鋼板狀構(gòu)件，長230 mm，寬210 mm、厚6 mm，如圖1所示。為了研究構(gòu)件局部結(jié)構(gòu)的變化對其位移振型和應(yīng)變振型的影響，對兩種結(jié)構(gòu)的板狀構(gòu)件進(jìn)行數(shù)值模態(tài)分析，其中板狀構(gòu)件1#的中心無通孔，稱為無孔構(gòu)件，另外一板狀構(gòu)件2#的中心有φ6 mm的通孔，稱為有孔構(gòu)件。由于數(shù)值模態(tài)分析的構(gòu)件三維幾何模型相對簡單，此處選擇直接在ANSYS經(jīng)典GUI界面中建立構(gòu)件的三維幾何模型。

圖1 數(shù)值分析構(gòu)件示意圖(mm)

Fig.1 Schematic of the workpiece used in the numerical analysis(mm)

材料設(shè)置：數(shù)值模態(tài)分析時(shí)需要用到45#鋼的密度ρ、彈性模量E、泊松比ν，具體數(shù)值分別為ρ=7 850 kg/m3、E=200 GPa、ν=0.28。

有限元模型：采用20節(jié)點(diǎn)的SOLID95單元對構(gòu)件劃分網(wǎng)格，建立其三維有限元模型。為了提高求解的精度，在試件的長度、寬度以及厚度三個(gè)方向設(shè)置單元?jiǎng)澐值臄?shù)目分別為55、45以及3，如圖1所示。然后采用掃掠(Volume Sweep)的方式對構(gòu)件劃分網(wǎng)格，可以得到規(guī)則的網(wǎng)格，也可以提高求解的精度。另外，SOLID95單元具有中間節(jié)點(diǎn)，這將使得求解的結(jié)果更為精確。

1.3 求解方法

采用ANSYS模態(tài)分析模塊中的Block Lanczos方法(分塊蘭索斯法)求解構(gòu)件的位移模態(tài)振型和應(yīng)變模態(tài)振型。Block Lanczos法特征值求解器是ANSYS模態(tài)分析模塊中的默認(rèn)求解器，采用的是Lanczos算法。這種算法不僅可以保證求解的精確性，更能提高求解的速度。

1.4 數(shù)值結(jié)果及分析

采用1.3小節(jié)的求解方法求解1.1小節(jié)建立的位移模態(tài)和應(yīng)變模態(tài)的控制方程可以獲得構(gòu)件的位移模態(tài)和應(yīng)變模態(tài)的數(shù)值結(jié)果。無孔構(gòu)件低階模態(tài)振型(模態(tài)頻率為572 Hz)如圖2所示，其中圖2(a)為其位移振型，圖2(b)為其應(yīng)變振型；有孔構(gòu)件低階模態(tài)振型(模態(tài)頻率為571 Hz)如圖3所示，其中圖3(a)為其位移振型，圖3(b)為其應(yīng)變振型。

(a) 位移振型

(b) 應(yīng)變振型圖2 無孔構(gòu)件低階模態(tài)振型(572 Hz)

Fig.2 Low order modal mode of the workpiece without the hole (572 Hz)

無孔構(gòu)件高階模態(tài)振型(模態(tài)頻率為3 231 Hz)如圖4所示，其中圖4(a)為其位移振型，圖4(b)為其應(yīng)變振型；有孔構(gòu)件高階模態(tài)振型(模態(tài)頻率為3 229 Hz)如圖5所示，其中圖5(a)為其位移振型，圖5(b)為其應(yīng)變振型。

(a) 位移振型

(b) 應(yīng)變振型圖3 有孔構(gòu)件低階模態(tài)振型(571 Hz)

Fig.3 Low order modal mode of the workpiece with the hole (571 Hz)

(a) 位移振型

(b) 應(yīng)變振型圖4 無孔構(gòu)件高階模態(tài)振型(3 231 Hz)Fig.4 High order modal mode of the workpiece without the hole(3 231 Hz)

(a) 位移振型

(b) 應(yīng)變振型圖5 有孔構(gòu)件高階模態(tài)振型(3 229 Hz)Fig.5 High order modal mode of the workpiece with the hole(3 229 Hz)

對比無孔構(gòu)件和有孔構(gòu)件的位移振型和應(yīng)變振型，可以發(fā)現(xiàn)構(gòu)件局部結(jié)構(gòu)的變化在其低階和高階應(yīng)變振型上都能夠得以體現(xiàn)，而低階和高階位移振型對于構(gòu)件局部結(jié)構(gòu)的變化都不敏感，同時(shí)從構(gòu)件的應(yīng)變振型上可以直觀地看出構(gòu)件各部分分布的動(dòng)應(yīng)變的大小，有利于直觀地評估構(gòu)件各部分殘余應(yīng)力的釋放效果，這說明了為什么從應(yīng)變振型的角度探討振動(dòng)時(shí)效激振頻率的確定方法。

掃頻法是目前常用的確定振動(dòng)時(shí)效激振頻率的方法，即對構(gòu)件進(jìn)行掃頻振動(dòng)，尋找構(gòu)件的共振頻率，然后在共振頻率下對構(gòu)件進(jìn)行振動(dòng)時(shí)效處理。然而這種方法確定的激振頻率并未考慮到構(gòu)件表面的殘余應(yīng)力分布狀態(tài)和共振頻率所對應(yīng)的應(yīng)變振型，這極大地降低了振動(dòng)時(shí)效的效果。從圖3(b)和圖5(b)的應(yīng)變振型分布規(guī)律可知，高階應(yīng)變振型具有較多的應(yīng)變峰值，更有利于消除殘余應(yīng)力分布狀態(tài)復(fù)雜的構(gòu)件的殘余應(yīng)力。假如構(gòu)件表面的峰值殘余應(yīng)力主要分布在圖3(b)所示的應(yīng)變峰值附近，此時(shí)選取圖3(b)應(yīng)變振型所對應(yīng)的模態(tài)頻率進(jìn)行振動(dòng)時(shí)效處理能夠獲得較為理想的振動(dòng)時(shí)效效果。若構(gòu)件表面的峰值殘余應(yīng)力分布在圖5(b)所示的兩個(gè)應(yīng)變峰值附近，此時(shí)選取圖5(b)應(yīng)變振型所對應(yīng)的模態(tài)頻率進(jìn)行振動(dòng)時(shí)效處理能夠獲得較為理想的振動(dòng)時(shí)效效果。從上面的分析，我們不難發(fā)現(xiàn)對于振動(dòng)時(shí)效的激振頻率的確定，不僅要考慮構(gòu)件表面的殘余應(yīng)力分布狀態(tài)，而且需要考慮構(gòu)件的應(yīng)變振型，兩者缺一不可。因此，為了提高振動(dòng)時(shí)效效果，對于激振頻率f的選擇應(yīng)以構(gòu)件表面的殘余應(yīng)力分布為依據(jù)，共振應(yīng)變函數(shù)ψr為準(zhǔn)則。

在此基礎(chǔ)上本文提出了基于應(yīng)變振型的振動(dòng)時(shí)效激振頻率的確定方法，其步驟主要包括：①獲取構(gòu)件表面的殘余應(yīng)力分布狀態(tài)；②建立構(gòu)件的有限元模型，對構(gòu)件進(jìn)行數(shù)值模態(tài)分析；③獲取構(gòu)件的各階應(yīng)變振型和模態(tài)頻率；④在步驟①和③的基礎(chǔ)上，初步確定振動(dòng)時(shí)效的激振頻率；⑤以數(shù)值模態(tài)分析確定的激振頻率為基礎(chǔ)，在其一定范圍內(nèi)對構(gòu)件進(jìn)行掃頻振動(dòng)，獲取振動(dòng)幅值較大的頻率作為振動(dòng)時(shí)效的激振頻率。通過本文提出的方法確定的激振頻率是以構(gòu)件表面的殘余應(yīng)力分布為依據(jù)，共振應(yīng)變函數(shù)ψr為準(zhǔn)則的，有利于獲得較為理想的振動(dòng)時(shí)效效果。

2 實(shí)驗(yàn)研究

2.1 實(shí)驗(yàn)材料

碳素結(jié)構(gòu)鋼被廣泛應(yīng)用于工程結(jié)構(gòu)中，這種材料在實(shí)際應(yīng)用中由于存在較大的殘余應(yīng)力而經(jīng)常產(chǎn)生損傷破壞，特別是在有些場合不允許采用熱時(shí)效對其進(jìn)行去應(yīng)力處理，因此有必要嘗試采用振動(dòng)時(shí)效對碳素結(jié)構(gòu)鋼的殘余應(yīng)力進(jìn)行消除研究。實(shí)驗(yàn)中選用45#優(yōu)質(zhì)碳素結(jié)構(gòu)鋼來制備實(shí)驗(yàn)試樣。

2.2 試樣制備

實(shí)驗(yàn)試樣(即1.2小節(jié)數(shù)值模態(tài)分析的試件)的尺寸為230 mm×210 mm×6 mm，如圖1所示。為了使得試樣的表面產(chǎn)生較大的拉伸殘余應(yīng)力，實(shí)驗(yàn)中采用激光表面處理工藝來產(chǎn)生殘余應(yīng)力。這種工藝是通過激光與試樣相互作用區(qū)域快速的升溫來改變試樣的微觀組織形態(tài)，試樣的表面產(chǎn)生熱損傷而產(chǎn)生初始拉伸殘余應(yīng)力。激光表面處理工藝在SLC-2030D龍門式CO2多功能激光加工機(jī)上展開，激光器的型號為Rofin SLAB DC030。原理圖如圖6所示。

圖6 激光表面處理工藝的原理圖Fig.6 Schematic of the laser surface processing

2.3 小孔法評估殘余應(yīng)力

試樣殘余應(yīng)力測點(diǎn)以及應(yīng)變花粘貼方向如圖7所示。實(shí)驗(yàn)中對試樣的兩個(gè)區(qū)域進(jìn)行激光表面處理，每個(gè)激光熱損傷區(qū)域的中心距離試樣的邊緣均為55 mm。實(shí)驗(yàn)中選取14個(gè)測點(diǎn)用于評價(jià)振動(dòng)前后試樣的殘余應(yīng)力，其中測點(diǎn)1#～7#被分為一組，用于測試振動(dòng)前試樣的殘余應(yīng)力，測點(diǎn)1′#～7′#被分為另一組，用于測試振動(dòng)后試樣的殘余應(yīng)力。每一組測點(diǎn)對稱的分布在垂直于垂直中心線的直線上，距離垂直中心線的距離分別為0、40 mm、60 mm、80 mm，每組測點(diǎn)距離試樣邊緣的距離均為80 mm。以測點(diǎn)4為例，其應(yīng)變花的粘貼方向如圖7所示，其余測點(diǎn)處應(yīng)變花的粘貼方向與測點(diǎn)4處應(yīng)變花的粘貼方向一致。通過應(yīng)變儀測試試樣鉆孔處的釋放應(yīng)變ε1、ε2、ε3，并將釋放應(yīng)變代入如下的方程就可以求解出試樣鉆孔處x軸向和y軸向的殘余應(yīng)力σx和σy

(11)

圖7 殘余應(yīng)力測點(diǎn)示意圖(mm)Fig.7 Schematic of the measured points for residual stress(mm)

2.4 振動(dòng)時(shí)效過程

振動(dòng)時(shí)效實(shí)驗(yàn)裝置主要包括：電磁式激振器、高頻振動(dòng)能量放大裝置、信號發(fā)生器、功率放大器、壓電式加速度傳感器、電荷放大器、示波器。首先，將試樣通過安裝孔裝夾在高頻振動(dòng)能量放大裝置上；然后，按照本文提出的方法對試樣進(jìn)行掃頻振動(dòng)，發(fā)現(xiàn)試樣在頻率717 Hz處有一個(gè)峰值較大的共振頻率點(diǎn)，以及在頻率3 311 Hz處有一個(gè)峰值較大的共振頻率點(diǎn)；最后，對其中一件試樣(稱為1#試樣)在頻率717 Hz處進(jìn)行振動(dòng)時(shí)效處理，對另一件試樣(稱為2#試樣)在頻率3 311 Hz處進(jìn)行振動(dòng)時(shí)效處理。

2.5 實(shí)驗(yàn)結(jié)果及分析

圖8為振動(dòng)時(shí)效前后1#試樣x軸向和y軸向的殘余應(yīng)力測試結(jié)果。圖9為振動(dòng)時(shí)效前后2#試樣x軸向和y軸向的殘余應(yīng)力測試結(jié)果。

(a) x軸向殘余應(yīng)力

(b) y軸向殘余應(yīng)力圖8 振動(dòng)時(shí)效前后1#試樣殘余應(yīng)力測試結(jié)果

Fig.8 Residual stress measuring results of 1# specimen before and after the vibratory stress relief

1#試樣和2#試樣振動(dòng)前殘余應(yīng)力的測試結(jié)果表明，試樣最大的殘余應(yīng)力出現(xiàn)在激光表面燒傷處理中心區(qū)域，隨著遠(yuǎn)離激光表面燒傷中心區(qū)域距離的增加，殘余應(yīng)力則逐漸降低。1#試樣和2#試樣振動(dòng)時(shí)效處理后殘余應(yīng)力峰值消除率均較大，而對于殘余應(yīng)力較小的地方(遠(yuǎn)離激光表面燒傷中心區(qū)域)，振動(dòng)時(shí)效處理后殘余應(yīng)力的消除率相對較小。由圖8可見，1#試樣經(jīng)過振動(dòng)時(shí)效處理后其x軸向和y軸向峰值殘余應(yīng)力的釋放量分別為106.8 MPa和131.8 MPa；由圖9可見，2#試樣經(jīng)過振動(dòng)時(shí)效處理后其x軸向和y軸向峰值殘余應(yīng)力的釋放量分別為183.6 MPa和226.1 MPa，表明根據(jù)殘余應(yīng)力分布狀態(tài)和應(yīng)變振型分布規(guī)律，合理的選擇振動(dòng)時(shí)效的激振頻率能夠獲得較為理想的振動(dòng)時(shí)效消除殘余應(yīng)力效果。

由圖8和圖9可見，2#試樣具有更好的振后殘余應(yīng)力均化分布效果，表明根據(jù)殘余應(yīng)力分布狀態(tài)和應(yīng)變振型分布規(guī)律，合理的選擇振動(dòng)時(shí)效的激振頻率，可以獲得較為理想的振動(dòng)時(shí)效均化殘余應(yīng)力效果。何聞等學(xué)者在文獻(xiàn)[3]中研究指出對試樣在較高的頻率下激振，試樣上分布較多的應(yīng)變峰值，并且動(dòng)應(yīng)變在試樣上的分布更加均勻，提高了試樣表面的殘余應(yīng)力均化效果；除此以外，通過1.4小節(jié)對試樣應(yīng)變振型的分析結(jié)果和2.5小節(jié)的殘余應(yīng)力測試結(jié)果可知，2#試樣的殘余應(yīng)力峰值分布在其應(yīng)變振型峰值的附近，所以2#試樣具有更好的殘余應(yīng)力均化效果。綜上所述，對于構(gòu)件進(jìn)行振動(dòng)時(shí)效處理時(shí)，需要根據(jù)構(gòu)件表面的殘余應(yīng)力分布狀態(tài)和應(yīng)變振型分布規(guī)律，合理選取振動(dòng)時(shí)效的激振頻率，能夠獲得較為理想的振動(dòng)時(shí)效消除和均化殘余應(yīng)力效果。

(a) x軸向殘余應(yīng)力

(b) y軸向殘余應(yīng)力圖9 振動(dòng)時(shí)效前后2#試樣殘余應(yīng)力測試結(jié)果

Fig.9 Residual stress measuring results of 2# specimen before and after the vibratory stress relief

3 結(jié) 論

(1)通過分析應(yīng)變振型和殘余應(yīng)力的分布規(guī)律，得到了基于應(yīng)變振型的振動(dòng)時(shí)效激振頻率的確定準(zhǔn)則，即為了提高振動(dòng)時(shí)效的效果，對于激振頻率f的選擇應(yīng)以構(gòu)件表面的殘余應(yīng)力分布為依據(jù)，共振應(yīng)變函數(shù)ψr為準(zhǔn)則。

(2)基于應(yīng)變振型的振動(dòng)時(shí)效激振頻率的確定方法，其步驟主要包括：①獲取構(gòu)件表面的殘余應(yīng)力分布狀態(tài)；②建立構(gòu)件的有限元模型，對構(gòu)件進(jìn)行數(shù)值模態(tài)分析；③獲取構(gòu)件的各階應(yīng)變振型和模態(tài)頻率；④在步驟①和③的基礎(chǔ)上，初步確定振動(dòng)時(shí)效的激振頻率；⑤以數(shù)值模態(tài)分析確定的激振頻率為基礎(chǔ)，在其一定范圍內(nèi)對構(gòu)件進(jìn)行掃頻振動(dòng)，獲取振動(dòng)幅值較大的頻率作為振動(dòng)時(shí)效的激振頻率。

(3)1#試樣經(jīng)過振動(dòng)時(shí)效處理后其x軸向和y軸向峰值殘余應(yīng)力的釋放量分別為106.8 MPa和131.8 MPa；2#試樣經(jīng)過振動(dòng)時(shí)效處理后其x軸向和y軸向峰值殘余應(yīng)力的釋放量分別為183.6 MPa和226.1 MPa，表明根據(jù)殘余應(yīng)力分布狀態(tài)和應(yīng)變振型分布規(guī)律，合理的選擇振動(dòng)時(shí)效的激振頻率能夠獲得較為理想的振動(dòng)時(shí)效效果。

參 考 文 獻(xiàn)

[1] WALKER C A, WADDELL A J, JOHNSTON D J. Vibratory stress relief-an investigation of the underlying process[J]. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, 1995, 209:51-58.

[2] KWOFIE S. Plasticity model for simulation, description and evaluation of vibratory stress relief[J]. Materials Science and Engineering A, 2009, 516:154-161.

[3] 蔣剛, 何聞, 鄭建毅. 高頻振動(dòng)時(shí)效的機(jī)理與實(shí)驗(yàn)研究[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2009, 43(7):1269-1272.

JIANG Gang, HE Wen, ZHENG Jiangyi. Mechanism and experimental research on high frequency vibratory stress relief[J]. Journal of Zhejiang University (Engineering Science), 2009, 43(7):1269-1272.

[4] 劉曉丹, 陶興華, 韓振強(qiáng). 振動(dòng)時(shí)效工藝在消除膨脹波紋管殘余應(yīng)力中的應(yīng)用[J]. 振動(dòng)與沖擊, 2015, 34(4):171-174.

LIU Xiaodan, TAO Xinghua, HAN Zhenqiang. Application of vibratory stress relief in relaxation of residual stress for expandable corrugated liners[J]. Journal of Vibration and Shock, 2015, 34(4):171-174.

[5] 韓衍昭, 劉愛敏. 振動(dòng)時(shí)效工藝參數(shù)選取的探討[J]. 熱加工工藝, 2012, 41(16):189-191.

HAN Yanzhao, LIU Aimin. Discussion on selection of technology parameters for vibration stress relief[J]. Hot Working Technology, 2012, 41(16):189-191.

[6] 沈華龍, 吳運(yùn)新, 郭俊康. 高強(qiáng)度鋁合金厚板振動(dòng)時(shí)效工藝的研究[J]. 振動(dòng)與沖擊, 2009, 28(8):191-194.

SHEN Hualong, WU Yunxin, GUO Junkang. VSR technology used in high intensity aluminum alloy thick plates[J]. Journal of Vibration and Shock, 2009, 28(8):191-194.

[7] 李德葆, 陸秋海. 實(shí)驗(yàn)?zāi)B(tài)分析及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2001.

[8] YAM L H, LEUNG T P, LI D B, et al. Theoretical and experimental study of modal strain analysis[J]. Journal of Sound and Vibration, 1996, 191(2):251-260.

[9] 徐麗, 易偉建. 應(yīng)變模態(tài)在框架結(jié)構(gòu)節(jié)點(diǎn)損傷診斷中的應(yīng)用研究[J]. 振動(dòng)、測試與診斷, 2006, 26(2):102-107.

XU Li, YI Weijian. Damage detection of frame structure joints using strain modal method[J]. Journal of Vibration Measurement & Diagnosis, 2006, 26(2):102-107.

[10] HE W, REN Y, CHEN C, et al. Analysis on technology of high-frequency vibratory stress relief[C]∥Proceedings of 2005 ASME Power Conference. Chicago, IL, 2005.

[11] HE Wen, GU Bangping, ZHENG Jianyi, et al. Research on high-frequency vibratory stress relief of small Cr12MoV quenched specimens[J]. Applied Mechanics and Materials, 2012, 157/158:1157-1161.

[12] SHALVANDI M, HOJJAT Y, ABDULLAH A, et al. Influence of ultrasonic stress relief on stainless steel 316 specimens: a comparison with thermal stress relief[J]. Materials and Design, 2013, 46:713-723.

[13] 王壬炎, 王時(shí)英. 超聲振動(dòng)時(shí)效的機(jī)理及實(shí)驗(yàn)研究[J]. 中國農(nóng)機(jī)化學(xué)報(bào), 2015, 36(2):59-62.

WANG Renyan, WANG Shiying. Mechanism analysis and experimental research of the ultrasonic vibration aging[J]. Journal of Chinese Agricultural Mechanization, 2015, 36(2):59-62.

[14] ZIENKIEWICZ O C, TAYLOR R L. 有限元方法基本原理(第1卷)[M]. 曾攀,等,譯. 北京: 清華大學(xué)出版社, 2008.

[15] American Society for Testing and Materials. ASTM E 837-08 Standard test method for determining residual stresses by the hole-drilling strain-gage method[S]. Pennsylvania: ASTM International, 2008.