This study, with the engineering background of the design of a stope involving a sublevel mining method in a certain underground metal mine, explored the application of stress-solving methods based on the complex variable function theory in actual engineering. Three mathematical calculation models based on the functions of a complex variable were established. Through triangle interpolation, mapping functions of a plane with a roadway section and a plane with the stope section were determined. An improved Schwarz alternating method was adopted to study the stability of the roadway and the influence of an adjacent roadway from the perspective of the stress field. In addition, in light of the distribution characteristics of a gangue in the stope, the design parameters of a pillar were optimized, with the pillar’s optimal dimensions determined. The results showed that when the magnitudes of two far-field principal stresses in the rock mass are relatively close, the distribution of the surrounding rock stress around the roadway is more uniform, and tensile stress is less likely to occur. The excavation of a neighboring roadway exacerbates to some extent the side stress of the other roadway, especially the compressive stress concentration on the side closer to the neighboring roadway. However, when the distance between the two roadways is significantly larger than the roadway size, this effect is not pronounced. In the engineering case studied in this research, the thickness of the pillar is approximately linearly positively correlated with the safety factor of the pillar approximately linearly negatively correlated with the recovery rate. Overall, this research explored the application of the complex variable function theory in an underground mine design, demonstrating its accuracy and practicality.