Cross-laminated timber (CLT) wall systems are composed of massive timber panels that are fastened together and to the horizontal elements (foundations or intermediate floors) with step joints and mechanical connections. Due to the high in-plane stiffness of CLT, the shear response of such systems depends strongly on the connections used. This paper proposes a numerical model capable of predicting the mechanical behavior and failure mechanisms of CLT wall systems. The wall and the element to which it is anchored are simulated using three-dimensional (3D) solid bodies, while the connections are modeled as nonlinear hysteretic springs. Typical racking tests of wall systems are reproduced by varying the assumptions used to schematize the behavior of the connections. Results are compared with test data published in the literature, and the differences are discussed. The influence of the boundary conditions (vertical load applied on top of the wall and friction at its base) and aspect ratio of the panel are investigated via a parametric numerical study. Finally, the performance of a wall system assembled with two CLT panels is analyzed, highlighting how the properties of the anchoring connections and vertical step joints affect the load-displacement response and energy dissipation.