A three-dimensional numerical avian model was developed to investigate the unsteady and turbulent aerodynamic performance of flapping wings for varying wingbeat frequencies and flow velocities (up to 12 Hz and 9 m s(-1)), corresponding to a reduced frequency range of k=0.22 to k=1.0 and a Reynolds number range of Re=16,000 to Re=50,000. The wings of the bird-inspired model consist of an elastic membrane. Simplifying the complicated locomotion kinematics to a sinusoidal wing rotation about two axes, the main features of dynamic avian flight were approximated. Numerical simulation techniques of fluid-structure interaction (FSI) providing a fully resolved flow field were applied to calculate the aerodynamic performance of the flapping elastic wings with the Reynolds averaged Navier-Stokes (RANS) approach. The results were used to characterize and describe the macroscopic flow configurations in terms of starting, stopping, trailing and bound vortices. For high reduced frequencies up to k=0.67 it was shown that the wake does not consist of individual vortex rings known as the discrete vortex ring gait. Rather, the wake is dominated by a chain of elliptical vortex rings on each wing. The structures are interlocked at the starting and stopping vortices, which are shed in pairs at the reversal points of the wingbeat cycle. For decreasing reduced frequency, the results indicate a transition to a continuous vortex gait. The upstroke becomes more aerodynamically active, leading to a consistent circulation of the bound vortex on the wing and a continuous spanwise shedding of small scale vortices. The formation of the vortices shed spanwise in pairs at the reversal points is reduced and the wake is dominated by the tip and root vortices, which form long drawn-out vortex structures.