Search

• Efficient solution procedure for non-linear eigen-value problem based on homotopy and quasi-Newton methods. • Numerical verification of the proposed formulation and testing of various solution variants. • Extensive analyses of results sensitivity to the model parameters. In this paper the numerical procedure to solve the non-linear eigen-value problem describing the free vibrations of the laminate plates with viscoelastic (VE) layers is presented. The derivation of that main equation based on the refined zig-zag theory and using fractional Zener material model to describe VE material is presented in the companion paper Lewandowski R, Litewka P. Wielentejczyk P. Free vibrations of laminate plates with viscoelastic layers using the refined zig-zag theory – Part 1. Theoretical background. The solution obtained using the continuation method provides the required dynamic characteristics of the plates – the frequency and the damping ratio. In the iterative part of each continuation step either the full Newton or the quasi-Newton iterative procedure are used. Several numerical examples are solved to verify accuracy, efficiency and versatility of the proposed formulation. Numerous parametric studies are carried out to determine the influence of various solution variants and combinations of material parameters onto the results describing the damping properties of laminate plates. [ABSTRACT FROM AUTHOR]

• The subspace iteration method is used to solve nonlinear eigenvalue problem. • The continuation method is used to solve the nonlinear Hermitian eigenproblem. • The method presented allows determination of only a part of eigenvalues. • The method was tested for structures with viscoelastic elements. • The models of VE elements are described with the help of fractional derivatives. The paper presents an extension of the subspace iteration method for application in systems with viscoelastic damping elements, which are described by both classical and fractional models. The presented method enables determination of only a certain number of eigenvalues and associated eigenvectors. This is very useful, especially when the considered problem has many degrees of freedom, in which case it would be very time-consuming or even impossible to determine all the eigenvalues. At the same time, from the engineering point of view, it is not necessary to know all the eigenvalues. In this paper, the subspace iteration method is used for the first time to solve the nonlinear eigenproblems which appear in the dynamic analysis of systems with viscoelastic damping elements. The proposed solution consists of assuming the number of eigenvalues to be determined, taking the initial point of iteration and then solving the nonlinear reduced eigenproblem with Hermitian matrices in each iterative loop. The solution to the reduced eigenproblem is obtained using the continuation method, and it is an additional novelty in the paper. The correctness and effectiveness of the proposed approach are illustrated with numerical examples. [ABSTRACT FROM AUTHOR]

• The refined zig-zag theory is used to solve the dynamic problem for composite beams with viscoelastic layers. • The fractional Zener model is used in analysis. • The dynamic characteristics of multi-layered viscoelastic beams are of interest. • The hierarchical FEM and the Laplace transform are used to derive the nonlinear eigenvalue problem. • The continuation method is adopted to solve the nonlinear eigenvalue problem. In the present paper the displacement based refined zig-zag theory is used, for the first time, to solve dynamic problems of composite and sandwich beams made of both elastic and VE layers. The linear constitutive relations are used to model physical properties of the VE layers. The rheological models can be in the classical or fractional versions, but the main focus is on the fractional Zener model. The dynamic characteristics of considered beams are evaluated. The hierarchical finite element and the Laplace transformation are used to derive the nonlinear eigenvalue problem and its solutions yield the required dynamic characteristics of beams. The continuation method is adopted to solve this problem numerically. Accuracy, effectiveness and efficiency of the proposed method are verified in several examples. The results of the presented calculations are summarized with some conclusions. [ABSTRACT FROM AUTHOR]