Rotating machinery fluid-solid coupling time domain calculation method and device, electronic equipment and medium

Abstract

The present disclosure relates to the field of fluid-structure coupling calculation of rotating machinery, and provides a fluid-structure coupling time domain calculation method and device for rotating machinery, an electronic device and a medium. The method comprises: establishing a mixed lubrication model and a structure deformation calculation model of the rotating machinery; solving the deformation and the flow field distribution by using an implicit separation method in a specified time step; a time step cycle step: detecting the flow field boundary distribution of the rotating machinery in real time, if no geometric interference occurs, keeping the deformation unchanged in a time period meeting the accuracy requirement, calculating the flow field distribution by using an explicit separation method, updating the deformation and the flow field distribution at a single time point after the time period, and if geometric interference occurs, triggering the implicit separation method to update the deformation distribution; calculating the system stress under the flow field distribution, and updating the motion state by using an implicit multi-step method; and repeating the time step cycle step until the entire load history is completed. The present disclosure improves the calculation efficiency and stability while ensuring the accuracy, and avoids the result oscillation divergence.

Description

Technical Field

[0001] This disclosure relates to the field of fluid-structure interaction calculation technology for rotating machinery, and particularly to a method, apparatus, electronic device, and medium for time-domain calculation of fluid-structure interaction in rotating machinery. Background Technology

[0002] In the field of rotating machinery, numerical solution methods for fluid-structure interaction (FSI) time-domain problems are mainly divided into two categories: strongly coupled methods and segregated methods. Their core objective is to coordinate the interaction between the fluid and the solid through iteration to ensure the convergence and stability of the solution.

[0003] The principle of strongly coupled methods is to combine fluid governing equations, such as the Navier-Stokes equations, with solid governing equations, such as structural dynamics equations, forming a single system of equations for unified solution. The advantages of strongly coupled methods are high accuracy and strong coupling, making them suitable for highly nonlinear or highly dynamic interaction problems. However, strongly coupled methods have extremely high computational complexity, requiring large-scale matrix operations, resulting in high memory and computing power demands and difficulties in engineering implementation. Therefore, this method is generally not used in general engineering calculations.

[0004] The principle of the partitioned method is to solve the fluid and solid domains separately, achieving coupling through iterative transfer of boundary conditions (such as displacement, force, or pressure). Partitioned methods include explicit partitioning, implicit partitioning, and hybrid methods. The explicit partitioned method works as follows: within each time step, the fluid domain is solved independently first, then fluid loads (such as pressure or shear force) are applied to the solid domain, the solid response is solved, and the fluid boundary conditions are updated. A single iteration moves to the next time step. The advantages of explicit partitioning are high computational efficiency and simple implementation. However, explicit partitioning has poor stability, requires extremely small time steps, and is suitable for weakly coupled or quasi-static problems (such as low-speed wind-induced vibration problems). The implicit partitioned method works as follows: within each time step, the fluid and solid solvers are iterated multiple times until the boundary conditions converge (e.g., the displacement residual is less than a threshold). The advantages of implicit partitioning are strong stability and suitability for strongly nonlinear problems (such as structural response problems under explosive impact). However, implicit separation methods are computationally expensive, and the number of iterations increases with the coupling strength. The principle of hybrid methods is to use implicit separation in locally dynamic regions (such as boundary layers) and explicit separation in other regions. Hybrid methods couple regions separately, but require an additional definition of the high-dynamic regions from the other regions. This means that different definitions can affect the solution process. Summary of the Invention

[0005] This disclosure aims to address at least one of the problems existing in the prior art by providing a method, apparatus, electronic device, and medium for time-domain computation of fluid-structure interaction in rotating machinery.

[0006] One aspect of this disclosure provides a time-domain calculation method for fluid-structure interaction (FSI) in rotating machinery, the method comprising:

[0007] A hybrid lubrication model is established to analyze and characterize the performance of rotating machinery under preset working conditions, and a structural deformation calculation model of the rotating machinery is obtained.

[0008] Within a specified time step, the deformation distribution and flow field distribution of the rotating machinery are iteratively solved using the implicit separation method based on the hybrid lubrication model and the structural deformation calculation model.

[0009] Time-step loop steps: The flow field boundary distribution of the rotating machinery is detected in real time. If the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered unchanged within a time period meeting the accuracy requirements. The flow field distribution of the rotating machinery is calculated using an explicit separation method, and at a single moment after the time period meeting the accuracy requirements, the deformation distribution and flow field distribution of the rotating machinery are updated using an implicit separation method. If the flow field boundary distribution shows geometric interference, the deformation distribution of the rotating machinery is updated using an implicit separation method. The system forces of the rotating machinery under the flow field distribution are calculated, and the motion state of the rotating machinery is updated using an implicit multi-step method based on the system forces.

[0010] Repeat the time step cycle until the entire load history of the rotating machinery is completed.

[0011] Optionally, the establishment of a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset operating conditions includes:

[0012] Based on the Reynolds equation for average flow rate considering surface morphology and the JFO cavitation boundary condition, a mass-conserving hydrodynamic lubrication model for the rotating machinery is established.

[0013] The hydrodynamic lubrication model is coupled with the micro-protrusion contact model of the rotating machinery and the rheological properties of the lubricating oil to establish the hybrid lubrication model.

[0014] Optionally, the step of iteratively solving the deformation distribution and flow field distribution of the rotating machinery using the implicit separation method based on the hybrid lubrication model and the structural deformation calculation model includes:

[0015] The elastic deformation is updated to the oil film thickness distribution according to the grid in the hybrid lubrication model;

[0016] The sum of fluid pressure and contact pressure of the rotating machinery is calculated based on the hybrid lubrication model.

[0017] The surface pressure values ​​of the rotating machinery are converted into a grid for calculating the elastic deformation of the rotating machinery;

[0018] The elastic deformation of the rotating machinery is solved based on the structural deformation calculation model.

[0019] The deformation distribution and flow field distribution of the rotating machinery are solved by using a cyclic iterative method.

[0020] Optionally, the system forces of the rotating machinery under the calculated flow field distribution include:

[0021] Calculate the pressure distribution of the rotating machinery based on the hybrid lubrication model;

[0022] Calculate the oil film force and friction force corresponding to the pressure distribution;

[0023] The forces acting on the system are calculated based on the pressure distribution, the oil film force, and the friction force.

[0024] Optionally, updating the motion state of the rotating machinery using an implicit multi-step method, based on the forces acting on the system, includes:

[0025] The motion state of the rotating machinery is updated using a modified extended backward differential formula based on adaptive adjustments of step size and order.

[0026] Optionally, the time period that meets the accuracy requirements is determined in advance based on time-domain calculation experiments of fluid-structure interaction in rotating machinery.

[0027] Another aspect of this disclosure provides a time-domain computational apparatus for fluid-structure interaction of rotating machinery, the apparatus comprising:

[0028] The fluid lubrication and solid deformation modeling module is used to establish a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset working conditions, and to obtain the structural deformation calculation model of the rotating machinery.

[0029] An initialization module is used to iteratively solve the deformation distribution and flow field distribution of the rotating machinery using an implicit separation method within a specified time step, based on the hybrid lubrication model and the structural deformation calculation model.

[0030] The time-step calculation module is used to execute the time-step loop steps: Real-time detection of the flow field boundary distribution of the rotating machinery; if the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered unchanged within a time period meeting the accuracy requirements; the flow field distribution of the rotating machinery is calculated using an explicit separation method; and at a single moment after the time period meeting the accuracy requirements, the deformation distribution and flow field distribution of the rotating machinery are updated using an implicit separation method; if the flow field boundary distribution shows geometric interference, the implicit separation method is triggered to update the deformation distribution of the rotating machinery; the system forces of the rotating machinery under the flow field distribution are calculated; and the motion state of the rotating machinery is updated using an implicit multi-step method based on the system forces.

[0031] The time-domain calculation module is used to repeatedly trigger the time-step calculation module to execute the time-step loop steps until the entire load history of the rotating machinery is completed.

[0032] Another aspect of this disclosure provides an electronic device comprising:

[0033] At least one processor; and,

[0034] A memory that is communicatively connected to at least one processor; wherein,

[0035] The memory stores instructions that can be executed by at least one processor, which enables the at least one processor to perform the rotating machinery fluid-structure interaction time-domain calculation method described above.

[0036] Another aspect of this disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the aforementioned method for time-domain computation of fluid-structure interaction in rotating machinery.

[0037] Another aspect of this disclosure provides a computer program product, including a computer program that, when executed by a processor, implements the aforementioned method for time-domain computation of fluid-structure interaction in rotating machinery.

[0038] Compared with existing technologies, this disclosure significantly reduces computation time while maintaining accuracy by dynamically switching between implicit and explicit separation methods, compared with existing implicit separation methods; compared with existing explicit separation methods, it effectively improves computational stability and avoids oscillation and divergence of computational results by dynamically switching between implicit and explicit separation methods; compared with existing hybrid methods, it does not require the introduction of additional control equations and achieves complementary advantages of implicit and explicit separation methods through dynamic coupling strategy. Attached Figure Description

[0039] One or more embodiments are illustrated by way of example with the corresponding pictures in the accompanying drawings. These illustrations do not constitute a limitation on the embodiments. Elements with the same reference numerals in the drawings are denoted as similar elements. Unless otherwise stated, the figures in the drawings are not to be limited by scale.

[0040] Figure 1 A flowchart of a time-domain calculation method for fluid-structure interaction in rotating machinery, provided as an embodiment of this disclosure;

[0041] Figure 2 A flowchart illustrating a time-domain calculation method for fluid-structure interaction in rotating machinery, provided as another embodiment of this disclosure;

[0042] Figure 3 A schematic diagram of the structure of a time-domain computing device for fluid-structure interaction of rotating machinery provided in another embodiment of this disclosure;

[0043] Figure 4 A schematic diagram of the structure of an electronic device provided in another embodiment of this disclosure. Detailed Implementation

[0044] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the various embodiments of this disclosure will be described in detail below with reference to the accompanying drawings. However, those skilled in the art will understand that many technical details have been provided in the various embodiments of this disclosure to facilitate a better understanding of the disclosure. However, the technical solutions claimed in this disclosure can be implemented even without these technical details and with various variations and modifications based on the following embodiments. The division of the various embodiments below is for ease of description and should not constitute any limitation on the specific implementation of this disclosure. The various embodiments can be combined with and referenced by each other without contradiction.

[0045] One embodiment of this disclosure relates to a time-domain computation method for fluid-structure interaction in rotating machinery, the process of which is as follows: Figure 1 As shown, it includes steps S110 to S140.

[0046] Step S110: Establish a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset working conditions, and obtain a structural deformation calculation model of the rotating machinery.

[0047] For example, rotating machinery here can be dynamically loaded sliding bearings. Under most common operating conditions, dynamically loaded sliding bearings are designed to operate in a fluid lubrication state, with no contact between the bearing and journal surfaces to avoid bearing wear. However, in special operating conditions, dynamically loaded sliding bearings may operate in a mixed lubrication or boundary lubrication state due to factors such as low speed, low lubricant viscosity, or heavy load. In this case, the minimum oil film thickness of the bearing will be on the order of magnitude of the bearing surface roughness peak height, resulting in micro-protrusion contact between the working surfaces, and these surface micro-protrusions will significantly affect lubrication analysis. For the special condition of bearing start-stop, the engine main bearing may begin operating before the oil supply pressure is established. In this case, the bearing will be in a lean lubrication state, which places higher demands on the lubrication model. Furthermore, under heavy load conditions, the local pressure of the bearing may be very high, and the resulting bearing deformation and changes in lubricant viscosity will have a non-negligible impact on lubrication analysis.

[0048] For example, in step S110, a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset operating conditions is established, including: establishing a mass-conserving hydrodynamic lubrication model for rotating machinery based on the average flow Reynolds equation considering surface morphology and JFO cavitation boundary conditions; and coupling the hydrodynamic lubrication model with the micro-protrusion contact model of rotating machinery and the rheological properties of lubricating oil to establish a hybrid lubrication model.

[0049] Specifically, the hydrodynamic lubrication model here can be expressed as:

[0050]

[0051] Where x represents the position in the x-direction of motion, y represents the position in the perpendicular direction of motion, and Θ represents the cavitation factor. ρ represents the density of the lubricating oil. avg φ represents the average density of the mixture within the gap. px φ represents the pressure-flow factor in the x-direction. py η represents the pressure-flow factor in the y-direction. h represents the height of the gap between the two surfaces, i.e., the oil film thickness. η represents the viscosity of the lubricant. p represents the average oil film pressure within the micro-element. U1 and U2 represent the movement velocities of the upper and lower surfaces, respectively. This represents the average gap height considering surface micro-protrusions. σ represents the combined roughness of the two surfaces. σ1 and σ2 represent the standard deviations of the heights of the micro-protrusions on the two rough surfaces, respectively. φ s This represents the shear flow factor. t represents time.

[0052] The micro-convexity contact model can be represented as: p ac =KE*F 5 / 2 (H). Among them, p ac Indicates the solid contact pressure. E* The combined elastic modulus of two contacting surfaces is determined based on the elastic moduli E1 and E2 of the materials of the two contacting surfaces and their Poisson's ratios υ1 and v2. The calculation method is as follows: F 5 / 2 (H) is a function of the dimensionless film thickness H, and H = h / σ. K represents a function of the topographic parameters of the two rough surfaces, and its expression is: β c and ζ c All of these are morphological parameters of rough surfaces.

[0053] The rheological properties of lubricating oil can be expressed using the viscosity-temperature formula and the viscosity-pressure formula. The viscosity-temperature formula can be expressed as: η0 is the reference viscosity. T a and T b These are reference temperatures, in °C. η1 represents the viscosity at temperature T.

[0054] The viscosity-pressure formula can be expressed as: Where z0 is the dimensionless viscosity coefficient.

[0055] By coupling the above hydrodynamic lubrication model with the micro-convexity contact model and the formulas related to the rheological properties of lubricating oil, a hybrid lubrication model can be obtained.

[0056] Step S120: Within a specified time step, based on the hybrid lubrication model and the structural deformation calculation model, the deformation distribution and flow field distribution of the rotating machinery are iteratively solved using the implicit separation method.

[0057] Specifically, the specified time step can be set according to actual needs.

[0058] For example, in step S120, based on the hybrid lubrication model and the structural deformation calculation model, the deformation distribution and flow field distribution of the rotating machinery are iteratively solved using the implicit separation method, including: updating the elastic deformation to the oil film thickness distribution according to the grid in the hybrid lubrication model; solving for the sum of fluid pressure and contact pressure of the rotating machinery according to the hybrid lubrication model; converting the surface pressure values ​​of the rotating machinery to the grid for calculating the elastic deformation of the rotating machinery; solving for the elastic deformation of the rotating machinery according to the structural deformation calculation model; and solving for the deformation distribution and flow field distribution of the rotating machinery using a cyclic iterative method.

[0059] Specifically, the fluid pressure can be calculated by discretizing the two-dimensional Reynolds equations using the finite difference method or the finite volume method, obtaining the corresponding discrete equation set. Then, given boundary conditions are applied to numerically solve the discrete equation set to obtain the fluid pressure. The two-dimensional Reynolds equations can be expressed as: Fluid pressure can be characterized by the parameter p, which represents the average oil film pressure within a micro-element.

[0060] Contact pressure can be expressed by the parameter p, which represents the contact pressure of a solid. ac To represent, i.e., p ac =KE*F 5 / 2 (H).

[0061] The formula for calculating elastic deformation can be expressed as: Among them, h e This represents the radial deformation of node (l,k) in rotating machinery. P(l,k) represents the pressure acting on node (l,k), including oil film pressure and solid contact pressure. dA is the area of ​​a single mesh. is the value of node (i,j,l,k) in the deformation influence coefficient matrix D. Its physical meaning is the radial deformation value produced at node (i,j) by a unit force at node (l,k). Here, l and i are the circumferential numbers of the nodes in the rotating machinery, k and j are the axial numbers of the nodes in the rotating machinery, and m and n represent the maximum values ​​of l and k, respectively.

[0062] When solving for the deformation and flow field distributions of rotating machinery using iterative methods, any existing iterative method can be employed. For example, a specific iterative method could be the Etkin-accelerated low-relaxation iterative method.

[0063] Step S130, Time Step Loop Step: Real-time detection of the flow field boundary distribution of the rotating machinery. If the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered to remain unchanged within a time period that meets the accuracy requirements. The flow field distribution of the rotating machinery is calculated using the explicit separation method, and the deformation distribution and flow field distribution of the rotating machinery are updated using the implicit separation method at a single moment after the time period that meets the accuracy requirements. If the flow field boundary distribution shows geometric interference, the implicit separation method is triggered to update the deformation distribution of the rotating machinery. The system forces of the rotating machinery under the flow field distribution are calculated, and the motion state of the rotating machinery is updated using the implicit multi-step method in combination with the system forces.

[0064] Specifically, the flow field boundary here refers to the physical boundary of the rotating machinery. For example, when the rotating machinery is a sliding bearing, the flow field boundary refers to the physical clearance between the journal and the bearing. If this physical clearance is negative, it indicates that geometric interference has occurred in the sliding bearing. Real-time detection of the flow field boundary distribution of the rotating machinery can be achieved using existing numerical simulation methods. To determine whether geometric interference has occurred, the clearance distribution of the rotating machinery can first be calculated using existing clearance distribution calculation methods, and then the presence of geometric interference can be determined based on this clearance distribution.

[0065] The specific process of calculating the flow field distribution of rotating machinery using the explicit separation method is actually the solution process of the mixed lubrication model. Specifically, the explicit separation method in the existing technology can be used for the solution, which will not be elaborated here.

[0066] When updating the deformation distribution and flow field distribution of the rotating machinery using the implicit separation method in step S130, the specific steps are the same as those in step S120, which iteratively solves the deformation distribution and flow field distribution of the rotating machinery using the implicit separation method. Therefore, they will not be repeated here.

[0067] For example, in step S130, calculating the system forces of the rotating machinery under the flow field distribution includes: calculating the pressure distribution of the rotating machinery according to the mixed lubrication model; calculating the oil film force and friction force corresponding to the pressure distribution; and calculating the system forces according to the pressure distribution, oil film force, and friction force.

[0068] Specifically, the formula for calculating oil film force can be expressed as:

[0069]

[0070] Among them, F x F y These represent the oil film forces in the x and y directions of the rotating machinery, respectively. L is the axial length of the rotating machinery, and R is the radius of the rotating machinery. θ represents an intermediate variable and... v represents the circumferential coordinate of the bearing in the rotating machinery, R b This indicates the bearing radius of rotating machinery.

[0071] Friction can be represented by solid contact friction, and its expression can be:

[0072]

[0073] Among them, F ac This represents the solid contact friction force. k0 is the correction factor for the Coulomb coefficient of friction, ξ0 is the ratio of classical Coulomb kinetic friction to the maximum static friction, α0 is the coefficient determining how quickly the friction coefficient decays from static friction to kinetic friction, σ0 is the coefficient determining how quickly k0 changes from -1 to +1, and μ... c is the Coulomb friction coefficient in solid contact, and u is the relative linear velocity between the two surfaces.

[0074] The system force is used to indicate the overall force situation of rotating machinery. After obtaining the oil film force and friction force, the oil film force and friction force can be vector-superimposed to obtain the system force of the rotating machinery.

[0075] For example, in step S130, the motion state of the rotating machinery is updated using an implicit multi-step method in combination with the system forces, including: updating the motion state of the rotating machinery using a modified extended backward differential formula (MEBDF) with adaptive adjustment based on step size and order.

[0076] Specifically, the general form of a k-th order MEBDF is:

[0077]

[0078] Here, "order k" means using k known motion states (n to n+k-1) of the rotating machinery to deduce its motion state at the (n+k)th time step. n+k Let y represent the motion state of the rotating machinery at time step n+k. n+k Let y represent dy / dt at time step n+k. n+j Let y represent the motion state of the rotating machinery at time step n+j. Indicated by, according to The calculated dy / dt. Indicates according to The calculated dy / dt. Where, and The motion states y of the rotating machinery at time step n+k and time step n+k+1 are predicted and obtained according to the standard backward differential formula (BDF) method. This represents the correction coefficient at time step k+1. Δt is the integration time step size. All are correction factors.

[0079] At each time step, y n+j (j=0,…,k-1) is known, and is first predicted using the standard backward differential formula (BDF) method. and The value of f is then obtained according to the definition of function f. and The value of y is finally obtained through iteration. n+k The value is corrected.

[0080] The function f can be defined as follows: Let the motion state of the rotating machinery be... Then we have: Where u represents the displacement of the rotating machinery, u represents the velocity of the rotating machinery, and T represents the transpose. Let t represent the acceleration of the rotating machinery, and let V represent time. The function V represents the process of solving the hybrid lubrication model based on the displacement and velocity of the rotating machinery to obtain the oil film force and contact force, and combining this with the external forces currently acting on the system to calculate the acceleration of the rotating machinery. At the initial time t0, the initial motion state of the rotating machinery is y0 = y(t0).

[0081] For example, the time period that meets the accuracy requirements is determined in advance based on the time-domain calculation experiment of fluid-structure interaction of rotating machinery.

[0082] Specifically, when conducting time-domain calculation experiments on fluid-structure interaction (FSI) of rotating machinery, the time-domain calculation method for FSI provided in this embodiment can be used within a specified time history. At the same time, the fully implicit separation method can be used to perform time-domain calculations of FSI of rotating machinery within the same time history. The calculation results of the two methods are compared. Under the premise of meeting the accuracy requirements, the largest possible time step in the time history is selected as the time period that meets the accuracy requirements, thereby taking into account both computational efficiency and accuracy requirements, and balancing computational speed and computational accuracy.

[0083] Step S140: Repeat the time step cycle until the entire load history of the rotating machinery is completed.

[0084] The time-domain computation method for fluid-structure interaction of rotating machinery provided in this disclosure, compared with the prior art, significantly shortens the computation time while ensuring accuracy by dynamically switching between implicit and explicit separation methods, compared with the existing implicit separation method; compared with the existing explicit separation method, it uses implicit separation to update the deformation distribution and flow field distribution of the rotating machinery for correction when geometric interference occurs, effectively improving computational stability and avoiding oscillation and divergence of computational results; compared with the existing hybrid method, it does not require the introduction of additional control equations, and achieves the complementary advantages of implicit and explicit separation methods through dynamic coupling strategy.

[0085] To enable those skilled in the art to better understand the above embodiments, a specific example is provided below.

[0086] Combined Figure 2 A time-domain calculation method for fluid-structure interaction in rotating machinery includes the following steps.

[0087] Set time t = 0, initialize calculation parameters and dynamic variables: establish a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset working conditions, and obtain the structural deformation calculation model of the rotating machinery. Set the output time step Δt, and define the dynamic variables of the rotating machinery. Where u represents the displacement of the rotating machinery. The speed of the rotating machinery is indicated by T, which stands for transpose.

[0088] The first two steps use the implicit separation method to calculate the deformation and flow field distribution: within the first two output time steps, i.e., when t < 2Δt, the implicit separation method is used to iteratively solve the deformation distribution and flow field distribution of the rotating machinery.

[0089] Time step loop: Specifically includes the following steps 1 to 4.

[0090] Step 1: Calculate the forces acting on the system under the flow field distribution. Using the implicit multi-step method, complete the system dynamics integration within the current Δt and update the dynamic variables y and time t.

[0091] Step 2: Determine whether the set output time has been reached. If yes, determine whether the termination time has been reached. If no, based on the current dynamic variable y, update the deformation distribution of the rotating machinery using the implicit separation method, record the current output time, and update the next output time.

[0092] Step 3: If the set output time has not been reached, then if t < 2Δt is not satisfied in the next Δt, keep the deformation of the rotating machinery unchanged, and use the explicit separation method to calculate the time domain change of the fluid field of the rotating machinery, i.e., the flow field distribution.

[0093] Step 4: Detect the physical space of the rotating machinery considering deformation. If the detection result shows that geometric interference occurs, the implicit separation method is triggered to update the deformation distribution of the rotating machinery and solve the flow field distribution of the rotating machinery. Return to Step 1 to perform iterative calculations until the end.

[0094] Repeat steps 1 through 4 until the entire load history of the rotating machinery is completed.

[0095] Another embodiment of this disclosure relates to a time-domain computing device for fluid-structure interaction in rotating machinery, such as... Figure 3 As shown, it includes a fluid lubrication and solid deformation modeling module 310, an initialization module 320, an intra-time step calculation module 330, and a time domain calculation module 340.

[0096] The fluid lubrication and solid deformation modeling module 310 is used to establish a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset working conditions, and to obtain a structural deformation calculation model of the rotating machinery.

[0097] The initialization module 320 is used to iteratively solve the deformation distribution and flow field distribution of rotating machinery using the implicit separation method within a specified time step, based on the hybrid lubrication model and the structural deformation calculation model.

[0098] The time-step calculation module 330 is used to execute the time-step loop steps: real-time detection of the flow field boundary distribution of the rotating machinery; if the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered to remain unchanged within a time period that meets the accuracy requirements; the flow field distribution of the rotating machinery is calculated using the explicit separation method; and at a single moment after the time period that meets the accuracy requirements, the deformation distribution and flow field distribution of the rotating machinery are updated using the implicit separation method; if the flow field boundary distribution shows geometric interference, the implicit separation method is triggered to update the deformation distribution of the rotating machinery; the system forces of the rotating machinery under the flow field distribution are calculated; and the motion state of the rotating machinery is updated using the implicit multi-step method based on the system forces.

[0099] The time-domain calculation module 340 is used to repeatedly trigger the time-step calculation module to execute the time-step loop steps until the entire load history of the rotating machinery is completed.

[0100] The specific implementation method of the rotating machinery fluid-structure interaction time-domain calculation device provided in this disclosure can be found in the rotating machinery fluid-structure interaction time-domain calculation method provided in this disclosure, and will not be repeated here.

[0101] The rotating machinery fluid-structure interaction time-domain computing device provided in this disclosure significantly shortens the computation time while ensuring accuracy by dynamically switching between implicit and explicit separation methods. When geometric interference occurs, the implicit separation method is used to update the deformation distribution and flow field distribution of the rotating machinery for correction, which effectively improves the computational stability and avoids oscillation and divergence of the computational results. It also eliminates the need to introduce additional control equations. The dynamic coupling strategy realizes the complementary advantages of implicit and explicit separation methods.

[0102] Another embodiment of this disclosure relates to an electronic device, such as Figure 4 As shown, it includes:

[0103] At least one processor 401; and,

[0104] Memory 402 is communicatively connected to at least one processor 401; wherein,

[0105] The memory 402 stores instructions that can be executed by at least one processor 401, which enables the at least one processor 401 to perform the rotating machinery fluid-structure interaction time-domain calculation method described in the above embodiments.

[0106] The memory and processor are connected via a bus, which can include any number of interconnecting buses and bridges, connecting various circuits of one or more processors and memories. The bus can also connect various other circuits, such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and will not be described further herein. The bus interface provides an interface between the bus and the transceiver. The transceiver can be a single element or multiple elements, such as multiple receivers and transmitters, providing a unit for communicating with various other devices over a transmission medium. Data processed by the processor is transmitted over the wireless medium via an antenna, which further receives data and transmits it to the processor.

[0107] The processor manages the bus and general processing, and also provides various functions, including timing, peripheral interfaces, voltage regulation, power management, and other control functions. Memory is used to store data used by the processor during operation.

[0108] Another embodiment of this disclosure relates to a computer-readable storage medium storing a computer program that, when executed by a processor, implements the time-domain calculation method for fluid-structure interaction of rotating machinery as described in the above embodiments.

[0109] That is, those skilled in the art will understand that all or part of the steps in the methods described in the above embodiments can be implemented by a program instructing related hardware. This program is stored in a storage medium and includes several instructions to cause a device (which may be a microcontroller, chip, etc.) or processor to execute all or part of the steps of the methods described in the various embodiments of this disclosure. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0110] Another embodiment of this disclosure relates to a computer program product, including a computer program that, when executed by a processor, implements the time-domain calculation method for fluid-structure interaction of rotating machinery as described in the above embodiments.

[0111] Those skilled in the art will understand that the above embodiments are specific implementations of this disclosure, and in practical applications, various changes can be made to them in form and detail without departing from the spirit and scope of this disclosure.

Claims

1. A time-domain calculation method for fluid-structure interaction in rotating machinery, characterized in that, The time-domain calculation method for fluid-structure interaction of rotating machinery includes: A hybrid lubrication model is established to analyze and characterize the performance of rotating machinery under preset working conditions, and a structural deformation calculation model of the rotating machinery is obtained. Within a specified time step, the deformation distribution and flow field distribution of the rotating machinery are iteratively solved using the implicit separation method based on the hybrid lubrication model and the structural deformation calculation model. Time-step loop steps: The flow field boundary distribution of the rotating machinery is detected in real time. If the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered unchanged within a time period meeting the accuracy requirements. The flow field distribution of the rotating machinery is calculated using an explicit separation method, and at a single moment after the time period meeting the accuracy requirements, the deformation distribution and flow field distribution of the rotating machinery are updated using an implicit separation method. If the flow field boundary distribution shows geometric interference, the deformation distribution of the rotating machinery is updated using an implicit separation method. The system forces of the rotating machinery under the flow field distribution are calculated, and the motion state of the rotating machinery is updated using an implicit multi-step method based on the system forces. Repeat the time step cycle until the entire load history of the rotating machinery is completed.

2. The time-domain calculation method for fluid-structure interaction in rotating machinery according to claim 1, characterized in that, The establishment of a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset operating conditions includes: Based on the Reynolds equation for average flow rate considering surface morphology and the JFO cavitation boundary condition, a mass-conserving hydrodynamic lubrication model for the rotating machinery is established. The hydrodynamic lubrication model is coupled with the micro-protrusion contact model of the rotating machinery and the rheological properties of the lubricating oil to establish the hybrid lubrication model.

3. The time-domain calculation method for fluid-structure interaction in rotating machinery according to claim 1, characterized in that, The step of iteratively solving the deformation distribution and flow field distribution of the rotating machinery using the implicit separation method based on the hybrid lubrication model and the structural deformation calculation model includes: The elastic deformation is updated to the oil film thickness distribution according to the grid in the hybrid lubrication model; The sum of fluid pressure and contact pressure of the rotating machinery is calculated based on the hybrid lubrication model. The surface pressure values ​​of the rotating machinery are converted into a grid for calculating the elastic deformation of the rotating machinery; The elastic deformation of the rotating machinery is solved based on the structural deformation calculation model. The deformation distribution and flow field distribution of the rotating machinery are solved by using a cyclic iterative method.

4. The time-domain calculation method for fluid-structure interaction in rotating machinery according to claim 1, characterized in that, The forces acting on the rotating machinery system under the calculated flow field distribution include: Calculate the pressure distribution of the rotating machinery based on the hybrid lubrication model; Calculate the oil film force and friction force corresponding to the pressure distribution; The forces acting on the system are calculated based on the pressure distribution, the oil film force, and the friction force.

5. The time-domain calculation method for fluid-structure interaction in rotating machinery according to claim 4, characterized in that, The method of updating the motion state of the rotating machinery by combining the forces acting on the system and using an implicit multi-step method includes: The motion state of the rotating machinery is updated using a modified extended backward differential formula based on adaptive adjustments of step size and order.

6. The time-domain calculation method for fluid-structure interaction of rotating machinery according to any one of claims 1 to 5, characterized in that, The time period that meets the accuracy requirements is determined in advance based on time-domain calculation experiments of fluid-structure interaction in rotating machinery.

7. A time-domain computational device for fluid-structure interaction of rotating machinery, characterized in that, The rotating machinery fluid-structure interaction time-domain computing device includes: The fluid lubrication and solid deformation modeling module is used to establish a hybrid lubrication model for analyzing and characterizing the performance of rotating machinery under preset working conditions, and to obtain the structural deformation calculation model of the rotating machinery. An initialization module is used to iteratively solve the deformation distribution and flow field distribution of the rotating machinery using an implicit separation method within a specified time step, based on the hybrid lubrication model and the structural deformation calculation model. The time-step calculation module is used to execute the time-step loop steps: Real-time detection of the flow field boundary distribution of the rotating machinery; if the flow field boundary distribution shows no geometric interference, the deformation of the rotating machinery is considered unchanged within a time period meeting the accuracy requirements; the flow field distribution of the rotating machinery is calculated using an explicit separation method; and at a single moment after the time period meeting the accuracy requirements, the deformation distribution and flow field distribution of the rotating machinery are updated using an implicit separation method; if the flow field boundary distribution shows geometric interference, the implicit separation method is triggered to update the deformation distribution of the rotating machinery; the system forces of the rotating machinery under the flow field distribution are calculated; and the motion state of the rotating machinery is updated using an implicit multi-step method based on the system forces. The time-domain calculation module is used to repeatedly trigger the calculation module within the time step to execute the time-step loop steps until the entire load history of the rotating machinery is completed.

8. An electronic device, characterized in that, include: At least one processor; as well as, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the rotating machinery fluid-structure interaction time-domain calculation method according to any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the time-domain calculation method for fluid-structure interaction of rotating machinery as described in any one of claims 1 to 6.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the time-domain calculation method for fluid-structure interaction of rotating machinery as described in any one of claims 1 to 6.

Images