Fibrous grain three-dimensional crystal modeling and analyzing method and device

Abstract

The invention discloses a fibrous grain three-dimensional crystal modeling and analyzing method, and belongs to the technical field of finite element analysis. The method comprises the following steps: acquiring a metal wire material fiber grain sample, processing to obtain a ctf file, and importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data; based on the grain size distribution data, constructing a three-dimensional fibrous grain model to obtain a geo file; importing the geo file into Gmsh software, and carrying out grid division to generate an inp file; on the basis of the grain orientation information, Euler angles with the same number as the three-dimensional fibrous grain models are screened by adopting an Euler angle screening method to obtain Euler angle data, and the Euler angle data passing verification are processed to obtain grain orientation parameters; the grain orientation parameters and the inp file are input into ABAQUS, and independent material attributes related to each grain are generated; finite element analysis is conducted on the fiber crystal grain sample according to the independent material attribute of each crystal grain, an analysis result is obtained, and the modeling accuracy of the fibrous crystal grains is improved through the method.

Description

Technical Field

[0001] The present application belongs to the field of finite element analysis technology, and in particular relates to a method and device for modeling and analyzing three-dimensional crystals of fibrous grains. Background Art

[0002] With the rapid development of computational materials science, a growing number of researchers are dedicated to quantitatively linking the microstructure and macroscopic properties of metal wire materials. In crystal plasticity simulations of metal wire materials, the fibrous grain structure (aspect ratio ≥ 3:1) serves as the primary carrier of plastic deformation, and its three-dimensional spatial distribution directly determines the material's anisotropic mechanical behavior. Therefore, to accurately predict key performance indicators such as tensile strength and fatigue life of metal wire materials, it is crucial to construct numerical models that incorporate realistic fibrous grain characteristics.

[0003] Currently, there are two main approaches for crystal plasticity modeling of wire materials: one is to construct a grain geometry model based on metallographic image processing. This method reconstructs grain outlines through binarization and connected region identification, but suffers from image distortion and grain boundary morphology distortion. The other is to build a polycrystalline model based on EBSD experimental data. While this method can preserve grain orientation information, it requires data conversion using multiple software packages, including MATLAB, Python, and ABAQUS, which can easily lead to mismatching between the grain spatial coordinates and orientation data. Existing methods often use regular polygonal grain approximations during modeling, failing to reproduce the axial extension characteristics and gradient grain boundary morphology of fibrous grains. This makes it difficult to accurately characterize the quantitative impact of the evolution of grain preferred orientation on the mechanical response of the material during the wire drawing process. Grain aspect ratios and orientation distributions vary significantly among different wire materials. Existing methods for modeling fibrous grain structures suffer from grain morphology distortion, low modeling efficiency and accuracy, and large data matching errors. Summary of the Invention

[0004] The present application aims to solve at least one of the technical problems existing in the prior art. To this end, the present application proposes a method and apparatus for three-dimensional crystal modeling and analysis of fibrous grains, which improves the accuracy of fibrous grain modeling.

[0005] In a first aspect, the present application provides a method for modeling and analyzing three-dimensional crystals of fibrous grains, the method comprising: Obtaining a fiber grain sample of a metal wire material, processing the fiber grain sample to obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; Importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample; Based on the grain size distribution data, a three-dimensional fibrous grain model is constructed to obtain a geo file including grain geometry information; Import the geo file into Gmsh software for grid division, and generate an inp file based on the divided data; Based on the grain orientation information, an Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters; Input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; Loads and boundary conditions are set based on the ABAQUS software, and finite element analysis is performed on the fiber grain sample according to the independent material properties of each grain to obtain analysis results.

[0006] According to one embodiment of the present application, the processing of the fiber grain sample to obtain a ctf file including grain information includes: The grain information of different cross sections of the fiber grain sample is obtained by using the EBSD probe of the scanning electron microscope, wherein the grain information includes grain orientation and size distribution, and EBSD raw data is obtained; The EBSD raw data is converted into a ctf file including grain information.

[0007] According to one embodiment of the present application, the step of importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample includes: Importing the ctf file into the MTEX ​​toolkit in MATLAB, extracting the lognormal distribution parameters of the cross-sectional grain size and the range of the longitudinal grain aspect ratio, and obtaining the grain orientation information and grain size distribution data of the fiber grain sample; The grain orientation information includes the orientation distribution of the longitudinal section, and the grain size distribution data includes the average grain size distribution data of the cross section and the aspect ratio data of the longitudinal section.

[0008] According to one embodiment of the present application, constructing a three-dimensional fibrous grain model based on the grain size distribution data to obtain a geo file including grain geometric information includes: Based on the grain size distribution data, a three-dimensional fibrous grain model was constructed using Neper software, wherein the aspect ratio of the longitudinal section of the three-dimensional fibrous grain model ranged from 3:1 to 15:1, the width-to-height ratio was 1:1, and the cross-sectional size of the three-dimensional fibrous grain model obeyed a lognormal distribution of μ=0.5-2 μm and σ=0.2-0.5; Where σ is the standard deviation of the cross-sectional data, and μ is the mean of the cross-sectional data.

[0009] According to one embodiment of the present application, the Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters, including: The grain ID sequence is generated based on the grain size using Neper software, and the EBSD raw data are arranged in descending order according to the grain ID sequence; Euler angle data are obtained by randomly selecting the same number of Euler angles as the three-dimensional fibrous grain model using a Monte Carlo algorithm. The Euler angle data file obtained by screening is expanded to a file format that conforms to the MTEX ​​interface and then imported into MTEX ​​to draw orientation difference maps and pole figures, which are compared and verified with the arranged EBSD raw data. The verified Euler angle data is processed to obtain the grain orientation parameters. If the verification fails, it is re-screened until it passes the verification.

[0010] According to one embodiment of the present application, the step of importing the geo file into the Gmsh software for meshing, and generating an inp file based on the meshed data, includes: Import the geo file into Gmsh software, set the maximum grid size to one fifth of the average grain size, and the minimum grid size to one twentieth of the average grain size; Based on the maximum grid size and the minimum grid size, the Delaunay algorithm is used to divide the data in the geo file into tetrahedral grids, and the inp file is generated based on the divided data.

[0011] According to one embodiment of the present application, setting loads and boundary conditions based on the ABAQUS software and performing finite element analysis on the fiber grain sample includes: Loading the UMAT subroutine in the ABAQUS software to define the constitutive model; Setting displacement loads and boundary conditions, wherein the displacement loads include unidirectional tensile loads and shear loads, the angle between the displacement load direction and the axial direction of the fibrous grains is 0°, 45°, or 90°, and the boundary conditions are to fix all degrees of freedom of the nodes at one end of the model and apply the displacement load to the nodes at the other end; Based on the displacement load and boundary conditions, a crystal plasticity finite element analysis is performed on the fiber grain sample to obtain an analysis result.

[0012] In a second aspect, the present application provides a device for modeling and analyzing three-dimensional crystals of fibrous grains, the device comprising: an acquisition module, configured to acquire a fiber grain sample of a metal wire material, process the fiber grain sample, and obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; A first processing module is used to import the ctf file into MATLAB to obtain the grain orientation information and grain size distribution data of the fiber grain sample; A second processing module is used to construct a three-dimensional fibrous grain model based on the grain size distribution data to obtain a geo file including grain geometric information; The third processing module is used to import the geo file into the Gmsh software for grid division, and generate an inp file based on the divided data; a fourth processing module, configured to, based on the grain orientation information, use an Euler angle screening method to screen Euler angles having the same number as the three-dimensional fibrous grain model to obtain Euler angle data, verify the Euler angle data, and obtain grain orientation parameters by processing the verified Euler angle data; a fifth processing module, configured to input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; The analysis module is used to set loads and boundary conditions based on the ABAQUS software, perform finite element analysis on the fiber grain sample according to the independent material properties of each grain, and obtain analysis results.

[0013] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the method for three-dimensional crystal modeling and analysis of fibrous grains as described in the first aspect above is implemented.

[0014] In a fourth aspect, the present application provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the three-dimensional crystal modeling and analysis method for fibrous grains as described in the first aspect above.

[0015] In a fifth aspect, the present application provides a chip comprising a processor and a communication interface, wherein the communication interface is coupled to the processor, and the processor is used to run programs or instructions to implement the fibrous grain three-dimensional crystal modeling and analysis method as described in the first aspect.

[0016] In a sixth aspect, the present application provides a computer program product, comprising a computer program, which, when executed by a processor, implements the fibrous grain three-dimensional crystal modeling and analysis method as described in the first aspect above.

[0017] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become obvious from the description below, or will be learned through practice of the present application.

[0018] The present invention provides a three-dimensional crystal modeling and analysis method for fibrous grains, which has the following advantages over the prior art: (1) The present invention processes a fiber grain sample of a metal wire material to obtain a ctf file and imports it into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample; constructs a three-dimensional fibrous grain model according to the grain size distribution data to obtain a geo file including grain geometry information; adopts an Euler angle screening method to screen the same number of Euler angles as the three-dimensional fibrous grain model according to the grain orientation information to obtain Euler angle data; imports the geo file into Gmsh software for meshing to generate an inp file; inputs the Euler angle data and the inp file into ABAQUS to generate independent material properties for each grain, and performs finite element analysis on the fiber grain sample based on the load and boundary conditions set by the ABAQUS software, which can more accurately simulate the mechanical behavior of metal materials under different loads and boundary conditions, effectively improves the simulation accuracy of the mechanical properties of metal materials, and realizes high-precision restoration and efficient simulation of fibrous grain structure through multi-software collaboration and automated processes, providing technical support for the prediction of mechanical properties of metal wire materials, has wide applicability, and can be applied to different working conditions.

[0019] (2) The present invention effectively improves the accuracy and consistency of the grain model by generating a grain ID sequence based on Neper software and arranging the EBSD raw data in descending order of grain size. Combining the Monte Carlo algorithm, the present invention randomly selects the same number of Euler angles as the three-dimensional fibrous grain model and verifies and screens out the Euler angle data that meets the conditions, which helps to more accurately establish the three-dimensional fibrous grain model. This reduces the errors and uncertainties in the model construction process, improves the efficiency of material property simulation, reduces calculation errors, and has higher authenticity.

[0020] (3) The present invention obtains grain information of different cross sections of the fiber grain sample by using the EBSD probe of the scanning electron microscope and converts it into a ctf file including the grain information, which can more accurately obtain the grain orientation and size distribution information, provide higher precision data support for subsequent three-dimensional modeling and finite element analysis, and can more accurately restore the actual morphology of the fibrous grains. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the description of the embodiments in conjunction with the following drawings, in which: Figure 1This is one of the flow charts of the fibrous grain three-dimensional crystal modeling and analysis method provided in the embodiments of the present application; Figure 2 This is the second flow chart of the fibrous grain three-dimensional crystal modeling and analysis method provided in the embodiments of the present application; Figure 3 : is the inverse pole figure of different cross sections of the tungsten filament grains provided in the embodiments of the present application; Figure 4 It is a model diagram provided in an embodiment of the present application that contains grain geometry information and material properties and is assigned boundary conditions; Figure 5 is a schematic diagram of the simulation results provided in the embodiments of the present application; Figure 6 Schematic diagram of the structure of the fibrous grain three-dimensional crystal modeling and analysis device provided in an embodiment of the present application; Figure 7 It is a structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0022] The following will be combined with the accompanying drawings in the embodiments of the present application to clearly describe the technical solutions in the embodiments of the present application. Obviously, the embodiments described are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field are within the scope of protection of this application.

[0023] The terms "first," "second," and the like in the specification and claims of this application are used to distinguish similar objects, and are not used to describe a specific order or precedence. It should be understood that the terms used in this manner are interchangeable where appropriate, so that the embodiments of this application can be implemented in an order other than that illustrated or described herein, and that the objects distinguished by "first," "second," and the like are generally of the same type, and do not limit the number of objects; for example, the first object can be one or more. In addition, the term "and / or" in the specification and claims refers to at least one of the connected objects, and the character " / " generally indicates that the objects connected are in an "or" relationship.

[0024] Below, in conjunction with the accompanying drawings, the fibrous grain three-dimensional crystal modeling and analysis method, fibrous grain three-dimensional crystal modeling and analysis device, electronic device and readable storage medium provided in the embodiments of the present application are described in detail through specific embodiments and their application scenarios.

[0025] The fibrous grain three-dimensional crystal modeling and analysis method can be applied to a terminal, and can be specifically executed by hardware or software in the terminal.

[0026] The terminal includes, but is not limited to, a portable communication device such as a mobile phone or tablet computer having a touch-sensitive surface (e.g., a touch screen display and / or a touch pad). It should also be understood that in some embodiments, the terminal may not be a portable communication device, but rather a desktop computer having a touch-sensitive surface (e.g., a touch screen display and / or a touch pad).

[0027] In the following embodiments, a terminal including a display and a touch-sensitive surface is described. However, it should be understood that the terminal may include one or more other physical user interface devices such as a physical keyboard, a mouse, and a joystick.

[0028] The embodiment of the present application provides a three-dimensional crystal modeling and analysis method for fibrous grains. The execution subject of the three-dimensional crystal modeling and analysis method for fibrous grains can be an electronic device or a functional module or functional entity in the electronic device that can implement the three-dimensional crystal modeling and analysis method for fibrous grains. The electronic devices mentioned in the embodiment of the present application include but are not limited to mobile phones, tablets, computers, cameras and wearable devices, etc. The three-dimensional crystal modeling and analysis method for fibrous grains provided in the embodiment of the present application is explained below using electronic devices as the execution subject as an example.

[0029] Figure 1 This is one of the flow charts of the fibrous grain three-dimensional crystal modeling and analysis method provided in the embodiments of the present application, such as Figure 1 As shown, the fibrous grain three-dimensional crystal modeling and analysis method includes: step 110, step 120, step 130, step 140, step 150, step 160 and step 170.

[0030] Step 110: Obtain a fiber grain sample of a metal wire material, process the fiber grain sample, and obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; In some embodiments, the processing of the fiber grain sample to obtain a ctf file including grain information includes: The grain information of different cross sections of the fiber grain sample is obtained by using the EBSD probe of the scanning electron microscope, wherein the grain information includes grain orientation and size distribution, and EBSD raw data is obtained; The EBSD raw data is converted into a ctf file including grain information.

[0031] It is easy to understand that a fiber grain sample of a metal wire material is obtained. The fiber grain sample includes multiple grains. The grain orientation information and size distribution data of the cross section and longitudinal section of the fiber grain sample are obtained through EBSD (Electron Backscatter Diffraction) of a scanning electron microscope to obtain the EBSD raw data. The EBSD raw data is converted into a ctf file and exported.

[0032] In this embodiment, by using the EBSD probe of the scanning electron microscope to obtain grain information of different sections of the fibrous grain sample and converting it into a ctf file including the grain information, the grain orientation and size distribution information can be obtained more accurately, providing higher-precision data support for subsequent three-dimensional modeling and finite element analysis, and more accurately restoring the actual morphology of the fibrous grains.

[0033] Step 120: import the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample; In some embodiments, the step of importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample includes: Importing the ctf file into the MTEX ​​toolkit in MATLAB, extracting the lognormal distribution parameters of the cross-sectional grain size and the range of the longitudinal grain aspect ratio, and obtaining the grain orientation information and grain size distribution data of the fiber grain sample; The grain orientation information includes the orientation distribution of the longitudinal section, and the grain size distribution data includes the average grain size distribution data of the cross section and the aspect ratio data of the longitudinal section.

[0034] The ctf format file was imported into MATLAB software, and the MTEX ​​toolkit was run in MATLAB to extract the logarithmic normal distribution parameters (mean μ, standard deviation σ) of the cross-sectional grain size and the measured distribution range of the longitudinal grain aspect ratio, thereby obtaining the grain orientation information and grain size distribution data of the fiber grain sample.

[0035] In this embodiment, by importing the ctf file into the MTEX ​​toolkit in MATLAB, the lognormal distribution parameters of the cross-sectional grain size and the range of the longitudinal grain aspect ratio are extracted, which reduces the manual operation process in the material orientation parameter definition process and improves the efficiency and accuracy of fiber grain sample modeling.

[0036] Step 130: construct a three-dimensional fibrous grain model based on the grain size distribution data to obtain a geo file including grain geometry information; It is easy to understand that according to the longitudinal section grain aspect ratio and cross-sectional average size distribution of the fiber grain sample, a geometric model is established in combination with the Neper software editing code. The code content includes but is not limited to: model size, grain size distribution, grain aspect ratio, and a geo file including grain geometric information is obtained.

[0037] It should be noted that when using Neper software to edit the code to establish the geometric model, the cross-sectional grain size (z) distribution function is set to the log-normal distribution function, and the longitudinal grain aspect ratio x:y value range is the actual aspect ratio distribution, and the aspect ratio y:z=1:1, so as to realize the construction of the three-dimensional geometric model of fibrous grains.

[0038] The calculation formula of the cross-sectional grain size distribution function is as follows:

[0039] in, is the cross-sectional grain size distribution function, is the cross-sectional grain size, σ is the standard deviation of the cross-sectional data, μ is the mean of the cross-sectional data, and they satisfy 0.5≤μ≤2, 0.2≤σ≤0.5.

[0040] Step 140: import the geo file into Gmsh software for meshing, and generate an inp file based on the meshed data; In some embodiments, the step of importing the geo file into the Gmsh software for meshing, and generating an inp file based on the meshed data, includes: Import the geo file into Gmsh software, set the maximum grid size to one fifth of the average grain size, and the minimum grid size to one twentieth of the average grain size; Based on the maximum grid size and the minimum grid size, the Delaunay algorithm is used to divide the data in the geo file into tetrahedral grids, and the inp file is generated based on the divided data.

[0041] It is easy to understand that when the geo file is imported into the Gmsh software for meshing, the maximum grid size is set to 1 / 5 of the average grain size, the minimum grid size is set to 1 / 20 of the average grain size, and the Delaunay grid algorithm is used to divide the data in the geo file into tetrahedral grids. The Delaunay grid division algorithm can reduce the total number of grids by 30%-50%, and the grid density in the grain boundary area is 2-3 times that of the inside of the grain.

[0042] In this example, the meshing accuracy and consistency of the grain model were effectively improved by importing the geo file into Gmsh software, setting the grid size based on the average grain size, and applying the Delaunay algorithm to tetrahedral meshing of the data in the geo file to generate an inp file. By controlling the ratio of the grid size to the grain size, the meshing quality was improved, the errors and uncertainties in the meshing process were reduced, the morphology and distribution of the fibrous grains could be controlled, and the efficiency and accuracy of fibrous grain modeling were improved.

[0043] Step 150: Based on the grain orientation information, an Euler angle screening method is used to screen the same number of Euler angles as the number of the three-dimensional fibrous grain models to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters; In some embodiments, the Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters, including: The grain ID sequence is generated based on the grain size using Neper software, and the EBSD raw data are arranged in descending order according to the grain ID sequence; Euler angle data are obtained by randomly selecting the same number of Euler angles as the three-dimensional fibrous grain model using a Monte Carlo algorithm. The Euler angle data file obtained by screening is expanded to a file format that conforms to the MTEX ​​interface and then imported into MTEX ​​to draw orientation difference maps and pole figures, which are compared and verified with the arranged EBSD raw data. The verified Euler angle data is processed to obtain the grain orientation parameters. If the verification fails, it is re-screened until it passes the verification.

[0044] It is easy to understand that based on the EBSD grain orientation information derived from MTEX, the Euler angle screening method is used to screen the Euler angle combination that matches the number of grains in the 3D fibrous grain model. The pole figures are then used to verify whether the orientation misorientation angle deviation and the pole figure intensity distribution error meet the preset thresholds. The Euler angle screening method includes the following steps: (1) Arrange the EBSD raw data in descending order of grain size according to the grain ID sequence generated by Neper; (2) The Monte Carlo algorithm is used to randomly select the Euler angle combination with the same number of grains as the three-dimensional fibrous grain model, expand the content of the screened Euler angle data file to a file format that conforms to the MTEX ​​interface, and then import it into MTEX ​​to draw the orientation difference map and pole figure, and compare and verify it with the original EBSD data (3) The Euler angle that has passed the verification is used as the Euler angle data. If the verification fails, return to (2) and re-screen until the verification error meets the requirements.

[0045] It should be noted that when comparing and verifying with the original EBSD data, the consistency judgment criteria are that the orientation difference angle deviation does not exceed 4° and the grain orientation intensity distribution error in the pole figure does not exceed 8%.

[0046] In this example, the accuracy and consistency of the grain model were effectively improved by generating a grain ID sequence based on Neper software and sorting the raw EBSD data in descending order of grain size. A Monte Carlo algorithm was then used to randomly select the same number of Euler angles as the three-dimensional fibrous grain model and perform verification. This screened out qualified Euler angle data, facilitating a more accurate construction of the three-dimensional fibrous grain model. This reduced errors and uncertainties in the model construction process, improved the efficiency of material property simulation, reduced computational errors, and achieved higher fidelity.

[0047] Step 160: Input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; Furthermore, the verified Euler angle data is combined with the inp file through a Python script, and the inp file is imported into ABAQUS to generate a geometric model. After assembly, the Python script is run to create materials and sections for different grains and assign material properties. The Python script assigns independent material properties to each grain based on the Euler angle data by calling the Material module and Section module of ABAQUS. Other parameters include crystal elastic modulus, Poisson's ratio and slip system parameters, and creates corresponding solid sections. The specific parameters can be adjusted accordingly according to the experimental data.

[0048] Step 170: Setting loads and boundary conditions based on the ABAQUS software, performing finite element analysis on the fiber grain sample according to the independent material properties of each grain, and obtaining analysis results.

[0049] Finally, the UMAT subroutine was loaded into ABAQUS, displacement loads were applied, and boundary conditions were set to perform crystal plasticity finite element analysis. The hardening model parameters in the UMAT subroutine were determined by fitting the true stress-strain curve of the wire material, with a fitting error of ≤5%. The UMAT subroutine path was specified through the material definition interface of ABAQUS.

[0050] For example, the angle between the direction of the applied displacement load and the axial direction of the fibrous grain is 0°, 45° or 90°, and the boundary condition is to fix all degrees of freedom at one end of the model.

[0051] According to the three-dimensional crystal modeling and analysis method of fibrous grains provided in the embodiment of the present application, a ctf file is obtained by processing a fibrous grain sample of a metal wire material and importing it into MATLAB to obtain grain orientation information and grain size distribution data of the fibrous grain sample; a three-dimensional fibrous grain model is constructed based on the grain size distribution data to obtain a geo file including grain geometric information; the Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model based on the grain orientation information to obtain Euler angle data; the geo file is imported into the Gmsh software for meshing to generate an inp file; the Euler angle data and the inp file are input into ABAQUS to generate independent material properties for each grain, and finite element analysis of the fiber grain sample is performed based on the load and boundary conditions set by the ABAQUS software, which can more accurately simulate the mechanical behavior of metal materials under different loads and boundary conditions, effectively improve the simulation accuracy of the mechanical properties of metal materials, and achieve high-precision restoration and efficient simulation of the fibrous grain structure through multi-software collaboration and automated processes, providing technical support for the prediction of mechanical properties of metal wire materials. It has wide applicability and can be applied to different working conditions.

[0052] In some embodiments, constructing a three-dimensional fibrous grain model based on the grain size distribution data to obtain a geo file including grain geometric information includes: Based on the grain size distribution data, a three-dimensional fibrous grain model was constructed using Neper software, wherein the aspect ratio of the longitudinal section of the three-dimensional fibrous grain model ranged from 3:1 to 15:1, the width-to-height ratio was 1:1, and the cross-sectional size of the three-dimensional fibrous grain model obeyed a lognormal distribution of μ=0.5-2 μm and σ=0.2-0.5; Where σ is the standard deviation of the cross-sectional data, and μ is the mean of the cross-sectional data.

[0053] In this example, Neper software was used to construct a three-dimensional fibrous grain model, incorporating the grain geometry and the distribution of longitudinal and cross-sectional dimensions to effectively improve the accuracy and consistency of the grain model. By controlling the model's longitudinal and cross-sectional aspect ratios, width-to-height ratios, and cross-sectional dimensions, the model's authenticity and rationality were enhanced, reducing errors and uncertainties during the model construction process.

[0054] In some embodiments, setting loads and boundary conditions based on the ABAQUS software and performing finite element analysis on the fiber grain sample includes: Loading the UMAT subroutine in the ABAQUS software to define the constitutive model; Setting displacement loads and boundary conditions, wherein the displacement loads include unidirectional tensile loads and shear loads, the angle between the displacement load direction and the axial direction of the fibrous grains is 0°, 45°, or 90°, and the boundary conditions are to fix all degrees of freedom of the nodes at one end of the model and apply the displacement load to the nodes at the other end; Based on the displacement load and boundary conditions, a crystal plasticity finite element analysis is performed on the fiber grain sample to obtain an analysis result.

[0055] Figure 2 This is the second flow chart of the fibrous grain three-dimensional crystal modeling and analysis method provided in the embodiment of the present application, such as Figure 2 As shown in the figure, the prediction of anisotropic mechanical behavior of metal wire materials with an aspect ratio ≥ 3:1 (such as tungsten wire and aluminum alloy wire) specifically includes the following steps: (1) For a given tungsten wire sample with a diameter of 0.39 mm, the EBSD probe of a scanning electron microscope is used to obtain the grain information of the cross section and longitudinal section of the sample and export it as a ctf file. Figure 3 is the inverse pole figure of different cross sections of the tungsten filament grains provided in the embodiment of the present application, such as Figure 3 As shown, the grain information includes grain orientation and size distribution; (2) Import the ctf format file containing the tungsten filament grain information into MATLAB software, and use the command to run the MTEX ​​toolkit to export the grain orientation and size distribution information of different cross sections; (3) Analyze the aspect ratio of the grains in the longitudinal section and the average size distribution of the cross section, and establish a geometric model by editing the code with Neper software. By setting the model size to a = 5 μm, b = 1 μm, c = 1 μm, the grain ratio to 11:1:1, the geometric mean of the cross-sectional grain size to 0.53126, and the geometric standard deviation to 0.353, the three-dimensional geometric model of the fibrous grains is constructed and then exported as a geo file; (4) Import the geo file containing the grain geometry information into the Gmsh software for meshing. Set the maximum mesh size to 1 / 5 of the average grain size and the minimum mesh size to 1 / 20 of the average grain size. Use the Delaunay mesh algorithm to complete the meshing. After the meshing is completed, export the inp file. (5) Using MATLAB code, the 3,421,257 EBSD orientation data exported by MTEX ​​were screened to find the Euler angle combinations with the same number of grains as the 849 grains generated by Neper; (6) Expand the content of the screened Euler angle file to a file format that conforms to the MTEX ​​interface and then import it into MTEX ​​to draw the orientation difference map and pole figure, and compare and verify them with the original EBSD data; (7) Edit the Python script based on the Euler angle data in (6), create materials and cross sections for different grains of the geometric model, and assign the cross sections to the corresponding component areas; (8) Import the inp file into ABAQUS to generate the geometric model. After assembly, run the Python script to create materials and cross sections for different grains and assign material properties; (9) Apply certain loads and boundary conditions to the model established in (8). Figure 4 This is a model diagram provided by the embodiment of the present application, which contains grain geometry information and material properties and is given boundary conditions. The UMAT subroutine path is set for simulation calculation, and finally the simulation results are extracted and analyzed. Figure 5 It is a schematic diagram of the simulation results provided in the embodiments of the present application.

[0056] In this example, by applying a uniaxial tensile or shear load and combining it with boundary conditions that fixate a node at one end of the model, the deformation behavior of a fiber grain sample under different loads was effectively simulated. This enables more accurate prediction of the mechanical response of fiber grain materials under actual working conditions, simulating the mechanical behavior of materials under different conditions, and providing a deeper understanding of the relationship between a material's microstructure and macroscopic properties, providing an important theoretical basis for material research, development, optimization, and application.

[0057] The fibrous grain three-dimensional crystal modeling and analysis method provided in the embodiments of this application can be executed by a fibrous grain three-dimensional crystal modeling and analysis device. In the embodiments of this application, the fibrous grain three-dimensional crystal modeling and analysis device executing the fibrous grain three-dimensional crystal modeling and analysis method is used as an example to illustrate the fibrous grain three-dimensional crystal modeling and analysis device provided in the embodiments of this application.

[0058] The present application also provides a three-dimensional crystal modeling and analysis device for fibrous grains, such as Figure 6 As shown, the fibrous grain three-dimensional crystal modeling and analysis device includes: an acquisition module 610, a first processing module 620, a second processing module 630, a third processing module 640, a fourth processing module 650, a fifth processing module 660 and an analysis module 670.

[0059] An acquisition module 610 is configured to acquire a fiber grain sample of a metal wire material and process the fiber grain sample to obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; A first processing module 620 is used to import the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample; A second processing module 630 is configured to construct a three-dimensional fibrous grain model based on the grain size distribution data, and obtain a geo file including grain geometry information; The third processing module 640 is used to import the geo file into the Gmsh software for grid division, and generate an inp file based on the divided data; A fourth processing module 650 is configured to, based on the grain orientation information, use an Euler angle screening method to screen the same number of Euler angles as the number of the three-dimensional fibrous grain model to obtain Euler angle data, verify the Euler angle data, and process the verified Euler angle data to obtain grain orientation parameters; A fifth processing module 660 is configured to input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; The analysis module 670 is used to set loads and boundary conditions based on the ABAQUS software, perform finite element analysis on the fiber grain sample according to the independent material properties of each grain, and obtain analysis results.

[0060] According to the three-dimensional crystal modeling and analysis method of fibrous grains provided in the embodiment of the present application, a ctf file is obtained by processing a fibrous grain sample of a metal wire material and importing it into MATLAB to obtain grain orientation information and grain size distribution data of the fibrous grain sample; a three-dimensional fibrous grain model is constructed based on the grain size distribution data to obtain a geo file including grain geometric information; the Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model based on the grain orientation information to obtain Euler angle data; the geo file is imported into the Gmsh software for meshing to generate an inp file; the Euler angle data and the inp file are input into ABAQUS to generate independent material properties for each grain, and finite element analysis of the fiber grain sample is performed based on the load and boundary conditions set by the ABAQUS software, which can more accurately simulate the mechanical behavior of metal materials under different loads and boundary conditions, effectively improve the simulation accuracy of the mechanical properties of metal materials, and achieve high-precision restoration and efficient simulation of the fibrous grain structure through multi-software collaboration and automated processes, providing technical support for the prediction of mechanical properties of metal wire materials. It has wide applicability and can be applied to different working conditions.

[0061] The fibrous grain three-dimensional crystal modeling and analysis device provided in the embodiment of the present application can achieve Figures 1 to 5 To avoid repetition, the various processes implemented in the embodiment of the fibrous grain three-dimensional crystal modeling and analysis method are not described here.

[0062] In some embodiments, as Figure 7As shown, an embodiment of the present application further provides an electronic device 700, comprising a processor 701, a memory 702, and a computer program stored in the memory 702 and executable on the processor 701. When the program is executed by the processor 701, each process of the above-mentioned embodiment of the fibrous grain three-dimensional crystal modeling and analysis method is implemented, and the same technical effect can be achieved. To avoid repetition, it will not be described here.

[0063] It should be noted that the electronic devices in the embodiments of the present application include the mobile electronic devices and non-mobile electronic devices mentioned above.

[0064] An embodiment of the present application also provides a non-transitory computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the various processes of the above-mentioned embodiment of the fibrous grain three-dimensional crystal modeling and analysis method, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0065] The processor is the processor in the electronic device in the above embodiment. The readable storage medium includes a computer-readable storage medium, such as a computer read-only memory (ROM), random access memory (RAM), a magnetic disk, or an optical disk.

[0066] An embodiment of the present application also provides a computer program product, including a computer program, which, when executed by a processor, implements the above-mentioned fibrous grain three-dimensional crystal modeling and analysis method.

[0067] The processor is the processor in the electronic device in the above embodiment. The readable storage medium includes a computer readable storage medium, such as a computer read-only memory ROM, a random access memory RAM, a magnetic disk or an optical disk.

[0068] An embodiment of the present application further provides a chip, which includes a processor and a communication interface, which are coupled to the processor. The processor is used to run programs or instructions to implement the various processes of the above-mentioned embodiment of the fibrous grain three-dimensional crystal modeling and analysis method, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0069] It should be understood that the chip mentioned in the embodiments of the present application can also be called a device-level chip, a device chip, a chip device, or an on-chip device chip, etc.

[0070] It should be noted that, in this article, the terms "comprise", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the statement "comprises a ..." does not exclude the presence of other identical elements in the process, method, article or device comprising the element. In addition, it should be noted that the scope of the methods and devices in the embodiments of the present application is not limited to performing functions in the order shown or discussed, and may also include performing functions in a substantially simultaneous manner or in the opposite order according to the functions involved. For example, the described method may be performed in an order different from that described, and various steps may also be added, omitted, or combined. In addition, the features described with reference to certain examples may be combined in other examples.

[0071] Through the description of the above embodiments, those skilled in the art will clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus the necessary general hardware platform. Of course, hardware can also be used, but in many cases the former is a more preferred embodiment. Based on this understanding, the technical solution of this application, or the portion that contributes to the prior art, can be embodied in the form of a computer software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, or optical disk) and includes a number of instructions for enabling a terminal (which can be a mobile phone, computer, server, or network device, etc.) to execute the fibrous grain three-dimensional crystal modeling and analysis methods of each embodiment of this application.

[0072] In the description of this application, "first feature" and "second feature" may include one or more such features.

[0073] In the description of this application, “plurality” means two or more.

[0074] The embodiments of the present application are described above in conjunction with the accompanying drawings, but the present application is not limited to the above-mentioned specific implementation methods. The above-mentioned specific implementation methods are merely illustrative and not restrictive. Under the guidance of this application, ordinary technicians in this field can also make many forms without departing from the purpose of this application and the scope of protection of the claims, all of which are within the protection of this application.

[0075] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "illustrative embodiments," "examples," "specific examples," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0076] Although the embodiments of the present application have been shown and described, those skilled in the art will appreciate that various changes, modifications, substitutions, and variations may be made to the embodiments without departing from the principles and intent of the present application, and that the scope of the present application is defined by the claims and their equivalents.

Claims

1. A three-dimensional crystal modeling and analysis method for fibrous grains, characterized in that: The method comprises: Obtaining a fiber grain sample of a metal wire material, processing the fiber grain sample to obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; Importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample; Based on the grain size distribution data, a three-dimensional fibrous grain model is constructed to obtain a geo file including grain geometry information; Import the geo file into Gmsh software for grid division, and generate an inp file based on the divided data; Based on the grain orientation information, an Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters; Input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; Loads and boundary conditions are set based on the ABAQUS software, and finite element analysis is performed on the fiber grain sample according to the independent material properties of each grain to obtain analysis results.

2. The fibrous grain three-dimensional crystal modeling and analysis method according to claim 1, characterized in that: The fiber grain sample is processed to obtain a ctf file including grain information, including: The grain information of different cross sections of the fiber grain sample is obtained by using the EBSD probe of the scanning electron microscope, wherein the grain information includes grain orientation and size distribution, and EBSD raw data is obtained; The EBSD raw data is converted into a ctf file including grain information.

3. The fibrous grain three-dimensional crystal modeling and analysis method according to claim 1, characterized in that: The step of importing the ctf file into MATLAB to obtain grain orientation information and grain size distribution data of the fiber grain sample includes: Importing the ctf file into the MTEX ​​toolkit in MATLAB, extracting the lognormal distribution parameters of the cross-sectional grain size and the range of the longitudinal grain aspect ratio, and obtaining the grain orientation information and grain size distribution data of the fiber grain sample; The grain orientation information includes the orientation distribution of the longitudinal section, and the grain size distribution data includes the average grain size distribution data of the cross section and the aspect ratio data of the longitudinal section.

4. The method for modeling and analyzing 3D crystals of fibrous grains according to claim 1, wherein: The method of constructing a three-dimensional fibrous grain model based on the grain size distribution data and obtaining a geo file including grain geometric information includes: Based on the grain size distribution data, a three-dimensional fibrous grain model was constructed using Neper software, wherein the aspect ratio of the longitudinal section of the three-dimensional fibrous grain model ranged from 3:1 to 15:1, the width-to-height ratio was 1:1, and the cross-sectional size of the three-dimensional fibrous grain model obeyed a lognormal distribution of μ=0.5-2 μm and σ=0.2-0.5; Where σ is the standard deviation of the cross-sectional data, and μ is the mean of the cross-sectional data.

5. The method for three-dimensional crystal modeling and analysis of fibrous grains according to claim 2, characterized in that: The Euler angle screening method is used to screen the same number of Euler angles as the three-dimensional fibrous grain model to obtain Euler angle data, the Euler angle data is verified, and the verified Euler angle data is processed to obtain grain orientation parameters, including: The grain ID sequence is generated based on the grain size using Neper software, and the EBSD raw data are arranged in descending order according to the grain ID sequence; Euler angle data are obtained by randomly selecting the same number of Euler angles as the three-dimensional fibrous grain model using a Monte Carlo algorithm. The Euler angle data file obtained by screening is expanded to a file format that conforms to the MTEX ​​interface and then imported into MTEX ​​to draw orientation difference maps and pole figures, which are compared and verified with the arranged EBSD raw data. The verified Euler angle data is processed to obtain the grain orientation parameters. If the verification fails, it is re-screened until it passes the verification.

6. The fibrous grain three-dimensional crystal modeling and analysis method according to claim 1, characterized in that: The step of importing the geo file into the Gmsh software for meshing and generating an inp file based on the meshed data includes: Import the geo file into Gmsh software, set the maximum grid size to one fifth of the average grain size, and the minimum grid size to one twentieth of the average grain size; Based on the maximum grid size and the minimum grid size, the Delaunay algorithm is used to divide the data in the geo file into tetrahedral grids, and the inp file is generated based on the divided data.

7. The method for modeling and analyzing fibrous grain three-dimensional crystals according to claim 1, characterized in that: The setting of loads and boundary conditions based on the ABAQUS software and performing finite element analysis on the fiber grain sample include: Loading the UMAT subroutine in the ABAQUS software to define the constitutive model; Setting displacement loads and boundary conditions, wherein the displacement loads include unidirectional tensile loads and shear loads, the angle between the displacement load direction and the axial direction of the fibrous grains is 0°, 45°, or 90°, and the boundary conditions are to fix all degrees of freedom of the nodes at one end of the model and apply the displacement load to the nodes at the other end; Based on the displacement load and boundary conditions, a crystal plasticity finite element analysis is performed on the fiber grain sample to obtain an analysis result.

8. A fibrous grain three-dimensional crystal modeling and analysis device, implemented using the fibrous grain three-dimensional crystal modeling and analysis method according to any one of claims 1 to 7, characterized in that: The device comprises: an acquisition module, configured to acquire a fiber grain sample of a metal wire material, process the fiber grain sample, and obtain a ctf file including grain information, wherein the fiber grain sample includes a plurality of grains; A first processing module is used to import the ctf file into MATLAB to obtain the grain orientation information and grain size distribution data of the fiber grain sample; A second processing module is used to construct a three-dimensional fibrous grain model based on the grain size distribution data to obtain a geo file including grain geometric information; The third processing module is used to import the geo file into the Gmsh software for grid division, and generate an inp file based on the divided data; a fourth processing module, configured to, based on the grain orientation information, use an Euler angle screening method to screen Euler angles having the same number as the three-dimensional fibrous grain model to obtain Euler angle data, verify the Euler angle data, and obtain grain orientation parameters by processing the verified Euler angle data; a fifth processing module, configured to input the grain orientation parameters and the inp file into ABAQUS to generate independent material properties for each grain; The analysis module is used to set loads and boundary conditions based on the ABAQUS software, perform finite element analysis on the fiber grain sample according to the independent material properties of each grain, and obtain analysis results.

9. An electronic 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 program, the fibrous grain three-dimensional crystal modeling and analysis method according to any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for modeling and analyzing fibrous grain three-dimensional crystals according to any one of claims 1 to 7 is implemented.

Images