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

Размещено на http://www.allbest.ru/

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

Л.Н. Панасюк,

Ю.Ш. Чубка,

В.С. Тюрина,

А.У-Б. Пошев

Аннотация

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

Л.Н. Панасюк, Ю.Ш. Чубка, В.С. Тюрина, А.У-Б. Пошев,

Академия строительства и архитектуры Донского государственного технического университета, Ростов-на-Дону.

В статье рассмотрен вопрос ползучести в строительной механике. Записаны системы нелинейных разрешающих уравнений для разных видов ядер ползучести. Рассмотрен алгоритм решения задачи нелинейной наследственности.

Ключевые слова: ползучесть, нелинейная наследственность, аппроксимации ядер ползучести, метод конечных элементов, ядра экспоненциального типа, метод наименьших квадратов, узлы коллокации.

Содержание статьи

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

В строительной механике одной из наиболее общих теорий описания процессов ползучести является наследственная линейная теория ползучести, называемая также теорией линейной наследственности Больцмана - Вольтера, которая в случае одномерного напряженного состояния записывается в форме (1). В (1) е - деформация, у- напряжение, E - мгновенный модуль упругости, K(t-ф) называется ядром ползучести разностного типа [1, 2]. Практическое использование механических моделей для решения прикладных задач невозможно без достоверной методики определения параметров, входящих в определяющие уравнения. Для наследственных моделей ползучести эти параметры, как правило, входят в структуру ядер наследственности и определяются по результатам обработки простых и легко воспроизводимых экспериментов.

(1)

Для описания процессов нелинейной теории наследственности в одномерном случае, предложенной Ю.Н. Работновым [2-3], используется следующее уравнение состояния:

(2)

В (2) ц - нелинейная функция связи между напряжениями и деформациями, полученная при мгновенном эксперименте (рис. 1), а ядро К определяется из эксперимента на линейную ползучесть. Значительным достоинством теории Ю.Н. Работнова является с одной стороны, достаточно простой состав экспериментов для определения параметров сложных процессов нелинейной наследственности, с другой - все преимущества теории наследственности по отношению к простым техническим гипотезам ползучести.

Рисунок 1. Мгновенная диаграмма деформирования материала

Ю.Н. Работновым использовано представление нелинейного функционала в виде ряда кратных интегралов Фреше [4]. В этом представлении ядра задавались произведением одинаковых функций различных аргументов.

С учетом обозначения (3) получены соотношения (4):

(3)

(4)

Из обращения ряда (4) следует (5).

(5)

В (5) под символом L*б понимается свертка функций L и у, где L - ядро ползучести. Обозначения в (5): е - деформации, у - напряжения одномерной теории наследственности, L* - интегральный оператор (оператор Вольтерры). При t=0 соотношения (5) сводится к зависимости нелинейной теории упругости, где ц(е) задает "мгновенную" диаграмму деформирования материала. В случае линейной зависимости в "быстром" эксперименте при t=0 имеем . При отсутствии эффекта ползучести (L?0) получаем соотношения одномерной нелинейной теории упругости:

.

Перепишем (4) в виде (6).

где (6)

На основе подхода Ю.Н. Работнова получим альтернативную форму соотношений наследственной упругопластичности. Для этого выполним разложение напряжений в ряд Фреше (7).

(7)

В (7) также предполагается выполнение гипотезы Ю.Н. Работнова - (8).

(8)

Используя обозначения (9), получаем физические соотношения в виде (10). После обращения (10) получим (11).

(9)

(10)

(11)

Здесь R - ядро релаксации, - модифицированная деформация. Функция f(у) (обратная к ?()) соответствует мгновенной диаграмме деформирования. Действительно, при t=0 выполняются соотношения, описывающие одну и ту же диаграмму мгновенного деформирования образца: ползучесть деформация строительная уравнение

.

Выполним переход от (10) к (12).

,

(12)

Итак, (4-6) соответствуют прямой форме Ю.Н. Работнова, а (9-11) альтернативной закона наследственной одномерной нелинейной теории упругости [5-7]. Далее ядра L и R будем записывать с индексом "единица" (L₁, R₁) для прямой формы, и с индексом "два" (L₂, R₂) для альтернативной формы.

Альтернативная форма закона позволяет получить разрешающие уравнения для случая сложного НДС, а также конечно-элементную реализацию (13).

(13)

В (23) интегрирование (1-R*)Дq можно выполнить численно. Например, интегрирование методом трапеций дает:

С учетом последнего равенства, уравнения равновесия метода конечных элементов аппроксимируются как (14).

(14)

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

1. Вначале определяются мгновенные компоненты НДС, относящиеся к моменту t=0. Здесь на условном бесконечно-малом отрезке времени при t=0 вязкие свойства не проявляются, и выполняется решение стационарной задачи.

2. Второй этап выполняется циклически до конечного времени счета. В пределах очередного шага модифицированные деформации, напряжения, модули для текущего момента времени. Строятся локальные матрицы жесткости и формируется глобальная система уравнений МКЭ (14) [8]. Решение этой системы дает вектор приращения модифицированных перемещений.

3. После достижения конечного времени счета выполняется контроль точности шага интегрирования Дt. Для чего проводится сгущение сетки временных конечных элементов (т.е. уменьшение Дt) до стабилизации решения.

Рассмотрена возможность аппроксимации ядер ползучести следующими известными и часто используемыми функциями:

Простейшее ядро Больцмана:

Усложнение ядра Больцмана (ядро Абеля):

Ядро Абеля описывает следующую зависимость при постоянном уровне напряжений:

Ядро Ю.Н. Работнова:

,

где Г() - гамма-функция.

Ядро Г.М. Слонимского:

.

Частный случай этого ядра:

Приведенные выше ядра были сингулярными, имеющими особенность в нулевой точке по времени.

Часто на практике используются ядра экспоненциального типа, не являющиеся сингулярными. В общем виде:

.

В частном случае часто используется ядро с одним слагаемым:

.

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

Итак, согласно теории Ю.Н. Работнова, ядро ползучести можно строить на простом эксперименте на линейную ползучесть [9, 10]. Для этого уровень напряжений на образец должен быть сравнительно невысоким. Таким, при котором не будет происходить ощутимого развития нелинейных эффектов, в частности, накопления необратимых деформаций в образце. Определить, при каком уровне воздействий не проявляются нелинейные эффекты, можно на простом эксперименте на "мгновенное" нагружение образца.

Рассмотрим основную зависимость (1) и выполним эксперимент при постоянном уровне нагрузки:

.

Тогда зависимость (1) (с учетом коммутативности операции сверки:

=)

примет вид (15):

(15)

Или, с учетом того, что

- мгновенное значение деформации при t=0 после приложения нагрузки , имеем:

(16)

Продифференцировав (16) по времени, получим:

. (17)

Выражение (17) задает явную связь между ядром ползучести и скоростью ползучести: функция ядра ползучести пропорциональна скорости роста деформаций при ползучести с коэффициентом пропорциональности , где - мгновенное значение деформации при t=0 после приложения нагрузки . Символом н(t) в (17) обозначена скорость нарастания деформаций ползучести.

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

Для этой аппроксимации в работе используется метод наименьших квадратов. Скорость ползучести в точке по времени приближается конечно-разностным выражением, в котором учитывается, что величина шага по времени при выполнении эксперимента может изменяться (см. рис. 2).

Рисунок 2. Конечно-разностная аппроксимация скорости ползучести

Для построения конечно-разностного выражения для определения скорости ползучести в момент времени t_k участок ползучести на интервале [t_k-1, t_k+1] аппроксимирован квадратной параболой с локальным временем в нулевой точке t_k:

(18)

Несложно убедиться в выполнении граничных условий в (18):

. (19)

В (18) и в (19) t - локальное время с нулем в t_k.

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

(20)

Подставим в (20) локальное t_k=0 и получим выражение для скорости ползучести в глобальный момент t_k:

(21)

Для определения параметров ядер ползучести по результатам лабораторных экспериментов используем МНК (метод наименьших квадратов). Этот математический подход основан на минимизации суммы квадратов отклонений некоторых функций от искомых переменных. Используется при решении переопределенных систем уравнений, а также для поиска решения нелинейных систем уравнений. В рассматриваемой задаче МНК можно использовать для приближения точечных значений эксперимента некоторой функцией.

Подход МНК записывают следующим образом:

В нашем случае имеется n значений некоторой переменной yи соответствующих переменных x. Задача заключается в аппроксимации взаимосвязи между y и x некоторой функцией f(x_t, b), известной с точностью до определяемых в процессе решения параметров b.

При решении задачи минимизации находим стационарные точки целевой функции, для чего находятся частные производные по неизвестным параметрам b. Система разрешающих уравнений получается после приравнивания производных к нулю:

(22)

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

Запишем системы нелинейных разрешающих уравнений для разных видов ядер ползучести:

Аппроксимация ядра Абеля. Определяются параметры (C, б) для приближения функции:

1. Aleksandrov A.F., Potapov V.D. Osnovy teorii uprugosti I plastichnosti [Fundamentals of the theory of elasticity and plasticity]. M.: Vysshaya shkola, 1990. 400 p.

2. Rabotnov Ju.N. Mehanika deformiruemogo tverdogo tela. [Mechanics of deformable solid]. M.: Nauka, 1979. 650 p.

3. Rabotnov Ju.N. Jelementy nasledstvennoj mehaniki tverdyh tel. [Elements of hereditary mechanics of solids]. M.: Nauka, 1977. 383 p.

4. Volterra V. Theory of Functionals and of Integral and Integro-Differential Equations. Dover Phoenix Editions, 1959. 304 p.

5. Vasil'kov G.V., Panasjuk L.N., Rogachkin P.L. Variacionnye postanovki zadach nasledstvennoj teorii techenija s izotropnym uprochneniem.Vychislitel'naja mehanika i modelirovanie raboty konstrukcij i sooruzhenij. [Variation tasking hereditary flow theory with isotropic hardening. Computational mechanics and modeling of structures and buildings]. Rostov n.D.: RGAS, 1992. 6 p.

6. Rabotnov Ju.N. Polzuchest' jelementov konstrukcij. [Creep of Structural Elements]. M.: Nauka, 1968.752 p.

7. Vasil'kov G.V., Panasjuk L.N., Selim Sh.I. Deformacionnaja teorija plastichnosti nasledstvennogo tipa. [Deformation theory of hereditary type of plasticity]. Rostov n.D.: RISI, 1991. 18 p. Dep. V. VINITI N 3792-V91.

8. Batht K.J. Finite Element Procedures. K.J. Baths. New Jersey: Prentice Hall, 1996. pp. 10.12.

9. Ju.Sh. Chubka, V.S. Tyurina, L.N. Panasjuk. Inћenernyj vestnik Dona (Rus), 2016, №3. URL: ivdon.ru/ru/magazine/archive/n3y2016/3742.

10. V.S. Tyurina, Ju.Sh. Chubka, L.N. Panasjuk. Inћenernyj vestnik Dona (Rus), 2016, №3. URL: ivdon.ru/ru/magazine/archive/n3y2016/3733.

Размещено на Allbest.ru