⚡ electrical engineering

A Physically Constrained Energy-Based Excess Pore Pressure Model for Saturated Sands with an Ordinary Differential Equation Interpretation

Deze studie herformuleert het energiegebaseerde model voor overtollige poriënwaterdruk in verzadigde schone zanden door het te interpreteren als een oplossing voor een Cauchy–Euler gewone differentiaalvergelijking, waardoor fysieke begrenzing wordt gegarandeerd, de kalibratie wordt teruggebracht tot één enkele parameter en een hoge nauwkeurigheid wordt bereikt over verschillende relatieve dichtheden.

Oorspronkelijke auteurs: Ferruh Gündoğan, Ayfer Erken

Gepubliceerd 2026-06-30

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Ferruh Gündoğan, Ayfer Erken

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (https://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

Het Grote Plaatje: Waarom storten zandkastelen in?

Stel je voor dat je op het strand een zandkasteel aan het bouwen bent. Als je er gewoon op staat, behoudt het zijn vorm. Maar als je er ritmisch op begint te springen, beginnen de zandkorrels te trillen. In een verzadigd zandkasteel (één die doordrenkt is met water) zorgt dit schudden ervoor dat de waterdruk tussen de korrels steeds harder gaat duwen.

Uiteindelijk wordt de waterdruk zo sterk dat de zandkorrels volledig uit elkaar worden geduwd. Het zand gedraagt zich niet langer als een vaste stof, maar begint zich te gedragen als een dikke vloeistof (zoals een milkshake). Dit wordt liquefactie (vervloeiing) genoemd, en dit is wat ervoor zorgt dat gebouwen wegzakken of kantelen tijdens aardbevingen.

Ingenieurs hebben een manier nodig om precies te voorspellen wanneer dit gebeurt. Ze gebruiken hiervoor een concept genaamd gedissipeerde energie—denk aan de totale hoeveelheid "gejakkel" of werk die het zand moet verrichten om zichzelf te herschikken. Hoe meer je het schudt, hoe meer energie het absorbeert, en hoe hoger de waterdruk wordt.

Het Oude Probleem: Een Gebroken Liniaal

Decennialang hebben ingenieurs een beroemde formule gebruikt (gemaakt door Berrill en Davis in 1985) om deze waterdruk te voorspellen. Het is als een simpele liniaal die zegt: "Hoe meer energie je erin stopt, hoe hoger de druk wordt."

Deze oude liniaal heeft echter twee grote gebreken:

Het heeft geen plafond: De formule blijft eeuwig doorgaan naar boven. In de echte wereld kan de waterdruk niet groter zijn dan het gewicht van de grond erboven (het kan niet boven de 100% uitkomen). De oude formule voorspelt soms onmogelijke getallen, zoals 150% druk, wat de natuurkunde breekt.
Het heeft verstrengelde knoppen: De formule heeft twee instellingen (parameters) die onafhankelijk van elkaar zouden moeten zijn, maar ze zitten eigenlijk aan elkaar vast. Als je de ene verandert, verandert de andere ook mee. Dit maakt het erg moeilijk om het model af te stemmen op verschillende soorten zand.

De Nieuwe Oplossing: Een Slim, Begrensd Model

De auteurs van dit artikel (Gündoğan en Erken) besloten deze liniaal te repareren. Ze keken naar gegevens uit echte laboratoriumtests op twee soorten zand: Ottawa F65 (een standaard, schoon zand dat in veel studies wordt gebruikt) en Akpınar-zand (een lokaal zand uit Istanbul).

In plaats van alleen een rechte lijn op een grafiek te tekenen, realiseerden zij zich dat de relatie tussen schudenergie en waterdruk meer lijkt op een reis in drie fasen:

Het Begin: Een langzame, zachte stijging.
Het Midden: Een snelle, steile klim.
Het Einde: Een vloeiende curve die precies afvlakt bij de limiet (100% druk).

Om dit te vangen, behandelden ze de wiskunde als een recept. Ze braken de formule af in twee delen:

De Groei-motor: Een basisterm die de snelle stijging afhandelt.
De Bemiddelaar: Een "verkeersregelaar" die de groei vertraagt naarmate de limiet wordt genaderd, zodat de 100%-plafond nooit wordt doorbroken.

Het "Synthetische" Geheim: Een Wiskundig Verhaal

Dit is het slimme deel van hun methode. Ze merkten op dat de vorm van hun nieuwe curve exact leek op de oplossingen van een specif kind type natuurkundige vergelijking genaamd een gewone differentiaalvergelijking (ODE).

Beschouw een ODE als een verhaal over hoe een systeem verandert in de loop van de tijd.

De Linkerkant van de vergelijking vertegenwoordigt het interne gedrag van het zand (hoe het weerstand biedt).
De Rechterkant vertegenwoordigt het externe schudden (de aardbeving).

Door hun model te herschrijven als deze vergelijking, konden ze strikte "verkeersregels" (fysische beperkingen) toepassen:

Regel 1: De druk moet altijd omhoog gaan (je kunt het zand niet "ont-schudden").
Regel 2: De druk moet exact stoppen met stijgen wanneer deze de 100% bereikt (liquefactie).
Regel 3: Op het moment dat de 100% wordt bereikt, moet de snelheid van de stijging nul zijn (het vlakt vloeiend af, zonder een grillige piek).

Door de wiskunde te dwingen deze regels te gehoorzamen, waren ze in staat om de "verstrengelde knoppen" te elimineren. Ze brachten het hele complexe model terug naar slechts één enkel getal (genoemd b) dat ze moeten vinden. De rest van de wiskunde wordt vastgesteld door de regels van de natuurkunde.

Wat Ze Hebben Gevonden

Ze hebben dit nieuwe "Eén-Getal-Model" getest tegen de labgegevens:

Nauwkeurigheid: Het kwam ongelooflijk goed overeen met de echte labresultaten (99% nauwkeurigheid).
Veiligheid: Het voorspelt nooit onmogelijke getallen (het blijft tussen 0% en 100%).
Eenvoud: Het is veel gemakkelijker te gebruiken dan het oude model omdat je slechts één instelling hoeft af te stemmen.

Het "Aha!"-moment: Twee Fasen van Falen

De wiskunde onthulde iets interessants over hoe zand faalt.

Vroege Fase: De eigen interne structuur van het zand (hoe dicht de korrels op elkaar zitten) bepaalt de druktoename. Dit is het "Homogene" deel van de vergelijking.
Late Fase: Zodra het zand heel dicht bij liquefactie komt (rond de 75% druk), verliest het zijn "geheugen". De interne structuur valt uiteen en de druktoename wordt bijna volledig gedreven door de externe schudenergie. Dit is het "Particuliere" deel van de vergelijking.

Samenvatting

De auteurs hebben een rommelige, onbegrensde en ingewikkelde formule voor het voorspellen van zandliquefactie genomen en deze opnieuw opgebouwd. Ze gebruikten een wiskundig "verhaal" (een vergelijking) om het model te dwingen zich te gedragen als de echte natuurkunde: vloeiend stijgen en precies stoppen bij de limiet. Het resultaat is een eenvoudiger, veiliger en nauwkeuriger instrument voor ingenieurs om te begrijpen wanneer zand tijdens een aardbeving in een vloeistof verandert.

Technische Samenvatting: Een fysiek beperkt energie-gebaseerd model voor overtollige poriëndruk in verzadigde zanden

Probleemstelling
Bodemliquefactie blijft een kritiek geotechnisch gevaar, dat traditioneel wordt beoordeeld via spanning-gebaseerde, vervormings-gebaseerde, poriëndruk-gebaseerde en energie-gebaseerde methoden. Onder deze methoden bieden Energie-gebaseerde Methoden (EBM) een duidelijk voordeel door de generatie van overtollige poriëndruk te relateren aan de cumulatieve gedissipeerde energie, waardoor ze onregelmatige seismische belasting kunnen accommoderen zonder de onzekerheid van equivalent-cyclus-conversies. De meest invloedrijke EBM-formulering is de Berrill–Davis (1985) machtswet ( $r_u = \alpha W^\beta$ ). Echter, deze studie identificeert twee fundamentele beperkingen in het Berrill–Davis-model:

Onbegrensde respons: De machtswet is wiskundig onbegrensd, waardoor de verhouding van de overtollige poriëndruk ( $r_u$ ) de eenheid kan overschrijden, wat fysiek onmogelijk is voor verzadigde zanden waar $r_u$ asymptotisch naar 1,0 zou moeten naderen bij liquefactie.
Parameterinterdependentie: De kalibratieparameters $\alpha$ en $\beta$ worden vaak als onafhankelijk behandeld, terwijl empirische studies aantonen dat zij wederzijds gecorreleerd zijn. Deze afhankelijkheid bemoeilijkt kalibratie, vermindert de fysieke interpreteerbaarheid en belemmert extrapolatie.

Methodologie
De auteurs stellen een geherformuleerd model voor uniforme, schone, verzadigde zanden voor die deze beperkingen aanpakt via een meerstaps afleidingsproces:

Data-analyse en decompositie: Gebruikmakend van cyclische triaxiale data van Ottawa F65 zand (LEAP-2017) en Akpınar zand, wordt de relatie tussen de genormaliseerde poriëndrukverhouding ( $r_u$ ) en de genormaliseerde gedissipeerde energie ( $W_N$ ) geanalyseerd op een log-log schaal. De auteurs observeren dat de data geen enkele rechte lijn volgen, maar bestaan uit ten minste drie lineaire segmenten.
Polynomiale benadering: De data worden benaderd door een zesde-orde polynoom zonder constante term. Deze functie wordt gedecomposeerd in een primaire groeiterm ( $W_N^b$ ) en een mediatiefunctie $g(W_N)$ .
Taylor-expansie en vereenvoudiging: De mediatiefunctie wordt verder geanalyseerd. Een Taylor-expansie van het exponentiële component onthult dat een derde-orde polynoom de datatrends adequaat representeert. Dit leidt tot een kandidaat-modelvorm: $r_u = W_N^b (a_0 + a_1 W_N + a_2 W_N^2 + a_3 W_N^3)$ .
Interpretatie via gewone differentiaalvergelijking (ODE): Door de wiskundige structuur te herkennen, wordt het model geherinterpreteerd als de oplossing van een niet-homogene Cauchy–Euler gewone differentiaalvergelijking. De linkerzijde vertegenwoordigt de interne respons van het systeem, terwijl de rechterzijde de externe belastingscondities vertegenwoordigt.
Fysieke beperkingen: Om onfysische fluctuaties en parameterinterdependentie te elimineren, wordt het model onderworpen aan specifieke randvoorwaarden bij het liquefactiepunt ( $W_N = 1, r_u = 1$ $W_{N} = 1, r_{u} = 1$ ):
- $r_u = 1$ bij $W_N = 1$ .
- Een horizontale raaklijn ( $r_u' = 0$ ) bij $W_N = 1$ , wat ervoor zorgt dat de curve afvlakt naarmate de limiet wordt bereikt.
- Monotone toename ( $r_u' > 0$ ) gedurende het gehele belastingsproces.
- Een specifieke relatie tussen coëfficiënten ( $a_2 = -3a_3$ ) afgeleid van het verdwijnen van de niet-homogene term bij liquefactie.

Deze beperkingen bepalen analytisch de coëfficiënten van de mediapolynoom in termen van de enkele parameter $b$ , waardoor het model wordt gereduceerd tot een vorm die uitsluitend wordt beheerst door $b$ .

Kernbijdragen

Fysiek beperkte formulering: Het voorgestelde model is geconstrueerd om inherent te voldoen aan $0 \le r_u \le 1$ , wat de onfysische overschatting van poriëndruk voorkomt die gebruikelijk is bij het Berrill–Davis-model.
Enkele kalibratieparameter: Door de mediacoëfficiënten analytisch te koppelen aan de groeiexponent $b$ , elimineert het model de noodzaak voor de onafhankelijke kalibratie van meerdere parameters, wat de interdependentie oplost die in eerdere formuleringen werd gevonden.
ODE-framework: De studie biedt een nieuwe interpretatie van de empirische poriëndrukrelatie als de oplossing van een niet-homogene Cauchy–Euler ODE. Dit framework scheidt de respons in een homogeen deel (dat de interne materiaaleigenschappen beheerst, dominant in de vroege stadia) en een particulier deel (dat de externe belastingscondities beheerst).
Validatie: Het model wordt gevalideerd tegen hoogwaardige cyclische triaxiale data voor Ottawa F65 en Akpınar zanden over een breed scala aan relatieve dichtheden (van los tot zeer dicht), waarbij een $R^2 \ge 0,99$ wordt behaald.

Resultaten

Het vereenvoudigde model (Eq. 20) reproduceert succesvol de gemeten poriëndrukgeschiedenissen voor zowel Ottawa F65 als Akpınar zanden.
Er is vastgesteld dat de parameter $b$ varieert met de relatieve dichtheid (bijv. $b \approx 0,36$ voor 51% relatieve dichtheid, $b \approx 0,43$ voor 80–91% relatieve dichtheid).
Het model legt de niet-lineaire evolutie van de poriëndruk vast, inclusief de initiële snelle stijging en de daaropvolgende afvlakking naarmate liquefactie wordt benaderd.
Een vergelijking met het oorspronkelijke Berrill–Davis-model laat zien dat de voorgestelde formulering een consistentere beschrijving biedt van de poriëndrukontwikkeling, met name nabij de eindconditie van volledige liquefactie, waar het oorspronkelijke model significante spreiding en onbegrensde groei vertoont.
De toepassing van het vereenvoudigde model op de resterende LEAP-2017 dataset bevestigt het vermogen om poriëndrukgeneratie over verschillende relatieve dichtheden en cyclische spanningverhoudingen te simuleren, mits $b \ge 0,27$ om lokale maxima te vermijden.

Significantie
Het artikel beweert dat dit werk een levensvatbaar, compact en fysiek consistent alternatief biedt voor de klassieke Berrill–Davis machtswet voor liquefactiestudies waarbij uniforme, schone zanden betrokken zijn. Door de respons te begrenzen en parameterinterdependentie te verwijderen via een ODE-gebaseerde afleiding, vereenvoudigt het model de kalibratie terwijl het een hoge nauwkeurigheid behoudt. De auteurs benadrukken dat hoewel de formulering empirisch is, de afleiding ervan uit regressie onder strikte fysieke beperkingen het in staat stelt het onderliggende mechanisme van liquefactie te vangen: een initiële fase die wordt beheerst door interne granulaire kenmerken (homogene oplossing) gevolgd door een fase waarin het materiaal zijn mechanische geheugen verliest en aanzienlijke energie vereist om volledige liquefactie te bereiken (particuliere oplossing). De reikwijdte is expliciet beperkt tot uniforme, schone, verzadigde zanden.