Головна Фізика та Механіка Сипкі матеріали — Піщаний насип та кут відкосу

🏜 Сипкі матеріали — Піщаний насип та кут відкосу

Симуляція сипких матеріалів: сипте зерна піску, спостерігайте зростання насипу, вимірюйте кут відкосу, відстежуйте лавинну динаміку. Клітинний автомат BTW.

Фізика та Механіка2DСередній60 FPS
granular-materials ↗ Відкрити окремо
DRAG · SCROLL · CLICK — керуйте прямо у вікні симуляції.

Фізика та рівняння

Клітинний автомат Бака-Танга-Візенфельда (BTW): якщо z[i] − z[i±1] > zc, то z[i] −= 1, z[i±1] += 1 (каскад перекидання → лавина). Кут відкосу θ = arctan(zc/діаметр_зерна). Значення zc регулюється повзунком від 1 до 8, а темп додавання — від 1 до 20 зерен за кадр.

Кут відкосу

Кут відкосу — найкрутіший кут, за якого насип сипкого матеріалу залишається стійким, зазвичай 25–40° для сухих матеріалів. Його визначає насамперед міжчастинкове тертя: кутасті частинки дають вищий кут через зчеплення, а волога через капілярні містки дозволяє ще крутіші схили.

Самоорганізована критичність

Насип природно еволюціонує до критичного схилу, за якого одне нове зерно може не спричинити нічого або запустити лавину, що охоплює всю систему. Розміри лавин підпорядковуються степеневому закону — та сама статистична ознака, яку спостерігають у землетрусах і лісових пожежах.

Про цю симуляцію

Сипкі матеріали — пісок, гравій, рис, фармацевтичні порошки — це скупчення дискретних частинок, що взаємодіють через тертя й контактні сили та поводяться так, як не поводиться жоден звичайний стан речовини: вони можуть текти, як рідини, витримувати навантаження, як тверді тіла, і стискатися, як гази. Ця симуляція сипле зерна на насип і відтворює отриману колективну динаміку за допомогою клітинного автомата піщаної купи Бака-Танга-Візенфельда (BTW), який лежить в основі проєктування силосів, бункерів і хвостосховищ.

🔬 Що це показує

Кожна колонка сітки зберігає висоту; зерна додаються по одному. Якщо схил між сусідніми колонками перевищує критичне значення zc, зерно перекидається вниз, що може спричинити каскад перекидань сусідів — лавину. Зерна, що досягають країв, випадають, тож насип самоорганізується до кута відкосу, балансуючи на межі стійкості.

🎮 Як користуватися

Оберіть режим — Конус (сипати в центр), Потік (сипати з лівого краю) або Дощ (сипати випадково) — і регулюйте Критичний схил zc та Темп додавання повзунками. Кнопка «Додати зерна» впорскує порцію, «Очистити» скидає насип до нуля, а клік по полотну сипле зерна прямо в колонку під курсором.

💡 Чи знали ви?

Розміри лавин у піщаному насипі підпорядковуються степеневому закону — типового розміру не існує. Одне додане зерно може не зробити нічого, а може спричинити обвал у масштабах усієї системи, тому ця модель стала засадничим прикладом самоорганізованої критичності — тієї самої статистики, яку бачать у землетрусах і лісових пожежах.

Поширені запитання

Чому сипкі матеріали поводяться інакше, ніж рідини та тверді тіла?

Рідини й тверді тіла перебувають у термодинамічній рівновазі — теплові флуктуації безперервно досліджують доступні мікростани. Зерна сипких матеріалів настільки великі, що теплова енергія kT нехтовно мала порівняно з гравітаційною потенційною енергією одного зерна — це співвідношення становить приблизно 10⁻¹² для піщинок 1 мм. Зерна не перебудовуються самовільно; вони залишаються в механічно стійкій конфігурації, доки не буде прикладено зовнішньої енергії (струшування, вібрація, зсув). Через це сипкі матеріали залежать від передісторії: той самий матеріал може бути щільним чи пухким, міцним чи слабким, залежно від того, як він був підготовлений.

Що таке кут відкосу і від чого він залежить?

Кут відкосу — це найкрутіший кут, за якого насип сипкого матеріалу залишається стійким і не сповзає, зазвичай 25–40° для сухих сипких матеріалів. Він визначається насамперед міжчастинковим тертям: вище тертя дозволяє крутіші схили. Форма частинок теж має значення — кутасті частинки мають вищий кут відкосу, ніж кулясті, оскільки їхнє зчеплення підвищує ефективне тертя. Волога створює капілярні містки між частинками, різко підвищуючи зчеплення. Динамічний кут відкосу (під час течії) трохи нижчий за статичний кут відкосу (у стані спокою), що спричиняє характерний цикл лавин.

Що таке ефект бразильського горіха?

Ефект бразильського горіха описує парадоксальне спостереження: коли суміш частинок різного розміру струшують вертикально, великі частинки піднімаються нагору, а дрібні опускаються вниз. Цей контрінтуїтивний результат зумовлений конвекцією: струшування створює висхідний конвективний потік у центрі й низхідний потік біля стінок. Великі частинки піднімаються конвективним потоком у центрі, але заширокі, щоб опуститися вузьким пристінним простором, і накопичуються нагорі. Просіювання теж робить свій внесок: коли шар сипкого матеріалу розширюється, дрібні частинки провалюються в проміжки під великими.

Що таке ефект Янсена і чому він важливий для конструювання силосів?

У високому силосі, заповненому сипким матеріалом, тиск біля основи не зростає безмежно з висотою заповнення (як у рідині), а насичується вже після глибини приблизно в один діаметр силосу. Цей ефект Янсена виникає через те, що тертя об стінки силосу перенаправляє частину ваги горизонтально на стінки, а не вниз. Конструктори силосів повинні враховувати насичення Янсена, щоб правильно розрахувати тиски на основу та стінки й уникнути катастрофічних обвалів силосу.

Як метод дискретних елементів (DEM) моделює сипкі матеріали?

DEM (розроблений Кунделлом і Страком 1979 року) подає кожну частинку як сферу і відстежує її положення, швидкість та орієнтацію. На кожному кроці за часом контактні сили обчислюються за моделями пружина-демпфер або Герца для нормальної сили та за кулонівським тертям для дотичної сили, а рух інтегрується за другим законом Ньютона (F = ma). DEM є стандартним інструментом для промислового моделювання течії сипких матеріалів, обробки фармацевтичних порошків і геомеханічного аналізу стійкості схилів.

Що контролює повзунок Критичний схил z_c у цій моделі піщаного насипу?

z_c (від 1 до 8) — це різниця висот, яку колонка піщаного насипу може витримати відносно сусідів, перш ніж перекинутися в цій моделі клітинного автомата, адаптованій із класичної піщаної купи Бака-Танга-Візенфельда. Щойно висота колонки перевищує сусідню на z_c, зерна перерозподіляються на цього сусіда, що може, своєю чергою, підштовхнути його за власний поріг — нижче значення z_c змушує насип перекидатися й лавинити легше, спричиняючи частіші, але зазвичай менші за розміром каскади.

Чому розміри лавин у цій симуляції підпорядковуються степеневому закону, а не групуються навколо одного типового розміру?

Це ознака самоорганізованої критичності: коли зерна додаються по одному, насип природно еволюціонує до критичного схилу, за якого одне нове зерно може не спричинити жодного обвалу або запустити лавину, що охоплює всю систему. Оскільки в правило перекидання не закладено характерного розміру лавини, розподіл розмірів лавин підпорядковується степеневому закону — та сама статистична ознака, яку спостерігають у реальних піщаних насипах, землетрусах і лісових пожежах.

Що змінює повзунок Темп додавання, і чи призводить швидше додавання зерен до більших лавин?

Темп додавання (від 1 до 20 зерен за кадр) визначає, скільки нових зерен додається до насипу щокадру анімації. Він прискорює заповнення насипу й досягнення ним критичного стану, але не змінює лежачу в основі степеневу статистику окремих лавин, оскільки кожна подія перекидання все одно розв'язується по одній взаємодії зерна з колонкою за тим самим правилом порогу z_c.

Що відбувається, коли я натискаю Очистити, і чому насип поводиться інакше одразу після цього?

Очистити скидає висоту кожної колонки сітки назад до нуля, видаляючи всі накопичені зерна. Одразу після цього насип перебуває далеко нижче критичного схилу всюди, тому перші додані зерна просто накопичуються, не викликаючи жодних лавин — характерна переривчаста, степенева поведінка лавин повертається лише тоді, коли достатньо зерен накопичиться, щоб висоти колонок знову наблизилися до порогу z_c.

Чи є ця модель піщаного насипу тим самим, що й фізична поведінка кута відкосу, описана вище?

Вони пов'язані, але відрізняються: кут відкосу описує максимальний стійкий кут схилу реального 3D-насипу сипкого матеріалу під дією гравітації й тертя, тоді як ця симуляція є спрощеною абстракцією клітинного автомата (піщана купа Бака-Танга-Візенфельда), де перекидання керується виключно цілочисельним правилом різниці висот z_c, а не безперервними фізичними силами. Клітинна модель відтворює ту саму якісну поведінку самоорганізованої критичності — періодичні лавини різних розмірів — не моделюючи фактичне тертя, форму й вологість, які визначають кут відкосу реального насипу.

Схожі симуляції