CC BY
54
References:
1. Balamirzoev A.G. Prediction of deformation bases waterworks on saline soils. A.G.Balamirzoev.- Rostov-na-Donu.: Izd-vo SKNC VSH, 2001.- 201p.
2. Balamirzoev A.G. Development of the theory and methods of suffusion deformations prediction when filtering reschinovatyh grounds waterworks. dissertation doctor of Technical Sciences (05.23.07): 25.05.06: utv.13.10.06 Balamirzoev Abdul Gadzhibalaevich.- Makhachkala, 2006.-409p.
3. Mamedov K. M., Ismailov F. M. A method for solving the convective diffusion equations for the case of filtration in porous media of limited capacity. Scientists universities notes MViSSO AzerbSSR Ser. X, 1973, pp. 41- 47.
4. The hydrodynamic and physico-chemical properties of rocks. Edited by N. N. Verigina. - M. Nedra, 1977- 271p.
УДК 539,41: 629.7.023
Булгаков А.И., Манапов Р.М., Муселемов Х.М., Устарханов О.М.
РАСЧЕТ ТРЕХСЛОЙНЫХ БАЛОК ПО МЕТОДУ КОНЕЧНЫХ ЭЛЕМЕН-ТОВС ПОМОЩЬЮ ПК ЛИРА
Bulgakov A. I., Manapov R. M., Muselemov Kh.M., Ustarkhanov O.M.
CALCULATION OF THREE-LAYER BEAMS BY THE FINITE ELEMENT-METHOD WITH THE HELP OF LIRA PC
Аннотация. В настоящее время широкой популярностью у инженеров и научных работников пользуются программные комплексы, в основе которых заложен метод конечных элементов (МКЭ). Основной концепцией МКЭ является непосредственная дискретизация рассчитываемой системы, которая расчленяется расчетной сеткой на конечные элементы. В статье рассмотрены вопросы сходимости метода конечных элементов МКЭ для пластинчатых элементов.
На основании сделанных предпосылок выполнен численный эксперимент расчета напряженно-деформированного состояния трехслойной балки и произведено сравнение результатов с теоретическими данными, полученными по результатам расчета трехслойных балок (ТБ) проведенными по методу изложенном в справочнике Кобелева В.Н.
Ключевые слова: метод конечных элементов, сходимость, трехслойные балки, узловые переменные, линейная теория упругости, программные комплексы.
Abstract. Currently very popular among engineers and scientists are the software systems, which are based on the finite element method (FEM).
The basic concept of the finite element method is a direct discretization of the calculated system that is divided by the calculated net into finite elements.
In the paper are considered the questions of convergence of the finite element method for plate elements.
Based on the premises made is carried out the numerical calculation experiment of the stress-strained state of three-layer beams and compared the results with theoretical data obtained from the results of the calculation of three-layer beam (TB) by the method described in the Kobelev V.N reference book
Key words: finite element method, convergence, three-layer beams, the nodal variables, the linear elasticity theory, software systems.
Введение. Трехслойные балки представляют собой сложную пространственную структуру, образованную верхним и нижним несущими слоями, а также сотовой конструкцией заполнителя. Соединение несущих слоев с сотовым заполнителем осуществляется с помощью клея (рис.1).
а)
б)
Рисунок 1 - Трехслойная балка а) общий вид; б) сотовый заполнитель
Конструкция трехслойной балки сама по себе является весьма жесткой при работе на плоский изгиб, что делает ее весьма привлекательной для применения в реальных конструкциях, особенно учитывая ее относительную легкость [2]. В настоящее время расчет трехслойных конструкций (ТК) ведется различ-
ными методами. Одним из наиболее распространенным является метод, изложенный в книге [1] под редакцией Кобелева В.Н.
Развитие компьютерных технологий привело к созданию большого количества программ расчета прочности и устойчивости конструкций, основанных на различных методах строительной механики и теории упругости которое можно использовать для расчета напряженно-деформированного состояния (НДС) трехслойных конструкций.
Революцию в строительной механике и теории упругости произвело появление метода конечных элементов (МКЭ). Неограниченные возможности реализации компьютерных моделей на основе МКЭ обусловили разработку новых методов математической физики для решения задач динамики, устойчивости, физической и геометрической нелинейности.
В настоящее время широкой популярностью у инженеров и научных работников пользуются программные комплексы, в основе которых заложен метод конечных элементов (МКЭ), однако лишь единицы из них ориентированы на строительную область [3]. Так в России и странах ближнего зарубежья наибольшей популярностью пользуются программные комплексы STARK и MicroFe (Россия), SCAD и LIRA (Украина) основанные на нормах и стандартах данного региона.
Постановка задачи. Исходя из вышеизложенных предпосылок, авторами сделан выбор в пользу программного комплекса LIRA (ПК ЛИРА).
Рассмотрим основные положения МКЭ реализованные в ПК ЛИРА. Выбор расчетной схемы для МКЭ в перемещениях легко поддается алгоритмизации, что делает этот метод привлекательным с точки зрения применения современных компьютерных технологий. В основе МКЭ лежит линейная теория упругости, уравнения равновесия для которой могут быть записаны в следующем виде:
Au = - BT(DBu) + f = 0, (1)
где B - матричный линейный дифференциальный оператор, с помощью которого вектор деформации s(u) выражается через вектор перемещений u, s(u) = Bu; D - матрица упругости, определяющая линейную связь между векторами напряжений Г и деформаций s, г = Ds (закон Гука); BT - транспонированная матрица.
Потенциальную энергию деформации обозначим
П (и) =1J (DBu)(Bu)dn =1J r(u)s(u)dn (2)
2 n 2 n
Работу внешних сил выразим формулой
W(u) =J f • udn, (3)
n
а полную энергию
J (и) = П (и) + №(и) (4)
Из формулы вытекает, что функция и является решением дифференциального уравнения (1). Функционал полной энергии вида (4) представляет собой функционал Лагранжа, условием минимума которого являются равенства принципа возможных перемещений - при любом возможном перемещении V сумма возможных работ внутренних и внешних сил равна нулю. Обозначим:
а(иу) = / фВи)(Ву^П = ¡о-(и)е(у^О - возможная работа внутренних
о о
сил;
(/, у) = №у) = 11 • уо - возможная работа внешних сил.
о
Тогда принцип возможных перемещений запишется в виде: а(иу) + (/у) = 0. (5)
Левая часть равенства (5) является производной функционала полной энергии (4). Задача решения уравнений равновесия (1) и вариационная задача минимизации функционала (4), т.е. нахождения функции и, удовлетворяющей (5), эквивалентны, т.е. имеют одно и то же решение.
Вариационная постановка задачи имеет определенные преимущества, которые вытекают из того, что порядок производных в (4) понижается в 2 раза. Это позволяет более удобно формулировать граничные условия, смягчить требования к базисным функциям и более просто представить разностные выражения.
Методы исследования. Метод конечных элементов относится к вариационным методам, т.к. позволяет отыскивать минимум функционала (4) на основе уравнений (5). Основной концепцией МКЭ является непосредственная дискретизация рассчитываемой системы, которая расчленяется расчетной сеткой на конечные элементы. Искомая функция перемещений по области системы
и( х), ХЕ О приближенно может быть принята в виде:
ь
ии(х) = 14 <Р1(х) (6)
I=1
где Ь - общее число узловых неизвестных.
Следует помнить, что в каждом узле может быть различное число неизвестных, т.к. на некоторые из них могут быть наложены связи.
128
Вестник Дагестанского государственного технического университета. Технические науки. №2 (41), 2016
На основе подстановки (6) в (5) задача определения непрерывной функции и (х) сводится к определению значений конечного числа узловых неизвестных ц, которые находят из системы уравнений:
а
аЪ
I (щ ) = (П (щ )-Ж (щ ))
Г\ ь
I ад Я}у} -2 ]=1
а
аЧг
а(иь, ип)- (/, иъ)
а
аЧ7
\\
V
ь
/, I Ч]ин ]=1 л
(7)
ь
= а(ин'П) +(/'П) = IЧ] != о
]=1 о о
при
I=\х...,ь
Решая систему (7) будем полагать, что и (х) удовлетворяет главным
граничным условиям. По найденным из (7) значениям на основе (6) определяется функция перемещений по области системы, а по ней на основе известных соотношений теории упругости и другие компоненты напряженно-деформированного состояния.
Обозначим:
= /(ПВф] )ВъОО=а(ф] ^ ); р = //ф1сШ=(/,щ ) (8) о о
Матрицу К с элементами К у называют матрицей жесткости или матрицей системы уравнений МКЭ, вектор Р с элементами Р - вектором нагрузок и вектором правых частей.
Обозначим через ц вектор узловых неизвестных, тогда уравнение (7) в матричном виде будет:
Кч + Р=0 (9)
Расчленение системы на конечные элементы, выполненные на первом этапе расчета, дает возможность представить возможные работы деформаций и внешних сил в виде сумм по отдельным элементам:
а(и,у) = 1 аг (и,у), (/,у)=1 (/,у)г .
г г
Такой подход позволяет составить элементы матрицы К и вектора Р из отдельных компонентов. Так, // -элемент матрицы К и I -элемент вектора Р определяются по формулам
Kj =1 Kjr; P=2 P ,
relJ rel
где r e ij, r e i - суммирование по всем элементам, содержащим l и J узловые неизвестные; Kljr, Plr - компоненты матрицы жесткости и вектора узловых сил r - конечного элемента, которые определяются аналогично (8):
Kjr = J (вVi )T DB(<Pj )dQr; Plr = J cpjfdo.r. (10)
nr nr
Это позволяет строить разрешающую систему уравнений МКЭ (7) на основе рассмотрения каждого отдельного конечного элемента.
После выбора системы базисных функций {(} процедура МКЭ представляется достаточно формализованной. Выбор {(} является самым ответственным этапом, так как определяет сходимость метода, точность решения задачи, разрешимость системы (7). Линейная независимость базисных функций проверяется достаточно легко и, как правило, выполняется для МКЭ автоматически.
На основе работы [4] можно сформулировать требования, которым должны удовлетворять функции Щ, чтобы обеспечить сходимость МКЭ:
- система базисных функций {(}должна принадлежать энергетическому
пространству НА дифференциального оператора задачи А . То есть представление разрешающей функции „ должно удовлетворять главным граничным условиям и обеспечивать существование по всей области О тех перемещений и их производных, которые входят в функционал (4);
- функции (р I должны быть линейно независимыми, для разрешимости системы (7);
- система базисных функций {(} должна быть полна в энергетическом пространстве оператора А .
Принадлежность к энергетическому пространству оператора А устанавливается существованием компонентов напряженно-деформированного состояния, которые входят в соответствующий функционал.
По области конечных элементов это требование, как правило, удовлетворяется автоматически, поэтому необходимо проверять неразрывность соответствующих компонентов только по линиям или поверхностям контактов конечного элемента.
Так как для изгибаемых плит и оболочек дифференциальный оператор А имеет четвертый порядок, а в функционал Лагранжа входят вторые производные по перемещениям, то необходимо обеспечить непрерывность, как перемещений, так и их первых производных по линиям контакта.
Обсуждение результатов. Основываясь на выше приведенных выкладках, было принято решение об использовании МКЭ при решении задачи расчета
трехслойных балок. В качестве расчетной схемы использовалась система вертикальных и горизонтальных пластин, со следующими характеристиками:
1. Несущие слои (горизонтальные пластины) - толщина 1мм, материал: алюминий Д16-АТ (Е=69000МПа, у=0.34, р=26.8кН/м3);
2. Сотовый заполнитель (вертикальные пластины) - толщина 0.12мм, материал: алюминий АМГ2-Н (Е=69000МПа, у=0.34, р=26.8кН/м3). Расчетная длина балки принята 70см. Высота балки - 6см. Ширина балки
варьировалась в зависимости от размера сот и составляла 4 полных соты. Ширина элементов, образующих соты составляла 5, 10, 15 и 25мм.
Условие опирания - жесткие заделки. Нагрузка на балку - равномерно распределенная, интенсивностью от 1,0 - 6,5кН/м .
Исходя из вида работы и напряженно-деформированного состояния элементов в качестве расчетных были приняты четырех- и трех- узловые конечные элементы пологих оболочек. Сопряжение элементов между собой принято жестким, что обеспечивает выполнение условия совместности перемещений и их производных входящих в функционал (4). Линейная независимость базисных функций проверяется для МКЭ автоматически. Полнота системы базисных функций также выполняется, так как порядок полинома, выражаемого линейными комбинациями функции (рш, р=т (р =2, порядок полинома; 2т =4 - порядок дифференциального оператора А ). Таким образом, сходимость решений доказана.
Результаты численного эксперимента приведены на рис. 2.
Рисунок 3- График прогибов полученных в работе В.Н. Кобелева [1]. £™х=0.389мм.
Рисунок 2- График прогибов полученных в результате расчета с применением МКЭ (ПК ЛИРА).
£пах=0.39мм.
Для сравнения приведены результаты расчета ТБ по методу изложенному в [1] (рис. 3).
Вывод. Сравнение результатов аналитического и численного методов, показывает достаточно хорошее совпадение вычисляемых параметров, что подтверждает правильность выбранной модели метода конечных элементов и принятый тип конечных элементов.
Библиографический список:
1. Кобелев В.Н. «Расчёт трёхслойных конструкций». М.:Машиностроение, 1984. 304с.
2. Панин В.Ф. «Конструкции с сотовым заполнителем». М.:Машиностроение, 1982. 152с.
3. А.С. Городецкий, И.Д. Евзеров. Компьютерные модели конструкций. Москва: Издательство Ассоциации строительных вузов, 2009.360с.
4. Михлин С.Г. Вариационные методы математической физики. - М.:Наука, 190. - 512 с.
References:
1. Kobelev V. N. «Calculation of sandwich structures». M.:Mashinostroenie, 1984. 304р.
2. Panin, V. F., «design of honeycomb». M.:Mashinostroit-nie, 1982. 152р.
3. A. S. Gorodetski, I. D. Evzerov. Computer models of designs. Moscow: Publishing house Association building universities, 2009.360 р.
4. Mikhlin S. G. Variational methods in mathematical physics. - M.:Nauka, 190. -512 p.
УДК 550.34; 622, 831 Магомедов Р.А.
ОЦЕНКА НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ МАССИВА ГОРНЫХ ПОРОД
Magomedov R.A.
EVALUATION OF STRESS-STRAIN STATE OF ROCK MASS
Аннотация. Предложен новый способ определения напряженно-деформированного состояния массива горных пород без пропуска отдельных его участков на основе совмещения наблюдений за динамикой пьезометрического уровня в специальных гидрогеологических скважинах заранее спланированной сети и изменением диаметра скважин измеряемого фотоупругими датчиками.
С началом формирования напряженно-деформированного состояния геологической среды, хаотичная структура гидрогеодинамического поля (ГГД-поля) переходит в аномальное состояние определенной формы, в соответствии с геометрической формой и размерами геологической структуры. Замкнутые линии равных аномальных значений пьезометрических уровней образуют кольцевые геометрические структуры, в центральных частях которых наблюда-
132
