Xreferat.com » Рефераты по промышленности и производству » Математическая модель процесса вытяжки трубчатой заготовки

Математическая модель процесса вытяжки трубчатой заготовки

1. Вариационные подходы к решению задач методом конечного элемента


Основная идея МКЭ основывается на замене некоторой непрерывной величины в пределах рассматриваемой области дискретной моделью, которая строится на множестве кусочно-непрерывных функций, определённых на конечном числе подобластей, называемых конечными элементами (КЭ). Неизвестная искомая величина в пределах каждого КЭ аппроксимируется, как правило, полиномиальной функцией заданного вида с учётом требования непрерывности на границах смежных КЭ. При этом выбор формы конечного элемента и вида выражения, аппроксимирующего действительный закон изменения исследуемой величины в пределах КЭ, является одним из наиболее ответственных моментов в общей процедуре МКЭ, от которого существенно зависит точность приближённого решения. Таким образом, непрерывная в пределах исследуемой области неизвестная величина (например, перемещение, скорость перемещения, напряжение, температура и т. д.) представляется через конечное число её дискретных значений в узлах элементов.

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


Математическая модель процесса вытяжки трубчатой заготовки,(7)


где AД - работа или мощность внутренних сил; AМ - работа или мощность, развиваемая массовыми силами; AВ - работа или мощность внешних сил.

Дальнейшая процедура МКЭ предусматривает представление выражения (2.1) в виде функционала значений, неизвестных только в узлах КЭ, и построение разрешающей системы уравнений путем минимизации J по всем узловым переменным:


Математическая модель процесса вытяжки трубчатой заготовки(7)


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

Способ получения разрешающих уравнений для КЭ, основанный на оптимизации функционала (1), является общепризнанным при теоретическом решении задач ОМД, поскольку вариационные принципы имеют наглядный физический смысл и достаточно строгое математическое обоснование.

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

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

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

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

Каждому из перечисленных вариационных принципов соответствует определённая форма МКЭ. Принципу минимума полной мощности (полной энергии) соответствует кинематический метод, принципу минимума дополнительной работы - метод напряжений, а вариационному принципу Рейсснера - смешанный метод.

При нагружении тела потенциальная энергия внешних сил изменяется. При этом внешние силы совершают работу. Потенциал внешних сил Q на возможных перемещениях δu численно равен работе этих сил:


Математическая модель процесса вытяжки трубчатой заготовки(7)


где Pi – поверхностные силы, S – площадь поверхности тела.

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


Математическая модель процесса вытяжки трубчатой заготовки(7)

где sij - компоненты тензора напряжения, eij - компоненты тензора деформации, V – объём тела.

Сумма энергии деформации и потенциала внешних сил равна полной потенциальной энергии:


Математическая модель процесса вытяжки трубчатой заготовки(7)


В соответствии с принципом возможных перемещений Лагранжа изменение полной потенциальной энергии на возможных перемещениях равняется нулю:


Математическая модель процесса вытяжки трубчатой заготовки(7)


При этом под возможными перемещениями du понимаются сколь угодно малые отклонения системы от положения равновесия, допускаемые наложёнными на систему связями. Из уравнения (6) следует, что в состоянии равновесия энергия П имеет стационарное значение. Можно показать, что в положении устойчивого равновесия этот экстремум соответствует минимуму.

С учётом изложённого вариационное уравнение Лагранжа для статической задачи имеет вид:


Математическая модель процесса вытяжки трубчатой заготовки(7)


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

2. Основные соотношения метода конечных элементов


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


Математическая модель процесса вытяжки трубчатой заготовки

Рис. 1. Конечный элемент в задаче осесимметричной деформации


Вектор перемещений узловых точек конечного элемента в случае осесимметричной деформации имеет вид:


Математическая модель процесса вытяжки трубчатой заготовки.(8)


Произвольная точка элемента получает перемещения ur и uz в направлении осей r и z. Поэтому матрица u имеет вид:

Математическая модель процесса вытяжки трубчатой заготовки.(9)


Узловые перемещения Математическая модель процесса вытяжки трубчатой заготовки и u связаны между собой матрицей аппроксимирующих функций N:


Математическая модель процесса вытяжки трубчатой заготовки(9’)


Наиболее распространен способ получения приближённых решений на основе использования вариационного уравнения по методу Релея - Ритца. Он заключается в том, что функции перемещений задаются в виде интерполяционного полинома. Если ограничиться полиномом первой степени, то эти функции будут иметь вид:


Математическая модель процесса вытяжки трубчатой заготовки(10)


Здесь ai - произвольные постоянные. При линейной аппроксимации стороны треугольника после деформирования элемента остаются прямыми.

Выразим ai через перемещения узлов элемента. В результате матрица N примет вид:


Математическая модель процесса вытяжки трубчатой заготовки(11)


S - площадь сечения элемента:

Математическая модель процесса вытяжки трубчатой заготовки,(12)


где ri, zi - координаты i-го узла в соответствующих осях.

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


Математическая модель процесса вытяжки трубчатой заготовки(13)


В условиях осесимметричной задачи тензор деформации второго ранга сводится к вектору:


Математическая модель процесса вытяжки трубчатой заготовки(14)


компоненты которого выражаются через производные перемещений по соответствующим координатам:


Математическая модель процесса вытяжки трубчатой заготовки.(15)


Связь между составляющими векторов деформации и перемещений можно представить одним матричным равенством:


Математическая модель процесса вытяжки трубчатой заготовки(16)


где B – матричный дифференциальный оператор:


Математическая модель процесса вытяжки трубчатой заготовки.(17)


Используя (16) и (17), можно выразить деформации через узловые перемещения

Математическая модель процесса вытяжки трубчатой заготовки.(18)

Матрица функций формы C для осесимметричной деформации:


Математическая модель процесса вытяжки трубчатой заготовки.(19)


Коэффициенты матрицы C зависят от координат r и z точки внутри элемента. Для треугольника с узлами в вершинах координаты r и z можно заменить средними по элементу значениями:


Математическая модель процесса вытяжки трубчатой заготовки(20)


Вектор напряжений s имеет вид:


Математическая модель процесса вытяжки трубчатой заготовки(21)

Выразим с помощью линейного закона, выражаемого матрицей жёсткости, напряжения через узловые перемещения


Математическая модель процесса вытяжки трубчатой заготовки,(21’)


где D – матрица материальных констант.

Потенциальная энергия деформации элемента с учётом (20) и (19)


Математическая модель процесса вытяжки трубчатой заготовки.(22)


Интеграл в выражении (2.22) есть матрица жёсткости выбранного элемента


Математическая модель процесса вытяжки трубчатой заготовки,(23)


Элементарный объём Математическая модель процесса вытяжки трубчатой заготовки. Поэтому матрица жёсткости элемента записывается следующим образом:


Математическая модель процесса вытяжки трубчатой заготовки,(24)


где S – площадь элемента.

С учётом проделанных преобразований уравнение равновесия элемента через узловые перемещения выражается в форме:


Математическая модель процесса вытяжки трубчатой заготовки(25)


где K - матрица жёсткости; P, Математическая модель процесса вытяжки трубчатой заготовки - векторы внешних сил и узловых перемещений, соответственно.

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


Математическая модель процесса вытяжки трубчатой заготовки(26)


Одним из способов решения задачи в приращениях является метод последовательных нагружений. Для квазистатической задачи приращения внешних сил DP вычисляются на шаге по времени Dt. При этом вектор внешних сил P в момент времени t равен:


Математическая модель процесса вытяжки трубчатой заготовки(27)


где n – шаг нагружения.

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


Математическая модель процесса вытяжки трубчатой заготовки(28)


где Математическая модель процесса вытяжки трубчатой заготовки- вектор приращений перемещений.


3. Представление матрицы жёсткости


В пределах упругости связь между приращениями напряжений и деформации выражается законом Гука. Согласно ему компоненты приращений деформации являются линейными функциями приращений напряжений. Пластическое состояние материала описывается теорией малых упругопластических деформации Ильюшина. Принимается теория изотропного упрочнения. Объёмная деформация в пластической зоне остается упругой и для нее выполняется объёмный закон Гука:


Математическая модель процесса вытяжки трубчатой заготовки,(29)


где q - относительное изменение объёма.

Модуль объёмного сжатия k для изотропного тела в случае осесимметричной деформации имеет вид:


Математическая модель процесса вытяжки трубчатой заготовки.(30)


Модуль сдвига G связан с модулем Юнга E и коэффициентом Пуассона n формулой:

в упругой области:


Математическая модель процесса вытяжки трубчатой заготовки(31)

в пластической:


Математическая модель процесса вытяжки трубчатой заготовки(32)

Здесь H – касательный модуль упрочнения. Коэффициент Ляме - l определяется формулой:


Математическая модель процесса вытяжки трубчатой заготовки(33)

Таким образом, матрица материальных констант D имеет вид:


Математическая модель процесса вытяжки трубчатой заготовки.(34)


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

Зная текущее состояние элемента, предел текучести, накопленную деформацию и приращения внешних сил, можно определить изменение напряжённо-деформированного состояния на шаге приращения перемещений Du и сил DР, используя для вычисления K по формуле Error: Reference source not found упругое или пластическое представление матрицы жёсткости.


4. Пластическая деформация


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

тело изотропно;

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

полные приращения составляющих деформации Deij складываются из приращений составляющих упругой деформации Deeij и пластической деформации Depij:


Математическая модель процесса вытяжки трубчатой заготовки;


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

Напряжённо-деформированное состояние элемента на i+1 шаге характеризуется интенсивностью деформации ei:


Математическая модель процесса вытяжки трубчатой заготовки(35)


где eij - компоненты тензора деформации.

Если интенсивность деформации какого - либо конечного элемента превысила текущий предел упругости по деформациям Математическая модель процесса вытяжки трубчатой заготовки, то этот элемент переходит из упругого в пластическое состояние. Если материал упрочняется при пластическом деформировании, то соответствующая пределу упругости деформация εе увеличивается на величину Deе (Рис. 7):


Математическая модель процесса вытяжки трубчатой заготовки(36)


Вычисление предела упругости по деформациям Математическая модель процесса вытяжки трубчатой заготовки, достигнутого на шаге k определяется суммированием:

Математическая модель процесса вытяжки трубчатой заготовки.(37)


Имеется в виду, что в упругой области предел упругости не изменяется, его приращения не вычисляются и равны нулю.

Накопленная пластическая деформация определяется разностью интенсивностей полной деформации ei и деформации ee, соответствующей пределу упругости:


Математическая модель процесса вытяжки трубчатой заготовки(38)


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


Математическая модель процесса вытяжки трубчатой заготовки,(39)


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

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


Математическая модель процесса вытяжки трубчатой заготовки,(40)

где eеp – деформация в области нелинейной упругости:


Математическая модель процесса вытяжки трубчатой заготовки.


Вектор приращений компонент тензора напряжения на шаге k в пластическом состоянии определяется по приращениям компонент деформации:


Математическая модель процесса вытяжки трубчатой заготовки.(41)


Вектор компонент напряжения на шаге k в упругом и пластическом состоянии суммируется по приращениям:


Математическая модель процесса вытяжки трубчатой заготовки.(42)


Интенсивность напряжений определяется по компонентам тензора напряжения sij:


Математическая модель процесса вытяжки трубчатой заготовки.(43)


Математическая модель процесса вытяжки трубчатой заготовки

Рис. 2. Изменение предела упругости по деформациям при упрочнении


Если интенсивность деформации уменьшилась:


Математическая модель процесса вытяжки трубчатой заготовки,(44)


то материал разгружается и переходит в упругое состояние. При нарушении неравенства (2.44) вновь происходит переход элемента в пластическое состояние.


5. Оценка повреждаемости заготовок


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

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

Если влиянием истории деформирования пренебречь, то можно использовать критерий Смирнова-Аляева:


Математическая модель процесса вытяжки трубчатой заготовки(45)


Либо, нормируя на единицу, получим меру повреждений y:

Математическая модель процесса вытяжки трубчатой заготовки(46)


где ep(h) - предельная деформация в момент появления первых трещин, обнаруживаемых визуально; h - показатель напряжённого состояния:


Математическая модель процесса вытяжки трубчатой заготовки(47)


s - среднее нормальное напряжение; si – интенсивность напряжений.

Для учёта влияния истории деформирования и использования соотношения (47) для простого нагружения, примем за меру повреждений y выражение (критерий Колмогорова):


Математическая модель процесса вытяжки трубчатой заготовки,(48)


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

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


6. Взаимодействие заготовки с инструментом


Заготовка представляет собой совокупность узлов, связанных между собой конечно-элементной сеткой. Далее “узлами” будут называться узлы конечно–элементной сетки заготовки. Поведение узлов описывается соотношениями МКЭ. Инструмент является абсолютно жёстким телом, ограниченным непроницаемыми для узлов заготовки границами. Границы инструмента аппроксимируются прямолинейными и радиусными участками, состыкованными друг с другом в единую поверхность. В дальнейшем под именем “граница” будет пониматься граница инструмента. Инструмент может быть неподвижным и подвижным. Для адекватного описания технологических процессов штамповочного производства математическая модель позволяет использовать один подвижный инструмент и несколько неподвижных (например, матрица и оправка).

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

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


Математическая модель процесса вытяжки трубчатой заготовки,(49)


где r – радиальная степень свободы; z – осевая; k – тангенс угла наклона границы.

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

Математическая модель процесса вытяжки трубчатой заготовки(50)


То есть в ячейки, соответствующие степеням свободы r и z данного узла, вносятся стоящие при них в уравнении (50) коэффициенты. Таким образом, удовлетворяется уравнение связи (49).


Математическая модель процесса вытяжки трубчатой заготовки

Рис. 3. Ограничения неподвижных границ


Подвижная граница перемещается в заданном направлении (вертикальном или горизонтальном, в зависимости от вида процесса) на величину DH на каждом шаге:


Математическая модель процесса вытяжки трубчатой заготовки(51)


где H – полный ход инструмента; i – число шагов решения.

Если после очередного перемещения границы какой-либо узел (или несколько) может оказаться “в теле” инструмента, то ему назначается принудительное перемещение на величину Dh:


Математическая модель процесса вытяжки трубчатой заготовкиилиМатематическая модель процесса вытяжки трубчатой заготовки(52)

Здесь R, Z – координаты подвижной границы; r, z – координаты узла.

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


Математическая модель процесса вытяжки трубчатой заготовкиилиМатематическая модель процесса вытяжки трубчатой заготовки(53)


Для системы уравнений это означает внесение в соответствующие данному узлу ячейки матрицы жёсткости коэффициентов 0 и 1 в зависимости от направления движения и в правую часть величины Dh.


Математическая модель процесса вытяжки трубчатой заготовки

Рис. 4 Ограничения подвижных границ


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


Математическая модель процесса вытяжки трубчатой заготовки


Рис. 5. Процедура уточнения положения узла при смещении его подвижным инструментом


7. Трение


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

Модуль напряжения трения определяется законом Кулона (54):


Математическая модель процесса вытяжки трубчатой заготовки,(54)


где sN – нормальное напряжение на границе инструмента; m – коэффициент трения скольжения; tK – касательное к границе напряжение .

Математическая модель процесса вытяжки трубчатой заготовки

Рис. 6. Контакт конечного элемента заготовки с инструментом


При пластическом течении касательное напряжение в элементах tK, контактирующих с внешними телами (матрица или пуансон), не должно превышать предел текучести при сдвиге tS:


Математическая модель процесса вытяжки трубчатой заготовки(55)


Для выполнения этого условия предложен следующий алгоритм расчёта. Для элемента, лежащего на границе, определяется касательное напряжение tK, по которому вычисляется узловая сила трения:


Математическая модель процесса вытяжки трубчатой заготовки(56)


где SK – площадь контакта элемента с границей инструмента. Если к узлу примыкают два элемента, лежащие на границе, то узловая сила трения получается суммированием сил, вычисленных по формуле (56). Для определения направления силы трения реализуется следующий алгоритм. Сначала выполняется шаг нагружения без учёта трения (FTP=0). Затем по результатам этого шага в каждом узле на границе контакта определяется нормальная узловая сила на границе контакта, приращение перемещения вдоль неё и, описанным

Если Вам нужна помощь с академической работой (курсовая, контрольная, диплом, реферат и т.д.), обратитесь к нашим специалистам. Более 90000 специалистов готовы Вам помочь.
Бесплатные корректировки и доработки. Бесплатная оценка стоимости работы.

Поможем написать работу на аналогичную тему

Получить выполненную работу или консультацию специалиста по вашему учебному проекту

Похожие рефераты: