Advances in Applied Mechanics T9 - (1966)
Advances in Applied Mechanics T9 - (1966)
P. G. DRAZIN
L. N. HOWARD
N. N. MOISEEV
A. A. PETROV
PIOTRPBRZYNA
R. M. ROSENBERG
ADVANCES IN
APPLIED MECHANICS
Edited by
G. G. CHERNYI W. OLSZAK
H. L. DRYDEN W. PRAGER
P. GERMAIN R. F. PROBSTEIN
L. HOWARTH H. ZIEGLER
Managing Editor
G. KUERTI
Case Institrzte of Technology, Cleveland, Ohio
VOLUME 9
1966
LIBRARY
OF CONGRESS CATALOG CARD 48-8503
NUMBER:
V
In preparing this ninth volume of the “Advances in Applied Mechanics”
we still had the invaluable aid of Dr. Dryden who, unfortunately, died on
December 2, 1966. His advice and judgment will be sorely missed in the
future.
Stability theory in two different fields of mechanics is the subject of two
contributions in this volume. The volume also contains a detailed report on
the sloshing of liquid in containers; the basic analysis of the sloshing problem
was set forth in Volume 8. A survey of the present theories of elastic-visco-
plastic behavior concludes the volume.
It is intended to publish at least the next volume of the “Advances” in
the form of several successive fascicles, in order to present the material as
rapidly as possible.
THE EDITORS
May, 1966
vii
Hydrodynamic Stability of Parallel Flow of Inviscid Fluid
BY P . G. DRAZIN AND L. N . HOWARD
Page
I . Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
I1 . Inertial Instability . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1. Eigenvalue Problem for Inertial Modes . . . . . . . . . . . . . . . . 3
2. General Stability Characteristics of Plane Parallel Flow . . . . . . . . 10
3 . The Initial-Value Problem and the Stability of Non-parallel Flow . . . . 22
4 . Stability Characteristics of Various Basic Flows . . . . . . . . . . . . 32
I11. Waves and Stability of Plane Parallel Flow of Inviscid Fluid under the Actions
of Various Force Fields . . . . . . . . . . . . . . . . . . . . . . . 43
1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
2. Internal Gravity Waves and Stability of a Fluid of Variable Density . . . 44
3. Sound Waves and Stability of Compressible Fluid . . . . . . . . . . . 47
4 . Planetary Waves and Stability in a Rotating System . . . . . . . . . 49
5 . Rossby Waves and Stability of Fluid in a Rotating System with Variable
Coriolis Parameter . . . . . . . . . . . . . . . . . . . . . . . . . 51
6. Magnetohydrodynamic Waves and Stability of an Electrically-Conducting
Fluid in a Magnetic Field . . . . . . . . . . . . . . . . . . . . . 62
IV . Heuristic Theory of Instability . . . . . . . . . . . . . . . . . . . . 54
1. Dimensional Analysis . . . . . . . . . . . . . . . . . . . . . . . 54
2. Physical Arguments . . . . . . . . . . . . . . . . . . . . . . . . 57
V. Instability of an Incompressible Fluid of Variable Density . . . . . . . . 80
1. General Stability Characteristics . . . . . . . . . . . . . . . . . . . 60
2 . Stability Characteristics of Various Basic Flows . . . . . . . . . . . . 68
VI . Stability of Other Parallel Flows . . . . . . . . . . . . . . . . . . . 79
1. Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
2 . The Semicircle Theorem for General Parallel Flow . . . . . . . . . . . 79
3. Inertial Instability of Axisymmetric Jets . . . . . . . . . . . . . . . 80
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
I . INTRODUCTION
Stability of parallel flow of inviscid fluid was first studied in the last
third of the nineteenth century. notably by Helmholtz [l]. Kelvin [2] and
Rayleigh [S] . They considered the inertial instability of a homogeneous
incompressible fluid. and its modification- Kelvin-Helmholtz instability-
when there is variation of density of the fluid transverse to the basic flow.
Subsequent authors have continued this work and gone on to consider
1
2 P. G. DRAZIN AND L. N. HOWARD
in some cases may be stabilized by the force field. We add a survey of results
and literature for each force field. The fields we consider are compressibility
in a fluid of variable temperature, buoyancy due to variations of density,
Coriolis force due to rotation of the system in which the parallel flow is
placed, variations of this Coriolis force transverse to the flow, and magneto-
hydrodynamic forces.
The similarity of all these problems is brought out in Section IV. Mech-
anisms of instability are discussed, and analyzed dimensionally to give some
general stability characteristics. Some of these dimensional arguments are
elaborated by physical ones.
Lack of space and time prevent our treating in detail the case of each
force field, so we have picked the single case of buoyancy due to variations
of density for detailed study in Section V. This case is as typical as any,
and has the advantages of practical importance and of advanced theoretical
development. We discuss this case in Section V much as we did inertial
instability in Section 11. Finally, in Section VI, we give some results on non-
planar parallel flows.
11. INERTIALINSTABILITY
I
lb1
A wide range of problems has been solved since, the theory has been extended
to viscous fluids, and applications of hydrodynamic stability are numerous.
This abundance of work makes a chronological account impractical, so we
develop the subject logically, referring to authors where appropriate. Strictly
speaking, a logical account should begin with the formulation of the problem
for an arbitrary, not necessarily parallel, basic flow. However, since the
vast majority of work in stability theory has been on the parallel flow case
and only a few general results are known otherwise, we shall defer our remarks
on non-parallel flow to Section 11.3, and begin here in the traditional manner.
We consider the stability of a basic plane parallel flow of inviscid fluid
with given velocity
fi+ = (@+(Y+)~O~O) cYi+ <Y+ <YaJ.
This is illustrated in Figure 1. We take the flow bounded by the two planes
y+ = ylS,yaI parallel to the flow. Each of these planes is either rigid or free,
so that either the normal velocity of the fluid is zero or the pressure con-
stant there. However, we allow one or both of the planes to be at infinity.
In this section we suppose the fluid is incompressible and homogeneous.
It is convenient to choose some velocity scale V of the basic flow zi,(y,)
and some length scale L of its transverse variations in order to introduce
dimensionless variables. Thus we define in the usual way the dimensionless
time, position vector, velocity, basic velocity and pressure as
t= t,V/L, r =r+/L u =u,/V, = ii,,/V,
ii = (6,0,0)
P =P*/PVa
respectively, where p is the density of the fluid. Now the Euler equations
of motion can be written
+
aulat U . PU = - vp.
Also the equations of continuity and incompressibility give
P*u=O.
The basic flow satisfies these equations and the boundary conditions
with uniform pressure fi. To study the instability of this flow one puts
+
u(r.4 = i(y) u'(r,t), $(r,t) = p'(r4 +
and neglects quadratic terms in the perturbations, denoted by primes.
Linearizing the equations of motion in this way for small perturbations,
we find
+ +
&#/at aa%#lax virtc2/dy = - apllax,
+
a q a t navllax = - apllay,
awllat + nawi/ax = - api/az,
+
atqax avl/ay awi/az = 0.+
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 6
for some functions i, 3, where a,y are real wave-numbers and c = c, ic, +
is a complex velocity. I t is to be understood that real parts are to be taken
to get physical quantities, this being permissible for the linear problem.
Thus each component travels as a wave in the direction of (a,O,y) with
+
phase speed ac,/(a2 y2)1/2and grows or dies away in time like exp (acit).
Therefore the wave is unstable when aci > O and stable when aci Q 0. It
is said to be neutrally stable when aci = 0.
This assumption that any disturbance can be represented by wave
components, according to the method of normal modes, serves to separate
the variables and reduce the linearized equations of motion from partial
to ordinary differential equations. They now become
for some complex constants A,, A,. Substitution into boundary conditions
(2.12) and elimination of A , , A , gives the eigenvalue relation
Then
V’ = Dq/Dt = aq/at + waq/ax = ia(w - c)V,
to first degree in the perturbation. But v = iav, so it follows that
(2.14) [d(W - 41 = 0
at y = yo.
Conditions (2.13), (2.14) can alternatively be proved mathematically.
We may take yo = 0 without loss of generality. Then stability equation
(2.11) has integral
[ ( w - c ) D v - ( D W ) ~ ]=? a2
1
~ (w - c)pldy.
-8
[ d ( w - 4Ih. = a2
i
-e
dy(w -
-e idyl(wl - c ) ~ , ,
The fundamentals of this section are due to Rayleigh [3]. For more
recent accounts, the book of Lin [9, Chaps. 4,8] and chapter of Stuart [lo]
are recommended. These two works, and Reid's survey [Ill, also place
the inviscid theory in its viscous context.
(2.16)
(2.17)
If c, > O it now follows that Dew must change sign at one or more points
in the field of flow. On assuming that Dew is continuous, there must be at
least one inflection point on (y l ,y2)and indeed an inflection point a t which
the velocity profile crosses its tangent, i.e. a relative maximum or minimum
in the basic oorticity i3 = - Dw. With the (weaker) assumption that i3
is only piecewise continuously differentiable (which has really been tacitly
made anyway in writing down the stability equation) we can still say that
6 must have a relative maximum or minimum.
In 1960 Fjmtoft [13, p. 261 proved the stronger necessary condition
for instability that (D%)(w - tos)< 0 somewhere in the field of flow, where
ys is a point at which D f vanishes, and where w, = w(y,). A proof comes
from the real part of equation (2.16),
Adding
J
YI
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 11
FIG2. (a) Stable: D*w < 0. (b) Stable: Dgw> 0. (c) Stable: D’w = 0 at ys, but Dw in-
creases where w < ws. (d) Possibly unstable: D*w = 0 at ys and (w - ws)Dsw< 0 elsewhere.
(2.18)
the minimum being for functions f that vanish a t the walls and have square-
integrable derivatives. Therefore a neutral eigensolution with positive
a = a, = (- As)1'2 exists if and only if I, < 0. There may be a finite number
of other eigensolutions for larger eigenvalues I, provided these are negative.
Also there may be other series of eigensolutions when c = w, for other
values of w, a t other points of inflection, and sometimes eigenfunctions with
real c not equal to the value of w at an inflection point, though these have
slightly singular behavior.
The existence of the neutral eigensolution with I, < 0 follows easily
when K ( y ) > n 2 / ( y 2- y J 2 everywhere on account of the well-known in-
equality,
i
Yn
Again, when K ( y ) > O over the field of flow and w vanishes a t the walls
but not between, trial of / = w shows that A,< 0.
Tollmien also demonstrated heuristically the existence of unstable waves,
whose limit as ci + 0 is the neutral s-wave above. This has been considered
alternatively by Lin [7, pp. 223-224; 9, pp. 122-1231 as follows. The
stability equation for Q, gives
(2.19) D*p,, + {I, - (D2w)/(w- w,)}p,, = 0.
Multiplying this by p, and subtracting vS times the equation for Q,we get
D(Q& - f@pls) - ( A - Is)Q% - ( o w { ( w - c)-l - (w - ws)-'}plqJ, = 0.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 13
YI I1
-
To find the unstable solutions near the neutral one, we take the limits
a +as, c w,, q +vs. Then
(w - wJ2
Yl Y.
(2.20)
For a known neutral eigensolution this gives eigenvalues c(a)near wsin the
complex c-plane and thence the associated eigenfunctions in the limit
ci -+ 0 f . If K ( y ) > 0 throughout the field of flow, then the imaginary
part of the denominator is positive, and it follows that (dci/daa)a=. as < 0,
with instability for a just less than a,.
Further, it can be proved as follows that there is instability only when
a < as. Suppose K ( y ) > O throughout the field of flow. Then, when cj # 0,
the real part of equation (2.16) plus (w, - c,)/ci times equation (2.17)
gives
14 P. G. DRAZIN AND L. N. HOWARD
Therefore
5
Ys
YI
(w - c)((DylS + a2(ollS) + ( D W ( q I 2- y * ( D 4 ( D q ) d y= 0.
The imaginaiy part of this gives
(2.21) ci l
YS
YI
+
~Dvla a2lqlady =
l !
(Dw)(vDv*- pl*Dq)dy.
YI Y1
Foote and Lin [18] noted that the average of the Reynolds shear stress
over a wavelength
5
errla
1
= - &If = - (a/2n) u'vfdx = -ia(qDy* - y*Drp)e*it.
4
n
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 16
-.
modes, that is, suppose we have an unstable mode for some a, and as a
say, the corresponding c approaches a real value: ci(a) 0. Then for this
-.
as,
neutral mode, (2.22) shows that D t must be zero everywhere, except possibly
at y = y,(w(y,) = c) where the limit of the right-hand side of (2.22) does
not exist. By consideration of the nature of the singularity in D t which
appears as c, --r 0 + it can be shown that at the “critical layer” y = yc (or
layers) z has (in the limit ci --+ 0 -1) a jump [t],of magnitude
ideas to the instability of wind whereby ocean waves are generated.) Taylor
went on to note that viscosity allowed momentum to be diffused from the
boundaries, and suggested that a given basic flow might thus be stable for
inviscid but unstable for viscous fluid. This suggestion has since been veri-
fied, for plane Poiseuille flow as an example.
Lin [7, pp. 226-2271 also has interpreted physically the mechanism of
inertial instability by consideration of the migration of vorticity. He regarded
the flow due to a neutrally-stable disturbance in Kelvin’s “cat’s-eye” di-
(2.24) +
(w - c)s{1DFl2 aalFla)dy = 0.
YI
1
YI
But
j.1
YI
UQdy = &dy,
YI r j .
Y,
dQdy =
Y1
(Q' + c,')Qdy.
0 2
1 YI
(W - wmin)(w - wmax)Qdy
-j
- { (ci2
YI
+ G') - (wmin + wmax)C, + wminwmax}Qdy#
the maximum and minimum being taken over the field of flow yl < y < y2.
Therefore
cia + cva - (wmin + ~max)cr+ w&wmx <0,
i.e. for unstable waves c lies in the semicircle
(2.26) {c, - +(wmin + wmax)}' + ci2 < {i(wmax - wmin))' (ci> 0)-
18 P. G. DRAZIN AND L. N. HOWARD
This shows that any eigenvalue c, real or complex, must lie in or on the
circle with center t(w,, +
wmin) and radius i(wmx - wmin).
Again, with G = q/(w. - c)Y2, the stability equation can be written
(2.26) D((w - c)DG} - (40% + U'(W - C) + )(Dw)*/(w- c)}G = 0.
This has an integral
j.
Y.
+ aalG12}+ &(D2w)IG/'+ ) ( D W ) -~ (c+)lC/(w
(w - c)(lDGI2 ~ -c)l*}dy=O,
9 8 = Hvpcy) + d- Y))
and the odd part
vo = H d Y ) - 9-4- Y ) }
of v are also eigenfunctions for the same c, a. This can be seen a t once from
the symmetric stability equation and boundary conditions. It can be shown
further that either v0 or 9' is identically zero. To show this, we multiply
the stability equation (2.11) for pe by yo, and subtract ye times the same
equation for vo. This gives
vPv0 -v p v . =Q
where w # c. Therefore
@hpo - q&pS = constant = value at wail
= 0.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW O F INVISCID FZUID 19
either a - 0 or a + a , -.
20 P. G . DRAZIN AND L. N. HOWARD
q(x,t) = F(y)eiacX- c t ) ,
The trivial solution appears as the first term in a power series expansion
of Q for small a. Heisenberg [S] found formally two solutions of the Ray-
leigh stability equation :
+
q+(y;a*,c)= (w - c ) { q i d y , c ) aaq,l(y,c)+ . . . + a*q~,d.y,c)+ . . .} (i= 1,2)
where
qlo(y',c)= 1, qeo(y,c)=
In these formulae the lower limit of integration is arbitrary, but may con-
veniently be taken as y l . The zeroth approximation for small a2 gives the
eigenvalue relation
(2.27)
1
YI
(w - c)-%y = 0.
F-constant x e-=Y as y +w
changing the limit of the eigenfunction F, which, like c a y , does not admit a
power series expansion in aB uniformly for y large.
In 1962 Drazin and Howard [S] considered long-wave disturbances of
unbounded and semi-bounded flows. The basis of their work for unbounded
flows is as follows. The stability equation (2.23)has two solutions F, (y;a,c)
defined by their asymptotic properties
F, WeFay, DF, N Fa e F a Y as y+ f 00.
These solutions are defined by the stability equation for given a,c,w(y),
and are in general independent. However, for an eigenfunction which
vanishes at y = f w,
F E K+F+ K-F- (- W < Y < W)
for some complex constants K,, the solutions F, being linearly dependent
when c is an eigenvalue corresponding to given a. Therefore the Wronskian
a general smoothly-varying shear layer behaves like the vortex sheet with
basic velocity
Indeed the limit (2.29)of c gives the exact result for this vortex sheet, as is
given in Section II.4.d. In fact there is also instability (but of smaller
growth) when w ( - 00) = ~ ( o o ) ,i.e. when the flow is of jet by@, the next
approximation then giving
as a 4 0 , where
at y = y,,,.
In discussing the eigenvalue relation for a general profile, bounded or
unbounded, we mentioned that c(a) may be a many-valued function. The
variational principle (2.18)suggested that there might be many values
of a for each c = ws and many values of c = w,. For symmetric profiles
we mentioned sinuous and varicose disturbances, when c is at least double
valued. We shall meet many-valued c in several examples of the Section 11.4,
finding that each branch of c(a) is well behaved and corresponds to a distinct
mode of instability. By continuity in a one might expect that for each
neutral eigensolution with c = ws, and therefore for approximately each
point of inflection, there is one mode of instability. Drazin and Howard [8]
considered unbounded flows and associated heuristically the neutral eigen-
solutions having c = w, for each zero ys of D'Jw with the small-a eigen-
solutions having c = w(y,,,) for each zero y,,, of Dw. However the general
problem defies oversimplification, and the modes have not been satis-
fact orily classified.
In the parallel flow case, the flow domain is taken as the strip
- do < <
< x < do, y1 y y,, and the boundary conditions are usually
taken to correspond to rigid walls at y = y1 and y,. In the general case,
it seems appropriate to take the flow domain to be some region R in the
plane or in space, and to allow the possibility of flow across the boundary
B of R. In this case, more attention has to be given to the boundary con-
ditions. We shall assume that R is fixed and that we have given in R a basic
flow U, which is steady, incompressible, and inviscid. We shall also
not consider any body forces. The stability problem is then formulated as
follows: given some initial conditions which differ from U, by a small
amount, in terms of some appropriate measure (e.g., the L , norm of the
difference), we find the time dependent flow U determined by these initial
conditions and some suitable boundary conditions which are satisfied by
U,. The flow is stable if U continues to differ from U, by a small amount
in terms of the selected measure. By “suitable boundary conditions” we
mean such as assure that the initial-value problem for the flow equations
will have a solution, and a unique one. The mathematical questions of
existence and uniqueness of the initial-value problem for the inviscid flow
equations with various boundary conditions do not appear to have received
as much study as they deserve, but it is not appropriate here to embark
on such a discussion. We shall give only some brief heuristic remarks. The
most familiar case is that the boundary B consists entirely of a rigid wall,
so that Us n = 0 on it. This condition probably is, in itself, a “suitable
boundary condition” in the above sense. More generally, if the flow crosses
B , we should expect to prescribe U n as a boundary condition on B , sub-
ject only to the requirement that its integral over B should vanish (other-
wise the cohtinuity equation alone would have no solution). However,
this condition alone is in general not sufficient to insure uniqueness of the
solution to the initial-value problem-consider for instance the plane
flows in the annulus 1 r < < 2 (polar coordinates) whose radial velocity
component is l / r and whose azimuthal component is ( l / r ) f ( r -
2 2t) where
f is a function which is zero for values of its argument 2 1, but is otherwise
arbitrary (it may be as smooth as desired). It is easily checked that these
velocity fields do satisfy the flow equations, they all have the same normal
component on the boundaries of the annulus, and are identical at t = 0.
In the example just given (two-dimensional and axisymmetric flow in
an annulus), uniqueness can be insured by prescribing, in addition to the
normal velocity component on the complete boundary, the tangential com-
ponent on the part (r = 1) of the boundary through which the fluid enters
the flow region, for all t > O . I t is clear that this additional information
is just sufficient to determine the function f . Note also that this example
shows that one may not in general prescribe the tangential component
where the fluid leaves the region (r = 2 ) .
24 P. G. DRAZIN AND L. N. HOWARD
We shall now show that these boundary conditions, namely the pres-
cription of the normal velocity component over the complete boundary
(in a manner consistent with the continuity equation) and of the tangential
velocity components over that part of the boundary where the flow is inward,
together with one additional assumption, are sufficient to insure the unique-
ness of the solution to the initial-value problem for general three-dimen-
sional incompressible inviscid flow.
Let U, be such a flow, in a region R with boundary B, and let
+
U = U, u be another. We assume that U and U, satisfy the same bound-
ary conditions, so that u * n = 0 on B and u = 0 on that part B, of B on
.
which U, n < 0. Write B, for the rest of B , on which U, n 2 0. Consider
+
now the deformation tensor D of U, , with components Dij = U&,j U , , .
<
Since Djj has zero trace, at least one of its eigenvalues is 0; let C(t) be
the supremum over R of the absolute value of the most negative eigenvalue
of D-we call C(t) the maximum shear of U,, and we assume that C(t) is
+
finite, initially and thereafter. Writing for the difference of the pressure
fields of the flows U and U, , divided by density, one obtains the following
equation for the “perturbation” u by subtracting the momentum equation
for U, from that for IT:
(2.31) Ut + U-Vu + U* m, + v+ = 0.
Note that the “perturbation” u is not necessarily small. Multiplying (2.31)
by u and using V oU = 0 we get
Integrating this over R , applying the divergence theorem to the last term,
and using the boundary conditions, we get
(2.34)
and thus
t
(2.36)
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 26
initial-value problem for equations (2.36) and (2.37), are the prescription
of Y (or U n) on B , and of R on B,. In fact, with these boundary conditions
and the assumptions that (U,(and R, remain finite, one can prove unique-
ness by a method similar to that used above for the case of the tangential
velocity components being given on B,.
We now consider the stability problem for plane flow, using this vorticity
boundary condition, to establish a result which may be regarded as giving
a generalization to non-parallel plane flow of Rayleigh’s inflection point
theorem. Let the basic steady flow be U, = vYo x k. Its vorticity 51,
= - AY, is constant along streamlines, and we shall write Ro = /(Yo),
though in some cases such a representation is not literally possible unless
f is regarded as multiple valued-different streamlines might carry the
same value of Yobut different values of a,. We are going to prove that
if f’(Y,)< 0 throughout the field of flow, then the flow is stable to two-
dimensional disturbances. As with Rayleigh’s theorem, this will be a suffi-
cient, but not necessary, condition for stability. The perturbation momen-
tum equation is:
(2.38) U$ + Qk x u + wk x U, + Vh = 0,
where h is the perturbation H - H , of the total head H = ilU21 + P.
From (2.38) we deduce:
(2.39) -”(’
atZ Iul2) + w u * (k x U,) + (uh) = 0.
(2.40) W# + u. v w + U * VR, = 0,
and from this we get:
(2.41)
Now
WU. VQ, = f’(Yo)wuoW, = f’w(u x k) (W,
x k) = /’COU* (k x U,).
If we now assume f’< 0 we can rewrite (2.41), using the fact that f’ is in-
dependent of t and constant along the basic streamlines, as:
Subtracting (2.42) from (2.39) and integrating over the flow region R , using
the divergence theorem and the boundary conditions U. n = 0 on B and
w = 0 on B, we get:
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 27
(2.44) <
0 2 / f ' ] d A 0.
Since /' < 0, this integral is positive definite, and it follows that the energy
of the perturbation, though it may possibly increase somewhat over its
initial value at the expense of the term J - w2/2f'dA,must remain bounded;
R
thus the flow is stable. It is interesting to note that the restriction to the
linear stability theory is not necessary for the (rather special) class of basic
flows which have f' constant and negative.
The relation of this result to Rayleigh's theorem is easy to see. For a
parallel flow with velocity profile .I&) we have
If there are no inflection points, w" is of one sign, say w f f >O. By adding
a suitable uniform translation if necessary which obviously does not affect
the stability properties of a parallel flow, we can assume that w >0 through-
out, and so f'c 0 and the flow is stable. (If w"< 0 we can add a suitable
uniform translation so that w < 0; but if w" changes sign this is not possible.)
If there is just one inflection point we can assume that w = 0 there, and
our result thus implies stability if w"/w >O; we thus also obtain Fjertoft's
extension of Rayleigh's theorem.
The above argument is not applicable if f' E 0, i.e. for constant vor-
ticity (in particular irrotational) flows. However such flows, like their
parallel prototype the plane Couette flow, are always stable, a t least with
the boundary conditions we have assumed here. For when VQo = 0, the
vorticity equation (2.40) shows that the perturbation vorticity w is constant
following particles. Since no new perturbation vorticity is brought in by
entering fluid particles, w cannot grow, and since the perturbation stream
function is determined from the Poisson equation A$ +
w = 0 with t,4 = 0
on B , 4 cannot grow either.
28 P. G. DRAZIN AND L. N. HOWARD
(;+.;)(%+!$)=a.
Therefore
-m
.l'
) da cos a x
@ ( x , y , ~=
m
2 bn(a)sin +nn(y+ 1 ) .
n-1
#(x,y,t) = da
-w i 2 #,,(aB+ +fi2na)cosech 2a
ml:
Evaluation of this solution for large t shows that @ = O(t-1) and therefore
that plane Couette flow is stable.
A more systematic approach to the initial-value problem, more suitable
for application to basic flows other than plane Couette, comes from use
of the Laplace transform with respect to time. This approach has been
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 20
developed for stability problems by Miles [27], Carrier and Chang [28],
Case [29. 301 and Dikii [31. 321. It has been reconciled with the method
of normal modes, both for inviscid and slightly viscous fluid, by Case [33]
and Lin [34].
First let us take Case's [29] solution of plane Couette flow, for comparison
with the other solution above. Let the Fourier transform with respect to
x and the Laplace transform with respect to t give
p(y,a,t) = 1 e- %,h'(x.y,t)dx,
i
Y ( y ; a , p )G e- ~'!P(y,a,t)dt.
0
Given the conditions that +' = 0 on the boundaries for all x,t and therefore
that
!P,P=O ( y = *I),
sinhall + v,,)sinhall - v )
a sinh 2a \d" .d
(2.50)
j.
PU(y;a,$)= G(y,yo)H(Yo,a)/($ iaw0)dro;
YI
+
where
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 31
(2.51)
such that
As in the case of plane Couette flow, we can now invert the transforms, etc.,
to find $'(x,y,t) as an explicit multiple integral of #'(x,y,O).
The crucial integral is the inversion of the Laplace transform,
(2.52)
We seek to evaluate this integral for large t to see whether the flow is stable.
When 1 is large, (ept(is large along C where Re p > 0. However, as we shall
see, cancellation of the large components of the integrand means that the
integral may not be large. The cancellation makes it advantageous to move
the contour close to the imaginary axis and thereby reduce the magnitude
of lePtI. Following the usual method of deforming the contour C, we must
study the singularities of Y, which is given by equation (2.50). Since Yl,!P,
satisfy equation (2.&51),their only singularities can be at the singularities
of this equation, namely where w ( y ) = ip/u = c. These singularities lie
only on the imaginary axis of the p-plane. Thus it is convenient to take
a new contour on the imaginary axis with indentations to the right of each
singularity or just to the right of the imaginary axis, say from E - ioo to
E + ioo for small E . To find P in terms of the integral along the new contour
we must know the other singularities of the integrand in the half-plane
Re $ > 0 ; these coming from zeros of the denominator of the Green's func-
tion G(y,y,), i.e. from the solutions of
W(Y1.Y2)
=0 (ReP>O).
This can be seen to be just the eigenvalue relation for the discrete spectrum
of normal modes with ci >O. [Of course, equation (2.51) is essentially the
Rayleigh stability equation with p = - iuc.] Supposing these eigenvalues
are already known, we may evaluate the residues a t these singularities of
the integrand in equation (2.52) and find
1
E+iW
exponentials
!P(r;u,t) = -
2ni Yb,a,+)ePtdt +
e-iw
The fastest growing component will be dominant for large t. It will come
from the discrete spectrum if there is one, being the term growing like
8' for the largest value of uc{. Case [29, p. 1481 and Dikii [33, p. 11801
32 P. G. DRAZIN AND L. N. HOWARD
have indicated that the integral above over the continuous spectrum decays
like l / t , so the discrete spectrum alone is associated with instability. Thus,
in seeking a criterion for instability, we h a y use the method of normal
modes and ignore the continuous spectrum.
4. Stability Characteristics of Variozcs Basic Flows
As will be seen, most of the basic flows with known stability charac-
zw
teristics are piecewise linear. However, the characteristics of many smoothly-
Y Y
FIG. 4. (a) Plane Couette flow: w = y. (a) Plane Poiseuille flow: w = 1 - ya.
(c) Sinusoidal flow: w = sin y. (d) Vortex sheet: w = VS(cd,2 y > 0). w = V,(O>
y2 - dJ. (e) Rectangular jet: w = I(lyl c I), w = O(lyl> 1). ( f ) Thin jet w = (d(y))l/*.
(9) (i) Symmetric jet in channel. (g) (ii) Antisymmetric shear layer in channel.
HYDRODYNAMIC STABILITY OF P A W E L FLOW OF INVISCID FLUID 33
(k)
(01 (PI
varying basic flows are partially known, and the advent of electronic com-
puters is allowing many to be found completely in numerical terms. Diagrams
of the basic flows accompany the examples below.
( D 2- a2)q= S(y - c ) ,
where 6 is the Dirac delta-function. This representation in terms of a general-
ized function is admissible because the typical component of wavenumber
a is really only one component in an integral. It gives a solution of the
continuous c-spectrum with eigenf unction
v= II
sinh a(c - 1)sinh a ( y
sinh a(c
a sinh 2u
+ 1)
+ 1)sinh cr(y - 1)
a sinh 2a
(- 1<Y<C)
(c <Y < 1)
for each value of a and for each value of c in the interval (- 1,l). The set
of eigenfunctions is complete so that an arbitrary initial disturbance of the
velocity field can be represented as a sum orintegral of them [cf.36,pp. 11-12].
In finding the s-solution we ignore the factor sin y (and thereby discard the
stable eigensolution corresponding to c = 0 in the continuous spectrum)
to get
vs = sin M Y - Yl)/(Y2- Y1))>
as = (1 - nznz/(yz- y1)2}1/2
for each positive integer 12 < (y, - yl)/n. I t follows that the flow is unstable
if (r, - yl) > n , but stable otherwise although the point of inflection lies
at y = 0 in the field of flow. This counter-example to the sufficiency of
Rayleigh's necessary condition for instability is due to Tollmien [14; see
also 7, pp. 219-2201,
When
Therefore
c = { V, coth ad, + V, coth adl
+ ilVa - V1l [coth ad1 coth ada]l/e)/{coth ad, + coth ad,}.
This gives complex c and instability for each a.
When d,,dg = m,V2 = 1 = - V,, this gives c = i. This is an example
of an antisymmetric flow with c, = 0. It also gives the limiting eigenvalue
(2.29) as a -P 0 of any smoothly-varying flow with w ( 00)= 1 = - w ( - cm).
The elevation of the material surface with mean level y = 0 is q =
+
If V,+ V , = V and dl,da = 00, one finds c = - V and q = (A Bt)e”(x- vl)
for arbitrary constants A , B . This has been described as the instability of
a flapping flag.
and
[
(ii) For an antisymmetric shear layer with
V +W Y - ib’) (ib‘ < y <b + i q
w= 2Vy/b’ (IY I< tb‘)
-V + nv(y + ib’) (+b’ 3 y 2 - b - gb‘)
Rayleigh gives [3, p. 3881
+ a sinh a(b + b‘)]*
cB= P {[(A - 2/b‘) sinh d sinh ab’
- aBsinhaab}/{assinh ab’sinh a(2b + b’)}.
Thus there is instability for some a only when - l i b < A < - l / b + 2ib’.
This result in some sense exemplifies Fjmtoft’s result for smoothly-varying
profiles represented in Figure 2.
38 P. G. DRAZIN AND L. N. HOWARD
0 (IYl>1)
a)@ -4 (1 > IYI > 4
(a > IY I)
the eigenvalue relation is [3, p. 3971
4(1 - - 2(1 - a)ac{2(1 - a)a ‘f e-*(l - e-2(1-a)a)}
+ {- 1 +2(1 - a)a + r2(1-a)a} e-*{1 - [I + 2(1 - a)a]e-2(1-a)a}
= 0,
where the upper sign is for the sinuous and the lower for the varicose mode.
When a = 1 we get an example of the rectangular jet (Section II.4.e).
When a = 0 we get a triangular jet, with
2aW +a(1- 2a - e-%)c + (a(1 + e-2”) -1 +e-a) =O
When
0 (IYl>1)
w = I u (-l<y<O)
1 (O< Y < 1)
the eigenvalue relation is [S, p. 2691
+ (1 - c ) tanha}{c*
(1 - C ) ~ { C % ~ tanha + (U- c)~}
+ (U- tanh a + (1 - c)%}{ca + (U- c)%tanh a} = 0.
C)~(C%
or
U{1 + i(ta)'/*+ ...}
for small a.
(k) S k a r Layer
When
= { rllrl (Irl>4
Yla (Irl<a)
the eigenvalue relation is [3, p. 3931
-2 4 *
c* = (4~*a*)-~{(1 - e-&}.
Only two special cases are of interest : a = 0 and a = 1. The former gives
Kelvin-Helmholtz instability of the vortex sheet (Section II.4.d). When
a = 1, c is pure imaginary for a< a, I 0.64 and real for a 2 a,. The
logarithmic growth rate acj is greatest when a = 0.4.
<
There are two modes unstable for each a. For a 4 tanh-1 (6l/*-- 2) both
unstable roots c are pure imaginary; for a >4 tanh-l (S1le - 1) both are
complex. Thus an antisymmetric profile does not neceSSarily give c, = 0.
for the varicose mode. For small a the sinuous mode has [8, p. 2811
c =a + i{#a - a* - i d l o g (24a-1) - &azni}l/*+ . , ,
and the varicose [8, p. 2791
c =1 + e(2/a)d(&a)21a(l+ gC(Wa)"(#n2a2)1/8 + O(a)}.
Numerical values of q,c, for some a have been given by Lessen and Fox
[38; cf. 81. We give results of a somewhat more complete recent calculation
in Table 1.
Stability Characteristics of other jet or wake profiles have been computed
by Hollingdale [39], Haurwitz and Panofsky [40], and Sato [41].
TABLE1. INSTABILITY
CHARACTERISTICS
OF THE B I C KL ~Y
JET
a cr ci Cr ci
-
TABLE2. INSTABILITY
CHARACTERISTICS
FOR UJ = 9 1 3 sech* y tanh y
a ci a 6, ci
0 0 0 1.o 0
0.1 0.270 0.05 0.944 0.071
0.2 0.362 0.1 0.903 0.103
0.3 0.426 0.2 0.831 0.133
0.6 0.495 0.3 0.775 0.139
0.7 0.612 0.4 0.732 0.130
1.o 0.472 0.6 0.700 0.116
1.5 0.317 0.6 0.678 0.098
1.8 0.199 0.7 0.660 0.077
2.0 0.114 0.8 0.648 0.062
6 0 0.9 0.641 0.044
42 P. G. DRAZIN AND L. N. HOWARD
The special case m = 0 of the example in Section I1.4.n gives the shear
layer
which has better known stability characteristics [S, p. 2811. The s-eigen-
solution is [43, 441
~d =1 - 1.786a + 1.52W + .. . .
Numerical calculations for this case have recently been given by Michalke
[&I.
Stability characteristics of some other smoothly-varying profiles of
shear-layer type have been computed by Hollingdale [39], Carrier [cf. 461
and Lessen and Fox [38]. Their results are similar to those for the hyperbolic
tangent shear layer.
When
Lin pointed out [47, p. 901 that the Rayleigh stability equation can be
transformed into the hypergeometric equation, a result used by Chiarulli
and Freeman 1471.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 43
1. Introrlzcction
The theory in this section consists chiefly of developments by Scandina-
vian meteorologistsin the 1920’s and 1930’s of Kelvin’s paper [2] of 1871. Their
work is collected in a book by Bjernes, Bjerknes, Solberg, and Bergeron [48]
of 1933. It may be more accessible to the reader in Chapter X of the book by
Godske. Bergeron, Bjerknes, and Bundgaard [49]. which also has a biblio-
graphy of relevant meteorological papers published by 1960. These me-
teorologists have considered the equations, boundary conditions, and eigen-
solutions for piecewise-constant velocity profiles under the influences of
combinations of density variation, compressibility, and rotation. Haunvitz
[SO] has considered the equations and boundary conditions under the same
combinations of force fields for smoothly-varying profiles also.
The very generality of these combinations and their meteorological
context has obscured some of the fluid dynamics and enabled other authors
to duplicate their work in ignorance of it. So we shall consider the force
fields separately in order to simplify the understanding and to compare
the effects of different force fields. Of course any combination of these
force fields can be considered, and. indeed, there are dozens of papers dealing
with various combinations. However these combinations need no special
techniques in their treatment, so we shall not describe then. If the reader
is interested in a physical problem of stability of parallel flow affected by some
combination, he may readily adapt the methods used for the fields separately.
Again, we have ignored the effect of surface tension on an interfacial
boundary condition. This may be physically important, but it is a simple
matter to apply the method of Kelvin [2] to cater for surface tension.
Treating one external force field at a time, and finding its effects on the
inertial instability discussed in the previous section. we first give the stability
equation and boundary conditions for two-dimensional wave-disturbances,
and comment on the validity or invalidity of Squire’s theorem. Then we find
the eigensolutions for two unbounded basic velocity profiles, that of static
equilibrium (w = 0), and that of a vortex sheet (w = y/lyl). It is found
that for static equilibrium a neutrally stable wave motion may occur for
most of the external force fields. For example, sound waves may occur
when the external “force” is due to compressibility of the fluid. Such a
wave motion is important both for its own sake and for its representation
of some stability characteristics of profiles with shear. The example of a
vortex sheet is also important, because a vortex sheet is the simplest flow
with shear, and because its stability characteristics represent those of any
smoothly-varying shear layer for long-wave disturbances. Finally we sum-
marize briefly the literature on stability problems for each force field.
44 P. G. DRAZIN AND L. N. HOWARD
Following Kelvin [2], many authors have studied the stability of parallel
horizontal basic flows of incompressible fluid with piecewise-constant
velocity and density distributions (i.e. with velocity and density uniform
in layers but varying from layer to layer) under the action of gravity.
Rayleigh [61] was the first to consider the stability of fluid of smoothly-
varying density distribution p,,,(y*) at rest, y* being the height. For the
particular density distribution p* = po exp (- By) (- bo < y* < bo, con-
stants Po, fi >O), he found neutrally-stable internal gravity waves of phase
velocity
(34 c* = a, = (gp/a,*)'lS,
on neglect of a&, i.e. on neglect of the variation of inertia due to the
variation of density but not of the buoyancy.
Taylor [62], Goldstein [63],and Haurwitz [60] consideredtwo-dimensional
disturbances of parallel horizontal flow of incompressible fluid under gravity
with smoothly-varying profiles of velocity and density. Their angysis
leads to the stability equation with dimensionless form,
(w - c)(D*- a*)p,- (0%)~ +
J y / ( w - c) - K{(w - c)Dy - ( h ) p , } = 0;
(3.2)
where the local Richardson number of the basic flow is
(3.3) =
Jcy) - ( g L B d p * / d Y * ) / W ,
8
In fact Squire’s theorem may be extended to this case [54], giving for
each three-dimensional wave making angle 0 with the basic flow a two-
dimensional wave of the same growth rate but effective Richardson number
Jcos*O and Froude number Fcos*0. Usually density variation with
dp,/dy, < 0 acts as a stabilizing agent, so we seek criteria of stability with
minima of J or F. Thus two-dimensional disturbances are the most unstable
usually, and anyway the stability characteristics of any three-dimensional
disturbance follow at once from knowledge of the characteristics of all
the two-dimensional disturbances,
The outstanding problem of this case is Kelvin-Helmholtz instability of
the basic vortex sheet with
(3.7)
Thus the flow is always unstable if there is heavy fluid above lighter
< A*). If there is heavy fluid below lighter, the flow is stable to those
waves with
(3.10)
for neutrally-stable internal waves at the interface of the fluids. Taylor [66]
recognized that the same analysis gives
The square-roots must be chosen with non-negative real parts in order that
the eigenfunction (3.13) is bounded at y = f 00. If a square-root is pure
imaginary. its sign must be chosen so that there is outward radiation of
energy at infinity; however, this occurs only for real c, in which case there
is stability anyway. It now follows from squaring up equation (3.14) that
(3.16) c=O and 8 2 1 ; c*/a*=l and a t = 00;
or
ca = - 1 + &as.
The second mode represents Rayleigh's internal gravity wave with I'= 0
and c,, = f a,. The first and second modes are isolated from one another,
and from the third mode, which is the only one that can give instability.
It can be seen that there is stability to all waves only when
(3.16) as 2 2, i.e. aa < +gp/Va.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW O F INVISCID FLUID 47
where p*P, is the pressure plus the centrifugal potential. Thus the vorticity
with respect to a non-rotating (inertial) frame is 2n plus the relative vorticity
O, = V x u*. Using the modified vorticity equation and the usual methods
of normal modes, Johnson [66] has found the stability equation of a three-
dimensional disturbance,
(W cos e - c ) % ( D-~ a9)cp - cos qrar
cos 8 - c ) ( D * w ) ~ ~
(3.24)
+ Ro-1 cos x(R,,-l + sin 8Dw)cp = 0 ;
cos x
where pl is defined by v' = iav(y)exp {i(ax + yz - act)}, x is the angle
between and the wave-number vector = (a,O,y),and the Rossby number is
a, = 2Q cos x/a*.
For the vortex sheet w = y/lyl (- ao< y < ao) the eigenfunction is
(3.28) 9, =
1+
(c - cos 8) exp (- ay(1 - a2/(c- cos 8)z}1/2)
+
(c cos 8) exp (ay(1 - a’J/(c cos 8)2]1/2)
and the eigenvalue relation is
(c - cos e)yi - ay(c - COS e ) y
(3.29)
+ (C + cos e)z(l - u ~ / (+c cos 8)z}1’z= 2a sin 8.
, I
When 8 = 0 this eigenvalue relation has the same form as that of (3.13),
(3.14) for a heterogeneous fluid in a non-rotating system. On squaring up
(3.29) etc. for general 8, we get the cubic equation in cz,
(3.30)
o =~ ( C Z ) 4 cosz ec6 + (8 cos4 e - a2(i + 3 cOSae)}c4
+ ( ~ ~ ~ S ~ ~ - ~ U ~ C O S ~ O ( ~+ - UC ~O }S C~ O~ )+ U ~ ~ ~ ~ ~ ~ ( ~ ~ - C O S
The first root is inadmissible. The latter pair give only a complex conjugate
pair of roots c when cos2 8 > 1/9, implying instability. When cos2 8 \< 1/9
there is stability provided a sin 8 2 0 but not otherwise.
The Scandinavian meteorologists [cf. 48,491 have considered the stability
characteristics of various piecewise constant velocity profiles and Johnson
[66] has treated the shear layer w = tanhy.
where
(3.36) as E pHo2/4nVapp,,
(3.371
. .
[{l- u4/(w - c)2} { ( w - c)Dp - (Dw)q}]= 0.
(2) When RM = 0 and the basic magnetic field is parallel to the flow,
the stability equation can be shown to be [70]
where
(3.41) N pHo2LL/4np,/V.
The boundary conditions are the same as with no magnetic field. There
is no progressive wave possible when w = 0, and the vortex sheet is unstable
for all values of H,, however large, although the magnetic field reduces
the instability [cf. 721.
1. Dimensional Analysis
The methods we shall use for.this problem can be readily applied to the
other stability problems of Section 111, which have a similar form. The
stability equation above shows that gravity occurs only in the product
- gD,p+/p+, for J ( y ) = - g(D,p,)La/p,Va. Therefore the eigenvalue
problem (3.12), (2.12) gives eigenvalues of the form
for the class of similar profiles w(y),p(y), where Jo* is the value of J + 3
-gD,p,/p, at any specified point y . Now dimensional analysis implies
that
(4.2) c = c(a,Jo)
where Jo = Jo+La/V~ acts as a characteristic value of J ( y ) . To solve a prob-
lem we find this relation explicitly, and, in particular, find the values Jo(a)
for which ci(a,Jo)= 0 but for which ci(a,Jo)> O nearby. These values
Jo(a)define the curve of neutral stability (or ~eutralcurve or stability boundary)
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 55
-
there is instability or not for any given wave and flow defined by the point
(a,Jo)*
Now consider the limit as Jo 0 for fixed a # 0. Then (some of) the
eigenvalues
(4.3) c - c(a,o),
which limit we suppose to exist and to equal the eigenvalues c(a) of the
Rayleigh stability equation (2.11), which is equation (4.1) with Jo 0.
It should be borne in mind that internal gravity waves exist for all Jo >O,
-.
ledge of c(a) from the theory of section 11 now tells us the behavior of
some of the branches of the eigenvalues c(a,Jo) as Jo 0, i.e. as the buoyancy
becomes small, as the velocity scale becomes large, or as the length scale
becomes small.
-
To consider mechanism (c) alone let us take the limit as w* + 0 for an
unbounded flow, i.e. as V 0 for fixed w ( y ) , p ( y ) , L, Jo*, u,,. In this limit
we suppose that c* tends to a function which is independent of w*, and
therefore of both L and V , because as w* vanishes its length scale and shape
cannot be relevant physically. Therefore c* is some function of a*, Joe
which has dimensions of velocity. This implies that
c*-kJi$/a*, as V 40,
where k is a (many-valued) dimensionless constant independent of wCy)
but dependent on pCy). But when V = 0 and the flow is unbounded it is
well known that there are internal gravity waves whose speeds do depend
on p b ) . These speeds will give k. For example, when
P* = Po* exp (- BY*),
it can be shown [Sl] that
cs = (gS)"%* :
therefore k = 1 if we choose Jo* = gfi. (Of course the arbitrary multi-
plicative constant in Joe affects k, because it is only the product k Ji:/+"
that
is determined physically.)
To determine mechanisms (a), (c) together let us again suppose that
the flow is unbounded. We can now let the length scale L of the velocity
profile tend to zero without altering the infinite domain of flow. Thus we
let L - 0 while a*, V , w ( y ) , g, p(y) are fixed. Then a = a,L + O and
Jo = - gL(Dp/p)o/Va+ 0, although Jo/a is fixed. Thus if we write
c = c(a,Jo/a)
66 P. G . DRAZIN AND L. N. HOWARD
(Y* > 0)
P*(Y*) + P*o (Y* = 0)
Q
'::[ (Y*<O)
-
as L + O if these limits exist.
Now let us review what happens in the limit as L - 0 for a profile of
-. -
from there being more than one mode of instability for a smoothly-varying
shear layer, yet only one for a vortex sheet; again the limits c, 0, a 0
may not be uniform.
2. Physical Argtcmnts
The mechanism of instability of a vortex sheet w* = Vy*/ly*I in a
compressible fluid at uniform temperature will now be described, essentially
in the way attributed to Ackeret [cf. 74, p. 2401. Consider a small irrota-
tional two-dimensional disturbance of the velocitv field in which the interface
between the two streams of speeds V , - V is distorted. Thus the interface
has small bends. If the streams are subsonic ( V < a*), then by continuity
the speed on the convex side of a bend has a small increase over its basic
value, and the flow on the concave side a small decrease. Now Bernoulli's
theorem for irrotational unsteady flow of barotropic inviscid fluid plausibly
suggests that the pressure decreases on the convex side and increases
on the concave side of the bend. This pressure difference induced across the
bend increases the curvature of the bend and thus causes instability of
the interface. By the theory of the Laval nozzle [cf. 74, 8 3.61 the speed
58 P. G. DRAZIN AND L. N. HOWARD
as y, -.
above its basic level is 11, =F,(y*) exp(ia,(x, - cat,)} where F, -FO*e-OL*IY*I
f 00. The density is j,-=below the jet. Close to the jet, within
a distance of order of magnitude L from the jet, i.e. much closer than a
wavelength, the long waves seem locally like a vertical translation of the
jet. Thus the jet oscillates like a string with form
110, = FO, exp {ia*(z* - c*t*)h
for q , is approximately constant on any vertical line near the jet.
In this motion the vertical mass-acceleration of the fluid on both sides
of the jet is in balance with buoyancy and the centrifugal force due to the
(small) curvature of the jet. The buoyancy comes from lifting fluid below
the jet a height 7, into space previously occupied by the lighter fluid above
the jet and vice versa. In this way the buoyancy gives rise to a pressure
disturbance across the jet
= g{p*-a, - p*m}qo*s
to first order for small q,.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 69
s
m
= a*2 P*W*%*dY*.
-m
Now over the effective width of the jet q+ = qo*, because a* is small; in
the distant regions where q* -+ 0 exponentially w* is small anyway. Therefore
this pressure difference
I
m
= a*%* P*w*2dY*
-m
Now /3, changes rapidly from its value at the origin to its values
at infinity, whereas q,, changes slowly like Fo,e-"*lY*l +"*(+* - c*t*) . There-
fore, for small a*, this expression for the mass-acceleration
m 0
= - a*cla{jkm + P*-m)qo*.
Finally the balance of pressure and mass-acceleration per unit area gives
1.e.
60 P. G. DRAZIN AND L. N. HC'YARD
as a* + O for fixed
g{P*-m - ij*m}/V8a*@*-c P*mh
g{A-m - P * m ) 2 a*
.jP*w*'~Y**
-m
V. INSTABILITY
OF AN FLUID
INCOMPRESSIBLE OF 1 ARIABLE DENSITY
. .
(6.3)
+ W-zn(n(I - n)(Bw)2- J ) ) H = 0.
(5.4)
Y, T We(l-*l(lDHI2 + a2lH12) + nWl-"(D%m)IHJ2
+ W-%{n(l -~ ) ( D w-)J))HI2dy
~ = 0.
This result of Howard [23] can lead to various properties according to the
value of n chosen.
When n = 1 we have
(5.6) +
[ ( D q ~ l ~oc21pl12 + W-'(D2w)(cp(' - W-2J1~\adr= 0.
YI
Therefore
(5.7) Dew = 2(w - c , ) J ( y ) / ( (-~ c,l2 + ci*}
somewhere in the field of flow. If D2w # 0 in the field of flow we further
have
(6.9) +
r ( w - c)*((DFI* a*lF(*)- JlFl*dy = 0.
Yl
This leads to the proof of the semicircle theorem, as in Section 11.2. The
extra term in the present case only strengthens the inequalities used provided
J 0 everywhere. Thus, when ci > O and J ( y ) 2 0 in the field of flow,
(6.10) {c, - a ( ~ & + Wna.)}~ + cis < { a ~ ~ m -a s wultn)}**
Howard proved this result [23] for a heterogeneous fluid originally. Even
when J < 0 somewhere, it follows that wulti,< c,< wmx. When J ( y ) 0 <
everywhere, equation (6.9) shows that no non-singular neutral mode can
exist, i.e. that either ci # 0 or c lies within the range of w and F is therefore
singular. However, when J ( y ) > 0 somewhere, it is possible that non-singular
neutral modes exist with c outside the range [w-,w-] of w ( y ) ; these
isolated neutral modes in fact occur as internal gravity waves.
When n = 4 we have
j.,
Y1
- c){lDHl* + a*lHl*}+ i(D2w)lHl*+ W-l(t(Dw)*- filHl*dy= 0.
(6.11)
(6.13)
Y,
P P
0 > - lDHlady = (aa
9,
+ {J - f(Dw)*}/lW(*)IHl*dy.
Therefore J C y ) < $(Dw)* somewhere in the field of flow. This gives Miles'
[77] sufficient condition of stability that J - t(Dw)*should be everywhere
non-negative. Further from inequality (6.13) we have
I j
Y,
(da*/dc)b= lim
e+a,
j. (2JW- 1--zn- w-ymY+
(0%)
P ~ ( 1 * -w a y ,
Yr YI
(6.16)
64 P. G. DRAZIN AND L. N. HOWARD
(6.18) lim
a+a#
i
YI
qp
w2(1-
i
d y = WF(1- n)Hsadr,
Yl
(5.20)
(6.21)
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 66
-
of unbounded flow for long waves. Their method is a natural generalization
-
of that for a homogeneous fluid with J = 0. The eigenfunction must be
such that p, constant x exp (Fay) as y f m, when DP 0 smoothly
-+
(6.22)
0= a(W-,
2
+ Wma) - 2J0+ 1
--m
{a(W2- Wm2)
This result in the limit as a --+ 0 agrees with that of (3.9)for Kelvin-Helmholtz
+
instability of a vortex sheet when (p-m - Pm)/(P-m Pm) << 1 , i.e. when
K ( y ) 0. For profiles of jet type with wim = 0, relation (6.22) gives
1
m
ca = Jo/a - aa +
( W z - c ~ ) ~ / 2Jo(Wa
W ~ - ca)/aWa
--m
(6.26)
+ Joa(l - 12)/a2W2dy+ . . . .
When w E 0 everywhere this relation in turn gives the speeds of the internal
gravity waves. For w(r) not identically zero, these internal gravity waves
are modified and become isolated stable modes. There are also a t least
two unstable modes for sufficiently small a,Jo. When Jo is of order aa,
equation (6.26) gives the sinuous mode with
(6.26) as a -0,
-m
there being stability to long waves when Jo is greater than this value.
as y -
It is physically realistic to have a model for which p ( r ) tends to constants
f 00. However, in many circumstances it is practical to suppose
p tends to zero or infinity. For example, one might suppose that P = po
exp (- By) for constant p > O to represent the atmospheric density in a
problem of instability that occurs effectively in a finite range of y ; in this
case J ( y ) = g p L Z / V ais independent of y . We have shown (3.1) that then
internal gravity waves have speed a = J1I2/a. I t can also be seen from the
stability equation (6.1) and boundary conditions (6.2) at infinity that
y - +
where L , = (1 - U ~ / W ; , ) ' / ~
provided
, that w(r) --+w,, smoothly as
f m. Now a direct generalization of the method of Drazin and Howard
[79] described above gives for small a and fixed a the eigenvalue relation
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 07
0 = L+wma +L-w-, + a
2
J J
- w 0
- L+L-Wm'
-m
5
0
1 - Wt-mjw'dy} + ... .
For a shear layer with I-, = - w , = - I, this relation in the first ap-
proximation gives the same stability characteristics (3.15) as a vortex sheet.
Higher approximations are obscured by the difficulty discussed at the end
of Section IV.2 for a compressible fluid. For a jet with w + = ~ 0, relation
(6.29) gives
(5.30) 0 = 2Lcs +a
5
--oo
(Wa - ca)(1 - L*c2/Wa)dy + ... ,
+
where L = (1 - U ~ / C ~ ) ) ' /In
~ . the first approximation for small a this
gives an unstable mode with
(5.31) c--li(a/u)
2
1
W
--m
wady as a -0.
In this special case with constant J and unbounded flow, we can find
a sufficient condition for stability. We have required that the real parts
of L , be non-negative in formula (6.29) in order that (p does not exponentially
increase at inwity. In fact the solution of tb initial-value problem must
die down as y + f 00, so isolated modes may have non-zero or even un-
bounded eigenfunctions at infinity, but any dense set of waves must have
eigenfunctions which tend to zero there. Therefore, when c is real, only
isolated waves may be not exponentially damped as y + f 00. It follows
that in the limit as ci --* 0 , L+ and L- are real and non-negative. Therefore
68 P. G. DRAZIN AND L. N. HOWARD
Again (5.36 i-iii) are eigensolutions, but it seems [83] that none is a
stability boundary.
where
-1 (Y< -1)
(6.42) W = l Y (-l<y<l)
1 (1< Y )
was first considered by Taylor [62, $31 and Goldstein [63,$53, 61. They
took essentially
i
8-m (ye-1)
(6.43) p= A (-l<y<l)
Pm ( 1 <Y )
and found the eigenvalue relation. When p,, = i@-= + p-) and
+
(p-- - pm)/@-aD pm)<< 1, there is instability if and only if
(6.44) 2a/(l + e-%) - 1< J0< 2 a / ( l - e-*) -1
where Jo -ggL@-, - &)/V*(p-, +
pm). When Jo = 0 this reduces
to Rayleigh’s result [cf. Section II.4.k] that there is instability if 0 < a < a, =k
0.64. For general values of Jo,a the stability boundaries etc. are shown in
Figure 6, after Goldstein. It can be seen that any wave unstable when
J o = 0 becomes stable when Jo is sufficiently large. However, other waves
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 71
r
are made unstable as Jo increases, some n a o w band of waves being unstable
however large Jo is.
Goldstein [63,5 61 further considered the shear layer with the continuous
density distribution
(Y< - 1)
(5.6) p= p-me-b(y+l) (- l < y < 1)
p- w (1< Y)l
where again j? is negligibly small except when multiplied by gravity.
Goldstein’s solution, involving Bessel functions, is complicated. In brief,
I I I
‘ I
it gives stability to all waves if and only if Jo g,8,La/Vn >, t. This condi-
tion is in agreement with Miles’ sufficient condition J/(Da)P>, for stability.
When a << 1, it can be shown from Goldstein’s work that is there is stability
when
J~ 2 - 2 - -as
4 - -a4
I6 -
(6.46)
3 9 46
...,
in agreement with relation (6.22). Recently Miles and Howard [a] have
clarified an obscure point in Goldstein’s paper and given some numerical
72 P. G. DRAZIN AND L. N. HOWARD
results for this example. The principal stability characteristics are shown
in Figure 6.
Holmboe [78] studied the model of Taylor and Goldst+ with density
distribution (6.43) by his method of symmetric waves. Holmboe also con-
sidered density distribution (6.40)with a single discontinuity at y = 0.
0.8
0.7 -
a
0.5 -
0.4
- UNSTABLE
0.3-
0.2-
JO
Then the eigenvalue relation [76, equation (7.6) essentially] can be shown
to be
(5.47) 4 a V - c a{(2a - 1 ) 2 - ,-*a + +
4aJ0} (Jo/a)(2a- 1 e-k)a = 0. +
I t follows that ca is complex only when
(2a - 1 +3 e - 9 - 2{2e-k(2a - 1 + e-h)]l/g
-
1 , 1 , 1 , 1 , I . I
0 0.2 0.4 0.6 0.8 1.0 1.2
JO
cr> 1)
and J,, gLs/V*,s<< 1, it can be shown that
(5.60) (2c2 + 1 - J o / a ) g- 4c2(1 - e-kl) = e-"(1 - Jo/a).
Howard [26] found stability with four real roots c when
(5.61)
JO
FIG. 8. Stability characteristics for the double shear layer w = y/!y)(lyl> l),
I= O(lyI < 1) with density = pa(’ + e ) ( y c - l), = po(lyl < I), p =
po(l--E)@> 1).
(g) Sicliley Jet
Howard and Drazin [83] have found various exact parts of the neutral
eigensolutions for
(6.52) w = sechly (- oo<y< oo),
with various density distributions.
When J ( y ) is constant, the following solutions may be verified.
(5.52i) Sinuous mode c = (6 + as)/16, J = aa(4- aa)(Q- aa)/225,
p = (sechsy - c)h(sechy)m (0< aa < 4)
where
k +
3(4 - aa)/2(6 aa), m -= 6aa/(6+ aa).
(6.62ii) Varicose mode c = (3 + aa)*/3(3+ 5x3,
J = aa(l - a*)(9- aa)(3+ a*)”9(3 + 6aP)%,
tp = tanh y(sech’y - c)*(sechy)” (0 <aa < 1)
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 76
0 0.05
JO
0.1 .K)9 .I27
q = (sechy)=(sechPy- +a)*- (1/2)a
Fig. 9a.
(0 < a < 2).
(6.63ii) Varicose mode
+
c = (3 aa)/3(1 a). + STABLE
Jo = a ( l - a)
* (3 - a)(3 +
aa)/Q(l a ) , + 1.5 -
q = tanh y(sech y)'
.(sechay - c)(l@)(*-a)(0 < a 1). < UNSTABLE
JO
FIG.10. Stability boundary of one mode for to = sech y tanh y, J = J,, sech*y ( - 00
< y < do): Ja = a9(9 - a9)/9.
,
It gives the stability boundary, shown in Figure ll(a).
a
0 UNSTABLE
0.I 0.2
J JO
FIG. Illa). Stability boundary for w = tanh y, J = constant (- m < y < a): J =
= aa(1 - a*) (b). Stability boundary for w = tanh y. J = Josech3 y(- ao< y < ao):
Jo = a(l - a).
STABLE
0 1 2 3 4 1 6
JO
FIG. 1 l(c). Stability boundary for ru = tanh y , J = 3J0 sech' y tanh9 y ( - 00 <y< 00) :
+ +
J o = a(a 3)/3 and J o = (a - l ) ( a 2)/3.
1. Discussion
Using (6.2) to eliminate v' and (6.3) to eliminate zb', (6.1) isreadily transformed
into:
This is of the same form s(w - c)aQ = 0 with Q > O as in the plane parallel
case, and the semicircle theorem thus follows immediately, as before.
+
at4;lat -I- vauzi/ax u;du/ar = - apilax,
aw;lat + uau,i/ax= - apt/&,
a q a t + uau;lax = - aptIraq,
a q a x + a(Yu;)lra + auci/ray = 0.
-
and G(0) = H ( 0 ) = 0; and thus we may use condition (6.9) at Y = 0 except
when n = 0. In fact it can be seen that in general equation (6.8) gives
G N constant x r"-' (n# 0) and G constant x Y (n = 0) as Y + 0.
If U or DU is discontinuous, at yo say, then the pressure must be con-
tinuous a t the material interface with mean position r = y,. It follows that
(6.10) [(U- c)D(YG)- (DU)(rC)]= 0 (r = Yo).
82 P. G. DRAZIN AND L. N. HOWARD
Note that when U is piecewise constant it is easier to work with the amplitude
@ = JGdr of the velocity potential rather than with G directly [W].
The eigenvalue problem ( 6 4 , (6.9) is essentially due to Rayleigh [9l].
It has prompted surprisingly little later work, in view of the scores of papers
on the analogous problem (2.11), (2.12) of plane parallel flow. Perhaps the
similarity of the two problems and their methods of solution has discouraged
duplication of work, perhaps the greater physical importance of plane parallel
flows has overshadowed that of round jets. However, there is one important
difference between the two eigenvalue problems, namely the essentially
three-dimensional nature of instability of a round jet. Experience of plane
parallel flows suggests that varicose instability (axisymmetric disturbances
with n = 0) of a round jet should be less than sinuous instability (n = l),
so it comes as no surprise to find that there is no analogue of Squire’s the-
orem. In fact, in a recent examination of non-axisymmetnc disturbances,
Batchelor and Gill [go] found that a certain jet is most unstable to the
mode t z = 1.
The eigenvalue problem ( 6 4 , (6.9) is symmetric in a and (- a), so we
can again take a 2 0 without loss of generality. There is also a symmetry
in G,c and G*,c* for the same a, so c is real for stability and complex for
instability. Again, we write c, >0 when there is instability, bearing in mind
the initial-value problem and the inviscid limit of the viscous problem
P21, C90, §23.
Rayleigh [91] found a necessary condition for instability, analogous to
there being a point of inflexion in the velocity profile of a plane parallel
basic flow. Essentially by multiplying the stability equation (6.8) by
&*/(U - c), integrating from Y , to Y,, and taking the imaginary part, he
found that
(6.12) ci lglzDQdr = 0,
TI
+
where g E rG/(U - c ) , Q = r(DU)/(na a%*). Therefore a necessary
condition for instability (ci >O) is that DQ = 0 somewhere in the field of
flow. This is equivalent to U having a point of inflection with respect to the
+
variable p = s(n* aara)/rdr= n8 log Y +
)a%*. This reduces to Ray-
leigh’s condition for plane parallel flow if one regards the round jet as being
plane parallel flow locally when y1,r2-+ 00 and Y, - y1 = y , - y1 is fixed.
The following general stability characteristics are due to Batchelor and
Gill [go], who give details of the proofs.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 83
(6.13) 1
Is
rl
( g ( W- c,)DQdr < 0.
Therefore, when ci > O . equation (6.12) gives
(6.14)
It can be shown that the Reynolds stress tensor, averaged over one
period 2n/a of x and one 2n/n of v, has orthogonal components
(6.16) tlr
~
(6.18)
-
- nc,(DU)lGla
(6.19)
2a(n2 + aPra)l/z{(U - c,), +cia}
The stability equation (6.8) gives
2c,ra[GlaDQ
(6.20) DW=-
(U- c,)' + Cj' *
Therefore, if
the solution (rC) of the stability equation will be monotonic and cannot
satisfy both boundary conditions. Thus a necessary condition for the exist-
ence of the singular neutral solution is that f i is not so large that
max {rDQ/(U,- U)}< 1. In fact this condition is quite restrictive.
Very few examples have been treated in- the literature. First we take
the exact solution of the Navier-Stokes equations for a viscous fluid as our
basic flow, namely Poiseuille flow in a pipe with
(6.23) U = Ara + B log r + C (rl < r < r2).
Rayleigh [91] investigated the stability of this basic flow in an inviscid
+ +
fluid. I t gives Q = (2Br4 C ) / ( n 2 a W ) , which varies monotonically
with r . Therefore the flow is always stable.
For the cylindrical vortex sheet,
(6.24)
Batchelor and Gill [go] used the velocity potential on each side of the
discontinuity to deduce that the eigenvalue is
(6.26) c = (1 + iLnl’z(a)}/(l + L ( a ) } ,
where &,(a) - K,,(a)I,,’(a)/K,,’(a)l,,(a)in terms of the modified Bessel
functions I,,,K, of the first and second kinds and their derivatives. This
flow is unstable for each pair of values %,a. As a 00 (i.e. as the radius
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 86
Refemaces
1. HELMHOLTZ, H.,tlber diskontinuierliche Flussigkeitsbewegungen, Monats. konigl.
fireuss. Akad. Wiss. Berlin 215-228 (1868); translation by F. Guthrie. On
discontinuous movements of fluids, Phil. Mag. (4) 86, 337-346 (1868) ; also
,.Wissenschaftliche Abhandlungen“. vol. 1, J. A. Barth, Leipzig. 1882.
2. KELVIN,W., The influence of wind on waves in water supposed frictionless. Phil.
Mag. (4) 42, 368-374 (1871); also ”Mathematical and Physical Papers,” vol. IV,
“Hydrodynamics and General Dynamics,” pp. 76-83, Cambridge Univ. Press, 1910.
3. RAYLEIGEL J. W. S., “The Theory of Sound” vol. 11, chap. XXI. Dover, New York,
1945 (reprint of 2nd ed. of 1894).
4. SQUIRE,H. B., On the stability of three-dimensional disturbances of viscous flow
between parallel walls, Proc. Roy. SOC.A, 142, 621-628 (1933).
5. LIN. C. C., “Hydrodynamic stability,” Proc. Symposia in Appl. Math., vol. V,
McGraw-Hill, New York, 1-18, 1954.
6. HEISENBERG,W., uber Stabilltiit und Turbulenz von Fliissigkeitsstr6men. Ann.
Phys., Leipzig (4) 74, 577-627 (1924); translated as Nat. Adv. Comm. Aero.,
Washington. Tech. Memo., no. 1291 (1951).
7. LIN, C. C., On the stability of two-dimensional parallel flows. Part 11: Stability
in an inviscid fluid, Quart. AHZ. Math. 5, 218-234 (1945).
8. DRAZIN,P. G. and HOWARD, L. N., The instability to long waves of unbounded
parallel inviscid flow, J. Fluid Mech. 14, 257-283 (1962).
9. LIN, C. C., “The Theory of Hydrodynamic Stability,” Cambridge Univ. Press,
London and New York (1955).
10. STUART, J . T., Hydrodynamic stability, i n “Laminar Boundary Layers,” chap. IX,
(L. Rosenhead, ed.), Oxford Univ. Press, London and.New York, 1983.
11. REID, W. H.. The stability of parallel flows, i n “Basic Developments in Fluid
Dynamics,” (M. Holt, ed.) Academic Press, New York, pp. 249-307, 1965.
12. RAYLEIGH, J . W. S., O n the stability, or instability, of certain fluid motions, PYOC.
Lond. Math. SOC.9, 57-70 (1880); also “Scientific Papers,’’ vol. I. pp. 474-487
Cambridge Univ. Press, London and New York. 1899.
13. FJBRTOFT. R., Application of integral theorems in deriving criteria of stability of
laminar flow and for the baroclinic circular vortex, Geofys. Publ. 17, no. 6
pp. 1-52 (1950).
14. TOLLMIBN, W.. Ein allgemeines Kriterium der Instabltiit Laminarer Geschwindig-
keitsverteilungen, Nachr. Ges. Wdss. Gottingen, Math.- Phys. Klasse, 60, 79-1 14
(1935); translated as Nat. Adv. Comm. Aero., Washington. Techn. Memo.. No. 792
(1936).
86 P. G. DRAZIN AND L. N. HOWARD
16. FRIEDRICHS, K. O., “Fluid Dynamics,” chap. IV. mimeographed lecture notes,
Brown University, Providence, Rhode Island, 1942.
16. HOWARD, L. N., The number of unstable modes in hydrodynamic stability problems.
J . de Mdcanique 8, 433-443 (1964).
17. ROSENBLUTH, M. N. and SIMON,A., Necessary and sufficient condition for the
stability of plane parallel inviscid flow. Phys. Fluids, 7, 667-668 (1964).
18. FOOTE,J . R. and LIN, C. C., Some recent investigations in the theory of hydro-
dynamic stability, Quart. Appl. Math. 8, 266-280 (1950).
19. TAYLOR,G. I., Eddy motion in the atmosphere, Phil. Trans. Roy. Soc. A, B1Q.
1-26 (1816).
20. LIGHTHILL, M. J., Physical interpretation of the mathematical theory of wave
generation by wind, J . Fluid Mech. 14, 386-398 (1962).
21. KELVIN,W., On a disturbing infinity in Lord Rayleigh’s solution for waves in a
plane vortex stratum, Nature 88, 45-46 (1880): also “Mathematical and Physical
Papers,” vol. IV., Cambridge Univ. Press, pp. 186-187, 1010.
22. GILL, A., A mechanism for instability of plane Couette flow and of Poiseuille
flow in a pipe, J . Ffuid Mech. 81, 603-611 (1966).
23. HOWARD, L. N.. Note on a paper of John W . Miles, J. Fluid Mech. 10, 609-612
(1961).
34. HBILAND,E.,On two-dimensional perturbation of linear flow, Geofys. Publ. 18,
NO. 9 pp. 1-12, (1863).
26. HOWARD.L. N., Neutral curves and stability boundaries in stratified flow, J.
Fluid Mech. 16, 333-342 (1963).
26. ORR,W. M’F.. The stability or instability of the steady motions of a perfect liquid
and of a viscous liquid. Part I : A perfect liquid, Proc. Roy. Irish Acad. %7 A,
9-68 (1907).
37. MILES, J. W., On the disturbed motion of a vortex sheet, J. Fluid Mech. 4, 638-662
(1968).
28. CARRIER, G. F. and CHANG.C. T., On an initial vaiue problem concerning Taylor
instability of incompressible fluids, Quart. ApPI. Madh. 18, 436-439 (1969).
20. CASE,K. M.,Stability of inviscid plane Couette flow, Phys. Ffuids 8, 143-148 (1960).
30. CASE, K.M.,Stability of an idealized atmosphere, I: Discussion of results, Phys.
Ffuids 8, 149-164 (1960).
31. DIKII,L. A., On the stability of plane parallel flows of an inhomogeneous fluid,
J. Appl. Math. Mech. 84. 367-369 (1960). This is English translation of Russian
paper in Prik. Math. Mekh. B4, 249-260 (1960).
32. DIKII.L.A.. The stability of plane-parallel flows of a n ideal fluid, Sou. Phys. Dohlady
6, 1179-1182 (1961). Translated from Russian in Dokkrdy Akad. Nauk. S S S R
186, 1068-1071 (1960).
33. CASE, K. M., Hydrodynamic stability and the inviscid limit, J . Fluid Mech. 10,
420-429 (1961).
34. LIN, C. C., Some mathematical problems in the theory of the stability of parallel
flows, 1. Fluid Mech. 10, 430-438 (1961).
36. ELIASSEN,A.. HeILAND, E.. and RIIS, E., Two-dimensional perturbation of a flow
with constant shear of a stratified fluid, Inst. Weather Climate Res., Oslo. hb-
lication No. 1 (1063).
36. SAVIC,P.,O n acoustically effective vortex motion in gaseous jets, Phil. Mag. ( 7 )
8g, 246-262 (1941).
37. SAVIC.P. and MURPHY, J. W., The symmetrical vortex street in sound sensitive
plane jets, Phil. Mag. (7). 84, 139-144 (1943).
38. L~SSEN, M. and FOX,J. A.. The stability of boundary layer type flows with infinite
boundary conditions, ,,60 Jahre Grenzschichtforschung,” (H. Wrtler and
W. Tollmien, eds.), Vieweg und Sohn, Braunschweig, pp. 122-126, 1965.
HYDRODYNAMIC STABILITY OF PARALLEL FLOW O F INVISCID FLUID 87
39. HOLLINGDALB, S., Stability and configuration of the wakes produced by solid
bodies moving through fluids, Phil. Mag. (7)99, 209- 57, (1940).
40. HAURWITZ, B. and PANOFSKY, H. A., Stability and meandering of the Gulf Stream,
Trans. Amer. Geophys. Union 81, 723-731 (1960).
41. SATO,H.. The stability and transition of a two-dimensional jet, J. Fluid Mech. 0 ,
53-80 (1960).
42. CURLE, N., On hydrodynamic stability of unlimited antisymmetrical velocity
profiles. Aero. Res. Council, London, unpublished Rept. No. 18664 (1956).
43. GARCIA,R. V., Barotropic waves in straight parallel flow with curved velocity
profiles, Tellus 8, 82-93 (1966).
44. CURLE,N.. Hydrodynamic stability of the laminar mixing region between parallel
streams, Aero. Res. Council, London, unpublished Rept. No. 18426 (1956).
45. MICHALKE.A., On the inviscid instability of the hyperbolic-tangent velocity profile.
J. Fluid Mech.. 19, 643-656 (1964).
46. ESCH,R. E., The instability of a shear layer between two parallel streams, J . Fluid
Mech. 8, 289-303 (1957).
47. CHIARULLI,P.and FREEMAN, J. C.. Stability of the boundary layer. Headquarters
Air Materiel Command, Dayton, unpublished Tech. Hept. No. F-TR/l197-1A
(1948).
48. BJERNBS. V., BJERKNES. J., BBRGERON, T.. and SOLBERG.H., .,Physikalische
Hydrodynamik,” Springer, Berlin, 1933. Translated into French as “Hydro-
dynamique Physique.” Presses Universitaires de France, Paris. 1934.
49. GODSKE,C. L., BLRGERON, T.. B~ERKNES, J. and BUNDGAARD, R. C., “Dynamic
Meteorology and Weather Forecasting,” Amer. Meteor. SOC., Boston, 1967.
60, HAURWITZ. B., Zur Theorie der Wellen in L u f t und Wasser, Veriiff. Geophys. ins;.
Llniv. Leipsig, 6, no. 1 (1931).
51. RAYLEIGH, J. W. S., Investigations on the character of the equilibrium of an in-
compressible fluid of variable density. Proc. Lond. M d h . SOC.14, 17&177 (1883).
Also in “Scientific Papers,” vol. 11, pp. 200-207, Cambridge Univ, Press,London
and New York, 1900.
62. TAYLOR, G.I., Effect of variation of density on the stability of superposed streams
of fluid, Proc. Roy. Soc. A, 183. 499-523 (1931).
53. GOLDSTEIN, S.. On the stability of superposed streams of fluid of different den-
sities, Proc. Roy. SOG.A. 182. 524-548 (1931).
64. YIH, C. S.. Stability of two-dimensional parallel flows for three-dimensional dis-
turbances, Quart. A p p l . Math., 12, 434-435 (1965).
65. TAYLOR,G. I., The instability of liquid surfaces when accelerated in a direction
perpendicular to their planes. I, Proc. Roy. SOC.A, W l , 192-194 (1950).
68. ALTERMAN,2.. Kelvin-Helmholtz instability in media of variable density, Phys.
Fluids, 4, 1177-1179 (1961).
57. MENKES,J., On the stability of a shear layer, J. Fluid Mech., 6, 518-622 (1969).
68. FEJER,J. A.. and MILES. J. W.. On the stability of a plane vortex sheet with respect
to three-dimensional disturbances, J. Fluid Mech., 16, 335-336 (1963).
59. LANDAU, L., Stability of tangential discontinuities in compressible fluid, Akad.
N a r k . S S S R . Comptes Rendus (Doklady), 44, 139-141 (1944).
60. HATANAKA, H., O n the stability of a surface of discontinuity in a compressible
fluid, J . SOC.Sci. Culture, Japan, e, 3-7 (1947).
81. LIEPMANN. H. W.and ROSHKO, A., ”Elements of Gas Dynamics,” Wiley, New York.
1957.
62. K ~ C H B M A ND., N , Sarungsbewegungen in einer Gasstr6mung mit Grenzschicht.
ZCit. angew. Math. Mech. 18. 207-222 (1938).
83. LIN, C. C., On the stability of the laminar mixing region between two parallel
streams in a gas, Naf.Adv. Comm. Aero., Washington, Tech. Note No. 2087 (1968).
88 P. G. DRAZIN AND L. N. HOWARD
64. ECKART. C., Extension of Howard’s circle theorem for adiabatic jets, Phys. Fluids
6, 1042-1047 (1963).
66. JOHNSON, J . A., The stability of a shearing motion in a rotating fluid, J. Fluid
Mech. 17, 337-362 (1963).
66. Kuo, H. L., Dynamic instability of two-dimensional non-divergent flow in a baro-
tropic atmosphere, J . Meteor., 6 , 105-122, (1940).
67. HOWARD. L. N. and DRAZIN,P. G., On instability of parallel flow of inviscid fluid
in a rotating system with variable Coriolis parameter, J . Math. and Phys. 48,
83-99. (1064).
68. LIPPS, F. B., The barotropic stability of the mean winds in the atmosphere, J.
Fluid Mech. 12. 397-407 (1962).
69. MICHAEL, D. H.. Stability of plane parallel flows of electrically conducting fluids,
Proc. Camb. Phil. Soc., 49, 166-168 (1953).
70. STUART, J . T., On the stability of viscous flow between parallel planes in the presence
of a coplanar magnetic field, Proc. Roy. SOC.A, 391. 189-206 (1964).
71. MICHAEL, D. H., The stability of a combined current and vortex sheet in a perfectly
conducting fluid, Proc. Camb. Phil. SOC.,61, 628432 (1956).
72. DRAZIN, P. G., Stability of parallel flow in a parallel magnetic field at small magnetic
Reynolds numbers, J . Fluid Mech., R, 130-142 (1960).
73. LOCK,R. C., The stability of flow of an electrically conducting fluid between parallel
planes under a transverse magnetic field, Proc. Roy. SOC.A, 388. 10.5-125 (1965).
74. LIEPMANN, H. W. and PUCKETT, A. E., “Introduction to Aerodynamics of a Com-
pressible Fluid,” Wiley, New York, 1947.
76. HOLYBOE,J., On the behavior of symmetric waves in stratified shear layers.
Geofys. Publ., 94, no. 2, 67-113 (1962).
76. SYNGE,J . L., The stability of heterogeneous liquids, Trans Roy. SOC. Canada, 27,
111, 1-18 (1933).
77. MILES, J. W., On the stability of heterogeneous shear flows, J. Fluid Mech., 10,
496-508 (1961).
78. MILES, J. W., On the stability of heterogeneous shear flows. Part 2. J . Fiuid Mech.,
16, 209-227 (1963).
79. DRAZIN, P. G. and HOWARD, L. N., Stability in a continuously stratified fluid,
Proc. Amer. SOL. Civil Engvs. E . M . Div. R7, 101-116 (1961). Also Trans A m r .
SOC. Civil h g r s . 19A. 849-864 (1963).
80. FJELDSTAD, J . E., Interne wellen, Geofys. Publ., 10, no. 6, pp. 1-35 (1933).
81. DAVIS,P. A . and PATTERSON, A. M., The creation and propagation of internal
waves. A literature survey. Unpublished tech. memo No. 56-2, Pacific Naval Lab..
Esquirnault, British Columbia (1956).
82. HeILAND, E., On the stability of Couette-flow of a fluid with parabolic gravitational
stability, Vitensk. Akad. Inst. Vaer-og Kiimasforskning. Publ. no. 1, (1954).
83. HOWARD, L. N. and DRAZIN, P. G., (to be published).
84. MILES, J. W. and HOWARD, L. N., Note on a heterogeneous shear flow, J. Fluid
Mech., 90, 331-336 (1964).
86. DRAZIN, P. G., The stability of a shear Layer in an unbounded heterogeneous in-
viscid fluid, J . Fluid Mech. 4, 214-224 (1958).
86. ESCH,H. E., Stability of the parallel flow of B fluid over a slightly heavier fluid,
J. Fluid Mech., 12. 192-208 (1962).
87. HOCKING, L. M . . The instability of a non-uniform vortex sheet, .I, Fluid Mech.,
l A , 177-186 (1964).
88. HOCKING,L. M., Note on the instabilitv of a non-uniform vortex sheet. J . Fluid
Mech., 31, 333-336 (1966).
89. SYNGE,J . L., Hydrodynamical stability, Semi-centenn. Publ. Amer. Math. Sac..
S, 227-269 (1038).
HYDRODYNAMIC STABILITY OF PARALLEL FLOW OF INVISCID FLUID 89
90. BATCHELOR. G. K. and GILL,A. E., Analysis of the stability of &symmetric jets,
J . Fluid Mech. 14. 629-551 (1962).
91. F~AYLBIGH,J. W. S., On the question of the stability of the flow of fluids, Phil.
Mag., 84. 59-70 ( 1 8 9 2 ) ; also “Scientific Papers”, vol. 111, Cambridge Univ.
Press, London and New York, pp. 676-584, 1902.
92. PRETSCH, J . , Uber die Stabilitat einer LaminarstrBmung in einem geraden Rohr
mit kreisf6rmigem Querschnitt, Z . angew. Math. Mech., el, 204-217 (1941).
The Calculation of Free Oscillations of a Liquid in a Motionless
Container
BY N. N. MOISEEV AND A. A. PETROV
I. INTRODUCTION
Only in the simplest cases one can manage to solve it analytically by the
Fourier method. These are the generally known oscillation problems of
a liquid in a container of cylindrical shape with a flat bottom (parallelepiped,
circular cylinder, cylinder with a circular ring or its sector as the base, and
elliptical cylinder). Apart from this, there are some papers which, in one
way or another manage to give solutions of the problem of liquid os-
cillations in a container of a certain special shape. These papers include
for example [2] and [3].
However, a number of engineering tasks require the skillful construction
of the spectrum of natural frequencies for containers of an almost com-
pletely arbitrary shape. Therefore a method of solution of these prob-
lems is required which can claim a good degree of universality in its
application. One of the possible methods satisfying this requirement is,
in the opinion of authors, the Ritz method. I t is based on the application
of Hamilton’s principle to the motion of an ideal incompressible liquid.
In paper [l] the applicability of the Ritz method has been proved for a
large class of problems.
The present paper can be considered as a continuation of this work.
Here, on the basis of the general results given in [l], we systematically
describe the application of the Ritz method to calculations of free oscilla-
tions of a liquid in containers of various shapes.
The present survey consists of two sections. In Section 1.1 we have
considered it necessary to give the basic results of the general theory. This
simplifies further reading of the material. In Sections 1.2 and 1.3 some
general thoughts of efficient realization of the Ritz scheme are given.
In the second section the solutions of practical problems of free oscil-
lations of liquid in containers of various shapes and the results of the cal-
culations are presented.
Finally, in the Appendix one of the methods of numerical solution of
an algebraic problem is discussed which arises as a result of the application
of the Ritz method.
A significant part of the results has been obtained by the authors of the
survey.
In addition, results are given which have been obtained by other authors
who solved problems by the variational method. However, here we have
used only the results of their calculations. As for the method, we have
adopted a uniform procedure for solving all the problems.
This survey considers only small oscillations of the liquid under the
action of mass forces. Many questions which by now have become actual
and to which much work has been devoted, are beyond the scope of this
paper. These are first and foremost the problems of liquid dynamics in
force fields of low intensity where surface tension begins to play a decisive
pa1 t, the dynamics of the liquid-containing body under these conditions,
and, finally, the dynamics of a body containing a viscous liquid.
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 93
P i
/-
x
FIG.1.
area) be C. Introduce the coordinate system Oxyz with the xy-plane parallel
to the free surface of the liquid; the z-axis is opposite in direction to
the vector g in Fig. 1.
As has been shown in [l], in the problems considered the liquid motion
can be considered as irrotational without loss of generality. Denote the
velocity potential by y(x,y,z,t) and the elevation of the free surface above
the plane S by c(x,y,t). If one assumes that the motion of the liquid is small,
then the function q~ is a solution of the following boundary problem [l]:
(1.1) Av=O in the region t,
_
av -- 0 on the surface C,
(1.2) an
on the surface S;
2. One calls free (or natural) oscillations the solution of the problem
(1.1)-(1.4) of the form
gJ(x,y,z,t)= @(x,y,z)cos d.
S(x,y,t)= Z(x,y)sin d,
where (I is the natural frequency. From the conditions (1.1)-(1.4) it follows
that the function 0 is a solution of the following eigenvalue problem:
(1.6) A@=O in the region 7,
-a@an
=o on the surface Z,
a@
~ = i Z Q , on the surface S.
(1.7) an
(1.9) ! I
F(@) = (I'@)'dz - 1. W d S .
s
where
(1.12)
(1.13) 2 p,am
m-1
-1 2 q,,,,,a,,,= 0,
m-1
12 .
= 1,2,. ..N.
96 N. N. MOISEEV AND A. A. PETROV
(1) If the functions x,, are harmonic and satisfy some of the boundary
conditions of the problem, one obviously will expect that the Ritz procedure
will rapidly converge. And one of the widely applicable ways of choosing the
set x,, consists in the construction of a set of linearly independent solutions
of Laplace's equation (1.6) in the region T . which satisfy the conditions
(1.6) as accurately as possible.
In Sections 11.1-11.4 solutions of a number of problems are given
which are based on the use of this idea.
In particular, if the problem (1.6)-(1.7) permits partial separation of
variables, i.e., if the solution of the problem is, let us say, of the form
To fix the ideas, assume that a container (the region t) has the shape
of a vertical cylinder whose height is h. Then, by substituting the expression
(1.16) into the functional (1.9) and by integrating with respect to z we
obtain the functional
(1.16)
s s
where
(1.17)
U(%) =
+
cosh o ( z h)
cosh o h
A=otanhoh.
Thus, for the case of a right cylinder with a flat bottom, we have a
variational problem for the functional (1.16), which is much simpler than
the initial one.
the eigenvalues by p,,,. Put the container, for which the solution of the
problem of free oscillations is required, wholly inside the region toof minimal
dimensions in such a way that the planes of areas S and So coincide (Fig. 2).
SO
I
r I
I
'\= //
\
20
FIG.2.
Clearly, from the regions tomentioned above there should be chosen one
whose shape is closest to the shape of the container, but here, certainly,
much arbitrariness remains. The set of the functions ($I~} can be taken
as coordinate functions by setting xrn= +,,,.
Denote the difference of volumes r0 and t by A t and the difference of
areas So and S by A S . The formulae (1.12)for the coefficientsof the system
(1.13) can now be rewritten in the form
$urn = j
Trn
v+n v+rndt - 1
AT
v+n v+rndt,
j v+nv+rndt = pnj + n + ~ s -
7a sa
This implies, because of the equivalence of the indices rc and m, that the
functions +,, are mutually orthogonal over the surface So. If we consider
that
5
sb
= 1,
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 99
Furthermore, from the same Green's formula and from the boundary
values of the functions $, it follows that
Hence, the coefficients of the system (1.13) can be represented in the form
convenient for caIculations
(1.18)
(1.19)
(1.20)
If the shape of a container is such that the free surfaces So and S co-
incide, then AS = 0 (cf. Fig. 2) and the equations (l.lS), (1.19) simplify to
qnm = Sum.
x = xl, Y =y d , z = ad.
FIG. 3.
where
and
where
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 101
2. Right circdar cylider. Denote the radius of the cylinder by R and the
height by H. Introduce the coordinate system as shown in Fig. 4, using
dimensionless variables
I I li
FIG.4.
where Jn(E) are Bessel functions of the first kind, KZ)is the mth root of
the equation J,'(K) = 0, and
The case n =0 is the exception; here the eigenvalues have multiplicity one.
3. Right cylinder with an anmdar base. Denote the radius of the inner
cylinder by R, and the radius of the outer cylinder by cR where c > 1 is a
102 N. N. MOISEEV AND A. A. PETROV
FIG.5 .
(1.24)
where
P,(r) = N,,'(K,,,("))J,,(K,,$"Y)
- J.'(K,,,(,,))N,,(K,,,(")Y),
J,,(E) and N,,(E) are Ressel functions of the first and second kind, K,,,('" is the
mth root of the equation
J:(K) .N,,'(CK)- N / ( K ) ' J:(CK) = 0,
and
m,n = 0,1,2,. . ..
The exception is again the case n = 0 with eigenvalues of multiplicity one.
Assume now that the region z, for which the problem of free oscillation
of a liquid is to be solved, is close to the region to(we assume here that the
planes of the free surfaces S and So coincide) for which the effective solution
* The main content of this section is discussed in [a].
F R E E LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 103
of a similar problem can be constructed. This means that for towe know
the system of the eigenvalues pmand the eigenfunctions #m. The closeness
of the regions t and to(S and So respectively) permits an efficient use of
perturbation theory.
In solving the problem for the region t by the Ritz method we shall
take the functions +m as a system of coordinate functions. It is assumed in
this case that the functions (Im can be continued outside the region toin
such a way that the system of these functions remains complete also in the
region t (unless z, completely encloses t).
Let us.now find the solution in the form
N
= 2
n=1
an+nq
This will lead us to the system (1.13) whose coefficients in this case will be
of the form
1 1
70 AT
But it has been shown in the previous section that for functions normalized
so that
the relations
hold.
Furthermore, since the
(1.25)
(1.26)
where E << 1 and the numbers l,,,,,and q,,,,, are of the order of unity.
104 N. N. MOISEEV AND A. A. PETROV
(1.30)
m
(1.31)
By substituting the series into system (1.29) and by equating the coefficients
of the same powers of 8 we obtain the following systems for the determina-
tion of the coefficients of the decomposition (1.30) and (1.31):
(1.32) (Pn - A('))U,(~) = 0, n = 1,2,. . ., N ,
The coefficient a:) has remained indeterminate. We shall find it later from the
conditions of normalization.
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 106
1
obtained in the form
v$sv$ndz - ys j+
$s$ssds
+c
N
@, = (1 + d ) ) * S
n=1
AT
Ps
AS
- Pn
- *us
U#S
These formulae can be simplified by using Green’s formula and the boundary
values of the function $,. Indeed, if
5
n#s
I , = ys - *s $dS.
z
Let us finally, determine the coefficient a$:’ from the normalization condition
I
s
asSdS = 1.
(1.37)
106 N. N. MOISEEV AND A. A. PETROV
111.APPLICATION
OF THE RITZMETHOD TO THE CALCULATIONOF NATURAL
OF A LIQUIDIN CONTAINERSOF VARIOUSSHAPES
OSCILLATIONS
FIG.6,
In the new coordinate system the axis of the cylinder will appear vertical
and the free surface of the liquid will be inclined at the angle ,8 with respect
to the plane xy (Fig. 7). Introduce the cylindrical coordinates
FIG.7.
b 0
I/ 5
2n
-h
1 r(r,e)
Here z(r,B) = x tan p = r cos 0 tan p is the equation of the free surfaces
h = H / R is the dimensionless average liquid depth and the number,
1 = aaR/g are the dimensionless eigenvalues of the problem.
Let us take as a set of coordinate functions those functions which des-
cribe the oscillations of the liquid in the vertical cylinder of depth h [see
(1.23) in Section 1.21:
1
+.1((0 = N,c", J,,(K,(")Y) sin 120
+
cosh K,(")(Z 12)
s,n = 1,2,. . .,
cosh ~,(")h '
These are harmonic functions and satisfy all the boundary conditions of
our problem except the condition at the free surface.
The solution of the variational problem (2.1) will be found in the form
N N
#= 2 aju)UJu)sin n6 + b , ( W j u ) cos ne,
s,n s,n
108 N. N. MOISEEV AND A. A. PETROV
where
p$"' 55 I
2n
= d8 rdr
0 0
1 rcosBtanD
-h
(VUJ") * VULm)sin n8 sin m8
nm
J J
0 0
'%n !
I v1 + tanaPU,c")(r,z(r,B))V,rm)(r,z(r,B))
cos n8 cos me&,
0 0
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 109
nm
- -U p )U,(m)cos ne sin me
1%
In our case
n 1
n 1
110 N. N. MOISEEV AND A. A. PETROV
N N
It is not difficult to show that for vertical cylinder axis (p = 0 ) ji$"") = fifi'"),
qLM) = &"') and the solutions of the systems coincide. This corresponds
to the fact that the eigenvalues have multiplicity two, and to each of them
belong two eigenfunctions and I,&(") (except for the case n = 0). If the
" t 3 5 7 9
N
Here are some results of the numerical calculations. The systems (2.2),
(2.3) have been solved by the escalator method. The number of equations
in the systems N was taken equal to 3,6,7,9. The dependence of the value
on the approximation number is shown in Fig. 8 for two values of the
angle f l : p = 30' (dotted curve) and fl = 80" (solid curve).
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 111
1.0
10 20 30 40 50 60
Po
FIG.9.
TABLE1
8"
Let the lateral surface of the container be a circular cone with the semi-
angle 8, and the bottom of the container a spherical cap of radius R.
Introduce the Cartesian coordinates O X Y Z with origin at the vertex of the
cone and with the axis Z directed downwards along the cone axis (the direc-
tion of the axis Z therefore coincides with the direction of g). We shall measure
the depth of the liquid in the container through the distance H from the
vertex of the cone to the free surface S (Fig. 10).
FIG.10.
Remember that the problem of free oscillations in the given region reduces
to the following eigenvalue-boundary problem :
(2.4) A@ = 0 in the region r,
(2.6)
a@ = 0
- for r = 1,
a@ - 0
-- for 6 = 8,,
aY a8
(2.6)
a@ = - A@ for z = k.
From the form of the boundary values (2.5), (2.6) one notices that the problem
admits the separation of variables @ = U(r,tl)V(q). Since the solution
must be periodical in q it is obvious that
Co6
I/= nq.
sin
I
i
FIG. 11.
where p is the value of the radius vector of a point of the free surface in polar
coordinates p,q on the free surface.
As can be seen from Fig. 11, the equation of the free surface has the form
Y cos 8 = h. On the free surface
p = h tan 6, dp = -de,
cosa 8
114 N. N. MOISEEV AND A. A. PETROV
and
au -
_
az
--au cos 8 - -
ar
aU sin
- ;;si;e)- .
0 = cos 6 - - -
a8 r
(y
By combining all these results the functional (2.10) can be reduced to the
form
.(
8,
au a ~ s i
-- n ~ ]
h ,o)
F(U)=
n
a8 Is r-h/care a sine de
c d e
(2.11)
sin 8
0
(2.13)
au
-= 0 for Y =1 and
au
-= 0 for 8=eW
ar ae
Thus, for the solution of the problem (2.11) by the Ritz method a system
of coordinate functions complete in the region 0 8 0, and satisfying < <
conditions (2.12), (2.13) must be constructed.
To do this, let us add to the conditions (2.12) and (2.13)
(2.14)
au
-= - oU for r = h,
ar
and consider the auxiliary eigenvalue problem (2.12)-(2.14) in the region
To (Fig. 11).
This problem can be easily solved by separation of variables, setting
U = a(r)p(O). By substitution of this expression in the conditions (2.12)-
(2.14) we obtain the following boundary problems for the functions a and B
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 116
(2.16) ye
d2a - 2y-da
~ - pa = 0 h< r < 1,
dY2 a7
da da
(2.16) --
-0 for r=l, - -- - w a for Y =h,
dY dr
_
dP -
-0 for 8 = 8,.
dE
Besides, for 8 = 0 the function p should be finite.
The solution of this problem is represented by the associated spherical
functions
Indices v, in this expression are chosen in such a way that for each integer n
110 N. N. MOISEEV AND A. A. PETROV
(2.23) u= c asUs
N
s= 1
0:
FIG.12.
N
FIG.13.
9X FIG.
O 14.
Y
The solution of the variational problem (2.11) in this case is again sought
for in the form given in (2.23) where the functions U s are calculated by
(2.26). The minimum conditions lead to the system (1.13) whose coefficients
are calculated by formulae similar to (2.24).
The results of calculations of the first natural frequency 1, of the single-
node oscillations (n = 1) are represented in Fig. 16, where the dependence
1.8
0.6
0.2
0:
FIG. 15.
of 4 on the angle 8 is given. Also, the experimental results taken from [8]
are represented by dots.
If one recalls the solution of the oscillation problem in the inclined cylinder
given in the previous section, one will agree that the oscillation problem
in inclined cones must be solved in a similar way. The set of functions
determined by (2.22) and (2.26) must be taken as coordinate functions.
FIG. 18.
aQ,
(2.27) --- 0 for r = 1,
ar
(2.28)
(2.29)
F(@)=
stc5: @-&-A
5
s
@ads
or
(2.30)
On the surface
In this formula
4
(2.33)
function containing the multiplicator cos nq and the other sin nq. The
exception is the case 4t = 0 , where the eigenvalues are simple. The axi-
symmetrical forms of liquid oscillations correspond to this eigenvalue.
The solution of the variational problem (2.33) is constructed by the Ritz
met hod.
FIG.17.
Um= vmPrnn(cos
6) m = 1,2,. . .
where Prn"(cos0) is the associated Legendre function of nth order, and to seek
the minimizing function U in the form
N
(2.34 u = C a,,,pPrn*(cose).
m = l
FIG. 18.
FIG. 19.
dots in the Figure. Fig. 19 shows the shape of the principal singIe-node
oscillation for h = 1.6.
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 123
FIG. 20.
vertically upwards (i,e. opposite to the vector g). The liquid depth H will
be measured from the plane X Y . Thus, H > 0 if the container is filled more
than half, otherwise H < 0 (see Fig. 20).
Introduce cylindrical coordinates R,6,X
X = X, Y = Rsin6, Z =- Rcos I?
a@
(2.37) -=0 for x = O and x = l ,
ax
(2.38)
a@
-= A@ for z = h,
az
124 N. N. MOISEEV AND A. A. PETROV
a@ -
_ -0 for s = 0.
(2.39)
as
As has been shown previously, the boundary problem of the kind (2.35)-
(2.39) is equivalent to the variational problem of minimization of the
functional (1.9). In the case under consideration this functional will have
the form
1 1
where
rl arccos [ - h/(r,+ 41
~ ( 6=) -
$0
FIG.21.
We shall again solve the problem (2.42) by the Ritz method. For the
realization of the Ritz method one should construct a complete set of co-
ordinate functions {x,,}.
For this purpose consider the following auxiliary boundary problem
for the region do:
(2.46) ax
-- for s = rl.
as - ox
126 N. N. MOISEEV AND A. A. PETROV
(2.47) ~-
dab
ass k2P = 0. dP
-= 0 for s = 0,
dP
-= o@ for s = Y,;
as as
(2.48)
(2.49) --
da - 0 for t=O and f = 1.
dE
The solution of the boundary problem (2.47) is of the form
cosh ks
(2.60) p=- cosh kr, ’ w = k tanh kr,,
where k is to be determined.
By transforming the variables, (2.48) can be reduced to Bessel’s equation
whose solution will be Bessel functions of imaginary order and argument.
However, since we need the set of coordinate functions for a numerical
solution of the variational problem it is quite sufficient for us to find an
approximate but simpler and more convenient solution of the problem
(2.48), (2.49).
For such a solution we will use the following ideas. Note that for rl-+oo
the boundary problem (2.48), (2.49) takes the form
where
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 127
(2.62) _
kv -
-
vz + 1/(. ). + (q)l
1
I
I----
2
\ ”
, v = 0,1,2,.. .
80 2
The physical meaning of the solution (2.51), (2.62) is quite clear. For
large values of the inner cylinder radius r, the wall curvature is small. The
entire motion is taking place in a layer of the order of a wavelength, i.e.,
of the order of 6, and the motion in the limit for r, 00 is similar to that
-+
in the space between vertical walls of infinite length. This motion is des-
cribed by a solution of the kind (2.51), (2.52).
Assume temporarily that the container is such that Y, is large and 811
is small, i.e., the length and the radius of the inner cylinder are larger than
the “width” 8.
It is physically clear that in this case the wavelength along the cylinder
axis will be large compared with the wavelength in the transverse direction.
Therefore, for not too large values of n, the number (r1/8,,)kcan be con-
sidered large compared with the number nndll. The equation (2.48) can then
be considered as an equation with the large parameter y = (rl/80)k which
makes it possible to apply asymptotic integration methods.
In [ l l ] the equation
where C, and C2 are arbitrary constants and yl(E) and y,($) linearly indepen-
dent solutions of the equation
where
.. .
(2.64)
Note that for rl -* do and small S/1 the formulae (2.53), (2.64) imply
(2.66)
where the functions wv(E) and the numbers k, are determined by (2.63)
and (2.64).
For h = 0 the expressions (2.56) give an approximate solution of liquid
oscillations in the region T,which is the more precise the greater the ratio 116.
The set of functions xv is complete in the region A , bounded by the
circumferential arcs and the radius 6 = 8,. The region A bounded by the
line z = h is completely immersed in the region A,, hence the set {xv}defined
by (2.65) is also complete in the region d and can be taken as coordinate-
function system in solving the variational problem (2.42)by the Ritz method,
setting
N
(2.66) u = v2- 0 avx..
F R E E LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 129
This set was used for the calculations carried out on a computer. As a result
the dependence of the eigenvalues on the parameters rl, It, and 118 has
been obtained and is shown in Figs. 22-24. In all computations n = 1 was
assumed.
h 0.0
h:-1.67
-= 5.0
FIG. 22.
These results indicate that for lhJ<Y, the solution converges rapidly. As
it was to be expected the convergence becomes worse if lhl> r,.
2.2
2.0
1.8
1.6
1.4
xt
1.2
1.0
0.8
0.6
0.4
0.2
FIG.25.
50
40
30
A,
20
10
5 10 15 :
-
I
S
FIG.24.
F R E E L I Q U I D OSCILLATIONS, MOTIONLESS CONTAINER 131
TABLE2
N
FIG.25.
TABLE3
3 4 5 6
FIG. 26.
coinciding with the cylinder axis, the plane YZ being one of the end faces,
the axis Z directed opposite t o the vector g. Again we shall measure the
liquid depth H from the plane XY as indicated in Fig. 26.
Here our dimensionless variables are
(2.67) (+) aw
a
w+4,,@=0 in the region t.
F R E E L I Q U I D OSCILLATIONS, MOTIONLESS CONTAINER 133
(2.58) -a@
=o
ax
for x = O and x = l ,
(2.59) -a@
=o
ar
for r = ro (r2= z2 + y2),
(2.60)
a@
-= A@ for z = h.
az
Using the Ritz method again, we construct the set of coordinate functions
by the procedure described in Section 1.2. For this purpose immerse the
container in a parallelepiped with length 1, height R H , and width 2d. +
For the set of coordinate functions we shall take that set which describes
the natural oscillations of the liquid in the parallelepiped. Again, it was shown
in Section 1.2 that these functions have the form
1
= N,c", cos mny
cosh K,,,(")(z ro)+
+("),,
cosh ~,,,(")(h + ro)'OS nnxJ
1 2m + 1 cosh R,,,(")(z + h)
&(") = N,o sin ~
(2.62)
We shall determine the solution of the variational problem (2.62) in the
form
N N
(2.63) @= 2 Urn(")+,,,(*) + 2 b,,,(")I?,,,("),
m.* m,n
where
differing from zero only in the case n = s. This means that the set (2.64)
decomposes into N independent sets which correspond to n = 0 , 1 , . . . , N .
The solution of the set corresponding to n = no describes oscillations
of the form @ = U(y,z)cos nonx which comprise no half-waves along the
axis x . Therefore we shall write from now on only one upper index 12 for all
the coefficients of the set (2.64).
It also follows without difficulties that the functions t,hk(") (at,h,,,(")/an)
and $?)(&(")/an) are even and the functions $k(")(a+,,,(')/a?8) are odd, both
with respect to 8.
Correspondingly, the functions +m(")+k("), $,,,(")$k(") are even, and the
functions t,h,,,(")$p)
are odd with respect to 8. It follows from here that
cly = 0 and PLi = 0; thus for each n the set (2.64) decomposes into two
independent sets
N N
(2.66) 2 ptiak(*) -A 2 &ak(nj = 0, m =0,1,. .. ,
k=O k=O
N N
[ K , , , W-
- K,,,(") cos (mmocos e) sinhcosh ~ (cos
~ O)]
+ yo)
~,,,(")(h
cos e
=I
4mn
+
sin [(m- k ) n sin €lo] sin [(m k)nsin O0]
for k # m,
2(m - k ) n -k 2(m k ) n +
136 N. N. MOISEEV AND A. A. PETROV
2k +1
2 ny sin ~ 2 ZYdY
sin 8.
+
sin [(2m 1)nsindo]
+ for k= m
=l
2(2m 1)n
FIG.27.
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 137
TABLE4
-\. 1 3 1 3
-
1
R
FIG. 28.
138 N. N. MOISEEV AND A. A . PETROV
The solutions of the set (2.66) describe the modes of oscillations that
are symmetric with respect to the plane y = 0 ;those of the set (2.67) describe
the modes asymmetric with respect to the same plane.
These sets were solved on an electronic computer. The results of the
computations are shown in Figs. 27 and 28. Fig. 27 shows the graphs of the
dependence of the dimensionless natural frequencies A = aaR/g on the liquid
depth for n = 0 (solid curves 1) and for n = 1, l / R = 1 (solid curves 2).
The approximate calculation method of the natural oscillations is des-
cribed in Section 1.3. The results of the calculations of eigenvalues by
formula (1.37) are shown in Fig. 27 (dotted curves). The typical values of
the coefficients of the series (2.63), obtained by the escalator method on a
computer and by formula (1.41) for various values of the variables, are
given in Table 4. From the comparison of all these results one can see that
the approximate method can be successfully applied in estimating the
natural oscillations of liquid in containers.
The dependence of natural frequencies on the relative length 1/R of the
container for n = 1 and h = - 0.8 (solid curves) and h = 0.2 (dotted
curves) is shown in Fig. 28.
Let the axis of the cylinder be parallel to the vector g and introduce
the coordinate system O X Y Z the plane X Y of which coincides with the plane
of the free surface of the liquid, the axis Z being opposite to the vector g.
FIG. 29.
Denote the cylinder radius by R, the liquid depth by H,the cover radius
by H,, and the free-surface radius by L (see Fig. 29). The dimensionless
variables are
x=-
X y = -Y z
R’ R’ R
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 139
a@ -
_ -0 for
(2.70) Y = 1,
a?-
(2.71) -a@
=o on the surface 2.
an
Here, the wetted part of the bottom and cover surfaces of the cylinder is
denoted by 27, and the dimensionless natural frequency by 1 = a2R/g.
The boundary problem (2.68)-(2.71) has the equivalent formulation
which requires the determination of the minimum of the functional
(2.72)
I
F(@) = (V@)%*- 1 @ V S .
1 S
I
Using again the Ritz method, let us use the method described in Sec-
tion 1.2 for the construction of the coordinate functions. Circumscribe
our cylinder with a flat-bottomed circular cylinder with radius R and
height H (Fig. 29). Take as coordinate functions the set which describes
the modes of oscillations in that cylinder. In Section 1.2 this set was
given as
(2.73)
where Js(~,,,(n)~)are Bessel functions of the first kind, the numbers K,,,(") are
the roots of the equation J,,'(K) = 0, and It = H / R the dimensionless depth
of the liquid.
140 N. N. MOISEEV AND A . A. PETROV
Similarly to what was done in Sections 11.2 and 11.3 one can show that
the spectrum of the problem under consideration has double multiplicity
(except for the case n = 0) so that two eigenfunctions, one with the factor
cos no, the other sin nf3 as factor, correspond to each eigenvalue. Thus the
solution of the variational problem (2.72) can be found for instance in the
form
0 1
where
1 = L/R, ro = R,/R.
If the liquid level is inside the lower spherical surface, then the region
under discussion should be put inside the cylinder of height H and of radius
k = L = //2RoH - H2. In the transition formula to the dimensionless
variables and also in the expressions for the function $2) and the coefficients
pz and q2, the value R should be replaced by 8.
After determination of the eigenvalues, they must be multiplied by
RIIBR,,H - H 2 .
The set of equations (1.13) was solved on an electronic computer
by the escalator method. The maximum number of the terms N of the sum
(2.74) has been taken equal to seven. We do not intend to give the graphs
of the dependence of the eigenvalues A on the number N , because they are
of the same kind as the corresponding graphs given in the previous sections.
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 141
TABLE6
1 2 3 4 5 6 7
10
6
Xi
5
0 1 2 3 4 5 6
h
FIG.30.
142 N. N. MOISEEV AND A. A. PETROV
We shall only note that, as was to be expected, the convergence of the process
became worse when the liquid level was either in the higher or in the lower
11
10
6
Xi
5
3
t
FIG. 31.
of the two spherical caps. However, even in these cases, for calculations
of the first three eigenvalues with a relative error less than 1%, it was suf-
ficient to take only seven terms in the sum (2.74). If the liquid level was
inside the straight cylindrical portion of the container, then for obtaining
the same accuracy three terms of the sum (2.74) were sufficient.
FREE L I Q U I D OSCILLATIONS, MOTIONLESS CONTAINER 143
FIG.33.
perpendicular to the torus axis, by R,, the inner and outer radii of the
free surface of the liquid by L, and L,, respectively. Introduce the coordinate
system O X Y Z whose axis Z coincides with the torus axis and is opposite
to the vector g, and the plane XY coincides with the plane of the free surface.
We shall measure the liquid depth H as it is shown in Fig. 32.
Introduce the dimensionless variables
X Y z=- z
x = - -R
,O y=q’ RO
and the cylindrical coordinates r, 8, z by the formulae
x=rcos8, y=zsin8, z=z.
(2.76) !
F(@)= (P@)'dt - iZ @'dS,
s I
to which we again shall apply the Ritz method. Let us construct the system of
coordinate functions by the procedure described in Section 1.2. Circumscribe
the region t with annular cylinders bounded by cylindrical surfaces of radii
R, and cR, = R, + 2 R , and by planes Z = 0 and Z = - H. As coordinate
functions we shall take the set of functions which represent the natural
oscillations in the region between two coaxial circular cylinders of height H.
The solution of this problem is given in Section 1.2 and has the form
It is not difficult to show that the spectrum of this problem is twice degenerate
for rt # 0, i.e., to each eigenvalue correspond two functions one of which
includes the factor cos no and the other sin d3. Thus the solution of the
variational problem (2.75) can be sought for in the form
FREE LIQUID OSCILLATIONS,MOTIONLESS CONTAINER 146
FIG.33. Key: Solid line ___ : a = 0.6; dashed line - - -: a = 0.4; open dots
represent experimental points when o! = 0.6 ; solid dots ( 0 )represent experimental
(0)
points when a = 0.4.
culated by (1.18) and (1.19). If in the case under discussion It > a, where
a = R/R,, then
146 N. N. MOISEEV AND A. A. PETROV
[$,(n)t,~)], = Ndr,
0 1,
n = l n = 2
h -\Y 1 2 1 2
If the liquid fills less than half of the torus then it is more convenient
to describe the region by cylindrical surfaces of radii R, = R, + R -
V2RH - H2 and c" R, = R, + +
R V2RH - H2. In this case 1, = 1 and
1, = c". In the transition to dimensionless coordinates in the expressions
for the functions $(,") and the coefficients py
and q t i the values R, and c
should be substituted for a,, and E. The eigenvalues obtained from the
solution of the set (1.13) should finally be multiplied then by Rol Po.
0.6
0.5
0.4
'4 0.3
0.2
0.1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
h
FIG.34
The system (1.13) has been solved by the escalator method on an elec-
tronic computer. It turned out that a t any values of the problem parameters
it was enough to take eight terms in the sum (2.77) in order to calculate
the first two natural frequencies with a relative error less than 1 per cent.
Figure 33 shows the dependence of the first two dimensionless natural
frequencies li = ai2R/g on the liquid depth h in two containers, the first
with a = 0.6 (solid curves) and the second with a = 0.4 (dotted curves).
This figure also contains the results of experimental measurements of
these frequencies for the container with a=0.6 (open circles) and for the
container with a = 0.4 (full circles).
In Table 6 the first two coefficients of the sum (2.77) for the case a = 0.4
are given.
The first natural frequency was also computed by the approximate
formulae (1.37) and (1.38). Figure 34 shows the dependence of A, on It,
computed by the escalator method on a computer (solid curve) and by the
approximate formula (1.37) (dotted curve). In this case also the ap-
proximate formula gives satisfactory results.
148 N. N. MOISEEV AND A. A. PETROV
APPENDIX
The application of the Ritz method reduces the problem of liquid os-
cillations in a container to the solution of the linear homogeneous equations
(1.13). If we introduce the matrices
p= ~pm\l:m=l# Q = IlqnmlIL1
and the vector x(a,,.. .,a,) then the system (1.13) can be replaced by the
matrix equation
(A.1) Px = AQx.
Since the problem is self-adjoint and kinetic and potential energies are
positive definite, the matrices P and Q are symmetrical and positive definite.
In the process of designing a scheme of calculation for the natural fre-
quencies and modes of free oscillations there have been tested various methods
to construct the solutions of the system (A.1). A. A. Petrov, Yu. V.
Pukhnachev and Yu. P. Popov have developed a special form of the escalator
method which makes it possible to improve sequentially the unknown
values. For the simplest case the matrix Q is unitary and the roots
of the equation (A.l) are simple this method has been discussed in various
textbooks (see, for example, [12]).
The modification of the escalator method for the case where Q is an
arbitrary, symmetrical, positive definite matrix is discussed below.*
The idea of the escalator method is as follows. Consider the equation
(A4 PRX = A Q R X .
The matrices P R and QR are composed of the R first elements in the R
first rows of the matrices P and Q. Assume that the solution of this equa-
tion is known. Denote its eigenvalues by &(R) and the eigenvectors corre-
sponding to these eigenvalues by x , ( ~ ) .
Further, compose the matrix P R + l by adding to the matrix P R the R
+
first elements of row and column number R 1 and the element $ R + l , R =a
located on the diagonal. Let us call this procedure “bordering of the matrix
PR.” Denote by p the row (column) of R elements which we have added
to the matrix. Similarly, border the matrix Q R by the row (column) q and
the element p located on the diagonal.
Consider the equation
(A.3) p R + l x = AQR+,l_x*
with
A , c R ) ( A ) = (p - Aq,xP),
Proceeding from the minimax properties of the eigenvalues (see, for in-
stance, Chapter I, $ 4 of [13]) it is easy to show that the eigenvalues of the
equations (A.2) and (A.3) are intertwined as follows:
(A.lO) 3y(R+1) < < &(R+l) < &(R) < . . . < 1 ~ \<( ~Ak",i".
)
160 N. N. MOISEEV A N D A. A. PETROV
I t is evident that the smaller the value (y1, the closer to A!? lies one of the
zeros of the escalator function A(R + l). The approximate relationship
2;'- p + u = - Y
can be improved:
this follows from the inequalities (A.lO). The inequalities &(R+l) <
&(R)
<
and k(R) 2:;') show that the multiplicity of this root may reach v
and, finally, Y + 1 if the given inequalities are replaced by equalities. Thus
in the transition from (A.2) to (A.3) the root multiplicity may either de-
crease by unity, remain unchanged, or increase by unity. N o other
possibilities exist.
Gmsider in more details the case when the root multiplicity decreases
by unity. To the eigenvalue I(R+l) of multiplicity Y - 1 we should put
in correspondence P - 1 eigenvectors. Suppose that all AI'R)(AI'R)) # 0,
FREE LIQUID OSCILLATIONS, MOTIONLESS CONTAINER 161
.
Y = 1,2,. ..v. We determine x,('+l) in the form of a linear combination
of the vectors xJR)with y = 0, s = 1,2,. , .,r + 1:
(p - A1(R)Q,X,(R+1)) = 0.
These conditions give us the set of equations
(QR+lX,(R+",X1(R+l)) =1 + tlltrl = 0,
( Q ~ + i d ~ + ' ) , ~ g ' ~ +=
~ ) 1) + t&i + t&,2 = 0,
.......................................
(p - A p q , x , c R + ' ) ) = A J R ) ( A + R ) ) + t,lA2(R)(Al(R)) + . . . + t,,A;yl(Ap) = 0.
Let now some of the AJR)(AJR)) be equal to zero. Then the corresponding
X,(~+I) is assumed to be equal to x i R ) , and the rest is constructed by the
procedure discussed above. But applying this procedure in the case where
all A J R ) ( A i R ) ) are equal to zero, we would be able to find for an (v - 1)-
fold eigenvalue v instead of v - 1 linearly independent eigenvectors, which
is impossible. However, one can prove the following statement : The vanishing
of all A1[R)(A,,,(R)), .
m = 1,2,. . ,v is achieved then and only then, when the
multiplicity of the root A;R+l) of the equation (A.3) is equal to or greater
than the multiplicity of the root &(R) of the equation (A.2).
The proof of this statement is rather complicated, and we do not present it
here. Thus, if the root multiplicity reqains, to each of the roots A,,,(R+l)=
m = 1.2,. . ..v we can put in coeppondence the vector x , ( ~ + ~ ) assuming
x J R + l ) = xJR), y = 0. In solving the escalator equation we may obtain
one more root equal to jzl(R) and'from formulae of the kind (A.S), (A.6)
the eigenvector corresponding to it. The root multiplicity, consequently,
will increase to v + 1.
Systematizing what was said ahqve one may explain the following
procedure of solving the escalator equation E R + l ( A ) = 0. We calculate
first
and on the basis of their closeness (or equality) to zero we decide on the
closeness (or coincidence, respectively) of zeroes of the escalator function and
162 N. N. MOISEEV AND A. A. PETROV
+
its poles, Partition, further, the semi-infiniteinterval [0,m) into R 1intervals
.
by the points {A$R),. . , i l R ( R ) } . Assume that one of the points at which
1yrJ is small (or equal to zero) lies between points at which (y,1 is, at least,
of the order of unity. Between these two points there are two eigenvalues
of the equation (A.3). One of them is close to ilf) and has already been
taken into account. If Iy,,J = 0, then, in order to find the other root we take the
whole interval [A!:!- l , A ! ~ ~ l ] and, simultaneously, from the last sum of the
formula (A.7), which determines the escalator function, we omit the fraction
with the index Y,. In case of coincidence of the root determined with z?
we notice an increase of the root multiplicity. If lyr,l # 0 then, for the
determination of the second root out of the two intervals [A!:i,,#:)] and
[A!P’,il!:il], we take either one or the other depending on whether the value
FRfl@)) is positive or negative. Similar reasoning can be camed out
in the case when between the points at which 1y.l is of the order of unity
there are some points at which lyrl is small or equal to zero.
We use the same procedure if the root multiplicity Y> 1, examining
the expression
CONCLUSION
References
12. FADDEEV,D. K., and FADDEEVA, V. N., “Computational Methods of Linear Algebra.”
Fizmatgiz, 1960. (An English translation is available published by Freeman,
San Francisco, 1963).
13. COURANT,R., and HILBERT, D., “Methods of Mathematical Physics.” Interscience,
New York. 1963.
Abbreviations:
BY R. M . ROSENBERG
.
Department of Mechanical Engineering. Division of Applied Mechanics. University of
California. Berkeley Califwnia
Page
I. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
I1 . Admissible Systems and Motions . . . . . . . . . . . . . . . . . . . . 168
1. The General System . . . . . . . . . . . . . . . . . . . . . . . . 158
2. Admissible Systems . . . . . . . . . . . . . . . . . . . . . . . . 159
3. Admissible Motions . . . . . . . . . . . . . . . . . . . . . . . . . 162
I11 .
The Trajectories in Configuration Space . . . . . . . . . . . . . . . . 163
1 . Transformation and Trajectories . . . . . . . . . . . . . . . . . . . 163
2. The Restricted Principle of Least Action . . . . . . . . . . . . . . . 166
IV . General Description of the Geometrical Method . . . . . . . . . . . . 167
V. Trajectories of Admissible Motion . . . . . . . . . . . . . . . . . . 169
1. General Properties . . . . . . . . . . . . . . . . . . . . . . . . . 169
2 . Admissible Trajectories of Autonomous Systems . . . . . . . . . . . . 170
VI . Normal-Mode Vibrations of Nonlinear Systems . . . . . . . . . . . . 173
1. A New Definition of Normal Modes . . . . . . . . . . . . . . . . . 173
2. Straight Modal Lines . . . . . . . . . . . . . . . . . . . . . . . . 173
3. Interpretation of Modal Lines . . . . . . . . . . . . . . . . . . . . 176
VII . Properties of Motions with a Rest Point . . . . . . . . . . . . . . . 176
1. The Transversals (or P-curves) . . . . . . . . . . . . . . . . . . . 176
.
2 The 2'-curves . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
VIII . Special Autonomous Systems . . . . . . . . . . . . . . . . . . . . 187
1. Smooth Systems . . . . . . . . . . . . . . . . . . . . . . . . . . 187
2. Uniform Systems . . . . . . . . . . . . . . . . . . . . . . . . . . 188
3. Sequential Anchored Systems . . . . . . . . . . . . . . . . . . . . 189
4. Homogeneous Systems . . . . . . . . . . . . . . . . . . . . . . . 190
6. Symmetric Systems . . . . . . . . . . . . . . . . . . . . . . . . 190
IX . Theateb-Functions . . . . . . . . . . . . . . . . . . . . . . . . . 191
1. Their Origin and Importance . . . . . . . . . . . . . . . . . . . . 191
2. The Inversions . . . . . . . . . . . . . . . . . . . . . . . . . . . 194
3. The Amp-Functions . . . . . . . . . . . . . . . . . . . . . . . . 196
4. The Sam- and Cam-Functions . . . . . . . . . . . . . . . . . . . . 197
X. Nonsimilar Normal-Mode Vibrations . . . . . . . . . . . . . . . . . 199
1. The Perturbation Potential . . . . . . . . . . . . . . . . . . . . . 199
2 . The Curved Modal Line . . . . . . . . . . . . . . . . . . . . . . . 200
3. Integrable Cases . . . . . . . . . . . . . . . . . . . . . . . . . . 201
4 . GeneralProperties . . . . . . . . . . . . . . . . . . . . . . . . . 202
XI . Exact Solutions to Steady-State Forced Vibrations . . . . . . . . . . . 203
1. Introductory Remarks . . . . . . . . . . . . . . . . . . . . . . . 203
2. Homogeneous Systems . . . . . . . . . . . . . . . . . . . . . . . 204
3. The Response Curves . . . . . . . . . . . . . . . . . . . . . . . . 210
166
166 R. M. ROSENBERG
I. INTRODUCTION
AND MOTIONS
11. ADMISSIBLESYSTEMS
where O,(w) is the negative of the potential energy V, stored in the spring
Sij terminating on the ith and jth mass points, w is the length change of
Si, and Fii is the force which has produced the length change. Evidently,
combining (2.1) and (2.2), one has
(2.3) F-(m)
$1 = - F;i(- 20).
Each mass point is connected to one, several, or all other mass points
+
by a spring; hence, the maximum number of springs is n(n 1)/2. If one
extremity of a spring terminates on an infinitely large mass point it is called
an anchor spring; all others are coq5ling springs. I t is not essential that all
+
n(n 1)/2 springs be present in the system. However, it is required that,
in the absence of additional constraints, no mass point can be displaced
without giving rise to an elastic force acting on every mass point. For in-
stance if
(2.4) u*=o, ( i = O ,... I - l ) ,
)
the rth mass would be decoupled from all its neighbors on the left, and the
+
system of n 1 mass points would be separated into two spring-mass
systems with no elastic coupling between them. Such a separation into
several mutually uncoupled systems is excluded. In consequence the mini-
mum number of springs is n.
Each mass point has a single degree of translational freedom ui in the
direction of the chain. In view of mo = 00, one has uo E 0; hence, the
system cannot accelerate in the rigid body mode.
Each mass point may be acted upon by a time-dependent force, called
an exciting force and denoted by fi(t). Exciting forces act in the direction
of the chain of mass points. It is assumed that the spring forces Fii and the
exciting forces f, are the only forces acting on the mass points.
2. Admissible Systems
In order to fit into the framework of vibrations, linear or nonlinear,
systems must have certain properties whose physical description is :
(i) in the absence of exciting forces, the system must possess a single
equilibrium configuration, and it can execute “free vibrations” about that
equilibrium configuration ;
160 R. M. ROSENBERG
(ii) in the presence of exciting forces, the system can execute so-called
“steadystate forced vibrations.”
Some of these conditions can be insured by endowing the potential
function of the entire system with certain properties. This potential
function is, evidently,
n-1 I(
O(0,.. .,O) = 0,
a .
-O(0,. .,O) = . .. = -a O(0,. . *,O) = 0.
8% aun
Equations (2.6) and (2.7) insure that there exists a unique equilibrium
configuration, and the origin of the ui has been so chosen that every ui
vanishes in the equilibrium configuration. Equation (2.8) is a direct con-
sequence of the symmetry properties of the individual springs, described
in (2.3).
A potential function 0 which satisfies all of the above properties will
be called admissible. Moreover, when 0 is admissible, the azctonomous system
i s admissible.
NONLINEAR VIBRATIONS 161
(2.10)
where T is a constant. Equation (2.9) says that the exciting forces are
periodic and, if there is more than one such force, they all have the same
period T.
Exciting forces /<(t)which satisfy the above conditions are called ad-
missible. Moreover, when both, 0 and the f,(t) are admissible, the nonauto-
nomous system is said to be admissible.
Among admissible systems, autonomous or not, we distinguish between
those that are linearizable, nonlinearizable, and weakly, or strongly non-
linear. In our system, these properties are associated with the spring forces;
hence, it may be best to define them in terms of these forces; i.e., in terms
of the partial derivatives of 0.
Because of the admissibility conditions on 0, these derivatives may
be written in the form
,. n
(2.11)
where the functions Pi vanish at the origin and are of degree higher than 1
and, because of (2.8),
(2.12)
-
aij = aii.
3. Admissible Motions
From the definition of admissible systems, it is evident that these may
be strongly nonlinear, or even nonlinearizable. In consequence, it is not
likely that one can find the general solutions of the equations of motion or,
the class of all motions of which the system is capable. However, it is
possible to deduce certain types of motion and, in particular, the motions
which are the most important from the point of view of technical application.
These are the motions which, in linear systems, are called the “natural free
vibratiqns in normal modes” if the system is autommo%s, and the “steady-
state forced vibrations” if the system is nmutonomous. The latter include
the phenomena of “resonance” and of the tuned dynamic vibration absorber.
Both, the free vibrations in normal modes, and the steady-state forced
vibrations are so-called “vibrations-in-unison”. An admissible system is
said to execute a vibration-in-unison if its motion satisfies all of the fol-
lowing properties [7] :
(i) There exists a T = corzstant such that
(2.15) ui(t) = z ~ i ( b + T), (i = 1,. . .,%).
In words, the motions of all masses are equiperiodic.
(ii) If t, is any instant of time, there exists a single to in
In words, during any time interval of a half period, the velocities of all
masses vanish precisely once or, all masses attain their maximum displace-
ment (in absolute value) precisely once during any interval of a half
period.
(iv) Let r (fixed) be any one of the i = 1,. . . ,n. Then every %&) and
ur(t) may be written, for a.ll t, in the form
(2.18)
NONLINEAR VIBRATIONS 163
where
systems, there aIways exists some domain [ x j , X j ] ,(i = 1,. . .,n),of the
configuration space in which every xi@) has an inverse
(3.8) t = ti( q) I
Hence, summing the equations (3.7) over N, and integrating them once over
the domain [xi,Xi] one finds
"
(3.9)
1 "
-2
2 i-1
=
i<* . .JJ + 12 + i=
u(x~, 21Fi(xjJX<)
where h is a constant of integration. and
(3.10)
i
F ~ ( x < , X=~ )fi(tj(xi))dxi.
xi
It is evident that (3.9) is an equation of the energy balance of the system,
and the quantities Fi(xi,XI),(i = 1,. . .,rc) represent the work done on the
unit mass p i n t by the force components fi(t) in a displacement [xi,Xi].
When these forces vanish identically, (3.9) reduces to the energy integral
of the autonomous system. The idea of utilizing the inversions (3.8) is also
employed in Rauscher's method [3].
The substitution of the first of (3.6) in (3.9) results in an equation in
xla which may be explicitly solved for that quantity. That equation together
with the system H and the second of (3.6) is sufficient to eliminate the
parameter t. Hence, the resulting equation is that of the trajectory of the
unit mass point. It turns out to be
r 1
(3.11)
where we have omitted the self-evident subscripts on the 6. The corresponding
equation for the autonomous system is
(3.12)
The equation (3.12), but restricted to two degrees of freedom, was
probably first given by Kauderer [lo]. An equation, equivalent to (3.12)
has been given earlier [ll];in it, the arc length s was used as a parameter.
166 R. M. ROSENBERG
The elimination of s which has resulted in the compact forms (3.11) and
(3.12) was first given by Kinney [12].
The systems of equations (3.11) and (3.12) are denoted as the systems M
because they are derivable from a principle akin to the principle of Mau-
pertuis, given below. Similarly, the notation H in (3.3) may recall Hamilton's
principle.
(3.14)
'I
and
(3.16)
is the square of the line element in the configuration space. Hence, equa-
tions (3.12) follow directly from the principle of Maupertuis.
Now, it is well known that the principle of Maupertuis is only applicable
to scleronomous conservative systems. However, if one defines an action
$1
one can show that (3.11) is immediately obtained from the variational
principle
(3.17) SA=O
and from (3.10). Thus one has the [13].
Restricted Princi~lcof Least Action: The natural motion of any system
.
with equations of motion (3.3), 0% a%y interval [xi,Xi], (i = 1,. .#n)in which
the imersions (3.8) exist, i s such that the action A, defined in (3.16), i s
a mircimzlm.
NONLINEAR VIBRATIONS 167
While it has only been shown here that the action is stationary, the
usual proof that the stationary value is minimum also succeeds here.
The above principle is restricted to the domain [xi,Xi] where the in-
versions (3.8) exist. In problems of dynamics in general, this domain is
not known u priori. However, because of (2.16), (2.17) and (2.18), that
domain is, for admissible motions, sufficiently large to contain at least the
closed intervals lying between the maximum amplitudes and the equilibrium
position, and the corresponding time interval is sufficiently large to cover
at least one quarter period. Since the systems M are used directly in the
methods to be presented, the restricted principle of least action is of con-
siderable usefulness in problems where admissible motions of admissible
systems are examined.
(44 u, = 2
i- 1
uijxj, (a = 1,. . ..n)
where the aii are either constants or functions oft. Under it, equations (3.3)be-
come linear while (3.11) and (3.12) remain nonlinear. The intrinsic nonlinearity
of the system M will perhaps be better appreciated if on recalls that its
integrals are Lissajous curves. These may have many points of self-crossing
and may even be compact on some domain of the configuration space. One
may conclude that admissible systems with linear springs are more easily
examined by means of the system H of equations of motion than by M.
The fortunate circumstance that the superposition theorem holds in this
case can, then, always be utilized to find the general motions and, if one
wants to work very hard, their trajectories in the configuration space [XI.
168 R. M. ROSENBERG
(4.4)
.. = UX,(Xl*%n(&
x1 * * * %&l)) + fl(&
an equation in x1 only.
If the system is autonomous, fl(t) 3 0, and (4.4) is of the form
(4.6) 21 = g(x1).
Then, for initial conditions (for instance)
(44 x1(0) = x,, ZJO) = 0 ,
and with the definition
(4.7)
V. TRAJECTORIES MOTION
OF ADMISSIBLE
1. General Properties
(5.3) xj' +
[ux, f i ( t ( 4 ) I = uxj + f j ( t ( x j ) ) s (I' = 2,. - 1%)
because
T(21,.. .,A?,,)= - 2
1
Xi'
2 i=l
and this is also the maximum equipotential surface because (6.8) coincides
with (6.6) when the kinetic energy vanishes. In order to call attention to
the fact that (5.8) defines the bounding surface when the system is auto-
nomous, it will be denoted as the La-surface.
From (2.5), (2.6), (2.7), and (5.8) it follows that:
(i) the La-surface is a closed, smooth surface surrounding the origin
of the configuration space, and star-shaped with respect to it [9].
From (2.8) one sees that
(ii) the La-surface is symmetric with respect to the origin of the configura-
tion space.
From the admissibility condition (2.6) and (6.8) it is evident that
+
(iii) the quantity U h > 0 a t the origin of the configuration space.
This quantity does not change sign in the open, finite domain D, bounded
by the La-surface, it vanishes on L,, and it is negative (at least near the
La-surface) in the infinite open domain bounded by the La-surface.
These properties lead immediately to
Theorem V-I1 : Every trajectory, admissible or not, defined by solutions
x I . ( x I ) ,(i = 2 , . . .,n) of (3.12), lies i n the closed, finite domain D, of the con-
figuration space, bounded by the La-surface defined i n (5.8).
It is not necessary that all trajectories actually attain the La-surJace.
However, if they do, this occurs necessarily at an instant t, when all a,(t,)
vanish. Trajectories which attain the L,-surface will be called trajectories
of motion with a rest point [19]. Not all motions with a rest point need be
admissible because the corresponding trajectories need not pass through
the origin. However, all admissible motions are motions with a rest point.
Next we state a property, first observed by Mawhin [6]. I t is
Theorem V-I11 : If 8 i s any trajectory, admissible or not, of an admissible,
autonomous system, the trajectory 8 ,symmetric to 8 with respect to the origiB,
is also a trajectory.
The proof follows immediately from (2.8) and its consequence
Finally, we have
(6.11)
But equations (6.11) are precisely the conditions for orthogonal intersec-
tion.
Special additional properties hold for similar vibrations-in-unison of
autonomous systems. They are given in [9].
Theorem V-VII: If there exists a straight line which intersects every equi-
potential szcrface of an admissible, autonomous system orthogonally, it i s a
trajectory.
(6.12) . +
U(X1,. .,xn) h* = 0. (O< h* < h)
and, in view of the properties of admissible potential functions, the equi-
potential surfaces are simple, closed, non-intersecting surfaces surrounding
the origin, and compact on D,.
It follows from (6.2) that, in the case of similar motion,
VI. NORMAL-MODE
VIBRATIONS SYSTEMS
OF NONLINEAR
(6.1) xi = rf& .. ..
(i = 1,. ,n)
and, specifically, they are [9]
j- 1
cos ej, e,, = n/2, (i= 1,. ..,n).
or in view of (6.7)
(6.9) 0, = 0, (i= 1,. . .,n - 1).
NONLINEAR VIBRATIONS 176
i
Y
X
- y=ct ; x
Suppose the system is linear. Then, the equipotential surfaces are given
by the concentric three-dimensional ellipsoids (with center at the origin)
(6.11) t(x2 + y2 + %a! = h*, o< h* < h,
and the vibration in the pth normal mode is given by
(6.12) x = Xp cos apt, JJ= Yp cos opt, z = 2, cos apt
where X,, Y,,Z p are constants.
It follows from these equations that the pth normal mode satisfies
Then, the modal line in the pth mode is the intersection of the planes
(6.14) P
y = cy+x, z = l$X.
(a1 (b)
FIG. 3. Modal line of nonlinear system with three degrees of freedom.
axes of the ellipsoid. The construction of a modal line, and its location
inside the bounding ellipsoid is shown in Figures 2(a) and (b).
If the problem is nonlinear, but the normal mode vibrations are similar,
equations (6.14) still apply, but, in that case, the bounding ovaloid is no
longer an ellipsoid.
If the problem is nonlinear] and the normal-mode vibrations are non-
similar, equations (6.14) are replaced by equations of the form
VII. PROPERTIES
OF MOTIONS WITH A REST POINT[19, 241
Motions with a rest point have been defined as those in which there
exists a time t, at which the velocity of all mass points of the physical system
vanishes. If the system is autonomous, the trajectories corresponding to
motions with a rest point are those which attain the Lasurface, and ad-
NONLINEAR VIBRATIONS 177
a
-O(U,V) = - auk - 6 V k + P(z4.V).
au
(7.2)
-o
a ( ~=, -~ ) tUk - avk + Q(.,.)
av
where k is a positive, odd integer, and P and Q are of degree higher than k.
The determinant
(7.3)
.. -
m2v =
a O(u,v),
av
178 R. M. ROSENBERG
(74
(7.7) 7-(2,$)- V ( x , y )= h.
The bounding surface (here a curve) is the locus of the points P(X,Y) for
which the kinetic energy vanishes when the energy level of the motion is k,
and the equation of the La-curve is
(7.8) + h = 0.
U(X,Y)
The equipotential curves are the functions
(7.9) Y = &)
which satisfy the equation
(7.10) +
U(x,qJ) k* = 0, OC h' <h
and they will be called E-crcrves. By differentiating (7.10), one finds that
E-curves satisfy the differential equation
(7.11)
(7.12) Y =e(4
which satisfy the differential equation
(7.13)
We observe that (7.13) is singular only at the origin, and the equation
of the trajectories (7.6) is singular only on the La-curve. Consequently, we
define as regular p o d s of the xy-plane all points, except the origin, which
lie in the open, finite domain D, bounded by the La-curve.
NONLINEAR VIBRATIONS 179
Theorem VII-I1 : The origin i s always a node for the P-curves of admissible,
autonomous systems having two degrees of freedom.
(7.14)
u, = XU,(O,O) + y ~ x y ( 0 , O+) . -
* I
u, = X V , y ( O , O ) + yU,,(O,O) + . . .*
(7.16) *=
dx
- xu, - 'pu,
XU, + pU,, i-
+ . ..
. ..
where the partial derivatives are evaluated at the origin. But, the integrals
of (7,11), and hence (7.16), are known; they are E-curves and they form a
continuum of simple, closed, nonintersecting curves surrounding the origin.
Then, it follows that the origin is a center for the E-curves and this requires,
by the Poincart! theory of singular points [l], that
(7.17)
(7.18)
These are two of the three conditions which are necessary for the origin to
be a node. But, in view of (7.16).
This is the third necessary condition, and (7.18) and (7.19) are also sufficient.
Hence, Theorem VII-I1 is satisfied in the lilrearizable case.
180 R. M. ROSENBERG
(7.20)
(7.21)
are introduced in the resulting equation. Then, one finds, instead of (7.17)
in the neighborhood of the origin
(7.22) _
dr] - - cvk + b(E - q ) k + . . .
at - - atk - b ( t - q)b + .. . *
The singularity at the origin of this equation may be discussed by the method
of Argemi and Sideriadks [22, 231. These authors consider the equation
(7.23)
(7.26)
where k(A) = g(1) - it/(& Now, this last equation has simple singular
points on the it-axis at the zeros Ri of k(it), and their discussion is accessible
to the Poincarh theory. Expansion of the numerator and denominator
of the right-hand side of (7.26) near the singular points gives
NONLINEAR VIBRATIONS 181
_ - h’(li)l + . . .
dl -
(7.27)
dx f(&)x + ,. ’
,
The authors show that only saddles and nodes can occur, and one has
a node if f(li)h’(&)> 0,
(7.28)
a saddle if f(A,)h’(l$)< 0.
The quantity t(&)h’(&)does not vanish, in general, because the zeros of
/ ( A ) do not coincide with those of h(1) when X # Y,and h’(1) does not
vanish at the zeros of h(1). Moreover, when the 1, are ordered, according
(a 1 (b)
FIG.4. Saddles and nodes in xl-plane and mapping on configuration space.
(7.29)
/(A) = - b[(2 - 8) + (1 - 4 s 1 v
h(A) = - b[A4 - aAs + P A - 13
where
C
a=2--
b'
(7.30)
It is now necessary to determine first the zeros & of h(l) and, next, ,the
sign of the product fh' at these zeros. Writing h(A) = 0 in the form
(7.31) A4 - 1 = (als - 8)l
one sees that the left-hand side is even in A and has two real zeros at 1 = f 1,
and the right-hand side is odd and has one or three real zeros. Hence, these
two functions have always at least two real intersections, or h(A) has always
at least two real zeros. One may assume without loss of generality*
that c >, a. Now, it follows from the second of (7.29) and from the assump-
tions on a, b, and c, that
h ( w ) = h(- 00) = - 00,
h(0) = b,
(7.32)
h(1) = - b(- a + 8) = - (C - a ) < 0.
h(- 1) = - h(1) = c - a 2 0.
* If, on the contrary G < a, we simply denote the displacement of mII by u and
that of m1 by u which restores the assumed relation.
NONLINEAR VIBRATIONS 183
one root, say 4, is positive and lies in the interval O < A,< 1. and & is
negative. Moreover, the slopes are
(7.33)
(7.34)
Since b > 0, a/b > 0, and 0 < 4 < 1, and since 4 < 0, one has from (7.36)
f(4)< 01
(7.36)
f(U< 0.
Combining (7.33) and (7.36),
(7.37)
184 R. M. ROSENBERG
Hence, the xl-plane has, at least, two singular points on the A-axis. One
lies on the positive branch between 0 and 1, and it is a node. The other
lies on the negative branch, and it is a saddle.
We shall assume that Al,z are the only zergs* of h ( l ) . Under this assump-
tion we have now demonstrated
Theorem VII-11: When k = 1.3 and h ( l ) has only two zeros, all trans-
versals of admissible, autonomous systems pass through the origin of the xy-
plane with common slope OOp.
To determine the properties of P-curves not near the origin, we determine
the locus of their inflection points, if any. By differentiating (7.13), one has
It follows from (7.39) that F-curves pass through, and are symmetric with
respect to, the origin.
Next we determine the locus of points such that the tangents to P-
curves at these points pass through the origin. Evidently, these are the
points for which
(7.40)
It can happen that k(L) has four real roots when R = 3. This has been shown
to be an exceptional case and is not treated here. [MI.
NONLINEAR VIBRATIONS 186
from the origin, that quadrant contains at least one G-curve. Consider
now the transversal issuing from P I . I t points initially towards the origin.
If Po and P I do not coincide, this particular transversal must have an
inflection point prior to its amval at the origin with slope 6;. In that
case, the transversal issuing from Po points initially above the origin, and
it must also have a point of inflection before arriving with slope 6; at the
origin. Hence, in general, there exist transversals with points of inflection,
and there exists at least one F-curve in the first quadrant. The transversals
will be endowed with a sense of direction by considering the points P ( X , Y )
which compose the Ldcurve as their points of issue. These directed trans-
versals in Figure 6 are supplied with arrow heads pointing toward the origin.
In this way, the transversals give not only the direction, but also the sense,
of the forces acting on the unit mass point of the pseudo-system.
2. The T-curves [19]
The trajectories corresponding to motions (at energy level h) of admissible,
autonomous systems fall into two general classes: those which attain the
La-surface, and those which do not. By definition, T-curves are those integral
curves of (3.12) which do intercept the La-surface. Hence, the T-curves
constitute the class of all trajectories which correspond to motion with a
186 R. M. ROSENBERG
(7.42)
(7.43)
1 Uxy
y”(X,Y) = -
3 1
But, by Property 1 and (7.13),
+ y‘(Uyy -
ux
UXZ
1-
- Y’UZY)
x.*
(74 y y x , ~=
) e y x , y ) = u,(x,Y)/u,(x,Y).
Consequently,
VIII. SPECIALAUTONOMOUS
SYSTEMS
1. Smooth Systems
It is evident from the .equations of motion (3.3) of the pseudo-system
that all differences between given, autonomous systems consist in differences
between their potential functions.
An autonomous system is said to be smooth if its potential function is
of the form
188 R. M. ROSENBERG
where the a$) are constants. Then, the spring force of a spring Sij between
the masses mi and mi is
where
One sees from (8.4) that uniform systems are those in which all springs are
equal, and all masses are equal (in which case one may, without loss of
generality, put the masses equal to unity).
For uniform systems, one can readily demonstrate an interesting and
useful property. Let us denote as the potential function of the associated,
linear system that whose potential function is [9]
where the values of a(l) are the same in (8.4) and (8.6). Then, one has
Theorem VIII-I : The modal lines of an autonomous, uniform system are
straight. Their directions coincide with those of the associated linear system.
One proof of this theorem is due to Mawhin [6] who shows that the modal
lines of (8.4) are straight, and that their direction is independent of m.
The importance of this theorem is evident. The modal lines, when U = 0,
are the principal axes of the n-dimensional ellipsoids
and the directions of these axes are readily found by means of linear eigen-
value theory. The time-historyof the motions can, then, be found subsequently
by simple quadratures.
The geometrical interpretation of this result is that the equations
+
where i = i 1.
Haughton [27], using methods similar to those of Mawhin [S]proves the
following remarkable
Theorem VIII-I1 : If an admissible. anchored, seqrcential system has straight
modal lines, their directims are the same as those of the associated, linear
system of potential function
(8.9)
If, in an anchored, sequential system, all masses are equal and all springs
are equal it is said to be anchored, sequential and d i m . In that case,
one can readily show that straight modal lines do, indeed, exist; the prob-
lem of normal mode vibrations can, then, be completely solved by con-
sidering a linear system and by performing one additional quadrature.
190 R. M. ROSENBERG
where k is a real number in O C k < do. I t derives its name from the fact
that the potential function is homogeneous in the xi of degree k 1.+
It is interesting to forni the expression for the spring forces. These
forces are given by
This equation shows that the spring forces are odd functions, and proportional
to the kth power, of the length-changes. Clearly, when one puts k = 1,
one recovers the linear system.
Physically, neither masses nor springs need be equal in homogeneous
systems. Instead, the springs are nonlinear “in the same way”; i.e., all
obey the same simple power law.
It is not difficult to construct a physical system for which k = 3 [21].
Consider an arrangement of masses and springs as in Figure 1. However, in
the present case, the translations ui are not in the direction of the chain,
but t t o d to it. Moreover, every spring is assumed l i w and its free
length is exactly equal to the distance between the masses interconnected
by that spring when the system is in the equilibrium configuration. Clearly,
the spring forces are always restoring, or 0 is negative, definite; moreover,
the spring forces reverse sign when the displacements do, or 0 is symmetric
with respect to the origin. However, infinitesimal displacements I(~ which
are small of the first order produce length-changes in the springs which
are small of the second order. Hence, for sufficiently small displacements,
k = 3.
By applying the criterion (6.3) one can readily prove [S]
Theorem VIII-111: The modal lines of admissible, autolumtous, homo-
g m u s systems are straight. Their direction i s defined by the sysiems of roots
of the transcedental eqwations (6.4).
6. Symmetric Systems
The physical properties of symmetric systems are the following: All masses
are equal (hence, assumed equal to unity, without loss of generality), all
spring forces are smooth (in the above sense), all anchor springs are equal,
and the spring forces in all coupling springs are polynomials of thesame degree.
This is a generalization of the “symmetric two-degree-of-freedom system”
1291 which was called “symmetric” because the system had literally physical
symmetry about its centerline. (It should be noted that Mawhin denotes
the spring between m, and m, as So, that between m, and mi as S j , and
that between m, and m,, as S,.)
If a symmetric system has only two degrees of freedom, one can readily
show [29] that straight modal lines exist, and they have the same inclination
as those of the associated linear system.
(9.4)
pj" = mil-.
We shall show that (9.3) has periodic solutions of period Tp,and it follows
from (9.1) that, then all x&), (i = 2,. , ,,n)are also periodic of the same
period. Hence, the resulting motion is a normal mode vibration.
We shall integrate (9.3) for the Cases I and 11 where, for
The conditions for Case I1 are ilzitial conditions because x, and are speci-
fied at the same instant of time. However, in Case I the velocity and dis-
placement are prescribed at different instants of time. Hence, in Case I,
it must be shown a posteriori that, if x,(O) = 0, then there exists an %,(O)
such that %,(to) = X,, (X,,t, > 0).
First, we change (9.3) into the canonical form by means of the trans-
formations [20]
(9.6) z = (c&8)l/2Xl"-1t; x, = t X 1 ; 1( = (k + 1)/2
and, because of the definition of n and the bounds on k,
< <
with 0 [ 1. The sign ambiguities in (9.11) and (9.12) are readily re-
solved with the result that the +
sign must be chosen in (9.11) and the
- sign in (9.12) [W].
The question whether there exists a real to>0 such that E ( T ~ )= 1
is now easily answered by observing that (9.11) exists when E = 1. In
fact, as observed by many authors 131, 32, 10, 331 under the change of
variable
(9.13) r& --s
(9.16)
2. The Inversions
From the above results it follows that the solutions of (9.8) are, for
case I
(9.16) z=
1 ’
(- )
1
Bp 2n’T ’
1
(9.17) t= &[B 1 1
-
with
(9.18) [* = €11”.
But, if € ( t is
) periodic, the solutions (9.16) and (9.17) cannot be single-
valued on an interval exceeding one-half period. Hence, it is desirable to
invert these solutions. The definitions used in the inversions, and the trans-
formations are summarized in the table below. Cumbersome sign distinction
are avoided if “n is regarded to behave like an odd integer.” This phrase
+
means that negative quantities, raised to the power fi p , where 9 is either
zero or an even integer, remain negative, and negative quantities, raised
+
to the power n q where q is an odd integer, become positive. The
absolute values of these powers are given by
(9.19) ,*+P = l,P+P
Case I Case I1
NONLINEAR VIBRATIONS 196
3. Thc Am9-Functions
When n = 1,
(9.21) amp %I+ = '%,a.
Particular values are given by
amp nu* = 4 2 ,
(9.22)
amp0 = 0,
where
(9.23)
du,
=1 when It-1, for u,=O
=oo when n < 1,)
for %=y+.
(9.26)
I
d
-(ampn%) = 0, for %= 0,
du,
d
-(ampn%)=O
dust
=1 for %=%+
=do when n t 1.
,
(9.27)
The sam- and cam-functions are periodic of the same period; i.e.,
(9.28)
+
sam (ml)= sam (nu1 4nu*),
cam ( 2 ~ =
~ cn
) (2-1/2,~a).
If ulo and U s o are those values of u1 and ugfor which amp nulo = amp nuaO,
(9.31) +
Sam2*(nulo) cam& (nuso)= 1,
thus generalizing a well-known identity of trigonometric functions. The
I
derivatives of these functions are
d
-
dU,O
[sam (nulo)]= cam" (nuno),
(9.32)
d
K~[cam (nuno)]= - sam" (nulo),
I
-
d4
[sam (aul)]= - n sam%- 1 ( n y ) ,
(9.33)
da
7[cam (*US)] = - n cam*"-' (w).
The sam- and cam-functions look as shown in Figure 8 and Figure 9 for
a variety of values of n in 9 n < <00.
198 R. M. ROSENBERG
For n = 1, the curves are the trigonometric sine and cosine curve,
respectively. For n = 4, they are parabolas, and for B = 00, they consist
of segments of straight lines.
It is evident that the solutions of (9.8)are for Case I and 11, respectively,
.f = Sam (m),
(9.34)
.f = cam (n7),
as one can readily verify by direct substitutions, and by making use of (9.33).
NONLINEAR VIBRATIONS 199
(10.4)
where the masses and spring parameters of the perturbed system are
(10.6)
*i = mi + &a& (i = 1,.. .,?$),
Gj = uj + && .
( j = 1,. . J ) .
While some of the ai and Pi may be zero, all are fixed, once chosen.
200 R. M. ROSENBERG
possesses, at some energy level h, a straight modal line, then the fierturbed system
moving at the same energy level h, possesses a wnique modal line in an &-neigh-
borhood of that of the parent system.
The proof of this theorem is lengthy and will not be reproduced here.
In view of this theorem, the modal line of the perturbed system may be
written as
(10.8) x&) = cixl + E&~(X,) + e*&(xl) + . . ., (i = 2 , . . .,n)
and this is valid everywhere except, possibly, on the La-surface where the
system M is singular.
Then, substituting (10.8) in (3.12); i.e., in the system MI and expanding
in powers of E , it turns out that must satisfy the system of linear, inhomo-
geneous, second-order equations
(10.9)
where
(10.10)
and a sirn..ar notation is used for oo* and for the partia. derivatives of
these quantities. Moreover, it turns out that hi must satisfy a system
of equations in which the left-hand side is identical to (10.9) while the right-
hand side becomes a function the tlj. Consequently, one can find &, if one
can solve (10.9), and we need concern ourselves only with (10.9); i.e., with
first-order perturbations.
NONLINEAR VIBRATIONS 201
3. Integrable Cases
(10.11)
(10.16) =5
a@,) h= 1
ch2 5
XI
0
(cj awe*
ax,
- w)
axj
dx,, (i = 2,. . .,%).
with two degrees of freedom [8], and by Kuo [18] for those with any finite
number of degrees of freedom, that the Ai = 0, (i = 2,. . .,n)as well, when
coupling is weak. The arguments leading to this result are lengthy and
not simple, even in systems with only two degrees of freedom, because the
system of equations M is singular on the La-surface. Therefore, we shall
not reproduce them here. However, the result leads directly to
Theorem X-11: The slope at the origin of the modal line of the Perturbed
system i s the same as that of its weakly coupled parent system, on which it
neighbors.
It follows that the modal line of the perturbed system is given, to first
order in E , by
(10.17) xj(xl) = cix, + e&(xl), (i = 2,. ..,n)
where
(10.18)
4. General Properties
Consider now (4.4) with f,(t) 0; i.e., the equation determining the
period of a normal-mode vibration. It is well known that the period T,
given in (4.9), is independent of X, if (4.4)is linear, but it does depend
on X, if (4.4) is nonlinear. Summarizing these results leads to an interesting
gradation of normal-mode vibrations. In it, we replace the words “modal
line” by “mode.”
(i) In linear systems, the normal-mode vibrations are always similar;
the normal-mode coordinates decouple the equations of motion for arbitrary
motions; mode and period of the normal-mode vibration are independent
of the energy level of the motion. Linear combinations of solutions are also
s o htions .
(ii) If the system is nonlinear, the normal-mode coordinates decouple
the equations of motion for that mode only. If the normal-mode vibration
is similar, the mode is independent of the energy level, but the period is
not. Linear combinations of solutions are no longer solutions.
(iii) If the system is nonlinear and the normal-mode vibration is non-
similar, mode as well as period depend on the energy level; otherwise, the
conclusions are the same as in (ii) above.
(11.1)
(11.3)
and that the Xiw-planes contain certain singular points. The number,
location, and character of the singular points of (11.3) can then be used to
determine the geometrical character of the response curves. In particular,
one can ascertain where (if a t all) resonance occurs, and at what frequencies
(if a t all) the phenomenon of “tuned dynamic vibration absorption” takes
place.
2. Homogeneom Systems
We shall show that the linear system under simple harmonic excitation
is merely a special case of the homogeneous system under the excitation
V1(t)lkwhere
(11.6) f i ( t ) = P,cam (Ha),
NONLINEAR VIBRATIONS 206
(11.8)
where
(11.9) D = (cX)~-*,
and c is a constant. Form (11.8) has a certain advantage over (11.5). In
(11.8), the quantity I plays a role similar to that of the circular frequency
of simple harmonic excitation, but it is not the frequency of excitation when
k > 1. Hence, I will be called the generalized circular frequency.
Then, the equations of motion of the physical system are
(11.10)
m,ii + aluk + a& - v ) k = P, camk (n~ / ( D A ~t )/ ,N )
m,i + a3vk- a2(u- v ) =~ 0,
those of the pseudo-system are
(11.11)
x+-=
.* ;L( Vm1
m, 7(
+=-=-=
m1 Vm1 VY_,)*=+camh(nEt),
( 11.12)
Simple substitution of
(11.13)
y = cx
into (11.11) and (11.12) (and making use of the formulas (9.33) for the
derivatives of cam-functions) shows that (11.13) are indeed the solutions,
and the constant c in the second of (11.13) is the same as that in (11.9).
But it follows from (11.13) that, in the case considered, x(t) and y(t)
are cam-functions of the same argument. Hence, the steady-state vibrations
are similar, or the trajectories satisfying (11.12) are straight lines.
If we define
(11.14) Y = cx.
the equations (11.11) and (11.12). with (11.13) and (11.14) substituted in
them, are three equations in the three unknowns X,Y and c, and they can
be solved explicitly for these quantities.
One finds
c = aJl*/{a,'Ik + (a, - ~ J h + W J ~ ) 1 / b } ,
(11.16) Xk = { P o m ~ / 2 [ a , + * - ?%a(~+~l/*A*)1P]k)/R,
~ / (as
Yk = { Poa,m,RI,}/R
where
NONLINEAR VIBRATIONS 207
(11.17)
(11.18)
x1 = x,cam (n YFt)
Xj(t)/Z,(t) cj = X j / X , , (i = 2 , . . .,m)
satisfy this system as well as the equations of the trajectory which are not
reproduced here. Now, there are m quantities Xi and m - 1 quantities
cj which constitute the unknowns, and there are m equations of motion
and m - 1 equations of the trajectory. Hence, the problem is determinate
WI.
It is now advantageous, as done in the linear forced vibration problem,
to introduce non-dimensional quantities as follows:
(11.19)
y = P/[a,/m,(k+1)q
where it should be remembered that X and P are nondimensional ampli-
tudes, and y is a nondimensional generalized frequency ratio. With the
introduction of these quantities, (11.15) become
where
208 R. M. ROSENBERG
One must now find the locus of the endpoints of the trajectories, i.e.,
the L-surface, which is, here, a curve. The L-curves are found as follows:
the first of (11.20)is solved for y as a function of E. This is substituted in
the second, thus eliminating y , and resulting in an equation of the form
(11.22) xk = qE).
If one substitutes in that equation
(11.23) E= PIR,
one finds an equation in x
and p ; this is the equation of one of the loci
of end-points of the trajectories. The second locus is found by substituting
x
- for 19; then, the equations of the two loci become [36]
(11.24) x k y + aa1(8- p ) k ( p+ p:2x\) - a z l a , , p l ~ ~ ~ *p = 0.
2
From these equations one can readily deduce the following properties of
the L-curves :
2
(11.26) ple = 2, aB1= 112, as2= 3, k =3
and for the associated linear problem in which all parameters have the
values given in (11.25),except that k = 1. The results in the 27-plane are
shown for these two examples in Figures lO(a) and (b), respectively.
For an understanding of the general theory, to be presented in the next
section, it is helpful to interpret these diagrams. In general, the trajec-
tories of steady-state forced vibrations are, here, segments of straight lines
NONLINEAR VIBRATIONS 209
( b)
FIG. 10. Trajectory of steady-state forced vibration in (a) nonlinear and (b) corre-
sponding linear systems.
210 R. M. ROSENBERG
that pass through the origin and terminate on the L-curves. Their slope
is a function of the frequency y and, in fact, their slope E(y) is given in the
first of (11.20).
Assume that y, i.e., essentially the frequency of excitation, is small,
so that E is less than the slope of the modal line in the first and third quadrant
(shown dotted). As y is increased, the trajectory begins to rotate in a counter-
clockwise direction and becomes longer, tending to infinity as the slope
E(y) approaches that of the modal line. As y is further increased, the trajec-
tory continues to rotate, but it becomes now shorter until E = 00. This
occurs where G ( y ) = 0 and, at this frequency, there is bounded 9-motion,
but no 3-motion. In other words, the +motion vanishes identically because
m2 acts as a tuned dynamic vibration absorber for m,. A further increase
in y will further rotate the trajectory, and the motion becomes out-of-phase
because the slope ?< 0. At first, the trajectory will lengthen until the mo-
tion becomes, again, unbounded when y = yz. However, further increases
(11.26)
where
(11.27) F ( y )= ~ G ’ o v w- G(YW’(Y),
and G and H are defined in (11.21).
Both of (11.26) are seen to have singular points on the y-axes a t the
zeros of H ( y ) . The first of them has additional singular points on the y-axis
at the zeros of G ( y ) .
Now, the generalized frequency ratios of normal-mode vibrations of the
autonomous system are, in general, functions of the amplitudes R and P,
and they tend to certain limits as the amplitudes tend to zero. It can be
shown [12, 131 that their limiting values coincide with the zeros of H ( y ) .
In other words, the equation
(11.28) H(Y) = 0
NONLINEAR VIBRATIONS 211
can be shown to play the part of the so-called “frequency equation” of the
autonomous system. Hence, both equations (11.26) have singular points
at the natural frequencies of normal-mode vibrations at vanishingly small
amplitudes.
Denote the zeros of H(y) by y,,. Then expanding the first of (11.26)
in the neighborhood of (K= 0, y = y,,) one finds
(11.29)
because
(11.30) fi(yn) = - G(y*)H’(yn).
Expanding the second of (11.26), one has directly
(11.31)
(11.33)
because
(a) ( b)
<
Theorem XI-11: The system (11.11)i s szcch that, for any k ifi 1 k < 00,
there exists a single nondimensional generalized fveqzrency ratio y* of the ex-
citing force swh that the mass m, acts as a tuned dynamic vibation absorber
for the mass ml that i s bcilrg excited.
Response curves in the 8 y and Py-planes for the values of the parameter
given in (11.26) are shown in Figures ll(a) and @), respectively.
It remains to map the Xy and Py-planes on the KQ and PQ-planes,
where the frequency Q of the forcing function is defined through (11.6).
It follows from (9.23) and (11.19) that this mapping is given by [%I
where T is the only period of f . Hence, the circular frequency of the periodic
force is
(12.7) w = 2n/T.
for U and t,4 and their partial derivatives, and using the definition of F(x,,X,)
given in (3.10), the equations of the tj turn out to be
216 R. M. ROSENBERG
(12.14)
But, this system of equations is of the form (lo.$), and questions of its
integration are discussed in the section on nonsimilar normal-mode vibrations,
and also by Kinney [12] for different examples.
Once the integrals #o)(x,) have been found the functions (12.11)
(12.16) xj(x1) = cjxl + E#')(x~), (i = 2,. . 1%)
are substituted in the first of the equations of motion; that equation can
then be integrated by quadratures for the conditions
3, = 0 when x1 = X,,
(12.16)
xl = XI when t=0.
(12.17) x, = X,(t,X,)
and its inverse
(12.18) t = tl(xl).
forced vibration pass through the origin and are monotone because they
are those of a vibration-in-unison; moreover, they must lie wholly in the
&-tubesbecause of (12.11). Consequently, the L-surface which is the locus
of the endpoints of these trajectories, also lies inside the &-tube.
The direction with which the trajectories pass through the origin is a
function of the frequency w of the exciting force as explained in connection
with Figure 10. If that frequency lies near that of a normal-mode vibration,
the trajectory lies near a modal line and is long. If not, the trajectory is,
generally, short, and the amplitudes of the forced vibration are of O(E).
In this way, one understands readily the occurrence of resonance.
It is also instructive to translate this information into the frequency-
amplitude planes, shown in Figure 14. The dotted lines, composed of the
w-axes and the backbone curves [29, 381, are those of every vibration-in-
unison of the autonomous system. In other words, if the autonomous system
does not movein normal-mode vibrations, it is at rest with respect to admissible
motions. (Because of Definition I, every vibration-in-unison of the auto-
nomous system is a normal-mode vibration.)
Hence, we may surround the w-axes and the backbone curves by &-tubes
as shown in the diagram, and the forced steady-state vibrations of the non-
autonomous neighboring system must be represented by response curves
lying inside these &-tubes. One such response curve, for a prescribed, small
force-amplitude is shown. It is evident, then, that the intersections of
218 R. M. ROSENBERG
w-axes and backbone curves are saddle points. In addition to these, there
may be star points as well, if the steady-state forced vibrations are similar.
If they are not, one can easily show that the only singular points in the
XIII. STABILITY
1. Definitions
Two definitions seem particularly useful in the examination of the
stability of vibrations-in-unison : “stability in the sense of Liapunov,”
or L-stability and “stability in the sense of Poincarb” or P-stability [39].
To define the former, let
x; = % i +. (.t ) , I
(13.1) (i = 1,. ..,n),
x;*(O) = xi, &*(O) = 0,
NONLINEAR VIBRATIONS 219
1
xi = .j*(X1),
(13.5) (i = 2, . . .,N),
Xi(&) = xi,
(13.6) xi = X i ( % , ) , (i = 2 , . . .,n)
denote any set of solutions of the system M. Then, the set xi*(xl) is called
P-stable if, for every E with 0 < E << X i , (i = 1,. . . ,n),there exists a 6(e)
<
with O < S(E) E such that
(13.7) Ixi(x1) - xi*(xI) I < E
whenever
(13.8) Ixi(x,) - Xi*(X,) I< 44.
The interpretations of these two definitions of stability are clear. Let
.
the space having coordinates (x,,. .,x,,,t) be called the event-space. Then
..
x,*(t), (i = 1,. ,N) is a curve in the event-space. If any solution xi@),
having a single point inside an &-tubeabout xi*(t) remains within it every-
where, the solution is L-stable. Similarly, xi*(xl) is a curve in the con-
figuration space. If any solution xi(xl), having a single point inside an
e-tube about xi*(xl) remains within it everywhere, the solution is P-stable.*
It is evident, therefore, that every L-stable solution is also P-stable, but
the converse need not be true. Moreover, every solution that is P-unstable,
is also L-unstable but, again, the converse need not be true.
(13.10) i'+B(t)E= 0
it turns out that the matrix B(t) is of the form
(13.11)
where B, and the Bi are square matrices with constant elements, and the
fi(t) are periodic scalar functions of time. Then, if a transformation matrix
T exists such that
(13.12) T-'B(t)T = d ( t )
NONLINEAR VIBRATIONS 221
3. Homogeneozcs Systems
It is simple to demonstrate [7]
Theorem XIII-11: The system of variatioml equations of admissible,
homogmeous systems with respect to any normal-mode vibratioH can dways
be decwpled.
To prove it, one shows that, in homogeneous systems, (13.11) can always
be reduced to
(13.13)
(13.14)
%ii = - b#(v) + A#(% - v )
where
(13.16) $(w) = wlwl*-l.
If the trajectory of a given normal mode is
(13.16) v = cu;
222 R. M. ROSENBERG
(13.17)
one can show that the variational equations have solutions that satisfy [41]
Bl = C,
(13.19)
B, = - ml/cm,.
The variational equations reduce to
(13.20) € +4 s $ ( q ) € =0
D 1 --- a
+ -((IA - c)Il - c(k--1,
m1 m1
(13.21)
But, the independent variable in (13.20) is t, not t. Now, one has from
(13.22)
(13.23)
(13.26)
where Dlg = &/E. It follows that the two equations (13.26) do not
contain the amplitude X.
The discussion of the characteristic exponents of the equations (13.26)
is difficult because they are Hill equations in which the periodic coefficients
axe cam-functions. The detailed stability properties of such Hill equations
are not known.
-1 0 1 2 3 4
However, one can gain some insight into the stability of normal-mode
vibrations of homogeneous systems by replacing the cam-function in (13.25)
by a cosine-function. Proceeding in this manner is, in fact, equivalent to
expanding the periodic coefficient into Fourier series and retaining only the
first term of the expansion. This is usually done in examining the details
of stability of vibratory motions of nonlinear systems [4]. It has the con-
sequence that equations (13.26) become Mathieu equations whose stability
properties are well known [a, 421.
224 R. M. ROSENBERG
(13.29)
k
B = - [a + A l l - cIR-1(1 - c)].
NONLINEAR VIBRATIONS 226
For any given system parameters] the formulas (13.29) can be utilized to
determine the location of P(S,&). Similar results are found, when
(13.30) B=B,=-m 11cm8'
In either case, it turns out that the point P(e,S) lies on a straight line in
the &plane which passes through the origin, which lies always in the right
half-plane, and whose slope is a function of k only. When K = 1, this slope
is zero so that the point lies then on the (positive) S-axis which consists
entirely of stable points. Thus, linear systems are always stable. When
k > 1, the slope of the straight line is positive, and when O < K < 1. the
slope is negative. The location of P(8,e) on the straight line is a function
of the system parameters [28].
(13.31)
(13.33) e(x*(o*r*(w= Q*
for U and its derivatives, it is easily shown that, in symmetric systems,
one has always
(13.34)
*
uxx = uy,
*
for both, the i and o-modes.
The coupled variational equations are
.. * *
E - uxxt - u x y q = 0,
(13.36) * *
q - UXYE- U,,q = 0.
Under the transformations
tl = 'I+ 5,
(13.36)
c, = q - E,
226 R. M. ROSENBERG
(13.311
(13.38) u = u ( L ) + U(N)
where
(13.39)
(13.40)
Note that the quantities arising from the linear part of the potential function
are the only ones independent of the mode.
It will be remembered from the earlier example that the nonlinear
vibrations x*(t),y*(t) will be replaced, in the vuriatiod eqaruhzs, by their
zeroth approximation; i.e., by simple harmonic functions. It is entirely
consistent with this approximation, to replace the frequencies of the normal-
mode vibrations by their Duffing approximations [4]. These are, in terms
of the quantities defined in (13.40)
t The superscripts (L)and (N)indicate the terms in U which give rise, respectively,
to linear and to nonlinear terms in the equations of motion.
NONLINEAR VIBRATIONS 227
If one utilizes
(13.42)
(13.46)
In (13.44) and (13.46), the + signs are used with the subscript 1, and the
- signs with the subscript 2.
Evidently then, (13.43) represents four equations, those for two modes,
and two equations for each mode. As in the earlier example, one of the
two equations in each mode defines points in the Strutt chart which lie always
on the boundary between stable and unstable regions of the &plane. Hence,
it is the other two equations which determine the stability; these are the
combinations (1,o) and (2,i).
In applying the foregoing formulas to the example of potential function
(13.38) and (13.39), it is convenient to introduce the nondimensional quanti-
ties
228 R. M. ROSENBERG
(13.47)
and a similar result in the o-mode. In both cases, the parameters &,8depend
on E ; i.e., on the amplitude of the motion.
FIG. 16. Stability of in-phase mode of nonlinear, symmetric system with two
degrees of freedom.
The results for the i-mode are shown in Figure 16 and are as follows:
(i) The point on the Strutt chart, deciding on the stability of the i-
mode lies on a straight line passing through the point (C#, and intersecting
the positive &axis at a point that depends on oc, only. This straight line
is called a stability chracteyistic.
(ii) The portion of the stability characteristic lying in the upper half-
plane applies yhen u8 > 0, that in the lower half-plane when a, < 0.
NONLINEAR VIBRATIONS 229
FIG. 17. Stability of out-of-phase mode of nonlinear, symmetric system with two
degrees of freedom.
The results for the o-mode are shown in Figure 17. They differ in some
respects from those for the i-mode.
(i) There exist, again, straight stability characteristics. These terminate
on a point lying on a straight line of unit slope, passing through the origin.
The point on this line. at which the stability characteristics terminate.
depends on a, only. Further, the stability characteristics intersect the
positive &axis at a point whose value depends on a, only. Thus, know-
ledge of a, and ol, determines the stability characteristic and the amplitude
of the motion determines the point on the stability characteristic. Items
(ii) and (iii), listed above for the i-mode also apply to the 0-mode. From
these results, knowledge of El, a,, a3 and the amplitude of the motion permits
the determination of its stability.
230 R. M. ROSENBERG
(13.61) 0 = V ( X ~ .,xn)
, e + E O ( X I , . . -,xn)j (14 IUl).
The equations of motion of the perturbed system are
(13.62) . + .
Y i = U,(X,,. .,xn) ~ o , ~ ( z ~.,x,)),
,. (i = 1,. ..,a),
and it has been shown in a previous section that this system has a normal-
mode vibration (in general nonsimilar) of the form
(13.63) g6(t)= xi*(t) + e&(t) + .. ., (i = 1,. ..#n),
Then, one can easily demonstrate [8].
Theorem XIII-111: The stability of the nomttal-mode Vibration (13.53)
of the Perturbed system i s the same as tlte stability of the similar n o d m o d e
vibratiolz (13.60) of the pareNt system.
For the proof, one forms the equations of the first variation of (13.62)
with respect to the solution (13.53). It turns out that these are precisely
the same equations as those of the first variation of the equations of motion
of the parent system with respect to the solutions (13.60). which proves the
theorem. Hence, if the variational equations of the parent system can be
decoupled, the stability problem of normal-mode vibrations of neighboring
systems is solved.
6. Forced Vibratiows
where use was made of notation (13.33) with respect to the second deriv-
atives of U and 4. Since ZJ is a quadratic form, the a W / a x i a x , are constants
for all i and k.
We Write (13.69) in the matrix form
(13.60) 2 + [B, + sB(l)]z= 0
where z is a column matrix, B, is a constant square matrix,
S
(13.81) B(1) =
8-
2 .. B,cosswt,
I$,.
It will be noticed that, here, we have not made the assumption that the
solutions (13.68) are to be treated as simple harmonic for purposes of a
stability analysis. In fact, that assumption is unnecessary here because the
approach through Mathieu equations is not used.
It is now supposed that B, and the B, do not satisfy Hsu's criterion
[40] and that, in consequence, (13.69) cannot be decoupled.
We introduce the transformation
(13.62) z = TE
g+(n+EzD,COSSd 1
s
(13.63) l=o
s-1
where A
2 is a diagonal matrix with elements
(13.64) ~ 1 < wg' < . . . < wn2,
2
and the
(13.66) D,= [d$'], (S = 1,. . . ,S)
are constant matrices.
Then, Hsu's results show that instability occurs when, and only when,
the following inequalities are satisfied :
(13.66)
Now, the functional relation between the frequency w of the exciting func-
tion and the amplitude X, of x,*(t) is known. Hence, the inequalities (13.66)
determine the amplitudes X , (if any) for which the motion is unstable.
Kinney [12] has applied this theory to the weakly nonlinear problem
of the dynamic vibration absorber. The analysis is too detailed and tedious
to be reproduced here because of the involved form of the matrix elements
dii'. His result for that problem is:
The portions of the response cwves (in the freqwncy-am#Littde planes)
which have negative slope cmespond to unstable motion. All other portions
of these cwves correspolzd to stable motion.
NONLINEAR VIBRATIONS 233
dx,
for all xl, is a modal line.
With the above definition we prove in this appendix the central
Theorem A-I: Every admissible azctonomotrs system possesses at kast
one simple trajectory.
Before presenting the proof of this theorem it may be helpful to sketch
its essential ideas; these are very similar to those used in the two-dimensional
problem [191.
We consider trajectories in configuration space which pass through (or
issue from) a point Po on the bounding surface; such trajectories are called
T-curves. We then demonstrate that every T-curve passes through an
infinity of isolated distinct points Qi(i = 1,2,.. .), corresponding to con-
figurations of stationary potential energy with respect to neighboring
points on T. The corresponding levels of potential energy are denoted by
hi, and the transit times through the Qj by t i . We show, further, that the
ti and hi are uniquely determined by the initial point P o , and both are
differentiable with respect to P o . Next, we demonstrate that there exists
at least one T = T*, issuing from a point Po* for which the ordered
Qi(i = 1,2,. . ,) lie alternately at the origin and at Po*. Finally, we show
that T* satisfies condition (iii) above ; this completes the proof.
(iv) Trajectories are the vectors q(t) which satisfy the equations of
motion,
a
miqi = - U(q).
%i
(i = 1,. .. ,n)
where the mi are real, positive constants. The equations of motion transform
under
-
(A4 x; = vmiqi
into
.. = -
X;
a U ( X )= U,i
ax,
or, in vector notation,
(A.3) P = VU.
With Kuo [18], we transform (A.3) into the system of 2n first-order equations
(A.9 ji = fib)
where
Yi = xi, yi+u = 2i I
fi=Yn+i, f,+n= UIi.
NONLINEAR VIBRATIONS 236
Since the system is admissible, the functions f j are of class C1, or the Jacobian
J is continuous, where
..
For proof of the second lemma one need only observe that trajectories
which attain the bounding surface intercept that surface orthogonally
(Theorem V-V). Hence, the initial conditions (position and slope) are
defined; then, in view of the uniqueness of (A.6), the trajectory must retrace
itself.
I t is now convenient to rewrite the energy integral in terms of the initial
energy when the initial point is not on the bounding surface:
(A.lO) - + tl141s= - U X ) + i l l q *
W)
where x and as functions of the initial values X and d are given in'(A.6)
because y and Y correspond, respectively, to ( x , i ) and (X,ft).
Differentia-
ting (A.lO) with respect to X and ft, we have:
-ax,
__ ... ax, _an,_...- an,,'
ax, ax, ax, ax,
ax,
- ax,
(A.ll)
ax,, axu
ax, ax, aa, an"
_-
ad, ax, a f t , aft,
and we shall demonstrate the existence and regularity of A(t), because these
are essential for the proof of the central theorem. For instance, the transpose
A* may be found from
(A.13)
with the initial condition A*(O) = I where I is the identity matrix. Since
J is of class C1, the existence, uniqueness, and continuity of A(t) is assured.
But J is a function of the configuration x ; hence the solution of (A.13)
depends on the initial configuration X, and we write it as
NONLINEAR VIBRATIONS 237
(A.14) A = A(t,X).
Let us consider A as a mapping in the phase space which, for fixed t and X
maps the point (- VU,$ onto the initial point (- Po U , x ) where voU
denotes grad U at X . Then, A has the important property, formulated in
Lemma A-111: For fixed t and X , the mapping A i s volume-preserving;
i.e., the determinant of A is unity. Conseq.uently, a unique inverse of A exists.
To prove this lemma, we make use of the Jacobi-Liouville formula
(A.16)
I
det IA*(t)I = exp tr J(s)ds
0
where J is defined in (A.6), and tr J denotes the trace of J . Now, the elements
in the main diagonal of J are all zero. Therefore, det A is unity and the
inverse A-l exists. Thus, we may write (A.ll) as
(A.17)
au - 0
--
at
(A.18)
238 R. M. ROSENBERG
for variations in transit time through a given energy level, and for variations
in the initial point. The existence of this differential is assured by the prop-
NONLINEAR VIBRATIONS 239
erties of U and the uniqueness of (A.6). Consider now one of the hj , say h, .
Its differential is the same as (A.20) except that U,= 0 by definition.
Hence,
(A.21)
Thus, the hi are continuous and differentiable with respect to the initial
point Po(X).
To prove the existence of a simple trajectory we extremize h j ( X ) subject
to the equation of constraint
(A.22) U(X) + h = 0.
In other words, we determine the stationary values of the potential energy
for a trajectory and we insure by the constraint equation (A.22) that this
trajectory is a T-curve. Thus, we define a function
(A.23) g ( X ) = hj(X) - ilU(X)
where il is a Lagrangian multiplier. Since we wish to deal here with the
phase space we may write, instead of (A.23)
(A.26)
(A.26)
(A.27)
(A.28)
240 R. M. ROSENBERG
The first of these states that U,= 0 at the origh since VU = 0 at the origin
only, and the second, that Ut = 0 on the bounding surface because that
surface is defined by R = 0. Hence, a T-curve T* exists which connects
an initial point on the bounding surface with the origin. To prove Theorem
A-I, it remains to show that the T*-curve is a simple trajectory.
Let Po* be the initial point of T* on the bounding surface. Let its points
of first, second, etc. tangency with equipotential surfaces be Q1,Q2,. .. .
Now, suppose that Q1 is not a t the origin. Then, in view of the above
demonstration, Q1 is on the bounding surface. In that case, T* moves at
first toward configurations of lower potential energy and then toward
configurations of higher potential energy. Thus, there is a point R between
Po* and Q1 at which U,= 0. But this contradicts the fact that Q1 is the
first point along T* where Ut= 0. Hence, Q1 is a t the origin. This com-
pletes the proof.
ACKNOWLEDGMENT
This work was supported, in part, by a grant from the National Science
Foundation.
Retmcnces
R. M. and Kuo, J .
8. ROSENBERG, K.,Nonsimilar normal mode vibrations of non-
linear systems having two degrees of freedom, Trans. ASME. Journ. Awl.
Mech., 81, E, 2, 283-290 (1964).
9. ROSENBERG, R. M. and Hsu. C. S., On the geometrization of normal vibrations
of nonlinear systems having many degrees of freedom, in Proc. IUTAM Symp.
on Nonl. Vibr.. Ukrainian Academy of Sciences, Kiev, 1. 380-415 (1963).
10. KAUDERER, H., Nichtlineare Mechanik, Springer, Gtjttingen, 1968.
11. ROSENBERG, R. M., The normal modes of nonlinear n-degree-of-freedom systems,
Trans. ASME. Journ. A w l . Mech., 29, E, 1, 7-14 (1962).
12. KINNEY,W. D., On the geometrization of the forced oscillations of nonlinear
systems having many degrees of freedom, Ph. D. Thesis, Univ. of Calif., Berkeley,
1905.
13. KINNEY.W. D. and ROSENBERG,
R. M., On steady-state forced vibrations in
strongly nonlinear systems having two degrees of freedom (Les vibrations
forcbs dans les systbmes nonlin6aires). Proc. CNRS I s t . Coll. 148, 351-373
(1965).
14. DARBOUX, G.. “Thbrie a n b r a l e des Surfaces,” vol. 2, pp. 438-509. Gauthier-Villar,
Paris, 1889.
16. J ., “Oscillatory Motions” (Translation of “Les Mouvements Vibratoires”),
HAAG,
vol. 1, pp. 36-39, Wadsworth Publ. Co., 1962.
16. ROSENBERG, R. M., La g6ometrie de la dynamique et les vibrations des systbmes
non-linbaires. L’Onde Ebctrique 48, 433, 425-434 (1963).
17. COURANT, R. and HILBERT,D., ,,Methoden der mathematischen Physik‘, vol. 1,
pp. 24&245. Springer, Berlin, 1931.
18. Kuo, J . K., On non-similar normal mode vibrations of nonlinear systems having
many degrees of freedom, Ph. D. Thesis, Univ. of Calif., Berkeley, 1964.
19. ROSENBERG, R. M., On the existence of normal mode vibrations of nonlinear systems
with two degrees of freedom, Quart. Appl. Math., 22, 3, 217-234 (1964).
20. ROSENBERG, R. M., The ateb(h)-functions and their properties, Quart. Appl. Math.
el, 6, 37-47 (1963).
21. ROSENBERG. R. M., “Les vibrations non-lineaires de systbmes B plusieurs degrbs
de libert6.” Centre de Rech. Phys., Marseille, 1984.
22. SIDERIAD~S, L. and ARGBMI,J.. Sur la singularit6 multiple des systhmes dyna-
miques plans, Comp. Rend., Acad. Scie. 268, 2037-2039, Gauthier-Villar, Paris
(1961).
23. ARGEMI,J. and SIDERIADkS, L., Sur la singularitd multiple dbfinie par les formes
homogbnes des systkmes dynamiques plans, Comp. Rend. 864, 4135-4136.
Paris (1962).
24. ROSENBERG, R. M., On motions with a rest point and the existence of normal
mode vibrations in nonlinear systems, in “Proc. Eleventh. Int. Congr. Appl.
Mech., Munich” 1964, (to appear).
25. MAHAREM, N., On homogeneous systems and the stability of their normal modes,
Ph. D. Thesis, Univ. of Calif., Berkeley, 1965.
26. CODDINGTON, E. A. and LEVINSON, N., “Theory of Ordinary Differential Equations,”
McGraw-Hill, New York, 1966.
27. HAUGHTON. K. E., Similar motion of multi-degree-of-freedom vibrating systems
with nonsymmetric springs, Ph. D. Thesis, Univ. of Calif., Berkeley, 1964.
28. ROSENBERG, R. M., On normal vibrations of a general class of nonlinear dual-mode
systems, Trans. ASME. Joum. Awl. Mech.. S8, E, 2, 275-283 (1961).
30. ROSBNBERG, R. M. and ATKINSON.C. P., On the natural modes and their stability
in nonlinear two-degree-of-freedom systems, Trans. ASME, Journ. Appl. Mech.,
B6, E, 3, 377-386 (1959).
242 R. M. ROSENBERG
243
244 PIOTR PERZYNA
1. PHYSICAL
FOUNDATIONS
G. I. Taylor [171], was able to perform tensile tests over a wide range
of average strain rates by firing bullets against an anvil bar held as a simple
beam between two tensile specimens. The other end of each specimen was
fastened to a ballistic pendulum. Since the velocity of the anvil after im-
pact decreased from the impact velocity to zero, he assumed a mean rate
of strain equal to 1/2(velocity of anvil)/(length of specimen). Applying
the law of conservation of energy and assuming that the kinetic energy
equals the plastic strain energy he determined for steel the variability of
the ratio of the dynamic to the static yield limit as a function of the mean
strain-rate, Fig. 1.
4:
9
*.
9
4 I I I 1 I I I I ,
FIG. 1. Effect of strain rate on dynamic yield in steel (G. I. Taylor [171]).
FIG.2, Effect of strain rate on the proportional limit and ultimate strength of a 0.22%
carbon steel (D. S . Clark and P. E. Duwez [39]).
FIG.3. Effect of stress on strain at constant strain rate of high purity aluminium (F. E.
Hauser, J. A. Simmons and J. E. Dorn [77]).
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 249
FIG.4. Super-purity iron based 0.21 per cent carbon steel curves: (a) strain rate against
strain for dynamic test; (b) stress against strain, A. Static, B. Dynamic; numbers
denote time in microseconds; impact velocity 198 in./sec. (J. Harding, E. 0.Wood
and J. D. Campbell [75]).
Stmin-per mi Strain-percent
FIG.6. Super-purity iron curves: (a) strain rate against strain for dynamic test; (b)
stress against strain, A. Static, B. Dynamic; numbers denote time in microseconds:
impact velocity 198 in./sec. ( J . Harding, E. 0. Wood and J. D. Campbell [76]).
with initial velocity of 400 inlsec, yielding occurs only in the neighbourhood
of end of the test piece under impact and is limited to about lOmm (Fig.6).
These tests showed that there is a strong dispersion due to rheologic effects.
The influence of the transversal motion on the phenomenon of propa-
gation of stress waves in a bar is discussed in a paper by H. Kolsky and
L. S. Douch [96].
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 251
The phenomenon of the delay time was studied in detail for mild steel
by J. E. Johnson, D. S. Wood, and D. S. Clark [87]. The delay time which
varies within the limits of 40 p e c to 1.5 msec is a function of the initial
impact stress varying from 80,000lb/ins to 5000 lb/in2, respectively. It was
also shown that for low carbon steel the delay effect in the dynamic com-
pressive test is of the same nature as in the static tensile test (Fig.7; cf.
also [28, 32, 38, 401).
o Cornpmsion impact
Rapid load tension
m-4 lo" 10-9 lo-' I
Delay time,sec
FIG.7. Stress versus log delay time, 73 "F (J. E. Johnson,D. S. Wood and D.S.Clark
~371).
FIG. 8. Static stress-strain curves obtained after dynamic loading. Curve A : annealed
specimen (not pulsed). Curves B, C, D: specimens pulsed at 277, 293 and 308 in./sec.,
respectively. Curves E, F, G, H : specimens impacted at 326, 340, 362 and 366 in./sec.,
respectively ( J . D. Campbell and C. J . Maiden [29]).
It is clear that in real soil the inelastic strain-rate tensor depends on the
time-history as well as on the path-history of the stress.
2-
/
/
a/’
0
co
.
€ ? - 0
’
if.-----
-
I
01 I I I I I I I I t
40-4 10-2 1 102 104
hte of shin (pemntqeper minute)
FIG.9. Relation between strength and rate of strain for (a) clay and (b) sand (dataof
D. W. Taylor, A. Casagrande and W. L. Shannon, from A. W. Skempton and A. W.
Bishop [167]).
where u is the nominal tensile stress and ef’ and if’are the nominal plastic
strain and strain rate, respectively. P. Ludwik [107], and H. Deutler
[68], on the basis of experimental investigations, and L. Prandtl [147],
on the basis of some theoretical considerations, proposed a logarithmic
+
expression for the function (see also A. L. Nadai [llS]).
In the analysis of the problem of propagation of plastic stress waves
in bars, the effect of strain rate was first taken into consideration by V. V.
Sokolovsky [161, 1621, and L. E. Malvern [110, 1111. The physical law
of L. E. Malvern can be expressed thus
t in psec
- 8
--
I I
I"
x (inches)
FIG.10. Strain distribution in a long bar (S. Rajnak, F. Hauser and J. E. Dorn [ISO]).
the actual stress and the stress computed from the static curve is the assump-
tion of L. E. Malvern [110]. This difference produces the strain rate in
agreement with the viscosity law. The elastic strain component is considered
to be independent of the strain rate. In agreement with these assumptions
266 PIOTR PERZYNA
and in contrast with the results of the theory assuming the static charac-
teristic as a basis for dynamic considerations, the front of a plastic wave
propagates with the same velocity as the front of an elastic wave.
The theoretical results obtained by L. E. Malvern show that in this
strain-rate dependent theory there is no region of constant strain in the
neighbourhood of the end subjected to the impact load. I t is worthwhile
to observe that in many experimental investigations indeed no such region
I 1 I I I
1 I I I
x (inches)
FIG.11. Strain rate distribution in a long bar (S.Rajnak. F. Hauserand J.E.Dorn [lSO)).
was found (see for instance the tests by J. D. Campbell [26], Fig. 6). A
detailed discussion of this problem can be found in the paper by H. G.
Hopkins [83]. S. Rajnak, F. E. Hauser. and J. E. Dorn [150], examined
theoretically the strain distribution in a long bar subjected to dynamic
load. As a basis for their study they took the physical relation of L. E.
Malvern (1.6), the function g(a, 8 ) being determined from experimental
investigations with aluminium. Their results confirm the validity of the
assumptions of the theory that accounts for the influence of the strain rate.
Fig. 10 represents the strain distribution along the bar, computed for various
times. Similarly, Fig. 11 represents the distribution of the strain rate.
Fig. 12 represents a comparison of theoretical and experimental results.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 257
x (inches)
FIG. 12. Final strain distribution in a long bar (S. Rajnak, F. Hauser and J . E. Dorn
[ 1501).
(1.7) +
dB = f(a,s,l)do g(o,E,B)dE + h(a,s,i)dt,
where the influence of strain acceleration is taken into account in addition
to that of strain rate.
Another interesting question is : what is the interpretation of the existing
physical relations in the light of the dislocation theory. This is the subject
matter of a paper by J. A. Simmons, F. Hauser, and J. E. Dorn [166].
A study by means of the dislocation scheme of the classical theory of K k r m h ,
Taylor, and Rakhmatulin and the viscoplastic theory in which the influence
of the strain rate is taken into consideration has shown the advantages of
the viscoplastic theory. It is hoped that some of the experimental tests
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 259
(1.9)
the same as in flow theory. Difficulties are encountered only if the varia-
bility of flow surfaces in the course of a deformation process of an elastic/
viscoplastic body is considered.
The distinction just introduced requires that the considerations con-
cerning the constitutive equations should be subdivided into two parts.
Thus, methods of description of elastic-viscoplastic materials will be dis-
cussed first in Secs. 2 to 6. Sections7 to 13 will be concerned with elastic/visco-
plastic materials and with the application of that scheme to the description
of the dynamic properties of plastic materials.
It will be assumed throughout the entire work that the displacement
gradients in the straining processes under consideration are infinitely small
(cf. the definition in the monograph by C. Truesdell and R. Toupin [183]).
11. CONSTITUTIVEEQUATIONS
Let us denote by cii the strain tensor and by uij the stress tensor in
Cartesian coordinates. For an elastic-viscoplastic body it will be assumed
that the strain tensor can be represented in the form of the sum
where E$, e,; 4j denote the elastic, viscous, and plastic strain components,
respectively.
Let us emphasize the fact that the equation (2.1) is not satisfied in
general. Its introduction will be treated as a fundamental simplifying
assumption.
The stress and strain deviator tensors used in what follows are defined
thus :
A. C. Eringen [68], and W.Noll and C. Truesdell [122]. They also contain
extensive lists of literature.
With the notations of [73] the physical relations between the tensors of
strain and stress in the viscoelastic region read
where p and K are the elastic constants of the material (shear and bulk
modulus).
Making use of this observation, the viscous strain components can be
separated off in a simple manner. From (2.4). (2.6) we find (cf. P. M.Naghdi
and S. A. Murch [llS])
t
262 PIOTR PERZYNA
A2
In the elastic case the plastic state produced will be represented by the
same point independently of the time in which the state represented by A is
reached, provided that the load path is the same. If the material is visco-
elastic the plastic state may be reached at different points A, or 4, say,
depending on the time in which the load path is through. This difference
is caused by the viscosity of the material and the dependence of the load
history on time. It is obvious that by passing through the path OA in the
same overall time but with different strain-rates at the same points of the path,
different yield limits will be obtained.
In order to describe the complicated problem of a viscoelastic material
becoming plastic the notion of flow wrface has been introduced by P. M.
Naghdi and S.A. Murch [118]:
(2.9) f = f (q,&GP) *
= / I ( E ~ ; ~which
) represents viscous effects and the strain-hardening para-
meter K. The latter depends, in turn, on the plastic strain and is defined in
the same manner as in the theory of plastic flow describing isotropic strain-
hardening of the material (cf. the definition (2.44)).
Let us now consider the time-Variability of a flow surface. The time-
derivative of the function f is
(2.10)
(2.11)
(2.16)
264 PIOTR PERZYNA
The expression (2.16) shows why the load criteria for the viscoplastic
case are changed, if compared to those known from the classical theory of
plasticity; the influence of the viscous properties of the material on the
yield condition is here responsible. There is also an influence of the strain-
rate on the yield condition, of which a detailed discussion will be given
later.
M. Reiner [152] discussed the problem of attainment of the plastic state
with a constant load, that is, the possibility of a material becoming plastic
in a creep process.
An original analysis of yield criteria or, more generally, criteria of ob-
taining a critical state was proposed by W. Olszak and 2. Bychawski [124,
1261. This idea is based on the observation that the fracture of some visco-
elastic bodies depends not only on the elastic energy but also on the rate a t
which the dissipation energy varies. Examples of bodies were analysed
where destruction occurs if the rate of variability of the dissipation energy
reaches a certain value. For such bodies the strain rate is decisive. For an
analysis of a sufficiently large class of real bodies the authors of [126] propose
the following condition for reaching the critical state* :
where W Eis the reduced elastic distortion energy and W , the rate of change
of the dissipation energy expressed in the dimension of energy. The quan-
tities w ~ , w E and k are certain constants characterizing the material.
In [126] several models are considered, the form of the criterion (2.16)
being analysed in detail.
* It should be observed, however, that the attainment of the critical state is under-
stood in [125] in a more general manner than the attainment of the plastic state. The
plastic state may constitute a particular case of the critical state.
FUNDAMENTAL PROBLEMS IN VISCOPLA STICITY 266
where t = 0 is the instant at which the external load increments are applied.
It is often more convenient to introduce two loading paths T i l ) , P i l )
and Ti2),Pie), which become different after t = 0 (Fig. 14). Then (2.18)
can be written in the form
(2.19)
(2.20)
S V V
266 PIOTR PERZYNA
f’ fk
FIG.14. Two paths of loading. Arrows join points of the same time.
Assuming that only homogeneous states of stress and strain are con-
sidered, we obtain from (2.19), on the basis of (2.20), the condition
(2.21)
i‘
0
[aij
(2) - a;{(1)3 [its’- 49at 2 0.
i *
0 t 4 tk yields
(2.23)
(2.24)
For a rapid load path the behaviour of the flow surface is exactly the
same as in the inviscid theory of flow for an ordinary path.
268 PIOTR PERZYNA
(2.28) fa = f(aii,&~I"),P('),K("))= 0
c----
-0
that the state A not only depends on the components of the stress tensor
but also on the quantities E $ @ ) , , ~ ( ~ ) , K (which
~), describe the history of the
loading process. The latter is characterized by the load path and is followed
by the varying state of the body at the time at which the point A is reached.
Let us now assume that a loading process takes place a t t > t, accompanied
for t , > t, by an increase of plastic strain ,4$ dt; unloading then starts so
that the stress state returns to the initial state a t t = tk. The application
of the inequality (2.26) to the load cycle at the instant t, t. <
t - t,, gives
For any rapid load cycle the flow surface is on the basis of (2.26) identical
with the instantaneous flow surface
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 269
(2.31)
- -
'a
(2.32)
Thus the inequality (2.29) for a rapid load path can be assumed in the same
form as in the inviscid theory of plastic flow. In particular, if ta is assumed
to be equal to zero and a$) is treated as an arbitrary point in the visco-
elastic region, identical with a,; we have
Let us again consider the load path in the stress space and assume that
the instantaneous stress point moves along the same path with various
velocities. The duration of each loading process will be different, and each
will involve different rheologic effects. The plastic state will therefore be
270 PIOTR PERZYNA
Assuming that the vector of plastic strain rate .4$ is directed along the
outer normal t o the instantaneous flow surface f = 0, P. M. Naghdi and
S. A. Murch [118] proposed the relation
(2.34)
where
(2.36)
(2.37)
(2.38) i$ = l o a/
if f < 0,
Ih(B(akl)) aai, if f = 0,
(2.40)
(2.41)
The function h plays now the role of the viscoplastic potential, and the
function f represents the yield criterion for the elastic-viscoplastic material.
In order to determine the coefficient N and the connection between functions
h and f certain conditions must be imposed on the constitutive equations
(2.41).
A possibility seems to exist to develop the theory of viscoplasticity in
a broader sense than that for a stable inelastic material. To do that, an
272 PIOTR PERZYNA
This work was not fully appreciated for a long time. It was only in the
years 1948-1960 that it attracted attention subsequent to the observations
of V. V. Sokolovsky [162], and later of L. E. Malvern [110], who showed
that the assumption of K. Hohenemser and W. Prager may be used as a
basis for description of certain dynamic properties of rate-sensitive materials.
This concerned, however, one-dimensional problems only (see Section I).
Further development of the idea of K. Hohenemser and W. Prager is
contained in the references [132-1341.
Owing to the assumption that the viscous properties of the materials
become manifest only after the passage to the plastic state and that these
properties are not essential in elastic regions, the basic concepts of the
description of viscoplastic properties differ from the methods explained
before. It will be assumed that the strain rate can be resolved into an elastic
and inelastic part
(2.42) j 6; + 6;.
ii=
The inelastic part of the strain rate, which is denoted by ig., represents
combined viscous and plastic effects.
Since the material has no viscous properties in the elastic region, the
choice of an adequate yield criterion will be much simpler than in the case
of an elastic-viscoplastic material. The initial yield condition, which will
be called the static yield criterion, will not differ from the known condition
of the inviscid theory of plasticity.
In order to keep our considerations sufficiently general we introduce a
static yield function in the form
(2.43)
where the function f(ui,&) depends on the state of stress uij and the state
of plastic strain4. The parameter K is defined by the expression
i.'j - -1 &j 1 --
+- 2Y s + yO(@(F))-
&j
aF I
2P E a0,
where the symbol (@(F))is defined as follows:
The function @(F) may be chosen to represent the results of tests on the
behavio~rof metals under dynamic loading.
It will be much more convenient to write the constitutive equations
(2.46) in the slightly different form
(2.47)
(2.48)
(2.49)
[Subsequent dynamic !w
The change of the yield surface during the deformation process is caused
by isotropic and anisotropic work-hardening effects, and by the influence of
the strain-rate effect.
(2.61)
where the function f(uH)depends now on the state of stress only. Thus
(2.62) f h j ) = fWiJwJa),
where J< denotes the first invariant of the stress tensor udf;Ja and Js are
the second and third invariants of the stress deviator sii. respectively. We can
write
276 PIOTR PERZYNA
where
(2.64) tij= sighj - #Ja &j.
(2.66)
If we then assume the Huber-Mises yield condition, i.e., that the function
f(c7q) = (Ja)"', we obtain by (2.47)*:
(2.69)
where
(2.62) u = #?) [+
1 m-1(;)] ,
where u = +(@) represents the static stress-strain relation for simple tension.
(2.63) F=--W J ~ )
1, c = const.
C
278 PIOTR PERZYNA
By assuming
(2.66)
(41.
and the dynamical yield criterion
The dependence of
Qi in Fig, 20.
n9 l/v
on as given by (2.68) is plotted for some function
It can be seen from (2.68) that the constitutive equations (2.67) lead
to a similar result as the inviscid theory of plasticity for isotropic work-
hardening material. But in the inviscid theory of plasticity for isotropic
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 279
where R, = Kv? denotes the radius of the static yield cylinder (Fig. 21).
In the inviscid theory of plasticity for isotropic work-hardening material
we have three possibilities according to whether Jg > O (loading), J2 = 0
(neutral loading), or j g < 0 (unloading).
FIG. 21. The cylindrical yield locus for elasticlvisco-(perfectly plastic) material.
“t
i
I f 1
c-const
c
&
(2.72) /(a&) =K or F = 0.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 281
(2.73)
(2.76)
(2.76)
I ae,acr,,
In the special case of elastic/visco-(perfectly plastic) material we obtain
by putting the result (2.68) into the relations (2.67)
(2.77)
(2.78)
of the incompressible perfectly plastic material [146]. These constitutive
equations hold only when the rate of deformation does not vanish.
B . A$plicatiort to Soils. As a particular case of the constitutive equations
(2.66), we shall study the elastic/visco-(perfectly plastic) soil for the following
static yield function*:
(2.79) F=
aJi + J;'' - 1,
k
where a is a constant describing the dilatation rate of the soil.
* For the discussion of the other cases see papers [127, 1281.
282 PIOTR PERZYNA
(2.81)
In the limiting case, y + oo, we obtain from (2.80) the known consti-
tutive equations for an elastic-perfectly plastic soil (see D. C. Drucker and
W. Prager [201]):
(2.83)
where
where the stress tensors u$) and ui:) and the corresponding strain-rate
tensors ij;) and St,)refer to two distinct mechanical paths. The first integral
is extended over a finite interval of time 0 - t - I,tk, the second term is the
value of the volume integral at tk minus the value at the instant of the
divergence of paths, t = 0. Necessarily, in all cases the velocity &@) through-
out the volume is the same as @) at t = 0,hence the second term of (2.88)
is positive or zero at all times tk. Therefore, the inequality
is much more significant (compare with (2.19)). It is obvious that the in-
clusion of dynamic terms does not affect the conclusions concerning the
stress-strain relations to be drawn as a result of neglecting these terms.
From the inequality (2.87),we reobtain the condition of orthogonality
of the inelastic strain-rate vector to the yield surface and that of convexity
of the subsequent dynamic loading surface. Both were assumed as basis
for introducing the general constitutive equations (2.47).
(2.89)
S=d/E+y*
[exp:(-- l)-lID
which, in slightly different form, was introduced by Malvern [llO].
204 PIOTR PERZYNA
We shall examine the types of function @(F) (2.88) in the light of exper-
imental data. At first we shall try to determine the constants y*, d, A,*
and B,* (a= 1,2,. ..,N) according to the results of D. S. Clark and P. E.
Duwez [39], so as to describe best the dependence of yield stress upon
strain rate, but we must remember that the strain rate in impact problems
accompanied by the wave propagation effect may change from zero to about
lo00 sec-l.
__-----
FIG.24. Comparison of the experimental data with the prediction of the power strain-
rate law.
0’ ‘ 4D’ ’ ’ I I I I I I 1 1 1 -
80I ’
a0 I
.(80 pw MI PBI) am
stmin mte (sBc.-’)
FIG.26. Comparison of the experimental data with the prediction of the linear strain-
rate law.
FIG.26. Comparison of the experimental data with the prediction of the exponential
strain-rate law.
286 PIOTR PERZYNA
The analogous comparisons with the linear law (2.91,) and the exponential
law (2.91,) are made in Fig. 25 for y* = 180-300 sec-l and in Fig. 26 for
y* = 80-240 sec-l.
To determine theconstants A,* and B,* (a = 1,2,. . ., N)wemay take the
data of D. S. Clark and P. E. Duwez [39], given L. as a function of a/ao- 1.
Hence we can write two systems of linear equations for A,* and B,*
..,N,
N
2 A,*[eXp(Xi)a - 11 = B ( X J .
a-1
i = 1,2,,
N
where Xi = ai/u0- 1 for some points of Q = o,,u,,. .., QN, and d(X,)
.
(i = 1,2,. .,N) are known from experimental data.
I
min mte(kc:')
m
I I
m s
I #
a '
m
I -
FIG.27. Approximation of the experimental curve by the functions (2.88,) and (2.88,).
it as a hypothesis. This hypothesis will permit the use of the same constants
in the general constitutive equations (e.g. in (2.68) or in (2.67)) as those
determined for one-dimensional states of stress.
TABLE2
7lme (psec.)
FI 28. Comparison of the experimental data for strain rate with the prediction of
the power law.
6=6/E+y* exp - -
o: ( J J
(2.93) N r ,
Time (psec.)
FIG.29. Comparison of the experimental data for strain rate with the predictions of
the exponential law.
The best agreement with experimental data can be observed for the power
function (P(F), (2.93,) and for the exponential function @(F), (2.93,). The
results for these functions are plotted in Figs. 28 and 29, respectively
(cf. [133]).
All theoretical curves show that in comparison with experimental data
the maximum value for the strain rate is shifted in time by about 15 psec,
which may be due to the delay time. Indeed these shifts have the order of
magnitude of the delay times observed in mild steel under stress of up to
100,OOO Ib/ins (see J. E. Johnson, D. S. Wood,and D. S. Clark [87]).
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 289
8
*
P
h,
9 6
9
4 4
0 P 4 6 8
Rillin(%)
FIG. 30. Comparison between the experimental and theoretical stress-strain curves
for mild steel.
0 4 6 8
Strain (“/o)
.FIG.31. Comparison between the experimental and theoretical stress-strain curves
for pure iron.
Consider now the inelastic part of the strain-rate tensor. For perfectly
plastic materials with Huber-Mises yield condition (2.94)gives
or equivalently
(2.97)
(2.98)
Since all experiments of interest were carried out under tension or compression,
(2.97)and (2.98)will be useful in order to give experimental evidence of
the present approach.
292 PIOTR PERZYNA
The most general case in which the quantities y and a, and the function
@ are temperature-variable has not been studied in the literature. The main
reason for this is that the next two cases considered, where only two quan-
tities vary with temperature, were found to give satisfactory correspondence
with experimental data over the entire range of temperature. Let us con-
sider in detail the following three particular cases:
A. a, = a,(@ and y = ~ ( 6and ) @ is independent of temperature. It was
shown that the viscosity of the plastic material proposed by Hohenemser
and Prager can be used to describe the strain-rate sensitivity of metals.
Since the viscosity constant y , similarly as the yield point a,, depends on
FIG.33. Compariaon of the theoretical predictions with the experimental data for
mild steel.
m-
m-
80-
go-
40-
Pure imn
1oQ-
Q (~~=4saoo-~ao(e-m)1~z
4?l-
8aT
80- "'~i77/45?1-Sl(s-l.R)
'\ n-10
FIG.36. Comparisonof the theoretical predictions with the experimental data for pureiron.
294 PIOTR PERZYNA
t .WmJ
Pure imn
arm-%
*s -
25-
Id-
29 -
22 -
w-
40 -
0.9-
08-
FIG. 37. Comparison of the theoretical predictions with the experimental data for
aluminium.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 296
(2.101)
The results of tests on aluminium, iron, and steel taken from [l, 34,
96, 109, 1121 and [149], have been elaborated for our purpose. Only the
range of temperature is considered, in which the material is metallurgically
“stable.” In order to use the results of tests the form of the function Q, has
to be assumed. The power function has been found to be preferable (see
(2.88,)), and the introduction of more involved function is unnecessary,
since the power function @ gives a very good agreement with all tests con-
sidered.
Circles and triangles in Figs. 33-37 are experimental points both from
static and dynamic tests. All static tests of u,, were approximated by either
power or exponential functions and the lower solid lines are plots of u,(0).
The upper solid lines represent u(i,0) for one or several values of strain
rate. It is seen, that by a suitable choice of y ( 0 ) very good agreement with
experimental data can be achieved. This is evidence of the correctness of
the formula (2.96),case A. In the same figures the dash and dash-dot lines
correspond to the hypotheses of Vishman, Zlatin, and Joffe [186] and Volo-
shenko [187], respectively. In both cases a definite deviation from the
experimental points is evident.
It should be pointed out here that the comparison with experiments has
been carried out under the assumption that the material is perfectly plastic,
hence the strain-hardening effect is disregarded according to (2.96). If we
include strain-hardening, (2.94) yields
(2.102)
or
(2.103)
for a one-dimensional state of stress. Eqs. (2.102) and (2.103) could not be
compared with experimental data because of lack of suitable tests. On the
other hand, an excellent agreement of (2.96) with all tests here presented by
no means denies the existence of strain-hardening properties when the rate
of strain is high. This is because of the arbitrariness in the choice of the
function y ( 0 ) which was taken to fit the experimental data. If we had in-
cluded the strain-hardening in the analysis, we would also have changed
the function y ( 0 ) . I t seems most probable that (2.102) and (2.103) might
give a better description of the dynamic response of certain metals in larger
ranges of variables.
Recent investigations indicate the strong interaction between yield
stress, strain rate and temperature in dynamic tests of various metals.
The present approach shows how a simple temperature modification of the
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 297
12. Linearization
(2.104)
(2.107)
(2.109)
(2.110)
s
During an A-process the tensor rijvanishes. Thus,
(2.111) +
ak(0ikvf 0jhvi)dV = (uikakvj + ujk&vi)dV
'5
V
=
2
(0&j + 0jkkki)dV = 0.
V
The states of stress and strain are homogeneous during this type of relaxation
test ; we therefore have
(2.112) i(u&' + a&i) = 0.
the relaxation process follows. Then from (2.67) and (2.112) we obtain a
system of six differential equations with respect to sij and akk
(2.115) $7 i k E'kr - 0.
u-a.
Consider now the simpler boundary value problem in which the loading
process is characterized by the surface tractions T i being prescribed on the
entire surface S, while during a relaxation process the surface velocities v j
vanish on the entire surface S. This kind of relaxation process we shall call
a B-process.
Because the states of stress and strain are homogeneous during a relaxa-
tion test, we now have conditions
(2.116) d..$7 -
- 0.
From (2.67) and (2.116) we obtain a system of five differential equations
for sii in the form
(2.117)
(2.118)
and it is worth noting that these conditions are valid for A-processes in
incompressible materials because then
(2.119) in = 0.
(2.120)
Jd0)+ nC
W
(2.124) L l ~ n + l ) ( t )- J ~ ( 0 l
=O
is the solution of the integral equation (2.121), and hence of the differential
equation (2.120). Of course, the solution (2.126) is valid only in the nonelastic
region J s >ka.
Full discussions of the solutions of the relaxation equations have been
presented in the papers [132, 128, 891.
t a
t b
Strain m@ d
c L
Temperature 6
t C
L -
strain
FIG. 98. Increasing strain rate i,or decreasing temperature 0 promotes fracture
(D. C. Drucker [SS]).
1. M a t h m t i c d Preliminaries
+
The operator a/at 1,) (a/az) in the ith equation represents differentiation
along the ith characteristic curve determined by
(3.4)
(3.5)
The solution of the initial value problern for the system (3.1) has been
studied by numerous authors. First existence and uniqueness were discussed
in full (see K. 0. Friedrichs [70], R. Courant and P. Lax [511, P. Hartman
and A. Wintner [76], P. Lax [103], A. Douglis [61] and R. Courant [54]).
Some theorems regarding existence and uniqueness may be obtained as
special cases of the same theorems proved for the system of first order partial
differential hyperbolic equations involving n independent variables (see
J. Schauder [154, 1551, M. Cinquini-Cibrario [36, 371 and R. Courant [54]).
In practice, to obtain the solution of the initial and boundary-value
problems for the system (3.1) the method of finite differences is used (cf.
R. Courant and K. 0. Friedrichs [50]). For a study of convergence of this
method and the conditions of its application see R. Courant, W. Isaacson and
M. Rees [52] and G. Prouse [l48].
Some modifications of the finite difference method, and its application
to the solution of the initial and boundary-value problems for quasi-linear
and semi-linear systems (3.1)were proposed by H. B. Keller and V. ThomCe
[93, 94, 174-1761, R. Courant [53, 541 investigated the possibility of ap-
plication of successive approximations to the solution of initial and boundary-
value problems for the system (3.1).
.xb?J(x)) =g ( x ~ . u ( x ~ q , ( x ) ) ~ ~ Y ( x , q ( ~ ) ) ) ,
Theorem 1 . If
1
I / ( ~ , Y # W X , % )- f ( X , y , ~ , ~ z , U 5 - UI
y ) L( (4 + 14, - u,[ + lay- zcyl).
(3.8) Ig(X,Wy) - g(x,u,u,)l 5 LI4 - + M ( 4 , - uyl,
(h(y,n,cx)- h(y,u,ux)I 5 L ( 4 - 211 + N(4, - 21x1 ;
where X, = max(xo,yo);
there exists a unique solution u(x,y) E CQ1 to the boundary value problem
(3.6), (3.7).
Proof. We shall discuss here only the basic elements of the proof,* which
we shall use to construct the solution of the nonlinear problem (3.6)-(3.7)
by means of the iteration method.
Consider the linear space CQ' of the functions u ( x , y ) . The norm in this
space we shall define in the following form:
(3.10) l SUPINXJJl
l l ~ ( X J J ) l=
Q
+ SUPl%(XJJ)
Q
I+ IUY(X>Y)I.
(3.11) fi(X,Y) =9 [ U ( X , Y ) I ,
where
W[.(x,y)l = *. + i
X*
g(s,.(s,~(s)),.y(s,~(S)))ds
(3.12) +1
Y
Y.
~(Zl.(~(Z),Z),.X(~(Z),Z))dZ + 5I
Y X'
Y* P(.,
~(s,z,.(s,z~,~,(s,z),.y(s,z))dsd%
+ j j~~s,z,.~s,z~,.,~s,z~,.~~s,
X* CPW
(3.16)
Note also that by assuming the Lipschitz condition (3.8) in the weaker form
5 -4
I f ( x , Y , u * M y )- f ( x , Y , ~ , w y ) lL((4, + (21, - uyl),
3.18) 1 5 Mp, - uyl,
Ig(x,%cY) - g ( x , w y )
/h(Y,%%)- h ( Y , W %I )5 Np, - %I
and using the Schauder theorem we may prove the existence of the solution
for our boundary-value problem.*
Consider now the special case of the boundary-value problem (3.7),
where we assume
(3.21)
where a,, q, Pi E CQO for i = 1,2. The same problem as (3.6), (3.22) was
posed by J. Conlan [46] who solved it by a polygon method.*
The mapping for the problem (3.6), (3.22) takes the form
(3.23)
* The Euler-Cauchy polygon method for proving the existence of a solution for
It has been shown (see [166-168, 991) that all classical problems (i.e. the
Cauchy, Darboux, Picard and Goursat problems) can be analogically treated
as special cases of the nonlinear problem (3.6)-(3.7). In such a procedure,
particular note must be taken of the restrictions on the functions which
describe the boundary conditions.
Note that all the results discussed here are valid for a system of second-
order semi-linear hyperbolic partial differential equations.*
C, A non-linear boundary-value problem for a linear hyperbolic equation.
Consider the normal form of the second order linear partial differential
equation of hyperbolic type in two independent variables :
(3.24) s(u) uxy + a ( x , ~ ) u+x +
b ( x ~ ~ ) U yC ( x , Y ) U = d(x,y).
+ + lWq)Vx*,q;x,y)dq,
(3.26)
where
4x,y) =S(x,y)
X*j Q(E)W,y*;x,y)@
Y*
where
310 PIOTR PERZYNA
KO = [IK~(~,x)I.IK~(~,~)I,IK~(€,~)~.~K~(~~Y)II.
8
(3.29) i = 1,3,6; i = 2,4,6; K = 7,9,11; 1 = 8,10,12;
(
2K0X0 1 + -
M*&")- + N* < 1;
there exists, throughout Q, a unique solution w(x,y) E CQ1to the boundary-
value problem (3.24),(3.7).
Proof. We shall show by successive approximations that the system
(3.27)possesses a unique continuous solution set a(%)
and n(y).
We define
the sequences {Q(,,](x)} and {II(,,)(x)} as follows
5 o(n
X rp(Z)
~(rn,(x)= ~ ( x ) + -~ ) ( E ) K ~ ( E , w + -*)(q)Ka(q~dq
X* Y*
1
x
(3.32)
The inequality
(3.36)
* That the problem (3.24), (3.7) has a unique solution can be shown easily.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 313
(3.37)
q@)
= d(x,y),
(3.39)
aO(Wx,v(x)) + .l(x)%(x,V(x)) + a,(x)Uy(X,v(XN = %(X)J
P ~ ( Y ) N $ ( Y ) A+ P i ( ~ ) u x ( $ @ )+
A P a @ ) % y ( # ( y ) , ~=
) us@),
zl(x*,y*) = u*.
The solution of the problem (3.39)has been given by C. Chu and J. B.
Diaz [36]. The case (x*,y*) = (0,O). was treated in the paper [S]by A. K.
Aziz and J. B. Diaz.
Using the solution (3.2)-(3.3).we can reduce the boundary-value problem
(3.16)to an equivalent system of linear functional-integral equations :
X 9(*)
J
X*
(3.40)
It is obvious that the problem (3.39) can be treated as a linearized form of the
problem (3.241, (3.7).
314 PIOTR PERZYNA
where
Theorem 2b. If
1~2(X*)/.l(x*)l -= 1, JBlb*)IP&.Y*)l<
1;
6. Y ( x * ) = y*, $@*) = x * ;
N= "0" [ l ~ l ( e ~ ~ ; ~ ) ~ , ~ ~ * ( e l ~ ; y ) ~ l ~
there exists, throughout Q,a unique solution rs(x,y) E CQ1to the boundary-
value problem (3.39).
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 316
Q(s)(x) =G*(x)
P(4
i
+ A ( ~ ) n ( f i - ~ ) ( p (+x ) )klO,r*;x)n(s-l)(5')dE
I.
J
Y*
(3.42)
for n = 1.2,3,. . ., and for arbitrary continuous functions Q(,(x) and n(&y).
Proceeding as in the proof of Theorem 2, we easily obtain from (3.42):
(3.44)
Theorem 2b is valid (see C , S. Chu and J. B. Diaz [36]). In that case, the
constants N and X, should satisfy the following inequality: SI?X,< 1.
G. Majcher [lo81 has considered the boundary-value problem consisting
of the partial differential equation (3.24) subject to the boundary conditions
(linear generalized F'icard problem)
,
u(x,O) = a&),
(3.45)
BoCr)MlbCr)A + A(r)4lb(r)d + A ( r ) N v ( l b b ) A = OaCy).
A. K.Aziz and J. B. Diaz [6] have shown that the problem (3.46)stated
by G. Majcher [lo81 may be treated as a particular case of the linear problem
(3.39).
316 PIOTR PERZY NA
Full discussion of the other particular cases of the linear problem (3.39)
has been presented in paper [a], where the history of the linear problem
and an extensive bibliography is also given.
All theorems discussed in sections B and C are valid “in the small.”
The inequalities (3.9), (3.34), (3.37), and (3.43) required to prove the existence
and uniqueness of the solutions are very strong and impose important rest-
rictions on the size of the rectangle Q and on the boundary conditions.
(3.46) = %(I$), U, = Ng = 0,
where u,,, uq, Ue are the spherical components of the displacement vector.
Assumption (3.46)(spherical symmetry) yields the following components
of the strain and stress tensors:
au
- a4
(3.47) &" =
C ' Epnp=Eee=-' I
0 --
l o o 0
P
-(K+Qp) 0 0 0 0
U= , A=
#ru - K 0 0 0 0
-1 0 0 0 0
0 0 0 0 0
(3.49)
B=
--
V
I
J
T h e characteristics (3.4) have the form
3p )
(3.60) I = const, I = I , f2 + const, where A = (
4p + 3 K 'Ie
dE,
1
- -au,,
2P
+-1
2P
au,, - (+ + vy.) at = 0,
(3.61)
+
- ( 4 ~ 3K)dv f3Ad~pp= 0,
respectively.
If a pressure exceeding the plasticity limit, p > p , is suddenly applied to
the surface of the cavity, then the characteristic line Y - Y, - & = 0 will
be a strong discontinuity. Along the discontinuity additional conditions
must be satisfied, which will be called the conditions of kinematic and
dynamic continuity. For the spherical discontinuity these conditions have
the form (see for instance G. Hopkins [82]):
(3.62)
(3.63)
respectively.
Since the strain rate 6 , at the discontinuity must be regarded as infinite,
the constitutive equations give the relation (see [llS])
(3.54)
da,
--
(3.66)
dr
- - qY,U")*
with the condition
(3.66)
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 319
where
(3.67)
To solve Eq. (3.66) one has to show that on the front of the shock Y = ro + At
the hardening parameter K = K(W,) appearing in (3.57) depends only on the
variables Y and u,,. By definition (2.44) we have for spherical waves
(3.68)
J J
0 0
(3.69) [Ew] =
1: 1= o ;
-
here [ ] denotes the jump of the quantity in brackets across the discontinuity.
Separating the component of the strain tensor E,, into the elastic part E:,
and the plastic part f and using (3.52) and (3.53) we can write
(3.61) 4 = uw/(E.ap)s
and after a straightforward calculation
(3.62) w,= 0.
Since on the front of the wave the plastic work equals zero, the problem
of the work-hardening plastic material is reduced to the perfectly plastic
material, the hardening parameter K being equal to the yield stress k.
In the case of cylindrical waves and plane waves the proof that along
the shock K(WJ = k is fully analogical and will be omitted.
Equation (3.65) with the condition (3.56) leads to the nonlinear Volterra
integral equation of the second kind
(3.63) 0,
1
= Po - ‘Y[E,u,(E)ldS.
VI
where Y(Y,U,,)
has now the form
(3.64)
320 PIOTR PERZYNA
FIG.39. Elastic and plastic regions in the t,r-plane for the case of discontinuous front
wave.
and
v(n) = - &)/@A),
(3.67)
= a!?/@P);
&$I Ew = 0.
The solution given by (3.66)-(3.67) is valid for Y < Y* (Fig. 39) where I*
satisfies the condition
(3.68) V
m= k.
The solution along the discontinuity for r >r* has the closed form
1 Y*
a,, = -
R*
-
Y '
v=--- 1
plR*
r*
1 '
(1-- ;;);*:I --
(3.69)
r*
e, = @ASR*)-' ;
, eW = 0,
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 32 1
where
(3.70) R* = 2p@Aak1/37-1.
ConsLucrnow the solution in the elastic/viscoplastic region, ..e., region
P+ in Fig. 39 which is bounded by the discontinuity r - ro - At = 0 and
the straight line I = yo. On the discontinuity the quantities b,, ow, u, E,,
and e, are determined by means of the solution (3.66)and (3.67), while
a,, is known on the line I = r,, from the boundary condition.
Denoting the intersection points of the characteristic net as shown in
Fig. 41 we approximate Eqs. (3.61) by the following difference equations:
322 PIOTR PERZYNA
(3.71)
+ 3K)
-( 4 ~ [v(i,m,N) -~ ( -1 l.m,n +1)I - 31[~r(l,m,n)
- c,p(I - 1,mn,n+1)1= 0.
In every point of the characteristic net we have a system of five algebraic
equations with respect to five unknowns a,,, aw, v , E,, and
The procedure of finite differences can be applied to the elastic/visco-
plastic region in the case shown in Fig. 40. In this case we know the quantities
a,,, aPp,v , err and ew at the point t = t, on the line r = Y, from the elastic
solution in the region E*.
vm
The boundary *'I of the elastic/viscoplastic region (see Figs. 39 and 40)
can be obtained by an approximate method using the condition = K(Y).
0 - AP 0 0 0
0 0 0 0
U= A= 0 0 0 0 ,
0 0 0 0
0 0 0 0
0 0 0 1
(3.76)
B=
- -V
r
where
a&,, = -
V
7
at,
(3.78)
1
derr- -(da,,
3K
+ da,, + da,,) + -dt = 0,
V
7
&-
33, or,- aw
P r 1
dr = 0.
(3.80) YJ(r,fJ,?)
=- 1-
a?,
2 7
+ 31 0 (m-
2wr
1).
and
The solution on the cylindrical discontinuity for r > Y* has the following
closed form:
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 326
z(;[ - y)+ y
q 2y@ - 1) - -
31 a,,
c1
- a,
r 1
(I-l.m.n+l)
Ar
(3.84) u, = @ ( I $ ) , u, = u, = 0.
(3.86)
(fr.
The characteristic lines (3.4) are
(3.91)
FUNDAMENTAL PROBLEMS IN VISCOPLASTLCITY 327
(3.03)
and
The strain tensor has only one component that is not identically zero, namely
(3.99) Eli = &aa(f,t),
U , A=
. I
-1 0 0
r
(3.100)
0
1
B = p 5 y @ [ 3 ( u 3 3 2- VK5 eK i 3 ) - l
(3.106) (ILL = Po
j
- v[S,aLa(S)l&
0
where
and
(3.107)
Using iterations we may write the solution of Eq. (3.105) in the form
where
v(n)
1 (*)
= -- (4- 1 (n)
(3.110) a**I ELL- __ U L L .
PAa
The condition for f * (see Fig. 39) is
(3.111) VJa(x*) = k-
The solution on the plane discontinuous wave for f > 2* is
330 PIOTR PERZYNA
(3.116) t = .r&.t).
- 13.
CascA.*Non-homogeneousmatmial; dinearficnctionQ)=@[(r/k(r))
It will be assumed that all functions describing the mechanical properties
of the material vary with the radius Y only. We then obtain the semi-linear
system of differential equations
(3.116) Uc + AU, + B = 0,
where
(3.117)
(3.118)
(3.119) h ( d = f [ru(*)/p(~)l”e
are real. If shearing tractions p(t) or shear strain rate w(t) exceeding the
plasticity limit, p , >$, or a,,>a,,,are suddenly applied to the surface
The problems which will be studied in the Case A are much more general than those
discussed in Refs. [ISO, 131, lS5].
332 PIOTR PERZYNA
I
r
where
(3.122)
and
(3.123)
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 333
The elastic/viscoplastic region 9*at Fig. 42 now takes the form of the
region 9 at Fig. 43. The region 9 is bounded by the characteristic y = 0
and the line x = y (not a characteristic).
The boundary conditions (3.118) now become
where
rY(n+l)(X,Y) = WY [.ccn,(XJ41*
(3.130)
We seek the solution of (3.129), which satisfies the following boundary
conditions
(3.131)
Thus the problem in the region 9 (see Fig. 43) for the linear function @
is reduced to the solution of the special case of the linear problem (3.39).
(3.133) b ( w )=-
3.0
2 { P o + &(x - r)l-l+ &Yofo/Q.
C ( X , y ) = 4&'[2y0 -k &(X - y)]-'.
The boundary conditions have a form similar to that of the previous case.
338 PIOTR PERZYNA
for spherical waves, cylindrical radial waves, and cylindrical shear waves,
and new coordinates
(3.136)
(where x,, = const) for plane waves, the equation of motion in the elastic
region may be written in the form
A
(3.138) Wxy = ~
x- Y
(ux - M y ) + ( x -6y)2 @*
The considered elastic regions in the x,y-plane are shown in Fig. 44.
Equation (3.136) is valid in the regions El and E,.
Consider first the problem in the region El. This region is bounded by the
characteristic y = 0 and the curve x = ~ ( y ) .On y = 0, we have u = 0 ;
on x = q(y) in the cases of spherical waves, shear waves, and plane waves
the linear condition
(3.137) ci(%x -~ y + caw =
) %s
where
must be satisfied.
Hence for all four types of waves the problem in region Elreduces to the
generalized Picard problem (see G. Majcher, [lOS]).
The solution of this generalized Picard problem may be expressed as
follows
(3.139)
i
U(%Y) = W , Y ;xo,q)Q(r])dr],
0
(3.140)
0
(3.143)
338 PIOTR PERZYNA
+c
oc
(3.144) RCy4) = J C r ( Y B l 7 ) n i l~ ( & J d D
with
Y
Thus in the case of the nonlinear condition (3.138). Eq. (3.140) is equivalent
to the equation
The Riemann function for Eq. (3.136) has the following closed form (see
R. Courant [64])
(3.149) Y(x,y;t,q)= ( x - Y)”(t - r)%- q)JF(- A- B;l,t),
where F(- 8,- B;l,() is the hypergeometric function and is determined
by the relation
(3.160)
and
(3.161)
The values of the constants A , 6, LO and p and the coefficients c,, cg and cg
are presented in Table 3. The case of a plane wave can be treated separately.
In this case the Riemann function V(x,y;{,q)= 1 and the solution is trivial.
In the region E, (Fig. 44) the problem for all four types of waves
may be reduced to a similar generalized Picard problem. But in this region
it is useful to apply the well-known Fourier transform method (see A. Kromm
[98] and G. Hopkins [82]).
6. Numerical Examples
We shall discuss here the numerical examples presented in the paper [136].
As a first example, the propagation of plane waves in the half space is
considered. In this case the geometric dispersion does not have such a great
influence on the solution as in the other cases. In the spherical problem the
convergence of the method of successive approximations, due to the geo-
metric dispersion, is so quick that the difference between the first and the
second iteration is practically insignificant, whereas in the case of plane
waves the first and the second iteration and for certain small regions of 2
even the next iterations can lead to different results (see curves 1-6 of
Fig. 46). Thus i t seems that in the case of plane waves the influence of the
work-hardening of the material may be studied more effectively.
For simplicity of the equations in practical applications, linear and
power-law forms of the function Q, and linear work-hardening of the material
have been assumed. Thus Q, has the form (cf. (2.88,))
(3.162)
FUNDAMENTAL PROBLEMS I N VISCOPLASTICITY 341
0 4 6 8
and
(b) for a power-law form of qZ
(3.166)
342 PIOTR PERZYNA
TABLE4
Young's Tangent
2.1 * l@k G cm-a 0.014062* loe k G cm-s
modulus E modulus
Poisson's
0.29 y1 = 45Osec-'
ratio u
Viscosity yp = 680 sec-'
Bulk constants y3 = 736sec-'
1.663- 106 kG cm-* y, = 7OOsec-I
modulus K
Shear
0.820- l@kG cm-*
modulus p
-4
I I I I
Vm
'OO 2 4 6 8 1 0 2 -
, Curves
h ~46. versus I . Comparison of the exact linear solution with the exact
power solution.
In Fig. 46 the curves V x / k versus L at the discontinuity are given. Full lines
represent the results of the linear solution, broken lines correspond to the
power-law solution. Material constants y1 and ys have been assumed here
the same as in the case of a perfectly plastic material.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 343
The object of the second example is to find the curve (see Fig. 39) dividing
the plastic and elastic region. The influence of such parameters as duration
time of loading, the shape and sign of the loading curve, constant of material
y and work-hardening parameter K on the character of the curve P* is also
discussed. The comparison of the results for work-hardening and perfectly
plastic material has been presented.
Computing the third iteration by (3.65) and (3.66),we obtain the fol-
lowing result a t the discontinuity:
(9) 17
o,, = nya(R*$o)
- i@y2+ $,(l - D )
{ + 150 + 18)
- A?y,R*Po [(R*p,)9(D3+ t3Dz
- D [(R*$o)z(D' + 3 0 + - R*p0(2D + 3) + 1 rz
1
+ R*p,,Dz [R*po (D +):
I
- I] r3 -
1
(R*@,)'D44}
--
2
1
R*$,(log .Iz} I* + L ' @ ~ , ( R * $ ~ ) ~logD % ~
1,
(3.168)
where
(3.167) D=- *"
A
(R*$, - 1)8.
344 PIOTR PERZYNA
Eqs. (3.67) are employed for the determination of the quantities a*. v ,
E,,and eCC at the discontinuity.
Condition (3.67), according to (3.152) and (3.153) together with (3.166),
may now be written as follows:
(3.168) R*&)(r*) = 1.
L
tm m
I
.uo
I I
*w
I I I I L-
.1m PaJ rh
FIG. 47. Curves a,, versus v / r o on the discontinuous wave.
In Fig. 47 we have two curves a,, versus r at the discontinuity for applied initial
pressures p , of 5800 Kg/cma and 3740 Kg/cma, all material data being assumed
as for the plane wave. The character of the curve r*/ro versus V x ( y 0 ) / k
FUNDAMENTAL PROBLEMS I N VISCOPLASTICITY 345
and
(3.162)
0
(3.163)
(3.164) 6; +2 4 = 0,
from (3.162), we obtain
1
a m - a,,
(3.165)
0
With the notation P(l,m,n) for the intersection point of the characteristic
net as in Fig. 41, the following recurrence formula for determination of the
parameter K for any fixed no has been employed:
346 PIOTR PERZYNA
49. Influence of the time duration on the character of boundary r*for linear form
FIG.
of preeeure.
Fx0.60. Influence of the time duration on the chsracter of boundary P for power form
of pressure.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 347
This way of obtaining the sth iteration has been chosen in the numerical
computations in order to satisfy the condition I K ( ~ + ' ) - &)I < 1. The
absolute value of an intensity a,, - oqqhas been taken because W , from the
definition of the plastic work, should be positive. The following pressure
functions have been assumed
FIG.51. Influence of the viscosity constant y on the shape of the curve r*.
In Figs. 49-51 and in Table 1 the results are presented. In Figs. 49 and 50
the curves r* for linear and quadratic form of P(r) and for various loading
periods have been plotted. I t is seen that if the period increases, the character
ofr* changes. The influence of the material constant y on the shape of the
r*
curve may be observed in Fig. 51. The difference is seen only at the
shock and in points lying close to it. From the results collected in Table 6
the small influence of the work-hardening effect is seen. However, it should
be taken into consideration that the hardening parameter K depends first
of all on the duration of the process of plastic deformation and on the amount
TABLE5 w
%
AV = 0.04 AY = 0.01
t = 1.34576- 10-7
-4833.33 - 1406.65 - 1407.28 - 1436.07
1.00
1.04 - 4635.99 - 1382.42
-4833.33
-4635.80 - 1382.66
-4833.33
-4642.19
- 1435.70
- 1409.03
-4833.33
-4642.05 - 1409.18 z
1.08 -4451.33 - 1354.87 -4451.05 - 1354.92 -4463.17 - 1379.60 -4462.98 -11379.63 8
1.12 -4278.79 - 1325.46 -4278.49 - 1325.41 -4295.68 - 1348.64 - 4295.47 - 1348.61
1.10 -4117.69 - 1295.19 -4117.38 - 1295.10 -4139.03 -1317.03 -4138.80 -11316.97 8
1.20 -3967.28 - 1264.78 -3966.96 - 1264.67 -3992.47 - 1285.38 -3992.24 - 1285.31
1.24 -3826.77 - 1234.71 -3826.46 - 1234.59 -3855.27 - 1254.12 - 3865.05 - 1254.04
1.28 -3695.41 - 1205.31 -3695.11 - 1205.19 -3726.73 - 1223.54 -3726.51 -1223.46
1.32 - 3572.49 - 1176.82 -3572.20 - 1176.70 -3606.17 - 1193.83 -3605.96 - 1193.76
1.36 - 3457.33 - 1149.37 - 3457.05 - 1149.25 -3492.97 -1165.12 - 3492.76 - 1165.04
1.40 -3349.31 - 1123.06 -3349.03 -1122.95 -3386.53 - 1137.47 -3386.34 - 1137.40
1.44 -3247.84 - 1097.95 -3247.57 - 1097.85 - 3286.33 - 11 10.92 - 3286.14 -1110.85
1.48 -3152.40 - 1074.07 -3152.14 - 1073.96 -3191.79 - 1085.45 -3191.60 - 1085.38
1.52 -3061.72 - 1051.10 -3061.41 -1051.00 -3102.50 - 1060.27 -3102.32 - 1060.20
t = 2.69152- 1 0 - 7
1.00 -3866.67 -541.45 -3866.67 -543.27 -3866.67 -604.02 -3866.67 -605.10
1.04 - 3691.25 - 575.89 -3690.98 - 576.67 -3698.83 - 633.44 - 3698.66 -633.90
1.08 -3527.76 -598.05 - 3527.39 - 598.34 - 3542.21 -651.28 - 3541.97 -651.41
TABLE5 (continued)
Ar = 0.04 AY = 0.01
-
- Po, perfectly - p,, work-hardening - Po. perfectly - pc work-hardening
plastic material material plastic material material
Qw Qw UY, aw QW =w Qw
IV. QUASI-STATIC
SOLUTIONS
1. Spherical Problem
a ~ ,
--+2
-I oprp=
a 8
(4.2) t )pH(&
~ , ~ ( a ,= u,,(b,t) = 0,
where a and b denote the inner and outer radii respectively and H ( t ) is the
Heaviside function.
[(:r-
The solution of the system (4.1) by means of Laplace transform leads
to the following formulae
y a k=
F - +(iy 11 -'[- p + 2k(l - e ~ p ( -m*t)) log-a
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 361
where
(4.6)
Quasi-static
01
FIG. 62. Curves vJu/k versus v/a for different time parameters (Ref. [194]).
The formulae for the components of the stress tensor a,,, uw and uee are
given in [191], where the influence of viscosity and time on the stress distri-
bution are also discussed. The interaction between elastic and inelastic
362 PIOTR PERZYNA
components of the stress tensor causes a certain levelling off in the distribution
of J 2 along the radius r , Fig. 52.
Since the deflections of the plate are supposed to be small a linear relation
between rate of curvature and rate of deflection is assumed:
&=--
a% 1 aw
(4.7) up=---
r8r
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 363
If the inertia effects are neglected the equation of equilibrium takes the
form
where F = V&/k - 1.
The principal idea of Eq. (4.9) is that the strain rate should be in general
a non-linear but uni-valued function of the excessive stresses above the
yield surface. Two particular types of the function @(F)will be considered,
namely, the linear function @(F)= F and the power function @(F)= Fd.
In the case of the power function the constitutive equations (4.9) give two
independent equations
(4.10)
(4.11)
The stresses a,, and aw can be evaluated from (4.10) and substituted into
(4.11). Taking account of the relations (4.6) the integration can be performed
for any integer exponent 6. This leads to the following formulae relating
rate of curvature k, and kp with the radial and circumferential bending
moments :
364 PIOTR PERZYNA
(4.12)
where the constant B is
(4.13) B = y/1/5h[(26 + 1)/26]d.
Note thatmo formal analom. between (4.10) and (4.12) exists. Equations
(4.12) involve a new term which vanishes only for the linear function, ie.,
for 6 = 1. This example shows that particular attention should be paid
in constructing a stress strain-rate relation appropriate for the generalized
quantities.
It is convenient to introduce dimensionless variables
M, M w
(4.14) m,=-M o '
Mo' P=B* v=- BR'
We have now five equations with five unknown functions m,, mvJv , k, and kv
However only three of them are of interest to us. These functions are m,,
m, I and v. After eliminationof k, and kpthe followingsystem of three ordinary
quasi-linear differential equations is obtained
(4.16)
FUNDAMENTAL PROBLEMS IN VIXOPLASTICITY 366
FIG.63. Curves m, and %versus i; for b = 1and several values of the loading paramebr
p’ (Ref. [193]).
358 PIOTR PERZYNA
24 t mrtmcp
6-3
FIG.
84. Curves m, and n"p, versus for d = 3 and several values of the loading parameter
p' (Ref. [l03]).
desirable to compare the moment and velocity distribution for the case of
non-linear and linear function @(F). This has been done in Figs. 67
and 68, where the solid line refers to the case 6 = 3 and the broken line
corresponds to the case. 6 = 1. Note that the moment distributions do not
c
P
FIG.56. The velocity fields for 6 = 1 and different values of the loading parameter p‘
(Ref. [193]).
FIG.68. The velocity fields for B =3 and different values of the loading parameter p’
(Ref. [193]).
TABLE6
P'
22-
10 -
ae-
06-
04-
02-
0
,
a2
, ,
0.4
I I
06
, I
ne
l l
1.0 F
-
FIG.57. Comparison between the distributions of the moments m, and % for non-linear
and linear function @(F) (Ref. 11981).
FIG.68. Comparison between the velocity distributions for non-linear and linear function
@(F) (Ref. [l93]).
tensor has a constant direction for each region while the others vanish. On
the other hand in the majority of cases for visco-plastic material there are
two or more non-vanishing components of the strain-rate tensor and the
direction of the strain-rate tensor is no longer constant. I t can be concluded
that, at least for simply supported circular plates, the linearized theory of
viscoplasticity due to W. Prager [146] may give a close qualitative and
quantitative assessment of the deformation under a given condition of
loading. Of course the above statement relates only to the linear function
@(F) for which the comparison was carried out.
The above procedure involves the assumption that the deflection of the
plate is small. It is straightforward to establish upper bounds on the values
of pressure and duration of impulse so as to remain within the limits of the
theory.
360 PIOTR PERZYNA
08 -
06
a4
a2 -
-
-
-Huber -Mises
---- Tmca
1 I
I
I
FIG.69. Comparison between the distributions of the moments m, and m, for the Huber-
Mises and Tresca yield functions (Ref. [193]).
Since the applied pressure p(t) was assumed to be constant the time
is formally eliminated from the solution. Although there are no restrictions
to extending of the solution over the range of variable pressure, a sufficiently
rapid time variation of the load intensities would require consideration of
inertia effects. The dynamics of a visco-plastic circular plate has been con-
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 361
sidered in Ref. [196]. The treatment makes the same assumptions as in the
quasi-static problem, but the transverse inertia motion is taken into ac-
count. It is shown that the dynamic flow of a viscoplastic circular plate
is described by an initial- and boundary-value problem for a quasi-linear
parabolic system.
W/BRp
FIG.60. Comparison between the velocity distributions for the Huber-Mises and Tresca
yield functions (Ref. [193]).
V. OTHERDYNAMICAL
PROBLEMS
a v u
r+2r=0.
We shall assume the power law @(F)= Fd [cf. (2.88)] and arbitrary
pressure +(t). The boundary and initial conditions have-now the fo&
(6.2) u,,(a,t) = - p(t). u,(b,t) = 0, u(r,O)= 0.
In this case the general solution of (6.1) given in Ref. [192] has the
following form
8
(6.4) v = (m3)7a(;) At),
where y ( t ) is determined by the differential equation
dY + Ay'ld
- = P(L),
(6.6)
dt
with initial condition y(0)= 0 and
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 303
Dynamic
t"
-\
'L
0 t
t-a I
I
I
I
t
I
I
I
I
1
I I I 5
9.I .
II 23 q4 25 r/a
FIG. 61. Dynamic curves VJTversus r/a for different time parameters (Ref. [lQ4]).
If the function @(F) is linear, i.e., 8 = 1, (6.6) has a closed-form solution
364 PIOTR PERZYNA
The time dependence of J&) for rectangular shape of the loading curve
and the linear function @(F) are plotted in Fig. 61,
In Ref. [192] the formulae for the components of stress a,,, a=, am,
strain E”,, eW, E80 and strain-rate tensors are given and the solutions for other
types of the function @(F) are discussed.
In order to compare the dynamic and the quasi-static solutions (cf.
Chapter IV, Section l), the same material and identical boundary conditions
in both cases should be assumed (see Ref. [194]). Consider the rigid/visco-
plastic material, the rectangular impulse and the linear function @(F) = F.
The corresponding solution of the quasi-static problem can be obtained
by putting into formulae (4.3)-(4.6) ,u -,00 and K --* 00 with the result
(6.10) v x / k=1 +
r] (g
and (6.9) for 8 = 1 and p(t) = const, yields
(6.12) V&.k = 1 + (- [ (-
1- - + 3 log 11 - exp (- ~ t ) ] ,
(6.13) zt(r,t)
1[7:(= -(); ]3 -1
yi(f+2vrlogz) b
FIG. 62. Comparison of the quasi-static with the dynamic solution (Ref. [1@4]).
FIG.63. Comparison of the quasi-static with the dynamic solution (Ref. [ 1941).
Experimental Summa
Cum C-SsOO(~ -1)7ax-'
C U M Of FhrkeS -
---
I I 1
c
9-8 a-g W" 1 8
strain mte ~,sec.?
FIG.64. Effect of rate of straining on yield stress ($5. R. Bodner and P. S. Symonds [20]).
rate on the yield stress was primarily responsible for the deviations between
theory and experiment. In theoretical treatment S. R. Rodner and P. S.
Symonds [20] have assumed that the influence of plastic strain rate E upon
yield stress 0 obeys the relation (2.89). The numerical values for y and 6
have been deduced from the experimental data of M. J. Manjoine [112] for
mild-steel specimens, and from experimental values collected by E. W. Parkes
[129] for various aluminium alloys (see Fig. 64). The predictions of the
rate-dependent, rigid-plastic theory are generally in satisfactory agreement
with test results for the final deformation (shape and magnitude), deformation
time, and strain time-history. On the other hand the application of an
overall strain-rate correction factor on the final deformation (cf. T. J. Mentel
[lla]) cannot be generally recommended since it may lead to serious errors.
The analysis in Ref. [20] depends on several special assumptions. In
particular, it has been assumed that the plastic region is finite but small
in length compared to the beam length and the deformations small enough
so that geometry changes could be ignored. Comparisons both within the
theory and with test results have shown that these assumptions are of
doubtful validity in the range of interest. They are not made in the analysis
given by T. C. T. Ting and P. S. Symonds [178] (see also T. C. T. Ting [189,
1901). This study shows that final plastic deformations are in good agreement
with those measured in tests performed by S. R. Bodner and P. S. Symonds
[201.
and H. Wolf [104]. These considerations have shown the very important
influence of the strain-rate effect on the deformation.
A rigid-(linearly viscoplastic) formulation of the theory for longitudinal
impact on a bar has been applied by V. V. Sokolovsky [lea] to several
problems of plane shear waves in a semi-infinite medium. Known solutions
of the heat equations were used.
T. C , T. Ting and P. S. Symonds [181] in the analysis of the longitudinal
impact on rigid/viscoplastic rods have used the linear stress-(strain rate)
law (cf. (2.93,))
(5.15) c = y*(a - ao).
Four cases have been solved : constant-velocity impact on a semi-infinite
bar, constant-velocity impact on a bar of finite length, impact of a finite
mass on a semi-infinite rod and impact of a finite mass on finite rod. In
every case the problem has been reduced to the solution of the heat equation.
G . I. Barenblatt and A. Y. Ishlinsky [6] have extended Taylor’s problem
to the viscoplastic material.
In Ref. [7] successive approximations have been suggested as a method
of analysis for the deformation of a strain-rate sensitive plastic cylinder
and the theoretical results have been compared with Whiffin’s experimental
data. Discussions of similar problems may be found in Refs, [loo, 101.
106, 150, 153, 190, 2001.
References
(Titles of Russian publications are translated).
1. ALDER, J. F., and PHILLIPS,V. A., The effect of strain rate and temperature on
the resistance of aluminium, copper and steel to compression, J . Inst. Metals 88.
80-86 (1954-55).
2. ALEXEIEV, N. A,, RAKHMATULIN, KH.A., and SAGAMONIAN, A. Y.,On the fundamental
equations for dynamics of soil, Zhuvlzal Prik. Mekh. Tekh. 8, 147 (1963).
3. ALTER,B. E. K., and CURTIS,C. W., Effects of strain rate on the propagation of
a plastic strain pulse along a lead bar, J . A p p l . Phys. 27, 1079-1085 (1956).
4. APPLEBY, E. J., and PRAGER, W., A problem in visco-plasticity, J . A p p l . Mech. 29,
381 - 384 (1962).
5. A m , A. K., and DIAZ,J. B., On a mixed boundary-value problem for linear hyper-
bolic partial differential equations in two independent variables, Avch. Rational
Mech. Anal. 10, 1-28 (1962).
6. BARENBLATT, G. I., and ISHLINSKY, A. Y., On the impact of a visco-plastic bar
on a rigid target (in Russian), Pvik. Mat. Mekh. 26, 497-602 (1962).
7. BEJDA,J., Analysis of deformation in a short visco-plastic cylinder striking a rigid
target, Avch. Mech. Stos. 16, 879-889 (1963).
8. BEJDA, J.. T h e wave problem of elastic/viscoplastic beams, J. MCan. (1966)
(forthcoming).
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 369
32. CAMPBELL, J. D.,and MARSH,K. J,, The effect of grain size on the delayed yielding
of mild steel, Philosofihical Magazine 7, 933- 952 (1962).
33. CHADWICK. P.,Cox, A. D., and HOPKINS,H. G., Mechanics of deep underground
explosions, Phil. Trans. Roy. SOC.London 2628A. 235-300 (1964).
34. CHIDDISTER, J . L., and MALVERN, L. E., Compression-impact testing of aluminium
at elevated temperature, Ex$. Mech. 8, 81-91 (1963).
35. CHU,S. C.,and DIAZ,J. B., Remarks on a mixed boundary-value problem for linear
hyperbolic partial differential equations in two independent variables, Arch. Ra-
tional Mech. Anal. 16, 187-196 (1964).
36. CINQUINI-CIBRARIO, M., Teoremi di unicitd per sistemi di equazioni a derivate
parziali in pih variabili indipendenti. Ann. di Mat., 48, 103-134 (1959).
37. CINQUINI-CIBRARIO, M.. Sistemi di equazioni a derivate parziali in pi& variabili
indipen’denti, Ann. di Mat. .44, 357-418 (1957).
38. CLARK, D.S.,and WOOD,D. S., The time delay for the initiation of plastic deforma-
tion at rapidly applied constant stress, Proc. A h . Soc. Testing Materials 40.
717-737 (1949).
39. CLARK, D. S.,and DUWEZ,P. E., The influence of strain rate on some tensile prop-
erties of steel, Proc. Amer. SOC.Testing Materials 60, 560-575 (1960).
40. CLARK,D. S.. The behavior of metals under dynamic loading, Trans Amer. SOC.
Metals 46, 34-52 (1954).
41. COLEMAN, B. D.,and NOLL,W., On the thermostatics of continuous media, Arch.
Rational Mech. Anal, 4, 97- 128 (1959).
42. COLEMAN, B. D., and NOLL.W., Foundations of linear viscoelasticity. Reviews oj
Modern Physics 81. 239-249 (1961).
43. COLEMAN, B. D.,Mechanical and thermodynamical admissibility of stress-strain
functions, Arch. Rational Mech. Anal. 9, 172-186 (1902).
44. COLEMAN,B. D.,and NOLL, W.,Material symmetry and thermostatic inequal-
ities in finite elastic deformations, Arch. Rational Mech. Anal. 16, 87-111
(1964).
45. CONLAN, J., The Cauchy problem and the mixed boundary-value problem for a
non-linear hyperbolic partial differential equation in two independent variables,
Arch. Rational Mech. Anal. 8, 355-380 (1959).
46. CONLAN, J., A generalized boundary value problem for uxy = f(r,y , u, uz, u,,),
Arch. Rational Mech. Anal. 18, 137-146 (1963).
47. COPSON,E. T., On the Riemann-Green function. Arch. Rational Mech. Anal. 1,
324-348 (1968).
48. COTTRELL, A. H., and BILBY,B. A., Dislocation theory of yielding and strain ageing
of iron, Pvoc. Phys. SOC.62A, 49 (1949).
49. COTTRLLL,A. H.,Deformation of solids at high rates of strain, an ‘eProc.of t h e
Conference on the properties of materials at high rates of strain, London 1957,”
pp. 1-12, Inst. Mech. Engrs., London (1957).
50. COURANT, R.,and FRIDRICHS, K. O., “Supersonic Flow and Shock Waves.” Inter-
science, New York. 1948.
61. COURANT, R., and LAX, P. D., Nonlinear partial differential equations with two
independent variables, Comm. Pure A$Pl. Math. 2, 255-273 (1949).
52. COURANT, R., ISAACSON, E., and REES,M., On the solution of non-linear hyperbolic
differential equations by finite differences, Comm. Pure Appl. Math. 6, 243- 255
(1962).
63. COURANT, R., Cauchy’s problem for hyperbolic quasi-linear systems of first order
partial differential equations in two independent variables, Comm. Pure Appl.
Math. 14, 257-2285 (1961).
64. COURANT, R., and HILBLRT,D., “Methods of Mathematical Physics,” vol. 2.
Interscience, New York, 1962.
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 37 1
55. COWPER, G. R., and SYMONDS, P. S., Strain hardening and strain rate effects in
the impact loading of cantilever beams, Technical Report No. 28, Brown
University, September 1957.
56. CRISTESCU, N., Some problems of the mechanics of extensible strings, stress wave
in anelastic solids, i n “I.U.T.A.M. Symposium, Brown University, April 1963”
(H. Kolsky and W. Prager. eds.). Springer-Verlag, Berlin. 1964, pp. 118- 132.
57. CRISTESCU. N . . Some dynamic problems of one-dimensional elastic/visco/plastic
bodies, Eleventh International Congress of Applied Mechanics, Munchen,
August 30 to September 5. 1964, (forthcoming).
68. DHUTLER,H., Experimentelle Untersuchungen hbex die Abhangigkeit der Zug-
spannungen von der Verformungsgeschwindigkeit, Phys. Z . 81,247- 259 (1932).
69. DIAZ,J . B., On an analogue of the Euler-Cauchy polygon method for the numerical
solution of uxy = f ( x , y , u, ux, u,,), Avch. Rational Mech. Anal. 1, 357-390(1958).
60. MAC DONALD, R. J., CARLSON, R. Id,, and LANKFORD, W. T.. The effects of strain
rate and temperature on the stress-strain relations of deep-drawing steel, Proc.
A m . Soc. Test. Materials 86, 704-120 (1956).
61. DoucrIs.A.,Someexistencetheoremsforhyperbolicsystemsof partialdifferential equa-
tions in two independent variables, Comm. Puve A p p l . Math. 5 , 119- 154 (1952).
62. DRUCKER, D. C., A definition of stable inelastic material, J . A p p l . Mech. 86, 101 - 106
(1969).
63. DRUCKER, D. C., On the role of experiment in the development of theory, in Proc.
Fourth U.S. National Congress Appl. Mech.. Berkeley, June 1962, Am. Soc. Mech.
Engrs., New York. 1963, pp. 18-33.
64. DRUCKER,D. C.. Stress-straiq-time relations and irreversible thermodynamics,
in “Proc. Inter, Symposium on Second-Order Effects in Elasticity, Plasticity
and Fluid Dynamics, Haifa, April 1962.” Pergamon Press, pp. 331-381.
65. DRUCKER. D. C., A continuum approach to the fracture of metals, i n “Fracture
of Solids”, (D. C. Drucker and J . J . Gilman, eds.), Wiley, New York 1983.
pp. 3-60.
66. DRUCKER, D. C., Some remarks on flow and fracture, Technical Report, Brown
University, May 1983.
67. DUWEZ,P. E., and CLARK,D. S., An experimental study of the propagation of
plastic deformation under conditions of longitudinal impact, Proc. A m . SOC.
Testing Malerials 47, 502-532 (1947).
68. ERINGEN.A. C.. “Nonlinear Theory of Continuous Media,” McGraw-Hill, New
York, 1982.
88. FREUDENTHAL, A. M.. and GEIRINGER, H., The mathematical theories of the inelastic
continuum, in “Encyclopedia of Physics.” vol. V1. Springer-Verlag, Berlin,
1958, pp. 229-443.
70. FRIEDRICHS, K. 0.. Nonlinear hyperbolic differential equations for functions of two
independent variables, American Journal of Mathematics, 70, 555- 589 (1948).
71. GOLDSMITH, W., ”Impact. The Theory and Physical Behaviour of Colliding Solids.”
Edward Arnold, London, 1960.
72. MAC GREGOR,C. W., and FISHER,J . C.. A velocity-modified temperature for the
plastic flow of metals, J . A p p l . Tech. 11. (1946).
73. GURTIN,M., and STERNBERG, E., On the linear theory of visco-elasticity, Arch.
Rational Mecfr. Anal. 11, 291-3356 (1962).
74. GIZATULINA, G. M., On the problem of deflection of a circular plate of visco-plastic
material, (in Russian) Zsledovaniya Upv. Plast. 3 (1964) Leningrad.
75. HARDING, J., WOOD, E. 0.. and CAMPBELL.J. D.. Tensile testing of materials
at impact rates of strain, J . Mech. Engrg. Sca. 2. 88-96 (1960).
76. HARTMAN, P., and WfrimER, A., On hyperbolic partial differential equations.
American Jouvnat o/ Mathematics, 74, 836-843 (1962).
372 PIOTR PERZYNA
77. HAUSER,F. E., SIMMONS, J . A., and DORN.J . E., Strain rate effects in plastic
wave propagation, Technical Report No. 3, University of California. June
1960.
78. HAWSER,F. E.. and WINTER. C. A., An experimental method for determining
stress-strain relations a t high strain-rates, Technical Report No. 4, University
of California, December 1960.
79. HILL, R., “The Mathematical Theory of Plasticity,” Oxford, 1950.
80. HOHENEMSER. K., and PRAGER,W., uber die Ansltze der Mechanik isotroper
Kontinua. Z A M M 19, 216-226 (1932).
81. HOPKINS.H. G., and WANG,A. J., Load-carrying capacities for circular plates of
perfectly plastic material with arbitrary yield condition, J. Mech. Phys. Solids
8, 117-129 (1954).
82. HOPKINS,H. G . . Dynamic expansion of spherical cavities in metals, Progvess i n
Solid Mechanics 1, 83- 164 (1960).
83. HOPKINS,H. G.. Dynamic anelastic deformations of metals, A#@. Mech. Rev. 14,
417-431 (1961).
84. HOPKINS,H. G., Non-linear stress-wave propagation in metals, i n “Progress in
Applied Mechanics - The Prager Anniversary Volume,” Macmillan, New York,
1963, pp. 55-71,
85. HOPKINS,H. G., Mechanical waves and strain-rate effects in metals, stress wave
in anelastic solids, Proc. I.U.T.A.M. Symposium, Brown University, April 1963,
(H.Kolsky and W. Prager, eds.). Springer-Verlag, Berlin, 1964. pp. 133-148.
86. JEFFREY,A.. The development of jump discontinuities in nonlinear hyperbolic
systems of equations in two independent variables, Arch. Rational Mech. Anal. 14,
27-37 (1963).
87. JOHNSON, J . E., WOOD,D. S., and CLARK,D. S., Delayed yielding in annealed low-
carbon steel under compression impact, Proc. A m . Soc. Tesling Materials 53,
765-764 (1953).
88. JOHNSON, J. E., WOOD, D. S.. and CLARK,D. S., Dynamic stress-strain relation for
annealed 2 s aluminium under compression impact, J . Appl. Mech. 90, 523- 529
(1953).
89. KALISKI,S., On certain equations of dynamics of an elastic-visco-plastic body.
The strain-hardening properties and the influence of strain rate. Bull. Acad.
Polon. Sci.; Sev. Sci. Techn. 11, 239 (1963).
90. KALISKI,S., NOWACKI, W. K., and WLODARCZYK, E., Propagation and reflection of
a spherical wave in an elastic-visco-plastic strain-hardening body, Pvoc. Vibr.
Probkms 6. 31-56 (1964).
91. KANTOROVICH. L. V., and AKILOV,G . P., “Functional analysis in normalized
spaces” (in Russian). Moscow 1969.
92. VON KARMAN, T., and DWWEZ. P.. The propagation of plastic deformation in solids.
in “6th International Congress for Applied Mechanics, Paris 1946’; also J . Appl.
Pbys. 21, 987-994 (1950).
93. KELLER,H. B., On the solution of semi-linear hyperbolic systems by unconditionally
stable difference methods, Comm. Pure Appl. Ma& 14. 447-471 (1961).
94. KELLER,H. B., and T H O M ~V., E , Unconditionally stable difference methods for
mixed problems for quasi-linear hyperbolic systems in two dimensions, Comm.
Pure A p p l . Math. 16, 63-73 (1002).
95. KOLSKY,H., and DOUCH, L. S., Experimental studies in plastic wave propagation,
J . Mcch. Phys. Solids 10, 195-223 (1962).
96. KRAFFT,J. M.,SULLIVAN, A. M.,and TIPPER,C. F.. The effect of static and dynamic
loading and temperature on the yield stress of iron and mild steel in compression,
Proc. Royal Soc. B H A , 114-127 (1964).
97. KRAFFT,J . M., and SULLIVAN, A. M., Effects of grain size and carbon content
FUNDAMENTAL PROBLEMS IN VISCOPLASTICITY 373
on the yield delay-time of mild steel. Trans. Amer. Soc. Mefak 61, 643-659
(1969) and discussion t o this paper by J. D. Campbell, pp. 669-6665.
98. KROYM.A.. Zur Ausbreitung von StoDwellen in Kreislochscheiben. Z A M M 38,
104 und 297 (1948).
99. KRZYZAfiSKI, M., “Second Order Partial Differential Equations” (in Polish).
Warsaw, vol. I (1957). vol. I1 (1962).
100. KUKUDJANOV, V. N., and NIKITIN,L. V., The propagation of waves in bars of
non-homogeneous elastic/visco/plastic material (in Russian), Isu. Akad. Nauk
S S S R , Mskh. Mash. 1, 53-59 (1960).
101. KUKUDJANOV, V. N., and NIKITIN,L. V., Impact of a bar with piecewise constant
yield limit on a rigid target (in Russian), mag. Zhur. 1, 177-183 (1961).
102. LANDAU, H. G., WHINER,J. H., and ZWICKY,E.E., Thermal stress in viscoelastic-
plastic plate with temperature-dependent yield stress. J . A p p l . Mech. 27.
297-302 (1960).
103. LAX, P. D., Nonlinear hyperbolic equations, Comm. Pure Appl. Math. 8, 231-258
(1963).
104. LEE, E. H., and WOLF,H., Plastic-wave propagation effects in high-speed testing,
J . A p p l . Mech. 18, 379-386 (1951).
105. LEE, E. H., Wave propagation in anelastic materials, Deformation and Flow of
Solids, i n “Colloquium Madrid 1955”, Springer-Verlag. Berlin, 1956, pp. 139- 136.
106. LUBLINER.J .. A generalized theory of strain-rate dependent plastic wave propaga-
tion in bars, J. Mech. Phys. Solids 12, 69-65 (1964).
107. LUDWIK,P., “Elemente der technologischen Mechanik.” Springer-Verlag, Berlin.
1909, pp. 44-47.
108. MAJCHER, G., Sur un problhme mixte pour 1’6quation du type hyperbolique,
Annales Polonici Mathematici, 6 , 12 1 - 139 (1958).
109. MAIDEN,C. J., and CAMPBELL, J. D., The static and dynamic strength of a carbon
steel at low temperature, Philosophical Magazine 8 , 872- 885 (1958).
110. MALVERN, L. E., The propagation of longitudinal waves of plastic deformation
in a bar of material exhibiting a strain-rate effect, J . Appl. Mech. 18, 203-208
(1951).
111. MALVBRN, L. E., Plastic wave propagation in a bar of material exhibiting a strain
rate effect, Qwr#.A w l . Math. 8, 406-411 (1951).
112. MANJOINE,M.. Influence of rate of strain and temperature on yield stresses of
mild steel, J. A M l . Mech. 11, 211-218 (1944).
113. MARSH,K. J., and CAMPBELL, J . D., The effect of strain rate on the delayed
yielding flow of mild steel, J . Mech. Phys. Solids 11, 48-69 (1963).
114. MENTEL,T. J,, The plastic deformation due to impact of cantilever beam with
an attached tip mass, J. Appl. Mech. B6, p. 515 (1958).
115. MIKLOWITZ, J., The initiation and propagation of the plastic zone in a tension
bar of mild steel as influenced by speed and stretching and rigidity of testing
machine, J. A p p l . Mech. 14, A-31-A-38 (1947).
116. NADAI,A. L., “Theory of Flow and Fracture of Solids.” McGraw-Hill, New York.
vol. 1, 1950, vol. 2, 1963.
117. NAGHDI,P. M., Stress-strain relations in plasticity and thermoplasticity, ar
“Plasticity” (E. H. Lee and P. S. Symonds, eds.), Pergamon Press, London,
1960, pp. 121-167.
118. NAGHDI,P. M., and MURCH,S. A., On the mechanical behavior of viscoelasticl
plastic solids, J. A$@. Mech. 80, 321-328 (1963).
119. NIKITIN,L. V., The propagation of elastic-visco-plastic waves in a thick-walled
tube (in Russian), Zzu. Vurov. Mash. 8, 14-23 (1958).
120. NIKITIN, L. V., The propagation of transversal elastic-visco-plastic waves in
beams and plates (in Russian), Zag. Sbw. 80, 31-46 (1960).
374 PIOTR PERZYNA
at the University of Texas, Parts I and 11, in “Plasticity” (E. M. Lee and
P. S. Symonds, eds.). Pergamon Press, London, 1960, pp. 453-487.
142. PRAGER. W., Mkanique des solides isotropes au deli du domaine blastique, Mlmo-
rial Sci., Math. 8 7 , Paris, (1937).
143. PRAGER,W.. Discontinuous fields of plastic stress and flow, Proc. of the Second
U.S. National Congress of Applied Mechanics. pp. 21-32 (1954).
144. PRAGER, W., Non-isothermal plastic deformation, Proc. Kon. Ned. Akad. Wet. B81.
176- 182 (1958).
145. PRAGER,W., Linearization in visco-plasticity, &err. Ing. Avchiu, 16, 152- 167
(1961).
146. PRAGER,W., “Introduction to Mechanics of Continua.” Ginn and Company.
Boston, 1961.
147. PRANDTL. L., Ein Cedankenmodell zur kinetischen Theorie der festen Korper,
Z A M M N, 85-106 (1928).
148. ~ O U S EG., , Sulla risoluzione del problema misto per le equazioni iperboliche non
lineari mediante le differenze finite, A n n . di Mat. 46, 313-341 (1958).
149. PUGH.H. L. D., CUANG,S. S.. and HOPKINS,B. E., Tensile properties of a high-
purity iron from - 1H6 “C to 200 “C a t two rates of strain, Phil. Mag. %, 753-768
(1963).
150. K A J N A K ,S., HAUSER, F. E..and DORN,J . E., Theoretical prediction of strain
distribution under impact loading, Technical Report No. 1, University of
California, June 1961.
151. RAKHMATULIN, H. A , , and SHAPIRO, G . S.. Propagation of a disturbance in a non-
linear elastic and non-elastic medium (in Russian), I z v . Akad. Nauk S S S H . .
Otdelenie tech. nauk 8 , 68-88 (1955).
158. REINER,M., Plastic yielding in anelasticity. ./. Mech. Phys. Solids 8. 265-2261
(1960).
153. RUBIN,R. J., Propagation of longitudinal deformation waves in a pre-stressed
rod of material exhibiting a strain-rate effect, ./. A p p l . P h y s . 86.. 528-5636
(1954).
1.54. SCHAUDBR, J ., Cauchysches Problem fiir partiell: Differentialgleichungen erster
Ordnung. Anwendung einiger sich auf die Absolutbetrlge der Lesungen be-
ziehenden Abschatzungen, Comentaviz Math. Helu. H. 263-283 (1937).
.,
155. S C H A U ~ E JR , Das Anfangswertproblem einer quasi-linearen hyperbolischen
Differentialgleichung zweiter Ordnung in beliebiger Anzahl von unabhlngigen
Veranderlichen, Fundamenta Mathematicae Y4, 213- 246 (1935).
156. SIMMONS, J. A,, HAUSER, F. E.,and DORN,J . E., Mathematical theories of plastic
deformation under impulsive loading. University of California, Publications
an Engineering 5, 177-230, (1962).
157. SKBMFTON, A. W., and BISHOP,A. W., “Building materials: their elasticity and
inelasticity” (M. Reiner, ed.). North-Holland Publ., Amsterdam, 1964,
Chapter X, Soils.
158. SKRIPXIN, V. A , , On the model of dynamical work-hardening continuous media
in the Mises-Reuss theory of plasticity (in Russian). Prak. Mat. Mekh. 87,
109-115 (1983).
159. SOKOLOV, D. D., On the problem of determination of the resistance of metals to
the plastic deformation in dependence on the strain rate and temperature,
(in Russian) Doklady Akademii Nauk S S S R 6I (1949).
160. SOKOLOV, D. D., “Strength of Metals in Plastic Deformations.” (in Russian),
Moscow. 1963.
161. SOKOLOVSKY, V. V., Propagation of elastic-viscoplastic waves in bars (in Russian),
Prikl. Mat. Mekh. 18, 261-280 (1948).
376 PIOTR PERZYNA
184. TRUESDELL, c., General and exact theory of waves in finite elastic strain, Arch.
Rational Mech. Anal. 8, 263-296 (1961).
185. TRUESDELL, C., and TOUPIN,R.. Static grounds for inequalities in finite strain
of elastic materials, Arch. Rational Mech. A n d . 18, 1-33 (1963).
186. VISHMAN, F. F., ZLATIN,N. A., and MODDE. V. S., Resistance to the deformation
of metals a t the rate of lo4 to 10-’ mm./sec., Zhurn. Tekhn. Fiz. 19 (1949).
187. VOLOSHENKO-KLITOVITSKY, Y. Y . , On the conformity of the law of yield limit
at low-temperatures loads, Izv. Akad. Nauk S S S R , Mekh. Mash. 1, 154-156
(1962).
188. WEHRLI,C.. und ZIEGLER,H., Einige mit dem Prinzip von der groBten Dissi-
-
pationsleistung vertrlgliche Stoffgleichungen, Z A M P 18, 372 393 (1962).
189. WHIFFIN,A. C.. The use of flat-ended projectiles for determining dynamic yield
stress, 11. Tests on various metallic materials, Proc. Royal Soc. 1R4A, 300-332
(1948).
190. WHITE,M. P., On the impact behavior of a material with a yield point, J . A p p l .
Mech. 16. 39-52 (1949).
191. WIERZBICKI, T., A thick-walled elasto-visco-plastic spherical container under
stress and displacement boundary conditions, Arch. Mech. Stos. 15. 297-308
(1963).
192. WIERZBICKI. T.. Impulsive loading of a spherical container with rigid-plastic and
strain-rate sensitive material, Avch. Mech. Slos. 15, 775- 790 (1963).
193. WIERZBICKI. T., Bending of rigid/visco-plastic circular plates, Arch. Mech. - 9 0 s .
16 (1964); Bull. Acad. Pol. Sci.. Se’r. sci. tech. 18, 611-618 (1964).
194. WIERZBICKI, T., On the impulsive loading of a spherical vessel, Bull. Acad. Pol.
Sci., Sbr. sci tech. 12, 217-224 (1964).
196. WIERZBICKI, T.. Dvnamics of visco/plastic circular plate, ‘4rch. Mech. Stos. 17,
851-869 (1965).
196. WINTNBR, A., On the conditions of validity of Riemann’s method of integration,
Quart. ApPl. Maths. 16, 94-98 (1957).
197. ZIEGLER,H., Uber ein Prinzip der groBten spezifischen Entropieproduktion und
seine Bedeutung fur die Rheologie. Rheologica Acta 2, 230-235 (1962).
198. ZIEGLER,H., Some extremum principles in irreversible thermodynamics with
application to continuum mechanics, Progress in Solid Mechanics 4. 91 - 193
(1963).
199. ZIEGLER,H.. Thermodynamic considerations in continuum mechanics. The 1964
Minta Martin Lecture, Massachusetts Institute of Technology.
200. ZVEREV,I. N., The propagation of perturbation in visco-elastic and visco-plastic
bars (in Russian), Prikl. Mat. Mekh. 14, 295-302 (1960).
201. DRUCKER, D. C.. and PRAGER,W., Soil mechanics and plastic analysis of limit
design, Quart. A p p L Math. 10, 157-165 (1952).
202. LEE, E. H., and TUPPER,S. J.. Plastic deformation in a steel cylinder striking
a rigid target, J. A p p l . Mech., el. 63-72 (1954).
Author Index
Numbers in parenthesesare reference numbers and are included to assist in Loeating referencesin which auto&
name8 are not mentioned in the test. Numbers in italics refer to pages on which the complete referenceare listed.
379
380 AUTHOR INDEX
Kinney. W. D., 166, 205(12). 207(12). Maiden, C. J.. 251. 252, 253. 290(109).
210(12, 13), 214(12), 216, 231, 232, 296(109), 369, 373
24 1 Malvern, L. E., 252, 254, 255, 273,
Klotter, K., 217(38), 242 277, 283, 290(34), 295(34), 296(39),
Kolsky, H., 249, 250, 2.57, 258, 372 370, 373
Krafft, J . M., 282, 258(97). 290, 296(96), Manjoine. M., 246, 247. 252. 296(112).
372 367, 373
Kromm, A., 340, 373 Marsh, K. J.. 246, 247. 249, 251(32),
Krzyzabski, M., 308, (99). 373 370, 373
Kuchemann, D., 49, 87 Menkes, J., 47, 87
Kukudjanov, V. N.. 368(100), 373 Mentel, T. J.. 365, 367, 373
Kuo, H. L.. 51(66), 52, 88 Michael, 1). H.. 52(60. 71). 53(71). 88
Kuo, J . K.. 163(8), 199(8, 18), 200. Michalke, A., 42, 87
201(8, 18), 202. 230(8, 43, 44). 233. Miishkes, A. D., 191(31), 193(31), 242
234, 241 Mikishev, G. N., 117(8), llE(8). 153
Miklowitz, J., 246, 373
I, Miles, J . W., 29, 48(68), 48, 62. 63
(77, 78), 71, 78, 86, 87, 88
Landau, H. G., 373
Minorski. N., 1.56. 179(1), 240
Landau, L., 48. 87, 191(33), 193(33).
Modde, V. S., 205(186), 377
242
Moiseev. N. N.. 91(1), 92(1), 93(l),
Lankford. W. T.. 252, 371
94(1), 95(1), 127(11), 753
Lax, P. I)., 303, 370, 373
Murch, S. A., 269, 261. 262, 267, 270,
Ledoux, P.. 158(.5), 240
373
Lee, E. ti., 245, 367, 368. 373, 377
Murphy, J. W.. 86
Lessen, M., 40, 42, 06
Levinson, N.. 220(26), 247
Liepmann, H. W.. 49, 57(74), 87. 88 N
Lifshitz, E. M., IYI(33). 193(33). 242
Nadai, A . L., 254, 373
Lighthill, M. J.. 15, 86
Naghdi. P. M.. 259. 261, 262. 267. 270.
Lin. C. C . , 6, 9(7, 9), 12, 14, 16, 20(7),
273, 281(117), 373
29. 36(7), 49, 58, 70(5), 71(5), 85. 86.
Nikitin, L. V., 318(110), 361, 368(100).
87
373, 374
Lipps, F. B.. 52, 88
Noll, W., 260, 261, 272, 370, 374
Lock, 1s. C., 53(73), 88
Nowacki, W. K., 316(90), 372
Loude, W.. 221. 222(41). 242
Lubliner. J., 250, 369(106), 373
Ludwik, P.. 264, 373 0
Olszak, W.. 264, 282(127, 128). 290(126),
M 300(128), 316(127), 374
MacDonald, K. J., 262, 377 Orr. W. M. F.. 28. 86
McDuff, J . N., 191(32). 193(32), 242
MacGregor. C. W., 371
MacLachlan, N. W., 223(42). 242
P
Maharem, N., 184(25), 241 Panofsky, H. A., 40, 87
Majcher, G., 313, 315, 337, 373 Parkes, E. W.. 365. 367. 374
Marwin, J., 158(6), 171, 176(6), 188, Patterson, A. M., 68, 88
189, 190, 240 Pearson, K., 193, 242
382 AUTHOR INDEX
A D
A-process def., 298 Delay-time, 246, 261
Admissible systems (NV), 158f.. 162 Difference equation, 321
Alfv6n waves, 63 Dislocation theory of crystalline material,
Autonomous system, admissible (NV). 258
160 Dissipation energy, 264
trajectories of, 170 Drucker’s postulate, 265
Amp-function, 196f. Duffing approximation, 228
Ateb-function, 191, I96 Dynamic test, 245
properties, influenced by temperature
and irradiation, 246, 262
B
B-process def., 299 E
Banach’s theorem, 305
Beta function, 193 E-curve def. (NV). 178
incomplete, 193 Elastic energy, 264
Bounding surface def. (NV), 169 Elastic-viscoplastic def., 269
Elastic/viscoplastic def., 259
Elastic/visco-(perfectly plastic) material,
C 277, 281
Elastic/viscoplastic region, boundary of,
Cam-function, 197, ff. 322
Cantilever-bepm specimens, 306f. e-tube (NV),217
Characteristics, 317, 321 Escalator method, 140, 147f.
Circular cylindrical container with apher- Exchange of stabilities, 19, 40, 65, 73
ical caps (FOL). 138
Coaxial cylinders, horizontal (FOL), 123
Configuration space (NV), 164 F
Conical container (FOL), 111 F-curve def. (NV), 184
Constitutive equations, 200ff. Finite differences, method of, 310
for elastic-viscoplastic material, 270 Flow surface, 202
for rate-sensitive plastic material, 272 convexity of, 266
Constitutive inequalities, 272 Fourier method (FOL), 97
Convexity of flow surface, 209 Fracture of time- and temperature-
Coordinate functions, choice of, 96 dependent materials, 300
Creep process, 264 Frank-Read dislocation source, 259
Critical layer (HS), 16f., 24, 03. 66 Free oscillations def. (FOL), 94
Critical state, 264 Frequency equation, generalized, (NV),
Cubical dilatation, plastic rate of, 282 21 1
Cylindrical container with horizontal axis approximate method of calculation,
(FOL). 132 102
384
SUBJECT INDEX 386
M
Magnetohydrodynamic waves, 52f. (1
Mathieu equation, 223 Quasi-static solutions. 36Off.
386 SUBJECT INDEX
R in r o t d n g system, 49
Rapid load path, 267 with variable Coriolis parameter, 61.
Rate-sensitive materials, 273 57
Rayleigh’s necessary condition for in- Stability of normal-mode vibration (NV),
stability, 2, lOf., 15, 26f., 36, 61, 66, 82 230
Rayleigh stability equation, 7, 9f., 12, 14, Stable elastic/viscoplastic material, dy-
17, 20, 31, 42. 44, 47, 50 namical condition for, 282
Rayleigh-Taylor instability, 46 Stable inelastic material, 264f.
Relaxation equation, 300 State depending on loading path, 244
Relaxation process for general states of on load history and on time, 246
stress, 298 Static yield function, 273
Reynolds stress, 14f., 42, 64, 83 Static properties, changed by previous
Riemann function def., 308 dynamic loading, 246, 251
right circular-cylindricalcontainer (FOL), Steady-state forced vibrations, 162, 204
101 Strain acceleration, influence of, 268
with annular base (FOL), 101 Strain-hardening effect, reduction of, 249
Ritz method, 93, 96f. Strain-rate sensitivity, 246, 266
Rossby waves, Slf., 67 of beams under impact, 366f.
of metals, 276
of plastic materials, 290, 291
S of yield limit. 246
Sam -f unction, 197ff. Strain rate and temperature, simultane-
Schauder’a theorem, 306 ous influence of, 291
Semicircle theorem, 17, 49, 63, 68, 79, 83 Stress space, 262
Semi-linear hyperbolic partial differential Stress-strain curve, dynamic, 248
equation, 303 Stress-strain relation. static 277
Shear waves, cylindrical, 316, 326ff. dynamic, 277
a,&-curves, dynamic, 248 Stress waves, numerical examples. 340
Simple trajectories, existence of, (NL), in elastic/visco-plastic soil, 316
233ff. propagation, 302
Soils, experimental results for, 263 solution in the elastic region, 336ff.
dilatation rate. of, 281 Strong discontinuity, 318
Sokolovsky-Malvern theory. 269 Strongly nonlinear (NV).161
Sound waves, 43. 47, 49 Strubel. method of, 231
Spherical container (FOL),118 Strutt chart, 224, 227
Spherical container, impulsive loading of
(PV),362 T
Spherical problem, quasistatic sol., 350 T-curves def. (NV).186f.
Squire’s theorem, 2, 6, 43, 46, 48. Mm.,82 Temperature-dependent materials, 29Of.
Stability boundary, 64, 63f., 69f.. 73. 76ff. Thermostatic and thermodynamic in-
Stability in the large (PV),266 equalities, 272
in the small, 265 Thick-walled spherical container, quasi-
Stability of compressible fluid, 47, 67 static sol., 360
of electrically-conducting fluid in a Time-dependence of stress and strain, 244
magnetic field, 62 Toroidally shaped container (FOL),143
of fluid of variable density, 44, 60, Transversal motion, influence of, 250
64, 68. 60 Transversals (= P-curves) def. (NV).
of non-parallel flow, 22 184f.
SUBJECT INDEX 387