Consider some sort of entity/object with generalized coordinate $q$ and generalized velocity $\dot{q}$. Now suppose that we know the form of the object's Lagrangian when it's in free motion, not interacting with anything. Denote that by $\mathcal{L}_{free}$. Now we create conditions for the object to interact with its surroundings, and the Lagrangian becomes $\mathcal{L}$. In classical mechanics, and even in electrodynamics, we have that $\mathcal{L}=\mathcal{L}_{free}-U$, where $U$ is a function of the generalized position, velocity and time that's usually called the potential, and its form depends only on the type of interaction. Okay. But why?
Can we reach that conclusion without starting from the equations of motion (i.e. Newton's Second Law)? That we can always describe interactions by adding some sort of extra term to the Lagrangian of free motion?
UPDATE 1: I've fiddled around with the concept a bit, and started to think about lagrangian multipliers; describing interactions as some sort of constraint $f=0$, and finding stationary points for the free action, $$S_{free}=\int_{t_1}^{t_2}\mathcal{L}_{free}\,\mathrm{d}t$$ subject to the constraint $f=0$, which is equivalent to finding the stationary points of a new functional $$S=\int_{t_1}^{t_2}\mathcal{L}_{free}\,\mathrm{d}t-\lambda f=\int_{t_1}^{t_2}\mathcal{L}_{free}-U\,\mathrm{d}t$$ where $U\equiv -\lambda\frac{df}{dt}$, and so we define a new Lagrangian $\mathcal{L}=\mathcal{L}_{free}-U$.
But this still isn't satisfactory, given that I had to assume that the free action is still 'stationary' despite knowing nothing about how interactions may affect that original Lagrangian.
UPDATE 2: I see someone mentioned coupling constants and whatnot. I've elaborated on that and got no satisfactory results. Suppose we can modulate the strength of our interaction through some constant $\alpha$. Then we can say that our Lagrangian is a perturbation of the original, of the form $$\mathcal{L}=\mathcal{L}_{free}-\alpha U,$$ for $U\equiv-\frac{\partial\mathcal{L}}{\partial\alpha}$. If we then absorb the 'correct' value of our constant (i.e. the value that gets us the right equations of motion) into $U$, we get that $\mathcal{L}=\mathcal{L}_{free}-U$. The problem with this approach is that some interactions can't be modulated to the proper intensity, e.g. very quick isentropic processes in a gas or anything with hysteretic behavior, and therefore this model is at least flawed.