CC BY
47
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2014. № 3 (36). С. 79—93
ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1340
УДК 539.376
СРАВНИТЕЛЬНЫЙ АНАЛИЗ ПРИБЛИЖЕННОГО АНАЛИТИЧЕСКОГО И КОНЕЧНО-ЭЛЕМЕНТНОГО РЕШЕНИЙ ДЛЯ НЕСООСНОЙ ТРУБЫ
В. П. Радченко, А. Д. Москалик, И. Е. Адеянов
Самарский государственный технический университет,
Россия, 443100, Самара, ул. Молодогвардейская, 244.
Аннотация
Рассмотрена краевая задача установившейся ползучести для несоосной толстостенной трубы, находящейся под внутренним давлением. Строится приближенное аналитическое решение данной задачи методом малого параметра до второго приближения включительно. Решение строится для плоского деформированного состояния. Используется гипотеза несжимаемости материала для деформаций ползучести. В качестве малого параметра используется величина смещения центров внутреннего и внешнего радиусов трубы. Основное внимание уделено вопросу сходимости полученного аналитического решения во втором приближении и оценке его погрешности. Отмечается, что проблема сходимости решена только для краевых задач в упругой области. Поэтому оценка погрешности в поставленной задаче решалась на основании сравнения приближенного аналитического решения с численным решением, построенным методом конечных элементов, для некоторых частных случаев. С учетом симметрии задачи конечно-элементная модель была построена для половины трубы. Количество конечных элементов — около 18 000. Вторая половина трубы заменялась граничными условиями с учетом симметрии задачи. Анализ аналитического и численного решений выполнен в зависимости от параметра нелинейности установившейся ползучести и параметра несоосности — отношения смещения центров внешнего и внутреннего диаметров к внешнему радиусу. Показано, что для труб с малыми значениями показателя установившейся ползучести (от 3 до 8) погрешность отклонения приближенного аналитического решения во втором приближении от численного решения вплоть до величины несоосности центров внутреннего и внешнего диаметров 0.1 составляет не более 9 %, а для труб с большим показателем нелинейности установившейся ползучести погрешность до 8 % наблюдается при параметре несоосности до 0.06. Приводятся результаты расчетов в табличной форме и в форме графиков. Даны рекомендации по использованию построенного приближенного аналитического решения в прикладных задачах.
© 2014 Самарский государственный технический университет.
Образец для цитирования: Радченко В. П., Москалик А. Д., Адеянов И. Е. Сравнительный анализ приближенного аналитического и конечно-элементного решений для несоосной трубы // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014. № 3(36). С. 7993. doi: 10.14498/vsgtu1340.
Сведения об авторах: Владимир Павлович Радченко (д.ф.-м.н., проф.; radch@samgtu.ru), заведующий кафедрой, каф. прикладной математики и информатики. Анна Давидовна Москалик (annmoskalik1@gmail.com; автор, ведущий переписку), аспирант, каф. прикладной математики и информатики. Игорь Евгеньевич Адеянов (к.т.н.; adigorev@gmail.com), доцент, каф. механики.
79
Радченко В. П., Москалик А. Д., Адеянов И. Е.
Ключевые слова: несоосная толстостенная труба, установившаяся ползучесть, приближенное аналитическое решение, метод малого параметра, второе приближение, конечно-элементная модель, численное решение, погрешность.
doi: http://dx.doi.org/10.14498/vsgtu1340
Введение. Разработка аналитических методов решения краевых задач ползучести для элементов конструкций с возмущенными границами вследствие физической нелинейности определяющих реологических соотношений представляет трудноразрешимую проблему. Один из подходов состоит в линеаризации граничных условий и реологических соотношений на основе метода малого параметра. Число работ в данном направлении необозримо. Однако, что касается внешних краевых задач реологии (с возмущенными границами), то здесь следует отметить ряд работ по устойчивости однородного растяжения полосы или цилиндра из вязкого материала, чувствительного к скорости деформирования, по отношению к малым возмущениям регулярной или произвольной формы свободных границ, решенных методом малого параметра [1-4]. С другой стороны, развиваются методы решения краевых задач с возмущенным по пространственным переменным полем реологических характеристик (внешние краевые задачи). Так, в работах [5-8] методом малого параметра построены аналитические решения для полей напряжений и скоростей деформаций вплоть до третьего приближения в стохастической краевой задаче установившейся ползучести толстостенной трубы под действием внутреннего давления, а в [9, 10] приведены решения аналогичной задачи ползучести для растягиваемой плоскости.
Постановка задачи установившейся ползучести с возмущенными границами методом малого параметра приведена в монографии Л. М. Качанова [11], где, в частности, для несоосной трубы построено решение в первом приближении. Имеются попытки решения этой задачи в работах [12,13].
Детально метод возмущений (малого параметра) для упругопластических тел изложен в монографии [14] и систематически развивался в научной школе Д. Д. Ивлева в работах его учеников [15-18 и др.] для различных условий пластичности, составных упругопластических тел, различных типов концентраторов и т. д.
К сожалению, практически во всех работах по решению внешних и внутренних краевых задач методом малого параметра для сред с реономными свойствами отсутствуют оценки сходимости рядов, поскольку это является сложной и неформализованной задачей. Здесь можно отметить только работу [19], в которой этот вопрос решен для микронеоднородной упругой среды.
В связи с вышеизложенным целью данной работы является сравнение приближенного аналитического решения задачи об установившейся ползучести толстостенной несоосной трубы, находящейся под внутренним давлением, методом малого параметра с учетом приближений до второго включительно и конечно-элементного решения данной задачи с помощью программного комплекса инженерного анализа ANSYS.
1. Постановка задачи. Рассматривается несоосная толстостенная труба с внутренним контуром радиуса r = a и смещенным на малую величину 5 центром внешнего контура радиуса r = b относительно центра внутренней окружности под действием внутреннего давления q (рис. 1).
80
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
Рис. 1. Схема несоосной трубы: 1 — внутренний контур трубы r = a; 2 — внешний смещенный контур трубы; 3 — внешний контур трубы r = b для осесимметричного случая
[Figure 1. The scheme of a misaligned tube: 1 — the inner contour of the tube; 2 — the outer displaced contour of the tube; 3 — the outer contour of the tube r = b for axisymmetric case;
5 is the small parameter]
В качестве малого параметра 5 принимается расстояние между центрами внешнего смещенного и внутреннего контуров трубы. При получении приближенного аналитического решения предполагается, что упругие деформации малы по сравнению с деформациями ползучести и ими можно пренебречь. С физической точки зрения это означает, что рассматриваются установившиеся поля скоростей деформаций ползучести и напряжений, т. е. деформацией ползучести, накопленной на первой стадии и вызванной перераспределением напряжений от упругого состояния до состояния установившейся ползучести, пренебрегаем. Разложение тензора напряжений aij, тензора скоростей деформаций ползучести еj и вектора скоростей перемещений Щ по малому параметру до членов второго порядка имеет вид
aij £ij Щ :
(0)
ij
+5ст(1) + 52 + O(53)
e°j + 5eij + 5С^ + O(5T
u0 + 5u(1) + 52щ(2) + O(53),
где индексы 0, 1 и 2 соответствуют нулевому, первому и второму приближениям.
Уравнение внешнего контура трубы с учетом возмущения 5 имеет вид (r cos в — 5)2 + r2 sin2 в = Ь2.
Раскладывая последнее соотношение в степенной ряд по параметру 5 и ограничиваясь членами второго порядка включительно, получаем
r = Ь + 5 cos в + 52(cos 2в — 1)/4b.
Задача решается в условиях плоского деформированного состояния:
ezz
0.
81
Радченко В. П., Москалик А. Д., Адеянов И. Е.
Предполагается несжимаемость материала для скоростей деформаций ползучести, что находит экспериментальное подтверждение [20,21]:
ётт + £$$ = 0.
Постановка задачи включает в себя уравнения равновесия
do'rr 1 d@rQ a rr aQQ
dr r 39 r ’
3UQQ 3UrQ n
d9 = Г dr rQ ,
которые линейны относительно компонент напряжений и, следовательно, выполняются для каждого приближения.
Аналогично, для каждого приближения выполняются уравнения совместности деформаций
£rr —
£tQ
dur dr ’
1 3uq ur
r d9 + r ’
1 /1 dur duQ
2 V r d9 + dr
r
В качестве определяющих используются соотношения теории установившейся ползучести со степенным законом
£ij
^ AarlSij,
где n, A — постоянные характеристики материала, Sij — aij — a^k/3 —девиа-тор напряжений, ae — интенсивность напряжений для случая плоской деформации.
2. Приближенное аналитическое решение задачи. Согласно [21] решение для нулевого приближения имеет вид
Дт (r) — Q[1 — (b/r)p],
aQ°>(r) = Q[1 — (1 — P)(b/r)p], (1)
ai°z>(r) — Q[ 1 — (1 + p)(b/r)p],
где
Q = Щгг1p = 2/n-
При этом стТ°> = 0 ввиду симметричности задачи для нулевого приближения.
Поскольку граница при r = а не возмущена и задано давление q, линеаризованное граничное условие на внутреннем радиусе трубы для последующих (после нулевого) приближений представимо в виде
a
(k>
rr
r=a
о,
a(k>
arQ
r=a
0,
82
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
где k = 1, 2 — номера приближений.
Линеаризованное граничное условие для первого приближения при r = b согласно [14] зависит от нулевого приближения:
a
(1)
rr
r=b
daffl
dr
cos 9,
arl)
r=b
Aa(0)
sin 9 ~b~
где Aa(0) = arr) — a^.
Для первого приближения согласно [22] выполняется
arr) (r, 9) = рГГ) (r) cos 9, a$ (r, 9)
p$(r)cos 9, ar^)
рГС (r) sin 9.
Функции рГГ) (r), P^0) (r), рГ^)(г) имеют следующий вид:
рГГ)(г) = 71 Cii+Ci2s(s—4p)r 1+Ci3
2L L W + 1
Р«(г) = ^[c„ + C12s2r-1 + Ci3 'w + 1
„-1,n v(v — 4P^-w-1
++C14
v(v + 4pw) -w-1
r-
+ +C14
w(w — 4P) r-v-1 v + 1 .
w(w + 4pv) -v-1" v + 1 -
1
Pr^r) = 2L [C12S2r 1 + Ci3v2r w 1 + C^wV
где
L = 3A ^bPQ)“-1, s = p — 2, v = s + ^ + 16p, w = s — ^ + 16p.
Здесь константы интегрирования C1i, i = 1, 2, 3, 4, определяются из граничных условий для первого приближения.
Линеаризованное граничное условие для второго приближения при r = b согласно [14] зависит от нулевого и первого приближений:
a(2) I -
°rQ lr=b
a
(2) | = rr r=b
1 [ Ap(1)
2 .
1 -
2 .
d
рП
b dr
dprr) 1 d2a(0)
dr
1 [
+ 2 1
Здесь Ap(1) = рГГ^ — /4с ■
dp
(4)
dr2
Aa(0) 11 ——J sin 29,
Aa(0) 2 (1)1
+ ~bP^\ cos29+
1 dar0r)
1 d2a-
dr
dr b2
1 dar0r) Aa(0)
+-------rA +
+ 2b dr +
2 dr2 2b dr 2b2
Для второго приближения согласно [23] выполняется:
1 (1) bPrC
a(r)(r, 9) = a^T(R)cos 29 + a^r(r),
4?^ 9) = aR(R) cos 29 + acC(r),
a^r, 9) = a^ (R) sin 29.
Решение зависит от функции R(r), имеющей вид
4
R(r) = ^[C2fc (r) + Pk ]mfc (r).
k=1
83
Радченко В. П., Москалик А. Д., Адеянов И. Е.
Здесь C2k(г) —известные функции радиуса г, ш1 = г(р+1)/2cos(hlnг),ш2 = = r(p+1)/2 sin(h ln г),шз = r(p-1)/2 cos(h ln г), m4 = r(p-1)/2 sin(h ln г), l = l(p) и h = h(p) —известные значения для конкретного материала; pk — константы интегрирования.
Использование первого и второго приближений метода возмущений позволяет определить напряжения и скорости деформаций ползучести в несоосной трубе. Исследование необходимости построения третьего и более высоких приближений связано с величиной погрешности первых двух приближений, аналитически оценить которую не представляется возможным. Поэтому оценка погрешности приближенного аналитического решения была выполнена на основании его сравнения с численным решением, построенным на основе метода конечных элементов.
3. Анализ приближенного аналитического решения. В качестве модельного примера рассмотрена труба с внутренним радиусом a = 0.115 м, внешним радиусом b = 0.15 м под действием внутреннего давления q = 22.07 МПа. В качестве материалов рассмотрены углеродистая сталь [24] и жаропрочный сплав ХН73МБТЮ (ЭИ698) [25], реологические характеристики которых следующие:
углеродистая сталь: n = 3.03, A = 9.04 ■ 10-9;
ХН73МБТЮ (ЭИ698): n = 10.96, A = 4.57 ■ 10-33.
В табл. 1 приведены безразмерные значения тангенциального напряжения на внешней границе трубы в первом ale(r, в) = a(0+1)(f, в)/а<(0в(г, в) при r = b/a + 5 cos в, в = п и во втором аЦ(¥, в) = a^+1+2)(f, в)/а^(т, в) при r = b/a + 5 cos в + 52(cos2в — 1)/4b, в = п приближениях в зависимости от безразмерной величины малого параметра 5 = b/a с шагом 0.01 для углеродистой стали и жаропрочного сплава ХН73МБТЮ (ЭИ698).
Из данных, приведенных в табл. 1, можно сделать вывод, что решение задачи о несоосной трубе имеет тенденцию к сходимости. Однако для теоретического вывода о скорости сходимости решения, построенного по методу малого параметра, данных о первых двух приближениях, вообще говоря, недостаточно. Поэтому следующим шагом исследования является численное решение поставленной задачи и сравнение данных приближенного аналитического и численного решений.
Таблица 1
Значения тангенциального напряжения для несоосной трубы [Values of tangential stresses for misaligned tube]
s 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Углеродистая сталь [Carbon Steel]
aee 1.0 1.04 1.09 1.13 1.17 1.21 1.25 1.29 1.33 1.37 1.40
aee 1.0 1.04 1.09 1.14 1.19 1.24 1.29 1.35 1.40 1.46 1.53
ХН73МБТЮ (ЭИ698) [KHN73MBTYU (EI698) Alloy]
aee 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
84
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
4. Конечно-элементная дискретная модель. Для построения численного решения разработана конечно-элементная модель толстостенной несоосной трубы, находящейся под внутренним давлением, с помощью программного комплекса ANSYS.
Решение выполняется двумя шагами:
1) упругое решение;
2) решение с учётом свойств ползучести материала за время 1 ■ 103 ч, поскольку за это время напряженное состояние выходит на стационарный режим, соответствующий стадии установившейся ползучести.
Модуль упругости E и плотность материала р принимаются постоянными для выбранного материала и температуры: углеродистая сталь: E = 1.56 ■ 105 МПа, р = 7630 кг/м3, T = 649°С;
ХН73МБТЮ (ЭИ698): E = 1.44 ■ 105 МПа, р = 7900 кг/м3, T = 775С.
В расчетах использовался плоский восьмиузловой элемент PLANE183, пригодный для моделирования плоского деформированного состояния и позволяющий учитывать большие деформации при установившейся ползучести. В силу симметрии задачи конечно-элементная модель строилась для половины несоосной трубы. Количество элементов для половины трубы составляет порядка 18 000. Отсутствующая часть трубы заменена условием симметричности по оси 0х (см. рис. 1).
Для оценки конечно-элементной модели на предварительном этапе решалась задача для осесимметричной трубы, находящейся в условиях установившейся ползучести под внутреннем давлением q. Решение задачи в такой постановке соответствует нулевому приближению (1) поставленной задачи, аналитическое решение которой хорошо известно [21].
При значениях времени t = 10 часов изменения напряжений становятся пренебрежимо малы, а деформации ползучести изменяются линейно, что позволяет говорить о состоянии установившейся ползучести материала при t ^ 10 часов. В качестве решения задачи в условиях установившейся ползучести использовались данные, соответствующие времени 103 часов.
5. Сравнение результатов приближенного аналитического и конечно-элементного решений. В работе проведена оценка погрешности приближенного аналитического и конечно-элементного решений задачи для осесимметричной и несоосной трубы с учетом второго приближения включительно на основе значений радиальных orr и тангенциальных напряжений в 15 равноотстоящих точках по координате r+
а ^ n ^ b + 5 cos 6 + 52(cos 26 — 1)/4, i = 1,2,..., 15
при 6 = 0 и 6 = п. Вычисление погрешности проведено по двум нормам:
Y^15
2.^=1
о.
(0+1+2)
- о:
£1=1 1о
и О
ans
\
у^15 Г (0+1+2) оans] 2\ 1/2
/Аг=1 J \
Z^i=1
5 roansl ‘ = 1
ш = r, в,
s
где О(0+1+2) = o:0:+1+2)(ri, 6), oans = °^(п,в) —расчетные значения для аналитического (два приближения) и численного решений соответственно. В табл. 2 приведены оценки различия результатов расчетов между приближенным аналитическим методом и конечно-элементным методом в процентах
85
ос
05
Таблица 2
Оценки различия результатов расчетов между приближенным аналитическим методом и конечно-элементным методом [Evaluation of the differences in the calculation results between the approximate analytical method and the finite element method]
s 0 0.02 0.04 0.06 0.08 0.10
Углеродистая сталь [Carbon Steel]
crrr(r, 0) 0.12/0.34 0.21/0.41 0.22/0.42 0.20/0.41 0.28/0.42 0.47/0.51
0>r(r,7r) 0.12/0.34 0.21/0.41 0.25/0.42 0.27/0.44 0.33/0.51 0,50/0,65
&ee(r, 0) 0.009/0.023 0.32/0.32 1.31/1.32 3.00/3.05 5.41/5.57 8.56/8.96
<?ee{r, л) 0.009/0.023 0.32/0.35 1.24/1.29 2.75/2.80 4.80/4.85 7.37/7.40
ХН73МБТЮ (ЭИ698) [KHN73MBTYU (Ё1698) Alloy]
(Jrr (^7 0) 0.17/0.38 0.18/0.40 0.34/0.44 1.10/1.07 — —
(7rr ((5, 7Г) 0.17/0.38 0.21/0.40 0.40/0.46 0.55/0.56 — —
°ее{5Л) 0.02/0.06 0.36/0.41 1.81/2.02 6.10/7.05 — —
(твв(5,Д 0.02/0.06 0.32/0.35 1.22/1.24 2.61/2.62 — —
Радченко В. П., Москалик А. Д., Адеянов И. E.
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
г = г/а
a
г = г/а
b
Рис. 2. Радиальные напряжения <rrr для несоосных труб из малоуглеродистой стали (а) и сплава ХН73МБТЮ (ЭИ698) (b) при в = 0, А = 0.06: 1- а(г0},2 —а(°r+1),3 —a(°r+1+2), 4 —
[Figure 2. Radial stresses arr for carbon steel (a) and KHN73MBTYU (EI698) Alloy (b) misaligned tubes, when в = 0, <1 = 0.06: 1 — ст(Г), 2 — стГГ+1), 3 — стГГ+1+2), 4 — o-Arns]
87
Радченко В. П., Москалик А. Д., Адеянов И. Е.
г = г/а
a
г = г/а
b
Рис. 3. Тангенциальные напряжения аее для несоосных труб из малоуглеродистой стали (а) и сплава ХН73МБТЮ (ЭИ698) (b) при в = п, <5
0.04: 1- а(Д, 2 — а(0+1),
3 — а
(0+1+2)
4 — аА
[Figure 3. Tangential stresses аее for carbon steel (a) and KHN73MBTYU (EI698) Alloy (b) misaligned tubes, when в = п, < = 0.04: 1 — ад0, 2 — а'д0е+1'), 3 — адд+1+2), 4 — aAens]
88
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
для малоуглеродистой стали и сплава ХН73МБТЮ (ЭИ698), используемых в качестве модельных материалов. Через дробную черту приведены оценки, вычисленные по двум нормам: s/a.
В качестве иллюстрации на рис. 2 приведены графики радиальной компоненты arr тензора напряжений для углеродистой стали (n = 3.03) и жаропрочного сплава ХН73МБТЮ (ЭИ698) (n = 10.96) при в = 0, 5 = 0.06. Верхний индекс у компонент тензора напряжений означает количество приближений, используемых для оценки величин. Индекс (Ans) означает конечно-элементное решение. На рис. 3 приведены графики для тангенциальной компоненты a$$ тензора напряжений при 5 = 0.04 и значении угла в = п, соответствующего максимальным значениям тангенциального напряжения a$$.
Анализ данных таблицы 2 и графиков на рис. 2, 3 позволяет сделать следующие выводы:
1) для труб с малыми значениями показателя установившейся ползучести материала (малоуглеродистая сталь) погрешность отклонения приближенного аналитического решения во втором приближении от численного решения по методу конечных элементов вплоть до величины несоосности центров внутреннего и внешнего диаметров 5 = 0.1 составляет не более 9 %; для труб с большим показателем нелинейности установившейся ползучести материала (сплав ХН73МБТЮ (ЭИ698)) для величины 0 ^ 5 ^ 0.06 — погрешность не более 8%, а для значений 5 > 0.06 — она существенно возрастает;
2) в расчетной практике для труб с небольшим показателем нелинейности установившейся ползучести материала (3 ^ n ^ 8) можно ограничиться вторым приближением вплоть до величины 5 = 0.1 (с погрешностью решения до 9%); для труб с показателем нелинейности установившейся ползучести материала n > 8 с погрешностью не более 9 % приближенное аналитическое решение во втором приближении можно использовать лишь до величины 5 = 0.06, а для величины 5 > 0.06, по всей видимости, необходимо иметь приближения более высокого порядка.
ORCID
Vladimir P. Radchenko: http://orcid.org/0000-0003-4168-9660 Anna Moskalik: http://orcid.org/0000-0001-7527-6237 Igor Adeyanov: http://orcid.org/0000-0002-3943-6173
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Hill R., Hutchinson J. W. Bifurcation phenomena in the plane tension test // J. Mech. Phys. solids, 1975. vol. 23. pp. 239-264. doi: 10.1016/0022-5096(75)90027-7.
2. Storen S., Rice J. R. Localized necking in thin sheets // J. Mech. Phys. solids, 1975. vol. 23, no. 6. pp. 421-441. doi: 10.1016/0022-5096(75)90004-6.
3. Hutchinson J. W., Neale K. W. Influence of strain-rate sensitivity on necking under uniaxial tension// Acta Metallurgica, 1977. vol. 25, no. 8. pp. 839-846. doi: 10.1016/0001-6160(77) 90168-7.
4. Келлер И. Э. Равновесные формы свободной границы при одноосном растяжении нелинейно-вязкой полосы// ПМТФ, 2010. Т. 51, №1. С. 117-124.
5. Радченко В. П., Попов Н. Н. Аналитическое решение стохастической краевой задачи установившейся ползучести для толстостенной трубы // ПММ, 2012. Т. 76, №6. С. 10231031.
89
Радченко В. П., Москалик А. Д., Адеянов И. Е.
6. Должковой А. А., Попов Н. Н. Решение нелинейной стохастической задачи ползучести для толстостенной трубы методом малого параметра // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2002. №16. С. 84-89. doi: 10.14498/vsgtu102.
7. Попов Н. Н., Исуткина В. Н. Построение аналитического решение двумерной стохастической задачи установившейся ползучести для толстостенной трубы // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2007. №2(15). С. 57-61. doi: 10.14498/ vsgtu535.
8. Должковой А. А., Попов Н. Н., Радченко В. П. Решение стохастической краевой задачи установившейся ползучести для толстостенной трубы методом малого параметра // ПМТФ, 2006. Т. 47, №1. С. 161-171.
9. Коваленко Л. В., Попов Н. Н., Радченко В. П. Решение плоской стохастической краевой задачи ползучести// ПММ, 2009. Т. 73, №6. С. 1009-1016.
10. Попов Н. Н., Самарин Ю. П. Исследование полей напряжений вблизи границы стохастически неоднородной полуплоскости при ползучести // ПМТФ, 1988. №1. С. 159-164.
11. Качанов Л. М. Теория ползучести. М.: Физматгиз, 1960. 455 с.
12. Радченко В. П., Башкинова Е. В. Решение краевых задач установившейся ползучести в полярных координатах методом возмущений // Вестн. Сам. гос. техн. ун-та. Сер. Техн. науки, 1998. №5. С. 86-91.
13. Башкинова Е. В. Решение краевой задачи установившейся ползучести для неосесимметричной толстостенной трубы // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2002. №16. С. 105-110. doi: 10.14498/vsgtu106.
14. Ивлев Д. Д., Ершов Л. В. Метод возмущений в теории упругопластического тела. М.: Наука, 1978. 208 с.
15. Кержаев А. П. Упругопластическое состояние тонкой кольцевой пластины при наличии трансляционной анизотропии при равномерном растяжении // Вестник Чувашского государственного педагогического университета им. И. Я. Яковлева. Серия: Механика предельного состояния, 2012. №2(12). С. 174-179.
16. Фоминых С. О. Упругопластическое состояние толстостенной трубы при взаимодействии различных видов пластической анизотропии // Вестник Чувашского государственного педагогического университета им. И. Я. Яковлева. Серия: Механика предельного состояния, 2011. №1(9). С. 201-2016.
17. Петров Н. И. О деформировании растягиваемой полосы, ослабленной пологими выточками // Вестник Чувашского государственного педагогического университета им. И. Я. Яковлева. Серия: Механика предельного состояния, 2014. №2(20). С. 36-45.
18. Никитин А. В., Тихонов С. В. Предельное состояние многослойной трансляционноанизотропной толстостенной трубы, находящейся под действием внутреннего давления // Вестник Чувашского государственного педагогического университета им. И. Я. Яковлева. Серия: Механика предельного состояния, 2014. №1(19). С. 88-94.
19. Кунташев П. А., Немировский Ю. В. О сходимости метода возмущений в задачах теории упругости// Изв. Акад. наук СССР. Мех. тверд. тела, 1985. №3. С. 75-78.
20. Никитенко А. Ф. Ползучесть и длительная прочность металлических материалов. Новосибирск: НГАСУ, 1997. 278 с.
21. Работнов Ю. Н. Механика деформируемого твердого тела. М.: Наука, 1979. 744 с.
22. Москалик А. Д. Анализ напряженно-деформированного состояния толстостенного несоосного цилиндра, находящегося под внутреннем давлением, в условиях установившейся ползучести методом малого параметра / Труды Девятой Всероссийской научной конференции с международным участием/ Матем. моделирование и краев. задачи. Самара: СамГТУ, 2013. С. 140-144.
23. Москалик А. Д. Применение метода возмущений к задаче о несоосной трубе в условиях установившейся ползучести// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2013. №4(33). С. 76-85. doi: 10.14498/vsgtu1290.
24. Малинин Н. Н. Прикладная теория пластичности и ползучести. М.: Машиностроение, 1975. 400 с.
90
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
25. Радченко В. П., Саушкин М. Н. Ползучесть и релаксация остаточных напряжений в упрочненных конструкциях. М.: Машиностроение-1, 2005. 226 с.
Поступила в редакцию 26/V/2014; в окончательном варианте — 18/VI/2014; принята в печать — 27/VIII/2014.
Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki. 2014. Issue 3 (36). Pp. 79—93 [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci. 2014. Issue 3 (36). Pp. 79—93]
ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1340
MSC: 74D10, 74G10
COMPARATIVE ANALYSIS OF THE APPROXIMATE ANALYTICAL AND FINITE ELEMENT SOLUTIONS FOR MISALIGNED TUBE
V. P. Radchenko, A. D. Moskalik, I. E. Adeyanov
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.
Abstract
The boundary value problem of steady-state creep for thick-walled misaligned tube under internal pressure was considered. The approximate analytical solution of this problem by method of small parameter including the second approach is under construction. The solution for the state of plane deformation is constructed. The hypothesis of incompressibility of material for creep strain is used. As a small parameter the misalignment of the centers of the inner and outer radii of the tube is used. The main attention to the convergence of the resulting analytical solution considering the second approximation and assessment of its error is paid. It is noted that the convergence problem is solved only for boundary value problems in the theory of elasticity. Therefore the error assessment in the problem is solved on the basis of a comparison of the approximate analytical solution with the numerical solution constructed on the finite element method, for some special cases. Considering the symmetry of the problem, the finite element model was built for the half tube. The number of finite elements is about 18,000. Considering the symmetry of the problem the second half of the tube is replaced by boundary conditions. Analysis of analytical and numerical solutions is executed depending on the steady-state creep nonlinearity parameter and misaligned parameter that is ratio of the misalignment of the centers of the outer and inner diameter to the outer radius. It is shown that the error
© 2014 Samara State Technical University.
How to cite Reference: Radchenko V. P., Moskalik A. D., Adeyanov I. E. Comparative Analysis of the Approximate Analytical and Finite Element Solutions for Misaligned Tube, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2014, no. 3(36), pp. 79-93. doi: 10.14498/vsgtu1340. (In Russian)
Authors Details: Vladimir P. Radchenko (Dr. Phys. & Math. Sci., Professor; radch@samgtu.ru), Head of Department, Dept. of Applied Mathematics & Computer Science. Anna D. Moskalik (annmoskalik1@gmail.com; Corresponding Author), Postgraduate Student, Dept. of Applied Mathematics & Computer Science. Igor’ E. Adeyanov (Cand. Techn. Sci.; adigorev@gmail.com), Associate Professor, Dept. of Mechanics.
91
Радченко В. П., Москалик А. Д., Адеянов И. Е.
of deviation of the approximate analytical solution in the second approximation from numerical solution until the misalignment value of the centers of the inner and outer diameters of 0.1 for the tubes with small exponent of the steady-state creep (3 to 8) is not more than 9 %, and error to 8 % for the tubes with a large exponent of the steady-state creep nonlinearity is observed in the misaligned parameter to 0.06. Results of computations are presented in tabular form and in the form of graphs. Recommendations for the use of the constructed approximate analytical solution in applied problems are given.
Keywords: thick-walled misaligned tube, steady-state creep, approximate analytical solution, small parameter method, second approximation, finite element method, numerical solution, solution error.
doi: http://dx.doi.org/10.14498/vsgtu1340 ORCID
Vladimir P. Radchenko: http://orcid.org/0000-0003-4168-9660 Anna Moskalik: http://orcid.org/0000-0001-7527-6237 Igor Adeyanov: http://orcid.org/0000-0002-3943-6173
REFERENCES
1. Hill R., Hutchinson J. W. Bifurcation phenomena in the plane tension test, J. Mech. Phys. solids, 1975, vol. 23, pp. 239-264. doi: 10.1016/0022-5096(75)90027-7.
2. Storen S., Rice J. R. Localized necking in thin sheets, J. Mech. Phys. solids, 1975, vol. 23, no. 6, pp. 421-441. doi: 10.1016/0022-5096(75)90004-6.
3. Hutchinson J. W., Neale K. W. Influence of strain-rate sensitivity on necking under uniaxial tension, Acta Metallurgica, 1977, vol. 25, no. 8, pp. 839-846. doi: 10.1016/0001-6160(77) 90168-7.
4. Keller I. E. Self-similar shapes of the free boundary of a nonlinear-viscous band under uniaxial tension, J. Appl. Mech. Tech. Phys., 2010, vol. 51, no. 1, pp. 99-105. doi: 10.1007/ s10808-010-0016-z.
5. Popov N. N., Radchenko V. P. 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. doi: 10.1016/j.jappmathmech.2013.02.011.
6. Dolzhkovoi A. A., Popov N. N. 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 (In Russian). doi: 10.14498/vsgtu102.
7. Popov N. N., Isutkina V. N. 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. 57-61 (In Russian). doi: 10.14498/ vsgtu535.
8. Dolzhkovoi A. A., Popov N. N., Radchenko V. P. Solution of the stochastic boundary-value problem of steady-state creep for a thick-walled tube using the small-parameter method, J. Appl. Mech. Tech. Phys., 2006, vol. 47, no. 1, pp. 134-142. doi: 10.1007/ s10808-006-0017-0.
9. Kovalenko L. V., Popov N. N., Radchenko V. P. Solution of the plane stochastic creep boundary value problem, J. Appl. Math. Mech., 2009, vol. 73, no. 6, pp. 727-733. doi: 10. 1016/j.jappmathmech.2010.01.013.
10. 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. doi: 10.1007/BF00909710.
11. Kachanov L. M. Teoriia polzuchesti [Creep theory]. Moscow, Fizmatgiz, 1960, 455 pp. (In Russian)
92
Сравнительный анализ приближенного аналитического и конечно-элементного решений ...
12. Radchenko V. P., Bashkinova E. V. Solution of the boundary value problems for steady creep in polar coordinates by the perturbation method, Vestn. Sam. gos. tekhn. un-ta. Ser. Tekhn. nauki, 1998, no. 5, pp. 86-91 (In Russian).
13. Bashkinova E. V. Solution of the value boundary problem of steady creep for non-axisymmetric thick-walled tube, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 16, 2002, pp. 105-110 (In Russian). doi: 10.14498/vsgtu106.
14. Ivlev D. D., Ershov L. V. Metod vozmushchenii v teorii uprugoplasticheskogo tela [Perturbation Method in the Theory of an Elastic-Plastic Body]. M., Nauka, 1978, 208 pp. (In Russian)
15. Kerzhaev A. P. Elastoplastic state of the thin annular plate in the presence of translational anisotropy under uniform tension, Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ia. Iakovleva. Seriia: Mekhanika predel’nogo sostoianiia, 2012, no. 2(12), pp. 174-179 (In Russian).
16. Fominykh S. O. Elastoplastic state of the thick-walled pipe by reacting the different types of plastic anisotropy, Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ia. Iakovleva. Seriia: Mekhanika predel’nogo sostoianiia, 2011, no. 1(9), pp. 201-2016 (In Russian).
17. Petrov N. I. On the deformation of stretched strip weakened sloping bores, Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ia. Iakovleva. Seriia: Mekhanika predel’nogo sostoianiia, 2014, no. 2(20), pp. 36-45 (In Russian).
18. Nikitin A. V., Tikhonov S. V. Limit condition anisotropic multilayer translationnally thick-walled pipes under internal pressure, Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ia. Iakovleva. Seriia: Mekhanika predel’nogo sostoianiia, 2014, no. 1(19), pp. 88-94 (In Russian).
19. Kuntashev P. A., Nemirovskii Yu. V. Convergence of the perturbation method in elastic problems, Izv. Akad. Nauk SSSR, Mekh. Tverd. Tela, 1985, no. 3, pp. 75-78 (In Russian).
20. Nikitenko A. F. Polzuchest’ i dlitel’naia prochnost’ metallicheskikh materialov [Creep and Long-Term Strength of Metal Materials]. Novosibirsk, Novosibirsk State University of Architecture and Civil Engineering, 1997, 278 pp. (In Russian)
21. Rabotnov Iu. N. Mekhanika deformiruemogo tverdogo tela [Mechanics of deformable solids]. Moscow, Nauka, 1979, 744 pp. (In Russian)
22. Moskalik A. D. Analysis of the stress-strain state of a thick-walled misalignment cylinder under internal pressure in steady-state creep by small parameter method, Proceedings of the Ninth All-Russian Scientific Conference with international participation. Part 1, Matem. Mod. Kraev. Zadachi. Samara, Samara State Technical Univ., 2013, pp. 140-144 (In Russian).
23. Moskalik A. D. The application of perturbation method to problem of misaligned tube in conditions of steady-state creep, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2013, no. 4(33), pp. 76-85 (In Russian). doi: 10.14498/vsgtu1290.
24. Malinin N. N. Prikladnaia teoriia plastichnosti i polzuchesti [Applied Theory of Plasticity and Creep]. Moscow, Mashinostroenie, 1975, 400 pp. (In Russian)
25. Radchenko V. P., Saushkin M. N. Polzuchest’ i relaksatsiia ostatochnykh napriazhenii v uprochnennykh konstruktsiiakh [Creep and Relaxation of Residual Stresses in Hardened Structures]. Moscow, Mashinostroenie-1, 2005, 226 pp. (In Russian)
Received 26/V/2014;
received in revised form 18/VI/2014;
accepted 27/VIII/2014.
93
