Desperate solve 3rd order differentail equation

  • Thread starter Thread starter Hugh_Struct
  • Start date Start date
  • Tags Tags
    Desperate
Hugh_Struct
Messages
2
Reaction score
0

Homework Statement



Hi,

I am trying to complete my MSc in Structural Engineering and being an engineer my maths sometimes let's me down. Please help me to solve this 3rd order differential equation...it relates to solving the later torsional buckling of a beam. Here goes!

So far i am at this point

(3) GIt dr/dx - EIw d3r/dx3 = M2 / (PI2EIz / L2) dr/dx

I assume there is a homogeneous solution to this where i could use the following boundary conditions

(v)o = (v)L = 0

(r)o = (r)L = 0

(d2r/dx2)0 = (d2r/dx2) = 0

Homework Equations



The problem begins with these two equations...

(1) EIz d2v/dv2 = -M r(X)

and

(2) GIt dr/dx - EIw d3r/dx3 = M dv/dx

Trial solutions of

(4) v(x) = w sin (PI x /L) and r(x) = r sin (PI x /L)

so differentiating (4) you can obtain

(5) v(x) = M / (PI2EIz / L2) r(x)

The Attempt at a Solution



Apparently substituting equation (5) into equation (2) and using the boundary conditions you can resolve the problem to find
3...below is the solution
M = Square root [ (PI2EIz / L2) * (GIt + PI2EIw / L2)]


I say apparently as i can only manage the substation part! I know you could also differentiate equation (2) to a 4th order differential and then substitute d^2v/dv^2 of equation (1) but i really need to know the other method.

I am SO sorry about the awful text...first time poster...long time reader!

Please help!

Thanks

Hugh

Much better to check this link out...page 2 and half of page 3 is what I'm after!

https://noppa.tkk.fi/noppa/kurssi/rak-54.3600/luennot/Rak-54_3600_lateral-torsional_buckling.pdf"
 
Last edited by a moderator:
Physics news on Phys.org
Hugh_Struct said:

Homework Statement



Hi,

I am trying to complete my MSc in Structural Engineering and being an engineer my maths sometimes let's me down. Please help me to solve this 3rd order differential equation...it relates to solving the later torsional buckling of a beam. Here goes!

So far i am at this point

(3) GIt dr/dx - EIw d3r/dx3 = M2 / (PI2EIz / L2) dr/dx

I assume there is a homogeneous solution to this where i could use the following boundary conditions

(v)o = (v)L = 0

(r)o = (r)L = 0

(d2r/dx2)0 = (d2r/dx2) = 0



Homework Equations



The problem begins with these two equations...

(1) EIz d2v/dv2 = -M r(X)

and

(2) GIt dr/dx - EIw d3r/dx3 = M dv/dx

Trial solutions of

(4) v(x) = w sin (PI x /L) and r(x) = r sin (PI x /L)

so differentiating (4) you can obtain

(5) v(x) = M / (PI2EIz / L2) r(x)


The Attempt at a Solution



Apparently substituting equation (5) into equation (2) and using the boundary conditions you can resolve the problem to find
3...below is the solution
M = Square root [ (PI2EIz / L2) * (GIt + PI2EIw / L2)]


I say apparently as i can only manage the substation part! I know you could also differentiate equation (2) to a 4th order differential and then substitute d^2v/dv^2 of equation (1) but i really need to know the other method.

I am SO sorry about the awful text...first time poster...long time reader!

Please help!

Thanks

Hugh

Much better to check this link out...page 2 and half of page 3 is what I'm after!

https://noppa.tkk.fi/noppa/kurssi/rak-54.3600/luennot/Rak-54_3600_lateral-torsional_buckling.pdf"


If GIt, EIw, M2 / (PI2EIz / L2) are constants, your equation is of the form A d3r/dx3 = B dr/dx, where A = GIt and B = GIt - M2 / (PI2EIz / L2) are constants. Letting y = dr/dx, you have a second-order ODE A d^2 y/dx^2 = B y, which has a sinusoidal solution if B > 0 (assuming A > 0) or an exponential solution if B < 0.

RGV
 
Last edited by a moderator:
Thanks RGV,

That makes sense I'm trying to use the 2nd O.D.E A d2 y/dx2 = B y together with the B.C's to solve the equation. I will post it if i can get it...I may be a while!

Yes all the other terms are constants.

Hugh
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top