Институт Государственного управления, Главный редактор -

)

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

5

33ТРГСУ412

Институт Государственного управления, Главный редактор -

)

5

33ТРГСУ412

Ростовский Государственный Строительный Университет

05.23.17 Строительная механика

Об алгоритмах корректировки математической модели слоистой конструкции на основе экспериментальных данных

Ляпин Александр Александрович, заведующий кафедрой информационных систем в cтроительстве

Мещеряков Иван Алексеевич, Аспирант

Аннотация

многослойный плоскопараллельный полупространство

Анализ моделей многослойных сред с плоскопараллельной слоистостью используются при оценке долговечности и прочности дорожных конструкций, аэродромных покрытий. В статье приводится расчетная модель многослойного полупространства, построенная на основе аналитико-численного подхода. Численный анализ выявил, что характеристикой, чувствительной к варьированию модулей упругости верхних слоев конструкции, является так называемая «чаша максимальных динамических прогибов». Эта же характеристика может быть получена и в результате эксперимента, что обуславливает необходимость корректировки математической модели. Для решения этой задачи использован метод наименьших квадратов с регуляризацией по А.Н. Тихонову. Исследование выявило недостатки этого метода, и была обоснована необходимость использования взвешенный метод наименьших квадратов.

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

Abstract

Of correction algorythms for the mathematical model of a multilayer structure based on experimental data

Lyapin Alexander Alexandrovich, Rostov State Univercity of Civil Engineering

Head of the department of information systems in construction;

Meshcheryakov Ivan Alexejevich, postgraduate student

The analysis of multilayer structures with flat-parallel foliation is used to estimate durability of pavement structures and airfield surfacings. The article contains calculation model for a multilayer structure based upon a analythical and numerical analysis. The so-called “maximum dynamic deflection cup” is the most susceptible characteristic for modifying density moduli of the structure's upper layers. This characteristic can be found by experimentation, thereby the algorythm for correcting of the mathematical model is needed. For this the version of leastsquares method with Tikhonov regularization is used. The studies had revealed certain defects of the mentioned method and substantiated the preferability of using the weighted least-squares method.

Keywords: Multilayer structures, correction of a mathematical model, regularization, weighted least-squares method

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

Рассматриваемая в работе область D представляет собой N -слойное упругое полупространство: D = D1?D2?...?DN, описываемое в цилиндрической системе координат (R,и,z) как (рисунок 1) :

- полупространство;

Упругие свойства сред в Dj , j = 01, ,...,N определяются плотностью с j , модулем упругости Ej и коэффициентом Пуассона н j .

Рис. 1. Область в цилиндрической системе координат

Движение среды в общем случае при отсутствии объемных сил определяется решением системы уравнений в частных производных [4] :

(1)

- перемещения в точке наблюдения r={R,и,z} в момент времени t .

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

Деформацию среды примем осесимметричной:

На границе области считаем заданными вектор-функции напряжений:

Система (1) с использованием закона Гука для линейно-упругого материала, связывающего компоненты тензора напряжений у( )j с компонентами тензора малых деформаций :

, (3)

может быть сведена к эквивалентной системе уравнений Ламе относительно функций перемещений точек среды:

0, (4)

, - приведенные частоты колебаний,

- скорости распространения продольных и поперечных волн в j - й среде. При наличии диссипации в среде, определяющейся вязкостью материала, коэффициенты Ламе становятся комплексными, что в свою очередь приводит к комплекснозначности приведенных частот. В этом случае для модуля упругости вязкоупругих слоев можно принять следующее соотношение [5]:

E? - длительный модуль упругости материала;

ф - температурно-зависимое время релаксации.

Так, например, для асфальтобетона можно принять: E? = 2.0e10Па, d =1.7 , k1=0.2, k2 =0.5, ф=0.1 с при 15 0 С.

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

На бесконечности потребуем выполнения условий излучения в форме принципа предельного поглощения [2].

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

Решив уравнения движения среды для установившегося состояния с временным множителем exp(?iщt), функцию u(r,t) получим в виде интеграла Фурье:

(5)

Здесь определяет импульс внешнего воздействия, P(щ) - его преобразование Фурье.

Разложим функцию внешней нагрузки по базису с наперед заданной степенью точности:

(6)

В качестве ортогональной системы базисных функций разложения выберем :

T - предельное время наблюдения за конструкцией, t* - время задержки [1].

Тогда для конечного набора стационарных решений U(r,щk )с использованием свойств преобразования Фурье легко показать выполнение:

(7)

для всех r из рассматриваемой области.

На основе численного анализа полученных соотношений выявлено, что характеристикой, чувствительной к варьированию модулей упругости верхних слоев конструкции, является изменение максимальной по времени амплитуды вертикального смещения max ; |u r;t | при удалении по поверхности конструкции от точки удара - «чаша максимальных динамических прогибов» [3]. Эта характеристика может быть получена и на основе проведения натурных измерений в результате установки на поверхности покрытия группы датчиков на разном удалении от точки удара. Изложенное определяет необходимость разработки численного алгоритма решения обратной задачи с целью корректировки расчетной модели. Принято решение использовать метод наименьших квадратов, как метод, обладающий большей численной устойчивостью.

Решение обратной задачи с применением метода наименьших квадратов с регуляризацией по А.Н. Тихонову приводится к следующей системе уравнений

Здесь U_T - расчетная чаша максимальных динамических прогибов (получаемая численно для заданных значений параметров), U_ЭКСП - заданная форма чаши (в общем случае - полученная в результате обработки экспериментальных данных).

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

где коэффициенты б>0 монотонно уменьшаются с ростом величины в i точке наблюдения, причем:

На рисунке 2 приведены результаты расчетов на модели слоистого полупространства чаши максимальных динамических прогибов поверхности конструкции, При расчетах величина динамического воздействия на поверхность конструкции определялась треугольным импульсом длительности 0.003 с и максимальным значением 2 МПа.

Институт Государственного управления, Главный редактор -

)

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

5

33ТРГСУ412

Институт Государственного управления, Главный редактор -

)

5

33ТРГСУ412

Рис. 2. Расчетная чаша максимальных динамических прогибов при изменении 0.25 ?R?1.5

Проектные параметры конструкции (при полном совпадении физических характеристик смежных слоев данные слои объединялись в один) приведены в таблице:

Таблица. Параметры дорожной одежды

В результате корректировки параметров модели по получаемым экспериментальным данным UЭКСП наиболее чувствительными оказывались параметры верхних слоев асфальтобетона и грунтового основания конструкции, диапазон изменения модулей упругости которых составлял 30-50 % от проектных.

Литература

1. Боев С.И., Селезнев М.Г. Об одном подходе в нестационарных задачах теории упругости // Изв. СКНЦ ВШ. Естеств. Науки. -1989. -№ 2. ?С.76?81.

2. Ворович И.И., Бабешко В.А. Динамические смешанные задачи теории упругости для неклассических областей. М.: Наука, 1989. -320 с.

3. Илиополов С.К., Селезнев М.Г., Углова Е.В. Динамика дорожных конструкций. - Ростов-на-Дону: Рост. Гос. Строит. Ун-т, 2002. - 258 с.

4. Новацкий В. Теория упругости. М.: Мир, 1975. -872 с.

5. H. Beckedahl, H. Hьrtgen, E. Straube Begleitende Forschung zur Einfьhrung des Falling Weigt Deflectometer (FWD) in der Bundesrepublik Deutschland // Forschung Straвenbau und Straвenverkehrstechnik, 1996. № 733.

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