Well, what you have here is a 3rd order equation in 3 variables. Not a pretty thing. A mathemation might have fun with it, but an engineer's approach is how can I get the correct answer in a reasonable way. (Apologies to my old math professors, if any are still around).
Looking at your problem definition, you really don't have 3 variables. You are using the golden ratio as a variable, but it is a constant -- a known number and thus it's square is also a known number. That gets you to a 3rd order equation in 2 variables, closer but still no cigar. If you'd allow that in a given situation you would know the thickness of available materials, you could make t a constant - and solve individually for different values of t. That's what I would do (actually, I might just set up a table in excel, which would be quicker than solving the equation).
But, for the solution to 3 decimal points:
Golden ratio = 1.618; squared = 2.618
a = shortest side
b= 1.618 *a
c= 2.618*a
t= thickness of material
V=(a-2t) * (1.618a - 2t) * (2.618a - 2t)
multiplying (first, the first two)
V= (1.618 a^2 - 5.236 at + 4T^2) * (2.618a - 2t)
finally:
V= 4.236 a^3 - 16.944 a^2*t + 20.944 at^2 -8t^3
I left out the combining terms steps (1 in each multiplication) because I assume you know how to do that.
So, now you just pick a value for the thickness of your material, adjust the constants for that value and you have a solvable cubic equation. You can then plug whatever value you want for V and solve for a.The solution to the cubic equation will give the 3 possible mathematical values for a, but only one will make sense.
You then calculate b and c by applying the golden ratio and its square. To check your answer, subtract double the thickness from each and multiply the three results together to check that this equals your target volume.
A couple of examples:
Let t=1
then
V= 4.236 a^3 - 16.944a^2 + 20.944a -8
pick a value for V and solve for a:
For V=100
a=4.249 ----- minus 2, equals 2.249
b=4.249*1.618=6.875 -------minus 2, equals 4.875
c=4.249*2.618=11.124-------minus 2, equals 9.124
Actual V = 2.249*4.875*9.124 = 100.034
If that isn't accurate enough for you, use more decimal points. This should work for all values of V
Another example:
Let t=.5
then
V= 4.236a^3 - 8.472 a^2 +5.236 a -1
For V=1,000
a=6.852 ----- minus 1, equals 5.852
b=6.852*1.618=11.087-------minus 1, equals 10.087
c=6.852*2.618=17.939-------minus 1, equals 16.939
Actual V = 5.852*10.087*16.939 = 999.893
|