УДК 539.376
МЕТОД РЕШЕНИЯ ЗАДАЧИ О ЧИСТОМ СДВИГЕ СТОХАСТИЧЕСКИ НЕОДНОРОДНОЙ ПЛОСКОСТИ В УСЛОВИЯХ УСТАНОВИВШЕЙСЯ ПОЛЗУЧЕСТИ
Н. Н. Попов, О. О. Чернова
Самарский государственный технический университет,
443100, Россия, Самара, ул. Молодогвардейская, 244.
E-mails: ponick25@gmail.com, chernova_olga@citydom.ru
Разработан аналитический метод решения нелинейной задачи установившейся ползучести при чистом сдвиге стохастически неоднородной плоскости на основе второго приближения метода малого параметра. Введено ограничение о малости упругих деформаций, которыми допустимо пренебречь. Стохастич-ностъ введена в определяющие соотношения ползучести, взятые в соответствии с нелинейной теорией вязкого течения, при помощи случайной однородной функции координат. С использованием метода 'разложения компонент тензора напряжений по малому параметру до членов второго порядка малости получена система дифференциальных уравнений в частных производных относительно напряжений для первого и второго приближений. Решение системы строилось путём введения функции напряжений. Вычислены математическое ожидание и дисперсии случайного поля напряжений. Проведён анализ результатов, полученных в первом и втором приближениях.
Ключевые слова: чистый сдвиг, метод малого параметра, установившаяся ползучесть, второе приближение, стохастическая задача.
Наиболее часто используемым методом решения стохастических краевых задач является метод малого параметра, который позволяет свести статистически нелинейную задачу к последовательности статистически линейных задач. Следует отметить, что данный подход связан с трудностями вычислительного характера, поэтому при решении нелинейных стохастических задач обычно ограничиваются первым (реже вторым) приближением метода малого параметра. Метод малого параметра детально разработан лишь для стохастических задач линейной теории упругости В. А. Ломакиным [1,2]. Его применение в теории ползучести сдерживается проблемами не только статистической, но и физической нелинейности.
В данной работе приводится решение нелинейной задачи установившейся ползучести при чистом сдвиге стохастически неоднородной плоскости на основе второго приближения метода малого параметра. Аналогичная задача в корреляционном приближении (соответствует первому приближению метода малого параметра) на основе метода спектрального представления рассматривалась в работе [3]. Вводится ограничение о малости упругих деформаций и считается, что ими допустимо пренебречь. Деформируемый материал считается стохастически неоднородным, так что компоненты тензоров напряжений и скоростей деформаций установившейся ползучести в декартовой ортогональной системе координат являются случайными функциями
Николай Николаевич Попов (к.ф.-м.н., доц.), доцент, каф. прикладной математики и информатики.
Ольга Олеговна Чернова, младший научный сотрудник, каф. прикладной математики и информатики.
координат х\ и ^2-
Определяющие соотношения ползучести принимаются в соответствии с нелинейной теорией вязкого течения (установившейся ползучести) в стохастической форме [4]
Ргу = - ^5цСГкк^ (1 + аи (Ж1,Ж2)) , (1)
где в2 = 0,5 (Забегу — — интенсивность напряжений, р^ — компоненты тензора деформаций, —компоненты тензора напряжений, 5^ — сим-
вол Кронекера, II (х\, Жг) — случайная однородная функция, описывающая возмущение реологических свойств материала с математическим ожиданием (17) = 0 и дисперсией (и2) = 1; с, п — постоянные материала, а —число, определяющее степень неоднородности материала. По повторяющимся индексам производится суммирование от 1 до 2. В дальнейшем предполагается, что а— малый параметр, по которому производится разложение в ряд (метод малого параметра).
К определяющим соотношениям ползучести (1) присоединяются уравнения равновесия для напряжений и условие совместности деформаций, сформулированное для скоростей деформаций ползучести:
)3 = Ц] = 1)2, (2)
^~гj^~klPjk■>il = О- (3)
Здесь Л=(_° д ^ —единичный антисимметричный псевдотензор.
Соотношения (1)-(3) задают стохастическую задачу ползучести, которая решается относительно напряжений на основе второго приближения метода малого параметра.
После подстановки определяющего соотношения (1) в уравнение совместности (2) и сокращения на величину (82) (га 5)/2 в работе [5] было получено следующее уравнение:
(Т1 1 11 3 / о \ 2 ,. Гтх 12 2 ^ тт\
^ 2 2 ^ + ^ ^ 2 8 8,22 + +
+ (п - 1) 82822а11,2 + (з2)2 аи,22^ (2(711 - 022) +
+ 2 ^—---------82822 (1 + аи) + («2) а11,2^ (2(711 — а‘2‘2) ,2 +
+ (з2) (1 + а17) (2(7ц — (722),22 +
(Т1 1 77^ 3 / о \ 2 г. гтЧ ТЬ 12 2 г-г\
+ ^ 2 2 + 2 5 5,п (1 + +
+ (п - 1) 82821а11,1 + (82)2аС/,11^ (2(722 - 0-ц) +
+ 2 ^(1 + СкС/) + (з2) а11,1^ (2(722 — С7ц) , 1 +
+ (в2)2 (1 + аП) (2(722 - он),и -— 3(п — 1) ^ - 82\822 (1 + <У.и) + 828212 (1 + а,и) +
+ 82822а11,1 + 82821а11,2^(712 - 6(82)2а11,12(712-
- 6 (^7^82821 (1 + аи) + (82)2аи, 1^ (712,2 —
- б(^7^82822 (1 +аи) + (82)2а11,2^ (712,1 “
- 6 (52)2(1 + аС/) (712,12 = 0. (4)
Для приближённого решения уравнения (4) используется метод разложения компонент тензора напряжений и интенсивности напряжений в2 по
малому параметру а до членов второго порядка малости:
"° + аа^ + а2сг2^, '-----0
„2 „2 , , „2
аЧ — аіі + а<тЬ + а2 (агз) — (5)
82 = 8^ + а8і + а2 82- (6)
В силу того, что в условии чистого сдвига детерминированные нормальные напряжения (7і! и а22 равны нулю, а касательное напряжение а^2 является постоянным, имеем
4 = 3(Т12; «1 = 6<Ті2а12’ «2 = (^їі)2 + {а22}2 ~ °11°22 + Кап)2 + 6а°12а\2.
Уравнение (4) путём подстановки в него представлений (5), (6) и прирав-
нивания коэффициентов при одинаковых степенях а приводится к системе двух уравнений:
(2(71! - °22),22 + (20-22 “ ^Іі) ,11 “ ^2,12 ~ бЬ^іг^іг = 6[/, 12(7?2, (7)
(2(711 — ^22) 522 + (2(7^2 — — ^а12,12 ~ , 12СГ12 —
= 6к(}81,181,2(7^12 + 6к (^1,12^2 + «1,1^12,2 + «1,2^12,1~
_ 5151,*2°~12 +§1,211,1(712 + 51,1и,2(Т01^\ -Л(51,22(2(7111-(722)+81,11(2(722-(7111) + 50 /
+281,2(2а\1-а12) ,2+^1,1(2а12-(7п),1) -и,22(2а11-а12)-^и,2(2а11-а12) ,2~
-и, 11 (2(722 -^и) -2С/, 1(2(722 ~(7п),1 + 611,12(7^ ~ 6С/С/, 12(7?2, (8)
где к = (п — 1)/(2в§), д = (п- 3)/(2в§).
Подставив в уравнения (7), (8) приближения «1, «2 интенсивности напряжений з2, получим систему дифференциальных уравнений в частных производных относительно (7^- И (7г2:
2°П,22 — °22,22 + 2СГ22,11 — ^ЬИ — = 6(7°2^, 12, (9)
^а11,22 а22,22 + ^а22,11 а11,11 ^Па12,12 — ‘^‘^^(1а12,1а12,2 {а12) +
+ ЗбА:СГ°2 (^а12,1а12,2 ~ а12а12,12 + а 12,2^,1^ 12 + сг12,1^?2СГ12) —
- 6ка°12 (а\2
,22(2СГ11 — СГ22) + ^12,и(^а22 ~ ^11) + ^а12,2(^а11,2 ~ °22,2) +
+ ^а12,1(^а22,1 ~ а11,1) ~ Г? 12) “ и,22(^11 ~ а22) ~ 2^,2(2сг{1?2 “ ^22,г)-- и,и(2(722 - °п) - 2С/,1(2сг22, 1 - ^11,1) + 6С7,12(712 - еии,12<т°12, (Ю)
где Г = (<7^)2 + ((7^2)2 - <7п<722 + З^)2'
К полученным уравнениям необходимо добавить уравнения равновесия для приближений напряжений
<„=0, 171 = 1,2. (11)
Если ввести функцию напряжений Р для приближений тензора напряжений с помощью формул
аг]1] = — ^1,у, = 5цАР2 — ^2,4 (12)
(А^т = -Рт,ц+-Рт,22 — оператор Лапласа), то вместо системы четырёх уравнений (9)—(11) относительно напряжений получим два линейных уравнения относительно И ^
^1,1111 + Р1,2222 + (3П — 1)^1,1122 = За^и,12, (13)
,1111 + ^2,2222 + (Зп — 1)^2,1122 = ^(-^Ъ Ю- (14)
Правая часть у?(^1,С7) уравнения зависит от производных функций и II. В силу громоздкости здесь она не представлена.
Однородную функцию II (Х1,Х2), с помощью которой даётся случайное поле возмущений реологических свойств материала, будем брать в виде [6]
и(Х1,Х2) = Ас08(Ш1Х1 + Ш2Х2 + <р), (15)
где А — центрированная случайная величина, ср — случайная величина, имеющая равномерное распределение на интервале (0;27г), причём (А^) = 0, и)1 и Ш2 — положительные неслучайные параметры.
Для удобства выкладок целесообразно перейти к функции комплексной переменной
и = Ае^Ш1Х1+Ш2Х2), (16)
где А = Аег{р. Функция II введена так, что Лей = 11.
Будем рассматривать уравнения (13), (14), в которых действительные
функции II, заменены комплексными II, причём Re.Fi =
Решение уравнения (13) имеет вид
Рг = А/е^Ш1Х1+Ш2Х2), (17)
где / — постоянная, найденная путём подстановки представлений (16), (17) в уравнение (13):
. ___________ЗШ1Ш2СГ12_______
ш\ + ш2 + (3 П — 1)(л)1и)2
Выделяя действительную часть имеем
= Г1е Р1 = А/ сов^Ж! + Ш2Х2 + ф) или с учётом формулы (15) — = /[/(ж ьж2).
Используя соотношение (12), решение задачи в первом приближении в напряжениях можно записать в виде
(1) ! 2 £тт 3аш1ш1а\2и
= а<тц = -аш2/и = , 4 4 —
6^2 (3Т1 — 1)бс,^Сс,2
(1) 1 2 tTT Заш\ш2а°12и
°22 = aa22 = -CX^ifU = -f-------------Г—--------- 2 2 ’ 18
+ ш2 + (3 п — 1)ш(ш2
m n 1 n п 3au}^u}ocr%U
*12 = <г12 + а<г12 = <г12 + оиши = <г12 - ^+i,, + (3n_lKll,r
Подставляя решение (17) в правую часть уравнения (14), можно получить
^2,1111 +-^2,2222 + (Зп — 1)^2,1122 = 3 А2Ъе2%<уШ1Х1+Ш2Х2\ (19)
где 6 = —2(36kqf2ujfuj2(a(l2yi+4:kf2ujiuj2(uj'f+uj2)2a(l2+6kf2ujiuj2(ujf+uj2)a(l2 + + 12kfu)fu)%)2a°12 + (u)j + u%)f - wiw2<7?2) •
Решение уравнения (19) можно искать в виде
р2 = А2де2^Ш1Х1+Ш2Х2\ (20)
д = const. После подстановки (20) в (19) для нахождения д получается алгебраическое уравнение, из которого следует
9 =
36
16 (cof + ш2 + (Зп — 1)ш2ш2)
Согласно (20) решение уравнения (19) определяется формулой
2^е2г(ш1х1+ш2Х2)
F2 =
16 (ujf + ш2 + (3 п — l)u)fu)2)
Выделяя действительную часть Р2 и используя соотношения (12), получим вторые приближения для компонент тензора напряжений:
<^11
ЗА2Ьш2 cos(2w^i + 20J2X2 + ф)
4 (uif + ш2 + (Зп — l)ojfoj2)
2 _ 3A2bujf cos(2u)iXi + 2ш2Х2 + (fi) , .
022 4(ujf + ш2 + (3 n — 1 )oj2oj2)
2
al2
3A2bwiUJ2 cos(2w^i + 20J2X2 + ф)
4 (cvf + uj2 + (3 n — l)ojjoj2)
Решение задачи (1)—(3) во втором приближении согласно (5), (18), (21) определяется формулами
(2) Заш\uj'^a^U 3A2bw% cos(2wia;i + 2^2X2 + <p)
°11 Ulf + U>2 + (з n — l)wfw2 4 (cvf + U>2 + (Зп — l)wfw2) '
(2) 3aojfoj2<Ti2U 3A2bw2 cos(2wiX\ + 2^2X2 + <p)
022 Ulf + U>2 + (3n — l)wf^2 4(w^ + + (3n — l)u)fu>2)
(2) о Заш2^^^ 3A2b(Vi(V2Cos(2u)iXi-h 2u)2X2 + <fi)
Gl2 Gl2 ojf + U>2 + (3n — 1)W2W2 4(wf + + (3n — 1)uj2uj2)
Найдём основные статистические характеристики случайного поля напряжений: математическое ожидание и дисперсии компонент тензора напряжений. В силу того, что по условию задачи (U) = 0, имеем (сг^) = (of,-)-
Покажем, что (of-) = 0. Для этого согласно (21) достаточно установить, что {A2 cos(2u)iXi + 2u)2X2 + ф) = 0. В силу независимости случайных величин Avup имеем (cos(2wia;i + 2Ш2Х2) +<р) = cos(2wiX\ + 2Ш2Х2) (cos Lp) + sin(2wia;i + + 2u)2X2){sm(p) = 0. Поскольку случайная величина Lp распределена равномерно на отрезке [0;27г], имеем (cos Lp) = 0 (см. [6]). Отсюда (of-) = 0. Таким образом, математические ожидания в первом и втором приближениях равны (7°-.
Дисперсии напряжений в первом приближении с учётом условия {U2) = 1 определяются формулами
D
(і)__ і ( і\2\ _ 9а2ш2шЦа°12)
И - / - / 4 4 , , 2 2\2’
[0^1 (ЗТІ — 1)Сс,^Сс,2^
(і) її 1 ч2Ч 9а2шІиі2 (^12)2
22 ~ ^ ~ (ш}+^ + (Зп-1)ш?ш|)2’
й<‘> = (И2)2) = ----------(,Т?2)2-------------2 ■
(ио\ + ио2 + (Зп — 1)ш2(х;|)
Дисперсии для второго приближения можно записать в виде
Ві? = В [°І] + °^В \аЬ] + 2«3 )• (23)
В дальнейшем предполагаем, что случайная величина А распределена по нормальному закону. Для нормально распределённой случайной величины все центральные моменты нечётных порядков равны нулю, а момент четвёртого порядка {АА) = 3{А2) = 6 [7]. Момент (<т^-сгг2-), входящий в (23), равен нулю, так как он пропорционален моменту третьего порядка (А3). Дисперсия .О[а?-], как показано в работе [8], пропорциональна величине 0,5(А4).
Таким образом, дисперсии напряжений во втором приближении согласно (22), (23) определяются формулами
(2) _ 9а2ш\ш1 (<т?2)2 27а4Ь2ш% (а°12)2
Ї-'ЛЛ — , . . ... О “Г
(ujf + LUrj + (3n — l)^2^!)2 16 + LOrj + (3n — l)^2^!)2
9a2cofco% ((Т12)2 27aAb2ujf
jj(2) _ ________~ ~~ - 1- ^ 1-4)_________|___________________________________
22 (wf + W2 + (Зп — l)^2^!)2 16 (uj\ + OJ2 + (3n — l)^2^!)2
jj(2) _ _____~~~ "1"^ 1^/______I_____________________________
12 (wf + + (3n — l)^2^!)2 16 (uj\ + Wg + (3n — l)^2^!)2
Анализ случайного поля напряжений был проведён при условии, что частоты флуктуаций микронеоднородностей и)\ и W2 относительно осей Х\ и Х2 равны между собой (00\ = U)2 = со). При этом условии случайное поле U, заданное разложением (15), можно считать близким к изотропному [1]. Дисперсии напряжений и в первом, и во втором приближениях будут равны между собой (-Оц = D22 = D12 = D) и не будут зависеть от со. Тогда формулы для вычисления коэффициента вариации будут иметь вид
= = d = = л/За /
<7° Зп + 1 ’ (7° 4(3п + 1)
где (7° = (7?2-
Значения коэффициента вариации d\ и с?2 (в процентах) как функции переменных а и п приведены в табл. 1. В столбце d\ представлены значения коэффициента вариации, вычисленные по первому приближению, в столбце с?2 — по второму приближению. Как видно из табл. 1, для материалов с высоким показателем нелинейности (п = 9) коэффициент вариации с?2 изменяется в пределах от 1,08% (а = 0,1) до 5,81% (а = 0,5). В случае низких показателей нелинейности, когда возможна полная физическая линеаризация закона ползучести (п = 1), разброс напряжений около среднего значения значительно больше: здесь значения коэффициента вариации заключены в пределах от 7,54% (а = 0,1) до 42,64% (а = 0,5).
В табл. 2 приведена относительная погрешность вычисления среднеквадратичного отклонения второго приближения относительно первого:
у/]№-у№)
5 =-------==-------100%.
Из табл. 2 видно, что имеется довольно обширная область значений параметров «пи, для которой погрешность 5 < 5% и первое (линейное) приближение вполне приемлимо для решения практических задач. Второе приближение необходимо использовать для материалов с большой степенью неоднородности (а = 0,4 ч- 0,5), особенно при малых значениях п.
9a2cofco^ (o"i2)2 27аАЪ2ио\ио2 {^и)2
Таблица 1
Значения коэффициента вариации d\ и с/2 (в процентах)
п 1 3 5 7 9
а di d2 di d2 di d2 di d2 di d2
0,1 7,50 7,54 3,00 3,01 1,88 1,88 1,36 1,37 1,07 1,08
0,2 15,00 15,35 6,00 6,09 3,75 3,80 2,73 2,77 2,14 2,17
0,3 22,50 23,66 9,00 9,29 5,63 5,80 4,09 4,22 3,21 3,32
0,4 30,00 32,69 12,00 12,67 7,50 7,91 5,45 5,75 4,29 4,52
0,5 37,50 42,64 15,00 16,29 9,38 10,17 6,82 7,40 5,36 5,81
Таблица 2
Относительная погрешность вычисления среднеквадратичного отклонения (в процентах)
n a 1 2 3 4 5 6 7 8 9
0,05 0,146 0,097 0,090 0,088 0,088 0,088 0,088 0,089 0,089
0,10 0,581 0,388 0,358 0,351 0,350 0,351 0,352 0,353 0,354
0,15 1,293 0,867 0,801 0,785 0,783 0,784 0,786 0,789 0,792
0,20 2,264 1,526 1,410 1,383 1,379 1,381 1,385 1,390 1,395
0,25 3,472 2,355 2,178 2,137 2,130 2,133 2,139 2,147 2,154
0,30 4,890 3,339 3,092 3,034 3,024 3,029 3,038 3,049 3,059
0,35 6,487 4,465 4,140 4,064 4,051 4,057 4,069 4,083 4,096
0,40 8,234 5,717 5,308 5,212 5,196 5,204 5,219 5,236 5,253
0,45 10,100 7,078 6,581 6,465 6,445 6,455 6,473 6,494 6,514
0,50 12,056 8,531 7,945 7,808 7,785 7,796 7,818 7,842 7,867
Работа выполнена при поддержке РФФИ (проект № 10—01 —00644—а).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Ломакин В. А. Статистические задачи механики твёрдых деформируемых тел. М.: Наука, 1970. 137 с. [Lomakin V. A. Statistical problems of the mechanics of deformable solids. Moscow: Nauka, 1970. 137 pp.]
2. Ломакин В. А. Проблемы механики структурно-неоднородных тел // Изв. АН СССР. МТТ, 1978. №6. С. 45-52. [Lorn,akin V.A. Problems in mechanics structurally nonuniform bodies// Izv. Akad. Nauk SSSR, Ser. Mekh. Tverd. Tela, 1978. no. 6. Pp. 45-52].
3. Попов H. П., Яшин М. А. Исследование случайных полей напряжений при чистом сдвиге стохастически неоднородной полуплоскости в условиях ползучести // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2010. №1(20). С. 104-110. [Popov N. N., Yashin М. A. The study of random fields of stress in pure shear stochastically inhomogeneous halfplane under creep // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. nauki, 2010. no. 1(20). Pp. 104-110].
4. Радченко В. П., Попов Н.Н. Стохастические характеристики полей напряжений и деформаций при установившейся ползучести стохастически неоднородной плоскости // Изв. вузов. Машиностроение, 2006. №2. С. 3-11. [Radchenko V. P., Popov N. N. Stochastic characteristics of stress and strain fields in steady-state creep of stochastically inhomogeneous plane// Izv. Vuzov. Mashinostroenie, 2006. no. 2. Pp. 3-11].
5. Попов H. П., Коваленко Л. В., Яшин М. А. Решение плоской нелинейной стохастической задачи ползучести методом спектральных представлений // Вестн. Сам. гос. техн. унта. Сер. Физ.-мат. науки, 2009. №2(19). С. 99-106. [Popov N. N.,Kovalenko L. V., Yashin M. A. Solution of plane nonlinear stochastic problem with spectral representation method // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. nauki, 2009. no. 2(19). Pp. 99-106].
6. Вентцель E. С., Овчаров Л. А. Прикладные задачи теории вероятностей. М.: Радио и связь, 1983. 416 с. [Wentzel Е. S., Ovcharov L. A. Applied Problems in Probability Theory. Moscow: Radio i Svyaz’, 1983. 416 pp.]
7. Пугачев В. С. Теория вероятностей и математическая статистика. М.: Физматлит, 2002. 496 с. [Pugachev V. S. Theory of Probability and Mathematical Statistics. Moscow: Fizmatlit, 2002. 496 pp.]
8. Попов H.H., Чернова О. О. Решение нелинейной задачи ползучести для стохастически неоднородной плоскости на основе второго приближения метода малого параметра // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2011. №4(25). С. 50-58. [Popov N. N., Chernova О. О. Solution of nonlinear creep problem for stochastically inhomogeneous
plane on the basis of the second approximation for small parameter method // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. nauki, 2011. no. 4(25). Pp. 50-58].
Поступила в редакцию 01/XI/2012; в окончательном варианте — 19/XI/2012.
MSC: 35Q74; 74Е35, 74К20
METHOD OF PURE SHEAR PROBLEM SOLVING FOR STOCHASTICALLY INHOMOGENEOUS PLANE IN A STEADY-STATE CREEP
N. N. Popov, O. O. Chernova
Samara State Technical University,
244, Molodogvardeyskaya St., Samara, 443100, Russia.
E-mails: ponick25@gmail.com, chernova_olga@citydom.ru
The analytical method for nonlinear problem of steady-state creep solving for pure shear of stochastically inhomogeneous plane on the basis of the second approximation method of small parameter was developed,. It is supposed that elastic deformations are insignificant and they can be neglected. Stochasticity was introduced into the determinative creep equation, which was taken in accordance with the nonlinear theory of viscous flow, through a homogeneous random, function of coordinates. By using the decomposition technique of stress tensor components in a small parameter to the members of the second order of smallness, partial differential system of the first and, the second approximation of stress was obtained,. This system was solved by the introduction of the stress function. The mathematical expectation and, variances of the random, stress field were calculated. The analysis of the results in the first and, second approximations was obtained,.
Key words: pure shear, small parameter method, steady-state creep, second approximation, stochastic problem.
Original article submitted 01/XI/2012; revision submitted 19/XI/2012.
Nikolay N. Popov (Ph.D. (Phys. & Math.)), Associate Professor, Dept, of Applied Mathematics & Computer Science.
Olga O. Chernova, Junior Research Scientist, Dept, of Applied Mathematics & Computer Science.