Aufgabenstellung
Bei diesem Randwertproblem wird ein Euler-Bernoulli-Balken (Elastizitätsmodul E , Flächenmoment 2-ten Grades I ) mit einer Streckenlast q0 im Bereich A-B belastet. In A ist der Balken fest eingespannt, wobei durch eine Einbau-Ungenauigkeit der Rand im Verhältnis 1:10 geneigt ist.
In B hält das verschiebliche Lager den Balken horizontal. Gegen über Punkt A hält das Lager in B den Balken in einem vertikalen Abstand von W zur Horizontalen.
Datei:Datei:Kit5.png.png Lageplan
Gesucht ist die Biegelinie des Balkens, der durch geometrische Randbedingungen vorverformt ist.
Lösung mit Maxima
In Kit4 finden Sie die Transformation der Bewegungsgleichung des Euler-Bernoulli-Balkens in die dimensionslose Form mit der allgemeinen Lösung
w
~
i
(
ξ
)
:=
μ
⋅
ξ
4
24
+
C
i
,
3
⋅
ξ
3
6
+
C
i
,
2
⋅
ξ
2
2
+
C
i
,
1
⋅
ξ
+
C
i
,
0
{\displaystyle \displaystyle {{\tilde {w}}_{i}}\left(\xi \right):={\frac {\mu \cdot {{\xi }^{4}}}{24}}+{\frac {{{C}_{i,3}}\cdot {{\xi }^{3}}}{6}}+{\frac {{{C}_{i,2}}\cdot {{\xi }^{2}}}{2}}+{{C}_{i,1}}\cdot \xi +{{C}_{i,0}}}
für die Bereiche i=1 (A-B ) und i=2 (B-C ).
Koordinaten Dazu gehören die neuen dimensionslose Koordinaten ξ1 und ξ2 :
Die Rand- und Übergangsbedingungen lauten dann mit ℓ2 = ℓ1 /2
w
~
1
(
0
)
=
0
w
~
1
′
(
0
)
=
1
/
10
{\displaystyle {\begin{array}{ccl}{\tilde {w}}_{1}(0)&=&0\\{\tilde {w}}'_{1}(0)&=&1/10\end{array}}}
,
w
~
1
(
1
)
=
w
~
2
(
0
)
w
~
1
′
(
1
)
=
w
~
2
′
(
0
)
w
~
1
″
(
1
)
=
w
~
2
″
(
0
)
w
~
1
‴
(
1
)
=
w
~
2
‴
(
0
)
{\displaystyle {\begin{array}{ccl}{\tilde {w}}_{1}(1)&=&{\tilde {w}}_{2}(0)\\{\tilde {w}}'_{1}(1)&=&{\tilde {w}}'_{2}(0)\\{\tilde {w}}''_{1}(1)&=&{\tilde {w}}''_{2}(0)\\{\tilde {w}}'''_{1}(1)&=&{\tilde {w}}'''_{2}(0)\end{array}}}
,
w
~
1
(
1
/
2
)
=
W
w
~
2
′
(
1
/
2
)
=
0
{\displaystyle {\begin{array}{ccl}{\tilde {w}}_{1}(1/2)&=&W\\{\tilde {w}}'_{2}(1/2)&=&0\end{array}}}
tmp
Damit es einfacher wird, lassen Tilde über w weg. Die Bewegungsgleichung für beide Bereiche ist dann
w
⁗
(
ξ
)
=
μ
{\displaystyle w''''(\xi )=\mu }
und deren allgemeine Lösung
w
1
(
ξ
)
:=
μ
⋅
ξ
4
24
+
C
1
,
3
⋅
ξ
3
6
+
C
1
,
2
⋅
ξ
2
2
+
C
1
,
1
⋅
ξ
+
C
1
,
0
w
2
(
ξ
)
:=
C
2
,
3
⋅
ξ
3
6
+
C
2
,
2
⋅
ξ
2
2
+
C
2
,
1
⋅
ξ
+
C
2
,
0
(hier ist
μ
=
0
)
{\displaystyle {\begin{array}{l}\displaystyle {{w}_{1}}\left(\xi \right):={\frac {\mu \cdot {{\xi }^{4}}}{24}}+{\frac {{{C}_{1,3}}\cdot {{\xi }^{3}}}{6}}+{\frac {{{C}_{1,2}}\cdot {{\xi }^{2}}}{2}}+{{C}_{1,1}}\cdot \xi +{{C}_{1,0}}\\\displaystyle {{w}_{2}}\left(\xi \right):={\frac {{{C}_{2,3}}\cdot {{\xi }^{3}}}{6}}+{\frac {{{C}_{2,2}}\cdot {{\xi }^{2}}}{2}}+{{C}_{2,1}}\cdot \xi +{{C}_{2,0}}{\text{ (hier ist }}\mu =0)\end{array}}}
Lösen der dimensionslosen Bewegungsgleichung
Text
tmp
Einsetzten der Lösung der Bewegungsgleichungen in die Randbedingungen liefert die Gleichungen
(
C
1
,
0
=
0
C
1
,
1
=
1
10
μ
24
+
C
1
,
3
6
+
C
1
,
2
2
+
C
1
,
1
+
C
1
,
0
=
C
2
,
0
μ
6
+
C
1
,
3
2
+
C
1
,
2
+
C
1
,
1
=
C
2
,
1
μ
2
+
C
1
,
3
+
C
1
,
2
=
C
2
,
2
μ
+
C
1
,
3
=
C
2
,
3
C
2
,
3
8
+
C
2
,
2
2
+
C
2
,
1
=
0
C
2
,
3
48
+
C
2
,
2
8
+
C
2
,
1
2
+
C
2
,
0
=
1
)
{\displaystyle {\begin{pmatrix}{{C}_{1,0}}=0\\{{C}_{1,1}}=\displaystyle {\frac {1}{10}}\\\displaystyle {\frac {\mu }{24}}+{\frac {{C}_{1,3}}{6}}+{\frac {{C}_{1,2}}{2}}+{{C}_{1,1}}+{{C}_{1,0}}={{C}_{2,0}}\\\displaystyle {\frac {\mu }{6}}+{\frac {{C}_{1,3}}{2}}+{{C}_{1,2}}+{{C}_{1,1}}={{C}_{2,1}}\\\displaystyle {\frac {\mu }{2}}+{{C}_{1,3}}+{{C}_{1,2}}={{C}_{2,2}}\\\mu +{{C}_{1,3}}={{C}_{2,3}}\\\displaystyle {\frac {{C}_{2,3}}{8}}+{\frac {{C}_{2,2}}{2}}+{{C}_{2,1}}=0\\\displaystyle {\frac {{C}_{2,3}}{48}}+{\frac {{C}_{2,2}}{8}}+{\frac {{C}_{2,1}}{2}}+{{C}_{2,0}}=1\end{pmatrix}}}
mit den Unbekannten
x
_
=
(
C
1
,
0
C
1
,
1
C
1
,
2
C
1
,
3
C
2
,
0
C
2
,
1
C
2
,
2
C
2
,
3
)
{\displaystyle {\underline {x}}={\begin{pmatrix}{{C}_{1,0}}\\{{C}_{1,1}}\\{{C}_{1,2}}\\{{C}_{1,3}}\\{{C}_{2,0}}\\{{C}_{2,1}}\\{{C}_{2,2}}\\{{C}_{2,3}}\end{pmatrix}}}
.
Formulation of Boundary Conditions
Text
tmp
Das Gleichungssystem
A
_
_
⋅
x
_
=
b
_
{\displaystyle {\underline {\underline {A}}}\cdot {\underline {x}}={\underline {b}}}
hat dabei die Koeffizientenmatrix
A
_
_
=
(
1
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
1
1
2
1
6
−
1
0
0
0
0
1
1
1
2
0
−
1
0
0
0
0
1
1
0
0
−
1
0
0
0
0
1
0
0
0
−
1
0
0
0
0
0
1
1
2
1
8
0
0
0
0
1
1
2
1
8
1
48
)
{\displaystyle {\underline {\underline {A}}}={\begin{pmatrix}1&0&0&0&0&0&0&0\\0&1&0&0&0&0&0&0\\1&1&{\frac {1}{2}}&{\frac {1}{6}}&-1&0&0&0\\0&1&1&{\frac {1}{2}}&0&-1&0&0\\0&0&1&1&0&0&-1&0\\0&0&0&1&0&0&0&-1\\0&0&0&0&0&1&{\frac {1}{2}}&{\frac {1}{8}}\\0&0&0&0&1&{\frac {1}{2}}&{\frac {1}{8}}&{\frac {1}{48}}\end{pmatrix}}}
sowie die rechte Seite
b
_
=
(
0
1
10
−
μ
24
−
μ
6
−
μ
2
−
μ
0
1
)
{\displaystyle {\underline {b}}={\begin{pmatrix}0\\{\frac {1}{10}}\\-{\frac {\mu }{24}}\\-{\frac {\mu }{6}}\\-{\frac {\mu }{2}}\\-\mu \\0\\1\end{pmatrix}}}
The Equations of Motion
Text
tmp
Wir erhalten
[
C
1
,
0
=
0
,
C
1
,
1
=
1
10
,
C
1
,
2
=
72
+
5
⋅
μ
30
,
C
1
,
3
=
−
444
+
95
⋅
μ
135
,
C
2
,
0
=
2436
+
25
⋅
μ
3240
,
C
2
,
1
=
−
5
⋅
μ
−
231
270
,
C
2
,
2
=
−
24
+
μ
27
,
C
2
,
3
=
40
⋅
μ
−
444
135
]
{\displaystyle [{{C}_{1,0}}=0,{{C}_{1,1}}={\frac {1}{10}},{{C}_{1,2}}={\frac {72+5\cdot \mu }{30}},{{C}_{1,3}}=-{\frac {444+95\cdot \mu }{135}},{{C}_{2,0}}={\frac {2436+25\cdot \mu }{3240}},{{C}_{2,1}}=-{\frac {5\cdot \mu -231}{270}},{{C}_{2,2}}=-{\frac {24+\mu }{27}},{{C}_{2,3}}={\frac {40\cdot \mu -444}{135}}]}
.===Solving===
Text
tmp
w(x)
Die Lösung plotten wir, z.B. für die Streckenlast μ=100 :
Die dimensionslose Auslenkung ist jetzt etwas größer als "1" - Sie müssen sich also jetzt überlegen, ob Ihre Theorie noch gültig ist!
Post-Processing
Text
Links
Literature