WO2002008793A1

WO2002008793A1 - Verfahren zur bestimmung der position eines sensorelementes

Info

Publication number: WO2002008793A1
Application number: PCT/CH2001/000431
Authority: WO
Inventors: Paul G. Seiler; Ralph K. Muench; Stefan Kirsch
Original assignee: Northern Digital Inc.
Priority date: 2000-07-26
Filing date: 2001-07-10
Publication date: 2002-01-31
Also published as: CA2441226A1; US20030200052A1; ATE312364T1; EP1303771B1; EP1303771A1; AU2001267254A1; CA2441226C; DE50108329D1; CN1330978C; US6836745B2; CN1459030A; JP2004505253A

Abstract

Die Erfindung betrifft zunächst ein Verfahren zur Bestimmung der Position eines Sensorelementes (300), mit Hilfe dessen ein von mindestens einer Feldgeneratoreinheit (200) ausgesendetes magnetisches Wechselfeld (210) gemessen wird, wobei aufgrund eines im Sensorelement (300) empfangenen Signals die Position des Sensorelementes (300) bestimmt wird. Erfindungsgemäss ist vorgesehen, dass, vorzugsweise in erster Näherung, Störfelder (410) berechnet werden, welche aufgrund von in elektrisch leitenden Objekten (400) erzeugten Wirbelströmen (420) entstehen. Des weiteren ist vorgesehen, dass die Position, welche ausgehend vom im Sensorelement (300) empfangenen Signal bestimmbar ist, aufgrund der berechneten Störfelder (410) korrigiert wird.

Description

Verfahren zur Bestimmung der Position eines Sensorelementes

Die vorliegende Erfindung betrifft ein Verfahren nach dem Oberbegriff des Patentanspruchs 1, eine Anwendung des Verfahrens, eine Vorrichtung zur Durchführung des Verfahrens sowie ein Computerprogrammprodukt .

Bei zahlreichen technischen und medizinischen Verfahren ist die genaue Kenntnis der Position eines bestimmten

Gegenstandes von entscheidender Bedeutung. Währenddem in der Medizin die Position von einzelnen Gewebeteilen - beispielsweise eines Tumors, der zur Zerstörung oder zur Wachstumsbegrenzung bestrahlt werden soll - bestimmt werden muss, ist die Positionserfassung zur Eingabe in ein Computersystem, beispielsweise für „Cyber Space™- Anwendungen, von allgemeiner Bedeutung. Eine solche Positionserfassungs- bzw. Positionseingabeeinheit wird in diesen Anwendungen auch etwa als dreidimensionale Maus bezeichnet.

Eine bekannte Vorrichtung bzw. ein bekanntes Verfahren zur Bestimmung der Position ist in der internationalen Patentanmeldung mit der Verδffentlichungsnummer WO 97/36192 der gleichen Anmelderin beschrieben. Gemäss der bekannten Lehre ist es vorgesehen, ein Wechselfeld mit Hilfe einer Feldgeneratoreinheit aufzubauen, wobei je nach der Anzahl Freiheitsgrade eines Sensorelementes, dessen Position bestimmt werden soll, mehrere Wechselfelder einander überlagert werden. Mit Hilfe einer Verarbeitungs- und

Steuereinheit, welche die Feldgeneratoreinheit einerseits steuert, anderseits die vom Sensorelement empfangenen Signale verarbeitet, wird die Position und gegebenenfalls die Lage der Sensoreinheit bestimmt. Der Inhalt der vorstehend genannten Veröffentlichung bildet diesbezüglich integrierender Bestandteil dieser Beschreibung.

Es hat sich gezeigt, dass bei einer magnetfeldbasierten Ortung, wie sie beispielsweise bei der bekannten Lehre gemäss WO 97/36192 angewandt wird, in benachbarten, elektrisch leitenden Objekten Wirbelstrδme erzeugt werden. Diese führen zu Verzerrungen des ursprünglichen magnetischen Wechselfeldes und damit zu systematischen

Fehlern. Das bedeutet, dass, wenn Position und Ausrichtung von Sensorelementen im verzerrten Wechselfeld so bestimmt werden, als wäre kein elektrisch leitendes Objekt zugegen, die gewonnenen Werte systematisch verfälscht sind.

Ein Verfahren zur Kompensation von Stδreffekten, welche durch leitende Objekte hervorgerufen werden, ist unter dem Namen „Distortion Mapping" bekannt. Dieses Verfahren wird beispielsweise im Aufsatz mit dem Titel „Calibration of Tracking Systems in a Surgical Environment" (Birkfellner et. al., IEEE Trans Med Imaging, Vol. 17(5), Seiten 737 bis 742, 1998) beschrieben. Beim bekannten Verfahren wird die Position und die Orientierung eines Sensorelementes ebenfalls mit Hilfe eines Positionsmesssystems, das auf einer magnetfeldbasierten Ortung beruht, vorgenommen, wobei zur Kompensation von Stδreffekten ein zweites Positionsmesssystems vorgesehen ist, welches durch elektrisch leitende Objekte nicht beeinflussbar ist. Die Differenz zwischen den mit den beiden Positionsmesssystemen ermittelten Positionen und Orientierungen werden in der

Folge zur Korrektur der mit Hilfe des magnetfeldbasierten Positionsmesssystems bestimmten Position und Orientierung verwendet. Das bekannte Verfahren weist jedoch den Nachteil auf, dass, zur Erreichung einer hohen Genauigkeit, die Positions- und Orientierungsdifferenz an möglichst vielen Punkten gemessen werden muss. Um weitere Punkte zu erhalten, muss zusätzlich ein aufwendiges Interpolationsverfahren angewendet werden. Der sehr hohe Aufwand wird insbesondere durch das folgende Beispiel deutlich: Soll ein Volumen von 1 m³ vermessen werden, wobei dies in den drei Achsen alle 10cm und in zehn verschiedenen Orientierungswinkeln geschehen soll, so erhält man 10' 000 Punkte. Darüber hinaus ist das erwähnte zweite Positionsmesssystem erforderlich.

Des weiteren ist ein Verfahren zur Kompensation von Stδreffekten bekannt, bei dem Magnetfelder durch gepulste Gleichstromfelder erzeugt werden, wobei die Kompensation von Wirbelstromeffekten dadurch vorgenommen werden, indem Magnetfeldmessungen erst nach dem Abklingen der im Messsignal enthaltenen WirbelStromanteile durchgeführt werden. Weiterführende Erläuterungen zum bekannten Verfahren können den Druckschriften US-5 453 686 und US-5 767 669 entnommen werden. Es hat sich gezeigt, dass die Genauigkeit der ermittelten Resultate ungenügend ist. Insbesondere ist die Kompensation unvollständig, wenn die Zerfallszeiten der WirbelStromanteile die Pulszeit zwischen zwei Gleichstrompulsen übersteigt. Zwar kann dem durch

Verlängerung der Pulszeit begegnet werden, doch führt dies zu einer unerwünschten geringeren Messrate. Des weiteren lässt sich die bekannte Kompensationsmethode nicht bei auf magnetischer Ortung basierenden Positionsmesssystemen anwenden, die magnetische Wechselfelder erzeugen.

Der vorliegenden Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren anzugeben, das eine verbesserte Bestimmung der Position und/oder der Lage eines Sensorelementes ermöglicht.

Diese Aufgabe wird durch die im kennzeichnenden Teil des Patentanspruchs 1 angegebenen Massnahmen gelöst. Vorteilhafte Ausgestaltungen der Erfindung, eine Anwendung des Verfahrens, eine Vorrichtung zur Durchführung des Verfahrens sowie ein Computerprogrammprodukt sind in weiteren Ansprüchen angegeben.

Mit dem erfindungsgemässen Verfahren wird es ermöglicht, den Einfluss von leitenden Objekten zu eliminieren, zumindest aber erheblich zu reduzieren. Des weiteren ist dieses Verfahren allgemeiner und genauer als die bekannten Verfahren. Schliesslich kann der geometrieabhängige Teil der Berechnungen im Sinne einer Systemkalibrierung vor dem eigentliche Einsatz des Positionsmesssystems vorgenommen werden.

Die Erfindung wird nachfolgend anhand von Zeichnungen beispielsweise näher erläutert. Dabei zeigen

Fig. 1 eine bekannte Anordnung, bestehend aus Feldgeneratoreinheit, Sensorelement und Verarbeitungs- und Steuereinheit, in schematischer Darstellung, mit einem elektrisch leitenden Objekt,

Fig. 2 ein elektrisch leitendes Objekt und

Fig. 3 ein Flussdiagramm mit einigen Verfahrensschritten des erfindungsgemässen Verfahrens. In Fig. 1 ist eine bekannte Anordnung, bestehend aus einer Feldgeneratoreinheit 200, einem Sensorelement 300 und einer Verarbeitungs- und Steuereinheit 100, dargestellt. Die Verarbeitungs- und Steuereinheit 100 ist jeweils über Leitungen mit der Feldgeneratoreinheit 200 einerseits und dem Sensorelement 300 anderseits verbunden. Während sich die Feldgeneratoreinheit 200 vorzugsweise an einem bekannten Ort befindet - was bedeutet, dass die Koordinaten x, y, z, inkl. der Ausrichtung im Koordinatensystem, bekannt sind - kann das Sensorelement 300 beliebig bewegt werden bzw. eine beliebige Position und Orientierung einnehmen. Es wird darauf hingewiesen, dass es denkbar ist, wie bereits aus der WO 97/36192 bekannt ist, dass das Sensorelement 300 ortsfest und die Feldgeneratoreinheit 200 frei, d.h. im Rahmen der zur Verfügung gestellten

Verbindungsleitung zur Verarbeitungs- und Steuereinheit 100, beweglich ist. Des weiteren ist auch ohne weiteres denkbar, dass die Verarbeitungs- und Steuereinheit 100 in mehreren Funktionseinheiten realisiert ist, wie zum Beispiel, dass die Steuereinheit zur Steuerung der

Feldgeneratoreinheit 200 in einem Funktionsblock und die Verarbeitungseinheit, in der die eigentliche Positionsberechnung der Position des Sensorelementes 300 vorgenommen wird, in einem anderen Funktionsblock realisiert ist. Diese Abwandlungen von der in Fig. 1 dargestellten Anordnung haben keinen Einfluss auf die Anwendbarkeit des erfindungsgemässen Verfahrens. Das Gleiche gilt auch für Ausführungsformen, bei denen mehrere Feldgeneratoren an unterschiedlichen Stellen vorgesehen sind, wie dies beispielsweise bei der Lehre gemäss WO 97/36192 der Fall ist.

Mit 400 ist in schematischer Darstellung ein elektrisch leitendes Objekt stellvertretend für diejenigen Gegenstände dargestellt, welche die magnetische Ortung des Sensorelementes 300 stört, indem im Objekt 400 Wirbelströme 420 erzeugt werden, aufgrund derer ein dem Wechselfeld 210 überlagerten Stδrfeld 410 entsteht.

Bevor auf das erfindungsgemässe Verfahren weiter beschrieben wird, werden im folgenden zunächst allgemeine Zusammenhänge bzw. Vorgehensweisen bei der magnetfeldbasierten Ortung erläutert.

Wie bereits erwähnt wurde, wird bei der magnetfeldbasierten Ortung, welche auch etwa als magnetische Ortung bezeichnet wird, die Position und/oder Ausrichtung von einem oder mehreren Sensorelementen 300 relativ zu einem oder mehreren Feldgeneratoreinheiten 200 bestimmt. Die Position r_s und

Orientierung n_s der Sensorelemente 5, kann durch Lösung des folgenden Gleichungssystems bestimmt werden, wenn man voraussetzt, dass Position r_Gj und Orientierung n_0J der Feldgeneratoreinheiten G_} bekannt sind:

F_i] = F(r_Sι,n_s ,r_G , n_G] ) (1)

Mit i ist das i-te Sensorelement und mit j die j-te Feldgeneratoreinheit gemeint. F ist dabei eine vom Magnetfeld abhängige Messfunktion von meistens einer

Komponente des Magnetfelds B(x,y,z,t) (z.B. die induzierte Spannung in einer Sensorspule) . F kann natürlich auch eine Funktion von vielen zusammengebauten Sensoren in einem Sensorelement sein, welches mehrere oder alle Komponenten gleichzeitig misst. Nach Art der Lösung dieses Gleichungssystems können magnetische Positionssysteme in zwei Klassen unterteilt werden:

I. Das Gleichungssystem wird invertiert, d.h. die Sensorelementpositionen können aus den gemessenen Magnetfeldern berechnet werden:

r_Sι = f_r(F_¥) und n_Sι = f_n(F_ϋ) (2>

Da die Inversion des Gleichungssystems nur in sehr speziellen Fällen möglich ist, kann man durch Näherung versuchen, die Feldgleichungen in eine invertierbare Form zu bringen.

II. Das Gleichungssystem wird durch Optimierung gelöst, d.h. die Sensorelementpositionen werden so lange variiert, bis die nach Gleichung 1 berechneten Werte F_y am besten mit den gemessenen Werten E_(ϊ übereinstimmen. Eine mögliche Methode wäre ein Chi²-Fit nach Levenberg- Marquardt. Dabei werden die Sensorpositionen r_s und n_s, so lange variiert bis

j K^ g )

minimal wird. Für weiterführende Angaben zur Methode nach Levenberg-Marquardt wird stellvertretend auf die Druckschrift mit dem Titel „Numerical Recipies in C" (W. H. Press, S. A. Teukolsky, W. T. Vetterling und B. P.

Flannery; Cambridge university Press; 1994) verwiesen. Es ist auch eine Kombination beider Lδsungsansätze möglich.

Da die vom Sensorelement Si gemessene Grδsse F₀ nur von der relativen Lage vom Sensorelement Si und der

Feldgeneratoreinheit G_j abhängt, sind die Rollen von .Sensorelement und Feldgeneratoreinheit in allen magnetischen Positionsmesssystemen austauschbar.

Werden zeitlich variierende magnetische Felder benützt, so erzeugen diese - wie bereits erwähnt wurde - in benachbarten, elektrisch leitenden Objekten 400 Wirbelstrδme 420. Diese führen zu Verzerrungen des ursprünglichen magnetischen Wechselfeldes 210 und damit zu systematischen Fehlern bei der Positionsbestimmung. Das bedeutet, dass, wenn Position und Ausrichtung von Sensorelementen im verzerrten Wechselfeld so bestimmt werden, als wäre kein elektrisch leitendes Objekt 400 zugegen, die gewonnenen Werte systematisch verfälscht sind.

Um das Verfahren der magnetischen Ortung auch in der Nachbarschaft von elektrisch leitenden Objekten 400 frei von durch diese hervorgerufenen Fehlern anwenden zu können, werden erfindungsgemass die Wechselfeldverzerrungen und ihr Einfluss auf die Bestimmung von Sensorelementposition und Sensorelementausrichtung bestimmt. Damit können die auftretenden systematischen Fehler korrigiert werden, womit die Genauigkeit der Position und/oder der Ausrichtung erheblich verbessert werden kann.

Bei einer vorgegebenen Messanordnung hat man feststellen können, dass ein Messfehler von 4cm mit dem erfindungsgemässen Verfahren auf weniger als 1,5mm reduziert werden konnte.

Grundsätzlich lassen sich diese Korrekturen auch mit Hilfe der Technik der finiten Elemente und den Gleichungen der Elektrodynamik finden. Die bevorzugte Ausführungsform des erfindungsgemässen Verfahrens zeichnet sich gegenüber der Methode der finiten Elemente weiter dadurch aus, dass eine massive Reduktion der Positionsberechung erreicht wurde, da vieles im Sinne einer Systemkaiibration im voraus berechnet werden kann.

Im folgenden wird das erfindungsgemässe Verfahren erläutert, wobei der Einfachheit halber zunächst davon ausgegangen wird, dass das Objekt 400 aus einer elektrisch leitenden Platte, d.h. aus einer flachen und begrenzten Fläche, besteht. Objekte 400, die in der Richtung einer gedachten Linie von der Feldgeneratoreinheit 200 zum Objekt 400 eine relevante Ausdehnung (Tiefe) besitzen, können ebenfalls mit dem erfindungsgemässen Verfahren behandelt werden. Für diesen Zweck wird die der Feldgeneratoreinheit 200 zugewandte Seite durch eine Vielflächen-Struktur angenähert. Dies ist zulässig, da die Wirbelstrδme 420 nur geringfügig in die Oberfläche eindringen. Die Tiefe eines dreidimensionalen Objektes 400 ist deshalb nicht relevant. Bei der mathematischen Bestimmung des Stδrfeldes 410 wird daher das Objekt 400 durch eine Vielflächenstruktur im vorstehend genannten Sinne angenähert.

Im folgenden werden Überlegungen betreffend die Berechnung von Feldverzerrungen bei einer leitende Platte dargestellt. Die daraus resultierenden Ergebnisse können analog bei allgemeineren Objektformen angestellt werden. Befindet sich in einem zeitlich veränderlichen Magnetfeld B₀(x,y,z,t) ein elektrisch leitendes Objekt 400, so werden Wirbelstrδme 410 (Fig. 1) in die Oberfläche des Objektes 400 induziert. Diese Wirbelstrδme 410 verursachen ein weiteres Magnetfeld B'(x,y,z,t) , welches dem ursprünglichen Feld B₀(x_}y,z,t) überlagert ist und zu einem resultierenden Feld Bj_teiXt »²»*) führt. B^_s(x,y,z,t) ist gegenüber dem Feld B₀(x,y,z,t) verzerrt. Um dieses verzerrte Feld berechnen zu können, ist es notwendig, das induzierte Wechselfeld B'(x,y,z,i) zu kennen. Es lässt sich aus dem Biot-Savart Gesetz der Elektrodynamik (Gleichung 4) ein Feld B^x_^y^i) berechnen, welches gut genug B'(x,y,z,t) beschreibt, falls der örtliche und zeitliche Verlauf der Wirbelstrδme 410 im Objekt 400 in N verschiedenen punktfδrmigen Stromelementen bekannt ist:

wobei P_R = (x,y,z) ein Punkt im Raum ist und der Vektor F vom

Stromelement zum Punkt P_R zeigt. Vorfaktoren können, falls notwendig, eingeführt werden und/oder können im Betrag des Faktors Δs(t) enthalten sein. Eine detailliertere Berechnung von B_x(x,y,z,t) durch Einführung von Längsstrδmen oder Flächenstrδmen etc. anstatt Punktstrδmen ist möglich.

Dadurch würde die Schreibweise von Gleichung 4 allerdings geändert werden. In den meisten Fällen kann die Gleichung 4 hingegen so übernommen werden, wie sie oben aufgeführt ist, falls N gross genug gewählt wird.

Die Feldverzerrungen werden also in zwei Schritten berechnet. Der erste Schritt ist die Bestimmung der Wirbelstrδme 410 und der zweite Schritt ist die Berechnung eines durch die Wirbelstrδme 410 erzeugten Stδrfelds

B_\(.x_>y_>z ) ι welches gut genug das Störfeld B'(x,y,z,t) beschreibt.

Fig. 2 zeigt das Objekt 400, das für die Bestimmung des Stδrfeldes 410 in eine Vielflächenstruktur, bestehend aus beliebig vielen Segmenten, aufgeteilt ist. Ausgehend von dieser Aufteilung und diversen weiteren Annahmen werden zunächst die Wirbelstrδme berechnet.

Die Wirbelstrδme fHessen an der Oberfläche des Objektes 400 mit einer für die Theorie "unwesentlichen"

Eindringtiefe. Um das oben erwähnte Stδrfeld B'(x,y,z,t) hinreichend genau zu berechnen, reicht es aus, den zeitlichen Stromverlauf in einigen Punkten auf der Oberfläche des Objektes 400 zu kennen. Die Anzahl der Punkte hängt von der verlangten Genauigkeit ab. Die Wirbelstrδme werden also in Punkten berechnet, welche auf oder in der Nähe der Oberfläche des Objektes liegen.

In einem ersten Schritt wird das Objekt in N beliebig geformte Segmente unterteilt, die sinnvollerweise (aber nicht notwendigerweise) das ganze Objekt abdecken. Die Segmente werden nachfolgend zur eindeutigen Unterscheidung mit S₍.{0 < z' < N-l} bezeichnet, wobei i als Index verwendet wird. In einem zweiten Schritt wird pro Segment ein Stützpunkt P, gewählt. Es ist sinnvoll aber nicht zwingend, gleich viele Stützpunkte wie Segmente zu definieren und diese eindeutig zu den Segmenten zu zuordnen. Im folgenden wird der Einfachheit halber davon ausgegangen, dass N Segmente S, mit je einem eindeutig zugeordneten Stützpunkt P_t definiert werden. Die Stromdichte l,(t) im Stützpunkt P_t eines jeden Segments S, wird mittels der folgenden Formel berechnet:

ΛΓ ϊ_i(t) = ∑-i!_j(t)mitj o i (5) =0

wobei ?_u(f) die Stromdichte des Wirbelstroms I_y(t) ist, welche durch die Flussänderung des Feldes von B₀(x,y,z,t) im Segment S,. verursacht wird und durch den Stützpunkt P₍ oder in der Umgebung des Stützpunktes P_; fliesst. Die Berechnung der einzelnen Wirbelstrδme I_tj(t) ist im nächsten Abschnitt beschrieben. Zunächst gilt:

mit λ als Richtungsvektor (oder eine dazu fast kolineare Richtung) der Stromlinie durch den Stützpunkt P_t im

Stützpunkt P₍ und mit A_s die Querschnittsfläche der Stromline ist, wobei:

A_s = π -r²(h) (5b) mit r = Radius der kreisrunden Querschnittsfläche; h = Eindringtiefe;

Sind nun die Stromdichten T^t) bekannt, so ist auch

als Störfeld verursacht durch B₀(x,y,z,t) berechenbar - dazu kann ζ(t) - -4(S,.) direkt in die Gleichung 4 eingesetzt werden, wobei -4(5,.) die Fläche des Segments S_t ist. B (x,y,z_tt) kann in den meisten Fällen als B'(x,y,z,t) betrachtet werden. Als Effekt zweiter Ordnung kann an dieser Stelle B_l(x,y,z,t) als Urfeld eingesetzt werden, um wiederum Wirbelströme für ein zweites Stδrfeld B₂(x,y,z,t) zu berechnen (Einfluss der Wirbelstrδme auf einander) , welches mit B₀(x,y,z,t) und B_{(x,y,z,t) überlagert ist. In zweiter Näherung wäre also B'(x,y,z,t) gleich der Summe aus B_x(x,y,z,t) und B₂ x,y,z,t) - diese iterative Vorgehensweise kann für Effekte beliebiger Ordnung fortgesetzt werden. Es ist jedoch in den meisten Anwendungen der Effekt erster Ordnung hinreichend genau.

Ein einzelner Wirbelstrom I_y(t) ist eine Stromlinie, die durch den Stützpunkt P_t fliesst und durch die zeitliche Flussänderung des Feldes durch den Stützpunkt P_} verursacht wird. Um I_ϋ t) zu berechnen sind die Induktivität L_s , sein Ohm¹ scher Widerstand R_(ϊ und die zeitliche Flussänderung dΦ . —- notwendig. Sind diese Grδssen bekannt, so ist I_y(i) durch die Lösung der Differentialgleichung

gegeben. In vielen Fällen mag B₀(x,y,z,t) zeitlich periodisch sein oder sogar harmonisch oszillieren, dies ist aber für die Gültigkeit der erfindungsgemässen Methode keine Notwendigkeit .

Die Induktivität L_ή und der Ohm'sehe Widerstand R_y sind durch die geometrische Form des Wirbelstroms I_g(t) gegeben dΦ . und die Flussänderung -—— durch das Feld B₀(x,y,z,t) an der dt

Stelle P_j zusammen mit der Fläche des Segments S_j . In den folgenden Schritten wird zuerst die Form des Wirbelstroms beschrieben, und anschliessend daraus die Induktivität L_tj und der Ohm' sehe Widerstand R_y berechnet ,

Man stelle sich vor, dass eine einzelne Magnetfeldlinie B durch eine kleine Fläche dA um einen Punkt P auf dem Objekt dringt. In diesem Fall wären die induzierten Stromlinien in der Nähe vom Punkt P kreisförmig und am Rande des Objektes 400 würden sie der Begrenzung folgen, also die Form der Umrandung des Objekts haben. Die Form eines beliebigen Wirbelstroms ist ein Strom entlang einer Höhenlinie einer Fläche, welche die Potentialgleichung

erfüllt, wobei <p(x,y) das Potential darstellt. Die Randbedingungen zur eindeutigen Lösung von Gleichung 7 (ermitteln von φ(x,y) ) sind aus Figur 2 zu entnehmen (nämlich φ„ am Rand des Objektes und φ_x im Punkt P_j, φ₀ <> φ_x) . Diese Potentialgleichung ist am besten numerisch zu lösen. Mit der Form des Wirbelstroms und der Eindringtiefe h des Stroms in das Material ist die einzelne Stromlinie als Leiterschleife mit einem ringförmigen

Materialquerschnitt mit dem Durchmesser der Eindringtiefe betrachtbar (andere sinnvolle

Materialquerschnittsgeometrien sind natürlich denkbar, ändern aber an den Berechnungen nichts wesentliches) .

Der Ohm¹ sehe Widerstand der Leiterschleife wird somit zu

mit R = Ohmscher Widerstand [ Ω ]

1 = Länge der Leiterschleife [ m ] h = Eindringtiefe des Stroms [ m ] p = spezifischer elektrischer Widerstand des Materials [ Ωm ]

Die Induktivität der Leiterschleife ist gegeben durch

2W

1 = (9)

und kann numerisch berechnet werden, wobei

B²

W = f -=-^dV (10)

J II.

Raum * 2^■*μ*(_{> die gespeicherte Energie im vom Leiter erzeugten Magnetfeld bedeutet, falls ein Strom i im Leiter fliesst. Es gibt "beliebig¹¹ viele andere Näherungsformein, welche Gleichungen (8) , (9) und (10) ersetzen könnten und ähnliche Resultate liefern.

Für die Berechnung des Flusses Φ,(t) aus dem Feld B₀(x,y,z,t) gilt:

Φ.( = 5₀(x,; ,z,t . (11)

wobei B₀(x,y,z,t) das ungestörte Feld an der Stelle P_} ist und A_j die Flächennormale des Segments S_j ist, mit dem Betrag der Fläche des betreffenden Segments. Formel 11 ist eine Näherungsformel für die allgemein gültige Formel

Φ _j(t) = ^B₀(x,y,z,t)dA (11a)

und darf angewendet werden, wenn B über A_j hinreichend homogen ist (so z.B. für kleine Flächen A_j) . An dieser Stelle soll erwähnt werden, dass die Fläche A_} des

Segments S_} nicht in allen Fällen vollständig innerhalb der

Leiterschleife I_y(t) liegt - Korrekturen könnten diesbezüglich angebracht werden, sind aber in der Regel nicht notwendig.

Die teilweise intensiven Berechnung der Induktivitäten L_i} und der Ohm' sehen Widerstände R_y können mittels der Formeln (8) , (9) und (10) im voraus berechnet werden, da diese allein von der Geometrie und des Materials des Objektes abhängig sind - dies im Sinne einer Systemkalibration. Einsetzen und Lösen der Formeln (4) und (5) , um B_γ(x,y,z,f) zu berechnen, können zu einem anderen Zeitpunkt erfolgen, insbesondere wenn das Erregerfeld B₀(x,y,z,t) bekannt ist. Iterationen, wie z.B. Einsetzen des

Feldes

zur Berechnung eines Feldes B₂ x,y,z,t) usw. , sind möglich. Solche Iterationen könnten aber auch im voraus erfolgen, indem man sie in die Induktivitäten L_tJ einbindet. Dies macht aber erst einen Sinn, wenn regelmässig Korrekturen höherer Ordnung benötigt werden.

In der Praxis werden bei der magnetischen Ortung Lage und Ausrichtung eines oder mehrerer Sensorelemente 300 (Fig. 1) in einem Magnetfeld bestimmt, das von einem oder mehreren Feldgeneratoreinheiten 200 erzeugt wird. In dem verwendeten Koordinatensystem ist die Position der Feldgeneratoreinheit 200 oder den Feldgeneratoreinheiten bekannt. Im Fall von magnetischen Wechselfeldern erzeugen benachbarte elektrisch leitende Objekte 400 Feldverzerrungen durch in den Objekten 400 induzierte Wirbelstrδme 420. Das erfindungsgemässe Verfahren zur Korrektur dieser Verzerrungen, dessen theoretische Grundlagen oben angegeben wurden, wird folgendermassen angewendet: Die Lage der elektrisch leitenden Objekte 400 im oben erwähnten Koordinatensystem ist bekannt oder wird durch Vermessung bestimmt. Die Objektkoordinaten werden in einem Computerprogramm, das zur Berechnung der Wirbelstrδme 420 und der daraus resultierenden Feldverzerrungen benutzt wird, derartig eingegeben, dass die in den oben angegebenen Formeln verwendeten Ortskoordinaten im durch die Feldgeneratoreinheit 200 definierten Koordinatensystem definiert sind. Mit dem Computerprogramm wird dann das durch die Wirbelstrδme 420 erzeugte Stδrfeld berechnet. Bei Berücksichtigung der Wirbelstrδme 420 ändert sich das Gleichungssystem 1 wie folgt:

P

Fy = F(r_Sι,n_Sι,r_Gj,n_Gj ) + ∑F£(r_Sι,n_Sι,r_G n_G]) (12)

*=ι

ι wobei F_y die durch den Wirbelstrδme 420 des Objektes k erzeugte Störung darstellt. P ist die Anzahl der Objekte. Wie diese Korrektur angewendet wird, hängt von der Art des magnetischen Positionsmesssystems ab.

I. Bei Systemen, die auf Gleichung 2 beruhen, werden die Messwerte iterativ korrigiert, d.h. man berechnet zuerst die ungestörte Lösung nach Gleichung 2. Mit der gefundenen Position des Sensorelementes 300 kann der Korrekturterm F¹ berechnet und von den Messungen F_y ^M abgezogen werden. Mit den korrigierten Messungen wird wieder eine Position berechnet. Dieser Algorithmus wird fortgesetzt, bis die Variation der berechneten Positionen unter gewissen Toleranzschwellen liegen.

II. Bei Systemen, die auf Gleichung 3 beruhen, muss der Lδsungsalgorithmus nicht geändert werden. In der Chi-Summe werden anstatt des Modells zur Berechnung der wirbelstromfreien Magnetfeldes F_y nach Gleichung 1 das

Magnetfeld mit Wirbelstromkorrekturen nach Gleichung 12 benutzt. III. Unter bestimmten Voraussetzungen kann es auch möglich sein, das GleichungsSystem 12 zu invertieren, was dann zu einer Lösung entsprechend Gleichung 2 führt.

Fig. 3 zeigt, in vereinfachter Darstellung, ein

Struktogramm eines nach dem erfindungsgemässen Verfahren arbeitenden Computerprogramms. Die einzelnen Verarbeitungsschritte wurden bereits anhand der Fig. 1 und 2 ausführlich erläutert.

Das erfindungsgemässe Verfahren kann auch vorzüglich für

Objekte mit Öffnungen (Löcher) angewendet werden, wobei die Anzahl L der Öffnungen beliebig sein kann. Dabei muss für die schon beschriebene Lδsungsmethode zunächst die Randbedingung der Potentialgleichung (7) an den Rändern der Öffnungen gleich dem Potential φ₀ am Rande des Objektes sein. Des weiteren kommen zu der beschriebenen Methode weitere N mal L (N=Anzahl Stützpunkte und L=Anzahl Offnungen) Stromlinien I_ik hinzu, die in der Summe (5) noch zu addieren sind (k wandert von 1 bis L) .

Die zusätzlichen Wirbelstromlinien I_ik sind einzeln analog zu den Wirbelstromlinien I_i:j berechenbar, d.h. Lösen der Potentialgleichung (7) für die Form des Stroms und berechnung der Induktivität und des Widerstands nach Gleichungen (8) und (9) . Bei der Potentialgleichung (7) ist allerdings zu beachten, dass die Randbedingung nicht "φj, im Punkt Pj" lautet, sondern "§_x am Rand der Öffnung k" . Bei grossen Öffnungen ist eventuell Formel 11a für die Berechnung des Flusses zu verwenden anstatt die Näherungsformel 11.

Einzelne Leiterschleifen können auch mit dieser Methode berechnet werden, da die vorher erwähnten Öffnungen beliebig nahe an der Umrandung der zu berechnenden Objekte ausgedehnt werden kann. Das einfachste Beispiel ist ein Kreisring, welches als Scheibe mit einer fast gleich grossen Öffnung betrachtet werden kann: In diesem Beispiel sind die Stromlinien I_i;i vernachlässigbar (die Stützpunkte könnten weggelassen werden) und es gibt nur ein I_ik dessen Form durch den Kreisring gegeben ist. Falls die Stützpunkte weggelassen werden, ist das Feld B_x durch das Linienintegral über die Gleichung 4 zu bestimmen.

Ein weiterer Aspekt besteht darin, dass Stδreinflüsse von unbekannten Objekten abgeschirmt werden, indem man zwischen Feldgeneratoreinheit und Objekt eine leitende Platte vorsieht, wobei deren Grδsse und Form als auch deren Position bekannt sind. Damit müssen zwar die Feldverzerrungen dieser Platte berücksichtigt werden, doch alle anderen elektrisch leitenden Objekte, welche sich auf der in Bezug auf die Feldgeneratoreinheit anderen Seite der Platte befinden, können wegen der Abschirmung unberücksichtigt bleiben.

Claims

Patentansprüche:

1. Verfahren zur Bestimmung der Position eines Sensorelementes (300) , mit Hilfe dessen ein von mindestens einer Generatoreinheit (200) ausgesendetes magnetisches Wechselfeld (210) gemessen wird, wobei aufgrund eines im Sensorelement (300) empfangenen Signals die Position des Sensorelementes (300) bestimmt wird, dadurch gekennzeichnet,

- dass, vorzugsweise in erster Näherung, Stδrfeider (410) berechnet werden, welche aufgrund von in elektrisch leitenden Objekten (400) erzeugten Wirbelströmen (420) entstehen, und - dass die Position, welche ausgehend vom im

Sensorelement (300) empfangenen Signal bestimmbar ist, aufgrund der berechneten Stδrfelder (410) korrigiert wird.

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet,

- dass die Wirbelstrδme (420) im Objekt (400) , ausgehend vom Wechselfeld (210) , berechnet werden und

- dass die Stδrfelder (410) , ausgehend von den berechneten Wirbelstrδme (420) berechnet werden.

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass zur weiteren Verbesserung der Positionsbestimmung mindestens eine weitere Iteration durchgeführt wird, indem

- weitere Wirbelstrδme im Objekt (400) , ausgehend von berechneten Stδrfeidern (410) , berechnet werden und

- weitere Stδrfelder, ausgehend von den weiteren Wirbelströmen, berechnet werden.

4. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, dass Position und Form der Objekte (400) bestimmt werden und dass widerstände nach der Gleichung

und Induktivitäten nach der Gleichung

2W

L =

im Objekt (400) berechnet werden.

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass zur weiteren Verbesserung der Induktivitäten eine weitere

Iteration durchgeführt wird, indem

- weitere Wirbelstrδme im Objekt (400) , ausgehend von berechneten Stδrfeidern (410) , berechnet werden,

- weiter Induktivitäten, ausgehend von den weiteren Wirbelströmen, bestimmt werden.

6. Verfahren nach Anspruch 4 und 5, dadurch gekennzeichnet, dass die Bestimmung der Objekte (400) , die Bestimmung der Widerstände und Induktivitäten im Objekt (400) im Sinne einer Systemkalibrierung im voraus, d.h. vor den das

Wechselfeld (210) berücksichtigenden Berechnungen, bestimmt werde .

7. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, dass zur Bestimmung der Wirbelstrδme wie folgt vorgegangen wird:

- die Objekte (400) werden in Segmente (S_£) und Stützstellen (Pi) aufgeteilt;

- in den Stützpunkten (Pi) wird die Stromdichte (ii) bestimmt;

- aus den Wirbelstrδ en (I_i:j) werden die Stromdichten (ii) bestimmt.

8. Anwendung des Verfahrens nach einem der Ansprüche 1 bis 7 für die magnetfeldorientierte Ortung bei „Cyber Spaceⁿ-

Anwendunge .

9. Vorrichtung zur Durchführung des Verfahrens nach Anspruch 1 bis 7, dadurch gekennzeichnet, dass mindestens eine Feldgeneratoreinheit (200) , mindestens ein Sensorelement (300) und eine Verarbeitungs- und Steuereinheit (100) vorgesehen sind, wobei die Feldgeneratoreinheit (200) und das Sensorelement (300) mit der Verarbeitungs- und Steuereinheit (100) verbunden sind.

10. Vorrichtung nach Anspruch 9, dadurch gekennzeichnet, dass mindestens ein elektrisch leitendes Objekt zum

Abschirmen der Feldgeneratoreinheit (200) vorgesehen sind.

11. Computerprogrammprodukt, das in den internen Speicher eines digitalen Computers geladen werden kann und Softwarecodeabschnitte umfasst, mit denen die Schritte gemäss den Ansprüchen 1 bis 6 ausgeführt werden, wenn das Produkt auf einem Computer läuft.