Skip to main content
1-Visitor
July 11, 2017
Solved

Problem with Odesolve in good Old MC15

  • July 11, 2017
  • 2 replies
  • 3628 views

Hi,

 

I want to calculate a deflection of a beam with Odesolve in MC15.

 

The deflection consists of two parts: y1(x) and y2(x).

 

I could solve the y1(x), but when trying to solve y2(x), i got an error-message:

"this variable is not defined"

You can see it near the end of the sheet.

11-07-2017 19-03-41.jpg

The boundary-conditions are: deflection in the bearings are zero, the bending moments as well.

 

What am i doing wrong here?

Thank you for your help.

 

ps.: The MC15 File can be uploaded only as a zipped file.

Best answer by Werner_E

You can't use the prime symbol (Ableitungsstrich) for the derivative of y1.

I would suggest to define a function for the derivative.Pic1.png

Pic2.png

 But I don't think that the result is correct. When I look at what the value of the first derivative of y1 at lo (6) is, I sometimes get an error (which I would have suspected at the end of the interval)Pic3.pngand sometimes it will (upon recalculation) evaluate to 0.Pic4.png

So I see two options, Either let the first odesolve block up to lo+0.1 (or something like lo*1.01) or you change the initial condition in your second solove block to y2*(lo)=y1_diff(lo-0.1) or something similar.

This gives a more plausible result:Pic5.png

 

 

 

 

2 replies

Werner_E25-Diamond IAnswer
25-Diamond I
July 11, 2017

You can't use the prime symbol (Ableitungsstrich) for the derivative of y1.

I would suggest to define a function for the derivative.Pic1.png

Pic2.png

 But I don't think that the result is correct. When I look at what the value of the first derivative of y1 at lo (6) is, I sometimes get an error (which I would have suspected at the end of the interval)Pic3.pngand sometimes it will (upon recalculation) evaluate to 0.Pic4.png

So I see two options, Either let the first odesolve block up to lo+0.1 (or something like lo*1.01) or you change the initial condition in your second solove block to y2*(lo)=y1_diff(lo-0.1) or something similar.

This gives a more plausible result:Pic5.png

 

 

 

 

vlehner1-VisitorAuthor
1-Visitor
July 11, 2017

Vielen Dank Werner,

Ich glaube das war der Knackpunkt mit dem definieren der Funktion für die erste Ableitung.

Auf die Schreibweise mit dem Ableitungsstrich wäre ich so schnell nicht gekommen.

 

Jetzt gilt es die Übergangsbedingungen sowie die Grenzen richtig zu finden/festzulegen.

 

Way to go!

vlehner1-VisitorAuthor
1-Visitor
July 11, 2017

...now it's done 🙂

THX Werner!

 

12-07-2017 00-43-40.jpg

23-Emerald I
July 12, 2017

Mathcad will have issues with Odesolve when the constraints aren't "initial conditions."

 

Attached is a finite difference solution.  You need a stronger beam!

vlehner1-VisitorAuthor
1-Visitor
July 12, 2017

Thanks Fred!