ON THE CONVEXITY OF POTENTIAL OF NONLINEAR ELASTIC MEDIA MODEL WITH TENSOR DAMAGE PARAMETER

Abstract


For solving problems of the theory of elasticity using nonlinear models, the question of the convexity of the potential used and the proof of the uniqueness of the solution are of decisive importance. The present work is devoted to the determination of conditions for the local strict convexity of the potential for the nonlinear elasticity model, which ensure the uniqueness of the solution to the problem in a sufficiently small neighborhood of the desired solution. The considered model is a generalization of a scalar nonlinear rheological model of brittle solid deformation for the case of the tensor damage parameter, the principal values of which describe the reduction in the cross-sectional area of the material in three orthogonal directions. An additional term of the second order in deformation in the elastic potential makes it possible to describe the dependence of the elastic moduli on the type of the stress-strain state, the dilatancy of the material under shear deformation, as well as the nonlinear deformation response even at low loads. The introduced damage tensor of the second rank makes it possible to describe the damage-induced anisotropy of the elastic properties of the material. Conditions of local strict convexity in the principal axes of the strain tensor are obtained in this work for the general case of misaligned strain and damage tensors. To illustrate the obtained convexity conditions, two special cases of the damage tensor type are considered: a transversely isotropic fractured medium with coaxial strain and damage tensors, and a transversely isotropic medium with obliquely oriented fracturing. For both cases, the dependence of the limiting values of damage on the parameter of the degree of anisotropy is shown. It is shown that in the case of weak damage anisotropy, the potential convexity conditions for the scalar damage parameter give a minorant estimate of the maximum allowable damage for various types of stress-strain state. For obliquely oriented fracturing, the dependences of the maximum permissible damage on the degree of anisotropy, the angle of inclination and the type of stress-strain state are plotted.

Full Text

Введение Нелинейным деформационным откликом на приложенное усилие обладает большинство горных пород [1-8], сыпучих сред [9-10] и волокнистых/зернистых композитов [11-14]. Для геоматериалов нелинейный отклик не зависит от пространственного масштаба и наблюдается вплоть до геологических структур километровых размеров [15-18]. Еще одной отличительной особенностью материалов такого типа является их разномодульность, в общем случае зависимость упругих свойств от вида напряженно-деформированного состояния. Для природных материалов зависимость свойств от вида напряженно-деформированного состояния обусловлена развитой дефектностью, сформированной как на стадии литогенеза, так и на стадии активного деформирования в тех или иных геологических условиях. В условиях сжатия/растяжения трещиноподобные дефекты, ориентированные нормально к действию максимального/минимального сжимающего напряжения, схлопываясь/открываясь, дают дополнительный вклад в общую деформацию материала. Для описания деформационного поведения нелинейных разномодульных сред используются различные подходы, которые можно разделить на два класса: модели с упругими модулями, зависящими от типа напряженно-деформированного состояния [19-26], модели с постулированием упругого потенциала особой формы [27-34]. В независимости от способа описания нелинейного деформационного отклика при решении краевых задач определяющее значение имеет вопрос устойчивости материала, которая определяется, в свою очередь, выпуклостью используемого потенциала, обеспечивающую единственность решения, и эллиптичностью уравнения движения. Игнорирование проверки обеспечения в моделях данного класса фундаментальных принципов (выпуклости потенциала, условия эллиптичности) приводит к проблемам сходимости при численной реализации решения конкретной краевой задачи [35; 36]. Настоящая работа посвящена исследованию условий локальной строгой выпуклости потенциала в неклассической модели нелинейной упругости с тензорным параметром поврежденности, обеспечивающих единственность решения упругой задачи в достаточно малой окрестности искомого решения (единственность в малом). Получены соотношения в общем случае несоосных тензора деформации и поврежденности. В качестве примера рассмотрен случай трансверсально-изотропной среды с трещиноватостью, ориентированной произвольно относительно эллипсоида действующих деформаций. 1. Нелинейная неклассическая модель деформирования хрупкого тела со скалярным параметром поврежденности В работе [37] на основе множества экспериментальных данных предложен общий вид упругой энергии в предположении, что это функция второго порядка инвариантов тензора деформации. Приведенные в работе выкладки обеспечили строгую математическую основу для формулировки энергии упругой деформации твердого тела, обсуждавшейся ранее в работах [32-33; 38-41]: (1) где - параметр Ламе, - дополнительный упругий модуль, определяющий степень нелинейности материала, - первый и второй инвариант тензора деформации соответственно. Отличительной особенностью соотношения (1) от традиционных моделей нелинейной упругости является то, что дополнительное слагаемое имеет тот же порядок, что и основные слагаемые. Традиционные модели нелинейной упругости базируются на включении в выражение для упругой энергии слагаемых более высоких порядков по деформации [49]. В исследовании [33] показано, что дополнительное слагаемое является следствием принятия гипотезы об избирательной активизации микротрещин при приложении нагрузки: активизируются только микротрещины, берега которых субпараллельны направлению действия максимального сжимающего напряжения. Дополнительное слагаемое второго порядка обеспечивает нелинейную связь компонент тензора напряжения и компонент тензора деформации и позволяет описать нелинейный отклик материала уже при малых нагрузках: (2) Упругие модули зависят от вида напряженно-деформированного состояния, определяемого параметром Параметр меняется от для всестороннего сжатия до при всестороннем растяжении; соответствует одноосному сжатию/растяжению, а - чистому сдвигу. В работах [32; 42] в модель введена зависимость упругих свойств от скалярного параметра поврежденности, описывающего плотность микротрещин, в виде где - материальный параметр, контролирующий переход от залечивания микротрещин к их росту. Параметр связан с упругими модулями и углом внутреннего трения соотношением: , (3) где и определяется по результатам традиционных испытаний материала по схеме Кармана. Устойчивость материала в процессе деформирования определяется двумя типами условий. Первый тип - условия выпуклости упругого потенциала, обеспечивающие единственность решения статической задачи теории упругости [47]. Формулировка и обобщение этих условий, в том числе для случая нелинейных моделей, представлены в работах [29; 32; 43; 48]. Второй тип условий - условия эллиптичности уравнения движения теории упругости [49]. Условия выпуклости и эллиптичности не всегда идентичны, особенно для нелинейной упругости [50]. Условия выпуклости являются более сильными, и выпуклость может быть потеряна раньше, чем эллиптичность. Настоящая работа посвящена исследованию условий выпуклости потенциала, обеспечивающих единственность решения статической задачи теории упругости. Необходимыми и достаточными условиями локальной строгой выпуклости потенциала (1) являются условия положительной определенности матрицы Гессе функции (4) В свою очередь матрица будет положительно определенной в случае положительности всех ее собственных значений [51]. В работах [32; 43] показано, что условия положительной определённости матрицы Гессе для потенциала вида (1) имеют вид , (5) где, , - эффективные упругие модули среды, a - параметр поврежденности. Условия (5) накладывают ограничения как на дополни- Рис. 1. Зависимость допустимых значений поврежденности α, определяемых условиями (5) локальной выпуклости потенциала, от вида напряженно-деформированного состояния Fig. 1. Dependence of the permissible values of damage α, determined by the conditions (5) of the local potential convexity, on the type of stress-strain state тельный упругий модуль , так и на допустимые максимальные (предельные) значения поврежденности a в зависимости от параметра вида напряженно-деформированного состояния x. Это обстоятельство является отличительной чертой нелинейных моделей. В моделях классической континуальной механики поврежденности разрушение материала происходит при , тогда как в нелинейных моделях максимальная поврежденность, характеризующая макроразрушение материала, определяется условием выпуклости потенциала. Нетрудно показать, что наиболее строгим условием для дополнительного упругого модуля является третье неравенство в (5). Минимальное допустимое значение определяется этим неравенством при и : (6) Для определения зависимости допустимого максимального значения поврежденности от параметра вида напряженно-деформированного состояния построим решения уравнений в случае , и , вычисленного из (6). Найденные зависимости представлены на рис. 1. Из рисунка видно, что при локальная выпуклость потенциала имеет место при . В случае максимальное допустимое значение поврежденности определяется неравенством . После точки - условием . При этом вторая критическая величина следует из равенства и имеет вид . (7) 2. Нелинейная модель деформирования хрупкого тела с тензорным параметром поврежденности Следуя классической теории анизотропной поврежденности [44-45] зададим параметр поврежденности материала в виде симметричного тензора второго ранга Главные направления этого тензора задают направления ортотропии материала, а главные значения определяют изменение эффективной площади поперечных сечений, перпендикулярных к каждой из осей ортотропии (рис. 2, а): , (8) где и - площади поперечного сечения неповрежденного и поврежденного материала в плоскости с нормалью а б Рис. 2. Представительный объем материала в поврежденном (а) и фиктивном неповрежденном состоянии (б) Fig. 2. Representative volume of material in damaged (a) and fictitious undamaged state (б) Для перехода от поврежденного состояния (характеризующегося несимметричными тензорами напряжения и деформации) к фиктивному неповрежденному состоянию примем гипотезу эквивалентности упругой энергии этих состояний [45]. При этом предположим, что ориентация тензора в поврежденной конфигурации не совпадает с ориентацией тензора деформации в конфигурации фиктивной неповрежденной (см. рис 2). Переход от одного состояния к другому подразумевает переход от несимметричных мер напряженного и деформированного состояния к симметричным. Введем симметризованный тензор фиктивной деформации, представляющий собой разность между деформацией в фиктивном неповрежденном состоянии и эффективной деформацией в поврежденном состоянии в виде , (9) характеризующий деформацию поврежденного материала. Инварианты тензора фиктивной деформации определяются как , и в частном случае в главных осях тензора деформации имеют вид (10) Обобщим потенциал, предложенный в [46], и представим энергию упругой деформации анизотропного поврежденного материала в виде (11) где , - параметры Ламе исходного неповрежденного изотропного материала, , параметры, характеризующие деградацию упругих модулей при росте поврежденности, - дополнительный модуль, ответственный за нелинейность при накоплении повреждений. Первое слагаемое в (11) представляет собой энергию деформирования исходного неповрежденного изотропного материала, второе слагаемое - энергию деформации материала с анизотропией, индуцированной поврежденностью, а третье слагаемое отвечает за нелинейность и разномодульность деформационного отклика материала, зависимость его свойств от вида напряженно-деформированного состояния. В случае нулевой поврежденности выражение (11) сводится к классической гуковской энергии упругого изотропного тела. Частный вид потенциала (11) для случая соосных тензоров деформации и поврежденности ранее был предложен в работе [46]. Дифференцируя потенциал (11) по компонентам тензора деформации получим выражения для компонент тензора напряжений Коши (следуя [46], примем ): (12) где , Из соотношения (12) следует особенность, присущая всем анизотропным моделям. При несоосных тензорах поврежденности и деформации тензор напряжения будет не соосен им обоим. Например, в случае диагонального тензора поврежденности при сдвиговом деформировании представительного объема (только ненулевых недиагональных компонентах) тензор напряжений будет иметь как недиагональные компоненты, связанные с характером приложенной деформации, так и диагональные компоненты, связанные с дилатансией за счет роста поврежденности. Необходимо также отметить, что при поврежденный материал будет ортотропным с упругими модулями, зависящими от компонент тензора поврежденности, вида напряженно-деформированного состояния и величины приложенных главных деформаций. В случае соотношение (11) редуцируется к нелинейной реологической модели (1) с модулями, линейно зависимыми от скалярного параметра поврежденности . Различие в степенях и является спецификой анизотропной формулировки и обусловлено использованием гипотезы эквивалентности энергии. Как было сказано ранее, для строгой выпуклости потенциала необходимо и достаточно, чтобы матрица Гессе была положительно определенной. Матрица Гессе для потенциала (11) в главных осях тензора деформации имеет вид: . Для положительной определенности матрица Гессе необходимо и достаточно выполнение следующих условий [51]: (13) В случае неравенства, совпадают и редуцируются к условию с Система неравенств (13) является необходимым и достаточным условием строгой локальной выпуклости потенциала (11) и единственности решения статической задачи нелинейной теории упругости. Рассмотрим несколько частных случаев вида тензора поврежденности и найдем предельные значения его компонент из условий (13). Трансверсально-изотропная трещиноватая среда. Широко распространенным типом анизотропии горных пород является трансверсальная изотропия, связанная с процессами осадконакопления и уплотнения материала. Анизотропия может быть вызвана не только слоистым строением среды, но и определенно ориентированной трещиноватостью, появившейся в процессе ее формирования. Рассмотрим представительный объем материала с ансамблем одинаково ориентированных микротрещин (рис. 3, а) - представительный объема трансверсально-изотропной среды. Трещиноватость описывается тензором поврежденности с главными значениями . Параметр определяет степень анизотропии поврежденности материала. Ориентация трещиноватости определяется эллипсоидом поврежденности (рис. 3, б). Рис. 3. Представительный объем трансверсально-изотропного трещиноватого материала (а) и соответствующей вид эллипсоида поврежденности (б) Fig. 3. Representative volume of transversely isotropic fractured material (a) and the corresponding view of the damage ellipsoid (б) Конфигурация действующих на образец механических воздействий описывается тензором деформации с главными значениями . Рассмотрим случай, когда тензоры деформации и поврежденности соосны. Для анализа условий строгой локальной выпуклости потенциала при различных видах напряженно-деформированного состояния будем варьировать в диапазоне определяя из выражения для параметра вида напряженно-деформированного состояния, который с заданным шагом изменяется от (случай всестороннего сжатия) до (случай всестороннего растяжения). Количество шагов для деформации и параметра совпадают. На рис. 4 представлены предельные значения поврежденности для различной степени ее анизотропии от величины параметра, характеризующего тип напряженно-деформированного состояния . Рис. 4. Зависимости допустимых значений поврежденности определяемых условиями (13), от вида напряженно-деформированного состояния, для параметра анизотропии поврежденности q = 0,95 (а) и q = 0,8 (б) Fig. 4. Dependences of the permissible values of damage determined by conditions (13), on the type of stress-strain state, for the parameter of damage anisotropy q = 0,95 (a) and q = 0,8 (б) Из данных рис. 4 видим, что независимо от величины пары условий и и дают тождественные результаты. При максимальное допустимое значение поврежденности равно единице, так же, как и в модели со скалярной поврежденностью. Начиная с допустимое значение поврежденности определяется условием При этом для малых степеней анизотропии есть еще одна критическая величина, начиная с которой допустимое значение поврежденности контролируется условием . Величина зависит от степени анизотропии поврежденности, при этом при любой величине выполняется условие . На рис. 5 представлена зависимость критической величины от степени анизотропии материала q. Видно, что в целом анизотропия поврежденности приводит к упрочнению материала (точка смещается в область больших значений При этом при слабой анизотропии критическая величина находится в области преимущественного сжатия материала, тогда как при величина асимптотически приближается к нулю, т.е. потеря устойчивости материала будет наблюдаться в условии чистого сдвига. Рис. 5. Зависимости критической величины от степени анизотропии поврежденности Fig. 5. Dependences of the critical value on the degree of damage anisotropy Сравнение полученных допустимых значений поврежденности для различной степени анизотропии с изотропным случаем (скалярная поврежденность) показывает, что условия строгой локальной выпуклости (5) дают минорантную оценку только в области (см. рис. 6). В случае обстановки преимущественного растяжения () допустимые значения параметра поврежденности определяются из условия системы неравенств (13). Рис. 6. Зависимости допустимых значений поврежденности от параметра x для случаев скалярной и анизотропной поврежденности Fig. 6. Dependences of the admissible values of damage on the parameter x for cases of scalar and anisotropic damage Трансверсально-изотропная среда с наклонно-ориентированной трещиноватостью. Рассмотрим обобщение предыдущей постановки на случай несоосных тензора поврежденности и тензора деформации. Как и ранее, примем, что трещиноватость описывается тензором поврежденности с главными значениями Ее наклон задается вращением эллипсоида поврежденности вокруг оси на угол между осью максимальной главной деформации и главной осью, соответствующей максимальной главной поврежденности (рис. 7). Будем работать в системе координат главных осей тензора деформации. Компоненты тензора поврежденности в этой системе координат будут определяться из ортогонального преобразования , где - ортогональный тензор, поворачивающих систему координат главных осей поврежденности в систему координат главных осей деформации. В случае и тензор поврежденности помимо диагональных будет иметь одну ненулевую недиагональную компоненту Рис. 7. Представительный объем трансверсально-изотропного материала с наклонно-ориентированной трещиноватостью (а) и соответствующей вид эллипсоида поврежденности (б) Fig. 7. Representative volume of transversely isotropic material with obliquely oriented fracturing (a) and the corresponding view of the damage ellipsoid (б) Используя выражения для компонент матрицы Гессе и условия ее положительной определенности (13), построим допустимые предельные значения поврежденности для различных углов ориентации трещиноватости b, степени ее анизотропии q и вида напряженно-деформированного состояния x Полученные зависимости для трех различных значений параметра анизотропии представлены на рис. 8. Из рисунка видно, что и вариация угла наклона эллипсоида поврежденности b, и степень анизотропии приводят к существенному изменению характера зависимости допустимой предельной поврежденности от параметра вида напряженно-деформированного состояния поврежденного материала. Так, при условия (5) для случая изотропной поврежденности являются минорантными, тогда как при уменьшении параметра q допустимая кривая для углов и лежит ниже кривой для изотропной поврежденности. Необходимо отметить, что при вариации угла наклона трещиноватости в материале максимальная воз- Рис. 8. Зависимости допустимых значений поврежденности , определяемых условиями (13), от вида напряженно-деформированного состояния и углов наклона для параметров анизотропии поврежденности q = 0,9 (а), q = 0,6 (б) и q = 0,3 (в) Fig. 8. Dependences of the permissible values of damage , determined by conditions (13), on the type of stress-strain state, and slope angles for the parameters of damage anisotropy q = 0,9 (a), q = 0,6 (б) and q = 0,3 (в) можная поврежденность, обеспечивающая его устойчивое поведение, изменяется немонотонно с изменением типа напряженно-деформированного состояния. Эта немонотонность связана с ролью трещиноватости при деформировании материала. Проиллюстрируем это на примере рис. 8, б. Так, при (случай трещиноватости, ориентированной перпендикулярно направлению максимальной главной деформации) наиболее сильные ограничения (относительно )имеют место для > 0, т.е. для случая преимущественного растяжения образца. Это вызвано тем, что трещины, раскрываясь, обеспечивают дополнительную деформацию материала, которая может расти лавинообразно. В случае (случай трещиноватости, ориентированной перпендикулярно направлению минимальной главной деформации) мы видим два интервала с наименьшими значениями допустимой поврежденности. Интервал соответствует случаю пропорционального или непропорционального многоосного растяжения, существующая трещиноватость в этом случае также дает вклад в общую деформацию. Второй интервал соответствует ситуации одноосного сжатия со свободной от напряжений боковой поверхностью или одноосного сжатия с принудительным растяжением в двух других ортогональных направлениях. Вертикальная трещиноватость в этом случае ориентирована нормально к направлению минимальной и/или промежуточной главной деформации, что способствует ее раскрытию и дополнительному деформационному вкладу. Роль ориентации трещиноватости также можно проследить на зависимости критической величины от угла наклона для различных степеней анизотропии поврежденности (рис. 9). Так, для слабой анизотропии () наблюдается разупрочнение материала при увеличении угла ориентации трещиноватости. Разупрочнение заключается в уменьшении критической величины x** при стремлении угла наклона к . В случае существенной анизотропии поврежденности (q < 0.5) можно выделить три характерных интервала изменения величины x**. Первый интервал к характеризуется потерей устойчивости материала при деформировании в условиях чистого сдвига. Второй интервал соответствует случаю, когда материал теряет устойчивость в случаях преимущественного растяжения. При (третий интервал) деформирование материала становится неустойчивым в условиях преимущественного сжатия (). Важно отметить, что в случае наклона трещиноватости субнормально действию минимальной главной деформации (угол близок к ) и высокой степени анизотропии критическая величина будет меньше . Рис. 9. Зависимости критической величины x** от угла наклона β для различных степеней анизотропии поврежденности Fig. 9. Dependences of the critical value x** on the inclination angle β for different degrees of damage anisotropy Заключение В данной работе осуществлен вывод условий локальной строгой выпуклости потенциала в неклассической модели нелинейной упругости с тензорным (второй ранг) параметром поврежденности. Получены условия выпуклости в главных осях тензора деформации для общего случая несоосных тензора деформации и поврежденности. Найденные условия обеспечивают единственность решения статической задачи теории упругости. Для иллюстрации полученных условий выпуклости рассмотрены два частных случая вида тензора поврежденности: трансверсально-изотропная трещиноватая среда с соосными тензорами деформации и поврежденности, трансверсально-изотропная среда с наклонно-ориентированной трещиноватостью. Обе конфигурации помимо поврежденности характеризуются также параметром степени анизотропии поврежденности, показывающем, насколько отличаются площади поперечных сечений в трех ортогональных направлениях, соответствующих главным осям тензора поврежденности. Для обоих случаев установлены наиболее строгие неравенства, показана зависимость предельных значений поврежденности от параметра степени анизотропии. Для наклонно-ориентированной трещиноватости построены зависимости предельной допустимой поврежденности от степени анизотропии, угла наклона и вида напряженно-деформированного состояния. Показано, что ориентация поврежденности в совокупности при разной степени ее анизотропии может приводить как к упрочнению материала, так и разупрочнению относительно случая изотропной трещиноватости.

About the authors

I. A Panteleev

Institute of continuous media mechanics UB RAS

V. Lyakhovsky

Geological survey of Israel

References

  1. Капустянский С.М. Анизотропия геоматериалов // Итоги науки и техники. Механика деформируемого твердого тела. - 1987. - Т. 18. - С. 53-113.
  2. Ставрогин А.Н., Протосеня А.Г. Прочность горных пород и устойчивость выработок на больших глубинах. - М.: Недра, 1985. - 271 с.
  3. Experimental investigation of the anisotropic mechanical behavior of phyllite under triaxial compression / Guowen Xu [et al.] // International Journal of Rock Mechanics and Mining Sciences. - 2018. - Vol. 104. - P. 100-112. doi: 10.1016/j.ijrmms.2018.02.017
  4. Динамическое деформирование мягких грунтовых сред: экспериментальные исследования и математическое моделирование / В.В. Баландин [и др.] // МТТ. - 2015. - № 3. - С. 69-77.
  5. Lockner D.A., Stanchits S.A. Undrained poroelastic response of sandstones to deviatoric stress change //j. Geophys. Res. - 2002. - Vol. 107. - P. B12. doi: 10.1029/2001JB001460
  6. Basaran C., Nie S. An irreversible thermodynamics theory for damage mechanics of solids. // Int. J. Damage Mach. - 2004. - Vol. 13. - P. 205-223. doi: 10.1177/1056789504041058
  7. Nonlinear and nonequilibrium dynamics in geomaterials / TenCate J.A. [et al.] // Phys. Rev. Lett. - 2004. - Vol. 93, № 6. - P. 065501. doi: 10.1103/PhysRevLett.93.065501
  8. Nonequilibrium and nonlinear dynamics in Berea and Fontainebleau sandstones: Low-strain regime / D. Pasqualini [et al.] //j. Geophys. Res. - 2007. - Vol. 112. - P. B01204. doi: 10.1029/2006JB004264
  9. Mechanical properties and non-homogeneous deformation of open-cell nickel foams: application of the mechanics of cellular solids and of porous materials / X. Badiche [et al.] // Materials Science and Engineering. - 2000. - Vol. 289. - P. 276-288. doi: 10.1016/S0921-5093 (00) 00898-4
  10. Краснощеков П.И., Федотов А.Ф. Упругие модули изотропных порошковых и пористых материалов // Вестн. Самар. гос. техн. ун-та. Сер. Физ.-мат. науки. - 2006. - Вып. 43. - С. 81-87. doi: 10.14498/vsgtu455
  11. Капитонов А.М., Редькин В.Е. Физико-механические свойства композиционных материалов. Упругие свойства. - Красноярск: Изд-во Сиб. федер. ун-та, 2013. - 532 с.
  12. Березин А.В. Влияние повреждений на деформационные и прочностные характеристики твердых тел. - М.: Наука, 1990. - 135 с.
  13. Бессонов Д.Е., Зезин Ю.П., Ломакин Е.В. Разно-сопротивляемость зернистых композитов на основе ненасыщенных полиэфиров // Изв. Сарат. ун-та. Сер. Математика. Механика. Информатика. - 2009. - Т. 9, вып. 4, ч. 2. - С. 9-13.
  14. Lomakin E.V., Fedulov B.N. Nonlinear anisotropic elasticity for laminate composites // Meccanica. - 2015. - Vol. 50, iss. 6. - P. 1527-1535. doi: 10.1007/s11012-015-0104-5
  15. Johnson P.A., Jia X. Nonlinear dynamics, granular media and dynamic earthquake triggering. Nature. - 2005. - Vol. 473. - P. 871-874. doi: 10.1038/nature04015
  16. Wu C., Peng Z., Ben-Zion Y. Non-linearity and temporal changes of fault zone site response associated with strong ground motion // Geophys. J.Int. - 2009. - Vol. 176, № 1. - P. 265-278. doi: 10.1111/j.1365-246X.2008.04005.x
  17. Water-level oscillations caused by volumetric and deviatoric dynamic strains / E. Shalev [et al.] // Geophys. J.Int. - 2016. - Vol. 204. - P. 841-851. doi: 10.1093/gji/ggv483
  18. Imaging and monitoring temporal changes of shallow seismic velocities at the Garner Valley near Anza, California, in relation to the M7.2 2010 El Mayor-Cucapah earthquake / L. Qin [et al.] //j. Geophys. Res. - 2020. - Vol. 125, № 1. - P. e2019JB018070. doi: 10.1029/2019JB018070
  19. Амбарцумян С.А., Хачатрян А.А. Основные уравнения теории упругости для материалов, разносопротивляющихся растяжению и сжатию // Изв. АН СССР. МТТ. - 1966. - № 2. - С. 44-53.
  20. Jones R.M. Stress-strain relations for materials with different moduli in tension and compression // AIAA Journal. - 1977. - Vol. 15, № 1. - P. 16-23. doi: 10.2514/3.7297
  21. Lomakin E.V., Rabotnov Yu.N. A theory of elasticity for an isotropic body with different moduli in tension and compression // Mechanics of Solids. - 1978. - Vol. 13, № 6. - P. 25-30.
  22. Medri G. A nonlinear elastic model for isotropic materials with different behavior in tension and compression // Journal of Engineering Materials and Technology. - 1982. - Vol. 104, № 1. - P. 26-28. doi: 10.1115/1.3225031
  23. A new computational framework for materials with different mechanical responses in tension and compression and its applications / Z. Du [et al.] // International Journal of Solids and Structures. - 2016. - Vol. 100. - P. 54-73. doi: 10.1016/j.ijsolstr.2016.07.009
  24. Vijayakumar K., Ashoka J.G. A bilinear constitutive model for isotropic bimodulus materials // Journal of Engineering Materials and Technology. - 1990. - Vol. 112, № 3. - P. 372-379. doi: 10.1115/1.2903341
  25. Ломакин Е.В. Механика сред с зависящими от вида напряженного состояния свойствами // Физическая мезомеханика. - 2007. - Т. 10, № 5. - С. 41-52.
  26. Пахомов Б.М. Вариант модели изотропного разно-модульного материала // Вестн. Моск. гос. техн. ун-та им. Н.Э. Баумана. Сер.: Машиностроение. - 2017. - № 6, Вып. 117. - С. 35-48. doi: 10.18698/0236-3941-2017-6-35-48
  27. Мясников В.П., Топалэ В.И. Моделирование сейсмической анизотропии в литосфере как разномодульном упругом теле // Изв. АН СССР. Физика Земли. - 1987. - № 5. - С. 22-30.
  28. Ломакин Е.В. Определяющие соотношения деформационной теории для дилатирующих сред // Изв. АН СССР. Механика твердого тела. - 1991. - № 6. - С. 66-75.
  29. Мясников В.П., Олейников А.И. Основы механики гетерогенно-сопротивляющихся сред. - Владивосток: Дальнаука, 2007. - 172 с.
  30. Пирогов С.А. Модель динамического деформирования и разрушения разносопротивляющихся материалов при интенсивных ударных воздействиях // Проблемы прочности и пластичности. - 2012. - Вып. 74. - С. 40-48.
  31. Fedorenko A.N., Fedulov B.N., Lomakin E.V. Buckling problem of composite thin-walled structures with properties dependent on loading types // PNRPU Mechanics Bulletin. - 2019. - № 3. - P. 104-111. doi: 10.15593/perm.mech/2019.3.11
  32. Lyakhovsky V., Ben-Zion Y., Agnon A. Distributed damage, faulting, and friction //j. Geophys. Res. - 1997b. - Vol. 102, № 27. - P. 635-649.
  33. Non-linear elastic behavior of damaged rocks / V. Lyakhovsky [et al.] // Geophys. J.Int. - 1997a. - Vol. 130. - P. 157-166. doi: 10.1111/j.1365-246X.1997.tb00995.x
  34. Dudko O.V., Lapteva A.A., Ragozina V.E. Nonstationary 1D Dynamics Problems for Heteromodular Elasticity with Piecewise-Linear Approximation of Boundary Conditions // PNRPU Mechanics Bulletin. - 2019. - № 4. - P. 37-47. doi: 10.15593/perm.mech/2019.4.04
  35. Convergence analysis of a finite element method based on different moduli in tension and compression / X. He [et al.] // International Journal of Solids and Structures. - 2009. - Vol. 46, № 20. - P. 3734-3740. doi: 10.1016/j.ijsolstr.2009.07.003
  36. An efficient finite element formulation for nonlinear analysis of clustered tensegrity / L. Zhang [et al.] // Engineering Computations. - 2016. - Vol. 33, № 1. - P. 252-273. doi: 10.1108/EC-08-2014-0168
  37. Hamiel Y., Lyakhovsky V., Ben-Zion Y. The elastic strain energy of damaged solids with applications to nonlinear deformation of crystalline rocks // Pure Appl. Geophys. - 2011. - Vol. 168. - P. 2199-2210. doi: 10.1007/S00024-011-0265-7.
  38. Ляховский В.А., Мясников В.П. О поведении упругой среды с микронарушениями // Изв. АН СССР. Физика Земли. - 1984. - № 10. - С. 71-75.
  39. Lyakhovsky V., Myasnikov V.P. Acoustics of theologically non-linear solids // PEPIAM. - 1988. - Vol. 50. - P. 60-64.
  40. Ляховский В.А., Мясников В.П. Поведение вязкоупругой среды с микронарушениями при растяжении и сдвиге // Изв. АН СССР. Физика Земли. - 1985. - №. 4. - С. 28-35.
  41. Ben-Zion Y., Lyakhovsky V. Analysis of aftershocks in a lithospheric model with seismogenic zone governed by damage rheology // Geophys.Int. J. - 2006. - Vol. 165. - P. 197-210. doi: 10.1111/j.1365-246X.2006.02878.x
  42. Agnon A., Lyakhovsky V. Damage distribution and localization during dyke intrusion // The physics and chemistry of dykes / eds. G. Baer, A. Heimann. - Rotterdam: Balkema, 1995. - P. 65-78.
  43. Олейников А.И., Могильников Е.В. Единственность решения краевых задач и устойчивость для разномодульного нелинейного материала // Дальневосточный математический журнал. - 2002. - Т. 3, № 2. - С. 242-253.
  44. Voyiadjis G.Z., Kattan P.I., Yousef M.A. Some basic issues of isotropic and anisotropic continuum damage mechanics // Handbook of Damage Mechanics. - Berlin: Springer, 2015. - P. 3-42. doi: 10.1007/978-1-4614-8968-9_1-1
  45. Zhang W., Cai Y. Continuum Damage Mechanics and Numerical Applications. - Berlin: Springer, 2010. - 1000 p. doi: 10.1007/978-3-642-04708-4
  46. Non-linear anisotropic damage rheology model: theory and experimental verification / I. Panteleev [et al.] // European journal of mechanics: A. Solids. - 2021. - Vol. 85. - P. 104085. doi: 10.1016/j.euromechsol.2020.104085
  47. Ekland I., Temam R. Convex Analysis and Variational Problems. - New York: Elsevier, 1976. - 390 p.
  48. Bazant Z.P., Cedolin L. Stability of Structures. Elastic, Inelastic, Fracture and Damage Theories. - New York: Oxford Univ. Press. 1991 - 984 p.
  49. Rudnicki J.W., Rice J.R. Conditions for the localization of deformation in pressure-sensitive, dilatant materials //j. Mech. Phys. Solids. - 1975. - Vol. 23. - P. 371-394.
  50. Schreyer H.L., Neilsen M.K. Analytical and numerical tests for loss of material stability // lnt. J. Num. Methods Eng. - 1996. - Vol. 39. - P. 1721-1736.
  51. Курош А.Г. Курс высшей алгебры. - М.: Наука, 1965. - 431 с.

Statistics

Views

Abstract - 21

PDF (Russian) - 25

Cited-By


PlumX


Copyright (c) 2022 Panteleev I.A., Lyakhovsky V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies