UHPC wet joint construction quality intelligent control method and system
By using frequency domain analysis and transient detection technology, construction vibration and environmental disturbance were separated, and the parameters of the vibrating equipment were adjusted. This solved the problem of misjudgment in the construction quality control of UHPC wet joints in bridge reconstruction and expansion projects, and achieved precise vibration operation and material uniformity.
Patent Information
- Application Number
- CN202511544272.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-10-28
- Publication Date
- 2025-11-28
- Estimated Expiration
- 2045-10-28
AI Technical Summary
In bridge reconstruction and expansion projects, existing technologies make it difficult to distinguish between construction vibration and environmental vibration during the construction of UHPC wet joints, leading to misjudgments in the control system and risks to quality control.
By using frequency domain analysis and transient detection technology, construction characteristics and environmental disturbance characteristics are separated, construction vibration or environmental vibration is identified, and the parameters of the vibrating equipment are adjusted according to the structural coupling degree and vibration stability to achieve precise control.
It improves the reliability of vibration signal analysis, ensures the accuracy of control commands, enhances the density and uniformity of UHPC material in the joint area, and guarantees the overall performance of the bridge structure.
Smart Images

Figure CN121028569A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of automatic control during construction, and more particularly, to a UHPC wet joint construction quality intelligent control method and system. BACKGROUND
[0002] In bridge reconstruction and expansion projects, in order to ensure the effective connection between the new and old structures, the wet joint technology of ultra-high performance concrete (UHPC) is often used. In order to ensure the pouring density and material uniformity of the joint, the key procedures such as vibration need to be controlled during construction. The existing technology usually relies on preset construction parameters and sensors to monitor the vibration state, aiming to maintain the stability of the construction operation through feedback information.
[0003] However, the above-mentioned existing technology has significant defects in actual application. Since the wet joint construction is often under the interference of complex environments such as existing traffic loads, the vibration signals collected by the monitoring system are actually mixed responses of construction operation and external environmental vibration, and the existing method cannot effectively distinguish the signal sources, which may cause the system to misjudge harmless environmental disturbances as abnormal construction states, and thus generate incorrect control instructions, not only interfering with the normal construction rhythm, but also directly affecting the final forming quality of the joint, constituting a potential risk in quality control. SUMMARY
[0004] In order to overcome the above-mentioned defects of the existing technology, the present application provides a UHPC wet joint construction quality intelligent control method and system to solve the problems raised in the background technology.
[0005] To achieve the above-mentioned purpose, the present application provides the following technical solutions:
[0006] A UHPC wet joint construction quality intelligent control method, comprising the following steps:
[0007] S1, collecting vibration signals of a plurality of preset positions of the wet joint structure;
[0008] S2, extracting a first type of frequency domain features related to construction operation and a second type of frequency domain features related to environmental disturbance from the vibration signals, respectively;
[0009] S3, identifying the vibration signals as construction vibration or environmental vibration based on the energy distribution of the first type of frequency domain features and the transient characteristics of the second type of frequency domain features;
[0010] S4, when the vibration signals are identified as construction vibration, calculating the energy values of the vibration signals of the plurality of preset positions and analyzing the decay gradient of the energy values with the increase of the distance between the preset positions and the vibration source, and evaluating the structure coupling degree according to the steepness of the decay gradient;
[0011] S5, based on the structural coupling degree, analyze the concentration of the first type of frequency domain features to evaluate the stability of the construction vibration, and adjust the control parameters of the vibrating equipment according to the stability evaluation result;
[0012] S6, execute the adjusted control parameters to control the vibrating equipment to vibrate the wet joint.
[0013] Further, the vibration signals of the wet joint structure at multiple preset positions are collected, including:
[0014] The vibration sensors are arranged at the preset positions at the new and old concrete combination interface of the wet joint structure and at the preset positions in the action area of the vibrating equipment, and the vibration signals are synchronously collected through the vibration sensors.
[0015] Further, from the vibration signals, the first type of frequency domain features related to the construction operation and the second type of frequency domain features related to the environmental disturbance are extracted, including:
[0016] The vibration signals are subjected to frequency spectrum analysis to obtain the frequency component distribution of the vibration signals;
[0017] Based on the frequency component distribution, the main band energy is calculated as the first type of frequency domain features;
[0018] Meanwhile, it is detected whether there is a transient pulse component in the vibration signal as the second type of frequency domain features.
[0019] Further, based on the energy distribution of the first type of frequency domain features and the transient characteristics of the second type of frequency domain features, the vibration signal is identified as construction vibration or environmental vibration, including:
[0020] It is judged whether the main band energy in the first type of frequency domain features exceeds the first threshold value and the transient pulse component amplitude in the second type of frequency domain features is lower than the second threshold value;
[0021] When the main band energy exceeds the first threshold value and the transient pulse component amplitude is lower than the second threshold value, the vibration signal is identified as construction vibration; otherwise, the vibration signal is identified as environmental vibration.
[0022] Further, when the vibration signal is identified as construction vibration, the energy values of the vibration signals at multiple preset positions are calculated and the decay gradient of the energy values with the increase of the distance between the preset positions and the vibration source is analyzed, and the structural coupling degree is evaluated according to the steepness of the decay gradient, including:
[0023] The energy value of the vibration signal at each preset position is calculated, and a relationship model of the energy value and the distance from the corresponding preset position to the vibration source is established;
[0024] The decay gradient of the energy value with the distance is determined through the relationship model;
[0025] The structure coupling degree is evaluated based on linearity of the decay gradient.
[0026] The higher the linearity is, the better the structure coupling degree is.
[0027] Further, the energy value of each preset position vibration signal is calculated, and the relationship model of the energy value and the distance from the corresponding preset position to the vibration source is established by: calculating the energy value by square integration of the vibration signal of each preset position, and linearly fitting the energy value and the distance from the corresponding preset position to the vibration source to establish the relationship model.
[0028] Further, the structure coupling degree is evaluated based on linearity of the decay gradient, including: calculating the residual sum of squares of the data points and the fitting straight line in the relationship model to evaluate the linearity, and the smaller the residual sum of squares is, the higher the structure coupling degree is.
[0029] Further, the stability of the construction vibration is evaluated by analyzing the concentration degree of the first type of frequency domain features based on the structure coupling degree, and the control parameters of the vibrating equipment are adjusted according to the stability evaluation result, including:
[0030] The stability evaluation index is established based on the structure coupling degree and the discrete degree of the energy distribution of the main frequency band in the first type of frequency domain features;
[0031] When the stability evaluation index is lower than the set threshold, the vibration intensity control parameter of the vibrating equipment is increased;
[0032] When the stability evaluation index is not lower than the set threshold, the vibration intensity control parameter of the vibrating equipment is maintained or reduced.
[0033] Further, the adjusted control parameters are executed to control the vibrating equipment to perform the vibrating operation on the wet joint, including:
[0034] The adjusted vibration intensity control parameter is output to the controller of the vibrating equipment;
[0035] The vibrating equipment is driven by the controller to perform the vibrating operation on the wet joint at the adjusted vibration intensity control parameter.
[0036] On the other hand, the present application provides a UHPC wet joint construction quality intelligent control system, including the following modules:
[0037] A signal acquisition module is used to acquire vibration signals of a plurality of preset positions of the wet joint structure;
[0038] A feature extraction module is used to extract a first type of frequency domain features related to construction operation and a second type of frequency domain features related to environmental disturbance from the vibration signals, respectively;
[0039] The vibration identification module is configured to identify the vibration signal as the construction vibration or the environmental vibration based on energy distribution of the first type of frequency domain feature and transient characteristics of the second type of frequency domain feature.
[0040] The coupling evaluation module is configured to calculate energy values of the vibration signals at the plurality of preset positions and analyze an attenuation gradient of the energy values with an increase of distances between the preset positions and the vibration source when the vibration signal is identified as the construction vibration, and evaluate the structural coupling degree according to a steepness of the attenuation gradient.
[0041] The parameter adjustment module is configured to analyze the concentration degree of the first type of frequency domain feature based on the structural coupling degree to evaluate the stability of the construction vibration, and adjust the control parameter of the vibrating equipment according to an evaluation result of the stability.
[0042] The operation execution module is configured to execute the adjusted control parameter to control the vibrating equipment to perform the vibrating operation on the wet joint.
[0043] Compared with the prior art, the present application has the following beneficial effects:
[0044] 1. By separating the construction features and the environmental disturbance features from the mixed vibration signal, the accurate perception and closed-loop control of the vibrating process are realized. First, the construction vibration and the environmental vibration are effectively distinguished through the frequency domain analysis and the transient detection technology, the misguidance of the environmental interference such as the driving load to the control system is excluded, the system can accurately identify the real construction state, on this basis, the structural coupling state of the new and old concrete interface is evaluated in real time through the energy propagation characteristics analysis of the construction vibration signal, and then the vibrating parameter is dynamically adjusted according to the coupling quality, the reliability of the vibration signal analysis is significantly improved, the misoperation easily produced by the traditional single threshold judgment is avoided, and the accurate generation of the control instruction is ensured.
[0045] 2. Through the collaborative analysis of the structural coupling degree and the concentration degree of the vibration features, a complete quality feedback closed loop is constructed, the interface combination quality can be judged according to the attenuation law of the vibration energy in the structure, and the vibrating intensity parameter is automatically optimized accordingly, the complete automatic chain from the state monitoring to the process adjustment is realized, the limitations of the traditional dependence on the fixed parameter or the artificial experience are overcome, the vibrating operation is always in the optimal working state, the compactness and uniformity of the UHPC material in the joint area are significantly improved, and a reliable guarantee is provided for the overall performance of the bridge structure. BRIEF DESCRIPTION OF DRAWINGS
[0046] Figure 1 A flowchart of the UHPC wet joint construction quality intelligent control method of the present application;
[0047] Figure 2 A structural schematic diagram of the UHPC wet joint construction quality intelligent control system of the present application. DETAILED DESCRIPTION
[0048] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all the other embodiments obtained by a person of ordinary skill in the art without creative work are within the protection scope of the present application.
[0049] Embodiment 1 Figure 1 A UHPC wet joint construction quality intelligent control method is given, which comprises the following steps:
[0050] S1, collecting vibration signals of multiple preset positions of the wet joint structure;
[0051] S2, extracting a first type of frequency domain features related to construction operation and a second type of frequency domain features related to environmental disturbance from the vibration signals;
[0052] S3, identifying the vibration signals as construction vibration or environmental vibration based on the energy distribution of the first type of frequency domain features and the transient characteristics of the second type of frequency domain features;
[0053] S4, when the vibration signals are identified as construction vibration, calculating the energy values of the vibration signals of the multiple preset positions and analyzing the attenuation gradient of the energy values with the increase of the distance between the preset positions and the vibration source, and evaluating the structure coupling degree according to the steepness of the attenuation gradient;
[0054] S5, analyzing the concentration of the first type of frequency domain features based on the structure coupling degree to evaluate the stability of the construction vibration, and adjusting the control parameters of the vibrating equipment according to the stability evaluation result;
[0055] S6, executing the adjusted control parameters to control the vibrating equipment to perform vibrating operation on the wet joint.
[0056] S1, collecting vibration signals of multiple preset positions of the wet joint structure, which is specifically implemented as:
[0057] During the construction of the wet joint structure, the collection of vibration signals first needs to determine a plurality of preset positions, which are selected based on the geometric characteristics and stress analysis of the wet joint structure to ensure that the vibration propagation characteristics can be fully captured. Specifically, the preset positions include the new and old concrete combination interface and the vibration equipment action area, which are determined by engineering drawings and field measurement, for example, the combination interface selects a point position at a certain interval from the interface edge, and the vibration equipment action area is arranged in a ring or linear array around the vibration rod insertion point. The vibration sensor uses a piezoelectric accelerometer, which works on the principle of using the characteristic of piezoelectric material to generate an electric charge signal when subjected to mechanical vibration. The sensor is fixed to the concrete surface by epoxy resin glue to ensure close contact with the structure to reduce signal loss. When arranging the sensor, the axis of the sensor should be consistent with the direction of vibration propagation to accurately measure the vibration acceleration value.
[0058] When arranging vibration sensors at the new and old concrete combination interface, first identify the specific boundary of the combination interface, which is determined by visual inspection or non-destructive testing methods such as ultrasonic scanning, then symmetrically arrange sensors on both sides of the interface, for example, select three point positions on each side of the interface, and evenly distribute the point positions according to the interface length to ensure that the entire interface area is covered. Before installation, the sensor is calibrated to check its frequency response range and sensitivity to ensure that it works linearly within the construction vibration frequency range. During installation, avoid placing the sensor in areas with dense steel bars or aggregates to prevent signal distortion, and use a protective cover to protect the sensor from construction debris.
[0059] When arranging vibration sensors at the vibration equipment action area, first determine the working radius and insertion depth of the vibration rod, and arrange the sensor array around the vibration point according to these parameters, for example, arrange a sensor every certain distance in the radial direction with the vibration point as the center, and the distance value is set based on vibration energy attenuation experimental data to ensure that the spatial distribution of vibration energy can be monitored. The type of sensor selected is consistent with that at the combination interface to ensure data comparability, and a magnetic base or clamp is used for quick installation to meet the needs of vibration equipment movement. After arrangement, verify the signal connectivity of all sensors, and check whether the output waveform is normal by applying a test vibration.
[0060] When synchronously collecting vibration signals through vibration sensors, a central control unit sends a synchronous trigger signal to all sensors. This synchronous signal is transmitted through wired or wireless communication protocols, such as using Ethernet or Bluetooth technology, to ensure that all sensors start data collection at the same time. The collection parameter settings include sampling frequency and sampling duration. The sampling frequency is set according to the Nyquist theorem, which is twice the expected highest frequency of vibration. The sampling duration covers the entire vibration operation period. During data collection, the signal quality is monitored in real time. If the signal is found to be saturated or the noise is too high, the sensor gain or filter settings are dynamically adjusted. The collected vibration signals are stored in the data recorder in the form of time series, which are used for feature extraction and analysis later.
[0061] S2, from the vibration signal, respectively extracting the first type of frequency domain features related to construction operation and the second type of frequency domain features related to environmental disturbance, which is specifically implemented as:
[0062] When extracting features from vibration signals, the collected vibration signals are first preprocessed to eliminate noise interference. The preprocessing method includes using a digital filter to smooth the original signal. The digital filter uses a low-pass filter design, and the cutoff frequency is set according to the typical frequency range of construction vibration and environmental disturbance. For example, the cutoff frequency can be set to a value higher than a certain percentage of the working frequency of the vibration equipment to ensure that the main frequency components of construction vibration are retained while high-frequency noise is suppressed. The preprocessed vibration signal is used for subsequent frequency domain feature extraction to ensure that the data quality meets the analysis requirements.
[0063] When performing spectral analysis on the vibration signal to obtain the frequency component distribution of the vibration signal, a frequency domain transformation method is used to convert the time domain vibration signal into a frequency domain representation. The frequency domain transformation method is based on the principle of discrete Fourier transform and obtains the frequency component distribution by calculating the amplitude of the vibration signal at different frequency points. In specific implementation, the vibration signal is divided into multiple time windows, and the length of each time window is set according to the duration of the vibration signal, for example, the length of the time window can cover several vibration periods to ensure that the frequency resolution is high enough. Then, the frequency domain transformation is applied to the signal in each time window to obtain the corresponding frequency amplitude spectrum, which represents the energy distribution of the vibration signal at each frequency component. The output of the frequency component distribution includes the frequency axis and the corresponding amplitude data, which are used for subsequent feature calculation.
[0064] When calculating the main band energy based on the frequency component distribution as the first type of frequency domain feature, the main band in the frequency component distribution is first identified, which is defined as the continuous frequency range with the highest amplitude. The identification method is achieved by finding the peak area in the frequency amplitude spectrum, for example, by setting an amplitude threshold, the frequency interval with an amplitude exceeding the threshold is determined as the main band, and the amplitude threshold is dynamically adjusted according to the overall amplitude level of the frequency component distribution, for example, the threshold can be set to a certain multiple of the average amplitude of the frequency amplitude spectrum. The main band energy is calculated by integrating the square of the amplitude of all frequency points in the main band to obtain the result, which is the main band energy value, representing the concentration of vibration energy caused by construction operations.
[0065] When detecting whether there is a transient pulse component in the vibration signal as the second type of frequency domain feature, a time domain analysis method is used to identify the sudden high amplitude event in the vibration signal, and the detection process includes calculating the envelope of the vibration signal, which is obtained by Hilbert transform or peak detection method, and then analyzing the trend of the envelope to locate the transient pulse. In specific implementation, a pulse detection threshold is set, which is calculated based on the standard deviation or average amplitude of the vibration signal, for example, the threshold can be set to the average amplitude of the signal plus several times the standard deviation, when the envelope exceeds the threshold, it is determined that there is a transient pulse component. The characteristic parameters of the transient pulse component include pulse amplitude, duration and occurrence time, which are used to distinguish the vibration caused by environmental disturbance.
[0066] When calculating the main band energy, it is necessary to ensure that the definition of the main band is consistent with the characteristics of the construction equipment, for example, the working frequency range of the vibrating equipment can be obtained by experimental measurement or equipment manual, and the range of the main band is set within the working frequency range to avoid mistakenly including other frequency components in the calculation. After the main band energy is calculated, normalization processing is performed to eliminate the influence of signal amplitude variation, and the normalization method divides the main band energy by the total energy of the entire frequency component distribution to obtain a relative energy value, which is used as the final output of the first type of frequency domain feature.
[0067] When detecting the transient pulse component, considering the diversity of environmental disturbances, the pulse detection threshold needs to be dynamically adjusted according to historical data or field conditions, for example, by analyzing the statistical characteristics of the vibration signal without construction operations to set a baseline threshold, and updating it in real time during the construction process according to the signal changes. The detected transient pulse component is further analyzed for its time characteristics, such as calculating the rise time and fall time of the pulse, to distinguish different types of disturbance sources, such as vehicle load or mechanical impact.
[0068] During spectral analysis, the accuracy of the frequency component distribution depends on the sampling frequency and signal length. The sampling frequency is set based on the highest frequency component of the vibration signal; for example, the sampling frequency should be at least twice the expected highest frequency to satisfy the Nyquist sampling theorem. The signal length is selected according to the analysis requirements, such as covering the entire vibration operation cycle to ensure completeness. The calculation results of the frequency component distribution are stored in the form of a spectrum or data list for easy feature extraction and comparison later.
[0069] The calculation of the dominant frequency band energy involves the selection of the frequency bandwidth. A bandwidth that is too narrow may miss important frequency components, while a bandwidth that is too wide may introduce noise. Therefore, the bandwidth is set based on the typical frequency width of construction vibration. For example, the main frequency range of vibration is determined through preliminary experiments, and the dominant frequency band is set within this range. During energy calculation, a numerical integration method is used to sum the frequency amplitude spectrum within the dominant frequency band. The integration result is expressed in energy units, such as square millimeters per square second, to maintain dimensional consistency.
[0070] The detection of transient pulse components must consider the influence of signal baseline drift. In the preprocessing stage, high-pass filtering or baseline correction methods are used to eliminate low-frequency drift, ensuring the accuracy of pulse detection. The detection algorithm uses a sliding window approach to analyze the vibration signal in real time. The window size is set according to the expected duration of the pulse; for example, the window length covers several sampling points to capture short-term pulse events. When the signal within the window exceeds a threshold, the pulse characteristics are recorded.
[0071] The entire feature extraction process ensures the independence of the first and second types of frequency domain features. The first type of frequency domain features focuses on the steady-state vibration characteristics of construction operations, while the second type targets transient events caused by environmental disturbances. These two types of features are processed in parallel to improve recognition efficiency. After extraction, the feature data is used for subsequent vibration signal classification and stability assessment, forming a complete quality control chain.
[0072] S3. Based on the energy distribution of the first type of frequency domain characteristics and the transient characteristics of the second type of frequency domain characteristics, the vibration signal is identified as construction vibration or environmental vibration. Specifically, this is implemented as follows:
[0073] In the process of vibration signal identification, the first type of frequency domain features and the second type of frequency domain features extracted from the vibration signal are first obtained. The first type of frequency domain features include the proportion of the main frequency band energy, which is calculated by dividing the main frequency band energy by the total energy of the vibration signal. The total energy is obtained by summing the squares of the amplitudes of all frequency points in the entire frequency component distribution. The proportion of the main frequency band energy characterizes the relative magnitude of the construction vibration energy in the total vibration energy. The second type of frequency domain features includes the amplitude of the transient pulse component, which is obtained by detecting sudden pulse events in the vibration signal and recording their maximum amplitude. A pulse event is defined as the process in which the signal amplitude rises sharply and recovers rapidly in a short period of time.
[0074] When determining whether the main frequency band energy proportion in the first type of frequency domain feature exceeds the first threshold value, the setting of the first threshold value is based on statistical analysis of historical construction vibration data. Specifically, the main frequency band energy proportion values of multiple construction vibration samples are collected, the statistical distribution characteristics of these values are calculated, such as calculating the mean and standard deviation, and the first threshold value is set to the mean plus several times the standard deviation to ensure that it can cover the energy concentration degree of typical construction vibration while avoiding misjudgment of environmental vibration. The first threshold value is expressed in percentage form, for example, the value range is between 60% and 80%, which is adjusted according to specific construction conditions and equipment types.
[0075] When determining whether the transient pulse component amplitude in the second type of frequency domain feature is lower than the second threshold value, the setting of the second threshold value is based on feature analysis of environmental disturbance vibration data. Pure environmental disturbance samples such as vibration signals caused by traffic load or wind-induced vibration are collected, the transient pulse component amplitudes in these samples are extracted, and their statistical characteristics are calculated, such as taking the maximum value or upper percentile of the amplitudes as a reference. The second threshold value is set to this reference value minus a certain margin to distinguish between construction vibration and environmental disturbance. The unit of the second threshold value is consistent with the amplitude unit of the vibration signal, for example, millimeters per second squared, and the specific value is dynamically adjusted according to the on-site environmental noise level.
[0076] When the main frequency band energy proportion exceeds the first threshold value and the transient pulse component amplitude is lower than the second threshold value, the vibration signal is identified as construction vibration. This judgment logic is based on the characteristics of construction vibration, which usually shows energy concentration in a specific frequency band and lack of sudden pulses. The identification process is realized through comparison circuits or program condition statements, such as using logical AND operation to check whether both conditions are met, and outputting the classification label of construction vibration when the conditions are met.
[0077] When the main frequency band energy proportion does not exceed the first threshold value or the transient pulse component amplitude does not fall below the second threshold value, the vibration signal is identified as environmental vibration. This identification considers that environmental disturbances often have random pulses and dispersed energy distribution. The identification result is stored in the form of binary flag or text description for subsequent quality control steps.
[0078] When setting the first threshold value, the influence of construction equipment type and concrete state needs to be considered. For example, different vibrating equipment has different working frequency ranges, resulting in differences in the reference value of the main frequency band energy proportion. Therefore, the first threshold value is determined through on-site calibration test. During calibration, vibration signals under normal construction conditions are recorded, the typical range of the main frequency band energy proportion is calculated, and the first threshold value is set near the lower limit of this range to ensure sensitivity.
[0079] When setting the second threshold, the spatio-temporal variation of environmental disturbance needs to be considered. For example, during a busy traffic period, the amplitude of environmental vibration is higher, and the second threshold is correspondingly increased. During a quiet night period, the second threshold is decreased. The dynamic adjustment of the second threshold is realized by real-time monitoring of the environmental vibration background value. For example, the vibration signal without construction is collected every certain time to update the statistical benchmark of the amplitude of the transient pulse component.
[0080] In the judgment process, the calculation of the energy proportion of the main frequency band needs to ensure the accuracy of the total energy. The total energy calculation includes all frequency components to avoid missing low or high frequency components. The frequency component distribution comes from the frequency spectrum analysis result. The parameters of the frequency spectrum analysis, such as frequency resolution and window function type, need to be consistent to ensure the comparability of the characteristics.
[0081] The detection of the amplitude of the transient pulse component needs to exclude the influence of signal noise. High-frequency noise is eliminated by pre-processing filtering. The pulse amplitude is the maximum value of the detected pulse event. The amplitude unit is consistent with the original vibration signal, such as acceleration unit, to maintain dimensional consistency.
[0082] The implementation of the recognition logic uses an iterative method, which processes each vibration signal segment independently. The length of the signal segment is set according to the vibration duration, such as covering several seconds of data in each segment, to ensure the timeliness of the recognition. The recognition result is associated with the time stamp, which is convenient for tracing the vibration event.
[0083] The entire recognition process ensures the cooperative use of the first type of frequency domain feature and the second type of frequency domain feature. The first type of frequency domain feature reflects the steady-state characteristics of vibration, and the second type of frequency domain feature reflects the transient characteristics of vibration. The dual judgment improves the classification accuracy. The construction vibration or environmental vibration label output by the recognition is used for subsequent steps such as structure coupling degree evaluation, forming a coherent quality control process.
[0084] S4, when the vibration signal is recognized as construction vibration, the energy values of the vibration signals at multiple preset positions are calculated and the attenuation gradient of the energy values with the increase of the distance between the preset positions and the vibration source is analyzed. The structure coupling degree is evaluated according to the steepness of the attenuation gradient, which is specifically implemented as:
[0085] When the vibration signal is recognized as construction vibration, first, the data segments recognized as construction vibration are selected from the collected vibration signals. These data segments correspond to the vibration signals at multiple preset positions. The preset positions are distributed on the wet joint structure, including the new and old concrete joint interface and the action area of the vibrating equipment. The vibration source position is determined through construction drawings or obtained by on-site measurement, such as the insertion point of the vibrating rod, whose coordinates are measured and recorded by a total station or GPS device. The distance from the preset position to the vibration source is obtained by calculating the straight-line distance between the two points. The distance unit uses meters or millimeters to maintain dimensional consistency.
[0086] When calculating the energy value of each preset position vibration signal, the square integral method is used to process the vibration signal. The square integral is realized by squaring the point-by-point data of the time series of the vibration signal and then numerically integrating the squared sequence in the time domain. The numerical integration uses the trapezoidal rule or the rectangular rule to accumulate the squared signal values within the sampling time interval to obtain the energy value. The unit of the energy value is related to the physical quantity of the vibration signal. For example, if the vibration signal is an acceleration signal, the energy unit can be square millimeters per square second. When calculating, the sampling frequency must be consistent to avoid errors.
[0087] When establishing the relationship model between the energy value and the distance from the corresponding preset position to the vibration source, the linear fitting method is used to correlate the energy value and the distance. Linear fitting is realized by the least squares principle, which finds a straight line that minimizes the sum of the squared deviations of the data points of the energy value and the distance from the straight line. The fitting process includes calculating the slope and intercept of the straight line. The slope represents the attenuation gradient, and the intercept represents the theoretical energy value at the vibration source. The relationship model is represented in the form of a linear equation, for example, the energy value equals the slope multiplied by the distance plus the intercept.
[0088] When determining the attenuation gradient of the energy value with respect to the distance through the relationship model, the attenuation gradient is directly obtained from the slope of the linear fitting. A negative slope value indicates that the energy decreases with increasing distance. The larger the absolute value of the slope, the faster the attenuation. The unit of the attenuation gradient is related to the units of energy and distance. For example, if the energy unit is square millimeters per square second and the distance unit is meters, the attenuation gradient unit is square millimeters per square second per meter.
[0089] When evaluating the structural coupling degree based on the linearity of the attenuation gradient, the linearity is evaluated by calculating the residual sum of squares of the data points and the fitted straight line in the relationship model. The residual sum of squares is the sum of the squares of the difference between the actual energy value of each data point and the predicted value of the fitted straight line. The smaller the residual sum of squares, the better the data points fit the straight line, the higher the linearity, and the better the structural coupling degree. The structural coupling degree represents the bonding quality of the new and old concrete interfaces of the wet joint. High linearity indicates that the vibration energy propagates uniformly in the structure and the interface is well coupled.
[0090] When calculating the energy value, the square integral needs to consider the duration of the vibration signal to ensure that the integral covers the entire construction vibration event, such as from the start to the end of the vibration, to avoid partial signal loss. After calculating the energy value, normalization processing is performed, such as dividing by the signal length, to eliminate the influence of time length and facilitate comparison of different positions.
[0091] In the linear fitting process, the distance data needs to be accurately measured, and the distance value is calculated from the preset position coordinates and the vibration source coordinates. The coordinate measurement uses engineering measuring instruments, with an accuracy of millimeters, to ensure the reliability of the distance data. If there are abnormal points, such as abnormally high or low energy values, during fitting, they can be removed through statistical methods such as the Rida criterion to improve the accuracy of fitting.
[0092] The calculation of the residual sum of squares includes calculating the residual for each data point of the preset position, which is the actual energy value minus the predicted value of the fitted straight line. Then, the sum of the squares of all residuals is calculated, and the value of the residual sum of squares is expressed in square units of energy, such as square millimeters per fourth power of seconds. When calculating, it is necessary to ensure that the number of data points is sufficient, such as at least three preset positions, to ensure statistical significance.
[0093] The evaluation of structural coupling degree is based on the size of the residual sum of squares, which is judged by setting a reference threshold. The threshold is determined through historical data or experiments, such as calculating the distribution of the residual sum of squares in known well-coupled wet joint samples, and taking the upper quantile of the distribution as the threshold. When the residual sum of squares is lower than the threshold, it is determined that the structural coupling degree is high.
[0094] The steepness of the decay gradient is reflected by the absolute value of the slope. A large absolute value of the slope indicates fast decay, which may be due to structural defects or interface discontinuity. However, in combination with the linearity evaluation, if the linearity is high, even if the decay is fast, it may also be well coupled, so a comprehensive judgment is needed.
[0095] When calculating the energy value, the vibration signal needs to be preprocessed, such as filtering to remove high-frequency noise, to ensure that the energy value accurately reflects the construction vibration characteristics. The preprocessing parameters are consistent with those in the feature extraction stage.
[0096] The applicability of linear fitting is verified by the correlation coefficient, which calculates the linear correlation between the energy value and the distance. A value close to 1 indicates a strong linear relationship, but in this embodiment, the residual sum of squares is mainly used as a linearity indicator.
[0097] The threshold setting of the residual sum of squares takes into account the construction conditions and material characteristics, such as concrete strength or vibration equipment type. The threshold is adjusted through on-site calibration tests, which are conducted on wet joints with known quality states, and the corresponding relationship between the residual sum of squares and the coupling degree is recorded.
[0098] Throughout the process, the data flow goes from the vibration signal to the energy value, then to the relationship model and the decay gradient, and finally outputs the structural coupling degree evaluation result, ensuring logical closed-loop between steps, consistent terminology, and all calculations based on physical dimensions correct and feasible.
[0099] S5, based on the structural coupling degree, analyze the concentration of the first type of frequency domain features to evaluate the stability of construction vibration, and adjust the control parameters of the vibrating device according to the stability evaluation result, which is specifically implemented as:
[0100] When the stability evaluation index is lower than the set threshold value, the vibration intensity control parameter of the vibrating device is increased, the set threshold value is determined by experimental calibration, when the experimental calibration, the value of the stability evaluation index is recorded under normal construction conditions, the statistical distribution of these values is calculated, for example, the lower percentile of the distribution is taken as the threshold value, the threshold value is represented by a numerical value, for example, fifty, when the stability evaluation index is lower than the threshold value, it is determined that the construction vibration is unstable, and the vibration intensity control parameter needs to be increased, the vibration intensity control parameter includes the vibration frequency or the vibration amplitude, the increase of the vibration intensity control parameter is realized by adjusting the driving voltage or current of the vibrating device, for example, the vibration frequency is increased by a certain percentage or the vibration amplitude is increased by a certain value, the adjustment amount is set based on the difference between the stability evaluation index and the threshold value, the larger the difference, the larger the adjustment amount, to ensure that the vibration effect is improved.
[0101] When the stability evaluation index is lower than the set threshold value, the vibration intensity control parameter of the vibrating device is increased, the set threshold value is determined by experimental calibration, when the experimental calibration, the value of the stability evaluation index is recorded under normal construction conditions, the statistical distribution of these values is calculated, for example, the lower percentile of the distribution is taken as the threshold value, the threshold value is represented by a numerical value, for example, fifty, when the stability evaluation index is lower than the threshold value, it is determined that the construction vibration is unstable, and the vibration intensity control parameter needs to be increased, the vibration intensity control parameter includes the vibration frequency or the vibration amplitude, the increase of the vibration intensity control parameter is realized by adjusting the driving voltage or current of the vibrating device, for example, the vibration frequency is increased by a certain percentage or the vibration amplitude is increased by a certain value, the adjustment amount is set based on the difference between the stability evaluation index and the threshold value, the larger the difference, the larger the adjustment amount, to ensure that the vibration effect is improved.
[0102] When the stability evaluation index is lower than the set threshold value, the vibration intensity control parameter of the vibrating device is increased, the set threshold value is determined by experimental calibration, when the experimental calibration, the value of the stability evaluation index is recorded under normal construction conditions, the statistical distribution of these values is calculated, for example, the lower percentile of the distribution is taken as the threshold value, the threshold value is represented by a numerical value, for example, fifty, when the stability evaluation index is lower than the threshold value, it is determined that the construction vibration is unstable, and the vibration intensity control parameter needs to be increased, the vibration intensity control parameter includes the vibration frequency or the vibration amplitude, the increase of the vibration intensity control parameter is realized by adjusting the driving voltage or current of the vibrating device, for example, the vibration frequency is increased by a certain percentage or the vibration amplitude is increased by a certain value, the adjustment amount is set based on the difference between the stability evaluation index and the threshold value, the larger the difference, the larger the adjustment amount, to ensure that the vibration effect is improved.
[0103] The vibration intensity control parameter of the vibrating device is maintained or reduced when the stability evaluation index is not lower than the set threshold, which means that the construction vibration is in a stable state, maintaining the vibration intensity control parameter means maintaining the current vibration frequency or vibration amplitude unchanged, reducing the vibration intensity control parameter means reducing the vibration frequency or vibration amplitude, and the reduction degree is determined based on the magnitude by which the stability evaluation index is higher than the threshold, the greater the magnitude, the more the reduction, so as to optimize energy consumption and avoid excessive vibration.
[0104] In calculating the dispersion degree of the main frequency band energy distribution, the main frequency band energy values of multiple vibration signal segments are used, the segments are segmented from the continuous vibration signal, and the segment length is set according to the vibration duration, for example, each segment covers several seconds of data to ensure statistical reliability. After the dispersion degree is calculated, normalization processing is performed, for example, dividing by the average value of the main frequency band energy to eliminate the influence of amplitude.
[0105] The establishment of the stability evaluation index needs to ensure the dimensional consistency of the structure coupling degree and the dispersion degree. The structure coupling degree is a dimensionless ratio, and the dispersion degree is in the unit of energy. Therefore, the dispersion degree is inverted and normalized before weighting, so that the two are in a comparable scale.
[0106] The dynamic adjustment of the set threshold considers the construction stage and material changes. For example, in the initial setting stage of the concrete, the threshold is set higher to strictly require stability, and in the final setting stage, the threshold is appropriately reduced. The adjustment is realized by real-time monitoring of construction parameters such as concrete slump or temperature.
[0107] The adjustment of the vibration intensity control parameter is realized through closed-loop control. The adjustment signal is sent to the controller of the vibrating device, and the controller changes the driving output according to the parameter value, for example, increases the voltage to increase the vibration frequency. The adjustment step is set according to the response characteristics of the device to avoid excessive fluctuations.
[0108] During the whole process, the calculation of the stability evaluation index and the threshold comparison are carried out in an iterative manner, and each construction cycle is updated once to ensure real-time performance. The adjustment result is fed back to the vibration operation to form a quality control cycle.
[0109] S6, execute the adjusted control parameter to control the vibrating device to vibrate the wet joint, and the specific implementation is:
[0110] When the adjusted control parameters are executed, the adjusted vibrating intensity control parameters are first obtained from the construction vibration stability evaluation results. The vibrating intensity control parameters include vibrating frequency and vibrating amplitude. These parameters are determined based on the comparison results of the stability evaluation index and the set threshold value. When the stability evaluation index is lower than the set threshold value, the vibrating intensity control parameters are correspondingly increased. When the stability evaluation index is not lower than the set threshold value, the vibrating intensity control parameters are maintained or reduced. The adjusted vibrating intensity control parameters exist in the form of digital signals or analog signals, for example, the vibrating frequency is in units of hertz, and the vibrating amplitude is in units of millimeters.
[0111] When the adjusted vibrating intensity control parameters are output to the controller of the vibrating device, industrial communication protocols are used for data transmission. The industrial communication protocols include Modbus or PROFIBUS. The transmission mode is wired or wireless communication. Shielded twisted pair wires are used for wired communication, and radio frequency transmission is used for wireless communication. The vibrating intensity control parameters are encoded before data transmission. The encoding methods include pulse code modulation or amplitude modulation, which ensures signal anti-interference in the transmission process. The controller is located inside the vibrating device or in an independent cabinet outside the vibrating device. The controller type is a programmable logic controller or an embedded controller, which has digital input and output interfaces and analog input and output interfaces. The vibrating intensity control parameters are transmitted to the controller through these interfaces.
[0112] When the vibrating device is driven by the controller to perform vibrating operation on the wet joint with the adjusted vibrating intensity control parameters, the controller receives the vibrating intensity control parameters, generates a driving signal through internal program processing, and transmits the driving signal to the actuator of the vibrating device. The driving signal includes a voltage signal or a current signal. The voltage signal ranges from zero to ten volts, and the current signal ranges from four to twenty milliamperes. The actuator includes an electric vibrator or a hydraulic vibrator. The electric vibrator adjusts the vibrating frequency by changing the input voltage, and the hydraulic vibrator controls the vibrating amplitude by adjusting the hydraulic oil flow. The actuator adjusts the vibrating intensity according to the driving signal, so that the vibrating device performs vibrating operation on the wet joint with the adjusted vibrating intensity control parameters.
[0113] When the vibrating intensity control parameters are output, the reasonableness of the parameter values needs to be verified. The verification methods include range check and logic check. The range check confirms that the vibrating frequency and the vibrating amplitude are within the allowable working range of the device, for example, the vibrating frequency ranges from fifty hertz to one hundred hertz, and the vibrating amplitude ranges from one millimeter to five millimeters. The logic check confirms that the combination of the vibrating frequency and the vibrating amplitude meets the device operation specification, avoiding overload or failure. After verification, the parameters are sent to the controller.
[0114] The controller drives the vibrating device in a closed-loop control mode to maintain the vibrating intensity stable. The closed-loop control monitors the actual output of the vibrating device in real time through sensors, for example, uses an encoder to measure the vibrating frequency and uses a displacement sensor to measure the vibrating amplitude. The actual output value is compared with the adjusted vibrating intensity control parameter. The driving signal is adjusted according to the deviation to ensure that the actual vibrating intensity is consistent with the set value.
[0115] During the vibrating operation, the controller records the execution parameters and the running state, including the vibrating start time, the duration and the vibrating intensity change curve. These data are stored in the controller memory or the external storage device for quality tracing and performance analysis.
[0116] The whole execution process ensures that the vibrating intensity control parameter is accurately implemented. The parameter transmission and the driving link set a fault protection mechanism. For example, when the communication is interrupted or the device is abnormal, the controller automatically switches to the safety mode. The safety mode uses the default vibrating intensity control parameter or stops the vibrating operation to prevent construction quality defects.
[0117] When the vibrating device vibrates the wet joint, the operation path and the coverage range are determined according to the construction scheme. The vibrating device moves along the length direction of the wet joint. The moving speed is coordinated and controlled with the vibrating intensity to ensure that the concrete is uniformly compacted. After the vibrating operation is completed, the device automatically stops and feeds back the completion signal to the control system.
[0118] Embodiment 2 Figure 2 A structural diagram of a UHPC wet joint construction quality intelligent control system is given. The UHPC wet joint construction quality intelligent control system comprises the following modules.
[0119] A signal acquisition module is used to acquire the vibration signals of multiple preset positions of the wet joint structure.
[0120] A feature extraction module is used to extract a first type of frequency domain feature related to the construction operation and a second type of frequency domain feature related to the environmental disturbance from the vibration signals.
[0121] A vibration identification module is used to identify the vibration signals as construction vibration or environmental vibration based on the energy distribution of the first type of frequency domain feature and the transient characteristics of the second type of frequency domain feature.
[0122] A coupling evaluation module is used to calculate the energy values of the vibration signals of multiple preset positions and analyze the decay gradient of the energy values with the increase of the distance between the preset positions and the vibration source when the vibration signals are identified as construction vibration. The structure coupling degree is evaluated according to the steepness of the decay gradient.
[0123] The parameter adjustment module is configured to analyze concentration of the first type of frequency domain features based on structural coupling degree to evaluate stability of the construction vibration, and adjust the control parameter of the vibrating device according to the stability evaluation result.
[0124] The operation execution module is configured to execute the adjusted control parameter to control the vibrating device to perform the vibrating operation on the wet joint.
[0125] The calculation in the embodiments is all the dimensionless numerical calculation, and the preset parameters and the threshold values in the calculation are set by the person skilled in the art according to the actual situation.
[0126] The above embodiments can be realized wholly or partially by software, hardware, firmware or any combination thereof. When realized by software, the above embodiments can be realized in the form of a computer program product wholly or partially.
[0127] Those skilled in the art can appreciate that the modules and algorithm steps of the examples described in conjunction with the embodiments disclosed herein can be realized by electronic hardware or a combination of computer software and electronic hardware. Whether the functions are realized by hardware or software depends on the specific application and the constraints of the technical solution. The person skilled in the art can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.
[0128] In addition, the functional modules in each of the embodiments of the present application can be integrated in one processing module, or each module can exist physically independently, or two or more modules can be integrated in one module.
[0129] In the several embodiments provided by the present application, it should be understood that the disclosed system, device and method can be implemented in other ways. For example, the device embodiments described above are only schematic. The division of the modules is only a logical function division. There can be another division manner in actual implementation, for example, a plurality of modules or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the modules shown or discussed can be indirect coupling or communication connection through some interface, device or module, and can be electrical, mechanical or other forms.
[0130] The above describes only the specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art can easily think of changes or replacements within the technical range disclosed by the present application, which should be covered in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.
[0131] Finally: the above only for the preferred embodiments of the present application, and not for limiting the present application, any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application, should be included in the scope of protection of the present application.
Claims
1. A method for intelligent control of construction quality of UHPC wet joints, characterized in that, Includes the following steps: S1. Collect vibration signals at multiple preset locations on the wet joint structure; S2. Extract the first type of frequency domain features related to construction operations and the second type of frequency domain features related to environmental disturbances from the vibration signal, respectively. S3. Based on the energy distribution of the first type of frequency domain characteristics and the transient characteristics of the second type of frequency domain characteristics, the vibration signal is identified as construction vibration or environmental vibration. S4. When the vibration signal is identified as construction vibration, calculate the energy value of the vibration signal at multiple preset locations and analyze the attenuation gradient of the energy value as the distance between the preset location and the vibration source increases. Evaluate the structural coupling degree based on the steepness of the attenuation gradient. S5. Based on the structural coupling degree analysis, the concentration of the first type of frequency domain characteristics is used to evaluate the stability of construction vibration, and the control parameters of the vibrating equipment are adjusted according to the stability evaluation results. S6. Execute the adjusted control parameters to control the vibrating equipment to vibrate the wet joint.
2. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, Vibration signals were collected at multiple preset locations on the wet joint structure, including: Vibration sensors are pre-positioned at the interface between the old and new concrete in the wet joint structure, and vibration sensors are also pre-positioned at the action area of the vibrating equipment. Vibration signals are collected synchronously through the vibration sensors.
3. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, From the vibration signal, the first type of frequency domain features related to construction operations and the second type of frequency domain features related to environmental disturbances are extracted, including: Spectral analysis is performed on the vibration signal to obtain the frequency component distribution of the vibration signal; The energy of the main frequency band is calculated based on the frequency component distribution as the first type of frequency domain feature; Simultaneously, the presence of transient pulse components in the vibration signal is detected as a second type of frequency domain feature.
4. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, Based on the energy distribution of the first type of frequency domain characteristics and the transient characteristics of the second type of frequency domain characteristics, vibration signals are identified as construction vibrations or environmental vibrations, including: Determine whether the proportion of the main frequency band energy in the first type of frequency domain features exceeds the first threshold and whether the amplitude of the transient pulse component in the second type of frequency domain features is lower than the second threshold. When the energy proportion of the main frequency band exceeds the first threshold and the amplitude of the transient pulse component is lower than the second threshold, the vibration signal is identified as construction vibration; otherwise, the vibration signal is identified as environmental vibration.
5. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, When vibration signals are identified as construction vibrations, the energy values of vibration signals at multiple preset locations are calculated, and the attenuation gradient of the energy value as the distance between the preset location and the vibration source increases is analyzed. The structural coupling degree is assessed based on the steepness of the attenuation gradient, including: Calculate the energy value of the vibration signal at each preset position and establish a relationship model between the energy value and the distance from the corresponding preset position to the vibration source. The attenuation gradient of energy value as a function of distance is determined using a relational model; Evaluate structural coupling based on the linearity of the decay gradient; A higher linearity indicates a better structural coupling.
6. The intelligent control method for construction quality of UHPC wet joints according to claim 5, characterized in that, The calculation of the energy value of the vibration signal at each preset position and the establishment of the relationship model between the energy value and the distance from the corresponding preset position to the vibration source include: calculating the energy value by performing square integral on the vibration signal at each preset position, and then linearly fitting the energy value with the distance from the corresponding preset position to the vibration source to establish the relationship model.
7. The intelligent control method for construction quality of UHPC wet joints according to claim 5, characterized in that, The linearity assessment of structural coupling based on decay gradients includes calculating the sum of squared residuals between the data points in the relational model and the fitted line to assess linearity; the smaller the sum of squared residuals, the higher the structural coupling.
8. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, The stability of construction vibration is assessed by analyzing the concentration of the first type of frequency domain characteristics based on structural coupling, and the control parameters of the vibratory compaction equipment are adjusted according to the stability assessment results, including: A stability evaluation index is established based on the structural coupling degree and the dispersion of the main frequency band energy distribution in the first type of frequency domain characteristics; When the stability assessment index is lower than the set threshold, increase the vibration intensity control parameter of the vibratory equipment. When the stability assessment index is not lower than the set threshold, maintain or reduce the vibration intensity control parameters of the vibratory equipment.
9. The intelligent control method for construction quality of UHPC wet joints according to claim 1, characterized in that, The adjusted control parameters are executed to control the vibratory equipment to vibrate the wet joints, including: The adjusted vibration intensity control parameters are output to the controller of the vibration equipment; The controller drives the vibratory equipment to perform vibration operation on the wet joints using the adjusted vibration intensity control parameters.
10. A UHPC wet joint construction quality intelligent control system, used to implement the UHPC wet joint construction quality intelligent control method according to any one of claims 1-9, characterized in that, Includes the following modules: The signal acquisition module is used to acquire vibration signals at multiple preset locations on the wet joint structure; The feature extraction module is used to extract, from the vibration signal, the first type of frequency domain features related to construction operations and the second type of frequency domain features related to environmental disturbances, respectively. The vibration identification module is used to identify vibration signals as construction vibration or environmental vibration based on the energy distribution of the first type of frequency domain characteristics and the transient characteristics of the second type of frequency domain characteristics. The coupling evaluation module is used to calculate the energy values of vibration signals at multiple preset locations and analyze the attenuation gradient of the energy values as the distance between the preset location and the vibration source increases when the vibration signal is identified as construction vibration. The structural coupling degree is evaluated based on the steepness of the attenuation gradient. The parameter adjustment module is used to evaluate the stability of construction vibration by analyzing the concentration of the first type of frequency domain characteristics based on the structural coupling degree, and to adjust the control parameters of the vibrating equipment according to the stability evaluation results. The operation execution module is used to execute the adjusted control parameters to control the vibrating equipment to vibrate the wet joint.
Citation Information
Patent Citations
New and old bridge splicing seam construction quality control method based on vibration coordination recognition
CN116818088A
Control model obtaining method and system of bridge prefabricated part reinforced concrete vibration scheme optimization and vibration quality real-time evaluation method based on energy analysis, controller and control method
CN119047170A
Stratum disturbance analysis method and system under shield construction coupling effect
CN120408222A
Big data analysis tunneling support quality evaluation method
CN120493109A
Geological disaster networking monitoring and early warning method
CN120564356A