© Copyright July 1998, January 2021

Web site: https://acwhyte.droppages.com

**Contents:**

1.1 About Kelvin Functions

1.2 Why the Fourth Derivative?

1.3 About Schleicher Functions

2.1 The Power Series

2.2 The First Derivative

2.3 The Second Derivative

2.4 The Third Derivative

2.5 The Fourth Derivative

2.6 The Asymptotic Series

1. INTRODUCTION

**KELGEN-1** is a Kelvin Function Generator program that will calculate and plot Kelvin
and Schleicher functions up to the 4th derivative

The program is public domain. You can download the program and documentation [0.6MB]. The README text file should be read. To run the program double click the file KELGEN (or KELGEN.EXE).

**1.1 ABOUT KELVIN FUNCTIONS.** Kelvin functions were defined
by Lord Kelvin (William Thomson) of Largs in his presidential address
to the Institution of Electrical Engineers in 1889 (ref. 1). Kelvin was
describing the ohmic resistance of a wire with a current density distribution
given by a relationship containing two functions *Ber* and *Bei*.
The origins of the functions stem from the paper by Oliver Heaviside in 1884
(ref. 2). In this paper the solution for oscillatory currents include two
functions defined as M and N which Kelvin subsequently defined as *Ber*
and *Bei* respectively. Kelvin defined these functions by three letters
for convenience in a similar way that we have *Sin*, *Cos* and
*Tan*. The functions stand for: Bessel real and Bessel imaginary.
F. W. Bessel being a German astronomer who was interested in the perturbation
of planets. He is famous for his work in certain types of equations and
functions known as Bessel equations and Bessel functions.

Closely associated with *Ber* and *Bei* are the complementary functions
*Ker* and *Kei* defined by Alexander Russell in his paper of 1909 (ref. 3).
The four functions *Ber*, *Bei*, *Ker*, and *Kei*, collectively
are known as the Kelvin or Thomson functions. The functions are related to the zero
order Bessel functions
*J _{o}*

where

*J _{o}*

*x *is the argument of the function

*i = *Ö-1

and it is customary to omit the zero order subscript when referring to the Kelvin functions.

Although Kelvin was interested in applying the functions to an electrical problem,
McLachlan (ref. 4) shows that there are a range of applications in acoustics,
electrical, fluid and mechanical engineering for which Bessel functions can provide a
solution. My own interest lies in mechanical engineering and in particular to
developing *AJAP1* and reference (6) to include a larger range of plate and shell
applications including: cylinders with variable thickness, shallow spherical shells,
circular flat plates on an elastic foundation, etc.

**1.2 WHY THE FOURTH DERIVATIVE?** Mechanical engineers are generally interested
in obtaining internal loads (forces and moments) and deformations (deflections and
rotations) of a structure (e.g. a cylinder, spherical shell , flat plate etc.) when
subjected to a given set of external loads and boundary conditions. A knowledge of the
internal loads allow the stresses in the structure to be calculated and these can be
assessed for safety, against a set of rules such as a pressure vessel code. In some
cases deformation is the limiting criterion rather than stress. e.g. bolted flanges
have been known to leak because the flange and hub were not stiff enough, even though
they met strength criteria.

As an example, if we consider a cylindrical water tank with variable wall thickness. This problem has been analysed by several investigators (ref. 8, 9 and 10). A solution for the rotation, moment and shear force requires successive differentiation starting from a knowledge of the deflection of the structure. The solution involves the Bessel equation of standard form with a parameter of unity.

h^{2}*d ^{2}*z

The solution of this equation is known and after some mathematics the radial deflection of the cylinder can be shown to be:

*y = [C _{1 }Ber'(*x

Where *Ber'(*x*), Bei'(*x*), Ker'(*x*)* and* Kei'(*x*)*
are the first derivatives of the Kelvin functions and x is
a real argument of the length *x.* i.e. x* = 2*lÖ*x*.
The terms *C** _{1 }*,

Thus it can be seen that the solution for the radial deflection contains the first
derivative of the Kelvin functions. The derivative *dy/dx* gives the rotation of
the cylinder and hence contains the second derivative of the Kelvin functions. The
meridional bending moment is obtained from *M _{x} = -Ddy^{2}/dx^{2}* and
hence contains the third derivative of the Kelvin functions. Finally the radial shear
force

The above approach, with some modifications, can be used to obtain solutions for a range of mechanical engineering problems including: shallow spherical shells, conical shells, circular plates on an elastic foundation, nozzles in pressure vessels and heat exchanger tube sheets. References (8), (9), (10), (11), (13), (14), (15) and (16) provide examples of the solution to these problems.

**1.3 ABOUT SCHLEICHER FUNCTIONS.** Some references e.g. (9), (10), (14) and (15)
make use of the Schleicher functions *Z _{1 }*,

*Z _{1 }(x) = Ber(x)*

*Z _{2 }(x) = -Bei(x)*

*Z _{3 }(x) = -(2/*p

*Z _{4 }(x) = -(2/*p

2. THE KELVIN FUNCTIONS

**2.1 THE POWER SERIES. **There are many text that quote the Kelvin functions
(or in the form of Schleicher functions) e.g. references (4), (5), (7), (9), (10), (13),
(14), (15), (16), (17) and (18). Reference (5) by Watson is comprehensive though very
mathematical. The text by McLachlan (ref. 4) is suitable for engineers and there are
several interesting practical problems discussed. References (17) and (18) are also
worth reading. The text of particular use to engineers interested in the stress
analysis of plates and shells are references (7), (8), (9), (10), (11), (13), (14),
(15) and (16). The notation and presentation used by the various references differ
and this text is no different in this respect, so it is necessary to be careful when
comparing the various text.

The functions are given by the following power series:

*Ber(x) = 1 - x ^{4}/(2^{4 } 2!^{2}) + x^{8}/(2^{8} 4!^{2}) - x^{12}/(2^{12 }6!^{2}) + - . . .*

*Bei(x) = x ^{2}/(2^{2} 1!^{2}) - x^{6}/(2^{6} 3!^{2}) + x^{10}/(2^{10} 5!^{2}) - + . . .*

*Kei(x) = -(*p*/4)Ber(x) - ln(*g*x/2)Bei(x) + x** ^{2}*f

*Ker(x) = (*p*/4)Bei(x) - ln(*g*x/2)Ber(x) - x** ^{4}*f

where *x* is a positive real argument

*ln* is *Loge*

g *= 1.781073* i.e. *ln(*g*) = 0.577216* is Euler's constant.

f*(1) = 1*

f*(2) = 1 + 1/2*

f*(3) = 1 + 1/2 + 1/3*

f*(4) = 1 + 1/2 + 1/3 + 1/4*

f*(5) = 1 + 1/2 + 1/3 + 1/4 + 1/5*

f*(6) = 1 + 1/2 + 1/3 + 1/4 + 1/5 + 1/6*

Thumbnail views of plots of: *Ber, Bei, Kei and Ker*. Click view for a larger plot.

**2.2 THE FIRST DERIVATIVE.** The first derivative can be obtained by a term by term differentiation of the power series in paragraph **2.1**.

*Ber'(x) = - x ^{3}/(2^{3} 1! 2!) + x^{7}/(2^{7} 3! 4!) - x^{11}/(2^{11} 5! 6!) + - . . .*

*Bei'(x) = x/(2 ^{1} 0! 1!) - x^{5}/(2^{5} 2! 3!) + x^{9}/(2^{9} 4! 5!) - + . . .*

*Kei'(x) = -(*p*/4)Ber'(x) - ln(*g*x/2)Bei'(x) - (1/x)Bei(x) *

*+ x*f*(1)/(2 ^{1} 0! 1!) - x^{5}*f

*Ker'(x) = (*p*/4)Bei'(x) - ln(*g*x/2)Ber'(x) - (1/x)Ber(x) *

*- x ^{3}*f

Thumbnail views of plots of: *Ber', Bei', Kei' and Ker'*. Click view for a larger plot.

**2.3 THE SECOND DERIVATIVE.** For the second derivative we can make use of the fact that the the functions satisfy a differential equation of the Bessel type.
Hence we can write:

Ber''(x) = -Bei(x) - (1/x)Ber'(x)

Bei''(x) = Ber(x) - (1/x)Bei'(x)

Kei''(x) = Ker(x) - (1/x)Kei'(x)

Ker''(x) = -Kei(x) - (1/x)Ker'(x)

Thumbnail views of plots of: *Ber'', Bei'', Kei'' and Ker''*. Click view for a larger plot.

**2.4 THE THIRD DERIVATIVE.** The third derivative is obtained by differentiating
the equations in paragraph **2.3**

Ber'''(x) = -Bei'(x) - (1/x)Ber''(x) + (1/x^{2})Ber'(x)

Bei'''(x) = Ber'(x) - (1/x)Bei''(x) + (1/x^{2})Bei'(x)

Kei'''(x) = Ker'(x) - (1/x)Kei''(x) + (1/x^{2})Kei'(x)

Ker'''(x) = -Kei'(x) - (1/x)Ker''(x) + (1/x^{2})Ker'(x)

Thumbnail views of plots of: *Ber''', Bei''', Kei''' and Ker'''*. Click view for a larger plot.

**2.5 THE FOURTH DERIVATIVE.** The fourth derivative is obtained by
differentiating the equations in paragraph **2.4**

Ber^{IV}(x) = -Bei''(x) - (1/x)Ber'''(x) + (2/x^{2})Ber''(x) - (2/x^{3})Ber'(x)

Bei^{IV}(x) = Ber''(x) - (1/x)Bei'''(x) + (2/x^{2})Bei''(x) - (2/x^{3})Bei'(x)

Kei^{IV}(x) = Ker''(x) - (1/x)Kei'''(x) + (2/x^{2})Kei''(x) - (2/x^{3})Kei'(x)

Ker^{IV}(x) = -Kei''(x) - (1/x)Ker'''(x) + (2/x^{2})Ker''(x) - (2/x^{3})Ker'(x)

Thumbnail views of plots of: *Ber ^{IV}, Bei^{IV}, Kei^{IV} and Ker^{IV}*. Click view for a larger plot.

**2.6 THE ASYMPTOTIC SERIES.** The power series discussed in paragraphs **2.1 **to** 2.5** are
applicable for small values of the argument *x*. With large values of the argument,
typically *x>6*, it is appropriate to use the asymptotic series to represent
the functions. There are several forms of presenting the asymptotic series ranging
from the most basic first approximation (references 8, 9, and 10) through to formulas
involving many series terms, e.g. Flügge (ref. 7). Two cases will be considered:
case 1 the most basic approximation and case 2 the more complete series due to Russell
(ref. 3).

**Case 1.** The most basic approximation generally considered sufficiently
accurate for most engineering purposes, e.g. references (8), (9) and (10), is to
use the following first term in the series.

*Ber(x)» {e ^{(x/}*Ö

*Bei(x) » {e ^{(x/}*Ö

*Kei(x) » -{*Ö*(*p/*2x)e ^{(-x/}*Ö

*Ker(x) » {*Ö*(*p/*2x)e ^{(-x/}*Ö

the first derivative becomes:

*Ber'(x)» {e ^{(x/}*Ö

*- {e ^{(x/}*Ö

*Bei'(x)» {e ^{(x/}*Ö

*- {e ^{(x/}*Ö

*Kei'(x)» {*Ö*(*p/*2x)e ^{(-x/}*Ö

*- {e ^{(x/}*Ö

*Ker'(x)» -{*Ö*(*p/*2x)e ^{(-x/}*Ö

*- {e ^{(x/}*Ö

It is customary to omit the third term (the one with *(2x)* in the denominator)
so the first derivative becomes:

*Ber'(x)» {e ^{(x/}*Ö

*Bei'(x)» {e ^{(x/}*Ö

*Kei'(x)» {*Ö*(*p/*2x)e ^{(-x/}*Ö

*Ker'(x)» -{*Ö*(*p/*2x)e ^{(-x/}*Ö

and this can be simplified to:

*Ber'(x)» {e ^{(x/}*Ö

*Bei'(x)» {e ^{(x/}*Ö

*Kei'(x)» {*Ö*(*p/*2x)e ^{(-x/}*Ö

*Ker'(x)» -{*Ö*(*p/*2x)e ^{(-x/}*Ö

**Case 2.** If the additional terms quoted by Russell (ref. 3) are included we
have a more complete asymptotic series as follows:

*Ber(x)= {e ^{(}*a

*Bei(x)= {e ^{(}*a

*Kei(x)= {*Ö*(*p/*2x)e ^{(}*e

*Ker(x)= {*Ö*(*p/*2x)e ^{(}*e

where

a = *x/*Ö*2 + 1/(8*Ö*2x) - 25/(384*Ö*2x ^{3}) -13/(128x^{4}) - . . .*

e = -*x/*Ö*2 - 1/(8*Ö*2x) + 25/(384*Ö*2x ^{3}) -13/(128x^{4}) + . . .*

b = *x/*Ö*2 - *p*/8 - 1/(8*Ö*2x) - 1/(16x ^{2}) -25/(384*Ö

y = -*x/*Ö*2 - *p*/8 + 1/(8*Ö*2x) - 1/(16x ^{2}) +25/(384*Ö

the first derivative becomes:

*Ber'(x)*=* **{*a*' - 1/(2x)}Ber(x) - {*b*'}Bei(x)*

*Bei'(x)*=* **{*b*'}Ber(x) + {*a*' - 1/(2x)}Bei(x)*

*Kei'(x)*=* ** {*y*'}Ker(x) + {*e*' - 1/(2x)}Kei(x)*

*Ker'(x)*=* **{*e*' - 1/(2x)}Ker(x) - {*y*'}Kei(x)*

where

a*'* = *1/*Ö*2 - 1/(8*Ö*2x ^{2}) + 75/(384*Ö

e*'* = -*1/*Ö*2 + 1/(8*Ö*2x ^{2}) - 75/(384*Ö

b*'* = *1/*Ö*2 + 1/(8*Ö*2x ^{2}) + 1/(8x^{3}) + 75/(384*Ö

y*'* = -*1/*Ö*2 - 1/(8*Ö*2x ^{2}) + 1/(8x^{3}) - 75/(384*Ö

The higher derivatives of the asymptotic series are calculated by the relationships
discussed previously in paragraphs **2.3** to **2.5**.

3. TEST RESULTS

In order to test the results the following steps have been taken:

(1) The results of the program KELGEN-1 have been compared against the quoted tables of results in references 8, 9, and 10. This was done by starting with a low number of terms in the power series and then increasing the number of terms until the results from KELGEN-1 and references 8, 9, and 10 were identical (or virtually identical) over the quoted range of argument X from 0 to 6. It was found that 10 terms in the power series were required to achieve this.

(2) A facility was incorporated into the program to plot the percentage difference (Mode Diff) between the results for terms in the power series and results for terms in the power series* minus one*. i.e. a plot of:

%Diff = (Result, for terms in series - Result, for terms in series-1)/Result, for terms in series

The %Diff between 10 and 9 terms in the series is so small that it plots as zero (or virtually zero) over the range of argument X = 0 to 6.

(3) A facility was incorporated into the program to plot the percentage error (Mode Error) between the asymptotic series and the power series. i.e. a plot of:

%Error = (Result, for asymptotic series - Result, for power series)/Result, for asymptotic series

By inspection it was found that Russell's asymptotic series gave a smoother transition from the power to the asymptotic series and the percentage error was much reduced compared to the basic simplified asymptotic series. Hence Russell's asymptotic series (the case 2 series) is recommended and made the program default for large values of the argument.

All these findings can be checked by the user by running the program, altering the defaults to suit and comparing results.

4. REFERENCES

1. Thomson, W. *Ether, Electricity, and Ponderable Matter.* Presidential Address to the Institution of Electrical Engineers*,* 10th January 1889 and addition of 13th May 1890. Mathematical & Physical Papers, Sir William Thomson, Volume III, Cambridge University Press 1890.

2. Heaviside, O. *The Inductance of Currents in Cores.* The Electrician, Volume XII, 3rd May 1884.

3. Russell, A. *The Effective Resistance and Inductance of a Concentric Main, and Methods of Computing the Ber and Bei and Allied Functions. *Philosophical Magazine and Journal of Science. 6 Series, Volume XVII, Number C, April 1909.

4. McLachlan, N.W. *Bessel Functions for Engineers.* 2nd Edition, Oxford Clarendon Press 1955.

5. Watson, G. N. *A Treatise on the Theory of Bessel Functions.* 2nd Edition, Cambridge 1962.

6. Whyte, A. C. *Basic Discontinuity Analysis of Multishell Axisymmetric Junctions.* Self published 1994 including the *BASIC* program *AJAP1* (Axisymmetric Junction Analysis Program Version 1).

7. Flügge, W. *Four Place Tables of Transcendental Functions.* Pergamon Press 1954.

8. Turner, C. E. *Introduction to Plate and Shell Theory.* Longmans 1965.

9. Hetényi, M. *Beams on Elastic Foundation.* The University of Michigan Press 1974.

10. Timoshenko, S.P. and Woinowsky-Krieger, S. *Theory of Plates and Shells.* 2nd Edition, McGraw- Hill 1959.

11. Flügge, W. *Stresses in Shells.* 2nd edition, Springer-Verlag 1973.

12. Schleicher, F. *Kreisplatten auf elastischer Unterlage.* Berlin 1926.

13. Jawad, M. H. and Farr, J. R. *Structural Analysis & Design of Process Equipment.* 2nd Edition, Wiley Interscience 1989.

14. Singh, K. P. and Soler, A. I. *Mechanical Design of Heat Exchangers*. Arcturus Publishers 1984.

15. Murray, N. W. and Stuart, D. G. *Behaviour of large taper hub flanges.* Symposium on Pressure Vessel Research Towards Better Design. I.Mech.E. 1961.

16. Leckie, F. A. and Penny, R. K. *Solutions for the Stresses at Nozzles in Pressure Vessels.* Welding Research Council Bulletin No. 90, 1963.

17. Bowman, F. *Introduction to Bessel Functions*. Dover Publications, 1958.

18. Farrell, O. J. and Ross, B. *Solved Problems: Gamma and Beta Functions, Legendre Polynomials, Bessel Functions*. MacMillan, 1963.

--Top-- |

--Technical-- |

--Home-- |