2016
октябрь-декабрь
№ 4 (33)
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
УДК 539.3
ТЕРМОНАПРЯЖЕННОЕ СОСТОЯНИЕ МНОГОСЛОЙНЫХ ПОЛИАРМИРОВАННЫХ ГЕЛИКОИДАЛЬНЫХ ОБОЛОЧЕК
2
Ю.В. Немировский1, А.И. Бабин2, Е.А. Сальский3
1 Институт теоретической и прикладной механики им. С.А. Христиановича Кузбасский государственный технический университет им. Т. Ф. Горбачева 3Кемеровский государственный университет
3
3
В общем случае оболочечная конструкция геликоидальной формы может состоять из т армированных слоев, различающихся составами связующего, материалами армирующих волокон и схемами армирования. Получены физические составляющие тензоров эффективных тангенциальных жесткостей и температурных напряжений для многослойного полиармированного композитного материала в системе координат, не связанной с микроструктурой материала.
При этом физические составляющие тензоров эффективных тангенциальных жестко-стей и температурных напряжений армированного материала зависят от модулей Юнга материалов арматуры и связующего, коэффициентов Пуассона арматуры и связующего, интенсивностей армирования и направления укладки волокон. Используемая структурная модель армированного материала позволяет после решения краевой задачи и определения осредненных напряжений находить напряжения в элементах композиции - арматуре и связующем материале.
С целью определения нелинейного закона распределения температуры по толщине многослойной полиармированной оболочки отрицательной гауссовой кривизны найдено аналитическое решение контактной задачи теплопроводности. Приведена замкнутая система уравнений статики многослойных геликоидальных оболочек с учетом внешних термосиловых полей в системе координат нормально связанной с отсчетной поверхностью. Порядок системы уравнений не зависит от количества слоев и схем армирования. Разработана методика расчета армированных оболочек отрицательной гауссовой кривизны, которая учитывает структурную неоднородность материала, а также возможность за счет варьирования внутренней структуры достигать наибольшей согласованности жесткостных полей конструкции с полями напряжений от внешних термосиловых воздействий.
Ключевые слова: полиармированная оболочка, прямой геликоид, температура, не связанная задача термоупругости, обобщенные кинематические гипотезы Тимошенко, закон Дюамеля-Неймана.
БОТ: 10.17212/1727-2769-2016-4-7-21
В современной строительной практике широкое применение находят линейчатые оболочки отрицательной гауссовой кривизны. В конструкциях современных многоэтажных гаражей используются сходни - пандусы, выполненные в виде прямой геликоидальной оболочки [1]. Элементы в виде геликоидальных оболочек могут найти широкое использование при оформлении перекрытий новых архитектурных ансамблей.
Исследование выполнено при финансовой поддержке гранта РФФИ № 15-01-00825.
Введение
© 2016 Ю.В. Немировский, А.И. Бабин, Е.А. Сальский
Одним из основных факторов риска для таких конструкций является возможность возникновения пожара в здании. Температура окружающей среды при пожаре может достигать 1000 °С и более в течение всего лишь десятков минут [2].
В этом случае тонкостенные элементы работают в условиях высокотемпературного неравномерного нестационарного нагрева, приводящего к возникновению тепловых напряжений, знание величины и характера которых необходимо для всестороннего анализа несущей способности конструкции. Тепловые напряжения как сами по себе, так и в сочетании с механическими напряжениями от внешних сил могут вызвать появление трещин в элементах конструкций из хрупких материалов, возникновение и развитие пластических деформаций, ведущих к полному или прогрессирующему разрушению композитных тонкостенных конструкций, их термовыпучивание.
В настоящее время геликоидные турбины разработчики пытаются приспособить под ветровые преобразователи энергии. Геликоидные турбины (Патент Австрии № 117749 от 10.05.1930) с винтовыми направляющими, которые закручивают поток воды в продольный вихрь и способствуют возникновению гироскопического эффекта, оказывающего действие на тело ротора турбины. Известны также геликоидные турбины с вертикальной осью вращения и лопастями, изогнутыми по винтовой линии (Патент США № 5451137).
С целью облегчения таких конструкций при их изготовлении можно использовать современные волокнистые композиты: боро-, органо-, углепластики.
Отметим специфические особенности оболочек из композитных материалов: резко выраженную анизотропию их деформативных свойств, ослабленное сопротивление поперечным деформациям, существенное различие механических и теп-лофизических характеристик слоев. Эти факторы имеют принципиальное значение [3-10] при определении полей напряжений тонкостенных слоистых элементов конструкций.
Строительная практика промышленно развитых государств подтверждает перспективность применения подобных материалов. Однако отсутствие надежных методик расчета армированных оболочек отрицательной гауссовой кривизны, которые бы учитывали структурную неоднородность, материала, а также возможность за счет варьирования внутренней структуры достигать наибольшей согласованности жесткостных полей конструкции с полями напряжений от внешних воздействий сдерживает проектирование подобных конструкций. Даже расчет изотропных геликоидальных оболочек строится на различных упрощающих предположениях, которые часто противоречат положениям теории оболочек и реальным условиям работы конструкции.
К примеру, весьма распространен подход, основанный на применении безмо-ментной теории оболочек [11, 12]. Однако авторы не замечают, что поскольку границами геликоидальной оболочки являются асимптотические линии, то это противоречит условиям применимости безмоментной теории [13] из-за возможности возникновения обобщенного краевого эффекта. С другой стороны, безмо-ментная краевая задача теории оболочек является статически определимой. Значит, для вычисления усилий и напряжений не используется физический закон деформирования, и, следовательно, усилия и напряжения никак не зависят от структуры и свойств материала оболочки. Все это делает неприемлемым использование безмоментной теории для анализа деформирования армированных оболочек отрицательной гауссовой кривизны, напряженное состояние которых может качественно изменяться за счет варьирования структуры материала оболочки [14].
В известных работах, посвященных расчету изотропных геликоидальных оболочек, на основе моментной теории также используются не всегда оправданные упрощающие предположения. Например, в [15, 16] авторы пренебрегают ради-
альной компонентой вектора смещений, что приводит к понижению порядка разрешающей системы дифференциальных уравнений. При этом не дается оценки погрешности данных допущений. Часто делается предположение о том, что оболочка имеет большое число витков, нагрузка и шаг геликоида не изменяются в окружном направлении и влиянием закрепления прямолинейных кромок можно пренебречь [25]. Эти предположения приводят к тому, что компоненты напряженно-деформированного состояния не зависят от окружной координаты и краевая задача описывается обыкновенными дифференциальными уравнениями. Подобный подход справедлив для оболочек с большим числом витков (больше десяти) [1], что в строительной практике встречается нечасто.
Экспериментальные исследования показали, что разрушения от поперечного сдвига часто лимитируют несущую способность армированных конструкций [17]. Учет обжатия влияет на напряженно-деформированное состояние конструкции значительно меньше, чем учет поперечного сдвига.
В связи с вышесказанным при выводе дифференциальных уравнений, описывающих напряженно-деформированное состояние геликоидальной оболочки, будем использовать обобщенные кинематические гипотезы Тимошенко [18], позволяющие учесть поперечные сдвиги, удовлетворяющие краевым условиям на лицевых поверхностях оболочки.
Постановка задачи
Рассмотрим прямую геликоидальную поверхность (см. рисунок), задаваемую радиусом-вектором:
г (г, Ф) = гёр+1 фёг, I = -2-,
где г, ф - криволинейная ортогональная система координат; t - шаг геликоида.
Компоненты основного метрического тензора, тензора кривизны, а также главные кривизны прямого геликоида определяются следующим образом [19, с. 291]:
2 2
йц = 1, ^12 = 0, ^22 = г +1 ;
Ьп = ь22 = ° Ь12 =-;
^Г 2 +12
Л11 Л22
= 0,
Л
12
г 2 +12
22
= -Ь22, К =--
(Г 2 +12)2
(1)
где а и g - дискриминанты метрических тензоров поверхности и пространства
соответственно; К - гауссова кривизна. Параметры Ламе имеют вид
А = 1; А2
= >/ г 2 +12
(2)
В общем случае оболочечная конструкция может состоять из т армированных слоев, различающихся составами связующего, материалами армирующих волокон и схемами армирования.
Поверхности раздела к-го и (к + 1)-го слоя описываются уравнениями:
к т
2 = 2к+1. 2к+1 = Ек], Ек] =2Н (к =1'т-1)-] =1 ]=1
Будем рассматривать оболочку достаточно тонкую, поэтому теплообменом с поверхностей (см. рисунок) г = Г(, г = Г2; ф = 0, ф = Дф можно пренебречь.
На поверхностях 2 = -Н , 2 = Н оболочки имеет место конвективный теплообмен (при граничных условиях первого, второго или третьего рода):
дТп
Л(33) -¿Г + а0,1 (Т(1) - Т(1) ) + 90,1 = 0, 2 = -Н;
д2
(3)
дТ(
Л(33) дт +а т, т+1 (Т(т+1) Т(т)) + 9т, т+1 = 0, 2 = Н ■
(4)
В (3), (4) температуры внешних сред и коэффициенты теплообмена обозначены через Т(0), Т(т+ц и а0,1, ат,т+1 соответственно; удельные мощности поверхностных источников: =-Н = 90,1, ч\2=Н = 9т,т+1 ■
На границах раздела армированных слоев выполняются условия идеального теплового контакта:
Т(к )12=2к+1 -= Т(к+1)
2=2к+1+ '
Л(к) ^ Л(33) д2
2=2к+1"
= Л(к+1) дТ(к+1) = Л(33) д2
(к = 1,2,-,т -1).
(5)
2=2к+1 +
В качестве первого приближения для неизвестных функций Т(к) примем нелинейный закон распределения температуры по толщине пакета, который следует из решения одномерных уравнений стационарной контактной задачи теплопроводности для оболочки, состоящей из т армированных слоев [20],
( |£л33 ^ ^
.а &
V У
а^к)=0 (к=1т),
(6)
где а и g - дискриминанты метрических тензоров поверхности и пространства соответственно; w(k) - удельная мощность объемных источников диссипативного характера.
Интегрируя уравнение (6), получаем 2т параметрическое семейство первообразных функций:
Т(к) (2) = А(к) + В(к)Э(2) - жТк) (2), к = 1
(7)
Ж2) =
24Ё
1п
К2 - н-4Е
К2 - н+4Е
ж
(к)
(2)=(л(33) Г1
а ^) ^
Ч g ^ ^
а д
а
а д.
Здесь К - гауссова кривизна; Н - средняя кривизна и Е - эйлерова разность [19].
После подстановки выражений (7) в граничные условия (3)-(5) однозначно определяются 2т коэффициентов А^), Б^) и находятся конечные выражения
для искомых функций:
Т{к) (г ) = Ф(к) -Sk_2 +
'к?)ф-(л(33) )-11'
^ I (зз)/ £1 3'(2к+1)
}-
<[(г) -9(Гк)]- ^ (г),
Оз
Ф(к) = & + Ф[_1 _ & (21)
к = 1, т
5 = 1 {(л(3+)1) )-1 [в( г, +1)-Э( г,)] ± ^ } + (^+1);
Ф = Р2 -1 + Р
9' (Гт+1) (Л(3тз)) )-1 £ ^^у - ^т (Гт+1)
(8)
3 =1 К« [9(г, +1) -9(г, )]; ] =1
Ф =
Рз - РОзО-1 +Ф
РК^Э'(7т+1) + Р2 [^т Ш-1 ] '
Радиальная, окружная и нормальная компоненты вектора смещений оболочки представляются в виде [10]
з
+Х г У1г(Ф), ™ = ^ Ф), (9)
г=1
иг (ф) = и0г (ф)
где иог, иоф - радиальная и окружная компоненты вектора смещений срединной поверхности оболочки; угг, угф - коэффициенты разложения в степенной ряд
смещений оболочки - функции координат срединной поверхности.
В дальнейшем предполагается, что оболочка находится под действием нормальной нагрузки д(г, ф). В этом случае на лицевых поверхностях г = ±И должны выполняться условия СТГ(ф)Г = 0 .
Компоненты тензора деформации выражаются через обобщенные смещения:
ди.
ег + (г3 - 3И2г)^^ - г
2 ) дузг д ю ;
дг
дг дг2
1 ди0ф 1 йА2 з 2 ч
еФ =--- +--- и 0г + (г - зИ 2)-,
Ф
А2 дф А2 йг
1 дУзф 1 йА2
---+---У зг
А2 дф А2 йг
- г
1 йА2 дю 1 д2ю 2 -+■
А2 йг дг А| дф2
I д«0ф 1 д«0г 1 А о , ^2 ч
дг а дф а -г
^3ф+_ 1 дуэг
дг А2 дф
1 -А2 а2 -г
"У 3ф
2 д2ю 1 -А2 -А2 -ю
А2 дг дф а2 -г -г - ф
- 2 Ьи ю| А2 |
(10)
^ = |(22 -Н2)тэг; V = |(22 -Н2)У3ф.
Дифференциальные уравнения равновесия и набор краевых условий на прямолинейных и винтовых кромках оболочки были получены с помощью принципа возможных перемещений.
Дифференциальные уравнения равновесия:
^(А2Тг) +—-НА* Тф= 0; А(А25) + НА* 5 +^ = 0; дг дф -г дг -г дф
(А2 Ыг ) + 2-
дг 2
д 2 М
гф
1 -А2 дм.
гф
1 д2 Мф д ( -А,
дг дф А2 -г дф А2 дф2 дг V -г
м ф|+
2А2
я
-5 = -А2 9;
12
(11)
д дН гф аА2
—(А2 Нг) +—-2 Нф- А2бг = 0; дг дф йг
аА2
дНф
— (А2 Нгф) + ^ Нгф +-А20ф = 0.
дг 2 гф йг гф дф 2 ф
Краевые условия на направляющих и прямолинейных кромках:
(Тф( г) - Тф( г) )(Мф(г ) - м0ф( г)
Г0; (5-Тг0(ф))(
5 - Тг(ф) ) (иф(г) - м0ф(г) ) = 0;
))
мф( г) - м
0
ф( г) ф(г)
• = 0; ((Нф(г ) Нг(ф)
дф(г)
)(у3ф(г) у03ф(г))
= 0;
( дм 2
гф + -А^м + 1 дмф - д0
дг А2 йг гф А2 дф
д дмгф аА2 0
(А2 мг) +—--г1 мф- А20° дг дф -г
(^ - V0) = 0 (при ф = 0, Дф); (12)
(V - м> ) = 0 (при г = г1, г2);
((гф - (ф) -у03г(ф)) =
Входящие в (11), (12) усилия и моменты определяются следующим образом:
Н Н
Тф(г) = | СТг(ф)-2; мг(ф) = | стг(ф)т-т; -Н -Н
1Н
— | /1(2К(ф)-2;
Н
г (ф)
Н
-Н
н н л Н
5 = I стгф-2; мгф= | стгфm-m; Нгф= -у | /1(2)стгф-2;
-Н -Н Н -Н (13)
Н
Qг(ф) = | У2 С2К(ф)-2
-Н
Внешние усилия и моменты, приложенные к краям оболочки, определяются по формулам:
Н 1 Н
Тф(г) = | Рг(ф)-2; Нг(ф) = Н2 1 /1(2)Рг(ф)-2;
(14)
-Н Н -Н
Н Н
м0(ф) = I 2Рг(ф)-Г; б0 =| Рz-т,
-Н -Н
где /1(2) = 2 - 3Н 2; /2(2) ; Рг, Рф, Р2 - компоненты распределенной
,3 2„. (= — -2
нагрузки, приложенной к границам оболочки.
Компоненты тензоров напряжений, деформаций и температура для к -го слоя связаны зависимостями [21]:
„(к) = А(к) Р(к) - с(к) Т(к)
°(ар) _Л;арХц)ь(Хц) (ар)
(15)
(к) = -(к) (к) (к) = ~(к) (к) Т(а3) = Р(ар)У(р3), У(а3) = 9(ар)Т(р3).
Здесь ||р(ар)|| = ||с(ар)| 1.
Физические составляющие А(арр), р тензоров эффективных тангенциальных упругих и температурных жесткостей А и с, тензора поперечных сдвиговых податливостей q для к-го армированного слоя определяются формулами [21]:
(к) (■ -юк 1)Е(к 1
1(ар^ц) 1 - v(k)2 Х
1 Ус
+ Ю(к ) Е (к ) (16) + юz Е(арХц), (16)
А
vc 2 ((^рц + 8ац8рХ)
((а р=ю2к) к;а р)+(1 -ю(2к) )хск) §ар,
2 (1 -ю2к) ))1 + vCk)),
9(ар) = ЕсСк) Ьар+юг А (ар), (17)
-(к) ) ^ +ю(к)г(к)
Е (к "с
10, если а Ф р,
§ар = и р
11, если а = р;
Е
Е(к) = -(к) + (1 --(к)) Ек + (1111) а С
) ЕСк) [-(к) уак) +(1 --(к) )Ск) ]2
-(к) (1 -у(к)2)к) +(1 --(к))(-V?)2))к) ,
(к) Еак)Е?) [-(к)уак) +(1 --(к))у»)] (1122) та(к) (1 -у(ак)2)+(1 --(к))(1 -уСк)2)
Е
Е
(к) 1_Е0М1_
(1212) 2 та(к) (1+уак))+(1 -та(к))))
(к) Е{к) Е (к)
(к) __дс да_
(2222) -(к) (1 -у(ак)2)к) +(1 --(к))(1 -у«2)
Е(к) = 0, Е(к) = 0;
г(к) = Г(11)
(1122) (1222) 2 (1+уСк) )(1+уак))
з(к) (1 + уСк))Е(к}+(1 --(к))( + у(к)))к)
г(к) 2 ^) (1 + у(ак)))) +(1 ))(1 + у»))#)
г(22) 2 Е (к) Е (к) ,
С ^а
к((й = [-(к )у(ак) +(1 -®(к) )уСк) *
-(к)ЕСк) (1 -у(к)2)к) +(1 --(к))Еак) (1 -у(ак)2)
к
(к) = (22)
(18)
(19)
(20)
-(к)Е(к) (1 - у(ак)2) + (1 - -(к)) Еак) (1 - уСк)2)
+-(к) (1 -у(ак))) +(1 -(к))(1 -уСк))хСк), (21)
-(к)ЕСк) (1 -у(ак)2)) + (1 --(к))Е(к) (1 -у(к)2)хСк)
-(к)Е(к) (1 - у(ак)2)+(1 - -(к))Е(к) (1 - у(к)2) ,
к (к) =к (к) =0 (12) (21) и'
В соотношениях (15)-(21) приняты следующие обозначения: ЕПк), у(П),
(к) Е(к)а(к)
у}„ ) = ——(к)--модули Юнга, коэффициенты Пуассона и температурного рас-
1 -уп
ширения для арматуры (п = а) и связующего (п = с) , -(к), —к) - интенсивности армирования в плоскости и по толщине к-го слоя.
Установленные уравнения модели однонаправленно армированного слоя достаточно просты и их использование в конкретных расчетах не вызывает никаких затруднений. В то же время эти уравнения при
Ea ^ Ec , v
обеспечивают непрерывный переход к уравнениям Дюамеля-Неймана однородного изотропного тела.
На практике оболочка нередко изготовлена из материалов однонаправленно армированных слоев, оси симметрии которых не совпадают с координатными линиями исходной поверхности. Такая ситуация имеет место, например, в перекрестно армированных (под углами к направлению Г = const) оболочках.
Преобразование физических компонентов эффективных механических (16), температурных (17) жесткостей и тензора поперечных сдвиговых податливостей осуществляется по известным формулам [22]:
A(k) - ) m(k)am(k)p m(km(k(к) _ д(к) aR
(™p5) _ (ap^)m " m mm Pm 5 ' _ 9(«Р)m ^'
(k) r.(k) a p C(*¿) _ c(ap)mX ml •
(22)
Здесь m(k^ _ m(k)2 _ cos у к, m(k)2 _-m(k)2 _ sin у к.
При этом физические составляющие тензоров эффективных тангенциальных жесткостей и температурных напряжений армированного материала зависят от модулей Юнга материалов арматуры и связующего, коэффициентов Пуассона арматуры и связующего, интенсивностей армирования и направления укладки волокон. Используемая структурная модель армированного материала позволяет после решения краевой задачи и определения осредненных напряжений находить напряжения в элементах композиции - арматуре и связующем [21-23].
Кроме того, можно учитывать неоднородность материала оболочки, находя для различных структур армирования слоев зависимость структурных параметров
-(k), —k), уk от координат срединной поверхности [24].
Заключение
Таким образом, расчет армированной геликоидальной оболочки сводится к решению краевой задачи, описываемой дифференциальными уравнениями равновесия (11), кинематическими соотношениями (9), физическими соотношениями (15), а также краевыми условиями на прямолинейных и винтовых кромках оболочки (12).
В отличие от [24], где подобные соотношения получены для оболочки, срединная поверхность которой отнесена к главным линиям кривизны, приведенные здесь уравнения получены для оболочки, отнесенной к асимптотическим линиям.
Если оболочка шарнирно оперта вдоль прямолинейных кромок, то краевые условия (12) при ф = 0, Лф примут вид
"0ф= ™ = Ъг = Я = Мф= Нф= 0. (23)
Условиям (23) можно полностью удовлетворить, если обобщенные смещения представить в виде тригонометрических рядов:
да да
w(r, ф) = Е ™к(г^п ькф; и0ф(г, ф) = Е м0ф(г^п ькф; к=1 к=1
иог (r, ф) = Е мог(г)с°8 ькф; уэг (г, ф) = еуЗг (г^ькф;
к= к= (24)
то
У3ф(г, Ф) = Е у3ф(г Ькф. к=1
В отличие от классической краевой задачи, описываемой дифференциальными уравнениями восьмого порядка [24], в рассматриваемой постановке порядок дифференциальных уравнений двенадцатый. Подставляя (24) в (10) с помощью (13), (15) получим следующие выражения для усилий и моментов:
то то
тг(ф) = Е тг(ф)к(г)с^ькф; ^ = Е Бк(г^пькф; к=1 к=1
(25)
мг(ф) = Емг(Ф)к(г^Ькф; МГф = ЕмГф(гЬкф; к=1 к=1
то то
нг(ф) = Е нг(Ф)к(г^пЬкф; нГф = Е нкГф (г)сов ЬкФ; к=1 к=1
то то
Qr = Е бгк(г^Ькф; 6Ф = Е Qфk(г)со§ЬкФ. к=1 к=1
Входящие в (25) функции Тгк(г), Тфк(г),..., Оф.(г) выражаются через обобщенные кинематические перемещения и компоненты тензора упругости.
Подставляя далее (25) в (11) и применяя процедуру Бубнова-Галеркина, получим бесчисленное множество систем обыкновенных дифференциальных уравнений двенадцатого порядка:
(ЛТк) + ЬкЯк -^ТТфк = 0; ОТ (Л>%) + -- ЬкТфк = 0;
аг аг аг аг
а2 ( а м ) 2Ь амГф Ьк Л нк Ьк м а (ал2 м ^ + — (лмгк) - -лу нгФ- А2 мФк - аг Ц2мФк J+
+ ^ = -Л Як; (26)
— (л2Нгк ) - ЬкНГф —Нфк - Л20^гк = 0;
а((Нф) + аЛ2Нк + ЬкНфк -Л20фк = 0,
где Як - коэффициенты разложения в ряд Фурье нормальной и распределенной нагрузки я(г, ф).
Применяя процедуру Бубнова-Галеркина к (12), получим также набор краевых условий на винтовых кромках для уравнений (26) при г = Г(, Г2 .
Таким образом, анализ напряженно-деформированного состояния армированной геликоидальной оболочки сводится к решению краевой задачи, описываемой множеством систем обыкновенных дифференциальных уравнений двенадцатого
порядка. Если же рассматривать оболочку в рамках классической теории Кирхгофа-Лява, то порядок разрешающих систем будет восьмой [24]. Возрастание порядка разрешающих систем дифференциальных уравнений вызывает значительное повышение интенсивности краевых эффектов по сравнению с классической теорией оболочек. Поэтому хорошо зарекомендовавшие себя численные методы решения краевых задач классической теории дискретной ортогонализации здесь малопригодны. Наиболее перспективным для решения поставленной краевой задачи представляется метод инвариантных погружений [21], разработанный для решения многоточечных краевых задач, описываемых дифференциальными уравнениями высоких порядков.
Учет поперечного сдвига необходим при анализе напряженно-деформированного состояния оболочек, изготовленных из материала армированного высокомодульными волокнами Ea » Ec, например боро- или углепластика. В то же время расчет оболочек, изготовленных из материалов с сопоставимыми механическими характеристиками арматуры и связующего, может проводиться в рамках классической теорий оболочек Кирхгофа-Лява.
ЛИТЕРАТУРА
1. Колтунов С.Я. К расчету напряженного состояния в конечных геликоидальных оболочках // Известия Академии наук СССР. Механика твердого тела. - 1980. - № 6. -C. 149-152.
2. Ройтман В.М. Инженерные решения по оценке огнестойкости проектируемых и реконструируемых зданий. - М.: Ассоциация «Пожарная безопасность и наука», 2001. -382 с.
3. Муц А.В. Выбор переменных проектирования при оптимизации последовательности укладки конструкций из слоистых композитов // Механика композитных материалов. -2016. - Т. 52, № 2. - C. 305-324.
4. Muc A., Ulatowska A. Design of plates with curved fiber format // Composite Structures. -2010. - Vol. 92, N 7. - P. 1728-1733.
5. Muc A., Muc-Wersgon' M. An evolution strategy in structural optimization problems for plates and shells // Composite Structures. - 2012. - Vol. 94, N 4. - P. 1461-1470.
6. Джагангиров А.А. Несущая способность трехслойной волокнистой композитной кольцевой пластинки, защемленной по краям // Механика композитных материалов. -2016. - Т. 52, № 2. - С. 385-392.
7. Немировский Ю.В, Мищенко А.В. Динамический расчет системы профилированных композитных стержней // Вычислительная механика сплошных сред. - 2015. - Т. 8, № 2. - C. 188-199.
8. Емельянов И.Г., Кузнецов А.В. Применение виртуальных элементов при определении напряженного состояния оболочек вращения // Вычислительная механика сплошных сред. - 2014. - Т. 7, № 1. - C. 245-252.
9. Баженов В.Г., Павленкова Е.В., Артемьева А.А. Численное решение обобщенных осесимметричных задач динамики упруго-пластических оболочек вращения при больших деформациях // Вычислительная механика сплошных сред. - 2012. - Т. 5, № 4. - C. 427-434.
10. Киреев И.В., Немировский Ю.В. Консервативный численный метод решения линейных краевых задач статики упругих оболочек вращения // Вычислительная механика сплошных сред. - 2012. - Т. 5, № 1. - C. 85-99.
11. Залесский В.Н. Двумерная безмоментная задача для оболочки в виде косого геликоида // Динамика и прочность машин. - Харьков, 1974. - Вып. 12. - С. 88-93.
12. Залесский В.Н. Двумерная безмоментная задача для геликоидальной оболочки в асимптотических линиях // Динамика и прочность машин. - Харьков, 1966. - Вып. 3. -С. 10-13.
13. Гольденвейзер А.Л. Теория упругих тонких оболочек. - М.: Наука, 1976. - 512 с.
14. Александров П.В., Немировский Ю.В. Напряженное состояние гиперболической оболочки вращения при неоднородном армировании // Известия вузов. Строительство и архитектура. - 1991. - № 10. - С. 27-32.
15. Котельникова А.П. Расчет пологой геликоидальной оболочки, нагруженной изгибающим моментом // Известия вузов. Машиностроение. - 1978. - № 12. - С. 4-9.
16. O'Mathuna D. Rotationally symmetric deformations in helicoidal shells // Journal of Mathematics and Physics. - 1963. - Vol. 42. - P. 85-111.
17. Прочность при изгибе стеклопластиков, полученных намоткой / С.Г. Абрамов, Н.С. Мезенцев, В.П. Николаев, В.Д. Попов // Проблемы прочности. - 1975. - № 10. -С. 62-64.
18. Немировский Ю.В., Резников Б.С. Прочность элементов конструкций из композитных материалов. - Новосибирск: Наука, 1986. - 166 с.
19. Каган В.Ф. Основы теории поверхностей в тензорном изложении. Ч. 1. - М.; Л.: ОГИЗ, 1947. - 216 с.
20. Бабин А. И., Немировский Ю.В. Температурные напряжения в коллекторах машин постоянного тока // Численные методы решения задач теории упругости и пластичности: труды XVII межреспубликанской конференции, Новосибирск, 3-5 июля 2001 г. / под ред. В.М. Фомина. - Новосибирск, 2001. - С. 12-27.
21. Андреев А.Н. Математическая модель термоупругого деформирования слоистых композитных оболочек и пластин // Известия АлтГУ. - 2014. - № 1-1 (81). - С. 19-21. -doi: 10.14258/izvasu(2014)1.1-02.
22. Андреев А.Н., Немировский Ю.В. Численный анализ напряженно-деформированного состояния слоистых оболочек вращения методом инвариантного погружения // Известия Академии наук Армянской ССР. Механика. - 1989. - Т. 42, № 1. - С. 9-19.
23. Немировский Ю.В. Об упруго-пластическом поведении армированного слоя // Прикладная механика и техническая физика. - 1969. - № 6. - С. 81-89.
24. Александров П.В., Немировский Ю.В. Напряженное состояние армированных геликоидальных оболочек // Известия вузов. Строительство и архитектура. - 1991. - № 9. -С. 18-24.
25. Ярошенко А.Р. Осесимметричная деформация винтовой оболочки с прямоугольным профилем // Динамика и прочность машин. - Харьков, 1974. - Вып. 12. - С. 3-9.
THERMOSTRESSED STATE OF THE HELICOIDAL MULTILAYER POLYREINFORCED SHELLS
Nemirovskii Yu.V.1, Babin A.I. 2, Sal'skii E.A. 3
1 Khristianovich Institute of Theoretical and Applied Mechanics, Novosibirsk, Russia 2Kuzbass State Technical University, Kemerovo, Russia 3Kemerovo State University, Kemerovo, Russia
In general, a helicoidal shell construction can be composed of m reinforced layers different in the matrix composition, reinforcing fiber materials and the reinforcement pattern. Physical components of effective tangential stiffness and temperature stress tensors of a multilayered composite reinforced in different directions were obtained in the coordinate system not related to the material microstructure.
In this case, physical components of tangential effective stiffness and heat stress tensors depend on the Young modulus, the Poisson ratio of the matrix and reinforcing phase, the reinforcement power and the reinforcement fiber direction. The used structural model of the reinforced material allows finding stresses in composite compounds, namely the reinforcing phase and matrix, after computing mean values by solving the boundary problem.
An analytical solution of the contact thermal conductivity problem was found with a view to obtain a non-linear temperature distribution law over the reinforced shell with a negative Gaussian curvature. A closed equation system of multilayered helicoidal shell statics with regard to external thermal fields in the coordinate system normally related to the reference surface is adduced. The order of the system does not depend on the number of layers or method of
reinforcement. A calculation method of reinforced shells with a negative Gaussian curvature which takes in account the material structural heterogeneity is developed. It makes it possible by varying an internal structure to obtain maximum coherence between construction stiffness fields and tension fields from external thermal and mechanical effects.
Keywords: multi-reinforced shell; straight helicoids; temperature; non-coupled thermoelastic problem; Timoshenko's generalized hypotheses, Duhamel-Neumann's law.
DOI: 10.17212/1727-2769-2016-4-7-21
REFERENCES
1. Koltunov S.Ya. K raschetu napryazhennogo sostoyaniya v konechnykh gelikoidal'nykh ob-olochkakh [Finite helicoidal shells stress condition computation]. Izvestiya Akademii nauk SSSR. Mekhanika tverdogo tela - Mechanics of Solids, 1980, no. 6, pp. 149-152. (In Russian)
2. Roitman V.M. Inzhenernye resheniya po otsenke ognestoikosti proektiruemykh i rekon-struiruemykh zdanii [Engineer solutions in designed and restorated construction's fire resistance estimation]. Moscow, Assotsiatsiya "Pozharnaya bezopasnost' i nauka" Publ., 2001. 382 p.
3. Muc A. Vybor peremennykh proektirovaniya pri optimizatsii posledovatel'nosti ukladki kon-struktsii iz sloistykh kompozitov [Choice of design variables in the stacking sequence optimization for laminated structures]. Mekhanika kompozitnykh materialov - Mechanics of Composite Materials, 2016, vol. 52, no. 2, pp. 305-324. (In Russian)
4. Muc A., Ulatowska A. Design of plates with curved fiber format. Composite Structures, 2010, vol. 92, no. 7, pp. 1728-1733.
5. Muc A., Muc-Wersgon' M. An evolution strategy in structural optimization problems for plates and shells. Composite Structures, 2012, vol. 94, no. 4, pp. 1461-1470.
6. Dzhagangirov A.A. Nesushchaya sposobnost' trekhsloinoi voloknistoi kompozitnoi kol'tsevoi plastinki, zashchemlennoi po krayam [Load-carrying capacity of a fiber-reinforced annular tree-layer composite plate clamped on its external and internal contours]. Mekhanika kompozitnykh materialov - Mechanics of Composite Materials, 2016, vol. 52, no. 2, pp. 385398. (In Russian)
7. Nemirovskii Yu.V., Mishchenko A.V. Dinamicheskii raschet sistemy profilirovannykh kompozitnykh sterzhnei [Composite rods with variable cross-section]. Vychislitel'naya mekhanika sploshnykh sred - Computational Continuum Mechanics, 2015, vol. 8, no. 2, pp. 188-199.
8. Emel'yanov I.G., Kuznetsov A.V. Primenenie virtual'nykh elementov pri opredelenii naprya-zhennogo sostoyaniya obolochek vrashcheniya [Virlual elements usage in determinating shells of revolution stress condition]. Vychislitel'naya mekhanika sploshnykh sred - Computational Continuum Mechanics, 2014, vol. 7, no. 1, pp. 245-252.
9. Bazhenov V.G., Pavlenkova E.V., Artem'eva A.A. Chislennoe reshenie obobshchennykh osesimmetrichnykh zadach dinamiki uprugo-plasticheskikh obolochek vrashcheniya pri bol'shikh deformatsiyakh [Elastoplastic shells of revolution dynamic general axial symmetric problems with high deformations numerical solution]. Vychislitel'naya mekhanika splosh-nykh sred - Computational Continuum Mechanics, 2012, vol. 5, no. 4, pp. 427-434.
10. Kireev I.V., Nemirovskii Yu.V. Konservativnyi chislennyi metod resheniya lineinykh kraevykh zadach statiki uprugikh obolochek vrashcheniya [Elastic shells of revolution linear boundary static problems conservative numeric solutions methods]. Vychislitel'naya mekhanika sploshnykh sred - Computational Continuum Mechanics, 2012, vol. 5, no. 1, pp. 85-99.
11. Zalesskii V.N. [2-dimensional skew helicoid shell momentless problem]. Dinamika i prochnost' mashin [Engine dynamics and durability]. Khar'kov, 1974, iss. 12, pp. 88-93. (In Russian)
12. Zalesskii V.N. [2-dimensional asymptotical line helicoidal shell momentless problem]. Dinamika i prochnost' mashin [Engine dynamics and durability]. Khar'kov, 1966, iss. 3, pp. 10-13. (In Russian)
13. Gol'denveizer A.L. Teoriya uprugikh tonkikh obolochek [Elastic thin shell theory]. Moscow, Nauka Publ., 1976. 512 p.
14. Aleksandrov P.V., Nemirovskii Yu.V. Napryazhennoe sostoyanie giperbolicheskoi oboloch-ki vrashcheniya pri neodnorodnom armirovanii [Hyperboloid nonhomogenously armored shell of revolution strained state]. Izvestiya vysshikh uchebnykh zavedenii. Stroitel'stvo i arkhitektura - News of Higher Educational Institutions. Construction and Architecture, 1991, no. 10, pp. 27-32.
15. Kotel'nikova A.P. Raschet pologoi gelikoidal'noi obolochki, nagruzhennoi izgibayushchim momentom [Stressed with bending moment flat hyperboloid computation]. Izvestiya vysshikh uchebnykh zavedenii. Mashinostroenie - Proceedings of Higher Educational Institutions. Маchine Building, 1978, no. 12, pp. 4-9.
16. O'Mathuna D. Rotationally symmetric deformations in helicoidal shells. Journal of Mathematics and Physics, 1963, vol. 42, pp. 85-111.
17. Abramov S.G., Mezentsev N.S, Nikolaev V.P., Popov V.D. Prochnost' pri izgibe stekloplas-tikov, poluchennykh namotkoi [Bending strength of glass-fiber-reinforced plastics produced by winding]. Problemy prochnosti - Strength of Materials, 1975, no. 10, pp. 62-64. (In Russian)
18. Nemirovskii Yu.V., Reznikov B.S. Prochnost' elementov konstruktsii iz kompozitnykh mate-rialov [Composite construction element resistance]. Novosibirsk, Nauka Publ., 1986. 166 p.
19. Kagan V.F. Osnovy teorii poverkhnostei v tenzornom izlozhenii. Ch. 1 [Surface theory bases in tensor terms. Pt. 1]. Moscow, Leningrad, OGIZ Publ., 1947. 216 p.
20. Babin A.I., Nemirovskii Yu.V. [Thermal stresses in direct current engine's collectors]. Chislen-nye metody resheniya zadach teorii uprugosti i plastichnosti: trudy XVII mezhrespublikanskoi konferentsii [Numerical methods in elasticity and plasticity problems: proceedings of the XVII inter-republican conference], Novosibirsk, 3-5 July 2001, pp. 12-27. (In Russian)
21. Andreev A.N. Matematicheskaya model' termouprugogo deformirovaniya sloistykh kompozitnykh obolochek i plastin [A mathematical model of thermoelastic deformation of layered composite shells and plates]. Izvestiya Altaiskogo gosudarstvennogo universiteta -The News of Altai State University, 2014, no. 1-1 (81), pp. 19-21. doi: 10.14258/ izvasu(2014)1.1-02
22. Andreev A.N., Nemirovskii Yu.V. Chislennyi analiz napryazhenno-deformirovannogo sos-toyaniya sloistykh obolochek vrashcheniya metodom invariantnogo pogruzheniya [Numerical analysis of rotation sandwich shells stress-strain state by invariant imbedding method]. Izvestiya Akademii nauk Armyanskoi SSR. Mekhanika, 1989, vol. 42, no. 1, pp. 9-19.
23. Nemirovskii Yu.V. Ob uprugo-plasticheskom povedenii armirovannogo sloya [Armored layer elastic-yelding behavior]. Prikladnaya mekhanika i tekhnicheskaya fizika - Journal of Applied Mechanics and Technical Physics, 1969, no. 6, pp. 81-89. (In Russian)
24. Aleksandrov P.V., Nemirovskii Yu.V. Napryazhennoe sostoyanie armirovannykh ge-likoidal'nykh obolochek [Armored helicoidal shells of rotation stress condition]. Izvestiya vysshikh uchebnykh zavedenii. Stroitel'stvo i arkhitektura - News of Higher Educational Institutions. Construction and Architecture, 1991, no. 9, pp. 18-24.
25. Yaroshenko A.R. [Axisymmetric deformation of the screw shell with a rectangular profile]. Dinamika i prochnost' mashin [Engine dynamics and durability]. Khar'kov, 1974, iss. 12, pp. 3-9.
СВЕДЕНИЯ ОБ АВТОРАХ
Немировский Юрий Владимирович - родился в 1936 году, д-р физ.-мат. наук, профессор, главный научный сотрудник, СО РАН, Институт теоретической и прикладной механики имени С .А. Христи-ановича. Область научных интересов: теория пластин и оболочек, теория термоупругости и термопрочности, теория теплопроводности и тепломассопереноса в композитных конструкциях, теория пластичности ползучести, динамики тонкостенных конструкций, механика композитных материалов и конструкций, теория оптимального проектирования. Опубликовано 455 научных работ. (Адрес: 630099, Россия, Новосибирск, Морской пр-т, 44. E-mail: nemiryury@mail.ru).
Nemirovskii YuriiVladimirovich (b. 1936) - Doctor of Sciences (Phys.&Math.), professor, chief research scientist, SB RAS, Khristia-novich Institute of Theoretical and Applied Mechanics, Novosibirsk, Russia. His research interests are currently focused on the plate and shell theory, the theory of thermal elasticity and thermal resistance, the thermal conduction and heat-and-mass transfer theory in composite structures, the thin-wall structure dynamics theory, mechanics of composites, and the optimal design theory. He is the author of 455 scientific papers. (Address: 44, Morskoi Prospekt, Novosibirsk, 630099, Russia. E-mail: nemiryury@mail.ru).
Бабин Анатолий Иванович - родился в 1953 году, старший преподаватель, Кузбасский государственный технический университет имени Т.Ф. Горбачева. Область научных интересов: связанная задача термоупругости, теория оболочек. Опубликовано 23 научные работы. (Адрес: 650055, Россия, Кемерово, Ворошилова, 24в. E-mail: anbabin@yandex.ru).
Babin Anatoly Ivanovich (b. 1953) - senior lecturer, Kuzbass State Technical University. His research interests are currently focused on the coupled thermal elasticity theory and the shell theory. He is the author of 23 scientific papers. (Address: 24c, Voroshilova St, Kemerovo, 650055, Russia. E-mail: anbabin@yandex.ru).
Сальский Евгений Алексеевич - родился в 1989 году, аспирант, Кемеровский государственный университет. Область научных интересов: теория оболочек, механика твердого тела. Опубликовано 3 научные работы. (Адрес: 650002, Россия, Кемерово, ул. Институтская, 12. Email: e_s_@mail.ru).
Sal'skii Evgeny Alekseevich (b. 1989) - postgraduate, Kemerovo State University. His research interests are currently focused on the shell theory and solid-state mechanics. He is the author of 3 scientific papers. (Address: 12, Institutskaia St, Kemerovo, 650002, Russia. E-mail: e_s_@mail.ru).
Статья поступила 09 января 2017 г.
Received January 9, 2017
To Reference:
Nemirovskii Yu.V., Babin A.I., Sal'skii E.A. Termonapryazhennoe sostoyanie mnogosloinykh poliarmirovannykh gelikoidal'nykh obolochek [Thermo-stressed state of the helicoidal multilayer polyreinforced shells]. Doklady Akademii nauk vysshei shkoly Rossiiskoi Federatsii -Proceedings of the Russian higher school Academy of sciences, 2016, no. 4 (33), pp. 7-21. doi: 10.17212/1727-2769-2016-4-7-21