METHOD FOR MEASURING SUPER-LARGE DEFORMATION OF PLANE
20230228562 · 2023-07-20
Assignee
Inventors
Cpc classification
International classification
Abstract
A method includes the steps of arranging mark points for image recognition on a plane of a test piece to be measured; recognizing and recording positions of two-dimensional Cartesian coordinates of each mark point of the test piece to be measured before and after each stretching; and determining a deformation gradient of each mark point and deformation measurement parameters of each mark point through a numerical method, where the deformation measurement parameters include a deformation gradient matrix, an elongation tensor matrix, a finite strain tensor matrix, an orthogonal tensor matrix, an angular tensor matrix, a rotation angle, and a curvature. According to the method, objective measurement of super-large deformation of the plane relating to rotation deformation is achieved.
Claims
1. A method for measuring a super-large deformation of a plane, comprising the steps of arranging enough mark points for an image recognition on a plane of a test piece to be measured; recognizing and recording positions of two-dimensional Cartesian coordinates of each mark point on the plane of the test piece to be measured before and after each stretching; and determining a deformation gradient of each mark point and deformation measurement parameters of each mark point by using a numerical method, wherein the deformation measurement parameters comprise: at least one of an elongation tensor matrix and a finite strain tensor matrix; and at least one selected from the group consisting of an orthogonal tensor matrix, an angular tensor matrix, a rotation angle, and a curvature.
2. The method according to claim 1, wherein the curvature of one mark point is calculated according to a variation of the rotation angle of an adjacent mark point.
3. The method according to claim 1, wherein a deformation gradient matrix [F] of each mark point is calculated, if a forward difference method is used, an instance is shown as follows: any mark point is taken as P.sub.1, an adjacent mark point of a mark point in an X-axis direction is P.sub.2, an adjacent mark point in a Y-axis direction is P.sub.4, the two-dimensional Cartesian coordinates before deformation are (X.sub.1, Y.sub.1), (X.sub.2, Y.sub.2) and (X.sub.4, Y.sub.4) respectively, for the three mark points (P.sub.1, P.sub.2, and P.sub.4), the two-dimensional Cartesian coordinates after deformation are (x.sub.1, y.sub.1), (x.sub.2, y.sub.2) and (x.sub.4, y.sub.4) respectively, and the deformation gradient matrix of the mark point P.sub.1 is:
4. The method according to claim 1, wherein formulas of
[U]=√{square root over ([F].sup.T.Math.[F])}
[V]=√{square root over ([F].Math.[F.sub.].sup.T)} are configured for determining a right elongation tensor matrix [U] and a left elongation tensor matrix [V] of each mark point.
5. The method according to claim 1, wherein formulas of
[H]=ln[U]
[h]=ln[V] are configured for determining the finite strain tensor matrix of each mark point, comprising a right strain tensor matrix [H] and a left strain tensor matrix [h].
6. The method according to claim 1, wherein a formula of
[R]=[F].Math.[U].sup.−1 is configured for determining the orthogonal tensor matrix [R] of each mark point.
7. The method according to claim 1, wherein a formula of
8. The method according to claim 2, wherein components of the curvature C of mark point P.sub.1 are C.sub.11 and C.sub.12 respectively, specific calculations are as follows:
9. The method according to claim 1, wherein all the mark points are uniformly distributed on one part or all parts of the plane to be measured.
10. The method according to claim 1, further comprising the step of performing a pre-stretching operation on the plane to be measured in a natural state to obtain a plane to be measured in a tensioned state as a whole in a stretching direction, wherein coordinate positions of mark points on the plane to be measured in such a state serve as initial positions of the mark points.
11. The method according to claim 1, wherein the mark points are almost circular, and a ratio of a mark point diameter to a mark point spacing is about 1:2 to 1:4.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0035]
[0036]
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
[0045]
DETAILED DESCRIPTION OF THE EMBODIMENTS
[0046] In the example, by means of a geometric parameter measurement process of super-large deformation of a plane of a material test piece, the measurement method in the present invention is illustratively described.
[0047] 1. Manufacturing of Mark Point Test Piece
[0048] According to the example of the present invention, the designed and manufactured material standard mark point test piece is of a dumbbell-shaped sheet structure as a whole, the periphery of the material standard mark point test piece is cut neatly, and circular mark points are printed in a middle effective deformation area so as to provide regular recognition mark points for digital image measurement. The specific manufacturing method includes the steps: firstly, designing a dumbbell-shaped test piece structure as shown in
[0049] 2. Test Piece Stretching and Acquisition of Digital Images
[0050] The manufactured mark point test piece is clamped on a tensile testing machine, and the test piece is stretched and loaded as shown in
[0051] After the pre-stretching is completed, the state is marked as an initial state, and a digital image in the initial state is shown in
[0052] 3 Recognition and Coding of Test Piece Mark Points
[0053] In this embodiment, centroid position coordinates of the mark points serve as coordinate values of the mark points. In order to obtain centroid position data of the mark points of the test piece, the digital image of the mark points of the test piece is processed by using a digital image processing technology, which may employ but not limited to the following steps:
[0054] (1) Image Preprocessing
[0055] Operations such as noise suppression and filtering restoration on an original digital image acquired in the stretching process of the test piece, and real information of the position change of the mark points in the deformation process of the test piece is reduced to the maximum extent. By using an image enhancement technology, the definition and contrast of contours of mark point edges of the test piece are further improved, thereby highlighting edge information of the mark points.
[0056] (2) Image Segmentation
[0057] The mark points on the test piece image are segmented from a test piece background, and mark point edge contour information is obtained so as to provide basic data for image recognition. The main process includes the steps of firstly reading a grayscale image of the test piece, then performing edge detection by using a differential operator such as Sobel, Prewitt and Roberts, next, performing a threshold normalization operation, and finally obtaining a binary image with only black and white colors.
[0058] (3) Image Stitching
[0059] Feature point extraction is performed on each image by using several mark point images in a certain stretching state, then matching is performed on feature points one by one, next, the matching image is copied to a specific position on another image, and finally, fusion processing is performed on an overlapping boundary to form a new image with natural transition.
[0060] (4) Image Recognition
[0061] For the image subjected to mark point edge information acquisition of the test piece after image segmentation, 8 connected domain numbering and marking are performed on all pixels belonging to the same pixel connected domain, thereby returning a pixel number matrix. According to a centroid algorithm, coordinate data of each mark point centroid on the test piece image is obtained by performing loop iterative computation on all mark point edge pixel matrices on the test piece image.
[0062] (5) Mark Point Ordering and Coding
[0063] Positions of a row and column of each mark point, and the total number of the mark points on the test piece image are numbered. For coordinate data in X and Y directions, associated matching and ordering are performed according to spatial distribution positions in the test piece, and summarization is performed to form coordinate data of the mark points, rows and columns of the positions, and an n×5 matrix of serial numbers of the total number of the mark points.
[0064] More methods and technologies for recording two-dimensional Cartesian coordinate positions before and after stretching of each mark point are well known to those skilled in the art, which are not repeated herein. However, no matter what kind of coordinate positions are used to determine the technical solutions, and none of these technical solutions exceeds the protection scope of the present invention.
[0065] 4. Calculation of Metric Parameters of super-large deformation of plane
[0066] Image recognition is performed on the mark points in the image according to the digital image after stretching of the test piece to obtain plane coordinate data of each mark point. Geometric metric parameters such as the deformation gradient [F], the right elongation tensor matrix [U], the left elongation tensor matrix [V], the orthogonal tensor matrix [R], the right elongation strain matrix [H], the left elongation strain matrix [h], the angular tensor matrix [A], the rotation angle a and the curvature C of each mark point in different stretching states are calculated. Detailed description will be made below with reference to specific calculation examples.
[0067] In the case, only 4 mark points, that is P.sub.1, P.sub.2, P.sub.3 and P.sub.4, of any area inside the test piece are selected, the coordinates before deformation are (X.sub.1, Y.sub.1), (X.sub.2, Y.sub.2), (X.sub.3, Y.sub.3) and (X.sub.4, Y.sub.4), and the coordinates after deformation are (x.sub.1, y.sub.1), (x.sub.2, y.sub.2), (x.sub.3, y.sub.3) and (x.sub.4, y.sub.4). Coordinate origins of all the mark points in the plane of the test piece before and after deformation are located at lower left corners, coordinate data of mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 before and after deformation are shown in Table 1, and distribution positions in the plane of the test piece are shown in
TABLE-US-00001 TABLE 1 Coordinate data of mark points before and after deformation Coordinates before Coordinates after Serial numbers deformation deformation of mark (mm) (mm) points X Y x y P.sub.1 71.25 −9.75 148.86 −6.9214 P.sub.2 71.75 −9.75 149.87 −6.9385 P.sub.3 71.75 −9.25 149.64 −6.5764 P.sub.4 71.25 −9.25 148.62 −6.5613
[0068] 4.1 Deformation Gradient
[0069] With the coordinate positions of 4 mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 in
[0070] (1) if mark point P.sub.1 is located at a lower left corner of a mark point area or an internal mark point area, a forward difference condition is satisfied, and the calculation expression is:
[0071] (2) if mark point P.sub.2 is located at a right boundary of the mark point area and a non-upper right corner, the deformation gradient [F.sub.2] of P.sub.2 needs to be calculated according to the backward difference in an X direction and calculated according to the forward difference in a Y direction, and the calculation expression is:
[0072] (3) if mark point P.sub.4 is located at a left boundary and a non-lower left corner, the deformation gradient [F.sub.4] of P.sub.4 is calculated according to the forward difference in the X direction and is calculated according to the backward difference in the Y direction, and the calculation expression is:
and
[0073] (4) when mark point P.sub.3 is located at an upper right corner of the test piece to be measured, the deformation gradient [F.sub.3] of P.sub.3 is calculated according to the backward difference in the X direction and the Y direction, and the calculation expression is:
[0074] Calculation data of the deformation gradient of mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 of the test piece is shown in Table 2, and a calculation method for the deformation gradient of the remaining mark points on the plane of the test piece are similar to the above four types of situations.
TABLE-US-00002 TABLE 2 Deformation gradient data of mark points Serial numbers of mark points F11 F12 F21 F22 P.sub.1 2.0200 −0.4800 −0.0342 0.7202 P.sub.2 2.0200 −0.4600 −0.0342 0.7242 P.sub.3 2.0400 −0.4600 −0.0302 0.7242 P.sub.4 2.0400 −0.4800 −0.0302 0.7202
[0075] 4.2 Elongation Tensor
[0076] According to a formula:
[0077] the right elongation tensor matrix [U] is calculated to obtain U.sub.11, U.sub.12, U.sub.21 and U.sub.22.
[0078] According to a formula:
[0079] the left elongation tensor matrix [V] is calculated to obtain V.sub.11, V.sub.12, V.sub.21 and V.sub.22.
[0080] With the right elongation tensor matrix [U] and the left elongation tensor matrix [V] of mark point P.sub.1 as instances, calculation methods for the remaining mark points are the same as the above calculation methods. The data of the elongation tensor matrix of P.sub.1 in the example is shown in Table 3.
TABLE-US-00003 TABLE 3 Various elements about [U] and [V] of elongation tensor matrix of mark point P.sub.1 U.sub.11 U.sub.12 U.sub.21 U.sub.22 V.sub.11 V.sub.12 V.sub.21 V.sub.22 1.9883 −0.3581 −0.3581 0.7880 2.0708 −0.1494 −0.1494 0.7054
[0081] 4.3 Finite Strain Tensor
[0082] According to a formula:
[0083] the right strain tensor matrix [H] is calculated to obtain H.sub.11, H.sub.12, H.sub.21 and H.sub.22.
[0084] According to a formula:
[0085] the left strain tensor matrix [h] is calculated to obtain h.sub.11, h.sub.12, h.sub.21 and h.sub.22.
[0086] In the example, with the right strain tensor matrix [H] and the left strain tensor matrix [h] of mark point P.sub.1 as instances, the calculation methods for the remaining mark points are the same as the above calculation methods. The data of the strain tensor matrix of P.sub.1 in the example is shown in Table 4.
TABLE-US-00004 TABLE 4 Various elements about [H] and [h] of strain tensor matrix of mark point P.sub.1 H.sub.11 H.sub.12 H.sub.21 H.sub.22 h.sub.11 h.sub.12 h.sub.21 h.sub.22 0.6575 −0.2839 −0.2839 −0.2939 0.7229 −0.1184 −0.1184 −0.3594
[0087] 4.4 Components of Angular Tensor Matrix, Rotation Angle and Curvature
[0088] According to a formula:
[0089] the orthogonal tensor matrix [R] is calculated to obtain R.sub.11, R.sub.12, R.sub.21 and R.sub.22.
[0090] According to a formula:
[0091] the angular tensor matrix [A] is calculated to obtain an independent component A.sub.12, and the rotation angle of the mark point is α=−A.sub.12. In the example, various elements of the orthogonal tensor matrix [R] and the angular tensor matrix [A] of mark point P.sub.1 are shown in
TABLE-US-00005 TABLE 5 Various elements of orthogonal tensor matrix [R] and angular tensor matrix [A] R.sub.11 R.sub.12 R.sub.21 R.sub.22 A.sub.11 A.sub.12 A.sub.21 A.sub.22 0.9870 −0.1606 0.1606 0.9870 0 −0.1613 0.1613 0
[0092] The curvature of mark point P.sub.1 is calculated by using the difference method, and the curvature is calculated according to the forward or backward difference scheme of the deformation gradient mentioned above. According to rotation angles α.sub.1, α.sub.2 and α.sub.4 corresponding to the three mark points P.sub.1, P.sub.2 and P.sub.4, the two components C.sub.11 and C.sub.12 of mark point P.sub.1 are calculated, that is,
[0093] Calculation results of the two curvature components C.sub.11 and C.sub.12 of mark point P.sub.1 are shown in Table 6.
TABLE-US-00006 TABLE 6 Rotation angle α and two curvature components of mark point P.sub.1 Rotation Adjacent mark points Rotation angle of adjacent angle Curvature C Coordinates before deformation (mm) mark points (rad) α (deg) (rad/mm) X.sub.2 X.sub.1 Y.sub.4 Y.sub.1 α.sub.2 α.sub.1 α.sub.4 α.sub.1 α.sub.1 C.sub.11 C.sub.12 71.75 71.25 −9.25 −9.75 −0.1539 −0.1613 −0.1615 −0.1613 −9.242 0.0148 0.0004
[0094] The state of an object before deformation is of a reference configuration, and the deformed state after being stressed is a current configuration. According to the example, geometric metric parameters characterizing the deformation are calculated by using the change of the two configurations, where
[0095] (1) the deformation gradient [F] describes the degree of deformation occurring near a mass point, that is, the change degree of a line element of a deformation body, the line element before and after deformation has both the tension and pressing change and the bending change, the product of the deformation gradient is decomposed into elongation tensor and orthogonal tensor, and the elongation tensor is divided into right elongation tensor [U] and left elongation tensor [V];
[0096] (2) the right elongation tensor [U] describes the tension and pressing change before rotation of the line element, and the left elongation tensor [V] describes the tension and pressing change after rotation of the line element;
[0097] (3) the orthogonal tensor [R] describes the bending change of the line element, that is, rotation deformation;
[0098] (4) the finite strain tensor is divided into right elongation strain tensor [H] and left elongation strain tensor [h], and is a natural logarithm of the elongation tensor so as to measure the tensile deformation degree of the line element;
[0099] (5) the angular tensor [A] is a natural logarithm of the orthogonal tensor [R] and is configured to determine an axis and a rotation angle of the rotation deformation of the line element;
[0100] (6) the rotation angle a is a vector expression of the angular tensor [A]; and
[0101] (7) the curvature C is configured to measure the degree of bending or rotation deformation of the line element.
[0102] What is described above shows that the present invention relates to the tensile deformation and the rotation deformation of the super-large deformation, and objective measurement of the super-large deformation of the plane is achieved. In the example, according to an innovative plane super-large deformation geometric measurement method, the mark point test piece is designed and manufactured, each mark point and coordinates thereof are recognized by means of images before and after stretching of the test piece so as to obtain the deformation measurement parameters such as the deformation gradient, the elongation tensor, the strain tensor, the orthogonal tensor, the rotation angle and the curvature of the super-large deformation of the materials, and the problems existing in the super-large deformation geometric measurement are solved in the aspect of principles and methods.
[0103] For graphic display of the deformation measurement parameters of large deformation of the plane, the example calculates a main value and a principal axis for each tensor matrix.
[0104] Main values and principal axes of the right elongation tensor matrix [U] are calculated to obtain and the main values λ.sub.u1 and λ.sub.u2, and the corresponding principal axes are p.sub.u11, p.sub.u12, p.sub.u21, and p.sub.u22 respectively. Main values and principal axes of the left elongation tensor matrix [V] are calculated to obtain the main values lambda λ.sub.v1 and λ.sub.v2, and the corresponding principal axes are p.sub.v11, p.sub.v12, p.sub.v21, and p.sub.v22 respectively.
[0105] In the example, with mark point P.sub.1 as an instance, data of the main values and the principal axes of the elongation matrix of P.sub.1 are shown in Table 4 respectively, and it may be found that λ.sub.u1=λ.sub.v1λ.sub.1 and λ.sub.u2=λ.sub.v2=λ.sub.2 in the table.
TABLE-US-00007 TABLE 7 Main values and corresponding principal axes of elongation tensor matrix of mark point P.sub.1 Main value λ Principal axis p.sub.u Principal axis p.sub.ν λ.sub.1 λ.sub.2 P.sub.u11 P.sub.u12 P.sub.u21 P.sub.u22 p.sub.ν11 P.sub.ν12 P.sub.ν21 P.sub.ν22 0.6892 2.0871 0.2658 0.9940 −0.9640 0.2658 0.1075 0.9940 −0.9940 0.1075
[0106] Main values and principal axes of the right strain tensor matrix [H] are calculated to obtain the main values λ.sub.H1 and λ.sub.H2, and the corresponding principal axes are p.sub.H11, p.sub.H12, p.sub.H21, and p.sub.H22 respectively. Main values and principal axes of the left strain tensor matrix [h] are calculated to obtain the main values λ.sub.1 and λ.sub.h2, and the corresponding principal axes are p.sub.h11, p.sub.h12, p.sub.h21, and p.sub.h22 respectively.
[0107] In the example, with mark point P.sub.1 as an instance, data of the main values and the principal axes of the strain matrix of P.sub.1 are shown in Table 6 respectively. Obviously, λ.sub.H1=λ.sub.h1=λ.sub.1 and λ.sub.H2=λ.sub.h2=λ.sub.2 in the table.
TABLE-US-00008 TABLE 8 Main values and corresponding principal axes of strain tensor matrix of mark point P.sub.1 Main value λ Principal axis p.sub.H Principal axis p.sub.h λ.sub.1 λ.sub.2 P.sub.H11 P.sub.H12 P.sub.H21 P.sub.H22 P.sub.h11 p.sub.h12 p.sub.h21 p.sub.h22 0.6278 1.7358 −0.2658 −0.9640 −0.9640 0.2658 −0.1075 −0.9940 −0.9940 0.1075
[0108] The vector diagram of the principal axes and the main values of the right elongation tensor [U] of mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 are shown in
[0109] The vector diagram of the principal axes and the main values of the right strain tensor [H] of mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 are shown in
[0110] (3) Characterization of Rotation Angle
[0111] The rotation angles a of mark points P.sub.1, P.sub.2, P.sub.3 and P.sub.4 are included angles between the principal axes of [U] and [V], and characterization of the rotation angles α is shown in
[0112] It must be noted that the method and technology for determining deformation measurement parameters with respect to a unidirectional tensile test piece involved in the present invention is suitable for the case in which a flat test piece is stretched in any stretching direction.