本文以铝合金试样作为研究对象,采用有限元仿真软件ABAQUS与裂纹扩展软件FRANC3D相结合的方法,进行铝合金的疲劳裂纹扩展试验仿真。首先,通过疲劳裂纹扩展试验及光学显微镜的裂纹形貌研究,获得疲劳裂纹扩展次数及裂纹扩展形貌;其次在有限元仿真平台ABAQUS中建立试验静力学模型;最后基于裂纹扩展软件FRANC3D预制初始裂纹,采用奇异性网格划分,并采用最大周向力准则进行裂纹扩展仿真,以及借助Paris公式计算剩余寿命。结果发现:基于ABAQUS和FRANC3D的铝合金裂纹扩展试验仿真与现实裂纹扩展试验结果有很好的一致性;并通过计算应力强度因子判别此裂纹扩展为I型(张开型)裂纹,从而得到裂纹更清晰有效的演化过程。<br/>Taking aluminum alloy specimens as the research object, the method of combining the finite ele-ment simulation software ABAQUS and the crack growth software FRANC3D is used to simulate the fatigue crack growth test of the aluminum alloy. First, the fatigue crack growth test and optical mi-croscope crack morphology research are conducted to obtain fatigue crack growth times and crack growth morphology; secondly, the experimental static model is established in the finite element simulation platform ABAQUS; finally, the initial crack is prefabricated based on the crack growth software FRANC3D, using singularity grid division, using the maximum circumferential force criterion for crack propagation simulation, and calculating the remaining life with the aid of the Paris formula. The results found that: the aluminum alloy crack growth test simulation based on ABAQUS and FRANC3D is in good agreement with the actual crack growth test results; and the crack growth is judged to be a type I (open type) crack by calculating the stress intensity factor, so as to obtain a clearer and more effective evolution process of the crack.
本文以铝合金试样作为研究对象,采用有限元仿真软件ABAQUS与裂纹扩展软件FRANC3D相结合的方法,进行铝合金的疲劳裂纹扩展试验仿真。首先,通过疲劳裂纹扩展试验及光学显微镜的裂纹形貌研究,获得疲劳裂纹扩展次数及裂纹扩展形貌;其次在有限元仿真平台ABAQUS中建立试验静力学模型;最后基于裂纹扩展软件FRANC3D预制初始裂纹,采用奇异性网格划分,并采用最大周向力准则进行裂纹扩展仿真,以及借助Paris公式计算剩余寿命。结果发现:基于ABAQUS和FRANC3D的铝合金裂纹扩展试验仿真与现实裂纹扩展试验结果有很好的一致性;并通过计算应力强度因子判别此裂纹扩展为I型(张开型)裂纹,从而得到裂纹更清晰有效的演化过程。
铝合金试验,ABAQUS,FRANC3D,应力强度因子,裂纹扩展特性
Jun Ma, Ying Cui
School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai
Received: Sep. 13th, 2021; accepted: Oct. 28th, 2021; published: Nov. 4th, 2021
Taking aluminum alloy specimens as the research object, the method of combining the finite element simulation software ABAQUS and the crack growth software FRANC3D is used to simulate the fatigue crack growth test of the aluminum alloy. First, the fatigue crack growth test and optical microscope crack morphology research are conducted to obtain fatigue crack growth times and crack growth morphology; secondly, the experimental static model is established in the finite element simulation platform ABAQUS; finally, the initial crack is prefabricated based on the crack growth software FRANC3D, using singularity grid division, using the maximum circumferential force criterion for crack propagation simulation, and calculating the remaining life with the aid of the Paris formula. The results found that: the aluminum alloy crack growth test simulation based on ABAQUS and FRANC3D is in good agreement with the actual crack growth test results; and the crack growth is judged to be a type I (open type) crack by calculating the stress intensity factor, so as to obtain a clearer and more effective evolution process of the crack.
Keywords:Aluminum Alloy Test, ABAQUS, FRANC3D, Stress Intensity Factor, Crack Propagation Characteristic
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着经济的发展,铝合金材料凭借它良好的可塑性、易加工等特点,在各行各业中得到了日益广泛的应用。探测铝合金材料裂纹扩展,找到在疲劳载荷作用下裂纹扩展过程的参数特性以及损伤断口的微观特性 [
目前国内外许多学者在铝合金疲劳裂纹试验研究方面开展了大量的工作,也提出了许多疲劳裂纹扩展模型来体现裂纹扩展过程中的分散性 [
试验样品为2219铝合金,其抗拉强度 σ b 为450 Mpa,屈服强度 σ s 为280 Mpa。试件采用杠铃式平板设计,根据美国试验材料学会(ASTM E8/E8M-15a)标准和金属材料拉伸试验室温试验方法(GB/T228. 1-2010)等符合国家标准的情况下进行制备。试件形状及尺寸如图1所示,试件厚度为4.5 mm。
图1. 试样的几何尺寸图
疲劳裂纹扩展试验在LETER疲劳试验机(如图2)上进行,试验载荷谱为1~10 kn的循环轴向加载,加载频率为30 Hz。如图3,设置了4组试验,在没有明显的试验误差下,记录相关参数平均值,获得铝合金试样的疲劳裂纹扩展次数为23,303次。
图2. 疲劳试验机
图3. 铝合金试样
在铝合金试样的断口附近取样,采用如图4所示的Quanta200环境扫描光镜,观察疲劳试样断口的微观形貌。从微观形貌看,疲劳断裂样品的断口分为裂纹源区、疲劳裂纹扩展区及快速断裂区,裂纹自裂纹源处呈放射状向芯部扩展,初始裂纹源为直径1 mm的半圆弧裂纹。
图4. 环境扫描光镜
图5为铝合金裂纹扩展模型的工作流程图,主要包括两大部分:ABAQUS有限元建模和FRANC3D裂纹扩展分析 [
图5. 疲劳裂纹模型工作流程图
采用ABAQUS有限元仿真软件,按照现场调研的试验方式,建立1:1的铝合金试样的有限元仿真模型,模型参数如表1所示。
模型名称 | 泊松比 | 弹性模量(MPa) | 密度(kg/m3) | 阻尼比 |
---|---|---|---|---|
铝合金试样 | 0.28 | 7.7 × 104 | 2.71 × 103 | 0.045 |
表1. 有限元模型主要结构参数
在模型的建立过程中,铝合金试样使用C3D8R单元建立,是因为采用8节点六面体线性减缩积分单元时,若裂纹的引入应起网格的扭曲变形,这种单元划分仍能保证分析的精度不会因此受到影响。
为模拟铝合金裂纹扩展试验的全过程,在load中赋予一端固定,一端10 kn的垂向载荷,有限元静力学模型如图6所示。
由图7静力学分析结果可知,应力的最大值出现在试样中部表面,最大值为27.78 Mpa,此位置最易产生裂纹,读取点坐标的位置作为初始裂纹萌生点。
为适应计算收敛需求,利用FRANC3D裂纹扩展分析软件,选取上图应力红色区域作为子模型并在应力最大值处插入如图8(a)所示半圆弧状裂纹,裂纹半径为0.5 mm,裂纹面平行于xoz平面,网格为四面体网格。
图6. 有限元静力学模型
图7. 静力学分析结果图
图8. (a) 裂纹面模型图;(b) 裂纹前缘图;(c) 试样裂纹位置及网格分布图
如图8(b),在插入裂纹后,裂纹前缘有明显的应力集中现象,应力值也显著提高,最高可达270.8 Mpa。由图8(c)所示含裂纹的子模型与全局模型的连接方式为Merger nodes连接,裂纹在应力比为0.1的循环载荷作用下发生裂纹扩展。裂纹在扩展过程中,该子模型的网格会不停重新划分,来满足裂纹的计算收敛性以及应力强度因子计算的精确性。
在完成初始裂纹的预制后,选用经典的Paris公式模型: d a / d n = c ( Δ k ) n ,建立试样的疲劳裂纹扩展速率模型。获取铝合金材料的断裂韧度Kic,应力强度因子门槛值∆Kth及疲劳参数n。以2219铝合金材料为例,查阅相关文献 [
给定裂纹扩展长度a = 0.2 mm为一个扩展步,读取ABAQUS里全局模型的工况,扩展到k值大于或等于断裂韧度Kic为止。如图9可知,0.5 mm的初始裂纹,在1 kn~10 kn~1 kn的循环载荷下,经过了11个扩展步,铝合金的裂纹扩展寿命为22,846次,对应的裂纹扩展长度为2.2 mm。从裂纹扩展速率 d a / d n 来看,也就是图中曲线的斜率来看,随着裂纹扩展的越长,曲线越陡,斜率越大,对应的裂纹扩展也就越快。
扩展速率模型 | 参数c | 参数n | 门槛值∆Kth (Mpa*√m) | 断裂韧度Kic (Mpa*√m) |
---|---|---|---|---|
Paris公式 | 1.29e−7 | 2.8 | 2.9 | 74.52 |
表2. 疲劳裂纹扩展Paris参数
图9. 疲劳裂纹扩展a-n曲线图
图10左为疲劳试样断口的微观形貌图,参照比例尺,明亮的半椭圆区域为初始疲劳裂纹:深度为0.5 mm,在裂纹扩展失稳阶段观察疲劳裂纹深度为2.9 mm,扩展了约为2.4 mm,疲劳扩展次数为23,303次。图10右为仿真模型的裂纹扩展面图,仿真的疲劳裂纹深度为2.7 mm,扩展了疲劳裂纹的长度为2.2 mm,计算的疲劳裂纹扩展次数为22,846次。并且对此发现,试验图的裂纹面和仿真结果一样都有明显的半圆弧水波状裂纹辉纹,是典型的疲劳载荷裂纹。因此仿真模型结果与试验断裂情况基本一致,从而验证了仿真模型的有效性。
图10. 试验与仿真裂纹面对比图
线弹性断裂力学理论认为,在外载荷的作用下,裂纹尖端会产生弹塑性应力场,而应力强度因子简称K因子,是反映裂纹尖端应力场强弱的唯一指标,应力、应变、位移等均可通过应力强度因子的大小表示。如图11所示,裂纹的基本类型有三种 [
图11. 裂纹分类图
然而在实际的工程问题中,受力状况往往都是很复杂的,不会是上述单一的一种受力情况,而是两种或两种以上基本型的组合,称为复合型裂纹。对于复合型裂纹,我们在FRANC3D中主要采用最大周向力准则来分析计算。
最大周向力断裂准则又叫最大应力准则,对于I-II-III复合型裂纹,其断裂准则为:
k I r ( θ ) = σ 0 2 π r = f I ( k I , k II , k III , θ )
式中:r为径向坐标; θ 为角坐标。
由于铝合金试样一直受轴向的循环载荷作用,边界条件及作用力的方向并未改变,因此在裂纹扩展中的径向坐标和角坐标并未变化,应力强度因子变化的趋势也不会改变。所以取了前四个step的应力强度因子图分析,图的横坐标是采用了量纲为1的裂纹前缘长度来表示应力强度因子所在的位置,即裂纹前缘点到裂纹边界点的弧长与整个裂纹前缘弧长的比值。三幅图中,我们可以明显看出图12(a) I型应力强度因子的数量级最高,这也就是说,I型裂纹是导致试样裂纹扩展的主要因素。并且I型裂纹应力强度因子图类似于开口向上的二次函数曲线,裂纹前缘边角的应力强度因子大于中间部分的应力强度因子,这是因为裂纹前缘呈半圆弧状且与所受横向应力方向相反,所以裂纹越深所受横向应力的影响越大,应力强度因子也就呈现了中间小两边大的情况。在图12(b) II型和图12(c) III型应力强度因子的大小都在0值附近浮动,在断裂力学中,应力强度因子如果小于0,表示对裂纹的扩展起抑制作用,并且相对于I型,II、III型应力强度因子的大小可以忽略不记。也就说II、III型应力强度因子对裂纹扩展不起作用,我们也就不予考虑。
从三幅应力强度因子图中,我们看出试验过程中,样件裂纹主要是I型扩展裂纹,I型扩展裂纹又是工程实际中最易扩展的裂纹。而引起I型裂纹扩展作用力方向垂直于裂纹前缘及裂纹面,正好与试验作用力的方式对应,符合工程实际。
图12. 3种应力强度因子图
通过ABAQUS和FRANC3D对铝合金试样裂纹试验的联合裂纹仿真分析,可以得到以下结论:
1) 通过铝合金试样试验断裂面与仿真预测裂纹断裂面的对比,发现二者裂纹扩展结果基本一致。验证了基于ABAQUS和FRANC3D联合仿真法来模拟疲劳裂纹扩展试验的合理性。
2) 在铝合金疲劳试验中,试样裂纹主要以I型扩展裂纹为主,应力强度因子k随着裂纹的扩展逐渐变大并且I型裂纹应力强度因子图类似于开口向上的二次函数曲线,裂纹前缘边角的应力强度因子大于中间部分的应力强度因子。结合Paris公式计算了疲劳裂纹扩展寿命,得到对应的a-n曲线图:发现随着裂纹扩展的越深,a-n曲线的斜率越大,也就是裂纹扩展的速度越快。
马 俊,崔 颖. 基于ABAQUS和FRANC3D的铝合金裂纹扩展试验研究Experimental Research on Crack Propagation of Aluminum Alloy Based on ABAQUS and FRANC3D[J]. 建模与仿真, 2021, 10(04): 1013-1022. https://doi.org/10.12677/MOS.2021.104101