neutral frequency of the beam) is considered in this example and
depicted in Fig. 3 [33,34]. A FORTRAN program is conducted and
three-node quadratic boundary elements are applied on the global
boundary of the beam. For all the numerical examples, plane stress
condition is considered.
The displacement u2 at point A computed by the standard
Newmark's algorithm are plotted in Fig. 4. In this figure the numerical
results for different values of the time step Δt are displayed.
It can be seen that all results become unstable: the smaller
the time step Δt, the earlier the instability begins. However,
simply increasing Δt does not solve the problem since larger time
steps lead to smoothed results. When using the standard Newmark's
algorithm, the most widely used combination of parameters
is δ¼1/2 and α¼1/4, which corresponds to the assumption
of a constant averaged acceleration. However, the resulting timestepping
scheme is not unconditionally stable and should therefore
be used with care