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.
Die Zustandsraumdarstellung besteht aus einer Steuergleichung
![{\displaystyle {\dot {z}}(t)=A\cdot z(t)+B\cdot u(t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/662562f476ba6706e13808b379f30da78c501eb9)
und einer Beobachtungsgleichung
![{\displaystyle y(t)=C\cdot z(t)+D\cdot u(t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f283691f15b7e0cf30ffb85e4916efbe03332d05)
Hierbei ist
der Zustandsvektor des Systems
,
entspricht dem Eingangsvektor und
dem Ausgangsvektor. Die Matrizen
und
sind die Matrizen des Filters. Die Systemmatrix
beschreibt die Dynamik des Prozesses, sie ist quadratisch; die übrigen Matrizen in der Regel nicht.
wird Eingangs- oder Steuermatrix genannt;
ist die Ausgangs- oder Beobachtungsmatrix und
ist die Durchgangsmatrix.
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):
![{\displaystyle {\dot {\mathbf {z} }}_{1}(t)=\mathbf {A} _{1}\cdot \mathbf {z} _{1}(t)+\mathbf {B} _{1}\cdot \mathbf {w} (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56fa61282e4ff3fed19cc85f4a5312ddeb34cbf0)
![{\displaystyle \mathbf {f} (t)=\mathbf {C} _{1}\cdot \mathbf {z} _{1}(t)+\mathbf {D} _{1}\cdot \mathbf {w} (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93658d10097949984ea6a7a33511016557e5bb26)
Die Systemzustandsvariable
hat keine physikalische Bedeutung. Die Systemmatrizen
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:
![{\displaystyle \mathbf {z} _{1d<k+1>}=\mathbf {A} _{1d}\cdot \mathbf {z} _{1d<k>}+\mathbf {B} _{1d}\cdot \mathbf {w} _{d<k>}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7db01943ff1e920935c49a604cb2743a8c03a2f)
![{\displaystyle \mathbf {f} _{d<k>}=\mathbf {C} _{1d}\cdot \mathbf {z} _{1d<k>}+\mathbf {D} _{1d}\cdot \mathbf {w} _{d<k>}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/917d78477d7aa17bba805306763461f5a83c98db)
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
die Größe (n x n),
die Größe (n x q),
die Größe (p x n) und
die Größe (p x q).
Wendet man die Systemdarstellung auf die Tragwerksberechnung an, dann ist die Bewegungsgleichung eines dynamisch belasteten Tragwerks dabei die Grundlage:
![{\displaystyle \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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a260722434d0fa2a94f77238d4064407cc1316bc)
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
![{\displaystyle {\dot {\mathbf {x} }}\left(t\right)={\dot {\mathbf {x} }}\left(t\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/61c51830572fe7ffd43bb7af9039c39a9e8fc98b)
erhält man die Zustandsraumdarstellung und damit ein Filter für das mechanische System der Tragwerksreaktion:
![{\displaystyle {\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e8ad679fce91518fdd224e643cc0cdfd394bd93b)
Durch Einführung von Filtermatrizen
und
erhält man
![{\displaystyle {\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c4694f4a7b90ccb8dcb6c28215100c01636f7235)
mit
![{\displaystyle \mathbf {z} _{2}(t)={\begin{pmatrix}\mathbf {x} (t)\\{\dot {\mathbf {x} }}(t)\end{pmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5f24fbe963c056dfd03c3098b23145ab7d253d0)
Diese Filtermatrizen
und
sind durch die Tragwerksstruktur direkt gegeben und können durch ein Finite-Elemente-Programm ermittelt werden. Ergänzt man eine Beobachtungsgleichung
![{\displaystyle \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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a44c471a516a34bc6d3318681454a39d1abe5ece)
mit
und
, dann erhält man ein System
für das Tragwerksfilter. Dieses Filter ist ein kontinuierliches Filter, das aus Differentialgleichungen besteht.
Die Filtergleichungen von Lastfilter
und Strukturfilter
können zu einem Gesamtfilter kombiniert werden. Bindeglied ist der Lastvektor
. Damit erhält man ein Filter, das sowohl Last wie Struktur beschreibt, und gaußsches weißes Rauschen
als Eingangssignal hat:
![{\displaystyle {\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b40e8855076de64976a767f151e0b8cf56b0d273)
![{\displaystyle \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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5bd7f356c123908c77c699f5d4368350bf93f3ba)
![{\displaystyle \Leftrightarrow {\dot {\mathbf {z} }}_{G}(t)=\mathbf {A} _{G}\cdot \mathbf {z} _{G}(t)+\mathbf {B} _{G}\cdot \mathbf {w} (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c33add360ec99dfb0e6a56c30866680649bfc0a2)
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.