激光强化后7050凹槽铝板的残余应力分布规律研究
来源:http://www.secondhandfilm.com??发布时间:2019-11-26 09:08??浏览量:返回列表
因具有高比强度、抗腐蚀性和疲劳性能, 7050-T7451铝合金是飞机结构主要用材之一。通常采取表面强化技术 (如机械喷丸强化、激光冲击强化等) 对飞机结构表面进行处理, 以提高其疲劳寿命和安全可靠性。与传统机械喷丸相比, 激光冲击能够获得较高的表层残余应力和较深的残余应力影响层
国内外学者对激光冲击强化技术进行了较为广泛的研究, 通过试验和数值模拟研究了激光冲击强化对残余应力分布的影响。激光冲击强化过程中, 激光参数对残余应力分布的影响极大, 虽然光斑中心的冲击压力较高, 但通过激光冲击试验发现, 单次冲击时光斑中心的残余压应力低于光斑边缘处所形成的残余压应力, 即形成“残余应力洞”现象, 当增加冲击次数时能够削弱这一现象并获得更加均匀的表面残余压应力
本文以7050-T7451凹槽铝板为研究对象, 测量其激光冲击强化残余应力场分布规律, 采用ABAQUS软件建立三维有限元仿真模型, 利用VDLOAD子程序施加冲击波压力载荷, 分析7050-T7451凹槽铝板在不同激光功率密度, 光斑直径和激光脉宽等条件下残余应力场的分布规律, 为飞机结构设计工作提供理论依据。
1 激光冲击强化试验
1.1 材料性能和试验件尺寸
本文中使用的试验件尺寸为80 mm×120 mm, 凹槽区域 (如图1所示) 尺寸为40 mm×20 mm, 对凹槽加工后进行处理, 以消除凹槽加工过程中产生的残余应力。激光冲击强化区域大于凹槽区域, 尺寸为48 mm×28 mm。材料是7050铝合金, 其主要化学成分见表1, 材料参数见表2。
表1 7050铝合金化学成分 导出到EXCEL
质量分数/%
Zn |
Mg | Cu | Cr | Zr | Si | Fe | Al |
5.7~6.7 |
1.9~2.6 | 2.0~2.6 | ≤0.04 | 0.08~0.15 | ≤0.12 | 0.15 | 其他 |
表2 7050铝合金材料参数 导出到EXCEL
密度/ (kg·m-3) |
弹性参数/MPa | 泊松比 | 屈服应力/MPa |
2 830 |
72 000 | 0.33 | 395 |

图1 试验件及其尺寸示意图 下载原图
1.2 激光冲击路径和光斑搭接情况
试验采用YS80-M165激光冲击强化设备, 为确保残余应力分布均匀, 经过多次测试, 选择X向搭接率为25%, Y向搭接率为35%, 冲击次数为2次, 波长为1 064 nm, 激光脉宽为20 ns, 激光束光斑直径为4 mm。激光冲击路径和光斑搭接情况如图2所示。

图2 激光冲击路径和光斑搭接示意图 下载原图
2 激光冲击强化后残余应力场分布数值分析
2.1 冲击波峰值压力
在数值分析中将激光诱导的冲击波简化为随时间变化的压力载荷, 然后直接作用于靶材表面。Fabbro
式中, α为内能与热能比值, 取0.1~0.2之间;Z为材料的折合声阻抗, 可由 (2) 式得到;I0为激光功率密度, 可由 (3) 式得到。
式中, Z1为靶材的声阻抗;Z2为约束层的声阻抗。
激光功率密度计算公式为
式中, γ为等离子体的绝热指数, 通常取0.7;τ为激光脉宽。
将 (2) 式和 (3) 式的计算结果代入 (1) 式可得激光冲击强化峰值压力为
2.2 有限元数值分析模型
为提高计算机运行速率、减少计算量, 建模时取试验件1/4进行建模, 如图3所示。在激光冲击强化区域将网格细化, 共有343 688个单元。模型的2个侧面设置固定边界条件, 模型的剖面设置为对称边界条件。

图3 有限元模型 下载原图
2.3 本构模型和VDLOAD子程序加载
激光冲击过程中, 激光诱导的等离子体冲击一般为几千MPa, 由于吸收涂层的气化时间极短, 因而在模拟中可完全忽略激光冲击时的热效应, 仅考虑材料在激光冲击波力效应下的动态响应
式中, A, B, C和n为材料常数, 具体数值如表3所示。
表3 Johnson-Cook模型材料参数
材料 |
A/MPa | B/MPa | C | n |
7050铝合金 |
300 | 400 | 0.013 | 0.13 |
为了实现激光强化过程中多点连续冲击压力的施加, 引入ABAQUS的VDLOAD子程序, 流程图如图4所示。

图4 VDLOAD加载流程图 下载原图
3 结果与讨论
3.1 试验与数值分析的表面残余应力对比
采用X射线衍射仪对冲击区的表面残余应力场进行测试, 试验测量误差为±17 MPa, 图5为试验测量和数值分析的激光冲击后表面残余应力 (σx) 分布对比。由试验测量值可以看出, 在原点 (凹槽中心点) 处残余压应力为256 MPa, 距原点8.5 mm, 9.5 mm, 10.5 mm, 12.0 mm, 13.0 mm和14.0 mm处 (沿表面方向如图3所示) 残余压应力分别为306 MPa, 268 MPa, 261 MPa, 207 MPa, 203 MPa和60 MPa, 与原点处相比残余压应力分别相差19.5%, 4.7%, 2.0%, -19.1%, -20.7%和-76.5%??梢钥闯? 激光冲击后表面残余压应力在凹槽的曲面处取得最大值306 MPa, 在冲击区域的边缘处表面残余压应力取得最小值60 MPa。

图5 表面残余应力的试验值和数值模拟值对比 下载原图
数值分析和试验测量的残余应力分布基本规律相同:从试验件凹槽中心到凹槽边缘, 残余应力逐渐下降, 在凹槽边缘底边处取最小值, 然后沿凹槽边缘上升而增大, 表面残余压应力在凹槽曲面底边处取最大值。实测表面残余应力最小值为-306 MPa, 数值模拟的表面残余应力最小值为-338 MPa, 两者相差10.4%;凹槽中心处的表面残余应力试验值为-256 MPa, 模拟值为-214 MPa, 两者相差16.4%;随着凹槽曲面上升, 压残余应力减小;远离冲击区域, 表面残余应力逐渐接近于零。
3.2 激光功率密度对残余应力场的影响
为研究激光功率密度对残余应力场的影响, 分别取0.84 GW/cm2, 1.95 GW/cm2, 3.06 GW/cm2, 4.18 GW/cm2和5.29 GW/cm2这5种激光功率密度, 其对应的峰值压力为1 052 MPa, 1 607 MPa, 2 015 MPa, 2352 MPa和2374 MPa, 其他参数不变进行数值分析。图6为不同功率密度激光冲击强化后沿表面方向和深度方向残余应力的分布情况 (其中沿表面和深度方向的起始点均为凹槽中心点) 。

图6 不同激光功率密度下残余应力场分布 下载原图
由图6a) 可以看出, 当激光功率密度取0.84, 1.95, 3.06, 4.18, 5.29 GW/cm2时, 凹槽中心点处表面残余压应力分别为32, 214, 230, 210, 195 MPa, 相邻功率密度间残余压应力增加的幅度分别为182, 16, -20, -15 MPa。随着激光功率密度的增加, 凹槽中心点处表面残余压应力先增加后减小, 当激光功率密度为3.06 GW/cm2时取得最大值230 MPa。这一现象的存在是由于材料中应力波的峰值压力小于Hugoniot弹性极限
3.3 光斑直径对残余应力场的影响
图7为不同光斑直径激光冲击强化后沿表面方向和深度方向残余应力场分布。沿表面方向的残余应力只对冲击区域内进行分析。

图7 不同光斑直径下残余应力场分布 下载原图
由图7a) 可以看出, 随着光斑直径从4 mm增加到6 mm, 凹槽中心点处表面残余压应力由214 MPa减小至30 MPa, 在凹槽的曲面处, 表面残余压应力取得最大值, 在冲击区域边缘处残余压应力逐渐减小。由图7b) 可以看出, 随着光斑直径的增加, 不仅最大残余压应力随之减小, 而且残余压应力影响层也逐渐减小。这种现象是因为残余压应力与冲击波峰值压力相关, 而光斑直径与冲击波峰值压力成反比, 如公式 (4) 所示, 当光斑直径从4 mm增至6 mm时, 对应的峰值压力从1 607 MPa减小到1 072 MPa。
3.4 激光脉宽对残余应力场的影响
图8为不同激光脉宽下激光冲击强化后沿表面方向和深度方向的残余应力场分布。

图8 不同激光脉宽下残余应力场分布 下载原图
由图8a) 可以看出, 激光脉宽分别为10, 20, 30, 40 ns时, 凹槽中心点处表面残余压应力分别为21, 214, 241, 288 MPa。这是因为与冲击波穿透约束层的时间间隔相比, 随着激光脉冲宽度的增加, 约束层受激光冲击后汽化物的扩散延长了激光冲击波的作用时间
4 结 论
1) 不论激光参数如何变化, 凹槽铝板的表面残余应力最小值都在凹槽曲面底边处取得, 凹槽曲率的变化影响着表面残余应力场分布。
2) 随着激光功率密度的增加, 冲击波峰值压力随之增加, 凹槽中心点处表面残余压应力先增加后减小, 当激光功率密度为3.06 GW/cm2时, 峰值压力为2 015 MPa, 此时凹槽中心点处表面残余压应力达到最大值230 MPa。
3) 光斑直径与冲击波峰值压力成反比, 随着光斑直径从4 mm增至6 mm, 峰值压力由1 607 MPa减小到1 072 MPa, 凹槽中心点处表面残余压应力由214 MPa减小至30 MPa。
4) 激光脉冲宽度对表面残余应力的影响与激光冲击波穿透汽化物层的时间间隔有关。随着激光脉宽从10 ns增至40 ns, 凹槽中心点处表面残余压应力由21 MPa增至288 MPa。
5) 对7050凹槽铝板激光冲击强化时, 激光功率密度、光斑直径和脉冲宽度分别取3.06 GW/cm2, 4 mm和40 ns时, 能够在凹槽表面产生最大的残余压应力。