DMUG-Archiv 2000

Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

Problem mit NDSolve

Hallo,

ich habe mal wieder ein Problem: 
Ich moechte numerisch die DGL fuer eine 1-dimensionale 
Wellengleichung loesen.
Fuer t=0 ist das System in Ruhe:
                        y[x,0]==0,
                        Derivative[0,1][y][x, 0]==0,
Ein Ende sei fest:
                        y[1,t]==0,

Wenn das andere Ende mit
                        y[0,t] ==Sin[2 Pi t]
peridisch bewegt wird, gibt es keine Probleme:

In[1]:=
welle1=
        NDSolve[
                {D[y[x, t],t, t] == D[y[x, t], x, x],
                        y[x,0]==0,
                        y[1,t]==0,
                        Derivative[0,1][y][x, 0]==0,
                        y[0,t]==Sin[2 Pi t]
                        }, 
                y, {x, 0, 1}, {t, 0, 5}
        ]

Out[1]=
{{y\[Rule]InterpolatingFunction[{{0,1.},{0.,5.}},"<>"]}}

Wenn ich aber eine periodische KRAFT 
                        Derivative[0,2][y][0,t]==Sin[2 Pi t]
auf ein Ende wirken lasse, findet Mathematica keine Loesung

In[2]:=
welle2=
        NDSolve[
                {D[y[x, t],t, t] == D[y[x, t], x, x],
                        y[x,0]==0,
                        y[1,t]==0,
                        Derivative[0,1][y][x, 0]==0,
                        Derivative[0,2][y][0,t]==Sin[2 Pi t]
                        }, 
                y, {x, 0,1}, {t, 0, 5}
        ]

sondern es gibt folgende Fehlermeldung:

NDSolve::bcnorm : Boundary conditions 
  {MethodOfLinesDump`dv[1, 1](1, 0)[t, 0] == Sin[2 Pi t]} 
  include non-normal derivatives

Eigentlich sollte das doch zu loesen sein. Sehe ich nur den Wald
vor lauter Baeumen nicht, oder ist das tatsaechlich so ein 
schwieriges Problem?


Mit freundlichen Gruessen

        Jochen Skupin


Antworten:
Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

DMUG DMUG-Archiv, http://www.mathematica.ch/archiv.html