Следует отметить, что при любых граничных концевых условиях, структура формулы (20) сохраняется; изменяется лишь содержание коэффициента A. Так, например, для жёстко заделанной по концам оболочки граничные условия имеют вид
2 |
|
x L 0 , |
|
|
x L 0 , |
(21) |
|
|
что равносильно отсутствию радиальных прогибов и углов поворота концевых сечений. Для оболочки со свободными краями концевые граничные условия соответствуют ра-
венствам
2 |
|
x L 0 , |
2 |
x L 0 , |
(22) |
|
то есть нормальные и касательные усилия должны быть равны нулю.
Достоверность полученных результатов подтверждается предельным переходом к оболочкам постоянной толщины известным из литературы [1]:
|
q |
kp |
EI |
n2 1 |
|
8ER3hK 3 |
A, K |
4 |
|
Dn4 (n2 1)2 |
, |
|
(23) |
|||||||||||||
|
|
R3 |
|
n4 (n2 1) |
|
|
4ER6h |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
где для оболочки со свободно опёртыми, заделанными и свободными краями |
|
|
||||||||||||||||||||||||
A |
2(sh2 KL cos2 KL) |
, |
|
A |
|
|
shKL |
sin 2KL |
|
|
, |
|
A |
shKL sin 2KL |
. |
|||||||||||
sh2KL sin 2KL |
|
|
|
2(sh2 KL sin2 KL) |
|
2(ch2 KL cos2 |
KL) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Для длинной оболочки n=2, A=1. Тогда из (20) имеем формулу, |
|
|
|
|||||||||||||||||||||||
|
|
|
|
q |
kp |
|
EI |
|
3 |
|
|
0,38 |
|
Eh2 |
|
h |
, |
|
|
|
|
|||||
|
|
|
|
|
R3 |
(1 2 )0,75 |
R |
|
R |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
полученную Л.Н. Балабухом и В.М. Марченко [3] на основе интегрального уравнения Фредгольма с симметричным ядром (EI=0)
q |
|
12 |
|
Eh |
h 1,5 |
||
|
|
|
|
|
. |
||
9 |
|
(1 2 )0,75 |
|
||||
kp |
|
|
|
R |
|||
|
Численные расчёты, выполненные на ЭВМ, показывают, что коэффициент A в значи- |
|||
тельной |
степени |
зависит от безразмерного параметра |
K2 L , относительной толщины |
|
h h1 h2 |
и длины |
L усиленного пояса. Для оболочки с жёстко защемлёнными и сво- |
||
бодными краями коэффициент A →1 при K2 L →∞. В то время, когда K2 L →0, величина ко- |
||||
эффициента A →∞. Это факт свидетельствует о том, что концевые диафрагмы сближаются с |
||||
уменьшением длины оболочки и дляK2 L 0 совпадают с силовым шпангоутом, поэтому |
||||
q |
. |
|
|
|
kp |
Для оболочки со свободными краями коэффициент A →0 при K2 L →0, так как в этом |
|||
случае концевые диафрагмы отсутствуют. Выявлена также зависимость A=A( h ) при фикси- |
||||
рованных значениях безразмерных параметров =0,2;0,6 |
и K2L 0,6;0,9;1,8 Безразмерная |
|||
толщина изменялась в пределах 0≤h ≤2,05. Анализ графических зависимостей показывает, |
||||
что для защемлённых и шарнирно-опёртых оболочек при прочных равных условиях, коэффициент A →1 с увеличением h . Для оболочек со свободными краями при отсутствии диафрагм 0<A≤1. Графики A=A( ) и K2L 0,6;0,9;1,8 при h =1,9 показывают рост A с увеличением параметра и для случая жёсткой заделки A= 4,25 при =0,5 и K2 L 0,6 .
Выводы. В работе предложен подход к решению задач устойчивости цилиндрических конструктивно-ортотропных оболочек, опирающийся на известный энергетический метод, представляющий собой синтез сопротивления материалов, строительной механики и вариационного исчисления. Решённая этим методом задача, касается одного из важнейших вопросов – устойчивости оболочек переменной жёсткости. Последние, являясь следствием тенденции к снижению массы, рациональному распределению материала и повышению эксплуатационных качеств конструкции, исследованы недостаточно. Есть основания полагать, что задачи термоупругости, колебаний и напряженно-деформированного состояния этих оболочек, могут быть решены рассмотренным здесь подходом. Достоверность полученных ре-
85
зультатов опирается на предельный переход к формулам, известным из литературы и сравнительный анализ численных расчётов.
Библиографический список
1.Кан С.Н. Строительная механика оболочек. - М.: Машиностроение, 1966. -508с.
2.Булатов С.Н. Устойчивость ступенчатой цилиндрической оболочки под действием радиальных давлений и осевых усилий. – М.: Строительная механика и расчет сооружений.
–М., 1970, №3.
3.Балабух Л.Н., Марченко В.М. Устойчивость тонкой оболочки, нагруженной местными усилиями. Труды ЦАГИ, 1952.
References
1.Caen S.N. Building mechanics of shells - Moscow: Machinery construction (Mashinostroenie), 1966. - 508p.
2.Bulatov S.N. Stability of stepped barrel shell under radial pressure and axial thrusts. - Moscow: Building mechanics and analysis of structure, 1970, №3
3.Balabukh L.N., Marchenko V.M. Stability of thin shell loaded with local load. - Researches of TsAGI, 1952
УДК 624.071.2 |
|
Воронежский государственный |
Voronezh State University of Architecture |
архитектурно-строительный университет |
and Civil Engineering |
Д-р техн. наук, проф. С.Н. Булатов |
Dr. tehn. sciences professor S.N. Bulatov |
Россия, г. Воронеж, тел. 8(4732)71-52-02 |
Russia, Voronezh, ph. 8(4732)71-52-02 |
С.Н. Булатов
ОБ ОДНОМ ПОДХОДЕ К ИССЛЕДОВАНИЮ ДЛИННЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК ПЕРЕМЕННОЙ ЖЁСТКОСТИ
Исследована устойчивость длинной конструктивно ортотропной цилиндрической оболочки переменной толщины при локальном нагружении боковым давлением. Получены достаточно простые расчетные зависимости для определения несущей способности тонкостенной конструкции.
Ключевые слова: цилиндрические оболочки, боковое давление, несущая способность.
S.N. Bulatov
ABOUT ONE GOING NEAR RESEARCH OF CYLINDRICAL SHELLS
OF VARIABLE INFLEXIBILITY
Stability long structurally orthotropic a cylindrical cover of a variable thickness is investigated at local loading by lateral pressure. Simple enough settlement dependences for definition of bearing ability of a thin-walled design are received.
Keywords: barrel shells, side pressure, load-carrying ability.
В работе предложен один из возможных подходов к определению критической радиальной нагрузки, приложенной к боковой поверхности конструктивно-ортотропной оболоч-
© Булатов С.Н., 2009
86
ки переменной толщины. Оболочка предполагается достаточно длинной, геометрически совершенной. Получены расчётные формулы и их частные случаи, известные из литературы. Последнее служит достаточным основанием для их практического применения в стройиндустрии, общем и химическом машиностроении, авиационной и ракетно-космической технике.
На рис. 1 приведена расчётная схема и геометрические размеры оболочки, нагруженной равномерно распределённым по поясу толщины h1 и длинны 2 , радиальным давлением интенсивности P. При P=Pкр конструкция разрушается – теряет устойчивость с образованием регулярно расположенных по окружности волн. Задача по определению Pкр решается энергетическим методом.
Рис. 1. Расчетная схема оболочки, нагруженной равномерно распределенным по поясу радиальным давлением
Следуя [1], потенциальную энергию деформации длинной оболочки запишем в виде
2 1 |
dx 2 2 dx 0, |
(1) |
0 |
|
|
где потенциальная энергия деформации первой Г1 и Г2 второй секций на единицу длины со-
ответственно равна |
D1 n2 |
1 2 |
|
2 |
|
ER2 h1 |
|
" 2 |
|
n2 |
1 |
|
2 |
|
|
||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
PKP |
|
R, |
(2) |
2R4 |
|
2n4 |
2R |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
ER h2 |
|
|
|
|
|
|
|||
|
|
|
D1 n |
|
1 |
2 |
|
|
" 2 R. |
|
|
(3) |
|||||||
|
|
|
|
|
|
2n4 |
|
|
|||||||||||
|
|
2 |
|
2R4 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь и ниже обозначения приняты по [2]. Функциональная неизвестная (x) не
терпит разрыв в силу симметричного расположения h1 относительно срединной поверхности оболочки (рис.1). В плоскости симметрии хоу предполагается отсутствие силового шпангоута (EI=0) и депланации контура поперечного сечения. Так что угол поворота при х=0
d |
|
x 0 |
0. |
(4) |
|
||||
dx |
|
|||
|
|
|
||
Условию (4) удовлетворяет функциональная неизвестная |
|
|||
fe kx sin kx cos kx , |
(5) |
|||
где f - амплитудное значение функции ,k - |
неизвестный коэффициент затухания. |
Про- |
||
дольные U, окружные V и радиальные W перемещения на краях длинной оболочки равны нулю. Поэтому Pкр не зависит от краевых граничных условий. С учётом (2) ÷ (5) выражение потенциальной энергии деформации (1) получает вид
6 1 t |
4 C 1 t |
2 |
|
2 |
1 |
P |
0 , |
(6) |
|||||||||||||||
|
|
|
|
1 |
1 |
|
3 |
|
|
|
|
|
|
|
|
1 |
êð |
|
|
|
|||
откуда искомая критическая нагрузка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Р |
|
|
|
в 1 t |
|
|
|
|
C 1 t |
|
|
|
|
1 |
|
, |
(7) |
||||||
кр |
|
1 |
1 |
|
2 |
2 |
|
1 |
|||||||||||||||
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
87
|
D1 n2 1 |
|
|
ER2 h1 K 4 |
|
|
|
|
1 |
|
|
|
1 |
2k |
2 |
cos 2k sin 2k , |
||
|
|
|
|
|
|
|
|
|
|
|||||||||
где в |
R3 |
, |
C |
|
, |
t1 1, |
t2 |
|
|
1 |
, 1 |
|
3 e |
|
||||
n4 n2 1 |
|
|
|
|||||||||||||||
|
|
|
||||||||||||||||
2 e 2k cos 2k sin 2k 2 .
Для изотропной оболочки переменной толщины t1 t2 t
P |
(1 t ) 4 C |
(1 t |
2 |
) |
(1 ) . |
(8) |
|||||||||||||||
kp |
|
|
|
1 |
|
3 |
|
|
|
|
|
|
|
|
1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Отсюда, как частные случаи, вытекают формулы для изотропной оболочки постоян- |
|||||||||||||||||||||
ной толщины |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D(n2 1 |
|
|
4 ER3hK |
4 |
|
|||||||||||||||
|
Pkp |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
(1 1) . |
(9) |
|||
|
|
R |
|
|
3 n |
4 |
(n |
2 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
1) |
|
|||||||||
С увеличением ширины нагружаемого пояса ( L) параметр |
k , а 1 0 , и |
||||||||||||||||||||
формула (9) принимает вид |
|
D n2 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Pкр |
|
|
4 ER3hK 4 |
|
||||||||||||||||
|
|
|
R3 |
|
|
|
|
3 |
n4 |
n2 |
1 |
. |
(10) |
||||||||
Для длинных оболочек n=2, поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
P |
|
|
3D |
|
ER3hK 4 |
|
|
|
(11) |
||||||||||
|
|
|
|
R3 |
|
|
|
36 |
|
|
|
|
|
||||||||
|
|
кр |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Минимальное значение Pкр имеет место при k=0. Тогда из (11) получим известную |
|||||||||||||||||||||
формулу Брайана [4] для нагруженной оболочки по всей поверхности |
|
||||||||||||||||||||
|
|
Pкр |
|
|
|
|
Eh3 |
|
|
|
|
|
|
(12) |
|||||||
|
|
4 1 2 R3 |
|
|
|
||||||||||||||||
Рассмотрим случай нагружения длинной оболочки погонной кольцевой нагрузкой q в плоскости симметрии xoy при отсутствии силового шпангоута (рис. 2).
Рис. 2. Длинная оболочка, нагруженная погонной кольцевой нагрузкой
Условие равенства работ внутренних и внешних сил для оболочки со ступенчатым изменением толщины записывается в виде
L
2 1 |
dx 2 2 dx 0 |
(13) |
0 |
|
|
где первые два слагаемых – потенциальная энергия деформации нагруженного кольцевого пояса ширины 2 и концевых секций оболочки, T-работа внешних сил для координаты x=0:
|
D |
2 |
|
2 |
hi |
|
|
|
|
|
|
i n2 |
1 |
2 |
ER |
" 2 |
R |
, (i=1,2) |
(14) |
||
2R4 |
|
|
||||||||
i |
|
|
2n4 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
88
|
|
|
|
|
T |
1 |
qкр R 0 H 0 Rd qкр |
n2 1 |
|
2 |
0 R . |
|
|
|
(15) |
||||||||||||||||||||
|
|
|
|
|
2 |
2R |
|
|
|
|
|
|
|||||||||||||||||||||||
С учётом (14) и (15) уравнение (13) получает вид |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
2 |
2 |
|
|
|
|
|
|
L |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
2Di n |
|
|
1 |
|
|
|
2 |
2 |
|
|
2ER |
h1 |
|
|
2 |
|
|
|
|
2 |
|
|
|
2 |
|
||||||||
qкр |
|
|
|
|
dx h |
|
|
|
|
|
|
|
" |
dx h " |
|
|
|
0 (16) |
|||||||||||||||||
|
R |
3 |
|
|
dx |
n |
4 |
n |
2 |
|
|
|
|
dx |
|
||||||||||||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Функцию деформации после потери устойчивости оболочки, удовлетворяющую граничному условию (4), примем в форме (5). Выполняя интегрирование в (16), получим расчётную формулу
q |
|
в 1 t |
C 1 t |
|
|
, |
в 1.5D n2 |
1 KR3 , |
C 2ER3 |
|
K 3 |
n4 n2 1 , (17) |
||
кр |
2 |
h |
||||||||||||
|
1 |
1 |
|
|
2 |
1 |
|
|
|
|
|
|||
|
|
|
|
|
h1 |
|
h2 , |
i |
и ti – коэффициенты; |
(i=1,2). |
|
|||
|
|
|
h |
|
|
|||||||||
Из формулы (17) вытекает ряд частных случаев. Например, для изотропной оболочки постоянной толщины имеем расчётную формулу
|
3 D1 n2 1 |
|
2ER3 |
|
K 3 |
|
||
|
|
h |
|
|||||
|
qкр 2 KR3 |
|
|
n4 n2 1 |
. |
(18) |
||
|
Коэффициент затухания K найдём |
из |
условия минимума |
критической нагруз- |
||||
ки qк |
K 0 |
|
|
|
|
|
|
|
|
K=[Dn4(n2-1)/4ER6h]0.5 . |
(19) |
||||||
Подстановка (19) в (18) с учётом n=2 для длинных оболочек приводится к известной [3] формуле
q |
|
|
12 |
|
Eh |
h 1.5 |
|
||
|
|
|
|
|
|
|
. |
(20) |
|
кр |
9 |
1 2 0.75 |
|
||||||
|
|
|
|
R |
|
||||
Рис. 3. Расчетная схема к решению задачи устойчивости, когда толщина оболочки изменяется вдоль образующей по линейному закону
Решим теперь задачу устойчивости оболочки, представленной на рис.3, когда толщина оболочки h=h(x) изменяется вдоль образующей по линейному закону. Условие равенства работ представимы в виде
|
|
2 dx T 0 |
(21) |
0 |
|
89