ОЦЕНКА КОЭФФИЦИЕНТА ПОПЕРЕЧНОЙ ДЕФОРМАЦИИ ПРИ ОСЕСИММЕТРИЧНОМ ДЕФОРМИРОВАНИИ СПЛОШНЫХ СРЕД

Полный текст

Введение. Механические константы массивных тел, в частности грунтовых массивов - начальный модуль упругости, начальный модуль сдвига, начальный коэффициент поперечной деформации и так далее - определяются на основе экспериментов. В результате обработки экспериментальных данных строятся всевозможные математические модели, описывающие механическое поведение массивных тел. Эти модели представляют собой, вообще говоря, совокупность двух соотношений. Первое соотношение устанавливает связь между объёмными напряжениями и объёмными деформациями: . (1) Второе соотношение устанавливает связь между сдвиговыми напряжениями и сдвиговыми деформациями: . (2) Здесь первый инвариант тензора напряжений; первый инвариант тензора деформаций; интенсивность касательных напряжений (величина, пропорциональная второму инварианту девиатора напряжений); интенсивность деформаций сдвига (величина, пропорциональная второму инварианту девиатора деформаций). Эксперименты подтверждают факт взаимного влияния инвариантов тензоров и девиаторов напряжений и деформаций на механическое поведение грунтов и грунтовых массивов. Это означает, что объёмное деформирование грунта определяется как гидростатическим давлением, так и сдвиговыми напряжениями; а изменение формы грунтовых массивов происходит не только при действии касательных напряжений, но и при гидростатическом давлении. Таким образом, в формулах (1), (2) и модуль объёмного расширения (сжатия) , и модуль сдвига , в самом общем случае должны зависеть как от объёмной, так и от сдвиговой деформации, то есть . Следовательно, формулы (1) и (2) устанавливают перекрёстные зависимости между первыми инвариантами тензоров и вторыми инвариантами девиаторов напряжений и деформаций. Связь между переменными модулями объёмного расширения (сжатия) K(?,Г), сдвига G(?,Г) и переменными модулями упругости E(?,Г) и коэффициентом поперечной деформации устанавливается известными соотношениями: ; (3) (4) Вид функций и определяется видом кривых объёмного (1) и сдвигового (2) деформирования. Вид функции определяется видом кривой ?x = ?x(?x) при одномерном деформировании Что касается функции , то её вид может быть установлен лишь с использованием соотношения (33) для конкретного вида деформирования сплошной среды (плоское, сферически-симметричное, осе-симметричное и так далее), для конкретной математической модели сплошной среды и конкретных численных значений констант, входящих в описание математической модели. Деформационная теория пластичности сыпучей среды. Одной из наиболее удачных математических моделей, описывающих механическое поведение грунтов и грунтовых массивов, без учёта длительного действия нагрузки, является деформационная теория пластичности сыпучей среды, предложенная профессором Г.А. Гениевым в 1974 году [1]. Данная модель иллюстрирует формоизменение и состоит из параллельно соединённых неупругого элемента и элемента сухого трения. Она позволяет описать некоторые закономерности деформирования грунтовых сред, в частности, влияние среднего напряжения на вид зависимостей между вторыми инвариантами девиаторов напряжений и деформаций, некоторые особенности процесса разгрузки, реализует возможность непосредственного перехода от зависимостей напряжения-деформации к условию предельного равновесия, предусматривает физическую нелинейность диаграмм работы материала, а также эффект дилатансии в области напряжённых состояний, близких к предельным. Переменный модуль объёмного расширения (сжатия) и переменный модуль сдвига при этом имеют вид: (3) Здесь - начальный коэффициент объёмного расширения (сжатия); - начальный модуль сдвига при чистом сдвиге; - предельная интенсивность деформаций сдвига; f - аналог коэффициента внутреннего трения; q - коэффициент дилатансии. Модель, представленная соотношениями (3) базируется на гипотезе о малости деформаций, в основе которой лежит принцип неизменности начальных размеров (принцип «затвердевания»), а потому описывает механическое поведение грунтовых массивов, являющихся физически-нелинейными по своей природе, без учёта геометрической нелинейности. Вместе с тем при статическом деформировании грунтовых массивов деформации, вообще говоря, уже нельзя считать малыми и построение разрешающих уравнений следует выполнять уже «после» деформации следуя терминологии В.В. Новожилова [2]. Для построения деформационных соотношений сплошных сред с учётом геометрической нелинейности было предложено несколько подходов [3], в частности, подход, основанный на допущении о независимости формы записи замыкающих уравнений (1) и (2) от вида рассматриваемой "геометрии" сплошной среды (гипотеза об эквивалентности формы записи замыкающих уравнений). В соответствии с данной гипотезой, основные физические соотношения геометрически и физически нелинейного аналога сплошной среды, будут иметь вид: . (4) Здесь K*(?*,Г*) - геометрически нелинейный аналог переменного модуля объёмного расширения (сжатия); G*(?*,Г*) - геометрически нелинейный аналог переменного модуля сдвига. Конкретный вид этих величин определяет геометрически нелинейную математическую модель сплошной среды, в частности геометрически нелинейным аналогом деформационной теории пластичности сыпучей среды будет математическая модель, описываемая соотношениями: (5) Здесь первый инвариант тензора обобщённых напряжений; первый инвариант тензора нелинейных деформаций; интенсивность обобщённых касательных напряжений (величина, пропорциональная второму инварианту девиатора обобщённых напряжений); интенсивность нелинейных деформаций сдвига (величина, пропорциональная второму инварианту девиатора нелинейных деформаций). Связь между переменными модулями объёмного расширения (сжатия) K*(?*,Г*), сдвига G*(?*,Г*) и переменными модулями упругости E*(?*,Г*) и коэффициентом поперечной деформации ?*(?*,Г*) устанавливаются соотношениями: (7) Осесимметричное деформирование. Рассмотрим построение кривых ?(?,Г) и ?*(?*,Г*) для массивного тела, находящегося в условиях осе-симметричного деформирования и описываемого деформационной теорией пластичности сыпучей среды в геометрически линейной и геометрически нелинейной постановке. Для определения функций ?(?,Г) и ?*(?*,Г*) рассмотрим осесимметричное деформирования бесконечного сплошного массива за пределами цилиндрической полости начального радиуса r0. Пусть внутренняя поверхность цилиндрической полости получает радиальное перемещение . Определение напряжённо-деформированного состояния бесконечного сплошного массива с цилиндрической полостью, находящегося в состоянии осе-симметричного деформирования, математически сводится к решению двухточечной краевой задачи для обыкновенного дифференциального уравнения второго порядка с переменными коэффициентами, если краевые условия зададим следующим образом: (8) Первое из соотношений (8) означает, что задано радиальное перемещение точек цилиндрической поверхности; второе соотношение говорит о том, что на расстоянии от оси цилиндра радиальные перемещения равны нулю. В случае геометрически линейной модели сплошной среды разрешающее уравнение в перемещениях имеет следующий вид [4]: , (9) где ; ; причём При этом компоненты деформации связаны с радиальными перемещениями соотношениями: . Для случая геометрически нелинейной модели сплошной среды разрешающее уравнение в перемещениях записывается в следующей форме [4]: . (10) Здесь ; причём При этом связь между компонентами нелинейных деформаций и перемещениями устанавливается соотношениями: . Решение краевых задач (9), (8) и (10), (8) выполнено методом конечных разностей на кусочно-неравномерной сетке с погрешностью в остаточном члене порядка 0(h2) для всех узлов сетки, включая граничные [5]. Кусочно- неравномерная сетка задавалась следующим образом: область решения - интервал делится на равных участков, а затем каждый участок делится ещё на участков, где . Такое деление обусловлено тем, что с левого края области решения [r0 - rn] характеристики напряжённо-деформированного состояния (напряжения, деформации и их инварианты) имеют значительный градиент, приводящий к необходимости задавать очень мелкую сетку. Система конечно-разностных уравнений и для задачи (9), (8) (без учёта геометрической нелинейности), и для задачи (10), (8) (с учётом геометрической нелинейности) имеет следующий вид: (11) причём для задачи (9), (8) , а для задачи (10), (8) (i = 2,…,n - 1). Коэффициент равен отношению расстояния между соседними узлами кусочно-неравномерной сетки: . По построению . Построение решения выполнялось методом итераций, причём на каждой итерации система конечно-разностных уравнений решалась методом прогонки. Доказательство устойчивости метода прогонки для системы (11) приведено в работе [6]. В качестве начального приближения решения принималось решение данной краевой задачи для физически и геометрически линейной модели сплошной среды (в силу его интегрируемости). Окончание процесса итераций контролировалось по перемещениям: , где номер итерации; погрешность вычислений. Описанный алгоритм решения двухточечной краевой задачи реализован в пакете программ "OSA", разработанном в среде программирования математического пакета "MATLAB". При выполнении расчётов параметры геометрии полупространства, с приведённой начальной плотностью , принимались следующими: начальный радиус цилиндрической полости , толщина сжимаемой толщи полупространства , количество узлов кусочно-неравномерной сетки , количество расчётных точек . Механические константы сплошной среды имели следующие значения: ; предельная интенсивность деформаций сдвига ; приведённый коэффициент дилатансии ; коэффициент трения f = 0,1; 0,5; 0,9. Вычисления выполнялись с точностью . Начальное смещение границы цилиндрической полости принималось равным . Это значение обеспечивало допредельную стадию работы материала на границе цилиндрической полости для принятых констант сплошной среды и параметров геометрии полупространства. На рис. 1 представлено решение сформулированной краевой задачи, то есть, показаны графики относительных радиальных перемещений ur/Гs, относительных радиальных и тангенциальных деформаций и напряжений , относительных инвариантных величин . Все кривые построены для коэффициента внутреннего трения . Характер изменения коэффициента поперечной деформации по толщине цилиндрического полупространства для различных значений коэффициента внутреннего трения показан на рис. 2. На рис. 3 показан характер изменения коэффициента бокового давления ?? = ?r по толщине цилиндрического полупространства для различных значений коэффициента внутреннего трения f. На рисунках пунктирные линии соответствует случаю геометрически линейной модели, сплошные линии - геометрически нелинейной модели сплошной среды. Выводы. 1. Численные исследования показывают, что величина коэффициента поперечной деформации существенно зависит как от механических ха- рактеристик материала полупространства, уровня развития напряжённо-деформированного состояния в рассматриваемой точке, так и от вида рассматриваемой математической модели сплошной среды (геометрически линейной или геометрически нелинейной). Максимального значения коэффициент поперечной деформации достигает в точках сплошной среды, примыкающих к краю цилиндрической полости, то есть там, где уровень напряжённо- деформированного состояния является максимальным. Для значительных уровней напряжённо-деформированного состояния коэффициент поперечной деформации существенно превышает предельное значение для линейно-упругого материала, равное 0,5. Для малых уровней напряжённо-деформированного состояния коэффициент поперечной деформации не превышает значения 0,5. Для геометрически нелинейной модели коэффициент поперечной деформации выше, чем для геометрически линейной модели сплошной среды, причём, чем меньше коэффициент внутреннего трения, тем эта разница больше. 2. Что касается коэффициента бокового давления, то он монотонно изменяется по глубине сжимаемой толщи по кривой выпуклостью вниз. Его наибольшие значения фиксируются в конце сжимаемой толщи. Наименьшие значения коэффициент бокового давления принимает ближе к краю цилиндрической полости. Для геометрически нелинейной модели коэффициент бокового давления выше, чем для геометрически линейной модели сплошной среды, причём эта разница тем больше, чем меньше коэффициент внутреннего трения. Заключение. Представленные в статье результаты могут быть использованы при анализе несущей способности и деформативности сплошных сред типа грунтовых оснований под здания и сооружения, механическое поведение которых описывается геометрически и физически нелинейными моделями.

Об авторах

СЕРГЕЙ ВАСИЛЬЕВИЧ БАКУШЕВ

ФГБОУ ВО "Пензенский государственный университет архитектуры и строительства", каф. "Механика"

Автор, ответственный за переписку.
Email: bakuchsv@mail.ru

БАКУШЕВ СЕРГЕЙ ВАСИЛЬЕВИЧ - доктор технических наук, профессор, профессор кафедры «Механика», ФГБОУ ВО «Пензенский государственный университет архитектуры и строительства». Научные интересы: геометрически и физически нелинейная механика сплошной среды.

440071, г. Пенза, ул. Лядова, дом 30, кв. 84

Список литературы