Airplane wheel cooling method and device

Abstract

The invention discloses an aircraft wheel cooling method and device, and the method comprises the steps: collecting the temperature of a brake disc and the air pressure of a hub in real time through a temperature and pressure sensor, calculating the temperature change rate and the pressure difference through an embedded processing module, and generating a thermodynamic initial parameter set; calculating a thermal stress distribution value of the brake disc based on the initial parameter set, and generating a cooling liquid flow and spraying angle instruction set in combination with a material threshold database; the opening degree of an annular pipeline piezoelectric valve is regulated according to an instruction set, an injection area is regulated according to flow velocity-surface area heat exchange efficiency, and waste heat is synchronously recovered to optimize power consumption; the thermal infrared imager obtains temperature field distribution, the initial parameter tolerance range is compared, the spraying angle and the flow instruction are corrected through a fuzzy PID algorithm, and temperature gradient errors are eliminated; and based on the temperature drop rate, the pressure stability and the waste heat recovery efficiency, an efficiency index is generated by using a neural network scoring model, and algorithm weights and piezoelectric valve response parameters are dynamically adjusted. The cooling efficiency is effectively improved, and the thermal stress risk is reduced.

Description

Technical Field

[0001] The present invention relates to aircraft wheel thermal management technology related to the field of aerospace engineering, combining technical fields such as intelligent algorithms in control science and engineering and waste heat recovery technology in energy and power engineering, and in particular to an aircraft wheel cooling method. Background Art

[0002] In modern aviation braking systems, the brake device undergoes a huge process of converting kinetic energy into thermal energy during the landing phase, and its thermal management efficiency directly affects flight safety and component life. Traditional cooling technology mostly uses passive air cooling structure or fixed flow liquid cooling solution, which has three major technical bottlenecks: first, the cooling response lags behind the transient temperature rise, resulting in uneven thermal stress distribution of the brake disc, which is easy to cause micro cracks; second, the fixed injection mode is difficult to adapt to the thermal load distribution of different braking intensities, resulting in low efficiency of coolant utilization; third, the braking energy recovery mechanism is missing, and the overall energy efficiency ratio of the system is insufficient. Existing technologies generally lack multi-physical field collaborative control capabilities, and cannot optimize the balance between cooling parameters and energy recovery in real time, resulting in the brake disc material being in a critical state of thermal fatigue for a long time. Especially for continuous take-off and landing conditions, the coupling effect of air pressure fluctuations and temperature gradients inside the wheel hub will accelerate the aging of seals and increase the probability of tire blowouts. According to statistics from the International Air Transport Association, 21% of unplanned maintenance events for landing gear are caused by thermal management failures. Therefore, it is urgent to develop intelligent thermal control technology to achieve precise and efficient operation of the wheel cooling system through deep integration of multimodal sensing, dynamic fluid control and energy circulation. Summary of the invention

[0003] A method and device for cooling an aircraft wheel, comprising: S1. The temperature and pressure sensors collect brake disc temperature and wheel hub air pressure in real time. The embedded processing module calculates the temperature change rate and pressure difference to generate the initial thermodynamic parameter set.

[0004] S2. Based on the initial parameter set, the thermal stress distribution value of the brake disc is solved, and the cooling instruction set of coolant flow and injection angle is generated in combination with the material threshold database.

[0005] S3. According to the cooling instruction set, the opening value of the piezoelectric valve of the annular pipeline is adjusted, and the injection area is adjusted according to the flow rate-surface area heat exchange efficiency. At the same time, the waste heat is recovered through the thermoelectric conversion device and the system power consumption balance value is optimized.

[0006] S4. The infrared thermal imager obtains the temperature field distribution, compares it with the initial parameter tolerance range, and corrects the injection angle and flow command through the fuzzy PID algorithm to eliminate the temperature gradient error.

[0007] S5. Based on the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process, use a neural network scoring model to generate an efficiency index, and dynamically adjust the algorithm weights and hardware response parameters.

[0008] An aircraft wheel cooling method as described above, wherein the semi-temperature and pressure sensor real-time collects the brake disc temperature and hub air pressure, and the embedded processing module calculates the temperature change rate and pressure difference to generate a set of initial thermodynamic parameters, including the following sub-steps: Perform filtering and denoising processing on the temperature and air pressure data collected by the temperature and pressure sensor, and extract the effective signal segment; Calculate the real-time change rate of the brake disc temperature as the temperature change rate according to the effective signal segment, and calculate the difference between the hub air pressure and the standard air pressure value as the pressure difference; Associate and store the temperature change rate, pressure difference, and time stamp, and generate a set of initial thermodynamic parameters including temperature, air pressure, temperature change rate, and pressure difference.

[0009] An aircraft wheel cooling method as described above, wherein, based on the set of initial parameters, solve the thermal stress distribution value of the brake disc, and combine with the material threshold database to generate a cooling instruction set for the coolant flow rate and injection angle, including the following sub-steps: Based on the set of initial thermodynamic parameters, simulate the thermal stress distribution of the brake disc by the finite element analysis method, and output the stress peak value and gradient value of each region; According to the temperature resistance threshold and compressive resistance threshold in the material threshold database, match the coolant demand level corresponding to the current stress distribution value; Combine the demand level with the heat exchange efficiency model to generate a cooling instruction set including coolant flow rate, injection angle, and duration.

[0010] An aircraft wheel cooling method as described above, wherein, according to the cooling instruction set, adjust the opening value of the piezoelectric valve of the annular pipeline, adjust the injection area according to the flow velocity-surface area heat exchange efficiency, and at the same time recover the waste heat through a thermoelectric conversion device and optimize the system power consumption balance value, including the following sub-steps: According to the flow rate value in the cooling instruction set, adjust the opening of the corresponding piezoelectric valve in the annular pipeline in proportion to control the coolant flow velocity; Based on the flow velocity-surface area heat exchange efficiency model, dynamically adjust the injection area to cover the thermal stress peak area; Convert the waste heat into electric energy through a thermoelectric conversion device, and calculate the system power consumption balance value in real time and optimize the energy distribution strategy.

[0011] An aircraft wheel cooling method as described above, wherein the infrared thermal imager obtains the temperature field distribution, compares the initial parameter tolerance range, and corrects the injection angle and flow rate instruction through the fuzzy PID algorithm to eliminate the temperature gradient error, including the following sub-steps: Obtain the temperature field distribution data of the brake disc through an infrared thermal imager, and extract the actual temperature gradient curve; Compare the actual temperature gradient with the tolerance range of the initial parameters, calculate the deviation amount and input it into the fuzzy PID controller; Adjust the coolant flow rate proportionality coefficient and the injection angle compensation value according to the deviation amount until the temperature gradient error is less than the preset threshold.

[0012] An aircraft wheel cooling method as described above, wherein comparing the actual temperature gradient with the tolerance range of the initial parameters, calculating the deviation amount and inputting it into the fuzzy PID controller includes the following sub-steps: When the temperature gradient deviation amount is greater than the upper tolerance limit, preferentially increase the flow rate proportionality coefficient; When the deviation amount is within the tolerance range but there are local over-limit conditions, preferentially adjust the injection angle compensation value; When the deviation amount continuously falls below the lower tolerance limit, reduce the flow rate proportionality coefficient and enable the waste heat recovery priority.

[0013] An aircraft wheel cooling method as described above, wherein based on the temperature decay rate, pressure stability coefficient and energy recovery efficiency value during the cooling process, using a neural network scoring model to generate an efficiency index, and dynamically adjusting the algorithm weights and hardware response parameters, including the following sub-steps: Real-time collect the temperature decay rate, pressure stability coefficient and energy recovery efficiency value during the cooling process; Input the above parameters into the pre-trained neural network scoring model to output the comprehensive efficiency index; Adjust the proportional-integral-derivative weights of the fuzzy PID algorithm according to the efficiency index, and optimize the response speed of the piezoelectric valve and the thermoelectric conversion efficiency parameters.

[0014] 8. An aircraft wheel cooling device, characterized in that it includes: Data acquisition and preprocessing module: The temperature and pressure sensors collect the brake disc temperature and hub air pressure in real time, and the embedded processing module calculates the temperature change rate and pressure difference to generate a set of initial thermodynamic parameters.

[0015] Thermal stress analysis and instruction generation module: Based on the set of initial parameters, solve the thermal stress distribution value of the brake disc, and combine with the material threshold database to generate a cooling instruction set for the coolant flow rate and injection angle.

[0016] Cooling execution and waste heat recovery module: According to the cooling instruction set, adjust the opening value of the piezoelectric valve in the annular pipeline, adjust the injection area according to the flow rate-surface area heat exchange efficiency, and at the same time recover waste heat through the thermoelectric conversion device and optimize the system power consumption balance value.

[0017] Temperature Field Feedback and Dynamic Correction Module: The infrared thermal imager acquires the temperature field distribution, compares it with the tolerance range of the initial parameters, and corrects the injection angle and flow rate commands through the fuzzy PID algorithm to eliminate the temperature gradient error.

[0018] Efficiency Evaluation and Adaptive Optimization Module: Based on the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process, a neural network scoring model is used to generate an efficiency index, and the algorithm weights and hardware response parameters are dynamically adjusted.

[0019] A computer storage medium, characterized by comprising: at least one memory and at least one processor; The memory is used to store one or more program instructions; The processor is used to run one or more program instructions to execute an aircraft wheel cooling method described in any one of the above.

[0020] The beneficial effects achieved by the present invention are as follows: In the present invention, the brake disc temperature and hub air pressure are real-time monitored by temperature and pressure sensors, and the thermodynamic parameters are dynamically analyzed in combination with the embedded processing module to establish a high-precision initial data set, overcoming the control error caused by data lag in the traditional method; based on the collaborative analysis of the thermal stress distribution calculation and the material threshold database, dynamic instructions for the coolant flow rate and injection angle are generated to achieve the precise adaptation of the cooling intensity and thermal load, optimizing the energy configuration while ensuring the material safety of the brake disc; the piezoelectric valve opening adjustment and dynamic matching technology of the injection area are used, combined with the waste heat recovery mechanism, to construct a cooling and energy regeneration collaborative system, significantly improving the heat dissipation efficiency and reducing the overall energy consumption of the system; through the closed-loop monitoring of the infrared thermal imager and the dynamic correction of the fuzzy PID algorithm, the temperature gradient fluctuation is effectively suppressed, ensuring the temperature uniformity of the brake disc surface and avoiding the deformation risk caused by local overheating; further introducing a neural network model, the cooling efficiency is intelligently evaluated and dynamically optimized by integrating multi-dimensional parameters, enhancing the adaptive ability of the system under extreme working conditions, forming an intelligent cooling system with continuous learning ability, and comprehensively improving the reliability, stability, and energy utilization efficiency of the braking system. Description of the Drawings

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments recorded in the present invention. For those of ordinary skill in the art, other drawings can also be obtained according to these drawings.

[0022] Figure 1 It is a flowchart of an aircraft wheel cooling method provided by an embodiment of the present application.

[0023] Figure 2It is a schematic diagram of an aircraft wheel cooling device provided by an embodiment of the present application. Specific embodiments

[0024] Combined with the accompanying drawings in the embodiments of the present invention below, the technical solutions in the embodiments of the present invention are clearly and completely described. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present invention.

[0025] Embodiment 1 As Figure 1 shown, an aircraft wheel cooling method according to an embodiment of the present application includes: Step S1: The temperature and pressure sensor collects the brake disc temperature and hub air pressure in real time, and the embedded processing module calculates the temperature change rate and pressure difference to generate a set of initial thermodynamic parameters, specifically including the following sub-steps: Step S11: Filter and denoise the temperature and pressure data collected by the temperature and pressure sensor, and extract the effective signal segment; The embedded system first performs composite noise reduction processing on the original data of the temperature and pressure sensor. For the random noise in the temperature signal, moving average filtering is used for time series smoothing, and high-frequency disturbances are suppressed by dynamically adjusting the time window; for the possible sudden outliers in the pressure signal, the median filtering algorithm is enabled to maintain the integrity of the signal waveform and eliminate transient interference terms at the same time. To further improve the signal quality, the system performs multi-resolution decomposition on the preprocessed data, uses specific wavelet basis functions to separate the noise components in different frequency bands, and weakens the residual interference through threshold processing methods. After the noise reduction is completed, the signal enters the effective segment recognition stage. Based on the change point detection algorithm, the statistical characteristics of the signal change trend are analyzed, and the transient interference during the sensor startup stage and the invalid data during system anomalies are excluded by combining the sliding hypothesis testing method. Finally, the stable signal interval with physical significance is locked.

[0026] Step S12: Calculate the real-time change rate of the brake disc temperature as the temperature change rate according to the effective signal segment, and calculate the difference between the hub air pressure and the standard air pressure value as the pressure difference; The calculation of the temperature change rate uses an improved differential and fitting hybrid algorithm. First, perform time series differential operation based on the temperature values of continuous sampling points. When it is detected that the temperature fluctuation exceeds the preset range, it automatically switches to the linear fitting mode: Among them, represents the optimal parameter value obtained by minimizing the objective function; represents finding the value that makes the subsequent expression take the minimum value; Indicates the total number of data points participating in the linear fitting; Indicates the discrete time index within the time window, with values ranging from , and is used to index the data points participating in the fitting; Indicates the time series data; Indicates the intercept term in the linear fitting model; Indicates the slope term in the linear fitting model; Indicates the time interval between adjacent time points.

[0027] Balance the calculation accuracy and real-time requirements by adaptively adjusting the time window range. In the differential pressure calculation section, the system obtains the vehicle load parameters in real time through the vehicle bus, dynamically matches the multi-dimensional pressure calibration curve corresponding to the wheel hub model, and at the same time introduces the air pressure compensation coefficient based on the ambient temperature to correct the original air pressure value to a comparable value under standard working conditions.

[0028] All intermediate calculation results need to pass a two-way verification mechanism: on the one hand, set the theoretical boundary values of parameter changes based on the thermodynamic properties of the material: Among them, Indicates the optimal parameter value obtained by minimizing the objective function; its value range satisfies ; Indicates The theoretical minimum value of; Indicates The theoretical maximum value of; Indicates the thermal conductivity of the material, reflecting the ability of the material to conduct heat; Indicates the area related to heat transfer; Indicates the density of the brake disc material; Indicates the specific heat capacity; Indicates the volume of the brake disc; Indicates the maximum heat.

[0029] On the other hand, construct a dynamic threshold model in combination with historical data, start the parallel review of the backup algorithm for abnormal calculation results, and ensure the physical rationality of parameter output and the reliability of the system.

[0030] Step S13: Correlate and store the temperature change rate, differential pressure, and time stamp to generate a thermodynamic initial parameter set including temperature, air pressure, temperature change rate, and differential pressure; The embedded system establishes a unified time reference through a hardware clock source and adopts an event-triggered mechanism in a multi-threaded processing architecture to achieve precise alignment of sensor data, calculation parameters, and timestamps. The system constructs a multi-dimensional data structure to integrate original observations, derived parameters, and environmental context information. The storage layer uses a combined mechanism of hierarchical buffering and persistent storage to maintain a rolling update of the latest parameter set in the real-time data area, and at the same time transfers structured data batches to non-volatile memory at fixed intervals.

[0031] The data encapsulation process introduces a cyclic redundancy check mechanism to ensure the integrity of information during transmission and establishes a fast retrieval channel through an address mapping table. To support subsequent fault diagnosis and life prediction functions, the system automatically performs data feature extraction and statistic calculation to form an enhanced thermodynamic parameter set containing information such as trend features and extreme value distributions.

[0032] Step S2: Based on the initial parameter set, solve the thermal stress distribution value of the brake disc, and combine it with the material threshold database to generate a cooling instruction set for coolant flow rate and injection angle, which specifically includes the following sub-steps: Step S21: Based on the initial thermodynamic parameter set, simulate the thermal stress distribution of the brake disc through the finite element analysis method, and output the stress peak and gradient values of each region; Construct a three-dimensional transient heat conduction model of the brake disc based on the initial thermodynamic parameter set, and its control equation is: Among them, represents the density of the brake disc material; represents the specific heat; represents the temperature the partial derivative with respect to time ; represents the Hamiltonian operator, which represents the divergence operation of a vector in three-dimensional space; represents the reference value of the thermal conductivity; represents the temperature coefficient; represents the temperature gradient; represents the braking friction coefficient, which reflects the friction characteristics between the friction surfaces during braking; represents the braking pressure; represents the contact speed, which is the relative movement speed between the brake disc and the friction component; represents the contact area.

[0033] Adopt an adaptive mesh refinement technique to locally encrypt the braking contact area and dynamically balance the calculation accuracy and real-time requirements. The setting of boundary conditions integrates real-time temperature field monitoring data and dynamic parameters of wheel hub air pressure, and characterizes the thermal conductivity and thermal expansion characteristics of the brake disc with temperature changes through a non-linear material model.

[0034] During the finite element solution process, an implicit time integration algorithm is introduced to iteratively solve the distribution of the thermo-mechanical coupling field within each calculation step, and a spatial distribution matrix containing equivalent stress, principal stress direction, and gradient change rate is output. For key regions, the system automatically extracts the stress peak and its spatial change rate to generate a two-dimensional thermal stress characteristic map marked in polar coordinates, providing a quantitative basis for subsequent decision-making.

[0035] Step S22: According to the temperature resistance threshold and compressive strength threshold in the material threshold database, match the coolant demand level corresponding to the current stress distribution value; The system has a built-in hierarchical material threshold library, and for each type of brake disc material, a dynamic performance atlas containing parameters such as phase change temperature, yield strength decay curve, and fatigue crack growth rate is established. In the stress matching process, spatial mapping technology is used to compare the stress distribution field output by the finite element with the three-dimensional failure surface in the material database point by point, and the over-threshold risk area is identified through tensor operations. The determination of the cooling demand level introduces a fuzzy logic algorithm, comprehensively considering characteristic parameters such as the over-limit ratio of the stress peak, the steepness of the gradient change, and the area of the high-temperature region, and outputting a decision vector containing multi-level instructions such as emergency cooling, maintaining heat dissipation, and preventive adjustment. At the same time, the system monitors the cumulative damage amount of the material performance in real time. When it is detected that the brake disc has experienced multiple thermal cycles, the aging correction coefficient is automatically loaded to dynamically compensate the threshold, ensuring that the determination logic conforms to the actual service state of the material.

[0036] Step S23: Combine the demand level with the heat exchange efficiency model to generate a cooling instruction set containing coolant flow rate, injection angle, and duration; The heat exchange efficiency model pre-sets the flow field distribution characteristics at different injection angles through computational fluid dynamics simulation, and establishes a non-linear relationship matrix of flow rate - coverage area - heat dissipation rate. The formulation of the cooling strategy uses a constrained optimization algorithm. On the premise of meeting the stress suppression requirements, the optimal combination of flow rate distribution and injection angle is solved with the minimum coolant consumption and the shortest response time as the objective function. For the multi-region collaborative cooling scenario, the system implements hierarchical injection control, preferentially performing high-precision directional cooling on the peak stress region, and at the same time reducing the overall temperature baseline through a wide-area coverage mode.

[0037] The encapsulation of the instruction set follows the vehicle bus protocol specification. In each control cycle, a structured instruction packet containing parameters such as flow rate setting value, nozzle deflection angle, and pulse modulation duty cycle is output, and dynamic feasibility verification is performed through the hardware-in-the-loop simulation module to ensure that the instructions are executable within the working boundary of the hydraulic actuator. The entire process implements a closed-loop feedback mechanism, and the cooling effect is evaluated in real time through the thermodynamic parameter set of the next cycle, and the response coefficient of the control strategy is dynamically adjusted.

[0038] Step S3: According to the cooling instruction set, adjust the opening value of the piezoelectric valve in the annular pipeline, adjust the spraying area based on the flow velocity - surface area heat exchange efficiency, and at the same time recover the waste heat through the thermoelectric conversion device and optimize the system power consumption balance value, which specifically includes the following sub - steps: Step S31: According to the flow value in the cooling instruction set, adjust the opening of the corresponding piezoelectric valve in the annular pipeline proportionally to control the flow rate of the coolant; Convert the target flow value in the cooling instruction into a driving voltage signal, and precisely adjust the bending deformation of the valve disc in the annular pipeline through the inverse piezoelectric effect of the piezoelectric ceramic. The control algorithm adopts a feed - forward - feedback composite strategy: where, The control quantity obtained after calculation by the feed - forward - feedback composite control algorithm; Represents the function operation rule for calculating the feed - forward term value according to the parameters in the brackets; Represents the viscosity coefficient related to the temperature ; Represents the pipe length; Represents the pipe diameter; Represents the target mass flow rate; Represents the proportional coefficient. In feedback control, it determines the strength of the proportional relationship between the controller output and the error signal; Represents the error value; Represents the integral coefficient. In feedback control, the integral link performs an integral operation on the error, Determines the strength of the integral action, mainly used to eliminate the steady - state error of the system; Represents the integral of the error, which is an integral operation of the error over time.

[0039] Based on the coolant viscosity - temperature characteristic curve, preset the reference opening value, and at the same time form a closed - loop correction by real - time detecting the actual flow rate through an ultrasonic flowmeter. Aiming at the hysteretic non - linear characteristic of the piezoelectric valve, the system loads a pre - calibrated hysteresis compensation model and adopts a dynamic parameter estimation method based on the Preisach operator to eliminate the phase deviation in the opening adjustment process. When multi - valve collaborative control is implemented, a decoupling algorithm is used to suppress the hydraulic coupling interference between pipelines through a pressure equalizer, ensuring that the flow rate distribution of each spraying unit and the instruction value error are controlled within the allowable range.

[0040] Step S32: Based on the flow velocity - surface area heat exchange efficiency model, dynamically adjust the spraying area to cover the peak area of thermal stress; The heat exchange efficiency model establishes a dimensionless correlation formula with the Reynolds number and Nusselt number as the core through pre-set nozzle flow field simulation data, and calculates in real time the convective heat transfer coefficient of the coolant jet on the brake disc surface at different flow rates. The dynamic region matching algorithm rasterizes the thermal stress distribution map output by the finite element method, uses image morphology methods to identify the geometric center and contour features of the high stress region, and maps the physical position to the injection unit array of the annular pipeline through coordinate transformation. The control strategy combines the greedy algorithm with predictive control, preferentially activates the nozzle group covering the peak stress region, and at the same time adjusts the opening gradient of adjacent nozzles in advance based on the heat diffusion rate prediction model to form a moving cooling window that follows the migration of the hot spot. During the execution process, the actual cooling effect is fed back through an infrared thermal imager, and the surface emissivity parameter in the model is dynamically updated to compensate for the heat transfer attenuation caused by brake dust pollution.

[0041] Step S33: Convert the waste heat into electrical energy through a thermoelectric conversion device, calculate the system power consumption balance value in real time, and optimize the energy distribution strategy. The thermoelectric conversion device adopts a segmented bismuth telluride-based thermoelectric module layout. According to the radial temperature gradient distribution characteristics of the brake disc, large-size modules are configured in series in the high-temperature area to increase the output voltage, and small modules are connected in parallel in the low-temperature area to reduce the internal resistance loss. The power management unit implements the maximum power point tracking algorithm, dynamically adjusts the load impedance matching through the perturbation observation method, and at the same time uses a super capacitor bank to suppress the power fluctuations caused by sudden changes in braking conditions. The system power consumption balance model introduces a dynamic weight factor, and automatically assigns the usage priority of the recovered electrical energy according to the vehicle power supply status - when the battery has a low power, it preferentially supplies the cooling system pump, otherwise the electrical energy is fed back to the vehicle power grid. The optimization algorithm solves the optimal solution of energy distribution in real time through the Lyapunov optimization framework, minimizes the net power consumption of the system while ensuring the cooling performance, and at the same time predicts the waste heat generation trend in the next few seconds through a Kalman filter, and pre-adjusts the working point of the thermoelectric module to improve the overall energy conversion efficiency.

[0042] Step S4: The infrared thermal imager obtains the temperature field distribution, compares it with the initial parameter tolerance range, and corrects the injection angle and flow rate commands through the fuzzy PID algorithm to eliminate the temperature gradient error, which specifically includes the following sub-steps: Step S41: Obtain the brake disc temperature field distribution data through the infrared thermal imager and extract the actual temperature gradient curve. The high-resolution infrared thermal imager uses non-uniformity correction technology to eliminate optical system aberrations and improves the temperature measurement accuracy through real-time calibration with a blackbody radiation source. During data acquisition, the rotation phase signal of the brake disc is synchronously obtained, and a motion compensation algorithm is used to eliminate the image blur caused by the centrifugal effect. The preprocessed temperature field matrix identifies the boundaries of the heat dissipation grooves through an edge detection algorithm, divides the radial and circumferential detection regions, calculates the modulus of the temperature gradient in each direction using a second-order partial differential equation, and generates a gradient intensity cloud map distributed in polar coordinates. For local high-temperature regions, the system expands the detection range through morphological dilation operations, combines the historical temperature distribution characteristics to identify the migration trajectory of abnormal hot spots, and constructs a spatio-temporal feature vector containing the gradient amplitude, direction, and change rate.

[0043] Step S42: Compare the actual temperature gradient with the initial parameter tolerance range, calculate the deviation amount and input it into the fuzzy PID controller; The system performs pattern matching between the measured gradient feature vector and the tolerance interval preset in the initial parameter set, and quantifies the deviation degree through the membership function. The rule base of the fuzzy PID controller is designed with a three-dimensional input structure: the gradient deviation amount, the deviation change rate, and the historical cumulative error are jointly used as the inference inputs. For three types of working conditions, the system implements different adjustment strategies: (1) When the temperature gradient deviation amount is greater than the upper tolerance limit, the flow rate proportionality coefficient is preferentially increased; When it is detected that the global gradient deviation breaks through the upper tolerance limit, a strong intervention mode is triggered, and an exponential weighting mechanism is introduced in the proportionality coefficient adjustment link to dynamically increase the flow rate adjustment weight according to the area ratio of the over-limit region.

[0044] (2) When the deviation amount is within the tolerance range but there are local over-limit situations, the injection angle compensation value is preferentially adjusted; For the working conditions with local over-limit but overall controllable, the system starts the spatial decoupling algorithm. Through the mapping relationship between the coordinates of the thermal imager and the nozzle array, the optimal injection angle compensation vector is calculated, and the geometric relationship between the injection coverage area and the hot spot position is described using the Jacobian matrix to implement directional compensation cooling.

[0045] (3) When the deviation amount continuously falls below the lower tolerance limit, the flow rate proportionality coefficient is decreased and the waste heat recovery priority is enabled.

[0046] When the system enters the low-deviation steady state, the energy efficiency optimization mode is activated. While reducing the basic flow rate, the working voltage of the thermoelectric conversion device is increased, and more waste heat is converted into storable electrical energy through dynamic impedance matching.

[0047] Step S43: Adjust the coolant flow rate proportionality coefficient and the injection angle compensation value according to the deviation amount until the temperature gradient error is less than the preset threshold; The adjustment of the flow rate proportionality coefficient adopts a variable universe fuzzy strategy, dynamically stretching and shrinking the adjustment range of the universe according to the real-time working conditions, allowing a larger adjustment step size in the stage of severe fluctuations to quickly suppress the deviation, and switching to the fine adjustment mode when approaching the steady state. The calculation of the jet angle compensation value introduces a prediction and correction mechanism. Based on the jet coverage model preset by computational fluid dynamics, the flow field distribution at different deflection angles is predicted, and then the nozzle pointing angle is corrected by the actual cooling effect feedback of the infrared thermal imager. The jet angle prediction is expressed by the following formula: Wherein, represents the predicted jet angle at time ; represents the current time; represents the integral variable, which is used to cumulatively calculate the relevant quantities within the time period from to during the integral operation; represents the time lag; represents the kernel function, which reflects the influence weight of different times on the prediction result of the current time ; represents the complementary error function; represents the temperature value at time ; represents the reference temperature, which is used as the reference benchmark value of the temperature; represents the thermal diffusion coefficient.

[0048] During the closed-loop control process, the system continuously monitors the attenuation rate of the gradient error, and automatically injects a damping factor to smooth the adjustment process when detecting the overshoot risk. After the parameter tuning is completed, the updated control instruction is sent to the annular pipeline actuator through the timestamp synchronization mechanism, and at the same time, the adjusted system state is written into the operation log of the non-volatile memory to provide the initial parameter optimization benchmark for the subsequent control cycle. The entire adjustment process implements a multiple protection mechanism. When the single adjustment amplitude exceeds the safety threshold or the continuous adjustment fails to meet the convergence standard, the fault diagnosis process is started and the backup control strategy is switched.

[0049] Step S5: Based on the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process, use the neural network scoring model to generate an efficiency index, and dynamically adjust the algorithm weight and hardware response parameters, which specifically includes the following sub-steps: Step S51: Collect the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process in real time; The calculation of the temperature decay rate adopts the weighted second derivative of the temperature change amount within the sliding window: Wherein, Indicates at the moment The second derivative of temperature after being weighted by a sliding window; Indicates the half-width of the sliding window, that is, the range of the number of data points considered forward and backward from the current moment ; Indicates the index variable for summation, which is used to traverse the terms corresponding to different time offsets within the sliding window, and its value range is from to , representing different time intervals relative to the current moment ; Indicates the standard deviation, which is a parameter in the weight function and controls the decay rate of the exponential part in the weight function; sinc(*) is the sinc function, which is a sine basis function; Indicates a parameter used in the sinc function and affects the shape of the weight function; Indicates the temperature value at the moment ; Indicates the temperature value at the moment , that is, relative to the current moment , the time offset is the temperature at the moment; Indicates the temperature value at the moment , that is, relative to the current moment , the time offset is the temperature at the moment; Indicates the time interval, which represents the time difference between two adjacent sampling moments; The high-frequency noise interference to trend judgment is eliminated by the Savitzky-Golay filter; the pressure stability coefficient is calculated from the real-time data of the annular pipeline pressure sensor array, the main frequency characteristics of the pressure fluctuation are extracted by the principal component analysis method, and the degree of normal distribution deviation of the pressure distribution is determined by combining the Kolmogorov-Smirnov test; the energy recovery efficiency value is obtained by the four-quadrant measurement method, the hot-end temperature difference, output voltage and load current of the thermoelectric conversion device are synchronously collected, and the real-time energy conversion efficiency is calculated through the dynamic impedance matching model. All the original data streams are synchronized at the millisecond level through the timestamp alignment mechanism and filtered by the rationality verification module to filter out transient outliers.

[0050] Step S52: Input the above parameters into the pre-trained neural network scoring model to output the comprehensive efficiency index; The pre-trained three-channel convolutional neural network receives the normalized time-series feature matrix. The first channel processes the third-order difference spectrogram of the temperature decay rate, the second channel analyzes the time-frequency distribution characteristics of the pressure stability coefficient, and the third channel analyzes the load characteristic curve of the energy recovery efficiency. The feature fusion layer uses the attention mechanism to dynamically allocate the contribution weights of each parameter. The fully connected layer outputs a three-dimensional scoring vector including the thermal management efficiency, system stability, and energy efficiency. The neural network model is continuously optimized through transfer learning. When deployed on the embedded side, the knowledge distillation technology is used to compress the model size. At the same time, an online verification module is loaded to monitor the distribution shift of the input data. When it is detected that the working condition exceeds the coverage range of the training data, it automatically switches to the analogical reasoning mode based on fuzzy similarity to ensure the reliability of the evaluation results.

[0051] Step S53: Adjust the proportional-integral-derivative weights of the fuzzy PID algorithm according to the efficiency index, and optimize the response speed of the piezoelectric valve and the parameters of the thermoelectric conversion efficiency. The adjustment of the fuzzy PID weights adopts a reinforcement learning framework. The efficiency index is used as the input of the reward function, and the membership function shapes of the fuzzy rule tables of the proportional, integral, and differential terms are updated online through the policy gradient method. The optimization of the response speed of the piezoelectric valve is achieved by adjusting the rising slope of the driving signal. The system establishes a dynamic response model of the piezoelectric ceramic, calculates the optimal step response parameters according to the current valve body temperature and driving history, and uses the pre-distortion compensation algorithm to offset the influence of the hysteresis effect. The improvement of the thermoelectric conversion efficiency implements dual-mode control. In the high-efficiency index interval, the aggressive search mode of the maximum power point tracking algorithm is enabled, and in the low-efficiency interval, it switches to the conservative tracking strategy based on the historical optimal parameters. All parameter adjustments are verified by the Lyapunov stability criterion, and a parameter version tree is established in the non-volatile memory to support fast rollback and fault tracing in case of system performance anomalies.

[0052] Embodiment 2 As Figure 2 shown, Embodiment 2 of the present application provides an aircraft wheel cooling device, including: Data acquisition and preprocessing module 21: The temperature and pressure sensors collect the brake disc temperature and the wheel hub air pressure in real time. The embedded processing module calculates the temperature change rate and the pressure difference, and generates an initial thermodynamic parameter set, including the following sub-modules: Data preprocessing sub-module 211: Filter and denoise the temperature and air pressure data collected by the temperature and pressure sensors, and extract the effective signal segment; The temperature and pressure sensor collects the original signals of the brake disc temperature and the wheel hub air pressure in real time through an embedded hardware interface, and the signals are transmitted to the preprocessing module after analog-to-digital conversion. For the noise interference in the signals, first, the sliding window median filtering algorithm is used to eliminate the sudden impulse noise, and the window length is adaptively adjusted according to the dynamic characteristics of the signals to balance the real-time performance and the denoising effect; the temperature signal is further processed by low-pass filtering to filter out the high-frequency noise introduced by electromagnetic interference or mechanical vibration, and the air pressure signal suppresses the short-term fluctuations through weighted moving average filtering to ensure the data smoothness. Subsequently, based on the preset effective range of physical quantities, the filtered data is screened by thresholds to eliminate the outliers beyond the reasonable range; at the same time, combined with the dynamic change rate verification mechanism, if the change rate of the temperature or air pressure of consecutive sampling points exceeds the safety threshold, it is determined as a sensor failure or transient interference, and this section of data is marked as invalid and an alarm is triggered until the signal returns to the stable range.

[0053] Feature calculation sub-module 212: Calculate the real-time change rate of the brake disc temperature as the temperature change rate according to the effective signal segment, and calculate the difference between the wheel hub air pressure and the standard air pressure value as the pressure difference; After the effective signal segment is input, perform time series difference operation based on the temperature values of consecutive sampling points, and automatically switch to the linear fitting mode when it is detected that the temperature fluctuation exceeds the preset range: Among them, represents the optimal parameter value obtained by minimizing the objective function; represents finding the value that makes the subsequent expression take the minimum value; represents the total number of data points participating in the linear fitting; represents the discrete time index within the time window, and the value is , used to index the data points participating in the fitting; represents the time series data; represents the intercept term in the linear fitting model; represents the slope term in the linear fitting model; represents the time interval between adjacent time points.

[0054] The pressure difference calculation adopts a dynamic reference correction strategy. The standard air pressure value is adjusted in real time according to the ambient temperature through a pre-calibrated temperature-pressure compensation table. The actual pressure difference is the difference between the current air pressure and the dynamic reference, ensuring that the calculation result reflects the air pressure deviation under the actual working conditions; if the current air pressure data is marked as invalid due to preprocessing, the air pressure value of the previous effective cycle is used as a temporary reference to maintain the calculation continuity and trigger the sensor health diagnosis synchronously. All feature calculation results need to pass the logical verification. If data conflicts are detected, the output is frozen and the historical valid values are called, and at the same time, the system monitoring module is reported for fault tracing.

[0055] Data storage and management sub-module 213: Associates and stores the temperature change rate, pressure difference, and timestamp, generating a set of initial thermodynamic parameters including temperature, air pressure, temperature change rate, and pressure difference; The characteristic calculation results are bound to a high-precision timestamp and encapsulated into a structured data unit containing temperature, air pressure, temperature change rate, and pressure difference, and real-time storage management is achieved through a circular buffer. The buffer uses a circular queue structure to cover the historical data of a fixed time window. When writing, the oldest data is automatically eliminated to prevent memory overflow. At the same time, the data associated with abnormal events is marked with priorities to ensure that key information is retained until the persistent storage stage. The storage module periodically batches and compresses the buffer data and transfers it to a non-volatile memory. The compression algorithm needs to balance efficiency and resource occupancy and adapt to the computing power of the embedded system; a check code is appended before data writing to ensure integrity. If the check fails during reading, an attempt is made to recover from the buffer or backup area, and an error log is recorded for subsequent analysis. When the storage capacity is lower than the safety threshold, the sampling frequency is adaptively reduced or an alarm is triggered to prioritize ensuring the storage reliability of core data, forming a closed-loop logic from real-time processing to long-term storage management.

[0056] Thermal stress analysis and instruction generation module 22: Based on the set of initial parameters, calculates the thermal stress distribution value of the brake disc, and combines with the material threshold database to generate a cooling instruction set for coolant flow rate and injection angle, including the following sub-modules: Thermodynamic modeling and simulation sub-module 221: Based on the set of initial thermodynamic parameters, simulates the thermal stress distribution of the brake disc through the finite element analysis method, and outputs the stress peak and gradient values of each region; Based on the set of initial thermodynamic parameters, constructs a three-dimensional transient heat conduction model of the brake disc, and its control equation is: Among them, represents the density of the brake disc material; represents the specific heat; represents the temperature partial derivative with respect to time ; represents the Hamiltonian operator, which represents the divergence operation of a vector in three-dimensional space; represents the reference value of the thermal conductivity; represents the temperature coefficient; represents the temperature gradient; represents the braking friction coefficient, reflecting the friction characteristics between the friction surfaces during braking; represents the braking pressure; represents the contact speed, the relative movement speed between the brake disc and the friction component; represents the contact area.

[0057] First, non-uniform grids are divided according to the geometric structure and material properties of the brake disc. The grid density is encrypted in potential areas with high temperature gradients to improve calculation accuracy. The boundary conditions are dynamically set based on real-time temperature field and air pressure data. The heat flux density is derived from the temperature change rate and friction power consumption model, and the convective heat transfer coefficient is corrected by the aerodynamic formula related to the wheel hub air pressure and pressure difference. The solver uses an implicit iterative algorithm to calculate the transient heat conduction equation, and outputs the stress peak and gradient distribution maps of each area of the brake disc at fixed time intervals, while recording the position and change trend of the maximum stress point. To balance the computational efficiency and real-time requirements, the model runs in a simplified version on an embedded system, and optimizes matrix operations through a hardware acceleration unit to ensure that the simulation results are output to the downstream module within a millisecond-level delay.

[0058] Decision rule matching sub-module 222: According to the temperature resistance threshold and compressive strength threshold in the material threshold database, match the coolant demand level corresponding to the current stress distribution value; The thermo-stress distribution values output by the simulation are matched with the material threshold database in real time. The database pre-stores parameters such as the temperature resistance threshold, compressive strength threshold, and fatigue life curve of the brake disc material. The matching logic is divided into multiple levels of judgment: First, locate the risk area according to the stress peak. If the peak exceeds the instantaneous tensile strength threshold of the material, the highest-level cooling demand is directly triggered; if the peak is within the safe range but the gradient value exceeds the fatigue damage critical value, the cooling level is dynamically evaluated according to the gradient amplitude and duration. The cooling demand levels are divided into four levels: emergency, high, medium, and low, and each level corresponds to a preset coolant flow rate reference value and angle adjustment range. An environment adaptive strategy is introduced during the matching process; if there is no current material matching item in the database, an interpolation estimation mode based on similar material properties is enabled, and a manual verification request is triggered to ensure the reliability of the decision.

[0059] Control strategy generation sub-module 223: Combine the demand level and the heat exchange efficiency model to generate a cooling instruction set including coolant flow rate, injection angle, and duration; Generate an executable cooling instruction set based on the cooling demand level and the heat exchange efficiency model. The heat exchange efficiency model integrates coolant physical properties parameters, pipeline characteristics, and injection dynamics, and real-time solves the optimal flow rate and angle combination through a non-linear equation set. The control instruction set includes the flow rate setting value, injection angle, and duration, and limits the frequent switching of instructions through an anti-saturation algorithm to avoid actuator overload. After the instruction is generated, a redundant verification mechanism is embedded: if the current actuator state conflicts with the instruction, switch to the backup strategy, and at the same time report to the system monitoring module for fault tracing to ensure the safety and executability of the control instruction.

[0060] Cooling Execution and Waste Heat Recovery Module 23: According to the cooling instruction set, adjust the opening value of the piezoelectric valves in the annular pipeline, adjust the spraying area based on the flow velocity-surface area heat exchange efficiency, and at the same time recover waste heat through the thermoelectric conversion device and optimize the system power consumption balance value, including the following sub-modules: Actuator Control Sub-module 231: According to the flow rate value in the cooling instruction set, proportionally adjust the opening of the corresponding piezoelectric valve in the annular pipeline to control the coolant flow velocity; The flow rate setting value in the cooling instruction set is converted into a piezoelectric valve opening adjustment signal through the embedded control unit. The piezoelectric valve adopts a closed-loop control strategy. Based on the deviation between the real-time feedback of the flow sensor and the instruction value, a feedforward-feedback composite strategy is adopted through the control algorithm: Among them, The control quantity obtained after calculation by the feedforward-feedback composite control algorithm; Represents the function operation rule for calculating the feedforward term value according to the parameters in the brackets; Represents the temperature Related viscosity coefficient; Represents the pipe length; Represents the pipe diameter; Represents the target mass flow rate; Represents the proportional coefficient. In feedback control, it determines the strength of the proportional relationship between the controller output and the error signal; Represents the error value; Represents the integral coefficient. In feedback control, the integral link integrates the error, Determines the strength of the integral action, mainly used to eliminate the steady-state error of the system; Represents the integral of the error, which is the integral operation of the error over time.

[0061] Make the actual coolant flow velocity accurately match the target value. The relationship between the piezoelectric valve opening and the flow rate is compensated by a pre-calibrated non-linear curve to ensure high-resolution adjustment in the low-flow section and fast response in the high-flow section. The valve control logic integrates a fault protection mechanism: if the flow sensor fails or the valve feedback is abnormal, switch to the open-loop mode, perform conservative control according to the preset voltage-flow mapping table, and at the same time trigger an alarm and record the fault code. When multiple valves work together, a time-division multiplexing strategy is adopted to avoid instantaneous current overload, and hardware interlock is used to ensure the timing safety of the valve actions.

[0062] Dynamic Optimization Sub-module 232: Based on the flow velocity-surface area heat exchange efficiency model, dynamically adjust the spraying area to cover the peak area of thermal stress; Based on the flow velocity - surface area heat exchange efficiency model, the optimal injection area is calculated in real time. The model inputs include the peak position of thermal stress, the current flow velocity of the coolant, and the nozzle layout parameters, and the outputs are the activation priorities and angle adjustment amounts of each nozzle. The optimization objective is to maximize the cooling efficiency of the thermal stress peak area under limited flow rate. The greedy algorithm is used to dynamically allocate the flow rate: nozzles covering high - stress areas are activated first, and the injection angles are fine - tuned according to the direction of the stress gradient. The nozzle angles are adjusted by stepper motors or piezoelectric actuators, and the angle resolution needs to match the spatial accuracy of the thermal stress distribution. The optimization process embeds inertia constraints to limit the mechanical wear caused by frequent switching of nozzle angles, and predicts the adjustment trend of the injection area through learning from historical data to reduce the real - time calculation load.

[0063] Energy management sub - module 233: Convert waste heat into electrical energy through a thermoelectric conversion device, calculate the power consumption balance value of the system in real time, and optimize the energy distribution strategy; The thermoelectric conversion device is arranged close to the high - temperature area of the brake disc and uses the Seebeck effect to convert waste heat into electrical energy. The conversion efficiency is calibrated in real time through the temperature difference, and the output voltage is stabilized by a converter and then stored in a supercapacitor bank as the auxiliary power supply of the system. The power consumption balance value is calculated by monitoring the ratio of the main power supply to the recovered electrical energy in real time, and the power consumption strategy of the cooling system is dynamically adjusted: during high waste - heat recovery periods, the coolant flow rate is increased to the upper limit to accelerate heat dissipation, and the redundant electrical energy is fed back to the main power supply; during low waste - heat periods, it switches to the low - power consumption mode to limit the energy consumption of non - critical loads. The energy distribution strategy is implemented through a fuzzy logic controller, comprehensively considering the remaining battery capacity, braking intensity, and thermal stress risk level, and giving priority to ensuring the stability of core functions. If the total power consumption of the system exceeds the safety threshold, a degradation protocol is triggered, and secondary functions are gradually shut down and an energy warning is reported.

[0064] Temperature field feedback and dynamic correction module 24: The infrared thermal imager obtains the temperature field distribution, compares it with the initial parameter tolerance range, and corrects the injection angle and flow rate commands through the fuzzy PID algorithm to eliminate the temperature gradient error, including the following sub - modules: Temperature field monitoring and feature extraction sub - module 241: Obtain the brake disc temperature field distribution data through the infrared thermal imager and extract the actual temperature gradient curve; The infrared thermal imager captures the distribution data of the surface temperature field of the brake disc in real time through high-frame-rate scanning. After the original data is corrected for non-uniformity and compensated for bad pixels, it is converted into a temperature matrix. The feature extraction algorithm is based on matrix gradient calculation, generating two-dimensional temperature gradient curves along the radial and circumferential directions of the brake disc, and locating high-temperature aggregation areas and gradient mutation points. To improve the anti-interference ability, the data fusion module synchronously accesses the discrete temperature values of the temperature and pressure sensors, and calibrates the local measurement error of the infrared thermal imager through the weighted average method. The extracted temperature gradient curve is bound to the time stamp and output to the downstream control module, while an anomaly detection mechanism is embedded: if it is detected that the temperature field distribution seriously deviates from the physical heat transfer model, it is determined that the sensor fails or there is optical contamination, and the self-cleaning process is triggered or switched to the backup sensor.

[0065] Intelligent control strategy sub-module 242: Compare the actual temperature gradient with the tolerance range of the initial parameters, calculate the deviation amount and input it into the fuzzy PID controller; The actual temperature gradient curve is compared with the tolerance range preset by the set of initial thermodynamic parameters, and the deviation amount calculates the global error and the area of the local over-limit region through spatial integration. The fuzzy PID controller dynamically adjusts the control parameters according to the amplitude, duration and spatial distribution characteristics of the deviation amount: (1), When the temperature gradient deviation amount is greater than the upper tolerance limit, the flow rate proportionality coefficient is preferentially increased; When the overall deviation of the temperature gradient exceeds the upper tolerance limit, the fuzzy rule base preferentially increases the flow rate proportionality coefficient, quickly suppresses the diffusion of the temperature gradient by increasing the coolant flow rate, and at the same time limits the maximum flow rate increase to prevent thermal shock.

[0066] (2), When the deviation amount is within the tolerance range but there are local over-limit situations, the injection angle compensation value is preferentially adjusted; If the deviation is within the global tolerance range but there are local hot spots, based on the coordinates of the high-temperature area and the injection coverage model, the nozzle angle compensation value is dynamically calculated, and the injection angle of the corresponding area is preferentially adjusted to enhance directional cooling, while maintaining the flow stability of other areas.

[0067] (3), When the deviation amount continuously falls below the lower tolerance limit, the flow rate proportionality coefficient is reduced and the priority of waste heat recovery is enabled.

[0068] When the deviation continuously falls below the lower tolerance limit, the flow rate proportionality coefficient is gradually reduced to reduce energy consumption, and the priority of the waste heat recovery subsystem is increased, and the redundant coolant flow is allocated to the high-efficiency working range of the thermoelectric conversion device. The control strategy incorporates hysteresis anti-chattering logic to avoid command oscillations caused by instantaneous noise.

[0069] Dynamic parameter adjustment sub-module 243: Adjust the coolant flow rate proportionality coefficient and the injection angle compensation value according to the deviation amount until the temperature gradient error is less than the preset threshold; According to the control parameter correction amount output by the fuzzy PID, the coolant flow rate proportionality coefficient and the jet angle compensation value are adjusted in real time. The flow rate adjustment adopts a progressive approximation strategy: based on the current flow rate, it is superimposed step by step according to the proportionality coefficient increment, and after each step of adjustment, a fixed period is delayed to wait for the temperature field response until the global gradient error converges to the preset threshold. The jet angle compensation value is converted into the stepping motor drive pulse number of each nozzle through a coordinate mapping table. The angle correction amount is strictly matched with the spatial position of the high-temperature area. Based on the jet coverage model preset by computational fluid dynamics, the flow field distribution under different deflection angles is predicted, and then the nozzle pointing angle is corrected through the actual cooling effect feedback of the infrared thermal imager. The jet angle prediction is expressed by the following formula: Where, represents the predicted jet angle at time ; represents the current moment; represents the integration variable, which is used to cumulatively calculate the relevant quantities in the time period from to during the integration operation; represents the time lag; represents the kernel function, which reflects the influence weight of different moments on the prediction result of the current moment ; represents the complementary error function; represents the temperature value at time ; represents the reference temperature, which is used as the reference benchmark value of the temperature; represents the thermal diffusivity.

[0070] The adjustment process embeds a saturation protection mechanism: the single - time flow rate adjustment amplitude does not exceed the system safety limit value, and the angle compensation range is limited by the mechanical structure limit. If the error is still not eliminated after the maximum number of adjustments, it is determined as an uncontrollable working condition, triggering the system - level degradation operation and reporting the fault code, and synchronously recording the temperature field data for off - line analysis.

[0071] The efficiency evaluation and adaptive optimization module 25: Based on the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process, a neural network scoring model is used to generate an efficiency index, and the algorithm weights and hardware response parameters are dynamically adjusted, including the following sub - modules: The data monitoring and acquisition sub - module 251: Real - time acquisition of the temperature decay rate, pressure stability coefficient, and energy recovery efficiency value during the cooling process; The system real - time acquires the key parameters during the cooling process through multi - source sensors: the temperature decay rate is calculated from the time - series data of high - precision temperature sensors, and the calculation of the temperature decay rate uses the weighted second - derivative of the temperature change amount within a sliding window: Among them, represents the second derivative of temperature weighted by a sliding window at time ; represents the half-width of the sliding window, that is, the range of the number of data points considered forward and backward from the current time ; represents the index variable for summation, which is used to traverse the terms corresponding to different time offsets within the sliding window, and its value range is from to , representing different time intervals relative to the current time ; represents the standard deviation, which is a parameter in the weight function and controls the decay rate of the exponential part in the weight function; sinc(*) is the sinc function, which is a sine basis function; represents a parameter used in the sinc function and affects the shape of the weight function; represents the temperature value at time ; represents the temperature value at time , that is, relative to the current time , the time offset is ; represents the temperature value at time , that is, relative to the current time , the time offset is ; represents the time interval, which is the time difference between two adjacent sampling times; The pressure stability coefficient is quantified by the ratio of the standard deviation to the mean of the wheel hub air pressure, and the pressure fluctuation intensity is dynamically evaluated by combining the range analysis within the time window; the energy recovery efficiency value is calculated in real time by the ratio of the output electric energy of the thermoelectric conversion device to the input power of the brake disc waste heat, and the input power is calculated by the product of the temperature difference and the heat flow rate. After all parameters are aligned with time stamps by the hardware synchronization unit, they are encapsulated into a structured data set, and outlier filtering and data normalization processing are performed before transmission to ensure that the data input to the downstream module conforms to the model input specification.

[0072] Performance evaluation and model calculation sub-module 252: Input the above parameters into the pre-trained neural network scoring model and output the comprehensive performance index; The pre-trained lightweight neural network model is deployed in the embedded inference engine. The input layer receives the normalized temperature decay rate, pressure stability coefficient, and energy recovery efficiency values. The hidden layer extracts non-linear features through the activation function, and the output layer generates a comprehensive efficiency index in the range of 0-1, representing the global performance level of the current cooling system. During the model training stage, historical operating condition data and expert scoring labels are used for supervised learning, and an online incremental learning mechanism is introduced: regularly fine-tune the model weights with new data at the edge side to adapt to long-term drift problems such as material aging and environmental changes. The inference process embeds dynamic confidence detection. If the input parameter combination exceeds the training data distribution range, the expert rule base is triggered to assist in scoring to avoid optimization failure caused by model misjudgment.

[0073] Parameter optimization and adaptive adjustment sub-module 253: Adjust the proportional-integral-derivative weights of the fuzzy PID algorithm according to the efficiency index, and optimize the response speed of the piezoelectric valve and the parameters of the thermoelectric conversion efficiency; Based on the high or low efficiency index, the gradient descent method is used to online adjust the weight distribution of the proportional, integral, and differential terms of the fuzzy PID controller: when the efficiency is low, the proportional coefficient is preferentially increased to accelerate the system response; when the efficiency is high, the integral term is strengthened to suppress the steady-state error. At the same time, the weight overshoot is restricted by the constraint conditions to avoid control oscillation caused by parameter mutation. The optimization of the piezoelectric valve response speed is achieved by dynamically reconstructing the voltage-opening non-linear mapping relationship of the drive circuit, switching to the low-gain mode in the high-efficiency interval to reduce power consumption, enabling the high-gain mode at low efficiency to improve the flow regulation accuracy, and embedding the dead zone compensation algorithm to eliminate the influence of mechanical hysteresis. The optimization of the thermoelectric conversion efficiency is achieved by dynamically adjusting the load impedance matching point of the TEG array through pulse width modulation, making it always approach the maximum power transfer state, and coordinating the coolant flow control to maintain the temperature difference within the high-efficiency working range of the thermoelectric material. The parameter adjustment process adopts a rolling horizon optimization strategy, and the optimization effect is re-evaluated based on the latest efficiency index every fixed period. If the efficiency does not improve after continuous adjustments for multiple times, it will automatically roll back to the historical stable parameter set and trigger system-level fault diagnosis. All parameter adjustment instructions are issued through the safety interface, strictly limited within the hardware physical limit, and the key operation records are logged with timestamps to support offline traceability and policy iteration verification.

[0074] Corresponding to the above embodiments, an embodiment of the present invention provides a computer storage medium, including: at least one memory and at least one processor; The memory is used to store one or more program instructions; The processor is used to run one or more program instructions to execute an aircraft wheel cooling method; Corresponding to the above embodiments, an embodiment of the present invention provides a computer-readable storage medium, and the computer storage medium contains one or more program instructions, and the one or more program instructions are used to be executed by the processor to execute an aircraft wheel cooling method.

[0075] An embodiment disclosed by the present invention provides a computer-readable storage medium. Computer program instructions are stored in the computer-readable storage medium. When the computer program instructions run on a computer, the computer is caused to execute the above-mentioned aircraft wheel cooling method.

[0076] In an embodiment of the present invention, the processor may be an integrated circuit chip with signal processing capabilities. The processor may be a general-purpose processor, a digital signal processor (DSP for short), an application specific integrated circuit (ASIC for short), a field programmable gate array (FPGA for short), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0077] It can implement or execute the various methods, steps and logic block diagrams disclosed in the embodiments of the present invention. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor. The steps of the method disclosed in combination with the embodiments of the present invention may be directly embodied as being executed and completed by a hardware decoding processor, or executed and completed by a combination of hardware and software modules in the decoding processor. The software module may be located in a random access memory, flash memory, read-only memory, programmable read-only memory, or electrically erasable programmable memory, register, a mature storage medium in the art. The processor reads the information in the storage medium and combines its hardware to complete the steps of the above method.

[0078] The storage medium may be a memory, for example, it may be a volatile memory or a non-volatile memory, or may include both a volatile and a non-volatile memory.

[0079] Among them, the non-volatile memory may be a read-only memory (ROM for short), a programmable read-only memory (PROM for short), an erasable programmable read-only memory (EPROM for short), an electrically erasable programmable read-only memory (EEPROM for short), or a flash memory.

[0080] The volatile memory may be a Random Access Memory (RAM) which serves as an external cache. By way of example but not limitation, many forms of RAM are available, such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDR SDRAM), Enhanced SDRAM (ESDRAM), Synchlink DRAM (SLDRAM), and Direct Rambus RAM (DRRAM).

[0081] The storage media described in the embodiments of the present invention are intended to include, but not limited to, these and any other suitable types of memory.

[0082] Those skilled in the art should be aware that, in one or more of the above examples, the functions described in the present invention can be implemented by a combination of hardware and software. When applying software, the corresponding functions can be stored in a computer-readable medium or transmitted as one or more instructions or codes on a computer-readable medium. The computer-readable medium includes computer storage media and communication media, where the communication media includes any medium that facilitates the transfer of a computer program from one place to another. The storage media can be any available medium accessible by a general or special purpose computer.

[0083] The above specific embodiments have further elaborated on the purpose, technical solutions, and beneficial effects of the present invention. It should be understood that the above is only the specific embodiments of the present invention and is not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, and improvements made on the basis of the technical solutions of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for cooling an aircraft wheel, characterized in that: include: S1. The temperature and pressure sensor collects the brake disc temperature and wheel hub air pressure in real time. The embedded processing module calculates the temperature change rate and pressure difference to generate the initial thermodynamic parameter set. S2. Based on the initial parameter set, calculate the thermal stress distribution value of the brake disc, and generate a cooling instruction set of coolant flow and injection angle in combination with the material threshold database; S3. According to the cooling instruction set, the opening value of the piezoelectric valve of the annular pipeline is adjusted, and the injection area is adjusted according to the flow rate-surface area heat exchange efficiency. At the same time, the waste heat is recovered through the thermoelectric conversion device and the system power consumption balance value is optimized; S4, infrared thermal imager obtains temperature field distribution, compares initial parameter tolerance range, corrects injection angle and flow command through fuzzy PID algorithm, and eliminates temperature gradient error; S5. Based on the temperature decay rate, pressure stability coefficient and energy recovery efficiency value during the cooling process, a neural network scoring model is used to generate an efficiency index, and the algorithm weights and hardware response parameters are dynamically adjusted.

2. The method for cooling an aircraft wheel according to claim 1, characterized in that: The temperature and pressure sensors collect brake disc temperature and wheel hub air pressure in real time. The embedded processing module calculates the temperature change rate and pressure difference to generate the initial thermodynamic parameter set, including the following sub-steps: Filter and denoise the temperature and air pressure data collected by the temperature and pressure sensors to extract the effective signal segments; The real-time change rate of the brake disc temperature is calculated as the temperature change rate according to the effective signal segment, and the difference between the wheel hub air pressure and the standard air pressure value is calculated as the pressure difference; The temperature change rate, pressure difference and timestamp are associated and stored to generate an initial thermodynamic parameter set including temperature, air pressure, temperature change rate and pressure difference.

3. The method for cooling an aircraft wheel according to claim 1, characterized in that: Based on the initial parameter set, the thermal stress distribution value of the brake disc is solved, and the cooling instruction set of coolant flow and injection angle is generated in combination with the material threshold database, including the following sub-steps: Based on the initial thermodynamic parameter set, the thermal stress distribution of the brake disc is simulated by finite element analysis, and the stress peak and gradient values ​​of each area are output; According to the temperature resistance threshold and pressure resistance threshold in the material threshold database, the coolant requirement level corresponding to the current stress distribution value is matched; The demand level and heat exchange efficiency model are combined to generate a cooling instruction set including coolant flow rate, injection angle and duration.

4. The method for cooling an aircraft wheel according to claim 1, characterized in that: According to the cooling instruction set, the opening value of the piezoelectric valve of the annular pipeline is adjusted, and the injection area is adjusted according to the flow rate-surface area heat exchange efficiency. At the same time, the waste heat is recovered through the thermoelectric conversion device and the system power consumption balance value is optimized, including the following sub-steps: According to the flow value in the cooling instruction set, the opening of the corresponding piezoelectric valve in the annular pipeline is proportionally adjusted to control the flow rate of the coolant; Based on the velocity-surface area heat exchange efficiency model, the injection area is dynamically adjusted to cover the thermal stress peak area; The waste heat is converted into electrical energy through a thermoelectric conversion device, the system power consumption balance value is calculated in real time, and the energy allocation strategy is optimized.

5. The method for cooling an aircraft wheel according to claim 1, characterized in that: The infrared thermal imager obtains the temperature field distribution, compares the initial parameter tolerance range, and corrects the injection angle and flow command through the fuzzy PID algorithm to eliminate the temperature gradient error, including the following sub-steps: The temperature field distribution data of the brake disc is obtained through an infrared thermal imager, and the actual temperature gradient curve is extracted; Compare the actual temperature gradient with the initial parameter tolerance range, calculate the deviation and input it into the fuzzy PID controller; The coolant flow rate proportional coefficient and the injection angle compensation value are adjusted according to the deviation until the temperature gradient error is less than the preset threshold.

6. The method for cooling an aircraft wheel according to claim 5, characterized in that: The actual temperature gradient is compared with the initial parameter tolerance range, the deviation is calculated and input into the fuzzy PID controller, including the following sub-steps: When the temperature gradient deviation is greater than the upper limit of the tolerance, the flow rate proportional coefficient is increased first; When the deviation is within the tolerance range but there is a local over-limit, the injection angle compensation value is adjusted first; When the deviation value is continuously lower than the lower tolerance limit, the flow rate proportional coefficient is reduced and the waste heat recovery priority is enabled.

7. The method for cooling an aircraft wheel according to claim 1, characterized in that: Based on the temperature decay rate, pressure stability coefficient and energy recovery efficiency value during the cooling process, the neural network scoring model is used to generate the performance index, dynamically adjust the algorithm weight and hardware response parameters, including the following sub-steps: Real-time collection of temperature decay rate, pressure stability coefficient and energy recovery efficiency value during cooling process; The above parameters are input into the pre-trained neural network scoring model to output the comprehensive performance index; The proportional-integral-derivative weights of the fuzzy PID algorithm are adjusted according to the efficiency index, and the piezoelectric valve response speed and thermoelectric conversion efficiency parameters are optimized.

8. An aircraft wheel cooling device, characterized in that: include: Data acquisition and preprocessing module: The temperature and pressure sensors collect brake disc temperature and wheel hub air pressure in real time, and the embedded processing module calculates the temperature change rate and pressure difference to generate the initial thermodynamic parameter set; Thermal stress analysis and instruction generation module: Based on the initial parameter set, it solves the thermal stress distribution value of the brake disc and generates a cooling instruction set for coolant flow and injection angle in combination with the material threshold database; Cooling execution and waste heat recovery module: According to the cooling instruction set, adjust the opening value of the piezoelectric valve in the annular pipeline, adjust the injection area according to the flow rate-surface area heat exchange efficiency, and recover waste heat through the thermoelectric conversion device to optimize the system power consumption balance value; Temperature field feedback and dynamic correction module: The infrared thermal imager obtains the temperature field distribution, compares the initial parameter tolerance range, and corrects the injection angle and flow command through the fuzzy PID algorithm to eliminate the temperature gradient error; Performance evaluation and adaptive optimization module: Based on the temperature decay rate, pressure stability coefficient and energy recovery efficiency value during the cooling process, the neural network scoring model is used to generate the performance index and dynamically adjust the algorithm weights and hardware response parameters.

9. A computer storage medium, characterized in that: include: at least one memory and at least one processor; A memory for storing one or more program instructions; A processor, configured to run one or more program instructions to execute an aircraft wheel cooling method according to any one of claims 1 to 7.

Images