CN112055871A - 通过耦合的图像相关性与机械建模来识别机械性能的计算机实现方法 - Google Patents
通过耦合的图像相关性与机械建模来识别机械性能的计算机实现方法 Download PDFInfo
- Publication number
- CN112055871A CN112055871A CN201880090358.4A CN201880090358A CN112055871A CN 112055871 A CN112055871 A CN 112055871A CN 201880090358 A CN201880090358 A CN 201880090358A CN 112055871 A CN112055871 A CN 112055871A
- Authority
- CN
- China
- Prior art keywords
- cal
- cin
- functional
- mechanical
- tot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000003384 imaging method Methods 0.000 claims abstract description 8
- 238000006073 displacement reaction Methods 0.000 claims description 38
- 238000005259 measurement Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 7
- ILJCBKGAZSIGIH-UHFFFAOYSA-N 2-oxoimidazolidine-1-carbothioamide Chemical compound NC(=S)N1CCNC1=O ILJCBKGAZSIGIH-UHFFFAOYSA-N 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims 1
- 238000001454 recorded image Methods 0.000 abstract 1
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000002091 elastography Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007431 microscopic evaluation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012014 optical coherence tomography Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
- G01B15/06—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring the deformation in a solid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/08—Testing mechanical properties
- G01M11/081—Testing mechanical properties by using a contact-less detection method, i.e. with a camera
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
- G01M5/005—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/001—Industrial image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/97—Determining parameters from multiple pictures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/607—Specific applications or type of materials strain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30164—Workpiece; Machine component
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Electromagnetism (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明的一般领域是识别经受机械应力的对象的机械参数的计算机实现的方法。根据本发明的方法包括:通过成像装置获取施加机械应力之前和期间拍摄的对象的图像的步骤;计算由于在所记录图像的建模基础上或在所述应力的理论机械建模基础上执行应力而带来的影响的三个步骤;定义等于两个模型之间差异的泛函的步骤;以及最小化所述泛函的最终步骤,以便实验模型尽可能接近理论机械模型。附加测量使得根据本发明的方法能够细化。
Description
本发明的一般领域是材料和结构科学领域,更确切地说是材料和结构在经受各种机械应力时的性能的知识领域。
主要的工业应用与材料的强度有关。然而,还存在诸如无损检测、某些生物学应用或医用弹性成像之类的其他应用领域。
一般地,测量是借助于包括成像设备的机械测试系统进行的。这些设备主要是在可见光波长或接近可见光波长下工作的视频摄像头或静态摄像头,或是诸如扫描电子显微镜或原子力显微镜之类的微观分析系统。断层扫描装置也可用于分析材料样品的整个体积。这种装置有各种类型。将非穷尽地提及X射线断层扫描装置、磁共振断层扫描装置或“MRI”和光学相干断层扫描装置。
这些成像设备的一般工作原理如下。例如,在应用确定的机械或热负载条件之前和期间,拍摄一定数量的表示样品的图像。这些图像可以是二维的、立体的或视频的,或甚至是体积图像。当然,如果期望识别具有物理尺寸的机械性能,那么对样品尺度的了解是必不可少的。
存在能够使用获得的图像来表征负载效果的各种方式。这些图像可以借助于数字图像相关性(CIN)来使用,这一方法也被缩写为“DIC”。将特别参考出版物“Imagecorrelation for shape,motion and deformation measurements:basic concepts,theory and applications”,MA Sutton,JJ Orteu,H Schreier,Springer(2009)。该方法包括在代表性运动学基础上分解样品的运动。这种基础可以是基于样品的几何网格的“有限元”类型的描述。将有负载情况下获得的多个图像与无负载情况下获得的那个或那些图像进行对比的分析使得可以测量应力情况下的样品的位移场UCIN(x),x表示样本点的坐标。图1示出了数字图像相关性方法表征的该一般原理。
同样的机械实验可以借助于有限元方法或借助于其他技术进行数值建模,以计算随后用UCAL(x)表示的位移场。
建模需要:
-样品的几何形状的知识;
-包括时间的负载条件的知识,;
-诸如区域边缘的位移或力的测量值之类的边界条件的知识,;
-样品或组成样品的各相的机械行为的规律的知识;
-针对多相介质的各相的精确位置的知识。
这些不同的元素可以经由统称为{pi}的参数集以数学方式描述,i是从1变化到n的索引,n是所讨论的参数数量。图2示出了通过建模表征的该一般原理。
当然,只有当位移场UCAL(x)与位移场UCIN(x)相同或非常接近时,建模才能真实地反映观测值。因此,参数{pi}的优化是借助于迭代循环来实现的,在迭代循环中修改该参数,直到两个位移场彼此对应到最佳可能的程度为止。图3示出了测量、建模和优化的整个过程,使得能够确定正确的参数。
这种方法有一些缺点。具体地说,图像总是有噪声的,并且DIC对这种噪声的或多或少的敏感程度取决于运动学基础的选择。在所研究区域的边缘处,这种影响尤其明显。建模的主要缺点是位移场对于某些参数的敏感性很低,使得这些参数难以确定。甚至可能在所有参数{pi}之间存在污染影响。最后,尽管存在测量缺陷,但在整个识别过程中,该位移场是不可更改的。
为了优化该方法,已经提出了各种方法。因此,名称为“Elastography imagingmodalities for characterising properties of tissue”的美国专利US 7257244描述了在弹性成像的限定背景下,在图像的相关性与所讨论样品的弹性性能的建模之间进行耦合直到获得图像之间的最佳重合的迭代方法。在该方法中,假定边界条件已知,并且在迭代期间不重新考虑边界条件。
J.Réthoré在Int.J.Num.Meth.Eng.84:631-660(2010)中发表的名称为“A fullyintegrated noise robust strategy for the identification of constitutive lawsfrom digital images”的出版物中提出了另一种方法。该方法通过引入要被最小化的泛函来关联识别与DIC,所述泛函通过强制实施两个位移场的识别来关联两个目标。识别部分通过使与平衡的偏差最小化来进行,所述与平衡的偏差是基于应用于位移的二阶微分算子的二次泛函。因此,它的最小化涉及关于位移U的四阶微分算子。
所描述的图像相关性是基于位移场的“有限元”类型离散化的全局方法。对于同一个位移场,一旦校正了当前位移场,参考图像与畸变图像之间的图像二次差以及与平衡的偏差最小化使得可以在通过DIC获得的测量值和识别结果之间找到折衷点,并且关于目标本构参数{p}对该折衷点进行优化。根据图像相关性直接且唯一地产生所传递的力不为零的边界条件。该方法经由图像的重合提供了对弹性且非线性的本构规律的识别。然而,应强调的是,边缘位移的测量受到噪声的严重影响,并且使用基于位移U的二阶导数的二次泛函(通过最小化四阶微分算子产生)导致对噪声的敏感度大大增加。用这种方法处理的示例清楚地证明在施加运动学的边缘附近的误差大大增加了。在比参考文献中使用的条件更为不利的情况下,这些边界条件甚至可能妨碍机械性能的正确确定。因此,可以观察到,即使所提议的DIC/识别耦合解决了上述一些缺点,但仍存在其他弱点,可能会限制定量识别能力。
更具体地,本发明的一个主题是一种识别经受机械应力的对象的、被称为“目标参数”的至少一个机械参数的计算机实现的方法,其特征在于,所述方法包括以下步骤:
-步骤1:通过成像装置获取应用所述机械应力之前和期间拍摄的所述对象的至少两个图像,并且测量所述对象的尺度因子;
-步骤2:计算第一泛函TCIN(UCIN),所述第一泛函TCIN(UCIN)对应于数字图像的相关性,所述数字图像的相关性取决于使用第一运动学基础表示的所述位移场UCIN,所述位移场UCIN是在有负载和无负载的情况下、在所述对象的所述图像之间的应力下、在所述对象的任意点处进行测量的;
-步骤3:计算所述对象的任意点处的计算位移场UCAL;
-步骤4:基于所述计算位移场UCAL计算第二泛函TCAL(UCAL,{p},{q}),所述计算位移场UCAL是使用第二运动学基础表示的,所述第二泛函对应于所述应力的机械模型的变分公式,所述变分公式取决于所述对象的所述几何形状、所施加的力、边界条件、至少所述目标参数{p}和预定机械参数{q};
-步骤5:计算二次范数形式的第三泛函TPAR(UCIN,UCAL),所述第三泛函TPAR(UCIN,UCAL)等于UCIN与UCAL之差;
-步骤6:将总泛函TTOT(UCIN,UCAL,{p},{q})关于UCIN、UCAL和{p}最小化,所述总泛函TTOT(UCIN,UCAL,{p},{q})至少包括以下项:
TTOT(UCIN,UCAL,{p},{q})=αTCIN(UCIN)+βTCAL(UCAL,{p},{q})+γTPAR(UCIN,UCAL)
-α、β和γ是三个非零权重系数,所述权重系数(α,β,γ)是根据与所述泛函中涉及的各种量相关联的不确定度和/或根据与泛函TTOT最小化相切的问题的条件数进行调整的。
有利的是,当对象的行为经受与时间有关的应力时,第二泛函TCAL(UCAL,{p},{q})取决于所确定的时间。
有利的是,该方法的步骤1包括力、时间或温度的附加测量值FMES,该方法的步骤3包括与附加测量值相对应的评估值FCAL,在步骤4之后是步骤4bis,步骤4bis是计算第四泛函TFOR(FCAL,FMES),第四泛函TFOR(FCAL,FMES)与这些量之间的二次偏差成比例,并且步骤5的总泛函TTOT(UCIN,UCAL,{p},{q})等于:
TTOT(UCIN,UCAL,{p},{q})=αTCIN(UCIN)+βTCAL(UCAL,{p},{q})+γTPAR(UCIN,UCAL)+χTFOR(FCAL({p},{q}),FMES)
χ是第四权重系数,权重系数χ是根据与泛函中涉及的各种量相关联的不确定度和/或根据与泛函TTOT的最小化相切的问题的条件数来调整的。
有利的是,总泛函TTOT的最小化是通过迭代方法实现的,这可能需要或不需要计算TTOT的梯度。
有利的是,第一运动学基础与第二运动学基础相同。
有利的是,测量不确定度是由具有收敛性的泛函TTOT的海塞函数通过可用的采集噪声的测量来估计的。
有利的是,第一运动学基础或第二运动学基础是在有限元网格上产生的。
本发明的另一主题是根据上述方法的适于识别至少一个机械参数的计算机设备。
本发明的另一主题是具有用于执行根据上述方法的方法的程序的计算机可读介质。
通过阅读以下不作限制的描述并且通过附图,本发明将得到更好的理解,并且其他优点将变得显而易见,其中:
已经描述的图1示出了使用称为CIN(DIC)的数字图像相关性方法进行监测来表示样本中的应力的各个步骤;
已经描述的图2示出了通过建模方法表示样本中应力的各个步骤;
已经描述的图3示出了通过根据现有技术的迭代方法优化目标参数{p}(诸如,机械性能、几何结构或边界条件)的各个步骤;
图4示出了通过根据本发明的迭代全局方法优化目标参数的各个步骤。
作为示例,图4示出了根据本发明的迭代全局方法优化目标参数的各个步骤。如果将根据本发明的方法与根据图3的现有技术的方法进行比较,则这两种方法之间的主要变化是运动学计算的位置,特别是关于边界条件的位置。在根据本发明的方法中,建模也介入到这些运动学值中,因此这些运动学值不再是常数,而是构成优化回路的一部分。
根据本发明的用于识别经受已知机械应力的试件或零件的本构材料的、被称为“目标参数”的至少一个机械参数的计算机实现的方法包括下述步骤。
第一步骤包括通过成像装置获取应用机械应力之前、期间和(对于某些应用)在施加机械应力之后拍摄的对象的数字图像,并且测量对象的尺度因子。事实上,如果想要确定具有物理尺寸的机械性能,那么对样品尺度的了解是必不可少的。
成像装置可以是单独使用或组合使用的、使得可以获得对象的至少一个图像的任何设备。通过所述装置获得的图像可以是在本领域技术人员熟知的各种波长范围内获得的光学图像。
其次,基于使用所选择的方法获得的图像来计算第一泛函TCIN。该第一泛函对应于基于位移场UCIN的数字图像的相关性,所述位移场UCIN是使用(例如,有限元网格上的)第一相关运动学基础来表示的。
传统上,该第一泛函是研究区域上的参考图像与位移场的一个或多个校正失真图像之间的平方差的和,但是可以选择诸如互相关或信息的联合熵之类的其他准则。
该方法的第一步骤可能包括诸如力、时间或温度测量值之类的附加测量值FMES。
第三步骤包括计算使用第二运动学基础表示的计算位移场UCAL,该第二泛函对应于应力的机械模型的变分公式,所述变分公式取决于对象的几何形状、所施加的力、边界条件、至少目标参数{p}和预定机械参数{q}。
实际上,同样的机械实验可以通过有限元方法或者通过其他技术来进行数值模拟,以计算随后用UCAL(x)表示的位移场。
第四步骤包括计算第二泛函TCAL(UCAL,{p},{q}),所述第二泛函TCAL(UCAL,{p},{q})取决于一种或更多种材料的机械行为、零件的几何结构、可能包括所施加的力和以目标{p}或预定{q}参数的形式聚集在一起的一个或多个所讨论时间的边界条件,以及计算目标参数的标称值集合。该方法的第四步骤可以包括对与附加测量值FMES相对应的力、时间或温度(如果它们可用)的评估值FCAL,从而可以制定与最后两个量的平方差成比例的附加泛函TFOR(FCAL,FMES),该附加泛函TFOR(FCAL,FMES)可能由测量值的方差的倒数进行加权。
第五步骤,以UCIN和UCAL之间的差的二次范数形式引入第三泛函TPAR(UCAL,UCIN)。这两个位移场与在目标参数被很好地识别且预定参数合适的情况下测量并计算的其他物理量以及所使用的模型合并在一起。
因此,所建议的识别原理是在最终步骤中,对这三个或四个泛函的权重和TTOT关于两个位移场UCAL和UCIN以及目标参数{p}进行最小化:
TTOT(UCIN,UCAL,{p},{q})=αTCIN(UCIN)+βTCAL(UCALV,{p},{q})+γTPAR(UCIN,UCAL)+χTFOR(FCAL({p},{q}),FMES)
在该方法不包括附加测量值的情况下,该泛函TTOT被简化为前三项。
应注意的是,如果应力范围、几何结构甚至图像的定义不合适,则问题可能仍然不适定的(ill-posed)。在这种情况下,并不是所有的目标参数都可以测量。对应于已识别参数与预期参数之间偏差的惩罚的Tikhonov正则化可能需要对问题进行数值求解。然后,应使用所获得的解自身的不确定度来判断该解,例如,通过考虑预先表征的测量噪声对最小化TTOT的影响而不考虑Tikhonov正则化来判断。
权重系数(α,β,γ,χ)使得能够根据与泛函中所涉及的量相关联的不确定度和/或根据与泛函TTOT最小化相切的问题的条件数来赋予各个项更大或更小的重要性。应当注意,任意选择的权重的任一者都可以设定为1。
如果不能直接获得机械模型的变分表达式,应注意的是,泛函TPAR关于UCAL的最小化可以例如在有限元代码中由线性弹性连接简单地表示,该线性弹性连接在各个节点处产生与UCAL与UCIN之间的偏差成比例的节点力。因此,对于使用当前专业计算机代码执行的有限元建模(其能够包括任意复杂的本构关系),所提出的TTOT公式在其关于UCAL的最小化中只需要在网格的各个节点处引入附加线弹性连接。所得到的解正是在UCIN、{p}和{q}固定的情况下使总泛函最小化的解。通过关于不同的未知数子集来交替进行最小化步骤,使得如果该问题是适定的,则可能达到目标最小化。
例如,可以经由连续线性化和校正,借助于牛顿-拉斐逊方法来执行泛函TTOT的最小化。
有利的是,具有收敛性的泛函TTOT的海赛函数(Hessian)使得如果例如在执行机械实验之前经由无应力情况下的重复采集获得了对采集噪声的测量值,则可以估计测量不确定度。特别地,可以理解,问题的适定性对应于严格正的特征值,并且在这种情况下,条件数对应于海赛函数的谱半径。否则,可以提出Tikhonov正则化。
有利的是,这些元素也使模型有效或无效。具体地:
-位移场UCIN使得能够估计图像相关性的剩余场,即,UCIN位移场的畸变和校正图像与参考图像之间的差;
-建模使得可以验证本构规律和平衡条件是否满足;
-可获得的任何附加测量值与从建模中得到的那些值进行比较;
-将两个位移场(一个位移场接近测量值UCIN,另一个位移场接近模型UCAL)组合在衡量这两种方法的一致性的同一个泛函中。
因此,所使用的各个泛函都提供了自己的验证。相反地,残差太大,无法与采集噪声信号模型或测量误差兼容,并且提供了关于如何丰富解释性模型或识别非预期测量伪影的指示。
有利的是,泛函的各个项的耦合使得能够补偿这个泛函或那个泛函的不适定性或病态性(ill-conditioned character)。例如,低对比度或照明不足的区域可能不允许仅使用泛函TCIN来测量UCIN。然后,泛函TPAR可以通过计算来补偿信息的缺乏。
对称地,当本构规律或几何非线性导致机械问题解的稳定性或唯一性丧失时,耦合泛函TPAR随即使得可以恢复问题的适定性,并且经由模型跟踪同一分支。
在收敛的情况下,针对位移场UCIN、不同计算出的位移场UCAL、模型FCAL的力、时间或温度或其他物理量的值、以及识别对象的目标参数{p}的估计,总泛函达到其最小值。这些参数可能是与一个或更多个相有关的材料的参数,或是对象的几何结构或其他量(例如边界条件)。
该方法由计算机实现,因此可以根据上述方法使计算机设备适于识别至少一个机械参数。
因此,根据本发明的方法基本上可以通过使用当前计算工具完全可以实现的数值计算装置来实现,但是在实践中可以分析更广泛的材料,或者使用更便宜的硬件(尤其是与给定最终结果的采集硬件的质量有关)。
Claims (10)
1.一种识别经受机械应力的对象的、被称为“目标参数”的至少一个机械参数的计算机实现的方法,其特征在于,所述方法包括以下步骤:
-步骤1:通过成像装置获取应用所述机械应力之前和期间拍摄的所述对象的至少两个图像,并且测量所述对象的尺度因子;
-步骤2:计算第一泛函TCIN(UCIN),所述第一泛函TCIN(UCIN)对应于数字图像的相关性,所述数字图像的相关性取决于使用第一运动学基础表示的位移场UCIN,所述位移场UCIN是在有负载和无负载的情况下、在所述对象的所述图像之间、在应力下的所述对象的任意点处进行测量的;
-步骤3:计算所述对象的任意点处的计算位移场UCAL;
-步骤4:基于所述计算位移场UCAL计算第二泛函TCAL(UCAL,{p},{q}),所述计算位移场UCAL是使用第二运动学基础表示的,所述第二泛函对应于所述应力的机械模型的变分公式,所述变分公式取决于所述对象的几何形状、所施加的力、边界条件、至少目标参数{p}和预定机械参数{q};
-步骤5:计算二次范数形式的第三泛函TPAR(UCIN,UCAL),所述第三泛函TPAR(UCIN,UCAL)等于UCIN与UCAL之差;
-步骤6:将总泛函TTOT(UCIN,UCAL,{p},{q})关于UCIN、UCAL和{p}最小化,所述总泛函TTOT(UCIN,UCAL,{p},{q})至少包括以下项:
TTOT(UCIN,UCAL,{p},{q})=αTCIN(UCIN)+βTCAL(UCAL,{p},{q})+γTPAR(UCIN,UCAL)
-α、β和γ是三个非零权重系数。
2.根据权利要求1所述的识别至少一个机械参数的计算机实现的方法,其特征在于,当所述对象的行为经受与时间相关的应力时,所述第二泛函TCAL(UCAL,{p},{q})取决于所确定的时间。
3.根据前述权利要求中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,所述方法的步骤1包括力、时间或温度的附加测量值FMES,所述方法的步骤3包括与所述附加测量值相对应的评估值FCAL,在步骤4之后是步骤4bis,所述步骤4bis是计算第四泛函TFOR(FCAL,FMES),所述第四泛函TFOR(FCAL,FMES)与这些量之间的二次偏差成比例,并且步骤5中的总泛函TTOT(UCIN,UCAL,{p},{q})等于:
TTOT(UCIN,UCAL,{p},{q})=αTCIN(UCIN)+βTCAL(UCAL,{p},{q})+γTPAR(UCIN,UCAL)+χTFOR(FCAL({p},{q}),FMES)
χ是第四权重系数,所述第四权重系数χ是根据与所述泛函中涉及的各个量相关联的不确定度和/或根据与所述泛函TTOT的所述最小化相切的问题的条件数来调整的。
4.根据前述权利要求中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,所述总泛函TTOT的所述最小化是通过迭代方法实现的。
5.根据权利要求1至3中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,所述总泛函TTOT的所述最小化是通过迭代方法实现的,从而需要计算TTOT的梯度。
6.根据前述权利要求中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,所述第一运动学基础与所述第二运动学基础相同。
7.根据前述权利要求中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,测量不确定度是由具有收敛性的所述泛函TTOT的海赛函数通过对可用的采集噪声的测量来估计的。
8.根据前述权利要求中的任一项所述的识别至少一个机械参数的计算机实现的方法,其特征在于,所述第一运动学基础或所述第二运动学基础是在有限元网格上产生的。
9.一种计算机设备,所述计算机设备适于识别如权利要求1至8中的任一项所述的至少一个机械参数。
10.一种计算机可读介质,所述计算机可读介质具有用于执行根据权利要求1至8中的任一项所述的方法的程序。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2018/054953 WO2019166085A1 (fr) | 2018-02-28 | 2018-02-28 | Procede mis en oeuvre par ordinateur pour l'identification de proprietes mecaniques par correlation d'images et modelisation mecanique couplees |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112055871A true CN112055871A (zh) | 2020-12-08 |
CN112055871B CN112055871B (zh) | 2024-10-11 |
Family
ID=61563381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880090358.4A Active CN112055871B (zh) | 2018-02-28 | 2018-02-28 | 通过耦合的图像相关性与机械建模来识别机械性能的计算机实现方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11215446B2 (zh) |
EP (1) | EP3759687B8 (zh) |
CN (1) | CN112055871B (zh) |
ES (1) | ES2911999T3 (zh) |
WO (1) | WO2019166085A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112597608B (zh) * | 2020-12-25 | 2023-03-14 | 清华大学 | 一种机械结合面切-法向耦合力学建模方法 |
CN113536570B (zh) * | 2021-07-16 | 2022-09-27 | 中国石油大学(北京) | 一种应力参数计算方法、装置、系统及计算机设备 |
CN115017770B (zh) * | 2022-06-10 | 2024-07-05 | 三峡大学 | 一种基于模态响应的结构损伤识别方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839699A (zh) * | 2010-05-14 | 2010-09-22 | 清华大学 | 一种测量微纳米金属互连线残余变形的方法 |
US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
CN103870690A (zh) * | 2014-03-14 | 2014-06-18 | 哈尔滨工业大学 | 一种基于有限元的超声理疗中软组织影响下的骨组织内部应力分布数值模拟方法 |
WO2014095760A1 (fr) * | 2012-12-20 | 2014-06-26 | European Aeronautic Defence And Space Company Eads France | Procédé pour l'évaluation de la performance métrologique d'une mesure d'un champ cinématique par corrélation d'images numériques |
US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
CN104603836A (zh) * | 2012-08-06 | 2015-05-06 | 范德比尔特大学 | 用于在图像引导手术期间校正数据变形的改进的方法 |
CN104700368A (zh) * | 2015-03-06 | 2015-06-10 | 南京大树智能科技股份有限公司 | 基于核函数的数字图像相关方法的位移场自适应平滑方法 |
CN105102921A (zh) * | 2012-12-20 | 2015-11-25 | 通用电气公司 | 用于监测蒸汽发生器中的操作条件的方法和系统 |
CN105765630A (zh) * | 2013-07-15 | 2016-07-13 | 空客集团有限公司 | 通过创建多个色度斑点图案测量受到应力的结构部件的形状、移动和/或变形的多尺度测量方法 |
CN106500612A (zh) * | 2016-10-12 | 2017-03-15 | 北京科技大学 | 材料构件高温变形的非接触测量系统及方法 |
CN106528982A (zh) * | 2016-10-26 | 2017-03-22 | 西安交通大学 | 一种有拉筋和围带的干摩擦阻尼失谐叶片的振动分析方法 |
CN107273590A (zh) * | 2017-05-31 | 2017-10-20 | 北京理工大学 | 一种用于复合材料本构方程参数确定的多目标优化方法 |
CN107726991A (zh) * | 2017-10-09 | 2018-02-23 | 北京航空航天大学 | 一种基于弱形式的数字图像相关应变场计算方案 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7257244B2 (en) | 2003-02-24 | 2007-08-14 | Vanderbilt University | Elastography imaging modalities for characterizing properties of tissue |
DE102007037726B4 (de) * | 2007-08-09 | 2010-07-08 | Lavision Gmbh | Verfahren zur berührungslosen Messung von Verformungen einer Oberfläche eines Messobjektes |
US8553233B2 (en) * | 2011-06-30 | 2013-10-08 | John W. Newman | Method and apparatus for the remote nondestructive evaluation of an object using shearography image scale calibration |
-
2018
- 2018-02-28 CN CN201880090358.4A patent/CN112055871B/zh active Active
- 2018-02-28 EP EP18708651.7A patent/EP3759687B8/fr active Active
- 2018-02-28 ES ES18708651T patent/ES2911999T3/es active Active
- 2018-02-28 WO PCT/EP2018/054953 patent/WO2019166085A1/fr unknown
- 2018-02-28 US US16/976,425 patent/US11215446B2/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839699A (zh) * | 2010-05-14 | 2010-09-22 | 清华大学 | 一种测量微纳米金属互连线残余变形的方法 |
US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
CN104603836A (zh) * | 2012-08-06 | 2015-05-06 | 范德比尔特大学 | 用于在图像引导手术期间校正数据变形的改进的方法 |
WO2014095760A1 (fr) * | 2012-12-20 | 2014-06-26 | European Aeronautic Defence And Space Company Eads France | Procédé pour l'évaluation de la performance métrologique d'une mesure d'un champ cinématique par corrélation d'images numériques |
CN105102921A (zh) * | 2012-12-20 | 2015-11-25 | 通用电气公司 | 用于监测蒸汽发生器中的操作条件的方法和系统 |
CN105765630A (zh) * | 2013-07-15 | 2016-07-13 | 空客集团有限公司 | 通过创建多个色度斑点图案测量受到应力的结构部件的形状、移动和/或变形的多尺度测量方法 |
US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
CN103870690A (zh) * | 2014-03-14 | 2014-06-18 | 哈尔滨工业大学 | 一种基于有限元的超声理疗中软组织影响下的骨组织内部应力分布数值模拟方法 |
CN104700368A (zh) * | 2015-03-06 | 2015-06-10 | 南京大树智能科技股份有限公司 | 基于核函数的数字图像相关方法的位移场自适应平滑方法 |
CN106500612A (zh) * | 2016-10-12 | 2017-03-15 | 北京科技大学 | 材料构件高温变形的非接触测量系统及方法 |
CN106528982A (zh) * | 2016-10-26 | 2017-03-22 | 西安交通大学 | 一种有拉筋和围带的干摩擦阻尼失谐叶片的振动分析方法 |
CN107273590A (zh) * | 2017-05-31 | 2017-10-20 | 北京理工大学 | 一种用于复合材料本构方程参数确定的多目标优化方法 |
CN107726991A (zh) * | 2017-10-09 | 2018-02-23 | 北京航空航天大学 | 一种基于弱形式的数字图像相关应变场计算方案 |
Non-Patent Citations (4)
Title |
---|
KHOO, S.-W等: "A Review of Surface Deformation and Strain Measurement Using Two-Dimensional Digital Image Correlation", 《METROLOGY AND MEASUREMENT SYSTEMS》, vol. 23, no. 3, 18 July 2016 (2016-07-18), pages 461 - 480 * |
M.BORNERT等: "Assessment of Digital Image Correlation Measurement Errors: Methodology and Results", 《EXPERIMENTAL MECHANICS》, vol. 49, 25 December 2008 (2008-12-25), pages 353 - 370 * |
代祥俊: "基于径向基函数的高精度变形测量方法研究", 《中国博士学位论文全文数据库 基础科学辑》, vol. 2017, no. 2, 15 February 2017 (2017-02-15), pages 004 - 22 * |
张立荣: "基于有限元法ZL50装载机铲斗强度分析与结构优化", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, vol. 2009, no. 1, 15 January 2009 (2009-01-15), pages 029 - 43 * |
Also Published As
Publication number | Publication date |
---|---|
EP3759687B1 (fr) | 2022-03-30 |
ES2911999T3 (es) | 2022-05-24 |
EP3759687A1 (fr) | 2021-01-06 |
US20210003389A1 (en) | 2021-01-07 |
US11215446B2 (en) | 2022-01-04 |
EP3759687B8 (fr) | 2022-06-01 |
WO2019166085A1 (fr) | 2019-09-06 |
CN112055871B (zh) | 2024-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sutton et al. | Full-field representation of discretely sampled surface deformation for displacement and strain analysis | |
CN112055871B (zh) | 通过耦合的图像相关性与机械建模来识别机械性能的计算机实现方法 | |
US10072924B2 (en) | System and method for quantifying deformation, disruption, and development in a sample | |
JP6503364B2 (ja) | 不織繊維複合布又は織繊維複合布を含む部材の検査装置及び方法 | |
Honarvar et al. | A comparison of direct and iterative finite element inversion techniques in dynamic elastography | |
Pan et al. | Some recent advances in digital volume correlation | |
Taillandier-Thomas et al. | Soft route to 4D tomography | |
JP2002532723A (ja) | ある容量内の電磁場特性を決定するための方法および装置 | |
Réthoré | Automatic crack tip detection and stress intensity factors estimation of curved cracks from digital images | |
US20180306691A1 (en) | Systems and Methods for Visualizing the Interior Morphology and Material Properties of Samples from Surface Deformations | |
Callaghan et al. | Spatial resolution and velocity field improvement of 4D‐flow MRI | |
van Dijk et al. | A global digital volume correlation algorithm based on higher-order finite elements: Implementation and evaluation | |
Villarraga-Gómez et al. | Effect of geometric magnification on dimensional measurements with a metrology-grade X-ray computed tomography system | |
Rahmani et al. | Regularized virtual fields method for mechanical properties identification of composite materials | |
Connesson et al. | Extension of the optimised virtual fields method to estimate viscoelastic material parameters from 3D dynamic displacement fields | |
Eberle et al. | A new method for estimating the bending stiffness curve of non-uniform Euler-Bernoulli beams using static deflection data | |
Albocher et al. | Adjoint-weighted equation for inverse problems of incompressible plane-stress elasticity | |
Galarce et al. | Displacement and pressure reconstruction from magnetic resonance elastography images: Application to an in silico brain model | |
JP2014021956A (ja) | 時間領域のモード解析手法,時間領域のモード解析プログラムおよび時間領域のモード解析プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
Boquet-Pujadas et al. | Reformulating optical flow to solve image-based inverse problems and quantify uncertainty | |
CN102193076B (zh) | 用于确定背景相位变化的方法和磁共振装置 | |
HajiRassouliha et al. | Subpixel measurement of living skin deformation using intrinsic features | |
Becker et al. | A robust finite element-based filter for digital image and volume correlation displacement data | |
Madani et al. | Identification of heterogeneous elastoplastic behaviors using the constitutive equation gap method | |
Bergers et al. | Enhanced global digital image correlation for accurate measurement of microbeam bending |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |