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

Abstract

The present invention discloses a simulation method and system based on multi-field coupling and dynamic mesh processing, belonging to the field of computational simulation. The method comprises the following steps: establishing and solving a coupled physical field model; discretizing the solution domain using finite element space, establishing and solving a set of coupled partial differential equations describing at least the temperature field, displacement field, and stress field, and obtaining numerical solutions for these physical fields in their current state; performing progressive element removal based on physical criteria; and finally performing mesh data structure updates and solver synchronization. The present invention can achieve progressive removal of metal materials when physical quantities reach critical conditions, and dynamically and robustly handle mesh changes during this 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 grid processing. Background Art

[0002] In many engineering applications, particularly in aerospace (e.g., thermal protection systems and engine components), metal materials are often exposed to coupled physical environments such as high temperature, high stress, and high-velocity airflow, causing them to undergo complex performance degradation and even mass loss processes (e.g., ablation and erosion). Accurately simulating these processes is crucial for component design, life prediction, and safety assessment.

[0003] Existing numerical simulation methods face the following challenges: 1) How to effectively deal with the dynamic evolution of the boundary shape and topology of the computational domain caused by material loss, which is difficult to apply to traditional fixed grids; 2) How to establish material removal criteria and grid processing strategies that both reflect physical mechanisms and ensure computational stability, avoiding, for example, numerical instability or the generation of suspended "island" units that do not conform to physical reality. Summary of the Invention

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

[0005] The object of the present invention is achieved through the following solutions:

[0006] A simulation method based on multi-field coupling and dynamic grid processing includes the following steps:

[0007] Establish and solve the coupled physical field model: perform finite element spatial discretization on the solution domain, establish and solve a set of coupled partial differential equations describing at least the temperature field, displacement, and stress field, and obtain the numerical solutions of these physical fields in the current state;

[0008] Then perform progressive unit removal based on physical criteria;

[0009] Finally, the grid data structure update and solver synchronization are performed.

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

[0011] Furthermore, the physical criterion-based progressive unit removal specifically includes the following sub-steps:

[0012] S11, status check and unit marking: based on the obtained numerical solution of the physical field, traverse the active units within the predetermined inspection range, and mark the units that meet the removal conditions as "units to be removed" according to at least one preset physical quantity threshold;

[0013] S12, identification and update of new exposed boundaries: For each unit E_del marked as a "to-be-removed unit", before it is physically removed, check all its active neighboring units N that are not marked as to-be-removed. If the shared internal surface F_internal between E_del and N is not currently an "exposed boundary", then update F_internal to the new "exposed boundary" and confirm it as a new eroded surface;

[0014] S13, physical removal of elements: deleting all elements currently marked as “elements to be removed” from the finite element mesh data structure;

[0015] S14, iterative erosion: repeat steps S11 to S13, focusing the inspection range on all active cells that are currently "exposed boundaries" until the iteration termination condition is met.

[0016] Furthermore, between the progressive element removal based on physical criteria and the execution of mesh data structure update and solver synchronization, the following steps are also included:

[0017] S21, after performing the progressive unit removal operation, starting from the predefined “stable boundary”, all active units connected to the “stable boundary” are identified through a graph search algorithm;

[0018] S22, all units that are active but not identified as connected are judged as “island units”;

[0019] S23, deleting all identified “island elements” from the finite element mesh data structure.

[0020] Furthermore, the execution of the grid data structure update is synchronized with the solver, specifically including the following sub-steps:

[0021] S31, updating connection relationships and states: after the unit removal is completed, modifying the grid data structure to reflect the changes in the connection relationships between nodes and units caused by the unit removal, as well as the updates of the activity states of the units and nodes;

[0022] S32, rebuilding adjacency and boundary information: recalculating or updating adjacency information between active units based on the updated connection relationship, and maintaining or rebuilding geometric entities and their topological information constituting all currently defined boundaries;

[0023] S33, synchronizing parallel environment data for parallel computing: synchronizing the geometry and topology information of the shared boundaries between the mesh parts held by each processor due to unit removal through inter-processor communication, to ensure the consistency of global mesh data;

[0024] S34, notifying the main solution framework: notifying the main finite element solution framework through an internal mechanism that the mesh data structure has changed to trigger the reinitialization of the internal state of the solver.

[0025] Furthermore, the coupled partial differential equations are solved by nonlinear iterative solution using the Newton method, wherein the linear equations generated in each Newton iteration are solved using a preconditioned iterative method.

[0026] Furthermore, in step S11, the status check and unit marking includes the following sub-steps:

[0027] If the state check and element marking is based on stress, the inspection scope is all active elements that are currently adjacent to at least one defined "exposure boundary"; and the stress index in each inspected element is also calculated.

[0028] Furthermore, in step S14, the iteration termination condition includes at least one of the following:

[0029] During the status check and unit marking steps of an iteration, no new units are marked as "units to be removed"; or

[0030] At the beginning of the iteration, the "set of units to be removed" is empty; or

[0031] The preset maximum number of iterations has been reached.

[0032] Furthermore, in step S21, the graph search algorithm includes breadth-first search (BFS), depth-first search (DFS), or connectivity analysis method.

[0033] A simulation system based on multi-field coupling and dynamic grid processing includes a computer device, wherein the computer device includes a processor and a memory, wherein a computer program is stored in the memory, and when the computer program is loaded by the processor, any of the methods described above is executed.

[0034] The beneficial effects of the present invention include:

[0035] (1) Simulating the gradual process: In particular, the stress-based iterative removal mechanism can simulate the physical process of material loss layer by layer rather than a one-time removal.

[0036] (2) Improved robustness: The island unit detection and removal mechanism significantly improves the numerical stability of large-scale unit removal simulations and the physical rationality of the results.

[0037] (3) Flexibility and applicability: It can be applied to threshold-based material failure criteria (temperature, stress, etc.) and can provide a basic framework for dynamic mesh processing for more complex aerospace ablation / erosion models. For ablation / erosion models, the influence of boundaries needs to be further expanded to consider the ablation effect caused by chemical reactions and the erosion effect caused by particle impact.

[0038] (4) Consideration of coupling effects: By executing within the coupled solution framework, the mutual influence between the evolution of the physical field and the change of the geometric shape is taken into account. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] 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 only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0040] Figure 1a is a temperature-based mesh deletion strategy;

[0041] Figure 1b is a stress-based mesh deletion strategy;

[0042] Figure 2a Delete schematic for temperature-based mesh;

[0043] Figure 2b Removed flow chart for stress-based mesh;

[0044] Figure 3 Time-varying calculation results for temperature-based grid deletion;

[0045] Figure 4a Temperature distribution diagram in time-varying calculation results after stress-based mesh deletion;

[0046] Figure 4b This is a stress distribution diagram from the time-varying calculation results of stress-based mesh deletion. vonmises_stress represents stress. DETAILED DESCRIPTION

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

[0048] 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 a coupled multi-physics field finite element simulation framework, and dynamically update the boundaries of the computational domain and maintain numerical stability in this process. It mainly solves the problems of how to automatically identify and mark failed units when simulating layer-by-layer material loss, how to update the exposed boundaries according to the removal process, how to avoid numerical islands caused by unit removal, and how to effectively integrate these processes with the main physical field solver.

[0049] Specifically, within a coupled multi-physics field finite element calculation framework, a simulation method based on multi-field coupling and dynamic mesh processing is provided, including the following steps:

[0050] Step (a): establish and solve the coupled physical field model.

[0051] 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.

[0052] 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:

[0053] 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:

[0054]

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

[0056]

[0057] 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.

[0058] 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:

[0059]

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

[0061]

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

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

[0064] Thermal strain: The deformation of the material under the influence of external heat flow and pressure, such as linear thermal expansion , where α is the coefficient of thermal expansion, Tref is the reference temperature, and I is the second-order unit tensor.

[0065] Strain decomposition: Total strain Calculated from the displacement gradient , and is decomposed into mechanical strain and thermal strain ,Right now .

[0066] Elastic stress: Mechanical strain produces stress through the elastic tensor C: .

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

[0068] Numerical solution strategy:

[0069] Time discretization: For transient problems, the time derivatives are discretized using time integration methods (e.g., implicit backward Euler method, Crank-Nicolson method, etc.).

[0070] Nonlinear solution: After discretization, the coupled system of equations usually forms a large nonlinear algebraic system of equations R(U,t)=0, where U is a global vector containing all unknown degrees of freedom of all nodes (such as temperature and displacement components). Use Newton's method for iterative solution:

[0071] At the kth iteration, the residual vector R(Uk) is calculated.

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

[0073] Solving linear systems Get the solution increment .

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

[0075] Repeat until the residual meets the convergence criterion.

[0076] Linear solution: The aforementioned linear systems are usually large and sparse, and are solved using preconditioned iterative methods such as the generalized minimum residual method (GMRES) or the biconjugate gradient stabilized method (BiCGStab), combined with efficient preconditioners such as algebraic multigrid (AMG) or incomplete LU decomposition.

[0077] Step (b) implements progressive element removal based on physical criteria.

[0078] 1) Status check and unit marking: Based on the physical field solution obtained in step (a), the active units within the predetermined inspection range are traversed and units that meet the removal conditions are marked as "units to be removed" based on at least one preset physical quantity threshold (for example, the unit stress index exceeds the stress threshold);

[0079] 2) Automatic identification and update of new exposed boundaries: For each cell E_del marked as a "to-be-removed cell", before it is physically removed, all its active neighboring cells N that are not marked as to-be-removed are checked. If the shared internal surface F_internal between E_del and N is not currently an "exposed boundary", F_internal is updated to the new "exposed boundary" and confirmed as the new eroded surface;

[0080] 3) Physical removal of elements: Delete all elements currently marked as "elements to be removed" from the finite element mesh data structure;

[0081] 4) Iterative erosion strategy: Repeat the stress-based process in step 1), step 2), and step 3), focusing the inspection scope on all active elements that are currently "exposed boundaries" until the iteration termination condition is met.

[0082] Specifically, at specific stages of the solution process (e.g., after each time step), the following physically-based element removal process is performed. Two methods are presented: one based on a temperature threshold and the other based on a stress threshold. The temperature-based removal method is used to test the algorithm's functionality and is not suitable for real-world applications. The stress-based removal method, on the other hand, can be used to simulate the progressive loss of real metal materials.

[0083] 1) Status check and marking ( Figure 1a Step B):

[0084] Based on temperature: Traverse the temperature of all cells on the eroded boundary. Calculate the temperature index of these cells (such as the average temperature of the cell, the maximum temperature of the integration point, etc.). If the temperature index exceeds the preset temperature threshold, the cell is marked as "to be removed". Figure 1a shown.

[0085] Stress-based: Traverse the stress of all elements on the eroded boundary. Calculate the stress index of these elements (such as the element average stress, the maximum von Mises stress at the integration point). If the stress index exceeds the preset stress threshold, the element is marked as "to be removed". Figure 1b shown.

[0086] 2) Boundary Update ( Figure 1a and Figure 1b Step C): For each cell E_del marked as "to be removed": traverse all its neighboring cells N. If neighbor N is active and not marked as "to be removed", add the internal shared surface F_internal between cell E_del and neighbor N to the definition of "exposed boundary", thereby confirming the new eroded surface. This step is completed before physically removing the cell to ensure the exposure of the new surface after the simulation material removal.

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

[0088] 4) Iterative erosion (for stress criteria, ( Figure 1bStep E, Step C, Step D)) Repeat the stress-based process in step 1) to step 3). That is, on the newly exposed boundary and in the corresponding unremoved elements on the original exposed boundary, check again whether the stress exceeds the limit, and mark, update the boundary, and remove it. This iterative process simulates the layer-by-layer erosion of the material until no new elements are marked in an iterative cycle or the preset maximum number of iterations is reached. For the temperature criterion, there is no such restriction. When the temperature of the eroded surface is higher than the threshold, it is directly deleted (such as Figure 2a As shown in the figure, the stress criterion is that due to the heat / force load on the material, the material may also have large stress inside. The material deletion conditions are often also satisfied by the internal elements, and the material loss often develops from the eroded surface. Therefore, the retrieval method is limited to layer-by-layer inspection and recursive mesh deletion (such as Figure 2b shown).

[0089] Step (c), island unit detection and processing (optional).

[0090] After performing the removal operation in step (b), starting from the predefined “stable boundary”, all active units connected to the “stable boundary” are identified through a graph search algorithm;

[0091] All units that are active but not identified as connected are judged as “island units”;

[0092] Remove all identified "island elements" from the finite element mesh data structure.

[0093] Specifically, it is recommended to add this connectivity analysis under the stress criterion, ( Figure 1b Step F): After confirming that no unit is marked as out of limit or the preset maximum number of iterations is reached, the "island" unit detection analysis begins, mainly performing connectivity analysis, including:

[0094] Define one or more "stable boundary" identifiers, which represent the boundaries of regions that are fixed or will not be completely removed during the process.

[0095] Start from all active cells adjacent to the "stability boundary".

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

[0097] Identify all cells in the grid that are active but not marked as reachable; these are called "island" cells.

[0098] Remove all identified "island" elements from the finite element mesh data structure.

[0099] Step (d) updates the grid data structure and synchronizes it with the solver. Specifically, the following steps are performed:

[0100] 1) Updating connectivity and status: After performing the unit removal in steps (b) and / or (c), modifying the mesh data structure to reflect the changes in connectivity between nodes and units resulting from the unit removal, as well as the updates to the activity status of the units and nodes;

[0101] 2) Reconstructing adjacency and boundary information: Based on the updated connectivity, recalculate or update the adjacency information between active units, and maintain or reconstruct the geometric entities and their topological information that constitute all currently defined boundaries (including newly generated "exposed boundaries");

[0102] 3) (If parallel computing) Synchronize parallel environment data: Through inter-processor communication, synchronize the geometry and topology information of the shared boundaries between the mesh parts held by each processor due to cell removal to ensure the consistency of global mesh data;

[0103] 4) Notify the main solution framework: Notify the main finite element solution framework through internal mechanisms, indicating that the mesh data structure has changed, to trigger the necessary reinitialization of the solver internal state (such as degree of freedom mapping, sparse matrix mode).

[0104] In order to more clearly illustrate the technical solution of the present invention, the present invention is described in detail below with reference to the accompanying drawings. Figure 1a~Figure 1b The flowcharts of two grid cell removal strategies based on different physical criteria in the method of the present invention are shown. Figure 1a The temperature-based mesh deletion strategy is illustrated. Within a time step, this strategy mainly performs a one-time check: after the coupled field calculation (the beginning of the time step), the unit temperature on the initial exposed boundary is checked (step B). If it exceeds the threshold, the adjacent boundaries are updated (step C) and the marked units are deleted (step D), and then the mesh is updated at the end of the time step. This represents a mode of removing units based on the current temperature field snapshot. This mesh deletion strategy is relatively simple. In actual working environments, the deletion of materials due to temperature is not in line with physical laws. On this basis, a stress-based material deletion strategy is simulated to simulate the progressive damage process of metal materials in a real thermal-mechanical coupling environment. Figure 1bThis diagram illustrates a stress-based mesh deletion strategy. This strategy embodies a more complex, iterative, layer-by-layer removal logic: After the coupled field calculation, the initial exposed boundary stresses are checked and marked (step B). If any elements are marked, the adjacent boundaries are updated (step C), and then the marked elements are deleted (step D). The key point is that after the deletion operation, the element stresses on all currently exposed boundaries (including newly exposed boundaries) are rechecked (step E). If new out-of-limit elements are found, the system returns to step C, forming an iterative cycle of steps C, D, E, and C, simulating the layer-by-layer erosion of the material until no new out-of-limit elements are found. After the entire iterative removal process is completed (regardless of whether any elements are removed), island detection and deletion are performed (step F) to ensure the physical connectivity of the mesh. Finally, the final mesh update is performed at the end of the time step. The core of this strategy is to simulate the stress-driven progressive physical erosion of the material by iteratively checking for newly exposed surfaces.

[0105] Figure 2a~Figure 2b The temperature-based mesh deletion strategy is schematically shown (e.g. Figure 1a Depending on the heat flux distribution (here showing the effect of non-uniform boundaries), the temperature exceeding the limit may cause single or multiple layers of elements to be removed at once, reflecting the instantaneous state of reaching the thermal threshold rather than advancing layer by layer. The corresponding calculation results are shown in Figure 3 As shown (the top is the heat flow boundary, the left and right sides are the adiabatic boundaries, and the bottom is constant temperature), the initial heat flow boundary given here is a function based on the coordinates. When a new exposed boundary appears, a new heat flow boundary is generated based on the coordinates of the new boundary for subsequent calculations. The same method can be extended to the pressure boundary. Subsequently, the heat flux distribution at the initial moment of transmission from the internal and external flow fields in the aerospace field to the interior of the structure is parameterized, and the newly generated boundary is parameterized based on the coordinates to update the heat flow boundary. Updating the boundary based on the coordinates can better conform to the actual environment faced by the new boundary, thereby simulating the loss and retreat of the material.

[0106] Figure 3 Schematic illustration of stress-based mesh deletion strategies (e.g. Figure 1b It emphasizes the iterative logic added within a time step: by checking layer by layer (corresponding to Figure 1b The C -> D -> E cycle), delete cells, update the mesh topology, determine the new exposure boundary through neighborhood search, and recheck whether the cells on these new boundaries exceed the limit, so as to achieve the effect of gradual removal of material from the initial exposure surface to the interior, which is closer to the physical erosion process. The corresponding calculation results are shown in Figure 2. Figure 4a~Figure 4bAs shown in the figure (the top is a constant heat flux and pressure boundary, the left and right sides are adiabatic boundaries, and the bottom is a 300K constant temperature and clamped boundaries). This mesh deletion strategy is mainly based on the consideration that material loss starts from the eroded surface and stress exceedance is not limited to the eroded surface. Therefore, it is necessary to check and delete the layers layer by layer starting from the top.

[0107] It should be noted that Figure 2a~Figure 2b and Figure 3 It is a numerical simulation process carried out under two different calculation boundaries.

[0108] It should be noted that in other embodiments of the present invention, it can be implemented by developing user-defined modules (such as User Objects or subroutines) based on a general finite element analysis software platform, including, for example:

[0109] Data structure: Use a collection to store the IDs of cells to be removed, removed, and newly exposed boundaries.

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

[0111] Module integration: Embed the element removal and island detection modules into the main solver's process control, for example, calling them after each time step and before updating the mesh.

[0112] Information transfer: Ensures that the physics solver has access to the latest mesh information and that the element removal module has access to the values ​​of the physics variables (temperature, stress, etc.) calculated at the current time step.

[0113] Connectivity algorithm: Graph search algorithm (BFS or DFS) is used for island detection.

[0114] The above implementation modes are all within the protection scope of the claims of the present invention.

[0115] The units involved in the embodiments of the present invention may be implemented in software or hardware, and the units described may also be provided in a processor. In some cases, the names of these units do not limit the units themselves.

[0116] According to one aspect of an embodiment of the present invention, a computer program product or computer program is provided, comprising computer instructions stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the computer device to perform the methods provided in the various optional implementations described above.

[0117] As another aspect, 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 independently and not incorporated into the electronic device. The computer-readable medium carries one or more programs, and when executed by the 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 grid processing, characterized in that: The following steps are involved: Establish and solve the coupled physical field model: perform finite element spatial discretization on the solution domain, establish and solve a set of coupled partial differential equations describing 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 unit removal based on physical criteria; Finally, the grid data structure update and solver synchronization are performed; The progressive unit removal based on physical principles specifically includes the following sub-steps: S11, status check and unit marking: based on the obtained numerical solution of the physical field, traverse the active units within the predetermined inspection range, and mark the units that meet the removal conditions as "units to be removed" according to at least one preset physical quantity threshold; S12, identification and update of new exposed boundaries: For each unit E_del marked as "to be removed", before it is physically removed, check all its active neighboring units N that are not marked as to be removed. If the shared internal surface F_internal between E_del and N is not currently an "exposed boundary", then update F_internal to the new "exposed boundary" and confirm it as the new eroded surface; S13, physical removal of elements: deleting all elements currently marked as "elements to be removed" from the finite element mesh data structure; S14, iterative erosion: repeat steps S11 to S13, focusing the inspection range on all active cells that are currently "exposed boundaries" until the iteration termination condition is met; Between the physical-based progressive element removal and the mesh data structure update and solver synchronization, the following steps are included: S21, after the progressive unit removal operation is completed, starting from the predefined "stable boundary", all active units connected to the "stable boundary" are identified through a graph search algorithm; S22, all units that are active but not identified as connected are determined as "island units"; S23, deleting all identified "island elements" from the finite element mesh data structure; In step S11, the status check and unit marking includes the following sub-steps: If the stress-based condition check and element marking is used, the check range is all active elements that are currently adjacent to at least one defined "exposure boundary"; the stress index within each checked element is also calculated.

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

3. The simulation method based on multi-field coupling and dynamic grid processing according to any one of claims 1 to 2, characterized in that: The execution of the grid data structure update and the solver synchronization specifically includes the following sub-steps: S31, updating connection relationships and states: after the unit removal is completed, modifying the grid data structure to reflect the changes in the connection relationships between nodes and units caused by the unit removal, as well as the updates of the activity states of the units and nodes; S32, rebuilding adjacency and boundary information: recalculating or updating adjacency information between active units based on the updated connection relationship, and maintaining or rebuilding geometric entities and their topological information constituting all currently defined boundaries; S33, synchronizing parallel environment data for parallel computing: synchronizing the geometry and topology information of the shared boundaries between the mesh parts held by each processor due to unit removal through inter-processor communication, to ensure the consistency of global mesh data; S34, notifying the main solution framework: notifying the main finite element solution framework through an internal mechanism that the mesh data structure has changed to trigger the reinitialization of the internal state of the solver.

4. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, characterized in that: The coupled partial differential equations are solved by nonlinear iterative Newton method, wherein the linear equations generated in each Newton iteration are solved by preconditioned iterative method.

5. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, characterized in that: In step S14, the iteration termination condition includes at least one of the following: During the status check and unit marking steps of an iteration, no new units are marked as "units to be removed"; or At the beginning of the iteration, the "set of units to be removed" is empty; or The preset maximum number of iterations has been reached.

6. The simulation method based on multi-field coupling and dynamic grid processing according to claim 1, characterized in that: In step S21, the graph search algorithm includes breadth-first search (BFS), depth-first search (DFS), or connectivity analysis method.

7. A simulation system based on multi-field coupling and dynamic grid processing, characterized in that: The method comprises a computer device, wherein the computer device comprises a processor and a memory, wherein a computer program is stored in the memory, and when the computer program is loaded by the processor, the method according to any one of claims 1 to 2 is executed.

Images