Method and device for analyzing honeycomb seal dynamic characteristics, storage medium and electronic device
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- SHENYANG BLOWER WORKS GROUP CORP
- Filing Date
- 2026-03-09
- Publication Date
- 2026-06-05
Smart Images
Figure CN122154548A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of sealing structure performance analysis technology, and in particular to a method and apparatus for analyzing the dynamic characteristics of a honeycomb seal, a storage medium, and an electronic device. Background Technology
[0002] In high-pressure (40 MPa) turbine units, the gas swirl-assisted honeycomb seal structure has become a key component to ensure stable rotor operation. The dynamic characteristic coefficient of the honeycomb seal with swirl-assisted structure is an important parameter for evaluating the stability of the rotor system. Therefore, it is crucial to accurately solve the dynamic characteristic coefficient of the seal.
[0003] However, the acquisition of dynamic characteristic coefficients (stiffness, damping) currently faces three major contradictions: (1) Insufficient accuracy: the commonly used control volume model (Bulk-Flow) acquisition method is fast but inaccurate; the commonly used Computational Fluid Dynamics (CFD) model acquisition method is accurate, but lacks experimental verification under high pressure conditions; (2) Too slow speed: traditional CFD modeling is complex and computationally intensive, and a single analysis under a single condition often takes several days, which cannot meet the needs of rapid iteration in engineering; (3) Too high cost: high pressure sealing test benches are scarce, and the test cycle is long and expensive. Therefore, there is an urgent need for a method that can shorten the calculation time of the dynamic characteristic coefficients of honeycomb seals from the 'day' level to the 'minute' level while ensuring the accuracy of high pressure conditions, so as to meet the engineering needs of rapid deployment. Summary of the Invention
[0004] In view of this, the present invention provides a method and apparatus for analyzing the dynamic characteristics of a honeycomb seal, a storage medium, and an electronic device. The main purpose is to solve the problem of how to improve the accuracy and speed of obtaining the dynamic characteristic coefficients of a gas-swirled honeycomb seal structure without increasing costs and while ensuring accuracy under high-pressure conditions.
[0005] According to one aspect of the present invention, a method for analyzing the dynamic characteristics of a honeycomb seal is provided, comprising: The honeycomb sealing fluid domain in the honeycomb sealing structure with gas spool is parametrically modeled to obtain the honeycomb sealing parametric fluid domain model. The efficiency of CFD analysis is optimized by combining the aforementioned parameterized fluid domain model of the honeycomb seal to obtain the optimized CFD analysis model; the efficiency optimization includes simplification of the single-frequency eddy model and simplification of the dynamic characteristic coefficients; A large number of computational samples were generated based on parameter combinations, and the sealing dynamic characteristics of each sample were calculated using an optimized CFD analysis model. By fitting and analyzing the sensitivity of the sealing dynamic characteristics results to different parameters, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb sealing structure is obtained.
[0006] Furthermore, the parametric modeling of the honeycomb sealing fluid domain in the gas-puffed honeycomb sealing structure to obtain a parametric fluid domain model for the honeycomb sealing includes: Obtain the pre-established complete swirl block structure model and analyze the influence of the complete swirl block structure model on the pre-swirl of the airflow inlet. Based on the aforementioned influence law, the pre-spinning parameter is used to replace the solid structure to indirectly reflect the inlet gas spinning effect; A parametric interface was built based on the modeling platform, and the key dimensions of the honeycomb seal were set as editable parameters. The pre-swirl parameter, which indirectly reflects the inlet gas swirl effect, is adapted in conjunction with the editable key dimensions of the honeycomb seal. This allows the model to be automatically updated and the corresponding inlet gas swirl effect to be loaded synchronously when the user edits the key dimensions of the honeycomb seal, ultimately resulting in a parameterized fluid domain model of the honeycomb seal that incorporates the inlet gas swirl effect.
[0007] Furthermore, prior to the efficiency optimization process of combining the parameterized fluid domain model of the cellular seal for CFD analysis, the method further includes: Verify the adaptability of the proposed honeycomb seal parameterized fluid domain model to single-frequency eddy currents; After successful verification, the pre-swirl parameters were used as inlet boundary conditions and incorporated into the flow field calculation of single-frequency eddy currents to ensure that the inlet gas swirl effect was accurately captured under single-frequency conditions.
[0008] Furthermore, the efficiency optimization process for CFD analysis based on the parameterized fluid domain model of the cellular seal, resulting in an optimized CFD analysis model, includes: The multi-frequency eddy model in the CFD analysis process is simplified into a single-frequency eddy model to reduce unnecessary frequency domain calculations; the single-frequency eddy model adopts the rotor single-frequency elliptical eddy equation.
[0009] Furthermore, the efficiency optimization process for CFD analysis based on the parameterized fluid domain model of the cellular seal, to obtain the optimized CFD analysis model, further includes: When the rotor is subjected to minute displacement and velocity disturbances, the relationship between the airflow force on the rotor and the disturbance displacement and velocity is linearized based on the single-frequency eddy current model as follows:
[0010] in, F x This is the excitation force in the x-direction of the airflow;F y This is the excitation force in the y-direction of the airflow; K Principal stiffness; k For cross stiffness; C Main damping; c For cross damping; X Displacement of rotor shaft center x direction; Y Displacement of rotor shaft center y direction; Rotor shaft disturbance speed x direction; Rotor shaft disturbance speed y direction; Performing a Fast Fourier Transform on the linearized relationship yields the following equations: The relationship between the change in sealing gas flow force in the frequency domain and the sealing dynamic characteristic coefficient and the displacement of small disturbances is:
[0011] in, ; It is the vortex angular velocity; K xx , K xy , K yx and K yy These are the four stiffness coefficients in the single-frequency eddy coulomb model; C xx , C xy , C yx and C yy These are the four damping coefficients in the single-frequency eddy current model; D x , D y ) and (Δ F x Δ F y These are the time-domain signals of the rotor eddy displacement and the frequency-domain signals corresponding to the time-domain values of the aerodynamic force changes, respectively. Utilizing the symmetry of the dynamic characteristic coefficients during rotor whirling around the axis, the eight dynamic characteristic coefficients in the single-frequency whirling model are analyzed. K xx , K xy , K yx , K yy , C xx ,C xy , C yx , C yy After simplification, the specific formula is as follows:
[0012]
[0013]
[0014]
[0015] The simplified dynamic characteristic coefficients are then substituted into the single-frequency eddy current formula for calculation. and The rotor motion states at two different times are shown below: when hour, ; when hour, ; Where a and b are the major and minor axes of the ellipse, respectively; , , and These are the optimized four sealing dynamic characteristic coefficients.
[0016] Furthermore, the step of analyzing the sensitivity of the sealing dynamic characteristic results to different parameters through fitting, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, includes: The influence of univariate parameters on the dynamic characteristics of honeycomb sealing was analyzed using the correction coefficient method, and the results of the univariate influence analysis were obtained. Regression analysis was used to analyze the influence of multivariate parameters on the dynamic characteristics of honeycomb sealing, and the results of the multivariate influence analysis were obtained. By integrating the results of the univariate influence analysis and the multivariate influence analysis, a unified quantitative relationship between parameters and sealing dynamic characteristics is formed.
[0017] Furthermore, the method also includes: Store the correspondence between the input parameters of all samples and the sealing dynamic characteristic coefficients to form a structured database; Enter the fitted polynomial coefficients and correction coefficients, and establish the connection between the quantization relationship and the database; A visual input interface is adopted so that engineers can input sealing structure parameters and operating condition parameters through the visual input interface; The built-in calling logic invokes the quantization relationship, and the cell sealing dynamic characteristic coefficient is calculated and output based on the quantization relationship.
[0018] According to another aspect of the present invention, an analytical apparatus for the dynamic characteristics of a honeycomb seal is provided, comprising: The parametric modeling module is used to perform parametric modeling of the honeycomb sealing fluid domain in a honeycomb sealing structure with gas spool, and obtain a parametric fluid domain model of the honeycomb sealing. The CFD optimization module is used to perform efficiency optimization processing on the CFD analysis in conjunction with the parameterized fluid domain model of the honeycomb seal, so as to obtain an optimized CFD analysis model; the efficiency optimization processing includes simplification processing of the single-frequency eddy model and simplification processing of the dynamic characteristic coefficients; The data preparation module is used to generate a large number of computational samples based on parameter combinations, and to calculate the sealing dynamic characteristics of each sample using an optimized CFD analysis model. The analysis and quantification module is used to analyze the sensitivity of the sealing dynamic characteristics results to different parameters through fitting analysis, and to establish a quantitative relationship between the parameters and the sealing dynamic characteristics, so as to obtain the influence law of different types of parameters on the stiffness and damping in the honeycomb sealing structure.
[0019] Furthermore, the parametric modeling module is also used for: Obtain the pre-established complete swirl block structure model and analyze the influence of the complete swirl block structure model on the pre-swirl of the airflow inlet. Based on the aforementioned influence law, the pre-spinning parameter is used to replace the solid structure to indirectly reflect the inlet gas spinning effect; A parametric interface was built based on the modeling platform, and the key dimensions of the honeycomb seal were set as editable parameters. The pre-swirl parameter, which indirectly reflects the inlet gas swirl effect, is adapted in conjunction with the editable key dimensions of the honeycomb seal. This allows the model to be automatically updated and the corresponding inlet gas swirl effect to be loaded synchronously when the user edits the key dimensions of the honeycomb seal, ultimately resulting in a parameterized fluid domain model of the honeycomb seal that incorporates the inlet gas swirl effect.
[0020] Furthermore, the device also includes a verification module, which is used for: Verify the adaptability of the proposed honeycomb seal parameterized fluid domain model to single-frequency eddy currents; After successful verification, the pre-swirl parameters were used as inlet boundary conditions and incorporated into the flow field calculation of single-frequency eddy currents to ensure that the inlet gas swirl effect was accurately captured under single-frequency conditions.
[0021] Furthermore, the CFD optimization module is also used for: The multi-frequency eddy model in the CFD analysis process is simplified into a single-frequency eddy model to reduce unnecessary frequency domain calculations; the single-frequency eddy model adopts the rotor single-frequency elliptical eddy equation.
[0022] Furthermore, the CFD optimization module is also used for: When the rotor is subjected to minute displacement and velocity disturbances, the relationship between the airflow force on the rotor and the disturbance displacement and velocity is linearized based on the single-frequency eddy current model as follows:
[0023] in, F x This is the excitation force in the x-direction of the airflow; F y This is the excitation force in the y-direction of the airflow; K Principal stiffness; k For cross stiffness; C Main damping; c For cross damping; X Displacement of rotor shaft center x direction; Y Displacement of rotor shaft center y direction; Rotor shaft disturbance speed x direction; Rotor shaft disturbance speed y direction; Performing a Fast Fourier Transform on the linearized relationship yields the following equations: The relationship between the change in sealing gas flow force in the frequency domain and the sealing dynamic characteristic coefficient and the displacement of small disturbances is:
[0024] in, ; It is the vortex angular velocity; K xx , K xy , K yx and K yy These are the four stiffness coefficients in the single-frequency eddy coulomb model; C xx , C xy , C yx and C yy These are the four damping coefficients in the single-frequency eddy current model; D x , D y ) and (Δ F x Δ F y These are the time-domain signals of the rotor eddy displacement and the frequency-domain signals corresponding to the time-domain values of the aerodynamic force changes, respectively. Utilizing the symmetry of the dynamic characteristic coefficients during rotor whirling around the axis, the eight dynamic characteristic coefficients in the single-frequency whirling model are analyzed. K xx , K xy , K yx , K yy , C xx , C xy , C yx , C yy After simplification, the specific formula is as follows:
[0025]
[0026]
[0027]
[0028] The simplified dynamic characteristic coefficients are then substituted into the single-frequency eddy current formula for calculation. and The rotor motion states at two different times are shown below: when hour, ; when hour, ; Where a and b are the major and minor axes of the ellipse, respectively; , , and These are the optimized four sealing dynamic characteristic coefficients.
[0029] Furthermore, the analysis and quantification module is also used for: The influence of univariate parameters on the dynamic characteristics of honeycomb sealing was analyzed using the correction coefficient method, and the results of the univariate influence analysis were obtained. Regression analysis was used to analyze the influence of multivariate parameters on the dynamic characteristics of honeycomb sealing, and the results of the multivariate influence analysis were obtained. By integrating the results of the univariate influence analysis and the multivariate influence analysis, a unified quantitative relationship between parameters and sealing dynamic characteristics is formed.
[0030] Furthermore, the device further includes a storage and execution module, the storage and execution module being used for: Store the correspondence between the input parameters of all samples and the sealing dynamic characteristic coefficients to form a structured database; Enter the fitted polynomial coefficients and correction coefficients, and establish the connection between the quantization relationship and the database; A visual input interface is adopted so that engineers can input sealing structure parameters and operating condition parameters through the visual input interface; The built-in calling logic invokes the quantization relationship, and the cell sealing dynamic characteristic coefficient is calculated and output based on the quantization relationship.
[0031] According to another aspect of the present invention, a storage medium is provided, wherein at least one executable instruction is stored therein, the executable instruction causing a processor to perform an operation corresponding to the above-described analysis method for the dynamic characteristics of cellular sealing.
[0032] According to another aspect of the present invention, an electronic device is provided, including a processor, a memory, a communication interface, and a communication bus, wherein the processor, the memory, and the communication interface communicate with each other through the communication bus; The memory is used to store at least one executable instruction that causes the processor to perform an operation corresponding to the analysis method of the dynamic characteristics of the cellular seal described above.
[0033] By employing the above-described technical solutions, the technical solutions provided by the embodiments of the present invention have at least the following advantages: This invention provides a method and apparatus for analyzing the dynamic characteristics of honeycomb seals, a storage medium, and an electronic device. Compared with existing technologies, this invention obtains a parameterized fluid domain model of the honeycomb seal by parametrically modeling the fluid domain in a honeycomb seal structure with gas swirl. Then, it performs efficiency optimization processing on the parameterized fluid domain model of the honeycomb seal to obtain an optimized CFD analysis model. The efficiency optimization processing includes simplification of the single-frequency eddy model and simplification of dynamic characteristic coefficients, realizing full-process optimization from modeling to analysis method. By conducting parameterized modeling of the honeycomb seal and optimizing and verifying the CFD analysis method, the computation time of the CFD method is significantly shortened without increasing costs and while ensuring accuracy under high-pressure conditions. This invention also generates a large number of calculation samples based on parameter combinations and uses an optimized CFD analysis model to calculate the sealing dynamic characteristics of each sample. By fitting and analyzing the sensitivity of the sealing dynamic characteristics to different parameters and establishing a quantitative relationship between parameters and sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping in the honeycomb sealing structure is obtained. This meets the needs of rapid engineering application of honeycomb sealing dynamic characteristics and provides a faster and more accurate basis for improving and optimizing the overall dynamic characteristics of the unit.
[0034] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description
[0035] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings: Figure 1 A flowchart illustrating an analysis method for the dynamic characteristics of a honeycomb seal provided in an embodiment of the present invention is shown. Figure 2 A flowchart illustrating another method for analyzing the dynamic characteristics of a honeycomb seal provided in an embodiment of the present invention is shown. Figure 3 A schematic diagram of the structure of an analysis device for the dynamic characteristics of a honeycomb seal provided in an embodiment of the present invention is shown; Figure 4 A schematic diagram of the structure of an electronic device provided by an embodiment of the present invention is shown. Detailed Implementation
[0036] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
[0037] Since existing analytical methods for obtaining the dynamic characteristics of honeycomb seals using Computational Fluid Dynamics (CFD) models lack experimental verification under high-pressure conditions, before developing the technical means of this invention, we first obtained experimental data of honeycomb seals at the highest pressure level from publicly available materials as a benchmark. We then used CFD methods to simulate and calculate the dynamic characteristics of the honeycomb seals under this high-pressure scenario. We compared the CFD prediction results with the experimental data, verifying parameters such as effective stiffness and effective damping, and verified the consistency using frequency domain response curves (excitation frequency-cpm dimension). The results show that the CFD analysis method under high pressure can accurately provide the dynamic characteristics of the honeycomb seals. This provides a reliable methodological foundation for subsequent full-process optimization and ensures the credibility of the calculation results in engineering applications.
[0038] This invention provides a method for analyzing the dynamic characteristics of honeycomb seals, such as... Figure 1 As shown, the method includes: 101. Parametric modeling of the honeycomb sealing fluid domain in the honeycomb sealing structure with gas spool is performed to obtain the parametric fluid domain model of the honeycomb sealing. In this embodiment of the invention, the current execution end performs parametric modeling of the honeycomb sealing fluid domain in the honeycomb sealing structure with gas rotor. The specific steps are as follows: (1) Obtain the complete swirl-blocking structure model established in advance, and analyze the influence of the complete swirl-blocking structure model on the pre-swirl of the airflow inlet; (2) Based on the influence law, the pre-spinning parameter is used to replace the solid structure, which indirectly reflects the inlet gas anti-spinning effect. In this embodiment, the tedious process of solid modeling of the honeycomb sealing structure with gas anti-spinning can be avoided by using the above steps (1) and (2), and the workload of geometric model construction can be reduced by indirectly reflecting the inlet gas anti-spinning effect. (3) A parameterized interface is built based on the modeling platform, and the key dimensions of the honeycomb seal are set as editable parameters; wherein, the key dimensions of the honeycomb seal include axial width, diameter, gap, etc., and the embodiments of the present invention do not make specific limitations.
[0039] (4) The pre-swirl parameter that indirectly reflects the inlet gas swirling effect is adapted in conjunction with the editable key dimensions of the honeycomb seal so that when the user edits the key dimensions of the honeycomb seal, the model is automatically updated and the corresponding inlet gas swirling effect is loaded synchronously, and finally the parameterized fluid domain model of the honeycomb seal that integrates the inlet gas swirling effect is obtained.
[0040] In this embodiment of the invention, the automatic model update and synchronous loading of the corresponding inlet gas swirling effect are mainly used as the inlet boundary conditions in the CFD analysis process. When the key dimensions of the honeycomb seal edited by the user change, the inlet boundary conditions of the CFD analysis will also change accordingly. In this embodiment, by replacing manual adjustment with parameterized driving through the above steps (3) and (4), the needs of multi-parameter iterative calculation in engineering can be adapted, which can greatly shorten the modeling time. At the same time, it combines the law of inlet gas swirling effect, which simplifies the CFD analysis without losing the core influencing factors.
[0041] 102. The efficiency of CFD analysis is optimized by combining the parameterized fluid domain model of the honeycomb seal to obtain the optimized CFD analysis model; the efficiency optimization includes simplification of the single-frequency eddy model and simplification of the dynamic characteristic coefficients; In this embodiment of the invention, in order to improve the overall efficiency of CFD analysis, the current execution end also performs efficiency optimization processing on the simplified honeycomb seal parameterized fluid domain model obtained in step 101. The specific optimization method is as follows: the multi-frequency eddy model in the CFD analysis process is simplified into a single-frequency eddy model to reduce unnecessary frequency domain calculations; the single-frequency eddy model adopts the rotor single-frequency elliptical eddy equation.
[0042] In this embodiment of the invention, the multi-frequency eddy current model in the original CFD analysis process is represented by the following formula:
[0043] in, x for x The eddy trajectory corresponding to the coordinates; y for y The coordinates correspond to the eddy trajectory; a and b are the major and minor axes of the ellipse, respectively; Indicates the vortex angular velocity. The eddy frequency is represented by N; the number of frequencies is represented by i; the frequency type is represented by t; and time is represented by t. In this embodiment, the analysis efficiency is improved by simplifying the model. The specific simplification process is as follows: (1) When the rotor is subjected to small displacement and velocity disturbances, the relationship between the airflow force on the rotor and the disturbance displacement and velocity is linearized based on the single-frequency eddy model as follows:
[0044] in, F x This is the excitation force in the x-direction of the airflow; F y This is the excitation force in the y-direction of the airflow; K Principal stiffness; k For cross stiffness; C Main damping; c For cross damping; X Displacement of rotor shaft center x direction; Y Displacement of rotor shaft center y direction; Rotor shaft disturbance speed x direction; Rotor shaft disturbance speed y direction; (2) By performing a fast Fourier transform on the linearized relationship, the relationship between the change in sealing gas flow force in the frequency domain and the sealing dynamic characteristic coefficient and the small disturbance displacement is obtained as follows:
[0045] in, ; It is the vortex angular velocity; K xx , K xy , K yx and K yy These are the four stiffness coefficients in the single-frequency eddy coulomb model; C xx , C xy , C yx and C yy These are the four damping coefficients in the single-frequency eddy current model; D x , D y ) and (Δ F x Δ F y These are the time-domain signals of the rotor eddy displacement and the frequency-domain signals corresponding to the time-domain values of the aerodynamic force changes, respectively. (3) Utilizing the symmetry of the dynamic characteristic coefficients when the rotor whirls around the axis, the eight dynamic characteristic coefficients in the single-frequency whirl model are analyzed. K xx , K xy , K yx , K yy , C xx , C xy , C yx , C yy After simplification, the specific formula is as follows:
[0046]
[0047]
[0048]
[0049] The simplified dynamic characteristic coefficients are then substituted into the single-frequency eddy current formula for calculation. and The rotor motion states at two different times are shown below: when hour, ; when hour, ; Where a and b are the major and minor axes of the ellipse, respectively; , , and These are the optimized four sealing dynamic characteristic coefficients.
[0050] In this embodiment of the invention, steps 101 and 102 achieve a dual optimization process of model simplification and coefficient simplification, which can significantly improve the efficiency of CFD analysis, solve the problem of long calculation time in traditional methods, reduce hardware resource consumption, and support rapid calculation of large-scale samples. For example, in engineering applications, users input key dimensions of the honeycomb seal and key simulation parameters, including medium, temperature, rotational speed, and sealing gap, through a parameterized interface. This embodiment of the invention does not impose specific limitations. Then, the model is automatically updated to a parameterized fluid domain model of the honeycomb seal adapted to the above parameters; then, the optimized CFD analysis is started, and rapid calculation is performed based on a single-frequency eddy current model, utilizing coefficient symmetry to obtain four core sealing dynamic characteristic coefficients.
[0051] It should be noted that, before the efficiency optimization process of CFD analysis performed by the current execution end in conjunction with the parameterized fluid domain model of the cellular seal in this embodiment, the method further includes: Verify the adaptability of the parameterized fluid domain model of the honeycomb seal to single-frequency eddy; for example, verify whether the core features retained in the parameterized fluid domain model of the honeycomb seal, such as the arrangement of honeycomb holes, sealing gap, and inlet and outlet flow channels, can accurately simulate the airflow state under single-frequency eddy, including the airflow response to rotor eddy and pressure distribution changes. The embodiments of the present invention do not make specific limitations.
[0052] After successful verification, the pre-swirl parameters are incorporated as inlet boundary conditions into the flow field calculation of single-frequency vortex, ensuring that the inlet gas swirl effect is accurately captured under single-frequency operating conditions. In this embodiment, the pre-swirl parameters are incorporated as inlet boundary conditions into the flow field calculation of single-frequency vortex, meaning that the gas flow force changes calculated by the single-frequency vortex model already include the regulating effect of the swirl structure on the inlet airflow.
[0053] 103. Generate a large number of computational samples based on parameter combinations, and use the optimized CFD analysis model to calculate the sealing dynamic characteristics of each sample group; In this embodiment of the invention, the current execution terminal also acquires a pre-built database of dynamic characteristics of honeycomb seals, combines the associated parameter information in the database, such as combining honeycomb seal structure parameters and operating condition parameters, including core parameters such as medium, pre-swirl ratio, temperature, axial width, rotational speed, diameter, inlet / outlet pressure, and sealing gap, and generates a large number of calculation samples based on the above parameters. In addition, the current execution terminal also uses the CFD analysis model optimized in steps 101 and 102 to calculate the sealing dynamic characteristic results of each group of samples, including calculating stiffness, damping coefficient, etc., which are not specifically limited in this embodiment of the invention.
[0054] It should be noted that the calculated sealing dynamic characteristics of each group of samples, together with the calculated samples, form the underlying database samples for subsequent fitting analysis.
[0055] 104. By fitting and analyzing the sensitivity of the sealing dynamic characteristics results to different parameters, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb sealing structure is obtained.
[0056] In this embodiment of the invention, the current execution end analyzes the sensitivity of the sealing dynamic characteristic results to different parameters through fitting analysis, and establishes a quantitative relationship between the parameters and the sealing dynamic characteristics. The specific process includes: (1) The influence of univariate parameters on the dynamic characteristics of honeycomb sealing was analyzed by the correction coefficient method, and the results of the univariate influence analysis were obtained; In this embodiment of the invention, the current execution end performs curve fitting on parameters such as different media, temperatures, rotational speeds, axial widths, sealing diameters, pre-swirl ratios, and sealing gaps to obtain the influence law of these parameters on the stiffness damping coefficients. First, the correction coefficient method is selected to analyze the influence of single-variable parameters on the dynamic characteristics of the honeycomb seal, obtaining the single-variable influence analysis results corresponding to each single variable. The correction coefficient method, also known as the calibration coefficient method, is a method of correcting coefficients to improve the accuracy, balance, and completeness of measurement data. It is also a coefficient added to the calculation formula to reflect the true performance as much as possible when deviations occur in data calculation and formula expression due to ideal versus reality, or reality versus survey results. In this embodiment of the invention, the correction coefficient method is used... β express.
[0057]
[0058] in, x Indicates a single-variable parameter. β 1 、β 2 、β 3 represents the correction factor.
[0059] (2) The influence of multivariate parameters on the dynamic characteristics of honeycomb sealing was analyzed by regression analysis, and the results of the multivariate influence analysis were obtained; In this embodiment of the invention, the current execution end also uses regression analysis to analyze the influence of multivariate parameters on the dynamic characteristics of the honeycomb seal, and obtains the results of the multivariate influence analysis; for example, variables such as temperature, rotational speed, axial width, sealing diameter, and pre-swirl ratio are determined as input variables, and the following nonlinear polynomial regression equation is used for regression analysis:
[0060] in, x 1 、x 2 、...、x n These represent the input variables such as temperature, rotational speed, axial width, sealing diameter, and pre-rotation ratio. β 1 、β 2 、...、β i , β i+1 ... β n 、β n+1 They represent the coefficients of the polynomial regression, y These are the predicted dynamic characteristic coefficients, i.e., the results of the multivariate influence analysis.
[0061] (3) Integrate the results of the univariate influence analysis and the multivariate influence analysis to form a unified quantitative relationship between parameters and sealing dynamic characteristics.
[0062] In this embodiment of the invention, the current execution end integrates the results of univariate influence analysis and multivariate influence analysis using the following formula to form a unified quantitative relationship between parameters and sealing dynamic characteristics:
[0063] in, x 1 to x n These are the individual variable parameters used in the correction coefficient method analysis, representing one of the variables such as temperature, rotational speed, axial width, sealing diameter, and pre-rotation ratio. x n+1 to x m The input parameters used in the regression analysis method represent temperature, rotation speed, axial width, sealing diameter, pre-rotation ratio, etc., and the embodiments of the present invention do not impose specific limitations.
[0064] It should be noted that in this embodiment, due to the large differences in the values of different parameters, the large parameters overwhelm the small parameters when the computer constructs the polynomial regression equation. Therefore, parameter normalization is required to reduce the differences between different parameters and improve the adaptability of the regression equation.
[0065] Furthermore, as a refinement and extension of the specific implementation of the above embodiments, and in order to improve the ability to rapidly apply these methods in engineering, another method for analyzing the dynamic characteristics of honeycomb seals is provided, such as... Figure 2 As shown, the method further includes: 201. Store the correspondence between the input parameters of all samples and the sealing dynamic characteristic coefficients to form a structured database; 202. Enter the fitted polynomial coefficients and correction coefficients, and establish the connection between the quantization relationship and the database; 203. A visual input interface is adopted so that engineers can input sealing structure parameters and operating condition parameters through the visual input interface; 204. The built-in calling logic is used to invoke the quantization relationship, and the dynamic characteristic coefficient of the honeycomb seal is calculated and output based on the quantization relationship.
[0066] In this embodiment of the invention, the current execution end realizes the rapid engineering application of the method of this application through the above steps 201 to 204. Without the need for professional CFD analysis, modeling and calculation knowledge, the dynamic characteristic coefficients under the target parameters can be quickly obtained through software. The calculation method, fitting law and sample data are integrated into a standardized tool, which is convenient for subsequent iterative upgrades and promotion.
[0067] This invention provides a method for analyzing the dynamic characteristics of honeycomb seals. Compared with existing technologies, this invention performs parametric modeling of the honeycomb seal fluid domain in a honeycomb seal structure with gas swirl, obtaining a parametric fluid domain model of the honeycomb seal. Then, it performs efficiency optimization processing on the parametric fluid domain model of the honeycomb seal to obtain an optimized CFD analysis model. The efficiency optimization processing includes simplification of the single-frequency eddy model and simplification of dynamic characteristic coefficients, realizing full-process optimization from modeling to analysis. By performing parametric modeling of the honeycomb seal and optimizing and verifying the CFD analysis method, the computation time of the CFD method is significantly shortened without increasing costs and ensuring accuracy under high-pressure conditions. This invention also generates a large number of calculation samples based on parameter combinations and uses the optimized CFD analysis model to calculate the sealing dynamic characteristic results of each sample. By fitting and analyzing the sensitivity of the sealing dynamic characteristic results to different parameters and establishing a quantitative relationship between parameters and sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb seal structure is obtained, meeting the needs of rapid engineering application of honeycomb seal dynamic characteristics and providing a faster and more accurate basis for improving and optimizing the overall dynamic characteristics of the unit.
[0068] As a response to the above Figure 1 The implementation of the method shown in this invention provides an analytical device for the dynamic characteristics of honeycomb seals, such as... Figure 3 As shown, the device includes: The parametric modeling module 31 is used to perform parametric modeling on the honeycomb sealing fluid domain in the honeycomb sealing structure with gas spool, and obtain the honeycomb sealing parametric fluid domain model. CFD optimization module 32 is used to perform efficiency optimization processing of CFD analysis in combination with the parameterized fluid domain model of the honeycomb seal, so as to obtain an optimized CFD analysis model; the efficiency optimization processing includes simplification processing of single-frequency eddy model and simplification processing of dynamic characteristic coefficients; The data preparation module 33 is used to generate a large number of computational samples based on parameter combinations, and to calculate the sealing dynamic characteristics of each sample using an optimized CFD analysis model. The analysis and quantification module 34 is used to analyze the sensitivity of the sealing dynamic characteristic results to different parameters through fitting analysis, and to establish a quantitative relationship between the parameters and the sealing dynamic characteristics, so as to obtain the influence law of different types of parameters on the stiffness and damping in the honeycomb sealing structure.
[0069] Furthermore, the parametric modeling module 31 is also used for: Obtain the pre-established complete swirl block structure model and analyze the influence of the complete swirl block structure model on the pre-swirl of the airflow inlet. Based on the aforementioned influence law, the pre-spinning parameter is used to replace the solid structure to indirectly reflect the inlet gas spinning effect; A parametric interface was built based on the modeling platform, and the key dimensions of the honeycomb seal were set as editable parameters. The pre-swirl parameter, which indirectly reflects the inlet gas swirl effect, is adapted in conjunction with the editable key dimensions of the honeycomb seal. This allows the model to be automatically updated and the corresponding inlet gas swirl effect to be loaded synchronously when the user edits the key dimensions of the honeycomb seal, ultimately resulting in a parameterized fluid domain model of the honeycomb seal that incorporates the inlet gas swirl effect.
[0070] Furthermore, the device also includes a verification module, which is used for: Verify the adaptability of the proposed honeycomb seal parameterized fluid domain model to single-frequency eddy currents; After successful verification, the pre-swirl parameters were used as inlet boundary conditions and incorporated into the flow field calculation of single-frequency eddy currents to ensure that the inlet gas swirl effect was accurately captured under single-frequency conditions.
[0071] Furthermore, the CFD optimization module 32 is also used for: The multi-frequency eddy model in the CFD analysis process is simplified into a single-frequency eddy model to reduce unnecessary frequency domain calculations; the single-frequency eddy model adopts the rotor single-frequency elliptical eddy equation.
[0072] Furthermore, the CFD optimization module 32 is also used for: When the rotor is subjected to minute displacement and velocity disturbances, the relationship between the airflow force on the rotor and the disturbance displacement and velocity is linearized based on the single-frequency eddy current model as follows:
[0073] in, F x This is the excitation force in the x-direction of the airflow; F y This is the excitation force in the y-direction of the airflow; K Principal stiffness; k For cross stiffness; C Main damping; c For cross damping; X Displacement of rotor shaft center x direction; Y Displacement of rotor shaft center y direction; Rotor shaft disturbance speed x direction; Rotor shaft disturbance speed y direction; Performing a Fast Fourier Transform on the linearized relationship yields the following equations: The relationship between the change in sealing gas flow force in the frequency domain and the sealing dynamic characteristic coefficient and the displacement of small disturbances is:
[0074] in, ; It is the vortex angular velocity; K xx , K xy , K yx and K yy These are the four stiffness coefficients in the single-frequency eddy coulomb model; C xx , C xy , C yx and C yy These are the four damping coefficients in the single-frequency eddy current model; D x , D y ) and (Δ F x Δ F y These are the time-domain signals of the rotor eddy displacement and the frequency-domain signals corresponding to the time-domain values of the aerodynamic force changes, respectively. Utilizing the symmetry of the dynamic characteristic coefficients during rotor whirling around the axis, the eight dynamic characteristic coefficients in the single-frequency whirling model are analyzed. K xx , K xy , K yx , K yy , C xx , C xy , C yx , C yy After simplification, the specific formula is as follows:
[0075]
[0076]
[0077]
[0078] The simplified dynamic characteristic coefficients are then substituted into the single-frequency eddy current formula for calculation. and The rotor motion states at two different times are shown below: when hour, ; when hour, ; Where a and b are the major and minor axes of the ellipse, respectively; , , and These are the optimized four sealing dynamic characteristic coefficients.
[0079] Furthermore, the analysis and quantification module 34 is also used for: The influence of univariate parameters on the dynamic characteristics of honeycomb sealing was analyzed using the correction coefficient method, and the results of the univariate influence analysis were obtained. Regression analysis was used to analyze the influence of multivariate parameters on the dynamic characteristics of honeycomb sealing, and the results of the multivariate influence analysis were obtained. By integrating the results of the univariate influence analysis and the multivariate influence analysis, a unified quantitative relationship between parameters and sealing dynamic characteristics is formed.
[0080] Furthermore, the device further includes a storage and execution module, the storage and execution module being used for: Store the correspondence between the input parameters of all samples and the sealing dynamic characteristic coefficients to form a structured database; Enter the fitted polynomial coefficients and correction coefficients, and establish the connection between the quantization relationship and the database; A visual input interface is adopted so that engineers can input sealing structure parameters and operating condition parameters through the visual input interface; The built-in calling logic invokes the quantization relationship, and the cell sealing dynamic characteristic coefficient is calculated and output based on the quantization relationship.
[0081] This invention provides an analysis device for the dynamic characteristics of honeycomb seals. Compared with existing technologies, this invention performs parametric modeling of the honeycomb seal fluid domain in a honeycomb seal structure with gas vortex, obtaining a parametric fluid domain model for the honeycomb seal. It then performs efficiency optimization processing on the parametric fluid domain model for CFD analysis, resulting in an optimized CFD analysis model. This efficiency optimization includes simplification of the single-frequency eddy model and simplification of dynamic characteristic coefficients, achieving full-process optimization from modeling to analysis methods. By performing parametric modeling of the honeycomb seal and optimizing and verifying the CFD analysis method, the computation time of the CFD method is significantly shortened without increasing costs and while ensuring accuracy under high-pressure conditions. Furthermore, this invention generates a large number of calculation samples based on parameter combinations and uses the optimized CFD analysis model to calculate the sealing dynamic characteristics of each sample. By fitting and analyzing the sensitivity of the sealing dynamic characteristics to different parameters and establishing a quantitative relationship between parameters and sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb seal structure is obtained. This meets the needs of rapid engineering application of honeycomb seal dynamic characteristics and provides a faster and more accurate basis for improving and optimizing the overall dynamic characteristics of the unit.
[0082] According to one embodiment of the present invention, a storage medium is provided, the storage medium storing at least one executable instruction, the computer-executable instruction being capable of executing the method for analyzing the dynamic characteristics of the cellular seal in any of the above method embodiments.
[0083] Figure 4 The diagram illustrates the structure of an electronic device according to an embodiment of the present invention. The specific embodiments of the present invention do not limit the specific implementation of the electronic device.
[0084] like Figure 4 As shown, the electronic device may include: a processor 402, a communications interface 404, a memory 406, and a communications bus 408.
[0085] The processor 402, communication interface 404, and memory 406 communicate with each other via communication bus 408.
[0086] Communication interface 404 is used to communicate with other network elements such as clients or other servers.
[0087] The processor 402 is used to execute program 410, which can specifically perform the relevant steps of the analysis method of the dynamic characteristics of the honeycomb seal described above.
[0088] Specifically, program 410 may include program code that includes computer operation instructions.
[0089] Processor 402 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement embodiments of the present invention. The electronic device includes one or more processors, which may be processors of the same type, such as one or more CPUs; or processors of different types, such as one or more CPUs and one or more ASICs.
[0090] Memory 406 is used to store program 410. Memory 406 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.
[0091] Specifically, program 410 can be used to cause processor 402 to perform the following operations: The honeycomb sealing fluid domain in the honeycomb sealing structure with gas spool is parametrically modeled to obtain the honeycomb sealing parametric fluid domain model. The efficiency of CFD analysis is optimized by combining the aforementioned parameterized fluid domain model of the honeycomb seal to obtain the optimized CFD analysis model; the efficiency optimization includes simplification of the single-frequency eddy model and simplification of the dynamic characteristic coefficients; A large number of computational samples were generated based on parameter combinations, and the sealing dynamic characteristics of each sample were calculated using an optimized CFD analysis model. By fitting and analyzing the sensitivity of the sealing dynamic characteristics results to different parameters, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb sealing structure is obtained.
[0092] It is obvious to those skilled in the art that the modules or steps of the present invention described above can be implemented using general-purpose computing devices. They can be centralized on a single computing device or distributed across a network of multiple computing devices. Optionally, they can be implemented using computer-executable program code, thereby storing them in a storage device for execution by a computing device. In some cases, the steps shown or described can be performed in a different order than those presented herein, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any particular combination of hardware and software.
[0093] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.
Claims
1. A method for analyzing the dynamic characteristics of honeycomb seals, characterized in that, include: The honeycomb sealing fluid domain in the honeycomb sealing structure with gas spool is parametrically modeled to obtain the honeycomb sealing parametric fluid domain model. The efficiency of CFD analysis is optimized by combining the aforementioned parameterized fluid domain model of the honeycomb seal to obtain the optimized CFD analysis model; the efficiency optimization includes simplification of the single-frequency eddy model and simplification of the dynamic characteristic coefficients; A large number of computational samples were generated based on parameter combinations, and the sealing dynamic characteristics of each sample were calculated using an optimized CFD analysis model. By fitting and analyzing the sensitivity of the sealing dynamic characteristics results to different parameters, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, the influence of different types of parameters on the stiffness and damping of the honeycomb sealing structure is obtained.
2. The method according to claim 1, characterized in that, The parameterized modeling of the honeycomb sealing fluid domain in the gas-swirled honeycomb sealing structure yields a parameterized fluid domain model for the honeycomb sealing, including: Obtain the pre-established complete swirl block structure model and analyze the influence of the complete swirl block structure model on the pre-swirl of the airflow inlet. Based on the aforementioned influence law, the pre-spinning parameter is used to replace the solid structure to indirectly reflect the inlet gas spinning effect; A parametric interface was built based on the modeling platform, and the key dimensions of the honeycomb seal were set as editable parameters. The pre-swirl parameter, which indirectly reflects the inlet gas swirl effect, is adapted in conjunction with the editable key dimensions of the honeycomb seal. This allows the model to be automatically updated and the corresponding inlet gas swirl effect to be loaded synchronously when the user edits the key dimensions of the honeycomb seal, ultimately resulting in a parameterized fluid domain model of the honeycomb seal that incorporates the inlet gas swirl effect.
3. The method according to claim 1, characterized in that, Before the efficiency optimization process of CFD analysis based on the parameterized fluid domain model of the cellular seal, the method further includes: Verify the adaptability of the described honeycomb seal parameterized fluid domain model to single-frequency eddy currents; After successful verification, the pre-swirl parameters were incorporated as inlet boundary conditions into the flow field calculation of single-frequency eddy, so that the inlet gas swirling effect could be accurately captured under single-frequency conditions.
4. The method according to claim 1, characterized in that, The efficiency optimization process for CFD analysis based on the parameterized fluid domain model of the honeycomb seal results in an optimized CFD analysis model, including: The multi-frequency eddy model in the CFD analysis process is simplified to a single-frequency eddy model to reduce unnecessary frequency domain calculations; the single-frequency eddy model adopts the rotor single-frequency elliptical eddy equation.
5. The method according to claim 4, characterized in that, The efficiency optimization process for CFD analysis based on the parameterized fluid domain model of the honeycomb seal, resulting in an optimized CFD analysis model, further includes: When the rotor is subjected to minute displacement and velocity disturbances, the relationship between the airflow force on the rotor and the disturbance displacement and velocity is linearized based on the single-frequency eddy current model as follows: in, F x This is the excitation force in the x-direction of the airflow; F y This is the excitation force in the y-direction of the airflow; K Principal stiffness; k For cross stiffness; C Main damping; c For cross damping; X Displacement of rotor shaft center x direction; Y Displacement of rotor shaft center y direction; Rotor shaft disturbance speed x direction; Rotor shaft disturbance speed y direction; Performing a Fast Fourier Transform on the linearized relationship yields the following equations: The relationship between the change in sealing gas flow force in the frequency domain and the sealing dynamic characteristic coefficient and the displacement of small disturbances is: in, ; It is the vortex angular velocity; K xx , K xy , K yx and K yy These are the four stiffness coefficients in the single-frequency eddy coulomb model; C xx , C xy , C yx and C yy These are the four damping coefficients in the single-frequency eddy current model; D x , D y ) and (Δ F x Δ F y These are the time-domain signals of the rotor eddy displacement and the frequency-domain signals corresponding to the time-domain values of the aerodynamic force changes, respectively. Utilizing the symmetry of the dynamic characteristic coefficients during rotor whirling around the axis, the eight dynamic characteristic coefficients in the single-frequency whirling model are analyzed. K xx , K xy , K yx , K yy , C xx , C xy , C yx , C yy After simplification, the specific formula is as follows: The simplified dynamic characteristic coefficients are then substituted into the single-frequency eddy current formula for calculation. and The rotor motion states at two different times are shown below: when hour, ; when hour, ; Where a and b are the major and minor axes of the ellipse, respectively; , , and These are the optimized four sealing dynamic characteristic coefficients.
6. The method according to claim 1, characterized in that, The process of fitting and analyzing the sensitivity of the sealing dynamic characteristics results to different parameters, and establishing a quantitative relationship between the parameters and the sealing dynamic characteristics, includes: The influence of univariate parameters on the dynamic characteristics of honeycomb sealing was analyzed using the correction coefficient method, and the results of the univariate influence analysis were obtained. Regression analysis was used to analyze the influence of multivariate parameters on the dynamic characteristics of honeycomb sealing, and the results of the multivariate influence analysis were obtained. By integrating the results of the univariate influence analysis and the multivariate influence analysis, a unified quantitative relationship between parameters and sealing dynamic characteristics is formed.
7. The method according to any one of claims 1 to 6, characterized in that, The method further includes: Store the correspondence between the input parameters of all samples and the sealing dynamic characteristic coefficients to form a structured database; Enter the fitted polynomial coefficients and correction coefficients, and establish the connection between the quantization relationship and the database; A visual input interface is adopted so that engineers can input sealing structure parameters and operating condition parameters through the visual input interface; The built-in calling logic invokes the quantization relationship, and the cell sealing dynamic characteristic coefficient is calculated and output based on the quantization relationship.
8. An analytical device for the dynamic characteristics of a honeycomb seal, characterized in that, include: The parametric modeling module is used to perform parametric modeling of the honeycomb sealing fluid domain in a honeycomb sealing structure with gas spool, and obtain a parametric fluid domain model of the honeycomb sealing. The CFD optimization module is used to perform efficiency optimization processing on the CFD analysis in conjunction with the parameterized fluid domain model of the honeycomb seal, so as to obtain an optimized CFD analysis model; the efficiency optimization processing includes simplification processing of the single-frequency eddy model and simplification processing of the dynamic characteristic coefficients; The data preparation module is used to generate a large number of computational samples based on parameter combinations, and to calculate the sealing dynamic characteristics of each sample using an optimized CFD analysis model. The analysis and quantification module is used to analyze the sensitivity of the sealing dynamic characteristics results to different parameters through fitting analysis, and to establish a quantitative relationship between the parameters and the sealing dynamic characteristics, so as to obtain the influence law of different types of parameters on the stiffness and damping in the honeycomb sealing structure.
9. A storage medium, characterized in that, The storage medium stores at least one executable instruction that performs the operation corresponding to the analysis method of the dynamic characteristics of the cellular seal as described in any one of claims 1-7.
10. An electronic device, characterized in that, It includes a processor, a memory, a communication interface, and a communication bus, wherein the processor, the memory, and the communication interface communicate with each other through the communication bus; The memory is used to store at least one executable instruction that causes the processor to perform an operation corresponding to the analysis method of the dynamic characteristics of the cellular seal as described in any one of claims 1-7.