Решение плоской задачи наследственной теории старения методом конечных элементов

Содержание

Слайд 2

Проблема В настоящее время накоплен значительный экспериментальный и теоретический материал по

Проблема
В настоящее время накоплен значительный экспериментальный и теоретический материал по ползучести

бетона. Вместе с тем известные подходы к расчету бетонных и железобетонных конструкций с учетом ползучести ориентированы главным образом на решение задач с относительно простой геометрией изделия и при неизменной внешней нагрузке. В литературе отсутствуют данные о конечноэлементном моделировании напряженно-деформированного состояния бетонных конструкций с учетом старения материала и истории нагружения.
В работе J.T. Boyle и J. Spence приведен численный алгоритм решения краевых задач ползучести (I и II стадии), основанный на использовании уравнения состояния , где σ – напряжение, t – время, Т – температура . Данный подход не учитывает изменение деформационных характеристик материала во времени, не описывает обратной ползучести при разгрузке и ориентирован на конструкции из металла и металлических сплавов.
Неявная схема численного интегрирования реализована ANSYSе.
Слайд 3

Цель исследований Разработка и программная реализация алгоритма метода конечных элементов для

Цель исследований
Разработка и программная реализация алгоритма метода конечных элементов для решения

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

Наследственные функции II рода для бетона Операторно-матричной форма физических соотношений ,

Наследственные функции II рода для бетона

Операторно-матричной форма физических соотношений

,

где векторы

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

,

,

- матрица модулей упругости материала (3×3),

- ядро релаксации (наследственная функция II рода),

- интегральный оператор,

τ – «возраст» материала в момент приложения нагрузки,

t – временная координата, отсчитываемая от момента времени τ.

Вид функции базируется на принятой механико-математической модели ползучести бетона и применяемой функции меры ползучести .

Слайд 5

Величину принимаем равной 2 сут. , (1) Теория старения Ядро релаксации

Величину принимаем равной 2 сут.

, (1)

Теория старения

Ядро релаксации

где функция меры

ползучести Н.Х. Арутюняна

,

закон изменения модуля упругости материала

константы материала, определяемые

из опытов на ползучесть при одноосной деформации

,

Достоинством представления наследственной функции в форме (1) является то, что после обработки этого выражения в среде символьного процессора системы Maple 12, получаем следующую, легко программируемую формулу:

Слайд 6

, где введены обозначения: (2)

,

где введены обозначения:

(2)

Слайд 7

Теория наследственного старения Функция меры ползучести С.В. Александровского . (3) Здесь

Теория наследственного старения

Функция меры ползучести С.В. Александровского

. (3)

Здесь обозначено:

;

- опытные

константы.

Выражение для ядра релаксации, полученное на основе функции (3) имеет вид:

, (4)

Слайд 8

После обработки выражения (4) символьным процессором системы Maple 12 и группировки

После обработки выражения (4) символьным процессором системы Maple 12 и группировки членов

получим следующую формулу для наследственной функции:

где

;

.

,

(5)

Слайд 9

где

где

Слайд 10

Графики функции для бетона различного «возраста» Теория старения (4) Теория наследственного

Графики функции для бетона различного «возраста»

Теория старения
(4)

Теория наследственного
старения

(5)

1 – τ = 2 сут; 2 - τ = 8 сут;
3 - τ = 14 сут; 4 - τ = 20 сут;
5 - τ = 30 сут; 6 - τ = 40 сут

1 – τ = 2 сут; 2 - τ = 4 сут;
3 - τ = 6 сут; 4 - τ = 10 сут;
5 - τ = 20 сут; 6 - τ = 30 сут

Константы материала:

Слайд 11

Конечноэлементная реализация плоской задачи теории ползучести Векторы-столбцы наследственных деформаций и перемещений

Конечноэлементная реализация плоской задачи
теории ползучести

Векторы-столбцы наследственных деформаций и перемещений

Отсюда

где

- матрица, образованная из функций формы.

В соответствии с принципом возможных перемещений имеем

,

где вариация энергии деформации

возможная работа внешних сил

.

(6)

Слайд 12

Здесь обозначено: - матрица, устанавливающая связь ; - объем, занимаемый КЭ;

Здесь обозначено: - матрица, устанавливающая связь

;

- объем, занимаемый КЭ;

поверхность КЭ, к

которой приложена
распределенная нагрузка;

- вектор-столбец объемной нагрузки;

- вектор-столбец распределенной нагрузки.

Из уравнения (6) следует матричное соотношение

,

где матрица жесткости КЭ

и вектор-столбец узловых сил КЭ .

Подставляя в выражение (7) получим

(7)

.

(8)

Слайд 13

Для вычисления интеграла воспользуемся численным методом, основанном на формуле трапеций. Разобьем

Для вычисления интеграла

воспользуемся численным методом, основанном на формуле трапеций. Разобьем

рассматриваемый временной интервал
на m равноотстоящих временных шагов , так чтобы
. Тогда, выражение (8) можно записать в форме

Или в компактном виде

,

где

В выражении (9) вектор-столбец соответствует упруго мгновенному решению задачи.

(9)

Слайд 14

Пример 1. Мгновенное нагружение и последующая ступенчатая разгрузка призматического бетонного образца

Пример 1. Мгновенное нагружение и последующая ступенчатая разгрузка призматического бетонного образца

с размерами 6×6×30 см (опыт А. Д. Росса ). “Возраст” бетона в момент нагружения образца τ = 28 сут, момент времени распалубки = 2сут. Конечноэлементная разбивка 2×8 (2 КЭ по ширине и 8 КЭ по высоте). Шаг интегрирования = 1 сут.

Числовые примеры

1 – теория упругой
наследственности;
2 – теория старения;
3 – теория наследственного старения

Слайд 15

Пример 1. Ползучесть однопролетной железобетонной балки, мгновенно нагруженной сосредоточенной силой посередине

Пример 1. Ползучесть однопролетной железобетонной балки, мгновенно нагруженной сосредоточенной силой посередине

пролета.

Диаметр арматуры 8 мм. Модуль упругости арматуры 2,1∙105 МПа. Разбивку балки на плоские КЭ выполняем сеткой 6×40 КЭ (6 КЭ по высоте и 40 КЭ по длине). Арматуру моделируем стержневыми КЭ (40 КЭ). Шаг интегрирования
=1 сут. Момент времени распалубки = 2сут.

Слайд 16

схема армирования I (арматура расположена в нижней части балки); схема армирования

схема армирования I
(арматура расположена в
нижней части балки);

схема армирования II
(симметричное
армирование)
1 – τ = 14 сут;
2 – τ = 28 сут

Теория упругой наследственности

Теория старения

Теория наследственного
старения

___

- - - -

Слайд 17

Выводы: В рамках теории старения и теории наследственного старения с помощью

Выводы:

В рамках теории старения и теории наследственного старения с помощью символьного

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