Non-woven cellulose fiber networks of low areal density are widely used in many industrial applications and consumer products. A discrete element method (DEM) modeling framework is advanced to simulate the formation of strongly anisotropic cellulose fiber network sheets in the dilute limit with simplified hydrodynamic and hydroelastic interactions. Our modeling accounts for in-plane fiber orientation and viscous drag indirectly by using theories developed by Niskanen (2018 Fundamentals of Papermaking, Trans. 9th Pulp and Paper Fundamental Research Symp. Cambridge, 1989 (FRC) pp 275–308) and Cox (1970 J. Fluid Mech. 44 791–810) respectively. Networks formed on a patterned and flat substrate are simulated for different fiber types, and their tensile response is used to assess the influence of the out-of-plane topographical pattern, specifically, on their stiffness and strength. Sheets with the same grammage and thickness, but composed with a higher fraction of softwood fiber (longer fibers with large diameter), have higher strength and higher strain to failure compared to sheets made from hardwood fibers (short fibers with small diameter). However, varying the fiber fraction produces only an insignificant variation in the initial sheet stiffness. The above simulation predictions are confirmed experimentally for sheets comprised of fibers with different ratios of Eucalyptus kraft and Northern Bleached Softwood Kraft fibers. Sheets with out-of-plane topography show an unsymmetric mass distribution, lower tensile stiffness, and lower tensile strength compared to those formed on a flat substrate. The additional fiber deformation modes activated by the out-of-plane topography, such as bending and twisting, explain these differences in the sheet mechanical characteristics.