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

Серия «Естественные науки» Устойчивость трёхслойной цилиндрической оболочки из композитного материала с ортотропным промежуточным слоем

к.т.н. доц. Мозгунов В.Н., к.т.н. доц. Серов М.В., к.т.н. проф. Александрова С.Г.,

к.т.н. доц. Зубарев А. А.

МГТУ им. Н.Э. Баумана, Университет машиностроения, НИТУ «МИСиС»

8(495)223-05-23, tm@mami.ru

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

Ключевые слова: устойчивость, трёхслойная цилиндрическая оболочка, ор-тотропный промежуточный слой.

Введение

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

Наибольшая часть публикаций, посвящённых проблеме статической устойчивости рассматриваемых оболочек, связана с упрощёнными моделями для промежуточного слоя [1]. Некоторые модели для заполнителя, как упругого основания представлены в работе [2]. В работе [3] изложены инженерные методы расчёта трёхслойных цилиндрических оболочек из композитных материалов, подкреплённых рёбрами жесткости и упругим цилиндром. Особый интерес представляют точные решения, полученные без каких-либо допущений о напряжённо-деформированном состоянии в промежуточном слое, так как эти решения позволяют проконтролировать точность вводимых допущений и уточнить границы применения технических теорий. Если промежуточный слой - ортотропный цилиндр, то решение можно построить, представляя перемещения в виде обобщённых степенных рядов. Такая теория развита в работе [4] и используется в данной статье.

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

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

Постановка задачи

Рассматривается устойчивость трёхслойной цилиндрической оболочки под действием внешнего гидростатического давления (рисунок 1 а). Материал наружного и внутреннего тонкостенных слоёв считается ортотропным.

Промежуточный слой представляет собой ортотропный цилиндр. Для обшивки исполь-

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

а

г I ' ' '

б)

h

Рисунок1. Расчётная схема трёхслойной цилиндрической оболочки с ортотропным промежуточным слоем: а) оболочка; б) фрагмент трёхслойного пакета

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

Исследования и результаты

1. Уравнения устойчивости трёхслойного пакета

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

4U+LVr + 43 W ± ч1а = о; L2U + L2V + LW ± чь=о;

LUi + L2V +

r('k

M) + T 0

L33 ^ Tia

H

да2

+ T 0

д2 V /Г дч-с

W + h гс

dß2

да dß

+

± Чгг = 0,

где: i = 1,2 - номера слоёв обшивки соответственно наружного и внутреннего;

Ui ,Vj ,Wi - проекции вектора полного перемещения точки срединной поверхности обшивки на оси а, ß и r; Ltj - дифференциальные операторы, содержащие коэффициенты жёсткости срединной поверхности слоёв обшивки на растяжение, сдвиг, изгиб и кручение [5].

Для дальнейшего решения уравнений устойчивости (1) необходимо знать закон распределения контактных поверхностных усилий qia , qiß , qr , а также величины погонных до-

критических усилий Т°а ,Tß безмоментного состояния трёхслойной оболочки.

2. Разрешающие уравнения для промежуточного слоя

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

qia = pia cos Ха cosn j;

4iß= Piß sinХаsinn j;

qr = pirsin Ха cosn j,

)

a

h

- pm

где: l = -j- ; m,n - параметры волнообразования в осевом и окружном направлениях; l - длина ортотропного цилиндра;

Pia , Рь, Pr - амплитудные значения контактной поверхностной нагрузки. Для ортотропного цилиндра уравнения в перемещениях имеют вид [4] :

LjjU + L12V + L13W = 0;

L21U + L22V + L2W = 0;

L31U + L32V + L33W = 0;

где: Lj - дифференциальные операторы; U,V,W - проекции вектора полного перемещения некоторой точки заполнителя на оси a, b и r . Представляя перемещения точек заполнителя при потере устойчивости в виде :

U = U ( r ) cos la cos n j;

V = V

W = W (r ) sin la cos n j, и вводя новую безразмерную координату:

Р = -

r - R

R

h h — < p < —

2r 2r

Rj + R2

где: R = —-2

2

среднии радиус ортотропного цилиндра,

h -

толщина заполнителя полу-

чим систему уравнений:

a-

d2

d

(1 + P)  W + (1 + P) dP

- ^22 - G 12 n2 - G13 x2 (1 + P )2 ïw +

+n

d

12 «22

x

( «12 + G12 )(1 + P ) d[p ~ G12 -d

V -U = 0;

n

( «21 + G12 )(1 + P ) —- + G12 + a22

w-\a.

(1+p)2 w+(1+P) dP -1

2

-«22 П - G23

x2 (1 + p )2 }v-( «23 + G23 ) nx (1 + p )U = 0;

x

2 d2

w + x ( «32 + G23 ) n (1 + P )V +

G

d2

d

(1+P)  &+(1+рУdp

-G3;n2 -x2 (1 + P)4u = 0.

Здесь: - компоненты симметричного тензора модулей упругости 1-го и 2-го ро-

да; р = г ф; с = 1

Введение переменной р позволяет выбрать решение системы дифференциальных

уравнений (5) при отсутствии особых точек в области || < в виде степенных рядов:

2R

w (p)=Z A rk;

k=0 ¥

V

k=0 ¥

U (p)=Z Ck pk.

Подставляя зависимости (6) в уравнения (5) и приравнивая члены при одинаковых степенях р , получим рекуррентные соотношения для определения коэффициентов Лк,Бк,Ск :

Ак = -[к (к -1)] (апЛк- +а12Лк_2 -авЛк_3 -амЛк_4 +а15Вк- +

+а16Вк-2 - а17Ск-1 а18Ск-2 - а19Ск-3 );

Вк =[к (к 1)] [а21 Лк-1 +а22Лк-2 а23Вк-1 а24Вк-2 + а25Вк-3 а26Вк-4 а27 (Ск-2 + Ск-3 )] ; (7) Ск = - [к (к - !)] 1 [(а31 Лк-1 + а32Лк-2 + а33Лк-3 + а34 (Вк-2 + Вк-3 ) + + а35^к-1 + а36Ск-2 - а37Ск-3 - а38Ск-4 )] ,

где: а - коэффициенты, состоящие из модулей упругости, параметров волнообразования и индекса суммирования. В соотношениях (7) коэффициенты Л0,Л1,В0,В1,С0,С1 остаются неизвестными и играют роль произвольных постоянных. Задаваясь для первого случая:

Л«= 1, Л1(1)= В01)= В1(1) = С01)= С1(1)= 0;

для шестого случая:

cj6) = i, a06) = A(6) = b06) = B|6) = C06) = o,

с помощью рекуррентных соотношений (7) находим шесть частных решений

W о) ,V(j) ,pj ] (j = 1,2,-6 ).

Тогда общее решение уравнений (5) запишется в виде:

6

W = ^ DjW(j) (р) sin la cos n j;

j=i 6

V = ^ DjV

j=i 6

U = ^ Dp(j)(р) cos la cos n j,

j=i

где: О (7 = 1,2,...6) - произвольные постоянные.

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

Затем определяются шесть частных решений о[7,о^'^р ,т(а'так, как это делалось выше для определения Ж7 ,У7 ,и 7.

Тогда общее решение для определения напряжений в ортотропном цилиндре запишется в виде:

1 6

Ттт:-г Z Di °(j) (Р) sin Xa cos nj

R (1 + Р) i=i

^я =

Z DjGp) (p) sin Xa cos nj;

R(1 + p) ] b

~ы\—ч ZDja¿ (p)sin Xa cos j; R (1+p) i=i

тф =

G

12

Z DjTj ( p ) sin Xa sin n j;

R

G

Tra = ^ Z DlTÍra) (p) COs Xa COs П j;

R U ] G

Tafi

Z DjT^ (p) cos Xa sin n j.

R

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

Рассмотрим алгоритм определения внешнего критического давления для трёхслойной цилиндрической оболочки.

3. Определение верхней критической нагрузки

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

Ui = U0i cos la cos n j;

Vi = V0i sin la sin n j;

Wi = W0i sin la cos n j, где: U0i,V0i,W0i - амплитудные значения перемещений несущих слоёв оболочки.

В силу тонкостенности несущих слоёв оболочки граничным условиям по перемещениям будем удовлетворять на срединных поверхностях:

U = U (pl);

V = V

W = w (p),

где: U(pi), V(pi), W(pi) - перемещения точек наружной и внутренней поверхностей заполнителя.

Подставляя в граничные условия (11) зависимости (10) и (8) при р = р получим значения произвольных постоянных Dj :

А U )..... .....U (б)(л ) -1 U 01

D. = V )..... .....V %,) V01

D6 W е > (p.)...... W ( 6 ) (.p 2) Wo.

Для определения контактных усилий входящих в уравнение (1) будем исходить из соотношений:

°г =

1

a =

Чга =- Тга (Рг )-' % =- Т ф (Рг )-'

Чгг =- ° г

где: тт(р), тгв (р), ог (р) - нормальные и касательные напряжения на наружной и внутренней поверхностях заполнителя. Подставляя в уравнение (13) зависимости (2) для контактной нагрузки, а также зависимости (9) для напряжений заполнителя при р = , получим выражения для амплитудных значений контактной нагрузки:

Р1а и01

Р1 в у01

II Р у.

Р2г К02

где: (г, у = 1,2...6) - элементы матрицы 6-го порядка, получаемой от перемножения

матрицы, составленной из частных решений для напряжений, на обратную матрицу, составленную из частных решений для перемещений. Подставляя в уравнение (1) зависимости (10) и (2) с учётом амплитудных значений контактных усилий (14) получаем уравнение устойчивости трёхслойной оболочки в матричном виде:

ф[1> - q

Фк1 рк1

— Р+3,1

- Q

к,1+3

ф(2) - Q

Фк1 рк+3,1+3

и

01

К

02

=Ф • и =0

где: Ф() (к,1 = 1,2,3) - элементы блочной матрицы, определяющиеся следующим образом:

Ф

-А I2 -Р —' Ф(0 = Ф(0 = (а + Р )—' Ф(0

1 га ^ 2 > ^12 ^21 V гаР га / ^ > ^13

Ф

А

"га р

Я

1;

Ф

(г) 22

п

-АР — - Р 12; Ф23

гр я2 га ' 23

Ф

32

\ га р га /

п 1

Ф () =

33

в. 14 + 2(В. р + 2РР)12— + В.р— + -

га V гар гр / т-,2 г р п4 п

/

Я 2 г р Я4 ' Я2

< V

п

А р - В. р12 + Вр —

гр гар гр ^2

Л

' У

Т012 - т0 п

га гр я 2 .

г

Нетривиальное решение уравнения (15) возможно, если определитель

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

710 ГТ10 "

1аТгр, позволяет получить квадратное уравнение относительно внешней нагрузки д :

а

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

Заключение

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

Предложенный алгоритм расчёта был апробирован на известных решениях задач устойчивости трёхслойных цилиндрических оболочек с лёгким заполнителем из пенопласта и с радиально-армированным заполнителем из стеклопластика, для которых ранее использовался модельный подход. Так верхнее критическое давление для оболочек с радиально-армированным заполнителем цгкр = 17,9 МПа достаточно хорошо согласуется с результатами

эксперимента = 14,6 МПа [6].

Литература

1. Сухинин С.Н. Прикладные задачи устойчивости многослойных композиционных оболочек. М.: Физматлит, 2010, 248 с.

2. Соломонов Ю.С. Георгиевский и др. Методы расчёта цилиндрических оболочек из композиционных материалов. М.: Физматлит, 2009, 262 с.

3. Иванов О.Н., Александрова С.Г., Мозгунов В.Н. Некоторые модели упругого анизотропного основания. - Машины и технологии переработки каучуков, полимеров и резиновых смесей. Ярославль, 1982, с. 106-111.

4. Гусь А.И. Устойчивость трёхмерных деформируемых тел. Киев: Наукова думка, 1971, 276 с.

5. Иванов О.Н., Мозгунов В.Н., Рогинский С.Л., Устойчивость трёхслойной радиально-армированной цилиндрической стеклопластиковой оболочки. В кН.: Совершенствование конструкций машин и аппаратов химических производств. М.: МИХМ, 1982. с. 23-27.

6. Иванов О.Н., Рогинский С.Л., Плешков Л.В., Мозгунов В.Н. Экспериментальное исследование устойчивости радиально-армированных цилиндрических оболочек из стеклопластика. В кН.: Стеклопластики и стекловолокно. М., 1982, вып.5, с. 6-9.