2.1.C
Gegeben
Symbole:
Drei Stäbe mit Querschnittsfläche \(A\). Stab 1 und 2 haben die Länge \(l\).
Stab 3 hat die Länge \(\tfrac l \beta\).
Zeigen Sie, dass für die 3 Stabkräfte gilt:
Berechnen Sie dazu die unvollständige Lösung mit 2·k + 2·s - r = 8 Gleichungen und Unbekannten. Gehen Sie wie folgt vor:
a) Gleichgewicht
Machen Sie einen Knotenschnitt am Knoten 1. Geben Sie zwei Gleichgewichtsbedingungen an.
b) Kinematik
Einheitsvektoren und Stabverlängerungen: Geben Sie die (genäherten) Stabverlängerungen \(\Delta l_1, \Delta l_2, \Delta l_3\) an abhängig von den Verschiebungsvektor-Komponenten \((u_{1x}, u_{1y})\). Verwenden Sie die Einheitsvektoren in Stabrichtung \(\boldsymbol e_1, \boldsymbol e_2, \boldsymbol e_3\).
Lösung
c) Elastizität
Geben Sie die Verlängerungen der Stäbe an abhängig von den Stab-Zugkräften und gegebenen Symbolen.
Lösung
d) Gleichungssystem und Lösung
Lösen Sie dazu das Gleichungssystem, und berechnen Sie die Stabkräfte.
Aus den acht Gleichungen lassen sich die acht Unbekannten \(S_1, S_2, S_3, u_{1x}, u_{1y}, \Delta l_1, \Delta l_2, \Delta l_3\) berechnen. Lösung für Stabkräfte:Lösung
e) Stabkräfte
Berechnen Sie die Stabkräfte für \(\beta=1\).
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 kg, m, s
from sympy import pi, S, pprint, symbols, cos, sin, Matrix, Eq, solve, sqrt, exp
# 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
S1, S2, S3 = symbols('S1 S2 S3')
u1x, u1y = symbols('u1x u1y')
dl1, dl2, dl3 = symbols('dl1 dl2 dl3')
# Given:
beta = symbols("beta")
# beta = 1
F, l, EA = symbols('F, l, EA')
if beta == 1:
pprint("\nSolution with β = 1:")
else:
pprint("\nSolution with β given as given symbol:")
# Solve for 8 unknowns:
unknowns = [S1, S2, S3, u1x, u1y, dl1, dl2, dl3]
s = sqrt(2)/2
eqns = [
Eq( S1 + S2*s , F ),
Eq( S3 + S2*s , 0 ),
Eq( dl1 , u1x ),
Eq( dl2 , s*(u1x + u1y) ),
Eq( dl3 , u1y ),
Eq( EA*dl1/l , S1 ),
Eq( EA*dl2/l , S2 ),
Eq( EA*dl3/l*beta , S3 ),
]
sol = solve(eqns, unknowns)
S1, S2, S3, u1x, u1y, dl1, dl2, dl3 = \
sol[S1],sol[S2],sol[S3],sol[u1x],sol[u1y],sol[dl1],sol[dl2],sol[dl3]
pprint("\nu1x:")
tmp = u1x
tmp = tmp.simplify()
pprint(tmp)
pprint("\nu1y:")
tmp = u1y
tmp = tmp.simplify()
pprint(tmp)
pprint("\ndl1, dl2, dl3:")
for dl in [dl1, dl2, dl3]:
tmp = dl
tmp = tmp.simplify()
pprint(tmp)
pprint("\nS1, S2, S3:")
for tmp in [S1, S2, S3]:
tmp = tmp.simplify()
pprint(tmp)
[ ]: einfügen.▶ drücken.