How to simulate a Double Pendulum with a Motor
I'm trying to find a specific formula for a variant of the 2D double pendulum experiment, but I'm afraid my physics skills are not up to the task. I'm hoping someone here can help me :). The situation is like in the double pendulum task, but with the following changes:
* There is no gravity
* The joint between the two rods (motor joint) cannot rotate freely, but an internal motor controls the angle between them. If the motor isn't moving, the two rods are basically a rigid body.
* I don't know the actual shape of the two rods, but I do know their individual masses, inertia and the location of the center of mass.
It looks something like this: (The rods look like rectangles here, but their actual shape is unknown)
α and β are the angles between the rods and the y-axis (one of them is negative here). Note that the joint, that fixes the upper joint in the picture to the world, (static joint) can rotate freely.
Now I'm not looking for the general formula to simulate this situation, but I'm interested in the following specific case: - First, both rods are at rest with the motor joint still at angle γ=β-α. 782
- Then, the motor spins a little bit and stops again. Once it stops, both rods will be at rest too, but the angles will be different: γ'=β'-α'.
The situation might look something like this now:
Note that the motor joint has a different angle, because the motor moved in the meantime. The static joint has a different angle as well, because the motor applies a force to both rods it is attached to. But before and after the motor movement, all velocities are zero.
Now, given α, β, γ' and all the mass information, how can I calculate α' (or β')?
It's fine if the formula assumes that the difference between γ and γ' is very small. I already tried to figure this out on my own with derivations of the double pendulum formulas on the internet, but I don't know how to deal with the motor joint in those formulas. I also implemented that situation in a 2D physics engine (box2d) and tried to deduce the formula in that simulation "experimentally", but couldn't figure it out. 782