Moment of inertia - inverted pendulum
Moment of inertia - inverted pendulum
(OP)
Hi,
I have been trying to understand the dynamics behind an inverted pendulum from this webpage:
h ttp://www. library.cm u.edu/ctms /ctms/simu link/examp les/pend/p endsim.htm
I don't get what the moment of inertia (I) of the rod is. I was thinking it should be about the end of the rod: m*(2L)^2, since L is defined as half the length of the rod, but the listed I value seems to be solved for a rod spinning about its center of gravity: 1/12*m*(2L)^2. Can someone explain this issue to me? Why should I be about the rod's cg?
Thanks,
- Paul
I have been trying to understand the dynamics behind an inverted pendulum from this webpage:
h
I don't get what the moment of inertia (I) of the rod is. I was thinking it should be about the end of the rod: m*(2L)^2, since L is defined as half the length of the rod, but the listed I value seems to be solved for a rod spinning about its center of gravity: 1/12*m*(2L)^2. Can someone explain this issue to me? Why should I be about the rod's cg?
Thanks,
- Paul





RE: Moment of inertia - inverted pendulum
Find the moment of inertia about the center of the rod. It is twice the moment of inertia of one half the rod about the center:
I1=2*int(mu*r^2*dr,r=0..L)
D
I1=2*[r^3/3], evaluated at r=L minus same thing at r=0:
I1 = 2/3 L^3*mu
substitute: mu=m/(2*L) into I1:
I1 := 1/3 L^2 m
substitute L=0.3,m=0.2 into I1:
I1 = 0.006
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
The particular choice of 2 coordinates x and theta as shown results in two equations (F = M*d^2x/dt^3 and T = I * d^2theta/dt^2) which are already decoupled. i.e. they are the modal variables. So it is very easy to solve for derivative of state variables. If you chose x and theta as angle from the bottom, they would not end up being decoupled, but with a very small amount of additional algebra you can still solve for derivatives of state variables to feed it into Matlab to get the same solution. (But you'd still get same solution)
The best choice of state variables to give the simplest, most direct and most intuitive path to solution can be a tricky thing when you have angles and displacements. For this simple problem, I think there is some intuition that suggests that a rigid body has 2 independent modes: one for translation of center of mass and one for rotation about center of mass, which leads to this selection without much fuss. I think maybe Lagrange can help provide insight about selection of variables and developing equation for more complicated systems. I studied it a few years ago and remember I was impressed by what you could do to develop equations of motion more simply. But I've forgotten most of it.
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
I'm not sure there is a big benefit to selecting center of bar for summing moments as opposed to bottom of bar for summing moments (maybe someone can comment?). Again you will come up with same solution either way (as long as you apply equations correctly.)
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
RE: Moment of inertia - inverted pendulum
Fnetx on M = M*d^2(x)/dt^2
Fnetx on m = m*d^2(xp)/dt^2
Fnety on m = m*d^2*(yp)/dt^2
Tnet on I = I*d^2*(theta)/dt^2
On first glance it looks like a 4DOF system. However there is one constraint which is that the positions of pendulum COM (xp,yp) must always remains a distance L from the position of mass M (x,0). Once you consider that constraint, I believe it turns out o be a 3DOF problem as we can see from the fact that if we know one variable theta, we can compute the two variables (xp,yp) and we still know the complete condition of the system (including the kinetic energy of each of the three inertia's m,M, I).
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
Once you consider that constraint, I believe it turns out to be a [B]3[B]DOF problem as we can see from the fact that if we know one variable theta, we can compute the two variables (xp,yp) and we still know the complete condition of the system (including the kinetic energy of each of the three inertia's m,M, I).
should have ben:
Once you consider that constraint, I believe it turns out to be a 2DOF problem as we can see from the fact that if we know two variables: X and theta (along with their derivatives), we can still compute (xp,yp) and we still know the complete condition of the system (including the kinetic energy of each of the three inertia's m,M, I).
In other words, it is a 2DOF system since 2 variables (and their derivatives) completely describe the system energy.
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
Mass moments of inertia are important for a problem like this.
The inertia should be taken about the axis of rotation, always. In this case it's the pin joint at the bottom attaching the rod to the cart. This joint constrains the rod to rotate only about this axis. It of course also translates due to the motion of the cart, which comprises the 2nd DOF.
The moment of inertia in this case is 1/3*mL^2 as seen here
http:
For more information about moments of inertia and their relations to dynamic systems there are a few dynamics texts books I would refer you to.
Fe
RE: Moment of inertia - inverted pendulum
The inverted pendulum system is a very good model to learn how to solve dynamical systems and control them. I remember it was the model I chose to learn nonlinear control theory on at uni.
Fe
RE: Moment of inertia - inverted pendulum
I looked it over as it is wiki so there is a possibility of mistakes. But it looks fine to me.
http://en.wikipedia.org/wiki/Inverted_pendulum
Fe
RE: Moment of inertia - inverted pendulum
=====================================
(2B)+(2B)' ?
RE: Moment of inertia - inverted pendulum
1/2Icgw^2 + 1/2mx'^2 +1/2my'^2
where x' and y' are the orthogonal coordinates of velocity at the cg
while 1/2 Iw^2 is the KE only when there is no translation of the pivot point.
RE: Moment of inertia - inverted pendulum
Thanks. You are correct. I guess my wording was not 100% correct. It does come out easier if you choose a pin joint when deriving the equation of motion as some of the forces applied to the rigid body don't create a moment there. This goes for most mechanisms.
Thanks again for the correction.
zekeman,
I agree totally when we are deriving the DE's using an energy method. Lagrange or virtual work ect.
But using the Newton-Euler method you can choose any point (and inertia), where the joint is the most efficient choice in many cases.
In the pendulum example we end up with the inertial part of the DE as torque=(angular accel)*(inertia about joint) in the DE which is intuitive and correct. Then the inertia is just Icg+m*l^2 by parallel axis theorem.
h
I'm sure we all agree, we just have different ways of thinking about it.
Fe
RE: Moment of inertia - inverted pendulum
FeX32, if I analyzed about the joint, would the DE simply be thetaddot = 1/(Icg + mL^2) * mgLsin(theta)? And everything else would remain the same in the link I showed?
RE: Moment of inertia - inverted pendulum
Summing moments about the joint you directly get
(Icg+mL^2)*thetadd+mgL*sin(theta)+mLxdd*cos(theta)=0
You then sum forces in the most principle direction. So yes the same as that link.
To solve the DE's I sugest using both Simulink and ODE45.
Fe
RE: Moment of inertia - inverted pendulum
(Icg+mL^2)*thetadd-mgL*sin(theta)+mLxdd*cos(theta)=0
Fe
RE: Moment of inertia - inverted pendulum
I agree with your comment,(18 Jul 11 0:07); equations both ways are equivalent. Guess my personal preference got in the way.
RE: Moment of inertia - inverted pendulum
Fe