CORRELI - LMT Cachan - ENS Cachan
CORRELI - LMT Cachan - ENS Cachan
CORRELI - LMT Cachan - ENS Cachan
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>CORRELI</strong> Q4 :<br />
A Software for “Finite-element”<br />
Displacement Field Measurements<br />
by Digital Image Correlation<br />
François HILD and Stéphane ROUX<br />
April 2008<br />
Internal report no. 269<br />
<strong>LMT</strong>-<strong>Cachan</strong><br />
<strong>ENS</strong> de <strong>Cachan</strong>/CNRS-UMR 8535/Université Paris 6/PRES UniverSud Paris<br />
61 avenue du Président Wilson, F-94235 <strong>Cachan</strong> Cedex, France<br />
Email: {francois.hild,stephane.roux}@lmt.ens-cachan.fr<br />
1
Abstract: Internal report no. 269<br />
This internal report completes the first ones (i.e., no. 230 and 254) on digital image correlation.<br />
The technique developed herein is based upon a multi-scale approach to determine “finite element”<br />
displacement fields by digital image correlation. The displacement field is first estimated<br />
on a coarse resolution image and progressively finer details are introduced in the analysis as the<br />
displacement is more and more securely and accurately determined. Such a scheme has been<br />
developed to increase the robustness, accuracy and reliability of the image matching algorithm.<br />
The details of the program are then presented. The procedure, <strong>CORRELI</strong> Q4 , is implemented<br />
in Matlab TM . The different steps are presented. The procedure is used on one example dealing<br />
with Portevin-Le Châtelier bands in an aluminum alloy, and with an artificially deformed<br />
picture of stone wool. Other published applications using the present code are listed.<br />
The software has been protected under the IDDN.FR.001.110004.000.S.P.2008.000.21000.<br />
Résumé : rapport interne n ◦ 269<br />
Ce rapport interne complète les précédents (i.e., n ◦ 230 et 254) sur la corrélation d’images<br />
numériques. La technique présentée ici est basée sur une détermination multi-échelles par<br />
corrélation d’images d’un champ de déplacement de type “élément fini”. Le champ de déplacement<br />
est d’abord déterminé sur une image de résolution grossière et des détails plus fins sont<br />
rajoutés au fur et à mesure que les évaluations sont obtenues de manière plus robuste et sure.<br />
L’algorithme développé ici a pour but d’augmenter la robustesse, la précision et la fiabilité de la<br />
technique de corrélation. Les détails du programme, <strong>CORRELI</strong> Q4 , implanté dans Matlab TM ,<br />
sont ensuite présentés. Ce programme est enfin utilisé dans un exemple correspondant à une<br />
bande de Portevin-Le Châtelier dans un alliage d’aluminium, et un autre correspondant à<br />
l’analyse d’une image artificiellement déformée d’une laine de roche. D’autres applications<br />
publiées utilisant le code de corrélation présenté ici sont listées.<br />
Le logiciel est protégé sous l’IDDN.FR.001.110004.000.S.P.2008.000.21000.<br />
2
1 Introduction<br />
The analysis of displacement fields from mechanical tests is a key ingredient to bridge the gap<br />
between experiments and simulations. Different optical techniques are used to achieve this<br />
goal [1]. Among them, digital image correlation (DIC) is appealing thanks to its versatility<br />
in terms of scale of observation ranging from nanoscopic to macroscopic observations with<br />
essentially the same type of analyses. Most developments based on correlation exploit mainly<br />
locally constant or linearly varying displacements [2].<br />
In Solid Mechanics, the measurement stage is only the first part of the analysis. The<br />
most important application is the subsequent extraction of mechanical properties, or quantitative<br />
evaluations of constitutive law parameters [3]. By having an identical description for the<br />
displacement field during the measurement stage and for the numerical simulation is the key<br />
for reducing the noise or uncertainty propagation in the identification chain. During the latter,<br />
there is usually a difference between the kinematic hypotheses made during the measurement<br />
and simulation stages. To avoid this source of noise, it is proposed to develop a DIC approach<br />
in which the measured displacement field is consistent with a finite element simulation. Consequently,<br />
the measurement mesh has also a mechanical meaning. Let us emphasize that in<br />
the present study, only the displacement field measurement is considered (i.e., it is a DIC technique),<br />
and no (finite-element) mechanical computation is performed. No constitutive law has<br />
been chosen, nor any identification performed. However the displacement evaluation is directly<br />
matched to a format ready to use for any further finite element modeling work.<br />
In the following, it is proposed to develop a Q4-DIC technique in which the displacements<br />
are assumed to be described by Q4P1-shape functions relevant to finite element simulations [4].<br />
The pattern-matching algorithm is based upon the conservation of the optical flow. Variational<br />
formulations are derived to solve this ill-posed problem. A spatial regularization was introduced<br />
by Horn and Schunck [5] and consists in a looking for smooth displacement solutions. The<br />
quadratic penalization is replaced by “smoother” ones based upon robust statistics [6, 7, 8]. In<br />
the present approach, the sought displacement field directly satisfies continuity. In its direct<br />
application, the conservation of the optical flow is a non-linear problem that is expressed in<br />
terms of the maximization of a correlation product when the sought displacement is piece-wise<br />
constant [9]. Other kinematic hypotheses are possible and a perturbation technique of the<br />
minimization of a quadratic error leads to a linear system as in finite element problems. To<br />
increase the measurable displacement range, a multi-scale setting is used as was proposed for<br />
a standard DIC algorithm [10].<br />
The report is organized as follows. Section 2 presents the general principles of a DIC<br />
approach. It is particularized to Q4P1-shape functions and it is thus referred to as Q4-DIC. In<br />
3
Section 3, all the details are given to run correli_q4.m as a Matlab TM file and discussed for an<br />
artificially deformed picture of stone wool. A picture of an aluminium alloy sample constitutes<br />
another test case discussed in Section 4 for the quantitative analysis of a Q4P1 kinematics as<br />
it offers a good illustration of a heterogeneous strain field (i.e., a localized band is observed<br />
in a tensile test). Section 5 briefly summarizes other applications with the same correlation<br />
technique, and Section 6 with extensions of Q4-DIC. Details can be found in the listed papers.<br />
2 Q4-Digital Image Correlation (Q4-DIC)<br />
In this section the principle of the perturbation approach is introduced. Let us underline that<br />
this approach applies to a wide class of functions, and is not confined to finite element shape<br />
functions. Other examples have been explored [11, 12, 13], using mechanically based functions,<br />
or using spectral decompositions of the displacement field [15, 14]. However, the discussion will<br />
be specialized to Q4P1-shape functions, which provide a versatile tool for the analysis of very<br />
different mechanical problems, ideally suited to finite element modeling.<br />
2.1 Principle of DIC with an arbitrary displacement basis<br />
Let us deal with two images, which characterize the original and deformed surface of a material<br />
subjected to a known loading. An image is a scalar function of the spatial coordinate that gives<br />
the gray level at each discrete point (or pixel) of coordinate x. The images of the reference<br />
and deformed states are respectively called f(x) and g(x). Let us introduce the displacement<br />
field u(x). This field allows one to relate the two images by requiring the conservation of the<br />
optical flow<br />
g(x) = f[x + u(x)] (1)<br />
Assuming that the reference image are differentiable, a Taylor expansion to the first order yields<br />
g(x) = f(x) + u(x).∇f(x) (2)<br />
Let us underline here that the differentiability of the original image is not simply an academic<br />
question, but we will come back to this point later on. The measurement of the displacement is<br />
an ill-posed problem. The displacement is only measurable along the direction of the intensity<br />
gradient. Consequently, additional hypotheses have to be proposed to solve the problem. For<br />
example, if one assumes a locally constant displacement (or velocity), a block matching procedure<br />
is found. It consists in maximizing the cross-correlation function [16, 9]. To estimate<br />
u, the quadratic difference between right and left members of Eq. (2) is integrated over the<br />
4
studied domain Ω and subsequently minimized<br />
∫∫<br />
η 2 = [u(x).∇f(x) + f(x) − g(x)] 2 dx (3)<br />
Ω<br />
The displacement field is decomposed over a set of functions Ψ n (x). Each component of the<br />
displacement field is treated in a similar manner, and thus only scalar functions ψ n (x) are<br />
introduced<br />
u(x) = ∑ a αn ψ n (x)e α (4)<br />
α,n<br />
The objective function is thus expressed as<br />
∫∫ [ 2<br />
∑<br />
η 2 = a αn ψ n (x)∇f(x).e α + f(x) − g(x)]<br />
dx (5)<br />
Ω α,n<br />
and hence its minimization leads to a linear system<br />
∑<br />
∫∫<br />
∫∫<br />
a βm [ψ m (x)ψ n (x)∂ α f(x)∂ β f(x)]dx = [g(x) − f(x)] ψ n (x)∂ α f(x)dx (6)<br />
β,m<br />
Ω<br />
Ω<br />
that is written in a compact form as<br />
Ma = b (7)<br />
where ∂ α f = ∇f.e α denotes the directional derivative. The matrix M and the vector b are<br />
directly read from Eq. (6)<br />
∫∫<br />
M αnβm = [ψ m (x)ψ n (x)∂ α f(x)∂ β f(x)]dx (8)<br />
Ω<br />
and<br />
∫∫<br />
b αn = [g(x) − f(x)] ψ n (x)∂ α f(x)dx (9)<br />
Ω<br />
Let us note that the role played by f and g is symmetric, and up to second order terms,<br />
exchanging those two functions will lead to simply exchanging the sign of the displacement.<br />
Thus in order to compensate for variations of the texture and to cancel the induced first order<br />
error in u, one substitutes f in the expression of the matrix M by the arithmetic average<br />
(f + g)/2. This symmetrization turns out to make the estimate of a more stable and accurate,<br />
although it requires more computation time associated with the computation and assembly of<br />
all elementary matrices and vectors.<br />
Last, the present development is similar to a Rayleigh-Ritz procedure frequently used in<br />
elastic analyses [4]. The only difference corresponds to the fact that the variational formulation<br />
is associated to the (linearized) conservation of the optical flow and not the principal of virtual<br />
work.<br />
5
2.2 Particular case: Q4P1-shape functions<br />
A large variety of functions Ψ may be considered. Among them, finite element shape functions<br />
are particularly attractive because of the interface they provide between the measurement of the<br />
displacement field and a numerical modeling of it based on a constitutive equation. Whatever<br />
the strategy chosen for the identification of the constitutive parameters, choosing an identical<br />
kinematic description suppresses spurious numerical noise at the comparison step. Moreover,<br />
since the image is naturally partitioned into pixels, it is appropriate to choose a square or<br />
rectangular shape for each element. This leads us to the choice of Q4-finite elements as the<br />
simplest basis. Each element is mapped onto the square [0, 1] 2 , where the four basic functions are<br />
(1−x)(1−y), x(1−y), (1−x)y and xy in a local (x, y) frame. The displacement decomposition<br />
(4) is therefore particularized to account for the previous shape functions of a finite element<br />
discretization. Each component of the displacement field is treated in a similar manner, and<br />
thus only scalar shape functions N n (x) are introduced to interpolate the displacement u e (x)<br />
in an element Ω e<br />
∑n e<br />
∑<br />
u e (x) = a e αnN n (x)e α (10)<br />
n=1 α<br />
where n e is the number of nodes (here n e = 4), and a e αn the unknown nodal displacements. The<br />
objective function is recast as<br />
η 2 = ∑ ∫∫ [ 2<br />
∑<br />
a e αnN n (x)∇f(x).e α + f(x) − g(x)]<br />
dx (11)<br />
e Ω e α,n<br />
and hence its minimization leads to a linear system (6) in which the matrix M is obtained<br />
from the assembly of the elementary matrices M e whose components read<br />
∫∫<br />
Mαnβm e = [N m (x)N n (x)∂ α f(x)∂ β f(x)]dx (12)<br />
Ω e<br />
and the vector b corresponds to the assembly of the elementary vectors b e such that<br />
∫∫<br />
b e αn = [g(x) − f(x)] N n (x)∂ α f(x)dx (13)<br />
Thus it is straightforward to compute for each element e the elementary contributions to M<br />
and b. The latter is assembled to form the global “mass” matrix M and “force” vector b,<br />
as in standard finite element problems [4]. The only difference is that the “mass” matrix and<br />
the “force” vector contain picture gradients in addition to the shape functions, and the “force”<br />
vector includes also picture differences. The matrix M is symmetric, positive (when the system<br />
is invertible) and sparse. These properties are exploited to solve the linear system efficiently.<br />
Last, the domain integrals involved in the expression of M e and b e require imperatively a pixel<br />
6
summation. The classical quadrature formulas (e.g., Gauss point) cannot be used because of<br />
the very irregular nature of the image texture. This latter property is crucial to obtain an<br />
accurate displacement evaluation.<br />
2.3 Sub-pixel interpolation<br />
In the previous subsection, the gradient ∇f(x) is used freely in the Taylor expansion leading<br />
to Eq. (2). However, f represents the texture of the initial image, discretized at the pixel<br />
level. Therefore, the definition of a gradient requires a slight digression. Previous works have<br />
underlined the importance of sub-pixel interpolation. In Ref. [17] a cubic spline was argued to<br />
be very convenient and precise. Here a different route is proposed, namely, a Fourier decomposition.<br />
The latter provides a C ∞ function that passes by all known values of f at integer<br />
coordinates. From such a mapping one easily defines an interpolated value of the gray level<br />
at any intermediate point. Moreover, one also exploits the same mapping for computing a<br />
gradient at any point. Finally, powerful Fast Fourier Transform (FFT) algorithms allow for a<br />
very rapid computation.<br />
There is however a weakness in this procedure related to the treatment of edges. Fourier<br />
transforms over a finite interval implicitly assume periodicity. Thus left-right or up-down<br />
differences induce spurious oscillations close to edges. To reduce edge effects, each element<br />
is enlarged to an integer power of two size, including a frame around each element. This<br />
enlarged element is only used for FFT purposes, and once gradients are estimated, the original<br />
element is cut out the enlarged zone, and thus the region where most of spurious oscillations are<br />
concentrated is omitted. Moreover, at present, an “edge-blurring” procedure is implemented,<br />
i.e., each border is replaced by the average of the pixel values of the original and opposite<br />
border ones. This again reduces the discontinuity across boundaries [10]. There exist a few<br />
alternative routes to limit or circumvent part of this artefact, namely, ad hoc windowing [10],<br />
neutral padding [18], symmetrization, or linear trend removal. Such options have not been<br />
tested.<br />
2.4 Multi-scale approach<br />
Even though a way of interpolating between gray level values at a sub-pixel scale was introduced<br />
above, the very use of a Taylor expansion requires that the displacement be small when<br />
compared with the correlation length of the texture. For a fine texture and a large initial<br />
displacement, this requirement appears as inappropriate to converge to a meaningful solution.<br />
Thus one may devise a generalization to arbitrarily expand the correlation length of the tex-<br />
7
ture. This is achieved through a coarse-graining step. Again many ways may be considered,<br />
such as a low pass filtering in Fourier or Wavelet spaces. A rather crude, but efficient way, is<br />
to resort to a simple coarse-graining in real space [10] obtained by forming super-pixels of size<br />
2 n × 2 n pixels, by averaging the gray levels of the pixels contained in each super-pixel.<br />
First one generates a set of coarse-grained pictures of f and g for super-pixels of size<br />
2 × 2 pixels, 4 × 4 pixels, 8 × 8 pixels and 16 × 16 pixels. Starting from the coarser scale, the<br />
displacement is evaluated using the above described procedure. This determination is iterated<br />
using a corrected image g where the previously determined displacement is used to correct for<br />
the image. These iterations are stopped when the total displacement no longer varies. At this<br />
point, one may estimate that a gross determination of the displacement has been obtained,<br />
and that only small displacement amplitudes remain unresolved. This lack of resolution is<br />
due to the fact that the small scale texture was filtered out. Thus finer scale images are used<br />
taking into account the previously estimated displacement to correct for the g image. Again<br />
the displacement evaluation is iterated up to convergence. This process is stopped once the<br />
displacement is stabilized at the finer scale resolution, i.e., dealing with the original images.<br />
Along the iterations, the “correction” of the deformed image by the previously determined<br />
displacement field are possible with different degrees of sophistication. For reasons of computation<br />
efficiency, only the most crude correction is performed in the present implementation,<br />
namely, each element is simply translated by the average displacement in the element. Integer<br />
rounded displacements are taken into account by a mere shift of coordinates, and sub-pixel<br />
translation is performed by a phase shift in Fourier space [19]. This is a very low cost correction<br />
since Fourier transforms are already required to compute gradients.<br />
At the present stage, the implementation is such that the same number of super-pixel is<br />
contained in each element. Thus as a finer resolution image is considered, the displacement is to<br />
be determined on a physically finer grid. The transfer of the displacement from one scale to the<br />
next one is performed using a linear interpolation, consistent with the Q4P1-shape functions<br />
that are used. This multi-resolution scheme is thus also a mesh refinement procedure which is<br />
performed uniformly (up to now) over the entire map.<br />
This multi-resolution scheme was previously implemented using an FFT-correlation approach<br />
to estimate the displacement field [10]. In this context, it leads to much more robust<br />
results. Large displacements and strains are measured using this algorithm, whereas a single<br />
scale procedure revealed to be severely limited. Similarly, using the present Q4-decomposition,<br />
this multi-resolution analysis revealed very precious to significantly increase the robustness and<br />
accuracy of the measurement.<br />
8
3 <strong>CORRELI</strong> Q4 : User’s guide<br />
The following section describes all the steps that can be followed when using the Q4-DIC code.<br />
First start Matlab TM : at least the version 5.3 is needed (or any newer one; versions 7.xx have<br />
been tested). Choose the directory in which the <strong>CORRELI</strong> files are put as the current directory.<br />
Type the command correli_q4 at the MATLAB prompt. The first step is to choose in a first<br />
menu (Fig. 1) to run a priori (performance) analyses, a Q4-DIC computation, to visualize a<br />
result, or to create a movie:<br />
• Texture: texture analysis;<br />
• Uncertainty: uncertainty analysis;<br />
• Resolution: resolution analysis;<br />
• Computation click: choice of the Region Of Interest (ROI) by mouse click;<br />
• Computation restart: get the ROI coordinates of a previous computation;<br />
• Computation data: give the coordinates (in pixels) of the ROI to analyze;<br />
• Visualization: visualize results of any previous computation.<br />
• Movies: generate movies (only active for versions greater than or equal to 7.0).<br />
¢¡¤£¦¥ §¨£¦¥©¨© <br />
¨<br />
§¨!"¡¤#%$&©'$(¥ §¨<br />
)*¥ +,#©¨ ¥ +%©'$(¥ §<br />
Figure 1: Menu to choose the type of analysis.<br />
9
3.1 A priori analyses<br />
The aim of the present section is to evaluate the a priori performances of the Q4-DIC technique<br />
applied to the picture that corresponds to the reference configuration of the experiment to be<br />
analyzed in Section 4. Figure 2 shows the texture used to measure displacement fields. It is<br />
obtained by spraying a white and black paint prior to the experiment. Another picture (of<br />
stone wool) will also be analyzed.<br />
¢ <br />
0.03<br />
<br />
<br />
<br />
0.025<br />
¤¥¡£¡<br />
¦§¡£¡<br />
Frequency<br />
0.02<br />
0.015<br />
0.01<br />
¥¡£¡<br />
0.005<br />
¨©¡£¡<br />
¡¥¡£¡<br />
¢¡£¡ ¤¥¡¥¡ ¦§¡£¡ ¨©¡£¡ ¡¥¡£¡ ¥¡¥¡<br />
0<br />
1 2 3 4 5<br />
Gray level<br />
x 10 4<br />
Figure 2: View of a reference picture and corresponding gray level histogram. The tension axis<br />
is vertical. The width of the sample is 30 mm.<br />
3.1.1 Texture characteristics<br />
The quality of the displacement measurement is primarily based on the quality of the image<br />
texture. Hence before discussing the result of the analysis, the characteristics of the texture<br />
are presented. The gray level was encoded on a 16-bit depth (even though the original depth<br />
was equal to 12 bits) in the image acquisition, and the true gray level dynamic range takes<br />
advantage of this encoding, as judged from the gray level histogram shown in Fig. 2. Such a<br />
histogram is a good indication of the global image quality to check for saturation problems.<br />
However, such a global characterization of the image is only of limited interest. It is mostly<br />
useful at the stage of image acquisition to set, say, the exposure time and / or the aperture.<br />
Many acquisition softwares offer such functions. However, since the actual dynamic range of<br />
gray levels is an important element to appreciate the quality of a picture, it is included here as<br />
a possible diagnostic tool of poor performance.<br />
What is more significant is the average of texture properties as estimated from sampling of<br />
sub-images in elements. This is more a characteristic of the patterns than of image acquisition.<br />
The point is to evaluate whether the sub-images carry enough information to allow for a proper<br />
10
analysis. Each element is characterized by its own gray level dynamic range, or its standard<br />
deviation of gray level. The latter quantity, averaged over all elements of a given size, and<br />
normalized by the maximum gray level used in the image, is shown in Fig. 3-a. Even for the<br />
smallest element sizes, this ratio is already as large as 0.06, and increases to about 0.13 for large<br />
element sizes. The higher the ratio, the smaller the detection threshold as shown in Eq. (2);<br />
the standard deviation being an indirect way of characterizing the sensitivity of the technique.<br />
One thus concludes that the gray level amplitude is large enough to allow for a good quality of<br />
the analysis even for element sizes as small as 4 pixels.<br />
0.13<br />
3.5<br />
Mean gray level fluctuation<br />
0.12<br />
0.11<br />
0.1<br />
0.09<br />
0.08<br />
0.07<br />
x (pixel)<br />
3.0<br />
2.5<br />
2.0<br />
1.5<br />
0.06<br />
2 4 8 16 32 64<br />
Element size (pixel)<br />
1.0<br />
2 4 8 16 32 64<br />
Element size (pixel)<br />
-a-<br />
-b-<br />
Figure 3: Fluctuation of gray level values averaged over elements of different sizes normalized<br />
by the gray level dynamic range of the image (a). Average of largest (+) and smallest (◦)<br />
correlation radii determined on elements of varying sizes (b).<br />
Another significant criterion is the correlation radius of the image texture. The latter is<br />
computed from a parabolic interpolation of the auto-correlation function at the origin. The<br />
inverse of the two eigenvalues of the curvature give an estimate of the two correlation radii, ξ 1<br />
and ξ 2 , shown in Fig. 3-b when averaged over all elements of a given size. The texture is rather<br />
isotropic (i.e., similar eigenvalues), and remains small (varying from 1-2 to about 3 pixels) for<br />
all element sizes. This indicates a very good texture quality that reveals small scale details<br />
even for small element sizes. If one wants at least one disk and its complementary surrounding<br />
to get a good estimate, the correlation radii should be less than one fourth of the element size.<br />
This is achieved for element sizes greater than 6 pixels in the present case.<br />
11
How to get these results?<br />
To show the difference with a natural texture, another picture will be considered. It corresponds<br />
to that of a stone wool sample [10]. Choose the Texture option in the first menu. Another<br />
menu appears in which the format of the pictures has to be given (Fig. 4):<br />
• Unknown format: it should one of the following formats .bmp, .CR2, .hbf, .hmf, .jpg,<br />
.png, .tif.<br />
• Image .bmp: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .bmp files.<br />
• Image .hbf or Image .hmf: these are HOLO3 (www.holo3.com) formats.<br />
• Image .jpg: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .jpg files.<br />
• Image .png: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .png files.<br />
• Canon EOS 350: this is a .CR2 raw file.<br />
• Image .tif: this is a classical 8-bit or 16-bit coded format. Make sure the pictures are<br />
stored as B/W .tif files.<br />
Figure 4: Menu to choose the format of the picture and then the reference picture.<br />
The reference image, which is not necessarily located in the same directory as <strong>CORRELI</strong> Q4 ,<br />
has to be chosen (Fig. 4). The next step consists in selecting the region of interest (ROI). Three<br />
options are possible:<br />
12
• Computation Click. Follow the instructions and click to choose the two end points of<br />
the ROI (Fig. 5).<br />
Figure 5: Choice of the Region Of Interest (ROI) by mouse click. Menu to choose the result<br />
file of a previous computation for which the same ROI will be considered.<br />
• Computation Restart. Indicate the results (.mat) file in which the ROI size is given<br />
(Fig. 5).<br />
• Computation Data. In the MATLAB Command Window, the user has to answer to four<br />
questions to choose the size of the ROI. The minimum and maximum values are given<br />
and they correspond to the image size:<br />
minimum horizontal coordinate 1
x<br />
Gray level histogram<br />
104<br />
2.5<br />
Region of interest<br />
2<br />
50<br />
Number of samples<br />
1.5<br />
1<br />
100<br />
150<br />
0.5<br />
200<br />
0<br />
120 130 140 150 160 170 180<br />
Gray level<br />
250<br />
50 100 150 200 250<br />
Figure 6: Histogram of the chosen ROI.<br />
believed that the measurement is not possible (i.e., there are not enough gradients to capture<br />
displacements). With this limit, it is concluded that about 30% of all 4-pixel elements do not<br />
meet this criterion. Such small size should therefore not be used. With 16-pixel and larger<br />
elements, this first criterion is always satisfied.<br />
100<br />
90<br />
80<br />
70<br />
Fluctuation criterion<br />
Unvalidated ZOI<br />
Validated ZOI<br />
0.04<br />
0.035<br />
0.03<br />
Mean fluctation<br />
Min. fluctation<br />
Limit<br />
ZOI percentage<br />
60<br />
50<br />
40<br />
30<br />
20<br />
Relative fluctuation<br />
0.025<br />
0.02<br />
0.015<br />
0.01<br />
10<br />
0.005<br />
0<br />
4 8 16 32 64 128<br />
ZOI size (pixels)<br />
0<br />
4 8 16 32 64 128<br />
ZOI size (pixels)<br />
Figure 7: Percentage of validated elements. Minimum and mean relative RMS fluctuations as<br />
functions of the element size in the analyzed picture of Fig. 6.<br />
Last, the principal correlation radii are shown in values expressed in pixels, and normalized<br />
by the element size (Fig. 8). A practical limit is chosen to be at most 25% of the element size.<br />
Above this value, it is believed that the measurement is not secure. With this limit, it is<br />
concluded that about 50% of all 8-pixel elements do not meet this criterion. Such small size<br />
should therefore not be used. With 16-pixel and larger elements, this second criterion is always<br />
satisfied. Last, let us note that the two correlation radii are significantly different. This result<br />
14
indicates that the texture is anisotropic (Fig. 6). The interested reader will find additional<br />
details in Ref. [18] concerning the determination of dominant orientations of a texture.<br />
ZOI percentage<br />
100<br />
90<br />
80<br />
70<br />
60<br />
50<br />
40<br />
30<br />
20<br />
Correlation radius criterion<br />
Unvalidated ZOI<br />
Validated ZOI<br />
Dimensionless correlation radii<br />
5<br />
4.5<br />
4<br />
3.5<br />
3<br />
2.5<br />
2<br />
1.5<br />
1<br />
R1mean<br />
R1max<br />
R2mean<br />
R2max<br />
Limit<br />
Correlation radii (pixels)<br />
20<br />
18<br />
16<br />
14<br />
12<br />
10<br />
8<br />
6<br />
4<br />
R1mean<br />
R1max<br />
R2mean<br />
R2max<br />
10<br />
0.5<br />
2<br />
0<br />
4 8 16 32 64 128<br />
ZOI size (pixels)<br />
0<br />
4 8 16 32 64 128<br />
ZOI size (pixels)<br />
0<br />
4 8 16 32 64 128<br />
ZOI size (pixels)<br />
Figure 8: Percentage of validated elements. Maximum and mean correlation radii as functions<br />
of the element size.<br />
With these two simple criteria, it is concluded that at least 16-pixel elements are to<br />
be chosen (i.e., both are satisfied simultaneously). It is worth noting that this first a priori<br />
analysis concerns only the texture itself. It is therefore a qualitative analysis. In the following,<br />
two quantitative analyses are proposed.<br />
3.1.2 Displacement uncertainty<br />
Prior to any computation, it is important to estimate the a priori performance of the approach<br />
on the actual texture of the image. If one changes the picture, one may not get exactly the<br />
same performance since it is related to the local details of the gray level distribution as shown<br />
in Section 2. This is performed by using the original image f only, and generating a translated<br />
image g by a prescribed amount u pre . Such an image is generated in Fourier space using a simple<br />
phase shift for each amplitude. This procedure implies a specific interpolation procedure for<br />
inter-pixel gray levels, to which one resorts systematically (see Section 2.3). The algorithm is<br />
then run on the pair of images (f, g), and the estimated displacement field u est (x) is measured.<br />
One is mainly interested in sub-pixel displacements, where the main origin of errors comes<br />
from inter-pixel interpolation. Therefore the prescribed displacement is chosen along the (1, 1)<br />
direction so as to maximize this interpolation sensitivity. To highlight this reference to the<br />
pixel scale, one refers to the x- (or y-) component of the displacement u pre ≡ u pre .e x varying<br />
from 0 to 1 pixel, rather than the Euclidian norm (varying from 0 to √ 2 pixel).<br />
The quality of the estimate is characterized by two indicators, namely, the systematic<br />
error, δ u = ‖〈u est 〉 − u pre ‖, and the standard uncertainty σ u = 〈‖u est − 〈u est 〉‖ 2 〉 1/2 . The<br />
change of these two indicators is shown in Fig. 9 as functions of the prescribed displacement<br />
amplitude for different element sizes l ranging from 4 to 128 pixels. Both quantities reach a<br />
15
maximum for one half pixel displacement, u pre = 0.5 pixel, and are approximately symmetric<br />
about this maximum. Integer valued displacements (in pixels) imply no interpolation and are<br />
exactly captured through the multi-scale procedure discussed above. This confirms that these<br />
errors are due to interpolation procedures. The results are shown in a semi-log scale to reveal<br />
the strong sensitivity to the element size, however a linear scale would show that both δ u and σ u<br />
follow approximately a linear increase with u pre from 0 to 0.5 pixel (and a symmetric decrease<br />
from 0.5 to 1 pixel).<br />
10 -1 0 0.2 0.4 0.6 0.8 1<br />
Mean displacement error (pixel)<br />
10 -2<br />
10 -3<br />
10 -4<br />
10 -5<br />
10 -6<br />
4<br />
8<br />
16<br />
32<br />
64<br />
128<br />
Standard uncertainty (pixel)<br />
10 0 0 0.2 0.4 0.6 0.8 1<br />
10 -1<br />
10 -2<br />
10 -3<br />
10 -4<br />
10 -5<br />
4<br />
8<br />
16<br />
32<br />
64<br />
128<br />
Prescribed displacement (pixel)<br />
Prescribed displacement (pixel)<br />
-a-<br />
-b-<br />
Figure 9: Mean error δ u and standard deviation σ u as a function of the prescribed displacement<br />
u pre for different element sizes l ranging from 4 to 128 pixels.<br />
To quantify the effect of the element size, the error and standard uncertainty, are averaged<br />
over u pre within the range [0, 1] as functions of the element size l. These data are shown in<br />
Fig. 10. A power-law decrease<br />
〈σ u 〉 = A 1+ζ l −ζ<br />
(14)<br />
〈δ u 〉 = B 1+υ l −υ<br />
for 8 ≤ l ≤ 128 pixels is usually observed as shown by a regression line on the graph. Both<br />
amplitudes are typically close to 1 pixel (more precisely A = 1.15 pixel and B = 1.07 pixel).<br />
The exponents are measured to be ζ = 1.96 and υ = 2.34. The data for l = 128 pixels seem<br />
to depart from the power-law trend with a tendency to saturate. These results quantify the<br />
trade-off the experimentalist has to face in the analysis of a displacement field, namely, either<br />
the measurement is accurate but estimated over a large zone, or it is spatially resolved but at<br />
the cost of a less accurate determination. This is a significant difference with classical finite<br />
element techniques for which convergence is achieved when the element size decreases. This<br />
is not the case when measurements are concerned. Let us however underline the following<br />
conclusions:<br />
16
©<br />
<br />
<br />
<br />
<br />
<br />
<br />
;<br />
<br />
<br />
<br />
<br />
:<br />
• Elements as small as l = 4 pixels may be used with an average error and standard<br />
uncertainty of the order of 0.1 pixel,<br />
• Systematic errors of the order of 10 −2 and 10 −3 pixel is reached for element sizes respectively<br />
equal to 8 and 16 pixels.<br />
• Standard uncertainties of the order of 2 × 10 −2 and 6 × 10 −3 pixel is reached for element<br />
sizes respectively equal to 8 and 16 pixels.<br />
• The systematic error in the determination of a displacement is such that evaluations will<br />
be “attracted” toward integer values. Correspondingly, transitions at half-integer pixel<br />
values for the displacement will appear as more abrupt. This phenomenon will be referred<br />
to as “integer locking” in the sequel, and will be discussed in detail. Let us underline<br />
that this spurious bias is revealed in this technique because the latter is used down to<br />
extremely small element sizes.<br />
¡ ¢ ¨<br />
¡£¢ ¨<br />
465 7 1'1"2-3<br />
3 465 7 89<br />
<br />
<br />
<br />
<br />
¡£¢§<br />
¡£¢§<br />
<<br />
¡£¢ ¦<br />
¡£¢ ¦<br />
¡£¢¥<br />
<br />
¡£¢ ¤<br />
¡£¢¥<br />
¡ ¨ ¡ §<br />
¡ ¨ ¡ §<br />
"!$# % &')( *+% ,- .<br />
"!$# % &'=( *+% ,- .<br />
/ >$/<br />
/ 0'/<br />
Figure 10: Average error 〈δ u 〉 and standard uncertainty 〈σ u 〉 as functions of the element size<br />
l. For the displacement uncertainty, the results obtained by Q4-DIC are compared with those<br />
obtained by FFT-DIC. The dashed lines correspond to power-law fits.<br />
For comparison purposes, the displacement uncertainties obtained with the present technique<br />
are compared with those of a standard FFT-DIC technique [10]. In that case, a weaker<br />
power-law decrease is observed with A = 1.00 pixel and ζ = 1.23 (Fig. 10b). This result shows<br />
that by using a continuous description of the displacement field, it enables for a decrease of the<br />
displacement uncertainty when the same element size is used. Conversely, for a given displacement<br />
uncertainty, the Q4-DIC algorithm allows one to reduce significantly the element size,<br />
17
thereby increasing the number of measurement points when compared to a classical FFT-DIC<br />
technique.<br />
How to get these results?<br />
The same instructions as presented for the texture analysis (Section 3.1.1) hold. We reproduce<br />
them for a more linear reading. Choose the Uncertainty option in the first menu. Another<br />
menu appears in which the format of the pictures is chosen (Fig. 4):<br />
• Unknown format: it should one of the following formats .bmp, .CR2, .hbf, .hmf, .jpg,<br />
.png, .tif.<br />
• Image .bmp: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .bmp files.<br />
• Image .hbf or Image .hmf: these are HOLO3 (www.holo3.com) formats.<br />
• Image .jpg: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .jpg files.<br />
• Image .png: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .png files.<br />
• Canon EOS 350: this is a .CR2 raw file.<br />
• Image .tif: this is a classical 8-bit or 16-bit coded format. Make sure the pictures are<br />
stored as B/W .tif files.<br />
The reference image, which is not necessarily located in the same directory as <strong>CORRELI</strong> Q4 ,<br />
has to be chosen (Fig. 4). The next step consists in selecting the region of interest (ROI). Three<br />
options are possible:<br />
• Computation Click. Follow the instructions and click to choose the two end points of<br />
the ROI (Fig. 5).<br />
• Computation Restart. Indicate the results (.mat) file in which the ROI size is given<br />
(Fig. 5).<br />
• Computation Data. In the MATLAB Command Window, the user has to answer to four<br />
questions to choose the size of the ROI. The minimum and maximum values are given<br />
and they correspond to the image size:<br />
18
minimum horizontal coordinate 1
`U X<br />
^ T<br />
R _<br />
]\\X<br />
\<br />
ZX<br />
[<br />
U VW<br />
YX<br />
<br />
¢¡¤£ ¥¡§¦ ¢¡©¨<br />
¢¡¤ <br />
K Y<br />
Q ZN<br />
T W<br />
X M<br />
VQ<br />
TO<br />
K S<br />
J K LM<br />
N OP<br />
¢¡¤£<br />
! #"#$!%'&)( *+-, .0/¤1325476<br />
¢¡¤£ ¥¡§¦ ¢¡©¨<br />
¢¡¤ <br />
!.0/¤132405'6 798§:=<br />
"!#$<br />
%'&)( *+-,<br />
<br />
[]\§^ _`_<br />
¢¡¤¦<br />
¢¡¤ ¨<br />
¢¡¤¦<br />
T U<br />
PS<br />
RQ<br />
SQ<br />
¢¡ <br />
¢¡ ¨<br />
Q R ST<br />
?@BACDFE¤GBHIJDFKLGNMOCP<br />
89;:=@?¤A;BCD=@EFAHG#I<br />
Figure 12: Mean displacement error and standard deviation as functions of the element size.<br />
3.1.3 Noise sensitivity<br />
Last, the effect of noise associated to the image acquisition (e.g., digitization, read-out noise,<br />
black current noise, photon noise [20]) on the displacement measurement is assessed. This analysis<br />
allows one to estimate the displacement resolution [21]. The reference image is corrupted<br />
by a Gaussian noise of zero mean and standard variation σ g ranging from 1 to 8 gray levels<br />
at each pixel with no spatial correlation. No displacement field is superimposed on the image,<br />
and the displacement field is then estimated. The standard deviation of the displacement field,<br />
σ u , is shown in Fig. 13 as a function of the noise amplitude σ g and for different element sizes<br />
ranging from 4 to 128 pixels. The quantity σ u is linear in the noise amplitude and inversely proportional<br />
to the element size. The latter properties are derived from the central limit theorem.<br />
A theoretical analysis of this problem is discussed in Refs. [11, 22], and leads to the<br />
following estimate of the standard deviation of the displacement field induced by a Gaussian<br />
white noise<br />
σ u =<br />
12√ 2σ g p<br />
7〈|∇f| 2 〉 1/2 l<br />
where p is the physical pixel size. For the present application, one computes 〈|∇f| 2 〉 1/2 ≈<br />
5340 pixel −1 , hence σ u ≈ 4.5 × 10 −4 σ g p/l. This theoretical expectation (neglecting the spatial<br />
correlation in the image texture) is consistent with the direct estimates shown in Fig. 13 (e.g.,<br />
for l = 4 pixels and σ g = 8 gray levels, the direct estimate is 1.2×10 −3 pixel to be compared with<br />
9 × 10 −4 pixel given by the above formula). In practice, with the used CCD camera, the noise<br />
level is given with a maximum range less than 3 gray levels. Consequently, the contribution of<br />
(15)<br />
20
-3<br />
1.5 x 10 Noise level (gray level)<br />
Standard uncertainty (pixel)<br />
1<br />
0.5<br />
0<br />
0 2 4 6 8<br />
4<br />
8<br />
16<br />
32<br />
64<br />
128<br />
Figure 13: Standard deviation of the displacement error versus noise amplitude for different<br />
element sizes (4, 8, 16, 32, 64 and 128 pixels) from top to bottom.<br />
image noise is negligibly small when compared to that induced by the sub-pixel interpolation.<br />
How to get these results?<br />
The same instructions as presented for the texture analysis (Section 3.1.1) hold. We reproduce<br />
them for a more linear reading. Choose the Resolution option in the first menu. Another<br />
menu appears in which the format of the pictures has to be chosen (Fig. 4):<br />
• Unknown format: it should one of the following formats .bmp, .CR2, .hbf, .hmf, .jpg,<br />
.png, .tif.<br />
• Image .bmp: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .bmp files.<br />
• Image .hbf or Image .hmf: these are HOLO3 (www.holo3.com) formats.<br />
• Image .jpg: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .jpg files.<br />
• Image .png: this is a classical 8-bit coded format. Make sure the pictures are stored as<br />
B/W .png files.<br />
• Canon EOS 350: this is a .CR2 raw file.<br />
• Image .tif: this is a classical 8-bit or 16-bit coded format. Make sure the pictures are<br />
stored as B/W .tif files.<br />
21
The reference image, which is not necessarily located in the same directory as <strong>CORRELI</strong> Q4 ,<br />
has to be chosen (Fig. 4). The next step consists in selecting the region of interest (ROI). Three<br />
options are possible:<br />
• Computation Click. Follow the instructions and click to choose the two end points of<br />
the ROI (Fig. 5).<br />
• Computation Restart. Indicate the results (.mat) file in which the ROI size is given<br />
(Fig. 5).<br />
• Computation Data. In the MATLAB Command Window, the user has to answer to four<br />
questions to choose the size of the ROI. The minimum and maximum values are given<br />
and they correspond to the image size:<br />
minimum horizontal coordinate 1
92<br />
+<br />
65<br />
0<br />
)*+<br />
*-,<br />
./+<br />
32<br />
,<br />
),+'<br />
,<br />
),+'<br />
¦¥£<br />
¦¥¤<br />
¦¥¡<br />
¦¥¢<br />
§<br />
4 7<br />
6:;5<br />
5<br />
67<br />
§¦ ¢<br />
. 8<br />
234<br />
-1<br />
7 1<br />
§¦ <br />
§¦ ¡<br />
/0<br />
%.<br />
2 3-+<br />
+4<br />
-'<br />
%&<br />
§¦ ©¨<br />
0.,1<br />
(*<br />
$%&<br />
(&' )<br />
8:9<br />
¡ ¢ £ ¤<br />
¥ §¦ ©¨<br />
74; 695<br />
¡ ¢ £ ¤<br />
¨©¦ !¦#"<br />
§! #"%$&§('<br />
Figure 14: Mean displacement error and standard deviation as functions of the noise level for<br />
different element sizes.<br />
4 5<br />
§¡©¤<br />
687<br />
/0<br />
-1<br />
%.<br />
-'<br />
%&<br />
§¡ ¨£<br />
9 :<br />
(*<br />
$%&<br />
(&' )<br />
¢¡£ ¤ ¤¥¡£ ¦<br />
©§¥ !©§#"<br />
Figure 15: Mean displacement error and standard deviation as functions of the noise level for<br />
different element sizes.<br />
low noise features. Otherwise, it might be desirable to increase the element size, or work harder<br />
to get pictures with a better quality (this is, unfortunately not always possible...).<br />
3.2 Computation<br />
Any of the following options has to be selected to run directly a computation:<br />
• Computation Click. Follow the instructions and click to choose the two end points of<br />
the ROI (Fig. 5).<br />
• Computation Restart. Indicate the results (.mat) file in which the ROI size is given<br />
(Fig. 5).<br />
• Computation Data. In the MATLAB Command Window, the user has to answer to four<br />
23
questions to choose the size of the ROI. The minimum and maximum values are given<br />
and they correspond to the image size:<br />
minimum horizontal coordinate 1
• The number of scales is the second most important parameter. As discussed above, the<br />
present algorithm is based upon a coarse-graining approach. When large displacements<br />
and / or strains are suspected to occur, it is desirable to increase the number of scales.<br />
The maximum number of scales depends on the size of the ROI in comparison to that<br />
of the elements (e.g., for a 512 × 512pixel-ROI, and 16-pixels elements, the maximum<br />
number of scales is 5, i.e., at least 2 × 2 elements are needed for the coarsest scale).<br />
• The number of iterations (the higher the number, the higher the accuracy and the computation<br />
time). This value is given to avoid endless iterations. This criterion should not<br />
be reached. In case it is reached a message warns the user:<br />
--------------------------- WARNING - WARNING -----------------------------<br />
Warning: No convergence: check results!!<br />
--------------------------- WARNING - WARNING -----------------------------<br />
• The number of images to analyze (the reference image is not included in that count);<br />
• When a sequence of more than one image is considered, the correlation can be performed<br />
either by considering always the same reference image for strains (in absolute value) less<br />
than or equal to 10% or by changing (or updating) the reference image. In the last case,<br />
the Image update Y/N has to be activated.<br />
• Last, it is possible to store the correlation results for any couple of analyzed pictures.<br />
This is made possible by activating the Independent calculation Y/N button.<br />
Validate the choice at the end (press the Validation button).<br />
The sequence of deformed pictures has to be selected. The user selects the whole sequence<br />
by using the same procedure as for the reference picture (Fig. 17).<br />
It is then possible to mask part of the ROI. A menu appears and seven different operations<br />
are possible:<br />
• Define exclusion polygon allows the user by mouse click to enter different exclusion<br />
regions. Follow the instructions on top of the picture that appears.<br />
• Remove polygon allows the user by mouse click to remove any of the existing exclusion<br />
polygon.<br />
25
• Define exclusion circle allows the user, for instance to mask a hole.<br />
instructions on top of the picture that appears.<br />
Follow the<br />
• Remove exclusion circle allows the user by mouse click to remove any of the existing<br />
exclusion circle.<br />
• Define inclusion circle allows the user to choose a circular ROI from a rectangular<br />
shape chosen before. This would be the case for the analysis of a Brazilian test [12].<br />
Follow the instructions on top of the picture that appears.<br />
• Remove inclusion circle allows the user by mouse click to remove any of the existing<br />
inclusion circle.<br />
• Redraw when the user is not happy with the (random) colors.<br />
• Exit to end the mask procedure.<br />
When the computation is completed, the result file has to be saved (Fig. 18). The extension<br />
is ‘.mat’ to be readable for a later visualization. The computation is now ended and the<br />
visualization stage starts. One can choose to perform another computation or to visualize any<br />
Figure 16: dialog box for choosing the correlation parameters.<br />
26
Figure 17: dialog box to choose the deformed picture(s).<br />
Figure 18: dialog box to choose the type of mask and to save the results<br />
results. Type the command correli_q4 at the MATLAB prompt. It can be noted that during<br />
the computations, different messages may appear. Some messages are only given to indicate<br />
that the computation is running normally.<br />
3.3 Visualization<br />
If Visualization is chosen, the result file (‘.mat’ extension) to be displayed has to be selected<br />
(Fig. 19).<br />
Figure 20 appears, in which different options can be chosen. In the present case, the<br />
vertical displacement field (expressed in pixels) is plotted for the stone wool texture for which<br />
an artificially deformed picture was created. A uniform nominal strain level of 0.25 was applied.<br />
This is an extreme case for the present correlation software that needs all the scales to capture<br />
properly the displacement field. The element size was equal to 16 pixels. It is worth noting<br />
that the maximum displacements are greater than 3 times the element size. Had the multiscale<br />
algorithm not been implemented, it would have been impossible to capture these levels (i.e., 5<br />
27
Figure 19: dialog box to read the results of a previous computation.<br />
scales were used in the present case). When a sequence is analyzed, any image can be selected<br />
Figure 20: Visualization of the measurement results. In the present case, the displacement field<br />
along the vertical direction.<br />
by using the two buttons bellow the Image No. message (Fig. 20).<br />
28
At the top left corner are options related to the strain measures:<br />
• infinitesimal: infinitesimal strains (i.e., symmetric part of displacement gradient, see<br />
Eq. (24));<br />
• nominal: nominal (Cauchy-Biot) strains (Eq. (58), when m = 1/2);<br />
• Green Lag.: Green-Lagrange strains (Eq. (58), when m = 1);<br />
• logarithmic: logarithmic (Hencky) strains (Eq. (58), when m → 0 + );<br />
• RdB (internal development);<br />
• Eigen value?: the eigen values of the selected strain measure are shown.<br />
It is worth noting that in the case of large strains (see details in the Appendix), the out-ofplane<br />
displacement is assumed to be small and is therefore neglected. Other assumptions can<br />
be made. They will have to be implemented on demand.<br />
Just below are options related to the type of component to visualize. The type of component<br />
to display has to be chosen, namely, in-plane strains, in-plane displacements, or out-ofplane<br />
rotation. The frame is always the same, namely, horizontal direction: 2, vertical direction:<br />
1. By choosing error, an error indicator (i.e., correlation residual |u(x).∇f(x) + f(x) − g(x)|<br />
normalized by the dynamic range of the picture) of the result is plotted. The closer to 0, the<br />
better the result (when no lighting variations occur, levels below few percentages are usually<br />
achieved).<br />
Two pictures are plotted. The left picture is always the very first reference picture. For<br />
the reference image, the chosen field is plotted. In the case of displacements, the average value<br />
is always in the middle of the scale. To change the scale, the dialog box in between the two<br />
images is to be used (Fig. 21).<br />
The number of contours can be changed. If 0 is set, a grayscale is used. If 11 is chosen,<br />
a fancy color scale appears, otherwise a conventional (i.e., hot) one is used. The amplitude of<br />
the displacement field is chosen with respect to the average value. For strains, two routes can<br />
be followed:<br />
• the maximum and minimum values are given by the user in the middle part of the dialog<br />
box;<br />
• the w or w/o mean button is activated and the average value corresponds to the middle<br />
of the scale, the range of which is chosen as the maximum strain.<br />
29
Figure 21: Contour and scale options. Mesh options. Amplification of the displacements.<br />
When the fill option is activated, the contours are filled on the reference image. The rigid<br />
body motion can be subtracted to the overall displacement by pushing the corresponding button<br />
(Rigid body motion Y/N). When the error indicator is plotted, there is no need to change any<br />
scale parameters, it is performed automatically. The w or w/o mean button can also be used.<br />
When selecting mesh, the undeformed and deformed meshes appear on the relevant images<br />
(Fig. 22). When selecting vector, the displacement vectors are shown (Fig. 21). It can be noted<br />
that both options can be used simultaneously. A slider enables the displacements to be amplified<br />
by a factor that can be chosen by the user (Fig. 21). When the amplification is greater than 1,<br />
the underlying image disappears.<br />
To further comment on the results on the artificially deformed picture, the corresponding<br />
(nominal) strain field is shown in Fig. 22. Since the w or w/o mean button is activated, the<br />
average strain can be read directly as the median value (i.e., 0.2396 ≈ 0.24 for a prescribed<br />
value of 0.25).<br />
3.4 Virtual Gauges<br />
Type the command gauge at the MATLAB prompt. This procedure allows for the computation<br />
of the average strain in a user-selected ROI. one needs to indicate the result file in which all<br />
the data needed are stored. Then the Region Of Interest (ROI) is selected by mouse. Click and<br />
maintain to select the ROI. If the ROI is larger than the mesh, then the size is automatically<br />
reset to the maximum size.<br />
In the MATLAB Command Window, average values of different strain tensors are given.<br />
30
Figure 22: Visualization of post-processed results. In the present case the normal (nominal)<br />
strains along the vertical direction.<br />
The corresponding eigen strains are also computed. The results can also be saved in an ASCII<br />
file. The computation is now ended. Type any command at the MATLAB prompt.<br />
4 Application to a tensile test<br />
In this section, an application of the previously proposed algorithm is carried out to analyze<br />
a tensile test performed on an aluminum alloy sample. In the plastic regime, the formation<br />
of localization strain bands is observed. The fact that for a given displacement uncertainty,<br />
smaller element sizes can be chosen in the present case (Q4-DIC) when compared with those<br />
of a standard FFT-DIC technique (Fig. 10b), enables one to better capture kinematic details<br />
in the localization band.<br />
31
4.1 Material and method<br />
The studied aluminum alloy is of type 5005 (i.e., more than 99 wt% of Al content and a small<br />
amount of Mg; these values were determined by electron probe micro-analysis). As shown in<br />
Fig. 2, the sample is a coated with sprayed black and white paints to create the random texture<br />
for the displacement field measurement. The sample size is 140 × 30 × 2 mm 3 . It is positioned<br />
within hydraulic grips of a 100 kN servo-hydraulic testing machine. Its alignment is checked<br />
with DIC measurements (i.e., no significant rotation of the sample is observed in the elastic<br />
domain). To have a first strain evaluation, an extensometer was used. Its pins are observed on<br />
the right edge of the sample (Fig. 2).<br />
An artificial light source is used to minimize gray level variations so that the conservation<br />
of the optical flow is considered as practically achieved. A CCD camera (12-bit digitization,<br />
noise less than 3 gray levels, resolution: 1024 × 1280 pixels) with a conventional zoom is<br />
positioned in front of the sample. In the present case, the physical size of one pixel is 25 µm.<br />
Two loading sequences are carried out. First, in the elastic domain, a controlled displacement<br />
rate of 5 µm/s is applied and pictures are taken for 12 µm-increments. Elastic properties may<br />
be identified [23]. This will not be discussed herein. Second, a controlled displacement rate of<br />
10 µm/s is applied to study strain localization and pictures are taken for 60 µm-increments.<br />
The following analysis of the displacement field is an increment between two image acquisitions<br />
in the “plastic” regime.<br />
Figure 23 shows the change of the average longitudinal strain with the number of pictures<br />
(or equivalently with time). This result was obtained by using the Q4-DIC analysis. Until the<br />
extensometer pins slipped (at about a 5% strain), the average strain measured by DIC and<br />
that given by the extensometer were close, even though the same surface was not analyzed.<br />
This response is typical of a Portevin-Le Châtelier (PLC) phenomenon or jerky flow [24].<br />
From a microscopic point of view, PLC effects are related to dynamic interactions between<br />
mobile dislocations and diffusing solute atoms [25, 26]. From a macroscopic perspective, it is<br />
related to a negative strain rate sensitivity that leads to localized bands that are simulated [27].<br />
Many experimental studies [28] however are based upon average strain measurements. There<br />
are also full-field displacement measurements performed by using, for instance, laser speckle<br />
interferometry [29]. Yet the spatial resolution did not allow for an analysis of the displacement<br />
field within the band. Additional insight is gained by using IR thermography [30].<br />
32
¢¡ ¨§<br />
¢¡ ¦©<br />
¢¡ ¨<br />
¢¡ ¨<br />
¢¡ <br />
<br />
¢¡ ¦¥<br />
¢¡ ¤£<br />
¥ © £ £¥ £©<br />
¤! "$#¨%'&)(#*,+¢&-%<br />
Figure 23: Mean strain for a region of interest of 1000 × 700 pixels as a function of the number<br />
of picture. The box shows the two pictures that are analyzed.<br />
4.2 Kinematic fields<br />
Let us now analyze the displacement field in between two states (0.3% mean strain apart,<br />
see Fig. 23). The same region of interest of size 1000 × 700 pixels is studied using the above<br />
method, with different element sizes ranging from 16 down to 4 pixels. Figures 24 and 25 show<br />
the resulting displacement fields (component U x and U y , respectively). Although the test is<br />
pure tension, the analysis reveals without ambiguity the presence of a localization band whose<br />
width is about 150 pixels, and across which the displacement discontinuity is about 2 pixels<br />
along the tension axis, and about 1 pixel perpendicular to it. Let us concentrate here on this<br />
single pair of images to validate the algorithm on a real experimental test and evaluate its<br />
performances.<br />
One notes that all element sizes may be used. As expected, the smallest element sizes are<br />
noisier, yet the agreement between all these determinations is excellent. Let us underline that<br />
FFT-DIC usually deals with element sizes equal or larger than 32 pixels, exceptionally 16 pixels<br />
for very favorable cases are used when locally constant displacement fields are sought. Using<br />
the Q4-DIC technique allows one to reduce the element size by a factor of 4, which means that<br />
the number of pixels in the element has been cut down by a factor of 16.<br />
Let us however note that one should be cautious about the fact that displacements have a<br />
tendency to be attracted toward integer values, especially for small element sizes. Therefore the<br />
direct evaluation of strains along the tension axis ε yy as obtained by the Q4P1-shape functions<br />
or equivalently as a simple finite difference is expected to be artificially increased at half-integer<br />
displacement components. Figure 26 shows such strain fields for 4 different element sizes from<br />
33
100<br />
100<br />
200<br />
200<br />
300<br />
300<br />
400<br />
400<br />
500<br />
500<br />
600<br />
600<br />
700<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-a-<br />
-b-<br />
¢¡ £¥¤ ¦¨§ © <br />
100<br />
200<br />
300<br />
400<br />
500<br />
<br />
100<br />
0.4<br />
200<br />
0.2<br />
300<br />
0<br />
400<br />
0.2<br />
500<br />
600<br />
0.4<br />
600<br />
700<br />
0.6<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-c-<br />
-d-<br />
100<br />
100<br />
200<br />
200<br />
300<br />
300<br />
400<br />
400<br />
500<br />
500<br />
600<br />
600<br />
700<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-e-<br />
-f-<br />
Figure 24: Map of U y displacement for different element sizes: (a) l = 16, (b) l = 12, (c)<br />
l = 10, (d) l = 8, (e) l = 6 and (f) l = 4 pixels. The physical size of one pixel is equal to<br />
25 µm.<br />
16 down to 8 pixels. For a size of 16 pixels, the localization band appears as a genuine zone<br />
of increased strains as compared to a “silent” (or elastic) background. For smaller element<br />
sizes, the edges of the shear band appear to concentrate still a higher strain. The same effect<br />
is apparent for element sizes 12, 10 and 8 pixels. The strain maps obtained for smaller element<br />
sizes are not shown, since the noise level becomes much higher and thus the measurement<br />
cannot be trusted. The same artefact of strain enhancement at the edges of the shear band is<br />
however observed.<br />
34
100<br />
100<br />
200<br />
200<br />
300<br />
300<br />
400<br />
400<br />
500<br />
500<br />
600<br />
600<br />
700<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-a-<br />
-b-<br />
£ ¥¦ ¨© <br />
100<br />
100<br />
¢¡¤£ ¥§¦ ¨§© <br />
200<br />
300<br />
400<br />
500<br />
600<br />
700<br />
2<br />
1.5<br />
1<br />
0.5<br />
0<br />
200<br />
300<br />
400<br />
500<br />
600<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-c-<br />
-d-<br />
100<br />
100<br />
200<br />
200<br />
300<br />
300<br />
400<br />
400<br />
500<br />
600<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
500<br />
600<br />
700<br />
800<br />
200 400 600<br />
-e-<br />
-f-<br />
800 1000<br />
Figure 25: Map of U x displacement for different element sizes: (a) l = 16, (b) l = 12, (c)<br />
l = 10, (d) l = 8, (e) l = 6 and (f) l = 4 pixels. The physical size of one pixel is equal to<br />
25 µm.<br />
4.3 Integer locking<br />
One notes on the previous figure that the U y -displacement is half-integer valued at the edge of<br />
the shear band. The larger strains at the edge of the band could therefore be interpreted as an<br />
artefact due to integer locking. Integer valued displacements being favored, an artificial gap is<br />
created for half-integer values displacement, and thus any gradient (finite difference operator)<br />
will underline this effect very markedly.<br />
To test this interpretation, the following test is proposed. An artificially translated image<br />
by 0.5 pixel is computed from the original one, using a fast Fourier transform, as the latter<br />
35
100<br />
100<br />
200<br />
200<br />
300<br />
300<br />
400<br />
400<br />
500<br />
500<br />
600<br />
600<br />
700<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-a-<br />
-b-<br />
y (pixel)<br />
ε xx<br />
x (pixel)<br />
100<br />
200<br />
300<br />
0.025<br />
0.02<br />
0.015<br />
100<br />
200<br />
300<br />
400<br />
500<br />
0.01<br />
0.005<br />
0<br />
400<br />
500<br />
600<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
0.005<br />
0.01<br />
0.015<br />
600<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
-c-<br />
-d-<br />
Figure 26: Map of the strain component ε xx for different element sizes: (a) l = 16, (b) l = 12,<br />
(c) l = 10 and (d) l = 8 pixels.<br />
provides a simple and numerically efficient way of interpolating the image at any arbitrary subpixel<br />
value. A genuine strain enhancement is thus expected to be identified at a fixed position<br />
in the reference image frame of coordinates, whereas a numerical artefact would be moved to<br />
a different location. Figures 27a and b show the U y displacement component starting from<br />
the original image or from the translated one (and where the 1/2 pixel translated has been<br />
corrected for). A good agreement is observed for the displacement field thus revealing a rather<br />
poor sensitivity to such a rigid translation. Figures 27c and d show the corresponding ε yy strain<br />
maps. On the latter set of figures, although high strain values tend to concentrate along two<br />
lines in both figures, the precise location of these bands is not stable. This is a signature of the<br />
integer locking phenomenon. Therefore the strain enhancement at the edge of the shear band<br />
is to be considered as an artefact.<br />
Let us underline that such a phenomenon results from the fact that the elements are<br />
reduced to a very small size, and still provide a very accurate determination of the displacement<br />
field, without much noise. Such a success encourages the user to decrease the size of the elements<br />
to very small values. By doing so, the determination of the displacement is much more prone<br />
36
U x (pixel)<br />
100<br />
200<br />
0<br />
100<br />
200<br />
300<br />
0.5<br />
300<br />
400<br />
500<br />
600<br />
1<br />
1.5<br />
400<br />
500<br />
600<br />
700<br />
2<br />
700<br />
800<br />
200 400 600<br />
800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-a-<br />
-b-<br />
y (pixel)<br />
ε xx<br />
x (pixel)<br />
100<br />
200<br />
0.04<br />
0.03<br />
100<br />
200<br />
300<br />
400<br />
500<br />
600<br />
0.02<br />
0.01<br />
0<br />
0.01<br />
0.02<br />
300<br />
400<br />
500<br />
600<br />
700<br />
0.03<br />
700<br />
800 0.04<br />
200 400 600 800 1000<br />
800<br />
200 400 600<br />
800 1000<br />
-c-<br />
-d-<br />
Figure 27: Map of U x displacement, estimated for a element size l = 12 pixels: (a) for the<br />
original image, (b) for an artificially translated image of 0.5 pixel. Maps of the corresponding<br />
normal strain component ε xx : (c) for the original image, (d) for an artificially translated image<br />
of 0.5 pixel. The physical size of one pixel is equal to 25 µm.<br />
to slight sub-pixel shifts (Fig. 10), here characterized as an attraction toward integer values,<br />
which appear as very significant upon differentiation (in the computation of the strain). This<br />
phenomenon should be identified before any further interpretation of the strain map to ensure<br />
its validity.<br />
4.4 Error maps<br />
As mentioned earlier, a very important output of the displacement measurement obtained from<br />
a minimization procedure is that the optimization functional provides not only a global quality<br />
factor of the determined field, but more importantly a spatial map of residuals, so that one<br />
may appreciate a specific problem that may be spatially localized.<br />
Figure 28 shows the different error maps obtained for different element sizes. This error<br />
is the remanent difference in gray levels that is still unexplained by the estimated displacement<br />
field. The first observation is that the error level does not vary very much with the element<br />
37
size. This is consistent with the fact that the displacement field is quite comparable for different<br />
element sizes. However, there is a slight increase in the error as the element size decreases. This<br />
is explained by the fact that the performance of the correlation algorithm degrades as the spatial<br />
resolution improves. This observation is in good agreement with the results to be expected from<br />
the analysis of Fig. 10.<br />
Figure 28: Map of the residual error η for different element sizes: (a) l = 16, (b) l = 12, (c)<br />
l = 10 (d) l = 8, (e) l = 6 and (f) l = 4 pixels.<br />
5 Applications of <strong>CORRELI</strong> Q4<br />
In the following, different references are given in which the software was used. The abstract is<br />
given for the reader to select the relevant ones. The name of the .pdf file is also given.<br />
38
5.1 Reference paper<br />
G. Besnard, F. Hild and S. Roux, “Finite-element” displacement fields analysis from digital<br />
images: Application to Portevin-Le Châtelier bands, Exp. Mech. 46 (2006) 789-803. EM2.pdf<br />
This is the paper in which the so-called Q4-DIC technique is discussed in details. A new<br />
methodology is proposed to estimate displacement fields from pairs of images (reference and<br />
strained) that evaluates continuous displacement fields. This approach is specialized to a finiteelement<br />
decomposition, therefore providing a natural interface with a numerical modeling of<br />
the mechanical behavior used for identification purposes. The method is illustrated with the<br />
analysis of Portevin-Le Châtelier bands in an aluminum alloy sample subjected to a tensile test.<br />
A significant progress with respect to classical digital image correlation techniques is observed<br />
in terms of spatial resolution and uncertainty.<br />
5.2 Identification of elastic properties<br />
F. Hild and S. Roux, Digital image correlation: from measurement to identification of elastic<br />
properties - A review, Strain 42 (2006) 69-80. Strain1.pdf<br />
The current state of the art of digital image correlation, where displacements can be determined<br />
for values less than one pixel, enables one to better characterize the behavior of materials<br />
and the response of structures to external loads. A general presentation of the extraction of<br />
displacement fields from pictures taken at different instants during an experiment is given.<br />
Different strategies can be followed to determine sub-pixel displacements. New identification<br />
procedures are then devised making use of full-field measurements. A priori or a posteriori<br />
routes can be followed. They are illustrated on the analysis of a Brazilian disk test.<br />
5.3 SIF measurements<br />
S. Roux and F. Hild, Stress intensity factor measurements from digital image correlation: postprocessing<br />
and integrated approaches, Int. J. Fract. 140 [1-4] (2006) 141-157. IJF4.pdf<br />
Digital image correlation is an appealing technique for studying crack propagation in<br />
brittle materials such as ceramics. A case study is discussed where the crack geometry, and<br />
the crack opening displacement are evaluated from image correlation by following two different<br />
measurement and identification routes. The displacement uncertainty can reach the nanometer<br />
range even though optical pictures are dealt with. The stress intensity factor is estimated with<br />
a 7% uncertainty in a complex loading set-up without having to resort to a numerical modeling<br />
of the experiment.<br />
39
R. Hamam, F. Hild and S. Roux, Stress intensity factor gauging by digital image correlation:<br />
Application in cyclic fatigue, Strain 43 (2007) 181-192. Strain2.pdf<br />
A fatigue crack in steel (CCT geometry) is studied via digital image correlation. The<br />
measurement of the stress intensity factor change during one cycle is performed using a decomposition<br />
of the displacement field onto a tailored set of elastic fields. The same analysis<br />
is performed using two different routes, namely, the first one consists in computing the displacement<br />
field using a general correlation technique providing the displacement field projected<br />
onto finite element shape functions, and then analyzing this displacement field in terms of the<br />
selected mechanically relevant fields. The second strategy, called integrated approach, directly<br />
estimates the amplitude of these elastic fields from the correlation of successive images. Both<br />
procedures give consistent results, and offer very good performances in the evaluation of the<br />
crack tip position (uncertainty of about 20 µm for a 14.5-mm crack), and stress intensity factors<br />
(uncertainty less than 1 MPa √ m).<br />
J. Réthoré, S. Roux and F. Hild, Noise-robust Stress Intensity Factor Determination from Kinematic<br />
Field Measurements, Eng. Fract. Mech. (2008) DOI: 10.1016/j.engfracmech.2007.04.018.<br />
EFM1.pdf<br />
Stress Intensity Factors are often estimated numerically from a given displacement field<br />
through an interaction integral formalism. The latter method makes use of a weight, the virtual<br />
crack extension field, which is under-constrained by first principles. Requiring a least noise<br />
sensitivity allows one to compute the optimal virtual crack extension. Mode I and mode II<br />
specialized fields are obtained and particularized for a given displacement functional basis. The<br />
method is applied to an experimental case study of a crack in a silicon carbide sample, whose<br />
displacement field is obtained by a digital image correlation technique. The optimization leads<br />
to a very significant uncertainty reduction up to a factor 100 of the non-optimized formulation.<br />
The proposed scheme reveals additional performances with respect to the integral domain choice<br />
and assumed crack tip geometry, which are shown to have a reduced influence.<br />
5.4 Analysis of elasto-plastic laws<br />
V. Tarigopula, O. S. Hopperstad, M. Langseth, A. H. Clausen, F. Hild, O.-G. Lademo and<br />
M. Eriksson, A Study of Large Plastic Deformations in Dual Phase Steel Using Digital Image<br />
Correlation and FE Analysis, Exp. Mech. 48 [2] (2008) 181-196. EXME3.pdf<br />
Large plastic deformation in sheets made of dual phase steel DP800 is studied experimentally<br />
and numerically. Shear testing is applied to obtain large plastic strains in sheet metals<br />
without strain localisation. In the experiments, full-field displacement measurements are car-<br />
40
ied out by means of digital image correlation, and based on these measurements the strain field<br />
of the deformed specimen is calculated. In the numerical analyses, an elastoplastic constitutive<br />
model with isotropic hardening and the Cockcroft-Latham fracture criterion is adopted to<br />
predict the observed behavior. The strain hardening parameters are obtained from a standard<br />
uniaxial tensile test for small and moderate strains, while the shear test is used to determine<br />
the strain hardening for large strains and to calibrate the fracture criterion. Finite Element<br />
(FE) calculations with shell and brick elements are performed using the non-linear FE code<br />
LS-DYNA. The local strains in the shear zone and the nominal shear stress-elongation characteristics<br />
obtained by experiments and FE simulations are compared, and, in general, good<br />
agreement is obtained. It is demonstrated how the strain hardening at large strains and the<br />
Cockcroft-Latham fracture criterion can be calibrated from the in-plane shear test with the aid<br />
of non-linear FE analyses.<br />
5.5 Identification of a damage law<br />
S. Roux and F. Hild, Digital Image Mechanical Identification (DIMI), Exp. Mech. DOI:<br />
10.1007/s11340-007-9103-3 (2008). EXME4.pdf<br />
A continuous pathway from digital images acquired during a mechanical test to quantitative<br />
identification of a constitutive law is presented herein based on displacement field analysis.<br />
From images, displacement fields are directly estimated within a finite element framework.<br />
From the latter, the application of the equilibrium gap method provides the means for rigidity<br />
field evaluation. In the present case, a reconditioned formulation is proposed for a better stability.<br />
Last, postulating a specific form of a damage law, a linear system is formed that gives a<br />
direct access to the (non-linear) damage growth law in one step. The two last procedures are<br />
presented, validated on an artificial case, and applied to the case of a biaxial tension of a composite<br />
sample driven up to failure. A quantitative estimate of the quality of the determination<br />
is proposed, and in the last application, it is shown that no more than 7% of the displacement<br />
field fluctuations are not accounted for by the determined damage law.<br />
5.6 Analyses of localized phenomena<br />
I. Elnasri, S. Pattofatto, H. Zhao, H. Tsitsiris, F. Hild and Y. Girard, Shock enhancement of<br />
cellular structures under impact loading: Part I Experiments, J. Mech. Phys. Solids 55 (2007)<br />
2652-2671. JMPS2.pdf<br />
This paper aims at showing experimental proof of the existence of a shock front in cellular<br />
structures under impact loading, especially at low critical impact velocities around 50 m/s.<br />
41
First, an original testing procedure using a large diameter Nylon Hopkinson bar is introduced.<br />
With this large diameter soft Hopkinson bar, tests under two different configurations (pressure<br />
bar behind/ahead of the supposed shock front) at the same impact speed are used to obtain the<br />
force/time histories behind and ahead of the assumed shock front within the cellular material<br />
specimen. Stress jumps (up to 60% of initial stress level) as well as shock front speed are measured<br />
for tests at 55 m/s on Alporas foams and nickel hollow sphere agglomerates, whereas no<br />
significant shock enhancement is observed for Cymat foams and 5056 aluminium honeycombs.<br />
The corresponding rate sensitivity of the studied cellular structures is also measured and it<br />
is proven that it is not responsible for the sharp strength enhancement. A photomechanical<br />
measurement of the shock front speed is also proposed to obtain a direct experimental proof.<br />
The displacement and strain fields during the test are obtained by correlating images shot with<br />
a high speed camera. The strain field measurements at different times show that the shock<br />
front discontinuity propagates and allows for the measurement of the propagation velocity. All<br />
the experimental evidences enable us to confirm the existence of a shock front enhancement<br />
even at quite low impact velocities for a number of studied materials.<br />
V. Tarigopula, O. S. Hopperstad, M. Langseth, A. H. Clausen and F. Hild, A study of localisation<br />
in dual phase high-strength steels under dynamic loading using digital image correlation<br />
and FE analysis, Int. J. Solids Struct. 45 [2] (2008) 601-619. IJSS5.pdf<br />
Tensile tests were conducted on dual-phase high-strength steel in a Split-Hopkinson Tension<br />
Bar at a strain-rate in the range of 150600/s and in a servo-hydraulic testing machine at a<br />
strain-rate between 10 −3 and 100/s. A novel specimen design was utilized for the Hopkinson bar<br />
tests of this sheet material. Digital image correlation was used together with high-speed photography<br />
to study strain localisation in the tensile specimens at high rates of strain. By using<br />
digital image correlation, it is possible to obtain in-plane displacement and strain fields during<br />
non-uniform deformation of the gauge section, and accordingly the strains associated with diffuse<br />
and localized necking may be determined. The full-field measurements in high strain-rate<br />
tests reveal that strain localisation started even before the maximum load was attained in the<br />
specimen. An elasto-viscoplastic constitutive model is used to predict the observed stressstrain<br />
behavior and strain localisation for the dual-phase steel. Numerical simulations of dynamic<br />
tensile tests were performed using the non-linear explicit FE code LS-DYNA. Simulations were<br />
done with shell (plane stress) and brick elements. Good correlation between experiments and<br />
numerical predictions was achieved, in terms of engineering stressstrain behavior, deformed geometry<br />
and strain fields. However, mesh density plays a role in the localisation of deformation<br />
in numerical simulations, particularly for the shell element analysis.<br />
42
6 Beyond <strong>CORRELI</strong> Q4<br />
In the following, different references are given in which extensions of the software were used.<br />
The abstract is given for the reader to select the relevant ones.<br />
6.1 XQ4-DIC<br />
J. Réthoré, F. Hild and S. Roux, Shear-band capturing using a multiscale extended digital<br />
image correlation technique, Comp. Meth. Appl. Mech. Eng. 196 [49-52] (2007) 5016-5030.<br />
CMAME2.pdf<br />
Finite elements have been used recently to solve the optical flow conservation principle<br />
invoked to determine displacement fields by digital image correlation. Inspired by these recent<br />
advances, and by the computational effort that has been accomplished during the past 10<br />
years for the simulation of discontinuities by the extended finite element method (XFEM), an<br />
extended correlation technique is introduced for capturing shear-band like discontinuities from<br />
images of real mechanical tests. Because of the specific information at hand (i.e., gray level<br />
images), this extended finite element approach is included in a non-linear multi-grid solver. The<br />
performances of the proposed approach are discussed on two examples, namely the analysis of<br />
a bolted assembly and the formation of Piobert-Lüders bands.<br />
J. Réthoré, G. Besnard, G. Vivier, F. Hild and S. Roux, Experimental investigation of localized<br />
phenomena using Digital Image Correlation, Phil. Mag. [submitted] (2008). PM1.pdf<br />
The present paper is dedicated to the use of enriched discretization schemes in the context<br />
of digital image correlation. The aim is to capture and evaluate strong or weak discontinuities<br />
of a displacement field directly from digital images. An analysis of different enrichment performances<br />
is provided. Two examples of strain localization illustrate the discussion, namely,<br />
a compression test on an HMX-based material and the Portevin-Le Châtelier band shown in<br />
Section 4.<br />
J. Réthoré, F. Hild and S. Roux, Extended digital image correlation with crack shape optimization,<br />
Int. J. Num. Meth. Eng. 73 [2] (2008) 248-272. IJNME2.pdf<br />
The methodology of eXtended finite element method is applied to the measurement of displacements<br />
through digital image correlation. An algorithm, initially based on a finite element<br />
decomposition of displacement fields, is extended to benefit from discontinuity and singular enrichments<br />
over a suited subset of elements. This allows one to measure irregular displacements<br />
encountered, say, in cracked solids, as demonstrated both in artificial examples and experi-<br />
43
mental case studies. Moreover, an optimization strategy for the support of the discontinuity<br />
enables one to adjust the crack path configuration to reduce the residual mismatch, and hence<br />
to be tailored automatically to a wavy or irregular crack path.<br />
6.2 Beam-DIC<br />
F. Hild, S. Roux, R. Gras, N. Guerrero, M. E. Marante and J. Flórez-López, Displacement<br />
Measurement Technique for Euler-Bernoulli Kinematics, Opt. Lasers Eng. [accepted] (2008).<br />
OLE1.pdf<br />
It is proposed to develop a digital image correlation procedure that is suitable for beams<br />
whose kinematics is described by an Euler-Bernoulli hypothesis. As a direct output, the degrees<br />
of freedom corresponding to flexural and axial loads are directly measured. The performance<br />
of the correlation algorithm is evaluated by using a picture of a cantilever beam experiment.<br />
One load level is analyzed with the present algorithm. The latter is validated by comparing the<br />
displacement field with that given by a finite element based correlation algorithm. It is also<br />
shown that a locally buckled zone is detectable with the present procedure.<br />
6.3 Stereo-correlation<br />
G. Besnard, F. Hild, J.-M. Lagrange, S. Roux and C. Voltz, Contact-free characterization of<br />
materials used in detonics experiments, Proceedings Europyro 2007 . EUROPYRO_07.pdf<br />
In detonics experiments, specimens are loaded through intense shock waves produced by<br />
explosives, and hence the simultaneous occurrence of large strains and high strain rates makes<br />
the analysis delicate. Attention is here focused on tantalum whose elasto-plastic properties are<br />
quite complex, possibly involving strain localization. In order to face this challenge, the choice<br />
has been made to use imaging techniques, both for quasi-static tests (few images per second<br />
at most) and dynamic experiments (up to one million frames per second). Displacement fields<br />
are analyzed through digital image correlation techniques applied to monovision for in-plane<br />
displacements, and stereocorrelation when the three dimensional components of the surface<br />
displacement field are sought. Quasi-static tests on a tantalum specimen give access to the<br />
elastic properties (Young’s modulus and Poisson’s ratio), but also reveal strain localization<br />
occurring in the plastic flow regime. Extension of this technique to stereocorrelation is underway<br />
in order to have access to microsecond time resolution, and sub-millimeter space resolution.<br />
44
6.4 C8-DIC<br />
S. Roux, F. Hild, P. Viot and D. Bernard, Three dimensional image correlation from X-Ray<br />
computed tomography of solid foam, Comp. Part A DOI: 10.1016/j.compositesa.2007.11.011<br />
[in press] (2008). CompPartA1.pdf<br />
A new methodology is proposed to estimate 3D displacement fields from pairs of images<br />
obtained from X-Ray Computed Micro Tomography (XCMT). Contrary to local approaches,<br />
a global approach is followed herein that evaluates continuous displacement fields. Although<br />
any displacement basis may be considered, the procedure is specialized to finite element shape<br />
functions. The method is illustrated with the analysis of a compression test on a polypropylene<br />
solid foam (independently studied in a companion paper). A good stability of the measured<br />
displacement field is obtained for cubic element sizes ranging from 16 voxels to 6 voxels. This<br />
last paper is a direct extension of Q4-DIC to three-dimensional displacement measurements in<br />
the bulk of a material.<br />
7 Conclusion<br />
A novel approach was developed to determine a displacement field based on the comparison<br />
of two digital images. The sought displacement field is decomposed onto a basis of continuous<br />
functions using Q4P1-shape functions as proposed in finite element methods. The latter corresponds<br />
to one of the simplest kinematic descriptions. It therefore allows for a compatibility of<br />
the kinematic hypotheses made during the measurement stage and the subsequent identification<br />
/ validation stage, for instance, by using conventional finite element techniques.<br />
The performance of the algorithm is tested on a reference image to evaluate the reliability<br />
of the estimation, which is shown to allow for either an excellent accuracy for homogeneous<br />
displacement fields, or for a very well resolved displacement field down to element sizes as small<br />
as 4×4 pixels. When compared with a standard FFT-DIC technique, Q4-DIC enables for a<br />
significant decrease of the displacement uncertainty when the same size of the elements is used.<br />
For the natural texture studied herein, it is shown that the main limitation is related to its<br />
poor dynamic range that makes it very sensitive to acquisition noise.<br />
Last, the displacement field is analyzed on an aluminium alloy sample where a localization<br />
band develops. In the present case, it corresponds to a Portevin-Le Châtelier effect. For element<br />
sizes ranging from 16 down to 4 pixels, the displacement field is shown to be reliably determined.<br />
As far as strains are concerned, spurious strain concentrations at the edge of the shear band<br />
is observed for small element sizes. This is attributed to integer locking of the displacement.<br />
Apart from this bias, the analysis has been shown to be operational and reliable, on a real<br />
45
experimental test with strain localization. The complete space-time kinematic analysis is now<br />
possible for the Portevin-Le Châtelier phenomenon even within the band thanks to the spatial<br />
resolution achieved by the Q4-DIC technique presented herein. It is deferred to a another<br />
publication in which an enriched kinematics is implemented [31].<br />
References<br />
[1] P. K. Rastogi, edts., Photomechanics, (Springer, Berlin (Germany), 2000), 77.<br />
[2] M. A. Sutton, S. R. McNeill, J. D. Helm and Y. J. Chao, Advances in Two-Dimensional and<br />
Three-Dimensional Computer Vision, in: Photomechanics, P. K. Rastogi, eds., (Springer,<br />
Berlin (Germany), 2000), 323-372.<br />
[3] A. Lagarde, edt., Advanced Optical Methods and Applications in Solid Mechanics,<br />
(Kluwer, Dordrecht (the Netherlands), 2000), 82.<br />
[4] O. C. Zienkievicz and R. L. Taylor, The Finite Element Method, (McGraw-Hill, London<br />
(UK), 4th edition, 1989).<br />
[5] B. K. P. Horn and B. G. Schunck, Determining optical flow, Artificial Intelligence 17<br />
(1981) 185-20.<br />
[6] P. J. Hubert, Robust Statistics, (Wiley, New York (USA), 1981).<br />
[7] M. Black, Robust Incremental Optical Flow, (PhD dissertation, Yale University, 1992).<br />
[8] J.-M. Odobez and P. Bouthemy, Robust multiresolution estimation of parametric motion<br />
models, J. Visual Comm. Image Repres. 6 (1995) 348-365.<br />
[9] M. A. Sutton, W. J. Wolters, W. H. Peters, W. F. Ranson and S. R. McNeill, Determination<br />
of Displacements Using an Improved Digital Correlation Method, Im. Vis. Comp. 1 [3]<br />
(1983) 133-139.<br />
[10] F. Hild, B. Raka, M. Baudequin, S. Roux and F. Cantelaube, Multi-Scale Displacement<br />
Field Measurements of Compressed Mineral Wool Samples by Digital Image Correlation,<br />
Appl. Optics IP 41 [32] (2002) 6815-6828.<br />
[11] S. Roux and F. Hild, Stress intensity factor measurements from digital image correlation:<br />
post-processing and integrated approaches, Int. J. Fract. 140 [1-4] (2006) 141-157.<br />
46
[12] F. Hild and S. Roux, Digital image correlation: from measurement to identification of<br />
elastic properties - A review, Strain 42 (2006) 69-80.<br />
[13] F. Hild, S. Roux, R. Gras, N. Guerrero, M. E. Marante and J. Flórez-López, Displacement<br />
Measurement Technique for Euler-Bernoulli Kinematics, Opt. Lasers Eng. [accepted]<br />
(2008).<br />
[14] B. Wagne, S. Roux and F. Hild, Spectral Approach to Displacement Evaluation From<br />
Image Analysis, Eur. Phys. J. AP 17 (2002) 247-252.<br />
[15] S. Roux, F. Hild and Y. Berthaud, Correlation Image Velocimetry: A Spectral Approach,<br />
Appl. Optics 41 [1] (2002) 108-115.<br />
[16] P. J. Burt, C. Yen and X. Xu, Local correlation measures for motion analysis: a comparative<br />
study, Proceedings IEEE Conf. on Pattern Recognition and Image Processing, (1982),<br />
269-274.<br />
[17] H. W. Schreier, J. R. Braasch and M. A. Sutton, Systematic errors in digital image correlation<br />
caused by intensity interpolation, Opt. Eng. 39 [11] (2000) 2915-2921.<br />
[18] S. Bergonnier, F. Hild and S. Roux, Local anisotropy analysis for non-smooth images,<br />
Patt. Recogn. 40 [2] (2007) 544-556.<br />
[19] J.-N. Périé, S. Calloch, C. Cluzel and F. Hild, Analysis of a Multiaxial Test on a C/C<br />
Composite by Using Digital Image Correlation and a Damage Model, Exp. Mech. 42 [3]<br />
(2002) 318-328.<br />
[20] G. Holst, CCD Arrays, Cameras and Displays, (SPIE Engineering Press, Washington DC<br />
(USA), 1998).<br />
[21] ISO, International Vocabulary of Basic and General Terms in Metrology (VIM), (International<br />
Organization for Standardization, Geneva (Switzerland), 1993).<br />
[22] G. Besnard, F. Hild and S. Roux, “Finite-element” displacement fields analysis from digital<br />
images: Application to Portevin-Le Châtelier bands, Exp. Mech. 46 (2006) 789-803.<br />
[23] G. Besnard, Corrélation d’images et identification en mécanique des solides, (University<br />
Paris 6, Intermediate MSc report 2005).<br />
[24] A. Portevin and F. Le Châtelier, C. R. Acad. Sci. Paris 176 (1923) 507.<br />
[25] A. H. Cottrell, A note on the Portevin-Le Châtelier effect, Phil. Mag. 44 (1953) 829-832.<br />
47
[26] P. G. McCormick, A model of the Portevin-Le Châtelier effect in substitutional alloys,<br />
Acta Metall. 20 (1972) 351-354.<br />
[27] A. H. Clausen, T. Borvik, O. S. Hopperstad and A. Benallal, Flow and fracture characteristics<br />
of aluminium alloy AA5083-H116 as function of strain rate, temperature and<br />
triaxiality, Mat. Sci. Eng. A364 (2004) 260-272.<br />
[28] D. Thevenet, M. Mliha-Touati and A. Zeghloul, Characteristics of the propagating deformation<br />
bands associated with the Portevin-Le Châtelier effect in an Al-Zn-Mg-Cu alloy,<br />
Mat. Sci. Eng. A291 (2000) 110-117.<br />
[29] R. Shabadi, S. Kumara, H. J. Roven and E. S. Dwarakadasa, Characterisation of PLC<br />
band parameters using laser speckle technique, Mat. Sci. Eng. A364 (2004) 140-150.<br />
[30] H. Louche, P. Vacher and R. Arrieux, Thermal observations associated with the Portevin-<br />
Le Châtelier effect in an Al-Mg alloy, Mat. Sci. Eng. A 404 (2005) 188-196.<br />
[31] J. Réthoré, G. Besnard, G. Vivier, F. Hild and S. Roux, Experimental investigation of<br />
localized phenomena using Digital Image Correlation, Phil. Mag. [submitted] (2008).<br />
[32] R. Hill, Aspects of Invariance in Solid Mechanics, Adv. Appl. Mech. 18 (1978) 1-75.<br />
48
Appendix: Strain measures<br />
From the displacement field, the strain field is usually sought. Different strain measures can be<br />
used to analyze a mechanical experiment.<br />
7.1 Infinitesimal strains<br />
A strain is defined in terms of relative variation of length of a line linking two material points<br />
of a deformable body. This intuitive definition will be completed in the following subsections.<br />
The displacement of a point P (x, y, z) is defined by three components (u, v, w) of a vector u in<br />
a cartesian coordinate system chosen in the undeformed state<br />
A second point Q in the vicinity of P moves to Q ′<br />
P (x, y, z) → P ′ (x + u, y + v, z + w) (16)<br />
Q(x + dx, y + dy, z + dz) →<br />
Q ′ (x + dx + u + du, y + dy + v + dv, z + dz + w + dw) (17)<br />
The quantities du, dv, dw are called relative displacements. When du, dv, dw are sufficiently<br />
small (i.e., infinitesimal) compared to the length scale dx, dy, dz, a first order approximation<br />
of the variations du, dv, dw can be written as<br />
du = ∂u<br />
∂x<br />
dv = ∂v<br />
∂x<br />
dw = ∂w<br />
∂x<br />
∂u ∂u<br />
dx + dy +<br />
∂y ∂z dz<br />
dx +<br />
∂v<br />
∂y<br />
The previous equation can be written in matrix form<br />
∂v<br />
dy + dz (18)<br />
∂z<br />
∂w ∂w<br />
dx + dy +<br />
∂y ∂z dz.<br />
du = gradu dx (19)<br />
Of the nine components obtained in Eqn. (18), they are grouped as follows for the normal<br />
strains<br />
for the shear strains<br />
ε xx = ∂u<br />
∂x , ε yy = ∂v<br />
∂y<br />
ε xy = ε yx = 1 ( ∂u<br />
2 ∂y + ∂v )<br />
∂x<br />
ε yz = ε zy = 1 ( ∂v<br />
2 ∂z + ∂w )<br />
∂y<br />
ε zx = ε xz = 1 ( ∂w<br />
2 ∂x + ∂u )<br />
∂z<br />
49<br />
, ε zz = ∂w<br />
∂z<br />
(20)<br />
(21)
and the infinitesimal rotations<br />
ω z = 1 (<br />
− ∂u<br />
2 ∂y + ∂v )<br />
∂x<br />
ω x = 1 (<br />
− ∂v<br />
2 ∂z + ∂w )<br />
∂y<br />
ω y = 1 (<br />
− ∂w<br />
2 ∂x + ∂u )<br />
∂z<br />
These results are related to the symmetric and antisymmetric decomposition of the gradient<br />
operator<br />
(22)<br />
gradu = grad s u + grad a u (23)<br />
where the infinitesimal strain tensor ε is defined by<br />
ε = grad s u = 1 (<br />
gradu + grad t u ) (24)<br />
2<br />
and the infinitesimal rotation tensor Ω is expressed as<br />
Ω = grad a u = 1 (<br />
gradu − grad t u ) (25)<br />
2<br />
It can be shown that an antisymmetric tensor can be related to a vector (here the rotation<br />
vector ω) by<br />
Examples:<br />
Ωdx = ω × dx (26)<br />
Let us first consider a rigid translation. The two points P and Q then move by an identical<br />
amount u<br />
P → P ′ (x + u, y + v, z + w)<br />
Q → Q ′ (x + dx + u, y + dy + v, z + dz + w) (27)<br />
so that the displacement gradient vanishes, so do the strains and the rotation vector.<br />
The second example deals with a small rigid rotation. The two points P and Q then move<br />
by different amounts u<br />
P → P ′ (x + u, y + v, z + w)<br />
Q → Q ′ (x + dx + u + du, y + dy + v + dv, z + dz + w + dw) (28)<br />
with<br />
du = −ω z dy + ω y dz<br />
dv = −ω x dz + ω z dx<br />
dw = −ω y dx + ω x dy (29)<br />
50
so that the displacement gradient has only an antisymmetric part (i.e., ω x , ω y , ω z are the components<br />
of the rotation vector ω), and the strains are vanishing.<br />
The third case is such that<br />
du = adx<br />
dv = bdy<br />
dw = cdz (30)<br />
The displacement gradient is symmetric (i.e., no rotation occurs). The strain tensor has non<br />
vanishing components associated to normal strains<br />
ε xx = a , ε yy = b , ε zz = c. (31)<br />
7.2 Finite transformation strains<br />
The aim of this part is to introduce the finite strains as an extension of the previous approach.<br />
The more difficult aspects will be discussed in the next part. Let us consider again a point<br />
P (x, y, z) experiencing a displacement u<br />
P → P ′ (x + u, y + v, z + w) (32)<br />
and a second point Q(x + dx, y + dy, z + dz) in the vicinity of P moves to Q ′<br />
Q → Q ′ (x + dx + u + du, y + dy + v + dv, z + dz + w + dw). (33)<br />
The theory of finite strain considers the length variation of a line joining the two considered<br />
points in a reference state (dl = P Q) and in a deformed state (dl ′ = P ′ Q ′ )<br />
(dl) 2 = dx 2 + dy 2 + dz 2<br />
(dl ′ ) 2 = (dx + du) 2 + (dy + dv) 2 + (dz + dw) 2 . (34)<br />
As a first approximation, Eqn. (18) is still used, and the length ratio becomes<br />
By noting that<br />
( ) dl<br />
′ 2<br />
=<br />
dl<br />
( dx<br />
dl + du<br />
dl<br />
du<br />
dl<br />
dv<br />
dl<br />
dw<br />
dl<br />
) 2<br />
+<br />
( dy<br />
dl + dv ) 2<br />
+<br />
dl<br />
= ∂u dx<br />
∂x dl + ∂u dy<br />
∂y dl + ∂u dz<br />
∂z dl<br />
= ∂v dx<br />
∂x dl + ∂v dy<br />
∂y dl + ∂v dz<br />
∂z dl<br />
= ∂w dx<br />
∂x dl + ∂w dy<br />
∂y dl + ∂w dz<br />
∂z dl ,<br />
( dz<br />
dl + dw<br />
dl<br />
) 2<br />
. (35)<br />
(36)<br />
51
Eqn. (35) can be rewritten as<br />
( ) dl<br />
′ 2<br />
=<br />
dl<br />
+<br />
+<br />
( dx<br />
dl + ∂u dx<br />
∂x dl + ∂u dy<br />
∂y dl + ∂u ) 2<br />
dz<br />
∂z dl<br />
( dy<br />
dl + ∂v dx<br />
∂x dl + ∂v dy<br />
∂y dl + ∂v ) 2<br />
dz<br />
(37)<br />
∂z dl<br />
( dz<br />
dl + ∂w dx<br />
∂x dl + ∂w dy<br />
∂y dl + ∂w ) 2<br />
dz<br />
.<br />
∂z dl<br />
The strains are expressed with respect to a reference state that is undeformed. Let a, b and c<br />
denote the cosine directors of dl<br />
so that, by definition<br />
a = dx<br />
dl<br />
, b = dy<br />
dl<br />
, c = dz<br />
dl , (38)<br />
( ) dl<br />
′ 2<br />
= (1 + 2E xx )a 2 + (1 + 2E yy )b 2 + (1 + 2E zz )c 2<br />
dl<br />
+ 4E xy ab + 4E yz bc + 4E zx ca, (39)<br />
with the normal strains<br />
E xx = ∂u<br />
∂x + 1 2<br />
E yy = ∂v<br />
∂y + 1 2<br />
E zz = ∂w<br />
∂z + 1 2<br />
[ (∂u ) 2 ( ) 2 ( ) ] 2 ∂u ∂u<br />
+ +<br />
∂x ∂y ∂z<br />
[ (∂v ) 2 ( ) 2 ( ) ] 2 ∂v ∂v<br />
+ +<br />
∂x ∂y ∂z<br />
[ (∂w ) 2 ( ) 2 ( ) ] 2 ∂w ∂w<br />
+ +<br />
∂x ∂y ∂z<br />
and the shear strains<br />
E xy = e yx = 1 ( ∂u<br />
2 ∂y + ∂v<br />
∂x + ∂u ∂u<br />
∂x ∂y + ∂v ∂v<br />
∂x ∂y + ∂w )<br />
∂w<br />
∂x ∂y<br />
E yz = e zy = 1 ( ∂v<br />
2 ∂z + ∂w<br />
∂y + ∂u ∂u<br />
∂y ∂z + ∂v ∂v<br />
∂y ∂z + ∂w )<br />
∂w<br />
∂y ∂z<br />
E zx = e xz = 1 ( ∂w<br />
2 ∂x + ∂u<br />
∂z + ∂u ∂u<br />
∂z ∂x + ∂v ∂v<br />
∂z ∂x + ∂w )<br />
∂w<br />
∂z ∂x<br />
(40)<br />
(41)<br />
corresponding to the Green-Lagrange strain tensor E<br />
E = 1 2<br />
(<br />
gradu + grad t u + grad t u gradu ) (42)<br />
52
The normal strains can be interpreted as<br />
[ (dx<br />
E xx = 1 )<br />
′ 2<br />
− 1]<br />
= 1 [<br />
(1 + µx ) 2 − 1 ]<br />
2 dx<br />
2<br />
[ (dy<br />
E yy = 1 )<br />
′ 2<br />
− 1]<br />
= 1 [<br />
(1 + µy ) 2 − 1 ] (43)<br />
2 dy<br />
2<br />
[ (dz<br />
E zz = 1 )<br />
′ 2<br />
− 1]<br />
= 1 [<br />
(1 + µx ) 2 − 1 ] ,<br />
2 dz<br />
2<br />
where µ x , µ y , µ z are the relative length variations (dx ′ − dx)/dx, (dy ′ − dy)/dy, (dz ′ − dz)/dz,<br />
respectively. Similarly, the shear strains are related to angle variations<br />
2E xy = √ √<br />
1 + 2E xx 1 + 2Eyy cos θ xy<br />
2E yz = √ √<br />
1 + 2E yy 1 + 2Ezz cos θ yz (44)<br />
2E zx = √ √<br />
1 + 2E zz 1 + 2Exx cos θ zx ,<br />
where θ xy , θ yz , θ zx are the distortions of the vectors dy ′ wrt. dx ′ , dz ′ wrt. dy ′ , dx ′ wrt. dz ′ ,<br />
respectively. (It is worth remembering that dx.dy = 0, dy.dz = 0, and dz.dx = 0.)<br />
By noting that<br />
−1 ≤ µ x , µ y , µ z ≤ +∞ , −π/2 ≤ θ xy , θ yz , θ zx ≤ π/2 (45)<br />
the following bounds can be obtained for the components of the Green-Lagrange strain components<br />
−1/2 ≤ E xx , E yy , E zz ≤ +∞ , −1/2 ≤ E xy , E yz , E zx ≤ 1/2. (46)<br />
7.3 Large transformation kinematics<br />
Let us consider two neighboring points P and Q of a body Ω in its reference state, and a<br />
corresponding frame R 0 . If a rigid body motion is applied, it corresponds to an isometry,<br />
namely the distance between the points P ′ and Q ′ after the transformation characterized by a<br />
tensor A<br />
P ′ Q ′ = A P Q (47)<br />
is identical to the initial one<br />
‖P ′ Q ′ ‖ 2 = P Q A t A P Q = ‖P Q‖ 2 , (48)<br />
so that<br />
A t A = I 2 (49)<br />
53
where I 2 is the second order unit tensor.<br />
expressed as<br />
or equivalently<br />
where the tensor<br />
dP ′ Q ′<br />
dt<br />
The first order derivative of the vector P ′ Q ′ is<br />
= V (Q ′ /R 0 ) − V (P ′ /R 0 ) (50)<br />
dP ′ Q ′<br />
= A<br />
dt<br />
˙ A t P ′ Q ′ (51)<br />
AA ˙ t is antisymmetric. As shown above, an antisymmetric tensor acting on<br />
a vector can also be written as the cross product of a vector, here the rotation vector ω of Ω<br />
with respect to R 0 , with P ′ Q ′ dP ′ Q ′<br />
= ω(Ω/R 0 ) × P ′ Q ′ . (52)<br />
dt<br />
Consequently, the vector field associated to a rigid body kinematics is such that<br />
V (Q ′ /R 0 ) = V (P ′ /R 0 ) + ω(Ω/R 0 ) × P ′ Q ′ (53)<br />
The motion of a body Ω is described by the function x ′ = x ′ (x, t) giving the position<br />
x ′ at a time t of a particle that occupied the position x prior to a deformation. At a fixed<br />
time, this function defines the deformation of the body between its reference configuration C<br />
and the current one C ′ . In Solid Mechanics, one chooses as reference configuration, the initial<br />
configuration C and the current configuration in the same frame. It is however preferable,<br />
to avoid any confusion, to distinguish these two frames, and in particular the Lagrangian<br />
coordinates x for C and the Eulerian coordinates x ′ for C ′ . For the sake of simplicity, the two<br />
systems are chosen orthogonal and normalized. As performed above, it is common to introduce<br />
the displacement vector u(x, t)<br />
x ′ (x, t) = x + u(x, t) (54)<br />
Contrary to infinitesimal motions, no particular hypotheses are made concerning u. The function<br />
x ′ (x, t) defines the global motion of the body. Locally, i.e., for two points P and Q<br />
separated by dx, one uses the deformation gradient tensor F relating an infinitesimal vector<br />
dx in the reference configuration (C) to dx ′ in the deformed configuration (C ′ )<br />
dx ′ = F dx (55)<br />
so that the tensor F can be related to the displacement gradient ∇u by<br />
F = I 2 + ∇u. (56)<br />
The deformation gradient tensor defines the local motion. To measure a strain, i.e., a shape<br />
change, it is necessary to eliminate the rotation. To characterize shape changes, one needs to<br />
54
determine length variations as well as angle variations. In both cases, the information is given<br />
by scalar products. Let us form the scalar product of dx ′ with itself<br />
dx ′ .dx ′ = dx.F t F .dx (57)<br />
where C = F t F denotes the right Cauchy-Green strain tensor. The Lagrangian strain measures<br />
are formed by using C. It is possible to define various strain measures. The use of one rather<br />
than another is a matter of choice or convenience. The only requirements are:<br />
• they vanish for any rigid body motion (i.e., C = I 2 for any rigid body motion),<br />
• usually, they are positive for an expansion and negative for a contraction.<br />
For Lagrangian measures, they can be expressed by using the strain tensors E m [32]<br />
⎧<br />
⎨<br />
E m =<br />
⎩<br />
1<br />
2m (Cm − 1) when m ≠ 0<br />
1<br />
2 when m → 0+ (58)<br />
When m = 1, the Green-Lagrange strain tensor is obtained (see Section 7.2), m = 1/2 is the<br />
Cauchy-Biot (or nominal) strain tensor and yields ∆L/L 0 for a uniaxial elongation, where L 0<br />
is a reference (or gauge) length and ∆L the length variation. The case m → 0 + corresponds to<br />
the logarithmic (or Hencky) strain tensor. The latter is the only additive strain measure. By<br />
definition, all these strains are equal to zero for a rigid translation (i.e., F = I 2 and C = I 2 )<br />
or a rigid rotation (i.e., F = R and C = I 2 ), where R is an orthogonal tensor. When the<br />
amplitude of the body motion is small as well as the strain gradient, all measures converge<br />
towards the infinitesimal strain tensor ε defined by<br />
ε = 1 2<br />
[<br />
∇u + (∇u)<br />
t ] . (59)<br />
55