Совершенствование схемы аппроксимации перемещений в стержневом конечном элементе для расчета подкрепленных пластин Текст научной статьи по специальности «Механика и машиностроение»

Численные методы расчета конструкций

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

И.Н. СЕРПИК*, д-р техн. наук, профессор

М.В. ШВЫРЯЕВ**, инженер

А.И. ТЮТЮННИКОВ***, инженер

* Брянская государственная инженерно-технологическая академия;

** Брянский государственный технический университет;

*** ООО «Научно-инженерная компания», Жуковский, Московская обл.

241007, РФ, Брянск, ул. В. Сафроновой, д. 75, 29, [email protected]

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

КЛЮЧЕВЫЕ СЛОВА: пластинчато-стержневые системы, метод конечных элементов, стержневой конечный элемент, условие совместности перемещений.

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

Рассмотрим стержень, имеющий постоянное по длине поперечное сечение с главными центральными осями Cz и Су (рис. 1). Ось Сz параллельна срединной плоскости П подкрепляемой обшивки. Полагаем, что стержень может быть подвергнут растяжению-сжатию, поперечному изгибу в двух главных плоскостях и чистому кручению. Центр изгиба совпадает с центром тяжести сечения. Считаем для стержня справедливыми классические гипотезы плоских сечений и ненадавливания волокон.

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

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

и° = иаЬ - еу07, (1)

где иаь - продольное перемещение точки прямой аЬ для рассматриваемого поперечного сечения стержня; ву - проекция вектора в на ось Су; 0г - угол поворота данного сечения относительно оси О2.

Рис. 1. Соединение стержня 1 с тонким листом 2: в - вектор эксцентриситета присоединения стержня к обшивке

Рис. 2. Плоский стержневой конечный элемент

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

к }={

е х

X г

ди'

дх д ^

дх

(2)

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

Соотношение (2) с учетом равенства (1) запишем таким образом:

{е в }=•

диаЬ

д 2 V

дх

- в

У дх2

д 2 V

дх

2

(3)

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

{5 в }42, V(1), 0^ иаь>, V(2), 0^2) }т,

(4)

где и«, v(г),0(/) - перемещения иаь, V и угол поворота 0г для узла / (/ = 1, 2).

Будем аппроксимировать перемещение иаь по линейному закону, а V - с помощью полинома третьей степени [1]. Тогда матрица деформаций [2] конеч-

ного элемента

к 1=

- V1 - вуФ1 - ву Ф2 V1 - ву Фз - ву Ф4 0 -ф1 -ф2 0 -ф3 -ф4

(5)

где

6 12 х

Ф1 =—2 + —Т" /2 /3

4 6х 6 12х

Ф2 = -7+тт; Фз = /2 -7т

2 6х

7+7

Ф4 = -- + -2"

Вектор обобщенных напряжений стержня

к }=

N М7

(6)

где N - продольная сила; М г - изгибающий момент относительно оси С7. Матрица упругости [2] конечного элемента

~ЕЛ 0

О ] =

0 Ей

(7)

где Е - модуль упругости материала; А - площадь поперечного сечения стержня; 1г - момент инерции поперечного сечения относительно оси С7.

Используя равенства (5)-(7) и зависимость для вычисления матрицы жесткости конечного элемента, занимающего объем V, [2]

К ]=|[Ве Г [О. 1Ве ,

(8)

V

получим

Ке ] =

/

I

0

(еУ / + п

- е /

у•/

12

Т\" у

)

-/ 0

0

еу/

I (еу / + п) 4^ / + п) еу/

у /

-(е;/+п) -

1 (еу-0

2 V у-1

6 (еу / + п)

2^ / + п)

- еу/

(еу / + п)

(еу / + п ) 0

Симметрично 12 (еу / + п)

Здесь / = ЕА/1; п = Е1г/1; I - длина конечного элемента. В соответствии с формулами (5)-(7) внутренние усилия N(0, М« ниях по узлам 1, 2 можно определить таким образом:

6

4(еу / + п)

(9)

в сече-

N(1) М« N(2)

где матрица

[5. ] =

О ] 00 [о. ]

М

1 6е^ I 12 6 I2 6е,

(2)

= [5е №е },

(10)

4е,

6е 2е

у у

I I

I

0 —

- 0 -— -

_ 1

-1 0 - -6-

4

1 2е

6

12

1 2

1

12 в_

7

2

1 4е,.

1 4

1

(11)

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

1

2

I

{5в }о = '"""

(12)

Г&Н

^ =\{5в} /,

где {5в}1= {м(1),в<1),в^,м(2),в^,0У2)}Т ; (13)

), в^х), ву) - значения проекции м перемещения узла г на ось Сг и углов поворота вх, в у поперечного сечения в этом узле относительно осей Сх и Су

(г =1, 2).

Тогда матрицу жесткости конечного элемента, отражающую все рассматриваемые деформации, можно записать следующим образом:

К 1о =

К1 о о К 11

(14)

где [Кв ]1 - матрица, связанная с деформациями изгиба в плоскости Охг и чистого кручения.

Матрица [Кв ] формируется по традиционной схеме [1] с учетом эксцентричного присоединения стержня к узлам на жестких консолях.

Пример

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

Рис. 3. Подкрепленная консольная пластина: 1 - стержень; 2 - обшивка

Задавалось: толщина листа к = 3 мм; сечение стержня - квадратная труба 60х60х3 мм; материал листа и стержня - сталь с модулем упругости £=2,06-105 МПа и коэффициентом Пуассона V = 0,3. На пластину действует сосредоточенная сила F = 1 кН.

N кН

-8 -10 -12 -14 -16 -18 -20

Н1

М2, Нм

750 700 650 600 550 500 450 400 350

Н1

• Предлагаемый конечный элемент —■— Nastran NX 7.5

--------------. ■-------------

Н2 Сетка

а)

Н3

Предлагаемый конечный элемент

Nastran NX 7.5

Работа обшивки моделировалась с помощью треугольного трехузлового конечного элемента, в котором мембранные перемещения аппроксимировались по линейному закону, а прогиб - на основе методики, представленной в работе [3]. Рассматривались равномерные сетки Hj (] = 1, 2, 3) конечных элементов разной густоты. На рис. 3 сплошными толстыми линиями показана сетка Н\, сплошными толстыми и тонкими линиями - сетка Н2, сплошными и штриховыми линиями -сетка Н3. Полученные решения сравнивались с результатами, установленными для тех же сеток с помощью пакета конечно-элементного анализа

NASTRAN Ж 7.5 (имеется лицензия на использование у Брянского государственного технического университета). В этом пакете треугольные трех-узловые конечные элементы пластин близки по точности к реализуемому нами подходу к описанию деформаций тонкой обшивки, а для моделирования стержней применяются классические конечные элементы на жестких консолях.

На рис. 4, 5 проиллюстрирована скорость сходимости рассматриваемых решений по мере сгущения сеток. Штриховыми линиями здесь отражены результаты, полученные с помощью пакета NASTRAN NX 7.5 для уточненной мелкой сетки Нь, дальнейшее сгущение которой не приводит к сколько-нибудь значительным изменениям результатов счета.

На рис. 6 приведены эпюры продольных сил и изгибающих моментов в стержне, найденные для сетки Н2. Кроме того, на этих схемах штриховыми линиями показаны графики изменения вдоль стержня внутренних усилий, полученные в программном комплексе NASTRAN NX 7.5 с помощью сетки Нь.

Н3

Н2 Сетка

б)

Рис. 4. Внутренние усилия в стержне по сечению заделки: а - продольная сила; б - изгибающий момент

vK

-1,8 -1,9 -2 -2,1 -2,2 -2,3 -2,4 -2,5

Предлагаемый конечный элемент

H1

H2

Сетка

H3

Рис. 5. Вертикальное перемещение стержня в узле К

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

13796

562

4583

354

К Я

lvlz> пм 144

а) б)

Рис. 6. Эпюры внутренних силовых факторов в стержне для сетки H2: а - расчет в пакете Nastran NX 7.5; б - расчет с помощью разработанного стержневого

конечного элемента Л и т е р а т у р а

1. Агапов, В.П. Метод конечных элементов в статике, динамике и устойчивости пространственных тонкостенных подкрепленных конструкций / В.П. Агапов. - М.: Изд. АСВ, 2000. - 152 с.

2. Zienkiewicz, O.C. The finite element method, Vol. 2: Solid mechanics/ O.C. Zienkiewicz, R.L. Taylor. - Oxford: Butterworth-Heinemann, 2000. - 459 pp.

3. Serpik, I.N. Development of a new finite element for plate and shell analysis by application of generalized approach to patch test / I.N. Serpik // Finite Elements in Analysis & Design. - 2010. - Vol. 46, № 11. - P. 1017-1030.

IMPROVEMENT OF AN APPROXIMATION SCHEME FOR ROD FINITE ELEMENT DISPLACEMENTS TO ANALYZE THE REINFORCED PLATES

Serpik I.N., Shviryaev M.V., Tyutyunnikov A.I.

The rod finite element providing displacement compatibility in case of eccentric connection a rod with a thin plate in the direction along the interface circuit is developed. This approach to the approximation of displacements allows to reduce the discontinuities of the rod internal forces on nodes. The results of the reinforced plate calculation showing a significant increasing of the solution accuracy in the comparison with the traditional procedure of discretization of plate-rod systems are presented.

KEY WORDS: plate-rod systems, finite element method, rod finite element, displacement compatibility condition.