МНОГОСЛОЙНЫЕ СВАРНЫЕ КОНСТРУКЦИИ И ТРУБЫ
АНАЛИЗ РАБОТЫ РУЛОНИРОВАННОЙ СТЕНКИ СОСУДА, ОСНОВАННЫЙ НА УСРЕДНЕНИИ СВОЙСТВ НАВИВКИ
В работах [1, 5] предложена схема рулонированной стенки сосуда, основанная на усреднении свойств навивки, показан путь идеализации, приводящий к схематизации стенки трехслойным цилиндром,, а также исходные уравнения и полученная с их помощью разрешающая система дифференциальных уравнений, записанная в нормальной форме. При этом жесткостные характеристики слоя, схематизирующего навивку, представлены в общем виде, чем предусмотрена возможность различных вариантов усреднения. В настоящей статье конкретизируется усреднение зависящей от микронеровностей контактной податливости между витками навивки и исследуется работа схемати-
Рис. 2. Усреднение свойств навивки проскальзывания (а) и контактной податливости (б).
Рис. 1. Схема идеализированной конструкции. |
вирующего ее слоя. Цель данного исследования — углубление понимания некоторых аспектов работы рулонированной стенки.
Путем, описанным в [1], приходим к идеализированной конструкции (рис. 1). Как видно, навивка представляет собой спираль Архимеда из пакета тонких листов (штриховые линии на рис. 1), толщина которого равна толщине h навиваемой полосы в реальной конструкции. Устремляя к нулю толщину тонких листов при сохранении общей толщины пакета h, заменяем всю навивку сплошным слоем с анизотропными свойствами. Таким образом, расчетная схема всей стенки представляется трехслойным цилиндром. Внутренний и наружный слои, соответствующие центральной трубе и кожуху, изотропны.
Возможность проскальзывания в навивке усложняет задачу о напряженно-деформированном состоянии описанного трехслойного цилиндра при изменении напряжений вдоль его длины, что имеет место в реальной конструкции из-за наличия в стенке кольцевых швов, связывающих отдельные обечайки между собой и с днищами сосуда. Поэтому на данном этапе исследования ограничиваемся более простой задачей — рассмотрением бесконечного цилиндра, нагруженного внутренним давлением р и осевой силой Z. Напряжения и деформации в нем не зависят от осевой координаты г.
При определении свойств слоя, схематизирующего навивку, усредняются проскальзывание между витками и жестностные характеристики.
Для усреднения проскальзывания рассмотрим элемент навивки, представляющий собой короткий отрезок, состоящий из двух половинок (по толщине) ее смежных витков. На рис. 2 показано их взаимное смещение в результате проскальзывания. При расслоении навивки в процессе схематизации на более тонкие листы это смещение распределяется ступенчато и в пределе переходит в непрерывную деформацию сдвига Yp,.
Усредненные жесткостные характеристики навивки зависят от микро - и макронеровностей контактирующих поверхностей. Макронеровности вызываются разнотолщинностью полосы и ее погнутостью. Различные виды этих дефектов не одинаково влияют на жесткостные характеристики многослоя в целом.
Исследования контактной жесткости, зависящей от микронеровностей (шероховатости), описаны в работе [2]. Применительно к многослойным сосудам исследования контактной податливости от сжимающих напряжений проводились в Иркутском НИИхиммаше [3, 4].
Усреднение жесткости выполняется следующим образом.
Рассмотрим элемент навивки, состоящий из двух половинок смежных витков и нагруженный нормальными сжимающими напряжениями q (рис. 2, б). Вследствие микронеровностей (шероховатости) контактирующих поверхностей в какой-то узкой зоне контакта деформируемость будет большей, чем в остальной части элемента. Для дополнительного по сравнению со сплошным материалом сближения б нагруженных граней показанного элемента в работах [3, 4] предлагается следующая эмпирическая зависимость от давления q
8 = Aqa, (1)
где А, а — постоянные величины, определяемые экспериментально на образцах из пакетов листов. Полная усредненная по толщине элемента деформация сжатия gp определяется формулой
ер = 1 + 4- <2>
Ей соответствует усредненный модуль упругости
Яр “V - <3>
Здесь Е — модуль упругости материала навивки, поперечная деформация и модули упругости в направлении навивки Es и вдоль оси цилиндра Ег остаются такими же как и для сплошного материала, т. е. Es = Ez = Е. Аналогичным образом можно получить и усредненную жесткость на сдвиг, характеризующуюся модулем сдвига Gps. Однако экспериментальными данными о контактной податливости сдвигу мы не располагаем.
Для описания напряженно-деформированного состояния в точке слоя, схематизирующего навивку, выделим из него в полярной системе координат г, 0 бесконечно малый элемент (рис. 3). Здесь же показаны координаты р, s, направленные по нормали и касательной к навивке. Угол а между осями этих систем определяется из уравнения
спирали Архимеда г = 0 соотношением
sin а — -=s. (4)
Y 4яага-)-Аа
На рис. 3 показаны также действующие по граням элемента радиальные ог и окружные ае напряжения. По площадкам, перпендикулярным осям р и s, кроме нормальных напряжений стр и aSf действуют и касательные напряжения tps, а в поперечном сечении слоя — только осевые напряжения ст2.
Выпишем все соотношения, характеризующие напряженно-деформированное состояние выделенного элемента. Ограничения, налагаемые на напряжения условиями контакта витков навивки, |
Рис. 3. Элемент слоя, схематизирующего навивку. |
(5) |
Ор <0, |
</, |
где / — коэффициент трения. При отсутствии проскальзывания удовлетворяется неравенство, при его наличии — равенство, которое записывается в виде уравнения тр5=±/аР, (6) |
где знак зависит от направления гР8. Дифференциальное уравнение равновесия |
<Эсг, |
<7- — Of. |
(7) |
= 0. |
+ |
дг ‘ г Соотношения между напряжениями в системах координат г, 9 и Pi 5 ог = ар cos2 а + os sin2 а + 2rps sin a cos а, otq = Op sin2 а + as cos2 a — 2xps sin a cos a, (8) Tr0 = (crs — Стр) sin a cos a + Tps (cos2 a — sin2 a) = 0. |
Соотношения между напряжениями и деформациями ер = - g д— (cts + tf*), д (Op + а2), |
(9) |
(ар + аа) = const, урз = |
+ ?ps- |
Е |
Е |
ps |
Здесь v — коэффициент Пуассона для материала навивки. Соотношения между деформациями в системах координат г, 0 и р, s
ег = ер cos2 a - f es sin2 a + yps sin a cos a,
e0 = єр sin'2 a - j - e„ cos'2 a — Yps sin a cos a, (10)
<y>r0 = 2 (es — 8P) sin a cos a + yps (cos2 a — sin2 a).
Уравнение совместности деформаций
а80 dr |
(11) |
= 0. |
+ |
во — є, |
Выражения деформаций ег, ее и Yre через радиальное и и окружное v перемещения
dv Y г0 — дт |
да дг |
Уравнения (6) — (12) дополняются интегральным условием равновесия по поперечному сечению всего трехслойного цилиндра
(13) |
2л Г azrdr = Z,
где г0, г3 — соответственно внутренний и наружный радиусы цилиндра.
С помощью соотношений (5) — (13) проведем анализ работы слоя, схематизирующего навивку. Рассмотрим этот слой при действии на него нагрузок в виде радиальных напряжений — стгі на внутренней и — аГ2 на наружной поверхностях. Знак минус указывает на то, что эти напряжения сжимающие.
Уравнения статики (6) — (8) позволяют определить все напряжения в слое, кроме о2, с точностью до одной постоянной интегрирования независимо от деформаций слоя, т. е. при наличии проскальзывания по всей его толщине слой является статически определимым. Для этого достаточно знать нагрузки на его поверхностях. Поскольку решения для напряжений содержат лишь одну постоянную интегрирования, то эти нагрузки не могут задаваться произвольно, между ними должно существовать определенное соотношение.
h
±2я 1~ CTp cjf> -' С | є Т ±2 Я/-Ї- |
Если пренебречь величинами порядка — по сравнению с единицей то можно получить простое решение указанных уравнений и представить выражения для напряжений в следующем виде:
(14)
t, jS = ± /ар » ± /Схе Здесь Сх — постоянная интегрирования.
Таким образом, отношение напряжений на границах слоя ~— должно находиться в пределах
Предельные случаи (равенства) соответствуют проскальзыванию по всей толщине слоя. Если выполняются неравенства, то проскальзывания может не быть или проскальзывает только часть слоя по тол-
о
щине. ото зависит от величины —- и упругих характеристик слоя.
rt
В этом случае слой в целом не будет статически определимым. Для той его части, в которой нет проскальзывания, уравнение (6) неприменимо. Ее напряженно-деформированное состояние описывается
системой уравнений (7) — (13) при Yps = 0. В случае проскальзывания по всей толщине слоя по найденным напряжениям (причем напряжение ст2 остается неизвестным) с помощью формул (9) и (10) можно найти выражения для деформаций ер, es, Yps и єг, ее, уиэ. В них будет входить неизвестная деформация проскальзывания Yps и неизвестная осевая деформация е2, являющаяся постоянной величиной.
Подставляя єг и еэ в уравнение совместности деформаций (И), получаем дифференциальное уравнение первого порядка относительно Yps. Его общее решение можно представить в виде суммы двух частных решений
Yps = Yps(1) (Or, <Т0, e2) + Yps2). (16)
Первое слагаемое — частное решение неоднородного уравнения, зависящее от напряжений и осевой деформации. Получить его аналитическим путем затруднительно, так как напряжения входят не только в правую часть, но и в коэффициенты уравнения. Выражаемая этим частным решением деформация проскальзывания дополняет упругую деформацию так, что в сумме они удовлетворяют условие совместности деформаций.
Второе слагаемое — частное решение уравнения при отсутствии напряжений, т. е. выражаемая им деформация проскальзывания может иметь место без нагружения слоя. Таким образом, при проскальзывании по всей толщине слой не только статически определим, но и геометрически изменяем.
В аналогичной выражению (16) форме можно представить решение для радиального перемещения и, т. е.
u = u(ii(ar, ае, е2) + и(2 (17)
а в общем решении для окружного перемещения и, получаемого интегрированием третьего из уравнений (12), появляется еще одно частное решение, т. е.
V — У(1> (стг, Сте, е2) + 1>(2) + v(3 (18)
Если пренебречь величинами порядка - у- по сравнению с единицей, то для Yps(2) можно получить простое приближенное выраже
ние, а именно:
, (19)
где С2 — постоянная интегрирования. Через него также просто выражаются решения и(2) и z/2). Эти выражения имеют вид
vW~-C2. (20)
Частное решение у(3) — это решение однородного уравнения, соответствующего (12). Оно выражает жесткий поворот слоя в целом и имеет вид
= С3Г, (21)
где Са — постоянная интегрирования.
Геометрической изменяемостью объясняется то обстоятельство, что слой не может воспринимать нагрузок, выходящих за пределы интервала, определяемого выражением (15). Если будет тенденция
уменьшить отношение ниже левого предела, то «сработает»
проскальзывание (19) и слой за счет него будет расширяться без уменьшения этого отношения. При тенденции к увеличению отноше - аг
ния —- выше правого предела проскальзывание будет происходить
в обратном направлении.
Таким образом, если в процессе нагружения проскальзывание
аг
происходит по всему слою навивки, то отношение -^-2- остается неизменным, т. е. имеет место простое нагружение слоя. В этом случае напряжения в нем пропорциональны нагрузке. Сложное нагружение
слоя, характеризующееся изменением отношения ——, может осу-
гі
ществляться в рамках соотношения (15). При этом слой разделяется на зоны с проскальзыванием и без проскальзывания. Граница между зонами смещается в процессе нагружения, проскальзывание может прекратиться совсем. Аналогичное явление имеет место и при простом нагружении, вызывающем проскальзывание не по всему слою. Это происходит вследствие нелинейной контактной податливости зоны, в которой отсутствует проскальзывание. При идеальном контакте [Ер = Е и Gps = 2(Е--) ) гРаниЧа межДУ зонами с проскальзыванием и без проскальзывания фиксирована, а напряжения и деформации пропорциональны нагрузке. В процессе разгрузки зона с проскальзыванием изменяется как при нелинейной контактной податливости, так и при идеальном контакте.
Следовательно, при сложном нагружении слоя, схематизирующего навивку, его напряженно-деформированное состояние зависит не только от величины нагрузки, но и от истории нагружения.
Нагрузки стГі и аг, — это напряжения взаимодействия среднего слоя трехслойного цилиндра с его внутренним и наружным слоями. Нелинейная жесткость среднего слоя приводит к тому, что даже при
Ягг
простом нагружении цилиндра в целом отношение —— изменяется
® г
в процессе нагружения и, таким образом, для среднего слоя нагружение является сложным. При расчете цилиндра на внутреннее давление приходится прослеживать весь процесс нагружения, пошагово увеличивая нагрузку. В случае идеального контакта задача упрощается, так как напряжения и деформации пропорциональны нагрузке, если, конечно, в процессе нагружения не допускается частичная разгрузка и в ненагруженном состоянии цилиндр не обладает предварительной напряженностью. При полной разгрузке цилиндр без предварительной напряженности возвращается в исходное ненапряженное состояние как при идеальном контакте, так и при нели
нейной контактной податливости в навивке, но напряженно-деформированное состояние не повторяет при этом пути, который оно проходит при нагружении, так как часть энергии расходуется на преодоление трения при проскальзывании. Предварительная напряженность может уменьшить или даже вовсе предотвратить проскальзывание.
Рис. 4. Распределение окружных напряжений по толщине стенки трехслойного цилиндра (п — число витков): |
1 — с учетом проскальзывания, контактной податливости и истории нагружения; 2 — с учетом проскальзывания; 3 — с учетом только контактной податливости; 4—решение задачи Ляме. |
При учете контактной податливости и истории нагружения расчеты трехслойного цилиндра, схематизирующего рулонированную стенку, весьма усложняются. Для их выполнения составлена программа на языке АЛГОЛ-бО, которая предусматривает численное интегрирование дифференциальных уравнений по методу Рунге — Кутта без внесения в них упрощений, использованных для получения приведенных выше приближенных решений. В процессе интегрирования удовлетворяются условия сопряжения слоев и отыскивается граница между зонами с проскальзыванием и без проскальзывания. Условия сопряжения слоев выражаются равенством их
радиальных напряжений стг, а также радиальных и и окружных v
перемещений, которое удовлетворяется также на границе зоны проскальзывания. Для ее отыскания используется то обстоятельство, что для обеих зон на их границе справедливы уравнения для обоих случаев — при наличии и отсутствии проскальзывания.
Проведенные расчеты показали, что проскальзывание при внутреннем давлении существенно зависит от коэффициента трения и жесткости наружного слоя, т. е. толщины кожуха. Причем его толщины, равной толщине навиваемой полосы, достаточно, чтобы проскальзывание заметно сказывалось на напряженно-деформированном состоянии цилиндра лишь при весьма малых значениях коэффициента трения — порядка 0,06 и менее. Если довести толщину кожуха до трех толщин навиваемой полосы, то при коэффициенте трения / ^ 0,05 проскальзывание практически не влияет на напряжения и деформации. Значительно большее влияние на них при таких значениях коэффициента трения оказывает контактная податливость, вычисленная по данным работ [3, 4]. При этом оказывается, что с ростом контактной податливости увеличивается и проскальзывание.
Результаты расчета трехслойного цилиндра, схематизирующего рулонированную стенку, изготовленную из стального прокатного листа, показаны на рис. 4. Размеры цилиндра следующие: г0= 400 мм; rt = 410; г2 = 470; г3 = 480 мм; средний слой его соответствует десяти виткам навивки толщиной h = 6 мм. Параметры, характеризующие жесткость такие: Е = 2,1 ■ 106 МПа; v = 0^; А = 0,0505;
d = 0,222 [3]. Модуль сдвига, Сра ввиду отсутствия соответствующих экспериментальных данных по контактной податливости сдвигу опре-
Е
делялся по формуле Gps = 0уг~г~ • Заметим, что величина Gps мало
Z (1 - J - V)
влияет на напряженно-деформированное состояние рассматриваемого
Е
цилиндра (результаты расчетов при принятом Gps и при Срз = 2 ц_ vy
°0 гт
практически совпадают). Результаты — приведены для р = 39 МПа
и / = 0,06. Спады кривых в районе восьми—десяти витков указывают на зоны проскальзывания. Как видно (рис. 4), учет контактной податливости и проскальзывания почти вдвое повышает напряжения
<Jq
в центральной трубе — = 10,5, по сравнению с решением задачи Ля-
°0 р г
ме — = 5,5, при атом основную роль играет контактная податливость.