Stable methods to solve the impedance matrix for radially inhomogeneous cylindrically anisotropic structures

Andrew N. Norris, Adam J. Nagy, Feruza A. Amirkulova

Research output: Contribution to journalArticlepeer-review

10 Scopus citations


A stable approach for integrating the impedance matrix in cylindrical, radial inhomogeneous structures is developed and studied. A Stroh-like system using the time-harmonic displacement-traction state vector is used to derive the Riccati matrix differential equation involving the impedance matrix. It is found that the resulting equation is stiff and leads to exponential instabilities. A stable scheme for integration is found in which a local expansion is performed by combining the matricant and impedance matrices. This method offers a stable solution for fully anisotropic materials, which was previously unavailable in the literature. Several approximation schemes for integrating the Riccati equation in cylindrical coordinates are considered: exponential, Magnus, Taylor series, Lagrange polynomials, with numerical examples indicating that the exponential scheme performs best. The impedance matrix is compared with solutions involving Buchwald potentials in which the material is limited to piecewise constant transverse isotropy. Lastly a scattering example is considered and compared with the literature.

Original languageEnglish (US)
Pages (from-to)2520-2531
Number of pages12
JournalJournal of Sound and Vibration
Issue number10
StatePublished - 2013

All Science Journal Classification (ASJC) codes

  • Condensed Matter Physics
  • Mechanics of Materials
  • Acoustics and Ultrasonics
  • Mechanical Engineering


Dive into the research topics of 'Stable methods to solve the impedance matrix for radially inhomogeneous cylindrically anisotropic structures'. Together they form a unique fingerprint.

Cite this