Inputs zu Lösungen

1

Literaturstudium

2

a Interpretation der Parameter: Angelegenheit der Physik oder der technischen Mechanik. Wichtig: Bei einer eindimensional gereichteten Bewegung kann z.B. m eine Masse bedeuten, d eine Dämpfungskonstante und k eine Federkonstante. Bei einer Drehbewegung wäre entsprechend Statt m das Massenträgheitsmoment I oder J zu setzen etc.
D wäre im ersten Fall eine auf die Feder und die Masse bezogenes Dämpfungsmass und (Subscript[ω, D])eine Kreisfrequenz, welche in der Lösung erscheint.

b

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[m y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->y0,y'[0]->b}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

(* Allgemeine Lösung bei freier Schwingung! *)

c

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[m y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->b}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

(* Allgemeine Lösung bei freier Schwingung! *)

d

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+1/2 y'[t]+1 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->1}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

Plot[u0[t],{t,0,15}];

e

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+1/2 y'[t]-1 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->1}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

u0[t]//Expand//N

Plot[u0[t],{t,0,10}];

(* ===> Explosion *)

f

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[m y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[DiracDelta[t] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

g

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+1/2 y'[t]+1 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[DiracDelta[t] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

U[s]//Apart

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

Plot[u0[t],{t,0,15}];

h

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[m y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[A Sin[ω t+ φ] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

i

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[m y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[A Sin[ω t] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

u0[t]//Expand

j

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[A Sin[ω t] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

k

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+d y'[t]+k y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->0,y'[0]->0,d->2 Sqrt[k]}

2. Rechte Seite transformieren

rechts=LaplaceTransform[A Sin[ω t] ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]/.d->2 Sqrt[k]

l, m

solv=Solve[ω^4+(d^2-2 k)ω^2+k^2==0,{ω}]//Simplify //Flatten

solv[[3]]

ωϵ1[ϵ_] := ω + ϵ /. 
        solv[[1]];  ωϵ1[ϵ]

ωϵ2[ϵ_]:=ω+ϵ/.solv[[2]]; ωϵ2[ϵ]

(* Fall mit den meisten positiven Anteilen unter Wurzeln *)

ωϵ2[0]//FullSimplify

ωϵ2[0]//FullSimplify//InputForm

(* (d*(-d + Sqrt[d^2 - 4*k]))/2 ist für positive k negativ *)

Solve[(d*(-d + Sqrt[d^2 - 4*k]))/2 + k == 0,{k}]

(* ωϵ2[0] ist daher nie positiv *)

ωϵ2[0]/Sqrt[d^2]//PowerExpand

ωϵ3[ϵ_]:=ω+ϵ/.solv[[3]]; ωϵ3[ϵ]

ωϵ4[ϵ_]:=ω+ϵ/.solv[[4]]; ωϵ4[ϵ]

n

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+ 2q y'[t]+ω0^2 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->1,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]/.d->2 Sqrt[k]

o1

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+ 2q y'[t]+ω0^2 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->1,y'[0]->0,q->1, ω0->1/2}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

Plot[u0[t],{t,0,10}];

o2

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+ 2q y'[t]+ω0^2 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->1,y'[0]->0,q->1, ω0->1}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

Plot[u0[t],{t,0,10}];

o3

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+ 2q y'[t]+ω0^2 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->1,y'[0]->0,q->1/2, ω0->1}

2. Rechte Seite transformieren

rechts=LaplaceTransform[0 ,t,s]

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U1[s]:=Y[s]/. solv;  U1[s]

U[s]:=(U1[s]//Apart)/.Sign[ω]->1

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

Plot[u0[t],{t,0,10}];

p

Eigene Experimente....

q

Remove["Global`*"]

1. Linke Seite transformieren, Anfangswerte anpassen

links = LaplaceTransform[1 y''[t]+ 2*1 y'[t]+2^2 y[t],t,s] /. {LaplaceTransform[y[t],t,s]->Y[s],y[0]->1,y'[0]->0}

2. Rechte Seite transformieren

rechts=LaplaceTransform[Sin[ω t] ,t,s]/.Sign[ω]->1

3. Gleichung links = rechts lösen

solv=Solve[links==rechts,{Y[s]}] // Flatten

4. Rücktransformation

U[s]:=Y[s]/. solv;  U[s]

u0[t_]:=InverseLaplaceTransform[U[s],s,t]//Simplify; u0[t]

5. Konstant schwingender Term isolieren, herausrechnen, Plot der Amplitude als Funtion der Frequenz

col1=Collect[Evaluate[u0[t]/.{Sin[ω t]->x,Sin[Sqrt[3] t]->x^2,Cos[ω t]->x^3,Cos[Sqrt[3] t]->x^4}],{x,x^2,x^3,x^4}];
col1

col2 =((col1)/.{x->Sin[ω t],x^2->Sin[Sqrt[3] t],x^3->Cos[ω t],x^4->Cos[Sqrt[3] t], Sqrt[ω^2]->ω})// Cancel

col3 = col2/. {E^(-t)->0, Sqrt[ω^2]->ω}

col4 = col3/. {Cos[ω t]->1, Sin[t ω]->0}

col5 = col3/. {Cos[ω t]->0, Sin[t ω]->1}

cω1 = col4/Sqrt[col4^2+col5^2]//Simplify

cω2 = col5/Sqrt[col4^2+col5^2]//Simplify

schwing = -(Sqrt[col4^2+col5^2](Sin[α] Cos[ω t]+Cos[α] Sin[t ω])//Simplify)/. α->ArcSin[cω1]

aml=-Sqrt[col4^2+col5^2]//Simplify

aml//InputForm

amplitude[ω_]:=1/Sqrt[16 - 4*ω^2 + ω^4]

Plot[amplitude[ω],{ω,0,10}];

Solve[D[amplitude[x],x]==0,{x}]

N[%]

6. Bemerkung: Die Amplitudenformel gibt dasselbe Resultat:

Sqrt[2^2-2*1^2]//N