CC BY
16
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
2017, Т. 159, кн. 4 С. 395-428
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 539.3
ОСЕСИММЕТРИЧНЫЕ ЗАДАЧИ О ГЕОМЕТРИЧЕСКИ НЕЛИНЕЙНОМ ДЕФОРМИРОВАНИИ И УСТОЙЧИВОСТИ ТРЕХСЛОЙНОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ С КОНТУРНЫМИ ПОДКРЕПЛЯЮЩИМИ СТЕРЖНЯМИ
И.Б. Бадриев1, М.В. Макаров1, В.Н. Паймушин1,2,
С.А. Холмогоров2
1 Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия 2Казанский национальный исследовательский технический университет имени А.Н. Туполева, г. Казань, 420111, Россия
Аннотация
Проведено численное исследование задачи о геометрически нелинейном осесимметричном деформировании трехслойной цилиндрической оболочки с трансверсальномягким заполнителем, подкрепленной в торцевых сечениях упругими стержнями. Для описания процесса деформирования использованы выведенные ранее уравнения уточненной геометрически нелинейной теории, позволяющие как изучить докритическое поведение оболочки, так и выявить все возможные формы потери устойчивости несущих слоев. Указанные уравнения основаны на введении в рассмотрение в качестве неизвестных контактных усилий взаимодействия внешних слоев с заполнителем, а также внешних слоев и заполнителя с подкрепляющими телами во всех точках поверхностей их сопряжения. Разработаны численные методы решения сформулированных задач. Они основаны на предварительном сведении исходных задач к системе интегро-алгебраических уравнений, при решении которых используется метод конечных сумм. Предложена методика исследования докритического и закритического геометрически нелинейного поведения оболочки при ее торцевом сжатии через контурные подкрепляющие стержни, согласно которой неустойчивые положения равновесия определяются методом продолжения решения по параметру при выборе в качестве параметра работы внешних сил. Предложен способ нахождения критической нагрузки (точки бифуркации), при достижении которой оболочка теряет устойчивость. Он основан на линеаризации исходной геометрически нелинейной задачи в окрестности её нелинейного решения с последующей формулировкой задачи на собственные значения с нелинейным вхождением параметра. Приведены результаты численных экспериментов. Проведен анализ результатов экспериментов.
Ключевые слова: трехслойная цилиндрическая оболочка, трансверсально-мягкий заполнитель, контурный подкрепляющий стержень, геометрическая нелинейность, контактные напряжения, осевое сжатие, осесимметричная деформация, метод конечных сумм, докритическое и закритическое поведение, точка бифуркации, линеаризованная задача, формы потери устойчивости
Введение
Трехслойные пластины и оболочки являются типовыми элементами конструкций того или иного назначения [1—5], в частности в судостроении [6-9], авиационной и ракетно-космической промышленностях [10-14]. Задачи механики их деформирования, как правило, формулируются без учёта их взаимодействия с контурными
Рис. 1. а) Схема соединения трехслойной оболочки со шпангоутом; б) схема подкрепляющего стержня
подкрепляющими стержнями. Результаты их исследований можно найти в работах [15-20], обширная библиография, посвященная таким задачам, содержится, в частности в обзорной статье [21].
Отдельное направление исследований в области механики трехслойных конструкций составляют задачи устойчивости, анализ которых дан в статье [22]. Отметим, что численные методы исследования процесса геометрически и физически нелинейного деформирования, включая вопросы выявления форм потери устойчивости (ФПУ) трехслойных пластин, освещены в статьях [23-28].
Трехслойные оболочки, составленные из двух жестких несущих слоев 1, 2 и расположенного между ними заполнителя 3, относящийся к классу трансверсальномягких, как правило, в реальных конструкциях соединяются с другими конструктивными элементами через шпангоуты 4 (рис. 1, а). Такое соединение зачастую выполняется так, что осевая сжимающая сила, приложенная к шпангоуту, в силу особенностей конструктивного исполнения трехслойной оболочки на контуре (рис. 1, б) главным образом передаётся на наружный несущий слой. Так как такие оболочки являются наиболее рациональными при их работе на изгиб и сжатие, одной из главных при их прочностном анализе является задача исследования возможных ФПУ.
В работах [29, 30] был проведен анализ конструктивных вариантов исполнения трехслойных элементов конструкций в виде пластин и оболочек с заполнителем, относящимся к классу трансверсально-мягких, и с внешними слоями, имеющими на внешнем контуре подкрепляющие стержни, предназначенные для обеспечения передачи нагрузки на несущие слои при взаимодействии с другими элементами конструкций. Для таких конструкций при малых деформациях и средних перемещениях была построена уточненная геометрически нелинейная теория, позволяющая описать процесс их докритического деформирования и выявить все возможные ФПУ несущих слоев и подкрепляющих элементов. Она основана на введении в рассмотрение в качестве неизвестных контактных усилий взаимодействия внешних слоев с заполнителем, а также внешних слоев и заполнителя с подкрепляющими стержнями во всех точках поверхностей их сопряжения. Для вывода основных уравнений равновесия, статических граничных условий для оболочки и подкрепляющих стержней, а также условий кинематического сопряжения внешних слоев с заполнителем, внешних слоев и заполнителя с подкрепляющими стержнями использовался предложенный ранее [31] обобщенный вариационный принцип Лагранжа. От всех известных вариантов предложенная в [29, 30] теория отличается высокой степенью точности и содержательности при введении в рассмотрение минимального количества неизвестных двумерных функций для оболочки, одномерных функций для подкрепляющих стержней и двумерных и одномерных контактных усилий взаимодействия элементов конструкции.
Согласно результатам, изложенным в [29, 30], для трехслойной оболочки, показанной на рис. 1, а, примыкающий к ней участок (рис. 1, б) при составлении расчётной схемы с большой степенью точности целесообразно рассматривать как подкрепляющий стержень, так как заштрихованный объем между двумя внешними слоями в реальных конструкциях всегда заполняется тем или иным жестким материалом. Существующие методы и известные пакеты прикладных программ для расчёта такого класса трехслойных пластин и оболочек с учётом реальных особенностей их конструктивного исполнения на контуре абсолютно не эффективны и практически не пригодны. В связи с этим настоящая работа посвящена разработке на основе уравнений, предложенных в [29, 30], численного метода решения осесимметричной геометрически нелинейной задачи о докритическом и закритиче-ском поведении трехслойной цилиндрической оболочки, относящейся к описанному классу конструкций, и нелинейной задачи на собственные значения, возникающей при линеаризации в окрестности решения исходной нелинейной задачи о докритическом деформировании, с целью выявления и исследования возможных осесимметричных ФПУ. В среде MatLab с целью численной реализации предложенных приближенных методов разработан комплекс программ и с его помощью проведены численные эксперименты. Отметим, что аналогичные исследования были ранее проведены в работе [32], в которой на базе уравнений [29, 30] была рассмотрена геометрически нелинейная задача о продольно-поперечном изгибе по цилиндрической форме трехслойной пластины с трансверсально-мягким заполнителем, подкрепленной в торцевых сечениях абсолютно твердыми телами, разработан численный метод её решения и проведено исследование докритического напряженно-деформированного состояния (НДС). 1
1. Постановка задачи
Будем считать, что внешние слои трехслойной цилиндрической оболочки имеют толщины 2h(k), а ортотропный трансверсально-мягкий заполнитель - толщину 2h (рис. 1, а). Пусть У С R2, У(к) С R2 - ограниченные области, занимаемые срединными плоскостями заполнителя, нижнего (к =1) и верхнего (к = 2) несущих слоев (здесь и далее индекс в круглых скобках соответствует верхнему и нижнему несущим слоям и к = 1, 2); z G [—h,h] и z(k) G [-h(k), h(k)] - координаты, отсчитываемые от У и У(д.) соответственно в направлении единичной нормали m к У, составляющий с единичными векторами ei и в2 выбранной системы координат правосторонний триэдр (рис. 2).
В работах [29, 30] для трехслойных оболочек описанного класса на основе обобщенного вариационного принципа Лагранжа [31] была предложена уточненная геометрически нелинейная модель деформирования и получены основные уравнения равновесия оболочки, подкрепляющих стержне и условия кинематического сопряжения жестких внешних слоев с заполнителем, а также внешних слоев и заполнителя с подкрепляющими стержнями. Для трехслойной цилиндрической оболочки, для которой параметры Ляме равны Ai = 1, A2 = R, а главные кривизны на У есть ki =0 и k2 = 1/R, описанные уравнения равновесия внешних слоев и кинематические условия их сопряжения с заполнителем будут иметь в осесимметричном случае следующий вид (здесь и далее 5(() = 1, 5(2) = —1):
dN(k + 5(к)Е3 dx 2h
w(2) — w(1)
+ *3
(к)
T(k) T 22
R
0,
(1)
dT
(k)
ii
+ 5(k)q1 = °
(2)
dx
__ ..+
X — x
+
Рис. 2. Трехслойная оболочка с контурными подкрепляющими стержнями
M1
X- «(2) _ H
(l)
W(1) — H(2)W(2) +
2hqi 2h3 d?qi
G13 3E3 dx2
0.
(3)
T
(k)
22 и обобщенные перере-
(k)
В уравнениях (1)—(3) тангенциальные усилия тЦ^
зывающие силы N(k) выражаются через неизвестные осевые перемещения и прогибы w(k) срединных поверхностей несущих слоев и касательные напряжения в заполнителе qi, постоянные по его толщине, по формулам
N((k) _ Q(k + H(k) qi,
Q(k) _ dM[kk Q 1 _ dx
+ IT1(1k)^(k)
M^k _ — D(i)
d2w(k) dx2 '
(k) _ (k)
Til _ B1l
du
(k)
(k) _ (k)
T 22 B22
■®[ ^+2( X)+er
dX + ^
(k)
(k) _
(k)
+ v21
dw(k)
dx
w(k')
R
(4)
2
где R - радиус срединной поверхности заполнителя (см. рис. 1, а), B^ _ _ 2h(k)E(k)/ (1 — v^v^i), i _ 1, 2, и < _ B(1) h2k)/3 - жесткость на растяжение-сжатие и изгибная жесткость к-го слоя, имеющего толщину 2h(k) и выполненного
из материала с модулями упругости E(k), E2k) и коэффициентами Пуассона v(2),
(k) 2 2
v21'; H(k) _ h(k) + h, 2h - толщина заполнителя, имеющего модуль упругости E3
в поперечном направлении и модуль поперечного сдвига G13; M^ - внутренний изгибающий момент в k-м слое; X(k), X3k) - компоненты поверхностной нагрузки, приведенной к срединной поверхности k-го слоя, Q 1k) - перерезывающие силы в k-м несущем слое без учёта касательных напряжений в заполнителе.
Если ввести в рассмотрение [29, 30] неизвестные контактные усилия взаимодействия внешних слоев и заполнителя с подкрепляющими стержнями (см. рис. 1, б, 2), то решения уравнений (1)-(4) при x 1 _ x± должны удовлетворять граничным условиям
M
(k)
11
L
(k) 11 >
dM(k)
dx
+Т
(k):,,(k)
11
+ H(k)q1 _ Q1 3 + hr1,
T(k)
Qli), q1 _ n, (5)
где Q(k, Q(3) и L(k ) - неизвестные погонные контактные усилия и изгибающие моменты, T1 - поверхностные контактные напряжения, формирующиеся в точках сопряжения подкрепляющего стержня с заполнителем в сечениях x 1 _ x± .
Предположим, что при нагружении оболочки осесимметричной нагрузкой в ней формируется НДС, являющееся симметричным относительно сечения х\ = 0. Тогда в этом сечении требуется удовлетворить условиям симметрии решений уравнений (1)—(4) вида
ш(к) = 0, Ni(k) = 0, u(k) =0, q1 = 0. (6)
При действии внешней погонной нагрузки, приложенной в точке A подкрепляющего стержня (рис. 1, а), его равновесие будет описываться следующей системой алгебраических уравнений:
f+ = - (Qii) + Qn) + T+ =0, (7)
fW = - (Q((3) + Q((23) + 2hTi + Bs W - Kz + T+ = 0, (8)
f+
J (p
L k=1
- Sik)Q[kiH(k)) +
+ b2+) (qH + Qi3 + 2hr^ + D
V
W
, Л____K _____
z R2 z R2
+ m+ = 0,
(9)
где [29, 30] Bs
E dQdz, Kz = EQdQdz, Dz = E Z2 dQdz — жест-
F F F
костные характеристики поперечного сечения подкрепляющего стержня, E - некоторый осредненный модуль упругости материалов подкрепляющего стержня в окружном направлении, T+, T+ - компоненты внешней погонной нагрузки T+ = T+n + T+m, приложенной в точке линии г a = Zah. + z^m (рис. 2), m+ = T+ za - T+ Za - погонный изгибающий момент.
Формирующиеся в подкрепляющих стержнях перемещения U, W и угол поворота V должны удовлетворять кинематическим условиям сопряжения стержней с несущими слоями
u(k) - (U - 6{k)H{k)v) =0 при 6Q[ki =0, (10)
w(k) - ^W + B(+)^ = 0 при SqU = 0, (11)
wk + v = 0 при SL(i) = 0, (12)
к которым необходимо добавить кинематическое условие
w(1) + w(2) + qi i - 2 (w + B2+)v) = 0, (13)
имеющее место при Sti = 0. Здесь и далее единицы после запятой в нижнем индексе означают производные по х .
Разрешающая система уравнений об определении вектор-функции неизвестных
X = (w
(1) w(2), U(1), U(2)
q Q(1) q(2) Q(k) Q(1) L(1) L(2) t U W ,q1,Q11 ,Q11 ,Q11 ,Q 13 ,L11 ,L11 ,T1,U, W, V)
состоит из пяти геометрически нелинейных дифференциальных уравнений равновесия цилиндра (1)-(4) с граничными условиями (5), трех линейных алгебраических уравнений равновесия подкрепляющего стержня (7)-(9) и четырёх кинематических условий сопряжения (10)—(13).
2. Сведение краевой задачи к интегро-алгебраическому виду
С целью нахождения численного решения сформулированной краевой задачи (1)—(13) об определении функции X методом конечных сумм [33-35] сведем её к системе интегро-алгебраических уравнений типа Вольтерра второго рода с дополнительными соотношениями для определения неизвестных констант интегрирования. Заметим, что в исходную дифференциальную задачу входят производные порядка 2п от искомых функций, в то время как интегральные уравнения будут содержать лишь п-е производные. Такая редукция сформулированной задачи проводится путём интегрирования уравнений (1)—(4), удовлетворения условиям (5), (6) и использования соотношений
u(k)(x)
du(k)
ds
■ ds,
X
qi(x)
X
0
dw(k)
dx
d2w(k) ds2
£
(k)(x) = w(k)
ds, w( ) (x) =
wa -
d2w(k)
ds2
ds d£,
(k) |
где wa = w|
0 x 0
+ . В результате относительно вектор-функции
xi
X(x) = (w(ll,w(ll,ul1l, и1д, qi,i,wiL>,w^>, Q\1, Qu ,Qi3 ,Qi3 ,Lii ,Lii ,тъ ua, v),
x € [0, x+] = [0, a], приходим к следующей системе интегро-алгебраической уравнений:
( ^ w(2) Q( 1 ) Q(2) Q(У Q(2) T( 1 ) T(2)
X
a
т (k) zi/f(k) | г E3
L11 - M11 16(k)
a £
X0
a
a n
w(a) - w(a1} -
s 0
d2w(2') d2w(1)
£
+ I Tnl ^dsdt + H(k) !f dscJt-
B.
(k)
22
a £
R
X0
V12
du(k) If f d2w(2)
ds I 2 ds2
0
dn
2i
ds 2 ds 2
a £
ds2
a n
dy dn
ds d£+
X0
+ rK?> -
d2w(k)
ds2
s0
dydn f ds d£
= -JJ X3k)dsd£, (14)
X0
Qff - Bk ( ^ I 1 ( / ^ dsf I R - jj ^ dsdi'jY
0 x 0
I S(k) I I ~ds dsd£ = 0, (15)
X 0
2h3 dq1 2h f f dq1
ЗЕ3 dx G13 J J ds
X0
£ 2
I GfH 57 d‘d( +//E - H
X0
k = 1
du(k) d2w(k')
I(k)^r~
ds d£ = 0, (16)
£
a
a
a
[ ^ dx — т1 = 0, J dx
0
(17)
S(k)Es
2h
B(k)
B22
R
(k)
k=i
a £
(k)
d2 w(k) ds2
dsdf
x 0
dx + Ql3) + hTi —
Vl2
du(k) If f d2 w(k)
ds2
ds + 2
dn
a n
+ r(W(ak) —
d2w(k)
ds2
s 0
dydn \ dx =
= — J x(k) dx, (18)
a
a
2
a
которую необходимо дополнить алгебраическими уравнениями равновесия стержня (7)—(9) и кинематическими условиями его сопряжения с оболочкой (10)—(12), записанными в интегральном виде
a
0
du(k) , . dx dx (U s(k)H(k)<f) =° (19)
w(k) — (w+b2+)v) =0, (20)
a , , [ d2w(k) , — dx + p = 0. dx (21)
0
Таким образом, для определения введенных в рассмотрение семнадцати неизвестных из X € Hх5 х Rx11 получена разрешающая система семнадцати интегро-алгебраический уравнений (14)—(21) и алгебраических уравнений (7)—(9), где H = = L2(0,a), a - полудлина пластины. Сформулированная краевая задача может быть записана в операторном виде
A(X) = F. (22)
3. Аппроксимация интегральных уравнений методом коллокаций по гауссовским узлам
Для аппроксимации полученных интегральных уравнений типа Вольтерра будем использовать предложенный в [33, 34] метод коллокаций по гауссовским узлам и способ построения интегрирующих матриц. Введем в рассмотрение интегральные операторы по формулам
x a a
J (f) = У f (f) df, J*(f) = y f (f) df, J3(f) = у f (f) df,
0 x 0
которые аппроксимируем конечномерными аналогами в виде интегрирующих матриц Ji, J2 , J3 соответственно (здесь и далее прямым шрифтом обозначены операторы, действующие в пространстве измеримых функций, готическим шрифтом -операторы, действующие в пространстве сеточных функций).
С этой целью на отрезке [0, а] введем сетку ш = {xi, i = 1, 2, ..., N}, соответствующую квадратурной формуле Гаусса
N
Mf) = 53 dif (Xi),
i=1
где {di}, {xi} - соответственно веса и узлы коллокаций, связанные с корнями полинома Лежандра степени N.
Обозначим через fi значение f в узле xi : fi = f (xi) и приблизим f на отрезке
N
[0, а] посредством интерполирующей функции f (x) » 53 fi li(x). В качестве базис-
i=1
ных выберем функции Лагранжа {li} по узлам {xj}. Таким образом, раскладывая функции li по полиномам Лежандра, строим интегрирующие матрицы Ji, J2. Введем в рассмотрение конечномерные операторы
Т = (ТШ Т (2) Т(1) Т (2) Т(1>) : H X, Jx12 H
Т1 = (I 1w ,т 1w ,Т 1и ,Т 1и ,Т 1q ): Hh Х R ^
X5 h ,
X5,
h,
(1) (2) (1) (2) XN X5 X12
Т2 = (Т 2w ,Т 2w, Т 2и ,Т 2и , U ): Hh Х R ^
Г = (Гт, rwi(), Г(2а) : HX5 Х RX12 ^ RX3,
Q = Q{\\ Q121), Q113), Q13, l11), L™) : H X5 Х RX12 ^ RX6 T = (Ms,Tz,TZ) : HX5 Х RX12 ^ RX2
по формулам
IЦХ (x) = L(1k) + D^wf, + H(k) J2 J1qn +
+ S(k)J2J1 Eh |w(2) - w(1) - J2J1 (w™ - wg
B.
(k)
22 т т
RJ2J1
>(k)v 22
V12U(k) + R (wik) - J2 J^v^l
(k)
(k)
I^X(x) = J2TRJ1wf, - -2R2 J2J1 (jWS) ;
+ S(k)J2J1q1,1;
(k) (k) (k)
1 1u X (x) = Q11 B11
E + R (w(k) - J2J1Wf
,1
T~) (k) 2
I (2kU)X (x) = - -R (J1wf) ;
11qX(x) = (JE + GG^J2J^ q11 + J2J1 - u<1 - H(1)w(11 - H(2)w(n) ;
ГтX(x) = J3q1,1 - T1;
гЦX(x) = 2)Тз{Ез [w™ - wR - J2J1 (wg - w^)] } +
(k)
B(k)
+ Q13 + hT1-----JR J3
TD (k) 2
r2klX(x) = --RJ3 (j^f) ;
+ (wR - J2J1w
(k)
V12U1 + r \ wa' - J2J1w11
q[\)X(x) — Xuf - (U - 5{k)H{k)v) ; Q^X(x) = wik) - (W + B(+)^ ; L[\)X (x) — Jw™ + p;
MsX (x)
E (4?
8{k)QuHW) + B(+) ( Q(13) + Q(123) +2hn
+ Dx
+ R2
Tz x (x) — Q13 + Qk3 + 2hTi +---—2—
Kz P
R2 ;
Kz W
R2 ;
Tc X (x) — Q(111) + Q™.
Тогда сформулированная геометрически нелинейная конечномерная задача может быть представлена в матричной форме
Ai X + A2(X) — F, (23)
где Ai — (I1, r,Q,T) : H£4 5 * * * x Rx12 ^ H£5 x Rx12 - линейный оператор, A2 — — (I2, 0x12) : Hx x Rx12 ^ Hx x Rx12 - нелинейный оператор; в общем случае правая часть F(T+,T+,m+ ,x3t^) принадлежит Hx5 x Rx12. Для задачи об осевом сжатии цилиндрической оболочки с эксцентриситетом для соответствующего компонента внешней погонной нагрузки и погонного момента, приложенных к подкрепляющему стержню, и для поперечной распределенной нагрузки на цилиндр положим TZ — m+ — X|k) — 0, тогда F — (Ox(5N+11),T+).
Таким образом, для функций, входящих в систему дифференциальных уравнений (1)-(3) со старшей производной порядка 2n, построена конечномерная схема (23) относительно n-й производной решения краевой задачи. После решения конечномерной задачи (23) решение исходной краевой задачи восстанавливается численным интегрированием при помощи полученных ранее интегрирующих матриц.
4. Итерационный метод и численные эксперименты
Для решения геометрически нелинейной проекционной схемы (23) будем ис-
пользовать следующий двухслойный итерационный процесс с опусканием нели-
нейности на нижний слой [36-39] c предобуславливателем, являющимся линейной
частью оператора разностной схемы
X (n+1) -A1----------
T
X1n)
+ (A1 + A2) X1n)
F,
(24)
где X(0) - заданное начальное приближение, t > 0 - итерационный параметр.
Докритическое НДС конструкции. Для численной реализации итерационного метода (24) решения задачи (23) об осевом сжатии оболочки силой T+, приложенной к подкрепляющему стержню с эксцентриситетом, разработан комплекс программ в среде Matlab. На его основе проводились численные расчёты при следующих значениях геометрических и упругих параметров оболочки и подкрепляющего стержня (см. рис. 1):
x+ — a — 50 см, 2h(1) — 2h(2) — 0.1 см, h — 1 см, G13 — 25 МПа,
E{k) — 103 • 103 МПа, E(2k) — 40 • 103 МПа, Е3 — 50 МПа, E — 20 • 103 МПа,
VV2 — 41 — 0.04, к — 1, 2, R — 52.5 см, za — Za — B(+) — h + h1.
Рис. 3. Прогибы точек срединных поверхностей несущих слоев w(k), см
X
Рис. 4. Осевые перемещения точек срединных поверхностей несущих слоев u(k), см
При выбранных значениях осредненного модуля E и геометрических параметров подкрепляющего стержня для жесткостных параметров Bs, Kz, Dz приняты значения Bs = 1.9708 • 107 Н, Kz = —61300 Н • м, Dz = 727.6401 Н • м2 . Число точек сетки N = 256, что обеспечивает приемлемую точность при использовании предложенного численного метода по всем параметрам НДС конструкции. Вычисления согласно (24) проводились до тех пор, пока норма невязки || F — (Ai + A2) X(n) || и разницы между итерациями Ц X(n+1) — X(n) ^ оставались больше заданной точности е = 5• 10—8 . В качестве нормы вектора g = (gi, g2,..., gm) выбиралась величина ||g|| = max { |gi\, |g2\,..., \gm\ }. Итерационный параметр подбирался эмпирически.
При заданном значении осевой сжимающей погонной нагрузки T+ = —50 кН/м в подкрепляющем стержне, а также в сечении его сопряжения с несущими слоями и заполнителем формируются параметры НДС, имеющие следующие значения:
Q(11) = —1.371 кН/м, q121) = —48.629 кН/м, Q(113) = —0.53221 кН/м,
Q(123) = —0.65079 кН/м, ьЦ = —4.8743 Н, ьЦ = —4.6653 Н,
U = —0.01209 см, W = 0.0092828 см, р = —0.0045715, т1 = —0.0375 МПа.
Результаты расчётов параметров НДС оболочки приведены на рис. 3-10.
На рис. 3 показаны графики изменения по длине оболочки прогибов точек срединных поверхностей несущих слоев. Видно, что они практически совпадают в силу малости деформаций поперечного обжатия заполнителя в докритическом состоянии конструкции. В окрестности соединения внешних слоев с подкрепляющим стержнем наблюдается сильная изменяемость функций прогибов, в силу условия р = 0 являющаяся следствием деформации выворачивания подкрепляющего стержня. Заметим, что прогиб верхнего несущего слоя в сечении Ж1 = a есть w(k)\x=x+ = W + b2+') см, в то время как прогиб точки O+ подкрепляющего стержня принимает значение W = 0.0010075 см.
На рис. 4 изображены осевые перемещения точек срединных поверхностей несущих слоев. Заметим, что в верхнем несущем слое их значения значительно больше, чем у нижнего слоя. Такая разница обусловлена передачей внешней сжимающей нагрузки главным образом на верхний несущий слой оболочки (Q Ц) = (2)
= —1.371 кН/м, Q11 = —48.629 кН/м в сечении Ж1 = a) в силу особенностей конструктивного исполнения оболочки и действия на неё внешней силы с эксцентриситетом и поворотом поперечного сечения подкрепляющего стержня против часовой стрелки на величину \р\. Из рис. 5 следует, что оболочка в целом находится в моментном состоянии (см. рис. 6) ввиду действия на неё сжимающей
о
-10 -20 -30 -40 -50
0 10 20 30 40 50
x
Рис. 5. Осевые мембранные нормальные Рис. 6. Изгибающие моменты в несущих
(к)
напряжения в несущих слоях ah
= T('l)/(2h(k)), МПа
слоях , Н
Рис. 7. Поперечные касательные напряжения в заполнителе ql, МПа
Рис. 8. Окружные мембранные напряжения в несущих слоях а(2) = T2(k)/(2h(k)), МПа
Рис. 9. Обобщенные перерезывающие силы Рис. 10. Перерезывающие силы в несущих
N(k\ кН/м
слоях Q(\ кН/м
нагрузки с эксцентриситетом, причем формирующееся в нижнем слое мембранное
(1)
нормальное напряжение аЦ в сечении Ж1 = а практически равно нулю, а значение тЦ + 2h.(2)CTn) в данном сечении - сжимающей нагрузке T((+) = —50 кН/м.
Рис. 11. Деформированное состояние срединной поверхности верхнего несущего слоя при Е = Е(2к) = 20 • 103 МПа
-50 0 50
Рис. 12. Деформированное состояние срединной поверхности верхнего несущего слоя при E = Е(2к) = 80 • 103 МПа
На рис. 7 показано изменение касательных напряжений в заполнителе. Видно, что в рассматриваемом случае нагружения и в силу особенностей конструктивного исполнения оболочки на её контуре они имеют большие градиенты в окрестности сечения х1 = а (в точках поверхности сопряжения заполнителя с подкрепляющим стержнем они должны принимать нулевое значение в случае отсутствия адгезионного соединения), где достигают максимального по модулю значения.
На рис. 8 показаны функции мембранных окружных напряжений в несущих слоях. Видно, что в окрестности сечения xi = а они являются положительными из-за значительных положительных радиальных перемещений точек осевой линии подкрепляющего стержня. По мере удаления от данного сечения эти напряжения переходят в сжимающие вследствие деформации выворачивания подкрепляющего стержня (см. рис. 8). В силу таких особенностей деформирования, судя по результатам, приведенным на рис. 9, 10, максимальное участие внешних слоев в восприятии поперечных касательных напряжений наблюдается лишь в окрестности сечения х1 = а.
С целью иллюстрации на рис. 11 изображен деформированный цилиндр в соответствии с найденными прогибами точек срединной поверхности верхнего несущего
Рис. 13. Зависимость погонной нагрузки Т^+^ от максимального по модулю значения прогиба w(2) верхнего слоя
слоя, где цвет характеризует значение мембранных нормальных напряжений в этом слое. Видно, что края цилиндра растягиваются в радиальном направлении. Это связано с тем, что модуль упругости в окружном направлении подкрепляющего стержня в два раза меньше, чем несущих слоев оболочки (E = E= 20-103 МПа).
На рис. 12 изображено деформированное состояние срединной поверхности верхнего несущего слоя при E = E= 60 • 103 МПа. Видно, что увеличение модуля упругости подкрепляющего стержня в окружном направлении в четыре раза приводит к значительному уменьшению деформации его выворачивания и, как следствие, к уменьшению окружных деформаций внешних слоев оболочки.
НДС конструкции в окрестности бифуркационного значения внешней нагрузки. На основе изложенного алгоритма путём табулирования по параметру внешней нагрузки отслеживался закон изменения внешней нагрузки в зави-
симости от максимального значения прогиба верхнего несущего слоя в окрестности места возможной потери устойчивости. Установленный закон в виде зависимости значения нагрузки Т^ от максимального по модулю значения прогиба w(2) верхнего слоя приведен на рис. 13. В окрестности некоторого бифуркационного значения нагрузки Т+ итерационный процесс (24) перестаёт сходиться. Результаты расчётов, получаемые до такого значения внешней нагрузки, дают информацию только о возможных устойчивых положениях равновесия. Другими словами, решая поставленную геометрически нелинейную задачу (23) с помощью итерационного метода (24), при табулировании по параметру погонной нагрузки Т^+) мы можем найти только докритические положения равновесия. Эксперименты показали, что итерационный процесс перестал сходиться при Т^+) = —460 кН/м и формировании следующих параметров НДС в сечении xi = a:
qW = —8.4973 кН/м, Qf! = —451.5027 кН/м, Q^ = 3.4434 кН/м,
Q(3 = —29.3836 кН/м, Ь{\] = 2.0716 Н, L^ = 38.732 Н,
U = —0.11963 см, W = 0.10614 см, р = —0.049425, п = 0.18416 МПа.
Законы изменения параметров НДС оболочки, найденные при указанном выше значении внешней нагрузки, приведены на рис. 14-17.
Сравнивая результаты, приведенные на рис. 3, 5, 7, 8, с результатами, приведенными на рис. 14-17, можно видеть как количественные, так и качественные изменения параметров НДС рассматриваемой конструкции по мере приближения значения внешней нагрузки к бифуркационному.
Рис. 14. Прогибы точек срединных по- Рис. 15. Поперечные касательные напря-верхностей несущих слоев w(k), см жения в заполнителе qi, МПа
Рис. 16. Осевые мембранные нормаль-
ад
ные напряжения в несущих слоях аЦ = = /(2h(k)), МПа 5
Рис. 17. Окружные мембранные напряжения в несущих слоях а(2) = Tk /(2h(k)), МПа
5. Линеаризация задачи в окрестности её нелинейного решения и постановка задачи на собственные значения
Главной целью настоящей работы является нахождение точки бифуркации задачи (22), определение соответствующего ей критического значения внешней нагрузки и форм потери устойчивости конструкции. В предлагаемом методе численного решения рассматриваемой задачи (22) точка бифуркации определяется как точка ветвления её решений. Один из наиболее употребительных способов решения указанной задачи поиска точки бифуркации состоит в линеаризации задачи в окрестности решения линейной задачи. В настоящей же работе предлагается проводить линеаризацию задачи (22) в окрестностях её нелинейных решений.
Основная идея определения точки бифуркации на основе решений однородных линеаризованных уравнений состоит в следующем. Предположим, что одна какая-то форма равновесия системы известна и нужно найти точку бифуркации этой формы равновесия. Для этого, не интересуясь поведением системы вдали от известной формы равновесия
X = (w
(1) w(2\u(1\u
(2) q Q(1) Q(2) Q(k) Q(1) L(1) Г(2) T TT W ,л) ,q1, Q11 ,Q11 ,Q11 13 , L11 , L11 , T1, UW
достаточно найти условия существования другой формы X + SX, отличной от исходной, но бесконечно к ней близкой
A(X + SX) = F,
(25)
Другими словами, точка, в окрестности которой существуют две формы равновесия X и SX (то есть A(X) = АДХ) = F), и является точкой бифуркации, где
SX = (Sw(i), Sw(2), Su(i), Su(2), Sqb SQ[\\ SqI) SQ[\\
SQ[l\ SL[\\ SL(2l\ Sti, SU, SW, Sy>)
есть малое приращение решения.
Сформулированная задача эквивалентна задаче о нахождении такой точки X, что дифференциал Фреше (см. [40]) D(X)SX оператора A(X) в точке X будет равен нулю: D(X)SX = 0. Покажем, что это действительно так, и получим вид линейного по SX оператора D(X) .
Имеем, что при фиксированном X и всех SX
A(X + SX) — A(X) = D(X)SX + R(X, SX), (26)
где D(X) - линейный непрерывный по SX оператор, lim R(X,SX)/||SX|| = 0,
||$X||—,0
D(X)SX - дифференциал Фреше оператора А в точке SX, а R(X, SX) - остаточный член этого дифференциала.
Следует также отметить, что задача о нахождении точки бифуркации, которая получается в результате выполнения соотношений (22) и (25) для одной и той же правой части, будет иметь вид
A(X + SX) - A(X) = 0, (27)
В силу равенства (26) и соотношения R(X, SX) ^ 0 задача (26) запишется в виде
D(X)SX = 0. (28)
Вместе с тем оператор А представим в виде А = Ai + А2, где Ai - линейный, А2 - нелинейный операторы. Тогда задача (26) после простых преобразований будет иметь вид Ai(SX)+A2(X+SX)- A2(X) = 0, где A2(X+SX) — A2(X) является дифференциалом Фреше оператора А2 в точке X. Обозначим его через D2(X)SX. Тогда задача (28) может быть сформулирована в виде Ai(SX) + D2(X)SX = 0 или
D(X)SX =(Ai + D2(X)) SX)=0, (29)
где D(X), D2(X) - линейные непрерывные по SX и нелинейные по X операторы. Применительно к сформулированной задаче (1)-(3), (5)—(13) задача (29) будет описываться следующими дифференциальными уравнениями:
(SN[k))A + (3 — 2k)c3 (dw(2) — Sw(i)) = 0, (30)
(ST1(k)),i + (3 — 2k) сз Sqi =0, (31)
Su(i) — Su(2) — H[k) (Sw(i)),i — H2k)(Sw(2))A +
+ 2hG-3Sqi — 2h(3Gi3) i(Sqi)jii = 0, (32)
SN[k)
D{ki)(Sw(k))Mi + B{k^w<:k)(Su(k)),i + uf(Sw
(k)
), i+
+ 1.5( w|k^2(Sw(k)), i + 2w(1k)(Sw(k)), i) + H(k)Sq
1
0,
S^f = Б{к) ((5п^)л + wf (Sw(k))A) , для которых при x = x± формулируются граничные условия
SMlf = sLf, Sqi - Sti, STlf = SQf, SN[k) = SQ(f + hSru (33)
а при x = 0 - условия
Su(k) = 0, SNf = 0, Su(k) = 0, Sqi = 0. (34)
Поскольку уравнения равновесия подкрепляющего стержня (7)-(9) и кинематические условия его сопряжения с несущими слоями (10)—(13) являются линейными алгебраическими уравнениями, их вид в задаче (29) останется прежним относительно вектор-функции неизвестных SX и они не будут зависеть от X. Задача (29), порожденная системой уравнений (30)—(34), (7)—(13) относительно SX, является однородной задачей с тривиальным решением. Поэтому она может быть сведена к эквивалентной спектральной задаче, а способ её решения будет основан на решении задачи на собственные значения. Для этого рассмотрим вспомогательную задачу
D^TX+^j SX = 0, (35)
которая получается в результате умножения в уравнении (29) X на Л/Т+. Здесь Л — число, характеризующее точку бифуркации в окрестности известной формы равновесия
X = (w
(1) w(2),u(1),u
(2) q Q(1) Q(2) Q(k) Q(1) Г(1) Г(2) t TT W ,o) ,q1,Q11 ,Q11 ,Q11 ,Q 13 , L11 , L11 ,T1,U, W,lf),
T+ — погонная нагрузка, соответствующая состоянию равновесия X. Заметим, что задача (29), сформулированная в виде (35), является задачей на собственные значения относительно пары Л , SX .
Итак, путём последовательного нагружения и решения геометрически нелинейной задачи (22), сформулированной в виде A(X) = F(Т+), табулируя по нагрузке Т+, находим множество форм равновесия {Xt} Э X. Решая при этом для каждой формы равновесия X (то есть на каждом шаге нагружения) задачу на собственные значения (35), необходимо следить за поведением и изменением минимального положительного собственного числа Л. В сформулированной таким образом задаче на собственные значения (35) условие существования другого состояния равновесия будет выполнено лишь в том случае, когда Л становится равным параметру нагружения Т+: Л = Т+, так как именно в этом случае выполняется равенство A(X + SX) = F(Т+) = F(Т+), а найденное значение Л будет соответствовать точке бифуркации. Выполнение условия Л = Т+ при решении задачи (35) эквивалентно решению задачи (29), а найденное значение Л является значением критической нагрузки, которое характеризует точку пересечения двух решений X и SX.
Представим левую часть задачи (35) в виде
4Т+) sx
А1 + ЛЛ1( Т+)) + л2л2(Т+) ISX,
где А1 — линейная часть оператора А задачи (22). Тогда задача (35) сводится к следующей квадратичной задаче на собственные значения для отыскания точки бифуркации:
^А1 + ЛЛ1
X
Т+
+ Л М Т+
SX = 0,
(36)
Для операторов Лх(X/T+ и м X т+ в соответствии c (30)-(34) с помощью используемого метода конечных сумм и вышеизложенного способа аппроксимации интегральных операторов строятся их конечномерные аналоги
Л1 = (K,Ow,5n+i2) : Hf5 х R*12 ^ Hf5 х Rf12, Л2 = (L,Ow,5n+i2) : Hf x R*12 ^ Hf х R*12,
в которых
K
X
T+
SX
( k^sx \
K^SX
kU1sx
K{2)SX ,
On,5n+12SX
kWHsx
V kWIsx
L
X
T+
SX
( L&SX I
lW)sx I,
\O3n+2,5n+12SXy
K(k)(X )SX = J2Bi1k1>(U'k> + v21R-1W(k))J1(SW(k)),1 + J2Bi1k1>((Su(k) ),1+
+ v21R-1(Swk - J2 J1(Sw(k))M))»! - в^^-1^^.]^»^)^
K(k)(X )SX = -B^wf-h (Sw(k)),u,
Kik}(X )SX = -B{2k2) v12R-1JstifJ1w<:k:l L(k)(X )SX = J2B[k1)l.5(S(k:))2J1(Sw(k)),11, X = (■U(1),-U(2),UM),U(1),U1), W(k) = w(k)/T+, U(k) = u(k)/T+, U1 = Q1/T+.
6. Метод продолжения по параметру — работе внешних сил
В настоящей работе центральное место занимает разработка метода определения полей напряжений, деформаций и перемещений в окрестности бифуркационного значения нагрузки в докритическом и закритическом положениях равновесия. Как было отмечено выше, при проведении расчётов путём табулирования по нагрузке T+ итерационный процесс (24) в окрестности некоторого значения T+, по косвенным причинам являющимся бифуркационным (рис. 13), перестаёт сходиться, поскольку, табулируя по указанному параметру нагружения, мы не можем выйти на ветвь решения, где приращение функции этого параметра меняет знак. Поэтому для поиска неустойчивых положений равновесия и исследования поведения оболочки в таких положениях предлагается алгоритм, основанный на глобально инкрементальной теории Лагранжа. Он является одним из вариантов реализации инкрементального алгоритма процесса продолжения решения по параметру нагрузки [41, 42], в соответствии с которым процесс деформирования представляется в виде реализации последовательности равновесных состояний при соответствующих уровнях нагружения. При том значении нагрузки, при котором итерационный процесс (24) перестает сходиться, в качестве нового параметра нагружения предлагается выбрать либо работу внешних сил At в силу её строгого возрастания, либо параметр кинематического нагружения U*. Вследствие этого требуется формулировка новых геометрически нелинейных задач относительно
уже функций приращений искомых неизвестных исходя из (1)—(3), (5)—(13). В случае, когда параметром продолжения является Ат , к системе дифференциальных уравнений (1)—(3), (5)—(13) добавляется дополнительное уравнение
At = T+U (37)
относительно нового дополнительного неизвестного T+ .
Таким образом, геометрически нелинейная задача о равновесии трехслойной цилиндрической оболочки, в которой параметром нагружения выбрана работа внешних сил At , относительно вектор-функции неизвестных
(T+
w(1), w(2),
,(1) 1 )
q 1, Q1V, Q 121), Q13 ,Q 123), L11, L1V, T1, U, W, ф
формулируется на основе дифференциальных уравнений равновесия оболочки (1)— (3) с краевыми условиями (5), (6), уравнений равновесия подкрепляющего стержня (7)—(9), кинематических условий сопряжения оболочки с подкрепляющим стержнем (10)—(13) и дополнительного уравнения (37), представляющего собой выражение для работы внешних сил. Заметим, что при её решении табулирование возможно как по работе внешних сил At , так и по параметру кинематического нагружения U*.
Для реализации описанной схемы метода продолжения по параметру для устойчивого равновесного положения конструкции в её докритическом состоянии, найденного путём решения задачи (23) с помощью итерационного процесса (24), на последнем шаге нагружения зафиксируем найденное решение
X = (T+ ,w(1),w(2)
,а) u(1)
q Q(1) Q(2) Q(1) Q(2) L(1) L(2) t U W ф) ,q1, Q11 , Q11 ,Q 13 ,Q 13 ,L11 ,L11 ,r1, U, W ,ф)-
(38)
Введем далее в рассмотрение оператор C : Hx5 х Rx 13 ^ R по формуле C(X) = T+U, где T+ - новая неизвестная скалярная величина. Тогда сформулированную геометрически нелинейную задачу (1)-(3), (5)—(13), (37) представим в следующем операторном виде:
B(X) = B1(X) + B2(X) = F, (39)
где B1(A1, 0) : Hх5 х Rx13 ^ Hх5 х Rx13 — линейный оператор, B2(A2, C) : Hх5 х х Rx 13 ^ Hx5 х Rx 13 — нелинейный оператор; вектор-функция правых частей F зависит только от скалярной величины работы внешних сил At .
Вектор-функцию неизвестных в неустойчивых положениях равновесия для некоторого значения работы внешних сил, равного Ат = At + ДАт, обозначим как
(T+,
?(1),W3(2),1(1),1(1),i1,Q?C1i), Q1V, Q^ Qg\ L[\\ L121),h, U, W, ф), (40)
где тильда над элементом р означает добавление приращения: p = p + Др. Тогда наряду с уравнением равновесия (39) относительно X можно записать еще одно уравнение относительно положения равновесия X
B(X) = B1(X) + B2(X )B1(X) + Bl(ДX) + B2(X) + B2(ДX) = F + ДF. (41)
Вычитая из уравнения (41) уравнение (39), получим следующую геометрически нелинейную задачу относительно вектор-функции приращений
Имеем, что B1(AX) + B2(X + AX) — B2(X) = Da(X)AX + R(X, AX), где Da(X)AX — линейный по AX и нелинейный по X оператор, R(X, AX) — нелинейный оператор.
Получим явный вид этих операторов. Подставив (40) в уравнения (1)—(3), приходим к системе дифференциальных уравнений
- D
(k)
11
d4w(k) d4 Aw(k) \
dx4
+
dx4
)+H(k){ ix + +(3 - 2k)c^(
,(2)
— w(1))+
+ (Aw(2) — Aw(1))) + X{3k) + Ax(k) + dx\ b(1)
du1k dAu(k)
+ —;----+
dx
dx
1 / dw(k) dAw(k) \2 (w(k) + Aw(k))
+ 2 V dx + dx ) +
dx
dx
R
(dw(k) dAw(k)\
dx dx
B(k)
B22 J ,,(k)
R
v12
du1k + dAu(k) 1 (dw(k) dAw(k)\
+ 2 у dx + dx J
2 i
dx
+
dx
+
w (k) +Aw(k) )
R J
d
dx
^ B
(k)
11
du1k dAu(k) 1 (dw1k dAw(k) \
dx
+
+
2i
+
+
dx 2 dx dx
+ X{1k) AX(k) + (3 — 2k)(q1 + Aq1) = 0,
' dw(k) dAw(k)\
H(k)[
k=1
(u(1) — u(2)) + (Au(1) — Au(2)) — ± hJ + djwL\
dx dx
+
2h 2h3 (d2q1 d2Aq1
+ G3(q1 + Aq1) — ЖЛ d? +
0
и выражению (T+ + AT+ )(U + AU) = At + AAt для работы внешних сил.
После простых преобразований и при учёте уравнений (1)—(3), которым удовлетворяют функции (38), приходим к геометрически нелинейной задаче (42), описывающейся системой пяти обыкновенных дифференциальных уравнений
dAN-
(k)
^ + Ax3k) +
dx
+ d (b<?
dx
dAu(k) dw1k dAw(k) 1(dAw(k) Д 2 и2ЦAw(k)
dx
+
dx
dx
+ n
2 dx
2
-) +
R
dw(k)
dx
+
+ b(1)
dw(k) dw(k) dAw(k) 1 / dw(k) у v21)Aw(k)
’ + ’ dx + 2 у dx J + R
dx + dx
dAw(k)
dx
+ csS(k)( Aw(2) — Aw(1)) —
B(k)
B22 Jv (k)
~R^12
dAu(k) dw1k dAw(k)
dx
+
dx dx
1
+ 2
fd2_Aw(k) dx2
ds
2i
a n
+R( Awak)—
d2Aw1k) dx2
dydn
(k)
d
+ Ax1k) + (3 — 2k)Aq1 + AL [BY1
dx dx
(k)
s 0
dw(k) dAw1k \
dx
dx
■
.. (i) . (2v tt dAw(k) 2h 2h3 d2Aq1
(Au(1) — Au(2)) —22 H(k) ——— + — Aq1 — —--------— = 0
k=1
dx
G13
ЗЕ3
+
- +
0, (43)
(44)
(45)
0
0
и алгебраическим уравнением
AT+ (U + AU) + T+AU = AAT (46)
относительно вектор-функции приращений
AX = (AT+, Aw(1), Aw(2), Au(1), Au(1), Aqb AQ(111), AQ(121),
AQ&, Aq12s), Ab[\\ AL(121), An, AU, AW, Aip),
зависящей от решения X. Здесь AT+ - новая неизвестная величина, а AAt — входной параметр задачи. Подставив (40)) в граничные условия (5) и (6), для системы уравнений (43)-(45) получим граничные условия при x = x±
(k) i D(k)
AN(k)+ B
11
dAu(k) dw(k) dAw(k)
+
dx
dx
dx
+ b(1)
dw(k)
dx
+
2 у dx
dw(k) dAw(k\ 1 ( dw(k\2: v^k_)w(k)
+ 1 ( dAw(k)\2+ v(k)Aw(k)
R
dw(k)
dx
■ +
dx
AM
(k)
li
AL'
(k), 11 ,
dx
Aq1
■ + -
2 dx At1
+
R
dAw(k)
AT1(1k) + B(1k1)
dx
dw(k) dAw(k)
= AQn + hAT1, (47)
dx
dx
AQ(1k1)
и при x = 0
AN1(k)+ ВЦ)
dAu(k) dw(k) dAw(k) 1
dx dx dx 2
dAw(k) A 2 Aw(k) dw(k)
dx R dx
+B
(k)
11
dw(k) dw(k) dAw(k) 1
dx dx dx 2
dw(k)\2 v^^w(k') dAw(k)
dx R dx
0, (48)
Au(k) = 0, Au(k) = 0, Aq1 = 0.
К составленной системе уравнений (43)-(46) необходимо добавить уравнения равновесия подкрепляющего стержня (7)-(9) и кинематические условия сопряжения подкрепляющего стержня с несущими слоями (10)—(13). Поскольку они являются линейными, вид этих уравнений относительно вектор-функции перемещений AX задачи (42) останется прежним и не будет зависеть от X.
Для приближенного решения сформулированной задачи (42) об определении вектор-функции приращений AX так же, как в разд. 1, использовался метод конечных сумм. В соответствии с этим методом дифференциальные уравнения равновесия оболочки (43)—(45) путём последовательного удовлетворения граничным условиям (47), (48) сведены к системе интегро-алгебраических уравнений
AL(1k1) + бЦ
d2w(k)
dx2
+ Л E1
+ V(k)^ZT
a £
2h
x 0
a n
Aw(2) — Aw(1) —
s 0
f d2Aw(2) 1 dx,2
* dydn + d.df + /[ (aT™ + Bk wLj d.)
dx dx2
0
x wi + (T(k1)+B k /d2 Aw(k) d8\ fd2w(k) AT (k) fd Aw(fc)
dx
dx dx2
0
dx2
-ds+Aт1k;)
dx2
ds
0
0
ds
df-
a £
a £
- Rff (A£’ + ВйЧ? ^f idAr- d")d-8d( = "/<49>
x0
x0
0
AQn -
AT
(k)
ii
+ в(!)
dw(k) f d2Aw(k) dx dx2
a £
ds \ + S,
(k)
d2
qi
dx2
ds d£ = 0,
X 0
(50)
2h3 dAqi 2h
ЗЕ3 dx + G13
a £
X0
d2qi
dx2
a £ 2
w7 £
X 0 k=i
Л dAu(k) d2Aw(k)
ivS(k) ^x Я(к) “d2-
ds d£
0, (51)
S
Е3
(k) 2h
d2Aw(2)
dx2
d2 Aw(i) dx2
ds^dt + Qk
+ hri-
1
R
a
0
(k) (k) (k
At22 + B22 vi2
(k) (k)
dw(k) f d2 Aw(k)
dx
dx2
d£ dx
a
j AX(3k)dx, 0
a
0
dqi
—— dx — Ti
dx
0,
(52)
(53)
dAu(k)
dx
dx — (AU — S(k)H(k)A<f) = 0,
Aw(k) — (AW + B(+)A^>) = 0,
a . .
( d2Aw(k)
dx2
■ dx + A^ = 0,
Af+ = — (AQii + AQii)) + AT+ = 0,
a
(54)
(55)
(56)
(57)
Af+
AL(1i) + AL12) — AQ(1i)H(i) + AQ(121)H(2) +
+ (AQi3 + AQi3 + 2hATi
AT+(U + AU) + AT+AU = AAT,
+ AT+ = 0,
(58)
(59)
где
AT
(k)
11
B
(k)
11
dAu(k) 1
dx +2
d d.) 2 + fl Aw?> —
a n
d2
qi
dx2
dy dy
s 0
AT
(k)_B(k)
22
22
■ {k)( dAu(k) + 1
ViA dd +2
£
d2Aw(k) dx2
ds
a n
+ R ( Aw(k) —
d2
qi
dx2
dy dy
s0
2
Решая сформулированную задачу (49)—(59), определим вектор-функцию приращений ДХ. Это позволяет в неустойчивых положениях равновесия оболочки определить решение X = X + ДХ, а следовательно, и все параметры закритиче-ского НДС.
Заметим, что применительно к сформулированной задаче явный вид введенных ранее в (42) операторов будет следующим:
Da(X)AU = (DA(X)Ax) : Н^ Х R^ ^ Н^ Х RXl3’
Ca(X)AX = AT+U + T+AU,
r(x, ax)
^2(AX)\
e
+
/rW)(x, ax)\
rW\x, ax)
exsn+ii
\ra(x, a x)J
: Hx5 x Rx13 -► Hx5 x Rx13,
RWk)(x, ax) = J2 B{k)-
3 dw(k) ( T d2Aw(k)Y
2 dx
J
dx2 у
, ra(x, ax) = t+au.
Формулы для операторов Da(X)AX и A2(AX) приведены в разд. 2 и 5.
В конечномерной операторной постановке сформулированная задача есть
Da (X)AX + R(X, AX) = AF,
(60)
где
AX = (AT+, Awg, Aw^, Au(1), Au(1), AqM, AQ Ц, AQ^,
AQ(3), AQ(3), AL(1), AL(2), An, AU, AW, A?)
является конечномерным аналогом вектор-функции неизвестных, X - решением задачи (23) о докритическом положении равновесия оболочки. Заметим, что вектор-функция правых частей AF зависит только от скалярной величины - приращения работы внешних сил AAt .
Поскольку оператор R(X, AX) является нелинейным относительно искомой неизвестной AX , для решения задачи (60) предложен итерационный процесс
AX(j+1 ) - AX(j) , s , s
Da(X)-----------------+ Da(X)AX(n) + R(X, AX(n)) = AF,
T
аналогичный (24).
Таким образом, табулируя сначала по нагрузке T+, приложенной в осевом направлении к подкрепляющему стержню, а затем по приращению работы внешних сил AAt , можно установить закон изменения внешней нагрузки T+ в зависимости от максимального значения прогиба верхнего несущего слоя w(2), показанной на рис. 18.
Реализованный метод продолжения решения по параметру позволяет проследить, как ведет себя оболочка в закритическом состоянии. На рис. 18-23 приведены результаты таких исследований для оболочки с указанными ранее параметрами. Их анализ показывает, что значение Т£ = 463.22 кН/м является критическим, так как дальнейшее увеличение прогибов происходит при уменьшающемся значении Т+ (см. рис. 18). При таком значении действующей нагрузки наиболее интенсивное волнообразование наблюдается у верхнего несущего слоя в области крепления подкрепляющего стержня с оболочкой (рис. 19). С целью иллюстрации на рис. 2023 для одного и того же значения нагрузки Т+ = 450 кН/м приведены графики
500
400
Е 300 z
+н^200
100
0
— -Геометрически нелинейное НДС ----Метод продолжения по параметру
0 0.1 0.2 0.3 0.4 0.5
||w(2)-w(2)(a)||-||w(2)(a)||
Рис. 18. Зависимость погонной нагрузки T+ от максимального по модулю значения прогиба w(2) верхнего слоя
Рис. 19. Деформированное состояние срединной поверхности верхнего несущего слоя в точке бифуркации
изменения прогиба верхнего слоя w(2), формирующихся в нем мембранных на-
- (1) (2)
пряжении а11, а22 и поперечного касательного напряжения в заполнителе q1 по длине оболочки в её докритическом (штриховая линия) и закритическом (сплошная линия) состояниях.
Анализ результатов численного решения линеаризованной задачи. В соответствии с изложенными выше алгоритмами, реализованными в среде Matlab, решая задачи (23), (60) при различных значениях нагрузки T+, мы находим такое её значение, когда определяемое путём решения задачи (36) собственное число становится равным значению торцевой нагрузки (то есть Л = T+). В такой постановке задачи были найдены точка бифуркации и соответствующее значение нагрузки Л = Т^ = 463.34 кН/м (рис. 24). Был проведен анализ поведения минимального
собственного значения Л задачи (36) при различных значениях нагрузки Т+, результаты которого приведены на рис. 24 в виде соответствующей зависимости. Видно, что решение линеаризованной задачи о возмущенном состоянии оболочки, определяемое в окрестности решения линейной задачи о докритическом НДС, приводит к собственным значениям Л > Т^ .
Рис. 20. Прогибы срединной поверхности верхнего несущего слоя w(2\ см
Рис. 21. Поперечные касательные напряжения в заполнителе qi, МПа
Рис. 22. Мембранные осевые напряжения в верхнем несущем слое ст® = T-)2/(2h(2)), МПа
С целью иллюстрации на рис. 25-28 изображены графики изменения компонент собственного вектора SX, соответствующих собственному значению Л = 463.34 задачи (36), по длине оболочки.
Видно, что в условиях торцевого нагружения подкрепляющего стержня реализовывается ФПУ с выпучиванием оболочки в зоне сопряжения подкрепляющего стержня с оболочкой. На рис. 25 приведены собственные функции прогибов несу-
Рис. 23. Мембранные окружные напряжения в верхнем несущем слое /(2h(2)),
МПа
Рис. 24. Зависимость собственного числа Л от значения нагрузки T+
х10-3
Рис. 25. Приращения прогибов точек срединных поверхностей внешних слоев при переходе в возмущенное состояние
Рис. 26. Приращения касательных напряжений в заполнителе при переходе в возмущенное состояние
щих слоёв, на котором видно, что в граничной области происходит волнообразование прогибов, причём ввиду действия на оболочку сжимающей нагрузки с эксцентриситетом наиболее интенсивное волнообразование наблюдается у верхнего несу-
Рис. 27. Приращения осевых мембранных напряжений в несущих слоях при переходе в возмущенное состояние
Рис. 28. Приращения окружных мембранных напряжений в несущих слоях при переходе в возмущенное состояние
щего слоя в области крепления подкрепляющего стержня с оболочкой. На рис. 2628 изображены собственные функции касательных напряжений в заполнителе, осевых и окружных мембранных напряжений в несущих слоях соответственно, которые в результате потери устойчивости также имеют большие градиенты изменений в области сопряжения с подкрепляющим стержнем.
Благодарности. Результаты исследований получены за счёт гранта Российского научного фонда (проект № 16-11-10299, разработка численных методов) и частично гранта РФФИ (проект № 16-08-00316, постановка задач и вывод основных уравнений) в рамках выполнения государственного задания Минобрнауки России № 9.5762.2017/ВУ (проект № 9.1395.2017/ПЧ, проведение численных экспериментов и анализ их результатов).
Литература
1. Александров А.Я., Брюккер Л.Э., Куршин Л.М., Прусаков А.П. Расчет трехслойных панелей. - М.: Оборонгиз, 1960. - 270 с.
2. Reddy J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis. -Boca Raton: CRC Press, 2004. - 831 p.
3. Badriev I.B., Paimushin V.N. Refined Models of Contact Interaction of a Thin Plate with Positioned on Both Sides Deformable Foundations // Lobachevskii J. Math. - 2017. -V. 38, No 5. - P. 779-793. - doi: 10.1134/S1995080217050055.
4. Birman V., Vo N. Wrinkling in sandwich structures with a functionally graded core // J. Appl. Mech.. - 2017. - V. 84, No 2. - Art. 021002, P. 1-8. - doi: 10.1115/1.4034990.
5. Zenkert D. An Introduction to Sandwich Construction. - London: Chameleon Press Ltd., 1995. - 277 p.
6. Crupi V., Epasto G., Guglielmino E. Comparison of aluminium sandwiches for lightweight ship structures: Honeycomb vs. foam // Marine Structures. - 2013. - V. 30. - P. 74-96. -doi: 10.1016/j.marstruc.2012.11.002.
7. Vasiliev V.V., Morozov E.V. Advanced Mechanics of Composite Materials and Structural Elements. - Elsevier, 2013. - 832 p.
8. Sutherland L., Soares G. Impact behavior of typical marine composite laminates // Composites: Part B. - 2006. - V. 37, No 2-3. - P. 89-100. - doi: 10.1016/j.compositesb.2005.09.001.
9. Johnson H.E., Louca L.A., Mouring S., Fallah A.S. Modelling impact damage in marine composite panels // Int. J. Impact Eng. - 2009. - V. 36, No 1. - P. 25-39. - doi: 10.1016/j.ijimpeng.2008.01.013.
10. Васильев В.В., Добряков А.А., Дудченко А. А. Основы проектирования и изготовления конструкций летательных аппаратов из композиционных материалов. - М.: МАИ, 1985. - 218 с.
11. Takeda N., Minakuchi S., Okabe Y. Smart composite sandwich structures for future aerospace application damage detection and suppression: A review // J. Solid Mechanics and Materials Engineering. - 2007. - V. 1. No 1. - P. 3-17. - doi: 10.1299/jmmp.1.3.
12. Minjing L., Zhanjun W. Application of composite honeycomb sandwich structure in aircraft // Sci. Technol. Rev. - 2016. - V. 34, No 8. - P. 21-25. (in Chinese)
13. Seibert H.F. Composite materials for aerospace applications // Bull. Materials Sci. -1999. - V. 22, No 3. - P. 657-664. - doi: 10.1007/BF02749982.
14. Crump D.A., Dulieu-Barton J.M., Savage J. The Manufacturing Procedure for Aerospace Secondary Sandwich Structure Panels // J. Sandwich Structures Materials. - 2010. -V. 12, No 4. - P. 421-447. - doi: 10.1177/1099636209104531.
15. Григолюк Э.И., Чулков П.П. Устойчивость и колебания трехслойных оболочек. -М.: Машиностроение, 1973. - 168 с.
16. Болотин В.В., Новичков Ю.Н. Механика многослойных конструкций. - М.: Машиностроение, 1980. - 375 с.
17. Григолюк Э.И., Коган Ф.А. Современное состояние теории многослойных оболочек // Прикл. механика. - 1972. - Т. 8, № 6. - С. 5-17.
18. Паймушин В.Н. Нелинейная теория среднего изгиба трехслойных оболочек с дефектами в виде участков непроклея // Прикл. механика. - 1987. - № 11. - С. 32-38.
19. Иванов В.А., Паймушин В.Н. Уточненная теория устойчивости трехслойных конструкций (нелинейные уравнения докритического равновесия оболочек с трансверсально-мягким заполнителем) // Изв. вузов. Матем. - 1994. - № 11. - С. 29-42.
20. Paimushin V.N., Bobrov S.N. Refined geometric nonlinear theory of sandwich shells with a transversely soft core of medium thickness for investigation of mixed buckling forms // Mechanics of Composite Materials. - 2000. - V. 36, No 1. - P. 59-66.
21. Noor A.K., Burton W.S., Bert Ch.W. Computational models for sandwich panels and shells // Appl. Mech. Rev. - 1996. - V. 49, No 3. - P. 155-199. - doi: 10.1115/1.3101923.
22. Паймушин В.Н. Теория устойчивости трехслойных пластин и оболочек (этапы развития, современное состояние и направления дальнейших исследований) // Изв. РАН. Механика твердого тела. - 2001. - № 2. - С. 148-162.
23. Badriev I.B., Makarov M.V., Paimushin V.N. Numerical Investigation of physically nonlinear problem of sandwich plate bending // Procedia Eng. - 2016. - V. 150. -P. 1050-1055. - doi: 10.1016/j.proeng.2016.07.213.
24. Badriev I.B., Banderov V.V., Makarov M.V. Mathematical Simulation of the problem of the pre-critical sandwich plate bending in geometrically nonlinear one dimensional formulation // IOP Conf. Ser.: Materials Science and Engineering. - 2017. - V. 208. -Art. 012002, P. 1-7. - doi: 10.1088/1757-899X/208/1/012002.
25. Badriev I.B., Garipova G.Z., Makarov M.V., Paymushin V.N. Numerical solution of the issue about geometrically nonlinear behavior of sandwich plate with transversal soft filler // Res. J. Appl. Sci. - 2015. - V. 10, No 8. - P. 428-435. - doi: 10.3923/rjasci.2015.428.435.
26. Badriev I.B., Makarov M.V., Paimushin V.N. Solvability of physically and geometrically nonlinear problem of the theory of sandwich plates with transversally-soft core // Russ. Math. - 2015. - V. 59, No 10. - P. 57-60. - doi: 10.3103/S1066369X15100072.
27. Вадриев И.Б., Гарипова Г.З., Макаров М.В., Паймушин В.Н., Хабибуллин Р.Ф. О решении физически нелинейных задач о равновесии трехслойных пластин с трансверсально-мягким заполнителем // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. -2015. - Т. 157, кн. 1. - С. 15-24.
28. Бадриев И.Б., Макаров М.В., Паймушин В.Н. Численное исследование физически нелинейной задачи о продольном изгибе трехслойной пластины с трансверсальномягким заполнителем // Вестн. Перм. нац. исслед. политехн. ун-та. Механика. -2017. - № 1. - С. 39-51. - doi: 10.15593/perm.mech/2017.1.03.
29. Paimushin V.N. theory of moderately large deflections of sandwich shells having a transversely soft core and reinforced along their contour // Mechanics of Composite Materials. -2017. - V. 53, No 1. - P. 1-16. - doi: 10.1007/s11029-017-9636-1.
30. Badriev I.B., Makarov M.V., Paimushin V.N. Contact statement of mechanical problems of reinforced on a contour sandwich plates with transversally-soft core // Russ. Math. -2017. - V. 61, No 1. - P. 69-75. - doi: 10.3103/S1066369X1701008X.
31. Паймушин В.Н. К вариационным методам решения пространственных задач сопряжения деформируемых тел // Докл. АН СССР. - 1983. - Т. 273, № 5. - С. 1083-1086.
32. Бадриев И.Б., Макаров М.В., Паймушин В.Н. Продольно-поперечный изгиб по цилиндрической форме трехслойной пластины, подкрепленной в торцевых сечениях абсолютно твердыми телами // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. -2017. - Т. 159, кн. 2. - С. 174-190.
33. Дартов Р.З., Паймушин В.Н. О методе интегрирующих матриц решения краевых задач для обыкновенных уравнений четвёртого порядка // Изв. вузов. Матем. -1996. - № 10. - С. 13-25.
34. Даутов Р.З., Карчевский М.М., Паймушин В.Н. К методу интегрирующих матриц для систем обыкновенных дифференциальных уравнений // Изв. вузов. Матем. -2003. - № 7. - С. 18-26.
35. Вахитов М.Б. Интегрирующие матрицы аппарат численного решения дифференциальных уравнений строительной механики // Изв. вузов. Авиац. техника. - 1966. -№ 3. - С. 50-61.
36. Badriev I.B., Makarov M.V., Paimushin V.N. Longitudinal and transverse bending by a cylindrical shape of the sandwich plate stiffened in the end sections by rigid bodies // IOP Conf. Ser Materials Science and Engineering. - 2016. - V. 158. - Art. 012011, P. 1-9. -doi: 10.1088/1757-899X/158/1/012011.
37. Бадриев И.Б., Макаров М.В., Паймушин В.Н. Геометрически нелинейная задача о продольно-поперечном изгибе трехслойной пластины с трансверсально-мягким заполнителем // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2016. - Т. 158, кн. 4. - С. 453-468.
38. Карчевский М.М. Итерационные схемы для уравнений с монотонными операторами // Изв. вузов. Матем. - 1971. - № 5. - С. 32-37.
39. Badriev I.B., Makarov M.V., Paimushin V.N. Mathematical Simulation of Nonlinear Problem of Three-point Composite Sample Bending Test // Procedia Eng. - 2016. -V. 150. - P. 1056-1062. doi: 10.1016/j.proeng.2016.07.214.
40. Вайнберг М.М. Вариационный метод и метод монотонных операторов. - М: Наука, 1972. - 416 с.
41. Григолюк Э.И., Шалашилин В.И. Проблемы нелинейного деформирования: Метод продолжения решения по параметру в нелинейных задачах механики твердого деформируемого тела. - М.: Наука, 1988. - 231 с.
42. Шалашилин В.И., Кузнецов Е.Б. Метод продолжения решения по параметру и наилучшая параметризация в прикладной математике и механике. - М.: Эдиториал УРСС, 1999. - 222 с.
Поступила в редакцию 11.04.17
Бадриев Ильдар Бурханович, доктор физико-математических наук, профессор кафедры вычислительной математики
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: ildar.badriev1@mail.ru
Макаров Максим Викторович, ассистент кафедры анализа данных и исследования операций
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: makarovmaksim@mail.ru
Паймушин Виталий Николаевич, доктор физико-математических наук, профессор кафедры прочности конструкций; главный научный сотрудник
Казанский национальный исследовательский технический университет имени
А.Н. Туполева
ул. К. Маркса, д. 10, г. Казань, 420111, Россия Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: vpajmushin@mail.ru
Холмогоров Сергей Андреевич, кандидат физико-математических наук, научный сотрудник
Казанский национальный исследовательский технический университет имени
А.Н. Туполева
ул. К. Маркса, д. 10, г. Казань, 420111, Россия E-mail: hkazan@yandex.ru
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA.
SERIYA FIZIKO-MATEMATICHESKIE NAUKI
(Proceedings of Kazan University. Physics and Mathematics Series)
2017, vol. 159, no. 4, pp. 395-428
The Axisymmetric Problems of Geometrically Nonlinear Deformation and Stability of a Sandwich Cylindrical Shell with Contour Reinforcing Beams
I.B. Badrieva* , M.V. Makarova** , V.N. Paimushina,b*** , S.A. Kholmogorovb****
aKazan Federal University, Kazan, 420008 Russia b Tupolev Kazan National Research Technical University, Kazan, 420111 Russia E-mail: *ildar.badriev1@mail.ru, **makarovmaksim@mail.ru, *** vpajmushin@mail.ru,
**** hkazan@yandex.ru
Received April 11, 2017 Abstract
A numerical investigation of the problem of geometrically nonlinear axisymmetric deformation of a sandwich cylindrical shell with a transversally soft core reinforced in the end sections by elastic rods has been carried out. To describe the process of deformation, we have used the previously derived equations of the refined geometrically nonlinear theory that allow to both study the subcritical behavior of the shell and to reveal all possible buckling forms of the carrier layers. These equations are based on the introduction, as unknown variables, of the contact forces of the interaction of the outer layers with the core, as well as of the outer layers and the filler with the reinforcing bodies at all points on the surfaces of their conjugation. Numerical methods for solving the formulated problems have been developed. They are based on the preliminary reduction of the original problems to a system of integro-algebraic equations, for the solving of which the finite sum method is used. A method has been proposed for investigating the subcritical and supercritical geometrically nonlinear behavior of a shell with its end compression through contour reinforcing rods, according to which unstable equilibrium positions are determined by the method of continuation of the solution with respect to the parameter when the external forces are selected as a parameter. A method has been proposed for finding the critical load (the bifurcation point) at which the shell buckling occurs. This method is based on the linearization of the initial geometrically nonlinear problem in the neighborhood of its nonlinear solution, followed by the formulation of the eigenvalue problem with a nonlinear presence of the parameter. The results of the numerical experiments have been discussed. The results of the experiments have been analyzed.
Keywords: sandwich cylindrical shell, transversally soft filler, contour reinforcing beam, geometric nonlinearity, contact stresses, axial compression, axisymmetric deformation, finite sum method, subcritical and supercritical behavior, bifurcation point, linearized problem, buckling forms
Acknowledgments. The study was supported by the Russian Science Foundation (project no. 16-11-10299) and in part by the Russian Foundation for Basic Research (project no. 16-08-00316) in the framework of the state task of the Ministry of Education of the Russian Federation (task no. 9.5762.2017/VU, project no. 9.1395.2017/PCh).
Figure Captions
Fig. 1. a) The scheme of connection of a sandwich shell with a frame; b) the scheme of the reinforcing beam.
Fig. 2. Sandwich plate shell with contour reinforcing beams.
Fig. 3. Deflections of the middle surfaces of the carrier layers w(k), cm.
(k)
Fig. 4. Axial displacements of the middle surfaces of the carrier layers щ , cm.
Fig. 5. Axial membrane normal stresses in the carrier layers u(1) = /(2h(k)), MPa.
Fig. 6. Bending moments in the carrier layers МЦ , N
Fig. 7. Transverse tangential stresses in the core q1, MPa.
Fig. 8. Circumferential membrane stresses in the carrier layers = T2(k)/(2h(k)), MPa.
Fig. 9. Generalized shear forces of the carrier layers N(k), kN/m.
Fig. 10. Shear forces of the carrier layers Q(k) , kN/m.
(k)
Fig. 11. The deformed state of the middle surface of the upper carrier layer at E = E2 =
= 20 • 103 MPa.
(k)
Fig. 12. The deformed state of the middle surface of the upper carrier layer at E = E2 =
= 80 • 103 MPa.
Fig. 13. The dependence of the running load T((+) on the maximum modulus of the deflection w(2) of the upper layer.
Fig. 14. Deflections of the middle surfaces of the carrier layers w(k), cm.
Fig. 15. Transverse tangential stresses in the core q1, MPa.
Fig. 16. Axial membrane normal stresses in the carrier layers o(1) = T(k)/(2h(k^), MPa.
Fig. 17. Circumferential membrane stresses in the carrier layers = Т2(^/(2h(k)), MPa.
Fig. 18. The dependence of the running load T(+^ on the maximum modulus of the deflection w(2) of the upper layer.
Fig. 19. The deformed state of the middle surface of the upper carrier layer at the bifurcation point
Fig. 20. Deflections of the middle surface of the upper carrier layer w(2) , cm.
Fig. 21. Transverse tangential stresses in the core q1, MPa.
Fig. 22. Membrane axial stresses in the upper carrier layer a(2) = T(2/(2h(2)), MPa.
Fig. 23. Circumferential membrane stresses in the upper carrier layer a22 = T^/(2h(2)), MPa.
Fig. 24. Dependence of the eigenvalue A on the load value T+.
Fig. 25. Increments in the deflections of the points of the middle surfaces of the outer layers at the transition to the disturbed state.
Fig. 26. Increments of the tangential stresses in the core at the transition to the disturbed state.
Fig. 27. Increments in axial membrane stresses in the carrier layers at the transition to the disturbed state.
Fig. 28. Increments of circumferential membrane stresses in the carrier layers at the transition to the disturbed state.
References
1. Aleksandrov A.Ya., Bryukker L.E., Kurshin L.M., Prusakov A.P. Calculation of Sandwich Panels. Moscow, Oborongiz, 1960. 270 p. (In Russian)
2. Reddy J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis. Boca Raton, CRC Press, 2004. 831 p.
3. Badriev I.B., Paimushin V.N. Refined models of contact interaction of a thin plate with positioned on both sides deformable foundations. Lobachevskii J. Math., 2017, vol. 38, no. 5, pp. 779-793. doi: 10.1134/S1995080217050055.
4. Birman V., Vo N. Wrinkling in sandwich structures with a functionally graded core. J. Appl. Mech., 2017, vol. 84, no. 2, art. 021002, pp. 1-8. doi:10.1115/1.4034990.
5. Zenkert D. An Introduction to Sandwich Construction. London, Chameleon Press Ltd., 1995. 277 p.
6. Crupi V., Epasto G., Guglielmino E. Comparison of aluminium sandwiches for lightweight ship structures: Honeycomb Vs. foam. Mar. Struct., 2013, vol. 30, pp. 74-96. doi: 10.1016/j.marstruc.2012.11.002.
7. Vasiliev V.V., Morozov E.V. Advanced Mechanics of Composite Materials and Structural Elements. Elsevier, 2013. 832 p.
8. Sutherland L., Soares G. Impact behavior of typical marine composite laminates. Composites, Part B, 2006, vol. 37, nos. 2-3, pp. 89-100. doi: 10.1016/j.compositesb.2005.09.001.
9. Johnson H.E., Louca L.A., Mouring S., Fallah A.S. Modelling impact damage in marine composite panels. Int. J. Impact Eng., 2009, vol. 36, no. 1, pp. 25-39. doi: 10.1016/j.ijimpeng.2008.01.013.
10. Vasil’ev V.V., Dobryakov A.A., Dudchenko A.A. Fundamentals of the Planning and Production of Aircraft Structures Made of Composite Materials. Moscow, MAI, 1985. 218 p. (In Russian)
11. Takeda N., Minakuchi S., Okabe Y. Smart composite sandwich structures for future aerospace application damage detection and suppression: A review. J. Solid Mech. Mater. Eng., 2007, vol. 1, no. 1, pp. 3-17. doi: 10.1299/jmmp.1.3.
12. Minjing L., Zhanjun W. Application of composite honeycomb sandwich structure in aircraft. Sci. Technol. Rev., 2016, vol. 34, no. 8, pp. 21-25.
13. Seibert H.F. Composite materials for aerospace applications. Bull. Mater. Sci., 1999, vol. 22, no. 3, pp. 657-664. doi: 10.1007/BF02749982.
14. Crump D.A., Dulieu-Barton J.M., Savage J. The manufacturing procedure for aerospace secondary sandwich structure panels. J. Sandwich Struct. Mater., 2010, vol. 12, no. 4, pp. 421-447. doi: 10.1177/1099636209104531.
15. Grigolyuk E.I., Chulkov P.P. Stability and Vibrations of Sandwich Shells. Moscow, Mashinostroenie, 1973. 168 p. (In Russian)
16. Bolotin V.V., Novichkov Y.N. Mechanics of Multilayered Structures. Moscow, Mashinos-troenie, 1980. 375 p. (In Russian)
17. Grigolyuk E.I., Kogan F.A. Present State of the Theory of Multilayered Shells. Prikl. Mekh., 1972, vol. 8. no. 6, pp. 5-17.(In Russian)
18. Paimushin V.N. Nonlinear theory of the central bending of three-layer shells with defects in the form of sections of bonding failure. Sov. Appl. Mech., 1987, vol. 23, no. 11, pp. 10381043. doi: 10.1007/BF00887186.
19. Ivanov V.A., Paimushin V.N. An improved theory of the stability of three-layer structures (non-linear equations for the subcritical equilibrium of shells with a transversely soft filler). Russ. Math., 1994, vol. 38, no. 11, pp. 26-39.
20. Paimushin V.N., Bobrov S.N. Refined geometric nonlinear theory of sandwich shells with a transversely soft core of medium thickness for investigation of mixed buckling forms. Mech. Compos. Mater., 2000, vol. 36, no. 1, pp. 59-66. doi: 10.1007/BF02681778.
21. Noor A. K., Burton W.S., Bert Ch.W. Computational models for sandwich panels and shells. Appl. Mech. Rev., 1996, vol. 49, no. 3, pp. 155-199. doi: 10.1115/1.3101923.
22. Paimushin V.N. The theory of stability of sandwich plates and shells (stages of development, current state and directions of further research). Izv. Ross. Akad. Nauk. Mekh. Tverd. Tela, 2001, no. 2, pp. 148-162. (In Russian)
23. Badriev I.B., Makarov M.V., Paimushin V.N. Numerical investigation of physically nonlinear problem of sandwich plate bending. Procedia Eng., 2016, vol. 150, pp. 1050-1055. doi: 10.1016/j.proeng.2016.07.213.
24. Badriev I.B., Banderov V.V., Makarov M.V. Mathematical simulation of the problem of the pre-critical sandwich plate bending in geometrically nonlinear one dimensional formulation. IOP Conf. Ser.: Mater. Sci. Eng., 2017, vol. 208, art. 012002, pp. 1-7. doi: 10.1088/1757-899X/208/1/012002.
25. Badriev I.B., Garipova G.Z., Makarov M.V., Paymushin V.N. Numerical solution of the issue about geometrically nonlinear behavior of sandwich plate with transversal soft filler. Res. J. Appl. Sci., 2015, vol. 10, no. 8, pp. 428-435. doi: 10.3923/rjasci.2015.428.435.
26. Badriev I.B., Makarov M.V., Paimushin V.N. Solvability of physically and geometrically nonlinear problem of the theory of sandwich plates with transversally-soft core. Russ. Math., 2015, vol. 59, no. 10, pp. 57-60. doi: 10.3103/S1066369X15100072.
27. Badriev I.B., Garipova G.Z., Makarov M.V., Paimushin V.N., Khabibullin R.F. Solving physically nonlinear equilibrium problems for sandwich plates with a transversally soft core. Lobachevskii J. Math., 2015, vol. 36, no. 4, pp. 474-481. doi: 10.1134/S1995080215040216.
28. Badriev I.B., Makarov M.V., Paimushin V.N. Numerical investigation of a physically nonlinear problem of longitudinal bending of sandwich plate with transversal-soft core. PNRPU Mech. Bull., 2017, no. 1, pp. 39-51, doi: 10.15593/perm.mech/2017.1.03.
29. Paimushin V.N. Theory of moderately large deflections of sandwich shells having a transversely soft core and reinforced along their contour. Mech. Compos. Mater., 2017, vol. 53, no. 1, pp. 1-16. doi: 10.1007/s11029-017-9636-1.
30. Badriev I.B., Makarov M.V., Paimushin V.N. Contact statement of mechanical problems of reinforced on a contour sandwich plates with transversally-soft core. Russ. Math., 2017, vol. 61, no. 1, pp. 69-75. doi: 10.3103/S1066369X1701008X.
31. Paimushin V.N. Variational methods for solving non-linear spatial problems of the joining of deformable bodies. Dokl. Akad. Nauk SSSR, 1983, vol. 273, no. 5, pp. 1083-1086. (In Russian)
32. Badriev I.B., Makarov M.V., Paimushin V.N. Longitudinal and transverse bending on the cylindrical shape of a sandwich plate reinforced with absolutely rigid bodies in the front sections. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matemati-cheskie Nauki, 2017, vol. 159, no. 2, pp. 174-190. (In Russian)
33. Dautov R.Z., Paimushin V.N. On the method of integrating matrices for the solution of boundary value problems for fourth-order ordinary equations. Russ. Math., 1996, vol. 40, no. 10, pp. 11-13.
34. Dautov R.Z., Karchevskii M.M., Paimushin V.N. On the method of integrating matrices for systems of ordinary differential equations. Russ. Math., 2003, vol. 47, no. 7, pp. 16-24.
35. Vakhitov M.B. Integrating matrices as a means of numerical solution of differential equations in structural dynamics. Izv. Vyssh. Uchebn. Zaved., Aviats. Tekh., 1966, no. 3, pp. 50-61. (In Russian)
36. Badriev I.B., Makarov M.V., Paimushin V.N. Longitudinal and transverse bending by a cylindrical shape of the sandwich plate stiffened in the end sections by rigid bodies. IOP Conf. Ser.: Mater. Sci. Eng., 2016, vol. 158, art. 012011, pp. 1-9. doi: 10.1088/1757-899X/158/1/012011.
37. Badriev I.B., Makarov M.V., Paimushin V.N. Geometrically nonlinear problem of longitudinal and transverse bending of a sandwich plate with transversally soft core. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2016, vol. 158, no. 4, pp. 453-468. (In Russian)
38. Karchevskii M.M. Iteration schemes for equations with monotone operators. Izv. Vyssh. Uchebn. Zaved., Mat., 1971, no. 5, pp. 32-37. (In Russian)
39. Badriev I.B., Makarov M.V., Paimushin V.N. Mathematical simulation of nonlinear problem of three-point composite sample bending test. Procedia Eng., 2016, vol. 150, pp. 10561062. doi: 10.1016/j.proeng.2016.07.214.
40. Vainberg M.M. Variational Method and Monotone Operators Method. Moscow, Nauka, 1972. 416 p.(In Russian)
41. Grigolyuk E.I., Shalashilin V.I. Problems of Nonlinear Deformation: Parameter Continuation Method in Nonlinear Problems of Solid Mechanics. Moscow, Nauka, 1988. 231 p.(In Russian)
42. Shalashilin V.I., Kuznetsov E.B. Parametric Continuation and Optimal Parametrization in Applied Mathematics and Mechanics. Dordrecht, Boston, London, Kluwer Acad. Publ., 2003. 236 p.
Для цитирования: Бадриев И.Б., Макаров М.В., Паймушин В.Н., Холмого-I ров С.А. Осесимметричные задачи о геометрически нелинейном деформировании ( и устойчивости трехслойной цилиндрической оболочки с контурными подкрепля-\ ющими стержнями // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2017. -
Т. 159, кн. 4. - С. 395-428.
For citation: Badriev I.B., Makarov M.V., Paimushin V.N., Kholmogorov S.A. / The axisymmetric problems of geometrically nonlinear deformation and stability of ( a sandwich cylindrical shell with contour reinforcing beams. Uchenye Zapiski Kazanskogo \ Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2017, vol. 159, no. 4, pp. 395-428.
(In Russian)
