Software: OptiY-Workflow - Einfache Toleranzkette

Autoren: Dr.-Ing. PD Thomas Nagel, Dr.-Ing. Alfred Kamusella

Dieser Beitrag wird gerade erstellt und ist noch unvollständig!

Aufgabenstellung

Am Beispiel einer einfachen Maßkette soll gezeigt werden, welche neuen Möglichkeiten die probabilistische Simulation im Vergleich zu klassischen, analytischen Methoden eröffnet:

Das Schlussmaß M₀ ergibt sich aus der Gesamt-Abmessung M₁ abzüglich der Teilmaße M₂ und M₃.

Vorgegeben:

• als Ausgangslösung stehen folgende Maße in der Maßkette:
  • M₁= 11,8-0,2
  • M₂= 1,3-0,1
  • M₃= 1,5±0,05
• diese Ausgangslösung besitzt also die Toleranzmittenmaße C_i mit den Toleranzen T_i:
  • C₁=11,70 mit T₁=0,2
  • C₂=  1,25 mit T₂=0,1
  • C₃=  1,50 mit T₃=0,1

Gesucht:

• das Nennmaß des Schlussmaßes N₀
• Toleranzmittenabmaß des Schlussmaßes E_c0
• Toleranz T₀ des Schlussmaßes
• eine Aussage, ob mit der ermittelten Toleranz T₀ die geforderte Toleranz T_{0_Max} des Schlussmaßes eingehalten wird
• günstigere Toleranzen für die Maßglieder der Toleranzkette, um die zulässige Toleranz des Schlussmaßes eventuell voll auszuschöpfen oder die Fertigung kostengünstiger zu gestalten.

Workflow-Modell der Maßkette

Nach dem Start von OptiY erscheint ein leerer Workflow, den wir nun nutzen, um darauf das Modell der Maßkette zu beschreiben. Dafür einleitend einige Erläuterungen zur Abbildung der toleranzbehafteten Maße auf die Entwurfsparameter (Nennwerte und Streuungen) von OptiY:

• Der zu einer OptiY-Streuung zugeordnete OptiY-Nennwert entspricht dem Toleranz-Mittenmaß C in der Begriffswelt der Maße. Man muss also beachten, dass ein OptiY-Nennwert nicht einem Nennmaß N_i innerhalb der Maßkette entspricht!
• Die OptiY-Streuung Toleranz entspricht dem Wert der Maß-Toleranz T um das Toleranz-Mittenmaß.
• Der aktuelle Wert eines Maßes M_i ist also immer die Summe aus dem Toleranzmittenmaß C_i und dem Istwert der Abweichung innerhalb der Streuung. Das wollen wir für jedes unserer drei Maße im Workflow beschreiben. Die Bedienung des Workflow-Editor wird in diesem Beispiel vorausgesetzt:

Im Beispiel werden alle Workflow-Elemente mit ausführlichen Bezeichnern versehen. Für den erfahrenen Nutzer genügen später auch die Kürzel Ci, Ti und Mi. Die einzelnen Größen sind immer ausführlich zu kommentieren und mit der richtigen Maßeinheit zu versehen.

• Toleranzmittenmaße (als Entwurfsparameter - Nennwerte):
  
  • Für die Toleranz-Analyse (Probabilistische Simulation einer Stichprobe) werden die Toleranzmittenmaße als Konstante angenommen.
  • Die aktuellen Istwerte des zugehörigen Maßes streut dann um diesen konstanten Wert.
• Toleranzen (als Entwurfsparameter - Streuungen):
  
  • Für die Toleranz-Analyse (Probabilistische Simulation einer Stichprobe) werden die Toleranzen als vorgegebene Konstante angenommen.
  • In Vorbereitung auf die Nachbildung der Maximum-Minimum-Methode wurde eine Gleichverteilung für den aktuellen Maßwert angenommen.
• unabhängige Maße (als Transfervariablen):
  
  • Transfer-Variablen bieten die Möglichkeit, ihren Wert über eine Formel zu berechnen.
  • Der aktuelle Ist-Wert eines jeden unabhängigen Maßes Mi berechnet man aus der Summe von Toleranzmittenmaß Ci und aktueller Maßabweichung im Rahmen des Toleranzbereiches Ti:
    
  • Den erforderlichen Formel-Ausdruck bearbeitet man im OptiY-Workflow mit dem Rechner:
    
• Schlussmaß (als Restriktion):
  Datei:Software OptiY-Workflow - Einfache Toleranzkette eigenschaft schlussmass.gif
  • Unter- und Obergrenze der Restriktionsgröße beschreiben den zulässigen Bereich für dieses Schlussmaß. Solange dafür noch keine Werte bekannt sind, sollte man die Grenzen so weit setzen, dass auftretende Werte des Schlussmaßes in jedem Fall zulässig sind (im Beispiel 0 bis 1 Meter).
  • Der Wert des Gewichtsfaktors hat vorläufig keine Bedeutung und sollte den Standardwert=1 behalten
  • Das Schlussmaß M0 wird über einen Formel-Ausdruck aus allen unabhängigen Maßen Mi berechnet:
    
  • Zum Bearbeiten dieses Ausdrucks benutzt man wieder den Rechner.

Maximum-Minimum-Methode

Diese Methode berücksichtigt bei der Berechnung der Schlussmaß-Toleranz die Worst-Case-Fälle und gewährleistet damit eine 100%-ige Einhaltung des Schlussmaßes.

Die klassische Maximum-Minimum-Methode berechnet folgende Größen des Schlussmaßes:

• das Nennmaß N₀ aus den unabhängigen Nennmaßen N_i
  • N₀ = -1 • (1,5 + 1,3 - 11,8) = 9,0
• das Toleranzmittenabmaß E_c0 aus den unabhängigen Toleranzmittenabmaßen E_ci
  • E_c0 = -1 • (0 - 0,05 + 0,1) = -0,05
• die Toleranz T₀ ergibt sich infolge der linearen Kette aus der Summe der Einzeltoleranzen T_i
  • T₀=0,2+0,1+0.1=0,4
• Damit ergibt sich das Schlussmaß M₀=8,95±0,2

Zum gleichen Ergebnis gelangt man mit der probabilistischen Simulation der Maßkette, wenn man jedes unabhängige Maß als gleichverteilt betrachtet. Die Elemente der Maßkette wurden beim Aufbau des OptiY-Workflows bereits richtig konfiguriert.

Die probabilistische Simulation der Maßkette ín Anlehnung an die Maximum-Minimum-Methode soll nun das erste Experiment sein, welches wir mit dem Workflow-Modell durchführen:

• Wir werden noch mehrere Experimente in diesem OptiY-Versuchsstand durchführen. Deshalb sollte man über das Kontextmenü (Rechter Mausklick) das Experiment Umbenennen in Maximum-Minimum.
• Als "Optimierungsverfahren" wählen wir die "Simulation".
  • Da Streuungen im Workflow-Modell existieren, handelt es sich hier um eine probabilistische Simulation unter Berücksichtigung der Toleranzen.
  • Praktisch entspricht dies der Simulation einer Stichprobe.
    
    

• Nach welchem Verfahren die probabilistische Simulation durchgeführt wird, konfiguriert man in der statistischen Versuchsplanung:
  • Für den einfachen linearen Zusammenhang ist eine Moment-Methode das genaueste und effizienteste Verfahren.
  • Es soll vorläufig als Ersatzfunktion M0=f(T1,T2,T3) eine Taylorreihe zweiter Ordnung genutzt werden.
  • Man sollte immer mit der Berücksichtigung von Interaktionen zwischen den Toleranzgrößen beginnen. Nur so kann man feststellen, ob wirklich Wechselwirkungen zwischen den Toleranzgrößen existieren.
• Mittels Analyse - Verteilungsdichte öffnen wir ein Grafikfenster und ziehen per Drag&Drop die Toleranzen und das Schlussmaß aus dem OptiY-Explorer in dieses Fenster.
• Nach Projekt - Start erfolgt die Simulation der Stichprobe und die Verteilungsdichten der gewählten Größen werden dargestellt. Leider kann man anhand der dargestellten Grenzwerte von Schlussmaß_M0 nur näherungsweise die tatsächlichen Grenzen abschätzen. OptiY schneidet Werte der Dichtefunktion ab, die kleiner als 1/100 ihres Maximalwertes sind:

• Mittels Analyse - Verteilungstabelle erhält man Stützstellen der Verteilungsfunktion aller Bewertungsgrößen für den Wertebereich von 0 bis 1 aufgelistet. In der ersten und letzten Zeile dieser Tabelle findet man die gesuchten Grenzen für das Schlussmaß:

• Bis auf geringe Rundungsfehler haben entspricht dies den Werten der klassischen Maximum-Minimum-Methode.

Das Ergebnis der Maximum-Minimum-Methode bedeutet:

Bei den vorgegebenen Fertigungsgenauigkeiten liegt mit einer Wahrscheinlichkeit von 100% das Schlussmaß innerhalb der Grenzen 8,95±0,2.

Wahrscheinlichkeitsbasierte Methode

Der Vorteil der probabilistischen Simulation ist die Anschaulichkeit der grafischen Ergebnisdarstellung:

• So gibt es bei den meisten Betrachtern sicher ein gewisses Aha-Erlebnis, dass trotz gleichverteilter Einzelmaße das Schlussmaß im Prinzip einer Normalverteilung entspricht.
• Es zeigt sich z.B. eine Häufung des Schlussmaßes in der Nähe des Toleranzmittenmaßes.
• Nur ein geringer Prozentsatz einer Stichprobe liegt mit dem Schlussmaß in der Nähe der Grenzwerte.

An dieser Erkenntnis setzt die wahrscheinlichkeitsbasierte Methode der Toleranzrechnung an:

• Alle unabhängigen Maße streuen um ihr Toleranzmittenmaß mit einer Verteilungsfunktion, welche von unterschiedlichsten Einflussgrößen abhängt.
• Im Normalfall oder bei mangelhafter Kenntnis der Verteilungsfunktionen kann man von einer Normalverteilung ausgehen.
• Die Verteilungsfunktion des Schlussmaßes ist abhängig von der Maßkette und den Verteilungsfunktionen der Einzelmaße.

Für lineare Ketten kann man die Toleranz T0 des Schlussmaßes auch noch analytisch berechnen:

• Der Streuungskoeffizient k_i der Einzelmaße ist abhängig von der Verteilungsfunktion (0,333 für Normalverteilung).
• Der Risikofaktor t für die Toleranz des Schlussmaßes ist abhängig davon, wieviel Prozent der Stichprobe im berechneten Toleranzbereich liegen sollen. Dabei wird die Studentverteilung zugrunde gelegt:
• Bei vorgebener Ausfallwahrscheinlichkeit (im Beispiel 0,3%, was T₀=6•σ entspricht), kann man über einen Web-Rechner der Uni Saarland nach Wahl des Ausfallbereiches, der Angabe großen Wertes df=1000000 und der Wahrscheinlichkeit des Ausfalls von z.B. 0.003 den exakten Wert von t=2.968 erhalten
  

• Die Toleranzbreite T₀ berechnet man mit diesen Werten:

[math]\displaystyle{ T_0=2{,}968 \cdot \sqrt{(0{,}333 \cdot 0{,}1 )^2+(0{,}333 \cdot 0{,}1)^2+(0{,}333 \cdot 0{,}2)^2}=0{,}2421 }[/math]

• Damit ergibt sich das Schlussmaß zu M₀= 8,95±0,12105 (normalverteilt).

Das bedeutet:

Bei gleichen Fertigungsgenauigkeiten wird eine reduzierte Schlussmaß-Toleranz erreicht! Es wird mit 99,7% Wahrscheinlichkeit im Beispiel ein Schlussmaß innerhalb der Grenzen 8,95 ± 0,12602 eingehalten.

Das werden wir nun in einem weiteren Experiment mit unserem Workflow-Modell anschaulich nachvollziehen:

• Dazu erstellen wir zuerst ein neues Experiment Wahrscheinlichkeit durch Duplizieren des Maximum-Minimum-Experiments (Kontextmenü der rechten Maustaste).
• Diese Kopie des Experiments erhält durch anschließendes Umbenennen den gewünschten Namen.
• Beim Kopieren werden der Workflow und alle Experiment-Einstellungen der Kopier-Vorlage übernommen. Die Darstellung der Ergebnisse muss man jedoch erneut mittels des Analyse-Menüs konfigurieren.
• Für alle Einzel-Toleranzen wählen wir anstatt der Gleichverteilung die Normalverteilung.
• Die anschließende Simulation ergibt eine Standard-Abweichung σ₀=0,0408248 für das Schlussmaß. Damit ergibt sich eine Toleranz T₀=6•σ=0,245 und das Schlussmaß zu M₀= 8,95±0,1225. Die mit dem Workflow-Modell ermittelte Toleranz des Schlussmaßes weicht im Beispiel um weniger als 1% (ca. 3 µm) vom Wert des analytischen Ansatzes ab.