DMUG-Archiv 2011

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

Re: Probleme mit NDSolve

Hallo Alexander,
eine Möglichkeit ist:

Werte = {m1 -> 2, m2 -> 2, g -> 9.81, l0 -> 49, d -> 10000, 
    \[Alpha] -> 30.21}; 
m = \[Phi] = 
    d/E^(2*\[Alpha]*(Sqrt[(x22[t] - x21[t])^2 + (y22[t] - y21[t])^
            2] - l0)) - 
     (2*d)/E^(\[Alpha]*(Sqrt[(x22[t] - x21[t])^2 + (y22[t] - y21[t])^
            2] - l0)); 
Fsx1 = D[\[Phi], x21[t]] /. Werte; 
Fsx2 = D[\[Phi], x22[t]] /. Werte; 
Fsy1 = D[\[Phi], y21[t]] /. Werte; 
Fsy2 = D[\[Phi], y22[t]] /. Werte; 
s = NDSolve @@ ({{Derivative[2][x21][t] == Fsx1/m1, x21[0] == 0, 
       Derivative[1][x21][0] == 30, m2*Derivative[2][x22][t] == 
        Fsx2, x22[0] == 50./Sqrt[2], Derivative[1][x22][0] == 0, 
       m1*Derivative[2][y21][t] == (-m1)*g + Fsy1, y21[0] == 0, 
       Derivative[1][y21][0] == 70, m2*Derivative[2][y22][t] == 
        (-m2)*g + Fsy2, y22[0] == 50/Sqrt[2], 
       Derivative[1][y22][0] == 0}, {x21[t], x22[t], y21[t], 
       y22[t]}, {t, 0, 10}} /. Werte); 
ParametricPlot[Evaluate[{{x21[t], y21[t]} /. s, 
    {x22[t], y22[t]} /. s}], {t, 0, 0.5}]

Gruss,

Rolf

On Sat, 15 Jan 2011 17:31:28 +0100, Alexander Knöttner
<alex.knoettner@XXXXXXX.at> wrote:
> Hallo Udo.
> 
> Tut mir Leid, auf diese Idee hätte ich natürlich kommen können, hätte
> auch
> ich mir das mühsame einfügen ersparen können.
> 
> Freut mich, dass sich jemand für mein Problem interessiert, Danke!
> 
> Liebe Grüße Alex
> 
> -----Ursprüngliche Nachricht-----
> Von: Udo und Susanne Krause [mailto:su.krause@XXXXXXX.ch] 
> Gesendet: Samstag, 15. Jänner 2011 17:23
> An: demug@XXXXXXX.ch; Alexander Knöttner
> Betreff: Re: Probleme mit NDSolve
> 
> Hallo Alexander,
> 
> könnten Sie nicht freundlicherweise das Notebook anhängen (Dateianhang 

> auswählen und die entsprechende Notebookdatei auswählen), dann müsste
> man  
> nicht alles abtippen am Sonnabend ...
> 
> Gruss
> Udo.
> 
> On Sat, 15 Jan 2011 15:11:25 +0100, Alexander Knöttner  
> <alex.knoettner@XXXXXXX.at> wrote:
> 
>> Werte={m1R2,m2R2,gR9.81,l0R40,dR10000,aR30.21}
>> m=f=d*-2*d*
>> Fsx1=D[f,x21[t]]/.Werte
>> Fsx2=D[f,x22[t]]/.Werte
>> Fsy1=D[f,y21[t]]/.Werte
>> Fsy2=D[f,y22[t]]/.Werte
>>
>
s=NDSolve[{x21''[t]==Fsx1/m1,x21[0]==0,x21'[0]Š30,m2*x22''[t]==Fsx2,x22[0]==
>
50/,x22'[0]==0,m1*y21''[t]==-m1*g+Fsy1,y21[0]==0,y21'[0]Š70,m2*y22''[t]==-m2
>
*g+Fsy2,y22[0]Š50/,y22'[0]==0},{x21[t],x22[t],y21[t],y22[t]},{t,0,10}]/.Wert
> e
>>
>
ParametricPlot[Evaluate[{{x21[t],y21[t]}/.s,{x22[t],y22[t]}/.s}],{t,0,0.5}]



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

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