A method and system for analyzing design parameter sensitivity of an anti-rotation plate

Abstract

The application provides a kind of anti-rotation plate design parameter sensitivity analysis method and system, method includes: selecting the geometric parameter of k anti-rotation plate as test factor, and based on central composite design test method generates the code table containing multiple anti-rotation plate geometric parameter combination;Based on the code table, generate multiple anti-rotation plate sealing structure, and carry out numerical simulation of steady-state rotor static eccentricity to each anti-rotation plate sealing structure under the same working condition, to obtain the index value of each anti-rotation plate sealing structure stability evaluation index;Based on stability evaluation index analysis anti-rotation plate parameter sensitivity to anti-rotation plate design, and determine the preferred anti-rotation plate design parameter combination according to the analysis result.The application is mainly based on test design DOE method to analyze the sensitivity of key geometric parameters and then obtain the preferred anti-rotation plate design parameter combination, which helps turbine designers to quickly and efficiently meet the anti-rotation plate design requirements in the engineering design process.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of performance analysis, in particular to a method and system for analyzing sensitivity of design parameters of anti-rotation plates. BACKGROUND

[0002] The annular seal structure is installed between the rotor and the stator in a turbine, and controls the leakage under high pressure difference. In steam turbines, gas turbines, pumps and compressors, the rotating seal is installed between the rotating part and the stationary part, and controls the leakage flow from the high pressure area to the low pressure area through the dynamic gap, which has a significant impact on the operating efficiency of the turbomachinery.

[0003] There are two aspects to evaluate the sealing performance of the annular seal structure: one is the sealing leakage characteristic, and the other is the sealing rotor dynamic characteristic. Compared with the damper seal structure such as the hole type seal and the bag type seal, the non-damper seal is more likely to cause insufficient rotor stability due to the unsteady airflow excitation force generated by the rotor whirl. Compared with the low inlet pre-rotation working condition seal, the cross stiffness increases significantly under the medium and high inlet pre-rotation working condition, and the airflow excitation force weakens the ability to suppress the forward whirl of the rotor. Therefore, the anti-rotation plate structure is usually arranged at the inlet of the non-damper seal to block or reverse guide the circumferential airflow to improve the stability of the sealing rotor.

[0004] Currently, the anti-rotation performance analysis of the seal anti-rotation plate is carried out by a large number of experiments and numerical calculations. However, there are three problems: the cost of a large number of experiments and numerical calculations is too high; the performance comparison of individual anti-rotation plates under specific conditions may not be suitable for new seal structures and working conditions; and the key geometric parameters affecting the performance of the anti-rotation plate cannot be determined. SUMMARY

[0005] Therefore, the technical problem to be solved by the present application is to overcome the defect of sensitivity analysis of key geometric parameters in the design process of the anti-rotation plate in the prior art, so as to provide a method and system for analyzing sensitivity of design parameters of anti-rotation plates.

[0006] The technical solution of the present application to solve the above technical problems is as follows:

[0007] In a first aspect, the present application provides a method for analyzing sensitivity of design parameters of anti-rotation plates, comprising the following steps:

[0008] Selecting k geometric parameters of the anti-rotation plate as test factors, and generating a code table containing various combinations of the geometric parameters of the anti-rotation plate based on a central composite design test method;

[0009] Generating a plurality of anti-rotation plate seal structures based on the code table, and performing numerical simulation of the steady-state rotor static eccentricity of each anti-rotation plate seal structure under the same working condition to obtain the index value of the stability evaluation index of each anti-rotation plate seal structure;

[0010] The stability evaluation index is analyzed based on the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design, and the optimal anti-rotation plate design parameter combination is determined according to the analysis result.

[0011] The anti-rotation plate design parameter sensitivity analysis method provided by the embodiment of the application can process the geometric parameters of the anti-rotation plate to obtain an encoding table, generate a plurality of anti-rotation plate sealing structures according to the encoding table, obtain the stability evaluation index of the sealing structure through data simulation, and analyze the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the index, so that the key geometric parameters of the anti-rotation plate stability enhancement and vibration suppression can be identified and the optimal anti-rotation plate design parameter combination can be obtained through a limited number of tests, which helps the turbine designers to quickly and efficiently meet the anti-rotation plate design requirements in the engineering design process.

[0012] Optionally, the geometric parameters include the anti-rotation plate length, the anti-rotation plate number, the anti-rotation plate stagger angle and the anti-rotation plate inclination angle.

[0013] The application selects a plurality of geometric parameters of the anti-rotation plate as the key geometric parameters to be investigated, analyzes the sensitivity of each geometric parameter to the anti-rotation plate design through subsequent test, and finally determines the key geometric parameter sensitivity of the anti-rotation plate design, so that the analysis range is more comprehensive and the obtained result is more reliable.

[0014] Optionally, the process of selecting the k geometric parameters of the anti-rotation plate as the test factors and generating the encoding table containing a plurality of anti-rotation plate geometric parameter combinations based on the central composite design test method includes: uniformly selecting 5 values for 4 geometric parameters in the required analysis value range; generating a four-factor five-level encoding table based on the parameter test of the central composite design method, and the encoding table gives a set of anti-rotation plate geometric parameter combinations for each item.

[0015] The embodiment of the application identifies the key parameters to obtain the optimal design by using the DOE module in the Isight parameter optimization, selects 4 geometric parameters as test factors by using the central composite design method, and uniformly divides the 4 geometric parameters into 5 levels to generate a four-factor five-level encoding table containing a plurality of anti-rotation plate geometric parameter combinations. The method expands the design space and obtains high-order information by adding the extreme points and the center points on the basis of the two-level factorial design, can provide sample data for the response surface approximation model, has the advantages of simple design, less test times, good predictability and the like, and can also ensure the reliability of the test result.

[0016] Optionally, the encoding table includes: 2 k +2k+1 kinds of anti-rotation plate geometric parameter combinations, corresponding to 2 kThe +2k+1 anti-whirl plate sealing design schemes include 2k k-factor two-level full-factor test schemes, 2k shaft point test schemes and one center point test scheme.

[0017] The application generates corresponding anti-whirl plate sealing design schemes by various anti-whirl plate geometric parameter combinations contained in the coding table, considers the influence degree of different parameters and parameter combinations on anti-whirl plate design, and the obtained sensitivity analysis result is more comprehensive.

[0018] Optionally, the stability evaluation index includes an average circumferential velocity of the downstream airflow of the anti-whirl plate and a tangential airflow excitation force of the rotor.

[0019] The application selects the two indexes of the average circumferential velocity of the airflow and the tangential airflow excitation force of the rotor to evaluate the stability of the design scheme, because the two indexes can accurately evaluate the stability of the sealing rotor, and the design scheme obtained by using the two indexes as evaluation indexes can better meet the anti-whirl plate design requirements in actual engineering design.

[0020] Optionally, the process of obtaining the index values of the stability evaluation indexes of each anti-whirl plate sealing structure includes: using a numerical simulation method of solving RANS equations to calculate the corresponding index values of the average circumferential velocity of the downstream airflow of the anti-whirl plate and the tangential airflow excitation force of the rotor of each anti-whirl plate sealing structure.

[0021] The application uses the method of solving RANS equations to obtain the stability evaluation indexes under the static eccentricity of the rotor, and the method of solving RANS equations is a main method for numerically simulating complex viscous flow fields at present, and the accuracy of the RANS equations after adding a turbulence model can be greatly improved in calculating lift, drag and moment. In the application, the calculation of the average circumferential velocity of the downstream airflow of the anti-whirl plate and the tangential airflow excitation force of the rotor by the RANS equations belongs to an example of calculating the drag by the RANS equations.

[0022] Optionally, the process of analyzing the sensitivity of the anti-whirl plate parameters to the anti-whirl plate design based on the stability evaluation indexes includes: normalizing the value range of the test factor to [-1, 1] to obtain a normalized coding table; performing quadratic polynomial fitting on the normalized coding table and the stability evaluation indexes of each sealing structure; obtaining the contribution rate and influence effect of each test factor to the stability evaluation indexes based on the fitting coefficients of the polynomials; and analyzing the sensitivity of the anti-whirl plate parameters to the anti-whirl plate design according to the obtained contribution rate and influence effect.

[0023] The embodiment of the present application performs quadratic polynomial fitting on the results of various anti-rotation plate geometric parameter combinations and corresponding design schemes by using commercial software Isight, and obtains the contribution rate and influence effect of each test factor on the stability evaluation index according to the results after fitting, and further analyzes the sensitivity of the anti-rotation plate design parameters. The whole data analysis grasps the changes and the correlation of the design parameters and the stability evaluation index after the test from the perspective of various parameters, so that the final analysis result is more reliable. In addition, the contribution rate and influence effect of each parameter of the anti-rotation plate on the stability evaluation index of the sealing design scheme are displayed through the Pateto chart and the main effect chart, so that the analysis process is more intuitive, and the complex data analysis work becomes simpler.

[0024] In a second aspect, the embodiment of the present application provides an anti-rotation plate design parameter sensitivity analysis system, the system comprises:

[0025] The encoding table generation module is configured to select k geometric parameters of the anti-rotation plate as test factors, and generate an encoding table containing various anti-rotation plate geometric parameter combinations based on a central composite design test method.

[0026] The evaluation index acquisition module is configured to generate a plurality of anti-rotation plate sealing structures based on the encoding table, and perform numerical simulation of the steady-state rotor static eccentricity of each anti-rotation plate sealing structure under the same working condition to obtain the index value of the stability evaluation index of each anti-rotation plate sealing structure.

[0027] The sensitivity analysis module is configured to analyze the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the stability evaluation index, and determine a relatively optimal anti-rotation plate design parameter combination according to the analysis result.

[0028] The anti-rotation plate design parameter sensitivity analysis system provided by the embodiment of the present application processes the geometric parameters of the anti-rotation plate by the central composite design method to obtain an encoding table, generates a plurality of anti-rotation plate sealing structures according to the encoding table, obtains the stability evaluation index of the sealing structure by data simulation, and analyzes the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the index, so that the key geometric parameters of the anti-rotation plate for stability enhancement and vibration suppression can be identified and a relatively optimal anti-rotation plate design parameter combination can be obtained through a limited number of tests, which helps turbine designers to quickly and efficiently meet the anti-rotation plate design requirements in the engineering design process.

[0029] In a third aspect, the embodiment of the present application provides a computer device, comprising a memory and a processor, which are communicatively connected with each other, and the memory stores computer instructions, and the processor executes the computer instructions to perform the method of the first aspect or any one of the optional embodiments of the first aspect.

[0030] In a fourth aspect, an embodiment of the present application provides a computer readable storage medium, which stores computer instructions for causing a computer to execute the method in the first aspect, or any optional implementation manner of the first aspect. BRIEF DESCRIPTION OF DRAWINGS

[0031] In order to more clearly illustrate the technical solutions of the embodiments of the present application or the prior art, the drawings needed to be used in the specific embodiments or prior art description will be briefly introduced. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0032] Figure 1 A flowchart of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0033] Figure 2 A static eccentricity diagram of a rotor of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0034] Figure 3 A circumferential velocity contribution rate Pateto diagram of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0035] Figure 4 A circumferential velocity main effect diagram of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0036] Figure 5 A tangential airflow excitation force contribution rate Pateto diagram of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0037] Figure 6 A tangential airflow excitation force main effect diagram of a design parameter sensitivity analysis method of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0038] Figure 7 A structural diagram of a design parameter sensitivity analysis system of an anti-rotation plate provided by an embodiment of the present application is shown in the figure.

[0039] Figure 8 A structural diagram of a computer device provided by an embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0040] In order to make the purposes, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0041] In addition, the technical features involved in the different embodiments of the present application described below can be combined with each other as long as there is no conflict.

[0042] The embodiments of the present application provide a design parameter sensitivity analysis method for anti-rotation plates, as shown in the figure, which specifically comprises the following steps: Figure 1

[0043] Step S1: Selecting k geometric parameters of the anti-rotation plate as test factors, and generating an encoding table containing various combinations of the geometric parameters of the anti-rotation plate based on a central composite design test method.

[0044] Specifically, in the embodiments of the present application, a labyrinth seal anti-rotation plate is used for example, which comes from the labyrinth seal test piece disclosed by Ertas et al. (Ertas BH, Delgado A, Vannini G. Rotordynamic force coefficients for three types of annular gas seals with inlet preswirl and high differential pressure ratio[J]. ASME Journal of Engineering for Gas Turbines and Power, 2012, 134(4): 042503.), and the number of sealing teeth is reduced from 14 to 6. The geometric parameters of the labyrinth seal are shown in the following table:

[0045] Geometric parameter name Value Seal tooth number 6.0 mm Seal gap 0.3 mm Chamber depth 4.01 mm Tooth tip included angle 30° Tooth spacing 5.0 mm Rotor diameter 170.0 mm

[0046] The embodiments of the present application use the DOE method in the parameter optimization of the commercial software Isight for test design, and more specifically, the central composite design method in the DOE method is selected for test design: a horizontal full-factor method enhanced by taking one center point and two additional angular points on the factor axis for each factor.

[0047] First, the factors are selected. In the structure of the labyrinth seal anti-rotation plate, the length (a) of the anti-rotation plate, the number (n) of the anti-rotation plates, the stagger angle (b) of the anti-rotation plates, the number (m) of the anti-rotation plate teeth, the tooth width (c) of the anti-rotation plate, the tooth height (d) of the anti-rotation plate, the tooth thickness (e) of the anti-rotation plate, the tooth tip radius (f) of the anti-rotation plate, the tooth root radius (g) of the anti-rotation plate, the tooth tip chamfer (h) of the anti-rotation plate, and the tooth root chamfer (i) of the anti-rotation plate are selected as the factors. ​Four parameters, including the tilt angle (θ) of the anti-spinning plates, were selected as test factors. Next, a range of values ​​was assigned to each parameter, and five levels were selected. For example, the range of the number of anti-spinning plates was 36-108, with five levels of 36, 54, 72, 90, and 108. Finally, a four-factor, five-level coding table was generated, as shown in the table below:

[0048] Horizontal number Anti-rotation plate length / mm Anti-rotation plate number Anti-rotation plate stagger angle / ° Anti-rotation plate tilt angle / ° 1 3 36 -30 -40 2 4 54 -15 -20 3 (reference) 5 72 0 0 4 6 90 15 20 5 7 108 30 40

[0049] Based on the coding table, 25 experimental schemes were generated using the central composite design method, including 16 four-factor, two-level full-factor experiments (experiment schemes 1-16), 1 central point experiment (experiment scheme 17), and 8 pivot point experiments (experiment schemes 18-25), as shown in the table below:

[0050]

[0051] Step S2: Generate multiple anti-spinning plate sealing structures based on the coding table, and perform a numerical simulation of steady-state rotor static eccentricity for each anti-spinning plate sealing structure under the same working conditions to obtain the index values ​​of the stability evaluation index of each anti-spinning plate sealing structure.

[0052] Specifically, in this embodiment of the invention, 25 anti-swirl plate sealing structures were generated according to the corresponding test schemes in the coding table. Two stability evaluation indices under rotor static eccentricity were obtained using the steady-state solution of the RANS equations in the CFX commercial software: the average circumferential velocity of the downstream airflow of the anti-swirl plate and the tangential airflow excitation force under rotor static eccentricity. The numerical calculations of different anti-swirl plate sealing structures used the same boundary conditions (inlet total pressure 0.69 MPa; inlet total temperature 14°C; inlet pre-swirl ratio 0.45; outlet static pressure 0.383 MPa; rotational speed 15000 rpm).

[0053] Figure 2 A schematic diagram of rotor static eccentricity is given. Rotor static eccentricity refers to the rotor center C not coinciding with the stator center O, and the rotor's vortex velocity around point O being 0 rpm. The computational mesh for the sealed eccentricity with anti-spinning plate structure is obtained by solving the displacement diffusion equation using dynamic meshing techniques based on a dynamic-static concentric mesh. Mesh deformation techniques are an important means of solving moving boundaries in numerical simulations. Mesh deformation techniques first specify the motion of nodes on the boundary, and then the motion of the remaining internal mesh nodes is obtained by solving the displacement diffusion equation, which is expressed as follows:

[0054]

[0055] In the formula: δ is the displacement of the grid node relative to the grid position at the previous moment; Γ dispThe grid stiffness is 0.5. The specific calculation method and the vortex model parameters are shown in the following table:

[0056] Item Value, attribute Solver ANSYS CFX Discrete format High-precision format Solving method Constant value method, dynamic mesh technology Turbulence model Standard k-ε turbulence model, improved wall function method Working medium Ideal air Vortex frequency / Hz 0 Eccentricity amplitude / mm delta = 10% C r ]]

[0057] Step S3: Sensitivity of the anti-rotation plate parameters to the anti-rotation plate design is analyzed based on the stability evaluation index, and a relatively optimal anti-rotation plate design parameter combination is determined according to the analysis result.

[0058] Specifically, in the embodiment of the present application, the value range of the test factor is normalized to [-1, 1], and after normalization, -1 represents a low level of 0%, 1 represents a high level of 100%, 0 represents a 50% level, and so on. Assuming that the maximum of a certain geometric parameter is x max and the minimum is x min , then after normalization:

[0059]

[0060] where i = 1, 2, …, 5, x i represents the parameter value before normalization, x i ' represents the parameter value after normalization.

[0061] All the measured stability evaluation indexes in the design scheme are brought into the commercial software Isight together with the normalized coding table to perform quadratic polynomial fitting. Assuming that the quadratic fitting polynomial of two factors is constructed as follows:

[0062] y = c0 + c1x1 + c2x2 + c3x1 2 + c4x2 2 + c5x1x2

[0063] Then the derivative is:

[0064] dy = c1dx1 + c2dx2 + 2c3x1dx1 + 2c4x2dx2 + c5d(x1x2)

[0065] Therefore, the main effects of the x1 and x2 terms are c1dx1 and c2dx2, respectively; the main effects of the x1 and x2 quadratic terms are 2c3x1dx1 and 2c4x2dx2, respectively; y is a selected response, corresponding to the various stability evaluation indexes in the embodiment of the present application: circumferential velocity or tangential airflow excitation force.

[0066] In the quadratic polynomial fitting process, R 2 is used to represent the fitting accuracy of the polynomial, and R 2 is defined as follows:

[0067]

[0068] where N is the number of sample points; yi is the index value of any stability evaluation index; is the response value after polynomial fitting; y mean is the mean value of response sample points. When R 2 is closer to 1, the fitting is more accurate.

[0069] The average circumferential velocity of the gas flow is fitted with a quadratic polynomial with the normalized encoding table. The average circumferential velocity of the gas flow downstream of the anti-whirl plate is usually used as an indirect index to judge the stability of the sealing rotor, and the smaller the circumferential velocity of the gas flow, the better the stability of the sealing rotor.

[0070] Figure 3 The Pateto diagram based on the average circumferential velocity downstream of the anti-whirl plate under the static eccentricity of the rotor is given, and the fitting accuracy R 2 of the quadratic polynomial fitting model is 0.909. In terms of positive effects, the top three contribution rates are the square of the stagger angle, the square of the length of the anti-whirl plate, and the interaction term (the product of the two) of the stagger angle and the tilt angle.

[0071] Figure 4 The main effect diagram based on the average circumferential velocity downstream of the anti-whirl plate under the static eccentricity of the rotor is given. The average circumferential velocity of the gas flow downstream of the anti-whirl plate gradually decreases with the increase of the tilt angle, indicating that the increase of the tilt angle of the anti-whirl plate is beneficial to enhancing the negative effect of the tilt angle to obtain a lower circumferential velocity of the gas flow. When the stagger angle of the anti-whirl plate increases, the average circumferential velocity of the gas flow downstream of the anti-whirl plate shows a trend of gradually decreasing first and then gradually increasing, and reaches the lowest value at the 50% level (i.e., the front edge of the vertical sealing teeth of the anti-whirl plate is arranged, ), which is because the arrangement of the front edge of the vertical sealing teeth of the anti-whirl plate can eliminate the positive effect of the square of the stagger angle and the positive effect of the interaction term of the stagger angle and the tilt angle. The average circumferential velocity of the gas flow downstream of the anti-whirl plate gradually decreases with the increase of the length of the anti-whirl plate, and reaches the lowest value at the 65% level (a=5, 6 mm), and then gradually increases with the increase of the length of the anti-whirl plate, i.e., the anti-whirl plate is too long and is not conducive to reducing the average circumferential velocity of the gas flow downstream of the anti-whirl plate. When the number of anti-whirl plates is at the level of 30%-74% (n=54, 72, 90), the average circumferential velocity of the gas flow downstream of the anti-whirl plate can achieve a relatively low value (<1 m / s). Based on the above analysis of the sensitivity of the circumferential velocity, the optimal combination of the design parameters of the anti-whirl plate is obtained as shown in the following table:

[0072]

[0073] The tangential excitation force of the rotor is fitted with a quadratic polynomial with the normalized encoding table. Figure 5 The Pateto diagram based on the tangential excitation force of the rotor under the static eccentricity of the rotor is given, and the fitting accuracy R 2is 0.865. When the evaluation index is changed from the average circumferential velocity of the gas flow downstream of the anti-whirl plate to the static eccentric rotor tangential gas flow excitation force, the anti-whirl plate parameter sensitivity changes obviously, and this change is because the circumferential pressure fluctuation size of the upstream anti-whirl plate area is considered in the calculation process of the rotor tangential gas flow excitation force. In terms of positive effects, the top three contribution rates are the square of the stagger angle, the interaction term of the stagger angle and the tilt angle, and the interaction term of the anti-whirl plate length and the number of anti-whirl plates.

[0074] Figure 6 The main effect diagram of the rotor tangential gas flow excitation force under static eccentricity is given, and the rotor tangential gas flow excitation force under static eccentricity gradually decreases with the increase of the tilt angle, indicating that it is also feasible to increase the negative effect by increasing the tilt angle of the anti-whirl plate to obtain a lower tangential gas flow excitation force. When the stagger angle of the anti-whirl plate increases, the rotor tangential gas flow excitation force under static eccentricity also presents a trend of first decreasing and then gradually increasing, and the minimum value is obtained at the 50% level (that is, the anti-whirl plate is arranged vertically in front of the leading edge of the sealing tooth, ). When the anti-whirl plate length or the number of anti-whirl plates increases to 75% level, the rotor tangential gas flow excitation force under static eccentricity gradually decreases; but when the anti-whirl plate length or the number of anti-whirl plates exceeds 75% level, the rotor tangential gas flow excitation force under static eccentricity hardly changes, which is caused by the decrease of the contribution rate of the square term of the anti-whirl plate length and the square term of the number of anti-whirl plates. Based on the above sensitivity analysis of the tangential gas flow excitation force, the optimal anti-whirl plate design parameter combination is as follows:

[0075]

[0076] The anti-whirl plate design parameter sensitivity analysis method provided in the embodiments of the present application can process the geometric parameters of the anti-whirl plate by the central composite design method to obtain an encoding table, generate a plurality of anti-whirl plate sealing structures according to the encoding table, obtain the stability evaluation index of the sealing structure by a data simulation method, and analyze the sensitivity of the anti-whirl plate parameters to the anti-whirl plate design based on the index, so that the key geometric parameters of the anti-whirl plate stability enhancement and vibration suppression can be identified and the optimal anti-whirl plate design parameter combination can be obtained through a limited number of tests, which is helpful for turbine designers to quickly and efficiently meet the anti-whirl plate design requirements in the engineering design process.

[0077] The embodiments of the present application provide an anti-whirl plate design parameter sensitivity analysis system, as shown in Figure 7 , the system comprises:

[0078] The encoding table generation module 1 is used to select k geometric parameters of the anti-whirl plate as test factors, and generate an encoding table containing a plurality of anti-whirl plate geometric parameter combinations based on the central composite design test method. For details, refer to the related description of step S1 in the method embodiment described above, which will not be repeated here.

[0079] The evaluation index obtaining module 2 is configured to generate a plurality of anti-whirl plate sealing structures based on the coding table, and perform numerical simulation of steady-state rotor static eccentricity of each anti-whirl plate sealing structure under the same working condition to obtain an index value of the stability evaluation index of each anti-whirl plate sealing structure. For details, refer to the related description of step S2 in the method embodiment, which will not be repeated here.

[0080] The sensitivity analysis module 3 is configured to analyze the sensitivity of the anti-whirl plate parameters to the anti-whirl plate design based on the stability evaluation index, and determine a relatively optimal combination of anti-whirl plate design parameters according to the analysis result. For details, refer to the related description of step S3 in the method embodiment, which will not be repeated here.

[0081] The anti-whirl plate design parameter sensitivity analysis system provided by the embodiment of the application can obtain a coding table by processing the geometric parameters of the anti-whirl plate through the central composite design method, generate a plurality of anti-whirl plate sealing structures according to the coding table, obtain the stability evaluation index of the sealing structure through data simulation, and analyze the sensitivity of the anti-whirl plate parameters to the anti-whirl plate design based on the index, so as to identify the key geometric parameters of the anti-whirl plate for stability enhancement and vibration suppression and obtain a relatively optimal combination of anti-whirl plate design parameters through a limited number of tests, which is helpful for turbine designers to quickly and efficiently meet the anti-whirl plate design requirements in the engineering design process.

[0082] Figure 8 The structure schematic diagram of the computer device in the embodiment of the application is shown, which includes a processor 901 and a memory 902, wherein the processor 901 and the memory 902 can be connected through a bus or other manners, Figure 8 For example, the connection through the bus is taken as an example.

[0083] The processor 901 can be a central processing unit (CPU). The processor 901 can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or a combination of the above.

[0084] The memory 902, as a non-transitory computer readable storage medium, can be used to store non-transitory software programs, non-transitory computer executable programs and modules, such as program instructions / modules corresponding to the methods in the above method embodiments. The processor 901 performs various functional applications and data processing of the processor by running the non-transitory software programs, instructions and modules stored in the memory 902, that is, implements the methods in the above method embodiments.

[0085] The memory 902 can include a program storage area and a data storage area, wherein the program storage area can store an operating system, application programs required by at least one function; and the data storage area can store data created by the processor 901 and the like. In addition, the memory 902 can include a high-speed random access memory, and can also include a non-transitory memory, such as at least one magnetic disk storage device, a flash memory device, or other non-transitory solid-state memory device. In some embodiments, the memory 902 can optionally include a memory disposed remotely with respect to the processor 901, which can be connected to the processor 901 through a network. Examples of the above network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof.

[0086] One or more modules are stored in the memory 902 and executed by the processor 901 to perform the methods in the above method embodiments.

[0087] The above computer device specific details can be understood in correspondence with the relevant description and effects of the corresponding method embodiments, which will not be described here.

[0088] Those skilled in the art can understand that all or part of the processes in the above method embodiments can be completed by a computer program instructing related hardware, and the implemented program can be stored in a computer readable storage medium. When the program is executed, it can include the processes of the above method embodiments. The storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk drive (HDD) or a solid-state drive (SSD), etc. The storage medium can also include a combination of the above types of memories.

[0089] Although the embodiments of the present application are described in conjunction with the accompanying drawings, various modifications and changes can be made by those skilled in the art without departing from the spirit and scope of the present application, and such modifications and changes fall within the scope defined by the appended claims.

Claims

1. A method for sensitivity analysis of anti-spinning plate design parameters, characterized in that, The method comprises the following steps: Select k The geometric parameters of each anti-spinning plate are used as test factors, and a coding table containing multiple combinations of anti-spinning plate geometric parameters is generated based on the central composite design test method. The coding table includes: 2 k +2 k +1 combination of anti-spin plate geometric parameters, corresponding to 2 k +2 k +1 anti-spinning plate sealing design scheme, which includes: 2 k kind k Two-level full factorial trial protocol, 2 k A test scheme for axis points and a test scheme for center points; Based on the coding table, a plurality of anti-rotation plate sealing structures are generated, and numerical simulation of steady-state rotor static eccentricity is performed on each anti-rotation plate sealing structure under the same working condition to obtain index values of stability evaluation indexes of each anti-rotation plate sealing structure, wherein the stability evaluation indexes comprise average circumferential velocity of anti-rotation plate downstream airflow and rotor tangential airflow excitation force; the process of obtaining the index values of the stability evaluation indexes of each anti-rotation plate sealing structure comprises: using a numerical simulation method of solving RANS equations to calculate corresponding index values of the average circumferential velocity of anti-rotation plate downstream airflow and the rotor tangential airflow excitation force of each anti-rotation plate sealing structure; Based on the stability evaluation indexes, the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design is analyzed, and a better anti-rotation plate design parameter combination is determined according to the analysis result, wherein the process of analyzing the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the stability evaluation indexes comprises: normalizing the value range of the test factor to [-1, 1] to obtain a normalized coding table; performing quadratic polynomial fitting on the normalized coding table and the index values of the stability evaluation indexes of each sealing structure; obtaining the contribution rate and influence effect of each test factor on each stability evaluation index based on the fitting coefficients of the polynomial; and analyzing the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design according to the obtained contribution rate and influence effect.

2. The anti-rotation plate design parameter sensitivity analysis method of claim 1, wherein, The geometric parameters comprise: anti-rotation plate length, anti-rotation plate number, anti-rotation plate stagger angle and anti-rotation plate inclination angle.

3. The anti-rotation plate design parameter sensitivity analysis method of claim 2, wherein, The selection k The process of taking the geometric parameters of the anti-rotation plate as test factors and generating an encoding table containing various combinations of the geometric parameters of the anti-rotation plate based on a central composite design test method includes: Five values are uniformly selected for each of the four geometric parameters in the value range to be analyzed; A four-factor five-level coding table is generated based on the central composite design method, and each item of the coding table gives a set of anti-rotation plate geometric parameter combinations.

4. A system for analyzing sensitivity of anti-rotation plate design parameters, comprising: The method comprises the following steps: The application discloses a method for generating a coding table for a screw protector. k The coding table comprises two k +2 k +1 kinds of screw protector geometric parameter combinations, and two k +2 k +1 kinds of screw protector sealing design schemes. k k The two-level full-factor test scheme, the two shaft point test schemes and the center point test scheme are used as test factors. k The center point test scheme is used as a test factor.​ The evaluation index acquisition module is configured to generate a plurality of anti-rotation plate sealing structures based on a coding table, and perform numerical simulation of steady-state rotor static eccentricity on each anti-rotation plate sealing structure under the same working condition to obtain index values of stability evaluation indexes of each anti-rotation plate sealing structure, wherein the stability evaluation indexes comprise average circumferential velocity of anti-rotation plate downstream airflow and rotor tangential airflow excitation force; the process of obtaining the index values of the stability evaluation indexes of each anti-rotation plate sealing structure comprises: using a numerical simulation method of solving RANS equations to calculate corresponding index values of the average circumferential velocity of anti-rotation plate downstream airflow and the rotor tangential airflow excitation force of each anti-rotation plate sealing structure; The sensitivity analysis module is configured to analyze the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the stability evaluation indexes, and determine a better anti-rotation plate design parameter combination according to the analysis result, wherein the process of analyzing the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design based on the stability evaluation indexes comprises: normalizing the value range of the test factor to [-1, 1] to obtain a normalized coding table; performing quadratic polynomial fitting on the normalized coding table and the index values of the stability evaluation indexes of each sealing structure; obtaining the contribution rate and influence effect of each test factor on each stability evaluation index based on the fitting coefficients of the polynomial; and analyzing the sensitivity of the anti-rotation plate parameters to the anti-rotation plate design according to the obtained contribution rate and influence effect.

5. An electronic device, comprising: The method comprises the following steps: A memory and a processor, which are connected in communication with each other, the memory storing computer instructions, and the processor executing the computer instructions to perform the anti-rotation plate design parameter sensitivity analysis method of any one of claims 1-3.

6. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions for causing the computer to perform the anti-rotation plate design parameter sensitivity analysis method of any one of claims 1-3.

Images