角点相关方法在材料力学实验教学中的应用.pdf
《角点相关方法在材料力学实验教学中的应用.pdf》由会员分享,可在线阅读,更多相关《角点相关方法在材料力学实验教学中的应用.pdf(6页珍藏版)》请在咨信网上搜索。
1、2 0 2 3年1 0月第3 9卷第5期皖西学院学报J o u r n a l o f W e s t A n h u i U n i v e r s i t yO c t.,2 0 2 3V o l.3 9 N o.5角点相关方法在材料力学实验教学中的应用周新贵1,2,郝英奇1,叶中豹1(1.安徽建筑大学 土木工程学院,安徽 合肥 2 3 0 6 0 1;2.中国科学技术大学 近代力学系,安徽 合肥 2 3 0 6 0 1)摘 要:材料力学实验是一门重要的工科专业基础实验课,对培养学生的实践创新能力具有重要的作用。基于自主改进的数字图像相关新算法-角点相关方法,建立了材料力学拉伸实验测量系统
2、,以学生为主体进行了铝合金的单轴拉伸实验,对比了角点相关方法和电测法的测量结果,结果表明,两种方法测试结果误差很小,且学生获得了丰富的全场位移信息,拓展了学生的知识面。实验表明,角点相关方法应用于拉伸实验在满足教学要求的前提下提高了测量精度。通过实验不仅调动了学生的学习积极性,也巩固了学生对力学知识的理解,激发了他们的科研兴趣,取得了很好的教学效果。关键词:角点相关方法;实验教学;拉伸实验;全场测量中图分类号:T B 3 0 1-4 2 文献标识码:A 文章编号:1 0 0 9-9 7 3 5(2 0 2 3)0 5-0 1 3 0-0 6*收稿日期:2 0 2 3-0 1-0 5基金项目:安
3、徽省教育厅质量工程项目(2 0 2 1 j y x m 0 2 7 9)。作者简介:周新贵(1 9 8 4-),女,江西鹰潭人,博士,讲师,研究方向:力学。材料力学实验是材料力学教学的一个重要实践性环节,也是我国高校工科学生重要的专业基础课之一1,2。实验课与理论课相互支撑、相互印证,在整个教学活动中,理论课为实验课提供实验原理依据,而实验现象、实验结果能验证理论知识点的正确性,培养学生“实践出真知”的学习研究态度。实验作为课程的重要组成部分,有助于学生对理论知识的理解和知识体系的构建,特别是对培养学生“理论与实践相结合”的能力具有重要作用。材料力学实验是在学习材料力学理论的前提下进行的巩固操
4、作,可以很好地培养学生的实验操作能力和工程协调意识,同时对学生创新能力的培养也有着非常重要的作用。现阶段材料力学实验以室内实验为主,内容可分为测定性实验和验证性实验,主要是为了呈现理论课程教学中的力学现象,或验证力学理论,起到巩固学生力学理论知识、启发学生独立思考的作用。实验项目一般包括金属的拉伸与压缩、扭转、剪切、拉压弹性模量测定、剪切弹性模量测定、纯弯曲梁正应力测定、梁的主应力测定。其中,拉伸实验是普通院校开设率最高的项目,传统的实验教学模式,该实验通常利用万能试验机和引伸计获得材料的应力-应变曲线,学生由此获得不同材料的强度极限、屈服极限、伸长率和断面收缩率等。目前,材料力学实验中材料的
5、变形主要采用电测法,即试件上粘贴的电阻应变片测出的表面应变来计算,该方法仅能获得试件单点的变形,无法观测到试件全场变形。随着材料力学有关问题的开展,如变形局部化以及更复杂的构件破裂过程等所需要的材料全表面变形过程,人们把目光转向了光测法3,4。光测法具有很多优势,并且在各个领域当中起到了非常重要的作用。与材料力学实验中常用的电测法相比,光学测量在观测形变方面有许多优点5,6:除了非接触式测量、自动化程度高、操作简单、测量精度高这些优点,能够获得全场信息更是一大优点。因此,本文基于自主改进的数字图像相关新算法-角点相关算法搭建了新的测量系统,把该系统应用到材料力学实验的材料拉伸实验中,让学生能观
6、察到拉伸过程中试件位移场的变化情况,使学生对金属材料轴向拉伸受力这一知识点得到巩固,同时可以激发学生的科研兴趣以及开阔学生的视野,增添实验课程学习的趣味性,以促进学生能力的培养。1 角点相关方法(C C M)角点相关方法(C C M)为自主改进的数字图像新算法,该方法结合角点提取和数字图像相关法,先通过H a r r i s算法检测出角点的整像素位移,再采用牛顿迭代法得到初始整像素角点变形后的亚像素位移,最终得到被测物体表面的位移场,利用潘兵提出的基于位移场局部最小二乘拟合计算全场应变7。经过模拟及实验验证,该方法得到的位移精度控制在0.0 1像素,应变的分辨率控制在1 0-4量级,且在计算表
7、面发生复杂破坏的模型同样适用。关于角点相关方法的详细介绍及精度评估见笔者文章8。1.1 H a r r i s角点检测原理角点是图像的主要局部特征,这些点类似于编码,但在保持着图形的主要特征的同时,可以用少量的数据代表更多的信息,从而使计算量得到减少,算法的运行速度得到提高。根据应用范围或者算法理论的不同,我们常用角点提取算法,主要可分成以下三类:第一种算法是基于边缘轮廓的提取算法,该算法对图像分割和边缘检测有重要的要求,具有算法复杂,步骤繁琐等缺点;第二种算法是基于模型匹配的角点检测算法,该算法的原理是确定一个给定的模板与图像的所有区域的相似和相关性,这种方法的各类角点需要大量模板来覆盖,在
8、实际中这显然很难做到,此外,多个模板会导致复杂性很大的匹配计算,所以在实际操作中,模板匹配使用的很少;三是基于图像灰度变化的角点检测,该算法是根据图像局部灰度的梯度变化来计算曲率,不需要事先对图像进行分割和边缘检测,其检测角点是通过考察目标像素在一定邻域范围内的灰度变化来完成的,H a r r i s算法就是这种类型的经典算法。1 9 8 8年,H a r r i s算法被提出,用于图像匹配9。角点的检测通常通过图像在不同方向上的灰度改变来确定,所以,假设图像上有一个可向不同方向移动的小窗口,如果经过的小窗口区域都没有灰度改变,则该区域没有角点。当窗口朝一个方向移动,灰度发生了很大变化,而在另
9、一个方向没有灰度改变,则可能是该区域的一条直线。当该区域内像素的灰度值在不同方向上发生了大幅改变时,将该位置的像素定义为角点。角点附近区域的自相关函数如下:Zu,v =x,yw(x,y)Ix+u,y+v -I(x,y)2(1)式中u,v分别为在x,y方向上的平移量,I(x+u,y+v)是(x,y)点平移后的图像灰度值,I(x,y)是(x,y)点在平移前的图像灰度值,w(x,y)是一个窗口函数,其可以是常数函数,也可以是高斯函数。自相关函数Zu,v 的泰勒公式为T(u,v)u,v x,yw(x,y)I2xIxIyIxIyI2yuv (2)其中,Ix、Iy分别为像素在x,y方向上的梯度。令1和2为
10、矩阵特征值,则有:R=12-k(1+2)2(3)k为常数项,其取值范围为0.0 4 0.0 6,角点处特征值1,2都较大。若某一像素的R值大于某一阈值,则该区域存在角点。H a r r i s角点检测算法具有易于实现、计算简单等优点,因此一直以来被广泛应用在各个领域当中1 0。但是,H a r r i s算法也存在如下缺点:一是由于所提取的角点精度都不高,是像素级别,因此角点位置会出现偏离现象和存在伪角点;二是不具有尺度不变性;三是工作效率低,运算量大,因为模板卷积运算过多1 1。针对H a r r i s算法存在的缺陷,国内外学者作了不少改进1 2,1 3。本文在整像素提取采用的是H a r
11、 r i s角点检测。在使用H a r r i s算法提取后,会有小部分点存在角点偏移,如图1所示,约有一个像素的坐标差距,同时整像素的精度也不能满足位移和应变计算的需求,因此要在整像素提取的基础上进行亚像素提取。图1 整像素角点提取示意图1.2 牛顿迭代法牛顿迭代法是数字图像相关方法中最常用亚像素搜索方法,它是通过不断迭代求出互相关函数的极值点,就可得出被测点的亚像素位移。图像子区位移模式选用一阶形函数:x=x0+u+uxx+uyyy=y0+v+vxx+vyy (4)该模式允许变形前后子区发生刚体位移和伸缩变形。其中u,v的初值为角点的整像素位移。选用131第5期周新贵,等:角点相关方法在材
12、料力学实验教学中的应用归一化最小平方距离相关函数:Cf,gp =Mx=-MMy=-M f x,y -fmMx=-MMy=-Mf x,y -fm 2-gx,y -gmMx=-MMy=-Mgx,y -gm 2 2(5)如果子区在变形前后相似度达到最高,则c值为最小,即梯度接近0:?f,gp =Cpi i=1,6=0(6)上式通过牛顿迭代法求解:?p =?p0 +?p0 p-p0 =0(7)整理得:p=p0-?p0?p0 (8)角点相关方法在亚像素提取阶段采用的是牛顿迭代法。1.3 全场应变测量本文采用的方法是由潘兵提出的基于位移场局部最小二乘拟合的全场应变测量方法。该方法的基本思路如下:运用完全二
13、维多项式,分片逐点拟合离散位移数据的局部区域,该离散位移数据以被计算点为中心,再用最小二乘法计算拟合多项式的系数。为了去除噪声,u,v位移场由二维一次多项式对位移场局部子域中的离散数据拟合而来,有:ux,y =a0+a1x+a2yv(x,y)=b0+b1x+b2y(9)用最小二乘法来求解拟合多项式系数a0,b2,u(x,y)、v(x,y)为离散的位移数据点。在N e w t o n-R a p h s o n算法中,假设(x,y)为参考图像子区中任一点,经变形后其对应于目标图像子区中的(x,y)点,则由形函数将两者联系起来:x=x+u+uxx+uyyy=y+v+vxx+vyy(1 0)子区中心
14、为u、v,x0、y0分别为在x和y方向上的位移,ux,uy,vx,vy为图像子区的位移梯度,点(x,y)到点(x0,y0)的距离为x、y,点(x0,y0)为参考图像子区中心点。要求到公式(1 0)中的6个参量,利用N e w t o n-R a p h s o n迭代公式对相关系数公式极值的求解。利用此方法既可以获取子区中心点(x0,y0)两个方向上位移,也可以得到相关的应变信息。当发生小变形时,C a u c h y应变分量可按式(1 1)计算得到:x=ux=a1y=vy=b2(1 1)x y=uy+vx=a2+b1 同理,当发生有限变形时,G r e e n应变分量可按式(1 2)计算得到
15、:Ex x=122ux+ux 2+vx 2 =12(2a1+a21+b21)Ey y=122vy+uy 2+vy 2 =12(2b2+a22+b22)(1 2)Ex y=12uy+vx+uxuy+vxvy =12(a2+b1+a1a2+b1b2)1.4 尺寸标定由于角点相关方法计算的位移信息是以像素为单位,为了让同学们更清楚地了解加载过程中试样表面的变形,需要将像素位移换算成以标准长度为单位计算的实际物理位移,因此,需要对所采集到的图像进行标定。在标定时,需要确定样品的实际物理尺寸与图像像素大小之间的比例关系,即放大倍数,单位为P i x e l/mm。本文研究实验所采用的试样上的网格为规则的
16、矩形,因此,角点相关算法计算位移与实际位移之间的比例关系可以通过网格尺寸来标定。假设实际网格中两个角点的距离为L,我们用程序对所研究的图像进行角点提取,可以计算出网格间两个角点的像素距离为L,则比例关系k为:k=L/L(1 3)假设s为使用角点相关算法计算得到试样表面某一点的位移,则该点的实际位移s 为:s=k*s(1 4)2 基于角点相关方法的拉伸实验教学探索在实际教学中,我们首先向学生介绍了拉伸实验231皖西学院学报第3 9卷的实验目的、实验设备以及实验原理,并详细讲解了角点相关方法的测试原理及其计算方法。具体实验操作流程如下。2.1 试样制备我们采用了铝合金标准试件进行拉伸试验,尺寸大小
17、如图2所示。单位:mm图2 试件尺寸目标区域网格点设置是角点相关法实现的基础,网格点的形状特征也直接影响着角点相关方法的测量精度和运算速度。目前设置网格点的方法有很多,实验中我们采取喷漆镀膜制点的方式,在试件表面待测区域设置特征网格点阵列。在设置网格点之前,首先需要对试件表面进行基本的处理工作。为了得到平整光滑的待测区域,教师指导学生先用锉刀、砂纸对试件待测区域表面进行打磨,去除表面毛刺划痕,避免标记点的凹凸不平。之后用酒精对试件表面杂物进行清洗,再用清水冲洗干净。我们制作的点阵为黑色背景,白色网格点,所以在试件表面清理完成之后用黑色哑光漆进行背景制作,喷漆过程需要保正目标区域均匀覆盖,厚度适
18、中。待背景黑漆晾干后再进行表面网格点的制作,网格点是通过镂空点阵模板制作而成,为了保证获得边界明显,特征清晰的网格点,镂空点阵模板需要孔洞大小一致,边界光滑,孔洞间距不能太近,镂空点阵模板整体厚度适中。实验中的镂空点阵模板尺寸为2 0 mm的铝质正方形模板,孔尺寸4 mm 4 mm,厚度1.5 mm。制点过程将镂空点阵模板放置于试件表面,尽可能与试件贴紧,以免中间存在缝隙。使用白色哑光漆透过孔洞喷洒制点,喷漆过程不宜过多也不能太少,喷漆过多会导致液态流动,致使网格点之间边界模糊;喷漆过少会造成标记点中颗粒离散,不利于自动网格法的特征提取和计算。待网格点白漆晾干,去除镂空点阵模板,获得待测网格点
19、阵列。待测区域全表面均布1 0行5列5 0个4 mm 4 mm相同尺寸的白色方格。应变片贴在试样背面以减少应变片引线对拍摄的干扰,竖向贴着的应变片用于测量轴向变形。试件如图3所示。实验之前,学生采用游标卡尺对试样的几何尺寸进行测量。图3 试件实物图2.2 实验系统角点相关法教学实验系统需准备:C C D相机一台、计算机一台、三脚架一套、补光光源一个,万能试验机,计算机软件部分包括图像采集系统、计算程序等模块。教师指导学生将试件安装至万能试验机夹持端,三脚架放置于试样前,实验拍摄过程中,为了防止离面位移对测量结果的影响,应尽量保证C C D相机测量角度与被测试件保持垂直,为此我们采用可多自由度调
20、节的三脚架固定相机,通过调节高低、俯仰等操作使拍摄角度满足我们的实验要求,继续调整直至待测区域位于图像中心。相机采图软件基于M F C对话框框架编写,能够实现拉伸过程图片的实时采集显示,也可以通过调节帧数来控制采图频率。调节相机光圈和曝光时间,并调节聚焦环,使补光光源对准试样,调整光源亮度及角度,从而保证试件表面的网格点图案的清晰成像。试验系统如图4所示。图4 实验系统2.3 设置参数与加载对万能试验机电脑控制界面设置以位移控制方式0.0 1 mm/s对试件拉伸,静态应变测试仪和C C D331第5期周新贵,等:角点相关方法在材料力学实验教学中的应用相机采集时间间隔均设置为1 s,对拉伸过程采
21、集的图片进行保存,各采集程序同步启动,以保证后续采集数据的匹配。2.4 角点相关方法变形场计算教师指导学生采用角点相关方法对所拍摄的图片进行计算,主要步骤包括:导入所有待计算图像,选择计算区域进行计算,最终获得位移场和应变场云图。2.5 数据处理分别用角点相关方法计算的应变值和应变仪所测值绘制应变-时间曲线,并对结果进行对比分析。2.6 撰写实验报告根据实验目的、步骤、实验结果和分析撰写实验报告。3 实验与教学效果分析同学们以应变片所测值为“真值”,把角点相关法计算值与应变仪测量值的差值作为误差,以每组数据中两条曲线上各对应的误差绝对值的平均值作为该组试验的误差代表值。图5为拉伸试验两种测量方
22、法得到的试件轴向应变-时间曲线。选取贴应变片网格测量结果进行显示。图5 两种结果对比由图5可以看出,角点相关法计算得到的结果与真实应变趋势一致,两者的误差平均值为1.4 91 0-4,且最大误差仅为5.0 31 0-4。试验初始阶段,两者误差较大,出现偏差的主要原因可能为试验机上下夹头的预夹角过程导致试件出现微小的上下滑动引起了测量误差。试验快结束阶段也出现误差较大情况,可能的原因为试件变形较大,应变片遭到破坏导致。由此表明,角点相关方法可以应用于材料力学金属材料拉伸实验中的应变测试,在教学要求得到满足的情况下也提高了测量精度。基于角点相关方法,同学们得到了试件典型时刻的轴向位移场,由图6可知
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 相关 方法 材料力学 实验教学 中的 应用
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。