Simulation method and system based on multi-field coupling and dynamic grid processing

Abstract

The invention discloses a simulation method and system based on multi-field coupling and dynamic grid processing, and belongs to the field of computational simulation, and the method comprises the steps: building and solving a coupling physical field model: carrying out finite element spatial discretization on a solving domain, building and solving a coupling partial differential equation set which at least comprises a temperature field, a displacement field and a stress field, and carrying out the finite element spatial discretization; obtaining numerical solutions of the physical fields in the current state; carrying out progressive unit removal based on a physical criterion; and finally, grid data structure updating and solver synchronization are executed. According to the method, the metal material can be gradually removed when the physical quantity reaches the critical condition, and the grid change is dynamically and robustly processed in the process.

Description

Technical Field

[0001] The present invention relates to the field of computational simulation, and more specifically, to a simulation method and system based on multi-field coupling and dynamic mesh processing. Background Art

[0002] In many engineering applications, especially in the aerospace field (such as thermal protection systems, engine components, etc.), metallic materials are often in a coupled physical environment of high temperature, high stress, and high-speed gas flow erosion, and undergo complex performance degradation and even mass loss processes (such as ablation, erosion). Accurately simulating these processes is crucial for component design, life prediction, and safety assessment.

[0003] Existing numerical simulation methods face challenges: 1) How to effectively handle the dynamic evolution of the computational domain boundary shape and topological structure caused by material loss, and traditional fixed meshes are difficult to apply; 2) How to establish a material removal criterion and a mesh processing strategy that can both reflect the physical mechanism and ensure computational stability, and avoid, for example, numerical instability or the generation of suspended "island" elements that do not conform to physical reality. Summary of the Invention

[0004] The purpose of the present invention is to overcome the deficiencies of the prior art, and provide a simulation method and system based on multi-field coupling and dynamic mesh processing, to achieve the progressive removal of metallic materials when physical quantities (such as stress) reach critical conditions, and dynamically and robustly handle mesh changes during this process.

[0005] The purpose of the present invention is achieved by the following solutions: A simulation method based on multi-field coupling and dynamic mesh processing, comprising the following steps: Establish and solve a coupled physical field model: perform finite element spatial discretization on the solution domain, establish and solve a coupled partial differential equation system describing at least a temperature field, displacement, and stress field, and obtain numerical solutions of these physical fields in the current state; Then perform progressive element removal based on physical criteria; Finally, perform mesh data structure update and solver synchronization.

[0006] Furthermore, the material constitutive relationship of the coupled physical field model includes the coupling effect of the temperature field on mechanical behavior.

[0007] Furthermore, the progressive element removal based on physical criteria specifically includes the following sub-steps: S11, state check and element marking: Based on the numerical solutions of the obtained physical fields, traverse the active elements within a predetermined inspection range, and mark the elements that meet the removal conditions as "elements to be removed" according to at least one preset physical quantity threshold; S12. Identification and Update of Newly Exposed Boundaries: For each cell E_del marked as "cell to be removed", before it is physically removed, check all its active neighbor cells N that are not marked as cells to be removed. If the shared internal face F_internal between E_del and N does not currently belong to the "exposed boundary", update F_internal to the new "exposed boundary" and confirm it as the newly eroded face. S13. Physical Removal of Cells: Delete all cells currently marked as "cells to be removed" from the finite element mesh data structure. S14. Iterative Erosion: Repeat steps S11 to S13, with the inspection scope focused on all active cells that are currently "exposed boundaries", until the iterative termination condition is met.

[0008] Furthermore, between the progressive cell removal based on physical criteria and the execution of mesh data structure update and solver synchronization, the following steps are also included: S21. After the progressive cell removal operation is completed, starting from the predefined "stable boundary", identify all active cells connected to the "stable boundary" through a graph search algorithm. S22. Determine all active cells that are not identified as connected as "isolated cells". S23. Delete all identified "isolated cells" from the finite element mesh data structure.

[0009] Furthermore, the execution of mesh data structure update and solver synchronization specifically includes the following sub-steps: S31. Update Connection Relationships and States: After the cell removal is completed, modify the mesh data structure to reflect the changes in the connection relationships between nodes and cells caused by the cell removal, as well as the update of the active states of cells and nodes. S32. Reconstruct Adjacency and Boundary Information: Based on the updated connection relationships, recalculate or update the adjacency information between active cells, and maintain or reconstruct the geometric entities and their topological information that constitute all currently defined boundaries. S33. Synchronize Parallel Environment Data for Parallel Computing: Through inter-processor communication, synchronize the geometric and topological information of the shared boundaries between the mesh parts held by each processor caused by the cell removal to ensure the consistency of the global mesh data. S34. Notify the Main Solving Framework: Through an internal mechanism, notify the main finite element solving framework that the mesh data structure has changed to trigger the re-initialization of the internal state of the solver.

[0010] Furthermore, the coupled partial differential equations are specifically solved using the Newton method for nonlinear iterative solution, where the linear equations generated in each Newton iteration are solved using the preconditioned iterative method.

[0011] Further, in step S11, the state check and cell marking include sub-steps: If it is a state check and cell marking based on stress, the inspection scope is all active cells currently adjacent to at least one defined "exposed boundary"; and the stress index within each inspected cell also needs to be calculated.

[0012] Further, in step S14, the iteration termination conditions include at least one of the following: In the state check and cell marking step of one round of iteration, no new cells are marked as "cells to be removed"; or At the beginning of the iteration, the "set of cells to be removed" is empty; or The preset maximum number of iterations is reached.

[0013] Further, in step S21, the graph search algorithm includes breadth-first search BFS, depth-first search DFS, or connectivity analysis method.

[0014] A simulation system based on multi-field coupling and dynamic mesh processing includes a computer device, where the computer device includes a processor and a memory, and a computer program is stored in the memory. When the computer program is loaded by the processor, the method described in any one of the above is executed.

[0015] The beneficial effects of the present invention include: (1) Simulating the progressive process: Especially the stress-based iterative removal mechanism can simulate the physical process of material loss layer by layer, rather than removing it all at once.

[0016] (2) Improving robustness: The island cell detection and removal mechanism significantly improves the numerical stability and physical rationality of the large-scale cell removal simulation.

[0017] (3) Flexibility and applicability: It can be applied to material failure criteria based on thresholds (such as temperature, stress, etc.), and can provide a basic framework for dynamic mesh processing for more complex aerospace ablation / erosion models. For ablation / erosion models, it is necessary to further expand the influence of the boundary, consider the ablation effect brought by chemical reactions, and the erosion effect brought by particle impact.

[0018] (4) Considering the coupling effect: By executing within the coupled solution framework, it ensures that the mutual influence between the physical field evolution and the geometric shape change is considered. Description of the Drawings

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings.

[0020] Figure 1a is the grid deletion strategy based on temperature; Figure 1b is the grid deletion strategy based on stress; Figure 2a is the schematic diagram of grid deletion based on temperature; Figure 2b is the flowchart of grid deletion based on stress; Figure 3 is the time-varying calculation result of grid deletion based on temperature; Figure 4a is the temperature distribution diagram in the time-varying calculation result of grid deletion based on stress; Figure 4b is the stress distribution diagram in the time-varying calculation result of grid deletion based on stress; where vonmises_stress represents stress. Detailed implementation manners

[0021] All the features disclosed in all the embodiments in this specification, or all the steps in the methods or processes implicitly disclosed, can be combined and / or extended, replaced in any way, except for mutually exclusive features and / or steps.

[0022] The present invention aims to provide a numerical method that can effectively handle the progressive removal of metal materials due to exceeding physical thresholds (such as stress) in the framework of coupled multi-physics finite element simulation, and dynamically update the calculation domain boundary and maintain numerical stability during this process. The main problems to be solved are how to automatically identify and mark failed elements, how to update the exposed boundary according to the removal process, how to avoid numerical islands caused by element removal, and how to effectively integrate these processes with the main physical field solver when simulating the layer-by-layer loss of materials.

[0023] Specifically, within a coupled multi-physics finite element calculation framework, a simulation method based on multi-field coupling and dynamic grid processing is provided, including the following steps: Step (a), establish and solve a coupled physical field model.

[0024] The solution domain is discretized in finite element space, and a set of coupled partial differential equations describing at least the temperature field and displacement / stress field is established and solved to obtain the numerical solution of these physical fields in the current state. Among them, the material constitutive relationship of the model includes the coupling effect of the temperature field on the mechanical behavior, such as the influence of thermal strain on the total strain and stress calculation.

[0025] Specifically, the finite element method (FEM) is used to spatially discretize the solution domain. This provides the foundation for subsequent element-based physical quantity calculations and element removal operations. For the physical phenomena that need to be simulated in the physical computational domain, a mathematical description is established, namely a system of partial differential equations (PDEs), which mainly include: Energy conservation (non-steady-state heat conduction): describes the spatiotemporal evolution of the temperature field T. Its strong form (not considering volume heat sources) is:

[0026] The corresponding weak form is: for all suitable test functions , solve T so that:

[0027] in, is the material density, c p is the material specific heat capacity, k is the material thermal conductivity, q n is through the boundary Γ q The normal heat flux density, n is the boundary Γ q The unit external normal vector, V is the infinitesimal volume, S is the infinitesimal area, is the entire computational domain. Formula (2) represents the transient term, diffusion term, and boundary heat flux from left to right.

[0028] Conservation of momentum (quasi-static force balance equation): describes the displacement field u and stress field of the material under load and temperature changes Since the inertia term is not considered, the equation is simplified to the force balance equation, and its strong form is:

[0029] The corresponding weak form is: for all suitable test functions , solve u so that:

[0030] where σ is the Cauchy stress tensor, is the volume force, is applied on the boundary Γ t The surface force on it.

[0031] Material constitutive theory and coupling: a bridge connecting different physical fields.

[0032] Thermal strain: When a structure is under the action of external heat flux and pressure, the material deforms. For example, linear thermal expansion , where α is the coefficient of thermal expansion, Tref is the reference temperature, and I is the second-order identity tensor.

[0033] Strain decomposition: The total strain is calculated from the displacement gradient and is decomposed into mechanical strain and thermal strain , that is .

[0034] Elastic stress: Mechanical strain generates stress through the elastic tensor C: .

[0035] Coupling relationship: Through the above relationships, changes in the temperature field T(x, t) will cause thermal strain, which in turn affects the mechanical strain and stress field σ, constituting the thermal driving force on the force.

[0036] Numerical solution strategy: Time discretization: For transient problems, time integration methods (e.g., implicit backward Euler method, Crank-Nicolson method, etc.) are used to discretize the time derivative.

[0037] Nonlinear solution: The coupled system of equations usually forms a large-scale nonlinear algebraic equation system R(U,t)=0 after discretization, where U is the global vector containing all the unknown degrees of freedom of the nodes (such as temperature, displacement components). The Newton method is used for iterative solution: At the k-th iteration, calculate the residual vector R(Uk).

[0038] Calculate the Jacobian matrix (tangent stiffness matrix) .

[0039] Solve the linear system to obtain the solution increment .

[0040] Update the solution , where λ is the line search step size.

[0041] Repeat until the residual satisfies the convergence criterion.

[0042] Linear solution: The above linear system is usually large-scale and sparse, and is solved using a preconditioned iterative method, such as the Generalized Minimum Residual method (GMRES) or the Bi-Conjugate Gradient Stabilized method (BiCGStab), and is combined with an efficient preconditioner, such as Algebraic Multigrid (AMG) or Incomplete LU factorization.

[0043] Step (b), implement progressive element removal based on physical criteria.

[0044] 1) State check and element marking: Based on the physical field solution obtained in step (a), traverse the active elements within a predetermined check range, and mark the elements that meet the removal conditions as "elements to be removed" according to at least one preset physical quantity threshold (for example, the element stress index exceeds the stress threshold); 2) Automatic identification and update of new exposed boundaries: For each element E_del marked as "element to be removed", before it is physically removed, check all its active neighbor elements N that are not marked as elements to be removed. If the shared internal face F_internal between E_del and N does not currently belong to the "exposed boundary", then update F_internal to the new "exposed boundary" and confirm it as the new eroded face; 3) Physical removal of elements: Delete all elements currently marked as "elements to be removed" from the finite element mesh data structure; 4) Iterative erosion strategy: Repeat steps 1) based on the stress process, 2) and 3), with the check range focused on all active elements that are currently "exposed boundaries" until the iterative termination condition is met.

[0045] Specifically, at a specific stage of the solution process (for example, after each time step), execute the following element removal process based on physical criteria. Here are two methods: one is the mesh deletion method based on the temperature threshold, and the other is the mesh deletion method based on the stress threshold. The mesh deletion method based on temperature is used to test the functional practicality of the algorithm and is not applicable to the actual environment; while the mesh deletion function based on the stress threshold can be used in the simulation process of the progressive loss of actual metal materials.

[0046] 1) State check and marking ( Figure 1a step B in ): Figure 1a Based on temperature: Traverse the temperatures of all elements on the eroded boundaries. Calculate the temperature indices of these elements (such as the average element temperature, the maximum temperature at the integration point, etc.). If the temperature index exceeds the preset temperature threshold, then mark the element as "to be removed". As

[0047] shown. Figure 1b Based on stress: Traverse the stresses of all elements on the eroded boundaries. Calculate the stress indices of these elements (such as the average element stress, the maximum von Mises stress at the integration point). If the stress index exceeds the preset stress threshold, then mark the element as "to be removed". As

[0048] shown. Figure 1a and Figure 1bStep C): For each unit E_del marked as "to be removed": Traverse all its neighboring units N. If the neighbor N is active and not marked as "to be removed", add the internal shared face F_internal between the unit E_del and the neighbor N to the definition of the "exposed boundary", thereby identifying the new eroded face. This step is completed before physically removing the unit, ensuring the exposure of the new surface after simulating the material removal.

[0049] 3) Physical removal of units ( Figure 1a and Figure 1b Step D): Delete all units marked as "to be removed" from the mesh data structure.

[0050] 4) Iterative erosion (for the stress criterion, ( Figure 1b Steps E, C, and D in)): Repeat the stress-based process in Step 1) to Step 3). That is, on the newly exposed boundary and in the corresponding non-removed units on the original exposed boundary, check again whether the stress exceeds the limit, and perform marking, boundary update, and removal. This iterative process simulates the layer-by-layer erosion of the material until no new units are marked in an iterative cycle or the preset maximum number of iterations is reached. For the temperature criterion, there is no such limitation. When the temperature of the eroded surface is higher than the threshold, it is directly deleted in its entirety (as Figure 2a shown). Due to the material being subjected to thermal / mechanical loads in the stress criterion, large stresses may also occur inside the material. The deletion condition of the material often also applies to the internal units, and the loss of the material often develops from the eroded surface. Therefore, the retrieval method is limited to layer-by-layer inspection and recursive deletion of the mesh (as Figure 2b shown).

[0051] Step (c), Detection and handling of isolated island units (optional).

[0052] After performing the removal operation in Step (b), starting from the predefined "stable boundary", identify all active units connected to the "stable boundary" through a graph search algorithm; Determine all active but unrecognized connected units as "isolated island units"; Delete all identified "isolated island units" from the finite element mesh data structure.

[0053] Specifically, it is recommended to add the connectivity analysis of this item under the stress criterion, ( Figure 1b Step F in)): After confirming that no unit exceeds the limit and is marked, or reaching the preset maximum number of iterations, start the detection and analysis of "isolated island" units, mainly for connectivity analysis, including: Define one or more "stable boundary" identifiers, which represent the boundaries of regions that are fixed or not completely removed in the actual process.

[0054] Starting from all active cells adjacent to the "stable boundary".

[0055] Use a graph search algorithm (such as breadth - first search BFS or depth - first search DFS) to traverse all active cells connected to each other through shared faces or nodes, and mark all reachable cells.

[0056] Identify all active cells in the grid that are not marked as reachable. These are the "island" cells.

[0057] Remove all identified "island" cells from the finite - element grid data structure.

[0058] Step (d), grid data structure update and solver synchronization. Specifically, perform the following steps: 1) Update connection relationships and states: After removing cells in steps (b) and / or (c), modify the grid data structure to reflect the changes in the connection relationships between nodes and cells caused by cell removal, as well as the update of the active states of cells and nodes; 2) Re - build adjacency and boundary information: Based on the updated connection relationships, recalculate or update the adjacency information between active cells, and maintain or re - build the geometric entities and their topological information that form all currently defined boundaries (including newly generated "exposed boundaries"); 3) (For parallel computing) Synchronize parallel environment data: Through inter - processor communication, synchronize the geometric and topological information of the shared boundaries between the parts of the grid held by each processor caused by cell removal to ensure the consistency of the global grid data; 4) Notify the main solution framework: Through an internal mechanism, notify the main finite - element solution framework that the grid data structure has changed to trigger the re - initialization of necessary solver internal states (such as degree - of - freedom mapping, sparse matrix pattern).

[0059] To more clearly illustrate the technical solution of the present invention, the present invention will be described in detail below with reference to the accompanying drawings. Figures 1a - 1b The flowcharts of two grid cell removal strategies based on different physical criteria in the method of the present invention are shown. Among them, Figure 1aIllustrates the temperature-based mesh deletion strategy. Within a time step, this strategy mainly performs a one-time check: after the coupled-field calculation (at the start of the time step), it checks the element temperatures on the initial exposed boundary (step B). If the temperature exceeds the threshold, it updates the adjacent boundary (step C) and deletes the marked elements (step D), and then enters the end-of-time-step phase for mesh update. This represents a pattern of element removal based on the current temperature field snapshot. This mesh deletion strategy is rather simplistic. In an actual working environment, material deletion due to temperature does not conform to physical laws. On this basis, a simulation of the stress-induced material deletion strategy is carried out to simulate the progressive damage process of metal materials in a real thermo-mechanical coupling environment. Figure 1b Illustrates the stress-based mesh deletion strategy. This strategy embodies a more complex iterative, layer-by-layer removal logic: after the coupled-field calculation, it checks and marks the stresses on the initial exposed boundary (step B); if any elements are marked, it updates the adjacent boundary (step C), and then deletes the marked elements (step D); the key lies in that after the deletion operation, it checks the element stresses on all current exposed boundaries (including the newly exposed boundaries) again (step E). If new over-limit elements are found, it returns to step C, forming an iterative loop of step C, step D, step E, step C to simulate the layer-by-layer erosion of the material until no new over-limit elements are found; after the entire iterative removal process ends (regardless of whether any elements are removed), it performs island detection and deletion (step F) to ensure the physical connectivity of the mesh; finally, it enters the end-of-time-step phase for the final mesh update. The core of this strategy is to simulate the stress-driven progressive physical erosion process of the material by iteratively checking the newly exposed surfaces.

[0060] Figures 2a - 2b Schematically shows the possible results of the temperature-based mesh deletion strategy (as shown in Figure 1a the process). According to the heat flux distribution (here, the non-uniform boundary influence is shown), the regions with temperature overrun may cause the removal of single or multiple layers of elements at once, reflecting the instantaneous state of reaching the heat threshold rather than a layer-by-layer progression. The corresponding calculation results are as shown in Figure 3 (the top corresponds to the heat flux boundary, the left and right sides are adiabatic boundaries, and the bottom is at a constant temperature). The initial heat flux boundary given here is a function of coordinates. When new exposed boundaries appear, new heat flux boundaries are generated based on the coordinates of the new boundaries for subsequent calculations. The same method can be extended to pressure boundaries. Subsequently, the initial heat flux distribution transferred from the internal and external flow fields in the aerospace field to the structure is parameterized, and the newly generated boundaries are updated with parameterized heat flux boundaries based on coordinates. Implementing boundary update based on coordinates can better conform to the actual environment faced by the new boundaries, thus simulating the retreat of material loss.

[0061] Figure 3Schematically shows the core logic and effect of the stress-based mesh deletion strategy (such as Figure 1b process). It emphasizes the iterative logic added within a time step: by checking layer by layer (corresponding to Figure 1b the C ->D ->E loop), deleting elements, performing mesh topology updates, determining newly exposed boundaries through neighborhood search, and then checking again whether the elements on these new boundaries exceed the limit, thus achieving the effect of progressive removal of materials from the initial exposed surface towards the interior, which is closer to the physical erosion process. The corresponding calculation results are as Figures 4a - 4b shown (corresponding to a constant heat flux and pressure boundary at the top, adiabatic boundaries on the left and right sides, a constant temperature of 300K and a fixed support boundary at the bottom). For this kind of mesh deletion strategy, it is mainly considered that the material loss starts from the eroded surface, and stress exceeding the limit does not only exist on the eroded surface, and it is necessary to check layer by layer and delete layer by layer starting from the top.

[0062] It should be noted that Figures 2a - 2b and Figure 3 are numerical simulation processes carried out under two different calculation boundaries.

[0063] It should be noted that in other embodiment scenarios of the present invention, it can be implemented on the basis of a general finite element analysis software platform by developing user-defined modules (such as User Objects or subroutines), for example, including: Data structure: Use sets to store the element IDs of the elements to be removed, already removed, and newly exposed boundaries.

[0064] Mesh interface: It is necessary to access and modify the interfaces provided by the finite element mesh library to obtain element adjacency relationships, boundary information, and perform element deletion, boundary addition / deletion operations.

[0065] Module integration: Embed the element removal and island detection modules into the process control of the main solver, for example, call after the solution of each time step is completed and before the mesh is updated.

[0066] Information transfer: Ensure that the physical field solver can access the latest mesh information, and the element removal module can obtain the physical field variable values (temperature, stress, etc.) calculated at the current time step.

[0067] Connectivity algorithm: Graph search algorithms (BFS or DFS) are used for island detection.

[0068] The above embodiments are all within the protection scope of the claims of the present invention.

[0069] The units involved in the embodiments of the present invention can be implemented in software or in hardware, and the described units can also be provided in a processor. Wherein, the names of these units do not, in some cases, constitute a limitation on the units themselves.

[0070] According to one aspect of the embodiments of the present invention, there is provided a computer program product or a computer program, the computer program product or the computer program including computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the computer device executes the methods provided in the above various alternative implementation manners.

[0071] As another aspect, the embodiments of the present invention further provide a computer-readable medium, which may be included in the electronic device described in the above embodiments; or may exist separately without being assembled into the electronic device. The above computer-readable medium carries one or more programs, and when the above one or more programs are executed by an electronic device, the electronic device implements the methods described in the above embodiments.

Claims

1. A simulation method based on multi-field coupling and dynamic mesh processing, characterized in that, The method includes the following steps: Establish and solve a coupled physical field model: discretize the solution domain in the finite element space, establish and solve a coupled partial differential equation system that describes at least the temperature field, displacement, and stress field, and obtain the numerical solutions of these physical fields in the current state; Then perform progressive element removal based on physical criteria; Finally, perform mesh data structure update and solver synchronization.

2. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, wherein The material constitutive relationship of the coupled physical field model includes the coupling effect of the temperature field on mechanical behavior.

3. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, wherein The progressive element removal based on physical criteria specifically includes the following sub-steps: S11, State check and element marking: Based on the numerical solutions of the obtained physical fields, traverse the active elements within a predetermined inspection range, and mark the elements that meet the removal conditions as "elements to be removed" according to at least one preset physical quantity threshold; S12, Identification and update of newly exposed boundaries: For each element E_del marked as "element to be removed", before it is physically removed, check all its active neighbor elements N that are not marked as elements to be removed. If the shared internal face F_internal between E_del and N does not currently belong to the "exposed boundary", then update F_internal to the new "exposed boundary" and confirm it as the new eroded face; S13, Physical removal of elements: Delete all elements currently marked as "elements to be removed" from the finite element mesh data structure; S14, Iterative erosion: Repeat steps S11 to S13, and focus the inspection range on all active elements that are currently "exposed boundaries" until the iterative termination condition is met.

4. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, wherein Between the progressive element removal based on physical criteria and the execution of mesh data structure update and solver synchronization, there is also a step: S21, After performing the progressive element removal operation, starting from the predefined "stable boundary", identify all active elements connected to the "stable boundary" through a graph search algorithm; S22, Determine all active elements that are not identified as connected as "island elements"; S23, Delete all identified "island elements" from the finite element mesh data structure.

5. The simulation method based on multi-field coupling and dynamic grid processing according to any one of claims 1 to 4, characterized in that The execution of mesh data structure update and solver synchronization specifically includes the following sub-steps: S31, Update connection relationship and state: After performing element removal, modify the mesh data structure to reflect the changes in the connection relationship between nodes and elements caused by element removal, as well as the update of the active states of elements and nodes; S32, Reconstruct adjacency and boundary information: Based on the updated connection relationship, recalculate or update the adjacency information between active elements, and maintain or reconstruct the geometric entities and their topological information that constitute all currently defined boundaries; S33, Synchronize parallel environment data for parallel computing: Through inter-processor communication, synchronize the geometric and topological information of the shared boundaries between the mesh parts held by each processor caused by element removal to ensure the consistency of the global mesh data; S34, Notify the main solution framework: Notify the main finite element solution framework through an internal mechanism that the mesh data structure has changed to trigger the re-initialization of the internal state of the solver.

6. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, characterized in that To solve the coupled partial differential equations, the Newton method is specifically used for non-linear iterative solution, and the linear equations generated in each Newton iteration are solved by the preconditioned iterative method.

7. The simulation method based on multi-field coupling and dynamic grid processing according to claim 3, characterized in that In step S11, the state check and element marking include sub-steps: If it is a state check and element marking based on stress, the inspection range is all active elements currently adjacent to at least one defined "exposed boundary"; the stress index within each inspected element also needs to be calculated.

8. The simulation method based on multi-field coupling and dynamic grid processing according to claim 3, characterized in that In step S14, the iteration termination conditions include at least one of the following: In the state check and element marking step of one round of iteration, no new element is marked as a "to-be-removed element"; or At the beginning of the iteration, the "set of to-be-removed elements" is empty; or The preset maximum number of iterations is reached.

9. The simulation method based on multi-field coupling and dynamic grid processing according to claim 4, wherein In step S21, the graph search algorithm includes breadth-first search BFS, depth-first search DFS, or a connectivity analysis method.

10. A simulation system based on multi-field coupling and dynamic grid processing, characterized in that It includes a computer device, the computer device includes a processor and a memory, and a computer program is stored in the memory. When the computer program is loaded by the processor, it executes the method according to any one of claims 1 to 4.

Images