Переменный изгиб трехслойного стержня со сжимаемым заполнителем в нейтронном потоке Текст научной статьи по специальности «Физика»

УДК 539.374

ПЕРЕМЕННЫЙ ИЗГИБ ТРЕХСЛОЙНОГО СТЕРЖНЯ СО СЖИМАЕМЫМ ЗАПОЛНИТЕЛЕМ В НЕЙТРОННОМ ПОТОКЕ

Э. И. Старовойтов1, Д. В. Леоненко2

1 Старовойтов Эдуард Иванович, доктор физико-математических наук, заведующий кафедрой строительной механики Белорусского государственного университета транспорта, 246653, Беларусь, Гомель, Кирова, 34, [email protected]

2Леоненко Денис Владимирович, доктор физико-математических наук, профессор кафедры строительной механики Белорусского государственного университета транспорта, 246653, Беларусь, Гомель, Кирова, 34, [email protected]

Рассмотрен переменный изгиб трехслойного упругопластического стержня со сжимаемым заполнителем в нейтронном потоке. Для описания кинематики несимметричного по толщине пакета приняты гипотезы ломаной линии: в тонких несущих слоях справедливы гипотезы Бернулли; в сжимаемом по толщине заполнителе выполняется гипотеза Тимошенко с линейной аппроксимацией перемещений по толщине слоя. Учитывается работа заполнителя в тангенциальном направлении. Физические соотношения связи напряжений и деформаций соответствуют теории малых упругопластических деформаций. Система дифференциальных уравнений равновесия получена вариационным методом. На границе предполагаются кинематические условия свободного опирания торцов стержня на неподвижные в пространстве жесткие опоры. Решение краевой задачи сведено к нахождению четырех искомых функций — прогибов и продольных перемещений срединных поверхностей несущих слоев. Аналитическое решение получено методом упругих решений с помощью теоремы о переменных нагружениях Москвитина. Проведен его численный анализ в случае непрерывных и локально распределенных нагрузок.

Ключевые слова: трехслойный стержень, пластичность, повторный изгиб, сжимаемый заполнитель, нейтронный поток.

10.18500/1816-9791-2017-17-2-196-208

ВВЕДЕНИЕ

В связи с широким применением в промышленности и строительстве неоднородных элементов конструкций при квазистатических и динамических нагрузках в различных физико-механических средах актуальной становится проблема создания адекватных математических моделей для описания их деформирования. В монографиях [1,2] рассмотрены различные математические модели квазистатического и динамического деформирования слоистых элементов конструкций, приведены постановки краевых задач, изложены методы их решения. В статьях [3-10] исследованы свободные и вынужденные колебания неоднородных оболочек, а также рассмотрены вопросы распространения волн и деформирования оболочек и пластин с термочувствительной толщиной. Публикации [11-15] посвящены теоретическому и экспериментальному исследованию параметров свободных колебаний трехслойных панелей и квазистатическому деформированию трехслойных упругих и упругопластических пластин при однократных и переменных нагрузках.

Здесь теория переменного нагружения Москвитина [15] применена к трехслойным упругопластическим стержням прямоугольного поперечного сечения со сжимаемым заполнителем, находящимся в нейтронном потоке.

1. НАГРУЖЕНИЕ ИЗ ЕСТЕСТВЕННОГО СОСТОЯНИЯ

Рассматривается несимметричный по толщине трехслойный стержень (рис. 1). Для изотропных несущих слоёв приняты гипотезы Бернулли. На границах контакта используются условия непрерывности перемещений. Материалы несущих слоев несжимаемы в поперечном направлении, в заполнителе учитывается его обжатие, деформации малые. Система координат х, у, £ связывается со срединной плоскостью заполнителя. На стержень действуют поверхностная нагрузка д(х) и подводится нейтронный поток плотностью в направлении, противоположном внешней нормали. Через ш^(х) и щ (х) обозначены прогибы и продольные перемещения срединных поверхностей несущих слоёв, Нк — толщина к-го слоя, Н3 = 2 (к = 1,2, 3 — номер слоя), Ь0 — ширина стержня.

Рис. 1. Расчетная схема трехслойного стержня Fig. 1. The design scheme of a three-layer rod

Радиационное облучение твердых тел сопровождается многочисленными эффектами, в результате которых возникает объемная деформация От, изменяются упругие и особенно пластические характеристики материала. Согласно известным экспериментальным данным [16], рост величины нейтронного потока

I (£) = & (1)

в пределах малых деформаций, как правило, приводит к радиационному упрочнению материала, вызванному ростом предела текучести — интенсивность потока, нейтрон/(м2с), £ — время), и возникновению объемной деформации, которую в линейном приближении можно считать прямо пропорционально потоку (см. [16]):

От = В1 (£), (2)

где В — константа, получаемая из опыта.

Пусть на границу тела, занимающего полупространство z ^ 0, параллельно оси z падают нейтроны с одинаковой интенсивностью Тогда интенсивность потока нейтронов, доходящих до плоскости параллельной границы z = const, будет (см. [16]):

<p(z) = , д = ano = a——р, (3)

где д — величина макроскопического эффективного сечения (1/см), a — эффективное сечение, отнесенное к одному ядру, n0 — число ядер в 1 см3, —0 — число Авогадро, р — плотность, — — атомный вес.

Для алюминия a = 0.21 • 10-24 см2, —0 = 6.02 2-1023 моль-1, плотность р = 2.7 г/см3, — = 27 а.е.м., подставляя в (3), получим д = 1.26 м-1.

Если в (1) интенсивность ^ = не зависит от времени, то к моменту t через сечение с координатой z пройдет интегральный поток

I (z,t) = ^o te-MZ. (4)

В реакторах имеет порядок 1017 — 1018 нейтрон/(м2с), а I0 достигает значений 1023 — 1027 нейтрон/м2, причем 0i достигает значений порядка 0.1. Следовательно, величина B может быть порядка 10-28 — 10-23 м2/нейтрон.

Предположим, что в процессе деформирования материалы несущих слоев стержня проявляют упругопластические свойства, заполнитель — нелинейно упругий. Для описания их деформирования с учетом объемной деформации (2) принимаем физические уравнения состояния теории малых упругопластических деформаций [17]:

45 = 2Gk(1 — ^<k)(в <ч, I))e® a(k) = Kk (3£<k) — BkI) (k = 1, 2, 3),

S3 = 2Gk (1 — Ш<3'(в <3>, I ))eZ?, 43) = 2G3(1 — ш(3)(в<3), I ))e<3, '

где s(k), e(k), a(k), £(k) — девиаторные и шаровые части тензоров напряжений и деформаций, sXz, eXz — тангенциальное напряжение и сдвиговая деформация в заполнителе, Gk, Kk — модули сдвига и объемного деформирования материала k-го слоя, ш(к) — универсальные функции нелинейности материалов несущих слоев (k = 1,2) и заполнителя (k = 3) при нагружении из естественного состояния, например, для сплава Д16Т:

(k) mi' 0, sik) < £yk)(I),

(в(к) 1) а Л l<") (I) Г* (kK (k)m ,,

—k 1 — , >4 (I) (6)

^(k) (4k),i ) =

4k)(I ) = вуо 1 + —(1 — exp(—CI))

1/2'

(к) (к) где £и — интенсивность деформаций, е\ — деформационный предел текучести

материала несущего слоя, I — интегральный нейтронный поток (4), Вк — константы материалов, получаемые экспериментально.

Перемещения в слоях и(к)(ж,^} и можно выразить через четыре иско-

мые функции (ж), и (ж), эд2(ж) и и2(ж) — прогибы и продольные перемещения срединных поверхностей несущих слоев:

в несущих слоях

¿1

«(1) = «1 - (г - с - у ) '1 ,х, «(2) = «2 - (г + с + у ) '2,Х,

ш(1) = ш1 (с < г < с + ), ш(2) = Ш2 (-с - ¿2 < г < -с);

в заполнителе

«

(3) = (1 + ^ (1 „1 + £ + Г1 - -

Ш

(3) =

с) \ 2 1

с2

-«2

4

2

г

1

2

= - 1 + - '1 + - 1 - - '2 (-с < г < с)

(7)

где г — координата рассматриваемого волокна, запятая в нижнем индексе обозначает операцию дифференцирования по следующей за ней координате.

Компоненты тензора деформаций следуют из (7) и соотношений Коши [2]. Внутренние усилия в слоях вводятся соотношениями

(8)

-к) = ^ а м<к) = б0/ ^ а °(3) = ь/ ^ а

Ь Ьк Ьз

-<3) = ь/,*? ам<? = а,,

Ьз Ьз

(к) (3) (3)

где оХх , оХг, о\г — компоненты тензора напряжений.

Уравнения равновесия в усилиях для рассматриваемого стержня получены вариационным методом Лагранжа:

Н - Р1,х = 0, Н + Р2,х = 0,

,хх + Н2 - Т1 ,х = Ь0 ^2 ,хх - Н2 - Т2 ,х = 0

Обобщенные внутренние усилия в (9) выражаются через усилия (8)

(9)

Н = Т =

2с ' 1 +

2с

Н2 =

—

(3)

(3)

2с +

—Г м

Р = + '

2

(3)

х + -Х1),

0(3) мХ3)

2с

2с ¿1

Р2 =

—

(3)

2

¿1

мХ^

2с ¿1

+ -Х2),

* = (с+^ -х1)-мх1)+^-х3)+4^мх3), (10)

Т = 1 + ^

М 0(3) м

(3)

хг

2с 2

2с

2

& = - с + -х2) - мХ2) - ^ -Х3) + ^ мХ3).

4

4с

В дальнейшем принимаются кинематические граничные условия свободного опи-рания торцов стержня на неподвижные в пространстве жесткие опоры. Тогда в поперечных сечениях х = 0,1 (1 — длина стержня) должны выполняться требования

Шк = «к ,Х = Шк ,ХХ = 0

(к = 1, 2).

(11)

Н

1

2

2

Используя соотношения (5), выделим в тензоре напряжений упругие (индекс е) и нелинейные (индекс ,) слагаемые, которые будут включать и радиационные добавки (к = 1, 2):

о£к) = оХХ)е - оХк^, оХХ)е = 2£кеХк) + 3Кк£(к) = К+42,

о(к)- = 4£(к),(к) + Ккв 1 о(3) = о(3)е - о(3)-

°ХХ 3 икйХХш + Кквк1 , хх оХХ °ХХ ,

2

о(3)е = к+£(3) + К- 4?, = 2 С3 (2£<3» - £<3))ш(3) + К3В3

о<? = о<3»е - о<Г, о(3)е = К+4? + К-42, о«3;- = 2С3(24? - е(3>),<3> + КВ3

о(3) = о(3)е - о(3)- о(3)е = Д3) о(3)- = о^ -(3), ,(3)

42

К+ = К к + 3 Ск, К- = К к - 2 £к.

Проведя подобную операцию с внутренними усилиями (8), получим

—Хк) = —Хк)е - —Хк)-, мХк) = мХк)е - мХк)- (к = 1, 2, 3), —(3) = —3)е - —(3)-, м]3) = м]3)е - м]3)-, мХ3) = м^3)е - мХ3)-, (12)

д(3) = 0(3)е - 0(3)-.

Линейные (индекс е) и нелинейные слагаемых (индекс ,) в (12) будут вычисляться по формулам типа (8) , в которых подынтегральные выражения принимаются из (12) с соответствующими индексами.

Подставив в обобщенные усилия (10) внутренние усилия (12), в которых деформации выражены через перемещения (7), и отправив полученное в уравнения (9), придем к системе нелинейных дифференциальных уравнений, описывающей деформирование упругопластического стержня в нейтронном потоке. Для решения этой системы применим метод упругих решений Ильюшина [17]. В этом случае получим следующую систему линейных дифференциальных уравнений в итерациях:

Й1 - а!«П - а4«П,хх - «5«П,хх + ,х + ^3,х - 2ав'П,ххх +

+«7 'П,ххх = Р + Р-П-1),

-Й1«1 + Й1«2 - 0-5«1 ,хх - а9«2 )ХХ - «10Ш1 ,х - а17Ш2 ,х - «6Ш1 5ххх +

+2а7шП ,ххх = ¿-Т 1) 5 (13)

-°2«1 ,х + «10«2 'Х + 2аб«1 ,ххх + °6«2,ххх + «11 'Г5хх - «12Ш2 ?хх +

+«15Ш^ххх - «16Ш11 ,хххх + 08- Ов= д + ^-П_1), -03«Т,х + °17«2,х - «7«1 ,ххх - 2°7«Т,ххх - «12'! )хх + °14'2 )хх_ «16,хххх + а^'^хххх - ав+ авшn =

где п — номер линейного приближения, «Ц, — линейные приближения искомых перемещений, коэффициенты

^3 ^3 / ¿Л К3- ^3 / М + К3-

= , а2 = 1 + --^, Й3 = — 1 + — +

2с 2 2с 2 2 2с 2

^ 2К+ с а4 = К+ + 3

«в =

1 "1 К+

«5 =

К+с

2с

3 7 3 7

2К+ с

ад = К+ "2 + —3^-• с

ав =

К+ с"1

а7 =

К+с^

аз

ап =

К— "1

а12 =

К3- ("1 + "2) , Сзс

4

+

а1з =

к+"2 + К++ с"2

12

6

а14 =

1 +

2с К-"2

а15 =

К+"1 + К3+с"2

12

6

а1в =

К+

2 2"1

12

¿1

К

а10 = — 1+--+

2 V 2с/ 2

^ А2 _

2с/ 6 '

м _ ^с

2с/ 6 ;

ас 2

1 + М _ Озс 2с 6

аз Л + "2.

а17 = 1" 1 + _

К

В(п —1) 7 (п —1)

дополнительных нагрузках рш \ \

(п—1) (п—1) ^

Яш , дш на первом шаге прибли-

(1)

жения нелинейные слагаемые принимаются равными нулю (шк ) = 0), а в дальнейшем вычисляются по результатам предыдущей итерации:

р

(п—1) _ 1 ( ттш(п—1)_ рш(п— 1)

= Ьо V 1 р •

1

Я

(п—1) _ 1 ( пш(п—1)

•■• = Ог ^ '

гтш(п—1) грш(п—1)

+ Н2 _ Т 1 •

"(п—1) = 1 (

Пи Ьо V

Нш(п—1) + р»ш(п—1)

(14)

1

дш

(п—1) _ _ ( пш(п—1)_ ттш(п—1)_ грш(п—1)

, — , I ^о •тт Но То •

Ог

Обобщенные усилия в (14) вычисляются по формулам типа (10), в которых нужно добавить в верхний индекс знак нелинейности ш и номер приближения (п _ 1). Граничные условия (11) должны выполняться на каждом шаге приближения.

Применение метода упругих решений позволяет на каждом шаге приближения задачу об упругопластическом изгибе рассматриваемого трехслойного стержня в нейтронном потоке поле сводить к соответствующей линейной задаче упругости с дополнительными нагрузками.

2. АНАЛИТИЧЕСКОЕ РЕШЕНИЕ

Решение системы дифференциальных уравнений (13) принимается в виде разложения в тригонометрические ряды, которые автоматически удовлетворяют граничным условиям (11):

Е Щт соэ (

т=1

ж

дшж

1

Е^Гт 81п(

/пшж

1

ж

Еи

п

ь2 = / ^ и2т т=1 ж

соэ

птж

2

Е ^т (

/пшж

т=1 т=1

где и™т, Щт, ЖТт — искомые амплитуды перемещений.

1

6

6

2

2

2

2

ж

ип =

1

1

=

1

Поперечная нагрузка и дополнительные усилия в слоях стержня также представляются в виде разложений в тригонометрические ряды:

Я = Я- 8т

ш=1

птх

2 С / птх \

= 1 я(х) 8т ^ —-—^ ё х,

я:-1 =

Е

ш=1

птх

Я^1

У:—

я:-1

(х) 8Ш

птх

ё х,

яп-1 = V Я^1

•У: / ^ а:—

ш=1

81П

птх

Я^1

•У:—

я:-1

. ,птх\ (х) 81п ( —;— ) а х,

к:-1 = Е к:- оо^ (

ш=1

кп-1 = - I к

п-1

, птх\ (х) 008 ( —:— ) а х,

р:-1 =

Е

ш=1

Р:ш 008

птх

рП-1 :ш

р:- 1

(х) 008

птх

ё х.

Пусть поверхностная нагрузка равномерно распределена в интервале а < х < 6,

тогда

Я—

2Яо

пт

пт6 пта 008 —--- 008

1

1

(15)

После подстановки перемещений (15), дополнительных усилий (15) и коэффициентов (15) в уравнения равновесия (13) получим систему линейных алгебраических уравнений для определения искомых амплитуд перемещений Ц—, Ц—, ИИ—, И2—:

_,п-1 :ш '

Ь1Ц— + ВД- + Ьз и— + 64^2— = р 62иГ- + 65и2ш + 66ИГ- - 67и2— = :

63 Ц— + ЬбЦ— + 68 и— + 6дИ2— = Яш + Я

64 Ц— - б7и2— + б9 и— + б1оИ2— = я

п-1 :ш :

(16)

7п-1 7:ш '

где коэффициенты 6^ зависят от параметра т и от температуры через коэффициенты ап, определенные в (13).

Решение системы (16) можно выписать в определителях, либо решать численно. Далее по формулам (15) вычисляются искомые функции. Перемещения в несущих слоях и заполнителе следуют из соотношений (7), деформации — из соотношений Коши, напряжения — из (5).

3. ПОВТОРНОЕ ЗНАКОПЕРЕМЕННОЕ НАГРУЖЕНИЕ

Пусть, начиная с момента £ = , осуществляется мгновенная разгрузка и повторное нагружение усилиями обратного знака я''. Они создадут в к-м слое стержня поле перемещений и^'(х), эд2(х), деформаций (х,г) и напряжений (х,г), г,^ = х,у. При этом будем предполагать, что за время разгрузки и последующего переменного нагружения воздействие нейтронного потока прекратилось, интегральный нейтронный поток зафиксирован и совпадает с уровнем к моменту начала разгрузки.

х

Я

—

о

2

1

1

1

о

х

2

1

о

X

о

2

1

1

1

о

Пределы текучести материалов , а"ук) в каждой точке стержня зависят, вообще говоря, от величин , а'Ц^ и Д(¿,£1), существовавших непосредственно перед началом разгрузки — конкретно от того, происходило ли в данной точке тела пластическое деформирование или нет. Следовательно, материалы слоев после нагружения из естественного состояния становятся деформационно анизотропными.

Введем параметры со звездочками, в которых величины с одним штрихом — напряжения, деформации и перемещения в стержне перед разгрузкой, двумя штрихами помечены аналогичные параметры в процессе второго полуцикла:

(к)* /(к) //(&) (к)* /Ш //Ш * / // /7 1 о оЛ

аЬ = аЬ - а Ь , 4 = - е Ь , Я = Я — Я (к = 1, 2, з),

*/// */// */// */ / ^ '

и1 = и1 — и15 = — и2 = и2 — и2, = — .

Для введенных напряжений и деформаций со звездочками примем физические уравнения состояния типа (5):

45* = 2^(1 — ^(к)*(4к)*,/1 ))4х*> а(к)* = 3Кк е(к)* (к = 1, 2, 3), 4?* = (1 — ^(3)*(г!3)*,/1 ))е!3;*, 43)* = 2^з(1 — ^(3)*(43)*,/1 М?*. (18)

Отличие соотношений (18) в том, что в связь шаровых частей тензоров интегральный нейтронный поток не входит.

Следует отметить, что выражения (7)-(11) и соотношения Коши в силу их линейности будут выполняться и для величин со звездочками. После проведения аналогичных преобразований получим систему нелинейных дифференциальных уравнений равновесия в итерациях для перемещений со звездочками, которая будет подобна си-

/10\ ТТ *(п — 1) 7 *(п — 1) *(п —1) *(п —1)

стеме (13). Дополнительные «внешние» нагрузки \ л/ \ я/ \ ; в ней вычисляются по формулам (14), в которых нелинейные составляющие будут иметь вид

жХк)/*(п-1) = -Ьо / ок4к)*(п-1)^(к)*(4к)*(п-1), /1

-

-I

3

мХк)/*(п-1) = -Ьо / ^£Хк)*(п-1)^(к)*(4к)*(п-1), /1)^,

3 Ьо / ^3(2^

жХ3)/*(п-1) = ±Ьо / ^3(2гХ3)*(п-1) — 43)*(п-1))^(3)*(43)*(п-1), /1

Ьз

Жр*^"= 3Ьо / ^31} — 43)*(п"*(43)*(п"1}, /1

Ьз

М^^-1} = 2 Ьо/ ) — 43М"-*(43)*(п"1}, /1

Ьз

м^з)/*(п-1) = 2Ьо J ^3£Х3)*(п-1)^(3)*(п-1)(£и3)*(п-1), д)^,

Ьз

д(3)/*(п-1) = 2Ьо J ^3е(3)*(п-1)^(3)*(п-1)(43)*(п-1), Д)^ (к = 1, 2).

Ьз

Входящая в них функция нелинейности со звездочкой ш(к)к(п' , 11), со-

гласно гипотезе Москвитина, может быть выражена через соответствующую функцию нелинейности при нагружении из естественного состояния (6):

ш(к)* (4^,11) = и(к) (£{ик>,11),

но с другими константами Л*к, а*к.

В случае граничных условий свободного опирания стержня по торцам на неподвижные в пространстве жесткие опоры краевая задача для величин со звездочками с точностью до обозначений совпадает с задачей о нагружении стержня из естественного состояния. Поэтому искомое аналитическое решение краевой задачи для величин со звездочками можно принять в тригонометрических рядах по типу (15). Входящие в них искомые амплитуды перемещений Ц™, , , определяются по рассмотренной ранее схеме. Искомые перемещения на втором полуцикле следуют из (17).

4. ЧИСЛЕННЫЕ РЕЗУЛЬТАТЫ

Численно исследованы перемещения стержня, слои которого составлены из материалов Д16Т-фторопласт-Д16Т. Относительные толщины слоев Н1 = к2 = 0.03, с = 0.09; упругие и радиационные параметры материалов приведены в [2]. Кроме этого принималось: В = 10-23 м2/нейтрон, что обеспечивает объемную деформацию в слоях рассматриваемого стержня 3-3.5 %; для дюралюминия — д = 1.26 см, для фторопласта д = 3.21 см. Воздействие нейтронного потока моделировалось увеличением предела текучести материала внешнего слоя на 20 %, предела физической нелинейности заполнителя — на 10 %.

На рис. 2 показано изменение прогибов первого слоя по середине стержня (х = 0,5) в зависимости от координат краев интервала нагрузки: 1 — упругий, 2', 2" — упругопластический при прямом и обратном нагружениях в нейтронном потоке.

0.04 W1 0.02

0.

-0.02 -0.04

J'1___

............

0

0.2

0.4 0.6

а / a

0.8

b 1.0

Wi 0.010

0.005

0

-0.005 -0.010 -0.015

2"

ss.

1

"1'....... ..........

0

0.15

0.30 0.45

б/b

0.60 а 0.75

Рис. 2. Изменение прогиба вдоль оси стержня в зависимости от длины интервала нагрузки Fig. 2. Change in deflection along the axis of the rod, depending on the length of the load

interval

Рис. 2, а иллюстрирует зависимость прогиба от длины интервала нагрузки b (a = 0); рис. 2, б — зависимость прогиба от координаты a левого края интервала нагрузки (b — a = 0.25).

Следует отметить, что за счет воздействия нейтронного потока на первом полуцикле нагружения происходит увеличение пределов текучести и жесткости материалов слоев. Это приводит к уменьшению упругопластического прогиба на 3-4 %. При повторном знакопеременном нагружении и равномерно распределенной и локальной прямоугольной нагрузками предел текучести материала несущих слоев изменяется только за счет силового упрочнения материала, так как воздействие нейтронного потока прекращается. Поэтому максимальное значение прогибов при повторных нагру-жениях в нейтронном потоке уменьшается на 1-2 %, однако в случае последующих циклических изменениях нагрузки это расхождение будет накапливаться.

ЗАКЛЮЧЕНИЕ

Предложенная методика исследования переменного деформирования трехслойных стержней со сжимаемым заполнителем может быть применена при любом количестве циклов нагружения. Она позволяет учитывать в инженерных расчетах воздействие нейтронного облучения и физическую нелинейность материалов слоев. Численные расчеты показали их существенное влияние на перемещения в стержне.

Благодарности. Работа выполнена при финансовой поддержке Российского научного фонда (проект № 14-49-00091-П).

Библиографический список

1. Болотин В. В., Новичков Ю. Н. Механика многослойных конструкций. М. : Машиностроение, 1980. 375 с.

2. Горшков А. Г., Старовойтов Э. И., Тарлаковский Д. В. Теория упругости и пластичности. М. : Физматлит, 2011. 416 с.

3. Kuznetsova E. L., Leonenko D. V., Starovoitov E. I. Natural Vibrations of Three-layer Circular Cylindrical Shells in an Elastic Medium // Mech. Solids. 2015. Vol. 50, № 3. P. 359-366. DOI: 10.3103/S0025654415030127.

4. Леоненко Д. В., Старовойтов Э. И. Импульсные воздействия на трехслойные круговые цилиндрические оболочки в упругой среде // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 1. С. 202-209. DOI: 10.18500/18169791-2015-15-2-202-209.

5. Мочалин А. А. Параметрические колебания неоднородной круговой цилиндрической оболочки переменной плотности при различных краевых условиях // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 2. С. 210-214. DOI: 10.18500/1816-9791-2015-15-2-210-215.

6. Тарлаковский Д. В., Федотенков Г. В. Пространственное нестационарное движение упругой сферической оболочки // Изв. РАН. МТТ. 2015. Т. 50, № 2. С. 118-128.

7. Tarlakovskii D. V.Fedotenkov G. V. Two-Dimensional Nonstationary Contact of Elastic Cylindrical or Spherical Shells //J. Machinery Manufacture and Reliability. 2014. Vol. 43, № 2. P. 145-152. DOI: 10.3103/S1052618814010178.

8. Kuznetsova E. L., Tarlakovskii D. V.Fedotenkov G. V. Propagation of Unsteady Waves in an Elastic Layer // Mech. Solids. 2011. Vol. 46, № 5. P. 779-787. DOI: 10.3103/S0025654411050128.

9. Fedotenkov G. V.Tarlakovskiy D. V. Analytic Investigation of Features of Stresses in Plane Nonstationary Contact Problems With Moving Boundaries //J. Math. Sci. 2009. Vol. 162, № 2. P. 246-253. DOI: 10.1007/s10958-009-9635-4.

10. Белосточный Г. Н., Русина Е. А. Оболочки и геометрически нерегулярные пластинки с термочувствительной толщиной // Докл. Рос. акад. естественных наук. Поволжское межрегион. отд-ние. 1999. № 1. С. 28-37.

11. Rabboh S., Bondok N., Mahmoud T., El Kholy H. The Effect of Functionally Graded Materials into the Sandwich Beam Dynamic Performance // Materials Sciences and Applications. 2013. Vol. 4. P. 751-760. DOI: 10.4236/msa.2013.411095.

12. Havaldar S., Sharma R. Experimental Investigation of Dynamic Characteristics of Multilayer PU Foam Sandwich Panels //J. Minerals and Materials Characterization and Engineering. 2013. Vol. 1, № 5. P. 201-206. DOI: 10.4236/jmmce.2013.15031.

13. Yang L., Harrysson O., West H., Cormier D. A. Comparison of Bending Properties for Cellular Core Sandwich Panels // Materials Sciences and Applications. 2013. Vol. 4, № 8. P. 471-477. DOI: 10.4236/msa.2013.48057.

14. Москвитин В. В. Циклическое нагружение элементов конструкций. М. : Наука, 1981. 344 с.

15. Москвитин В. В., Старовойтов Э. И. К исследованию напряженно-деформированного состояния двухслойных металлополимерных пластин при циклических нагружениях // Изв. АН СССР. МТТ. 1986. № 1. С. 116-121.

16. Ильюшин А. А., Огибалов П. М. Упругопластические деформации полых цилиндров. М. : Изд-во Моск. ун-та, 1960. 224 с.

17. Ильюшин А. А. Пластичность. Ч. 1. Упругопластические деформации. М. : Гостехиздат, 1948. 376 с.

Образец для цитирования:

Старовойтов Э. И., Леоненко Д. В. Переменный изгиб трехслойного стержня со сжимаемым заполнителем в нейтронном потоке // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2017. Т. 17, вып. 2. С. 196-208. 001: 10.18500/1816-9791-2017-172-196-208.

Variable Bending of a Three-layer Rod with a Compressed Filler in the Neutron Flux

E. I. Starovoitov1, D. V. Leonenko2

1 Eduard I. Starovoitov, ORCID: 0000-0002-2550-5377, Belarusian State University of Transport, 34, Kirova str., Gomel, Belarus, 246653, [email protected]

2Denis V. Leonenko, ORCID: 0000-0001-8003-9279, Belarusian State University of Transport, 34, Kirova str., Gomel, Belarus, 246653, [email protected]

The present paper considers variable bending of a three-layer elastoplastic bar with a compressible filler in the neutron flux. To describe kinematic properties of an asymmetric through thickness pack we have accepted the hypotheses of a broken line as follows: Bernoulli's hypothesis is true in the thin bearing layers; Timoshenko's hypothesis is true in the compressible through thickness filler with a linear approximation of displacements through the layer thickness. The filler's work is taken into account in the tangential direction. The physical stress-strain relations correspond to the theory of small elastoplastic deformations. By the variational method a system of differential equilibrium equations has been derived. The kinematic conditions of simply supported faces of the bar on the immovable in space rigid bases are presumed on the boundary. The solution of the boundary problem is reduced to the search for four functions, namely: deflections and longitudinal displacements of the medial surfaces of the bearing layers. An analytical solution has been derived by the method of elastic solutions by using the Moskvitin's theorem about variable loadings. Its numerical analysis has been performed for the case of the uniform distribution of the continuous and local loads.

Key words: three-layer rod, plastic, variable bending, compressible filler, neutron flux.

Acknowledgements: This work was supported by the Russian Science Foundation (project no. 14-49-00091

References

1. Bolotin V. V., Novichkov Yu. N. Mekhanika mnogosloinykh konstruktsii [Mechanics of Multilayer Structures]. Moscow, Mashinostroenie, 1980. 375 p. (in Russian).

2. Gorshkov A. G., Starovoitov E. I., Tarlakovskii D. V. Teoriia uprugosti i plastichnosti [Foundations of the Theory of Elasticity and Plasticity]. Moscow, Fizmatlit, 2011. 416 p. (in Russian).

3. Kuznetsova E. L., Leonenko D. V., Starovoitov E. I. Natural vibrations of three-layer circular cylindrical shells in an elastic medium. Mech. Solids, 2015, vol. 50, no. 3, pp. 359— 366. DOI: 10.3103/S0025654415030127.

4. Leonenko D. V., Starovoitov E. I. Impulsive action on the three-layered circular cylindrical shells in elastic media. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2015, vol. 15, iss. 1, pp. 202-209. DOI: 10.18500/1816-9791-2015-15-2-202-209 (in Russian).

5. Mochalin A. A. Parametric Oscillations of a Non-uniform Circular Cylindrical Shell of Variable Density under Different Boundary Conditions. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2015, vol. 15, iss. 2, pp. 210-214. DOI: 10.18500/1816-9791-201515-2-210-215 (in Russian).

6. Tarlakovskii D. V., Fedotenkov G. V. Nonstationary 3D Motion of an Elastic Spherical Shell. Mech. Solids, 2015, vol. 50, no. 2, pp. 118-128. DOI: 10.3103/S0025654415020107.

7. Tarlakovskii D. V., Fedotenkov G. V. Two-Dimensional Nonstationary Contact of Elastic Cylindrical or Spherical Shells. Journal of Machinery Manufacture and Reliability, 2014, vol. 43, no. 2, pp. 145-152. DOI: 10.3103/S1052618814010178.

8. Kuznetsova E. L., Tarlakovskii D. V., Fedotenkov G. V. Propagation of Unsteady Waves in an Elastic Layer. Mech. Solids, 2011, vol. 46, no. 5, pp. 779-787. DOI: 10.3103/S0025654411050128.

9. Fedotenkov G. V., Tarlakovskiy D. V. Analytic Investigation of Features of Stresses in Plane Nonstationary Contact Problems With Moving Boundaries. J. Math. Sci., 2009, vol. 162, no. 2, pp. 246-253. DOI: 10.1007/s10958-009-9635-4.

10. Belostochny G. N., Rusina E. A. Obolochki and Geometrically Irregular Plates with a Temperature-Sensitive Thickness. Dokl. Ross. Akad. estestvennykh nauk. Povolzhskoe mezhregional'noe otdelenie , 1999, no. 1, pp. 28-37 (in Russian).

11. Rabboh S., Bondok N., Mahmoud T., El Kholy H. The Effect of Functionally Graded Materials into the Sandwich Beam Dynamic Performance. Materials Sciences and Applications, 2013, vol. 4, pp. 751-760. DOI: 10.4236/msa.2013.411095.

12. Havaldar S., Sharma R. Experimental Investigation of Dynamic Characteristics of Multilayer PU Foam Sandwich Panels. Journal of Minerals and Materials Characterization and Engineering, 2013, vol. 1, no. 5, pp. 201-206. DOI: 10.4236/jmmce.2013.15031.

13. Yang L., Harrysson O., West H., Cormier D. A. Comparison of Bending Properties for Cellular Core Sandwich Panels. Materials Sciences and Applications, 2013, vol. 4, no. 8, pp. 471-477. DOI: 10.4236/msa.2013.48057.

14. Moskvitin V. V. Tsiklicheskoe nagruzhenie elementov konstruktsii [Cyclic Loading of Elements of Designs]. Moscow, Nauka, 1981. 344 p. (in Russian).

15. Moskvitin V. V., Starovoitov E. I. To the Study of the Stress-Strain State of Two-Layer Metal-Polymer Plates under Cyclic Loading. Izv. AN SSSR. Mekhanika tverdogo tela, 1986, no. 1, pp. 116-121 (in Russian).

16. Il'iushin A. A., Ogibalov P. M. Uprugoplasticheskie deformatsii polykh tsilindrov [Elastic-Plastic Deformations of Hollow Cylinders]. Moscow, Moscow Univ. Press, 1960. 224 p. (in Russian).

17. Il'iushin A. A. Plastichnost'. Ch. 1. Uprugoplasticheskie deformatsii [Plastic. Ch. 1. Elastic-Plastic Deformation]. Moscow, Gostekhizdat, 1948. 376 p. (in Russian).

Cite this article as:

Starovoitov E. I., Leonenko D. V. Variable Bending of a Three-layer Rod with a Compressed Filler in the Neutron Flux. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2017, vol. 17, iss. 2, pp. 196-208 (in Russian). DOI: 10.18500/1816-9791-2017-17-2-196-208.