Method for controlling heating temperature of multiple heating rods
By dynamically calculating and updating the heating control amount of multiple heating rods in the annular fuel heating and blasting test of the coolant loss accident, the problem of difficult to control the heating rate of multiple heating rods with high accuracy in the prior art is solved, and the real simulation of the thermodynamic conditions of the coolant loss accident is achieved, supporting technological development and design optimization.
Patent Information
- Application Number
- CN202510358395.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-03-25
- Publication Date
- 2025-06-13
AI Technical Summary
In the ring fuel heating and blasting test of coolant loss accident, it is difficult for the prior art to control the high-precision heating rate of multiple heating rods simultaneously, making it difficult to truly simulate the thermodynamic conditions of coolant loss accidents.
A method for controlling the heating temperature of a plurality of heating rods is provided, by determining the expected heating target temperature and target heating rate of the heating rod, calculating the current temperature and heating rate deviation, and then determining the heating control amount, and dynamically update the heating control amount through temperature closed loop feedback.
The synchronous control of the heating rate of multiple heating rods is achieved, ensuring that the heating rod is heated accurately and stably according to the target heating rate, successfully reproducing the thermodynamic conditions of the coolant loss accident, supporting the development of coolant loss accident technology and the optimization of fuel component design.
Smart Images

Figure CN120152071A_ABST
Abstract
Description
Technical Field
[0001] Embodiments of the present application relate to the technical field of temperature control, and more particularly to a method for controlling the heating temperature of multiple heating rods. Background Art
[0002] The statements herein merely provide background information related to the present application and do not necessarily constitute prior art.
[0003] The coolant loss accident (LOCA) annular fuel heating burst test is a destructive test used to evaluate the bulging and bursting characteristics of nuclear fuel cladding during a coolant loss accident, so as to provide safety data for fuel element design. In the coolant loss accident annular fuel heating burst test, it is necessary to simultaneously implement a constant, stable, and precise heating rate control for multiple heating rods until the fuel element heats up and bursts, so as to reproduce the thermodynamic conditions of the coolant loss accident. Summary of the Invention
[0004] A brief overview of the present application is given below to provide a basic understanding of certain aspects of the present application. It should be understood that this overview is not an exhaustive overview of the present application. It is not intended to identify the key or important parts of the present application, nor is it intended to limit the scope of the present application. Its purpose is merely to present certain concepts in a simplified form as a prelude to a more detailed description to follow.
[0005] Embodiments of the present application provide a method for controlling the heating temperature of multiple heating rods, which includes the following steps: S10: Determine the expected heating target temperature and the target heating rate of multiple heating rods; S20: Determine the current temperature of multiple heating rods and the temperature at the previous moment; S30: Determine the current heating rate of multiple heating rods according to the current temperature and the temperature at the previous moment of multiple heating rods determined in step S20; S40: Determine the heating rate deviation of multiple heating rods according to the target heating rate in step S10 and the current heating rate in step S30; S50: Determine the heating control amount of multiple heating rods according to the heating rate deviation; S60: Heat multiple heating rods according to the heating control amount; S70: Repeat steps S20 - S60 until multiple heating rods are heated to the expected heating target temperature.
[0006] In the embodiments of the present application, according to the heating rate deviation of multiple heating rods, the heating control amount of the multiple heating rods is determined, so as to realize the synchronous control of the heating rates of the multiple heating rods, and ensure that the multiple heating rods can all be heated at the target heating rate. At the same time, by repeating steps S20 - S60, temperature closed-loop feedback can be realized to dynamically update the corresponding heating control amount according to the real-time temperatures of the multiple heating rods, and further ensure that the multiple heating rods can always accurately and stably maintain the target heating rate during the entire heating process. This method makes up for the deficiency in the prior art that the heating rates of multiple heating rods cannot be simultaneously controlled with high precision, successfully reproduces the thermodynamic conditions of the loss-of-coolant accident, and thus provides technical support for the development of the loss-of-coolant accident technology and the optimization of fuel element design. BRIEF DESCRIPTION OF THE DRAWINGS
[0007] Other objects and advantages of the present application will be apparent from the following description of the embodiments of the present application with reference to the accompanying drawings, and will help to provide a comprehensive understanding of the present application.
[0008] Figure 1 is a flowchart of a method for controlling the heating temperature of multiple heating rods according to an embodiment of the present application.
[0009] Figure 2 is a layout schematic diagram of 49 heating rods, and the embodiments of the present application act on the 49 heating rods.
[0010] Figure 3 is a fitting curve graph of the heating rate of the heating rod, and the embodiments of the present application act on multiple heating rods.
[0011] It should be noted that the drawings are not necessarily drawn to scale, but are only shown in a schematic manner that does not affect the understanding of the reader. DETAILED DESCRIPTION OF THE EMBODIMENTS
[0012] In the following, exemplary embodiments of the present application will be described in conjunction with the accompanying drawings. For the sake of clarity and conciseness, not all features of the actual implementation manners are described in the specification. However, it should be understood that many implementation-specific decisions must be made during the development of any such actual implementation manner in order to achieve the specific goals of the developer, for example, to comply with those system- and business-related constraints, and these constraints may vary with different implementation manners. In addition, it should also be understood that although the development work may be very complex and time-consuming, for those skilled in the art who benefit from the content of the present application, such development work is only a routine task.
[0013] Here, it should also be noted that, in order to avoid obscuring the present application due to unnecessary details, only the device structures and / or processing steps closely related to the solution according to the present application are shown in the drawings, while other details less related to the present application are omitted.
[0014] The inventors of the present application have found that in the coolant loss accident annular fuel heating explosion test, it is necessary to simultaneously implement a constant, stable, and precise heating rate control for multiple heating rods to truly simulate the thermodynamic conditions of the coolant loss accident. However, due to the lack of a method in the prior art that can simultaneously perform high-precision heating rate control for multiple heating rods, it is difficult to truly simulate the thermodynamic conditions of the coolant loss accident, thereby restricting the development of the coolant loss accident technology and the optimization of fuel element design.
[0015] Based on this, an embodiment of the present application provides a method for controlling the heating temperature of multiple heating rods. Figure 1 It is a flowchart of the method for controlling the heating temperature of multiple heating rods according to an embodiment of the present application. As Figure 1 shown, it includes the following steps: S10: Determine the expected heating target temperature and the target heating rate of multiple heating rods; S20: Determine the current temperature of multiple heating rods and the temperature at the previous moment; S30: Determine the current heating rate of multiple heating rods according to the current temperature and the temperature at the previous moment of multiple heating rods determined in step S20; S40: Determine the heating rate deviation of multiple heating rods according to the target heating rate in step S10 and the current heating rate in step S30; S50: Determine the heating control amount of multiple heating rods according to the heating rate deviation; S60: Heat multiple heating rods according to the heating control amount; S70: Repeat steps S20 - S60 until multiple heating rods are heated to the expected heating target temperature.
[0016] In an embodiment of the present application, according to the heating rate deviation of multiple heating rods, the heating control amount of multiple heating rods is determined, thereby realizing synchronous control of the heating rates of multiple heating rods and ensuring that multiple heating rods can all be heated at the target heating rate. At the same time, by repeating steps S20 - S60, temperature closed-loop feedback can be realized to dynamically update the corresponding heating control amount according to the real-time temperatures of multiple heating rods, and further ensure that multiple heating rods can always accurately and stably maintain the target heating rate throughout the heating process. This method makes up for the deficiency in the prior art that the heating rates of multiple heating rods cannot be simultaneously controlled with high precision, successfully reproduces the thermodynamic conditions of the coolant loss accident, and thus provides technical support for the development of the coolant loss accident technology and the optimization of fuel element design.
[0017] In some embodiments, the power of the heating rod can be 13KW.
[0018] In some embodiments, the step S50 further includes the following steps: S51: determining the adjusted temperatures of a plurality of heating rods according to the heating rate deviation; S52: converting the adjusted temperatures of the plurality of heating rods into adjusted voltages of the plurality of heating rods; S53: determining the current voltages of the plurality of heating rods; S54: determining the heating control amounts of the plurality of heating rods according to the adjusted voltages in step S52 and the current voltages in step S53.
[0019] In the embodiments of the present application, by converting the adjusted temperatures of a plurality of heating rods into corresponding adjusted voltages, further comparing the current voltages and the adjusted voltages, and determining the amount of the voltage actually controlled according to the comparison result of the two, thus, the control of the voltage not only considers the amount to be adjusted but also considers the current voltage value, realizing the feedback control of the voltage. In the foregoing steps, the feedback control of the temperature is realized, thus realizing the temperature-voltage dual feedback control of a plurality of heating rods. Moreover, since the adjusted temperatures determined in the temperature feedback process can accurately reflect the heating requirements of the plurality of heating rods, and the voltage control amounts of the actuators for heating the heating rods are determined in the voltage feedback process, the temperature and heating rate of the plurality of heating rods can be accurately controlled, thereby ensuring that the plurality of heating rods can accurately and stably rise to the expected heating target temperature at the target heating rate.
[0020] In some embodiments, the adjusted temperatures of the plurality of heating rods are realized by adjusting the voltage, that is, converting the adjusted temperature into the adjusted voltage.
[0021] In some embodiments, the temperature feedback proportional integral derivative (PID) control method can be used to determine the adjusted temperatures of the plurality of heating rods in step S51, and the voltage feedback proportional integral derivative control method can be used to determine the heating control amounts of the plurality of heating rods in step S54.
[0022] In some embodiments, the step S51 further includes the following steps: S511: Determine the proportional parameter of temperature regulation, the integral parameter of temperature regulation, and the derivative parameter of temperature regulation according to the heating rate deviation; S512: Determine the number of adjacent heating rods of each heating rod; S513: Determine the current temperatures of each heating rod and its adjacent heating rods according to the current temperatures of multiple heating rods; S514: Determine the current heating rates of each heating rod and its adjacent heating rods according to the current heating rates of multiple heating rods; S515: Determine the ratio of the heating rate of each heating rod and its adjacent heating rods to the target heating rate according to the current rates of each heating rod and its adjacent heating rods in step S514 and the target heating rate in step S10; S516: Adjust the proportional parameter of temperature regulation, the integral parameter of temperature regulation, and the derivative parameter of temperature regulation according to the number in step S512, the current temperatures of each heating rod and its adjacent heating rods in step S513, the current rates of each heating rod and its adjacent heating rods in step S514, and the ratio in step S515, and obtain the adjusted proportional parameter of temperature regulation, the adjusted integral parameter of temperature regulation, and the adjusted derivative parameter of temperature regulation; S517: Determine the adjusted temperatures of multiple heating rods according to the adjusted proportional parameter of temperature regulation, the adjusted integral parameter of temperature regulation, the adjusted derivative parameter of temperature regulation in step S516, the heating rate deviation in step S40, and the current temperatures of multiple heating rods in step S20.
[0023] In the embodiments of the present application, first, by determining the number of adjacent heating rods of each heating rod and combining the current temperature and current heating rate of itself and its adjacent heating rods, the influence of adjacent heating rods on the temperature and heating rate of this heating rod is fully considered. At the same time, according to the ratio of the heating rate of each heating rod and its adjacent heating rods to the target heating rate, the proportional parameter, integral parameter, and derivative parameter of temperature regulation are adjusted, thereby effectively improving the accuracy of the proportional parameter, integral parameter, and derivative parameter of temperature regulation, and further improving the control accuracy of the adjusted temperatures of multiple heating rods, ensuring that multiple heating rods can reach the target heating rate stably and accurately.
[0024] In some embodiments, the adjusted temperatures of multiple heating rods can be determined according to the adjusted proportional parameter, integral parameter, derivative parameter, heating rate deviation, and current temperature of multiple heating rods, and in combination with the PID operation function in the temperature feedback proportional integral derivative control method.
[0025] In some embodiments, the proportional parameter, integral parameter, and derivative parameter of temperature regulation before adjustment can be determined according to the heating rate deviation of multiple heating rods and in combination with the experience of those skilled in the art.
[0026] In some embodiments, the proportional parameter, integral parameter, and derivative parameter of the temperature regulation before adjustment can be the proportional parameter, integral parameter, and derivative parameter of the temperature regulation in the previous cycle of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation after adjustment.
[0027] In some embodiments, the cycle can be a proportional integral derivative control cycle. Exemplarily, the cycle length of the proportional integral derivative control cycle can be set to 1 s.
[0028] In some embodiments, the proportional parameter, integral parameter, and derivative parameter of the temperature regulation after adjustment have a strong coupling with the current heating rate and current temperature of each heating rod.
[0029] In some embodiments, the action intervals of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation can be divided according to the decision of the fuzzy control algorithm. And the proportional parameter, integral parameter, and derivative parameter of the temperature regulation are adjusted accordingly according to the interval where they are located, so as to accurately determine the adjustment amplitude of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation. While improving the adjustment accuracy and adjustment efficiency of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation, it effectively avoids the oscillation phenomenon caused by the overshoot of the temperature feedback system, and significantly improves the stability of the temperature feedback system. Exemplarily, the action intervals can be divided into an allowable zone, a fine-tuning zone, and an overshoot zone.
[0030] In some embodiments, when the proportional parameter, integral parameter, and derivative parameter of the temperature regulation enter the allowable zone, there is no need to adjust the proportional parameter, integral parameter, and derivative parameter of the temperature regulation. When the proportional parameter, integral parameter, and derivative parameter of the temperature regulation enter the fine-tuning zone, it is necessary to fine-tune the proportional parameter, integral parameter, and derivative parameter of the temperature regulation. When the proportional parameter, integral parameter, and derivative parameter of the temperature regulation enter the overshoot zone, it is necessary to strengthen the adjustment of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation.
[0031] In some embodiments, in step S516: when the ratio of the heating rate of each heating rod in step S515 and its adjacent heating rod to the target heating rate is less than 10%, the proportional parameter of the temperature regulation before adjustment, the integral parameter of the temperature regulation before adjustment, and the derivative parameter of the temperature regulation before adjustment and the proportional parameter of the temperature regulation after adjustment, the integral parameter of the temperature regulation after adjustment, and the derivative parameter of the temperature regulation after adjustment satisfy the following relational expressions:
[0032] P i△n =P i△o .
[0033] I i△n =I i△o .
[0034] D i△n = D i△o 。
[0035] Among them, P i△o represents the proportional parameter of temperature regulation before adjustment; I i△o represents the integral parameter of temperature regulation before adjustment; D i△o represents the differential parameter of temperature regulation before adjustment; P i△n represents the proportional parameter of temperature regulation after adjustment; I i△n represents the integral parameter of temperature regulation after adjustment; D i△n represents the differential parameter of temperature regulation after adjustment.
[0036] In the embodiments of the present application, when the ratio of the heating rate of each heating rod and its adjacent heating rod to the target heating rate is less than 10%, there is no need to adjust the proportional parameter, integral parameter, and differential parameter of temperature regulation, thereby while ensuring the adjustment accuracy of the proportional parameter, integral parameter, and differential parameter of temperature regulation, streamlining unnecessary steps and improving the adjustment efficiency of the proportional parameter, integral parameter, and differential parameter of temperature regulation.
[0037] In some embodiments, in step S516: when there is a ratio greater than 10% and less than 20% among the ratios of the heating rates of each heating rod and its adjacent heating rod in step S515 to the target heating rate, the proportional parameter of temperature regulation before adjustment, the integral parameter of temperature regulation before adjustment, and the differential parameter of temperature regulation before adjustment and the proportional parameter of temperature regulation after adjustment, the integral parameter of temperature regulation after adjustment, and the differential parameter of temperature regulation after adjustment conform to the following relational expressions:
[0038]
[0039] Among them, P i△o represents the proportional parameter of temperature regulation before adjustment; I i△o represents the integral parameter of temperature regulation before adjustment; D i△o represents the differential parameter of temperature regulation before adjustment; P i△n represents the proportional parameter of temperature regulation after adjustment; I i△n represents the integral parameter of temperature regulation after adjustment; D i△n represents the differential parameter of temperature regulation after adjustment; α represents the temperature difference fine-tuning factor coefficient, determined by experiments; β represents the temperature change rate difference fine-tuning factor coefficient, determined by experiments; T i represents the current temperature of the i-th heating rod; △T i and △T i / s iAll represent the adjustment parameters of the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment before adjustment of the i-th heating rod.
[0040] In the embodiments of the present application, when the ratio of the heating rate of each heating rod and its adjacent heating rod to the target heating rate enters the range of 10%-20%, the fine-tuning factor coefficient of the temperature difference between adjacent heating rods and the fine-tuning factor coefficient of the difference in temperature change rate are introduced, and the fuzzy control is used for the fine-tuning of the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment for the fine-tuning factor coefficient of the temperature difference between adjacent heating rods and the fine-tuning factor coefficient of the difference in temperature change rate, which significantly improves the adjustment accuracy and reliability of the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment, and ensures the stability and control accuracy of the temperature feedback system.
[0041] In some embodiments, the fine-tuning factor coefficient of the temperature difference and the fine-tuning factor coefficient of the difference in temperature change rate can be determined according to the actual arrangement and distribution of multiple heating rods to consider the influence of the heating rod distribution on the temperature distribution and heating rate distribution.
[0042] In some embodiments, in step S516: when the ratio of the heating rate of each heating rod and its adjacent heating rod in step S515 to the target heating rate is greater than 20%, the proportional parameter of the temperature adjustment before adjustment, the integral parameter of the temperature adjustment before adjustment, and the derivative parameter of the temperature adjustment before adjustment and the proportional parameter of the temperature adjustment after adjustment, the integral parameter of the temperature adjustment after adjustment, and the derivative parameter of the temperature adjustment after adjustment satisfy the following relationship:
[0043]
[0044] Among them, P i△o represents the proportional parameter of the temperature adjustment before adjustment; I i△o represents the integral parameter of the temperature adjustment before adjustment; D i△o represents the derivative parameter of the temperature adjustment before adjustment; P i△n represents the proportional parameter of the temperature adjustment after adjustment; I i△n represents the integral parameter of the temperature adjustment after adjustment; D i△n represents the derivative parameter of the temperature adjustment after adjustment; γ represents the overshoot factor coefficient of the temperature difference, determined by experiments; δ represents the overshoot factor coefficient of the difference in temperature change rate, determined by experiments; ε is the oscillation coefficient of the adjustment system, determined by experiments; T i represents the current temperature of the i-th heating rod; △T i and △T i / s iAll represent the adjustment parameters of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation before adjustment of the i-th heating rod.
[0045] In the embodiments of the present application, when the ratio of the heating rate of each heating rod and its adjacent heating rod to the target heating rate enters the range of 20%-100%, a temperature difference overshoot factor coefficient, a temperature change rate difference overshoot factor coefficient, and a regulation system oscillation coefficient are introduced, and fuzzy control is used for the temperature difference overshoot factor coefficient and the temperature change rate difference overshoot factor coefficient to strengthen the adjustment of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation, significantly improving the adjustment accuracy and reliability of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation after adjustment, and ensuring the stability and control accuracy of the temperature feedback system.
[0046] In some embodiments, the temperature difference overshoot factor coefficient and the temperature change rate difference overshoot factor coefficient can be determined according to the actual arrangement and distribution of multiple heating rods, so as to consider the influence of the distribution of heating rods on the temperature distribution and heating rate distribution.
[0047] In some embodiments, the regulation system oscillation coefficient is introduced when the adjustment amplitude of the proportional parameter of the temperature regulation increases, so as to suppress the oscillation of the temperature feedback system, and then realize the rapid stability of the temperature feedback system.
[0048] In some embodiments, in step S516: △T i is determined by the following method:
[0049]
[0050] wherein, △T i represents the adjustment parameter of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation before adjustment of the i-th heating rod; T i△1 、T i△2 and T i△n all represent the current temperature difference between the i-th heating rod and its adjacent heating rods; n represents the number of adjacent heating rods to the i-th heating rod.
[0051] In the embodiments of the present application, by introducing the current temperature difference between a plurality of heating rods and their adjacent heating rods, the adjustment of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation of each heating rod fully considers the influence brought by the temperature difference between this heating rod and the adjacent heating rods. Furthermore, the adjustment accuracy and reliability of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation after adjustment of this heating rod are significantly improved. At the same time, due to fully considering the influence of the temperature difference between adjacent heating rods, the identity when multiple heating rods are heating simultaneously is effectively solved, ensuring that multiple heating rods can all accurately heat continuously at a constant and stable target heating rate until the expected heating target temperature is reached.
[0052] In some embodiments, in step S516: △T i / s i is determined by the following method:
[0053]
[0054] wherein, △T i / s i represents the adjustment parameter of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation before adjustment of the i-th heating rod, and T i / s i△1 、T i / s i△2 and T i / s i△n all represent the current heating rate difference between the i-th heating rod and its adjacent heating rods; n represents the number of adjacent heating rods to the i-th heating rod.
[0055] In the embodiments of the present application, by introducing the current heating rate difference between a plurality of heating rods and their adjacent heating rods, the adjustment of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation of each heating rod fully considers the influence brought by the heating rate difference between this heating rod and the adjacent heating rods. Furthermore, the adjustment accuracy of the proportional parameter, integral parameter, and derivative parameter of the temperature regulation after adjustment of this heating rod is significantly improved. At the same time, due to fully considering the influence of the heating rate difference between adjacent heating rods, the identity when multiple heating rods are heating simultaneously is effectively solved, ensuring that multiple heating rods can all accurately heat continuously at a constant and stable target heating rate until the expected heating target temperature is reached.
[0056] In some embodiments, in step S516: T i△1 、T i△2 、T i△n are determined by the following method:
[0057] T i△1 =Ti+1 -T i 。
[0058] T i△2 =T i -T i-1 。
[0059] T i△n =T i -T n 。
[0060] Among them, T i+1 、T i-1 and T n all represent the current temperature of the heating rods adjacent to the i-th heating rod; T i represents the current temperature of the i-th heating rod.
[0061] In the embodiments of the present application, through the above expressions, the current temperature differences between each heating rod and all its adjacent heating rods can be accurately determined in sequence, thereby providing calculation data for determining the adjustment ratio parameter, integral parameter, and differential parameter of the temperature adjustment before adjustment. Specifically, the (i + 1)-th heating rod and the (i - 1)-th heating rod are both adjacent to the i-th heating rod, and T i+1 represents the current temperature of the (i + 1)-th heating rod, T i-1 represents the current temperature of the (i - 1)-th heating rod, and T n represents the current temperature of the heating rods other than the (i + 1)-th heating rod and the (i - 1)-th heating rod that are adjacent to the i-th heating rod.
[0062] Exemplarily, the heating rods are arranged in a 7*7 array and numbered from left to right and from top to bottom. When i = 18, T i represents the current temperature of the 18th heating rod, and T i+1 can represent the current temperature of the 19th heating rod adjacent to the right of the 18th heating rod, and T i-1 can represent the current temperature of the 17th heating rod adjacent to the left of the 18th heating rod, and T n can represent the current temperature of the 11th heating rod adjacent to the upper side of the 18th heating rod and the current temperature of the 25th heating rod adjacent to the lower side of the 18th heating rod.
[0063] In some embodiments, T i△n can be determined by first performing an integral calculation on the difference in the temperature rise change rate of each heating rod and its adjacent heating rods, and then taking the average of the integral results.
[0064] In some embodiments, T i△n can be determined by the following method:
[0065]
[0066] T i△n =T iE / 10.
[0067] Among them, T iE represents the integral of the current temperature difference between the ith heating rod and its adjacent heating rod; T i△n represents the current temperature difference between the ith heating rod and its adjacent heating rod; n represents the number of heating rods adjacent to the ith heating rod; k represents the summation variable; T i -T i-1 Represents the current temperature difference between the i-th heating rod and its adjacent heating rod.
[0068] In some embodiments, in step S516: T i / s i△1 , T i / s i△2 and T i / s i△n Determined by:
[0069] T i / s i△1 =T i+1 / s i+1 -T i / s i .
[0070] T i / s i△2 =T i / s i -T i-1 / s i-1 .
[0071] T i / s i△n =T i / s i -T n / s n .
[0072] Among them, T i+1 / s i+1 , T i-1 / s i-1 and T n / s n Both represent the current heating rate of the heating rod adjacent to the i-th heating rod.
[0073] In the embodiments of the present application, through the above expression, the current heating rate difference between each heating rod and all the heating rods adjacent to it can be accurately determined in sequence, thereby providing calculation data for determining the adjustment parameters of the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment before adjustment. Specifically, the (i + 1)-th heating rod and the (i - 1)-th heating rod are both adjacent to the i-th heating rod, T i+1 / s i+1 represents the current heating rate of the (i + 1)-th heating rod, T i-1 / s i-1 represents the current heating rate of the (i - 1)-th heating rod, T n / s n represents the current heating rate of the heating rods other than the (i + 1)-th heating rod and the (i - 1)-th heating rod that are adjacent to the i-th heating rod.
[0074] Exemplarily, the heating rods are arranged in a 7×7 array and numbered from left to right and top to bottom. When i = 18, T i / s i represents the current heating rate of the 18th heating rod, T i+1 / s i+1 can represent the current heating rate of the 19th heating rod adjacent to the right of the 18th heating rod, T i-1 / s i-1 can represent the current heating rate of the 17th heating rod adjacent to the left of the 18th heating rod, T n / s n can represent the current heating rate of the 11th heating rod adjacent to the top of the 18th heating rod and the current heating rate of the 25th heating rod adjacent to the bottom of the 18th heating rod.
[0075] In some embodiments, T i / s i△n can be determined by first performing an integral calculation on the difference in the rate of change of temperature rise between each heating rod and its adjacent heating rods, and then taking the average value.
[0076] In some embodiments, T i / s i△n can be determined by the following method:
[0077]
[0078] T i / s i△n = T i / s iE / 10.
[0079] Wherein, T i / s iErepresents the integral of the current heating rate difference between the i-th heating rod and its adjacent heating rod, T i / s i△n represents the current heating rate difference between the i-th heating rod and its adjacent heating rod; n represents the number of heating rods adjacent to the i-th heating rod; k represents the summation variable; T i / s i -T i-1 / s i-1 represents the current heating rate difference between the i-th heating rod and its adjacent heating rod.
[0080] In some embodiments, the step S54 further includes the following steps: S541: setting the proportional parameter of voltage regulation, the integral parameter of voltage regulation, and the differential parameter of voltage regulation; S542: determining the heating control amounts of multiple heating rods according to the adjusted voltage in step S52, the current voltage in step S53, and the proportional parameter of voltage regulation, the integral parameter of voltage regulation, and the differential parameter of voltage regulation in step S541.
[0081] In some embodiments, the proportional parameter of voltage regulation, the integral parameter of voltage regulation, and the differential parameter of voltage regulation can be determined according to the experience of those skilled in the art.
[0082] In some embodiments, in step S54: the current voltage, the adjusted voltage, and the heating control amounts of multiple heating rods satisfy the following expression:
[0083] U out =(U ref -U fed )·f(P, I, D).
[0084] Wherein, U ref represents the adjusted voltage of multiple heating rods; U out represents the heating control amounts of multiple heating rods; U fed represents the current voltage of multiple heating rods; P represents the proportional parameter of voltage regulation; I represents the integral parameter of voltage regulation; D represents the differential parameter of voltage regulation; f(P, I, D) represents a function involving P, I, and D.
[0085] In the embodiments of the present application, the above expression defines the relationship between the heating control amounts of multiple heating rods and the current voltage, the adjusted voltage, and the proportional parameter, integral parameter, and differential parameter of voltage regulation of multiple heating rods, so as to ensure that the heating control amounts of multiple heating rods determined by the above expression can accurately reflect the heating requirements of multiple heating rods, and further enable multiple heating rods to all heat up to the expected heating target temperature at a constant, stable, and accurate target heating rate.
[0086] In some embodiments, f(P, I, D) can be the PID operation function in the voltage feedback proportional integral derivative control method.
[0087] The following uses specific examples to further illustrate the process of using the method in this application to control the heating temperature of multiple heating rods.
[0088] Figure 2 is a layout schematic diagram of 49 heating rods, and the embodiments of this application act on these 49 heating rods. As Figure 2 shown, 49 heating rods with a power of 13 KW are arranged in a 7*7 array, and the 49 heating rods are numbered sequentially from left to right and from top to bottom.
[0089] The expected heating target temperature is set to 500 °C, the target heating rate is set to 7 °C / s, and the proportional integral derivative control period is set to 1 s.
[0090] Determine the current temperature of each heating rod, the temperature of the previous moment of each heating rod, and the current heating rate.
[0091] According to the current heating rate and the target heating rate, determine the heating rate deviation of each heating rod, and determine the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment before adjustment according to the heating rate deviation of each heating rod.
[0092] Determine the number of adjacent heating rods of each heating rod. As Figure 2 shown, the number of adjacent heating rods of each heating rod can be 2, 3, or 4. Specifically, the number of adjacent heating rods of the heating rods numbered 1, 7, 43, and 49 is 2; the number of adjacent heating rods of the heating rods numbered 2, 3, 4, 5, 6, 8, 14, 15, 21, 22, 28, 29, 35, 36, 42, 44, 45, 46, 47, and 48 is 3; the number of adjacent heating rods of the heating rods numbered 9, 10, 11, 12, 13, 16, 17, 18, 19, 20, 23, 24, 25, 26, 27, 30, 31, 32, 33, 34, 37, 38, 39, 40, and 41 is 4.
[0093] Determine the current temperature and the current heating rate of each heating rod and its adjacent heating rods, and determine the ratio of the heating rate of each heating rod and its adjacent heating rods to the target heating rate according to the current heating rate of each heating rod and its adjacent heating rods.
[0094] When the ratio of the heating rate of the heating rod and its adjacent heating rods to the target heating rate is less than 10%, there is no need to adjust the proportional parameter, integral parameter, and derivative parameter of the temperature adjustment of this heating rod.
[0095] When there is a ratio greater than 10% and less than 20% among the heating rate ratios of a heating rod and its adjacent heating rods to the target heating rate, the proportional parameter, integral parameter, and derivative parameter of the temperature regulation of this heating rod are finely adjusted according to the following expression to obtain the adjusted proportional parameter, integral parameter, and derivative parameter of the temperature regulation.
[0096]
[0097] When there is a ratio greater than 20% among the heating rate ratios of a heating rod and its adjacent heating rods to the target heating rate, the proportional parameter, integral parameter, and derivative parameter of the temperature regulation of this heating rod are strengthened and adjusted according to the following expression to obtain the adjusted proportional parameter, integral parameter, and derivative parameter of the temperature regulation.
[0098]
[0099] Determine the adjusted temperature of each heating rod according to the current temperature, heating rate deviation, adjusted proportional parameter, integral parameter, and derivative parameter of the temperature regulation of each heating rod.
[0100] Convert the adjusted temperature of each heating rod into an adjusted voltage.
[0101] Determine the proportional parameter, integral parameter, and derivative parameter of the voltage regulation of each heating rod. Determine the current voltage and adjusted voltage of each heating rod.
[0102] According to the current voltage, adjusted voltage, proportional parameter, integral parameter, derivative parameter of the voltage regulation of each heating rod and the following expression, determine the heating control amount of each heating rod.
[0103] U out =(U ref -U fed )·f(P, I, D).
[0104] Repeat the steps from determining the current temperature of the 49 heating rods to determining the heating control amount of each heating rod until the current temperature of the 49 heating rods reaches the set expected heating target temperature of 500 °C.
[0105] The stability and accuracy of controlling the heating temperature of multiple heating rods using the method of this application will be described below.
[0106] According to the current temperature of the 49 heating rods determined above, their corresponding number of repetitions, and the proportional-integral-derivative control period, plot a heating rate curve graph as shown in Figure 3 shown. Figure 3It is a fitting curve graph of the heating rate of the heating rod, and the embodiments of the present application act on the plurality of heating rods. Among them, the abscissa is time (s); the ordinate is temperature (°C).
[0107] As Figure 3 shown, the change trend of the heating rate curve is approximately a straight line slanting upwards. According to the current temperature of the plurality of heating rods, their corresponding number of repetitions, and the proportional integral differential control period, the heating rate curve is fitted to obtain the relationship formula and correlation coefficient between temperature and time: y = 7.1305x + 52.807, R2 = 0.9971. Among them, the independent variable x represents time (s), the dependent variable y represents temperature (°C), and the slope 7.1305 represents the heating rate (°C / s). Since the correlation coefficient R2 is approximately 1, it indicates that the above relationship formula between temperature and time is highly consistent with the heating rate curve, and further indicates that the above relationship formula between temperature and time can accurately reflect the temperature change trend of the plurality of heating rods. From the above relationship formula between temperature and time, it can be seen that the heating rate of 7.1305 °C / s is close to the target heating rate of 7 °C / s, indicating that the method for controlling the heating temperature of the plurality of heating rods provided by the present application can implement a constant, stable, and accurate heating rate control for the plurality of heating rods simultaneously, and enable the plurality of heating rods to be heated to the expected heating target temperature according to the target heating rate.
[0108] For the embodiments of the present application, it should also be noted that, without conflict, the embodiments of the present application and the features in the embodiments can be combined with each other to obtain new embodiments.
[0109] The above is only the specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. The protection scope of the present application shall be subject to the protection scope of the claims.
Claims
1. A method for controlling the heating temperature of a plurality of heating rods, characterized in that: It includes the following steps: S10: Determine the estimated heating target temperature and target heating rate of the plurality of heating rods; S20: Determine the current temperature of the plurality of heating rods and the temperature at the previous moment; S30: determining a current heating rate of the plurality of heating rods according to the current temperature of the plurality of heating rods determined in step S20 and the temperature at the previous moment; S40: determining a heating rate deviation of the plurality of heating rods according to the target heating rate of step S10 and the current heating rate of step S30; S50: determining heating control amounts of the plurality of heating rods according to the temperature rise rate deviation; S60: heating the plurality of heating rods according to the heating control amount; S70: Repeat steps S20-S60 until the plurality of heating rods are heated to the estimated heating target temperature.
2. The method according to claim 1, characterized in that The step S50 also includes the following steps: S51: determining the adjustment temperatures of the plurality of heating rods according to the temperature rise rate deviation; S52: converting the adjusted temperatures of the plurality of heating rods into adjusted voltages of the plurality of heating rods; S53: Determine the current voltages of the plurality of heating rods; S54: Determine the heating control amount of the plurality of heating rods according to the adjusted voltage in step S52 and the current voltage in step S53.
3. The method according to claim 2, characterized in that The step S51 also includes the following steps: S511: determining a proportional parameter of temperature regulation, an integral parameter of temperature regulation, and a differential parameter of temperature regulation according to the heating rate deviation; S512: Determine the number of adjacent heating rods of each heating rod; S513: determining the current temperature of each heating rod and its adjacent heating rods according to the current temperatures of the plurality of heating rods; S514: determining a current heating rate of each heating rod and its adjacent heating rods according to the current heating rates of the plurality of heating rods; S515: determining a ratio of the heating rate of each heating rod and its adjacent heating rod to the target heating rate according to the current rate of each heating rod and its adjacent heating rod in step S514 and the target heating rate in step S10; S516: According to the quantity in step S512, the current temperature of each heating rod and its adjacent heating rods in step S513, the current rate of each heating rod and its adjacent heating rods in step S514, and the ratio in step S515, adjust the proportional parameter of the temperature regulation, the integral parameter of the temperature regulation, and the differential parameter of the temperature regulation, and obtain the adjusted proportional parameter of the temperature regulation, the adjusted integral parameter of the temperature regulation, and the adjusted differential parameter of the temperature regulation; S517: Determine the adjusted temperatures of the multiple heating rods according to the adjusted proportional parameter of the temperature adjustment in step S516, the adjusted integral parameter of the temperature adjustment, the adjusted differential parameter of the temperature adjustment, the heating rate deviation in step S40, and the current temperature of the multiple heating rods in step S20.
4. The method according to claim 3, characterized in that: In step S516: When the ratio of the heating rate of each heating rod and its adjacent heating rod in step S515 to the target heating rate is less than 10%, the proportional parameter of the temperature regulation before adjustment, the integral parameter of the temperature regulation before adjustment, and the differential parameter of the temperature regulation before adjustment and the proportional parameter of the temperature regulation after adjustment, the integral parameter after adjustment, and the differential parameter after adjustment meet the following relationship: P i△n =P i△o ; I i△n =I i△o ; D i△n =D i△o ; Among them, P i△o I represents the proportional parameter of the temperature regulation before the adjustment; i△o represents the integral parameter of the temperature regulation before the adjustment; D i△o represents the differential parameter of the temperature regulation before the adjustment; P i△n I represents the proportional parameter of the adjusted temperature regulation; i△n represents the integral parameter of the adjusted temperature regulation; D i△n represents the differential parameter of the adjusted temperature regulation.
5. The method according to claim 3, characterized in that: In step S516: When the ratio of the heating rate of each heating rod and its adjacent heating rod to the target heating rate in step S515 is greater than 10% and less than 20%, the proportional parameter of the temperature regulation before adjustment, the integral parameter of the temperature regulation before adjustment, and the differential parameter of the temperature regulation before adjustment and the proportional parameter of the temperature regulation after adjustment, the integral parameter of the temperature regulation after adjustment, and the differential parameter of the temperature regulation after adjustment meet the following relationship: Among them, P i△o I represents the proportional parameter of the temperature regulation before the adjustment; i△o represents the integral parameter of the temperature regulation before the adjustment; D i△o represents the differential parameter of the temperature regulation before the adjustment; P i△n I represents the proportional parameter of the adjusted temperature regulation; i△n represents the integral parameter of the adjusted temperature regulation; D i△n represents the differential parameter of the adjusted temperature regulation; α represents the temperature difference fine-tuning factor coefficient, which is determined by experiment; β represents the temperature change rate difference fine-tuning factor coefficient, which is determined by experiment; T i Indicates the current temperature of the i-th heating rod; △T i With △T i / s i All represent adjustment parameters of the proportional parameter of the temperature adjustment before adjustment, the integral parameter of the temperature adjustment before adjustment, and the differential parameter of the temperature adjustment before adjustment of the i-th heating rod.
6. The method according to claim 3, characterized in that In step S516: When the ratio of the heating rate of each heating rod and its adjacent heating rod in step S515 to the target heating rate is greater than 20%, the proportional parameter of the temperature regulation before adjustment, the integral parameter of the temperature regulation before adjustment, and the differential parameter of the temperature regulation before adjustment and the proportional parameter of the temperature regulation after adjustment, the integral parameter of the temperature regulation after adjustment, and the differential parameter of the temperature regulation after adjustment meet the following relationship: Among them, P i△o I represents the proportional parameter of the temperature regulation before the adjustment; i△o represents the integral parameter of the temperature regulation before the adjustment; D i△o represents the differential parameter of the temperature regulation before the adjustment; P i△n I represents the proportional parameter of the adjusted temperature regulation; i△n represents the integral parameter of the adjusted temperature regulation; D i△n represents the differential parameter of the adjusted temperature regulation; γ represents the temperature difference overshoot factor coefficient, which is determined by experiment; δ represents the temperature change rate difference overshoot factor coefficient, which is determined by experiment; ε is the regulation system oscillation coefficient, which is determined by experiment; T i Indicates the current temperature of the i-th heating rod; △T i With △T i / s i All represent adjustment parameters of the proportional parameter of the temperature adjustment before adjustment, the integral parameter of the temperature adjustment before adjustment, and the differential parameter of the temperature adjustment before adjustment of the i-th heating rod.
7. The method according to claim 5 or 6, characterized in that: In step S516: △T i Determined by: Among them, △T i represents the adjustment parameters of the proportional parameter of the temperature adjustment before adjustment, the integral parameter of the temperature adjustment before adjustment and the differential parameter of the temperature adjustment before adjustment of the i-th heating rod; T i△1 、T i△2 and T i△n Both represent the current temperature difference between the i-th heating rod and its adjacent heating rod; n represents the number of heating rods adjacent to the i-th heating rod.
8. The method according to claim 5 or 6, characterized in that: In step S516: △T i / s i Determined by: Among them, △T i / s i represents the adjustment parameter of the proportional parameter of the temperature adjustment before adjustment, the integral parameter of the temperature adjustment before adjustment and the differential parameter of the temperature adjustment before adjustment of the i-th heating rod, T i / s i△1 、T i / s i△2 and T i / s i△n Both represent the current heating rate difference between the i-th heating rod and its adjacent heating rod; n represents the number of heating rods adjacent to the i-th heating rod.
9. The method according to claim 7, characterized in that: In step S516: T i△1 , T i△2 , T i△n Determined by: T i△1 =T i+1 -T i ; T i△2 =T i -T i-1 ; T i△n =T i -T n ; Among them, T i+1 , T i-1 and T n Both represent the current temperature of the heating rod adjacent to the i-th heating rod; T i Indicates the current temperature of the i-th heating rod.
10. The method according to claim 8, characterized in that In step S516: T i / s i△1 , T i / s i△2 and T i / s i△n Determined by: T i / s i△1 =T i+1 / s i+1 -T i / s i ; T i / s i△2 =T i / s i -T i-1 / s i-1 ; T i / s i△n =T i / s i -T n / s n ; Among them, the T i+1 / s i+1 , T i-1 / s i-1 and T n / s n Both represent the current heating rate of the heating rod adjacent to the i-th heating rod.
11. The method according to any one of claims 2 to 10, characterized in that: In step S54: The current voltages, the adjustment voltages, and the heating control amounts of the plurality of heating rods conform to the following expressions: U out =(U ref -U fed )·f(P,I,D); Among them, U ref represents the adjustment voltage of the plurality of heating rods, U out represents the heating control amount of the plurality of heating rods; U fed represents the current voltage of the plurality of heating rods; P represents the proportional parameter of voltage regulation; I represents the integral parameter of voltage regulation; D represents the differential parameter of voltage regulation; f(P, I, D) represents a function involving P, I and D.