A Mesoscopic Damage Numerical Simulation Method, Device and Equipment for a Magnesium Matrix Composite

Abstract

The present invention discloses a mesoscopic damage numerical simulation method, device and equipment for a magnesium-based composite material. A magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes is established. The magnesium-based composite material model includes a carbon nanotube model and a magnesium-based composite material matrix model; cohesive elements for replacing the interface phase of the magnesium-based composite material are established on the surface of the carbon nanotube model to obtain a carbon nanotube model with an interface; periodic boundary conditions and load boundary conditions are applied to the magnesium-based composite material RVE grid model, and contact properties and constraint boundary conditions between the carbon nanotube model with an interface and the magnesium-based composite material matrix model are given to obtain a magnesium-based composite material RVE grid model including boundary conditions; the material properties and boundary conditions of the magnesium-based composite material are input into the mesoscopic damage numerical simulation analysis model of the magnesium-based composite material, and the mesoscopic damage numerical simulation results are calculated. The present invention can solve the problem of insufficient research on the mesoscopic damage mechanism of magnesium-based composite materials.

Description

Technical Field

[0001] The present invention belongs to the technical field of finite element analysis of composite materials, and particularly relates to a mesoscopic damage numerical simulation method, device and equipment for magnesium-based composite materials. Background Art

[0002] Magnesium-based composite materials have the advantages of low density, high specific strength, good recyclability, etc., and are widely used in new energy vehicles, manned spacecraft, medical devices, electronic packaging and other fields. However, the service environment of magnesium-based composite materials is harsh and the load conditions are variable, resulting in complex damage and failure forms of composite materials, which are difficult to study.

[0003] At present, for the research on the mesoscopic damage mechanism of composite materials, for example, Patent Application No. CN113971354A provides a finite element method for studying the bending fracture damage of nanocomposites. However, it only targets the two-dimensional model of composite materials and the reinforcement is particles, and it cannot effectively apply to the three-dimensional real structure and damage mechanism of composite materials, and is not applicable to nanoscale fiber reinforcements. Due to the complex structure of composite materials and the existence of the interface between the reinforcement and the matrix, the structure of composite materials changes greatly, and the interface damage form has a great influence on the failure of composite materials. However, the current exploration of the mesoscopic damage mechanism of magnesium-based composite materials is not clear. Studying the mesoscopic damage mechanism of magnesium-based composite materials is crucial for the design and performance regulation of magnesium-based composite materials. There is an urgent need for a mesoscopic damage numerical simulation method for magnesium-based composite materials to solve this problem. Summary of the Invention

[0004] Aiming at the problems existing in the prior art, the present invention provides a mesoscopic damage numerical simulation method, device and equipment for magnesium-based composite materials, which can solve the problem of insufficient research on the mesoscopic damage mechanism of magnesium-based composite materials and improve the prediction accuracy of the numerical simulation method.

[0005] In order to solve the above technical problems, the present invention is realized through the following technical solutions:

[0006] A mesoscopic damage numerical simulation method for magnesium-based composite materials, comprising:

[0007] Establish a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes, where the magnesium-based composite material model includes a carbon nanotube model and a magnesium-based composite material matrix model;

[0008] Establish cohesive elements on the surface of the carbon nanotube model to replace the interface phase of the magnesium-based composite material, and obtain a carbon nanotube model with an interface;

[0009] Perform Boolean operations and assembly on the carbon nanotube model with an interface and the magnesium-based composite material matrix model to obtain a magnesium-based composite material RVE model;

[0010] Correspondingly assign material properties to the matrix phase, interface phase, and reinforcement phase in the RVE model of the magnesium-based composite material, and assign a matrix damage constitutive model and a cohesive constitutive criterion to the RVE model of the magnesium-based composite material to obtain an RVE model of the magnesium-based composite material containing material properties;

[0011] Perform mesh division on the RVE model of the magnesium-based composite material containing material properties to obtain an RVE mesh model of the magnesium-based composite material;

[0012] Apply periodic boundary conditions and load boundary conditions to the RVE mesh model of the magnesium-based composite material, and assign contact properties and constraint boundary conditions between the carbon nanotube model with an interface and the matrix model of the magnesium-based composite material to obtain an RVE mesh model of the magnesium-based composite material containing boundary conditions;

[0013] Assign a finite element analysis mode to the RVE mesh model of the magnesium-based composite material containing boundary conditions to obtain a mesoscopic damage numerical simulation analysis model of the magnesium-based composite material;

[0014] Input the material properties and boundary conditions of the magnesium-based composite material obtained from experiments into the mesoscopic damage numerical simulation analysis model of the magnesium-based composite material, and calculate the mesoscopic damage numerical simulation results, where the numerical simulation results include a stress field, a strain field, and a damage field.

[0015] Furthermore, the establishment of the model of the magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes includes:

[0016] Write a script for generating a model of the magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes based on the random adsorption sequence algorithm;

[0017] Assign the volume ratio of carbon nanotubes to the script,

[0018] Assign the type of the magnesium-based composite material model to be a 3D type RVE to the script;

[0019] Assign the basic dimensions of the magnesium-based composite material model to the script, where the basic dimensions include the side length dimension of the matrix model of the magnesium-based composite material and the length and diameter of the carbon nanotube model;

[0020] Run the script after assigning parameters to generate a model of the magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes.

[0021] Furthermore, use the Abaqus finite element simulation software to run the script after assigning parameters to generate a model of the magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes.

[0022] Furthermore, the thickness of the cohesive element used to replace the interface phase of the magnesium matrix composite material is 0.

[0023] Furthermore, in the representative volume element (RVE) model of the magnesium matrix composite material containing material properties, the material properties include:

[0024] The material properties, failure mode, damage judgment criterion, and damage evolution of the magnesium matrix composite matrix model;

[0025] The material properties, bilinear separation criterion, and parameters of the cohesive element of the interface phase. The parameters of the cohesive element include the interface phase thickness, surface force, and interface separation parameter;

[0026] The material properties and damage mode of the carbon nanotube model.

[0027] Furthermore, the meshing of the RVE model of the magnesium matrix composite material containing material properties is as follows:

[0028] The magnesium matrix composite matrix model and the carbon nanotube model in the RVE model of the magnesium matrix composite material containing material properties are respectively meshed using three-node tetrahedral elements;

[0029] The interface phase cohesive element in the RVE model of the magnesium matrix composite material containing material properties is meshed using six-node three-dimensional cohesive elements.

[0030] Furthermore, applying periodic boundary conditions to the RVE grid model of the magnesium matrix composite material includes:

[0031] Writing a script based on the Python language to apply periodic boundary conditions to the RVE grid model of the magnesium matrix composite material;

[0032] The script couples all the grid nodes on the opposite faces, opposite edges, and opposite vertices of the RVE grid model of the magnesium matrix composite material to complete the application of periodic boundary conditions;

[0033] Applying load boundary conditions to the RVE grid model of the magnesium matrix composite material includes:

[0034] Applying load boundary conditions to the RVE grid model of the magnesium matrix composite material using a surface load on the upper surface of the Y-axis. The load direction is the positive direction of the Y-axis, and the load mode is a smooth analysis step.

[0035] Furthermore, the finite element analysis mode is the dynamic explicit analysis mode.

[0036] A mesoscopic damage numerical simulation device for magnesium matrix composite materials, comprising:

[0037] The first establishment module is used to establish a magnesium-based composite material model with a random distribution of three-dimensional carbon nanotubes. The magnesium-based composite material model includes a carbon nanotube model and a magnesium-based composite material matrix model;

[0038] The second establishment module is used to establish cohesive elements on the surface of the carbon nanotube model to replace the interface phase of the magnesium-based composite material, obtaining a carbon nanotube model with an interface;

[0039] The assembly module is used to perform Boolean operations and assembly on the carbon nanotube model with an interface and the magnesium-based composite material matrix model, obtaining a representative volume element (RVE) model of the magnesium-based composite material;

[0040] The first assignment module is used to respectively assign corresponding material properties to the matrix phase, interface phase, and reinforcement phase in the magnesium-based composite material RVE model, and assign a matrix damage constitutive model and a cohesive constitutive criterion to the magnesium-based composite material RVE model, obtaining a magnesium-based composite material RVE model containing material properties;

[0041] The mesh generation module is used to perform mesh generation on the magnesium-based composite material RVE model containing material properties, obtaining a magnesium-based composite material RVE mesh model;

[0042] The second assignment module is used to apply periodic boundary conditions and load boundary conditions to the magnesium-based composite material RVE mesh model, and assign contact properties and constraint boundary conditions between the carbon nanotube model with an interface and the magnesium-based composite material matrix model, obtaining a magnesium-based composite material RVE mesh model containing boundary conditions;

[0043] The third assignment module is used to assign a finite element analysis mode to the magnesium-based composite material RVE mesh model containing boundary conditions, obtaining a mesoscopic damage numerical simulation analysis model of the magnesium-based composite material;

[0044] The calculation module is used to input the material properties and boundary conditions of the magnesium-based composite material obtained from experiments into the mesoscopic damage numerical simulation analysis model of the magnesium-based composite material, and calculate the mesoscopic damage numerical simulation results. The numerical simulation results include a stress field, a strain field, and a damage field.

[0045] A device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of a mesoscopic damage numerical simulation method of a magnesium-based composite material are implemented.

[0046] Compared with the prior art, the present invention has at least the following beneficial effects:

[0047] The present invention provides a mesoscopic damage numerical simulation method for magnesium matrix composites. The established magnesium matrix composite model with randomly distributed three-dimensional carbon nanotubes is more realistic and reliable, and can accurately explore the microscopic deformation behavior of the matrix composite compared with the traditional single-cell RVE model; the cohesive unit is used to replace the interface phase of the magnesium matrix composite, and the initial damage occurrence time of the magnesium matrix composite can be well predicted by regulating the interface bonding strength; the periodic boundary conditions and load boundary conditions are applied to the RVE grid model of the magnesium matrix composite, and the contact properties and constraint boundary conditions between the carbon nanotube model with interface and the magnesium matrix composite matrix model are given to obtain the RVE grid model of the magnesium matrix composite containing boundary conditions, so that the model can better reflect the actual damage mechanism of the magnesium matrix composite. The method provided by the present invention provides a reference for the mesoscopic scale finite element analysis of magnesium matrix composites, and can provide reliable theoretical guidance for the design, performance regulation and damage behavior monitoring of composites in practical applications.

[0048] In order to make the above objects, features and advantages of the present invention more obvious and understandable, the following specifically enumerates preferred embodiments and, in conjunction with the accompanying drawings, the detailed description is as follows. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions in the specific embodiments of the present invention, the following will briefly introduce the drawings required for the description of the specific embodiments. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0050] Figure 1 It is a flowchart of a mesoscopic damage numerical simulation method for a magnesium matrix composite of the present invention.

[0051] Figure 2 It is a schematic diagram of a three-dimensional RVE model and a grid model established by using finite element software in the present invention, where:

[0052] (a) Schematic diagram of the three-dimensional model of the magnesium matrix composite;

[0053] (b) Schematic diagram of the grid model of the magnesium matrix composite.

[0054] Figure 3 It is the constitutive model of the magnesium matrix composite of the present invention, where:

[0055] (a) Constitutive model of the matrix material of the magnesium matrix composite;

[0056] (b) Cohesive constitutive model of the interface phase of the magnesium matrix composite.

[0057] Figure 4 It is a schematic diagram of the simulation results and experimental results of the present invention.

[0058] Figure 5 This is the schematic diagram of the final fracture of the magnesium matrix composite material RVE model in the present invention, where:

[0059] (a) Tensile stress of the magnesium matrix composite material;

[0060] (b) von Mises stress of the interface phase;

[0061] (c) Tensile strain of the magnesium matrix composite material;

[0062] (d) Tensile stress of carbon nanotubes in the magnesium matrix composite material.

[0063] Figure 6 This is the schematic diagram of the fracture mechanism of the magnesium matrix composite material in the present invention. Specific embodiments

[0064] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0065] As a specific embodiment of the present invention, in combination with Figure 1 shown in the figure, a mesoscopic damage numerical simulation method for a magnesium matrix composite material specifically includes the following steps:

[0066] S1. Establish a magnesium matrix composite material model with randomly distributed three-dimensional carbon nanotubes. The magnesium matrix composite material model includes a carbon nanotube model and a magnesium matrix composite material matrix model.

[0067] Optionally, the establishment of the magnesium matrix composite material model with randomly distributed three-dimensional carbon nanotubes is as follows:

[0068] S1.1. Write a script for generating a magnesium matrix composite material model with randomly distributed three-dimensional carbon nanotubes based on the random adsorption sequence algorithm;

[0069] S1.2. Assign the volume ratio of carbon nanotubes to the script,

[0070] S1.3. Assign the type of the magnesium matrix composite material model as a 3D type RVE to the script;

[0071] S1.4. Assign the basic dimensions of the magnesium matrix composite material model to the script. The basic dimensions include the side length dimension of the magnesium matrix composite material matrix model and the length and diameter of the carbon nanotube model;

[0072] S1.5. Run the script with the assigned parameters to generate a model of magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes.

[0073] Preferably, use Abaqus finite element simulation software to run the script with the assigned parameters to generate a model of magnesium-based composite material with randomly distributed three-dimensional carbon nanotubes.

[0074] S2. Establish cohesive elements on the surface of the carbon nanotube model to replace the interface phase of the magnesium-based composite material, obtaining a carbon nanotube model with an interface.

[0075] Preferably, the thickness of the cohesive elements used to replace the interface phase of the magnesium-based composite material is 0.

[0076] S3. Perform Boolean operations and assembly on the carbon nanotube model with an interface and the magnesium-based composite material matrix model to obtain a representative volume element (RVE) model of the magnesium-based composite material.

[0077] Specifically, in the magnesium-based composite material RVE model containing material properties, the material properties specifically include the following:

[0078] The material properties, failure modes, damage judgment criteria, and damage evolution of the magnesium-based composite material matrix model;

[0079] The material properties, bilinear separation criterion, and parameters of the cohesive elements of the interface phase. The parameters of the cohesive elements include the interface phase thickness, surface force, and interface separation parameter;

[0080] The material properties and damage modes of the carbon nanotube model.

[0081] S4. Assign corresponding material properties to the matrix phase, interface phase, and reinforcement phase in the magnesium-based composite material RVE model, and assign a matrix damage constitutive model and a cohesive constitutive criterion to the magnesium-based composite material RVE model to obtain a magnesium-based composite material RVE model containing material properties.

[0082] S5. Mesh the magnesium-based composite material RVE model containing material properties to obtain a meshed model of the magnesium-based composite material RVE.

[0083] Optionally, the meshing of the magnesium-based composite material RVE model containing material properties is specifically as follows:

[0084] Use three-node tetrahedral elements (C3D4) to mesh the magnesium-based composite material matrix model and the carbon nanotube model in the magnesium-based composite material RVE model containing material properties respectively;

[0085] The cohesive elements in the interface phase of the magnesium matrix composite RVE model containing material properties are meshed using six-node three-dimensional cohesive elements (COH3D6).

[0086] S6. Periodic boundary conditions and load boundary conditions are applied to the magnesium matrix composite RVE mesh model, and contact properties and constraint boundary conditions are assigned between the carbon nanotube model with an interface and the magnesium matrix composite matrix model to obtain a magnesium matrix composite RVE mesh model containing boundary conditions.

[0087] Optionally, the application of periodic boundary conditions to the magnesium matrix composite RVE mesh model is as follows:

[0088] A script for applying periodic boundary conditions to the magnesium matrix composite RVE mesh model is written based on the Python language.

[0089] All mesh nodes on the opposite faces, opposite edges, and opposite vertices of the magnesium matrix composite RVE mesh model are coupled by the script to complete the application of periodic boundary conditions.

[0090] A three-dimensional RVE model of fiber randomly distributed reinforced composites is implemented using a script written in Python. With periodic structures, periodic meshes, and periodic boundary conditions, the model can better reflect the actual damage mechanism of the composite material.

[0091] Optionally, the application of load boundary conditions to the magnesium matrix composite RVE mesh model is as follows:

[0092] Load boundary conditions are applied to the magnesium matrix composite RVE mesh model using a surface load on the upper surface of the Y-axis, with the load direction being the positive Y-axis and the load mode being a smooth analysis step.

[0093] S7. A finite element analysis mode is assigned to the magnesium matrix composite RVE mesh model containing boundary conditions to obtain a mesoscopic damage numerical simulation analysis model of the magnesium matrix composite.

[0094] Preferably, the finite element analysis mode is a dynamic explicit analysis mode.

[0095] S8. The material properties and boundary conditions of the experimentally obtained magnesium matrix composite are input into the mesoscopic damage numerical simulation analysis model of the magnesium matrix composite, and mesoscopic damage numerical simulation results are calculated. The numerical simulation results include a stress field, a strain field, and a damage field.

[0096] Specifically, the damage parameters and material properties of the magnesium matrix composite obtained through uniaxial tensile experiments of the magnesium matrix composite, and the boundary conditions are the load conditions and constraint conditions of the uniaxial tensile experiment.

[0097] The following introduces a mesoscopic damage numerical simulation method for a magnesium-based composite material provided by the present invention in combination with a specific implementation case. The specific steps are as follows:

[0098] (1) Based on the random sequential adsorption algorithm, write a script for generating a three-dimensional model of a magnesium-based composite material with randomly distributed carbon nanotubes to realize a three-dimensional model of a composite material reinforced by randomly distributed carbon nanotubes (as shown in Figure 2 Figure a), and the specific steps are as follows.

[0099] ① In the Abaqus finite element simulation software, assign parameters to the script. The parameters include the volume ratio of the reinforcement of 0.4237 vol.%, which is equivalent to the addition amount of carbon nanotubes in the experimental material of 0.5% (wt.%), and the type of the magnesium-based composite material model is 3D type RVE;

[0100] ② Assign the basic dimensions of the model: the size of the RVE three-dimensional model is 20 μm × 20 μm × 20 μm; the diameter of the reinforcement model is 0.3779 μm and the length is 3.779 μm;

[0101] ③ Run the script with the assigned parameters to generate a three-dimensional model of a magnesium-based composite material with randomly distributed carbon nanotubes.

[0102] (2) RVE modeling: Establish cohesive elements (0-thickness) on the interface between the carbon nanotube model and the magnesium-based composite material matrix model to simulate the interface bonding. Perform a Boolean operation and assembly on the carbon nanotube model with the interface and the magnesium-based composite material matrix model to obtain a magnesium-based composite material RVE model;

[0103] (3) Assign material properties to the matrix phase, interface phase, and fiber reinforcement in the magnesium-based composite material RVE model established in step (1), and select appropriate damage evolution criteria. The interface phase adopts a bilinear cohesive model, and the BK criterion is used as the failure criterion for the interface region. Define the traction force and fracture energy of the cohesive element, as shown in Figure 3 Figure. Obtain a magnesium-based composite material RVE model containing material properties.

[0104] The specific parameters are as follows:

[0105] Table 1

[0106] Material parameters of Mg–6Zn alloy and CNTs.

[0107]

[0108] Table 2

[0109] Cohesive zone model parameters.

[0110]

[0111] (4) Mesh the magnesium matrix composite RVE model containing material properties in step (3). Use the finite element method based on four-node linear tetrahedral elements (C3D4) to mesh the continuum (magnesium matrix composite matrix model and carbon nanotube model). Use four-node three-dimensional cohesive elements (COH3D6) to establish the interface layer cohesive element mesh for establishing the interface phase. To verify the reliability of the simulation results, a large number of mesh elements are set around the carbon nanotubes, as shown in Figure 2 Figure b. The RVE model is established using 246,935 elements and 44,148 nodes respectively. Among them, there are 241,473 units of C3D4 and 5462 units of COH3D6.

[0112] (5) For the magnesium matrix composite RVE mesh model obtained in step (4), apply periodic boundary conditions using a script, couple all mesh nodes on the opposite faces, opposite edges, and opposite vertices of the model, so that the RVE model can reflect the stress and strain of the macroscopic model.

[0113] According to the actual working conditions and loads, apply load boundary conditions to the magnesium matrix composite RVE mesh model obtained in step (4). Apply surface pressure loading on the upper surface of the Y-axis, and the load is applied through a smooth analysis step. Obtain the magnesium matrix composite RVE mesh model containing boundary conditions.

[0114] (6) For the magnesium matrix composite RVE mesh model containing boundary conditions obtained in step (5), assign the dynamic explicit analysis mode, and assign the total analysis step length to be 2.5 s. Obtain the mesoscopic damage numerical simulation analysis model of the magnesium matrix composite.

[0115] (7) For the mesoscopic damage numerical simulation analysis model of the magnesium matrix composite obtained in step (6), create a job and submit it to the computer for analysis. Using the calculation function of the finite element software, the stress-strain curve of the composite material under the same working conditions is output. Comparing with the experimental results, good consistency is found, verifying the effectiveness of this finite element simulation method, as shown in Figure 4 Figure.

[0116] (8) Analyze the fracture situation of the composite material model obtained in step (7), and explore the damage initiation position, damage mode, and damage evolution of the composite material. Explore the effects of the introduction of carbon nanotube reinforcements and interface phases on the mechanical properties, load transfer, deformation behavior, and damage mechanism of magnesium matrix composites, as follows:

[0117] The simulation and experimental results show that, as Figure 5 shown, the initial damage of the composite material occurs near the carbon nanotubes. As the tensile load increases, the stress concentration at the ends of the carbon nanotubes increases, and interfacial debonding is the main form of crack propagation in the carbon nanotube / magnesium composite material. The carbon nanotubes can blunt the crack tip and change the propagation direction of the main crack, thus delaying fracture. The carbon nanotubes play a dominant role in controlling the fracture mode of the composite material. As Figure 6 shown, the pull-out and bridging of the carbon nanotubes are the main strengthening mechanisms of the composite material, which is beneficial to improving the strength and toughness of the composite material.

[0118] In summary: The mesoscopic damage numerical simulation method of a magnesium-based composite material according to the present invention records the damage and fracture responses of the carbon nanotube / magnesium composite material under uniaxial tension by combining experiments and numerical simulations. In order to conduct mesoscopic mechanical simulations, a typical unit cell model composed of a random distribution model of carbon nanotubes and cohesive units in the carbon nanotube / magnesium interface phase is constructed to describe the packing characteristics of the fibers. It is found that local interface failure occurs in the area around the carbon nanotubes. When the crack extends to the carbon nanotubes, the crack will deflect along the wall surface of the carbon nanotubes, which will delay the crack propagation. Among the foreseeable results, the carbon nanotubes are the main barriers to preventing crack propagation. At the same time, some carbon nanotubes exposed on the fracture surface endow it with the functions of nanofiber pull-out and bridging. The experimental results are highly consistent with the simulation calculation results, verifying the reliability of the numerical simulation method involved in the present invention.

[0119] In an embodiment of the present invention, a computer device is provided. The computer device includes a processor and a memory. The memory is used to store a computer program. The computer program includes program instructions. The processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, and is suitable for implementing one or more instructions. Specifically, it is suitable for loading and executing one or more instructions to implement the corresponding method flow or corresponding function. The processor described in the embodiment of the present invention can be used to implement the operations of a mesoscopic damage numerical simulation method of a magnesium-based composite material.

[0120] In one embodiment of the present invention, when a mesoscopic damage numerical simulation method for a magnesium-based composite material is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above-mentioned embodiment methods of the present invention, it can also be completed by instructing relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above-mentioned various method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file, or some intermediate form, etc. The computer-readable storage medium includes permanent and non-permanent, removable and non-removable media, and information storage can be achieved by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data.

[0121] The computer storage medium can be any available medium or data storage device accessible by a computer, including but not limited to magnetic memory (such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc.), optical memory (such as CDs, DVDs, BDs, HVDs, etc.), and semiconductor memory (such as ROM, EPROM, EEPROM, non-volatile memory (NANDFLASH), solid-state drives (SSD)), etc.

[0122] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-readable storage media (including but not limited to disk memory, CD-ROM, optical memory, etc.) containing computer-usable program code.

[0123] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products of the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of the processes and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for implementing the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0124] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction device that implements the functions specified in one or more of the processes Figure 1 and / or boxes Figure 1 specified in one or more of the boxes.

[0125] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one or more of the processes Figure 1 and / or boxes Figure 1 specified in one or more of the boxes.

[0126] Finally, it should be noted that the above-described embodiments are only specific embodiments of the present invention, used to illustrate the technical solutions of the present invention, rather than limiting it. The protection scope of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: any person skilled in the art within the technical scope disclosed by the present invention can still modify the technical solutions described in the foregoing embodiments or easily conceive of changes, or perform equivalent replacements on some of the technical features; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A mesoscopic damage numerical simulation method for a magnesium-based composite material, characterized in that, it includes: Establish a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes. The magnesium-based composite material model includes a carbon nanotube model and a magnesium-based composite material matrix model; Establish cohesive elements on the surface of the carbon nanotube model to replace the interface phase of the magnesium-based composite material, and obtain a carbon nanotube model with an interface; Perform a Boolean operation and assembly on the carbon nanotube model with an interface and the magnesium-based composite material matrix model to obtain a representative volume element (RVE) model of the magnesium-based composite material; Assign corresponding material properties to the matrix phase, interface phase, and reinforcement phase in the RVE model of the magnesium-based composite material, and assign a matrix damage constitutive model and a cohesive constitutive criterion to the RVE model of the magnesium-based composite material to obtain an RVE model of the magnesium-based composite material containing material properties; Perform mesh division on the RVE model of the magnesium-based composite material containing material properties to obtain an RVE mesh model of the magnesium-based composite material; Apply periodic boundary conditions and load boundary conditions to the RVE mesh model of the magnesium-based composite material, and assign contact properties and constraint boundary conditions between the carbon nanotube model with an interface and the magnesium-based composite material matrix model to obtain an RVE mesh model of the magnesium-based composite material containing boundary conditions; Assign a finite element analysis mode to the RVE mesh model of the magnesium-based composite material containing boundary conditions to obtain a mesoscopic damage numerical simulation analysis model of the magnesium-based composite material; Input the material properties and boundary conditions of the magnesium-based composite material obtained from experiments into the mesoscopic damage numerical simulation analysis model of the magnesium-based composite material, and calculate the mesoscopic damage numerical simulation results. The numerical simulation results include a stress field, a strain field, and a damage field.

2. The mesoscopic damage numerical simulation method for a magnesium-based composite material according to claim 1, characterized in that, The establishment of the magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes includes: Write a script for generating a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes based on the random adsorption sequence algorithm; Assign the volume ratio of carbon nanotubes to the script; Assign the type of the magnesium-based composite material model as a 3D type RVE to the script; Assign the basic dimensions of the magnesium-based composite material model to the script. The basic dimensions include the side length dimension of the magnesium-based composite material matrix model and the length and diameter of the carbon nanotube model; Run the script with parameters assigned to generate a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes.

3. The mesoscopic damage numerical simulation method for a magnesium-based composite material according to claim 2, characterized in that, Use Abaqus finite element simulation software to run the script with parameters assigned to generate a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes.

4. The mesoscopic damage numerical simulation method for a magnesium-based composite material according to claim 1, characterized in that, The thickness of the cohesive element used to replace the interface phase of the magnesium-based composite material is 0.

5. The mesoscopic damage numerical simulation method for a magnesium-based composite material according to claim 1, characterized in that, In the magnesium-based composite material RVE model containing material properties, the material properties include: The material properties, failure modes, damage judgment criteria, and damage evolution of the magnesium-based composite material matrix model; The material properties of the interface phase, the bilinear separation criterion, and the parameters of the cohesive element. The parameters of the cohesive element include the interface phase thickness, the surface force, and the interface separation parameter; The material properties and damage modes of the carbon nanotube model.

6. The mesoscopic damage numerical simulation method of a magnesium-based composite material according to claim 1, characterized in that the mesh generation of the magnesium-based composite material RVE model containing material properties is specifically as follows: The three-node tetrahedral element is used to perform mesh generation on the magnesium-based composite material matrix model and the carbon nanotube model in the magnesium-based composite material RVE model containing material properties respectively; The six-node three-dimensional bonding element is used to perform mesh generation on the interface phase cohesive element in the magnesium-based composite material RVE model containing material properties.

7. The mesoscopic damage numerical simulation method of a magnesium-based composite material according to claim 1, characterized in that the application of periodic boundary conditions to the magnesium-based composite material RVE mesh model includes: Writing a script based on the Python language for applying periodic boundary conditions to the magnesium-based composite material RVE mesh model; The script couples all mesh nodes on the opposite faces, opposite edges, and opposite vertices of the magnesium-based composite material RVE mesh model to complete the application of periodic boundary conditions; The application of load boundary conditions to the magnesium-based composite material RVE mesh model includes: Applying load boundary conditions to the magnesium-based composite material RVE mesh model by using a surface load on the upper surface of the Y axis. The load direction is the positive direction of the Y axis, and the load mode is a smooth analysis step.

8. The mesoscopic damage numerical simulation method of a magnesium-based composite material according to claim 1, characterized in that the finite element analysis mode is the dynamic explicit analysis mode.

9. A mesoscopic damage numerical simulation device for a magnesium-based composite material, characterized in that it includes: The first establishment module is used to establish a magnesium-based composite material model with randomly distributed three-dimensional carbon nanotubes. The magnesium-based composite material model includes a carbon nanotube model and a magnesium-based composite material matrix model; The second establishment module is used to establish cohesive elements on the surface of the carbon nanotube model to replace the interface phase of the magnesium-based composite material, and obtain a carbon nanotube model with an interface; The assembly module is used to perform Boolean operations and assembly on the carbon nanotube model with an interface and the magnesium-based composite material matrix model to obtain a magnesium-based composite material RVE model; The first assignment module is used to assign corresponding material properties to the matrix phase, interface phase, and reinforcement phase in the magnesium-based composite material RVE model respectively, and assign a matrix damage constitutive model and a cohesive constitutive criterion to the magnesium-based composite material RVE model to obtain a magnesium-based composite material RVE model containing material properties; The mesh generation module is used to perform mesh generation on the magnesium-based composite material RVE model containing material properties to obtain a magnesium-based composite material RVE mesh model; A second assignment module, configured to apply periodic boundary conditions and load boundary conditions to the magnesium matrix composite RVE mesh model, and assign contact properties and constraint boundary conditions between the carbon nanotube model with an interface and the magnesium matrix composite matrix model, so as to obtain a magnesium matrix composite RVE mesh model including boundary conditions; A third assignment module, configured to assign a finite element analysis mode to the magnesium matrix composite RVE mesh model including boundary conditions, so as to obtain a mesoscopic damage numerical simulation analysis model of the magnesium matrix composite; A calculation module, configured to input the material properties and boundary conditions of the magnesium matrix composite obtained by experiments into the mesoscopic damage numerical simulation analysis model of the magnesium matrix composite, and calculate a mesoscopic damage numerical simulation result, where the numerical simulation result includes a stress field, a strain field, and a damage field.

10. A device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein, when the processor executes the computer program, the steps of a mesoscopic damage numerical simulation method of a magnesium matrix composite as described in any one of claims 1 to 8 are implemented.

Images