In order to investigate effects of anisotropic fiber packing on stresses in composites, a Volume Integral Equation Method is applied to calculate the elastostatic field in an unbounded isotropic elastic medium containing multiple orthotropic inclusions subject to remote loading, and a Mixed Volume and Boundary Integral Equation Method is introduced for the solution of elastostatic problems in unbounded isotropic materials containing multiple anisotropic inclusions as well as one void under uniform remote loading. A detailed analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out for square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively. Also, an analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out, when it is assumed that a void is replaced with one inclusion adjacent to the central inclusion of square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively, due to manufacturing and/or service induced defects. The effects of random orthotropic fiber packing on stresses at the interface between the isotropic matrix and the central orthotropic inclusion are compared with the influences of square and hexagonal orthotropic fiber packing on stresses. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with multiple orthotropic inclusions and one void, it will be established that these new methods are very accurate and effective for investigating effects of general anisotropic fiber packing on stresses in composites.