Need to solve the differential equation for beam deflection and get the following plot

11 次查看（过去 30 天）
Haris Hameed
Haris Hameed 2019-11-5
评论： Haris Hameed 2019-11-5
Hi all
i am trying to solve the beam defelction equation and get the plot (as shown in image),
can someone guide me further how to code all this, i developed some of it but cant proceede further
thanks for your time
clc
clear
for F=0:0.1:0.5 %force
%F=0;
d=0.6; %diameter mm
l=sqrt(72); %length mm
A=pi/4*d^2; %cross section dia
E=2460; %youngs modulus in MPa
mu=0.3;
G=E*2*(1+mu);
I=1/4*pi*(d/2)^4;
k=1/1.1; %shape factor
i=0;
w(1)=0;
th=35.264*pi/180;
a=F/(k^2*G*A);
b=F/(E*I);
th=32*pi/180;
for x=0:0.2:l
i=i+1;
%w(i)=-a*x+b*x^2/4*(2-l);
w(i)=-(a*x)+(b*x^3/6)-(b*x^2*l/4);
h(i)=w(i)*cos(th);
%w(1,43)=0;
end
x=0:0.2:l;
plot(x+h,w+h)
hold on
plot(-x+h,w+h)
end
  3 个评论
显示 1更早的评论
David Wilson
David Wilson 2019-11-5
Are you sure you have the correct equations? Typically for beam problems you have a 4th order differential equation (not 1st as you have written), and therefore you need 4 boundary conditions, and to solve it you need bvp4c or equivalent.

请先登录，再进行评论。

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

采纳的回答

Richard Brown
Richard Brown 2019-11-5
I think the problem you're having is not fully figuring out your solution before diving in and coding it. You have a differential equation of the form
where are constants. This equation is valid on the domain . It's easy enough to solve by integrating directly
You can figure out what is from the boundary condition (all the other terms go away when ). If you think of it this way and define the various constants in terms of the physical parameters, it's a bit easier to see what's going on - all of the parameters "get in the way" and make the expressions look overcomplicated.
So that's the mathematical part, which is pretty straightforward. The problem is relating it to your picture -- what are we looking at? What is the actual geometry of a problem? Where is the force being applied? The picture doesn't look like a cantilevered beam.
  1 个评论
Haris Hameed
Haris Hameed 2019-11-5
its like four beams joined together, force is applied at top, the deformed and undeformed shapes are shown, (dotted as undeformed). i thought to model one beam and then mirror it to get this shape. the beam deflects in y direction but how to apply (move) it in x direction.
prob2.jpg

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by