  Laminated glass elements, which consist of stiff elastic glass layers
connected with a compliant viscoelastic polymer foil, exhibit geometrically
non-linear and time/temperature-sensitive behavior. In computational modeling,
the viscoelastic effects are often neglected or a detailed continuum
formulation typically based on the volumetric-deviatoric elastic-viscoelastic
split is used for the interlayer. Four layerwise beam theories are introduced
in this paper, which differ in the non-linear beam formulation at the layer
level (von K\'{a}rm\'{a}n/Reissner) and in constitutive assumptions for the
interlayer (a viscoelastic solid with the time-independent bulk modulus/Poisson
ratio). We perform detailed verification and validation studies at different
temperatures and compare the accuracy of the selected formulation with
simplified elastic solutions used in practice. We show that all the four
formulations predict very similar responses. Therefore, our suggestion is to
use the most straightforward formulation that combines the von K\'{a}rm\'{a}n
model with the assumption of time-independent Poisson ratio. The simplified
elastic model mostly provides a response in satisfactory agreement with full
viscoelastic solutions. However, it can lead to unsafe or inaccurate
predictions for rapid changes of loading. These findings provide a suitable
basis for extensions towards laminated plates and glass layer fracture, owing
to the modular format of layerwise theories.
