A reasonable theoretical model of the stress arch (SA) is the key to explore its morphological evolution, but the existing theoretical models of the SA have issues, such as simple parameters and disregarding hard roof (HF). To study the influence of motion-induced fracture of the HF on morphological evolution of stress arch (MESA), the dynamic relationship between the fracturing motion of the HF and evolutionary expansion of the SA was analyzed and a calculation model of coal/rock mass interfacial stress was introduced for theoretical SA model optimization. Hence, a recursive algorithm was used to derive the morphological evolution equation of the stress arch (MEESA) under the influence of the HF, and Flac3D was used to simulate the MESA. Results indicate that the SA is formed in the overlying strata after mining, and SA morphology is controlled by the HF and experiences dynamic evolution with the fracturing motion of the HF. Before the SA expands to HF, it presents arch-shaped free expansion. After SA expands to HF, HF has the shielding effect on longitudinal expansion of the SA, thereby SA presents semi-arch-shaped transverse expansion, and HF has a constraining effect on the SA and SA goes through longitudinal expansion rapidly after HF is fractured, therefore SA presents calabash-shaped expansion. MEESA fits in well with the numerical experimental result. Hence, MEESA can be used to describe SA morphologies accurately. The conclusions in this study can provide a theoretical reference for mine pressure control.