Multi-scale modeling and simulation prediction method and device for two-phase composite material

Abstract

The invention discloses a multi-scale modeling and simulation prediction method and device for a two-phase composite material, and the method comprises the following steps: S1, carrying out the molecular dynamics simulation of the penetration process of a solution in the composite material based on molecular dynamics, extracting the microscopic information in the simulation process, and obtaining a microstructure model of the composite material; s2, mapping the microstructure model to a finite element grid, and reserving structures between particles and interfaces to obtain a finite element model; s3, inserting a cohesive zone unit into an interface of the finite element model, and modeling an interface adhesion and separation phenomenon between a matrix and particles of the composite material based on the cohesive zone unit; and S4, carrying out analogue simulation on the whole damage evolution process based on the finite element model to realize mechanical simulation analysis and prediction. According to the method, composite material modeling is carried out based on a real microstructure, and more accurate mechanical simulation analysis is carried out based on the modeling.

Description

Technical Field

[0001] This invention relates to the field of multi-scale three-dimensional modeling technology for two-phase composite materials, and particularly to a multi-scale modeling and simulation prediction method and apparatus for two-phase composite materials. Background Technology

[0002] In dual-phase metallic composites, the mechanical behavior typically exhibits a significant transition from ductile fracture to brittle fracture with increasing brittle phase volume fraction, accompanied by changes in microstructure and composition. This transition in ductile-brittle material is crucial for fracture resistance, damage tolerance, and overall structural reliability, posing a significant challenge to the reliable application of industrial components. However, accurately predicting its initiation and evolution remains a highly challenging task, primarily due to the complex structures and interactions among phase morphology, interfacial bonding, and multi-scale deformation localization.

[0003] Currently, methods for modeling and predicting the mechanical properties of two-phase metal composites often rely on idealized microstructures. While this idealization simplifies the problem to some extent, it ignores the complex microstructural inhomogeneities present in real materials, such as particle size distribution, the diversity of phase interface morphology, and internal defects. These microstructural characteristics have a crucial impact on the mechanical properties of materials, especially their ductile-brittle transition behavior. Traditional modeling and prediction methods have significant limitations in accurately capturing and predicting the initiation and evolution processes.

[0004] Therefore, developing a multi-scale 3D modeling method that can take into account both microstructural details and macroscopic mechanical response, and effectively describe interface behavior, is of great significance for a deeper understanding of the ductile-brittle transition mechanism of biphase metal composites. Summary of the Invention

[0005] In view of this, it is necessary to provide a multi-scale modeling and simulation prediction method and device for two-phase composite materials to effectively solve the technical problem that current modeling of two-phase composite materials does not accurately consider their complex microstructure, thus affecting the simulation analysis of their mechanical properties.

[0006] This invention provides a multi-scale modeling and simulation prediction method for two-phase composite materials, comprising the following steps: Step S1: Based on molecular dynamics, perform molecular dynamics simulation on the solution infiltration process in the composite material, extract microscopic information during the simulation process, and obtain the microstructure model of the composite material. Step S2: Map the microstructure model to a finite element mesh, and retain the structure between particles and interfaces to obtain a finite element model. Discretize the finite element model and refine the local mesh. Step S3: Insert cohesive zone elements into the interface of the finite element model, and model the interface adhesion and delamination phenomenon between the matrix and particles of the composite material based on the cohesive zone elements. At the same time, the bilinear traction separation law is used to describe the cohesive behavior, and the simulation of the crack initiation and propagation process at the interface is completed. Step S4: Simulate the entire damage evolution process based on the finite element model to achieve mechanical simulation analysis and prediction.

[0007] Preferably, step S1 specifically comprises: Based on the initial microstructure state of the composite material, initial model parameters are set, including atom type, atom number, and positional distribution. An initial model of the composite material is constructed based on the initial model parameters. Set the temperature, pressure, and boundary condition parameters for the solution permeation process simulation, set the simulation time step and total simulation time, and perform molecular dynamics simulation based on the initial model in the set simulation environment; After the simulation, the microscopic information of atoms during the solution infiltration process is extracted. Based on the microscopic information, the solution infiltration path, solution infiltration rate and microstructure evolution characteristics are analyzed, and then the microstructure model of the composite material is obtained.

[0008] Preferably, in step S2, the microstructure model is mapped to a finite element mesh, and the structure between particles and interfaces is preserved to obtain the finite element model, specifically as follows: The microstructure model is transformed onto a finite element mesh using a mapping algorithm, and the microscopic information of the microstructure is transformed onto the nodes and elements of the finite element mesh to obtain a finite element model that matches the microstructure model. During the mapping process, particle shape, particle size, particle distribution, interface roughness, and interface curvature features are preserved. The finite element model is then checked and corrected.

[0009] Preferably, in step S2, the finite element model is discretized and local mesh refinement is performed, specifically as follows: The finite element model is discretized using four-node linear tetrahedral elements, and local mesh refinement is performed at particle corners and triple interface connections.

[0010] Preferably, step S3 specifically comprises: Geometric analysis is performed on the interface of the composite material to identify the interface location, shape, and connection method, which serve as the geometric information of the cohesive zone unit. Parameters related to the bilinear traction separation law are set, and boundary conditions and load steps are set to simulate the interaction between the matrix and particles under different working conditions. Based on the relevant parameters and different working conditions, the cohesive zone unit is inserted into the interface to simulate the interface behavior.

[0011] Preferably, step S4 specifically comprises: A vertical displacement-controlled tensile load is applied to the finite element model to simulate the stress under actual working conditions. Symmetrical boundary conditions are applied to the lateral surface to simulate plane strain conditions. The entire damage evolution process is simulated in the simulation environment, including matrix plasticity, ductile fracture, particle cracking, and interface delamination. Then, mechanical simulation analysis is performed on the mechanical response and failure mechanism of the composite material.

[0012] Preferably, the mechanical simulation analysis of the composite material's mechanical response and failure mechanism in step S4 specifically includes: During the simulated damage evolution process, the stress and strain distribution of the model are monitored, and the monitoring data are analyzed to obtain the damage initiation location, development path and final failure mode. Based on the analysis results, it is determined whether the stress concentration area matches the actual vulnerable part and whether the strain change conforms to the damage evolution law.

[0013] Preferably, step S4 further includes: Before conducting simulation analysis, the material parameters of the model are set according to the physical properties of the composite material to ensure that the simulation results reflect the mechanical properties of the material; when applying load, the load size and loading rate are determined according to the actual working conditions.

[0014] Preferably, step S4 further includes: The simulation results are post-processed to visualize the evolution process.

[0015] The present invention also provides a multi-scale modeling and simulation prediction device for two-phase composite materials, including a memory and a processor. The memory stores a computer program, and when the computer program is executed by the processor, it implements the multi-scale modeling and simulation prediction method for the two-phase composite materials.

[0016] Compared with existing technologies, the advantages of this invention are as follows: This invention abandons the idealized model assumptions commonly used in traditional research, and constructs a model based on the microstructural characteristics of actual materials. This allows for the accurate capture of the true morphology and distribution of the matrix and particles in composite materials, as well as the complex interface conditions between them. This highly realistic modeling method makes the simulation results closer to the mechanical behavior of actual materials, providing a reliable basis for accurately predicting the performance of two-phase composite materials under different working conditions, and greatly improving the accuracy and reliability of the predictions. Secondly, mechanical simulation of the three-dimensional model is performed by applying displacement-controlled tensile loads in the vertical direction and symmetrical boundary conditions on the lateral surfaces to accurately simulate plane strain conditions. During the simulation, the entire damage evolution process can be fully presented, covering various situations such as matrix plasticity, ductile fracture, particle cracking, and interface delamination. With the help of this simulation process, the mechanical response and failure mechanisms of two-phase composite materials can be analyzed in depth and comprehensively. Attached Figure Description

[0017] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 is a flowchart of an embodiment of a multi-scale modeling and simulation prediction method for two-phase composite materials provided by the present invention; Figure 2 yes Figure 1 Microstructure and particle size distribution of Cu-W composite materials with increasing tungsten content in the illustrated embodiments; Figure 3 yes Figure 1 The initial and final microstructures of the three Cu-W composite materials W70, W80 and W90 in the illustrated embodiments are shown. Figure 4 yes Figure 1 Microstructure reconstruction and finite element mesh generation diagram of Cu-W composite material in the illustrated embodiment; Figure 5 yes Figure 1 The finite element mesh diagram generated by interpolation in the illustrated embodiment; Figure 6 yes Figure 1 The deformation, strain, and stress field evolution diagrams of the W60 Cu-W composite material under tensile load in the illustrated embodiment are shown. Figure 7 yes Figure 1 The evolution of deformation, strain, and stress field of the W90 Cu-W composite material under tensile load in the illustrated embodiment is shown in the figure. Figure 8a yes Figure 1The diagram shows the transformation of the tensile mechanical properties of Cu-W composite materials in the illustrated embodiment from ductile to brittle-like behavior as the W content increases. Figure 8b yes Figure 1 The diagram shows the evolution of the ultimate tensile strength and elongation at break of the Cu-W composite material with increasing W content in the illustrated embodiment.

[0018] Figure 9a yes Figure 1 Typical stress-strain curves of Cu-W composite materials at different tungsten mass fractions in the embodiments shown; Figure 9b yes Figure 1 The graph shows the functional relationship between the tensile strength and elongation of the Cu-W composite material in the illustrated embodiment and the tungsten content. Detailed Implementation

[0019] Preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, which form part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, but are not intended to limit the scope of the present invention.

[0020] Example 1 Before introducing the embodiments of the present invention, the inventive concept of the method will be explained first.

[0021] Currently, in the 3D modeling of two-phase composite materials, such as Cu-W, the traditional method is to use an idealized or simple geometric model as its microscopic model. However, the problem now is: In the process of three-dimensional modeling of Cu-W composite materials, the molecular dynamics calculation is too small to directly reflect the mechanical properties of the material at the macroscopic scale. On the other hand, macroscopic finite element simulation cannot accurately capture the microstructural features and interactions inside the material, resulting in a large deviation between the simulation results and the actual experimental data. This makes it impossible to accurately predict the mechanical behavior and performance of the two-phase composite material under actual working conditions.

[0022] Within the scope of mechanical analysis of the microstructure and properties of Cu-W composite materials, the prediction of the mechanical properties of Cu-W composite materials is involved.

[0023] While molecular dynamics (MD) simulations can reveal the mechanisms of atomic motion and interaction at the microscopic level of materials, providing atomic-level details for understanding material properties, their computational scale is limited. They are usually only applicable to systems at the nanometer to micrometer scale and are difficult to directly relate to the mechanical behavior of materials at the macroscopic scale.

[0024] The Finite Element Method (FEM) is an important engineering method for data analysis. It can discretize a complex continuum into multiple simple finite elements, solve the mechanical and physical properties of each element through mathematical equations, and then obtain an approximate solution to the overall problem.

[0025] In view of this, please refer to Figure 1 This embodiment provides a multi-scale modeling and simulation prediction method for two-phase composite materials, including the following steps: Step S1: Based on molecular dynamics, perform molecular dynamics simulation on the solution infiltration process in the composite material, extract microscopic information during the simulation process, and obtain the microstructure model of the composite material. Step S2: Map the microstructure model to a finite element mesh, and retain the structure between particles and interfaces to obtain a finite element model. Discretize the finite element model and refine the local mesh. Step S3: Insert cohesive zone elements into the interface of the finite element model, and model the interface adhesion and delamination phenomenon between the matrix and particles of the composite material based on the cohesive zone elements. At the same time, the bilinear traction separation law is used to describe the cohesive behavior, and the simulation of the crack initiation and propagation process at the interface is completed. Step S4: Simulate the entire damage evolution process based on the finite element model to achieve mechanical simulation analysis and prediction.

[0026] Specifically, this embodiment uses LAMMPS software to simulate the melt penetration process in Cu-W composite materials; ABAQUS software is used to map the real microstructure obtained from the molecular dynamics simulation of the melt penetration process onto a finite element mesh, preserving the complex morphology of the W particles and the Cu-W interface; four-node linear tetrahedral elements C3D8 are used for discretization, and local mesh refinement is implemented in high stress concentration areas such as particle corners and triple interfaces; cohesive zone elements COH3D8 inserted at all Cu-W interfaces are used to model the interface debonding phenomenon between the Cu matrix and W particles, and the bilinear traction separation law is used to describe the cohesive behavior, thereby simulating the crack initiation and propagation process at the interface; a displacement-controlled tensile load is applied in the vertical direction, while symmetrical boundary conditions are applied on the lateral surface to simulate plane strain conditions; the simulation fully presents the entire damage evolution process, including matrix plasticity, ductile fracture, particle cracking, and interface delamination, thus enabling a comprehensive analysis of mechanical response and failure mechanisms.

[0027] To better illustrate the method of this invention, Cu-W composite material is used as a specific two-phase composite material to provide a detailed explanation of the multi-scale three-dimensional modeling and mechanical simulation prediction process.

[0028] like Figure 2 The diagram shows the microstructure and particle size distribution of the Cu-W composite material, providing a basis for the subsequent analysis of the microstructure of the Cu-W composite material.

[0029] As a preferred embodiment, step S1 will first be described in detail, such as... Figure 3 As shown, the initial and final microstructures of three Cu-W composite materials were analyzed using molecular dynamics methods and LAMMPS software to simulate the melt infiltration process in the Cu-W composite materials. Specifically, an initial model of the Cu-W composite material was constructed, setting reasonable atomic types, quantities, and initial positional distributions based on the actual material conditions to accurately reflect the initial microstructure state of the composite material. Next, the boundary condition parameters such as temperature and pressure required for the simulation were determined, and the simulation time step and total simulation time were set. The selection of the time step balanced computational efficiency and simulation accuracy, while the total simulation time ensured sufficient observation of the dynamic changes in the melt infiltration process. Then, LAMMPS software was run for molecular dynamics simulation, monitoring changes in key physical quantities such as energy, temperature, and pressure in real time during the simulation to ensure the stability and reliability of the simulation process. After the simulation, information such as the position and velocity of each atom during the melt infiltration process was extracted for subsequent analysis of the melt infiltration path, velocity, and microstructure evolution characteristics.

[0030] After completing step S1, proceed to step S2, as follows: Figure 4 As shown, the microstructure of Cu-W composite materials is remodeled from the atomic scale to the continuous medium scale: 2. Reconstruction of the microstructure and generation of the finite element mesh. The real microstructure obtained from molecular dynamics simulations of the melt infiltration process is mapped to the finite element mesh using ABAQUS software. A finite element model matching the molecular dynamics simulation results is created in ABAQUS software to ensure that the model accurately reflects the geometric features of the real microstructure. Through mapping algorithms or tools, information such as atomic positions and particle morphology obtained from molecular dynamics simulations is accurately converted to the nodes and elements of the finite element mesh, achieving a realistic reproduction of the microstructure. During the mapping process, the complex morphology of the W particles and the Cu-W interface is preserved, including the shape, size, and distribution of the particles, as well as the roughness and curvature of the interface, to ensure the accuracy of subsequent mechanical simulation predictions. Necessary checks and corrections are performed on the mapped finite element model, such as eliminating mesh distortion and optimizing mesh quality, to improve the computational efficiency and stability of the model. Figure 5 As shown, the three-dimensional finite element mesh models of W60 and W90 composite materials after mapping are obtained. Next, in step S2, the discretization process is performed. Four-node linear tetrahedral elements (C3D8) are used to discretize the mapped finite element model, and local mesh refinement is implemented in high-stress concentration areas such as particle corners and triple interface connections. This allows for more accurate simulation of the mechanical behavior of these key areas, improving the accuracy of the simulation results.

[0031] Next, in step S3, the interface debonding phenomenon between the Cu matrix and W particles is modeled using cohesive zone elements (COH3D8) inserted at all Cu-W interfaces. The cohesion behavior is described using the bilinear traction-separation law, thereby simulating the crack initiation and propagation process at the interface. Geometric analysis is performed on the Cu-W interfaces to identify their location, shape, and connectivity, providing geometric information for inserting cohesive zone elements. In ABAQUS software, cohesive zone elements (COH3D8) are inserted at all Cu-W interfaces based on the analysis results to simulate interface behavior. Key parameters of the bilinear traction-separation law are determined based on experimental data or theory; these parameters affect the simulation effect of cohesion behavior. Boundary conditions and load steps are set to simulate the interaction between the Cu matrix and W particles under different operating conditions.

[0032] Finally, in step S4, a vertical displacement-controlled tensile load is applied to the model to ensure the accuracy and stability of the load application. Symmetrical boundary conditions are appropriately applied to the lateral surfaces to accurately simulate plane strain conditions, providing a reliable foundation for subsequent analysis. The vertical displacement-controlled tensile load is applied to the model to ensure the accuracy and stability of the load application, simulating the stress conditions under actual working conditions. The simulation presents a detailed view of the complete damage evolution process, focusing not only on the plastic changes of the matrix but also comprehensively capturing phenomena such as ductile fracture, particle crack generation, and interface delamination, thereby conducting an in-depth and comprehensive analysis of mechanical response and failure mechanisms. Figure 6 , Figure 7 As shown, Figure 6 The stress-strain diagrams are for W60 and W90 composite materials. Figure 7 The diagram shows the structural distribution of W60 and W90 composite materials.

[0033] In step S4, mechanical simulation analysis is performed on the mechanical response and failure mechanism of composite materials. Specifically, during the simulated damage evolution process, the stress and strain distribution of the model are monitored in real time. By analyzing the data, the damage initiation location, development path and final failure mode are understood. For example, it is observed whether the stress concentration area matches the actual vulnerable part and whether the strain change conforms to the damage evolution law.

[0034] The physical models involved in performing mechanical simulations: in For the gradient of the deformation tensor, This is a deformation tensor, applicable to large deformations, where X is the material coordinate, representing its position within the reference configuration. These are spatial coordinates, representing the position within the current configuration. Let be the displacement vector, representing the displacement from the reference configuration to the current configuration. Let I be the second-order unit tensor. The gradient operator is the reference coordinate X.

[0035] in Represents the stress tensor. For volume forces, For density, For acceleration; For strain tensor, For volume, For variational operators, For surface traction force; Let be the strain energy density function. Shear modulus For the right Cauchy-Green tensor, It is a volume ratio. This is the first Lamé constant.

[0036] in For shape functions, Let i be the displacement component of the i-th node. Here is the stiffness matrix. This is the global nodal displacement vector. The residual vector of the k-th iteration , .

[0037] Specifically, before performing mechanical simulation analysis in step S4, the material parameters of the model, such as elastic modulus, Poisson's ratio, and yield strength, need to be accurately set according to the actual physical properties of the Cu-W composite material to ensure that the simulation results reflect the mechanical properties of the material. When applying loads, in addition to vertical displacement tensile loads and lateral symmetric boundary conditions, the load magnitude and loading rate should be determined according to the actual working conditions, as the loading rate will affect the damage evolution and should be set reasonably to simulate the material response under different loading conditions.

[0038] Specifically, after the simulation in step S4 is completed, the simulation results are post-processed, and the damage evolution process is presented as images or animations using visualization technology, providing a basis for subsequent mechanical response and failure mechanism analysis. The plastic deformation area of ​​the matrix, the initiation point and propagation direction of ductile fracture, the location of particle cracks, and the interface delamination can be seen intuitively.

[0039] like Figure 8a , Figure 8b The figures show the evolution of tensile mechanical properties of Cu-W composite materials, specifically the transition from ductile to brittle behavior via stress-strain curves and the evolution of tensile mechanical properties via ultimate tensile strength and elongation at break. With increasing tungsten content, the material's plasticity gradually decreases, and the failure mode becomes more abrupt. The W60 composite exhibits typical ductile behavior, with an ultimate tensile strength of approximately 510 MPa and an elongation at break of approximately 15%, while also possessing a clear yield point and a significant strain hardening region. With increasing tungsten content, the deformation capacity after yielding gradually weakens. The W70 sample has an ultimate tensile strength of approximately 580 MPa and an elongation at break of approximately 9%, while the W80 sample has an ultimate tensile strength of approximately 625 MPa and an elongation at break of approximately 7%, maintaining moderate ductility, but the stress drop after peak load is steeper. The W90 composite shows significantly reduced strength and ductility, with an ultimate tensile strength of approximately 590 MPa and an elongation at break of less than 5%. Its stress-strain curve exhibits limited strain hardening, and the load-bearing capacity decreases sharply after yielding, indicating a clear transition towards brittle fracture behavior.

[0040] like Figure 9a , Figure 9b The figures show the stress-strain diagram and the functional relationship between tensile strength and elongation as a function of tungsten content for the Cu-W composite material obtained after simulation. The stress-strain curve of W60 exhibits a significant work hardening stage, followed by necking at approximately 8% strain and eventual failure, which is in perfect agreement with the experimental trend. The calculated values ​​for W60 and W90 are in high agreement with the measured values, verifying the accuracy of the multi-scale model in capturing the trade-off between strength and ductility.

[0041] Example 2 This embodiment provides a multi-scale modeling and simulation prediction device for two-phase composite materials, including a memory and a processor. The memory stores a computer program, and when the computer program is executed by the processor, it implements the multi-scale modeling and simulation prediction method for two-phase composite materials described in Embodiment 1.

[0042] The multi-scale modeling and simulation prediction device for two-phase composite materials provided in this embodiment is used to realize the multi-scale modeling and simulation prediction method for two-phase composite materials. Therefore, the multi-scale modeling and simulation prediction device for two-phase composite materials also has the technical effects of the value method, which will not be described in detail here.

[0043] The above description is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in the present invention should be included within the scope of the present invention.

Claims

1. A multi-scale modeling and simulation prediction method for two-phase composite materials, characterized in that, Includes the following steps: Step S1: Based on molecular dynamics, perform molecular dynamics simulation on the solution infiltration process in the composite material, extract microscopic information during the simulation process, and obtain the microstructure model of the composite material. Step S2: Map the microstructure model to a finite element mesh, and retain the structure between particles and interfaces to obtain a finite element model. Discretize the finite element model and refine the local mesh. Step S3: Insert cohesive zone elements into the interface of the finite element model, and model the interface adhesion and delamination phenomenon between the matrix and particles of the composite material based on the cohesive zone elements. At the same time, the bilinear traction separation law is used to describe the cohesive behavior, and the simulation of the crack initiation and propagation process at the interface is completed. Step S4: Simulate the entire damage evolution process based on the finite element model to achieve mechanical simulation analysis and prediction.

2. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, Step S1 specifically involves: Based on the initial microstructure state of the composite material, initial model parameters are set, including atom type, atom number, and positional distribution. An initial model of the composite material is constructed based on the initial model parameters. Set the temperature, pressure, and boundary condition parameters for the solution permeation process simulation, set the simulation time step and total simulation time, and perform molecular dynamics simulation based on the initial model in the set simulation environment; After the simulation, the microscopic information of atoms during the solution infiltration process is extracted. Based on the microscopic information, the solution infiltration path, solution infiltration rate and microstructure evolution characteristics are analyzed, and then the microstructure model of the composite material is obtained.

3. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, In step S2, the microstructure model is mapped to a finite element mesh, and the structure between particles and interfaces is preserved to obtain the finite element model, specifically as follows: The microstructure model is transformed onto a finite element mesh using a mapping algorithm, and the microscopic information of the microstructure is transformed onto the nodes and elements of the finite element mesh to obtain a finite element model that matches the microstructure model. During the mapping process, particle shape, particle size, particle distribution, interface roughness, and interface curvature features are preserved. The finite element model is then checked and corrected.

4. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, In step S2, the finite element model is discretized and local mesh refinement is performed, specifically as follows: The finite element model is discretized using four-node linear tetrahedral elements, and local mesh refinement is performed at particle corners and triple interface connections.

5. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, Step S3 specifically involves: Geometric analysis is performed on the interface of the composite material to identify the interface location, shape, and connection method, which serve as the geometric information of the cohesive zone unit. Parameters related to the bilinear traction separation law are set, and boundary conditions and load steps are set to simulate the interaction between the matrix and particles under different working conditions. Based on the relevant parameters and different working conditions, the cohesive zone unit is inserted into the interface to simulate the interface behavior.

6. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, Step S4 specifically involves: A vertical displacement-controlled tensile load is applied to the finite element model to simulate the stress under actual working conditions. Symmetrical boundary conditions are applied to the lateral surface to simulate plane strain conditions. The entire damage evolution process is simulated in the simulation environment, including matrix plasticity, ductile fracture, particle cracking, and interface delamination. Then, mechanical simulation analysis is performed on the mechanical response and failure mechanism of the composite material.

7. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 6, characterized in that, In step S4, mechanical simulation analysis is performed on the mechanical response and failure mechanism of the composite material, specifically as follows: During the simulated damage evolution process, the stress and strain distribution of the model are monitored, and the monitoring data are analyzed to obtain the damage initiation location, development path and final failure mode. Based on the analysis results, it is determined whether the stress concentration area matches the actual vulnerable part and whether the strain change conforms to the damage evolution law.

8. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, Step S4 further includes: Before conducting simulation analysis, the material parameters of the model are set according to the physical properties of the composite material to ensure that the simulation results reflect the mechanical properties of the material; when applying load, the load size and loading rate are determined according to the actual working conditions.

9. The multi-scale modeling and simulation prediction method for two-phase composite materials according to claim 1, characterized in that, Step S4 further includes: The simulation results are post-processed to visualize the evolution process.

10. A multi-scale modeling and simulation prediction device for two-phase composite materials, characterized in that, It includes a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, it implements the multi-scale modeling and simulation prediction method for two-phase composite materials as described in any one of claims 1-9.

Images