LAUSR

Abstract Variable Angle Tow (VAT) composite panels, whereby fibre angle distributions vary continuously in-plane, have more scope for tuning structural properties than traditional straight fibre laminated composite materials. Currently, the geometrically non-linear response of such structures is typically modelled using path-following methods implemented in commercial finite element analysis. However the high computational cost given from both the incremental-iterative procedure and fine mesh required for accurate modelling of the variable trajectory of fibre direction, can be excessively time consuming. Driven by these shortcomings, we use a fast multi-modal Koiter asymptotic method for investigating the nonlinear buckling behaviour of VAT cylindrical panels in compression. In doing so, this paper provides new insight into the multi-modal description and nonlinear buckling mode interactions which are responsible for the highly nonlinear behaviour of cylindrical panels in compression with a particular emphasis on the nonlinear components in the mathematical asymptotic description. As such, through an extensive parametric virtual testing programme we construct, for each panel under investigation, an asymptotic solution projecting the equilibrium equations in the subspace of its first 30 buckling modes. Then, at several points of the resulting equilibrium path, the percentage contribution of each of these modes to the solution is measured. On the basis of the percentage obtained, those modes which give the largest contribution to the solution are identified. The main novelty of this work is that these results are used as a priori information for re-projecting the equilibrium equations into a new modal subspace to which only the buckling modes with the largest participation belong. Moreover, for the first time, we observe that in a multi-modal Koiter-inspired description the modal shapes giving the largest contribution to the asymptotic solution have the same degree of symmetry. In this way, for a generic panel under consideration, we show that it is possible to obtain a more computationally-efficient asymptotic description than current approaches. Finally, the new computed equilibrium paths are compared to benchmark results using the commercial finite element software ABAQUS and the computational advantages given by using a priori information in the asymptotic expansion are commented upon.