Volume coefficient matrix rapid solving method for residual stress field inference and application

Abstract

The invention provides a volume coefficient matrix rapid solving method for residual stress field inference and application, and relates to the technical field of material stress measurement and characterization. The method comprises the following steps: S1, calculating an inverse matrix of a non-constraint node sub-stiffness matrix of an initial state of a part; s2, calculating a non-constraint node sub-stiffness matrix variable quantity caused by part material removal in each process step; s3, calculating an inverse matrix of a non-constraint node sub-stiffness matrix of the part in each process step; s4, calculating a volume coefficient matrix of the part in each process step; and S5, assembling the volume coefficient matrixes of the parts in each process step into a volume coefficient matrix. The method solves residual stress in residual stress field inference by utilizing a volume coefficient matrix. According to the method, the calculation efficiency of the volume coefficient matrix is remarkably improved, and the problem of low calculation efficiency of the volume coefficient matrix in the prior art is effectively solved.

Description

Technical Field

[0001] The present invention relates to the technical field of material stress measurement and characterization, and in particular to a method for quickly solving a volume coefficient matrix for residual stress field inference and its application. Background Art

[0002] Against the backdrop of rapid advancements in aerospace technology, new-generation aircraft structural components are becoming increasingly complex and sophisticated, while simultaneously requiring ever-shorter manufacturing cycles. This places higher demands on machining quality. However, the workpiece undergoes multiple steps during the manufacturing process, including heat treatment and machining, which inevitably generate residual stresses in the blank. Residual stresses significantly impact machining accuracy, material strength, dimensional stability, and fatigue resistance, and can even cause deformation or failure. Therefore, accurately assessing the residual stress field distribution and optimizing the machining process have become critical steps in improving part machining quality and a key research focus in the current manufacturing industry.

[0003] After searching existing technical documents, it was found that patent document CN114154364A discloses a method for inferring the initial residual stress field based on deformation force. This method uses the influence of the unbalanced stress field in the non-cutting state under the equivalent clamping constraint of the deformation force on the deformation trend. According to the principle of virtual work, the mechanical relationship between a series of deformation forces and residual stress fields monitored in multiple non-cutting states is established, namely (in, represents the volume coefficient matrix; represents the residual stress, represents the deformation force). This method is an effective means of inferring the residual stress field. However, the calculation of the volume coefficient matrix, a key parameter for inferring residual stress from deformation force, relies entirely on finite element simulation, resulting in low computational efficiency and difficulty meeting the requirements of engineering applications.

[0004] Therefore, there is an urgent need for an innovative method that can significantly improve the efficiency of volume coefficient matrix calculation to break through the limitations of existing technologies. Summary of the Invention

[0005] The purpose of the present invention is to provide a method and application for quickly solving the volume coefficient matrix for residual stress field inference, which can solve the problem of low efficiency in volume coefficient matrix calculation in the prior art and reduce calculation time.

[0006] To achieve the above object, the present invention provides a method for quickly solving the volume coefficient matrix for residual stress field inference, comprising the following steps: Step S1: Calculate the inverse matrix of the non-constrained node sub-stiffness matrix of the initial state of the part; Step S2: Calculate the change in the sub-stiffness matrix of the unconstrained nodes caused by the material removal of the parts in each process step; Step S3: Calculate the inverse matrix of the unconstrained node sub-stiffness matrix of each step part; Step S4: Calculate the volume coefficient matrix of the parts in each process step; Step S5: Assemble the volume coefficient matrices of the parts in each process step into a volume coefficient matrix.

[0007] Preferably, in step S3, the specific expression of the inverse matrix of the unconstrained node sub-stiffness matrix of each step part is: ; ; in, Indicates the The overall stiffness matrix of the work step; Indicates the Unconstrained node sub-stiffness matrix of the work step parts; Indicates the Unconstrained-constrained node sub-stiffness matrix of the work step part; Indicates the Constrained-unconstrained node sub-stiffness matrix of the work step parts; Indicates the The constraint node sub-stiffness matrix of the work step parts; represents the identity matrix; Indicates the The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the The change in the stiffness matrix of the unconstrained nodes caused by the material removal of the work step parts; Indicates the number of work steps.

[0008] Preferably, in step S3, Based on the idea of ​​block matrix, through elementary row transformation, Convert to upper triangular matrix , the specific expression is: ; in, express The number of rows of the block matrix composed of non-zero elements in .

[0009] Accordingly, through the conversion The same elementary row transformation will Convert to , the specific expression is: ; ; ; ; in, Indicates the first The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the first -1 Inverse matrix of the non-constrained node stiffness matrix of the step part; express The block sub-matrix of ; express The block sub-matrix of ; Represents the row block index, ; Represents the column block index, ; The number of rows is , The number of rows and columns is based on Perform calculations; represents transpose; Finally, through the conversion The opposite elementary row transformation will Transform again to get .

[0010] Preferably, in step S4, the specific expression of the volume coefficient matrix of the parts in each process step is: ; in, Indicates the The sub-geometry matrix corresponding to the work step constraint node; Indicates the The sub-geometry matrix corresponding to the non-constrained nodes of the work step; Indicates the Volume coefficient matrix of process step parts.

[0011] Preferably, in step S5, the specific expression of assembling the volume coefficient matrix of the parts in each process step into the volume coefficient matrix is: ; in, represents the volume coefficient matrix.

[0012] The present invention also provides an application of the above-mentioned method for quickly solving the product coefficient matrix in residual stress field inference, using the volume coefficient matrix to solve the residual stress. The specific expression is: ; in, represents residual stress; Indicates deformation force.

[0013] Therefore, the present invention adopts the above-mentioned method and application for quickly solving the volume coefficient matrix for residual stress field inference, and the beneficial technical effects are as follows: The present invention reduces the dimension of matrix inversion in the process of solving the volume coefficient matrix through block matrix inversion and iterative update algorithm, significantly improves the calculation efficiency of the volume coefficient matrix, and effectively solves the problem of low calculation efficiency of the volume coefficient matrix in the prior art. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] Figure 1 It is the parts model drawing; Figure 2 It is the finite element model diagram of some working steps; among them, Figure 2 (a) is the finite element model of the second step; Figure 2 (b) is the finite element model of the 5th step; Figure 2 (c) is the finite element model of the 10th step; Figure 2 (d) is the finite element model of the 16th step; Figure 3 Comparison chart of residual stress results; Figure 4 The present invention provides a flow chart of a method for quickly solving a volume coefficient matrix for residual stress field inference and its application. DETAILED DESCRIPTION

[0015] The technical solution of the present invention is further described below with reference to the accompanying drawings and embodiments.

[0016] Unless otherwise defined, technical or scientific terms used in the present invention shall have the same meaning as commonly understood by one of ordinary skill in the art to which the present invention belongs.

[0017] Example 1

[0018] like Figure 1 As shown in the figure, the part selected in this embodiment is a rectangular blank of aluminum alloy with a length of 32 mm, a width of 16 mm, and a height of 20 mm, and its material density is 2830 kg / m 3 , Young's modulus is 71.7 GPa, and Poisson's ratio is 0.33.

[0019] Three fixed constraint points were selected at the bottom of the part and fixed constraints were applied to simulate fixed clamping. Spring element points in the z direction were set at the four corners of the part to simulate the constraints of the deformation force monitoring device on the workpiece and obtain deformation force data.

[0020] The part is divided into 2mm grids and 16 areas are selected as material removal areas. A material removal operation is performed in each process step. The material removal process of each process step is simulated by the function of ABAQUS finite element software. After completing the boundary condition setting, the finite element model file in inp format is exported. The finite element model of some process steps is as follows Figure 2 As shown. Among them, Figure 2 (a) is the finite element model of the second step; Figure 2 (b) is the finite element model of the 5th step; Figure 2 (c) is the finite element model of the 10th step; Figure 2 (d) in the figure is the finite element model of the 16th process step.

[0021] Based on the finite element model, the following data were extracted: .

[0022] in, Indicates the The overall stiffness matrix of the work step; Indicates the The sub-geometry matrix corresponding to the work step constraint node; Indicates the The sub-geometry matrix corresponding to the non-constrained nodes of the work step; , Indicates the number of work steps.

[0023] like Figure 4 As shown, the present invention proposes a method for quickly solving the volume coefficient matrix for residual stress field inference, comprising the following steps: Step S1: According to the node constraints, the overall stiffness matrix is ​​divided into four sub-stiffness matrices. The specific expression is: ; in, Indicates the Unconstrained node sub-stiffness matrix of the work step parts; Indicates the Unconstrained-constrained node sub-stiffness matrix of the work step part; Indicates the Constrained-unconstrained node sub-stiffness matrix of the work step parts; Indicates the The constraint node sub-stiffness matrix of the work step part.

[0024] The stiffness matrix of the unconstrained nodes in the initial state of the part , calculate the inverse matrix of the unconstrained node stiffness matrix of the initial state of the part .

[0025] Step S2: Calculate the change in the unconstrained node sub-stiffness matrix caused by the material removal of each step. The specific expression is: ; in, Indicates the The change in the stiffness matrix of the unconstrained node caused by the material removal of the work step part.

[0026] In this embodiment, .

[0027] Step S3: Calculate the inverse matrix of the unconstrained node sub-stiffness matrix of each step part. The specific expression is: ; in, represents the identity matrix; Indicates the The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the The inverse matrix of the unconstrained node sub-stiffness matrix of the work step part.

[0028] Based on the idea of ​​block matrix, through elementary row transformation, Convert to upper triangular matrix , the specific expression is: ; in, express The number of rows of the block matrix composed of non-zero elements in .

[0029] Accordingly, through the conversion The same elementary row transformation will Convert to , the specific expression is: ; ; ; ; in, Indicates the first The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the first -1 Inverse matrix of the non-constrained node stiffness matrix of the step part; express The block sub-matrix of ; express The block sub-matrix of ; Represents the row block index, ; Represents the column block index, ; The number of rows is , The number of rows and columns is based on Perform calculations; Indicates transpose.

[0030] Finally, through the conversion The opposite elementary row transformation will Transform again to get .

[0031] Step S4: Calculate the volume coefficient matrix of the parts in each process step. The specific expression is: ; in, Indicates the The sub-geometry matrix corresponding to the work step constraint node; Indicates the The sub-geometry matrix corresponding to the non-constrained nodes of the work step; Indicates the Volume coefficient matrix of process step parts.

[0032] Step S5: Assemble the volume coefficient matrices of the parts in each process step into a volume coefficient matrix. The specific expression is: ; in, represents the volume coefficient matrix.

[0033] The volume coefficient matrix was calculated using a computer with an Intel(R) Core(TM) i7-9700F processor and 16GB of memory. The traditional inversion method took 817.5 seconds to calculate the volume coefficient matrix, while the method of the present invention took only 433 seconds, with a 47% improvement in computational efficiency. This strongly demonstrates that the method of the present invention has significant advantages in improving the efficiency of volume coefficient matrix calculation.

[0034] Monitor the deformation force data of parts in real time during the processing of parts, and solve the residual stress based on the solved volume coefficient matrix The specific expression is: ; in, represents residual stress; Indicates deformation force.

[0035] The part is divided into 10 layers along the thickness direction, and residual stress (including residual stress in the x-direction and residual stress in the y-direction) is preset in each layer. The residual stress of the part is solved based on the method of the present invention, and the obtained residual stress result is compared with the preset residual stress. The comparison results are as follows: Figure 3 The mean absolute error of the residual stress inference in the x-direction is 4.3 MPa, and the mean absolute error of the residual stress inference in the y-direction is 2.6 MPa. These comparisons demonstrate the excellent prediction accuracy of the proposed method, verifying its effectiveness and reliability in residual stress field inference.

[0036] Therefore, the present invention adopts the above-mentioned method and application of rapid solution of volume coefficient matrix for residual stress field inference. Through an innovative matrix operation solution method, while ensuring the calculation accuracy, it significantly improves the calculation efficiency of the volume coefficient matrix, and effectively solves the problem of low calculation efficiency of the volume coefficient matrix in the existing technology.

[0037] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit the same. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that they can still modify or replace the technical solutions of the present invention with equivalents, and these modifications or equivalent replacements cannot cause the modified technical solutions to deviate from the spirit and scope of the technical solutions of the present invention.

Claims

1. A method for quickly solving the volume coefficient matrix for residual stress field inference, characterized in that: The following steps are involved: Step S1: Calculate the inverse matrix of the non-constrained node sub-stiffness matrix of the initial state of the part; Step S2: Calculate the change in the sub-stiffness matrix of the unconstrained nodes caused by the material removal of the parts in each process step; Step S3: Calculate the inverse matrix of the unconstrained node sub-stiffness matrix of each step part; Step S4: Calculate the volume coefficient matrix of the parts in each process step; Step S5: Assemble the volume coefficient matrices of the parts in each process step into a volume coefficient matrix.

2. A method for rapidly solving a volume coefficient matrix for residual stress field inference according to claim 1, characterized in that: In step S3, the specific expression of the inverse matrix of the unconstrained node sub-stiffness matrix of each step part is: ； ； in, Indicates the The overall stiffness matrix of the work step; Indicates the Unconstrained node sub-stiffness matrix of the work step parts; Indicates the Unconstrained-constrained node sub-stiffness matrix of the work step part; Indicates the Constrained-unconstrained node sub-stiffness matrix of the work step parts; Indicates the The constraint node sub-stiffness matrix of the work step parts; represents the identity matrix; Indicates the The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the The change in the stiffness matrix of the unconstrained nodes caused by the material removal of the work step parts; Indicates the number of work steps.

3. A method for rapidly solving the volume coefficient matrix for residual stress field inference according to claim 2, characterized in that: In step S3, Based on the idea of ​​block matrix, through elementary row transformation, Convert to upper triangular matrix , the specific expression is: ； in, express The number of rows of the block matrix composed of non-zero elements in ; Accordingly, through the conversion The same elementary row transformation will Convert to , the specific expression is: ； ； ； ； in, Indicates the first The inverse matrix of the non-constrained node stiffness matrix of the work step parts; Indicates the first -1 Inverse matrix of the non-constrained node stiffness matrix of the step part; express The block sub-matrix of ; express The block sub-matrix of ; Represents the row block index, ; Represents the column block index, ; The number of rows is , The number of rows and columns is based on Perform calculations; represents transpose; Finally, through the conversion The opposite elementary row transformation will Transform again to get .

4. A method for rapidly solving a volume coefficient matrix for residual stress field inference according to claim 3, characterized in that: In step S4, the specific expression of the volume coefficient matrix of the parts in each process step is: ； in, Indicates the The sub-geometry matrix corresponding to the work step constraint node; Indicates the The sub-geometry matrix corresponding to the non-constrained nodes of the work step; Indicates the Volume coefficient matrix of process step parts.

5. A method for rapidly solving a volume coefficient matrix for residual stress field inference according to claim 4, characterized in that: In step S5, the volume coefficient matrix of each step part is assembled into a specific expression of the volume coefficient matrix: ； in, represents the volume coefficient matrix.

6. An application of the method for quickly solving the volume coefficient matrix according to any one of claims 1 to 5 in residual stress field inference, characterized in that: The volume coefficient matrix is ​​used to solve the residual stress. The specific expression is: ； in, represents residual stress; Indicates deformation force.

Images