Message Boards Message Boards

0
|
17229 Views
|
4 Replies
|
2 Total Likes
View groups...
Share
Share this post:

How to solve a higher order partial differential equation with boundaries

Posted 10 years ago
Attachments:
POSTED BY: ihsan faseel
4 Replies

Looks difficult, even setting all constants equal to 1 one gets

In[8]:= DSolve[{D[y[x, t], {t, 2}] + D[y[x, t], {x, 4}] + D[y[x, t], {x, 4}, t] == 0,
  (y[x, t] /. x -> 0) == 0, (D[y[x, t], x] /. x -> 0) == 0, (D[y[x, t], {x, 2}] /. x -> L) == 0, 
  (D[y[x, t], {x, 3}] /. x -> L) == 0}, y[x, t], {x, t}]
Out[8]= (* actually the input *)

so give L a value and try NDSolve:

In[9]:= With[{L = 5},
 NDSolve[{D[y[x, t], {t, 2}] + D[y[x, t], {x, 4}] + D[y[x, t], {x, 4}, t] == 0,
   (y[x, t] /. x -> 0) == 0, (D[y[x, t], x] /. x -> 0) == 0, (D[y[x, t], {x, 2}] /. x -> L) == 0, 
   (D[y[x, t], {x, 3}] /. x -> L) == 0}, y[x, t], {x, t}]
 ]

During evaluation of In[9]:= NDSolve::underdet: There are more dependent variables, {y[x,t],(y^(0,1))[x,t],(y^(0,2))[x,t]}, than equations, so the system is underdetermined. >>
Out[9]= (* actually the input *)

first job you should work on is to get a numerical solution with unit constants. Seemingly the problem is incorrect formulated.

POSTED BY: Udo Krause
Posted 10 years ago

sir, actually "E","I" ,"P","A","C" are constants. like : E-Youngs modulus. I-moment of inertia. A-area. C-damping coefficient.

So what should i do sir.Should i define all these with there values.Actually i want to get a solution containing these constants as itself.

POSTED BY: ihsan faseel
Posted 10 years ago

"EI" is not defined, as "PA" and "CI" are not defined either.

POSTED BY: Simon Tyran

It appears that there are more unknowns than equations.

POSTED BY: S M Blinder
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract