Do you ever work on a problem that you know ought to be simple, ought to have a straight forward answer, and yet it just does not seem to work out? That is where I am right now.
For about a year, off and on, I have been working on a particular mechanics problem, not a particularly esoteric problem at all, that I should be able to work in a matter of a few hours and that would be the end of it. It is the sort of problem that an undergraduate ME should be able to put together and have a satisfactory solution in perhaps a couple of days. So what is the great difficulty?
Well, I am insisting on checking it, and that is where the rub comes. I can get numbers without too much difficulty. It is just a matter of solving Newton’s Second Law equations for a system of three bodies simultaneously, a problem that produces a system of seven linear equations in this case. The unknowns are six internal forces and the system acceleration. With the computer, this is all pretty straight forward, and crunching the numbers is not that big a problem. My program runs in much less than a second, so time is not a problem. There are some kinematic calculations to be done before the force solution is made which makes the use of the computer very advantageous, but the calculations fairly fly. I have it programed both in Maple 8 and also in BASIC and I get the same results in both.
But then there is that matter of checking it. Newton’s Second Law is pretty nearly infallible for mechanical systems, but the person applying it is not. So it is important to check the result. I can substitute the solutions back into the expressions for Newton’s Second Law, and they are satisfied, but all that says is that I got a good numerical solution. It does not say that I applied Newton’s Second Law correctly. What else can I do?
One of the happy alternatives to Newton’s Second Law for dyanical systems is the Lagrange Equation of motion, an energy based approach to determining the motion of a dynamic system. In many cases, it is significantly easier to apply than Newton’s Second Law and less prone to errors. Generally speaking, all that is required is to write the system kinetic coenergy, take a few derivatives, write the virtual work of the external forces, and viola’, the equation of motion appears! From this expression, the system acceleration is quickly solved and this provides a check on the Newton’s Law formulation.
Except when it does not! In the problem at hand, for a year or so, I have not been able to reconcile these two. I know that I am getting good numerical solutions to the Newton’s Law equations, so that is not the problem. In some cases, the accelerations determined by the two approaches are fairly close to each other; in other cases, they are significantly different.
This is truly a “head-banging” situation. I trust both of these approaches, and when they are in disagreement, it means one of them must be incorrect (it really means that I have incorrectly applied one of them). It looks like it is time to lay this problem aside again for a few weeks and hope for new inspiration. What to do?