63 results about "Solid mechanics" patented technology

The invention discloses a fluid-solid coupling numerical simulation method and device based on an embedded discrete fracture model. A structured grid is adopted to divide a reservoir matrix, and a discrete fracture model is embedded to construct an oil and gas reservoir geological grid. Establishing a transmission function, and formulating a porous medium fluid seepage control equation according to the transmission function. The equation is solved at each time node, a pore pressure field is obtained, and an equivalent pore pressure field is obtained through weighted averaging and serves as a boundary condition of the solid mechanical control equation. Expanding an oil and gas reservoir geological grid model, combining a crack-rock constitutive model, an effective stress principle and a stress balance equation, constructing a solid mechanical control equation, and solving and simulating rock deformation by adopting a finite element method and a Newton-Rafson method. Based on the constitutive relation of the composite material, matrix and crack seepage attributes are updated, the interaction between a seepage field and a stress field is considered in an iterative coupling mode, and multi-field coupling simulation is achieved. The problem that an existing method cannot be applied to a real oil and gas reservoir is solved.

The invention provides a rock mass high-temperature deformation simulation method and system based on deep learning near-field dynamics, and relates to the field of rock mass thermal-mechanical coupling analysis and research, and the method comprises the steps: obtaining the parameter information of a to-be-simulated target rock mass, including the geometric shape of the rock mass and the physical attribute of the rock mass; dispersing the target rock mass into a plurality of material points according to the parameter information, and constructing a thermal-mechanical coupling model of the target rock mass; based on the initialized thermal-mechanical coupling model, time step iterative simulation of a solid mechanical field is conducted, and finally the crack propagation process and temperature distribution of the target rock mass are obtained; according to the method, a solid field near-field dynamics theory based on strain energy density is introduced, a fluid field near-field dynamics equation is combined, temperature distribution in a rock medium is predicted through a physical information neural network, a deep learning near-field dynamics framework of rock mass high-temperature deformation damage is formed, and therefore rock mass high-temperature deformation is efficiently and accurately simulated.

The invention discloses a 2.5 D Chiplet structure physical field coupling characteristic analysis method. The method specifically comprises the following steps: step 1, adding a solid heat transfer physical field interface, a solid mechanical physical field interface and a thermal expansion physical field interface to a 2.5 D Chiplet structure; 2, boundary conditions are added to a solid heat transfer physical field interface, and the temperature distribution condition of the 2.5 D Chiplet structure is obtained; step 3, solving the deformation and stress distribution conditions of the 2.5 D Chiplet structure; 4, the solid heat transfer physical field interface and the solid mechanical physical field interface are coupled through the expansion physical field interface; and 5, carrying out grid division on the 2.5 D Chiplet structure to obtain the temperature, stress and deformation distribution conditions of the structure. The method can accurately predict the distribution of the temperature field and the force field on the packaging chip, and facilitates the thermal analysis of the packaging structure.

The invention discloses a method and a system for simulating influence of interface pressure of a prefabricated cable intermediate joint on ultrasonic non-linear parameters. The simulation method comprises the following steps: establishing a two-dimensional simulation solid mechanics geometric model; setting material parameters of the solid mechanical geometric model; setting boundary conditions of the solid mechanical geometric model; a solid mechanical geometric model and setting of material parameters and boundary conditions are utilized, and a steady-state mechanical analysis mode is adopted to conduct contact stress simulation under different stress conditions; importing the solid mechanical geometric model subjected to stress simulation into a sound field coupling model, coupling a solid mechanical physical field and an electrostatic physical field in a sound field, and performing ultrasonic signal simulation; and in combination with the actual state of the prefabricated cable intermediate joint, carrying out simulation calculation on the influence of different interface pressures on ultrasonic non-linear parameters. The accurate value of the interfacial pressure can be obtained by analyzing the nonlinear parameters of the echo signals.

The invention belongs to the technical field of laser and substance interaction, and relates to a simulation method and system for a temperature field and a stress field of a short-pulse laser irradiation material. The simulation method comprises the following steps: determining simulation parameters including short pulse laser irradiation parameters and material attribute parameters; constructing a thermal-mechanical coupling transient model of the material, and determining boundary conditions of a heat transfer physical field and a solid mechanics physical field of the model; meshing division is carried out on the thermal coupling transient model, and solver parameters are configured; the method specifically comprises the following steps: inputting a short pulse laser pulse width and a short pulse laser frequency into a pre-trained linear neural network model, and calculating to obtain solver parameters; and solving the thermal coupling transient model by applying a solver to obtain temperature field and stress field distribution conditions of the material at different irradiation moments. According to the simulation method, the complex transient physical process of the material under the action of the short pulse laser can be comprehensively and accurately captured.

The invention discloses a simulation method and system for progressive damage and failure removal of a composite material, and belongs to the field of computer-aided engineering, and the method comprises the following steps: S1, coupling solving: solving a partial differential control equation set describing heat transfer and solid mechanics, and considering a coupling effect caused by thermal strain and the influence of temperature on material attributes, obtaining the temperature field and displacement field distribution of the current time step or iteration step; s2, material state calculation: on the basis of the temperature and the displacement field obtained in the coupling solution step, calculating the material state of an integral point of each unit in the model by using a predefined composite material constitutive model; and S3, unit processing and grid updating: after the material state calculation is completed, executing unit processing and grid updating. According to the invention, the final failure behavior of the structure can be predicted more comprehensively and accurately.

The invention discloses a horizontal salt cavity circulating hydrogen storage and release whole process representation numerical simulation method which comprises the following steps: constructing an initial stratum model which comprises a mudstone upper layer, a rock salt interlayer and a mudstone lower layer; an ellipsoid is embedded in the initial stratum model, and a one-fourth stratum model with simplified calculation is obtained; adding a physical field module to the 1 / 4 stratum model to obtain an optimal stratum model; wherein the multi-physics field module comprises a solid mechanics physics field module, a Darcy law physics field module and a porous medium heat transfer physics field module; the optimal stratum model is subjected to grid division, the initial ground stress of the divided optimal stratum model is calculated firstly, then multi-field coupling data after the initial ground stress is added are calculated, and the multi-field coupling data comprise a temperature field, a stress field and a displacement field. The design capacity of the reservoir can be effectively increased, and the stable leakproofness of the reservoir is ensured by improving the structural form of the gas storage.

The invention discloses a free vibration semi-analytical calculation method for a liquid filling functional gradient pipeline on an elastic foundation, and relates to the technical field of solid mechanics and fluid-solid coupling vibration analysis. The method comprises the following steps: acquiring materials, geometry and boundary parameters of a pipeline, a fluid and a foundation, establishing a fluid-solid coupling model based on a three-dimensional elastic theory, simulating foundation support by adopting a two-parameter Pasternak model, and deducing a control equation and boundary conditions of a pipeline structural domain and a fluid domain; discretizing a pipeline reference surface and a fluid boundary through a two-dimensional high-order spectrum unit, and converting an equation into a standard coordinate system in combination with proportional boundary coordinate transformation; proportional boundary finite element control equations of a structural domain and a fluid domain are deduced respectively and converted into a first-order differential equation set, then a fine integration method and an eigenvalue decomposition method are adopted for solving, and then based on a fluid-solid coupling condition coupling equation, the equation is converted into a generalized eigenvalue problem for solving the inherent frequency and the vibration mode. According to the method, only precision, efficiency and discrete boundaries are considered, and support is provided for pipeline design optimization and safety evaluation.

The invention discloses a fracturing optimization method and system for preventing pressure channeling of a deep shale reservoir and a medium, and belongs to the technical field of shale oil and gas geological exploitation. Comprising the following steps: acquiring different-scale fracture characteristics and ground stress distribution rules of a deep shale reservoir; different-scale fracture mechanical parameters and seepage parameter characteristics of the rock sample are obtained; building a deep shale reservoir multi-physics coupling fracturing model by using a damage mechanics theory according to deep shale reservoir matrix characteristics; the method comprises the following steps: setting parameters according to solid mechanics, porous media and Darcy law based on a deep shale reservoir and different-scale fracture characteristics, and establishing a platform fracturing numerical model; on the basis of a platform fracturing numerical model, an optimal fracturing chart influenced by multi-scale fracture development is obtained through orthogonal experimental simulation, so that a fracturing parameter combination enabling the effective fracturing volume to be maximum and avoiding inter-well pressure channeling is determined; according to the fracturing optimization method for preventing pressure channeling of the deep shale reservoir, efficient development of deep shale oil and gas is effectively guaranteed.

The invention provides a molten pool dynamic simulation model and system for simulating laser drilling based on a level-set method, and the model comprises a two-dimensional molten pool simulation model which comprises heat transfer-laminar flow (two-phase flow) and level-set physical field coupling; on the basis of the confirmed implementation process parameters and material physical property parameters, a molten pool interface simulation metallographic phase under given process operation is compared with an actual metallographic phase, and a heat transfer mathematical model in the laser drilling process is obtained through simulation confirmation; and the three-dimensional thermal-mechanical coupling model is constructed based on a heat transfer mathematical model and further coupled with a solid mechanical physical field, so that temperature distribution and stress distribution at the welding printing position in the laser drilling process are obtained through simulation. The technical problems that in an existing laser drilling welding technology, parameter debugging depends on experience, dynamic monitoring of a molten pool is difficult, and thermal stress prediction means are lacked can be effectively solved, and powerful theoretical support and data basis are provided for optimizing the welding technology and improving the welding quality.

The invention provides an LBM large-deformation fluid-solid interface stabilization method and system based on energy flux regulation, and belongs to the technical field of computational fluid mechanics and solid mechanics coupling numerical simulation. Identifying a fluid-solid boundary area, and constructing a buffer layer area; solving fluid in a buffer layer area, and embedding an energy flux regulation mechanism; applying a load to the structure through fluid, and calculating a structure response; fluid is fed back, and the dynamic boundary is updated; and iterating repeatedly, and stopping and outputting after a set termination condition is reached. Active regulation and control of the fluid-structure interface energy exchange process are achieved, the numerical stability and physical precision of fluid-structure interaction simulation under the large deformation condition are remarkably improved, and the method has the advantages of being simple in structure, high in expandability and low in calculation overhead; the method is suitable for complex coupling scenes such as flexible structural body propulsion, airfoil pitching motion and high Reynolds number structural response analysis, and has wide engineering application value.

The present application belongs to the technical field of laser and material interaction, and relates to a kind of optical element laser irradiation damage prediction method, system and terminal.The damage prediction method includes: establishing the three-dimensional geometric model of optical element, defining the material attribute parameter and environmental parameter of optical element, and constructing laser energy deposition function;With laser energy deposition function as laser heat source, a three-dimensional transient model of solid heat transfer and solid mechanics bidirectional coupling is constructed, and the boundary conditions of solid heat transfer and solid mechanics are set;Dynamic damage evolution equation is constructed, the laser irradiation surface of optical element is meshed, the solver is configured based on condition criterion, the current damage variable distribution is obtained by solving, the material attribute parameter is updated based on the current damage variable distribution, the damage variable distribution is continuously calculated after updating, until the preset termination condition is reached, and the damage prediction result is obtained.The present application can comprehensively and accurately predict the complex transient physical process of material under the action of high-repetition-rate short-pulse laser.

The embodiment of the invention provides a method and device for obtaining thermodynamic characteristics of a fan tower based on COMSOL, and relates to the technical field of data processing, and the method specifically comprises the steps: constructing a three-dimensional transient model set for the fan tower based on COMSOL; the three-dimensional transient model set comprises a solid heat transfer model, a solid mechanical model and a surface-to-surface radiation model; acquiring a simulation parameter set; according to the simulation parameter set, constructing a target geometric model corresponding to the wind turbine generator tower; determining an initial value and a boundary condition corresponding to each model in the three-dimensional transient model set; and after grid division is carried out on the target geometric model, solving calculation is carried out on each model in combination with a material attribute parameter set, an environment parameter set and the initial value and the boundary condition corresponding to each model so as to obtain the temperature, the stress and the strain of each node corresponding to each grid. According to the invention, the tower drum under sunlight can be simulated, and thermodynamic characteristics can be obtained.

The present invention belongs to the technical field of the interaction between laser and matter, and relates to a simulation method and system for the temperature field and stress field of a material irradiated by short-pulse laser. The simulation method includes: determining simulation parameters, including short-pulse laser irradiation parameters and material property parameters; constructing a thermo-mechanical coupling transient model of the material, and determining the boundary conditions of the heat transfer physical field and the solid mechanics physical field of the model; performing grid division on the thermo-mechanical coupling transient model, and configuring solver parameters; specifically: inputting the short-pulse laser pulse width and the short-pulse laser frequency into a pre-trained linear neural network model, and calculating to obtain the solver parameters; applying the solver to solve the thermo-mechanical coupling transient model to obtain the temperature field and stress field distribution of the material at different irradiation times. The simulation method of the present invention can comprehensively and accurately capture the complex transient physical process of the material under the action of short-pulse laser.

The application belongs to but is not limited to the technical field of electronic core industry, and discloses a polycrystalline NCM electrode particle cracking and pulverization prediction method and system, which comprises the following steps: constructing an anisotropic polycrystalline fracture phase field model, including a grain boundary phase field equation, a solid mechanics equation, a Li + diffusion equation and a fracture phase field equation; constructing a two-dimensional geometry of a polycrystalline NCM electrode particle, setting required parameters of physical fields, carrying out grid division and setting initial conditions and boundary conditions; using a finite element simulation software to perform numerical solution, and obtaining cracking and pulverization distribution of the polycrystalline NCM electrode particle. + The application considers bidirectional coupling of Li diffusion and mechanical stress, cracking and pulverization and mechanical stress, and uses grain boundary phase field calculation to correct the critical energy release rate, so that the application can accurately reflect the mechanism of cracking and pulverization of the polycrystalline NCM electrode particle, and improve the accuracy of cracking and pulverization prediction.

This invention proposes a method and system for simulating the vibrational displacement distribution of magnetic particles based on multiphysics coupling, comprising: Step S1: constructing a three-dimensional geometric model containing biological tissue and magnetic particles; Step S2: setting material parameters for each region; Step S3: configuring multiphysics; Step S4: performing mesh generation; Step S5: setting solver parameters; Step S6: using the finite element method to solve the force of the magnetic field on the magnetic particles in Maxwell's equations and describing the magnetic field distribution in the magnetic particle region; Step S7: using the magnetic force obtained in Step S6 as a load, inputting it into the solid mechanics field module, and solving the vibrational displacement distribution in the magnetic particle region based on the wave equation. Furthermore, the system of this invention includes: a three-dimensional geometric model group, a material library, a mesh module, and a multiphysics module. This invention accurately simulates the vibrational displacement distribution of magnetic particles within biological tissue under the action of an alternating magnetic field, revealing the coupling law between magnetic fields and solid mechanics.

The invention discloses a radial point interpolation finite element method for linear elastic solid harmonic response and application, and belongs to the technical field of solid mechanics and engineering numerical simulation. The method comprises the following steps: discretizing a linear elastic solid problem domain into a standard finite element grid, deriving a nodal function through a radial basis function, selecting a local enhanced interpolation basis function to construct a newly added degree of freedom and complete local numerical approximation, and constructing a global interpolation format of a displacement vector field in an element in combination with unit decomposition characteristics; based on the interpolation format and a Galerkin weak form of an elastic dynamic problem, applying a displacement boundary condition to obtain a system forced vibration equation; acquiring steady-state dynamic response and frequency domain response characteristics under different frequency harmonic excitation; according to the method, high-order interpolation can be realized without increasing nodes, the numerical precision is improved, the computing resource and time consumption is reduced and a more reliable numerical solution is provided under the condition that the grid density is relatively small, and the method has important application value in engineering practice.

The invention discloses a two-dimensional simulation method for a bulk acoustic wave delay line based on COMSOL Muliphysics. The method comprises the following steps: S1, establishing a two-dimensional geometric model of the bulk acoustic wave delay line by adopting a geometric module in C0MS0L Multiphysics simulation software; s2, performing region division on the geometric model, and setting material parameters; s3, setting parameters of a solid mechanics physical field module and an electrostatic physical field, performing grid division on each region, setting an event physical field, and correspondingly setting the parameters of each physical field module; s4, setting multi-physical field coupling; and S5, calculating and solving to obtain a simulation result of the bulk acoustic wave delay line. The vibration mode and the S parameter are determined through analog simulation of the bulk acoustic wave delay line, and the method can shorten the time cost for actual use and reduce the manufacturing cost.

The application discloses a shield machine main drive sealing performance prediction method and device under actual complex service conditions, and relates to the field of mixed lubrication. The method takes a rotary lip seal as a research object, is based on a two-dimensional axisymmetric finite element model of a shield machine main drive rotary lip seal, and combines solid mechanics, lubricating oil film fluid mechanics, sealing lip surface micro-convex body contact mechanics, sealing lip surface normal and tangential elastic deformation, and thermodynamic analysis to quantitatively analyze data representing sealing performance of the shield machine main drive seal under actual complex service conditions, such as viscous friction and friction torque under different pressure differentials. The application can predict the sealing performance of the shield machine main drive seal under actual complex service conditions.

The invention provides a rock mass chemical corrosion simulation method and system based on neural network near-field dynamics, and relates to the technical field of fractured rock mass chemical corrosion simulation, and the method comprises the steps: building a solid mechanical field model based on a solid field near-field dynamics theory of a Hoek-Brown failure criterion, and building a chemical field model based on a neural network near-field dynamics theory; in a solid field, average principal stress among material points is calculated in real time based on a Hoek-Brown failure criterion, a nonlinear shear tensile failure relation in a stress state is obtained, in a chemical field, a concentration distribution prediction model is incorporated into the chemical field, and in an iteration time step of each solid mechanical field, a nonlinear shear tensile failure relation is obtained. And predicting the concentration distribution of the substance points in real time by adopting a concentration distribution prediction model, calculating a corresponding stress state, updating damage parameters, feeding back the updated damage parameters to the concentration distribution prediction model, and iteratively predicting the concentration distribution of the substance points in the next time step until an iteration termination condition is reached, thereby realizing a stress-corrosion coupling simulation calculation process.

The invention particularly relates to a solid mechanical calculation method and device based on double-grid dynamic weight optimization, and the method comprises the steps: judging whether a solid mechanical calculation demand in a plane plate tension process exists or not; if so, establishing a solid mechanics problem in the surface plate tension process, and establishing a coarse mesh discrete model and a fine mesh discrete model; establishing a first nonlinear mapping relation and a second nonlinear mapping relation based on a preset deep neural network and the two models; and constructing a composite loss function, and based on a preset training strategy, performing staged cooperative training on the composite loss function according to the first nonlinear mapping relationship and the second nonlinear mapping relationship until a preset training round is reached or the trained composite loss function meets a preset precision requirement, and solving according to the trained composite loss function to obtain a solution of a solid mechanics problem displacement field. Therefore, the problem that the convergence speed of a fine mesh model is low in the prior art is solved, and the training efficiency is improved.

The embodiment of the present application provides a method for analyzing the evolution of defects caused by fracture and electromigration based on the phase field method, determining the parameter information of the defective specimen and the multi-physical field load information, and establishing a simplified model; wherein, the multi-physical field load information includes: the load and boundary conditions of the force-thermal-electric field; based on the simplified model, a defect evolution phase field model is constructed; the defect evolution phase field model includes an electromigration phase field model coupled with the force-thermal-electric field and a thermoelastic fracture phase field model; the defect evolution phase field model is solved. In this solution, for each time step in the electromigration phase field evolution process, the electrostatic model, the electrothermal model, the thermoelastic solid mechanics model, and the thermoelastic fracture phase field model are solved for the defect morphology, so as to classify the crack penetrating the defective specimen into the final defect morphology and obtain the complete process of the evolution of the defect morphology over time.

A method for measuring high-strain-rate dynamic hardness of a material surface based on AFM knocking belongs to the technical field of solid mechanical property testing and high-speed machining, and comprises the following steps: S1, enabling a material to generate plastic deformation at a high strain rate by an AFM probe; s2, accurately measuring an applied load and a material deformation amount; and S3, performing post-processing analysis to calculate the dynamic hardness of the material. Compared with an existing method for testing the mechanical property of the material under the high strain rate, the method has the advantages that the strain rate of plastic deformation is increased on the basis of ensuring nanoscale deformation of the material, and the dynamic mechanical property response condition of the surface of the material is accurately measured. Therefore, the test method provided by the invention promotes the development of the high-strain-rate mechanical property of the surface of the test material under the nanoscale.

A simulation analysis device and analysis method for powder pressing, belonging to the field of powder metallurgy. To solve the problem of inaccurate numerical model of powder pressing. The present invention applies Comsol Multiphysics to establish a two-dimensional axisymmetric space dimension, and selects solid mechanics to establish a geometric model of powder pressing, respectively setting a definition part, a solid mechanics 1 part, and a solid mechanics 2 part; setting the material properties of the metal powder part and the mold part, and meshing the geometric model; setting a study-two-stage pressing, including setting and modifying the model configuration of the research step, and auxiliary scanning; setting a study-single-stage pressing, including establishing a research setting, modifying the model configuration of the research step, and auxiliary scanning; performing solution calculations to obtain simulation results of powder pressing, and processing each of the simulation results. The present invention provides theoretical guidance and technical support for the application of powder pressing.

The invention provides a material deformation behavior high-robustness prediction method and system based on a bionic intelligent constitutive structure, and belongs to the technical field of solid mechanics and computational intelligence cross, and the method comprises the steps: obtaining a multi-modal data sample under multiple materials and multiple scenes; the method comprises the following steps: constructing an initial BP neural network constitutive model, taking a multi-modal data sample as input, taking stress as output, simulating a stem cell hierarchical differentiation mechanism by using a bionic hierarchical architecture to obtain an optimized genetic algorithm, optimizing an initial weight and an initial threshold of the initial BP neural network constitutive model by using the optimized genetic algorithm to obtain a transition weight and a transition threshold, and constructing the stem cell hierarchical differentiation mechanism based on the transition weight and the transition threshold. Optimizing the transition weight and the transition threshold through a gradient descent algorithm to obtain an improved BP neural network constitutive model; and predicting the deformation behavior of the to-be-predicted material based on the improved BP neural network constitutive model. According to the method, the deformation behaviors of multiple materials in multiple scenes can be quickly and accurately predicted.

The application discloses a kind of deep shale reservoirs anti-channeling fracturing optimization method, system and medium, belong to shale oil and gas geological exploitation technical field;Including the following steps: obtaining the different scale fracture characteristics and ground stress distribution law of deep shale reservoir;Different scale fracture mechanics parameters and percolation parameter characteristics of rock sample are obtained;By using damage mechanics theory to construct deep shale reservoir multi-physical field coupling fracturing model of matrix characteristics of deep shale reservoir;Based on deep shale reservoir and different scale fracture characteristics, according to solid mechanics, porous medium and Darcy's law, the parameter is set, and the platform pressure numerical model is established;Based on platform pressure numerical model, the optimal fracturing chart affected by multi-scale fracture development is obtained by orthogonal experiment simulation, to determine the maximum effective fracturing volume and avoid interwell channeling of fracturing parameter combination;The application proposes a kind of deep shale reservoir anti-channeling fracturing optimization method, effectively ensures the efficient development of deep shale oil and gas.

This invention discloses a method and system for simulating progressive damage and failure removal in composite materials, belonging to the field of computer-aided engineering. The method includes the following steps: S1, coupled solution: solving the partial differential governing equations describing heat transfer and solid mechanics, considering the coupling effects caused by thermal strain and the influence of temperature on material properties to obtain the temperature field and displacement field distribution of the current time step or iteration step; S2, material state calculation: based on the temperature and displacement fields obtained in the coupled solution step, calculating the material state for each integral point of the model using a predefined composite material constitutive model; S3, element processing and mesh updating: performing element processing and mesh updating after the material state calculation is completed. This invention can more comprehensively and accurately predict the final failure behavior of the structure.

This disclosure provides a method and apparatus for iterative coupled hybrid experiments on fluid-structure interaction boundaries, applicable to fluid mechanics, solid mechanics, numerical simulation, and engineering fields. The experimental method includes: decomposing an initial fluid-structure interaction structure into multiple substructures, including physical and numerical substructures; determining the influencing factors of the experimental convergence index based on the response error between the physical and numerical substructures, where the response error characterizes the error between the input response information of the physical substructure and the output response information of the numerical substructure, and the influencing factors include the iterative algorithm and the substructures; updating the input response information, output response information, and substructures based on an algorithm optimization strategy and parameter update strategy using the iterative algorithm when the response error does not meet a preset error threshold, thereby obtaining the target structure; and determining the hydrodynamic pressure time history and target motion results based on the target numerical substructure, the target physical substructure, and the preset response flow.

The present application discloses a fluid-solid coupling numerical simulation method and device based on an embedded discrete fracture model. The reservoir matrix is ​​divided by structured grids, and a discrete fracture model is embedded to construct a geological grid for oil and gas reservoirs. A transfer function is established, and the porous media fluid seepage control equation is formulated accordingly. The equation is solved at each time node to obtain the pore pressure field, and the equivalent pore pressure field is obtained by weighted average as the boundary condition of the solid mechanics control equation. The geological grid model of the oil and gas reservoir is expanded, and the solid mechanics control equation is constructed by combining the fracture-rock constitutive model, the effective stress principle and the stress equilibrium equation. The finite element method and the Newton-Raphson method are used to solve and simulate rock deformation. Based on the constitutive relationship of the composite material, the matrix and fracture seepage properties are updated, and the interaction between the seepage field and the stress field is considered through iterative coupling to achieve multi-field coupling simulation. The problem that the existing methods cannot be applied to real oil and gas reservoirs is solved.

The invention discloses an acoustic modeling method based on a Voronoi layered double-phase microstructure, which comprises the following steps: generating a Voronoi crystal grain structure according to the average crystal grain size of a polycrystalline metal material; the method comprises the following steps: constructing a layered double-phase microstructure geometric model by setting interlayer spacing and an area proportion coefficient; different elastic matrixes are given to the grain matrix and the layered second-phase microstructure through a Bond transformation method; according to the elastic matrix, a solid mechanics physical field, an ultrasonic excitation source and boundary conditions are set through joint simulation software, grid division is conducted, and an acoustic simulation model used for simulating ultrasonic propagation is established; and according to the acoustic simulation model, setting a solver for calculation so as to obtain and output the propagation characteristics of the ultrasonic waves in the layered double-phase microstructure. The method disclosed by the invention is of great significance to research on the propagation rule of the ultrasonic waves in the polycrystalline metal material with the layered double-phase microstructure.