95
Original Article
Normal Impact of A Rigid Cone-shaped against A Viscoelastic Plate on Viscoelastic Foundation
Duong Tuan Manh
*VNU University of Engineering and Technology, 144 Xuan Thuy, Cau Giay, Hanoi, Vietnam Received 10 July 2021
Revised 18 July 2021; Accepted 18 July 2021
Abstract: This paper considers the problem of a low-velocity normal impact of a rigid cone-shaped upon a viscoelastic plate. The contact force is defined by the modified Hertz's contact law.
Approximate solutions of the system of nonlinear integro-differential equations for the contact force and local indentation have been obtained.
Keywords: Normal impact, viscoelastic plate, Kelvin-Voigt model.
1. Introduction *
The problems with the analysis of the impact of two bodies have a widespread application in various fields of science and technology. Because these problems belong to dynamic contact interaction, their solutions are connected with severe mathematical and calculation difficulties. Some approaches and methods have been suggested in articles by Abrate [1], Rossikhin and Shitikova [2].
Rossikhin et al., [3] investigated the collision of two viscoelastic shells, viscoelastic features of which are described by the standard linear solid model. During the impact process, the local bearing of the materials of the colliding viscoelastic shells is taken into account, the solution in the contact domain is found via the modified Hertz contact theory involving the operator representation of viscoelastic analogs of Young’s modulus and Poisson’s ratio.
In the paper, we will consider a special but very important for engineering practice case when a rigid cone-shaped impacts a viscoelastic plate.
________
* Corresponding author.
E-mail address: iam.mr.manh@gmail.com https//doi.org/10.25073/2588-1124/vnumap.4661
2. Problem Formulation
Let us consider the problem on a transverse impact of a rigid cone-shaped upon a viscoelastic Kirchhoff-Love plate when the viscoelastic features of the surrounding medium are described by a Kelvin-Voigt model. In this case, the equations of motion of a rigid cone-shaped of mass ms and the viscoelastic rectangular plate with the dimensions a and b and of thickness h have, respectively, the form
Figure 1. Scheme of the normal impact of a rigid cone-shaped against a plate.
,
( ) ( )
s
t
m z t = − P
(1)2
( , , ) ( , , ) ( , ) ( )
2 2
( , , )
( , , ) ( , , )
w T G
a b
D w x y t hw x y t L w F P t x y w x y t
K w x y t C K w x y t
t
+ − = − −
− − +
(2)
1
21
( , ) F 2 L w w
E = −
, (3)where ρ is the density of plate, w(x,y,t) is the displacements of an arbitrary point in the plate in the z- direction, P(t) is the contact force, F is Airy function stress,
is the Laplace operator,2 x a
−
,2 x b
−
is the Dirac delta function, d and z(t) are the penetration depths, θ is the angle between the?
h d
? z
Cone-shaped
Viscoelastic plate
Foundation
Damping coefficient Spring foundation Shear foundation
plane and cone,
K
w,K
G,C
T are the Winkler, Pasternak and Damper modulus parameters, respectively, and1
tE E
E
= +
( )
3
12 1
2D Eh
=
−
Utilizing the KelvinVoigt model with the bulk modulus K remaining constant during the process of mechanical loading of this material [4], it is found [3] that
( ) ( )
2
1
1 E E
tL
E L
= + −
−
,where E, , are the relaxed Young’s modulus, Poisson’s ratio and viscosity, respectively,
1
E
=
−=
is the creep time, and( )
( ) ( )
( ) ( ) ( ) ( ) ( )
0
1 ;
2 1 2 1
; ; 2 exp ex '
1 2 1 2 p '
4
v v
tL t t dt
−
− −
− + − −
= = = = −
Equations (1), (2) and (3) are subjected to the following initial conditions:
(0)
0, (0) 0, (0) 0
z = V z = P =
(4)( , ,0) 0, ( , ,0) 0
w x y = w x y =
(5)(0, , ) 0, ( , , ) ( ,0, ) 0, ( , , ) 0
w y t = w a y t = w x t = w x b t =
(6)(0, , ) ( , , ) (0, , ) ( , , ) 0
x x y y
M y t = M a y t = M y t = M b y t =
(7)1; 2;
( , , ) ( , , ) ( , , ) 0 : 0,
x y xy
0,
x a
N x y t p N x y t p N x y t at
y b
= = = − =
=
, (8)where all four edges of a rectangular plate are simply supported and immovable, V0 is the initial velocity of the collision.
Integrating twice Eq. (1) yields
0 0
( ) 1 ( ')( ') '
t
s
z t P t t t dt V t
= − m − +
(9)In the case of the collision of a rigid cone-shaped against a plate (Figure 1), the solution in the contact domain could be found using Hertz’ theory. Thus, the contact force is defined as [5]
( )
2 2
2 ( )
( ) ,
1 tan
E z t
P t = −
(10)Replacing Eq. (1) into Eq. (10) yields
( )
2 2
2 ( )
( ) 0
1 tan
s
E z t m z t
+ =
−
(11)The approximate solution of the equation (3) satisfying the boundary conditions (6), (7) and (8) can be written as
w( , , ) sin m x sin n y x y t f
a b
=
(12)2 2 2 2
2 1 2
2 2
2 2
cos cos
32 E 1
ta m x b n y p y 2 p x 2
F f
E b a a b
= + + +
+
, (13)where f is function of time t, m, n are odd numbers.
Replacing Eq. (13) into Eq. (2) yields
( )
2 2 4 4 4 4 4
2 2 4 4
2 2
2 2 2 2
2 1 2 2
2 2 2 2
3 2
2
sin sin 2 sin sin
sin sin 4 cos2
32
sin si 4 2
12
c 1
n o
2 1
1 s
3
t
t
E m x n y m n m n m x n y
f f
a b a b a b a b
m m x n y E n y
f f n a p
a a b E b
n m x n y E m x
f f m b
b a b E a
h h
−
−
+ +
− −
+
+ +
−
−
−
−
+ +
2 2
2 2
2
2
2 sin sin sin sin
( ) 2
( sin sin sin sin )
w T
G
p m x n y
a b
P t x y K C
m n
K n
m x n y
f f
a b a b
m x y m x n y
a b a
f a b b
+
=
+
− − − −
−
(14)
With due account for orthogonality of sines on the segments 0 ≤ x ≤ a; 0 ≤ y≤ b, we obtain:
( )
2 2 4 4 4 4 4
2 2 4 4
2 2
2 2 2 2
2 1
2 2
2 2 2 2
2 2 3
2 sin sin 2
sin sin
sin sin 4 cos2
32
sin sin 4 co 2
1
1
s 2
32 1
1
t
t
E m x n y m n m n
f a b a b a b
m x n y
f a b
m m
E h
x n y E n y
f f n a p
a a b
m h
E b
n m x n y m x
f f b p
b a b E a
−
−
+ +
+
−
− +
+
+
−
−
0 0
2 2
2 2
2
sin sin
s
0
in sin sin sin
sin n
( ) 2 2
( sin sin si )
a b
w T
G
m x n y
a b
m x n y m x n y
f f
a b a b
m x n y m x n y
a b
dydx
a b
P t x y K C
m n
K f a b a b
=
− −
+ +
− +
+
The result of which could be written in the form of
( ) ( ) ( ) ( )
3 21 2 3 4 5 6
( )
P t = A L − L f + A L − L f + A f + A f + A f + A f f
, (15)where
2 2 4
1 2 2
2 2 4
2 2 2
2 2 2 2 2 2
3 2 1 2 2 2 2
2 2 4 2 2 4
5 3 3
2 2 4 2
2
2 6
2 3
3
4
4
3 3
7 ;
48
48 4
4 4 4
4
64 64
64 tan
64 2
T
w G
m n ab
A E
a b
m n ab ab
A a b
m n ab ab ab m n
A p p
a b a b
ab
m n b m n a
A E E
a b
m n b m n h
h C
K K
A h
A a
a b
A
= +
= + +
= + + + +
+
= +
=
=
=
(16) The system of nonlinear integro-differential equations describing the dynamic behavior of two bodies during their collision takes the form of
( )
( ) ( )
( ) ( )
7
2
3 2
1 2 3 5 6
2
4
( )
( ) ( ) 0
1
( ) 0
( ) s
s
s
P
f m z t
m z t E z t
m z t f f
t A
A L L f A L L f A f A f A A
+
+ =
−
=
+ − + + + =
−
− +
(17)
3. Approximate Calculations
First we substitute the approximate solution of the equation (11) to form
( ) ( )
2( ) ( )
Hz t = N z t + M z t z t
, (18)where
( ) ( ) ( ) ( )
7
7 s
s
N A E L L
m
M A L L
m
= − −
= − −
0
( )
k( )
k
z t a t
=
= (19)
Replacing Eq. (19) into Eqs. (18) yields
' '' ''
'' ''' '' '''
0
0 0 0 0 0 0 0
' '' ''
'' ''' '' '''
0 8
9
0 0
0 0 0
( ) ( ( ''') ( ''')) ''' '' '
( ( ''') ( ''')) ''' '' '
t t t t
t t t t
k k k
k k k
t t t t
t t t t
k k
k k
a t V t e e A t e e A t dt dt dt
e e B t e e B t dt
A
A dt dt
− −
= = =
− −
= =
= −
−
−
−
, (20)
where
7 8
7 9
s
s
A A E m A A
m
=
=
Note that
( ) ( ) ( )
1 1
0 0
1 1
! !
!
i k
t k
t x k k i t
i k
i
e e x dx k t k e
k i
− − −
+ +
=
− −
= −
−
(21)The first two of them have the form
4 3
8 8 9
2 2
8 3 3 9 2 2 8 4 4
2
9 3 3 8 5 5 9
0
2 0
1 0
2 2 2
0 0
2 2 2
0 0 0
(
1 1 1 1 1 1
12
1 )
( )
1 1
6
1 6
3 3
1 1
4 1 4
1
o
t t
a t V t
V A
a t V
V t V t
e
A A
t t
A
e e
V t V
A A
A A A V
− − −
− − + − − −
− − + −
=
=
+ −
− − +
− −
4 48 5 5 9
2 2
0 0 4 4
1 1 1 1
6 4
t t
A V A
e
V
−
−
− − + −
(22)
We obtain the solution
4 3
8 8 9
2 2
8 3 3 9 2 2 8
2 0 2
0
2 2 2
0
4 4
9 8 5 9
0 0
2 2 2
3 0 5 0
3 4
1 1 1 1 1 1
( ) 12 3 3
1 1 1 1 1 1
6
1 1
4 6 4
t o
t t
A A A
z
t
t V t t
V t V
V t V
e e e
V
A A A
V
A A e
V t A
− − − −
− − + − − −
− − + −
+ −
− − + − − −
4
8 5 5 9
2 2
0 4 4
0
1 1 1 1
6 4
t
AV AV
− − + −
(23)
In the experiments of the collision problems, the collision time is about 10-3s, so the series (19) converges very quickly. Here, we only need the first two terms to approximate the value of the solution.
Another method such as balancing the coefficients of the series can be found in [3].
The obtained equation by replacing Eq. (23) into Eq. (15) can be used for nonlinear dynamical analysis of the collision.
4. Conclusion
In the present paper, the problem on the normal impact of a rigid cone-shaped upon a viscoelastic plate has been studied using the damping features of the impactor modeled by the linear Kelvin-Voigt model. An approximate analytical solution has also been found.
Acknowledgments
This work was supported by VNU University of Engineering and Technology under Project CN.19.05.
References
[1] S. Abrate, Impact on Composite Structures, Cambridge University Press, Cambridge, 1998.
[2] Y. A. Rossikhin, M. V. Shitikova, Analysis of Two Colliding Fractionally Damped Spherical Shells in Modelling Blunt Human Head Impacts, Central European Journal of Physics, Vol. 11, 2013, pp. 760-778, https://doi.org/10.2478/s11534-013-0194-4.
[3] Y. A. Rossikhin, M. V. Shitikova, D. T. Manh, Modelling of the Collision of Two Viscoelastic Spherical Shells, Mechanics of Time-Dependent Materials, Vol. 20, 2016, pp. 481-509, https://doi.org/10.1007/s11043-016-9308-x.
[4] Y. N. Rabotnov, Creep of Structure Elements, Naka, Moscow, 1966.
[5] I. N. Sneddon, The Relation Between Load and Penetration in The Axisymmetric Boussinesq Problem for A Punch of Arbitrary Profile, International Journal of Engineering Science, Vol. 3, 1965, pp. 47-57, https://doi.org/10.1016/0020-7225(65)90019-4.