Finding a Variable through iterating. This one is killing me! 1

[maloneb]

Hey guys, I'm relatively new to mathcad and have hit a wall. All the data needed is attached to this message.

The formulae are lifted straight out of DNV-OS-F101 for the Local Buckling - Collapse of off shore Pipelines.

The Problem is that the data provided must be used to find a value of t_nom by setting the value of P_c = P_coll.

In Excel, the program runs with an estimated value of t_nom and iterates until the value of P_c > P_coll, therefore giving a suitable value of t_nom.

Is there any way of doing this in mathcad?

Thanks guys, I'm really stuck here.

 

[jghrist]

Generally, this type of problem can be solved using a solve block. Please save the workbook in version 11 or 13 and attach.

 

[maloneb]

Here is the workbook in version 11.

 

[IRstuff]

In order to get anything to work at all, your units must balance, which they don't. I've added what I think you're trying to do:

http://files.engineering.com/getfil...3dc-45af-86c6-5b00c4c85e68&file=Collapse2.mcd

TTFN

FAQ731-376
Chinese prisoner wins Nobel Peace Prize

 

[maloneb]

To IRstuff: I haven't altered the equation at all though, it's pulled straight out of a DNV document, so you would have thought they'd balance?!?!

In the document you've presented, t_nom is still undefined earlier in the equation, and Find(t_nom) is also undefined, how do I solve this?

I really am a complete newbie with mathcad, I'm talking baby steps...

 

[Occupant]

I have done a simple test and as you can see, it would be difficult for Matcad to find a solution. However, that is assuming I put the correct units in where you neglected to do so.

 

[jghrist]

To IRstuff: I haven't altered the equation at all though, it's pulled straight out of a DNV document, so you would have thought they'd balance?!?!

Often sources show equations, particularly empirically derived ones, with constants that require particular units. You either need to use those particular units and not explicitly show them in MathCAD or adjust the equations for the units used.

 

[IRstuff]

The OP has two variables that he uses that ostensibly should have units, E and sigma.y; assuming both have pascals as units allows the equations to be balanced. However, the RHS expression is extremely flat, while the other side is quite steep, making it difficult for the solver to find a solution.

As a general rule, one should graph the expressions to see whether they behave as expected, and whether a solution is plausibly achievable. The equations, as written have graphical solution ~0.5183226, which seems absurdly small. There also is an anomalous D.0/t.nom expression, where all the other instances are t.nom/D.0. However, changing that one instance makes little difference to the problem as a whole.

http://files.engineering.com/getfil...751-42be-8ed0-c874bcf5ac47&file=Collapse2.mcd

TTFN

FAQ731-376
Chinese prisoner wins Nobel Peace Prize

 

[jghrist]

0.518m doesn't seem absurdly small (it is meters, not mm) if this is pipeline wall thickness. If D.o = 350mm is the pipeline outside diameter, a 518mm wall thickness is impossibly large. If you use Minerr instead of Find, you get a solution 0.5183226*m.

 

[IRstuff]

You're right, I'd forgotten that the quickplot defaulted to meters, and yes, Minerr gets around the issues with Find.

TTFN

FAQ731-376
Chinese prisoner wins Nobel Peace Prize

 

[IRstuff]

File with Minerr attached:

http://files.engineering.com/getfil...3bd-4cab-85e8-32a4b893fc8d&file=Collapse3.mcd

TTFN

FAQ731-376
Chinese prisoner wins Nobel Peace Prize

 

[jghrist]

The problem might be the units. The modulus of elasticity of steel is about 2·10^11 Pa, not 207000Pa. If you use 2.07·10^11Pa, you get a t.nom = 5.16mm

 

[IRstuff]

The OP seems to have lost interest, but thanks for the new info; it's only just so much math for me, since it's well outside of my purview. Half a millimeter seems a bit thin, to me, but it feels much more plausible than half a meter.

TTFN

FAQ731-376
Chinese prisoner wins Nobel Peace Prize

 

[Occupant]

Just for the record, as IRstuff said, the OP has lost interest a long time ago: If you change the Modulus of Elasticity from whatever it was to steel then you probably have to look at Sigma.y too. If you use Pascal as unit, I think, it refers to soft butter.

 

[maloneb]

Hey guys, sorry, I didn't lose interest in this post, but I've been away for the bank holiday weekend. I'm just looking through everything now and will reply properly shortly.

 

[maloneb]

I didn't realise I'd missed the units for E and sigma.y, my bad. The units are MPa.

5.16mm is extremely close to the desired answer, a hell of a lot closer than I had been getting.

Just to clarify, is Minerr a pre-built in function to Mathcad?

Thanks for all the help guys

 

[Occupant]

Yes, Minerr is build in, but I did not use it. With the new units I came to 5.5mm - check it out.

 

[maloneb]

Occupant,

Do you have the file you have used? Could you post it please?

 

[Occupant]

Well, I didn't actually solve it, it is an approximation, based on the fact that the equation has to zero out - but you'll see.

 

[jghrist]

Actually, there are three solutions, -0.579mm, +0.569mm, and 5.51mm. Minerr gives them all, depending on your guess value.