A micromechanical model was developed to account for the particular microtexture of upper bainite in low alloy steels, i.e. the non-random spatial distribution of variants within a given former austenite grain. A self-consistent scheme and an Eshelby approach, considering both transformation shape strain and viscoplastic strain as eigenstrains, was applied to estimate coupling between parent austenite and two or more bainite variants without any applied stress. Model predictions concerning self-accommodation between variants are sensitive to the plane of the first “lattice invariant shear” in the crystallographic model used to determine the shape strain. No obvious effect of the constitutive equations of phases and of the other model parameters was found.