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

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

Литература

1. Авернин А.Н., Батыршин И.З. Нечеткие множества в моделях управления и искусственного интеллекта - М.: Мир, 1986.

2. Михайлец В.Б. Оценка объектов с помощью пакета нечеткой логики Fuzzy программы Matlab.

3. Пащенко Д.С., Солдатов В.В., Левиков Д.А. Методы интеллектуального и нечеткого управления проектами автоматизации сельскохозяйственных предприятий

4. Развитие московской канализации (сборник научно-технических статей), ОАО «Можайский полографический комбинат», 2003г., 326 стр.

5. Комиссаров Ю.А., Гордеев Л.С., Нгуен Суан Нгуен Анализ и синтез систем водообеспе-чения химических производств, Москва «Химия» 2002г, 496 стр.

6. Яковлев С.В., Карелин Я.А., Ласков Ю.М., Калицун В.И. Водоотведение и очистка сточных вод, М., Стройиздат 1996г., 592 стр.

7. Дьяконов В.П. Matlab 6/6.1/6.5+SIMULINK 4/5 основы применения, Москва Слон - Пресс 2002, 768стр.

8. Саймон Хайкин Нейронные сети полный курс, Вильямс 2008г., 1104стр.

9. Круглов В.В., Дли М.И., Голунов Р.Ю. Нечеткая логика и искусственные нейронные сети. М: Наука. Физматлит - 2001 г

10. Медведев В.С., Потемкин В.Г., Нейронные сети Matlab 6, Москва, Диалог МИФИ 2002г., 490 стр.

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

д.т.н. проф. Сухомлинов Л.Г., к.т.н. доц. Петров В.К.

МГТУ «МАМИ»

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

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

При постановке задачи на шаге используем следующие декартовы прямоугольные системы координат: глобальную систему координат X,Y,Z при рассмотрении всей дискретной модели в целом и локальные системы координат x,y,z, плоскости x,y которых совпадают с плоскостями треугольных элементов в начале шага нагружения. Матрицы перехода от ло-

кальных координат к глобальным и наоборот формируем традиционным для метода конечных элементов образом.

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

Ч

В каждой точке соответствующей поверхности определяем внешнюю нормаль п (в

проекциях на оси глобальной системы координат) и орты е1, е2, лежащие в плоскости, касательной к поверхности инструмента.

Для вектора перемещений узла дискретной модели оболочки на шаге нагружения используем обозначение и . Связь между проекциями (Ых, у, Ыг) и

(их, Ыу , и2) вектора и на

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

1

— р

8 8 О ху

Обозначим через ^, уу, 2 компоненты тензора конечных деформаций Грина-Лагранжа, вычисляемые относительно конфигурации модели в начале шага нагружения. В

рамках принятого предположения о малости приращения деформаций на шаге параметры ррр

хх, уу, ху для каждого элемента (учитывая, что координаты х, у - прямоугольные) с точностью до пренебрежения этими параметрами по сравнению с единицей можно рассматривать как малые относительные удлинения и сдвиги. На основе известных соотношений нели-

„ 8 8 8 нейной теории упругости для параметров ^, у, х записываем:

8 = V +1 (Х2 + V2 + V2) / ч

хх хх 2^ хх ух гх )

8 = V + V + V V + V V + V V

ху ху ух хх ху ух уу гх гу

где

V =дЫа V

ах ^ а Ят,

дх , ду , (а = х, у, 2).

Производные, входящие в записи (3), выражаем через узловые перемещения с учетом принятого линейного закона распределения перемещений по каждому элементу.

с*

Малое относительное удлинение гг по толщине элемента на шаге нагружения определяем, используя условие несжимаемости

г = -8

гг хх уу

8г = -8х - 8уу

Соотношения (2) - (4) позволяют выразить приращения деформаций ^, у, ху, гг в каждом элементе модели через узловые перемещения на шаге нагружения в проекциях на оси глобальной системы координат Х,У,2.

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

8 8 8 8

щие малые, но конечные компоненты деформации хх, ^, ууу, гг на шаге. В результате

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

хх , , х , гг в элементах дискретной модели оболочки на шаге нагружения

СТ =

2 Я + 2 ст

3 1 + 2Я ^

[(1 + Я + Явуу]

1 Я + 2 ст

ст =---е

ху 3 1 + 2 Я е х

где

ст =

е =

1

1 + 2Я V 3

у 2

4

стхх -ст

) +ст2х +ст;. +:

I

ехх -е

уу

2(1 + 2 Я )Ху ]

)2 + (е - Яе22 )2 + ((егг-е„ )2 + ^

хУ

Здесь Я - параметр нормальной анизотропии материала, а - эквивалентное напряжение, 8 - эквивалентное приращение деформации. Между параметрами а и 8 имеет место зависимость, описываемая соотношениями вида

а = ф

е* = е+е,

где величины

и

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

конце шага нагружения, а ф(е) есть функция упрочнения материала, устанавливаемая экспериментально.

Принятые физические соотношения (5) - (7) с использованием геометрических соотношений (2) - (3), выраженных в терминах узловых перемещений, приводят к явной зависимости между узловыми перемещениями и напряжениями в элементах модели на шаге нагру-жения.

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

Е

Суммирование в левой части равенства (8) осуществляется по всем элементам модели,

а суммирование в правой части - по всем узлам модели. Величины п и Л представляют со-

1*

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

п = п

Величины

р р р 77

х, 7 , 2 представляют собой проекции вектора узловой силы [ на оси глобальной системы координат.

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

о С, О,

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

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

р р р $

жесткого инструмента, проекции х, 7 , 2 соответствующей узловой силы р заранее неизвестны. Поэтому упомянутые уравнения равновесия модели в этом случае необходимо до-

3

2

*

полнить контактными условиями.

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

Р =-№ I

Здесь ц - коэффициент трения,

Ай„

перемещ ение узла относительно поверхности ин-

Р = р п

струмента (в тангенциальном направлении), п

нормальная составляющая силы

Ре Р РпП - тангенциальная составляющая силы р (сила трения).

С учетом того, что

Ре = РеА + Ре,2 е2 Айе = АйеЛ + Айе,2 ё2

Ре,1 = Р  е1 Ре,2 = Р  е2 Айе,1 = Ай е ' е1 Айе,2 = <

1 1 ? ?

векторное равенство (10) переписывается в виде двух равенств в проекциях на заданные тан-

е  е2

генциальные направления 1, 2 в соответствующей точке поверхности инструмента.

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

+ йх, У + йУ, 1 + й1 )= 0

Решение сформулированной задачи о деформации дискретной модели оболочки на шаге нагружения строится (как и в работе [1]) на основе линеаризующих итерационных процессов. При этом линеаризация геометрических соотношений, уравнений равновесия и условий типа (12) выполняется по схеме Ньютона. Линеаризация физических соотношений (5) и условий контактного взаимодействия (10) осуществляется по схеме переменных параметров.

В качестве примера рассмотрим процесс формовки квадратной закрепленной по контуру заготовки из листового алюминия квадратным пуансоном. Этот процесс был экспериментально исследован в работе [3].

Схема формовки показана на рис. 1. Листовая заготовка толщиной Ь0 = 0,5 мм закреплена по контуру размером 80x80 мм (ь х ь). Материал листа изотропный с характеристикой

МПа. Квадратный пуансон имеет размеры 52x52 мм ( а х а ), радиус

упрочнения

ст = 180ё°,285

его скругленных кромок

г = 4

мм. Формовка осуществлялась в условиях смазки, при кото-

рой коэффициент трения был оценен величиной ц = 0,02, а также в отсутствие смазки, когда

Ц = 0,2.

защемление

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

Помимо экспериментальных данных в работе [3] имеются и результаты численных расчетов на основе безмоментной упругопластической инкрементальной конечноэлементной модели.

Раздел 3. Естественные науки. Схема дискретизации четверти заготовки представлена на рис. 2 (аналогичная сетка использована и в работе [3]). Полученная в результате численного моделирования форма оболочки показана на рис. 3.

Рис. 2. Схема дискретизации четверти заготовки

11=16 ми

Рис. 3. Форма оболочки

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

15 Р, [кН]

10

0

^ =0,2

—0— Эксперимент при р=0,02 —— Эксперимент при р-0,2 Результаты расчета на основе модели [3] Результаты расчета на основе предлагаемой модели

=0,2

-Л

у \

^=0,02

... :: ——" =0,02

О

10

15 20

У, [мм]

Рис. 4. Зависимость силы на пуансоне от его перемещения при формовке оболочки

и = 16 мм

Таблица изозначений еи

1 -0,310

2 -0,279

3 -0,284

4 -0,217

5 -0,186

6 -0,155

7 -0,124

8 -0,093

9 -0,062

10 -0,031

Рис. 5. Распределение толщинных деформаций в отформованной оболочке

в случае ц = 0,02

На рис. 5 представлена картина распределения логарифмических толщинных деформаций ь вдоль четверти заготовки при перемещении пуансона и = 16 мм для случая ц = 0,02. Видно, что при формовке квадратным пуансоном максимальное утонение оболочки имеет место в зоне скругленного угла пуансона (в диагональном сечении заготовки). Это полностью согласуется с тем, что наблюдалось в эксперименте (эксперимент показал, что

= -0,3 ёь = -0,2 „ ~

ь - на диагонали и ь ' - на стороне квадратной области; расчет соответственно

sh =-0,31 sh =-0,19, дает h и h ' (см. рис.5)).

На рис. 4 зафиксирован момент потери несущей способности (разрыв) оболочки по падению силовой характеристики процесса (при ц = 0,2) в связи с катастрофическим утонением в зоне угла пуансона.

Выводы

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

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

Литература

1. Сухомлинов Л.Г., Энгельсберг В.К. Конечноэлементная система автоматизированного расчета напряженно-деформированного состояния тонких оболочек в процессах осесим-метричного формоизменения под действием жестких штампов // Известия вузов. Машиностроение, 1989. - N.3. - с.66-71.

2. Сухомлинов Л.Г., Генин Е.В. Численное решение задач о больших пластических деформациях тонких неосесимметричных оболочек под действием заданных нагрузок //Известия вузов. Машиностроение. - 1990. - N.1. - с.16-21.

3. Nakamachi E. Finite element modelling of the punch press forming of thin elastic plates // Proceedings of the NUMMIFORM'86 Conference. - Gothenburg. - 25-29 August. - 1986. -pp.333338.