550 likes | 669 Vues
Corrélation d'images numériques: Stratégies de régularisation et enjeux d'identification. Stéphane Roux, François Hild LMT, ENS-Cachan. Atelier « Problèmes Inverses », Nancy, 7 Juin 2011. Image 2. Image 1. Relative displacement field ?. Image 2. Image 1. Deformed image. Reference image.
E N D
Corrélation d'images numériques: Stratégies de régularisation et enjeux d'identification Stéphane Roux, François Hild LMT, ENS-Cachan Atelier « Problèmes Inverses », Nancy, 7 Juin 2011
Image 2 Image 1 Relative displacement field ?
Image 2 Image 1
Deformed image Reference image Relative displacement field ?
Displacement field Uy Deformed image Reference image
Displacement fields are nice, but … Can we get more ?
Stress intensity Factor, Crack geometry Image 2 Image 1
Damage field Deformed image Reference image
Constitutive law Deformed image Reference image
Outline • A brief introduction to “global DIC” • Mechanical identification • Regularization
From texture to displacements DIC in a nutshell
Digital Image Correlation • Images (gray levels) indexed by time t • Texture conservation (passive tracers) (hypothesisthatcanberelaxed if needed)
Problem to solve • Weak formulation: Minimize wrt u where the residual is Provides a spatiallyresolvedqualityfield of the proposed solution
Solution • The problemisintrinsicallyill-posed and highly non-linear ! • A specificstrategy has to bedesigned for accurate and robust convergence • It impacts on the choice of the kinematic basis
Global DIC • Decompose the soughtdisplacementfield on a suited basis providing a naturalregularization • Yn: • FEM shapefunction, X-FEM, … • Elastic solutions, Numericallycomputedfields, Beamkinematics…
The benefit of C0 regularization ZOI size / Element size (pixels) Key parameter = (# pixels)/(# dof)
Example: T3-DIC* Pixel size = 67 mm *[Leclerc et al., 2009, LNCS 5496 pp. 161-171]
Example: T3-DIC Ux (pixel) 0.46 0.28 0.11 -0.06 -0.23 [H. Leclerc]
Example: T3-DIC Uy (pixel) 0.54 0.35 0.15 -0.04 -0.24
Example: T3-DIC Residual 28 21 14 7 0 Mean residual = 3 % dynamic range
The real challenge • For solidmechanics application, the actual challenge is • not to get the displacementfields, but rather • to identify the constitutive law (stress/strain relation) • The simplest case islinearelasticity
Plane elasticity • A potential formulation canbeadoptedshowingthat the displacementfieldcanbewrittengenerically in the complex plane as whereF and Y are arbitraryholomorphicfunctions • mis the shearmodulus, • kis a dimensionlesselastic constant (related to Poisson’s ratio)
Plane elasticity • It suffices to introduce a basis of test functions for F(z) and Y(z) and considerthat and are independent • Direct evaluation of 1/m and k/m
Validatedexamples • Brazilian compression test • Cracks
Example 1:Brazilian compression test • Integrated approach: decomposition of the displacement field over 4 fields (rigid body motion + analytical solution)
Integrated approach Identifiedproperties for the polycarbonate m 880 MPa n 0.45 In good agreement withliterature data
Need for coupling to modelling • Elasticity (or incremental non-linearbehavior) • FEM
Dialog DIC/FEA modeling • Local elastic identification R. Gras, Comptest2011
More generalframework • Inhomogeneouselasticsolid • Non-linear constitutive law • Plasticity • Damage • Non-linearelasticity
Mechanical regularization • The displacement field should be such that or in FEM language for interior nodes. This can be used to help DIC
Integrated DIC • Reach smaller scale H. Leclerc et al., Lect. Notes Comp. Sci. 5496, 161-171, (2009)
Tikhonov type regularization • Minimization of • Regularizationisneutralwith respect to rigid body motion • How should one chooseA ?
Spectral analysis • For a test displacementfield Regularization log(||.||2) DIC Cross-over scale log(k)
Boundaries • The equilibrium gap functionalisoperativeonly for interiornodes or free boundaries • Atboundaries, information maybelacking • Introduce an additionalregularizationterm(e.g. ) • Extendelasticbehavioroutside the DIC analyzedregion
Regularizationat voxel scale • An example in 3D for a modest size 243 voxels
Voxel scale DVC 1 voxel 5.1 µm Displacementnorm (voxels) Vertical displacement (voxels) H. Leclerc et al., Exp. Mech. (2011)
Identification • As a post-processingstep, a damage lawcanbeidentifiedfrom the minimization ofwhereU has been measured and Kisknown • Manyunknowns !
Validation < 5.3 %
Constitutive law State potential (isotropic damage) State laws Dissipated power Thermodynamic consistency Growth law ~ equivalentscalarstrain
Use of a homogeneous constitutive law • Postulating a homogeneouslaw, damage is no longer a twodimensionalfield of unknowns, but a (non-linear) function of the maximum strainexperienced by an element of volume.
Damage growth law • Identifiedform or truncation
Identified damage image 11 log10(1-D)