1.3.E
Gegeben
Symbole:
Größen:
a) Kraftvektoren
Berechnen Sie die \((x,y,z)\)-Komponenten von \(\boldsymbol F_1\) und \(\boldsymbol F_2\)
für die gegebenen Symbole und
für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\):
Komponenten von \(\boldsymbol F_1\): Komponenten von \(\boldsymbol F_2\): Für die gegebenen Größen:Lösung
b) Vektor von A nach P
Berechnen Sie die \((x,y,z)\)-Komponenten des Vektors \(\boldsymbol r\) von A nach P:
für die gegebenen Symbole und
für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\):
Für die gegebenen Symbole: Für die gegebenen Größen:Lösung
c) Momentvektor
Die Kräfte \(\boldsymbol F_1\) und \(\boldsymbol F_2\) erzeugen ein Moment \(\boldsymbol M_A\) um A.
Berechnen Sie die \((x,y,z)\)-Komponenten von \(\boldsymbol M_A\). Und zwar für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\). Verwenden Sie Vektorrechnung und das Kreuzprodukt.
Dies berechnen mit Komponenten liefert: Einsetzen der gegebenen Größen liefert:Lösung
d) Moment um die z-Achse
Berechnen Sie dasselbe - aber ohne Verwendung von Vektorrechnung oder Kreuzprodukt. Runden Sie auf Rundestellenwert \(0{,}01\). Zeigen Sie, dass:
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, 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
# Printing:
def printout(x, unit=1, prec="nil", sub_list=[]):
x = x.subs(sub_list)
x /= unit
try:
x = x.applyfunc(simplify)
except:
tmp = x.simplify()
if tmp != None:
x = tmp
if prec != "nil":
x = iso_round(x, prec)
pprint(x)
return
# ---
beta, a, b, c = symbols("beta, a, b, c")
F1, F2 = symbols("F1, F2", positive=True)
F1x, F1y = symbols("F_{1x}, F_{1y}", positive=True)
prec = 2
sub_list=[
(beta , 60 *deg ),
(a , S(25)/100 *m ),
(b , S(125)/1000 *m ),
(c , S(3)/10 *m ),
(F1x , 4 *N ),
(F1y , 3 *N ),
(F2 , 6 *N ),
]
pprint("\nF1 / N:")
F1 = Matrix([F1x,F1y, 0])
printout(F1, unit=N, prec="nil", sub_list=sub_list)
pprint("\nF2 / N:")
F2x = F2*cos(beta)
F2y = -F2*sin(beta)
F2 = Matrix([F2x,F2y, 0])
printout(F2, unit=N, prec=prec, sub_list=sub_list)
pprint("\nr / m:")
r = Matrix([b+c, a, 0])
printout(r, unit=m, prec=prec, sub_list=sub_list)
pprint("\nMA / Nm:")
F = F1 + F2
MA = r.cross(F)
printout(MA, unit=N*m, prec=prec, sub_list=sub_list)
pprint("\nMAz / Nm:")
MAz = - a*(F1x + F2x) + (b + c)*(F1y + F2y)
printout(MAz, unit=N*m, prec=prec, sub_list=sub_list)
[ ]: einfügen.▶ drücken.