Multi-scale finite element method and system based on sub-model nesting and medium
Through a multi-scale finite element method based on sub-model nesting, the problem that the prior art is difficult to simulate the impact of the microstructure of materials under complex loads is solved, and an accurate analysis of the stress, deformation and damage evolution of large-size structures is achieved.
Patent Information
- Application Number
- CN202511013256.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-07-23
- Publication Date
- 2025-08-22
- Estimated Expiration
- 2045-07-23
AI Technical Summary
Existing macro-microscopic multi-scale simulations are difficult to simulate the impact of material microstructure on stress, deformation or damage evolution under complex loads, especially key areas or areas of concern for large-size structures.
A multi-scale finite element method based on sub-model nesting is adopted. Through step-by-step nesting sub-models, combined with Voronoi method and scanning electron microscopy and other technologies, the grains are divided and the grain orientation characteristics are imparted to achieve simulation transfer from macroscopic scale to microscopic scale, and the crystal plastic finite element model is used for analysis.
The simulation and analysis of stress, deformation and damage evolution of microstructures in the areas of concerned areas of materials or structural parts under complex loads is realized, especially the study of deformation and damage evolution behavior of large-size structures.
Smart Images

Figure CN120526907A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the field of materials science and engineering, and in particular to a multi-scale finite element method, system and medium based on sub-model nesting. Background Art
[0002] From a macroscopic perspective, high-stress areas or geometrically discontinuous areas of a material or structure are often weak points for damage generation or failure. From a microscopic perspective, the generation of damage is closely related to the anisotropy and local stress concentration of polycrystalline materials. While existing macro-micro multiscale simulations can reflect the effects of grain size, orientation, and other factors on stress, deformation, or damage evolution at the microscopic scale, they generally struggle to simulate the effects of material microstructure on stress / deformation / damage evolution under complex loads. This is especially true for multiscale simulations that consider the effects of microstructure in critical or relevant areas of large-scale structures under complex loads. Therefore, achieving multiscale simulations of stress, deformation, or damage evolution under complex loads that consider the effects of material microstructure is an urgent problem that needs to be solved. Summary of the Invention
[0003] The present invention provides a multi-scale finite element method, system and medium based on sub-model nesting to solve the problem that existing methods are difficult to simulate the influence of material microstructure on stress, deformation or damage evolution under complex loads.
[0004] To achieve the above objectives, in a first aspect, the present invention relates to a multi-scale finite element method based on sub-model nesting for multi-scale finite element simulation of materials or structural parts, comprising: Step 1: Establish a finite element simulation model of a preset macro-scale material or structural component, determine a pre-extracted first-level sub-model region of the region of interest; set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract a first-level sub-model based on the pre-extracted first-level sub-model region, and record the displacement boundary conditions of the first-level sub-model, wherein the macro-scale is millimeter-level or above; Step 2 uses the recorded displacement boundary conditions of the first-level submodel to determine a pre-extracted second-level submodel region of the region of interest, performs finite element simulation and mesh-independent verification, extracts a second-level submodel based on the pre-extracted second-level submodel region, and records the displacement boundary conditions of the second-level submodel; Step 3 verifies whether the second-level sub-model meets the preset micro-scale simulation. If so, further sub-model extraction is stopped; Step 4 is to divide the grains and assign grain orientation characteristics based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods for the sub-model that can perform microscale simulation, and use the recorded boundary displacement conditions of the sub-model that can perform microscale simulation as the boundary conditions of the sub-model that can perform microscale simulation to obtain a crystal plasticity finite element simulation model, wherein the microscale is sub-millimeter, micron, or nanometer.
[0005] Preferably, step 3 further includes: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th level sub-model to continue the finite element simulation of the i+1-th level sub-model until the extracted i+1-th level sub-model meets the preset micro-scale simulation. Then stop the simulation and obtain the i+1-th level sub-model that can perform the micro-scale simulation, where i is an integer greater than or equal to 2.
[0006] Preferably, the step of establishing a finite element simulation model of a macro-scale material or structural component and determining a pre-extracted first-level sub-model region of a region of interest is as follows: A three-dimensional global model of the material or structural component is established in Abaqus, and meshing is performed using preset reduced integration units. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structural component is performed. The pre-extracted first-level sub-model area of the area of interest is determined through Abaqus pre-analysis.
[0007] Preferably, the region of interest is a preset position or a preset area.
[0008] To achieve the above-mentioned object, in a second aspect, the present invention relates to a multi-scale finite element system based on sub-model nesting, which is used for multi-scale finite element simulation of materials or structural parts, comprising: A first-level sub-model simulation unit is used to establish a finite element simulation model of a material or structural component at a macro scale, determine a pre-extracted first-level sub-model region of a region of interest, set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract a first-level sub-model based on the extracted first-level sub-model region, and record the displacement boundary conditions of the first-level sub-model, wherein the macro scale is millimeter level or above; a second-level sub-model simulation unit, configured to use the recorded displacement boundary conditions of the first-level sub-model to determine a pre-extracted second-level sub-model region of the region of interest, perform finite element simulation of the second-level sub-model, perform finite element simulation and mesh-independent verification, extract a second-level sub-model based on the pre-extracted second-level sub-model region, and record the displacement boundary conditions of the second-level sub-model; a verification unit, configured to verify whether the second-level sub-model satisfies a preset micro-scale simulation, and if so, stop further sub-model extraction; A micro-simulation unit is used to divide grains and assign grain orientation characteristics based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods for a sub-model that can perform micro-scale simulation, and use the recorded boundary displacement conditions of the sub-model that can perform micro-scale simulation as the boundary conditions of the sub-model that can perform micro-scale simulation to obtain a crystal plasticity finite element simulation model, wherein the micro-scale is sub-millimeter, micron, or nanometer.
[0009] Preferably, the verification unit is further used to: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th level sub-model to continue the finite element simulation of the i+1-th level sub-model until the extracted i+1-th level sub-model meets the preset micro-scale simulation. Then stop the simulation and obtain the i+1-th level sub-model that can be used for micro-scale simulation, where i is an integer greater than or equal to 2.
[0010] Preferably, the step of establishing a finite element simulation model of a macro-scale material or structural component and determining a pre-extracted first-level sub-model region of a region of interest is as follows: A three-dimensional global model of the material or structural component is established in Abaqus, and meshing is performed using preset reduced integration units. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structural component is performed. The pre-extracted first-level sub-model area of the area of interest is determined through Abaqus pre-analysis.
[0011] Preferably, the region of interest is a preset position or a preset area.
[0012] To achieve the above objectives, the third aspect of the present invention also relates to a computer-readable storage medium, in which instructions are stored, and when the instructions are run, the above-mentioned multi-scale finite element method based on sub-model nesting is executed.
[0013] The present invention relates to a multi-scale finite element method, system, and medium based on sub-model nesting, which has the following beneficial effects compared to the prior art: This invention is used to simulate and analyze the effects of the microstructure of a material or structural component's region or location of interest on stress, deformation, or damage evolution under complex loads. By nesting sub-models in a hierarchical manner, the complex load action can be transferred from the macroscopic scale of the material or structure to the microscopic region of the region or location of interest, enabling targeted simulation and analysis of the stress, deformation, and damage in the microscopic region. Combined with grain plasticity finite element analysis, this method can simulate and analyze the stress, deformation, and damage evolution of a material or structural component's region or location of interest under complex loads, taking into account the influence of grain size, grain orientation, and other factors. This invention provides a method for studying the deformation and damage evolution behavior of materials and structures, particularly large-scale structures, under complex loads. BRIEF DESCRIPTION OF THE DRAWINGS
[0014] Figure 1 This is a flow chart of a multi-scale finite element method based on sub-model nesting in Example 1.
[0015] Figure 2 It is a geometric model of a reinforced shell structure and a diagram of positions of interest of Example 1 of a multi-scale finite element method based on nested sub-models in Example 1.
[0016] Figure 3 It is a diagram of boundary conditions and applied loads of a reinforced shell structure of Example 1 of a multi-scale finite element method based on nested sub-models in Example 1.
[0017] Figure 4 It is a Mises stress cloud diagram in the global model of Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0018] Figure 5 It is a displacement cloud map in the global model of Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0019] Figure 6 This is a mesh division diagram of a global model and sub-models at different levels of Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0020] Figure 7 This is a displacement comparison diagram of the global model sub-model area and the sub-model simulation in Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0021] Figure 8 This is a grain morphology and mesh division diagram of a 10 μm×10 μm×10 μm sub-model crystal plasticity finite element model of Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0022] Figure 9This is the stress distribution result of the macroscopic finite element simulation of the 10 μm×10 μm×10 μm sub-model area of Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0023] Figure 10 This is a stress distribution result diagram of the crystal plasticity finite element simulation of a 10 μm×10 μm×10 μm sub-model area in Example 1 of a multi-scale finite element method based on sub-model nesting in Example 1.
[0024] Figure 11 It is a structural diagram of a multi-scale finite element system based on sub-model nesting in Example 2 of the present invention. DETAILED DESCRIPTION
[0025] The present invention will be further described in detail below with reference to the accompanying drawings and examples. It will be understood that the specific embodiments described herein are intended only to illustrate the present invention and are not intended to limit the present invention. It should also be noted that, for ease of description, the accompanying drawings only illustrate portions relevant to the present invention, not all structures.
[0026] Example 1
[0027] A multiscale finite element method based on nested submodels, see Figure 1 As shown, the present invention provides a multi-scale finite element method based on sub-model nesting, which is used for multi-scale finite element simulation of materials or structural parts, including the following steps: S101 to S104.
[0028] S101: Establish a finite element simulation model of a macro-scale material or structural part, determine the pre-extracted first-level sub-model area of the area of interest; set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract the first-level sub-model based on the pre-extracted first-level sub-model area, and record the displacement boundary conditions of the first-level sub-model, where the macro-scale is millimeters or above.
[0029] Macroscale: Focuses on the overall behavior and properties of materials or structures, including the strength, stress-strain relationship, etc. Here, the index is above the millimeter level, such as 10 -3 Common observation methods include uniaxial / multiaxial experiments, on-site stress testing, and full-field displacement observations.
[0030] Microscopic scale: refers to the study of very small structures or phenomena within materials or structures, such as crystal structure, crystal deformation, etc. The research objects include crystal deformation, micro cracks, etc., and the order of magnitude is at the nanometer, micrometer or sub-millimeter level, for example, at 10 -8 meters to 10 -3Commonly used observation methods include scanning electron microscopy (SEM), atomic force microscopy (AFM), nanoindentation experiments, etc. In this embodiment, a finite element simulation model of a macro-scale material or structural component may be established based on a target element, for example, the target element is stress.
[0031] In this embodiment, the region of interest may be a preset position or a preset area, and the preset position is the position of interest.
[0032] In this embodiment, a finite element simulation model of a macro-scale material or structural component is established to determine the pre-extracted first-level sub-model region of the region of interest, specifically: A 3D global model of the material or structure is established in Abaqus, and meshing is performed using preset reduced integration elements. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structure is performed. Pre-analysis in Abaqus determines the pre-extracted first-level sub-model area of the area of interest.
[0033] S102: Using the recorded displacement boundary conditions of the first-level sub-model, determine the pre-extracted second-level sub-model area of the area of interest, perform finite element simulation and mesh-independent verification, extract the second-level sub-model based on the pre-extracted second-level sub-model area, and record the displacement boundary conditions of the second-level sub-model.
[0034] Specifically, the displacement boundary conditions of the recorded first-level sub-model are used to determine the pre-extracted second-level sub-model area of the area of interest, finite element simulation and mesh-independent verification of the second-level sub-model are carried out, the second-level sub-model is extracted based on the pre-extracted second-level sub-model area, and the displacement boundary conditions of the second-level sub-model are recorded.
[0035] S103: Verify whether the second-level sub-model meets the preset micro-scale simulation. If so, stop further sub-model extraction.
[0036] In this embodiment, S103 further includes: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th sub-model to continue the finite element simulation of the i+1-th sub-model until the extracted i+1-th sub-model meets the micro-scale simulation. Then stop the simulation and obtain the i+1-th sub-model that can be used for micro-scale simulation, where i is an integer greater than or equal to 2, and the micro-scale is sub-millimeter, micron or nanometer.
[0037] S104: For the sub-model capable of micro-scale simulation, based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods, the grains are divided and given grain orientation characteristics, and the boundary displacement conditions extracted from the upper-level sub-model of the sub-model capable of micro-scale simulation are recorded as the boundary conditions of the sub-model capable of micro-scale simulation, to obtain a crystal plasticity finite element simulation model.
[0038] The resulting crystal plasticity finite element simulation model can be used to simulate and analyze the stress, deformation, and damage evolution of materials or structural parts under complex loads in areas or locations of interest, taking into account the influence of grain size, grain orientation characteristics, etc.
[0039] In order to better illustrate the solution of the present invention, Example 1 is given below. Figure 2-10 As shown, the following steps are included: Example 1: The following is a focus on the titanium alloy reinforced shell structure under high hydrostatic pressure ( Figure 2 ) is taken as an example to further illustrate.
[0040] A 3D global model was established in Abaqus, and the sub-model area was determined by preliminary analysis. The mesh was divided using a 20-node hexahedral quadratic reduced integration element (C3D20R). The boundary conditions and applied loads were as follows: Figure 3 Based on the linear elastic constitutive relationship, finite element simulation is carried out in combination with the material elastic modulus, density and Poisson's ratio. Figure 4 The Mises stress distribution cloud diagram of the global model of the reinforced shell structure is given. Figure 5 The displacement distribution cloud diagram of the global model of the reinforced shell structure is given. Combined with the finite element simulation results, the Figure 2 The area of interest is located in the 10 mm × 10 mm × 7 mm sub-model ( Figure 6 b) and verify the independence of the sub-model mesh. Record the displacement boundary conditions of the sub-model. Figure 6 In the figure, (a) global model, (b) 10 mm×10 mm×7 mm sub-model, (c) 1 mm×1 mm×1 mm sub-model, (d) 100 μm×100μm×100 μm sub-model, and (e) 10 μm×10 μm×10 μm sub-model.
[0041] Similarly, extract Figure 2 The next level sub-model of the area where the position of interest is located is used to carry out finite element simulation of the sub-model until the extracted sub-model meets the micro-scale simulation requirements. Figure 4 and Figure 5 The size and mesh information of sub-models at different levels are given. Figure 7 Given Figure 6 Part (e) shows the displacement comparison between the sub-model region of the global model and the sub-model simulation when the sub-model is extracted from the global model. It can be seen that the displacement contours of the sub-model region of the global model and the sub-model simulation are well consistent, and the error in node displacements extracted under the same path is less than 1.1%.
[0042] Development Figure 6 Finite element simulation of crystal plasticity in the 10 μm×10 μm×10 μm sub-model region shown in part (e). The Voronoi method is used here to construct the polycrystalline microstructure morphology of the sub-model region, considering only the α phase, and randomly generating the orientation parameters of each grain based on Matlab. The obtained grain morphology and mesh division are shown in the figure below. Figure 8 In the actual simulation process, information such as grain morphology and crystallographic orientation should be obtained through scanning electron microscopy, electron backscatter diffraction and other technologies, and then a crystal plasticity finite element model reflecting the actual material microstructure characteristics should be constructed. Figure 6 The Mises stresses given by the macroscopic finite element simulation and crystal plasticity finite element simulation of the sub-model region shown in part (e) are as follows: Figure 9 and Figure 10 As shown. It can be seen that due to the influence of grain morphology, crystallographic orientation, etc., different from Figure 9 The macro finite element method was used to obtain Figure 6 The Mises stress in the sub-model region shown in part (e) is Figure 10 The Mises stress in the given sub-model area changes dramatically, which realizes the influence and analysis of the material microstructure on the stress of the structural parts under local complex loads (such as the loading of the shell structure under hydrostatic pressure, multi-axial load, etc.).
[0043] Example 2 A multi-scale finite element system based on nested sub-models, such as Figure 11 As shown, the multi-scale finite element simulation for materials or structural parts is implemented by electronic equipment hardware with a central processing unit, which can be implemented by a personal computer, smart terminal, local area network, server, etc. In this embodiment, please refer to Figure 11 , including a first-level sub-model simulation unit 61, a second-level sub-model simulation unit 62, a verification unit 63 and a microscopic simulation unit 64.
[0044] A first-level sub-model simulation unit 61 is configured to establish a finite element simulation model of a material or structural component at a macroscopic scale, determine a pre-extracted first-level sub-model region of a region of interest, set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract a first-level sub-model based on the pre-extracted first-level sub-model region, and record displacement boundary conditions of the first-level sub-model, wherein the macroscopic scale is millimeter-level or above; A second-level submodel simulation unit 62 is configured to use the recorded displacement boundary conditions of the first-level submodel to determine a pre-extracted second-level submodel region of the region of interest, perform finite element simulation and mesh-independent verification, extract a second-level submodel based on the pre-extracted second-level submodel region, and record the displacement boundary conditions of the second-level submodel; A verification unit 63 is used to verify whether the second-level sub-model meets the preset micro-scale simulation, and if so, stop further sub-model extraction; The micro-simulation unit 64 is used to divide the grains and assign grain orientation characteristics based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods for the sub-model that can be simulated at the microscale, and use the recorded boundary displacement conditions of the sub-model that can be simulated at the microscale as the boundary conditions of the sub-model that can be simulated at the microscale to obtain a crystal plasticity finite element simulation model, wherein the microscale is sub-millimeter, micrometer, or nanometer.
[0045] In this embodiment, the verification unit 63 is further configured to: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th sub-model and continue to cyclically carry out the finite element simulation of the i+1-th sub-model until the extracted i+1-th sub-model meets the preset micro-scale simulation. Then stop the simulation and obtain the i+1-th sub-model that can be used for micro-scale simulation, where i is an integer greater than or equal to 2.
[0046] In this embodiment, a finite element simulation model of a macro-scale material or structural component is established to determine the pre-extracted first-level sub-model region of the region of interest, specifically: A 3D global model of the material or structure is established in Abaqus, and meshing is performed using preset reduced integration elements. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structure is performed. Pre-analysis in Abaqus determines the pre-extracted first-level sub-model area of the area of interest.
[0047] In this embodiment, the region of interest is a preset position or a preset region, a preset position and a position of interest.
[0048] The implementation process and effect of the multi-scale finite element system based on sub-model nesting in this embodiment are the same as those of the multi-scale finite element method based on sub-model nesting described in Example 1, and will not be repeated here.
[0049] Example 3 The present invention relates to a computer-readable storage medium, which stores instructions. When the instructions are executed, a multi-scale finite element method based on sub-model nesting is implemented. The implementation process method and effect thereof at runtime are the same as those of the multi-scale finite element method based on sub-model nesting described in Example 1, and will not be repeated here.
[0050] It should be noted that, in this document, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or apparatus comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or apparatus comprising the element.
[0051] The above are only preferred embodiments of the present invention and are not intended to limit the patent scope of the present invention. Any equivalent structure or equivalent process transformation made using the contents of the present invention description and drawings, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present invention.
Claims
1. A multi-scale finite element method based on sub-model nesting, characterized in that: Used for multi-scale finite element simulation of materials or structures, including: Step 1: Establish a finite element simulation model of a macro-scale material or structural component, determine a pre-extracted first-level sub-model region of the region of interest; set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract a first-level sub-model based on the pre-extracted first-level sub-model region, and record the displacement boundary conditions of the first-level sub-model, wherein the macro-scale is millimeter-level or above; Step 2 uses the recorded displacement boundary conditions of the first-level submodel to determine a pre-extracted second-level submodel region of the region of interest, performs finite element simulation and mesh-independent verification, extracts a second-level submodel based on the pre-extracted second-level submodel region, and records the displacement boundary conditions of the second-level submodel; Step 3 verifies whether the second-level sub-model meets the preset micro-scale simulation. If so, further sub-model extraction is stopped; Step 4 is to divide the grains and assign grain orientation characteristics based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods for the sub-model that can perform preset microscale simulation, and use the recorded boundary displacement conditions of the sub-model that can perform microscale simulation as the boundary conditions of the sub-model that can perform microscale simulation to obtain a crystal plasticity finite element simulation model, wherein the microscale is sub-millimeter, micron, or nanometer.
2. The multi-scale finite element method based on sub-model nesting according to claim 1, characterized in that: Step 3 also includes: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th level sub-model to continue the finite element simulation of the i+1-th level sub-model until the extracted i+1-th level sub-model meets the preset micro-scale simulation. Then stop the simulation and obtain the i+1-th level sub-model that can perform the micro-scale simulation, where i is an integer greater than or equal to 2.
3. The multi-scale finite element method based on sub-model nesting according to claim 1, characterized in that: Establish a finite element simulation model of the macro-scale material or structural component and determine the pre-extraction first-level sub-model area of the area of interest, specifically: A three-dimensional global model of the material or structural component is established in Abaqus, and meshing is performed using preset reduced integration units. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structural component is performed. The pre-extracted first-level sub-model area of the area of interest is determined through Abaqus pre-analysis.
4. The multi-scale finite element method based on sub-model nesting according to claim 1, characterized in that: The region of interest is a preset position or a preset area.
5. A multi-scale finite element system based on sub-model nesting, characterized in that: Used for multi-scale finite element simulation of materials or structures, including: A first-level sub-model simulation unit is configured to establish a finite element simulation model of a material or structural component of a preset macroscopic scale, determine a pre-extracted first-level sub-model region of a region of interest, set boundary conditions and load conditions, perform finite element simulation and mesh-independent verification, extract a first-level sub-model based on the pre-extracted first-level sub-model region, and record the displacement boundary conditions of the first-level sub-model, wherein the macroscopic scale is on the order of millimeters or above; a second-level sub-model simulation unit, configured to use the recorded displacement boundary conditions of the first-level sub-model to determine a pre-extracted second-level sub-model region of the region of interest, perform finite element simulation and mesh-independent verification, extract a second-level sub-model based on the pre-extracted second-level sub-model region, and record the displacement boundary conditions of the second-level sub-model; a verification unit, configured to verify whether the second-level sub-model satisfies a preset micro-scale simulation, and if so, stop further sub-model extraction; A micro-simulation unit is used to divide grains and assign grain orientation characteristics based on one or more of the Voronoi method, scanning electron microscopy, and electron backscatter diffraction methods for a sub-model that can be simulated at a microscale, and use the recorded boundary displacement conditions of the sub-model that can be simulated at a microscale as the boundary conditions of the sub-model that can be simulated at a microscale to obtain a crystal plasticity finite element simulation model, wherein the microscale is sub-millimeter, micrometer, or nanometer.
6. The multi-scale finite element system based on sub-model nesting according to claim 5, characterized in that: The verification unit is further configured to: Verify whether the second-level sub-model meets the preset micro-scale simulation. If not, use the recorded displacement boundary conditions of the i-th level sub-model to continue the finite element simulation of the i+1-th level sub-model until the extracted i+1-th level sub-model meets the preset micro-scale simulation. Then stop the simulation and obtain the i+1-th level sub-model that can perform the micro-scale simulation, where i is an integer greater than or equal to 2.
7. The multi-scale finite element system based on sub-model nesting according to claim 5, characterized in that: The method of establishing a finite element simulation model of a macro-scale material or structural component and determining a pre-extracted first-level sub-model region of a region of interest is as follows: A three-dimensional global model of the material or structural component is established in Abaqus, and meshing is performed using preset reduced integration units. Based on the linear elastic constitutive relationship, combined with the material elastic modulus, density, and Poisson's ratio, finite element simulation of the material or structural component is performed. The pre-extracted first-level sub-model area of the area of interest is determined through Abaqus pre-analysis.
8. The multi-scale finite element system based on sub-model nesting according to claim 5, characterized in that: The region of interest is a preset position or a preset area.
9. A computer-readable storage medium, characterized in that: The storage medium stores instructions, which, when executed, execute a multi-scale finite element method based on sub-model nesting according to any one of claims 1 to 4.
Citation Information
Patent Citations
Fatigue damage and service life evaluation method based on Abaqus platform
CN110222439A
Crystal plasticity finite element simulation method of grain boundary reinforced metal matrix composite material
CN116776660A
Metal material compression creep constitutive model construction method
CN120260762A
