Mechanical behaviour of 3D-printed composite parts is affected by the volume fraction, aspect ratio and type of fibre reinforcement. Although in the literature experimental approaches have been used to characterise the effects of the above factors on the mechanical properties of 3D printed parts, time and cost of the manufacturing process as well as the uncertainty associated with a large number of experimental techniques are the key issues. This study aims to address these challenges by developing a methodology based on a multi-scale Finite Element (FE) analysis of representative volume element (RVE) of 3D printed composite parts to predict the effective orthotropic properties. To account for the effects of fibre features, RVEs were modelled considering variables of volume fraction, aspect ratios and type of short fibres. To study the main and interaction effects of the above variables on the mechanical properties of 3D printed composite parts, a structured approach based on the Design of Experiments is used. The FE stress analysis of the RVE provides an understanding about the potential failure modes such as interfacial debonding between fibres and matrix, interlayer and intralayer delamination that may occur in load-bearing 3D printed composite parts. The FE computed mechanical properties are validated against experimental data through a series of mechanical testing of flexure, Iosipescu, and short beam shear which were conducted in conjunction with the Digital Image Correlation technique. As a result, certainty is obtained in using the proposed approach for a fast iterative design of 3D printed composite parts prior to industrial applications.