Reduction method for the non‐linear analysis of symmetric anisotropic panels

Reduction method and computational procedures are presented for reducing the size of the analysis model and the number of degrees of freedom used in predicting the non-linear response of symmetric anisotropic panels. The two key elements of the method are (a) operator splitting, or decomposition of the characteristic arrays of the finite element model into sums of orthotropic and non-orthotropic contributions, (b) application of a reduction method through the successive use of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate a small number of global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the classical Rayleigh-Ritz technique. The global approximation vectors are selected to be those commonly used in single (or multiple) parameter perturbation techniques, namely a non-linear solution corresponding to zero non-orthotropic arrays and a number of its derivatives with respect to an anisotropic tracing parameter (and possibly, to a load or arc-length parameter in the solution space). The size of the analysis model used in generating the global approximation vectors is identical to that of the corresponding orthotropic structure. The effectiveness of the proposed reduction method is demonstrated by means of a numerical example, and its potential for solving quasi-symmetric non-linear problems of anisotropic panels is discussed.