One way is to make a function. Change n
and m
to pick the mode of vibration using the slider. Move the time
slider to change the time. The x,y
are hardcoded ranges. The plot range is fixed in order to better capture the motion scale.
Change as needed.

code
ClearAll[n, m, x, y, t, f];
f[n_, m_, t_, x_, y_] =
FullSimplify[(t Cos[(m^2 Pi^2)/9 + (n^2 Pi^2)/9] +
t Sin[(m^2 Pi^2)/9 + (n^2 Pi^2)/9]) Sin[(n Pi x)/3] Sin[(m Pi y)/
3]]
Manipulate[
Module[{x, y},
Plot3D[f[n, m, t, x, y], {x, 0, 3}, {y, 0, 3},
PerformanceGoal -> "Quality",
PlotRange -> {Automatic, Automatic, {-5, 5}},
PlotLabel -> Row[{"time ", t}]
]
],
{{t, 0, "time"}, 0, 10, .1, Appearance -> "Labeled"},
{{n, 2, "n"}, 1, 10, 1, Appearance -> "Labeled"},
{{m, 3, "m"}, 1, 10, 1, Appearance -> "Labeled"},
TrackedSymbols :> {t, n, m}
]