Classification and Examples of Differential Equations and Their Applications, 2019 PDF
Classification and Examples of Differential Equations and Their Applications, 2019 PDF
of Differential Equations
and their Applications
Mathematics and Physics for Science and Technology
Series Editor: L.M.B.C. Campos
Director of the Center for Aeronautical
and Space Science and Technology
Lisbon University
Volume IV
Ordinary Differential
Equations with Applications
to Trajectories and Oscillations
Book 9
By
L.M.B.C. Campos
Director of the Center for Aeronautical
and Space Science and Technology
Lisbon University
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 334 87-2742
This book contains information obtained from authentic and highly regarded sources.
Reasonable efforts have been made to publish reliable data and information, but the
author and publisher cannot assume responsibility for the validity of all materials or the
consequences of their use. The authors and publishers have attempted to trace the copyright
holders of all material reproduced in this publication and apologize to copyright holders if
permission to publish in this form has not been obtained. If any copyright material has not
been acknowledged, please write and let us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted,
reproduced, transmitted, or utilized in any form by any electronic, mechanical, or other means,
now known or hereafter invented, including photocopying, microfilming, and recording, or in
any information storage or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www
.copyright.com (http://www.copyright.com/) or contact the Copyright Clearance Center,
Inc. (CCC), 222 Rosewood Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit
organization that provides licenses and registration for a variety of users. For organizations
that have been granted a photocopy license by the CCC, a separate system of payment has been
arranged.
Contents����������������������������������������������������������������������������������������������������������������� vii
List of Classifications and Tables�������������������������������������������������������������������������xi
Preface������������������������������������������������������������������������������������������������������������������� xxi
Acknowledgments��������������������������������������������������������������������������������������������� xxv
About the Author��������������������������������������������������������������������������������������������� xxvii
Physical Quantities������������������������������������������������������������������������������������������� xxix
vii
viii Contents
Classification 10.1�������������������������������������������������������������������������������������������������94
Classification 10.2����������������������������������������������������������������������������������������������159
References�����������������������������������������������������������������������������������������������������������203
Index���������������������������������������������������������������������������������������������������������������������209
List of Classifications and Tables
Classifications
C 10.1 500 Standards of Single and Simultaneous Systems
of Ordinary Differential Equations����������������������������������������������������� 94
A. Linear Differential Equation with Constant Coefficients�������������������� 95
A. a. Unforced: Standards I–VI��������������������������������������������������������� 95
A. b. Forced—Direct Method: Standards VII–XII��������������������������� 95
A. c. Forced—Method of Inverse Polynomial of Ordinary
Derivatives: Standards XIII–XX������������������������������������������� 96
A. d. Complete Integrals�������������������������������������������������������������������� 97
B. Linear Differential Equation with Homogeneous Coefficients���������� 97
B. a. Unforced: Standards XXI–XXIV����������������������������������������������� 97
B. b. Forced—Direct Method: Standards XXV–XXX���������������������� 97
B. c. Forced—Method of Inverse Polynomial of Homogeneous
Derivatives: Standards XXXI–XXXVII�������������������������������� 98
B. d. Complete Integrals��������������������������������������������������������������������� 99
C. Linear Finite Difference Equation with Constant Coefficients����������� 99
C. a. Unforced: Standards XXXVII–XL��������������������������������������������� 99
C. b. Forced: Standards XLI–XLIV�������������������������������������������������� 100
D. Method of Variation of Parameters: Standard XLV��������������������������� 100
E. Influence or Green’s Function: Standard XLVI����������������������������������� 101
F. Fourier Series-—Discrete Spectrum: Standard XLVII������������������������� 101
G. Fourier Integral—Continuous Spectrum: Standard XLVIII�������������� 101
H. Laplace Transform—Initial Conditions: Standard XLIX������������������� 102
I. First-Order Differential Equations��������������������������������������������������������� 102
I. a. Solvable by Quadratures: Standards L–LIV�������������������������� 102
I. b. Ricatti: Standards LV–LXIII������������������������������������������������������ 103
I. c. Homogeneous Differential Equation: Standards
LXVI–LXV����������������������������������������������������������������������������� 104
J. Differentials����������������������������������������������������������������������������������������������� 104
J. a. First-Order Differential in Two Variables: Standards
LXVI–LXXIII������������������������������������������������������������������������� 104
J. b. First-Order Differential in Three Variables: Standards
LXXIV–LXXVI���������������������������������������������������������������������� 106
J. c. Three-Dimensional Potentials: Standards
LXXVII–LXXXI��������������������������������������������������������������������� 106
xi
xii List of Classifications and Tables
Tables
T 10.1 Combinations of Boundary Conditions for a Circular Plate
with a Circular Hole������������������������������������������������������������������������������ 49
T 10.2 Constants in the Shape of Neutral Surface of Heavy
Circular Plate with a Circular Hole����������������������������������������������������� 51
T 10.3 Boundary Conditions for the Bending of a Plate������������������������������ 61
T 10.4 General Solution (10.350a, b) of Coupled System of
Finite Difference Equations (10.346a, b)���������������������������������������������� 81
Preface
xxi
xxii Preface
Summary of Volume IV
Following the three volumes on complex, transcendental, and generalized
functions, this fourth volume is the first volume on boundary and initial
value problems, and concerns ordinary differential equations with a vari-
ety of applications, including trajectories and oscillations. The differential
equations are classified in 500 standards indicating their properties, meth-
ods of solution, and examples, and include: linear and non-linear ordinary
differential equations and simultaneous systems; in particular, linear equa-
tions with constant coefficients having a characteristic polynomial, and
linear equations with variable coefficients leading to 23 special functions.
The applications consist of 500 problems with physical formulation, math-
ematical solutions, and interpretation of results; the problems include linear
and non-linear oscillators with one or several degrees-of-freedom, with or
without damping, and forcing including ordinary, parametric, multiple, and
non-linear resonance. Other applications include trajectories, orbits, electro-
mechanical systems, and chains with deterministic and chaotic behavior;
also included are deformations and buckling of elastic beams, membranes,
and plates. Included in the contents are:
The fourth volume of the series justifies renewing some of the acknowl-
edgments also made in the first three volumes, to those who contributed
more directly to the final form of the volume: Ms. Ana Moura, L. Sousa, and
S. Pernadas for help with the manuscripts; Mr. J. Coelho for all the drawings;
and at last, but not least, to my wife as my companion in preparing this work.
xxv
About the Author
xxvii
xxviii About the Author
xxix
xxx Physical Quantities
t — time: 2.1.1
— time of reception: 2.1.11
t0 — time of emission: 2.1.11
u — longitudinal displacement of a bar: 6.4.1
u — displacement vector: 6.4.1, 6.8.5
v — velocity: 2.1.1
v — velocity vector: 5.2.1
w — group velocity: N6.16
x — independent variable in ordinary differential equation: 1.1.1
— Cartesian coordinate: 2.1.1
x — position vector of observer: 2.1.11
y — dependent variable in ordinary differential equation: 1.1.1
— Cartesian coordinate: 2.1
y — position vector of source: 2.1.11
z — Cartesian coordinate: 6.1.1
— specific impedance: N7.47
N — unit outer normal to a surface: 5.5.1
— turning moment: 6.7.6
N n — normal turning moment: 6.8.16
N n — augmented normal turning moment: 6.9.5
N x , N y — turning moments: 6.7.6
N rs — undamped modal matrix: 8.2.14
P — reduced pressure perturbation spectrum: N7.23
P2 N — dispersion polynomial: 8.2.2
Pij — dispersion matrix: 8.2.2
Q — heat: 3.9.11
— heat source: 6.4.2
— concentrated torque: E10.10.1
Qr — transformation matrix: 8.2.7
R — electrical resistance: 2.1.6, 4.7.1, 8.9.2
— perfect gas constant: 5.5.17
— radius of curvature: 6.1.1
— reflection coefficient: N7.47
S — entropy: 3.9.11
— dynamo parameter: 4.7.3
— strain: 6.4.1
— surface adsorption coefficient: N7.48
Sa — row of components of the strain tensor: 6.8.7
Sij — strain tensor: 6.5.1, 6.8.8
— scattering matrix: N7.53
T — temperature: 3.9.11
— thrust: 4.8.13
— tangential tension along an elastic string: 6.1.1
— transmission coefficient: N7.51
T — effective axial tension: 6.3.12
T — stress vector: 6.6.2, 6.8.3
Ta — row of components of the stress tensor: 6.8.7
Tij — stress tensor: 6.5.1, 6.8.3
o
Tij — residual stresses: 6.8.6
U — internal energy: 3.9.11
Physical Quantities xxxv
y ′′′ + y ′′ − 2 y = x 2 + x + 1 + e 3 x + e x + sinh ( 2 x )
(10.1)
+ e −2 x cosh x + e 2 x sin x + e − x cos x + e x sin x
and obtain the complete integral. The linear third-order ordinary differential
equation (10.1) has characteristic polynomial (10.2a–c) of ordinary derivatives
(10.2d):
( )
P3 ( D ) = D3 + D2 − 2 = ( D − 1) D2 + 2 D + 2 = ( D − 1) ( D + 1 − i ) ( D + 1 + i ) ,
(10.2a–c)
d
D≡ : P3′ ( D ) = 3D2 + 2 D = D ( 3D + 2 ) ; (10.2d, e)
dx
the roots (10.3a) of (10.2c) specify the general integral (10.3b) of the unforced
equation:
1
2 Differential Equations and Applications
1 1 1
(
P3 ( D )
x2 + x + 1 = − )
2 1 − D2 ( 1 + D )/2
(
x2 + x + 1 )
(10.4)
1 D2 x2 x
= − 1 +
2 2
( ) (
+ O D3 x 2 + x + 1 = − )− − 1.
2 2
1 x2 x
y ( x) − A x −
x
( B cos x + C sin x ) + + +1
2 2
e3x x ex 1 e2x e −2 x 1 e − x e −3 x
= + + − + +
P3 ( 3 ) P3′ ( 1) 2 P3 ( 2 ) P3 ( −2 ) 2 P3 ( −1) P3 ( −3 )
because; (i) the exponential exp ( 3x ) is non-resonant (1.129a, b); (ii) the
exponential exp(x) is resonant (1.134a, b); (iii) the hyperbolic sine sinh(2x)
is non-resonant (1.147a, b); (iv) the product of exponential exp ( −2x ) and
hyperbolic cosine cosh(x) is not resonant (1.165a, b); (v) the product of the
exponential exp(2x) and circular sine sin(x) is also non-resonant (1.173a, b);
(vi) the product of the exponential exp ( −2x ) by the circular cosine cos x is
resonant (1.174a–c); (vii) the product of the exponential exp(x) by the circu-
lar cosine cos(x) is non-resonant (1.176a, b). The characteristic polynomial
(10.2b) and its first-order derivative (10.2e) evaluate the coefficients (10.6a–i)
in (10.5):
1 1 1 1 1 1 1 1
= , = , = , = − , (10.6a–d)
P3 ( 3 ) 34 P3′ ( 1) 5 P3 ( 2 ) 10 P3 ( −2 ) 6
1 1 1 1 1 1 1 − 5i
=− , =− , = = , (10.6e–g)
P3 ( −1) 2 P3 ( −3 ) 20 P3 ( 2 + i ) 3 ( 1 + 5 i ) 78
1 1 1− 2i 1 1 1+ i
=− =− , = =− , (10.6h, i)
P3′ ( i − 1) 2 (1 + 2i ) 10 P3 ( 1 + i ) 4 ( i − 1) 8
Examples 10.1 to 10.20 3
x2 + x + 2 e3x
y ( x ) = A e x + e − x (B cos x + C sin x) − +
2 34
x e x e 2 x e −2 x e − x e −3 x e 2 x
+
5
+
20
+
12
−
4
−
40
+
78
( sin x − 5 cos x ) (10.7)
x e− x ex
−
10
( cos x + 2 sin x ) − ( sin x + cos x ) ,
8
d
δ≡x : P3 ( δ ) = δ ( δ − 1)( δ − 2 ) + 4 δ ( δ − 1) + 2δ − 2
dx (10.9a–d)
(
= ( δ − 1) δ ( δ + 2 ) + 2 = ( δ − 1) δ + 2 δ + 2 ,
2
)
coincides with (10.9d) ≡ (10.2c). Thus, the substitutions (10.10a, b), which imply
(10.10c, d):
y n+ 3 + y n+ 2 − 2 y n = 3n + 2 n− 1 + ( −2 )
n− 1
+
( − )n + ( −3 )n , (10.12)
2 2
and obtain the complete solution. The characteristic polynomial of forward
finite differences (10.13a) of (10.12) is (10.13b) ≡ (10.9d) ≡ (10.2c):
∆y n = y n+ 1 : P3 ( ∆ ) = ∆ 3 + ∆ 2 − 2; y ( x ) ↔ y n , y ( k ) ( x ) ↔ y n+ k , e ax ↔ a n ,
(10.13a–e)
the changes (10.13c–e) imply:
{ }
n
e( ↔ ( a + i b) =
a + ib ) x 1/2
exp i arc tan ( b/a )
n
a2 + b2
(10.13f, g)
2 n/2
= a + b 2
exp i n arc tan ( b/a ) ,
for example:
e − x ± i x ↔ ( −1 ± i ) = ( )
n n
2 e ± i 3π/4 = 2 n/2 e ± i 3π n/4
(10.13h, i)
= 2 n/2 cos ( 3 π n/4 ) ± i sin ( 3 π n/4 ) ,
and also:
( i a ) ± ( − i a ) = an e iπ/2 ± e − iπ/2 ( ) ( )
n n n n
e iax ± e − iax
cos, sin ( ax ) = ↔
{2, 2i} {2, 2i} {2, 2i}
(10.13j–m)
e inπ/2 ± e − inπ/2
= an = a n cos, sin ( nπ/2 ) .
{2, 2i}
Use of the relations (10.13a–c, h, i) transforms: (i) the differential equation
(10.1) with the fourth, sixth, and seventh terms on the right-hand side (r.h.s.)
to the finite difference equation (10.12); (ii) the corresponding terms of the
complete integral (10.7) to solution (10.14):
2 n− 2 ( −2 ) ( − ) − ( −3 ) . n− 2 n n
πn π n 3
n
y n = A + 2 n/2 B cos + C sin + + − −
4 4 34 5 3 4 40
(10.14)
Examples 10.1 to 10.20 5
The solution (10.14) [complete integral (10.11)] of the linear finite difference
(differential) equation (10.12) [(10.8)]: (i) have been obtained using the trans-
formations (10.13a–m) [(10.10a–d)] applied to the equivalent linear differ-
ential equation with constant coefficients (10.1) and its complete integral
(10.7); (ii) can be confirmed using directly the methods in the section(s) 1.9
(1.7–1.8).
k k a 2 −2 λt k a 2 −2 λt
α = 0: Φ ( t ) = x ( t ) = cos 2 ( ωt ) = 1 + cos ( 2ωt ) ,
2
e e
2 2 4
(10.15a, b)
1 k a 2 −2 λt
Ev ( t ) = m x ( t ) = ω sin ( ωt ) + λ cos ( ωt )
2 2
e
2 2 ω 20
k a 2 −2 λt 2
= e ω sin 2 ( ωt ) + λ 2 cos 2 ( ωt ) + 2 λω cos ( ωt ) sin ( ωt ) (10.15c)
2 ω 20
k a 2 −2 λt 2 λω λ2 − ω2
= e 1 + 2 sin ( 2ωt ) + cos ( 2ωt ) ,
4 ω 0 ω 2
0
k a 2 −2 λt λω λ2
E ( t ) = Φ ( t ) + Ev ( t ) = e 1 + 2 sin ( 2ωt ) + 2 cos ( 2ωt ) . (10.15d)
2 ω 0 ω 0
6 Differential Equations and Applications
The exact average values over a period starting at an arbitrary time t and
ending at t + τ1 are given by:
Φ ( t ) , Ev ( t ) , E ( t )
k a2 2 λω λ2 − ω2 2 λω 2λ 2
= I0 + I+ , I0 + I− + I+ , 2 I0 + 2 I− + 2 I+ ,
4 ω02
ω02
ω0 ω0
(10.16a–c)
t + 2 π/ω
t+ τ
ω e −2 λt ω −2 λt
I0 =
1
τ ∫ t
e −2 λt
dt =
2 π −2 λ t
=
4π λ
e (
1 − e −4 πλ/ω , (10.17a) )
1 t+ τ
ω t + 2 π/ω
I± =
τ ∫t
e −2 λt cos, sin ( 2ω t ) dt =
2π
Re, Im
∫
t
e −2 λt + i 2 ω t dt
=−
ω
4π
Re, Im
λ − iω
e {
1 −2( λ− iω )t t + 2 π/ω
t } (10.17b)
ω λ + iω −2( λ− iω )t
=
4π
Re, Im 2
λ + ω
2
e
1 − e −4 πλ/ω ;
( )
the latter splits into:
ω 1 − e −4 πλ/ω −2 λt
{I + , I − } = e
4π ω 20 (10.18a, b)
{
× λ cos ( 2ω t ) − ω sin ( 2ω t ) , λ sin ( 2ω t ) + ω cos ( 2 ω t ) . }
Substituting (10.17a; 10.18a, b) in (10.16a–c) specifies the exact average over a
period (10.19a) of, respectively. (10.19b/c/d) the potential (10.15b), kinetic (10.15c).
and total energies (10.15d) of a subcritically damped harmonic oscillator:
k a 2 −2 λt
( t, t + τ ) ≡ ( t, t + 2π/ω ) : Φ ( t ) , Ev ( t ) , E ( t ) =
16π
e (
1 − e −4 πλ/ω )
ω λω ω2 ω ωλ
× + 2 cos ( 2ω t ) − 2 sin ( 2ω t ) , + 2 cos ( 2ω t )
λ ω 0 ω 0 λ ω0
ω2 ω 2ω λ
+ sin ( 2ωt ) , 2 + 2 cos ( 2ω t ) .
ω 20 λ ω0
(10.19a–d)
The total energy decays for strong (E10.2.1) or weak (E10.2.2) damping.
Examples 10.1 to 10.20 7
2µ
Ψ ( t ) = − µ x ( t ) = − Ev ( t ) = − 4 λ Ev ( t ) ; (10.20a–c)
2
m
thus, the average of the dissipation over a period (10.19a) is given (10.19c) by
(10.21b):
k a 2 −2 λt
Ψ ( t ) = − 4 λ Ev ( t ) = −
4π
e (
1 − e 4 π λ/ω )
(10.21a, b)
ω λ2 ω2λ
× ω + 2 cos ( 2 ω t ) + 2 sin ( 2ω t ) .
ω0 ω0
The time derivative of the total energy (10.15d) is balanced by the dissipa-
tion (10.20a) at any time in the absence of external forces, that is, A = 0 in
(2.12); thus, the same relation applies to the time derivative of the average
of the total energy over a period (10.19d) that is given (10.20a–c; 10.21a, b) by
(10.22a):
d E ( t ) k a 2 −2 λt
dt
=
4π
e (
1 − e 4 πλ/ω )
(10.22a, b)
ω λ2 ω2 λ
× ω + 2 cos ( 2ω t ) + 2 cos ( 2ω t ) = − ψ ( t ) ,
ω0 ω0
minus the dissipation (10.21b), both averaged over a period (10.22b). In the
case of weak damping (10.23a), substitution of (10.23b) gives (10.23c):
4π λ 4π λ
λ 2 << ω 20 ∼ ω 2 : 1 − e −4 πλ/ω ∼ ∼ . (10.23a–c)
ω ω0
k a 2 −2 λt 4πλ ω 0 λ
Φ (t ) = e + cos ( 2ω 0 t ) − sin ( 2ω 0 t )
16π ω0 λ ω0
(10.24a)
k a 2 −2 λt λ λ2
= e 1 − sin ( 2ω 0 t ) + O 2 ,
4 ω0 ω0
8 Differential Equations and Applications
k a 2 −2 λt 4πλ ω 0 λ
Ev ( t ) = e + cos ( 2ω 0 t ) + sin ( 2ω 0 t )
16π ω 0 λ ω 0
(10.24b)
k a 2 −2 λt λ λ2
= e 1 + sin ( 2ω 0 t ) + O 2 ,
4 ω0 ω0
k a 2 −2 λt 4πλ ω 0 k a 2 −2 λt
λ 2 << ω 20 : E (t ) = e 2 = e . (10.24c, d)
16π ω0 λ 2
The further assumption (10.25a) that is met by (10.25b) leads to the equiparti-
tion of potential energy and kinetic energy (10.25c–e):
ω0 k a 2 −2 λt E (t )
sin ( 2ω 0 t ) << : ω 0 >> λ , e = Φ ( t ) = Ev (T ) = . (10.25a–e)
λ 4 2
Thus, for (problem 10) a subcritically damped oscillator (2.106a–c) omitting the
phase (10.15a): (i)(ii) the instantaneous value (average over a period) are given by
(10.15b/c/d) [(10.19b/c/d)] for the potential/kinetic/total energies and by (10.20a–c)
[(10.21a, b)] for the dissipation; (iii) in the case (problem 11) of weak damping (10.24c)
the total energy (10.24d) decays at twice the damping rate; (iv) with the further
assumption (10.25a) there is equipartition (10.25c–e) of potential and kinetic ener-
gies; (v) the assumption (10.25a) is met by (10.25b), which implies negligible damp-
ing, except in the exponential term exp( −2 λ t ), which is much smaller than unity for
2 t > 1/λ but is close to unity for short time t << 2/λ .
k k f a2 2 m f a2 2
Φ (t ) = x* ( t ) = t sin 2 ( ω 0 t ) = t sin 2 ( ω 0 t ) , (10.26a)
2
2
2 8ω 0 8
m m f a2
sin ( ω 0 t ) + ω 0 t cos ( ω 0 t )
2
Ev ( t ) = x * ( t ) =
2
2
2 8ω0
(10.26b)
m f a2 2
sin ( ω 0 t ) + ω 0 t sin ( 2 ω 0 t ) + ( ω 0 t ) cos 2 ( ω 0 t ) ,
2
=
8 ω 20
m 2 2 t 1
E ( t ) = Φ ( t ) + Ev ( t ) = f a t + sin ( 2 ω 0 t ) + 2 sin 2 ( ω 0 t ) , (10.26c)
8 ω0 ω0
Examples 10.1 to 10.20 9
m f a2
A ( t ) = f ( t ) x * ( t ) = cos ( ω 0 t ) sin ( ω 0 t ) + ω 0 t cos ( ω 0 t )
2 ω0 (10.26d)
m f a2
= sin ( 2 ω 0 t ) + 2 ω 0 t cos 2 ( ω 0 t ) .
4ω 0
The averages over a period starting at time t and ending at time t + τ 0 are
given by:
m f a2
Φ ( t ) , Ev ( t ) , E ( t ) , A ( t ) =
8 ω 20 (10.27a–d)
1
2{ 1
× J − , + J 0 + J + , + J 0 + J + + J − , 4 ω 0 J1 ,
2 }
involving besides (2.116c; 2.117a, b) the integrals:
t +τ 0
{J , J , J } = 1
0 1 ±
τ0 t ∫ {
ω 0 t sin ( 2ω 0 t ) , ω 0 t cos 2 ( ω 0 t ) , ω 20 t 2 cos 2 , sin 2 ( ω 0 t ) dt , }
(10.28a–d)
to be evaluated exactly.
The first integral (10.28a) is evaluated by a single integration by parts:
t +τ 0
1 d
J0 = −
2τ0 ∫t
t cos ( 2ω 0 t ) dt
dt
1 ω 0 t + 2 π/ω 0 (10.29)
∫
t +τ 0
= − t cos ( 2ω 0 t ) t + cos ( 2ω 0 t ) dt
2τ0 4π t
1 1 t +τ 0 1
= − cos ( 2ω 0 t ) + sin ( 2ω 0 t ) t = − cos ( 2ω 0 t ) .
2 8π 2
{ }
t +τ 0
ω 0 2 t +τ0 1 d
= t
4τ 0 t
+
4τ 0 ∫ t
t sin ( 2ω 0 t ) dt
dt
ω0 1 1 t +τ0 (10.30)
∫
t +τ 0
= ( t + τ ) 2
− t 2
+
t sin ( 2 ω t )
− sin ( 2 ω t ) dt
4τ 0 4τ 0
0 0 0
t
4τ 0 t
ω 1 1
= 0 ( 2 t + τ 0 ) + sin ( 2ω 0 t ) +
t +τ 0
cos ( 2ω 0 t ) t
4 4 8τ 0 ω 0
ω t π 1
= 0 + + sin ( 2ω 0 t ) ,
2 2 4
so that only the remaining pair is (10.28c, d).
10 Differential Equations and Applications
t +τ 0
1
J± =
2τ0 ∫t
( ω 0t )2 1 ± cos ( 2ω 0t ) dt
{ }
t +τ 0
ω 20 3 t +τ0 1 d
=
6τ 0
t
t
±
4τ 0 ω 0 ∫t
( ω 0 t )2 sin ( 2ω 0 t ) dt
dt
t +τ 0
ω 20 1 1
( ) ∫
t +τ 0
− t 3 ±
3
= t + τ0 ( ω 0 t ) sin ( 2ω 0 t )
2
ω 0 t sin ( 2ω 0 t ) dt
6τ 0 8π t 2τ0 t
ω 20 1 2 1
6
= (
3 t 2 + 3 t τ 0 + τ 02 ±
8π
)
ω 0 τ 0 ( 2 t + τ 0 ) sin ( 2ω 0 t ) ± cos ( 2ω 0 t ) ;
4
(10.31)
the last term on the r.h.s. of (10.31) is J 0 in (10.28a), which is specified by (10.29),
so that no more integrations by parts are needed:
ω 20 t 2 2π 2 ω 0t π 1
J± = + π ω0 t + ± + sin ( 2ω 0 t ) ± cos ( 2ω 0 t ) , (10.32a, b)
2 3 2 2 4
1 π ω 0t 1
J 0 = − cos ( 2ω 0 t ) , J1 = + + sin ( 2ω 0 t ) , (10.33a, b)
2 2 2 4
m f a2
Φ ( t ) , Ev ( t ) , E ( t ) , , A ( t ) =
16 ω 20
4π 2 1
{
× ω 20 t 2 + 2 π ω 0 t +
3
− ( π + ω 0 t ) sin ( 2ω 0 t ) − cos ( 2ω 0 t ) ,
2
4π 2 1
1 + ω 20 t 2 + 2 π ω 0 t + + ( π + ω 0 t ) sin ( 2ω 0 t ) − cos ( 2ω 0 t ) ,
3 2
8π 2
− cos ( 2ω 0 t ) , 4 π ω 0 + 4ω 20 t + 2ω 0 sin ( 2ω 0 t ) } .
1 + 2ω 20 t 2 + 4π ω 0 t +
3
(10.34a–d)
Examples 10.1 to 10.20 11
The activity of the applied forces (10.34d) equals the rate of change of the
total energy (10.34c):
d d f a2 m
E ( t ) = 2ω 20 t 2 + 4πω 0 t − cos ( 2ω 0 t )
dt dt 16ω 20 (10.35)
2
fa m
= 4 ω 0 t + 4πω 0 + 2ω 0 sin ( 2ω 0 t ) = A ( t ) .
2
16ω 20
If has been proved that an undamped harmonic oscillator (problem 29) corre-
sponds to: (i) the displacement (2.176b) and velocity (2.176c); (ii) the potential energy
(2.177a) ≡ (10.26a), the kinetic-energy (2.177b) ≡ (10.26b), and the total energy
(2.130) ≡ (10.26c), and the activity (2.180a) ≡ (10.26d) of the applied force; (iii) their
exact average over the period from t to t + τ 0 is given exactly by (10.34a–d); (iv) the
terms growing fastest with time t are:
fat
2
E (t ) m f a2 t d
Φ ( t ) = Ev ( t ) = m = , A ( t ) = = E ( t ) , (10.36a–c)
4 2 4 dt
in agreement with (2.178a–c; 2.179b) valid for amplitude varying slowly over a wave
period; (v) within the latter approximation (iv) there is equipartition of the potential
(10.36a) and kinetic (10.36b) energies, whose sum is the total energy (10.36c), grow-
ing like the square of time t due to the work of the applied forces growing linearly
in time.
∂n at
P2 ( D ) = D2 + ω 20 : f0 t n = f0 lim
a→ 0 ∂an
( )
e , (10.39a, b)
and the forcing to (10.37b) ≡ (10.39b). The corresponding response is (problem 39)
obtained using the inverse of the characteristic polynomial:
−1 ∞
D2 n ( − )m D 2 m t n
∑ω
fa n fa fa
λ = 0: x* ( t ) = t = 2 1 + 2
t = 2
P2 ( D ) ω0 ω0 ω0 m= 0
2m
0
(10.40a, b)
≤ n/2
( − )m
∑ω
fa n!
= 2 t n− 2 m .
ω0 m= 0
2m
0 ( n − 2 m)!
fa 1
x* ( t ) = Re t e iω at = Re f a e iω at t
P
2 ( D ) P2 ( D + iω a )
(10.41)
f a e iω at f a e iω at
= Re t = Re 2 t .
( D + iω a ) + ω 0
2 2 2
ω0 − ωa + 2 i ωaD + D
2
This splits into: (i) the non-resonant (10.42a, b) case (problem 40):
fa 2 iω a D + D2
−1
λ = 0, ω 0 ≠ ω a : x* ( t ) = Re e iω at 1 + t
2 2
ω0 − ωa ω 20 − ω 2a
fa 2 i ωaD 2
= 22
ω − ωa
Re e iω at 1 − ω 2 − ω 2 + O D t
( )
0 0 a
(10.42a–c)
fa 2 iω a
= Re e iω at t − 2
2
ω − ωa
0
2
ω 0 − ω 2a
fa 2 ωa
=
ω − ω 2a
2 t cos ( ω a t ) + ω 2 − ω 2 sin ( ω a t ) ,
0 0 a
(problem 41) the resonant case (10.43a, b) does not hold (10.42c) and (10.41)
leads to (10.43c):
f a e iω at f e iω at D
−1
x* ( t ) = Re
a
λ = 0: ω 0 ≡ ω a : t = Re 1 + 2 i ω t
D ( 2 i ω 0 + D ) 2 i ω 0 D 0
i f a e iω at D i f a e iω at 1
= − Re 1 − 2 i ω t = − Re t − 2 i ω
2 ω 0 D 0 2 ω 0 D 0
i f a iω at t 2 t f a t 2 fa t
= − Re e − = sin ( ω 0 t ) + cos ( ω 0 t ) ,
2 ω 0 2 2 i ω 0 4 ω 0 4 ω 20
(10.43a–c)
showing that the amplitude of the displacement is a quadratic function
of time t (10.43c) instead of linear for the forcing (10.37c) and the highest
( )
order term O t 2 in (10.43c) is out-of-phase by π/2 relative to the forcing
(10.37c), that is, O ( t ). Setting (10.44b) in (10.42c) leads to (10.44d):
fa t
λ = 0 , ω a = 0 , n = 1: x* ( t ) = , (10.44a–d)
ω 20
which is the same result (10.44d) ≡ (10.40b) for (10.44c). It has been shown that
the undamped (10.37a) harmonic (10.37b) oscillator with forcings (10.37b) [(10.37c)]
has displacements (10.40a, b) [(10.42a, c)/(10.43a, c)], respectively, in the [non-
resonant/resonant cases]. The non-resonant (resonant) cases differ in that the ampli-
tude of the displacement is a function of time t with the same exponent (exponent one
unit higher) than the forcing and the highest powers are in-phase (out-of-phase by π/2).
fa f a e iω at
x* ( t ) = Re t e iω at = Re t
P
2 ( D ) P
2 ( D + iω a )
f a e iω at
= Re 2 t
ω 0 + 2 λ ( D + iω a ) + ( D + iω a )
2
f a e iω at
= Re 2 t
ω 0 − ω a + 2 iλω a + 2 ( λ + iω a ) D + D
2 2
f a e iω at 2 ( λ + iω a ) D + D2
−1
= Re 2 2 1 + 2 2 t
ω 0 − ω a + 2 iλω a ω 0 − ω a + 2 iλω a (10.47)
2 ( λ + iω a ) D
ω 20 − ω 2a − 2 iλω a
= Re f a e iω at 1− 2 t
( )
2
ω 20 − ω 2a + 4λ 2 ω 2a ω 0 − ω 2a + 2 iλω a
= fa t
(ω 2
0 )
− ω 2a cos ( ω a t ) + 2 λω a sin ( ω a t )
(ω )
2 2
0 − ω 2a + 4λ 2 ω 2a
−
(
ω 20 − ω
2 fa
)
2 2
+ 4λ ω 2 2
2 {
(
λ ω 20 + ω 2a )
2
(
− 4 ω 2a ω 2a + λ 2
)
a
a
(
cos ( ω a t ) − ω a ω 20 − ω 20
)
2
+ 4λ 2 ω 20 sin ( ω a t ) .
}
In the last term on the r.h.s. of (10.47) was used:
( λ + iω a ) (ω 20 − ω 2a − 2iλω a ) = ( λ + iω a ) ω 20 − ω 2a( ) (
− 4λ 2 ω 2a − 4iλω a ω 20 − ω 2a )
2 2
(
= λ ω 20 − ω 2a ) (
− 4λ 2 ω 2a + 4ω 2a ω 20 − ω 2a )
2
(
+ i ω a ω 20 − ω 2a ) (
− 4λ 2 ω 2a − 4λ 2 ω 20 − ω 2a )
2
(
= λ ω 20 + ω 2a ) (
− 4ω 2a ω 2a + λ 2 )
2
+ iω a ω 20 − ω 2a − 4λ 2 ω 20 . ( )
2
(10.48)
Examples 10.1 to 10.20 15
fa t f a cos ( ω 0 t ) sin ( ω 0 t )
ωa = ω0 : x* ( t ) = sin ( ω 0 t ) + + . (10.49a, b)
2 λ ω0 2 λω 0 ω 0 λ
y C cos x
LII: y ′ = − + cos x , y ( x) = + sin x + , (10.51a, b)
x x x
y y2 x
LIII, LXIV: y′ = + , y ( x ) = , (10.52a, b)
x x2 C − log x
y 2 − x2
LXIV: y′ = , y ( x ) = x C − 2 log x , (10.53a, b)
xy
−1
1 1 ex x
eξ
LV:
x
2
y ′ = y − y − ; y ( x) = + 2 C +
x x ∫ dξ , (10.54a, b)
ξ 2
and use the appropriate method to obtain the respective general integrals
(10.50b–10.54b), where C is an arbitrary constant.
The first differential equation (10.50a) is (subsection 3.2.1), separable (standard L)
and thus, (10.55a) integrable by quadratures leading to (10.55b) ≡ (10.50b):
dy
∫ y = ∫ log x dx, log y = log C + x log x − x. (10.55a, b)
16 Differential Equations and Applications
1 x
dξ 1 (10.56a–c)
P=− ,
x
Q = cos x: X ( x ) = exp −
∫ ξ
= exp ( − log x ) = ,
x
x x
∫ ∫
Y ( x) = ξ cos ξ dξ = x sin x − sin ξ dξ = x sin x + cos x , (10.56d)
1 1 x dξ 1 x
dξ
P = , Q = 2 : X 2 ( x ) = exp −
x x ξ x ∫
= , Y2 ( x ) = −
ξ
= − log x ;
∫
(10.57a–d)
1 1
n = 2: = ( C − log x ) . (10.58a, b)
y ( x) x
y/x
dv C x
f ( v ) = v + v2 : x = e C exp
∫ v 2
= e exp − ,
y
(10.59a, b)
1 y /x
y2
∫ v dv = exp − 2 , (10.60a, b)
f ( v) = v − : e − C/2 x = exp −
v 2x
Examples 10.1 to 10.20 17
which when solved for y gives (10.60b). The fifth differential equation (10.54a)
is of the Riccati type (standard LV), with (3.58) coefficients (10.61a–c) and a
particular integral is (10.61d) leading (3.63a–d) to (10.61e–g):
1 1
Q = − , P = 1, R = − 1 : g1 ( x ) = ,
x x
(10.61a–g)
x
2 ex 1 x
eξ
Z ( x ) = exp
∫ 1 −
ξ
dξ =
x 2
=
X ( x)
, Y ( x) =
∫ ξ2
dξ ;
0 = dy − ( Py + Q ) dx = A dx + B dy:
∂B ∂ A (10.62a–g)
{ A , B} = {− Py − Q, 1 } , Ω≡ −
∂x ∂y
= − P ≠ 0,
the condition (3.183a) is met (10.63a, b) and thus, there exists (3.183b) an
integrating factor (10.62c, d) depending only on x whose differential is
(10.62e):
Ω x
1
−
B
= −P ≡ f ( x) : λ ( x ) = exp −
∫ P ( ξ ) dξ ≡
X ( x)
, dλ = − P λ dx.
(10.63a–e)
The integrating factor (10.63c) is the same as for the linear unforced first-
order differential equation (3.19a) because it does not depend on Q. The
general integral depends on Q and is obtained multiplying the inexact
18 Differential Equations and Applications
differential (10.62a) by the integrating factor (10.62b) and using its differential
(10.63c) leading to an exact differential (10.64a, b):
0 = λ dy − λ ( Py + Q ) dx = λ dy + y dλ − λ Q dx
(10.64a, b)
x
= d λ y −
∫ λ ( ξ ) Q ( ξ ) dξ ,
C = λ ( x) y( x) −
∫ λ ( ξ ) Q ( ξ ) dξ
(10.65a, b)
y ( x) x Q (ξ)
=
X ( x)
−
∫ X (ξ)
dξ ,
LXVI: (x 2
)
+ y 2 dx + y ( y + 2 x ) dy = 0 , x 3 + y 3 + 3 x y 2 = C, (10.66a, b)
x x2
LXVI : dx − dy = 0, x 2 = 2 C y 2 , (10.67a, b)
y2 y3
2xy
LXX: dx + ( x − y ) dy = 0, 3 x 2 y − y 3 = C , (10.71a, b)
x+y
2x
LXXII: 3 dx − dy = 0 , y 2 = C x 3 , (10.73a, b)
y
Examples 10.1 to 10.20 19
∂X ∂Y
{X , Y } = {x 2 + y 2 , y 2 + 2 x y } : = 2y = ,
∂y ∂y
(10.74a–c)
x3 + y 3
0 = x 2 dx + y 2 dy + y 2 dx + 2 x y dy = d + x y2 ,
3
x x2 ∂X 2 x ∂Y x dx x 2 x2
{X , Y } =
2
, − 3 : =− 3 = , 0 = 2 − 3 dy = d 2 ,
y y ∂y y ∂x y y 2y
(10.75a–c)
∂X ∂Y
{X , Y } = {cos x cos y , − sin x sin y} : = − cos x sin y =
∂y ∂x (10.76a–c)
0 = cos x cos y dx − sin x sin y dy = d ( sin x cos y ) ,
1
{X , Y } = {2 x y , − x 2 } : Ω = − 4 x, λ ( y ) = : (10.78a–c)
y2
2 x y dx − x 2 dy 2 x x2 x2
0= 2
= dx − 2 dy = d ; (10.78d)
y y y y
20 Differential Equations and Applications
2xy 2 x2
{X , Y } = , x − y : Ω = 1− , λ ( x , y ) = x + y: (10.79a–c)
x + y ( x + y )2
2xy y2
0 = (x + y)
dx + ( x − y ) dy = 2 x y dy + x 2 dy − y 2 dy = d x 2 y − ;
x+ y 3
(10.79d)
{X , Y } = {3 x y , 2 x 2 } : Ω = x , λ ( x y ) = x y: (10.80a–c)
( ) (
0 = x y 3 x y dx + 2 x 2 dy = 3 x 2 y 2 dx + 2 x 3 y dy = d x 3 y 2 ; (10.80d) )
2x 2 x2
{X , Y } = 3, − : Ω=− , λ ( x , y ) = 2 : (10.81a–c)
y y y
x2 2x 3 x2 2x3 x3
0= 3 dx − dy = dx − dy = d y 2 ; (10.81d)
y2 y y2 y3
Ω 2 y
dξ 1
∫
(10.82a, b)
= − = g ( y ): λ ( y ) = exp −2 = .
X y ξ y2
The inexact differential (10.71a) with coefficients (10.79a, b) and curl (10.79c)
meets (standard LXX) the condition (3.196a) for the existence (10.83a) of an
integrating factor (10.83b) ≡ (10.79d) depending on x + y:
Ω 1 x+ y
dξ
∫
= = h ( x + y ): λ ( x + y ) = exp = x + y. (10.83a, b)
X −Y x+ y ξ
Examples 10.1 to 10.20 21
The inexact differential (10.72a) with coefficients (10.80a, b) and curl (10.80c)
meets (standard LXXI) the condition (3.201a) for the existence (10.84a) of an
integrating factor (10.84b) ≡ (10.80d) depending on xy:
Ω 1 xy
dξ
∫
(10.84a, b)
= = j ( x y ): λ ( x y ) = exp = x y.
X x −Yy xy ξ
The inexact differential (10.73a) with coefficients (10.81a, b) and curl (10.81c)
meets (standard LXXII) the condition (3.206a) for the existence (10.85a) of an
integrating factor (1–0.85b) ≡ (10.81d) depending on y/x:
x2 Ω x y y y/x
dξ x 2
Xx+Yy
= −2 = k :
y x
λ = exp −2
x ∫ =
ξ y2
. (10.85a, b)
The general integrals of the five inexact differentials in two variables are
obtained using an integrating factor that depends only on: (i) one vari-
able (10.70a; 10.82a, b); (ii–iv) the sum (10.71a; 10.83a, b), the product (10.72a;
10.84a, b), and the ratio (10.73a; 10.85a, b) of two variables; (v) the two variables
(10.70a; 10.77c) in a different way from the preceding (ii–iv).
dx dy
= vx ( x , y ) = x y , = vy ( x , y ) = y 2 − x 2 , (10.86a, b)
dt dt
and determine (problem 77) the paths near the stagnation point.
The paths are specified by the ratio of (10.87a), leading to a first-order ordi-
nary differential equation of degree two, which is not included in the cases
of first degree in section 4.1:
dy y 2 − x 2 y x y 1
= = − = f : f ( v) = v − ; (10.87a, b)
dx xy x y x v
22 Differential Equations and Applications
(
y 2 = x 2 C 2 − 2 log x . )
The integral (10.87d) curves (Figure 10.1). (i) are symmetric relative to
the y-axis (10.88b) and exist in the interval with boundaries (10.88a) and
(10.88 c, d); (ii) pass through the origin (10.88e) and cross (10.88f) the vertical
line x = 1 at y = ± c:
C2
0 ≤ x ( y ) = x ( − y ) ≤ exp ≡ xmax , y ( 0 ) = 0, y (1) = ± C , (10.88a–f)
2
y2 dy
lim = lim ( C − log x ) = ∞ , lim = 0,
x→ 0 x 2 x→ 0 x→± y dx
(10.88g–j)
dy 1
lim = ±C ,
x→ 1 dx C
y
√2
1 C = √2
1/√2
1 C=1
C=
√2
1
4
√e √e e x
–1/√2
–1
–√2
FIGURE 10.1
An example of singular point of a family of curves or streamlines is the point through which
all streamlines pass with the same tangent direction.
Examples 10.1 to 10.20 23
(iii) the slope is infinite at the origin (10.88g, h) where the tangent is ver-
tical, is zero (10.88i) at the crossing of the diagonal of quadrants and is
(10.88j) at the crossing of the vertical line x = 1. For example the integral
curve: (a) C = 1 crosses the x-axis at the ends of the interval 0 ≤ x ≤ e = 1.649
and has zero slope at x = 1, y = ±1; (b) C 2 crosses the x-axis at the ends
of the interval 0 ≤ x ≤ e = 2.718 and has slope ± 2 1/ 2 at the points
x = 1, y = ± 2 ; (c) C = 1/ 2 crosses the x-axis at the ends of the interval
0 ≤ x ≤ 4 e = 1.284 and has slope ±1/ 2 2 with opposite sign to the case
(b) at x = 1, y = ±1/ 2 .
dx dy
= vx ( x , y ) = − 2 x y , = vy ( x , y ) = y 2 − x 2 , (10.89a, b)
dt dt
dy 1 x 2 − y 2 1 x y y 1 1
= = − = f , f ( v) = − v . (10.90a, b)
dx 2 xy 2 y x x 2v
1 y/x
v 1 y2 C
log x −
3
log C = 2
∫ 1− 3v
dv = − log 1 − 3 2 ; x 2 − 3 y 2 = .
3 x x
(10.91a, b)
C
x x →∞ x→∞
( x →−0
1
3 )(
C
x x→ 0
)
= lim x − y 3 x + y 3 , ∞ = lim x 2 − = lim y 2 ,
0 = lim
(10.91c–f)
along: (i) the straight lines (10.91c, d) with slopes ±1/ 3 making angles ± π/6
with the x-axis; (ii) along the y-axis (10.91e, f) since y → ±∞ as x → −0. Thus the
asymptotes divide the plane into six angular sectors for ϕ = ± π/6, ± π/2, ± 5π/6
that are the only curves passing through the origin (Figure 10.2).
24 Differential Equations and Applications
π/2
5π/ 6 π/6
7π/ 6 11π / 6
3π / 2
FIGURE 10.2
A different singular point from that in Figure 10.1 is the point through which only six stream-
lines pass, which are the asymptotes of all other streamlines.
( ) ( )
∂ y 2 − x 2 / ∂ x = − 2 x = ∂ ( −2 x y )/y , ∂( −2 x y )/ ∂ x + ∂ y 2 − x 2 / ∂ y = 0.
(10.93a, b)
The latter velocity field (10.89a, b) corresponds to the potential flow (section
I.12.19) of complex conjugate velocity:
( )
f ′ = vx − i vy = − 2 x y − i y 2 − x 2 = i ( x + i y ) = i z 2 ;
2
(10.94a, b)
Examples 10.1 to 10.20 25
the stream function Ψ, which is the imaginary part of the complex potential:
C i z3 i ( x + i y )3 x3
2
= Im f = Im = Im = − x y + . (10.95a–d)
3 3 3 3
y ′3 C + p 2 /2 − p 3 /3
C: y = y ′ 2 x + ; x ( p; C) = , (10.97a–c)
3 (1 − p )−2
y x2
CII: x = + y ′, y ( x ; C ) = C ( x − C ) ; y ( x ) = , (10.98a–c)
y′ 4
CIII: y ′ 3 − y ′ 2 y − 4 x 2 y ′ + 4 x 2 y = 0, (y − x 2
)( )( )
+ C y + x 2 + C y − C e x = 0,
(10.99a, b)
and obtain the general integrals (10.96b–10.99b) and the special integrals
when they exist (10.96c–10.98c). The roman numerals at left identify the stan-
dard in the classification of differential equations in the recollection 10.1.
1 1 ∂y 1 1
h( p) = : y ( x; C) = C x + : = x − 2 = 0, C = , y ( x) = 2 x ,
p C ∂C C x
(10.100a–c)
26 Differential Equations and Applications
y
y = C x + 1/C
y = 2 √−x y = 2 √x
envelope
1/C
x
−b −a a b
FIGURE 10.3
The cusped parabola is the envelope of a family of straight lines whose ordinate at the origin
is the inverse of the slope.
p3 2 1
g ( p ) = p2 , h( p) =
3
: ϕ ( p ) = exp
∫ 1 − p = (1 − p) 2 , (10.101a–d)
x ( p ; C) =
1
(1 − p ) 2 { ∫
C + p ( 1 − p ) dp = } C + p 2 /2 − p 3 /3
( 1 − p )2
. (10.102)
y dp
x= + p ≡ g ( y, p ): = 0; p = C, y ( x ; C ) = C ( x − C ) , (10.103a–d)
p dy
Examples 10.1 to 10.20 27
y = x 2 /4
y = Cx − C 2
− C2
FIGURE 10.4
A family of straight lines distinct from that in Figure 10.3 is the case of ordinate at the origin
equal to minus the square of the slope that has a smooth parabola as an envelope.
hence (10.103c) and the general integral is (10.103d) = (10.98b). The C-discriminant
(10.104a) leads (10.104b) to the special integral (10.104c) ≡ (10.98c):
∂y x x2
0= = x − 2C, C= , y ( x) = . (10.104a–c)
∂C 2 4
The integral curves are straight lines (10.103d) of slope C and ordinate −C 2 at the
origin, tangent to the smooth parabola (10.104c), which is their envelope (Figure 10.4).
( )
0 = y ′ 3 − y ′ 2 y − 4 x 2 y ′ + 4 x 2 y = ( y ′ − y ) y ′ 2 − 4 x 2 = ( y ′ − y )( y ′ − 2 x )( y ′ + 2 x ) ;
(10.105)
there are three solutions (10.106a–c) of (10.105):
y1− 3 ( x ) = C e x , ± x 2 − C , 0 = y − y1 ( x ) y − y 2 ( x ) y − y 3 ( x ) , (10.106a–d)
1
CXI: y ′′ = , y ( x ) = C2 + C1 x + x log x − x , (10.107a, b)
x
x2
I, CX: y ′′ = y ′ + x: y ( x ) = C2 + C1 e x − − x − 1, (10.108a, b)
2
1
3 y ( x ) = C2 + 2 ( C1 + x )
3/2
CVIII: y ′′ = , , (10.110a, b)
y′
3 2
CVIII: y ′′ =
2
y , x ( y ) = C2 +
∫ (C + y )
1
3 −1/2
dy , (10.111a, b)
2y
CVIII: y ′′ =
3 y′
: x ( y ) = C2 +
∫ (C + y ) 1
3 −1/3
dy , (10.112a, b)
x ( y ) = C2 +
∫ 2 (C sinh y )
−1/2
CVIII: y ′′ = − cosh y: 1 dy , (10.113a, b)
y′ dy
CVIII: y ′′ =
y
: y ( x ) = C2 +
∫ C + log y ,
1
(10.114a, b)
1
15 y ( x ) = C3 + C2 x + 4 ( C1 + x )
5/2
CXIV: y ′′′ = : , (10.115a, b)
2 y ′′
CXVI: y ′′′ = − 2 y ′′ y ′: x ( y ) = C3 +
∫ C
− 1/2
2 + log ( C1 + x ) dy , (10.116a, b)
y ′′ x2
CXIII: y ′′′ = 2
: y ( x ) = C3 + C2 x + − C1−3 ( 1 + C1 x ) log ( 1 + C1 x ) − 1 ,
x 2C1
(10.117a, b)
N −2
x a+ N −1
CXIII, L: x y( N ) ( x ) = a y( N − 1) ( x ) : y ( x ) = ∑
n= 0
Cn+ 1 x n + CN
( a + 1)( a + N − 1)
,
(10.118a, b)
Examples 10.1 to 10.20 29
1 e x /2
2
CXVII: y ′′ + x +
x
y ′ + 2 y = 0: y ( x ) = e
− x 2/2
C2 + C1
∫ x
dx ,
(10.119a, b)
CXVIII: x y ′′ y ± y ′y + y ′ 2 x = 0: y + ( x ) = C2 + C1 log x ,
(10.120a–c)
y − ( x ) = C2 x 2 − C1 ,
(
CXIX: y ′′ y − y ′ 2 + 6 x y 2 = 0: y ( x ) = C2 exp x C1 − x 2 , ) (10.121a, b)
The second-order linear differential equation (10.108a) misses both the inde-
pendent and dependent variables (standard CX) and is also linear with
constant coefficients (1.197a–c) and hence has: (i) characteristic polynomial
(10.124a) with (standard I) distinct roots zero and unity corresponding to the
general integral in the first two terms on the r.h.s. of (10.124b):
1
0 = P2 ( D ) = D2 − D = D ( D − 1) , y ( x ) = C2 + C1 e x + x;
D ( D − 1)
(10.124a, b)
(ii) the last term on the r.h.s. of (10.124b) is (standard XIII) the particular inte-
gral (10.124c–e):
1 1 x2 x2 x2
D ( D − 1)
x=
D−1 2
= − 1 + D + D2 + O D2
2
=−
2
( )
− x − 1;
(10.124c–e)
30 Differential Equations and Applications
(iii) sum of (10.124b) and (10.124e) is the complete integral (10.108b). The lin-
ear second-order differential equation with variable coefficients (10.109a) ≡
(10.125a) is a separable first-order differential equation (10.125b) in y ′:
y ′′ d d
= − tan x: ( log y ′ ) = log ( cos x ) , (10.125a–d)
y′ dx dx
y ′ 2 = 2 ( x + C1 ) , y ′ = 2 ( x + C1 )
1/2
y ′dy ′ = dx , ; (10.126a–c)
(ii) a further integration of (10.126c) leads to the general integral (10.110b). The
second-order differential equation (10.111a) is (standard CVIII) also equiva-
lent (5.138a) to a first-order equation (10.127a) for y ′:
3 2
( )
−1/2
y ′dy ′ = y dy , y ′ 2 = y 3 + C1 , dx = y 3 + C1 dy , (10.127a–c)
2
( )
−1/3
3 y ′ 2 dy ′ = 2 y d y: y ′ 3 = y 2 + C1 , dx = y 2 + C1 dy , (10.128a–c)
leading to (10.128b, c) and hence, to the general integral (10.112b). The differ-
ential equation (10.113a) is (standard CVIII) reducible to first-order and leads
(5.138a) to (10.129a) whose integration yields (10.129b, c):
and hence to the general integral (10.113b). The non-linear differential equa-
tion (10.114a) omits the independent variable (standard CVIII) and leads to a
first-order differential equation (10.130a) of degree two:
dy ′ 1 dy
0 = y′ − : y ′ = 0: y = C0 ; dy ′ = , (10.130a–d)
dy y y
which has: (i) a special integral (10.130b) involving only one constant of inte-
gration (10.130c); (ii) the second factor (10.130d) in (10.130a) leads (10.131a) to
the general integral (10.131b) ≡ (10.114b):
dy
y ′ = C1 + log y , dx = . (10.131a, b)
C1 + log y
The differential equations reducible to the first order may be of the second
(higher) order [E10.9.2 (E10.9.3)].
2 y ′′ dy ′′ = dx , ( y ′′ )2 = x + C1 , y ′′ = x + C1 , (10.132a–c)
dy ′′ dy ′′
y′ + 2 y ′′ 2 = 0: y′ = 0 ; y ( x ) = C0 ; = − dy , (10.133a–d)
dy 2 y ′′ 2
and hence, has: (i) a special integral (10.133b) involving only one constant of
integration (10.133c); (ii) the second factor in (10.133a) leads to (10.133d), which
may be integrated as (10.134a) ≡ (10.134b), implying (10.134c) by a second
integration:
1 dy
y + C1 = : 2 y ′dy ′ = ; y ′ 2 = C2 + log ( C1 + y ) , (10.134a–c)
2 y ′′ y + C1
and yielding the general integral (10.116b) by a third integration. The third-
order differential equation (10.117a) is (standard CXVI) non-linear with
32 Differential Equations and Applications
dy ′′ dx 1 1 x 1 1
= : = C1 + , y ′′ = = 1− ; (10.135a–d)
y ′′ 2 x 2 y ′′ x 1 + C1 x C1 1 + C1 x
x
y ′ = C2 + − C1−2 log ( 1 + C1 x ) ; (10.136)
C1
(iii) a third integration leads to the general integral (10.117b). The differen-
tial equation (10.118a) of order N can be solved (standard CXIII) as a linear
unforced differential equation with constant coefficients or as a first-order
separable (standard L) equation (10.137b) for the derivative (10.137a) of order
N − 1, which may be: (i) integrated once (10.137c, d):
d 1 d d 1 1
0= +
dx x dx
+ x y =
+
dx x
( y ′ + x y ) = y ′′ + x + y′ + 2 y ;
x
(10.138a, b)
since the operators have variable coefficients they generally do not commute,
that is, (10.138a, b) is distinct from:
d d 1 d y
0= + x + y = + x y ′ +
dx dx x dx x
(10.139a, b)
1 1
= y ′′ + x + y ′ + y 1 − 2 ,
x x
Examples 10.1 to 10.20 33
with commutator:
d 1 d d d 1 1
+ + x − + x + y = 1 + 2 y . (10.140)
dx x dx dx
dx x x
The second-order ordinary differential equation (10.119a): (i) splits into two
first-order equations (10.141a, b):
z dz dx C1
z = y ′ + x y: z′ + = 0; + = 0, z= ; (10.141a–d)
x z x x
(ii) the equation (10.141b) ≡ (10.141c) has the integral (10.141d); (iii) substitu-
tion of (10.141d) in (10.141a) leads to a linear differential equation (10.142a) ≡
(10.142b) ≡ (3.27) with coefficients (10.142c; 10.141e):
C1
y′ = −x y + = P y + Q: P ( x ) = − x ,
x
2
(10.142a–f)
x
x2 x
e ξ /2
X ( x) = exp −
∫ ξ dξ = exp −
2
Y ( x ) = C1
∫ ξ
dξ ,
0 = x y ′′y + y ′ 2 x + y ′y = ( x y ′y )′ , ( )
C1 = 2 x y ′y = x y 2 ′ , (10.143a, b)
( y )′ = Cx ,
2 1
y 2 = C1 log x + C2 , (10.144a, b)
34 Differential Equations and Applications
x
( )
0 = x y ′′y + y ′ 2 x + y ′y − 2 y ′y = x y ′y − y 2 ′ , C1 = x y ′y − y 2 = y 2 ′ − y 2 ,
2
( )
(10.145a, b)
leading to: (i) the first integral (10.145b); (ii) the latter is a non-linear first-order
differential equation, which becomes linear (10.145b) ≡ (10.146b) for the vari-
able (10.146a):
z ≡ y2 :
x
2
z ′ − z = C1 ; z′ = P z + Q, { P , Q} ≡ { 2 2 C1
,
x x }
. (10.146a–e)
X x = exp 2
x x
2 C1 C
( ) ∫ dξ = x 2 ,
ξ
Y ( x) =
∫ ξ3
dξ = − 21 ; (10.147a–c)
x
z ( x ) = X ( x ) C2 + Y ( x ) = x 2 C2 − C1 = y ( x ) ,
2
(10.147e, f)
2
y′ y ′′ y ′′ y ′
v≡ , = v′ + v2 : 0= − + 6 x = v ′ + 6 x ; (10.148a–d)
y y y y
y′
−3 x 2 = v − C1 = − C1 , − x 3 = log y − C1 x − log C2 , (10.149a, b)
y
Examples 10.1 to 10.20 35
dw dw
y ′ x ≡ w, x 2 y ′′ = w − 1 : 0 = x 2 y ′′ + x 2 y ′ 2 + y ′x = w + w ,
dy dy
(10.150a–d)
which is a first-order equation (10.150d) for the dependent variable
(10.150a). The special integral solution w = 0 of (10.150d) is y = const by
(10.150a) and satisfies (10.122a). The general integral solution of (10.150d)
is (10.151a):
dw dy C
= − w: C1 e − y = w = x , 0 = 1 dx − e y dy = C1 log x − e y + C2 ,
dy dx x
(10.151a–c)
leading to (10.151b) and hence (10.150c), to the general integral (10.151d) ≡
(10.122b).
in the following cases: (a) a clamped beam (10.153a–d) with (Figure 10.5a), a
concentrated torque Q in the middle (10.153e); (b) a pinned beam (10.154a–d)
with (Figure 10.5b), a concentrated force F in the middle (10.154e); (c) a
clamped-pinned beam (10.155a–d) under (Figure 10.5c) its own weight
(10.155e); (d) a cantilever beam (10.156a–d) subject (Figure 10.5d) to a shear
stress (10.156e) increasing linearly from the clamped end:
ζ ( 0 ) = 0 = ζ′ ( 0 ) ; ζ ( L ) = 0 = ζ′ ( L ) : f ( x ) = Q δ ′ ( x − L/2 ) , (10.153a–e)
ζ ( 0 ) = 0 = ζ′ ( 0 ) ; ζ ( L ) = 0 = ζ′′ ( L ) : f ( x ) = ρ g , (10.155a–e)
b
F
d
qx
FIGURE 10.5
The linear bending of a beam under traction can be considered for 4 × 4 = 16 combina-
tions of support and loads, of which four are illustrated: (i) for a concentrated torque
(a) [transverse force (b)] or a constant (c) [linearly increasing (d)] shear stress; (ii) for
supports clamped (a) [pinned (b)] at both ends or clamped at one end and pinned (c) [free
(d)] at the other end.
d2 d
T d
f ( x ) = 0, p 2 =
0 = ζ′′′′ − p 2ζ′′ = 2
> 0: − p + p ζ ( x ).
EI
dx dx dx
(10.157a–d)
Q
− = E I ζ1′′( 0 ) ; ζ1 ( 0 ) = 0 = ζ1′ ( 0 ) , (10.159a–c)
4
which specifies the curvature, and the coordinate and slope are zero (10.159b, c).
The shape of the beam is skew-symmetric relative to the middle (10.160a)
L L
ζ1 ( x ) = − ζ1 ( L − x ) , ζ1 = − ζ1 = 0, (10.160a, b)
2 2
{ } {
A + C = 0 = B + p D: ζ1 ( x ) = A 1 − cosh ( px ) + D sinh ( px ) − px , }
(10.161a–c)
38 Differential Equations and Applications
from four in (10.158) to two in (10.161c). The remaining two constants of inte-
gration are determined substituting (10.161c) in (10.159a) [(10.160b)] leading to
[(10.162a) [(10.162c)]:
Q/4 = EI A p 2 = AT , { } {
D sinh ( pL/2 ) − pL/2 = A cosh ( pL/2 ) − 1 , }
(10.162a–c)
and using (10.157b) in (10.162b). Solving (10.162b, c) for A and D, and substitut-
ing into (10.161c) yields:
L 4T cosh ( pL/2 ) − 1
0≤x≤ : ζ1 ( x ) = 1 − cosh ( px ) + sinh ( px ) − px ,
2 Q sinh ( pL/2 ) − pL/2
(10.163a, b)
the shape (10.163b) of the first half (10.163a) of the elastica of the beam under
traction (10.157b).
If the traction force is weak in the sense (10.164a), then (III.7.6b) ≡ (10.146b)
and (III.7.7b) ≡ (10.146c):
6 6
pL T L2 p2 x2 p4 x4
( p x ) ≤ 2 = 2 E I << 1:
6
cosh ( px ) − 1 =
2
+
24
(
+ O p6 x6 )
p3 x3 p5 x5
sinh ( px ) − px =
6
+
120
(
+ O p7 x7 , )
(10.164a–c)
imply (10.164d–f):
cosh ( pL/2 ) − 1
−1
3 p2 x2 p2 x2
= 1 +
sinh ( pL/2 ) − pL/2 p x 12
1 +
20
(
+ O p4 x4
)
3 1 1 3 px
p x
1 + p2 x2 −
12 20
+ O p4 x4 =
= +
p x 10
+ O p3 x3 ; ( ) ( )
(10.164d–f)
p 2 x 2 px 6 pL
pL
4
4T p2 x2
( px )4 ≤
1 +
12
− << 1:
2
+ ;
3 pL 20 Q
ζ1 ( x ) = −
2
(10.165a–c)
Examples 10.1 to 10.20 39
2
4EI Q x2 2x Tx L
T 2 << 2 : ζ1 ( x ) = − 1 − L + 12EI x − , (10.166a, b)
L 8EI 5
specifies the shape of the elástica (10.166b) for weak axial traction (10.166a).
The beam equation (10.157c) simplifies to the bar equation (10.167b) in the
absence of axial tension (10.167a):
T = 0: ζ′′′′
0 = 0; ζ0 ( x ) = A + B x + C x 2 + D x 3 , (10.167a–c)
Q
0 = ζ01 ( 0 ) = A, 0 = ζ′01 ( 0 ) = B , − = ζ′′01 ( 0 ) = 2 C ,
4EI
(10.168a–d)
L L
2
DL
0 = ζ01 = C + 2 ;
2 4
Q 2C Q Q x2 2x
C=− , D=− = : ζ10 ( x ) = − 1− , (10.169a–c)
8EI L 4EIL 8EI L
which coincides with (10.166b) ≡ (III.4.220b) for zero tension (10.167a). Thus,
the (problem 222) linear bending (10.152b) of a uniform beam (6.10a, b) subject
to a traction (10.152a) clamped at both ends (10.153a–d) under a concentrated
torque (10.153e) leads to: (i) the transverse displacement of the first (10.163a, b)
[second(10.170a, b)] half of the elastica, which is skew-symmetric (10.160a):
L 4T
≤ x ≤ L: ζ1 ( x ) = cosh p ( L − x ) − 1
2 Q
(10.170a, b)
cosh ( pL/2 ) − 1
− sinh p ( L − x ) − p ( L − x ) ,
sinh ( pL/2 ) − pL/2
40 Differential Equations and Applications
confirming that the displacement is zero at the clamped ends (10.153a, c) and in the
middle (10.160b); (ii) the slope (10.171b) is symmetric (10.171a):
Qp cosh ( pL/2 ) − 1
ζ1′ ( L − x ) = ζ1′ ( x ) = − sinh ( px ) + cosh ( px ) − 1 ,
4T sinh ( pL/2 ) − pL/2
(10.171a, b)
confirming it vanishes at the clamped end (10.153b, d); (iii) the curvature (10.172c)
is skew-symmetric (10.172b) and specifies the bending moment (10.172a) and
implies that M1 ( 0 ) = −Q/4 in agreement (10.159a) with the bending moment at
the support:
The shape of the elastica simplifies to (10.169c) for a bar in the absence of axial tension
(10.167a); the lowest-order correction for a beam in the case of weak axial traction
(10.166a) is (10.166b).
P L
− 2 ( 0 ) − T ζ′ ( 0 ) ,
= E I ζ′′′ ζ′2 = 0, (10.173a, b)
2 2
P pL
− = − T B, 0 = B + p D cosh , (10.176a, b)
2 2
the two remaining constants of integration, and hence:
L 2T 1 pL
0≤x≤ : ζ2 ( x ) = x − sinh ( px ) sech , (10.177a, b)
2 P p 2
6 −1
pL p2 x2 p4 x4
( px ) ≤ << 1: sech ( px ) = 1 +
6
+
2 2 24
2
p2 x2 p4 x4 p2 x2 p2 x2 5p4 x4
= 1 − − + = 1− + ;
2 24 2 2 24
(10.178a–d)
2T p2 x2 p 2 L2 5 p 4 L4
ζ2 ( x ) = 1 − 1 + 1 − +
Px 6 8 384
(10.179a)
2 2 2 2 4 4 4 2 2
p x p L 5p L p Lx
=− + − + ,
6 8 384 48
P x 2 3L2 T L2 2 5L2
2
4E I
T 2 << 2 : ζ02 ( x ) = − x − − x − 8 , (10.179b, c)
L 12 EI 4 8EI
42 Differential Equations and Applications
under weak traction (10.179b). In the absence of axial tension (10.167a) the
beam (10.157c) becomes a bar (10.167b) and the constants in the shape (10.167c)
of the elastica are determined by the boundary conditions (10.175a, b;
10.173a, b) ≡ (10.180a–d):
P
0 = ζ02 ( 0 ) = A, 0 = ζ′′02 ( 0 ) = 2C , − 02 ( 0 ) = 6 D,
= ζ′′′
2EI
(10.180a–d)
L 3DL2
0 = ζ′02 = B + ;
2 4
P 3DL2 P L2 P 2 3L2
D=− , B=− = : ζ20 ( x ) = − x − , (10.181a–c)
12EI 4 16EI 12 EI 4
which coincides with (10.179c) ≡ (III.4.114b) for zero axial tension (10.167a).
In conclusion, (problem 223) the linear bending (10.152b) of a uniform beam
(6.10a, b) pinned at both ends (10.154a–d) under traction (10.152a) with a concen-
trated force (10.154e) at the middle (Figure 10.5b) leads to: (i) the transverse dis-
placement of the first (10.177a, b) [second (10.182a, b)] half of the elastica, which is
symmetric (10.175a):
L 2T 1 pL
≤ x ≤ L: ζ2 ( x ) = L − x − sinh p ( L − x ) sech , (10.182a, b)
2 P p 2
confirming that it vanishes at both ends (10.154a, c); (ii) the slope (10.183b):
P pL
ζ′2 ( x ) = − ζ′2 ( L − x ) =
2T 1 − cosh ( p x ) sech 2 , (10.183a, b)
pL
M2 ( x ) = − E I ζ′′2 ( x ) = − E I ζ′′2 ( L − x ) = − T sinh ( px ) sech , (10.184a–c)
2
Examples 10.1 to 10.20 43
which vanishes at the pinned ends; (iv) the maximum deflection (slope in modulus) is:
L PL 2 pL
ζ2 max = ζ2 = 1− tanh , (10.185)
2 4T pL 2
P pL
ζ′2 max = ζ′2 ( 0 ) = 1 − sech = − ζ′2 ( L ) , (10.186a, b)
2T 2
at the middle (10.185) [at the supports (10.186a, b)]. In the absence of tension (10.167a)
the shape of the elastica of the bar is (10.181a–c) and the lowest-order correction for
a beam under weak traction (10.179b) is (10.179c).
E I ζ′′′′
3 − T ζ ′′
3 = ρ g ; (10.187)
d ρ g −2 ρ g x2
( )
−1
D≡ : ζ3* ( x ) = E I D4 − T D2 ρg = − D ( 1) = − . (10.188a–d)
dx T 2T
ρ g x2
ζ3 ( x ) = ζ ( x ) + ζ3* ( x ) = A + B x − + C cosh ( px ) + D sinh ( px ) .
2T
(10.189a, b)
In the case (c) of a beam (Figure 10.5c) clamped (10.155a, b) ≡ (10.190a, b) at one
end substitution of (10.189b) yields (10.161a, b), so that only two constants of
integration remain in (10.190c):
ρ g x2
ζ3 ( 0 ) = 0 = ζ3′ ( 0 ) : { } {
ζ3 ( x ) = A 1 − cosh ( px ) + D sinh ( px ) − px − } 2T
.
(10.190a–c)
44 Differential Equations and Applications
1 − cosh ( pL )
sinh ( pL ) − pL A ρ g L2 /2 T
ζ3 ( L ) = 0 = ζ3′′( L ) :
= ;
sinh ( pL ) D ρ g/p 2 T
− cosh ( pL )
(10.191a–c)
p L > 0: { }
∆ ≡ sinh ( pL ) − pL cosh ( pL ) = cosh ( pL ) tanh ( pL ) − pL ≠ 0,
(10.192a–c)
A sinh ( p L ) pL − sinh ( pL )
ρ g L2 /2 T ∆
= ; (10.193)
D cosh ( p L ) 1 − cosh ( pL ) ρ g/ p 2 T ∆
ρg ρ g x4
T = 0: ζ′′′′
03 = , ζ03 ( 0 ) = + A + B x + C x 2 + D x 3 , (10.194a–c)
EI 24 E I
whose complete integral (10.194c) consists of the sum of: (i) a particular
integral of the unforced equation (10.194b), which is the first term; (ii) the
general integral (10.167c) of the unforced equation (10.167b). The boundary
conditions for clamping (10.190a, b) at one end, and pin-joint (10.191a, b) at
the other end determine the constants of integration in (10.194c), namely
(10.195a–c):
1 L C ρ g L2 1
A = 0 = B, =− , (10.195a–c)
2 6 L D 24 E I 12
Examples 10.1 to 10.20 45
ρ g L2 5ρ g L ρ g x 2 x 2 5 x L L2
C= , D=− ; ζ03 ( x ) = − + ,
16 E I 48 E I 8 E I 3 6 2
(10.196a–c)
the shape of the bar is given by (10.194c) with (10.195a, b; 10.196a, b) leading
to (10.196c) ≡ (III.4.64).
In conclusion, (problem 224) the linear bending (10.152b) of a heavy (10.155e)
uniform (6.10a, b) beam under traction (10.152a) clamped (10.155a, b) at one end
and pinned (10.155c, d) at the other end leads to: (i) the transverse displacement of
the elastica (10.190c) with (10.193) confirming that the displacement is zero at the
supports (10.190a; 10.191a); (ii) the slope (10.197):
{
ζ3′ ( x ) = p − A sinh ( px ) + D cosh ( px ) − 1 −} ρgx
T
; (10.197)
ρgEI
M3 ( x ) = − E I ζ3′′( x ) = T A cosh ( px ) − D sinh ( px ) + ; (10.198a, b)
T
leading to the shear stress (10.199c), in agreement with (10.199d) ≡ (10.155e). The
bending momenta at the clamped support is (10.200a):
ρgEI
;M3 ( 0 ) = T A +
F3 ( 0 ) = pT D, F3 ( L ) = F3 ( 0 ) + ρ g L ,
T
(10.200a–c)
and the transverse force at the two supports (10.200b, c) differs by the weight of the
beam.
E I ζ′′′′
4 − T ζ ′′
4 = q x , (10.201)
46 Differential Equations and Applications
d q −2 q x3
( )
−1
D≡ : ζ 4* ( x ) = E I D4 − T D2 qx = − D ( x) = − . (10.202a–d)
dx T 6T
Adding (10.202d) to the general integral (10.158) of the unforced case (10.157c)
leads (10.202e) to (10.202f):
q x3
ζ4 ( x ) = ζ ( x ) + ζ 4* ( x ) = A + B x − + C cosh ( px ) + D sinh ( px ) , (10.202e, f)
6T
q x3
ζ ( 0 ) = 0 = ζ′ ( 0 ) : { } {
ζ4 ( x ) = A 1 − cosh ( p x ) + D sinh ( p x ) − p x − } 6T
.
(10.203a–c)
q q p 2 L2
ζ4′′′ ( L ) − p 2 ζ4′ ( L ) = 0 : 0 = − + p3D + , (10.204a, b)
T 2T
q L2
ζ4′ ( L ) = 0: A sinh ( pL ) − D cosh ( pL ) − 1 = − , (10.205a, b)
2Tp
q q L2 q L2
p3D = − , A = D tanh ( p L ) − csch ( p L ) − csch ( p L ) ;
T 2E I 2Tp
(10.206a, b)
In the case of a bar in the absence of traction (10.207a), the equation (10.201)
is replaced by (10.207b):
qx q x5
T = 0: ζ′′′′
04 = , ζ04 ( x ) = + A + B x + C x 2 + D x 3 , (10.207a–c)
EI 120 E I
whose complete integral is (10.207c) the sum of (10.158) with the particular
integral of (10.207b). The constants of integration are determined by the
boundary conditions at the: (i) clamped end (10.156a, b) leading to (10.195a, b);
(ii) at the free end (10.156c, d) leading to (10.208a, b):
q L2 q L3 q L3
04 ( L ) = 6 D +
0 = ζ′′′ , ′′ ( L ) = 2 C + 6 DL +
0 = ζ04 = 2C − ;
2EI 6E I 3E I
(10.208a–c)
substitution of (10.208a) in (10.208b) specifies the remaining constant of inte-
gration (10.208c); substitution of (10.195a, b; 10.208a, c) in (10.207c) yields:
q x2
ζ04 ( x ) =
120 E I
{ }
x 3 − 10 x L2 + 20 L3 , (10.209)
q x2
ζ4′ ( x ) = p D cosh ( px ) − 1 − A sinh ( px ) − , (10.210)
2T
confirming that it vanishes at the clamped end (10.156b); (iii) the curvature (10.211b),
which specifies the bending moment (10.211a):
q xEI
M4 ( x ) = − E I ζ4′′( x ) = T A cosh ( p x ) − D sinh ( p x ) + ; (10.211a, b)
T
q E I q x2
+ F4 ( x ) = E I ζ′′′ ( x ) − T ζ′ ( x ) = pT D −
, f 4 ( x ) = F4′ ( x ) = q x ,
T 2
(10.212a–d)
48 Differential Equations and Applications
qE I
M4 ( 0 ) = T A, F4 ( 0 ) = pT D − . (10.213a, b)
T
Eh3 ρ g r4
D= : ζ ( r ) = + C1 + C2 log r + ( C3 − C4 ) r 2 + C4 r 2 log r ,
(
12 1 − σ 2 ) 64D
(10.214a, b)
dζ ρ g r 3 C2
ζ′ ( r ) ≡ = + + 2 C3 r − C4 r + 2 C4 r log r , (10.215a, b)
dr 16D r
Examples 10.1 to 10.20 49
d 2ζ 3 ρ g r 2 C2
ζ′′ ( r ) ≡ = − 2 + 2 C3 + C4 + 2 C4 log r , (10.216a, b)
dr 2 16 D r
σ ρ g r2
Mr ( r ) = − D ζ′′ + ζ′ = − (3 + σ )
r 16
(10.217a, b)
C2
+ D ( 1 − σ ) 2 − C4 − D ( 1 + σ ) C3 + C4 log r ,
r
D ′ ρ g r 4 C4 D
N r ( r ) = − ( r ζ′ )′ = − − , (10.218a, b)
r 2 r
the stress couple (10.217a, b), and the turning moment (10.218a, b).
The linear bending of a heavy circular plate with a concentric (10.219a) circular
hole (Figure 10.6) leads (Table 10.2) to (problem 272) eight cases (Table 10.1):
b ≤ r ≤ a;
(10.219a–e)
I: ζ ( a ) = 0 = ζ′ ( a ) , ζ ( b ) = − d , ζ′ ( b ) = 0 ,
IV: ζ ( a ) = 0 = Mr ( a ) , Mr ( b ) = 0 = N r ( b ) , (10.222a–d)
V: ζ ( a ) = 0 = Mr ( a ) , ζ ( b ) = − d , Mr ( b ) = 0, (10.223a–d)
TABLE 10.1
Combinations of Boundary Conditions for a Circular Plate with a Circular Hole
Boundary Outer
Support Clamped Supported Free
Inner Clamped case I, Figure 10.6a VI ≡ reverse case II VII ≡ reverse case III
Supported case II, Figure 10.6b case V, Figure 10.6e VIII ≡ reverse case IV
Free case III, Figure 10.6c case IV, Figure 10.6d statically unstable
Note: There are 3 × 3 = 9 combinations of clamped/pinned/free boundaries for the weak bend-
ing of a circular elastic plate with a circular hole; of these eight are static cases, consisting
of five distinct cases (Table 10.2 and Figures 10.5a–c) and three interchanges.
50 Differential Equations and Applications
e
b
c
a
FIGURE 10.6
The linear bending of a circular elastic plate with a concerntric circular hole (Table 10.1) leads
to five distinct cases (Table 10.2), namely: (a) the case I: clamped at both boundaries; (b) the
case II: clamped at the outer boundary and supported at the inner boundary (and vice-versa
for the case VI); (c) the case III: clamped at the outer boundary and free at the inner boundary
(and vice-versa for the case VII); (d) the case IV: supported at the outer boundary and free at the
inner boundary (and vice-versa for the case VIII): (e) the case V supported at both the outer and
inner boundaries. The case IX free at both boundaries is not static and is omitted.
Examples 10.1 to 10.20 51
TABLE 10.2
Constants in the Shape of Neutral Surface of Heavy Circular Plate with a Circular Hole
Case I II III IV V
Figure 10.6a 10.6b 10.6c 10.6d 10.6e
Outer Clamped Clamped Clamped Supported Supported
Boundary
ζ ( a) = 0 ζ ( a) = 0 ζ ( a) = 0 ζ ( a) = 0 ζ ( a) = 0
ζ′ ( a ) = 0 ζ′ ( a ) = 0 ζ′ ( a ) = 0 Mr ( a ) = 0 Mr ( a ) = 0
Inner Clamped Supported Free Free Supported
Boundary
ζ (b ) = − d ζ (b ) = − d Nr (b ) = 0 Nr (b ) = 0 ζ ( a) = − d
ζ′ ( b ) = 0 Mr ( b ) = 0 Mr ( b ) = 0 Mr ( b ) = 0 Mr ( b ) = 0
Note: For the distinct cases (Figures 10.5a–e) of a circular heavy elastic plate with a concentric
circular hole (Table 10.1) the boundary conditions are indicated.
namely: (I) clamped at both boundaries (Figure 10.6a) with (10.219b–e) a height
difference d; (II) clamped at the outer boundary (10.220a, b) and supported at the
inner boundary (10.220c, d) (Figure 10.6b): (III) clamped at the outer boundary
(10.221a, b) and free at the inner boundary (10.221c, d) (Figure 10.6c); (IV) sup-
ported at the outer boundary (10.222a, b) and free at the inner boundary (10.222c,
d) (Figure 10.6d); (V) supported both at the outer boundary (10.223a, b) and the
inner boundary (10.223c, d) (Figure 10.6e); (VI, VII, VIII) the cases II, III, IV of
unequal boundaries can be reversed, respectively, for (VI) supported outside and
clamped inside, (VII) clamped inside and free outside and (VIII) supported inside
and free outside. The remaining combination IX of 3 × 3 = 9 is free both at
the inner boundary and the outer boundary, which is not statically stable
because there is no support.
∂2 ζ ∂2 ζ ∂2 ζ
Txy = 0: Txx + Tyy ( x ) = ρ , (10.224a, b)
∂x2 ∂y 2 ∂t 2
52 Differential Equations and Applications
T xy = 0 T yy = const > 0
T xx = const > 0 T xx
0
x
T yy
a
+ +
+
−
−
+
−
+ −
n=2 n=3 n=4
b c d
ζ
e
n<n
cut-off evanescent
x
cut-on
n>n
propagating
FIGURE 10.7
A homogeneous elastic membrane held at two parallel boundaries and (a) like an infi-
nite strip in the orthogonal direction, subject to constant normal tractions and no shear
stress, has transverse modes of oscillation of all orders, for example, (b) two, (c) three, and
(d) four. The modes below (above) a transition value n are (e) cut-on (cut-off) or propagating
(evanescent).
Examples 10.1 to 10.20 53
dT yy / dx > 0
T xy = 0
T xx = const > 0 T xx
0
x
a
T yy
ζ T yy
x
x
+
−
c
+
−
FIGURE 10.8
In Figure 10.6ª, the stress normal to the boundary (a) increases monotonically with dis-
tance along the boundary, each transverse mode has a turning point, such that before
(after) (b) the modes are oscillating (decaying), corresponding to (c) propagating
(evanescent) waves.
54 Differential Equations and Applications
where: (i) the mass density per unit area ρ is constant because the membrane
is homogeneous; (ii) there is no shear stress (10.224a); (iii)(iv) the stress paral-
lel (normal) to the sides is constant (non-uniform or constant). The case of
constant tension is considered first (E10.12.1) leading to propagating (evanes-
cent) that is cut-on (cut-off) modes. The case of non-uniform motion leads to
turning points (E10.12.2), as in the ray theory (notes 5.10–5.19).
nπ y
ζ ( x , y = 0, t ) = 0 = ζ ( x , y = a, t ) : ζn ( x , y , t ) = e iωt sin X ( x ) , (10.225a–c)
a
2 1/2
2 2 2 2 π n
d X n /dx + k X n = 0, Txx kn = ρω − Tyy , (10.226a, b)
n
a
ωa ρ
n< ≡ n: X n± ( x ) = C± exp ( ± i kn x ) , (10.227a, b)
π Tyy
nπ y
n < n: ζn± ( x , y , t ) = C± sin exp i ( ω t ± kn x ) , (10.228a, b)
a
nπ y
= C± sin exp i kn ( un t ± x ) , (10.228c)
a
Examples 10.1 to 10.20 55
2 −1/2
Txx ω Tyy π n
u0 = : un = = u0 1 − , (10.229a–c)
ρ kn ρ ω a
where (10.229a) is the speed (10.229b) of the zero-order mode. The modes of order
higher (10.230a) than (10.227a) are (Figure 10.7e) cut-off or evanescent modes
(10.230b; 10.225c) ≡ (10.230c):
nπy
n > n: X n± ( x ) = C± exp ( ± kn x ) , ζn± ( x , y , t ) = C± sin exp ( ± kn x ) cos ( ωt ) ,
a
(10.230a–c)
2 x > x : propagation,
aω
kn ( x ) = 0: Tyy ( x ) = ρ : (10.231a–d)
nπ x < x: evanescence,
and the traction takes the value (10.131b). Before the turning point (10.231c),
there are vibrations like (10.228a–c; 10.229a, b) cut-on or propagating modes
(Figure 10.8b, c); beyond the turning point, the modes decay like cut-off or evanes-
cent modes (10.230a–c). In a strict sense, the wavenumber k only exists (10.226b)
when it is constant, that is for uniform traction (Figure 10.7a–e), corresponding to
sinusoidal vibrations; in the case of non-uniform tension the wavenumber (10.226b)
is not constant:
1/2
ρω 2 Tyy ( x ) nπ
2
kn ( x ) = − , (10.232)
Txx Txx a
56 Differential Equations and Applications
nπy x
ζn ( x , y , t ) = C± kn ( x )
∫ kn ( ξ ) dξ ,
−1/2
x > x: sin exp i ω t ±
a
(10.233a, b)
nπy x
x < x : ζn ( x , y , t ) = C± kn ( x )
∫ kn ( ξ ) dξ ,
−1/2
sin cos ( ωt ) exp ±
a
(10.234a, b)
when (7.431a–g) the local wavelength is much smaller than the lengthscale of varia-
tion of the stress:
2
2π
λ n ( x ) 2 =
k ( x )
2 2
(
≤ L = Tyy ∂Tyy /dx ) −2
. (10.235)
n
If the ray approximation (notes 5.16–5.17, 7.35–7.27) does not hold, then the
notion of wavenumber or wavelength loses its meaning, even locally, and exact
global solutions must be sought to describe accurately vibrations or waves.
The vibrations of a membrane held are similar to quasi-one-dimensional
waves in non-uniform medium (notes 7.28–7.45), so analogous methods
apply.
{X , Y , Z} = {y 2 − z2 , z2 − x 2 , x 2 − y 2 } ;
(10.236a, b)
x 2 ( dx ) − y 2 ( dy ) + ( dz ) − 2 x dx dz = 0.
2 2 2
dx dy dz
X + Y + Z = 0 = x 2 X + y 2Y + z 2 Z ; = = ; (10.237a–c)
X Y Z
Examples 10.1 to 10.20 57
the curve tangent to the vector field (10.236a) satisfies (7.12a–c) ≡ (10.237c)
leading to (10.238a, b):
dx + dy + dz = 0 = x 2 dx + y 2 dy + z 2 dz: x + y + z = C1 , x 3 + y 3 + z 3 = C2 ,
(10.238a–d)
0 = A ( dx ) + B ( dy ) + C ( dz ) + 2 D dx dy + 2 E dx dz + 2 F dy dz ; (10.239)
2 2 2
{ A , B , C , D , E , F } = { x 2 , − y 2 , 1 , 0 , − x , 0 } . (10.240a–f)
A BC + 2 DE F − A F 2 − BE 2 − C D2 = − x 2 y 2 + x 2 y 2 = 0 ; (10.240g)
( )
2
0 = ( C dz + E dx + F dy ) −
2
E 2 − C A dx + F 2 − C B dy
(10.241a–c)
= ( dz − x dx ) − ( y dy ) = ( dz − x dx − y dy )( dz − x dx + y dy ) .
2 2
x 2 + y 2 = 2 z + C1 , x 2 − y 2 = 2 z + C2 , (10.242a, b)
LXXIV: ( )
y 2 dx + 2 x y + z 2 dy + 2 y z dz = 0 ; x y 2 + y z 2 = C ; (10.243a, b)
LXXV: y z dx + 2 x z dy + 3 x y dz = 0 ; x y 2 z 3 = C ; (10.244a, b)
LXXVI: 2 y dx − x dy + z dz = 0 , y = C1 x: C1 x 2 + z 2 = C2 ; (10.245a–c)
LXXXII: ( ) ( )
y 2 dx + 2 x y + z 2 dy + 2 y z + u2 dz + 2 z ud u = 0:
(10.246a, b)
x y 2 + y z 2 + zu2 = C ;
y dx + x dy + z dz − z du = 0, y = C1 x , u = C2 z: 2 C1 x 2 − C2 z 2 = − z 2 .
(10.248a–d)
( )
0 = X .dx = y 2 dx + 2 x y dy + z 2 dy + 2 y z dz = d y 2 x + z 2 y , (10.249d)
( )
λ = y z 2 : 0 = λ X . dx = y 2 z 3 dx + 2 x y z 3 dy + 3 x y 2 z 2 dz = d x y 2 z 3 , (10.250e, f)
the integrating factor (10.250e) leads to the solution (10.250f) ≡ (10.244b). Thus,
the family of surfaces (10.244b) is orthogonal to the vector field (10.250b). The differ-
ential (10.245a) ≡ (10.251a, b) is (7.29 a–d) inexact (10.251c) and non-integrable
(10.251d):
0 = X .dx , X = { 2 y , − x , z} ; ∇ ∧ X = {0, 0 − 3} , X . ∇ ∧ X = −3 z ; ( )
(10.251a–d)
1
y = C1 x: 0 = 2 y dx − x dy + z dz = C1 x dx + z dz =
2
( )
d C1 x 2 + z 2 , (10.251e, f)
( ) ( ) ( ) (
0 = y 2 dx + 2 xy dy + z 2 dy + 2 yz dz + u2 dz + 2 zudu = d y 2 x + z 2 y + u2 z . )
(10.252c, d)
λ = yz 2 u3 : λ ( y z u dx + 2 x z u dy + 3 x y udz + 4 x y z du )
(
= y 2 z 3 u4 dx + 2 x y z 3 u4 dy + 3 x y 2 z 2 u4 dz + 4 x y 2 z 3 u3 du = d x y 2 z 3 u4 , )
(10.253c–e)
60 Differential Equations and Applications
0 = X . dx : X = { y , x , z , − z } , (10.254a, b)
0 = C1 x dx + x d ( C1 x ) + z dz − z d ( C2 z ) = 2C1 x dx + ( 1 − C2 ) z dz
z2 (10.254c–e)
= d C1 x 2 + ( 1 − C2 ) ,
2
TABLE 10.3
Boundary Conditions for the Bending of a Plate
Arbitrary Regular
Shape of Boundary Curve Circular Rectangular
Figure 6.24a 6.24b 6.24c
Equation (6.467e) (6.467c) (6.467a)
Slope (6.464a) (6.464c) (6.464e)
Stress Couple* (10.261)/(10.280) (10.283b)/(6.471d) (10.283d)/(6.471b)
Turning Moment* (10.262)/(10.281) (6.467d)/(6.467d) (6.467b)/(6.467b)
Clamped* (6.474a, b, d; 10.263b) (6.473a, b, d; 10.284)/ (6.472a–d)/(6.472a–d)
(10.271a–e) (6.474a–d) (6.473a–d)
Figure 6.23a
Supported or Pinned* (6.477a, b, d; 10.285)/ (6.476a, b, d; 10.276)/ (6.475a–d)/(6.475a–d)
(10.273a–e) (6.477a–d) (6.476a–d)
Figure 6.23b, c
Free Boundary* (6.480a, b; 10.278; (6.479a, b, d; 10.276)/ (6.478a, b, d; 10.282)/
(10.277a–e) 10.279)/ (6.479a–d) (6.478a–d)
Figure 6.23d (6.480a–d)
* for a plate made of a pseudo-isotropic orthotropic (isotropic) material before (after) the dash:
* first bending stiffness (6.413c, d) [(6.579a; 6.581a)];
* second bending stiffness (10.255a) [(10.255b)];
Note: The boundary conditions for the linear bending of an elastic plate lead to 3 × 3 × 2 = 18
cases combining: (i) three types of boundary: clamped, pinned, and free; (ii) three shapes:
arbitrary closed regular curve, circular, and straight; (iii) two types of elastic material:
isotropic and pseudo-isotropic orthotropic.
h3
C44 = C44 → D ( 1 − σ ) = const:
12
where the first bending stiffness (6.593a; 6.595a) [(6.413c, d)] for a pseudo-
isotropic orthotropic (isotropic) plate need not be constant, and the second
stiffness (10.255b) [(10.255a)] is constant (6.454a). The coefficient of the first
(second) term of (10.255c) is in agreement with (6.609b) [(6.609c)] the turn-
ing moment (6.430a) [the stress couple (6.422a–c; 6.423a–c)] with additional
contributions coming from the third term on the r.h.s. of (10.255c). The latter
62 Differential Equations and Applications
involves the deformation vector I given by (6.448d) whose normal compo-
nent (6.451e) appears in (10.255c). Designating by θ the angle of the normal
with the x-axis (Figure 6.24a) leads to the unit normal (6.462a) [tangent
(6.462c)] vectors and associated derivatives (6.462b) [(6.462d)]. Using (6.462a),
the normal component of the deformation vector (6.448d) is given by:
( ) (
I n = I . n = cos θ ∂ xy ζ ∂ y ( δζ ) − ∂ yy ζ ∂ x ( δζ ) )
(10.256)
( )
+ sin θ ∂ xy ζ ∂ x ( δζ ) − ( ∂ xx ζ ) ∂ y ( δζ ) ,
for substitution in the last term of the r.h.s. of (10.255c), leading to an integra-
tion by parts along the boundary (E10.15.2).
∫∂c
C44 I n ds =
∫ ∂c
{ ( ) (
C44 2 sin θ cos θ ∂ xy ζ − sin 2 θ ( ∂ xx ζ ) − cos 2 θ ∂ yy ζ δ ( ∂n ζ ) )
( ) ( )
+ cos θ sin θ ∂ yy ζ − ∂ xx ζ + cos 2 θ − sin 2 θ ∂ xy ζ ∂ s ( δζ ) ds , }
(10.257)
sin ( 2θ )
X≡
2
( )
∂ yy ζ − ∂ xx ζ + cos ( 2θ ) ∂ xy ζ, (10.258)
∫ ∂c
C44 X ∂ s ( δζ ) ds +
∫
∂c
(
δζ ∂ s C44 X ds = ) ∫ ∂ (C
∂c
s 44 )
X δζ ds = C44 X δζ = 0,
(10.259a–d)
and then in (10.255c) leads to (6.609a) ≡ (10.260), the elastic energy per unit
length of the boundary:
dEˆ
{ ( ) ( )
E d = d = − ∂n D ∇ 2ζ + C44 ∂ s sin θ cos θ ∂ xx ζ − ∂ yy ζ − cos ( 2θ ) ∂ xy ζ δζ,
ds }
{ }
+ D ∇ 2ζ + C44 sin ( 2θ ) ∂ xy ζ − sin 2 θ ∂ xx ζ − cos 2 θ ∂ yy ζ δ ( ∂n ζ ) ,
(10.260)
( ) ( )
N n = ∂n D ∇ 2ζ − C44 ∂ s sin θ cos θ ∂ xx ζ − ∂ yy ζ − cos ( 2θ ) ∂ xy ζ . (10.262)
In the case of the nonlinear strong bending of a plate with in-plane stresses, the
augmented turning moment is given by (6.682a–c). If the bending moment D is con-
stant: (i) it can be taken out of the derivative in the first term on the r.h.s. of (10.262);
(ii) from (10.255b) the Poisson ratio is also constant. If (problem 303) the plate is
supported (Figure 6.23b) or pinned (Figure 6.23c) the vanishing of the displacement
on the boundary (10.263a) simplifies the stress couple (10.261) [turning moment
(10.262)] to (10.263b) [(10.263c, d)]:
dθ dθ
ζ = 0: Mn = − D ∂nm ζ + ∂n ζ + C44 ∂n ζ . (10.263a, b)
ds ds
Nn dθ
D = const:
D
( )
= ∂n ∇ 2 ζ = ∂nnn ζ +
ds
∂nn ζ . (10.263c, d)
If (problem 304) the plate is clamped (Figure 6.23a) the conditions (10.264a, b)
simplify the stress couple to (10.264c) and leaves the turning moment unchanged
(10.263d) ≡ (10.264d):
Nn dθ
ζ = 0 = ∂n : Mn = − D ∂nn ζ, = ∂nnn ζ + ∂nn ζ . (10.264a–d)
D ds
64 Differential Equations and Applications
The relations (10.261, 10.262; 10.263a–d; 10.264a–d) hold for a plate with any shape
(Figure 6.24a) provided that the boundary is a closed regular curve. The proof of
(10.263a–d) is made next (E10.15.4) before considering the three main types
of boundary conditions (E10.15.6–E.10.15.8) for pseudo-isotropic orthotropic
and isotropic plates (E10.15.5).
dθ
∂ xy ζ = ( cos θ ∂n − sin θ ∂ s ) sin θ ∂n ζ = cos θ sin θ ∂nn ζ − ∂n ζ , (10.266a)
ds
dθ
∂ yy ζ = ( sin θ ∂n + cos θ ∂ s ) sin θ ∂n ζ = sin 2 θ ∂nn ζ + cos 2 θ ∂n ζ , (10.266b)
ds
dθ
∂ xx ζ = ( cos θ ∂n − sin θ ∂ s ) cos θ ∂n ζ = cos 2 θ ∂nn ζ + sin 2 θ ∂n ζ , (10.266c)
ds
where (10.267a) was used since θ is the direction of the normal, hence, does
not vary along it, implying (10.267b), which is the inverse (10.267c) of the
radius of curvature (6.461c):
dθ 1 dθ
∂n θ = 0,
= ; ζ = 0: ∇ 2ζ = ∂ xx ζ + ∂ yy ζ = ∂ nn ζ +
∂s θ ≡ ∂n ζ ,
ds R ds
(10.267a–f)
from (10.266b, c) follows that the Laplacian (10.267e) along the boundary
(10.267d) is given by (10.267f) in terms of normal derivatives in agreement
with the term in curved brackets in (6.474d). These results apply to pseudo-
isotropic orthotropic plates, which include as the particular case the isotro-
pic plate (E10.15.5).
Examples 10.1 to 10.20 65
dθ
( )
N n = ∂n D ∇ 2ζ = ∂n D ∂nn ζ +
ds
∂ n ζ ;
(10.268a, b)
(
= − 2 sin 2 θ cos 2 θ + sin 4 θ + cos 4 θ ) ddsθ ∂ ζ (10.269)
n
dθ dθ
( )
2
= − sin 2 θ + cos 2 θ ∂n ζ = − ∂n ζ .
ds ds
Mn dθ dθ
− = ∇ 2ζ − ( 1 − σ ) ∂n ζ = ∂nn ζ + σ ∂n ζ , (10.270a, b)
D ds ds
where the Lapacian (10.267f) was used in (10.270b). The boundary conditions are
considered next for the three main cases of (i) clamped plate (E10.15.6), (ii)
supported or pinned plate (E10.15.7), and (iii) free boundary (E10.15.8).
s ( θ ) ∈C1 ( ∂C ) : ζ = 0 = ∂n ζ , Mn ≠ 0 ≠ N n , (10.271a–e)
66 Differential Equations and Applications
1 1
(1
Mr = − D ∂rr ζ + ∂r ζ − C44 ∂r ζ = − D ∂rr ζ − D + C44 ∂r ζ .
r r r
) (10.272a, b)
s ( θ ) ∈C 1 ( ∂C ) : ζ = 0 ≠ ∂ n ζ, Mn = 0 ≠ N n . (10.273a–e)
C44 dθ dθ
χ≡ : ∂nn ζ + ∂n ζ = χ ∂n ζ ; (10.274a, b)
D ds ds
dθ
χ = 1 − σ: ∂nn ζ + σ ∂n ζ = 0. (10.275a, b)
ds
In the case of a supported or pinned circular plate, made of an isotropic mate-
rial, the boundary conditions are (6.476a–d), where only (6.476c) is changed
for a pseudo-isotropic orthotropic material to (10.276):
1 χ
∂rr ζ + ∂r ζ = ∂r ζ . (10.276)
r r
s ( θ ) ∈C 1 ( ∂C ) : ζ ≠ 0 ≠ ∂n ζ: Mn = 0 = N n . (10.277a–e)
The vanishing of the stress couple (10.261) [turning moment (10.262)] at the
free boundary lead to the first (10.278) [second (10.279)] boundary conditions
for a pseudo-isotropic orthotropic plate involving the parameter (6.606a–c):
( ) ( )
χ −1 ∂n ∇ 2ζ = ∂ s sin θ cos θ ∂ xx ζ − ∂ yy ζ − cos ( 2θ ) ∂ xy ζ . (10.279)
In the case (10.275a) of an isotropic plate, the general stress couple (10.261)
[turning moment (10.262)] simplifies to (10.280) [(10.281)]:
{ }
Mn = − D ∇ 2ζ − ( 1 − σ ) sin ( 2θ ) ∂ yy ζ − sin 2 θ ∂ yy ζ − cos 2 θ ∂ yy ζ , (10.280)
( ) ( )
N n = ∂n D ∇ 2ζ − D ( 1 − σ ) ∂ s sin θ cos θ ∂ xx ζ − ∂ yy ζ − cos ( 2θ ) ∂ yy ζ . (10.281)
∂ yy ζ + ( 1 − χ ) ∂ xx ζ = 0. (10.282)
material, in particular isotropic; (iii) the three main boundary conditions are clamped
(Figure 6.23a), supported (Figure 6.23b), or pinned (Figure 6.23c) or free bound-
ary (Figure 6.23d). The boundary conditions (6.480a–d) involve the displacement,
the slope, and the normal components of the stress couple and the turning moment.
Normal component means radial for a circular plate and along the y-axis for the
straight side along the x-axis of a rectangular plate. For arbitrary/circular/rectangu-
lar plate the boundary is given, respectively, by (6.467e, c, a) and the slope by (6.464a,
c, e). The normal component of the stress couple (turning moment) is given by
(10.261) [(10.262)] for a pseudo-isotropic orthotropic plate and by (10.280) [(10.281)]
for an isotropic plate. The stress couple for a pseudo-isotropic orthotropic (isotropic)
plate simplifies to (10.283b) [(6.471d)] for a circular plate (10.283a) and to (10.283d)
[(6.471b)] for a straight boundary (10.273c):
Mr 1− χ My
r = const:
= ∂rr ζ + − ; y = 0: − = ∂ yy ζ + ( 1 − χ ) ∂ xx ζ .
D r D
(10.283a–d)
The turning moment for a circular (6.467d) [rectangular (6.467b)] plate is the same
for the isotropic and the pseudo-isotropic orthotropic plate.
The general form of the stress general couple (10.261) [turning moment (10.262)]
simplifies: (i) to (10.263b)-[(10.263d)] for supported or pinned boundary (10.263a, c);
(ii) to (10.264c) [(10.264d)] for clamped boundary (10.264a, b). Thus, the boundary con-
ditions for arbitrary/circular/straight boundaries, respectively, of a pseudo-isotropic
orthotropic (isotropic) plate are: (i) for clamped (6.474a, b, d; 10.263b)/(6.473a, b, d;
10.284)/(6.472a–d) [(6.474a–d)/(6.473a–d)/(6.472a–d)]:
1 1
Mr = − D ∂rr ζ + ∂r ζ − C44 ∂r ζ ; (10.284)
r r
1− χ
∂rr ζ + ∂r ζ = 0 ; (10.285)
r
f D 2
( ) ( ) ( )(
= ∇ ζ − ∂ yy Θ ( ∂ xx ζ ) − ( ∂ xx Θ ) ∂ yy ζ + 2 ∂ xy Θ ∂ xy ζ , (10.286)
h h
)
replacing the first bending stiffness for an isotropic plate (6.413d) ≡ (10.287a)
by that for a pseudo-isotropic orthotropic plate (6.593a; 6.595a) ≡ (10.287b):
E h3 h3 C13
← D → C11 − C13 C . (10.287a, b)
(
12 1 − σ 2
)
12 23
Concerning the second complementary equation (6.657) of the pair: (i) the
derivation from (6.659a, b) to the first equality in (6.666a–c) is independent
of the type of material; (ii) the elimination of the in-plane displacements is
made (10.288), as in (6.667a), retaining the strains from the first equality in
(6.666a–c) leading to:
1
2
2
( )
∂ yy ( ∂ x ζ ) + ∂ xx ( ∂ x ζ ) − 2 ∂ xy ( ∂ x ζ ) ∂ y ζ = ∂ yy Sxx + ∂ xx Syy + 2 ∂ xy Sxy ;
2
(10.288)
( ∂ ζ) − ( ∂ ζ)( ∂ ζ) = ∂ ( A
xy
2
xx yy yy 11 ) (
Txx + A12 Tyy + ∂ xx A12 Txx + A22 Tyy )
(10.289)
(
+ 2 ∂ xy A44 Txy ,)
and on the r.h.s. the compliance matrix (6.583) is used for an orthotropic
material, bearing in mind that there are only in-plane stresses (6.312c–e);
(iv) substituting (6.655e) in (10.289) it follows that (problem 318) the comple-
mentary equation relating the transverse displacement to the stress function for the
strong bending of an orthotropic plate (10.290b):
assuming that the components of the compliance matrix are constant (10.290a).
The first (10.286) [second (10.290a, b)] equation of the coupled pair is valid
for a pseudo-isotropic (general) orthotropic plate. In the complementary
70 Differential Equations and Applications
σ yx 1
{ A22 , A11 , 2 ( A12 + A44 ) } = E1 ,
1
Ex
,− 2 +
Ex Gyz
. (10.291a–c)
y
are needed, namely: (i)(ii) three (one) in the first (10.287b) ≡ (10.293b) [second
(6.595b) ≡ (10.293a)] bending stiffness:
C44 =
h3Gyz
, D=
h3
1 − σ yz σ zy −
(
Ex σ zx + σ zy σ yx
2
)
; (10.293a, b)
6 12 C0 Ey Ez Ey σ zy + σ zx σ xy
(iii)(iv) five (six) in the first (6.596a) ≡ (10.294a) [second (6.596b) ≡ (10.294b)]
conditions for a pseudo-isotropic orthotropic material:
( ) (
Ez 1 − σ xy σ yx Ey ( 1 − σ xz σ zx ) − Ex 1 − σ yz σ zy )
(10.294a)
( ) ( )
2 2
= Ey σ zy + σ zx σ xy − Ex σ zx + σ zy σ yx ,
E
( ) (
Ez 1 − σ xy σ yx 1 − σ xz σ zx − x σ yx + σ yz σ zx
Ey
)
(10.294b)
( ) ( ) (
= σ xy + σ xz σ zy Ey σ zy + σ zx σ xy − Ex σ yz + σ yz σ zx . )
In conclusion, (problem 319) the transverse displacement and stress function in
the strong non-linear deflection of an elastic plate, made of a pseudo-isotropic ortho-
tropic material, satisfy the coupled non-linear fourth-order partial differential equa-
tions (10.286; 10.290a, b) where: (i) the first bending stiffness (10.287b) is given
by (10.293b); (ii) the second bending stiffness that appears in the boundary condi-
tions (10.261; 10.262) is given by (10.293a). Four (seven) terms of the compliance
(6.583) [stiffness (6.586)] matrix appear (10.291a–c) [(10.292a–g)] in the balance
Examples 10.1 to 10.20 71
equations (i-ii) and the boundary conditions (ii-iii). The three Young moduli Ei , six
Poisson ratio σ ij, and three shear moduli Gij satisfy (6.584b) [(6.587a–c)] in the com-
pliance (stiffness) matrix for an orthotropic material. In addition, (problem 320) a
pseudo-isotropic orthotropic material satisfies two additional relations (6.596a, b) ≡
(10.294a, b).
y ′′z + y z ′ = 0, z ′′ + z y = 0, (10.295a, b)
dy dz yv
= y′ ≡ u , = z ′ ≡ v , u′ = y ′′ = − , v ′ = z ′′ = − z y . (10.296a–d)
dx dx z
Thus, the order of the system cannot exceed four, but could be less.
Eliminating for y leads to:
yv ′
2
y ′y v y 2 v ′ y 2 v z ′ yv
y ′′′y = − y = − − + = y ′y ′′ + y 3 + = y ′y ′′ + y 3 + y ′′ 2 ,
z z z z2 z
(10.297a–d)
which is of the third-order. Thus, the system could be of order three or order
four. To find out which is the case, the system is eliminated for z:
z ′z ′′
z ′′′ = − ( z y )′ = − z ′y − y ′z = − y ′z, (10.298a–c)
z
72 Differential Equations and Applications
z ′z ′′ ′ z ′′ ′
z 2 z ′′′′ = z 2 − z 2 z ′y ′ − z 3 y ′′ = z ′′′ z ′z + z ′′ 2 z − z ′′z ′ 2 + z 2 z ′ + z 2 z ′y ,
z z
y ′′ + c z ′′ + α y + β z = 0 = y ′ + z ′ + γ y + δ z , (10.300a, b)
P2 ( D ) = λ ( D − µ ) ,
2
c = 1: (10.302a, b)
D2 + α c D2 + β y ( x)
= 0 , (10.303)
D + γ D + δ z ( x)
( ) (
P3 ( D ) = D2 + α ( D + δ ) − c D2 + β ( D + γ ) )
(10.304a, b)
= ( 1 − c ) D3 + ( δ − c γ ) D2 + ( α − β ) D + α δ − β γ ,
c = 0: P3 ( D ) = D3 + δ D2 + ( α − β ) D + α δ − β γ , (10.305a, b)
c = 1: P2 ( D ) = ( δ − γ ) D2 + ( α − β ) D + α δ − β γ , (10.306a, b)
c = 0: ( )
P3 ( D ) = D3 − 1 ( D + 2 ) = D3 + 2 D2 − D − 2 , (10.307a–c)
δ = 2: β − α = 1, β γ − α δ = 2: β = 1 + α, γ = 2, (10.308a–e)
y ′′ + α y + ( 1 + α ) z = 0 = y ′ + z ′ + 2 ( y + z ) . (10.309a, b)
The roots (10.310a) of the characteristic polynomial (10.301b) lead to the gen-
eral integral (10.310b):
c = 1: P2 ( D ) = λ D2 − 2 λ µ D + µ 2 λ , (10.311a, b)
δ − γ = λ, β − α = 2 λ µ, α δ − β γ = µ 2 λ , (10.312a–c)
α δ − β γ = µ2 λ =
( 2 µ λ )2 = (β − α )2 , (10.313a–d)
4λ 4 (δ − γ )
β−α
( β − α )2 = 4 ( δ − γ ) ( α δ − β γ ) : λ = δ− γ, 2µ = , (10.314a–c)
δ−γ
74 Differential Equations and Applications
in which case the parameters ( λ , µ ) are given by (10.312a, b) ≡ (10.314b, c). The
differential system (10.300a, b; 10.302a) ≡ (10.315a, b):
y ′′ + z ′′ + α y + β z = 0 = y ′ + z ′ + γ y + δ z, y ( x ) = e µx ( D1 + D2 x ) .
(10.315a–c)
has the general integral (10.315c) where D1− 2 are arbitrary constants.
In the case (10.309a, b) substitution of (10.310b) in (10.309a) ≡ (10.316a) leads
to the second general integral (10.316b):
( 1 + α ) z ( x ) = − y ′′ ( x ) − α y ( x ) = − ( α + 4 ) C1 e −2 x − ( 1 + α ) ( C2 cosh x + C3 sinh x ) .
(10.316a, b)
z ( x ) = e µx ( E1 + E2 x ) : D2 ( µ + γ ) + E2 ( µ + δ ) = 0,
(10.317a–c)
D1 ( µ + γ ) + D2 + E1 ( µ + δ ) + E2 = 0.
µ+γ µ+γ γ −δ
E2 = − D2 , E1 = − D1 + D2 , (10.318a, b)
µ+δ µ+δ ( + δ )2
µ
µ+δ γ −δ
z ( x ) = − e µx D1 + D2 x − D2 . (10.319)
µ+γ (µ + γ ) (µ + δ )
In conclusion, the coupled system of linear differential equations with constant coef-
ficients (10.300a, b) for (10.301a) [(10.302a)]: (i) simplifies to (10.309a, b) [(10.315a,
b)]; (ii) has the cubic (10.301b) [quadratic (10.302b)] characteristic polynomial if
the four constants satisfy (10.308a, d, e)[(10.314a–c)]; (iii) has the general integral
(10.310b; 10.316b) [(10.315c; 10.319)] where C1,2,3 ( D1,2 ) are arbitrary constants.
[E10.18.1 (E10.18.2)]; (iii) linear finite difference equations with constant coef-
ficients (E10.18.3–E10.18.6).
y ′ + z′ − y = e 2 x + x, y ′′ + z ′ − y + z = e x + x 2 , (10.320a, b)
D − 1 D y ( x)
D≡
d
: = e 2 x 1 + e x 0 + x , (10.321a, b)
dx D − 1 D + 1
2
z ( x) 0
2
1 x
P3 ( D ) =
D − 1 D ( )
= D2 − 1 ( 1 − D ) = − ( D − 1) ( D + 1)
2
D2 − 1 D + 1 = − D3 + D2 + D − 1,
(10.322a–d)
with double (single) root 1(−1) leading to the general integral (10.323a) for y(x),
where (A, B, C) are arbitrary constants:
y ( x) = ( A + B x) ex + C e− x ; z′ ( x ) = y ( x ) − y ′ ( x ) , (10.323a, b)
the general integral (10.323b) for z follows from (10.320a) without forcing and
substitution of (10.323a) leads to (10.323c):
z′ ( x ) = − B e x + 2 C e − x , z ( x) = − B ex − 2 C e− x , (10.323c, d)
y ( x) bx
b + 1 −b 1
= lim e
z ( x ) b→ 2 P3 ( b ) 1 − b 2 b − 1 0
(10.324a–d)
e 2x b+1 e 2x
3 −1
= lim =− =e
2x
.
P3 ( 2 ) b→ 2 1 − b 2 3 −3 1
76 Differential Equations and Applications
The forcing by the second term on the r.h.s. of (10.321b), which is a doubly
resonant exponential (7.122a, b) is given by (10.325a–d):
y ( x) b + 1 −b 0
= lim 1 ∂ 2
2
bx
e
z ( x ) b→1 P3′′( b ) ∂b
1 − b b − 1
2
1
1 ∂2 bx −b
= lim e
b→ 1 2 − 6 b ∂b 2 b − 1
(10.325a–d)
1 −b −1
= − lim e bx x 2 + 2x
4 b→1 b − 1 1
e x x( x + 2)
= .
4 −2 x
The forcing by the last term on the r.h.s. of (10.321b) leads to particular inte-
grals (10.326a–d) using the inverse matrix of polynomials of derivatives:
y ( x) D + 1 − D x 1− x
= 1 = −
1
z ( x ) P3 ( D ) 1 − D D − 1 x 2
2
1 − D − D2 + O D3( ) 3 x − x 2
1 − x x
= − 1 + D + 2 D2 + O D3 ( )
= 2
3 x − x 2 x − x + 1
.
(10.326a–d)
Adding the general integral (10.323a, d) of the unforced system on the l.h.s. of
(10.320a, b) ≡ (10.321a, b) to the particular integrals (10.324d; 10.325d; 10.326d)
corresponding to the three forcing terms on the r.h.s. leads to:
x( x + 2) x
y ( x) = ( A + B x) ex + C e− x − e2x + e + x , (10.327a)
4
x ex
z ( x) = − B ex − 2 C e− x + e2x − + x 2 − x + 1, (10.327b)
2
as the complete integral (10.327a, b) of the coupled forced system of linear differential
equations with constant coefficients (10.320a, b) ≡ (10.321a, b).
Examples 10.1 to 10.20 77
d δ − 1 δ y (x)
δ≡x : = x 2 1 + x 0 + 1 log x ,
dx δ ( δ − 1) + δ − 1 δ + 1 z (x) 0 1 2
(10.329a, b)
C x
y ( x ) = ( A + B log x ) x + − x 2 + log x ( 2 + log x ) + log x , (10.330a)
x 4
2C x
z ( x) = − B x − + x 2 − log x + log 2 x − log x + 1, (10.330b)
x 2
as the complete integral (10.330a, b) of the forced coupled linear system of differential
equations with homogeneous coefficients (10.328a, b) ≡ (10.329a, b).
∆ − 1 ∆ yn
∆ { y n , zn } ≡ { y n + 1 , zn + 1 } : 2 = 0 , (10.332a, b)
∆ − 1 ∆ + 1 zn
e ax ↔ xa ↔ an , xn ↔ log n x ↔ n , (10.333a–d)
lead from the first three (two) terms of (10.327a) ≡ (10.330a) [(10.327b) ≡
(10.330b)] to (10.334a) [(10.334b)]:
yn = A + B n + C ( − ) , zn = − B − 2 C ( − ) ,
n n
(10.334a, b)
as the general solution (10.334a, b) of the coupled unforced system of linear finite
difference equations with constant coefficients (10.331a, b).
∆ − 1 −1 y 1
∆ + 1 ∆
= b cos ( 2 φ ) = Re b e
z 1
i 2 φ
(
. ) (10.336a, b)
R2 ( ∆ ) = Det Rm , r ( ∆ ) = ∆ ( ∆ − 1) + ∆ + 1 = ∆ 2 + 1; (10.338a–c)
Examples 10.1 to 10.20 79
∆ 1 1 1 + ∆
Rr , m ( ∆ ) Gm = = . (10.339a, b)
−∆ − 1 ∆ − 1 1 −2
( )
R2 b e i 2 φ = b 2 e 4 iφ + 1 = b 2 cos ( 4φ ) + 1 + i b 2 sin ( 4φ ) , (10.341)
1 b 2 cos ( 4φ ) + 1 − i b 2 sin ( 4φ ) b 2 e − i 4φ + 1
b 2 e i 4φ ≠ −1:
= = ,
R2 b e i2φ
( )
b 2 cos ( 4φ ) + 1 + b 4 sin 2 ( 4φ ) 1 + b + 2 b cos ( 4φ )
2 4 2
(10.342a–c)
y
b e i 2 φ ≠ ± i : = Re
( )(
b e i 2 φ 1 + b 2 e − i 4φ ) 1 + b e i 2 φ
1 + b + 2 b cos ( 4φ )
4 2
z − 2
e i 2 φ + b e i 2 ( + 1)φ + b 2 e i 2 ( − 2 )φ + b 3 e i 2 ( − 1)φ
b
= × Re ,
1 + b 4 + 2 b 2 cos ( 4φ ) −2 e i 2 φ − 2 b 2 e i 2( − 2 )φ
(10.343a–c)
π cos ( 2 φ ) + b 2 cos 2 ( − 2 ) φ
(b, φ) ≠ 1, ± : z = − 2 b , (10.344a, b)
4 1 + b 4 + 2 b 2 cos ( 4φ )
c ≡ b e i 2 φ = ± i: ( )
y + 1 − y − z = Re c = y + 1 + z + 1 + y , (10.345a–c)
π
( ) (
Re c = Re e ± i π /2 : ) y + 1 − y − z = cos
2
= y + 1 + z + 1 + y .
(10.346a–c)
∆ − 1 −1 y 1
= c , (10.347a)
∆ + 1 ∆ z 1
has for (10.347b) the solutions (10.347c, d), where (10.337b; 10.338c) were used:
y 1 ∂ 1 + c
c ≠ ± i: = Re lim c −
c →± i R2′ ( c ) ∂ c
z
2
log c c + c + 1 π 1 ± i
= Re lim = Re c ± i π /2
c→± i 2 c −2 4 −2
Thus, the system (10.345b, c) [(10.346b, c)] has particular solution (10.347c, d) ≡
(10.348a–c) [(10.347e, c) ≡ (10.349a–c):
1+ c 2 c
c ≠ ± i: y = c , z = − ; (10.348a–c)
1 + c2 1 + c2
Examples 10.1 to 10.20 81
π π
y = A cos + B sin , (10.350a)
2 2
π π
z = ( B − A ) cos − ( A + B ) sin , (10.350b)
2 2
whose value is indicated in Table 10.4, where (A,B) are arbitrary constants deter-
mined by initial values. The complete solution adds to the general solution (10.350a,
b), the particular solutions (10.344a–c), simplifying to (10.348a–c) for (10.345a–c)
and (10.349a–c) for (10.349a, b). The general solution (10.350a, b) is obtained as
follows: (i) the characteristic polynomial (10.338c) has roots (10.351a, b) lead-
ing to the solutions (10.351c, d):
π π
∆ = ± i = e ± iπ/2 : (e )
± i π/2
= e ± i π/2 = cos ± i sin , (10.351a–d)
2 2
TABLE 10.4
General Solution (10.350a, b) of Coupled System of Finite Difference Equations
(10.346a, b)
= 4n 4n + 1 4n + 2 4n + 3
π 1 0 –1 0
cos
2
π 0 1 0 –1
sin
2
y A B –A –B
z B–A A–B A–B A+B
Note: The general integral (10.350a, b) of the coupled system of linear finite difference equa-
tions with constant coefficients (10.334a, b) leads to four cases, for , a multiple of 4, or
adding 1, 2, or 3.
82 Differential Equations and Applications
of which a real linear combination with arbitrary constants (A, B) was chosen
in the r.h.s. (10.350a); (ii) substituting (10.350a) in the l.h.s. of (10.335a) leads
to (10.351e, f):
π π π π π π
z = y + 1 − y = A cos + − cos + B sin + − sin
2 2 2 2 2 2
π π π π
= A − sin − cos + B cos − sin ,
2 2 2 2
(10.351e–g)
∂ H dx ∂ H dv ∂H µ v + f ∂H
H > 0 > H = + =v − ; (10.353a–d)
∂ x dt ∂ v dt ∂x m ∂v
2
a a2
β > 0: + b β > 0 → − 2 ab < + b 2β, (10.354a–c)
β β
mα 2 m x
m 2
2
v + f 2 = (1 − α ) v2 +
β = 1:
2
H ≥E−
0
f ( ξ ) dξ − (
2
f α. ) ∫
(10.355a–c)
The middle term is compared with the last term on the r.h.s. of (10.355c)
using L’Hôpital rule (I.19.35):
f ( x )
2
2f′ f
lim x = lim = 2 f ′ ( 0 ) = 2 lim x k ( x ) ′ = 2 k ( 0 ) > 0,
f (10.356a–d)
∫ f ( ξ ) dξ
x→ 0 x→ 0 x→ 0
where k(0) is the resilience of the spring at zero displacement, as for a lin-
ear restoring force. It follows that there exists a constant (10.357a) such that
(10.357b) holds:
mα c x x
( 1 − α ) v 2 + 1 −
m
f ( x ) ≤ c
∫ f ( ξ ) dξ ; H ≥
f ( ξ ) dξ ,
∫
2
c > 0:
0 2 2 0
(10.357a–c)
2 x
α < 1, :
mc
H > 0; f ( x) > 0 <
∫
0
f ( ξ ) dξ . (10.358a–e)
dg dg dg dx
g ≡ , g′ = : g = = g ′ v , (10.359a–d)
dt dx dx dt
84 Differential Equations and Applications
H = m v ( v + α f ) + f v + α m f ′ v 2 = − ( µ v + f )( v + α f ) + v ( f + α m f ′ v )
= − µ v2 + α m f ′ v2 − α f 2 − α µ f v .
(10.360a–c)
The inequality (10.354c) is applied to the last term on the r.h.s. of (10.359c)
leading to (10.360d):
µα 2 f2
H ≤ − µ v 2 + α m v 2 f ′ − α f 2 + v β + . (10.360d)
2 β
Introducing the upper bound (10.361a) for f ′ ( x ), which must be finite, leads
to (10.361b, c):
α mB αβ µ
B = max f ′ ( x ) :
H ≤ −µ v 2 1 − − − α f 2 1 − < 0.
µ 2 2β
(10.361a–c)
2 1 µ
β> ; α< < , (10.362a–c)
µ mB β 1 + mB
+
µ 2
(ii) a small enough α so that (10.362b), the first term on the r.h.s. of (10.361b),
is also negative; (iii) the value (10.362c) of α is compatible with (10.358a) since
it suffices to comply with the smallest of the two upper bounds. Thus, α
meeting (10.358a; 10.361c) is sufficient to ensure the conditions (10.353a, b) of
asymptotic stability, quod erat demonstrandum (QED).
which is in the cases when the generalized circular and hyperbolic cosine
and sine (subsections 9.4.12–9.4.20) do not apply directly.
The analysis is distinct for: (i)(ii) a simple m = −1 (double m = −2) pole
[E10.20.1(E10.20.2)], for which the origin is a regular singularity leading to
regular integrals (elementary power solutions); (iii)(iv) a pole of third m = −3
(higher m = −4 − 5....) order [E10.20.3(E10.20.4–E10.20.5)], for which the origin
is an irregular singularity and irregular integrals are needed.
m = − 1: x y ′′ ± y = 0, y ( x) = ∑x
k=0
k+a
c k ( a ) , (10.364a–d)
∞ ∞
∞ (10.365a–c)
= xa ∑n= 0
x n ( n + a )( n + a + 1) cn+ 1 ( a ) ± cn ( a ) .
From (10.365c) follows the recurrence formula for the coefficients (10.366a):
( n + a )( n + a − 1) cn ( a ) = cn− 1 ( a ) , n = 0: a± = 0 , 1, (10.366a–d)
and for (10.366b) the indices (3.366c, d). The higher index (10.366d) leads
(10.366a) to the coefficients (10.367a–c):
cn− 1 ( 1) c0 ( 1) c0 ( 1)
cn ( 1) = = () = ()
n n
. (10.367a, b)
( n + 1) n ( n + 1) (
n 2
n − 1) 2
...2 2
( n + 1)! n !
∞
( )n x n+ 1
c0 ( 1) = 1: y1± ( x ) = ∑
n= 0
n ! ( n + 1)!
, (10.368a, b)
Since the indices (10.366c, d) differ by an integer the second particular inte-
gral linearly independent from (10.368b) is a regular integral of the second
kind, second type (9.412) specified by:
∂
∞
y 0± ( x ) = lim
a→ 0 ∂ a ∑
a x a + n cn ( a ) , (10.369)
n= 0
cn ( a ) =
( )n c 0 ( a ) . (10.370)
( n + a )( n + a − 1)2 ... ( a + 1)2 a
∂ a
∞
( x )n
c0 ( a ) = 1: y 0± ( x ) = lim
a→ 0
x a +
∂a ∑ ,
( n + a )( n + a − 1)2 ... ( a + 1)2
n= 1
(10.371a, b)
leads to:
∞
( x )n
y ( x ) = lim x a 1 + a log x +
±
0 ∑ ( n + a ) ( n + a − 1)2 ... ( a + 1)2
a→ 0
n= 1
1 2 2
log x − − ... −
n+ a n+ a−1 a + 1
∞
( x )n log x − 1 − 2 ... − 2
= 1+ ∑
n= 1
n ! ( n − 1)! n n−1
∞
( x )n log x − ψ 1 + n − ψ n + 2 ψ 1
= 1+ ∑
n= 1
n ! ( n − 1)!
( ) ( ) ( )
∞
( x )n ψ 2 + m + ψ 1 + m − 2 ψ 1 ,
= 1 + log x y ±
1 ( x) − ∑
m= 0
(
m ! ( m + 1)!
) ( ) ( )
(10.372a–d)
Thus, the general integral (10.373b) of the generalized hyperbolic (circular) differen-
tial equation (10.364b)[(10.364c)] with parameter (10.364a) in (10.363a–c) is (stan-
dard CCLXXXV) a linear combination (10.373b) with arbitrary constants ( C0 , C1 ),
valid in the finite complex x-plane (10.373a) of: (i) a regular integral of the first kind
(10.368b); (ii) a regular integral of the second kind type two (10.372d) consisting
(10.374a) of the sum of: (ii-1) a preliminary function equal to unity; (ii-2) a logarithm
multiplying a function of the first kind:
∞
( x )n ψ 1 + n + ψ n ,
y 0± ( x ) = log x y1± ( x ) + y 0± ( x ) : y 0± ( x ) = 1 − ∑
n= 1
(
n ! ( n + 1)!
) ( )
(10.374a, b)
m = − 2: x 2 y ′′ y = 0, (10.375a–c)
y ( x) = xa : 0 = a ( a − 1) 1 = P2± ( a ) , 2 a± = 1 ± 1 ± 4 . (10.376a–c)
The upper sign leads to the roots (10.377a) and the general integral (10.377b):
5 5
2 a± = 1 ± 5: y ( x ) = x 1/2 C0 cosh log x + C1 sinh log x ,
2 2
(10.377a, b)
3 3
2 a± = 1 ± i 3: y ( x ) = x 1/2 C0 cos log x + C1 sin log x ,
2 2
(10.378a, b)
88 Differential Equations and Applications
and the lower sign leads to the roots (10.378a) and the general integral (10.378b).
Thus, the generalized hyperbolic (circular) differential equation (10.375b) [(10.375c)]
with a double pole (10.375a) in (10.363b) [(10.363c)] has (standard CCLXXXVI)
general integral (9.377b) [(9.378b)] where (C0 ,C1 ) are arbitrary constants. The cases
of (10.363a–c) remaining are poles of order higher than the second leading to
an irregular singularity at the origin and irregular integrals (E10.20.3).
has an irregular singularity at the origin (10.379c) [(10.379d)] and hence, the
solution has an essential singularity, which is involving unending descending
powers x − n. This suggests as change of variable the inversion relative to the ori-
gin (9.553a, b) ≡ (10.380a, b) leading (9.555c–e) the differential equation (10.380c):
1
ξ= , w ( ξ ) = y ( x ): 0 = ξ w ′′ + 2 w ′ ± ξ p w. (10.380a–d)
x
In the case of a triple pole (10.381a, b), the differential equation (10.381c, d) has
a regular singular at infinity:
m = −3; p = 0: ξ w ′′ + 2 w ′ ± w = 0, (10.381a–d)
∞ ∞
∑n= 0
ξ n+ a en ( a ) = w ( ξ ) = y ( x ) = ∑x
n= 0
− n− a
e n ( a ) , (10.382a, b)
∞ ∞
∞ (10.383a–c)
= ξa ∑ n= 0
ξ n ( n + a + 1)( n + a + 2 ) en+ 1 ( a ) ± en ( a ) ,
Examples 10.1 to 10.20 89
setting (10.384b) leads to the índices (10.384c, d). The recurrence formula
(10.384a) may be applied iteratively (10.385b, c) up to the zero order coefficient
to which may be given the value unity (10.385a):
en − 1 ( a ) () n
e1 ( a ) = 1: en ( a ) =
= .
( n + a )( n + a + 1) ( n + a + 1)( n + a )2 ... ( 2 + a )2 ( 1 + a )
(10.385a–c)
∞ ∞
∑ n(!( n +) 1)! ,
−n
∑x
x
y ±
0 ( x) = −n
en ( 0 ) = (10.386a, b)
n= 0 n= 0
∂
∞
w−±1 ( ξ ) = lim
a→−1
( a + 1) ξ a
∂a
∑ξ
n= 0
n
en ( a )
∞
( ξ )m
= lim
∂
a→ − 1 ∂ a
ξ
a
a + 1 +
∑
m= 1
( m + a + 1)( m + a ) ... ( a + 2 )
2 2
∞
( ξ )m
= lim ξ a 1 + ( a + 1) log ξ +
a→ − 1
∑ ( m + a + 1)( m + a )2 ... ( a + 2 )2
m= 1
1 2 2
log x − − − ... −
m+ a+1 m+2 a + 2
∞
( ξ )m
= ξ −1 + ξ −1 ∑
m= 1
m ! ( m − 1)!
1
log ξ − −
2
m m−1
− ... − 2
∞
( ξ )n
=ξ +−1
∑n= 0
log ξ − ψ ( 2 + n) − ψ ( 1 + n) + 2 ψ ( 1) ,
n ! ( n + 1)!
(10.387a–e)
90 Differential Equations and Applications
m = − 3: x 3 y ′′ ± y = 0, (10.388a–c)
has the general integral (standard CCLXXXVII) valid outside the origin (9.389a),
specified by a linear combination (9.389b) with arbitrary constants ( C0 , C1 ):
of asymptotic regular integrals of the first kind (10.386b) and the second kind
(10.387e) ≡ (10.390a):
∞
( x )− n ψ 2 + n + ψ 1 + n ,
y −1 ( x ) = x − log x y 0 ( x ) − y 2 ( x ) ; y 2 ( x ) =
± ± ± ±
∑n= 0
(
n ! ( n + 1)!
) ( )
(10.390a, b)
the asymptotic regular integral of the second kind (10.390a) is of the second type and
consists of: (i) a preliminary function x; (ii) the integral of the first kind (10.386b)
multiplied by a logarithm; (iii) the complementary function (10.390b). The only
remaining case of solution of the generalized circular (hyperbolic differen-
tial equation (10.363b) [(10.363c)] is the case of a pole of order more than three
that is solvable in terms of the generalized circular (hyperbolic) functions of
the inverse variable (E10.20.4).
m = − 4, −5,... : y ′′ ± x m y = 0, (10.391a–c)
leads for the inverse variable (10.379b; 10.380a–d) to the differential equation
(10.392a):
0 = ξ w ′′ + 2 w ′ ± ξ − m− 3 w = ( ξ w )′′ ± ξ − m− 4 ξ w , (10.392a, b)
1
− m − 4 = 0, 1, 2 ,... ∈|N 0 : ξ w (ξ) = y( x), (10.393a, b)
x
Examples 10.1 to 10.20 91
1 1
y + ( x ) = x B0 cos ; − m − 4 + B1 sin ; − m − 4 , (10.396)
x x
1 1
y − ( x ) = x E0 cosh ; − m − 4 + E1 sinh ; − m − 4 . (10.397)
x x
1 1
x 4 y ′′ + y = 0: y ( x ) = x B0 cos + B1 sin , (10.398a, b)
x x
1 1
x 4 y ′′ − y = 0: y ( x ) = x E0 cosh + E1 sinh . (10.399a, b)
x x
with arbitrary constants ( C0 , C1 ) of the two Airy functions of the inverse variable 1/x
multiplied by the variable x:
1 1
x 5 y ′′ + y = 0: y ( x ) = x C0 Ai + C1 Bi . (10.400a, b)
x x
Conclusion 10
Two examples of integral curves or streamlines of a flow: (Figure 10.1) all
passing through the origin; (Figure 10.2) not passing through the origin,
except for the six asymptotes. The cusped (smooth) parabola is the envelope
of a one-parameter family of straight lines [Figure 10.3 (10.4)]. A beam, which
is a bar under axial tension, can be subject to 16 possible combinations of
loads and supports, of which four are illustrated (Figures 10.5a–d), using
four types of (i) support: clamped, pinned, clamped-pinned or cantilever;
(ii) loads: concentrated torque or force or uniform or linearly increasing shear
stress. The linear bending of a circular plate with a circular hole allows for
eight combinations (Table 10.1), of which five are distinct (Table 10.2), namely:
(Figure 10.6a/e): clamped/supported at both boundaries; (Figure 10.6b/c)
clamped at the outer boundary and supported/free at the inner boundary;
(Figure 10.6d) supported at the outer boundary and free at the inner bound-
ary. For an elastic membrane without shear stress and under uniform ten-
sion parallel to the supports and with uniform (non-uniform) traction at the
supports [Figure 10.7a (10.8a)], there are modes of all orders n = 1, 2 ,..., for
Examples 10.1 to 10.20 93
94
Classification 10.1 95
to CCI. The fifth group relating to the special functions is too extensive to
cover in a single volume; the methods of solution of linear differential equa-
tions with variable coefficients are presented in detail and are illustrated
with generalized Bessel and other special functions. Classification 10.1 with
500 standards on the methods of solution of differential equations is comple-
mented by Classification 10.2, with detailed application to 500 physical and
engineering problems.
∑
dn y
An = B ( x ) . (10.401)
n= 0
dx n
0 = PN ( a ) = ∑ A a . (10.402)
n= 0
n
n
B ( a ) = B e x . (10.403)
B( x) = ∑B
m= 0
m x m ; (10.409)
B ( x ) = e bx ∑B
m= 0
m x m ; (10.411)
∑
dn y
An x n = B ( x ) , (10.417)
n= 0
dx n
0 = PN ( a ) = ∑A a ( a − 1)...( a − n + 1) , (10.418)
n= 0
n
B ( x ) = x a ; (10.419)
98 Differential Equations and Applications
B( x) ≠ 0
B( x) = ∑B
m= 0
m log m x ;(10.425)
∑A y
k=0
n n+ k = Bn . (10.432)
Bn = b n ; (10.433)
∑
d dn y
L y ( x) ≡ An ( x ) = B ( x ) , (10.437)
dx n= 0
dx n
d
L y m ( x ) = 0 , W ( y1 , y 2 ,..., y N ) ≠ 0: y ( x ) =
dx
∑
m= 1
Cm ( x ) y m ( x ) ,
(10.438a–c)
with the coefficients not constant, but rather functions of the independent
variable x, to be determined.
Classification 10.1 101
d
L G ( x ; ξ ) = δ ( x − ξ ) , (10.439)
dx
b
B ∈C ( a, b ) : y ( x) =
∫ a
B ( ξ ) G ( x ; ξ ) d ξ . (10.440a, b)
N N
πmx
∑ ∑ B exp i
dn y
B ∈C ( − L, + L ) : An = B( x) ≡ m . (10.441a, b)
n= 0
dx n n= 0
L
∑A
dn y +∞
B ∈L1 ∩ F (−∞ , + ∞):
n= 0
n
dx n
= B( x) ≡
∫ −∞
B ( k ) e ikx dk. (10.442a–c)
102 Differential Equations and Applications
∑
dn y ∞
B ∈L1 (0, L) ∩ V (L , + ∞):
n= 0
An
dx n
= B( x) =
∫
0
B ( s ) e − sx ds. (10.443a–c)
dy
y′ ≡ : F ( x , y ; y ′ ) = 0 (10.444a, b)
dx
dy X ( x )
X , Y ∈E (|R ) : = ; (10.445a, b)
dx Y ( y )
y ′ = P ( x ) y 2 + P ( x ) y + Q ( x ) . (10.450)
u′′ + A u′ + B u = 0. (10.451)
ab R′ ( x )
y′ = + A + B − y + R ( x ) y ; (10.452)
2
R ( x) R ( x)
a2 R′ ( x )
y′ = + 2 a − y + R ( x ) y ; (10.453)
2
R ( x) R ( x)
ab a + b − 1 R′ ( x )
y′ = + − y + R ( x ) y ; (10.454)
2
x R ( x)
2
x R ( x)
a2 2 a R′ ( x )
y′ = + − y + R ( x ) y . (10.455)
2
x R ( x) x
2
R ( x)
104 Differential Equations and Applications
λ
a, b ∈|R: y ′ = a x2λ + y + b y 2 . (10.456a–c)
x
y
y ∈B (|R ) : y ′ = f ; (10.457a, b)
x
y
y′ = a + b . (10.458)
x
X = {X , Y } : 0 = X ( x , y ) dx + Y ( x , y ) dy = X . dx . (10.459a–c)
∂Y ∂X
Ω ( x, y ) ≡ − = 0: X ( x , y ) dx + Y ( x , y ) dy = dΦ ( x , y ) ; (10.460a, b)
∂x ∂y
Ω ( x , y ) ≠ 0: X ( x , y ) dx + Y ( x , y ) dy = λ ( x , y ) dΦ ( x , y ) ; (10.461a, b)
Ω ( x, y ) x
Y ( x, y )
≡ f ( x ): X ( x , y ) dx + Y ( x , y ) dy = exp
∫ f ( ξ ) d ξ dΦ ( x , y ) ;
(10.462a, b)
Classification 10.1 105
Ω ( x, y ) y
X ( x, y )
≡ g ( y ): X ( x , y ) dx + Y ( x , y ) dy = exp
∫ g ( η) d η dΦ ( x , y ) ;
(10.463a, b)
Ω ( x, y )
≡ h ( x + y ): X ( x , y ) dx + Y ( x , y ) dy
X ( x, y ) − Y ( x, y )
(10.464a, b)
x+ y
= exp
∫ h ( z ) dz dΦ ( x , y ) ;
Ω ( x, y )
≡ j ( xy ) : X ( x , y ) dx + Y ( x , y ) dy
x X ( x, y ) − y Y ( x, y )
(10.465a, b)
x+ y
= exp
∫ j ( u) du dΦ ( x , y ) ;
x2 Ω ( x, y ) y
≡ k : X ( x , y ) dx + Y ( x , y ) dy
x X ( x, y ) − y Y ( x, y ) x
(10.466a, b)
x+ y
= exp
∫ k ( v ) dv dΦ ( x , y ) ;
λ = x a − 1 y b − 1 : a y d x + b x d y + x α− a y β− b ( α y d x + β x dy ) =
1
λ
(
d x a y b + xα y β .)
(10.467a, b)
106 Differential Equations and Applications
X = {X , Y , Z} : 0 = X ( x , y , z ) dx + Y ( x , y , z ) dy + Z ( x , y , z ) dz = X . dx .
(10.468a, b)
n, m = 1,..., N : 0= ∑ X (x
n= 0
n m ) dxn . (10.477a, b)
Ωmn =
∂X n ∂X n
−
∂ xm ∂ xm
=0 ⇔ ∑ X (x
n= 1
n m ) dxn = dΦ ( xm ) ; (10.478a, b)
H mns ≠ 0 ⇔ ∑ X (x
n= 0
n m ) dxn ≠ λ ( xm ) dΦ ( xm ) . (10.480a, b)
0= ∑
n= 0
X n ( xm ) dxn → 0= ∑ Y (y
m= 1
m m ) dyn . (10.481a, b)
108 Differential Equations and Applications
y y
X ( x , y ) = x q X 1, , Y ( x , y ) = x q Y 1, ; (10.482a, b)
x x
∂Y ∂X X ( x , y ) dx + Y ( x , y ) dy
≠ : dΦ ( x , y ) = ; (10.483a, b)
∂x ∂y x X ( x, y ) + y Y ( x, y )
∂X ∂Y
q ≠ −1; ≡ : x X ( x , y ) + y Y ( x , y ) = C. (10.484a–c)
∂y ∂x
x x
X n ( x1 ,...., xN ) = ( x1 ) X 1, 2 ,..., N ;
q
m, n = 1,..., N : (10.485a, b)
x1 x1
N
N
∂X n ∂X n
≠
∂ xm ∂ xm
: ∑n= 1
X n ( xm ) dxn =
n= 1
∑
X n ( xm ) xn dΦ ( xm ) ;
(10.486a, b)
q ≠ − 1;
∂X n ∂X n
=
∂ xm ∂ xm
: ∑ X (x
n= 1
n m ) xn = C. (10.487a–c)
0= ∑
i1 ,.., iM = 1
X i1 ...i dx i ... dx iM . (10.488)
M
Classification 10.1 109
0 = A ( x , y ) ( dx ) + B ( x , y )( dy ) + D ( x , y ) dx dy . (10.489)
2 2
0 = A ( dx ) + B ( dy ) + C ( dz ) + 2 D dx dy + 2E dx dz + 2 F dy dz , (10.490a)
2 2 2
A BC + 2 DE F = A F 2 + BE 2 + C D2 . (10.490b)
0 = P ( x , y ) C 2 + Q ( x , y ) C + R ( x , y ) . (10.492)
110 Differential Equations and Applications
dy dy
h ∈D ( R ) : y = x + h ; (10.493a, b)
dx dx
*standard C. D’Alembert differential equation (10.494c) involving (10.494a, b)
two differentiable functions (subsections 5.3.3–5.3.4; E10.8.1):
dy dy
g , h ∈D ( R ) : y = xg + h . (10.494a–c)
dx dx
dy
( )
h ∈D |R 2 : y = h x, ;
dx
(10.495a, b)
dy
( )
g ∈D |R 2 : x = g y , .
dx
(10.496a, b)
M m
0= ∑
m= 0
dy
Am . (10.497)
dx
Classification 10.1 111
dΦ = ∑X
m= 1
n dxn . (10.498)
Ψ ( X 1 , x2 ,...., xN ) = Φ ( x1 , x2 ,...., xN ) − x1 , X 1 : dΨ = − x1 d X 1 + ∑X
m= 2
m dxm .
(10.499a, b)
dy dY dY
0 = F x; y, = F ;X − Y , X ; (10.500a, b)
dx dX dX
*standard CVI. Special case:
dy dy dY
0 = Φ x − y − x Ψ = Φ( Υ ) − Ψ (X ) . (10.501a, b)
dx dx dX
dy d2 y
y′ ≡ , y ′′ ≡ 2 : 0 = F ( x ; y , y ′ , y ′′ ) . (10.502a–c)
dx dx
dy ′
0 = F x ; y , y ′, y ′ . (10.503)
dy
dy ′
0 = F ( y , y ′ , y ′′ ) = F y , y ′ , y ′ . (10.504)
dy
112 Differential Equations and Applications
*standard CIX. Omitting both the independent variable and slope (subsec-
tion 5.6.3):
0 = F ( y , y ′′ ) ; (10.505)
z ≡ y ′: 0 = F ( x , y ′ , y ′′ ) = F ( x ; z, z ′ ) ; (10.506a, b)
*standard CXI. Missing both the slope and the dependent variable (subsec-
tion 5.6.4; E10.9.1):
0 = F ( x ; y ′′ ) . (10.507)
dN y
y( N ) ( x ) ≡
dx N
(
: 0 = F x ; y ; y ′ , y ′′ ,....y( N ) . ) (10.508a, b)
( )
0 = F x , y( N ) ; (10.509)
z = y( N − 1) ( x ) : ( )
0 = F x , y( N − 1) , y( N ) = F ( x ; z, z ′ ) ; (10.510a–c)
w = y( N − 1) ( x ) : ( )
0 = F y( N − 1) , y( N ) = F ( w , w ′ ) . (10.511a–c)
u = y( N − 2 ) ( x ) : ( )
0 = F y( N − 2 ) , y( N ) = F ( u, u′′ ) ; (10.512a–c)
z ≡ y( N − P) ( x ) : ( ) ( )
0 = F x , y( N − P) ,..., y( N ) = F x ; z,..., z( P) . (10.513a–c)
Classification 10.1 113
d d
PN QM y ( x ) = 0 ; (10.514)
dx dx
(
0 = F x , y , y ′ ,..., y( N ) = ) dp
dx P {(
G x ; y , y ′ ,..., y ( )
N−p
) }. (10.515)
L. f. Higher-Order Homogeneous Differential Equations (section 5.9):
*standard CXIX. Homogeneous differential equation of the first kind order
N (subsection 5.9.1; E10.9.6):
y y ′ y ′′ y( N )
0 = F , y ′ , y ′′ ,..., y( N ) = G x , , ,..., ; (10.516a, b)
x x y y
( )
0 = F y , y ′x , y ′′x 2 ,..., y( N ) x N . (10.517)
dy m
m = 1,..., M: = Ym ( x ; y1 ,...., y M ) ; (10.518a, b)
dx
m, s = 1,..., M ; N m , s ∈|N : ( )
0 = Fs x ; y m , y m′ , y m′′ ,..., y m( m , s ) ; (10.519a–c)
N
114 Differential Equations and Applications
dxn
m, s = 1,..., N : = X n ( xm ) ; (10.520a, b)
dt
M Nm , r
∑∑
dk y
m, s = 1,..., M ; N m , r ∉|N : Am , r , k ( x ) = Bm ( x ) , (10.521a–c)
m, r =1 k = 0
dx k
M Nm , r
∑∑
dk y
m = 1,..., M ; N m , r ∉|N : Am , r , k = Bm ( x ) ; (10.522a–c)
m, r =1 k = 0
dx k
M Nm , r
∑∑
dk y
m = 1,..., M ; N m , r ∉|N : Am , r , k x k = Bm ( x ) . (10.523a–c)
m, r =1 k = 0
dx k
Classification 10.1 115
Nm , r
PN ( a ) = Det { Pm , r ( a )} , Pm , r ( a ) ≡ ∑A
k=0
m, r , k a k . (10.524a, b)
Pm ( x ) = Gm e bx ; (10.525)
Nm , r
QN ( a ) = Det (Qm , r ) , Qm , r ( a ) ≡ ∑A
k=0
m, r , k a ( a − 1) ... ( a − k + 1) . (10.529a, b)
Bm ( x ) = Gm x b ; (10.530)
M Nm , r
m = 1,..., M ; N m , r ∈|N : Bm , = ∑∑ A
m, =1 k = 0
m, r , k y r , k + ; (10.534a–c)
M Nm , r
∆ k yr , = yr , + k : Bm , = ∑∑ A
m, =1 k = 0
m, r , k ∆ k y r , k + ; (10.535a, b)
118 Differential Equations and Applications
Nm , r
RN ( ∆ ) = Det {Rm , r ( a )} , Rm , r = ∑A
k=0
m, r , k ∆ k . (10.536a, b)
*standard CLXIV. The roots of the characteristic polynomial specify the nat-
ural sequences (subsection 7.8.3);
*standard CLXV. Natural sequences corresponding to a single root of the
characteristic polynomial (subsection 7.8.4);
*standard CLXVI. Natural sequences corresponding to a multiple root of the
characteristic polynomial (subsection 7.8.4);
*standard CLXVII. General solution for characteristic polynomial with all
roots distinct (subsection 7.8.5);
*standard CLXVIII. General solution for characteristic polynomial with dis-
tinct multiple roots (subsection 7.8.5);
*standard CLXIX. Arbitrary constants determined from starting conditions
(subsections 7.8.5 and 7.8.7);
*standard CLXX. Coefficients of natural sequences determined from com-
patibility conditions (subsections 7.8.5–7.8.6);
*standard CLXXI. Diagonal system of natural sequences corresponding to
single roots of the characteristic polynomial (subsections 7.8.8 and 7.8.10);
*standard CLXXII. Lower triangular system of natural sequences corre-
sponding to multiple roots of the characteristic polynomial (subsection 7.8.8);
*standard CLXXIII. Block-lower triangular system of natural integrals corre-
sponding to distinct multiple roots of the characteristic polynomial (subsec-
tion 7.8.9).
Bm , = Gm b ; (10.537)
N M
∑n= 0
∂2 F
∂ xn2
+ ∑A
m= 0
m
∂2 F
∂t m
= 0. (10.541)
∇ ≡ 2
∑
n= 1
∂2 F
∂ xn2
: ∇ 2 F + k 2 F = 0. (10.542a, b)
M +∞
k ≡ 2
∑A
m= 0
m (−i ω ) m
: F ( xn , t ) =
∫ F ( x , t ) e
n
− iωt
dω . (10.543a, b)
−∞
( )
y ( x ) = J ν ( x ) , Yν ( x ) ; ν2 = a 2 + ( N − 1/ 2 ) , x 2 y ′′ + ( N − 1) x y ′ + x 2 − a 2 y = 0 ;
2
(10.544a–c)
dy
= f ( x, y ; λ ) , y ( x0 ) = y 0 + ε: y = y ( x ; y 0 , λ , ε ) . (10.546a–c)
dx
dy m
= fm ( x , y n ; λ ) , y m ( x0 ) = y m ,0 + ε m y = y ( x ; y 0, m , λ , ε m ) .
m, n = 1,...., M:
dx
(10.547a–d)
dN y
n = 0,..., N − 1:
dN y
dx N
= F x ; dx N ; λ , ( )
y( n) ( x0 ) = y(0n) + ε n , y = y x ; y(0n) , λ , ε n .
(10.548a–d)
*standard CCXIV. Existence and unicity of solution (10.548d) of the N-th order
differential equation (10.548a, b) without parameter λ, with boundary condi-
tions (10.548c) without ε n = 0 perturbations (subsections 9.1.18 and 9.1.20);
*standard CCXV. Robustness theorem for the solution (10.548d) of the non-
linear N-th order differential equation (10.548a, b) without the parameter λ,
and with small pertubations ε n of the (10.548c) boundary condition (subsec-
tions 9.1.18 and 9.1.20);
*standard CCXVI. Uniformity theorem for the solution (10.548d) of the non-
linear N-th order differential equation (10.548a, b) with parameter λ, and
without ε n = 0 perturbation of (10.548c) the boundary conditions (subsections
9.1.18 and 9.1.20).
Classification 10.1 123
∑
dn y
An ( x ; λ ) = B ( x ; λ ) ; m = 0 ;..., N − 1:
n= 0
dx n (10.549a–d)
y( m) ( x0 ) = y(0m) + ε n , ( )
y = y x ; y(0m) , λ , ε m .
d2 d
m 2 + µ ( x ) + k ( x ) x ( t ) = 0 ; (10.553)
dt dt
d2 d
m 2 + µ ( t ) + k ( t ) x ( t ) = 0 ; (10.554)
dt dt
d2
2 + ω 0 1 + 2 h cos ( ω e t ) x ( t ) = 0. (10.555)
2
dt
Classification 10.1 125
∑
dy n
m, n = 1,..., M ; Amn ( x ) ≤ B < ∞: = Anm ( x ) y m ( x ) , (10.556a–c)
dx m= 1
*standard CCXXXV. Existence of at least one eigenvalue for the linear auton-
omous system of differential equations (10.556a, c) with bounded (10.556b)
coefficients (subsection 9.2.13);
*standard CCXXXVI. Existence of at least one eigenvalue (standard CCXXXV)
or equivalently existence of exponential positive upper (negative lower)
bounds (standards CCXXXIII–CCXXXIV) for the solution and first N – 1
derivatives of a linear differential equation (10.557c) of order N with bounded
(10.557a, b) coefficients (subsection 9.2.14):
∑
dn y
m = 0,..., N ; Amn ( x ) ≤ B < ∞: Anm ( x ) = 0. (10.557a–c)
m= 1
dx n
dy n
= Fn ( xm ( t )) , x0 = x ( 0 ) , x = x ( t ; x0 ) , (10.558d–f)
dx
*standard CCXXXIX. Upper bound for the solution (10.558f) of the autono-
mous differential system (10.558a–d) with the initial (10.558e) condition (sub-
sections 9.2.16–9.2.20);
126 Differential Equations and Applications
dN y dy d 2 y d N −1 y
(
F ∈C 2 |R N − 1 : ) dx N
= F y , ,
dx dx 2 ,....,
dx N − 1
. (10.559a, b)
n = 0,...., N ; An ( t ) = An ( t + τ ) : ∑
m= 0
An ( x )
dn x
dt n
= 0. (10.560a–c)
d2
I (t ) = I (t + τ ): 2 + I ( t ) x ( t ) = 0. (10.561a, b)
dt
∑
dn y
n = 0,...., N : An ∈A (|C ) : An ( x ) = 0. (10.562a–c)
n= 0
dx n
y ( x) = xa ∑
β= 0
gβ logβ x ∑c
n =− ∞
n ,β x n , (10.563)
involving the index a, its multiplicity α, the logarithmic factors gβ and the
coefficients cn,β of the Laurent series (subsection 9.4.6);
*standard CCLXV. Classification of analytic, regular, and irregular integrals
of an unforced linear first-order differential equation in the case of the point
at infinity that is excluded from the standard CCLXII (subsection 9.4.7);
*standard CCLXVI. Upper bound in modulus for the exponential of a com-
plex matrix whose elements are bounded in modulus (subsection 9.4.8);
*standard CCLXVII. – Existence of a solution bounded in modulus for a linear
first-order autonomous system with integrable coefficients (subsection 9.4.8);
*standard CCLXVIII. Analytic, regular, and irregular integrals of a linear
autonomous system of differential equations (subsection 9.4.8).
y ( x) = ∑
n= 0
cn x n = ∑
n= 0
x n ( n)
n!
y ( 0 ) ; (10.563a, b)
Classification 10.1 129
*standard CCLX. Particular case of the standard CCLIX for a linear second-
order differential equation in invariant form (10.564b) with analytic (10.564a)
coefficient (subsection 9.4.9):
∞
d2
I ( x) = ∑I
k=0
k xk : 2 − I ( x ) y ( x ) = 0 ;
dx
(10.564a, b)
*standard CCLXI. First method of solution of the standard CCLIX, with the
standard CCLX as example, specifying the MacLaurin series (10.563b) as the
analytic integral (subsection 9.4.10);
*standard CCLXII. Second alternative and equivalent method of solution to
the standard CCLXI specifying the recurrence formula for the coefficients cn
of the (10.563a) analytic integral (subsection 9.4.11).
y ′′ − x y = 0: y ( x ) = c0 Ai ( x ) + c1 Bi ( x ) ; (10.566a, b)
∑ a x
dy m
m, n = 1,..., N : Pm , n ∈A (|C ) : = m, n
+ pm , n ( x ) y n ( x ) .
dx n= 1
(10.569a–d)
x < R: y ( x) = xa ∑x
n= 0
n
cn ( a ) , (10.570a, b)
∑x
dn y
n = 0,..., N ; pn ∈A (|C ) : n
pn ( x ) = 0. (10.571a–c)
n= 0
dx n
d2 y dy
p , q ∈A (|C ) : x2 2
+ x p ( x) + q ( x ) y ( x ) = 0. (10.572a–c)
dx dx
0 = a ( a − 1) + a p ( 0 ) + q ( 0 ) = ( a − a+ )( a − a− ) ; (10.573a, b)
′′ + P ( x ) y1,2
y1,2 ′ ( x ) + Q ( x ) y1,2 ( x ) = 0, W ( y1 , y 2 ) = y1 y 2′ − y 2 y1′ ,
(10.574a–d)
x
W ( x ) = W ( x0 ) exp
∫ x0
p ( ξ ) dξ ;
µ
y ( x ) = J νµ ( x ) , Yνµ ( x ) , H µ(1,2
,ν ( x ):
)
( )
x 2 y ′′ + x 1 − x 2 y ′ + x 2 − ν2 y = 0.
2
(10.575a–e)
Classification 10.1 133
( )
y ( x ) = J ν ( x ) , Yν ( x ) , H ν(1,2 ) ( x ) : x 2 y ′′ + x y ′′ + y ′ + x 2 − ν2 y = 0 ; (10.576a–e)
*standard CCCXVI. In the further particular case of zero order ν = 0 the inte-
grals (10.577a–d) of the original Bessel differential equation (10.577e) of order
zero (subsections 9.5.22);
y ( x ) = J 0 ( x ) , Y0 ( x ) , H 0(1,2 ) ( x ) : x y ′′ + y ′ + x y = 0 ; (10.577a–e)
y ( x ) = jn ( x ) , y n ( x ) , hn(1,2 ) ( x ) : x 2 y ′′ + 2 x y ′ + x 2 − n ( n + 1) y = 0, (10.578a–e)
d2 y dy
P, Q ∈A ( 0 < x − x0 < R ) : 2
+ x P ( x) + Q ( x ) y = 0; (10.579a, b)
dx dx
*standard CCXXIV. Example of the existence of at most one regular integral
near an irregular singularity (subsection 9.6.2).
M ∞
1
Ω =
x ∑m= 1
Ωm
xm
:
1
y ( x ) = exp Ω x a
x
∑x c ( a) ;
m= 0
n
n (10.580a, b)
136 Differential Equations and Applications
*standard CCCXXVI. Analytic, regular, and irregular integrals near the point
at infinity of a linear second-order differential equation (subsection 9.6.4);
*standard CCCXXVII. Comparison of an asymptotic series with an asymp-
totic expansion (subsection 9.6.5);
*standard CCCXXVIII. Properties of asymptotic expansions including (i) the
sum, (ii) the Cauchy rule for the product, (iii) term-by-term differentiation,
and (iv) integration (subsection 9.6.5).
+∞
y ( x ) = ( x − x0 )
a
∑ ( x − x ) c ( a) ; (10.581)
n =−∞
0
n
n
d2 z
∞
+ b0 +
dθ2 ∑b cos ( 2 θ) u(θ) = 0. (10.582)
=1
dy
∞
d2 y
0 = x2
dx 2
+x − b0 +
dx
∑b ( x
n= 1
n
n
)
+ x − n y = 0. (10.583)
*standard CCCLXV. Relation between the gamma (standard CCCII) and beta
(standard CCCLXIII) functions (subsection 9.8.3);
*standard CCCLXVI. Evaluation of integrals in terms of beta (standards
CCCLIII–CCCLV) or gamma (standards CCCII–CCCIII) functions (subsec-
tion 9.8.4).
∑( pn + qn x) dxn = 0 ;
dn y
y ( x) =
∫
L
v ( s ) e sx dx:
n= 0
(10.584a, b)
µ η 2 µν
η = x2 , ( )
y ( x ) = x± ν Ψ ± x2 : η Ψ ′′± + 1 ± ν −
Ψ ′± +
4 8
Ψ = 0;
(10.585a–c)
∑
dn y
n= 0
An ( x )
dx n
= 0: y ( x) =
∫ k ( x, s) v ( s) ds ; (10.586a–b)
Γ
y ( x ) = F (α ; γ ; x ) , G (α ; γ ; x ): x y ′′ + ( γ − x ) y ′ − α y = 0; (10.588a–c)
142 Differential Equations and Applications
N M
∏
n= 1
d
x
dx
− an − x
∏
n= 1
d
x
dx
− bm y ( x ) = 0 ; (10.589)
Classification 10.1 143
x 3 y ′′′ + ( 3 − x − α 1 − α 2 − α 3 ) x 2 y ′′
+ 1 − α 1 − α 2 − α 3 + α 1α 2 + α 1α 3 + α 2 α 3 − x ( 1 − β1 − β 2 ) x y ′ (10.590)
− ( α 1α 2 α 3 + β 1 β 2 x ) y = 0 ;
x 3 y ′′′ + ( 3 − β1 − β 2 − β 3 ) x − 1 x 2 y ′′
+ ( 1 − β1 − β 2 − β 3 + β1β 2 + β1β 3 + β 2β 3 ) x y ′ + α 1α 2 − 1 x y ′ (10.591)
− (β1β 2β 3 x − α 1α 2 ) y = 0 ;
y ( x ) = F ( α , β ; γ ; x ) , G ( α , β ; γ ; x ) : x ( 1 − x ) y ′′ + γ − ( α + β + 1) x y ′ − αβ y = 0,
(10.592a–c)
*standard CDIV. Derivation of the Gaussian hypergeometric differential equa-
tion (10.592c) from the Gaussian hypergeometric series (10.592a) (note 9.14);
*standard CDV. General integral of the Gaussian hypergeometric differential
equation (10.592c) as a linear combination of Gaussian hypergeometric func-
tion of the first kind (10.592a) when γ is not an integer (note 9.15);
*standard CDVI. General integral of the Gaussian hypergeometric differ-
ential equation (10.592c) in the case of γ-positive integer, which is excluded
from the standard CDV, when a Gaussian hypergeometric functions of the
second kind (10.592b) appears (note 9.16);
*standard CDVII. Gaussian hypergeometric function of the second kind
(10.592b) with γ = n, a positive integer, which appears in this case in the solu-
tion of the Gaussian hypergeometric differential equation (10.592a) (note 9.17);
*standard CDVIII. General integral of the Gaussian hypergeometric differen-
tial equation (10.592c) in the case when γ is zero or a negative integer, which
is excluded from the standards CDV and CDVI, and unlike (like) the stan-
dard CDV(CDVI) involves a Gaussian hypergeometric function of the second
kind (10.592b) (note 9.18);
*standard CDIX. Gaussian hypergeometric function of the second kind
(10.592b), which appears in the general integral of the Gaussian hypergeomet-
ric differential equation (10.592c) when γ is zero or a positive integer (note 9.18);
*standard CDX. Particular case of the Gaussian hypergeometric function of
the second kind (10.592b) for γ = 1 (note 9.19);
*standard CDXI. Simplification of the Gaussian hypergeometric function of
the second kind (10.592b) when neither α nor β are negative integers (note 9.19);
*standard CDXII. Simplification of the Gaussian hypergeometric functions
of the first kind (10.592a) or second kind (10.592b) when α or β is a negative
integer (note 9.19);
*standard CDXIII. Summary of the general integrals of the Gaussian hyper-
geometric differential equation (10.592c) involving Gaussian hypergeometric
functions of the first kind (10.592a) or second kind (10.592b) for all cases of
values of the parameters ( α , β, γ ) (note 9.19).
y ( x ) = Pν ( x ) , Qν ( x ) : (1 − x ) y ′′ − 2 x y ′ + ν ( ν + 1) y = 0 ; (10.593a–c)
2
Classification 10.1 145
(10.594a–d)
y = 0,
y ( x ) = Pνµ ( x ) , Rνµ ( x ) :
(1 − x ) y ′′ − 2 x y ′ + ν ( ν + 1) − 1 −µx
2
2
(10.595a–c)
(1 − x ) y ′′ − ( 2 + λ ) x y ′ + ν ( ν + 1) − 1 −µ x
2
y ( x ) = Aνµ, λ ( x ) , Bνµ, λ ( x ) : 2
2 y = 0;
(10.596a–c)
y ( x ) = Pν , λ ( x ) , Rν , λ ( x ) : (1 − x ) y ′′ − ( 2 + λ ) x y ′ + ν ( ν + 1) y = 0 ; (10.597a–c)
2
y ( x ) = I νa ,b ( x ) : (1 − x ) y ′′ + b − a − ( a + b + 2 ) x y ′ + ν ( ν + a + b + 1) y = 0.
2
(10.598a, b)
y ( x ) = Tν ( x ) , Sν ( x ) : (1 − x ) y ′′ − x y ′ + ν
2 2
y = 0 ; (10.599a–c)
y = ( x ) = U ν ( x ) , Vν ( x ) : (1 − x ) y ′′ − 3 x y ′ + ν ( ν + 2 ) y = 0 ; (10.600a–c)
2
y ( x ) = Lν ( x ) , N ν ( x ) : x y ′′ + ( 1 − x ) y ′ + ν y = 0. (10.601a–c)
y ( x ) = Lµν ( x ) : x y ′′ + ( 1 + µ − x ) y ′ + ν y = 0. (10.602a, b)
y ( x ) = H ν ( x ) , Iν ( x ): y ′′ − 2 x y ′ + 2 ν y = 0. (10.603a–c)
w = ( θ ) = H ( b1 ,....., bM ; cos θ ) :
d2 w
dθ2
+ ∑b
m= 1
m cos ( m θ ) = 0. (10.604a, b)
y ( x ) = M ( a, q ; x ) : (1 − x ) y ′′ − x y ′ − x y ′ + ( a + 2 q − 4 q x ) y = 0 ;
2 2
(10.605a, b)
*standard CDLVII. Relation between the Mathieu function (10.605a) and the
extended Gaussian hypergeometric functions (10.606a) (note 9.35);
*standard CDLVIII. General integral of the extended Gaussian hypergeo-
metric differential equation (10.606b) with upper ( α , β ), lower γ, and asymp-
totic ( A1 ,..., AM ) parameters as a linear combination of extended Gaussian
hypergeometric functions of the first kind (10.606a) for γ not an integer
(notes 9.35–9.36):
y ( x ) = F ( α , β ; γ ; A1 ,..., AM ; x ) :
M
(10.606a, b)
x ( 1 − x ) y ′′ + γ − ( α + β + 1) x y ′ − α β +
∑m= 1
Am x m y = 0.
M
y ( x ) = F ( α , β ; γ ; A1 ,..., AM ; x ) : x y ′′ + ( γ − x ) y ′ − α +
∑A
m= 1
m x m y = 0;
(10.607a, b)
150 Differential Equations and Applications
y ( x ) = J ± ν ( A1 ,..., AM ; B1 ,..., BM ; x ) :
L
M
(10.608a, b)
x 2 y ′′ + x 1 +
∑
= 2
B x y ′ + x 2 − ν2 +
∑
m= 1
Am x m y = 0.
x
y ′′ + P ( x ) y ′ + Q ( x ) y = 0, R ( x ) = exp
∫ p ( ξ ) d ξ :
(10.609a, b)
P′ P2
( R y ′ )′ + RQ y = 0 ; z = y R: z ′′ + Q −
2
+
4
z = 0. (10.610a–c)
y ( x ) = F (α ; γ ; x ): (e −x
)
x γ y ′ ′ − α e − x x γ y = 0, (10.611a, b)
1 α γ γ γ2
z ( x ) = e − x/2 x γ /2 F ( α ; γ ; x ) : z ′′ − + − − + z = 0 ; (10.612a, b)
4 x 2x 2x
2
4 x 2
1 k m2 1
y ( x ) = Wk , m ( x ) : y ′′ − − + 2 − y = 0 ; (10.613a, b)
4 x x 4 x 2
Classification 10.1 151
*standard CDLXVI. Relation between the Whittaker (10.613a) and the gener-
alized Bessel (10.575a) function (note 9.41);
*standard CDLXVII. Self-adjoint (10.614a, b) [invariant (10.615a, b)] form of the
generalized Bessel differential equation (10.575e) (note 9.41):
′
µ x2 µ x2
y ( x ) = J νµ ( x ) : x exp −
y
4
′ 2 2
+ x − ν exp −
(
4
y = 0 , ) (10.614a, b)
−µ x 2 µ ν2 1 µ µ x2
z ( x ) = x exp Jν ( x ): z ′′ + 1 − 2 + + 1− z = 0.
2 x 4x 2
2 8
(10.615a, b)
′
y ( x ) = F (α , β; γ ; x ): x (1 − x )
γ y ′ − α β x γ ( 1 − x )
α+β+ 1− γ 1+α+β− γ
y = 0;
(10.616a, b)
(1+α+β− γ )/2
z ( x ) = x γ /2 ( 1 − x ) F (α, β; γ ; x ):
γ ( 2 − γ ) 1 − ( α + β − γ )2 γ ( 1 + α + β − γ ) − 2 α β (10.617a, b)
z ′′ + + + z = 0.
4 x
2
4 1 − x2 ( )
2 x (1 − x )
′ µ2
(
1 − x2 ) (
y ′ + 1 − x 2 )
1/2 +λ λ− 1/2
y ( x ) = Pνµ, λ ( x ) : ν ( ν + 1) − 1 − x 2 y = 0 ;
(10.618a, b)
( )
1/4 +λ/2
z ( x ) = 1 − x2 Pνµ, λ ( x ) :
′
(
y ( x ) = P ν , λ ( x ) : 1 − x 2 ) (
y ′ + ν ( ν + 1 + λ ) 1 − x 2 )
1/2 +λ λ− 1/2
y = 0 ; (10.620a, b)
152 Differential Equations and Applications
ν ( ν + 1 + λ ) x 2 ( 1 + λ/2 )( 3 + λ/2 )
( )
1/4 +λ/2
z ( x ) = 1 − x2 P ν , λ ( x ) : z ′′ + + z = 0.
1 − x2
1 − x 2 2
(
)
(10.621a, b)
y ( x ) = I νa ,b ( x ) : { (1 − x ) 2 1+ a/2 + b/2
exp ( b − a ) arg tanh x y ′ }′
( )
a/2 + b/2
+ 1 − x2 exp ( b − a ) arg tanh x y = 0 ;
(10.622a, b)
b−a
( )
1/2 + a/4 + b/4
z ( x ) = 1 − x2 exp arg tanh x I νa ,b ( x ) :
2
b − a 2 − 2x b − a 4 + a + b − x2 a + b + 2 a + b + 6
z ′′
−
( ) ( )( ) ( )( )
(10.623a, b)
4 1− x 2 2
( )
a + b + 2 − ν ( ν + a + b + 1)
+ z = 0.
1 − x2
y ( x ) = Tν ( x ) : ( ′
)
1 − x 2 y ′ + ν2 1 − x 2 y = 0 ; (10.624a, b)
1/2 + ν2 5 x 2
( )
1/4
z ( x ) = 1 − x2 Tν ( x ) : z ′′ + − z = 0 ; (10.625a, b)
1− x
2
4 1 − x2
′′
(
1 − x2 ) y ′ + ν ( ν + 2 )( 1 − x ) y = 0 ;
3/2
y ( x ) = U ν ( x ):
3/2
(10.626a, b)
3/2 + ν ( ν + 2 ) 21/9
( )
3/4
z ( x ) = 1 − x2 U ν ( x ): z ′′ + − z = 0 ; (10.627a, b)
1 − x2 1 − x 2 ( )
Classification 10.1 153
(e )
′
y ( x ) = Lµν ( x ) : − x 1+µ
x y ′ + ν e − x x 1+µ y = 0 ; (10.628a, b)
1 + µ + ν 1 − µ2
z ( x ) = e − x/2 x 1/2 +µ/2 Lµν ( x ) : z ′′ + + z = 0. (10.629a, b)
x 2 x2
( ) ( )
′
y ( x ) = H ν ( x ): exp − x 2 y ′ + 2 ν exp − x 2 y = 0 ;
(10.630a, b)
x2
z ( x ) = exp − H ν ( x ) :
2
( )
z ′′ + 1 + 2 ν − 4 x 2 z = 0. (10.631a, b)
( ) 1
y ( x ) = x a J νµ c x b : x 2 y ′′ + 1 − 2 a − µ b c 2 x 2 b x y ′
2
(10.632a, b)
2 2 2 2 2 2b µa
+ a − b ν + b c x 1 + y = 0.
2 b
µ x2
y ( x ) = J νµ ( i x ) : x 2 y ′′ + 1 +
2
(
x y ′ − x 2 + ν2 y = 0. ) (10.633a, b)
µx
y ( x ) = J νµ ( x ):
x 2 y ′′ + 1 −
4
1 2
4
2
(
x y ′ + x − ν y = 0. ) (10.634a, b)
154 Differential Equations and Applications
(
y ( x ) = J 2µν 2 i x : ) ( )
x 2 y ′′ + ( 1 + µ x ) x y ′ − x + ν2 y = 0. (10.635a, b)
( )
y ( x ) = J νµ x 2 : ( ) ( )
x 2 y ′′ + 1 − µ x 4 x y ′ + 4 x 4 − ν2 y = 0. (10.636a, b)
µ ν 2ν µ
( )
y ( x ) = x ν J νµ x ν : x 2 y ′′ + 1 − 2 ν −
2
x x y ′ + ν2 1 − ν2 1 − x 2 ν 1 + y = 0.
2
(10.637a, b)
µ x2 µ ν
y ( x ) = x ν J νµ ( x ) : x y ′′ + 1 − 2 ν − y ′ + x2 1 + y = 0. (10.638a, b)
2 2
µq 3 µ
( )
y ( x ) = x J pµ− 1/2 i x q : x 2 y ′′ +
2
x y ′ − p ( p − 1) + q 1 + x 2 y = 0.
4
(10.639a, b)
µq 2 µ
y ( x ) = x 1/2 + p J1/2
µ
(
+p i x q : )
x y ′′ − 2 p −
2
2
1
x y ′ − q 1 + p + x y = 0.
2
(10.640a, b)
q p µ q 2 p−1 µ 2 p−2
y ( x ) = x J1/
µ
( 2 p ) i p x : y ′′ + 2 p x y′ − q 1 + x y = 0. (10.641a, b)
4 p
Classification 10.1 155
µ 2 3/2 µ 2 µ
y ( x ) = x J1/3 x : y ′′ − x y ′ + 1 + x y = 0. (10.642a, b)
3 3 6
2
y ( x ) = x J1/3 x 3/2 = c1 A i ( − x ) + c2 B i ( − x ) : y ′′ + x y = 0. (10.643a–c)
3
f ′ g ′′ g ′ µ g ′g
y ( x ) = f ( x ) J ν ( g ( x )) :
y ′′ − 2 + − − y′
f g′ g 2
(10.644a–c)
f ′ f ′ g ′′ g ′ µ g ′g f ′′ g′2
+ 2 + − − − + g ′ 2 − ν2 2 y = 0.
f f g′ g 2 f g
g ( x) µ µ g ′g
y ( x) = J ν ( g ( x )) : y ′′ + y′
g′ ( x) 2
g ′′′ 3 g ′′ 2 1 2 g′
2
µ
+ + g′2 −
4g′ 2
+
− ν
2
+ ( )
g ′′g − g ′ 2 y = 0.
2g 4 g 4
(10.645a, b)
g ′′ g ′ µ g ′g
y ( x ) = g ( x ) J νµ ( g ( x )) : + ( 2 a − 1) −
a
y ′′ + y′
g′ g 2
(10.646a, b)
g′2
µa
(
+ a 2 − ν2 ) g2
+ g′2 1 −
y = 0.
2
156 Differential Equations and Applications
( )
f ( x ) = x a exp α xβ , g ( x ) = c xb , (10.647a, b)
µ b c 2 2b
( ) ( )
y ( x ) = x a exp α xβ J νµ c xb : x 2 y ′′ + 1 − 2 a − 2 α β xβ −
2
x x y′
µ a 2b µ b c 2 2b
+ a 2 − b 2 ν2 + b 2 c 2 1 + x + α β ( 2 a − β ) xβ + α β xβ α β xβ + x y = 0,
2b 2
(10.647c, d)
µ b c 2 2b
( ) ( )
y ( x ) = x a exp α x b J νµ c x b : x 2 y ′′ + 1 − 2 a − 2 α b x b −
2
x x y′
µ b c2
( )
+ a 2 − b 2 ν2 + α b ( 2 a − b ) x b + b 2 c 2 + α 2 x 2 b +
2
( )
a + α b x b x 2 b y = 0.
(10.648a, b)
µ b c 2 2b
( )
y ( x ) = x a exp ( ± i c x ) J νµ c x b : x 2 y ′′ + 1 − 2 a 2 i b c x b −
2
x x y′
µ b c2
+ a 2 − b 2 ν2 ± i b c ( 2 a − b ) x b +
2
( )
a ± i c b x b x 2 b y = 0.
(10.648c, d)
µ b c 2 2 bx
( )
y ( x ) = x − a J νµ c e b x : x 2 y ′′ + 2 a +
2
x e x y′
(10.649a, b)
µ a
+ a ( a − 1) − b 2 ν2 − c 2 e 2 bx 1 + 2
x y = 0.
2 b
Classification 10.1 157
µ x 2 2/x µ x 2/x 2
( )
y ( x ) = x J νµ e 1/ x : x 4 y ′′ −
2
e y ′ + 1 +
2
e − ν y = 0. (10.650a, b)
1 f ′ µx
y ( x ) = f ( x ) J νµ ( x ) : y ′′ + − 2 + y′
x f 2
(10.651a, b)
ν2 f ′ f ′ 1 µ x f ′′
+ 1 − 2 + 2 − − − y = 0.
x f f x 2 f
( )
a
*standard CDXCVII. Particular case of (10.651a, b) with f = x 2 − 1 leading
to (10.652a, b):
µ x2 2
( ) ( ) ( ) ( )
a 2
y ( x ) = x 2 − 1 J νµ ( x ) : x 2 x 2 − 1 y ′′ + ( 1 − 4 a ) x 2 − 1 + x − 1 ≠ x x2 − 1 y ′
2
( )( ) ( )
+ x 2 − ν2 x 2 − 1 + 4 a x 2 − µ a x 4 x 2 − 1 y = 0.
2
(10.652a, b)
µ x2
y ( x ) = sec x J νµ ( x ) : x 2 y ′′ + 1 − 2 x tan x + x y′
2
(10.653a, b)
2 µx 2
− ν + x tan x 1 + y = 0.
2
158 Differential Equations and Applications
µ x2
y ( x ) = csc x J νµ ( x ) : x 2 y ′′ + 1 − 2 x cot x + x y′
2
(10.654a, b)
2 µ x2
− ν + x cot x 1 + y = 0.
2
x
∫ j ( ξ ) d ξ = y ( x ) J νµ ( x ) :
−1
y ( x ) = exp
µ x2 µ x3 j
x 2 y ′′ + 1 − 2 x j +
2 ( )
x y ′ + x 2 1 + j 2 − j ′ − x j − ν2 −
2
y = 0.
(10.655a–c)
Classification 10.2
159
160 Differential Equations and Applications
d2 x dx
2
+ 2λ + ω 20 x = f ( t ) . (10.656)
dt dt
m x + µ x + k x = Fm ( x ) ; (10.657)
1
L q + R q + q = Fe ( x ) . (10.658)
c
x + ( λ − ν ) x + ω 20 x = 0. (10.659)
*problem 33. Distinct applied and oscillation (10.662c) frequencies and dis-
tinct applied decay and natural damping (subsection 2.8.4): ω ≠ ω a , λ ≠ λ a ;
1/2
ω = ω 20 − λ 2 . (10.662c)
+∞
nπt
f ∈F ( 0, T ) : x + 2 λ x + ω 20 x = ∑ f cos
n =−∞
n
T
≡ f ( t ) . (10.664a, b)
G + 2 λ G + ω 20 G = δ ( t ) . (10.666)
b2
x + 2 λ x + ω 20 x = f a exp − 2 (10.667)
T
*problem 69. Undamped forcing by Gaussian random noise with mean µ and
variance σ 2 (note 2.13):
( t − µ )
2
fa
x + 2 λ x + ω 20 x = exp − . (10.668)
σ 2π 2 σ 2
dy c x + d y
ad ≠ bc: = . (10.669a, b)
dx a x + b y
*problem 70. Contact point of inflexion of paths tangent to the principal line
(subsection 4.1.4) for:
∆ ≡ ( a − d ) + 4 c b = 0 ; (10.670)
2
*problem 71. Saddle point through which pass two principal lines, to which
the paths are asymptotically tangent (subsection 4.1.6): ∆ > 0, b c > a d;
*problem 72. Nodal point through which pass all paths tangent to a princi-
pal line (subsection 4.1.7): ∆ > 0 in (10.670) and:
ad − bc
ϑ=4 ≠ 1, 0. (10.671)
(a + d + ∆ )
2
*problem 73. Focal point through which pass all straight paths (subsec-
tion 4.1.7): ϑ = 1;
*problem 74. Simple point through which passes one of the set of parallel
straight paths (subsection 4.1.7): ϑ = 0;
*problem 75. Common center of elliptical paths (subsection 4.1.8):
∆ < 0, a + d = 0;
*problem 76. Asymptotic point as the limit of spiral paths (subsection 4.1.9):
∆ < 0, a + d ≠ 0.
Classification 10.2 165
B. b. Paths Near a Stagnation Point of the Second Degree (example 10.7):
*problem 77. Paths tangent to the y-axis at the origin (E10.7.1):
dy y 2 − x 2
= ; (10.672)
dx xy
*problem 78. Paths with three asymptotes crossing at the origin with equal
angles (E10.7.2):
dy x 2 − y 2
2 = . (10.673)
dx xy
dx dy
= a x + b y, = c x + d y . (10.674a, b)
dt dt
a+d d2 d
λ=− , ω 20 = a d − b c: 2 − 2λ + ω 20 x , y ( t ) = 0, (10.675a–c)
2 dt dt
with (19a) damping λ > 0 (amplification λ < 0) and (19b) natural frequency
ω 20 > 0 (growth rate ω 20 < 0) in (10.675c);
*problem 80. Critical damping corresponding to (problem 70) to a contact
point of inflection (subsection 4.2.4): λ = ω 0;
*problem 81. Supercritical damping corresponding to (problem 72) to a nodal
point (subsection 4.2.5): λ > ω 0 , 0 < r s ≠ 0:
r , s = − λ ± λ 2 − ω 20 ; (10.676a, b)
d dx
m ( t ) + k ( t ) x ( t ) = 0. (10.678)
dt dt
d2 x
+ ω 20 1 + 2 h cos ( ω e t ) x ( t ) = 0. (10.679)
dt 2
2
ωe t 2 ω0
w (θ) ≡ x (t ) , θ ≡ , a≡ , q ≡ − h a : (10.680a–d)
2 ω e
d2 w
+ a − 2 q cos ( 2θ ) w ( θ ) = 0. (10.680e)
dθ2
C. c. Stability/Instability/Indifference of
Solutions (subsections 4.3.3–4.3.6):
*problem 92. For linear unforced second-order differential equation (10.681e)
with (10.681a, b) periodic coefficients (subsections 4.3.3–4.3.5):
n = 0, 1, 2: An ( t ) = An ( t + τ ) , ∑ n= 0
An ( t )
dx x
dt n
= 0; (10.681a–c)
Classification 10.2 167
*problem 93. For undamped harmonic oscillator (10.682b) with (10.682a) peri-
odic “natural frequency” (subsection 4.3.6):
d2 x
ω (t ) = ω (t + τ ): + ω ( t ) x ( t ) = 0.
2
(10.682a, b)
dt 2
m±1
m, n = 1, 2,...: ω e± =
m
( ω 0 + ε ) , (10.683a–c)
including excitation frequency equal to the natural frequency
( n = 1, m = 2 : ω +
)
e = ω0 .
2
ω h d2 x dx
λ 2 << b 2 ≡ 0 : + 2λ + ω 20 1 + 2 h cos ( ω e t ) x ( t ) = 0. (10.684a, b)
2 dt 2
dt
g d2θ dθ g 8p
λ 2 << , θ2 << 1: + 2λ + 1+ cos ( ω e t ) θ ( t ) = 0. (10.685a–c)
L dt 2 dt L L
2
d2θ dΦ 1 dx
m = F ( x) = − : m + Φ ( x ) = E . (10.686a, b)
dt 2 dx 2 dt
1 2 β
Φ( x) = k x 1 + x 2 , (10.687)
2 2
Classification 10.2 169
1 2 2 1
Φ( x) = k x 1 + ε α x + ε 2β x 2 . (10.688)
2 3 2
d2 x dx
m
dt 2
+µ
dt
( )
+ k x 1 + ε α x + ε 2 β x 2 = Fa cos (β ω a t ) . (10.690)
170 Differential Equations and Applications
*problem 117. Non-linear resonance with applied frequency close to the nat-
ural frequency ω a = ω 0 + ε leading to the existence of amplitude jumps and
hysterisis loop (subsections 4.6.1–4.6.3);
*problem 118. Problem 117 for subharmonic non-linear resonance with
applied frequency close to half of the natural frequency 2 ω a = ω 0 + ε
leading to smaller amplitude jumps and weaker hysteresis (subsections
4.6.4–4.6.5);
*problem 119. Non-linear resonance with applied frequency ω a = 2 ( ω 0 + ε )
close to the first harmonic of the natural frequency leading to smooth ampli-
tude variation, or amplitude jumps or suppressions of oscillations (subsec-
tions 4.6.6–4.6.7);
*problem 120. Non-linear resonance with applied frequency ω e = 3 ( ω 0 + ε )
close to the second harmonic of the natural frequency leading to no oscilla-
tion, or suppression of oscillations or a threshold amplitude for the existence
of oscillations (subsections 4.6.8–4.6.9);
*problem 121. Non-linear resonance at rational multiples of the natural fre-
quency, that is, submultiples of harmonics of the natural frequencym (sub-
section 4.6.10).
dJ ZΩ
L +RJ = ± J . (10.691)
dt 2π
2 π R << Z Ω : L J = − R J ± S J 2 . (10.692a, b)
d2 J 1 d J dJ R M
− ± b J bJ J + = 0. (10.693)
2
dt J d t dt L a L
*problem 126. Without driving torque M = 0 the electric current and angular
velocity decay to zero (subsection 4.7.5);
*problem 127. Balance of driving and Lorentz force torques M Ω = R J 2 leads
to constant electric current and angular velocity (subsection 4.7.6).
dv±
= b v± − a ( v± ) . (10.694)
2
m
dt
*problem 128. For upper sign in (10.694) linear and quadratic damping lead-
ing to decay;
*problem 129. For lower sign in (10.694) linear amplification and quadratic
damping leading to a steady state.
g + ω 02 sin θ = 0.
ω0 = : θ (10.695a, b)
L
µ > 0: + µ θ θ + m g sin θ = 0.
mLθ (10.696a, b)
x = f ( x ; q ) . (10.697)
d2 x dx
+ γ h( x) + ω 20 x = 0. (10.699)
dt 2 dt
*problem 140. Weak γ << 1 damping (subsection 4.8.6);
*problem 141. Strong γ >> 1 damping (subsection 4.8.7).
x3 d2 x ∂Φ
Φ ( x; a) = − a x: 2
=− = a − x 2 ; (10.700a, b)
3 dt ∂x
x 4 a x2 d2 x ∂Φ
Φ ( x; a) = − + − b x: =− = x 2 − a x + b . (10.701a, b)
4 2 dt 2 ∂x
Classification 10.2 173
dx dy
= ax − b x y , = − c y + d x y . (10.702a, b)
dt dt
∂y ∂v
{u, v} ≡
dx dy 1
, ∈C ( D ) : + ≠ 0. (10.703a, b)
dt dt ∂x ∂y
dx dy
u, v ∈C 2 (|R ) : = u ( x, y ) , = v ( x , y ) ; (10.704a–c)
dt dt
dx dv
= v, = −2 λ v − ω 20 x . (10.705a, b)
dt dt
d2 y dy
A2 ( x ) 2
+ A1 ( x ) + A0 ( x ) y ( x ) = B ( x ) . (10.706)
dx dx
d2 y dy
2
+ P ( x) + Q ( x ) y = f ( x ) . (10.707)
dx dx
*problem 162. Idem (10.708a) the general integral (10.709c) of the unforced
(10.709a, b) equation (notes 5.4–5.5):
*problem 163. Idem (10.710a) the complete integral (10.710c) of the forced
(10.710b) equation (notes 5.4–5.5):
du x
dx
+ u2 + P u + Q = F ( x ) exp −
∫ u ( ξ ) d ξ ; (10.711)
du dy
R ( x) + Q ( x ) y = T ( x ) ; (10.712)
dx dx
d2 v
+ I ( x ) v ( x ) = G ( x ) ; (10.713)
dx 2
d2 Ψ
+ α 2 ( x − xr ) Ψ ( x ) = 0. (10.714)
dx 2
∂2 Φ ∂2 Φ ∂2 Φ 1 ∂2 Φ
+ + = . (10.715)
c ( x ) ∂t
2
∂x2 ∂y 2 ∂z2 2
d 2 Ψ ω 2
2
+ 2 − k Ψ ( x ) = 0. (10.716)
c ( x )
2
dx
176 Differential Equations and Applications
d 4ζ d 2ζ
EI − T = 0. (10.717)
dx 4 dx 2
*problem 184. Clamped at both ends (subsection 6.1.5);
*problem 185. Pinned at both ends (subsection 6.1.6);
*problem 186. Clamped at one end, pinned at the other end (subsection 6.1.7);
*problem 187. Non-linear and linear boundary conditions at a free end (sub-
section 6.1.8);
*problem 188. Cantilever beam, which is clamped at one end and free at the
other end (subsection 6.1.9);
*problem 189. Comparison of critical buckling loads for four diferent combi-
nations of supports (subsection 6.1.9).
′ ′
dζζ′ 2 3 2 ′ ζ′ 2
ζ′ ≡
1 − ; ζ′ 3 << 1:
ζ ′′ 1 − ζ ′ + p 2
ζ ′ 1 − = 0.
dx 2 2
2
(10.718a–c)
Classification 10.2 177
0 ≤ a ≤ L: k ζ ( a ) = E I ζ′′′ ( a ) − T ζ′ ( a ) ; (10.719a, b)
*problem 195. Cantilever beam with a point translational spring at the free
end (subsection 6.2.2);
*problem 196. Clamped beam with a point translational spring in the middle
(subsection 6.2.3);
*problem 197. Pinned beam with a point translational spring in the middle
(subsection 6.2.4);
*problem 198. Clamped-pinned beam with a point translational spring in the
middle (subsection 6.2.5).
0 ≤ a ≤ L: − E I ζ′′ ( a ) = k ζ′ ( a ) ; (10.720a, b)
*problem 200. Cantilever beam with a point rotary spring at the free end
(subsection 6.2.7);
*problem 201. Clamped beam with a point rotary spring in the middle (sub-
section 6.2.8);
*problem 202. Pinned beam with a point rotary spring in the middle (subsec-
tion 6.2.9);
*problem 203. Clamped-pinned beam with a point rotary in the middle (sub-
section 6.2.10).
*problem 206. Buckling loads of all orders for a cantilever beam with a point
rotary spring at the free end (subsections 6.2.14–6.2.15);
*problem 207. Problem 206 with a point translational spring (subsections
6.2.14 and 6.2.16);
*problem 208. Buckled shape of a cantilever beam with a translational or
rotary spring at the free end (subsections 6.2.14–6.2.17).
d 4ζ d 2ζ
EI 4
− T 2 + k ′ζ = 0. (10.721)
dx dx
*problem 210. Pinned beam under traction T > 0 without k ′ = 0, spring sup-
port (subsection 6.3.3);
*problem 211. Pinned bar with attractive springs k ′ > 0 and T = 0, no axial
tension (subsection 6.3.4);
*problem 212. Pinned bar supported on repulsive springs k ′ < 0 without
T = 0, axial tension (subsection 6.3.5);
*problem 213. Pinned beam with balance T 2 = 4 E I k ′ > 0 of axial compression
T < 0 and attractive spring k ′ > 0 (subsection 6.3.6);
*problem 214. Pinned beam with balance T 2 = 4E I k ′ > 0 of axial traction T > 0
and attractive spring k ′ > 0 (subsection 6.3.7);
*problem 215. Pinned beam with predominance T 2 > 4 E I k ′ of axial traction
T > 0 and k ′ > 0 attractive springs (subsection 6.3.8);
*problem 216. Pinned beam with predominance T 2 > 4 E I k ′ of axial compres-
sion T < 0 and k ′ > 0 attractive springs (subsection 6.3.9);
*problem 217. Pinned beam with predominance T 2 > 4 E I k ′ of axial tension
over k ′ < 0 repulsive springs (subsection 6.3.10);
*problem 218. Pinned beam with predominance 0 < T 2 < 4 E I k ′ of axial atrac-
tive springs k ′ > 0 over axial tension (subsection 6.3.11);
*problem 219. Extension of the problems 185 and 210–218 from pinned to can-
tilever, clamped and clamped-pinned beams using an effective axial tension
(subsection 6.3.12);
*problem 220. To all problems 185 and 210–219, apply a combined tension-
spring buckling load that does not affect the shape of the buckled elastica
(subsection 6.3.13);
*problem 221. Extension of problem 219 to (10.722) to bending under own
weight (subsection 6.3.14):
d 4ζ d 4ζ
EI 4
− T 4 + k ′ ζ = ρ g . (10.722)
dx dx
Classification 10.2 179
d 4ζ d 2ζ
T > 0: EI 4
−T = f ( x ) . (10.723a, b)
dx dx 2
L
ζ ( 0 ) = 0 = ζ′ ( 0 ) , ζ ( L ) = 0 = ζ′ ( L ) : f ( x ) = Q δ ′ x − ; (10.724a–e)
2
L
ζ ( 0 ) = 0 = ζ′′ ( 0 ) , ζ ( L ) = 0 = ζ′′ ( L ) : f ( x) = F δ′ x − ; (10.725a–e)
2
ζ ( 0 ) = 0 = ζ′ ( 0 ) , ζ ( L ) = 0 = ζ′′ ( L ) : f ( x ) = ρ g ; (10.726a–e)
x2
F = const: E ( x ) = E0 + ( E1 − E0 ) . (10.729a, b)
L2
∂4 Θ ∂4 Θ ∂4 Θ
( )
Θ∈D 4 |R 2 : 0=
∂x 4
+
∂y 4
+ 2 2 2 ;
∂x ∂y
(10.731a, b)
*problem 236. Strains and stresses due to the drilling of a circular hole in an
infinite thin elastic plate (subsection 6.5.6);
*problem 237. Stress concentration near a circular hole in an infinite plate
with normal stresses at infinity (subsections 6.5.7–6.5.8);
*problem 238. In-plane deformations, strains, and stresses in a circular disk
subject to a uniform pressure along the periphery (subsection 6.5.9).
∂ ∂ζ
Tαβ Nβ δζα = 0: 0=
∂ xα Tαβ ∂ x + f ( xγ ) ; (10.732a, b)
β
∂2 ∂2 ∂2
T11 + T22 + 2 T12 G ( x1 , x2 ) = δ ( x1 ) δ ( x2 ) ; (10.733)
∂ x1 ∂ x1 ∂ x2 ∂ x2 ∂ x1 ∂ x2
∂2 ∂2 ∂2
T11 2 + T22 2 + 2 T12 ζ ( x1 , x2 ) = − f ( x1 , x2 ) . (10.734)
∂ x1 ∂ x2 ∂ x1 ∂ x2
*problem 251. Turning moments and linear weak deflection of an isotropic elas-
tic plane (10.735b) with constant (10.735a) bending stiffness (subsection 6.7.6):
Eh3 f ( x, y ) ∂4 ζ ∂4 ζ ∂4 ζ
D≡ = const: = + + 2 . (10.735a, b)
( )
12 1 − σ 2 D ∂x 4 ∂y 4 ∂x2 ∂y 2
O. d. W
eak Linear Bending of Isotropic Circular
Plates (subsections 6.7.15–6.7.19):
*problem 266. Concentrated force at the center and clamped at the boundary
(subsections 6.7.15–6.7.16);
Classification 10.2 183
( )
f = ∇ 2 D ∇ 2ζ − h ∂α (Tαβ ∂β ζ ) . (10.736)
*problem 274. Comparison of elastic stressed plate and beam under tension
(subsection 6.8.1);
*problem 275. Balance equation for the transverse displacement of a uniform
stressed plate and beam (subsection 6.8.1);
*problem 276. Extreme cases of transverse deflection of a membrane and
weak bending of an elastic plate (subsection 6.8.2).
dvi
fi + ∂ j Tij = ρ ai = ρ ; (10.737a, b)
dt
*problem 278. Newton law (1686) in the absence of stresses for mass inde-
pendent of time (subsection 6.8.3);
*problem 279. Balance of moments of forces and symmetry of the stress ten-
sor (subsection 6.8.4);
184 Differential Equations and Applications
*problem 280. Work of the volume and surface forces (subsection 6.8.5);
*problem 281. Internal energy, including the work of deformation of the
stresses on the strains (subsection 6.8.6).
D≡
h3
C −
( C13 )2 = h3 C − ( C23 )2 . (10.738)
11 22
12 C33 12 C33
dEˆ
E d = d = − N x δζ − Mn δ ( ∂n ζ ) ; (10.739)
ds
*problem 301. Vanishing of the total elastic energy along the boundary of a (i)
isotropic or (ii) pseudo-isotropic orthotropic plate (E10.15.1–E10.15.2);
*problem 302. Stress couple and turning moment for a pseudo-isotropic
orthotropic plate of arbitrary shape (E10.15.3);
*problem 303. simplification of the stress couple and turning moment for a
pinned plate of arbitrary shape (E10.15.3–E10.15.4);
*problem 304. Problem 303 for a clamped plate (E10.15.3–E10.15.4);
*problem 305. Comparison of (a) stress couple and (b) turning moment for
(i) clamped or (ii) pinned plate of an (α) isotropic and (β) pseudo-isotropic
orthotropic material (E10.15.5);
*problem 306. Boundary conditions for elastic plates: (i) of arbitrary, circular,
or rectangular shape; (ii) with clamped, pinned, or free boundary; (iii) made of
isotropic or pseudo-isotropic orthotropic elastic material (E10.15.6–E10.15.9).
( ) ( ) ( )
D∇ 2ζ − f − h ∂ yy Θ ( ∂ xx ζ ) + ( ∂ xx Θ ) ∂ yy ζ − 2 ∂ yy Θ ∂ yy ζ , (10.740a)
( ) ( )( )
∇ 4 Θ = E ∂ yy ζ − ∂ xy ζ ∂ yy ζ ; (10.740b)
2
∇ 4ζ = − f , ∇ 4Θ = 0 ; (10.741a, b)
*problem 315. Exact strain tensor due to transverse and in-plane displace-
ments of a plate (subsection 6.9.2);
*problem 316. Using the Hooke law and stress function to eliminate the in-
plane displacements and keep the transverse displacement in (10.740b) (sub-
section 6.9.3).
*problem 322. Boundary condition for the in-plane stress vector and dis-
placement (subsection 6.9.4);
*problem 323. Principle of virtual work, including elastic energies of bend-
ing, deflection, and deformation (subsection 6.9.5);
*problem 324. Surface integral in problem 323 leading to (10.740a, b) the bal-
ance equations (subsections (subsection 6.9.5);
*problem 325. Boundary integral in problem 323 leading to the boundary
conditions (subsection 6.9.5);
*problem 326. Second boundary condition from problem 325 involving the
augmented turning moment (subsection 6.9.5);
*problem 327. Conditions for clamped, pinned, and free boundaries in the
strong bending of an elastic plate (subsection 6.9.5).
{ }
•
∂ζ • ∂ζ • 2 −1/2 ′
ζ′ ≡ ,ζ≡ : 1 − T ζ′ 1 + ζ′
ρ ζ = f a ( x , t ) ; (10.742a–c)
∂x ∂t
*problem 340. Non-linear wave equation (10.743) for the transverse displace-
ment ζ of an elastic membrane with mass density ρ1 and transverse applied
force f a per unit area and isotropic tangential tension T, which are all func-
tions of position and time (note 6.4):
{ } = f ( x, y, t). (10.743)
•
•
ρ2 ζ − ∇. T 1 + ( ∇ζ.∇ζ )
−1/2
a
*problem 341. Linear torsional (10.744a) wave equation (10.744b) for the angle
of rotation φ along a straight elastic rod with moment of inertia I and tor-
sional stiffness of the cross-section C and applied axial torque or moment
Ma , which are all functions of position and time (note 6.5):
•
I φ• − C φ′ ′ = M x , t ;
( )
•
τ ≡ φ: a( ) (10.744a, b)
*problem 342. Linear wave equation (10.745) for the longitudinal displace-
ment u of a straight elastic rod with mass density per unit length ρ1, Young
modulus E and applied longitudinal force f a , all depending on position and
time (note 6.6).
•
ρ u• − E u′ ′ = F x , t ; (10.745)
1 ( ) a( )
*problem 343. Non-linear wave equation (10.746) for the transverse displace-
ment ζ of an elastic beam with mass density per unit length ρ1, bending
stiffness B, tangential tension T, and applied shear stress f a , which are all
functions of position and time (note 6.7):
*problem 344. Linear wave equation (10.747) for the transverse displacement
ζ of an elastic stressed plate with mass density ρ2 and transverse force f a per
unit area in-plane stresses Tαβ and bending stiffness D, which are all func-
tions of position and time (note 6.8):
•
•
2 2
( )
ρ2 ζ − ∂α (Tαβ ∂β ζ ) + ∇ D ∇ ζ = f a ( x1 , x2 , t ) . (10.747)
*problem 348. Particular case of problem 346 for the classical one-dimensional
bending wave equation (10.750) for the wave variable ζ with damping χ,
translational ω t , and rotational ω r natural frequencies, elastic wave speed,
ce and dispersion stiffness parameter b (note 6.10):
•• •
ζ+ 2 χζ+ ω t2 ζ + ω 2r ζ′ − ce2 ζ′′ + b ζ′′′′ = ρ1−1 f a ( x , t ) ; (10.750)
*problem 349. Particular case of problem 347 for the classic two-dimensional
bi-dimensional bending wave equation (10.751) for the wave variable ζ, as in
the one-dimensional case (10.750) with two wave speeds c1,2 in the principal
stress directions, a vector ω r1,2 natural frequency for the rotary spring and a
dispersion stiffness parameter b (note 6.10):
•• •
ζ+ 2 χζ+ ω t2 ζ + ω 2r 1 ∂1 ζ + ω 2r 2 ∂2 ζ − c12 ∂11 ζ
(10.751)
( ∂11 + ∂22 ) fa ( x, y , t ) .
2 2 2 −1
− c ∂22 ζ + b
2 ζ=ρ 1
190 Differential Equations and Applications
∂Φ ∂Φ ∂2 Φ ∂2 Φ ∂2 Φ ∂N + M Φ
0 = F x , t ; Φ, , , 2 , 2 , 2
,..., N M ; (10.753)
∂ x ∂t ∂ x ∂t ∂ x ∂t ∂ x ∂t
N M
∂N + M Φ
∑∑A
n= 0 m= 0
nm ( x, t )
∂ x N ∂t M
= B ( x , t ) ; (10.754)
*problem 371. Wave equation (10.756d) for the torsion τ, which is the rate of
rotation along the axis (10.756a), in the linear-free torsional oscillations of a
straight elastic rod with moment of inertia I and torsional stiffness C of the
cross-section that (10.756b, c) do not depend on time (note 7.10):
∂φ ∂I ∂C ∂2 τ 1 ∂ ∂τ
τ≡ ; =0= : − C ( x ) = 0 ; (10.756a–d)
∂x ∂t ∂t ∂t 2 I ( x ) ∂ x ∂x
∂σ ∂E ∂2 u 1 ∂ ∂u
=0= : − E ( x ) = 0 ; (10.757a–c)
∂t ∂t ∂t 2
σ ( x ) ∂x ∂x
∂b ∂h ∂2 ζ g ∂ ∂ζ
ρ, g = const ; =0= : − h ( x ) b ( x ) = 0 ; (10.758a–e)
∂t ∂t ∂t 2
b ( x ) ∂x ∂x
*problem 374. Wave equation (10.759e) for the electric field E of electro-
magnetic waves in a medium with di-electric permittivity ε and mag-
netic permeability µ independent of time (10.759c, d) in the absence of
electric charges and currents (10.759a, b) involving the speed of light in
vacuo c0 (notes 7.14–7.15):
∂ε ∂µ ∂2 E c 2 ∂ 1 ∂E
q = 0 = J; =0= : − 0 = 0 ; (10.759a–e)
∂t ∂t ∂t 2
ε ( x ) ∂x µ ( x ) ∂x
*problem 375. Wave equation (10.760c) for the acoustic pressure perturbation
p of free-linear (10.760a) sound waves, which is with velocity perturbation
v small compared with the sound speed c s , in a horn, which is a duct with
variable cross-sectional area A, which does not (10.760b) depend on time
(notes 7.16–7.17):
2
v ∂A ∂2 p c2 ∂ ∂p
c << 1, ∂t = 0: − s A ( x ) = 0. (10.760a–c)
s ∂t 2
A ( x ) ∂ x ∂x
Classification 10.2 193
∂2 Φ 2 ∂ Φ 1 ∂Φ
2
− c ( x ) 2 + = 0, (10.761a)
∂t 2
∂ t L ( x ) ∂t
∂2 Φ ∂ 2 ∂Ψ Ψ
− c ( x ) + = 0, (10.761b)
∂ x L ( x )
2
∂t ∂ x
*problem 392. Classical wave equation for sound in a duct of constant cross-
section (note 7.25);
*problem 393. Ray approximation for sound propagation in a horn with a
wavelength short compared with the lengthscale of variations in cross-sec-
tion (note 7.26).
*problem 424. Energy balance equation the work of the external forces minus
the dissipation by friction forces against the total energy, consisting of kinetic
and potential energies (subsections 8.1.5–8.1.6);
*problem 425. Particular case of decoupled oscillators with (i) separate equa-
tions of motion and (ii)(iii) potential (dissipation) functions that involve only
squares of position (velocity) and no cross-products (subsections 8.1.5–8.1.6);
*problem 426. Generalized equation of motion with explicit accelerations xr ,
involving damping λ rs and oscillations ω 2rs matrices, and reduced external
forces fr ( t ) (subsection 8.2.1):
N
r = 1,..., N :
d 2 xr
dt 2
+ ∑ 2 λ
s=1
rs
dxs
dt
+ ω 2rs xs = fr ( t ) .
(10.762a, b)
Fr = 0: ( 2
0 = Det δ rs ξ + 2 λ rs ξ + ω 2
rs ) = ∏ ( ξ − ξ ) ; (10.763a–c)
n
n= 1
d 2 q d 2 q
= 1,..., N : 2
+ 2 λ 2 + ω 2 q = g ( t ) ; (10.764a, b)
dt dt
*problem 438. Solution of the equations of motion (10.762a, b) for the linear
forced damped multidimensional oscillator via (10.764a, b) modal coordi-
nates (subsections 8.2.11–8.2.12);
*problem 439. Modal matrix for the sinusoidal forcing of the (10.762a, b) lin-
ear damped multidimensional oscillator (subsection 8.2.13);
*problem 440. Sinusoidal forcing of (10.762a, b) the undamped λ rs = 0 linear
multidimensional oscillator (subsection 8.2.14);
*problem 441. Sinusoidal forcing, beats, and resonance of an undamped
linear multidimensional oscillator, in terms of modal coordinates (subsec-
tion 8.2.15);
*problem 442. Sinusoidal forcing of a linear multidimensional damped oscil-
lator, in terms of modal coordinates (subsection 8.2.16);
*problem 443. Diagonalization to a sum of squares of the potential and
kinetic energies and dissipation function of a linear multidimensional
oscillator (10.762a, b), using modal coordinates (subsection 8.2.16);
*problem 444. Three (four) cases of free (forced) oscillations of (10.764a, b) the
linear multidimensional oscillator (subsection 8.2.16).
x 1 λ 11 λ 12 x 1 ω 12 ω 12 2 x1 f1 ( t )
+ 2 + 2 = . (10.765)
x2 λ 21 λ 22 x 2 ω 21 ω 22
2
x2 f2 ( t )
*problem 446. Modal frequencies for the undamped free oscillations of two
coupled electrical circuits analogous to the two-dimensional mechanical
oscillator in problem 445 (subsection 8.3.2);
*problem 447. Coupled vertical translational and rotational oscillations of a
two-wheeled vehicle with distinct springs and dampers and choice of the
moment of inertia and position of the center of mass (subsection 8.3.3);
*problem 448. Analogies among the variables, parameters, and forcing of the
coupled electrical (mechanical) circuits in the problem(s) 446 (445 and 447),
which satisfy a similar system of coupled linear differential equations (109)
with constant coefficients (subsection 8.3.4);
*problem 449. Analogous modal matrices, consisting of mass, damping and
resilience matrices, and analogous forcing vectors for problem(s) 446 (445
and 447) of electrical (mechanical) linear two-dimensional oscillators (sub-
section 8.3.5).
*problem 459. Explicit natural frequencies and dampings for the free decou-
pled translational and rotational oscillations of a homogeneous rod, as a par-
ticular case of problem 458 (subsection 8.4.10);
*problem 460. Four exact relations satisfied by the modal frequencies and
dampings of problem 458, for strong damping and coupling of translational
and rotational motions (subsection 8.4.11);
*problem 461. Weak modal dampings and unchanged modal frequencies for
problem 458 ≡ 460, in the case of weak damping (subsection 8.4.12).
*problem 472. Modal dampings of the vibration absorber in the case of weak
damping (subsection 8.6.4);
*problem 473. Operation of the vibration absorber, including application of
the initial conditions (subsections 8.6.5–8.6.6).
203
204 Differential Equations and Applications
1815 Cauchy, A. A. “Mémoire sur les Fonctions qui ne Peuvent Obtenir que Deux
Valeurs.” Journal de l'École Polytechnique 17, 29–107 (Oeuvres, série 2, 1, 91–169,
Gaultier-Villars, 1905, Paris).
1818 Fourier, J. B. Theorie Analytique de la Chaleur. Reprinted Dover, New York, 1953.
1824 Bessel, F. W. “Untersuchung der Theils der Planetarischen Storungen
Welchear aus des Bewegung der Sonne Entsteht.” Berliner Abhandlungen.
1827 Cauchy, A. L. Exercices Mathematiques 2, 159 (Oeuvres 7, 198).
1836 Kummer, E. E. “Ueber die Hypergeometrische Reihe.” Journal für die reine und
angewandte Mathematik 83, 127–172.
1838 Airy, G. B. “On the Intensity of Light Near a Caustic.” Transactions of the
Cambridge Philosophical Society 6, 379–402.
1850 Malmsten, C. J. “Théoremes sur l’Équation Differentielle.” Cambridge and
Dublin Mathematical Journal 5, 180–182.
1853 Riemann, B. “Üeber die Darstellbarkeit einer Function durch eine
Trignometrische Reihe.” Gescimmelte Mathematischan Werke 227–271, reprinted
Dover, New York.
1856 Weiertrass, K. T. W. “Theorie des Abel’schen Functionen.” Journal für
Mathematik 52, 285–379.
1859 Chebychev, P. L. “Sur les Questions de Minima, qui se Ratachent à la
Representation Approximative des Fonctions.” Mémoires de l'Académie
Scientique de St. Petersbourg 7, 199–291 (Oeuvres 1, 271–378).
1859 Jacobi, C. G. J. “Untersuchungen über die Differentialglechung der
Hypergeometrischen Reihe.” Journal für die reine und angewandte Mathematik
56, 149–165.
1863 Maxwell, J. C. Treatise of Electricity and Magnetism. Clarendan Press, Oxford.
1864 Hermite, C. “Sur un Nouveau Development en Série de Fonctions.” Comptes
Rendus de l’Académie des Sciences 58, 93–100, 266, 273 (Oeuvres 2, 293–308,
Gauthiers Villars, 1908).
1864 Lipshitz, R. “De Explicatione per Series Trignometricas Instituenda
Functionum Unices Variablis Arbitrarium, et Praecipue Earum, qua per
Variablis Spatium Finitum Valorum Maximum et Minimum Numerarum
Habent Infinitum Disquisition.” Journal für die reine und angewandte Mathematik
63, 296–308.
1866 Fuchs, J. L. “Zur Theorie der Linearen Differentialgleichungen mit
Verändlichen Koefflizienten.” Journal für die reine und angewandte Mathematik
66, 121–160; 68, 354–385.
1867 Newmann, C. G. Theorie der Bessel’schen Funktionen, Teubner, Leipzig.
1868 Lommel, E. C. J. “Zur Theorie der Bessel’schen Functionen.” Mathematische
Annalen 14, 510–536.
1869 Hankel, H. “Die Cylinderfunktionen Erster and Zweiter Art.” Mathematische
Annalen 1, 467–501.
1873 Frobenius, G. “Ueber die Integration des Linearen Differentialgleichungen
durch Reihen.” Journal für die reine und angewandte Mathematik 76, 214–235.
1873 Mathieu, E. Cours de Physique Mathematique. Gautiers-Villars, Paris.
∞
1879 Laguerre, E. “Sur l’integrale ∫
x
ξ −1 e −ξ d ξ.” Bulletin de la Societé Mathematique
Acoustic rays, see Optics and acoustic Asymptotic stability of damped non-
rays linear oscillator, 82–84
Acoustics of horns, 193–195 negative-definite stability function,
exact elementary solutions, 194 83–84
Gaussian horn and displacement positive-definite stability function,
amplifier, 194 83
power law horns and cylindrical/ Axial traction, 179–180
spherical waves, 194–195
properties, 193–194 Balance equations for strong bending of
reflection and transmission at isotropic elastic plate, 186
interfaces, 195 Balance of forcing, 10–11
Adjoint kernel operator, 141 Banded diagonal system, 126
Airy differential equation, 175 Bars, 159, 182
Airy functions of inverse variable, 91–92 absence of traction, 47
Amplification, fast, 10–11 elástica of, 42–43
Amplified oscillator, 161 equation, 39
Amplitude, 168 pinned, 178
decays, 55 straight, 179–180
displacement is a quadratic function, Beams
13 bending under traction, see Bending
excitation, 167–168 of beam under traction
excitation frequency and, 166–167 cantilever with linearly increasing
Analogies among six types of waves, shear stress, 45–48
193 clamped with concentrated torque,
Analytic coefficients, 127–130; see also 37–40
Linear differential equations elastica under traction, 36–37
with analytic coefficients heavy clamped-pinned, 43–45
Analytic integrals, 90–91 linear bending under traction, 35–36,
Arbitrary order, 126 36
Associated Laguerre differential pinned with concentrated force,
equation, 147–148 40–43
Associated Legendre differential Beats, 199
equation, 145 Bending of beam under traction, 35–48
Asymptotic integrals cantilever beam with a linearly
analytic and high-order poles, 90–91 increasing shear stress, 45–48
for cylindrical Bessel differential clamped beam with concentrated
equation, 136 torque, 37–40
for generalized Bessel differential elastica of a beam under traction,
equation, 136–137 36–37
for spherical Bessel differential heavy clamped-pinned beam, 43–45
equation, 136 linear bending, 35–36, 36, 179
triple pole and, 88–90 pinned beam with concentrated
Asymptotic point, 165, 166 force, 40–43
209
210 Index