This is a clarified follow-up to a previously solved question. The new objective is to compute the soil wedge area, weight, and centroid (CGx, CGy) behind a retaining wall using the mathematical approach described in the attached Mathcad files, and input data is already explained in the attachment. I am using Mathcad Prime 11. Any help would be appreciated.
Thanks @terryhendicott ,
I very much appreciate sharing your solutions, which shows that you are familiar with the problem I need to resolve. The example you shared is still does not give what I would need, and I do not know how to tweak it!! My MathCAD skills are still not as good as you may think!
I will need to calculate the area(s) and centroid(s) (CGx, CGy) of each shape bounded by the vertical plane and the inclined planes bounded by the dashed lines incrementally. The challenges to do this is to identify the coordinates of the intersecting points with the soil, then calculate the area(s), whether they are triangles or polygon, and their centroids (x & y) for each shape and so on.
This is the "trial wedge method". I appreciate sharing your problem, but I am not sure how to tailor it to get it to do the trial wedge.
For a given angle ß (running from 90° down to 34°), which polygon are you looking for to calculate area and centroid?
The left or the right in this picture:
Is there a guarantee that the first segment in xy has an intersection with the vertical line at B1?
This means that the second x-value in xy is always larger than B1.
EDIT: Your reply to Terry seems to make clear that you have the left version in mind where the vertical line at B1 is a border.
The last question is not that important as we can run through xy until we come across a point to the right of B1,
P.S.: I think it was a mistake that you marked this post as the correct answer. You should remove that mark as others may ignore a thread they think is already answered.
Thanks @werner ...
1) The left.
2) I am not clear on the second question, regarding intersection with B1, but I will answer both scenarios I thought you might mean
a) If I assume that you mean the vertical line, which is referenced as "Pressure plane", this is where the first segment starts, so the vertical line is one side of the segment. The intersection I mean will be with the soil profile which is the topmost line.
b) If you mean β₁ instead (the angle), what I am showing is read as βᵢ "with subscript letter i not 1"
as the iteration progresses, the wedge continuously increases in size, as you show on the left, and each new wedge has its own updated centroid.
The angle βᵢ is not intended to represent a single fixed plane, as I shown (I apologize!). Rather, it is a dynamic parameter representing each incremental dashed plane (β₁, β₂, β₃, β₄, and so on).
The dashed lines I am showing at these locations are just schematic for illustration purposes only. In practice, the dashed planes represent approximately 1,000 individual planes, each associated with a unique corresponding β angle value, starting from the vertical line (pressure plane), and rotating until it with the solid green line. I understand it was confusing,
As shown in my file, the angle β is generated using an equation that computes about 1,000 incremental angles ranging from 90° to 34°, with each step equal to 1/1000 of the total angular range.
I realize the dashed lines I added, for illustration, may have caused some confusion, my apologies!. The blue solid line shown happens to align with one of the planes in my model, but the dashed planes are not fixed at those locations. They are virtual and schematic only, generated dynamically during the iteration process.
I’ve also added a couple of notes to the original sketch for clarity. (attached). Please let me know if this addresses your questions, and I truly appreciate your help.
What I meant with my second question was if the soil polygon could also look that way:
Here the vertical line at x=B which is the left border of the polygon of interest does not intersect the first segment of the soil polygon.
But it doesn't matter anymore. The solution I am about to post should be able to deal with this situation as well.
OK, based on the informations so far here is my attempt.
I provide a function "calculate" which, given the angle beta and the soil polygon as input returns area, center of gravity, weight and also the polygon so we can plot it. Its not completely self-contained as it relies on variables B1 and gamma to be defined in the worksheet.
Of course this function "calculate" can also be used to create a table of values. I have chosen 1120 steps because this will make for a nice step value of 0.05 degree but of course that value can be freely chosen.
Because Prime does not allow to change the units in a matrix with mixed units i decided to create a table with dimensionless values and add a header line which specifies the units.
Hope this helps
Prime 11 sheet attached
Here is the modified worksheet.
I also added a way to create single vectors for all the desired quantities instead of a compact table as it seems that's what you actually are looking for.
You may have to check the logic in your function lsurLL because this function returns 0 feet for any value xb provided - at least when used with LL index 1 and 3. Or is this intentional?
Always check your functions by evaluating them using some test data!
You may also test them with a vector of sample values if you use vectorization:
