Filter (Strukturanalyse)

In der Strukturdynamik wird als Filter ein Modell bezeichnet, das im Zustandsraum dargestellt wird. Diese Darstellung eines Filters als ein System ist auch in der Elektrotechnik zur Beschreibung von Schwingungssystemen üblich.

Zustandsraumdarstellung

Die Zustandsraumdarstellung besteht aus einer Steuergleichung

$\dot z(t)= A \cdot z(t)+ B \cdot u(t)$

und einer Beobachtungsgleichung

$y(t)= C \cdot z(t)+ D \cdot u(t)$

Hierbei ist z(t) der Zustandsvektor des Systems $S \left( { A, B, C, D} \right)$, u(t) entspricht dem Eingangsvektor und y(t) dem Ausgangsvektor. Die Matrizen A,B,C und D sind die Matrizen des Filters. Die Systemmatrix A beschreibt die Dynamik des Prozesses, Sie ist quadratisch; die übrigen Matrizen in der Regel nicht. B wird Eingangs- oder Steuermatrix genannt; C ist die Ausgangs- oder Beobachtungsmatrix und D ist die Durchgangsmatrix.

Formfilter (Lastfilter)

Ein Formfilter ist die Systemdarstellung einer stochastischen dynamischen Last. Dabei wird am Systemeingang gaußisches weißes Rauschen w(t) angenommen. Am Systemausgang ergibt sich der gesuchte Lastvektor f(t):

$\dot{\mathbf z}_1(t) = \mathbf A_1 \cdot \mathbf z_1(t) + \mathbf B_1 \cdot \mathbf w(t)$

$\mathbf f(t) = \mathbf C_1 \cdot \mathbf z_1(t) + \mathbf D_1 \cdot \mathbf w(t)$

Die Systemzustandsvariable $\mathbf z_1(t)$ hat keine physikalische Bedeutung. Die Systemmatrizen $\mathbf A_1,\mathbf B_1,\mathbf C_1,\mathbf D_1$ sind unbekannt und müssen bestimmt werden. Dieser Vorgang wird als Identifikation bezeichnet.

Das oben dargestellte Formfilter ist ein kontinuierliches Filter, das aus Differentialgleichungen besteht. Für die Identifikation der Systemmatrizen ist es günstiger, eine diskrete Filterdarstellung zu wählen:

$\mathbf z_{1d <k+1>}=\mathbf A_{1d}\cdot \mathbf z_{1d <k>}+\mathbf B_{1d}\cdot \mathbf w_{d <k>}$

$\mathbf f_{d <k>}=\mathbf C_{1d}\cdot \mathbf z_{1d <k>}+\mathbf D_{1d}\cdot \mathbf w_{d <k>}$

Die Matrizen des diskreten Systems sind unterschiedlich von denen des kontinuierlichen Systems. Sie können jedoch mit einer diskret-kontinuierlichen Transformation umgerechnet werden.

Die Ordnung des Systems sei mit n bezeichnet, die Anzahl der Lastkanäle (Ausgangskanäle) mit p und die Anzahl der Eingangskanäle mit q. Dann hat $\mathbf A_1$ die Größe (n x n), $\mathbf B_1$ die Größe (n x q), $\mathbf C_1$ die Größe (p x n) und $\mathbf D_1$ die Größe (p x q).

Tragwerksfilter (Strukturfilter)

Wendet man die Systemdarstellung auf die Tragwerksberechnung an, dann ist die Bewegungsgleichung eines dynamisch belasteten Tragwerks dabei die Grundlage:

$\mathbf M \cdot \ddot{\mathbf x}(t) + \mathbf D \cdot \dot{\mathbf x}(t) + \mathbf K \cdot \mathbf x(t) = \mathbf T \cdot \mathbf f(t)$

Hier ist T eine Transformationsmatrix, die die Freiheitsgrade des Tragwerks mit den Kanälen der Last f(t) verbindet. M, D und K sind die Masse-, Dämpfungs- und Steifigkeitsmatrix des Tragwerks. Die Systemzustandsvariable x(t) entspricht den Verschiebungsgrößen des Tragwerks.

Durch einfache Umformung und Hinzufügen der Identitätsgleichung

$\dot{\mathbf x} \left( t \right) = \dot{\mathbf x} \left( t \right)$

erhält man die Zustandsraumdarstellung und damit ein Filter für das mechanische System der Tragwerksreaktion:

$\begin{pmatrix} \mathbf{\dot x}(t) \\ \ddot{\mathbf x}(t) \end{pmatrix} = \begin{pmatrix} \mathbf 0 & \mathbf I \\ {-\mathbf M^{-1}\mathbf K} & {-\mathbf M^{-1}\mathbf D} \end{pmatrix} \cdot \begin{pmatrix} \mathbf x(t) \\ \dot{\mathbf x}(t) \end{pmatrix} + \begin{pmatrix} \mathbf 0 \\ {-\mathbf M^{-1}\mathbf T} \end{pmatrix} \cdot \mathbf f(t)$

Durch Einführung von Filtermatrizen $\mathbf A_2$ und $\mathbf B_2$ erhält man

$\dot{\mathbf z}_2 \left( t \right) = \mathbf A_2 \cdot \mathbf z_2 \left( t \right) + \mathbf B_2 \cdot \mathbf f\left( t \right)$

mit

$\mathbf z_2(t) = \begin{pmatrix} \mathbf x(t) \\ \dot{\mathbf x}(t) \end{pmatrix}$

Diese Filtermatrizen $\mathbf A_2$ und $\mathbf B_2$ sind durch die Tragwerksstruktur direkt gegeben und können durch ein Finite Elemente-Programm ermittelt werden. Ergänzt man eine Beobachtungsgleichung

$\mathbf z_2 \left( t \right) = \mathbf C_2 \cdot \mathbf z_2 \left( t \right) + \mathbf D_2 \cdot \mathbf f \left( t \right)$

mit $\mathbf C_2=\mathbf 1$ und $\mathbf D_2=\mathbf 0$, dann erhält man ein System $\mathbf S_2 \left( {\mathbf A_2,\mathbf B_2,\mathbf C_2,\mathbf D_2} \right)$ für das Tragwerksfilter. Dieses Filter ist ein kontinuierliches Filter, das aus Differentialgleichungen besteht.

Verbindung von Lastfilter und Strukturfilter zu einem Gesamtfilter

Die Filtergleichungen von Lastfilter $\mathbf S_1 \left( {\mathbf A_1,\mathbf B_1,\mathbf C_1,\mathbf D_1} \right)$ und Strukturfilter $\mathbf S_2 \left( {\mathbf A_2,\mathbf B_2,\mathbf C_2,\mathbf D_2} \right)$ können zu einem Gesamtfilter kombiniert werden. Bindeglied ist der Lastvektor $\mathbf f(t)$. Damit erhält man ein Filter, das sowohl Last wie Struktur beschreibt, und gaußsches weißes Rauschen $\mathbf w(t)$ als Eingangssignal hat:

$\begin{pmatrix} \dot{\mathbf x}(t) \\ \ddot{\mathbf x}(t) \\ \dot{\mathbf z}_1(t) \end{pmatrix} = \begin{pmatrix} \mathbf 0 & \mathbf I & \mathbf 0 \\ {-\mathbf M^{-1}\mathbf K} & {-\mathbf M^{-1}\mathbf D} & {\mathbf M^{-1}\mathbf T \mathbf C_1} \\ \mathbf 0 & \mathbf 0 & {\mathbf A_1} \end{pmatrix} \cdot \begin{pmatrix} \mathbf x(t) \\ \dot{\mathbf x}(t) \\ \mathbf z_1(t) \end{pmatrix} + \begin{pmatrix} \mathbf 0 \\ \mathbf 0 \\ {\mathbf B_1} \end{pmatrix} \cdot \mathbf w(t)$

$\Leftrightarrow \begin{pmatrix} \dot{\mathbf z}_2(t) \\ \dot{\mathbf z}_1(t) \end{pmatrix} = \begin{pmatrix} \mathbf A_2 & {\mathbf B_2 \mathbf C_1} \\ \mathbf 0 & \mathbf A_1 \end{pmatrix} \cdot \begin{pmatrix} \mathbf z_2(t) \\ \mathbf z_1(t) \end{pmatrix} + \begin{pmatrix} \mathbf 0 \\ \mathbf B_1 \end{pmatrix} \cdot \mathbf f(t)$

$\Leftrightarrow \dot{\mathbf z}_G(t)=\mathbf A_G\cdot \mathbf z_G(t)+\mathbf B_G\cdot \mathbf w(t)$

Dieses Filter wird bei der Kovarianzanalyse in der Strukturdynamik zur Untersuchung von Problemen mit stochastischen Lasten verwendet, um statistische Kennwerte (Varianzen und Kovarianzen) der Tragwerksantwort zu erhalten.