Consider a damped, driven harmonic oscillator:
$$m\frac{d^2x}{dt^2} + \beta\frac{dx}{dt} + kx = F(t) $$
I want to write this equation for time-varying $m$, $\beta$ and $k$ which are the mass, viscous damping coefficient and spring coefficient respectively.
In some papers, for the time-varying version of this equation the authors simply write:
$$m(t)\frac{d^2x}{dt^2} + \beta(t)\frac{dx}{dt} + k(t)x = F(t) $$
however this seems not quite right to me because according to Newton's second law, force equals to the time-derivative of momentum, then for time-varying mass:
$$F= \frac{d}{dt}\left(m\frac{dx}{dt}\right) = m\frac{d^2x}{dt^2} + \frac{dm}{dt}\frac{dx}{dt}$$
So there should be an extra term in the equation:
$$m(t)\frac{d^2x}{dt^2} + \left(\beta(t)+\frac{dm}{dt}\right)\frac{dx}{dt} + k(t)x = F(t) $$
I am not sure whether there should be extra terms for time-varying $\beta$ and $k$ too.
In short, I would like to know what time-dependent harmonic oscillator equation looks like.