(Ю. Г. Коноплев и др.) - Динамическая устойчивость упругих пластин и оболочек - учебное пособие-Казанский ун-т (2012)
(Ю. Г. Коноплев и др.) - Динамическая устойчивость упругих пластин и оболочек - учебное пособие-Казанский ун-т (2012)
(Ю. Г. Коноплев и др.) - Динамическая устойчивость упругих пластин и оболочек - учебное пособие-Казанский ун-т (2012)
УПРУГИХ ПЛАСТИН И
ОБОЛОЧЕК
Учебное пособие
Казанский университет
2012
УДК 539.3
Д46
Авторы:
Ю.Г. Коноплев, Л.У. Бахтиева, В.И. Митряйкин,
Ф.Х. Тазюков
Рецензенты:
доктор физ.-мат. наук Ю.П. Артюхин (К(П)ФУ)
доктор физ.-мат. наук В.Н. Паймушин (КНИТУ-КАИ)
ISBN 978-5-905787-30-0
3
Учебное пособие состоит из четырех частей.
В первой части дается общая постановка задачи, предлагается
алгоритм ее решения.
Во второй части учебного пособия исследуется
линеаризованное уравнение динамики оболочки, формулируются
начальные условия и критерий динамической устойчивости. На
примере круговой цилиндрической оболочки, находящейся под
действием внешнего давления, линейно изменяющегося во
времени, описывается метод решения задачи, позволяющий
получить формулы для расчета критических нагрузок, волновых
чисел, коэффициентов динамичности и критических импульсов.
Третья часть посвящена решению ряда конкретных задач.
Рассмотрены тонкие изотропные оболочки и пластины
постоянной толщины, различным образом закрепленные по
торцам и подверженные нагрузкам, изменяющимся во времени по
некоторым законам (линейному, квадратичному и др.). Для всех
задач получены расчетные формулы, позволяющие делать выводы
о влиянии различных физических и геометрических параметров на
динамическую устойчивость. Для уточнения констант, входящих
в полученные формулы, рекомендовано проведение
экспериментальных исследований.
В четвертой части учебного пособия описана методика
проведения испытаний оболочек и пластин на динамическую
устойчивость. Приведены результаты некоторых экспериментов,
сделаны выводы о практической применимости полученных
аналитических зависимостей.
Предложенные рекомендации могут быть использованы при
расчёте на устойчивость элементов конструкций судов,
летательных аппаратов, различных емкостей, резервуаров и
других инженерных сооружений, находящихся под действием
быстро изменяющихся во времени нагрузок.
Учебное пособие предназначено для студентов и аспирантов,
специализирующихся по методам расчета конструкций, а также
для работников научно-исследовательских институтов и
конструкторских бюро.
4
ПРИНЯТЫЕ ОБОЗНАЧЕНИЯ
х, у, z – декартовы координаты,
t – время,
h – толщина оболочки (пластины),
Е – модуль упругости,
ν – коэффициент Пуассона,
ρ – плотность материала оболочки,
E
V = – скорость звука в материале, ,
ρ
Е h3
D=
(
12 1 − ν 2 ) – изгибная жёсткость,
∂ 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. ПОСТАНОВКА ЗАДАЧИ
6
f) полученные зависимости дополнить требованием, чтобы в
предельном случае (при медленном нагружении) они совпадали с
соответствующими формулами для статических задач.
4. Рассчитать критические нагрузки, волновые числа,
коэффициенты динамичности, критические импульсы.
Результаты получить в виде формул, справедливых при
произвольных параметрах оболочек и пластин.
5. Уточнить экспериментальным путем коэффициенты,
входящие в полученные структурные зависимости.
⎡⎛ ∂ 2 w ⎞ 2 ∂ 2 w ∂ 2 w ⎤ Еh ∂ 2 w
∇ 4 F = Еh ⎢⎜⎜ ⎟⎟ − 2 ⎥−
⎢⎣⎝ ∂x ∂у ⎠ ∂x ∂у 2 ⎥ R ∂x 2 ,
⎦
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
F = f (t ) sin w = ξ (t ) sin
m1 x ny m1 x ny
sin , sin , (5)
R R R R
mπR
m1 =
l
8
d 2ξ
+ λ (t ) ξ = 0 ,
dt 2
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 .
9
ξʹ (1) = 0, (8)
W = W0 + ξ , (9)
W0
Z ′′ + λ (τ ) Z = λ (τ ) , Z = W/h.
h
EhW0 T20 P0 R 2t
= = − P0 t , откуда W0 = − , (10)
R2 R Eh
10
Задача, таким образом, сводится к решению уравнения (11) с
начальными условиями (12) и дополнительным условием (13). Для
решения используется процедура метода Бубнова.
(
Z = A τ 2 + Bτ α , )
где А – амплитуда прогиба, α – константа, определяемая в
процессе решения задачи ( α ≠ 2 ). Из условия (13) следует, что
B = – 2/α , тогда
⎛ 2 ⎞
Z = A ⎜τ 2 − τ α ⎟ . (14)
⎝ α ⎠
1
⎛ 1 ⎞ α −2
τ1 = ⎜ ⎟ .
⎝ α −1 ⎠
(15)
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)
⎡1 2 ⎤ 5 ⎡1 2 ⎤
с4 = τ 14 ⎢ − ⎥, с5 = τ 1 ⎢ − ⎥.
⎣ 4 α (α − 1) (α + 2) ⎦ ⎣ 5 α (α − 1) (α + 3) ⎦
А с2 + Р 0 с4 = А с3 + Р 0 с5.
12
D 2 β 2 m14 b1ρ h R
Р0 tк = (n + 6 )+ 2 2 , (18)
R3 n n tк
Таблица 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
где
− 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)
b2 D 2 β 2 m14 b3 ρ h R
Р0 tк = (n + 6 )+ 2 2 , (20)
R3 n n tк
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)
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Па
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
из которых находим
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
tк tк
P0 t к2
I к = ∫ Р(t )dt = ∫ Р0 t dt =
0 0
2 .
c2 h2 E
Iк = ⋅
2 4 1 −ν 2 RV . (27)
18
3. РЕШЕНИЕ НЕКОТОРЫХ ЗАДАЧ
кm1 x ny кm1 x ny
F = f (t ) sin sin , w = ξ (t ) sin sin ,
R R R R
19
где к – коэффициент, характеризующий способ закрепления
торцов оболочки.
Таблица 3
Результаты расчетов для различных видов закрепления
Шарнирное опирание
1 240 240 490 670
20
динамической критической нагрузки при высоких скоростях
нагружения не зависит от характера закрепления торцов оболочки,
т. к. величина m1 в эту формулу не входит. Коэффициент γ в
таблице 3 равен отношению критического давления при заданном
виде закрепления к критическому давлению при шарнирном
опирании.
⎡ β 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
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].
∂2w
D∇ 4 w − T10wxx − T20wyy − 2S10wyy + ρ h = 0,
∂t 2
mπх nπу
w = ξ (t ) sin sin
a b .
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к
D bρ b
q0 tк = (m12 b 2 + 2) + 12 2 , b= .
hb 2
m1 tк π (28)
2D 2 2
q0 tк = (m1 b + 1) . (29)
hb 2
π2
2
4D ⎛h⎞
q = 2 = q* E ⎜ ⎟ ,
к
q* =
3 (1 −ν 2 ) , m = 1.
ст (30)
b h ⎝b⎠
b1
K д4 − K д3 − = 0, (31)
S
25
2
⎛ EVπ h3 ⎞
где K д = q0tk к
q , S = q ⎜⎜ *3
⎟
4 ⎟ .
ст
⎝ q0 b ⎠
2
⎛a⎞
m = (2 K д − 1) ⎜ ⎟ .
2
(32)
⎝b⎠
26
mπx 2 πy
w = ξ (t ) sin sin
a b .
λ (τ ) =
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)
2D ⎛ 2 2 4 ⎞
q0 t к = ⎜ т1 b + ⎟ . (36)
b 2h ⎝ 3⎠
4
16 3 128 2 b1 ⎛ 4 ⎞
K д4 − Kд − Kд − ⎜ ⎟ = 0 , (37)
21 441 S ⎝7⎠
27
которое решается численно.
Таблица 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
mπx
w = ξ (t ) sin y.
a
28
3.5. Круговая усеченная коническая оболочка
при внешнем давлении
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
Р(t ) r tgγ
Т1 = − , T20 = − P(t ) r tgγ .
2
29
r0
γ
R0
r
r1
φ
R1
Рисунок 3
30
μ – параметр, при выборе которого следует пользоваться
рекомендациями работы [1], а именно:
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
31
⎛ 2 η12 ⋅ m24 ⎞ b ρ hr
Р ( tк ) =
D
⎜⎜ n1 θ + ⎟⎟ + 12 2 1
r1 θ tgγ
3
⎝ n1 ⎠ θ n1 t к tgγ .
6 (40)
⎛ 2 η 2 m4 ⎞
Р ( tк ) =
2D
⎜⎜ n1 θ − 1 6 2 ⎟⎟ (41)
r1 θ tgγ
3
⎝ n ⎠
(1 − 3а)⋅ (1 − а) 2
= Р0 а (42)
32
2
1 P 0 η12 m24 b1 ρ h r110 Р02 θ tg 2γ
а= , откуда n18 = = . (43)
Р0 θ 4 D 3θ 3
1 3
4 Dη1 2 m2 θ 4
Рст = ⋅ ,
r13θ tgγ
3
4
3
Таблица 5
Зависимость параметра μ от z
θ ст θ η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)
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
n13 = P1 θ . (48)
34
3 1
Р t3 c 2 E h2 θ 2
Iк = 1 к = 3 . (50)
3 3 1 − ν 2 V r1 θ tgγ
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γ
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)
θ
d 2ξ ⎛ P(t ) ⎞
+ ω 2 ⎜⎜1 + ⎟ ξ = 0,
dt 2
⎝ Pст ⎟⎠ (57)
d 2ξ ⎛ P(t ) ⎞
+ ⎜ 1 + ⎟ ξ = 0.
dτ 2 ⎜⎝ Рст ⎟⎠ (58)
36
Анализ уравнения (58) показывает, что структурная формула
для коэффициента динамичности устанавливает связь между
динамической неустойчивостью, колебаниями и устойчивостью
при статических нагрузках. При этом оказывается, что изменение
формы, размеров оболочки, условий закрепления и т. п. влияет на
процесс динамической потери устойчивости лишь постольку,
поскольку при этом изменяется спектр собственных частот и
критические нагрузки в статике. Иначе говоря, влияние
конкретных свойств геометрии оболочки на её устойчивость при
динамических нагрузках в достаточной мере учитывается такими
её интегральными характеристиками, как ω и Рст. Поэтому
формула для коэффициента динамичности может быть записана в
виде
37
Р(t)
P(t)
J к∧
I к∧
0 tк t
Рисунок 4
q0 t
Рисунок 5
38
Для определения критического времени tк достаточно найти
∧
величину критического импульса I к для линейного закона
нагружения, построить график функции Р(t) и подобрать
графически величину tк таким образом, чтобы площадь
∧
криволинейной трапеции оказалась равной I к .
Tr 0 = Tφ 0 = −q0 t h , Trφ = 0 ,
∂2w
DΔΔw + q0t Δw + ρh = 0, (60)
∂t 2
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,
D b1 ρ R 2
q0 tк = α n2,m + .
R 2 h α n2,m tк2 (62)
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
⎡ r ⎛ r ⎞n ⎤
wmn = J n ( β n, m ) ⎢ − ⎜ ⎟ ⎥,
⎣⎢ R ⎝ R ⎠ ⎦⎥
r J n (r ) − (1 −ν ) J n +1 (r ) = 0,
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 ⎠
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 ,
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
к
( )
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
,
44
таком угле возможен, так как при достаточно высоких скоростях
нагружения справедлива формула [11]
⎢ ln ⎜⎝ R2 ⎟⎠ (R1 + R2 )⎥
1/ 4
1/ 2
⎡ ⎛ R1 ⎞ 2 2 ⎤
⎛ q Eh ⎞
q0 tк = 3,96 ⎜ 0 ⎟ ⎢ ⎥ .
⎝ V ⎠ ⎢ R12 − R22 ⎥
⎣ ⎦
2 Trφ 2 ∂ ⎛ 1 ∂w ⎞ ∂2w
DΔ Δw − 2 r0 ⎜ ⎟ + ρ h 2 = 0,
r ∂r ⎜⎝ r ∂φ ⎟⎠ ∂t
πz
w = ξ (t ) sin sin n (φ − kz),
z1
45
Интегрирование по Бубнову приводит к следующему
характеристическому уравнению [11]
( )
+ α −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
α 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
4/3 n 9 11 14 19
2 n 4 5 6 8
4 n 2 2 3 4
⎛ d 6V d 4V d 2V ⎞
( EI )ш ⎜⎜ 6 + 2 4 + 2 ⎟⎟ = − R 4 T * ,
⎝ dφ dφ dφ ⎠
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
Решение уравнения следует искать в виде
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
51
3.12. Оболочка с вырезами
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 (φ )ζ (τ ).
∂φ
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 ⎦
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φ.
⎠
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
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 ,
⎛ Рфэ ⎞ ⎛ I фэ ⎞
⎜ − 1⎟ ⎜ − 1⎟ = 1
⎜Р ⎟⎜ ⎟ ,
⎝ ст ⎠ ⎝ I к ⎠
э
с асимптотами Рст и Iк, где Рф - давление во фронте волны, I фэ −
импульс давления. В качестве Рст может быть принято
критическое значение при статическом внешнем давлении,
полученное экспериментально, или близкая к нему теоретическая
величина [17].
Построенная на основе экспериментальных данных кривая
устойчивости позволяет определить асимптотическое значение Iк,
т. е. величину критического импульса давления при высоких
скоростях нагружения. Эта величина сравнивается со значением,
рассчитанным по формуле (55).
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 мм. Модуль
упругости определялся динамическим способом.
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
синхронизации, ДД – датчик давления, ДП – датчик перемещения,
К – коммутатор, ПД – преобразователь формы импульса, У –
усилитель, ФА – фотокамера, ∑ - сумматор аналоговых сигналов.
До начала испытаний проводилась тарировка всего
измерительного комплекса по всем каналам и со всеми
коэффициентами усиления. Датчик давления градуировался на
специально изготовленной тарировочной трубе.
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
Рисунок 8
2γ=600
I фэ КПа·с z1=1,260
Рст.=0,020 МПа
0,08
0,04
Рст
Iк = 0,0310 КПа·с
Рисунок 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
ром приведены осциллограммы процесса потери устойчивости и
построенные по ним эпюры прогибов для различных моментов
времени: при квазистатическом нагружении (а), и для
ступенчатого импульса (б). При более интенсивном нагружении
указанный эффект исчезал и вмятина начинала развиваться
симметрично, захватывая сразу почти всю поверхность панели,
так что изгибные деформации концентрировались возле
прямолинейных кромок.
Вид остаточных вмятин после испытания панели на
динамическую нагрузку может ввести в заблуждение
относительно формы потери устойчивости. Как правило, после
интенсивного нагружения у всех испытанных панелей более или
менее отчётливо наблюдалось по две вмятины. Однако
осциллограммы прогибов свидетельствуют об образовании в
процессе потери устойчивости одной обширной вмятины с весьма
пологим центром и крутыми склонами. Причём это наблюдалось
даже у тех панелей, которые при статическом нагружении теряли
устойчивость с образованием двух вмятин.
Очевидно, иллюзия образования двух вмятин после
интенсивного нагружения связана с уже отмеченной
концентрацией изгибных напряжений, а, следовательно, и
остаточных деформаций вблизи прямолинейных кромок.
73
полученного по зависимости (55) можно получить значение
коэффициента c0. Как показал эксперимент, с учетом коррекции
это значение необходимо брать в пределах 6,03 < c0 < 6,29.
Для цилиндрических панелей на основе экспериментальных
данных построены формулы для коэффициента динамичности K
при линейном и квадратичном законах нагружения, а также для
критической длительности импульса tк при ступенчатом законе
нагружения.
При расчёте панелей на динамическую устойчивость
необходимо, прежде всего, определить – расчётом, по справочным
данным или экспериментально – низшую собственную частоту
панели ω и критическое давление при статическом нагружении
Рст. Затем, задавшись величиной скорости нарастания нагрузки q0,
определить соответствующее значение параметра α и по
приведенным на рисунках 17-19 эмпирическим формулам
вычислить значение коэффициента динамичности.
Величина динамической критической нагрузки определяется
как
Р ин = K Рст .
−0 , 36
1,1 ⎛ P ⎞
tk = ⎜⎜ ⎟⎟
ω ⎝ Pст ⎠ .
74
Рисунок 20 ( а) Рисунок 20( б)
а:Δt = 0,5 мс б: Δt = 0,25 мс
Рисунок 20
75
ЛИТЕРАТУРА
76
8. Дилекторский Б.Д. Экспериментальное исследование
устойчивости цилиндрических оболочек при статическом
внешнем давлении и гидроударе. [Текст] / Б.Д. Дилекторский,
Ю.Г. Коноплев, А.А. Саченков, А.Л. Хавкин // Теория и методы
исследования оболочек сложной формы: Межвузовский сборник.
– Казань: изд-во КАИ. – 1978. – С. 38-42.
77
16. Саченков А.В. Об устойчивости кольцевых пластин и
пологой усеченной сферической оболочки при сдвиге. [Текст] /
А.В. Саченков //Исслед. по теории пластин и оболочек. – 1965. –
Вып. 3. – С. 35-42.
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
79
Учебное издание
ДИНАМИЧЕСКАЯ УСТОЙЧИВОСТЬ
УПРУГИХ ПЛАСТИН
И ОБОЛОЧЕК
Учебное пособие
80