数字图像相关方法若干关键问题研究进展

本文报道了作者及其所在课题组近期在数字图像相关(DIC)测量方法上取得的重要研究进展。主要包括:(1)通过对DIC方法中反向组合高斯-牛顿算法的理论误差分析,提出了新的理论误差公式,进一步证明了反向组合高斯-牛顿算法在提高计算速度和测量精度方面的综合优势;(2)采用提出的理论误差公式,发展了数字散斑场的优化及制作方法,保证了测量结果的一致性和正确性;(3)基于相机阵列和图像拼接技术,发展了超分辨率数字图像相关方法,大大提高了DIC测量方法的应变测量分辨率;(4)提出了大视场条件下的三维系统标定方法,实现了三维测量系统的外参实时标定和多相机测量系统中相机位姿的自动矫正;(5)研制了便携式原位三维测量仪和多尺度DIC测量系统,实现了三维实时数字图像相关测量,进一步满足了DIC方法在工业在线检测和医学领域中的应用需求。     

A Method of Detecting the Onset of Localized Necking Based on Surface Geometry Measurements

To more accurately capture the onset of localized necking and obtain necking limit strains, this paper proposes a method of detecting the onset of localized necking in the Marciniak test (i.e. under in-plane deformation). The method is merely based on the measured surface geometry of the test specimen using digital image correlation (DIC) techniques. It was inspired by the observation of a sudden increase of the surface curvature obtained from 2D curvature fits along the direction across the surface of the sheet. This increase of the surface curvature is detected just before the dimples, which form the final localized neck, become obvious in the DIC measurements. The appearance of this signal is explained by a neck expansion theory defined by propagation of the instability along the direction of the neck, which is a physical behavior of materials.     

On the Use of NURBS Functions for Displacement Derivatives Measurement by Digital Image Correlation

In this paper, we propose to investigate the potential improvement of using Non-Uniform Rational B-Spline (NURBS) functions for displacement measurements by digital image correlation (DIC). The aim is at improving the performance of DIC to capture with low uncertainty and low noise levels not only the displacement field but also its derivatives. Indeed, when the displacement field is used to feed constitutive law identification procedures, displacement derivatives are required and thus may be measured with robustness. Two examples illustrate the potential of NURBS for DIC: a compressive test on a wood sample and a bending test on a steel beam. For the latter, beam kinematics are adopted and NURBS are used in order to capture the variation of the curvature (second derivative of the displacement) along the beam axis. For these two examples, an error study based on a decomposition of the error into the correlation error and the interpolation error, is carried out and shows the great potential of NURBS functions for DIC.     

Effects of the curvature radius of a cylindrical shell on the stress intensity factors of surface cracks

The effects of the curvature radius of a cylindrical shell on stress intensity factors are investigated in circumferential (inner and outer) semielliptical surface cracks in a cylindrical shell. What is new in this paper is to have given: (1) The stress intensity factors for surface cracks in a cylindircal shell are determined by photoelastic technique. (2) By a special method photoelastic slices are handled for obtaining a clear caustic curve, and the stress intensity factors for surface cracks in a cylindrical shell are determined by the caustic method. (3) An approximate equation of curvature correction factor F_c is proposed. (4) Effects of the curvature radius R of a cylindrical shell on the stress intensity factors of surface cracks are obtained. The results of this paper are in fair agreement with already existing analytical results. The approximate equation of curvature correction factor F_c can be widely used for engineering purposes.     

Thin film/substrate systems featuring arbitrary film thickness and misfit strain distributions. Part II: Experimental validation of the non-local stress/curvature relations

The classical Stoney formula relating local equibiaxial film stress to local equibiaxial substrate curvature is not well equipped to handle realistic cases where the film misfit strain, the plate system curvature, and the film thickness and resulting film stress vary with in-plane position. In Part I of this work we have extended the Stoney formula to cover arbitrarily non-uniform film thickness for a thin film/substrate system subject to non-uniform, isotropic misfit strains. The film stresses are found to depend non-locally on system curvatures. In Part II we have designed a demanding experiment whose purpose is to validate the new analysis for the case of radially symmetric deformations. To achieve this, a circular film island with sharp edges and a radially variable, but known, thickness is deposited on the wafer center. The plate system’s curvatures and the film stress distribution are independently measured by using white beam and monochromatic X-ray microdiffraction (μXRD) measurements, respectively. The measured stress field (from monochromatic μXRD) is compared to the predictions of various stress/curvature analyses, all of which have the white beam μXRD measurements as input. The results reveal the shortcomings of the “local” Stoney approach and validate the accuracy of the new “non-local” relation, most notably near the film island edges where stress concentrations dominate.     

扁化效应对弯曲段套管抗挤强度的影响对比分析

针对弯曲段套管抗挤强度随井眼曲率变化规律问题，利用材料力学理论将弯曲段套管等价为均布载荷作用下的弯曲梁模型，根据拉梅厚壁筒理论和强度理论得到套管抗挤强度公式．借鉴扁化变形在圆管类结构中的研究思路，将其拓展到套管抗挤强度分析，对比分析是否考虑扁化效应下四种规格套管不同井眼曲率下的抗外挤强度．研究表明，套管轴向应力随井眼曲率基本呈线性变化；套管抗挤强度随井眼曲率的增加非线性下降．井眼曲率为2o/30 m时，是否考虑套管扁化变形得到的5 1/2"×10.54 mm P110套管抗挤强度误差为12.8 %；井眼曲率为12o/30 m时，误差46.9 %．井眼曲率为2o/30 m时，4 1/2"×9.65 mm TP140V套管抗挤强度误差11.9 %；井眼曲率为12o/30 m时，误差44.7 %．可见，考虑扁化变形后套管抗挤强度小于未考虑扁化变形情况；井眼曲率不大时，两种算法的差距较小，随着井眼曲率的增加，套管扁化变形对套管抗挤强度的影响不能忽略．     

数字图像相关方法中匹配及过匹配形函数的误差分析

基于图像子区的数字图像相关方法需采用合适的形函数来近似目标图像子区的真实变形。由于实际测量时目标子区的局部变形往往是未知的,实际采用的不同阶次(零阶、一阶和二阶)的形函数不可避免地产生误匹配(欠匹配和过匹配)问题,从而引入位移测量的系统或随机误差。尽管由欠匹配形函数引起的系统误差已被充分认识,由过匹配形函数引起的位移误差仍缺少理论解释。本文首先推导出采用一阶和二阶形函数的数字图像相关方法的随机误差理论公式,随后采用一系列数值实验验证了理论公式的准确性。结果显示:过匹配形函数不会引入额外的系统误差,但会增加随机误差,且二阶形函数的随机误差是一阶形函数的二倍。考虑到由欠匹配一阶形函数引入的系统误差往往远大于过匹配二阶形函数的随机误差,因此在未能确知变形的情况下,推荐使用二阶形函数。     

Recent Progress in Digital Image Correlation

In this paper, we report the following important progress recently made in the basic theory and practical implementation of digital image correlation (DIC) for deformation measurement. First, we answer a basic but confusing question to the users of DIC: what is a good speckle pattern for DIC? We present a simple, easy-to-compute yet effective global parameter, called mean intensity gradient, for quality assessment of the entire speckle pattern. Second, we provide an overview of various correlation criteria used in DIC for evaluating the similarity of the reference and deformed subsets, and demonstrate the equivalence of three robust and most widely used correlation criteria, i.e., a zero-mean normalized cross-correlation (ZNCC) criterion, a zero-mean normalized sum of squared difference (ZNSSD) criterion and a parametric zero-mean normalized sum of squared difference (PSSD_ab) criterion with two additional unknown parameters, which elegantly unifies these correlation criteria for subset-based pattern matching. Third, we describe an iterative least squares (ILS) algorithm for accurate subpixel motion detection, which is proved to be equivalent to the existing Newton–Raphson algorithm, but the principle and implementation of ILS algorithm is more straightforward and easier. Finally, to overcome the two limitations of existing subset-based DIC technique, we introduce a robust and generally applicable reliability-guided DIC technique, in which the calculation path is guided by the ZNCC coefficients of computed points, to determine the genuine full-field deformation of an object with complex shape.     

Cross Polarization for Improved Digital Image Correlation

Digital image correlation (DIC) is a surface deformation measurement technique for which accuracy and precision are sensitive to image quality. This work presents cross polarization, the use of orthogonal linear polarizers on light source(s) and camera(s), as an effective method for improving optical DIC measurements. The benefits of cross polarization are characterized through quantitative and statistical comparisons from two experiments: rigid body translation of a flat sample and uniaxial tension of a superelastic shape-memory alloy (SMA). In both experiments, cross polarization eliminated saturated pixels that degrade DIC measurements, and increased image contrast, which enabled higher spatial precision by using smaller subsets. Subset sizes are usually optimized for correlation confidence interval (typically with subsets of 21×21 px or larger), but can be decreased to achieve the highest possible spatial precision at the expense of increased correlation confidence intervals. Smaller subset sizes (such as 9×9 px) require better images to maintain correlation within error thresholds. By comparing DIC results from a uniaxial SMA tension test with unpolarized and cross-polarized images, we show that for 9×9 px subsets, the loss of valid DIC data points was reduced almost ten-fold with cross polarization. The only disadvantage we see to cross polarization is the decrease in specimen illumination due to transmission losses through the polarizers, which can easily be accommodated with sufficiently intense light sources. With the installation of relatively inexpensive linear polarizing filters, an optimum optical DIC setup can provide even better DIC measurements by delivering images without saturated pixels and with higher contrast for increased DIC spatial precision.     

The Geometry of Barotropic Flow

In this paper we construct a new noninvariant Riemannian metric on the semidirect product of the diffeomorphism group of a manifold and the space of positive functions on that manifold, which has the property that certain geodesics give the equations of barotropic fluid mechanics. We compute a formula for its curvature, analyze the sign of the curvature, and determine directly the growth of Jacobi fields in a few special cases.     

损伤拉索的等效弹性模量及其参数分析

损伤拉索会出现线形松弛、应力水平降低的情况,必然会影响拉索的等效弹性模量。本文首先引入损伤程度、位置及范围3个参数,用以描述拉索损伤形态的特征,建立损伤拉索索力和线形计算公式,采用数值方法计算了损伤拉索弦向等效弹性模量精确数值,并和经典的等效弹性模量公式的计算结果进行了比较分析,分析了考虑损伤时两种不同计算方法结果的误差。计算表明,对于500m弦向长度以内的损伤拉索,拉索的弦向长度Lc越大,倾角越小,等效弹性模量的损失越大,并且应用割线模量公式计算的误差也越大,当Lc=500m时,损伤拉索相对误差值在2.5%~4.5%之间。弦向应变越小,等效弹性模量损失越大,弦向应变在[0.001,0.004]内,应用割线模量公式计算的相对误差小于3.5%。损伤程度及损伤范围对引用等效弹性模量公式的误差影响较大,倾角对等效弹模公式相对误差的影响也不容忽视。弦向长度、弦向应变、倾角和损伤程度参数都是通过改变拉索的松弛程度进而影响等效弹性模量的数值以及公式的误差。     

Grid Method for Microscale Discontinuous Deformation Measurement

The objective of this paper is to explore both grid method and Digital Image Correlation (DIC) technique for microscale and discontinuous displacement measurements, such as those associated with crack tips. First, the principle of the grid method is revisited. The grid method and DIC technique are then applied to computer generated images to calculate the displacement field around crack tips. Finally, the grid method is applied to actual experimental images of fracture tests which are conducted inside a Scanning Electron Microscope (SEM) chamber. A new technique is developed to generate microscale pattern that is suitable for both grid method and DIC technique. The displacement fields calculated from grid method are compared with those from DIC technique to identify the strengths and weaknesses of each technique for the microscale and discontinuous displacement measurements. It has been determined that grid method can obtain data closer to the discontinuity than DIC; however, DIC produces smoother displacement fields at the far field. Using this new pattern generation technique, both grid method and DIC technique can be applied to the fracture test at the microscale to complement with each other to achieve the best experiment results.     

A synthetic Schlieren method for the measurement of the topography of a liquid interface

An optical method for the measurement of the instantaneous topography of the interface between two transparent fluids, named free-surface synthetic Schlieren (FS-SS), is characterised. This method is based on the analysis of the refracted image of a random dot pattern visualized through the interface. The apparent displacement field between the refracted image and a reference image obtained when the surface is flat is determined using a digital image correlation (DIC) algorithm. A numerical integration of this displacement field, based on a least square inversion of the gradient operator, is used for the reconstruction of the instantaneous surface height, allowing for an excellent spatial resolution with a low computational cost. The main limitation of the method, namely the ray crossing (caustics) due to strong curvature and/or large surface-pattern distance, is discussed. Validation experiments using a transparent solid model with a wavy surface or plane waves at a water–air interface are presented, and some additional time-resolved measurements of circular waves generated by a water drop impact are discussed.     

Drone-Based StereoDIC: System Development,Experimental Validation and Infrastructure Application

Background

Digital Image Correlation (DIC) is widely used for remote and non-destructive structural health evaluation of infrastructure. Current DIC applications are limited to relatively small areas of structures and require the use of stationary stereo vision camera systems that are not easy to transfer and deploy in remote areas.

Objective

The enclosed work describes the development and validation of an Unmanned Aircraft System (UAS, commonly known as drone) with an onboard stereo-vision system capable of acquiring, storing and transmitting images for analysis to obtain full-field, three-dimensional displacement and strain measurements.

Methods

The UAS equipped with a StereoDIC system has been developed and tested in the lab. The drone system, named DroneDIC, autonomously hovers in front of a prestressed railroad tie under pressure and DIC data are collected. A stationary DIC system is used in parallel to collect data for the railroad tie. We compare the data to validate the readings from the DroneDIC system.

Results

We present the analysis of the results obtained by both systems. Our study shows that the results we obtain from the DroneDIC system are similar to the ones gathered from the stationary DIC system.

Conclusions

This work serves as a proof of concept for the successful integration of DIC and drone technologies into the DroneDIC system. DroneDIC combines the high accuracy inspection capabilities of traditional stationary DIC systems with the mobility offered by drone platforms. This is a major step towards autonomous DIC inspection in portions of a structure where access is difficult via conventional methods.

     

Transient pressure effects in the evolution equation for premixed flame fronts

A nonlinear evolution equation for a scalar field G(x, t) is derived, whose level surface G ₀=const. represents the interface of a thin premixed flame propagating in a flow field. The derivation is an extended version of an equation already proposed by Markstein [1]. It was reconsidered by Williams [2] as a basis for theoretical and numerical analysis and takes, in addition to flame curvature and flame stretch time variations of the bulk pressure, heat loss and nonconstant transport coefficients into account. The equation is an extension of earlier analyses where a flame evolution equation was derived for slightly wrinkled flames such that the front can be described by a single-valued function of a normal coordinate. That formulation excluded situations where the mean flame front has an arbitrary shape in space. Here the more general situation is analysed by using a two-length-scale asymptotic analysis. The leading-order solution of this analysis is equivalent to the equation originally derived by Markstein [1]. In addition to nonconstant properties and heat-loss effects, that had already been considered by Clavin and Nicoli [3], the influence of transient changes of the bulk pressure is analysed. All these effects are combined into a unified formulation which will serve as a basis for a new flamelet concept for premixed turbulent combustion.     

Calculating the Thrust Characteristics of Nozzles from Their Outlet Cross-Section Parameters

Within the framework of boundary-layer theory, simple formulas are presented for finding the viscous loss from the gas parameters in the nozzle outlet cross-section. The analysis is performed for ordinary Laval nozzles and spike nozzles. It is found that for nozzles with large expansion ratios the viscous loss is almost independent of the outlet cross-section parameters and is determined only by the parameter values on the nozzle contour. The effect of the longitudinal nozzle curvature on this loss is investigated. It is shown that the viscous losses calculated from the nozzle outlet parameters and by integrating along the nozzle contour with account for the longitudinal curvature fully coincide.     

Towards Criteria Characterizing the Metrological Performance of Full-field Measurement Techniques

Users of full-field measurement methods like Digital Image Correlation (DIC) often aim to perform measurements with the best trade-off between spatial resolution, bias and measurement resolution. Whenever two full-field methods are compared, it is essential that these criteria are taken into consideration. Recently a metrological efficiency indicator for full-field measurements has been proposed and discussed. This indicator combines measurement resolution and spatial resolution. It has been shown to be invariant to the subset size in the case of Local DIC. The goal of this article is to discuss a method, which determines both the spatial and the measurement resolutions for a given bias for two different DIC methods, in order to obtain the metrological efficiency indicator for each of these methods. The benefit of this indicator is that it does not depend on setting parameters such as the subset size, which are chosen by the user. As such, it can be considered as intrinsic to each technique, thus enabling fair comparison. Local DIC and triangular finite element based Global DIC will be the subject of this investigation. With this setting, their respective subset and triangular element sizes will be related to the spatial resolution of both methods for a given acceptable bias. By using the metrological efficiency indicator, the performance of the two methods will be compared and discussed to a new level of detail. Generally speaking, the indicator shows that the metrological performance of both methods is similar, confirming their popularity. However, it will be shown that, depending on the choice of what an acceptable bias is, one of the method may be preferred to another. The results show that for the specific DIC versions used in the study, for cases for which a significant bias is acceptable, Local DIC outperforms Global DIC, while the opposite is true in the case for which the bias requirements are more stringent. Finally, the quadratic versions of both DIC versions are shown to significantly outperform their respective linear versions.

     

DIC Challenge: Developing Images and Guidelines for Evaluating Accuracy and Resolution of 2D Analyses

With the rapid spread in use of Digital Image Correlation (DIC) globally, it is important there be some standard methods of verifying and validating DIC codes. To this end, the DIC Challenge board was formed and is maintained under the auspices of the Society for Experimental Mechanics (SEM) and the international DIC society (iDICs). The goal of the DIC Board and the 2D–DIC Challenge is to supply a set of well-vetted sample images and a set of analysis guidelines for standardized reporting of 2D–DIC results from these sample images, as well as for comparing the inherent accuracy of different approaches and for providing users with a means of assessing their proper implementation. This document will outline the goals of the challenge, describe the image sets that are available, and give a comparison between 12 commercial and academic 2D–DIC codes using two of the challenge image sets.     

考虑范德华力曲率效应的双壁碳纳米管外压屈曲

针对双壁碳纳米管外压屈曲问题,研究了层间范德华力的曲率效应对临界外压的影响。应用弹性双层圆柱壳模型,考虑层间范德华力不仅与层间距有关而且与挠度曲率的变化有关,导出了外压屈曲临界压力解析公式。计算得出在不同半径、不同长细比下,外压屈曲临界压力的数值结果,并与经典壳的结果和忽略范德华力曲率效应的结果做了比较。结果显示,对于小半径的双壁碳纳米管曲率效应对外压屈曲有效明显的影响。     

钢管高强混凝土叠合柱弯矩—曲率关系的合成模型

本文采用现有高强混凝土及钢管混凝土本构关系,运用合成法编制了计算程序,对钢管高强混凝土叠合柱的截面弯矩－曲率关系进行了全过程分析,并讨论了主要参数对上述关系的影响;经大量回归分析,建立了钢管高强混凝土叠合柱截机弯矩－曲率关系的三折线模型,从而为叠合柱结构的抗震分析奠定了基础。