CN106291708A - A kind of method and device revising data - Google Patents

Abstract

The present invention provides a data correction method and device, the method comprising: based on the exposed length of the test base on the surface of the rock and soil body, establishing a finite element structure model corresponding to the test base at each length; using the base The position node where the sensor is placed is used as the simulated monitoring node; the preset external excitation load is applied to the bottom node group on the corresponding finite element structure model, and the vibration data of each bottom node and the vibration data of each simulated monitoring point are obtained; according to each The vibration data of the bottom node and the vibration data of each simulated monitoring point are fitted with each data correction function; the vibration data measured by the base sensor is corrected according to each data correction function. In this way, according to each data correction function, the vibration data measured by the base sensor can be corrected, and the vibration data of the actual measuring point on the rock and soil mass can be obtained by solving, thereby improving the accuracy and accuracy of the data, and more accurately reflecting the Vibration characteristics of the actual measuring point on the rock and soil mass.

Description

Method and device for correcting data

technical field

The invention relates to the technical field of geotechnical testing, in particular to a method and device for correcting data.

Background technique

In the field vibration test of rock and soil, the vibration signal of rock and soil is generally collected by using the test sensor.

In the prior art, the fixing method of the vibration test sensor on the measuring point is generally to use gypsum and 502 glue to directly fix the sensor on the rock and soil mass test site to be tested. This method has four disadvantages: first, the fixation effect is not good for the loose rock-soil medium on the surface; The additional workload of the test; third, for special environments, such as tunnels containing accumulated water, short-term or perennial water accumulation on the surface (such as depressions, streams) and other environments that are not suitable for sensor work, the side walls of tunnels and rock and soil In the environment where it is difficult to place the sensor, the normal operation of the sensor is blocked; fourth, for the vibration signal at a certain depth of the rock and soil that needs to be tested, if the sensor is directly fixed on the surface of the measuring point without excavating to the measuring point, the measured The vibration signal cannot truly reflect the signal at the measuring point. Moreover, in the process of collecting vibration data generated by blasting operations, due to the dynamic response of the base structure under the action of blasting vibration, there is a certain deviation between the measured blasting vibration data and the actual monitoring object position, which affects the accuracy of the data. accuracy.

Based on this, there is an urgent need for a data correction method and device to solve the above-mentioned problems in the prior art.

Contents of the invention

Aiming at the problems existing in the prior art, the embodiment of the present invention provides a method and device for correcting data, so as to solve the problem in the prior art that when the vibration signal of the rock and soil body is measured by the sensor fixed on the base, due to blasting There is a certain deviation between the vibration data and the actual monitoring object position, which leads to technical problems that the accuracy of the data cannot be guaranteed.

The invention provides a method for correcting data, the method comprising

Based on the exposed length of the test base on the surface of rock and soil, establish the corresponding finite element structure model of the test base at each length;

Using the position node where the sensor is placed on the base as the simulated monitoring node;

Apply the preset external excitation load to the bottom node group on the corresponding finite element structure model to obtain vibration data of each bottom node and vibration data of each simulated monitoring point;

Fit each data correction function according to the vibration data of each bottom node and the vibration data of each simulated monitoring point;

The vibration data measured by the base sensor is corrected according to each data correction function; wherein, the data correction function specifically includes: a frequency correction function and a vibration velocity amplitude correction function.

In the above scheme, the establishment of the finite element structural model corresponding to the test base at each height specifically includes:

According to the actual size of each part of the base, use the dynamic finite element software ANSYS/LS-DYNA to establish the corresponding finite element structure model of the test base at each height; wherein,

The base includes:

a fixed abutment, one end of which is connected to the movable plate;

connect the high rod, one end of the connected high rod is connected with the other end of the fixed abutment;

a sensor, the sensor is installed on the fixed base;

A drill rod, one end of the drill rod is connected with the other end of the connecting rod.

In the above solution, when the length of the test base exposed in the rock and soil mass is 25 cm, the frequency correction function is specifically: f _X =0.489f _XM ^1.1951 , f _Y =0.489f _YM ^1.1951 and f _Z =1.6157 f _ZM ^0.9354 ;

The vibration velocity amplitude correction function is specifically: v _X =0.4185V _XM ^1.1538 , v _Y =0.4898V _YM ^0.961 and v _Z =0.9958V _ZM ^0.9982 ; wherein,

The f _x is the actual main vibration frequency of the measuring point in the X direction, the f _XM is the main vibration frequency of the sensor data in the X direction; the f _y is the actual main vibration of the measuring point in the Y direction Frequency, the f _YM is the main vibration frequency of the sensor data in the Y direction; the f _z is the actual main vibration frequency of the measuring point in the Z direction, and the f _ZM is the actual measuring point in the Z direction Main vibration frequency;

The v _x is the actual vibration velocity amplitude of the measuring point in the X direction, and the V _XM is the vibration velocity amplitude of the sensor data in the X direction; the f _y is the actual vibration velocity of the measuring point in the Y direction Vibration velocity amplitude, the V _YM is the vibration velocity amplitude of the sensor data in the Y direction; the f _z is the actual vibration velocity amplitude of the measuring point in the Z direction, and the V _ZM is the sensor data The vibration velocity amplitude of the data in the Z direction.

In the above solution, when the length of the test base exposed in the rock and soil mass is 30 cm, the frequency correction function is specifically: f _X =0.7272f _XM ^1.1332 , f _Y =2.3405f _YM ^0.8134 and f _Z =1.3762 f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.2181V _XM ^1.6255 , v _Y =0.3509V _YM ^1.2087 and v _Z =0.9953V _ZM ^0.9961 .

In the above scheme, when the length of the test base exposed in the rock and soil mass is 35m, the frequency correction function is specifically: f _X =2.2699f _XM ^0.8643 , f _Y =2.1075f _YM ^0.8949 and f _Z =0.8069 f _ZM ^1.1271 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0902V _XM ^2.2975 , v _Y =0.2412V _YM ^1.559 and v _Z =0.9968V _ZM ^0.9934 .

In the above solution, when the length of the test base exposed in the rock and soil mass is 40m, the frequency correction function is specifically: f _X =4.8804f _XM ^0.7062 , f _Y =3.2278f _YM ^0.7639 and f _Z =1.3762 f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0665V _XM ^2.6674 , v _Y =0.2127V _YM ^1.743 and v _Z =0.9996V _ZM ^0.9906 .

In the above scheme, when the length of the test base exposed in the rock and soil mass is 45m, the frequency correction function is specifically: f _X =12.601f _XM ^0.435 , f _Y =3.5085f _YM ^0.749 and f _Z =1.3762 f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0794V _XM ^2.9451 , v _Y =0.1738V _YM ^2.0956 and v _Z =1.0032V _ZM ^0.9848 .

In the above solution, when the length of the test base exposed in the rock and soil mass is 50m, the frequency correction function is specifically: f _X =10.822f _XM ^0.4837 , f _Y =3.3109f _YM ^0.7379 and f _Z =1.3762 f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0539V _XM ^3.7508 , v _Y =0.1628V _YM ^2.4078 and v _Z =1.0044V _ZM ^0.9824 .

In the above solution, when the length of the test base exposed in the rock and soil mass is 55m, the frequency correction function is specifically: f _X =10.068f _XM ^0.4822 , f _Y =4.7793f _YM ^0.6468 and f _Y =4.7793 f _YM ^0.6468 ;

The vibration velocity amplitude correction function is specifically: v _X =0.1545V _XM ^3.3643 , v _Y =0.352V _YM ^1.9237 and v _Z =0.989V _ZM ^0.984 .

The present invention also provides a device for correcting data, the device comprising:

Establishing a unit, which is used to establish a finite element structural model corresponding to the test base at each length based on the exposed length of the test base on the surface of the rock and soil body;

An applying unit, the applying unit is used to apply a preset external excitation load to the bottom node group on the corresponding finite element structural model, and obtain vibration data of each bottom node and vibration data of each simulated monitoring point;

A fitting unit, the fitting unit is used to fit each data correction function according to the vibration data of each bottom node and the vibration data of each simulated monitoring point;

A correction unit, the correction unit is used to correct the vibration data measured by the base sensor according to each data correction function; wherein, the data correction function specifically includes: a frequency correction function and a vibration velocity amplitude correction function; the simulation The monitoring node is a location node where sensors are placed on the base.

The present invention provides a data correction method and device, the method comprising: based on the exposed length of the test base on the surface of the rock and soil body, establishing a finite element structure model corresponding to the test base at each length; using the base The position node where the sensor is placed is used as the simulated monitoring node; the preset external excitation load is applied to the bottom node group on the corresponding finite element structure model, and the vibration data of each bottom node and the vibration data of each simulated monitoring point are obtained; according to each bottom node The vibration data and the vibration data of each simulated monitoring point are fitted with each data correction function; according to each data correction function, the vibration data of different depths of the measured rock and soil mass are corrected; wherein, the data correction function specifically includes: a frequency correction function and a vibration velocity Amplitude correction function. In this way, the vibration data measured by the base sensor can be corrected according to each data correction function, and the vibration data of the actual measuring point on the rock and soil mass can be obtained by solving, thereby improving the precision and accuracy of the data, and more accurately reflecting The vibration characteristics of the actual measuring points on the rock and soil.

Description of drawings

FIG. 1 is a schematic flowchart of a method for correcting data provided by Embodiment 1 of the present invention;

FIG. 2 is a schematic diagram of the overall structure of the base provided by Embodiment 1 of the present invention;

Fig. 3 is the top view of the movable plate provided by the embodiment of the present invention;

Fig. 4 is a top view of the fixed base provided by the embodiment of the present invention;

Fig. 5 is the side view of the connecting rod provided by the embodiment of the present invention;

Fig. 6 is a top view of the connecting pole provided by the embodiment of the present invention;

Fig. 7 is the side view of the drilling rod provided by the embodiment of the present invention;

Fig. 8 is the top view of the drilling rod provided by the embodiment of the present invention;

Fig. 9 is a top view of the height-connecting cylinder provided by the embodiment of the present invention;

Fig. 10 is a top view of the force applying component provided by the embodiment of the present invention;

FIG. 11 is a schematic structural diagram of a device for correcting data provided by Embodiment 2 of the present invention;

Fig. 12 is a schematic diagram of the base provided by the third embodiment of the present invention fixed on the weathered rock.

detailed description

In order to improve the precision and accuracy of the test data, so that the test data can more accurately reflect the vibration characteristics of the actual measuring point, the present invention provides a method and device for data correction. The method includes: The length exposed in the soil, establish the finite element structure model corresponding to the test base at each height; use the position node where the sensor is placed on the base as the simulated monitoring node; apply the preset external excitation load to the corresponding finite element structure On the bottom node group on the model, the vibration data of each bottom node and the vibration data of each simulated monitoring point are obtained; each data correction function is fitted according to the vibration data of each bottom node and each simulated monitoring point; according to each data correction function, the base The vibration data measured by the sensor is corrected; wherein, the data correction function specifically includes: a frequency correction function and a vibration velocity amplitude correction function.

The technical solution of the present invention will be further described in detail below with reference to the drawings and specific embodiments.

Embodiment one

This embodiment provides a method for data correction, as shown in Figure 1, the method includes the following steps:

Step 110, based on the exposed length of the test base on the surface of the rock and soil body, a finite element structure model corresponding to the test base at each length is established.

In this step, according to the actual size of each component of the base, use the dynamic finite element software ANSYS/LS-DYNA to establish the finite element structure model corresponding to the test base at each height; the materials in the model are divided into Lagrange nets by Soild164 elements grid. In order to ensure a uniform grid, all grid sizes are divided into 0.2cm, and the calculation uses the cm-g-us unit system. Among them, the base hardware structure materials are all made of Q235 ordinary carbon structural steel, the material constitutive model is MAT_ELASTIC, and its specific parameters are as follows: density 7.85g/cm3, elastic modulus 200GPa, Poisson's ratio 0.3, tensile strength 420MPa .

The exposed length of the base in the rock and soil mass may include: 25cm, 30cm, 35cm, 40cm, 45cm, 50cm, 55cm.

Wherein, as shown in Figure 2, the base includes a fixed base, a movable plate 1, a high rod 2, a sensor 3, and a drill rod 4; wherein,

One end of the fixed abutment is connected with the movable plate 1; one end of the connecting rod 2 is connected with the other end of the fixed abutment; the sensor 3 is installed on the fixed abutment for testing rock and soil The vibration signal of the body; one end of the drilling rod 4 is connected with the other end of the extension rod 2 for screwing the base into the rock and soil body. The sensor 3 is a vibration sensor.

Specifically, referring to Fig. 2 , the fixed abutment includes: a fixed plate 5, a base block 6 and a height connecting tube 7; the fixed plate 5 is arranged with a first screw along a diagonal line, and the first screw includes two , one corner of the fixed plate 5 is connected with one corner of the movable plate 1 by a first screw, and the other corner of the fixed plate 5 is connected with the other corner of the movable plate 1 by a first screw; the base block 6 One end of the base block 6 is connected to the other end of the fixing plate 5 by welding, and the other end of the base block 6 is connected to one end of the connecting rod 2 by threads. One end of the drill rod 4 is connected to the other end of the connecting rod 2 through the connecting tube 7 .

Wherein, referring to Fig. 3, a slot is provided on the movable plate 1 for dredging the wiring of the sensor 3; in the present embodiment, the diameter of the slot is 2 cm, and the thickness is 0.5 cm; the movable plate The height of the center enclosure of 1 is 0.5cm, and the thickness is 0.2cm. The movable plate 1 is also provided with a wing wall with an inner side length of 6cm to limit the displacement of the sensor 3. In other embodiments, the diameter and thickness of the slot; the center height, thickness and side length of the wing wall of the movable panel 1 can be set according to actual needs.

Referring to FIG. 5 , a slot is opened on the fixed base for dredging the wiring of the sensor 3 ; the diameter of the slot is 2 cm, and the thickness is 0.5 cm.

6 and 7, one end of the connecting rod 2 is also provided with three screw holes and threads. In this embodiment, the diameter of the thread is 2.5 cm, and the outer diameter of the connecting rod 2 is 5 cm. Scale marks are also provided on the connecting rod 2 . In other embodiments, the thread diameter and outer diameter of the extension rod 2 can be set according to actual needs.

8 and 9, one end of the drill rod 4 is also provided with three screw holes and threads. In this embodiment, the length of the drill rod 4 is 25 cm, and the tapered part at the other end is 10 cm. Scale marks are also arranged on the drill rod 4 . In other embodiments, the length of the drill rod 4 and the length of the tapered portion can be set according to actual needs.

Referring to Fig. 10, in the present embodiment, the length of the connecting tube 7 is 5 cm, the inner diameter is 1.5 cm, and the outer diameter is 5 cm; in other embodiments, the length, inner diameter and outer diameter of the drill rod connecting tube 7 are It can be set according to actual needs.

Here, in order to measure the vibration signal of the rock-soil mass under different coating thicknesses or ponding depths, the connecting rods 2 may include a plurality, and the lengths of each connecting rod 2 are different.

Step 111, taking the position node where the sensor is placed on the base as the simulated monitoring node.

In this step, in the established finite element structure model, the position node where the sensor is placed on the base is used as the simulated monitoring node to more realistically simulate the actual test environment, so that the simulated monitoring node data can truly reflect the actual test environment. sensor data.

Step 112, applying preset external excitation loads to the bottom node groups on the corresponding finite element structural model, and obtaining vibration data of each bottom node and vibration data of each simulated monitoring point.

In this step, in order to truly reflect the vibration data of rock and soil in the actual test environment, the preset external excitation load is applied to the bottom node group on the corresponding finite element structural model, and the vibration data and Vibration data of each simulated monitoring point.

Among them, the external excitation load specifically includes: 20Hz, 30Hz, 40Hz, 50Hz, 60Hz, 80Hz, 100Hz main vibration frequency; 1cm/s, 2cm/s, 3cm/s, 4cm/s, 5cm/s, 6cm/s, The vibration velocity amplitude of 7cm/s. The main vibration frequency and vibration velocity amplitude are combined to form forty-nine sets of excitation loads, which are applied to the corresponding finite element structural model of each exposed length, for example, 20Hz and 1cm/s, 2cm/s, 3cm/s respectively , 4cm/s, 5cm/s, 6cm/s, and 7cm/s form seven sets of excitation loads, and apply this set of excitation loads to the finite element structure model with an exposed length of 25cm, and so on, a total of forty-nine sets of excitation loads load.

Step 113, fitting each data correction function according to the vibration data of each bottom node and the vibration data of each simulated monitoring point.

In this step, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 25cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points . According to the least squares method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 25cm is obtained. The data correction function specifically includes: a frequency correction function and a vibration velocity amplitude correction function, specifically as follows:

The frequency correction function is specifically: f _X =0.489f _XM ^1.1951 , f _Y =0.489f _YM ^1.1951 and f _Z =1.6157f _ZM ^0.9354 ;

The vibration velocity amplitude correction function is specifically: v _X =0.4185V _XM ^1.1538 , v _Y =0.4898V _YM ^0.961 and v _Z =0.9958V _ZM ^0.9982 ; wherein,

The f _x is the actual main vibration frequency of the measuring point in the X direction, the f _XM is the main vibration frequency of the sensor data in the X direction; the f _y is the actual main vibration of the measuring point in the Y direction Frequency, the f _YM is the main vibration frequency of the sensor data in the Y direction; the f _z is the actual main vibration frequency of the measuring point in the Z direction, and the f _ZM is the actual measuring point in the Z direction Main vibration frequency;

The v _x is the actual vibration velocity amplitude of the measuring point in the X direction, and the V _XM is the vibration velocity amplitude of the sensor data in the X direction; the f _y is the actual vibration velocity of the measuring point in the Y direction Vibration velocity amplitude, the V _YM is the vibration velocity amplitude of the sensor data in the Y direction; the f _z is the actual vibration velocity amplitude of the measuring point in the Z direction, and the V _ZM is the sensor data The vibration velocity amplitude of the data in the Z direction.

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 30 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 30cm is obtained, as follows:

When the length of the test base exposed in the rock and soil mass is 30 cm, the frequency correction function is specifically: f _X =0.7272f _XM ^1.1332 , f _Y =2.3405f _YM ^0.8134 and f _Z =1.3762f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.2181V _XM ^1.6255 , v _Y =0.3509V _YM ^1.2087 and v _Z =0.9953V _ZM ^0.9961 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 35 cm, 49 sets of calculation results can be obtained, including the bottom node vibration data and the simulated monitoring point vibration data. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 35cm is obtained, as follows:

The frequency correction function is specifically: f _X =2.2699f _XM ^0.8643 , f _Y =2.1075f _YM ^0.8949 and f _Z =0.8069f _ZM ^1.1271 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0902V _XM ^2.2975 , v _Y =0.2412V _YM ^1.559 and v _Z =0.9968V _ZM ^0.9934 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 40 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 40cm is obtained, as follows:

The frequency correction function is specifically: f _X =4.8804f _XM ^0.7062 , f _Y =3.2278f _YM ^0.7639 and f _Z =1.3762f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0665V _XM ^2.6674 , v _Y =0.2127V _YM ^1.743 and v _Z =0.9996V _ZM ^0.9906 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 45 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 45cm is obtained, as follows:

The frequency correction function is specifically: f _X =12.601f _XM ^0.435 , f _Y =3.5085f _YM ^0.749 and f _Z =1.3762f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0794V _XM ^2.9451 , v _Y =0.1738V _YM ^2.0956 and v _Z =1.0032V _ZM ^0.9848 .

Correspondingly, when a set of excitation loads with a main vibration frequency of 80 Hz and a vibration velocity amplitude of 6 cm/s are applied to the finite element structural model with an exposed length of 50 cm, 49 sets of calculation results can be obtained, including Bottom node vibration data and simulated monitoring point vibration data. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 50cm is obtained, as follows:

The frequency correction function is specifically: f _X =10.822f _XM ^0.4837 , f _Y =3.3109f _YM ^0.7379 and f _Z =1.3762f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0539V _XM ^3.7508 , v _Y =0.1628V _YM ^2.4078 and v _Z =1.0044V _ZM ^0.9824 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 50 cm, 49 sets of calculation results can be obtained, including the bottom node vibration data and the simulated monitoring point vibration data. According to the least square method, the vibration data of the internal nodes and the vibration data of the simulated monitoring points are fitted, and the correction function when the exposed length is 55cm is obtained, as follows:

The frequency correction function is specifically: f _X =10.068f _XM ^0.4822 , f _Y =4.7793f _YM ^0.6468 and f _Y =4.7793f _YM ^0.6468 ;

The vibration velocity amplitude correction function is specifically: v _X =0.1545V _XM ^3.3643 , v _Y =0.352V _YM ^1.9237 and v _Z =0.989V _ZM ^0.984 .

Step 114, correcting the vibration data measured by the base sensor according to various data correction functions.

In this step, after each data correction function is obtained, the vibration data measured by the base sensor can be corrected according to each data correction function.

The data correction method provided in this embodiment is based on the exposed length of the test base on the surface of the rock and soil body, and establishes the corresponding finite element structural model of the test base at each length, and according to the vibration data of each bottom node and each The vibration data of the simulated monitoring points are fitted to each data correction function; according to each data correction function, the vibration data measured by the base sensor can be corrected, and the vibration data of the actual measurement points on the rock and soil can be obtained by solving, thereby improving the accuracy of the data. Accuracy and accuracy can more accurately reflect the vibration characteristics of the actual measuring point.

Embodiment two

Corresponding to Embodiment 1, this embodiment provides a device for correcting data. As shown in FIG. 11 , the device includes: an establishment unit 101, an application unit 102, a fitting unit 103, and a correction unit 104; wherein,

The establishment unit 101 is configured to establish a finite element structure model corresponding to the test base at each length based on the length of the test base exposed on the surface of the rock and soil body.

Specifically, the establishment unit 101 uses the dynamic finite element software ANSYS/LS-DYNA to establish the corresponding finite element structural model of the test base at each height according to the actual size of each component of the base; the materials in the model all use the Soild164 unit Divided into a Lagrange grid. In order to ensure a uniform grid, all grid sizes are divided into 0.2cm, and the calculation uses the cm-g-us unit system. Among them, the base hardware structure materials are all made of Q235 ordinary carbon structural steel, the material constitutive model is MAT_ELASTIC, and its specific parameters are as follows: density 7.85g/cm3, elastic modulus 200GPa, Poisson's ratio 0.3, tensile strength 420MPa .

The exposed length of the base in the rock and soil mass may include: 25cm, 30cm, 35cm, 40cm, 45cm, 50cm, 55cm.

Wherein, as shown in Figure 2, the base includes a fixed base, a movable plate 1, a high rod 2, a sensor 3, and a drill rod 4; wherein,

One end of the fixed abutment is connected with the movable plate 1; one end of the connecting rod 2 is connected with the other end of the fixed abutment; the sensor 3 is installed on the fixed abutment for testing rock and soil The vibration signal of the body; one end of the drilling rod 4 is connected with the other end of the extension rod 2 for screwing the base into the rock and soil body. The sensor 3 is a vibration sensor.

Specifically, referring to Fig. 2 , the fixed abutment includes: a fixed plate 5, a base block 6 and a height connecting tube 7; the fixed plate 5 is arranged with a first screw along a diagonal line, and the first screw includes two , one corner of the fixed plate 5 is connected with one corner of the movable plate 1 by a first screw, and the other corner of the fixed plate 5 is connected with the other corner of the movable plate 1 by a first screw; the base block 6 One end of the base block 6 is connected to the other end of the fixing plate 5 by welding, and the other end of the base block 6 is connected to one end of the connecting rod 2 by threads. One end of the drill rod 4 is connected to the other end of the connecting rod 2 through the connecting tube 7 .

Wherein, referring to Fig. 3, a slot is provided on the movable plate 1 for dredging the wiring of the sensor 3; in the present embodiment, the diameter of the slot is 2 cm, and the thickness is 0.5 cm; the movable plate The height of the center enclosure of 1 is 0.5cm, and the thickness is 0.2cm. The movable plate 1 is also provided with a wing wall with an inner side length of 6cm to limit the displacement of the sensor 3. In other embodiments, the diameter and thickness of the slot; the center height, thickness and side length of the wing wall of the movable panel 1 can be set according to actual needs.

Referring to FIG. 5 , a slot is opened on the fixed base for dredging the wiring of the sensor 3 ; the diameter of the slot is 2 cm, and the thickness is 0.5 cm.

6 and 7, one end of the connecting rod 2 is also provided with three screw holes and threads. In this embodiment, the diameter of the thread is 2.5 cm, and the outer diameter of the connecting rod 2 is 5 cm. Scale marks are also provided on the connecting rod 2 . In other embodiments, the thread diameter and outer diameter of the extension rod 2 can be set according to actual needs.

8 and 9, one end of the drill rod 4 is also provided with three screw holes and threads. In this embodiment, the length of the drill rod 4 is 25 cm, and the tapered part at the other end is 10 cm. Scale marks are also arranged on the drill rod 4 . In other embodiments, the length of the drill rod 4 and the length of the tapered portion can be set according to actual needs.

Referring to Fig. 10, in the present embodiment, the length of the connecting tube 7 is 5 cm, the inner diameter is 1.5 cm, and the outer diameter is 5 cm; in other embodiments, the length, inner diameter and outer diameter of the drill rod connecting tube 7 are It can be set according to actual needs.

Here, in order to measure the vibration signal of the rock-soil mass under different coating thicknesses or ponding depths, the connecting rods 2 may include a plurality, and the lengths of each connecting rod 2 are different.

After the finite element structural model of the base is established, the applying unit 102 is used to apply the preset external excitation load to the bottom node group on the corresponding finite element structural model, and obtain the vibration data of each bottom node and Vibration data of each simulated monitoring point;

Specifically, the external excitation load specifically includes: main vibration frequencies of 20Hz, 30Hz, 40Hz, 50Hz, 60Hz, 80Hz, 100Hz; 1cm/s, 2cm/s, 3cm/s, 4cm/s, 5cm/s, 6cm/s , 7cm/s vibration velocity amplitude. The main vibration frequency and vibration velocity amplitude are combined to form forty-nine sets of excitation loads, which are applied to the corresponding finite element structural model of each exposed length, for example, 20Hz and 1cm/s, 2cm/s, 3cm/s respectively , 4cm/s, 5cm/s, 6cm/s, and 7cm/s form seven sets of excitation loads, and apply this set of excitation loads to the finite element structure model with an exposed length of 25cm, and so on, a total of forty-nine sets of excitation loads load.

The application unit 102 uses the position node where the sensor is placed on the base as a simulated monitoring node to more realistically simulate the actual test environment, so that the simulated monitoring node data can truly reflect the sensor data in the actual test environment. In order to truly reflect the vibration data of rock and soil in the actual test environment, the preset external excitation load is applied to the bottom node group on the corresponding finite element structural model, and the vibration data of each bottom node and each simulated monitoring point are obtained vibration data.

After the applying unit 102 acquires the vibration data of each bottom node and each simulated monitoring point, the fitting unit 103 is used to fit each data correction function according to each bottom node vibration data and each simulated monitoring point vibration data.

Specifically, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 25 cm, 49 sets of calculation results can be obtained, including vibration data of bottom nodes and vibration data of simulated monitoring points. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 25 cm. The data correction function specifically includes: a frequency correction function and a vibration velocity The amplitude correction function is as follows:

The frequency correction function is specifically: f _X =0.489f _XM ^1.1951 , f _Y =0.489f _YM ^1.1951 and f _Z =1.6157f _ZM ^0.9354 ;

The vibration velocity amplitude correction function is specifically: v _X =0.4185V _XM ^1.1538 , v _Y =0.4898V _YM ^0.961 and v _Z =0.9958V _ZM ^0.9982 ; wherein,

The f _x is the actual main vibration frequency of the measuring point in the X direction, the f _XM is the main vibration frequency of the sensor data in the X direction; the f _y is the actual main vibration of the measuring point in the Y direction Frequency, the f _YM is the main vibration frequency of the sensor data in the Y direction; the f _z is the actual main vibration frequency of the measuring point in the Z direction, and the f _ZM is the actual measuring point in the Z direction Main vibration frequency;

The v _x is the actual vibration velocity amplitude of the measuring point in the X direction, and the V _XM is the vibration velocity amplitude of the sensor data in the X direction; the f _y is the actual vibration velocity of the measuring point in the Y direction Vibration velocity amplitude, the V _YM is the vibration velocity amplitude of the sensor data in the Y direction; the f _z is the actual vibration velocity amplitude of the measuring point in the Z direction, and the V _ZM is the sensor data The vibration velocity amplitude of the data in the Z direction.

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 30 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 30 cm, as follows:

When the length of the test base exposed in the rock and soil mass is 30 cm, the frequency correction function is specifically: f _X =0.7272f _XM ^1.1332 , f _Y =2.3405f _YM ^0.8134 and f _Z =1.3762f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.2181V _XM ^1.6255 , v _Y =0.3509V _YM ^1.2087 and v _Z =0.9953V _ZM ^0.9961 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 35 cm, 49 sets of calculation results can be obtained, including the bottom node vibration data and the simulated monitoring point vibration data. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 35 cm, as follows:

The frequency correction function is specifically: f _X =2.2699f _XM ^0.8643 , f _Y =2.1075f _YM ^0.8949 and f _Z =0.8069f _ZM ^1.1271 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0902V _XM ^2.2975 , v _Y =0.2412V _YM ^1.559 and v _Z =0.9968V _ZM ^0.9934 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 40 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 40 cm, as follows:

The frequency correction function is specifically: f _X =4.8804f _XM ^0.7062 , f _Y =3.2278f _YM ^0.7639 and f _Z =1.3762f _ZM ^0.9894 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0665V _XM ^2.6674 , v _Y =0.2127V _YM ^1.743 and v _Z =0.9996V _ZM ^0.9906 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 45 cm, 49 sets of calculation results can be obtained, including the vibration data of the bottom nodes and the vibration data of the simulated monitoring points. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 45 cm, as follows:

The frequency correction function is specifically: f _X =12.601f _XM ^0.435 , f _Y =3.5085f _YM ^0.749 and f _Z =1.3762f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0794V _XM ^2.9451 , v _Y =0.1738V _YM ^2.0956 and v _Z =1.0032V _ZM ^0.9848 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 50 cm, 49 sets of calculation results can be obtained, including the bottom node vibration data and the simulated monitoring point vibration data. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 50 cm, as follows:

The frequency correction function is specifically: f _X =10.822f _XM ^0.4837 , f _Y =3.3109f _YM ^0.7379 and f _Z =1.3762f _ZM ^0.9849 ;

The vibration velocity amplitude correction function is specifically: v _X =0.0539V _XM ^3.7508 , v _Y =0.1628V _YM ^2.4078 and v _Z =1.0044V _ZM ^0.9824 .

Correspondingly, when forty-nine sets of excitation loads are respectively applied to the finite element structural model with an exposed length of 50 cm, 49 sets of calculation results can be obtained, including the bottom node vibration data and the simulated monitoring point vibration data. The fitting unit 103 fits the vibration data of the internal nodes and the vibration data of the simulated monitoring points according to the least square method, and obtains a correction function when the exposed length is 55 cm, as follows:

The frequency correction function is specifically: f _X =10.068f _XM ^0.4822 , f _Y =4.7793f _YM ^0.6468 and f _Y =4.7793f _YM ^0.6468 ;

The vibration velocity amplitude correction function is specifically: v _X =0.1545V _XM ^3.3643 , v _Y =0.352V _YM ^1.9237 and v _Z =0.989V _ZM ^0.984 .

After each data correction function is obtained, the correction unit 104 is configured to correct the vibration data measured by the base sensor according to each data correction function.

In practical applications, the establishment unit 101, the application unit 102, the fitting unit 103, and the correction unit 104 may be composed of a central processing unit (CPU, Central Processing Unit), a digital signal processor (DSP, DigitalSignal Processor), Realized by programmable logic array (FPGA, Field Programmable Gate Array) and micro control unit (MCU, Micro Controller Unit).

The data correction device provided in this embodiment is based on the exposed length of the test base on the surface of the rock and soil body, and establishes the corresponding finite element structure model of the test base at each length, and according to the vibration data of each bottom node and each The vibration data of the simulated monitoring point is fitted with each data correction function; according to each data correction function, the vibration data measured by the base sensor can be corrected, and the vibration data of the actual measurement point on the rock and soil can be obtained by solving, thereby improving the accuracy of the data And accuracy, it can more accurately reflect the vibration characteristics of the actual measuring point.

Embodiment Three

In practical applications, the base of Embodiment 1 can be used to test rock and soil bodies in different environments to obtain vibration data, and the data correction method provided in Embodiment 1 and the data correction device provided in Embodiment 2 can be used to correct the data. ,details as follows:

When the rock and soil body is weathered rock, the sensor is not easy to fix. As shown in FIG. 12 , a first sensor 121 is arranged on the base to monitor the vibration signal. In order to verify the accuracy of the correction algorithm, the upper weathered layer of the rock is excavated to form a trench, and the second sensor 122 is arranged, and the measured vibration signal is compared with the vibration amplitude and frequency corrected by the correction algorithm. During actual monitoring, the exposed length of the base is 55cm. The final vibration monitoring data obtained are shown in Table 1:

Table 1

The data measured by the first sensor 121 are corrected using the data correction function with a base exposed length of 55 cm as shown in Table 2:

Table 2

It can be seen from Table 2 that the maximum error between the corrected vibration velocity of the first sensor 121 and the vibration velocity measured by the second sensor 122 is 2.81%, and the maximum frequency error is 4.02%, indicating the accuracy of the data correction method Still good.

When the rock and soil body is loose soil, the sensor is not easy to fix, and the base is used to monitor the ground vibration in the loose soil environment. According to the actual environment on site, the exposed height of the base is selected as 25cm. After obtaining the monitoring data (monitoring vibration velocity and monitoring frequency) of the first sensor 121, the data correction function that the exposed length of the base is 25cm is used to measure the first sensor 121. After the data is corrected, the corrected data is shown in Table 3:

table 3

When using the base to monitor the rock and soil under the water accumulation environment, according to the actual environment on site, the exposed height of the base is selected as 30cm, and after obtaining the monitoring data (monitoring vibration velocity and monitoring frequency) of the first sensor 121, use After the data measured by the first sensor 121 is corrected by the data correction function with a base exposed length of 30 cm, the corrected data are shown in Table 4:

Table 4

In this embodiment, after the data measured by the first sensor 121 is corrected by using the data correction method provided in the first embodiment, the vibration characteristics of the actual measuring point can be more accurately reflected.

The above description is only a preferred embodiment of the present invention, and is not used to limit the protection scope of the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention shall be included in the within the protection scope of the present invention.

Claims (10)

1. A method of modifying data, the method comprising

Establishing a finite element structure model corresponding to the test base under each length based on the exposed length of the test base on the surface of the rock-soil body;

taking a position node of the base where the sensor is placed as a simulation monitoring node;

applying a preset external excitation load to the bottom node group on the corresponding finite element structure model to obtain vibration data of each bottom node and vibration data of each simulation monitoring point;

fitting each data correction function according to each bottom node vibration data and each simulation monitoring point vibration data;

correcting the vibration data measured by the base sensor according to each data correction function; wherein the data modification function specifically includes: a frequency correction function and a vibration velocity amplitude correction function.

2. The method of claim 1, wherein establishing the finite element structure model corresponding to the test base at each height comprises:

establishing a finite element structure model corresponding to the test base under each height by using dynamic finite element software ANSYS/LS-DYNA according to the actual size of each component of the base; wherein,

the base includes:

one end of the fixed base platform is connected with the movable plate;

one end of the height connecting rod is connected with the other end of the fixed base station;

a sensor mounted on the fixed base;

and one end of the drill rod is connected with the other end of the heightening rod.

3. The method according to claim 2, characterized in that when the length of the test base exposed in the geotechnical body is 25cm, the frequency correction function is in particular: f. of_X＝0.489f_XM ^1.1951、f_Y＝0.489f_YM ^1.1951And f_Z＝1.6157f_ZM ^0.9354；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.4185V_XM ^1.1538、v_Y＝0.4898V_YM ^0.961And v_Z＝0.9958V_ZM ^0.9982(ii) a Wherein,

f is_xFor measuring the actual principal vibration frequency of the point in the X direction, said f_XMFor the sensor data in the X directionA primary vibration frequency; f is_yFor measuring the actual principal vibration frequency in the Y direction, said f_YMThe main vibration frequency of the sensor data in the Y direction; f is_zFor measuring the actual principal vibration frequency in the Z direction, said f_ZMThe main vibration frequency of the actual measuring point in the Z direction is obtained;

v is_xFor the actual vibration velocity amplitude of the measuring point in the X direction, said V_XMA vibration velocity amplitude in the X direction is taken as the sensor data; f is_yFor the actual vibration velocity amplitude of the measured point in the Y direction, said V_YMThe vibration speed amplitude of the sensor data in the Y direction is obtained; f is_zFor the actual vibration velocity amplitude of the measuring point in the Z direction, said V_ZMThe vibration velocity amplitude in the Z direction is the sensor data.

4. The method according to claim 3, characterized in that when the length of the test base exposed in the geotechnical body is 30cm, the frequency correction function is in particular: f. of_X＝0.7272f_XM ^1.1332、f_Y＝2.3405f_YM ^0.8134And f_Z＝1.3762f_ZM ^0.9894；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.2181V_XM ^1.6255、v_Y＝0.3509V_YM ^1.2087And v_Z＝0.9953V_ZM ^0.9961。

5. The method according to claim 3, characterized in that when the length of the test base exposed in the rock-soil mass is 35m, the frequency correction function is in particular: f. of_X＝2.2699f_XM ^0.8643、f_Y＝2.1075f_YM ^0.8949And f_Z＝0.8069f_ZM ^1.1271；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.0902V_XM ^2.2975、v_Y＝0.2412V_YM ^1.559And v_Z＝0.9968V_ZM ^0.9934。

6. The method according to claim 3, characterized in that when the length of the test base exposed in the rock-soil mass is 40m, the frequency correction function is in particular: f. of_X＝4.8804f_XM ^0.7062、f_Y＝3.2278f_YM ^0.7639And f_Z＝1.3762f_ZM ^0.9894；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.0665V_XM ^2.6674、v_Y＝0.2127V_YM ^1.743And v_Z＝0.9996V_ZM ^0.9906。

7. The method according to claim 3, characterized in that when the length of the test base exposed in the rock-soil mass is 45m, the frequency correction function is in particular: f. of_X＝12.601f_XM ^0.435、f_Y＝3.5085f_YM ^0.749And f_Z＝1.3762f_ZM ^0.9849；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.0794V_XM ^2.9451、v_Y＝0.1738V_YM ^2.0956And v_Z＝1.0032V_ZM ^0.9848。

8. The method according to claim 3, characterized in that when the length of the test base exposed in the rock-soil mass is 50m, the frequency correction function is in particular: f. of_X＝10.822f_XM ^0.4837、f_Y＝3.3109f_YM ^0.7379And f_Z＝1.3762f_ZM ^0.9849；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.0539V_XM ^3.7508、v_Y＝0.1628V_YM ^2.4078And v_Z＝1.0044V_ZM ^0.9824。

9. The method according to claim 3, characterized in that when the test base is exposed in the geotechnical body for a length of 55m, the frequency correction function is in particular: f. of_X＝10.068f_XM ^0.4822、f_Y＝4.7793f_YM ^0.6468And f_Y＝4.7793f_YM ^0.6468；

The vibration velocity amplitude correction function is specifically as follows: v. of_X＝0.1545V_XM ^3.3643、v_Y＝0.352V_YM ^1.9237And v_Z＝0.989V_ZM ^0.984。

10. An apparatus for modifying data, the apparatus comprising:

the establishing unit is used for establishing a finite element structure model corresponding to the test base under each length based on the exposed length of the test base on the surface of the rock-soil body;

the application unit is used for applying a preset external excitation load to the bottom node group on the corresponding finite element structure model to obtain vibration data of each bottom node and vibration data of each simulation monitoring point;

the fitting unit is used for fitting each data correction function according to each bottom node vibration data and each simulation monitoring point vibration data;

the correction unit is used for correcting the vibration data measured by the base sensor according to each data correction function; wherein the data modification function specifically includes: a frequency correction function and a vibration velocity amplitude correction function; the simulation monitoring node is a position node of the base for placing the sensor.