Maximum of graph

Forum|Forum|10 years ago
August 31, 2015
3 replies
7674 views

In the attached prime worksheet I am graphing the shear and moment in a beam. Is there a simple (maybe 1 or 2 lines) for determining the minimum or maximum values on either graph? The three lines I have for determining the maximum is convoluted and inaccurate (the maximum is 8kip*ft, not 5.9 kip*ft).

1_306-beam-checks-mcdx.zip

Other

Best answer by Werner_E

Thanks for providing the picture.

Its a Prime specific problem and has to do with the different handling of units introduced in Prime. Unfortunately a lot of routines, especially numeric iteration algorithms like maximize, have not been adopted and tend to fail when used with units This was a problem in some cases in MC15, too, but not to that extent. Things got worse, even if the new unit handling is an improvement in general.

While numeric integration and numeric algorithms like the one used in maximize have its precision limits and can not be expected to get the 100% correct result 8 kip ft, the result given (approx. 6 kip ft) is not acceptable at all.

There are workarounds I could think of and depending on the situation they all have their benefits.

Best choice in your case seems to be to change the system variable TOL before the definition of the numeric integration.

You can do this globally in the calc ribbon or in the sheet (which makes for better doc) and it doesn't help to set it too low - results won't get better:

BTW, even though Mathcad seems not to care, you should not use a variable of integration as a limit in the same integral. Thtas the reason I changed you x to xi.

One other way to copy with the problem is to completety do without units, but this would mean to pass on one of the outstandig features of Mathcad.

A compromise is to get rid of the units temporarily by defining a unitless dummy function to maximize. You can of course additionally play around with TOL.

WE

R

RichardJ

19-Tanzanite

I can't read your worksheet because it's in Prime 3.1 format, but look up the maximize and minimize functions.

W

Werner_E

25-Diamond I

I can't read Prime 3.1 format either.

Do you graph a function or do you graph one vector against the other?

in the latter case you may look up the functions localmax and localmin.

L

LucMeekes

23-Emerald IV

You have functions to plot the data: here's the end result (it's a PITA to get Prime sheet data intio this forum). There's programming involved in getting Vd(x), so I can't evaluate due to Prime Express.

It should be possible to find the x-value of the extreme(s) for M(x) by finding the root(s) of V(x), because that's the derivative of M(x):

With a guess value xg:=0, find an extreme of M(x) with xM:=root(V(xg),xg), then the maximum of M(x) is found as M(xM).

If the root is not found, change the guess value.

For V(x) you can simply read the extremes from the graph.

To do it numerically you procede likewise as with M(x), but you need to find the roots of the derivative of V(x).

xg:=0 xV:=root(dd(V(xg),xg),xg). {Where dd(F(x),x) is the derivative operator to take the derivative of f(x) wrt x.}

Changing the guess value xg will allow to find other extremes.

Success!
Luc

W

Werner_E

Answer