基于dic的非接触式全场应变测量系统设计毕业设计论文

《基于dic的非接触式全场应变测量系统设计毕业设计论文》由会员分享，可在线阅读，更多相关《基于dic的非接触式全场应变测量系统设计毕业设计论文（95页珍藏版）》请在装配图网上搜索。

1、本科毕业设计（论文）基于dic的非接触式全场应变测量系统设计毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作 者 签 名： 日 期： 指导教师签名： 日期： 使用授权说明本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的

2、印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名： 日 期： 学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名： 日期： 年 月 日学位论文版权使用授权书本学位论文作者完全了解学校有关保留

3、、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期： 年 月 日导师签名： 日期： 年 月 日指导教师评阅书指导教师评价：一、撰写（设计）过程1、学生在论文（设计）过程中的治学态度、工作精神 优 良 中 及格 不及格2、学生掌握专业知识、技能的扎实程度 优 良 中 及格 不及格3、学生综合运用所学知识和专业技能分析和解决问题的能力 优 良 中 及格 不及格4、研究方法的科学性；技

4、术线路的可行性；设计方案的合理性 优 良 中 及格 不及格5、完成毕业论文（设计）期间的出勤情况 优 良 中 及格 不及格二、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？ 优 良 中 及格 不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？ 优 良 中 及格 不及格三、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义 优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？ 优 良 中 及格 不及格3、论文（设计说明书）所体现的整体水平 优 良 中 及格 不及格建议成绩： 优 良 中 及格 不及格（在所选等级前的内画“”）指导教师： （签

5、名） 单位： （盖章）年 月 日评阅教师评阅书评阅教师评价：一、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？ 优 良 中 及格 不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？ 优 良 中 及格 不及格二、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义 优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？ 优 良 中 及格 不及格3、论文（设计说明书）所体现的整体水平 优 良 中 及格 不及格建议成绩： 优 良 中 及格 不及格（在所选等级前的内画“”）评阅教师： （签名） 单位： （盖章）年 月 日 教研室（或答辩小组）及教学系

6、意见教研室（或答辩小组）评价：一、答辩过程1、毕业论文（设计）的基本要点和见解的叙述情况 优 良 中 及格 不及格2、对答辩问题的反应、理解、表达情况 优 良 中 及格 不及格3、学生答辩过程中的精神状态 优 良 中 及格 不及格二、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？ 优 良 中 及格 不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？ 优 良 中 及格 不及格三、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义 优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？ 优 良 中 及格 不及格3、论文（设计说明书）所体现的整体

7、水平 优 良 中 及格 不及格评定成绩： 优 良 中 及格 不及格（在所选等级前的内画“”）教研室主任（或答辩小组组长）： （签名）年 月 日教学系意见：系主任： （签名）年 月 日燕山大学毕业设计（论文）任务书学院： 电气工程学院 系级教学单位：自动化仪表系 学号091203021094学生姓名专 业班 级工业自动化仪表题目题目名称基于DIC的非接触式全场应变测量系统设计题目性质1.理工类：工程设计 （ ）；工程技术实验研究型（ ）；理论研究型（ ）；计算机软件型（ ）；综合型（ ）2.管理类（ ）；3.外语类（ ）；4.艺术类（ ）题目类型1.毕业设计（ ） 2.论文（ ）题目来源科研课题

8、（ ） 生产实际（ ）自选题目（ ） 主要内容1、了解数字图像相关法应变测量原理；2、设计基于DIC的非接触式应变测量系统；3、完成全场应变测量实验，分析测量系统的性能；4、撰写毕业论文基本要求1、 方案设计合理、查阅文献充分；2、 理论分析正确、论证严密；3、 设计、仿真模型可靠，结果与理论相符；4、 毕业设计论文符合撰写规范、要求；参考资料1、牛永强, 物体内部三维位移场测量算法研究. 中国科学技术大学硕士学位论文,2010.2. 冈萨雷斯. 数字图像处理. 北京:电子工业出版社,2011周 次第14周第58周第912周第1316周第1718周应完成的内容查阅资料阅读文献方案论证总体设计理

9、论分析模型研究软件设计仿真调试撰写论文准备答辩指导教师：林洪彬职称：讲师 2012年12月11日系级教学单位审批： 年 月 日摘要摘要近几年来，用数字图像相关法测量变形已成为世人瞩目的测量方法，并且应用越来越广泛，自从2007年以来，国际上基于数字图像相关的方法（Digital Image Correlation,DIC）测量全场应变的技术已经被广泛应用于各种学科的研究，例如材料力学，断裂力学，微观纳米应变测量，各种新型材料测量等。这种测量方法具有非接触性、应用广泛、精度较高、全场测量、数据采集简单、测量环境要求不高、易于实现自动化等优点。伴随着各种高分辨率显微镜的出现，这种测量技术可以测量微

10、米甚至纳米的变形，这样使其拥有了更加广泛的发展空间。本论文凭借国内外在此领域的研究成果作为基础，针对这种方法的测量精度的问题，开展了仿真与实际实验研究，重点分析各种实验因素和外界的干扰对数字图像相关法的计算精度的影响。仿真研究部分，用数值方法模拟从而生成散斑图样，分析了理想状态下散斑与相关计算对精度的影响。借助仿真分析的结果作为基础，开展实际实验的验证，获的了得取最佳计算精度的实际参考值。针对各种产生误差的原因进行分析，并且提出了消除误差或者抑制误差的方法和措施。关键词数字图像相关法；应变和位移测量；测量精度；力传感 I 燕山大学本科生毕业设计（论文）AbstractIn recent yea

11、rs，with the digital image correlation method to measure the deformation has become a measurement method of world attention ，and more and more extensive application，since 2007，method based on digital image correlation（digital image correlation ，dic）of strain measurement technology has been widely use

12、d in many disciplines，such as mechanics of materials，fracture materials，micro nano strain measurement ，various new materials measurement etc，this method of measurement is non-contact，wide application，high precision，full-filled measurement，data collection is simple，measurement environmental requireme

13、nts is not high， easy to realize automation etc.accompanied by a variety of high resolution microscopy，this measurement technique can be the deformation measurement of micron and nanometer，so it has more development space. The diploma by the domestic and foreign research results in this field as a b

14、asis for this method, the measurement accuracy of the problem, research on Simulation and actual experiment, influence the precision of calculating focuses on analysis of various experimental factors and external disturbance on digital image correlation method. The simulation study, simulation to ge

15、nerate the speckle pattern by using numerical method, analyzed the ideal speckle and related effects on the precision of calculation. By the simulation analysis result as the basis, to carry out verification experiments, the optimal calculation precision in practical reference value. To analyze the

16、causes for various errors, and puts forward the methods and measures to eliminate the error or the suppression error.Keywords digital image correlation method; strain and displacement measurement; measuring accuracy; force sensorIII 目 录摘要IAbstractII第1章 绪论11.1 背景和意义11.2 国内外研究现状及趋势21.3 本文的研究内容4第2章 数字图

17、像相关法基本理52.1 数字图像相关法基本原理52.1.1 二维数字图像相关基本原理52.1.2 物体面内变形的表征62.1.3 相关标准72.2数字图像相关搜索方法92.3亚像素位移定位法112.3.1 灰度插值法122.3.2 相关系数拟合法132.4基于位移场的应变估计142.5 本章小结16第3章 基于仿真数字散斑图像的变形测量分析173.1 模拟散斑图的生成原理173.2于仿真散斑图的数字图像相关方法变形测量精度析193.2.1计算窗口大小对计算精度的影响193.2.2 散斑颗粒大小对计算精度的影响213.2.3 散斑颗粒数量对计算精度的影响233.2.4相关函数对计算精度的影响25

18、3.3.本章小结28第4章 数字图像相关变形测量实验研究及误差析294.1 随机散斑生成294.2 单向拉伸试验304.3 实验数据处理314.4 实验影响因素分析324.5 本章小结34结论35参考文献36致谢38附录139附录245附录352附录459附录56539第1章 绪论 第1章 绪论1.1 背景和意义实自上世纪七十年代，光电子技术和数字图像技术的出现，特别是近年来CCD（Charge-coupled Device）摄像机、计算机软硬件及数字图象处理技术的飞速发展，使得数字图像相关(Digital ImageCorrelation ,DIC)法应运而生并得到广泛应用。该方法最初是在上

19、世纪八十年代初由美国南卡罗莱纳州大学（University of South Carolina）的Peter和Ranson提出的。该大学的 Sutton 教授是最早从事该方法的研究人员之一，并一直致力于数字图像相关方法的研究和应用，发表了一系列重要论文，推动了数字图像相关方法的发展，被公认为该研究领域的权威学者。数字图像相关法是一种基于现代数字图像处理和分析技术的新型先进光测技术，它通过分析变形前后物体表面的数字图像来获取被测物体表面的变形信息。数字图像相关法的本身特性决定了它有一些特殊优势：（1）实验设备和实验过程相对简单。被测物体表面的散斑可以是人工喷制，也可以直接以物体表面的自然纹理作为

20、实验标记；另外仅需要单个或两个固定的CCD，来拍摄被测物体表面变形前后的数字图像。通过计算拍摄过程中得到的一组数字图像即可得到物体表面的变形信息。（2）对光源光源和隔震条件要求不高。直接使用白光或者自然光作为照明光源，不需要激光光源和隔振台，这样对测量环境的要求非常低，非常容易实现现场测量。（3）测量过程易于实现自动化。整个实验过程拍摄的都是数字图像，不要传统的胶片记录，省去了繁琐的显影、定影操作，充分发挥了计算机在数字图像处理中的优势。（4）适用测量范围广泛。将数字图像相关方法与不同空间分辨率的图像采集设备结合，来实现对宏观、微观尺度变形的测量。因为数字图像相关方法处理的数据对象是数字化的图

21、像，所以在满足一定的前提条件下，各种高分辨率的数字化图像采集设备所获取的数字图像均可以作为数字图像相关法的处理对象。例如利用各种显微成像设备所获得的数字图像可以可以利用数字图像相关方法来分析处理，从而实现微观变形的测量。同样对高速动态摄影设备所获得的序列数字图像进行处理分析，则可以实现对物体高速瞬时变形的测量。另外，将数字图像相关方法与 X 射线断层扫描成像设备（X-Computed Tomography，X-CT）或同步辐射光 CT 相结合，还可以将数字图像相关法应用与对不透明物体内部三维位移场的测量，这是任何一种其它光学测量方法都难以实现的。由此可见，数字图像相关方法可以很容易的与高空间分

22、辨率和高时间分辨率的图像采集设备相结合，从而实现多种条件下物体变形的测量。将这些先进的图像采集设备与数字图像相关方法结合，不仅直接提高了数字图像相关方法的测量精度，更拓展了其适用范围和应用领域。对数字图像相关方法的研究，不仅具有重要的理论意义，而且更具有非常重要的实际意义。首先，该方法具有许多独特的优势，可以弥补其它测量方法的一些不足；其次，随着高空间分辨率和时间分辨率的图像采集设备的不断出现，该方法可以很方便的应用于新的研究领域。1.2 国内外研究现状及趋势20世纪80年代，来自美国的W.H.Peters、W.F.Ranson和来自日本的I.Yamaguchi同时提出了数字图像相关的算法。W

23、.H.Peters和W.F.Ranson通过对变形前后的灰度图像进行相关迭代可以找到相关系数的极值，从而得到物体的位移值。I.Yamaguchi则用测量物体变形前后光强的相关峰值从而推导出物体的位移。近几十年来，国内外的许多学者对数字图像相关法进行了更加深入的研究，并且取得了不小的成果。由于数字图像相关法更加依赖于计算机的计算速度，所以许多学者就如果提高计算机的计算速度进行了探索研究。Peters在1983年用数字图像相关法进行了刚体测量研究。在同一年，sutton用粗细相关结合的方法对相关搜索进行了进一步的完善和改进，大大提高了计算机的计算速度。随后，T.C.Chu 和Ranson完善了数字

24、图像相关法的测量理论。1989年，Russell用数字图像相关法进行了复合材料的裂纹测量，得出了裂纹二维全场应变。同年，Bruck用插值重建Newton-Raphson迭代算法解决了许多相应的问题，使得数字图像相关法有了很大程度的提高和改进，理论更加趋于完善。1993年，来自加拿大的陆华教授用统计学的方法对数字图像结果进行了误差分析，并且提出了较小误差的办法。1995年，Coburn利用陆华教授的方法对陶瓷损伤进行了实验，实现了陶瓷的无损检测。1996年，Sutton把数字图像相关法应用到了高温测量的领域，同时利用相关的条件和设备把数字图像相关的方法发展到了三维领域。1997年，H.Lu用数字

25、图像相关法对圆柱的变形前后的应变进行了比较和相关的测试。1998年，endroux,Smith改进了精度研究的方法。同年，Y.J.Chao利用此种方法冲击满负荷下裂纹的扩展情况。Rand James L利用数字图像相关方法分析散斑运动，从而得到相关位移，并且通过位移得到应变，做出应变-应力曲线。近几年来，在广大学者的共同努力下，数字图像相关法有了进一步的发展。2001年，Chevalie对橡胶材料的单双轴拉伸进行了数字相关分析，可以得到应力-应变关系。2004年，Yamaguchi对粗糙物体表面也可以进行数字图像相关分析的能力。在国内，数字图像相关方法也得到了广泛的关注，各种数字相关的学术研究

26、也如火如荼的开展起来了。1989年，高建新首次在我国进行了数字相关分析研究，从理论上分析了数字图像相关法，并且进行了物体刚体位移实验测量，但是由于当时实验条件的影响，实验效果不是很理想，精度也不是很高。1992年，李喜德用亚像素的方法对散斑图像进行检测和重建。1993年，李宝琛教授对裂纹金属的损伤进行了全场应变研究，并且得到了比较令人满意的结果。1994年，芮嘉白博士提出了很好的搜索改进方案，例如十字搜索法，在节约时间的基础上也提高了计算的精度。1995年，高建新在生物力学领域提出了多用途数字图像相关的测量系统。1998年，姜锦虎对可能对数字图像相关法造成误差的因素进行了研究。2002年，马少

27、鹏根据岩石材料的变形改进了数字图像相关法。2004年，樊雪松对正交异性材料进行数字相关分析，得到应力应变曲线。相关搜索和相关迭代式数字图像相关法的两种测量方法。相关搜索只能得到位移信息，而应变可以通过位移求导获得。相关迭代可以同时得到位移和应变，是全场测量的方法。在当今时代，随着科技的越来越发达，数字图像相关法慢慢的朝着旧材料到新型材料，普通环境到恶劣环境，静态到动态，微观到宏观的趋势不断发展。1.3 本文的研究内容本文的研究内容可分为三部分。一、对数字图像相关法进行了解并熟悉，通过物体内变形的表征，相关标准，搜索方法和亚像素插值法对数字图像相关法进行研究。二、通过给定位移和变形值，利用数值模

28、拟方法生成仿真散斑图。利用数字图像相关软件计算得出不同散斑颗粒大小和颗粒数量对计算精度的影响，并且通过不同计算窗口的选择确定权衡计算精度和计算速度最佳窗口大小。最后通过模拟不同干扰，比较各相关函数的抗干扰性能。三、利用仿真分析的结果进行实验研究和分析。对试件受力拉伸后变形场进行了测量，同时分析了实验因素对图像采集质量和计算精度的影响。获得了影响计算精度的主要因素，并提出相应的消除误差或抑制误差的方法和有效措施，为后续的实验提高精度提供了可靠的经验。第2章 数字图像相关法基本原理 第2章 数字图像相关法基本理数字图像相关法是依据计算机技术，视觉技术和光电技术的发展而发展起来的一门技术。其原理非常

29、简单，就是通过图像匹配的方法分析试件表面变形前后的散斑图像，对图像进行相关处理得到相应的位移场和应变场。由于该方法对光源和环境的要求都比较低，近年来得到了广泛的应用。并且所需散斑可以是自然纹理，也可以是人工散斑，非常容易实现。本章主要是介绍数字图像相关法的基本原理，以及数字图像相关法的两种搜索方法-搜索法和迭代法，并且介绍了亚像素插值的方法。2.1 数字图像相关法基本原理当光束照射到物体表面时，反射回来的场是一种粒子结构，这个结构式随机的，叫做散斑场。散斑场于物体结构相对应，当物体表面发生变化时，散斑场也跟着发生相应的变化。所以，可以通过匹配散斑图像的方式来获得相应的位移和应变。数字图像相关法

30、的操作过程是先由ccd摄像机进行拍摄变形前后的散斑图，经过A/D转换得到相应的灰度场，然后对两幅灰度图像进行相关运算，找到相应的位移和应变的信息。因为散斑分布是随机的，所以每一点和它周围的散斑是不一样的，我们在相关运算过程中，可以将变形前和变形后的散斑图像分割成很多网格，每一个网格就是一个相应的子集。这样，我们就可以这个子集为载体，分析物体的相应的位移信息，将所有的子集进行计算，就可以得到相应的位移场。在数字图像相关算法中，我们将变形前后的两幅散斑图分别设为F(x ,y)和G(x ,y)。数字图像相关基本思想是在F(x ,y)中找到一个子区，通过子集中的灰度信息，按照一定的搜索方法在变形后的图

31、像G(x ,y)进行相关计算，找到与样本子集相对应的区域，通过分析子集中的位置和形状变化，可以得到物体在该点的位移和应变信息。2.1.1 二维数字图像相关基本原理如图2.1所示，在参考图像的中心点p(x0,y0)周围取（2N+1）（2N+1）像素的子区，用一定的搜索方法在变形后的图像中搜索与参考图像相关系数最大或者是最小的目标子区，该子区必须以p(x0,y0)为中心，从而确定该中心在x和y方向的位移分量分别为u和v，选择正方形为子区是因为正方形比单一的像素点灰度值变换更多，很容易和其他子区进行区分。通常在实际的计算中，我们将参考图像划分为许多虚拟网格的形式，可以通过计算网格的位移来得到物体表面

32、的全场的位移信息。网格节点距离一般设置2-10个像素。YX变形前子区变形后子区vuunvn图2-1 变形前后子区的几何形状2.1.2 物体面内变形的表征为了更好地找到前后图像的匹配程度，我们需要一个相关的标准，匹配的目的就是找到相关系数的最大值，如果这个峰值确定了，那么相应的位置也就确定了。在变形后的图形中子区肯定会发生一定的变化，因此需要对图像进行连续性观察。如图2-1所示，变形后的图像中p点附近的q在变形后的p点附近。因此，参考图像中的各个点可以与变形后的图像子区一一对应，对应所用的函数为： (2-1)式中，a0(Xi,YI),b0(Xi,YI)分别是x和y方向的函数。如果在参考子区和变形

33、子区只存在刚体平移，即子集中各点的位移一样，则可以用零阶形函数来描述： (2-2)显然零阶形函数只能满足简单的子区变形。假设变形前后的目标子区相对于变形前的参考图像子区不但其中心位置发生移动，而且其形状也产生改变，就需要引入形函数来描述这种对应关系。一阶形函数可以描述子区的平移、旋转、剪切、法向应变以及它们的组合。(2-3)在一阶的基础上，lu提出了更加复杂的二阶函数： (2-4)在公式（2-3）-（2-4）中，x=xi-x0, y=yi-y0,u,v是参考子区中点p(x0 y0)。在x和y方向的位移分量是ux uy vx vy参考图像子区的一阶位移梯度,uxx uyy uxy vxx vyy

34、 vxy是参考子区的二阶位移梯度。2.1.3 相关标准如前面所述，参考子区和变形后子区之间的相似度需要一个预先定义的相关标准来进行分析相关标准来进行分析。以P(x0 ,y0)为中心取子区F，假设为(2 N + 1)(2 N+ 1)个像素。则F就记录了P(x0 ,y0)点周围随机分布的斑点灰度值的信息，由统计学定义为二维样本空间。物体移动后，原来子区F处的斑点，就位于子区G处相应位置，斑点间一一对应，这是另一样本空间，概率与统计认为两样本空间完全相关，相关系数为1。若有变形，相关系数会相应下降，但通过位移导数项对子区G尺寸上的修正，相关系数仍然可以达到1。若取一个非相应位置处的子集G ，则G 与

35、F就不完全相关，相关系数就小于1。因此可由两子集的相关性来确定某点P的位移及其导数。相关函数应该满足的以下要求：(1)可操作性(2)可靠性(3)抗干扰性(4)较小的计算量。相关公式是相关运算的关键，先前的文献中提出的与数字图像相关的最常用的相关公式有十多种，这些相关标准大致可分为两大类，互相关准则和平方和准则，如表2-1和表2-2所示。表2-1常用互相关准则互相关准则 定义互相关（cc） Ccc=归一化互相关（ncc） Cncc=零均值归一化互相关函数（zncc）Czncc=表2-2常用平方和准则平方和准则 定义平方和（ssd） Cssd=归一化平方和（nssd） Cnssd=零均值归一化平方

36、和(zssd) Czssd=表2-1和表2-2中，； （2-5）； （2-6）； （2-7）从上述两表可以看出，互相关准则和平方和准则是相关联的。例如 ZNCC 准则可由ZNSSD 准则推导而来,CZNSSD(P)=21- CZNSSD(P),具体的推导可参考文献。类似的NCC准则也可以从NSSD推导而来，CNSSD(P)=21- CNCC(P)，同时也可以看出，当目标子区存在一个线性平移时，利用ZNCC和ZNSSD相关准则计算出来的相关值均不变。因此，ZNCC和 ZNSSD的抗噪声干扰性能最好，同时对光强的补偿和线性放大不敏感。类似的，NCC和NSSD对光强的线性放大不敏感但对光强的补偿敏感

37、。CC和SSD对光强的补偿和线性放大均比较敏感。2.2数字图像相关搜索方法在数字图像相关方法的求解中，相关搜索的技巧经历了从简单到复杂，从大计算量到相对小的计算量，从采用经典的数学理论到采用近现代数学理论的发展过程。相关算法早期采用的方法是粗细搜索法，此方法编程简单，易于实现，其不足是搜索速度慢。针对这个问题，人们提出了十字搜索法，该方法在提高收敛速度方面有所改善。对于较多采用的牛顿迭代法，其相对于粗细搜索法而言，在很大程度上提高了搜索效率。相关算法有很多种，其中典型算法主要有：粗细搜索法、双参数法、牛顿迭代法和梯度算法等。1.双参数法用相关函数进行6个变量(迭代计算时, 首先改变u 和v值,

38、 其余4个变量保持不变。当迭代找到相关系数C的极值时，下一步变化另两个参数和，其余4个变量保持不变，迭代求得对应相关系数极值的和。再接下来变化最后两个参数和，计算对应极值相关系数时的参数值。完成一轮迭代后，再缩小参数的变动范围，重复上述迭代过程，直到满足迭代的终止条件。此时得到的变量（）便为所求的变形参量。2. 粗细搜索法每次用1个像素的步长来搜索计算区域的互相关系数的极值，便得到整像素位移值。如果将搜索的步长改为0.1或0.01像素便可分别获得0.1像素或 0.01像素精度的亚像素位移值。实际使用中通常首先对变形前后的数字图像进行逐点相关计算，找出使相关系数C取极值的位置，这是对真实位移的粗

39、搜索，所得位移值即为整像素位移值。然后是细搜索，在粗搜索的基础上，对子区分别进行亚像素重建，在重建的图像上再进行相关搜索，找出使相关系数C取极值的亚像素值。由于亚像素位置的灰度值必须用插值来获得，这就决定了这种搜索方法必然是很耗时的。3.牛顿迭代法在相关计算中，找到真实位移的条件是相关系数C能够达到极值，即将真实位移及其导数问题转化为数值最优问题。牛顿迭代法基于逐渐改善初值修正项的计算来得到相关系数的极值，在解决最优问题上是比较好的方法。设需要分析的6个参量为： （2-8）式中， t表示转置，牛顿迭代的基本过程为：（1）预估，选定参数的初始值.（2）迭代，迭代所用的公式为 (2-9)式中，i代

40、表迭代次数；为相关系数 在点的 Hessian 矩阵，即二阶导数矩阵(3) 控制，满足时终止迭代，表示所设定的允许误差。(4) 梯度搜索法用f表示变形前的图像，g表示变形后的图像，则有 （2-10）其中(u0,v0)为像素点 (x,y)的位移。为了确定亚像素位移，假设在原图像中以像素点(x,y)为中心的子区。对中的第i个点(xi,yi)的位移可表示为：； （2-11）其中；。代入（2-7）式可得 (2-12)为求解方程中的6个未知量，引入，并设迭代过程：,其中 k 表示迭代次数, 其余各项意义分别为： （2-13） （2-14）； (2-15)； (2-16)这样通过迭代计算，即可得到点处的位

41、移及其梯度值。上述迭代过程的初值，可以采用十字搜索法寻找一个整像素位移点(U , V)，对于所感兴趣区域中的各测量点，分别重复上述迭代过程，即可得到相应的位移场。梯度法在亚像素搜索精度和速度上有明显优势。特别适合于微小变形的测量。2.3亚像素位移定位法由于散斑图记录的是离散的灰度信息，数字图像相关法处理的是数字的图像（最小单位为像素），在相关搜索的时候窗口的平移也只能以像素为单位进行，因此相关搜索所能获得的位移只能是像素的整数倍。然而在实际应用中，位移值一般不会恰好为整像素，而且由于CCD 摄像机的像素有限，整像素位移定位精度在精密测量中远远不够。在数字图像相关方法中，通常采用亚像素定位技术提

42、高测量精度。采用亚像素插值法能够有效地对散斑图像进行相关运算。2.3.1 灰度插值法亚像素位移求解方法主要分为基于灰度插值（或拟合）的亚像素定位法和基于相关系数插值（或拟合）的亚像素定位法。前者一般采用粗细结合的方法，先以像素为最小单位进行位移搜索，得到的位移值是整像素值；再对离散的灰度进行插值或拟合，然后以0.1像素或0.01像素为最小单位重复上述过程，得到亚像素位移部分。插值方法主要有双线性插值和双三次插值等。而后者主要是假设相关函数主峰的分布符合某种模型，对整像素相关搜索结果及其周围相邻8个点组成的相关系数矩阵进行拟合或插值，得到一个连续曲面，然后求该曲面的极值点作为亚像素位移的求解结果

43、。常用的插值方法有高斯曲面插值、抛物面插值和梯度插值等，常用的曲面拟合方法有高斯拟合和二次多项式拟合等。1.双线性灰度插值法一种常用的最简单的插值方法是双线性插值法，其插值函数为： （2-17）对于四条边皆平行于坐标轴的任何矩形平面，有唯一的用于矩形顶点插值的二元线性多项式。假定我们要在一个矩形网格的四个顶点中间一点( x , y )进行插值，并设点(x ,y)由 四 条 边 都 平 行 于 坐 标 轴 的 矩 形 包 围 。此 矩 形 的 顶 点 坐 标 是(x0 ,y0)、(x0 , y1)、( x1 ,y0)、( x1 , y1)，二元线性插值的系数就是由矩形的4个顶点来确定，其函数值分

44、别为f(x0 ,y0)、f(x0 , y1)、f( x1 ,y0)、f( x1 , y1)，把每一个顶点的坐标代入方程，则系数满足如下方程: （2-18） 联立上述四个方程求解得： （2-19） 如果要得到更好的效果，可以采用更高阶的插值算法。双三次样条插值有着优良的性质。从数学上看，双三次样条函数近似于分段的二元六次多项式，在节点处具有一阶和二阶连续导数，光滑程度较高，保证了插值函数的二阶导数连续性。这样双三次样条函数所得到的灰度导数也将有很好的连续性。2. 双三次样条灰度插值法和双线性插值方法类似，通过对 x 方向和 y 方向分别进行一次三次样条插值就可以获得双三次样条插值。用双三次样条插

45、值方式对点(x ,y)进行插值，可以表示为 （2-20） 2.3.2 相关系数拟合法由于相关系数矩阵在以最大值为中心的一个单峰区域上通常近似地满足高斯分布，因此可以通过拟合的方法得到该区域的解析曲面函数，取曲面极值点为目标的亚像素位置。一般采用的拟合方法有高斯拟合和二元多项式拟合，对于相关系数曲面比较平缓的情况，通常采用二元多项式拟合。拟合函数为 （2-21）通常取3 3的拟合窗口，则有 (2-22) 利用最小二乘法，可以求得上述方程组中的6个系数。在拟合曲面的极值点处，应满足： (2-23)由此可以求解得位移 (2-24)2.4基于位移场的应变估计根据2.2和2.3节描述的搜索方法可以获得亚

46、像素精度的位移场，然而在许多材料机械性能测试和结构应力分析中，应变场更能反映材料的性能。应变是位移的数值微分，也就是说如果位移的梯度可以通过算法（如 Newton-Raphson 算法、遗传算法等）进行算出的话，那么相应地就可以获得应变场。值得注意的是，由于运用算法计算位移梯度时不可避免地会引入误差，这使得这种数值微分计算应变的方法仅适用于局部应变大于0.01的场合。另外尽管在理论上应变和位移可以用数值微分的方式来描述，但是微分过程可能会放大位移场中所包含的噪声。因此如果直接进行数值微分，计算出来的应变结果显然是不可靠的。例如，位移的计算误差为 0 .02像素，网格划分的步长为5个像素，那么用

47、前向差分的方法计算得到的应变误差为= (0.02 +0.02)/5 = 8000，用中心差分方法得到的应变误差为前向差分法的一半也即4000。这些误差不可避免地会掩盖被测物体的许多真实应变信息。通过先对位移场进行平滑再进行差分的方法可以提高应变计算的精度。基于这一考虑，Sutton等提出了一种基于罚函数有限元方法（FEM）进行位移场平滑，再用微分进行应变计算的理论。这种方法被Shi等应用于计算电子封装中的热变形。近年来姚学锋等对FEM方法进行了改进。另外Wang和Tong分别提出了薄板样条平滑方法以及其他3种平滑算法来消除位移场中的噪声。平滑操作后明显去除了位移场中所包含的噪声，同时也提高了后

48、续应变计算的精度。然而用FEM方法或薄板样条平滑方法进行平滑的过程都是冗长而复杂的。实际中最常用的应变估计方法是Wattrisse和潘兵等提出的逐点局部最小二乘拟合法。为了分析薄板样件拉伸过程中的局部应变情况，Wattrisse提出了一种基于数字图像相关（DIC）方法的从离散和噪声位移场利用局部最小二乘法来计算应变。为了获得位移计算边缘点的应变，需要在图像边缘对位移场进行连续的延伸。潘兵等也提出过类似的理论，让这种方法更适用于图像边缘、孔、裂纹和其他不连续区域。对原始的离散位移数据用逐点最小二乘拟合的方法来求解应变，由于拟合的过程包含了噪声的去除，因此所得到的应变精度较直接差分会显著提高。该方

49、法的基本思想就是用完全二维多项式对离散位移数据的局部子区（含(2M + 1) (2M+ 1)个均匀分布数据点的计算窗口，如图 2-2 所示）进行分片逐点拟合，利用最小二乘法得到拟合多项式的系数，由拟合多项式的系数就很容易得到拟合区域中间点的对应值和各阶导数，作为该点的平滑和差分的结果。接着，局部子区移动到下一个待计算点，计算新的数据子区的中心点平滑和差分的结果。为了更好地滤除离散位移数据中的噪声和得到较好的平滑效果，考虑拟合函数为二维一次多项式。则拟合的位移场为： （2-25）这里 x, y = M ,M为局部位移场中各数据点的局部坐标；a0,b2 为待求的拟合多项 式 系 数 ；u ( x,

50、 y ), v ( x, y )为离散位移数据点。得到拟合多项式的系数a0, b2后，小变形情况下的应变分量即可按下式进行计算：； （2-26）对于边界点、孔洞、裂纹附件区域由于局部位移数据所包含的有效数据点要小于(2 M + 1)(2 M+ 1)个，在计算系数矩阵 X 和位移向量 u ( x, y ), v ( x, y )时可以忽略这些无效的数据点。实际上由于待求系数向量只有三个未知数，因此只要应变计算窗口中的有效数据点大于三个即可。由于通常的应变计算窗口都取在1111 - 2121点之间，因此以上处理方式是合理可靠的。2.5 本章小结本章介绍了已有的二维数字图像相关方法的基本原理及各种相

51、关搜索方法，比较了现有相关系数表达式的特性，考虑了相关算法对计算精度和速度的影响，对比了不同相关搜索方法的优缺点。为本文后续仿真和实际实验中相关系数和搜索方法的选择提供了依据。结果可归纳为以下几点：（1）数字图像相关最常用的相关函数为互相关准则和平方和准则，并且它们是相关联的。ZNCC和ZNSSD的抗噪声干扰性能最好，同时对光强的补偿和线性放大不敏感。类似的NCC和NSSD对光强的线性放大不敏感但对光强的补偿敏感CC和SSD对光强的补偿和线性放大均比较敏感。（2）位移搜索方法中双参数法需要进行多次的迭代。粗细搜索法易于操作，但是计算量大且耗时。牛顿迭代和梯度算法在计算效率上都有较大的提高，但是

52、迭代初值的选择和梯度的计算对计算的效率有至关重要的影响。（3）对于亚像素重构，主要有插值和拟合两种方法。插值点选择越多，计算量越大，精度越高，一般折衷取五点插值。相关系数拟合通常是利用相关系数矩阵进行拟合得到一个连续曲面，对于相关系数曲面比较平缓的情况通常可以采用二元多项式拟合。第3章 基于仿真数字散斑图像的变形测量分析 第3章 基于仿真数字散斑图像的变形测量分析3.1 模拟散斑图的生成原理 本文采用Peng Zhou提出的算法，用若干随机分布高斯光斑强度迭代的方法模拟变形前后的试件表面的光强分布。首先由计算机模拟生成一幅散斑图作为变形前的散斑图，然后对于图像中的每一点，根据已知位移和应变找到

53、变形后其所对应的位置，将该点的灰度值作为该位置处的灰度值。假定以高斯函数模拟散斑颗粒的光强分布，且图像背景光强均匀则变形前后图像灰度分布函数可分别用式(3-1)和式(3-2)表示： （3-1）（3-2）式中，模拟散斑颗粒的数目 散斑颗粒的大小i0高斯光斑的中心光强，通常设为 1变形前后散斑颗粒的位置由均匀随机分布函数产生，也可以根据具体需要调整随机分布函数，而可根据如下关系得到： (3-3)位移分量可以按照如下关系得到： (3-4)式中,刚体位移位移的一阶导数，即梯度项若要引入非均匀应变，则式（3-4）中的位移梯度项需要满足如下关系： (3-5a) （3-5b) (3-5c) (3-5d)模拟

54、散斑图的算法流程如图 3-1 所示。确定散斑参数（散斑颗粒数目，散斑大小，散斑变形参数）确定变形前散斑颗粒的位置生成基准图像确定变形后散斑颗粒数的位置生成目标图像结束开始图3-1散斑图生成算法流程通过对(3-1)和(3-2)式表达的函数进行离散采样，便能够得到一幅计算机仿真散斑图如图 3-2所示，并可以对其任意平移、转动或施加变形。散斑数目和颗粒大小本节不做过多阐述，所以数量和大小都是随机的，在以下会对散斑数量和大小有进一步的说明。散斑的获取为进行仿真实验奠定了必要的基础。 图3-2（a）初始图像 图3-2（b) 变形后的图像 3.2于仿真散斑图的数字图像相关方法变形测量精度析根据误差分析理论

55、，亚像素位移测量算法的计算误差由均值误差和标准差组成。假设平移散斑图所预加的真实位移为d，而由某种亚像素位移测量算法计算得到的第 i 点(i = 1,2, ,N)的位移为di ，那么该算法的均值误差可定义为: (3-6)这里，为N点位移计算结果的算术平均值，标准差可相应定义为: (3-7)3.2.1计算窗口大小对计算精度的影响计算窗口（子区）大小是数字图像相关中的一个关键的参数。在整像素位移搜索时，一方面，计算窗口越小计算量就越小、计算的效率也越高，但是计算出来的(x ,y )相邻各点的相关系数会有较大变化。另一方面，窗口较大可以补偿、降低各种噪声的影响，计算出来的相邻各点相关系数变化趋于平缓，但是相应地会增加计算量。为了研究计算窗口大小对计算精度的影响，用计算机模拟的方法生成一幅散斑图片，模拟散斑数1500，模拟散斑4 pixel，散斑图像大小428162pixels，变形后的图片相对于初始图片仅在水平方向上平移5 pixel。如果选择不同的子区大小来计算图像中心点的位移（采用互相关函数），计算区域内各点亚像素位移均值、均方差和计算总耗时如图3-3所示。（a） （b） 图3-3 子区大小对计算精度的影响从图中数据可以看出，子区大小在1111 - 4141p