Skip to main content
21-Topaz II
March 28, 2014
Solved

Model with simple frustrum - volume calculation wrong

  • March 28, 2014
  • 1 reply
  • 7746 views

I make a lot of models that involve pyramidal frustrums (transition of a rectangular section to another rectangular section, linearly). Rather than use the measured value calculated by Creo, I thought I'd have a go at calculating the volume directly from the dimensions.

When I check my calculations, the volume calculated directly from the dimensions differs from the Creo analysis result.

I minimized the accuracy value, but still get the exact same result.

Here are my simple values

Base rectangle: 5 X 3 -> Base area: 15

Top rectangle: 4 X 2 -> Top area: 8

Height: 1

Frustrum of a pyramid volume is

Volume = ( Area Base + Area Top + SQRT ( Area Base * Area Top ) ) * Height / 3 = 11.3182

Creo, on the other hand, gives a volume of 11.3333.

Anyone have any insights into this?

This is simple shape, and kind of scares me about what might be happening with more complicated geometry.

    Best answer by dschenken

    I think Creo is right.

    I'm not sure where the formula came from but ...

    L = l_base - ((l_base-l_top)/delta_z) * Z

    W = w_base - ((w_base-w_top)/delta_z) * Z

    A = L * W =

    l_base*w_base

    - (l_base*((w_base-w_top)/delta_z) + w_base*((l_base-l_top)/delta_z)) * Z

    + ((l_base-l_top)/delta_z)*((w_base-w_top)/delta_z)*Z*Z

    Integrating:

    V =

    l_base*w_base*Z

    - (l_base*((w_base-w_top)/delta_z) + w_base*((l_base-l_top)/delta_z)) * Z*Z/2

    + ((l_base-l_top)/delta_z)*((w_base-w_top)/delta_z)*Z*Z*Z/3

    Plugging it all in I get 11.33...

    1 reply

    dschenken1-VisitorAnswer
    1-Visitor
    March 28, 2014

    I think Creo is right.

    I'm not sure where the formula came from but ...

    L = l_base - ((l_base-l_top)/delta_z) * Z

    W = w_base - ((w_base-w_top)/delta_z) * Z

    A = L * W =

    l_base*w_base

    - (l_base*((w_base-w_top)/delta_z) + w_base*((l_base-l_top)/delta_z)) * Z

    + ((l_base-l_top)/delta_z)*((w_base-w_top)/delta_z)*Z*Z

    Integrating:

    V =

    l_base*w_base*Z

    - (l_base*((w_base-w_top)/delta_z) + w_base*((l_base-l_top)/delta_z)) * Z*Z/2

    + ((l_base-l_top)/delta_z)*((w_base-w_top)/delta_z)*Z*Z*Z/3

    Plugging it all in I get 11.33...

    1-Visitor
    March 29, 2014

    I should be able to simplify the integral to be the same as the original formula, but so far no luck. I looked at the Wikipedia page that showed the derivation of the original formula, but it doesn't yield terms I can yet convert.

    I get 11.3182 from the original formula, but don't see the reason for the discrepancy.

    KenFarley21-Topaz IIAuthor
    21-Topaz II
    March 29, 2014

    I got the original formula from the CRC math tables. A handy book from olden times which has all sorts of nice answers to unusual geometry questions. It's an odd conundrum. One would, reasonably, expect any sound mathematical method to yield the same results.

    I can see where your integration comes from, and how it gets the number Creo gets. But still, I wonder if I'm missing something crucial about the frustrum...