На правах рекламы:
ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Авторитетность издания

ВАК - К1
RSCI, ядро РИНЦ

Добавить в закладки

Следующий номер на сайте

4
Ожидается:
09 Декабря 2024

Алгоритм и программа расчета напряженно-деформированного состояния песчаных грунтов при циклическом нагружении

The algorithm and program for calculating the sands stress-strain state under cyclic loading
Статья опубликована в выпуске журнала № 2 за 2013 год. [ на стр. 254-257 ]
Аннотация:На основе современной упруго-пластической модели Дафалиаса–Манзари, неявного метода Эйлера и итерационной схемы Ньютона–Рафсона разработана численная модель расчета деформаций несвязного грунта. Выполнена численная реализация разработанной модели на ЭВМ. Проведено моделирование поведения песчаных грунтов при монотонном и циклическом нагружениях. Сравнение результатов расчета с экспериментальными данными, полученными для песка различной плотности, показало, что численная модель достоверно описывает поведение грунта как при монотонном, так и при циклическом нагружении. Результаты теоретического прогноза напряженно- деформированного состояния песчаных грунтов при монотонном и циклическом нагружении показали эффектив- ность вычислительного алгоритма и позволяют рекомендовать его для конечно-элементных расчетов при решении геотехнических задач.
Abstract:A numerical model of non-cohesive soil deformation calculating is developed On the basis of modern elasticplastic Dafalias-Manzari model, implicit Euler method and the iterative Newton-Raphson scheme. The numerical implementation of the model on the computer is executed. The simulation of the sand behavior at a monotonic and cyclic loading is carried out. The calculation data and experimental information for sand with various compactness were compared. The computational model definitely describes sand behavior with monotonic and cyclic loading. The results of the theoretical prediction of the sands stress-strain state under monotonic and cyclic loading show the efficiency of the computational algorithm. They allow recommending this prediction for finite-element analysis for geotechnical problems solving.
Авторы: Миронов В.А. () - Тверской государственный технический университет, г. Тверь, Россия, доктор технических наук, Софьин О.Е. (sofjin@mail.ru) - Тверской государственный технический университет (доцент ), г. Тверь, Россия, кандидат технических наук
Ключевые слова: циклическое нагружение., моделирование, упруго-пластичность, несвязные грунты
Keywords: cyclic loading, modeling, elasto-plasticity, non-cohesive soils
Количество просмотров: 17158
Версия для печати
Выпуск в формате PDF (7.68Мб)
Скачать обложку в формате PDF (1.35Мб)

Размер шрифта:       Шрифт:

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

Упруго-пластические модели, как правило, формулируются в скоростях, используя законы теории пластического течения. Численное интегрирование этих моделей имеет решающее значение для успешного моделирования геотехни- ческих задач. Схема, показывающая место под- программы интегрирования упруго-пластической модели неявным методом Эйлера внутри программы метода конечных элементов, приведена на рисунке 1. Основной целью подпрограммы является вычисление приращений напряжений и параметров упрочнений грунта Ds и Dq. Здесь используются обозначения D=n+1–n, где n+1 и n соответствуют значениям функции, вычисленным в моменты времени tn+1 и tn. Кроме того, подпрограмма вычисляет соответствующие касательные модули грунта , которые будут применяться в последующих итерациях. При использовании касательных модулей нелинейный конечно-элементный алгоритм имеет асимптотическую скорость сходимости, что является одним из ключевых элементов для эффективного геотехнического анализа.

Для моделирования механического поведения песчаных грунтов наиболее популярной является двухповерхностная модель Дафалиаса и Манзари [1]. Она разработана на основе концепции критических состояний и позволяет описать наблюдаемые в опыте эффекты упрочнения и дилатансии плотных и рыхлых песков.

Согласно модели Дафалиаса–Манзари коэффициент пористости в критическом состоянии грунта определяется уравнением

 где ec, ref – начальное значение критического коэффициента пористости; p – давление; pat – атмосферное давление; lс, x – константы грунта.

Скорость полной деформации упруго-пластического грунта складывается из упругой и пластической частей:  где  – соответственно тензоры скоростей полной, упругой и пластической деформации.

Упругие свойства грунта описываются соотношениями Гука:  где  – тензор-девиатор скоростей деформации и скорость объемной деформации;  – тензор-девиатор скоростей напряжений и скорость среднего напряжения; G, K – модули сдвига и объемной деформации.

Модули G и K определяются выражениями:

 

где G0 – константа грунта; n – коэффициент Пуассона.

Поверхность текучести задается в форме конуса вращения с уравнением  где ; m – угловой коэффициент, принятый константой грунта; a – безразмерный тензор-девиатор, определяющий положение оси конуса в пространстве главных напряжений.

Принимается неассоциированный закон пластического течения. Нормаль к поверхности текучести  определяется выражением  где  – девиаторная часть тензора, определяющего направление нормали к поверхности текучести в p-плоскости;  – тензор второго ранга; dij – символ Кронекера.

Для учета влияния третьего инварианта тензора напряженного состояния используется функция  где  – отношение угловых коэффициентов при растяжении и сжатии; q – угол вида напряженного состояния Лоде.

Линия, параллельная нормали n, пересекает предельную, критическую и дилатансионную поверхности в трех точках (рис. 2) c соответствующими тензорами кинематического упрочнения грунта:

где  na – константы грунта;  – параметр упрочнения.

Скорости пластической деформации определяются тензором

 

где

Коэффициент дилатансии грунта

 

Изменение Ad зависит от изменения тензора структурной дилатансии z:

 где А0 – константа грунта.

Закон кинематического упрочнения принимается в виде  с функцией  где ain – начальные значения тензора a при инициировании нового процесса нагружения;  где h0, ch – константы грунта.

Закон изменения тензора структурной дилатансии записывается в виде  где cz, zmax – константы грунта.

Численная реализация упруго-пластической модели была выполнена на основе неявного метода Эйлера и итерационной схемы Ньютона–Рафсо­на [2, 3].

Алгоритм получения нелинейного решения заключается в следующем:

1. Задаются начальные значения напряжений и параметров упрочнения:

2. Проверяются значения функции текучести и невязок:  

Если  и , то итерации заканчиваются, иначе – переход к п. 3.

3. Вычисляются приращения напряжений и параметров упрочнения:

4. Вычисляются новые значения напряжений и параметров упрочнения:

  

Далее изменяется номер итерации k←k+1 и выполняется переход к п. 2.

В каждой итерации k уравнения невязки линеаризуются, чтобы найти более близкое решение с  При последующей итерации k+1 значения  находятся путем перехода от  к  с использованием упругих модулей , законов пластического течения и упрочнения. Так как поверхность текучести выпуклая, применяемая итерационная схема (рис. 3) гарантирует сходимость решения.

На основе разработанного алгоритма составлена вычислительная программа для ЭВМ на языке Visual Fortran. Программа апробирована на результатах известных опытов по трехосному сжатию изотропно-уплотненных образцов песчаного грунта в стабилометре [4].

Характеристики грунта, при которых выполнялись расчеты, представлены в таблице, где они разделены на группы в зависимости от описываемых ими свойств [5].

Константы модели

Свойство

Характеристики

Упругость

G0

125

n

0,05

Критическое состояние

M

1,25

c

0,712

ec,ref

0,934

l

0,019

x

0,7

Дилатансия

nd

2,1

A0

0,704

nb

1,25

Кинематическое упрочнение

h0

7,05

ch

0,968

Структурная дилатансия

zmax

2,0

cz

600

На рисунке 4 приведены расчетные кривые, которые сопоставлены с экспериментальными данными, полученными для песка различной плотности при монотонном нагружении с последующей разгрузкой. Как видим, результаты расчета достоверно описывают поведение грунта.

На рисунке 5 представлен теоретический прогноз деформации плотного песка при циклическом нагружении с постоянным p. Очевидно, что поведение со значительным разрыхлением будет наблюдаться и для реальных образцов плотного грунта.

Полученные результаты теоретического прогноза напряженно-деформированного состояния песчаных грунтов при монотонном и циклическом нагружении показали эффективность вычислительного алгоритма и позволяют рекомендовать его для конечно-элементных расчетов в геотехническом проектировании.

Литература

1.     Dafalias Y.F., Manzari M.T., Journ. of Engineering Mechanics, 2004, no. 130 (6), pp. 622–634.

2.     Zienkiewicz O.C., Taylor R.L., The Finite Element Method, Solid Mechanics, Butterworth-Heinemann, 2000, Vol. 24, 59 p.

3.     Миронов В.А., Софьин О.Е. Упруговязкопластическое деформирование анизотропных оснований // Вестн. МГСУ, 2011. № 7. С. 570–576.

4.     Verdugo R., Ishihara K., Soils and Foundations, 1996, no. 36 (2), pp. 81–91.

5.     Taiebat M., Jeremic B., Dafalias Y.F., Kaynia A.M., Cheng Z., Soil Dynamics and Earthquake Engineering, 2010, no. 30, pp. 236–257.

References

1.  Dafalias Y.F., Manzari M.T.,  Journ. of Engineering Me-chanics, 2004, no. 130 (6), pp. 622–634.

2.  Zienkiewicz O.C., Taylor R.L.,  The Finite Element Meth-od,  Solid Mechanics, Butterworth-Heinemann, Vol.  24, 2000,
459 p.

3.  Mironov  V.A.,  Sofin  O.E.,  Vestnik MGSU  [Proc. of  Mos-cow State Construction Univ.], no. 7, 2011, pp. 570–576.

4.  Verdugo R., Ishihara K., The steady state of sandy soils. Soils and Foundations, 1996, no. 36(2), pp. 81 –91.

5.  Taiebat M., Jeremic  B., Dafalias Y.F., Kaynia A.M., Cheng Z.,  Soil Dynamics and Earthquake Engineering, 2010, no. 30, pp. 236–257.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=3505&lang=&lang=&like=1
Версия для печати
Выпуск в формате PDF (7.68Мб)
Скачать обложку в формате PDF (1.35Мб)
Статья опубликована в выпуске журнала № 2 за 2013 год. [ на стр. 254-257 ]

Возможно, Вас заинтересуют следующие статьи схожих тематик: