Gelöste Aufgaben/UEBA

Aufgabenstellung

Wie man die Sensoren eines Smartphones für Messungen einsetzt, zeigen wir hier.

Hier soll der Elastizitätsmodul des Stab-Materials mit Hilfe eines Biegeversuchs bestimmt werden. Der Stab der Länge ℓ hat einen quadratischen Querschnitt der Höhe h und die Masse m. Für die Messung nutzen Sie als Sensor ihr Smartphone, das Ihnen die statische Beschleunigung a_y, a_z in y- und z-Richtung angibt.

Sie haben Ihre Standard-Lösungen für den Euler-Bernoulli-Balken zu Hause vergessen und müssen sich mit einer Näherungslösung für den Balken nach dem Verfahren von Rayleigh-Ritz zufrieden geben.

Lageplan

Für den Versuch wird der Balken wird durch seine Gewichtskraft belastet. Er ist in A fest eingespannt, die Masse des Smartphones sei vernachlässigbar.

Gesucht ist die Lösung mit dem Ansatz von Rayleigh-Ritz und zwei Trial-Functions.

Gegeben: ℓ=3 m, m=28 kg, g=9.81 m/s², h=10 cm,

                 a_y = 0.9793658 m/s², a_z= 9.7609909 m/s²

Lösung mit Maxima

Smartphones haben eine Reihe von Sensoren eingebaut - und die man kann als Ingenieur prima gebrauchen ...

Ein Erklärungsvideo zu Beschleunigungs-Sensoren finden Sie hier.

tmp

Das Rayleigh-Ritz-Verfahren ist hier ein bisschen unter der experimentellen Bestimmung von Parametern des Systems "vergraben".

Der Lösungsweg ist aber der gleiche wie gewohnt.

Header

Text

1+1

tmp

Die Streckenlast auf den Balken ist

$q_0=\frac{mg}{\ell }$

und wir arbeiten mit der dimensionslosen Länge

$\xi =\frac{x}{\ell }$.

System-Parameters

Text

1+1

tmp

Formfunctions

Als Ansatzfunktion für die Näherungslösung verwenden wir

$\begin{array}{lll}w(x)=& & W\cdot \underset{=:{\varphi }_1}{\underbrace{({\xi }^2\cdot (3-2\cdot \xi ))}}\\ & +& \mathrm{\Phi }\cdot \underset{=:{\varphi }_2}{\underbrace{({\xi }^2\ell \cdot (\xi -1))}}\end{array}$

.

Koordinaten

Die Trial-Functions ϕ_i erfüllen dabei die Bedingungen

$w(\ell )=W$

und

$w^{\prime }(\ell )=\mathrm{\Phi }$

Formfunctions

Text

1+1

tmp

Für den Ritz-Ansatz brauchen wir die Terme der Potentiellen Energie des Systems

$U=\mathrm{\Pi }-A$,

wobei Π die Formänderungsenergie ist und A die Arbeitsfunktion der äußeren Last q₀.

Es ist

$\begin{array}{ll}\mathrm{\Pi }& =\frac{1}{2}\underset{\ell }{\int }EI(w^{″}{)}^{2}dx\\ A& =\underset{\ell }{\int }q(x)\cdot w(x)dx\text{ mit }q(x)=\frac{mg}{\ell }\end{array}$.

Π und A setzen wir in U ein und schreiben die skalare Gleichung in Matrizenform an. Dabei müssen wir

$\frac{d{\varphi }_i}{x}=\frac{d{\varphi }_i}{\xi }\cdot \underset{=\frac{1}{\ell }}{\underbrace{\frac{d\xi }{x}}}$

beachten. Mit den gesuchten Größen

$\underset{\_}{Q}=\left(\begin{array}{c}W\\ \mathrm{\Phi }\end{array}\right)$

erhalten wir für das Potential

$U=\frac{1}{2}EI\cdot {\underset{\_}{Q}}^T\cdot \underset{\_}{\underset{\_}{K}}\cdot \underset{\_}{Q}-{\underset{\_}{Q}}^T\cdot \underset{\_}{P}$

wobei für die Koeffizienten der Matrix K

$k_{i,j}=\frac{EI}{{\ell }^3}\cdot \underset{0}{\overset{1}{\int }}{{\varphi }_i}^{″}(\xi )\cdot {{\varphi }_j}^{″}(\xi )d\xi \text{ Achtung! Hier mit }(.{)}^{\prime }=\frac{d(.)}{d\xi }$

gilt und die rechten Seite P

$p_i=q_0\ell \cdot \underset{0}{\overset{1}{\int }}{\varphi }_i(\xi )d\xi$.

✔ Matlab© und Polynome:

Alle verwendeten Funktionen in diesem Beispiel sind Polynome. Matlab© bietet für die Manipulation von Polynomen eine Klasse von Functions an, die hier sehr hilfreich sind: https://de.mathworks.com/help/matlab/polynomials.htm

Potential Energy

Text

1+1

tmp

Die Gleichgewichtsbedingungen

$\frac{dU}{d{Q}_i}\stackrel{!}{=}0$

sind erfüllt, wenn

$\frac{EI}{{\ell }^3}\cdot \left(\begin{array}{cc}12& -6\ell \\ -6\ell & 4{\ell }^2\end{array}\right)\cdot \left(\begin{array}{c}W\\ \mathrm{\Phi }\end{array}\right)=mg\cdot \left(\begin{array}{c}\frac{1}{2}\\ \frac{\ell }{12}\end{array}\right)$.

Dieses lineare Gleichungssystem hat die Lösung

$\left(\begin{array}{c}W\\ \mathrm{\Phi }\end{array}\right)=\frac{mg{\ell }^3}{EI}\left(\begin{array}{c}\frac{1}{8}\\ \frac{1}{6\ell }\end{array}\right)$.

Zum Vergleich schauen wir jetzt bei der analytischen Lösung des Problems nach:

$\begin{array}{lll}\hat{w}& =w_{max}& =\frac{{\ell }^4}{8EI}\cdot q_0\\ \hat{\phi }& =w{\text{'}}_B& =\frac{{\ell }^3}{6EI}\cdot q_0\end{array}$.

Und sehen: unsere Lösung nach Rayleigh-Ritz stimmt am Rand exakt mit der analytischen Lösung überein, es ist

$\begin{array}{ccc}W& =& \hat{w}\\ \mathrm{\Phi }& =& \hat{\varphi }\end{array}$!

Solving

Text

1+1

tmp

Lösungsverläufe für w, φ, M, Q

Die Näherungslösung für Auslenkung w(x), Kippwinkel ϕ(x), Moment M(x) und Querkraft Q(x) können wir jetzt auftragen. Dafür verwenden wir die weiteren Referenzgrößen

$\begin{array}{ll}{M}_{ref}& =mg\ell \\ Q_{ref}& =mg\end{array}$

und tragen die Lösungen dimensionslos auf:

Für unsere Fragestellung schreiben wir die Lösung um als

$\begin{array}{ll}EI& =\frac{mg{\ell }^2}{6\mathrm{\Phi }}\\ W& =\frac{3\mathrm{\Phi }\ell }{4}\end{array}$

und bestimmen uns den Drehwinkel aus den gemessenen statischen Beschleunigungen (vgl. Bild oben). Dabei ist

$a_y=g\cdot \sin (\mathrm{\Phi })$.

🧨 Warum nehmen wir nicht die z-Komponente der Beschleunigung?:

Die az-Komponente können wir nicht gebrauchen weil ${\frac{\partial a_z}{\partial \mathrm{\Phi }}|}_{\mathrm{\Phi }=0}=0$

Da für kleine Winkel sin(Φ) ≈ Φ gilt, ist

$\mathrm{\Phi }=0.1$

Wir berechnen die gesuchte Materialkonstante mit

$I=\frac{h^4}{12}$

zu

$E=878.976.000\frac{N}{m^2}$.

Post-Processing

Text

1+1

---

Links

• ...

Literature

• ...

Lageplan