CN118094067B - Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer - Google Patents
Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer Download PDFInfo
- Publication number
- CN118094067B CN118094067B CN202410521484.8A CN202410521484A CN118094067B CN 118094067 B CN118094067 B CN 118094067B CN 202410521484 A CN202410521484 A CN 202410521484A CN 118094067 B CN118094067 B CN 118094067B
- Authority
- CN
- China
- Prior art keywords
- water
- stable
- signal
- subbase
- acceleration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005056 compaction Methods 0.000 title claims abstract description 49
- 238000004364 calculation method Methods 0.000 title claims description 22
- 230000001133 acceleration Effects 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims abstract description 59
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000009467 reduction Effects 0.000 claims abstract description 6
- 230000005284 excitation Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 11
- 230000010354 integration Effects 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 16
- 208000002740 Muscle Rigidity Diseases 0.000 claims 7
- 239000011800 void material Substances 0.000 claims 4
- 230000005279 excitation period Effects 0.000 claims 2
- 238000005259 measurement Methods 0.000 claims 2
- 239000010410 layer Substances 0.000 abstract description 50
- 239000002356 single layer Substances 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000010276 construction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000013441 quality evaluation Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Road Paving Machines (AREA)
Abstract
本发明公开了水稳底基层智能压实综合刚度系数指标计算方法，包括如下步骤：采集振动轮质心加速度信号与压路机位置信号，加速度信号降噪处理，二次积分法求解水稳底基层位移信号，利用振动轮‑水稳底基层脱空段加速度信号拟合求解激振力信号，利用激振力信号与振动轮‑水稳底基层接触段加速度信号求解振动轮‑水稳底基层接触力信号，基于余弦函数拟合求解接触力信号幅值、相位、平均值，将接触力信号代入三自由度多层弹塑性动力学模型理论位移方程中，拟合求解水稳底基层各弹塑性力学参数，计算水稳底基层综合刚度系数指标。本发明克服了现有智能压实刚度指标对水稳底基层采用弹性假定与单层假定。
The present invention discloses a method for calculating a comprehensive stiffness coefficient index of a water-stable base layer for intelligent compaction, comprising the following steps: collecting a vibration wheel mass center acceleration signal and a roller position signal, performing noise reduction processing on the acceleration signal, solving a displacement signal of the water-stable base layer by a quadratic integral method, solving an exciting force signal by fitting an acceleration signal of a debonding section of the vibration wheel-water-stable base layer, solving a contact force signal of the vibration wheel-water-stable base layer by using an exciting force signal and an acceleration signal of a contact section of the vibration wheel-water-stable base layer, solving an amplitude, a phase, and an average value of the contact force signal by fitting based on a cosine function, substituting the contact force signal into a theoretical displacement equation of a three-degree-of-freedom multilayer elastic-plastic dynamic model, solving various elastic-plastic mechanical parameters of the water-stable base layer by fitting, and calculating a comprehensive stiffness coefficient index of the water-stable base layer. The present invention overcomes the problem that the existing intelligent compaction stiffness index adopts elastic assumptions and single-layer assumptions for the water-stable base layer.
Description
技术领域Technical Field
本发明属于水稳底基层智能压实检测技术领域，更具体地说，涉及一种水稳底基层智能压实综合刚度系数指标计算方法。The present invention belongs to the technical field of intelligent compaction detection of water-stable base layers, and more specifically, relates to a method for calculating a comprehensive stiffness coefficient index of intelligent compaction of a water-stable base layer.
背景技术Background technique
压实是水稳底基层施工的重要环节，充分压实能提高水稳底基层的强度、稳定性，对于保障路基路面结构的稳定性、耐久性具有重要意义。当前水稳底基层施工主要采用灌砂法、环刀法等传统压实度测量方法，通过随机抽样检测评价压实质量，存在随机性、滞后性、破坏性的不足。针对上述问题，研究人员提出了智能压实技术，通过在压路机上安装加速度传感器、定位装置，实时采集振动轮加速度数据与压路机位置信息，计算智能压实控制指标并与位置信息对应，实现了全部工作区域压实质量的实时无损检测。当前最常用的水稳底基层智能压实指标为刚度系数指标，其计算方法为取单个激振周期内振动轮/水稳底基层接触力-水稳底基层位移滞回曲线中静平衡点与最大位移点的割线斜率或速度为零点的切线斜率。然而，刚度系数指标对水稳底基层采用弹性假定与单层假定，未消除水稳底基层塑性性质与下层路基刚度特性对水稳底基层动力学响应的影响，导致刚度系数无法反映水稳底基层自身的弹塑性刚度特性。因此该指标仅适用于压实完成阶段单层路基的压实质量评价，对水稳底基层压实质量评价的准确性较低。此外，该指标取割线斜率或切线斜率的计算方法仅能利用个别数据点，导致指标量值易受压路机工作功率不稳定、振动轮机械噪声等现场偶然因素干扰，在实践中稳定性较低。针对上述问题，本技术领域需要一种考虑水稳底基层弹塑性性质与下层路基刚度特性的水稳底基层智能压实刚度系数指标，同时在实践中具有较高稳定性，以实现水稳底基层压实质量的准确表征。Compaction is an important part of the construction of water-stable subbase. Sufficient compaction can improve the strength and stability of the water-stable subbase, which is of great significance for ensuring the stability and durability of the roadbed and pavement structure. At present, the construction of water-stable subbase mainly adopts traditional compaction measurement methods such as sand filling method and ring knife method. The compaction quality is evaluated by random sampling detection, which has the shortcomings of randomness, hysteresis and destructiveness. In response to the above problems, researchers have proposed intelligent compaction technology. By installing acceleration sensors and positioning devices on the roller, the acceleration data of the vibrating wheel and the position information of the roller are collected in real time, and the intelligent compaction control index is calculated and corresponds to the position information, realizing the real-time non-destructive detection of the compaction quality of the entire working area. The most commonly used intelligent compaction index of the water-stable subbase is the stiffness coefficient index, which is calculated by taking the secant slope of the static equilibrium point and the maximum displacement point in the hysteresis curve of the vibrating wheel/water-stable subbase contact force-water-stable subbase displacement in a single excitation cycle or the tangent slope of the zero velocity point. However, the stiffness coefficient index adopts elastic assumption and single-layer assumption for the water-stable subbase, and does not eliminate the influence of the plastic properties of the water-stable subbase and the stiffness characteristics of the lower roadbed on the dynamic response of the water-stable subbase, resulting in the stiffness coefficient being unable to reflect the elastic-plastic stiffness characteristics of the water-stable subbase itself. Therefore, this index is only applicable to the compaction quality evaluation of a single-layer roadbed at the compaction completion stage, and the accuracy of the compaction quality evaluation of the water-stable subbase is low. In addition, the calculation method of the secant slope or tangent slope of this index can only use individual data points, resulting in the index value being easily disturbed by accidental factors on site such as unstable working power of the roller and mechanical noise of the vibrating wheel, and having low stability in practice. In response to the above problems, the technical field needs a water-stable subbase intelligent compaction stiffness coefficient index that takes into account the elastic-plastic properties of the water-stable subbase and the stiffness characteristics of the lower roadbed, and at the same time has high stability in practice, so as to achieve accurate characterization of the compaction quality of the water-stable subbase.
因此，亟需一种新的考虑水稳底基层弹塑性性质与下层路基刚度特性的水稳底基层智能压实综合刚度系数指标计算与现场应用方法。Therefore, a new method for calculating and applying the comprehensive stiffness coefficient index of the water-stable base layer by intelligent compaction that takes into account the elastic-plastic properties of the water-stable base layer and the stiffness characteristics of the underlying roadbed is urgently needed.
发明内容Summary of the invention
本发明提出了一种水稳底基层智能压实综合刚度系数指标计算与现场应用方法。The invention proposes a method for calculating and applying a comprehensive stiffness coefficient index of a water-stable base layer through intelligent compaction.
为了解决上述技术问题至少之一，根据本发明的一方面，提供了一种水稳底基层智能压实综合刚度系数指标计算与现场应用方法，包括如下步骤：In order to solve at least one of the above technical problems, according to one aspect of the present invention, a method for calculating and applying a comprehensive stiffness coefficient index of a water-stable subbase layer by intelligent compaction is provided, comprising the following steps:
步骤1：在振动轮中心轴上安装加速度传感器，用于采集振动压实过程中的振动轮加速度信号；Step 1: Install an acceleration sensor on the central axis of the vibration wheel to collect the acceleration signal of the vibration wheel during the vibration compaction process;
步骤2：在压路机驾驶舱顶部安装GPS定位装置，用于采集压路机位置信息，与智能压实刚度指标对应反映压实质量分布状况；Step 2: Install a GPS positioning device on the top of the roller cab to collect the roller position information, which corresponds to the intelligent compaction stiffness index to reflect the compaction quality distribution;
步骤3：去除原始加速度信号中的高频分量，对原始加速度信号进行降噪处理；Step 3: Remove the high-frequency components in the original acceleration signal and perform noise reduction on the original acceleration signal;
步骤4：采用平均值法、二次积分法将加速度信号转换为位移信号；Step 4: Use the average value method and quadratic integration method to convert the acceleration signal into a displacement signal;
步骤5：将振动轮-水稳底基层脱空段加速度信号代入水稳底基层振动压实动力学模型脱空段理论控制方程中，求解振动轮激振力信号；Step 5: Substitute the acceleration signal of the vibration wheel-water-stable base layer voiding section into the theoretical control equation of the voiding section of the water-stable base layer vibration compaction dynamics model to solve the vibration wheel exciting force signal;
步骤6：将激振力信号与振动轮-水稳底基层脱空段加速度信号代入水稳底基层振动压实动力学模型接触段理论控制方程中，求解振动轮-水稳底基层接触力信号；Step 6: Substitute the exciting force signal and the acceleration signal of the de-airing section of the vibration wheel-water-stable base into the contact section theoretical control equation of the vibration compaction dynamics model of the water-stable base to solve the contact force signal of the vibration wheel-water-stable base;
步骤7：利用余弦函数拟合振动轮-水稳底基层接触力信号，求解接触力信号的幅值、相位、平均值，拟合公式如下：Step 7: Use the cosine function to fit the contact force signal of the vibration wheel-water-stable base layer, and solve the amplitude, phase, and average value of the contact force signal. The fitting formula is as follows:
（5） (5)
式（5）中 为接触力信号，为振动轮-水稳底基层接触力信号平均值， 为振动轮-水稳底基层接触力信号幅值， 为振动轮-水稳底基层接触力信号初始相位，为激振频率；根据上式拟合振动轮-水稳底基层接触力信号，求解 、 与 ；In formula (5) is the contact force signal, is the average value of the contact force signal between the vibration wheel and the water-stable base layer, is the contact force signal amplitude between vibration wheel and water-stable base layer, is the initial phase of the contact force signal between the vibration wheel and the water-stable base layer, is the excitation frequency; according to the above formula, the contact force signal between the vibration wheel and the water-stable base is fitted to solve , and ;
步骤8：将接触力信号代入基于三自由度多层弹塑性动力学模型的水稳底基层理论位移方程中，拟合水稳底基层实测位移信号，求解水稳底基层各弹塑性力学参数；Step 8: Substitute the contact force signal into the theoretical displacement equation of the water-stable base layer based on the three-degree-of-freedom multi-layer elastic-plastic dynamic model, fit the measured displacement signal of the water-stable base layer, and solve the elastic-plastic mechanical parameters of the water-stable base layer;
步骤9：通过提取振动轮-水稳底基层接触力幅值与水稳底基层位移幅值之比，建立水稳底基层综合刚度系数kcs与各弹塑性力学参数的数学关系式，即的计算公式，如下所示：Step 9: By extracting the ratio of the vibration wheel-water-stable base contact force amplitude to the water-stable base displacement amplitude, a mathematical relationship between the comprehensive stiffness coefficient kcs of the water-stable base and various elastic-plastic mechanical parameters is established, that is, The calculation formula is as follows:
（7） (7)
式（7）中 为水稳底基层弹性刚度系数， 为水稳底基层弹性粘度系数， 为水稳底基层塑性刚度系数， 为激振频率。In formula (7) is the elastic stiffness coefficient of the water-stable subbase, is the elastic viscosity coefficient of the water-stable base layer, is the plastic stiffness coefficient of the water-stable base layer, is the excitation frequency.
进一步的，步骤1中，加速度传感器的采样频率不低于1000Hz，量程不低于100m/s2。Furthermore, in step 1, the sampling frequency of the acceleration sensor is not less than 1000 Hz, and the measuring range is not less than 100 m/s2.
进一步的，步骤2中，GPS定位装置采用RTK实时动态差分技术，在驾驶舱顶部安装移动站，在路边视野开阔处安装基站，利用基站定位信息实时修正移动站定位信息。Furthermore, in step 2, the GPS positioning device uses RTK real-time dynamic differential technology, installs a mobile station on the top of the cockpit, installs a base station at a roadside with a wide field of view, and uses the base station positioning information to correct the mobile station positioning information in real time.
进一步的，步骤3中，加速度信号的降噪处理先采用快速傅里叶变换方法将加速度时域信号转换成频域信号，去除频率高于140Hz的高频分量，再采用快速傅里叶逆变换方法将降噪后的频域信号转换为时域信号。Furthermore, in step 3, the acceleration signal is denoised by first converting the acceleration time domain signal into a frequency domain signal using the fast Fourier transform method, removing high-frequency components with a frequency higher than 140 Hz, and then converting the denoised frequency domain signal into a time domain signal using the inverse fast Fourier transform method.
进一步的，步骤4中，平均值法用于修正加速度信号，基于计算区间内加速度信号均值为0的原则，选取1s计算区间，去除加速度信号的直流分量，计算公式如下：Furthermore, in step 4, the average value method is used to correct the acceleration signal. Based on the principle that the average value of the acceleration signal in the calculation interval is 0, a 1s calculation interval is selected to remove the DC component of the acceleration signal. The calculation formula is as follows:
（1） (1)
式（1）中表示加速度直流分量，表示加速度信号值，表示计算区间内采样点的数目，表示去除直流分量后的加速度信号值。In formula (1), represents the DC component of acceleration, Indicates the acceleration signal value, Indicates the number of sampling points in the calculation interval, Indicates the acceleration signal value after removing the DC component.
二次积分法用于获取速度信号与位移信号，积分后基于单个激振周期内积分面积为0的原则去除直流分量，再将各周期内的速度信号与位移信号分别按顺序拼接为完整信号；直流分量的计算公式如下：The quadratic integration method is used to obtain the velocity signal and displacement signal. After integration, the DC component is removed based on the principle that the integral area in a single excitation cycle is 0, and then the velocity signal and displacement signal in each cycle are spliced into a complete signal in sequence; the calculation formula of the DC component is as follows:
（2） (2)
式（2）中 为直流分量， 为激振频率， 为采样频率， 为采样点数，为速度或位移信号值， 为截取采样点后单个激振周期剩余区段的比例长度。In formula (2), is the DC component, is the excitation frequency, is the sampling frequency, is the number of sampling points, is the velocity or displacement signal value, It is the proportional length of the remaining section of a single excitation cycle after the sampling point is intercepted.
进一步的，步骤5中，求解振动轮激振力信号的求解公式如下：Furthermore, in step 5, the solution formula for solving the vibration wheel exciting force signal is as follows:
（3） (3)
式（3）中 为激振力幅值， 为激振频率， 为激振力初始相位， 为振动轮质量， 为机架质量， 为加速度信号， 为重力加速度；选用振动轮-水稳底基层脱空段加速度波形，根据上式拟合激振力信号，求解 、 与 。In formula (3) is the exciting force amplitude, is the excitation frequency, is the initial phase of the exciting force, is the mass of the vibration wheel, is the rack mass, is the acceleration signal, is the acceleration of gravity; the acceleration waveform of the vibration wheel-water-stable subgrade debonding section is selected, and the exciting force signal is fitted according to the above formula to solve , and .
进一步的，步骤6中，求解振动轮-水稳底基层接触力信号的求解公式如下：Furthermore, in step 6, the solution formula for the contact force signal between the vibration wheel and the water-stable base layer is as follows:
（4）。 (4).
进一步的，步骤8中，水稳底基层理论位移方程如下：Furthermore, in step 8, the theoretical displacement equation of the water-stable subbase is as follows:
（6） (6)
式（6）中 为水稳底基层位移幅值， 为位移相对于接触力的相位滞后， 为水稳底基层弹性刚度系数， 为水稳底基层弹性粘度系数， 为水稳底基层塑性刚度系数， 为使水稳底基层发生塑性变形的振动轮-水稳底基层接触力阈值， 为振动轮影响深度范围内的下层路基弹性刚度系数；根据上式拟合水稳底基层实测位移信号，求解 、 、 、 。In formula (6) is the displacement amplitude of the water-stable base layer, is the phase lag of displacement relative to contact force, is the elastic stiffness coefficient of the water-stable subbase, is the elastic viscosity coefficient of the water-stable base layer, is the plastic stiffness coefficient of the water-stable base layer, The contact force threshold between the vibration wheel and the water-stable base layer for plastic deformation of the water-stable base layer is: is the elastic stiffness coefficient of the lower roadbed within the depth range affected by the vibration wheel; according to the above formula, the measured displacement signal of the water-stable subgrade is fitted to solve , , , .
根据本发明的另一方面，提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明的水稳底基层智能压实综合刚度系数指标计算与现场应用方法中的步骤。According to another aspect of the present invention, a computer-readable storage medium is provided, on which a computer program is stored, and when the program is executed by a processor, the steps of the method for calculating and applying the comprehensive stiffness coefficient index of the water-stable base layer intelligently compacted according to the present invention are implemented.
根据本发明的又一方面，提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本发明的水稳底基层智能压实综合刚度系数指标计算与现场应用方法中的步骤。According to another aspect of the present invention, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the method for calculating and applying the comprehensive stiffness coefficient index of the intelligent compaction of the water-stable base layer of the present invention are implemented.
与现有的技术相比较，本发明的上述方法的有益效果为：Compared with the prior art, the above method of the present invention has the following beneficial effects:
（1）本发明提出了考虑水稳底基层弹塑性性质与下层路基刚度特性的综合刚度系数指标，克服了现有智能压实刚度指标对水稳底基层采用弹性假定与单层假定，导致计算得到的水稳底基层刚度特性与实际刚度特性偏差较大，表征压实质量的准确性较低的不足；(1) The present invention proposes a comprehensive stiffness coefficient index that takes into account the elastic-plastic properties of the water-stable subbase and the stiffness characteristics of the lower roadbed, which overcomes the shortcomings of the existing intelligent compaction stiffness index that adopts elastic assumptions and single-layer assumptions for the water-stable subbase, resulting in a large deviation between the calculated stiffness characteristics of the water-stable subbase and the actual stiffness characteristics, and a low accuracy in characterizing the compaction quality;
（2）本发明改进了现有刚度系数指标计算方法，对单个激振周期内的有效数据进行最小二乘法拟合，提高了综合刚度系数指标的稳定性，克服了指标仅计算个别数据点的切线或割线斜率，导致指标量值易受压路机工作功率不稳定、振动轮机械噪声等现场偶然因素干扰的不足；(2) The present invention improves the existing stiffness coefficient index calculation method, performs least squares fitting on the effective data within a single excitation cycle, improves the stability of the comprehensive stiffness coefficient index, and overcomes the shortcomings of the index that only calculates the tangent or secant slope of individual data points, resulting in the index value being easily disturbed by accidental factors on site such as unstable working power of the roller and mechanical noise of the vibrating wheel;
（3）本发明提出了振动轮-水稳底基层接触力信号获取方法，打破了当前压路机设备缺乏振动轮相位监测技术，导致振动轮-水稳底基层接触力信号无法测量的技术壁垒。(3) The present invention proposes a method for obtaining the contact force signal of a vibrating wheel-water-stable base layer, which breaks the technical barrier that the current roller equipment lacks a vibrating wheel phase monitoring technology, resulting in the inability to measure the contact force signal of the vibrating wheel-water-stable base layer.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本发明实施例的技术方案，下面将对实施例的附图作简单地介绍，显而易见地，下面描述中的附图仅仅涉及本发明的一些实施例，而非对本发明的限制。In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings of the embodiments will be briefly introduced below. Obviously, the drawings in the following description only relate to some embodiments of the present invention, but are not intended to limit the present invention.
图1是本发明计算模型示意图。FIG1 is a schematic diagram of a calculation model of the present invention.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。To make the purpose, technical solution and advantages of the embodiment of the present invention clearer, the technical solution of the embodiment of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiment of the present invention. Obviously, the described embodiment is only a part of the embodiment of the present invention, not all of the embodiments.
除非另作定义，此处使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。Unless otherwise defined, technical or scientific terms used herein shall have the same meaning as commonly understood by one of ordinary skill in the art to which the invention belongs.
实施例1：如图1所示，本发明提出了一种水稳底基层智能压实综合刚度系数指标计算与现场应用方法，具体步骤如下：Embodiment 1: As shown in FIG1 , the present invention proposes a method for calculating and applying a comprehensive stiffness coefficient index of a water-stable base layer through intelligent compaction, and the specific steps are as follows:
步骤1：在振动轮中心轴上安装加速度传感器，用于采集振动压实过程中的振动轮加速度信号。加速度传感器的采样频率不低于1000Hz，量程不低于100m/s2。Step 1: Install an acceleration sensor on the central axis of the vibration wheel to collect the acceleration signal of the vibration wheel during vibration compaction. The sampling frequency of the acceleration sensor shall not be less than 1000Hz, and the measuring range shall not be less than 100m/s2.
步骤2：在压路机驾驶舱顶部安装GPS定位装置，用于采集压路机位置信息，与智能压实刚度指标对应反映压实质量分布状况。Step 2: Install a GPS positioning device on the top of the roller cab to collect the roller position information, which corresponds to the intelligent compaction stiffness index to reflect the compaction quality distribution status.
GPS定位装置采用RTK实时动态差分技术，在驾驶舱顶部安装移动站，在路边视野开阔处安装基站，利用基站定位信息实时修正移动站定位信息。The GPS positioning device uses RTK real-time dynamic differential technology. A mobile station is installed on the top of the cockpit, and a base station is installed on the roadside with a wide field of view. The base station positioning information is used to correct the mobile station positioning information in real time.
步骤3：去除原始加速度信号中的高频分量，对原始加速度信号进行降噪处理。Step 3: Remove the high-frequency components in the original acceleration signal and perform noise reduction on the original acceleration signal.
加速度信号的降噪处理先采用快速傅里叶变换方法将加速度时域信号转换成频域信号，去除频率高于140Hz的高频分量，再采用快速傅里叶逆变换方法将降噪后的频域信号转换为时域信号。The denoising process of the acceleration signal first converts the acceleration time domain signal into a frequency domain signal using the fast Fourier transform method, removes the high-frequency components with a frequency higher than 140 Hz, and then converts the denoised frequency domain signal into a time domain signal using the inverse fast Fourier transform method.
步骤4：采用平均值法、二次积分法将加速度信号转换为位移信号。平均值法用于修正加速度信号，基于计算区间内加速度信号均值为0的原则，选取1s计算区间，去除加速度信号的直流分量，计算公式如下：Step 4: Use the average method and quadratic integration method to convert the acceleration signal into a displacement signal. The average method is used to correct the acceleration signal. Based on the principle that the average value of the acceleration signal in the calculation interval is 0, a 1s calculation interval is selected to remove the DC component of the acceleration signal. The calculation formula is as follows:
（1） (1)
式（1）中表示加速度直流分量，表示加速度信号值，表示计算区间内采样点的数目，表示去除直流分量后的加速度信号值。In formula (1), represents the DC component of acceleration, Indicates the acceleration signal value, Indicates the number of sampling points in the calculation interval, Indicates the acceleration signal value after removing the DC component.
二次积分法用于获取速度信号与位移信号，积分后基于单个激振周期内积分面积为0的原则去除直流分量，再将各周期内的速度信号与位移信号分别按顺序拼接为完整信号。直流分量的计算公式如下：The quadratic integration method is used to obtain the velocity signal and displacement signal. After integration, the DC component is removed based on the principle that the integral area in a single excitation cycle is 0, and then the velocity signal and displacement signal in each cycle are spliced into a complete signal in sequence. The calculation formula of the DC component is as follows:
（2） (2)
式（2）中 为直流分量， 为激振频率， 为采样频率， 为采样点数， 为速度或位移信号值， 为截取采样点后单个激振周期剩余区段的比例长度。In formula (2), is the DC component, is the excitation frequency, is the sampling frequency, is the number of sampling points, is the velocity or displacement signal value, It is the proportional length of the remaining section of a single excitation cycle after the sampling point is intercepted.
步骤5：将振动轮-水稳底基层脱空段加速度信号代入水稳底基层振动压实动力学模型脱空段理论控制方程中，求解振动轮激振力信号，求解公式如下：Step 5: Substitute the acceleration signal of the vibration wheel-water-stable base debonding section into the theoretical control equation of the debonding section of the water-stable base vibration compaction dynamics model to solve the vibration wheel exciting force signal. The solution formula is as follows:
（3） (3)
式（3）中 为激振力幅值， 为激振频率， 为激振力初始相位， 为振动轮质量， 为机架质量， 为加速度信号， 为重力加速度。选用振动轮-水稳底基层脱空段加速度波形，根据上式拟合激振力信号，求解 、 与 。In formula (3) is the exciting force amplitude, is the excitation frequency, is the initial phase of the exciting force, is the mass of the vibration wheel, is the rack mass, is the acceleration signal, is the acceleration of gravity. The acceleration waveform of the vibration wheel-water-stable subbase debonding section is selected, and the exciting force signal is fitted according to the above formula to solve , and .
步骤6：将激振力信号与振动轮-水稳底基层脱空段加速度信号代入水稳底基层振动压实动力学模型接触段理论控制方程中，求解振动轮-水稳底基层接触力信号，求解公式如下：Step 6: Substitute the exciting force signal and the acceleration signal of the vibration wheel-water-stable base debonding section into the contact section theoretical control equation of the water-stable base vibration compaction dynamics model to solve the vibration wheel-water-stable base contact force signal. The solution formula is as follows:
（4）。 (4).
步骤7：利用余弦函数拟合振动轮-水稳底基层接触力信号，求解接触力信号的幅值、相位、平均值，拟合公式如下：Step 7: Use the cosine function to fit the contact force signal of the vibration wheel-water-stable base layer, and solve the amplitude, phase, and average value of the contact force signal. The fitting formula is as follows:
（5） (5)
式（5）中为接触力信号， 为振动轮-水稳底基层接触力平均值， 为振动轮-水稳底基层接触力幅值， 为振动轮-水稳底基层接触力初始相位，为激振频率；根据上式拟合振动轮-水稳底基层接触力信号，求解 、 与 。In formula (5) is the contact force signal, is the average contact force between the vibration wheel and the water-stable base layer, is the contact force amplitude between vibration wheel and water-stable subgrade, is the initial phase of the contact force between the vibration wheel and the water-stable subgrade, is the excitation frequency; according to the above formula, the contact force signal between the vibration wheel and the water-stable base is fitted to solve , and .
步骤8：将接触力信号代入基于三自由度多层弹塑性动力学模型的水稳底基层理论位移方程中，拟合水稳底基层实测位移信号，求解水稳底基层各弹塑性力学参数，水稳底基层理论位移方程如下：Step 8: Substitute the contact force signal into the theoretical displacement equation of the water-stable base layer based on the three-degree-of-freedom multi-layer elastic-plastic dynamic model, fit the measured displacement signal of the water-stable base layer, and solve the elastic-plastic mechanical parameters of the water-stable base layer. The theoretical displacement equation of the water-stable base layer is as follows:
（6） (6)
式（6）中 为水稳底基层位移幅值， 为位移相对于接触力的相位滞后， 为水稳底基层弹性刚度系数， 为水稳底基层弹性粘度系数， 为水稳底基层塑性刚度系数， 为使水稳底基层发生塑性变形的振动轮-水稳底基层接触力阈值， 为振动轮影响深度范围内的下层路基弹性刚度系数。根据上式拟合水稳底基层实测位移信号，求解、 、 、 。In formula (6) is the displacement amplitude of the water-stable base layer, is the phase lag of displacement relative to contact force, is the elastic stiffness coefficient of the water-stable subbase, is the elastic viscosity coefficient of the water-stable base layer, is the plastic stiffness coefficient of the water-stable base layer, is the contact force threshold between the vibrating wheel and the water-stable subbase that causes plastic deformation of the water-stable subbase, and is the elastic stiffness coefficient of the lower roadbed within the depth range of the vibrating wheel. According to the above formula, the measured displacement signal of the water-stable subbase is fitted to solve , , , .
步骤9：通过提取振动轮-水稳底基层接触力幅值与水稳底基层位移幅值之比，建立水稳底基层综合刚度系数kcs与各弹塑性力学参数的数学关系式，即kcs的计算公式，如下所示：Step 9: By extracting the ratio of the vibration wheel-water-stable base contact force amplitude to the water-stable base displacement amplitude, a mathematical relationship between the comprehensive stiffness coefficient kcs of the water-stable base and various elastic-plastic mechanical parameters is established, that is, the calculation formula of kcs, as shown below:
（7）。 (7).
本发明依托沿江高速项目，在现场进行综合刚度系数的准确性验证。选用压路机型号为徐工XS265。施工层位为水稳碎石底基层，松铺厚度30cm。选取两个压实条带，每个条带选取四个灌砂法测点，每一遍碾压完成后依次选取一个测点测量压实度，用于建立各智能压实控制指标（CMV、刚度系数、综合刚度系数）与压实度间相关关系。The present invention relies on the Yanjiang Expressway project to conduct comprehensive stiffness coefficient on site. The accuracy of the test was verified. The roller model selected was XCMG XS265. The construction layer was a water-stable crushed stone subbase with a loose thickness of 30 cm. Two compaction strips were selected, and four sand filling method measuring points were selected for each strip. After each rolling, one measuring point was selected in turn to measure the compaction degree, which was used to establish the correlation between each intelligent compaction control index (CMV, stiffness coefficient, comprehensive stiffness coefficient) and the compaction degree.
综合刚度系数与压实度的相关性（R2=0.88）显著优于CMV、刚度系数（R2=0.68、R2=0.08），表明本发明提出的综合刚度系数指标在表征水稳底基层压实质量方面具有更高准确性与稳定性。Comprehensive stiffness coefficient The correlation with compaction degree (R2=0.88) is significantly better than CMV and stiffness coefficient (R2=0.68, R2=0.08), indicating that the comprehensive stiffness coefficient proposed in this invention is The index has higher accuracy and stability in characterizing the compaction quality of water-stable base.
实施例2：本实施例的计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现实施例1的水稳底基层智能压实综合刚度系数指标计算与现场应用方法中的步骤。Embodiment 2: The computer-readable storage medium of this embodiment stores a computer program thereon, which, when executed by a processor, implements the steps in the method for calculating and applying the comprehensive stiffness coefficient index of the water-stable base layer intelligently compacted in embodiment 1.
本实施例的计算机可读存储介质可以是终端的内部存储单元，例如终端的硬盘或内存；本实施例的计算机可读存储介质也可以是所述终端的外部存储设备，例如终端上配备的插接式硬盘，智能存储卡，安全数字卡，闪存卡等；进一步地，计算机可读存储介质还可以既包括终端的内部存储单元也包括外部存储设备。The computer-readable storage medium of this embodiment may be an internal storage unit of the terminal, such as a hard disk or memory of the terminal; the computer-readable storage medium of this embodiment may also be an external storage device of the terminal, such as a plug-in hard disk, a smart memory card, a secure digital card, a flash memory card, etc. equipped on the terminal; further, the computer-readable storage medium may also include both an internal storage unit of the terminal and an external storage device.
本实施例的计算机可读存储介质用于存储计算机程序以及终端所需的其他程序和数据，计算机可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。The computer-readable storage medium of this embodiment is used to store computer programs and other programs and data required by the terminal. The computer-readable storage medium can also be used to temporarily store data that has been output or is to be output.
实施例3：本实施例的计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现实施例1的水稳底基层智能压实综合刚度系数指标计算与现场应用方法中的步骤。Example 3: The computer device of this example includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for calculating and applying the comprehensive stiffness coefficient index of the water-stable base layer intelligently compacted in Example 1 are implemented.
本实施例中，处理器可以是中央处理单元，还可以是其他通用处理器、数字信号处理器、专用集成电路、现成可编程门阵列或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等，通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等；存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据，存储器的一部分还可以包括非易失性随机存取存储器，例如，存储器还可以存储设备类型的信息。In this embodiment, the processor may be a central processing unit, or other general-purpose processors, digital signal processors, application-specific integrated circuits, readily available programmable gate arrays or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The memory may include a read-only memory and a random access memory, and provide instructions and data to the processor. A portion of the memory may also include a non-volatile random access memory. For example, the memory may also store information about the device type.
本领域内的技术人员应明白，实施例公开的内容可提供为方法、系统、或计算机程序产品。因此，本方案可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且，本方案可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器和光学存储器等）上实施的计算机程序产品的形式。Those skilled in the art will appreciate that the disclosed contents of the embodiments may be provided as methods, systems, or computer program products. Therefore, the present solution may take the form of hardware embodiments, software embodiments, or embodiments combining software and hardware. Moreover, the present solution may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program codes.
本方案是参照根据本方案实施例的方法、和计算机程序产品的流程图和/或方框图来描述的，应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合；可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present scheme is described with reference to the method according to the embodiment of the present scheme, and the flowchart and/or block diagram of the computer program product. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the processes and/or boxes in the flowchart and/or block diagram, can be implemented by computer program instructions; these computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate a device for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-OnlyMemory，ROM）或随机存储记忆体（RandomAccessMemory，RAM）等。Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.
本发明所述实例仅仅是对本发明的优选实施方式进行描述，并非对本发明构思和范围进行限定，在不脱离本发明设计思想的前提下，本领域工程技术人员对本发明的技术方案作出的各种变形和改进，均应落入本发明的保护范围。The examples described in the present invention are merely descriptions of the preferred implementation modes of the present invention, and are not intended to limit the concept and scope of the present invention. Without departing from the design concept of the present invention, various modifications and improvements made to the technical solutions of the present invention by engineers and technicians in this field should all fall within the protection scope of the present invention.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410521484.8A CN118094067B (en) | 2024-04-28 | 2024-04-28 | Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410521484.8A CN118094067B (en) | 2024-04-28 | 2024-04-28 | Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118094067A CN118094067A (en) | 2024-05-28 |
| CN118094067B true CN118094067B (en) | 2024-07-05 |
Family
ID=91148070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410521484.8A Active CN118094067B (en) | 2024-04-28 | 2024-04-28 | Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118094067B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111122087A (en) * | 2020-01-06 | 2020-05-08 | 山东大学 | A system and method for measuring stiffness coefficient and viscous damping coefficient of compacted soil |
| CN114021066A (en) * | 2021-10-14 | 2022-02-08 | 东南大学 | Measurement and Calculation Method of Low Variability Modulus Smart Compaction Metrics |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114117615B (en) * | 2021-12-02 | 2022-09-20 | 中国科学院武汉岩土力学研究所 | Method and device for determining performance of roadbed of highway section and processing equipment |
| CN114839970A (en) * | 2022-04-11 | 2022-08-02 | 中铁三局集团第五工程有限公司 | Intelligent paving method and system for asphalt water-stable base layer and asphalt lower surface layer |
-
2024
- 2024-04-28 CN CN202410521484.8A patent/CN118094067B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111122087A (en) * | 2020-01-06 | 2020-05-08 | 山东大学 | A system and method for measuring stiffness coefficient and viscous damping coefficient of compacted soil |
| CN114021066A (en) * | 2021-10-14 | 2022-02-08 | 东南大学 | Measurement and Calculation Method of Low Variability Modulus Smart Compaction Metrics |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118094067A (en) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Feng et al. | A kNN algorithm for locating and quantifying stiffness loss in a bridge from the forced vibration due to a truck crossing at low speed | |
| CN109858156B (en) | Vehicle and structure information simultaneous identification method based on axle coupling vibration | |
| CN104120644B (en) | A kind of detection method of road-surface evenness based on Gravity accelerometer | |
| WO2020119305A1 (en) | Real-time detection algorithm for vibration information of vibration roller based on acceleration magnitude | |
| CN102506986B (en) | Test system and method for mode and vibration of self-supporting tower and large-span power transmission tower | |
| CN110532657B (en) | Pier structure state evaluation method based on variable speed vehicle excitation and wavelet packet analysis | |
| Quqa et al. | Automatic identification of dense damage-sensitive features in civil infrastructure using sparse sensor networks | |
| CN104408917B (en) | Based on smart mobile phone acceleration transducer real-time traffic method of estimation | |
| CN104008644B (en) | A kind of traffic noise on urban roads measuring method based on Gradient Descent | |
| CN105588599B (en) | The adaptive corrective method of vehicle-mounted mobile measuring system vibration displacement error | |
| CN110631786A (en) | A rapid evaluation method for the bearing capacity of beam bridges based on parking vibration response | |
| CN104215421A (en) | Quick bridge impact coefficient determination method | |
| CN104750830B (en) | The cycle method for digging of time series data | |
| CN118094067B (en) | Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer | |
| CN108860011A (en) | A kind of overload of vehicle recognition methods and system | |
| Zhang et al. | The state-of-the-art on time-frequency signal processing techniques for high-resolution representation of nonlinear systems in engineering | |
| Yang et al. | Bridge damage identification based on synchronous statistical moment theory of vehicle–bridge interaction | |
| CN105915594A (en) | Device for monitoring earth rockfill and asphalt concrete dam material rigidity in real time in rolling processes | |
| CN111353238A (en) | Pier scouring depth identification method based on vehicle sensing | |
| Erdogan et al. | Friction coefficient measurement for autonomous winter road maintenance | |
| CN110826807A (en) | Method for rapidly predicting dynamic resilience modulus of roadbed filler in seasonal frozen region | |
| CN119332746A (en) | A method and system for on-site detection of backfill foundation compaction quality | |
| JP6809172B2 (en) | Estimating method, estimation program and estimation device | |
| Bruwer et al. | Comparison of GPS and MEMS support for smartphone-based driver behavior monitoring | |
| CN114004041B (en) | Measurement and calculation method of intelligent compaction vibration modulus index to avoid phase detection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |