To predict the Elasto-Plastic Behaviors (EPBs) of aligned inclusions reinforced compos-ites, this paper develops an interpolative Mori-Tanaka/Double-Inclusion (MT-DI) homogenization model with the secant formulation, and gives the numerical implementation algorithms of the devel-oped MT-DI model with the secant formulation. The Finite Element (FE) homogenization method is implemented to provide the''exact"EPBs of the composites and thus validate the MT-DI model with the secant formulation. The MT-DI model with the 2nd-order secant formulation is validated to provide the more accurate predictions, while the MT-DI model with the 1st-order secant formu-lation always gives the stiffer predictions. The results show that using the macro-stress and macro-strain as the inputs, the MT-DI model with the secant formulation gives the identical predictions. The predictions of the MT-DI model with the secant formulation vary between those of the MT and DI models with the secant formulation. Meanwhile, the MT-DI model with the secant formulation does not predict the accurate EPBs for the phases of the composites.