УДК 539.376
ПОСТРОЕНИЕ АНАЛИТИЧЕСКОГО РЕШЕНИЯ ПЛОСКОЙ СТОХАСТИЧЕСКОЙ НЕЛИНЕЙНОЙ КРАЕВОЙ ЗАДАЧИ УСТАНОВИВШЕЙСЯ ПОЛЗУЧЕСТИ С УЧЁТОМ ГРАНИЧНЫХ ЭФФЕКТОВ
Л. В. Коваленко, Н. Н. Попов
Самарский государственный технический университет,
443100, Самара, ул. Молодогвардейская, 244.
E-mails: flytitmouse@mail.ru; ponick25@gmail.com
Приводится решение нелинейной стохастической краевой задачи ползучести тонкой пластины при плоском напряженном состоянии при условии, что упругие деформации малы и ими допустимо пренебречь. Определяющее соотношение ползучести, взятое в соответствии с нелинейной теорией вязкого течения, сформулировано в стохастической форме. При помощи метода малого параметра нелинейная стохастическая задача сводится к системе трех линейных уравнений в частных производных относительно флуктуаций тензора напряжений.
Эта система при помощи перехода к функции напряжений была сведена к одному дифференциальному уравнению, решение которого представлено в виде суммы двух рядов. Первый ряд задает решение вдали от границы тела без учета краевых эффектов, второй ряд представляет собой решение в пограничном слое, его члены быстро затухают по мере удаления от границы пластины. На основе полученного решения проведен статистический анализ случайных полей напряжений вблизи границы пластины.
Ключевые слова: установившаяся ползучесть, стохастически неоднородная пластина, случайное поле напряжений, метод малого параметра, краевой эффект.
Рассматривается стохастическая краевая задача о нелинейной ползучести тонкой пластины при плоском напряженном состоянии при условии, что упругие деформации малы и ими допустимо пренебречь. Материал пластины считается стохастически неоднородным, так что тензоры напряжений и скоростей деформаций являются случайными функциями координат xi и Х2 в декартовой ортогональной системе координат.
Нелинейная стохастическая краевая задача ползучести состоит из дифференциальных уравнений равновесия для напряжений
Людмила Викторовна Коваленко (к.ф.-м.н.), ассистент, каф. прикладной математики и информатики. Николай Николаевич Попов (к.ф.-м.н., доц.), доцент, каф. прикладной математики и информатики.
aij,j = 0, i = 1, 2; j = 1, 2, условия совместности деформаций
ЛІ7 Лк1 Pjk,il = ° детерминированных граничных условий
(1)
(2)
®ij nj IГ — qi
(3)
и определяющих соотношений ползучести, взятых в соответствии с нелинейной теорией вязкого течения в стохастической форме [1]
В формулах (1)—(4) о^ — компоненты тензора напряжений; — компоненты тензора деформаций; Л^ —единичный антисимметричный псевдотензор; 5гз — символ Кронекера; и (ж1, х2) — случайная однородная функция, описывающая флуктуации реологических свойств материала с математиче-
териала. Постоянная а определяет степень неоднородности материала и для реальных материалов может изменяться в пределах от 0,05 до 0,5 [2]. Точка обозначает дифференцирование по времени. По повторяющимся индексам производится суммирование от 1 до 2.
Путём разложения компонент тензора напряжений о^ по малому параметру а
нелинейная стохастическая задача (1), (2), (4) может быть приведена к системе линейных дифференциальных уравнений в частных производных:
В силу трудностей математического характера при решении системы (5), как и в работах [3-5], будем ограничиваться первым приближением метода малого параметра, которое справедливо для слабо неоднородных сред. Путём введения функции напряжения ¥ по формулам (при т = 1)
вместо системы уравнений (5) можно получить единственное дифференциальное уравнение относительно ¥:
(2 + q^) + 2 (2 + qhl2) 1122 + (2 + q/f) ¥,2222 —
= —a (liU,22 + I2U,ii) (7)
ским ожиданием (U) — 0 и дисперсией (U2) — 1; c, n, a — постоянные ма-
m=i
(5)
— (n — !)/2Sq; г ' —некоторая величина, зависящая от решений на предыдущих шагах. При выводе системы (5) предполагалось, что a(Q) — const, a(2 — 0. Краевые условия для системы (5) с учётом (3) имеют вид aj nj |г— 0.
(6)
с граничными условиями
(^,22^1 — ^,12П2)|Г = 0, ( — ^,12^1 + Р, 11 п2)1Г = 0. (8)
Пусть однородная функция и(^1,^2), описывающая стохастические свойства материала, является почти периодической быстро осциллирующей функцией координат [6]:
U = ^ Ak cos (wckx1 + wdkx2 + pk),
k=i
где ш — большой параметр, имеющий размерность, обратную длине; Ck, dk — безразмерные величины порядка единицы; Ak — центрированные одинаково распределённые случайные величины, рk — случайные величины, распределенные равномерно на интервале (0;2п), причём все величины Ak и pk независимы.
Для удобства выкладок целесообразно перейти к функции комплексного переменного:
ГО
U — ^2 Ak exp(iw(ckxi + dkX2)), Ak — Ak exp ipk. (9)
k=i
Функция U введена так, что Re U — U. Тогда решение краевой задачи (7), (8), в котором U и F заменены величинами U и F (Re F — F), можно представить в виде
F = ^(vk + Wk), (10)
k=i
где Ук — частное решение уравнения (7), полученное при замене функции и к-тым членом разложения (9), а Wk — решение соответствующего (7) однородного уравнения, удовлетворяющее на границе Г условиям:
(vk,22n1 — vk, 12n2)|r = (wk,22n1 — wk, 12n2) |г, (vk, 11n2 — vk, 12n1)|r = (wk, 11n2 — wk, 12n1) 1 p.
(11)
Разыскивая функции у к в виде
Ук = /к ехр(ш(0кХ1 + dкЖ2)), (12)
для /к можно получить следующее выражение:
/ = _______Ак (11 ^к + 12с|)________
к ш2(2(ск + ¿к)2 + я(1^к + кск)) ’
Ряд £^=1 Wk задаёт при ш ^ 1 решение типа пограничного слоя, быстро затухающее по мере удаления вглубь области тела. Построим решение Wk
типа пограничного слоя в области Х2 ^ —Ь вблизи границы Х2 = —Ь тела. Делая замену
Wk = дк (¿)ехр[гш(ск Х1 — ¿к Ь)], Ь = ш(х2 + Ь), (13)
из однородного уравнения, соответствующего (7), для функции дк(Ь) можно получить
(2 + Ч1\) — 2ск(2 + ~-дТ + ск (2 + ?12)дк = 0.
Решение этого уравнения, когда все корни т'к соответствующего характеристического уравнения простые, определяется формулой
дк(Ь) = ^2 Ска ехр (гкЬ), (14)
8=1
где Сак — произвольные постоянные.
Из четырёх корней характеристического уравнения два корня т3к и т\ имеют положительные действительные части. Так как при Х2 ^ то краевой эффект должен затухать, то две постоянные, отвечающие этим корням, равны нулю. Для нахождения двух других констант используются граничные условия
-дк
-Ь
9k |t=0 — fk, —^Г
— idk fk,
t=0
полученные из условий (11) с использованием выражений (12), (13).
В качестве примера рассмотрим ползучесть стохастически неоднородной полуплоскости X2 ^ 0, находящейся в условиях плоского напряженного состояния. Пусть к границе полуплоскости X2 — 0 приложены нагрузки
а22 Ц=0 — а2? — const, ai2 |x2=0 — °,
а напряжение ац удовлетворяет условию макроскопической однородности {ац) — а^0 — const, которое соответствует приложению нагрузки на бесконечности (xi —^ ±то).
Решение уравнения (14) типа пограничного слоя при условии а^ — а22 имеет вид
9k(t) — erk‘ - (15)
k k - k k r2 ri r2 ri
где
\ ±ai + v/ai + a2 rk,2 — -Ck(Bi ± iB2), Bi,2 — ^----------------------------,
2 + qlil2 y/2q\li — ¿21
ai — 2 + ql2 , a2 — 2 + ql2 '
гг (0) (0) (о)
При условии аЦ = a22 = a(0) корни характеристического уравнения являются кратными. В этом случае решение задачи приведено в [7] и здесь не рассматривается.
Подставляя выражения (12), (13), (15) в соотношение (10) и выделяя действительную часть F, можно получить
ОО
a
F = ReF = V' AkMk cos (ckwxi + Pk)
Ш2
cos dk wx2—
k=l
e-Ck В1ШХ2
(B2 cos CkВ2ШХ2 + Bi sin CkВ2ШХ2)
CkBi^
В 2 Ck
B2
/ g—CkBi^X2 4
— sin (ckwx1 + pk ) ( sin dkwx2--------------dk sin ckB2wx2 )
V B2 Ck )
где Мк = fкш2/Ак.
Компоненты тензора напряжений можно найти из выражения (15), используя формулы (6). В силу громоздкости они здесь не выписаны.
Вычисление дисперсий случайного поля напряжений 2) про-
изводилось при условии, что все величины Ск и dk равны 1. При этом условии случайное поле и, заданное разложением (9), можно считать близким к изотропному [6]. С учётом условий, наложенных на случайные величины Ак и Рк, и равенства (и2) = 1 дисперсии случайного поля напряжений определяются следующими выражениями:
D11 =Л
, g—Biux2 ч 2
^- cos wx2 +-----B-----(B2 + B|)(B2 cos B2Wx2 — Bi sin B2Wx2)j +
g—В1ШХ2 x 2i
g " ~2 B2
g— i 2
+ (^sin dkwx2 +-----B---((B2 + B|) sin B2Wx2 — 2B1B2 cos B2Wx2)j
D22 = Л
, g—В1ШХ2 4 2
(cos wx2------- ---(B2 cos B2wx2 + B1 sin B2wx2)) +
B2
g—Biux2 л 2!
g— i 2 2 + (^sin wx2 +---B------sin B2Ux2)
D12 = Л
, g—Biux2 \ 2
^— sin wx2 +----B-----(B2 + B|) sin B2Wx2) +
, g—BiuiX2 >2i
+ (cos wx2 +------- ----(B1 sin B2wx2 — B2 cos B2ux2))
B2
где
Л
B2
8a2s0(l1 +12)2
(1652 + (п - 1)(11 + Ь)2)2'
На основе полученных аналитических решений проведен статистический анализ случайного поля напряжений, в результате которого установлено, что
вблизи границы полуплоскости существует узкий пограничный слой, в котором разброс напряжений намного больше, чем для глубинных слоев.
На рис. 1 приведен типичный график нормированных дисперсий напряжений (х2)/Б^ (то) в зависимости от безразмерной координаты шх2.
Здесь Н = /^и ~ параметр нагружения.
С ростом шх2 дисперсии довольно быстро приближаются к постоянным значениям, совпадающим с их значениями для неограниченной среды. Концентрация напряжения, возникающая на границе полуплоскости за счёт неоднородности материала, вычисляемая по формуле
Р = / Бц(0)/Бц (то) = \! 1 + а\ + 2а\ + ,
может быть в два и более раз больше, чем на бесконечности, в зависимости от параметра нагружения и степени неоднородности материала.
В качестве второго примера рассматривалась ползучесть стохастически неоднородной полосы с быстро осциллирующими реологическими свойствами. Считается, что бесконечная стохастически неоднородная полоса (пластина) —то < XI < то, —Ь ^ Х2 ^ Ь растягивается вдоль оси XI постоянными напряжениями ац = а0, которые приложены на бесконечности (х1 ^ ±то), а её границы Х2 = ±Ь свободны от напряжений:
а22 |ж2=±Ь = 0, а12 |Ж2=±Ь = 0.
На рис. 2 представлены зависимости нормированных дисперсий = = (шх2)/В^(то) от безразмерной координаты шх2 при степени нелиней-
ности установившейся ползучести п = 3. Видно, что на границах полосы
О 2 4 6 8 шх2
Рис. 1. Нормированные дисперсии Йу при п = 3, Н = 2: 1 — О01, 2 — Д02, 3 — Д02
Рис. 2. Графики нормированных дисперсий D0 при n = 3: 1 — 2 — D02, 3 — В®2
Х2 = ±шЪ разброс напряжений стц имеет наибольшее значение, здесь D01=1,5.
Дисперсии D02, D02 принимают максимальные значения в пограничном слое.
Вдали от границ полосы дисперсии принимают те же значения, что и для
неограниченной плоскости.
Таким образом, флуктуации напряжений в пограничном слое играют существенную роль, и их необходимо учитывать при решении вопроса о надёжности элементов конструкций по критерию длительной прочности.
Работа выполнена при поддержке РФФИ (проект № 10-01-00644-a).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кузнецов В. А. Ползучесть стохастически неоднородных сред в условиях плоского напряжённого состояния/ В сб.: Математическая физика: Сб. научн. трудов. Куйбышев: КПтИ, 1977. С. 69-74. [Kuznetsov V. A. Creep of stochastically nonuniform media under conditions of a plane stress state/ In: Mathematical Physics (collected scientific papers). Kuibyshev: KPtI, 1977. Pp. 69-74].
2. Попов Н. Н., Самарин Ю. П. Пространственная задача стационарной ползучести стохастически неоднородной среды// ПМТФ, 1985. №2. С. 150-155; англ. пер.: Popov N.N., Samarin Yu. P. Spatial problem of stationary creep of a stochastically inhomogeneous medium// J. Appl. Mech. Tech. Phys., 1985. Vol. 26, no. 2. Pp. 296-301.
3. Ломакин В. А., Шейнин В. И. Концентрация напряжений на границе случайнонеоднородного упругого тела// Изв. АИ СССР. МТТ, 1974. №2. С. 65-70. [Lomakin V. A., Sheinin V. I. Stress concentration at the boundary of a randomly inhomogeneous elastic body // Izv. AN SSSR. MTT, 1974. no. 2. Pp. 65-70].
4. Наумов В. Н. Напряженное состояние случайно-неоднородного упругого полупространства// Изв. АН СССР. МТТ, 1976. №2. С. 58-63. [Naumov V. N. Stressed state of a randomly inhomogeneous elastic half-space// Izv. AN SSSR. MTT, 1976. no. 2. Pp. 58-63].
5. Попов Н. Н., Коваленко Л. В. Поля напряжений на границе стохастически неоднородной полуплоскости при ползучести // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2006. №42. С. 61-66. [Popov N.N., Kovalenko L. V. Stress fields on the boundary of a stochastically inhomogeneous half-plane with creep // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2006. no. 42. Pp. 61-66].
6. Ломакин В. А. Статистические задачи механики твёрдых деформируемых тел. М.: Наука, 1970. 139 с. [Lomakin V. A. Statistical problems of the mechanics of solid deformable bodies. Moscow: Nauka, 1970. 139 pp.]
7. Попов Н. Н., Самарин Ю. П. Исследование полей напряжений вблизи границы стохастически неоднородной полуплоскости при ползучести // ПМТФ, 1988. №1. С. 159-164; англ. пер.: Popov N. N., Samarin Yu. P. Stress fields close to the boundary of a stochastically inhomogeneous half-plane during creep// J. Appl. Mech. Tech. Phys., 1988. Vol. 29, no. 1. Pp. 149-154.
Поступила в редакцию 16/XII/2010; в окончательном варианте — 14/II/2011.
MSC: 35Q74; 74E35, 74K20
CONSTRUCTION OF ANALITICAL SOLUTION OF 2D STOCHASTICALLY NONLINEAR BOUNDARY VALUE PROBLEM OF STEADY CREEP STATE WITH RESPECT TO THE BOUNDARY EFFECTS
L. V. Kovalenko, N. N. Popov
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russia.
E-mails: flytitmouse@mail.ru; ponick25@gmail.com
The solution of nonlinear stochastically boundary value problem of creep of a thin plate under plane stress is developed,. It is supposed that elastic deformations are insignificant and they can be neglected.. Determining equation of creep is taken in accordance with nonlinear theory of viscous flow and is formulated in a stochastic form. By applying the method of small parameter nonlinear stochastic problem reduces to a system of three linear partial differential equations, which is solved about fluctuations of the stress tensor. This system with transition to the stress function has been reduced to a single differential equation solution of which is represented as a sum of two series. The first row gives the solution away from the boundary of the body without boundary effects, the second row represents the solution boundary layer, its members quickly fade as the distance increases from plate boundary. Based on this solution, the statistical analysis random stress fields near the boundary of the plate was taken.
Key words: steady creep, stochastically inhomogeneous plate, random field of stresses, perturbation theory, boundary effect.
Original article submitted 16/XII/2010; revision submitted 14/II/2011.
Ludmila V. Kovalenko (Ph. D. (Phys. & Math.)), Assistant, Dept. of Applied Mathematics & Computer Science. Nikolay N. Popov (Ph. D. (Phys. & Math.)), Associate Professor, Dept. of Applied Mathematics & Computer Science.