| Abstract: | A finite element method based on the virtual work principle to determine the steady state response of frams in free or forced periodic vibration is introduced. The axial and flexural deformations are coupled by mean of the induced axial force along the element. The spatial discretization of the deformations is achieved by the usual finite element method and the time discretization by Fourier coefficients of the nodal displacements. No unconventional element matrices are needed. After applying the harmonic balance method, a set of non-linear algebraic equations of the Fourier coefficients is obtained. These equations are solved by the Newtonian iteration method in terms of the Fourier coefficient increments. Nodal damping can easily be included by a diagonal damping matrix. The direct numerical determination of the Fourier coefficient increments is difficult owing to the presence of peaks, loops and discontinuities of slope along the amplitude-frequency response curves. Parametric construction of the response curves using the phase difference between the response and excitation is recommended to provide more points during the rapid change of the phase (i.e. at resonance). For undamped natural vibration, the method of selective coefficients adopted. Numerical examples on the Duffing equation, a hinged–hinged beam, a clamped–hinged beam, a ring and a frame are given. For reasonably accurate results, it is shown that the number of finite elements must be sufficient to predict at least the linear mode at the frequency of interest and the number of harmones considered must satisfy the conditions of completeness and balanceability, which are discussed in detail. |
