Научная статья на тему 'Неустановившееся движение жёсткопластического осесимметричноного тела по внутренней поверхности конуса (обжатие)'

Неустановившееся движение жёсткопластического осесимметричноного тела по внутренней поверхности конуса (обжатие) Текст научной статьи по специальности «Математика»

CC BY
92
29
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕУСТАНОВИВШЕЕ РАВНОВЕСНОЕ ДВИЖЕНИЕ / ЖЁСТКОПЛАСТИЧЕСКОЕ ОСЕСИММЕТРИЧНОЕ ТЕЛО / КОНИЧЕСКАЯ МАТРИЦА / ЖЁСТКОПЛАСТИЧЕСКИЙ МАТЕРИАЛ / ПРОЦЕСС ОБЖАТИЯ / ДИСКРЕТНЫЕ УРАВНЕНИЯ

Аннотация научной статьи по математике, автор научной работы — Маркин Алексей Александрович, Лыу Туан Ань

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Неустановившееся движение жёсткопластического осесимметричноного тела по внутренней поверхности конуса (обжатие)»

Известия Тульского государственного университета Естественные науки. 2014. Вып. 3. С. 89-97 Механика

УДК 539.3

Неустановившееся движение жёсткопластического осесимметричноного тела по внутренней поверхности конуса

(обжатие) *

А. А. Маркин, Лыу Туан Ань

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

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

Введение

В работе [1] дано описание установившегося процесса движения тонкого тела по конической поверхности в рамках жёсткопластической модели Треска. В статьях [2, 3] получена общая система уравнений, описывающая неустановивщееся движение осесимметричных тел по поверхностям в эйлеровых координатах при условии пластичности мизеса. В статье [4] определены характеристики напряженно-деформированного состояния заготовки в процесс её раздачи при нестационарном течении по конической матрице.

В данной статье рассматривается нестационарное течение по конической матрице в процессе обжатия.

Постановка задачи

Пусть заготовка в начальном состоянии находится внутри конической поверхности, имеет постоянную толщину Ло. Положения верхнего и нижнего

* Работа выполнена при финансовой поддержке Министерства образования и науки РФ (госзадание № 467).

торцев заготовки в начальном состоянии задаются координатами (0)5*1 и тогда 1о = (0)52 - (0)^1 — начальная длина заготовки. Нижний торец движется к точке О вверху со скоростью VI, а верхний торец не нагружен.

Рис. 1. Течение материала по внутренней поверхности конуса

Представление системы уравнений в безразмерном виде

В работе [4] система уравнений, описывающая движение заготовки, приведена к следующему безразмерному виду:

dj Л sin y ( 1 1 dh \ ( 2п \

- Y - 1Л +1 — íWy + 2t , (1)

dx S cos (y + VS h dx J \ 3 /

dV = Л у cos Y (2) dx S cos (y + dh hdV hTr

дт = даХ + SSV (3)

dV = - дЛ

дх дт dS (x, т)

(4)

= -V (x,t), (5)

дт

где х = (50 - (1)50) /10 ; т = ((2)50 - (2)5) /10 .

Таким образом, имеем систему пяти нелинейных уравнений (1)—(5) в частных производных относительно пяти неизвестных. Сформулируем граничные условия для случая, когда задано движение нижнего торца (2)50, а верхний торец свободен от напряжений, тогда

V (x) lx=1 = -1,

(6)

п

Y (x)|x=0 = о

I x=0 н0 "11 lx=0

Принимаем следующие начальные условия:

°"iilx=o = 0-

h|T=o = ho ; dX

= Л|т=0 = 1 ; S |т=0 = So = x + a.

(7)

(8

т=0

Подставим начальные условия в уравнение (3). В результате получим

дТс

dx

sin Y0

x + a cos (y0 + )

1 / 2п

—— tg Y + x + a V 3

(9)

Уравнение (9) с граничным условием

п

7о (х)|х=о = 3

позволяет найти распределение угла вида 70 (ж) в начальный момент.

Из (2) с учетом (8) следует уравнение, определяющие поле скоростей Ро (ж) в начальный момент:

dV) V0

cos Y0

(10)

dx x + a cos (Y0 + 2t)' Граничное условие V0|x=1 = -1.

Если из (9), (10) найдены функции Y0 (x) и V0 (x), то из (2) и (3) находим производные по «времени» от кинематических характеристик в начальный момент:

дл

дТ

dh дТ

т=0

т=0

. , dV0 h0

= h^--+ — V0

dx S0

dV0; dS

dx ' дт

= -V0 (x)

т=0

Дискретизация системы (1)—(8) по переменной х

Для уточнения представления производных от скорости по координате используем неявную схему дискретизации, тогда

дx

1 / Vi+1 - Vi-1

2 I Ax

(11)

Использование неявной схемы для дискретизации приводит к системе нелинейных уравнений относительно узловых перемещений поэтому будем определять 7^ (ж) из уравнения (1) по явной схеме. В результате

1

x=x

сочетания явной и неявной схем система (1)—(5) сводится к следующей системе обыкновенных дифференциальных уравнений по «времени» относительно узловых величин:

АгДж sin(Yi_i) / АгДж hi+1 - hi-Л / 2п\

Yi = Yi-1 cos (7i_i + f ) Л^T + MYi-1 + '

(12)

V V = 2А*ДхУ cos Yi (13)

Vi+1 - Vi-1 = ~ Vi cos (Yi + f) ' (13)

^ = -Vi, (14)

ar

dAi 1 V+1 - Vi_1

> (15)

dr 2 Дх

dhi hi Vi+1 - Vi_1 hi

df = Ai 2Дх + SiVi- (16)

Система (12)-(16) интегрируется при следующих краевых условиях: Yo (т) = 3 — угол вида в начальной точке х = 0 соответствует отсутствию напряжения ац на верхнем торце;

РП (т) = -1 — граница детали х = 1 движется с постоянной единичной скоростью.

Начальные условия для системы (12)-(16) имеют вид: hi|T=o = ho — толщина детали в начальный момент постояна; Ai|T=o = 1 — деталь в начальный момент не деформирована; Si|T=o = Soi = Дх.г — текущие координаты в начальный момент совпадают с начальными.

Дискретизация системы (12)—(16) по «времени»

Разобъем интервал по «времени» на т отрезков Дт. Значение неизвестных узловых функций в момент т = Дтк обозначим Y; Л^;

Здесь к = 1,2,..., т. Уравнения (12), (13) запишем в «момент» к — 1. Это позволит при к = 1 учитывать начальные условия. Дискретная форма уравнений (12), (13) принимает следующий вид:

_ _ Лк—МАж 8Ш 7к-М-1

7к— 1 ,г 7к— 1,г— 1 0 / , 2п \

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ьк-1,г 008 (7к—1,г—1 + 2п) Лк;—1,гАж . 1,г+1 — 1,г— 1 \ , / . 2п

+ --!--1,г— 1 + "Г" ,

^к—1,г 2Лк;— 1,г ) \ 3 / (17

П .

7к—1,0 _ 3; £о,г _ Аж.г + а;

Ло,г _ Ло _ ооп81; _ ; ^,„+1 _ Лк,га;

2АхЛк— М 008 7к_1,г

^к—1,г+1 - ^к—1,г—1 _ —-7-^к—

^к—1,г 008 (7к— М +

Рк—1,„ _ -1; Рк—1,„+1 _ -1; Ло,г _ 1; £о,г _ Аж.г + а.

(18)

Отметим, что в уравнения (17) и (18) входят начальные и краевые условия.

Для интегрирования эволюционных уравнений (14)—(16) используем явную схему Эйлера. В результате получаем:

_ ^к— 1,г — Атрк— 1,г,

£о,г _ Аж.г + а, 1 Ат

Лк,г _ Лк—1,г — 2 д^ 1,г+1 — ^к—1,г—0 ,

Ло,г _ 1

, , , , л /^—1,1+1 — ^к—1,г—1 , ^к—1,г

_ Лк—1,г + Л:—1,гАт -—----+

(19)

(20)

Ло г _ Ло _ 00П81 .

(21)

Последовательность решения системы уравнений (17)—(20)

Решение строится на последовательности шагов по «времени». Для каждого значения к _ 1,2,...,т находятся соответстующие значения искомых функций в узловых пространственных точках ж. Для построения решения в начальный момент принимаем к _ 1. В этот момент уравнение (17) принимает с учетом начальных условий следующий вид:

Дх йШ Y0,г-l Дх / 2п\

Y0,г = Y0,г-1 - д-Г";--т-- "Т-—- Y0,г-1 + "Г"

' ' Дх.г + а сой ^0,г-1 + 2т) Дхд + а \ 3 /

п

Yo,o

3

(22)

Из (22) определяем последовательно Y0,г для г = 1,2,..., п. Зная Yo,г, определяем узловые значения скорости при к = 1 из (18):

2Дх сой Y0,' (Дх.г + а) сой + 2т)

^0,г+1 — = —-;---—!-2ПГ

, . .. , (23)

^0,«. = —1; ^0,га+1 = - 1

Для определения Т/0,г необходимо решить систему линейных уравнений (23). Определив У0,г, переходим к определению £1>г; Л1>г; А1>г. Из (19)-(21) при к = 1 получим

51,» = 50,г - ДтУ0,г = Дх.г + а - Дт^г, (24)

1 Дт

А1,г = 1 - ^ ДХ (Чг+1 - ) , (25)

Л1,г = Л0,г + Л0,гДт —^-т-—^--+ —^ . (26)

\ 2АМДх 51,»)

После определения 51>г; Л1>г; А1>г полагаем к = 2 и переходим снова к решению уравнений (17) и (18), а затем к (19) и (20). Далее цикл продолжается до заданого к = т момента «времени».

На рис. 2-5 приведены зависимости характеристик напряженного-деформированного состояния от координаты X для последовательных значений параметра т. При вычислении Дх принималось равным 10-3 и Дт = 3Дх. В качестве лагранжевой переменный принимается координата X = 1 - х с началом отсчета в точке А.

Рис. 2. Закон изменения толщины заготовки при различных шагах

«времени»

Рис. 3. Закон изменения напряжения 7ц при различных шагах «времени»

Рис. 4. Закон изменения напряжения <22 при различных шагах «времени»

Из графиков, показанных на рис. 2 следует, что толщина заготовки с ростом «времени» увеличивается, причем толщина верхнего торца увеличивается больше чем нижнего. Из графиков, показанных на рис. 3 и рис. 4 следуют, что напряжения <ц и <22 сжимающие. При этом с ростом «времени» они изменяются незначительно.

Рис. 5. Закон изменения нормального давления при различных шагах

«времени»

Из графиков, показанных на рис. 5, следует, что с ростом «времени» нормальное давление увеличивается, причём в окрестности нижнего торца меньше, чем в окрестности верхнего торца.

Список литературы

1. Ильюшин А.А. Труды (1946-1966). Т. 2. Пластичность. М.: Физматлит, 2004. 480 с.

2. Маркин А.А., Лыу Туан Ань. Движение тонкого жёсткопластического тела по поверхности с осевой симметрией // Изв. ТулГУ. Естественные науки. 2012. Вып. 3. С. 93-101.

3. Маркин А.А., Лыу Туан Ань. Движение тонкого жёсткопластического тела по конической и тороидальной поверхностям // Изв. ТулГУ. Естественные науки. 2013. Вып. 3. С. 154-156.

4. Маркин А.А., Лыу Туан Ань . Нестационарное течение жёсткопластического осесимметричноно тела по конической поверхности // Изв. ТулГУ. Естественные науки. 2014. Вып. 2. С. 164-173.

Маркин Алексей Александрович (markin-nikram@yandex.ru), д.ф.-м.н., профессор, зав. кафедрой, кафедра математического моделирования, Тульский государственный университет.

Лыу Туан Ань (myname_is_luutuananh@yahoo.com), аспирант, кафедра математического моделирования, Тульский государственный университет.

The nonstationary motion of a rigid-plastic axisymmetric body on a inner surface of a cone (compression)

A. A. Markin, Luu Tuan Anh

Abstract. Based on system of equations describing nonstationary motion of axisymmetric body, in which the body material is rigid-plastic incompressible, to view the change characteristics of the strain state and the angle of view stressed during transition of material on the conical matrix.

Keywords: velocity field gradient, velocity deformation tensor, equilibrium motion, rigid-plastic materials, conical matrix, transitional motion, loop-system.

Markin Alexey (markin-nikram@yandex.ru), doctor of physical and mathematical sciences, professor, head of department, department of mathematical modelling, Tula State University.

Lyu Tuan An' (myname_is_luutuananh@yahoo.com), postgraduate student, department of mathematical modelling, Tula State University.

Поступила 09.06.2014

i Надоели баннеры? Вы всегда можете отключить рекламу.