1.2.G
Gegeben
Symbole:
Größen:
Die Kraft \(\boldsymbol R\) wirkt am Teilchen P.
Zerlegen Sie die Kraft \(\boldsymbol R\) in:
die Kraft \(\boldsymbol F_1\) entlang der Richtung \(\boldsymbol e_1\)
die Kraft \(\boldsymbol F_2\) entlang der Richtung \(\boldsymbol e_2\) , sodass:
a) Vektor-Komponenten
Sei:
\(F_1\): Komponente von \(\boldsymbol F_1\) in Richtung des Einheitsvektors \(\boldsymbol e_1\), sodass \(\boldsymbol F_1 = F_1 \boldsymbol e_1\)
\(F_2\): Komponente von \(\boldsymbol F_2\) in Richtung des Einheitsvektors \(\boldsymbol e_2\), sodass \(\boldsymbol F_2 = F_2 \boldsymbol e_2\)
Die \((x,y)\)-Komponenten der Kräfte sind dann:
mit den Abkürzungen:
Stellen Sie damit und aus der Bedingung, dass \(\boldsymbol F_1 + \boldsymbol F_2 = \boldsymbol R\) zwei Gleichungen für \(F_1\) und \(F_2\) auf. Und berechnen Sie dann die zwei Komponenten \(F_1\) und \(F_2\) aus diesen zwei Gleichungen. Und zwar:
für die gegebenen Symbole und
für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\).
Dies ist ein lineares Gleichungssystem für \(F_1\) und \(F_2\). Lösung: Der Nenner darf nicht gleich null sein. Die zwei Gleichungen in einer Matrix-Gleichung zusammengefasst: Berechnung der Inversen der Matrix \(A\) nach dieser Formel: Lösung für \(F_1\) und \(F_2\): Einsetzen der gegebenen Größen führt auf:Lösung
Lösungsweg
b) Grüne Komponenten
Berechnen Sie die \((x,y)\)-Komponenten der Kräfte , und zwar für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\). Zeigen Sie, dass:
Einsetzen der gegebenen Größen führt auf das, was zu zeigen ist.Lösung
c) Vektor-Addition
Prüfen Sie Ihr Ergebnis zeichnerisch mit einem Kräfteparallelogramm, in das Sie \(\boldsymbol F_1\) und \(\boldsymbol F_2\) und \(\boldsymbol R\) einzeichnen.
Source Code eines Programms dazu: Programm ausführen? Dazu: Copy: Source Code in die Zwischenablage kopieren. Paste: Source Code ins Eingabefeld hinter Play: Knopf SymPy
from sympy.physics.units import kg, m, s
from sympy import pi, S, pprint, symbols, cos, sin, Eq, solve, Matrix
# Units:
N = kg*m/s**2
k, M, G = 10**3, 10**6, 10**9
kN = k*N
Pa = N/m**2
MPa, GPa = M*Pa, G*Pa
deg = pi/180
half = S(1)/2
# Rounding:
def iso_round(x, pv):
try:
x = float(x)
tmp = round(x, pv)
except TypeError:
x = x.evalf()
from functools import partial
func = partial(round, ndigits = pv)
tmp = x.applyfunc(func)
return tmp
R, p1, p2 = symbols("R, p1, p2", real=True)
prec = 2
sub_list = [
( R, S(9)/2*N ) ,
( p1, 45 *deg ) ,
( p2, 15 *deg ) ,
]
# Trig:
C1, S1 = cos(p1), sin(p1)
C2, S2 = cos(p2), sin(p2)
# Shortcuts:
c1, c2, s1, s2 = symbols("c1, c2, s1, s2")
sub_list_shortcuts = [
(C1, c1),
(S1, s1),
(C2, c2),
(S2, s2),
]
F1, F2 = symbols("F1, F2", real=True)
unknowns = [F1, F2]
eq1 = Eq(F1*C1 + F2*C2, 0)
eq2 = Eq(F1*S1 + F2*S2, -R)
eqns = [eq1, eq2]
sol = solve(eqns, unknowns)
F1, F2 = sol[F1], sol[F2]
pprint("\nF1:")
tmp = F1
pprint(tmp.subs(sub_list_shortcuts))
pprint("\nF1 / N:")
tmp = tmp.subs(sub_list)
tmp /= N
tmp = iso_round(tmp, prec)
pprint(tmp)
pprint("\nF2:")
tmp = F2
pprint(tmp.subs(sub_list_shortcuts))
pprint("\nF2 / N:")
tmp = tmp.subs(sub_list)
tmp /= N
tmp = iso_round(tmp, prec)
pprint(tmp)
F1v = Matrix([F1*C1, F1*S1])
F2v = Matrix([F2*C2, F2*S2])
pprint("\n(F1x, F1y):")
tmp = F1v
pprint(tmp.subs(sub_list_shortcuts))
pprint("\n(F1x, F1y) / N:")
tmp = tmp.subs(sub_list)
tmp /= N
tmp = iso_round(tmp, prec)
pprint(tmp)
pprint("\n(F2x, F2y):")
tmp = F2v
pprint(tmp.subs(sub_list_shortcuts))
pprint("\n(F2x, F2y) / N:")
tmp = tmp.subs(sub_list)
tmp /= N
tmp = iso_round(tmp, prec)
pprint(tmp)
[ ]: einfügen.▶ drücken.