I. Підпірна стіна — від тиску ґрунту до коефіцієнта запасу
🧱Підпірна стіна — тиск ґрунту за Ренкіним та стійкість
Регулюйте висоту стіни, кут внутрішнього тертя ґрунту та навантаження на поверхні; спостерігайте за силою тиску, коефіцієнтами запасу та анімованим руйнуванням.
Ґрунт — це не суцільний моноліт. Насипте його купою, і він прагне осісти до свого природного кута укосу; підпірна стіна існує саме для того, щоб утримувати цю масу на місці. Перше запитання, яке ставить інженер: наскільки сильно тисне ґрунт? Теорія тиску ґрунту Ренкіна 1857 року відповідає на це одним коефіцієнтом, що виводиться з кута внутрішнього тертя ґрунту.
Активний коефіцієнт K_a зменшується зі зростанням кута тертя — пухкий, слизький ґрунт тисне сильніше, ніж щільний, гравелистий. Оскільки тиск лінійно зростає з глибиною, його інтегрування по висоті стіни дає трикутний розподіл, рівнодійна сила якого діє на третині висоти від основи. Саме ця єдина сила та точка її прикладання визначають усе інше.
K_a = tan^2(45 - phi/2) (активний коефіцієнт тиску ґрунту)
P_a = 0.5 * K_a * gamma * H^2 (рівнодійна сила тиску на метр стіни)
діє на висоті H/3 від основиОсь філософія проєктування, яка дивує новачків: ви обчислюєте не один запас міцності, а кілька, бо стіна може зруйнуватися по-справжньому різними способами. Вона може обертатися навколо носка (перекидання), може ковзати вперед цілком (ковзання) або може розчавити ґрунт під одним із країв (несуча здатність). Кожен вид руйнування має власну діаграму вільного тіла й власний коефіцієнт запасу, і їх перевіряють окремо, бо проходження одного жодним чином не гарантує проходження іншого. Стіна, достатньо важка, щоб не ковзати, усе ще може бути занадто вузькою, щоб протистояти перекиданню.
Для кожного з них прийнято різні порогові значення: як правило, коефіцієнт не менше 2,0 проти перекидання і не менше 1,5 проти ковзання. Симуляція також відстежує ексцентриситет рівнодійної на основі — якщо вона виходить за межі середньої третини, частина основи переходить у розтяг, якого не витримують ні кладка, ні ґрунт.
FS_overturning = M_resisting / M_overturning >= 2.0
FS_sliding = (mu * N) / P_a >= 1.5
eccentricity e = B/2 - M_net/N; keep e <= B/6 (середня третина)Спробуйте: Зменшіть кут тертя ґрунту з 35° до 20° і подивіться, як різко зростають K_a та сила тиску, — а потім подивіться, який коефіцієнт запасу відмовляє першим. За більшості геометрій стіна перекидається раніше, ніж ковзає, тому основу плити зазвичай розширюють з боку п'яти.
II. Гемодіаліз — як мембрана замінює нирку
🩸Гемодіаліз — кліренс, Kt/V та зниження сечовини
Проведіть сеанс діалізу: встановіть кліренс, напрямок потоку та тривалість і спостерігайте, як рівень сечовини в крові знижується до адекватного Kt/V.
Здорова нирка фільтрує весь об'єм крові багато разів на добу, видаляючи сечовину та інші азотисті відходи. Коли нирки відмовляють, цю роботу виконує діалізатор за допомогою пучка порожнистих напівпроникних волокон: кров тече всередині волокон, чиста діалізуюча рідина тече ззовні, а сечовина дифундує крізь мембрану вниз за градієнтом концентрації.
Найважливіше конструктивне рішення — пускати обидві рідини протитечією, тобто в протилежних напрямках. Якби вони текли в один бік, концентрації вирівнялися б десь на середині шляху й дифузія зупинилася б. За протитечії кров завжди зустрічає діалізат, який трохи чистіший за неї саму, тож градієнт зберігається по всій довжині волокна, а кліренс максимізується. Це той самий трюк, який еволюційно виробили зябра риб і нирки ссавців.
Клінічно видалення відходів описується кінетикою одного пулу: сечовину розглядають так, ніби вона розчинена в одному добре перемішаному об'ємі, і її концентрація спадає експоненційно в міру того, як діалізатор її очищує.
C(t) = C0 * exp(-K * t / V)
K = кліренс діалізатора (мл/хв)
t = тривалість сеансу
V = об'єм розподілу сечовини (≈ загальна вода організму)Два числа підказують клініцисту, чи був сеанс достатньо ефективним. Безрозмірний показник Kt/V об'єднує кліренс, час і розмір пацієнта в один індекс адекватності; широко застосовується цільове значення не менше 1,2 за сеанс. Коефіцієнт зниження сечовини — простіший прикроватний показник — це просто частка падіння сечовини в крові від початку до кінця сеансу. Симуляція пов'язує три регулятори — кліренс, час і об'єм — з обома показниками, тож можна побачити, чому більшому пацієнту потрібен або швидший діалізатор, або довший сеанс для досягнення тієї самої адекватності.
Kt/V = -ln(C_end/C0) (спрощена модель одного пулу); ціль >= 1,2
URR = (C0 - C_end) / C0 * 100%Спробуйте: Зафіксуйте кліренс і скоротіть сеанс — Kt/V опуститься нижче 1,2, і панель позначить сеанс як неадекватний. Тепер збільшіть об'єм розподілу (більший пацієнт) за того самого часу й подивіться, як адекватність знову падає: доза залежить від розміру тіла, а не лише від налаштувань апарата.
III. Світлова крива фотосинтезу — закон лімітувальних факторів
🌱Фотосинтез — світлова крива, C3 проти C4
Змінюйте інтенсивність світла й спостерігайте, як чиста асиміляція окреслює криву крізь точки компенсації та насичення.
Освітіть листок сильніше, і він фіксуватиме більше вуглецю — але лише до певної межі. Світлова крива має три режими. За слабкого освітлення швидкість зростає майже лінійно, обмежена суто кількістю фотонів, що надходять; її нахил — це квантовий вихід. Зі зростанням освітленості крива вигинається й насичується, обмежена вже подальшою біохімією, а не фотонами. А в самому низу, нижче точки компенсації, дихання вивільняє більше CO₂, ніж фіксує фотосинтез, тож листок стає чистим джерелом вуглекислого газу.
Симуляція моделює чисту асиміляцію як неправильну (нерівносторонню) гіперболу мінус стала складова темнового дихання, позначаючи як точку компенсації (чиста швидкість дорівнює нулю), так і точку насичення (крива вирівнюється).
A_net = (phi*I * A_max) / (phi*I + A_max) - R_d
phi = квантовий вихід (моль CO2 / моль фотонів)
I = падаючий потік фотонів
A_max = швидкість асиміляції за насичення світлом
R_d = темнове дихання
Точка компенсації: A_net = 0; Насичення: dA/dI -> 0Це наочна ілюстрація закону лімітувальних факторів Блекмана 1905 року: коли процес залежить від кількох вхідних чинників, швидкість визначає той чинник, якого найменше. На світанку обмежувальним чинником є світло; опівдні світла достатньо, і вузьким місцем стають CO₂ або температура. Посилення й без того рясного чинника нічого не дає — треба піднімати саме дефіцитний.
Пресети C3 та C4 показують, чому це має агрономічне значення. Рослини C4 — кукурудза, цукрова тростина, багато тропічних злаків — концентрують CO₂ навколо свого фіксувального ферменту, пригнічуючи марнотратне фотодихання, яке дошкуляє рослинам C3 за яскравого й спекотного клімату. Результат — вище плато насичення й краща ефективність використання води саме там, де культури C3 страждають.
Спробуйте: Встановіть пресет C3 і підвищте температуру: плато провисає, оскільки дається взнаки фотодихання. Перемкніться на C4 — і плато за високої освітленості тримається твердо: наочна ілюстрація того, чому тропічне сільське господарство спирається на культури C4.
IV. Зрошення та водний баланс ґрунту — модель відра
💧Зрошення та ґрунтова волога — баланс, MAD і дефіцит
Проведіть сезон: опади й зрошення наповнюють кореневу зону, евапотранспірація її виснажує, а тригер MAD планує наступний полив.
Для інженера зі зрошення коренева зона культури — це відро. Вода надходить — опади та зрошення; вода витікає — евапотранспірація, глибинне дренування повз корені та поверхневий стік. Зміна запасу вологи в ґрунті за будь-який проміжок часу — це просто різниця між ними. Планування зрошення — це мистецтво наповнювати відро раніше, ніж культура пересохне, не переповнюючи його й не марнуючи воду на дренаж.
delta_S = P + I - ET - D - R (водний баланс ґрунту)
P,I = опади, зрошення; D,R = дренаж, поверхневий стік
ET = ET0 * Kc (еталонна ET x коефіцієнт культури)У відра є дві важливі позначки. Польова вологоємність — це скільки води утримує ґрунт через день-два після просочення, коли надлишок уже стік під дією сили тяжіння. Точка в'янення — це рівень, нижче якого корені більше не можуть подолати утримання ґрунтом залишкової води. Різниця між ними, помножена на глибину коріння, — це загальна доступна вода, яку культура справді може використати.
З поливом не чекають до точки в'янення — культура зазнала б стресу задовго до цього. Натомість частка допустимого виснаження (MAD) встановлює комфортний тригер, зазвичай близько половини доступної води, за досягнення якого планується наступне зрошення. Сам попит культури — це зумовлена погодою еталонна евапотранспірація, масштабована коефіцієнтом культури, який зростає й спадає протягом вегетаційного сезону.
TAW = (theta_FC - theta_WP) * root_depth (загальна доступна вода)
Планувати зрошення, коли виснаження >= MAD * TAW (MAD ≈ 0,5)
Дефіцитне зрошення: подавати менше повного попиту, щоб підвищити ефективність використання водиСимуляція також досліджує дефіцитне зрошення, коли аграрій свідомо подає менше повного попиту. Урожайність падає, але часто набагато менше, ніж заощаджена вода, — тож ефективність використання води, тобто врожай на одиницю води, може зростати. У регіонах з дефіцитом води саме цей компроміс визначає всю гру.
Спробуйте: Встановіть піщаний ґрунт (низька польова вологоємність) і глибококореневу культуру, потім глинистий ґрунт, і порівняйте, як часто кожен із них ініціює зрошення. Піщані ґрунти утримують мало води й вимагають частих невеликих поливів; глинисті буферизують культуру протягом тривалих посушливих періодів.
V. Ель-Ніньйо (ENSO) — квазіперіодичне кліматичне коливання
🌊Ель-Ніньйо — осцилятор ENSO із запізненням
Нахиліть термоклин Тихого океану й спостерігайте, як океанічний індекс Ель-Ніньйо коливається нерегулярними циклами тривалістю 2–7 років.
Ель-Ніньйо — це не погодне явище, а повільна перебудова тропічної частини Тихого океану, яка хвилями поширюється на клімат усього світу. Вона виникає через тісний зв'язок між океаном і атмосферою, вперше описаний Якобом Б'єркнесом: теплі температури поверхні моря на сході послаблюють східні пасати, слабші вітри дають термокліну сплощитися, а сплощений термоклин виносить на поверхню ще тепліші води. Залишений сам на себе, цей позитивний зворотний зв'язок призвів би до неконтрольованого зростання.
Що зупиняє цей процес — і перетворює неконтрольоване зростання на коливання — так це запізнення. Повільні океанічні хвилі, викликані тими самими змінами вітру, перетинають Тихоокеанський басейн, відбиваються від його західної межі й повертаються через кілька місяців, несучи протилежний сигнал. Модель осцилятора із запізненням описує це одним диференціальним рівнянням із запізненням: локальний член зростання, запізніле від'ємне зворотне зв'язування та кубічний член, що обмежує амплітуду.
dT/dt = a*T(t) - b*T(t - tau) - eps*T(t)^3
a = зв'язане зростання за Б'єркнесом (позитивний зворотний зв'язок)
b = запізніле хвильове зворотне зв'язування (від'ємне), що надходить після затримки tau
eps = кубічне насичення, що обмежує амплітудуВзаємодія зростання, запізнення й насичення породжує самопідтримне коливання з періодом приблизно від двох до семи років. Але воно квазіперіодичне, а не годинникове: рівняння перебуває близько до хаотичних режимів, і невеликі збурення від погодного шуму зсувають час і силу кожного циклу. Саме ця чутливість пояснює, чому прогнози ENSO на понад рік залишаються складними і чому «бар'єр весняної передбачуваності» щороку засмучує сезонних прогнозистів.
На практиці події оголошують за допомогою океанічного індексу Ель-Ніньйо — ковзного середнього аномалій температури поверхні моря у визначеному регіоні Тихого океану. Якщо він достатньо довго перевищує +0,5°C — це Ель-Ніньйо; якщо опускається нижче −0,5°C — це Ла-Нінья. Симуляція нахиляє термоклин у реальному часі й відстежує проходження індексу через цей поріг.
Спробуйте: Подовжіть затримку tau, і цикл розтягнеться до верхньої межі діапазону 2–7 років; скоротіть її, і події почнуть відбуватися густо й часто. Додайте трохи шуму й подивіться, як колись регулярне коливання стає по-справжньому нерегулярним — саме тому реальний ENSO не підкоряється простому календарю.
Спільна нитка
П'ять дуже різних галузей — геотехніка, нефрологія, фізіологія рослин, агрономія та кліматична динаміка — і в кожній той самий інтелектуальний прийом: взяти фізичний механізм, записати його як рівняння, а потім прочитати рішення з результату. Коефіцієнт запасу, Kt/V, лімітувальний фактор, тригер виснаження, перетин індексом порогу. Рівняння — це легка частина; знати, за яким числом стежити, — ось у чому інженерія. Відкрийте будь-яку з п'яти симуляцій, і відповідь змінюється при кожному повороті регуляторів.