A new beam finite element based on a new discrete-layer laminated beam theory with sublaminate first-order zig-zag kinematic assumptions is presented and assessed for thick and thin laminated beams. The model allows a laminate to be represented as an assemblage of sublaminates in order to increase the model refinement through the thickness, when needed. Within each sublaminate, discrete-layer effects are accounted for via a modified form of DiSciuva's linear zig-zag laminate kinematics, in which continuity of interfacial transverse shear stresses is satisfied identically. In the computational model, each finite element represents one sublaminate. The finite element is developed with the topology of a four-noded rectangle, allowing the thickness of the beam to be discretized into several elements, or sublaminates, if necessary, to improve accuracy. Each node has three engineering degrees of freedom, two translations and one rotation. Thus, this element can be conveniently implemented into general purpose finite-element codes. The element stiffness coefficients are integrated exactly, yet the element exhibits no shear locking due to the use of a consistent interdependent interpolation scheme. Numerical performance of the current element is investigated for an arbitrarily layered beam, a symmetrically layered beam and a sandwich beam with low and high aspect ratios. The comparisons of numerical results with elasticity solutions show that the element is very accurate and robust. (C) 1997 Elsevier Science Ltd.