Wall surface wear simulation method and device independent of grid size, equipment and medium

Abstract

The present application relates to the technical field of numerical simulation, and particularly relates to a wall wear simulation method and device independent of grid size, equipment and medium, comprising: obtaining a model file of a structure and a background grid independent of a facet element size of the model file; obtaining, through simulation, a collision contact point between a granular particle and a wall of the structure and a wall wear volume generated by the collision contact point due to collision contact; accumulating, according to the wall wear volume generated by the collision contact point due to collision contact, a wall wear volume of a first grid node closest to the collision contact point in the background grid; and performing wear diffusion on the wall of the structure according to a latest accumulated wall wear volume of a second grid node corresponding to a sampling point in the background grid, to obtain a distribution cloud picture of the wear depth on the wall grid. Through the present application, wall wear simulation independent of the facet element size of the model file can be performed, and the accuracy and engineering applicability of wear prediction are improved.

Description

Technical Field

[0001] This invention relates to the field of numerical simulation technology, and more specifically, to a method, apparatus, equipment, and medium for simulating wall wear independent of mesh size. Background Technology

[0002] Numerical simulation is an important tool for studying the interaction between granular particles and structures, allowing us to obtain physical quantities that are impossible to obtain experimentally. For simulating wall wear, the Discrete Element Method (DEM) is one of the most important numerical methods. DEM uses spherical or non-spherical discrete elements to represent granular particles and employs mesh-based CAD model files (such as STL / OBJ format) to describe the geometry of the structure. The mesh elements in the model file themselves have specific dimensions.

[0003] When using the traditional DEM method for wear prediction, the wear-causing surface unit is usually located first based on the collision contact point between the particles and the structure wall. Then, the wear volume generated by a single collision is divided by the area of ​​the surface unit to obtain the average wear depth of the surface unit.

[0004] However, traditional methods have significant limitations. The wear volume generated by a collision is only recorded on a single patch element at the point of impact. When the patch element size of the structure model changes, even if the total wear volume remains constant, the calculated average wear depth will also change, leading to simulation results that do not match the actual physical wear distribution. Summary of the Invention

[0005] In view of this, the purpose of this application is to provide a method, apparatus, device and medium for simulating wall wear that is independent of mesh size. This method can simulate wall wear in a way that is independent of the size of the surface elements in the solid wall boundary model file, thereby eliminating the influence of the surface size of the structure on the simulation results and improving the accuracy and engineering applicability of wear prediction.

[0006] In a first aspect, embodiments of this application provide a method for simulating wall wear that is independent of mesh size, the method comprising:

[0007] Obtain the solid wall boundary model file of the structure, and a background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch elements in the solid wall boundary model file;

[0008] The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, as well as the wall wear volume caused by the collision contact point.

[0009] Based on the wall wear volume generated by the collision at the collision point, the wall wear volume of the first grid node in the background grid closest to the collision point is accumulated.

[0010] Continue with the collision simulation at the next simulation time step; and based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, perform wear diffusion on the wall of the structure to obtain a distribution cloud map of the wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

[0011] In one possible implementation, a background mesh surrounding the solid wall boundary of the structure, independent of the patch element size of the solid wall boundary model file, is obtained through the following steps:

[0012] The node spacing of the background mesh is calculated based on the preset minimum size of the granular particles.

[0013] Based on the geometric dimensions of the solid wall boundary model file, determine the minimum geometric bounding box of the solid wall boundary model;

[0014] The background mesh is obtained by dividing the interior of the geometric bounding box into background meshes based on the node spacing.

[0015] In one possible implementation, the discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the wall wear volume caused by the collision contact point, including:

[0016] The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the sliding distance of the granular particles on the wall surface of the structure, as well as the interaction force between the granular particles and the structure.

[0017] The volume of wall wear caused by the collision at the point of impact is calculated based on the sliding distance of the granular particles on the wall surface of the structure and the interaction force between the granular particles and the structure.

[0018] In one possible implementation, the first mesh node closest to the collision point is determined using the following formula:

[0019] ;

[0020] in, The coordinates of the first grid node are... The X-axis coordinate of the first grid node is... The Y-axis coordinate of the first grid node. The Z-axis coordinate of the first grid node is... This indicates that the result is rounded to the nearest integer. The X-axis coordinates of the collision contact points. The Y-axis coordinate of the collision contact point. The Z-axis coordinate of the collision contact point. The starting point of the background grid X-axis coordinates The starting point of the background grid Y-axis coordinate, The starting point of the background grid Z-axis coordinate, This represents the node spacing corresponding to the background grid.

[0021] In one possible implementation, the step of performing wear diffusion on the wall surface of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid to obtain a distribution cloud map of the wear depth on the wall grid includes:

[0022] Determine whether the current simulation time step is the preset output time step for the wall wear simulation;

[0023] If so, then based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, wear diffusion is carried out on the wall of the structure to obtain a distribution cloud map of wear depth on the wall grid.

[0024] In one possible implementation, the step of performing wear diffusion on the wall surface of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid to obtain a distribution cloud map of the wear depth on the wall grid includes:

[0025] For each second grid node, the wear impact value of the second grid node on the corresponding sampling point is calculated based on the distance between the second grid node and the corresponding sampling point;

[0026] Using the wear impact value of each second grid node on the corresponding sampling point as a weight, the latest accumulated wall wear volume of all second grid nodes is weighted and summed to obtain the wear depth at the sampling point in the distribution cloud map.

[0027] In one possible implementation, calculating the wear impact value of the second grid node on the corresponding sampling point based on the distance between the second grid node and the corresponding sampling point includes:

[0028] Substitute the distance between the second grid node and the corresponding sampling point into the diffusion function to calculate the wear influence of the second grid node on the corresponding sampling point.

[0029] The diffusion function satisfies the following condition:

[0030] ;

[0031] ;

[0032] ;

[0033] in, This represents the wear impact value when the distance between the second grid node and the corresponding sampling point is 0. This is the preset maximum wear impact value; The distance between the second grid node and the corresponding sampling point is The wear impact value at that time; R is the preset wear diffusion radius; This represents the distance between the second grid node and the corresponding sampling point. The azimuth angle of the second grid node relative to the corresponding sampling point.

[0034] Secondly, embodiments of this application also provide a wall wear simulation device independent of mesh size, the device comprising:

[0035] The acquisition module is used to acquire the solid wall boundary model file of the structure, and the background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch unit of the solid wall boundary model file.

[0036] The calculation module is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file using the discrete element method, to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, and the wall wear volume caused by the collision contact point.

[0037] The calculation module is also used to accumulate the wall wear volume of the first grid node in the background grid that is closest to the collision contact point based on the wall wear volume generated by the collision contact point in this collision contact.

[0038] The calculation module is also used to continue the collision simulation for the next simulation time step; and to perform wear diffusion on the wall of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, so as to obtain the distribution cloud map of wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

[0039] In one possible implementation, the acquisition module is specifically configured to acquire a background mesh surrounding the solid wall boundary of the structure that is independent of the patch element size of the solid wall boundary model file through the following steps:

[0040] The node spacing of the background mesh is calculated based on the preset minimum size of the granular particles.

[0041] Based on the geometric dimensions of the solid wall boundary model file, determine the minimum geometric bounding box of the solid wall boundary model;

[0042] The background mesh is obtained by dividing the interior of the geometric bounding box into background meshes based on the node spacing.

[0043] In one possible implementation, the calculation module is specifically used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file using the discrete element method, to obtain the sliding distance of the loose particles on the wall surface of the structure, and the interaction force between the loose particles and the structure; based on the sliding distance of the loose particles on the wall surface of the structure and the interaction force between the loose particles and the structure, the module calculates the wall wear volume generated by the collision contact point due to the collision contact.

[0044] In one possible implementation, the calculation module is specifically used to determine the first mesh node closest to the collision contact point using the following formula:

[0045] ;

[0046] in, The coordinates of the first grid node are... The X-axis coordinate of the first grid node is... The Y-axis coordinate of the first grid node. The Z-axis coordinate of the first grid node is... This indicates that the result is rounded to the nearest integer. The X-axis coordinates of the collision contact points. The Y-axis coordinate of the collision contact point. The Z-axis coordinate of the collision contact point. The starting point of the background grid X-axis coordinates The starting point of the background grid Y-axis coordinate, The starting point of the background grid Z-axis coordinate, This represents the node spacing corresponding to the background grid.

[0047] In one possible implementation, the calculation module is specifically used to determine whether the current simulation time step is the preset output time step of the wall wear simulation; if so, then based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, wear diffusion is performed on the wall of the structure to obtain a distribution cloud map of the wear depth on the wall grid.

[0048] In one possible implementation, the calculation module is specifically used to calculate the wear influence value of the second grid node on the corresponding sampling point based on the distance between the second grid node and the corresponding sampling point for each second grid node; and to use the wear influence value of each second grid node on the corresponding sampling point as a weight to perform a weighted summation of the latest accumulated wall wear volume of all second grid nodes to obtain the wear depth at the sampling point in the distribution cloud map.

[0049] In one possible implementation, the calculation module is specifically used to substitute the distance between the second grid node and the corresponding sampling point into the diffusion function to calculate the wear influence value of the second grid node on the corresponding sampling point.

[0050] The diffusion function satisfies the following condition:

[0051] ;

[0052] ;

[0053] ;

[0054] in, This represents the wear impact value when the distance between the second grid node and the corresponding sampling point is 0. This is the preset maximum wear impact value; The distance between the second grid node and the corresponding sampling point is The wear impact value at that time; R is the preset wear diffusion radius; This represents the distance between the second grid node and the corresponding sampling point. The azimuth angle of the second grid node relative to the corresponding sampling point.

[0055] Thirdly, embodiments of this application also provide an electronic device, including: a processor, a storage medium, and a bus, wherein the storage medium stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor communicates with the storage medium via the bus, and the processor executes the machine-readable instructions to perform the steps of the wall wear simulation method independent of mesh size as described in any of the first aspects.

[0056] Fourthly, embodiments of this application also provide a computer-readable storage medium storing a computer program, which, when executed by a processor, performs the steps of the wall wear simulation method independent of mesh size as described in any of the first aspects.

[0057] This application provides a method, apparatus, device, and medium for simulating wall wear independent of mesh size. The method includes: acquiring a solid wall boundary model file of a structure and a background mesh of the structure independent of the patch element size of the model file; performing collision simulation on the solid wall boundary model to obtain the collision contact points between the granular particles and the wall of the structure at the current simulation time step, and the wall wear volume generated by the collision contact points; accumulating the wall wear volume of the first mesh node in the background mesh closest to the collision contact point based on the wall wear volume generated by the collision contact points; continuing the collision simulation; and performing wear diffusion on the wall of the structure based on the latest accumulated wall wear volume of the second mesh node corresponding to the sampling point in the background mesh to obtain a distribution cloud map of the wear depth on the wall mesh. This application provides a way to simulate wall wear independent of the patch element size of the model file, improving the accuracy and engineering applicability of wear prediction. Attached Figure Description

[0058] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0059] Figure 1 A schematic flowchart of a wall wear simulation method independent of mesh size provided in an embodiment of this application is shown;

[0060] Figure 2 This illustration shows a schematic diagram of the collision process between granular particles and the wall surface of a structure, as provided in an embodiment of this application.

[0061] Figure 3 This illustration shows a schematic diagram of the diffusion function provided in this application embodiment when the preset wear diffusion radius is 2;

[0062] Figure 4 This illustration shows a schematic diagram of wear diffusion at sampling points on the solid boundary of a structure, provided in an embodiment of this application.

[0063] Figure 5This diagram illustrates a comparison of simulation results between the traditional wall wear simulation method provided in this application and the wall wear simulation method provided in this application.

[0064] Figure 6 A schematic diagram of a wall wear simulation device independent of mesh size provided in an embodiment of this application is shown;

[0065] Figure 7 A schematic diagram of the structure of an electronic device provided in an embodiment of this application is shown. Detailed Implementation

[0066] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. It should be understood that the accompanying drawings in this application are for illustrative and descriptive purposes only and are not intended to limit the scope of protection of this application. Furthermore, it should be understood that the schematic drawings are not drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of this application. It should be understood that the operations in the flowcharts may not be implemented in sequence, and steps without logical contextual relationships may be reversed or implemented simultaneously. In addition, those skilled in the art, guided by the content of this application, may add one or more other operations to the flowcharts, or remove one or more operations from the flowcharts.

[0067] Furthermore, the described embodiments are merely some, not all, of the embodiments of this application. The components of the embodiments of this application described and illustrated herein can typically be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0068] To enable those skilled in the art to utilize the content of this application, and in conjunction with the specific application scenario of "numerical simulation technology," the following implementation methods are provided. For those skilled in the art, the general principles defined herein can be applied to other embodiments and application scenarios without departing from the spirit and scope of this application. Although this application is primarily described within the "numerical simulation technology field," it should be understood that this is merely an exemplary embodiment.

[0069] It should be noted that the term "comprising" will be used in the embodiments of this application to indicate the presence of the features declared thereafter, but does not exclude the addition of other features.

[0070] The following is a detailed description of a wall wear simulation method independent of mesh size provided in the embodiments of this application.

[0071] Reference Figure 1 The diagram shown is a flowchart illustrating a wall wear simulation method independent of mesh size provided in an embodiment of this application. The exemplary steps of this embodiment are described below:

[0072] S101. Obtain the solid wall boundary model file of the structure, and the background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch elements in the solid wall boundary model file.

[0073] In this embodiment, the surface of the structure is generally made of metal, such as an excavator or a water turbine (a type of power machinery that converts the energy (kinetic or potential energy) of water flow into rotational mechanical energy). The solid boundary of the structure refers to the structural wall surface that is eroded, collided with, or worn by loose particles (such as stones). The solid boundary model file of the structure is a general CAD drawing file such as .stl or .obj, which is also the geometric model file of the structure. It is used to characterize the solid wall boundary of the structure in subsequent discrete element simulations and is the object of collision and wear of loose particles. The background mesh is obtained by dividing the solid boundary model in the solid boundary model file into a background mesh. Specifically, the background mesh surrounding the solid boundary of the structure, which is independent of the size of the surface element in the solid boundary model file, is obtained through the following steps:

[0074] Step 1: Calculate the node spacing of the background mesh based on the preset minimum size of the granular particles.

[0075] In this embodiment, the node spacing of the background mesh is independent of the geometric dimensions of the structure model file. The node spacing can be set. ;in, This is the preset size factor for the background mesh (a positive value below 1.0 can be used to ensure calculation accuracy). The preset minimum size of the granular particles for collision simulation with structures.

[0076] Step 2: Determine the minimum geometric bounding box of the solid wall boundary model based on the geometric dimensions of the solid wall boundary model file.

[0077] In this embodiment, the solid wall boundary model file of the structure is imported, and the minimum geometric bounding box of the solid wall boundary model is defined according to the coordinate value range of the model file to ensure that the solid wall boundary of the structure is surrounded by the nodes of the background mesh.

[0078] Step 3: Divide the interior of the geometric bounding box into a background mesh according to the node spacing to obtain the background mesh.

[0079] S102. The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, as well as the wall wear volume caused by the collision contact point.

[0080] In the embodiments of this application, reference is made to Figure 2 The diagram illustrates the collision process between granular particles and the wall of a structure, as provided in this embodiment. The granular particles come into contact with the wall (i.e., the solid boundary), and the red dots represent the collision contact points. A discrete element method (DEM) is used to simulate the collision of the solid boundary model in the solid boundary model file. The interaction between the granular particles and the solid boundary (i.e., the wall) of the structure during their movement is calculated, and the coordinates of the collision contact points are recorded. ;in, The X-axis coordinates of the collision contact points. The Y-axis coordinate of the collision contact point. The Z-axis coordinate of the collision contact point.

[0081] The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the wall wear volume caused by the collision contact point, including:

[0082] Step 1: Use the discrete element method to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the sliding distance of the granular particles on the wall surface of the structure, as well as the interaction force between the granular particles and the structure.

[0083] Step 2: Based on the sliding distance of the granular particles on the wall surface of the structure and the interaction force between the granular particles and the structure, calculate the volume of wall wear caused by the collision at the point of contact.

[0084] In this embodiment, the wear volume generated by the collision is calculated based on the wear model used in the simulation. The wear model describes the cutting and erosion effects on the wall surface when granular particles interact with it, resulting in a certain degree of volume loss on the structure's wall surface. The sliding distance of the granular particles on the structure's wall surface can be substituted into the wear model below to calculate the wall wear volume generated by the collision contact point:

[0085] ;

[0086] in, The volume of wall wear caused by the collision at the point of contact. The preset wear coefficient; The interaction force between the granular particles and the structure is the force of their collision. The sliding distance of the granular particles on the wall surface of the structure; The material hardness of the structure is given. The wear coefficient k and the material hardness H are known coefficients, while the interaction force F and the sliding distance s can be calculated through DEM simulation. Therefore, a definite wear volume can be obtained at each simulation time step.

[0087] S103. Based on the wall wear volume generated by the collision at the point of contact, accumulate the wall wear volume of the first grid node in the background grid that is closest to the point of contact.

[0088] In this embodiment of the application, firstly, based on the coordinates of the collision contact point obtained above and the known starting point of the background mesh... The coordinates and the node spacing of the background mesh The first mesh node closest to the collision point is determined. Specifically, the first mesh node closest to the collision point (i.e., ...) is determined using the following formula. Figure 2 (The nearest grid node marked with an "X")

[0089] ;

[0090] in, The coordinates of the first grid node are... The X-axis coordinate of the first grid node is... The Y-axis coordinate of the first grid node. The Z-axis coordinate of the first grid node is... This indicates that the result is rounded to the nearest integer. The X-axis coordinates of the collision contact points. The Y-axis coordinate of the collision contact point. The Z-axis coordinate of the collision contact point. The starting point of the background grid X-axis coordinates The starting point of the background grid Y-axis coordinate, The starting point of the background grid Z-axis coordinate, This represents the node spacing corresponding to the background mesh. There is a certain difference between the coordinates of the first mesh node and the coordinates of the collision contact point, but this difference is negligible when dx is sufficiently small.

[0091] Then, based on the wall wear volume generated by the collision at the contact point, the wall wear volume of the first mesh node is calculated. Accumulate:

[0092] .

[0093] S104. Continue the collision simulation for the next simulation time step; and based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, perform wear diffusion on the wall of the structure to obtain the distribution cloud map of the wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

[0094] In this embodiment, the DEM simulation time continues to advance. When the collision simulation reaches a certain time step, the calculation results need to be saved and visualized; this time step is called the output time step. Therefore, it is determined whether the current simulation time step is the preset output time step for the wall wear simulation. If so, wear diffusion is performed on the wall of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, resulting in a wear depth distribution cloud map on the wall grid. Each output time step outputs a distribution cloud map, and each distribution cloud map contains the wear depth of at least one sampling point.

[0095] The wall mesh is determined by the solid wall boundary model file. The wall mesh is equal to the patch element of the solid wall boundary model file and is related to the patch element size.

[0096] Here, since wear volume and wear depth have a linear superposition characteristic, in the optional embodiment provided in this application, the diffusion of wear volume is not processed at each simulation time step, but only at the output time step, thereby significantly improving computational efficiency; specifically, based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, wear diffusion is performed on the wall of the structure to obtain a distribution cloud map of the wear depth on the wall grid, including:

[0097] Step 1: For each second grid node, calculate the wear impact value of the second grid node on the corresponding sampling point based on the distance between the second grid node and the corresponding sampling point.

[0098] In this embodiment, the distance between the second grid node and the corresponding sampling point is substituted into the diffusion function to calculate the wear influence of the second grid node on the corresponding sampling point. The function value of the diffusion function characterizes the attenuation weight of the wear effect generated by the collision between the granular particles and the wall of the structure as it diffuses from the collision point along the wall towards the surrounding area.

[0099] Here, the diffusion function is used to post-process the wear calculation results. Physically, when granular particles collide with the wall of a structure, the wear depth should be greatest at the collision point, and then gradually decrease towards the surrounding area. Therefore, the diffusion function satisfies the following condition:

[0100] ;

[0101] ;

[0102] ;

[0103] in, This represents the wear impact value when the distance between the second grid node and the corresponding sampling point is 0. This is the preset maximum wear impact value; The distance between the second grid node and the corresponding sampling point is The wear impact value at that time; R is the preset wear diffusion radius; This represents the distance between the second grid node and the corresponding sampling point. The azimuth angle of the second grid node relative to the corresponding sampling point.

[0104] In the diffusion function condition, the first equation indicates that the diffusion function at... The maximum value is obtained when = 0, and the second expression represents... = Its impact is 0 at that time, and It must be monotonically decreasing. The third equation indicates that the integral of the wear depth along the plane equals the volume of wear, satisfying volume conservation. In principle, any function satisfying the above equation... All of these can be used as diffusion functions. This embodiment provides one of the most commonly used diffusion functions:

[0105] ;

[0106] Reference Figure 3 The diagram shown is a schematic of the diffusion function provided in this application embodiment when the preset wear diffusion radius is 2. Figure 4 This is a schematic diagram illustrating wear diffusion at sampling points on the solid boundary of a structure, provided in an embodiment of this application; from Figure 3 and Figure 4 As can be seen, the diffusion function has the largest value at the collision point and gradually decreases as the diffusion radius increases. When outputting the post-processing results, it is only necessary to know the coordinates p of a sampling point on the solid boundary of the structure, and then iterate through all second grid nodes within a distance R around that sampling point (such as...). Figure 4 As shown in the figure, the distance r between the sampling point and the second grid node is calculated, and the wear influence value of the second grid node on the corresponding sampling point is calculated.

[0107] Step 2: Using the wear impact value of each second grid node on the corresponding sampling point as a weight, perform a weighted summation of the latest accumulated wall wear volume of all second grid nodes to obtain the wear depth at the sampling point in the distribution cloud map.

[0108] In this embodiment, the wear influence value of each second grid node on the corresponding sampling point is used as a weight by the following formula, and the latest accumulated wall wear volume of all second grid nodes is weighted and summed to obtain the wear depth at the sampling point in the distribution cloud map:

[0109] ;

[0110] in, The wear depth at the sampling point. The number of nodes in the second grid. The function value is the degree of wear impact of the second grid node on the corresponding sampling point. The distance between the m-th second grid node and its corresponding sampling point. Let m be the coordinates of the second grid node. Let x be the x-coordinate of the m-th second grid node. Let be the Y-axis coordinate of the m-th second grid node. Let be the Z-axis coordinate of the m-th second grid node.

[0111] Here, the selection of sampling points can be arbitrary in principle, but in general, the selected sampling points are the center coordinates or node coordinates of the mesh patch elements in the model file, so that the distribution cloud map of wear depth can be obtained through simple visualization tools.

[0112] Reference Figure 5 The diagram shows a comparison of simulation results between the traditional wall wear simulation method provided in this application and the wall wear simulation method provided in this application. It can be seen that the traditional method can only record collisions at the surface element where the collision point is located. When the size of the surface element of the structure changes, although the wear volume remains the same, the average wear depth will change accordingly. In contrast, the method of this application is independent of the mesh surface size of the structure's geometry file. It only requires setting the wear diffusion radius R to obtain a wear depth independent of size, which is more consistent with the real physical process and has better simulation prediction capabilities.

[0113] Based on the same inventive concept, this application also provides a wall wear simulation device that is independent of grid size, corresponding to the wall wear simulation method that is independent of grid size. Since the principle of the device in this application is similar to the wall wear simulation method that is independent of grid size described above in this application, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be described again.

[0114] Reference Figure 6The diagram shown is a schematic of a wall wear simulation device independent of mesh size provided in an embodiment of this application. The device includes:

[0115] The acquisition module 601 is used to acquire the solid wall boundary model file of the structure, and the background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch unit of the solid wall boundary model file.

[0116] The calculation module 602 is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file using the discrete element method, to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, and the wall wear volume caused by the collision contact point.

[0117] The calculation module 602 is further configured to accumulate the wall wear volume of the first grid node in the background grid that is closest to the collision contact point based on the wall wear volume generated by the collision contact point in this collision.

[0118] The calculation module 602 is also used to continue the collision simulation in the next simulation time step; and to perform wear diffusion on the wall of the structure according to the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, so as to obtain the distribution cloud map of wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

[0119] This application provides a wall wear simulation device that is independent of mesh size. This device can simulate wall wear in a way that is independent of the size of the patch cells in the model file, thereby improving the accuracy of wear prediction and its engineering applicability.

[0120] like Figure 7 As shown in the embodiment of this application, an electronic device 700 includes a processor 701, a memory 702, and a bus. The memory 702 stores machine-readable instructions that can be executed by the processor 701. When the electronic device is running, the processor 701 communicates with the memory 702 via the bus. The processor 701 executes the machine-readable instructions to perform the steps of the wall wear simulation method that is independent of the grid size described above.

[0121] Specifically, the memory 702 and processor 701 mentioned above can be general-purpose memory and processor, without any specific limitations. When the processor 701 runs the computer program stored in the memory 702, it can execute the above-mentioned wall wear simulation method that is independent of the mesh size.

[0122] Corresponding to the above-described wall wear simulation method independent of mesh size, this application embodiment also provides a computer-readable storage medium storing a computer program, which, when run by a processor, executes the steps of the above-described wall wear simulation method independent of mesh size.

[0123] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems and devices described above can be referred to the corresponding processes in the method embodiments, and will not be repeated here. In the several embodiments provided in this application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of modules is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed mutual coupling or direct coupling or communication connection can be through some communication interfaces; the indirect coupling or communication connection of devices or modules can be electrical, mechanical, or other forms.

[0124] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0125] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0126] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, ROM, RAM, magnetic disks, or optical disks.

[0127] The above are merely specific embodiments of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for simulating wall wear independent of mesh size, characterized in that, The method includes: Obtain the solid wall boundary model file of the structure, and a background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch elements in the solid wall boundary model file; The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, as well as the wall wear volume caused by the collision contact point. Based on the wall wear volume generated by the collision at the collision point, the wall wear volume of the first grid node in the background grid closest to the collision point is accumulated. Continue with the collision simulation at the next simulation time step; and based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, perform wear diffusion on the wall of the structure to obtain a distribution cloud map of the wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

2. The wall wear simulation method independent of mesh size according to claim 1, characterized in that, The following steps are used to obtain a background mesh surrounding the solid wall boundary of the structure, which is independent of the patch element size of the solid wall boundary model file: The node spacing of the background mesh is calculated based on the preset minimum size of the granular particles. Based on the geometric dimensions of the solid wall boundary model file, determine the minimum geometric bounding box of the solid wall boundary model; The background mesh is obtained by dividing the interior of the geometric bounding box into background meshes based on the node spacing.

3. The wall wear simulation method independent of mesh size according to claim 1, characterized in that, The solid wall boundary model in the solid wall boundary model file is subjected to collision simulation using the discrete element method to obtain the wall wear volume caused by the collision contact point, including: The discrete element method is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file to obtain the sliding distance of the granular particles on the wall surface of the structure, as well as the interaction force between the granular particles and the structure. The volume of wall wear caused by the collision at the point of impact is calculated based on the sliding distance of the granular particles on the wall surface of the structure and the interaction force between the granular particles and the structure.

4. The wall wear simulation method independent of mesh size according to claim 1, characterized in that, The first grid node closest to the collision point is determined using the following formula: ； in, The coordinates of the first grid node are... The X-axis coordinate of the first grid node is... The Y-axis coordinate of the first grid node. The Z-axis coordinate of the first grid node is... This indicates that the result is rounded to the nearest integer. The X-axis coordinates of the collision contact points. The Y-axis coordinate of the collision contact point. The Z-axis coordinate of the collision contact point. The starting point of the background grid X-axis coordinates The starting point of the background grid Y-axis coordinate, The starting point of the background grid Z-axis coordinate, This represents the node spacing corresponding to the background grid.

5. The wall wear simulation method independent of mesh size according to claim 1, characterized in that, The step of performing wear diffusion on the wall surface of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid to obtain a distribution cloud map of the wear depth on the wall grid includes: Determine whether the current simulation time step is the preset output time step for the wall wear simulation; If so, then based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, wear diffusion is carried out on the wall of the structure to obtain a distribution cloud map of wear depth on the wall grid.

6. The wall wear simulation method independent of mesh size according to claim 1 or 5, characterized in that, The step of performing wear diffusion on the wall surface of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid to obtain a distribution cloud map of the wear depth on the wall grid includes: For each second grid node, the wear impact value of the second grid node on the corresponding sampling point is calculated based on the distance between the second grid node and the corresponding sampling point; Using the wear impact value of each second grid node on the corresponding sampling point as a weight, the latest accumulated wall wear volume of all second grid nodes is weighted and summed to obtain the wear depth at the sampling point in the distribution cloud map.

7. The wall wear simulation method independent of mesh size according to claim 6, characterized in that, The step of calculating the wear impact value of the second grid node on the corresponding sampling point based on the distance between the second grid node and the corresponding sampling point includes: Substitute the distance between the second grid node and the corresponding sampling point into the diffusion function to calculate the wear influence of the second grid node on the corresponding sampling point. The diffusion function satisfies the following condition: ； ； ； in, This represents the wear impact value when the distance between the second grid node and the corresponding sampling point is 0. This is the preset maximum wear impact value; The distance between the second grid node and the corresponding sampling point is The wear impact value at that time; R is the preset wear diffusion radius; This represents the distance between the second grid node and the corresponding sampling point. The azimuth angle of the second grid node relative to the corresponding sampling point.

8. A wall wear simulation device independent of grid size, characterized in that, The device includes: The acquisition module is used to acquire the solid wall boundary model file of the structure, and the background mesh surrounding the solid wall boundary of the structure that is independent of the size of the patch unit of the solid wall boundary model file. The calculation module is used to perform collision simulation on the solid wall boundary model in the solid wall boundary model file using the discrete element method, to obtain the collision contact point between the granular particles and the wall of the structure at the current simulation time step, and the wall wear volume caused by the collision contact point. The calculation module is also used to accumulate the wall wear volume of the first grid node in the background grid that is closest to the collision contact point based on the wall wear volume generated by the collision contact point in this collision contact. The calculation module is also used to continue the collision simulation for the next simulation time step; and to perform wear diffusion on the wall of the structure based on the latest accumulated wall wear volume of the second grid node corresponding to the sampling point in the background grid, so as to obtain the distribution cloud map of wear depth on the wall grid; the distance between the sampling point and the corresponding second grid node is less than or equal to the preset wear diffusion radius.

9. An electronic device, characterized in that, include: The device includes a processor, a storage medium, and a bus, wherein the storage medium stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor communicates with the storage medium via the bus, and the processor executes the machine-readable instructions to perform the steps of the mesh-size-independent wall wear simulation method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the steps of the wall wear simulation method independent of mesh size as described in any one of claims 1 to 7.

Images