В работе изложена процедура численного описания поведения упруго-хрупкопластичных материалов. Отмечены основные трудности, возникающие при численном моделировании поведения таких материалов под действием нагрузок, и методы их решения. Кратко описаны некоторые алгоритмы, как известные в литературе, так и разработанные автором, с использованием которых был решен ряд задач. Представлены некоторые результаты расчетов локализации деформации и разрушения геоматериалов.

The processes of localized deformation band formation and crack generation in brittle-plastic materials are simulated numerically in various loading conditions. The patterns of cracking and strain localized band for a sandstone specimen in compression are obtained. To describe inelastic behavior of pressure-sensitive materials, the mathematical tools and concepts of plasticity theory are employed. In the computations the Nikolaevskii model with a non-associated flow rule is used. In the course of “plastic” deformation damages are accumulated. Macrocrack opening and material fracture are described explicitly. For this purpose the node-splitting technique is used and the free-surface conditions are given on newly formed boundaries. This provides self-acting account of crack generation for the entire calculation area in the course of deformation. The use of the considered approaches in describing the behavior of a rock specimen as an ideal homogeneous medium allows a good agreement of fracture patterns as well as of various types of stress-strain dependences with the experimental findings. In order to describe in detail the behavior peculiarities of heterogeneous porous and cracked materials, it is necessary to consider their structural features and crack growth.

В работе представлен анализ нелинейных особенностей поведения горных пород, которые наблюдаются на диаграммах нагружения, но в рамках традиционных моделей упругопластических сред, как правило, игнорируются. Рассмотрены начальный этап деформирования и разгрузка образцов горных пород. Нелинейность поведения на данных этапах нагружения интерпретируется с позиции частичного закрытия трещин, которые образовались в ходе деформировании за пределом упругости или за предыдущую историю нагружения. Предложены феноменологические соотношения, позволяющие учесть указанные нелинейные особенности при численном моделировании. Рассмотрен этап запредельного деформирования, соответствующий стадии локализации деформации и образования магистральных трещин, и предложены поправки для более корректного определения параметров моделей.

Моделирование процессов деформации и разрушения угля выполнено на мезо- и макроуровнях. На мезоуровне численно исследован процесс разрушения мезообъема угля как композиционного материала, в котором явно учтены значимые структурные элементы, включая крупные поры. По известным физико-механическим параметрам компонентов получены усредненные параметры для заданного состава и пористости. На основе полученных усредненных характеристик рассмотрен процесс резания угля на макроуровне, с использованием определяющих соотношений, интегрально учитывающих внутреннее трение и дилатансию в деформируемом материале. Особенности напряженно-деформированного состояния вблизи кромки режущего инструмента воспроизводились при изучении поведения мезообъемов с целью оценки выхода пылевых частиц. Полученные данные о распределении пылевых частиц по размерам и массе хорошо согласуются с экспериментальными данными.

На основе молекулярно-динамического моделирования исследовалось влияние анизотропии динамического нагружения на процессы пластической деформации и разрушения в кристаллите вблизи свободной поверхности. Показано, что существую пороговые величины деформации. которые определяют раскрытие затравочной трещины и ее последующую остановку. Из результатов моделирования следует, что при нагружении в направлении [100] распространение трещины прекращается, после того как вблизи ее вершины начинает формироваться разупорядоченная область.