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

Технология и технологические

машины

ЛАРИН Сергей Николаевич

кандидат технических наук, доцент

УДК 621.983; 539.374

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

из анизотропного материала

— 1 в режиме кратковременной ползучести

С.Н. Ларин

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

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

The article presents the results of theoretical investigations of stressed and defomed states, power circumstances, geometric sizes of an original piece and extreme defamation levels in the process of pneumatic faming of hemispherical parts from anisotropic material under a short-duration creeping mode. The influence of mechanical properties anisotropy, loading conditions and geometric sizes of an original piece on the investigated parameters of this process has been established.

Keywords: anisotropy, mathematical model, high strength materials, hemispherical parts, pneumatic forming, creeping, damageability, failure.

1 Работа выполнена по ведомственной целевой программе «Развитие научного потенциала высшей школы (2009—2011 гг.)», грантам РФФИ, гранту Президента РФ для поддержки молодых ученых — кандидатов наук и по государственным контрактам в рамках Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009—2013 гг.

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

Изотермическое формоизменение куполообразных деталей газом из листовых высокопрочных алюминиевых и титановых сплавов имеет значительные преимущества перед традиционными методами обработки и весьма перспективно при использовании его в промышленности [1—4].

Основные предположения

Рассмотрим деформирование круглой листовой заготовки радиусом R0 и толщиной h0 свободным выпучиванием в режиме вязкого течения материала под действием избыточного давления газа p = p0 + apt"p в сферическую матрицу (рис. 1). Здесь p0, ap ,np — константы на-гружения.

По внешнему контуру заготовка закреплена. Материал заготовки принимается трансвер-сально-изотропным с коэффициентом анизотропии R; напряженное состояние оболочки — плоским, т. е. напряжение, перпендикулярное плоскости листа, равно нулю (az = 0). Рассматривается деформирование в меридиональной плоскости оболочки, как мембраны. В силу симметрии механических свойств материала относительно оси заготовки и характера действия внешних сил меридиональные, окружные и нормальные к срединной поверхности заготовки напряжения и скорости деформаций являются главными.

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

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

Рис. 1. Схема к расчету деформированного состояния срединной поверхности заготовки в меридиональной плоскости

поверхности и сечения оболочки конической поверхностью р(, перпендикулярной дуге меридиана, определяют по формуле

Pm = Pt = P =

н2 + r02

2H

где Н — высота купола в данный момент времени деформирования.

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

£c =

Stc

2HH H2 + R,

-; £c =

mc

2HH

h 2 + r2 ' £ cc .

h h'

10 1 -"0

Здесь Н = с1Н / йГ, Иг = й/ / Л.

По контуру заготовка закреплена (точка к), т. е. скорость деформации вдоль контура равна нулю, £,Ск = 0, и в соответствии с ассоциированным законом течения

£ л - 0; а ,к -

Яа

тк

1 + Я

' £ гк £

тк

где Я — коэффициент нормальной анизотропии при вязком течении материала.

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

£с -

т

0

9 а

С£а]а; £ С -/г

С08 9

81п а

£с -

г

Здесь 9 — текущий угол между вертикальной осью симметрии заготовки и радиус-вектором, определяющим положение точки в сечении срединной поверхности диагональной плоскостью; а - (1а / &.

При деформации оболочки принималось, что на каждом этапе деформирования имеет место радиальное течение точки срединной поверхности в меридиональной плоскости относительно нового центра в момент , + &, т. е. в направлении 9 + (9.

Связь между углом а и временем деформирования ,, когда задана функциональная связь Н - Н(,), устанавливается следующим образом:

Н

а - 2агс1£ -

Я

Толщина оболочки в куполе срединной поверхности оболочки (9 - 0) определяют по выражению

/г - Нп

1+

Н2

ЯI

Изменения толщины оболочки от времени деформирования , в месте ее закрепления (9 - а) оценивают по формуле

Н

/г - кп

Я

1+

Н2

\

Я -

агс1§

Н_

Я

Напряженное и деформированное состояние оболочки

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

ат , аг Р. РР,

Рт

+ р, г; а т* 2г

Решая их совместно, с учетом того, что Рт - Р,, находим

РР

а т - а ' - •

Эквивалентные скорость деформации £се

и напряжение а е в вершине купола (точка с) и в точке закрепления оболочки по контуру (точка к) для анизотропного материала вычисляют соответственно по выражениям:

л/3

а „„ -

2л/2 + Я

1/2

£с -

£ ек

' ек

'2 (2 + Я)(Я + 1)' _

3 2Я + 1 | ^тк;

1/2

а

2Я + 1

2

тк

Силовые режимы и предельные возможности формоизменения

Рассмотрим медленное изотермическое деформирование оболочки из материала, для которого справедливы уравнения состояния энергетической теории ползучести и повреждаемости [2]:

£ С =

B(aе / a ,o)П

(i-»a)m

Ancp

где B, n, m — константы материала, зависящие от температуры испытаний; »А — повреждаемость материала при вязкой деформации по энергетической модели разрушения; Ар —

удельная работа разрушения при вязком течении материала; ш А = d»A / dt; £ си a е — эквивалентные скорость деформации и напряжение соответственно; a — эквивалентное напряжение, разделяющее вязкое и вязкопластическое течение материала.

Величину удельной работы разрушения А^р

при вязком течении анизотропного материала определяют по выражению

Апср = D(b0 + b1 cos a + b2 cos p + b3 cos y),

где D, b0, b1, b2, b3 — константы материала; a — среднее напряжение, a = (a 1 + a2 + a3)/3; a 1 ,a2, a3 — главные напряжения; a, в, у — углы ориентации первой главной оси напряжений a1 относительно главных осей анизотропии x , y и z соответственно.

Поскольку давление p в каждый момент деформирования равномерно распределено по поверхности оболочки, то будем находить его величину в вершине купола оболочки (точка с).

Подставив в первое из уравнений (12) входящие в него величины a е и £ С, определяемые по формулам (10), с учетом соотношений (1), (4) и (9), получим

pndt =

П+1

a*0 (1-шAc)m 22n+2(2 + R)T Hn+1hndH

n+1

~2

3 2 B

n+1

Толщина оболочки И определяется по выражению (6).

Найдем величину накопления повреждаемости . Подставив во второе уравнение состояния (12) выражения (10) с учетом (1), (4) и (9), получим

p

1+

H

2 V

R

0 у

h0 Anp

н.

Это уравнение удобно использовать, если на-гружение такое, что p = const.

Если подставить первое уравнение состояния во второе, то получим другую форму уравнения для нахождения повреждаемости:

' е 0

n+1

a; b 1/n

Уравнение (15) удобно использовать при £cc = £c1 = const. В последнем случае интегрирование уравнения (15) приводит к выражению вида

»Ac =1-

n — m

n+1

АПр Bl/n

n

n—m

.

Время разрушения t* определяется из усло-

вия »Ac = 1:

t=

АПр B 1/nn

n

a е0

n+1

Давление р, необходимое для реализации условий деформирования, вычисляют по формуле

p(t ) =

aе0(1 — »Ac )m/n 22(2 + R)1/2Hhfr z 1/я

л/3£1/и

Зависимость = ^) находят согласно соотношению (16), а Н = Н(/) может быть определена из уравнения

Н2 + я2

ln:

t=

242 + R H 02 + R02

s

£ c1

Предельную высоту купола Н* найдем по уравнению (19) при t = I*.

Задание функции Н = Н({) позволяет определить = (t) из выражений (15) или (16), а функцию р = р^) вычисляют по формуле (13).

2

c

»Ac =

c

»Ac =

n

Аналогичным образом выполнены исследования напряженного и деформированного состояний заготовки в точке закрепления оболочки (точка к), также получены основные уравнения и соотношения для решения поставленной задачи в предположении, что поведение материала подчиняется уравнениям кинетической теории ползучести и повреждаемости при известном законе давления от времени Р- Р (,)и при постоянной эквивалентной скорости деформации в куполе заготовки £е1.

Разработаны алгоритм расчета силовых и деформационных параметров исследуемого технологического процесса и программное обеспечение для ЭВМ.

Анализ результатов

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

Расчеты выполнены для титанового сплава ВТ6С при температуре Т - 860 ° С, поведение которого описывается энергетической теорией ползучести и повреждаемости, и для титанового сплава ВТ14 при температуре Т - 950 ° С, поведение которого подчиняется кинетической

теории ползучести и повреждаемости. Механические характеристики этих материалов при формоизменении в условиях вязкого течения материала приведены в работе [3].

Зависимости изменения величин давления газа Р, относительных величин толщины заготовки в куполе Нс - Нс / /г0 ив месте ее закрепления кк - кк / к0, высоты куполообразной заготовки Н' - Н / Я0 от времени деформирования , для титанового сплава ВТ6 (Т - 860 °С) при постоянной величине эквивалентной скорости деформации в куполе заготовки £е1 представлены на рис. 2. Экспериментальные данные на рисунке обозначены точками.

Из анализа результатов расчетов и графических зависимостей следует, что с ростом времени деформирования , до определенного предела резкое увеличивается относительная высота заготовки Н' и уменьшается относительная толщина заготовки в куполе Нс ив месте ее закрепления кк. Дальнейшее увеличение времени деформирования , приводит к плавному изменению исследуемых величин. В момент времени ,, близком к разрушению заготовки, происходит резкое изменение относительных величин Н', Нс и Нк. Это связано с интенсивным ростом

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

Установлено, что изменение относительной толщины в куполе заготовки Нс происходит бо-

Рис. 2. Зависимости изменения Р, Н' и к в рассматриваемых точках заготовки от времени деформиролвания

, для титанового сплава ВТ6С

лее интенсивно по сравнению с изменением относительной толщины в месте ее закрепления Ик. С ростом времени деформирования t эта разница увеличивается и может достигать 50%.

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

Сопоставление теоретических и экспериментальных данных по относительной толщине в куполе заготовки Нс и месте ее закрепления Ик, а также относительной высоте заготовки Н'

указывает на удовлетворительное их согласование (до 10%).

Установлено, что разрушение заготовки при изотермическом деформировании происходит в куполе детали, где наблюдается максимальное утонение заготовки.

Зависимости изменения времени разрушения t*, относительной высоты Н' и толщины в куполе заготовки И* в момент разрушения, определенных по величине накопленных микроповреждений при ш^ = 1, от величины постоянной эквивалентной скорости деформации

в куполе заготовки ^е1 и коэффициента анизотропии Я представлены на рис. 3 и 4 соответственно. Показано, что увеличение параметров закона нагружения ар, пр и величины постоянной эквивалентной скорости деформации в куполе заготовки ^е1 приводит к уменьшению времени разрушения и и относительной высоты заготовки Н'*, а также к увеличению относительной толщины в куполе заготовки И*.

Из анализа рис. 4 следует, что коэффициент нормальной анизотропии Я оказывает существенное влияние на величину времени разрушения и относительные величины Н'*, И*. С ростом коэффициента анизотропии Я относительное значение И* резко увеличивается, а время разрушения t* и относительная высота заготовки Н'* резко уменьшаются. Установлено, что не учет анизотропии механических свойств заготовки при анализе процесса изотермического формоизменения сферической оболочки дает погрешность в оценки времени разрушения t* порядка 35%, а относительной высоты Н'* и толщины в куполе заготовки И* в момент разрушения — 15%.

Анализ результатов расчетов показывает, что предельные возможности формоизменения в режиме вязкого течения материала, поведение которого подчиняется кинетической теории ползучести и повреждаемости (сплав ВТ14), не зависят от условий нагружения заготовки.

Рис. 3. Зависимости изменения t, и Н, И, в куполе заготовки от ^е для титанового сплава ВТ6

Рис. 4. Зависимости изменения и Н',, к, от К (кинетическая теория; К0 - 300; Р0 - 0,013 МПа ;

а -4-10"3 МПа/ е"Р; пр - 0,6)

Установлена существенная зависимость времени разрушения от параметров нагружения ар, пр и величины постоянной эквивалентной скорости деформации £е1. Увеличение параметра нагружения ар от значения 0,2-10"3 МПа / с"Р

до 1,4 • 10-3 МПа / с"Р и пр от значения 0,5 до 0,8

при фиксированных других параметрах приводит к уменьшению времени разрушения и в 1,8 раза.

Показано существенное влияние геометрических размеров заготовки на величину времени разрушения . Установлено, что увеличение относительной величины радиуса заготовки К0 - К0 / к0 от 200 до 800 приводит к уменьшению времени разрушения более чем в 4 раза.

Литература

1. Романов К.И. Механика горячего формоизменения металлов. М.: Машиностроение, 1993. 240 с.

2. Яковлев С.П., Яковлев С.С., Андрейченко В.А. Обработка давлением анизотропных материалов. Кишинев: Квант, 1997. 332 с.

3. Изотермическое деформирование высокопрочных анизотропных материалов / С.С. Яковлев и др. М.: Машиностроение, 2004. 427с.

4. Изотермическая пневмоформовка анизотропных высокопрочных листовых материалов / С.С. Яковлев, С.Н. Ларин и др. М.: Машиностроение, 2009. 352 с.

Статья поступила в редакцию 10.09.2011 г