Вікіпедія

Метод матеріальної точки

Метод матеріальної точки (ММТ, англ. material point method, MPM) — це числовий метод, який використовується для моделювання поведінки твердих тіл, рідин, газів та будь-якого іншого суцільного матеріалу. Зокрема, це стійкий метод дискретизації простору для моделювання багатофазових взаємодій (тверде тіло-рідина-газ). У ММТ неперервне тіло описується кількома невеликими елементами Лагранжа, які називаються «матеріальними точками». Ці матеріальні точки оточені фоновою сіткою, яка використовується для обчислення різних характеристик, зокрема, таких як градієнт деформації. На відміну від інших методів на основі сітки, таких як метод скінченних елементів, метод скінченних об'ємів або скінченних різниць, ММТ не є методом на основі сітки, а натомість класифікується як безсітковий або неперервний метод на основі частинок, прикладами якого є гідродинаміка згладжених частинок та перидинаміка. Попри наявність фонової сітки, ММТ не стикається з недоліками методів на основі сітки (висока деформація заплутування, похибки адвекції тощо), що робить його перспективним і потужним інструментом у обчислювальній механіці.

ММТ був запропонований на початку 1990-х років професорами Деборою Л. Сульскі, Чженем Ченом та Говардом Л. Шрайером із університету Нью-Мексико, як розширення подібного методу, відомого як FLIP (подальше розширення методу частинок в комірці[en], англ. Particle-in-cell, PIC) для обчислювальної динаміки твердого тіла. Після цієї початкової розробки ММТ вже отримав подальший розвиток як у національних лабораторіях, як в університеті Нью-Мексико, так і в університеті штату Орегон, університеті Юти та інших закладах США та світу. Останнім часом кількість інституцій, які досліджують ММТ, зростає завдяки популярності та обізнаності з різних джерел, зокрема, таким як використання ММТ у діснеївському фільмі Крижане серце.

Алгоритм ред.

Моделювання ММТ складається з наступних етапів:

(До фази інтеграції часу)

  1. Ініціалізація сітки та матеріальних точок.
    1. Геометрія розбивається на сукупність матеріальних точок, кожна зі своїми властивостями матеріалу та початковими умовами (швидкість, напруження, температура тощо).
    2. Сітка, яка використовується лише для надання місця для розрахунків градієнта, зазвичай створюється для покриття достатньо великої області, щоб заповнити очікуваний обсяг обчислювальної області, необхідної для моделювання.

(Під час фази інтеграції часу — явне формулювання)

  1. Величини матеріальних точок екстраполюються на вузли сітки.
    1. Маса матеріальної точки ( ), імпульс ( ), напруження ( ), і зовнішні сили ( ) екстраполюються на вузли в кутах комірок, всередині яких знаходяться матеріальні точки. Найчастіше це робиться за допомогою стандартних функцій лінійної форми ( ), які так само використовуються в МСЕ.
    2. Сітка використовує значення матеріальної точки для створення мас ( ), швидкості ( ), внутрішніх та зовнішніх векторів сил ( ,  ) для вузлів:
       
       
       
       
  2. На сітці розв'язуються рівняння руху.
    1. 2-й закон Ньютона розв'язується, щоб отримати прискорення у вузлі ( )
       
    2. У вузлах обчислюються нові швидкості ( ).
       
  3. Похідні терміни екстраполюються назад до матеріальних точок
    1. Прискорення матеріальної точки ( ), градієнт деформації ( ) (або швидкість деформації ( ) залежно від того, яка теорія деформацій використовується) екстраполюється з навколишніх вузлів за допомогою функцій, аналогічним тим, що раніше були застосовані до ( ).
       
       
    2. Змінні в матеріальних точках: положення, швидкості, деформації, напруження тощо, потім оновлюються цими швидкостями залежно від вибраної схеми інтегрування та відповідної моделі побудови.
  4. Скидання сітки.
    Тепер, коли матеріальні точки повністю оновлено на наступному часовому кроці, сітка скидається, щоб дозволити почати наступний часовий крок.

Історія PIC та ММТ ред.

PIC спочатку був задуманий для розв'язання задач динаміки рідини і розроблений Харлоу[en] в Лос-Аламосській національній лабораторії в 1957 році[1]. Одним із перших кодів PIC була програма Fluid-Implicit Particle (FLIP), яка була створена Брекбілом у 1986 році[2] і відтоді постійно розвивалася. До 1990-х років метод PIC використовувався в основному в гідродинаміці.

Мотивовані потребою кращого моделювання проблем проникнення в твердій динаміці, Сулскі, Чен і Шрайер почали в 1993 році переформулювати PIC і розробляти ММТ при фінансовій підтримці Sandia National Laboratories[3]. Початковий ММТ був потім додатково розширений Барденхагеном та іншими, щоб додати у модель фрикційний контакт[4], що дозволило моделювати зернистий потік[5], і за Нейрном включити явні тріщини[6] та поширення тріщини (відоме як CRAMP).

Нещодавно реалізація ММТ, заснована на мікрополярному континуумі Коссера[7], була використана для моделювання зернистого потоку з високим зсувом, такого як скидання силосу. Використання ММТ було розширено в геотехнічній інженерії завдяки нещодавній розробці квазістатичного, неявного вирішувача ММТ, який забезпечує чисельно стабільний аналіз задач великої деформації в механіці ґрунтів[8].

Щорічні семінари з використання ММТ проводяться в різних місцях у Сполучених Штатах. П'ятий семінар ММТ відбувся в Університеті штату Орегон, в Корваллісі, штат Орегон, 2 і 3 квітня 2009 року.

Застосування PIC та ММТ ред.

Використання методу PIC або ММТ можна розділити на дві великі категорії: по-перше, існує багато застосувань, що включають динаміку рідини, фізику плазми, магнітогідродинаміку та багатофазні застосунки. Друга категорія застосунків включає задачі з механіки твердого тіла.

Динаміка рідини та багатофазне моделювання ред.

Метод PIC був використаний для моделювання широкого спектра взаємодій рідини і твердого тіла, включаючи динаміку морського льоду[9], проникнення в біологічні м'які тканини[10], фрагментацію наповнених газом балонів[11], розсіювання атмосферних забруднювачів[12], багатомасштабне моделювання, що поєднує молекулярну динаміку з ММТ[13][14] та взаємодією рідина-мембрана[15]. Крім того, код FLIP на основі PIC був застосований в інструментах магнітогідродинаміки та обробки плазми, а також для моделювання в астрофізиці та течії по вільній поверхні[16].

В результаті спільних зусиль математичного факультету Каліфорнійського університету в Лос-Анджелесі та студії Walt Disney Animation Studios ММТ було успішно використано для імітації снігу в комп'ютерному анімаційному фільмі Крижане серце 2013 року[17][18][19].

Механіка твердого тіла ред.

ММТ також широко використовується в механіці твердого тіла для імітації удару, проникнення, зіткнення та відскоку, а також поширення тріщини[20][21]. ММТ також став широко використовуваним методом у галузі механіки ґрунтів: його використовували для моделювання зернистого потоку, випробування на швидкість чутливих глин[22], зсувів[23][24][25], силосного розвантаження, забивання паль, випробування на падіння конуса[23][26][27][28], наповнення ковша та руйнування матеріалу; а також моделювання розподілу напружень у ґрунті[29], ущільнення та зміцнення. Зараз він використовується в проблемах механіки деревини, таких як моделювання поперечного стиснення на рівні клітин, включаючи контакт клітинної стінки[30]. Ця робота також отримала нагороду Джорджа Марра як наукова стаття року від Товариства деревознавства та технологій[31].

Класифікація кодів PIC та ММТ ред.

ММТ в контексті чисельних методів ред.

Однією з підмножин числових методів є безсіткові методи[en], які визначаються як методи, для яких «попередньо визначена сітка не потрібна, принаймні в інтерполяції змінних поля». В ідеалі такий метод не використовує сітку «у процесі розв'язання задачі, що обумовлена диференціальними рівняннями в часткових похідних, у заданій довільній області, з урахуванням усіх видів граничних умов», хоча наявні методи не є ідеальними та не вдаються до хоча б один із цих підходів. Безсіткові методи, які також іноді називають методами частинок, мають «загальну особливість, що історія змінних, які описують стан, простежується в точках (частинках), які не пов'язані з жодним елементом сітки, спотворення якої є джерелом чисельних ускладнень». Як видно з цих різних інтерпретацій, деякі вчені вважають ММТ безсітковим методом, а інші — ні. Проте всі погоджуються, що ММТ є методом частинок.

Довільні методи Лагранжа-Ейлера (Arbitrary Lagrangian Eulerian, ALE) утворюють ще одну підмножину числових методів, яка включає ММТ. Чисто лагранжеві[en] методи використовують структуру, в якій простір розбивається на початкові підоб'єми, шляхи потоків яких потім відображаються в часі. Чисто ейлерівські[en] методи, з іншого боку, використовують структуру, в якій описується рух матеріалу відносно сітки, яка залишається нерухомою в просторі протягом усього розрахунку. Як видно з назви, методи ALE поєднують лагранжеву та ейлерову системи відліку.

Підкласифікація ММТ та PIC ред.

Методи PIC можуть бути засновані або на колокації сильної форми, або на дискретизації слабкої форми основного диференціального рівняння в частинних похідних (ДРЧП). Методи, засновані на сильній формі, правильно називаються PIC методами скінченного об'єму. Ті, що ґрунтуються на слабкій дискретизації форми ДРЧП, можна назвати PIC або ММТ.

Розв'язувачі ММТ можуть моделювати задачі в одному, двох або трьох просторових вимірах, а також можуть моделювати осесиметричні задачі. ММТ може бути реалізований для розв'язування квазістатичних або динамічних рівнянь руху залежно від типу задачі, яку потрібно моделювати. Декілька версій ММТ включають метод узагальненої інтерполяційної точки[32]; метод інтерполяції в області конвекційних частинок[33]; метод інтерполяції найменших квадратів конвекційних частинок[34].

Інтеграція за часом, що використовується для ММТ, може бути як явною[en], так і неявною. Перевагою неявної інтеграції є гарантована стабільність навіть для великих часових кроків. З іншого боку, явна інтеграція виконується набагато швидше і простіше у реалізації.

Переваги ред.

Порівняння з МСЕ ред.

На відміну від методу скінченних елементів (МСЕ), ММТ не вимагає періодичних кроків перемішування та перетворення змінних стану, і тому краще підходить для моделювання великих деформацій матеріалу. У ММТ усю інформацію про стан обчислення зберігають частинки, а не точки сітки. Таким чином, коли сітка повертається у вихідне положення після кожного циклу обчислення, числова помилка не виникає, і не потрібен алгоритм повторної сітки.

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

Порівняння з чистими методами частинок ред.

Оскільки в ММТ вузли залишаються фіксованими на звичайній сітці, обчислення градієнтів є тривіальним.

У моделюванні з двома або більше фазами досить легко виявити контакт між сутностями, оскільки частинки можуть взаємодіяти через сітку з іншими частинками в тому самому тілі, з іншими твердими тілами та з рідинами.

Недоліки ММТ ред.

ММТ є дорожчим з точки зору зберігання даних, ніж інші методи, оскільки ММТ використовує сітку, а також дані частинок. ММТ є дорожчим з точки зору обчислень, ніж МСЕ, оскільки сітку потрібно скинути в кінці кожного кроку обчислення ММТ та повторно ініціалізувати на початку наступного кроку. У ММТ можуть виникнути помилкові коливання, коли частинки перетинають межі сітки в ММТ, хоча цей ефект можна мінімізувати за допомогою узагальнених методів інтерполяції (GIMP). У ММТ, як і в МСЕ, розмір і орієнтація сітки можуть впливати на результати обчислень: наприклад, у ММТ, як відомо, місце надриву особливо чутливе до точності сітки. Однією з проблем стабільності в ММТ, яка не виникає в МСЕ, є помилки перетину комірок і помилки порожнього простору[35], оскільки кількість точок інтегрування (матеріальних точок) не залишається постійною в комірці.

Примітки ред.

  1. ↑ Johnson, N. L. (1996). The legacy and future of CFD at Los Alamos. Proceedings of the 1996 Canadian CFD Conference (англ.). OSTI 244662.
  2. ↑ Brackbill, J. U.; Ruppel, H. M. (1986). FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions. Journal of Computational Physics. 65 (2): 314—343. Bibcode:1986JCoPh..65..314B. doi:10.1016/0021-9991(86)90211-1. ISSN 0021-9991.
  3. ↑ Sulsky, D.; Chen, Z.; Schreyer, H. L. (1994). A particle method for history-dependent materials. Computer Methods in Applied Mechanics and Engineering. 118 (1): 179—196. doi:10.1016/0045-7825(94)90112-0. ISSN 0045-7825. Архів оригіналу за 6 лютого 2022. Процитовано 6 квітня 2022.
  4. ↑ Bardenhagen, S. G.; Brackbill, J. U.; Sulsky, D. L. (1998). Shear deformation in granular materials (англ.). doi:10.2172/329539. OSTI 329539.
  5. ↑ Więckowski, Zdzisław; Youn, Sung-Kie; Yeon, Jeoung-Heum (1999). A particle-in-cell solution to the silo discharging problem. International Journal for Numerical Methods in Engineering (англ.). 45 (9): 1203—1225. Bibcode:1999IJNME..45.1203W. doi:10.1002/(SICI)1097-0207(19990730)45:9<1203::AID-NME626>3.0.CO;2-C. ISSN 1097-0207.
  6. ↑ Nairn, J. A. (2003). Material Point Method Calculations with Explicit Cracks. Computer Modeling in Engineering & Sciences. 4 (6): 649—664. doi:10.3970/cmes.2003.004.649.
  7. ↑ Coetzee, Corne J. (2004). The modelling of granular flow using the particle-in-cell method (Дипломна робота PhD) (англ.). Stellenbosch : University of Stellenbosch. Архів оригіналу за 18 червня 2019. Процитовано 29 червня 2022.
  8. ↑ Beuth, L., Coetzee, C.J., Bonnier, P. and van den Berg, P. «Formulation and validation of a quasi-static material point method.» In 10th International Symposium on Numerical Methods in Geomechanics, 2007.
  9. ↑ Wang, R.-X; Ji, S.-Y.; Shen, Hung Tao; Yue, Q.-J. (2005). Modified PIC method for sea ice dynamics. China Ocean Engineering. 19: 457—468 — через ResearchGate.
  10. ↑ Ionescu, I., Guilkey, J., Berzins, M., Kirby, R., and Weiss, J. «Computational simulation of penetrating trauma in biological soft tissues using MPM [Архівовано 15 квітня 2021 у Wayback Machine.].»
  11. ↑ Banerjee, Biswajit (2012). Material point method simulations of fragmenting cylinders. ResearchGate (англ.). arXiv:1201.2439. Bibcode:2012arXiv1201.2439B. Процитовано 18 червня 2019.
  12. ↑ Patankar, N. A.; Joseph, D. D. (2001). Lagrangian numerical simulation of particulate flows. International Journal of Multiphase Flow. 27 (10): 1685—1706. doi:10.1016/S0301-9322(01)00025-8. ISSN 0301-9322.
  13. ↑ Lu, H.; Daphalapurkar, N. P.; Wang, B.; Roy, S.; Komanduri, R. (2006). Multiscale simulation from atomistic to continuum – coupling molecular dynamics (MD) with the material point method (MPM). Philosophical Magazine. 86 (20): 2971—2994. Bibcode:2006PMag...86.2971L. doi:10.1080/14786430600625578. ISSN 1478-6435. S2CID 137383632.
  14. ↑ Ma, Jin (2006). Multiscale Simulation Using the Generalized Interpolation Material Point Method, Discrete Dislocations and Molecular Dynamics (Дипломна робота PhD). Oklahoma State University. Архів оригіналу за 19 вересня 2015. Процитовано 6 квітня 2022.
  15. ↑ York, Allen R.; Sulsky, Deborah; Schreyer, Howard L. (2000). Fluid–membrane interaction based on the material point method. International Journal for Numerical Methods in Engineering (англ.). 48 (6): 901—924. Bibcode:2000IJNME..48..901Y. doi:10.1002/(SICI)1097-0207(20000630)48:6<901::AID-NME910>3.0.CO;2-T. ISSN 1097-0207.
  16. ↑ Liu, Wing Kam; Li, Shaofan (2002). Meshfree and particle methods and their applications. Applied Mechanics Reviews. 55 (1): 1—34. Bibcode:2002ApMRv..55....1L. doi:10.1115/1.1431547. ISSN 0003-6900.
  17. ↑ Marquez, Letisia (27 лютого 2014). UCLA's mathematicians bring snow to life for Disney's "Frozen". UCLA Today. Архів оригіналу за 10 March 2014. Процитовано 6 березня 2014.
  18. ↑ Alexey Stomakhin; Craig Schroeder; Lawrence Chai; Joseph Teran; Andrew Selle (August 2013). A material point method for snow simulation (PDF). Walt Disney Animation Studios. Архів оригіналу (PDF) за 24 March 2014. Процитовано 6 березня 2014.
  19. ↑ Making of Disney's Frozen: A Material Point Method For Snow Simulation. CG Meetup. 21 листопада 2013. Архів оригіналу за 7 січня 2014. Процитовано 18 січня 2014.
  20. ↑ Karuppiah, Venkatesh (2004). Implementation of irregular mesh in MPM for simulation of mixed mode crack opening in tension (Дипломна робота Master's). Oklahoma State University.
  21. ↑ Daphalapurkar, Nitin P.; Lu, Hongbing; Coker, Demir; Komanduri, Ranga (1 січня 2007). Simulation of dynamic crack growth using the generalized interpolation material point (GIMP) method. International Journal of Fracture (англ.). 143 (1): 79—102. doi:10.1007/s10704-007-9051-z. ISSN 1573-2673. S2CID 20013793.
  22. ↑ Tran, Quoc-Anh; Solowski, Wojciech; Thakur, Vikas; Karstunen, Minna (2017). Modelling of the Quickness Test of Sensitive Clays Using the Generalized Interpolation Material Point Method. Landslides in Sensitive Clays. Advances in Natural and Technological Hazards Research. 46: 323—326. doi:10.1007/978-3-319-56487-6_29. ISBN 978-3-319-56486-9. Архів оригіналу за 6 квітня 2022. Процитовано 6 квітня 2022.
  23. ↑ а б Tran, Quoc-Anh; Solowski, Wojciech (2019). Generalized Interpolation Material Point Method modelling of large deformation problems including strain-rate effects – Application to penetration and progressive failure problems. Computers and Geotechnics. 106 (1): 249—265. doi:10.1016/j.compgeo.2018.10.020. Архів оригіналу за 6 квітня 2022. Процитовано 6 квітня 2022.
  24. ↑ Llano-Serna, Marcelo A.; Farias, Márcio M.; Pedroso, Dorival M. (2016). An assessment of the material point method for modelling large scale run-out processes in landslides. Landslides (англ.). 13 (5): 1057—1066. doi:10.1007/s10346-015-0664-4. ISSN 1612-510X.
  25. ↑ Llano Serna, Marcelo Alejandro; Muniz-de Farias, Márcio; Martínez-Carvajal, Hernán Eduardo (21 грудня 2015). Numerical modelling of Alto Verde landslide using the material point method. DYNA. 82 (194): 150—159. doi:10.15446/dyna.v82n194.48179. ISSN 2346-2183. Архів оригіналу за 18 серпня 2020. Процитовано 6 квітня 2022.
  26. ↑ Tran, Quoc-Anh; Solowski, Wojciech; Karstunen, Minna; Korkiala-Tanttua, Leena (2017). Modelling of Fall-cone Tests with Strain-rate Effects. Procedia Engineering. 175: 293—301. doi:10.1016/j.proeng.2017.01.029. Архів оригіналу за 6 квітня 2022. Процитовано 6 квітня 2022.
  27. ↑ Llano-Serna, M.A.; Farias, M.M.; Pedroso, D.M.; Williams, David J.; Sheng, D. (2016). Simulations of Fall Cone Test in Soil Mechanics Using the Material Point Method. Applied Mechanics and Materials. 846: 336—341. doi:10.4028/www.scientific.net/AMM.846.336. ISSN 1662-7482.
  28. ↑ Llano-Serna, M; Farias, M (3 червня 2014), Hicks, Michael; Brinkgreve, Ronald; Rohe, Alexander (ред.), Use of generalized material point method (GIMP) to simulate shallow wedge penetration, Numerical Methods in Geotechnical Engineering (англ.), CRC Press: 259—264, doi:10.1201/b17017-48, ISBN 9781138001466
  29. ↑ Llano-Serna, M.A.; Farias, M.M. (2016). Validación numérica, teórica y experimental del método del punto material para resolver problemas geotécnicos. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería (ісп.). 32 (2): 110—115. doi:10.1016/j.rimni.2015.02.008.
  30. ↑ Nairn, John A. (2007). Numerical Simulations of Transverse Compression and Densification in Wood. Wood and Fiber Science (амер.). 38 (4): 576—591. ISSN 0735-6161. Архів оригіналу за 5 липня 2020. Процитовано 6 квітня 2022.
  31. ↑ Society of Wood Science and Technology: George Marra Award Recipients. 2007. Архів оригіналу за 23 вересня 2007. Процитовано 18 червня 2019.
  32. ↑ Bardenhagen, S. G.; Kober, E. M. (2004). The Generalized Interpolation Material Point Method. Computer Modeling in Engineering & Sciences. 5: 477—496. doi:10.3970/cmes.2004.005.477. Архів оригіналу за 19 січня 2022. Процитовано 6 квітня 2022.
  33. ↑ Sadeghirad, A.; Brannon, R. M.; Burghardt, J. (2011). A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations. International Journals for Numerical Methods in Engineering. 86 (12): 1435—1456. doi:10.1002/nme.3110. Архів оригіналу за 19 січня 2022. Процитовано 6 квітня 2022.
  34. ↑ Tran, Quoc-Anh; Solowski, Wojciech; Berzins, Martin; Gulkey, James (2020). A convected particle least square interpolation material point method. International Journals for Numerical Methods in Engineering. 121 (6): 1068—1100. doi:10.1002/nme.6257. Архів оригіналу за 6 квітня 2022. Процитовано 6 квітня 2022.
  35. ↑ Tran, Quoc-Anh; Solowski, Wojciech (2017). Temporal and null‐space filter for the material point method. International Journal for Numerical Methods in Engineering. 120 (3): 328—360. doi:10.1002/nme.6138. Архів оригіналу за 6 квітня 2022. Процитовано 6 квітня 2022.

Посилання ред.

Отримано з https://uk.wikipedia.org/w/index.php?title=Метод_матеріальної_точки&oldid=36381504