Hydraulic synchronous pushing construction control system based on multi-source information fusion

Abstract

The present application relates to high-end equipment manufacturing technology field, especially to a kind of hydraulic synchronous pushing construction control system based on multi-source information fusion, the system includes acquisition module, deviation module, dangerous module, instability module, hierarchical module, execution module and adjustment module.The present application introduces wind-induced buffeting spectrum analysis canyon wind field frequency domain characteristics, according to the difference between the main frequency of buffeting and natural frequency, dynamically adjusts the length of sliding window, suppresses the interference of wind-induced buffeting noise to attitude change rate, makes dangerous trend judgment accurately distinguish instantaneous buffeting and real deterioration, reduces false alarm rate, improves early warning lead time, and provides high-fidelity inclination rate input for instability determination, realizes the leap of canyon wind field pushing control from passive response to active accurate filtering, effectively solves the technical problems that wind-induced buffeting noise and real dangerous trend are difficult to distinguish due to fixed window calculation inclination rate in the prior art.

Description

Technical Field

[0001] This invention relates to the field of high-end equipment manufacturing technology, and in particular to a hydraulic synchronous jacking construction control system based on multi-source information fusion. Background Technology

[0002] The jacking construction of long-span bridges in strong wind environments such as canyons, river crossings, and valley crossings presents a key technical challenge in mountainous transportation infrastructure construction and bridge engineering in complex terrain. These bridge sites typically feature complex wind fields, strong gusts, significant turbulence, and variable wind directions. During the jacking process, the bridge beams are constantly exposed to dynamic wind loads, and wind-induced buffeting and overturning moments pose severe challenges to the beam's attitude stability, structural safety, and jacking synchronization accuracy. Simultaneously, the transient characteristics of wind speed and direction in canyon wind fields make attitude monitoring signals during jacking construction susceptible to high-frequency noise interference, and the risk of coupled resonance between the structure's natural frequencies and wind-induced excitation frequencies further exacerbates the uncertainty of instability. How to achieve multi-source information fusion perception, accurate judgment, and adaptive control during the jacking process under strong wind interference has become a critical technical bottleneck that urgently needs to be overcome in this field.

[0003] Chinese Patent Publication No. CN114578731A discloses a control system for a jacking device with monitoring function, comprising: a main computer for processing signals sent from various sensors and outputting control signals after processing; the main computer includes a CPU motherboard, a display module, an expansion module, and a storage module; the output terminal of the CPU motherboard is electrically connected to the input terminal of the display module; the CPU motherboard is bidirectionally electrically connected to the expansion module and the storage module; the expansion module is bidirectionally connected to a bus; multiple sub-controllers are connected to the main computer via a bus; each sub-controller includes a sensor acquisition module and a pump station control module, and both the sensor acquisition module and the pump station drive module are bidirectionally electrically connected to the bus; an actuator, the input terminal of which is electrically connected to the output terminal of the pump station control module; and a detection mechanism, the output terminal of which is electrically connected to the input terminal of the sensor acquisition module.

[0004] Therefore, the existing technology has the following problems: it relies on the independent threshold judgment of a single sensor signal, which makes it impossible to comprehensively identify the synchronization deviation of multi-parameter coupling, and it is easy to cause missed judgments or false judgments; it relies on a fixed sliding window or instantaneous value to calculate the tilt change rate, which makes it impossible to distinguish between wind-induced flutter noise and the actual attitude deterioration, and it is easy to cause frequent false alarms; it relies on static preset thresholds for anomaly judgment, which makes it unable to adapt to environmental changes during construction, and it is easy to miss the alarm during the strong disturbance period and false alarm during the weak disturbance period. Summary of the Invention

[0005] To address this, the present invention provides a hydraulic synchronous jacking construction control system based on multi-source information fusion, which overcomes the technical problem in the prior art where wind-induced buffeting noise is difficult to distinguish from the actual dangerous trend due to the fixed window calculation of the tilt angle change rate.

[0006] To achieve the above objectives, the present invention provides a hydraulic synchronous jacking construction control system based on multi-source information fusion, comprising: The data acquisition module is used to acquire in real time the maximum displacement difference between each jacking point, the pressure at each jacking point, the inclination angle, stress coefficient and natural frequency of the jacked beam, and the wind speed in the jacking construction area during the canyon wind field jacking construction. The deviation module is used to determine whether there is a deviation in the pushed beam based on the comparison result of the deviation index and the preset deviation threshold, wherein the deviation index is determined based on the displacement difference and the pressure. The danger module is used to determine whether the pushed beam is in danger when there is a deviation, based on the rate of change of the tilt angle and the stress coefficient within the sliding window. The length of the sliding window is adjusted according to the wind-induced flutter spectrum and the natural frequency. The wind-induced flutter spectrum is determined according to the wind speed and the natural frequency. The instability module is used to determine whether the pushed beam is about to become unstable when a danger occurs, based on the stability index and its changing trend. The stability index is determined based on the pressure, the tilt angle and the correction torque, and the correction torque is determined based on the wind speed. The grading module is used to determine the execution action of the jacking construction based on the stability index and the trend of change when instability is imminent. An execution module, which outputs execution commands to the hydraulic execution unit for jacking construction based on the execution action; The adjustment module is used to adjust the preset deviation threshold based on the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and the time correlation between the deviation judgment result and the danger judgment result.

[0007] Furthermore, the deviation module includes: The deviation index calculation submodule is used to calculate the deviation index based on the maximum displacement difference and the pressure at each jacking point; The threshold comparison submodule is used to compare the deviation index with a preset deviation threshold to determine whether a deviation exists.

[0008] Furthermore, hazardous modules include: The wind-induced buffeting spectrum calculation submodule is used to calculate the wind-induced buffeting spectrum based on the time-series data of the wind speed and extract the dominant frequency of the buffeting. A sliding window dynamic adjustment submodule is used to adjust the length of the sliding window according to the difference between the dithering main frequency and the natural frequency; The tilt angle change rate calculation submodule is used to perform differential calculation on the tilt angle within the sliding window to obtain the tilt angle change rate. The danger trend determination submodule is used to compare the tilt angle change rate and the stress coefficient with their respective preset tilt angle change rate threshold and preset stress coefficient threshold. When the tilt angle change rate is greater than the preset tilt angle change rate threshold and the stress coefficient is greater than the preset stress coefficient threshold, it is determined that the pushed beam is in danger.

[0009] Furthermore, the sliding window dynamic adjustment submodule includes: A frequency comparison unit is used to obtain a first frequency difference based on the difference between the dithering main frequency and the natural frequency. The window length calculation unit is used to extend the length of the sliding window according to the first frequency difference and the preset resonance threshold when the first frequency difference is less than the preset resonance threshold; and to restore the length of the sliding window to the preset window length when the first frequency difference is greater than or equal to the preset resonance threshold.

[0010] Furthermore, the instability module includes: The static load moment calculation submodule is used to calculate the total jacking force on the outer side and the total jacking force on the inner side according to the pressure, and to calculate the stabilizing moment and the static load overturning moment in combination with the tilt angle. The corrected moment calculation submodule is used to calculate the wind-induced overturning moment based on the wind speed; The stability index calculation submodule is used to calculate the stability index based on the stabilizing moment, the static overturning moment, and the wind-induced overturning moment. The trend calculation submodule is used to determine the trend based on the stability index sequence within a preset stable period; The instability trend determination submodule is used to determine whether the pushed beam is about to become unstable based on the stability index and the change trend.

[0011] Furthermore, the hierarchical module includes: The risk level determination submodule is used to determine the current risk level based on the stability index and the trend of change. The risk level includes low risk, medium risk and high risk. The lateral correction execution submodule is used to select lateral active correction when the risk level is medium risk.

[0012] Furthermore, the lateral correction execution submodule includes: Lateral offset calculation unit, which is used to calculate the current lateral offset based on the tilt angle and beam geometry parameters; The correction force determination unit is used to determine the required lateral correction force based on the lateral offset and the preset correction gain coefficient. The lateral push execution unit is used to control the auxiliary jacks arranged on the side of the beam to apply the lateral correction force until the updated lateral offset is restored to the preset allowable range.

[0013] Furthermore, the execution module includes: The instruction generation submodule is used to generate corresponding hydraulic control instructions based on the lateral active correction action output by the lateral correction execution submodule. The instruction output submodule is used to output the hydraulic control instructions to the hydraulic actuator of the jacking operation.

[0014] Furthermore, the adjustment modules include: The frequency deviation calculation submodule is used to calculate the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and obtain the second frequency difference. The time correlation calculation submodule is used to determine the time correlation based on the number of deviation events and the number of dangerous events. The number of deviation events is the total number of times the deviation judgment result indicates that there is a deviation within the preset observation period, and the number of dangerous events is the number of times the danger judgment result indicates that there is a danger within the preset correlation window after the deviation judgment result indicates that there is a deviation. The threshold adjustment submodule is used to adjust the preset deviation threshold based on whether the second frequency difference is less than the preset resonance threshold and whether the time correlation is greater than the preset correlation threshold.

[0015] Furthermore, the threshold adjustment submodule includes: A resonance risk adjustment unit is used to determine a first adjustment range based on the second frequency difference and the preset resonance threshold when the second frequency difference is less than the preset resonance threshold. A statistical correlation adjustment unit is used to determine a second adjustment range based on the difference between the time correlation degree and the preset correlation threshold when the time correlation degree is greater than the preset correlation threshold. The integrated adjustment unit is used to adjust the preset deviation threshold according to the first adjustment range and the second adjustment range.

[0016] Compared with existing technologies, this method achieves joint evaluation of position synchronization and load balance on a unified scale by normalizing the maximum displacement difference and the pressure at each jacking point and then fusing them into a single dimensionless deviation index using Euclidean distance. This overcomes the shortcomings of existing technologies, which suffer from fragmented multi-parameter coupling risks due to the different dimensions of displacement and pressure and separate threshold judgments. Furthermore, it introduces wind-induced buffeting spectrum analysis to analyze the frequency domain characteristics of the canyon wind field in real time. The sliding window length of the tilt angle change rate is dynamically adjusted based on the difference between the dominant buffeting frequency and the natural frequency of the beam, suppressing the interference of wind-induced buffeting noise on the attitude change rate extraction from the frequency domain root. This enables the danger trend judgment to distinguish between wind-induced instantaneous buffeting and actual attitude deterioration, solving the problem of difficulty in distinguishing noise and danger caused by the existing fixed-window calculation of the tilt angle change rate. Finally, it employs logic and judgment based on the tilt angle change rate and stress coefficient, and integrates the stability of static load and wind load. By employing a fixed-index dynamic evaluation and a dual-driven threshold adaptive correction mechanism based on the resonance frequency difference and the statistical correlation between deviation and hazard events, a physically progressive causal chain from deviation perception and hazard identification to instability early warning is constructed. This overcomes the limitations of single-parameter independent threshold judgment and static threshold setting, achieving adaptive matching of control parameters with the construction environment and structural state. It also solves the problems of missed alarms during periods of strong disturbance and false alarms during periods of weak disturbance by using fixed thresholds. By integrating frequency domain resonance analysis, time domain adaptive filtering, and static-dynamic coupling models, a technological leap from passive response to proactive and precise prevention and control is achieved. This significantly improves the safety, adaptability, and control accuracy of hydraulic synchronous jacking construction in complex environments such as canyon wind fields, and effectively solves the technical problem in existing technologies where wind-induced flutter noise is difficult to distinguish from the actual hazard trend due to the fixed window calculation of the tilt angle change rate. Attached Figure Description

[0017] Figure 1 This is a schematic diagram of the hydraulic synchronous jacking construction control system based on multi-source information fusion in this embodiment; Figure 2 This is a logic diagram for determining whether the pushed beam is in danger, as shown in this embodiment. Figure 3 This is a logic diagram for determining whether the pushed beam is about to become unstable in this embodiment. Figure 4 This is a flowchart for adjusting the preset deviation threshold in this embodiment. Detailed Implementation

[0018] To make the objectives and advantages of the present invention clearer, the present invention will be further described below with reference to embodiments; it should be understood that the specific embodiments described herein are merely for explaining the present invention and are not intended to limit the present invention.

[0019] Preferred embodiments of the present invention will now be described with reference to the accompanying drawings. Those skilled in the art should understand that these embodiments are merely illustrative of the technical principles of the present invention and are not intended to limit the scope of protection of the present invention.

[0020] Please see Figure 1 As shown, it is a schematic diagram of the hydraulic synchronous jacking construction control system based on multi-source information fusion in this embodiment.

[0021] This embodiment provides a hydraulic synchronous jacking construction control system based on multi-source information fusion, including: The data acquisition module is used to acquire in real time the maximum displacement difference between each jacking point, the pressure at each jacking point, the inclination angle, stress coefficient and natural frequency of the jacked beam, and the wind speed in the jacking construction area during the canyon wind field jacking construction. The deviation module is used to determine whether there is a deviation in the pushed beam based on the comparison result of the deviation index and the preset deviation threshold, wherein the deviation index is determined based on the displacement difference and the pressure. The danger module is used to determine whether the pushed beam is in danger when there is a deviation, based on the rate of change of the tilt angle and the stress coefficient within the sliding window. The length of the sliding window is adjusted according to the wind-induced flutter spectrum and the natural frequency. The wind-induced flutter spectrum is determined according to the wind speed and the natural frequency. The instability module is used to determine whether the pushed beam is about to become unstable when a danger occurs, based on the stability index and its changing trend. The stability index is determined based on the pressure, the tilt angle and the correction torque, and the correction torque is determined based on the wind speed. The grading module is used to determine the execution action of the jacking construction based on the stability index and the trend of change when instability is imminent. An execution module, which outputs execution commands to the hydraulic execution unit for jacking construction based on the execution action; The adjustment module is used to adjust the preset deviation threshold based on the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and the time correlation between the deviation judgment result and the danger judgment result.

[0022] In this embodiment, the following parameters are predetermined and stored in the control system based on bridge design drawings, structural calculation sheets, jack factory parameters, material properties, and engineering specifications: Geometric support parameters include the outer lever arm, inner lever arm, beam width, and center of gravity height. The outer lever arm is the horizontal distance from the outer jack thrust line to the overturning rotation center, i.e., the outermost support point of the beam. The inner lever arm is the horizontal distance from the inner jack thrust line to the same rotation center. The center of gravity height is obtained through finite element calculations or design values ​​based on the mass distribution of each beam segment. Beam physical property parameters include self-weight, material thermal expansion coefficient, and material density of each segment. Jack characteristic parameters include the effective piston area and geometric parameters such as the position of the overturning rotation center. These parameters are calibrated and entered before the jacking operation and participate in real-time control calculations as fixed constants. Those skilled in the art can obtain these parameters based on conventional engineering practices.

[0023] In this embodiment, the preset jacking step length refers to the travel distance of a single horizontal jack in a single jacking cycle. It is determined based on the maximum stroke of the jack and the allowable range of beam deformation under stress. In this embodiment, the preset step length is stored in the control system according to the selected jack model. It is used to calculate the theoretical target stroke of each jacking point in real time and compare it with the actual displacement to obtain the maximum displacement difference. This provides an accurate displacement reference for synchronous deviation control, supports unequal step length differential speed correction in curved jacking, reduces the complexity of real-time calculation, and prevents the beam from being over-pushed or under-pushed, thus ensuring construction accuracy and progress.

[0024] In this embodiment, the current radius of curvature is obtained by pre-setting the design alignment and querying the cumulative travel distance: Before the jacking construction, the horizontal curve parameters of the bridge are entered into the control system. These parameters include the starting mileage, radius of curvature, and transition curve parameters for each alignment segment (where the transition curve parameter is the spiral curve parameter A, which can be directly read from the bridge design drawings or based on the circular curve radius R and the total length L of the transition curve). Calculation; During the jacking process, the current linear segment is determined in real time based on the total jacking stroke accumulated by the displacement sensor, and the design radius of curvature is calculated according to the corresponding formula based on the linear segment: for circular curve segments, a constant radius is used; for transition curve segments, the spiral curve formula R=A is used. 2 The calculation is / L, where L is the current arc length calculated from the starting point of the transition curve.

[0025] A wire-type displacement sensor is installed at the end of the piston rod of each jack to collect the piston stroke of each jack in real time. This stroke is the actual stroke of each jacking point. The theoretical target stroke of each point is calculated in real time based on the preset jacking step length, the current radius of curvature, and the beam width. The maximum absolute value of the difference between the actual stroke and the corresponding theoretical target stroke is taken as the maximum displacement difference. A resistance strain gauge pressure sensor is installed at the oil inlet of each jack's hydraulic cylinder. The current signal is converted into the pressure value of each jacking point. This pressure value is the final pressure obtained for each jacking point. Then, the pressure is converted into the real-time jacking stroke of each jacking point based on the effective area of ​​each jack's piston. Force; MEMS biaxial tilt sensors are installed at key sections at the bottom of the beam, directly outputting the tilt angle via a digital interface; fiber optic strain gauges are attached to key stress sections of the beam, and the wavelength drift is read by a fiber optic demodulator and converted into measured axial strain. Simultaneously, a temperature sensor measures the surface temperature to calculate the theoretical thermal expansion strain, and the ratio of the measured strain to the theoretical strain is used as the stress coefficient; a MEMS accelerometer is deployed at the mid-span of the beam to continuously collect vibration time history data, and the peak frequency of the power spectrum is extracted as the natural frequency using fast Fourier transform; an ultrasonic anemometer is installed on the top of the temporary pier, and the wind speed is calculated in real time based on the ultrasonic time-of-flight method. All sensor data are transmitted to the central controller via fieldbus, and each parameter, except for the natural frequency, is continuously acquired at an update rate of no less than 10Hz.

[0026] It should be noted that under conditions of insufficient environmental excitation, such as no wind or light wind, the beam's vibration response is weak, which may lead to unstable real-time natural frequency extraction results. To ensure the continuity of the system's judgment, this embodiment adopts the following strategy: when the peak amplitude of the power spectrum of the acceleration time history data is lower than a preset energy threshold, it is determined that the environmental excitation is insufficient. In this case, the real-time extracted frequency value is no longer used, but the preset theoretical natural frequency is directly used as the current natural frequency. This theoretical natural frequency is obtained through finite element calculation based on the beam design parameters. This strategy ensures that stable natural frequency input can be obtained for judgments such as dynamic adjustment of the sliding window and adaptive deviation threshold under various wind conditions.

[0027] The preset deviation threshold is a critical value used to determine whether the synchronization deviation index exceeds the limit; its physical meaning is the maximum permissible degree of synchronization deviation. The value is determined by comprehensively considering the accuracy requirements of the jacking construction, the structural safety margin, and the sensor noise level. During the system calibration phase, finite element simulation is used with the stability index decreasing to a preset simulation safety reference value as the boundary condition. The synchronization deviation is gradually increased until the stability index reaches this reference value, and the deviation index at this point is recorded as the initial reference for the preset deviation threshold. After correction through on-site calibration tests, in this embodiment, the preset deviation threshold is initially set to 0.8. This calibration method ensures that the preset deviation threshold matches the actual safety margin of the structure, preventing the deviation judgment from being too sensitive, leading to frequent warnings, or too sluggish, causing risk accumulation, thereby fundamentally ensuring the safety and continuity of the jacking construction.

[0028] The preset simulation safety reference value refers to the critical reference value used in finite element simulation to calibrate the preset deviation threshold. Its physical meaning is the minimum allowable overturning safety margin. The value is determined by comprehensively referencing the requirement in the "Crane Design Code" that the overturning stability coefficient should not be less than 1.2, as well as the dynamic characteristics of jacking construction. In this embodiment, it is set to 1.2. Its beneficial effect is that it provides a clear engineering boundary for calibrating the preset deviation threshold, ensuring that the deviation judgment matches the actual safety margin of the structure, balancing early warning sensitivity and safety.

[0029] The preset observation period refers to a fixed time window used to adjust the difference between the peak frequency and natural frequency of the wind-induced buffeting spectrum, as well as the time correlation between the deviation judgment result and the hazard judgment result. The value is determined by comprehensively considering the time-varying characteristics of the dominant frequency of the canyon wind field buffeting and the progress rhythm of the jacking construction. In this embodiment, based on the measured data of the canyon wind field and the system response requirements, the preset observation period is set to 30 seconds, balancing statistical reliability and adaptive response speed. This allows the system to adjust its judgment sensitivity promptly according to changes in wind field conditions or construction disturbance intensity, while avoiding frequent threshold changes due to short-term random fluctuations. Those skilled in the art can routinely adjust the period within the range of 15 to 60 seconds based on the actual wind field change rate and construction progress.

[0030] The acquisition module acquires multi-source heterogeneous parameters such as maximum displacement difference, pressure, tilt angle, stress coefficient, natural frequency, and wind speed in real time during the canyon wind-driven jacking construction. The deviation module integrates the displacement difference and pressure into a deviation index and compares it with a preset threshold, achieving a joint assessment of position synchronization and load balance. The hazard module adaptively adjusts the sliding window length based on the wind-induced buffeting spectrum to suppress high-frequency noise and uses the logic of tilt angle change rate and stress coefficient to determine the hazard trend, effectively distinguishing between wind-induced buffeting interference and actual attitude deterioration. The instability module calculates the stability index based on pressure, tilt angle, and wind-induced correction torque and monitors its changing trend, incorporating static load and wind load into the dynamic assessment of overturning safety margin, identifying instability risks in advance from three dimensions: current value, direction of change, and rate of deterioration. The grading module takes differentiated actions based on the stability index and its trend, and the adjustment module further adaptively corrects the deviation threshold based on the difference between the buffeting main frequency and the natural frequency and the time correlation between deviation and hazard, forming a closed-loop optimization. The entire system follows a physical progressive causal chain from synchronization deviation and dangerous trends to instability warning. It deeply couples time-domain filtering, frequency-domain resonance analysis, static and wind load superposition, statistical learning and physical mechanism driving, achieving a leap from passive response to active and precise prevention and control. It significantly improves the safety, adaptability and control accuracy of jacking construction in complex environments such as canyon wind fields, and effectively solves the technical problem in the existing technology that it is difficult to distinguish between wind-induced flutter noise and real dangerous trends due to the fixed window calculation of the tilt angle change rate.

[0031] Specifically, the deviation module includes: The deviation index calculation submodule is used to calculate the deviation index based on the maximum displacement difference and the pressure at each jacking point; The threshold comparison submodule is used to compare the deviation index with a preset deviation threshold to determine whether a deviation exists.

[0032] The preset maximum allowable displacement difference refers to the maximum permissible value of the displacement difference between each jacking point and the theoretical target during the jacking construction. Its physical meaning is the control limit of the synchronous deviation. The value is determined by comprehensively considering the accuracy requirements of the jacking construction, the safety margin of the beam structure, and the noise level of the displacement sensor. In this embodiment, based on engineering experience and finite element simulation results, it is initially set to 10mm. This parameter is used to normalize the real-time maximum displacement difference, converting it into a dimensionless displacement deviation component, thereby integrating it with the pressure deviation on the same scale.

[0033] The theoretical target pressure at each jacking point is calculated as follows: Let the beam's self-weight be G, the horizontal distance from the center of gravity to the inner jack be a, and the horizontal distance to the outer jack be b. The distance between the inner and outer jacks is D = a + b. Using the inner jack as the fulcrum, from the moment balance, we get Fout × D = G × a, thus obtaining the theoretical support reaction force of the outer jack: Fout = G × a / D. Using the outer jack as the fulcrum, we obtain the theoretical support reaction force of the inner jack: Fin = G × b / D. Dividing the support reaction force by the effective area of ​​the corresponding jack piston yields the theoretical target pressure. The lever arms a and b are fixed design values, obtained from the bridge design drawings and pre-stored in the control system, and do not change with the radius of curvature. This calculation only considers self-weight and does not involve wind load or lateral movement. The obtained theoretical target pressure serves as the benchmark for pressure normalization in the deviation module.

[0034] The preset maximum allowable pressure deviation refers to the absolute value of the maximum permissible difference between the real-time pressure at each jacking point and the corresponding theoretical target pressure. The value is determined by comprehensively considering the rated working pressure of the jack, the control accuracy of the hydraulic system, the safety margin of the beam's stress, and the measurement error of the pressure sensor. In this embodiment, through finite element simulation and on-site trial jacking calibration verification, the preset maximum allowable pressure deviation is initially set to 10% of the rated pressure. This provides a standardized and normalized benchmark for the pressure deviation, allowing the pressure deviation component and the displacement deviation component to merge on the same dimensionless scale, avoiding the dominance of a single parameter's deviation index due to dimensional differences. In the deviation index calculation submodule, the real-time maximum displacement difference is divided by the preset maximum allowable displacement difference to obtain the normalized displacement deviation. For each jacking point, the absolute value of the difference between its real-time pressure and the corresponding theoretical target pressure is calculated, and the maximum value among all points is taken as the maximum pressure deviation. The maximum pressure deviation is divided by the preset maximum allowable pressure deviation to obtain the normalized pressure deviation. Based on the normalized displacement deviation and the normalized pressure deviation, the deviation index is calculated using the Euclidean distance formula. When the deviation index exceeds the preset deviation threshold, a synchronization deviation is determined to exist; otherwise, no synchronization deviation is determined, the current jacking state is maintained, and subsequent dangerous trend judgments are not triggered.

[0035] By employing Euclidean distance to fuse normalized displacement deviation and normalized pressure deviation into a deviation index, joint evaluation of position synchronization and load balance under a unified scale is achieved. This overcomes the problem of fragmented multi-parameter coupling risks caused by the different dimensions of displacement and pressure and separate threshold judgments in traditional methods. Euclidean distance fusion can simultaneously reflect two abnormal modes: single-parameter over-limit and dual-parameter synchronous degradation. This ensures that the deviation index is equally sensitive to hidden conditions such as displacement not exceeding limits but pressure imbalance or pressure not exceeding limits but displacement deviation. By using preset deviation thresholds for binarization judgment, the next level of danger trend judgment is triggered only when the index exceeds the threshold, forming an early warning and progressively advancing control chain. This significantly reduces computational load and false alarm rate while ensuring control sensitivity, providing a reliable and self-consistent input starting point for subsequent danger trend identification and instability judgment.

[0036] Please see Figure 2 As shown, it is the logic diagram for determining whether the pushed beam is in danger in this embodiment.

[0037] Specifically, the hazardous modules include: The wind-induced buffeting spectrum calculation submodule is used to calculate the wind-induced buffeting spectrum based on the time-series data of the wind speed and extract the dominant frequency of the buffeting. A sliding window dynamic adjustment submodule is used to adjust the length of the sliding window according to the difference between the dithering main frequency and the natural frequency; The tilt angle change rate calculation submodule is used to perform differential calculation on the tilt angle within the sliding window to obtain the tilt angle change rate. The danger trend determination submodule is used to compare the tilt angle change rate and the stress coefficient with their respective preset tilt angle change rate threshold and preset stress coefficient threshold. When the tilt angle change rate is greater than the preset tilt angle change rate threshold and the stress coefficient is greater than the preset stress coefficient threshold, it is determined that the pushed beam is in danger.

[0038] The wind-induced buffeting spectrum calculation submodule performs a Fast Fourier Transform (FFT) on continuously collected wind speed time-series data within a preset wind-induced time window to obtain the power spectral density curve of the wind-induced buffeting, and extracts the frequency component with the largest amplitude in the curve as the dominant frequency of the buffeting. The preset wind-induced time window refers to the duration of the data segment extracted when performing a FFT on the continuously collected wind speed time-series data; its physical meaning is the time coverage range of a single spectral analysis. It depends on the frequency resolution requirements and the real-time response requirements. In this embodiment, considering the typical period of wind speed change rate in the canyon wind field (10-20 seconds) and the natural frequency of the beam (0.2-2Hz), the preset wind-induced time window is set to 30 seconds, with a frequency resolution of approximately 0.033Hz. This accurately captures resonance risk, balancing spectral resolution and real-time response, ensuring that the extraction of the dominant frequency of the buffeting meets the accuracy requirements for resonance detection while also being able to update quickly during sudden changes in wind conditions. Those skilled in the art can adjust the time window within the range of 20 to 60 seconds according to the actual wind field characteristics and structural dynamic characteristics. In this embodiment, a segmented processing method with a 50% overlap rate is adopted, and the dithering main frequency is updated every 15 seconds to achieve quasi-continuous output.

[0039] The tilt angle change rate calculation submodule performs numerical difference processing on the continuously acquired tilt angle time-series data within the current sliding window. Specifically, it divides the tilt angle difference between adjacent sampling points within the window by the sampling time interval to obtain the instantaneous angular velocity sequence. Then, it uses the least squares method to fit the linear trend of this sequence, and the resulting slope is the tilt angle change rate within the window. In this embodiment, the tilt sensor sampling frequency is 20Hz, and the initial length of the sliding window is 5 seconds. The calculation method remains unchanged even when the window is dynamically extended. Using linear regression instead of simple average difference can effectively suppress the influence of high-frequency noise on the change rate estimation, making the extraction of attitude deterioration trends more robust.

[0040] The preset tilt angle change rate threshold is a critical value used to determine whether the tilt angle change rate is abnormal. The value is determined by comprehensively considering the tilt angle fluctuation amplitude caused by road surface irregularities, minor settlement of temporary piers, etc., and the structural safety margin during normal jacking. In this embodiment, based on finite element simulation and on-site trial jacking data, the preset tilt angle change rate threshold is set to 0.05° / s. This provides a quantitative limit for the attitude deterioration index in the judgment of dangerous trends, avoiding frequent false alarms caused by normal fluctuations, while ensuring that a timely warning is triggered when the attitude change rate truly exceeds the safe range. Those skilled in the art can routinely adjust the threshold within the range of 0.02° / s to 0.1° / s according to the actual beam stiffness, jacking speed, and monitoring noise level.

[0041] The preset stress coefficient threshold is a critical value used to determine whether the stress coefficient is abnormal. The value is determined by comprehensively considering the elastic limit of the steel and the safety margin for temperature effects. In this embodiment, based on the material properties and structural analysis of Q345 steel, the preset stress coefficient threshold is set to 1.2. This provides a quantitative benchmark for stress anomaly indicators in the determination of dangerous trends, effectively eliminating interference from normal thermal expansion and strain caused by temperature changes. It ensures that the warning is triggered only when the structure experiences additional constraint stress or stress imbalance, conforming to the underlying physical logic that stress anomalies must be separated from thermal effects. Those skilled in the art can routinely adjust the threshold within the range of 1.1 to 1.5 according to the actual steel grade and structural constraint conditions.

[0042] It should be noted that when the rate of change of inclination angle exceeds a preset threshold and the stress coefficient exceeds a preset threshold, the hazard trend determination submodule determines that the pushed beam is in danger and triggers the instability module to determine instability. When only the rate of change of inclination angle exceeds a preset threshold and the stress coefficient is less than or equal to a preset threshold, or when only the stress coefficient exceeds a preset threshold and the rate of change of inclination angle is less than or equal to a preset threshold, the hazard trend determination submodule determines the situation to be in a state of alert and does not trigger the instability module.

[0043] In the alert state, the system intensifies monitoring, including but not limited to dynamically shortening the preset observation period of the adjustment module based on the proximity of the current wind-induced buffeting dominant frequency to the natural frequency. The shortened period should be no less than the minimum data window length required for calculating the tilt change rate in the danger module, in order to more quickly reassess whether to escalate to a danger state. If the alert state continues for more than the preset alert duration without escalation, the system automatically restores normal monitoring parameters; if both parameters exceed the limit simultaneously during the alert state, it is immediately escalated to danger and instability judgment is triggered.

[0044] In this embodiment, the preset attention duration is set to 10 seconds. This duration is determined based on the typical evolution time of 10 to 30 seconds from a single-parameter anomaly to a two-parameter anomaly and the judgment cycle of the dangerous module. It can effectively capture the trend of danger escalation and avoid the system being in an enhanced monitoring state for a long time.

[0045] The wind-induced buffeting spectrum calculation submodule extracts the dominant buffeting frequency from wind speed time-series data and compares it with the natural frequency of the beam. The sliding window dynamic adjustment submodule adaptively adjusts the window length for calculating the tilt angle change rate based on the frequency difference, thus solving the problem of time-varying signal-to-noise ratio of tilt angle signals in canyon wind fields from the root of the frequency domain. At the same time, the tilt angle change rate and stress coefficient are determined by logical integration, ensuring that the danger warning is triggered only when the attitude continues to deteriorate and the structural stress is abnormal. This effectively isolates false alarms caused by wind-induced buffeting spikes or temperature transients of a single parameter, achieving a unified approach to accurate extraction of danger trends and suppression of false alarms.

[0046] Specifically, the sliding window dynamic adjustment submodule includes: A frequency comparison unit is used to obtain a first frequency difference based on the difference between the dithering main frequency and the natural frequency. The window length calculation unit is used to extend the length of the sliding window according to the first frequency difference and the preset resonance threshold when the first frequency difference is less than the preset resonance threshold; and to restore the length of the sliding window to the preset window length when the first frequency difference is greater than or equal to the preset resonance threshold.

[0047] When the first frequency difference is less than a preset resonance threshold, the window length calculation unit extends the length of the sliding window according to the ratio of the first frequency difference to the preset resonance threshold, i.e., sliding window length = preset window length × [1 + preset smoothing coefficient × (preset resonance threshold - first frequency difference)]. When the first frequency difference is greater than or equal to the preset resonance threshold, the sliding window length returns to the preset window length. In this embodiment, the sliding window length is dynamically updated every 0.5 seconds to ensure a rapid response to changes in wind conditions.

[0048] A preset smoothing coefficient is used to linearly map the deviation between the first frequency difference and the preset resonance threshold to the extension factor of the window length. Its value is determined to ensure that, within the typical wind field fluttering frequency fluctuation range, the window length linearly increases as the frequency difference decreases. In this embodiment, the preset window length is 5 seconds, the preset resonance threshold is 0.2 Hz, and the preset smoothing coefficient is set to 5, with units of Hz. -1 When the first frequency difference decreases from 0.2Hz to 0, the window length increases linearly from 5 seconds to 10 seconds. This value is less than the preset maximum window length of 30 seconds, so the actual range of window length variation is limited to between 5 seconds and 10 seconds, achieving a smooth transition.

[0049] It should be noted that, to prevent the sliding window length from being excessively extended due to an excessively small first frequency difference, thus affecting the real-time performance of the attitude change rate calculation, a maximum allowable value for the sliding window length is set in this embodiment. When the sliding window length calculated by the formula exceeds the preset maximum window length, it is limited to the preset maximum window length. The setting of the preset maximum window length takes into account the response speed of the beam's attitude change, the timeliness of the hazard trend determination, and the frequency resolution requirements required for frequency domain analysis. This ensures that wind-induced flutter noise can be effectively suppressed under extremely high resonance risk, while avoiding lag in the update of the tilt angle change rate due to an excessively long window, thus ensuring the real-time performance of the hazard trend determination. In this embodiment, the preset maximum window length is set to 30 seconds.

[0050] The preset resonance threshold is a critical difference used to determine whether the buffeting dominant frequency and the natural frequency of the beam are close. The value is determined by comprehensively considering the typical bandwidth of the buffeting energy concentration frequency band in the canyon wind field (0.1 to 0.3 Hz) and the sensitivity of the beam structure's damping ratio to resonance amplification. In this embodiment, based on structural dynamics characteristics and finite element simulation analysis, the preset resonance threshold is set to 0.2 Hz. This provides a clear frequency domain trigger condition for adjusting the sliding window length, ensuring that the window extension operation is only initiated within the dangerous frequency band related to the structural dynamic characteristics. This maintains the stability of the attitude change rate calculation while effectively suppressing wind-induced noise.

[0051] The preset window length refers to the base value of the time window used for calculating the tilt change rate when no resonance risk is triggered. The value is determined by comprehensively considering the response speed of the beam's attitude change during the jacking construction and the real-time requirements of the control system. In this embodiment, based on the tilt sensor sampling frequency of 20Hz and the frequency distribution characteristics of the valley wind field buffeting, the initial length of the sliding window is set to 5 seconds, corresponding to 100 sampling points. This initial value can effectively suppress random noise under normal wind conditions, while ensuring that the calculation delay of the attitude change rate does not exceed 5 seconds, meeting the timeliness requirement for dangerous trend judgment.

[0052] The sliding window length is continuously adjusted by the ratio of the first frequency difference to the preset resonance threshold. When the first frequency difference is less than the threshold, the window length is extended proportionally, and the extension factor is inversely proportional to the frequency proximity, thereby achieving adaptive smoothing of wind-induced flutter noise. When the difference is greater than or equal to the threshold, the window returns to the preset length to avoid unnecessary response delay and to match the filtering strength with the resonance risk in real time. This ensures both rapid response under normal operating conditions and sufficient noise suppression in the resonance danger frequency band.

[0053] Please see Figure 3 As shown, it is the logic diagram for determining whether the pushed beam is about to become unstable in this embodiment.

[0054] Specifically, the instability module includes: The static load moment calculation submodule is used to calculate the total jacking force on the outer side and the total jacking force on the inner side according to the pressure, and to calculate the stabilizing moment and the static load overturning moment in combination with the tilt angle. The corrected moment calculation submodule is used to calculate the wind-induced overturning moment based on the wind speed; The stability index calculation submodule is used to calculate the stability index based on the stabilizing moment, the static overturning moment, and the wind-induced overturning moment. The trend calculation submodule is used to determine the trend based on the stability index sequence within a preset stable period; The instability trend determination submodule is used to determine whether the pushed beam is about to become unstable based on the stability index and the change trend.

[0055] The stabilizing torque is calculated using the following formula: Mstable = Fouttotal × Louteq + G × (B / 2), where Mstable is the stabilizing torque, Fouttotal is the total jacking force obtained by multiplying the pressure values ​​measured by the pressure sensors of each outer jack by the effective area of ​​their respective pistons, Louteq is the equivalent force arm of the outer jack calculated using the following formula, G is the weight of the beam, and B is the width of the beam.

[0056] Where i is the index variable of the jack, N is the number of outer jacks, Fout,i is the jacking force of the i-th jack, which is calculated by multiplying the pressure value measured by the pressure sensor at that location by the effective area of ​​the corresponding piston, and Lout,i is the lever arm value of the i-th outer jack, which is determined in the design stage according to the actual arrangement of the jacks. It is a fixed constant, obtained from the bridge design drawings and pre-stored in the control system.

[0057] The static overturning moment is calculated according to the following formula: Moverturn = Final × Lineq + G × h × tanθ, where Moverturn is the static overturning moment, Final is the total inner thrust obtained by multiplying the pressure values ​​measured by the pressure sensors of each inner jack by the effective area of ​​their respective pistons, Lineq is the inner equivalent force arm, the calculation of which is the same as that of the outer equivalent force arm, and will not be repeated here; G is the weight of the beam, h is the height of the beam center, and θ is the tilt angle.

[0058] The corrected moment calculation submodule calculates the additional overturning moment generated by the wind load on the jacked beam based on the real-time wind speed, obtaining the wind-induced overturning moment, which serves as a correction term for the static load overturning moment. The calculation formula is as follows:

[0059] Where Mwind is the wind-induced overturning moment, ρ is the air density (1.225 kg / m³ in this embodiment), Cd is the wind load shape coefficient (1.3 according to the "Specifications for Wind-Resistant Design of Highway Bridges" in this embodiment), A is the windward area of ​​the beam, which is pre-stored in the system according to the beam design dimensions; v is the real-time wind speed; hwind is the vertical distance from the centroid of the windward area in the beam design drawings to the overturning rotation center, which is a fixed constant, obtained from the bridge design documents and entered into the control system. Used to convert rice from 0.05 Nb to 1 kNb.

[0060] The stability index is calculated in real time as follows: the total overturning moment is obtained by summing the static overturning moment and the wind-induced overturning moment, and the stability index is obtained by dividing the stability moment by the total overturning moment.

[0061] Understandably, the stability index calculation formula in this embodiment adopts the following simplified assumption: the wind load is only included in the denominator in the rectangular form of the wind-induced overturning force it generates, without separately correcting the change in the outer jacking force caused by the wind load in the stabilizing moment. This is because the total outer jacking force is measured in real time by pressure sensors, naturally including the actual impact of the wind load on the support reaction force. Therefore, the stabilizing moment itself implicitly reflects this coupling effect; further correction of the stabilizing moment would lead to redundant calculations. This simplification is biased towards safety: when the wind load causes a decrease in the outer jacking force, the stabilizing moment automatically decreases, while the denominator increases the wind-induced overturning moment, resulting in a slightly larger decrease in the stability index than the actual physical process, thus providing an early warning from the system; when the wind load causes an increase in the outer jacking force, the stabilizing moment automatically increases, but the denominator still increases the wind-induced overturning moment, and the stability index may be slightly higher than the actual safety margin. However, considering that strong winds mainly produce an overturning effect during jacking construction, this error is within an acceptable range.

[0062] The preset stabilization duration refers to the length of the time window used to determine the changing trend of the continuous stability index sequence. The value is determined by comprehensively considering the stability index update frequency and the actual response speed of the push attitude change. In this embodiment, the preset stabilization duration is set to 10 seconds, corresponding to 100 sampling points. This effectively suppresses high-frequency fluctuations through linear regression and promptly captures the continuous downward trend of the stability index under wind-induced disturbances.

[0063] The trend calculation submodule performs linear regression on the stability index sequence continuously collected within a preset stable period and uses the slope of the fitted line as the trend. When the slope is negative, it indicates that the stability index is continuously decreasing and the safety margin is deteriorating. The absolute value of the slope reflects the rate of deterioration.

[0064] The preset safety threshold is a critical value used to determine whether the stability index has fallen to the danger zone. The value is determined by comprehensively referencing the requirements for the anti-overturning stability coefficient in the "Crane Design Code" and the dynamic characteristics of canyon wind-driven jacking construction. In this embodiment, based on finite element simulation and on-site trial jacking calibration, the preset safety threshold is set to 1.2. This ensures that an emergency response is triggered only when the safety margin is truly insufficient and continues to deteriorate, avoiding missed alarms due to an excessively low threshold or frequent false shutdowns due to an excessively high threshold. Those skilled in the art can routinely adjust the threshold within the range of 1.1 to 1.5 according to the specific bridge span, beam stiffness, and wind load conditions.

[0065] The instability trend determination submodule comprehensively determines whether the pushed beam is about to become unstable based on the comparison between the current stability index and the preset safety threshold, as well as the sign and absolute value of the changing trend. Specifically, when the stability index is lower than the preset safety threshold and the changing trend is negative, the beam is determined to be about to become unstable; when the stability index is lower than the preset safety threshold but the changing trend is non-negative, only a safety margin insufficient warning is issued, and it is not determined to be about to become unstable, allowing the system to maintain the current state and strengthen monitoring; when the stability index is higher than the preset safety threshold but the changing trend is negative and its absolute value exceeds the preset rapid deterioration threshold, it indicates that although the safety margin has not yet reached the critical value, it is decaying at a relatively fast rate. To prevent sudden instability, it is also determined to be about to become unstable, and the lateral active correction in the hierarchical module is triggered in advance to suppress the deterioration trend. In this embodiment, the preset rapid deterioration threshold is set to 0.1 / s, and those skilled in the art can routinely adjust it within the range of 0.05 / s to 0.2 / s according to the actual dynamic characteristics of the beam and the safety redundancy.

[0066] The static load moment calculation submodule calculates the stabilizing moment and static overturning moment based on real-time jacking force and tilt angle, accurately quantifying the fundamental contribution of self-weight and jacking force to the overturning margin from a static perspective. The modified moment calculation submodule introduces real-time wind speed to calculate wind-induced overturning moment, superimposing the dynamic load of the canyon wind field onto the total overturning moment, making up for the deficiency of traditional methods in quantifying the impact of environmental disturbances. The stability index calculation submodule uses the ratio of the stabilizing moment to the total overturning moment as the stability index, making this index directly correspond to the overturning safety margin in a physical sense. The trend calculation submodule extracts the time-domain trend of the stability index through linear regression within a sliding window, accurately capturing the decay rate of the safety margin. The instability trend judgment submodule integrates the current stability index and its rate of change, adopting multiple criteria such as low margin and continuous deterioration and high margin but rapid decay, realizing a substantial upgrade from static threshold warning to dynamic trend warning, significantly improving the lead time and accuracy of instability warning under strong disturbance environments such as canyon wind fields, and providing a reliable decision-making basis for graded control.

[0067] Specifically, the hierarchical module includes: The risk level determination submodule is used to determine the current risk level based on the stability index and the trend of change. The risk level includes low risk, medium risk and high risk. The lateral correction execution submodule is used to select lateral active correction when the risk level is medium risk.

[0068] The grading module classifies the current state into three levels—low risk, medium risk, and high risk—based on the comparison between the stability index and a preset safety threshold, as well as the sign and absolute value of the changing trend: Low risk corresponds to a stability index greater than the preset safety threshold and a non-negative changing trend; medium risk corresponds to a stability index greater than the preset safety threshold but a negative changing trend and an absolute value greater than the preset rapid deterioration threshold, or a stability index less than or equal to the preset safety threshold and a non-negative changing trend; high risk corresponds to a stability index less than or equal to the preset safety threshold and a negative changing trend. When a medium risk determination is established, the lateral correction execution submodule initiates lateral active correction to prevent the risk from escalating to high risk and triggering an emergency shutdown. The definitions and values ​​of the preset safety threshold and the preset rapid deterioration threshold are the same as described in the aforementioned instability module and will not be repeated here.

[0069] The grading module determines the current risk level based on the stability index and its changing trend, and executes corresponding actions for different risk levels: when the risk level is low, no additional actions are triggered, the system maintains the current jacking state and continues monitoring; when the risk level is medium, the lateral correction execution submodule initiates active lateral correction, controlling the auxiliary jacks arranged on the side of the beam to apply correction force to limit lateral displacement; when the risk level is high, an emergency stop is executed, all hydraulic valves are immediately locked, and an audible and visual alarm is issued. The risk level determination submodule maps the stability index and its changing trend into three levels of risk: low, medium, and high. It initiates proactive lateral correction only at medium risk, maintains monitoring at low risk, and executes an emergency shutdown at high risk, achieving risk-driven differentiated control. This avoids construction interruptions caused by accidental corrections at low risk, ensures timely suppression of lateral displacement at medium risk to prevent escalation, and reserves a bottom line for decisive shutdown at high risk. Thus, it maximizes the continuity and efficiency of the jacking operation while ensuring structural safety.

[0070] Specifically, the lateral correction execution submodule includes: Lateral offset calculation unit, which is used to calculate the current lateral offset based on the tilt angle and beam geometry parameters; The correction force determination unit is used to determine the required lateral correction force based on the lateral offset and the preset correction gain coefficient. The lateral push execution unit is used to control the auxiliary jacks arranged on the side of the beam to apply the lateral correction force until the updated lateral offset is restored to the preset allowable range.

[0071] The lateral offset calculation unit calculates the current lateral offset based on the real-time tilt angle θ and the center of gravity height h and rotation center position in the beam's geometric parameters, using the formula Δ=h×tanθ, where Δ is the horizontal offset distance of the beam's center of gravity relative to the overturning rotation center. When the tilt angle is small, tanθ≈θ, and a linear approximation Δ=h×θ can also be used to simplify the calculation.

[0072] The preset correction gain coefficient refers to the proportional coefficient that converts the lateral offset into the required lateral correction force. The value is determined by comprehensively considering the lateral stiffness of the beam, the response characteristics of the hydraulic system, and the safety requirements of the correction. In this embodiment, based on finite element simulation and on-site debugging, the preset correction gain coefficient is set to 0.5 kN / mm, establishing a linear mapping relationship between the offset and the correction force. This allows the correction force to adapt to the offset, ensuring gentle correction for small offsets and rapid response for large offsets, while avoiding over-adjustment or under-adjustment due to fixed force output. Those skilled in the art can adjust it within the range of 0.2 kN / mm to 1.0 kN / mm according to the actual beam stiffness and hydraulic system characteristics.

[0073] The preset dead zone threshold refers to the allowable range of no lateral offset. The value is determined by comprehensively considering the accuracy of the tilt sensor and the structural rigidity, to avoid frequent start-stop of the correction unit due to sensor noise or minor wind-induced vibrations. In this embodiment, the preset dead zone threshold is set to 2mm, effectively filtering out interference from minor offsets and sensor noise on the correction action, and preventing increased energy consumption and component wear caused by frequent start-stop of the hydraulic system. Those skilled in the art can adjust it within the range of 1mm to 5mm according to the actual sensor accuracy and structural rigidity.

[0074] The correction force determination unit employs a proportional control strategy with a dead zone to determine the lateral correction force. When the absolute value of the current lateral offset does not exceed a preset dead zone threshold, it is determined to be an ineffective offset, and no correction force is output. When the lateral offset exceeds the preset dead zone threshold, the magnitude of the correction force is equal to the value of the lateral offset exceeding the dead zone multiplied by a preset correction gain coefficient, and the direction of the correction force is opposite to the direction of the lateral offset. This ensures that small offsets will not trigger invalid correction, and once the offset exceeds the allowable range, the correction force increases linearly with the offset, achieving smooth and adaptive active correction.

[0075] It should be noted that, to achieve hysteresis control of the offset between the preset dead zone threshold and the preset allowable range, the lateral correction in this embodiment employs a dual-threshold hysteresis comparator: the start threshold is set to the preset dead zone threshold, and the stop threshold is set to the preset allowable range. When the absolute value of the lateral offset is greater than the start threshold, correction is initiated and a correction force is output; when correction is in progress and the absolute value of the lateral offset is less than the stop threshold, correction is stopped; when the absolute value of the lateral offset is between the stop threshold and the start threshold, the current correction state remains unchanged, i.e., if correction is in progress, the correction force continues to be output; if correction is not in progress, no output is output. This hysteresis control mechanism effectively avoids frequent start-stop cycles caused by signal fluctuations near a single threshold, ensuring the smoothness of the correction action and the service life of the actuator.

[0076] The preset allowable range refers to the target range to which the lateral offset should recover after correction. The value is determined by comprehensively considering the beam's lateral stiffness, construction accuracy requirements, and sensor noise levels, ensuring that the offset is below this value and thus considered corrected, avoiding over-correction. In this embodiment, based on engineering experience and finite element analysis, the preset allowable range is set to 1mm, which, together with the preset dead zone threshold of 2mm, forms a hysteresis interval. This ensures that the current correction state remains unchanged when the offset is between 1mm and 2mm, effectively suppressing frequent starts and stops and hydraulic system oscillations caused by measurement noise or minor fluctuations. Those skilled in the art can routinely adjust the range to 0.5mm to 2mm according to actual accuracy requirements, and the preset allowable range should be less than the preset dead zone threshold.

[0077] By filtering out minute offsets and sensor noise through a preset dead zone threshold, the frequency of invalid actions of the hydraulic actuator is reduced. The preset correction gain coefficient realizes a linear proportional relationship between the correction force and the effective offset, enabling the correction force to be adaptively adjusted as needed. The preset dead zone threshold is greater than the preset allowable range, avoiding critical oscillations caused by signal fluctuations near a single threshold, thus improving the stability of the control system. While ensuring that the lateral offset is effectively suppressed, the number of start-stop cycles of the hydraulic system is significantly reduced, energy consumption and mechanical wear are reduced, and the smoothness of the correction action and the service life of the actuator are improved.

[0078] Specifically, the execution module includes: The instruction generation submodule is used to generate corresponding hydraulic control instructions based on the lateral active correction action output by the lateral correction execution submodule. The instruction output submodule is used to output the hydraulic control instructions to the hydraulic actuator of the jacking operation.

[0079] The instruction generation submodule, based on the lateral correction force determined by the lateral active correction action, maps the correction force to the corresponding proportional valve drive current value or PWM duty cycle signal according to the current-force characteristic curve of the hydraulic proportional valve, thereby generating hydraulic control instructions. The instruction output submodule sends the hydraulic control instructions in real time to the hydraulic actuator of the corresponding lateral auxiliary jack in the jacking construction through a fieldbus, such as CAN bus, Profinet or EtherCAT, so as to drive it to output the required correction force.

[0080] By precisely mapping the correction force requirement in the lateral active correction action to a proportional valve drive current or PWM duty cycle signal, and using the command output submodule to send control commands to the hydraulic actuator unit without delay via a real-time industrial fieldbus, it is ensured that the lateral auxiliary jack can output correction force accurately and on demand. This guarantees the accuracy of control commands and the real-time transmission, avoiding correction lag or over-adjustment caused by signal conversion errors or communication delays, thus providing a reliable execution guarantee for the rapid response and stable control of the entire system.

[0081] Please see Figure 4 As shown, it is a flowchart of adjusting the preset deviation threshold in this embodiment.

[0082] Specifically, the adjustment modules include: The frequency deviation calculation submodule is used to calculate the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and obtain the second frequency difference. The time correlation calculation submodule is used to determine the time correlation based on the number of deviation events and the number of dangerous events. The number of deviation events is the total number of times the deviation judgment result indicates that there is a deviation within the preset observation period, and the number of dangerous events is the number of times the danger judgment result indicates that there is a danger within the preset correlation window after the deviation judgment result indicates that there is a deviation. The threshold adjustment submodule is used to adjust the preset deviation threshold based on whether the second frequency difference is less than the preset resonance threshold and whether the time correlation is greater than the preset correlation threshold.

[0083] At the end of each preset observation period, the frequency deviation calculation submodule obtains the peak frequency of the wind-induced buffeting spectrum and the natural frequency of the pushed beam that was last updated in that period, respectively, and subtracts the two directly. The difference is recorded as the second frequency difference.

[0084] The preset correlation window refers to the time range used to determine whether a dangerous event has occurred after a deviation event. Its value is determined by comprehensively considering the sliding window length of the tilt change rate in the danger module and the calculation time of the stable change rate, combined with the physical evolution lag time from deviation to danger. In this embodiment, the sliding window length of the tilt change rate and the calculation time of the stable change rate are 5 seconds and 10 seconds, respectively. The physical evolution lag time from deviation to danger is typically 10–30 seconds. Therefore, the preset correlation window is set to 30 seconds to cover most processes from synchronous deviation to dangerous trends, while avoiding the introduction of irrelevant noise due to an excessively long window. This ensures that the deviation threshold adjustment can truly reflect the correlation between the deviation event and subsequent danger.

[0085] At the end of each preset observation period, the frequency deviation calculation submodule obtains the peak frequency of the most recently updated wind-induced flutter spectrum and the natural frequency of the most recently updated pushed beam within the preset observation period, and calculates the absolute value of the difference between the two as the second frequency difference.

[0086] The time correlation calculation submodule calculates the total number of deviation events within each preset observation period, and counts the number of dangerous events within the preset correlation window after each deviation event. The ratio of the number of dangerous events to the number of deviation events is then calculated, and the resulting ratio is the time correlation.

[0087] The preset resonance threshold is a critical difference used to determine whether the dominant frequency of wind-induced buffeting is close to the natural frequency of the beam. The value is determined by comprehensively considering the typical bandwidth of the concentrated buffeting energy in the canyon wind field (0.1–0.3 Hz) and the sensitivity of the beam structure's damping ratio to resonance amplification. In this embodiment, based on structural dynamics characteristics and finite element simulation analysis, the preset resonance threshold is set to 0.2 Hz. This ensures that the adjustment of the sliding window length and the correction of the deviation threshold only activate within the dangerous frequency band where the wind-induced excitation is sufficiently close to the structure's natural frequency, avoiding frequent triggering of adaptive adjustment due to meaningless frequency fluctuations. In this embodiment, the preset resonance threshold used by the adjustment module and the preset resonance threshold used by the danger module are the same value; however, depending on the actual working conditions, they can also be calibrated independently.

[0088] The preset correlation threshold refers to the critical probability value used to determine whether the time correlation degree has reached the deviation threshold adjustment condition. In this embodiment, based on the statistical analysis of historical data of similar bridge jacking construction, when the proportion of dangers occurring after a deviation event exceeds 60%, the leniency of the deviation threshold has already posed a significant risk to construction safety. At the same time, by combining finite element simulation and Monte Carlo method, risk-false alarm trade-off analysis is carried out under different thresholds, and it is determined that setting the preset correlation threshold to 0.6 can balance the early warning sensitivity and stability.

[0089] The frequency deviation calculation submodule monitors in real time the proximity of the wind-induced buffeting frequency to the natural frequency of the beam. The time correlation calculation submodule calculates the conditional probabilities of deviation events and hazardous events. Then, the threshold adjustment submodule dynamically adjusts the preset deviation threshold based on a dual criterion of frequency domain resonance risk and time domain causality. When the buffeting frequency approaches the natural frequency, the resonance risk is considered to have increased, and the deviation threshold is proactively lowered to tighten the deviation tolerance and improve early warning sensitivity. When deviation events frequently induce subsequent hazards, the current deviation threshold is considered too lenient, and the threshold is also lowered to strengthen forward warning. This dual-condition fusion mechanism overcomes the limitations of a single criterion, enabling rapid response to sudden wind changes and adaptive adjustment based on long-term construction data. This ensures that deviation judgment remains in an optimally sensitive state, effectively solving the technical problems of missed alarms during periods of strong disturbance and false alarms during periods of weak disturbance with fixed thresholds. It achieves adaptive matching of control parameters with the construction environment and structural state.

[0090] Specifically, the threshold adjustment submodule includes: A resonance risk adjustment unit is used to determine a first adjustment range based on the second frequency difference and the preset resonance threshold when the second frequency difference is less than the preset resonance threshold. A statistical correlation adjustment unit is used to determine a second adjustment range based on the difference between the time correlation degree and the preset correlation threshold when the time correlation degree is greater than the preset correlation threshold. The integrated adjustment unit is used to adjust the preset deviation threshold according to the first adjustment range and the second adjustment range.

[0091] The preset maximum adjustment range refers to the maximum allowable change ratio for a single threshold adjustment. The value is determined by comprehensively considering system stability requirements and adaptive response speed. In this embodiment, based on finite element simulation and on-site debugging, the preset maximum adjustment range is set to 0.1. This provides sufficient adjustment sensitivity while ensuring system stability, allowing the threshold to converge smoothly and quickly to the optimal value, while avoiding control jitter caused by over-adjustment.

[0092] When the resonance risk adjustment unit determines that the second frequency difference is less than the preset resonance threshold, the first adjustment amplitude is calculated using the following formula:

[0093] Where a1 is the first adjustment amplitude, amax is the preset maximum adjustment amplitude, Δf is the second frequency difference, and Δfth is the preset resonance threshold.

[0094] When the correlation degree in the statistical correlation adjustment unit is determined to be greater than the preset correlation threshold, the second adjustment range is calculated using the following formula:

[0095] Where a2 is the second adjustment range, amax is the preset maximum adjustment range, R is the time correlation degree, and Rth is the preset correlation threshold.

[0096] It is understandable that as R approaches 1, deviation events almost inevitably lead to danger, and the adjustment range is larger; when R ≤ Rth, a2 = 0.

[0097] The preset minimum allowable deviation threshold refers to the lowest limit at which the preset deviation threshold allowed by the system can be lowered. In this embodiment, the preset minimum allowable deviation threshold is 0.5 times the initial threshold, which is 0.4. The setting of this lower limit value includes the noise amplitude of the tilt sensor measurement, the response dead zone of the hydraulic actuator, and the structural safety margin, ensuring that the threshold is not lower than the minimum value required to maintain basic anti-interference capability. While adaptively lowering the threshold based on resonance risk and statistical correlation to improve early warning sensitivity, a lower limit protection is set for threshold adjustment to avoid excessive lowering that would cause the system to become too sensitive and generate frequent false alarms, thereby achieving a balance between sensitivity and stability.

[0098] The comprehensive adjustment unit weights and merges the first adjustment range and the second adjustment range to obtain the total adjustment range. The weight of the resonance risk adjustment range is set to 0.6, and the weight of the statistical correlation adjustment range is set to 0.4 to highlight the priority of resonance risk. The updated preset deviation threshold is calculated according to the following formula: new threshold = original threshold × (1 - total adjustment range). At the same time, a lower limit protection is set to ensure that the new threshold is not lower than the preset minimum allowable deviation threshold.

[0099] The first adjustment amplitude is determined based on the proximity of the fluttering dominant frequency to the natural frequency; the smaller the frequency difference, the higher the resonance risk and the larger the adjustment amplitude, thus enhancing the sensitivity driven by the physical mechanism. The second adjustment amplitude is determined based on the conditional probability of deviation events and dangerous events; the higher the correlation, the more significant the leniency of the threshold and the larger the adjustment amplitude, thus achieving optimization and correction driven by historical data. A weighted fusion strategy is used to integrate the two to uniformly update the preset deviation threshold. By organically combining the structural dynamics resonance principle with statistical learning laws, the adjustment of the deviation threshold can respond quickly to sudden changes in wind conditions and continuously optimize based on long-term construction data, effectively avoiding the limitations of a single criterion and achieving a dynamic balance between control sensitivity and environmental adaptability. This overcomes the shortcomings of fixed thresholds, which are prone to missed detections during periods of strong disturbance and false detections during periods of weak disturbance.

[0100] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A hydraulic synchronous jacking construction control system based on multi-source information fusion, characterized in that, include: The data acquisition module is used to acquire in real time the maximum displacement difference between each jacking point, the pressure at each jacking point, the inclination angle, stress coefficient and natural frequency of the jacked beam, and the wind speed in the jacking construction area during the canyon wind field jacking construction. The deviation module is used to determine whether there is a deviation in the pushed beam based on the comparison result of the deviation index and the preset deviation threshold, wherein the deviation index is determined based on the displacement difference and the pressure. The danger module is used to determine whether the pushed beam is in danger when there is a deviation, based on the rate of change of the tilt angle and the stress coefficient within the sliding window. The length of the sliding window is adjusted according to the wind-induced flutter spectrum and the natural frequency. The wind-induced flutter spectrum is determined according to the wind speed and the natural frequency. The instability module is used to determine whether the pushed beam is about to become unstable when a danger occurs, based on the stability index and its changing trend. The stability index is determined based on the pressure, the tilt angle and the correction torque, and the correction torque is determined based on the wind speed. A grading module is used to determine the execution action of the jacking construction based on the stability index and the trend of change when instability is imminent. An execution module, which outputs execution commands to the hydraulic execution unit for jacking construction based on the execution action; The adjustment module is used to adjust the preset deviation threshold based on the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and the time correlation between the deviation judgment result and the danger judgment result.

2. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 1, characterized in that, The deviation module includes: The deviation index calculation submodule is used to calculate the deviation index based on the maximum displacement difference and the pressure at each jacking point; The threshold comparison submodule is used to compare the deviation index with a preset deviation threshold to determine whether a deviation exists.

3. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 2, characterized in that, Hazardous modules include: The wind-induced buffeting spectrum calculation submodule is used to calculate the wind-induced buffeting spectrum based on the time-series data of the wind speed and extract the dominant frequency of the buffeting. A sliding window dynamic adjustment submodule is used to adjust the length of the sliding window according to the difference between the dithering main frequency and the natural frequency; The tilt angle change rate calculation submodule is used to perform differential calculation on the tilt angle within the sliding window to obtain the tilt angle change rate. The danger trend determination submodule is used to compare the tilt angle change rate and the stress coefficient with their respective preset tilt angle change rate threshold and preset stress coefficient threshold. When the tilt angle change rate is greater than the preset tilt angle change rate threshold and the stress coefficient is greater than the preset stress coefficient threshold, it is determined that the pushed beam is in danger.

4. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 3, characterized in that, The sliding window dynamic adjustment submodule includes: A frequency comparison unit is used to obtain a first frequency difference based on the difference between the dithering main frequency and the natural frequency. The window length calculation unit is used to extend the length of the sliding window according to the first frequency difference and the preset resonance threshold when the first frequency difference is less than the preset resonance threshold; and to restore the length of the sliding window to the preset window length when the first frequency difference is greater than or equal to the preset resonance threshold.

5. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 1, characterized in that, The instability module includes: The static load moment calculation submodule is used to calculate the total jacking force on the outer side and the total jacking force on the inner side according to the pressure, and to calculate the stabilizing moment and the static load overturning moment in combination with the tilt angle. The corrected moment calculation submodule is used to calculate the wind-induced overturning moment based on the wind speed; The stability index calculation submodule is used to calculate the stability index based on the stabilizing moment, the static overturning moment, and the wind-induced overturning moment. The trend calculation submodule is used to determine the trend based on the stability index sequence within a preset stable period; The instability trend determination submodule is used to determine whether the pushed beam is about to become unstable based on the stability index and the change trend.

6. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 1, characterized in that, The hierarchical module includes: The risk level determination submodule is used to determine the current risk level based on the stability index and the trend of change. The risk level includes low risk, medium risk and high risk. The lateral correction execution submodule is used to select lateral active correction when the risk level is medium risk.

7. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 6, characterized in that, The lateral correction execution submodule includes: Lateral offset calculation unit, which is used to calculate the current lateral offset based on the inclination angle and beam geometry parameters; The correction force determination unit is used to determine the required lateral correction force based on the lateral offset and the preset correction gain coefficient. The lateral push execution unit is used to control the auxiliary jacks arranged on the side of the beam to apply the lateral correction force until the updated lateral offset is restored to the preset allowable range.

8. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 1, characterized in that, The execution module includes: The instruction generation submodule is used to generate corresponding hydraulic control instructions based on the lateral active correction action output by the lateral correction execution submodule. The instruction output submodule is used to output the hydraulic control instructions to the hydraulic actuator of the jacking operation.

9. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 1, characterized in that, The adjustment module includes: The frequency deviation calculation submodule is used to calculate the difference between the peak frequency and the natural frequency of the wind-induced flutter spectrum within a preset observation period, and obtain the second frequency difference. The time correlation calculation submodule is used to determine the time correlation based on the number of deviation events and the number of dangerous events. The number of deviation events is the total number of times the deviation judgment result indicates that there is a deviation within the preset observation period, and the number of dangerous events is the number of times the danger judgment result indicates that there is a danger within the preset correlation window after the deviation judgment result indicates that there is a deviation. The threshold adjustment submodule is used to adjust the preset deviation threshold based on whether the second frequency difference is less than the preset resonance threshold and whether the time correlation is greater than the preset correlation threshold.

10. The hydraulic synchronous jacking construction control system based on multi-source information fusion according to claim 9, characterized in that, The threshold adjustment submodule includes: A resonance risk adjustment unit is used to determine a first adjustment range based on the second frequency difference and the preset resonance threshold when the second frequency difference is less than the preset resonance threshold. A statistical correlation adjustment unit is used to determine a second adjustment range based on the difference between the time correlation degree and the preset correlation threshold when the time correlation degree is greater than the preset correlation threshold. The integrated adjustment unit is used to adjust the preset deviation threshold according to the first adjustment range and the second adjustment range.

Images