Full Text

Введение В работе построена методика численного исследования изотропных материалов с использованием левого тензора Коши–Грина, для которых физические соотношения задаются с помощью упругого потенциала. В первой части рассматривается кинематика среды, которая описывается с помощью левого тензора Коши–Грина, тензора пространственного градиента скорости и тензора деформации скорости. Далее получено разрешающее уравнение в скоростях напряжений Коши–Эйлера. Во второй части построены определяющие соотношения и уравнения пластического течения с использованием уравнения предельного состояния, ассоциированного закона течения и потенциала упругой деформации, из которого получены линеаризованные физические соотношения для обобщенной производной Трузделла напряжений. В третьей части описан метод проецирования напряжений на поверхность текучести. В четвертой – алгоритм решения нелинейной задачи. Используется метод последовательных нагружений и процедура проецирования напряжений на поверхность текучести с итерационным уточнением напряженно-деформированного состояния (НДС). В пятой части рассмотрен пример построения определяющих соотношений для материала II порядка. В последней части приводятся результаты решения задач. Численная реализация основана на методе конечных элементов. 1. Кинематика среды. Уравнение в скоростях напряжений В качестве тензоров, описывающих деформацию и скорость деформации, используются тензор градиента деформации , мера деформации Фингера , тензор пространственного градиента скорости , тензор деформации скорости , где – компоненты вектора скорости ; – компоненты радиуса-вектора; – диадное произведение орт декартовой системы координат. Напряженное состояние описывается с помощью тензора истинных напряжений , определенного в актуальном состоянии. Для решения задач с учетом физической нелинейности (в первую очередь задач пластического деформирования) получили распространение формулировки разрешающих вариационных уравнений в скоростях деформаций и напряжений. Такие уравнения могут быть получены дифференцированием по времени уравнения принципа виртуальных мощностей в актуальной конфигурации , (1) где – текущий объем; – часть поверхности, на которой заданы усилия; , – векторы объемных и поверхностных сил соответственно. После линеаризации (1) получим уравнение в скоростях напряжений (2) где – относительное изменение объема. 2. Определяющие соотношения Определяющие соотношения получим, используя потенциальную энергию упругой деформации, которую определяет функция W. В качестве аргументов функции W примем компоненты тензора меры деформации Фингера, т.е. (3) тогда тензор напряжений Коши–Эйлера будет выражаться в следующем виде: (4) Для изотропного материала, свойства которого не зависят от направления, функция удельной потенциальной энергии деформации должна зависеть лишь от инвариантов соответствующих тензоров: (5) тогда тензор напряжений можно выразить как (6) Линеаризуя соотношение (4), получим выражение для скорости напряжений Коши–Эйлера (7) и (8) где введено обозначение (9) В результате получаем физическое соотношение упругого деформирования для производной Трусделла в виде линейного уравнения (10) где – производная Трусделла тензора напряжений . Моделирование упругопластических деформаций основано на аддитивном представлении полной деформации скорости, т.е. , где – упругая составляющая, а – пластическая. Предполагается справедливость ассоциированного закона течения (11) где ‑ скорость пластических деформаций; – функция текучести. Рассмотрим в качестве критерия упругого деформирования условие Губера–Мизеса, которое для изотропной среды допускает обобщение в виде (12) где – интенсивность напряжений; – функция упрочнения; – параметр упрочнения. Тогда, используя (12), пластическую деформацию скорости можно записать следующим образом: . (13) 3. Метод проецирования напряжений на поверхность текучести По известным параметрам k-го состояния определим (k + 1)-е по следующей формуле: где – приращение параметра (времени), определяющее переход от предыдущего состояния к последующему. Запишем последнее соотношение в следующем виде: (14) Здесь (15) – так называемый тензор «пробных» напряжений. Уравнение (14) определяет снос напряжений на поверхность текучести. 4. Общий алгоритм решения Процесс деформирования представляется в виде последовательности равновесных состояний. Переход от предыдущего состояния к последующему происходит путем приращения нагрузки. Разрешающее уравнение на k-м шаге нагружения будет иметь вид (16) Так как решаются задачи квазистатические, то можно перейти от скоростей к приращениям, например (и принять ). Решая уравнение (16), получим вектор перемещений , который определяет конфигурацию на следующем шаге нагружения: , (17) и напряженное состояние . (18) Далее для учета пластических деформаций применяем метод проецирования (14). В результате использования метода «проецирование напряжения на поверхность текучести» полученное напряженное состояние не удовлетворяет разрешающей системе уравнений. Поэтому воспользуемся итерационным уточнением НДС. Эта итерационная процедура основана на введении в разрешающее уравнение вариации мощности «дополнительных напряжений» на возможных деформациях скорости, где дополнительные напряжения определяются как разность истинных и пробных напряжений, определяемых по формуле (15) Итоговое уравнение для m-й итерации на k-м шаге нагружения в матричном виде , (19) где – вектор дополнительных напряжений. 5. Физические соотношения и разрешающее уравнение для материала II порядка Рассмотрим пример построения физических соотношений для следующего потенциала упругих деформаций: (20) где – параметры Ляме. Получим выражение для скорости изменения напряжений Коши–Эйлера (21) где (22) Подставив физические соотношения (22) в уравнение в скоростях напряжений (2), получим разрешающее уравнение (23) 6. Численные примеры В качестве базового в настоящей работе используется произвольный восьмиузловой конечный элемент. При вычислении интегралов используется схема численного интегрирования. Окончательно получается система линейных алгебраических уравнений на k-м шаге нагружения: (24) 6.1. Упругопластическое растяжение круглого стержня Рассмотрим задачу растяжения круглого стержня со следующими параметрами: R = 6,413 мм, R0 = 0,982R мм, L = 26,667 мм. Отметим, что для конкретизации места образования шейки в центре стержня задается снижение радиуса на 1,8 %. Критерием пластичности служит условие Губера–Мизеса (12). Функция нелинейного изотропного упрочнения имеет вид: . Параметры материала: МПа, , МПа, МПа, , . Рис. 1. Конечно-элементная дискретизация половины стержня Конечно-элементная дискретизация половины стержня представлена на рис. 1. На торце задается перемещение. Так как задача обладает тремя плоскостями симметрии, то рассматривалась восьмая часть стержня с соответствующими граничными условиями. На рис. 2 показана интенсивность пластических деформаций для конечного положения и диаграмма «сила – перемещение». а б Рис. 2. Интенсивность пластических деформаций (а); диаграмма «сила – перемещение торца» (б); сплошная кривая – решение по описанной методике; ♦ – решение; ○ – решение 6.2. Упругопластическое деформирование прямоугольной плиты Исследуем напряженно-деформированное состояние плиты из материала второго порядка под действием давления 110 МПа. Сторона плиты 20 см, толщина h = 1 см. Механические свойства: модуль упругости Е = 206900 МПа, коэффициент Пуассона µ = 0,29. Верхние ребра не имеют вертикального смещения. В силу симметрии была рассмотрена четверть плиты. Критерием пластичности служит условие Губера–Мизеса (12). Функция нелинейного изотропного упрочнения имеет вид со следующими характеристиками: МПа, , . На рис. 3 показано распределение интенсивности пластических деформаций, на рис. 4 приведен график зависимости вертикального перемещения центральной точки плиты от давления. Рис. 3. Распределение интенсивности пластических деформаций Рис. 4. График зависимости вертикального перемещения центральной точки плиты от давления Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов № 12-01-00955, 12-01-97026, 12-01-31212.

About the authors

Ruslan Lavrentievich Davydov

Kazan Federal University

Email: ruslan.davydov@mail.ru
18, Kremlyovskaya st., 420008, Kazan, Russian Federation postgraduate, Department of Theoretical Mechanics, Kazan Federal University

Lenar Usmanovich Sultanov

Kazan Federal University

Email: Lenar.Sultanov@kpfu.ru
18, Kremlyovskaya st., 420008, Kazan, Russian Federation Ph. D. in Physical and Mathematical Sciences, Ass. Professor, Department of Theoretical Mechanics, Kazan Federal University

References

Statistics

Views

Abstract - 194

PDF (Russian) - 217