mA = 1; mB = 3;
h = 1;
l = 1;
time = 0.6;
\[Theta]0 = ArcTan[h/l];
s = NDSolve[{h/(Sin[\[Theta][t]])^2 \[Theta]'[t] == Sqrt[(
20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/(
mA + mB*(Cos[\[Theta][t]])^2)], \[Theta][
0] == \[Theta]0}, \[Theta], {t, 0, time}];
\[Theta] = \[Theta] /. s[[1]];
Plot[{Sqrt[(20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/(
mA + mB*(Cos[\[Theta][t]])^2)],
Sqrt[(20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/(
mA + mB*(Cos[\[Theta][t]])^2)]*Cos[\[Theta][t]]}, {t, 0, time}]
Clear["Global`*"]

这是一个具有实际物理背景的方程，照理说，肯定图像的呀

首先，不是没图像，而是图像几乎就是θ=0，然后，你在h=1时观察到的那个非零解其实是数值误差。你提升下WorkingPrecision再解试试。
如果这个方程的解不应如此，那一定是你的方程错了，请再仔细查查。