(Ю. Г. Коноплев и др.) - Динамическая устойчивость упругих пластин и оболочек - учебное пособие-Казанский ун-т (2012)

Скачать как pdf или txt
Скачать как pdf или txt
Вы находитесь на странице: 1из 80

ДИНАМИЧЕСКАЯ УСТОЙЧИВОСТЬ

УПРУГИХ ПЛАСТИН И
ОБОЛОЧЕК

Учебное пособие

Казанский университет
2012
УДК 539.3
Д46

Печатается по решению Редакционно-издательского совета


ФГАОУВПО «Казанский (Приволжский) федеральный университет»

Авторы:
Ю.Г. Коноплев, Л.У. Бахтиева, В.И. Митряйкин,
Ф.Х. Тазюков

Рецензенты:
доктор физ.-мат. наук Ю.П. Артюхин (К(П)ФУ)
доктор физ.-мат. наук В.Н. Паймушин (КНИТУ-КАИ)

Д46 Динамическая устойчивость упругих пластин и оболочек: учеб.


пособие / Л.У. Бахтиева, Ю.Г. Коноплев, В.И. Митряйкин, Ф.Х. Тазюков. –
Казань: Казан. ун-т, 2012. – 80 с.

ISBN 978-5-905787-30-0

Предназначено для студентов и аспирантов, специализирую-


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

ISBN 978-5-905787-30-0 © Коллектив авторов, 2012.


© Казанский университет, 2012.
ВВЕДЕНИЕ

Тонкостенные оболочки и пластины являются важными


составными частями многих конструкций: это корпуса
всевозможных машин и судов, фюзеляжи и крылья летательных
аппаратов, различные емкости, резервуары и другие инженерные
сооружения.
Одной из основных причин разрушения оболочечных
конструкций является потеря устойчивости исходной формы
равновесия. Разрушение при этом происходит практически
мгновенно, что может привести к тяжелым последствиям, поэтому
расчет конструкций на устойчивость очень ответственен. Вместе
с тем он и очень сложен, особенно когда внешние нагрузки
быстро изменяются во времени [6, 11, 20].
Получение достоверных данных о критических нагрузках
требует исследования нелинейных уравнений, для их решения
применяются в основном численные методы, на основе которых
разрабатываются программные комплексы [5, 9, 14, 22].
Однако в расчётной практике придается большое значение
разработке и совершенствованию простых аналитических методов
решения конкретных задач, наглядно отражающих влияние
отдельных геометрических и физических параметров на
устойчивость конструкции. Именно такой эффективный подход к
задачам устойчивости обеспечивает теоретико-
экспериментальный метод, который позволяет получить решение
в аналитическом виде, а затем скорректировать полученные
формулы с помощью экспериментальных данных [8, 10, 17, 19].
Цель учебного пособия – изложить основные положения
теоретико-экспериментального метода применительно к решению
задач устойчивости цилиндрических и конических оболочек,
прямоугольных и круглых пластин, а также пологих панелей под
действием динамических нагрузок; наглядно продемонстрировать
достоинства этого метода при решении целого ряда практически
важных задач.

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

4
ПРИНЯТЫЕ ОБОЗНАЧЕНИЯ

х, у, z – декартовы координаты,
t – время,
h – толщина оболочки (пластины),
Е – модуль упругости,
ν – коэффициент Пуассона,
ρ – плотность материала оболочки,
E
V = – скорость звука в материале, ,
ρ
Е h3
D=
(
12 1 − ν 2 ) – изгибная жёсткость,

u, v, w – перемещения точек срединной поверхности,


Т1 , Т2 , S – усилия в срединной поверхности,
F – функция усилий,
φ, r – криволинейные координаты,
γ – угол конусности оболочки,
φ1 = φ sin γ,
R0 , R1 – радиусы верхнего и нижнего оснований конуса,
Н – высота конуса,
R – радиус цилиндрической оболочки,
l – длина цилиндрической оболочки,
d= a/b; а, b – стороны прямоугольной пластинки (рисунок 1),
m – число полуволн вдоль образующей,
n – число волн по окружности,
ω – частота собственных колебаний,
∂4 ∂4 ∂4
∇4 = + 2 +
∂х 4 ∂х 2 ∂у 2 ∂у 4 ,

∂ 4 2 ∂3 1 ∂ 2 1 ∂ 2 ∂ 4 2 ∂3 4 ∂2 1 ∂4
∇∇ = 4+
2 2
− + + − + +
∂r r ∂r 3 r 2 ∂r 2 r 3 ∂r r 2 ∂φ12∂r 2 r 3 ∂φ12∂r r 4 ∂φ12 r 4 ∂φ14 .

5
1. ПОСТАНОВКА ЗАДАЧИ

Рассматриваются круговые цилиндрические оболочки под


действием динамического внешнего давления при различных
способах закрепления торцов, прямоугольные и круглые пластины
под действием динамического сжатия или динамического сдвига,
конические оболочки под действием динамического внешнего
давления (закон изменения давления – линейный, квадратичный и
т. п.), пологие цилиндрические панели под действием поперечной
нагрузки. Для расчета таких конструкций на устойчивость
предлагается:
1. Использовать линеаризованные уравнения динамики
тонких оболочек, полученные из общих нелинейных уравнений
путем варьирования (с учетом только поперечных сил инерции).
2. Проинтегрировать уравнения по координатам с помощью
процедуры Бубнова, в результате получить одно обыкновенное
дифференциальное уравнение второго порядка.
3. Проинтегрировать уравнение по времени, при этом:
a) задать начальные условия и сформулировать критерий
динамической устойчивости,
b) разбить временной интервал на два участка (в
зависимости от направления сил инерции),
c) ввести в рассмотрение полный прогиб оболочки, что
позволит получить нетривиальное решение для нулевых
начальных условий,
d) выбрать аппроксимирующую функцию прогиба,
удовлетворяющую начальным условиям и критерию
динамической устойчивости,
e) применить процедуру Бубнова,

6
f) полученные зависимости дополнить требованием, чтобы в
предельном случае (при медленном нагружении) они совпадали с
соответствующими формулами для статических задач.
4. Рассчитать критические нагрузки, волновые числа,
коэффициенты динамичности, критические импульсы.
Результаты получить в виде формул, справедливых при
произвольных параметрах оболочек и пластин.
5. Уточнить экспериментальным путем коэффициенты,
входящие в полученные структурные зависимости.

2. МЕТОД РЕШЕНИЯ ЗАДАЧИ

Метод решения описан для задачи устойчивости круговой


цилиндрической оболочки под действием линейно
изменяющегося во времени внешнего поперечного давления.

2.1. Вывод разрешающего уравнения

Исходные нелинейные уравнения динамики цилиндрической


оболочки приведены в [6]

∂2w ∂2w ∂2w 1 ∂2F ∂2w


D ∇4w = T + T + 2 S + + P − ρ h , (1)
∂x 2 ∂y 2 ∂x∂y R ∂x 2 ∂t 2
1 2

⎡⎛ ∂ 2 w ⎞ 2 ∂ 2 w ∂ 2 w ⎤ Еh ∂ 2 w
∇ 4 F = Еh ⎢⎜⎜ ⎟⎟ − 2 ⎥−
⎢⎣⎝ ∂x ∂у ⎠ ∂x ∂у 2 ⎥ R ∂x 2 ,

где Р = Р0 t – динамическое внешнее давление, Р0 – скорость


нагружения.

7
После применения операции варьирования и допущения о
безмоментности первоначального напряжения система (1)
приобретает вид

1 ∂2 w
D ∇4 w − F xx − T10 w xх − T20 w уу − 2S 0 w ху + ρ h 2 = 0 , (2)
R ∂t
Eh
∇4 F = w xx
R

где w , F - вариации функций прогиба и усилий, Т10 , Т20 , S –


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

∂T10 ∂S 0 ∂S 0 ∂T20 T20


+ =0, + = 0, = − P. (3)
∂x ∂y ∂x ∂y R

Для рассматриваемой задачи решение системы (3) имеет вид

Т10 = S0 = 0, Т20 = – РR = – Р0 R t. (4)

К системе (2) применяется процедура метода Бубнова. При


этом для случая шарнирного опирания аппроксимирующие
функции выбираются в виде

F = f (t ) sin w = ξ (t ) sin
m1 x ny m1 x ny
sin , sin , (5)
R R R R

mπR
m1 =
l

После интегрирования система (2) сводится к одному


дифференциальному уравнению

8
d 2ξ
+ λ (t ) ξ = 0 ,
dt 2

или после введения безразмерного параметра времени

τ = t t (tк – критическое время, соответствующее моменту


k

хлопка), к уравнению вида

d 2ξ
+ λ (τ ) ξ = 0 , 0 ≤ τ ≤ 1, (6)
dτ 2

tк2 ⎡ ⎛ 4 β 2 m14 ⎞⎤ Eh R 2
λ (τ ) = ⎢ 0− P R 3 2
n t τ + D ⎜
⎜ n + ⎟ ⎥ , β 2
=
ρ h R4 n 4 ⎟⎠⎥⎦
к
⎢⎣ ⎝ D .

Замечание. При выводе формулы (6) сделано допущение:


m1 << n, что справедливо при достаточно высоких скоростях
нагружения.

2.2. Начальные условия и критерий динамической


устойчивости

Поскольку функция ξ(τ) является вариацией полного прогиба


оболочки, начальные условия для неё имеют вид

ξ (0) = 0, ξʹ (0) = 0. (7)

В задачах динамики необходимо также сформулировать


условие на правом краю интервала времени, т. е. критерий
динамической устойчивости. Из работы А.В. Саченкова [18]
следует, что это условие должно быть выбрано в виде

9
ξʹ (1) = 0, (8)

т. е. критическое время должно соответствовать обращению в


нуль скорости характерных точек оболочки.

2.3. Преобразование разрешающего уравнения

Задача Коши (6) – (7) имеет лишь тривиальное решение.


Поэтому введем в рассмотрение амплитуду полного прогиба
оболочки

W = W0 + ξ , (9)

где W0 – амплитуда докритического (безмоментного) прогиба.


Подставляя в (6) ξ = W – W0, приводим уравнение к виду

W0
Z ′′ + λ (τ ) Z = λ (τ ) , Z = W/h.
h

Из третьего уравнения системы (3) следует

EhW0 T20 P0 R 2t
= = − P0 t , откуда W0 = − , (10)
R2 R Eh

тогда уравнение для Z можно записать в виде

Z + λ(τ) Z = – λ(τ) P 0 τ, P 0 = P0 R2 tк /Eh2. (11)

Из условий (7) и (8) с учетом (9) следуют

Z (0) = Zʹ (0) = 0 – начальные условия, (12)

Zʹ (1) = 0 – критерий устойчивости. (13)

10
Задача, таким образом, сводится к решению уравнения (11) с
начальными условиями (12) и дополнительным условием (13). Для
решения используется процедура метода Бубнова.

2.4. Выбор аппроксимирующей функции

Аппроксимирующая функция прогиба выбирается в


соответствии с начальными условиями (12) в виде

(
Z = A τ 2 + Bτ α , )
где А – амплитуда прогиба, α – константа, определяемая в
процессе решения задачи ( α ≠ 2 ). Из условия (13) следует, что
B = – 2/α , тогда

⎛ 2 ⎞
Z = A ⎜τ 2 − τ α ⎟ . (14)
⎝ α ⎠

Абсцисса точки перегиба для кривой (14) определяется


равенством

1
⎛ 1 ⎞ α −2
τ1 = ⎜ ⎟ .
⎝ α −1 ⎠
(15)

2.5. Решение уравнения

Уравнение (11) необходимо интегрировать по двум участкам:


0≤ τ ≤ τ1 и τ1 ≤ τ ≤ 1, так как при переходе через точку τ1
меняется тип уравнения (коэффициент λ(τ) меняет знак).
В результате интегрирования на первом участке (0, τ1)
получается уравнение

11
( ) (
А с1 + с2 В − с3 А + Р 0 с4 В − с5 А = 0 , ) (16)

t к2 t2 β 2 m1
4

где А= ⋅ P0 R 3 n 2 t к , В = к 4 D ( n4 + )
ρ hR 4
ρ hR n4 ,

2τ 13 (2α 4 − 7α 3 + 7α 2 − 8α + 12)
с1 = ,
3α (α 2 − 1) (2α − 1)

τ 15 (2α 6 + 3α 5 − 9α 4 − 39α 3 + 23α 2 + 40α + 60)


с2 =
5α 2 (α − 1) 2 (α + 3) (2α + 1) ,

τ 16 (α 6 + 3α 5 − 5α 4 − 27α 3 + 4α 2 + 36α + 48)


с3 = ,
6α 2 (α − 1) 2 (α + 4) (α + 1)

⎡1 2 ⎤ 5 ⎡1 2 ⎤
с4 = τ 14 ⎢ − ⎥, с5 = τ 1 ⎢ − ⎥.
⎣ 4 α (α − 1) (α + 2) ⎦ ⎣ 5 α (α − 1) (α + 3) ⎦

В уравнении (16) при tк → ∞ и Р0 → 0 должен


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

А с2 + Р 0 с4 = А с3 + Р 0 с5.

Из этого условия находим амплитуду

А= Р0 (с5 – с4) / (с2 – с3). (17)

Уравнение (16) с учётом (17) преобразуется к виду

12
D 2 β 2 m14 b1ρ h R
Р0 tк = (n + 6 )+ 2 2 , (18)
R3 n n tк

где b1 = с1 (с4 – с5) / (с3 с4 – с2 с5) – величина, зависящая от


параметра α. Численные значения этого коэффициента для
различных значений α приведены в таблице 1.

Таблица 1
Зависимости параметров τ1, b1, b2 и b3 от коэффициента α

α τ1 b1 b2 b3
1,5 0,250 1825,82 5,31 535,25
1,6 0,279 1306,06 1,88 88,30
1,7 0,305 998,68 1,60 53,67
1,8 0,328 800,58 1,50 41,05
1,9 0,349 664,59 1,44 34,60
3,0 0,500 232,73 1,26 21,95
4,0 0,577 154,55 1,21 21,98
5,0 0,630 120,19 1,18 23,20
10,0 0,760 67,48 1,11 31,51
20,0 0,849 45,19 1,07 47,18
30,0 0,887 37,74 1,05 61,12
40,0 0,908 33,89 1,04 72,75
50,0 0,922 31,50 1,04 86,17
60,0 0,932 29,87 1,03 97,83
70,0 0,940 28,67 1,03 109,10

Результат интегрирования по второму участку (τ1, 1) имеет


вид
13
А (а1 + а2 В - а3 А )+ Р0 (а4 В - а5 А )= 0, (19)

где

− 8α 3 + 32α 2 − 8α − 12 2α 4 + 7α 3 − 37α 2 + 60
а1 = − с1 , а 2 = − с2 ,
3 (α + 1)(2α − 1) 5α 2 (α + 3)(2α + 1)

α 4 + 5α 3 − 20α 2 − 12α + 48 α 2 + 2α − В
а3 = − с а = − с4 ,
6α 2 (α + 1)(α + 4) , 4α (α + 2)
3 4

α 2 + 3α − 10
а5 = − с5 .
5α (α + 3)

С учетом формулы (17) уравнение (19) преобразуется к виду

b2 D 2 β 2 m14 b3 ρ h R
Р0 tк = (n + 6 )+ 2 2 , (20)
R3 n n tк

где коэффициенты b2 и b3 зависят от α (таблица 1).


Эксперименты показывают, что волновое число m для
оболочек средней длины можно принять равным единице, а число
n существенно зависит от скорости нагружения. Для его
определения к уравнениям (18) и (20) необходимо добавить
условие минимума нагрузки по волновому числу n, которое на
основании (18) можно записать в виде

2 D 2 β 2 m14
Р0 t к = 3 ( n − 6 ) . (21)
R n

14
Задача, таким образом, сводится к решению системы
алгебраических уравнений (18), (20), (21) относительно
неизвестных tк , n и α.
После преобразований получается уравнение


(1 – а )2 (1–3 а ) = P0 а , (22)

β 2 m14 ↔
2b1 + b3 − b1b2 ~ b ρ h R10 P 2
где а= = , P0 = 1 3 2 40 .
b1b2 + 2b1 + 3b2 n8 4 D β m1 (23)

Величина α входит в уравнение (22) сложным образом.


Поэтому решить его можно лишь численно.
Найденное значение α подставляется в (23) для определения
волнового числа n, а затем в (21) для определения критического
времени tк (или критической нагрузки Рк=P0tk ). В таблице 2
приведены примеры расчета оболочки с параметрами: R = 0,09 м,
ρ = 3,04*103 нсек2/м4, h = 8*10-4 м, l = 0,2 м, Е = 7,56*1010 ПА для
шарнирного опирания и для заделки, результаты сравниваются с
экспериментальными и теоретическими данными, имеющимися в
литературе.

2.6. Получение расчётных формул

Расчёты показывают (таблица 2), что при значениях параметра


Р 0 ≤ 10 (что соответствует небольшим скоростям нагружения Р0)
решение динамической задачи мало отличается от статического
решения. При Р 0 >10 (большие скорости нагружения) становится
возможным упростить уравнения (18), (20), (21) и получить
решение в виде аналитических формул. Из таблицы 2 видно, что с

15
Таблица 2
Расчёты для цилиндрических оболочек

Граничные
условия Шарнирное
Заделка Эксп. [12] Теор. рез. раб. [6]
Скорость опирание
нагружения
Р0 α n P0 tк·10 -5Па α n P0 tк·10 -5Па n P0 tк·10 -5Па n P0 tк·10 -5Па

0 75 6 2,4 510 7 3,6 6 2,5 6 2,4

200 14 6 2,4 312 7 3,6 – – – –


500 2,1 6 2,4 73 7 3,6 – – – –
600 1,6 6 3,3 31 7 4,6 – – – –
1000 1,59 7 4,9 14 9 7,7 – – – –
2000 1,58 8 6,7 2 9 8,0 7 6,5 7 5,2

4700 1,58 9 8,9 1,59 10 10 8 8,1 9 7,4

6500 1,58 10 10,6 1,58 11 11 9 10 10 8,8

16
ростом скорости нагружения параметр α уменьшается, его
значение стабилизируется и становится ≈ 1,6. При этом b1(α)
≈ 1300 (таблица 1). Для решения задач, таким образом, нет
необходимости находить α из уравнения (22), а можно
воспользоваться только уравнениями (18) и (21), где коэффициент
b1 известен. Пренебрегая в указанных уравнениях величиной
β 2 m14
по сравнению с единицей (что при высоких скоростях
n8
нагружения оправдано), получаем для определения величин Р0 tк
и n следующие зависимости

D n 2 b1 ρ hR 2 Dn 2
Р 0 tк = + 2 2 , Р0 t к = , (24)
R3 n tк R3

из которых находим

b1 ρ hR10 P02 R10 P02


n8 = = 432 (1 − ν 2 3
) b
V 2 E 2 h8 ,
1
4D3
1/ 4
2 D ⎛ b ρ hR10 P02 ⎞
1/ 2
c h ⎛ EP0 R ⎞
P0 t к = 3 ⎜⎜ 1 ⎟⎟ = ⎜ ⎟ , 25)
R ⎝ 4D 3 ⎠ 4
1 −ν 2 R⎝ V ⎠

где значение константы с = 4,5 при b1 = 1300. Выражение (25)


совпадает с аналогичной зависимостью, полученной в работе [12].
Формула для коэффициента динамичности Kд на основании
(25) и известной формулы для статической критической нагрузки
[6] будет

17
1
Р t ⎛ Р0 R Rl 2 ⎞ 2 4 Dβ 1/ 2 m
К ∂ = 0 к к ≈ 1,2 с 1 − ν 2 ⎜⎜ ⋅ 3 ⎟⎟ , Р кст = 3 / 4 3 1 . (26)
Рст ⎝ EV h ⎠ 3 R

Можно получить также формулу для критического импульса


Iк, связанного с нагрузкой соотношением

tк tк
P0 t к2
I к = ∫ Р(t )dt = ∫ Р0 t dt =
0 0
2 .

С учетом (25) формула для критического импульса имеет вид

c2 h2 E
Iк = ⋅
2 4 1 −ν 2 RV . (27)

В работе [11] к решению системы уравнений (1) предложен


другой поход, основанный на их предварительном
интегрировании по времени. При этом результаты расчетов для
оболочки с теми же параметрами мало отличаются от данных
таблицы 2.
Однако следует отметить, что входящая в уравнение (24)
величина b1, а следовательно, и константа c в формулах (25)–(27)
зависят от метода решения уравнения (11), от выбора
аппроксимирующей функции (14) и от критерия динамической
устойчивости.
Поэтому значение константы рекомендуется уточнить
экспериментальным путём. Для уточнения константны можно
также воспользоваться экспериментальными данными,
имеющимися в литературе [4, 12].

18
3. РЕШЕНИЕ НЕКОТОРЫХ ЗАДАЧ

3.1. Общий алгоритм решения

Для расчёта оболочки или пластины на динамическую


устойчивость необходимо:
a) вывести линеаризованные уравнения динамики оболочки
(пластины),
b) проинтегрировать уравнения методом Бубнова по
пространственным координатам и получить уравнение вида (6),
где коэффициент λ (t) зависит от конкретной задачи,
c) перейти от уравнения (6) к уравнению (11),
d) получить уравнение типа (18), не вычисляя конкретных
значений коэффициента b1,
e) записать условие минимума критической нагрузки по
волновым числам типа уравнения (21),
f) решить совместно уравнения (18), (21) и получить
расчётные формулы для основных динамических характеристик,
g) уточнить экспериментальным путём константы, входящие
в формулы.

3.2. Цилиндрическая оболочка при внешнем


давлении (другие случаи закрепления торцов)

Линеаризованные уравнения динамики имеют вид (2).


Аппроксимирующие функции выбираются согласно выводам
работы [7] в виде

кm1 x ny кm1 x ny
F = f (t ) sin sin , w = ξ (t ) sin sin ,
R R R R

19
где к – коэффициент, характеризующий способ закрепления
торцов оболочки.

Таблица 3
Результаты расчетов для различных видов закрепления

Граничные Граничные Критическое давление (кПа)


условия при условия при к кПа кПа
2 ⋅ 105
кПа
х=0 х=l Р0=0 104 105
с с с

Шарнирное опирание
1 240 240 490 670

357 358 588 697


Жесткое защемление
1,5 γ=1,5 1,49 1,2 1,04

1,25 298 299 515 677


w=wхх=Т1=0 w=wхх=и=0
γ=1,5 1,25 1,05 1,01

0,5 119 149 485 670


Свободный край
γ=0,5 0,62 0,99 1,00

wxxx=wxx=T1 0,6 142 163 470 670


w=wх=и=0
=S=0 γ=0,6 0,68 0,96 1,00

179 191 475 663


wххх=wx=и=S=0 0,75
γ=0,7 0,8 0,97 0,991

В таблице 3 приведены значения к для различных видов


закрепления. Для нахождения динамических характеристик
достаточно во всех уравнениях пунктов 2.5 и 2.6 заменить m1 на
km1. В частности, формула (25) показывает, что значение

20
динамической критической нагрузки при высоких скоростях
нагружения не зависит от характера закрепления торцов оболочки,
т. к. величина m1 в эту формулу не входит. Коэффициент γ в
таблице 3 равен отношению критического давления при заданном
виде закрепления к критическому давлению при шарнирном
опирании.

3.3. Влияние предварительного статического


нагружения

В случае, когда шарнирно опертая цилиндрическая оболочка


предварительно нагружена статическим внешним давлением q и
статическим осевым сжатием T, характеристическое уравнение
(18) принимает вид

⎡ β 2 m14 ⎤ b1ρ h R 4
( )
2
( Р0 tк + q) n 2 R 3 + Tm12 R 2 = D ⎢ m12 + n 2 + 2⎥
+ .
⎢⎣ (
m12 + n 2 )⎥⎦ tк2

Результаты вычислений критического внешнего давления для


предварительно нагруженной оболочки, а также
экспериментальные данные получены в [8] и отображены на
рисунке 1, где Pk – динамическая критическая нагрузка в случае

ненагруженной оболочки, P* – динамическая критическая


нагрузка в случае нагруженной оболочки. Данные рисунка 1
получены для оболочки, предварительно нагруженной
статическим внешним давлением q, не превышающим величину
критического статического давления qСТ. Пунктирная линия
соответствуют скорости нагружения P0 = 490 МПа/c, сплошная
линия – скорости нагружения P0 = 880 МПа/c.

21
Р* /Рк
1,0

0,75

0,5

0,25

0,0 q/qСТ
0,0 0,25 0,5 0,75 1,0

Рисунок 1

22
Анализ полученных результатов показывает, что
предварительное нагружение статическим внешним давлением
снижает величину динамической критической нагрузки, однако
это влияние уменьшается с возрастанием скорости нагружения.
Аналогичные результаты получаются, если оболочка
предварительно нагружена статическим осевым сжатием [8].

3.4. Прямоугольная пластина при динамическом


сжатии

Пусть прямоугольная пластина равномерно сжата в одном


направлении линейно изменяющимися во времени силами
(рисунок 2). Линеаризованное уравнение динамики в данном
случае имеет вид [6]

∂2w
D∇ 4 w − T10wxx − T20wyy − 2S10wyy + ρ h = 0,
∂t 2

где Т20 = S10 = 0, Т10/h = - q0 t, q0 – скорость нагружения.

а) В случае шарнирного опирания всех сторон пластины


аппроксимирующую функцию прогиба можно выбрать в виде

mπх nπу
w = ξ (t ) sin sin
a b .

Интегрируя уравнение по пространственным координатам


методом Бубнова, приходим к уравнению (6), где

tк2 ⎡ ⎤
λ (τ ) =
D
− m12 q0tкτ + (m12 + n1 ) 2 ⎥ , m1 = mπ a, n1 = nπ b .
2

ρ⎣ ⎢ h ⎦

23
х
q=q0 t

a
b

q=q0 t

Рисунок 2

24
Уравнение типа (18) будет иметь вид

D 2 bρ
m12 q0 tк = (m1 + n1 ) 2 + 12 .
2

h tк

В работе [6] указано, что при одностороннем сжатии пластина


теряет устойчивость с образованием одной полуволны вдоль у,
поэтому n = 1 (n1 = π/b, n1 << m1), и уравнение может быть
записано в виде

D bρ b
q0 tк = (m12 b 2 + 2) + 12 2 , b= .
hb 2
m1 tк π (28)

Условие минимума нагрузки по волновому числу m1 примет


вид

2D 2 2
q0 tк = (m1 b + 1) . (29)
hb 2

Известное решение статической задачи приведено в [6]

π2
2
4D ⎛h⎞
q = 2 = q* E ⎜ ⎟ ,
к
q* =
3 (1 −ν 2 ) , m = 1.
ст (30)
b h ⎝b⎠

Решая совместно уравнения (28) и (29) с учётом (30), можно


придти к приближенному уравнению для коэффициента
динамичности

b1
K д4 − K д3 − = 0, (31)
S

25
2
⎛ EVπ h3 ⎞
где K д = q0tk к
q , S = q ⎜⎜ *3

4 ⎟ .
ст
⎝ q0 b ⎠

Численное решение уравнения (31) приводит к результатам,


совпадающим с данными работы [6]
при S = 1: Kд = 5, т = 3,
при S = 10: Kд = 3, т = 2,

число т определяется из (29) с учётом (30), т.е. из условия

2
⎛a⎞
m = (2 K д − 1) ⎜ ⎟ .
2
(32)
⎝b⎠

При достаточно высоких скоростях нагружения (S < 1)


b1=1300 и из (31) следует
1/ 2
⎛ q b b3 ⎞
1/ 4
⎛b ⎞
(1 −ν ) 27b1
4
Kд ≈ ⎜ 1 ⎟ = с1 ⎜⎜ 0 ⋅ 3 ⎟⎟ 2 3/ 4
, с1 = ≈ 1,5 . (33)
⎝S⎠ ⎝ EV h ⎠ π2

По найденному значению Kд с помощью (32) определяется


волновое число т, а затем с помощью (29) динамическая
критическая нагрузка q0 tк .
Значение константы с1 в формуле (33) рекомендуется
уточнить экспериментально или с помощью данных, имеющихся в
литературе [6].

б) В случае, когда пластинка жёстко защемлена вдоль


незагруженных сторон, а нагруженные стороны закреплены
шарнирно, аппроксимирующая функция прогиба может иметь вид

26
mπx 2 πy
w = ξ (t ) sin sin
a b .

Коэффициент λ(τ) в уравнении (6) будет

λ (τ ) =
tк2
3ρ b 4
⎡ D
( ⎤
)
⎢⎣− 3m1 b q0tкτ + h 3m1 b + 8m1 b + 16 ⎥⎦ ,
2 4 4 4 2 2

m1=mπ/a,
что приводит к характеристическому уравнению

D ⎛ 4 4 8 2 2 16 ⎞ b1ρ b 4
m12 b 4 q0 tк = ⎜ m1 b + m1 b + ⎟ + 2 . (34)
h⎝ 3 3⎠ tк

Решение статической задачи:

2
⎛h⎞
2
4 ⎛a⎞ 7D 7
т =
2
⎜ ⎟ q ≈ 2 = q* E ⎜ ⎟ .
к
ст
3⎝b⎠ , b h 4 ⎝b ⎠ (35)

К уравнению (34) добавляется условие минимума нагрузки по


волновому числу т1

2D ⎛ 2 2 4 ⎞
q0 t к = ⎜ т1 b + ⎟ . (36)
b 2h ⎝ 3⎠

Совместное решение уравнений (34) и (36) с учётом (35)


приводит к уравнению для коэффициента динамичности

4
16 3 128 2 b1 ⎛ 4 ⎞
K д4 − Kд − Kд − ⎜ ⎟ = 0 , (37)
21 441 S ⎝7⎠

27
которое решается численно.

При высоких скоростях нагружения (S ≤ 1) можно положить


1/ 4
4⎛b ⎞
Kд = ⎜ 1 ⎟ . (38)
7⎝ S ⎠

Сравнивая выражения (38) и (33) с учётом (35) и (30), можно


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

Таблица 4
Результаты численных расчётов для прямоугольной пластины
Скорость
нагружения
q0 ·10-5 Па/с 0 100 1000 2000
Условие
закрепления
Шарнирное q0tк ·10-5 Па 2,81 7,45 21,6 30,24
опирание т 1 2 4 5
Жёсткое q0tк ·10-5 Па 4,92 8,26 22,2 30,4
защемление т 2 2 4 5
Свободный край q0tк ·10-5 Па 1,00 6,44 20,4 28,81
т 1 2 4 5

В таблице 4 приведены результаты численных расчётов


пластинки с параметрами: а = 0,1 м, b = 0,2 м, V = 5 * 103 м/сек,
Е = 7,56 * 1010Па, h = 2 * 10-4м, ν = 0,3 для различных условий
закрепления сторон, при этом в случае, когда сторона у = b
свободна, а остальные стороны закреплены шарнирно,
аппроксимирующая функция прогиба выбирается в виде

mπx
w = ξ (t ) sin y.
a
28
3.5. Круговая усеченная коническая оболочка
при внешнем давлении

Линеаризованные уравнения динамики оболочки (рисунок 3)


имеют вид [1]

1 ⎛ ∂2 w ⎞ Ehctg 2γ ∂ 2 F ∂2 w
2 2 2
⎜ 2

D∇ ∇ w + 4 ctg γ ⎜ 2 2 + ctg γ ⋅ w ⎟ − − −Т10 2 −
r ⎝ ∂φ1 ⎠ r ∂r 2 ∂r

⎛ 1 ∂w 1 ∂2 w w⎞ ∂2 w
− Т10 ⎜⎜ + 2 + сtg 2
γ ⎟ + ρ h = 0,
⎝ r ∂r r ∂φ1
2
r 2 ⎟⎠ ∂t 2

∂2 w
r∇ 2∇ 2 F + = 0. (39)
∂r 2

Область применимости уравнений ограничивается


неравенством n ≥ ctg
4
1
2
γ, где n1 = n /sin γ.

Докритические усилия можно определить из соотношений

Р(t ) r tgγ
Т1 = − , T20 = − P(t ) r tgγ .
2

К системе (39) применяется процедура Бубнова, при этом для


свободно опертой оболочки аппроксимирующие функции можно
выбрать в виде

F = f (t ) e μ z cos n1φ1 sin m1 z , w = ξ (t ) e μ z cos n1φ1 sin m1 z ,

где введены обозначения

z = ln r/r0; z1 = ln r1 /r0; m1 = mπ/z,

29
r0

γ
R0
r
r1

φ
R1

Рисунок 3

30
μ – параметр, при выборе которого следует пользоваться
рекомендациями работы [1], а именно:

при z1 < 2,7 μ = 1,2,


при 2,7 ≤ z1 ≤ 3,5 μ = 1,6,
при z1 > 3,5 μ = 2.

После интегрирования система (39) сводится к уравнению (6),


коэффициент λ (τ) имеет вид

t к2 ⎡ ⎛ 4 η12 m24 ⎞⎤
=
λ (τ) ρ h r 4 ⎢ − P ( t кτ ) r1
3 2
n1 θ tg γ + D ⎜⎜ n1 θ + 4 ⎟⎟⎥ ,
1 ⎣ ⎝ n1 ⎠⎦

где приняты следующие обозначения

Е r12 ctg 2γ
η2 = m24 = (m12 + μ 2 )
2
, ,
D

η12 =
η 2 (1 − e −2 μ z 1
) (m + (μ + 1) ) (μ + 1)
2 2

(1 − e ) (m + μ ) μ ,
1
( )
− 2 μ +1 z1 2 2
1

θ=
(1 − e − ( 2 μ +1) z1
) (m + (μ + 1) ) (μ + 1)
2 2

(1 − e ) (m + (μ + 0.5) ) (μ + 0.5) ,
1
− 2 ( μ +1) z1 2 2
1

θ=
(1 − e − 2 ( μ −1) z1
) (m + (μ + 1) ) (μ + 1)
2 2

(1 − e ) (m + (μ − 1) ) (μ − 1) .
1
− 2 ( μ +1) z1 2 2
1

Уравнение типа (18) имеет вид

31
⎛ 2 η12 ⋅ m24 ⎞ b ρ hr
Р ( tк ) =
D
⎜⎜ n1 θ + ⎟⎟ + 12 2 1
r1 θ tgγ
3
⎝ n1 ⎠ θ n1 t к tgγ .
6 (40)

Условие минимума нагрузки по волновому числу n1


записывается в виде

⎛ 2 η 2 m4 ⎞
Р ( tк ) =
2D
⎜⎜ n1 θ − 1 6 2 ⎟⎟ (41)
r1 θ tgγ
3
⎝ n ⎠

Совместное решение уравнений (40) и (41) позволяет


находить значение критических нагрузок и волновых чисел для
любого закона изменения давления во времени.

а) Линейный закон изменения давления (Р = Р0 t)

Левые части уравнений (40) и (41) приобретают вид


Р(tк) = Р0 tк, и задача сводится к кубическому уравнению

(1 − 3а)⋅ (1 − а) 2
= Р0 а (42)

где а = (η1 m2 ) (n1 θ ) , P 0 = (b1ρ h r1 P0 θ tg γ ) (4D η1 m2 θ ) .


2 4 8 10 2 2 23 2 4 2

Уравнение (42) решается численно или при помощи формулы


Кардана. По найденному значению a определяется волновое
число n1, а затем из формулы (41) – динамическая критическая
нагрузка.
При достаточно высоких скоростях нагружения ( Р 0 >10)
решение уравнения (42) упрощается и можно положить

32
2
1 P 0 η12 m24 b1 ρ h r110 Р02 θ tg 2γ
а= , откуда n18 = = . (43)
Р0 θ 4 D 3θ 3

Для динамической критической нагрузки получается формула


1
2 D n12θ ⎡ E 2 h 4 P02θ ⎤ 4 b
Pдин = 3 = c2 ⎢ ⎥ , c 2 = 4 1 ≈ 4 . (44)
r1 θ tgγ 2 2
(
⎢⎣ 1 −ν V r1 θ tg γ ⎥⎦
2 2 2
) 3

Чтобы получить расчётную формулу для коэффициента


динамичности, необходимо в уравнении (40) перейти к пределу
при tк → ∞, Р0 → 0 и найти решение статической задачи

1 3
4 Dη1 2 m2 θ 4
Рст = ⋅ ,
r13θ tgγ
3
4
3

где параметр μ для величин θ , θ , m2 и η1 выбирается согласно


таблице 5.

Таблица 5
Зависимость параметра μ от z

z 0,18 0,25 0,5 0,75 1 1,25 1,5 2 2,5 3


μ=μСТ 2,05 2,1 2,3 2,45 2,5 2,6 2,65 2,7 2,8 4,0

Формула для коэффициента динамичности имеет вид


1
⎛ r14 Р0 z12 tg 2γ ⎞2
Р t
K∂ = 0 к = 3,6 с2 f∂ (z1 ) 1 −ν 2
Рст
( )1
2
⎜⎜ 3
⎟⎟ , (45)
⎝ V Eh ⎠

θ ст θ η1ст ⎛ μ 2z12 ⎞
где f ∂ (z1 ) = , f ст = ⎜1 + 2 ⎟⎟ .
θ ст θ f ст ηст ⎜⎝ π ⎠

33
Величины с индексом «ст» соответствуют статическому
значению параметра μ (таблица 5).
Формула для критического импульса получается в виде
1
Р t2 c E h2 θ 2
Iк = 0 к = 2
2 2 1 −ν 2 V r1 θ tgγ . (46)

Входящая в формулы (44), (45), (46) константа с2 уточняется


экспериментальным путём.

б) Квадратичный закон изменения давления (Р = Р1 t 2)

Левые части уравнений (40) и (41) имеют вид Р(tk) = Р1 tk2.


В этом случае для определения волнового числа получается
уравнение

n θ − P1 n + η m = 0 , P 1
8 5 2 4
=
( b Р r ρ h θ tgγ )
1 1 1
7
1
2

1 1 1 2 , (47)
D

решая которое при достаточно высоких скоростях нагружения


( P1 ≥ 100) с незначительной погрешностью находим

n13 = P1 θ . (48)

Тогда формула для динамической критической нагрузки будет


1
⎡ E 2 h 4 P1 θ ⎤ 3 b1
Pдин = P t = c3 ⎢
2
⎥ , c3 = 2
( )
3
12 . (49)
1 k 2
⎢⎣ 1 −ν 2 V 2 r12 θ tg 2γ ⎥⎦

Формула для критического импульса получается в виде

34
3 1
Р t3 c 2 E h2 θ 2
Iк = 1 к = 3 . (50)
3 3 1 − ν 2 V r1 θ tgγ

в) Кубический закон изменения давления (Р = Р2 t3)

Уравнение для волнового числа имеет вид

5
P2 r19 b1 2 (ρ h ) 2 θ tgγ
3
⎛ η 2m4 ⎞ 2 3

n ⎜⎜θ + 1 8 2 ⎟⎟
8
1 = P2 = 5 , (51)
⎝ n1 ⎠ D 2

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


получить

5
n18 = P 2 / θ 2
, (52)
1
⎡ E 3 θ h 6 P2 ⎤ 4 3
= P2 t к = c4 ⎢
3
⎥ ⎛ b1 ⎞ 8
P∂ин. c = 2 ⎜ ⎟ . (53)
(⎣ )
⎢ 1 −ν 2 2 V 3 r 3 θ 3 tg 3γ ⎥ , 4
3
1 ⎦ ⎝ 12 ⎠

Формула для критического импульса принимает вид

4 1
P t4 c 3 E h2 θ 2
Iк = 2 к = 4 (54)
4 4 1 −ν 2 V r1 θ tgγ

Сравнение формул (46), (50) и (54) показывает, что при


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

35
c0 f (z1 ) E h 2
1
E h 2θ 2
I к = c0 = , (55)
1 − ν 2 V r1 θ tgγ 1 − ν 2 V r1 tgγ

θ
где f (z1 ) = , (56)
θ

с0 – константа, подлежащая определению экспериментальным


путем.

3.6. Пологая цилиндрическая панель под действием


поперечной нагрузки

Структурные формулы для задач динамической устойчивости


могут быть получены и другим путём. Покажем это на примере
задачи динамической устойчивости пологих цилиндрических
панелей.
Уравнение движения (6) в этом случае имеет вид [3]

d 2ξ ⎛ P(t ) ⎞
+ ω 2 ⎜⎜1 + ⎟ ξ = 0,
dt 2
⎝ Pст ⎟⎠ (57)

где ω – собственная частота оболочки, соответствующая


рассматриваемой форме потери устойчивости, Р(t) – параметр
нагрузки, Рст – критическое значение нагрузки при статическом
нагружении.
После перехода к безразмерной переменной τ = tω,
уравнение принимает вид

d 2ξ ⎛ P(t ) ⎞
+ ⎜ 1 + ⎟ ξ = 0.
dτ 2 ⎜⎝ Рст ⎟⎠ (58)

36
Анализ уравнения (58) показывает, что структурная формула
для коэффициента динамичности устанавливает связь между
динамической неустойчивостью, колебаниями и устойчивостью
при статических нагрузках. При этом оказывается, что изменение
формы, размеров оболочки, условий закрепления и т. п. влияет на
процесс динамической потери устойчивости лишь постольку,
поскольку при этом изменяется спектр собственных частот и
критические нагрузки в статике. Иначе говоря, влияние
конкретных свойств геометрии оболочки на её устойчивость при
динамических нагрузках в достаточной мере учитывается такими
её интегральными характеристиками, как ω и Рст. Поэтому
формула для коэффициента динамичности может быть записана в
виде

K = Ф(ω, Рст), (59)

где Ф – функция, определяемая путём обработки эксперимен-


тальных данных.

3.7. Произвольный закон нагружения

Для приближенных расчётов при произвольном законе


нагружения предлагается графический метод определения
критического времени tк, т. е. момента потери устойчивости.
Поскольку геометрический смысл величины критического
импульса – это площадь криволинейной трапеции, ограниченной
осью t, прямыми t = 0 и t = tк и кривой зависимости нагрузки
от времени Р(t) (рисунок 4), то на основе предыдущих
рассуждений делается вывод, что для импульсов различной
формы эти площади примерно одинаковы.

37
Р(t)

P(t)

J к∧
I к∧

0 tк t

Рисунок 4

q0 t

Рисунок 5

38
Для определения критического времени tк достаточно найти

величину критического импульса I к для линейного закона
нагружения, построить график функции Р(t) и подобрать
графически величину tк таким образом, чтобы площадь

криволинейной трапеции оказалась равной I к .

3.8. Круглая пластина при динамическом сжатии

В случае, когда сплошная круглая пластина радиуса R


сжимается усилиями q0 t, равномерно распределенными по
контуру (рисунок 5), имеем

Tr 0 = Tφ 0 = −q0 t h , Trφ = 0 ,

где r, φ − полярные координаты.

Линеаризованное уравнение динамики в рассматриваемом


случае имеет вид [11]

∂2w
DΔΔw + q0t Δw + ρh = 0, (60)
∂t 2

где Δ − оператор Лапласа в полярной системе координат.


Аппроксимирующая функция прогиба выбирается согласно
рекомендациям работы [11]

w = ξ (t ) wmn (r ) cos nφ , (61)

где функция wmn (r) зависит от граничных условий.

39
а) Если пластина защемлена по контуру, то должны
выполняться условия

d wmn
wmn = =0 при r = R,
dr

согласно которым можно положить

⎡ r ⎛ r ⎞n ⎤
wmn = J n (α n +1, m ) ⎢ − ⎜ ⎟ ⎥
⎢⎣ R ⎝ R ⎠ ⎥⎦ ,
где Jn (r) − функция Бесселя n-го порядка, α n+1, m − m-й корень
уравнения

J n +1 (r ) = 0,

m − число полуволн вдоль радиуса пластины, α 1, 1 = 3,83.


Интегрирование уравнения (60) по Бубнову приводит к
характеристическому уравнению, аналогичному уравнению (18)

D b1 ρ R 2
q0 tк = α n2,m + .
R 2 h α n2,m tк2 (62)

В предельном случае статического нагружения (tk → ∞)


справедлива формула

D D
qстk = α 21,1 2
= 14,68 2 . (63)
Rh Rh

40
Минимизируя в (62) критическое время по волновому числу
α n,m и учитывая соотношение (63), получаем при b1 = 1300
формулы для критической нагрузки, коэффициента
динамичности и критического импульса

1/ 2 2 1/ 2
⎛ q Eh ⎞ ⎛R⎞ ⎛q h⎞
q0 tк = 4,71⎜ 0 ⎟ , K д = 3,5 ⎜ ⎟ ⎜ 0 ⎟ .
⎝ V ⎠ ⎝ h ⎠ ⎝ EV ⎠

⎛ Eh ⎞
I k = 11,1⎜ ⎟.
⎝V ⎠

б) Если контур пластины оперт, должны выполняться


условия

d 2 wmn ν dwmn
wmn = + =0 при r = R,
dr 2 r dr

и функция wmn (r) может быть выбрана в виде

⎡ r ⎛ r ⎞n ⎤
wmn = J n ( β n, m ) ⎢ − ⎜ ⎟ ⎥,
⎣⎢ R ⎝ R ⎠ ⎦⎥

где β n,m − корень уравнения

r J n (r ) − (1 −ν ) J n +1 (r ) = 0,

β 0,1 = 2,047, ν = 0,3.

41
Интегрирование уравнения (1) по Бубнову приводит в этом
случае к характеристическому уравнению

D b1 ρ h R 2 ( β n,m − 2 − 2ν )
2 2
q0 tк = β n2,m + ,
R2 β n2,m tк2 ( β n2,m − 1 +ν 2 )

из которого при достаточно больших скоростях нагружения


получаем

1/ 2
D D ⎛ q Eh ⎞
q =β
k
ст
2
= 4,19 2 , q0 tк = 3,96 ⎜ 0 ⎟ ,
⎝ V ⎠
0,1 2
Rh Rh
⎛ Eh ⎞
I k = 77,84 ⎜ ⎟.
⎝V ⎠

Анализ полученных формул показывает, что при высоких


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

3.9. Кольцевой сектор при равномерном сжатии

Рассматриваемая задача (рисунок 6) возникает в процессе


производства конических оболочек методом «свертывания».
Для самоуравновешенности внешних сил необходимо на
контуре дуги радиуса R2 приложить сжимающие усилия

R2 − R1
qr = qφ , qφ = q0 t.
R2

42
qr

R1

qφ R2

2φ 0

Рисунок 6

43
При условии шарнирного опирания на радиальных кромках
φ = ± φ0 аппроксимирующая функция прогиба может быть
выбрана в виде

w = ξ (t ) r sin n1 φ , n1 = nπ / φ0 ,

что, как показано в [11], приводит после интегрирования по


времени и координатам к характеристическому уравнению

1
2
[
q0 tк (R1 − R2 ) n12 (R1 + R2 ) + (R1 − R2 ) = ]
( ) ⎞ + b ρ (R − R ).
4 4
n1 − 1 ln ⎛⎜ 1
D 2 2 R
= ⎟
1 1 2
h ⎝ R2 ⎠ t2
к

Минимизируя в этом уравнении критическое время tk по


волновому числу n1, получаем

( )
4 D n12 − 1 ln ⎛⎜ 1 ⎞⎟
R
⎝ R2 ⎠ .
q0 tk =
(
h R1 − R22
2
)
Для критической статической нагрузки получаем формулу

( 2
)
2 D n12 − 1 ln ⎛⎜ 1 ⎞⎟
R
⎝ R2 ⎠
q стk = , n1 = π
[
h (R1 − R2 ) n1 (R1 + R2 ) + (R1 − R2 )
2
]
2φ0
,

из которой для пластины с углом 2φ0 = π следует вывод о


невозможности перехода в изогнутое состояние под равномерной
статической нагрузкой. В динамике процесс «свертывания» при

44
таком угле возможен, так как при достаточно высоких скоростях
нагружения справедлива формула [11]

⎢ ln ⎜⎝ R2 ⎟⎠ (R1 + R2 )⎥
1/ 4

1/ 2
⎡ ⎛ R1 ⎞ 2 2 ⎤
⎛ q Eh ⎞
q0 tк = 3,96 ⎜ 0 ⎟ ⎢ ⎥ .
⎝ V ⎠ ⎢ R12 − R22 ⎥
⎣ ⎦

3.10. Кольцевая пластина при динамическом сдвиге

В случае, когда сдвигающее усилие линейно зависит от


времени, линеаризованное уравнение динамики имеет вид [11]

2 Trφ 2 ∂ ⎛ 1 ∂w ⎞ ∂2w
DΔ Δw − 2 r0 ⎜ ⎟ + ρ h 2 = 0,
r ∂r ⎜⎝ r ∂φ ⎟⎠ ∂t

где Trφ = T0 t − погонное сдвигающее усилие в сечении r = r0

( r0, r1 − внутренний и внешний радиусы пластины), T0 −


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

πz
w = ξ (t ) sin sin n (φ − kz),
z1

где z = ln (r/r0), z1 = ln (r1/r0), n − число полуволн вдоль


окружности, k − тангенс угла наклона складки к радиусу кольца.

45
Интегрирование по Бубнову приводит к следующему
характеристическому уравнению [11]

n12 k Sk = n14 (k 2 + 1) + 6n12 k 2 + 1 + 2n12 + α −2 +


2

( )
+ α −1 2 + 2n12 k 2 − 2n12 + S02 Sk−2 , (64)

где введены обозначения

1/ 2
⎛ z1 ⎞
( )
4
2 T0 tk r02 T0 ⎛ r0 ⎞ 4 z1
Sk = , S = ⎜ ⎟ e − 1 ⎜⎜ 2 ⎟ ,
2 ⎟
α 2D + π
0
EV ⎝ h ⎠ ⎝ 1
z ⎠
n12 = n 2α − 2 , α = π / z1.
Минимизируя в (64) критическое усилие Sk по волновому
числу n1 и параметру k, придем к уравнениям

( )
n12k Sk = 2 n14 k 2 + 1 + 6 n12k 2 + 2 n12α −1 k 2 − 1 ,
2
( )
= 4 n (k + 1) + 12 n k + 2 n12α −1.
4
n12k Sk 4
1
2 2 2
1

Решая полученную систему уравнений, получим зависимости


усилия Sk и волнового числа n1 от параметра k

α 2 +1 k
Sk = 8 ⋅ 2 ⋅ ,
α 1− k2

n = n1 α =
2 2 2 (3k 2 − 1)α 2 + (k 2 + 1)
.
1− k4

46
Подставляя полученные зависимости в уравнение (64), придем
к уравнению для определения параметра k

( ) ( ) (
k 4 2k 2 − 1 + α −2 2k 2 + k 4 − 1 k 2 + α −4k 4 = S02 1 − k 8 , )
S0 α 2
S0 =
(
16 1 + α 2
,
)
которое решается численно.

Таблица 6
Результаты численных расчетов

r1/r0 S0 0 1 3 10

Tk 11,34 13,72 18,58 29,87

4/3 n 9 11 14 19

k2 0,50 0,56 0,65 0,76

Tk 11,41 13,82 18,98 29,97

2 n 4 5 6 8

k2 0,49 0,55 0,64 0,75

Tk 11,88 14,36 19,91 32,01

4 n 2 2 3 4

k2 0,45 0,52 0,62 0,74

В таблице 6 приведены результаты расчетов критических


усилий, волновых чисел и угла наклона вмятин для пластин с
различным отношением внешнего радиуса кольца к внутреннему
при некоторых скоростях нагружения. При малых скоростях
47
нагружения полученные результаты совпадают с
соответствующими значениями для статической задачи, которая
подробно исследовалась в [16]. При малых значениях отношения
r1/r0 ( ln (r1/r0) < 2 ) результаты совпадают с соответствующими
значениями для задачи устойчивости прямоугольной пластины
при динамическом сдвиге.

3.11. Подкрепленная цилиндрическая оболочка

Для решения задач устойчивости подкрепленных оболочек


существует два подхода. В первом случае подкрепленную
оболочку рассматривают как однородную ортотропную. Другой
подход заключается в том, что отдельно интегрируются уравнения
устойчивости гладкой оболочки и подкрепления, а затем
выполняются условия стыковки. В отличие от первого, второй
подход не получил широкого распространения при решении задач
динамической устойчивости. Однако с его помощью можно
получить практически важные результаты.
Пусть круговая цилиндрическая оболочка подкреплена в
середине шпангоутом и нагружена динамическим внешним
давлением. Для расчета такой оболочки можно воспользоваться
методикой, предложенной в [2], где сделаны следующие
предположения:
1. Поведение обшивки описывается уравнениями
полубезмоментной теории.
2. Центры тяжести поперечных сечений шпангоута лежат в
срединной поверхности обшивки, оси шпангоута нерастяжимы.
3. В решении учитывается только изгибная жесткость
шпангоута в своей плоскости.
4. Ширина шпангоута пренебрежимо мала по сравнению с
длиной оболочки.
48
Изгиб кругового шпангоута с нерастяжимой осью под
действием распределенной касательной нагрузки интенсивности
T* описывается уравнением

⎛ d 6V d 4V d 2V ⎞
( EI )ш ⎜⎜ 6 + 2 4 + 2 ⎟⎟ = − R 4 T * ,
⎝ dφ dφ dφ ⎠

где (EI)Ш − изгибная жесткость шпангоута, V − окружные


перемещения точек осевой линии шпангоута. Передаваемая на
шпангоут касательная нагрузка равна разности сдвигающих сил на
участках обшивки, прилегающих слева S_ и справа S+ к
шпангоуту

T * = S− − S+ .
Потеря устойчивости обшивки описывается уравнением

∂ 4Φ D ⎛ ∂ 8Φ ∂ 6Φ ∂ 4Φ ⎞
Eh + ⎜ + 2 + ⎟+
∂x 4 R 6 ⎜⎝ ∂φ 8 ∂φ 6 ∂φ 4 ⎟⎠
P0 t ⎡ ∂ 6Φ ∂ 4Φ ⎤ ρ h ∂ 6Φ
+ 3 ⎢
+ 4⎥+ 2 = 0, φ = y .
R ⎣ ∂φ 6
∂φ ⎦ R ∂φ ∂ t
4 2 R

Здесь Ф − функция перемещений, определяющая


перемещения по формулам

∂Φ ∂Φ 1 ∂ 2Φ
u= , v=− , w = .
∂x R ∂φ 2 R ∂φ 2
Решение уравнения следует искать в виде

Φ = ς (τ ) A(x )cos nφ.

49
После интегрирования по времени приходим к уравнению

⎛λ⎞
4

A − ⎜ ⎟ A = 0,
IV

⎝R⎠

λ 4=
( ) ⎡P t − ( n
R n4 n2 − 1 2
)
− 1 D b1 ρ h R ⎤
− 2
Eh


0 k
R3 ( ) ⎥.
n − 1 tk2 ⎦
Выражая касательную нагрузку, передаваемую на шпангоут,
через функцию перемещений, можно найти перемещения точек
оси шпангоута [11]. При этом должны выполняться следующие
четыре условия:
a) равенство окружных перемещений участков обшивки,
прилегающих к шпангоуту слева и справа;
b) равенство осевых перемещений;
c) равенство осевых сил в обшивке слева и справа;
d) равенство окружных перемещений точек оси шпангоута
соответствующим окружным перемещениям обшивки.
Выполнение всех стыковочных условий приводит к системе
линейных однородных уравнений [11], равенство нулю
определителя которой дает характеристическое уравнение для
определения параметра λ.
Если оболочка свободна оперта по обоим торцам, то
уравнение для параметра λ распадается на два, одно из которых
соответствует местной потере устойчивости, когда обшивка
теряет устойчивость, а шпангоут сохраняет круговую форму
(критическое давление Рм не зависит от жесткости шпангоута).
Второе уравнение соответствует общей потере устойчивости,
когда оболочка теряет устойчивость вместе с подкреплением.
Критическое давление Роб = P0 tk находится из уравнения

50
P0 tк =
E h λ4
+
D n2 −1 (
+ 2
)
b1 ρ h RP02
(
R n4 n2 −1 ) R3 n − 1 (P0tk )
2
(
,
)
в котором значение параметра λ зависит от жесткости шпангоута
и волнового числа n.
В таблице 7 приведены результаты расчетов [11] отношения
критической нагрузки общей потери устойчивости к критической
нагрузке местной потери устойчивости в зависимости от скорости
нагружения для оболочки с параметрами: R/l=1, R/h=500,
5 4
P*=936 P0R /VEh .
Во второй строке таблицы 7 приведены значения критической
жесткости шпангоута (EI)k. Если жесткость шпангоута меньше
критической, то происходит общая потеря устойчивости.
Увеличение жесткости шпангоута приводит к повышению
критической нагрузки. При достижении жесткости шпангоута
критического значения происходит смена форм устойчивости с
общей на местную. Дальнейшее повышение жесткости шпангоута
не влияет на критическую нагрузку.

Таблица 7
Расчеты для подкрепленной оболочки

P* ∙ 10 – 4 0 1 2 4 8

(EI / Dl)k 1,50 1,10 0,76 0,42 0,16

Pоб / Pм 2 1,46 1,25 1,10 1,03

Расчеты показывают, что при увеличении скорости


нагружения критическая жесткость шпангоута уменьшается, а
отношение Pоб / Pм стремится к единице.

51
3.12. Оболочка с вырезами

Рассмотрим круговую цилиндрическую оболочку,


ослабленную системой прямоугольных вырезов, расположенных
на равном расстоянии друг от друга вдоль образующей. Пусть на
оболочку действует динамическое внешнее давление.
Учитывая, что при высоких скоростях нагружения
критическая нагрузка слабо зависит от граничных условий и
предполагая, что при потере устойчивости кривизна оболочки
изменяется лишь в окружном направлении, будем судить об
устойчивости оболочки по устойчивости кольца, мысленно
выделенного двумя перпендикулярными к оси оболочки
сечениями. Ширина b кольца равна расстоянию между центрами
вырезов по образующей. Линеаризованное уравнение динамики в
этом случае имеет вид [11]

4 ∂ w
2
EI 2
∇1 w + P0 t R ∇1w + ρ hR
3
= 0, (65)
b ∂t 2
∂2
∇1 = 2 + 1, t = tkτ , w = w (φ )ζ (τ ).
∂φ

К решению уравнения (65) применим подход, предложенный в


[11], согласно которому сначала проинтегрируем его по времени

EI 2 ρ h R4w
∇1 w + P0 tk R3 ∇1w + = 0.
b tk2
Полученное уравнение с переменными коэффициентами
можно решить методом последовательных приближений [15].
Нулевое приближение выберем в виде

52
⎡ P t R3b ρ hR 4b ⎤
⎥ × ∑ Ai ri e , w = ∑ Ai ri e .
ri φ k ri φ
w = − ⎢1 + 0 k + (k )

⎣ EI EI tk2 ⎦

Здесь r1 = –r3 , r2 = –r4 – корни характеристического уравнения

r 4 + a2 r 2 + a0 = 0, (66)
2π 2π
1 ⎛ P0tk R 3b ⎞ 1 ⎛ P0tk R 3b ρ hR 4b ⎞
a0 =
2π ∫0 ⎜⎝
⎜ 2 +
EI ⎠
⎟⎟ dφ , a2 =
2π ∫0 ⎜⎜⎝ EI + EI tk2 ⎟⎟ dφ.

Функция w(φ ) должна удовлетворять условию симметрии.


Если оболочка теряет устойчивость по n волнам, то для
определения коэффициентов Ai имеем условия

dw d 2w
π = 0, = 0.
dφ φ =
n dφ 2 φ =
π
2n

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


придем к характеристическому уравнению вида

n4 − a2n2 + a0 = 0.
Если вырезы отсутствуют, то характеристическое уравнение
имеет вид

bρ h R 4
( 2
)
EI n − 1 − P0 tk R b n − 1 +
2 3
( 2
) = 0.
tk2
Минимизируя в этом уравнении критическую нагрузку по
волновому числу, получаем

53
1/ 4
⎛ ρ EI P02 ⎞ 3EI
Pд = P0 tk = 7,2 ⎜⎜ 2
⎟⎟ , Pстk = 3 .
⎝ R ⎠ Rb

Пусть оболочка перфорирована квадратными вырезами a∗a


(a << R), расположенными симметрично (k1 – количество вырезов
вдоль окружности). Тогда окружность рассматриваемого кольца
можно представить в виде нескольких участков с различной
жесткостью. Участок с вырезом характеризуется жесткостью
(
EI1 = EI 1 − a
b
).
Коэффициенты в уравнении (66) примут вид

P0tk R 3b P0tk R 3b ρ hR 4b
a0 = 2 + , a2 = −1 + + ,
EIδ EIδ EI tk2
−1
⎡ a 2 k1 ⎤
δ =⎢ + 1⎥ .
⎣ 2π R (b − a ) ⎦

После преобразований получаем значения критических


нагрузок (статической и динамической) для оболочки с вырезами

Pст = Pст δ , Pд = Pд δ 1/ 2 ,

т.е. отношение критической нагрузки для оболочки с вырезами к


критической нагрузке для оболочки без вырезов в динамике равно
корню квадратному из аналогичного отношения в статике. В [11]
показано, что, в отличие от внешнего давления, при осевом
сжатии влияние малых вырезов в статике и динамике
практически одинаково. В связи с приближенной постановкой
задачи эти выводы носят лишь качественный характер и требуют
уточнения с помощью эксперимента.
54
4. МЕТОДИКА ЭКСПЕРИМЕНТА И ЕГО РЕЗУЛЬТАТЫ

4.1. Общие положения

Эксперимент проводится с целью уточнения констант,


входящих в расчётные формулы. Методика эксперимента описана
для круговой конической оболочки и пологой цилиндрической
панели.
Границы устойчивости при динамическом нагружении
оболочек описываются уравнением гиперболы [4]

⎛ Рфэ ⎞ ⎛ I фэ ⎞
⎜ − 1⎟ ⎜ − 1⎟ = 1
⎜Р ⎟⎜ ⎟ ,
⎝ ст ⎠ ⎝ I к ⎠

э
с асимптотами Рст и Iк, где Рф - давление во фронте волны, I фэ −
импульс давления. В качестве Рст может быть принято
критическое значение при статическом внешнем давлении,
полученное экспериментально, или близкая к нему теоретическая
величина [17].
Построенная на основе экспериментальных данных кривая
устойчивости позволяет определить асимптотическое значение Iк,
т. е. величину критического импульса давления при высоких
скоростях нагружения. Эта величина сравнивается со значением,
рассчитанным по формуле (55).

4.2. Требования к проведению испытаний

Установление величины критического импульса связано с


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

4.3. Форма, размеры и методы отбора образцов

Анализ работ, в которых представлены результаты


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

56
2) материал должен допускать использование сравнительно
простого технологического процесса изготовления моделей (это
обстоятельство особенно важно при массовых испытаниях);
3) сортамент материала должен быть промышленного
изготовления в диапазоне толщины 0,1-0,5 мм.
На основании изложенного, для изготовления оболочек
использовался материал 20Х, для изготовления панелей −
материал Х18Н9Т.
Геометрические параметры оболочек выбирались в диапазоне
размеров, удобных для экспериментальных исследований и были
обусловлены геометрическими характеристиками ударной трубы.
Все образцы имели радиус нижнего основания R = 80 мм.
Толщина листа материала оболочек определялась с помощью
вертикального оптиметра ИЗВ-2. Для этого вырезались полоски
из листа материала вдоль и поперек проката размером 400*40 мм
на различных участках, в центре которых проводились замеры
толщины. Значение толщины определялось как среднее от трех
размеров в каждой точке.
Заготовки для оболочек вырезались из листовой стали по
специальным шаблонам, что позволило выдерживать
необходимую точность геометрических размеров
изготавливаемых образцов. Поскольку в дальнейшем измерялись
начальные погиби в различных точках поверхности оболочки, то
на поле заготовки наносилась координатная сетка. Это позволяло
не только отыскивать необходимые точки, но и контролировать
изготовление оболочек.
Оболочки изготавливались путем свертывания на специально
изготовленных приспособлениях, которые представляют собой
жёсткий конус с продольным сквозным пазом в зоне,
соответствующей сварному шву оболочки. Паз предназначен для
57
подвода сварочных электродов. Наружная поверхность конуса
тщательно полировалась до размера, соответствующего
внутреннему диаметру оболочки. Заготовка обтягивалась по
поверхности конуса, обматывалась резиновым жгутом и
подвергалась точечной сварке на машине ТКМ-7. Съём
изготовленных образцов производился взаимным смещением
оболочки и жесткого конуса при помощи специального упора.
Последней операцией технологического процесса являлось
шлифование торцов оболочек. При этой операции устраняются
погрешности и дефекты резки, и других предшествующих
операций. Шлифование торцов производилось на поверочной
плите, обтянутой наждачной бумагой. Цилиндрические панели
прямоугольной формы в плане, изготавливались из листовой
стали толщиной 0,12-0,22 мм. Параметр кривизны панелей
варьировался в пределах от 40 до 120 и удлинение от 1 до 25.
Панели крепились с помощью эпоксидного клея и точечной
сварки к жёстким оправкам, изготовленным из точеной стальной
трубы диаметром 216 мм с толщиной стенки 6 мм. Модуль
упругости определялся динамическим способом.

4.4. Испытательное оборудование

Экспериментальные исследования при статическом и


динамическом нагружении внешним давлением проводились на
универсальной экспериментальной установке.
Установка объединяет в себе необходимые нагрузочные и
регистрирующие устройства и приборы с постоянной
электромеханической связью и позволяет проводить испытания
образцов и в статике, и в динамике при одинаковых граничных
условиях. Принципиальная схема установки с измерительной
аппаратурой приведена на рисунке 7, она состоит из камер
58
высокого 9 и низкого 4 давления. В качестве рабочей среды от
компрессора 10 в камеру высокого давления через редуктор
попадает воздух, давление которого контролируется по манометру
8 с ценой деления 0,01 МПа. Ударная волна создается внезапным
разрывом диафрагмы, разделяющей камеры высокого и низкого
давления.
При разрыве диафрагмы производится быстрый ввод газа в
канал камеры низкого давления, где формируется ударная волна,
движение которой направлено в сторону среза ударной трубы. На
срезе ударной трубы при помощи специальных приспособлений
крепятся испытуемые образцы. Для обеспечения возможности
более широкого варьирования законом нагружения изготовлена
специальная приставка, позволяющая преобразовать ступенчатый
закон нагружения в линейный и квадратичный.
Измерительный комплекс приборов настроен для измерения
скорости движения падающей волны, определения давления во
фронте ударной волны, измерения величин деформации в
определенных точках оболочки, изучения характера
волнообразования, перемещений (прогибов) в трех точках панели
и осреднённого прогиба в среднем сечении панели. Скорость
движения ударной волны замерялась по временному интервалу
прохождения волной базового расстояния l = 1800 мм между
двумя датчиками 12, идентичными по своим характеристикам.
При прохождении фронта ударной волны контактные датчики
замыкают цепи подачи сигналов соответственно на запуск и
включение частотомера 43-33 (3), который работает в режиме
счёта времени.
Для замера давления во фронте ударной волны
использовались кварцевые пьезометрические датчики давления 17
типа РНД 5/18(2).

59
Σ С1-74 ФА

С1-18 ФА
3 4 5 6 7 8
К С8-9Д

У У У У БС 43-33 БП

ДП
ДП
ДП
1 БА
ДП
9
17
ФА ПМ-1
16
2 15 14 13 12 11 10

Рисунок 7

60
Один устанавливался на расстоянии 5 мм от поверхности
оболочки, другой был вмонтирован в корпус обтекателя.
Сигнал от датчика давления подавался на экран
пьезоизмерительного устройства РМ-1 (1) и с него проводилось
фотографирование. Было разработано устройство для временного
и дистанционного управления затвором фотоаппарата. Кроме
того, для измерения давления применялись датчики ёмкостного
типа с чувствительной мембраной диаметром 5 мм (1-я
резонансная частота 60 кгц). Процесс выпучивания оболочек при
динамическом нагружении исследовался с помощью скоростной
киносъёмочной камеры СКС-1 (14). Съёмка проводилась при
искусственном освещении с частотой 4000 кадр/с.
Для регистрации динамических деформаций оболочек
применялись проволочные тензометры сопротивления с базой
10 мм, наклеенные в сечении, которое соответствует, как было
установлено в предварительных опытах, наибольшему прогибу.
Запись показаний электрических тензометров сопротивления,
усиленных прибором ТОПАЗ-03 (13) фиксировалась электронным
двухлучевым запоминающим осциллографом С1-20 (11).
Регистрация прогиба в отдельных точках позволяет судить о
форме потери устойчивости, а кривая осреднённого прогиба
используется для определения момента потери устойчивости.
Для регистрации прогибов применялись бесконтактные
ёмкостные датчики, работающие на принципе частотной
модуляции [13]. Высокая чувствительность датчиков позволила
использовать их с электродами малых размеров и тем самым
добиться приемлемой локализации зон измерения. В качестве
регистрирующей аппаратуры использовались серийные
осциллографы С8-9А, С1-18 и С1-74. На рисунке 7 приняты
следующие условные обозначения: БД – блок датчиков, БС – блок

61
синхронизации, ДД – датчик давления, ДП – датчик перемещения,
К – коммутатор, ПД – преобразователь формы импульса, У –
усилитель, ФА – фотокамера, ∑ - сумматор аналоговых сигналов.
До начала испытаний проводилась тарировка всего
измерительного комплекса по всем каналам и со всеми
коэффициентами усиления. Датчик давления градуировался на
специально изготовленной тарировочной трубе.

4.5. Проведение испытаний и обработка


результатов

При определении границ устойчивости для каждой серии


конических оболочек вначале задавались заведомо безопасные
э э
величины давления во фронте Рф и импульса I ф . Далее оболочку
подвергали воздействию динамических нагрузок с возрастающей
интенсивностью до тех пор, пока не происходило выпучивание
стенок, регистрируемое с помощью скоростной киносъёмки и
визуально.
Согласно принятому критерию, критическая динамическая
нагрузка считается достигнутой в тот момент, когда стрела
остаточного прогиба становится равной примерно толщине
оболочки. Подобным прогибам, судя по данным опытов,
соответствует достижение приведёнными напряжениями предела
текучести материала.
При экспериментах импульс нагрузки I фэ изменялся от
э
0,2 до 2 КПа С, давление во фронте Рф в пределах 0,1-3 МПа.
На рисунках 8 – 11 представлены границы устойчивости в
плоскости «давление-импульс» для каждой серии оболочек.
Сплошной линией показана граница, разделяющая области
устойчивости и неустойчивости.

62
На рисунках 8 – 10 показана зависимость критического
импульса нагрузки от угла конусности при постоянном
отношении r1 / r0 = 3,55. Влияние параметра z1 = ln r1 / r0 на
величину критического импульса нагрузки при постоянном
значении 2γ = 40 0 показано на рисунке 11.
Разброс в показаниях Iк при определенном законе
нагружения для данной серии не превышал 5-10%. Стабильность
результатов при динамических испытания выше, чем это
наблюдалось при испытаниях подобных образцов на статическую
нагрузку, причём с увеличением величины давления во фронте и
уменьшением величины импульса нагрузки разброс в показаниях
уменьшается.
Как показали отдельные экспериментальные исследования,
наибольшее влияние на разброс опытных значений критических
нагрузок оказывают начальные неправильности. При
динамическом нагружении, в отличие от статического, за счёт
резкого возрастания нагрузки и инерционности деформируемой
оболочки образование выпучины в месте, где имеется начальная
погибь, очевидно, запаздывает и процесс волнообразования
происходит по всей поверхности образца. Поэтому начальные
неправильности при динамической потере устойчивости
оказывают меньшее влияние на значение критической нагрузки,
чем это наблюдается при статических испытаниях. Причём с
увеличением интенсивности нагружения это влияние
уменьшается.
То, что процесс волнообразования при потере устойчивости
происходит практически одновременно по всей поверхности
оболочки, подтверждают результаты скоростной киносъёмки
процесса выпучивания при динамическом нагружении.
Установлено, что процесс выпучивания длился 13·10-4с.

63
I фэ КПа·с
2γ=400
z1=1,260

Рст.=0,013 МПа
0,08

0,04

Рст Iк = 0,0297 КПа·с

0,04 0,08 0,12 Рфэ МПа

Рисунок 8

2γ=600
I фэ КПа·с z1=1,260

Рст.=0,020 МПа
0,08

0,04

Рст
Iк = 0,0310 КПа·с

0,04 0,05 0,12 Рфэ МПа

Рисунок 9

64
I фэ КПа·с
2γ=400
z1=0,157, Рст3=0,056 МПа
Рст3
z1=0,397, Рcт2=0,029 МПа
0,08
z1=1,045, Р cт1=0,022 МПа

Рст2
Iк1 = 0,0503 КПа∙с
0,04
Iк2 = 0,0384 КПа∙с
Рст1
Iк3 = 0,0311 КПа·с

0
0,04 0,08 0,12 Рфэ МПа

Рисунок 10

I фэ КПа·с
2γ=800
z1=1,260

Рст=0,02 МПа
0,08

0,04
э
Рст Iк = 0,0319 КПа·с

0
0,04 0,08 0,12 Рфэ МПа

Рисунок 11

65
Через 4·10-5с на оболочке наметилось 10 равномерно
расположенных вмятин. В дальнейшем вмятины увеличивались в
плане, на границах вмятин появлялись резко выраженные рёбра.
Следует отметить, что при воздействии динамической
нагрузки, близкой по величине к критической, происходила
потеря устойчивости, как в упругой, так и в пластической зоне.
На первоначальных кадрах были видны упругие выпучины,
которые перестраивались при дальнейшем нагружении.
Как показывают экспериментальные исследования, характер
волнообразования при динамическом нагружении конических
оболочек внешним давлением существенно отличается от
статического. При динамическом нагружении для всех серий
оболочек проявляются более высокие формы потери
устойчивости. С увеличением интенсивности нагружения число
волн в окружном направлении увеличивается. Выпучины в плане
имеют форму вытянутого вдоль образующей конуса
шестиугольника с более резко очерченным контуром и меньшей,
чем в статике, глубиной.
Следует отметить, что проявление более высоких форм потери
устойчивости наблюдается при воздействии нагрузок в 1,5-2,8
больших, чем при статическом нагружении. Для случаев
воздействия динамической нагрузкой, по величине близкой к
статической (квазистатической) характер волнообразования и
число волн при динамическом и статическом нагружении
одинаковы. Экспериментальные исследования конических
оболочек с различными значениями параметра r1 / r0
показывают, что с увеличением усечённости конуса
(уменьшением r1 / r0) границы устойчивости расширяются.
На рисунке 12 представлены осциллограммы для
иллюстрации характера зависимости амплитуды деформаций от

66
времени при различной интенсивности динамической нагрузки.
При этом выбраны показания пары тензорезисторов, которые
находились в области с наибольшим остаточным прогибом
(тензорезистор 1 наклеен вдоль образующей, тензорезистор 2 – в
кольцевом направлении).
Вначале оболочку нагружали безопасной нагрузкой, после
воздействия импульса нагрузки остаточных деформаций не
возникало. При воздействии критического импульса происходит
потеря устойчивости оболочки «хлопком». Деформации меняют
свой знак и имеют значительные величины.
Испытание цилиндрических панелей проводилось следующим
образом:
1) После закрепления панели на оправке с помощью
лекальной линейки контролировалось отсутствие искажений
цилиндрической формы вдоль образующей.
2) Для каждой панели экспериментально определялись
собственные частоты до пятой гармоники и формы свободных
колебаний (использовался метод перемещающейся опоры [21]).
Это испытание служило дополнительным элементом контроля на
отсутствие начальных неправильностей и дефектов монтажа
панели.
3) Оправка с панелью закреплялась в камере нагружающего
устройства, и проводилось её испытание на потерю устойчивости
с регистрацией на фотоплёнку осциллограмм давления и
прогибов. «Визуальный» контроль потери устойчивости
осуществлялся с помощью запоминающего осциллографа С8-9А.
Определение критической нагрузки Рст также оказалось
возможным проводить на описанной установке путём создания
предельно низкого темпа нарастания давления. Поскольку время
от начала нагружения до потери устойчивости составляло при

67
Рисунок 12

68
этом 5-6 мс, что примерно на порядок больше периода
собственных колебаний панели, то динамические эффекты
практически не проявлялись, то есть нагружение носило
квазистатический характер.
На рисунках 13-16 приведены характерные осциллограммы
для основных типов нагружения:
рис. 13 – квазистатическое нагружение, 2 мс/дел.;
рис. 14 – нагружение ступенчатым импульсом, 0,5 мс/дел.;
рис. 15 – линейный закон нагружения, 0,2 мс/дел.;
рис. 16 – квадратичный закон нагружения, 0,2 мс/дел.
Результаты проведённых экспериментов отражены на
рисунках 17-19.
Полученные экспериментальные данные позволяют сделать
некоторые выводы качественного характера:
1) Получено подтверждение высказанного ранее
предположения о том, что геометрия оболочки в достаточной мере
учитывается интегральными характеристиками ω и Рст.
Действительно, экспериментальные точки, особенно на рисунках
17 и 18, довольно хорошо группируются возле общих кривых,
несмотря на то, что они получены для панелей различного
удлинения и с различными значениями параметра кривизны.
2) При квазистатическом нагружении развитие вмятины
начинается, как правило, по несимметричной форме. При этом
отдельные точки панели совершают колебательные движения со
сдвигом по фазе, так что центр вмятины как бы перетекает по
поверхности панели от одного края к другому и обратно. Этот же
эффект наблюдался и при нагружении ступенчатым импульсом в
том случае, когда давление лишь незначительно превосходило
Рст. Отмеченное явление иллюстрируется на рисунке 20, на кото-

69
Рисунок 13

Рисунок 14

70
Рисунок 15

Рисунок 16

71
Рисунок 17 Рисунок 18

Рисунок 19

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

4.6. Определение констант и функций

Для определения константы с0 в зависимости (55)


используются границы устойчивости в координатах "давление-
импульс", построенные с учётом экспериментальных результатов.
Используя эти границы, можно экспериментально определить для
каждой серии оболочек значение Iк. Сравнивая эти величины со
значениями критического импульса для данной серии,

73
полученного по зависимости (55) можно получить значение
коэффициента c0. Как показал эксперимент, с учетом коррекции
это значение необходимо брать в пределах 6,03 < c0 < 6,29.
Для цилиндрических панелей на основе экспериментальных
данных построены формулы для коэффициента динамичности K
при линейном и квадратичном законах нагружения, а также для
критической длительности импульса tк при ступенчатом законе
нагружения.
При расчёте панелей на динамическую устойчивость
необходимо, прежде всего, определить – расчётом, по справочным
данным или экспериментально – низшую собственную частоту
панели ω и критическое давление при статическом нагружении
Рст. Затем, задавшись величиной скорости нарастания нагрузки q0,
определить соответствующее значение параметра α и по
приведенным на рисунках 17-19 эмпирическим формулам
вычислить значение коэффициента динамичности.
Величина динамической критической нагрузки определяется
как

Р ин = K Рст .

В случае ступенчатого импульса вместо скорости нагружения


следует задать величину одного из параметров импульса – Р
(амплитуду) или tк (длительность). Тогда критическое значение
другого параметра можно определить по эмпирической формуле

−0 , 36
1,1 ⎛ P ⎞
tk = ⎜⎜ ⎟⎟
ω ⎝ Pст ⎠ .

74
Рисунок 20 ( а) Рисунок 20( б)
а:Δt = 0,5 мс б: Δt = 0,25 мс

Рисунок 20

75
ЛИТЕРАТУРА

1. Агеносов Л.Г. Устойчивость и свободные колебания


тонких цилиндрических и конических оболочек кругового сечения
при различных краевых условиях. [Текст] / Л.Г. Агеносов, А.В.
Саченков // Исслед. по теории пластин и оболочек. – 1964. –
Вып. 2. – С. 111-126.

2. Алфутов Н.А. Основы расчета на устойчивость упругих


систем. [Текст]: монография / Н.А. Алфутов. – Москва,
Машиностроение. – 1978. – 312 с.

3. Амиро И.Я. К определению критических значений быстро


возрастающих во времени сжимающих сил. [Текст] / И.Я. Амиро
// Прикладная механика. – 1979. – Т. 15, № 5. – С. 54-60.

4. Андерсон Д. Динамическое выпучивание цилиндрических


оболочек под действием нестационарного бокового давления.
[Текст] / Д. Андерсон, Х. Линдберг // Ракетная техника и
космонавтика, – 1968. – Т. 6, № 4. – С. 332-347.

5. Баженов В.Г. Экспериментально-теоретическое исследо-


вание нелинейного деформирования и потери устойчивости
оболочек вращения при изгибе. [Текст] / В.Г. Баженов, А.И.
Кибец, М.В. Петров, Т.Г. Федорова, Д.В. Шошин, А.А. Артемьева
// Проблемы прочности и пластичности. – 2010. – Вып. 72. –
С. 80-84.

6. Вольмир А.С. Нелинейная динамика пластинок и


оболочек. [Текст]: монография / А.С. Вольмир. – Москва, Наука.
– 1972. – 432 с.

7. Григолюк Э.И. Устойчивость оболочек. [Текст]:


монография / Э.И. Григолюк, В.В. Кабанов. – Москва, Наука. –
1978. – 359 с.

76
8. Дилекторский Б.Д. Экспериментальное исследование
устойчивости цилиндрических оболочек при статическом
внешнем давлении и гидроударе. [Текст] / Б.Д. Дилекторский,
Ю.Г. Коноплев, А.А. Саченков, А.Л. Хавкин // Теория и методы
исследования оболочек сложной формы: Межвузовский сборник.
– Казань: изд-во КАИ. – 1978. – С. 38-42.

9. Карпов В.В. Программный комплекс исследования


устойчивости оболочек. [Текст]: монография / В.В. Карпов,
Д.А. Баранова, Р.Т. Беркалиев. – СПб: СПбГАСУ. – 2009. – 104 с.

10. Коноплев Ю.Г. Теоретико-экспериментальный метод в


задачах устойчивости цилиндрических оболочек эллиптического
сечения. [Текст] / Ю.Г. Коноплев, А.А. Саченков. // Исслед. по
теории пластин и оболочек. – 1984. – Вып. 17, № 1. – С. 135-
152.

11. Коноплев Ю.Г. Устойчивость упругих пластин и оболочек


при нестационарных воздействиях. [Текст]: монография /
Ю.Г. Коноплев, Ф.Х. Тазюков. – Казань: изд-во Казанск. ун-та. –
1994. – 124 с.

12. Минеев В.Е. Исследование устойчивости замкнутых


цилиндрических оболочек при динамическом действии
всестороннего давления. [Текст] / В.Е. Минеев // Исслед. по
теории пластин и оболочек. – 1970. – Вып. 6-7. – C. 596-623.

13. Михлин Б.З. Высокочастотные емкостные и индуктивные


датчики. [Текст]: монография / Б.З. Михлин. – Москва,
Госэнергоиздат. – 1960. – 87 с.

14. Пикуль В.В. К теории устойчивости оболочек. [Текст] /


В.В. Пикуль // Вестник ДВО РАН. – 2006. – № 4. – С. 81-86.

15. Преображенский И.И. Устойчивость и колебания


пластинок и оболочек с отверстиями. [Текст]: монография / И.И.
Преображенский. – Москва, Машиностроение. – 1981. – 192 с.

77
16. Саченков А.В. Об устойчивости кольцевых пластин и
пологой усеченной сферической оболочки при сдвиге. [Текст] /
А.В. Саченков //Исслед. по теории пластин и оболочек. – 1965. –
Вып. 3. – С. 35-42.

17. Саченков А.В. Теоретико-экспериментальный метод


исследования устойчивости пластин и оболочек. [Текст] /
А.В. Саченков //Исслед. по теории пластин и оболочек. – 1970. –
Вып. 6-7. – С. 391-433.

18. Саченков А.В. Динамический критерий устойчивости


пластин и оболочек. [Текст] / А.В. Саченков // Исслед. по теории
пластин и оболочек. – 1976. – Вып. 12. – С. 281-293.

19. Саченков А.В. Решение динамических задач устойчивости


оболочек теоретико-экспериментальным методом. [Текст] /
А.В. Саченков // Прикладные проблемы прочности и
пластичности. – 1977. – № 7. – С. 137-141.

20. Саченков А.А. Устойчивость круговой цилиндрической


оболочки при совместном действии осевой сжимающей силы
и гидравлического удара. [Текст] / А.А. Саченков // Исслед. по
теории пластин и оболочек. – 1984. – Вып. 17, № 2. – С. 65-71.

21. Шалабанов А.К. Способ измерения форм колебаний


твердых деформируемых тел. [Текст]: Авторское свидетельство №
654860. / А.К. Шалабанов. // Бюллетень изобретений. – 1979. – №
12. – С. 36.

22. Якушев В.Л. Нелинейная деформация и устойчивость


тонких оболочек. [Текст]: монография / В.Л. Якушев: РАН. –
Москва: Наука. – 2004. – 274 с.

78
ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ ........................................................................................................................ 3
ПРИНЯТЫЕ ОБОЗНАЧЕНИЯ ..................................................................................... 5
1. ПОСТАНОВКА ЗАДАЧИ........................................................................................ 6
2. МЕТОД РЕШЕНИЯ ЗАДАЧИ ................................................................................ 7
2.1. Вывод разрешающего уравнения ............................................................. 7
2.2. Начальные условия и критерий динамической устойчивости......... 9
2.3. Преобразование разрешающего уравнения ......................................... 10
2.4. Выбор аппроксимирующей функции .................................................... 11
2.5. Решение уравнения ..................................................................................... 11
2.6. Получение расчётных формул ................................................................ 15
3. РЕШЕНИЕ НЕКОТОРЫХ ЗАДАЧ ..................................................................... 19
3.1. Общий алгоритм решения ........................................................................ 19
3.2. Цилиндрическая оболочка при внешнем давлении
(другие случаи закрепления торцов) ................................................................ 19
3.3. Влияние предварительного статического нагружения .................... 21
3.4. Прямоугольная пластина при динамическом сжатии ...................... 23
3.5. Круговая усеченная коническая оболочка
при внешнем давлении ......................................................................................... 28
3.6. Пологая цилиндрическая панель под действием поперечной
нагрузки .................................................................................................................... 36
3.7. Произвольный закон нагружения........................................................... 37
3.8. Круглая пластина при динамическом сжатии .................................... 39
3.9. Кольцевой сектор при равномерном сжатии ..................................... 42
3.10. Кольцевая пластина при динамическом сдвиге .................................. 45
3.11. Подкрепленная цилиндрическая оболочка .......................................... 48
3.12. Оболочка с вырезами .................................................................................. 52

4. МЕТОДИКА ЭКСПЕРИМЕНТА И ЕГО РЕЗУЛЬТАТЫ ............................. 55


4.1. Общие положения ....................................................................................... 55
4.2. Требования к проведению испытаний .................................................. 55
4.3. Форма, размеры и методы отбора образцов ....................................... 56
4.4. Испытательное оборудование ................................................................. 58
4.5. Проведение испытаний и обработка результатов ............................. 62
4.6. Определение констант и функций ......................................................... 73
ЛИТЕРАТУРА................................................................................................................. 76

79
Учебное издание

Коноплев Юрий Геннадьевич


Бахтиева Ляля Узбековна,
Митряйкин Виктор Иванович,
Тазюков Фэрид Хоснутдинович

ДИНАМИЧЕСКАЯ УСТОЙЧИВОСТЬ
УПРУГИХ ПЛАСТИН
И ОБОЛОЧЕК

Учебное пособие

Отпечатано в полном соответствии с представленным оригинал-макетом


_________________________________________________________________

Подписано в печать 24.04.2012.


Форм. 60х84 1/16. Гарнитура «Таймс». Печать ризографическая.
Печ. л. 5. Тирах 150. Заказ 113.

Лаборатория оперативной полиграфии Издательства КФУ


420045, Казань, Кр. Позиция, 2а
Тел. 233-72-12

80