2.4.4.G
Gegeben
Symbole:
Größen:
Profil ist durch Querkraft \(Q_z\) (Schnittkraft in \(z\)-Richtung) beansprucht.
Berechnen Sie die Schubspannungen aufgrund Querkraft bei A und B. Gehen Sie wie folgt vor:
a) Flächenträgheitsmoment für Größen
Zeigen Sie, dass:
Berechnung wie üblich führt auf: Rundung liefert die Behauptung. Sei: \(S_{yA}\): Statisches Moment der Teilfläche rechts von A und \(S_{yB}\): Statisches Moment der Teilfläche unterhalb und rechts von B. Für die gegebenen Symbole: Berechnen Sie: Statische Momente: Berechnen Sie diese Statischen Momente in \(\mathrm{cm}^3\) (Kubikzentimeter) und gerundet auf Rundestellenwert \(0{,}1\) für die gegebenen Größen. Zeigen Sie, dass: Einsetzen liefert die Behauptung. Für diese Größen: Berechnen Sie die Schubspannungen bei A und B mit: 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 Lösung
b) Statische Momente für Symbole
Lösung
c) Statische Momente für Größen
Lösung
d) Schubspannung
SymPy
from sympy.physics.units import kg, m, s
from sympy import pi, S, pprint, symbols, cos, sin, Matrix, Eq, solve, sqrt, exp, simplify
# 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
cm = m/100
mm = m/1000
# 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
# ---
b, h = symbols("b, h")
tg, ts = symbols("t_g, t_s")
Q = symbols("Q")
sub_list=[
( b, 45 *mm ),
( h, 80 *mm ),
( tg, 8 *mm ),
( ts, 6 *mm ),
( Q, 53 *kN ),
]
pprint("\nIyy / cm⁴:")
Is = ts *h**3/12
Ig = (b-ts)*tg**3/12 + ((h - tg)/2)**2*(b-ts)*tg
Iyy = Is + 2*Ig
printout(Iyy, unit=cm**4, prec="nil", sub_list=sub_list)
printout(Iyy, unit=cm**4, prec=1, sub_list=sub_list)
A1 = tg *(b-ts)
A2 = ts * h/2
z1 = h/2 - tg/2
z2 = h/4
pprint("\n(SA, SB) / cm³:")
SA = z1 * A1
SB = z1*A1 + z2*A2
printout(SA, unit=cm**3, prec=1, sub_list=sub_list)
printout(SB, unit=cm**3, prec=1, sub_list=sub_list)
pprint("\n(τ_A, τ_B) / MPa:")
tauA = - SA * Q / Iyy / tg
tauB = - SB * Q / Iyy / ts
printout(tauA, unit=MPa, prec=1, sub_list=sub_list)
printout(tauB, unit=MPa, prec=1, sub_list=sub_list)
[ ]: einfügen.▶ drücken.