Der Sporn der Läuferin sitzt an einer bestimmten Stelle: an der Unterseite des Fersenbeins, im Sehnenansatz, dort wo die Fußsohle vorn an der Ferse beginnt. Auf dem Schema aus Kapitel 1 ist er ein kleiner knöcherner Fortsatz, genau da. Zwei Theorien streiten seit Jahrzehnten, warum er dort wächst — die eine sagt Zug, die andere Druck. Dieses Kapitel entscheidet den Streit nicht. Es macht ihn rechenbar: Wir setzen die geschlossene Umbauschleife aus Kapitel 17 einmal auf einen reinen Zug an der Enthese und einmal auf reinen Bodendruck an, beide unter demselben Maßstab, und sehen nach, wo jede Hypothese ihren Knochen baut und welches Muster sie legt. Am Ende steht eine Frage, die das Buch stellt und nicht mehr beantwortet: was eine Einlage an diesem Muster ändern würde.
Lernziele
Nach diesem Kapitel kannst du
die Zughypothese und die Druckhypothese des Fersensporns benennen und sagen, warum die Literatur streitet (Kumai und Benjamin 2002 deuten histologische Befunde als Druckhinweis, Zwirner und Mitarbeitende 2021 wieder als Zughinweis — der Streit ist offen),
erklären, wie ein Modell zwei Hypothesen vergleichbar macht: beide unter einem gemeinsamen Referenzreiz rechnen, Aufbau-Ort und Aufbau-Muster nebeneinanderlegen — und dabei sagen, welche Lesart jeder Hypothese das Modell überhaupt auf die Probe stellt,
den Befund des Modells lesen: an welchem Ort jede der beiden Hypothesen aufbaut, wie man prüft, ob dieser Ort am Netz hängt, und warum der Zugaufbau am Sporn-Ort kein Beleg für die Zughypothese ist,
die zwei Modellgrenzen deuten (Sättigung mit Fenster-Kriterium; Lagernähe an der Talus-Kante, dicht am Druck-Cluster) und benennen, welcher Baustein einem Modell fehlt, das die Wirkung eines Eingriffs — einer Einlage — ausrechnen soll.
Das Rechenwerk lädt das offline vorberechnete Datenpaket der Hypothesen-Läufe und legt die Zeichenhilfen bereit. Es rechnet selbst nichts und deutet nichts — aufklappen lohnt nur, wer einer Zahl oder einem Namen nachgehen will.
Das Rechenwerk dieses Kapitels — zum Nachlesen aufklappen
# --- Nur Zeichen- und Ladehilfen. Die EINE Rechenquelle der Hypothesen-Läufe ist# programme/kap18/hypothesen.py; jedes gezeichnete Feld stammt aus dem OFFLINE# vorberechneten Datenpaket hypothesen_daten.json (Kap.-12/16-Muster). Diese Zelle# konfiguriert nur und deutet nichts.import jsonimport pathlibimport sysimport numpy as npimport matplotlib.pyplot as pltimport matplotlib.patheffects as pfximport matplotlib.tri as mtrifrom matplotlib.colors import LogNormfrom matplotlib.patches import Circlesys.path.insert(0, str(pathlib.Path("programme/gemeinsam").resolve()))sys.path.insert(0, str(pathlib.Path("programme/kap14").resolve()))sys.path.insert(0, str(pathlib.Path("programme/kap17").resolve()))sys.path.insert(0, str(pathlib.Path("programme/kap18").resolve()))import hypothesen as H18komma = H18.kommaPAKET = json.loads(pathlib.Path("programme/kap18/hypothesen_daten.json").read_text())NETZE = {s: json.loads(pathlib.Path("netze/netz_%s.json"% s).read_text())for s in ("grob", "mittel", "fein")}KONTUR = np.array(json.loads( pathlib.Path("../geometrie/kalkaneus_kontur.json").read_text())["punkte"]) *1000.0# Die beiden Netzstufen der Vergleichsrechnung (Geometrie + Cluster); der# Netzstabilitäts-Regler zeichnet sie, der Fließtext zitiert ihre Kennzahlen.REFNETZ = json.loads(pathlib.Path("netze/referenznetz.json").read_text())["stufen"]def wandert(lastfall):"""Wie weit das Cluster-Zentrum von der feinen zur feineren Stufe rückt (mm).""" a = REFNETZ["fein"][lastfall]["zentrum_mm"] b = REFNETZ["feiner"][lastfall]["zentrum_mm"]return ((a[0] - b[0]) **2+ (a[1] - b[1]) **2) **0.5# Feste Dichtezunahme-Skala: inferno, 0–0,80 g/cm³. Das ist NICHT die# Standard-Dichteskala (magma 0,2–1,0) — die Dichtezunahme Δρ ist eine andere Größe,# deshalb eine eigene, an jeder Abbildung benannte Skala.ZU_MIN, ZU_MAX =0.0, PAKET["parameter"]["zunahme_vmax"]SPORN = PAKET["parameter"]["sporn_ort_mm"]MISCH = PAKET["regler_mischung"]["werte"] # die α-Stufen des Mischungs-Reglersdef ort(punkt):"""Ein Ort als „(x / y)“ in Millimetern mit deutschen Kommas."""return"(%s / %s)"% (komma(punkt[0], 1), komma(punkt[1], 1))def abstand_zum_ansatz(punkt):"""Abstand eines Ortes zum Faszienansatz (= realer Sporn-Ort) in Millimetern."""return ((punkt[0] - SPORN[0]) **2+ (punkt[1] - SPORN[1]) **2) **0.5def zunahme_bild(ax, netz, feld, titel, zentrum_mm=None, netz_kanten=False, schrift=7):"""Ein Dichtezunahme-Feld über einem Buchnetz, feste inferno-Skala 0–0,80 g/cm³, mit Ort, Orientierung und 10-mm-Maßstab. ``zentrum_mm`` markiert das Cluster-Zentrum (Kreis). Der Titel steht als echte Achsen-Überschrift (nicht als Text im Feld, damit die gemeinsame Colorbar ihn nicht schneidet); der 10-mm-Maßstab liegt auf EIGENER Höhe UNTER der Orientierungszeile (Kapitel-14-Entzerrung). ``netz_kanten`` legt die Dreieckskanten über das Feld (Kapitel-17-Konvention: weiß, dünn, halbdurchsichtig). Nur für die Netzstabilitäts-Reihe: Dort ist das Netz selbst der Gegenstand, und ohne Kanten sähen alle drei Stufen gleich aus.""" kn = np.asarray(netz["knoten"]) *1000.0 tri = mtri.Triangulation(kn[:, 0], kn[:, 1], netz["dreiecke"]) tp = ax.tripcolor(tri, facecolors=np.asarray(feld), cmap="inferno", vmin=ZU_MIN, vmax=ZU_MAX)if netz_kanten: ax.triplot(tri, color="white", lw=0.18, alpha=0.35)if zentrum_mm isnotNone: ax.plot(zentrum_mm[0], zentrum_mm[1], "o", mfc="none", mec="#38bdf8", mew=2.4, ms=15) ax.set_title(titel, fontsize=schrift +1, fontweight="bold", pad=3)# Orientierungszeile kompakt (hinten/vorn/oben) — kurze Labels, damit sie sich in# gestapelten schmalen Panels nicht überlappen; Ferse/Zehen stehen in der Caption. ax.text(0, -5.0, "hinten", fontsize=schrift, color="#666") ax.text(66, -5.0, "vorn", fontsize=schrift, color="#666", ha="right") ax.annotate("oben", xy=(-6, 24), xytext=(-6, 12), fontsize=schrift, color="#444", ha="center", arrowprops=dict(arrowstyle="->", color="#444", lw=1.1))# 10-mm-Maßstab auf EIGENER Höhe, klar unter der Orientierungszeile. ax.plot([28, 38], [-12.5, -12.5], color="#333", lw=2.0, solid_capstyle="butt") ax.text(33, -15.0, "10 mm", ha="center", fontsize=schrift, color="#333") ax.set_xlim(-10, 70); ax.set_ylim(-17, 44); ax.set_aspect("equal"); ax.axis("off")return tpdef gemeinsame_colorbar(fig, tp, achsen, laenge=52, schrift=7):# Colorbar VOLLSTÄNDIG rechts außerhalb der Panelspalte (constrained_layout der# Figur reserviert den Platz; die Panel-Titel sind echte set_title und werden# nicht geschnitten). ``laenge`` ist das Seitenverhältnis des Farbbalkens: Ohne# die Vorgabe blieb er bei der Standardlänge und stand nur neben dem MITTLEREN# Panel, als gehörte er allein zu diesem. cb = fig.colorbar(tp, ax=achsen, fraction=0.028, pad=0.03, location="right", aspect=laenge) cb.set_label("Dichtezunahme Δρ (g/cm³), inferno 0–0,80", fontsize=schrift) cb.ax.tick_params(labelsize=schrift)
18.1 Der Streit: zieht die Faszie, oder drückt der Boden?
Die naheliegende Erklärung des Fersensporns steht in vielen Lehrbüchern: Die Plantarfaszie spannt sich von den Zehen zur Ferse; wo sie am Fersenbein ansetzt, zieht sie bei jedem Schritt. Der Knochen, so die Idee, antwortet auf diesen Zug und wächst ihm entgegen — heraus kommt der Sporn. Das ist die Zughypothese, und sie klingt zwingend. Nur ist sie in der Forschung umstritten.
WarnungEine publizierte Gegenposition zur Zugdeutung
Verbreitete Vermutung:„Die Plantarfaszie zieht am Fersenbein, der Knochen wächst dem Zug entgegen — daher der Sporn.“ Das ist die Zughypothese, und sie liegt nahe: Der Sporn sitzt am Faszienansatz, und Zug ist die Kraft, die man dort am eigenen Fuß spürt.
Eine Studie, die anders deutet — Beobachtung und Deutung getrennt. Kumai und Benjamin (Journal of Rheumatology 2002, Band 29, Heft 9, Seiten 1957–1964) untersuchten radiologisch und histologisch die Fersen von 17 älteren Leichen — also verstorbenen Menschen, keine Belastungen am lebenden Fuß. Beobachtet haben sie unter anderem, dass Sporne auch an der tiefen Seite der Faszie auftreten, dass das Gewebe degenerative Veränderungen zeigt und dass größere Sporne innen senkrecht verlaufende Knochenbälkchen (vertikale Trabekel) tragen. Gedeutet haben die Autoren das als Hinweis auf wiederholte vertikale Druck- statt Zugbelastung. Wichtig: Die vertikalen Trabekel sind eine Struktur im Sporn — sie sagen nichts über die Ausrichtung des ganzen Sporns, und die Studie misst keine Belastung und weist keinen Ursache-Wirkungs-Zusammenhang nach.
Die Gegenposition steht ebenfalls in der Literatur. Zwirner und Mitarbeitende (Scientific Reports 2021) kommen morphologisch ausdrücklich wieder zur Zugdeutung und betonen, dass die Trabekelorientierung allein den Streit nicht entscheidet. Der Streit ist also offen — beide Deutungen sind publiziert.
Was dieses Kapitel daraus macht: Es entscheidet den Streit nicht. Es rechnet beide Hypothesen mit derselben Schleife und legt Ort, Muster und Eingriffs-Antwort nebeneinander. Der reifere Nutzen eines Modells liegt nicht darin, die Wahrheit zu behaupten, sondern zu zeigen, welche Frage die Hypothesen überhaupt trennt — und welche nicht.
18.2 Wie man zwei Hypothesen vergleichbar macht
Zwei Hypothesen sind nur dann vergleichbar, wenn man sie an denselben Maßstab hält. Der Reiz, auf den der Mechanostat aus Kapitel 16 antwortet, ist die Abweichung vom Referenzreiz \(S_0\) — dem Sollwert, bei dem weder auf- noch abgebaut wird. Gäbe man dem Zug- und dem Druckszenario je einen eigenen\(S_0\) (jeweils das Mittel seines eigenen Reizfelds), dann rechnete man beide auf ihr eigenes Gleichgewicht herunter und hätte die Lastunterschiede gerade wegnormalisiert — der Reiz wächst quadratisch mit der Last, und genau diese Größenunterschiede sind das Interessante. Deshalb gilt hier, wie in der Vergleichsrechnung, ein gemeinsamer Referenzreiz für beide: das flächengewichtete Mittel der beiden Start-Reizfelder. Auf dem feinen Buchnetz ist das \(S_0 \approx 4\,8\) J/kg — ein Modellwert dieser gröberen Rechnung, nahe dem Referenzreiz der bildbasierten Kontur (4,84 J/kg).
Das Modell dieses Kapitels ist der Aufbau der numerisch geprüften Vergleichsrechnung, auf die Buchnetze übertragen: eine feste kortikale Schale am Rand (17 GPa, sie wird nicht umgebaut), ein spongiöser Kern und — das ist neu gegenüber Kapitel 17 — eine weiche Enthese am vorderen unteren Abhang des Fersenbeins, dem realen Sporn-Ort, mit niedriger Startdichte (0,20 g/cm³). Dort, im faserknorpeligen Ansatzgewebe, entsteht der Enthesophyt beim realen Menschen durch Verknöcherung; im Modell steht dafür allein die Dichtezunahme dieser weichen Region. Das ist damit der vierte Knochen, den dieses Buch rechnet — Tabelle 13.3 in Kapitel 13 stellt alle vier nebeneinander: die homogen kortikalen Buchnetze der Kapitel 10 bis 15, das gleichmäßige Dichtefeld der Kapitel 16 und 17, den Aufbau dieses Kapitels (Schale, Kern, weiche Enthese) und die Vergleichsrechnung auf der bildbasierten Kontur. Wer Zahlen aus zwei Zeilen dieser Tabelle nebeneinanderlegt, vergleicht zwei verschiedene Knochen. Der Aufbau-Ort schließlich ist keine Einzelzelle, sondern eine definierte Größe: das flächengewichtete Zentrum des größten zusammenhängenden Clusters der obersten zehn Prozent der Dichtezunahme (die Definition rechnen wir gleich von Hand nach).
18.2.1 Welchen „Druck“ dieses Kapitel prüft
Ein Wort noch zum zweiten Lastfall, bevor der Regler etwas zeigt, das leicht mehr zu beweisen scheint, als es beweist. „Druck“ heißt in diesem Kapitel: die Bodenlast der Standphase, verteilt über die Auftrittsfläche, bei fest gelagerter Talus-Gelenkfläche — genau die Last, die in Kapitel 12 eingeführt wurde. Geprüft wird damit die Lesart: Baut der Lastpfad der Standphase den Knochen dort auf, wo der Sporn sitzt?
Kumai und Benjamin meinen mit Druck etwas Engeres: eine wiederholte örtliche Zusammendrückung des Ansatzgewebes selbst, also eine Kraft, die unmittelbar an der Spornregion senkrecht auf die Enthese drückt. Diese zweite Lesart zu prüfen hieße, eine eigene Last anzusetzen — eine Normallast auf die Enthese-Facette oder auf die Auftrittszone unmittelbar darunter — und sie gegen den Zug zu stellen. Das rechnet dieses Buch nicht: Woher diese Last käme und wie groß sie wäre, hinge daran, wie der Fuß den Boden berührt, und die Berührung selbst gibt jedes Modell dieses Buches vor, statt sie auszurechnen (Kapitel 1). Wo im Folgenden „der Druck baut hier“ steht, ist deshalb immer der Bodendruck der Standphase gemeint, nicht die lokale Kompression der Literatur. Das ist kein Nebensatz, sondern die Grenze des Vergleichs, den dieses Kapitel anstellt.
Bevor wir umschalten, eine Vorhersage.
WichtigVorhersage-Punkt
Bevor du die Animation startest: Erst baut das Modell unter reinem Zug an der Enthese auf, dann schaltest du auf reinen Druck von unten um. Wandert der Aufbau-Ort weit, bleibt er an derselben Stelle, oder verschwindet er? Und woran willst du die beiden Hypothesen dann noch unterscheiden, wenn nicht am Ort? Lege dich fest, dann spiele ab.
Deutung der Animation. Zwei Hypothesen, zwei ganz verschiedene Bilder. Unter Zug verknöchert das weiche Ansatzgewebe — ein heller Sättigungsfleck wächst vorn unten am Abhang, genau in der Spornregion. Unter Druck entsteht etwas anderes: eine Säule zieht sich den Lastpfad hinauf, von der Ferse zur fest gelagerten Talus-Facette oben, während die Enthese fast unberührt bleibt. Der Aufbau-Ort wandert also nicht ein Stück weit — er sitzt bei den beiden Hypothesen an zwei völlig verschiedenen Zonen. Das ist der erste Teil der Antwort auf die Vorhersagefrage: Der Ort trennt die beiden Hypothesen. Der zweite Teil — was dieser Ortsbefund beweist und was nicht — ist die eigentliche Lehre dieses Kapitels.
Der Regler legt die beiden Endzustände direkt nebeneinander, damit man die Muster und die Aufbau-Orte messen kann.
WichtigVorhersage-Punkt
Bevor du umschaltest: Schätze, wo jedes der beiden Felder auf diesem feinen Buchnetz seinen stärksten Aufbau hat. Liegen die beiden Aufbau-Orte nah beieinander oder weit auseinander? Und — wichtiger — was unterscheidet die beiden Felder außer dem Ort? Lege dich fest, dann schalte um.
Deutung. Auf dem Buchnetz sind die beiden Aufbau-Orte weit voneinander entfernt: Der Zug-Cluster sitzt bei rund (29 / 12,2) mm, der Druck-Cluster bei rund (17,6 / 24,4) mm oben an der Talus-Säule — gut 17 mm auseinander. Beim Zug ist mit dieser Zahl Vorsicht geboten, und der Kasten „Zwei Ortszahlen, zwei Fragen“ weiter unten sagt genau warum: Sie ist die Mitte eines langgezogenen Bandes, nicht dessen hellste Stelle. Der stärkste Zug-Aufbau liegt bei (23,6 / 5,5) mm, also nur 2,17 mm vom realen Sporn-Ort — das ist die Spornregion. Man könnte einwenden, das sei nur die Buchnetz-Demo. Genau hier setzt die Vergleichsrechnung an — und sie bestätigt es.
18.3 Der Befund: der Druck baut woanders
Dass der Zug am Sporn-Ort aufbaut, konnte nicht ausbleiben. Der Zugursprung wurde dorthin gelegt, und eine Zuglast erzeugt an ihrem eigenen Ansatz die höchste Spannung; die weiche Enthese sitzt ebenfalls dort, und ohne sie gäbe es am Sporn-Ort gar keinen Aufbau. Der Sporn-Ort steckt also zweimal in den Annahmen dieser Rechnung, bevor sie beginnt. Der Befund steckt im anderen Fall: Der Bodendruck baut nicht am Sporn-Ort, sondern oben, wo sein Lastpfad an der Talus-Facette endet — in der Vergleichsrechnung 17 Millimeter vom Zugaufbau entfernt, und auf einem feineren Netz ebenso. Unter den Annahmen dieses Modells passt der Bodendruck räumlich nicht zum Sporn. Was das beweist und was nicht, sagt der Kasten weiter unten; und ob dieser Druckfall überhaupt die Druckhypothese der Literatur trifft, ist die Frage, die der Abschnitt „Welchen ‚Druck‘ dieses Kapitel prüft“ offen benannt hat.
Die numerisch geprüfte Vergleichsrechnung ist denselben Test auf der feinen bildbasierten Kontur gefahren, mit derselben Methode: gemeinsamer Referenzreiz, Cluster-Zentrum als Aufbau-Ort, Faszienursprung am realen Sporn-Ort. Ihre Zahlen stehen einmal vollständig im Kasten; spätere Nennungen im Kapitel runden sie.
HinweisDie Zahlen der Vergleichsrechnung
Größe
Wert
was sie sagt
Cluster-Zentrum Zug
(21,8 / 6,1) mm
der Zug baut in der Spornregion
Cluster-Zentrum Druck
(18,3 / 23,0) mm
der Druck baut oben an der Talus-Facette
Ortsabstand der beiden
17,33 mm, auf dem feineren Netz 17,58 mm
die Orte liegen weit auseinander, und das bleibt so
Schwelle
2,06 mm (zwei Elementkanten)
darunter wäre ein Ortsunterschied nicht auflösbar
Vollfeld-Differenz Zug ↔︎ Druck
0,746 (relative \(L_2\))
auch die ganzen Felder sehen verschieden aus
gemeinsames \(S_0\)
4,84 J/kg
der Maßstab, unter dem beide gerechnet wurden
Alle Werte gehören zur Rechnung auf der bildbasierten Kontur — es sind keine Buchnetz-Zahlen.
Das Zug-Cluster liegt in der Spornregion, dicht vor den Lastfacetten, dort wo der echte Sporn sitzt. Das Druck-Cluster liegt oben an der fest gelagerten Talus-Facette. Sie sind 17,33 mm auseinander (auf dem noch feineren Netz 17,58 mm), weit über der Schwelle von 2,06 mm — und das bestätigt, was schon das Buchnetz zeigte.
Auch das Muster des ganzen Feldes trennt die beiden Fälle: Die Vollfeld-Differenz der beiden Zunahmefelder beträgt 0,746 (relative \(L_2\), Zug gegen Druck innerhalb derselben Rechnung). Der Zug sättigt am Abhang; der Druck baut die Lastpfad-Säule hinauf zur Talus-Facette. Ort und Muster ziehen an einem Strang.
HinweisZwei Modellrechnungen, ein Ortsbefund — was er beweist und was nicht
Beide Zahlenfamilien sagen dasselbe. Die Buchnetz-Demo dieses Kapitels (die gröbere Rechnung, wie schon in Kapitel 13 und Kapitel 17) legt die beiden Aufbau-Orte rund 17 mm auseinander: Zug in der Spornregion, Druck oben an der Talus-Säule. Die Vergleichsrechnung (bildbasierte Kontur, eigenes feineres Netz, mit voll aufgelöster weicher Enthese) findet 17,33 mm, mit demselben Bild. Dass zwei verschieden feine Rechnungen (anderes Netz, anderes \(S_0\)) hier übereinstimmen, macht den Ortsbefund robust.
Und doch ist Vorsicht geboten, gleich dreifach. Erstens ist der Zugaufbau am Sporn-Ort kein unabhängiger Beleg der Zughypothese, aus dem oben genannten Grund: Das Zug-Cluster liegt nur 3,5 mm vom Faszienansatz entfernt, den das Modell selbst dorthin gesetzt hat; diese Vorhersage konnte gar nicht scheitern. Tragfähig ist die andere Zahl: Der Druckfall baut 20,5 mm vom Ansatz entfernt auf, oben an der Talus-Säule — dass die beiden Fälle im Ort so weit auseinanderfallen, war nicht gesetzt. Das Modell zeigt damit, dass der Zugfall räumlich zum Sporn passt, nicht dass Zug die Ursache ist. Zweitens sitzt das Druck-Cluster dicht an der fest eingespannten Talus-Gelenkfläche (Modellgrenze zwei) — Lagernähe kann den Aufbau dort überzeichnen. Drittens prüft der Druckfall die Bodenlast der Standphase, nicht die örtliche Kompression des Ansatzgewebes, von der Kumai und Benjamin sprechen. Ihre klinische Gegenposition wird durch dieses Modell also nicht widerlegt; sie bleibt die reale Referenz. Was der Ortsbefund leistet, ist didaktisch: Er macht zwei Lastfälle, über die man sonst nur reden kann, am Ort und am Muster unterscheidbar.
Damit steht der Sporn zum dritten Mal im Buch an einer Stelle, an der man leicht mehr mitnimmt, als dasteht — nach dem ersten Regler in Kapitel 1 und der Einlage in Kapitel 12. Also noch einmal der Kasten: Wer sagt hier eigentlich was?
HinweisDas Modell sagt · die Studienlage sagt · dein Arzt sagt
Das Modell sagt: Unter reinem Zug an der Enthese verdichtet sich die Spornregion, unter reinem Bodendruck die Säule oben an der Talus-Facette, rund 17 Millimeter davon entfernt. Das ist eine Aussage über dieses Modell — über seine Kontur, seine angenommenen Einzellasten, seine weiche Enthese am Sporn-Ort und sein fest eingespanntes Talus-Lager.
Die Studienlage sagt: Ob Zug oder Druck den Sporn entstehen lässt, ist bis heute strittig; Kumai und Benjamin deuten 2002 histologische Befunde als Druckhinweis, Zwirner und Mitarbeitende 2021 morphologische wieder als Zughinweis. Ein Sporn findet sich auch bei Menschen ohne jede Beschwerde, und ob er selbst schmerzt, ist nicht geklärt. Keine dieser Arbeiten misst eine Belastung im lebenden Fuß.
Dein Arzt sagt: Was an einem bestimmten Fuß los ist und was ihm hilft — Pause, Übungen, Einlage —, entscheidet eine Untersuchung, keine Rechnung. Dieses Buch gibt dazu keine Empfehlung und ersetzt keine.
Wie robust der Druck-Ort ist, zeigt eine kleine Probe am gemeinsamen Referenzreiz \(S_0\). Der ist die zentrale Stellschraube eines fairen Vergleichs — gäbe man beiden Hypothesen je einen eigenen, hätte man die Lastunterschiede wegnormalisiert. Gibt man dem Buchnetz-Druck statt seines eigenen \(S_0 \approx 4\,8\) J/kg den das \(S_0\) der Vergleichsrechnung von 4,8 J/kg — sonst nichts —, dann bleibt das Druck-Cluster praktisch an derselben Stelle: von rund (17,6 / 24,4) mm zu rund (17,5 / 24,4) mm — beide Male oben an der Talus-Facette. Im geprüften Bereich (rund 4,1 bis 4,8 J/kg) hängt der neue Druck-Ort also nicht an der Wahl des \(S_0\); er folgt den Lastpfaden zur Talus-Facette, nicht dem gewählten Reizmaßstab. Dass \(S_0\) für beide Hypothesen derselbe sein muss, bleibt trotzdem die Regel des fairen Vergleichs.
Der reale Sporn-Ort auf der Buchnetz-Kontur ist der Enthese-Konturpunkt (24,07 / 3,35) mm — derselbe Faszien-Ansatzpunkt wie in Kapitel 1. Jetzt kommt eine Frage, die harmloser klingt, als sie ist: Wie weit ist der Aufbau von diesem Punkt entfernt? Sie hat zwei Antworten, und sie unterscheiden sich um das Fünffache. Welche man bekommt, hängt davon ab, womit man den Aufbau-Ort misst.
HinweisZwei Ortszahlen, zwei Fragen
Der Zug-Aufbau ist kein runder Fleck. Er ist ein langgezogenes Band, das an der Enthese am hellsten ist und von dort nach vorn-oben ausläuft. Für so eine Form geben zwei vernünftige Ortsmaße zwei verschiedene Antworten:
Der Bandschwerpunkt — das flächengewichtete Zentrum des größten zusammenhängenden Clusters, die Definition der Vergleichsrechnung, die wir gleich von Hand nachrechnen. Er beantwortet: Wo liegt die Mitte der Aufbauzone?
Der Spitzenort — der flächengewichtete Schwerpunkt aller Elemente, deren Zunahme mindestens 95 % des Maximums erreicht. Er beantwortet: Wo baut das Modell am stärksten auf?
Bei einem runden Fleck fallen beide zusammen. Bei einem Band nicht: Der Schwerpunkt liegt in der Mitte, die Spitze an einem Ende. Und weil ein feineres Netz das Band länger und schlanker auflöst, wandert der Bandschwerpunkt mit der Netzfeinheit, der Spitzenort nicht:
Buchnetz
Bandschwerpunkt
Spitzenort
grob
6,2 mm
1,86 mm
mittel
5,23 mm
2,51 mm
fein
10,13 mm
2,17 mm
(Abstand zum realen Sporn-Ort.) Die Frage „wo baut das Modell am stärksten auf?” hat also eine netzstabile Antwort, die Frage „wo liegt die Mitte der Zone?” nicht.
Damit steht fest, welche Zahl wo gebraucht wird, und beide bleiben im Kapitel. Geht es um die Nähe zum Sporn — also darum, ob der Zugaufbau dorthin fällt, wo der echte Sporn sitzt —, zählt der Spitzenort; er steht deshalb überall dort im Vordergrund, wo dieses Kapitel den Zugaufbau mit dem Sporn-Ort vergleicht. Geht es um den Abstand zwischen den beiden Aufbau-Zonen — also darum, ob Zug und Druck überhaupt auseinanderliegen —, zählt der Bandschwerpunkt: Er ist die Größe, mit der die Vergleichsrechnung misst, und beide Zonen werden mit demselben Maß gemessen. Der Mischungsregler und die Netzstabilität weiter unten arbeiten deshalb mit ihm, und sie sagen es an Ort und Stelle noch einmal.
Abbildung 18.1 legt beide Zahlen neben den realen Sporn-Ort.
Code
g = PAKET["haupt"]fig, ax = plt.subplots(figsize=(6.4, 4.6))# Konturumriss.kx =list(KONTUR[:, 0]) + [KONTUR[0, 0]]ky =list(KONTUR[:, 1]) + [KONTUR[0, 1]]ax.plot(kx, ky, color="#8a6d3b", lw=1.6)ax.fill(kx, ky, color="#f0e6d2", alpha=0.5)# Der Bandschwerpunkt auf allen drei Netzstufen — seine Wanderung ist die Aussage.# Die gepunktete Kette zeigt sie; grob und mittel liegen nur 1,0 mm auseinander, ihre# Marken berührten sich bei direkter Beschriftung. Deshalb je Stufe eine eigene# Markenform und die Namen in der Legende (oben links, im leeren Bildteil).kette = [st["zug_zentrum_mm"] for st in PAKET["regler_netz"]["stufen"]]ax.plot([p[0] for p in kette], [p[1] for p in kette], ":", color="#1f6feb", lw=1.2)form = {"grob": ("o", 9), "mittel": ("^", 9), "fein": ("D", 8)}for st in PAKET["regler_netz"]["stufen"]: zs = st["zug_zentrum_mm"] marke, groesse = form[st["stufe"]] ax.plot(zs[0], zs[1], marke, mfc="none", mec="#1f6feb", mew=1.8, ms=groesse, label="Bandschwerpunkt, Netz %s"% st["stufe"])# Der Spitzenort — die netzstabile Zahl. Er liegt nur 2,17 mm vom Sporn-Ort, also# müssen beide Marken klein bleiben, damit sie sich nicht überdecken.sz = g["zug"]["spitze_mm"]ax.plot(sz[0], sz[1], "s", mfc="none", mec="#d62728", mew=2.2, ms=9, label="Spitzenort (stärkster Aufbau)")ax.plot(SPORN[0], SPORN[1], "+", color="#2ca02c", mew=2.6, ms=10, label="realer Sporn-Ort (Enthese)")ax.annotate("realer Sporn-Ort (Enthese)", xy=(SPORN[0] +0.7, SPORN[1] -0.6), xytext=(33.0, -5.0), fontsize=8, color="#2ca02c", va="center", arrowprops=dict(arrowstyle="->", color="#2ca02c", lw=1.1, connectionstyle="arc3,rad=-0.15"))ax.annotate("Spitzenort\n(stärkster Aufbau)", xy=(sz[0] -1.0, sz[1] +0.4), xytext=(1.0, 0.5), fontsize=8, color="#d62728", va="center", arrowprops=dict(arrowstyle="->", color="#d62728", lw=1.1, connectionstyle="arc3,rad=0.2"))ax.legend(loc="upper left", fontsize=7.5, framealpha=0.92, handlelength=1.6, borderpad=0.5, labelspacing=0.45)ax.text(2, -9.5, "hinten (Ferse)", fontsize=7, color="#666")ax.text(64, -9.5, "vorn (Zehen)", fontsize=7, color="#666", ha="right")ax.annotate("oben", xy=(-6, 24), xytext=(-6, 12), fontsize=7, color="#444", ha="center", arrowprops=dict(arrowstyle="->", color="#444", lw=1.1))# 10-mm-Maßstab auf eigener Höhe unter der Orientierungszeile.ax.plot([48, 58], [-14.0, -14.0], color="#333", lw=2.0)ax.text(53, -16.4, "10 mm", ha="center", fontsize=7)ax.set_xlim(-10, 70); ax.set_ylim(-19, 46); ax.set_aspect("equal"); ax.axis("off")plt.tight_layout(); plt.show()
Abbildung 18.1: Die Buchnetz-Kontur mit dem realen Sporn-Ort (grünes Kreuz, Enthese-Konturpunkt 24,07/3,35 mm am vorderen unteren Abhang, wie in Kapitel 1) und den beiden Ortszahlen des Zug-Aufbaus. Der Spitzenort (rotes Quadrat, stärkster Aufbau) liegt gut 2 mm vom realen Sporn-Ort und behält diese Lage auf jedem Buchnetz. Der Bandschwerpunkt (blaue Kette aus Kreis, Dreieck und Raute für grob, mittel und fein) liegt weiter weg und wandert mit der Netzfeinheit — er misst die Mitte eines langgezogenen Bandes, nicht dessen stärkste Stelle; grob und mittel liegen dabei nur 1,0 mm auseinander, deshalb stehen die Stufennamen in der Legende statt am Marker. Das Druck-Cluster sitzt weit oben an der Talus-Facette (hier nicht eingezeichnet — vergleiche das Zug- und Druckfeld am Umschalter oben), so auch in der Vergleichsrechnung. Sagittalschnitt, hinten links, vorn rechts, oben ↑; Maßstab 10 mm.
Deutung. Der stärkste Zug-Aufbau liegt 2,17 mm vom realen Sporn-Ort — die Zughypothese baut den Knochen dort, wo er tatsächlich sitzt, und diese Zahl bleibt über alle drei Buchnetze zwischen 1,9 und 2,5 mm. Die Druckhypothese baut ihn dagegen woanders, oben an der Talus-Facette. Der Ort trennt die beiden Hypothesen also.
Das Bild zeigt zugleich, warum man die Frage genau stellen muss. Die blauen Kreise sind derselbe Aufbau, nur anders vermessen: der Bandschwerpunkt, auf drei Netzstufen. Er liegt zwischen 5,2 und 10,1 mm vom Sporn und wandert dabei quer über das Fersenbein. Wer diese Zahl für den Aufbau-Ort hielte, müsste sagen: „Das Ergebnis hängt vom Netz ab.” Es hängt aber nicht das Ergebnis vom Netz ab, sondern die Kennzahl vom Ergebnis — sie misst die Mitte einer Zone, die auf feineren Netzen länger wird. Eine netzabhängige Zahl ist nicht immer ein netzabhängiges Ergebnis; manchmal ist sie nur die falsche Frage.
Wie diese Nähe zu lesen ist, steht oben im Befund-Abschnitt: Sie ist kein Beleg, sondern die Folge einer Annahme. Neben dem Ort bleibt deshalb die zweite, robustere Trennfrage wichtig — welches Muster jeder Fall legt. Ihr geht der Rest dieses Kapitels nach.
18.4 Das Cluster-Zentrum von Hand: was ein „Aufbau-Ort“ überhaupt ist
Damit die Ortszahl kein Orakel bleibt, rechnen wir ihre Definition an einem winzigen Beispiel nach — und zwar die des Bandschwerpunkts, denn er ist die Zahl, die die Vergleichsrechnung vorgibt und mit der der Ortsabstand zwischen den Hypothesen gemessen wird. Er ist nicht das eine Element mit der größten Zunahme (ein Einzelelement wäre netzabhängiges Rauschen; deshalb mittelt auch der Spitzenort über alle Elemente ab 95 % des Maximums, nicht über eines), sondern das flächengewichtete Zentrum des größten zusammenhängenden Clusters der stärksten Zunahmen. Drei Schritte stecken darin: eine Schwelle setzen (das oberste Zehntel der Zunahmen), die Elemente oberhalb der Schwelle zu zusammenhängenden Gruppen verbinden (über gemeinsame Kanten), und vom größten dieser Cluster das flächengewichtete Mittel der Schwerpunkte bilden. Die folgende Zelle führt genau das auf einem Hand-Gitter aus zwölf Elementen aus; das Buch prüft bei jedem Bau, dass dieser sichtbare Lehrcode und der schnelle Löser hinter den Reglern Ziffer für Ziffer dasselbe rechnen.
Das Cluster-Zentrum von Hand, Schritt für Schritt
# Das Cluster-Zentrum von Hand, auf einem 3x2-Hand-Gitter (zwölf Elemente). Erst die# Schwelle (oberstes Zehntel der Zunahmen), dann die groesste zusammenhaengende# Gruppe ueber die Nachbarliste, dann das flaechengewichtete Zentrum. Die Zunahmen# sind ein frei erfundenes Demo-Feld (keine Kanon- oder Messwerte).import hypothesen as werkzunahme = [0.02, 0.03, 0.04, 0.60, 0.05, 0.04,0.01, 0.02, 0.03, 0.62, 0.03, 0.02] # zwei helle Nachbarn: Elem. 3 und 9schwer = werk.HANDGITTER_SCHWERnachbarn = werk.HANDGITTER_NACHBARNflaeche = werk.HANDGITTER_FLAECHEzentrum, schwelle, groesse = werk.handnetz_cluster_python(zunahme, schwer, nachbarn, flaeche)
Tabelle 18.1: Das Cluster-Zentrum auf dem Hand-Gitter aus zwölf Elementen.
Größe
Wert
Schwelle (oberstes Zehntel)
0,545
größter Cluster
2 Elemente
Cluster-Zentrum
(1,5 / 1)
Deutung unter der Tabelle. Beim obersten Zehntel (die Schwelle liegt hier bei 0,545) bleiben von den zwölf Elementen nur die zwei hellsten übrig — Element 3 und Element 9, die über eine gemeinsame Kante zusammenhängen und so einen Cluster bilden. Ihr flächengewichtetes Zentrum liegt bei (1,5 / 1,0), genau zwischen den beiden hellen Elementen (das eine bei \(y = 0{,}5\), das andere bei \(y = 1{,}5\)) — nicht auf dem einen hellsten. Genau diese Rechnung, über die Hunderte Elemente des Buchnetzes ausgeführt, liefert die Cluster-Zentren der Bilder oben.
18.5 Die Mischung: ein Ort dazwischen, zwei Orte, oder springt es?
Reiner Zug und reiner Druck sind die beiden Extreme. Ein echter Schritt ist eine Mischung aus beidem. Was macht der Aufbau-Ort, wenn man den Reiz stufenweise vom reinen Zug zum reinen Druck dreht? Das lässt sich raten — oder rechnen. Der Regler mischt die beiden Lasten linear, \(f(\alpha) = (1-\alpha)\,f_{\text{Zug}} +
\alpha\,f_{\text{Druck}}\), und rechnet den Umbau unter der gemischten Last frisch (der Reiz wächst quadratisch mit der Last, also ist das keine Überblendung zweier fertiger Bilder, sondern jedes Mal eine neue Gleichgewichtsrechnung).
WichtigVorhersage-Punkt
Bevor du den Schieber bewegst: Wenn du von reinem Zug (\(\alpha = 0\)) zu reinem Druck (\(\alpha = 1\)) mischst — entsteht ein Aufbau-Ort dazwischen, der langsam von der Enthese zur Talus-Säule wandert? Bleiben es zwei Orte? Oder bleibt der Ort erst an der Enthese und springt dann irgendwann? Lege dich fest, dann schiebe.
Deutung. Die Antwort ist keine der drei glatten Erwartungen — und sie ist unruhiger, als man hofft. Der Aufbau-Ort wandert nicht stetig von der Spornregion zur Säule; er springt, und er springt sogar hin und her. Gemessen wird hier, wie angekündigt, mit dem Bandschwerpunkt: Es geht um den Wechsel zwischen zwei Zonen, und beide Zonen müssen mit demselben Maß gemessen werden. Bei reinem Zug (\(\alpha = 0\)) meldet er (29 / 12,2) mm — die Mitte des langen Sättigungsbandes, nicht dessen hellste Stelle. Bei \(\alpha\) = 0,25 rutscht sie hinunter an die Enthese, (21,6 / 6,2) mm. Bei \(\alpha\) = 0,5 sitzt sie oben auf der Säule, (25,2 / 20,4) mm. Bei \(\alpha\) = 0,75 ist sie wieder unten, (21,7 / 4,5) mm — und erst bei reinem Druck steht sie fest oben, (17,6 / 24,4) mm.
Dieses Hin und Her ist keine Panne der Rechnung, sondern die Kennzahl, die sich selbst entlarvt. Es gibt in Wahrheit zwei Aufbau-Zonen — das weiche Ansatzgewebe unten, die Lastpfad-Säule oben —, und die Ortszahl meldet stets nur die gerade größere. Solange ein hoher Zuganteil das Ansatzgewebe belastet, verknöchert es kräftig und stellt den größten Cluster; sobald der Druck überwiegt, gewinnt die Säule. Weil beide Zonen über weite Strecken ähnlich groß sind, kippt die Rangfolge mehrfach, und die eine Zahl hüpft mit. Das ist der Vorbehalt zur Ortszahl in seiner schärfsten Form: Auch wenn der Ort die beiden Reinformen sauber trennt, unterschlägt eine einzige Zahl für „den“ Aufbau-Ort, dass das Feld zwei Zonen hat — und wer sie für eine Wanderung liest, liest ein Ranglisten-Ergebnis als Bewegung. Bemerkenswert nebenbei: Das reine Zugfeld läuft schnell in die Sättigung — was direkt zur ersten Modellgrenze führt.
18.5.1 Den Mischungsgrad selbst drehen: warum es springt
Der Regler oben rechnet auf dem Buchnetz und zeigt das Ergebnis; warum der Ort springt statt zu wandern, bleibt darin verborgen. Zwei Zutaten stecken dahinter. Die erste macht aus dem Wechsel überhaupt einen Sprung: Es gibt zwei getrennte Reizzonen — unten die weiche Enthese, oben die Talus-Säule —, und der Aufbau-Ort ist immer nur das Zentrum der größeren. Kippt die Rangfolge, springt das Zentrum von der einen Zone zur anderen; ein gleitender Ort dazwischen kann so gar nicht entstehen. Die zweite Zutat entscheidet, wann der Wechsel kommt: Der Reiz wächst quadratisch mit der Last (halbierst du eine Last, viertelt sich der Reiz, den sie erzeugt), also trägt in der Mischung \(f(\alpha) = (1-\alpha)\,f_{\text{Zug}} + \alpha\,f_{\text{Druck}}\) der Zug mit \((1-\alpha)^2\), der Druck mit \(\alpha^2\) zum Reiz bei. Weil der Zugreiz von einem hohen Wert startet, bleibt er lange der größere — der Wechsel kommt erst spät. Der folgende Block macht beide Zutaten nachrechenbar, auf einem bewusst idealisierten Zwei-Zonen-Hand-Gitter: zwei feste Reiz-Zonen, quadratisch gemischt. Die energetische Überlagerung der beiden Lastfälle (den Kreuzterm) lässt er absichtlich weg, um die Mechanik pur zu zeigen — er bildet also nicht die frische Gleichgewichtsrechnung des Buchnetzes nach, sondern deren Prinzip. Er sucht das Cluster-Zentrum mit genau der Definition, die der Lehrcode oben eingeführt hat, und meldet die beiden stärksten Zonenreize und den Ort.
WichtigVorhersage-Punkt
Bevor du \(\alpha\) setzt: Du drehst den Mischungsgrad von 0 (reiner Zug) auf 1 (reiner Druck). Wandert das Cluster-Zentrum bei jedem Schritt ein Stück von unten nach oben — oder bleibt es unten und springt erst spät nach oben? Und bei welchem \(\alpha\) kippt es? Lege dich fest, dann probiere 0, 0,25, 0,50, 0,75, 0,87, 0,88 und 1,0 — die beiden mittleren Werte klammern den Kipppunkt ein.
Editierbare Zelle. Der folgende Block ist zum Selbstprobieren: Er mischt den Reiz mit deinem \(\alpha\) und zeichnet das Gitter samt Aufbau-Ort — der blaue Kreis wandert mit. Er trägt keinen neuen Gedanken; wer ihn überspringt, verliert nichts.
Code
fig, achsen = plt.subplots(3, 1, figsize=(4.6, 5.4), layout="constrained")tp = H18.misch_handgitter_bild(fig, achsen, [0.0, 0.88, 1.0])cb = fig.colorbar(tp, ax=list(achsen), fraction=0.046, pad=0.03, ticks=[0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 8.0], label="Demo-Reiz (willkürliche Einheiten), inferno logarithmisch 0,1–8,0")# Klarschrift statt 10^-1 / 10^0: die Leserin liest hier Zahlen, keine Exponenten.cb.ax.set_yticklabels([komma(w, 1) for w in (0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 8.0)])cb.ax.tick_params(labelsize=7)plt.show()
Abbildung 18.2: Das α-Reizfeld auf dem 3×2-Hand-Gitter bei drei Mischungen untereinander: α = 0 (reiner Zug), α = 0,88 (knapp über dem Kipppunkt) und α = 1 (reiner Druck). Der blaue Kreis ist das Cluster-Zentrum. Es bleibt bis α = 0,87 unten an der Enthese und springt bei α = 0,88 nach oben auf die Talus-Säule: ein Sprung, weil der Aufbau-Ort stets nur die größere von zwei getrennten Zonen meldet, und ein später Sprung, weil der Reiz quadratisch mit der Last wächst und der hohe Zuganteil bis zuletzt dominiert. Der Demo-Reiz überstreicht zwei Zehnerpotenzen (8,00 unten bei α = 0 gegen 0,20 oben bei α = 1), deshalb eine logarithmische inferno-Skala 0,1–8,0 statt der linearen des übrigen Kapitels; am Kipppunkt trennen die beiden Zonen weniger als zwei Prozent, weshalb ihr Zahlenwert zusätzlich im Bild steht. Idealisiertes Zwei-Zonen-Hand-Gitter (ohne energetische Überlagerung der Lastfälle), willkürliche Einheiten, keine Kanon- oder Buchnetzwerte; die untere Zellreihe ist die Enthese-Zone, die obere die Talus-Säulen-Zone — ein abstraktes Gitter, kein Fersenschnitt.
Deutung. Das ist die Antwort auf die Vorhersagefrage — und sie deckt sich qualitativ mit dem Buchnetz oben: Der Aufbau-Ort wandert nicht, er springt. Beide Zutaten liest man an den Zahlen ab. Dass es überhaupt ein Sprung ist und kein Gleiten, liegt an den zwei getrennten Zonen: Der Ort ist stets das Zentrum der größeren, und sobald die Rangfolge kippt, springt er mit ihr von unten nach oben. Der Kipppunkt liegt hier bei \(\alpha = 0{,}879\): Bei \(\alpha = 0{,}87\) ist der Enthese-Reiz mit 0,211 noch größer als der Säulen-Reiz mit 0,191, bei \(\alpha = 0{,}88\) ist es umgekehrt (0,193 gegen 0,195) — im Bild oben stehen beide Zahlen, denn zwei Prozent Unterschied zeigt keine Farbe. Wann er kippt, bestimmt die quadratische Gewichtung: Der Enthese-Reiz fällt mit \((1-\alpha)^2\), der Säulen-Reiz steigt mit \(\alpha^2\); weil der Zug von einem sehr hohen Wert startet, kippt es erst kurz vor dem reinen Druck. Auf dem Buchnetz oben, mit seinen zwei wirklich gerechneten Zonen, kippt die Rangfolge dagegen mehrfach — die dortigen Zonen sind einander viel ähnlicher als in diesem überzeichneten Demo. Der Mechanismus ist derselbe, nur die Kipppunkte sind andere. Und wieder zeigt sich: Es sind zwei Zonen, und die eine Ortszahl meldet immer nur die gerade größere.
18.6 Modellgrenze eins — das Zugfeld sättigt
Das reine Zugfeld läuft rasch in die Sättigung: Die weiche Enthese erreicht an ihrer hellsten Stelle die obere Dichtegrenze \(\rho_{\max} = 1{,}0\) g/cm³ und kann dort nicht weiter aufbauen. Danach kommt der Lauf zur Ruhe — auf jedem der drei Buchnetze, nur unterschiedlich schnell: 748 Iterationen auf dem groben, 1319 auf dem mittleren und 1165 auf dem feinen. Dass die feineren Netze länger brauchen, ist erwartbar: Je kleiner die Elemente, desto kleiner die Dichteänderung, die ein einzelnes von ihnen je Iteration beiträgt, und desto länger dauert es, bis die Summe unter die Abbruchschranke fällt.
WarnungEin Befund, der keiner ist: „wird nicht stationär“
Hier lauert eine Falle, in die man mit einer einzigen zu kleinen Zahl im eigenen Programm tappt. Jeder Lauf bricht irgendwann ab — spätestens bei einem festgelegten Iterationsbudget. Setzt man es zu knapp, hören die Läufe auf, bevor sie ausgelaufen sind, und in der Ausgabe steht dann ein Feld, das sich noch ändert. Wer das nicht prüft, deutet es als Eigenschaft des Modells: „Auf dem feinen Netz kriecht das Feld weiter und wird nicht stationär.”
Nachrechnen entlarvt es. Der Zuglauf braucht auf dem mittleren Netz 1319, auf dem feinen 1165 Iterationen, der Drucklauf auf dem feinen 1375. Mit einem Budget von 900 wären alle drei abgebrochen worden — und das Abbrechen hätte wie ein Befund über das Netz ausgesehen. Dieses Kapitel rechnet deshalb mit 3000; alle Läufe laufen aus.
Wer eine Rechnung abbricht, muss das Abbrechen vom Ergebnis unterscheiden. „Wird nicht stationär” und „war noch nicht fertig” sehen in der Ausgabe gleich aus — der Unterschied liegt allein darin, ob man noch einmal mit mehr Budget gerechnet hat.
Für Läufe, die tatsächlich nicht auslaufen, gibt es das Fenster-Kriterium (aus der Vergleichsrechnung übernommen): Ein Lauf gilt als auswertbar, sobald sein Cluster-Zentrum über die letzten 100 Iterationen um weniger als eine Elementkante wandert. Hier wandert es 0 mm — es wird gar nicht gebraucht. Abbildung 18.3 zeigt, wie die Enthese in die Sättigung läuft.
Code
s = PAKET["saettigung"]fig, ax = plt.subplots(figsize=(6.4, 3.6))ax.plot(s["iterationen"], s["enthese_mittel"], "-", color="#1f6feb", lw=1.8, label="flächengew. Enthese-Dichte (g/cm³)")ax.plot(s["iterationen"], [s["rho_enthese"] + d for d in s["dichte_max"]], "-", color="#e8710a", lw=1.8, label="Startdichte + stärkste Zunahme (g/cm³)")ax.axhline(s["rho_max"], ls="--", color="#888", lw=1.0)ax.text(s["iterationen"][-1], s["rho_max"] +0.01, "ρ_max = 1,0 g/cm³", ha="right", va="bottom", fontsize=8, color="#555")ax.set_xlabel("Iteration (Rechenschritt, keine Zeit)")ax.set_ylabel("Dichte (g/cm³)")ax.set_ylim(0.15, 1.08)ax.legend(fontsize=8, loc="center right")ax.grid(True, alpha=0.3)plt.tight_layout(); plt.show()
Abbildung 18.3: Der Zuglauf des feinen Buchnetzes über die Iterationen: die flächengewichtete mittlere Dichte der Enthese-Elemente (blau) und die stärkste Dichtezunahme im Feld (orange). Die stärkste Zunahme läuft rasch in die Sättigung (Δρ = 0,80, also von der Enthese-Startdichte 0,20 bis zur Obergrenze 1,00 g/cm³) und bleibt dort; die flächengewichtete Enthese-Dichte steigt langsamer, weil nicht alle Enthese-Elemente gleich stark belastet sind. Nach der Sättigung kommt der Lauf zur Ruhe — auf jedem Buchnetz, das feine braucht dafür nur länger. Die Achse zählt Iterationen (Rechenschritte), keine Wochen — eine Iteration ist keine Zeit.
Deutung. Die stärkste Zunahme im Feld erreicht schon nach wenigen Dutzend Iterationen die Obergrenze und bleibt dort — die Enthese ist an ihrer hellsten Stelle gesättigt. Die flächengewichtete Enthese-Dichte steigt gemächlicher, weil die Enthese-Region auch Elemente enthält, die weniger Last sehen. Sobald oben nichts mehr wachsen kann, zählt für die Aussage nur noch der Ort. Das ist keine Panne, sondern eine Modellgrenze mit einer klaren Lehre: Iteration ist keine Zeit. Das Modell sagt, wo der Sporn wächst, nicht wie lange er dazu braucht — die Zahl der Iterationen bis zur Sättigung ist ein Rechenweg, keine Aussage über Wochen oder Monate.
18.7 Modellgrenze zwei — die Druck-Säule strebt zur Talus-Kante
Die zweite Grenze betrifft das Druckfeld. Seine Säule strebt nach oben zur Talus-Facette — und die ist im Modell fest eingespannt (das Lager, das die Gelenkreaktion trägt). Ein fest eingespanntes Randstück konzentriert die Spannung, wie schon Kapitel 13 an der Lagerkante gezeigt hat: Es ist eine Lagerkantensingularität (keine Kerbe — die Kontur hat dort keinen einspringenden Winkel, sondern eine harte Randbedingung). Der Aufbau, den das Druckfeld am oberen Ende der Säule zeigt, kann durch diese Lagernähe überzeichnet sein. Abbildung 18.4 zeichnet die Lagerkante ein und misst den Abstand nach.
Code
g = PAKET["haupt"]netz = NETZE[PAKET["haupt"]["stufe"]]fig, ax = plt.subplots(figsize=(6.0, 4.9), layout="constrained")tp = zunahme_bild(ax, netz, g["druck"]["zunahme"], "Druck-Aufbau + Talus-Lager", zentrum_mm=g["druck"]["zentrum_mm"])kn = np.asarray(netz["knoten"]) *1000.0talus = np.asarray(netz["rand_talus"])tp_pts = kn[talus]ordnung = np.argsort(tp_pts[:, 0])# Die Lagerkante liegt AUF dem Netzrand: nach innen schwarzes Feld, nach außen weißer# Seitengrund. Eine reine weiße Linie verschwindet dort im Papier — deshalb ein# dunkler Saum unter der weißen Linie, damit sie nach beiden Seiten trägt.ax.plot(tp_pts[ordnung, 0], tp_pts[ordnung, 1], color="white", lw=2.2, solid_capstyle="round", path_effects=[pfx.Stroke(linewidth=4.8, foreground="#1f2937"), pfx.Normal()])ax.annotate("Talus-Lagerkante\n(fest eingespannt)", xy=(tp_pts[ordnung][len(ordnung)//2, 0], tp_pts[ordnung][len(ordnung)//2, 1]), xytext=(44, 38), fontsize=8, color="#222", arrowprops=dict(arrowstyle="->", color="#222", lw=1.0))fig.colorbar(tp, ax=ax, fraction=0.030, pad=0.03, label="Dichtezunahme Δρ (g/cm³), inferno 0–0,80")plt.show()# Abstand des Druck-Cluster-Zentrums zur nächsten Talus-Lagerfacette (berechnet).zentrum = np.array(g["druck"]["zentrum_mm"])facetten = netz["facetten_talus"]def _abstand_segment(p, a, b): ab = b - a tt = np.clip(np.dot(p - a, ab) / np.dot(ab, ab), 0.0, 1.0)returnfloat(np.linalg.norm(p - (a + tt * ab)))d_lager =min(_abstand_segment(zentrum, kn[int(f[0])], kn[int(f[1])]) for f in facetten)print("Abstand Druck-Cluster-Zentrum ↔ nächste Talus-Lagerfacette: %s mm"% komma(d_lager, 1))print("→ das Cluster-Zentrum klebt NICHT an der Lagerkante, sitzt ihr aber nahe:")print(" die Lagernähe bleibt ein Deutungshinweis (mögliche Überzeichnung, Kapitel 13).")
Abbildung 18.4: Das Druck-Aufbaufeld des feinen Buchnetzes (inferno 0–0,80 g/cm³) mit der fest eingespannten Talus-Lagerkante (weiße Linie mit dunklem Saum — sie liegt genau auf dem Netzrand und wäre als reine weiße Linie gegen den weißen Seitengrund nicht zu sehen) und dem Cluster-Zentrum (blauer Kreis). Die Lastpfad-Säule strebt nach oben zur Lagerkante; deren feste Einspannung konzentriert die Spannung (Lagerkantensingularität wie in Kapitel 13, keine Kerbe) und kann den Aufbau am oberen Säulenende überzeichnen — ein Deutungshinweis, kein Messfehler. Das Cluster-Zentrum selbst liegt aber rund 8,8 mm UNTER der Lagerkante (kein Cluster-Element berührt sie): Die Säule klebt nicht an der Kante. Sagittalschnitt, hinten links, vorn rechts, oben ↑; Maßstab 10 mm.
Abstand Druck-Cluster-Zentrum ↔ nächste Talus-Lagerfacette: 8,8 mm
→ das Cluster-Zentrum klebt NICHT an der Lagerkante, sitzt ihr aber nahe:
die Lagernähe bleibt ein Deutungshinweis (mögliche Überzeichnung, Kapitel 13).
Deutung. Die Druck-Säule strebt nach oben zur Lagerkante, und weil das Talus-Lager fest eingespannt ist, konzentriert es dort die Spannung — der Aufbau am oberen Säulenende kann also von der Lagernähe überzeichnet sein, so wie die Kerbspitze in Kapitel 13 unrealistisch hohe Spannungen zeigte. Das ist der Deutungshinweis, den auch die Vergleichsrechnung nennt, und er ist beim Ortsbefund besonders ernst zu nehmen: Der Druck-Ort sitzt gerade in dieser Lagerzone. Wichtig für die saubere Deutung der Ortszahl ist zweierlei. Erstens liegt das Cluster-Zentrum rund 8,8 mm unter der Kante, und kein Cluster-Element berührt das Lager — die Säule klebt also nicht an der Kante, es gibt kein unmittelbar an der Lagerkante lokalisiertes Cluster. Zweitens bleibt trotzdem offen, ob die feste Einspannung den Aufbau in dieser Zone verstärkt; deshalb ist der Druck-Ort ein Deutungshinweis, keine klinische Aussage. Ein Modell, das seine eigenen Verdachtsstellen ausweist und nachmisst, ist brauchbarer als eines, das eine glatte Zahl verspricht.
18.8 Netzstabilität: hält der Befund — und was macht das Buchnetz?
Ein Befund taugt nur, wenn er nicht am Netz hängt. Die Vergleichsrechnung ist deshalb auf zwei Netzstufen gefahren worden. Der folgende Regler zeigt allein sie (offline gerechnete Vergleichsdaten, wie die Netzleiter aus Kapitel 13) — keine Buchnetz-Zahlen.
Der Regler zeigt links die Übersicht mit zwei Lupenausschnitten und rechts, was in den Lupen steckt: das echte Netz um jedes der beiden Cluster. Gefüllt sind die Cluster-Elemente, der Ring ist ihr flächengewichtetes Zentrum — wieder der Bandschwerpunkt, denn gemessen wird der Abstand zwischen den beiden Zonen. Der Umschalter tauscht die Netzstufe: 3146 Dreiecke mit 1,03 mm Kante gegen 5320 mit 0,79 mm.
WichtigVorhersage-Punkt
Bevor du umschaltest: Auf dem feineren Netz sind die Dreiecke kleiner, und das Cluster wird aus mehr Elementen bestehen. Was macht der Aufbau-Ort dabei — springt er, wandert er langsam davon, oder bleibt er, wo er ist? Und wie weit dürfte er wandern, ohne dass man von einem netzabhängigen Ergebnis sprechen müsste? Lege dich fest, dann schalte um.
Deutung. In den Lupen ändert sich sichtbar viel und in der Sache nichts. Die Dreiecke schrumpfen, das Zug-Cluster wächst von 129 auf 214 Elemente, das Druck-Cluster von 147 auf 254 — es ist also nicht dasselbe Cluster, sondern eines aus feineren Bausteinen. Sein Zentrum aber rückt nur um 0,17 mm (Zug) beziehungsweise 0,16 mm (Druck) weiter, und das ist der Vergleich, auf den es ankommt: weniger als eine einzige Elementkante. Kleiner als ein Element kann ein Netz nichts auflösen; eine Verschiebung unterhalb dieser Größe ist Rechenrauschen, keine Ortsänderung. Der Abstand der beiden Cluster bleibt entsprechend bei rund 17 mm (17,33 → 17,58 mm), weit über der Schwelle von zwei Elementkanten. Die belastbare Aussage steht damit fest — der Zug baut in der Spornregion, und der Bodendruck baut woanders, oben an der Talus-Facette.
Das ist die Bedeutung des Wortes „netzstabil”: nicht, dass zwei Zahlen zufällig gleich aussehen, sondern dass sich das Ergebnis nicht ändert, während man die Bausteine austauscht, aus denen es gerechnet wird.
Und die anderen Buchnetze? Auch dort trennt der Befund die Hypothesen auf jeder Netzstufe — und in derselben Größenordnung wie die Vergleichsrechnung. Abbildung 18.5 zeigt das Druckfeld auf drei Buchnetzen, und zwar mit eingezeichneten Dreieckskanten: Von 585 über 1416 auf 4186 Dreiecke wird die Auflösung mehr als versiebenfacht — man sieht der Reihe also an, worüber sie eine Aussage macht. Das Ergebnis ändert sich dabei kaum: Die Lastpfad-Säule bleibt auf jedem Netz oben, und der Buchnetz-Ortsabstand bleibt rund 17 mm (grob 16,8, mittel 17,3, fein 16,7 mm). Er stimmt in der Größenordnung mit der Vergleichsrechnung überein — es ist also nicht bloß die Eigenheit einer Modellfamilie: Beide Rechnungen, die gröbere Buchnetz-Demo und die feinere Vergleichsrechnung, trennen den Druck (oben an der Talus-Facette) vom Zug (unten in der vorderen Tuberhälfte) um rund 17 mm.
Code
stufen = PAKET["regler_netz"]["stufen"]# Die Panels sind seitenverhältnistreu (80 mm breit, 61 mm hoch); Höhe und Breite# bleiben auf das Verhältnis der drei gestapelten Panels abgestimmt, sonst wird der# Überschuss zu leerem Rand statt zu größeren Knochen. Die Figur ist bewusst groß:# Quarto leitet die Anzeigebreite aus der Zollangabe ab, und in dieser Reihe muss# man einzelne Dreieckskanten unterscheiden können.fig, achsen = plt.subplots(len(stufen), 1, figsize=(7.4, 13.0), sharex=True, sharey=True, layout="constrained")for ax, s inzip(achsen, stufen): tp = zunahme_bild(ax, NETZE[s["stufe"]], s["druck_zunahme"],# Titel bewusst kurz: Die Dreieckszahl gehört zum sichtbaren# Netz, die Knotenzahl steht in der Ausgabe unter dem Bild.# Ein längerer Titel liefe in die Spalte der Farbleiste."%s · %d Dreiecke · Ortsabstand %s mm"% (s["stufe"], len(NETZE[s["stufe"]]["dreiecke"]), komma(s["abstand_mm"], 1)), zentrum_mm=s["druck_zentrum_mm"], netz_kanten=True, schrift=11)gemeinsame_colorbar(fig, tp, list(achsen), schrift=11)plt.show()print("Zwei netzstabile, übereinstimmende Ortszahlen — beide Familien trennen den Ort:")print(" Vergleichsrechnung (belastbar): Abstand %s → %s mm (fein → feiner)"% (komma(PAKET["referenz"]["abstand_mm"], 2), komma(PAKET["referenz"]["abstand_feiner_mm"], 2)))for s in stufen:print(" Buchnetz %-6s%4d Kn.: Abstand %s mm"% (s["stufe"], s["n_knoten"], komma(s["abstand_mm"], 1)))
Abbildung 18.5: Buchnetz-Beobachtung (nicht die Vergleichsrechnung): das Druck-Aufbaufeld auf dem groben, mittleren und feinen Buchnetz (inferno 0–0,80 g/cm³), mit Druck-Cluster-Zentrum (blauer Kreis) und Buchnetz-Ortsabstand je Stufe. Die weißen Linien sind die Dreieckskanten des jeweiligen Netzes — von 585 groben Dreiecken oben bis zu 4186 feinen unten wächst die Zahl auf mehr als das Siebenfache, und man sieht der Reihe damit an, worüber sie eine Aussage macht. Der Aufbau bleibt trotzdem derselbe: Die Säule steht auf jedem Netz oben, der Ortsabstand bleibt rund 17 mm und stimmt mit der Vergleichsrechnung (17,33 mm) überein — netzstabil, und beide Modellfamilien erzählen dasselbe. Sagittalschnitt, hinten links, vorn rechts, oben ↑; Maßstab 10 mm. Buchnetz-Rechnung.
Zwei netzstabile, übereinstimmende Ortszahlen — beide Familien trennen den Ort:
Vergleichsrechnung (belastbar): Abstand 17,33 → 17,58 mm (fein → feiner)
Buchnetz grob 330 Kn.: Abstand 16,8 mm
Buchnetz mittel 768 Kn.: Abstand 17,3 mm
Buchnetz fein 2199 Kn.: Abstand 16,7 mm
Deutung. Zwei sauber gerechnete Modelle, und ihr Ortsabstand stimmt überein: rund 17 mm in der Vergleichsrechnung wie auf dem Buchnetz, auf jeder Netzstufe. Netzstabil ist die Trennung der beiden Lastfälle — der Druck-Ort bleibt oben an der Talus-Facette, der Zug unten in der vorderen Tuberhälfte; die genaue Zug-Ortszahl rückt auf dem feinsten Buchnetz etwas vor, ihr Abstand zum Druck bleibt. Diese Trennung ist robust. Und dennoch trägt sie den Streit nicht allein — sie kann es gar nicht, denn der Zug-Ursprung wurde an den Sporn-Ort gelegt, die Lagernähe des Druck-Clusters bleibt ein Deutungsvorbehalt, und der Druckfall prüft die Bodenlast, nicht die örtliche Kompression der Literatur. Deshalb bleibt neben dem Ort das Muster die zweite Trennfrage: Sättigungsfleck am Abhang gegen Lastpfad-Säule.
Und die dritte, die eigentlich interessante, bleibt eine Frage: was eine Einlage an diesem Muster ändern würde. Sie zu beantworten hieße, in die Reizverteilung einzugreifen — und dafür müsste die Auflagefläche aus der Rechnung kommen statt aus der Vorgabe. Dieses Buch stellt die Frage scharf; es beantwortet sie nicht.
18.9 Die Blackbox geöffnet
In Kapitel 1 stand ein Kasten mit vier Reglern. Du hast an ihnen gedreht, ein Bild ist erschienen, und warum es so aussah, wusste niemand außer dem Programm. Genau das war die Verabredung: erst bedienen, dann verstehen. Jetzt ist das Buch zu Ende, und der Kasten steht noch einmal da — derselbe Regler, dasselbe Vorführmodell. Was sich geändert hat, ist, dass jeder seiner Knöpfe inzwischen einen Namen hat, ein Kapitel und eine Rechnung.
WichtigVorhersage-Punkt
Bevor du noch einmal drehst: Drei der vier Regler greifen an einer Größe an, die dieses Buch inzwischen wirklich ausrechnet; einer greift an einer Größe an, die jedes Modell dieses Buches vorgibt. Welcher ist der vierte — und woran hättest du das schon in Kapitel 1 merken können? Lege dich fest, dann dreh.
Deutung. Der vierte ist die Einlage. Die drei anderen Regler ändern Lasten, und Lasten sind seit Kapitel 2 die gegebene Größe dieses Buches; die Einlage dagegen ändert, wie der Fuß den Boden berührt — und die Berührung ist die eine Größe, die kein Modell dieses Buches ausrechnet. Der Reihe nach:
Und das Bild selbst, das die vier Regler erzeugen, ist inzwischen ebenfalls kein Zauberstück mehr: Es ist ein Feld über einem Netz aus Dreiecken (Kapitel 10), das aus einem Gleichungssystem kommt (Kapitel 5) und aus dem Gleichgewicht in jedem Punkt (Kapitel 8, Kapitel 9); seine Farbe hat eine Skala, eine Einheit und eine Orientierung (Kapitel 11, Kapitel 15); und ob man ihm trauen darf, entscheidet nicht seine Schönheit, sondern eine Netzstudie (Kapitel 13). Ein Unterschied bleibt: Das Vorführmodell aus Kapitel 1 gibt den Ort der Verdichtung noch vor. Dieses Kapitel hat ihn ausgerechnet.
Vier Fragen, die du jetzt selbst stellen kannst. Wer eine Rechnung lesen kann, kann auch sagen, was ihr fehlt. Vier Lücken dieses Buches sind so groß, dass sie eigene Aufgabenklassen sind — und du kannst sie inzwischen genau benennen:
Die Berührung mit dem Boden. Wie groß eine Auflagefläche ist, hängt davon ab, wie stark man drückt; wer eine Einlage rechnen will, muss Fläche und Druck gleichzeitig finden. Das ist ein Kontaktproblem, und es ist die offene Frage dieses Buches.
Die dritte Dimension. Alles hier ist ein Sagittalschnitt — eine Scheibe von 33 mm Dicke. Ein Fersenbein ist ein Körper, und ein Sporn hat auch eine Breite. Was aus der Ebene heraus geschieht, sagt dieses Buch nicht.
Der Vergleich mit Messdaten.Kapitel 13 hat Verifikation und Validierung getrennt: Geprüft ist hier, dass richtig gerechnet wird, nicht dass das Richtige gerechnet wird. Ein Vergleich mit gemessenen Dichteverteilungen realer Fersenbeine steht aus — und wäre der nächste Schritt.
Die Richtung. Der Knochen dieses Buches ist isotropes Material: in jeder Richtung gleich steif. Die Trabekel aus Kapitel 15 haben aber eine Richtung, und ein Umbaugesetz, das nur die Dichte ändert, kann sie nicht drehen. Ein Modell, das auch die Richtung umbaut, ist eine eigene Klasse.
Keine dieser vier Fragen ist eine Ausrede. Sie sind der Ertrag: Am Anfang stand ein Kasten, den man nur bedienen konnte. Am Ende steht ein Kasten, dessen Grenzen man benennen kann — und das ist der Unterschied zwischen einem Regler und einem Modell.
Übungen
Ü 18.1 (Verstehen). Erkläre in einem Satz, warum man beide Hypothesen unter einem gemeinsamen Referenzreiz \(S_0\) rechnen muss und nicht unter je einem eigenen.
HinweisMusterlösung zu Ü 18.1
Mit je einem eigenen \(S_0\) (dem Mittel des jeweiligen Reizfelds) rechnete man beide Szenarien auf ihr eigenes Gleichgewicht herunter und normalisierte damit gerade die Lastunterschiede weg — der Reiz wächst quadratisch mit der Last, sodass ein starker und ein schwacher Lastfall nach ihrer je eigenen Normierung gleich aussähen. Nur ein gemeinsamer Maßstab macht die beiden Hypothesen überhaupt vergleichbar; sonst vergliche man zwei verschieden geeichte Waagen.
Ü 18.2 (Verändern). Öffne den Zug/Druck-Umschalter. Lies für beide Felder ab, wo der stärkste Aufbau sitzt, und schätze den Ortsabstand auf dem Buchnetz. Vergleiche ihn mit den 17,33 mm der Vergleichsrechnung. Was bedeutet es, dass beide Rechnungen denselben Ortsabstand finden — und was folgt daraus nicht?
HinweisMusterlösung zu Ü 18.2
Auf dem feinen Buchnetz sitzt der stärkste Zug-Aufbau am vorderen unteren Abhang, gut zwei Millimeter vom realen Sporn-Ort (rund 23,6 / 5,5 mm), der Druck-Aufbau oben an der Talus-Säule (rund 17,6 / 24,4 mm) — der Ortsabstand der beiden Bandschwerpunkte ist rund 17 mm. Die Vergleichsrechnung misst ganz ähnlich 17,33 mm, mit demselben Bild (Zug unten in der Spornregion, Druck oben). Dass zwei verschieden feine Rechnungen übereinstimmen, macht den Ortsbefund robust: Der Ort trennt die beiden Lastfälle. Daraus folgt aber nicht, dass Zug die Ursache des Sporns ist: Eine Zuglast baut an ihrem eigenen Ansatz auf — das ist Mechanik und keine Entdeckung; das Druck-Cluster sitzt lagernah (mögliche Überzeichnung); und der gerechnete Druck ist die Bodenlast der Standphase, nicht die örtliche Kompression des Ansatzgewebes, von der Kumai und Benjamin sprechen. Ihre klinische Gegenposition bleibt unwiderlegt. Die Lehre: Der Ort trennt, aber die tragende Aussage ist die über den Druck — er baut woanders.
Ü 18.3 (Übertragen). Das reine Zugfeld läuft in die Sättigung. Begründe, woran das liegt, und sage, was das für eine Aussage über die Dauer des Sporn-Wachstums bedeutet.
HinweisMusterlösung zu Ü 18.3
Die weiche Enthese läuft unter reinem Zug an ihrer hellsten Stelle in die obere Dichtegrenze \(\rho_{\max} = 1{,}0\) g/cm³ (Sättigung) und kann dort nicht weiter aufbauen. Danach kommt der Lauf auf jedem Buchnetz zur Ruhe; die feineren brauchen nur mehr Iterationen dafür (748 / 1319 / 1165), weil jedes einzelne kleine Element je Iteration weniger beiträgt. Für die Dauer bedeutet das: gar nichts. Eine Iteration ist kein Tag und keine Woche — das Modell sagt, wo der Sporn wächst, nicht wie lange er dazu braucht. Die Zahl der Iterationen bis zur Sättigung ist ein Rechenweg, keine biologische Zeitspanne.
Roter Faden
Wo kam das schon vor, wo kommt es wieder?Zurück: Die getrennten Zug- und Druckfelder sind die Fortsetzung von Kapitel 15 (dort die Hauptspannungen, hier der Umbau, den sie treiben); die geschlossene Schleife und der Reizfilter stammen aus Kapitel 17, und mit ihnen die Warnung, einer Ortszahl nicht blind zu trauen, ohne sie über die Netzstufen geprüft zu haben; die Totzone, die den Umbau zur Ruhe bringt, ist die aus Kapitel 16; die Lagerkante, die den Druckaufbau überzeichnen kann, ist die Singularität aus Kapitel 13; und der Widerlegungskasten als Idee stand schon in Kapitel 1; und die vier Regler, die das Buch eröffnet haben, stehen im Abschnitt „Die Blackbox geöffnet“ noch einmal da, jeder mit dem Kapitel, das ihn erklärt hat. Vor: Von hier führte der Weg zu einem Eingriff in genau diese Reizverteilung — zu der Frage, was eine Einlage ändert. Sie verlangt eine Auflagefläche, die aus der Rechnung kommt statt aus der Vorgabe, und damit ein Kontaktproblem. Das liegt jenseits dieses Buches; benennen kannst du es jetzt.
Teil VI abgeschlossen — die Studiums-Tabelle
Damit schließt Teil VI, „Der Knochen antwortet“ — und mit ihm die Studiums-Tabelle des Vorworts. Teil VI hat sie um zwei Zeilen erweitert, und beide sind eingelöst: Die Totzone im Mechanostat (Kapitel 16) ist als Totband eines Regelkreises und nichtlineare Rückkopplung erschienen; das „nochmal rechnen, bis sich nichts mehr ändert“ (Kapitel 17) ist die Fixpunktiteration mit ihrer Konvergenz. Dieses Kapitel hat beide auf den Forschungsstreit angewandt — die Totzone und die Sättigung machen den Umbau auswertbar, die Iteration liefert den Aufbau-Ort.
Weil dies das letzte Teilende ist, steht in der dritten Zeile der Rest: die elf Zeilen aus Teil I bis V, die dieses Kapitel nicht mehr einführt, sondern alle zugleich benutzt.
Im Buch als …
Im Studium heißt es …
Stand nach Teil VI
Totzone im Mechanostat
Totband eines Regelkreises, nichtlineare Rückkopplung
eingeführt in Kapitel 16, hier auf Zug und Druck angewandt
„nochmal rechnen, bis sich nichts mehr ändert“
Fixpunktiteration, Konvergenz
eingeführt in Kapitel 17, hier bis zur Sättigung und zum Fenster-Kriterium geführt
Kraftpfeil · Tabelle · „alle Knoten gleichzeitig“ · Steigung · Fläche · Gleichgewicht in jedem Punkt · „die Gleichung abschwächen“ · Netz aus Dreiecken · festhalten, ziehen, betten · Zug und Druck trennen · „das Netz war zu grob“
Vektor · Matrix · lineares Gleichungssystem · Ableitung · Integral · partielle Differentialgleichung · schwache Form · Diskretisierung · Randbedingungen · Eigenwerte und Hauptachsen · Konvergenzstudie
in Teil I bis V eingelöst, in diesem Kapitel alle elf zugleich im Einsatz
Hier endet der Weg dieses Buches. Er hat bei einer Läuferin mit Fersenschmerz begonnen und führt zu einem Modell, das sagen kann, wo es rechnet, wie genau, und was es nicht weiß. Die Frage des Buchtitels — warum der Knochen wächst, wo er wächst — hat in diesem Kapitel ein Bild bekommen: Abbildung 18.1 legt den Ort des stärksten berechneten Aufbaus neben den realen Sporn-Ort, und die beiden liegen 2,17 mm auseinander. Das ist die Antwort, die dieses Buch geben kann: Der Knochen wächst dort, wo der Reiz am größten ist — und das Modell zeigt, wo es diesen Reiz sieht, unter welchen Annahmen und mit welchem Vorbehalt. Die Frage, die als nächste käme — was ändert eine Einlage? —, verlangt eine Auflagefläche, die das Modell selbst ausrechnet. Sie ist gestellt und bleibt stehen; und wer bis hierher gelesen hat, kann sagen, warum sie offen ist.
Was dieses Kapitel NICHT tut
Es entscheidet den Streit nicht: Weder Zug noch Druck wird zur Ursache des Sporns erklärt. Dass der Ort die beiden Lastfälle trennt (Zug in der Spornregion, Druck oben), ist kein unabhängiger Beleg der Zughypothese: Dass eine Zuglast an ihrem eigenen Ansatz die höchste Spannung erzeugt, ist Mechanik und keine Entdeckung — der Zugaufbau dort konnte gar nicht ausbleiben; und das Druck-Cluster sitzt lagernah, was seinen Ort überzeichnen kann. Vor allem aber prüft es die Druckhypothese nicht in der Lesart ihrer Urheber: Gerechnet ist der Bodendruck der Standphase, nicht die örtliche Kompression des Ansatzgewebes. Die klinische Gegenposition von Kumai und Benjamin (der Sporn als Antwort auf vertikale Druck- statt Zugbelastung) wird durch dieses Modell also nicht widerlegt — sie bleibt die reale Referenz. Es macht keine klinische Aussage über den Sporn eines realen Menschen; die Orte, Abstände und Dichten sind Modellwerte unter dokumentierten Annahmen (Geometrie aus der bildbasierten Kontur, angenommene reine Einzellasten, weiche Enthese als Materialregion). Es kennt keine biologische Umbauzeit: Eine Iteration ist keine Zeit — das Modell sagt wo, nicht wann. Es prüft nicht die in der Literatur diskutierte Mikrofraktur-Reparatur-Hypothese (der Weg „Mikroschaden → Reparatur → Anbau“ ist im Modell nicht abgebildet). Und es lässt keinen Rand wandern: Der Sporn entsteht hier als Dichtezunahme eines weichen Ansatzgewebes, nicht als Verschiebung der Knochenkontur. Ein Modell, das die Kontur selbst wandern lässt, ist eine eigene Aufgabenklasse und in diesem Buch nicht enthalten. Es kennt auch keine Verknöcherungsschwelle: Im Modell gibt es allein eine Dichte zwischen 0,20 und 1,00 g/cm³ und eine Totzone um den Referenzreiz, innerhalb derer nichts geschieht — es rechnet Dichte um, nicht Gewebe. Ab welcher Belastung aus faserknorpeligem Ansatzgewebe tatsächlich Knochen wird, ab wann die Verknöcherung also einsetzt, ist damit nicht beantwortet; die Totzone ist eine Schwelle des Umbaugesetzes und keine biologische. Schließlich sind die Buchnetz-Zahlen dieses Kapitels (Ortsabstand ~17 mm, \(S_0 \approx 4\,8\) J/kg, die Aufbau-Orte) Modellwerte dieser Buchnetz-Rechnung; die belastbaren Hypothesentest-Kennzahlen (Ort 17,33 mm, Vollfeld-Differenz 0,746, \(S_0\) 4,84 J/kg) stammen aus der bildbasierten Vergleichsrechnung und werden hier zitiert, nicht neu erzeugt. Dass beide \(S_0\) nur um rund zwei Prozent auseinanderliegen, ist ein Nebenbefund und kein Beleg: Es sind weiterhin zwei verschiedene Modelle mit zwei verschiedenen Konturen.