An acoustic modeling method based on voronoi lamellar dual-phase microstructure

CN121659648BActive Publication Date: 2026-08-11INNER MONGOLIA UNIV OF SCI & TECH
View PDF 0 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2025-12-04
Publication Date
2026-08-11

AI Technical Summary

Technical Problem

而在金属材料的实际微观组织结构中,存在双相的微观结构,在声学模型构建中并未考虑双相微观组织结构的建立

Benefits of technology

本发明通过融合Voronoi晶粒结构与层状第二相建模,突破了传统单相模型的局限,构建了与真实微观组织高度一致的声学模型;利用Comsol与MATLAB联合仿真实现了参数化建模流程,显著提升了建模效率与可控性;该模型能精确模拟超声波在双相材料中的传播行为,为揭示材质劣化与声学响应关联机制、发展无损评价方法提供了有效的技术手段。

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN121659648B_ABST
    Figure CN121659648B_ABST
Patent Text Reader

Abstract

This invention discloses an acoustic modeling method based on a Voronoi layered two-phase microstructure, comprising: generating a Voronoi grain structure according to the average grain size of the polycrystalline metal material; constructing a geometric model of the layered two-phase microstructure by setting the interlayer spacing and area ratio coefficients; assigning different elastic matrices to the grain matrix and the layered second-phase microstructure respectively using the Bond transformation method; establishing an acoustic simulation model for simulating ultrasonic wave propagation by setting the solid mechanics physical field, ultrasonic excitation source, and boundary conditions and performing mesh generation using co-simulation software based on the elastic matrices; and setting a solver to perform calculations based on the acoustic simulation model to obtain and output the propagation characteristics of ultrasonic waves in the layered two-phase microstructure. This invention is of great significance for studying the propagation law of ultrasonic waves in polycrystalline metal materials with a layered two-phase microstructure.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention belongs to the field of non-destructive evaluation technology of material degradation damage, and particularly relates to an acoustic modeling method based on Voronoi layered two-phase microstructure. Background Technology

[0002] High-temperature pressure-bearing components are prone to material degradation damage during service, with common metallic materials including carbon steel and low-alloy steel. Material degradation causes changes in the internal microstructure of metallic materials, including material structure, grain size, and grain shape. Changes in the microstructure lead to degradation of the material's mechanical properties, and in severe cases, can cause serious safety accidents. The aforementioned metallic materials all possess a two-phase microstructure, and changes in this microstructure significantly affect the extraction of ultrasonic characteristic parameters. For example, coarse grain size, carbide precipitation, and decomposition of the second phase can cause changes in backscattered signals and reduced echo amplitude during ultrasonic propagation. Extracting ultrasonic characteristic parameters allows for the evaluation of material degradation damage in metallic materials, ensuring safe equipment operation. Therefore, establishing an acoustic model that accurately describes the layered two-phase microstructure is of great significance for studying the propagation laws of ultrasound in polycrystalline metals and developing acoustic evaluation methods for material degradation damage.

[0003] Currently, acoustic models of polycrystalline metallic materials mostly use Voronoi diagrams to establish grain structures, and are largely based on single-phase microstructures, with corresponding physical fields added for simulation studies. However, the actual microstructure of metallic materials contains two-phase microstructures, which are not considered in the construction of acoustic models.

[0004] Therefore, there is an urgent need to propose an acoustic modeling method based on the Voronoi layered two-phase microstructure. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention proposes an acoustic modeling method based on Voronoi layered two-phase microstructure. This method introduces a layered microstructure within the grains, building upon the Voronoi grain structure. This method is of great significance for studying the propagation laws of ultrasound in polycrystalline metal materials with layered two-phase microstructures and for developing acoustic evaluation methods for material degradation damage.

[0006] This invention provides an acoustic modeling method based on Voronoi layered two-phase microstructure, comprising the following steps: Based on the average grain size of polycrystalline metal materials, a Voronoi grain structure is generated; Based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficient, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model. Based on the described layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure respectively by the Bond transformation method. Based on the elasticity matrix, a solid mechanical physical field, ultrasonic excitation source, and boundary conditions are set using co-simulation software, and a mesh is generated to establish an acoustic simulation model for simulating ultrasonic wave propagation. Based on the acoustic simulation model, a solver is set up to perform calculations to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure.

[0007] Optionally, a Voronoi grain structure is generated based on the average grain size of the polycrystalline metal material, specifically including: Metallographic testing was performed on polycrystalline metal material samples to obtain the average grain size. The Voronoi function in MATLAB is used to generate a grain structure with a preset grain diameter, and the area distribution of the generated grains is calculated.

[0008] Optionally, based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficients, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model, specifically including: Based on the set interlayer spacing parameters and area ratio coefficient, obtain the logical judgment conditions used to distinguish the matrix from the second phase; Based on the aforementioned logical judgment conditions, the geometric coordinates inside the grain are determined to obtain the spatial distribution of the matrix and the second phase; Based on the spatial distribution of the matrix and the second phase, a layered microstructure is generated inside the grain to construct the final layered two-phase microstructure geometric model.

[0009] Optionally, based on the layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure respectively using the Bond transformation method, specifically including: Based on the Euler angles corresponding to each grain, obtain the rotation matrix of each grain in the global coordinate system; Based on the rotation matrix and the known elastic stiffness matrix in the crystallographic coordinate system, the anisotropic elastic matrix in the global coordinate system is obtained by calculation using the Bond transformation method. Based on the anisotropic elastic matrix, material properties are assigned to the grain matrix and the layered second phase microstructure in the model, respectively.

[0010] Optionally, based on the elasticity matrix, a solid mechanical physical field, ultrasonic excitation source, and boundary conditions are set using co-simulation software, and a mesh is generated to establish an acoustic simulation model for simulating ultrasonic wave propagation, specifically including: Based on the elasticity matrix and material property parameters, select the solid mechanics physical field and assign values ​​to the material properties in the co-simulation software; Based on the selected physical field and ultrasonic excitation frequency, the geometric model of the layered biphasic microstructure is divided into free triangular meshes, and the maximum unit size is set to be between one-twentieth and one-tenth of the ultrasonic wave length. Based on the requirements of ultrasonic simulation, a force source along the normal direction is set at the top of the model as the ultrasonic excitation source. Based on the model boundary constraints, set free boundary conditions at the bottom of the model and symmetrical boundary conditions on the left and right sides of the model. Based on the mesh size and ultrasonic propagation characteristics, a simulation time step matching the mesh is set and transient study parameters are configured to establish an acoustic simulation model for simulating ultrasonic wave propagation.

[0011] Optionally, based on the acoustic simulation model, a solver is set to perform calculations to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure, specifically including: Based on the physical field settings and mesh generation of the acoustic simulation model, select the MUMPS solver and set the transient study parameters to start the model calculation; Based on the time-domain data calculated from the model, probes are set in the model to monitor the displacement response at specific locations; Based on the displacement response data, a time-domain waveform diagram of ultrasonic wave propagation is plotted and output using a one-dimensional plotting group; Based on the displacement response data, the peak value of the frequency domain variance of the ultrasonic backscatter signal is calculated and output through frequency domain analysis.

[0012] This invention also proposes an acoustic modeling system based on Voronoi layered two-phase microstructure for implementing the method, comprising: The structure generation module is used to generate Voronoi grain structures based on the average grain size of polycrystalline metal materials. The model building module is used to generate a layered second-phase microstructure inside the grain based on the Voronoi grain structure by setting the interlayer spacing and area ratio coefficient, so as to construct a layered two-phase microstructure geometric model. The material assignment module is used to assign different elastic matrices to the grain matrix and the layered second-phase microstructure respectively according to the geometric model of the layered two-phase microstructure using the Bond transformation method. The simulation modeling module is used to set the solid mechanical physical field, ultrasonic excitation source and boundary conditions and perform mesh generation using co-simulation software based on the elastic matrix, and to establish an acoustic simulation model for simulating ultrasonic wave propagation. The propagation output module is used to set up a solver to perform calculations based on the acoustic simulation model, so as to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure.

[0013] The present invention also proposes a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method.

[0014] The present invention also proposes a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method.

[0015] The present invention also proposes a computer program product, including a computer program that, when executed by a processor, implements the steps of the method.

[0016] Compared with the prior art, the present invention has the following advantages and technical effects: This invention overcomes the limitations of traditional single-phase models by integrating Voronoi grain structure and layered second-phase modeling, constructing an acoustic model that is highly consistent with the actual microstructure. The parametric modeling process is realized through joint simulation using Comsol and MATLAB, significantly improving modeling efficiency and controllability. This model can accurately simulate the propagation behavior of ultrasonic waves in two-phase materials, providing an effective technical means for revealing the correlation mechanism between material degradation and acoustic response, and for developing non-destructive evaluation methods. Attached Figure Description

[0017] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings: Figure 1 This is a schematic diagram of the method flow according to an embodiment of the present invention; Figure 2 This is a schematic diagram of the Voronoi grain structure model according to an embodiment of the present invention; Figure 3 This is a schematic diagram of the layered microstructure model inside the grains according to an embodiment of the present invention; Figure 4 This is a schematic diagram illustrating the excitation source and boundary conditions for the grains in an embodiment of the present invention; Figure 5 This is a schematic diagram of the output waveform at a certain point on the top of the grain model in an embodiment of the present invention; Figure 6This is a schematic diagram of the peak variance of the backscattered radio frequency domain in the ultrasound of a layered biphasic microstructure according to an embodiment of the present invention. Detailed Implementation

[0018] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0019] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0020] Example 1 This embodiment proposes an acoustic modeling method based on Voronoi layered two-phase microstructure. Addressing the two-phase microstructure problem in polycrystalline metal materials, it combines the Voronoi algorithm to construct a geometric model of the layered two-phase microstructure, achieving an accurate description of the microstructure in material degradation damage. This avoids the limitation of traditional Voronoi diagrams, which can only construct acoustic models with single-phase random grain orientations. This provides a solution for studying the propagation laws of ultrasound in polycrystalline metal materials and developing acoustic evaluation methods for material degradation damage. Specifically, it includes the following steps: Based on the average grain size of polycrystalline metal materials, a Voronoi grain structure is generated; Based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficient, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model. Based on the described layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure respectively by the Bond transformation method. Based on the elasticity matrix, a solid mechanical physical field, ultrasonic excitation source, and boundary conditions are set using co-simulation software, and a mesh is generated to establish an acoustic simulation model for simulating ultrasonic wave propagation. Based on the acoustic simulation model, a solver is set up to perform calculations to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure.

[0021] As a specific implementation method, the acoustic modeling method based on Voronoi layered two-phase microstructure proposed in this embodiment involves metallographic testing of polycrystalline metal material samples and calculation of the average grain size; generating a geometric model with a specified grain size using the Voronoi function in MATLAB software; generating the layered microstructure inside the grain by setting the interlayer spacing and area ratio coefficients; calculating the elastic matrix of the grain and layered structure using the Bond transform method, and assigning the elastic matrix to the grain and layered microstructure; and using the Comsol with MATLAB co-simulation software, selecting the solid mechanics physical field, setting the corresponding boundary conditions, and completing the establishment of the acoustic model.

[0022] It is feasible to generate a Voronoi grain structure based on the average grain size of the polycrystalline metal material, specifically including: Metallographic testing was performed on polycrystalline metal samples to obtain the average grain size; the Voronoi function in MATLAB was used to generate grain structures with preset grain diameters, and the area distribution of the generated grains was calculated.

[0023] Specifically, metallographic testing was performed on polycrystalline metal samples using an optical microscope. Grains were randomly intercepted using straight lines. If the total length of the measurement line is L and the number of grains intercepted is N, then for continuously distributed single-phase austenite grains, the average intercept length is L1 = L / N. L1 is the intercept length of a three-dimensional object; it refers to the average intercept length within the object when randomly intercepting it. When the number of measurements is sufficiently large, the intercept length L2 of a two-dimensional section equals L1. The Voronoi function in MATLAB was used to generate grain structures with a specified grain diameter, and the area distribution of the generated grains was calculated to determine whether it follows a normal distribution.

[0024] Implementable, based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficients, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model, specifically including: Based on the set interlayer spacing parameters and area ratio coefficients, logical judgment conditions for distinguishing the matrix and the second phase are obtained; based on the logical judgment conditions, the geometric coordinates inside the grain are determined to obtain the spatial distribution of the matrix and the second phase; based on the spatial distribution of the matrix and the second phase, a layered microstructure is generated inside the grain to construct the final layered two-phase microstructure geometric model.

[0025] Specifically, the interlayer spacing d and the area ratio coefficient a are set, the inclination angle of the straight line is θ, and the matrix and the second phase are determined by the logical expression mod(|sin(θ)*x+con(θ)*y|,d) / d≥a, thus generating the layered microstructure inside the grain.

[0026] Implementable, based on the aforementioned layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure respectively using the Bond transformation method, specifically including: Based on the Euler angles corresponding to each grain, obtain the rotation matrix of each grain in the global coordinate system; based on the rotation matrix and the known elastic stiffness matrix in the crystallographic coordinate system, calculate the anisotropic elastic matrix in the global coordinate system using the Bond transformation method; based on the anisotropic elastic matrix, assign material property values ​​to the grain matrix and the layered second phase microstructure in the model, respectively.

[0027] Specifically, in this embodiment, the elastic constants of the grains are set in the material module, and the anisotropic elastic matrix of the grains and layered structure in the global coordinate system is calculated using the Bond transformation method, C'=RCR. T C is the Voigt elastic stiffness matrix of a single grain in the crystallographic coordinate system, and R is the rotation matrix corresponding to each grain, using Euler angles as the rotation matrix a(θ1, φ1, ..., ...) ) can be calculated: (1) Based on this, the rotation matrix R(θ1, φ1, ...) corresponding to each grain is calculated. (2) The anisotropic elastic matrix of grains and layered structures in the global coordinate system was calculated using the Bond transformation method. (3) In the formula, C is the Voigt elastic stiffness matrix of a single grain in the crystallographic coordinate system, and C' is the elastic stiffness matrix of the grain after rotation.

[0028] Implementably, based on the elasticity matrix, a co-simulation software is used to set the solid mechanical physical field, ultrasonic excitation source, and boundary conditions, and to perform mesh generation to establish an acoustic simulation model for simulating ultrasonic wave propagation. Specifically, this includes: Based on the elastic matrix and material property parameters, a solid mechanics physical field is selected in the co-simulation software, and the material properties are assigned values. According to the selected physical field and the ultrasonic excitation frequency, the layered two-phase microstructure geometric model is divided into free triangular meshes, with the maximum element size set between one-twentieth and one-tenth of the ultrasonic wave length. Based on the ultrasonic simulation requirements, a force source along the normal direction is set at the top of the model as the ultrasonic excitation source. According to the model boundary constraints, free boundary conditions are set at the bottom of the model, and symmetrical boundary conditions are set on the left and right sides of the model. Based on the mesh size and ultrasonic propagation characteristics, a simulation time step matching the mesh is set, and transient research parameters are configured to establish an acoustic simulation model for simulating ultrasonic wave propagation.

[0029] Specifically, using the Comsol with MATLAB co-simulation software, a solid mechanics physical field was selected, and rotation matrices a(θ1, φ1, ..., ...) were set for different grains. Anisotropic elasticity matrix was set. Material properties and excitation frequency parameters were set, and a free triangular mesh was generated for the layered two-phase microstructure model. The maximum element size was set between 1 / 20 and 1 / 10 of the wavelength. A force source was set at the top, the bottom was set as a free boundary condition, and the left and right sides were set as symmetrical boundary conditions. A simulation step size matching the mesh size was set. A transient study was set up in the study, and the simulation time step was set for simulation calculation.

[0030] Implementably, based on the acoustic simulation model, a solver is set up to perform calculations to obtain and output the propagation characteristics of ultrasound waves in the layered biphasic microstructure, specifically including: Based on the physical field settings and mesh generation of the acoustic simulation model, the MUMPS solver is selected and transient study parameters are set to start the model calculation; based on the time-domain data obtained from the model calculation, probes are set in the model to monitor the displacement response at specific locations; based on the displacement response data, the time-domain waveform of ultrasonic wave propagation is plotted and output using a one-dimensional plotting group; based on the displacement response data, the frequency domain variance peak value of the ultrasonic backscatter signal is calculated and output through frequency domain analysis.

[0031] Specifically, the MUMPS solver in Comsol software was used to calculate the layered two-phase microstructure model, and the results were saved. Probes were added to the results to monitor the displacement of points in the simulation. One-dimensional and two-dimensional plotting groups were set up to plot the displacement of ultrasonic simulation points, and the relevant data were exported. The construction of this acoustic simulation model provides support for the testing of high-temperature pressure-bearing components that have experienced material degradation damage during service. It can provide technical support for the acoustic detection of material degradation damage samples and enable effective research on the propagation laws under different damage levels.

[0032] This embodiment establishes an effective geometric model in the acoustic simulation by constructing a grain model based on the Voronoi algorithm and using interlayer spacing and area ratio coefficients to establish the microstructure of the second phase within the grain. The use of Comsol with MATLAB co-simulation software simplifies the modeling process and enables parametric modeling of different simulation models. The simulation model allows for better study of the propagation of ultrasound in layered two-phase microstructures, achieving acoustic evaluation of material degradation and damage.

[0033] On the other hand, this embodiment also proposes an acoustic modeling system based on Voronoi layered two-phase microstructure for implementing the method, including: The structure generation module is used to generate Voronoi grain structures based on the average grain size of polycrystalline metal materials. The model building module is used to generate a layered second-phase microstructure inside the grain based on the Voronoi grain structure by setting the interlayer spacing and area ratio coefficient, so as to construct a layered two-phase microstructure geometric model. The material assignment module is used to assign different elastic matrices to the grain matrix and the layered second-phase microstructure respectively according to the geometric model of the layered two-phase microstructure using the Bond transformation method. The simulation modeling module is used to set the solid mechanical physical field, ultrasonic excitation source and boundary conditions and perform mesh generation using co-simulation software based on the elastic matrix, and to establish an acoustic simulation model for simulating ultrasonic wave propagation. The propagation output module is used to set up a solver to perform calculations based on the acoustic simulation model, so as to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure.

[0034] Example 2 Figure 1 This is a flowchart of an acoustic modeling method based on Voronoi layered two-phase microstructure in this embodiment. The specific implementation method is as follows: (1) The Voronoi function in MATLAB software is used to generate a grain structure of a specified size. A seed point is randomly arranged within a small square with side length a, and control coefficients are used to adjust the grain shape so that the overall grain area follows a normal distribution. The schematic diagram of the generated grain model is shown below. Figure 2 As shown.

[0035] (2) Set the interlayer spacing and area ratio coefficients, and generate the layered microstructure within the grain using logical expressions. The set interlayer spacing and area ratio coefficients are used to distinguish between grains and layered microstructures. A layered microstructure with a fixed spacing or area ratio is generated within the selected grain, as shown in the following diagram. Figure 3 As shown.

[0036] (3) Calculate the elasticity matrices of grains and layered structures using the Bond transform method. The elasticity matrices are calculated using the elastic constants of different phases from existing literature. The Bond transform method is used to calculate the elasticity matrices of grains and layered structures, and values ​​are assigned to different grains and layered structures using MATLAB software.

[0037] (4) Using Comsol with MATLAB co-simulation software, boundary conditions were set for simulation analysis. The solid mechanics physical field was invoked through the co-simulation software, applying a force source along the negative y-axis. Free boundary conditions were set at the bottom, and symmetrical boundary conditions were set on the left and right sides. The simulation time was set to T=4μs, the excitation frequency f=10MHz, and the maximum mesh element size was λ / 10. Figure 4 As shown.

[0038] (5) Set up the MUMPS solver to calculate the model, extract the displacement of a point at a vertex, and the result is as follows: Figure 5 As shown, save the corresponding calculation results.

[0039] (6) Repeating steps (2) to (5) yields acoustic simulation models with different interlayer spacing and area ratios. Extracting different acoustic parameters allows for the evaluation of the layered two-phase microstructure model. The extracted backscattered signal frequency domain variance peak value is shown below. Figure 6 As shown, the acoustic model of the layered two-phase microstructure of polycrystalline metal materials proposed in this embodiment can be used for the study of ultrasonic propagation mechanisms and the evaluation of material degradation and damage.

[0040] Example 3 This embodiment also discloses a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method described in Embodiment 1.

[0041] Example 4 This embodiment also discloses a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the steps of the method described in Embodiment 1.

[0042] Example 5 This embodiment also discloses a computer program product, including a computer program that, when executed by a processor, implements the steps of the method described in Embodiment 1.

[0043] The above are merely preferred embodiments of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. An acoustic modeling method based on Voronoi lamellar dual-phase microstructure, characterized in that, Includes the following steps: Based on the average grain size of polycrystalline metal materials, a Voronoi grain structure is generated; Based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficient, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model. Based on the described layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure respectively by the Bond transformation method. Based on the elasticity matrix, a solid mechanical physical field, ultrasonic excitation source, and boundary conditions are set using co-simulation software, and a mesh is generated to establish an acoustic simulation model for simulating ultrasonic wave propagation. Based on the acoustic simulation model, a solver is set up to perform calculations to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure. Based on the Voronoi grain structure, by setting the interlayer spacing and area ratio coefficients, a layered second-phase microstructure is generated inside the grain to construct a layered two-phase microstructure geometric model, specifically including: Based on the set interlayer spacing parameters and area ratio coefficient, obtain the logical judgment conditions used to distinguish the matrix from the second phase; Set the interlayer spacing d and the area ratio coefficient a, and the straight line tilt angle θ. Use the logical expression mod(|sin(θ)*x+con(θ)*y|,d) / d≥a to determine the matrix and the second phase, and generate the layered microstructure inside the grain. Here, x represents the abscissa value of a point inside the grain on the two-dimensional plane, and y represents the ordinate value of a point inside the grain on the two-dimensional plane. Based on the aforementioned logical judgment conditions, the geometric coordinates inside the grain are determined to obtain the spatial distribution of the matrix and the second phase; Based on the spatial distribution of the matrix and the second phase, a layered microstructure is generated inside the grain to construct the final layered two-phase microstructure geometric model.

2. The acoustic modeling method of claim 1, wherein, Based on the average grain size of polycrystalline metal materials, a Voronoi grain structure is generated, specifically including: Metallographic testing was performed on polycrystalline metal material samples to obtain the average grain size. The Voronoi function in MATLAB is used to generate a grain structure with a preset grain diameter, and the area distribution of the generated grains is calculated.

3. The acoustic modeling method of claim 1, wherein, Based on the described layered two-phase microstructure geometric model, different elastic matrices are assigned to the grain matrix and the layered second-phase microstructure using the Bond transformation method, specifically including: Based on the Euler angles corresponding to each grain, obtain the rotation matrix of each grain in the global coordinate system; Based on the rotation matrix and the known elastic stiffness matrix in the crystallographic coordinate system, the anisotropic elastic matrix in the global coordinate system is obtained by calculation using the Bond transformation method. Based on the anisotropic elastic matrix, material properties are assigned to the grain matrix and the layered second phase microstructure in the model, respectively.

4. The acoustic modeling method of claim 1, wherein, Based on the elasticity matrix, a solid mechanical physical field, ultrasonic excitation source, and boundary conditions are set using co-simulation software, and a mesh is generated to establish an acoustic simulation model for simulating ultrasonic wave propagation. Specifically, this includes: Based on the elasticity matrix and material property parameters, select the solid mechanics physical field and assign values ​​to the material properties in the co-simulation software; Based on the selected physical field and ultrasonic excitation frequency, the geometric model of the layered biphasic microstructure is divided into free triangular meshes, and the maximum unit size is set to be between one-twentieth and one-tenth of the ultrasonic wave length. Based on the requirements of ultrasonic simulation, a force source along the normal direction is set at the top of the model as the ultrasonic excitation source. Based on the model boundary constraints, set free boundary conditions at the bottom of the model and symmetrical boundary conditions on the left and right sides of the model. Based on the mesh size and ultrasonic propagation characteristics, a simulation time step matching the mesh is set and transient study parameters are configured to establish an acoustic simulation model for simulating ultrasonic wave propagation.

5. The acoustic modeling method of claim 1, wherein, Based on the acoustic simulation model, a solver is set up to perform calculations to obtain and output the propagation characteristics of ultrasound waves in the layered biphasic microstructure, specifically including: Based on the physical field settings and mesh generation of the acoustic simulation model, select the MUMPS solver and set the transient study parameters to start the model calculation; Based on the time-domain data calculated from the model, probes are set in the model to monitor the displacement response at specific locations; Based on the displacement response, a time-domain waveform diagram of ultrasonic wave propagation is plotted and output using a one-dimensional plotting group; Based on the displacement response, the peak value of the frequency domain variance of the ultrasonic backscatter signal is calculated and output through frequency domain analysis.

6. An acoustic modeling system based on Voronoi lamellar dual phase microstructure, characterized in that, For implementing the method according to any one of claims 1-5, comprising: The structure generation module is used to generate Voronoi grain structures based on the average grain size of polycrystalline metal materials. The model building module is used to generate a layered second-phase microstructure inside the grain based on the Voronoi grain structure by setting the interlayer spacing and area ratio coefficient, so as to construct a layered two-phase microstructure geometric model. The material assignment module is used to assign different elastic matrices to the grain matrix and the layered second-phase microstructure respectively according to the geometric model of the layered two-phase microstructure using the Bond transformation method. The simulation modeling module is used to set the solid mechanical physical field, ultrasonic excitation source and boundary conditions and perform mesh generation using co-simulation software based on the elastic matrix, and to establish an acoustic simulation model for simulating ultrasonic wave propagation. The propagation output module is used to set up a solver to perform calculations based on the acoustic simulation model, so as to obtain and output the propagation characteristics of ultrasonic waves in the layered biphasic microstructure.

7. A computer device comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1-5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1-5.

9. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1-5.