Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

DETERMINATION OF TENSION IN ZONE SOYEDENENY OF COVERS OF ROTATION FROM DIVERSE MATERIALS ON THE BASIS OF A METHOD OF FINAL ELEMENTS

Nikolaev A.P. 1 Kiselev A.P. 1 Gureeva N.А. 1 Kiseleva R.Z. 1 Leonteva V.V. 1
1 Volgograd State Agricultural University
1499 KB
For definition strained the deformed condition in crossing zones of arbitrary loaded c shells of rotation from diverse materials on the basis of a method of finite elements earlier developed volume six-sided finite element with nodal unknown in the form of displacements and their derivatives is used. For the finite elements adjacent to border of a joint of shells of rotation, relation between nodal unknown of one cover taken for main, and nodal unknown of other shells adjoining main are received. On the basis of the received relations transformations of matrixes of rigidity and vectors of nodal loadings of the finite elements adjoining border of a joint of shells from diverse materials are executed. Based on the analysis of the calculation results can be concluded about the correctness of the algorithm for determining the stress-strain state in the areas sochlineniya shells of revolution under arbitrary loading.
finite element method
аrbitrarily loaded shells of revolution
surround hexagonal finite element
the nodal unknowns
the conditions at the border crossing membranes
heterogeneous material
1. Golovanov A.I., Tuleneva O.N., Shigabudinov A.F. Finite Element Method in statics and dynamics of thin-walled structures. M. Fizmatlit, 2006. 392 p.
2. Bath K.Y. Finite element methods. M. Fizmatlit, 2010. 1022 р.
3. Nikolaev A.P., KlochkovYu.V., Kiselev A.P., Gureeva N.A. Calculation of shells on the basis of FEM in two-dimensional formulation. Volgograd, 2009. 194 p.
4. Kiselev A.P. Stock approximation displacement fields surround hexagonal finite element Sci-Tech. journal «Structural Mechanics engineering constructions and buildings» no. 1, People’s Friendship University, Moscow, 2007.
5. Kiselev A.P. Volumetric finite element in the form of a triangular prism with the first derivatives nodal displacements / A.P. Kiselev, A.P. Mykolaiv / Math. Universities, Ser. «Construction». 2006. no. 1. рр. 13–18.
6. Gureeva N.A. Octagonal mixed finite element formulation based on the Reissner functional. MBTU Bauman Education News : Machinery, M.: no. 5, рp. 23–28.
7. Gureeva N.A. Calculation of multilayer cladding surround finite element / N.A. Gureeva, A.P. Kiselev, R.Z. Kiseleva News VSTU. Volgograd, 2010. no. 4. рр. 125–128 .
8. Kiselev A.P. Calculation of multilayer shells of revolution and plates surround the finite element / A.P. Kiselev, N.A. Gureeva, R.Z. Kiseleva Math. Universities, Ser. «Construction». 2010. no. 1. рр. 106–112.
9. Kiselev A.P. Using three-dimensional finite element calculations of strength of sandwich panels / A.P. Kiselev, N.A. Gureeva R.Z. Kiseleva Structural Mechanics engineering constructions and buildings. 2009. no. 4. рр. 37–40.
10. Sedov A.I. Continuum Mechanics. M.: «Science», 1976, v. 1, 535 p.

Из-за сложности решения дифференциальных уравнений, описывающих деформированное состояние оболочек вращения, большое распространение получили численные методы определения их напряженно-деформированного состояния. Среди численных методов особое место занимает метод конечных элементов (МКЭ) в различных формулировках: в формулировке метода перемещений разрабатывались конечные элементы в двумерной постановке [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.