Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

SOLUTION OF THREE-DIMENSIONAL DEFORMATION PROBLEMS FOR HETEROGENEOUS AREAS USING THE VARIABLE SEPARATION METHOD BASED ON A VARIATIONAL FORMULATION

Gorshkov A.V. 1 Spevak L.F. 1
1 Institute of Engineering Science RAS
The paper deals with the three-dimensional deformation modeling for heterogeneous bodies. The modeling based on the variational principle of virtual velocities and stresses. It is proposed that the deforming domain consists of sub-domains with different mechanical properties, and that the material of each sub-domain is isotropic. Conjugation conditions at the boundaries between the sub-domains are specified. The calculation of deformation process parameters on the finite time interval is fulfilled using the Lagrangian coordinates. The isochronous variational problem is solved for an arbitrary time moment. As a result the system of ordinary differential equations is obtained for the determining of time-dependent parameters. The system is solved numerically for given initial conditions. As an example, the calculations of the layered perform setting. The displacement representation used allows to simulate the barreling in the process of deformation. The conjugation conditions proposed results in the decreasing the order of differential equation system.
variational principle
deformation
stress
mathematical simulation

Математическое моделирование процессов деформирования материалов (ковки, прокатки, резания и др.) с целью их усовершенствования сопровождается расчетом полей деформаций и напряжений в обрабатываемом изделии. Задача расчета деформирования однородных областей с помощью вариационного принципа скоростей и напряжений исследовалась В.Л. Колмогоровым и его учениками в работах [4, 5, 6]. В [1 ,2] построены варианты вариационного принципа для неизтормических процессов и для случая несимметричного тензора напряжений. В работе [3] приведен пример расчета напряженного состояния цилиндра с каналом.

Рассмотрим задачу о деформировании кусочно-неоднородной заготовки, занимающей объем V и состоящей из K подобластей gorskiv1.wmf. Решение задачи проведем в лагранжевых переменных, так как рассматриваются большие деформации. В каждой подобласти введем локальную сопутствующую систему координат. Будем обозначать лагранжевы локальные переменные трехмерным вектором gorskiv2.wmf с компонентами gorskiv3.wmf. Текущие координаты точек подобласти относительно неподвижной системы отсчета будем обозначать gorskiv4.wmf. Используем суммирование по умолчанию по повторяющимся индексам. Символ в скобках показывает номер подобласти. Суммирование по нему указывается явно. Поэтому номер подобласти располагаем где удобнее – вверху или внизу.

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

gorskiv5.wmf (1)

Здесь gorskiv6.wmf – интенсивность скоростей сдвига, gorskiv7.wmf – интенсивность касательных напряжений, gorskiv8.wmf компоненты девиатора тензора скоростей деформации, gorskiv9.wmf – компоненты девиатора тензора напряжений, gorskiv10.wmf – первый инвариант тензора скоростей деформаций, gorskiv11.wmf – первый инвариант тензора напряжений, m – номер подобласти.

Предполагается, что поле скоростей gorskiv12.wmf и поле напряжений gorskiv13.wmf в пределах каждой подобласти непрерывны.

На внешней поверхности S области V заданы граничные условия:

на части поверхности gorskiv14.wmf заданы скорости gorskiv15.wmf;

на части поверхности Sf заданы внешние усилия gorskiv16.wmf, где gorskiv17.wmf – компоненты вектора нормали к поверхности Sf ;

на части поверхности Ss задана сила трения gorskiv18.wmf, gorskiv19.wmf, gorskiv20.wmf – скорость скольжения, gorskiv21.wmf – единичный вектор в направлении скорости скольжения.

Звездочкой отмечены заданные функции.

Для поверхностей gorskiv22.wmf выполняются соотношения: gorskiv23.wmf и gorskiv24.wmf, gorskiv25.wmf.

В случае неоднородной области на поверхностях контакта подобластей рассматриваются два варианта. Первый: прилипание – скорости непрерывны, для усилий выполняются условия равновесия. Пусть соприкасаются подобласти с номерами p и q. Тогда выполняются соотношения

gorskiv26.wmf,

gorskiv27.wmf. (2)

Здесь gorskiv28.wmf – компонента скорости по нормали к поверхности подобласти, gorskiv29.wmf – касательная компонента скорости, gorskiv30.wmf – компонента усилия на нормаль к поверхности подобласти, gorskiv31.wmf – компонента усилия, касательная к поверхности подобласти gorskiv32.wmf.

Второй – скольжение поверхностей. Поверхность скольжения подобластей gorskiv33.wmf и gorskiv34.wmf обозначим gorskiv35.wmf. В этом случае нормальные компоненты скоростей совпадают, так как тело сплошное, для нормальных компонент усилий выполняются условия равновесия:

gorskiv36.wmf, gorskiv37.wmf. (3)

Разность касательных компонент скоростей равна скорости скольжения

gorskiv38.wmf, (4)

(номер в скобках совпадает с номером поверхности) разность касательных усилий определяется силой трения

gorskiv39.wmf, (5)

где gorskiv40.wmf единичный вектор, направленный параллельно скорости скольжения gorskiv41.wmf.

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

gorskiv42.wmf (6)

где

gorskiv43.wmf

gorskiv44.wmf

gorskiv45.wmf.

Суммирование в gorskiv46.wmf ведется по парам скользящих поверхностей p и q.

Приближенное решение в произвольный момент времени t по предложенному методу следует искать с помощью принципа (6) в виде линейной комбинации функций координат

gorskiv47.wmf (7)

Здесь gorskiv48.wmf, gorskiv49.wmf система линейно-независимых функций лагранжевых координат. Функции выбраны так, что gorskiv50.wmf и gorskiv51.wmf являются виртуальными. Коэффициенты gorskiv52.wmf, gorskiv53.wmf – варьируемые при фиксированном t функции времени; (в правой части (6) по повторяющимся индексам i, j суммирование не производится).

После подстановки (7) в (6) и варьирования, уравнение (6) превратится в систему обыкновенных дифференциальных уравнений относительно gorskiv54.wmf, которую дополняем условиями непрерывности скоростей и равновесия (2) или (3, 4) и условиями трения (5).

Рассмотрим классическую задачу обработки металлов давлением – ковку параллелепипеда плоскими штампами, которая рассматривалась неоднократно в аналогичной постановке [1–4]. Схема задачи показана на рис. 1. Показаны также глобальные оси координат и локальные оси для третьей области.

Для простоты будем считать, что заготовка состоит из прямоугольных слоев размера gorskiv55.wmf постоянной толщины gorskiv56.wmf (gorskiv57.wmf не зависит от y1, y2, но может меняться от слоя к слою). Предполагается, что в процессе обжатия все слои остаются прямоугольными, но могут деформироваться в усеченные пирамиды.

Пусть нижний штамп неподвижен, а верхний перемещается поступательно вниз со скоростью gorskiv58.wmf. Предполагается, что на контакте с инструментом превалирует зона скольжения Ss; зона прилипания SV, примыкающая к центру нижнего и верхнего оснований, практически отсутствует (gorskiv59.wmf). Предположим, что материал параллелепипеда обладает известными реономными свойствами, на контактной поверхности со штампами действует некоторый известный закон трения. Боковая поверхность параллелепипеда – это поверхность типа Sf, на которой gorskiv60.wmf.

gor1.tif

Рис. 1. Схема напряженно-деформированного состояния при осадке многослойного параллелепипеда: 1 – верхний штамп; 2 – заготовка; 3 – нижний штамп

Расчеты проводились для модели материала, для которой функции (1) имеют вид

gorskiv61.wmf, gorskiv62.wmf. (8)

(бингамовский пластик).

Закон трения выбран по Зибелю [1,2]

gorskiv63.wmf. (9)

Очевидно, что закон трения (9) не зависит от скорости скольжения и не имеет обратной функции. В выражениях (8) и (9) gorskiv64.wmf известные величины: предел текучести при чистом сдвиге, коэффициенты «вязкости» и трения соответственно.

Будем предполагать, что физические координаты точек элемента имеют вид

gorskiv65.wmf (10)

Такое представление перемещений позволит учесть возникающую в процессе деформации бочкообразность изделия. Вектор перемещений gorskiv66.wmf связан с вектором gorskiv67.wmf соотношениями

gorskiv68.wmf

После подстановки (7, 8, 10) в (6) и варьирования получим систему относительно <<gorskiv69.wmf>> переменных. Часть переменных связаны и их можно исключить из системы, понизив порядок.

Считаем, что перемешивания материала между слоями заготовки нет, и поэтому уравнение неразрывности выполняется для каждого слоя. В лагранжевых переменных уравнение неразрывности примет вид

gorskiv70.wmf,

gorskiv71.wmf (11)

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

gorskiv72.wmf

Отсюда получаем связь переменных gorskiv73.wmf и gorskiv74.wmf

gorskiv75.wmf(12)

При деформировании заготовки слои не отделяются друг от друга. Поэтому должно выполняться равенство, связывающее переменные gorskiv76.wmf и gorskiv77.wmf

gorskiv78.wmf, gorskiv79.wmf,

и gorskiv81.wmf (13)

Используя соотношения (12, 13) исключим из системы уравнений переменные

gorskiv82.wmf, gorskiv83.wmf, gorskiv84.wmf и gorskiv85.wmf.

Чтобы исключить переменные gorskiv86.wmf введем дополнительные переменные gorskiv87.wmf. Используя условие неразрывности (11) исключим из системы переменные gorskiv88.wmf, gorskiv89.wmf. Переменная gorskiv90.wmf сохраняется. В итоге вектор переменных будет содержать следующие неизвестные функции

gorskiv91.wmf,

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

gorskiv93.wmf

Здесь gorskiv94.wmf – матрица коэффициентов при вторых производных, блочная. Блоки трех-диагональные.

В случае однородной среды типа (8) можно ввести относительное время по формуле

gorskiv95.wmf.

Очевидно, что при такой замене в модели будет только один существенный параметр – относительный коэффициент вязкости m. В случае неоднородной среды относительное время выберем из соотношения

gorskiv96.wmf,

где gorskiv97.wmf – максимальное значение gorskiv98.wmf, а gorskiv99.wmf – минимальное значение gorskiv100.wmf.

gor2.wmf

Рис. 2. Изменение параметров

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

gorskiv101.wmf (нижнее основание неподвижно) gorskiv102.wmf

gorskiv103.wmf

gorskiv104.wmf,

где gorskiv105.wmf – заданная скорость бойка в момент касания.

Относительные значения параметров

gorskiv106.wmf, gorskiv107.wmf, gorskiv108.wmf.

Результаты решения представлены на рис. 2. Используется относительное время.

Такое соотношение параметров можно объяснить следующим образом: gorskiv109.wmf меняется меньше, чем gorskiv110.wmf, так как влияет сила трения на поверхности контакта с инструментом; изменение gorskiv111.wmf существенно меньше, так как сказываются силы инерции.

На рис. 3 тонкой линией показана исходная форма заготовки, жирной – форма после деформации. Пунктиром показано разделение заготовки на слои.

gor3.wmf

Рис 3. Изменение формы заготовки к моменту времени t=0,05 по сравнению с начальной формой

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