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
- stable
- water
- subbase
- signal
- 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 52
- 238000004364 calculation method Methods 0.000 title claims abstract description 30
- 230000001133 acceleration Effects 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 53
- 238000006073 displacement reaction Methods 0.000 claims abstract description 34
- 239000011800 void material Substances 0.000 claims abstract description 13
- 230000009467 reduction Effects 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 50
- 208000002740 Muscle Rigidity Diseases 0.000 claims description 34
- 238000003860 storage Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 12
- 230000010354 integration Effects 0.000 claims description 12
- 230000005279 excitation period Effects 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 4
- 239000002356 single layer Substances 0.000 abstract description 4
- 239000010410 layer Substances 0.000 description 46
- 238000010586 diagram Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 3
- 238000001514 detection method Methods 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
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000004519 manufacturing process 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
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000003068 static effect Effects 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 invention discloses a calculation method of an intelligent compaction comprehensive rigidity coefficient index of a water-stable subbase, which comprises the following steps: collecting a vibration wheel centroid acceleration signal and a road roller position signal, conducting acceleration signal noise reduction processing, solving a water-stable subbase displacement signal by a quadratic integral method, solving an exciting force signal by utilizing a vibration wheel-water-stable subbase void section acceleration signal fitting, solving a vibration wheel-water-stable subbase contact force signal by utilizing the exciting force signal and the vibration wheel-water-stable subbase contact section acceleration signal, solving a contact force signal amplitude value, a phase value and an average value based on cosine function fitting, substituting the contact force signal into a three-degree-of-freedom multilayer elastoplastic dynamic model theoretical displacement equation, solving each elastoplastic mechanical parameter of the water-stable subbase by fitting, and calculating a water-stable subbase comprehensive stiffness coefficient index. The invention overcomes the defect that the existing intelligent compaction stiffness index adopts elasticity assumption and single-layer assumption for the water-stable subbase.
Description
Technical Field
The invention belongs to the technical field of intelligent compaction detection of a water-stable subbase layer, and particularly relates to a method for calculating an intelligent compaction comprehensive rigidity coefficient index of the water-stable subbase layer.
Background
Compaction is an important link of construction of the water-stable subbase, and full compaction can improve the strength and stability of the water-stable subbase, so that the method has important significance for guaranteeing the stability and durability of the roadbed pavement structure. The current construction of the water-stable subbase mainly adopts traditional compaction degree measuring methods such as a sand filling method, a ring cutter method and the like, and the compaction quality is evaluated through random sampling detection, so that the defects of randomness, hysteresis and destructiveness exist. According to the intelligent compaction technology, researchers provide an intelligent compaction technology, acceleration data of the vibrating wheel and position information of the road roller are collected in real time through installing an acceleration sensor and a positioning device on the road roller, intelligent compaction control indexes are calculated and correspond to the position information, and real-time nondestructive testing of compaction quality of all working areas is achieved. The most commonly used intelligent compaction index of the water-stable subbase layer at present is a rigidity coefficient index, and the calculation method is to take the secant slope of a static balance point and a maximum displacement point or the tangent slope of which the speed is zero point in a vibration wheel/water-stable subbase layer contact force-water-stable subbase layer displacement hysteresis curve in a single excitation period. However, the stiffness coefficient index adopts elasticity assumption and single-layer assumption on the water-stable subbase layer, so that the influence of the plastic property of the water-stable subbase layer and the stiffness characteristic of the lower roadbed on the dynamic response of the water-stable subbase layer is not eliminated, and the stiffness coefficient cannot reflect the elastoplastic stiffness characteristic of the water-stable subbase layer. Therefore, the index is only suitable for the compaction quality evaluation of the single-layer roadbed at the compaction completion stage, and the accuracy of the compaction quality evaluation of the water-stable subbase layer is low. In addition, the calculation method of the index taking secant slope or tangent slope can only use individual data points, so that the index magnitude is easily interfered by site accidental factors such as unstable working power of the road roller, vibration turbine noise and the like, and the stability is lower in practice. Aiming at the problems, the technical field needs an intelligent compaction stiffness coefficient index of a water-stable subbase considering the elastoplasticity property of the water-stable subbase and the stiffness characteristic of a lower roadbed, and has higher stability in practice so as to realize accurate representation of compaction quality of the water-stable subbase.
Therefore, a new calculation and field application method for the intelligent compaction comprehensive rigidity coefficient index of the water-stable subbase layer considering the elastoplasticity property of the water-stable subbase layer and the rigidity property of the subbase bed is needed.
Disclosure of Invention
The invention provides a water-stable subbase intelligent compaction comprehensive rigidity coefficient index calculation and field application method.
In order to solve at least one of the above technical problems, according to an aspect of the present invention, there is provided a method for calculating and applying a water-stable subbase intelligent compaction comprehensive stiffness coefficient index, including the steps of:
Step 1: an acceleration sensor is arranged on the central shaft of the vibrating wheel and is used for collecting the acceleration signal of the vibrating wheel in the vibrating compaction process;
step 2: a GPS positioning device is arranged at the top of the cockpit of the road roller and used for collecting the position information of the road roller and reflecting the compaction quality distribution condition corresponding to the intelligent compaction stiffness index;
step 3: removing high-frequency components in the original acceleration signals, and carrying out noise reduction treatment on the original acceleration signals;
step 4: converting the acceleration signal into a displacement signal by adopting an average value method and a secondary integration method;
Step 5: substituting the acceleration signal of the shaking wheel-water stable subbase void section into a theoretical control equation of a shaking compaction dynamic model void section of the water stable subbase, and solving a shaking force signal of the shaking wheel;
step 6: substituting the exciting force signal and the acceleration signal of the shaking wheel-water stable subbase void section into a theoretical control equation of a contact section of a shaking compaction dynamic model of the water stable subbase, and solving the shaking wheel-water stable subbase contact force signal;
step 7: fitting a vibration wheel-water stable subbase contact force signal by using a cosine function, and solving the amplitude, phase and average value of the contact force signal, wherein the fitting formula is as follows:
（5）
in (5) In order to contact the force signal,For the mean value of the oscillating wheel-water stable sub-layer contact force signals,For the vibration wheel-water stable sub-layer contact force signal amplitude,For the initial phase of the vibrating wheel-water stable sub-layer contact force signal,Is the excitation frequency; solving according to the fit vibration wheel-water stable subbase contact force signal 、And (3) with ；
Step 8: substituting the contact force signal into a water-stable subbase theoretical displacement equation based on a three-degree-of-freedom multilayer elastoplastic dynamic model, fitting the water-stable subbase actual measurement displacement signal, and solving each elastoplastic mechanical parameter of the water-stable subbase;
step 9: by extracting the ratio of the contact force amplitude of the vibrating wheel-water stable subbase layer to the displacement amplitude of the water stable subbase layer, the mathematical relationship between the comprehensive rigidity coefficient kcs of the water stable subbase layer and each elastoplastic mechanical parameter is established, namely The calculation formula of (2) is as follows:
（7）
in (7) Is the elastic rigidity coefficient of the water-stable subbase layer,Is the elastic viscosity coefficient of the water-stable subbing layer,Is the plastic rigidity coefficient of the water-stable subbase layer,Is the excitation frequency.
Further, in the step 1, the sampling frequency of the acceleration sensor is not lower than 1000Hz, and the measuring range is not lower than 100m/s2.
In step 2, the GPS positioning device adopts RTK real-time dynamic differential technology, a mobile station is installed at the top of the cockpit, a base station is installed at the wide roadside view, and the positioning information of the mobile station is corrected in real time by using the positioning information of the base station.
In step 3, the noise reduction process of the acceleration signal firstly adopts a fast fourier transform method to convert the acceleration time domain signal into a frequency domain signal, removes high frequency components with frequency higher than 140Hz, and then adopts an inverse fast fourier transform method to convert the frequency domain signal after noise reduction into a time domain signal.
Further, in step 4, the average value method is used for correcting the acceleration signal, and based on the principle that the average value of the acceleration signal in the calculation interval is 0, a 1s calculation interval is selected, the direct current component of the acceleration signal is removed, and the calculation formula is as follows:
（1）
In the formula (1) The direct current component of the acceleration is indicated,The value of the acceleration signal is indicated,Representing the number of sampling points within the computation interval,Indicating the acceleration signal value after removal of the dc component.
The secondary integration method is used for obtaining a speed signal and a displacement signal, removing a direct current component based on the principle that the integration area is 0 in a single excitation period after integration, and splicing the speed signal and the displacement signal in each period into complete signals in sequence respectively; the calculation formula of the direct current component is as follows:
（2）
In the formula (2) As a direct current component of the power supply,For the excitation frequency to be the same,For the sampling frequency to be the same,For the number of sample points,In order to be a velocity or displacement signal value,The proportional length of the rest section of the single excitation period after the sampling point is intercepted.
In step 5, a solution formula for solving the exciting force signal of the vibrating wheel is as follows:
（3）
In (3) For the amplitude of the exciting force,For the excitation frequency to be the same,For the initial phase of the exciting force,For the mass of the vibrating wheel,Is the quality of the machine frame,As the acceleration signal, the acceleration signal is,Gravitational acceleration; selecting vibration wheel-water stable subbase void section acceleration waveform, fitting exciting force signal according to the above formula, and solving 、 And (3) with。
Further, in step 6, the solution formula for solving the vibration wheel-water stable subbase contact force signal is as follows:
（4）。
Further, in step 8, the theoretical displacement equation of the water-stable subbase layer is as follows:
（6）
In (6) Is the displacement amplitude of the water-stable subbase layer,In order to shift the phase lag with respect to the contact force,Is the elastic rigidity coefficient of the water-stable subbase layer,Is the elastic viscosity coefficient of the water-stable subbing layer,Is the plastic rigidity coefficient of the water-stable subbase layer,To plastically deform the water-stable underlayment,The elastic rigidity coefficient of the lower roadbed in the depth range is influenced by the vibrating wheel; according to the above-mentioned fitting water-stable subbase actual measurement displacement signal, solving、 、 、 。
According to another aspect of the present invention, there is provided a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of the water-stable under-layer intelligent compaction comprehensive stiffness coefficient index calculation and field application method of the present invention.
According to a further aspect of the invention there is provided a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, said processor implementing the steps in the water-stable under-layer intelligent compaction comprehensive stiffness index calculation and field application method of the invention when said program is executed.
Compared with the prior art, the method has the beneficial effects that:
(1) The invention provides a comprehensive rigidity coefficient index taking the elastoplasticity property of the water-stable subbase layer and the rigidity property of the lower roadbed into consideration, and overcomes the defects that the rigidity property of the water-stable subbase layer obtained by calculation has larger deviation from the actual rigidity property and the accuracy of representing compaction quality is lower because the existing intelligent compaction rigidity index adopts elasticity assumption and single-layer assumption for the water-stable subbase layer;
(2) The invention improves the existing rigidity coefficient index calculation method, carries out least square fitting on the effective data in a single excitation period, improves the stability of the comprehensive rigidity coefficient index, and overcomes the defect that the index magnitude is easily interfered by site accidental factors such as unstable working power of a road roller, vibration turbine noise and the like because the index only calculates the tangential line or secant slope of individual data points;
(3) The invention provides a method for acquiring a vibration wheel-water stability subbase contact force signal, which breaks through the technical barrier that the vibration wheel-water stability subbase contact force signal cannot be measured due to the lack of a vibration wheel phase monitoring technology of the current road roller equipment.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following brief description of the drawings of the embodiments will make it apparent that the drawings in the following description relate only to some embodiments of the present invention and are not limiting of the present invention.
FIG. 1 is a schematic representation of a computational model of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. It will be apparent that the described embodiments are some, but not all, embodiments of the invention.
Unless defined otherwise, technical or scientific terms used herein should be given the ordinary meaning as understood by one of ordinary skill in the art to which this invention belongs.
Example 1: as shown in fig. 1, the invention provides a water-stable subbase intelligent compaction comprehensive rigidity coefficient index calculation and field application method, which comprises the following specific steps:
Step 1: and an acceleration sensor is arranged on the central shaft of the vibrating wheel and is used for collecting the acceleration signal of the vibrating wheel in the vibrating compaction process. The sampling frequency of the acceleration sensor is not lower than 1000Hz, and the measuring range is not lower than 100m/s2.
Step 2: and a GPS positioning device is arranged at the top of the cockpit of the road roller and used for collecting the position information of the road roller and reflecting the compaction quality distribution condition corresponding to the intelligent compaction stiffness index.
The GPS positioning device adopts an RTK real-time dynamic differential technology, a mobile station is installed at the top of a cockpit, a base station is installed at the wide roadside view, and the positioning information of the mobile station is corrected in real time by utilizing the positioning information of the base station.
Step 3: and removing high-frequency components in the original acceleration signal, and carrying out noise reduction treatment on the original acceleration signal.
The noise reduction processing of the acceleration signal firstly adopts a fast Fourier transform method to convert the acceleration time domain signal into a frequency domain signal, removes high-frequency components with the frequency higher than 140Hz, and then adopts an inverse fast Fourier transform method to convert the frequency domain signal after noise reduction into the time domain signal.
Step 4: and converting the acceleration signal into a displacement signal by adopting an average value method and a quadratic integration method. The average value method is used for correcting 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, the direct current component of the acceleration signal is removed, and the calculation formula is as follows:
（1）
In the formula (1) The direct current component of the acceleration is indicated,The value of the acceleration signal is indicated,Representing the number of sampling points within the computation interval,Indicating the acceleration signal value after removal of the dc component.
The secondary integration method is used for obtaining a speed signal and a displacement signal, removing a direct current component based on the principle that the integration area is 0 in a single excitation period after integration, and then splicing the speed signal and the displacement signal in each period into complete signals in sequence. The calculation formula of the direct current component is as follows:
（2）
In the formula (2) As a direct current component of the power supply,For the excitation frequency to be the same,For the sampling frequency to be the same,For the number of sample points,In order to be a velocity or displacement signal value,The proportional length of the rest section of the single excitation period after the sampling point is intercepted.
Step 5: substituting the acceleration signal of the shaking wheel-water stable subbase void section into a theoretical control equation of a shaking compaction dynamic model void section of the water stable subbase, solving a shaking force signal of the shaking wheel, and solving the following formula:
（3）
In (3) For the amplitude of the exciting force,For the excitation frequency to be the same,For the initial phase of the exciting force,For the mass of the vibrating wheel,Is the quality of the machine frame,As the acceleration signal, the acceleration signal is,Gravitational acceleration. Selecting vibration wheel-water stable subbase void section acceleration waveform, fitting exciting force signal according to the above formula, and solving、And (3) with。
Step 6: substituting the exciting force signal and the vibrating wheel-water stable subbase void section acceleration signal into a water stable subbase vibration compaction dynamics model contact section theoretical control equation, and solving the vibrating wheel-water stable subbase contact force signal, wherein the solving equation is as follows:
（4）。
step 7: fitting a vibration wheel-water stable subbase contact force signal by using a cosine function, and solving the amplitude, phase and average value of the contact force signal, wherein the fitting formula is as follows:
（5）
In (5) In order to contact the force signal,For the vibration wheel-water stable sub-layer contact force average,For the vibration wheel-water stable sub-layer contact force amplitude,For the initial phase of the vibrating wheel-water stable sub-layer contact force,Is the excitation frequency; solving according to the fit vibration wheel-water stable subbase contact force signal 、 And (3) with 。
Step 8: substituting the contact force signal into a water-stable subbase theoretical displacement equation based on a three-degree-of-freedom multilayer elastoplastic dynamic model, fitting the water-stable subbase actual measurement displacement signal, and solving each elastoplastic mechanical parameter of the water-stable subbase, wherein the water-stable subbase theoretical displacement equation is as follows:
（6）
In (6) Is the displacement amplitude of the water-stable subbase layer,In order to shift the phase lag with respect to the contact force,Is the elastic rigidity coefficient of the water-stable subbase layer,Is the elastic viscosity coefficient of the water-stable subbing layer,Is the plastic rigidity coefficient of the water-stable subbase layer,The vibration wheel-water stable subbase contact force threshold value for plastically deforming the water stable subbase is the elastic rigidity coefficient of the subbase roadbed in the depth range of the vibration wheel. According to the above-mentioned fitting water-stable subbase actual measurement displacement signal, solving、 、 、 。
Step 9: by extracting the ratio of the contact force amplitude of the vibrating wheel-water stable subbase layer to the displacement amplitude of the water stable subbase layer, a mathematical relation formula of the comprehensive rigidity coefficient kcs of the water stable subbase layer and each elastoplastic mechanical parameter, namely a calculation formula of kcs is established, and the mathematical relation formula is as follows:
（7）。
the invention relies on the high-speed project along the river to carry out the comprehensive rigidity coefficient on site Is verified by the accuracy of (2). The model number of the road roller is Xu Gong XS265. The construction layer is a water-stabilized macadam subbase layer, and the loose pavement thickness is 30cm. Two compaction strips are selected, four sand filling method measuring points are selected for each strip, one measuring point is sequentially selected after each rolling is completed, and the compaction degree is measured and is used for establishing the correlation between each intelligent compaction control index (CMV, stiffness coefficient and comprehensive stiffness coefficient) and the compaction degree.
Coefficient of integrated stiffnessThe correlation with the compactness (R2=0.88) is obviously better than CMV, the rigidity coefficient (R2=0.68, R2=0.08), which shows that the comprehensive rigidity coefficient provided by the inventionThe index has higher accuracy and stability in the aspect of representing the lamination quality of the water-stable base.
Example 2: the computer readable storage medium of this embodiment has stored thereon a computer program which, when executed by a processor, performs the steps of the water stable under-layer intelligent compaction comprehensive stiffness coefficient index calculation and field application method of embodiment 1.
The computer readable storage medium of the present embodiment may be an internal storage unit of the terminal, for example, a hard disk or a memory of the terminal; the computer readable storage medium of the present embodiment may also be an external storage device of the terminal, for example, a plug-in hard disk, a smart memory card, a secure digital card, a flash memory card, etc. provided 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 the present embodiment is used to store a computer program and other programs and data required for a terminal, and the computer-readable storage medium may also be used to temporarily store data that has been output or is to be output.
Example 3: the computer device of this embodiment includes a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps in the water-stable under-layer intelligent compaction comprehensive stiffness coefficient index calculation and field application method of embodiment 1 when the program is executed.
In this embodiment, the processor may be a central processing unit, or may be other general purpose processors, digital signal processors, application specific integrated circuits, off-the-shelf programmable gate arrays or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like, where the general purpose processor may be a microprocessor or the processor may also be any conventional processor, or the like; the memory may include read only memory and random access memory, and provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory, e.g., the memory may also store information of the device type.
It will be appreciated by those skilled in the art that the embodiment(s) disclosure may be provided as a method, system, or computer program product. Thus, the present approach may take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present aspects may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, magnetic disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present aspects are described with reference to flowchart illustrations and/or block diagrams of methods and computer program products according to embodiments of the invention, it being understood that each flowchart illustration and/or block diagram illustration, and combinations of flowcharts and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions; these computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Those skilled in the art will appreciate that implementing all or part of the above-described methods in accordance with the embodiments may be accomplished by way of a computer program stored on a computer readable storage medium, which when executed may comprise the steps of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-only memory (ROM), a random access memory (RandomAccessMemory, RAM), or the like.
The examples of the present invention are merely for describing the preferred embodiments of the present invention, and are not intended to limit the spirit and scope of the present invention, and those skilled in the art should make various changes and modifications to the technical solution of the present invention without departing from the spirit of the present invention.
Claims (10)
1. The method for calculating the intelligent compaction comprehensive rigidity coefficient index of the water-stable subbase layer is characterized by comprising the following steps of:
Step 1: an acceleration sensor is arranged on the central shaft of the vibrating wheel and is used for collecting the acceleration signal of the vibrating wheel in the vibrating compaction process;
step 2: a GPS positioning device is arranged at the top of the cockpit of the road roller and used for collecting the position information of the road roller and reflecting the compaction quality distribution condition corresponding to the intelligent compaction stiffness index;
step 3: removing high-frequency components in the original acceleration signals, and carrying out noise reduction treatment on the original acceleration signals;
step 4: converting the acceleration signal into a displacement signal by adopting an average value method and a secondary integration method;
Step 5: substituting the acceleration signal of the shaking wheel-water stable subbase void section into a theoretical control equation of a shaking compaction dynamic model void section of the water stable subbase, and solving a shaking force signal of the shaking wheel;
step 6: substituting the exciting force signal and the acceleration signal of the shaking wheel-water stable subbase void section into a theoretical control equation of a contact section of a shaking compaction dynamic model of the water stable subbase, and solving the shaking wheel-water stable subbase contact force signal;
step 7: fitting a vibration wheel-water stable subbase contact force signal by using a cosine function, and solving the amplitude, phase and average value of the contact force signal, wherein the fitting formula is as follows:
（5）
In (5) In order to contact the force signal,For the mean value of the oscillating wheel-water stable sub-layer contact force signals,For the vibration wheel-water stable sub-layer contact force signal amplitude,For the initial phase of the vibrating wheel-water stable sub-layer contact force signal,Is the excitation frequency; solving according to the fit vibration wheel-water stable subbase contact force signal、And (3) with；
Step 8: substituting the contact force signal into a water-stable subbase theoretical displacement equation based on a three-degree-of-freedom multilayer elastoplastic dynamic model, fitting the water-stable subbase actual measurement displacement signal, and solving each elastoplastic mechanical parameter of the water-stable subbase;
Step 9: by extracting the ratio of the contact force amplitude of the vibrating wheel-water stable subbase layer to the displacement amplitude of the water stable subbase layer, the comprehensive rigidity coefficient of the water stable subbase layer is built Mathematical relation with each elastoplastic mechanical parameter, namely a calculation formula, is as follows:
（7）
in (7) Is the elastic rigidity coefficient of the water-stable subbase layer,Is the elastic viscosity coefficient of the water-stable subbing layer,Is the plastic rigidity coefficient of the water-stable subbase layer,Is the excitation frequency.
2. The method according to claim 1, wherein in step 1, the sampling frequency of the acceleration sensor is not lower than 1000Hz, and the measuring range is not lower than 100m/s2.
3. The method of claim 2, wherein in step 2, the GPS positioning device uses RTK real-time dynamic differential technology, installs a mobile station on top of the cockpit, installs a base station at the wide roadside view, and corrects the mobile station positioning information in real time using the base station positioning information.
4. The method according to claim 1, wherein in step 3, the noise reduction processing of the acceleration signal converts the acceleration time domain signal into the frequency domain signal by using a fast fourier transform method, removes high frequency components with a frequency higher than 140Hz, and converts the frequency domain signal after noise reduction into the time domain signal by using an inverse fast fourier transform method.
5. The method of claim 4, wherein in step 4, an average value method is used for correcting the acceleration signal, a 1s calculation interval is selected based on the principle that the average value of the acceleration signal in the calculation interval is0, the direct current component of the acceleration signal is removed, and the calculation formula is as follows:
（1）
In the formula (1) The direct current component of the acceleration is indicated,The value of the acceleration signal is indicated,Representing the number of sampling points within the computation interval,Representing the acceleration signal value after the DC component is removed;
the secondary integration method is used for obtaining a speed signal and a displacement signal, removing a direct current component based on the principle that the integration area is 0 in a single excitation period after integration, and splicing the speed signal and the displacement signal in each period into complete signals in sequence respectively; the calculation formula of the direct current component is as follows:
（2）
In the formula (2) As a direct current component of the power supply,For the excitation frequency to be the same,For the sampling frequency to be the same,For the number of sample points,In order to be a velocity or displacement signal value,The proportional length of the rest section of the single excitation period after the sampling point is intercepted.
6. The method of claim 5, wherein in step 5, the solution formula for solving the excitation force signal of the vibrating wheel is as follows:
(3)
In (3) For the amplitude of the exciting force,For the excitation frequency to be the same,For the initial phase of the exciting force,For the mass of the vibrating wheel,Is the quality of the machine frame,As the acceleration signal, the acceleration signal is,Gravitational acceleration; selecting vibration wheel-water stable subbase void section acceleration waveform, fitting exciting force signal according to the above formula, and solving、And (3) with。
7. The method of claim 6, wherein in step 6, the oscillating wheel-water stable sub-layer contact force signal is solvedThe solution formula of (2) is as follows:
（4）。
8. the method of claim 7, wherein in step 8, the water stable underlayment theoretical displacement equation is as follows:
(6）
In (6) Is the displacement amplitude of the water-stable subbase layer,In order to shift the phase lag with respect to the contact force,Is the elastic rigidity coefficient of the water-stable subbase layer,Is the elastic viscosity coefficient of the water-stable subbing layer,Is the plastic rigidity coefficient of the water-stable subbase layer,To plastically deform the water-stable underlayment,The elastic rigidity coefficient of the lower roadbed in the depth range is influenced by the vibrating wheel; according to the above-mentioned fitting water-stable subbase actual measurement displacement signal, solving、、、。
9. A computer-readable storage medium having stored thereon a computer program, characterized by: the program, when executed by a processor, implements the steps in the water-stable subbase intelligent compaction comprehensive stiffness coefficient index calculation and field application method as set forth in any one of claims 1-8.
10. A computer device comprising a memory, a processor and a computer program stored on the memory and operable on the processor, wherein the processor, when executing the program, performs the steps of the method for calculating and applying the index of the intelligent compaction integrated stiffness coefficient of the water-stable subbase according to any one of claims 1-8.
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 | 山东大学 | System and method for measuring rigidity coefficient and viscous damping coefficient of compacted soil body |
| CN114021066A (en) * | 2021-10-14 | 2022-02-08 | 东南大学 | Measuring and calculating method of low-variability modulus type intelligent compaction metering index |
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 | 山东大学 | System and method for measuring rigidity coefficient and viscous damping coefficient of compacted soil body |
| CN114021066A (en) * | 2021-10-14 | 2022-02-08 | 东南大学 | Measuring and calculating method of low-variability modulus type intelligent compaction metering index |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118094067A (en) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103411607B (en) | Pedestrian's step-size estimation and dead reckoning method | |
| CN100432681C (en) | Alternating Current Frequency Monitoring Method | |
| Park et al. | Displacement estimation using multimetric data fusion | |
| Nitsche et al. | Comparison of machine learning methods for evaluating pavement roughness based on vehicle response | |
| CN109357822A (en) | A kind of quick test and evaluation method of bridge changed based on Vehicle-Bridge Coupling System time-varying dynamic characteristic | |
| CN109858156B (en) | Vehicle and structure information simultaneous identification method based on axle coupling vibration | |
| CN110619106B (en) | Bridge damage positioning method and quantification method thereof | |
| CN109916485B (en) | Dynamic vehicle weighing method and device | |
| CN104268883A (en) | Time-frequency spectrum curve extracting method based on edge detection | |
| CN109791047A (en) | Mileage and velocity estimation | |
| CN105426587A (en) | Method for completing pavement condition acquisition and monitoring based on smart phone | |
| CN111191696B (en) | Deep learning algorithm-based steel bar layering method and system | |
| CN104697436A (en) | Fourier series-based circular induction synchronizer error model analysis method | |
| CN106786561A (en) | A kind of Low-frequency Oscillation Modal Parameters discrimination method based on adaptive Kalman filter | |
| CN107729592A (en) | Traced back the Time variable structure Modal Parameters Identification of track based on broad sense subspace | |
| CN111693283B (en) | Bearing unsteady state fault feature non-rotating speed extraction method | |
| CN104063543A (en) | Wheel-rail combination roughness identification method for rail traffic | |
| CN105388012A (en) | Order tracking method based on nonlinear frequency modulation wavelet transformation | |
| CN101975633A (en) | Method for measuring energizing force of engine by continuous refinement analytical Fourier transform method | |
| CN103217695A (en) | Capturing method of L1 and L2C signals for global positioning system (GPS) | |
| CN104698325A (en) | Method for determining low-frequency oscillating and mandatory oscillating through negative damping mechanism of power system | |
| CN102955004B (en) | Subway tunnel segment service performance detection method based on wave velocity determination | |
| CN102269803B (en) | Method for correcting low-frequency components in discrete spectrum based on time delay | |
| CN112229355A (en) | Bridge movement detection method and system based on signal processing and trailer response | |
| CN118094067B (en) | Calculation method for intelligent compaction comprehensive rigidity coefficient index of water-stable subbase layer |
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 |