Valve structure simulation method and device based on finite element analysis and readable storage medium

Abstract

The invention discloses a valve structure simulation method and device based on finite element analysis and a readable storage medium, and belongs to the technical field of valves, the simulation method mainly comprises a multi-physics field coupling method, control equations of a plurality of physics fields are simultaneous to obtain a control equation set, and the control equation set comprises a physics field parameter array corresponding to the physics field; the control equation set comprises a fluid mechanics equation, a damp-heat conduction equation and an elastic mechanics equation, a unified equation set is formed, and then discretization processing is carried out; performing discretization to convert a continuous partial differential equation into an algebraic equation set; for a transient problem, a display or implicit time integration scheme is adopted to perform time discretization on a simultaneous equation. An apparatus includes a memory, a processor, and a computer program. A terminal program for executing the method is stored in the storage medium. Factors such as fluid flow, damp heat conduction and structural mechanics can be considered at the same time, the simulation real-time performance is improved by integrating an advanced method and an optimization strategy, and the precision is improved while the working efficiency is ensured.

Description

Technical Field

[0001] The present application belongs to the field of valve technology, and in particular relates to a valve structure simulation method and device based on finite element analysis, and a readable storage medium. Background Art

[0002] Valves are control components in fluid delivery systems, with functions such as shutoff, regulation, diversion, backflow prevention, pressure stabilization, diversion, or overflow pressure relief. Valves are devices used to control the direction, pressure, and flow of fluids in fluid systems. They are devices that allow the flow of media (liquids, gases, powders) within piping and equipment to flow or stop and control their flow. Valves achieve functions such as diversion, shutoff, throttling, non-return, diversion, or overflow pressure relief by changing the cross-section of the passage and the direction of the medium flow. An electric valve is a device that uses an electric actuator to drive the valve for switching or regulation. The electric valve converts electrical energy into mechanical energy and drives the valve to open or close through a motor, thereby achieving fluid control.

[0003] The valve industry is developing toward intelligent, automated, high-precision, and high-efficiency features. For example, servo valves are high-precision control valves. As a key component of a servo system's precision hydraulic or pneumatic control, servo valves are used to precisely control fluid pressure and flow. They offer advantages such as fast response, high output power, and high control accuracy. Their performance directly impacts the reliability and safety of servo system operations, and they are widely used in fields such as petroleum, chemical engineering, metallurgy, shipbuilding, aviation, and aerospace. Because the operating principle of servo valves involves the interaction of multiple physical fields (such as fluid mechanics, hygrothermal dynamics, and structural mechanics), multi-physics simulation is crucial for optimizing servo valve design and performance.

[0004] In modern engineering design of servo valves, many systems require consideration of the interactions between multiple physical phenomena, such as fluid flow, heat transfer, and structural mechanics. Traditional single-physics simulations cannot accurately capture these complex interactions, resulting in inaccurate or inefficient design results. While existing multi-physics coupling simulation methods exist, the simplified assumptions currently employed can lead to simulation results that are inconsistent with actual conditions. Other approaches comprehensively consider multiple physics fields, but often only consider fluid-structure coupling or thermo-structure coupling. This inaccurate simulation results leads to insufficient precision in valve structure design, limiting improvements in accuracy, especially in complex engineering applications such as aerospace and intelligent equipment, which require high precision. Furthermore, current valve structure design suffers from the following drawbacks: high computational cost, as multi-physics coupling typically requires extensive computing resources, making it difficult to meet these requirements through manual calculations or experimental testing; difficulty in convergence, as strong coupling between different physical fields can lead to numerical instability and increased solution difficulty; and low computational efficiency, as extensive simulation calculations are performed on the valve structure during design, a task currently limited in efficiency. Therefore, developing a simulation method that can effectively handle multi-physics coupling is crucial for improving design accuracy and efficiency. Summary of the Invention

[0005] The present application aims to at least to some extent solve the technical problem of inaccurate simulation results of valve structures. To this end, the present application provides a valve structure simulation method and device based on finite element analysis, and a readable storage medium, which can simultaneously consider factors affecting the valve structure such as fluid flow, heat and moisture conduction, and structural mechanics. By establishing a multi-physical field coupling relationship for simulation, advanced methods and optimization strategies are integrated to improve the real-time performance of multi-physical field coupling simulation, thereby ensuring work efficiency while improving accuracy.

[0006] In a first aspect, an embodiment of the present application provides a valve structure simulation method based on finite element analysis, which includes:

[0007] Designing a geometric model according to structural parameters of the valve structure and dividing the geometric model into multiple components;

[0008] Define the material parameters of each component and associate the material parameters with each component;

[0009] Finite element simulation model software is used to mesh the geometric model to form mesh units, and all mesh units are used as simulation analysis objects;

[0010] Establish a coupled model of the effects of fluid flow, moisture transfer, heat transfer, and mechanical deformation on valve structure; set boundary conditions and initial conditions;

[0011] Import the coupling model into the simulation model software to perform simulation analysis of multi-physics field coupling;

[0012] Based on the simulation results, mesh elements with failure risks in the geometric model are identified and the design is optimized.

[0013] In some embodiments, the coupling model includes:

[0014] The continuity equation of fluid is used to describe the motion state of fluid flow;

[0015] The heat conduction equation is used to describe the temperature changes of each component caused by the ambient heat source and the component's own Joule heating;

[0016] The moisture diffusion equation is used to describe the humidity changes of each component caused by ambient humidity and its own moisture diffusion;

[0017] Dynamic equations are used to describe the changes in displacement, velocity, acceleration and stress of each component of the valve structure itself.

[0018] In some embodiments, the dynamic equation is described based on a multi-degree-of-freedom system as the relationship between the displacement vector, velocity vector, acceleration vector and external force vector, wherein the coefficients of the displacement vector, velocity vector and acceleration vector are the mass matrix, damping matrix and stiffness matrix respectively.

[0019] In some embodiments, when performing simulation analysis of multi-physical field coupling, the equations corresponding to the coupling model are solved to obtain simulation results of temperature distribution, humidity distribution, stress distribution, and displacement distribution.

[0020] In some embodiments, when solving the equations corresponding to the coupled model:

[0021] Introducing the energy conservation equation, mass conservation equation, momentum conservation equation, and coupling mechanism of multi-physics field coupling;

[0022] Discretize the continuous multi-physics field into a finite number of physics field units corresponding to the grid units;

[0023] The discretized algebraic equations for each physics element are solved to obtain an approximate solution.

[0024] In some embodiments, introducing a coupling mechanism includes: establishing a relationship between temperature, humidity and stress through an equilibrium equation, a displacement equation and a thermal stress-strain equation; establishing a relationship between stress changes caused by temperature and stress changes caused by humidity through thermal stress and wet stress; establishing a relationship for introducing a phase change effect by adding a latent heat of vaporization term; and establishing a relationship between flow fluid pressure, density and flow velocity through a continuity equation and a mass conservation equation.

[0025] In some embodiments, solving the discretized algebraic equations specifically includes:

[0026] By introducing test functions, the partial differential equations of energy conservation equation, mass conservation equation and momentum conservation equation are transformed into integral forms to form weak form expressions respectively;

[0027] Substitute the temperature field corresponding to each grid cell expressed by the interpolation function into the weak form corresponding to the energy conservation equation, substitute the humidity field into the weak form corresponding to the mass conservation equation, and substitute the displacement field into the weak form corresponding to the momentum conservation equation to obtain the discretized equations respectively;

[0028] Use interpolation functions to represent the temperature, humidity, and displacement of each grid cell;

[0029] The discretized equations are solved by coupled simulation.

[0030] In some embodiments, when performing coupled simulation to solve each discretized equation:

[0031] The momentum conservation equation is discretized using the Newmark method, and then coupled with the discretized energy conservation equation and the discretized mass conservation equation for solution.

[0032] In some embodiments, when performing coupled simulation solutions, a coupling mechanism is introduced, including: the relationship between thermal stress and wet stress, the relationship between temperature and latent heat of vaporization, the fluid mass equation, and the fluid continuity equation.

[0033] In a second aspect, an embodiment of the present application provides a device for valve structure simulation based on finite element analysis, comprising a memory, a processor, and a computer program stored in the memory and executable in the processor, characterized in that the computer program comprises steps corresponding to the above-mentioned valve structure simulation method based on finite element analysis.

[0034] In a third aspect, an embodiment of the present application provides a readable storage medium, which stores a terminal program, and is characterized in that when the terminal program is executed, it performs the valve structure simulation method based on finite element analysis as described above.

[0035] It can be seen from the above technical solution that the beneficial effects of this application are:

[0036] 1. The method of this application designs a geometric model based on the structural parameters of the valve structure and subdivides it into multiple components. This can accurately restore the actual structure of the valve. The material parameters of each component are defined and associated, so that the model can reflect the actual material properties, laying the foundation for subsequent precise analysis. Valve seats and valve cores made of different materials have different mechanical properties and thermal conductivity. This allows accurate parameter setting to simulate the interaction under actual working conditions. By meshing and using them as simulation objects, it can improve computational efficiency and accuracy, ensuring the accuracy of the description of complex valve structures. A coupled model of fluid flow, moisture transfer, heat transfer, and mechanical deformation is established, comprehensively considering factors such as fluid flow, moisture heat transfer, and structural mechanics that affect the valve structure. Finally, the coupled model is imported into simulation model software for simulation solution analysis. The simulation is carried out by establishing multi-physics field coupling relationships, integrating advanced methods and optimization strategies. The valve structure generates heat during operation, humidity changes may affect material properties, and mechanical deformation interacts with other factors. This multi-physics field coupling analysis can simulate the valve operating state in a more real-time manner and obtain more accurate results, thereby improving the real-time performance of the multi-physics field coupling simulation, ensuring work efficiency and improving accuracy. By determining the simulation results, grid cells with failure risks can be identified, accurately displaying problems in the valve structure design, facilitating adjustments and ensuring the accuracy of the valve structure.

[0037] 2. The device of the present application can efficiently execute the structural simulation method based on finite element analysis through a computer program stored in the memory and running on the processor, and can perform complex numerical calculations more quickly, which greatly improves the calculation efficiency of the simulation and saves the time cost of valve structure design.

[0038] 3. The readable storage medium of the present application can conveniently store and transfer the above-mentioned valve structure simulation method based on finite element analysis in the form of software, which is convenient for the development, design and use of the valve structure. The readable storage medium can execute the above-mentioned valve structure simulation method based on finite element analysis on different devices, which is convenient for transfer when the equipment needs to be replaced, and can meet the use requirements of different devices. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for the description of the embodiments one by one. Obviously, the drawings described below are some embodiments of the present application. For those skilled in the art, other embodiments and drawings can be obtained based on these drawings without inventive work. The block diagrams shown in the drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities can be implemented in software form, or in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or microcontroller devices. The flowcharts shown in the drawings are only illustrative and do not necessarily include all content and operations / steps, nor do they have to be executed in the order described. For example, some operations / steps can be decomposed, while some operations / steps can be combined or partially combined. Therefore, the actual execution order may change according to actual circumstances.

[0040] Figure 1 A schematic diagram showing the steps of an embodiment of a valve structure simulation method based on finite element analysis of the present invention is shown;

[0041] Figure 2 A schematic diagram showing the relationship between the coupled model of multiple physical fields;

[0042] Figure 3 A schematic diagram of a deformation simulation model of a geometric model in a valve structure simulation method based on finite element analysis of the present invention is shown. DETAILED DESCRIPTION

[0043] The following will be combined with the accompanying drawings corresponding to the specific embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. The following detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the present application for protection. The described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, various different configurations can be arranged and designed. All other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0044] The present application is described below with reference to specific embodiments and with reference to the accompanying drawings:

[0045] Please refer to Figure 1 The first aspect of the present application provides an embodiment of a valve structure simulation method based on finite element analysis, including:

[0046] S1. Design a geometric model based on the structural parameters of the valve structure, and divide the geometric model into multiple components. The multiple components are divided according to the various combination relationships or functional relationships of the valve structure. The valve structures involved in the embodiments of the present application include various types, including both mechanical valves and electric valves. Mechanical valves include ball valves, butterfly valves, check valves, etc. Electric valves include solenoid valves, servo valves, etc. Servo valves include electro-hydraulic, pneumatic or hydraulic types. Electro-hydraulic types include proportional valves and jet tube valves. The geometric model of the valve structure is preliminarily designed based on the performance requirements, such as dimensional requirements, material requirements, and various physical parameter requirements. Physical parameters such as vibration performance, high and low temperature, acceleration, etc. are determined according to the product application scenario requirements, and load specifications and design specifications must also be considered.

[0047] S2. Define the material parameters of each component and associate the material parameters with each component. This will associate the material parameters with the different structures of the geometric model of the valve structure. For example, the valve core uses cemented carbide and the valve sleeve uses stainless steel. It also includes material parameters such as hardness and elastic modulus, and associates each component with the material parameters. This ensures that the designed geometric model can meet the design requirements and can also be used to evaluate the performance of the materials of each component in complex environments.

[0048] S3. Use finite element simulation model software to perform finite element meshing on the geometric model to form mesh units, and use all mesh units as simulation analysis objects; meshing is a common method in existing simulation analysis, and this embodiment also adopts a conventional meshing method. Through meshing, complex geometric shapes can be divided into multiple simple units, thereby achieving approximate solutions to complex geometric shapes.

[0049] S4. Establish a coupling model of the effects of fluid flow, humidity transfer, heat transfer and mechanical deformation on the valve structure; set boundary conditions and initial conditions; introduce the effects of different environmental information on the valve structure. By establishing the above coupling model, the important parameters affecting the valve structure can be taken into account, the establishment of the coupling model of multiple physical fields can be completed, and an accurate simulation model of the valve structure can be obtained more comprehensively. Through the boundary conditions and initial conditions, the constraints of the coupling model can be determined to ensure the uniqueness and stability of the solution and simulate the stress and deformation conditions in the actual engineering environment.

[0050] S5. Import the coupling model into the simulation model software to perform simulation analysis of multi-physical field coupling; different physical fields have their own control equations, and these physical fields are interconnected. By importing the coupling model into the simulation model software, each physical field and geometric model can be linked. This process is a conventional simulation method. With the help of the powerful computing function of the simulation model software, the distribution and relationship of each physical field can be quickly calculated, thereby providing a basis for design. The finite element simulation model software can use conventional analysis software such as ANSYS, SolidWorks, etc.

[0051] S6. Based on the simulation results, identify mesh units in the geometric model that are at risk of failure and optimize the design. The simulation model software can quickly calculate the simulation results, and based on the simulation results, judge the mesh units in the geometric model that are at risk of failure. Adjust the geometric model positions corresponding to the mesh units at risk of failure, first mark the mesh units at risk of failure as weak mesh units, then readjust the design of the geometric model, optimize, and return to the division of finite element mesh units to continue. Then repeat the above steps S5 and S6 until each mesh unit of the final geometric model meets the design standards. The standards may include various parameters such as stress, safety factor, displacement, deformation, fatigue, etc., and ultimately no weak mesh units exist. Weak mesh units refer to areas in the valve structure that are prone to structural failure and are key areas that should be optimized. Generally, simulations judge from several points: high stress value areas, i.e., stress concentration areas, areas with excessive deformation, fatigue hot spots, and areas prone to fatigue damage. If weak mesh units exist, the valve structure will be deformed or damaged.

[0052] The method of the present application designs a geometric model based on the structural parameters of the valve structure and subdivides multiple components, which can accurately restore the actual structure of the valve. The material parameters of each component are defined and associated, so that the model can reflect the actual material properties, laying the foundation for subsequent accurate analysis. The valve seat and valve core of different materials have different mechanical properties and thermal conductivity. This allows accurate parameter setting to simulate the interaction under actual working conditions. By meshing and using them as simulation objects, it can improve computational efficiency and accuracy, ensuring the accuracy of the description of complex valve structures. A coupled model of fluid flow, moisture transfer, heat transfer and mechanical deformation is established, comprehensively considering factors such as fluid flow, moisture heat transfer and structural mechanics that affect the valve structure. Finally, the coupled model is imported into the simulation model software for simulation solution analysis. By establishing a multi-physics coupling relationship for simulation, advanced methods and optimization strategies are integrated. The valve structure generates heat during operation, humidity changes may affect material properties, and mechanical deformation interacts with other factors. Through this multi-physics coupling analysis, the working state of the valve can be simulated more in real time, and more accurate results can be obtained, thereby improving the real-time performance of the multi-physics coupling simulation, ensuring work efficiency and improving accuracy. By determining the simulation results, grid cells with failure risks can be identified, accurately displaying problems in the valve structure design, facilitating adjustments and ensuring the accuracy of the valve structure.

[0053] At the same time, by determining the parameters of the optimal valve structure under different environments, this application can conduct in-depth analysis of the valve structure and evaluate the structural performance, and then optimize and improve the design and performance of the valve structure to ensure that the valve structure meets the performance requirements in different application scenarios; this application also comprehensively considers the mechanical properties of the valve structure and various environmental factors, and can design and simulate more accurately, improve the design efficiency of the valve structure, and reduce the development cost of the valve structure.

[0054] In some embodiments, the finite element meshing of the geometric model specifically includes: 1. Adaptive mesh refinement: dynamically adjusting the mesh density based on the local physical field strength to ensure high resolution in key areas. Dynamic adjustment of the mesh is achieved by the following methods: local encryption: increasing the mesh resolution in areas where the solution changes dramatically; local coarsening: reducing the mesh resolution in areas where the solution changes slowly; hierarchical meshing: dividing the computational domain into multiple levels of meshes, each with a different resolution. The key to adaptive mesh refinement is how to determine which areas need to be encrypted or coarsened, which is usually based on error estimation or the gradient of physical quantities. 2. Refinement criteria: error estimation, which determines which areas require finer meshes by estimating local errors (estimated errors by comparing solutions on different meshes); physical field characteristics, for areas that are sensitive to specific physical phenomena, such as high gradient areas in fluids, stress concentration points in structures, etc., finer meshes should be used; user-defined criteria, manually specifying certain areas for refinement based on engineering experience and design requirements. Through the above-mentioned mesh refinement, the mesh density is dynamically adjusted based on the local physical field strength to ensure high resolution in key areas; at the same time, it can support multi-scale modeling, using coarse and fine meshes simultaneously in the same model, balancing computational cost and accuracy requirements.

[0055] Please refer to Figure 2 In some embodiments, the coupling model includes: the continuity equation of the fluid, and the heat conduction equation, the moisture diffusion equation, and the dynamic equation. The fluid continuity equation is used to describe the motion state of the fluid flow; the heat conduction equation is used to describe the temperature changes of each component caused by the environmental heat source and its own Joule heat generation; the moisture diffusion equation is used to describe the humidity changes of each component caused by the environmental humidity and its own moisture diffusion; the dynamic equation is used to describe the changes in the displacement, velocity, acceleration and stress of each component of the valve structure. This coupling model involves the interaction between fluid flow, humidity transmission, heat transfer and structural deformation, forming a coupled model of flow, moisture, heat and force. In order to describe the coupling relationship of these physical fields, it is usually necessary to solve a set of coupled partial differential equations. The numerical solution method based on the finite element method is an effective tool for dealing with such complex problems.

[0056] In some embodiments, the dynamic equation is described based on a multi-degree-of-freedom system as the relationship between the displacement vector, velocity vector, acceleration vector, and external force vector, where the coefficients of the displacement vector, velocity vector, and acceleration vector are the mass matrix, damping matrix, and stiffness matrix, respectively. Vibration, shock, and acceleration simulation of valve structures are often used to analyze the response characteristics of the structure under dynamic loads. Based on the dynamic equation, the changes in physical quantities such as displacement, velocity, acceleration, and stress of the system are described. The motion equation based on the multi-degree-of-freedom system (MDOF) can be described as:

[0057]

[0058] Where M is the mass matrix; C is the damping matrix; K is the stiffness matrix; X(t) is the displacement vector; and are the second-order derivative and first-order derivative of the displacement vector respectively; F(t) is the external force vector.

[0059] When there is no external excitation, F(t) = 0; when there is vibration excitation, F(t) = F0cos(ω0t); when there is impact excitation, When acceleration is applied, F(t)=at.

[0060] Where ω0 is the excitation frequency; F0 is the reference value of the external force vector, which is used as the excitation amplitude; a is the acceleration excitation value; and t is time.

[0061] In some embodiments, when performing simulation analysis of multi-physics field coupling, the equations corresponding to the coupled model are solved to obtain simulation results for temperature, humidity, stress, and displacement distributions. Specifically, the equations in the coupled model are solved simultaneously to ultimately calculate the distribution of various physical quantities on each component. For example, using a servo valve as an example, simulation analysis using the coupled model ensures the performance of the valve structure's geometric model under different operating conditions during the design phase, thereby determining whether the servo valve can meet the system's flow and pressure control requirements.

[0062] In some embodiments, when solving the equations corresponding to the coupled model, a finite element numerical solution is used, including:

[0063] S51. Introduce the energy conservation equation, mass conservation equation, momentum conservation equation, and coupling mechanism for multi-physics field coupling. The coupling mechanism includes the relationship between the four physical fields of flow, moisture, heat, and force. The thermal simulation calculation of the servo valve uses the heat conduction equation to describe the temperature changes of the component itself caused by external heat sources and the Joule heating generated by the component itself. The specific expression is as follows:

[0064]

[0065] Where T is the temperature of each structural component itself; k is the thermal conductivity of the material; ρ2 is the material density of each component; C p is the specific heat capacity of the material at normal pressure; Q is the external heat source; is the temperature gradient; T0 is the constant temperature at the boundary; and t is time.

[0066] The wet simulation calculation of the servo valve uses the moisture diffusion equation to describe the humidity changes of the component itself caused by external humidity and its own moisture diffusion.

[0067]

[0068] Where ω is the humidity of each component; D is the moisture diffusion coefficient; S 湿 is the moisture source term; t is time; is the humidity gradient.

[0069] Because the servo valve itself exists in a hot and humid environment, its material properties are affected by temperature and humidity, leading to changes in their inherent parameters. Mechanical analysis is required to analyze the physical structural characteristics of each component. To describe the coupling relationships between these physical fields, the corresponding governing equations are derived, and appropriate coupling terms are introduced. The following is a derivation of the basic theoretical formulas for hot and humid coupling simulation, covering conservation of energy, mass, momentum, and coupling effects.

[0070] The basic equation of fluid flow is:

[0071]

[0072] Where V is the velocity field of the fluid; q is the flow rate; k f is the thermal conductivity of the material; μ is the viscosity of the fluid; and P is the fluid pressure.

[0073] The continuity equation of the fluid ensures the conservation of mass of the fluid and is:

[0074]

[0075] Where V is the velocity field of the fluid; q is the flow rate, which can be used as a fluid source term; Φ is the porosity; t is time; and ρ1 is the fluid density.

[0076] Next, we construct the basic equations of structural mechanics. The strain tensor ε can be derived from the displacement gradient as follows:

[0077]

[0078] Where T is the absolute temperature; u is the displacement vector; is the displacement gradient.

[0079] The constitutive relationship between stress and strain is:

[0080] σ=C:ε

[0081] Where C is the elastic stiffness tensor; ε is the strain tensor; and σ is the stress.

[0082] For the coupled flow, moisture, heat and force problem, the energy conservation equation is as follows:

[0083]

[0084] Where k(T,ω,P) is the thermal conductivity, which is affected by temperature T, humidity w, and pressure P; Q int is the internal heat source term; Q latent is the latent heat source term, which represents the heat released or absorbed due to phase change per unit volume; V is the velocity field of the fluid; ρ1 is the fluid density; is the temperature gradient; P is the fluid pressure; ω is the humidity of each component itself; C p is the specific heat capacity of the material at normal pressure; t is time.

[0085] When moisture undergoes a phase change, it releases or absorbs energy that affects the energy balance and temperature distribution of the system.

[0086]

[0087] Among them, Q latent is the latent heat source term, which represents the heat released or absorbed due to phase change per unit volume; L is the latent heat of vaporization of water; ω is the humidity of each component itself; and t is time.

[0088] For the coupled flow, moisture, heat and force problem, the mass conservation equation is:

[0089]

[0090] Where D(T,ω,P) is the moisture diffusion coefficient, which is affected by temperature T and humidity ω; V is the fluid velocity; P is the fluid pressure; t is time; S 湿 is the moisture source term.

[0091] Based on the Arrhenius model and combining the nonlinear effects of temperature and humidity, the following formulas for the thermal conductivity k(T, ω, P) and the diffusion coefficient D(T, ω, P) were obtained. Their reliability was confirmed through simulation calculations and experimental verification.

[0092]

[0093] Among them, E a is the activation energy; represents the effect of temperature on thermal conductivity; R is the gas constant, 8.314 J / (mol·K); n and m are empirical coefficients, usually 1 or 2; α 敏is the pressure sensitivity coefficient; T is the temperature; ω is the humidity of each component itself; P is the fluid pressure; D0 is the reference diffusion coefficient; k0 is the reference thermal conductivity.

[0094] For the coupled problem of fluid, moisture, heat and force, the momentum conservation equation is:

[0095]

[0096] Among them, b is the body force; α is the thermal expansion coefficient; β is the hygroscopic expansion coefficient; γ is the pressure sensitivity coefficient; I is the unit tensor; C is the elastic stiffness tensor; ε is the strain tensor; T is the temperature; ω is the humidity of each component itself; ρ is the density, which can be either the fluid density ρ1 or the material density of each component ρ2; u is the displacement; t is time.

[0097] S52. Discretize the continuous multi-physics field into a finite number of physics field units corresponding to the grid units.

[0098] In order to numerically solve these complex coupling problems, the finite element method (FEM) is usually used to discretize the continuous physical field into a finite number of physical field units. These physical field units correspond to grid units. By simulating the physical field units, the physical field corresponding to the geometric model corresponding to each grid unit can be simulated, thereby achieving comprehensive coupling of multiple physical fields. Discrete multiple physical fields into a finite number of physical field units, and different solution methods such as the finite element method, finite difference method, finite volume method, and boundary element method can be used to discretize multiple physical fields. For the valve structure of this application with a more complex geometry and the coupling problem of multiple physical fields, the finite element method is used for analysis.

[0099] S53. Solve the discretized algebraic equations of each physical field unit to obtain an approximate solution.

[0100] The discretized algebraic equations for each physical field unit are nonlinear partial differential equations. To solve these nonlinear partial differential equations, weak forms and the Galerkin method are typically used. The weak form is obtained by multiplying the original partial differential equation by an appropriate test function and integrating it over the solution domain. The weak form does not require the equation to be strictly satisfied at every point, but rather satisfied in an integral sense, which reduces the requirement for smoothness of the solution. The weak form facilitates the conversion of partial differential equations into integral equations, which can then be solved, simplifying the calculations. The Galerkin method is a method for solving the weak form. By selecting basis functions, the integral equation is converted into a system of algebraic equations, providing a more suitable solution for the numerical solution of multi-physics field problems.

[0101] In some embodiments, solving the discretized algebraic equations specifically includes:

[0102] S531. By introducing test functions, the partial differential equations of the energy conservation equation, mass conservation equation, and momentum conservation equation are converted into integral forms to form weak form expressions respectively.

[0103] The test functions introduced are: δT, δω and δu, which transform the partial differential equation into an integral form, including:

[0104] Weak form of the energy conservation equation:

[0105]

[0106] Where Ω is the integration area; ρ is the density, which can be either the fluid density ρ1 or the material density of each component ρ2; C p is the specific heat capacity of the material at normal pressure; T is the ambient temperature; t is time; V is the velocity field of the fluid; ω is the humidity of each component itself; P is the fluid pressure; Q int is the internal heat source term; Q latent is the latent heat source term; is the temperature convection term due to fluid flow.

[0107] Weak form of the mass conservation equation:

[0108]

[0109] Where ω is the humidity of each component itself; t is time; D(T, ω, P) is the moisture diffusion coefficient; is the gradient of the humidity process quantity; is the moisture convection term caused by fluid flow; S 湿 is the moisture source term.

[0110] Weak form of the momentum conservation equation:

[0111]

[0112] Where C is the elastic stiffness tensor; ε is the strain tensor; T is the temperature; ω is the humidity of each component; ρ is the density, which can be either the fluid density ρ1 or the material density ρ2 of each component; u is the displacement; t is the time; α 热 is the thermal expansion coefficient; β is the hygroscopic expansion coefficient; γ is the pressure sensitivity coefficient; I is the unit tensor; and b is the body force.

[0113] S532. Substitute the temperature field corresponding to each grid cell expressed by the interpolation function into the weak form corresponding to the energy conservation equation, substitute the humidity field into the weak form corresponding to the mass conservation equation, and substitute the displacement field into the weak form corresponding to the momentum conservation equation to obtain the discretized equations respectively.

[0114] The temperature field is expressed by the interpolation function T(x) and substituted into the weak form to obtain the discretized energy conservation equation:

[0115]

[0116] Where ρ is the density, which can be either the fluid density ρ1 or the material density of each component ρ2; C p is the specific heat capacity of the material at normal pressure; N i 、N j is the shape function at a node position, i and j are integer indices representing nodes; V is the velocity field of the fluid; k(T, ω, P) is the thermal conductivity of the material; Q int is the internal heat source term; Q latent is the latent heat source term; M is the mass matrix; C 阻 is the damping matrix; K is the stiffness matrix; T is the ambient temperature.

[0117] The humidity field is expressed by the interpolation function ω(x) and substituted into the weak form to obtain the discretized mass conservation equation:

[0118]

[0119] Among them, N i 、N j is the shape function at a certain node position; ω is the humidity of each component itself; S 湿 is the moisture source term; M is the mass matrix; K is the stiffness matrix.

[0120] The displacement field u(x) is expressed as an interpolation function and substituted into the weak form to obtain the discretized momentum conservation equation:

[0121]

[0122]

[0123] Where ρ is the density, which can be either the fluid density ρ1 or the material density ρ2 of each component; N i 、N j is the shape function at a certain node position; C 阻 is the damping matrix; u is the displacement vector; b is the body force; M is the mass matrix; K is the stiffness matrix; and t is time.

[0124] S533. Use an interpolation function to represent the temperature, humidity, and displacement of each grid cell.

[0125] Use linear or nonlinear shape functions to interpolate the physical field. Specifically, we assume that we have a finite element mesh with n nodes, and the unknown quantity at each node (such as temperature T i 、Humidityi and displacement u i ) can be expressed by the interpolation function:

[0126]

[0127] Among them, N i (x) is the shape function, satisfying N i (x) = δ ij (Kronecker delta); x represents the node position, i and j are integer indices representing the node.

[0128] S534. Perform coupled simulation to solve each discretized equation.

[0129] When performing coupled simulation solutions on discretized equations, the main problem is dealing with a system of equations consisting of multiple interrelated equations. A variety of solutions can be used, such as direct simultaneous solution, iterative solution, etc. For smaller-scale and highly coupled systems of equations, all discretized equations can be combined to form a large system of equations for solution; for large-scale systems of equations, especially when the coupling between equations is not very tight, or the coupling relationship can be simplified by appropriate methods, an iterative solution method can be used. This application uses the above-mentioned intelligent coupling algorithm, a new iterative solver, combined with preconditioning technology and relaxation factor optimization to accelerate the convergence process.

[0130] In some embodiments, when performing coupled simulation to solve each discretized equation:

[0131] The Newmark method is used to discretize the momentum conservation equation, and then coupled simulation is performed with the discretized energy conservation equation and the discretized mass conservation equation to solve it. The Newmark method is a commonly used implicit time integration method, which is widely used in the numerical solution of structural dynamics and fluid-solid coupling problems. It is suitable for second-order ordinary differential equations (such as the momentum conservation equation) and has good numerical stability and accuracy. In the fluid-humidity-heat-force coupling problem, the Newmark method can be used to discretize the momentum conservation equation, and at the same time, combined with the discretization methods of other physical fields (such as temperature field and humidity field), to realize multi-physical field coupling simulation. Specifically:

[0132] When discretizing the momentum conservation equation, the Newmark method is used to discretize the momentum conservation equation and obtain the time step t n+1 The acceleration, displacement and velocity at .

[0133]

[0134] Among them, u n 、 They represent displacement, velocity, and acceleration respectively; Δt is the time step; β and γ are Newmark method parameters, usually β = 0.25 and γ = 0.5.

[0135] Get time step t n+1 The acceleration equation at is:

[0136]

[0137] Where M is the mass matrix; K is the stiffness matrix; C 阻 is the damping matrix; is the acceleration; Δt is the time step; γ is the pressure sensitivity coefficient; u n is unique; β is the coefficient of expansion; F n+1 is the external force vector at time step n+1.

[0138] Then discretize the energy conservation equation and mass conservation equation.

[0139] The discretization format of the energy conservation equation is:

[0140]

[0141] Where ρ is the density, which can be either the fluid density ρ1 or the material density of each component ρ2; C p is the specific heat capacity of the material at normal pressure; T n is the temperature at time step n; Δt is the time step; V is the velocity field; k(T n+1 ,ω n+1 , P n+1 ) is the thermal conductivity of the material; q int is the internal heat source term; Q latent is the latent heat source term.

[0142] The discretization format of the mass conservation equation is:

[0143]

[0144] Among them, ω n is the humidity in n time steps; Δt is the time step; D(T n+1 ,ω n+1 , P n+1 ) is the moisture diffusion coefficient at time step n+1; S 湿 is the moisture source term.

[0145] In order to solve the above governing equations, appropriate boundary conditions and initial conditions need to be specified. Initial conditions include the initial temperature, humidity, displacement, and fluid pressure distribution inside the material. Common boundary conditions include:

[0146] Temperature boundary condition: specifies the temperature or heat flux on the material surface; for example, the ambient temperature and the material surface temperature are 30°C.

[0147] Humidity boundary condition: Specifies the humidity or moisture flux on the material surface; for example, the relative humidity of the environment and the material surface is 85%.

[0148] Mechanical boundary conditions: specify the displacement or stress on the material surface; for example, use the fixed connection parts of the servo valve as fixed constraint positions; set the wall as a no-slip wall condition.

[0149] Fluid boundary conditions: specify the fluid inlet pressure, outlet pressure or flow rate, such as the fluid inlet pressure is 7MPa and the fluid outlet pressure is 0MPa.

[0150] The loading conditions such as acceleration, impact and vibration are converted into action load mode for loading. For example, the overall acceleration is 50g, and the directions are respectively acting in the X, Y and Z axes.

[0151] The above discretized momentum conservation equation is coupled with the discretized energy conservation equation and the discretized mass conservation equation for simulation and solution. After the above arrangement, we get:

[0152]

[0153] Where: M is the mass matrix; K is the stiffness matrix; C 阻 is the damping matrix; γ is the pressure sensitivity coefficient; Δt is the time step; β is the expansion coefficient;

[0154] In some embodiments, when performing coupled simulation solutions, a coupling mechanism is introduced, including: the relationship between thermal stress and wet stress, the relationship between temperature and latent heat of vaporization, the fluid mass equation, and the fluid continuity equation.

[0155] In coupled hygrothermal-mechanical problems, changes in temperature, humidity, and stress affect each other and jointly affect the mechanical properties of the material. Therefore, coupling terms are needed to describe these interactions.

[0156] Mechanical analysis is mainly based on elastic mechanics, and its related equations include equilibrium equations, displacement equations, and thermal stress-thermal strain relationships as follows:

[0157]

[0158] Where ρ is the density, which can be either the fluid density ρ1 or the material density ρ2 of each component; μ is the damping coefficient; is the displacement vector; is the applied stress; It is an elastic relationship. and is a differential operator, It is thermal strain or wet strain caused by changes in temperature or humidity; is the stress vector; It is strain; It's heat stress.

[0159] Based on the law of conservation of energy and the coupled hygrothermal model, changes in temperature and humidity interact, generating thermal and hygroscopic stresses within the material, which in turn affect its mechanical properties. For example, rising temperature can increase the material's hygroscopicity, while rising humidity can cause the material to swell, generating hygroscopic stresses. These stresses, combined with mechanical stresses, create a complex coupling effect. Thermal and hygroscopic stresses can establish a coupling mechanism through hygrothermal stresses:

[0160] σ thermo-hygral =α 热 TI+βωI

[0161] Among them, σ thermo-hygral is the moist-thermal coupling stress; α 热 is the coefficient of thermal expansion; β is the coefficient of hygroscopic expansion; I is the unit tensor; T is the ambient temperature; ω is the humidity of each component itself.

[0162] When the moisture in the material undergoes a phase change (such as evaporation or condensation), it is accompanied by the release or absorption of latent heat, which affects the distribution of the temperature field. The phase change effect can be established by introducing the latent heat term to establish a coupling mechanism:

[0163]

[0164] Among them, Q latent is the latent heat source term; L is the latent heat of vaporization of water; ω is the humidity of each component itself; and t is time.

[0165] Fluid flow affects the deformation of solid structures, and vice versa; fluid-solid interaction can be realized through the fluid pressure term P and the solid displacement term u in the momentum conservation equation.

[0166] When analyzing the flow path of a servo valve, it is necessary to solve the governing equation of fluid flow. Based on the law of mass conservation, the continuity equation can be expressed as:

[0167]

[0168] in, ρ1 represents the velocity and density of the fluid respectively; t is time.

[0169]

[0170] Where ρ1 is the fluid density; P is the fluid pressure, υ is the fluid viscosity, represents the vector acceleration due to gravity; is the displacement vector; t is the time.

[0171] In some embodiments, parametric optimization aims to achieve specific goals (such as minimizing stress and maximizing efficiency) by optimizing design parameters in multi-physics coupled simulations. Genetic Algorithm (GA) is a commonly used global optimization method suitable for finding the optimal combination of design parameters for complex problems. For servo valve design, it is necessary to optimize its geometric parameters (such as flow channel width and valve body thickness) so that the servo valve's performance indicators (such as flow rate, pressure loss, structural stress, etc.) are optimized under given operating conditions.

[0172] Genetic algorithms imitate natural selection and genetic mechanisms, searching for the optimal solution by simulating the evolutionary process. Its main steps include: Initializing the population: randomly generating a set of initial solutions (individuals), each individual representing a possible design scheme. Evaluating fitness: calculating the fitness value of each individual based on a predefined objective function. Selection: selecting individuals to enter the next generation based on fitness values. Crossover: generating new individuals by exchanging some genetic information. Mutation: randomly changing certain genes to increase population diversity. Termination condition: stopping iteration when certain conditions are met (such as reaching the maximum number of iterations or finding a satisfactory solution). The specific implementation steps are as follows:

[0173] First, define an objective function to evaluate the pros and cons of each design solution. For servo valve design, optimize its design parameters, including flow channel width and valve body thickness. The goal can be to comprehensively consider multiple performance indicators, such as minimizing pressure loss and maximizing flow rate. Then, create an initial population, with each individual consisting of a set of design parameters. Individuals are selected for the next generation based on their fitness values. New individuals are generated by exchanging some genetic information. Some genes are randomly changed to increase population diversity. Each of these steps is integrated into a main loop until the termination condition is met. Using the genetic algorithm and parameter optimization method described above, the optimal combination of design parameters can be automatically found.

[0174] In the above step S6, when determining the mesh cells with failure risk in the geometric model, the mesh cells with failure risk are marked as weak mesh cells. The determination of whether a mesh cell is a weak mesh cell is usually based on thresholds of key parameters such as stress, strain or deformation.

[0175] Based on the commonly used Von Mises stress and safety factor M of valves s , displacement and deformation, and fatigue analysis to determine which mesh elements are considered weak, including:

[0176] 1. Von Mises stress: If the Von Mises stress σ of a grid cell VM Exceeds the yield strength of the material σ y , then the unit is considered a weak mesh unit. For example, according to the material properties of a valve, the yield strength is 1550MPa. The stress of the mesh unit is compared with the yield strength. If the stress is less than the yield strength, the design is qualified. Otherwise, it is marked as a weak mesh unit.

[0177] 2. Safety margin (M s ): The ratio of design stress to actual working stress, used to measure the safety margin of the structure.

[0178] Normally, M s Should be greater than 1; for critical components and high-risk applications, M is recommended s At least 1.2 to 1.35; if the M of some units s If the value is smaller than this range, these cells can be regarded as weak mesh cells.

[0179] Ensure the safety margin M of the servo valve s At least above 1.2.

[0180]

[0181] Where [σ] is the allowable stress, MPa. The strength limit is used as the benchmark for brittle materials, and the yield limit is used as the benchmark for plastic materials. Its value is 1550 MPa. f is the safety factor, dimensionless. The safety factor for quasi-static loads is generally 1.5, and the safety factor for vibration loads is generally greater than or equal to 1.2. Its value is 1.2. max is the maximum stress value, MPa; to meet the safety margin of the servo valve, the maximum stress value of the grid unit should be less than 1076.39 MPa.

[0182] 3. Displacement and deformation: Set reasonable displacement limits based on design requirements and operating conditions. For example, in a precision servo system, any displacement exceeding the micron level may be unacceptable.

[0183] Through finite element analysis, the deformation of each component of the servo valve under the action of the coupled stress field is calculated and compared with the design displacement field of each part obtained according to the design standards. If it is less than, the design of the part is qualified, otherwise it is unqualified.

[0184] 4. Fatigue Analysis: For components such as servo valves that are frequently subjected to cyclic loads, fatigue damage must also be considered. If the predicted service life is lower than the expected service life, there is a risk of fatigue failure, and the relevant mesh elements should be considered weak.

[0185] For this servo valve, the design life is 107 The SN curve (stress-life curve) is used to evaluate the fatigue life of the material under cyclic load. If the predicted fatigue life is lower than the design requirement, it is considered that there is a risk of fatigue failure and is marked as a weak mesh element.

[0186] Determining which mesh elements in a servo valve are weak depends on the specific engineering requirements, material properties, and application environment. By combining the aforementioned methods with specialized knowledge and conducting detailed analysis, we can effectively identify and improve mesh elements at risk of failure in the design, allowing for further design optimization.

[0187] As in the embodiment of the present application, the deformation of the servo valve structure under the action of comprehensive stress is calculated by finite element analysis, and the maximum displacement X is obtained. max , compared with the design displacement X of each structural component, if it is less than, the structural component design is qualified, otherwise the design is unqualified, and the safety margin of each structural component is at least greater than 1.2. Figure 3 As shown in the figure, the final calculation results of the servo valve body are as follows: the horizontal forward displacement of the valve body is 0.0049mm; the horizontal lateral displacement is 0.0012mm; the vertical downward displacement is 0.0019mm; and the maximum stress peak on the valve body is 67.078MPa. Based on its material properties (the yield stress of aluminum alloy is 260MPa), the safety margin of the maximum stress can be calculated to be 3, further verifying the reliability of the servo valve assembly body in terms of strength and rigidity. Under some abnormal operating conditions, such as extreme high and low temperatures, random vibration, and large acceleration, the servo valve will have certain limitations. Therefore, when designing the servo valve structure, sufficient margin must be left for each servo valve component to ensure the structural strength and stability of the servo system under abnormal operating conditions.

[0188] In a second aspect of the present application, an embodiment is provided of a device for simulating a valve structure based on finite element analysis, comprising a memory, a processor, and a computer program stored in the memory and executable in the processor. The computer program comprises the steps corresponding to the above-mentioned method for simulating a valve structure based on finite element analysis. When the processor executes the computer program, the steps in the above-mentioned method embodiment are implemented, for example Figure 1 Alternatively, when a processor executes the computer program, the functions of the modules / units in the embodiment of the apparatus for valve structure simulation based on finite element analysis are realized. It should be understood that the apparatus of each embodiment of the present application can be realized based on a memory and a processor, wherein each memory is used to store a computer program for executing the above-mentioned method of the present application, and the processor executes the computer program, so that the apparatus for valve structure simulation based on finite element analysis realizes the method of each embodiment described above.

[0189] This application uses a computer program stored in a memory and running on a processor to efficiently execute a structural simulation method based on finite element analysis, and can perform complex numerical calculations more quickly, thereby greatly improving the computational efficiency of the simulation and saving the time cost of valve structure design.

[0190] In some embodiments, the above-mentioned device for valve structure simulation based on finite element analysis can be a desktop computer, a notebook computer, an industrial computer, a PDA, a tablet computer or other mobile terminal, as well as a computer device such as a cloud server, and is not limited to which operating system is installed. The device for valve structure simulation based on finite element analysis may include, but is not limited to: a processor and a memory. Those skilled in the art will understand that the example of the device for valve structure simulation based on finite element analysis does not constitute a limitation on the device for valve structure simulation based on finite element analysis, and may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the device for valve structure simulation based on finite element analysis may also include input devices, output devices, network access devices, buses, etc.

[0191] The third embodiment of the present application provides a readable storage medium, which stores a terminal program. When the terminal program is executed, it performs the valve structure simulation method based on finite element analysis as described above. If the integrated module / unit of the valve structure simulation method based on finite element analysis is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a readable storage medium. Based on this understanding, the present application implements all or part of the process of the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a readable storage medium. When the computer program is executed by a processor, the steps of each method in the above-mentioned embodiment can be implemented. When the computer program is executed by the processor, the specific implementation of each step and the technical effect produced are the same as those in the above-mentioned method embodiment. For the sake of brief description, for parts not mentioned in this embodiment, reference can be made to the corresponding content in the above-mentioned method embodiment.

[0192] The present application uses a readable storage medium to conveniently store and transfer the above-mentioned valve structure simulation method based on finite element analysis in the form of software, which is convenient for the development, design and use of the valve structure. The readable storage medium can execute the above-mentioned valve structure simulation method based on finite element analysis on different devices, which is convenient for transfer when the equipment needs to be replaced, and can meet the usage requirements of different devices.

[0193] Regarding the specific implementation of this application, it should be noted that:

[0194] In the description of the present application, the embodiments of the present application can be implemented with electronic hardware, computer program products, or a combination of computer software and electronic hardware, and are described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of processes and / or boxes in the flowchart and / or block diagram, can be implemented by instructions. The above-mentioned computer program can be loaded on a computer and stored in a computer-readable storage medium that can guide a computer or other programmable data processing device to work in a specific manner, so that the instructions of the computer program execute a series of operating steps to produce computer-implemented processing, produce a manufactured product including an instruction device, and implement it in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 The computer-readable storage medium may include any entity or device capable of carrying computer program code, such as a USB flash drive, a mobile hard drive, a magnetic disk, an optical disk, a computer memory, a read-only memory, a random access memory, an electrical carrier signal, a telecommunications signal, a software distribution medium, or other recording medium. Examples include a hard disk, a memory, an optical storage device, a plug-in hard disk, a memory card, a secure digital card, a flash memory card, or at least one magnetic disk storage device, a flash memory device, or other volatile solid-state storage device.

[0195] In the description of this application, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, apparatus, or readable storage medium comprising a series of elements includes not only those elements but also other elements not explicitly listed that are consistent with the concept of this application, or elements inherent to such process, method, apparatus, or readable storage medium. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional elements in the process, method, apparatus, or readable storage medium comprising the element.

[0196] In the description of this application, reference to the terms "some embodiments", "optional embodiments", "examples", "specific examples", "optional examples" or "optional embodiments" and the like means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of this application, but does not mean that these embodiments illustrate and describe all possible forms of the invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification.

[0197] In addition, the technical solutions between the various embodiments can be combined with each other, but this must be based on the fact that they can be implemented by ordinary technicians in this field. Although the embodiments of the present application have been shown and described, these embodiments can be changed, modified, replaced, and modified in various ways without departing from the principles and purpose of the present application. Ordinary technicians in this field can understand that various other specific changes and combinations of embodiments made based on the technical inspiration disclosed in this application without departing from the essence of the present application are still within the scope of protection defined by the claims of the present invention and their equivalent technical solutions.

[0198] At the same time, the technical solutions between the various embodiments can be combined with each other, but they must be based on the ability of ordinary technicians in this field to implement them. When the combination of technical solutions is mutually contradictory or cannot be implemented, it should be deemed that such a combination of technical solutions does not exist and is not within the scope of protection required by this application.

Claims

1. A valve structure simulation method based on finite element analysis, characterized in that: include: designing a geometric model according to structural parameters of the valve structure, and dividing the geometric model into multiple components; defining material parameters of each of the components, and associating the material parameters with each of the components; Using finite element simulation model software to mesh the geometric model to form mesh units, and using all mesh units as simulation analysis objects; Establish a coupled model of the effects of fluid flow, moisture transfer, heat transfer, and mechanical deformation on valve structure; set boundary conditions and initial conditions; Importing the coupling model into the simulation model software to perform simulation analysis of multi-physics field coupling; Based on the simulation results, mesh elements in the geometric model with failure risks are identified and the design is optimized.

2. The valve structure simulation method based on finite element analysis according to claim 1, characterized in that: The coupling model includes: The continuity equation of fluid is used to describe the motion state of fluid flow; The heat conduction equation is used to describe the temperature changes of each component caused by the ambient heat source and the component's own Joule heating; The moisture diffusion equation is used to describe the humidity changes of each component caused by ambient humidity and its own moisture diffusion; The dynamic equations are used to describe the changes in displacement, velocity, acceleration and stress of each component of the valve structure itself.

3. The valve structure simulation method based on finite element analysis according to claim 2, characterized in that: The dynamic equation is described as the relationship between the displacement vector, velocity vector, acceleration vector and external force vector based on a multi-degree-of-freedom system, wherein the coefficients of the displacement vector, velocity vector and acceleration vector are the mass matrix, damping matrix and stiffness matrix respectively.

4. The valve structure simulation method based on finite element analysis according to claim 2, characterized in that: When performing the simulation analysis of multi-physical field coupling, the equations corresponding to the coupling model are solved to obtain simulation results of temperature distribution, humidity distribution, stress distribution and displacement distribution.

5. The valve structure simulation method based on finite element analysis according to claim 4, characterized in that: When solving the equations corresponding to the coupling model: Introducing the energy conservation equation, mass conservation equation, momentum conservation equation, and coupling mechanism of multi-physics field coupling; Discretizing the continuous multi-physics field into a finite number of physics field units corresponding to the grid units; The discretized algebraic equations of each of the physical field elements are solved to obtain an approximate solution.

6. The valve structure simulation method based on finite element analysis according to claim 5, characterized in that: The introduction of the coupling mechanism includes: establishing a relationship between temperature, humidity and stress through the equilibrium equation, displacement equation and thermal stress-strain equation; establishing a relationship between stress change caused by temperature and stress change caused by humidity through thermal stress and wet stress; establishing a relationship for introducing phase change effect by adding latent heat of vaporization term; establishing a relationship between flow fluid pressure, density and flow velocity through continuity equation and mass conservation equation.

7. The valve structure simulation method based on finite element analysis according to claim 5, characterized in that: Solving the discretized algebraic equations specifically includes: By introducing test functions, the partial differential equations of energy conservation equation, mass conservation equation and momentum conservation equation are transformed into integral forms to form weak form expressions respectively; Substitute the temperature field corresponding to each grid unit expressed by the interpolation function into the weak form corresponding to the energy conservation equation, substitute the humidity field into the weak form corresponding to the mass conservation equation, and substitute the displacement field into the weak form corresponding to the momentum conservation equation, to obtain the discretized equations respectively; Using an interpolation function to represent the temperature, humidity and displacement of each grid cell; A coupled simulation is performed to solve the discretized equations.

8. The valve structure simulation method based on finite element analysis according to claim 7, characterized in that: When performing coupled simulation to solve the discretized equations: The momentum conservation equation is discretized using the Newmark method, and then coupled with the discretized energy conservation equation and the discretized mass conservation equation for solution.

9. The valve structure simulation method based on finite element analysis according to claim 7, characterized in that: When performing the coupled simulation solution, a coupling mechanism is introduced, including: the relationship between thermal stress and wet stress, the relationship between temperature and latent heat of vaporization, the fluid mass equation, and the fluid continuity equation.

10. A device for simulating valve structure based on finite element analysis, comprising a memory, a processor, and a computer program stored in the memory and executable in the processor, characterized in that: The computer program includes steps corresponding to the valve structure simulation method based on finite element analysis as described in any one of claims 1 to 9.

11. A readable storage medium storing a terminal program, characterized in that: When the terminal program is executed, the valve structure simulation method based on finite element analysis as described in any one of claims 1 to 9 is performed.

Images