www.volsu.ru
DOI: https://doi.org/10.15688/mpcmjvolsu.2017.5.5
УДК 539.3 ББК 22.251
НАПРЯЖЕННО ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ТОЛСТОСТЕННОЙ ТРУБЫ ИЗ УПРУГОГО КОМПОЗИЦИОННОГО В СРЕДНЕМ ИЗОТРОПНОГО МАТЕРИАЛА
Александр Степанович Кравчук
Доктор физико-математических наук, доцент, профессор кафедры био- и наномеханики,
Белорусский государственный университет
просп. Независимости, 4, 220030 г. Минск, Белоруссия
Анжелика Ивановна Кравчук
Кандидат физико-математических наук,
доцент кафедры веб-технологий и компьютерного моделирования, Белорусский государственный университет [email protected]
просп. Независимости, 4, 220030 г. Минск, Белоруссия
Аннотация. Впервые решена краевая задача для твердого композиционного тела без использования нелокальных гипотез о малости объема композиционного материала, для которого устанавливаются эффективные характеристики. Предполагается, что разброс коэффициентов Пуассона около среднего значения мал, и поэтому можно использовать условие постоянства коэффициентов Пуассона и равенство его среднему значению для всех компонент композиционного материала. В связи с особенностями постановки задачи в случае композитного материала нет возможности отдельно рассматривать плоское напряженное состояние и плоскую деформацию поперечного сечения трубы, и оба этих состояния будут участвовать в оценке напряженно-деформированного состояния изучаемого объекта согласно стандартным гипотезам Фойгта и Рейсса. Получено приближение Хилла для средних по представительному объему напряжений и деформаций. В силу того, что напряжения и перемещения как на внутрен-^ ней, так и на внешней границе постоянны, то и их средние значения по участку границы сЗ любой площади постоянны и равны исходным значениям. В связи с этим с методической точки зрения для задачи Ляме для трубы показано, что решения, построенные с <1 использованием гипотез Фойгта и Рейса, самодостаточны и не требуют использования
И
дополнительных предположений о малости элементов усреднения внутри упругого тела.
м
^ Установлено, что решение по напряжениям зависит от средних значений модулей упругости по Фойгту и Рейссу, а по деформациям определяется только средними значения-^ ми модуля упругости по Рейссу. Получены формулы, определяющие напряженно-де-и формированное состояние композиционной в среднем изотропной плоскости с отвер-^ стием. Эти решения могут быть применены в качестве оценочных значений напряже-^ ний и деформаций в поперечном сечении грунта вокруг свай при, например, бурозабив-© ном способе их погружения в вечной мерзлоте.
Ключевые слова: композиционный структурно неоднородный материал, дискретная случайная величина, усреднение, эффективные деформационные характеристики, гипотеза Фойгта, гипотеза Рейсса, приближение Хилла.
Введение
Решение задачи Ляме для толстостенного цилиндра нашло широкое применение в инженерии и строительстве [3]. Кроме того, в настоящее время данная задача находит широкое применение в расчетах и оптимизации конструкций радиационно-тепловых экранов ядерных реакторов, их тепловых и биологических защит и т. д. [1]. Определение напряженного состояния трубы под внутренним давлением является неотъемлемой частью методик анализа дефектности трубопроводной системы [2].
Однако, несмотря на значительный интерес к этой задаче, ее решение для композиционного материала даже в случае упругости его компонент отсутствует. В частности, она может найти применение при расчете НДС труб из керамических материалов.
До настоящей статьи все решения краевых задач для композиционных тел выполнялись с явным использованием гипотезы о малости представительного объема композиционного тела по сравнению с его размерами. В данной статье впервые эта гипотеза отсутствует, вместо этого авторами явно применены гипотезы Фойгта и Рейсса к решению двух отдельных оценочных краевых задач и последующим усреднением полученных величин напряжений и деформаций.
Основные уравнения и гипотезы
Предполагается, что композиционный материал состоит из п компонент, к-я компонента которого (k = 1, п) имеет различные модули упругости Ек и одинаковые коэффициенты Пуассона
1 п
V = —У Vk (где Уь - коэффициенты Пуассона компонент). Относительная объемная доля каж-п ^^
ПЬ=1 п
дого материала ук (0 < ук < 1, ^ ук = 1). Будем обозначать через R1 и R2 радиусы внутренней и
к=1
внешней границ сечения цилиндра соответственно (рис. 1). При решении задачи будем использовать результаты для толстостенного цилиндра, изложенные в [3, с. 90-93].
Рис. 1. Толстостенный цилиндр
Гипотеза об однородности деформаций. Плоская деформация
В данном разделе проведем первый этап решения поставленной задачи, заключающийся в использовании гипотезы Фойгта о том, что при простейшем нагружении в материале имеет место однородная деформация, то есть для любого к-го слоя (к = 1, п) [4, с. 4]:
8гг,к = 8гг, 8ее,к = 8ее, 8ге,к = 8ге . (1)
Кроме того, из независимости решения от угла е следует, что 8ге = 0 .
В полярной системе координат в случае плоской осесимметричной деформации уравнение совместности деформаций имеет вид [3, с. 90]:
д8ее 8гг — 8ее = 0 (2)
дг г
а уравнения состояния определяются обобщенным законом Гука [3, с. 32]
-((1 ,к -^ее,к ),
s.... = ll 1 -vlc .-va„
Ek
see = ^ЧО-vKe,k-v°rr ,k 1 ^
Ek
где а ггк, сеек - нормальные компоненты напряжений к-го слоя (к = 1, п).
Непосредственно из (3) можно сразу же получить связь между средними значениями напряжений и деформаций по реализации модулей упругости композита. Поскольку известно, что относительная объемная доля каждого материала с модулем упругости Ек равна ук (к = 1, п), то из (3) следует очевидное равенство:
(Е)г-8гг =(1 + у|1 -уНаХ-у.(аее)г 1
{Е)у .8ее=(1 + у)((1 -у)-(аее)г-у.<а„)Д (4)
где
п п п
{Е)) =Хук • Ек , (агг)) =Еук -аггк , (аее)) =Еук 'аее,к .
к=1 к=1 к=1
Учитывая, что при гипотезе Фойгта 8гг = (8г^)) и 8ее = (8ее) , то для решения поставленной оценочной задачи остается расшифровать, что представляют собой (аг^) и (аее) .
1 дфк д 2фк
Далее, подставляя (3) в (2), используя выражения агг к =--, аее к =-получим уравне-
' г дг ' дг _
ние для определения функции напряжений Эри фк к-го концентрического слоя (к = 1, п) [3, с. 90-91 ]:
д 3Фк +1д 2Фк -1 =0 дг3 г дг2 г2 дг
решением которого является функция фк = С к 1п(г) + С2 к г 2 + С Таким образом, напряжения в к-м слое (к = 1, п) примут вид:
3,k-
Ci k Ci k
a rr, k = + 2 C 2, k, aee, k = + 2 C 2, k. (5)
r r
Из (1) и (3) следует, что для любых k и т (k, т = 1,п) выполнено:
1 + у
Ек
1 + у
((1 -V)
гг,k
Ьоо ="
Ей
((1 - ^сее„
- v)Сддк -vo
* )=^((1 -V)«
Ет
, к )=^«1 -V)«
е,т),
г,т )'
(6)
Решая уравнения (6) относительно сгг к и «ее к, будем иметь:
=Ек. = Ек.
«гг,к = Е"®гг,т , «ее,к = Е"СТев,т , т т
(7)
и с учетом (5) получим соотношения для коэффициентов функции ф к:
С = Ек с с = Ек с
С1,к _ Е с1,т , с 2,к _ Е с 2,т .
Ет Ет
(8)
Таким образом, нормальные напряжения на внешней и внутренней границах примут вид:
С1,1 , о Е
Е_
1
«ГГ,1 = ^Г + ^ С2,п , «гг,п ^Си— + 2 С2,и .
г2 Е
Е\ г
(9)
При граничных условиях сгг1| = ся и сггп| д = с^ из (9) получим:
Откуда следует
Я 2 • Е -с = Е • С + 2• Я 2 • Е •С
Я22 - Е1 - СЯ2 = Еп - С1,1 + 2 -Я22 - Е1 - С2,п .
С1,1 =
Е,
с Я--к-с
Я Е
Ч2 • Я22
^ Я22 - Я12
Я12 • ЕЕ" -С Я - Я22 -С я С, =-^
2,п
2
2•(Я12 - Я22)
(10)
Из (5), (8), (9) и (10) получаем
с„,, =
Ек Ек с Я--- •с Я
Я1 е Я2
Я
V Е1
• Я22 1 + 2 ' Ех
2 Ек 2 Ек
-■ с„ - Я, —- •
• —— • (
-'Я1 2
СТЕ
Я22 - Я12 г2
2 •(Я12 - Я22)
2 _ г. 2 Ек
С аа ь —
Ек Ек с Я--- •с Я
Я1 Е Я2
V Е1
\ 2 2 Я. —— •сЯ -Я2
^ Я12 • Я22 1 + 1 Е1 Я1 2
.-. г
Я22 - Я12 г2
2•Я -Я22)
(11)
Очевидно, поскольку и первым, и последним слоем могут быть любые материалы со своими модулями упругости Е{ и Е. (i = 1, п , ] = 1, п), то (11) необходимо переписать в виде:
Ь гг
Я
Я
п
r,k,i, j
Et E, g R--• G R
R e R2
E
j
R12 • R22 1
R22 - R12 r2
R 2 • ^ • . - R Ek
. -_ . /
+2
E
E
2 •( Ri2 - R22)
e,t,i, j
EkEk
k GR---
R1 E R2
Л R2 • R22 1
E
j
R22 - R,2 r2
(12)
+2
R,2--— • G R — R92--— • G R
1 e Ri E R2
2•(R12 -r22)
+
G
R
+
G
Далее, домножая на (12) у^ (^ = 1,п) и суммируя по k, получаем:
<■<■,',j ¡у
= ( E)v
^ - ч 1 R12 • R22 1
E. E.
R
R22 - R12 r2
■2 X El
2 • ZRl E.
-R
2 ^
! ^ E
2•(R12 -R22)
(Gee,i,j)y =4 E)y
^Gr1 GR2 1 R12 • R22 1
E. E
v ' j у
R22 - R12 r2
■2 < E\
R 2 • G R1 - R 2 • G r2
1 ' E. 2 ' E
2 •(R12 - R22) •
Домножая результат на у 1, суммируя по индексу I, далее домножая на у у и суммируя по индексу у, окончательно определяем и :
{E)v
(E).
(СТ R1 -СТ R2 )
± + 2
2 7 R2 - R,2 r2
R12 ^R1 - R22 ^R2
2•(RR12 -R22)
Er
>у = № -(-R1 -R2 )
RÜ ± + 2
2/R2 -R,2 r2
R12 • 0R1 - R22 ^R2
2 (R112 - R22)
(13)
где
Er =
in ^
ZZX Ex
vX=1 xy
-1
Из (4) с учетом (13) для гипотезы Фойгта получаем:
(в > = ^
\ rr/y
(вее) у =
(Er (1 + v)
(Er
(Cr Cr2 ) R22 - R12 r2
- (C^-CR2 )
R12 • R22 1 + 2•(!-2•v) R12 ^OR -R22^ '
2^2 - R22)
R12 • R22 1 ~ (l ~ ) R12 • ^R R2 •or
r^R27+2 (1 -2-v> 2R -R22)
(14)
Гипотеза об однородности напряжений. Плоское напряженное состояние
В данном разделе проведем второй этап решения поставленной задачи, заключающийся в использовании гипотезы Рейсса о том, что при простейшем нагружении в материале имеет место однородное напряженное состояние, то есть для любого ^го слоя (k = 1, п) [4, с. 5]:
агг,к =агг, аее,к =аее, аге,к =аге. (15)
Кроме того, из независимости решения от угла е следует, что ог0 = 0 . В полярной системе координат в случае плоского осесимметричного напряженного состояния уравнение совместности деформаций в терминах напряжений имеет вид [3, с. 90]:
даее „да гг
-у--и + у
-(1 + у)
( агг - а- ^
ее
V
дг дг
а уравнения состояния определяются обобщенным законом Гука [3, с. 32]:
-(8гг к + у8ее,к ),
= 0, (16)
Ек ,
агг = ^-2 (8ггк +У8й
1 -у
Ек -(8 + у8 ) (17)
\,ьее,^ уьггкЬ
аее = ,-218ее,к +у8
1 -у
где 8гг к, 8ее,к - нормальные компоненты деформаций к-го слоя (к = 1,п).
Разделив оба уравнения (17) на Ек, домножив их на ук (к = 1, п) и суммируя по к из (17), получаем:
1
Е "гг = 1-уг ((8гг)д + у'(8ее)я ),
Е "ее = ^ (^е)я +у'(8гг)я ),
я
1 ^ * +у./8 \ ) (18)
Л,2 \\-ее/я + ч \ гг/я'
— V
( п Л 1 п п
где {Е)я = Х^т- , (8гЛя =^Ук-8гг,к , (8ее)д =^Ук-8
„ 1ЕК
Так как при гипотезе Рейсса агг = (агг)я и аее = (аее)я, то для решения второй оценочной задачи остается расшифровать, что представляют собой (8гг)я и (8ее)д.
а2.
Подставляя в (17) выражения агг =1 —, аее = д-у, получим уравнение для определения
гг г дг дг
функции напряжений Эри ф [3, с. 90-91]:
д 3ф 1 д 2ф 1 дф Л
—Т +--2--Т^Т = 0, (19)
дг3 г дг 2 г 2 дг
решением которого является функция ф(г) = С1п(г) + С2 г2 + С3. Таким образом, напряжения имеют вид:
С С
агг = -| + 2С2, аее = 1 + 2С2. (20)
г г
При тех же, что и в первой задаче, граничных условиях агг | д = а^ и агг | ^ = а^ из (20)
получим [3, с. 92]:
С = (а я,-ая2 , (21)
Я2 - Я1
12
.2 л 2
" = 2 .(я12 - 2Я22) '
Из (15) и (17) следует, что для любых к и т (к,т = 1,п) выполнено:
= Jlv ^"--k +VS00-k) = YT^ (srr,m + vsee,m
E^L )_ ) (22)
CTee =~TE(see,t ) = -E7 (see,„ + vs„v„ )■
1 — v 1 — v
Решая уравнения (22) относительно srr,£ и see k, будем иметь:
E E
_ m _ m
S rr ,k = S rr ,m ' S66,i = ~E~ S66,m . (23)
Ek Ek
Кроме того, из (21) и (22) можно получить:
1 + V( ) R2 • R22 1 1 -V R>r -R22ar2
s„t=-(cr -cR -2T-- +---7^-2-(24)
E^ R RjR22 - R12 r2 Ek (R12 - 2R22) ' ' '
1 + v(c )R12• R22 1 + 1-v R12aR1 -R22cr2
Ek VR CrJR22 - R12 r2 Ek ' (R12 - 2R22) •
Домножив оба уравнения (24) на ук (к = 1, п) и суммируя по к, получаем искомые средние выражения для деформаций (егг^ и ^8ее):
/ \= 1+у( - ) Я* • К22 _1 +1—у Я а«! -Я2 аЛ2 (25)
= Е)к ^ а^Я22 - Я12 г2 + {Е)к " (я/- 2Я22) ' (25)
/8 \ = 1+у(а а ) *12 • *22 1 + 1 - у я - *2
* а* Я/ - Я12 г2 Е)К (я12 - 2Я22) .
Далее с помощью (18) вычисляем (агг)я и (аее)*:
/ \ =( - ) *12 • *22 _1 + 1-у я12ая22а*2
\агг/я 1а«1 ^ / „ 2 „ 2 „2 + ^ „
R22 - R12 r2 1 + v R - 2R22)
/ \ =-( - )R12 • R22 _! + R>R - R22aR2 GrJR22 - R12 r2 R2- 2R22) •
(26)
Решение краевой задачи для трубы из композиционного материала
Очевидно, необходимо предположить, что искомые средние деформации (ег^, ^8ее) и напряжения (агг), (аее^ определяются выражениями:
(8„) = а.(8гг)) +(1 -а).(8гг)к,
(8е^ = а^8е^) + (1 -а).(8ее) я,
К) = а-(а„)) +(1 -а).(а„)я, (27)
(аее) = а • (аее)) + (1 - а> (аее)*,
где а (0 <а< 1) - вещественная константа. Определяя среднее интегральное по возможным значениям а, можно получить из (27) простейшее выражение, называемое в теории композиционных материалов приближением Хилла:
(О =1 • [Юу + Юп), (Вее) = 1 • ((Вее)V + (еее)п),
2
(28)
Ю = 2 • ((агг )у + (аЛ ), (аее ) = 2 • ((аее )у + (аее )п ).
Устремляя П2 ^ да при ап = 0 , можно перейти к рассмотрению композиционной в среднем изотропной плоскости с отверстием, тогда для нее из (13), (14) и (25), (26) можно получить:
а
Е, Е.
я
•а
• г ч _ •
• г2 , (аее)V
Еу
•а
Ш
*1 „2
(агг)д = ап • ^Г, (аее)д = _^ • ТТ.
Д2
п' г2
п' г2
(29)
/ \ =/ \ =(1±^) / \ =/ \ =_(1±^)
\Вгг/у \ гг/п = Е 'СЯ1 • г2 , у =\Вее/я = ф ' г2 •
Подставляя (29) в (28), получаем, что эффективные деформации по Хиллу в композиционной плоскости с отверстием совпадают со средними по Фойгту и Рейссу, а эффективные напряжения зависят от отношения средних модулей упругости по Фойгту и Рейссу (рис. 2):
в„ =
(1+у) _ Я
Е/
п I \ (1 + у) а*1 • а
пи
п • г2
= 2'
1
. Е
а
г1 • „2 , (аее) = 2
1
п' г2
' Е 1
Е)
л
п У
а
• гт _ • -
п,2
п г2
(30)
Рис. 2. Непрерывное распределение эффективных напряжений (аее) по (30) в плоскости из композиционного линейноупругого материала с отверстием радиуса Е1 = 1 м в зависимости от концентраций компонент двухкомпонентного материала с давлением внутри отверстия
р = _ап = 106 Па с модулями упругости компонент Е1 = -1010 Па, Е2 = -109 Па при у1 = у, у2 = 1 — у : 1 - у = 0; 2 - у = 1/2; 3 - у = 1
2
2
1
Выводы
Впервые решена краевая задача для твердого композиционного тела без использования нелокальных гипотез о малости объема композиционного материала, для которого устанавливаются эффективные характеристики.
В связи с особенностями постановки задачи в случае композитного материала нет возможности отдельно рассматривать плоское напряженное состояние и плоскую деформацию поперечного сечения трубы, и оба этих состояния будут участвовать в оценке напряженно-деформированного состояния изучаемого объекта согласно стандартным гипотезам Фойгта и Рейсса.
Получено приближение Хилла для средних по представительному объему напряжений и деформаций.
В силу того что и напряжения и перемещения как на внутренней, так и на внешней границе постоянны, то и их средние значения по участку границы любой площади постоянны и равны исходным значениям. В связи с этим с методической точки зрения для задачи Ляме для трубы показано, что решения, построенные с использованием гипотез Фойгта и Рейса, самодостаточны и не требуют использования дополнительных предположений о малости элементов усреднения внутри упругого тела.
Установлено, что решение по напряжениям зависит от средних значений модулей упругости по Фойгту и Рейссу, а по деформациям определяется только средними значениями модуля упругости по Рейссу.
Получены формулы, определяющие напряженно-деформированное состояние композиционной в среднем изотропной плоскости с отверстием. Эти решения могут быть применены в качестве оценочных значений напряжении и деформаций в поперечном сечении грунта вокруг свай при, например, бурозабивном способе их погружения в вечной мерзлоте. В этом случае эффективный модуль упругости замерзшего грунта будет определяться как средний по Рейссу модуль упругости в отдельности каждой из компонент грунта и льда.
Методическое изложение решения задачи позволит как инженерам, так и научным работникам проводить собственные исследования в этом направлении.
СПИСОК ЛИТЕРА ТУРЫ
1. Андреев, В. И. Оптимизация по прочности толстостенных оболочек / В. И. Андреев, И. А. Потехин. -М. : Изд-во МГСУ, 2011. - 86 с.
2. Ватульян, А. О. Об определении внутреннего давления в цилиндре по данным акустического зондирования / А. О. Ватульян, В. В. Дударев // Дефектоскопия. - 2014. - N° 4. - С. 52-60.
3. Жемочкин, Б. Н. Теория упругости / Б. Н. Жемочкин. - М. : Госстройиздат, 1957. - 256 с.
4. Тарасюк, И. А. Сужение «вилки» Фойгта - Рейсса в теории упругих структурно неоднородных в среднем изотропных композиционных тел без применения вариационных принципов / И. А. Тарасюк, А. С. Кравчук // APRIORI. Серия: Естественные и технические науки. - 2014. - № 3. - Электрон. текстовые дан. - Режим доступа: http://apriori-joumal.ru/seria2/3-2014/Tarasyuk-Kravchuk.pdf (дата обращения: 16.11.2014). - Загл. с экрана.
REFERENCES
1. Andreev V.I., Potekhin I.A. Optimizatsiyapoprochnosti tolstostennykh obolochek [Optimization on the Strength of Thick-Walled Shells]. Moscow, MGSU Publ., 2011. 86 p.
2. Vatulyan A.O., Dudarev V.V. Ob opredelenii vnutrennego davleniya v tsilindre po dannym akusticheskogo zondirovaniya [On Measuring the Internal Pressure in the Cylinder (Based on the Data of Acoustic Sounding]. Defektoskopiya, 2014, no. 4, pp. 52-60.
3. Zhemochkin B.N. Teoriya uprugosti [Theory of Elasticity]. Moscow, Gosstroyizdat Publ., 1957. 256 p.
4. Tarasyuk I.A., Kravchuk A.S. Suzhenie «vilki» Foygta - Reyssa v teorii uprugikh strukturno neodnorodnykh v srednem izotropnykh kompozitsionnykh tel bez primeneniya variatsionnykh printsipov [Narrowing of Voigt-Reuss "Range" in the Theory of Elasticity of Structural Heterogeneous, Average Isotropic, Composite Bodies without Application of Variational Principles]. APRIORI. Seriya: Estestvennye i tehnicheskie nauki, 2014, no. 3. URL: http://apriori-journal.ru/seria2/3-2014/ Tarasyuk-Kravchuk.pdf. (accessed November 16, 2014).
STRESS-STRAIN STATE OF THICK-WALLED PIPE OF AN ELASTIC COMPOSITE IN AVERAGE ISOTROPIC MATERIAL
Aleksander Stepanovich Kravchuk
Doctor of Physical and Mathematical Sciences, Associate Professor, Professor of Department of Bio- and Nanomechanics, Belarus State University [email protected]
Prosp. Nezavisimosti, 4, 220030 Minsk, Belarus
Anzhelika Ivanovna Kravchuk
Candidate of Physical and Mathematical Sciences, Associate Professor,
Department of Web Technologies and Computer Simulation,
Belarus State University
Prosp. Nezavisimosti, 4, 220030 Minsk, Belarus
Abstract. For the first time the boundary value problem is solved for a composite solid body without the application of non-local hypotheses about the smallness of representative volume of composite material, for which the effective characteristics are set. It is assumed that the variation of Poisson's ratios for all the components of the composite material about the mean value is small, and therefore it is possible to use the condition of constancy of Poisson's coefficients and its equality to the average value of all components of the composite material. Due to the peculiarity of this problem, it is not possible in the case of the composite material to separately consider the plane stress and the plane strain state of the pipe, because both these states participate in the evaluation of the stress-strain state of the studied object, according to the standard hypothesis of Voigt and Reuss. An approximation for Hill average value of stresses and strains for representative volume are obtained. Due to the fact that stress and displacement both on the inner and outer border are constant, their average values at any section of the boundary area are constant and equal to their original values. In this connection, from the methodological point of view to the Lame problem, it is shown that solutions built using Voigt and Reuss hypotheses are self-sufficient and do not require the use of additional assumptions about the smallness of the averaging elements inside the elastic body. It was established that a solution for the stresses depend on the Voigt and Reuss mean values of the elastic moduli, and the deformations are determined only by the Reuss mean values of elasticity modulus. The formulas that determine the stress-strain state in composite average isotropic plane with the hole is obtained. These solutions can be used as estimates of stresses and strains in the cross section of soil around the pile using, for example, brown-pile method of their immersion in the permafrost conditions.
Key words: composite heterogenic material, discrete random variable, averaging, effective deformation characteristics, Voigt hypothesis, Reuss hypothesis, Hill approximation.