CN112100702B - Additive material small crack propagation numerical simulation method considering microstructure - Google Patents
Additive material small crack propagation numerical simulation method considering microstructure Download PDFInfo
- Publication number
- CN112100702B CN112100702B CN202010940897.1A CN202010940897A CN112100702B CN 112100702 B CN112100702 B CN 112100702B CN 202010940897 A CN202010940897 A CN 202010940897A CN 112100702 B CN112100702 B CN 112100702B
- Authority
- CN
- China
- Prior art keywords
- crystal
- crack propagation
- model
- microstructure
- crack
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y50/00—Data acquisition or data processing for additive manufacturing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/10—Additive manufacturing, e.g. 3D printing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- Materials Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to an additive material small crack propagation numerical simulation method considering a microstructure, which comprises the following steps of: (1) establishing a flat piece model according to the size of a 'bone rod type' crack propagation test piece of the additive manufacturing material, partitioning the model according to whether crystal grains are introduced or not, and respectively endowing isotropic and anisotropic material properties; (2) carrying out grid division on the flat piece model, and carrying out grid local encryption on a crack expansion area in the middle of the model; (3) exporting a model geometry and grid attribute file, and combining mathematical characterization of MATLAB software on the crystal morphology and the crystal orientation of the additive manufacturing material to complete microstructure modeling; (4) the crack is prefabricated for the model, the expansion finite element method is adopted to simulate the small crack expansion behavior of the material, and the parameter value in the FATIGUE fracture criterion is determined according to the actual condition of the material; (5) and (3) simulating the small crack propagation behavior of the material based on finite element software ABAQUS to obtain a crack propagation path and a crack propagation rate.
Description
Technical Field
The invention discloses a numerical simulation method for small crack propagation of an additive material by considering a microstructure, which is a numerical simulation method for small crack propagation of a special microstructure of an additive manufacturing material, including crystal morphology and crystal orientation influence, and belongs to the technical field of aerospace engines.
Background
The additive manufacturing is different from the traditional manufacturing process, the materials are formed by stacking from bottom to top layer by layer, and the additive manufacturing has the characteristics of quick processing, low production cost, short manufacturing period, high material utilization rate and the like, and has remarkable advantages in the aspect of processing complex structures. The structure can be further lightened and the cost can be reduced when the structure is applied to an aircraft engine.
The rapid melting and cooling solidification in the additive manufacturing process enable the formed material to have metallurgical defects such as unfused, air holes and inclusions, the crack initiation can be accelerated, and the crack propagation accounts for more than 90% of the fatigue life of the structure in three stages of fatigue failure. Research also shows that the structural fatigue crack propagation life is mostly consumed in a small crack stage which is greatly influenced by the microstructure, and the propagation mechanism of the structural fatigue crack propagation life is obviously different from that of the conventional long crack. Therefore, in the fatigue failure process of the additive manufacturing structural component, the influence of the microstructure must be considered, and the small crack propagation stage of the additive manufacturing structural component is researched intensively.
However, the additive manufacturing material is affected by high temperature gradient during processing, and the microstructure thereof is greatly different from that of the traditional manufacturing and forming. For example, coarse columnar primary β grains growing along the direction of the deposit exist in the microstructure of the selected area laser melting TC4, α + β lamellar structure is inside the grains, and the material will exhibit significant anisotropy. Therefore, it becomes important how to consider the effect of a particular microstructure on the propagation of small cracks.
At present, the research on the small crack propagation of the additive manufacturing material is mainly based on tests, and the tests have the problems of high cost, difficulty in operation, low process repeatability and the like, so that the research on the numerical simulation method for the small crack propagation of the additive manufacturing material considering the microstructure is innovative.
Disclosure of Invention
The technical scheme of the invention is as follows: the method can consider the influence of the special microstructure of the additive manufacturing material, including crystal morphology and crystal orientation, realize effective simulation of the small crack propagation behavior of the additive manufacturing material, and lay a foundation for accurate prediction of the subsequent fatigue life of the additive manufacturing material.
The technical scheme of the invention is as follows: aiming at a special microstructure of an additive manufacturing material, an MATLAB script program is compiled through an INP data file of commercial finite element software ABAQUS to realize the characterization of the crystal morphology and the crystal orientation of the material, and a propagation finite element method is adopted to simulate the propagation behavior of a small crack to obtain a crack propagation path and a crack propagation rate.
The method comprises the following implementation steps:
the method comprises the steps of firstly, establishing a flat plate model by only keeping the minimum section part of the stepped shape of a test piece and combining the Shengweinan principle to improve the calculation efficiency on the basis of the size of a 'bone-rod-shaped' crack propagation test piece made of an additive manufacturing material. In order to further reduce the number of grid cells and reduce the calculation amount, a model area is divided according to whether crystal grains are introduced, only a crack propagation area in the middle part of the model is defined as an area considering a microstructure, and anisotropic material properties are given to the crack propagation area; the remaining regions do not take into account microstructural effects, i.e. no grains need to be introduced, giving isotropic material properties.
And secondly, carrying out grid division on the established flat piece model, wherein the linear reduction integration unit is more accurate in displacement solving result, so that the unit type is an 8-node hexahedron linear reduction integration unit (C3D 8R). The model has high requirement on the calculation precision of the middle imported grain region, so that the grid units in the region need to be encrypted, the number of seed points on each side of the grid units is increased, a local precise grid is established, and the solving precision of the numerical simulation method on the small crack propagation behavior of the material is ensured; in addition, a cell transition region is arranged around the region, so that the cell size is transited from dense to sparse.
And thirdly, deriving an INP data file in ABAQUS finite element software based on the flat piece model after meshing, extracting three-dimensional coordinates of mesh unit nodes and serial numbers corresponding to 8 nodes of each unit, generating a corresponding number of crystal core points in an area needing to be led in the crystal grains according to information such as real crystal grain size and crystal morphology of the additive manufacturing material, and partitioning each mesh unit by using a Voronoi algorithm to finish crystal morphology characterization. And (2) writing an MATLAB script program according to the corresponding position coordinates and the strength of each position in the additive manufacturing material EBSD pole figure, wherein the program comprises the following steps:
the Euler angle method is selected to represent a sufficient number of randomly generated crystal orientations passing through the Euler angleφ、The rotating XYZ sample coordinate system can be expressed as g in hkl crystal coordinate system, which expression is as follows:
wherein three crystal axes of the hkl crystal coordinate system are firstly rotated around the Z axis from the XYZ sample coordinate systemAngle, the rotated coordinate system being rotated again about the X-axisAngle, finally rotation about Z axisThe angle is obtained. Then, the hkl crystal coordinate system can be expressed as g in the XYZ sample coordinate system-1I.e. the crystal orientation can be expressed as g-1。
For the representation of the position coordinates of a certain crystal direction in each projection plane in the XYZ sample coordinate system, taking the XOY plane as the projection plane as an example, the expression of the position coordinates (x, y) is as follows:
wherein, alpha is the included angle between the crystal direction and the Z axis, beta is the included angle between the crystal direction and the X axis, and gamma is the included angle between the crystal direction and the Y axis.
And screening out crystal orientations which meet the positions and the texture strength of the EBSD pole figures of the material according to the number of crystal grains required by the middle introduced crystal grain area, and endowing the crystal orientations to the crystal grains to finish the crystal orientation characterization. And importing the grid unit partition information obtained by running an MATLAB script program and the screened crystal orientation information into an INP file so as to update the established basic model.
Fourthly, prefabricating cracks for a flat piece model with introduced crystal grains, simulating the small crack propagation behavior of the material by adopting a finite element propagation method, wherein the analysis step adopts a direct cycle analysis step, the fracture criterion adopts a FATIGUE criterion, and the conditions for starting the crack propagation of an enrichment unit in a certain layer are as follows:
ΔG=Gmax-Gmin
wherein G ismax、GminThe strain energy release rates corresponding to the maximum load and the minimum load during cyclic loading, c1、c2The number of cycles N required for the crack to start propagating is determined as a material parameter. Taking the value c, irrespective of crack initiation, due to the presence of pre-cracks1=0、c 20. The propagation rate of the layer of enrichment units after the crack begins to propagate is as follows:
Gthresh<Gmax<Gpl
wherein G isthresh、GplRespectively satisfying the upper and lower limits of the strain energy release rate of Paris area, c3、c4Determining the propagation rate of the crack after the crack begins to propagate as a material parameterc3、c4The value is calculated by C, m in the Paris formula of the material, and the expression is as follows:
wherein E is the elastic modulus of the material, and mu is the Poisson's ratio.
And fifthly, setting load and boundary conditions according to actual conditions, simulating the small crack propagation behavior of the additive manufacturing material by adopting an expansion finite element method of ABAQUS finite element software based on the established model, and acquiring a small crack propagation path and a crack propagation rate in a visualization module.
Compared with the prior art, the invention has the advantages that:
according to the method, the influence of the special microstructure of the additive manufacturing material is introduced into the numerical simulation of the small crack propagation of the material, and the microstructure modeling is completed through the characterization of the crystal morphology and the crystal orientation of the material, so that the simulation of the small crack propagation behavior of the additive manufacturing material is more accurate, and the related research deficiency is made up.
Drawings
FIG. 1 is a flow chart of an embodiment of a method for simulating a small crack propagation value of an additive material considering a microstructure according to the present invention;
FIG. 2 is a schematic diagram of a finite element model of additive material small crack propagation according to an embodiment;
FIG. 3(a) is a simulation result of a small crack propagation path of an additive material according to an embodiment;
fig. 3(b) is a simulation result of a small crack propagation rate of an additive material according to an embodiment.
Detailed Description
The technical scheme of the additive material small crack propagation numerical simulation method considering the microstructure is further explained below with reference to the accompanying drawings. The material to be investigated in this example was SLM (selective laser melting) TC4 alloy.
As shown in fig. 1, the specific implementation process of the present invention is as follows:
in the first step, a flat plate piece basic model with the size of 4.5mm by 1.5mm by 0.05mm is established for improving the calculation efficiency according to the size of an SLM TC4 material 'bone bar type' crack propagation test piece. The width of the basic model only keeps the minimum section part of the test piece in the step shape, namely 1.5 mm; the Saint Vietnam principle is combined, and the length of the model is only kept 4.5 mm; the model is changed from a 3-dimensional model to a 2.5-dimensional model, namely, the thickness of the model is reduced, and the dimension in the direction is only kept 0.05 mm. In order to further reduce the number of grid cells and the calculation amount, the model area is divided according to whether the crystal grains are introduced, only the area 1.0mm x 0.2mm in the middle of the model is defined as the area considering the microstructure, and anisotropic material properties are given to the area, as shown in table 1; the remaining regions are not considered for microstructural effects, i.e. no grains need to be introduced, giving isotropic material properties, as shown in table 2.
TABLE 1 SLM TC4 Material Anisotropic elastic constant
TABLE 2 SLM TC4 Material Isotropic Material Properties
And secondly, carrying out grid division on the established flat piece model, wherein the linear reduction integration unit is more accurate in displacement solving result, so that the unit type is an 8-node hexahedron linear reduction integration unit (C3D 8R). The model has high requirement on the calculation accuracy of the middle imported grain region, so that the grid cells in the region need to be encrypted, the number of seed points on each side of the grid cells is increased, a local precise grid is established, the cell size is finally determined to be 0.005mm through simulation verification, and the calculation amount is reduced as much as possible on the premise of ensuring the solving accuracy of the numerical simulation method on the small crack propagation of the material; in addition, a cell transition region is arranged around the region, so that the cell size is transited from dense to sparse, and 12524 grid cells are formed by combining models.
And thirdly, deriving an INP data file of ABAQUS finite element software based on the flat piece model after the meshing is finished, and extracting the three-dimensional coordinates of the mesh unit nodes and the serial numbers corresponding to 8 nodes of each unit. The microstructure of the SLM TC4 material has columnar original beta grains growing along the direction of a deposited layer, the grains are in an alpha cluster-like structure, and the same alpha cluster has the same or similar crystal orientation, so that the alpha cluster is used as a basic unit for characterization of crystal morphology. Based on the above, a corresponding number of crystal nucleus points are generated in the region of 1.0mm x 0.2mm in the middle of the crystal grains to be introduced, 8000 grid units are partitioned by using a Voronoi algorithm, and the crystal morphology characterization is completed. Writing an MATLAB script program by taking the corresponding position coordinates and the intensity of each position in the SLM TC4 material EBSD pole figure as the basis, wherein the calculation steps of the program comprise:
the Euler angle method is selected to represent a sufficient number of randomly generated crystal orientations passing through the Euler angleφ、The rotating XYZ sample coordinate system can be expressed as g in hkl crystal coordinate system, with the expression:
wherein three crystal axes of the hkl crystal coordinate system are firstly rotated around the Z axis from the XYZ sample coordinate systemAngle, the rotated coordinate system being rotated again about the X-axisAngle, finally rotation about Z axisThe angle is obtained. Then, the hkl crystal coordinate system can be expressed as g in the XYZ sample coordinate system-1I.e. the crystal orientation can be expressed as g-1。
Regarding the representation of the position coordinates of a certain crystal orientation in each projection plane in the XYZ sample coordinate system, taking the XOY plane as the projection plane as an example, the expression of the position coordinates (x, y) is as follows:
wherein, alpha is the included angle between the crystal direction and the Z axis, beta is the included angle between the crystal direction and the X axis, and gamma is the included angle between the crystal direction and the Y axis.
And screening out crystal orientations which meet the positions and the texture strength of the EBSD pole figure of the SLM TC4 material according to the number of crystal grains required by the area of the middle introduced crystal grains, and endowing the crystal orientations to the crystal grains to finish the crystal orientation characterization. And importing the grid unit partition information obtained by running an MATLAB script program and the screened crystal orientation information into an INP data file so as to update the established basic model.
The model was created as shown in figure 2. The middle part displays the whole situation of the basic model of the flat piece; the right side of the figure shows the area after the crystal grains are led in the middle in an enlarged manner, and unit areas with different colors and gray levels represent different crystal grains; the left side of the figure shows the crystal orientation of the model after screening by the MATLAB script program.
And fourthly, prefabricating a crack with the length of 0.1mm for a flat piece model with the introduced crystal grains, simulating the small crack propagation behavior of the material by adopting a finite element expansion method, wherein the analysis step adopts a direct cyclic analysis step, the increment of the analysis step is 0.01, and the initial value of the Fourier series is 24. The fracture criterion is FATIGUE criterion, and the condition that the crack of a certain layer of enrichment unit begins to expand is as follows:
ΔG=Gmax-Gmin
wherein G ismax、GminThe strain energy release rates corresponding to the maximum load and the minimum load during cyclic loading, c1、c2The number of cycles N required for the crack to start propagating is determined as a material parameter. Taking the value c, irrespective of crack initiation, due to the presence of pre-cracks1=0、c 20. The layer is rich inThe propagation rate after the initiation of the propagation of the unit crack was as follows:
Gthresh<Gmax<Gpl
wherein G isthresh、GplRespectively satisfying the upper and lower limits of the strain energy release rate of Paris area, c3、c4Determining the propagation rate of the crack after the crack begins to propagate as a material parameterc3、c4The value is calculated by C, m in the Paris formula of the material, and the expression is as follows:
wherein E is the elastic modulus of the material, and mu is the Poisson's ratio. The SLM TC4 material Paris formula has C corresponding to a value of 3.311 × 10-8M corresponds to a value of 2.99, E115 GPa, μ 0.324, and E' 128.488GPa in a plane strain state. From this, it can find3=1.495,c4=4.70663×10-5。
Fifthly, setting boundaries and load conditions according to actual conditions, wherein the boundary conditions of the model are that the left side is hinged, the right side is circularly loaded, and the maximum stress sigma of the circular load ismax380MPa, stress ratio R of 0.06 and loading frequency of 1 Hz. Based on the established model, the ABAQUS finite element software expansion finite element method is adopted to simulate the small crack expansion behavior of the SLM TC4 material, and a small crack expansion path and a crack expansion rate are obtained in a visualization module. The simulation results are shown in fig. 3.
The above examples are provided only for the purpose of describing the present invention, and are not intended to limit the scope of the present invention. The scope of the invention is defined by the appended claims. Various equivalent substitutions and modifications can be made without departing from the spirit and principles of the invention, and are intended to be within the scope of the invention.
Claims (3)
1. A method for simulating a small crack propagation value of an additive material by considering a microstructure is characterized by comprising the following steps:
step (1): establishing a flat piece model according to the size of a 'bone rod type' crack propagation test piece of the additive manufacturing material, partitioning the model according to whether crystal grains are introduced or not, and respectively endowing isotropic and anisotropic material properties; the bone-bar-type crack propagation test piece is a test piece used in an in-situ fatigue test of an additive manufacturing material, and a flat piece model is established based on the Saint-Venn principle on the basis of the size of the test piece, so that the calculation efficiency can be improved on the premise of ensuring the accuracy of a result compared with a complete bone-bar-type crack propagation model; the step of partitioning the model according to whether the grains are introduced or not refers to that the microstructure modeling is only carried out in the crack expansion area in the middle of the model, the influence of the microstructure is not considered in the rest areas, and compared with the mode that the grains are introduced into the whole model, the subsequent calculation amount can be reduced; respectively endowing material properties means endowing the region with the crystal grains introduced in the middle with anisotropic material properties, and endowing the rest regions with isotropic material properties;
step (2): carrying out mesh division on the flat piece model established in the step (1), and carrying out local mesh unit encryption on the middle introduced grain region; the grid division refers to directly distributing and dividing the seeds on each side of the flat piece model to establish an integral coarse grid, wherein the unit type is an 8-node hexahedron linear reduction integral unit C3D 8R; the grid cell encryption is to increase the number of seed points on each side of a grain region led into the middle of the grain region, establish a local precise grid, and set a cell transition region around the region to ensure that the cell size is transited from dense to sparse so as to ensure the solving precision of the numerical simulation method on the small crack propagation behavior of the material;
and (3): deriving an INP file in ABAQUS finite element software based on the flat piece model after the gridding division is completed, and characterizing the crystal morphology and the crystal orientation of the additive manufacturing material by using MATLAB software to complete the microstructure modeling; the step of exporting the INP file in the ABAQUS finite element software is to extract the three-dimensional coordinates of the grid unit nodes and the serial numbers corresponding to 8 nodes of each unit; the characterization of the crystal morphology of the additive manufacturing material by using MATLAB software refers to the generation of a corresponding number of crystal nucleus points in an area in which crystal grains need to be introduced according to the actual crystal grain size, crystal morphology and material texture information of the additive manufacturing material, and the division of grid units of the additive manufacturing material by using a Voronoi algorithm, wherein the units in the same area form one crystal grain; the characterization of the crystal orientation of the additive manufacturing material by using MATLAB software refers to extracting position coordinates and intensity of each position in an electron back scattering diffraction EBSD polar diagram, screening out qualified orientation from a sufficient number of randomly generated crystal orientations on the basis of the position coordinates and the intensity, and endowing the qualified orientation to each crystal grain; the step of finishing the microstructure modeling is to run an MATLAB script program, introduce the obtained grid unit partition information and the crystal orientation information into an INP file and update the established basic model according to the grid unit partition information and the crystal orientation information;
and (4): prefabricating cracks for the updated flat piece model in the step (3), simulating the small crack propagation behavior of the material by adopting a finite element propagation method, and determining each relevant parameter value in the FATIGUE fracture criterion according to the actual condition of the material; the updated flat piece model is a model after the microstructure modeling is completed, and comprises the representation of the crystal morphology and the crystal orientation in the middle area; the prefabricated crack refers to an initial crack which needs to be inserted into a flat piece model for a certain length when a FATIGUE fracture criterion is used;
and (5): simulating the small crack propagation behavior of the additive manufacturing material based on ABAQUS finite element software to obtain a crack propagation path and a crack propagation rate; the crack propagation path is acquired by directly displaying the options of PHILSM and STATUS SXFEM in ABAQUS software field variable output in a visualization module; the crack growth rate is obtained by plotting the crack growth rate da/dN obtained by the secant method on the basis of the simulation result against the crack length a.
2. The method for simulating the small crack propagation value of the additive material considering the microstructure according to claim 1, wherein:
the characterization of the material crystal morphology in the step (3) is based on the real microstructure of the material, and an MATLAB script program is further compiled by extracting the grid unit information in the INP data file, and specifically comprises the following steps:
the Euler angle method is selected to represent a sufficient number of randomly generated crystal orientations passing through the Euler angleφ、The rotating XYZ sample coordinate system can be expressed as g in hkl crystal coordinate system, with the expression:
wherein three crystal axes of the hkl crystal coordinate system are firstly rotated around the Z axis from the XYZ sample coordinate systemAngle, the rotated coordinate system being rotated again about the X-axisAngle, finally, rotation about Z axisObtaining an angle; then, the hkl crystal coordinate system is expressed as g in the XYZ sample coordinate system-1I.e. the crystal orientation is expressed as g-1;
Regarding the representation of the position coordinates of a certain crystal orientation in each projection plane in the XYZ sample coordinate system, taking the XOY plane as the projection plane as an example, the expression of the position coordinates (x, y) is as follows:
wherein, alpha is the included angle between the crystal direction and the Z axis, beta is the included angle between the crystal direction and the X axis, and gamma is the included angle between the crystal direction and the Y axis;
and screening out the crystal orientation which meets the positions of the EBSD pole figure and the texture intensity of the material according to the number of the crystal grains required by the middle introduced crystal grain area, and endowing the crystal orientation to each crystal grain to finish the crystal orientation characterization.
3. The method for simulating the small crack propagation numerical value of the additive material considering the microstructure according to claim 1, wherein:
the characterization of the material crystal orientation in the step (3) is based on the material texture of the EBSD pole figure, and the random orientation is screened by writing an MATLAB script program, so as to obtain the crystal orientation according with the position coordinates and the strength at each position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010940897.1A CN112100702B (en) | 2020-09-09 | 2020-09-09 | Additive material small crack propagation numerical simulation method considering microstructure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010940897.1A CN112100702B (en) | 2020-09-09 | 2020-09-09 | Additive material small crack propagation numerical simulation method considering microstructure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100702A CN112100702A (en) | 2020-12-18 |
CN112100702B true CN112100702B (en) | 2022-06-24 |
Family
ID=73751241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010940897.1A Active CN112100702B (en) | 2020-09-09 | 2020-09-09 | Additive material small crack propagation numerical simulation method considering microstructure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100702B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378432B (en) * | 2021-05-28 | 2022-10-21 | 南京理工大学 | Numerical simulation method for crack propagation on RPV pipe pitting pit based on propagation finite element |
CN114186461A (en) * | 2021-12-08 | 2022-03-15 | 北京理工大学 | Multi-axial fatigue life prediction method |
CN114065553B (en) * | 2021-12-31 | 2024-07-26 | 华侨大学 | Method for establishing anisotropic monocrystalline material molecular dynamics model |
CN114638131A (en) * | 2022-01-28 | 2022-06-17 | 南京航空航天大学 | Crystal plasticity finite element modeling method based on electron back scattering diffraction analysis and Gmesh software |
CN114626263B (en) * | 2022-02-25 | 2024-06-11 | 北京航空航天大学 | High-temperature alloy material short crack propagation numerical simulation method based on crystal plasticity |
CN117150858B (en) * | 2023-09-06 | 2024-03-26 | 天津大学 | Crack-containing continuous fiber reinforced resin matrix composite three-dimensional finite element modeling method |
CN117854652B (en) * | 2024-01-16 | 2024-07-23 | 华中科技大学 | Modeling method and system for weld microstructure crystal plasticity finite element model |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2858410B1 (en) * | 2003-07-28 | 2005-09-23 | Electricite De France | METHOD FOR DETERMINING CONSTRAINTS, DEFORMATIONS, DAMAGE TO PARTS CONSISTING OF SOLID MATERIAL |
CN108629147B (en) * | 2018-05-02 | 2021-07-09 | 西安电子科技大学 | Polycrystalline geometric modeling method |
CN109725123B (en) * | 2019-02-28 | 2020-05-22 | 北京航空航天大学 | Crack propagation life determination method considering shot peening strengthening surface layer grain refinement |
CN110926942B (en) * | 2019-11-19 | 2022-01-07 | 西安理工大学 | Numerical analysis method for rolling contact fatigue cracks of rails in ABAQUS |
-
2020
- 2020-09-09 CN CN202010940897.1A patent/CN112100702B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN112100702A (en) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112100702B (en) | Additive material small crack propagation numerical simulation method considering microstructure | |
KR102313343B1 (en) | Systems and methods for predicting structure and properties of atomic elements and alloy materials thereof | |
CN105373683B (en) | A kind of prediction technique of 20CrMnTiH steel thermal deformation process Microstructural Evolution rule | |
Romanova et al. | A method of step-by-step packing and its application in generating 3D microstructures of polycrystalline and composite materials | |
CN106649964B (en) | GA-ELM algorithm-based aluminum alloy die casting grain size prediction method | |
CN103257214A (en) | Material property distribution determination for fatigue life calculation using dendrite arm spacing and porosity-based models | |
Liu et al. | A numerical approach to simulate 3D crack propagation in turbine blades | |
Zhao et al. | Fatigue crack propagation within Al-Cu-Mg single crystals based on crystal plasticity and XFEM combined with cohesive zone model | |
Han et al. | Temperature-dependent fatigue response of a Fe44Mn36Co10Cr10 high entropy alloy: A coupled in-situ electron microscopy study and crystal plasticity simulation | |
Ji et al. | Grain scale modelling and parameter calibration methods in crystal plasticity finite element researches: A short review | |
Rollett et al. | Modeling polycrystalline microstructures in 3D | |
Zhu et al. | Plastic deformation mechanism and defect patterning under nanoindentation in medium entropy alloy CoCrNi | |
He et al. | Effects of secondary orientation on fatigue crack initiation in a single crystal superalloy | |
Xu et al. | Temperature-dependent, multi-mechanism crystal plasticity reveals the deformation and failure behaviour of multi-principal element alloys | |
Klippel | Constitutive equations for simulation of metal cutting with meshless methods on GPU | |
CN108256246A (en) | Design method and device of porous material based on Unity3D | |
CN115007820B (en) | Ultrasonic continuous casting method and system for reducing wall thickness sensitivity of spherical graphite tissue | |
CN112163356A (en) | Method for finite element modeling of additive manufacturing microstructure ultrasonic detection | |
Xiaoguang et al. | A multi-scale framework for life reduction assessment of turbine blade caused by microstructural degradation | |
Wang et al. | In-situ dwell-fatigue fracture experiment and CPFE simulation of SLM AlSi10Mg alloy at high temperature | |
Dai | A Study of Solidification Structure Evolution during Investment Casting of Ni-based Superalloy for Aero-Engine Turbine Blades. | |
Li et al. | Orientation effect on twin formation mechanisms of tungsten under uniaxial tension/compression | |
Moore et al. | A nonconformal nonlocal approach to calculating statistical spread in fatigue indicator parameters for polycrystals | |
Meng et al. | The microstructural evolution and mechanical response of laser direct energy deposition Inconel 718 alloy based on simulation and experimental methods | |
CN115906714B (en) | Plate PECVD reaction chamber airflow simulation method and simulation equipment thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |