Design method for conversion between cold state and hot state of rotor and stator

Abstract

The invention discloses a rotor and stator cold and hot state conversion design method which comprises the following steps: determining a rotor and stator part needing to be optimized, and obtaining a real hot state runner line; establishing a cold-state finite element model of the rotor and stator part, and calculating the cold-state finite element model by using a solid mechanics calculation method; obtaining a calculation thermal state finite element model and a calculation thermal state runner line according to a calculation result; moving each node on the calculation thermal-state flow channel line to the real thermal-state flow channel line to obtain an estimated thermal-state finite element model; and calculating the deviation absolute value of each node corresponding to the estimated thermal state finite element model and the calculated thermal state finite element model, performing threshold value judgment on the deviation absolute value, and when the deviation absolute value of each node is smaller than a threshold value, outputting the cold state finite element model as a final cold state finite element model. According to the method, by means of the mode from the false cold state to the hot state to the true cold state and by means of the design platformization mode, the cold state model of the rotor and stator component can be obtained, and the development period is remarkably shortened.

Description

Technical Field

[0001] The present invention relates to the field of aero-engines, and particularly to a stator-rotor design method. Background Art

[0002] In the process of aero-engine design, after fully considering the differences and influences of the hot and cold flow passages of the stator-rotor, the hot-cold state conversion of the stator-rotor is an essential link. For the existing hot-cold state conversion of aero-engine stator-rotor, since the stator-rotor is a newly designed part, when designing the component performance test piece, existing parts cannot be borrowed, and the costs of time, labor, manufacturing, etc. are too high, and the economic benefits are not high. Summary of the Invention

[0003] An object of the present invention is to provide a design method for the hot-cold state conversion of a stator-rotor.

[0004] The design method for the hot-cold state conversion of the stator-rotor to achieve the above object includes the following steps: S1. Determine the stator-rotor part to be optimized, and obtain the true hot-state flow path line; S2. Establish a cold-state finite element model of the stator-rotor part, apply loads to the cold-state finite element model, assign parameters, and use the computational solid mechanics method to calculate the cold-state finite element model; S3. Obtain a calculated hot-state finite element model and a calculated hot-state flow path line according to the calculation results; S4. Move each node on the calculated hot-state flow path line to the true hot-state flow path line to obtain a predicted hot-state finite element model; S5. Calculate the absolute value of the deviation of the corresponding nodes of the predicted hot-state finite element model and the calculated hot-state finite element model, and perform a threshold judgment on the absolute value of the deviation. When the absolute value of the deviation of each node is less than the threshold, output the cold-state finite element model as the final cold-state finite element model; when there is at least one absolute value of the deviation greater than the threshold, repeat steps S2 to S5, and iterate the cold-state finite element model and the calculated hot-state finite element model until the absolute value of the deviation of each node is less than the threshold.

[0005] In one or more embodiments, the iteration process of step S5 includes the following steps: superimpose the deformation amount of each node in the current hot-state finite element model on each node in the current cold-state finite element model to obtain a transitional cold-state finite element model; use computational solid mechanics to calculate the transitional cold-state finite element model to obtain a transitional calculated hot-state finite element model and a transitional hot-state flow path line; move the transitional hot-state flow path line to the true hot-state flow path line to obtain a transitional predicted hot-state finite element model; calculate the absolute value of the deviation of the corresponding nodes of the transitional predicted hot-state finite element model and the transitional calculated hot-state finite element model.

[0006] In one or more embodiments, in step S4, directly assign the coordinate values of the true hot-state flow path line to each node on the calculated hot-state flow path line.

[0007] In one or more embodiments, the method further includes S6. calculating a cold-state geometric model according to the final cold-state finite element model of the stator-rotor part to be optimized.

[0008] In one or more embodiments, when establishing the cold-state finite element model of the stator-rotor part in step S2, a component model and a flow field region associated with the stator-rotor part are also established.

[0009] In one or more embodiments, in step S2, constraint conditions are applied to the cold-state finite element model and loads are applied.

[0010] In one or more embodiments, the loads include gas force distribution loads, rotational angular velocity loads, and temperature loads.

[0011] In one or more embodiments, the threshold is 0.001 mm or 0.0001 mm.

[0012] In one or more embodiments, the stator-rotor includes one or more of a rotor disk, rotor blades, a vortex reducer, a casing, stator blades, a rectifier, and a shaft.

[0013] Another object of the present invention is to provide a storage medium storing a computer program, and when the computer program runs, it executes the above-mentioned stator-rotor cold and hot state conversion design method.

[0014] The above-mentioned stator-rotor cold and hot state conversion design method can utilize the existing real hot-state flow path line information, and in combination with the known overall engine information, after determining the specific stator-rotor part to be optimized, the initially established cold-state finite element model is used as a false cold-state model for hot-state conversion, and finally a true cold-state model is obtained, and the true hot-state flow path is used for correction. The above method can accurately and quickly obtain the cold-state model, and can borrow the integrated concept of platform design, fully borrow existing components, and on the basis of achieving the purpose of quickly designing the stator-rotor and shortening the development cycle, further reduce costs and improve economic benefits. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The above and other features, properties, and advantages of the present invention will become more apparent from the following description in conjunction with the drawings and embodiments, where:

[0016] Figure 1 is a flowchart of the stator-rotor cold and hot state conversion design method;

[0017] Figure 2 is a schematic diagram of the iterative principle of the transition model;

[0018] Figure 3It is a flowchart of a specific embodiment of the design method for the cold-hot state conversion of the rotor-stator. Detailed implementation mode

[0019] The present invention will be further described below in conjunction with specific embodiments and drawings. More details are set forth in the following description to facilitate a full understanding of the present invention. However, the present invention is obviously capable of being implemented in many other ways different from this description. Those skilled in the art can make similar generalizations and deductions according to the actual application situation without departing from the connotation of the present invention. Therefore, the protection scope of the present invention should not be limited by the content of this specific embodiment.

[0020] It should be noted that these and subsequent other drawings are only examples, and they are not drawn under the condition of equal proportion, and should not be used to limit the actual required protection scope of the present invention.

[0021] The rotor refers to the assembly formed by all rotating parts in the engine, generally including the rotor disk, rotor blades, vortex reducer, etc. The stator refers to the assembly formed by all non-rotating parts in the engine, generally including the casing, stator blades, rectifier, etc.

[0022] Cold-hot state conversion refers to the process of converting the rotor-stator from the cold state to the hot state, or the process of converting the rotor-stator from the hot state to the hot state. The cold-state rotor-stator refers to the rotor-stator in the state under the standard temperature environment, and the hot-state rotor-stator refers to the rotor-stator in the design point state defined by aerodynamics, and can also be considered as the rotor-stator in the working state. In the design stage of the component, accurate cold-state rotor-stator data needs to be given in order to obtain rotor-stator components with qualified dimensions in the subsequent processing link.

[0023] In the traditional cold-hot state conversion method of the rotor-stator, generally only hot-state data can be obtained, and it has relatively high time, labor, and manufacturing costs to convert it into a cold-state model.

[0024] The cold-hot state conversion method of the rotor-stator in this application is as Figure 1As shown, it includes the following steps: S1. Determine the stator-rotor part to be optimized and obtain the real hot-state runner line; S2. Establish a cold-state finite element model of the stator-rotor part, apply loads to the cold-state finite element model, assign parameters, and use the computational solid mechanics method to calculate the cold-state finite element model; S3. Obtain the calculated hot-state finite element model and the calculated hot-state runner line according to the calculation results; S4. Move each node on the calculated hot-state runner line to the real hot-state runner line to obtain the predicted hot-state finite element model; S5. Calculate the absolute value of the deviation of the corresponding nodes of the predicted hot-state finite element model and the calculated hot-state finite element model, and perform a threshold judgment on the absolute value of the deviation. When the absolute value of the deviation of each node is less than the threshold, output the cold-state finite element model as the final cold-state finite element model; when there is at least one absolute value of the deviation greater than the threshold, repeat steps S2 to S5, and iterate the cold-state finite element model and the calculated hot-state finite element model until the absolute value of the deviation of each node is less than the threshold.

[0025] Further, please understand in combination with Figure 2 the schematic diagram shown.

[0026] In step S1, determine the stator-rotor part to be optimized, and obtain other numerous data such as the real hot-state runner line, part material, and part dimensions according to the existing test data or empirical data. For example, if the parts to be optimized are the 7th and 8th stage rotor blades in an eight-stage compressor, then the data of the 7th and 8th stage rotor disks and the data of the 1st to 6th stage rotors can be directly obtained from the database of the existing engine design platform. In this way, some parts can be optimized based on the existing parts.

[0027] In step S2, establish a cold-state finite element model of the stator-rotor part according to the data obtained in S1. For example, establish a cold-state finite element model of the 7th and 8th stage rotor blades, as Figure 2 shown by the cold-state finite element model 310. The cold-state finite element model 310 is the initially established false cold-state model.

[0028] Preferably, when establishing the cold-state finite element model of the stator-rotor part, also establish the component model and the flow field region associated with the stator-rotor part. For example, also establish the 1st to 6th stage rotor models at the same time to achieve global calculation and improve the calculation accuracy.

[0029] Apply to the cold-state finite element model 310 loads including but not limited to gas force distribution loads, rotational angular velocity loads, and temperature loads, and apply constraint conditions. Assign characteristic parameters such as material density and thermal conductivity, carry out finite element analysis, and use the computational solid mechanics method to calculate the cold-state finite element model, as Figure 2 shown in step C.

[0030] Through calculation, the hot-state finite element model in step S3 can be obtained, as Figure 2The shown hot-state finite element model 311, and the calculation results also include the calculation of the hot-state runner line. There are also a series of nodes on the calculated hot-state runner line. Each node of the calculated hot-state finite element model 311 will undergo corresponding deformation.

[0031] Continue with step S4. Move each node on the calculated hot-state runner line to the actual hot-state runner line to obtain the estimated hot-state finite element model 312. For example, as shown in step N, directly assign the coordinate values of the actual hot-state runner line to each node on the calculated hot-state runner line, while the node positions in other locations remain unchanged. The new model obtained is the estimated hot-state finite element model 312.

[0032] Continue with step S5. Calculate the absolute value of the deviation of the corresponding nodes of the estimated hot-state finite element model 312 and the calculated hot-state finite element model 311, and perform a threshold judgment on the absolute value of the deviation, as Figure 2 shown in step M.

[0033] It can be understood that the node numbers between the hot-state models are in a one-to-one correspondence, but the data carried by each node is different. For example, the coordinate of a certain node P in the estimated hot-state finite element model 312 is (X2, Y2), and the coordinate of the corresponding numbered node P' in the calculated hot-state finite element model 311 is (X1, Y1). Subtracting the two, there is an absolute value of deviation for each node. The absolute value of deviation for nodes with the same coordinate data is 0, and the absolute values of deviation for nodes with different data vary.

[0034] When the absolute value of the deviation of each node is less than the threshold, output the cold-state finite element model as the final cold-state finite element model. For example, when the threshold is set to 0.001, when the absolute value of the deviation of each node is less than 0.001, it is considered that the accuracy of the estimated hot-state finite element model 312 and the calculated hot-state finite element model 311 obtained at this time is acceptable, and it is considered that the cold-state finite element model 310 deduced from the calculated hot-state finite element model 311 is acceptable and output as the final cold-state finite element model.

[0035] The output final cold-state finite element model is calculated to obtain the desired cold-state geometric model. This cold-state geometric model can be delivered to the processing link for the actual processing of components to obtain cold-state components with appropriate dimensions. When the cold-state components with these appropriate dimensions are in the working state, they can be deformed into the desired hot-state state as expected to obtain good aerodynamic performance.

[0036] When there is at least one absolute value of deviation greater than the threshold, repeat the above steps to iterate the cold-state finite element model and the calculated hot-state finite element model. For example, when the maximum value of the absolute value of deviation is greater than the threshold, such as when Figure 2 the absolute value of deviation obtained in step M is greater than the threshold, it is necessary to perform Figure 2Perform update iterations on the cold-state finite element model according to the steps A and R shown.

[0037] The iteration process includes the following steps.

[0038] Superimpose the deformation amounts of each node in the current hot-state finite element model onto each node in the current cold-state finite element model to obtain a transitional cold-state finite element model. That is, superimpose Figure 2 the deformation amounts of each node in the calculated hot-state finite element model 311 obtained in the current step shown onto the cold-state finite element model 310 of the current step to obtain a transitional cold-state finite element model 320.

[0039] Use computational solid mechanics to calculate the transitional cold-state finite element model 320, that is, repeat step C, to obtain a transitional calculated hot-state finite element model 321 and a transitional hot-state runner line.

[0040] Move the transitional hot-state runner line to the real hot-state runner line, that is, repeat step N, to obtain a transitional predicted hot-state finite element model 322.

[0041] Calculate the absolute values of the deviations of the corresponding nodes of the transitional predicted hot-state finite element model 322 and the transitional calculated hot-state finite element model 321, that is, repeat step M.

[0042] Compare the absolute values of the deviations between the transitional predicted hot-state finite element model 322 and the transitional calculated hot-state finite element model 321 with the threshold again. If the absolute values of the deviations of each node are all less than the threshold, output the transitional cold-state finite element model 320 in the current state as the final cold-state finite element model; if there is at least one absolute value of the deviation greater than the threshold, continue the iteration.

[0043] The process of continuing the iteration is as follows: Superimpose the deformation amounts of each node in the transitional calculated hot-state finite element model 321 onto the transitional cold-state finite element model 320 to obtain a secondary transitional cold-state finite element model 330. Use computational solid mechanics to calculate the secondary transitional cold-state finite element model 330, that is, repeat step C, to obtain a secondary transitional calculated hot-state finite element model 331 and a secondary transitional hot-state runner line. Move the secondary transitional hot-state runner line to the real hot-state runner line, that is, repeat step N, to obtain a secondary transitional predicted hot-state finite element model 332. Calculate the absolute values of the deviations of the corresponding nodes of the secondary transitional predicted hot-state finite element model 332 and the secondary transitional calculated hot-state finite element model 331, that is, repeat step M.

[0044] Compare the absolute values of the deviations between the secondary transition prediction hot-state finite element model 332 and the secondary transition calculation hot-state finite element model 331 with the threshold value again. If the absolute values of the deviations at each node are all less than the threshold value, the transition cold-state finite element model 330 in the current secondary state is used as the final cold-state finite element model for output, and the iteration is ended as shown in step O; if there is at least one absolute value of the deviation greater than the threshold value, continue the iteration according to the above steps.

[0045] Thus, through the above method, based on the existing components, a hypothetical stator-rotor cold-state finite element model can be obtained. By performing finite element simulation, an initial hypothetical hot-state stator-rotor model can be obtained. After being corrected by the actual hot-state flow path, the initial hot-state stator-rotor model can be obtained. Through finite element iteration, the final stator-rotor cold-state model can be obtained, which is a process of obtaining from a false cold state to a hot state and then to a true cold state. This method can make the best use of the existing component performance test piece components for cold-hot state conversion design, and quickly obtain a suitable cold-state model through the iteration of steps. Especially for the component performance test pieces, core engines, verification engines, etc. of aeroengines, it will greatly reduce the costs of time, labor, machining and manufacturing, etc., and thus improve the economic benefits.

[0046] This application uses specific terms to describe the embodiments of this application. For example, "one embodiment", "an embodiment", and / or "some embodiments" mean a certain feature, structure, or characteristic related to at least one embodiment of this application. Therefore, it should be emphasized and noted that the "one embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more at different positions in this specification does not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of this application can be appropriately combined.

[0047] Although the present invention is disclosed above with preferred embodiments, it is not used to limit the present invention. Any person skilled in the art can make possible changes and modifications without departing from the spirit and scope of the present invention. Therefore, any modification, equivalent change, and modification made to the above embodiments based on the technical essence of the present invention without departing from the content of the technical solution of the present invention all fall within the protection scope defined by the claims of the present invention.

Claims

1. A design method for hot-and-cold state conversion of a rotor and stator, characterized in that: The steps include: S1. Determine the rotor and stator parts that need to be optimized and obtain the real hot flow path line; S2. Establish a cold finite element model of the rotor-stator part, apply load to the cold finite element model, assign parameters, and calculate the cold finite element model using computational solid mechanics methods; S3. Obtaining a calculated hot finite element model and a calculated hot flow path line according to the calculation results; S4. moving each node on the calculated hot state flow path line to the real hot state flow path line to obtain an estimated hot state finite element model; S5. Calculate the absolute value of the deviation of each corresponding node of the estimated hot finite element model and the calculated hot finite element model, perform threshold judgment on the absolute value of the deviation, and when the absolute value of the deviation of each node is less than the threshold, output the cold finite element model as the final cold finite element model; When there is at least one deviation whose absolute value is greater than the threshold, steps S2 to S5 are repeated to iterate the cold finite element model and the calculated hot finite element model until the deviation absolute values ​​of all nodes are less than the threshold.

2. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: The iterative process of step S5 includes the following steps: Superimposing the deformation of each node in the current hot finite element model to each node in the current cold finite element model to obtain a transition cold finite element model; Using computational solid mechanics to calculate the transition cold state finite element model to obtain a transition calculation hot state finite element model and a transition hot state flow path line; Moving the transition hot state flow path line to the real hot state flow path line to obtain a transition estimated hot state finite element model; The absolute value of the deviation of each corresponding node of the transition estimation hot state finite element model and the transition calculation hot state finite element model is calculated.

3. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: In step S4, each node on the calculated hot-state flow path line is directly assigned a coordinate value of the real hot-state flow path line.

4. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: The method further comprises S6. calculating and obtaining a cold geometric model according to the final cold finite element model of the rotor-stator part to be optimized.

5. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: When the cold finite element model of the rotor-stator part is established in step S2, the component model and flow field area associated with the rotor-stator part are also established.

6. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: In step S2, constraints and loads are applied to the cold finite element model.

7. The rotor-stator cold-hot state conversion design method according to claim 6 is characterized in that: The loads include gas force distribution loads, rotational angular velocity loads and temperature loads.

8. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: The threshold value is 0.001 mm or 0.0001 mm.

9. The rotor-stator cold-hot state conversion design method according to claim 1 is characterized in that: The rotor and stator include one or more of a rotor disk, rotor blades, a vortex reducer, a casing, stator blades, a rectifier, and a shaft.

10. A storage medium storing a computer program, characterized in that: When the computer program is running, the rotor-stator cold-hot state conversion design method as described in any one of claims 1-9 is executed.

Images