The calculation of accurate arbitrary-order sensitivities of eigenvalues and eigenvectors is crucial for structural analysis applications, including topology optimization, system identification, finite element model updating, damage detection, and fault diagnosis. Current approaches to obtaining sensitivities for eigenvalues and eigenvectors lack generality, are complicated to implement, prone to numerical errors, and are computationally expensive. In this work, a novel methodology is introduced that uses hypercomplex automatic differentiation (HYPAD) and semi-analytical expressions to obtain arbitrary-order sensitivities for eigenfrequency problems. The new methodology exhibits no sign of truncation nor subtractive cancellation errors regardless of the order of the sensitivity, it is general, and can obtain any high-order sensitivities with the simplicity of first-order computations. A numerical example is presented to verify the accuracy of the method, where the free vibration of a homogeneous cantilever beam is studied. For this problem, up to third-order sensitivities of the eigenvalues and eigenvectors with respect to the material and geometrical parameters were obtained, considering the cases of close and distinct eigenvalues. The results were verified using analytical equations, showing excellent agreement for the eigenvalues and the eigenvectors. The new method promises to facilitate the computation of sensitivities for eigenfrequency problems into routine practice and commercial software.

A derivative-based uncertainty quantification (UQ) method called HYPAD-UQ that utilizes sensitivities from a computational model was developed to approximate the statistical moments and Sobol' indices of the model output. Hypercomplex automatic differentiation (HYPAD) was used as a means to obtain accurate high-order partial derivatives from computational models such as finite element analyses. These sensitivities are used to construct a surrogate model of the output using a Taylor series expansion and subsequently used to estimate statistical moments (mean, variance, skewness, and kurtosis) and Sobol' indices using algebraic expansions. The uncertainty in a transient linear heat transfer analysis was quantified with HYPAD-UQ using first-order through seventh-order partial derivatives with respect to seven random variables encompassing material properties, geometry, and boundary conditions. Random sampling of the analytical solution and the regression-based stochastic perturbation finite element method were also conducted to compare accuracy and computational cost. The results indicate that HYPAD-UQ has superior accuracy for the same computational effort compared to the regression-based stochastic perturbation finite element method. Sensitivities calculated with HYPAD can allow higher-order Taylor series expansions to be an effective and practical UQ method.