К вопросу о разработке математической модели изотермического стесненного деформирования анизотропной листовой заготовки в прямоугольную матрицу Текст научной статьи по специальности «Физика»

УДК 621.983; 539.374

К ВОПРОСУ О РАЗРАБОТКЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ИЗОТЕРМИЧЕСКОГО СТЕСНЕННОГО ДЕФОРМИРОВАНИЯ АНИЗОТРОПНОЙ ЛИСТОВОЙ ЗАГОТОВКИ В ПРЯМОУГОЛЬНУЮ МАТРИЦУ

С.Н. Ларин, В.И. Платонов

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

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

Прямоугольные элементы, являющиеся частями многослойных листовых конструкций, получают посредством пневмоформовки внутренних листов с наружными листами до полного их прилегания. Проанализируем ту стадию, на которой реализуется стесненное деформирование. В работе будет исследовано формирование элементов в углах оболочки в плоскостях уог и хог. Учтем, что Ь > а > Н1. Принимаем, что материал заготовки - анизотропен; учитываем, что главные оси анизотропии совпадают с осями координат х, у, г. Большая сторона заготовки перпендикулярна направлению прокатки х. Известны давление Р1, высота оболочки Щ, накопленная повреждаемость Юд и распределение толщины оболочки = (ф) в момент ? = ?1, где ф - угол, характеризующий положение точки на угловом элементе заготовки. После контакта вершины купола с обшивкой принимаем, что осуществляется то состояние, при котором толщина оболочки меняется одинаково в каждой точке оболочки от начальных размеров при ? = ^, а форма деформируемой угловой части оболочки в плоскости хог сохраняет форму части окружности, а в плоскости уог -сначала форму части эллипса с последующим переходом в форму части окружности.

Разобьем стесненное деформирование на два этапа. На первом этапе в плоскостях уог и хог реализуется плоский участок в районе вершины

купола до того времени, когда 51 = 51* = Ь - Щ и Sз = 53* = а - Щ соответственно. На втором этапе реализуется симметричное формообразование оболочки относительно новых осей симметрии О2О2 и О3 О3 с образова-

183

нием симметрично плоских участков в угловой части оболочки; при этом форма деформируемой свободной угловой части в указанных выше плоскостях имеет форму части окружности (рис. 1 и 2). Примем, что на первом этапе в плоскости уог форма эллипса сохраняется, его полуось ОС не изменяется, оставаясь равной Н1, а полуось Od изменяется от размера Ь до размера Н1, после чего реализуется второй этап второй стадии деформирования с изменением формы оболочки на часть окружности.

Рис. 1. Формообразование угловых Рис. 2. Формообразование угловых элементов в плоскости уо2 элементов в плоскости хог

Рассмотрим возможные ситуации на втором этапе деформирования:

1) 51 £ 51*, £3 £ £3* ;

2) 51 £ 51*,53 > 53*(54 > 0);

3) 51 > 51* (52 > 0), 53 > 53* (54 > 0).

Рассмотрим два деформированных состояния на первом этапе в плоскости уог: первое с длиной средней линии свободной поверхности оболочки в виде эллипса Ьо и длиной участка контакта £1 и второе с длиной средней линии свободной поверхности оболочки Ь и длиной участка контакта 51 + dS1. При переходе из первого состояния во второе приращение меридиональной деформации определится так:

184

йе

р [1,5(6 - 51 - й51 + И1 ]+ 51 + й51

4

У

Р

4

1,5(6 - + И1) -^1 (Ь - 51)И1 ]- 51 1,5(6 - 51 + И1) (Ь - 51) И1 ] + 51

р 4

р 4

1,5(Ь - 51 + И1) (Ь - 51)И1 ] - 51

(1)

Запишем выражение, которое приближенно позволяет определить периметр эллипса

Ь = р[1,5(Ь + И1) -■^ЬИ~1], (2)

где Ь и И - полуоси у получаемого эллипса.

В выражении (1) приведем к подобию члены, и воспользовавшись

разложенной

в ряд функций л/Т-Х = 1 - —

йе у =

2

0,5 +1 2Ц

Х, получим

И

Ь - 5

1)

й51

[1,5(Ь - 51 + И1) (Ь - 51) И1 ] + - 51

р

(3)

Скорость деформации в данном случае оценивается по соотноше-

нию

0,5 + ^Щ/(Ь - 51)

й5

1

(4)

'у 1,5(Ь - 51 + И1 (Ь - 51)И1 + 451 / р й '

Для дальнейшего решения задачи выделим и оценим два близких состояния заготовки на первом этапе рассматриваемой стадии формовки в

/ / о

плоскости х о^ , которая параллельна хог с условием у = 5^ первое состояние - радиус нейтрального слоя р и длина зоны контакта 5з и второе - когда радиус нейтрального слоя р + йр и длина зоны контакта 5з + й5з. Примем, что 5з < а - И.

Приращение окружной деформации в условиях перехода из первого во второе состояние определяется в следующем виде:

р Хйа + йр Ха + й5з

а

где tg — ;

И

a - 53

йе Х =■

а = 2агctg- И1

р ха + 53

(5)

а-5

йа = 2

И1й53

3

(а-5з)2 + И2 '

(а - 53)2 + И2 (а - 53)й53

р х =-—^-йр х = --

2 И

1

И

1

Подставив эти выражения в уравнение (5), получим

' a - Н Л 1--- аг^

2 Н

de1

H

1

a - S

3

[(а - S3)2 + H12]arctg—H— + H1S3

^3.

(7)

а - £

3

В случае равенства длины поверхности контакта Sз скорость де-

формации в искомом направлении находится по формуле

' а - S3 Н Л 1--- аг^

2 Н

X

Н

1

а - S

3

dS

3

х

[(а - S3)2 + Щ2]аг^—Н— + H1S3 Л

(8)

а - ^

Можно отметить, что в одно и то же время с образованием зоны контакта Sl на поверхности уог реализуется образование зоны контакта Sз в плоскости хог. Связь между данными поверхностями не известна, что значительно затрудняет анализ процесса.

Учтем, что Sl > $1* = Ь - Щ или S2 > 0, тогда образуется зона контакта параллельно оси О2О2. Исследуем два похожих состояния заготовки в углах на втором этапе рассматриваемой стадии. Приращение деформации в меридиональном направлении во время перехода из первого состояния во второе находится по формуле

р

deCy

dp у — + 2dS2

р

(9)

р у- + 2$ 2 + $

2

1*

Из рис. 2 получим

Р у = Щ - $2, dp у = -dS2 .

Уравнение (9) имеет вид

(10)

г

de

р

2

v 2

dS2

у

у

2

р

2

р

(11)

$2 + Н12 + $1*

и скорость деформации будет определяться как

X у =■

2

р

2

-5

2 -

р

2

р

2.

$ 2 + Н12 + ^

(12)

В/ /

случае, при котором в плоскости х о^ , параллельной хог при у = $1, $1* = Ь - Н\, $3 > а - Н1 или > 0, приращение окружной деформации найдем по формуле

р

dp х — + 2dS4

de

2

х

р

(13)

Р х 2 + ^4 + $3*

Из рис. 2 определим

Рх = Н -$4, dpх =-dS4.

(14)

В результате подстановки (13) и (14) в выражение (11) получим

dS4

deл

' 2л 2

v

у

v

2-р

2 ,

г

р

$4 + Н12 + $3*

(15)

X

2-р

v 2

у

х

/

р

^4.

2

v 2

р

(16)

$ 4 + Н12 + $3*

Когда > 0 Й > $1* = Ь - Н1), > 0 ^ > $3* = а - Н1), значения Xх и Ху находятся по выражениям (12) и (16). Используя уравнение равновесия безмоментной оболочки

о

у

о х

р. к

(17)

Ру Рх

и выражения для определения отношений скоростей деформаций Xсх и Ху,

определим компоненты напряжений

о

у

рру к

1 + Xх*х (1 + ) уХхХу Ру

о х =

V XyЯy (1 + Ях ) + ^х*х*у Рх ;

X хКх (1 + Ку ) + ^ уЯхКу

X уЯу (1 + Ях ) + X хКхК

о

у

(18) (19)

хху

Здесь

Р х =

(а - $3)2 + Н2

2Н

если $3 £ $3* = а - Щ.

1

Рх = Н1 - $4, если > $3*($4 > 0),

\3/2

Р у = (Ь - ^)2 Н12

2

■ + •

у

2

Н4 (Ь - 51 )4 ,

187

если $1 £ $1* = Ь - Н1,

/

г

ру = И1 - 52, если 51 > 51*= Ь - И1(52 > 0).

Запишем выражение для получения отношений скоростей деформаций, предварительно подставив уравнения для получения напряжений в определенной точке 51 = 5{, границы контакта в момент времени t на поверхности уо2:

Xх = А(1). (20)

х у

Запишем выражения для оценки величины радиусов кривизны в плоскостях уо2 и х 012 границы контактных зон, используя полученные ранее выражения:

р = (Ь-51)2. р = (а-5з)2 + И12 , рх = (а-5з)2 + И1 (21) р у И1 ' р Х 2И1 ' р у 2(Ь - 51) , ( )

если 51 £ 51* и 5з £ 5з*

Допустим, что известно значение 51 = 5l(t) с условием t > ^ и, зна-

й51

чит, —-, где ^ - время, когда вершина изделия достигнет обшивки. В Ш

итоге запишем следующую формулу для оценки 5з :

й53 = а($ Ш51. (22)

Л Р2 (53) А

Данное выражение возможно решить методом последовательных приближений численно путем перехода к конечно-разностному уравнению

5з (tй) = 5з (tй-1) + А^) /1((51:\ [51 (tй) - 51 (^-1)]. (23)

^2(53, tn-1)

Само решение возможно с момента времени t2 = tl + Дt. В первом приближении примем, что А^) = A(tl). Значение А(^) известно на основании решенной ранее задачи о свободной формовке. В следующих приближениях определяются значения р Х по выражению (21) и А^п) по

уравнению (22), 53) - по формуле (12). Значение радиуса кривизны р у

определим по (21).

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

если 51 £ 51* и 53 £ 53*.

В случае выполнения условий 51 £ 51*, 53 > 53* (54 > 0) величины ру и р Х будут определяться как

Р

у

= (Ь - ^ Н1 '

Р х = Н1 - $4. -

Р у

Р х_ (Н1 - $4)Н1

(Ь -

(24)

с учетом

В итоге

X у = ВД)

X х =

' 2-р ^

v 2

у

/

р

2

v 2

р

4.

(25)

Я 4 + Н1 2 + $3*

Xх = ).

X у

(26)

Для оценки $4 и ^ при заданном $1 = ) по аналогии с представленным выше образцом получим

^ = А1(: dS±. (27)

Данная формула может быть решена методом последовательных приближений:

$4 (:й) = Я4(:й-1) + 4 (:п) ^^ )] [^ (:й) - $1 (:п-1)]. (28)

^^ 4(:п-1)]

В случае, при котором $1 > Sl*(S2 > 0), $3 > Sз*(S4 > 0), радиусы кривизны Рх и Ру будут определяться при заданном значении $2:

Р х = Н1 - $4; Ру = Н1 - $2, -Рх = . (29)

' Ру Н - Б2

Скорости деформаций на поверхности в точке $2 могут быть определены следующим образом:

X у = ^^ ; X х = ад^. (30)

В итоге

Xх = л2($).

x у

(31)

Выразим уравнение, с помощью которого в дальнейшем можно будет найти $4 и ^ при известной $2 = $2 (:):

^ 4 = а2 (: ) .F3(S 2)

^ (32)

^ ^4(з4)

Это уравнение также может быть решено методом последовательных приближений:

Р3^3(:п)]

Я4 (:п ) = $4 (:п-1) + А2(:п-1)

^ 4(:п-1)]

189

&2(:п) - $2(:п-1)]. (33)

Решение последовательных приближений осуществляется аналогично описанному ранее.

Зная значения X Х и Xу в точке 51 и используя условие несжимаемости, найдем (учитывая, что 51 £ 51* и 53 £ 53*)

= -^1(51) Ш51 - ^(53) ^ (34)

НА А А

й53 Л, ^1(51) й51

При условии —- = А^)—--- получим

А (53) А

-—= -[ A(t) + Цад)—1. (35)

Н ш ш

Определим из представленного выше выражения изменение толщины за время Дt = tn - tn-1:

ДН(п) = Н(^-1)[1 + А^ Щ(50 А^) 51(^), (36)

где Н(tn_г) - толщина заготовки в искомой точке, когда t = tn-1, которая по своей траектории движения попала в точку 5^п) с толщиной купола

Н[ 51 (^)] = Н(^-1) {1 - [ А^) +ВД (51)[51 (^) - 51 (^-1)]}. (3 7) Когда 51 < 51*, 53 > 53* (54 > 0),

Н[ 5^п)] = Н(^-1){1 - [4(0 + ВД(51)[ЗДп) - 5^-1)]}. (38) Когда 51 > 51*(52 > 0), 53 > 53*(54 > 0),

Н[52(tn)] = Н(^-1){1 - [А2^) + ВД^)^) - 52^-1)]} . (39) Рассмотрим теперь случай, когда формовка происходит в плоскости хо2. Решение произведем аналогично представленному ранее. Задача сведется к оценке отношения скоростей деформации в искомой точке 53 = 53, а также поверхности контакта в момент t в Х02.

При условии, что 53 = 53 ^), получим следующее уравнение для оценки 51 :

Ш51 = 1 ^2(53) Ш53

А А^) ^1(51) А По аналогии с предыдущими получим ( при 53 £ 53*, 51 £ 51*)

(40)

51 (^) = 51 (^-1) + [53 (^) - 53 (^-1)]. (41)

А<^) 51(tn-1)]

На каждом последовательном приближении уточняются величины р у, А(^), Ж^).

Если выполняются неравенства 53 > 53*(54 > 0), 51 £ 51*, радиусы кривизны р Х и ру вычисляются по формулам

Р х = Н1 - $4; Р у

= (Ь -

Н1

(42)

X

= 4(:).

у

Воспользуемся зависимостью $4 = S4(t), для того чтобы найти

1 ^4( $ 4) (8 4

(43)

Л А1(:) $1)

Представленную выше зависимость решим методом последовательных приближений:

ЗДп ) = Sl(tn-1) +

1 Ж:п%4(:п) - $4(:п-1)].

(44)

4(:) ад^-1)]

При $3 > Sз*(S4 > 0), $1 > ^^ > 0) величины радиусов кривизны

Р х = Н1 - $4; Р у = Н1 - $2

и

X х = А2(: ).

у

X

Воспользуемся зависимостью $4 = $4 (:), найдем $2, ^: $ 2(:п) = $ 2(:п-1) +

1 ^^)] ^(Ь) - $4(:п-1)]

(45)

Ы:) ^з[ $ 2(:п-1)] Представленную выше зависимость решим методом последовательных приближений.

Получим из условия несжимаемости формулу (при $3 £ $3*,

$1 £ $1*)

I¡к = -^з)-^dS1

к ( ( (

(46)

Учитывая, что ^ ($1)

1 (к к &

¡Б1 = 1 ^

А(:)

найдем

1 +

1

А(:).

F2(Sз)

(47)

Используя конечные разности, выразим формулу для оценки изменения толщины за время А::

Ак = -к(:п-1)

1 +

1

А(:)

) А$з

(48)

где к(:п-1) - толщина оболочки в точке при : = :п-1, которая по траектории попала в точку $3 (:п) - границу контакта заполнителя и обшивки. Толщину в точке Sз(tn) определим по выражению

191

1 +

h[ s3(t„)] = h(t„ _i)|i - L a(í )

При S3 > S3*(S4 > 0), S1 < s1*

1

F2(S3)[S3(tn) - S3(tn-i)]

h[S 4(tn)] = h(tn-i)

1 -

1 +

1

4(t)

При S3 > S3*(S4 > 0), S1 > S1*(S2 > 0)

F4(S4)[S4(tn) - S4(tn-1)]

h[S 4(tn)] = h(tn-1)

1-

1 +

Ä2(t)

F4(S4)[S4(tn ) - S4(tn-1)]

Работа выполнена в рамках грантов РФФИ № № 16-48-710016 и 1608-00020 и гранта администрации Тульской области.

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

1. Ларин С.Н. Изотермическое деформирование элементов листовых конструкций цилиндрического и прямоугольного сечения в режиме кратковременной ползучести // Извести вузов. Машиностроение. 2011. №11. С. 44-50.

2. Ларин С.Н. Пневмоформовка ячеистых панелей из анизотропного материала // Известия Тульского государственного университета. Технические науки. 2010. Вып. 3. С. 51-61.

3. Математическая модель свободного изотермического деформирования анизотропной листовой заготовки в прямоугольную матрицу / С.Н. Ларин [и др.] // Известия Тульского государственного университета. Механика деформированного твердого тела и обработка металлов давлением. 2003. Вып. 1. С. 3-13.

4. Яковлев С.С., Ларин С.Н., Трегубов В.И. Изотермическая пнев-моформовка элементов ячеистых многослойных листовых конструкций из анизотропных высокопрочных материалов в режиме ползучести / под ред. С.С. Яковлева. Тула: Изд-во ТулГУ, 2011. 173 с.

Ларин Сергей Николаевич, д-р техн. наук, проф., mpf-tulaaramhler.ru, Россия, Тула, Тульский государственный университет,

Платонов Валерий Иванович, канд. техн. наук, доц., mpf-tulaaramhler.ru, Россия, Тула, Тульский государственный университет

TO THE PROBLEM OF ISOTHERMAL DEVELOPMENT OFMÄTHEMÄ TICÄL MODELS CRAMPED DEFORMATION ANISOTROPIC SLAB IN A RECTANGULAR MATRIX

S.N. Larin, V.I. Platonov 192

>

1

>

The multilayer sheet structures rectangular elements get pnevmoformovkoy sheets previously rigidly connected on a path from the outside-governmental sheets. In this connection the urgent problem considered in the paper, andput-yuschayasya in assessing the impact of process parameters on the steady flow of pro-process. The results of the modeling process isothermal uneasy th deformation anisotropic slab into a rectangular matrix to-torye allow in the future to assess the strength and defect under different conditions of deformation-tions.

Key words: pnevmoformovka constrained deformation, stress, strain, force.

Larin Sergey Nikolaevich, doctor of technical sciences, professor, mpf-tulaarambler. ru, Russia, Tula, Tula State University,

Platonov Valeriy Ivanovich, candidate of technical sciences, associate professor, mpf-tulaa rambler. ru, Russia, Tula, Tula State University

УДК 621.983

МЕТРОЛОГИЧЕСКИЙ МОНИТОРИНГ ИЗДЕЛИЙ, ПОЛУЧАЕМЫХ

РОТАЦИОННОЙ ВЫТЯЖКОЙ

В.И. Трегубов, Д.И. Благовещенский, Г. А. Нуждин, М.В. Ларина

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

Ключевые слова: качество, точность, ротационная вытяжка, эксперимент.

Одним из важнейших требований, предъявляемых к тонкостенным цилиндрическим деталям из малоуглеродистой стали 10, легированных сталей 30ХМА, 10ГН, 12ХЗГНМФБА и алюминиевого сплава АМГ6, является обеспечение заданных формы и размеров. Обеспечение толщины стенки с заданными отклонениями является одним из наиболее важных технических требований при производстве деталей ответственного назначения. Известно, что колебания толщины стенки в готовых деталях зависят в значительной степени от разностенности исходных заготовок, точности используемого оборудования и инструмента. В процессе деформирования исходной заготовки при ротационной вытяжке с утонением стенки происходит изменением величины разностенности, как правило, в сторону уменьшения. Однако в отдельных случаях (при использовании заготовок с малой исходной разностенностью) разностенность в готовой детали может превышать разностенность исходной заготовки, что связано с достижимой

193