Multi-mode fretting wear calculation method, device and system and storage medium

Abstract

The invention discloses a multi-mode fretting wear calculation method, device and system and a storage medium. The multi-mode fretting wear calculation method comprises the steps that S1, dynamic node data is obtained, and the wear amount is calculated; s2, calculating a geometric self-adaptive wear direction; s3, performing node classification and index optimization and adaptive adjustment of an increment step; and S4, carrying out million-level node parallel computing by adopting a block storage and matrix compression technology, and controlling the grid scanning times through KMESHSWEEP parameters. By adopting the technical scheme of the invention, the technical problems of poor geometric adaptability, low algorithm convergence and limited calculation scale in fretting wear simulation are solved.

Description

Technical Field

[0001] The present invention belongs to the field of fretting wear technology, and in particular relates to a multi-mode fretting wear calculation method and device, system, and storage medium. Background Art

[0002] Fretting wear is a surface damage phenomenon caused by the reciprocating motion of mechanical contact parts with small amplitudes. It is common in key parts such as fastening connections and bearing fits in the fields of aerospace, rail transportation, precision machinery, etc. Long-term fretting can lead to material delamination, crack initiation and expansion, significantly reducing the service life of the structure. Traditional experimental research methods have long cycles and high costs, and it is difficult to observe the wear evolution process in real time. Numerical simulation technology provides an effective means to predict wear morphology and life assessment. However, existing finite element simulation methods face significant challenges in simulating fretting wear: first, the dynamic changes in the wear direction under complex contact conditions are difficult to describe accurately, especially the local stress concentration effects in geometric mutation areas such as corners and edges are often simplified; second, there is a lack of efficient algorithms for the coupling mechanism of slip increment and contact pressure, resulting in poor convergence of multi-physics field iterations; third, the existing user subroutines are insufficient in adaptive mesh update capabilities for large-scale node wear, making it difficult to balance calculation accuracy and efficiency. Summary of the Invention

[0003] The technical problem to be solved by the present invention is to provide a multi-mode fretting wear calculation method and device, system and storage medium.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] A multi-mode fretting wear calculation method, comprising:

[0006] Step S1, obtaining dynamic node data and calculating wear amount;

[0007] Step S2, calculating the geometrically adaptive wear direction;

[0008] Step S3: node classification and index optimization and incremental step adaptive adjustment;

[0009] Step S4: Use block storage and matrix compression technology to perform parallel computing on millions of nodes, and control the number of grid scans through the KMESHSWEEP parameter.

[0010] Preferably, in step S1, the node contact pressure, slip increment and spatial coordinates are obtained in real time by calling the built-in subroutines GETVRMAVGATNODE and GETVRN, and the local instantaneous wear amount is calculated in combination with the improved Archard wear model.

[0011] Preferably, in step S1, the Fortran common storage area is used to define global arrays OLDSLIPZ and TEMPSLIPZ to store the node historical slip data.

[0012] 4. The multi-mode micro-wear calculation method according to claim 3 is characterized in that in step S1 and step S2, the coordinate difference between the node and the adjacent master node in the three-dimensional space is calculated in real time by introducing the master-slave node relative displacement method, the wear direction of the corner area is determined by combining the vector decomposition technology, and the local coordinate system transformation is used to achieve accurate mapping of the global wear amount to the local deformation direction.

[0013] Preferably, in step S3, key contact elements are screened based on the node number range; the time step is dynamically controlled by the PNEWDT parameter, and when the local wear amount changes suddenly or the iteration does not converge, the step size is automatically reduced and the calculation is restarted.

[0014] The present invention also provides a multi-mode fretting wear calculation device, characterized by comprising:

[0015] The first processing module is used to obtain dynamic node data and calculate wear amount;

[0016] A second processing module is used to calculate the geometric adaptive wear direction;

[0017] The third processing module is used for node classification and index optimization and incremental step adaptive adjustment;

[0018] The fourth processing module is used to perform parallel computing on millions of nodes using block storage and matrix compression technology, while controlling the number of grid scans through the KMESHSWEEP parameter.

[0019] The present invention also provides a multi-mode fretting wear calculation system, comprising: a memory and a processor, wherein the memory stores a computer program executed by the processor, and the computer program executes a multi-mode fretting wear calculation method when executed by the processor.

[0020] The present invention also provides a storage medium, on which a computer program is stored. When the computer program is run, the multi-mode fretting wear calculation method is executed.

[0021] This invention effectively solves the technical problems of poor geometric adaptability, low algorithm convergence, and limited computational scale in fretting wear simulation, providing a high-precision and high-efficiency solution for wear life prediction of complex mechanical systems. It has the following technical effects:

[0022] 1. The present invention introduces the relative displacement method of master and slave nodes to calculate the coordinate difference between a node and an adjacent master node in three-dimensional space in real time. It combines vector decomposition technology to accurately determine the wear direction in the corner area, and uses local coordinate system transformation (through the direction cosine matrix ALOCAL) to achieve accurate mapping of global wear to local deformation direction, significantly improving the calculation accuracy of wear vectors under complex geometric boundary conditions.

[0023] 2. Based on the improved Archard wear model, the present invention dynamically obtains the node-level contact pressure (CPRESS) and slip increment (CSLIP) by calling the GETVRMAVGATNODE subroutine, and uses the common storage area (COMMON block) to efficiently store historical slip data (OLDSLIPZ, TEMPSLIPZ) to avoid repeated access to hard disk data; at the same time, combined with the incremental step convergence control parameter (PNEWDT), the time step is dynamically adjusted to ensure the stability and computational efficiency of multi-physics field iteration.

[0024] 3. The present invention proposes a node classification strategy: three-dimensional vector updates are used for corner area nodes (such as NODE 333100-333200), while the wear of internal nodes of the mating surface is only accumulated along the normal direction (ULOCAL(NDIM)), which significantly reduces the calculation dimension; at the same time, the JELEMLIST index and JELEMTYPE screening mechanism are optimized to reduce redundant unit traversal, realize efficient grid update of millions of nodes, and provide reliable support for engineering-level micro-wear simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0026] Figure 1 This is a flow chart of a multi-mode fretting wear calculation method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0027] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0028] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0029] Example 1:

[0030] like Figure 1 As shown, an embodiment of the present invention provides a multi-mode fretting wear calculation method, comprising:

[0031] Step S1, obtaining dynamic node data and calculating wear amount;

[0032] Step S2, calculating the geometrically adaptive wear direction;

[0033] Step S3: node classification and index optimization and incremental step adaptive adjustment;

[0034] Step S4: Use block storage and matrix compression technology to perform parallel computing on millions of nodes, and control the number of grid scans through the KMESHSWEEP parameter.

[0035] As an implementation method of an embodiment of the present invention, in step S1, by calling the built-in subroutines GETVRMAVGATNODE and GETVRN, the node contact pressure (CPRESS), slip increment (CSLIP) and spatial coordinate (COORD) are obtained in real time, and the local instantaneous wear volume (SURFV) is calculated in combination with the improved Archard wear model. The expression is:

[0036] SURFV=K*CPRESS*ΔCSLIP

[0037] Where K is the material wear coefficient, and ΔCSLIP is dynamically updated through the historical slip arrays (OLDSLIPZ and TEMPSLIPZ).

[0038] In step S1, the Fortran common storage area (COMMON block) is used to define global arrays OLDSLIPZ and TEMPSLIPZ to store node historical slip data, avoiding frequent reading and writing of the hard disk and improving data access efficiency.

[0039] As an implementation of an embodiment of the present invention, in step S2, for the geometric mutation area (such as NODE333100-333200), the master-slave node relative displacement method is used to correct the three-dimensional vector of the corner area, specifically:

[0040] Define the master node (Master = NODE + 1000) and obtain its coordinates (cmaster_x, cmaster_y, cmaster_z);

[0041] The space vector difference between the calculation node and the master node is combined with the wear volume SURFV to generate the wear vector in the global coordinate system (WVGLOBAL);

[0042] The global wear vector is mapped to the local coordinate system (WVLOCAL) through the direction cosine matrix ALOCAL, and the node displacement ULOCAL is updated to achieve three-dimensional direction adaptive correction.

[0043] For nodes in non-corner areas, the wear amount is directly accumulated along the normal direction of the local coordinate system (ULOCAL(NDIM)) to reduce the calculation dimension and improve efficiency.

[0044] As an implementation method of an embodiment of the present invention, in step S3, key contact elements are screened based on the node number range (such as JELEMLIST and JELEMTYPE) to reduce unnecessary element traversal; the maximum number of elements is predefined in combination with the NELEMMAX parameter to avoid performance loss caused by dynamic memory allocation.

[0045] The time step is dynamically controlled by the PNEWDT parameter. When the local wear amount changes suddenly or the iteration does not converge, the step size is automatically reduced and the calculation is restarted to ensure numerical stability.

[0046] As an implementation method of an embodiment of the present invention, in step S4, block storage (JGVBLOCK) and matrix compression technology are used to support efficient parallel computing of millions of nodes (such as DIMENSION OLDSLIPZ (500000)), and the number of grid scans is controlled by the KMESHSWEEP parameter to balance the calculation accuracy and resource consumption.

[0047] The present invention dynamically obtains node contact pressure, slip increment and coordinate information, and combines it with the Archard wear model to achieve real-time calculation of local wear. A node classification strategy is innovatively introduced, and the master-slave node relative displacement method is used for the corner area to accurately calculate the three-dimensional wear vector, and the global wear is mapped to the local coordinate system through coordinate transformation; for the internal nodes of the mating surface, the calculation is simplified based on the normal cumulative slip, which significantly improves the adaptability to complex geometric boundary conditions. In addition, the access efficiency of historical variables is optimized through the common storage area (COMMON block), and combined with the incremental step convergence control (PNEWDT parameter), the numerical stability problem in the large-scale grid update process is effectively solved, providing an efficient and reliable simulation tool for the quantitative prediction and life optimization of micro-wear in engineering.

[0048] Example 2:

[0049] An embodiment of the present invention further provides a multi-mode fretting wear calculation device, comprising:

[0050] The first processing module is used to obtain dynamic node data and calculate wear amount;

[0051] A second processing module is used to calculate the geometric adaptive wear direction;

[0052] The third processing module is used for node classification and index optimization and incremental step adaptive adjustment;

[0053] The fourth processing module is used to perform parallel computing on millions of nodes using block storage and matrix compression technology, while controlling the number of grid scans through the KMESHSWEEP parameter.

[0054] Example 3:

[0055] An embodiment of the present invention further provides a multi-mode fretting wear calculation system, comprising: a memory and a processor, wherein the memory stores a computer program executed by the processor, and the computer program executes a multi-mode fretting wear calculation method when executed by the processor.

[0056] Example 4:

[0057] An embodiment of the present invention further provides a storage medium having a computer program stored thereon. The computer program executes the multi-mode fretting wear calculation method when running.

[0058] The embodiments described above are merely descriptions of preferred embodiments of the present invention and are not intended to limit the scope of the present invention. Without departing from the spirit of the present invention, various modifications and improvements made to the technical solutions of the present invention by persons skilled in the art should fall within the scope of protection defined by the claims of the present invention.

Claims

1. A multi-mode fretting wear calculation method, characterized in that: include: Step S1, obtaining dynamic node data and calculating wear amount; Step S2, calculating the geometric adaptive wear direction; Step S3: node classification and index optimization and incremental step adaptive adjustment; Step S4: Use block storage and matrix compression technology to perform parallel computing on millions of nodes, and control the number of grid scans through the KMESHSWEEP parameter.

2. The multi-mode fretting wear calculation method according to claim 1, wherein: In step S1, the built-in subroutines GETVRMAVGATNODE and GETVRN are called to obtain the node contact pressure, slip increment and spatial coordinates in real time, and the local instantaneous wear amount is calculated in combination with the improved Archard wear model.

3. The multi-mode fretting wear calculation method according to claim 2, wherein: In step S1, the Fortran common storage area is used to define global arrays OLDSLIPZ and TEMPSLIPZ to store the node historical slip data.

4. The multi-mode fretting wear calculation method according to claim 3, wherein: In step S1 and step S2, by introducing the relative displacement method of master and slave nodes, the coordinate difference between the node and the adjacent master node in the three-dimensional space is calculated in real time, the wear direction of the corner area is determined by combining the vector decomposition technology, and the local coordinate system transformation is used to realize the accurate mapping of the global wear amount to the local deformation direction.

5. The multi-mode fretting wear calculation method according to claim 4, wherein: In step S3, key contact elements are selected based on the node number range; the time step is dynamically controlled by the PNEWDT parameter. When the local wear amount changes suddenly or the iteration does not converge, the step size is automatically reduced and the calculation is restarted.

6. A multi-mode fretting wear calculation device, characterized in that: include: The first processing module is used to obtain dynamic node data and calculate wear amount; A second processing module is used to calculate the geometric adaptive wear direction; The third processing module is used for node classification and index optimization and incremental step adaptive adjustment; The fourth processing module is used to perform parallel computing on millions of nodes using block storage and matrix compression technology, while controlling the number of grid scans through the KMESHSWEEP parameter.

7. A multi-mode fretting wear calculation system, characterized in that: include: A memory and a processor, wherein the memory stores a computer program executed by the processor, and when the computer program is executed by the processor, the multi-mode fretting wear calculation method according to any one of claims 1 to 5 is executed.

8. A storage medium, characterized in that: The storage medium stores a computer program, which, when running, executes the multi-mode fretting wear calculation method according to any one of claims 1 to 5.

Images