ZfP in Anwendung, Entwicklung und Forschung
Print | | | |
| DGZfP-JAHRESTAGUNG 2003 ZfP in Anwendung, Entwicklung und Forschung | |||
| Start > Beiträge > Vorträge > Grundlagen der Ultraschallprüfung: | Print |
| Experimentelle Verifizierung von Modell-vorstellungen zur Schallausbreitung in anisotropen MaterialienH. Ernst, BAMV. Munikoti, BAM H. Wüstenberg, BAM Kontakt: Ernst Hardy Dipl.-Ing. 1. EinführungDie Notwendigkeit ein konkretes Verständnis über die teilweise sehr komplexen Wellenausbreitungsvorgänge in akustisch anisotropen Werkstoffen zu entwickeln, führte zu einer Reihe von Modellansätzen die sehr genaue Aussagen zu diesen Schallausbreitungsphänomenen machen. Die Berechnung der Schallausbreitung erfolgt dabei nach analytischen (z.B. [1], [2], [3] und [4]), approximativen (z.B. [5], [6], [7] und [8]) oder numerischen (z.B. [9] und [10]) Ansätzen. Besonders aufschlussreiche Erkenntnisse können aus den Modellierungen gezogen werden, wenn Wellenfelder in der Einschallebene der Prüfköpfe zu festen Zeitpunkten dargestellt werden. Anhand dieser anschaulichen Wellenfeldbilder kann das Ausbreitungsverhalten elastischer Wellen in anisotropen Werkstoffen gut studiert werden. Die messtechnische Verifizierung der Aussagen erfolgte bis jetzt nur durch die Messung einzelner Teilaspekte der Wellenausbreitung wie zum Beispiel der Richtungsabhängigkeit der Schallgeschwindigkeit und der radialen Schalldruckverteilung der quasi-Longitudinalwelle (qL). Mit den Ergebnissen dieser konventionellen Messmethoden war es bisher jedoch nicht möglich, die in zahlreichen Modellansätzen berechnete asymmetrische Form der Wellenfelder sowie das ungewöhnliche Fokussierungs- bzw. Divergenzverhalten speziell der quasi-Transversalwelle (qT1) in einer anschaulichen Form darzustellen und somit zu verifizieren. Diese Arbeit hat sich daher zum Ziel gesetzt, die Schallausbreitung in akustisch anisotropen Medien messtechnisch möglichst vollständig nach Ort- und Zeit zu erfassen, so dass eine anschauliche Beschreibung der Schallausbreitungsvorgänge durch Momentaufnahmen der Wellenfelder möglich ist. Dazu wird eine Visualisierungstechnik [11], [12] eingesetzt, mit der es möglich ist Momentaufnahmen der Wellenfeldausbreitung auf Festkörperoberflächen zu erstellen. In dieser Arbeit werden einige der wichtigsten modellierten Schallausbreitungsphänomene in akustisch anisotropen Medien mit gemessenen und anschaulich visualisierten Darstellungen verglichen und bestätigt. 2. Probenmaterial2.1. Si-Einkristall (kubisch)
Die Untersuchungen an dem einkristallinen Silizium (p-Typ, Bor dotiert) wurden an der in Abbildung 1 dargestellten Probengeometrie durchgeführt. Eine der drei <100>-Symmetrie-achsen des kubischen Kristallsystems ist entlang der Zylinderachse orientiert.
2.2. Austenitischer Stahlguss X 5 CrNi 18 10 (transversal isotrop)
Die elastischen Konstanten des austenitischen Stahlgusses wurden durch Schallgeschwindigkeitsmessungen bestimmt.
3. Modellierung der SchallausbreitungZiel war es das komplexe Ausbreitungsverhalten der von den Prüfkopfaperturen emittierten Wellenfelder in der Einschallebene der Prüfköpfe in Form von Momentaufnahmen anschaulich darzustellen. Die Berechnungen sollen häufig auftretende Anwendungsfälle bei der Werkstoffprüfung wie gekrümmte Grenzflächen und Schrägeinschallung einschließen. Die Berechnung der prüfkopfgenerierten transienten Felder in anisotropen Materialien erfolgt nach dem Huygenschen Prinzip durch eine linienförmige Aneinanderreihung von Punktquellen und Superposition der einzelnen Punktquellenfelder. Abbildung 3 zeigt die geometrische Konfiguration von zehn Punktquellen Q, welche die Schwingeroberfläche ersetzten.
Zur Berechnung der komplexen Teilchenverschiebung u an einem beliebigen Aufpunkt P(x,z) innerhalb des anisotropen Materials werden die von den Quellpunkten Q emittierten Wellen entsprechend ihrer Laufzeit und Richtwirkung zum Punkt P(x,z) aufsummiert. Dabei muss beachtet werden, dass der Energiefluss jeder durch die Richtung des Teilchenverschiebungsvektors v (Polarisation) gekennzeichneten Wellenart in einem anisotropen Material mit Gruppengeschwindigkeit erfolgt. Deswegen muss bei den Berechnungen anstelle der Phasengeschwindigkeit die Gruppengeschwindigkeit cg verwendet werden [13]. Desweiteren muss beachtet werden, dass die Punktrichtwirkungen x in anisotropen Materialien eine Funktion der jeweiligen kristallographischen Orientierung sind und im allgemeinen eine asymmetrische Form haben. Da die Berücksichtigung dieser richtungsabhängigen Punktrichtwirkungen eine genauere Berechnung der Schallausbreitung ermöglicht, wurden sie mit Hilfe des Reziprozitätstheorems [14], [15] berechnet. Um die transienten Felder möglichst genau zu berechnen, wurde die Zeitfunktion der in den Experimenten verwendeten Prüfköpfe sowie die Impulsantwort der abtastenden elektrodynamischen Sonde in den Simulationsrechnungen berücksichtigt. Zur Berechnung der transienten Wellenfelder wurden die zeitharmonischen Teilchenverschiebungsfelder für diskrete Komponenten fi des Frequenzspektrums des verwendeten Prüfkopfes superpositioniert. Durch die a priori eingeführte Verzögerungszeit t0 wurde eine Laufzeitverschiebung erzeugt, durch die der Schwerpunkt des modellierten Impulses im Zeitbereich verschoben werden konnte. Die komplexe Teilchenverschiebung u im Aufpunkt P(x,z) des anisotropen Materials zum Zeitpunkt t0 bei Anregung durch einen Stoßwellenprüfkopf auf einer freien Oberfläche ergibt sich durch Superposition von k Punktquellenfeldern,
wobei rm der Abstand zwischen der Punktquelle Qm und dem Aufpunkt P(x,z) innerhalb des anisotropen Materials (siehe Abb.: 3), fi die i-te diskrete Frequenzkomponente des Spektrums des verwendeten Prüfkopfes und Ai die zugehörige Amplitude, cg die Gruppengeschwindigkeit der jeweiligen Wellenart I (Longitudinal oder Transversal) und xd,I die Punktrichtwirkung mit d = z für normale und d = x für tangentiale Anregung ist. Die Anzahl der Punktquellen k ergibt sich aus dem Durchmesser der Prüfkopfapertur, da der Abstand der einzelnen Punktquellen kleiner sein muss als die Hälfte der kleinsten Wellenlänge der erzeugten Wellenart um das Sampling-Theorem einzuhalten. Abbildung 4 (links) zeigt den Realteil des komplexen Teilchenverschiebungsfeldes einer qL-Welle zum Zeitpunkt t0=4ms bei impulsartiger Anregung durch einen Normalprüfkopf mit einem Durchmesser von 12mm in einem -22,5° um die Y-Achse rotierten kubischem Silizium-Einkristall. Einschallebene ist die (010)-Ebene des Silizium Einkristalls. Im rechten Teil von Abbildung 4 ist das asymmetrische Schallfeld des Normalprüfkopfes (Betrag der komplexen Teilchenverschiebung u) dargestellt.
Die Implementierung des Einflusses der zylindrischen Bauteilkrümmung und des Einschallwinkels ßPl in der Plexiglasvorlaufstrecke erfolgte auf Basis einer Näherung, welche davon ausgeht, dass sich in den kurzen Plexiglasvorlaufstrecken ebene Wellen ausbreiten. Aus der Geometrie der Vorlaufstrecke ergeben sich die neuen Koordinaten der einzelnen Punktquellen auf der Grenzfläche. Die Berücksichtigung der unterschiedlichen Schallaufwege zwischen den Punktquellen der Prüfkopfapertur Qm und den neuen Punktquellen auf der Grenzfläche erfolgt durch eine zeitverzögerte Anregung. Desweiteren wurde ein einfaches Strahlenmodell entwickelt, mit dem die Wellenfrontausbreitung durch eine gekrümmte Grenzfläche zwischen isotropem Vorlaufmaterial und anisotropem Material berechnet werden kann. Voraussetzung für diesen Ansatz ist die Modellierung der Wellenfront im Vorlaufmaterial (Plexiglas) mittels Punktquellensynthese. Das Modell postuliert, dass jeder Punkt der Wellenfront Ursprung einer senkrecht zu ihr orientierten ebenen Welle ist. Die Richtung der Strahlen repräsentieren die Richtung der Wellenausbreitung. Es muss darauf hingewiesen werden, dass die beiden Modellansätze nur für die Berechnung der Wellenausbreitung in den {100}-Symmetrieebenen des kubischen Kristallsystems (Silizium) beziehungsweise für die Meridianebene des transversal isotropen Kristallsystems (austenitischer Stahlguss) gelten. Für generelle Ausbreitungsprobleme in welchen die Einschallebene nicht in einer Symmetrieebene des Kristallsystems liegt muss auf die in der Einführung erwähnten numerischen und approximativen Verfahren zurückgegriffen werden. 4. MessanordnungEine schematische Darstellung der Messanordnung zur Visualisierung der Wellenausbreitung in opaken Festkörpern [11], [12] zeigt Abbildung 5. Die Messanordnung besteht aus einem konventionellen Ultraschallmessgerät, einem Sendeprüfkopf, einer elektrodynamischen Messsonde, einem Manipulatorsystem und einem Messrechner mit Transientenrecorder zur A/D Wandlung der HF-Signale.
Um die räumliche Ausbreitung der Wellenfelder zu visualisieren wurde zunächst der relevante Bereich auf der Probeoberfläche mäanderförmig abgerastert. Das von der Sonde empfangene Zeitsignal (A-Bild) wurde digitalisiert und anschließend im Rechner abgespeichert. Im linken Teil von Abbildung 6 ist der Aufbau einer auf diese Weise erstellten Messdatenmatrix schematisch dargestellt. Die vertikalen Schnittebilder durch die Messdatenmatrix ergeben B-Bilder in denen der zeitliche Signalverlauf über der Abtastlinie dargestellt ist.
Ein horizontales Schnittbild durch die Messdatenmatrix zu einem definierten Zeitpunkt t ergibt ein C-Bild, in welchem die Amplituden über dem Ort abgebildet sind. In diesen C-Bildern ist die räumliche Ausbreitung der Ultraschallwellenfronten des vom Sendeprüfkopf emittierten Impulses zum Zeitpunkt t abgebildet. Durch die definierte Änderung des Zeitpunktes kann somit die zeitliche und räumliche Ausbreitung des Ultraschallimpulses visualisiert werden. Dargestellt in Abbildung 6 rechts ist das C-Bild mit dem Wellenfeld eines senkrecht in feinkörnigem isotropen ferritischem Stahl einschallenden Transversalwellenprüfkopfes (V153/34504/30/1MHz/0,5") nach einer Schalllaufzeit von 5 ms. 5. Vergleich von Messung und ModellierungAbbildung 7 zeigt für drei Zeitpunkte (1ms, 4ms und 7ms) gemessene und berechnete Schnappschüsse der Wellenausbreitung der qT1-Welle für Einschallung in die kristallographische [100]-, [210]- und [110]-Richtung des Silizium-Einkristalls. In den Prinzipskizzen über den Spalten mit den zugehörigen Wellenfeldern ist die kristallographische Orientierung des Si-Einkristalls relativ zur Einschallrichtung schematisch dargestellt. Deutlich zu erkennen sind die unterschiedlichen Formen und Geschwindigkeiten der Wellenfelder für die verschiedenen Einschallrichtungen in Silizium-Einkristall.
In Abbildung 8 sind die visualisierten Wellenfelder der qL-Welle für Einschallwinkel relativ zur Stengelkristallitorientierung von 0°, 22,5° und 45° den Berechnungen basierend auf dem Strahlenmodell und der Punktquellensynthese gegenübergestellt. Die Schalllaufzeit war für alle Momentaufnahmen der Wellenausbreitung konstant bei 7,6 ms. Abbildung 8 verdeutlicht in anschaulicher Weise, wie drastisch die Ausbreitung der Wellenfelder in dem transversal isotropem Stahlguss von der Orientierung der Stengelkristallite bestimmt wird.
In Abbildung 9 sind die gemessenen und berechneten Wellenfelder der qT1-Welle für Einschallwinkel von 0° und -10° relativ zur Stengelkristallitorientierung dargestellt. Die Schnappschüsse der Wellenausbreitung entsprachen einer Laufzeit von jeweils 8 ms im Stahlguss. In Abbildung 9 a ist für eine Einschallung in Richtung der Stengelkristallite die ungewöhnliche sichelförmige Form des Wellenfeldes zu erkennen, die schon von der Einschallung in die [100]-Richtung des Si-Einkristalls (siehe Abbildung 7) bekannt ist und auch hier auf die konkave Krümmung der Slownesskurve in dieser Richtung zurückzuführen ist. Wie empfindlich die qT1-Welle bereits auf geringe Änderungen der Stengelkristallitorientierung reagiert zeigt Abbildung 9 b. Ein Rotation der Stengelkristallitorientierung um lediglich 10° entgegen dem Uhrzeigersinn resultiert in eine deutliche Ablenkung des Wellenfeldes der qT1-Welle.
In Abbildung 10 sind gemessene und berechnete Momentaufnahmen der Ausbreitung der qT1-Wellenfelder in austenitischem Stahlguss gegenübergestellt, wobei die qT1-Wellen durch schrägen Einfall von Longitudinalwellen durch Brechung erzeugt wurden. Eingeschallt wurde mit einem konventionellem Winkelprüfkopf mit einem Keilwinkel von ßPl=50° welcher in isotropem ferritischem Stahl eine Transversalwelle unter 70° erzeugen würde. Die Einschallebene war die Meridianebene, wobei die Stengelkristallite senkrecht zur Koppelfläche orientiert waren. Über den Momentaufnahmen der Wellenausbreitung in Abbildung 10 ist eine Schema der Messanordnung mit der Orientierung der Stengelkristallite relativ zur Einschallrichtung sowie dem Messbereich auf der Probenoberfläche abgebildet. Abbildung 10 zeigt, wie ungewöhnlich die Wellenfeldausbreitung in dem transversal isotropem Stahlguss sein kann. Bei Einschallung entsprechend der dargestellten Messanordnung werden im austenitischen Stahlguss zwei Transversalwellen erzeugt, welche sich mit unterschiedlicher Geschwindigkeit in verschiedenen Richtungen ausbreiten. Die Ursache dieses Phänomens lässt sich anhand des Konzeptes der Slownesskurve der ebenen Welle ableiten. Infolge der konkaven Krümmung werden zwei physikalisch ausbreitungsfähige qT1-Wellen erzeugt. Besonders ungewöhnlich ist das Ausbreitungsverhalten der zweiten qT1-Welle. Das Wellenfeld der zweiten qT1-Welle breitet sich mit größerem Einschallwinkel und deutlich höherer Geschwindigkeit als die erste qT1-Welle aus, wobei die Wellenfronten beziehungsweise der Wellenvektor in Richtung der Grenzfläche zum Plexiglas zeigt. Die starke Zerklüftung der Wellenfronten ist eine Folge der Inhomogenität des Untersuchungsmaterials und der dadurch hervorgerufenen Streuung an den Korngrenzen.
6. Zusammenfassung
1 Literaturverzeichnis
| |
|
| |
