The shape and extent of the plastic zone can effectively reflect the degree of deformation and failure of the surrounding rock. By applying complex variable functions and the theory of elasticity, a mechanical model for the surrounding rock of a rectangular roadway is established. The analytical solution for the stress in the surrounding rock of a rectangular roadway under a non-uniform stress field is derived, and the spatial stress distribution pattern in the surrounding rock of the rectangular roadway is obtained. Furthermore, using the Mohr-Coulomb criterion, the boundary equation of the plastic zone in the surrounding rock of a rectangular roadway is derived, and the morphological characteristics of the plastic zone are obtained. The results indicate that the plastic zones on both sides of the rectangular roadway are arc-shaped. The asymmetry of the stress environment has a significant impact on the plastic zones in the roof and floor of the rectangular roadway, causing the shape of the plastic zones to resemble a “petal” pattern. The lateral pressure coefficient and the height-to-width ratio have a significant impact on the stress distribution and the morphological characteristics of the plastic zone in the surrounding rock of the rectangular roadway. By validating the computational accuracy of the analytical solution through numerical simulations and conducting an engineering applicability analysis based on a coal mine roadway, it has been demonstrated that the analytical calculation method is both rational and possesses good engineering applicability. This provides reliable theoretical guidance for the design calculations and construction planning of similar roadway projects.