A radiation simulation method, system and computer device for heating turbine blades

Abstract

The present invention discloses a radiation simulation method, system and computer equipment for heating turbine blades, which relates to the field of radiation simulation technology, including the steps of: obtaining the radiation flux of a known radiation receiving surface; ... N A quartz lamp array consisting of 10 lamp tubes, and obtaining the j The relative position information of each measuring point is obtained; the quartz lamp array is matrix-converted, and the quartz lamp array after matrix conversion is processed to obtain a linear equation group; the density of the radiation flux is input into the linear equation group to obtain the power of each lamp in the quartz lamp array; by changing the power of each lamp, the heat flux distribution on the radiation receiving surface is changed to simulate the radiation. The present invention changes the heat flux distribution on the radiation surface by changing the power of the quartz lamp array. The operation is simple and convenient to meet the timeliness and rapidity of the radiation heat flux distribution requirements, and the heat flux distribution on the radiation surface can be changed more quickly as required in the experiment.

Description

Technical Field

[0001] The present invention relates to the field of radiation simulation technology, and in particular to a radiation simulation method, system and computer equipment for heating turbine blades. Background Art

[0002] In recent years, in the heating simulation process of turbine blades or similar small test pieces, aerodynamic heating and thermal radiation heating are the two main methods for thermal performance testing in the aerospace field. The aerodynamic heating method simulates the aerodynamic heating environment through a high-temperature and high-speed wind tunnel and a gas heating device. The electromagnetic induction heating method is a non-contact heating technology that can heat the heated object without contacting it. It has many advantages such as high heating efficiency, low energy consumption, small size, local heating, and no pollution, making it have a very wide range of applications. However, according to existing research, the preparation cycle for tail flame heating of the gas test bench is long, involves many resources, and has a high cost of use. The strong electromagnetic field generated by the sensor during electromagnetic induction heating will interfere with nearby sensors and vibration table systems, which is not conducive to later testing; in addition, the structure and thermal boundary conditions of the turbine blade are very complex, and the shape of the induction coil is difficult to determine, making it difficult to simulate the actual thermal load on the overall area scale of the turbine blade.

[0003] Thermal radiation heating has competitive advantages such as long heating time, high heating power, and multi-zone control. It has become an effective and widely used method for comprehensive thermal performance testing. Thermal radiation heaters include graphite element heaters, solar furnaces, quartz lamps, and xenon lamps. Quartz lamp heaters have the characteristics of high thermal efficiency, long life, small thermal inertia, and easy electronic control. They are very suitable for simulating transient pneumatic heating temperatures. At the same time, they are small in size and high in power, which makes it easy to splice according to the heating shape of different test products, and can ensure a larger heating area while achieving precise control of the local temperature of small-volume specimens.

[0004] The disadvantage of the prior art is that the traditional optimization method using the position of each lamp in the lamp array as the optimization parameter can no longer meet the demand for rapidly changing the radiant heat flux on the surface of the turbine blade. The radiant heat flux of the quartz lamp array on the surface of the turbine blade is related to the position, height, arrangement, power and other parameters of the quartz lamp array. It takes a certain amount of time to adjust the geometric parameters of the lamp array, such as position and height, and the lamp array also needs to be rearranged. In non-uniform aerodynamic thermal tests, it is usually necessary to quickly change the radiant heat flux distribution on the test surface in a short time. The traditional method faces the disadvantages of being time-consuming and labor-intensive. Summary of the invention

[0005] The purpose of the present invention is to address the deficiencies of the above-mentioned prior art and to provide a radiation simulation method, system and computer equipment for heating turbine blades, so as to solve the problem in the prior art that adjusting the geometric parameters of the lamp array, such as position and height, requires a certain amount of time for debugging, and the lamp array also needs to be rearranged. In non-uniform aerodynamic thermal tests, it is usually necessary to quickly change the radiation heat flux distribution of the test surface in a short period of time. Traditional methods face the problems of being time-consuming and labor-intensive.

[0006] The present invention specifically provides the following technical solution: a radiation simulation method for heating turbine blades, comprising the following steps:

[0007] obtaining a radiation flux density of a known radiation receiving surface of a heated turbine blade;

[0008] use N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The relative position information of each measuring point;

[0009] Performing matrix transformation on the quartz lamp array and the relative position information, obtaining the radiation heat flow relationship of the measuring point according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and obtaining a power linear equation group based on the radiation heat flow relationship;

[0010] Inputting the radiation flux density into a power linear equation group to obtain the power of each lamp in the quartz lamp array;

[0011] By adjusting the power of each lamp tube, the heat flux distribution on the radiation receiving surface is changed to perform radiation simulation of heating turbine blades; wherein the heat flux distribution is a manifestation of the radiation flux density.

[0012] Preferably, the step of obtaining the radiation flux density of the known radiation receiving surface of the heating turbine blade comprises the following steps:

[0013] Get the heat power, surface area and solid angle of the light along the propagation direction of the lamp, and get the radiation flux and radiation flux density. The specific expression is:

[0014] .

[0015] in, I l is the radiation flux on the radiation receiving surface; n is the thermal power of the lamp; A n is the surface area of ​​the lamp tube; is the solid angle along the radiation propagation direction, is the radiation flux density; T is the unit time, is the solid angle along the direction of radiation propagation per unit time.

[0016] Preferably, the use N The support lamps are used to construct a quartz lamp array. The specific expression is:

[0017] .

[0018] in, H j The first j The power of the quartz lamp; q ij H j The first j Quartz lamp in the i The heat flow at each measuring point, q ij For the j Support lamp for the first i The relative position information of each measuring point; b i For the i Radiative heat flux at each measuring point.

[0019] Preferably, the step of obtaining each lamp tube for heating the turbine blades j The relative position information of each measuring point is expressed as follows:

[0020] .

[0021] in, q s It is the radiation heat flux of a single row of quartz lamp array received by a certain point on the radiation receiving surface; p is a point on the radiation receiving surface; Q 0 is the thermal power of the lamp; H p The distance between the lamp and the radiation receiving surface p Distance in the vertical direction; L It is the actual heated length of the lamp filament; , α p for p of α value, β p for p of β value, γ p for p of γ value, y The distance from a point on the radiation receiving surface to the lamp tube y The distance of the direction, xThe distance from a point on the radiation receiving surface to the center of the lamp tube x The distance of the direction, N is the number of lamps.

[0022] Preferably, the quartz lamp array and relative position information are transformed into a matrix, and the specific expression is:

[0023] .

[0024] .

[0025] .

[0026] Among them, q is the relative position information matrix, H is the power of the quartz lamp in the quartz lamp array, b is M The determinant of the radiation heat flux at each measuring point.

[0027] Preferably, the radiation heat flow relationship of the measuring point is obtained based on the quartz lamp array matrix obtained after the matrix transformation and the relative position information matrix, and the power linear equation group is obtained based on the radiation heat flow relationship. The specific expression is:

[0028] .

[0029] .

[0030] Among them, B is the coefficient matrix; c is the transpose of b multiplied by q on the left .

[0031] Preferably, if the power linear equations , and the coefficient matrix B is full rank, then the power linear equations have a unique solution, through Get the power of each lamp in the quartz lamp array.

[0032] Preferably, after performing the radiation simulation of heating the turbine blades, the following step is further included: obtaining a relative error based on the distance between the lateral heat flux distribution and each point on the target normal distribution curve, and the relationship between the distance and the target heat flux value. The specific expression is:

[0033] .

[0034] in, is the relative error, is the lateral heat flow distribution, are the coordinates of each point on the target normal distribution curve.

[0035] The present invention provides a radiation simulation system for heating turbine blades, comprising:

[0036] An acquisition module, used for acquiring the radiation flux and radiation flux density of a known radiation receiving surface of a heating turbine blade;

[0037] Location information acquisition module, used to use N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The relative position information of each measuring point;

[0038] A linear equation group acquisition module is used to perform matrix transformation on the quartz lamp array and the relative position information, and obtain the radiation heat flow relationship of the measuring point according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and obtain the power linear equation group based on the radiation heat flow relationship;

[0039] A power acquisition module, used for inputting the radiation flux density into a power linear equation group to obtain the power of each lamp in the quartz lamp array;

[0040] The simulation module is used to change the heat flux distribution on the radiation receiving surface by adjusting the power of each lamp tube, so as to perform radiation simulation of heating turbine blades; wherein the heat flux distribution is a manifestation of the radiation flux density.

[0041] The present invention provides a computer device, comprising a memory and a processor, wherein a program is stored in the memory, and when the program is executed by the processor, the processor executes the steps of a radiation simulation method for heating turbine blades.

[0042] Compared with the prior art, the present invention has the following significant advantages:

[0043] The present invention constructs a quartz lamp array to obtain the relative position information of each lamp tube in the quartz lamp array with respect to the measuring point, and obtains the radiation heat flow relationship of the measuring point through the relative position information in the matrix form of the quartz lamp array, and obtains the power linear equation based on the radiation heat flow relationship. It is not necessary to change or adjust the geometric parameters of the quartz lamp array, and the power of each lamp tube can be obtained through the relationship between the power linear equation and the radiation flux density. Moreover, it is only necessary to change the power of the quartz lamp array to change the heat flow distribution on the radiation surface. The operation is simple and convenient, so as to meet the timeliness and rapidity of the radiation heat flow distribution requirements, and the heat flow distribution on the radiation surface can be changed more quickly as required in the experiment. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 It is a comparison diagram between the target distribution curve and the inverse problem result curve (Monte Carlo) in the present invention;

[0045] Figure 2 is a graph showing the relative error between the solution to the inverse problem and the target heat flow in the present invention;

[0046] Figure 3It is the overall flow chart of the present invention. DETAILED DESCRIPTION

[0047] The following is a clear and complete description of the technical solutions of the embodiments of the present invention in conjunction with the drawings in the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0048] The present invention further explores and improves the method of non-partitioned simulation of non-uniform aerodynamic heating to meet the needs of rapid design. In essence, it is to solve the radiation inverse problem of the power of each lamp in the quartz lamp array with a known heat flux distribution on a test surface. Based on the power fast linear analysis method, the power model of the blade lattice model is analyzed.

[0049] like Figure 3 As shown, an embodiment of the present application provides a radiation simulation method for heating turbine blades, comprising the following steps:

[0050] Step S1: Acquire the radiation flux and radiation flux density of a known radiation receiving surface of a heated turbine blade.

[0051] Get the heat power, surface area and solid angle of the light along the propagation direction of the lamp, and get the radiation flux and radiation flux density. The specific expression is:

[0052] (1)

[0053] in, I l is the radiation flux on the radiation receiving surface; n is the thermal power of the lamp; A n is the surface area of ​​the lamp tube; is the solid angle along the radiation propagation direction, is the radiation flux density; T is the unit time, is the solid angle along the direction of radiation propagation per unit time.

[0054] Step S2: Utilization N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The relative position information of each measuring point.

[0055] Specifically:

[0056] Note that the radiation flux density and power in equation (1) of the radiation flux of the quartz lamp on the radiation test surface are linearly related, so there is an inverse solution to this direct problem, that is, the power of each tube in the quartz lamp array can be solved by knowing the distribution of the radiation flux density of the radiation receiving surface. N The quartz lamp array composed of lamp tubes has a total of M test point data, then based on formula (1):

[0057] (2)

[0058] Among them, in the formula, H j The first j The power of a quartz lamp; q ij H j The first j Quartz lamp in the i The heat flow at each measuring point, q ij For the j Support lamp for the first i The relative position information of each measuring point; b i For the i Radiative heat flux at each measuring point.

[0059] Step S3: Matrix transformation is performed on the quartz lamp array and the relative position information, and the radiation heat flow relationship of the measuring point is obtained according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and the power linear equation group is obtained based on the radiation heat flow relationship.

[0060] Among them, obtain each lamp for the first j The relative position information of the measuring points, q ij Contains j Support lamp for the first i The relative position information of each measuring point is expressed as follows:

[0061] (3)

[0062] Among them, the coefficient q s contains the relative position information of the jth lamp to the i-th measuring point, p is a point on the radiation receiving surface; Q 0 is the thermal power of the lamp; H p The distance between the lamp and the radiation receiving surface p Distance in the vertical direction; L It is the actual heated length of the lamp filament; , α p forp of α value, β p for p of β value, γ p for p of γ value, y The distance from a point on the radiation receiving surface to the lamp tube y The distance of the direction, x The distance from a point on the radiation receiving surface to the center of the lamp tube x The distance of the direction, N is the number of lamps.

[0063] Converting formula (2) into matrix form is:

[0064] (4)

[0065] Where:

[0066] (5)

[0067] (6)

[0068] (7)

[0069] Among them, q is the relative position information matrix, H is the power of the quartz lamp in the quartz lamp array, and b is the determinant of the radiation heat flux at M measuring points.

[0070] Step S4: input the density of the radiation flux into the power linear equation group to obtain the power of each lamp in the quartz lamp array.

[0071] In order to facilitate the solution, we multiply both sides of equation (4) by q T Then we have:

[0072] (8).

[0073] This formula is transformed into the general form of solving the power linear equation system, where B is the coefficient matrix; c is the transpose of b multiplied by q on the left .

[0074] In formula (8), if And the coefficient matrix B is full rank, then the linear equation system has a unique solution, and the power of each lamp in the quartz lamp array can be obtained by solving equation (4).

[0075] Step S5: by adjusting the power of each lamp tube, the heat flux distribution on the radiation receiving surface is changed to perform radiation simulation of heating turbine blades; wherein the heat flux distribution is a manifestation of the radiation flux density.

[0076] Inverse problem analysis of radiation heat flow on three-dimensional blade surface:

[0077] Then, the surface of the blade with slightly complex geometric conditions was tested. Seven lamps with random powers (375, 614, 299.99, 569.99, 190, 290, 439W) were arranged along the leaf height at the leading edge of the back of the blade, and the heat flow at a series of points on the back of the blade was obtained as the ideal heat flow point. Based on the fast linear analysis method for quartz lamp power optimization, the ideal power was basically consistent with the initial setting power, proving the theoretical feasibility of this method.

[0078] Finally, the feasibility of this method in practical application is verified. A Gaussian distribution heat flux is given to verify this method. A lamp array is arranged with five horizontal lamps and two vertical lamps. The height of the horizontal lamp row from the test surface is 40 mm, and the spacing is 24 mm. The height of the horizontal lamp row from the test surface is 80 mm, and the spacing is 22.5 mm. Input into b of formula (4): Figure 1 The target normal distribution curve value in the equation is obtained, and the geometric parameter matrix of each lamp is obtained and then the equation group is solved. The power of the five lamps is 70.14, 209.13, 372, 211.88, 70.84, 49.98, and 50W respectively. The power obtained by the inverse problem analysis algorithm is modeled in Tracepro software, and the heat flux distribution of the test surface obtained by Monte Carlo simulation under this power is obtained. The comparison between the two is as follows: Figure 1 shown.

[0079] To further explore its accuracy, the difference between the lateral heat flux distribution in the numerical simulation method and each point on the target normal distribution curve is divided by the target heat flux value to obtain the relative error:

[0080] (9)

[0081] in, is the relative error, is the lateral heat flow distribution, are the coordinates of each point on the target normal distribution curve.

[0082] The relative error curves of the two are as follows: Figure 2 shown.

[0083] From the comparison in the above figure and the relative error curve, it can be seen that the lamp array power obtained by the power rapid analysis method is substituted into Tracepro for simulation. Due to the overflow of light at the edge of the receiving surface and insufficient superposition, the obtained heat flux distribution has a large error at the edge of the receiving surface, while the overlap with the target heat flux curve in most areas of the center is high, and the relative error is controlled within 5%. It can be considered that the power rapid analysis method can achieve the target Gaussian distribution of heat flux and meet the requirements of speed and accuracy in non-uniform aerodynamic heating tests.

[0084] Based on the above analysis, for the ideal heat flow on the blade surface with slightly complex geometric conditions, the quartz lamp array fast linear analysis method can obtain better power matching, has high engineering applicability, and can realize the rapid design of heat flow distribution.

[0085] Based on the above method, the present invention provides a radiation simulation system for heating turbine blades, including: an acquisition module, a position information acquisition module, a linear equation group acquisition module, a digital power acquisition module and a simulation module.

[0086] The acquisition module is used to obtain the radiation flux and radiation flux density of the known radiation receiving surface of the heating turbine blade; the position information acquisition module is used to use N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The linear equation group acquisition module is used to perform matrix transformation on the quartz lamp array and the relative position information, and obtain the radiation heat flow relationship of the measuring points according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and obtain the power linear equation group based on the radiation heat flow relationship; the power acquisition module is used to input the radiation flux density into the power linear equation group to obtain the power of each lamp in the quartz lamp array; the simulation module is used to change the heat flux distribution on the radiation receiving surface by adjusting the power of each lamp, and perform radiation simulation of heating turbine blades; wherein the heat flux distribution is the embodiment of the radiation flux density.

[0087] The present invention also provides a computer device, including a memory and a processor. The memory stores a program. When the program is executed by the processor, the processor executes the steps of a radiation simulation method for heating turbine blades.

[0088] According to the disclosed embodiments, a computing device may communicate with one or more external devices (e.g., keyboards, pointing devices, Bluetooth communications, etc.), or with any device (e.g., routers, modems, etc.) that enables a computing device to communicate with one or more other computing devices.

[0089] The above content is a further detailed description of the present invention in combination with a specific preferred embodiment. For technicians in the technical field to which the present invention belongs, several simple deductions or substitutions can be made without departing from the concept of the present invention, which should be regarded as belonging to the protection scope of the present invention.

Claims

1. A radiation simulation method for heating turbine blades, characterized in that: The steps include: obtaining a radiation flux density of a known radiation receiving surface of a heated turbine blade; use N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The relative position information of each measuring point; Performing matrix transformation on the quartz lamp array and the relative position information, obtaining the radiation heat flow relationship of the measuring point according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and obtaining a power linear equation group based on the radiation heat flow relationship; Inputting the radiation flux density into a power linear equation group to obtain the power of each lamp in the quartz lamp array; By adjusting the power of each lamp tube, the heat flux distribution on the radiation receiving surface is changed to perform radiation simulation of heating turbine blades; wherein the heat flux distribution is a manifestation of the radiation flux density.

2. A method for simulating the radiation of heating turbine blades according to claim 1, characterized in that: The step of obtaining the radiation flux density of the known radiation receiving surface of the heating turbine blade specifically includes the following steps: Get the heat power, surface area and solid angle of the light along the propagation direction of the lamp, and get the radiation flux and radiation flux density. The specific expression is: in, I l is the radiation flux on the radiation receiving surface; n is the thermal power of the lamp; A n is the surface area of ​​the lamp tube; is the solid angle along the radiation propagation direction, is the radiation flux density; T is the unit time, is the solid angle along the direction of radiation propagation per unit time.

3. A radiation simulation method for heating turbine blades according to claim 1, characterized in that: The utilization N The support lamps are used to construct a quartz lamp array. The specific expression is: in, H j The first j The power of a quartz lamp; q ij H j The first j Quartz lamp in the i The heat flow at each measuring point, q ij For the j Support lamp for the first i The relative position information of each measuring point; b i For the i Radiative heat flux at each measuring point.

4. A radiation simulation method for heating turbine blades according to claim 1, characterized in that: The method of obtaining each lamp tube for heating the turbine blade j The relative position information of each measuring point is expressed as follows: in, q s It is the radiation heat flux of a single row of quartz lamp array received by a certain point on the radiation receiving surface; p is a point on the radiation receiving surface; Q 0 is the thermal power of the lamp; H p The distance between the lamp and the radiation receiving surface p Distance in the vertical direction; L It is the actual heated length of the lamp filament; , α p for p of α value, β p for p of β value, γ p for p of γ value, y The distance from a point on the radiation receiving surface to the lamp tube y The distance of the direction, x The distance from a point on the radiation receiving surface to the center of the lamp tube x The distance of the direction, N is the number of lamps.

5. A method for simulating the radiation of heating turbine blades according to claim 1, characterized in that: The quartz lamp array and relative position information are transformed into a matrix, and the specific expression is: Among them, q is the relative position information matrix, H is the power of the quartz lamp in the quartz lamp array, b is M The determinant of the radiation heat flux at each measuring point.

6. A radiation simulation method for heating turbine blades according to claim 5, characterized in that: The radiation heat flow relationship of the measuring point is obtained according to the quartz lamp array matrix obtained after the matrix transformation and the relative position information matrix, and the power linear equation group is obtained based on the radiation heat flow relationship. The specific expression is: Among them, B is the coefficient matrix; c is the transpose of b multiplied by q on the left .

7. A radiation simulation method for heating turbine blades according to claim 6, characterized in that: If the power linear equations , and the coefficient matrix B is full rank, then the power linear equations have a unique solution, through Get the power of each lamp in the quartz lamp array.

8. A radiation simulation method for heating turbine blades according to claim 1, characterized in that: After the radiation simulation of heating turbine blades is performed, the following steps are further included: obtaining a relative error according to the distance between each point on the lateral heat flux distribution and the target normal distribution curve, and the relationship between the distance and the target heat flux value. The specific expression is: in, is the relative error, is the lateral heat flow distribution, are the coordinates of each point on the target normal distribution curve.

9. A radiation simulation system for heating turbine blades, characterized in that: include: An acquisition module, used for acquiring the radiation flux and radiation flux density of a known radiation receiving surface of a heating turbine blade; Location information acquisition module, used to use N The quartz lamp array is constructed by using the lamp tubes, and the efficiency of each lamp tube for heating the turbine blade is obtained. j The relative position information of each measuring point; A linear equation group acquisition module is used to perform matrix transformation on the quartz lamp array and the relative position information, and obtain the radiation heat flow relationship of the measuring point according to the quartz lamp array matrix and the relative position information matrix obtained after the matrix transformation, and obtain the power linear equation group based on the radiation heat flow relationship; A power acquisition module, used for inputting the radiation flux density into a power linear equation group to obtain the power of each lamp in the quartz lamp array; The simulation module is used to change the heat flux distribution on the radiation receiving surface by adjusting the power of each lamp tube, so as to perform radiation simulation of heating turbine blades; wherein the heat flux distribution is a manifestation of the radiation flux density.

10. A computer device, characterized in that: It comprises a memory and a processor, wherein a program is stored in the memory, and when the program is executed by the processor, the processor executes the steps of a radiation simulation method for heating turbine blades as claimed in any one of claims 1 to 8.

Images