Устойчивость пологих ортотропных оболочек двоякой кривизны при шарнирно-подвижном закреплении контура
Автор: Каменев И.В., Семенов А.А.
Статья в выпуске: 2, 2018 года.
Бесплатный доступ
Оболочечные конструкции часто применяются в разных областях техники, и их исследование важно для многих прикладных задач. Для исключения концентрации напряжений вблизи контура, особенно в угловых точках оболочки, используется шарнирно-подвижное закрепление контура конструкции. В данной работе рассматриваются пологие оболочки двоякой кривизны, квадратные в плане, выполненные из ортотропных материалов и закрепленные по контуру шарнирно-подвижно. Математическая модель основывается на гипотезах теории оболочек Тимошенко - Рейснера, учитывающей поперечные сдвиги, и представлена в виде системы уравнений в смешанной форме. Также учитывается геометрическая нелинейность. Для решения системы дифференциальных уравнений применяется метод Бубнова-Галеркина, что позволяет свести задачу к решению системы нелинейных алгебраических уравнений. Показана сходимость метода при увеличении количества слагаемых аппроксимации. Полученная система является нелинейной и решается методом Ньютона. Разработанный алгоритм реализован в среде аналитических вычислений Maple 2017. Проводится верификация предложенного алгоритма посредством сравнения результатов расчета тестовой задачи с результатом, полученным другими авторами. Совмещение графика зависимости «нагрузка-прогиб» показало хорошую согласованность данных. Проводится анализ устойчивости трех вариантов пологих оболочечных конструкций двоякой кривизны. По каждой из них получены результаты расчета для четырех вариантов ортотропных материалов. На оболочки действует внешняя равномерно-распределенная поперечная нагрузка, закрепление контура - шарнирно-подвижное. Для всех исследованных конструкций приводятся значения критических нагрузок потери устойчивости, значения наибольшего прогиба, соответствующего данным нагрузкам, а также графики зависимости «нагрузка-прогиб». Сделаны выводы о напряженно-деформированном состоянии рассматриваемых оболочек.
Оболочки, математическая модель, компьютерное моделирование, уравнения в смешанной форме, шарнирно-подвижное закрепление, модель тимошенко-рейснера, устойчивость, напряженно-деформированное состояние, метод бубнова-галеркина, метод ньютона
Короткий адрес: https://sciup.org/146281860
IDR: 146281860 | DOI: 10.15593/perm.mech/2018.2.04
Текст научной статьи Устойчивость пологих ортотропных оболочек двоякой кривизны при шарнирно-подвижном закреплении контура
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2018PNRPU MECHANICS BULLETIN
Оболочечные конструкции широко используются в разных областях техники, и их исследование важно для многих прикладных задач: например, подобные конструкции применяются в строительстве и машиностроении [1, 2]. Оболочечные конструкции подвергаются воздействию различных нагрузок [2–14]. Так, в работах [3–6] рассматривается влияние температуры. В работах [7, 8] проводится анализ устойчивости конструкций под воздействием осевого сжатия, причем в работе [7] также проводятся расчеты устойчивости под воздействием точечной нагрузки. В [2, 9–11] рассматривается воздействие статической поперечной равномерно-распределенной нагрузки. В работах [12–14] исследуется процесс деформирования конструкций, находящихся под воздействием динамической нагрузки, причем в [13] также исследуется состояние конструкции под воздействием осевого сжатия.
Конструкции разной геометрии деформируются по-разному, поэтому тип рассматриваемых конструкций также имеет большое значение для исследования их напряженно-деформированного состояния. В большинстве рассмотренных работ [3, 4, 7, 12, 16] исследуются цилиндрические оболочки, в [13, 14] – конические, в [16] исследуются тороидальные, в [2, 17, 19] – пологие двоякой кривизны, в [20] – прямоугольные пластины. В работе [10] приводится общая модель деформирова- ния для оболочек канонической формы (тороидальные, цилиндрические, конические, пологие двоякой кривизны). Исследование устойчивости таких конструкций является актуальной задачей.
Также актуален вопрос облегчения подобных конструкций, для чего применяются современные композиционные материалы, такие как стеклопластики, боро-пластики, органопластики, углепластики и текстолиты, а также полимеры [17, 18, 21, 22]. Их использование позволяет существенно снизить вес конструкции.
Приведенные материалы обладают свойством ортотропии, т.е. его физические свойства различны по взаимно перпендикулярным направлениям. Благодаря этому они превосходят традиционные материалы и сплавы по своим механическим и физическим свойствам. Так, композиционные материалы обладают большой удельной прочностью [19] и стойкостью к воздействиям высоких температур и вибрационным нагрузкам [17].
Применение материалов более сложной структуры приводит к необходимости решения задач уточнения и усовершенствования математической модели их деформирования и выбору устойчивого и точного алгоритма ее исследования [23].
Для исследования оболочечных конструкций применяются различные модели. Так, могут быть использованы математические модели в форме функционала полной потенциальной энергии деформации [2, 23]. Также, что особенно удобно в случае шарнирно- подвижно опертых по контуру конструкций, используются уравнения в смешанной форме [10, 12, 24–28]. Для решения задач устойчивости могут использоваться различные численные методы, такие как метод конечных элементов [2, 20, 29] и метод Ритца [23], а также метод Бубнова–Галеркина [12–15].
Для исключения концентрации напряжений вблизи контура, особенно в угловых точках оболочки, используется шарнирно-подвижное закрепление контура конструкции [30]. Подобное закрепление уменьшает вероятность потери прочности конструкции, что особенно важно, если она изготовлена из ортотропного материала1.
1. Математическая модель деформирования пологой оболочки
(MF (Ф) + k,F2 (Ф)) - £ (F ( Ф) 01 + F3 (Ф)02) -
~( F 2 ( Ф ) 0 2 + F 3 ( Ф ) 0 1 ) + Q + ^ + q = 0; д у d x д у
d dx
—
d
+ 8y
—F ( Ф )-- ^2- f ( ф ) 1 E^h 2V ’ Exh 1V J
EhF (Ф)-
+
^ f ( ф ) I = E 2 h 2 JJ
Будем рассматривать пологие оболочки двоякой кривизны, квадратные в плане. Срединную поверхность оболочки толщиной h примем за координатную. Оси x , y ортогональной системы координат направлены по
= -(x22 -X1X2 + kxX2 + kyX);
df E1 h3 л . . и 4h\. I л A (X1 +^21X2) +Г I 10^12 X12 I Qx = 0;
dx (12(1 - Ц12Ц21) J dу (12
d f E2h3 / , aY d Г h3
2(X?+^i?Xi) +G7X17 - Qv = 0, dy (12 (1 -Ц12Ц21)( 1 )J dx (12 12Л12 J
линиям главных кривизн оболочки. Ось z ортогональна срединной поверхности и направлена в сторону вогну-
тости.
где E , E 2, Pj2, ц21 - модули упругости и коэффициенты
Пуассона материала; G , G , G – модули сдвига в
Параметры Ляме и параметры кривизны пологой оболочки двоякой кривизны соответственно равны плоскостях xOy, x Оz, у Оz соответственно; Ох = dW dx ’
A = 1, B = 1, к = —, к xRyR
SW
O2 =-^^ [10]; d у
W ( x , у ) - перемещение точки средин
Оболочка по контуру закреплена шарнирноподвижно, находится под действием равномерно-распределенной поперечной нагрузки q (рис. 1).
ной поверхности оболочки вдоль оси z ; Xv Х2 — функции изменения кривизн, Хх2 — функция изменения кру- чения,

Рис. 1. Схематичное изображение пологой оболочки Fig. 1. Schematic representation of a shallow shell
d^x д^ у 1 fdTx d^ у I
Х 1 = ; Х 2 = д ; Х 12 = ~I д+ I . (2)
д x д у 2 ^ д у д x J
Здесь функции Т , Т - углы поворота отрезка нормали к срединной поверхности в сечениях x O z , у О z соответственно;
Qx = khG 13 ( Т x -0 , ) ; Qy = khG 23 ( Т у -0 2 ) , к = 5.
В соотношениях (1) F ( Ф ) , F 2 ( Ф ) , F 3 ( Ф ) - обозначения усилий через функцию напряжений в срединной поверхности оболочки Ф ( x , у ) :
В качестве математической модели деформирования оболочечной конструкции воспользуемся уравнениями в смешанной форме, полученными в работе [10] для оболочек общего вида (пологих двоякой кривизны, цилиндрических, конических, тороидальных и др.), а также позволяющими учесть наличие подкрепления ребрами жесткости. После некоторых упрощений для гладких пологих оболочек двоякой кривизны получим
F 1 ( Ф ) =
д 2 Ф _ "i? ;
F 2 ( Ф ) =-гт; F 3 ( ф ) - - . (3)
d x d x д у
Здесь также используется упрощенная форма функций изменения кривизн и кручения, соответствующая модели Кирхгофа–Лява [10]:
„ SO, . d02 . 1 <50, d02I
X1 = "15 ; X2 = "H ; X12 = I ""5" + ""5 I .
d x d x 2 усу d x J
Описанные выше уравнения представляют собой систему уравнений равновесия в смешанной форме для оболочек из ортотропного материала с учетом геометрической нелинейности и поперечных сдвигов. Неизвестными функциями здесь являются W(x,у), Ф(x,у), — x (x, у ), —у (x, у ) .
2. Алгоритм решения уравнений математической модели
Для решения полученных уравнений относительно неизвестных функций W ( x , у ) , Ф ( x , у ) , V Д x , у ) ,
— ( x , у ) применим метод Бубнова-Галеркина. Тогда неизвестные функции будут иметь вид
N
W (x, у ) = £ W (k) X 1( k) Y 1( к);
к = 1
N
Ф ( x , у ) = £ ф ( к ) X 2 ( к ) Y 2 ( к ) ;
к = 1
N
V x ( x , у ) = £ — x ( к ) X 3 ( к ) Y 3 ( к);
к = 1
N
Vу (x,у) = ! —у (к)X4(к)Y4(к).
к = 1
Здесь W ( к ) , Ф ( к ) , — Д к ) , — ( к ) - неизвестные числовые параметры; X 1 ( к ) - X 4 ( к ) , Y 1 ( к ) - Y 4 ( к ) -некоторые аппроксимирующие функции.
Аппроксимирующие функции X 1 ( к ) - X 4 ( к ) ,
Y 1 ( к ) - Y 4 ( к ) выберем исходя из краевых условий при шарнирно-подвижном закреплении контура оболочки:
при x = 0, x = a
аппроксимирующие функции (8), в систему уравнений в смешанной форме (1). Обозначим первое уравнение L 1 ( W , Ф , V x , V у ) , второе - L 2 ( W , Ф , V x , V у ) , третье - L 3 ( W , Ф , V x , V у ) , четвертое - L 4 ( W , Ф , V x , V у ) . Затем строим систему F из 4 N уравнений относительно неизвестных числовых параметров W ( к ) , Ф ( к ) , V x ( к ) , V у ( к ) :
Гг ал ,
JJ L 1 ( W , Ф , V x , V у ) X 1 ( к ) Y 1 ( к ) dxdy
L 00
ab
N
- к = 1
N
F =
И L 2 ( W , Ф , V x , V у ) X 2 ( к ) Y 2 ( к ) dxdy
ab
- к = 1
N
= о
JJ L ( W , Ф , V, , V ) X 3 ( к ) Y 3 ( к ) dxdy
L00 у - 1 = 1
r abb у I N
J J L 4 ( W , Ф , V x , V у ) X 4 ( к ) Y 4 ( к ) dxdy
Система (9) нелинейная в силу нелинейности самой рассматриваемой модели. Для ее решения применим метод Ньютона.
В соответствии с методом Ньютона составим матрицу Якоби F' данной системы
—у = 0, -2- = 0, W = 0, —x = 0, Vу = 0; (6) дx дуд при у = 0, у = b а 2ф а2Ф-—
= 0, = 0, W = 0, —у = 0, —,= 0.(7)
д x 2 д у 2 д у
Тогда в качестве аппроксимирующих можно принять следующие тригонометрические функции:
, _ Г ( 2 к - 1 ) п A , _ Г ( 2 к - 1 ) п A
X 1(к) = sin -----— x ; Y 1(к) = sin
I a J I bJ
, _ Г ( 2 к - 1 ) п A , _ Г ( 2 к - 1 ) п A
X2(к) = sin -----— x ; Y2(к) = sin ab
" . J XJ
, _ Г ( 2 к - 1 ) п A , _ Г ( 2 к - 1 ) п A
X3(к) = cos -----— x ; Y3(к) = sin
I a J I bJ
, _ Г ( 2 к - 1 ) п A , _ Г ( 2 к - 1 ) п A
X4(к) = sin -----— x ; Y4(к) = cos ab
В соответствии с методом Бубнова–Галеркина подставляем разложения искомых функций (5), используя
ГдF1( k )1 д W ( I ) к. ,
Г —> , Ч пА'
д W ( I ) к , "-F3 ( к ) 1 N д W ( I ) t
^F4 ( к ) 1 N l- w ( I ) - к ,.
rfFk ^ ) 1
1дФ ( I ) - к , I
№(к ) 1 N
1-ф < I ) - к ,.
I N
1-ф ( I ) - к ,.
ГдЕ^^ к ) 1 N [-Ф ( I ) - к ,.
r^Fka 1
-v x (I) к , г _ -,/V
- f2 ( к ) l-V x ( I ) - к ,.
- f 3 ( к ) " N -— x ( i )!,. ~aFW 1 N L-v x ( i ) - к.
'-F1 ( к )"
-V у ( I )
- F 2( к )~ |
N |
(10) |
-—, ( I ) |
||
= 0 |
||
- F 3( к )~ |
N |
|
-—, ( I ) |
к , I = 1 |
|
- F 1 ( к )~ -—, ( I ) |
N к , 1 = 1 J |
Наконец, выбирается начальное приближение { W 0 ( I ) , Ф 0 ( I ) , V x 0 ( I ) , V у 0 ( I ) , I = 1.. N } и решается система относительно вектора поправок { A W ( I ) , АФ , ( I ) , AV x 1 ( I ) , AV у 1 ( I ) , I = 1.. N } :
Г A
La w ( I ) - |
||
F ' ( W 0 , Ф 0 , — x 0 , — у 0 ) • |
L Aф 1 ( I ) - |
= - F ( W 0 , Ф 0 , V x 0 , — у 0 ) |
L AV x 1 ( I ) -
IAV у 1 ( I ) -J
Отсюда находятся значения на следующем шаге:
W 1 ( I ) = W 0 ( I ) + A W 1 ( I ) , Ф 1 ( I ) = Ф 0 ( I ) + AФ 1 ( I ) , V x 1 ( I ) = — x ( I ) + AV x ( I ) , V у 1 ( I ) = V у ( I ) + AV у ( I ) ,
I = 1... N (12)
Приведенный итерационный процесс продолжается до тех пор, пока не будет выполнено условие
N
J Z ( A W ( I ) +АФ k ( I ) +AV xk ( I ) + av yk ( I ) )
V I =1 , — < 10 - 3 . (13)
N
J Z ( W k 2 ( I ) +ф 2 ( I ) +V xk ( I ) +V yk ( I ) )
v I =1
Для улучшения скорости сходимости метода рекомендуется использовать в качестве начального приближения на очередном шаге значения, полученные на шаге предыдущем.
Полученные коэффициенты W ( k ) , Ф ( k ) , ¥Д k ) , V ( k ) подставляются в разложение (5), таким образом находятся аппроксимации неизвестных функций при заданной нагрузке q .
При исследовании устойчивости оболочек решается геометрически нелинейная задача и строится график зависимости «нагрузка–прогиб» в некоторой точке оболочки, например в ее центре. Анализируются экстремумы этого графика, и таким образом, находятся верхние и нижние критические нагрузки. При этих нагрузках «хлопком» происходит переход на новое равновесное состояние. По сути дела для нахождения критических нагрузок применяется критерий Ляпунова, когда малому изменению входного параметра (нагрузки) соответствует существенное изменение выходного параметра (прогиба).
-
3. Верификация модели
Прежде чем приступить к исследованию устойчивости оболочечных конструкций, необходимо провести верификацию математической модели для проверки ее адекватности. Сравнение проводится с результатами, полученными в [4] и [7].
В рассмотренных примерах исследовалась цилиндрическая оболочка, выполненная из алюминия, а также шарнирно-подвижно закрепленная по контуру. Цилиндрическую оболочку можно рассматривать как частный случай пологой оболочки двоякой кривизны, но с нулевым коэффициентом кривизны для одной из сторон: kx = 0 или k y = 0, т.е. один из радиусов равняется бесконечности.
Параметры конструкции геометрические:
h = 0,01 м, a = 0,2м, b = 0,2003м, kx = 0 м - 1, k y = 1 м - 1;
физические:
Ev = E2 = 7 • 10 4 МПа, ц12 = ц21 = 0,3,
Gn = 26923,08 МПа, G 13 = G 23 = 76923,08 МПа.
Исследование устойчивости приведенной конструкции проводилось в среде аналитических вычислений Maple с удержанием N = 16 членов разложения в методе Бубнова–Галеркина (5). На рис. 2 приводятся графики зависимости «нагрузка–безразмерный прогиб в центре» для данной цилиндрической панели, полученные в [4] и [7], а также по предложенной в данной работе методике.

Рис. 2. Сравнительный анализ расчета конструкции
Fig. 2. Comparative analysis of the structure calculation
Как видно из представленного графика, все три кривые достаточно близко лежат друг к другу. Таким образом, можно сделать вывод об адекватности разработанной модели и применимости к ней приведенного алгоритма.
-
4. Расчетные параметры
Рассматривались конструкции трех разных размеров (табл. 1) из четырех разных материалов (табл. 2). Таким образом, всего было исследовано 12 конструкций (табл. 3)
Таблица 1
Варианты геометрических характеристик
Variants of geometric characteristics
Table 1
Характеристика |
Вариант геометрии |
||
1 |
2 |
3 |
|
h , м |
0,09 |
0,09 |
0,09 |
а = b , м |
5,4 |
10,8 |
18 |
R 1 = R 2, м |
20,25 |
40,05 |
45,27 |
Таблица 2
Варианты физических характеристик материала
Table 2
Variants of the physical characteristics of the material
Характеристика |
Углепластик M60J/Epoxy [8] |
Графит AS/3501/Epoxy [7] |
E-Glass/Epoxy [7] |
Стеклопластик Т10/UPE22-27/ Epoxy [8] |
Е 1 , МПа |
330000 |
138000 |
60700 |
29400 |
Е 2 , МПа |
59000 |
8960 |
24800 |
17800 |
G 12 = G 13 = G 23 , МПа |
3900 |
7100 |
12000 |
3010 |
µ 12 |
0,320 |
0,300 |
0,230 |
0,123 |
µ 21 |
0,057 |
0,019 |
0,094 |
0,074 |
Показатель анизотропии |
5,59 |
15,4 |
2,44 |
1,65 |
Таблица 3
Варианты рассматриваемых конструкций
Table 3
Variants of the constructions under consideration
Вариант конструкции |
Материал |
Вариант геометрических характеристик |
1 |
Углепластик M60J/Epoxy |
1 |
2 |
Графит AS/3501/Epoxy |
1 |
3 |
E-Glass/Epoxy |
1 |
4 |
Стеклопластик Т10/UPE22-27/Epoxy |
1 |
5 |
Углепластик M60J/Epoxy |
2 |
6 |
Графит AS/3501/Epoxy |
2 |
7 |
E-Glass/Epoxy |
2 |
8 |
Стеклопластик Т10/UPE22-27/Epoxy |
2 |
9 |
Углепластик M60J/Epoxy |
3 |
10 |
Графит AS/3501/Epoxy |
3 |
11 |
E-Glass/Epoxy |
3 |
12 |
Стеклопластик Т10/UPE22-27/Epoxy |
3 |
Для определения оптимального числа членов разложения проведем сравнение получаемых данных при N = 4,9,16,25. Будем рассматривать конструкции 9, 10, 11 и 12, т.е. вариант геометрических характеристик 3, выполненный из всех четырех приведенных материалов. Результаты для варианта конструкции 9 приведены на рис. 3, для варианта 10 – на рис. 4, для варианта 11 – на рис. 5, для варианта 12 – на рис. 6.
Обозначение W i соответствует прогибу в центре конструкции при i = N членах разложения, а W^ -наибольшему прогибу оболочки.
Для оболочки варианта 10 также проведем сравнение значений нормальных напряжений с , с и интенсивности напряжений с . Анализ значений будем проводить в центре конструкции, так как на краях (в том числе в угловых точках) значения всех компонентов напряжений будут равны нулю. Полученные результаты приведены на рис. 7–9 и в табл. 4.

Рис. 4. Сходимость результатов для варианта конструкции 10
Fig. 4. Convergence of results for the 10th variant of the construction

Рис. 3. Сходимость результатов для варианта конструкции 9
Fig. 3. Convergence of results for the 9th variant of the construction

Рис. 5. Сходимость результатов для варианта конструкции 11
Fig. 5. Convergence of results for the 11th variant of the construction

Рис. 6. Сходимость результатов для варианта конструкции 12

Рис. 9. Сходимость интенсивности напряжений с, для варианта конструкции 10
Fig. 9. Convergence of the intensity strain с for the 10th variant of the construction

Рис. 7. Сходимость напряжений с для варианта конструкции 10
Fig. 7. Convergence of ах strains for the 10th variant of the construction
Таблица 4
Значения напряжений в центре оболочки варианта 10 в момент потери устойчивости
Table 4
The strain values corresponding to the critical buckling load in the center of the 10th variant of the construction
N |
σ х , МПа |
σ у , МПа |
σ i , МПа |
4 |
99,38 |
–20,38 |
110,98 |
9 |
–28,55 |
–27,98 |
28,27 |
16 |
–0,48 |
–23,32 |
23,08 |
25 |
–13,01 |
–26,03 |
22,55 |

Рис. 8. Сходимость напряжений с y для варианта конструкции 10
Fig. 8. Convergence of с strains for the 10th variant of the construction
Полученные значения напряжений для оболочки варианта 10 показывают хорошую сходимость для σ у и σ i , для σ х сходимость есть, но выражена слабее.
Как можно заметить из представленных в данном разделе графиков и таблиц, расчеты на устойчивость при N = 16 и N = 25 дают близкий результат вне зависимости от материала, из которого изготовлена конструкция. Таким образом, все дальнейшие расчеты будут проводиться при удержании N = 16 членах разложения в методе Бубнова–Галеркина (5).
-
5. Исследование устойчивости
Проводится расчет устойчивости 12 конструкций, приведенных в табл. 3, с геометрическими характеристиками, указанными в табл. 1, и физическими характеристиками, указанными в табл. 2.
На рис. 10 приведены графики «нагрузка–прогиб» для варианта геометрических параметров 1 для всех четырех материалов. Индекс W k соответствует наибольшему прогибу конструкции из материала k , а индекс W k – прогибу в центре. Для рассматриваемой конструкции наибольший прогиб равен прогибу в центре оболочки: W ^ax = W* .

Рис. 10. График зависимости «нагрузка-прогиб» для варианта геометрических параметров 1
Fig. 10. Graph of the "load-deflection" dependence for the 1st variant of geometric parameters
На рис. 11 приведены графики «нагрузка-прогиб» для варианта геометрических параметров 2 для всех четырех материалов.
На рис. 12 приведены графики «нагрузка-прогиб» для варианта геометрических параметров 3 для всех четырех материалов.

Рис. 12. График зависимости «нагрузка-прогиб» для варианта геометрических параметров 3
Fig. 12. Graph of the "load-deflection" dependence for the 3rd variant of geometric parameters

Рис. 11. График зависимости «нагрузка-прогиб» для варианта геометрических параметров 2
Fig. 11. Graph of the "load-deflection" dependence for the 2nd variant of geometric parameters
Анализируя полученные результаты, можно заметить, что для конструкций с вариантом геометрических параметров 1 прогиб в центре совпадает с максимальным для всех четырех материалов. Это связано с отношением длины сторон оболочки к ее толщине. Для остальных двух вариантов геометрических параметров при той же толщине размеры в плане существенно больше, чем у первого варианта. Отсюда возникает различие между значениями прогиба в центре конструкции и ее максимальным прогибом, особенно вблизи критической нагрузки потери устойчивости. Для вариантов конструкций 5–8 (вариант геометрических параметров 2) это характерно только для стекловолокна E-Glass/Epoxy, а для вариантов конструкций 9–12 (вариант геометрических параметров 3) – уже для всех исследуемых материалов.
Все значения критических нагрузок q , при которых оболочки теряют устойчивость, сведены в табл. 5, а значение наибольшего прогиба, соответствующего этой нагрузке, – в табл. 6.
Таблица 5
Значения критических нагрузок потери устойчивости qkr
Table 5
The values of the critical loads of buckling qkr
Вариант геометрии |
q kr , Мпа |
|||
Углепластик M60J/Epoxy |
Графит AS/3501/Epoxy |
E-Glass/Epoxy |
Стеклопластик Т10/UPE22-27/ Epoxy |
|
1 |
0,575 |
0,205 |
0,265 |
0,105 |
2 |
0,081 |
0,059 |
0,120 |
0,036 |
3 |
0,098 |
0,067 |
0,070 |
0,030 |
Таблица 6
Значения наибольшего прогиба W kr , соответствующего потере устойчивости
The values of the maximum deflection W kr corresponding to the critical buckling load
Table 6
Вариант геометрии |
kr max , м |
|||
Углепластик M60J/Epoxy |
Графит AS/3501/Epoxy |
E-Glass/Epoxy |
Стеклопластик Т10/UPE22-27/ Epoxy |
|
1 |
0,300 |
0,130 |
0,115 |
0,135 |
2 |
0,261 |
0,155 |
0,115 |
0,159 |
3 |
0,328 |
0,143 |
0,091 |
0,109 |

Рис. 13. Критические нагрузки потери устойчивости в зависимости от выбранного материала
Fig. 13. The critical buckling loads, depending on a chosen material

Рис. 14. Значения максимального прогиба, соответствующего потере устойчивости, в зависимости от выбранного материала
Fig. 14. The values of the maximum deflection corresponding to the critical buckling load, depending on a selected material
Для лучшего восприятия данных в табл. 5, 6 те же самые значения приведены на рис. 13 и 14 соответственно. На следующих графиках индекс W i обозначает максимальный прогиб оболочки из материала i , соответствующий потере устойчивости, то есть при нагрузке q = q i .
Как видно из представленных графиков и таблиц, наибольшее значение критической нагрузки потери устойчиво- сти из рассмотренных материалов дают углепластик M60J/Epoxy и стекловолокно E-Glass/Epoxy. При этом оболочки, выполненные из стекловолокна E-Glass/Epoxy, к моменту достижения критической нагрузки потери устойчивости имеют наименьший максимальный прогиб среди всех четырех рассматриваемых материалов, в то время как конструкции из углепластика – наибольший.
Таким образом, можно говорить о том, что наилучшим выбором среди четырех рассмотренных материалов будет E-Glass/Epoxy в силу высокой устойчивости и большой несущей способности выполненных из него конструкций. Однако заметим, что углепластик M60J/Epoxy в данном случае был выбран однонаправленный, в то время как стеклопластик имеет волокна, уложенные в обоих направлениях. Укладка волокон углепластика в обоих направлениях может сделать его наиболее оптимальным материалом.
Заключение
В ходе исследования рассматривались пологие оболочки двоякой кривизны, прямоугольные в плане. Оболочки закреплены по контуру шарнирно-подвижно и выполнены из ортотропного материала. По результатам исследования были получены следующие выводы и результаты:
Список литературы Устойчивость пологих ортотропных оболочек двоякой кривизны при шарнирно-подвижном закреплении контура
- Кривошапко С.Н. О возможностях оболочечных сооружений в современной архитектуре и строительстве//Строительная механика инженерных конструкций и сооружений. -2013. -№1. -С. 51-56.
- Соловей Н.А., Кривенко О.П., Малыгина О.А. Конечно-элементные модели исследования нелинейного деформирования оболочек ступенчато-переменной толщины с отверстиями, каналами и выемками//Инженерно-строительный жур. -2015. -№ 1(53). -С. 56-69. DOI: 10.5862/MCE.53.6
- Спасская М.В., Трещев А.А. Термоупругое деформирование цилиндрической оболочки из анизотропного разносопротивляющегося материала//Вестн. Чуваш. гос. пед. ун-та им. И.Я. Яковлева. Сер. Механика предельного состояния. -2015. -№ 1. -С. 65-74.
- Duc N.D., Tung H.V. Nonlinear response of pressure-loaded functionally graded cylindrical panels with temperature effects//Composite Structures. -2010. -Vol. 92 -P. 1664-1672. DOI: 10.1016/j.compstruct.2009.11.033
- Гумерова Х.С. Влияние деформации поперечного сдвига на устойчивость ортотропной термочувствительной цилиндрической оболочки//Вестн. Казан. технол. ун-та. -2017. -Т. 20, № 7. -С. 91-92.