Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

Computational Technology Reviews
ISSN 20448430 Computational Technology Reviews
Volume 2, 2010 Historic, Recent and Ongoing Applications of the WittrickWilliams Algorithm
F.W. Williams and D. Kennedy
Cardiff School of Engineering, Cardiff University, United Kingdom F.W. Williams, D. Kennedy, "Historic, Recent and Ongoing Applications of the WittrickWilliams Algorithm", Computational Technology Reviews, vol. 2, pp. 223246, 2010. doi:10.4203/ctr.2.10
Keywords: eigenvalues, vibration, buckling, natural frequencies, buckling load factors, modes, frames, plate assemblies, dynamic stiffness matrix, transcendental eigenproblems, WittrickWilliams algorithm.
Summary
Eigenvalues (i.e. natural frequencies or buckling load factors) of structures are usually found by the finite element method, which involves a linear eigenproblem for which many solution methods exist. This paper covers the much harder transcendental eigenproblem which occurs when finite element discretisation errors are avoided by using the exact member equations yielded by solution of their governing differential equations. It outlines the WittrickWilliams (WW) algorithm, which was developed to obtain the eigenvalues of such transcendental eigenproblems, summarising its development history, its many applications and the development of associated mode finding methods. Key underlying concepts are indicated.
Working at a trial value of the eigenparameter, the WW algorithm computes the value of J, the number of eigenvalues of the structure between zero and the trial value, by summing the 'sign count' property of the transcendental overall stiffness matrix of the structure with J_{m} for all members of the structure, where J_{m} is the number of clamped/clamped member eigenvalues between zero and the trial value. An overview is given of member stiffness matrices, and associated values of J_{m}, for a wide range of member types. These include: BernoulliEuler beams and beamcolumns; Timoshenko beams and beamcolumns; and plates within prismatic plate assemblies which are isotropic or anisotropic (to cover both metal and laminated composite plates) and which may include outofplane shear deformation. Extensions of the algorithm include: cases for which the transcendental overall stiffness matrix is complex and Hermitian instead of being real and symmetric; the use of Lagrangian multipliers; very fast solutions for structures which are repetitive along one, two or three Cartesian directions; rotationally repetitive structures with members along the axis of rotational periodicity; multilevel substructuring; modes computed to the near machine accuracy to which the eigenvalues have always been found; design optimisation of stiffened wing panels using the WW algorithm to obtain sensitivities to design changes; modal density plots; and locating pass bands and stop bands. Some nonstructural areas of application of the WW algorithm are indicated, including using a structural analogy in mathematics both to enable second and fourth order SturmLiouville problems to be solved competitively and also to obtain unique results for effectively infinite level trees of SturmLiouville equations. purchase the fulltext of this paper (price £20)
go to the previous paper 
