CN104809362B - It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations - Google Patents

It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations Download PDF

Info

Publication number
CN104809362B
CN104809362B CN201510267349.6A CN201510267349A CN104809362B CN 104809362 B CN104809362 B CN 104809362B CN 201510267349 A CN201510267349 A CN 201510267349A CN 104809362 B CN104809362 B CN 104809362B
Authority
CN
China
Prior art keywords
sub
district
deformation
speckle pattern
strain
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.)
Expired - Fee Related
Application number
CN201510267349.6A
Other languages
Chinese (zh)
Other versions
CN104809362A (en
Inventor
王学滨
马冰
杜亚志
冯威武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201510267349.6A priority Critical patent/CN104809362B/en
Publication of CN104809362A publication Critical patent/CN104809362A/en
Application granted granted Critical
Publication of CN104809362B publication Critical patent/CN104809362B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a kind of speckle pattern preparation method comprising heterogeneous deformation based on numerical computation technology, including:Using non-classical numerical computation technology, the strain of each unit and the displacement of each node are obtained;Establish unit in numerical model or surround the certain area of node and associated with sub-district in Digital Image Correlation Method;The position of each pixel in sub-district before deformation and gray value are mapped in the sub-district after deformation.Essence of the invention is that:Result using the result of non-classical numerical computation technology as Digital Image Correlation Method, realize and undeformed speckle pattern is mapped on the grid after deformation.The present invention realizes Overall Steps using computer programming, can make to high-accuracy high-efficiency rate comprising non-homogeneous or even Deformation Localization speckle pattern, moreover, the limitation of type not loaded.The speckle pattern that the present invention makes is mainly used in field of scientific study, for example, for examining computational accuracy of Digital Image Correlation Method etc. under the conditions of non-homogeneous or even Deformation Localization.

Description

It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations
Technical field
The present invention relates to the strain measurement technique in engineering detecting, particularly one kind to include non-homogeneous or even Deformation Localization Speckle pattern high accuracy, high efficiency computer programming preparation method.
Background technology
Digital Image Correlation Method is that the two images (speckle pattern) of the body surface to being gathered before and after deformation carry out correlation Processing, to realize the measurement of object strain field and displacement field.This method light path is simple, and measuring environment requirement is simple, can be in the wild Using.Speckle can be laser-formed or artificial speckle or some natural textures.Greatly into cosmic space Galaxy moves, as low as in the measurement of nano material deformation, the shadow of all visible this method application.
Speckle pattern is a kind of medium or carrier, and the deformation or motion of object are comprised on a series of speckle patterns of shooting. Digital Image Correlation Method obtains the deformation of object and the characteristics of motion by related operation.The calculating of Digital Image Correlation Method Accuracy test be unable to do without the speckle pattern for including deformation or motion.In the lot of documents delivered, Digital Image Correlation Method Computational accuracy is examined and operated generally according to following manner:First, to a certain amount of deflection of speckle pattern pre-add or amount of exercise; Then, then with Digital Image Correlation Method calculated, differed by the deflection or amount of exercise that compare result of calculation and pre-add Degree, evaluate the precision of Digital Image Correlation Method.In addition, during related operation, some parameters, example are generally required Such as, the size of sub-district, spacing etc., influence research of these parameters to result of calculation are also required to the speckle for including deformation or motion Figure.It should be pointed out that including translation and rotation to the basic exercise of speckle pattern pre-add, the basic deformation of speckle pattern pre-add is included drawing Stretch, compress or shear.These accentuation are typically uniform, and this is suitable to the digitized map under the conditions of homogeneous deformation or rigid motion Computational accuracy inspection or parameter sensitivity Journal of Sex Research as correlation technique.But under loaded conditions, the deformation of object is not always Uniformly.By taking rock-soil material uniaxial compression as an example, with the increase of longitudinal strain, its deformation field is from being uniformly distributed to non-homogeneous point Cloth changes, until there is deformation localization phenomenon, macroscopic cracking hereafter just occurs.Shear deformation localization region is referred to as shearing Band.The distribution of shear band internal strain is very uneven.Nearest research is found, outside the shear band of sand material, strain Distribution and the uneven (numeral of Wang Xuebin, Du Yazhi, the mountain uniaxial compressions damp sand sample part of Pan one and overall volume strain Image correlation method observes Geotechnical Engineering journals, 2014,36 (9):1648-1656;Wang Xuebin, Du Yazhi, Pan Yishan, wait bases In Digital Image Correlation Method etc. different water cut sand sample shear band observation rock-soil mechanics, 2015,36 (3) under strain rate: 625-632).In order to accurately measure the strain field under the conditions of non-homogeneous or even Deformation Localization using Digital Image Correlation Method And displacement field, the computational accuracy of Digital Image Correlation Method on this condition is just have to be understood that, and the acquisition of this precision is not from Open and include non-homogeneous or even Deformation Localization speckle pattern.At present, the system of non-homogeneous or even Deformation Localization speckle pattern is included Make the also rare report of method.Some simple functions of some existing preparation method Main Basiss are made under the conditions of heterogeneous deformation Speckle pattern, without generality and wide applicability.Some documents are included using the theoretical answer making of shear band deformation field (Wang Xuebin, Du Yazhi, the mountain of Pan one consider that single order is related to the digital picture of Second Order Displacements gradient to the speckle pattern of Deformation Localization Comparison engineering mechanics of the method in shear band measurement, 2013,30 (7):282-287), simple shear condition, nothing are only applicable to Legal system makees the speckle pattern for including heterogeneous deformation outside shear band.In fact, the theoretical answer of existing shear band deformation field is very It is limited, due to the complex nature of the problem, under the relatively simple load-up condition such as compression shear, direct shearing, Biaxial Compression, shearing The acquisition of theoretical answer with interior deformation field is also extremely difficult.
The invention provides a kind of high accuracy comprising non-homogeneous or even Deformation Localization speckle pattern, high efficiency computer Program preparation method.Essence of the invention is that:Using the result of non-classical numerical computation technology as Digital Image Correlation Method Result, realize by undeformed speckle pattern be mapped to deformation after grid on.This method is due to introducing non-classical numerical value meter Calculation technology, it is no longer necessary to any theoretical answer about shear band deformation field, various simple load can be made using this method And include non-homogeneous or even Deformation Localization speckle pattern under the conditions of combined load.
The content of the invention
Under the conditions of non-homogeneous or even Deformation Localization, it is used for the inspection of Digital Image Correlation Method computational accuracy to make With the speckle pattern of calculating parameter sensitivity study, heterogeneous deformation is included based on numerical computation technology the invention provides a kind of Speckle pattern preparation method, the result of non-classical numerical computation technology pass through foundation by the result as Digital Image Correlation Method Unit or surround the certain area of node in numerical model and associated with sub-district in Digital Image Correlation Method, obtain comprising it is non- Even or even Deformation Localization speckle pattern.The present invention realizes Overall Steps using computer programming, substantially increases comprising non-equal The precision and efficiency that even or even Deformation Localization speckle pattern makes, moreover, the limitation of type not loaded, is carried suitable for simple Lotus type (for example, stretching, compression and shearing etc.), it is also applied for combined load type (for example, draw cut with compression shear etc.).
To solve the above problems, the invention provides a kind of non-homogeneous speckle pattern making side based on numerical computation technology Method, it is characterised in that including:
Using non-classical numerical computation technology, the strain of each unit and the displacement of each node are obtained;
Establish unit in numerical model or surround the certain area of node and associated with sub-district in Digital Image Correlation Method;
The position of each pixel in sub-district before deformation and gray value are mapped in the sub-district after deformation.
It is described to be further using the strain of non-classical numerical computation technology acquisition each unit and the displacement of each node:
The second order of Internal length and strain gradient is introduced in the yield function of the elastic and plastic strain softening constitutive model , two-dimension numerical calculation model is created, obtains the strain of each unit and the displacement of each node.
Further, wherein,
The two-dimension numerical calculation model is made up of some square shaped cells, is connected between unit by node, unit chi It is very little to be equal to sub-district size.
Further, wherein,
The yield function f has following form:
Wherein, F (σ) is yield stress index,To soften rule,For non local softening parameter, can be written as follows Form:
Wherein, k is local softening parameter, and c is the material parameter proportional to material internal length parameter,It is general to draw Laplacian operater, under two-dimensional condition, its expression formula is:
Further, wherein,
The parameter relevant with length is introduced in the yield function, is graded numerical computations using finite element or finite difference When technology is solved, the width for the band that localizes will be now uniquely determined, and (i.e. grid is quick for the different phenomenon of grid difference result of calculation Perception) it will disappear, the uniqueness of numerical result and objective performance are protected.
Sub-district in the certain area established unit in numerical model or surround node and Digital Image Correlation Method Associate and be further:
The displacement of sub-district central point is obtained, is either directly obtained by numerical computation technology or is saved by 4 of 1 unit The displacement result of point takes average acquisition;The strain of sub-district is obtained, is either directly obtained or by some lists by numerical computation technology The strain result of member takes average acquisition, or the displacement field directly obtained by numerical technique to be obtained by centered Finite Difference Methods.
Further, wherein,
The displacement result of 4 nodes by 1 unit takes the specific formula of the displacement of average acquisition sub-district central point For:
Wherein, the displacement that u, v are respectively sub-district central point on both horizontally and vertically.
Further, wherein,
It is described take average acquisition sub-district by the strain result of some units the specific formula of strain be:
Wherein, εx、εy、γxyRespectively horizontal line strain, vertical line strain and the shear strain of sub-district, n are sub-district center Number of unit around point.When the central point of sub-district is located at the inside of computation model, n=4;When the central point of sub-district is positioned at meter When calculating on the border of model, n=1 or 2.The central point that n=1 represents sub-district takes 4 jiaos in computation model, in this special case, Unit around sub-district central point only has 1, and the strain of sub-district is the strain of this unit.
Further, wherein,
It is described to be by the specific formula of centered Finite Difference Methods acquisition strain field by displacement field:
Wherein, Δ is the distance between node, and (i, j) represents the i-th row, jth row.
It is described by the position of each pixel in sub-district before deformation and gray value be mapped to deformation after sub-district on be further:
The position of each pixel and gray value in sub-district are obtained before deformation, obtains the position of each pixel in sub-district after deformation, will The gray value of each pixel is transformed in the sub-district after deformation in sub-district before deformation.
Further, wherein,
The concrete methods of realizing that the gray value by each pixel in sub-district before deformation is transformed in the sub-district after deformation can To be classified as two classes:Direct mapping and back mapping.
Further, wherein,
The Direct mapping is by being circulated successively to each pixel in sub-district before deformation, obtaining its sub-district after deformation Interior position, specific formula are:
Wherein, Q (x, y) is the coordinate of any pixel central point in sub-district before deformation, Q'(x', y') for the point after deformation Coordinate in sub-district, u, v are respectively sub-district center point P (x0, y0) on both horizontally and vertically displacement (need to obtain the metric system with Conversion relation between pixel, and then the metric system displacement is converted into pixel displacement), u=x'0–x0, v=y'0–y0, Δ x, Δ y Any pixel central point Q (x, y) to sub-district center point P (x in sub-district before respectively deforming0, y0) horizontal and vertical distance (figure 1-a and Fig. 1-b).Various partial derivatives describe the deformation of sub-district,If calculate Position be located exactly on the central point of a certain pixel, obtain the gray value of the pixel, in most cases, calculate The position gone out will fall in sub-pixel location, and now, row interpolation is entered in reply gray scale field, by the gray scale of each pixel in sub-district before deformation Value is transformed in the sub-district after deformation.
Further, wherein,
The back mapping is by being circulated successively to each pixel in sub-district after deformation, obtaining its sub-district before being deformed Interior position.The formula of back mapping is:
Wherein, Q'(x', y') for after deformation in sub-district any pixel central point coordinate, Q (x, y) be the point before being deformed Coordinate in sub-district.
Further, wherein the back mapping is before proceeding, it is necessary to calculate the coordinate of 4 angle points of sub-district after deformation.
It is of the present invention a kind of non-homogeneous speckle pattern preparation method, non-classical numerical value to be included based on numerical computation technology The result of computing technique is by the result as Digital Image Correlation Method, by establishing unit in numerical model or surrounding node Certain area associates with sub-district in Digital Image Correlation Method, and acquisition includes non-homogeneous or even Deformation Localization speckle pattern. The present invention realizes Overall Steps using computer programming, substantially increases comprising non-homogeneous or even Deformation Localization speckle pattern system The precision and efficiency of work, moreover, the limitation of type not loaded, suitable for simple load type (for example, stretching, compressing and cutting Cut), it is also applied for combined load type (for example, draw cut with compression shear etc.).
Brief description of the drawings
Fig. 1-a are the schematic diagram of sub-district before the deformation in Digital Image Correlation Method, and Fig. 1-b are the signal of sub-district after deformation Figure;
In figure, 1 is the image before being deformed in Digital Image Correlation Method, and 2 be a sub-district in image, P before deforming (x0, y0) for the central point of sub-district before deformation, Q (x, y) is the central point of any pixel in sub-district before deformation, and 3 be digital picture phase Image after being deformed in the method for pass, 4 be a sub-district in image, P'(x' after deforming0, y'0) for deformation after sub-district center Point, Q'(x', y') for any point in sub-district after deformation, Δ x, Δ y be respectively any pixel central point Q in sub-district before deforming (x, Y) to sub-district center point P (x0, y0) horizontal and vertical distance, u, v are respectively sub-district center point P (x0, y0) in horizontal and vertical Displacement on direction;
Fig. 2 is the schematic diagram calculation of the embodiment of the present invention one;
In figure, 5 be mathematical calculation model, and 6 be a unit in mathematical calculation model, its size and sub-district size phase Together, 7 be unit in numerical model node, the arrow at node represents the displacement of node, and 8 be the center of unit in numerical model Point;
Fig. 3-a are grid chart before the deformation of numerical model, and Fig. 3-b and Fig. 3-c are different speckle patterns;
Grid chart after the deformation of numerical model when Fig. 4-a are longitudinal strain=0.025, Fig. 4-b be longitudinal strain= Grid chart after the deformation of numerical model when 0.0375, grid chart after the deformation of numerical model when Fig. 4-c are longitudinal strain=0.05, Grid chart after the deformation of numerical model when Fig. 4-d are longitudinal strain=0.0625;
Fig. 5-a, Fig. 5-b, Fig. 5-c, Fig. 5-d are respectively to transform to Fig. 3-b in Fig. 4-a, Fig. 4-b, Fig. 4-c, Fig. 4-d Result figure on grid;
Fig. 6-a, Fig. 6-b, Fig. 6-c, Fig. 6-d are respectively to transform to Fig. 3-c in Fig. 4-a, Fig. 4-b, Fig. 4-c, Fig. 4-d Result figure on grid;
Fig. 7 is the schematic diagram calculation of the embodiment of the present invention two;
In figure, 9 be the certain area that node is surrounded in mathematical calculation model, and its size is identical with sub-district size.
Embodiment
Technical scheme is described in detail below in conjunction with accompanying drawing.
Illustrated using two kinds of specific embodiments, the difference of two schemes is Digital Image Correlation Method Sub-district 2 is different from the corresponding region in numerical model 5 before deformation in (its general principle is shown in Fig. 1-a and Fig. 1-b), embodiment one A unit 6 before deformation in middle Digital Image Correlation Method in the corresponding numerical model of sub-district 2, the center of sub-district 2 before deformation The central point 8 of unit 6 in the corresponding numerical models of point P (x, y), sub-district before the deformation in embodiment two in Digital Image Correlation Method The certain area 9 of encirclement node in 2 corresponding numerical models, before deformation in the corresponding numerical model of the center point P (x, y) of sub-district 2 One node 7 of unit 6.
As shown in Fig. 2 be the schematic diagram calculation of the embodiment of the present invention one, including step:
Step 1:Using non-classical numerical computation technology, the strain of each unit and the displacement of each node are obtained.
The step for implement it is specific as follows:First, in the yield function (formula (1-3)) of the elastic and plastic strain softening constitutive model The middle second order term for introducing Internal length and strain gradient, to avoid the mesh dependence of numerical result, it is ensured that numerical value The objectivity and uniqueness of result of calculation;Then, two-dimension numerical calculation model 5 (Fig. 3-a) is created, the square shaped cells length of side is 5mm, the parameter in boundary condition, loading environment and the constitutive model of computation model is given, obtain the strain of each unit 6 and each section The displacement (Fig. 4-a, Fig. 4-b, Fig. 4-c and Fig. 4-d) of point 7.
Step 2:Establish the association of sub-district in unit and Digital Image Correlation Method in numerical model.
The step for specific implementation it is as follows:First, according to the displacement of each node 7 in numerical model 5, obtained according to formula (4) The displacement of the central point 8 of each unit 6, the strain of each unit are obtained by the first step;Then, select consistent with computation model size A second best in quality speckle pattern, as the speckle pattern 1 before deformation, speckle pattern can be according to relevant algorithm manual manufacture, also can be right Rock or other materials surface with native texture are taken pictures, by intercept a part be used as speckle pattern (Fig. 3-b and Fig. 3- C), the size of sub-district 2 is equal to the size of unit 6, is 11 × 11 pixels;Finally, the central point 8 of unit is appointed as deforming preceding son The center point P (x, y) in area 2, the size of unit 6 are appointed as the size of sub-district 2 before deforming, and the displacement of unit center point 8 is appointed as The displacement of the center point P (x, y) of sub-district 2 before deformation, the strain of unit 6 are appointed as the strain of sub-district 2 before deforming.
Step 3:The gray value of each pixel in sub-district before deformation is mapped in the sub-district after deformation.
The step for embodiment it is as follows:First, calculated according to the displacement of each node 7 of unit in numerical model 5 The position of 4 angle points of sub-district after deformation;Then, the gray value of each pixel Q (x, y) in the sub-district 2 in image before deformation 1 is converted In sub-district 4 after to deformation in image 3.
The specific implementation method of said process has two kinds:Direct mapping and back mapping.Direct mapping is to scan in order Each pixel Q (x, y) before deformation in sub-district 2, its position in sub-district 4 after deformation is then calculated according to formula (7), passed through Interpolation transforms to the gray value of each pixel Q (x, y) in the sub-district 2 in image before deformation 1 in the sub-district 4 after deformation in image 3. Now, the number of the pixel in unit length, i.e. conversion relation between the metric system and pixel need to be obtained.If the position calculated It is located exactly on the central point of a certain pixel, obtains the gray value of the pixel, in most cases, the position calculated Putting to fall in sub-pixel location, and now, row interpolation is entered in reply gray scale field.Back mapping is sub-district 4 after scanning distortion in order Interior each pixel Q'(x', y'), its position in sub-district 2 before being deformed is calculated according to formula (8), schemed before being deformed by interpolation Transformed to as the gray value of each pixel Q (x, y) in the sub-district 2 in 1 in the sub-district 4 after deformation in image 3, as a result such as Fig. 5-a, figure 5-b, Fig. 5-c, Fig. 5-d, Fig. 6-a, Fig. 6-b, Fig. 6-c, shown in Fig. 6-d.
As shown in fig. 7, be the schematic diagram calculation of the embodiment of the present invention two, including step:
Step 1:It is identical with the step 1 of embodiment one.
Step 2:The association of sub-district in the certain area for surrounding node and Digital Image Correlation Method is established in numerical model.
The step for specific implementation it is as follows:First, according to the strain for the element number and each unit for surrounding any node, root The strain in the region 9 for surrounding any node is obtained according to formula (5), the displacement of each node is obtained by the first step;Secondly, selection and meter The consistent the second best in quality speckle pattern of moulded dimension is calculated, as the speckle pattern 1 before deformation;Finally, the position of node 7 is referred to It is set to the center point P (x, y) of sub-district 2 before deforming, the size for surrounding the region 9 of node is appointed as the size of sub-district 2, node 7 Displacement is appointed as the displacement of the center point P (x, y) of sub-district 2, and the strain for surrounding the region 9 of any node is appointed as the strain of sub-district 2. Another method for obtaining the strain in the region for surrounding any node is to carry out the equation of the ecentre to the displacement field of node according to formula (6) Point, so, it is not necessary to which the strain of some units to surrounding any node is averaged.
Step 3:It is identical with the step 3 of embodiment one.
It should be pointed out that because the unit number for surrounding node on border may be 1 or 2, so, node on model boundary is surrounded Region size may be surround model internal node region size half or a quarter, this can make mapping process Become somewhat complicated.In order to avoid drawbacks described above, in mapping process, it may be considered that only to the node inside those encirclement models Region operated, the node on model boundary is not operated.So, the speckle pattern after the deformation after mapping will be than becoming The small circle of speckle pattern before shape.This processing mode also has another advantage:Centered Finite Difference Methods are easily implemented, otherwise, it is necessary to To carrying out centered difference again after the positional displacement interpolation of the node on model boundary.
By contrast, embodiment one is most preferred embodiment, can guarantee that numerical model is intactly mapped to speckle pattern, and It is and easy to implement.
It is of the present invention it is a kind of non-homogeneous speckle pattern preparation method is included based on numerical computation technology, with non-classical number It is worth result of the result of computing technique as Digital Image Correlation Method, by establishing unit in numerical model or surrounding node Certain area associates with sub-district in Digital Image Correlation Method, and acquisition includes non-homogeneous or even Deformation Localization speckle pattern, The main application of the speckle pattern is field of scientific study:Under the conditions of non-homogeneous or even Deformation Localization, digital picture is related The computational accuracy of method is examined and calculating parameter sensitivity study, in addition, it may also be used for the purpose of decoration or image display.This hair It is bright to realize Overall Steps using computer programming, substantially increase what is made comprising non-homogeneous or even Deformation Localization speckle pattern Precision and efficiency, moreover, the limitation of type not loaded, is widely used.

Claims (3)

  1. A kind of 1. speckle pattern preparation method comprising heterogeneous deformation based on numerical computation technology, it is characterised in that including:
    Using numerical computation technology, by introducing Internal length in the yield function of the elastic and plastic strain softening constitutive model With the second order term of strain gradient, the two-dimension numerical calculation model being made up of some square shaped cells is created, passes through section between unit Point is connected, and unit size is equal to sub-district size, obtains the strain of each unit and the displacement of each node;
    Establish unit in two-dimension numerical calculation model or surround the certain area of node and sub-district in Digital Image Correlation Method Association;
    The position of each pixel in sub-district before deformation and gray value are mapped in the sub-district after deformation.
  2. A kind of 2. speckle pattern making side comprising heterogeneous deformation based on numerical computation technology according to claim 1 Method, it is characterised in that described to establish the certain area and digital picture phase of unit or encirclement node in two-dimension numerical calculation model The association of sub-district, is further divided into the method for pass:
    The displacement of sub-district central point is obtained, is either directly obtained by numerical computation technology or 4 nodes by 1 unit Displacement result takes average acquisition;The strain of sub-district is obtained, is either directly obtained or by some units by numerical computation technology Strain result takes average acquisition, or the displacement field directly obtained by numerical computation technology to be obtained by centered Finite Difference Methods.
  3. A kind of 3. speckle pattern making side comprising heterogeneous deformation based on numerical computation technology according to claim 1 Method, it is characterised in that it is described that the position of each pixel in sub-district before deformation and gray value are mapped in the sub-district after deformation, enter one Step is divided into:
    The position of each pixel and gray value in sub-district are obtained before deformation, the position of each pixel in sub-district after deformation is obtained, will deform The gray value of each pixel is transformed in the sub-district after deformation in preceding sub-district.
CN201510267349.6A 2015-05-22 2015-05-22 It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations Expired - Fee Related CN104809362B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510267349.6A CN104809362B (en) 2015-05-22 2015-05-22 It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510267349.6A CN104809362B (en) 2015-05-22 2015-05-22 It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations

Publications (2)

Publication Number Publication Date
CN104809362A CN104809362A (en) 2015-07-29
CN104809362B true CN104809362B (en) 2018-04-06

Family

ID=53694178

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510267349.6A Expired - Fee Related CN104809362B (en) 2015-05-22 2015-05-22 It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations

Country Status (1)

Country Link
CN (1) CN104809362B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105157594B (en) * 2015-09-05 2018-02-23 辽宁工程技术大学 A kind of Digital Image Correlation Method based on half son area split plot design
CN105654444B (en) * 2016-01-27 2018-07-13 辽宁工程技术大学 A kind of restored method of warped digital pattern
CN106152964B (en) * 2016-07-19 2018-08-07 辽宁工程技术大学 A kind of measuring method of plane planted agent's variability field based on speed linearity fitting
CN106767481B (en) * 2016-12-30 2018-12-28 辽宁工程技术大学 A kind of half sub-district related optical measurement method of strain localization band internal strain field
CN110188759A (en) * 2019-06-21 2019-08-30 江苏开放大学(江苏城市职业学院) One kind strain field sub-district dynamic selection method in loading by means of digital image correlation method
CN110926356B (en) * 2019-12-05 2021-11-02 武汉科技大学 Speckle preparation method for high-temperature deformation test of shaped refractory material
CN111369421A (en) * 2020-02-25 2020-07-03 东南大学 Image compression-based GPU-CPU joint calculation digital image method
CN111369549B (en) * 2020-03-10 2022-11-08 北京大学 Digital image deformation characterization method and device, electronic equipment and medium
CN112819962B (en) * 2021-02-10 2024-04-16 南京航空航天大学 Non-uniform grid division and local grid density method in digital image correlation
CN114689644B (en) * 2022-03-29 2023-03-31 清华大学 High-temperature environment parameter measuring method and device
CN116642750B (en) * 2023-07-24 2023-10-20 长江三峡集团实业发展(北京)有限公司 Rock strain localization starting time prediction method, device and equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564856A (en) * 2012-01-09 2012-07-11 西安交通大学 M integral measurement method based on plastic multi-defect material relevant to digital image
CN103149086A (en) * 2013-01-15 2013-06-12 华南理工大学 Device and method for measuring toughness damage evolution in process of metal uniaxial tension
CN103808567A (en) * 2014-02-28 2014-05-21 中国石油大学(华东) Mechanical property testing device and mechanical property testing method for soldered joint

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564856A (en) * 2012-01-09 2012-07-11 西安交通大学 M integral measurement method based on plastic multi-defect material relevant to digital image
CN103149086A (en) * 2013-01-15 2013-06-12 华南理工大学 Device and method for measuring toughness damage evolution in process of metal uniaxial tension
CN103808567A (en) * 2014-02-28 2014-05-21 中国石油大学(华东) Mechanical property testing device and mechanical property testing method for soldered joint

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于微粒子群优化算法的数字散斑图像相关方法;梁智锦等;《激光技术》;20140930;第38卷(第5期);摘要、第604页第1.2-1.3节、第605页第2.1节 *

Also Published As

Publication number Publication date
CN104809362A (en) 2015-07-29

Similar Documents

Publication Publication Date Title
CN104809362B (en) It is a kind of to include heterogeneous deformation speckle pattern preparation method based on numerical computations
US10048226B2 (en) Imaging method and apparatus based on magnetic flux leakage testing
CN110599528B (en) Unsupervised three-dimensional medical image registration method and system based on neural network
CN106600683B (en) A kind of adaptive method for reconstructing of grid model towards Bone CT sequence image
US8803943B2 (en) Formation apparatus using digital image correlation
CN108038902A (en) A kind of high-precision three-dimensional method for reconstructing and system towards depth camera
CN109360246B (en) Stereoscopic vision three-dimensional displacement measurement method based on synchronous subarea search
CN105678757B (en) A kind of ohject displacement measuring method
CN111340944B (en) Single-image human body three-dimensional reconstruction method based on implicit function and human body template
CN103076347B (en) Measurement method for mechanical injury of brittle material based on in-situ X-ray tomography
Li et al. Measuring unsaturated soil deformations during triaxial testing using a photogrammetry-based method
CN105654483A (en) Three-dimensional point cloud full-automatic registration method
CN107271460A (en) A kind of porous material internal moisture saturation degree changes the quantitatively characterizing method of spatial distribution
Wang et al. Mesh-based digital image correlation method using non-uniform elements for measuring displacement fields with high gradient
CN108871197B (en) Displacement field measuring method, device, equipment and storage medium for material surface
CN104239711B (en) A kind of determination method of research of joints roughness coefficient of rock
Dufour et al. 3D surface measurements with isogeometric stereocorrelation—application to complex shapes
JP6681221B2 (en) Structure analysis device, structure analysis method, and three-dimensional woven fiber material manufacturing method
CN101980304A (en) Three-dimensional digital volume image distortion measuring method
CN109584368B (en) Method and device for constructing three-dimensional structure of biological sample
CN107240103A (en) Boundary processing method in a kind of digital volume related algorithm split based on image
CN103116888A (en) Method for solving intrinsic parameters of cameras by plane triangles
CN106815871B (en) Modeling method of scanning electron microscope imaging system
CN109829939B (en) Method for narrowing search range of multi-view image matching same-name image points
CN113344983B (en) Multi-point cloud registration method based on Ping Miandian cloud segmentation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180406

Termination date: 20200522

CF01 Termination of patent right due to non-payment of annual fee