We consider and analyze applying a spectral inverse iteration algorithm and its subspace iteration variant for computing eigenpairs of an elliptic operator with random coefficients. With these iterative algorithms the solution is sought from a finite dimensional space formed as the tensor product of the approximation space for the underlying stochastic function space, and the approximation space for the underlying spatial function space. Sparse polynomial approximation is employed to obtain the first one, while classical finite elements are employed to obtain the latter. An error analysis is presented for the asymptotic convergence of the spectral inverse iteration to the smallest eigenvalue and the associated eigenvector of the problem. A series of detailed numerical experiments supports the conclusions of this analysis.
1 January 2019 in Computer Methods in Applied Mechanics and Engineering
In this paper, we give an overview of the issues one must consider when designing methods for vibration based health monitoring systems for perforated thin shells especially in relation to frequency response analysis. In particular, we allow either the material parameters or the structure or both to be random. The numerical experiments are computed using the standard high order finite element method with stochastic collocation for the cases with random material and Monte Carlo for those with damaged or random structures. The results display a wide range of responses over the experimental configurations. In perforated shell structures, the internal boundary layers can play an important role especially when damage is allowed within the penetration patterns. The computational methodology advocated here can be used to build statistical databases that are necessary for development of probabilistic damage identification methods.
Numerical simulation of thin solids remains one of the challenges in computational mechanics. The 3D elasticity problems of shells of revolution are dimensionally reduced in different ways depending on the symmetries of the configurations resulting in corresponding 2D models. In this paper, we solve the multiparametric free vibration of complex shell configurations under uncertainty using stochastic collocation with the p‐version of finite element method and apply the collocation approach to frequency response analysis. In numerical examples, the sources of uncertainty are related to material parameters and geometry representing manufacturing imperfections. All stochastic collocation results have been verified with Monte Carlo methods.
Frequency response analysis under uncertainty is computationally expensive. Low-rank approximation techniques can significantly reduce the solution times. Thin perforated cylinders, as with all shells, have specific features affecting the approximation error. There exists a rich thickness-dependent boundary layer structure, leading to local features becoming dominant as the thickness tends to zero. Related to boundary layers, there is also a connection between eigenmodes and the perforation patterns. The Krylov subspace approach for proportionally damped systems with uncertain Young’s modulus is compared with the full system, and via numerical experiments, it is shown that the relative accuracy of the low-rank approximation of perforated shells measured in energy depends on the dimensionless thickness. In the context of frequency response analysis, it then becomes possible that, at some critical thicknesses, the most energetic response within the observed frequency range is not identified correctly. The reference structure used in the experiments is a trommel screen with a non-regular perforation pattern with two different perforation zones. The low-rank approximation scheme is shown to be feasible in computational asymptotic analysis of trommel designs when the proportional damping model is used.