Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

ОПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ В ЗОНЕ СОЕДИНЕНИЯ ОБОЛОЧЕК ВРАЩЕНИЯ ИЗ РАЗНОРОДНЫХ МАТЕРИАЛОВ НА ОСНОВЕ МКЭ

Николаев А.П. 1 Киселёв А.П. 1 Гуреева Н.А. 1 Киселёва Р.З. 1 Леонтьева В.В. 1
1 ФГБОУ ВПО «Волгоградский государственный аграрный университет»
Для определения напряженно-деформированного состояния в зонах пересечения произвольно нагруженных оболочек вращения из разнородных материалов на основе метода конечных элементов используется ранее разработанный объемный шестигранный конечный элемент с узловыми неизвестными в виде перемещений и их производных. Для конечных элементов, примыкающих к границе сочленения оболочек вращения, получены соотношения между узловыми неизвестными одной оболочки, принятой за основную, и узловыми неизвестными другой оболочки, примыкающей к основной. На основе полученных соотношений выполнены преобразования матриц жесткости и векторов узловых нагрузок конечных элементов, примыкающих к границе сочленения оболочек из разнородных материалов. На основе анализа результатов расчета можно сделать вывод о корректности алгоритма определения напряженно-деформированного состояния в зонах сочленения оболочек вращения при произвольном нагружении.
МКЭ
произвольно нагруженные оболочки
объёмный шестигранный конечный элемент
узловые неизвестные
условия на границе пересечения оболочек
разнородный материал
1. Голованов А.И., Тюленева О.Н., Шигабудинов А.Ф. Метод конечных элементов в статике и динамике тонкостенных конструкций. – М.:Физматлит, 2006. – 392 с.
2. Бате К.Ю. Методы конечных элементов. – М.: Физматлит, 2010. – 1022 с.
3. Николаев А.П., Клочков Ю.В., Киселев А.П., Гуреева Н.А. Расчет оболочек на основе МКЭ в двумерной постановке. – Волгоград, 2009. – 194 с.
4. Киселев А.П. Векторная аппроксимация полей перемещений объемного шестигранного конечного элемента // Строительная механика инженерных конструкций и сооружений. – 2007. – № 1.
5. Киселев А.П. Объемный конечный элемент в виде треугольной призмы с первыми производными узловых перемещений / А.П. Киселев, А.П. Николаев // Изв. Вузов, сер. «Строительство». – 2006. – № 1. – С. 13–18.
6. Гуреева Н.А. Восьмиугольный конечный элемент в смешанной формулировке на основе функционала Рейснера // МБТУ им.Баумана, Известия вузов: Машиностроение. – № 5. – С. 23–28.
7. Гуреева Н.А. Расчет многослойной оболочки с использованием объемного конечного элемента / Н.А. Гуреева, А.П. Киселев, Р.З. Киселева // Известия ВолгГТУ. – Волгоград, 2010. – № 4. – С. 125–128.
8. Киселев А.П. Расчет многослойных оболочек вращения и пластин с использованием объёмного конечного элемента / А.П. Киселев, Н.А. Гуреева, Р.З. Киселёва // Изв. Вузов, сер. «Строительство». – 2010. – № 1. – С. 106–112.
9. Киселёв А.П. Использование трёхмерных конечных элементов в расчётах прочности многослойных панелей / А.П. Киселев, Н.А. Гуреева, Р.З. Киселёва // Строительная механика инженерных конструкций и сооружений. – 2009. – № 4. – С. 37–40.
10. Седов А.И. Механика сплошной среды. – М:. Наука, 1976. – т. 1,535.

Из-за сложности решения дифференциальных уравнений, описывающих деформированное состояние оболочек вращения, большое распространение получили численные методы определения их напряженно-деформированного состояния. Среди численных методов особое место занимает метод конечных элементов (МКЭ) в различных формулировках: в формулировке метода перемещений разрабатывались конечные элементы в двумерной постановке [1, 2, 3] и в трехмерной постановке [4, 5]; в смешанной формулировке использовались объемные конечные элементы [6]. Объемные конечные элементы в формулировке метода перемещений успешно использовались для расчета слоистых конструкций [7, 8, 9].

В настоящей работе для расчета произвольно нагруженной оболочки вращения в координатной системе s, θ, ζ используется шестигранный восьмиузловой конечный элемент с узлами i, j, k, l на нижней грани по координате ζ и узлами m, n, p, h по верхней грани [4].

Используемая в настоящей работе матрица жесткости объёмного шестигранного конечного элемента формируется на основе равенства работ внешних и внутренних сил [4, 5] и представляется выражением

nikol01.wmf (1)

где nikol02.wmf – вектор узловых неизвестных в криволинейной системе координат s, θ, ζ;

nikol03.wmf

(ϖ = i, j, k, l, m, n, p, h) [K] – матрица жесткости элемента в глобальной системе координат; {f} – вектор узловых нагрузок элемента в глобальной системе координат.

1. Геометрия оболочки вращения в узловой точке. Положение произвольной точки М срединной поверхности произвольно нагруженной оболочки вращения в декартовой системе координат xoz определяется радиус-вектором (рис. 1)

nikol04.wmf (2)

где r = r(x) – радиус вращения точки М относительно оси ox; nikol05.wmf – орты декартовой системы координат; θ – угол, отсчитываемый от вертикального диаметра против часовой стрелки.

pic_10.tif

Рис. 1. Перемещение точки в результате деформирования оболочки из положения Mζ в положение Mζ*

Векторы локального базиса точки М определяются выражениями

nikol06.wmf

nikol07.wmf (3)

где nikol08.wmf – производная радиуса вращения по дуге меридиана s.

Соотношения (3) можно представить в матричном виде

nikol09.wmf

где nikol10.wmf (4)

Производные векторов локального базиса определяются дифференцированием (3) и с учетом (4) представляются в матричном виде

nikol11.wmf (5)

где nikol12.wmf

Радиус-вектор произвольной точки оболочки Mζ, отстоящей на расстоянии ζ от срединной поверхности, можно представить выражением

nikol13.wmf (6)

Базисные векторы точки Mζ определяются дифференцированием (5)

nikol14.wmf

nikol15.wmf (7)

nikol16.wmf

Произвольная точка Mζ оболочки под действием заданной нагрузки займет положение Mζ*, которое определяется вектором nikol17.wmf с компонентами в базисе точки M срединной поверхности

nikol18.wmf (8)

Производные вектора перемещения по координатам s, r, θ, ζ с учётом (5) имеют вид

nikol19.wmf

nikol20.wmf

nikol21.wmf (9)

Деформации в точке Mζ* определяются выражениями

nikol22.wmf (10)

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

nikol23.wmf (11)

где nikol24.wmf – вектор – строка компонент деформаций в узловой точке оболочки.

Связь между напряжениями и деформациями определяется соотношениями механики сплошной среды [10]

nikol25.wmf (12)

где λ, μ – параметры Ламе;

nikol26.wmf

– первый инвариант тензора напряжений; gmn, gmn – ковариантные и контравариантные компоненты метрического тензора.

Зависимость (12) для узловой точки можно представить в виде

nikol27.wmf (13)

где nikol28.wmf.

2. Преобразование узловых величин в точке на грани сочленения оболочек из разнородных материалов. Рассматриваются две произвольно нагруженные оболочки вращения в декартовых системах координат xyz и nikol29.wmf. Связь между ортами этих систем считается известной (рис. 2)

nikol30.wmf (14)

где nikol31.wmf

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

nikol33.wmf (15)

nikol34.wmf (16)

pic_11.tif pic_12.tif

Рис. 2. Оболочки вращения в декартовых системах координат xoz и nikol32.wmf

Для этого в узловой точке на заданной поверхности пересечения оболочек используется ортогональный базис nikol35.wmf, определенный через векторы базиса декартовых координат

nikol36.wmf (17)

где nikol37.wmf

С использованием (4), (7) можно по (17) получить матричное соотношение

nikol38.wmf (18)

где базисы

nikol39.wmf nikol40.wmf

nikol41.wmf

относятся соответственно к основной и примыкающей оболочкам.

Векторы nikol42.wmf лежат в плоскости грани пересечения, а вектор nikol43.wmf нормален к поверхности пересечения оболочек.

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

nikol44.wmf (19)

nikol45.wmf (20)

Между векторами (19) и (20) записывается матричная зависимость

nikol46.wmf (21)

где [I] – матрица, на главной диагонали которой элементы равны единице.

Для определения соотношений между компонентами векторов (20), (16) используются следующие условия.

1. Условие о равенстве векторов перемещений в базисах nikol47.wmf и nikol48.wmf

nikol49.wmf

откуда с использованием (18) определяются компоненты nikol50.wmf через nikol51.wmf

nikol52.wmf (22)

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

nikol53.wmf (23)

2. Используется выражение производной скаляра a по направлению nikol54.wmf в криволинейной ортогональной системе координат

nikol55.wmf(24)

На основании (24) можно записать выражения

nikol56.wmf

nikol57.wmf

nikol58.wmf (25)

Для примыкающей оболочки можно записать аналогичные соотношения

nikol59.wmf

nikol60.wmf

nikol61.wmf (26)

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

nikol62.wmf (27)

где nikol63.wmf

Используя соотношения (18), из (27) можно сформировать матричное выражение

nikol64.wmf (28)

где nikol65.wmf

Принимая во внимание (13) и (15), можно выразить напряжения в базисе nikol66.wmf через перемещения базиса nikol67.wmf

nikol68.wmf (29)

Используя базис nikol69.wmf примыкающей оболочки, можно получить матричное соотношение

nikol70.wmf (30)

На основании выражений (29), (30), (25), (26) формируются матричные соотношения

nikol71.wmf (31)

Из условия равенства nikol72.wmf и nikol73.wmf получается

nikol74.wmf (32)

С использованием (32) формируется матрица преобразования [Т] для матрицы жесткости и вектора узловых нагрузок граничного конечного элемента примыкающей оболочки

nikol75.wmf (33)

Пример № 1. Определялось напряженнo-деформированное состояние цилиндра со сферическим днищем, находящегося под действием внутреннего давления интенсивности q (рис. 3). Цилиндр и днище выполнены из разнородных материалов.

pic_13.tif

Рис. 3. Цилиндр со сферическим днищем под действием внутреннего давления интенсивности q

Были приняты следующие исходные данные: l1 = 0,2 м, l2 = 0,1 м, l3 = 0,09 м, q = 10 Н, h = 0,0005 м, b = 0,05 м, E = 2·105 МПа, ν = 0,3, E′ = 2·106 МПа, ν′ = 0,25.

Конструкция разбивалась на 10 конечных элементов по толщине, на 100 элементов по длине цилиндра и на 50 по дуге круговой оболочки.

По полученным результатам построена эпюра нормальных напряжений σxx (рис. 4) в сечении 1–1 (рис. 3). Условие равновесия по силам (ΣX = 0) правой части оболочки от сечения 1-1 выполняется с погрешностью δ = 0,6 %.

На основе анализа результатов выполненного примера расчета можно сделать вывод о корректности алгоритма определения напряженно-деформированного состояния в зонах сочленения оболочек вращения на основе разработанного конечного элемента [4].

pic_14.wmf

Рис. 4. Эпюра нормальных напряжений σxx в сечении 1–1 цилиндра со сферическим днищем

Рецензенты:

Голованов В.К., д.т.н., профессор кафедры «Начертательная геометрия и графика» ВГТУ, г. Волгоград;

Кукса Л.В., д.т.н., профессор, заведующий кафедрой «Сопротивление материалов», ВолГАСУ, г. Волгоград.

Работа поступила в редакцию 26.02.2014.


Библиографическая ссылка

Николаев А.П., Киселёв А.П., Гуреева Н.А., Киселёва Р.З., Леонтьева В.В. ОПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ В ЗОНЕ СОЕДИНЕНИЯ ОБОЛОЧЕК ВРАЩЕНИЯ ИЗ РАЗНОРОДНЫХ МАТЕРИАЛОВ НА ОСНОВЕ МКЭ // Фундаментальные исследования. – 2014. – № 5-2. – С. 269-275;
URL: http://www.fundamental-research.ru/ru/article/view?id=33863 (дата обращения: 13.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074