УДК 620.179.101
ЦИКЛИЧЕСКОЕ СЖАТИЕ МНОГОСЛОЙНОГО МНОГОПРОЛЁТНОГО
ГОФРИРОВАННОГО ПАКЕТА
© 2013 И. Д. Эскин, Р. И. Алкеев, В. И. Иващенко
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Методом Бубнова - Галеркина выполнено решение задачи о циклическом сжатии многослойного многопролётного гофрированного пакета с учётом поворота сечений пакета, лежащих в вершинах гофров. Получены соотношения, описывающие любые процессы загрузки в поле упругогистерезисных петель пакета.
Многослойный многопролётный гофрированный пакет, циклическое сжатие, метод Бубнова -Галеркина.
Имея решение для многослойного одиночного гофра [2], можно построить решение для многослойного многопролётного гофрированного пакета.
Условные обозначения
1 = 1, 3, 5, 7, ... - порядковые номера процессов нагружения, 1 = 2, 4, 6, 8, ... - порядковые номера процессов разгрузки,
] = 1, 2, 3, ..., т - номер проскользнувшего гофра и номер этапа загрузки пакета, т - число гофров в пакете, п - число лент в пакете, (*) - индекс, указывающий, что значения параметров взяты в конце этапа загрузки пакета (для одиночного гофра в конце процесса загрузки), Ъ - толщина ленты пакета, Ь - ширина пакета, /- начальный выгиб гофра ленты, Уо -текущий начальный выгиб гофра ленты,
Г - шаг гофра ленты,
I = ЬЪ /12 - момент инерции поперечного сечения ленты,
Е - модуль упругости материала лент пакета,
ц - коэффициент трения скольжения на контактных поверхностях пакета, к = 1, 2, 3, ..., т - порядковый номер гофра,
Р1 - текущее значение циклической силы,
действующей на одиночный гофр на 1-й
загрузке,
*
Р I - значение циклической силы, действующей на одиночный гофр в конце 1-й загрузки,
х - координата текущего сечения гофра, у. - текущее приращение деформации к-го гофра на. этапах 1-й загрузки в сечении
х,
у.. - текущее приращение деформации на
] этапах 1-й загрузки в серединах пролетов пакета при х = Г/2,
у. - текущее значение деформации на .
этапах 1-й загрузки в серединах пролетов пакета при х = 1/2,
А у. - текущее значение деформации на
.-м этапе 1-й загрузки в серединах пролетов пакета при х = 1/2, УШ]. - текущее значение остаточного выгиба к-го гофра на .-м этапе 1-й загрузки в сечении х,
уУ0. - остаточный выгиб на .-м этапе 1-й
загрузки в серединах пролётов пакета при х = Г/2,
Рк.1 - сжимающая сила, действующая на к-й гофр на.-м этапе 1-й загрузки, Г. - сжимающая сила, действующая на пакет на.-м этапе 1-й загрузки, 5(х - Г/2) - импульсивная функция Дирака, которая определяется следующим образом:
5 = 0 при г/'2 < x < г/2, д=<х> при x = t / 2,
|8ах =1,
£ = Ък - площадь поперечного сечения гофра,
Ыц - продольная сила, растягивающая или сжимающая одиночный гофр на первом этапе 1-й разгрузки или 1-й повторной нагрузки,
Ыу - продольная сила, растягивающая или сжимающая к-й гофр на у-м этапе 1-й повторной нагрузки в случае, когда у < к.
Первичное нагружение. При первичном сжатии пакета гофры распрямляются и точки контакта смещаются симметрично от середины. При этом на каждую у-ю половину гофра, считая от концов пакета, действуют поперечные силы Pjl/2 и продольные силы:
где к=0,1,2,.. ,,т - текущий номер вершины гофра, и в случае, когда т1 < т, - в виде:
' ' (3)
~ ~ рк ~ р рк у = у сое-= у — cos
У 1к ./10 2 члл
т1
г
2
т1
/X У X У1 X УК
л. Р
-2
Рис.1. Схема нагружения
многопролетного гофра:
а) нумерация гофров и их вершин,
б) первый полугофр (¡=1)
(
т
i=у-1
I р+
р
Л
(1)
Распределение углов поворота сечений лент, лежащих в вершинах гофров, описывается некоторой плавной функцией, симметрично убывающей от концов пакета к его середине. Причём можно предположить, что углы поворота концов лент определяются соотношением (6) из [2]. Справедливо также считать, что при достаточно большом числе гофров, начиная с да1-й вершины, сечения лент, лежащие в вершинах гофров, не поворачиваются. Здесь также номера вершинам присваиваются симметрично средней вершине пакета, начиная с концевых, которым присвоен номер «0» (рис. 1).
С учётом высказанных соображений функцию распределения углов поворота сечений, лежащих в вершинах гофров, в случае, когда т1 > т, запишем в виде:
у
~ рк у cos —
* 10 2к
р
=у—
рк 2т
(2)
Примем, что углы поворота сечений вершин с номерами, большими ть равны нулю. Силу, действующую в процессе первичного нагружения на гофр, составленный из двух первых симметрично расположенных половинок (рис. 1), найдём по (28) из [2] при =1 и кэ=1,5. Определим силу, действующую на любой у-й гофр, составленный из у-х половинок, номера которых содержатся в интервале 2 < у < т1. Дифференциальное уравнение, описывающее изгиб у-го гофра, запишется в виде:
+ т(
ах
Рп + ... + р1
' ^ = Р/ ( х -1 ]
ах Р} ^ 2 0
где
/ (, 2рх
У01=-211 - У1.
(4)
(5)
Решение уравнения ищем в виде:
а < (, 2рх } . рх
У = ~ТI1 - СМ"Г | + Ау
(6)
2
1=1
г
где й, и А, - постоянные, определяемые из граничных условий:
' /ЛЧ ' , ч ~ Р Р (, -1)
У! (0) = -У! (0 = у -С08-^ 7
У i (л 2лх 1
У = -Т I1 - cos"T
(12)
t
2
mi
t
У1121 = У г
(7)
Функция (6) удовлетворяет также граничным условиям:
У! (0)= у! )= 0,
1 I 2.
У1 1-1 = 0.
(8)
Найдя dj и Aj из (7), окончательно получим
У1 (1 — cos-)| 1 — cos- I -Г
(9)
У =—1 (1 - cos——-:
Л 2 2 m1
В результате получим значение силы, действующей на j-й гофр, в виде:
2Со~! - У - /) (13)
р __4=1_1___4 '
1 2- /)
]= т1+1, т1+2,...,т.
Соотношения (28) из [2], (10), (13) позволяют последовательно найти силу, действующую на первый, второй, третий и так до последнего т-го гофра. Силу, действующую на пакет, найдем как равнодействующую этих сил
т1 т
=Рп+1: рл +1 р..
~ Р (j -1) . РХ
+ У cos— -- sin -
j=2
j'=m1+1
2m
в виде:
t
Решая уравнение методом Галерки-на, получим
F1 =
1,5С „у 1
4 - ^ f f - У1
t I Зж
Pj1
2Соу (1 - cosР(j-1)] + 4 С0Л1 2m1 J 3 -2т* &«! ж t k=1 у ж (j -1) У cos^- ^ 2m, 1
2 +1 p 2 t у ( ж (j "1)1 f У 1 cos -'- I- f Г \ 2m1 0 _ 4 у ж ж (j-1) + тУ t °os У 3 Л1 2 m1
2С 0 Aj + up- (f - у 1 Bj P*
j=2 2+тж 2t
У11 — p2mT|-f
4 У1 mp ж (j -1)
+ -cos — -
3t 2m1
j-1 Р2 uL рР k=1 1 У1 p (j -1)' cos^- 2 m1 0 - f 2 1 » 2С 0у 1 - up (у 1 - flpH (14) + l ' k=1 ,
2+1 тР2 2 t " у f Р (j - 1) ГЛ 2m1 - f 4 у ж ж (j -1)' j=m1+1 2+mp- У 1 - f) + — Uy cos v ; 2t v 1 ' t 2m, 3 1 m1 где
(10)
где Со определяется по (33) из [2].
Определим силу, действующую на любой j-й гофр, составленный из двух симметрично расположенных j-х половинок, номер которых содержится в интервале т1 + 1 < j < т.
Дифференциальное уравнение, описывающее изгиб гофра в этом случае, запишется в виде:
A = 1 + 1--11 cos
Аj 1 Зж 1
ж (j -1)
Bj =1 -l1 - ЗРIcos
2 m1 Р(j -1)
2
(15)
m1
EJ
^ + '(P11+-+Pm'+" + p~+yj^ " j iX - 5
p1 d у01
(11)
j= m1+1, m1+2,...,m,
где у01 определяется из (5).
Решение (11) ищем в виде, аналогичном (8) из [2]:
Поскольку силы Р^ отличаются незначительно, упростим соотношения (10), (13) и (14), заменив в них, как и в [1], продольную силу, действующую на j-ю половину гофра, приближённым значением. Тогда для силы, действующей на j-й гофр, с номером, расположенным в интервале 2 < j < т1, получим
р С 0~ А,. (16)
]=2,3,...,ть
+
+
+
Для силы, действующей на j-й гофр с номером в интервале т1+1<]<т:
С„Л (17)
Рп =
1 +
тк
4/
(2 , - 1)(у 1 -I)
Сила, действующая на гофрированный пакет в процессе первичного нагру-жения, определяется из соотношения:
1,5 с„~ ,
^ 1 = ■
4 - К ( 8К - ~ / 13к *1
С „ А, ~!
т,
+ И-2
]=21 + ^К2
2] -1 4
с „~\
1 В] - I)
(18)
] - + ^ (2] -1)(~! -1)
Значение параметра т1 определено экспериментально (результаты экспериментального расчетного исследования приведены в [3]) и равно т1=2. Подставив это значение т1 в (18), для силы, действующей на пакет в процессе первичного нагружения, окончательно получим
„ 1,5 + А2 + ._1_ Г ~
F1--г.-" +—:—2-+ .-2- С„у1,
4 - ^ (31 - 1 + ^ ~ В2 -1) ]-31 + ^ (2, -1)(~. -1)
(19)
где значения А2 и В2 получим из (15) при 7=2.
Разгрузка. На первом этапе разгрузки все гофры деформируются без смещения их вершин и поворота сечений, лежащих в вершинах. Здесь так же, как и на последующих этапах разгрузки, принимаем, что пакет деформируется, как одна лента с жёсткостью, равной сумме жёст-костей отдельных лент. Возможность использования этого допущения для многослойного многопролетного гофрированного пакета подтверждена хорошим совпадением экспериментальных и теоретических результатов. Кроме того, модель многослойного многопролётного гофрированного пакета, у которого в каждом пролёте все слои проскальзывают относительно друг друга и опорных плит одновременно, но проскальзывание постепенно распространяется от крайних пролётов к среднему, является вполне приемлемой.
На основании теоретических результатов, полученных для многослойного однопро-лётного гофра, можно утверждать, что вышеуказанное допущение будет справедливо и для этой модели многослойного многопролётного гофрированного пакета.
Текущее значение приращения деформации гофров на первом этапе разгрузки запишется в виде (8) из [2].
Продольная сила, действующая на первый гофр, составленный из двух крайних полугофров, определяется из соотношения (13) из [2]. Продольную силу, действующую на первом этапе разгрузки на й гофр, составленный из двух симметрично расположенных ]-х половинок, определим по (13) из [2], предварительно определив удлинение упругой линии гофра по (8) из [2].
Тогда для гофров, номер ] которых содержится в интервале 2 < ] < т1, продольная сила запишется в виде:
м - 3/С„ у1
1 V ,12 2,2
к И
Для гофров с номером ], содержащимся в интервале т1 + 1 <] < т:
3/ С„~!
1 - У1 В,
(2„)
N,12 -
2 , 2 к И
1 - У1
(21)
Дифференциальное уравнение, описывающее изгиб первого гофра на первом этапе разгрузки, запишется в виде:
Жх
+т Р,
Ж У1 Жх
112 7 2 ах
11 = ( Р1
12'
Р
:> (х - 2
(22)
Для любого ]-го гофра запишется в виде:
а - м ау
+т
Жх
' * •
Р
. +Р (=1 ^
ах
ж2 У.2
Жх
= ( Р,12-р
1х - 2
(23)
]=2,3,4,...,т,
где у12 определяется по (8) из [2], N112 - по (13) из [2], N,12 при 2 <] < т1 - по (2„) и при т1 + 1 < ] < т - по (21). Остаточный выгиб первого гофра на первом этапе разгрузки:
+
+
+
+
И. ~ У, 2ж) -- . ж
-.V - у Л1 -со§— 1+у^-,
У 0112 У011 + У12 2 У
остаточный выгиб /-го гофра:
(24)
У,
0 /12 2
У 1
1 - - i-V- у 1;
2^1
2рх 1 - сое-i +
у/ . Р(/ -1) . Рх + У --вт—,
2ш1 '
2 < / <
(25)
Г =
- -12 {Со
2
2 л2
(т -Ш1 ^ - + ¿'(V - у' В
¿1т -1 + ~! рЛ
(30)
где силы Р*/1 определяются при у = у! из соотношений: (28) из [2] при /=1 и 1=1; (16) при 2 < / < т1 и (17) при т1+1 < / < т.
Учитывая, что т!=2, окончательно получим
у0/12 2 VУ1 f У
т1 +1 < / < т.
2ях 1 - сов-
12 Л /
(26)
Уравнения (22) и (23) решим методом Галеркина. Правые и левые части уравнений умножим на 1-со^ 2лхЛ и проинтегрируем их в пределах от 0 до t.
В результате, пренебрегая членами второго порядка малости, найдём силы, действующие на первом этапе разгрузки на гофры, составленные из симметрично расположенных полугофров, в виде:
на первый гофр
= • _ ~ 1 Р112 Р11 у121С0
1+(f - 3Р у
2t
-У^^РпГ,
(27)
на любой у-й гофр
Р/12 Р;1 У12 ]С0
2 < / < т„
'+(f - У
2t
Р.. + ...+ Р/
(28)
Р/12 Р/1
1+^(V - /Г - Ж [р;1+...++...+Р'и+Р 2Л V 0 2t 11 Гт11 г1-11 2
т1 +1 < / < т.
(29)
Г и = Г,
- у12 {с 0 2
т+ф г - 2)/-+^f - у*) +(f -
^ Мт -1+#4
(31)
Первый этап разгрузки закончится тогда, когда продольные силы преодолеют силы трения на крайних вершинах пакета и произойдёт проскальзывание этих вершин. Определим силу, действующую на пакет в конце первого этапа разгрузки.
Приращение деформации гофров в конце первого этапа разгрузки можно определять двумя путями: из условия
N =т Рп+Р2)
(32)
и условия сшивки процессов разгрузки первого гофра при неподвижных и подвижных концевых сечениях
4 +
тр (8
Т v 3р г - у
р - -12с1+V - ¿у;)2 -
(33)
Силу, действующую на пакет на первом этапе разгрузки, найдём как равнодействующую этих сил в виде:
Из (33), учитывая, что деформация пакета в конце первого этапа разгрузки
у/гу;--;,' (34)
и отбросив величины второго порядка малости, найдём
3
/=1
1
2
у 1
^ 1 1 1,5 С оу 1 - Р и 2 [4+иж (38ж / - у ;)
1,5 с о+р : МЖ~-{ с о [1 + 2к2 (/ - 3ж у *)] 2 _ Р • 1 Р11 21 | 2 4 + ^ ( 38ж / - у1)
(35)
Подставив (35) в (31), определим силу, действующую на пакет в конце первого этапа разгрузки. Деформацию пакета в конце первого этапа разгрузки найдём, подставив (34) в (35). При дальнейшей разгрузке проскальзывания постепенно будут распространяться от крайних гофров к среднему гофру.
Так как т1=2, для полного рассмотрения процесса разгрузки достаточно рассмотреть ещё второй этап, когда первый гофр уже проскользнул, j-й этап, когда проскользнул j - 1-й гофр, составленный из расположенных симметрично относительно середины пакета полугофров (номер j содержится в интервале т), и последний этап разгрузки, когда все гофры пакета проскользнули.
Рассмотрим второй этап разгрузки. Силу, действующую на первый гофр, составленный из симметрично расположенных полугофров, найдём по соотношению
(28) из [2], подставив в него вместо у значение деформации пакета на втором этапе разгрузки
У 22 = У *2 "ДУ 22.
(36)
^ 0 1 У12 -Д У 22
„ иж ( 8 , У* 4 у
4 + "—I—f - у + Д у t 13ж У12 У 2
Дифференциальное уравнение, описывающее деформацию любого Ъ-го гофра (1 < k < m), запишется в виде:
74 у 72
-ы йУ
йх
йх
^ Р-+^) ^
Р: 22 - \х - 2
(37)
где Ык22 определяется из соотношения (20), : = 2 и из (21), если 2 < : < т, при подстановке в них у22 вместо у12 и
2211 - ^
2жх
У 22 = - У2
Решив уравнение (37) методом Га-леркина, получим силу, действующую на второй гофр:
Р222 = Р21 -(3У1*2 +Д У 2
(38)
С0 ^ + ¥
иж
2t
р + Р21
11
и на любой Ъ-й гоф
Ръ22 = Рй -((У~12 + Д У„КСо
у.2+д у 22л / - у. в2м г- у, д
1) 2
р:
1+
ъ_
2 И
f - Ух^'1 f - Уа + Д У 2:
иж- (р;+...+р:„+р 2t и 11 2
Определим силы, действующие на каждый из т - 1 непроскользнувших гофров.
Р 22
2 < : < т. (39)
Силу, действующую на пакет на втором этапе разгрузки, найдём, просуммировав все усилия, действующие на отдельные гофры:
-Со (Уа+Д У 22Jim
2к
+(У12+Д У
12 У 220 2t
\1ж
Уп+Д У ^ - У* Д2 J+1- У* в2) +(т - ^ - У*]/-уи+Д У
т
+ 1Р,
§| т - У +11 Р*1 -1Р*
+
x
2
+
:=2
Приращение деформации в конце второго этапа разгрузки, когда произойдёт проскальзывание второго гофра, может быть найдено из условия
1*11
2 т Р*22 + N 222\ = т
* л
р* + Р21 + Р222 , Р122 ри 2
2
2
с 0 а2 i ^ "а у 2
1" Зтр
\г
= Р21 "1У12 + А У
у,2 "а у 220 в2 " / С
где N222 определяется из (20), Р*222 - из (38) при А~22 = А~ *22 , Р*11 и Р*21 -соответственно по (28) из [2] и (16) при
= ~ V
Приращение можно также найти из условия сшивки процессов разгрузки второго гофра при неподвижных и подвижных концевых сечениях в конце второго этапа разгрузки:
1+^ И ~12+А / - ~ В2 )+[ / - ~ в.
* Л 2
2/
Р + Р
Р11 2
(41)
Способ определения стоящей в пра- Приращение деформации пакета,
вой части (41) силы, действующей на вто- найденное из условия (41) с точностью до рой гофр при проскальзывающих конце- членов второго порядка малости и при вых сечениях, поясним ниже при рас- условии малости величины смотрении последнего этапа разгрузки.
запишется в виде:
у .2+ау
А У 22
Со Р [1 - 3Т' в, - /)]
3тР * Со + 4/ В2 Р21 [1-3Т' в, - /)] С. 2?1 / - ^)2. 2 тр 2 Р+Р))
-Ух=
(42)
Подставив (42) в (40) и учитывая малость величины ()2, найдём си-
12+а~ и)
лу в конце второго этапа разгрузки. Соответствующую деформацию пакета ~*
У 22
найдём, подставив (42) в (36).
Рассмотрим 7-й этап разгрузки (3 < 7 < т). Определим силы, действующие на проскользнувшие гофры. Силу, действующую на первый гофр, найдём по (28) из
[2], подставив вместо ~ значение дефор-
У1
мации пакета на 7-м этапе разгрузки, равное
^2 = ~^ "А2. (43) Сила, действующая на второй гофр,
равна
Со(1„-А~„] (44)
Р2 /2
1"
4/
У^г "А У/2 I Вг " /
Сила, действующая на любой к-й гофр (3 < к < j - 1), равна
р = с 0 а. (~ -а ~/2) (45)
рк?2
-(2к - -а - /
4/ 4 7-1,2 22
Способ определения сил (44) и (45) изложен ниже при рассмотрении последнего этапа разгрузки.
Силы, действующие на (т - 2 + 1) непроскользнуших гофров, определяются аналогично силам, действовавшим на втором этапе на 3-й, 4-й, ..., да-й непросколь-знувшие гофры. Приведём значения этих сил
Рк2=Р4УМ2-АУ
х|С
3<к<т
трР1
(46)
*
*
X
1=1
Сила, действующая на пакет на _]-м этапе разгрузки, запишется в виде:
^ су,-, "ауу
у20
15
а
ищ 8
4+ 2 уу-12 +%
1 т
у-1
1
4
уу-12~АУуМ-/
31- -1)1
V у-12 •'у2
у -АУу -/
(47)
су ~ауУу+{1+к ^ -д+у+у+тр{уу-1,2+%
у-1 * т ( 1 ' *
(т-у +1)!рй +1т-к +1 рй
4=1 к=7ч
+!р,
к=у
N 2=т
3 £ у £ т
Р Р
Р + Р +... + Р* + + - Р
1 11 1 21 1 у-1,1 2 2 у
Ру-
Приу=3 в выражении (47) Е1 = 0. 3-й этап закончится, когда проскользнет у-й гофр.
Приращение деформации в конце у-
го этапа может быть определено из усло- или из условия сшивки процессов раз-
вия преодоления продольной силой Ыу грузки у-го гофра с неподвижными и подвижными концевыми сечениями в конце
сил трения, приложенных к у-й вершине
с 0 ( Уу-1,2-а у *
у-го этапа разгрузки
1 - ^г(2 у- <у *-1,2-а у ^ -/
Р,1 "(. уу-1,2 +а у,2 л с 0 ^ i / - у1 л / - у1 - уу-1,2 +а у,2
2* ^(рр*)
(48)
2г 2
а У, 2 =
Отбросив члены второго порядка малости, из (48) найдём Л Р - У *-1,2 А4]-С0^*.1,2 3 <у < т,
^ <2 у - </+У *-1,2-А У *)+1
• л, + А5
(49)
где
лз =
1 - ^ (2, -1)
У У-1,2 /
л4 = С 0
2 К
1 + -у 1 ][/ - у 1 + Уу-1,2
* *
тр
2г
2 / * n
Р У-1 .
Р1р:
2 =1
л =
5 „,2^0
С 0 Уу-1,2 I/ - У1
2 К у-1-2
1 - М- (2 у -1)1 у ^ - /
4г
у-1,2
тр
4г
(2 у -1) Р*1 - С 0-
Силу, действующую на пакет, и деформацию пакета в конце ]-го этапа разгрузки (3 < ] < т) найдём, подставив (49) в (47) и в (43).
На следующем этапе разгрузки многослойный многопролётный гофрированный пакет деформируется при подвижных концевых сечениях пролётов.
Дифференциальные уравнения, описывающие деформацию каждого гофра пакета, будут отличаться от соответствующих дифференциальных уравнений последнего этапа первичной нагрузки только знаком перед членом, содержащим силы трения. Способы решения этих уравнений и упрощения конечных выражений
■=1
3
для действующих сил аналогичны описанным при рассмотрении последнего этапа первичной нагрузки.
Поэтому ограничимся лишь приведением значений действующих сил.
Сила, действующая на первый гофр, будет равна
р 1.5СД (50)
Р 12
4 + ~'зР * - ~
Сила, действующая на 2-й гофр, запишется в виде:
СА У (51)
Р ="
22
1 -3Ж (~ В. - /
4/ ^
где А. и В2 определяются из соотношений (15) при ] = 2.
Сила, действующая на любой к-й гофр (3 < к < т), будет равна
С„~2 (52)
рк 2 = ■
1 - ^Т (2к -<~ 2 - *
Сила, действующая на пакет на последнем этапе разгрузки, запишется в виде:
^ 2 = С 0~,
1,5
А2
4 +
тр (^
2/ I Зр
У
4/ 1У 2 В2
-/
к=з 1 - ЦР
4/
(2к -1
У2
(53)
и деформация пакета ~2 в сечениях х=//2 может принимать последовательно все значения в интервале 0 < ~2 < _~т2.
Повторная нагрузка. При полной
разгрузке пакета (У2 =0) повторное нагружение его осуществляется по процессу первичного нагружения и зависимость между действующей силой и деформацией пакета при этом нагружении определяется соотношением (19).
При неполной разгрузке (У2 > 0) на первом этапе повторной нагрузки все гофры будут деформироваться без смещения и поворота сечений, лежащих в вершинах. При дальнейшем нагружении, как и в процессе разгрузки, будет происходить проскальзывание вершин гофров, постепенно распространяющееся от крайних гофров к среднему. На последнем этапе повторного нагружения пакет будет деформироваться со смещением всех вершин (за исключением средней) и зависимость между действующей силой и деформацией пакета будет описываться соотношением (19).
Анализ соотношений, полученных при рассмотрении первичного нагружения и разгрузки, позволяет получить все необ-
ходимые соотношение, описывающие любую ьую загрузку пакета. Соотношения, описывающие повторную нагрузку, при
неполной разгрузке (У2 > 0) могут быть получены из этих соотношений как частный случай.
|-ая загрузка пакета. Здесь возможно несколько случаев загрузки пакета.
Случай 1. Первый гофр не проскользнул в ^1-м процессе загрузки.
В этом случае ья загрузка идёт по первому этапу ^1-й загрузки и при достижении действующей силой значения F*i -2 пакет при дальнейшем приращении силы (приращение того же знака) будет загружаться по последнему этапу /-2-й загрузки.
Случай 2. Первый гофр проскользнул в /-1-й загрузке.
На первом этапе /-й загрузки гофры деформируются без смещения вершин и приращение деформации в текущем сечении и деформацию гофров в сечениях х = 1/2 определим из соотношений:
(54)
У,
и
(-1)+
1У1/ (л 2Рх —' i 1 - ^-
2 i /
Уи У1-1
-мл
(55)
+
+
где у- - деформация пакета в конце 1-й загрузки.
Продольная сила, действующая на первый гофр, будет равна
^-н-рс^ун- (56)
Силу, действующую на первый гофр, найдём, решив уравнение
,4 ,2 ; ,2
- N +(-1), т Р^ ^ =
=( Р11 - о ( х - 1
2
ах
в виде
р=л,+у, 1(-1),+1со
1+-
3 (-
1}
(3Р У-
у У11
ах
- N+Ц, = (рш - р^(х -1
ах
ах
2
(59)
где
р / -,\,+1
N 21, = (-1)+1
3С 01 21 2 р п
У1
Ум " 7+У к 2 (1 - 1 С08"
Я*
N 2-1=(-1) С у,-1
р п
8 а р , _ — i сов — к 2 v 3р0 4
у,-1 - 7+У--11 -—1 со§-
(63)
При вычислении R*11 _ 1 из (62) сле-
дует иметь в виду, что
( п* А
Я/к 2 Ч-Г
Р*К +
Р,
2к 2
!к 2 2
(64)
Учитывая (62), (63), (64), можно
R*11 _ 1 записать в следующем виде:
я2-1-(-Г
Л ^ р1к2 2
3С01 ^
(57)
Р 1
41 Р1-11
!х-1)+1у, у 1.1- 7 - .У к2 v1 - ¿0со§ р
(58)
Уравнение, описывающее деформацию второго гофра, запишется в виде:
(65)
р2 Н 4!+1ч i/ -1 ' 2 v 3р) 4
Тогда, решив (59), получим силу, действующую на второй гофр на первом этапе 1-й загрузки, в виде:
р.-рл-ыс (■+щ} [£-v-у^ к
я! р2 , |
f_у,_1 +Ук2с084-(-1) У ( 1х! у р . р
у =---1 1-С08- 1+1) со^вг^,
-Г021 2 V 1 ) Л. 4 4
(60)
у*. _ деформация второго гофра в конце к.-й загрузки (1 < к. < , — 1), на которой последний раз второй гофр проскользнул (принято, что когда второй гофр отслоился на к.-й загрузке, то его деформация в отслоенном состоянии определялась соотношением (9)),
(66)
Определим силу, действующую на /-й гофр (2 < / < т). Дифференциальное уравнение, описывающее деформацию у-го гофра на первом этапе -й загрузки, получим из (59), заменив в нём индекс «2» на «/»:
у1, а2у0/1 * а2у1, (Р р у
ах 2<<т,
(67)
где
^-(^"й у.( У *-1 -/
V -у >(-1гу, (1-
У 0/1 1
Я/ 1-1 = Я/(-2 N/-1
1рх
(68)
(61)
2, - 1 _ продольная сила, действующая на второй гофр в конце 1—1-ой загрузки.
Запишем Я*., - 1 в виде реккурентной зависимости
Ян-1 = Ян-2 - N1 -1'
¡=к+2,к+3,..., (62)
а приращение продольной силы на 1— 1-й загрузке 1 определим из (61) соответствующей заменой индексов в виде:
Якг(-Г
Р*к+ -.+Р*-1,к.+
Р
/к/
/-1к 1
1 < к/ <1 -1 ,
N *1-1=(-1)^ у ц у : - V)
(69)
(70)
т
1
т
Тогда силу, действующую на У-й гофр, найдём в виде:
Р,=;+(-1)+1у1, С
1+-
2 к
- ж. я I 2t Ку'-1{
2 < ] < т .
(71)
Сила, действующая на пакет на первом этапе 1-й загрузки, имеет вид
р=(-Гу С
т+-
3ж 1 - 71 +
У-1- 7 - У:2(1-3Ж ]°084
ж
(72)
+(т-2)| /
_ж
2t
(-1 и р-+§;
-ер;«.
у=1
Первый (или любой у-й) этап закончится либо когда приращение деформации сменит знак (тогда закончится ¡-я загрузка), либо когда проскользнёт первый (или любой ;-й) гофр. Случай смены знака приращения не нуждается в дополнительных пояснениях.
Приращение деформации гофров в конце первого этапа ¡-й загрузки в случае проскальзывания первого гофра найдём, отбрасывая члены второго порядка малости, из условия, аналогичного (33), в виде:
У , = "
1 4
(73)
иж
4t
где
= 1 -(-1)
р;„ -(-1) "с .+-а,
2
иж
4t
У,-1
- /
а, = (-1)'+1 со
3 ( 8
1+13ж У,-1- 7
+иж. р
4t Р1'-1'
Силу, действующую на пакет в конце первого этапа ¡-й загрузки, найдём, подставив (73) в (72).
Случай 3. В ¡-1-й загрузке проскользнуло ; гофров (2 <; < т).
В этом случае к; = ¡-1 и для любого г-го гофра из; гофров, проскользнувших в процессе ¡-1-ой загрузки,
£-1=(-1)м(£ р^
; А
Рп-1
1 < г < ;, 2 < ; < т.
(74)
При г=1 сумма в (74) равна нулю. Для любого 5-го гофра ( < я < т) из т-;,
не проскользнувших в процессе ¡-1-ой загрузки гофров, величина R*si-1 определяется из реккурентного соотношения, аналогичного (62):
Ли = -Е N
г=1
¡=к 2+2, к 2+3-...,
(75)
где
Р
(76)
Рассмотрим первый этап ¡-й загрузки. На этом этапе гофры деформируются без смещения вершин и приращение деформации в текущем сечении и деформация гофров в сечениях х = 1/2 определяются из соотношений (54) и (55).
Приращение продольной силы и сила, действующая на первый гофр, определяются из соотношений (56) и (58). Приращение продольной силы и сила, действующая на второй гофр, определяются соответственно соотношениями (61) и (66), в которых учтено, что кг = г - 1. Приращение продольной силы, действующее на любой из т - 2 остальных гофров, определяется из соотношения (7о), а сила - из (71).
При определении силы, действующей на любой г-й гофр из числа проскользнувших на ¡-1-ой загрузке, в (71) вместо Я*1-1 подставляется Я*п-1 из (74), а при определении силы, действующей на любой я-й гофр (; < я < т) из числа т-;, не проскользнувших в процессе ¡-1-й загрузки, величина Я*,,- берётся из (75).
2
3
2
3
+
2
2
и
+
2
г=1
2
2
2
Сила, действующая на пакет на первом этапе 1-й загрузки, определяется соотношением
Р И-таС 0
2 к
^.У,1 -/1 +
У^ В2 -/] +(т - < у:- /
п-1 я -1
1 п=1 я = ,+1
+: р,
]=1
(77)
Приращение деформации гофров в из (55) при замене индекса «1» на индекс
конце первого этапа ьй загрузки определяется из (73). Силу, действующую на па-
«2».
Силу, действующую на пакет, с учё-
кет в конце первого этапа 1-й загрузки, том (40), (74), (75) и малости у2г по срав найдём, подставив (73) в (77).
Рассмотрим второй этап 1-й загрузки. Деформация пакета Уъ определяется
1С
У,
-1 + (-1)+1~2
4+НГ ^ (| / - 1 -(-1)+1^
нению с У 2, или максимальной деформацией пакета на к-й загрузке найдём в виде:
(78)
^гуь
т-1+-
2 к
в-/) +(т - Чу,- -/
р
: к„-1 +: к,р,-г
Приращение деформации на втором этапе -й загрузки найдем из условия, аналогичного (41), отбросив члены, содержащие (у*1, + Ау*2,)2:
[р-1+(-1)'+1з~1-С0 А2~ *
А У 2, =-
(-1)1 А2 С0 - ^ В2 р'„-1 - 3т4Р 1 ~ В - / 0А9
(79)
где
А9 = С 0
1+2к21 / - У -в
Р- к
2? =2,,-1'
1 -
(-1) ^ 1, В2 - /
Подставив в (78) и (55) (в этом выражении предварительно индекс «1» заменим на индекс «2»), приращение деформации в конце второго этапа, равное
У = У +а У , (80)
У 2, У 1, У 2,
и (79), найдём соответственно силу, действующую на пакет и деформацию пакета в конце второго этапа ,-й загрузки.
Рассмотрим теперь любой к-й этап ьй загрузки (3 < к < J, 2 < j < т). Сила, действующая на пакет, запишется в виде:
3
т +
+
3
]=2
2
2
= С0
Ук
1+(-1)+1 А Ук,
1,5
А2
Ч^Й /-Ук.-(-1)*А?к) 1-(-1)3тР2[(£.Л(-«*Ал)в,-/
,1-(-1»т? е--<1л(-1)"а я-/
+|р>(-1Г(У1-1,+ау.]с (т- к+Т+2к21/ -УУ-^- ~-1+(-1)(ук-"+А
(
р I хк-1 + -1
2 1 г=к я=,+1
+
+
1
Деформацию пакета найдём из соотношения
~.г£,+(-1Г^, (82)
Приращение деформации за к-й этап /-й загрузки найдём из условия, аналогичного (48), пренебрегая членами, содержащими Ду2ы :
0 [ Р 'и-1 + (-1)' + к-1,' ^ ]- С * _и
(-1)/+1 С
0 + ~к-1,'С 0 2^' 1 - ~-1
-(-1)'+1 - ^12+^(2к -1)(-1)+1 р : 3 < к < I, 3 < I < т , (83)
где
2
1 - т!г (2к -1)
Ук-1,'
-/
1 + -
2 И
р =С0~г
'У-г
1" У^ + (-1)~ к
- Р- Р ■
21 г=к ''-1'
1,5
А.
2
=1 - ^ 2 - #к-„+(-1)'
У.-и -1
Подставив (82) и (83) в (81), найдём силу в конце к-го этапа '-й загрузки.
При достижении на j+1-м этапе '-й загрузки (1 < j < т) состояния, когда у]+11=У*1—2, дальнейшее загружение пакета при сохранении знака приращения деформации будет происходить по последнему этапу i-2-И загрузки.
В случае, когда на 1—1-й загрузке проскользнули все гофры ( = т), в соотношениях (77), (78), (81) следует принять R*si — 1 = 0 и подставить в них j = т.
На последнем этапе 1-й загрузки в этом случае пакет деформируется при всех подвижных концевых сечениях пролётов. Сила, действующая на пакет на последнем этапе 1-й загрузки, запишется в виде:
А 2_+ (84)
4 +
(- 1)
3 р
/ -
у
(-1)
3 \ХР
4 /
У В 2 - 1
I = 3
4/ ^ I - ЧУ - 1
_ (-1)
где деформация пакета принимает все значения в интервале
(-гх < (-1)'+1у' < (-1) '+1у;.
Отметим, что при чётных '=2, 4, 6,... соотношения (54) .... (84) будут описывать разгрузочные процессы, а при нечётных — нагрузочные.
Например, при '=1 соотношение (84) совпадает с соотношением (19), описывающем силу, действующую на пакет при первом нагружении.
При '=2, определяя Я*Гг-1 из (74) и учитывая, что ]=т и Р*„-_1=0, из соотношений (54) . (84) найдём ранее полученные соотношения, описывающие первую разгрузку, а при ' = 3 — соотношения, описывающие первую повторную нагрузку.
Библиографический список
1. Кондрашов, Н.С. Упруго-фрикционные характеристики демпферов с гофрированными лентами [Текст] / Н.С. Кондрашов // Вибрационная прочность и надежность двигателей и систем летательных аппаратов. — Вып. 51. Куйбышев: КуАИ, 1972. — С.45-62.
2. Эскин, И.Д. Циклическое сжатие многослойного однопролетного пакета. [Текст] / И.Д. Эскин, Р.И. Алкеев, В.И. Иващенко // Вестник Самар. гос. аэрокос. ун-та. — 2012. — Вып. 3 (34), Ч. 2. — С.97-106.
3. Эскин, И.Д. Упрощенная модель циклического сжатия многослойного многопролетного гофрированного пакета. Экспериментальное и расчетное исследование. [Текст] / И.Д. Эскин, Р.И. Алкеев, В.И. Иващенко // Вестник Самар. гос. аэрокос. ун-та. — 2013. — №1.
А10
10
3
А" = С0
+
8
CYCLIC COM PRESSION OF A MULTILAYER MULTISPAN LAMINATED PACKAGE OF CORRUGATED SHEETS
©2013 I. D. Eskin, R. I. Alkeev, V. I. Ivashchenko
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
The task on cyclic compression of a multilayer multispan package of corrugated sheets is solved by the Bubnov-Galerkin method taking into account the rotation of the package sections located at the vertices of corrugations. The relations describing any loading processes in the field of elastic - hysteresis loops of the package are obtained.
Multilayer multispan package of corrugated sheets, cyclic compression, Bubnov-Galerkin method.
Информация об авторах
Эскин Изольд Давидович, кандидат технических наук, старший научный сотрудник, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: eski-nid [email protected]. Область научных интересов: конструкционное демпфирование, динамика и прочность.
Алкеев Руслан Иванович, старший преподаватель кафедры инженерной графики, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: конструкционное демпфирование, динамика и прочность.
Иващенко Владимир Иванович, кандидат технических наук, доцент, заведующий кафедрой инженерной графики, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: конструкционное демпфирование, динамика и прочность.
Eskin Izold Davidovich, candidate of technical science, senior researcher, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: eskinid [email protected]. Area of research: structural damping, dynamics and strength.
Alkeev Ruslan Ivanovich, assistant of the engineering drawing department, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: [email protected]. Area of research: structural damping, dynamics and strength.
Ivashchenko Vladimir Ivanovich, candidate of technical science; associate professor, head of the engineering drawing department, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: [email protected]. Area of research: structural damping, dynamics and strength.