close

Вход

Забыли?

вход по аккаунту

?

Устойчивость и колебания цилиндрической оболочки переменной толщины с криволинейным краем.

код для вставкиСкачать
2003
ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА
Сер. 1
Вып. 3
(№ 17)
МЕХАНИКА
УДК 539.3
Л. С. Елисеева, С. Б. Филиппов
УСТОЙЧИВОСТЬ И КОЛЕБАНИЯ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ
ПЕРЕМЕННОЙ ТОЛЩИНЫ С КРИВОЛИНЕЙНЫМ КРАЕМ∗
1. Введение. Асимптотическим методом решаются задачи о потере устойчивости
под действием равномерного внешнего давления и низкочастотных колебаниях тонких
цилиндрических оболочек переменной толщины. Характерной особенностью рассматриваемых краевых задач является локализация форм потери устойчивости и форм
колебаний вблизи образующих цилиндрических оболочек.
Предполагается, что цилиндрическая оболочка с криволинейным краем замкнута
в окружном направлении. Выберем за единицу длины ее радиус R и введем на срединной поверхности оболочки ортогональную систему безразмерных криволинейных
координат s ∈ [0, l(ϕ)], ϕ ∈ (−π, π], где l(ϕ) — длина образующей. В частности, при
l(ϕ) = lc − tg β cos ϕ оболочка имеет косой край с углом наклона β (рис. 1).
Рис. 1. Цилиндрическая оболочка с косым краем.
Предположим, что безразмерная толщина оболочки h зависит только от координаты
ϕ:
h = h0 g(ϕ),
где h0 — малый параметр, g ∼ 1.
∗ Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант 01-01-00327).
c Л. С. Елисеева, С. Б. Филиппов, 2003
84
Безразмерная система уравнений, описывающая колебания и устойчивость цилиндрической оболочки, имеет вид [1]:
∂2Φ
+ λZ = 0,
∂s2
∂2w
= 0.
ε4 Δ(g −1 ΔΦ) +
∂s2
ε4 Δ(g 3 Δw) −
(1)
Здесь w(s, ϕ) — нормальный прогиб, Φ(s, ϕ) — функция напряжений,
Δ=
∂2
∂2
+
,
2
∂s
∂ϕ2
ε8 =
h20
,
12(1 − ν 2 )
ν — коэффициент Пуассона.
В задаче о потере устойчивости безмоментного напряженного состояния под действием равномерного внешнего давления p
Z = ε2
∂2w
,
∂ϕ2
λ=
p
,
Eh0 ε6
где E — модуль Юнга.
В задаче о свободных колебаниях
Z = −gw,
λ=
ρR2 ω 2
,
ε4 E
где ω — частота колебаний, ρ — плотность материала.
Предполагается, что на краях оболочки s = 0 и s = l(ϕ) задано по четыре однородных граничных условия.
В работах [1], [2] предложен эффективный метод построения асимптотических разложений решений системы (1) для задачи о потере устойчивости. В качестве примеров
в [2] рассмотрены случаи l = lc , h = h(ϕ) и l = l(ϕ), h = h0 . В первом случае формы
потери устойчивости локализуются в окрестности образующей, на которой функция
h(ϕ) имеет строгий минимум. Во втором случае локализация форм происходит вблизи наиболее длинной образующей цилиндрической оболочки. Использование того же
асимптотического метода в задаче о низкочастотных колебаниях описано в статье [1].
В данной работе метод [1], [2] применяется в случае l = l(ϕ), h = h(ϕ), т.е. исследуется влияние на локализацию форм потери устойчивости и форм колебаний одновременного действия двух факторов: криволинейного края и переменной толщины
оболочки.
2. Асимптотические разложения решений. При потере устойчивости под действием внешнего давления и низкочастотных колебаниях напряженное состояние цилиндрической оболочки складывается из основного полубезмоментного состояния и
простого краевого эффекта. Ограничимся определением λ с точностью до величин
O(ε2 ). В этом случае достаточно построить основное состояние, удовлетворив на каждом краю оболочки только двум главным граничным условиям [2].
Решение системы (1) ищем в виде
ϕ
∞
i
w(s, ϕ, ε) = exp
q(ϕ) dϕ
εk wk (s, ϕ),
ε ϕ0
k=0
λ=
∞
εk λk .
(2)
k=0
85
Выполнение условий
Im q(ϕ0 ) = 0,
Im
dq
(ϕ0 ) > 0
dϕ
обеспечивает экспоненциальное затухание решения при удалении от наиболее слабой
образующей ϕ = ϕ0 . Функция Φ имеет такое же асимптотическое разложение как и w.
Подстановка (2) в (1) и граничные условия дает последовательность уравнений для
определения неизвестных функций q(ϕ), wk (s, ϕ), Φk (s, ϕ) и чисел λk . Система уравнений нулевого приближения имеет вид
d2 Φ0
− q 4 g 3 w0 + λ0 N w0 = 0,
ds2
d2 w0
q4
+ Φ0 = 0,
2
ds
g
(3)
где N = q 2 в задаче устойчивости, N = g в случае колебаний.
Исключив из системы (3) функцию Φ0 , получим уравнение
d4 w0
− α4 w0 = 0,
ds4
α4 = λ0
q4 N
− g2 q8 ,
g
(4)
совпадающее по форме с уравнением колебаний балки.
При шарнирном опирании краев оболочки граничные условия для уравнения (4)
имеют вид
d2 w0
w0 =
= 0, s = 0, s = l(ϕ).
ds2
В этом случае
πn
, n = 1, 2, . . .
w0n = sin αn s, αn =
l
В общем случае решение уравнения (4) можно искать в виде линейной комбинации
балочных функций, причем αn (ϕ) = zn /l(ϕ), где числа zn зависят только от граничных
условий (0 < z1 < z2 < . . .).
Из второй формулы (4) следует, что величина λ0 является функцией параметров q
и ϕ:
q4 g3
α4 g
.
(5)
λ0 = f (q, ϕ) = 4 +
q N
N
Следуя [2], в качестве нулевого приближения для собственного значения λ возьмем
λ0 = min f (q, ϕ) = f (q0 , ϕ0 ).
q,ϕ
Числа q0 и ϕ0 являются корнями системы уравнений
λq =
∂f
= 0,
∂q
λϕ =
∂f
= 0.
∂ϕ
(6)
Поправка первого приближения определяется по формуле [2]:
m + 1
λ1 =
λqq λϕϕ − λ2qϕ , m = 0, 1, 2 . . .
2
где частные производные
(7)
∂2f
,
∂q 2
(8)
λqq =
86
λϕϕ =
∂2f
,
∂ϕ2
λqϕ =
∂2f
∂q∂ϕ
вычисляются при q = q0 , ϕ = ϕ0 . Образующая ϕ = ϕ0 является наиболее слабой, если
q0 , ϕ0 — решение системы (6), для которого λqq λϕϕ − λ2qϕ > 0.
3. Закон изменения толщины оболочки. Предположим, что поперечное сечение
оболочки представляет собой круг радиуса R = 1 с центром в точке O, из которого
вырезан круг радиуса r < 1 (рис. 2).
Рис. 2. Поперечное сечение оболочки.
Если расстояние e между центрами кругов не равно нулю, то оболочка имеет переменную безразмерную толщину h(ϕ). Справедливы равенства
a2 + e2 − 2ae cos ϕ = r2 ,
h(ϕ) = 1 − a,
(9)
определяющие точный закон изменения толщины оболочки. Ввиду того, что в данной
работе находятся приближенные решения краевых задач, целесообразно воспользоваться более простым приближенным выражением для функции h(ϕ).
Принимая во внимание равенство (9), получаем, что наименьшее h0 и наибольшее
hm значения толщины оболочки определяются по формулам
h0 = h(0) = 1 − r − e,
hm = h(π) = 1 − r + e.
Следовательно,
r = 1 − (γ + 1)h0 ,
e = γh0 ,
γ = (η − 1)/2,
η = hm /h0 .
(10)
Функцию a ищем в виде
a = a0 + a1 h 0 + . . .
(11)
Подстановка (10) и (11) в (9), дает равенства
a0 = 1,
a1 = γ cos ϕ − γ − 1.
С учетом второй формулы (9) получаем
h(ϕ) " h0 g(ϕ),
g(ϕ) = 1 + γ(1 − cos ϕ).
(12)
4. Устойчивость. Рассмотрим задачу об устойчивости безмоментного напряженного состояния цилиндрической оболочки переменной толщины с криволинейным краем
под действием равномерного внешнего бокового давления. Для определения верхнего критического давления следует найти наименьшее значение параметра нагрузки λ.
Этому значению соответствуют n = 1, m = 0. Формула (5) для задачи устойчивости
принимает вид
α4 g
f (q, ϕ) = 6 + q 2 g 3 .
q
87
Система уравнений (6) эквивалентна системе
q8 =
3α4
,
g2
2l
5g =
,
l
g
(13)
где l = dl/dϕ, g = dg/dϕ. При выводе системы (13) использовано тождество α /α =
−l /l. С учетом формул (8) и (13) получаем
q 2 g 4g 2 α
3
2
2
+ 10gg + 73g
λqq = 16g , λqϕ = 24qg g , λϕϕ =
.
(14)
3
α
Подставив (14) в (7) и воспользовавшись равенством α /α = 2l2 /l2 − l /l, найдем
поправку первого приближения:
15g 2
Λ1
10g 4l
3
, Λ1 =
−
+ 2 .
λ1 = 2qg
3
g
l
g
В качестве примера рассмотрим задачу о потере устойчивости оболочки с косым
краем, толщина которой меняется по закону (12). Второе уравнение системы (13) в
этом случае принимает вид
A(ϕ) sin ϕ = 0,
A(ϕ) = 5γ/g − 2t/l,
(15)
где
g = 1 + γ(1 − cos ϕ),
l = lc − t cos ϕ,
t = tg β,
lc > t,
β — угол наклона косого края (см. рис. 1). Множество корней уравнения (15), лежащих
в интервале (−π, π], состоит из чисел 0, π и корней уравнения A(ϕ) = 0, эквивалентного
уравнению
5γlc − 2t(1 + γ)
cos ϕ = c(γ), c(γ) =
,
3γt
которое при
2t
2t
γ1 < γ < γ2 , γ1 =
, γ2 =
(16)
5lc + t
5lc − 5t
в интервале (−π, π) имеет два решения
ϕ = ϕ∗ ,
ϕ = −ϕ∗ ,
ϕ∗ = arccos c(γ).
(17)
Следовательно, при γ γ1 и γ γ2 уравнение (15) имеет два корня: 0 и π, а при
γ1 < γ < γ2 — четыре корня: 0, π, ϕ∗ , −ϕ∗ . Корню ϕ0 соответствует наиболее слабая
образующая, если
Λ1 (ϕ0 ) = 2A(ϕ0 ) cos ϕ0 + 15γ 2 sin2 ϕ0 /g 2 (ϕ0 ) > 0.
В случае 0 γ < γ1 справедливы неравенства
Λ1 (0) = 2A(0) < 0,
Λ1 (π) = −2A(π) > 0,
поэтому ϕ0 = π, т.е. наиболее слабой является нижняя образующая.
Если выполнено неравенство (16), то
Λ1 (0) < 0,
88
Λ1 (π) < 0,
Λ1 (±ϕ∗ ) = 15γ 2 sin2 ϕ∗ /g 2 (ϕ∗ ) > 0,
и на оболочке имеются две наиболее слабых образующих: ϕ0 = ±ϕ∗ .
При γ > γ2 наиболее слабой образующей соответствует ϕ0 = 0, так как
Λ1 (0) > 0,
Λ1 (π) < 0.
Ввиду того, что λqq λϕϕ − λ2qϕ = 0 при γ = γk , k = 1, 2, полученными асимптотическими формулами можно пользоваться только при γ = γk .
Предположим, что длина оси оболочки lc и угол наклона косого края β фиксированы, а отношение максимальной толщины оболочки к минимальной η = hm /h0 = 1 + 2γ
увеличивается. Как и в случае оболочки постоянной толщины (η = 1), при η < η1 =
1+2γ1 влияние косого края приводит к локализации формы потери устойчивости вблизи нижней, наиболее длинной образующей. Дальнейшее увеличение величины η вызывает появление на оболочке двух наиболее слабых образующих, которые постепенно
поднимаются вверх, до тех пор пока отношение η не достигнет значения η2 = 1 + 2γ2 . В
этот момент две наиболее слабые образующие сливаются друг с другом, превращаясь
в образующую ϕ0 = 0. При последующем увеличении отношения η образующая ϕ0 = 0
остается наиболее слабой. Это свидетельствует о том, что при η > η2 переменность
толщины оказывает определяющее воздействие на положение наиболее слабой образующей, так как образующая ϕ0 = 0 является наиболее слабой и для оболочки с прямыми
краями (β = 0).
В таблице 1 для шарнирной опертой оболочки с параметрами lc = 3, h0 = 0.01,
β = π/4, ν = 0.3 приведены значения q0 , ϕ0 , λ0 , λ1 и λ = λ0 + ελ1 при различных
значениях η. В рассматриваемом случае η1 = 1.25, η2 = 1.4.
Таблица 1
η
q0
ϕ0
λ0
λ1
λ
1.00
1.10
1.20
1.30
1.35
1.45
1.50
1.02
0.99
0.97
1.11
1.27
1.44
1.44
3.14
3.14
3.14
1.68
1.02
0.00
0.00
1.34
1.74
2.17
2.60
2.72
2.76
2.76
1.17
1.13
0.79
1.01
0.99
0.83
1.17
1.65
2.01
2.36
2.84
2.95
2.95
3.03
Параметр нагрузки λ увеличивается с увеличением η. При η > η2 параметры q0 и
λ0 не зависят от η, и увеличение λ происходит за счет роста поправки λ1 .
5. Колебания. Рассмотрим свободные колебания цилиндрической оболочки переменной толщины с криволинейным краем. В этом случае формула (5) приобретает вид
f (q, ϕ) =
α4
+ q4 g2 ,
q4
а система уравнений (6) эквивалентна системе
q4 =
α2
,
g
2l
g
= .
l
g
Из формул (8) и (18) следует, что
λqq = 32g 2 q 2 ,
λqϕ = 16q 3 gg ,
λϕϕ = q 4
(18)
4g 2 α
+ 2gg + 5g 2 .
α
(19)
89
Поправку первого приближения находим с помощью формул (7) и (19):
λ1 = 2(m + 1)q 3 g 2
2Λ1 ,
Λ1 =
4l
g 2
2g −
− 2.
g
l
g
Определим положения наиболее слабых образующих на цилиндрической оболочке
с косым краем, толщина которой меняется по закону (12). Значения угла ϕ, которым
могут соответствовать слабые образующие, являются корнями второго уравнения системы (18). В рассматриваемом случае это уравнение принимает вид
A(ϕ) sin ϕ = 0,
A(ϕ) = γ/g − 2t/l.
(20)
Пусть ϕ0 ∈ (−π, π) — корень уравнения A(ϕ) = 0. Тогда
Λ1 (ϕ0 ) = 2A(ϕ0 ) cos ϕ0 − γ 2 sin2 ϕ0 /g 2 (ϕ0 ) = −γ 2 sin2 ϕ0 /g 2 (ϕ0 ) < 0.
Следовательно, образующая ϕ = ϕ0 не является наиболее слабой. Наиболее слабыми
могут быть только образующие ϕ = 0 и ϕ = π, причем первая из них будет наиболее
слабой при
(21)
Λ1 (0) = 2A(0) > 0,
а вторая — при
Λ1 (π) = −2A(π) > 0.
Введем обозначения:
γ1 =
2t
,
lc − t
γ2 =
(22)
2t
.
3lc − t
Неравенство (21) справедливо при γ > γ1 , а неравенство (22) при 0 γ < γ2 в случае
lc > 3t и при любых γ 0 в случае lc 3t.
Зафиксировав длину оси оболочки lc и угол наклона косого края β, будем увеличивать отношение η максимальной толщины оболочки к минимальной. При 1 η <
η1 = 1 + 2γ1 формы колебаний локализованы вблизи наиболее длинной образующей
ϕ = π. Дальнейшее увеличение η приводит к появлению на оболочке второй наиболее
слабой образующей ϕ = 0. Если lc > 3 tg β, то при достижении η величины η2 = 1 + 2γ2
образующая ϕ = π перестает быть наиболее слабой, и при η > η2 на оболочке имеется
только одна наиболее слабая образующая ϕ = 0. В противном случае при всех η > η1
на оболочке имеются две наиболее слабые образующие ϕ = π и ϕ = 0.
Алгоритмы приближенного решения краевых задач теории колебаний и теории
устойчивости практически совпадают, однако решения этих задач при l(ϕ) = lc −t cos ϕ,
h(ϕ) = h0 [1+γ(1−cos ϕ)] различаются не только количественно, но и качественно. В задаче устойчивости при изменении параметра γ в интервале (γ1 , γ2 ) две наиболее слабые
образующие ϕ = ϕ0 и ϕ = −ϕ0 перемещаются по поверхности оболочки (0 < ϕ0 < π). В
задаче колебаний тоже имеется область изменения γ, которой соответствуют две наиболее слабые образующие ϕ = 0 и ϕ = π, однако положение этих образующих не зависит
от γ.
Найдем приближенные значения наименьшего параметра частоты λ для шарнирно
опертой оболочки с параметрами lc = 3, h0 = 0.01, β = π/4, ν = 0.3. В таблице 2
приведены значения q0 , ϕ0 , λ0 , λ1 и λ = λ0 + ελ1 при различных значениях η.
В рассматриваемом примере η1 = 3. Из тех значений λ, которые при η1 > 3 соответствуют двум наиболее слабым образующим ϕ = 0 и ϕ = π, в таблице 2 приведено
90
Таблица 2
η
q0
ϕ0
λ0
λ1
λ
1.00
1.50
2.00
2.50
3.50
4.00
4.50
0.89
0.80
0.74
0.70
0.65
0.63
1.25
3.14
3.14
3.14
3.14
3.14
3.14
0.00
1.23
1.85
2.47
3.08
4.32
4.93
4.93
1.97
2.67
3.31
3.91
5.04
5.57
6.82
1.69
2.48
3.24
4.00
5.50
6.24
6.53
наименьшее. При η = 4 образующим ϕ = 0 и ϕ = π соответствуют одинаковые наименьшие значения λ0 и λ1 , однако параметры q0 , характеризующие скорость затухания
форм колебаний, различаются. В таблице приведено значение q0 = 0.63 для образующей
ϕ = π. Для образующей ϕ = 0 эта величина в два раза больше, поэтому в окрестности
ϕ = 0 форма колебаний затухает быстрее, чем вблизи ϕ = π.
Summary
Eliseeva L. S., Filippov S. B. Buckling and vibrations of a cylindrical shell of variable thickness
with a curvilinear edge.
Buckling under the uniform external pressure and low-frequency vibrations of a thin circular
closed cylindrical shell are considered. The shell has variable thickness and a curvilinear edge.
Asypmtotic expansions of the localized buckling and vibration modes are derived and analуzed.
Approximate explicit formulae for the critical external pressure and the fundamental vibration
frequency are obtained.
Литература
1. Товстик П. Е. Некоторые задачи устойчивости цилиндрических и конических оболочек
// Прикл. математика и механика. 1983. Т. 47. Вып. 5. С. 815–822.
2. Товстик П. Е. Устойчивость тонких оболочек: асимптотические методы. М., 1995.
Статья поступила в редакцию 26 декабря 2002 г.
91
Документ
Категория
Без категории
Просмотров
7
Размер файла
235 Кб
Теги
крае, толщины, оболочка, цилиндрическом, устойчивость, колебания, переменных, криволинейных
1/--страниц
Пожаловаться на содержимое документа