Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 99
PROCEEDINGS OF THE ELEVENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping
Paper 100

Solving Homogeneous Trees of Sturm-Liouville Equations using an Infinite Order Determinant Method

A. Watson1, W.P. Howson2, C. James1 and C. Williams1

1Department of Aeronautical and Automotive Engineering, Loughborough University, United Kingdom
2Cardiff School of Engineering, Cardiff University, United Kingdom

Full Bibliographic Reference for this paper
A. Watson, W.P. Howson, C. James, C. Williams, "Solving Homogeneous Trees of Sturm-Liouville Equations using an Infinite Order Determinant Method", in B.H.V. Topping, (Editor), "Proceedings of the Eleventh International Conference on Computational Structures Technology", Civil-Comp Press, Stirlingshire, UK, Paper 100, 2012. doi:10.4203/ccp.99.100
Keywords: exact dynamic stiffness matrix, Wittrick-Williams algorithm, infinite order determinant, graph theory, Sturm-Liouville operator.

Summary
Quantum graph problems occur in many disciplines of science and engineering and they can be solved by viewing the problem as a structural engineering one. The Sturm-Liouville operator acting on a tree is an example of a quantum graph and the structural engineering analogy is the axial vibration of an assembly of bars with identical topology. Previous work has demonstrated that a numerical procedure can be used to provide insight into the spectral behaviour of the free Laplacian in one dimension acting on a metric tree. In the case of an infinitely long tree, the eigenvalues form bands of continuous spectra with an eigenvalue of infinite multiplicity in the middle of the gaps. In this paper, as in previous work, only finite length trees are considered.

The required eigenvalues are found using the well-established Wittrick-Williams algorithm in conjunction with a dynamic stiffness approach. It is known that the plot of the determinant of the stiffness matrix versus a trial eigenparameter has several poles and that the eigenvalues of the system can occur when the determinant of the stiffness matrix is either zero or infinite. For a finite tree there are no bands of continuous spectra, however all eigenvalues fall within the lower and upper bounds of each band and the bands are repetitive.

The determinant plot starts from a positive value at zero on the abscissa and after crossing the abscissa one or more times tends to infinity as the trial eigenparameter tends to the clamped/clamped member eigenvalue. An established procedure transforms this plot into what is known as the 'infinite order determinant', in which there are no poles and all the original eigenvalues now correspond to the determinant being zero. These properties are used in this paper to show that although the transformed plot tends to zero for large values of trial eigenparameter, a small modification can produce a plot that repeats at regular intervals and corresponds to the band gap structure of the spectrum for an infinitely long tree. From the results of finite length trees it is then possible to estimate the edges of the spectrum for the infinite system. Furthermore, it is known that the eigenvalues of finite trees are obtainable from a family of stepped bars. This offers the possibility of significantly reducing the order of the problem to be solved and leads to a very efficient way of estimating the spectral limits.

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £65 +P&P)