Next Article in Journal
Uniqueness Criteria for Solving a Time Nonlocal Problem for a High-Order Differential Operator Equation l(·)—A with a Wave Operator with Displacement
Next Article in Special Issue
New Tribological Aspects in the Micro-Areas of the Symmetric Rolling-Sliding Contact
Previous Article in Journal
Using K-Means Clustering in Python with Periodic Boundary Conditions
Previous Article in Special Issue
On the Vibrations of a Rigid Solid Hung by Kinematic Chains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Response of Vibratory Piling Machines for Ground Foundations

by
Adrian Mihai Goanță
1,
Polidor Bratu
1,2 and
Nicușor Drăgan
1,*
1
Research Center for Mechanics of Machines and Technological Equipment—MECMET, Faculty of Engineering and Agronomy in Braila, “Dunărea de Jos” University of Galati, 800008 Galati, Romania
2
The Institute of Solid Mechanics of the Romanian Academy, 021652 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Submission received: 6 May 2022 / Revised: 10 June 2022 / Accepted: 12 June 2022 / Published: 14 June 2022
(This article belongs to the Special Issue Dynamic Systems and Mechanics)

Abstract

:
Vibrating technological equipment for the introduction of piles and columns into the ground of construction foundations (named vibratory piling machines) is crucial in the process of building stable and resilient foundations for civil engineering, hydrotechnical construction, special construction (e.g., military constructions), bridges, roads and industrial platforms. During the works carried out by the construction companies in various geographical areas of Romania, particularities of the dynamic technological regimes influenced by the nature of the land were identified at the deep introduction of the construction elements in the form of piles or circular (tubular) columns. The results of applied research, rheological modeling and optimization of vibrating equipment, highlight the need for an analytical approach that takes into account the parametric variations of the elastic and damping characteristics of some categories of soils on the depth of piles or foundation columns. In this context, the paper presents the calculation model with the dynamic response for the vibrating equipment of insertion with disturbing forces of 200–1250 kN for piles or columns with lengths of 10–30 m. The novelty of the research study consists in the linear rheological model, which was adopted in the form of a Maxwell–Voigt–Kelvin schematic of the type (E-V)–(E|V), with a discrete variation in four values for stiffness and damping of the soil, as the piles or columns vibrate and advance in the ground foundation. Practical experience of the authors in the field of using vibrogenerators for the introduction of piles in various types of ground foundations led to the adoption of the rheological model with variable damping coefficients depending on the depth of penetration into the soil. The curves of the dissipated power confirm the experimental data obtained in situ, in accordance with the rheological indoor tests of the different types of soil foundations.

1. Introduction

The capability level of the vibrating machines in this specialized class is evaluated according to the size of the functional parameters in the working mode, namely: vibration amplitudes, transmitted force to the soil (ground foundation), dynamic insulation in the machine and dissipated energy [1,2].
In Romania, vibrating machines with dynamic action are used, whose operating regime is correlated with the nature of the soil [3,4]. Thus, several distinct classes of foundation soils have been identified, consisting mainly of the following categories: loess-sensitive soils, contractile soils (with volume variations), fine and uniform sands with the possibility of liquefaction, clayey and dusty soils and sandy clay soils [5,6].
An equivalent linear calculation model, which takes into account a complex rheological behavior, is adopted in this paper. Also, for theoretical reasons, numerical analysis is performed using a linear visco-elastic approach [7,8].
According to the Federal Highway Administration (FHWA), the relationships between maximum soil particle velocity and pile peak distance are based on the nominal energy of the vibrogenerator for two types of foundation soil: Class II (soils with standard penetration test values 5–15 strokes per 0.3 m) and Class III (hard soils with standard penetration test values of 15–50 strokes per 0.3 m). In general, there is an increase in the kinematic parameters of the cell (displacements, speeds, accelerations) as the rigidity of the ground and the energy of the vibrogenerator increase [9]. The classification of piles can be made according to the method of installation, through which piles are usually classified into two categories: non-displaced piles and displaced piles. The first category of piles requires the prior removal of the ground, while the trained prefabricated piles (the second category) are characterized by the appearance of ground movements to allow the piles to enter the ground, without the need to remove the ground beforehand [10].
The interaction between piles and foundation soils was studied for the past five decades [11]. The scientific and technical literature presents various models of the interaction between the foundation soil and the prefabricated piles introduced by shocks and/or vibrations; analysis and calculation models are of three types: empirical/experimental models, theoretical models and engineering models [12,13,14,15].
The dynamic interaction between the pile and the foundation soil was not sufficiently studied due to the complexity of the vibrogenerator–pile–foundation soil model [16]. The studied dynamic models of the interactions of vibrogenerator–pile–foundation soil are theoretical models [17,18,19], numerical simulations models [20] and experimental models (in situ or indoor tests) [21,22,23]. Numerical simulation and, especially, analytical modeling of the movement of vibrating piles in saturated soil is quite difficult [24]. Although open-ended piles (circular or rectangular tube type) are common in engineering practice, studies focused mainly on closed-ended piles [25].
Experimental tests are often used to establish numerical simulation models and to validate theoretical models. Some experimental data on the vibration introduction of prefabricated piles into the foundation field are available in the literature. The experimental tests, in the laboratory and in situ, took into account the two types of piles: experiments with closed-ended piles and experiments with open-ended piles [26,27].
For vibratory machines in this technological class, the Maxwell–Voigt–Kelvin rheological model, denoted (E-V)–(E|V), with dynamic force generators (vibrating motors), was realized, so that between the physical–mechanical and terrain-modeling parameters, the rigidity k and the damping coefficient c, the amplitude of the steady-state forced vibrations A and the transmitted force to the ground Q could be determined by an analytical calculation relation [28,29,30].
In this way, the calculation relations and the families of curves were determined, which eloquently reflect the variations of the dynamic parameters for the technological regimes defined by the pulsation ω and the amplitude A of the technological vibrations [31,32].
This paper is organized as follows: Section 2 presents a linear dynamic model of the vibratory piling machine in interaction with the soil foundation; the dynamic analysis, in steady-state forced vibration condition of this model, consists in determination of the amplitude displacements, dynamic forces transmitted to the ground, dynamic transmissibility coefficient and dissipated energy per harmonic cycle function of damping characteristic of the soil and the functional parameters of the vibratory piling machine. Section 3 presents a case study of a vibratory piling machine. It shows graphs of the variation in the dynamic and energetic parameters of the vibrogenerator–pile–foundation soil system considered as an (E-V)–(E|V) rheological model, function of damping ratio ζ and the relative pulsation Ω . The plotted parameters of the systems are the amplitudes of the displacements of the mobile mass of the vibrating machine, the amplitudes of the displacements and of the deflections of the components of the rheological model, the amplitudes of the dynamic transmitted force to the ground, the dynamic transmissibility coefficient and the dissipated energy per cycle. Section 4 emphasizes the importance of technological parameter Ω on the analysis of the dynamic characteristics of the vibrogenerator–pile–foundation soil system. The final concluding remarks are presented in Section 5.

2. Linear Dynamic Model of the Vibratory Piling Machine

The linear dynamical system with unidirectional vertical motion is shown in Figure 1 [32,33]. The serial structure of the system consists of mass m, the rheological Maxwell model type (E-V) and the rheological model Voigt–Kelvin type (E|V). The dynamic model is disturbed by the vertical harmonic force F ( t ) = F 0 sin ω t , the coordinates being the vertical displacements: x = x ( t ) of the mass m, y = y ( t ) in the Voigt–Kelvin model and z = z ( t ) of the point of connection of the elements of the Maxwell model.
.
The rheological characteristics (stiffness, damping) of the linear models are ( k ,   c ) for the Voigt–Kelvin (E|V) model and ( N k ,   S c ) for the Maxwell (E-V) model, where S and N are positive multiplier numbers and the relation between the two is S = λ N , where λ is also a positive multiplier number.

2.1. Dynamic Response of the Linear Rheological System

For the dynamic system shown in Figure 1, the governing equations (differential equations of motion), expressed in complex numbers, can be written as follows [34,35],
{ m x ˜ ¨ + c S ( x ˜ ˙ z ˜ ˙ ) + k x ˜ = F 0 sin ω t k y ˜ + c y ˜ ˙ k N ( z ˜ y ˜ ) = 0 k N ( z ˜ y ˜ ) = c S ( x ˜ ˙ z ˜ ˙ )
The steady-state forced vibration responses, in terms of displacements, velocities and accelerations, are
x ˜ = X ˜ e j ω t   d d t   x ˜ ˙ = j ω X ˜ e j ω t   d d t   x ˜ ¨ = ω 2 X ˜ e j ω t  
y ˜ = Y ˜ e j ω t   d d t   y ˜ ˙ = j ω Y ˜ e j ω t   d d t   y ˜ ¨ = ω 2 Y ˜ e j ω t
z ˜ = Z ˜ e j ω t   d d t   z ˜ ˙ = j ω Z ˜ e j ω t   d d t   z ˜ ¨ = ω 2 Z ˜ e j ω t   ,
which, replaced in (1), lead to the linear algebraic system in complex quantities X ˜ ,   Y ˜ ,   Z ˜ in the form
{ ( m ω 2 + j c ω S ) X ˜ j c ω S Z ˜ = F 0 [ k ( 1 + N ) + j c ω ] Y ˜ k N Z ˜ = 0 j c ω S X ˜ k N Y ˜ + [ k N + j c ω S ] Z ˜ = 0
The complex determinant D ˜ of the linear algebraic system matrix is obtained as follows,
D ˜ = d e t [ m ω 2 + j c ω S 0 j c ω S 0 k ( 1 + N ) + j c ω k N j c ω S k N k N + j c ω S ]   ,
from where
D ˜ = ( m ω 2 + j c ω S ) { [ k ( 1 + N ) + j c ω ] [ k N + j c ω S ] k 2 N 2 } j c ω S { j c ω S [ k ( 1 + N ) + j c ω ] }
Using the notations α = m ω 2 , β = c ω , δ = k N , γ = k ( 1 + N ) and S = λ N , the expression of the determinant D ˜ becomes
D ˜ = [ α δ ( γ δ ) β 2 S ( α + δ ) ] + j β [ α ( γ S + δ ) + δ S ( γ δ ) ]
Returning to the structural parameters ( m , c , k , ω ) of the system, the determinant D ˜ can be written function of the variables ( c , ω ) as follows,
D ˜ ( c , ω ) = [ m ω 2 k 2 N + c 2 ω 2 λ N m ω 2 c 2 ω 2 k λ N 2 ]                                      + j c ω [ m ω 2 k ( N + λ N + λ N 2 ) + k 2 λ N 2 ]
As a function of the dimensionless parameters of the system, ζ = c 2 m k —damping ratio and Ω = ω ω n —relative pulsation, the determinant D ˜ can be written
D ˜ ( ζ ,   Ω ) = N k 3 [ Ω 2 ( 2 ζ Ω ) 2 λ N + Ω 2 ( 2 ζ Ω ) 2 λ ] + j ( 2 ζ Ω ) N k 3 [ λ N Ω 2 ( 1 + λ + λ N ) ]
or
D ˜ ( ζ ,   Ω ) = N k 3 [ C + j D ]
The complex determinant D ˜ X of unknown variable X ˜ can be calculated as follows,
D ˜ X = d e t [ F 0 0 j c ω S 0 k ( 1 + N ) + j c ω k N 0 k N k N + j c ω S ] ,
or
D ˜ X ( c , ω ) = F 0 [ k 2 N c 2 ω 2 λ N ] + j c ω k N F 0 [ 1 + λ + λ N ]
As a function of the dimensionless parameters ( ζ ,   Ω ) , the determinant D ˜ X can be written as follows,
D ˜ X ( ζ ,   Ω ) = F 0 k 2 N { [ 1 ( 2 ζ Ω ) 2 λ ] + j [ ( 2 ζ Ω ) ( 1 + λ + λ N ) ] } ,
or
D ˜ X ( ζ ,   Ω ) = F 0 k 2 N [ A + j B ]
The complex determinant D ˜ Y of unknown variable Y ˜ can be calculated as
D ˜ Y ( c , ω ) = d e t [ m ω 2 + j c ω S F 0 j c ω S 0 0 k N j c ω S 0 k N + j c ω S ] = j F 0 c ω k λ N 2
Or the function of variables ( ζ ,   Ω ) :
D ˜ Y ( ζ ,   Ω ) = j F 0 k 2 λ N 2 ( 2 ζ Ω )
With the notation E = ( 2 ζ Ω ) , the determinant D ˜ Y becomes
D ˜ Y ( ζ ,   Ω ) = F 0 k 2 λ N 2 [ j E ]
The complex determinant D ˜ Z of unknown variable Z ˜ can be calculated as follows,
D ˜ Z = d e t [ m ω 2 + j c ω S 0 F 0 0 k ( 1 + N ) + j c ω 0 j c ω S k N 0 ] ,
from where
D ˜ Z ( c , ω ) = F 0 c ω λ N [ c ω + j k ( 1 + N ) ]
Function of the dimensionless parameters ( ζ ,   Ω ) , the determinant D ˜ Z can be written
D ˜ Z ( ζ ,   Ω ) = F 0 k 2 λ N ( 2 ζ Ω ) [ ( 2 ζ Ω ) + j ( 1 + N ) ]
or, using the notations a = 2 ζ Ω and b = 1 + N :
D ˜ Z ( ζ ,   Ω ) = F 0 k 2 λ N ( 2 ζ Ω ) [ a + j b ]

2.1.1. Forced Vibration: Amplitude Displacement x ( t ) in Steady-State Condition

The harmonic dynamic response in complex displacement X ˜ is
X ˜ = D ˜ X ( ζ ,   Ω ) D ˜ ( ζ ,   Ω ) = F 0 k 2 N [ A + j B ] N k 3 [ C + j D ] = F 0 k A + j B C + j D
where it uses the notations
{ A = 1 ( 2 ζ Ω ) 2 λ ;   B = ( 2 ζ Ω ) ( 1 + λ + λ N ) C = Ω 2 ( 2 ζ Ω ) 2 λ N + Ω 2 ( 2 ζ Ω ) 2 λ                           D = ( 2 ζ Ω ) [ λ N Ω 2 ( 1 + λ + λ N ) ]                                
From Equation (13) the complex displacement X ˜ is obtained as follows,
X ˜ = F 0 k 1 C 2 + D 2 [ ( A C + B D ) + j ( B C A D ) ]
with the amplitude of steady-state forced vibration | X ˜ | = X 0 obtained by
X 0 = F 0 k A 2 + B 2 C 2 + D 2
As a function of the dimensionless parameters ( ζ ,   Ω ) , the amplitude X 0 can be written as follows,
X 0 ( ζ ,   Ω ) = F 0 k [ 1 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 ( 1 + λ + λ N ) 2 [ Ω 2 + ( 2 ζ Ω ) 2 λ N Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2

2.1.2. Forced Vibration: Amplitude Displacement y ( t ) in Steady-State Condition

The harmonic dynamic response in complex displacement Y ˜ is
Y ˜ = D ˜ Y ( ζ ,   Ω ) D ˜ ( ζ ,   Ω ) = F 0 k λ N j E C + j D = F 0 k λ N 1 C 2 + D 2 ( D E j C E )
where E, C and D are obtained by Equations (10a) and (13a).
The amplitude of steady-state forced vibration | Y ˜ | = Y 0 can be calculated as follows,
Y 0 = F 0 k λ N E C 2 + D 2
Using the expressions from (10a) and (13a) for the terms E, C and D, the amplitude Y 0 can be written
Y 0 ( ζ ,   Ω ) = F 0 k λ N 2 ζ Ω C 2 + D 2   ,
or
Y 0 ( ζ ,   Ω ) = F 0 k λ N 2 ζ Ω [ Ω 2 + ( 2 ζ Ω ) 2 λ N Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2

2.1.3. Forced Vibration: Amplitude Displacement z ( t ) in Steady-State Condition

Using the notations from Equations (12a) and (13a), the harmonic dynamic response in complex displacement Z ˜ is
Z ˜ = D ˜ Z ( ζ ,   Ω ) D ˜ ( ζ ,   Ω ) = F 0 k λ ( 2 ζ Ω ) a + j b C + j D   ,  
or
Z ˜ = F 0 k λ ( 2 ζ Ω ) 1 C 2 + D 2 [ ( a C + b D ) + j ( b C a D ) ]
The amplitude of steady-state forced vibration | Z ˜ | = Z 0 can be calculated as follows,
Z 0 = F 0 k λ ( 2 ζ Ω ) a 2 + b 2 C 2 + D 2   ,
from where
Z 0 = F 0 k λ ( 2 ζ Ω ) ( 2 ζ Ω ) 2 + ( 1 + N ) 2 C 2 + D 2   ,
or
Z 0 = F 0 k λ ( 2 ζ Ω ) ( 2 ζ Ω ) 2 + ( 1 + N ) 2 [ Ω 2 + ( 2 ζ Ω ) 2 λ N Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2

2.2. Dynamic Force Transmitted to the Ground Foundation

The dynamic force Q ˜ is transmitted to the ground foundation through the spring k and the damper c in the Voigt–Kelvin rheological model, thus, the force Q ˜ can be written as a sum of elastic force of the spring and viscous force of the damper
Q ˜ = k y ˜ + c y ˜ ˙ = ( k + j c ω ) Y ˜ e j ω t   ,  
or, using displacement Y ˜ , obtained by (16), and its derivate
Q ˜ = F 0 λ N E C 2 + D 2 [ ( D E C ) + j ( D E + C ) ]   ,
where terms E, C and D use the expressions from (10a) and (13a).
The amplitude | Q ˜ | = Q 0 of the dynamic force transmitted to the ground foundation is
| Q ˜ | = Q 0 = F 0 λ N E 1 + E 2 C 2 + D 2 = F 0 λ N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 C 2 + D 2   ,
or
Q 0 ( ζ ,   Ω ) = F 0 λ N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 [ Ω 2 + ( 2 ζ Ω ) 2 λ N Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2

2.3. Dynamic Transmissibility Coefficient

The dynamic coefficient of transmissibility of the force Q ˜ = Q 0 e j ( ω t + φ ) in relation to the harmonic perturbing force F ˜ = F 0 e j ω t is defined as the ratio between the two amplitudes of the forces, as [36,37]
T ( ζ ,   Ω ) = Q 0 ( ζ ,   Ω ) F 0 = λ N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 C 2 + D 2   ,
where :
{ C = Ω 2 ( 2 ζ Ω ) 2 λ N + Ω 2 ( 2 ζ Ω ) 2 λ D = ( 2 ζ Ω ) [ λ N Ω 2 ( 1 + λ + λ N ) ]

2.4. Deflections of the Viscous Dampers in Maxwell–Voigt–Kelvin Model

2.4.1. Deflection of Simple Viscous Damper in Voigt–Kelvin Model

The deflection of the simple viscous damper in the Voigt–Kelvin model (with damping coefficient c) is
y ˜ = Y ˜ e j ω t = Y 0 e j φ e j ω t = Y 0 e j ( ω t + φ )   ,
where:
φ is the phase between the displacements y ˜ and z ˜
Y 0 is the amplitude obtained by (18).

2.4.2. Deflection of Viscous Damper in Maxwell Model

The deflection of the viscous damper in the Maxwell model (with damping coefficient cS) is
u ˜ = x ˜ z ˜ = ( X ˜ Z ˜ ) e j ω t ,
where X ˜ and Z ˜ are obtained by (14) and (20), respectively.
The deflection of the viscous damper in the Maxwell model can be written as
u ˜ = U ˜ e j ω t   ,
where:
U ˜ = X ˜ Z ˜ = F 0 k [ A C + B D C 2 + D 2 + j B C A D C 2 + D 2 ] ( 2 ζ Ω ) λ F 0 k [ a C + b D C 2 + D 2 + j b C a D C 2 + D 2 ]   ,
The deflection U ˜ can be written as complex number
U ˜ = F 0 k 1 C 2 + D 2 [ G + j H ]   ,
where the expression of the amplitude | U ˜ | = U 0 is calculated as follows
U 0 2 = F 0 2 k 2 1 [ C 2 + D 2 ] 2 [ G 2 + H 2 ]   ,
or
U 0 = F 0 k G 2 + H 2 C 2 + D 2   ,
with the real constants:
{ G = ( A C + B D ) + ( 2 ζ Ω ) λ ( a C + b D ) H = ( B C A D ) + ( 2 ζ Ω ) λ ( b C a D )
By substituting the expressions for the constants G and H in Equation (34), the square of the amplitude U 0 of the deflection of the viscous damper in the Maxwell model is obtained as
U 0 2 = F 0 2 k 2 ( C 2 + D 2 ) [ ( A 2 + B 2 ) + ( 2 ζ Ω ) 2 λ 2 ( a 2 + b 2 ) + 2 ( 2 ζ Ω ) λ ( a A + b B ) ] ( C 2 + D 2 ) 2   ,
or
U 0 2 = F 0 2 k 2 L   ,
where
L = 1 C 2 + D 2 [ ( A 2 + B 2 ) + ( 2 ζ Ω ) 2 λ 2 ( a 2 + b 2 ) + 2 ( 2 ζ Ω ) λ ( a A + b B ) ]
The amplitude U 0 result is
U 0 = F 0 k L = F 0 k ( A 2 + B 2 ) + ( 2 ζ Ω ) 2 λ 2 ( a 2 + b 2 ) + 2 ( 2 ζ Ω ) λ ( a A + b B ) C 2 + D 2   ,
where a = 2 ζ Ω and b = 1 + N and the real constants A, B, C, D are obtained by (13a).

2.5. Dissipated Energy in the Viscous Dampers in Maxwell–Voigt–Kelvin Model

The dissipated energy per cycle in the viscous damper in the Voigt–Kelvin model (with damping coefficient c) is
W d V K = π c ω Y 0 2   ,
or, taking into consideration that c ω = ( 2 ζ Ω ) k and the expression of the amplitude Y 0 obtained by (18)
W d V K = π F 0 2 k λ 2 N 2 ( 2 ζ Ω ) 3 C 2 + D 2
The dissipated energy per cycle in the viscous damper in the Maxwell model (with damping coefficient cS = cλN) is
W d M = π ( c λ N ) ω U 0 2   ,
or, taking into consideration the square expression of the amplitude U 0 obtained by (18)
W d M = π F 0 2 k ( 2 ζ Ω ) λ N L
The total dissipated energy per cycle is the sum of dissipated energies in each type of viscous damper
W d = W d V K + W d M   ,
or
W d ( ζ , Ω ) = π F 0 2 k [ λ 2 N 2 ( 2 ζ Ω ) 3 C 2 + D 2 + λ N ( 2 ζ Ω ) L ]   ,
where the real constants C, D and L are obtained by (13a) and (38a).

3. Dynamic Parameters Analysis of a Vibratory Piling Machine: Case Study

The analysis of the dynamic parameters was performed for a calculus model of a vibratory piling machine with the following inertial and rheological characteristics [33]:
  • the mobile mass of the equipment (incl. pile mass) m = 10 4   kg
  • two cases of perturbing force:
    (1) harmonic force F ( t ) = F 0 sin ω t
    (2) harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ω t ;
  • the amplitude of harmonic perturbing force (piling force) F 0 = 1250   kN ;
  • the static moment of the dynamic unbalanced masses (for harmonic inertial force) m 0 r = 5   kg   m ;
  • discrete variable stiffness k = ( 1 ÷ 5 ) × 10 5   kN / m ;
  • discrete variable viscous damping c = ( 2 ÷ 15 ) × 10 2   kNs / m ;
  • discrete multiplier numbers N = 1 ÷ 10 and S = 1 ÷ 10 ;
  • perturbing force pulsation (steady-state regime) ω = ( 200 ÷ 500 )   rad / s ;
  • damping ratio of the ground foundation ζ = 0.15 ÷ 0.32 .

3.1. Amplitude of Steady-State Forced Vibration of Mobile Mass of the Equipment

Figure 2 shows the variation in steady-state vibration amplitudes of the mass m function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250   kN and multiplier numbers N = 10 , S = 10 for the Maxwell model. The curves of the amplitude X 0 ( ζ , Ω ) were plotted using Equation (15a).
Figure 3 shows the variation in steady-state vibration amplitudes of the mass m function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of inertial harmonic perturbing force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) and multiplier numbers N = 10 , S = 10 for the Maxwell model. The curves of the amplitude X 0 ( ζ , Ω ) were plotted using Equation (15a), with F 0 k = m 0 r m Ω 2 .

3.2. Amplitude of the Harmonic Deflection in Voigt–Kelvin Model

Figure 4 shows the variation in harmonic deflection in the Voigt–Kelvin model function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250   kN and multiplier numbers N = 10 , S = 10 for the Maxwell model. The curves of the deflection amplitude Y 0 ( ζ , Ω ) were plotted using Equation (18a).
Figure 5 shows the variation in harmonic deflection of the Voigt–Kelvin model function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of inertial harmonic perturbing force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) and multiplier numbers N = 10 , S = 10 for the Maxwell model. The curves of the deflection amplitude Y 0 ( ζ , Ω ) were plotted using Equation (18a), with F 0 k = m 0 r m Ω 2 .

3.3. Amplitude of the Harmonic Variation in the Displacement in Maxwell Model

Figure 6 shows the variation in harmonic displacement of the connection point CP in the Maxwell model, function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250 kN and multiplier numbers N = 10 , S = 10 . The curves of the displacement amplitude Z 0 ( ζ , Ω ) were plotted using Equation (22a).
Figure 7 shows the variation in harmonic displacement of the connection point CP in the Maxwell model, function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of inertial harmonic perturbing force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) and multiplier numbers N = 10 , S = 10 . The curves of the displacement amplitude Z 0 ( ζ , Ω ) were plotted using Equation (22a), with F 0 k = m 0 r m Ω 2 .

3.4. Dynamic Coefficient of Transmissibility

Figure 8 shows the variation in the coefficient of transmissibility T function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250   kN and multiplier numbers N = 10 , S = 10 for the Maxwell model. The curves of the coefficient of transmissibility T ( ζ , Ω ) were plotted using Equation (27).

3.5. Amplitude of the Harmonic Deflection in Maxwell Model

Figure 9a,b show the variation in harmonic deflection in the Maxwell model function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250   kN and multiplier numbers N = 10 , S = 10 . The curves of the deflection amplitude U 0 ( ζ , Ω ) were plotted using Equation (39).
Figure 10 shows the variation in harmonic deflection in the Maxwell model function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) and multiplier numbers N = 10 , S = 10 . The curves of the deflection amplitude U 0 ( ζ , Ω ) were plotted using Equation (39), with F 0 k = m 0 r m Ω 2 .

3.6. Dissipated Energy in the Viscous Dampers of Maxwell–Voigt–Kelvin Model

Figure 11 shows the variation in the total dissipated energy per cycle in the dampers of the Maxwell–Voigt–Kelvin model coefficient function of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of harmonic perturbing force F ( t ) = F 0 sin ( ω t ) with constant amplitude F 0 = 1250 kN and multiplier numbers N = 10 , S = 10 . The curves of the dissipated energy per cycle W d ( ζ , Ω ) were plotted using Equation (40).
Figure 12 shows the variation in the total dissipated energy per cycle in the dampers of the Maxwell–Voigt–Kelvin model coefficient of relative pulsation Ω , at four discrete/different values of damping ratio ζ, in the case of inertial harmonic perturbing force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) and multiplier numbers N = 10 , S = 10 . The curves of the dissipated energy per cycle W d ( ζ , Ω ) were plotted using Equation (40), with F 0 k = m 0 r m Ω 2 .

4. Discussion

(1) Hypothetically, in the case of a very high damping coefficient in the Maxwell rheological model, i.e., c S or λ , the dynamic rigidity of the damper increases (theoretical c ω S ) and the Maxwell model becomes a Hooke rheological model [32,38]. In this case, the Maxwell–Voigt–Kelvin rheological model (see Figure 1) becomes a Hooke–Voigt–Kelvin rheological model as seen in Figure 13. Thus, the amplitude of the dynamic harmonic force transmitted to the ground foundation is
Q 0 H V K = lim λ Q 0 M V K = lim λ F 0 λ N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 [ Ω 2 ( 2 ζ Ω ) 2 λ N + Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2   ,
or
Q 0 H V K = lim λ F 0 N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 [ Ω 2 λ ( 2 ζ Ω ) 2 N + Ω 2 ( 2 ζ Ω ) 2 ] 2 + ( 2 ζ Ω ) 2 [ N Ω 2 ( 1 λ + 1 + N ) ] 2
After the limit calculation, the amplitude of dynamic harmonic transmitted force becomes
Q 0 H V K = F 0 N ( 2 ζ Ω ) 1 + ( 2 ζ Ω ) 2 ( 2 ζ Ω ) 4 ( Ω 2 N ) 2 + ( 2 ζ Ω ) 2 [ N Ω 2 ( 1 + N ) ] 2   ,
or
Q 0 H V K = F 0 N 1 + ( 2 ζ Ω ) 2 ( 2 ζ Ω ) 2 ( Ω 2 N ) 2 + [ N Ω 2 ( 1 + N ) ] 2  
(2) Under the same hypothesis as in (1), the dynamic coefficient of the transmissibility of Hooke–Voigt–Kelvin rheological model is:
T H V K ( ζ , Ω ) = Q 0 H V K F 0 = N 1 + ( 2 ζ Ω ) 2 ( 2 ζ Ω ) 2 ( Ω 2 N ) 2 + [ N Ω 2 ( 1 + N ) ] 2   ,
(3) Considering the Hooke–Voigt–Kelvin rheological model according to (1), the amplitude of steady-state forced vibration of mobile mass m of the piling machine can be written
X 0 H V K = lim λ X 0 M V K = F 0 k lim λ [ 1 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 ( 1 + λ + λ N ) 2 [ Ω 2 + ( 2 ζ Ω ) 2 λ N Ω 2 ( 2 ζ Ω ) 2 λ ] 2 + ( 2 ζ Ω ) 2 [ λ N Ω 2 ( 1 + λ + λ N ) ] 2   ,
or
X 0 H V K = F 0 k lim λ [ 1 λ ( 2 ζ Ω ) 2 ] 2 + ( 2 ζ Ω ) 2 [ 1 λ + 1 + N ] 2 [ Ω 2 λ ( 2 ζ Ω ) 2 N + Ω 2 ( 2 ζ Ω ) 2 ] 2 + ( 2 ζ Ω ) 2 [ N Ω 2 ( 1 λ + 1 + N ) ] 2
X 0 H V K = F 0 k ( 2 ζ Ω ) 2 + ( 1 + N ) 2 ( 2 ζ Ω ) 2 ( Ω 2 N ) 2 + [ N Ω 2 ( 1 + N ) ] 2

5. Conclusions

Using a dynamic model and the hypothesis of linear behavior in stabilized dynamic regime (steady-state forced vibration), and based on the calculation relations, the parametric quantities defining the operation of the vibratory piling machines (equipment for the introduction of piles/columns in the ground foundation) were evaluated. Thus, the conclusions in the case of perturbing harmonic forces F ( t ) = F 0 sin ( ω t ) and F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) can be summarized as follows:
(a) The amplitude of the technological vibrations in a post-resonance regime remains relatively constant at the variation of the relative pulsation Ω 2 ;
(b) The variation in the maximum displacements and deflections, as internal parameters of the system, expressed by X 0 ( ζ , Ω ) , Y 0 ( ζ , Ω ) , Z 0 ( ζ , Ω ) and U 0 ( ζ , Ω ) in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 and Figure 9 and Figure 10, highlights the internal behavior of the ground foundation in the regime of steady-state harmonic forced vibrations with the maximum force P = F 0 = 2000   kN and F 0 ( ω ) = ( m 0 r ω 2 ) ,   respectively ;
(c) The amplitude of the dynamic force Q 0 ( ζ ,   Ω ) transmitted to the ground foundation, function of the parameters of the piling machine and the rheological parameters of the soil, can be calculated with Equation (26);
(d) The ability to transmit the technological vibrations at the top of the pilot/column to the ground foundation is given by Equation (27), with the parametric variation in the coefficient of transmissibility T from Figure 8;
(e) The dissipated energy during the process of vibration introduction of the pilot/column into the soil can be calculated with Equation (40), with the parametric variations from Figure 11 and Figure 12;
(f) From the curves of the amplitude of steady-state forced vibration of mobile mass m (Figure 2 and Figure 3), it can be concluded that the stabilized technological regime for the vibratory piling machine is in post-resonance, for Ω 2 (where the amplitude of forced vibration is the same, regardless of the value of damping ratio ζ ). The post-resonance operating condition Ω 2 can be ensured by increasing the total moving mass m (depending on the mass of the pilot/column, as well as the elasticity of the foundation soil) or by increasing the pulsation ω of the harmonic force, or both;
(g) The use of the vibratory piling machine in the resonance stabilized regime ( Ω 1 regardless of the size of the structural damping ratio ζ of the ground) could lead to larger amplitudes of the steady-state forced vibration: 40 70   mm in case of harmonic force F 0 sin ( ω t ) see Figure 2, 9 16   mm in case of inertial harmonic force ( m 0 r ω 2 ) sin ( ω t ) see Figure 3. However, although the resonant operating mode could lead to higher forward speeds of the pilot/column into the foundation ground (hence higher productivity of the technological process), with the change of the damping ratio ζ of the foundation soil, the amplitude of the forced vibrations varies significantly, the process of piling through vibrations becoming technologically unstable and, therefore, uncontrollable. According to the case studies presented in Figure 2 and Figure 3, the amplitude of the forced vibrations in the resonance mode varies significantly depending on the damping ratio ζ , with the linear functions A res ( ζ ) = 78.9 121.3 ζ   [ mm ] in the case of harmonic force F 0 sin ( ω t ) and A ω res ( ζ ) = 21.8 39.4 ζ   [ mm ] in the case of inertial harmonic force ( m 0 r ω 2 ) sin ( ω t ) ;
(h) Considering a technological regime Ω 2 for the vibratory piling equipment, the transmissibility coefficient of the force is T 0.5 regardless of the value of damping ratio ζ (Figure 8), meaning that the transmitted force to the pile/column is, at maximum, half of the dynamic force of the vibrating motor. In the case of a post-resonance regime with Ω 2 , the higher the damping ratio ζ , the higher the transmissibility coefficient T, due to the forces transmitted by the dampers of the Maxwell and Voigt–Kelvin rheological elements, especially for high deformation speeds. Despite the fact that the operation of the technological equipment is technologically unstable in the resonance work mode, high transmissibility coefficients can be obtained, all the more so as the damping coefficient of the foundation ground decreases. As Figure 8 demonstrates, the transmissibility coefficient at the resonance varies inversely with the damping ratio ζ of the foundation ground according to the law T res ( ζ ) = 4 6.07 ζ for 0.15 ζ 0.32 ;
(i) The curves in Figure 11 and Figure 12 show that the dissipated energies per cycle are maximum at resonance and the higher the value, the lower the damping ratio ζ (due to the higher values of the damping deflections at the lower values of the damping coefficients). For the technological regime with Ω 2 , the values of dissipated energy are 50 100   kJ / cycle for harmonic force F ( t ) = F 0 sin ( ω t ) and 10 20 kJ / cycle for inertial harmonic force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) .
(j) According to Figure 11, for a hypothetical stabilized resonance work regime, the theoretical value of dissipated energy in the case of harmonic perturbing force is 1300 2600   kJ / cycle , meaning a dissipated power P 40 , 000 80 , 000   kW for a nominal pulsation of ω = 200   rad / s ( f = 31.8   Hz or n = 1909   min 1 ). According to Figure 12, in the same hypotheses as above, the value of dissipated energy in the case of inertial harmonic perturbing force is 77 116   kJ / cycle , meaning a dissipated power P 2500 3700   kW . These huge values of dissipated powers lead to the conclusion that resonant operating modes, even when stable, require very high installed powers of the drive vibro-motors;
(k) For a nominal pulsation of ω = 200   rad / s , a vibromotor with rotating unbalanced masses and a technological regime with Ω 5 , the dissipated power in damping elements is about P 60 150   kW or P 80 200   HP . Taking into consideration the mechanical efficiency of the vibratory piling machines used and the in situ experimental test values (for the engines nominal power of P N 200 400   HP ), the complex linear Hooke–Voigt–Kelvin rheological model with variable damping coefficient can be used for the dynamic analysis of the vibrogenerator–pile–foundation soil system;
(l) Regardless of the type of disturbing force generated by the vibratory piling machine, the harmonic force F 0 sin ( ω t ) and inertial harmonic force ( m 0 r ω 2 ) sin ( ω t ) , engineering practice on site has shown that the pre-resonant and resonant operating modes are unstable and uncontrollable, therefore, not recommended from a technological perspective.
The dynamic calculation model presented in the paper can be used to evaluate the technological working capacity of a vibratory piling machine or equipment used for inserting piles/columns in various types of soils for the construction of foundations.
The limitations of the presented models are as follows:
  • the rheological model has linear elements (elasticity, damping) and cannot describe the possible nonlinear behavior of the foundation soil;
  • the dynamic model of the system considers the pilot/column as a rigid solid body with infinite stiffness;
  • the pilot’s advancement in the foundation ground does not take into account any possible obstacles in the foundation ground (e.g., very hard rocks) or irregularities (e.g., cavities, liquefied soil, groundwater).

Author Contributions

Conceptualization, P.B., N.D.; methodology, A.M.G. and P.B.; software, A.M.G. and N.D.; validation, A.M.G., P.B. and N.D.; formal analysis, N.D. and A.M.G.; investigation, A.M.G. and N.D.; data curation, N.D.; writing—original draft preparation, P.B. and N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tonciu, O. Parametric analysis of vibrating equipment for the insertion of piles into the ground. ACTA Tech. Napoc. Appl. Math. Mech. Eng. 2020, 63, 95–100. [Google Scholar]
  2. Forssblad, L. Vibratory Soil and Rock Fill Compaction; Robert Olsson Tryckeri AB: Stockholm, Sweden, 1981. [Google Scholar]
  3. Capatana, G.F. Dynamic Behaviour of Complex Interaction Between Vibratory Drum Equipment and Natural Terrain Based on Rheological Evaluations. In Proceedings of the 10th HSTAM International Congress on Mechanics Chania, Crete, Greece, 25–27 May 2013. [Google Scholar]
  4. Capatana, G.F. Dynamic Simulation of the Vibratory Roller-Terrain Interaction Using an Elastoplastic Approach. In Annals of Dunărea de Jos University of Galaţi, XIV, Mechanical Engineering; Galati University Press: Galati, Romania, 2013. [Google Scholar]
  5. Lambe, T.W.; Whitman, R.V. Soil Mechanics; John Wiley&Sons: New York, NY, USA, 1969. [Google Scholar]
  6. Bejan, S. Dynamic Response Analysis of the Road System Compaction According to the Forced Vibration Mode. Rom. J. Acoust. Vib. 2014, XI, 164–166. [Google Scholar]
  7. Ishibashi, I.; Zhang, X. Unified dynamic shear moduli and damping ratios of sand and clay. Soils Found. 1993, 33, 182–191. [Google Scholar] [CrossRef] [Green Version]
  8. Colaço, A.; Ferreira, M.A.; Costa, P.A. Empirical, Experimental and Numerical Prediction of Ground-Borne Vibrations Induced by Impact Pile Driving. Vibration 2022, 5, 80–95. [Google Scholar] [CrossRef]
  9. Bay, J.A. A Summary of the Research on Pile Driving Vibrations. In Proceedings of the Pile Driving Contractor’s Association, 7th Annual Winter Round, Atlanta, GA, USA, 21–22 February 2003. [Google Scholar]
  10. Fernandes, M.M. Analysis and Design of Geotechnical Structures, 1st ed.; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  11. Hunt, C.E.; Pestana, J.M.; Bray, J.D.; Riemer, M. Effect of pile driving on static and dynamic properties of soft clay. J. Geotech. Geoenviron. Eng. 2002, 128, 13–24. [Google Scholar] [CrossRef]
  12. Wang, S.; Zhu, S. Global Vibration Intensity Assessment Based on Vibration Source Localization on Construction Sites: Application to Vibratory Sheet Piling. Appl. Sci. 2022, 12, 1946. [Google Scholar] [CrossRef]
  13. Deckner, F.; Viking, K.; Hintze, S. Ground vibrations due to pile and sheet pile driving: Prediction models of today. In Proceedings of the 22nd European Young Geotechnical Engineers Conference, Gothenburg, Sweden, 26–29 August 2012; pp. 107–112. [Google Scholar]
  14. Kim, D.S.; Lee, J.S. Propagation and attenuation characteristics of various ground vibrations. Soil Dyn. Earthq. Eng. 2000, 19, 115–126. [Google Scholar] [CrossRef]
  15. Dai, D.; Zhang, Y.; Zhang, Y.; Wang, Z.; Li, Z. Kinematic Response of End-Bearing Piles under the Excitation of Vertical P-Waves Considering the Construction Effect. Appl. Sci. 2022, 12, 3468. [Google Scholar] [CrossRef]
  16. Li, Q.; Li, X.; Wen, M.; Hu, L.; Duan, W.; Li, J. Dynamic Responses of a Pile with a Cap under the Freezing and Thawing Processes of a Saturated Porous Media Considering Slippage between Pile and Soil. Appl. Sci. 2022, 12, 4214. [Google Scholar] [CrossRef]
  17. Varghese, R.; Boominathan, A.; Banerjee, S. Stiffness and load sharing characteristics of piled raft foundations subjected to dynamic loads. Soil Dyn. Earthq. Eng. 2020, 133, 106117. [Google Scholar] [CrossRef]
  18. Liming, Q.; Xuanming, D.; Changjie, Z.; Chongrong, W.; Guangwei, C. Numerical and test study on vertical vibration characteristics of pile group in slope soil topography. Earthq. Eng. Eng. Vib. 2021, 20, 377–390. [Google Scholar] [CrossRef]
  19. Yang, X.; Wang, L.; Wu, W.; Liu, H.; Jiang, G.; Wang, K.; Mei, G. Vertical Dynamic Impedance of a Viscoelastic Pile in Arbitrarily Layered Soil Based on the Fictitious Soil Pile Model. Energies 2022, 15, 2087. [Google Scholar] [CrossRef]
  20. Kwon, S.Y.; Yoo, M. Evaluation of Dynamic Soil-Pile-Structure Interactive Behavior in Dry Sand by 3D Numerical Simulation. Appl. Sci. 2019, 9, 2612. [Google Scholar] [CrossRef] [Green Version]
  21. Halder, P.; Manna, B. Large scale model testing to investigate the influence of granular cushion layer on the performance of disconnected piled raft system. Acta Geotech. 2021, 16, 1597–1614. [Google Scholar] [CrossRef]
  22. Rashidifar, M.A.; Rashidifar, A.A.; Abertavi, A. Nonlinear characteristics of the pile soil system under vertical vibration. Univers. J. Eng. Sci. 2016, 4, 59–65. [Google Scholar] [CrossRef]
  23. Fattah, M.Y.; Zabar, B.S.; Mustafa, F.S. Effect of saturation on response of a single pile embedded in saturated sandy soil to vertical vibration. Eur. J. Environ. Civ. Eng. 2020, 24, 381–400. [Google Scholar] [CrossRef]
  24. Machacek, J.; Staubach, P.; Tafili, M.; Zachert, H.; Wichtmann, T. Investigation of three sophisticated constitutive soil models: From numerical formulations to element tests and the analysis of vibratory pile driving tests. Comput. Geotech. 2021, 138, 104276. [Google Scholar] [CrossRef]
  25. Wei, J.; Wang, W.; Wu, J. Hydro-Mechanically Coupled Numerical Modelling on Vibratory Open-Ended Pile Driving in Saturated Sand. Appl. Sci. 2022, 12, 4527. [Google Scholar] [CrossRef]
  26. Vogelsang, J.; Huber, G.; Triantafyllidis, T. Experimental Investigation of Vibratory Pile Driving in Saturated Sand; Springer International Publishing: Cham, Switzerland, 2017; pp. 101–123. [Google Scholar]
  27. O’Neill, M.W.; Vipulanandan, C.; Wong, D. Laboratory modeling of vibro-driven piles. J. Geotech. Eng. 1990, 116, 1190–1209. [Google Scholar] [CrossRef]
  28. Spanu (Stefan), G.C.; Capatana, G.F.; Potirniche, A. Analysis of the Transmissibility Ratio and the Isolation Degree of the Vibration for 1DOF Mechanical Systems with Zener Viscous Damping Model. Synth. Theor. Appl. Mech. 2019, 10, 43–56. [Google Scholar]
  29. Capatana, G.F.; Potirniche, A.; Spanu (Stefan), G.C. The Dynamics of the 1 DOF Mechanical Systems with Viscous Damping Zener Model. Amplitude Factor Analysis. Synth. Theor. Appl. Mech. 2019, 10, 25–34. [Google Scholar]
  30. Spanu (Stefan), G.C.; Capatana, G.F.; Potirniche, A. Comparative Analysis of the Dynamic Parameters of Mechanical Systems 1DOF with Complex Reological Behavior Elastomeric Bearings. Synth. Theor. Appl. Mech. 2019, 10, 225–230. [Google Scholar]
  31. Bratu, P.; Dobrescu, C. Evaluation of the Dissipated Energy in Vicinity of the Resonance, depending on the Nature of Dynamic Excitation. Rom. J. Acoust. Vib. 2019, 16, 66–71. [Google Scholar]
  32. Dobrescu, C.F. Dynamic Response of the Newton Voigt–Kelvin Modelled Linear Viscoelastic Systems at Harmonic Actions. Symmetry 2020, 12, 1571. [Google Scholar] [CrossRef]
  33. Bratu, P. Hysteretic Loops in Correlation with the Maximum Dissipated Energy, for Linear Dynamic Systems. Symmetry 2019, 11, 315. [Google Scholar] [CrossRef] [Green Version]
  34. Rao, M. Mechanical Vibrations; Addison-Wesley Pub. Co.: Boston, MA, USA, 1986. [Google Scholar]
  35. Radeș, M. Mechanical Vibrations; Editura Printech: București, Romania, 2006. [Google Scholar]
  36. Vasile, O. Active Vibration Control for Viscoelastic Amortization Systems under the Action of Inertial Forces. Rom. J. Acoust. Vib. 2017, 14, 54–58. [Google Scholar]
  37. Sireteanu, T.; Giuclea, M.; Mitu, A.M. An analitical approach for approximation of experimental hysteretic by Bouc-Wen model. Procc. Rom. Acad. Ser. A 2009, 10, 43–54. [Google Scholar]
  38. Muscă (Anghelache), D.G.; Potîrniche, A.; Căpățână, G.F. Analyse of dynamic operation mode for drilling machine equipped with drill for pits digging. In Proceedings of the ModTech International Conference 2020, Modern Technologies in Industrial Engineering Iasi, Mamaia, Romania, 23–27 June 2020. [Google Scholar]
Figure 1. Scheme of the linear dynamic system Maxwell–Voigt–Kelvin (E-V)–(E|V) perturbed by a harmonic force F ( t ) = F 0 sin ( ω t ) .
Figure 1. Scheme of the linear dynamic system Maxwell–Voigt–Kelvin (E-V)–(E|V) perturbed by a harmonic force F ( t ) = F 0 sin ( ω t ) .
Symmetry 14 01238 g001
Figure 2. The amplitude of steady-state vibration of mass m in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Figure 2. The amplitude of steady-state vibration of mass m in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Symmetry 14 01238 g002
Figure 3. The amplitude of steady-state vibration of mass m in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Figure 3. The amplitude of steady-state vibration of mass m in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Symmetry 14 01238 g003
Figure 4. The amplitude of the harmonic variation in the deflection of Voigt–Kelvin model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Figure 4. The amplitude of the harmonic variation in the deflection of Voigt–Kelvin model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Symmetry 14 01238 g004
Figure 5. The amplitude of the harmonic variation in the deflection of Voigt–Kelvin model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Figure 5. The amplitude of the harmonic variation in the deflection of Voigt–Kelvin model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Symmetry 14 01238 g005
Figure 6. The amplitude of the harmonic variation in the displacement of the point of connection in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Figure 6. The amplitude of the harmonic variation in the displacement of the point of connection in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Symmetry 14 01238 g006
Figure 7. The amplitude of the harmonic variation in the displacement of the point of connection in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Figure 7. The amplitude of the harmonic variation in the displacement of the point of connection in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Symmetry 14 01238 g007
Figure 8. Dynamic coefficient of transmissibility T function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 2000   kN .
Figure 8. Dynamic coefficient of transmissibility T function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 2000   kN .
Symmetry 14 01238 g008
Figure 9. The amplitude of the harmonic variation in the deflection of the damper in Maxwell model in dynamic regime—variation in amplitude function of ζ and   Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN ; (a) Ω [ 0 ,   10 ] ; (b) detail for Ω [ 0 ,   0.1 ] .
Figure 9. The amplitude of the harmonic variation in the deflection of the damper in Maxwell model in dynamic regime—variation in amplitude function of ζ and   Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN ; (a) Ω [ 0 ,   10 ] ; (b) detail for Ω [ 0 ,   0.1 ] .
Symmetry 14 01238 g009
Figure 10. The amplitude of the harmonic variation in the deflection of the damper in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Figure 10. The amplitude of the harmonic variation in the deflection of the damper in Maxwell model in dynamic regime—variation in amplitude function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Symmetry 14 01238 g010
Figure 11. The dissipated energy per cycle W d   [ J ] function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Figure 11. The dissipated energy per cycle W d   [ J ] function of ζ and Ω ; harmonic force F ( t ) = F 0 sin ( ω t ) , F 0 = 1250   kN .
Symmetry 14 01238 g011
Figure 12. The dissipated energy per cycle W d   [ J ] function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Figure 12. The dissipated energy per cycle W d   [ J ] function of ζ and Ω ; harmonic inertial force F ( t ) = ( m 0 r ω 2 ) sin ( ω t ) = m 0 r m k Ω 2 sin ( ω t ) .
Symmetry 14 01238 g012
Figure 13. Scheme of the Hooke–Voigt–Kelvin rheological model (E)–(E|V) perturbed by a harmonic force F ( t ) = F 0 sin ( ω t ) .
Figure 13. Scheme of the Hooke–Voigt–Kelvin rheological model (E)–(E|V) perturbed by a harmonic force F ( t ) = F 0 sin ( ω t ) .
Symmetry 14 01238 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Goanță, A.M.; Bratu, P.; Drăgan, N. Dynamic Response of Vibratory Piling Machines for Ground Foundations. Symmetry 2022, 14, 1238. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14061238

AMA Style

Goanță AM, Bratu P, Drăgan N. Dynamic Response of Vibratory Piling Machines for Ground Foundations. Symmetry. 2022; 14(6):1238. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14061238

Chicago/Turabian Style

Goanță, Adrian Mihai, Polidor Bratu, and Nicușor Drăgan. 2022. "Dynamic Response of Vibratory Piling Machines for Ground Foundations" Symmetry 14, no. 6: 1238. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14061238

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