DE102019210336A1 - Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen - Google Patents

Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen Download PDF

Info

Publication number
DE102019210336A1
DE102019210336A1 DE102019210336.7A DE102019210336A DE102019210336A1 DE 102019210336 A1 DE102019210336 A1 DE 102019210336A1 DE 102019210336 A DE102019210336 A DE 102019210336A DE 102019210336 A1 DE102019210336 A1 DE 102019210336A1
Authority
DE
Germany
Prior art keywords
strain
bearing
sensor
estimation method
sensor probe
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
DE102019210336.7A
Other languages
English (en)
Inventor
Defeng Lang
Hendrik Anne Mol
Stijn van Eesbeek
Henricus van Genuchten
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SKF AB
Original Assignee
SKF AB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SKF AB filed Critical SKF AB
Priority to DE102019210336.7A priority Critical patent/DE102019210336A1/de
Priority to US16/916,235 priority patent/US11428590B2/en
Priority to CN202010654974.7A priority patent/CN112284584A/zh
Publication of DE102019210336A1 publication Critical patent/DE102019210336A1/de
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • G01L1/24Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0009Force sensors associated with a bearing
    • G01L5/0019Force sensors associated with a bearing by using strain gages, piezoelectric, piezo-resistive or other ohmic-resistance based sensors
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C19/00Bearings with rolling contact, for exclusively rotary movement
    • F16C19/02Bearings with rolling contact, for exclusively rotary movement with bearing balls essentially of the same size in one or more circular rows
    • F16C19/14Bearings with rolling contact, for exclusively rotary movement with bearing balls essentially of the same size in one or more circular rows for both radial and axial load
    • F16C19/16Bearings with rolling contact, for exclusively rotary movement with bearing balls essentially of the same size in one or more circular rows for both radial and axial load with a single row of balls
    • F16C19/163Bearings with rolling contact, for exclusively rotary movement with bearing balls essentially of the same size in one or more circular rows for both radial and axial load with a single row of balls with angular contact
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C19/00Bearings with rolling contact, for exclusively rotary movement
    • F16C19/52Bearings with rolling contact, for exclusively rotary movement with devices affected by abnormal or undesired conditions
    • F16C19/522Bearings with rolling contact, for exclusively rotary movement with devices affected by abnormal or undesired conditions related to load on the bearing, e.g. bearings with load sensors or means to protect the bearing against overload
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C41/00Other accessories, e.g. devices integrated in the bearing not relating to the bearing function as such
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C43/00Assembling bearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C2233/00Monitoring condition, e.g. temperature, load, vibration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0009Force sensors associated with a bearing
    • G01L5/0023Force sensors associated with a bearing by using magnetic sensors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

Belastungsschätzverfahren zum Schätzen einer Lagerbelastung, das in einer elektronischen Steuereinheit (24) einer Lageranordnung implementiert ist, die ein Lager (12) hat, das einen ersten Ring (14), einen zweiten Ring (16) und eine Reihe von Wälzkörpern (18) aufweist, die dazwischen angeordnet sind, wobei einer der Ringe (14) mit zumindest einem Dehnungssensorfühler (20a, 20b, 20c) ausgestattet ist, der dazu ausgebildet ist, Dehnungssignale an die elektronische Steuereinheit (24) zu übertragen.
Das Verfahren weist auf:
- Extrahieren einer Forminformation der Kurvenform des Dehnungssignals unter Verwendung von harmonischen Komponenten (u, v) des Dehnungssignals; und
- Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion der harmonischen Komponenten (u, v) des Dehnungssignals.

Description

  • Die vorliegende Erfindung bezieht sich auf das Gebiet der Lagerbelastungsschätzverfahren.
  • Um Lagerbelastungen zu erfassen, ist es bekannt, die Wälzkörperbelastung aus dem Dehnungssignal an der Sensorposition für jeden vorbeilaufenden Wälzkörper zu berechnen. In einem zweiten Schritt wird die Lagerbelastung aus den Wälzkörperbelastungen rekonstruiert.
  • Die Dehnung hängt von der Größe und Richtung der Wälzkörperbelastung und der Position und Orientierung des Bereichs, an dem der Sensor angeordnet ist, bezüglich der Wälzkörperbelastungsvektoren ab.
  • Andere Lagerbelastungsschätzverfahren zum Schätzen der Belastung berechnen eine Messung einer Signalamplitude. Bekannte Rechenverfahren verwenden Spitze-zu-Spitze-Verfahren (PKP), Effektivwertverfahren oder einen Durchschnitt eines absoluten Dehnungswerts.
  • Es ist auch bekannt, eine Belastung, die auf ein statisches, nicht-rotierendes Lager wirkt, mit einem Feld von Dehnungssensoren zu messen, die in den festen stationären Ring eines Lagers integriert sind.
  • Es kann Bezug auf die Dokumente US2013211742 - A1 , US 7,444,888 - B2 und WO 2011/128047 - A1 genommen werden, die Sensorlager, die mit Dehnungssensoren ausgestattet sind, zusammen mit Verfahren zum Bestimmen eines Belastungsvektors beschreiben, der auf ein Lager mit Wälzkörpern wirkt. Üblicherweise weisen die gemessenen Dehnungssignale für jeden Sensor eine Reihe von scharfen Spitzen auf, die den vorbeilaufenden Wälzkörpern entsprechen. Die Belastungsdetektion kann ein Analysieren, beispielsweise Mitteln über eine bestimmte Zeitdauer, des Signals über zumindest manche der Signalspitzen aufweisen, um eine verlässliche Information zu erhalten.
  • Jedoch können bekannte Verfahren eine falsche Rekonstruktion der aktuellen Gesamtlagerbelastung haben, da sich das gemessene Dehnungssignal aufgrund der Veränderung des Kontaktwinkels verändert. In anderen Worten verändert die Abweichung des Kontaktwinkels der Wälzkörper die Transferfunktion der Lagerbelastung und des Dehnungssignals.
  • Deswegen können Verfahren, die eine Spitze-zu-Spitze-Dehnung verwenden, um die Lagerbelastung zu berechnen, eine Schätzungsunsicherheit enthalten.
  • Um die Genauigkeit von Lagerbelastungsschätzverfahren zu verbessern, wurden globale und lokale Dehnungen verwendet, um die Lagerbelastungen zu schätzen.
  • Jedoch entfernt keines dieser Verfahren Fehler aufgrund einer Wälzkörperkontaktwinkelschwankung.
  • Es gibt somit einen Bedarf, Fehler aufgrund einer Wälzkörperkontaktwinkelschwankung in Lagerbelastungsschätzverfahren zu entfernen.
  • Das Ziel der vorliegenden Erfindung ist es dementsprechend, ein Lagerbelastungsschätzverfahren mit einer verbesserten Genauigkeit und Berücksichtigung der Schwankungen des Wälzkörperkontaktwinkels bereitzustellen.
  • Es ist eine besondere Aufgabe der vorliegenden Erfindung, ein Belastungsschätzverfahren zum Schätzen einer Lagerbelastung bereitzustellen, das in einer elektronischen Steuereinheit einer Lageranordnung implementiert ist, die ein Lager hat, das einen ersten Ring, einen zweiten Ring und eine Reihe von Wälzkörpern aufweist, die dazwischen angeordnet sind, wobei einer der Ringe mit zumindest einem Sensorfühler ausgestattet ist, der dazu ausgebildet ist, Dehnungssignale an die elektronische Steuereinheit zu übertragen.
  • Das Verfahren weist das Extrahieren einer Forminformation der Kurvenform des Dehnungssignals unter Verwendung von harmonischen Komponenten des Dehnungssignals auf.
  • Das Verfahren weist des Weiteren das Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion der harmonischen Komponenten des Dehnungssignals auf.
  • Belastungsschätzungen werden somit unter Verwendung der Kontaktwinkelparameter, die die Kontaktwinkelschwankung darstellen, berechnet. Der Kontaktwinkel wird aus den harmonischen Komponenten in dem Dehnungssignal geschätzt.
  • In einer Ausführungsform weist das Belastungsschätzverfahren auf:
    • - Extrahieren der Forminformation der Kurvenform des Dehnungssignals unter Verwendung von zumindest einer ersten und einer zweiten harmonischen Komponente des Dehnungssignals;
    • - Formulieren der ersten und der zweiten harmonischen Komponente der Wälzkörperkurvenform, die vorher extrahiert wurde, in Skalaren; und
    • - Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion des berechneten Verhältnisses der ersten und der zweiten Harmonischen, der Sensorfühlerpositionen und von Kalibrierungskonstanten.
  • Das Verhältnis der Energie der ersten und der zweiten harmonischen Komponente beschreibt die Beziehung zu der Kontaktwinkelschwankung.
  • In einer anderen Ausführungsform verwendet das Belastungsschätzverfahren mehr als einen Dehnungssensorfühler pro Wälzkörperzwischenraum in der Laufrichtung, beispielsweise zwei Dehnungssensorfühler, wodurch im Ortsraum statt im Zeitraum abgetastet wird. Die Sensorfühler sind umfänglich regelmäßig beabstandet.
  • Das Verfahren weist auf:
    • - Extrahieren der Forminformation der Kurvenform des Dehnungssignals unter Verwendung von zumindest einer ersten und einer zweiten harmonischen Komponente der Dehnungssignale von den regelmäßig beabstandeten Sensorfühlern;
    • - Formulieren der ersten und der zweiten harmonischen Komponente der Wälzkörperkurvenform, die vorher extrahiert wurde, in Skalaren; und
    • - Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion des berechneten Verhältnisses der ersten und der zweiten Harmonischen, der Sensorfühlerpositionen und der Kalibrierungskonstanten.
  • Beispielsweise verwendet das Belastungsschätzverfahren zwei Reihen von zumindest einem Sensorfühler pro Reihe zum Messen der Dehnung abhängig von dem Kontaktwinkel und Berechnen der Harmonischen der zwei Reihen, wodurch die Genauigkeit der Belastungsrekonstruktion an der Position des Wälzkörpers weiter verbessert wird.
  • Vorteilhafterweise verwendet das Belastungsschätzverfahren mehr als einen Dehnungssensorfühler pro Wälzkörperzwischenraum, der an einem der Ringe angeordnet ist, beispielsweise drei, fünf oder sieben, was die Messung von spezifischen harmonischen Komponenten des Dehnungssignals aufgrund der Wälzkörperbelastungsvektorgröße und Richtung mit korrekter Zurückweisung von Biegemoden niedriger Ordnung, unter denen der Ring auch beansprucht ist, wodurch die Genauigkeit der Belastungsrekonstruktion erhöht wird.
  • Jeder Dehnungssensorfühler kann beispielsweise einen direkt abgelagerten Dünnschichtdehnungsmessstreifen, der eine isolierende Beschichtung zwischen dem Ring und einem Schätzgitter verwendet, und ein Gitter aufweisen, das durch Rasterdruck und Laserglühen der metallischen Schicht des Sensorfühlers in ein funktionelles Gitter abgelagert wurde.
  • Andere Arten von Dehnungssensorfühlern können verwendet werden, beispielsweise ein Dickschichtdehnungsmessstreifen aus Stahl, der dünne Platten oder Strukturen aufweist, die an dem Ring an zwei Linien unter Verwendung von Schweißen (Laser, wie insbesondere in WO2011006523 A1 und WO2009056334 A1 beschrieben ist) angebracht und derart platziert und dimensioniert sind, dass sie aussagekräftig die erste und die zweite Harmonische des Dehnungsunterschieds an der Oberfläche der zwei Schweißlinien beobachten.
  • Alternativ weist jeder Dehnungssensorfühler einen glasfaserbasierten Dehnungsmessstreifen auf, der an einen der Ringe oder in einer Nut, die in dem Ring hergestellt ist, geklebt ist.
  • In einer anderen Ausführungsform weisen die Sensorfühler zumindest einen Weggeber in zumindest einer Reihe auf, um die relative Bewegung in radialer und umfänglicher Richtung aufgrund der Wälzkörperbelastung, die unter dem Hohlraum oder der Nut oder der Tasche, die anderweitig durch einen Dehnungssensor besetzt ist, durchläuft, zu beobachten. Die Bewegung wird in harmonischen Begriffen der Belastung und der Richtung der Belastung dargestellt und das dargestellte Berechnungsverfahren gilt für die Schätzung der Belastung.
  • Die vorliegende Erfindung und ihre Vorteile werden durch ein Studium der detaillierten Beschreibung von spezifischen Ausführungsformen besser verstanden werden, die als nicht-beschränkende Beispiele angegeben sind und durch die angehängten Zeichnungen dargestellt sind, in denen:
    • 1 und 2 schematische Halbschnittansichten einer Lageranordnung sind, die Lagerkräfte und Winkel zeigen;
    • 3 ein Flussdiagramm ist, das ein Lagerbelastungsschätzverfahren gemäß einer Ausführungsform der Erfindung behandelt; und
    • 4 verschiedene Graphen darstellt, die die Beziehung zwischen einer Lagerbelastungsschätzung unter Verwendung von verschiedenen Schätzverfahren und dem Kontaktwinkel der Wälzkörper zeigen.
  • Zunächst wird Bezug auf 1 und 2 genommen, die eine Ausführungsform einer Lageranordnung 10 mit einem Lager 12 darstellen, dessen Belastung zur Erfassung durch ein Belastungsschätzverfahren gemäß der Erfindung vorgesehen ist.
  • Das Lager 12 weist einen feststehenden Außenring 14 als einen ersten Ring, einen rotierenden Innenring 16 als einen zweiten Ring, eine Reihe von Wälzkörpern 18, wie beispielsweise Kugeln, die dazwischen angeordnet sind, und einen Käfig (nicht gezeigt) auf.
  • Das Lager 12 kann von der Kugellagerart sein.
  • Wie dargestellt ist, ist der Außenring 14 des Lagers 12 mit Dehnungssensorfühlern 20a, 20b, 20c ausgestattet. Die Dehnungssensorfühler 20a bis 20c können Dehnungsmessstreifen, beispielsweise Mikromessungen, sein, die an der radial äußeren Fläche des Außenrings 14 oder in einer definierten Nut oder Tasche (nicht gezeigt) befestigt sind.
  • Die Erfindung ist nicht auf dieses Beispiel beschränkt. In der Tat kann der Innen- oder Außenring des Lagers mit einer unterschiedlichen Anzahl von Dehnungssensorfühlern, beispielsweise zumindest einem Dehnungssensorfühler, ausgestattet sein.
  • Die Dehnungssensorfühler 20a bis 20c sind mit einer elektronischen Steuereinheit, ECU 24, verbunden, die dazu ausgebildet ist, Signale, die durch die Dehnungssensorfühler 20a-20c erhalten werden, wenn das Lager 12 läuft oder feststehend ist, zu empfangen und zu verarbeiten und eine Lagerbelastung zu schätzen.
  • Ein Belastungsschätzverfahren 30 ist in der elektronischen Steuereinheit 24 implementiert. Das Belastungsschätzverfahren 30 ist auf alle Lagerschätzvorrichtungen anwendbar, die zumindest einen Dehnungssensor haben.
  • Das Belastungsschätzverfahren 30 ist dazu ausgebildet, zusätzliche Signaleigenschaften, wie beispielsweise die Forminformation der Kurvenformen, zu verwenden, um die Genauigkeit der Kontaktwinkelschätzung zu verbessern.
  • In der Tat verwenden bekannte Verfahren, die Spitze-zu-Spitze-Verfahren (PKP) oder Effektivwertverfahren (RMS) verwenden, nur eine Amplitudenmessung des Signals und sind somit nicht genau, da sie keine Veränderungen des Kontaktwinkels berücksichtigen.
  • Der Kontaktwinkel hat einen Einfluss auf die Transferfunktion von der Wälzkörperbelastung zu dem Sensorsignal. Die zwei Grundsignaleigenschaften des Ausgangs der Transferfunktion, die durch den Kontaktwinkel beeinflusst werden, sind die Verstärkung und Form der Kurvenform, die an dem Sensor wahrgenommen wird, wenn ein Wälzkörper vorbeiläuft.
  • Das Ziel des Belastungsschätzverfahrens 30 ist es, die Abhängigkeit der Transferfunktion von dem Kontaktwinkel zu reduzieren.
  • Das Belastungsschätzverfahren 30 verwendet Kontaktwinkelparameter u, v, die die Skalardarstellungen der ersten und der zweiten Harmonischen sind, als Eingang, um Lagerbelastungen zu schätzen. Wie weiter beschrieben wird, werden diese Kontaktwinkelparameter u, v mit Kalibrierungsparametern kombiniert, um eine Schätzung von axialen und radialen Belastungen auf das Lager zu berechnen.
  • Die axiale Belastung FA auf das Lager gleicht der Summe der axialen Komponenten der Z Wälzkörperbelastungsvektoren Fc,j, wobei j der Index des Wälzkörpers ist. Axiale Komponenten können unter Verwendung des Kontaktwinkels aj der Kontaktbelastung gemäß der Gleichung unten geschrieben werden: F A = j = 1 Z s i n α j . F c , j
    Figure DE102019210336A1_0001
  • Die radiale Belastung in x- und y-Richtung kann als die Summe der x- und y-Projektionen der radialen Komponenten aller Wälzkörperbelastungsvektoren geschrieben werden. Sie sind durch die Winkelposition δi der Wälzkörper wie folgt definiert: F x = j = 1 Z c o s δ j . c o s α j . F c , j
    Figure DE102019210336A1_0002
    F y = j = 1 Z s i n δ j . c o s α j . F c , j
    Figure DE102019210336A1_0003
  • Die axiale Projektion X wird unter Verwendung des Kontaktwinkels γi an der Dehnungssensorposition βi durchgeführt.
  • Die axiale Belastungsschätzung F̂A kann unter Verwendung einer Funktion c wie folgt geschrieben werden: F ^ A = i = 1 N c . s i n γ i . u i j = 1 Z s i n α j . F c , j
    Figure DE102019210336A1_0004
  • Diese Annahme vermeidet die explizite Berechnung der Rollenbelastungen und Wälzkörperpositionen.
  • Gleichermaßen können Schätzungen der radialen Belastung in x- und y-Richtung unter Verwendung der Winkelposition der Dehnungssensoren βi wie folgt definiert werden: F ^ x = i = 1 N c . c o s β i . c o s γ i . u i
    Figure DE102019210336A1_0005
    F ^ y = i = 1 N c . s i n β i . c o s γ i . u i
    Figure DE102019210336A1_0006
  • In den Gleichungen 4 bis 6 ist ui ein Kontaktwinkelparameter, der aus Messungen abgeleitet ist. Die Sensorpositionen βi werden als bekannt angenommen und c stellt eine beliebige konstante Abbildung Dehnung-zu-Last dar.
  • Das Belastungsschätzverfahrens 30 weist einen ersten Schritt 32 auf, in dem eine Forminformation der Kurvenform unter Verwendung von zumindest einer ersten und einer zweiten harmonischen Komponente u, v des Signals extrahiert wird. U und v sind Kontaktwinkelparameter. Die Kombination dieser Parameter wird verwendet, um eine Information über die Amplitude sowie über die Form der Kurvenform bereitzustellen, da sich die Form mit Kontaktwinkelschwankungen ändert.
  • Die Forminformation ist als eine Energie der sinusförmigen Harmonischen quantifiziert. Mehrere Verfahren können verwendet werden, um die Harmonischen von dem Signal zu trennen, wie beispielsweise eine Kalman-Filterung, Vektorprojektion, Fouriertransformation. Diese Verfahren sind bekannt und werden nicht weiter beschrieben.
  • Bei Block 34 sind die erste und die zweite harmonische Komponente u, v der Kurvenform des Wälzkörpers 18, die bei Block 32 extrahiert wurde, in Skalaren ausgedrückt. Das Energieverhältnis der ersten und der zweiten harmonischen Komponente u, v ist somit berechnet. Das Verhältnis stellt die Kontaktwinkelschwankung dar.
  • Bei Block 36 wird der Kontaktwinkel des Wälzkörper 18 abhängig von dem bei Block 34 berechneten Verhältnis wie folgt geschätzt: c . s i n γ i = C ( u , v ) = c u i + c v i v i u i
    Figure DE102019210336A1_0007
  • Bei Block 38 wird die Lagerbelastungsschätzung als eine Polynomfunktion der ersten und der zweiten harmonischen Energieskalare berechnet.
  • Die axiale Belastungsschätzung wird durch Substituieren der Gleichung 7 in Gleichung 4 wie folgt berechnet: F ^ A = i = 1 N c u i . u i + c v i . v i
    Figure DE102019210336A1_0008
  • Die Konstanten cui und cvi sind abhängig von Dehnungsparametern unter Verwendung von Kalibrierungskonstanten k geschrieben: c u i = k u i 1 + k u i 2 . u i
    Figure DE102019210336A1_0009
    und c v i = k v i 1 + k v i 2 . v i
    Figure DE102019210336A1_0010
  • Indem die Gleichungen 9a und 9b in Gleichung 8 substituiert werden, wird die axiale Belastungsschätzung für eine einzelne Messung wie folgt geschrieben: F ^ A = i = 1 N k u i 1 . u i + k u i 2 . u i 2 + k v i 1 . v i + k v i 2 . v i 2
    Figure DE102019210336A1_0011
    F ^ A = [ U 1   U 2   V 1   V 2 ] = [ k A u 1 k A u 2 k A u 3 k A u 4 ]
    Figure DE102019210336A1_0012
    wobei:
    • Un, Vn Matrizen sind, die u und v harmonische Komponenten enthalten; und
    • KA ein 4 Nx1 Spaltenvektor ist.
  • Die Erfindung ist nicht auf die Zahl 4 und die Verwendung von zwei harmonischen Komponenten beschränkt. In der Tat kann die Belastungsschätzung unter Verwendung von höheren Ordnungen und/oder Harmonischen ohne Verlust der Gültigkeit der Berechnung berechnet werden.
  • Die Polynomfunktion der axialen Belastungsschätzung schließt somit den Kontaktwinkel aus, was die Genauigkeit des Belastungsschätzverfahrens verbessert.
  • Mit einer ähnlichen Beweisführung kann die Radialbelastungsschätzung basierend auf Dehnungsparametern und Kalibrierungskonstanten berechnet werden.
  • In einer allgemeinen Weise werden Belastungsschätzungen unter Verwendung der Kontaktwinkelparameter ui, vi, die die Kontaktwinkelschwankung darstellen, der Sensorpositionen βi und der Kalibrierungskonstanten k berechnet. Der Kontaktwinkel γi wird nicht direkt berechnet, sondern implizit aus dem Dehnungssignal geschätzt.
  • Dank des Verfahrens wird die Lagerbelastung mit einer höheren Genauigkeit im Vergleich zu bekannten Verfahren geschätzt, wie in 3 gesehen werden kann.
  • Dank des Belastungsschätzverfahrens werden Signalkomponenten verwendet, um den Einfluss des Kontaktwinkels aus der Gleichung für die Schätzung der Belastung herauszubekommen.
  • Harmonische Komponenten werden somit als eine Kontaktwinkelschwankung verwendet, um die Genauigkeit der Belastungsschätzung zu verbessern.
  • 4 stellt vier Graphen G1, G2, G3, G4 dar, die die Beziehung zwischen der Lagerbelastung und dem Kontaktwinkel der Wälzkörper darstellen.
  • Der Graph G1 in durchgezogener Linie stellt die Beziehung zwischen der wahren Lagerbelastung und dem Lagerkörperkontaktwinkel dar.
  • Der Graph G2 in gestrichelten Linien stellt die Beziehung zwischen der Lagerbelastung, die gemäß dem Verfahren vor von 2 geschätzt ist, und dem Lagerkörperkontaktwinkel dar.
  • Der Graph G3 in gestrichelten Linien stellt die Beziehung zwischen der Lagerbelastung, die gemäß einem Verfahren, das Spitze-zu-Spitze verwendet, geschätzt ist, und dem Lagerkörperkontaktwinkel dar.
  • Der Graph G4 in gepunkteten Linien stellt die Beziehung zwischen der Lagerbelastung, die gemäß einem Verfahren, das einen Effektivwert verwendet, geschätzt ist, und dem Lagerkörperkontaktwinkel dar.
  • Gemäß 4 kann gesehen werden, dass die geschätzte Lagerbelastung des Verfahrens gemäß der Erfindung besonders nahe an der realen Lagerbelastung ist. Der Einfluss der Kontaktwinkelschwankung ist somit aus dem Lagerbelastungsschätzverfahren gemäß der Erfindung entfernt.
  • ZITATE ENTHALTEN IN DER BESCHREIBUNG
  • Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
  • Zitierte Patentliteratur
    • US 2013211742 A1 [0006]
    • US 7444888 B2 [0006]
    • WO 2011128047 A1 [0006]

Claims (9)

  1. Belastungsschätzverfahren zum Schätzen einer Lagerbelastung, das in einer elektronischen Steuereinheit (24) einer Lageranordnung implementiert ist, die ein Lager (12) hat, das einen ersten Ring (14), einen zweiten Ring (16) und eine Reihe von Wälzkörpern (18) aufweist, die dazwischen angeordnet sind, wobei einer der Ringe (14) mit zumindest einem Dehnungssensorfühler (20a, 20b, 20c) ausgestattet ist, der dazu ausgebildet ist, Dehnungssignale an die elektronische Steuereinheit (24) zu übertragen, wobei das Verfahren aufweist: - Extrahieren einer Forminformation der Kurvenform des Dehnungssignals unter Verwendung von harmonischen Komponenten (u, v) des Dehnungssignals; und - Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion der harmonischen Komponenten (u, v) des Dehnungssignals.
  2. Belastungsschätzverfahren gemäß Anspruch 1, welches aufweist: - Extrahieren der Forminformation der Kurvenform des Dehnungssignals unter Verwendung einer ersten und einer zweiten harmonischen Komponente (u, v) des Dehnungssignals; - Formulieren der ersten und der zweiten harmonischen Komponente (u, v) der Kurvenform der Wälzkörper (18), die vorher extrahiert wurde, in Skalaren; und - Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion des berechneten Verhältnisses der ersten und der zweiten Harmonischen, der Sensorfühlerpositionen (βi) und von Kalibrierungskonstanten (k).
  3. Belastungsschätzverfahren gemäß Anspruch 1, welches mehr als einen Dehnungssensorfühler pro Wälzkörperzwischenraum in der Laufrichtung verwendet, wobei die Sensorfühler umfänglich regelmäßig beabstandet sind und wobei das Verfahren aufweist: - Extrahieren der Forminformation der Kurvenform des Dehnungssignals unter Verwendung von zumindest einer ersten und einer zweiten harmonischen Komponente (u, v) der Dehnungssignale von den regelmäßig beabstandeten Sensorfühlern; - Formulieren der ersten und der zweiten harmonischen Komponente (u, v) der Kurvenform der Wälzkörper (18), die vorher extrahiert wurde, in Skalaren; und - Berechnen einer Lagerbelastungsschätzung als eine Polynomfunktion des berechneten Verhältnisses der ersten und der zweiten Harmonischen, der Sensorfühlerpositionen (βi) und der Kalibrierungskonstanten (k).
  4. Belastungsschätzverfahren gemäß einem der vorhergehenden Ansprüche, welches zwei Reihen von zumindest einem Sensorfühler pro Reihe zum Messen der Dehnung abhängig von dem Kontaktwinkel und Berechnen der Harmonischen der zwei Reihen aufweist.
  5. Belastungsschätzverfahren gemäß einem der vorhergehenden Ansprüche, welches eine ungerade Anzahl von Dehnungssensoren, beispielsweise drei, fünf oder sieben Dehnungssensoren pro Wälzkörperzwischenraum, aufweist, wobei die Dehnungssensoren an einem der Ringe (14) angeordnet sind.
  6. Belastungsschätzverfahren gemäß einem der vorhergehenden Ansprüche, wobei jeder Dehnungssensorfühler kann beispielsweise einen direkt abgelagerten Dünnschichtdehnungsmessstreifen, der eine isolierende Beschichtung zwischen dem Ring und einem Schätzgitter verwendet, und ein Gitter aufweist, das durch Rasterdruck und Laserglühen der metallischen Schicht des Sensorfühlers in ein funktionelles Gitter abgelagert wurde.
  7. Belastungsschätzverfahren gemäß einem der Ansprüche 1 bis 5, wobei jeder Dehnungssensorfühler einen glasfaserbasierten Dehnungsmessstreifen aufweist, der an einen der Ringe (14) geklebt ist.
  8. Belastungsschätzverfahren gemäß einem der Ansprüche 1 bis 4, wobei die Sensorfühler zumindest einen Weggeber in zumindest einer Reihe aufweisen.
  9. Lageranordnung (10), welche ein Lager (12) aufweist, das einen ersten Ring (14), einen zweiten Ring (16) und eine Reihe von Wälzkörpern (18) hat, die dazwischen angeordnet sind, wobei einer der Ringe (14) mit zumindest einem Dehnungssensorfühler (20a, 20b, 20c) ausgestattet ist, wobei die Lageranordnung (10) des Weiteren eine elektronische Steuereinheit (24) zum Empfangen und Verarbeiten von Signalen aufweist, die durch den Dehnungssensorfühler (20a, 20b, 20c) erhalten werden, wobei die elektronische Steuereinheit (24) dazu ausgebildet ist, das Verfahren gemäß einem der vorhergehenden Ansprüche anzuwenden.
DE102019210336.7A 2019-07-12 2019-07-12 Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen Pending DE102019210336A1 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102019210336.7A DE102019210336A1 (de) 2019-07-12 2019-07-12 Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen
US16/916,235 US11428590B2 (en) 2019-07-12 2020-06-30 Method for estimating a bearing load using strain parameters to account for contact angle variation
CN202010654974.7A CN112284584A (zh) 2019-07-12 2020-07-09 用于使用考虑接触角变化的应变参数估算轴承载荷的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102019210336.7A DE102019210336A1 (de) 2019-07-12 2019-07-12 Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen

Publications (1)

Publication Number Publication Date
DE102019210336A1 true DE102019210336A1 (de) 2021-01-14

Family

ID=74092047

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102019210336.7A Pending DE102019210336A1 (de) 2019-07-12 2019-07-12 Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen

Country Status (3)

Country Link
US (1) US11428590B2 (de)
CN (1) CN112284584A (de)
DE (1) DE102019210336A1 (de)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102021211493A1 (de) 2021-10-12 2023-04-13 Aktiebolaget Skf Vorrichtung und Verfahren zum Abschätzen einer Last in einem Lager
CN114397279B (zh) * 2022-01-19 2023-07-18 天津大学 任意应变状态下二维材料和应变物体的应变状态检测方法

Family Cites Families (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6095879A (en) * 1998-05-14 2000-08-01 Mccarthy; Peter T. Methods for creating consistent large scale blade deflections
JP2004061151A (ja) * 2002-07-25 2004-02-26 Nsk Ltd 軸受装置の接触角測定方法及び装置
CN100442041C (zh) * 2003-05-22 2008-12-10 日本精工株式会社 用于滚动轴承单元的载荷测量装置以及载荷测量滚动轴承单元
NL1023948C2 (nl) 2003-07-18 2005-01-19 Skf Ab Werkwijze en sensoropstelling voor belastingsmeting op een lager met rolelementen.
JP4543643B2 (ja) * 2003-09-12 2010-09-15 日本精工株式会社 転がり軸受ユニットの荷重測定装置
NL1024372C2 (nl) * 2003-09-24 2005-03-29 Skf Ab Werkwijze en sensoropstelling voor belastingmeting op een lager met rollend element gebaseerd op modale vervorming.
ES2335135T3 (es) * 2005-05-30 2010-03-22 Interfleet Technology Ab Metodo y sistema para determinar una pluralidad de componentes de carga sobre rueda.
BE1017141A3 (nl) * 2006-05-17 2008-03-04 Punch Powertrain N V Continu variabele transmissie voor motorvoertuigen.
WO2009056159A1 (en) 2007-11-02 2009-05-07 Aktiebolaget Skf Combination of bearing component and sensor
US9120280B2 (en) * 2009-04-30 2015-09-01 Michelin Recherche Et Technique S.A. Tire uniformity through identification of process signatures
CN102472689B (zh) * 2009-06-30 2015-01-14 米其林研究和技术股份有限公司 在沿轮胎胎圈的不同轨道位置通过烧蚀进行的均匀性校正
CN102498371B (zh) 2009-07-17 2015-07-01 Skf私人有限公司 应变传感器
WO2011128047A1 (en) 2010-04-12 2011-10-20 Aktiebolaget Skf Load on object derived from natural-mode series modelling
CN103038618B (zh) 2010-06-07 2015-02-25 Skf公司 感应机器中轴承上负荷的设备、负荷决定系统及方法
CN102022426B (zh) * 2010-11-26 2012-12-05 天津大学 改善滚动轴承运转性能的方法
DE112012004163B4 (de) * 2011-10-06 2019-06-27 Kabushiki Kaisha Kobe Seiko Sho Reifenuniformitätstestvorrichtung und Reifenuniformitätstestverfahren
EP2838714B1 (de) * 2012-04-19 2018-12-26 Compagnie Générale des Etablissements Michelin Gleichförmigkeitserhöhung durch identifikation diskreter effekte
BR112014026461A2 (pt) * 2012-04-24 2017-06-27 Skf Ab módulo para determinar uma característica de operação de um mancal
CN102705451B (zh) * 2012-05-21 2015-05-27 深圳职业技术学院 消减谐波齿轮传动偏心误差的方法及自动调心波发生器
US8687206B2 (en) * 2012-05-29 2014-04-01 United Technologies Corporation Optical detection of airfoil axial position with NSMS
CN102788695B (zh) * 2012-07-18 2014-12-10 南京航空航天大学 一种滚动轴承磨损的识别方法
WO2014051613A1 (en) * 2012-09-28 2014-04-03 Compagnie Generale Des Etablissements Michelin Improvement of the uniformity of a tire using tooling signature analysis
US20210182298A1 (en) * 2012-12-06 2021-06-17 Autonomous Id Canada Inc. Gait-based biometric data analysis system
JP6195768B2 (ja) * 2013-09-19 2017-09-13 Ntn株式会社 センサ付車輪用軸受のキャリブレーション方法
DE112013007631T5 (de) * 2013-11-25 2016-09-22 Aktiebolaget Skf Vorrichtung und Verfahren der Lagerüberwachung
US10508960B2 (en) * 2013-12-20 2019-12-17 Aktiebolaget Skf Load determining system for a rolling element bearing
GB2531258A (en) * 2014-10-13 2016-04-20 Skf Ab Method and data processing device for detecting a load distribution in a roller bearing
DE112015005624T5 (de) 2014-12-16 2017-09-07 Aktiebolaget Skf Lastbestimmungssystem für ein Wälzlager
GB2542422A (en) 2015-09-21 2017-03-22 Skf Ab Bearing including at least two strain sensor probes
JP2017072438A (ja) * 2015-10-06 2017-04-13 Ntn株式会社 非破壊検査装置および軸受の製造方法
US9677972B2 (en) * 2015-10-26 2017-06-13 Commercial Time Sharing Inc. System and method for characterizing tire uniformity machines
DE102016010549A1 (de) * 2016-08-31 2018-03-01 Sensodrive Gmbh Drehmomentsensor mit Nebenschlussspeiche
DE102016010551B3 (de) * 2016-08-31 2018-02-08 Sensodrive Gmbh Drehmomentsensor mit radialelastischer Momentübertragung
CN107239590A (zh) * 2017-04-12 2017-10-10 上海大学 一种利用有限元建模计算柔性轴承疲劳寿命的方法
WO2018208176A1 (en) * 2017-05-12 2018-11-15 Epps Iain A bearing monitoring method and system
CN107345858B (zh) * 2017-07-25 2018-08-14 安徽大学 一种列车轴承轨边信号冲击成分快速提取方法
US20210033528A1 (en) * 2018-02-19 2021-02-04 Shimadzu Corporation Gas absorption spectroscopic device
CN109342061B (zh) * 2018-09-03 2020-12-15 上海大学 内力耗散力矩平衡式谐波传动双柔性构件性能试验装置
CN109799090B (zh) * 2019-01-08 2020-09-18 长安大学 采用频带3分区的经验小波变换的轴承特征频率提取方法
CN109708891B (zh) * 2019-01-30 2021-07-20 华南理工大学 一种柔性椭圆轴承滚道故障诊断方法
CN109883702B (zh) * 2019-03-05 2020-04-28 西安交通大学 一种基于时频域统计特征的电机轴承故障诊断方法

Also Published As

Publication number Publication date
CN112284584A (zh) 2021-01-29
US20210010883A1 (en) 2021-01-14
US11428590B2 (en) 2022-08-30

Similar Documents

Publication Publication Date Title
DE19937203A1 (de) Überwachung der Lebensdauer und der Belastung von Lagern
DE602004005227T2 (de) Verfahren und sensoranordnung zur lastmessung an einem walzelementlager auf der basis von modelldeformation
EP1944582A1 (de) Verfahren zur bestimmuing einer einflussgrosse auf die exzentrizitat in einem wineklmesser
DE10058424B4 (de) Vorrichtung und Verfahren zur Messung von Kräften und Momenten, die auf in aerodynamischen Windkanälen getestete Modelle wirken
DE10001813C2 (de) Meßsystem zur Ermittlung von Luftdaten eines Luftfahrzeuges sowie ein Verfahren zur Bestimmung der Luftdaten
EP2278262B1 (de) Verfahren zur Ermittlung der Form eines Kurbelzapfens einer Kurbelwelle
DE102018200047A1 (de) Sensorrolle
DE102019210336A1 (de) Verfahren zum Schätzen einer Lagerbelastung unter Verwendung von Spannungsparametern, um eine Kontaktwinkelschwankung zu berücksichtigen
DE112017002650T5 (de) Vorrichtung und Verfahren zur Vorhersage des Ermüdungszustands von Wälzanlagen
EP2439508A1 (de) Messradsatz für Schienenfahrzeuge
DE112020001460T5 (de) Verfahren zum Beziehen eines Kontaktwinkels eines Schrägkugellagers und Verfahren zum Herstellen einer Radlagervorrichtung
DE102017213806A1 (de) Kalibration von Fahrzeugsensoren
DE102011080282A1 (de) Verfahren und Messvorrichtung zur Untersuchung eines magnetischen Werkstücks
DE2326046C2 (de) Verfahren und Vorrichtung zur Qualitätsprüfung von Reifen, insbesondere von Kraftfahrzeugreifen
EP0457200B1 (de) Verfahren zur Messung von Winkeln und Winkelkennlinien mit Hilfe eines Kreisels
DE102018124685A1 (de) Drehmomentmessvorrichtung für ein Spannungswellengetriebe
EP1493006B1 (de) Wälzlager mit Sensoren
DE102016220195A1 (de) Verfahren und Messanordnung zur Erfassung von Schlupf bei Wälzlagern
DE112021002133T5 (de) Anzeichen-Detektionsvorrichtung und Anzeichen-Detektionsverfahren
DE102022116167A1 (de) Sensorvorrichtung
DE102017120756A1 (de) Sensorlagereinheit mit Wälzlager und Messring
DE102021211493A1 (de) Vorrichtung und Verfahren zum Abschätzen einer Last in einem Lager
EP1219921B1 (de) Verfahren zur Bestimmung der Abweichung des Pixelortes der Pixel mindestens einer Bildaufnahmematrix von der Sollposition
DE202014104737U1 (de) Wälzlageranordnung
DE102021124940A1 (de) Verfahren zur Ermittlung der bruchmechanischen Beanspruchung mindestens eines Risses bei Belastung eines Bauteils, Computerprogrammprodukt und Vorrichtung zur Durchführung des Verfahrens, sowie deren Verwendung