2. 航空工业西安飞机工业(集团)有限责任公司,西安 710089
2. AVIC Xi’an Aircraft Industry (Group) Company LTD., Xi’an 710089
超声喷丸成形(Ultrasonic peen forming)技术是利用超声发生器产生的强烈冲击波经冲击介质(弹丸或撞针)作用在金属材料表面,具有振动频率高、加载作用强和冲击载荷上升前沿时间短等优点,属于高能高应变速率变形[1],能够在金属表面产生纳米结构层和残余压应力场,可以显著地提高金属零件的机械性能、疲劳寿命和耐腐蚀性能[2-4]。
近年来,国内外学者在相关理论和实验的基础上,对超声喷丸进行了诸多研究。KAWANO S等人[5]为了提高核反应堆零件的表面层性能,利用直径为3 mm的不锈钢弹丸对316 L不锈钢和600合金进行超声喷丸试验。YIN F等人[6]研究了AISI-1018钢经过超声喷丸发生剧烈塑性变形的表面形态,利用非接触式3D轮廓测量仪对喷丸后的工件进行测量,探讨了弹丸直径和喷丸时间对表面形态的影响。YANG X等人[7]采用基于内德勒曼(GTN)的韧性损伤本构方程的有限元模型研究了超声喷丸对AISI 304不锈钢产生的塑性损伤问题,发现随着弹丸冲击速度的增加,受喷材料表面残余压应力值随之增加,且残余压应力层深度也有了明显的增大。史学刚等人[8]研究了超声喷丸覆盖率分别为100%~300%时2024铝合金的疲劳性能,发现喷完产生的表面纳米层和残余压应力使其疲劳寿命提高了5倍以上。
随着航空器件向着大规模、大尺寸、高精度和轻量化发展,需要满足长寿命、低耗能、低维修、高质量的要求,喷丸成形技术存在的的局限性越发明显起来,如工艺参数多且不易控制、工作效率低、易发生球面变形、精确成形能力有待提高,喷完后零件表面粗糙度值有待降低。尽管超声喷丸能够在一定程度上改善以上状况,但是在自由状态下,喷丸成形过程中也易发生球面变形,而预应力喷丸成形能够在一定程度上改变喷丸的球面变形趋势,并且还能够提高零件的成形性能[9]。BARRETT和TODD[10]试验研究了弹性预应力喷丸成形技术对残余应力场分布的影响,发现预应力喷丸成形后的残余压应力较传统自由喷丸成形有了大幅度的提高。MIAO H Y[11]建立了喷丸饱和度、表面覆盖率和粗糙度与喷丸时间之间的定量关系,并且试验研究了喷丸速度和喷丸时间对残余应力的影响以及喷丸过程中预弯力矩对窄板和方板弧高值的影响。XIE L等[12]研究了预应力喷丸对(TiB+TiC)/Ti-6Al-4V复合材料表面性能的影响,通过沃伊特方法得到的衍射峰积分宽度分析了晶粒尺寸和表面变形层微观应变。HU Y等[13]通过在方形试验件的单方向上施加预弯曲率来实现预应力激光喷丸成形,建立了一种基于本征应变的新模型来模拟预应力对弯曲变形和残余应力场分布的影响。
文中结合现行壁板超声喷丸成形和预应力喷丸成形的相关研究,从理论角度分析探讨了预应力喷丸成形后残余应力场的分布以及成形曲率的大小。通过试验研究,从预应力大小、成形轨迹和偏置距离以及壁板厚度等方面阐明预应力超声喷丸的成形特性,为研制开发机翼壁板超声喷丸成形工艺奠定理论与技术基础。
1 预应力超声喷丸成形理论预应力喷丸时,受喷材料表面处于拉应力状态下接受撞针或弹丸的撞击。与自由状态下相比,在相同的撞击能量下,这种拉应力有助于受喷表层材料的延伸,同时加深残余压应力层的深度以及残余压应力值的大小。以在板料宽度方向(弦向)上施加弹性预弯为例,分析喷丸过程的残余应力以及成形曲率半径的大小,如图1所示。
板料的x方向上施加弹性预弯后,如图1(a)所示。板料y方向上的预弯力矩可以忽略不计,x方向单位宽度上的预弯力矩可表示为:
式中:E为弹性模量;I为单位宽度的惯性矩,大小为I=h3/12,h为板料厚度。
板料的弯曲变形会在其内部产生相应的应力场,x方向上由预弯力矩产生的应力大小为:
喷丸过程中,板料上表面处于拉应力状态下接受撞针或弹丸的撞击。喷丸后,板料表层发生剧烈塑性变形并在材料内部诱导产生了应力
喷丸处理后,在解除约束条件之前,在应力
由于施加的是弹性预弯,因而从预弯夹具上取下受喷板料以后,板料受到的预弯力矩和延伸力以及由它们引起的应力有一部分消失,在这里分别引入回弹系数αi和βi,释放约束条件后板料内部的延伸力和弯曲力矩可表示为[16]:
而式(3)中的
忽略
式中:μ表示板料的泊松比;δε表示喷丸后残余压应力深度。
将Rx除以Ry,可得:
由于
由以上理论分析可以看出:预应力超声喷丸成形后,板料内部的残余应力增大,预弯方向上的成形曲率半径减小,与预弯垂直方向上的成形曲率半径增大,这种变化随着预弯曲率(预弯力矩)的增加而加剧。因此,预应力超声喷丸能够有效地克服自由状态喷丸成形时板料的球面变形趋势,并提高其机械性能、疲劳寿命和耐腐蚀性能。
2 预应力超声喷丸成形试验 2.1 超声喷丸设备超声喷丸成形试验采用南京航空航天大学自主研制的数控超声喷丸装置进行,其工作原理是将交流电转换成频率为20 kHz左右的高能超声波,然后通过换能器将高频超声振动经变幅杆转换成振幅为20~50 μm的机械振动,带动撞针将产生的机械能通过撞击传递到工件表面,从而使板料发生塑性变形,图2为超声喷丸装置示意图。
2.2 预弯装置和试验件如图3所示,将试验件紧固在预弯模具的弯曲曲面上,通过控制预弯模具弯曲曲面的曲率半径来实现不同的预弯曲率,从而使得试验件上表层获得不同程度的弯曲应力。由于预应力喷丸成形需要施加的是弹性预弯,因而试验件表面受到的拉应力不能超过材料的屈服极限。结合Von-Mises屈服准则和弹塑性理论,弹性预应力的预弯曲率半径大小必须满足如下不等式:
式中:ν为泊松比;E为弹性模量;h为试验件厚度,mm;σs为材料屈服强度。
试验所用材料为2024-T351铝合金,该材料采用固溶处理加自然时效处理的热处理方式。其化学成分和力学性能分别见表1和表2。
试验件长宽尺寸为:120 mm×50 mm,厚度分别为1.5、2 .0、2.5、3.0和3.5 mm,由方程(10)可计算得到3.5 mm试验件的弹性极限预弯曲率半径为363 mm,为了保证试验板件均为弹性预弯,预弯曲率半径分别取400、600和1 200 mm。
(w/%) | |||||
Element | Si | Fe | Cu | Mn | Mg |
Content | 0.5 | 0.5 | 3.8−4.9 | 0.3−0.9 | 1.2−1.8 |
Element | Cr | Zn | Ti | Al | |
Content | 0.1 | 0.25 | 0.15 | Bal. |
Parameters | Values |
Tensile strength, σb / MPa | 470 |
Yield strength, σs / MPa | 325 |
Elongation, δ/% | 10 |
Fatigue strength, σ−1 / MPa | 105 |
Hardness HB / (N·mm−1) | 120 |
Elastic modulus, E / GPa | 68 |
Density, ρ / (kg· m−3) | 2 770 |
采用控制变量法进行预应力超声喷丸成形试验研究,重点研究预弯力矩,成形轨迹,偏置距离,壁板厚度等参数对2024-T351铝合金成形曲率的影响,偏置距离与成形轨迹示意图如图4。在试验件的中心选取尺寸为80 mm×50 mm作为喷丸区域,选用单撞针工具头,3 mm直径的撞针、变幅杆振幅调为24 μm,机床进给速度为3 000 mm/min。
具体试验研究的安排如下:
(1) 成形轨迹是影响超声喷丸成形的重要因素,它对板料成形的外形轮廓具有非常重要的影响。采用图4所示两种成形轨迹,在预弯曲率半径Rp为+∞、1 200、600、400 mm的模具上,偏置距离d=1 mm,对厚度h为2.5 mm的试验件进行超声喷丸成形试验,研究弹性预应力对超声喷丸成形曲率的影响以及在有无预应力状态下成形轨迹对成形曲率的影响。
(2) Rp为+∞、1 200、600、400 mm,选取成形效果较优的成形轨迹,拟定偏置距离d分别为0.6、0.8、1、1.2和1.4 mm,对h为2.5 mm的试验件进行超声喷丸试验,研究预应力状态下偏置距离对超声喷丸成形性能的影响。
(3) 采用同样的成形轨迹,d=1 mm,对h为1.5、2、2.5、3和3.5 mm,在Rp为+∞、1 200、600和400 mm的条件下进行超声喷丸成形试验,研究壁板厚度与预应力超声喷丸成形之间的作用规律。
为了进一步研究弹性预应力对超声喷丸成形综合性能的影响,对在不同预弯曲率半径下喷丸后的试验件进行残余应力和表面粗糙度测量。残余应力测量设备采用Proto公司制造的腐蚀设备和残余应力测试仪,采用同倾固定Ψ加摆动法测量,Cr靶材kα辐射,测试晶面为(311),实验前对材料表面进行5~10 s的腐蚀,去除表面氧化膜,以确保试验测量的准确性。试验件表面粗糙度通过Surftest SJ-310进行测量,探针方向垂直喷丸轨迹,每个试验件测量8个点。
3 结果与讨论 3.1 弧高值和成形曲率超声喷丸试验后壁板成形曲率半径R的结果通过采用测量弧高值AHp然后根据其几何关系计算得出,测量所用弧高仪如图5所示,其测头部位是由上下移动的测针和位移传感器连接组成,通过表面光滑平整的校准块校零,两侧滑块可以调整跨距。成形曲率半径的表达式如下:
由于AHp远小于R,所以式(12)可以简写成:
式中,L表示测量跨距,弦向跨距为50 mm,展向跨距为100 mm。
表3列出了4种预弯曲率半径下的预弯力矩、表面最大拉应力以及采用图4(b)所示的弦向成形轨迹(IACD)对厚度为2.5 mm壁板进行超声喷丸后弦向和展向的弧高值平均值以及对应的成形曲率半径。预弯力矩和表面最大拉应力可分别由式(1)和式(2)计算得出。从表中数据可以看出,随着预弯曲率半径的减小,壁板展向的弧高值AHx随之增大,预弯曲率半径为400 mm时,AHx为2.754 mm,较自由状态而言,弧高值增大了150.6%;而弦向的弧高值AHy随着预弯曲率半径的减小而减小,预弯曲率半径为400 mm时的弧高值比自由状态减小了38.1%。
Radius of prebend, Rp / mm | Prebending moment, Mpre / (N·mm) | Surface stress, σst / MPa | Resulting arc height / mm | Forming radius of curvature / mm | ||
Span wise | Chord wise | Span wise | Chord wise | |||
+∞ | 0 | 0 | 1.099 | 0.126 | 1 137 | 2 470 |
1 200 | 73.8 | 70.8 | 1587 | 0.106 | 788 | 2 947 |
600 | 147.6 | 141.6 | 2.219 | 0.088 | 563 | 3 551 |
400 | 221.4 | 212.4 | 2.754 | 0.078 | 454 | 4 020 |
将表3中的弧高值和成形曲率半径数据转化为曲线,如图6所示。从图6(a)可以看出:壁板展向(弦向)弧高值随着预弯力矩的增加而增加(减小),且几乎呈线性变化。由图6(b)可知,预弯力矩Mpre与成形曲率半径R之间的关系和其与弧高值的关系不尽相同。随预弯力矩的增大,喷丸后壁板展向成形曲率半径Rx先快速减小后逐渐趋于平缓,而弦向成形曲率半径Ry随预弯力矩的增大几乎呈线性增大。当预弯曲率Rp为400 mm时,Rx为454 mm,Rp为800时,Rx为563 mm,Rp为1 200 mm时,Rx为788 mm;相比自由状态成形曲率半径1 137 mm,成形曲率半径分别减少了60.1%、50.5%、30.7%。而壁板弦向上的成形曲率Ry远远大于Rx,并随着预弯力矩的增加而线性增大,预弯曲率半径Rp为400 mm时,与自由状态相比其成形曲率半径增大了62.8%。
就同一试验件而言,其弦向和展向成形曲率半径的比值随着预弯曲率的增加而增大,当预弯曲率为0时,Ry与Rx的比值为2.17,随着预弯曲率的增加,Ry与Rx的比值最大达到了8.82,表明了预应力超声喷丸能够大幅度减小球面变形趋势。
成形轨迹是超声喷丸的重要参数,对板料成形的外形轮廓具有非常重要的影响。图7对比分析了图4两种成形轨迹,在不同预弯条件下对厚度为2.5 mm壁板进行超声喷丸后得到的成形曲率。从图中可以看出,两种成形轨迹得到的Rx和Ry随预弯力矩的变化趋势基本一致,即Rx随预弯力矩的增大先快速而后平缓减小,Ry随预弯力矩的增大近乎线性增大。另一方面,弦向喷丸(IACD)得到的Rx比展向喷丸(IASD)要小,而Ry比展向喷丸要大。显然,选择弦向成形轨迹对壁板进行预应力超声喷丸不仅能够进一步提高其成形曲率,还能更好地改善其球面变形趋势。
偏置距离对于板料的成形量以及诱导产生的残余应力场的大小和分布具有非常重要的影响。上述分析结果表明,IACD较IASD能进一步提高壁板的成形曲率半径,并在更大程度上消除球面变形的趋势。因此,在IACD的基础上,分析了偏置距离d对壁板展向成形曲率的影响,结果如图8所示。可以看出,Rx随偏置距离的减小而减小,这是因为随着偏置距离的减小,壁板单位面积受到撞针冲击的次数越多,传递到壁板中的冲击能量也就越多,从而使得Rx越来越小。当Rp=+∞,d=1.4 mm时,Rx为1 482 mm,而在Rp=400 mm,d=0.6 mm时,Rx为347 mm,较Rp=+∞减小了76.6%。
可见,在预应力状态下减小偏置距离可以大幅度减小壁板的成形曲率半径。这可能是由于预应力的存在,超声喷丸后壁板内部诱导产生了更大的残余压应力,而偏置距离的减小使得残余压应力进一步加大,导致壁板的弯曲变形程度加大。
图9分析了不同厚度壁板在预弯曲状态下超声喷丸成形后的展向成形曲率半径。由图9可知,无论是在自由状态下还是预弯状态下,成形曲率半径都随着壁板厚度的增加而增加,且增加幅度越来越大。自由状态下,2.5和3.5 mm壁板的成形曲率半径由1.5 mm壁板的897 mm分别增大到1 141和2 070 mm,分别增大了27.2%和130.7%;预弯曲率半径为400 mm时,2.5和3.5 mm壁板较1.5 mm的成形曲率半径由387 mm增大为455 mm和804 mm,分别增加17.6%和107.8%。可见,预应力状态下,壁板成形曲率半径随其厚度增加的增加幅度要比自由状态小,且随着预弯曲率半径的减小,增加幅度越来越小。与自由状态相比较而言,预弯曲率半径为400 mm时,1.5 mm壁板的成形曲率半径减少了56.8%,2.5和3.5 mm分别减少了60.1%和61.1%。可见预应力状态与自由状态相比,壁板成形曲率半径的减少幅度并不随着其厚度的增加而减小。因此,利用预应力超声喷丸成形工艺可以提高更大厚度零件的成形曲率。
3.2 残余应力材料表面层的残余压应力场不仅能够有效的阻止材料内裂纹的萌生,还能通过残余压应力场的作用增加裂纹的闭合效应减小疲劳裂纹的扩展。超声喷丸利用在材料内部引入残余压应力场来改善整个工件的应力分布状况达到提高疲劳寿命的目的,已有研究表明,引入的残余压应力场越大,试样的疲劳寿命越高[17]。为更好地描述喷丸后残余压应力场的特征,将其归纳为5个特征参量,即上表面残余压应力值σtsrs、下表面残余压应力值σbsrs最大残余压应力σmrs、最大残余压应力值深度Zm和残余压应力层深度Z0。
图10比较了厚度为2.5 mm的壁板在不同预弯力矩下超声喷丸后上表面和下表面的残余应力值的大小,可以看出:无论是否处于预应力状态进行超声喷丸成形,壁板的上下表面均产生了残余压应力,且σtsrs>σbsrs;随着预弯力矩的增大,σtsrs、σbsrs均随之呈现出近似线性增大的趋势,该结果与HU Y X[13]等人得出的结论基本一致。
图11分析了Rp=400、600、1 200 mm以及+∞条件下超声喷丸后2.5 mm壁板残余应力场分布。从图可以看出,与自由相比较,预应力喷丸在壁板内部诱导产生的σtsrs、σmrs、Zm以及Z0均更大,而且预应力喷丸在壁板内部产生的残余拉应力也比自由状态更大。随着预弯曲率半径减小,壁板的上表面残余压应力值(σtsrs)随之增大,预弯曲率半径为400 mm的σtsrs(229 MPa)较自由状态(153 MPa)增大了49.7%,σmrs、Zm以及Z0也随着预弯曲率半径的减小而增大。可见,预应力超声喷丸能够进一步提高受喷零件的残余压应力场,从而提高受喷零件的抗疲劳性能。
3.3 表面粗糙度表面粗糙度直接影响着材料零部件的耐腐蚀性、耐磨性、疲劳强度及导流、导电等性能。在传统的喷丸成形中,弹丸流高速撞击材料的表面,产生大量的弹坑,使得受喷材料的表面非常粗糙,严重影响其表面质量。因此,对预应力超声喷丸成形后板料表面粗糙度的研究具有非常重要的意义。表4为不同预弯条件下,原材料以及受喷件表面粗糙度主要参数Ra的测量值。可以看出:超声喷丸后,试验件的表面粗糙度Ra的平均值有所增加,由原材料的0.562 μm增加至1.903 μm,且随着预弯曲率半径的减小,Ra平均值先略微增加然后减小,这是因为预弯产生的弹性变形使得喷丸过程中发生的塑性变形更加充分和均匀,在试样件表面形成的更低的波峰,导致粗糙度值略微下降。可见,预应力超声喷丸成形不仅产生了更大的残余压应力场,提高了受喷零件的成形曲率,而且并没有使其表面粗糙度值变大。
但是,相对于原材料而言,预应力超声喷丸后零件的表面粗糙度值有了大幅提高,而表面粗糙度的增大,往往会造成零件表面的损伤,应力集中也会明显增加,对疲劳性能非常不利。在实际应用中,为了降低喷丸成形件的表面粗糙度值,通常采用打磨进行处理,打磨后喷丸所产生弹坑的波峰大幅度降低,从而实现减小表面粗糙度值的效果。此外,对成形件采用小压力、大弹丸的工艺参数进行喷丸强化不仅可以减小成形件的表面粗糙度值,还可以进一步提高其表面显微硬度和残余应力场。
Radius of prebend, Rp / mm | Surface roughness, Ra / μm | ||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | Average | |
+∞ | 1.029 | 1.153 | 1.066 | 0.986 | 1.095 | 1.092 | 1.161 | 1.158 | 1.093 |
1 200 | 1.124 | 1.146 | 1.153 | 1.209 | 1.056 | 1.193 | 0.997 | 1.161 | 1.130 |
600 | 1.035 | 0.892 | 1.019 | 0.934 | 0.995 | 1.037 | 1.109 | 1.114 | 1.017 |
400 | 0.979 | 0.954 | 1.024 | 0.923 | 0.823 | 0.792 | 1.04 | 0.892 | 0.928 |
Raw materials | 0.587 | 0.509 | 0.624 | 0.563 | 0.517 | 0.569 | 0.521 | 0.605 | 0.562 |
文中结合预应力喷丸成形和超声喷丸成形两种工艺方法,对2024-T351铝合金壁板采用试验研究分析了不同工艺参数对成形性能的影响。基于分析结果,得出了以下结论:
(1) 无论是否处于预应力状态,超声喷丸成形壁板的展向成形曲率半径远小于弦向;且随着预弯力矩的增大,展向成形曲率半径随之减小,而弦向成形曲率半径随之增大,大幅度地减小了壁板的球面变形趋势。
(2) 喷丸成形轨迹和偏置距离是影响预应力超声喷丸的重要因素。较展向成形轨迹而言,弦向成形轨迹能够进一步提高(降低)其展向(弦向)的成形曲率,减小球面变形趋势。随着偏置距离的减小,壁板展向的成形曲率半径随之减小,且其减小幅度随着预弯曲率半径的减小而减小。
(3) 预应力超声喷丸能够实现更大厚度的板材成形,且与自由状态相比较而言,预应力状态下壁板成形曲率半径的减小幅度并不随其厚度的增加而减小。
(4) 超声喷丸成形在壁板的上下表面均产出了残余压应力,且上表面的残余压应力值较下表面的大。随着预弯曲率半径的减小,表面残余压应力值(σtsrs、σbsrs)、最大残余压应力值(σmrs)、最大残余压应力值深度(Zm)、残余压应力深度(Z0)均随之增大。
(5) 预应力超声喷丸成形不仅能够诱导产生更大更深的残余应力场,提高其成形曲率,且成形后试验件表面的粗糙度并未随着预弯曲率半径的减小而增加,表明了预应力超声喷丸能够获得更好的表面层性能。
[1] | SEKKAL A C, LANGLADE C, VANNES A B. A micro/macro impact test at controlled energy for erosion and phase-transformation simulation[J]. Tribology Letters, 2003, 15(3): 265-274. |
点击浏览原文 | |
[2] | MORDYUK B N,PROKOPENKO G I. Ultrasonic impact peening for the surface properties’ management[J]. Journal of Sound & Vibration, 2007, 308(3): 855-866. |
[3] | FAN Z, XU H, LI D, et al. Surface nano crystallization of 35# type carbon steel induced by ultrasonic impact treatment (UIT)[J]. Procedia Engineering, 2012, 27: 1718-1722. |
点击浏览原文 | |
[4] | HACINI L, LE N V, BOCHER P. Evaluation of Residual Stresses Induced by Robotized Hammer Peening by the Contour Method[J]. Experimental Mechanics, 2009, 49(6): 775-783. |
点击浏览原文 | |
[5] | KAWANO S, KAWAGISHI A, SUEZONO N, et al. Development of ultrasonic shot peening technique for reactor components to improve structural integrity against stress corrosion cracking[C]. International Conference on Nuclear Engineering, 2009: 543-547. |
[6] | YIN F, HUA L, WANG X, et al. Numerical modeling and experimental approach for surface morphology evaluation during ultrasonic shot peening[J]. Computational Materials Science, 2014, 92(5): 28-35. |
[7] | YANG X, ZHOU J, LING X. Study on plastic damage of AISI 304 stainless steel induced by ultrasonic impact treatment[J]. Materials & Design, 2012, 36: 477-481. |
[8] |
史学刚, 鲁世红, 张炜. 铝合金超声波喷丸成形制件表面完整性研究[J]. 中国机械工程, 2013, 24(22): 3100-3104.
SHI X G, LU S H, ZHANG W. Study on surface integrity of aluminum alloy ultrasonic shot peening forming workpieces[J].China Mechanical Engineering, 2013, 24(22): 3100-3104(in Chinese). |
点击浏览原文 | |
[9] | HACKEL L A, HALPIN J M, HARRIS F B. Pre-loading of components during laser peenforming: EP, US 6670578 B2[P]. 2002. |
[10] | BARRETT C F, TODD R. Investigation of the effects of elastic pre-stressing technique on magnitude of compressive residual stress induced by shot peen forming of thick aluminum plates[C]. Proceedings of the 2nd International Conferences on shot peening,1984: 15-21. |
[11] | MIAO H Y, DEMERS D, LAROSE S, et al. Experimental study of shot peening and stress peen forming[J]. Journal of Materials Processing Technology, 2010, 210(15): 2089-2102. |
点击浏览原文 | |
[12] | XIE L, JIANG C, LU W, et al. Effect of stress peening on surface layer characteristics of (TiB+TiC)/Ti–6Al–4V composite[J]. Materials & Design, 2012, 33: 64-68. |
[13] | HU Y, LI Z, YU X, et al. Effect of elastic prestress on the laser peen forming of aluminum alloy 2024-T351: Experiments and eigenstrain-based modeling[J]. Journal of Materials Processing Technology, 2015, 221: 214-224. |
点击浏览原文 | |
[14] | AL-HASSANI S T S. Mechanical aspects of residual stress development in shot peening[J]. Hasscmi, 1981. |
[15] | MIAO H Y, LAROSE S, PERRON C, et al. Numerical simulation of the stress peen forming process and experimental validation[J]. Advances in Engineering Software, 2011, 42(11): 963-975. |
点击浏览原文 | |
[16] | HASHMI S, BATALHA G F, VAN TYNE C J, et al. Comprehensive materials processing[M]. Elsevier, 2014: 325-326. |
[17] | ZÁVODSKÁ D, GUAGLIANO M, BOKŮVKA O, et al. Fatigue resistance of low alloy steel after shot peening[J]. Materials Today Proceedings, 2016, 3(4):1220-1225. |
点击浏览原文 |