Hi folks,
My problem is a simple one - calculate the vertical overburden stress in a soil mass. The simplified equation is just [ unit weight of soil ] * [ depth ]. However, in a layered system with a different unit weight for each layer, the calculation is best formed as the sum of the vertical stress increments. To me that looks like a stepwise function constructed as a series of IF-ELSEIF-ELSE statements as shown in my example below. It works but there is surely a more clever way.
Lets say the bottom of layer 1 is z1 and the stress is p1, and the stress at some mid-depth in layer 2, say z2, is delta_p2, such that the stress at the point in question is p(z2) = p1 + delta_p2. What is the most efficient way of structing that calculation as a function depth (z) so that it can be scaled to any number of layers. My best idea so far is an "If" function to evaluate the stress at the bottom of each layer plus the portion of the n-th layer at the target depth. Its cumbersome. In a 4 layer system, the function has 4 steps.
Obviously depth z and unit weight gamma are defined. They are each single column matrices. For depth, z_top is the depth to the top of each layer, and the unit weight is define similarly. In this case z is negative to indicate depth below ground surface.
As it is now, I suspect the program never ends, because it calls itself repeatedly, recursively.
I would have thought so to that this would happen if no local function is defined which is the case for z>z.top[0.
But both MC15 and P11 throw an error in this case:
The error message in MC15 ("Type "2" can not be loaded.") is a bit strange but both seem to refuse to refer recursively to the main/worksheet function if a first pass showed that a function with that name is defined in one of the if branches.
Both programs detect the infinite recursion in case no local function of that name is defined anywhere in the program
Its also remarkable that the first error message from MC15 came in German (the language set in Windows) even though my language settings for Mathcad are to be in English.
So I guess that the program made by @AV_14389071 works OK for him as long as he does not provide a z value larger than z.top_0.
But nonetheless I agree that he should replace the local sigma function definitions for a simple variable or simply do without and use the return statement (which also means that he could do without the if_else and just use a couple of ifs).
Speaking of recursion, we can use it to make the function definition a little bit smaller which can be especially useful if there would be even more layers:
Or even shorter using a for-loop which works for any number of layers
stacking minus infinity on top of ztop was necessary because Prime unfortunately does not use short curcuit evaluation of combined Boolean expressions.
If you want to get rid of the if/else statements and have a more flexible way of expanding the set of layers, you may want to try this:
Define:
and, just for the purpose of illustration:
then this function:
will give you gamma values for each value of z:
Of course you can replace/redefine gamma.s with any other vector of expressions, just make sure there are as many as there are comparisons.
Your example would be covered by:
There may be a more clever method to construct the vector with the gamma values, but it's late...
Success!
Luc
Here are some ways to achieve what you are looking for.
All three approaches work, regardless of the number of layers, provided that the number of elements in ztop and gamma is the same, and that the values in ztop are in descending order.
All three functions are defined in a way so they return a zero value in case z is larger than the first (largest) value in ztop.
I already posted the solution using recursion:
Here an approach inspired by Terry using linear interpolation
And here a version using the "match" function as suggested by Stuart
EDIT: Added a further variant. Probable more conventional, using a for-loop and taking advantage of the fact the program execution is stopped when a "return" statement is encountered (otherwise we had to use a while-loop).
Comparison of the various methods
Prime 11 worksheet attached
