刘海淙１,李湘萍２

(１．厦门航空有限公司飞机维修工程部,厦门３６１００６;

２．中国民航大学中欧航空工程师学院,天津３００３００)

摘　要:建立了基于三维数字图像相关(３DＧDIC)的非接触光学测试系统,对２A１２铝合金弹性

常数进行了测试,在此基础上分析了DIC子区和步长的大小对材料弹性常数测试结果的影响,并

探讨了２A１２铝合金在拉伸载荷作用下的变形演化规律.结果表明:３组２A１２铝合金试样测得的

弹性模量和泊松比分别为６６．３２,６７．５１,６３．０８GPa和０．３３３,０．３３６,０．３４７;改变子区大小将引起图

像噪音的变化,对测试结果的影响较改变步长的影响更为明显,步长则主要作用于DIC算法的运

行时间;２A１２铝合金试样在整个拉伸过程中依次经过弹性变形、均匀塑性流动以及应变局部化三

个阶段后发生断裂.

关键词:２A１２铝合金;弹性常数;数字图像相关;子区;步长;灵敏度分析

中图分类号:O３４８．１;TG１４６．２ 文献标志码:A 文章编号:１００１Ｇ４０１２(２０１８)０８Ｇ０５５７Ｇ０６

DigitalImageCorrelationTestMethodandItsParameterSensitivityAnalysisforElasticConstantsof２A１２Aluminum Alloy

LIUHaicong１,LIXiangping２

(１．AircraftMaintenanceEngineeringDepartment,XiamenAirlines,Xiamen３６１００６,China;

２．SinoＧEuropeanInstituteofAviationEngineering,CivilAviationUniversityofChina,Tianjin３００３００,China)

Abstract:NonＧcontactopticalmeasurementsystem basedonthreeＧdimensionaldigitalimagecorrelation

(３DＧDIC)wasestablishedtotesttheelasticconstantsof２A１２aluminumalloy．TheeffectsofdifferentsizesofDIC

parameters,i．e．subsetsizeandstepsize,ontheelasticconstanttestresultsofthematerialwereinvestigated,and

thedeformationevolutionlawof２A１２aluminumalloyundertheactionoftensileloadwasalsodiscussed．The

resultsshowthattheelasticmodulusandPoissonratioofthreegroupsof２A１２aluminumalloyspecimenswere

６６．３２,６７．５１,６３．０８GPaand０．３３３,０．３３６,０．３４７,respectively．Changingthesubsetsizecausedthechangeof

imagenoise,subsetsizehadmoreinfluenceonthetestresultscomparedtostepsize,andstepsizemainlycontrolled

thecomputationaltimeofDICalgorithm．The２A１２aluminumalloyspecimenfracturedafterthreestagesofelastic

deformation,uniformplasticflowandstrainlocalization．

Keywords:２A１２ aluminum alloy;elasticconstant;digitalimagecorrelation;subsetsize;stepsize;sensitivityanalysis

方向[５Ｇ６].白晓虹[７]运用二维数字图像相关(TwoＧDimensionalDIC,２DＧDIC)技术,以DP５９０双相钢和３０４不锈钢为试验材料,测得各自的塑性应变比r 值以及全场r 值分布;胡斌等人[８]通过在圆柱体橡胶试样单轴压缩试验中采集变形图像,得到小应变范围内橡胶的应力Ｇ应变曲线,计算出材料的弹性模量.ZHU 等人[９]利用三维数字图像相(ThreeDimensionalDIC,３DＧDIC)测试系统并搭载电子万能试验机,研究了低碳钢棒材在拉伸载荷作用下的全场塑性变形演化过程,并通过３DＧDIC 还原的试样三维形貌获得了瞬时净截面积,进而绘制出了材料的真应力Ｇ真应变曲线.

２A１２铝合金具有高强度、高韧性和可进行热处理强化等优良性能,成为了航空工业中应用最广泛的铝合金之一[１０Ｇ１１].目前,对该材料拉伸弹性常数的测试研究多数采用接触式的测试方法,例如采用引伸计和电阻式应变片,获得标距段的平均工程应变或某些离散位置的点应变.传统测试手段具有以下两点不足:①应变片一般只满足单点、单向测量,难以追踪试样在整个拉伸过程中的变形情况,以及材料在非均匀变化出现后的应变局部化现象;②引伸计有易打滑、

易损伤试件以及自重和夹持力会引起试件附加变形等缺点[１２],给材料应变测量带来一定的误差.笔者利用先进的DIC 光学测试系统和分析软件,对２A１２铝合金在单轴拉伸载荷作用下的全场纵向应变和横向应变进行了试验测量,通过拟合相关曲线获得材料的弹性模量及泊松比;在此基础上,围绕DIC算法运行过程所涉及的两个重要参数子区和步长,对其不同尺寸对测试结果的影响开展了灵敏度分析,为其合理取值提供参考;最后利用DIC全场测量的优越性,对试样表面随拉伸过程的变形演化规律进行了探讨.

１　数字图像相关法及其关键参数

图１为利用３DＧDIC技术测量试样表面三维位移的示意图.首先在左相机拍摄的参考图像中确定分析区域(AreaofInterest,AOI).由于单一像素点的灰度值不具有唯一性,通常不可能在变形图像中识别出与参考图像相同的单个像素点[１３].因此,为了精确辨识测试点P(x０,y０,z０)的空间位置,软件建立了一个以该点为中心并包含(２M ＋１)×(２M ＋１)像素的正方形子区(Subset),通过相关性计算可追踪右相机参考图像中的对应子区.根据预先标定的相机内外部参数可获得P 点的三维坐标.同样地,在变形后左、右相机拍摄的数字图像中也可追踪出该图像子区的相应位置,从而确定变形后中心点的三维坐标P１(x１,y１,z１).二者坐标之差即为所求的三维位移(u,v,w ),其中u 和v 表示面内位移,w 表示离面位移.最后对位移场采用有限应变方程进行差分计算,求得物体表面相应的应变场,即横向应变εxx 、纵向应变εyy 以及剪切应变εxy

在数字图像匹配过程中,DIC技术采用空域迭代的非线性优化算法,如NewtonＧRaphson算法和效率更高的GaussＧNewton算法,结合高精度亚像素插值格式[１４],求解相关函数中用来描述变形子区精确位置和形貌特征的位移向量及位移梯度.DIC技术是全场全过程追踪试样表面变形信息的光学测试方法.

当采集数量较多的变形图像时,如疲劳试验,图像匹配算法将重复运行上百万次,这就对DIC的运行效率提出了相当高的要求.为此,通过设定步长(Step)大小,人为控制子区相关性计算所间隔的像素点数量,可以在一定程度上缩短算法运行时间.由此可知,子区和步长的取值将直接影响DIC测试精度.

２　试样制备与试验方法

２．１　试样制备

T２２８．１－２０１０«金属材料拉伸试验第１部分:室温试验方法»确定如下:厚度为１ mm,总长度为１６０mm,平行段长度为８０mm,宽度为１２．５mm,形状为标准哑铃型.采用黑白哑光漆人工喷制黑白相间的随机散斑,作为DIC追踪试样表面变形的信息载体.

２．２　试验方法

试验装置包括Instron５９８２电子万能试验机,最大载荷为１００kN,同时搭载２台像素分辨率为３２４８×４８７２pixel２的AVTProsilicaGE４９００型相机,采用焦距为１０５mm 的RodenstockRodagon镜头,共同构成３DＧDIC 光学变形测量系统.表１为相机通过圆点标靶标定后获得的内、外部参数.其中:相机光心位置和焦距长度的单位均为像素;角度α,β,γ 分别表示光轴的转角、相机夹角以及相机２相对相机１的水平倾角;Tx ,Ty ,Tz 是２台相机的

间隔距离,分别为－２２０．３５１,０．３１３,４９．９３４mm,基线距离为２２５．９３６mm.由表１可知,２台相机基本被固定在同一水平面,相距２２５．９３６mm,夹角约为２０°.

试样下端固定,上端沿机架竖直向上缓慢移动,并处在相机的视场中央,如图２ 所示.试验开始时,２台相机同时开启,拍摄试样整个拉伸过程的数字图像,采集频率为０．５ 帧??s－１.电子万能试验机按位移方式加载,速率为１mm??min－１,数据采样频率为１点??s－１,分别记录时间、位移、载荷和拉伸应力４ 种信息.在分析软件中,首先划分DIC 计算的分析区域(图２),几何尺寸为１２×７０ mm２,对应分辨率为３１０ × １ ８１０ pixel２, 尺度因子λ ＝２５．８５pixel??mm－１.用于相关性计算的子区和步长初始大小设置为３５×３５pixel２ 和７pixel.

图２　２A１２铝合的３DＧDIC位移测量示意图

Fig.２　Schematicdiagramof３DＧDICdisplacementmeasurement

of２A１２aluminumalloy

３　试验结果与讨论

３．１　２A１２铝合金弹性常数测试

每个采样时刻下AOI纵向工程应变均值εDIC yy 和横向工程应变均值εxDxIC分别为

式中:N 表示DIC测试点数量,(εyy )i 和(εxx )i 分别表示第i 测试点的纵向工程应变和横向工程应变.

图３　２A１２铝合金试样在拉伸应力下的工程应力Ｇ应变曲线

Fig.３　EngineeringstressＧstraincurvesof２A１２aluminumalloy

specimensundertensilestress

工程应力取自电子万能试验机,为瞬时载荷除以初始截面积.由此绘制试样在拉伸载荷作用下的工程应力Ｇ应变曲线,如图３所示.对材料弹性阶段进行线性拟合,相关系数用R２ 表示,直线斜率即为材料的弹性模量E.由图３可见,２A１２铝合金的工程应力Ｇ应变曲线主要分为两个部分,即应力迅速增加的弹性阶段和大范围塑性变形阶段,并且材料没有发生明显的屈服现象.试样的泊松比υ 为拟合横向应变Ｇ纵向应变曲

线所在直线的斜率,即

图４为２A１２铝合金试样在拉伸应力下的横向

应变Ｇ纵向应变曲线.最终利用DIC技术测得３组２A１２铝合金试样的弹性模量、泊松比、抗拉强度及断裂应变,见表２.由表２可知,３号试样的力学性能参数与另外两组试样的差别较大,其弹性模量较１号试样和２ 号试样分别减小了３．２４ GPa 和４．４３GPa.可能是组成３号试样的合金成分含量与另外两个试样的不同,减弱了晶体结构的原子结合

力.１号试样和２号试样的平均弹性模量(６６．９２±０．８４)GPa、平均泊松比０．３３５±０．００２均与文献[１５]给出的结果非常接近,表明该试验环境下DIC方法测得的应变结果是有效的.