Full Text

Введение Известно, что макроскопические упругие, реологические и прочностные свойства хрупких материалов чувствительны к скорости деформирования и могут существенно отличаться от соответствующих характеристик механического поведения, определенных в «квазистатических» условиях нагружения [1-9]. У хрупких материалов существенная зависимость характеристик прочностных и реологических свойств от скорости деформации наблюдается в диапазоне от 102 с-1 до 107 с-1. При этом эффективные упругие модули хрупких материалов проявляют заметные изменения при скоростях деформации выше 103 с-1. Распространенным инструментом, используемым для детального изучения особенностей динамического механического поведения материалов, включая динамику разрушения, является численное моделирование в лагранжевых координатах с использованием явных схем интегрирования. Наиболее известным «лагранжевым» численным методом, применяемым для решения таких задач, является метод конечных элементов (МКЭ) [9-13]. Хорошо развитый математический аппарат этого метода позволяет реализовывать сложные реологические модели механического отклика материалов и решать сложные задачи механики деформируемого твердого тела. Отметим, что при использовании этого метода для решения динамических задач возникают и хорошо известные сложности, связанные с моделированием множественного разрушения, сопровождающегося интенсивным массопереносом и перемешиванием. Помимо МКЭ, в настоящее время для моделирования динамических процессов множественного разрушения и переноса фрагментов активно применяются численные методы дискретного подхода к описанию среды, в частности методы дискретных элементов (МДЭ) [14-18]. Несмотря на известные преимущества при моделировании разрушения, МДЭ обладают рядом ограничений, к которым относятся, как правило, более низкий порядок точности по пространству в сравнении с традиционно используемыми конечно-элементными реализациями и недостаточное развитие математического формализма, ограниченного использованием квазистатических механических моделей материалов. Преодоление некоторых важных ограничений этих методов и объединение их преимуществ возможно, например, путем создания комбинированных (МКЭ/МДЭ) методов [19, 20]. Создание подобных «гибридных» методов является возможным вследствие общности базовых принципов реализации моделей упругости, пластичности и разрушения в МКЭ и МДЭ, использующих явные схемы интегрирования. Поэтому уточнение существующих и построение новых реологических моделей, учитывающих особенности динамического отклика хрупких материалов, является актуальным для развития широкого класса «лагранжевых» численных методов. Неупругое поведение хрупких материалов связано главным образом с накоплением микроповреждений (образованием пор и трещин различного масштаба, их ростом или залечиванием, объединением, коллапсом и т.д.) [21-25]. Вовлечение дефектов кристаллической решетки в процесс деформирования таких материалов имеет место только в условиях высоких давлений и температур [24, 26-28]. Поэтому для описания динамического неупругого поведения хрупких материалов применяются модифицированные модели пластичности, которые учитывают высокую чувствительность параметров неупругого отклика к величине давления, а также сложный характер связи между сдвиговой и объемной пластическими деформациями (неассоциированные законы течения). При этом особенности используемого пластического потенциала и предельной поверхности определяют область применения модели. Например, в работе [29] предложена динамическая реологическая модель для моделирования асфальтовых бетонов, основанная на неассоциированном законе пластического течения с предельной поверхностью типа Вермира (Vermeer). Также широкое распространение получили такие модели динамического поведения хрупких изотропных материалов, как модель Тейлора-Чена-Кусзмаула [30, 31], Холмквиста-Джонсона-Кука [32], «Continuous Surface Cap» (CSC) модель [33], модель бетона Karagozian и Case (KCC) [34], модель Риделя-Тома-Хермайера (RHT) [35] и ряд других [36]. Эти модели включают в себя значительное количество модельных параметров и при их корректном определении обеспечивают адекватное описание особенностей динамического деформирования и разрушения хрупких материалов, включая чувствительность свойств к давлению и скорости нагружения. Среди динамических моделей отметим модели, основанные на использовании модифицированных критериев Друкера-Прагера для описания неупругой деформации и разрушения. В частности, одна из таких моделей (Drucker-Prager Cap), модифицированная для описания зависимости поведения хрупких материалов от скорости деформации, описана в работе [37]. Она учитывает роль избыточного давления, приводящего к коллапсу пор [38], размерный эффект [39, 40] и влияние скорости деформации на величину сцепления [41]. Различные модификации модели пластичности Друкера-Прагера широко применяются для моделирования динамического неупругого отклика хрупких материалов [42-44]. Отметим также, что наряду с «изотропными» моделями активно развиваются динамические модели ортотропных хрупких материалов, учитывающие анизотропию упругих и неупругих (в том числе прочностных) механических свойств [9, 45]. В настоящей работе ограничимся рассмотрением локально изотропных материалов, в которых механические свойства элементарных объемов полагаются ориентационно независимыми. В описанных моделях применяется традиционный подход, учитывающий чувствительность характеристик неупругого отклика материала к динамике изменения локального напряженно-деформированного состояния путем введения зависимостей параметров моделей пластичности и разрушения от скорости деформации. В то же время скорость деформации (в отличие от физической массовой скорости) является техническим параметром, характеризующим усредненную по объему скорость изменения размеров образца или его фрагмента. Как правило, при экспериментальном определении зависимостей параметров применяемой модели пластичности или разрушения от скорости деформации используются интегральные значения для всего образца. При этом локальные (например, в области формирования разрушающей трещины) значения могут существенно отличаться от интегральной величины. Альтернативным способом описания поведения материалов в условиях динамического нагружения является подход, развитый в работах таких авторов, как Журков, Регель, Петров, Морозов [46-49]. В рамках данного подхода в качестве основного параметра динамического деформационного отклика материала используется характерное время деградации материала T. Параметр T является физическим и, как правило, характеризует процесс зарождения и развития магистральной трещины (разрушение) в образце или фрагменте материала или формирование системы микроповреждений, обусловливающих макроскопически неупругое поведение этого фрагмента. Очевидно, что временной параметр T является масштабно зависимым, причем иерархия соответствующих масштабов связана с иерархией элементов внутренней структуры материала [50-52]. В работах Журкова этот параметр связывается с термофлуктуационными эффектами в кристаллической решетке и определяется температурой и приложенной нагрузкой. В более поздних работах Петрова и Морозова в качестве характерного времени деградации или разрушения используется интегральный физический параметр с размерностью времени, имеющий более общий смысл и называемый инкубационным временем. Так, в работе [53] предлагается трактовка параметра T как характерного времени релаксации напряжений. Авторы отмечают, что для различных процессов релаксации (пластическая деформация, разрушение) и различных материалов механизмы релаксации напряжений могут качественно различаться, однако это не является препятствием для описания этих процессов на основе общего подхода. В данной работе кинетический подход к построению модели динамического механического поведения материалов реализован на примере модели пластичности Николаевского (неассоциированный закон пластического течения с критерием пластичности в форме Мизеса-Шлейхера) и критерия разрушения Друкера-Прагера [54, 55]. Первоначальная реализация этой модели применима для решения квазистатических задач. В настоящей работе впервые предложена динамическая модификация данной модели, позволяющая учитывать особенности деформирования хрупких материалов при скоростях деформации до ~103 c-1. 1. Описание модели квазистатического механического поведения хрупких материалов Макроскопическое описание неупругой деформации и разрушения хрупких гетерогенных материалов, как правило, осуществляется с использованием моделей механического поведения, включающих первый, второй и третий инварианты тензора напряжений. Это связано с особенностями описания неупругой деформации, обусловленной возникновением и развитием повреждений структуры (трещин различного масштаба, пор). Рост концентрации и размеров повреждений обусловливает эффект дилатансии, т.е. возрастание объемной неупругой деформации. В данной работе рассматривается одна из популярных моделей неупругого деформационного отклика хрупких материалов - модель Николаевского (неассоциированный закон пластического течения). Реализации модели в рамках различных лагранжевых численных методов успешно используются для изучения особенностей деформирования хрупких гетерогенных материалов и сред [24, 54-57]. В основе этой модели лежит неассоциированный закон пластического течения с критерием пластичности в форме Мизеса-Шлейхера , где J1 и J2 - первый инвариант тензора напряжений и второй инвариант девиатора тензора напряжений; w - безразмерный параметр, пропорциональный коэффициенту внутреннего трения a (a = 3w); smean - среднее напряжение; seq - интенсивность напряжений. Условием наступления предельного состояния является выполнение неравенства (1) где Y - сцепление. В модели Николаевского постулируется линейная связь между объемной и сдвиговой скоростью пластической деформации (2) где - первый инвариант тензора скорости пластической деформации; - второй инвариант девиатора тензора скорости пластической деформации; L - коэффициент дилатансии. Для реализации модели Николаевского в рамках численных методов, использующих явные схемы интегрирования, широко применяется алгоритм Уилкинса (Wilkins) [58]. В рамках алгоритма Уилкинса решение упругопластической задачи разбивается на два этапа: 1) Решение упругой задачи в приращениях. Примеры численной реализации упругой задачи в рамках МДЭ можно найти в более ранних работах авторов [54, 59]. Рис. 1. Схематичное описание работы алгоритма Уилкинса Fig. 1. Scheme of the Wilkins algorithm 2) Приведение компонент девиатора тензора напряжений к поверхности текучести в случае выполнения неравенства (1) на текущем шаге интегрирования уравнений движения (рис. 1). Математически приведение можно выразить в терминах девиаторов напряжений: , (3) где - компоненты девиатора тензора напряжений (исходные значения, полученные после решения упругой задачи); D¢ab - приведенные (скорректированные) компоненты девиатора тензора напряжений (здесь и далее символом «¢» обозначены скорректированные значения соответствующих величин), M£1 - коэффициент приведения (коррекции) напряжений. Таким образом, в области упругого отклика моделируемого элементарного объема хрупкого материала в качестве определяющего соотношения используется линейный закон Гука. Для описания релаксации напряжений за пределом упругости применяется алгоритм Уилкинса [58]. В рамках модели Николаевского приведение тензора напряжений к соответствующей точке предельной поверхности осуществляется с использованием следующего выражения: , (4) где ; - символ Кронекера. Приведенные («исправленные») значения инвариантов тензора напряжений имеют вид [55] (5) где G - модуль сдвига; K - модуль объемного сжатия. Отметим, что в приведенной модели величина упругих модулей элементарного объема материала полагается постоянной и не зависящей от величины накопленной неупругой деформации (отражающей степень поврежденности элементарного объема). В настоящей работе в качестве критерия локального разрушения используется двухпараметрический критерий Друкера-Прагера, учитывающий характерную для хрупких материалов чувствительность сдвиговой прочности к гидростатическому давлению и имеющий тот же структурный вид, что и критерий Мизеса-Шлейхера: , (6) где sc - предел прочности на сжатие; , st - предел прочности на растяжение; и среднее напряжение и интенсивность напряжений. 2. Описание модели динамического механического поведения хрупких материалов Главной идеей кинетической модели динамического неупругого поведения хрупких материалов, развиваемой в представленной работе, является учет конечного времени процесса релаксации локальных напряжений, обусловленного возникновением повреждений различного масштаба. Описанная выше статическая модель адекватно описывает неупругий отклик материала при нагружении с малыми скоростями деформации, когда инерционными эффектами можно пренебречь (квазистатический режим). В рамках квазистатического приближения фактически полагается, что релаксация локальных напряжений и связанное с этим перераспределение напряжений в образце осуществляются за пренебрежимо малое время (при численном моделировании - за один шаг интегрирования по времени). Однако при динамическом нагружении вклад инерционных членов в перераспределение кинетической и потенциальной энергии в объеме образца становится значимым или даже определяющим. Поэтому в динамическом режиме нагружения характеристики интегрального отклика образцов оказываются зависимыми от скорости деформирования. Сказанное обосновывает необходимость учета конечного времени релаксации локальных напряжений, поскольку его величина будет определять динамику перераспределения напряжений. В современных кинетических моделях неупругости и прочности хрупких материалов, как правило, применяются интегральные критерии разрушения, основанные на вычислении приращения импульса нагружения за некоторое время t, называемое временем инкубации разрушения [49, 60]. Параметр t интерпретируется как временной период, за который формируется разрушающая трещина при постоянной нагрузке, равной статическому пределу прочности материала Такой интегральный критерий разрушения можно записать в виде (7) где - скалярный силовой параметр напряженного состояния; t - время. В различных реализациях критерия в качестве силового параметра используется максимальное главное напряжение либо комбинация инвариантов тензора напряжений. Очевидно, что выполнение критерия (7) возможно в случае, если приложенная нагрузка достигнет или превысит величину статического предела прочности материала и не опустится ниже в течение некоторого конечного промежутка времени. При этом, если приложенная нагрузка на протяжении времени t остается постоянной и равна статическому пределу прочности материала, разрушение произойдет за время t (характерное время инкубации разрушения). Экспериментальное определение величины параметра t является сложной задачей. Более того, его количественное значение зависит от пространственного масштаба рассматриваемого образца [49]. Поэтому в данной работе предлагается упрощенная реализация динамического критерия разрушения (7): , (8) где - динамическое значение предела прочности материала, H(t0 + t - t) - функция Хевисайда. Такая формулировка критерия (7) подразумевает, что при достижении скалярным силовым параметром s статического предела прочности () в момент времени t0 начинается процесс разрушения. Длительность этого процесса определяется динамикой изменения параметра s при t > t0. К моменту формирования разрушающей трещины (при t = tfr) величина параметра s составит некоторое значение ³ (конкретное значение зависит от динамики изменения напряженного состояния в период времени (tfr - t0)). Предположим, что за время формирования трещины скорость деформации не претерпевает существенных изменений. В этом случае скорость роста параметра s за время (tfr - t0) также не претерпит существенных изменений. В рамках данного предположения время разрушения (tfr - t0), локальная скорость деформации и динамический предел прочности однозначно связаны между собой. При этом значение любого параметра из тройки (, (tfr - t0), ) может быть однозначно определено, если известны значения двух других параметров. В случае одноосного сжатия с постоянной скоростью время разрушения может быть вычислено по формуле , (9) где и - статическое и динамическое значения предела прочности материала при одноосном сжатии; E - модуль Юнга; - скорость осевой деформации; Tf - время разрушения (в общем случае отличное от времени инкубации). Аналогичное соотношение можно записать и для одноосного динамического растяжения с постоянной скоростью. Отметим, что к настоящему времени накоплен достаточно большой массив данных о зависимостях динамической прочности образцов хрупких материалов от скорости деформации при одноосном сжатии и растяжении [1-8]. Наличие этих данных позволяет получить зависимости Tf () в указанных условиях нагружения. Важно отметить, что эти зависимости являются нелинейными и монотонно убывающими (таким образом, функция Tf () является однозначной). Пусть для некоторого заданного вида напряженного состояния (в общем случае отличного от одноосного нагружения) известна зависимость Tf (), полученная при постоянной скорости деформирования. Положим, как и для случая одноосного нагружения, что эта зависимость является монотонно убывающей. Тогда сделанное выше допущение постоянства скорости роста локальных напряжений в течение временного интервала Tf позволяет предложить следующую численную реализацию динамического критерия разрушения (8). В рамках данной реализации полагаем, что разрушение материала происходит не в момент t0 достижения статического предела прочности s = , а спустя некоторое время, в течение которого в рассматриваемой локальной области материала протекают внутренние переходные процессы (инкубация разрушения). Результатом этих процессов станет образование новых свободных поверхностей (разрушение материала), приводящее к падению сопротивления внешнему нагружению. При этом в течение времени разрушения Tf рассматриваемая локальная область сохраняет целостность. Как отмечалось выше, сделанное приближение постоянства скорости деформации на протяжении Tf обеспечивает единственное значение динамической прочности, отвечающее текущей скорости изменения скалярного силового параметра s, и, следовательно, единственное значение Tf (рис. 2). Тогда для определения момента окончания подготовки разрушения в каждый момент времени t > t0 после достижения статического предела прочности производится проверка неравенства , (10) Рис. 2. Схема динамического разрушения Fig. 2. Scheme of dynamic fracture где s(t) - текущее локальное значение силового параметра; Tf = t - t0, а значение определяется из «калибровочной» зависимости Tf(). При выполнении условия (10) (в момент времени tfr) динамическая прочность материала считается достигнутой, процесс инкубации завершенным, и происходит локальное разрушение материала. Соответствующие величины Tf = tfr - t0 и рассматриваются как локальные значения времени разрушения и динамической прочности. Отметим, что предложенная схема численной реализации динамического критерия (8) является общей для силовых параметров различного структурного вида. В настоящей работе в качестве силового параметра s используется предложенная Друкером и Прагером комбинация первых двух инвариантов тензора напряжений. При этом формулировка условия (10) принимает следующий вид: , (11) где = (t - t0) - динамический предел прочности на сжатие, = (t - t0) - динамический предел прочности на растяжение. Применение критерия (11) возможно в случае, если известны зависимости (Tf) и (Tf). Экспериментальные данные свидетельствуют о том, что чувствительными к скорости нагружения являются не только прочностные, но и неупругие свойства хрупких материалов, такие как сцепление Y, коэффициент внутреннего трения a и коэффициент дилатансии L [5]. Это связано с тем, что неупругий отклик хрупких материалов связан с возникновением системы внутренних несплошностей (повреждений и микротрещин), формирование которых длится конечное время. Для адекватного описания неупругого деформирования хрупких материалов при динамическом нагружении необходимо учитывать зависимости величин параметров модели от времени формирования системы несплошностей. Поскольку неупругий отклик материала связан с релаксацией локальных напряжений, время формирования системы микроповреждений далее будем называть временем релаксации напряжений или сокращенно - временем релаксации Tr. Как и в случае динамического разрушения, конкретное значение времени релаксации зависит от динамики изменения напряженного состояния. В частном случае постоянства напряженного состояния в точке, отвечающей пределу текучести, время релаксации является кинетической константой материала, аналогично времени инкубации разрушения t. Так как в рассматриваемом случае механизмы неупругой деформации аналогичны механизмам локального разрушения, модель динамического неупругого отклика хрупких материалов также может быть построена на основе формализма кинетической теории прочности. В рамках данной модели предполагается, что скорость деформации и скорость изменения параметра напряженного состояния s остаются постоянными в течение времени релаксации (конкретный вид параметра s может быть отличным от того же в динамическом критерии разрушения). В этом случае динамический предел текучести является однозначной монотонно убывающей функцией времени релаксации. Для частного случая одноосного нагружения вид данной зависимости можно получить из уравнения, аналогичного (9), с использованием пределов текучести вместо пределов прочности и модулем вместо E. В общем случае форма записи динамического критерия пластичности идентична (8): , (12) где - динамический предел текучести. Конкретный вид этого критерия определяется применяемой моделью пластичности. Как и в модели разрушения, полагаем, что релаксация напряжений происходит не в момент t0 достижения статического предела текучести s = , а спустя определенное время (время релаксации). Алгоритм Уилкинса заключается в решении упругой задачи в приращениях на каждом шаге интегрирования уравнений движения и приведении вычисленных таким образом напряжений к соответствующей точке предельной поверхности согласно соотношениям (4)-(5) при выполнении неравенства F ³ Yst. В численной реализации статической модели приведение осуществляется за шаг интегрирования по времени. В динамической модели напряжения приводятся к поверхности текучести постепенно, в течение времени релаксации. Полагаем, что для рассматриваемого материала известна зависимость величины сцепления от времени релаксации Ydyn(Tr), которая является монотонно убывающей функцией. При выполнении неравенства F ³ Yst на данном шаге интегрирования текущее значение параметра F рассматривается как текущее (динамическое) значение предела упругости на сдвиг: Ydyn = F. При этом использование «калибровочной зависимости» Ydyn(Tr) позволяет определить текущее (отвечающее текущей динамике изменения напряжений) значение времени релаксации Tr = Tr(F). Полагая, что релаксация напряжений (т.е. их приведение к точке Yst предельной поверхности) на протяжении периода времени Tr происходит по линейному закону, уменьшение величины параметра F за шаг по времени Dt можно записать на основе линейной интерполяции: . (13) Приращение DF используется в соотношениях (4)-(5) динамической реализации модели пластичности Николаевского вместо «статической» разницы (F-Y). Текущие значения других параметров модели пластичности (wdyn и Ldyn) в соотношениях (4)-(5) определяются на основе «калибровочных» зависимостей wdyn(Tr) и Ldyn(Tr) с использованием текущего значения времени релаксации Tr. Здесь предполагается, что кривые wdyn(Tr) и Ldyn(Tr) известны и имеют монотонный характер. Важно отметить, что калибровочные зависимости Ydyn(Tr), wdyn(Tr) и Ldyn(Tr) должны соответствовать друг другу, т.е. конкретные значения Ydyn, wdyn и Ldyn при некотором Tr должны отвечать одному и тому же состоянию материала (т.е., быть полученными при одних и тех же условиях нагружения). В общем случае величина сцепления является функцией не только времени релаксации (или скорости деформации), но и функцией деформации (Ydyn = Ydyn(εpl, Tr)). Так, в модели Николаевского функция Y(ε) определяет упрочнение материала за пределом упругости. Известно, что при динамическом нагружении материала функция упрочнения претерпевает значительные изменения. А именно наряду с ростом предельных напряжений (пределов упругости, прочности) наблюдается снижение доли пластической деформации материала и, как следствие, увеличение наклона функции упрочнения. Поэтому для адекватного моделирования деформационного отклика материала при динамических воздействиях необходимо учитывать снижение доли пластической деформации при высоких скоростях нагружения. В данной работе предлагается учитывать снижение доли пластической деформации путем введения функции εpl(Tr) с использованием того же временного параметра Tr, что и для остальных параметров модели. При этом очевидно, что функция εpl(Tr) должна быть монотонно возрастающей, то есть стремиться к верхнему пределу при стремлении к квазистатическому режиму нагружения. Как уже отмечалось ранее, упругие свойства материала также чувствительны к скорости деформирования, однако этот эффект проявляется при достаточно высоких скоростях деформации (>103 c-1) и в рамках данной работы не рассматривается. Следует отметить, что предложенная модель динамического поведения хрупких материалов фактически представляет собой обобщение статических моделей неупругого деформирования и разрушения материала, учитывающее конечное (хотя и малое) время релаксационных процессов. Это обеспечивает возможность моделирования как квазистатических, так и динамических процессов в рамках единого формализма. 3. Определение зависимостей неупругих/прочностных параметров модели хрупких материалов от времени релаксации/разрушения Предложенная модификация модели Николаевского позволяет описывать механическое поведение хрупких гетерогенных локально изотропных материалов в условиях динамического нагружения. Однако для решения конкретных задач необходимо знать вид зависимости каждого из параметров модели от времени релаксации. В идеальном случае эта информация должна быть получена на основе обобщения экспериментов по динамическому деформированию образцов рассматриваемого материала требуемого масштаба при использовании различных видов напряженного состояния. Как уже отмечалось, эксперименты по определению времени разрушения и времени релаксации для различных видов напряженного состояния и различных материалов весьма сложны. Это обусловливает нехватку данных о зависимости прочностных и реологических характеристик материалов от времени разрушения и релаксации. С другой стороны, в последние десятилетия накоплена большая база экспериментальных данных о зависимости свойств хрупких материалов от скорости деформации в условиях одноосного деформирования и при некоторых других видах нагружения. Анализ этих данных показывает, что, несмотря на различие в абсолютных значениях величин динамических характеристик, зависимости приведенных динамических значений этих параметров (полученных нормированием на соответствующие статические значения) от скорости деформации могут быть с хорошей точностью аппроксимированы едиными зависимостями для широкого класса хрупких материалов [1-5]. Используя единые зависимости свойств материалов от скорости деформации, можно получить (например, по формуле (9)) требуемые зависимости параметров модели от времени релаксации и разрушения для конкретного хрупкого материала. В данной работе для получения таких зависимостей использованы экспериментальные данные по динамическому нагружению хрупких материалов в условиях одноосного сжатия и растяжения, собранные и обобщенные в работе [5]. Из экспериментальных данных хорошо известны зависимости пределов прочности на сжатие и растяжение от скорости деформации для широкого спектра хрупких материалов. В настоящей работе использовались обобщенные в [5] зависимости пределов прочности на сжатие и растяжение от скорости деформации (рис. 3, а, кривые 1 и 2 соответственно): , (14) , (15) где - скорость деформации в единицах с-1; = 1 с-1 - нормировочный коэффициент. Рис. 3. Зависимости пределов прочности при сжатии и растяжении (а) и сцепления, коэффициента внутреннего терния и пластической деформации (б) от скорости деформации для бетонов и горных пород Fig. 3. Dependencies of the compressive and tensile strengths (a) and cohesion, internal friction coefficient and plastic strain (b) on strain rate for concretes and rocks Эти данные о динамической прочности хрупких материалов можно представить в терминах характерного времени релаксации напряжений (в данном случае - времени разрушения). Так, для определения характерного времени разрушения конкретного материала (например, песчаника E = 16 ГПа, = 70МПа, = 31,5 МПа) в условиях одноосного сжатия использовалась формула (9) . (16) Полученная зависимость является убывающей и носит ярко выраженный нелинейный характер. При малых значениях времени разрушения (соответствующих высоким скоростям деформаций) она быстро убывает, а при достаточно больших значениях Tf (соответствующих квазистатическим скоростям деформации) асимптотически стремится к единице. Построенный набор точек , соответствующий интервалу скоростей деформации 10-3 c-1 - 103 c-1, не удается удовлетворительно аппроксимировать аналитической функцией простого вида. Наилучшая аппроксимация возможна в случае разделения зависимости на два участка. Область малых значений Tf аппроксимировалась логарифмической зависимостью, а область высоких значений Tf (пологий участок) аппроксимировалась линейной функцией (17) где время разрушения Tf приводится в секундах, Tf 0 = 0,00175 с. Аналогичным образом для песчаника была определена зависимость , также аппроксимируемая двумя функциями в области малых и больших значений Tf: (18) Выражения (17)-(18) получены для песчаника на основе единых кривых зависимости прочности от скорости деформации. Для обобщения (17)-(18) на класс хрупких материалов в качестве аргумента функций должны использоваться «масштабированные» параметры в (17) и в (18), где , , - предел прочности на сжатие для «опорного» материала; - предел прочности на растяжение для «опорного» материала; Eref - модуль Юнга опорного материала; Tfref - время разрушения для «опорного» материала в соответствующем виде нагружения. В данном случае в качестве опорного материала выступает песчаник. Для моделирования зависимости величины сцепления от скорости деформации также использовались сведенные в работе [5] экспериментальные данные для различных хрупких материалов, с хорошей точность аппроксимируемые зависимостью (см. рис. 3, б, кривая 1) (19) где = 317 с-1 - нормировочный коэффициент. Переход от скорости деформации к времени релаксации для произвольного хрупкого материала осуществлялся с помощью формулы, аналогичной (9): , (20) . (21) Так как полагается, что исходная зависимость описывает широкий класс хрупких материалов, функция является инвариантной частью в выражении для определения времени релаксации Tr. Множитель по своей сути является масштабным коэффициентом, определяющимся «статическими» характеристиками материала. В первом приближении зависимость можно аппроксимировать линейной функцией вида , (22) числовые коэффициенты в которой имеют размерность времени. Определяемая из (20) зависимость используется в формуле (13) для определения приращения параметра Ф при численной реализации динамической модели пластичности. Остальные параметры модели (a, L, e pl) также являются функциями скорости деформации, а следовательно, и времени релаксации Tr. В частности, известно, что зависимость приведенной величины коэффициента внутреннего трения от скорости деформации является единой для широкого класса хрупких материалов [5] и может быть аппроксимирована функцией (см. рис. 3, б, кривая 2): , (23) где - скорость деформации в единицах с-1 и использован тот же нормировочный множитель в показателе экспоненты, что и в (19). Используя определенную на основе (20) связь между Tr и , можно получить зависимость приведенного коэффициента внутреннего трения от времени релаксации для «опорного» материала. Расширить область применения этой зависимости на класс хрупких материалов можно путем введения безразмерного материального коэффициента: . (24) Тогда выражение для нахождения времени релаксации для произвольного хрупкого материала может быть записано в виде (25) где (26) Зная зависимость приведенного коэффициента внутреннего трения от времени релаксации для опорного материала и используя коэффициент R, можно построить «единую» зависимость приведенного коэффициента внутреннего трения от времени релаксации (в качестве «опорного» материала использовался песчаник (Yst)ref = 15,4 МПа, Eref = 16 ГПа): , (27) где время разрушения приводится в секундах, а Tr0 = 6,7×10-6 c. Зависимость коэффициента дилатансии от времени релаксации может быть получена аналогичным образом. Однако ввиду отсутствия исчерпывающей информации о зависимости коэффициента дилатансии от скорости деформации в представленных далее тестовых расчетах использована зависимость величины коэффициента дилатансии от времени релаксации того же вида, что и для коэффициента внутреннего трения: . (28) Учет зависимости величины вклада пластической составляющей в полную деформацию от в представленных ниже расчетах осуществлялся путем введения экспоненциального закона снижения этого вклада с ростом (см. рис. 3, б, кривая 3): (29) где - скорость деформации в единицах с-1 и использован тот же нормировочный множитель в показателе экспоненты, что и в (19). Как и в случае коэффициента внутреннего трения, была получена зависимость для «опорного» материала и обобщена путем введения материального коэффициента R: (30) где время разрушения приводится в секундах и используется тот же нормировочный коэффициент, что и в выражениях (27)-(28). 4. Верификация динамической модели Для демонстрации корректности и адекватности предложенной модели динамического деформирования хрупких материалов был проведен ряд численных тестов с использованием одного из представителей класса лагранжевых численных методов, а именно метода подвижных клеточных автоматов (англ. Movable cellular automaton method - MCA [54, 61, 62]). Метод подвижных клеточных автоматов является активно развивающимся представителем группы МДЭ, основными особенностями которого являются приближение однородного распределения деформаций и напряжений в объеме элемента и многочастичная формулировка межэлементного взаимодействия, полученная на основе использования модели клеточных автоматов Виннера-Розенблюта [63, 64]. Моделировались испытания на одноосное сжатие и растяжение с постоянной скоростью двумерных образцов прямоугольной формы размером 6×9 мм. Скорость деформирования варьировалась в пределах от 10-3 с-1 (что полагалось соответствующим квазистатическому нагружению) до 103 с-1. В качестве предела прочности образца принималась максимальная величина регистрируемого сопротивления нагружению образца, отнесенная к площади торцевой поверхности. Рассматривались образцы двух модельных материалов (таблица), механические параметры которых близки к характерным механическим свойствам песчаника и бетона нормальной прочности. Механические характеристики тестовых образцов бетона и песчаника Mechanical properties of concrete and sandstone samples for tests Образцы Коэффициент Пуассона Плотность, кг/м3 E, ГПа syst, МПа scst, МПа stst, МПа a L Бетон 0,194 4660 38,6 21,8 45,8 21 0,63 0 Песчаник 0,28 2200 16 40 70 31,5 0,57 0,36 На рис. 4 приведены расчетные диаграммы одноосного сжатия модельных образцов бетона и песчаника при различных скоростях деформации. Представленные диаграммы получены в приближении плоского напряженного состояния. Рисунки демонстрируют корректное изменение интегральных свойств материала с ростом скорости деформирования, а именно увеличение пределов упругости и прочности, а также уменьшение пластической составляющей деформации (увеличение коэффициента деформационного упрочнения). На основе анализа диаграмм одноосного сжатия и растяжения образцов рассматриваемых материалов были получены зависимости пределов прочности на сжатие и растяжение от скорости деформации (рис. 5). Можно видеть, что результаты численного моделирования хорошо согласуются с экспериментальными данными [5], которые использовались для получения зависимостей времени релаксации от параметров модели (Tr(Y), , ). Это подтверждает адекватность развитой динамической численной модели неупругого деформирования и разрушения хрупких материалов. Для наглядной демонстрации необходимости учета изменения прочностных и неупругих характеристик материалов с ростом скорости деформирования было проведено моделирование испытаний образцов бетона на одноосное сжатие, в котором зависимости параметров модели от скорости деформирования не учитывались. Скорость деформации образцов варьировалась в пределах от 10-3 с-1 до 103 с-1. Результаты, полученные в ходе численных экспериментов, приведены на рис. 6. Хорошо видно, что при использовании «квазистатических» значений параметров модели изменение прочностных характеристик образцов с ростом скорости деформации является несущественным. Только в области высоких скоростей деформации имеет место небольшое увеличение прочности образцов. Данный эффект является следствием инерционности образцов, которая начинает играть значимую роль при высокоскоростном деформировании материала. Таким образом, результаты моделирования демонстрируют важность и наглядно подтверждают необходимость учета конечного времени формирования повреждений и трещин различного масштаба в хрупких материалах при решении задач динамического нагружения. Рис. 4. Диаграммы одноосного сжатия образцов бетона (а) и песчаника (б) при скоростях деформации : 10-3 с-1 - 1; 10 с-1 - 2; 100 с-1 - 3 Fig. 4. Uniaxial compression loading diagrams of concrete (a) and sandstone (b) samples at various strain rates : 10-3 s-1 - 1; 10 s-1 - 2; 100 s-1 - 3 Рис. 5. Зависимости приведенного предела прочности на сжатие (а) и растяжение (б) образцов бетона и песчаника от скорости деформации: - бетон; - песчаник. Пунктирной линией показана аппроксимация экспериментальных данных [5] Fig. 5. Numerical tests results: dependencies of specified compressive strength (a) and specified tensile strength (b) of the concrete and sandstone samples on the strain rate: - concrete; - sandstone. Dash line is the experimental data approximation [5] Рис. 6. Зависимость предела прочности на сжатие от скорости деформации, полученная без учета зависимости свойств материала от скорости деформации, для образцов бетона и песчаника: - бетон; - песчаник. Пунктирной линией показана аппроксимация экспериментальных данных [5] Fig. 6. Numerical tests results with no dependence of model parameters on strain rate: dependence of specified compressive strength of the concrete and sandstone samples on the strain rate: - concrete; - sandstone. Dash line is the experimental data approximation [5] Динамика и итоговая картина разрушения являются важными результатами моделирования деформирования и разрушения хрупких материалов. При моделировании испытаний на одноосное сжатие в модельных образцах бетона и песчаника наблюдалась классическая картина разрушения: при «квазистатических» скоростях сжатия образцы разрушаются путем формирования магистральных трещин, от которых развиваются вторичные трещины. С ростом скорости деформации усиливается степень фрагментации образцов, а также появляются эффекты откола небольших фрагментов с боковых поверхностей. Подобные закономерности и особенности разрушения наблюдаются и в испытаниях реальных образцов хрупких материалов [36, 65]. Полученные результаты свидетельствуют, что развитая в работе динамическая модель может применяться для изучения зависимостей величины интегральной энергии разрушения «квазихрупких» систем и распределения фрагментов по размерам от амплитуды и скорости нагружения. Заключение Основной идеей развитого подхода к построению моделей динамического механического поведения хрупких материалов является учет конечного времени зарождения несплошностей различного масштаба (времени релаксации Tr и времени разрушения Tf) в материале на основе кинетической теории прочности. В отличие от традиционных динамических моделей, здесь основным управляющим параметром, определяющим динамику отклика материала, является не локальная скорость деформации, а физический параметр - время зарождения несплошностей. Отметим, что при этом не требуется строгого разделения задач деформирования материала (и соответствующих математических моделей) на статическую и динамическую и тем самым обеспечивается возможность корректного перехода от режима квазистатического нагружения в высокоскоростной динамический режим и наоборот. Более того, хотя развитый подход реализован на примере конкретных моделей пластичности и разрушения хрупких материалов, он естественным образом может быть применен к различным моделям неупругого отклика и разрушения, реализуемым в рамках различных лагранжевых численных методов (включая не только методы конечных и дискретных элементов, но и методы конечных разностей). Дальнейшее развитие этого подхода может быть связано с учетом зависимостей величин упругих модулей материала от степени поврежденности (что позволит более точно описывать упругий отклик материала при относительно больших неупругих деформациях) и анизотропии локальных механических свойств. Поскольку управляющие параметры модели (характерные времена зарождения несплошностей различного масштаба) определяются масштабно-зависимой величиной времени инкубации [53], они также являются масштабно-зависимыми. Поэтому при экспериментальном определении ключевых зависимостей параметров численной модели от времени зарождения несплошностей необходимо учитывать масштабный фактор. В частности, характерный размер образцов, использующихся для определения зависимостей прочностных и реологических свойств от скорости деформации (или времени зарождения несплошностей), должен быть сопоставим с размером конечного/дискретного элемента, поведение которого будет определяться этими зависимостями. Развитый в работе подход к построению динамических моделей является актуальным для решения нового класса прикладных задач, связанных с природными и техногенными динамическими воздействиями на конструкции из искусственных строительных материалов, включая бетоны, элементы конструкций из керамики и природные каменные материалы. В частности, реализация предложенного подхода в рамках метода частиц позволяет решать задачи прогнозирования времени разрушения материалов и конструкций в зависимости от амплитуды и скорости нагружения с учетом структурных особенностей рассматриваемых материалов. При этом на основе решения обратной задачи предложенный подход дает возможность проводить оценки характерного размера элементов внутренней структуры материала, определяющих разрушение на рассматриваемом масштабном уровне.

About the authors

A S Grigoriev

Institute of Strength Physics and Materials Science

E V Shilko

Institute of Strength Physics and Materials Science; National Research Tomsk State University

V A Skripnyak

National Research Tomsk State University

A G Chernyavsky

S.P. Korolev Rocket and Space Corporation “Energia”

S G Psakhie

Institute of Strength Physics and Materials Science; National Research Tomsk State University

References

Statistics

Views

Abstract - 290

PDF (Russian) - 111