变量替换老出问题，没办法推导出式子

\[Phi]1 = D[\[Phi][x, y, t], x, x, t]
\[Phi]2 = D[\[Phi][x, y, t], x, x, x]
\[Phi]3 = D[\[Phi][x, y, t], x, x, y]
\[Phi]4 = D[\[Phi][x, y, t], y, y, t]
\[Phi]5 = D[\[Phi][x, y, t], y, y, x]
\[Phi]6 = D[\[Phi][x, y, t], y, y, y]
\[Phi]7 = D[\[Phi][x, y, t], t]
\[Phi]8 = D[\[Phi][x, y, t], x]
\[Phi]9 = D[\[Phi][x, y, t], y]
\[Phi] = \[Psi]0[y, t] + \[Psi]
\[Phi]1 + \[Phi]4 - \[Lambda]^2 \[Phi]7 - \[Phi]9*\[Phi]2 + \[Phi]9*\
\[Phi]5 + \[Lambda]^2*\[Phi]9*\[Phi]8 + \[Phi]8*\[Phi]3 + \[Phi]8*\
\[Phi]6 - \[Lambda]^2*\[Phi]8*\[Phi]9 + \[Beta]*\[Phi]8 == 0
TraditionalForm[\[Phi]1 + \[Phi]4 - \[Lambda]^2*\[Phi]7 - \[Phi]9*\
\[Phi]2 + \[Phi]9*\[Phi]5 + \[Lambda]^2 \[Phi]9*\[Phi]8 + \[Phi]8*\
\[Phi]3 + \[Phi]8*\[Phi]6 - \[Lambda]^2*\[Phi]8*\[Phi]9 + \[Beta]*\
\[Phi]8 == 0]
\[Psi] = \[Psi][\[Xi], y, \[Tau]]
\[Psi]0 = \[Psi]0[y, \[Tau]]
\[Phi] = \[Psi]0[y, t] + \[Psi][\[Xi], y, \[Tau]]
expr = \[Phi]1 + \[Phi]4 - \[Lambda]^2*\[Phi]7 - \[Phi]9*\[Phi]2 + \
\[Phi]9*\[Phi]5 + \[Lambda]^2*\[Phi]9*\[Phi]8 + \[Phi]8*\[Phi]3 + \
\[Phi]8*\[Phi]6 - \[Lambda]^2*\[Phi]8*\[Phi]9 + \[Beta]*\[Phi]8 == 0
expr /. \[Xi] -> \[Epsilon]^(
1/2) (x - Subscript[c, 0] t) \[Tau] -> \[Epsilon]^(3/2) t
这个执行以后会有TerminatedEvaluation["RecursionLimit"]我这个是哪里错了呢求大佬指点