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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2013. № 4 (33). С. 76—85

УДК 539.376

ПРИМЕНЕНИЕ МЕТОДА ВОЗМУЩЕНИЙ К ЗАДАЧЕ О НЕСООСНОЙ ТРУБЕ В УСЛОВИЯХ УСТАНОВИВШЕЙСЯ ПОЛЗУЧЕСТИ

А. Д. Москалик

Самарский государственный технический университет,

Россия, 443100, Самара, ул. Молодогвардейская, 244.

E-mail: [email protected]

Рассмотрена задача определения напряжённо-деформированного состояния толстостенной несоосной трубы, находящейся под внутренним давлением на стадии установившейся ползучести. Проведена линеаризация задачи с использованием метода возмущений. Построено второе приближение данной задачи. Проанализировано влияние несоосности трубы на напряжённо-деформированное состояние с учётом второго приближения.

Ключевые слова: установившаяся ползучесть, метод возмущений, несоосная толстостенная труба, второе приближение.

Введение. Разработка аналитических методов решения краевых задач ползучести для элементов конструкций с возмущенными границами вследствие физической нелинейности определяющих реологических соотношений представляет собой трудноразрешимую проблему. Один из подходов состоит в линеаризации граничных условий и реологических соотношений на основе метода малого параметра. Постановка задачи установившейся ползучести с возмущенными границами методом малого параметра приведена в монографии Л. М. Качанова [1], где, в частности, для несоосной трубы построено решение в первом приближении. Имеются единичные попытки решения этой задачи в работах [2,3]. Хорошо разработаны методы решения стохастических краевых задач для соосной трубы, но с возмущённым по пространственным переменным полем реологических характеристик. Здесь в работах В. П. Радченко и Н. Н. Попова с соавторами [4-6] построены решения вплоть до третьего приближения.

Детально метод возмущений (малого параметра) для упругопластических тел изложен в монографии [7] и систематически развивался в научной школе Д. Д. Ивлева в работах его учеников [8-10 и др.] на случай различных условий пластичности, составных упругопластических тел, различных типов концентраторов и т. д.

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

1. Постановка задачи. Рассматривается несоосная толстостенная труба с несмещенным внутренним контуром радиуса r = а и смещенным на малую величину 5 центром внешнего контура радиуса r = b относительно центра внутренней окружности (рис. 1).

Анна Давидовна Москалик, аспирант, каф. прикладной математики и информатики.

76

Применение метода возмущений к задаче о несоосной трубе...

В качестве малого параметра 5 принимается расстояние между центрами внешнего и внутреннего контуров трубы. Предполагается, что упругие деформации малы по сравнению с деформациями ползучести и ими можно пренебречь. С физической точки зрения это означает, что рассматриваются установившиеся поля скоростей деформаций ползучести и напряжений, т. е. деформацией ползучести, накопленной на первой стадии и вызванной перераспределением напряжений от упругого состояния до состояния установившейся ползучести, пренебрегаем. Разложение тензора напряжений jj, тензора скоростей деформаций ползучести iij и вектора скоростей перемещений Ui по малому параметру до членов второго порядка имеет вид

Рис. 1. Схема несоосной трубы: 1 — внутренний контур трубы r = а; 2 — внешний контур трубы; 3 — внешний контур трубы r = b для осесимметричного случая

jij

£ij Ui -■

CTij0)+5CTij1}+5H2) +O(53)

+54j + 52 ^if + °(53)

u0 + 5u(1) + 52u(2) + O(53).

Уравнение внешнего контура трубы с учетом возмущения 5 имеет вид

(r cos 0 — 5)2 + r2 sin2 0 = b2.

Раскладывая последнее соотношение в степенной ряд по параметру 5 и оставляя члены второго порядка включительно, получаем

r = b + 5 cos 0 + 52(cos 20 — 1)/4b.

Задача решается в условиях плоского деформированного состояния (ezz = 0). Предполагается несжимаемость материала для скоростей деформаций ползучести, что находит экспериментальное подтверждение [11,12]:

£rr + £вв = °. (1)

Постановка задачи включает в себя уравнения равновесия

d(jrr

dr

1 djre jrr — jee djee

r дв r ’ дв

dJre r dr

2jrb i

(2)

которые линейны относительно компонент напряжений и, следовательно, выполняются для каждого приближения.

77

А. Д. Москалик

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

диг , 1 дйо ur 1 (1 dur див ив \

£rr = ——, £вв = + _, £гв = ^ ^ ~~слТГ +—сГ.-----Г" )• (3)

dr

r дв r

2 V r дв дг

r

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

= 3 Aan-1S ■

£ij = о AOe Sij ,

(4)

где n, A — постоянные характеристики материала, Sij = Oj — 3&kk — деви-

нтенс ажен

^3 Г/

•>ij — ^ij 3'

атор напряжений, ое — интенсивность напряжений, определяемая в случае плоской деформации выражением

Ое

2 L

/ \ 2 /. \211/2

\arr — авв) + \4агв) •

Разложение ое по малому параметру S позволяет определить аП 1, представленное в виде степенного ряда:

а

n-1

^У3|Ла(°)| ^

n1

1 + S

(n — 1)Ло(1) |Ло(0)|

+

+ £2 n - 1 ((n — 2) [Ла(1)] 2+2Ла(2)Ла(0) + 4 [Д)]2 (5)

[Ло(0)]

2

где для удобства записи введены обозначения Ло(к) = оД — o^, k = 0,1, 2 номера приближений.

Согласно [12] решение для нулевого приближения имеет вид Ла(0) = —Qp(b/r)p, Д = 0,

при этом ввиду симметричности задачи для нулевого приближения

q

Q=

(b/a)p — 1’

p = 2/n.

Для первого приближения

Ла(1) = Лр(1) (r) cos в, = рГ1 (r) sin в,

где функции Лр(1)(г), p^J(r) согласно [13] имеют следующий вид:

(6)

Лр(1)(г) = — — C12psr 1 + C13pvr w 1 + C14pwr - 1

рГв) (r)

2Г

L '

1Г

2L

C12s2 r-1 + C13v2r-w-1 + CMw2r-v-1

78

Применение метода возмущений к задаче о несоосной трубе...

Здесь

/ /3 )n-1

L = 3Л( — &Q) ,

\ n )

s + д/s2 + 16 p s — л/ s2 + 16 p

s = p — 2, v =---------------------, w =--------------------.

Значения коэффициентов C12, C13, C14 в настоящей работе не приводятся в силу их громоздкости.

Используя (5), (6), линеаризованные определяющие соотношения (4) для второго приближения можно записать в следующем виде:

дст(2) = J^(2) _ a (7)

ДK(r) t'rr a, [ '

K (r) = Lrs, a = AR cos 2 9 + Лф,

AR

Лф

в =

1

2 Да(0) 1

2 Да(0)

n — 1 2

n — 1 2

[Др(1)]2+*[+]2) , [ДР(1)]2 — s[p‘i)]2),

a(2)

arO

BR sin 29

1 ,(2)

K(r) rO

в,

BR

n — 1

2Д0М

(8)

Линеаризованные граничные условия для второго приближения при r = b согласно [7] зависят от нулевого и первого приближений:

a(2)| =

a rO lr=b 2

a,

(2) |

1

rr 'r=b 2 1 dr

Др(1) d_

b dr

dprr 1 d2a>

— 2

PrO -

(0) rr

Да(0)

1 da^T

1

+ 2L

dpr1 dr

2

2b dr

1 d2arr

2 dr2

+

sin 29,

Да(0) b2

da„

1

2b dr

+ 2рГ0) cos 29+ Да(0) 1 (1)'

-----------рГо)

+

2b2

b

(9)

Краевые условия на внутреннем контуре во втором приближении (при r имеют вид

a

(2)

rO

r=a

о,

ar(2r)

r=a

о.

= а) (10)

Уравнения (1)-(3), (7), (8) с граничными условиями (9), (10) образуют краевую задачу для нахождения 2-го приближения в напряжениях и дальнейшего полного построения решения поставленной краевой задачи с учётом второго приближения.

2. Решение задачи. С учётом вида граничных условий (9) и решения для первого приближения [3,13] вводится предположение, что скорость радиаль-

(2)

ных перемещений ur является суммой двух составляющих, одна из которых

зависит от радиуса r и угла 9, а вторая — только от радиуса r; скорость тан-

(2)

генциальных перемещений щ зависит и от радиуса r, и от угла 9:

U-2) = UR cos 29 + Uf, U O2) = UR sin 29,

(11)

79

А. Д. Москалик

где -R = -R(r), -ф = -ф (r), -R = -R(r) — неизвестные, подлежащие определению функции.

Представление для скоростей перемещений (11) позволяет выполнить условие несжимаемости материала (1) тождественно. Для этого необходимо потребовать выполнения следующих равенств:

д- R

дг

2 R

+ _ -r + r

- R r

д-ф -ф

о, -Д + — = о.

dr r

(12)

Первое уравнение из (12) тождественно выполняется путём введения функции скоростей перемещений Z(r, в) = R(r) sin 20 такой, что

-R cos 2 в

1 dZ

r дв,

-R sin 2в

д(

дr ’

(13)

Второе уравнение из (12) позволяет определить составляющую скоростей перемещений, не зависящую от угла в:

- ф = C/г

и, следовательно,

efr = —C/r2. (14)

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

,(2) = _ д(2) = 1 д^_ 1 _0Д + d-f -(2) = 1 /_ 1 Д + Д _ 1 д( \

rr вв r2 дв r двдr dr , гв A r2 дв2 дr2 rOr)

Так как Z(r, в) = R(r) sin2в, выражения для скоростей деформаций ползучести примут вид

егГ2) = — efg = {—2B!r 1 + 2Rr 2) cos 2в +

d-

ф

dr

Д = 1 (R" — R'r-1 + 4Rr-2) sin2в

(15)

Полученные формулы для скоростей деформаций (15) используются в уравнениях (7), (8):

Да(2) = [ 4Р (—R'r~p+1 + Rr~p) — Ar . L

cos 2в

Д = туг \R!'r~p+2 — R'r~p+1 + 4Rr~p 2L -

-^pCr~p — Аф = ДaR cos20 + Даф, (16)

L

sin2в = aR sin2в.

При составлении основного уравнения для нахождения функции R(r) необходимо предварительно продифференцировать по в первое уравнение (2), по r второе уравнение (2), а также найти вторую производную по в и по r

80

Применение метода возмущений к задаче о несоосной трубе...

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

- 2nR

dr2

+ 3-

d-R , д ДaR

dr

-Д- aR

re +2 —— + 4 +2

ДaR

dr

r

0

r

(17)

Опираясь в уравнении (17) на обозначения, введённые в (16), получаем основное уравнение для нахождения функции R(r):

RIV + (6 - 2p)R'"r-1 + (p2 - 20p + 11)R"r-2+

+ (15p2 — 20p + 5)RV-3 + (16 — 12p2 + 8p)Rr-4 = 2Lrp-3Y (r), (18)

где

d2BR dBR dAR BR Ar

Y (r) — r —-t~2—+ 3—-+ 2—-----+ 4---+ 2 .

dr2 dr dr r r

Использование степенного представления R(r) — rv позволяет получить линейное однородное дифференциальное уравнение четвёртого порядка для нахождения собственных значений

v4 — 2pv3 + (p2 — 14p + 4)v2 + (14p2 — 4p)v + (16 — 12p2 + 8p) — 0,

решением которого являются корни

p Jp2 + 28p — 8 ± 4^61^ _ 36p — 12

v1’2’3’4 = 2 ±---------------2----------------.

Прямым численным анализом определено, что при показателе нелинейности n — 3 ^ 20 однородное уравнение имеет комплексные корни:

vi,2,3,4 — p/2 ± 1/2 ± ih.

Здесь l — l(p) и h — h(p) —известные значения для конкретного материала. Воспользуемся тригонометрическим представлением полученного решения согласно [14]:

R(r) — C21r(p+l)/2 cos(h ln r) + C22r(p+l)/2 sin(h ln r)+

+ C23r(p-l)/2 cos(h ln r) + C24r(p-l)/2 sin(h ln r),

где C21, C22, C23, C24 — константы интегрирования.

Использование метода вариации произвольных постоянных [14] позволяет получить общее решение неоднородного дифференциального уравнения в виде

4

R(r) — ^[C2fc (r) + pk ]mfc (r). (19)

k=1

81

А. Д. Москалик

Здесь ш1 = r(p+1)/2 cos(hlnr),m2 = r(p+1)/2 sin(hlnr),m3 = r(p 1)/2 cos(hlnr), ш4 = r(p-l)/2 sin(h ln r); pk — константы интегрирования.

Подставляя полученное решение для функции R(r) (19) и для скорости деформаций ползучести Д (14) в уравнения равновесия (2) и обозначая

daR 1 BR Ar

= 1 (-R"r-p+1 + (4p + 1)R'r-p - (4p + 4)Rr-p-1) - 2 — - —, dr L ' r r

da% 2p i A^

~A~ = ~fCr p + —, dr L r

получаем

даГг

daR

rr- cos 2d +

dafr

dr dr dr

Интегрируя по r, приходим к выражению для радиальной составляющей тензора напряжений

afr = (aRr + KR) cos 2d + (aff + K^),

где KR и K^ — константы интегрирования.

Составляющая тензора напряжений ст^ определяется по второму соотношению (16).

Подставляя полученное решение для а$ и ст^ в граничные условия (9), (10), определяем константы интегрирования p1, p2, p3, p4, C, KR, K^. По-

(2)

скольку уравнение для нахождения arr распадается на два уравнения по признаку наличия зависимости от угла d, получаем шесть уравнений и семь неизвестных констант интегрирования, что позволяет без нарушения общности решения положить KR = 0.

Использование полученного решения позволяет определить напряжения и скорости деформаций ползучести в трубе с учётом второго приближения.

Анализ результатов. В качестве модельного примера рассмотрена труба с внутренним радиусом a = 0,115 м, внешним радиусом b = 0,15 м и с характеристиками материала n = 10,96, A = 9,58 ■ 10-23 под действием внутреннего давления q = 22,07 МПа.

В работе детально исследованы зависимости arr (r,d), а$в (r,d) и err (r,d) от приведённого радиуса r = r/a.

На рис. 2-4 соответственно приведены графики для радиальной и тангенциальной компонент тензора напряжений и радиальной компоненты тензора скоростей деформаций (поскольку £$$ = — err), где верхний индекс у компонент тензора напряжений и тензора скоростей деформаций означает количество приближений, используемых для оценки величин. Графики построены при 5 = 0,1 при значении угла d = п, соответствующем максимальным значениям тангенциального напряжения . В таблице приведены значения ®вв = а(0+1)/0дв при учёте первого приближения на внешней границе трубы

при r = b/a + 5 cos d и значения аЦ = а(0+1+2)/a(e°J при учёте второго приближения на внешней границе трубы при r = b/a + 5 cos d + 52 (cos 2d — 1)/4b

82

Применение метода возмущений к задаче о несоосной трубе...

Рис. 2. Радиальные напряжения при n = 10,96,0 = п,4 = 0,1: 1 —

„(0) 2 — „(0+1) 3 — „(0+1+2)

Рис. 3. Тангенциальные напряжения при n = 10,96,0 = п,4

1

40,2 — „<°+1),3-

(0+1+2)

ее

0,1:

Рис. 4.

Скорости радиальных деформаций при n 4 = 0,1: 1— 40, 2-

t(0+1) t-rr ?

3—

• (0+1+2) Ьгг

10,96, 0

п,

Значения тангенциальных напряжений

4 0 0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,10

„ее 1,0 1,04 1,09 1,13 1,17 1,22 1,26 1,30 1,34 1,38 1,42

„ее 1,0 1,04 1,10 1,14 1,19 1,24 1,29 1,34 1,40 1,46 1,52

83

А. Д. Москалик

при n = 10,96, в = п, вычисленные с шагом 0,01 по величине 5. Из данных, приведённых в таблице, можно сделать вывод, что решение задачи о несоосной трубе имеет тенденцию к сходимости. Однако для изучения скорости сходимости решения нужны дополнительные исследования.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Л. М. Качанов, Теория ползучести. М.: Физматгиз, 1960. 455 с. [L. M. Kachanov, Creep Theory. Moscow: Fizmatgiz, 1960. 455 pp.]

2. В. П. Радченко, Е. В. Башкинова, “Решение краевых задач установившейся ползучести в полярных координатах методом возмущений” // Вестн. Сам. гос. техн. ун-та. Сер. Технические науки, 1998. №5. С. 86-91. [V. P. Radchenko, E. V. Bashkinova, “Solution of the value boundary problems for steady creep in polar coordinates by the perturbation method” // Vestn. Samar. Gos. Tekhn. Univ. Ser. Tekhnicheskie Nauki, 1998. no. 5. Pp. 8691].

3. Е. В. Башкинова, “Решение краевой задачи установившейся ползучести для неосесимметричной толстостенной трубы”// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2002. № 16. С. 105-110. [E. V. Bashkinova, “Solution of the value boundary problem of steady creep for non-axisymmetric thick-walled tube” // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2002. no. 16. Pp. 105-110].

4. В. П. Радченко, Н. Н. Попов, “Аналитическое решение стохастической краевой задачи установившейся ползучести для толстостенной трубы”// ПММ, 2012. Т. 76, №6. С. 1023-1031; англ. пер.: N. N. Popov, V. P. Radchenko, “Analytical solution of the stochastic steady-state creep boundary value problem for a thick-walled tube” // J. Appl. Math. Mech., 2012. Vol.76, no. 6. Pp. 738-744.

5. А. А. Должковой, Н. Н. Попов, “Решение нелинейной стохастической задачи ползучести для толстостенной трубы методом малого параметра” // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2002. №16. С. 84-89. [A. A. Dolzhkovoy, N. N. Popov, “Solution of the nonlinear stochastic creep problem for a thick-walled tube by method of small parameter” // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2002. no. 16. Pp. 84-89].

6. Н. Н. Попов, В. Н. Исуткина, “Построение аналитического решения двумерной стохастической задачи установившейся ползучести для толстостенной трубы” // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2007. №2(15). С. 90-94. [N. N. Popov, V. N. Isutkina, “Construction of an Analytical Solution of a Two-Dimensional Stochastic Problem of the Steady Creep for a Thick-Walled Pipe” // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2007. no. 2(15). Pp. 90-94].

7. Д. Д. Ивлев, Л. В. Ершов, Метод возмущений в теории упругопластического тела. М.: Наука, 1978. 208 с. [D. D. Ivlev, L. V. Ershov, Method of Perturbations in the Theory of an Elastic-Plastic Body. Moscow: Nauka, 1978. 208 pp.]

8. А. П. Кержаев, “Упругопластическое состояние тонкой кольцевой пластины при наличии трансляционной анизотропии при равномерном растяжении” // Вестник Чуваш. гос. пед. ун-та им. И. Я. Яковлева. Серия: Механика предельного состояния, 2012. №2(12). С. 174-179. [A. P. Kerzhaev, “Elastoplastic state of a thin ring plate at transmitting anisotropy at uniform stretching” // Vestnik of I. Yakovlev Chuvash State Pedagogical University, Line: Mechanics of definable state, 2012. no. 2(12). Pp. 174-179].

9. С. О. Фоминых, “Упругоидеальнопластическое состояние анизотропной трубы” // Вестник Чуваш. гос. пед. ун-та им. И. Я. Яковлева. Серия: Механика предельного состояния, 2010. №2(8). С. 623-627. [S. O. Fominykh, “Elastoidealplastic condition of anisotropic pipe” // Vestnik of I. Yakovlev Chuvash State Pedagogical University, Line: Mechanics of definable state, 2010. no. 2(8). Pp. 623-627].

10. С. О. Фоминых, “Упругопластическое состояние толстостенной трубы при взаимодействии различных видов пластической анизотропии” // Вестник Чуваш. гос. пед. ун-та им. И. Я. Яковлева. Серия: Механика предельного состояния, 2011. №1(9). С. 201-206.

84

Применение метода возмущений к задаче о несоосной трубе...

[S. O. Fominykh, “Elastoplastic state of a thick-walled pipe at the interaction of different types of plastic anisotropy” // Vestnik of I. Yakovlev Chuvash State Pedagogical University, Line: Mechanics of definable state, 2011. no. 1(9). Pp. 201-206].

11. А. Ф. Никитенко, Ползучесть и длительная прочность металлических материалов. Новосибирск: НГАСУ, 1997. 278 с. [A. F. Nikitenko, Creep and Creep-Rupture Strength of Metallic Materials. Novosibirsk: Novosib. Gos. Arkhit.-Stroit. Univ., 1997. 278 pp.]

12. Ю. Н. Работное, Механика деформируемого твердого тела. М.: Наука, 1979. 744 с. [Yu. N. Rabotnov, Mechanics of Deformable Solids. Moscow: Nauka, 1979. 744 pp.]

13. А. Д. Москалик, “Анализ напряжённо-деформированного состояния толстостенного несоосного цилиндра, находящегося под внутреннем давлением, в условиях установившейся ползучести методом малого параметра” / В сб.: Труды девятой Всероссийской научной конференции с международным участием.Часть 1: Математические модели механики, прочности и надёжности элементов конструкций / Матем. моделирование и краев. задачи. Самара: СамГТУ, 2013. С. 140-144. [A. D. Moskalik, “Analysis of the stress-strain state of a thick-walled misaligned cylinder in conditions of steady-state creep under internal pressure by method of small parameter” / In: Proceedings of the Ninth AllRussian Scientific Conference with international participation. Part 1 / Matem. Mod. Kraev. Zadachi. Samara: Samara State Technical Univ., 2013. Pp. 140-144].

14. Л. Э. Эльсгольц, Дифференциальные уравнения и вариационное исчисление. М.: Наука, 1969. 424 с. [L. E. Elsgolts, Differential Equations and the Calculus of Variations. Moscow: Nauka, 1969. 424 pp.]

Поступила в редакцию 14/IX/2013; в окончательном варианте — 26/XI/2013.

MSC: 74D10, 74G10

THE APPLICATION OF PERTURBATION METHOD TO PROBLEM OF MISALIGNED TUBE IN CONDITIONS OF STEADY-STATE CREEP

A. D. Moskalik

Samara State Technical University,

244, Molodogvardeyskaya st., Samara, 443100, Russia.

E-mail: [email protected]

The problem of determining the stress-strain state of the thick-walled misaligned tube under internal pressure on steady-state creep is considered. The task linearization with the perturbation method is carried out. The second approximation of this problem is constructed. The effect of misalignment of the tube on the stress-strain state considering the second approximation is analyzed.

Keywords: steady-state creep, pertubation method, thick-walled misaligned tube, second approximation.

Original article submitted 14/IX/2013; revision submitted 26/XI/2013.

Anna D. Moskalik, Postgraduate Student, Dept. of Applied Mathematics & Computer Science.

85