We investigate the prolonged plastic deformation of aluminum containing θ’ phase with a multistage approach combining molecular dynamics (MD), continuum modeling (CM) and discrete dislocation dynamics (DDD). The time of performed MD calculations is sufficient for about a hundred dislocation–precipitate interactions. With this number of interactions, the inclusion of θ’ is not only cut, but also scattered into individual copper atoms in an aluminum matrix. Damage to the crystal structure of inclusion and activation of the cross-slip of dislocation segments cause a decrease in acting stresses in the MD system. The rate of this effect depends on θ’ diameter and occurs faster for small inclusions. The effect of decreasing the resistance of precipitate is further introduced into the dislocation–precipitate interaction CM by reducing the precipitate effective diameter with an increase in the number of interactions. A model of dislocation–precipitate interaction accounting for the softening of inclusions is further implemented into DDD. Dependences of flow stress in aluminum with θ’ phases on volume fraction and typical diameter of precipitates are obtained. Manifestation of inclusion softening is possible in such an alloy, which leads to the flow stress decrease during deformation. The range of volume fractions and typical diameters of θ’ phases corresponding to the possible decrease in flow stress is distinguished.