A mathematical model for Nordic skiing

Oorspronkelijke auteurs: Jane Shaw MacDonald, Rafael Ordoñez Cardales, John M. Stockie

Gepubliceerd 2026-05-19

📖 5 min leestijd🧠 Diepgaand

Oorspronkelijke auteurs: Jane Shaw MacDonald, Rafael Ordoñez Cardales, John M. Stockie

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je voor dat je probeert te voorspellen hoe snel een langlaufer precies een wedstrijd zal voltooien. Het gaat niet alleen om de kracht van hun benen; het is een complexe dans tussen hun spieren, de vorm van het met sneeuw bedekte parcours, zwaartekracht, wind en zelfs hoe ze bochten nemen.

Dit artikel is als een wiskundig receptenboek voor het simuleren van die wedstrijd. De auteurs, wiskundigen en wetenschappers, bouwden een computerprogramma dat fungeert als een "virtuele langlaufer" om te zien hoe verschillende factoren het resultaat beïnvloeden. Hier is hoe ze dit deden, uitgelegd in eenvoudige termen:

1. De kaart tekenen (Het parcours)

Echte ski-parcours zijn geen perfecte rechte lijnen; het zijn kronkelige, hobbelige 3D-paden. Meestal hebben we slechts enkele verspreide GPS-punten (zoals stippen op een kaart) om het parcours te beschrijven.

Het probleem: Als je die stippen gewoon met rechte lijnen verbindt (zoals een kind stippen op een pagina verbindt), ziet het pad er gekarteld en onrealistisch uit. Als je probeert het glad te strijken met standaard wiskundige krommen, ontstaan er soms "spookheuvels" of dalen die in werkelijkheid niet bestaan (zoals een wiebelige tekening).
De oplossing: De auteurs gebruikten een speciaal type wiskundige gladmaking genaamd een Hermite-spline. Denk hierbij aan een flexibele liniaal die perfect door de GPS-punten buigt zonder nep-heuveltjes te creëren. Het creëert een glad, realistisch pad waarover hun virtuele langlaufer kan reizen.

2. De fysica van de virtuele langlaufer (De motor)

Zodra het pad is getekend, plaatsen ze een "virtuele langlaufer" erop. Deze langlaufer wordt geleid door de wetten van de natuurkunde (de wetten van Newton), die de auteurs omzetten in een reeks vergelijkingen.

De krachten: De langlaufer wordt geduwd en getrokken door vier hoofdcomponenten:
1. Spierkracht: De langlaufer duwt zich vooruit. Het model gaat ervan uit dat ze het hardst duwen bij het beklimmen van hellingen (langzaam) en meer coaten bij het afdalen (snel).
2. Zwaartekracht: De zwaartekracht trekt hen de hellingen af (waardoor ze versnellen) en houdt hen tegen bij het beklimmen van hellingen (waardoor ze vertragen).
3. Wrijving: De sneeuw wrijft tegen de ski's, waardoor ze vertragen.
4. Luchtwrijving: De lucht duwt tegen hen aan, vooral wanneer ze snel gaan.
De wiskunde: Ze losten deze vergelijkingen op met een geavanceerde rekenmachine (een computergenerator) die zijn snelheid aanpast om het antwoord exact goed te krijgen, zelfs wanneer het terrein lastig wordt.

3. De 3D-twist (Bochten nemen en remmen)

De meeste eerdere modellen keken alleen naar de wedstrijd van de zijkant (2D), alsof je naar een film op een plat scherm kijkt. Maar echt skiën gebeurt in 3D.

De nieuwe functie: De auteurs voegden de mogelijkheid toe voor de langlaufer om links en rechts te draaien. Wanneer een langlaufer scherp een bocht neemt tijdens een afdaling, moet hij remmen om niet van het parcours te vliegen.
De analogie: Stel je voor dat je met een auto om een scherpe bocht rijdt. Als je te snel gaat, slip je. De langlaufer moet "schuiven" of "stappen" om te vertragen. Het model berekent deze "remkracht". Ze ontdekten dat hoe een langlaufer een bocht neemt, enkele seconden kan toevoegen aan of aftrekken van hun totale tijd – een enorm verschil in een wedstrijd waar winnaars vaak door fracties van een seconde worden gescheiden.

4. Het model testen

Het team testte hun virtuele langlaufer tegen real-world data:

De "baseline"-test: Ze draaiden een simulatie op een parcours van 4,2 km en vergeleken dit met echte wedstrijdtijden. Hun model was ongelooflijk nauwkeurig en kwam binnen enkele seconden overeen met de echte resultaten.
De "elite"-test: Ze simuleerden een wedstrijd van 15 km met 36 verschillende echte atleten. Door de instelling "spierkracht" in hun computer aan te passen, konden ze perfect de finishtijden van trage skiërs, snelle skiërs en zelfs de wedstrijdwinner matchen.
De vermoeidheidsfactor: Ze merkten op dat echte skiërs aan het einde van een lange wedstrijd vertragen omdat ze moe worden. Hun basismodel hield hier geen rekening mee, dus ze toonden aan hoe je een "vermoeidheidsschakelaar" kunt toevoegen om de virtuele langlaufer langzamer te laten worden naarmate de wedstrijd vordert.

Waarom dit belangrijk is

De auteurs zeggen dat dit niet alleen voor sportliefhebbers is. Ze ontwierpen dit artikel om te laten zien dat wiskunde die je op de universiteit leert (zoals calculus en computerprogrammeren) echte, rommelige problemen kan oplossen.

Het bewijst dat het gebruik van een gladder, nauwkeuriger kaart (de spline) betere resultaten oplevert dan het gebruik van een gekartelde, eenvoudige kaart.
Het toont aan dat 3D-effecten (zoals bochten nemen en remmen) cruciaal zijn voor het begrijpen hoe topsporters winnen.
Het biedt gratis, open-source computercode die coaches, wetenschappers en studenten kunnen gebruiken om te experimenteren met verschillende wedstrijdstategieën.

Kortom, het artikel bouwt een digitale tweeling van een langlaufer. Het neemt een ruwe kaart, past de wetten van de natuurkunde toe en simuleert een wedstrijd zo nauwkeurig dat het ons helpt de kleine details te begrijpen – zoals hoe een langlaufer een bocht neemt – die het verschil kunnen betekenen tussen goud en zilver.

Technische Samenvatting: Een Wiskundig Model voor Noordse Ski

Probleemstelling
Noordse (langlauf)skiën presenteert een complex dynamisch systeem waarbij een atleet een driedimensionale ruimtekromme navigeert, gekenmerkt door frequente veranderingen in hoogte, richting en kromming. In tegenstelling tot alpineskiën, waar de zwaartekracht voornamelijk de beweging aandrijft, moeten Noordse skiërs voortstuwende krachten genereren om voorwaartse beweging op vlakke stukken en hellingen te handhaven, terwijl ze weerstand van sneeuwfrictie en aerodynamische drag managen. Bestaande modellen, zoals het door Moxnes, Sandbakk en Hausken (MSH) voorgestelde, hebben gewone differentiaalvergelijkingen (ODE's) gebruikt om de dynamiek van skiërs op tweedimensionale (2D) parcoursen te beschrijven. Deze modellen vertrouwen echter vaak op stuksgewijs lineaire benaderingen van de parcoursgeometrie, die de gladde kromming van echte banen niet vastleggen en numerieke artefacten introduceren. Bovendien negeren standaard 2D-modellen de kritieke 3D-effecten van bochten, specifiek de remkrachten die nodig zijn om strakke afdalingbochten veilig te nemen. Dit artikel adresseert de behoefte aan een realistischer, hoogwaardig wiskundig model dat gladde parcoursinterpolatie, rigoureuze krachtenbalans en 3D-bochtdynamica integreert, terwijl het toegankelijk blijft voor bachelorstudenten in calculus en wetenschappelijk rekenen.

Methodologie
De auteurs ontwikkelen een uitgebreid modelleringkader gebaseerd op Newtoniaanse mechanica, geformuleerd als een stelsel niet-lineaire gewone differentiaalvergelijkingen (ODE's). De methodologie verloopt via verschillende sleutelfasen:

Parcoursgeometrie en Interpolatie:
- Het artikel gaat voorbij aan de stuksgewijs lineaire representaties die in eerdere werken werden gebruikt. In plaats daarvan wordt Cubische Hermite-spline-interpolatie (specifiek de makima-variant in MATLAB) ingezet om gladde, vormbehoudende krommen te construeren uit schaarse GPS-data of 2D-hoogteprofielen.
- Deze aanpak vermijdt het "Runge-fenomeen" (schijnbare oscillaties) dat geassocieerd is met standaard kubische splines bij het aanpassen van onregelmatig gespatieerde data.
- Het parcours wordt geparametriseerd door geprojecteerde booglengte ( $\xi$ ) om meetkundige grootheden zoals helling en kromming nauwkeurig weer te geven.
- Het model onderscheidt tussen 2D-simulaties (alleen hoogteprofiel) en 3D-simulaties (volledige ruimtelijke coördinaten $x, y, z$ ), waarbij hellingshoeken ( $\theta$ ) en azimuthhoeken ( $\phi$ ) worden berekend uit spline-afgeleiden.
Besturende Vergelijkingen (Dynamica):
- 2D-model: De skiër wordt behandeld als een puntmassa. Het stelsel lost snelheid ( $v$ ) en geprojecteerde afstand ( $\xi$ ) op met behulp van twee gekoppelde ODE's:
  $v' = \frac{P(v)}{mv} - g \sin \theta - \mu g \cos \theta - \beta v^2$
  $\xi' = v \cos \theta$
  Hierbij omvatten de krachten voortstuwing ( $F_p$ ), zwaartekracht ( $F_g$ ), sneeuwfrictie ( $F_s$ ) en aerodynamische drag ( $F_d$ ). De voortstuwingskracht is afgeleid van een snelheidsafhankelijke vermogensfunctie $P(v)$ , en drag wordt gemodelleerd als kwadratisch in snelheid met een coëfficiënt die schakelt tussen "rechtop" en "tuck"-posities op basis van snelheidsdrempels.
- 3D-extensie: Het model integreert een remkracht ( $F_b$ ) om rekening te houden met middelpuntvliedende versnelling in strakke bochten. Deze kracht is evenredig met de middelpuntvliedende versnelling ( $\kappa v^2$ ) en wordt alleen geactiveerd wanneer de skiër zich op een afdaling bevindt ( $\theta < 0$ ) en de versnelling een specifieke drempel overschrijdt. De remcoëfficiënt ( $\gamma$ ) varieert op basis van de gebruikte techniek (schuifbocht versus stapbocht).
Numerieke Implementatie:
- De ODE-systemen worden opgelost met ode113 van MATLAB, een hoog-orde, variabele-stap, variabele-orde solver (Adams-Bashforth-Moulton) met adaptieve tijdstappen en foutcontrole.
- Booglengte-integralen voor afstand en kromming worden benaderd met de regel van Simpson (vierde-orde nauwkeurigheid).
- De code zorgt ervoor dat de skiër blijft vastgezet op het geïnterpoleerde spline-pad, waardoor "drift"-fouten worden geëlimineerd die veelvoorkomend zijn in uitgebreide ODE-systemen die de verticale positie onafhankelijk oplossen.

Belangrijkste Bijdragen

Gladde Geometrische Representatie: Het artikel demonstreert dat het vervangen van stuksgewijs lineaire parcoursbenaderingen door Hermite-splines realistischere simulaties oplevert. Waar MSH kubische splines als minder nauwkeurig vond vanwege oscillaties, toont dit werk aan dat vormbehoudende Hermite-splines deze artefacten elimineren terwijl ze een gladdere, nauwkeurigere representatie van het parcours bieden.
3D-Remdynamica: De auteurs introduceren een nieuw 3D-component in Noordse ski-modellen: een remkrachtmodel afhankelijk van parcourskromming en skiersnelheid. Dit maakt het mogelijk om strategische remtechnieken (schuifbocht versus stapbocht) op steile, strakke afdalingen te simuleren, een factor die eerder werd overgeslagen in 2D-modellen.
Validatie en Kalibratie: Het model wordt gevalideerd tegen:
- De MSH-basisgegevens (4,2 km parcours), waarbij nauwe overeenkomst wordt getoond in finishtijden en snelheidsprofielen.
- Experimentele data van Welde et al. (15 km wedstrijd), waarbij succesvol wordt gedifferentieerd tussen "langzame", "snelle" en "winnaar"-skiërs door de maximale vermogensparameter ( $P_{max}$ ) aan te passen.
- Werkelijke FIS-gecertificeerde parcoursen (Ole- en Stephanie-parcoursen in Toblach, Italië), waarbij wordt geverifieerd dat het model voldoet aan homologatiecriteria (bijv. hellingslimieten, parcourslengte).
Educatief Kader: Het artikel structureert het model expliciet om elementaire concepten uit het bacheloronderwijs in calculus (vectorcalculus, afgeleiden, integralen) en wetenschappelijk rekenen (ODE-solvers, interpolatie) te benutten, dienend als brug tussen theoretische wiskunde en toegepaste sportwetenschap.

Resultaten

2D-simulaties: Op het MSH-parcours produceerde het Hermite-spline-model een finishtijd van 823 seconden, wat ligt tussen het resultaat van de lineaire spline (813 s) en het oscillerende kubische spline-resultaat (geschat 835 s). Het model slaagde erin de snelheidsprofielen en krachtenbalansen die in de basisstudie werden waargenomen, succesvol te repliceren.
Parametergevoeligheid: Het model bevestigde dat de massa van de skiër en het vermogensoutput kritieke differentiatoren zijn. Door $P_{max}$ aan te passen, kon het model de prestaties van elite-skiërs (Winnaar: $P_{max} = 434$ W) versus langzamere concurrenten ( $P_{max} = 346$ W) matchen.
Impact van 3D-remmen: Simulaties op het 4,2 km Ole-parcours toonden aan dat het toepassen van remkrachten de finishtijden met 10–30 seconden verhoogde, afhankelijk van de techniek en de drempel. De "schuifbocht"-techniek (hogere remcoëfficiënt) resulteerde in langzamere tijden dan de "stapbocht". Op het krommere "Stephanie"-parcours was het tijdsverlies door remmen nog prominenter (tot +63 seconden), wat de strategische belangrijkheid van bochtbeheer benadrukt.
Vermoeidheidsmodellering: Het artikel merkt op dat het constante-vermogensmodel de "positieve pacing" (vertragen tijdens de wedstrijd) die in echte wedstrijden wordt waargenomen, niet volledig kan vastleggen. Een eenvoudige modificatie waarbij $P_{max}$ in de loop van de tijd afneemt, werd voorgesteld om spiervermoeidheid na te bootsen, hoewel er geen volledig fysiologisch model werd ontwikkeld.

Betekenis en Claims
De auteurs positioneren dit werk als een demonstratie van hoe elementaire wiskundige concepten kunnen worden toegepast om wereldse problemen in de sportwetenschap op te lossen. Zij claimen dat hun model biedt:

Nieuwe Inzichten: Specifiek de kwantificering van hoe remtechnieken op strakke afdalingbochten de wedstrijdtijden significant beïnvloeden, een factor die vaak wordt over het hoofd gezien in 2D-analyses.
Educatieve Waarde: Het model dient als een concreet voorbeeld voor bachelorcursussen in differentiaalvergelijkingen en numerieke methoden, waarbij studenten wordt getoond hoe ze van natuurwetten naar computationele oplossingen kunnen gaan.
Praktische Nut: De open-source MATLAB-code stelt sportwetenschappers en coaches in staat om wedstrijdstrategieën te testen en parcoursgeometrieën te analyseren zonder complexe biomechanische modellen van scratch te hoeven ontwikkelen.

Het artikel blijft bescheiden in zijn omvang, erkennend dat het een "voorlopige stap" is naar een werkelijk realistisch model. Het sluit expliciet complexe metabole processen, gedetailleerde spierbiomechanica en geavanceerde vermoeidheidsmodellen uit, en richt zich in plaats daarvan op de balans van krachten en de geometrische beperkingen van het parcours. De auteurs suggereren dat toekomstig werk "vooruitkijkende" remstrategieën en meer verfijnde pacing-modellen op basis van experimentele data zou kunnen integreren.