Trying to solve a system of 2 ordinary differential equations with boundary conditions
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi there, I am trying to solve a system of ordinary differential equations with boundary conditions. These equations are part of my thesis and govern thick cylinders under internal pressure. I have checked my equations with the general theroy as suggested by Gabriel Lame for thick walled cylinders. The codes simplify as follows:
A,B,C and D are the functions of E and poisson's ratio.
the boundary conditions are t=1, y=0 and at t=R, y=-P. (As in the case of cylinder with internal pressure only)
I have tried to solve it in matlab but get an [ empty sym ] in return. I have limited experience with matlab so would appreciate any help possible.
My code:
syms E nu a11= (1/E); a12= (-nu/E); a33= (1/E); A= (a12)/(a11+a12); B= ((a33)-((2*a12^2)/(a11+a12))); C= (a11)/(a11^2-a12^2); D= (2*a12+a11)/(a11+a12);
syms x(t) y(t) eqns = [diff(x,t)==A*(x/t)+B*y, diff(y,t)==-C*(x/t^2)-D*y]; cond = [y(0) == 0, y(2)==1]; withSimplifications = dsolve(eqns, cond)
0 Kommentare
Antworten (4)
Torsten
am 22 Nov. 2016
I doubt that your system can be solved analytically using dsolve.
Try a numerical solution with bvp4c.
If you start at t=0, problems will arise because t is in the denominator in both equations.
Furthermore, are you sure that you don't have to prescribe any boundary condition for x ?
Best wishes
Torsten.
Ketav Majumdar
am 23 Nov. 2016
1 Kommentar
Karan Gill
am 23 Nov. 2016
Is this example in the documentation helpful? https://www.mathworks.com/help/symbolic/dsolve.html#bve5nw0-1_1
Siehe auch
Kategorien
Mehr zu Symbolic Math Toolbox finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!