想将红色方框区域的表达式中的sigemaZ用其他参数表达出来，请问用solve（红色框选部分，sigemaZ）来求解为什么会一直运行中呢？或者有其他可以求解的方式么？

\[Epsilon]z=(-r*cos\[Alpha])/R;
\[Sigma]c=(p/t+(cos\[Alpha] *\[Sigma]z)/(R-y))*r;
\[Epsilon]t=(\[Sigma]z+\[Sigma]c)*\[Epsilon]z/(\[Sigma]c-2\[Sigma]z);
\[Epsilon]c=(2\[Sigma]c-\[Sigma]z)*\[Epsilon]z/(2\[Sigma]z-\[Sigma]c);
\[Sigma]z^2-\[Sigma]z*\[Sigma]c+\[Sigma]c^2==K^2*((4/9)*((2*r^2*(cos\[Alpha])^2+R^2*(\[Epsilon]t^2+\[Epsilon]c^2-\[Epsilon]t*\[Epsilon]c))/R^2))^n
Solve[\[Sigma]z^2-r \[Sigma]z (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))+r^2 (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))^2==(4/9)^n K^2 ((2 cos\[Alpha]^2 r^2+R^2 ((cos\[Alpha]^2 r^2 (\[Sigma]z+r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y)))^2)/(R^2 (-2 \[Sigma]z+r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y)))^2)-(cos\[Alpha]^2 r^2 (\[Sigma]z+r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))) (-\[Sigma]z+2 r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))))/(R^2 (2 \[Sigma]z-r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))) (-2 \[Sigma]z+r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y))))+(cos\[Alpha]^2 r^2 (-\[Sigma]z+2 r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y)))^2)/(R^2 (2 \[Sigma]z-r (p/t+(cos\[Alpha] \[Sigma]z)/(R-y)))^2)))/R^2)^n,\[Sigma]z]