For the past 3 years our group has been developing numerical models of ballistic impact into layered fibrous systems consisting of orthogonal UD tapes as well as plain-weave fabrics. The projectiles have included right-circular cylinders (RCCs) and spherical shapes. Viscoelastic interactions and slip between both the projectile and the target and between target layers and individual yarns have been modeled, as well varying degrees of yarn crimp in both the UD layers and the fabrics. Besides exploring new phenomena resulting from layer-to-layer interactions, these models have been useful for validating previous analytical results in the literature, where various simplifying assumptions were necessary. It has been possible to check various scaling relationships, that are difficult and costly to investigate experimentally. In the past, the numerical models have assumed an axisymmetric framework whereby the numerical solution of only one quarter of the plane was necessary. In this paper we report on our recent work in modeling projectile impact into the full plane, thus making it possible for impact to occur at an arbitrary angle to the target as well as for multiple, spaced impacts to occur at the same or different times. We can also investigate boundary effects on the target (clamped, free) where projectile impact occurs close to one boundary or near a corner, as well as the effect of limited target dimensions especially when the target has relatively low areal density. In this case projectile velocities also tend to be low, so that tension waves have sufficient time to reflect from the boundary back to the point of impact, thereby affecting the local strains and back-face deformation. The paper will discuss recent results on the various aspects above.