Computer-implemented method for identifying mechanical properties by coupled correlation of images and mechanical modelling

Patent metadata

Field	Value
Publication number	US-11215446-B2
Application number	US-201816976425-A
Country	US
Kind code	B2
Filing date	Feb 28, 2018
Priority date	Feb 28, 2018
Publication date	Jan 4, 2022
Grant date	Jan 4, 2022

Abstract

A computer-implemented method for identifying mechanical parameters of an object subjected to mechanical stress is provided. The method comprises a step of acquiring, by an imaging means, images of the object taken before and during the application of the mechanical stress, three steps of calculating the effects due to the stress carried out either on the basis of the modeling of the recorded images or on the basis of a theoretical mechanical modeling of the stress, a step of defining a functional equal to the difference between the two models and a last step of minimizing said functional so that the experimental model is as close as possible to the theoretical mechanical model. Additional measurements make it possible to refine the method.

First claim

The invention claimed is:

1. A computer-implemented method for identifying at least one mechanical parameter called a “target parameter” of an object subjected to a mechanical stress, wherein said computer-implemented method comprises the following steps:
Step 1: acquiring, by an imaging means, at least two digital images of the object taken before and during an application of the mechanical stress and measuring a scale factor of the object;
Step 2: calculating with a computer a first functional T CIN (U CIN ) corresponding to a correlation of the at least two digital images depending on a displacement field U CIN represented using a first kinematic base, said displacement field U CIN being measured at any point of the object under stress between the at least two digital images of the object under load and without load;
Step 3: calculating with the computer a calculated displacement field U CAL at any point of the object;
Step 4: calculating with the computer a second functional T CAL (U CAL ,{p},{q}) on a basis of the calculated displacement field U CAL represented using a second kinematic base, this second functional corresponding to a variational formulation of a mechanical model of the stress depending on a geometry of the object, forces applied, boundary conditions, at least a target parameter {p} and predetermined mechanical parameters {q};
Step 5: calculating with the computer a third functional T PAR (U CIN ,U CAL ) in a form of a quadratic norm, equal to a difference between U CIN and U CAL ; and
Step 6: minimizing with the computer, with respect to U CIN , U CAL and {p}, a total functional T TOT (U CIN ,U CAL ,{p},{q}) comprising at least the terms:
T TOT ( U CIN ,U CAL ,{p},{q })=α T CIN ( U CIN )+β T CAL ( U CAL ,{p},{q })γ T PAR ( U CIN ,U CAL )

α,β and γ being three non-zero weighting coefficient,
wherein the step of minimizing with the computer identifies the “target parameter” that relates to at least one of the following: a mechanical property of a constitutive material of the object, one or more phases of the object, a geometry of the object, and/or quantities of the object; and
wherein the imaging means comprises at least one of the following: a video camera, a still camera, a scanning electron microscope, an atomic force microscope, a tomography apparatus, an X-ray tomography apparatus, a magnetic resonance tomography apparatus, and/or an optical coherence tomography apparatus.

2. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein when a behavior of the object is subjected to a time-dependent stress, the second functional T CAL (U CAL , {p}, {q}) is dependent on determined times.

3. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein step 1 of the method comprises additional measurements F MES of forces, times or temperatures, step 3 of the method comprises evaluations F CAL corresponding to said additional measurements, step 4 is followed by a step 4bis of calculating a fourth functional T FOR (F CAL ,F MES ) proportional to a quadratic deviation between these quantities and a total functional T TOT (U CIN ,U CAL ,{p},{q}) of step 5 is equal to:
T TOT ( U CIN ,U CAL ,{p},{q })=α T CIN ( U CIN )+β T CAL ( U CAL ,{p},{q })+γ T PAR ( U CIN ,U CAL )+χT FOR (F CAL ({p},{q}), F MES )

χ being a fourth weighting coefficient, said weighting coefficient χ is adjusted according to uncertainties associated with various quantities involved in the functionals, and/or according to a condition number of a problem tangent to the minimization of the functional T TOT .

4. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein the minimization of the total functional T TOT is carried out by iterative method.

5. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein the minimization of the total functional T TOT is carried out by iterative method, requiring the calculation of a gradient of T TOT .

6. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein the first kinematic base is identical to the second kinematic base.

7. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein a measurement uncertainty is estimated by a Hessian of the functional T TOT with convergence by having a measurement of acquisition noise available.

8. The computer-implemented method for identifying at least one mechanical parameter as claimed in claim 1 , wherein the first kinematic base or the second kinematic base is produced on a finite element mesh.

9. A computer device suitable for identifying at least one mechanical parameter as claimed in claim 1 .

10. A computer-readable medium with program for executing the method as claimed in claim 1 .

11. A mechanical test system configured to identify at least one mechanical parameter called a “target parameter” of an object subjected to a mechanical stress, wherein said the mechanical test system comprises:
an imaging device configured to acquire at least two digital images of the object taken before and during an application of the mechanical stress and measuring a scale factor of the object;
a computer configured to calculate a first functional T CIN (U CIN ) corresponding to a correlation of the at least two digital images depending on a displacement field U CIN represented using a first kinematic base, said displacement field U CIN being measured at any point of the object under stress between the at least two digital images of the object under load and without load;
the computer further configured to calculate a calculated displacement field U CAL at any point of the object;
the computer further configured to calculate a second functional T CAL (U CAL ,{p},{q}) on a basis of the calculated displacement field U CAL represented using a second kinematic base, this second functional corresponding to a variational formulation of a mechanical model of the stress depending on a geometry of the object, forces applied, boundary conditions, at least a target parameter {p} and predetermined mechanical parameters {q};
the computer further configured to calculate a third functional T PAR (U CIN ,U CAL ) in a form of a quadratic norm, equal to a difference between U CIN and U CAL ; and
the computer further configured to calculate a minimization of, with respect to U CIN , U CAL and {p}, a total functional T TOT (U CIN ,U CAL ,{p},{q}) comprising at least the terms:
T TOT ( U CIN ,U CAL ,{p},{q })=α T CIN ( U CIN )β T CAL ( U CAL ,{p},{q })γ T PAR ( U CIN ,U CAL )

α,β and γ being three non-zero weighting coefficients,
wherein the step of minimizing with the computer identifies the “target parameter” that relates to at least one of the following: a mechanical property of a constitutive material of the object, one or more phases of the object, a geometry of the object, and/or quantities of the object; and
wherein the imaging device comprises at least one of the following: a video camera, a still camera, a scanning electron microscope, an atomic force microscope, a tomography apparatus, an X-ray tomography apparatus, a magnetic resonance tomography apparatus, and/or an optical coherence tomography apparatus.

12. A non-transitory computer-readable medium having instructions for execution by a computer for identifying at least one mechanical parameter called a “target param

Assignees

• Centre Nat Rech Scient
• Ecole Normale Superieure Cachan

Inventors

• Hild François
• Leclerc Hugo
• Roux Stéphane

Classifications

• G06T7/32 Primary

  using correlation-based methods · CPC title

• G01M5/005

  by means of external apparatus, e.g. test benches or portable test systems · CPC title

• G01N3/08

  by applying steady tensile or compressive forces (G01N3/28 takes precedence) · CPC title

• G06F2111/10

  Numerical modelling · CPC title

• G06T2207/10012

  Stereo images · CPC title