CC BY
2Электронный журнал «Труды МАИ». Выпуск № 67
www.mai.ru/science/trudy/
УДК 539.3
К расчету напряжений в цилиндрических сосудах при несимметричном гидростатическом давлении и нагреве Нерубайло Б. В.1*, Ву Ба Зуи1, Зайцев В. М.2**
1 Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3, 125993, Россия Прогресстех, 2-я Хуторская улица, 38а, Москва, 12728 7, Россия
*e-mail: borisn@km.ru ** e-mail: good-for-you@mail.ru
Аннотация
Рассматриваются краевые задачи о напряженно-деформированном состоянии круговых цилиндрических оболочек с различными условиями закрепления на краях при действии гидростатического давления, соответствующего горизонтальному расположению сосуда, частично заполненного жидкостью, и температуры, постоянной вдоль оси и кусочно -постоянной по контуру.
Ключевые слова: цилиндрическая оболочка, гидростатическое давление, сосуд, основное напряженное состояние, краевой эффект, температура
Введение
Определение напряженно-деформированного состояния цилиндрических оболочек при действии гидростатического давления, постоянного по длине, что соответствует горизонтальному расположению сосудов, частично заполненных жидкостью, представляет интерес для многих отраслей техники [1-3]. Например, обечайка топливных баков, предназначенных для размещения компонентов жидкого топлива (окислителя, горючего), в полете нагружена в том числе внутренним избыточным давлением, складывающимся из давления наддува и гидростатического давления. Кроме того, часть обечайки сосуда, свободная от жидкости, может нагреваться иногда до значительной температуры (250...300о С)[3]. Причем,по длине оболочки ее можно принять постоянной, а по контуру на границе сухой и смоченной частей она резко изменяется [1], что может вызвать появление значительных температурных напряжений.
В случае шарнирного закрепления оболочки некоторые результаты расчета напряжений при гидростатическом давлении представлены в монографии В.З.Власова [4], а вопрос о напряженном состоянии при появлении температурного поля не рассматривался.
Постановка задачи
Применение в работе [4] полубезмоментной теории в сочетании с методом двойных тригонометрических рядов существенно сужает возможности изучения оболочек с различными краевыми условиями, так какв случае краевых условий, отличных от шарнирного закрепления, невозможно найти изгибающие моменты и, следовательно, изгибные напряжения на краях, так как там обращается в нуль нормальное перемещение,
связанноес изгибающим моментом недифференциальной зависимостью, ибо равенство нулю перемещения автоматически влечет за собой отсутствие изгибающего момента.
1. Разрешающие дифференциальные уравнения общей теории оболочек. При действии произвольно распределенных по поверхности оболочки нормальной нагрузки р(а, /) и температурного поля Т(а, /), постоянного по толщине, если принять за основу гипотезы Кирхгофа-Лява, исходные уравнения общей теории круговой цилиндрической оболочки могут быть сведены к следующим дифференциальным уравнениям в частных производных восьмого порядка относительно разрешающих функций Ф(а, /), Ф (а, /) [5]:
ЬФ(а,/) = Я4 Б — р(а, /) ,
ЬФ (а, /) = ЯТ(а, /3),
с
(1)
( д4 д4 ^ 1-"2 я4
где Ь = У2У2(У2 +1)2 -2(1 — у)
уда4 да2/2 ,
V2 +
1 -уу_
с2 да4
V2 =--л---, с2 = к 2 /12Я2; а,/ — продольная и окружная безразмерные координаты;
да д/
Я, И - радиус и толщина оболочки; Е, V - модуль упругости и коэффициент Пуассона, а
г
Ек3
коэффициент линейного температурного расширения, Б = ^^-— - цилиндрическая
жесткость.
Все искомые факторы - перемещения, усилия, моменты и т.д. - связаны с разрешающими функциями Ф(а, /), Ф (а, /) дифференциальными соотношениями [5,7].
Представим внешнюю поверхностную радиальную нагрузку и температурув виде следующих разложений:
р(а,/) = р0в(а)^в„ совП//; (2)
п
п=0
T(a,() = T0 в(а)^вп cosn(,
n
n=0
где p0 = —yR, у - плотность жидкости, p0, T0 - амплитудные значения давленияи температуры; 6(a) - функция распределения давления и температуры вдоль образующей; вп, вп - коэффициенты Фурье разложения нагрузки и температуры в ряды Фурье. Закон изменения давления по контуру описывается функцией 6(() [4] icos( — cos(0, p\<(
6(Р) = io, ' <3)
вкоторой Р0 - половина угла заполнения сосуда жидкостью(Рис.1)
.Рис 1. Распеределение кусочно-косинусоидального давления При принятом здесь горизонтальном расположении сосуда функция распределения нагрузки по длине в(а) = 1, а по контуру в(Р) представляется рядом Фурье
ад со
вр) = j вп cos пр = в0+ в cos Р + j в„ cos иР, (4)
и=О
и=2
где вп - коэффициенты ряда Фурье
2 po 2
6o = — Í (cos Р — cos Po)d( = — (sin (o — (o cos (o);
-77" ^ -77"
(5)
2 P 2
6n = — I (cos ( — cos P )cos n(dp = — ni n
sin(n — 1)P sin(n + 1)P sin n(0
' cos Pq
n
2(и -1) 2(и +1)
Здесь первые два члена ряда (п=0, п=1) соответствуют деформированию сосуда без изменения формы поперечного кругового сечения, что считается известным из курса
3
сопротивления материалов.Кроме того, при п=0, п > 1 следует учесть возникающий при этом «классический», осесимметричный (при п=0) и несимметричный (при п > 1) краевые эффекты [7].Прип > 2 имеет место основное напряженное состояние [4,6].
Температурное поле в (2) описывается кусочно-постоянной функцией в окружном направлении оболочки и будет представлено в п.3.
2.Задача о действии гидростатического давления. Решение разрешающего уравнения общей теории оболочек (1) ищем в виде
Ф(a,р) = ^Фп (a)cosnp.
(6)
Делая подстановку разложений (2), (6) в уравнение (1), получим обыкновенное дифференциальное уравнение относительно амплитудного значения разрешающей функции
Ф п («)
ЬФ „ (а) = р0 Я4 01вп в(а),
(7)
где оператор L определяется следующим образом:
L =
da2
- n2 +1
d2
J
da2
-n
- 2 (1 -v)
d
2 С d4
1 -v2 d4
da2
da4
-n
+ -
c2 da
4 '
В соответствии с представлением внешней нагрузки и разрешающей функции разложим в одинарные тригонометрические ряды перемещения, усилия и моменты
u
(a, P) = ^Un (a)cosnP; v(a,P) = ^K (a)sinnP; w(a,P) = ^Wn (a)cosnP;
(8)
T (a, P) = Yjm (a) cos nP; S, (a, P) = £ Sm (a) sin nP; Gt (a, p) = £ Gm (a) cos np, (i = 1;2)
n=0 n=1 n=0
Амплитудные значения Un (a),..., Gin (a) выражаются через ф(a) в результате
подстановки (6),(8) в соотношения, связывающие искомые факторы с разрешающей функцией [7]. Например,
Wn (a) =
^d2 2^
da
- n
Фп (a),
(9)
G2n (a) = -
D_
R2
d2 d2 v4- - n2+1x4- - n2)2
da
da
Фп (a).
Восьмой порядок основного разрешающего дифференциального уравнения (1) и, соответственно, (7) требует постановки на торцах оболочки восьми граничных условий, записанных относительно кинематических или статических факторов.Так что процедура
n=0
J
J
n=0
n=1
n=0
2
отыскания постоянных из системы восьми алгебраических уравнений в общем случае громоздка и таит в себе большую вероятность ошибки. Этого в значительной мере можно избежать, если применить метод асимптотического синтеза напряжённого состояния (МАС). Разрешающая функция Ф„(а) может быть представлена приближённо [6, 7]
i Ф°° (а) + Фк (а) (n < n ) Ф (a)«J n () () ( ) (10)
n ( ) | ФП (а) (n > n ) ( )
(индексы «о», «k», «u» относятся к основному состоянию, краевому эффекту и изгибному состоянию, соответственно).
Однако при рассматриваемой здесь нагрузке, постоянной по длине и имеющей, как
правило,большой угол заполнения Р0 >>у/h / R функция Ф (a) с достаточной точностью может быть определена следующим образом:
Ф ( a V п v '
Ф° (а) при n > 2 п (11) Фk (а) при n > 0
В соответствии с (11) могут быть записаны и искомых факторы w(a, Р) = w° (а, Р) + wk (а, Р);
Т(а, Р) = (а, Р); Т2(а, Р) = Tf (а, Р); (12)
Ох(а, Р) = Gk (а, Р) + -^2°(а, Р); в2(а, Р) = ^(а, Р) + vG^ (а, Р). Полное напряжённо-деформированное состояние, таким образом, складывается из основного состояния и краевого эффекта на основе (11, 12). Причём, каждое из них находится обособленно путем решения соответствующей краевой задачи.
Рассмотрим в качестве примера цилиндрическую оболочку, у которой один край защемленный (а=0), а другой шарнирно закрепленный (а= ах), так что краевые условия при использовании общей теории оболочек формулируются следующим образом: а=0: u(^P)= ^а,Р)=^а,Р)= Wa (а, Р =0 а=а : w(а,P)=V(а,P)= T(а,Р) = G(а,Р) = 0(13)
А в соответствии с идеей применяемого здесь МАС краевые условия (13) должны быть расчленены на следующие две группы [7]: Тангенциальные краевые условия u(o,P)= V (а,Р) =0, ( а=0)
^а,Р)= Тх(а, Р) = 0, (а=а) (14)
Нетангенциальныекраевые условия ^ (а,Р)= м/а (а, Р) =0, (а=0)
^ (а,Р) = 0,(а, Р) = 0, ( а=а) (15)
Если перейти к амплитудным значениям искомых факторов, получим вместо (14),(15) следующую запись краевых условий: и (а) = V (а) = 0, ( а=0)
V» = Т1п (а) = 0, (а=а) (14а)
Ж (а) = Ж' (а) = 0 (а=0) п па
Жп (а) = Оы (а) = 0, ( а=а) (15а)
В качестве подтверждения возможности получения решения на основе (11), (12), приведем числовую информацию, относящуюся к случаю нагружения оболочки, имеющей относительную толщину h/R= 1/100, двумя диаметрально противоположными нормальными нагрузками, действующими по отрезкам контура, каждый длиной 0,5R, примерно равной принятой в [4].
Максимальные значения кольцевого изгибающего момента 02п / Р, имеющие место в середине нагруженных отрезков, представлены в табл.1 при различных значениях номеров гармоник п , на основе различных приближенных уравнений.
Таблица 1
Численные значения кольцевого момента на основе различных уравнений теории оболочек
Номер Общая Основное Краевой Основное Изгибное
гармоники теория оболочек состояние эффект состояние плюс состояние
п краевой эффект
1 0,0170 0,0137 0,0036 0,0173 0,0986
2 0,0144 0,0108 0,0031 0,0139 0,0430
3 0,0118 0,0084 0,0025 0,0109 0,0225
4 0,0080 0,0058 0,0017 0,0075 0,0115
5 0,0040 0,0030 0.0009 0,0039 0,0048
40 0,0001 0,0001 0,0002 0,0003 0,0001
Как видно из табл.1, основное состояние в сумме с краевым эффектом приводят к численным результатам, близким к полученным по общей теории.В то же время изгибное состояние при высоких номерах является ничтожным, что открывает возможность исключить его из рассмотрения в (10) и распространить в получающихся решения(11, 12) суммирование до П ^да [7].
Запишем соответствующие разрешающие уравнения и соотношения, связывающие амплитудные значения искомых факторов и разрешающей функции для основного состояния и краевого эффекта [7].
Рассмотрим вначале основное состояние в (11),(12), которое описывается уравнениями полубезмоментной теории оболочек. Для разрешающей функции Ф (а) получаем
обыкновенное дифференциальное уравнение путём соответствующего упрощения [7] уравнения (7)
( d4 ^ п К2 с2 -- + 4¡/п Фп (а) = —— 9п9(а), 4^и4 =-^п4(п2-1)2; п > 2.(16) Амплитудные значения
ч -а ) Ек 1 -V
перемещений, усилий и моментов связаны с разрешающей
функцией Ф (а) посредством следующих дифференциальных зависимостей:
ип(а) = п2Ф»; Уп(а) = пъФп(а); ВД = пФп(а);
Т71 77/ 7~ч
Т1п(а) = —п2Ф»; Бя(а) = -—п Ф»; С2п(а) = —п\п2-1)Фп(а), (17) К К К
Воспользуемся решением на основе метода начальных параметров [7] для произвольно распределенной по длине нагрузки и запишем его для рассматриваемого здесь случая, положив в нема = а /2. Оно представлено в табл.2. Таблица 2
Решение дифференциального уравнения основного напряженного состояния при действии гидростатического давления
Т1п (о) ^ (0) К(о) и*(0) -0пук2
Тп (а) Ф2 п - (ф1 + Фз) 2 и2 Ф п М(Фз -Ф1) П2 ^ , ч 9 2 ф4 (а)
Sя (а) и ™ (Фз п -Ф1) Ф2 2 а3 -£-( Ф1 п2 + Фз) 2 и2 — Ф4 п [Фх (а) + Фз(а)] 2Ап
Vn (а) п 2,2 ф4 2 п2 4, з(Фз -Ф1) Ф2 п 2U (Ф1 + Фз) 3 fr & - Ф2 (а)] 4Un
U» 1 (Ф1 + Фз) п - 2,2 Ф4 и £(Фз п -Ф1) Ф2 2 -"-г [Фз (а) - Ф(а)] 4Un
Примечание к табл.2: Tln(0), S„(0), V*(0), U*(0) - начальные параметры; Ф1п (а) = ch juna sin juna; Ф 2n (а) = chuecos ¡лпа; Ф3и (а) = shuna sin ипа; Ф 4n (а) = shuna cos una
(индексы "n" у функций Ф1я (а),Ф2я(а),Ф3я (а), Ф4и (а) и коэффициента u опущены).
Амплитудные значения перемещений С/я (а) , V (а) с помощью табл. 2 записываются по формулам
V (а) = Rv»; Un (а) = RKO).
Искомые факторы основного состояния теперь определяются, принимая во внимание найденные из решения краевой задачи амплитудные значения Un (а), Vn (а), T1n (а) Sn (а)
и разложения (8).
В качестве примера приведем окончательные выражения для нормального усилия, нормального перемещения и кольцевого изгибающего момента в оболочке с жестко защемленными краями
T (а, Р) =1 ¿ {-2Ф2 (аТ + [ф (а) + Фз (а)]$, + 2Ф4 (а)} cos nfi,
p0R 4 n=2 U2
E 1 R n n4 - -
(а, Р) =1RЁ -4^n{-2Ф4(а)Т1о - [Фз (а) - Ф1 (а)]й + 2[1 - Ф2 (а)]} cos nfi,
PoR 8 h n=2 u
G2 (а, Р) =1 ¿ -А-{-2Ф4(а)То - [Фз (а) - ф (а^ + 2[1 - Ф2(а)]} cosnfi, P00R 2 n^ n -1
С1 (а )To = Ф4 (а1)Да) - [Фз (а) - Ф (а1 )]В(а),
С-1(а1)5г0 = 2Ф4(а1)В(а1) + [Фз(^!> + ф^МаД (19)
Да) = 1 - ф (а), 2В(а) = ф (а) + ф (а), С(а) = Ф42 (а) +1 [Ф32 (а) - ф2 (а
На основе полубезмоментной теории получены численные результаты, иллюстрирующие напряженно - деформированное состояние оболочки в зависимости от краевых условий, для угла заполнения сосуда жидкостью 2 Р0 = 66 град (Рис. 2-4). В качестве объекта исследования
принята оболочка с геометрическими размерами и углом заполнения, принятыми в монографии [4]Х=4000 см, R=320 см, h=0,6 см. Представленные на рис. 2-4 результаты расчетов позволяют судить о весьма существенном влиянии граничных условий на величину и характер напряженно-деформированного состояния оболочки.Например, величина продольного усилия, и, следовательно, напряжения при наличии жестко защемленного края в 1,5-2 раза превышает напряжения по сравнению с условиями шарнирного закрепления. Причем, это относится к сжимающим напряжениям, при которых возможна потеря устойчивости. Окружное напряжение и нормальное перемещение в средней части оболочки меньше, чем при шарнирном закреплении.
Рис 2. Изменение радиальногоперемешениявдоль нулевой образующей оболочки с различными граничными условиям
О 0,25 0,5 0,75
Рис 3. Изменение нормального продольного усили вдоль нулевой образующей оболочки с различными граничными условиям
Рис 4. Изменение радиального перемешения, нормального продольного усили, концевого изгибающего момента вдоль нулевой образующей оболочки с различными граничными условиям
1. Шарнирные закрепления краев
2. Шарнирный и жестко защемленный края
3. Жестко защемленные краев.
Перейдем к рассмотрению напряженного состояния типа краевого эффекта, включающего не только осесимметричную составляющую (n=0), как это обычно принято при расчете сосудов давления, но и несимметричные (n > 1) формы[7, 8].
Разрешающее дифференциальное уравнение при этом записывается в виде [7]
д 4 w r 4
+ 4^ w = —Р(а, P); 4г,А = 12(1 -v2)(R / h)2; w = w(a, P) (15) da D
Здесь нормальное перемещение w = w(a, P), являющееся функцией двух координат, идентифицируется с осесимметричным краевым эффектом[7] , то есть может быть представлено в виде
n n
w(a, P) = ^ wn (a) cos nP = w(a)^ cos nP
n=0 n=0 (18)
Для нагрузки p(a, P) имеют место приведенные ранее разложения (2),(10). Принимая их во внимание, в случае действия гидростатического давления получим обыкновенное дифференциальное уравнение и силовые факторы
f ^ 4 Л ^ 4
—- + V w„ (a) = pо — 0 0(a) ^ a a J D
Eh D .. D ...
T2n(a)=—wn(a); Gin(a) = wn(a); g2 =vGi; Qm(a) = w„(a)- (19)
Решение по методу начальных параметров получается путем адаптации записанного в [7] решения для случая действующего гидростатического давления. Далее, если длина оболочки
L>>2 L = S^Rh , (20)
где L0 - зона затухания краевого эффекта, целесообразно строить приближенное решение у
интересующего нас края оболочки, рассматривая ее как полубесконечную.
Пусть справедливо сильное неравенство (20), и нас интересует жестко заделанный край,
тогда реальная оболочка моделируется полубесконечной с краевыми условиями из (15а)
a=0: W (a) = W' (a) = 0 , (21)
n na
которые позволяют найти две произвольные постоянные интегрирования уравнения (1 9), принимая во внимание,и получить все искомые факторы краевого эффекта по известной процедуре.
w(a, P) = yRRR- [1 - (p(a) -X(a)] 0(P); T2 (a, P) = yR2 [1 - <p(a) -X(a)] 0(P); (22)
Eh
Gi(a,P) = yR'hW(a)-*(a)]0(P); G2(a,P) = vG,(atP); 0(fi) = £0n cosnp.
6(1 -V ) n=0
По расчетам и в соответствии с физическим смыслом задачи в нижней точке сосуда нормальное перемещение, направлено в сторону увеличения радиуса, кольцевое усилие растягивающее, а изгибающие моменты в заделке создают сжимающие напряжения на внешней поверхности оболочки. Кольцевое усилие на краю равно нулю, а вдали от края принимает максимальное значение, равное найденному по безмоментной теории. Нормальные напряжения на внутренней и наружной поверхностях оболочки находятся по известным формулам [9]..
Найденное здесь напряженно-деформированное состояние следует дополнять балочным напряжением (a,P) = a(a)cosP, (a) = M(a)/ W от возникающего в оболочке, как в
балке, изгибающего момента от гидростатической нагрузки (M - изгибающий момент, W -момент сопротивления сечения оболочки), и прогибом.
З.Термоупругая задача при кусочно-постоянной вдоль контура температуре.При
появлении температурного поля в оболочке(Рис.5)напряженное состояние строится аналогично изложенному вышеслучаю гидростатического давления.
Рис 5
Для безразмерной функции в (Р) изменения температуры Т(а,Р) = Т0в(а)в(Р) по оболочке
в(Р) = 1, р|<р0; в(Р) = 0, р\>р0
при начале координат, помещенном в
диаметральной плоскости, проходящей через середину нагретой области,получаем следующее разложение [7] в ряд Фурье:
в (Р)=А +1 у ^соз пр. (23)
2 п п=1 п
Для разрешающей функции основного напряженного состояния Фи (а) получаем обыкновенное дифференциальное уравнение [7]
Г ^4 Л 1 г2
4 с _ „4/„2 1Ч2
2
*4 4^4 Фп (а)^Т1-Л70 в(а)в 4^4 = -СТп4(п2-1)2; п > 2. (24) 1 -у 1 -V
У') " " ' 1 -у 1 — у"
Амплитудные значения перемещений, усилий и моментов связаны с разрешающей функцией Фя (а) посредством следующих дифференциальных зависимостей:
и (а) = п2Ф»; V» = п3Ф»; ВД = п'Фп(а);
ЕЪ „ ЕЪ -
Т1 = —п Фп(а)- — а^ в„ в(а);
Л 1 — V
ЕЙ - П -
^ (а) = -—п Ф»; С2п (а) = — п\п2-1)Фп (а), (25) К к
При использовании метода начальных параметров решение (24), с учетом (25) представлено в табл.3.
Таблица 3
Решение дифференциального уравнения основного напряженного состояния при действии температурного поля
Т1п(0) Ьп(0) К(а) и* (а) агЕ1г Т ' 0 в 1 -V п
ТтМ Ктт(а) Ктз(а) Ма) Кти(а) Ф2(а-£) -1
Кзт(а) Кзз(а) Ма) Кзи(а) 3
У*(а) Кут(а) Ма) Ма) Куи(а) п 1ТЧ Ф\(а-£\) 2№п
и* (а) Кит(а) Киз(а) Киу(а) Кии(а) -Цф (а-%1) + Фъ(а-^)] 2№п
Построенное в форме табл.3 решение легко реализуется для конкретного случая температурного поля и краевых условий.
Пусть оболочка на обоих концах имеет следующие условия жесткого защемления: и (а, Р) = V (а, Р) = 0, а = 0,а=а,
(26)
ип (а) = V (а) = 0, а = 0,а=а •
Тогда на основании табл.3 и краевых условий (26) получаем продольное усилие, нормальное перемещение и кольцевой изгибающий момент в виде:
1—V 1 " 1
Т (а, Р) = - ^ —{2Ф2 (а)4 + И (а) + Фъ (а)] +
а(ЕНТ л п=2 п
+ 2[ф (а) — 1]} Бт пР0 соб пР;
(27)
1 — V
——^(а, Р) = —-1- £ -Пг{— 2Ф4(а)7—0 + [Ф3(а) — Ф^а)]^ — а ОД 2л п=2
— 2ф(а)}БтпР собпР;
С2 (а, Р) =--аЕТ0 Т" £—{—2Ф4 (а)Т—0 + [Фз (а) — Ф— (а)Ё, —
Л3(1 — V ) 1 — V 6 п=2 П
— 2Ф4(а)}БтпР0 соб пР,
гдеобозначено:
~ (а—)Ф4(а—) + 2Ж(а—)[Фз(а—) — Ф—(а—)] . _ф4(а—)N(а—) — [Фз(а—) + ФМ)Шах)
Т 10 1 5 ^0 1
Ф42 (а—) + -[Фз2(а!) — Ф12(а1)] Ф^а) + — Ф^(аг)]
Обратимся к построению напряженно-деформированного состояния краевого эффекта при неосесимметричном температурном поле Т (а, Р).
Дифференциальное уравнение относительно разрешающей функции (а,Р), описывающее такое состояние, и искомые факторы имеют вид:
^д4 1—v2^ 1 "2
чда4' + с2
w(а, Р) = — аЯТ(а, Р). (28) с
ЕЬ ЕЬ
Т2 (а, Р) = — -—— Ца, Р) — —а ЯТ (а, Р);
( V ) V (29)
в1(а,Р) = — д2^Р); 02(а,Р) = vGl(а,Р). Я да
Здесь w = w(a,P), т.е. является функцией двух координат:
п п
w(a, P) = У wn (a) cos nP = w(a)^ cos nP
n=0 n=0
Для w(a) и силовых факторов получаем обыкновенное дифференциальное уравнение fd4 ^ I-"2
+ 4^4 w(a) = - l—VatRTaer^e(a) (30)
У daA j
с2
ГТ1 , . ЕЙ к/, ЕЙ ч ^ ✓ ч Б ё ^ (а) ^ . .
Т2(а) = ---(а)---а/0в„в(а); Ц(а) = -— , ( ); 02(а) = уЦ(а).(31)
(1 -у) К 1 -V К аа
Полагаем, что рассматриваемые оболочки имеют длину, удовлетворяющую условию (20),
то есть их края не влияют друг на друга, тогда оболочка может быть рассчитана как
полубесконечная, для которой
w(a) = С((а) + с2^(а) -а{КТ0 ;
™'(а) = л{- [((а) + Х(а)] С + [((а) - Х(а)] С}.
При краевых условиях w = 0, ' = 0, (а = 0) на основе (31) получаем систему двух
алгебраических уравнений
c- atRT0 = 0; - atRT0 + c = 0
из которой находим произвольные постоянные интегрирования: С1 = С2 =аДТ0. (33)
На основании (30), с учетом (3),(33) получаемследующие выражения для искомых факторовтермоупругой задачи:
w(a) = [((а) + х(а) -1] а(КТ0; ЕЙ
Т(а) = -¡г w(a) = -[1 - ((а) - ^(а)] аЕиТ; (34)
к
Г( \ п 2 — "Г л л/з(1 -у2) Й2
°!(а) = 2г1 w (а) =1—2—~аЕТ;
К 1 -у 6
^ (а) = уОх (а).
С помощью (34) находим максимальное продольное и кольцевое изгибные напряжения.
6G1(a) _ 73(1 -v2)
(a) = — о = ± л-- atET0; а2(a) = v°"1(a)-
1 h2 1 -v
4. Численный пример. Вычислим термоупругие напряжения для конкретного материалацилиндрической оболочки при гипотетическом уровне кусочно-постоянного по контуру и постоянного по длине оболочки температурного поля[3]. Пусть ^ = 100град, тогда при v = 0,3, at «16-106 град1; E « 2-104 daH / мм2 , что соответствует стали марки 1Х18Н10Т[10],получаемизгибное продольное напряжение ох (a) « 1,815 atET0 « 58 daH / мм2,
котороепочти в два раза превышает предел текучести материала <J02~ 24,8 daH / мм2 [10].
Кольцевое изгибное напряжение: а2(a) = vay(a) «17,4 daH/мм2. Кольцевое усилие, и, следовательно, напряжение, постоянное по толщине, является сжимающим: а2 (a) = T2 (a)/ h = -[1 -y(a) -^(a)]a, ET0. В заделке оно равно нулю.
Из сравнения видно, что изгибное напряжение больше кольцевого, постоянного по толщине оболочки и может достигать в реальных случаях [3] значительной величины. Библиографический список
1. Кан С.Н., Свердлов И.В. Расчет на прочность самолета. - М.: Машиностроение, 1966. 519 с.
2. Новиков В.Н., Авхимович Б.М., Вейтин В.Е. Основы устройства и конструирования летательных аппаратов. М.: Машиностроение, 1991. 368 с.
3. Моссаковский В.И., Макаренков А.Г., Никитин П.И., Саввин Ю.И., Спиридонов И.Н. Прочность ракетных конструкций. - М.: Высшая школа, 1990. 359 с.
4. Власов В.З. Тонкостенные пространственные системы. - М.: Госстройиздат, 1958. 502 с.
5. Образцов И.Ф., Нерубайло Б.В. О методах синтеза напряженного состояния в теории оболочек // ДАН СССР, 1983,№ 3, с. 54-56.
6. Образцов И.Ф., Нерубайло Б.В., Андрианов И.В. Асимптотические методы в строительной механике тонкостенных конструкций. - М.: Машиностроение, 1991.-416 с.
7. Нерубайло Б.В. Локальные задачи прочности цилиндрических оболочек.- М.: Машиностроение, 1983. 248 с.
8. Чернина В.С. Статика тонкостенных оболочек вращения. - М. Наука, 1966. 455 с.
9. Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки.- М., Физматгиз, 1963. 635 с.
10. Безухов Н.И., Бажанов В.Л., Гольденблат И.И., Николаенко Н.А., Синюков А.М. Расчеты на прочность, устойчивость и колебания в условиях высоких температур(Под редакцией И.И. Гольденблата).- М., Машиностроение, 1965. 567с.
