УДК 630*378 С.В. Посыпанов
Северный (Арктический) федеральный университет
Посыпанов Сергей Валентинович родился в 1961 г., окончил в 1983 г. Архангельский лесотехнический институт, кандидат технических наук, доцент кафедры водного транспорта леса и гидравлики Северного (Арктического) федерального университета. Имеет около 30 печатных работ в области лесотранспортного использования малых и средних рек, информационного обеспечения лесопромышленного производства. Тел.: 8(8182) 21-61-50
КОМБИНИРОВАННЫЙ МЕТОД РАСЧЕТА ПАКЕТА КРУГЛЫХ ЛЕСОМАТЕРИАЛОВ, УЛОЖЕННОГО НА ГОРИЗОНТАЛЬНОМ ОСНОВАНИИ
На основе эластиковой теории получены аналитические зависимости для определения натяжения в обвязках пакета круглых лесоматериаов, уложенного на горизонтальном основании.
Ключевые слова: лесосплав, пакет, пучок, обвязка, эластика.
Для транспортировки древесины по водным путям без потерь, минимизации воздействия на водную среду весьма существенна правильная оценка прочностных характеристик лесо-транспортных единиц, в частности пакетов круглых лесоматериалов или лесосплавных пучков. В принятой ныне трактовке под прочностью лесосплавного пучка понимается способность обвязок обеспечивать сохранение его формы под действием сил распора бревен. Она характеризуется прочностью обвязок [4].
Наибольшую нагрузку в состоянии покоя обвязки испытывают, когда пакет лесоматериалов уложен на твердое основание, например на плотбище. Этот случай будем рассматривать как расчетный.
Согласно эластиковой теории [2] совокупность комплектующих пакет бревен уподобляется некой несжимаемой жидкости, а сам пакет - части матерчатой невесомой бесконечно длинной цилиндрической оболочки, наполненной под давлением этой жидкостью.
Поперечный профиль пакета бревен, лежащего на плоской горизонтальной поверхности (рис. 1), в данном случае может быть описан параметрическими уравнениями в осях хОу:
* = \НГ> (!)
У = ЯД, (2)
где у - ордината, расстояние по вертикали от любой точки кривой профиля до уровня жидкости в пьезометрической трубке на расчетной схеме; х - абсцисса; Н - характерная высота.
Рис. 1. Расчетная схема пакета бревен, лежащего на плоской горизонтальной поверхности
/ = /(0, ф) = 2|(2 - sm2 0) ^(0) - F(0, ф) _- 2 f (0) - Е(®, <р)
А = А(0, ф) - л/l" sin2 © sin2q>; (4) 0 - модулярный угол;
(3)
К (0), Е (0) - полные эллиптические интегралы; Т7 (0, ф), Е (0, ф) - неполные эллиптические интегралы; ф - угол, изменяющийся вдоль кривой так, как отмечено на рис. 1.
Для заданной кривой, называемой бесперегибной эластикой, характерная высота Я и модулярный угол 0 постоянны. Таким образом, в рассматриваемом случае параметром является угол ф.
Из уравнений (1) и (2) можно получить выражения для определения основных геометрических характери-
стик пакета. Ограничимся лишь теми, которые понадобятся нам в данном случае.
Высота пакета
я =Н -И = (1-соъ®)Н = 2Н ът2 —. (5)
2
Ширина пакета
В =2х 1
в
^Я/(0,45°). (6)
Произведение ординаты точки линии (у) и ее кривизны (р)
yp = -H2sm2® = const. (7) 4
Площадь поперечного сечения пакета
0 = Я2[(2- sin2 ®)К-2Е]. (8)
На основании формул (5) - (7) запишем коэффициент формы пакета:
С =
Яп 21
2 ©
sin ©)|Х - F(©,45 )] - 2[Е - Е(&,45 )]cosecz
(9)
Как видно из выражения (9), коэффициент С является только функцией от модулярного угла 0. Между этими величинами существует однозначная зависимость. Можно также сказать, что угол 0 является функцией от С. Каждая из этих величин однозначно определяет форму эластики, а характерная высота Н - размеры кривой.
Недостаток выражений, полученных благодаря эластиковой теории, заключается в том, что они не позво-
ляют получить в явном виде зависимость одних общеупотребимых параметров пакета от других. Приходится задаваться различными значениями 0 и ф и находить соответствующие параметры, например высоту и ширину пакета.
В данном случае представляет
интерес безразмерная величина
*
р Я .
п
В соответствии с выражением (5)
h = Н р (1-C0S 0)Я 1-cos©
Эта величина показывает относительную высоту или толщину слоя «пучковой жидкости» в пьезометре, во сколько раз она больше высоты пучка. Согласно выражению (10) кр зависит от модулярного угла 0 или коэффици-
H
2H sin
2
2sin2
0
(10)
2 2 ента формы С, т. е. только от формы пучка.
Для практических целей зависимость кр от 0 удобнее заменить зависимостью кр от С, что несложно сделать, имея в виду выражение (9). Зави-
1
симость кр=у(С) можно представить в табличном виде (табл. 1).
Таблица 1
вторая - относительное полуотверстие эластики:
©,град С hp ©,град С hp
45 1,16 3,41 75 1,75 1,35
60 1,33 2,00 80 2,04 1,21
65 1,43 1,73 85 2,59 1,10
70 1,56 1,52 89 4,04 1,02
Ъп
1
[(2 - sin2 &)К - 2Е]. (12)
Для речных пакетов, т. е. при С = 2.. .3 величина Н превышает высоту пакета Нп не более чем на 21 %. При этом зависимость hp от С в указанном диапазоне близка к линейной. Для промежуточных С значения hp можно с достаточной точностью получить в результате интерполяции.
В связи с тем, что при сплотке пучков на берегу зачастую приходится в качестве опоры использовать соседний готовый пучок, интерес могут представлять такие безразмерные величины, как с/Н и n/H (рис. 1). Первая из них - это относительная высота наиболее выступающей боковой точки пучка:
п
Н 2'
Учитывая, что верхняя часть эластики с обеспечением достаточной точности может быть заменена полуэллипсом, а нижняя двумя четвертями эллипсов [2], можно утверждать, что величины кЕ, Ьп в совокупности с С дают достаточно полное представление о форме эластики. Размеры ее определяются величиной Н.
Для полноты картины рассмотрим еще одну безразмерную величину - относительную площадь поперечного сечения пакета юН. С учетом выражения (8)
юя = —- = (2-sin2 ©Ж-2^. (13)
Я
к = — = 1 — • Е Н
2-
sin9 ©
2
(11)
Согласно формулам (11) - (13) рассматриваемые безразмерные величины зависят только от модулярного угла 0. Для практики лесосплава их лучше заменить зависимостями от коэффициента формы С, как мы это уже делали для кр. Они представлены в табличной форме (табл. 2).
Таблица 2
©,град С hE Ьп с/Нп n/Нп юд
45 1,16 0,134 0,040 0,457 0,136 0,080
60 1,33 0,209 0,137 0,419 0,273 0,273
65 1,43 0,232 0,197 0,402 0,341 0,394
70 1,56 0,253 0,280 0,384 0,426 0,561
75 1,75 0,270 0,400 0,364 0,540 0,801
80 2,04 0,282 0,584 0,342 0,707 1,168
85 2,59 0,290 0,918 0,318 1,005 1,835
89 4,04 0,293 1,717 0,298 1,748 3,435
Для речных пучков при С = 2.. .3 величина кЕ меняется очень незначительно. Высота наиболее выступающей боковой точки пучка составляет чуть менее 30 % от Н, которая, в свою очередь, в указанном диапазоне С немного превышает высоту самого пучка. Для озерных пучков при С <1,5 разброс кЕ хотя и несколько больше,
но тоже невелик. Здесь Н существенно больше высоты пучка. Полагаем, что интерес представляют также отношения с и п к высоте пучка Нп, которые также приведены в табл. 2.
Зависимости относительного полуотверстия эластики от коэффициента формы пакета близки к линейным. Для этого случая мы подобрали
эмпирические выражения (при С = = 1,16.. .4,04)
Ъп = 0,586С - 0,630; (14)
— = 0,550С-0,446. (15) И
п
Зависимость относительной площади от коэффициента формы пакета также близка к линейной. Подобрано эмпирическое выражение (при С = = 1,16.4,04)
юя = 1,172С - 1,261. (16)
Геометрические параметры пакетов, получаемые с помощью эласти-ковой теории, достаточно точны, чего нельзя сказать о силовых характеристиках [6]. Дело в том, что эластиковая теория не учитывает силы внутреннего трения в пакете, трение обвязок о бревна. Кроме того, не следует забывать об уже отмеченных неудобствах, связанных с невозможностью в явном виде выразить одни общеупотребимые параметры пакета через другие. В связи с этим, приняв во внимание полученную зависимость кр=у(С), рассмотрим силовые характеристики пакета, использовав другой подход к решению задачи.
Расчетная схема сил, действующих на пакет, уложенный на горизонтальном основании, представлена на рис. 2.
Совокупность бревен пакета будем считать сыпучей средой, которая сжимается под воздействием сил тяжести (равнодействующая О) и реакции основания (равнодействующая Я). В результате вертикального сжатия сыпучей среды возникают горизонтальные распорные усилия. Вертикальное давление р в определенной точке наиболее высокого продольного сечения пучка может быть определено по формуле
Р = 1дЩко,
(17)
где уд - объемная масса бревен в пакете; г| - коэффициент полнодревесно-сти пакета;
g - ускорение свободного падения; кО - расстояния по вертикали от точки О (см. рис. 2).
/1
\
Т/
7777777777777^777777777777
Вп
77
Т77777777777Т 6
Рис. 2. Расчетная схема сил, действующих на пакет: а, б - схемы соответственно внешних и внутренних сил
Горизонтальное давление q (см. эпюру на рис. 2) пропорционально
вертикальному р\
q = pi а, (18)
где А-а - коэффициент активного давления.
Горизонтальное давление в верхней точке пакета
ch = la4SXAH-Hu)- (19)
в нижней
= (20)
Сила горизонтальных давлений для рассматриваемого сечения (равнодействующая распорных усилий)
^ = 0,5(д1+д2)Яп4 (21)
где L - длина бревен в пакете.
Это усилие уравновешивается реакциями в верхних и нижних ветвях обвязок Яв и Ян.
Сумма моментов относительно точки А:
1 2
ян -РН *-= 0; (22)
н п п С1Х + С12 относительно точки К:
¥Н
2 1
3-3--кн =0 (23)
в п
я,+я2
Выразим Ян и Яв из формул (22), (23), вместо , д2 и ^ подставим правые части выражений (19) - (21) и после преобразований получим
Дн=0,5удЛ^а1Яп(Я-^Яп); (24)
2
Дв=0,5удт^а1Яп(Я--Яп). (25)
Измерить высоту Н весьма проблематично. Вспомним введенную нами относительную высоту кр, равную отношению Н к Нп. С учетом этого отношения после несложных преобразований перепишем выражения (24) и
(25):
9 1
Ян=0,5у д^КШ;(Ир--у, (26)
2
Яв=0,5удЛ^а1Яп2(/7р--). (27)
Определить величину кр можно ранее рекомендованным способом.
Коэффициент активного давления для сыпучей среды может быть установлен из выражения
= (28)
Натяжения распределяются на обе обвязки пакета. Подбор сечений обвязок при известном их натяжении отражен в соответствующей нормативно-справочной литературе, не вызывает вопросов и поэтому здесь не рассматривается. Как видно из формул
(26), (27), большее натяжение у пакета, лежащего на твердом основании, возникает в нижних ветвях обвязок.
В настоящее время для определения натяжения в обвязках лесосплавного пучка, размещенного на горизонтальном основании, в нормативно-справочной и другой специальной литературе [1, 3, 5, 6] рекомендуется эмпирическая формула
Н = 0,15 пщ. (29)
С-0,83
В упомянутых литературных источниках не уточняется, для какой части обвязки рекомендуется эта формула.
Для сопоставления эмпирической и аналитических формул преоб-р азуем выражения (26) и (27). Учитывая, что масса пакета может быть оп-р еделена по формуле
т- — В Н Ьг\у ,
4 П п д
а ширина пакета
Вп = СНШ
запишем выражения (26) и (27) в виде
2 к(К-Ь
-1Щ-,
К =
3
Сп
я
С.л
.
Теперь эмпирическая и аналитические формулы отличаются выражениями, расположенными перед произведением Назовем их удельными натяжениями обвязок - натяжениями на единицу веса пакета. Сравним эти величины, представив зависимости удельных натяжений обвязок от коэффициента формы пакета в графическом виде (рис. 3).
По графикам видно, что удельные натяжения в нижних ветвях обвязок, полученные аналитическим путем и по эмпирической формуле при С > 1,25, достаточно хорошо совпадают. При меньших значениях С результаты отличаются существенно, что говорит об
g 0,7
со
S о-»-'-'-'-'-'
>» 1,0 1,5 2,0 2,5 3,0 3,5 4,0
Коэффициент формы пакета
Рис. 3. Графики зависимостей удельного натяжения обвязок от коэффициента формы пакета: 1 - данные, полученные эмпирически; 2, 3 - то же аналитически для нижних и верхних ветвей обвязок
ограниченности применения эмпирической формулы. Эта формула только выглядит более простой, расчеты по ней не проще, так как не только С, но обычно и m устанавливают вычислительным путем. Кроме того, для дальнейшего анализа прочности удержания бревен в пакете необходимо знать натяжение не только в нижних, но ив верхних ветвях обвязок.
Преимущество полученных аналитических формул перед эмпирическими состоит и в том, что они нагляднее отражают физическую сущность явления, показывают зависимость натяжений в обвязках от изначально определяющих их величин. Натяжения в обвязках пропорциональны длине бревен, квадрату высоты пакета, его полнодревесности, объемной массе бревен и ускорению свободного падения. Они уменьшаются с увеличением угла внутреннего трения бревен. Зависимость от коэффициента формы пакета сказывается в выражениях (26) и (27) через относительную высоту к1,.
Для речных пучков с коэффициентом формы 2.3 эта зависимость хотя и довольно значима, но, вероятно, в меньшей степени, чем принято считать. Так, при увеличении коэффициента формы от 2 до 3 при прочих равных условиях натяжение в нижних ветвях обвязок уменьшится на 18, в верхних - на 36 %.
СПИСОК ЛИТЕРАТУРЫ
1. Водный транспорт леса: справочник / под ред. В.И. Патякина. М.: Лесн. пром-сть, 1973. 404 с.
2. Воробьев А.Г. О расчете по эластиковой теории пучковых плотов для случая нахождения их на суше // Лесоинж. дело. 1958. № 3. С. 40. (Науч. докл. высш. шк.).
3. Инструкция по эксплуатации такелажа на лесосплаве. М.: ВНИПИЭИЛеспром, 1980. 134 с.
4. Плоты (конструкция, эксплуатация, технология) / М.Н. Фоминцев [и др.]; под ред. М.Н. Фоминцева. М.: Лесн. пром-сть, 1978. 216 с.
5. Правила (технические условия) сплотки, формирования и оснастки плотов для буксировки в Северодвинском бассейне. Архангельск, 1989. 65 с.
6. Реутов Ю.М. Расчеты пучков (пакетов) круглых лесоматериалов. М.: Лесн. пром-сть, 1975.152 с.
Поступила 24.03.09
S.V. Posypanov
Northern (Arctic) Federal University
Combined Calculating Method of Round Timber Pile Stacked on Horizontal Foundation
The analytical dependencies are received based on the elastic theory for determination of tension in binding of round timber piles stacked on horizontal foundation.
Keywords: timber rafting, pile, bundle, binding , elastics.