cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Help us improve the PTC Community by taking this short Community Survey! X

Problem solving second order ODE because of too many initial conditions

Sdruman
4-Participant

Problem solving second order ODE because of too many initial conditions

I am having issues solving a second order ODE because i need to force 3 boundary conditions, 2 for the starting point and one for the last one. I am forced to do this because otherwise the solution is incorrect so i can just obtain a solution with the correct derivative in 0 or the function with the correct value in 0. Of course if i enable all the three conditions mathcad says to me i have too many initial conditions but the problem is they are equally fundamental to obtain the correct solution. I will attach a picture also. Do you know how can i force a third condition solving a second order ODE? Or in general how to force much more conditions than the ODE order? thanks in advance.ode.png

1 ACCEPTED SOLUTION

Accepted Solutions
Sdruman
4-Participant
(To:Sdruman)

I solved the problem by myself.

I got the correct solution changing how large the interval is.

So basically changing it you obtain the only solution that is both balanced and with correct slips (it is a static problem). Basically i am changing the embedment length of the bar till i have the correct solution. You can do the same also changing the axial load at the beginning while length is fixed.

So basically if you have issues like these in future just try to change the boundary value or the length of the interval and look what it happens, you will see that odesolve solution change and you can find the one tha satisfy you. I would love to find a way to automatize this process using programming but seems like i am not able to put given and odesolve inside it. I will try maybe in future. Thanks everyone.

Immagine.png

View solution in original post

8 REPLIES 8
LucMeekes
23-Emerald III
(To:Sdruman)

Please attach the worksheet, not just a picture.

(For my convenience, save it as Mathcad 11 .mcd format; don't mind the warnings.)

 

Luc

Sdruman
4-Participant
(To:LucMeekes)

here you are

LucMeekes
23-Emerald III
(To:Sdruman)

I got this far:

LM_20190116_DEsolve.png

Success!

Luc

Sdruman
4-Participant
(To:LucMeekes)

So basically i will never be able to get the solution i want right? I mean that one that had both 0 value for derivative and funtion in zero plus the defined value at the end of the border. Because as i said i am forced in giving 3 conditions for this kind of problem otherwise the solution i got from ODE solve is not correct.


@Sdruman wrote:

So basically i will never be able to get the solution i want right?


I guess it means that the solution you are looking for simply does not exist!

Sdruman
4-Participant
(To:Sdruman)

I solved the problem by myself.

I got the correct solution changing how large the interval is.

So basically changing it you obtain the only solution that is both balanced and with correct slips (it is a static problem). Basically i am changing the embedment length of the bar till i have the correct solution. You can do the same also changing the axial load at the beginning while length is fixed.

So basically if you have issues like these in future just try to change the boundary value or the length of the interval and look what it happens, you will see that odesolve solution change and you can find the one tha satisfy you. I would love to find a way to automatize this process using programming but seems like i am not able to put given and odesolve inside it. I will try maybe in future. Thanks everyone.

Immagine.png

Being able to change the interval length of course changes the rules of the game and gives you one additional degree of freedom and so a solution exists.

Here is a way to let Mathcad find the optimal value for L.an. You have to turn the odesolve-block into a function of L.an and later use a normal solve block to find the value of L so, that sl(0)=0 (the condition not considered by odesolve.

 

B.png

Sdruman
4-Participant
(To:Werner_E)

Thanks my friend you saved my life 😄

Top Tags