How can I solve this question? And the ODE equation made me confused.

Question

Afifah Khairur 2020-6-29

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/556408-how-can-i-solve-this-question-and-the-ode-equation-made-me-confused

评论： Ameer Hamza 2020-6-29

WhatsApp Image 2020-06-29 at 00.26.08.jpeg

Hi, I have a problem regarding my assignment using MATLAB. This is my first time using this MATLAB. Can anybody help me with this problem ?

A simply supported beam of length is loaded by a distributed load and a tensile axial force as shown in the figure. The deflection of the beam, is determined from the solution of the following ODE:

(d^2 y)/(dx^2 )=1/EI [1+(dy/dx)^2 ]^(3/2) [1/6 q_0 (Lx-x^3/L)+Ty]

with the boundary conditions y(0)=0 and y(L)=0 .

EI=1.2×(10)^7 Nm^2 is the flexural rigidity, q_0= 30x(10)^3 N/m, and T=20x(10)^3 N .

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ameer Hamza 2020-6-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/556408-how-can-i-solve-this-question-and-the-ode-equation-made-me-confused#answer_458473

编辑：Ameer Hamza 2020-6-29

See bvp4c() for boundary condition problems. Something like this

EI = 1.3e7;
q0 = 30e3;
T = 20e3;
L = 4;
xmesh = linspace(0, L, 100);
odeFun = @(x, y) [y(2);
                  1./EI*(1+y(2).^2).^(3/2)*(1/6*q0*(L*x-x.^3/L)+T*y(1))];
              
bcFun = @(yl, yr) [yl(1)-0;  % y(0) = 0
                   yr(1)-0]; % y(L) = 0;
               
guess = bvpinit(xmesh, [0; 0]);
    
y_sol = bvp4c(odeFun, bcFun, guess);
plot(y_sol.x, y_sol.y);
legend({'$y$', '$\dot{y}$'}, 'Location', 'best', 'FontSize', 16, 'Interpreter', 'latex');