# The Exact Solutions for a Point Mass Moving along a Stretched String on a Winkler Foundation

Authors: Q. Gao; J. Zhang; H. W. Zhang; W. X. Zhong
Journal: Shock and Vibration
 (2014)

---

## Abstract

This paper derives the exact solutions for a point mass moving along a stretched infinite string on a Winkler foundation at a constant velocity. The solutions for the contact force between the string and the mass are derived and then the displacement responses of the string can be obtained easily. The solutions cover infinite string subjected to a moving mass at subsonic, sonic, or supersonic velocities. When time tends to infinity, the asymptotical solutions for the contact force between the mass and the string and for the displacement of the contact point are derived. The formulas derived are shown to be correct by comparison with the semianalytical method.


---

## Body

## 1. Introduction

The dynamic response of a structure due to moving forces or masses is of importance in engineering applications, such as the vehicle and track interaction, the train and bridge interaction, and the pantograph and catenary interaction [1–3]. For such problems, a lot of analytical and numerical research has been achieved.By using the Laplace transformation, an analytical solution for a stretched infinite string subjected to a moving force at a constant speed was derived by Kanninen and Florence [4], and a set of three solutions was obtained for subsonic, sonic, and supersonic velocities. A stretched infinite string excited by a moving force with varying speed was investigated by Flaherty Jr. [5]. Using the Laplace transformation and Fourier transformation, a stretched semi-infinite string subjected to a moving force with a constant acceleration was considered by Sagartz and Forrestal [6]. Using the Laplace transformation, the dynamic response of an inhomogeneous stretched infinite string on an elastic foundation subjected to a moving force at a constant speed was investigated by Wolfert et al. [7]. An integral equation of the response of a stretched infinite string on an elastic foundation subjected to a moving force was obtained by Gavrilov [8], and an asymptotical solution was also presented. Applying the Fourier transformation, the dynamic response of a stretched infinite string on an elastically supported membrane subjected to a moving force was investigated by Dieterman and Kononov [9]. The dynamic response of two parallel stretched infinite strings on an elastically supported membrane subjected to two moving forces was obtained by Kononov and Dieterman [10]. A closed-form solution for the dynamic response of a stretched finite double-string which is connected by a Winkler elastic layer subjected to a moving harmonic force at a constant speed was presented by Rusin et al. [11].An analytical solution for a stretched finite string subjected to a mass moving at a constant speed was presented by Smith [12], but the mass of the string was neglected in this model. A numerical solution for a stretched semi-infinite string subjected to a constantly accelerating moving mass was derived by Rodeman et al. [13], and an asymptotical solution was also obtained when the effect of the inertia of the moving mass was small. An integral equation for the response of a stretched finite string subjected to a moving mass at a constant speed was obtained by Yang et al. [14], and then they used a numerical integration method to solve this equation.A space-time method for numerical solution of a stretched finite string subjected to a moving mass at a constant speed was proposed by Bajer and Dyniewicz [15], and then by using a combination of the Fourier integral transformation and time integration methods, a semi-analytical method for the same problem was proposed by Dyniewicz and Bajer [16]. A moving coordinate system was introduced by Kruse et al. [17] to obtain the eigenfrequencies of a stretched infinite string on a viscoelastic foundation subjected to a two-mass oscillator at a constant speed. The Lagrange multiplier method was adopted by Lee et al. [18] to describe the dynamic responses of a coupled moving mass-stretched beam with a separation between each mass. A stretched infinite string on an elastic foundation subjected to a constant accelerating moving mass was considered by Gavrilov and Indeitsev [19], and an integral equation to describe the interaction force between the mass and the string was given and solved numerically. Many numerical methods and semianalytical solutions to solve the problems in the vibrations of structures subjected to moving loads or masses were presented by Bajer and Dyniewicz [3].A stretched string on a Winkler foundation subjected to a moving point mass is a typical and important model to simulate the dynamic interaction between a structure and a moving mass. The exact solutions can be used as models for physical experiments and as benchmarks for testing numerical methods, and so forth. However, to the authors' knowledge, there is no exact solution for the problem of a stretched string on a Winkler foundation subjected to a moving point mass. In this paper, the exact solution for a stretched infinite string on a Winkler foundation subjected to a moving point mass is considered. The solutions cover the situation of the mass moving at subsonic, sonic, or supersonic velocities. The asymptotical solutions for the contact force between the mass and the string and for the displacement of the contact point are derived.


## 2. The Governing Equation for the Moving Mass System

As shown in Figure1, let us consider a point mass moving along an infinite string on a Winkler foundation. The line mass density and tensile force of the string are ρ and T, respectively, the coefficient of elasticity of the foundation is k, and the mass of the moving mass is m. The mass moves in the positive x direction starting from the original point O at a constant speed V.Figure 1
A point mass moving along a stretched string on a Winkler foundation.Assume that the contact force between the mass and the string is denoted byp
(
t
), and then the interaction between the mass and the string can be described as in Figure 2. Therefore, the motion equation of the mass is given by
(1)
m
y
¨
=
-
p
(
t
)
-
m
g
,
y
(
0
)
=
y
˙
(
0
)
=
p
(
0
)
=
0
,
where y is the displacement of the mass and g denotes the acceleration due to gravity. The motion equation and the initial conditions of the string are
(2)
ρ
∂
2
u
∂
t
2
=
T
∂
2
u
∂
x
2
-
k
u
+
p
(
t
)
δ
(
x
-
V
t
)
,
(3)
u
(
x
,
0
)
=
∂
u
∂
t
(
x
,
0
)
=
0
,
where δ denotes the Dirac delta function. The continuous condition for the mass and the string is
(4)
y
(
t
)
=
u
(
V
t
,
t
)
.
Therefore, (1) and (4) yield
(5)
m
d
2
u
(
V
t
,
t
)
d
t
2
=
-
p
(
t
)
-
m
g
,
p
(
0
)
=
0
.
Then, combining (2) with (5) gives the governing equation for the string subjected to a mass moving at a constant speed
(6)
ρ
∂
2
u
∂
t
2
=
T
∂
2
u
∂
x
2
-
k
u
-
m
g
δ
(
x
-
V
t
)
-
m
[
∂
2
u
(
x
,
t
)
∂
2
t
+
2
V
∂
2
u
(
x
,
t
)
∂
x
∂
t
+
V
2
∂
2
u
(
x
,
t
)
∂
x
2
]
x
=
V
t
in which the terms 2
m
V
∂
2
u
(
x
,
t
)
/
∂
x
∂
t
|
x
=
V
t and m
V
2
∂
2
u
(
x
,
t
)
/
∂
x
2
|
x
=
V
t are the Coriolis force and centripetal force, respectively.Figure 2
The interaction between the mass and the string.Equation (2) can also be written as
(7)
∂
2
u
∂
t
2
=
c
2
∂
2
u
∂
x
2
-
κ
2
u
+
ρ
-
1
p
(
t
)
δ
(
x
-
V
t
)
,
where
(8)
c
=
T
ρ
,
κ
=
k
ρ
in which c is the velocity of the wave. For different engineering applications, the value of c can be different. For example, if a track is considered, the magnitude of c is about 10
3
m
/
s, while if the contact wire of the pantograph-catenary system is considered, the magnitude of c is about 10
2
m
/
s. The solution for (2) and (3) can be given by [20, Section 7.4]
(9)
u
(
x
,
t
)
=
1
2
ρ
c
∫
0
t
p
(
τ
)
d
τ
∫
x
-
c
(
t
-
τ
)
x
+
c
(
t
-
τ
)
δ
(
σ
-
V
τ
)
vvvvvvvvvvvvvv
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
σ
)
2
]
d
σ
,
where J
0 denotes the Bessel function of the first kind of order 0, and
(10)
α
=
κ
c
.
When V
<
c, (9) can be written as
(11)
u
(
x
,
t
)
=
1
2
ρ
c
{
0
,
-
c
t
>
x
∫
0
(
x
+
c
t
)
/
(
c
+
V
)
p
(
τ
)
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
c
t
≤
x
≤
V
t
∫
0
(
x
-
c
t
)
/
(
V
-
c
)
p
(
τ
)
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
V
t
≤
x
≤
c
t
0
,
c
t
<
x
.
When V
>
c, (9) can be written as
(12)
u
(
x
,
t
)
=
1
2
ρ
c
{
0
,
-
c
t
>
x
∫
0
(
x
+
c
t
)
/
(
c
+
V
)
p
(
τ
)
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
c
t
≤
x
≤
c
t
∫
(
x
-
c
t
)
/
(
V
-
c
)
(
x
+
c
t
)
/
(
c
+
V
)
p
(
τ
)
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
c
t
≤
x
≤
V
t
0
,
V
t
<
x
and when V
=
c, (9) can be written as
(13)
u
(
x
,
t
)
=
1
2
ρ
c
{
0
,
-
c
t
>
x
∫
0
(
x
+
c
t
)
/
2
c
p
(
τ
)
×
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
c
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
c
t
≤
x
≤
c
t
0
,
c
t
<
x
.Combining (5) with (11)–(13) gives the governing equation of the contact force when the mass moves at the subsonic, supersonic, or sonic velocities, respectively. In the following sections, the analytical solution for the contact force is derived.


## 3. The Solutions When the Mass Moves at Subsonic Velocity

WhenV
<
c, using (11) the following equation can be obtained as
(14)
d
2
u
(
V
t
,
t
)
d
t
2
=
1
2
ρ
c
[
p
˙
(
t
)
+
β
∫
0
t
{
-
β
J
0
[
β
(
t
-
τ
)
]
+
J
1
[
β
(
t
-
τ
)
]
t
-
τ
}
p
(
τ
)
d
τ
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
+
J
1
[
β
(
t
-
τ
)
]
t
-
τ
}
p
(
τ
)
d
τ
]
,
where J
1 denotes the Bessel function of the first kind of order 1, and
(15)
β
=
α
c
2
-
V
2
.
Then, according to (5) and (14), the following differential equation can be obtained:
(16)
p
˙
(
t
)
+
β
∫
0
t
{
-
β
J
0
[
β
(
t
-
τ
)
]
+
J
1
[
β
(
t
-
τ
)
]
t
-
τ
}
p
(
τ
)
d
τ
=
-
2
ρ
c
m
p
(
t
)
-
2
ρ
c
g
.
Equation (16) gives the governing differential equation for contact force p
(
t
) when the mass moves at subsonic velocity; that is, V
<
c.Performing the Laplace transformation to (16) gives
(17)
(
m
2
ρ
c
s
2
s
2
+
β
2
+
1
)
p
~
(
s
)
-
m
2
ρ
c
p
(
0
)
+
m
g
s
=
0
,
where p
~
(
s
) denotes the Laplace transformation of p
(
t
). Because p
(
0
)
=
0, so
(18)
p
~
(
s
)
=
-
2
m
g
ρ
c
s
2
+
β
2
s
(
m
s
2
+
2
ρ
c
s
2
+
β
2
)
.
Equation (18) can also be written as
(19)
p
~
(
s
)
=
m
g
[
-
1
s
+
m
s
1
m
(
s
2
+
β
2
)
+
2
ρ
c
s
2
+
β
2
-
m
β
2
]
.
Therefore, the inverse Laplace transformation of (19) gives
(20)
p
(
t
)

=

-
m
g
+
m
2
g
L
-
1
[
s
m
(
s
2
+
β
2
)
+
2
ρ
c
s
2
+
β
2
-
m
β
2
]
,
where L
-
1
(
#
) denotes the inverse Laplace transformation of function #. According to the following property of Laplace transformation:
(21)
L
-
1
[
s
f
~
(
s
)
]
=
d
f
(
t
)
d
t
,
if
f
(
+
0
)
=
0
,
we have
(22)
L
-
1
[
s
1
m
(
s
2
+
β
2
)
+
2
ρ
c
s
2
+
β
2
-
m
β
2
]
=
d
d
t
L
-
1
[
1
m
(
s
2
+
β
2
)
+
2
ρ
c
s
2
+
β
2
-
m
β
2
]
.
Then, by using the following two equations [21, equation (5) in Section 5.1] and [22, equation (9a)]
(23)
L
-
1
[
f
~
(
s
2
+
a
2
)
]
=
f
(
t
)
-
a
∫
0
t
f
(
t
2
-
τ
2
)
J
1
(
a
τ
)
d
τ
,
L
-
1
[
1
m
s
2
+
2
ρ
c
s
-
m
β
2
]
=
1
ρ
2
c
2
+
m
2
β
2
exp
⁡
(
-
ρ
c
m
t
)
sinh
⁡
(
ρ
2
c
2
+
m
2
β
2
m
t
)
.
Equation (22) gives
(24)
L
-
1
[
1
m
(
s
2
+
β
2
)
+
2
ρ
c
s
2
+
β
2
-
m
β
2
]
=
1
σ
[
∫
0
t
exp
⁡
(
-
ρ
c
m
t
2
-
τ
2
)
exp
⁡
(
-
ρ
c
t
m
)
sinh
⁡
(
σ
t
m
)
-
β
∫
0
t
exp
⁡
(
-
ρ
c
m
t
2
-
τ
2
)
×
sinh
⁡
(
σ
m
t
2
-
τ
2
)
J
1
(
β
τ
)
d
τ
∫
0
t
exp
⁡
(
-
ρ
c
m
t
2
-
τ
2
)
]
,
where
(25)
σ
=
ρ
2
c
2
+
m
2
β
2
.
Finally, (20), (22), and (24) give
(26)
p
(
t
)
=
m
g
2
σ
{
+
(
σ
-
ρ
c
)
×
{
J
1
(
β
τ
)
exp
⁡
[
(
σ
-
ρ
c
)
m
t
2

-

τ
2
]
t
2
-
τ
2
exp
⁡
[
(
σ
-
ρ
c
)
m
t
]
-
β
t
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
(
σ
-
ρ
c
)
m
t
2

-

τ
2
]
t
2
-
τ
2
d
τ
}
+
(
σ
+
ρ
c
)
×
{
J
1
(
β
τ
)
exp
⁡
[
(
σ
-
ρ
c
)
m
t
2

-

τ
2
]
t
2
-
τ
2
exp
⁡
[
-
(
σ
+
ρ
c
)
m
t
]
-
β
t
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
(
σ
+
ρ
c
)
m
t
2

-

τ
2
]
t
2
-
τ
2
d
τ
}
}
-
m
g
.Equation (25) shows that σ
-
ρ
c
>
0 and so the term exp
⁡
[
(
σ
-
ρ
c
)
t
/
m
] in (26) tends to infinity when time tends to infinity, which would cause numerical problem in computation. By using
(27)
exp
⁡
[
(
σ
-
ρ
c
)
m
t
]
=
2
cosh
⁡
[
(
σ
-
ρ
c
)
m
t
]

-

exp
⁡
[
-
(
σ
-
ρ
c
)
m
t
]
.
Equation (26) can be written as
(28)
p
(
t
)
=
m
g
2
σ
{
+
(
σ
-
ρ
c
)
×
{
exp
⁡
[
(
σ
-
ρ
c
)
m
t
]
-
2
β
t
×
∫
0
t
J
1
(
β
τ
)
cosh
⁡
[
m
-
1
(
σ

-

ρ
c
)
t
2

-

τ
2
]
t
2
-
τ
2
d
τ
}
+
(
σ
+
ρ
c
)
exp
⁡
[
-
(
σ
+
ρ
c
)
m
t
]
+
(
σ
-
ρ
c
)
β
t
×
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
(
σ
-
ρ
c
)
m
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
-
(
σ
+
ρ
c
)
β
t
×
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
(
σ
+
ρ
c
)
m
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
}
-
m
g
.
Since [22, equation (3) in Section 6.737 and equation (1) in Section 8.464]
(29)
2
t
β
∫
0
t
J
1
(
β
τ
)
cosh
⁡
[
m
-
1
(
σ
-
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
=
t
β
π
J
0.5
[
t
2
(
β
2
-
m
-
2
(
σ
-
ρ
c
)
2
-
i
m
-
1
(
σ
-
ρ
c
)
)
]
×
J
0.5
[
t
2
(
β
2
-
m
-
2
(
σ
-
ρ
c
)
2
+
i
m
-
1
(
σ
-
ρ
c
)
)
]
=
exp
⁡
[
t
m
-
1
(
σ
-
ρ
c
)
]
+
exp
⁡
[
-
t
m
-
1
(
σ
-
ρ
c
)
]
-
2
cos
⁡
(
ω
t
)
,
where i
=
-
1 and
(30)
ω
=
2
ρ
c
(
σ
-
ρ
c
)
m
,
therefore (28) and (29) give
(31)
p
(
t
)
=
m
g
[
σ
-
ρ
c
σ
cos
⁡
(
ω
t
)
-
1
]
+
g
(
t
)
,
where
(32)
g
(
t
)
=
m
g
2
σ
{
+
(
σ
-
ρ
c
)
β
t
×
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
m
-
1
(
σ
-
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
-
(
σ
+
ρ
c
)
β
t
×
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
m
-
1
(
σ
+
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
+
(
σ
+
ρ
c
)
exp
⁡
(
-
σ
+
ρ
c
m
t
)
-
(
σ
-
ρ
c
)
exp
⁡
(
-
(
σ
-
ρ
c
)
m
t
)
}
.
Equation (31) gives the analytical solution for the contact force between the mass and the string when the mass moves at subsonic velocities. Once the contact force is obtained, the displacement of the string can be obtained using (11).Since the string is infinite, it is worthwhile to see what happens to the contact force when the time tends to infinity. By using variable substitution, the following equation can be obtained as(33)
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
m
-
1
(
σ
-
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
=
∫
0
π
/
2
J
1
[
β
t
cos
⁡
(
θ
)
]
exp
⁡
[
-
m
-
1
(
σ
-
ρ
c
)
t
sin
(
θ
)
]
d
θ
.
By using |
J
1
(
x
)
|
≤
1
/
x, x
≥
0,β
>
0 and m
-
1
(
σ
-
ρ
c
)
>
0, we have
(34)
|
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
γ
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
|
≤
∫
0
π
/
2
exp
⁡
[
-
γ
t
sin
(
θ
)
]
β
t
cos
⁡
(
θ
)
d
θ
=
φ
(
t
)
,
where
(35)
φ
(
t
)
=
1
β
t
[
-
4
γ
t
5
-
3
π
3
/
2
2
1
/
4
L
3
/
4
(
γ
t
)
4
(
γ
t
)
3
/
4
Γ
(
3
/
4
)
+
π
3
/
2
(
2
γ
t
)
1
/
4
2
Γ
(
3
/
4
)
[
I
-
1
/
4
(
γ
t
)
-
L
7
/
4
(
γ
t
)
]
]
,
γ
=
m
-
1
(
σ
-
ρ
c
)
,
in which I
-
1
/
4
(
γ
t
) denotes the modified Bessel function of the first kind of order -
1
/
4, L denotes the modified Struve function, and Γ denotes the Gamma function. By using Maple software, it is easy to prove that
(36)
lim
⁡
t
→
∞
[
t
φ
(
t
)
]
=
0
.
Therefore, (33), (34), and (36) give
(37)
lim
⁡
t
→
∞
|
t
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
m
-
1
(
σ
-
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
|
=
0
.
In a similar way, we can prove that
(38)
lim
⁡
t
→
∞
|
t
∫
0
t
J
1
(
β
τ
)
exp
⁡
[
-
m
-
1
(
σ
+
ρ
c
)
t
2
-
τ
2
]
t
2
-
τ
2
d
τ
|
=
0
.
By using (32), (37), and (38), we can obtain
(39)
lim
⁡
t
→
∞
g
(
t
)
=
0
.
Then, (31) and (39) give
(40)
lim
⁡
t
→
∞
p
(
t
)
=
m
g
[
σ
-
ρ
c
σ
cos
⁡
(
ω
t
)
-
1
]
.
Equation (40) shows that when the mass moves at subsonic velocities and when the time tends to infinity, there exists an asymptotical solution for the contact force.Furthermore, we consider the asymptotical solution for the displacement of the contact point between the string and the mass. According to (11), when x
=
V
t, the displacement of the contact point is given by
(41)
u
(
V
t
,
t
)
=
1
2
ρ
c
∫
0
t
p
(
τ
)
J
0
[
β
(
t
-
τ
)
]
d
τ
.
A combination of (31) and (41) gives
(42)
u
(
V
t
,
t
)
=
1
2
ρ
c
{
m
g
(
σ
-
ρ
c
)
σ
∫
0
t
[
g
(
τ
)
-
m
g
]
J
0
[
β
(
t
-
τ
)
]
d
τ
+
m
g
(
σ
-
ρ
c
)
σ
∫
0
t
cos
⁡
[
ω
(
t
-
τ
)
]
J
0
(
β
τ
)
d
τ
}
.Firstly, the first term of the right-hand of (42) is considered, which is the convolution of g
(
t
) and J
0
(
β
t
). The limit of the convolution of two functions has the following property [23]; that is, let f
,
h
[
0
,
∞
)
→
R be locally integrable and assume that
(43)
f
∞
=
lim
⁡
t
→
∞
f
(
t
)
,
∫
0
∞
|
g
(
t
)
|
d
t
<
∞
,
then
(44)
lim
⁡
t
→
∞
∫
0
t
h
(
τ
)
f
(
t
-
τ
)
d
τ
=
f
∞
·
∫
0
∞
h
(
τ
)
d
τ
.
By using (39) and (44) and the following equation:
(45)
∫
0
∞
J
0
(
β
τ
)
d
τ
=
1
β
,
we can obtain
(46)
lim
⁡
t
→
∞
∫
0
t
[
g
(
τ
)
-
m
g
]
J
0
[
β
(
t
-
τ
)
]
d
τ
=
-
m
g
β
.The second term of the right-hand of (42) can be written as
(47)
∫
0
t
cos
⁡
[
ω
(
t
-
τ
)
]
J
0
(
β
τ
)
d
τ
=
cos
⁡
(
ω
t
)
∫
0
t
cos
⁡
(
ω
τ
)
J
0
(
β
τ
)
d
τ
+
sin
(
ω
t
)
∫
0
t
sin
(
ω
τ
)
J
0
(
β
τ
)
d
τ
.
Since ω can also be written as ω
=
β
2
-
m
-
2
(
σ
-
ρ
c
)
2, so β
>
ω
>
0, and then the table of integrals gives [22, equations (7) and (8) In Section 6.671]
(48)
lim
⁡
t
→
∞
∫
0
t
cos
⁡
(
ω
τ
)
J
0
(
β
τ
)
d
τ
=
1
β
2
-
ω
2
,
lim
⁡
t
→
∞
∫
0
t
sin
(
ω
τ
)
J
0
(
β
τ
)
d
τ
=
0
,
and (47) and (48) give
(49)
lim
⁡
t
→
∞
∫
0
t
cos
⁡
[
ω
(
t
-
τ
)
]
J
0
(
β
τ
)
d
τ
=
cos
⁡
(
ω
t
)
β
2
-
ω
2
.A combination of (42), (46), and (49) gives
(50)
lim
⁡
t
→
∞
u
(
V
t
,
t
)
=
m
g
2
ρ
c
[
(
σ
-
ρ
c
)
β
2
-
ω
2
cos
⁡
(
ω
t
)
-
1
β
]
.
Therefore, when the mass moves at subsonic velocities and when the time tends to infinity, there exists an asymptotical solution for the displacement of the contact point.


## 4. The Solutions When the Mass Moves at Supersonic Velocity

WhenV
>
c, (12) gives
(51)
d
2
u
(
V
t
,
t
)
d
t
2
=
0
.
Then, (5) gives
(52)
0
=
-
p
(
t
)
-
m
g
,
which yields
(53)
p
(
t
)
=
-
m
g
.
Therefore, when V
>
c, (12) and (53) give
(54)
u
(
x
,
t
)
=
-
m
g
2
ρ
c
{
0
,
-
c
t
≥
x
∫
0
(
x
+
c
t
)
/
(
c
+
V
)
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
c
t
<
x
≤
c
t
∫
(
x
-
c
t
)
/
(
V
-
c
)
(
x
+
c
t
)
/
(
c
+
V
)
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
x
-
V
τ
)
2
]
d
τ
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
c
t
≤
x
<
V
t
0
,
V
t
≤
x
.
Equations (53) and (54) give the analytical contact force and displacement solution of the string, respectively.Next, the displacement solution for the string at special locations is considered. Whenx
=
V
t, (54) gives
(55)
u
(
V
t
,
t
)
=
0
.
So, the displacement of the contact point is always zeros when the mass moves at supersonic velocities.Equations (53) and (55) can be explained as follows. For the case of the supersonic speeds, the velocity of the moving mass is greater than the wave speed, so the string in front of the mass is always still and so the contact force between the mass and the string is -
m
g.Whenx
=
c
t, (54) gives
(56)
u
(
c
t
,
t
)
=
-
m
g
2
ρ
c
∫
0
2
c
t
/
(
c
+
V
)
J
0
[
α
c
2
(
t
-
τ
)
2
-
(
c
t
-
V
τ
)
2
]
d
τ
.
By using variable substitution, (56) can be rewritten as
(57)
u
(
c
t
,
t
)
=
-
m
g
t
2
ρ
(
c
+
V
)
∫
0
π
J
0
[
α
c
t
V
-
c
c
+
V
sin
(
τ
)
]
sin
(
τ
)
d
τ
.
Then, by using the following equation [22, equation (9) in Section 6.681]:
(58)
∫
0
π
J
υ
[
2
a
sin
(
y
)
]
sin
(
2
μ
y
)
d
y
=
π
sin
(
μ
π
)
J
υ
-
μ
(
a
)
J
υ
+
μ
(
a
)
,
we have
(59)
u
(
c
t
,
t
)

=

-
m
g
2
ρ
π
t
c
+
V
J
-
0.5

[
α
c
t
2
V
-
c
c
+
V
]
J
+
0.5
[
α
c
t
2
V
-
c
c
+
V
]
.
Furthermore, using the following equation:
(60)
J
0.5
(
z
)
=
sin
(
z
)
2
π
z
,
J
-
0.5
(
z
)
=
cos
⁡
(
z
)
2
π
z
,
(59) becomes
(61)
u
(
c
t
,
t
)
=
-
m
g
ρ
α
c
V
2
-
c
2
sin
(
α
c
t
V
-
c
c
+
V
)
.


## 5. The Solutions When the Mass Moves at Sonic Velocity

WhenV
=
c, (13) gives
(62)
d
2
u
(
V
t
,
t
)
d
t
2
=
1
2
ρ
c
p
˙
(
t
)
,
and then (5) gives
(63)
p
˙
(
t
)
=
-
2
ρ
c
m
p
(
t
)
-
2
ρ
c
g
.
The solution of (63) can be given by
(64)
p
(
t
)
=
m
g
[
exp
⁡
(
-
2
ρ
c
m
t
)
-
1
]
.
Therefore, when the mass moves at sonic velocity, the contact force is given by a very simple formula.Substituting (64) into (13) and performing some variable substitutions give the following equation:
(65)
u
(
x
,
t
)
=
m
g
2
ρ
c
{
0
,
-
c
t
>
x
[
c
2
t
2
-
x
2
J
1
(
α
c
2
t
2
-
x
2
)
c
α
(
c
t
-
x
)
x
+
c
t
2
c
∫
0
1
exp
⁡
(
-
ρ
(
x
+
c
t
)
m
τ
)
×
J
0
[
α
c
2
t
2
-
x
2
1
-
τ
]
d
τ
-
c
2
t
2
-
x
2
J
1
(
α
c
2
t
2
-
x
2
)
c
α
(
c
t
-
x
)
]
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
c
t
≤
x
≤
c
t
0
,
c
t
<
x
.
Then, using (65), it is easy to obtain the following two equations:
(66)
lim
⁡
x
→
c
t
-
u
(
x
,
t
)
=
g
(
2
ρ
c
m
)
-
2
[
1
-
exp
⁡
(
-
2
ρ
c
m
t
)
-
2
ρ
c
m
t
]
,
lim
⁡
x
→
c
t
+
u
(
x
,
t
)
=
0
.
Therefore, the displacement of the contact point is discontinuous when the mass moves at sonic velocity.


## 6. Numerical Results

In this section, some numerical results for a point mass moving along an infinite string on a Winkler foundation are presented. If we letτ
=
c
t, then (2) becomes
(67)
∂
2
u
∂
τ
2
=
∂
2
u
∂
x
2
-
k
T
u
+
1
T
p
(
t
)
δ
(
x
-
V
c
τ
)
.
Equation (67) shows that the key parameter is the ratio of the velocity of the mass to the wave speed instead of their absolute values. Therefore, we pay our attention to the different values of the ratio of the velocity of the mass to the wave speed and the following parameters are used. The line mass density and tensile force of the string are ρ
=
1
(
kg
/
m
) and T
=
1
(
N
), respectively; the coefficient of elasticity of the foundation is k
=
10
(
N
/
m
2
); the mass of the moving mass is m
=
2
(
kg
); and the acceleration due to gravity is g
=
10
(
m
/
s
2
). The mass moves in the positive x direction starting from the original point O at a constant speed V; see Figure 1. For the above parameters, the velocity of the wave in the string is 1
(
m
/
s
).To demonstrate the validity of the proposed formula, a comparison with the semianalytical method based on the Fourier analysis is presented. The semianalytical method based on the Fourier analysis [3, Section 3.1.1] reduces the partial differential equation to an ordinary differential equation by applying the Fourier analysis, and the ordinary differential equation is then solved numerically. The accuracy of the semianalytical method depends on both the number of Fourier series and numerical method used to solve the ordinary differential equation. The semianalytical results are obtained for the string with finite length 20 (m) and the mass moving at velocity V
=
0.4
(
m
/
s
) from the middle point of the string to three-quarter length of the sting. In this situation, the semianalytical results are not affected by the string end conditions. The contact force is shown in Figure 3, in which the black, blue, and olive lines correspond to the results obtained by using the semianalytical method with 800, 2000, and 3200 Fourier series, respectively, and the red line corresponds to the results obtained using the proposed formula. Figure 3 shows that as the number of Fourier series increases, the results obtained using the semianalytical method become closer to those obtained from the proposed method.Figure 3
The contact force obtained from the proposed formula and the semianalytical method with different numbers of Fourier series whenV
=
0.4.The contact forces between the mass and the string forV
=
0.4, 0.6, 0.8, 1, 2, and 3
(
m
/
s
) are given in Figure 4, which shows that the contact force for supersonic velocities is a constant m
g (as shown in Section 4); the contact force for sonic velocity decays exponentially and tends to m
g, while the contact force for subsonic velocities is oscillating. As shown in Section 3, for subsonic velocities, there exists an asymptotical solution for the contact force when the time tends to infinity. Figure 5 gives the difference d
p between the analytical and asymptotical contact forces and shows that the analytical contact force tends to the asymptotical solution when time tends to infinity.Figure 4
The contact force between the mass and the string.Figure 5
The difference between the analytical and asymptotical contact forces.The displacements of the contact point forV
=
0.4, 0.6, 0.8, 1, 2, and 3
(
m
/
s
) are given in Figure 6, which shows that the displacement of the contact point for supersonic velocities is zero, the displacement of the contact point for sonic velocity increases exponentially, while the displacement of the contact point for subsonic velocities is oscillating. As shown in Section 3, for subsonic velocities, the displacement of the contact point has an asymptotical solution when time tends to infinity. Figure 7 gives the difference d
u
c between the analytical and asymptotical displacements of the contact point and shows that the analytical displacement of the contact point tends to the asymptotical solution when time tends to infinity. For supersonic velocities (i.e., V
=
2 and 3), the displacement of the string at x
=
c
t is given in Figure 8.Figure 6
The displacement of the contact point.Figure 7
The difference between the analytical and asymptotical displacements of the contact point.Figure 8
The displacement of the string atx
=
c
t when the mass moves at supersonic velocities.The displacements of the string when the mass moves at subsonic velocities (i.e.,V
=
0.4, 0.6 and 0.8) are given in Figure 9, and the displacements of the string when the mass moves at sonic (i.e., V
=
1) and at supersonic velocities (i.e. V
=
2 and 3) are given in Figure 10.The displacement of the string when the mass moves at subsonic velocities (a)V
=
0.4; (b) V
=
0.6; and (c) V
=
0.8.
(a)
(b)
(c)The displacement of the string when the mass moves at sonic or supersonic velocities (a)V
=
1; (b) V
=
2; and (c) V
=
3.
(a)
(b)
(c)


## 7. Conclusions

This paper has presented exact solutions for a point mass moving along a stretched string on a Winkler foundation. The exact solution of the contact force between the string and the mass has been derived, and then the displacement response of the string is given. The solutions cover infinite string subjected to a moving mass at subsonic, sonic, or supersonic velocities. When the mass moves at supersonic velocities, the contact force between the mass and the string isp
(
t
)
=
-
m
g and the displacement of the contact point is always zero. When the mass moves at sonic velocity, the contact force between the mass and the string is p
(
t
)
=
m
g
[
exp
⁡
(
-
2
ρ
c
t
/
m
)
-
1
] and the displacement of the contact point is discontinued. When the mass moves at subsonic velocities, the contact force and the displacement of the contact point are complicated, while their asymptotical solution when time tends to infinity can be given as lim
⁡
t
→
∞
⁡
p
(
t
)
=
m
g
[
(
(
σ
-
ρ
c
)
/
σ
)
cos
⁡
(
ω
t
)
-
1
] and lim
⁡
t
→
∞
u
(
V
t
,
t
)
=
(
m
g
/
2
ρ
c
)
[
(
(
σ
-
ρ
c
)
/
β
2
-
ω
2
)
cos
⁡
(
ω
t
) − 1
/
β
], respectively.



---

*Source: 136149-2014-06-30.xml*