0
|
11459 Views
|
6 Replies
|
3 Total Likes
View groups...
Share
GROUPS:

# Particular Solution of 4th degree differential equation (column buckling)

Posted 9 years ago
 I'm trying to use Mathematica to solve the Euler Column buckling problem (http://www.continuummechanics.org/cm/columnbuckling.html), but its not working out. Mathematica only shows me the zero solution. Is there any way of getting the no zero displacements solutions? For now, what I have is basically this: equation = EIy''[x] + Py[x] == 0 DSolve[{equation, y == 0, y[L] == 0}, y[x], x] Out= {{y[x] -> 0}} Thanks, Luís Valarinho
6 Replies
Sort By:
Posted 9 years ago
 Hello Kay!Thanks for the suggestion. What I did actually was using the nomenclature that we use in Civil Engineering But you're right, it can lead to mistakes. Thanks for the other suggestions too. I was actually wondering if it is possible to do it in one step. We cannot ask Mathematica to give all the possible solutions? Or adding an assumption that C should be different than zero? (something that is not necessarily truth).Thanks, Luís Valarinho
Posted 9 years ago
 I recommend not using capital letters for your constants (or at least give them unique names) because for example I is the symbol for Sqrt[-1] and E is the euler number Here is how I would solve it: In:= Clear["Global*"] In:= sol = DSolve[{y''[x] + ps/(es is) y[x] == 0, y == 0}, y[x], x] Out= {{y[x] -> C Sin[(Sqrt[ps] x)/(Sqrt[es] Sqrt[is])]}} In:= Solve[Sqrt[ps] l/Sqrt[es is] == n Pi, ps] Out= {{ps -> (es is n^2 \[Pi]^2)/l^2}} `
Posted 9 years ago
 Ok, thanks!
Posted 9 years ago
 The solution given is correct but it isn't the only solution. Yes, you'll have to handle the boundary condition at L separately.
Posted 9 years ago
 So, if I understood Mathematica doesn't solve the equation automatically. The last boundary condition will have to be me adding it, right?
Posted 9 years ago
 Replace y[L] ==0 with y' == yp0 and then combine y[L] ==0 with the solution you get.