In this research, we describe parallel numerical simulation method to compute the stress distribution of three-dimensional composite materials with fiber-reinforcements. The fibers are longitudinally surrounded by matrix, and distributed in a cube arrangement. This multi-phase polymer is modeled using hexahedral finite element discretization. The effects of the computation based on parallel CG solver with a preconditioned and the arrangement of reinforcements are discussed.