1.4.1.G
Gegeben
Symbole:
Größen:
Kreissektor mit konstanter Dicke \(t\) und konstanter Dichte \(\rho\).
Hinweis
\(\left(x,y\right)\)-Koordinaten des Ortsvektors zum Schwerpunkt S einer Scheibe mit konstanter Dicke und Dichte in Polarkoordinaten:
wobei \(A=\int {\mathsf d}A\) der Flächeninhalt der Scheibe ist.
a) Flächeninhalt und Masse
Berechnen Sie die Fläche \(A\) der Scheibe und ihre Masse \(m\).
Für die gegebenen Symbole: Berechnen Sie:
Für die gegebenen Größen: Berechnen Sie \(A\) in \(\mathrm{m}^2\) (Quadratmeter) und \(m\) in \(\mathrm{kg}\) (Kilogramm). Runden Sie auf Rundestellenwert \(0{,}01\). Zeigen Sie, dass:
Für die gegebenen Symbole: Für die gegebenen Größen:Lösung
b) Schwerpunktposition
Berechnen Sie symbolisch die Koordinaten \(\left(x_m, y_m \right)\) des Schwerpunkts S der Scheibe.
Für die gegebenen Symbole: Zeigen Sie, dass:
Für die gegebenen Größen: Runden Sie auf Rundestellenwert \(0{,}01\), und zeigen Sie, dass:
Für die gegebenen Symbole: Für die gegebenen Größen und gerundet auf Rundestellenwert \(0{,}01\):Lösung
c) Zusatzmasse
Die Masse der Scheibe ist:
Es wird eine Punktmasse der Masse \(M=2m\) hinzugefügt.
Für die gegebenen Symbole: Berechnen Sie die Koordinaten \(\left(x_M, y_M\right)\) des Schwerpunkts der Punktmasse, sodass der Schwerpunkt des Verbunds der beiden Massen im Koordinatenursprung O liegt:
Für die gegebenen Größen: Runden Sie auf Rundestellenwert \(0{,}01\). Zeigen Sie, dass:
Lösung
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 *
from sympy import *
# Units:
(mm, cm) = ( m/1000, m/100 )
kN = 10**3*newton
Pa = newton/m**2
MPa = 10**6*Pa
GPa = 10**9*Pa
deg = pi/180
half = S(1)/2
# Rounding:
def iso_round(x, pv):
try:
x = float(x)
tmp = round(x, pv)
except TypeError:
from functools import partial
func = partial(round, ndigits = pv)
tmp = x.applyfunc(func)
return tmp
def printout(x, unit=1, prec="nil", sub_list=[]):
x = x.subs(sub_list)
x /= unit
try:
x = x.applyfunc(simplify)
except:
x = x.simplify()
if prec != "undefined":
x = iso_round(x, prec)
pprint(x)
return
# ---
var("R, t, rho")
sub_list=[
(R, 1 *m),
(t, 1 *m/10),
(rho, 1000 *kg/m**3),
]
prec = 2
A = pi * R*R / 6
mass = rho * A * t
pprint("\nA:")
pprint(A)
pprint("\nm:")
pprint(mass)
pprint("\nA / m²:")
printout(A, unit=m*m, prec=prec, sub_list=sub_list)
pprint("\nm / kg:")
printout(mass, unit=kg, prec=prec, sub_list=sub_list)
pprint("\nxm / m:")
xm = R/pi *sqrt(3)
printout(xm, unit=m, prec=prec, sub_list=sub_list)
pprint("\nym / m:")
ym = R/pi
printout(ym, unit=m, prec=prec, sub_list=sub_list)
pprint("\nxM / m:")
xM = -xm/2
printout(xM, unit=m, prec=prec, sub_list=sub_list)
pprint("\nyM / m:")
yM = -ym/2
printout(yM, unit=m, prec=prec, sub_list=sub_list)
[ ]: einfügen.▶ drücken.