Аналитические решения квазистатических задач термоупругости с переменными физическими свойствами среды

Тип работы:
Реферат
Предмет:
Механика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Вестн. Сам. гос. техн. ун-та. Сер. Физ. -мат. науки. 2014. № 2 (35). С. 130−135
УДК 517. 958:539. 3(4)
АНАЛИТИЧЕСКИЕ РЕШЕНИЯ КВАЗИСТАТИЧЕСКИХ ЗАДАЧ ТЕРМОУПРУГОСТИ С ПЕРЕМЕННЫМИ ФИЗИЧЕСКИМИ СВОЙСТВАМИ СРЕДЫ
В. А. Кудинов, А. Э. Кузнецова, А. В. Еремин, Е. В. Котова
Самарский государственный технический университет,
Россия, 443 100, Самара, ул. Молодогвардейская, 244.
Используя ортогональной метод Бубнова-Галёркина, на основе разработанного метода построения систем координатных функций, в любом приближении точно удовлетворяющих неоднородным граничным условиям, получено высокой точности приближённое аналитическое решение нелинейной квазистатической задачи термоупругости для бесконечного полого цилиндра с переменными по '-радиальной координате физическими свойствами. Математическая постановка задачи включает нелинейные уравнения относительно искомой функции перемещения и неоднородные граничные условия. Решение разыскивается в таком виде, чтобы оно заранее точно удовлетворяло граничным условиям задачи. Точное выполнение граничных условий осуществляется благодаря использованию координатных функций особой конструкции. Неизвестные коэффициенты решения находятся путём составления невязки исходного дифференциального уравнения и выполнения требования ортогональности невязки ко всем координатным функциям. Отсюда относительно неизвестных коэффициентов решения получается система алгебраических линейных уравнений, число которых равно числу приближений принятого решения. Показано, что с увеличением числа приближений точность решения существенно возрастает. Так, уже в девятом приближении невязка исходного дифференциального уравнения равна нулю практически во всем диапазоне изменения пространственной переменной. Максимальная невязка в шестом приближении составляет е = 5 • 10−4.
Ключевые слова: термоупругость, переменные физические свойства среды, аналитическое решение, координатные функции, ортогональный метод Бубнова- Галёркина.
При скорости изменения температуры, значительно меньшей скорости распространения звука в материале, температурное поле практически не зависит от вызываемых им деформаций. В этом случае можно не учитывать зависящее от деформации слагаемое во взаимосвязанном уравнении тепло-
ISSN: 2310−7081 (online), 1991−8615 (print) — doi: http: //dx. doi. org/10. 14 498/vsgtu1219 © 2014 Самарский государственный технический университет.
Образец цитирования: В. А. К у д ин о в, А. Э. К у з не цов а, А. В. Е ре м ин, Е. В. К о т о в а, «Аналитические решения квазистатических задач термоупругости с переменными физическими свойствами среды» // Вестн. Сам. гос. техн. ун-та. Сер. Физ. -мат. науки, 2014. № 2(35). С. 130−135. doi: 10. 14 498/vsgtu1219.
Сведения об авторах: Василий Александрович Кудинов (д.ф. -м.н., проф.), заведующий кафедрой, каф. теоретических основ теплотехники и гидродинамики. Анастасия Эдуардовна Кузнецова, аспирант, каф. теоретических основ теплотехники и гидромеханики. Антон Владимирович Еремин (к.т.н.), старший преподаватель, каф. теоретических основ теплотехники и гидромеханики. Евгения Валериевна Котова (к.т.н.), старший преподаватель, каф. теоретических основ теплотехники и гидромеханики.
E-mail addresses: totig@yandex. ru (V.A. Kudinov), kuznetsovaae@rambler. ru (A.E. Kuznetsova), a.v. eremin@list. ru (A.V. Eremin, Corresponding author), larginaevgenya@mail. ru (E.V. Kotova)
130
Аналитические решения квазистатических задач термоупругости …
проводности и инерционные члены в уравнениях равновесия.
Математические постановки задач, в которых не учитываются указанные составляющие термоупругости, называются квазистатическими. В данном случае время теряет смысл переменной и становится параметром, а напряжения определяются исходя из температурного поля в данный конкретный момент времени. В случае учёта указанных факторов задачи термоупругости называются динамическими. Математические постановки таких задач включают сложные неоднородные гиперболические уравнения, методы решения которых пока ещё недостаточно разработаны.
Исходная математическая постановка квазистатической задачи термоупругости включает: уравнения равновесия- уравнения совместности деформаций- соотношения между деформациями и радиальными перемещениями (геометрические уравнения) — формулы закона Гука (физические уравнения). Если температура изменяется лишь в радиальном направлении (плоская осесимметричная задача), краевую задачу термоупругости можно свести к задаче о плоском напряжённом состоянии (тонкий круглый диск) или к задаче о плоской деформации (бесконечный цилиндр). При этом исходная математическая постановка задачи, включающая систему дифференциальных уравнений, приводится к одному уравнению относительно искомой функции напряжения или перемещения с соответствующими граничными условиями.
Рассмотрим задачу термоупругости для длинного полого цилиндра в случае его плоской деформации, обусловленной плоским осесимметричным температурным полем. Математическая постановка задачи относительно искомой функции перемещения с переменными по радиальной координате физическими свойствами среды в данном случае имеет вид [1−4]
d
dr
«dU Er-) + dr
v dE E 1 + v d
----------U = -----r-(aET),
1 — v dr r) 1 — v dr
(1 — v)+ vU — (1 + v) aT dr r
r=r 1, r=r 2
0,
(1)
(2)
где U = U® — радиальное перемещение- r — радиальная координата- E = = E® — модуль упругости- a = a® — коэффициент линейного расширения- v — коэффициент Пуассона- T = T® — температура.
Согласно граничному условию (2), радиальные напряжения ar ® в точках r = ri и r = r2 равны нулю.
При известных перемещениях радиальные и окружные напряжения определяются по формулам
?Х r

E
(1 + v)(1 — 2v) E
(1 + v)(1 — 2v)
(1 — v) dU + vU — (1 + v) aT~|, dr r
'- dU. U, _i
v- + (1 — v)-----(1 + v) aT ,
dr r
(3)
(4)
где ar = ar ® — радиальное напряжение- ® — окружное напряжение.
Трудности получения точного аналитического решения задачи (1), (2) связаны с её нелинейностью. Рассмотрим способ нахождения её приближённого аналитического решения, согласно которому решение задачи (1), (2) при-
131
В. А. Кудинов, А. Э. Кузнецова, А. В. Еремин, Е. В. Котова
нимается в виде [4]
П
и (г) = Ф (г) + ^ qk & lt-fk ®, (5)
k=2
где Ф (г) = Ai + А2Г — функция, удовлетворяющая неоднородным граничным условиям (2) — A1, A2 — неизвестные коэффициенты- ^k® = 1 + B1r + B2rk — координатные функции, удовлетворяющие однородным граничным условиям (2) (при равенстве нулю слагаемых, содержащих температуру) — B1, B2 — неизвестные постоянные- n — число приближений.
Для определения постоянных Ai и A2 используются неоднородные граничные условия (2). Подставляя функцию Ф (г) в (2), относительно постоянных Ai и A2 получаем следующую систему двух алгебраических линейных уравнений:
A2(1 — v) + v Ai + AГ1 — (1 + v) a (ri)T (ri) = 0, ri
A2(1 — v) + vAi + Ar2 — (1 + v) a (r2)T (r2) = 0. r2
Из решения системы уравнений (6) находим
Ai
A2
rir2(1 + v)(a (ri)T (ri) — a (r2)T (r2)) v (ri — r2)
(1 + v)(ria (ri)T (ri) — r2a (r2)T (r2)) ri — r2 '
Для определения постоянных Bi и B2 используются граничные условия (2). Положив последний член соотношения (2) равным нулю, относительно неизвестных Bi и B2 получаем следующую систему двух алгебраических линейных уравнений:
(Bi + fcB2rk-^ (1 — v) + (Bi + kB2rk-1)(1 — v) +
v (1 + Biri + B2rk) ri
v (1 + Bir2 + B2rk) r2
0,
0.
(7)
Из решения системы уравнений (7) находим
B = _ v (rf — rk) B = _________v (r2 — ri)_________
1 r1r2(Гг1 — rj0−1)' 2 r1r2(Гт1 — rk-1)(k — kv + v)
После нахождения неизвестных постоянных Ai, A2, Bi, B2 соотношение (5) будет точно удовлетворять граничным условиям (2). Неизвестные коэффициенты qk, k = 2, 3,…, n, находятся из выполнения уравнения (1). Для этого составляется его невязка и требуется ортогональность невязки ко всем координатным функциям ipj®, j = 2, 3,…, n:
dU / dU (v dE
— Er- ±--------
dr dr 1 — v dr
E) U — r-^- (aET)
r 1 v dr
(r)dr = 0. (8)
132
Аналитические решения квазистатических задач термоупругости …
г, mm
Рис. 1. Эпюры радиальных напряжений (расчёт по формуле (3)): 1 -третье приближение- 2 — шестое приближение- 3 — девятое приближение
[Figure 1. The diagrams of the radial stresses (calculation by formula (3)): 1 — the third approximation- 2 — the sixth approximation- 3 — the ninth approximation]
Рис. 2. Эпюры окружных напряжений (расчет по формуле (4)): 1-третье приближение- 2 — шестое приближение- 3 — девятое приближение
[Figure 2. The diagrams of the circumferential stresses (calculation by formula (4)): 1 — the third approximation- 2 — the sixth approximation- 3 — the ninth approximation]
Рис. 3. Невязка уравнения (1): 1 — шестое приближение- 2 — девятое приближение [Figure 3. The disparity of equation (1): 1 — the sixth approximation- 2 — the ninth approximation]
133
В. А. Кудинов, А. Э. Кузнецова, А. В. Еремин, Е. В. Котова
Подставляя (5) в (8), после определения интегралов относительно неизвестных коэффициентов qk будем иметь систему n — 2 алгебраических линейных уравнений. После определения из решения этой системы неизвестных коэффициентов qk, к = 2, 3,…, n, приближённое аналитическое решение задачи (1), (2) находится из (5).
В качестве конкретного примера найдём решение задачи термоупругости для полого цилиндра при следующих исходных данных: п = 14 мм- Г2 = = 39 мм- v = 0. 2- E = 19.5 ¦ 109 кг/м2- a = 10 ¦ 10−6 К-1. Распределение температуры по толщине слоя цилиндра принималось в виде
T® = 170 — 6592r + 13 8308r2 — 1. 7402 ¦ 106r3.
Результаты расчётов представлены на рис. 1, 2. Их анализ позволяет сделать заключение о том, что как для радиальных, так и для окружных напряжений отмечается незначительное различие результатов шестого и девятого приближений. Следовательно, уже в девятом приближении полученное решение практически совпадает с точным. Оценка невязки уравнения (1) для шестого и девятого приближений подтверждает данное заключение (см. рис. 3). Как видно из рис. 3, невязка уравнения (1) в девятом приближении практически во всем диапазоне изменения переменной r близка к нулю.
Заключение. Разработана методика получения высокоточных приближенных аналитических решений квазистатических нелинейных задач термоупругости (плоская деформация, плоское напряжённое состояние) с переменными физическими свойствами среды. Решение, найденное с помощью ортогонального метода Бубнова-Галёркина при использовании полученной в данной работе системы координатных функций в любом приближении точно удовлетворяющих неоднородным краевым условиям задачи, представляет быстро сходящийся ряд. Так, уже в девятом приближении полученное решение практически совпадает с точным, что подтверждается оценкой невязки основного дифференциального уравнения краевой задачи. Отметим, что ввиду нелинейности этого уравнения точные аналитические решения подобных задач в настоящее время не получены.
Работа выполнена при поддержке Минобрнауки России в рамках базовой части государственного задания ФГБОУ ВПО «СамГТУ». Код проекта: 1273.
This work was supported by the Russian Ministry of Education and Science within the base portion of the state task to Samara State Technical University, project no. 1273.
СПИСОК ЛИТЕРАТУРЫ/ REFERENCES
1. B. A. Boley, J. H. Weiner, Theory of thermal stresses, New York, John Wiley, 1960, xvi+586 pp.- Б. Боли, Дж. Уэйнер, Теория температурных напряжений. М.: Мир, 1964. 517 с.
2. А. Д. Коваленко, Введение в термоупругость. Киев: Наукова думка, 1965. 202 с. [A. D. Kovalenko, Vvedenie v termouprugost' [Introduction to Thermoelasticity], Kiev, Naukova Dumka, 1965, 202 pp. (In Russian)]
3. S. P. Timoshenko, J. Goodyear, Theory of Elasticity, New York, McGraw-Hill, 1970- С. П. Тимошенко, Дж. Гудьер, Теория упругости. М.: Наука, 1979. 560 с.
4. В. А. Кудинов, Э. М. Карташов, В. В. Калашников, Аналитические решения задач тепломассопереноса и термоупругости для многослойных конструкций. М.: Высшая школа, 2005. 430 с. [V. A. Kudinov, E. M. Kartashov, V. V. Kalashnikov, Analiticheskie resheniya zadach teplomassoperenosa i termouprugosti dlya mnogosloynykh konstruktsiy
134
Аналитические решения квазистатических задач термоупругости …
[Analytical problem solving heat and mass transfer and thermoelasticity for multilayer structures], Moscow, Vysshaya shkola, 2005, 430 pp. (In Russian)]
Поступила в редакцию 01/IV/2013- в окончательном варианте — 05/XII/2013- принята в печать — 19/III/2014.
MSC: 74F05- 74B20
ANALYTICAL SOLUTIONS OF THE QUASISTATIC THERMOELASTICITY TASK WITH VARIABLE PHYSICAL PROPERTIES OF A MEDIUM
V. A. Kudinov, A. E. Kuznetsova, A. V. Eremin, E. V. Kotova
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443 100, Russian Federation.
The high-precision approximate analytic solution of the nonlinear quasi-static problem of thermoelasticity for an infinite hollow cylinder with variable along the radial coordinate physical properties is obtained, using the orthogonal Bubnov-Galyorkin method developed by the construction of systems of coordinate functions exactly satisfying inhomogeneous boundary conditions in any approximation. The mathematical formulation includes non-linear equations for the unknown function of displacement and inhomogeneous boundary conditions. The desired solution is supposed to precisely satisfy the boundary conditions in advance. The exact fulfillment of the boundary conditions is achieved using the coordinate functions of special design. The unknown coefficients are found by constructing the disparity of original differential equation, that should be orthogonal to all the coordinate functions. Hence, the unknown coefficients of solution yields a system of linear algebraic equations, which number is equal to the number of approximations of the solution. It is shown that the solution accuracy increases substantially with increasing the number of approximations. Thus, already in the ninth approximation the disparity of original differential equation is zero almost the entire range of the spatial variable. The maximum disparity in the sixth approximation is e = 5 • 10−4.
Keywords: thermoelasticity problem, variable physical properties of a medium, analytical solution, coordinate functions, Bubnov-Galyorkin orthogonal method.
Received 01/IV/2013-
received in revised form 05/XII/2013-
accepted 19/III/2014.
ISSN: 2310−7081 (online), 1991−8615 (print) — doi: http: //dx. doi. org/10. 14 498/vsgtu1219 © 2014 Samara State Technical University.
Citation: V. A. Kudinov, A. E. Kuznetsova, A. V. Eremin, E. V. Kotova, «Analytical Solutions of the Quasistatic Thermoelasticity Task with Variable Physical Properties of a Medium», Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz. -Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & amp- Math. Sci. ], 2014, no. 2 (35), pp. 130−135. doi: 10. 14 498/vsgtu1219. (In Russian) Authors Details: Vasiliy A. Kudinov (Dr. Phys. & amp- Math. Sci.), Head of Dept., Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics. Anastasiya E. Kuznetsova, Postgraduate Student, Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics. Anton V. Eremin, Assistant (Cand. Techn. Sci.), Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics. Eugeneya V. Kotova, Assistant (Cand. Techn. Sci.), Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics.
E-mail addresses: totig@yandex. ru (V.A. Kudinov), kuznetsovaae@rambler. ru (A.E. Kuznetsova), a.v. eremin@list. ru (A.V. Eremin, Corresponding author), larginaevgenya@mail. ru (E.V. Kotova)
135

ПоказатьСвернуть
Заполнить форму текущей работой