Next Article in Journal
A Secure Video Steganography Based on the Intra-Prediction Mode (IPM) for H264
Next Article in Special Issue
Wind Turbine Gearbox Condition Monitoring Based on Class of Support Vector Regression Models and Residual Analysis
Previous Article in Journal
Process-Driven and Flow-Based Processing of Industrial Sensor Data
Previous Article in Special Issue
A Novel Infrared Thermography Sensing Approach for Rapid, Quantitative Assessment of Damage in Aircraft Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Monitoring of Underground Structures in Multi-Layer Media by Dynamic Methods

1
Department of Information Systems in Construction, Faculty of IT-Systems and Technologies, Don State Technical University, Gagarin, 1, Rostov-on-Don 344000, Russia
2
Department of Transport Systems, Faculty of Roads and Transport Systems, Don State Technical University, Gagarin, 1, Rostov-on-Don 344000, Russia
3
Department of Life Safety and Environmental Protection, Faculty of Life Safety and Environmental Engineering, Don State Technical University, Gagarin, 1, Rostov-on-Don 344000, Russia
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(18), 5241; https://doi.org/10.3390/s20185241
Submission received: 31 July 2020 / Revised: 4 September 2020 / Accepted: 10 September 2020 / Published: 14 September 2020
(This article belongs to the Special Issue Sensors for NDT Diagnostics and Health Monitoring)

Abstract

:
The actual problem of structural monitoring and modeling of dynamic response from buried building is considered in the framework of arbitrary dynamic load. The results can be used for designing underground transport constructions, crossings, buried reservoirs and foundations. In existing methods, the system of sensors that register the response to a dynamic action does not allow for effective interpretation of the signal without understanding the dynamic features and resonance phenomena. The analytical and numerical solution of the problem of the dynamics of a buried object in a layered medium is considered. A multilayer half-space is a set of rigidly interconnected layers characterized by elastic properties. At a distance, an arbitrary dynamic load acts on the half-space, which causes oscillations in the embedded structure, and the sensor system registers the response. The problem of assessing the dynamic stress-strain state (DSSS) is solved using Fourier transforms with the principle of limiting absorption. As an example, the behavior of an embedded massive structure of an underground pedestrian crossing under the influence of a dynamic surface source on a multilayer medium is considered, as well as instrumental support of the sensor system. The solution in the form of stress, strain and displacement fields is obtained and compared with the experimental data. The frequency-dependent characteristics of the system are determined and the possibility of determining the DSSS by a shock pulse is shown.

1. Introduction

One of the critical tasks in the industry is diagnostics of the state and dynamic structural monitoring of constructions in a multilayer media to ensure their reliability and safe operation [1]. At the same time, the features of wave propagation in layered media, multiple reflections of disturbances from the boundaries of layers make such problems very difficult to interpret. Instrumental support of dynamic tests of buried objects is inevitably associated with the solution of dynamic problems of wave propagation and assessment of the stress-strain state of structures under various influences.
Problems of the dynamics of layered media have long attracted the attention of researchers. Many concepts of strength have changed with the advent of laminates. Damascus steels used to make weapons have excellent toughness and cutting properties due to their layered structure. Bulletproof vests are made of composite laminates and revolutionized personal firearm protection.
The geophysical structure of the Earth’s surface makes it possible to formulate the problems of the dynamics of buildings and structures lying on layered media under seismic influences. The dynamic behavior of buildings and structural elements under seismic influences are considered in the works of Terenzi et al. [2], Di Lorenzo et al. [3], and Sabino et al. [4].
At the same time, it is being improved as an instrumental base for digital measurements (Li Zhu et al. [5], Xiao-Wei Ye et al. [6]), and methods of signal processing and analysis. Anaissi et al. [7] presented a new algorithm for detecting and assessing damage in bridge structures. This technique uses tensor analysis to combine data and feature extraction, and additionally uses one-class support vector machines to detect structural damage. The indirect boundary element method was proposed by Ba et al. [8] to study three-dimensional (3D) dynamic responses of a two-dimensional (2D) hill in a layered elastic half-space subject to obliquely incident waves. The proposed method is applied to construct a scattered field in closed areas.
The spatial problem (3D) of the dynamic interaction of the soil structure was investigated by Hana et al. [9] using a model of a system with several degrees of freedom on a cylindrical rigid foundation embedded in a layered fluid-saturated poroelastic half-space using the indirect boundary element method (IBEM).
An efficient method for calculating layered composite materials was proposed in [10,11] based on the contact layer model.
The problems of non-destructive testing (NDT) and monitoring of the state of structures of engineering structures, such as oil and gas pipelines, rails, aircraft components, and adhesive joints, are considered in the work of Lugovtsova [12]. Modeling and analysis of wave propagation and their interaction with damage due to dispersion and multimodal nature was carried out on the basis of the finite element method with a scaled boundary. The results obtained made it possible to analyze the propagation of waves in multilayer structures and are important for the further development of non-destructive testing of engineering structures consisting of several layers.
The dynamics of layered structures is effectively used in the analysis of the stability of landslide slopes. In the works of Zhan et al. [13], the behavior of landslide slopes under the action of dynamic loads was considered. Experimental studies by Zou et al. [14] made it possible to simulate the behavior of bridge elements with different peak accelerations. Evaluation of the behavior of earthen dams was performed by Cai et al. [15]. The design and analysis of the destruction of sand structures was carried out by Zhou [16] to study their stability under dynamic loads. Using nonlinear and machine learning methods, i.e., grouped data processing (GMDH) and multivariate adaptive regression spline (MARS), two models were applied and developed to predict slope deformation in earth dams for different types of earthquakes. An analytical solution was presented by Guo [17] based on a limit equilibrium analysis.
The influence of geological factors of the coal seam on the initiation pressure and crack propagation under dynamic impacts was studied by Cao et al. [18]. As a result of theoretical analysis, a model was created for calculating pressure on a multilayer coal seam based on the fracture criterion in the theory of maximum tensile stress.
Analysis of the dynamic response of heterogeneous layered media with the aim of developing a method of NDT and determining their mechanical and strength characteristics was carried out in [19,20]. A comprehensive diagnostic method was proposed in [21] based on the analysis of the unsteady response of a multilayer structure, on the surface of which a shock load acts. During the solution procedure, various methods of genetic algorithms are used which have shown good results. The NDT method (Bertocci et al. [22]), based on scanning acoustic microscopy, for damage detection becomes very useful for detecting defects in multilayer structures with a thickness of several microns, providing a low time investment (limitation for other NDT methods) and quantitative analysis based on measurements.
Some problems of nonstationary dynamics are solved by direct methods using numerical analysis. Kuvat et al. [23] applied methods for determining the dynamic properties of mixtures of sand and bitumen, which can be used as damping materials in the foundations of building structures. In cyclic tests, the influence of the content and properties of bitumen on the ability to absorb energy was studied.
The design of buildings and structures under dynamic loads from moving vehicles is considered in the works of Zhou et al. [24] and Ai et al. [25]. Sun et al. [26] noted that high-performance roads require a reliable tool that can calculate the structural response caused by moving vehicles. Models based on the method of spectral elements are proposed for efficiently predicting the three-dimensional dynamic response of layered systems under a moving load. The study of the dynamic characteristics of a layered poroelastic soil lying under an uneven pavement exposed to traffic loads was carried out by Lyu et al. [27]. The pavement layer is modeled as a thin Kirchhoff plate, while the underlying base and substrate are treated as a saturated two-layer poroelastic medium. A numerical verification of the proposed models allows one to accurately predict the dynamic response of layered systems caused by a moving load.
Algorithms for artificial intelligence and creation of a neural network were developed in the articles by Beskopylny et al. [28] and Lyapin [29] based on direct algorithms for solving (Lyapin et al. [30]). The data of the response of the time domain of the multilevel structure to the diagnostic shock load are used as input information. This approach provides a tool for restoring the elastic properties of layers according to the corresponding thickness coordinates.
The above analysis of the literature shows that many problems of the dynamics of multilayer media are very relevant, are of great practical importance, and are solved either numerically or experimentally. At the same time, the dynamic features of wave reflections from the boundaries of media and the resonance properties of objects can be determined with great difficulty from the numerical or experimental analysis. In addition, in the process of designing a sensor system, frequency-dependent characteristics of the environment appear, leading to nonlinear responses and complicating data analysis. Thus, the purpose of this article is to construct an analytical method that allows simulating dynamic effects in layered media, as well as obtaining fundamental solutions, identifying resonance zones and features of the dynamic response of a structure using the example of underground structures.

2. Analytical Approach to the Problem of the Underground Object Dynamics in a Layered Medium

Let us consider the area occupied by a linear elastic (viscoelastic) medium, which is a multilayer half-space, under the surface of which an underground structure of arbitrary shape is located (Figure 1).
A multilayer half-space consists of a set of layers D = D 1 D 2 D N D 0 where D 1 = { x > 0 ; y ( , + ) } is occupied by the half-space. D j = { x ( x j , x j 1 ) ; y ( , + ) } , x j = i = 2 j h i ; is jth layer (j = 2, …, N); and D 0 is the area occupied by an underground structure.
The physical properties of the medium are described by the density ρ j and propagation velocities of transverse and longitudinal waves V S j , V P j .
The conditions for joining dissimilar layers are considered rigid with the requirement of continuity of the vectors of displacements and stresses when passing through the interface.
Without limiting generalization, rigid connection conditions can be changed to homogeneous frictionless contact conditions including continuity of normal displacements and absence of tangential stresses as well as Coulomb friction.
On the surface of the medium outside the structure in a limited area Ω , a system of forces q ( t ) oscillating with a frequency ω is set.
The wave fields for the displacements of the foundation points and the layered medium are described by integral representations (1) obtained on the basis of the dynamic reciprocity theorem
u ( 0 ) ( r 0 ) = γ T * ( r 0 , r ) n ( r ) u ( 0 ) ( r ) d s + Λ U * ( r 0 , r ) τ ( r ) d s + + Ω U * ( r 0 , r ) | x = 0 q ( y ) d y ;
u ( r 0 ) = γ T * * ( r 0 , r ) n ( r ) u ( r ) d s + Λ U * * ( r 0 , r ) τ ( r ) d s + + Ω U * * ( r 0 , r ) | x = 0 q ( y ) d y .
The components of the displacement vector u ( 0 ) ( r 0 ) at the boundary γ of the final body (structure) and the stress vector τ ( r ) in the area Λ of joining the structure with the layered foundation are unknown. Functions describing stresses and displacements at a point r for a homogeneous plane with a source r 0 are selected as fundamental solutions T * ( r 0 , r ) , U * ( r 0 , r ) for the region D 0 .
For a layered foundation, the matrices of fundamental solutions T * * ( r 0 , r ) , U * * ( r 0 , r ) are implemented using the superposition principle. This construction of a solution for a multilayer medium is based on the derivation of constitutive relations for one layer with stress vectors specified on its faces [26].
Let, in the local coordinate system for the jth layer ( x , y ) : x ( 0 , h j ) , y ( , ) , the amplitude functions of displacements under the action of a lumped source r 0 have the form
U ( * * ) ( r 0 , r ) = { U 1 ( j ) ( r 0 , r ) , U 2 ( j ) ( r 0 , r ) } T = { U x ( j ) ( r 0 , r ) , U y ( j ) ( r 0 , r ) } T
The functions U k ( j ) ( r 0 , r ) satisfying the equations of motion, according to the proposed method are sought in the form
U k ( j ) = U k ( j , 1 ) + U k ( j , 2 ) + U k ( j , 3 )
In a similar consideration of the half-plane, we assume U k ( j , 1 ) 0 .
Here, U k ( j , n ) ( x , y ) , n = 1 , 2 , the terms of this representation, are solutions of the Lamé equations for a homogeneous half-plane satisfying the boundary conditions
t ( j , 1 ) ( 0 , y ) = μ j X ( j , 1 ) ( y ) ,             t ( j , 2 ) ( h j , y ) = μ j X ( j , 2 ) ( y )
We represent the displacement vector U ( j , 1 ) ( x , y ) in the form of a Fourier integral through the transformants of the stress vector X ( j , 1 ) ( y )
U ( j , 1 ) ( x , y ) = 1 2 π Γ P ( j , 1 ) ( x , α ) X ˜ ( j , 1 ) ( α ) e i α y d α
The contour Γ is determined by the principle of limiting absorption: it bypasses the positive poles of the integrand from below, negative from above, and coincides with the real axis in the rest of the part.
Matrix components P ( j , 1 ) are given in Appendix A and depend on relative frequencies θ j 1 = ω V p j ; θ j 2 = ω V S j , where V p j , V s j are the wave speeds for corresponding medium.
To take into account the viscoelastic properties of the medium, the loss tangents are introduced for longitudinal and shear waves tg γ p ( j ) , tg γ s ( j ) , equal to ratio of imaginary parts of relative frequencies squares to real ones tg γ p ( j ) = Im ( θ j 1 2 ) Re ( θ j 1 2 ) , tg γ s ( j ) = Im ( θ j 2 2 ) Re ( θ j 2 2 ) . In this case, contour Γ equals the real axis completely.
Similar to formula (5), the displacements for the half-plane x h j through functions X ˜ ( j , 2 ) ( α ) are determined, where the elements P n m ( j , 2 ) ( x , α ) satisfy the relations
P n m ( j , 2 ) ( x , α ) = ( 1 ) δ n m P n m ( j , 1 ) ( h j x , α ) ,   n , m = 1 , 2
δ n m is the Kronecker symbol.
Determining the stress state of the layer as the sum of the corresponding solutions for two half-planes, we obtain
{ σ x , τ x y } ( j , 1 ) T ( x , y ) = 1 2 π Γ W ( j , 1 ) ( x , α ) X ˜ ( j , 1 ) ( α ) e i α y d α
where W ( j , 1 ) ( x , α ) are given in Appendix B.
For the second group of terms, we find
W n m ( j , 2 ) ( x , α ) = ( 1 ) δ n m + 1 W n m ( j , 1 ) ( h j x , α )
The functions U k ( j , 3 ) determine displacements in a homogeneous plane with the parameters of the layer under consideration from the action of a concentrated source of oscillations p ( r 0 ) = { p 1 ( x 0 , y 0 ) , p 2 ( x 0 , y 0 ) } in the form of a set of cylindrical waves and correspond to the components of the matrix U * ( r 0 , r )
U k ( j , 3 ) ( x 0 , y 0 , x , y ) = l = 1 2 U k l ( j , 3 ) ( x 0 , y 0 , x , y ) p l ( x 0 , y 0 ) ; k = 1 , 2
Using the decomposition formulas, they can be written in the Fourier transform form
U ˜ k l ( j , 3 ) ( x 0 , y 0 , x , α ) = F y [ U k l ( j , 3 ) ] = + U k l ( j , 3 ) ( x 0 , y 0 , x , y ) exp ( i α y ) d x U k l ( j , 3 ) ( x 0 , y 0 , x , y ) = F α 1 [ U ˜ k l ( j , 3 ) ] = 1 2 π Γ U ˜ k l ( j , 3 ) ( x 0 , y 0 , x , α ) exp ( i α y ) d x k , l = 1 , 2
The form of functions U ˜ k l ( j , 3 ) ( x 0 , y 0 , x , α ) is given in Appendix A.
Following the same logics for fundamental solutions for stresses, we get:
T k n ( j , 3 ) ( x 0 , y 0 , x , y ) = l = 1 2 T l k n ( j , 3 ) ( x 0 , y 0 , x , y ) p l ( x 0 , y 0 ) ; k , n = 1 , 2
T l k n ( j , 3 ) ( x 0 , y 0 , x , y ) in Fourier transformation form are given in Appendix B.
The introduced Fourier transforms of the stress functions X ˜ ( j , k ) ( α ) of representations (2), (7) are unknown and must be determined from the conditions for joining the dissimilar components of the layered half-plane with each other. Satisfying the equalities of the components of the vectors of displacements and stresses when passing through the interfaces between the media in the Fourier transforms, we obtain a system of linear algebraic equations 4 N + 2 with unknowns
A ( α ) X ˜ ( α ) = B ( α )
where X ˜ ( α ) is the general vector of unknown stresses for a multilayer structure.
In the frictionless case and under coupling conditions, the view of matrix A ( α ) will change. This leads to change in dispersion properties of construction and difference in properties for wave’s refraction and transition through partition boundary which is also presented in [31,32]. In particular, for road related constructions, the effects of smooth contact lead to decreasing in structural properties and physically to redistribution of vibration energy to higher frequency area especially in side zones around loading area. Such effect can be used for methods of road coatings quality inspection.
The fundamental solutions obtained in this way have the important property of the absence of stresses on the day surface x = x N . This excludes from participation in representations (1) unknown displacements on this part of the boundary of the region D .
Letting the point r 0 be further to the boundary of the area occupied by the foundation, we obtain a system of boundary integral equations for the vectors u ( 0 ) and τ
1 2 u ( 0 ) ( r 0 ) γ T * ( r 0 , r ) n ( r ) u ( 0 ) ( r ) d s Λ U * ( r 0 , r ) τ ( r ) d s = = Ω U * ( r 0 , r ) | x = 0 q ( y ) d y , r 0 γ ,
1 2 u ( 0 ) ( r 0 ) + γ T * * ( r 0 , r ) n ( r ) u ( 0 ) ( r ) d s Λ U * * ( r 0 , r ) τ ( r ) d s = = Ω U * * ( r 0 , r ) | x = 0 q ( y ) d y , r 0 Λ .
The system of Equations (13) is solved by the method of linear or constant boundary elements with subsequent spline approximation of the obtained distributions of unknown functions γ and Λ . The centers of boundary elements and docking nodes are chosen as collocation points.
Using the solution to system (13) and subsequent numerical differentiation, the authors consider the concentration of dynamic stresses at the boundary of the foundation depending on its depth and the ratio of the elastic parameters of the foundation layers.
As a computational example, Figure 2 represents σ x normal stresses distribution on lower edge of buried by 0.4 m in triple layered foundation of 0.5 m height. The soil foundation structure is supposed to be normal with increasing of wave speeds values along the layers thickness. The dynamic load is created by surface point sources, moved from foundation by 2.5 m and with the level of 100 kN. The plots represent cases of source oscillations with frequencies ω = 200 and 300 rad/s. The figure shows that stresses distribution can have areas of equal sign, which leads to bending and rotation of the foundation.

3. Experimental Methods

3.1. Experimental Technique

Experimental studies were carried out at an underground pedestrian crossing located in the city center of Rostov-on-Don, Russia. The total volume of the structure is up to 1200 m3, the depth of the footing is up to 3.8 m, and the depth of occurrence under the carriageway above the upper beams is 0.5 m. The view of the pedestrian crossing is shown in Figure 3.
Vibration measuring 12-channel device VK-12 (Figure 4) is developed and implemented on the basis of the E14–440 module manufactured by the company “L-Card” the level of noise (micro seismic vibrations) and transmission of signals in digital form for further processing in a computer (Table 1).
During experimental studies, the accelerometers sensors were fixed in the center of the ceilings using anchor bolts. This type of fastening provides the best match between the sensor and the overlap.
During experimental studies, the level of vibration velocities from the following sources was investigated:
(1)
passage of a KamAZ vehicle weighing 20.5 tons at a speed of 30 km/h;
(2)
passage of a KamAZ vehicle weighing 20.5 tons at a speed of 50 km/h;
(3)
impact on the floor.
Programs for signal processing are based on the transformation of data arrays, and the construction of amplitude–frequency characteristics (AFC) on the basis of the obtained AVX by applying the fast Fourier transform. For a more correct result, for each series of measurements, a correlation analysis was carried out, followed by averaging and smoothing of the results.
Amplitude–frequency characteristics obtained from ABX using fast Fourier transform reflect the integral values of the amplitudes for each of the frequencies, in a given interval of the original signal.
It should also be noted that, when choosing the method of vibrations excitation for evaluating the mechanical properties of surface layers and the structural state of the shallow laying object, it is possible to use effectively such diagnostic installations as falling weight of large mass HWD.

3.2. Results Processing Technique

In accordance with GOST R 52892–2007, the vibration velocities were obtained using the integral Fourier transform to the recorded amplitude–time characteristics of the vibration accelerations. Namely, considering the properties of the Fourier transform, as well as that A ( t ) = V ˙ ( t ) , the following dependence is valid
V ¯ ( ω ) = A ¯ ( ω ) i ω
Here, ω is the circular frequency of oscillations, i is the imaginary unit, A(t) and V(t) are the accelerations and velocities of oscillations, V ¯ ( ω ) , A ¯ ( ω ) are the Fourier images of the velocities and accelerations of oscillations, respectively, and
V ¯ ( ω ) = + V ( t ) exp ( i ω t ) d t
A ¯ ( ω ) = + A ( t ) exp ( i ω t ) d t
Thus, having registered vibration accelerations, the transition to vibration velocities was carried out by the formula
V ( t ) = 1 2 π + A ¯ ( ω ) i ω exp ( i ω t ) d ω
For the experimental determination of the frequency of natural bending vibrations of the beam, we used the results of recording the vibration accelerations under impact.
The shock impact corresponds to a Gaussian pulse [31]
p ( t ) = P exp ( a 2 t 2 )
Here, P is the peak value of the shock load and 2 / a is the duration of the pulse. The spectrum of a Gaussian pulse is determined by the formula
S ( ω ) = P + exp ( a 2 t 2 ) exp ( i ω t ) d t = P π a exp ( ( ω 2 a ) 2 )
where ω = 2 π f is the circular vibration frequency and f is the linear vibration frequency.
With a short exposure time, the spectrum of a Gaussian pulse is close to a constant over a wide frequency range.
Under external influence on the structure, the spectrum of the recorded response depends on the spectrum of external influence and the spectral characteristics of the structure. This dependence is expressed by the formula
S 0 = S B S K
where S0 is the response spectrum, SB is the spectrum of external influences, and SK is the structure’s own spectrum. Then, the spectrum of the structure is determined by the formula
S K = S 0 S B
Thus, it follows from formula (21) that under impact action the response spectrum will differ very little from the structure’s own spectrum.
Since the direct registration of the dynamic response of the structural elements of the transition to an arbitrary dynamic impact is difficult due to the large flow of pedestrians and vehicles, the calculated response spectra from the reference load can be used to assess the impact of an arbitrary non-stationary load on the structural elements. In this study, the impact was taken as a reference effect when a KamAZ vehicle with a mass of 20.5 t passed at speeds of 30 and 50 km/h.
To determine the dynamic response to an arbitrary dynamic action, it is necessary to do the following:
Perform test registration of the response to the dynamic impact under investigation in an arbitrary place (for example, on a road with low traffic volume).
Register the response to the reference action in the same place.
Get the design response by the formula
V ( t ) = 1 2 π + V ¯ N S ( ω ) V ¯ E S ( ω ) V ¯ E ( ω ) exp ( i ω t ) d ω
where V ¯ E S ( ω ) is the spectrum of response rates to the reference impact, V ¯ N S ( ω ) is the spectrum of response rates to the investigated non-stationary impact, and V ¯ E ( ω ) is the spectrum of response rates to the reference impact recorded on the structural elements of the pedestrian crossing.

4. Results and Discussion

As an example of the application of the proposed calculation method, let us consider the problem of determining the dynamic stresses on the inner boundary of a buried object caused by loads on the day surface of a multilayer medium. The problem is relevant when analyzing the strength and performance of underground structures. One of the approaches to its solution is to establish links between the maximum (peak) voltage values and the maximum values of vibration velocities recorded by sensors at measurement points. When considering the vibrations of beams (slabs) as elements of ground structures, usually GOST R 52892–2007 uses a ratio of the form
σ E ρ 3 G t o t G k n ν max
where, in accordance with GOST R 52892–2007, E is dynamic modulus of elasticity of the beam; ρ is the density of the beam material; G t o t is the total weight of the beam, including dead weight G and surface load; k n is modal coefficient corresponding to the fundamental mode of vibration, having a value in the range from 1 to 1.33; and v max is the maximum value of the vibration velocity along the entire length of the beam
However, the presence of a layered medium, into which a buried object is immersed, does not allow applying this formula directly to underground structures and requires taking into account the processes of wave propagation in inhomogeneous media.
Note also that, when waves propagate in homogeneous media, the relationship between the maximum values of stresses and velocities is preserved. Thus, for the simplest case of wave propagation from a surface lumped oscillating ω source of shear in the z direction of the axis, the displacements have the form
u z ( x , y ) = i 2 H 0 ( 1 ) ( θ x 2 + y 2 )
where θ = ω μ / ρ is the reduced vibration frequency, H 0 ( 1 ) is the Hankel function of the first kind, and μ is the shear modulus of the material. Hence, it is easy to show that, for all frequencies with distance r = x 2 + y 2 > > 1 from the source, σ z r v max μ ρ .
Thus, in the case of an object immersed in a layered medium, it can be established that
σ y E ρ C ( E i , ρ i , h i ) ν max
The dimensionless characteristic function C ( E i , ρ i , h i ) depends on the material properties of the components of the layered medium and the thicknesses of the layers h i , i = 1 , , N , primarily located above the object.
As examples of calculations (Figure 5), a problem is considered in which an object is located at a depth in Layer 3 and experiences dynamic effects q(t). Figure 6 and Figure 7 show the dependences C ( E i , ρ i , h i ) on the parameters of the problem.

4.1. The Underground Object Is Buried in the D3 Layer and the D4 Layer Is Missing

Consider the problem of the dynamics of an underground object located according to the scheme in Figure 5, when the D4 layer is absent. The properties of the media (parameters with index 0 refer to an underground facility) are as follows: E0 = 2 e4 MPa, E1 = 100 MPa, E2 = 80 MPa, ρ 0 = 2400 kg/m3, ρ 1 = 2000 kg/m3, and ρ 2 = 1900 kg/m3. Layer thicknesses are h 1 = , h 2 = 4 m, and h 3 = 4 m. Since there is no D4 layer, h4 = 0. Object dimensions are a = 6 m and b = 2.5 m. The wall thickness is 0.3 m. Depth of the object is h 0 = 1.2 m.
Figure 6 shows the results of calculating the characteristic function C ( E i , ρ i , h i ) depending on the elastic modulus of the layer in which the object is immersed. A wide range of elastic moduli is considered, from E = 1 × 105 Pa (that is, silty soils) to E = 1 × 1010 Pa (granite and harder). Figure 5 shows that soft layers and very hard layers increase the characteristic function C ( E i , ρ i , h i ) , and, in the zone of materials corresponding to sandy loams and loams, the characteristic function is close to 1. This makes it possible to effectively use dependence (25) in a wide range of experimental data.

4.2. Underground Object Buried in D3 Layer and D4 Layer Available

Medium properties are as follows: E0 = 2 e4 MPa, E1 = 100 MPa, E2 = 80 MPa, E3 = 200 MPa, ρ 0 = 2400 kg/m3, ρ 1 = 2000 kg/m3, ρ 2 = 1900 kg/m3, and ρ 3 = 2000 kg/m3. Layer thicknesses are h 1 = , h 2 = 4 m, h 3 = 4 m, and h 4 = 1 m. Object dimensions are a = 6 m and b = 2.5 m. Wall thickness is 0.3 m. Depth of the object is h 0 = 0.2 m.
Figure 7 shows the dependence of the characteristic function on the modulus of elasticity of the surface layer E4; the layer height is conventionally assumed to be 1 m.
Figure 7 shows that, in soft layers, when E4 of the surface layer D4 is less than the elastic modulus of the underlying layer E3, that is, E4 << E3, the characteristic function increases significantly, which is due to resonance phenomena in a layered medium. The upper layer resonates and makes a significant nonlinear contribution to the general picture of the dynamics of a layered medium. If the modulus of elasticity of the near-surface layer E4 is close to the modulus of the underlying layer E3, the characteristic function is close to 1. In the case when the rigidity of the surface layer E4 is greater than the underlying E3, that is, E4 >> E3, the harder layer dampens all resonances, being a kind of low-frequency filter in the system, and the characteristic function drops to 0.3.
Figure 8 shows the dependence of the characteristic function on the thickness of the surface layer. Asphalt concrete with a modulus of elasticity E4 = 1600 MPa was taken as a layer material.
Figure 8 shows that for the accepted value of the modulus of elasticity E4 = 1600 MPa in a wide range of layer thickness D4; the characteristic function is in the range from 0.8 to 1.2, which makes it possible to effectively apply dependence (25) in practical measurements of the stress–strain state.
Figure 9 shows the dependences of the amplitude–frequency characteristics of the layered medium. Here, 1 is a half-space; 2 is a softer near-surface layer E4 < E3; and 3 is a harder near-surface layer E4 > E3. P is the load, a is the coordinate of the point of deepening of the object, and E2 is the elastic modulus of the half-space. The dependence is plotted in dimensionless quantities.
Figure 9 shows that the softer near-surface layer E4 < E3 resonates and, due to multiple reflections from the layer boundaries, gives a change in the vibration amplitude. This effect differs from the results obtained by Pryakhina et al. [33], since, in our case, a massive object differs from a solid inclusion. A tougher layer works differently. Only in the low-frequency region there are still resonance frequencies corresponding to bending vibrations; in the rest of the range, resonances are damped by the upper hard layer. The results obtained are in good agreement with those of Zhou et al. [24], who obtained a solution for a cylindrical region in a layered medium.
Figure 10a,b shows the amplitude–time and amplitude–frequency dependences under shock loading.
Figure 10b shows that, under shock loading, the most energetic resonances lie in the low-frequency region and in the region of about 43 Hz, and then the frequency response fades. Note that the presence of low-frequency maxima of the amplitude–frequency characteristic is explained by the presence of a massive object immersed in a layered medium as a whole. In the case of significant rigidity of the underlying half-space, such maxima can be significant. In the limiting case of a rigidly clamped layer, the presence of isolated low-frequency resonances for a surface object was considered in the works of I.I. Vorovich [34,35].
Figure 11a,b shows the amplitude–time and amplitude–frequency dependences when a KamAZ car is driven at a speed of 50 km/h.
Figure 11 shows that the energy resonances lie in the low-frequency region. Resonance frequencies at 10, 17, and 19 Hz are due to the frequency of external influences when driving a car. Accordingly, the response of the structure changes its vibration spectrum.
Figure 12 shows an example of calculating the overlap spectrum using the amplitude–frequency characteristic of the vibration velocities of the overlap under impact and the Gaussian pulse spectrum. The figure clearly shows a peak around 42 Hz. It should be noted that, when passing vehicles on the amplitude–frequency characteristics, resonance peaks were not observed at frequencies above 20 Hz. It follows from this that the frequency range of the impact of vehicles lies below the natural frequencies of bending vibrations of floors.

5. Conclusions

The article discusses the problem and proposes a method for monitoring underground structures of arbitrary shape, lying in layered media, under various external dynamic influences. An analytical solution to the problem of the dynamics of layered media with buried objects of arbitrary shape is obtained. On the example of an underground passage, a practical example of the estimation of resonance frequencies and characteristic function for various characteristics and geometry of layers, is considered.
The instrumental analysis of the stress–strain state is complicated by the fact that at the characteristic function, which depends on the properties of the medium and has a nonlinear character. The frequency ranges in which the characteristic function is close to 1 are identified and recommendations are given for processing sensor signals for buried objects.
It is shown that the softer near-surface layer E4 < E3 resonates and, due to multiple reflections from the layer boundaries, gives a change in the vibration amplitude. A tougher layer works differently. Only in the low-frequency region there are still resonance frequencies corresponding to bending vibrations; in the rest of the range, resonances are damped by the upper hard layer.
As part of the research methodology development, it is proposed to use it in solving of inverse structural mechanics problems for reconstruction of buried object properties by measuring and analyzing vibration responses on a medium surface and an internal object geometry, as well as in the development of linear consolidation theory methods for layered water-saturated soils.

Author Contributions

Conceptualization, A.B., A.L., and B.M.; methodology, A.B., and A.L.; software, A.L.; validation, A.B., A.L., and B.M.; formal analysis, A.B., A.L.; investigation, A.B., A.L.; resources, A.B., A.L.; data curation, A.B., A.L., and B.M.; writing–original draft preparation, A.B.; writing–review and editing, A.B., A.L. and B.M.; visualization, A.B.; supervision, A.B., A.L., and B.M.; project administration, A.B., A.L., and B.M.; funding acquisition, B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Russian Foundation for Base Research (project 18-01-00715-a).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Components of displacements matrix:
P 11 ( j , 1 ) ( x , α ) = σ j 1 Δ R { ζ j 2 e j 1 + 2 α 2 e j 2 }
P 12 ( j , 1 ) ( x , α ) = i α Δ R { 2 σ j 1 σ j 2 e j 1 + ζ j 2 e j 2 }
P 21 ( j , 1 ) ( x , α ) = i α Δ R { ζ j 2 e j 1 + 2 σ j 1 σ j 2 e j 2 }
P 22 ( j , 1 ) ( x , α ) = σ j 2 Δ R { 2 α 2 e j 1 ζ j 2 e j 2 }
Δ R = 4 α 2 σ j 1 σ j 2 + ς j 4 ,   e j k = e σ j k x
ζ j 2 = α 2 + σ j 2 2 , σ j k 2 = α 2 θ j k 2
U ˜ 11 ( j , 3 ) ( x 0 , y 0 , x , α ) = A j [ σ j 1 E j 1 + α 2 E j 2 / σ j 2 ]
U ˜ 12 ( j , 3 ) ( x 0 , y 0 , x , α ) = U ˜ 21 ( j , 3 ) ( x 0 , y 0 , x , α ) = i α A j [ E j 1 E j 2 ] sign ( x x 0 )
U ˜ 22 ( j , 3 ) ( x 0 , y 0 , x , α ) = A j [ σ j 2 E j 2 + α 2 E j 1 / σ j 1 ]
E j k = exp ( σ j k | x 0 x | ) , k = 1 , 2 ; A j = exp ( i α y 0 ) 2 θ j 2 2

Appendix B

Components of stresses matrix:
W 11 ( j , 1 ) ( x , α ) = 1 Δ R { ζ j 4 e j 1 4 σ j 1 σ j 2 α 2 e j 2 }
W 12 ( j , 1 ) ( x , α ) = 2 i α σ j 2 ζ j 2 Δ R { e j 1 e j 2 } ,
W 21 ( j , 1 ) ( x , α ) = 2 i α σ j 1 ζ j 2 Δ R { e j 1 e j 2 }
W 22 ( j , 1 ) ( x , α ) = 1 Δ R { 4 σ j 1 σ j 2 α 2 e j 1 + ζ j 4 e j 2 }
T ˜ 111 ( j , 3 ) ( x 0 , y 0 , x , α ) = A j [ ς j 2 E j 1 2 α 2 E j 2 ] sign ( x x 0 )
T ˜ 211 ( j , 3 ) ( x 0 , y 0 , x , α ) = i α A j σ j 1 [ ς j 2 E j 1 2 σ j 1 σ j 2 E j 2 ]
( T ˜ 112 ( j , 3 ) ( x 0 , y 0 , x , α ) = i α A j σ j 2 [ 2 σ j 1 σ j 2 E j 1 ς j 2 E j 2 ]
T ˜ 212 ( j , 3 ) ( x 0 , y 0 , x , α ) = A j [ 2 α 2 E j 1 + ς j 2 E j 2 ] sign ( x x 0 )
T ˜ 122 ( j , 3 ) ( x 0 , y 0 , x , α ) = A j [ η j 2 E j 1 + 2 α 2 E j 2 ] sign ( x x 0 )
T ˜ 222 ( j , 3 ) ( x 0 , y 0 , x , α ) = i α A j σ j 1 [ ς j 2 E j 1 + σ j 2 ( 2 θ j 2 2 θ j 1 2 ( 1 σ j 2 σ j 1 ) ) E j 2 ]
η j 2 = 2 α 2 + λ j μ j θ j 1 2

References

  1. Beskopylny, A.; Lyapin, A.; Andreev, V. Layered structures mechanical properties assessment by dynamic tests. MATEC Web Conf. 2017, 117, 00018. [Google Scholar] [CrossRef]
  2. Terenzi, G.; Bazzani, C.; Costoli, I.; Sorace, S.; Spinelli, P. Advanced Seismic Retrofit of a Mixed R/C-Steel Structure. Buildings 2019, 9, 241. [Google Scholar] [CrossRef] [Green Version]
  3. Di Lorenzo, G.; De Martino, A. Earthquake Response of Cold-Formed Steel-Based Building Systems: An Overview of the Current State of the Art. Buildings 2019, 9, 228. [Google Scholar] [CrossRef] [Green Version]
  4. Sabino, A.; Mannella, A.; de Leo, A.M. Seismic Response of a Structure Equipped with an External Viscous Damping System. Buildings 2020, 10, 19. [Google Scholar] [CrossRef] [Green Version]
  5. Zhu, L.; Fu, Y.; Chow, R.; Spencer, B.F.; Park, J.W.; Mechitov, K. Development of a High-Sensitivity Wireless Accelerometer for Structural Health Monitoring. Sensors 2018, 18, 262. [Google Scholar] [CrossRef]
  6. Ye, X.W.; Su, Y.H.; Xi, P.S. Statistical Analysis of Stress Signals from Bridge Monitoring by FBG System. Sensors 2018, 18, 491. [Google Scholar] [CrossRef] [Green Version]
  7. Anaissi, A.; Makki Alamdari, M.; Rakotoarivelo, T.; Khoa, N.L.D. A Tensor-Based Structural Damage Identification and Severity Assessment. Sensors 2018, 18, 111. [Google Scholar] [CrossRef] [Green Version]
  8. Ba, Z.; Wang, Y.; Liang, J.; Liu, X. 3D dynamic responses of a 2D hill in a layered half-space subjected to obliquely incident plane P-, SV- and SH-waves. Eng. Anal. Bound. Elem. 2019, 105, 129–145. [Google Scholar] [CrossRef]
  9. Hana, B.; Lianga, J.; Fub, J.; Liua, R. 3D dynamic soil-structure interaction in layered, fluid-saturated, poroelastic half-space. Soil Dyn. Earthq. Eng. 2019, 120, 113–126. [Google Scholar] [CrossRef]
  10. Andreev, V.; Tsybin, N.; Turusov, R. Layered composite and contact layer. Effective modulus of elasticity. E3 S Web Conf. 2019, 97, 04071. [Google Scholar] [CrossRef]
  11. Andreev, V.; Turusov, R.; Tsybin, N. Layered composite and contact layer. Normal separation and transversal strength. MATEC Web Conf. 2018, 251, 04066. [Google Scholar] [CrossRef]
  12. Lugovtsova, Y.; Bulling, J.; Boller, C.; Prager, J. Analysis of Guided Wave Propagation in a Multi-Layered Structure in View of Structural Health Monitoring. Appl. Sci. 2019, 9, 4600. [Google Scholar] [CrossRef] [Green Version]
  13. Zhan, Z.; Qi, S. Numerical Study on Dynamic Response of a Horizontal Layered-Structure Rock Slope under a Normally Incident Sv Wave. Appl. Sci. 2017, 7, 716. [Google Scholar] [CrossRef] [Green Version]
  14. Zou, Z.; Lei, D.; Jiang, G.; Luo, B.; Chang, S.; Hou, C. Experimental Study of Bridge Foundation Reinforced with Front and Back Rows of Anti-Slide Piles on Gravel Soil Slope under El Centro Waves. Appl. Sci. 2020, 10, 3108. [Google Scholar] [CrossRef]
  15. Cai, M.; Koopialipoor, M.; Armaghani, D.J.; Thai Pham, B. Evaluating Slope Deformation of Earth Dams Due to Earthquake Shaking Using MARS and GMDH Techniques. Appl. Sci. 2020, 10, 1486. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, Z.; Lei, J.; Shi, S.; Liu, T. Seismic Response of Aeolian Sand High Embankment Slopes in Shaking Table Tests. Appl. Sci. 2019, 9, 1677. [Google Scholar] [CrossRef] [Green Version]
  17. Guo, S.; Qi, S.; Yang, G.; Zhang, S.; Saroglou, C. An Analytical Solution for Block Toppling Failure of Rock Slopes during an Earthquake. Appl. Sci. 2017, 7, 1008. [Google Scholar] [CrossRef]
  18. Cao, S.; Li, X.; Zhou, Z.; Wang, Y.; Ding, H. Investigation of Hydraulic Fracturing Crack Propagation Behavior in Multi-Layered Coal Seams. Appl. Sci. 2020, 10, 1153. [Google Scholar] [CrossRef] [Green Version]
  19. Beskopylny, A.; Onishkov, N.; Korotkin, V. Assessment of the Fatigue Durability of the Rolling Contact. In Proceedings of the International Scientific Conference Energy Management of Municipal Transportation Facilities and Transport EMMFT 2017, Khabarovsk, Russia, 10–13 April 2017; pp. 184–191. [Google Scholar] [CrossRef]
  20. Beskopylny, A.; Kadomtseva, E.; Strelnikov, G.; Morgun, L.; Berdnik, Y.; Morgun, V. Model of heterogeneous reinforced fiber foam concrete in bending. IOP Conf. Ser. Mater. Sci. Eng. 2018, 365, 032023. [Google Scholar] [CrossRef]
  21. Beskopylny, A.; Lyapin, A.; Kadomtsev, M.; Veremeenko, A. Complex method of defects diagnostics in underground structures. MATEC Web Conf. 2018, 146, 02013. [Google Scholar] [CrossRef] [Green Version]
  22. Bertocci, F.; Grandoni, A.; Djuric-Rissner, T. Scanning Acoustic Microscopy (SAM): A Robust Method for Defect Detection during the Manufacturing Process of Ultrasound Probes for Medical Imaging. Sensors 2019, 19, 4868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kuvat, A.; Sadoglu, E. Dynamic properties of sand-bitumen mixtures as a geotechnical seismic isolation material. Soil Dyn. Earthq. Eng. 2020, 132, 106043. [Google Scholar] [CrossRef]
  24. Zhou, S.; He, C.; Guo, P.; Di, H.; Xiao, J. Three-Dimensional Analytical Model for Coupled Track-Tunnel-Soil System in a Multilayered Half-Space. J. Eng. Mech. 2020, 146, 04019121. [Google Scholar] [CrossRef]
  25. Ai, Z.; Ren, G. Dynamic analysis of a transversely isotropic multilayered half-plane subjected to a moving load. Soil Dyn. Earthq. Eng. 2016, 83, 162–166. [Google Scholar] [CrossRef]
  26. Sun, Z.; Kasbergen, C.; Skarpas, A.; Anupam, K.; vanDalen, K.; Erkens, S. Dynamic analysis of layered systems under a moving harmonic rectangular load based on the spectral element method. Int. J. Solids Struct. 2019, 180–181, 45–61. [Google Scholar] [CrossRef]
  27. Lyu, Z.; Qian, J.; Shi, Z.; Gao, Q. Dynamic responses of layered poroelastic ground under moving traffic loads considering effects of pavement roughness. Soil Dyn. Earthq. Eng. 2020, 130, 105996. [Google Scholar] [CrossRef]
  28. Beskopylny, A.; Lyapin, A.; Anysz, H.; Meskhi, B.; Veremeenko, A.; Mozgovoy, A. Artificial Neural Networks in Classification of Steel Grades Based on Non-Destructive Tests. Materials 2020, 13, 2445. [Google Scholar] [CrossRef]
  29. Lyapin, A. On the layered structures material properties definition. IOP Conf. Ser. Mater. Sci. Eng. 2019, 698, 066015. [Google Scholar] [CrossRef]
  30. Shatilov, Y.; Lyapin, A. Vibration-Based Damage Detection Techniques for Health Monitoring of Construction of a Multi-Storey Building. Mater. Sci. Forum 2018, 931, 178–183. [Google Scholar] [CrossRef]
  31. Tiraturyan, A.N.; Uglova, E.V.; Lyapin, A.A. Studying the energy distribution of the dynamic influences of road transport on the layers of nonrigid pavements. PNRPU Mech. Bull. 2017, 2, 178–194. [Google Scholar] [CrossRef]
  32. Liu, J. A dynamic modelling method of a rotor-roller bearing-housing system with a localized fault including the additional excitation zone. J. Sound Vib. 2020, 469, 115144. [Google Scholar] [CrossRef]
  33. Pryakhina, O.D.; Smirnova, A.V. On the study of waveguide properties of a stack of elastic layers with a set of rigid inclusions. Mech. Solids 2009, 44, 380–388. [Google Scholar] [CrossRef]
  34. Vorovich, I.I.; Pryakhina, O.D.; Tukodova, O.M.; Frejgejt, M.R. An approach to solving the dynamic problems for laminated electroelastic and anisotropic media. J. Appl. Math. Mech. 1995, 59, 625–634. [Google Scholar] [CrossRef]
  35. Vorovich, I.I.; Babeshko, V.A.; Pryakhina, O.D. Dynamics of Massive Bodies and Resonance Phenomena in Deformable Media; Nauch. Mir: Moscow, Russian, 1999. [Google Scholar]
Figure 1. Layered half-space scheme with an arbitrary underground structure.
Figure 1. Layered half-space scheme with an arbitrary underground structure.
Sensors 20 05241 g001
Figure 2. Normal stresses on lower edge of foundation.
Figure 2. Normal stresses on lower edge of foundation.
Sensors 20 05241 g002
Figure 3. Underground pedestrian crossing.
Figure 3. Underground pedestrian crossing.
Sensors 20 05241 g003
Figure 4. Vibration measuring device (a); and scheme of the sensors installation (b).
Figure 4. Vibration measuring device (a); and scheme of the sensors installation (b).
Sensors 20 05241 g004aSensors 20 05241 g004b
Figure 5. Model for calculating the dynamic response of an object D0 in a multilayer medium D1, D2, D3 and D4.
Figure 5. Model for calculating the dynamic response of an object D0 in a multilayer medium D1, D2, D3 and D4.
Sensors 20 05241 g005
Figure 6. Dependence C ( E i , ρ i , h i ) on the modulus of elasticity of the immersion layer E3, Pa.
Figure 6. Dependence C ( E i , ρ i , h i ) on the modulus of elasticity of the immersion layer E3, Pa.
Sensors 20 05241 g006
Figure 7. Dependence C ( E i , ρ i , h i ) on the modulus of elasticity of the immersion layer E3, Pa, h 4 = 1 .
Figure 7. Dependence C ( E i , ρ i , h i ) on the modulus of elasticity of the immersion layer E3, Pa, h 4 = 1 .
Sensors 20 05241 g007
Figure 8. Dependence C ( E i , ρ i , h i ) on thickness h4 of the surface layer E4 at E 4 = 1600 MPa (asphalt concrete).
Figure 8. Dependence C ( E i , ρ i , h i ) on thickness h4 of the surface layer E4 at E 4 = 1600 MPa (asphalt concrete).
Sensors 20 05241 g008
Figure 9. Dependence of the amplitude–frequency characteristics of the layered medium: 1, half-space; 2, softer near-surface layer E4 < E3; 3, harder near-surface layer E4 > E3.
Figure 9. Dependence of the amplitude–frequency characteristics of the layered medium: 1, half-space; 2, softer near-surface layer E4 < E3; 3, harder near-surface layer E4 > E3.
Sensors 20 05241 g009
Figure 10. Amplitude–time (a) and amplitude–frequency (b) characteristics of vibration velocities, impact load.
Figure 10. Amplitude–time (a) and amplitude–frequency (b) characteristics of vibration velocities, impact load.
Sensors 20 05241 g010
Figure 11. Amplitude–time (a) and amplitude–frequency (b) characteristics of vibration velocities, passage of a Kamaz car at a speed of 50 km/h.
Figure 11. Amplitude–time (a) and amplitude–frequency (b) characteristics of vibration velocities, passage of a Kamaz car at a speed of 50 km/h.
Sensors 20 05241 g011
Figure 12. Response spectra of underpass overlap: 1, theoretical; 2, experimental.
Figure 12. Response spectra of underpass overlap: 1, theoretical; 2, experimental.
Sensors 20 05241 g012
Table 1. Main characteristics of the vibration measuring device.
Table 1. Main characteristics of the vibration measuring device.
1Dynamic range of input levels50 db
2Bit depth analog-to-digital Converter12 bit
3Power consumption of the IV-1 device <250 mA
4Complex measurement error<5%
5Operating temperature range0 °C … +55 °C
6Conversion rate>2500 mW/m−1 s2

Share and Cite

MDPI and ACS Style

Lyapin, A.; Beskopylny, A.; Meskhi, B. Structural Monitoring of Underground Structures in Multi-Layer Media by Dynamic Methods. Sensors 2020, 20, 5241. https://doi.org/10.3390/s20185241

AMA Style

Lyapin A, Beskopylny A, Meskhi B. Structural Monitoring of Underground Structures in Multi-Layer Media by Dynamic Methods. Sensors. 2020; 20(18):5241. https://doi.org/10.3390/s20185241

Chicago/Turabian Style

Lyapin, Alexandr, Alexey Beskopylny, and Besarion Meskhi. 2020. "Structural Monitoring of Underground Structures in Multi-Layer Media by Dynamic Methods" Sensors 20, no. 18: 5241. https://doi.org/10.3390/s20185241

APA Style

Lyapin, A., Beskopylny, A., & Meskhi, B. (2020). Structural Monitoring of Underground Structures in Multi-Layer Media by Dynamic Methods. Sensors, 20(18), 5241. https://doi.org/10.3390/s20185241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop