Solve for dimensions of rectangular prism given inertia


by krenzke
Tags: inertia, inverse, iterative
krenzke
krenzke is offline
#1
Dec25-10, 09:42 PM
P: 2
I am trying to solve something of an inverse problem for inertia, and having a real tough time with it. Any suggestions would be helpful.

For a rectangular prism the mass and moments of inertia are:
m = rho*A*B*C
Ixx = m/12*(B^2 + C^2)
Iyy = m/12*(A^2 + C^2)
Izz = m/12*(A^2 + B^2)

Where m is mass, rho density, A width, B length, C height, and Ixx, Iyy, Izz are the principal moments of inertia. My problem involves a hollow rectangular prism, with wall thickness t. So the equations then become:

mOuter = rho*A*B*C
mInner = rho*(A-2*t)*(B-2*t)*(C-2*t)
m = mOuter - mInner
Ixx = mOuter/12*(B^2 + C^2) - mInner/12*((B-2*t)^2 + (C-2*t)^2)
Iyy = mOuter/12*(A^2 + C^2) - mInner/12*((A-2*t)^2 + (C-2*t)^2)
Izz = mOuter/12*(A^2 + B^2) - mInner/12*((A-2*t)^2 + (B-2*t)^2)

To see how these equations work imagine taking a rectangular prism and subtracting out a smaller rectangular prism from the inside. Given length (A), width (B), height(C), thickness (t), and density (rho) I can easily solve for the mass and moments of inertia. What I would like to do is given m, rho, Ixx, Iyy, and Izz solve for A,B,C, and t. I tried doing the hand-calc to solve the 4 equations for the 4 unknowns and quickly ran into a seemingly intractable 15th-degree polynomial in t with tons of unknown coefficients. If anyone has any ideas on how to solve this I would greatly appreciate it. It doesn't have to be a closed form solution (though that would be best), iterative or other methods would work as well.

Any ideas?
Phys.Org News Partner Physics news on Phys.org
Physicists consider implications of recent revelations about the universe's first light
Vacuum ultraviolet lamp of the future created in Japan
Grasp of SQUIDs dynamics facilitates eavesdropping
Filip Larsen
Filip Larsen is online now
#2
Dec26-10, 09:52 AM
PF Gold
Filip Larsen's Avatar
P: 946
Welcome to PF!

If you want to pursue a numerical solution coded by yourself (as opposed to using a numerical solver in, say, Matlab or Mathematica), one of the standard tools in that toolbox would be to try solve it using Newton-Raphson [1] on your system of equations. If you want to try that, you should probably write it up as a system of equations of 4 unknowns and 4 knowns, as density is just a scale factor, and note that you can break symmetry in A, B and C by restricting the solution to [itex]A \geq B \geq C > 2t > 0[/itex] and make an initial guess that is asymmetric, like for instance [itex] B = A/2, C = A/4, t = A/16[/itex] which should give the initial guess [itex]A^3 = \frac{512}{43}\frac{m}{\rho}[/itex].

[1] http://en.wikipedia.org/wiki/Newton'...s_of_equations
krenzke
krenzke is offline
#3
Dec27-10, 12:51 AM
P: 2
Thanks for the reply! I finally did try a newton-rhapson method coded in c++, and it seems to work pretty well. I've been using the eigenvalues of the inertia tensor as the initial guess, with some fraction of the smallest value as the thickness. I have noticed that it will periodically converge to a non-physical value (i.e. A,B,C or t less than zero), but in this case taking the absolute value of the non-physical answer and putting that back into the iterative solver seems to come up with a different (and realistic) solution. Not sure if there is a better way to ensure that the solution lands in the positive quadrant...perhaps a boundary function of some sort?

Filip Larsen
Filip Larsen is online now
#4
Dec27-10, 02:58 AM
PF Gold
Filip Larsen's Avatar
P: 946

Solve for dimensions of rectangular prism given inertia


Quote Quote by krenzke View Post
I have noticed that it will periodically converge to a non-physical value (i.e. A,B,C or t less than zero), but in this case taking the absolute value of the non-physical answer and putting that back into the iterative solver seems to come up with a different (and realistic) solution. Not sure if there is a better way to ensure that the solution lands in the positive quadrant...perhaps a boundary function of some sort?
I'm probably just stating the obvious here, but anyway:

If the non-physical solutions pop up after a large jump in state (due to a small-magnitude Jacobian giving rise to a large state change) you could try limit the magnitude of the state change.

If the non-physical solution appears during a normal "slow" convergence sequence, then your initial guess must belongs to the attraction basin of that solution and a different initial guess is (as you say) obviously needed, probably one where the value for A, B and C are more equal in magnitude. If non-physical solutions are a real problem a practical "work-around" may be to map out the attraction basins of different initial guesses to see if there is a class or pattern of initial guesses that always converge to the physical solution and then stick to those.

If I recall correctly, there should also be a fair bit of convergence theory which may help you if you want to analyse convergence in a more theoretical way, although I'm not sure how well it will apply to your non-linear system.


Register to reply

Related Discussions
Angle of diagonal of a rectangular prism Differential Geometry 0
moment of inertia of a rectangular prism Classical Physics 15
Quantum - Electron in an infinite rectangular prism well Advanced Physics Homework 2
Moment of inertia for rectangular prism General Engineering 2
Related Rates (rectangular prism) Calculus & Beyond Homework 4