Identification of damage fields using kinematic measurementsIdentification de champs d'endommagement à l'aide de mesures cinématiques

It is proposed to determine damage fields on the surface of a material by using only full-field displacement measurements. A finite element approach is developed for which the unknown damage field is assumed to be piece-wise constant. Two examples are discussed in which the displacement field is either obtained by finite element analyses or measured by digital image correlation. To cite this article: D. Claire et al., C. R. Mecanique 330 (2002) 729–734.     

Inverse method to determine elastic constants using a circular disk and moiré interferometry

We propose an inverse method, using a circular disk in diametrical compression, for the simultaneous determination of two elastic constants,E andv, from a single displacement map. Moiré interferometry combined with the phase-shifting technique provides a full-field displacement field. An overdeterministic approach using the least-squares method is implemented to fit the experimentally determined displacements to the theoretical solution. An implementation guideline is provided, considering the effects of accidental rigid-body motions, random noise and imperfect position of the origin. Accuracy and repeatability of the proposed method are verified experimentally.     

Investigation of displacement fields in an abrasive waterjet drilling process: Part 2. Numerical analysis

In this paper, the displacement fields associated with the abrasive waterjet (AWJ) drilling process were simulated using the finite element method. A threedimensional finite element model was established, and justifiable pressure loads were used in the numerical model to simulate the AWJ drilling process. It was assumed that the pressure load in the AWJ could be resolved into three components, such as impact jet pressure, shear and normal pressure. The effect of these three pressure loads and their magnitudes on the surface displacement were investigated as a function of the jet penetration depth through numerical modeling. Using the hybrid experimental-numerical stress analysis approach, the transient state of stress and strain associated with the notch crest of the jet-induced hole at the impingement zone of the target material during AWJ piercing can be modeled numerically. It was found that the shear contributed the most in shaping the displacement contour patterns and that the jet pressure did not play a dominant role in determining theu field displacement. The jet pressure and shear had the most effect on thev field displacement contour pattern. It was demonstrated that the principal stresses at the bottom of the cavity increase as the depth of the hole increases.     

A justification of the von Kármán equations

The method of asymptotic expansions, with the thickness as the parameter, is applied to the nonlinear, three-dimensional, equations for the equilibrium of a special class of elastic plates under suitable loads. It is shown that the leading term of the expansion is the solution of a system of equations equivalent to those of von Kármán. The existence of solutions of this system is established. It is also shown that the displacement and stress corresponding to the leading term of the expansion have the specific form generally assumed in the usual derivations of the von Kármán equations; in particular, the displacement field is of Kirchhoff-Love type. This approach also clarifies the nature of admissible boundary conditions for both the von Kármán equations and the three-dimensional model from which these equations are obtained. A careful discussion of the limitations of this approach is given in the conclusion.     

数字图像相关分析法增量位移场测试技术

利用位移场的连续性,对亚像素位移场的算法进行了一些改进,设计了一套分步计算位移场、应变场的测量计算方法,较好地解决了数字图像相关分析法计算精度和效率.采用增量位移场叠加的方法计算大应变位移场,采用局部平面拟合的方法计算应变场.通过对高分子材料拉伸试验位移场的测量和结果标定,说明该方法具有较强的实用性和计算精度.同时,由于避免了对亚像素点的搜索,大大提高了计算效率.     

High-magnification moiré interferometer for crack tip analysis of steels

A high-magnification moiré interferometer, particularly suitable for near-tip field analysis in cracked materials, is described. It has a submillimeter field of view, a high-resolution image sensor (1.4 million pixels), X-Y-Z translation stage and an optical fiber light delivery system. These features enable the microscope head to observe the crack tip while the specimen is loaded in a standard tensile test machine. Automated fringe pattern analysis, using temporal phase shifting and spatial phase unwrapping, enables thex ory displacement component to be measured and the corresponding in-plane strain component computed. The displacement placement accuracy is better than 40 nm, and the effective strain gage dimension is ∼ 25 μm. Furthermore, the interferometer has a built-in white light microscope that allows the observation of the specimen granular microstructure in exact registration with the displacement field. The interferometer has hence been employed to investigate the near-tip fields of a precracked stainless steel specimen under load. The influence of the grain boundaries on the measured displacement fields was relatively minor. The near-tip strain field shows a significant asymmetrical behavior despite pure mode lloading conditions.     

Template matching for improved accuracy in molecular tagging velocimetry

In 2D molecular tagging velocimetry (MTV), tags are written into a fluid flow with a laser grid and imaged at discrete times. These images are analyzed to calculate Lagrangian displacement vectors, often by direct cross correlation. The cross correlation method is inherited from particle imaging velocimetry, where the correlated images contain a random pattern of particles. A template matching method is presented here which takes advantage of the known geometry of laser written tag grids in MTV to achieve better accuracy. Grid intersections are explicitly located in each image by correlation with a template with several linear and rotational degrees of freedom. The template is a continuous mathematical function, so the correlation may be optimized at arbitrary sub-pixel resolution. The template is smooth at the spatial scale of the image noise, so random error is substantially suppressed. Under typical experimental conditions at low imaging resolution, displacement uncertainty is reduced by a factor of 5 compared to the direct cross correlation method. Due to the rotational degrees of freedom, displacement uncertainty is insensitive to highly deformed grids, thus permitting longer delay times and increasing the relative accuracy and dynamic range of the measurement. In addition, measured rotational displacements yield velocity gradients which improve the fidelity of interpolated velocity maps.     

A finite element implementation of the stress gradient theory

In this contribution, a finite element implementation of the stress gradient theory is proposed. The implementation relies on a reformulation of the governing set of partial differential equations in terms of one primary tensor-valued field variable of third order, the so-called generalised displacement field. Whereas the volumetric part of the generalised displacement field is closely related to the classic displacement field, the deviatoric part can be interpreted in terms of micro-displacements. The associated weak formulation moreover stipulates boundary conditions in terms of the normal projection of the generalised displacement field or of the (complete) stress tensor. A detailed study of representative boundary value problems of stress gradient elasticity shows the applicability of the proposed formulation. In particular, the finite element implementation is validated based on the analytical solutions for a cylindrical bar under tension and torsion derived by means of Bessel functions. In both tension and torsion cases, a smaller is softer size effect is evidenced in striking contrast to the corresponding strain gradient elasticity solutions.

     

Fracture prediction for piezoelectric ceramics based on the electric field saturation concept

The electrical field saturation model is applied to the fracture prediction of piezoelectric materials containing electrically impermeable cracks. This model is analogously similar to the electric displacement saturation model that available in the literature. An electrical field saturation strip near the crack front is introduced in the analytical model. The stress intensity factor K and the energy release rate G are obtained in closed-form. It is found that fracture predictions based on K and G criteria are identical. Fracture predictions based on the electric field saturation model and the electric displacement model are also found to be the same.     

Measuring stress intensity factors with a camera: Integrated digital image correlation (I-DIC)Mesure de facteurs d'intensité des contraintes avec une caméra : corrélation d'images numériques intégrée (CINI)

A novel ‘integrated’ approach coupling image correlation and elastic displacement field identification provides a powerful and accurate tool to evaluate mode I and II stress intensity factors. This technique is applied to silicon carbide subjected to a sandwiched three-point bend test, using digital pictures obtained in optical microscopy where the pixel physical scale is about 2 μm. A crack whose maximum opening is 500 nm can be detected and its geometry identified. The toughness is determined well within a 10%uncertainty. To cite this article: F. Hild, S. Roux, C. R. Mecanique 334 (2006).     

The application of Moire interferometry in the measurement of displacement field and strain field at notch-tip and crack-tip

This paper presents the application of Moire interferometry in measuring the displacement and strain field at notch-tip and crack-tip before and after crack propagation. The experiment is carried out using a three point bending beam with a notch. TheN _x andN _y fringe patterns representing displacement field, and the ΔN _x/Δx and ΔN _γ/Δγ fringe patterns representing the strain field are obtained. The sensitivity of the measured displacement is 0.417 μm per fringe order. The displacement and strain distribution along the sectionx=0 have been worked out according toN _x andN _γ fringe patterns. The project supported by Chinese Academy of Sciences and National Natural Science Foundation of China     

Free discontinuity finite element models in two-dimensions for in-plane crack problems

Two different free discontinuity finite element models for studying crack initiation and propagation in 2D elastic problems are presented. Minimization of an energy functional, composed of bulk and surface terms, is adopted to search for the displacement field and the crack pattern. Adaptive triangulations and embedded or r-adaptive discontinuities are employed. Cracks are allowed to nucleate, propagate, and branch. In order to eliminate rank-deficiency and perform local minimization, a vanishing viscosity regularization of the discrete Euler–Lagrange equations is enforced. Converge properties of the proposed models are examined using arguments of the Γ-convergence theory. Numerical results for an in-plane crack kinking problem illustrate the main operational features of the free discontinuity approach.     

应用数字图像相关方法测量含缺陷试样的全场变形

利用数字图像相关方法测量表面带孔洞、裂纹、缺口等缺陷试样的全场变形是许多实际测量任务中经常遇到的问题。就此问题，本文阐述了一种先对要避免计算的缺陷区域进行标记，在随后进行的相关计算中直接避免这些标记区域的方法。在已计算得到全场位移的情况下，文中提出了基于局部位移场最小二乘拟合的方法来计算区域边界、孔洞、裂纹或缺口附近等区域应变。最后对单侧边带半圆缺口试样的单向疲劳拉伸实验的计算结果充分显示本文方法的有效性和可靠性。     

Full-Field Surface 3D Shape and Displacement Measurements Using an Unfocused Plenoptic Camera

Full-field surface 3D shape and displacement measurements using a single commercial unfocused plenoptic camera (Lytro Illum) are reported in this work. Before measurements, the unfocused plenoptic camera is calibrated with two consecutive steps, including lateral calibration and depth calibration. Each raw image of a checkerboard pattern recorded by Lytro Illum is first extracted to an array of sub-aperture images (SAIs), and the center sub-aperture images (CSAIs) at diverse poses are used for lateral calibration to determine intrinsic and extrinsic parameters. The parallax maps between the CSAI and the remaining SAIs at each pose are then determined for depth parameters estimation using depth calibration. Furthermore, a newly developed physical-based depth distortion model is established to correct the serious distortion of the depth field. To realize shape and deformation measurements, the raw images of a test sample with speckle patterns premade on its surface are captured by Lytro Illum and extracted to arrays of SAIs. The parallax maps between the CSAI and the target SAIs are obtained using subset-based digital image correlation. Based on the pre-computed intrinsic and depth parameters and the disparity map, the full-field surface 3D shape and displacement of a test object are finally determined. The effectiveness and accuracy of the proposed approach are evaluated by a set of experiments involving the shape reconstruction of a cylinder, in-plane and out-of-plane displacement measurements of a flat plate and 3D full-field displacement measurements of a cantilever beam. The preliminary results indicate that the proposed method is expected to become a novel approach for full-field surface 3D shape and displacement measurements.     

Dynamics analysis of slender vehicle influenced by random factors

In this paper, the aeroelastic problems of slender vehicles under the influence of random factors and thrust are studied. An aeroelastic dynamic model of a free-free Euler–Bernoulli beam considering thrust and aerodynamic forces is established based on Hamilton’s principle of nonconservative systems. On this basis, considering the influence of random factors, the elastic modulus and viscous drag are regarded as one-dimensional continuous stationary random fields and discretized. The stochastic finite element method is used to solve the dynamic model, and the results are compared with the Monte Carlo simulation results. Then, the influence of the correlation of the random field on the elastic displacement is further analyzed. The following simulation results are obtained: (1) the stochastic factor analysis model established in this paper can reflect the statistical characteristics of aeroelastic response well; (2) the stronger the correlation of the random field is, the greater the expectation of elastic displacement, but as the correlation increases, the expectation tends to be constant; and (3) it is necessary to choose the discrete length of the random field reasonably, and the discrete length depends on the correlation characteristics of the random field studied.     

Mechanically Regularized FE DIC for Heterogeneous Materials

In situ tensile tests in a scanning electron microscope (SEM) have been conducted on a 8-layer 5-harness satin carbon fibre and epoxy matrix composite to observe the first stages of damage at the scale of fibres and matrix. A speckle pattern based on a suspension of alumina particles was applied onto the surface of the specimen to facilitate the use of digital image correlation (DIC). Local and finite element (FE) DIC are compared on pictures acquired during the tensile tests, with and without a speckle pattern. FE DIC with mechanical regularization was found to be the only approach able to measure displacement fields at a fine enough resolution in both cases. This method, initially created for homogeneous materials, was then adapted to heterogeneous materials. First, a microstructure consistent mesh was created and used for correlation purposes. Second, the difference between the mechanical properties of the constituents is taken into account in the mechanical regularization. Last, the accuracy of the method is analysed. The adaptation presented herein was proved to be able to measure displacement fields in the matrix between fibres with an error of 10 nm (a fifth of a pixel) and to detect the initiation of the first damage mechanisms by means of the mechanical residuals.

     

扩展数字图像相关方法中裂尖位移函数的表征研究

探寻适用于扩展数字图像相关方法(X-DIC)的裂纹尖端位移函数,对于提高该方法在裂尖测试精度方面具有重要意义。本文基于断裂力学裂纹尖端位移场函数,进行主导项影响分析,并探究裂尖位移函数中各项及其组合项对位移场表征的贡献程度及对测试精度的影响。通过对I型、II型及混合型裂纹进行测试,发现各项及其组合项在不同裂纹类型测试中的影响有着较大区别,从而得出在X-DIC测试计算中影响裂尖位移场表征的主导项。最后,应用上述研究成果对某航空超硬铝合金的单边裂纹进行拉伸测试。     

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.     

Asymptotic Justification of the Kirchhoff–Love Assumptions for a Linearly Elastic Clamped Shell

The displacement vector of a linearly elastic shell can be computed by using the two-dimensional Koiter's model, based on the a priori Kirchhoff–Love assumptions. These hypotheses imply that the displacement of any point of the shell is an affine function of the transverse variable x ₃. The term independent of x ₃ of this approximation is equal to the displacement vector of the two-dimensional Koiter's model. The term linear in x ₃ depends on the infinitesimal rotation vector of the normal. After an appropriate scaling, we estimate here the difference between the three-dimensional displacement and this affine vector field in the case of shells clamped along their entire lateral face. Besides, in the case of shells with uniformly elliptic middle surface, taking into account the term depending of the rotation of the normal allows to improve the asymptotic estimate between the three-dimensionnal displacement and Koiter's bidimensional displacement. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the static and dynamic equilibrium of concentrated loads in linear acoustics and linear elasticity

The displacement field in an unbounded linear elastic fluid subjected to a time-dependent point force is obtained by using integral transform techniques. Differentiation of the displacement field yields the pressure field. It is shown that the pressure on the surface of a spherical ball B _r of radius r centered at the point where the force is applied is statically equivalent in the limit as r0 to only one-third of the force. The remaining two-thirds are carried by the inertia terms. It is also shown, by an independent reasoning, that a point force cannot be carried in static equilibrium by a linear elastic fluid.The displacement field corresponding to an unbounded isotropic linear-elastic solid subjected to a time-dependent point force (the Stokes solution) is also obtained by using integral transform techniques. As is well-known, the tractions of the Stokes solution on the surface of a spherical ball B _r are statically equivalent in the limit as r0 to the force itself; consequently, the inertia terms do not contribute to the dynamic equilibrium of B _r. The contrast between the response of a fluid and that of an isotropic solid under the action of a point force is discussed.