Resonance frequency prediction method, electronic equipment, storage medium and device

Abstract

The invention discloses a resonant frequency prediction method, electronic equipment, a storage medium and a device. The method comprises the following steps: constructing a low-frequency dynamic elastic modulus measurement equipment digital model; performing mesh generation on the digital model; establishing a control equation based on multi-field constraints; loading the control equation to the digital model after mesh generation; loading the boundary constraint condition and the contact boundary condition to the digital model loaded with the control equation; and solving the digital model loaded with the boundary constraint condition and the contact boundary condition, and obtaining the resonance frequency and the corresponding resonance mode of the digital model. According to the method, the control equation, the boundary constraint condition and the contact boundary condition are loaded to the digital model, the working condition of the low-frequency dynamic elastic modulus testing equipment is simulated, the resonant frequency of the low-frequency dynamic elastic modulus testing equipment can be accurately predicted, and therefore the specific influence of the resonant frequency on the measurement result can be effectively judged; and a theoretical basis is provided for constructing a low-frequency elastic modulus test state with a wider frequency band.

Description

Technical Field

[0001] The present invention belongs to the field of geophysical technologies, and more specifically, relates to a resonance frequency prediction method, an electronic device, a storage medium, and a device. Background Art

[0002] In the field of geophysics, especially in the measurement of seismic wave dispersion and attenuation, the stress-strain measurement method is regarded as the most commonly used and crucial measurement means. The measurement device developed based on the stress-strain method is called a low-frequency dynamic elastic modulus measurement device. To accurately obtain data, developing and improving the relevant low-frequency dynamic elastic modulus measurement device is a task that cannot be ignored. However, although low-frequency devices have been developed in the past, many devices are adversely affected by resonance in practical applications, resulting in contaminated and distorted data. Especially when the measurement frequency exceeds a certain range, such as 200 - 400 Hz, the device will generate inaccurate data, thus unable to provide effective constraints for the measurement of seismic wave dispersion and attenuation. Sun et al. (2018) conducted a resonance frequency analysis and prediction on the low-frequency dynamic elastic modulus measurement device developed by Batzle et al., and preliminarily predicted the resonance frequency and resonance mode. However, their research did not consider relevant factors such as prestress and ceramic vibration sources, which limited the prediction accuracy of the resonance frequency. Therefore, in such a background, it is particularly important to construct a method that can directly predict the resonance frequency of low-frequency devices and consider key factors such as prestress and ceramic vibration sources. This can not only help researchers optimize the performance of the device before manufacturing, but also provide more accurate and reliable measurement data for the entire field of geophysics. Considering this demand, we propose this patent, aiming to fill the gap in the existing technology and contribute to the continuous development of the field of geophysics. The information disclosed in the background art section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as an admission or any form of implication that this information constitutes the prior art known to those skilled in the art. Summary of the Invention

[0003] The object of the present invention is to propose a resonance frequency prediction method, an electronic device, a storage medium, and a device, so as to accurately predict the resonance frequency of a low-frequency dynamic elastic modulus measurement device, and thus be able to effectively judge the specific influence of the resonance frequency on the measurement result.

[0004] To achieve the above object, the present invention proposes a resonance frequency prediction method, an electronic device, a storage medium, and a device.

[0005] According to the first aspect of the present invention, a resonance frequency prediction method is proposed, including:

[0006] Build a digital model of the low-frequency dynamic elastic modulus measurement device;

[0007] Perform mesh generation on the digital model;

[0008] Establish a control equation based on multi-field constraints;

[0009] Load the control equation into the digital model that has completed the mesh generation;

[0010] Load the boundary constraint conditions and the contact boundary conditions into the digital model that has already loaded the control equation;

[0011] Solve the digital model that has loaded the boundary constraint conditions and the contact boundary conditions to obtain the resonance frequency and the corresponding resonance mode of the digital model.

[0012] Optionally, the multi-field constraints include:

[0013] Solid mechanics field constraint, fluid field constraint, and electric field constraint.

[0014] Optionally, the control equation includes:

[0015]

[0016]

[0017]

[0018]

[0019] where λ is the eigenvalue, u s is the displacement of the solid material, is the divergence operator, σ is the stress tensor, e is the coupling tensor matrix of the piezoelectric material, ε 0 is the air dielectric constant, ε r is the relative dielectric parameter of the material, V is the voltage, C is the stiffness tensor, : is the tensor operator, is the gradient operator, α is the Biot coefficient, ρ b is the density of the solid material, ρ f is the fluid density, ρ c (λ) is the coupling density, specifically τ is the pore tortuosity of the rock sample, φ is the porosity, η is the fluid viscosity, κ is the permeability of the rock sample.

[0020] Optionally, loading the boundary constraint conditions includes:

[0021] Apply a specified displacement u s = disp;

[0022] The fixed boundary condition is loaded at the bottom of the device, and the displacement u s = 0;

[0023] The boundary condition of the rock sample is loaded with an undrained boundary condition:

[0024] The voltages loaded at both ends of the piezoelectric vibration source are respectively: V = 0, V = V0;

[0025] Among them, disp refers to the displacement, and V0 is the high voltage value.

[0026] Optionally, the contact boundary condition includes:

[0027]

[0028] Among them, p n is the penalty function factor, g n is the contact gap function, p 0 is the pressure when the contact gap is 0, If means if, else means or, that is, if g n <p o / p n , then T n = -p n g n + p 0 , otherwise, T n = 0.

[0029] Optionally, the finite element method is used to solve the digital model loaded with the boundary constraint condition and the contact boundary condition.

[0030] Optionally, the finite element method includes:

[0031] ARPACK.

[0032] According to the second aspect of the present invention, a resonance frequency prediction device is proposed, including:

[0033] A construction module for constructing a digital model of a low-frequency dynamic elastic modulus measurement device;

[0034] A mesh generation module for performing mesh generation on the digital model;

[0035] An establishment module for establishing a control equation based on multi-field constraints;

[0036] A first loading module for loading the control equation into the digital model that has completed the mesh generation;

[0037] A second loading module, configured to load boundary constraint conditions and contact boundary conditions into the digital model that has loaded the control equation;

[0038] A solving module, configured to solve the digital model that has loaded the boundary constraint conditions and the contact boundary conditions, and obtain the resonance frequency and the corresponding resonance mode of the digital model.

[0039] According to a third aspect of the present invention, an electronic device is provided, and the electronic device includes:

[0040] At least one processor; and,

[0041] A memory communicatively connected to the at least one processor; wherein,

[0042] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the resonance frequency prediction method according to any one of the first aspects.

[0043] According to a fourth aspect of the present invention, a non-transitory computer-readable storage medium is provided, and the non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the resonance frequency prediction method according to any one of the first aspects.

[0044] The beneficial effects of the present invention are as follows: Based on the multi-field coupling theory, a control equation is established and loaded into the digital model of the low-frequency dynamic elastic modulus measurement device, and then the boundary constraint conditions and the contact boundary conditions are loaded into the digital model of the low-frequency dynamic elastic modulus measurement device to simulate the working conditions of the low-frequency dynamic elastic modulus test device, and further obtain the resonance frequency and the corresponding resonance mode of the low-frequency dynamic elastic modulus measurement device; the present invention can accurately predict the resonance frequency of the low-frequency dynamic elastic modulus measurement device, so as to effectively judge the specific influence of the resonance frequency on the measurement result, provide a theoretical basis for constructing a low-frequency elastic modulus test state with a wider frequency band, provide strong technical support for manufacturing a more efficient and stable low-frequency dynamic elastic modulus measurement device, optimize the device design scheme, greatly save the device construction cost, and promote the technological progress in the field of geophysics, especially in the measurement of seismic wave dispersion and attenuation.

[0045] The system of the present invention has other characteristics and advantages, which will be obvious from the accompanying drawings incorporated herein and the subsequent detailed description, or will be described in detail in the accompanying drawings incorporated herein and the subsequent detailed description, and these drawings and detailed description are used together to explain the specific principles of the present invention. Description of the Drawings

[0046] The above and other objects, features, and advantages of the present invention will become more apparent by describing the exemplary embodiments of the present invention in more detail with reference to the accompanying drawings. In the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0047] Figure 1 A flowchart showing the steps of a resonance frequency prediction method according to the present invention is presented.

[0048] Figure 2 a, Figure 2 b, Figure 2 c and Figure 2 d respectively show schematic diagrams of the digital model of the low-frequency dynamic elastic modulus measurement device, the digital model of the low-frequency dynamic elastic modulus measurement device after mesh division, the mode of the digital model of the low-frequency dynamic elastic modulus measurement device at non-resonant frequencies, and the mode of the digital model of the low-frequency dynamic elastic modulus measurement device at resonant frequencies according to Embodiment 2 of the present invention.

[0049] Figure 3 A schematic diagram showing the variation of the solid displacement at the center point of the digital model of the low-frequency dynamic elastic modulus measurement device with frequency according to Embodiment 2 of the present invention is presented.

[0050] Figure 4 a, 4b, 4c, and 4d respectively show schematic diagrams of the Young's modulus, Young's modulus attenuation, Poisson's ratio of the Young's modulus, and Poisson phase difference of the Young's modulus of tempered glass according to Embodiment 2 of the present invention.

[0051] Figure 5 A schematic diagram of a resonance frequency prediction device according to Embodiment 3 of the present invention is presented. Detailed implementation manners

[0052] The present invention will be described in more detail below with reference to the accompanying drawings. Although the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. Instead, these embodiments are provided to make the present invention more thorough and complete, and to fully convey the scope of the present invention to those skilled in the art.

[0053] As Figure 1 shown, a resonance frequency prediction method according to the present invention includes:

[0054] Construct a digital model of a low-frequency dynamic elastic modulus measurement device;

[0055] Perform mesh division on the digital model;

[0056] Establish a control equation based on multi-field constraints;

[0057] Load the control equations into the digital model with the mesh generation completed;

[0058] Load the boundary constraint conditions and contact boundary conditions into the digital model with the control equations already loaded;

[0059] Solve the digital model with the boundary constraint conditions and contact boundary conditions loaded to obtain the resonance frequency and the corresponding resonance modes of the digital model.

[0060] Specifically, the present invention first constructs a digital model of the low-frequency dynamic elastic modulus measurement device according to the device size and materials, where the material parameters in the model are obtained through experimental measurements and the device's factory parameters; then, according to the characteristics of the digital model of the low-frequency dynamic elastic modulus measurement device, technical methods such as hexahedron meshing, tetrahedron meshing, and scanning are used to perform mesh generation on the digital model of the low-frequency dynamic elastic modulus measurement device; based on multi-field constraints, control equations are established. The physical laws of the low-frequency dynamic elastic modulus measurement device are mainly constrained by three fields, namely the solid mechanics field constraint, the fluid field, and the electric field. Among them, the solid mechanics field uses the stress-strain equation to constrain the solid materials in the device, including but not limited to the device's skeleton, rock samples, etc.; the fluid field mainly uses Darcy's law to describe the fluid flow inside the rock samples, and the electric field mainly uses the electrostatic field to describe the working electric field of the piezoelectric ceramic vibration source; the coupling between the electric field and the solid mechanics field is described by the piezoelectric effect, and the coupling between the fluid field and the solid mechanics field is described by the Biot effect; boundary constraint conditions are constructed according to the working state of the low-frequency dynamic elastic modulus measurement device, mainly including fixed displacement (simulating axial pressure), fixed constraint (simulating a rigid base), a boundary pressure gradient of 0 (simulating an undrained boundary condition), and a specified voltage for the piezoelectric material (simulating the power supply condition of the vibration source); in addition to various boundary constraint conditions, contact boundary conditions are used at each contact position; the boundary constraint conditions and contact boundary conditions are loaded into the digital model of the low-frequency dynamic elastic modulus measurement device with the control equations already loaded, and using the finite element method, the digital model of the low-frequency dynamic elastic modulus measurement device with the control equations, boundary constraint conditions, and contact boundary conditions loaded is solved to obtain the resonance frequency and the corresponding modes of the low-frequency dynamic elastic modulus measurement device. Through the present invention, the resonance frequency of the low-frequency dynamic elastic modulus measurement device can be accurately predicted, thereby effectively judging the specific influence of the resonance frequency on the measurement results, providing a theoretical basis for constructing a wider frequency band of low-frequency elastic modulus test states, providing strong technical support for manufacturing more efficient and stable low-frequency dynamic elastic modulus measurement devices, optimizing the device design scheme, and greatly saving the device construction cost.

[0061] In one example, the multi-field constraints include:

[0062] Solid mechanics field constraint, fluid field constraint, and electric field constraint.

[0063] Specifically, the physical laws of the low-frequency dynamic elastic modulus measurement device are mainly constrained by three fields, namely the solid mechanics field, the fluid field, and the electric field. Among them, the solid mechanics field uses the stress-strain equation to constrain the solid materials in the device, including but not limited to the device's skeleton, rock samples, etc.; the fluid field mainly uses Darcy's law to describe the fluid flow inside the rock samples, and the electric field mainly uses the electrostatic field to describe the working electric field of the piezoelectric ceramic vibration source; the coupling between the electric field and the solid mechanics field is described by the piezoelectric effect, and the coupling between the fluid field and the solid mechanics field is described by the Biot effect.

[0064] In one example, the control equations include:

[0065]

[0066]

[0067]

[0068]

[0069] where λ is the eigenvalue, u s is the displacement of the solid material, is the divergence operator, σ is the stress tensor, e is the coupling tensor matrix of the piezoelectric material, ε 0 is the air dielectric coefficient, ε r is the relative dielectric parameter of the material, V is the voltage, C is the stiffness tensor, : is the tensor operator, is the gradient operator, α is the Biot coefficient, ρ b is the density of the solid material, ρ f is the fluid density, ρ c (λ) is the coupling density, specifically τ is the pore tortuosity of the rock sample, φ is the porosity, η is the fluid viscosity, and κ is the permeability of the rock sample.

[0070] In one example, the loading boundary constraint conditions include:

[0071] Apply a specified displacement to the top of the piston column under axial compression: u s = disp;

[0072] Apply a fixed boundary condition to the bottom of the device, and set the displacement u s = 0;

[0073] Apply an undrained boundary condition to the boundary of the rock sample:

[0074] Apply voltages to both ends of the piezoelectric vibration source, which are respectively: V = 0, V = V0;

[0075] Wherein, disp refers to the displacement, and V0 is the high-voltage value.

[0076] In one example, the contact boundary conditions include:

[0077]

[0078] Wherein, p n is the penalty function factor, g n is the contact gap function, p 0 is the pressure when the contact gap is 0, If means if, else means or, that is, if g n <p 0 / p n , then T n =-p n g n +p 0 , otherwise, T n =0.

[0079] Specifically, the contact boundary conditions at each contact position of the digital model of the low-frequency dynamic elastic modulus measurement device of the present invention are The present invention can calculate the frictional force on the contact boundary according to the contact pressure T n as μT n , where μ is the friction coefficient.

[0080] In one example, the finite element method is used to solve the digital model with the loaded boundary constraints and contact boundary conditions.

[0081] In one example, the finite element method includes:

[0082] ARPACK.

[0083] Specifically, the finite element method adopted by the present invention is ARPACK, and other finite element methods can also be used to solve the digital model of the present invention, such as the Subspace method, the Block Lanczos method, the Power Dynamics method, the Reduced Householder method, the Damped method, and the Unsysmmetric method.

[0084] The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments, but it is not a limitation of the present invention. It should be noted that, without conflict, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0085] Embodiment 1

[0086] This embodiment provides a resonance frequency prediction method, including:

[0087] Construct a digital model of the low-frequency dynamic elastic modulus measurement device according to the device size and materials, where the material parameters in the model are obtained through experimental measurements and the device's factory parameters; according to the characteristics of the digital model of the low-frequency dynamic elastic modulus measurement device, use hexahedral meshes, tetrahedral meshes, and scanning and other technical methods to perform mesh division on the digital model of the low-frequency dynamic elastic modulus measurement device; establish a control equation based on multi-field constraints. The physical laws of the low-frequency dynamic elastic modulus measurement device are mainly constrained by three fields, namely the solid mechanics field constraint, the fluid field, and the electric field. Among them, the solid mechanics field uses the stress-strain equation to constrain the solid materials in the device, including but not limited to the device's skeleton, rock samples, etc.; the fluid field mainly uses Darcy's law to describe the fluid flow inside the rock sample, and the electric field mainly uses the electrostatic field to describe the working electric field of the piezoelectric ceramic vibration source; the coupling between the electric field and the solid mechanics field is described by the piezoelectric effect, and the coupling between the fluid field and the solid mechanics field is described by the Biot effect; the control equation includes:

[0088]

[0089]

[0090]

[0091]

[0092] where λ is the eigenvalue, u s is the displacement of the solid material, is the divergence operator, σ is the stress tensor, e is the coupling tensor matrix of the piezoelectric material, ε 0 is the air dielectric coefficient, ε r is the relative dielectric parameter of the material, V is the voltage, C is the stiffness tensor, : is the tensor operator, is the gradient operator, α is the Biot coefficient, ρ b is the density of the solid material, ρ f is the fluid density, ρ c (λ) is the coupling density, specifically τ is the pore tortuosity of the rock sample, φ is the porosity, η is the fluid viscosity, κ is the permeability of the rock sample;

[0093] Construct boundary constraint conditions according to the working state of the low-frequency dynamic elastic modulus measurement device, mainly including fixed displacement (simulating axial pressure), fixed constraint (simulating a rigid base), boundary pressure gradient of 0 (simulating an undrained boundary condition), and specifying the voltage of the piezoelectric material (simulating the power supply condition of the vibration source); apply a specified displacement to the top of the piston column of the axial pressure: u s = disp; apply a fixed boundary condition to the bottom of the device, and set the displacement u s= 0; The boundary condition of the rock sample is loaded with an undrained boundary condition: The voltages applied to both ends of the piezoelectric source are: V = 0, V = V0; where disp refers to the displacement, and V0 is the high voltage value; in addition to various boundary constraint conditions, contact boundary conditions are used at each contact position: Among them, p n is the penalty function factor, g n is the contact gap function, p 0 is the pressure when the contact gap is 0, If means if, else means or, that is, if g n <p 0 / p n , then T n = -p n g n + p 0 , otherwise, T n = 0; According to the contact pressure T n calculate the frictional force on the contact boundary, which is μT n , where μ is the friction coefficient; use the ARPACK finite element method to solve the digital model of the low-frequency dynamic elastic modulus measurement device with the loaded control equations, boundary constraint conditions, and contact boundary conditions to obtain the resonance frequency and the corresponding mode.

[0094] Example 2

[0095] This example provides a resonance frequency prediction method, including:

[0096] Construct a geometric model. According to the device dimensions, construct a digital model of the low-frequency dynamic elastic modulus measurement device. The model is as shown in Figure 2 a; The material parameters in the model are obtained through experimental measurements and factory parameters. The physical properties of solid materials are shown in Table 1, and piezoelectric materials are shown in Table 2;

[0097] Mesh generation. According to the characteristics of the geometric model, use hexahedral meshes, tetrahedral meshes, and scanning and other technical methods to mesh the device geometric model. The meshed model is as shown in Figure 2 b;

[0098] Control equations. The physical laws of the low-frequency dynamic elastic modulus measurement device are mainly constrained by three fields, namely the solid mechanics field constraint, the fluid field, and the electric field. Among them, the solid mechanics field uses the stress-strain equation to constrain the solid materials in the device, including but not limited to the framework of the device, rock samples, etc. The fluid field mainly uses Darcy's law to describe the fluid flow inside the rock sample, and the electric field mainly uses the electrostatic field to describe the working electric field of the piezoelectric ceramic source. The coupling between the electric field and the solid mechanics field is described by the piezoelectric effect, and the coupling between the fluid field and the solid mechanics field is described by the Biot effect;

[0099]

[0100]

[0101]

[0102]

[0103] where λ is the eigenvalue, λ = -iω, and u s is the displacement of the solid material, is the divergence operator, σ is the stress tensor, e is the coupling tensor matrix of the piezoelectric material, and ε 0 is the air dielectric constant, and ε r is the relative dielectric parameter of the material, V is the voltage, C is the stiffness tensor, and : is the tensor operator, is the gradient operator. α is the Biot coefficient. The density of the solid material is ρ b , and ρ f is the fluid density, and the coupled density ρ c (λ) is closely related to the pore tortuosity τ of the rock sample, the fluid density ρ f , the porosity φ, the fluid viscosity η, and the permeability κ of the rock sample, specifically M is the bulk modulus of the pores, and its specific definition is where φ is the porosity and K f is the bulk modulus of the fluid;

[0104] Load the governing equation (1) into all materials including the steel fixture material ( Figure 2 b dark gray area), the aluminum material ( Figure 2 b silver gray area), the sandstone sample ( Figure 2 b light yellow area), the piezoelectric material ( Figure 2 b red area), and the butyl rubber ( Figure 2 b green area), load the governing equation (2) into the piezoelectric ceramic area ( Figure 2 b red area), load the governing equation (3) into the sample area ( Figure 2 b light yellow area), and constrain the fluid flow; the modes of the digital model of the low-frequency dynamic elastic modulus measurement device at non-resonant frequencies are as shown in Figure 2 c, and the modes of the digital model of the low-frequency dynamic elastic modulus measurement device at resonant frequencies are as shown in Figure 2 d;

[0105] Boundary conditions: Based on the working state of the low-frequency dynamic elastic modulus measurement device, boundary constraint conditions are constructed, mainly including fixed displacement (simulating axial compression), fixed constraint (simulating a rigid base), a boundary pressure gradient of 0 (simulating an undrained boundary condition), and a specified voltage for the piezoelectric material (simulating the power supply condition of the seismic source). Specifically:

[0106] Apply a specified displacement (prestress boundary condition) to the top of the piston column under axial compression: u s = disp, where disp is the specified displacement;

[0107] Apply a fixed boundary condition to the bottom of the device, making the displacement u s = 0;

[0108] Set the boundary condition of the rock sample to an undrained boundary condition:

[0109] The voltages at both ends of the piezoelectric seismic source are: V = 0; V = V0, where V0 is the high voltage value.

[0110] In addition to various constrained boundary conditions, contact boundary conditions are used at each contact position,

[0111]

[0112] where p n is the penalty function factor, g n is the contact gap function, and p 0 is the pressure when the contact gap is 0. If means if, that is, if g n < p 0 / p n , T n = - p n g n + p 0 , otherwise, T n = 0.

[0113] Calculate the frictional force on the contact surface according to the contact pressure, which is μT n , where μ is the friction coefficient;

[0114] The contact positions of each material are contact boundary conditions, the friction coefficient is 0.1, the fixed boundary condition is loaded to the bottom of the gray steel base, the specified displacement is loaded to the top of the steel material, the specified displacement is 1.5 microns, the fluid boundary condition is loaded around the yellow sample, the voltage V0 = 600V is loaded to the bottom of the piezoelectric seismic source, the top is grounded, V = 0V;

[0115] Solve: Using the finite element method, solve for the resonance frequency and the corresponding mode according to the control equation and boundary constraint conditions, and evaluate the variation of the device displacement with frequency. The solution method is ARPACK; Figure 3The solid displacement of the model center point position is shown as a function of frequency. It can be seen that the resonance frequency of the device in this embodiment is above 1000 Hz. To verify the accuracy of the evaluation, the device in this embodiment was used to measure tempered glass, and the Young's modulus of the tempered glass is as Figure 4 shown in a, the attenuation of the Young's modulus is as Figure 4 shown in b, the Poisson's ratio of the Young's modulus is as Figure 4 shown in c, and the Poisson's phase difference of the Young's modulus is as Figure 4 shown in d. It can be seen that when the frequency is around 1000 Hz, unstable test results appear for all four elastic parameters, which indicates that the prediction result of the present invention matches the actual test result, proving the effectiveness of the present invention.

[0116] Table 1 Solid Material Parameters

[0117]

[0118] Table 2 Piezoelectric Material Parameters

[0119]

[0120] Embodiment 3

[0121] As Figure 5 shown, this embodiment provides a resonance frequency prediction device, including:

[0122] A construction module for constructing a digital model of a low-frequency dynamic elastic modulus measurement device;

[0123] A mesh generation module for performing mesh generation on the digital model;

[0124] An establishment module for establishing a control equation based on multi-field constraints;

[0125] A first loading module for loading the control equation into the digital model that has completed mesh generation;

[0126] A second loading module for loading boundary constraint conditions and contact boundary conditions into the digital model that has been loaded with the control equation;

[0127] A solution module for solving the digital model that has been loaded with boundary constraint conditions and contact boundary conditions to obtain the resonance frequency and corresponding resonance mode of the digital model.

[0128] Embodiment 4

[0129] This embodiment provides an electronic device, which includes:

[0130] At least one processor; and,

[0131] A memory communicatively connected to the at least one processor; wherein,

[0132] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the resonance frequency prediction method in Embodiment 1.

[0133] An electronic device according to an embodiment of the present disclosure includes a memory and a processor, and the memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0134] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In an embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0135] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain good user experience effects, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included in the protection scope of the present disclosure.

[0136] For a detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details will not be repeated here.

[0137] Embodiment 5

[0138] This embodiment provides a non-transitory computer-readable storage medium, and the non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the resonance frequency prediction method in Embodiment 1.

[0139] A computer-readable storage medium according to an embodiment of the present disclosure stores non-transitory computer-readable instructions thereon. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the methods of the foregoing embodiments of the present disclosure are executed.

[0140] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (such as CD-ROM and DVD), magneto-optical storage media (such as MO), magnetic storage media (such as magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (such as memory card), and media with built-in ROM (such as ROM cartridge).

[0141] The embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

[0142] The embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A resonance frequency prediction method, characterized in that, it includes: Construct a digital model of a low-frequency dynamic elastic modulus measurement device; Perform mesh division on the digital model; Establish control equations based on multi-field constraints; Load the control equations into the digital model that has completed the mesh division; Load boundary constraint conditions and contact boundary conditions into the digital model that has been loaded with the control equations; Solve the digital model loaded with the boundary constraint conditions and the contact boundary conditions to obtain the resonance frequency and the corresponding resonance mode of the digital model.

2. The resonance frequency prediction method according to claim 1, characterized in that, the multi-field constraints include: Solid mechanics field constraint, fluid field constraint and electric field constraint.

3. The resonance frequency prediction method according to claim 2, characterized in that, the control equations include: where λ is the eigenvalue, u s is the displacement of the solid material, is the divergence operator, σ is the stress tensor, e is the coupling tensor matrix of the piezoelectric material, ε 0 is the air dielectric constant, ε r is the relative dielectric parameter of the material, V is the voltage, C is the stiffness tensor, : is the tensor operator, is the gradient operator, α is the Biot coefficient, ρ b is the density of the solid material, ρ f is the fluid density, ρ c ρ(λ) is the coupling density, specifically τ is the pore tortuosity of the rock sample, φ is the porosity, η is the fluid viscosity, κ is the permeability of the rock sample.

4. The resonance frequency prediction method according to claim 1, characterized in that, loading the boundary constraint conditions includes: Apply a specified displacement to the top of the axially compressed piston column: u s = disp; Fixed boundary conditions are applied at the bottom of the device, setting the displacement u s = 0; Loading of boundary conditions for rock samples - Undrained boundary conditions: The voltages loaded at both ends of the piezoelectric vibration source are respectively: V = 0, V = V0; where disp refers to the displacement, and V0 is the high voltage value.

5. The resonance frequency prediction method according to claim 1, characterized in that, the contact boundary conditions include: Among them, p n is the penalty function factor, g n is the contact clearance function, p 0 is the pressure when the contact clearance is 0. If means if, else means or, that is, if g n < p 0 / p n , then T n = - p n g n + p 0 , otherwise, T n = 0.

6. The resonance frequency prediction method according to claim 1, characterized in that, The finite element method is used to solve the digital model loaded with the boundary constraint conditions and the contact boundary conditions.

7. The resonance frequency prediction method according to claim 6, characterized in that, the finite element method includes: ARPACK.

8. An electronic device, characterized in that, the electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the resonance frequency prediction method according to any one of claims 1-7.

9. A non-transitory computer-readable storage medium, characterized in that, This non-transitory computer-readable storage medium stores computer instructions, and these computer instructions are used to cause a computer to execute the resonance frequency prediction method according to any one of claims 1-7.

10. A resonance frequency prediction device, characterized in that, it includes: A construction module for constructing a digital model of a low-frequency dynamic elastic modulus measurement device; A mesh division module for performing mesh division on the digital model; An establishment module for establishing control equations based on multi-field constraints; A first loading module for loading the control equations into the digital model that has completed the mesh division; A second loading module for loading boundary constraint conditions and contact boundary conditions into the digital model that has been loaded with the control equations; A solving module for solving the digital model loaded with the boundary constraint conditions and the contact boundary conditions to obtain the resonance frequency and the corresponding resonance mode of the digital model.

Images