Математическое моделирование и анализ контактного взаимодействия в ортотропных средах на основе сглаживания коэффициентов жесткости

На правах рукописи

НЕКРАСОВА Наталия Николаевна

Математическое моделирование и анализ

контактного взаимодействия В ОРТОТРОПНЫХ СРЕДАХ на основе сглаживания коэффициентов жесткости

Специальность: 05.13.18 – Математическое моделирование,

численные методы и

комплексы программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук

Воронеж – 2011

Работа выполнена в ГОУ ВПО «Воронежский государственный технический университет»

Научный руководитель доктор технических наук, профессор

Алейников Сергей Михайлович

Научный консультант доктор технических наук, профессор

Бурковский Виктор Леонидович

Официальные оппоненты: доктор физико-математических наук,

профессор

Агранович Юрий Яковлевич;

кандидат технических наук, доцент

Крючкова Ирина Николаевна

Ведущая организация ГОУ ВПО «Липецкий государственный

технический университет»

Защита состоится 16 июня 2011 г. в 1230 часов в конференц-зале на заседании диссертационного совета Д 212.037.01 ГОУ ВПО «Воронежский государственный технический университет» по адресу: 394026, г. Воронеж, Московский просп., 14.

С диссертацией можно ознакомиться в научно-технической библиотеке ГОУ ВПО «Воронежский государственный технический университет».

Автореферат разослан 16 мая 2011 года.

Ученый секретарь

диссертационного совета Барабанов В.Ф.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Работа выполнена в соответствии с одним из научных направлений ГОУ ВПО "Воронежский государственный технический университет" "Вычислительные системы и программно–аппаратные комплексы" (ГБ 2010.18).

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

В соответствии с данной целью в работе поставлены и решены следующие задачи:

анализ существующих подходов к численному моделированию пространственных конструкций на упругом основании;

разработка модели и алгоритма численного анализа влияния неоднородности сжимаемой толщи основания на напряженное состояние плиты;

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

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

разработка алгоритма численного анализа состояния плит переменной геометрии, взаимодействующих с упругими основаниями;

разработка программного обеспечения модели анализа контактного взаимодействия плит с упругим основанием.

Методы исследования. В работе использованы методы теории моделирования, методы математической физики, вычислительной математики, разделы теории сопротивления материалов, объекто-ориентированного программирования.

Результаты соответствуют следующим пунктам паспорта специальности: п. 3 ("Разработка, обоснование и тестирование эффективных вычислительных методов с применением современных компьютерных технологий"); п. 4 ("Реализация эффективных численных методов и алгоритмов в виде комплексов проблемно-ориентированных программ для проведения вычислительного эксперимента"); п. 5 ("Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента").

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

модель формализованного описания напряженно-деформированного состояния ортотропных плит, отличающаяся возможностью анализа объекта моделирования в условиях пространственно-неоднородных оснований;

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

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

аналитические контактные модели взаимодействия ортотропных плит переменной геометрии, обеспечивающие учет специфических особенностей упругих неклассических оснований;

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

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

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

Реализация и внедрение результатов работы. Теоретические результаты диссертационной работы внедрены в учебный процесс на кафедре "Автоматика и информатика в технических системах" ГОУ ВПО "Воронежский государственный технический университет" в рамках дисциплины "Моделирование систем".

Апробация работы. Основные положения диссертационной работы докладывались на следующих конференциях, семинарах и совещаниях: XVI Международной конференции по теории оболочек и пластин (Н.Новгород, 1994); 47-й научно-технической конференции по проблемам архитектуры и строительных наук (Воронеж, 1994); II Украинской научно-технической конференции по механике грунтов и фундаментостроению (Полтава, 1995); научных семинарах кафедры высшей математики ГОУ ВПО " Воронежский государственный архитектурно-строительный университет" (1996 – 2010); всероссийской конференции "Новые технологии в научных исследованиях, проектировании, управлении, производстве" (Воронеж, 2010); а также на научных семинарах кафедры "Автоматика и информатика в технических системах" ГОУ ВПО "Воронежский государственный технический университет" (2005–2010).

Публикации. Основные положения диссертационной работы опубликованы в 15 научных работах, в том числе 3 – в изданиях, рекомендованных ВАК РФ. В работах, опубликованных в соавторстве и приведенных в конце автореферата, лично соискателю принадлежат: формализованное описание напряженно-деформированного состояния ортотропных плит в условиях пространственно-неоднородных оснований [4,5,6,7,9,10], аналитическая модель контактного взаимодействия ортотропных плит с упругим основанием, отличающаяся учетом свойств материала плиты и смешанных граничных условий [8], алгоритм численного анализа, позволяющий осуществлять расчет напряженно-деформированного состояния плит на основе сглаживания границ между областями с различными жесткостными характеристиками [1,2,3,13,14], структура программного комплекса анализа напряженно-деформированного состояния плит переменной и постоянной жесткости [11,15].

Структура и объем работы. Диссертационная работа состоит из введения, 4 глав, списка литературы из 87 наименований, основных выводов и приложений. Основная часть работы изложена на 120 страницах, содержит 12 таблиц и 26 рисунков.

ОСНОВНОЕ содержание работы

Во введении обоснована актуальность темы и практическая значимость работы с указанием цели и задач исследования. Отмечена научная новизна и практическая значимость полученных результатов.

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

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

Проблема расчета плит переменной жесткости на упругих основаниях возникает в задачах оптимального проектирования фундаментов, ориентированных на выполнение требований технологического характера: снижения их материалоемкости, обеспечения надежности, долговечности и экономичности. Это также связано с проектированием плит, изменение жесткостей (толщин) которых обусловливается расположением устанавливаемого технологического оборудования или значительным изменением нагрузки на отдельных участках многосвязных плит. В соответствии с механическими свойствами, плиты переменной жесткости относятся к трем типам характеризующимися следующими особенностями: материал плиты однородный, но при этом изменяется ее толщина; толщина плиты постоянна, но ее жесткость изменяется в результате армирования, изменчивости физических свойств материала, частичных пластических деформаций или трещинообразования; одновременно изменяется толщина и жесткость. В любом из этих трех случаев эффективная жесткость может описываться ступенчатым или непрерывным законами. Наиболее широкие исследования в данной области выполнены для круглых пластин переменной толщины. Это работы А.Д. Коваленко, В.А. Гришина, Е.Б. Кореневой, А.Г. Юрьева, Н.А. Смоляго и др.

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

На основе проведенного анализа определены цель и основные задачи исследования.

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

(1)

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

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

, (2)

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

Граничные условия на контуре плиты имеют вид:

1. защемление: (3)
2. шарнирное опирание: (4)
3. свободный край: (5)

где n – внешняя нормаль к контуру плиты ; изгибающий и крутящий моменты; перерезывающая сила; производная по нормали; производная по дуге контура плиты.

Таким образом, математическая формулировка рассматриваемой задачи сводится к совместному решению дифференциального (1) и интегрального (2) уравнений с граничными условиями на краях контура плиты (3-5). Если выразить из уравнения (1) и подставить в (2), то придем к сложному интегро-дифференциальному уравнению:

.

Решение сформулированной пространственной контактной задачи при любой модели упругого основания представляет значительные трудности, а для прямоугольных плит точные решения к настоящему времени практически отсутствуют. В диссертационной работе решение интегрального уравнения (2) относительно контактных давлений считается известным и является одним из входящих данных для разрабатываемого программного модуля. Численно-аналитическое решение уравнение (2) заимствовано из работ Л.В. Канторовича, В.И. Крылова, Ш.Э. Мекиладзе, С.М. Алейникова.

Затем, записывая левую часть уравнения (1), а также уравнений (3-5) в конечных разностях и подставляя интенсивность реактивного давления в правую часть (1), определяем прогибы, изгибающие моменты, перерезывающие силы в плите. При этом используется конечно-разностная схема типа сквозного счета, когда граница раздела областей с различными жесткостными характеристиками не выделяется явно. Для жесткостей, входящих в выражения коэффициентов конечно-разностных уравнений в узлах сетки, реализован принцип сглаживания.

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

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

Уравнение (1) во внутреннем узле конечно-разностной сетки имеет вид:

, (6)

.

Рис. 1. Конечно-разностная сетка, нанесенная на плиту

Частные производные в каждом слагаемом левой части уравнения (6) представлены в развернутом виде, на основе известных конечно-разностных соотношений для производных второго порядка (рис. 2):

;

;

;

;

.

Таким образом, уравнение (6) в узле окончательно принимает следующий вид:

, (7)

где разностные коэффициенты определяются по следующим формулам:

;

;

;

;

;

;

;

;

;

; ; ; .

Здесь - отношение квадратов шагов конечно-разностной сетки в направлениях и , соответственно.

Для удовлетворения второго порядка точности, значения жесткостей в сеточных узлах в работе определены согласно трем типам аппроксимации: 1) как среднее арифметическое значение жесткостей в двух соседних узлах (в рассматриваемом узле сетки и узле, расположенного слева от данного узла); 2) значения жесткостей в узле определяются как полусумма величин, обратных к среднему арифметическому значений жесткостей в соседних узлах; 3) значения жесткостей в сеточных узлах берем равными среднему арифметическому значению жесткостей в полуцелых узлах (промежуточные узлы, расположенные в серединах отрезков между узлами основной сетки).

а) б)

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

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

, (8)

здесь [Н] - матрица дифференциального оператора относительно вектора прогибов {W}; {q} и {p} - векторы внешней нормальной нагрузки и реактивного давления, рассматриваемые в узлах конечно-разностной сетки.

Система (8) решается методом Гаусса с исключением главного элемента по строке. В практических расчетах для матрицы системы (8) нами оценивалось её число обусловленности, которое показывает меру близости результатов к вырожденности. Проведенная оценка обусловленности на основе многочисленных расчетов показала, что решение системы уравнений методом исключения Гаусса дает вполне приемлемые результаты в качестве первого приближения.

С целью улучшения полученного решения оно в дальнейшем уточнялось итерационным методом. Результаты проведенных расчетов показали, что после 3-5 итераций решение, как правило, сходится в пределах машинной точности .

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

Для оценки сходимости численных решений и иллюстрации возможностей предложенной методики в работе приведены примеры расчета изгиба плит. Например, рассмотрена квадратная железобетонная плита переменной толщины со стороной под действием равномерно распределенной нагрузки. Толщина плиты скачкообразно меняется относительно диагонали квадрата (рис. 3, а) и по окружности (рис. 3, б), центр которой совпадает с центром плиты. Объем плиты при любом законе изменения толщины остается постоянным и равен . Для учета свойств ортотропии плиты полагаем Характеристики основания: Проведен сравнительный анализ типов аппроксимации жесткостей в сеточных узлах для плит с различной границей раздела жесткостей (толщин), различными граничными условиями на контуре и упругими основаниями.

а) б)

Рис. 3. Области плиты со скачкообразным изменением жесткостей (толщины): а) относительно диагонали; б) по окружности с радиусом R=4 м

Для защемленной по контору плиты, контактирующей с упругим полупространством и скачкообразным изменением толщины плиты в два раза (рис. 3,а в области 1 , а в области 2 и на самой диагонали , максимальный прогиб увеличится по сравнению с плитой постоянной толщины почти в полтора раза. Как показывают расчеты из таблицы, для рассматриваемой функции изменения толщины все три типа осреднения приводят к качественно подобному распределению прогибов; количественное отличие не превышает 5%, что вполне удовлетворяет заданной точности расчетов.

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

Максимальные прогибы защемленной плиты

	Тип аппроксимации
	1	2	3
	0,2005	0,2148	0,2050

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

Рис. 4. Изолинии равных прогибов W, см защемленной плиты с круговой линией раздела функции толщины на упругом полупространстве (первый тип аппроксимации жесткостей): а) h1=1м, h2=1,5м; б) h1=2,5478м, h2=1м

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

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

Для исследования влияния на напряженно-деформированное состояние плиты переменной жесткости упругого основания проведены расчеты для шарнирно опертой изотропной плиты, лежащей на упругом клине . Граница раздела ступенчатого изменения толщин расположена параллельно оси OY, изменяется в пределах . Толщину той части плиты, которая расположена ближе к ребру клина обозначим за , а толщину оставшейся части - . Будем предполагать, что . На рис. 7, а построены изолинии прогибов для случая, когда .

Рис. 5. Изолинии равных прогибов W, см защемленной плиты с круговой линией раздела функции толщины на упругом полупространстве (второй тип аппроксимации жесткостей): а) h1=1м, h2=1,5м; б) h1=2,5478м, h2=1м

Рис. 6. Изолинии равных прогибов W, см защемленной плиты с круговой линией раздела функции толщины на упругом полупространстве (третий тип аппроксимации жесткостей): а) h1=1м, h2=1,5м; б) h1=2,5478м, h2=1м

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

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

Рис. 7. Изолинии прогибов W, см шарнирно опертой плиты с кусочно-постоянной функцией толщины на упругом клине при: а); б)

Рис. 8. Изолинии прогибов W, см ортотропной шарнирно опертой плиты а) упругое полупространство; , ; б) упругий клин; ,

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

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

Предельные значения прогибов, полученные в рамках вычислительного эксперимента, совпадают в пределах 1,5 – 3%, что свидетельствует о достаточном качестве результата, который обеспечивает разработанный программный комплекс.

Рис. 9. График сходимости значений прогибов, по базовой модели (ряд 1)

и по эталонной модели (ряд 2)

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. С позиции системной методологии осуществлен анализ существующих подходов к численному моделированию пространственных конструкций на упругом основании.

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

3. Разработаны модель и алгоритм численного анализа влияния неоднородности сжимаемой толщи основания на напряженное состояние плиты.

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

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

6. Разработан алгоритм численного анализа состояния плит переменной геометрии, взаимодействующих с упругими основаниями.

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

Основные результаты диссертации опубликованы в следующих работах:

Публикации в изданиях, рекомендованных ВАК РФ

1. Бурковский В.Л. Численная реализация модели цилиндрического изгиба ортотропных плит переменной жесткости на упругом основании/В.Л. Бурковский, Н.Н. Некрасова// Вестник Воронежского государственного технического университета. 2009. Т.5. №10. С. 101-107.

2. Бурковский В.Л. Математическое моделирование влияния механических свойств основания на изгиб ортотропных плит с переменным коэффициентом жесткости/ В.Л. Бурковский, Н.Н. Некрасова// Вестник Воронежского государственного технического университета. 2009. Т.5. №10. С. 117-122.

3. Бурковский В.Л. Анализ адекватности модели изгиба фундаментных плит на основе конечно-элементного моделирования/ В.Л. Бурковский, Н.Н. Некрасова, В.М. Флавианов// Вестник Воронежского государственного технического университета. 2010. Т.6. №9.С. 15-17.

Статьи и материалы конференций

4. Алейников С.М. Численный расчет ортотропных плит на упругом основании переменной толщины/С.М. Алейников, Н.Н. Некрасова // Расчет и проектирование оснований и фундаментов в сложных инженерно-геологических условиях: сб. науч. тр. Воронеж: ВИСИ, 1992. С. 22-31.

5. Алейников С.М. Сходимость численных решений при расчете плит на упругом полупространстве методом сеток/С.М. Алейников, М.Д. Гончаров, Н.Н. Некрасова// Современные методы статического и динамического расчета сооружений и конструкций: сб. науч. тр. Воронеж: ВИСИ, 1993. Вып. 2. С.61- 67.

6. Некрасова Н.Н. Численный метод расчета плит на упругом полупространстве/Н.Н. Некрасова// Материалы 47-й науч.-техн. конф. Воронеж: ВИСИ, 1994. С.18-21.

7. Алейников С.М. Пространственная контактная задача изгиба пластин на упругом неклассическом основании/С.М. Алейников, М.Д. Гончаров, Н.Н. Некрасова//Труды XVI международной конференции по теории оболочек и пластин. Н-Новгород: Нижегор. ун-т, 1994. Т.3. С. 8-13.

8. Алейников С.М. Расчет ортотропных фундаментных плит при смешанных граничных условиях/С.М. Алейников, М.Д. Гончаров, Н.Н. Некрасова // Труды Украинской научно-технической конференции по механике грунтов и фундаментостроению. Ч.2. Эффективные фундаменты, изготавливаемые с выемкой грунта. Полтава: Полтав. технич. универ-т, 1995. С.79-83.

9. Алейников С.М. Изгиб ортотропных фундаментных плит, расположенных на упругих неклассических основаниях/ С.М. Алейников, Н.Н. Некрасова//Известия ВУЗов: Строительство. 1996. № 9. С.55-60.

10. Aleynikov S.M., Nekrasova N.N. Contact problem for ortotropic faundation stabs with consideration on oleformation peculiarities of spatial and nonhvomogeneous bases //Studia Geotehnica et Mechanica, V.20, №1-2, p.63-104, 1998.

11. Алейников С.М. Программа "ORTOPLIT" расчёта ортотропных фундаментных плит переменной жесткости на упругих неклассических основаниях/С.М. Алейников, Н.Н. Некрасова // Компьютерные учебные программы и информации. М. 2002.

12. Алейников С.М. Изгиб ортотропных плит переменной жесткости на упругих неклассических основаниях/ С.М. Алейников, Н.Н. Некрасова// Тезисы доклада 61-й научно–технической конференции. Новосибирск, 2004. С. 10.

13. Бурковский В.Л. Математическое моделирование изгиба фундаментных плит переменной жесткости на упругом основании/ В.Л. Бурковский, Н.Н. Некрасова //Новые технологии в научных исследованиях, проектировании, управлении, производстве: сб. тез. всерос. конф. Воронеж: ВГТУ, 2010. С. 30.

14. Бурковский В.Л. Конечно-разностное моделирование пространственного контакта ортотропных плит переменной жесткости и упругого основания/ В.Л. Бурковский, Н.Н. Некрасова //Новые технологии в научных ис-

следованиях, проектировании, управлении, производстве: сб. тез. всерос. конф. Воронеж: ВГТУ, 2010. С. 9.

15. Бурковский В.Л., Некрасова Н.Н. Программный модуль «SIMORT» расчета ортотропных плит с переменным коэффициентом жесткости на упругих основаниях. ГОСФАП, № 5020110087 от 10.01. 2011 г.

Подписано в печать 13.05. 2011.

Формат 60x84/16. Бумага для множительных аппаратов.

Усл. печ. л. 1,0. Тираж 85 экз. Зак. №

ГОУ ВПО «Воронежский государственный технический университет»

394026 Воронеж, Московский просп., 14