Three-dimensional quartz lamp temperature control heating virtual test method

Abstract

The invention provides a three-dimensional quartz lamp temperature control heating virtual test method, and the method comprises the steps: building the three-dimensional heat flow distribution of a quartz lamp heater for the radiation heating of the surface of a test piece based on a radiation angle coefficient, and obtaining the nominal heat flow of the quartz lamp heater through a control feedback method; and a quartz lamp heater temperature control heating process in a ground test is simulated, so that efficient virtual thermal test simulation is realized.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of structural heat transfer, and relates to a three-dimensional quartz lamp temperature control heating virtual test method. BACKGROUND

[0002] In ground thermal tests such as aircraft thermal environment simulation tests, static-thermal combined tests and thermal vibration tests, quartz lamp heaters are usually used for temperature control heating according to a temperature control curve. Virtual thermal test simulation is carried out in the test scheme preparation stage, which is an important means to evaluate the feasibility and effectiveness of the test scheme.

[0003] In the ground test using the quartz lamp heater for temperature control heating, the area far from the edge of the heater has a relatively uniform heat flux, which is an effective heating area. However, when there are multiple temperature zones on the surface of the test piece, the area located at the edge or boundary of the temperature zone is inevitably close to the edge of the heater of the temperature zone itself and will be affected by the heater of the adjacent temperature zone. Therefore, in the virtual thermal test simulation, the three-dimensional heat flux distribution of the quartz lamp heater radiating and heating the surface of the test piece must be simulated. If a three-dimensional cavity radiation method is directly used for simulation, the radiation heat flux must be calculated in real time through a high-dimensional sparse matrix in the heat transfer process, and the calculation efficiency is too low. SUMMARY

[0004] The present application aims to at least solve one of the problems in the prior art or related art.

[0005] To this end, the present application provides a three-dimensional quartz lamp temperature control heating virtual test method, and the purpose of the present application is to simulate the three-dimensional heat flux distribution of the quartz lamp temperature control heating in the ground test, and to realize efficient and high-precision virtual thermal test simulation.

[0006] The technical solution of the present application is as follows:

[0007] According to one aspect, a three-dimensional quartz lamp temperature control heating virtual test method is provided, which comprises:

[0008] Step 1, constructing a finite element model of the test piece and each temperature zone heater quartz lamp tube, defining a heating surface of the test piece, and setting an initial temperature;

[0009] Step 2, defining cavity radiation, comprising:

[0010] Taking the finite element model of the test piece and each temperature zone heater quartz lamp tube constructed in step 1 as input, defining the cavity radiation in which the outer surface of each temperature zone heater quartz lamp tube and the heating surface of the test piece participate, and the cavity radiation in which only the heating surface of the test piece participates;

[0011] Step 3, calculating the radiation angle coefficient, comprising:

[0012] Using the finite element models of the test piece and the quartz lamps of the heaters in each temperature zone constructed in step 1, and the cavity radiation groups defined in step 2 as input, calculate the sum of the radiation angle coefficients at each location on the heating surface of the test piece under each cavity radiation group;

[0013] Step 4: Based on the sum of the radiation angle coefficients at each location on the heating surface of the test piece under the radiation of each group of cavities obtained in step 3, calculate the three-dimensional radiation distribution coefficient of the heater in each temperature zone;

[0014] Step 5: Initializing the heating control, including:

[0015] Using the finite element model of the test piece constructed in step 1 as input, select the integration points on the heating surface as the temperature control points of the heaters in each temperature zone, and record the coordinates of the temperature control points. The initial heat flux of the heating surface is set to 0, and the control sensitivity parameter Г of the heaters in each temperature zone is set. j and temperature control curve;

[0016] Step 6: Finite element heat transfer calculation, including:

[0017] Using the finite element model of the test piece constructed in step 1 and the temperature control point determined in step 5 as input, perform a heat transfer calculation for one incremental step based on the current temperature of the test piece and the heat flow boundary of the heating surface, and record the temperature T at each unit of the heating surface. i , Temperature control point temperature of heater in each temperature zone T j and time increment Δt;

[0018] Step 7: Update the nominal heat flow of each zone heater, including:

[0019] The control sensitivity parameter Г set in step 5 j , temperature control curve, and the temperature control point temperature T obtained in step 6 j As input, calculate the heat flow Q of the heater in temperature zone j new，j ;

[0020] Step 8: Based on the obtained three-dimensional radiation distribution coefficient α of the heater in each temperature zone ji And the nominal heat flow Q of the heater in each temperature zone new，j , calculate the actual heat flux at each location on the heating surface;

[0021] Step 9: Using the finite element model of the test piece constructed in step 1 and the actual heat flow at each location on the heating surface as input, update the heat flow boundaries at each location on the heating surface in the finite element model.

[0022] Furthermore, in step 1, the quartz lamp tube is modeled using shell elements, the number of circumferential elements is not less than 8, the thickness is 1 mm and is offset inward.

[0023] Furthermore, the step 4 specifically includes:

[0024] The sum of the radiation angle coefficients of each part of the test piece heating surface under the radiation of each group of cavities is input, and the three-dimensional radiation distribution coefficient of the heater quartz lamp in the temperature zone j at the heating surface unit i is α ji The calculation method is as follows:

[0025] α ji = μ ji - μ i

[0026] Wherein, μ ji is the sum of the radiation angle coefficients of the heating surface unit i under the radiation of the heater quartz lamp in the temperature zone j and the test piece heating surface; μ i is the sum of the radiation angle coefficients of the heating surface unit i under the radiation of the test piece heating surface only.

[0027] Further, the heater heat flow Q new，j of the temperature zone j is calculated by the following formula:

[0028] Q new,j = Q old,j + Γ j (T con,j -T j ) Δt

[0029] Wherein, T con，j is the control temperature of the temperature zone j control curve corresponding to the current time, and Q old，j is the heater heat flow of the temperature zone j before updating.

[0030] Further, in step seven, if Q new，j is less than 0, Q new，j is updated to 0.

[0031] Further, in step eight, the actual heat flow Q i at the heating surface unit i is calculated by the following formula:

[0032]

[0033] Wherein, λ i is the emissivity at the heating surface unit i, and σ is the Boltzmann constant.

[0034] Further, the method further comprises:

[0035] Step ten, if the current time reaches the calculation termination time, the calculation is terminated, and if the current time does not reach the calculation termination time, the next incremental step heat transfer calculation is performed.

[0036] According to another aspect, a computer device is provided, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the method described above when executing the computer program.

[0037] The technical solution described above uses a three-dimensional radiation distribution coefficient based on a radiation angle coefficient to calculate the surface radiation heating of a quartz lamp heater on a test piece, and combines a control feedback method, so that the calculation precision of a virtual thermal test simulation on a region close to the edge of a temperature zone and a boundary region is greatly improved without increasing the calculation amount of finite element heat transfer analysis, and the method can be used for ground thermal test simulation of aircraft thermal environment simulation tests, static-thermal combined tests and thermal vibration tests, and is a new efficient and practical engineering heat transfer simulation method. BRIEF DESCRIPTION OF DRAWINGS

[0038] The accompanying drawings, which constitute a part of the specification, illustrate embodiments of the present application and serve to explain the principles of the present application. It is apparent that the accompanying drawings below only show 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.

[0039] Figure 1 A specific implementation flowchart of the present application is shown;

[0040] Figure 2 A finite element model schematic diagram of a test piece and each temperature zone heater quartz lamp tube according to a specific embodiment of the present application is shown;

[0041] Figure 3 A cavity radiation schematic diagram of an outer surface of a single temperature zone heater quartz lamp tube and a test piece heating surface according to a specific embodiment of the present application is shown;

[0042] Figure 4 A radiation angle coefficient and a schematic diagram of each part of a test piece heating surface under cavity radiation of an outer surface of a single temperature zone heater quartz lamp tube and a test piece heating surface according to a specific embodiment of the present application are shown;

[0043] Figure 5 A test piece surface temperature distribution schematic diagram according to a specific embodiment of the present application is shown.

[0044] Among the above drawings, the following reference signs are included:

[0045] 1, temperature zone 1 quartz lamp tube; 2, temperature zone 2 quartz lamp tube; 3, temperature zone 3 quartz lamp tube; 4, test piece. DETAILED DESCRIPTION

[0046] It should be noted that, in the absence of conflict, the embodiments in this application and the features in the embodiments can be combined with each other. The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. The following description of at least one exemplary embodiment is actually only illustrative and is in no way intended to limit the present invention and its application or use. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0047] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.

[0048] Unless otherwise specifically stated, the relative arrangement of the parts and steps, the numerical expressions and the numerical values ​​set forth in these embodiments do not limit the scope of the present invention. At the same time, it should be understood that, for ease of description, the sizes of the various parts shown in the drawings are not drawn according to the actual proportional relationship. The techniques, methods and equipment known to those of ordinary skill in the relevant art may not be discussed in detail, but where appropriate, the techniques, methods and equipment should be considered as part of the authorization specification. In all examples shown and discussed here, any specific values ​​should be interpreted as being merely exemplary and not as limiting. Therefore, other examples of the exemplary embodiments may have different values. It should be noted that similar numbers and letters represent similar items in the following figures, and therefore, once an item is defined in one figure, it does not need to be further discussed in subsequent figures.

[0049] like Figures 1-5 As shown, in one embodiment of the present invention, a three-dimensional quartz lamp temperature control heating virtual test method is provided, and the test method specifically includes:

[0050] Step 1: Modeling Initialization

[0051] First, construct the finite element model of the test piece and the quartz lamps of the heaters in each temperature zone, define the heating surface of the test piece, and set the initial temperature;

[0052] Preferably, the quartz lamp is modeled using shell elements, the number of circumferential elements should be no less than 8, the thickness should be 1 mm and the elements should be offset inwards.

[0053] For example, the quartz lamp tube 1 of the temperature zone 1, the quartz lamp tube 2 of the temperature zone 2, the quartz lamp tube 3 of the temperature zone 3 and the test piece 4 are shown as follows. Figure 2

[0054] Step two, define the cavity radiation

[0055] With the finite element model input of the test piece and the quartz lamp tube of each temperature zone heater constructed in step one, the cavity radiation participated by the outer surface of the quartz lamp tube of each temperature zone heater and the heating surface of the test piece and the cavity radiation only participated by the heating surface of the test piece are defined respectively;

[0056] Among them, Figure 3 A schematic diagram of the cavity radiation participated by the outer surface of the single temperature zone heater quartz lamp tube and the heating surface of the test piece is shown according to the specific embodiment of the present application.

[0057] Step three, calculate the radiation angle coefficient

[0058] With the finite element model of the test piece and the quartz lamp tube of each temperature zone heater constructed in step one and the groups of cavity radiation defined in step two as input, the radiation angle coefficient sum of each place of the heating surface of the test piece under each group of cavity radiation is calculated.

[0059] Among them, Figure 4 A schematic diagram of the radiation angle coefficient sum of each place of the heating surface of the test piece under the cavity radiation participated by the outer surface of the single temperature zone heater quartz lamp tube and the heating surface of the test piece is shown according to the specific embodiment of the present application.

[0060] In addition, the specific calculation process of the radiation angle coefficient sum of each place of the heating surface of the test piece is a conventional technical means, which will not be described in detail here.

[0061] Step four, calculate the three-dimensional radiation distribution coefficient of each temperature zone heater

[0062] With the radiation angle coefficient sum of each place of the heating surface of the test piece under each group of cavity radiation obtained in step three as input, for the quartz lamp tube of the temperature zone j heater, the three-dimensional radiation distribution coefficient α ji of the test piece heating surface unit i is calculated as follows.

[0063] α ji = μ ji - μ i

[0064] Among them, μ ji is the radiation angle coefficient sum of the heating surface unit i participated by the cavity radiation of the quartz lamp tube of the temperature zone j heater and the heating surface of the test piece; μ i is the radiation angle coefficient sum of the heating surface unit i only participated by the cavity radiation of the heating surface of the test piece.

[0065] Step five, heating control initialization ​

[0066] With the finite element model of the test piece constructed in step one as input, select integral points on the heating surface as the control temperature points of the heaters in each temperature zone, and record the coordinates of the control temperature points, wherein set the initial heat flux of the heating surface to 0, and set the control sensitivity parameters of the heaters in each temperature zone j and the control temperature curve;

[0067] Step six, finite element heat transfer calculation.

[0068] With the finite element model of the test piece constructed in step one and the control temperature points determined in step five as input, perform a heat transfer calculation of an incremental step with the current temperature of the test piece and the heat flux boundary of the heating surface, and record the temperature T i of each element on the heating surface, the temperature T j of the control temperature points of the heaters in each temperature zone, and the time increment At;

[0069] Step seven, update the nominal heat flux of the heaters in each temperature zone.

[0070] With the control sensitivity parameters Γ j , the control temperature curve, and the temperature T j of the control temperature points obtained in step six as input, calculate the heat flux Q new，j of the heaters in temperature zone j as follows:

[0071] Q new,j = Q old,j + Γ j (T con,j -T j )Δt

[0072] Wherein T con，j is the control temperature of the control temperature curve corresponding to the current time in temperature zone j;

[0073] In this step, if Q new，j obtained is less than 0, update Q new，j to 0;

[0074] Step eight, calculate the actual heat flux of each place on the heating surface

[0075] With the three-dimensional radiation distribution coefficient a ji of each heater in each temperature zone obtained in step four and the nominal heat flux Q new，j of each heater in each temperature zone obtained in step seven as input, calculate the actual heat flux Q i of element i on the heating surface as follows:

[0076]

[0077] Wherein λ i is the emissivity of element i on the heating surface, and σ is the Boltzmann constant;

[0078] Step nine, updating the heating surface heat flow boundary

[0079] With the finite element model of the test piece constructed in step one and the actual heat flow of each part of the heating surface obtained in step eight as inputs, the heat flow boundary of each part of the heating surface in the finite element model is updated;

[0080] Step ten, calculation termination determination

[0081] If the current time reaches the calculation termination time, the calculation is terminated, and if the current time does not reach the calculation termination time, the next incremental step heat transfer calculation is performed.

[0082] That is, in consideration of improving the calculation accuracy and efficiency of virtual heat test simulation, the embodiment calculates the three-dimensional radiation distribution coefficient of the radiation heating of the quartz lamp heater on the surface of the test piece based on the radiation angle coefficient, obtains the nominal heat flow of the quartz lamp heater by using the control feedback method, and combines the two to obtain the radiation heating heat flow of the quartz lamp heater on the surface of the test piece, simulates the temperature control heating process of the quartz lamp heater in the ground test, and then realizes efficient virtual heat test simulation.

[0083] By using the above configuration method, the three-dimensional radiation distribution coefficient of the radiation heating of the quartz lamp heater on the surface of the test piece is calculated based on the radiation angle coefficient, and the control feedback method is combined, so that the calculation accuracy of the virtual heat test simulation for the regions close to the edge and the junction of the temperature zone is greatly improved without increasing the calculation amount of the finite element heat transfer analysis, and the method can be used for ground heat test simulation such as aircraft thermal environment simulation test, static-thermal combined test and thermal vibration test, and is a new efficient and practical engineering heat transfer simulation method.

[0084] According to another embodiment, a computer device is provided, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, and the processor implements the above method when executing the computer program.

[0085] Features described and / or shown with respect to one embodiment above can be used in one or more other embodiments in the same or similar manner, and / or in combination with or in place of features in other embodiments.

[0086] It should be emphasized that the term "comprises / comprising" as used herein indicates the presence of the stated features, integers, steps or components, but does not preclude the presence or addition of one or more other features, integers, steps, components or groups thereof.

[0087] The above methods of the present invention can be implemented by hardware, or by a combination of hardware and software. The present invention relates to a computer-readable program that, when executed by a logic component, enables the logic component to implement the apparatus or components described above, or to implement the various methods or steps described above. The present invention also relates to a storage medium for storing the above program, such as a hard disk, magnetic disk, optical disk, DVD, flash memory, etc.

[0088] The many features and advantages of these embodiments are apparent from this detailed description, and thus, the appended claims are intended to cover all such features and advantages of these embodiments that fall within the true spirit and scope thereof. Furthermore, since numerous modifications and changes will readily occur to those skilled in the art, it is not intended that the embodiments of the invention be limited to the exact construction and operation illustrated and described, but rather that all suitable modifications and equivalents be covered within the scope thereof.

[0089] Parts of the present invention that are not described in detail are well known to those skilled in the art.

Claims

1. A three-dimensional quartz lamp temperature control heating virtual test method, characterized in that: The test method includes: Step 1: Construct a finite element model of the test piece and the quartz lamps of the heaters in each temperature zone, define the heating surface of the test piece, and set the initial temperature; Step 2: Define cavity radiation, including: Using the finite element models of the test piece and the quartz lamps of each temperature zone heater constructed in step 1 as input, define the cavity radiation in which both the outer surface of the quartz lamps of each temperature zone heater and the heating surface of the test piece participate, as well as the cavity radiation in which only the heating surface of the test piece participates; Step 3: Calculate the radiation angle coefficient, including: Using the finite element models of the test piece and the quartz lamps of the heaters in each temperature zone constructed in step 1, and the cavity radiation groups defined in step 2 as input, calculate the sum of the radiation angle coefficients at each location on the heating surface of the test piece under each cavity radiation group; Step 4: Based on the sum of the radiation angle coefficients at each location on the heating surface of the test piece under the radiation of each group of cavities obtained in step 3, calculate the three-dimensional radiation distribution coefficient of the heater in each temperature zone; Step 5: Initializing the heating control, including: Using the finite element model of the test piece constructed in step 1 as input, select the integration points on the heating surface as the temperature control points of the heaters in each temperature zone, and record the coordinates of the temperature control points. The initial heat flux of the heating surface is set to 0, and the control sensitivity parameter Г of the heaters in each temperature zone is set. j and temperature control curve; Step 6: Finite element heat transfer calculation, including: Using the finite element model of the test piece constructed in step 1 and the temperature control point determined in step 5 as input, perform a heat transfer calculation for one incremental step based on the current temperature of the test piece and the heat flow boundary of the heating surface, and record the temperature T at each unit of the heating surface. i , Temperature control point temperature of each temperature zone heater T j and time increment Δt; Step 7: Update the nominal heat flow of each zone heater, including: The control sensitivity parameter Г set in step 5 j , temperature control curve, and the temperature control point temperature T obtained in step 6 j As input, calculate the heat flow Q of the heater in temperature zone j new，j ; Step 8: Based on the obtained three-dimensional radiation distribution coefficient α of the heater in each temperature zone ji And the nominal heat flow Q of the heater in each temperature zone new，j , calculate the actual heat flux at each location on the heating surface; Step 9: Using the finite element model of the test piece constructed in step 1 and the actual heat flow at each location on the heating surface as input, update the heat flow boundaries at each location on the heating surface in the finite element model.

2. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 1, characterized in that: In step 1, the quartz lamp tube is modeled using shell elements, with the number of circumferential elements being no less than 8, the thickness being 1 mm, and being offset inward.

3. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 1, characterized in that: The step 4 specifically includes: The sum of the radiation angle coefficients at each location on the heating surface of the test piece under the radiation of each group of cavities is input. For the quartz lamp of the heater in temperature zone j, its three-dimensional radiation distribution coefficient α at unit i on the heating surface of the test piece is ji , which is calculated as follows: a ji =μ ji -m i Among them, μ ji The sum of the radiation angle coefficients at heating surface unit i for the quartz lamp of the heater in temperature zone j and the heating surface of the test piece participating in the cavity radiation; μ i Since only the heating surface of the test piece participates in the cavity radiation, the radiation angle coefficient at the heating surface unit i is and.

4. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 3, characterized in that: The heat flow Q of the heater in temperature zone j is calculated by the following formula new，j : Q new,j =Q old,j +Γ j (T con,j -T j )Δt Among them, T con，j is the control temperature of the temperature zone j corresponding to the current moment, Q old，j To update the heat flow of the heater in the front temperature zone j.

5. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 4, characterized in that: In step 7, if Q is obtained new，j If Q is less than 0, new，j Updated to 0.

6. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 4, characterized in that: In step 8, the actual heat flux Q at the heating surface unit i is i Calculated by the following formula: Among them, λ i is the emissivity of heating surface unit i, and σ is the Boltzmann constant.

7. A three-dimensional quartz lamp temperature control heating virtual test method according to claim 1, characterized in that: The method further comprises: Step 10: If the current moment reaches the calculation termination time, the calculation is terminated. If the current moment does not reach the calculation termination time, return to step 6 and perform the next incremental step heat transfer calculation.

8. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method according to claims 1 to 7 is implemented.

Images