Моделирование процессов упругопластического деформирования материалов по многозвенным кусочно-ломаным прямолинейным траекториям
- Авторы: Зубчанинов ВГ1, Алексеев АА1, Гультяев ВИ1
- Учреждения:
- Тверской государственный технический университет
- Выпуск: № 3 (2017)
- Страницы: 203-215
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/143
- DOI: https://doi.org/10.15593/perm.mech/2017.3.12
- Цитировать
Аннотация
Рассматривается математическое моделирование процессов упругопластического деформирования материала сталь 45 в векторном пространстве А.А. Ильюшина по многозвенным неаналитическим траекториям вида плоских многозвенных ломаных траекторий, состоящем из четырех прямолинейных кусочно-ломаных участков при его одновременном растяжении либо сжатии с кручением. При моделировании использовалась математическая модель теории процессов для плоских траекторий с уточненными аппроксимациями функционалов процесса, содержащая все необходимые параметры сложного нагружения для плоских траекторий. Основные уравнения математической модели приведены к задаче Коши, для численного решения которых и получения расчетных результатов использовался метод Рунге-Кутты четвертого порядка точности. Приведена методика определения материальных параметров аппроксимаций в выражениях функционалов на основе обработки экспериментальных данных по базовым траекториям типа «смещенного веера» двузвенных траекторий. Экспериментальные исследования проведены на автоматизированном испытательном комплексе СН-ЭВМ в лаборатории механических испытаний Тверского государственного технического университета. В качестве образцов для экспериментальных исследований были использованы тонкостенные цилиндрические оболочки из стали 45 в состоянии поставки. Численные теоретические расчеты сравниваются с экспериментальными данными для оценки достоверности математической модели теории процессов при сложном нагружении для данного класса траекторий деформирования. Предлагаемая математическая модель теории процессов дала результаты, хорошо соответствующие данным эксперимента, как по скалярным, так и векторным свойствам материала. Это подтверждает достаточную для практических задач достоверность расчетных данных и точность построенных аппроксимаций функционалов процессов используемой математической модели теории процессов.
Полный текст
1. Основные положения Тензоры напряжений и деформаций, характеризующие напряженно-деформированное состояние (НДС) точки тела, можно разложить на шаровые тензоры и девиаторы [1-3]: , , (1) где - символ Кронекера; , , , , (2) - модули шаровых тензоров и девиаторов соответственно; , , , , (3) - компоненты девиаторов и направляющих тензоров напряжений и деформаций соответственно. При простом пропорциональном нагружении направляющие тензоры напряжений и деформаций и определяющие соотношения теории малых упругопластических деформаций, предложенной в [1], с учетом упругости объемной деформации имеют вид , , , (4) где - модуль объемной упругости; - пластический модуль сдвига; - универсальная единой диаграммы упрочнения материалов Роша и Эйхингера при простом нагружении, определяющая только скалярные свойства материалов для произвольного НДС. При сложном нагружении направляющие тензоры и учет в определяющих соотношениях только скалярных свойств материалов недостаточен. В работах А.А. Ильюшина [2-4] и В.Г. Зубчанинова [11] в линейном совмещенном евклидовом пространстве с ортонормированным неподвижным базисом было использовано представление тензоров и в виде векторов: (5) где , - векторы напряжений и деформаций в одномерном подпространстве объемного растяжения и сжатия с гидростатической осью, характеризуемой единичным вектором ; , - векторы напряжений и деформаций формоизменения в пятимерном девиаторном подпространстве . Координаты векторов связаны с компонентами тензоров и девиаторов взаимно однозначными преобразованиями [2-6] (6) Модули векторов в пятимерном подпространстве равны модулям девиаторов напряжений и деформаций соответственно: , , , . (7) В работе [3] А.А. Ильюшин представил основное соотношение постулата изотропии, которое в дальнейшем в работах [2, 7-9] было дополнено соотношением для отображения влияния векторных свойств материалов. Общие определяющие соотношения теории процессов получены в работах [5, 6]. Они отражают связь между векторами напряжений и деформаций формоизменения с учетом скалярных и векторных свойств материалов. Для случая плоских траекторий без учета Э2-эффекта определяющие соотношения в скалярной форме имеют вид [5, 6] (8) где , - функционалы процесса деформирования, зависящие от параметров сложного нагружения; s - длины дуги траектории деформирования, ее кривизны и углов излома . Угол сближения характеризует направление по отношению к касательной к траектории деформирования в каждой ее точке. Этот угол отражает влияние векторных свойств материала на процесс деформирования. 2. Математическая модель теории процессов в плоских задачах В математической модели в дополнение к определяющим соотношениям (8) используются универсальные аппроксимации функционалов [5] (9) (10) где - универсальная функция нагружения Одквиста-Ильюшина для процессов, близких к простым, без учета их истории; - приращение длины дуги траектории после ее излома в некоторой точке K; - разница в точках излома между значениями универсальной функции Одквиста-Ильюшина и расчетным значением модуля вектора напряжений ; ; индекс «нолик» у пластического модуля сдвига соответствует значению в точке излома траектории; ; (11) - функция, описывающая после излома траектории скалярный нырок напряжений при сложной разгрузке и последующем вторичном пластическом деформировании, и ее производная; ; (12) - функция, учитывающая ориентацию вектора напряжений в процессе деформирования, и ее значение в точке излома при значении угла сближения каждого из участков неаналитической траектории; - материальные параметры материала, для определения которых ниже предложена методика на основе обработки экспериментальных данных для каждого материала. Для аппроксимации универсальной функции упрочнения Одквиста-Ильюшина при простом нагружении использовались выражения [5] (13) Здесь ; - предел текучести при растяжении; - граница участков диаграммы деформирования, разделяющая упругую часть диаграммы и площадку текучести () от участка самоупрочнения материала (); - упругий модуль сдвига; - экспериментально определяемые параметры материала из опытов на простое нагружение по типу «центрального веера» прямолинейных траекторий [5]. При заданных начальных условиях для координат вектора деформаций и начальных значениях угла основные уравнения (8) с конкретизированными функционалами (9), (10) приводятся к задаче Коши, для численного решения которой и определения координат вектора напряжений и угла сближения использовался метод Рунге-Кутты четвертого порядка точности в программном приложении MathWorks MatLab. Ранее данная математическая модель была использована [10] для описания процессов деформирования только по двузвенным ломаным прямолинейным траекториям. Другие заслуживающие внимания идеи построения определяющих соотношений и проведения базовых экспериментальных исследований представлены в работах [12-30]. 3. Методика определения материальных параметров аппроксимаций Для определения материальных параметров использован базовый опыт траекторий «смещенного веера» для кусочно-ломаной прямолинейной траектории с углом излома на 135 градусов. Общая диаграмма прослеживания процесса деформирования такой траектории в сравнении с универсальной функцией упрочнения Одквиста-Ильюшина представлена на рис. 1. Рис. 1. К определению материальных параметров Fig. 1. To the definition of material parameters Разность между значениями и действительными значениями модуля напряжений на диаграмме при определяется так: (14) В точке минимума нырка напряжений на диаграмме имеем , следовательно, . Из (14) с учетом (11) имеем . (15) В точке при из (14) с учетом (11) получаем . (16) Возьмем отношение , где после преобразований может быть получена функция , откуда находится параметр b, выражаемый через W-функцию Ламберта . (17) Для вычисления параметра при вычисленном значении по (17) предлагается использовать выражение (16), откуда (18) Для вычисления параметра используется значение производной в точке излома траектории при , (19) где из (13) . При в точке K с учетом (12) имеем при упругой разгрузке для получается , откуда Тогда с учетом (18) (20) При найденном значении по (18) и (20) при различных показателях степени может быть найден спектр значений параметров , , выбор которых определяется наилучшим соответствием данным базового эксперимента. 4. Результаты экспериментального исследования и математического моделирования Экспериментальное исследование было выполнено на автоматизированном испытательном комплексе СН-ЭВМ имени А.А. Ильюшина в лаборатории механических испытаний кафедры «Сопротивление материалов, теории упругости и пластичности» Тверского государственного технического университета. В качестве образцов для экспериментальных исследований были использованы тонкостенные цилиндрические оболочки из стали 45 в состоянии поставки, имеющие в рабочей части: длину мм, толщину мм и диаметр срединной поверхности мм. Начальная изотропия материала образцов с достаточной степенью точности была подтверждена в опытах на простое нагружение (растяжение, сжатие и кручение), где после обработки этих диаграмм были приняты следующие значения материальных параметров для стали 3 в аппроксимации (13): , , , , , , . Для оценки достоверности расчетных результатов по предлагаемой математической модели в пространстве деформаций девиаторного подпространства реализована нетривиальная неаналитическая траектория по типу «песочные часы», состоящая из четырех прямолинейных кусочно-ломаных участков (рис. 2). На первом прямолинейном участке реализовывалось пропорциональное кручение по координате до значения ; на втором участке при изломе траектории на угол реализовывалось комбинированное растяжение и кручение до значений ; на третьем участке при изломе траектории на угол при осуществлялось кручение до на четвертом участке при изломе вновь на угол реализовывалось комбинированное сжатие и кручение с приходом в начало координат. На рис. 3 приведен отклик на реализованную траекторию деформирования в виде траектории нагружения в плоскости совмещенного девиаторного подпространства . На рис. 4 приведены результаты расчета и экспериментальные данные для диаграмм прослеживания процессов сложного нагружения , характеризующих скалярные свойства материалов в зависимости от длины дуги сложного деформирования, на рис. 5 - диаграммы для углов сближения , характеризующие векторные свойства материалов. На рис. 6, 7 приведены локальные диаграммы деформирования растяжения-сжатия по компонентам и чистого сдвига по компонентам . На рис. 2-7 цифрами 1, 2, 3, 4 обозначены точки начала соответствующих участков реализованной кусочно-ломаной траектории. Экспериментальные данные на рисунках отмечены точками; модельные расчетные данные отмечены синим цветом. В расчете на всех участках траектории использовались значения материальных параметров, найденные по (17), (18), (20) при показателе степени , которые были равны соответственно: , МПа, . Показатель степени на втором участке принимался равным 4. На последующих участках его значение подбиралось из условия соответствия глубины «нырка напряжений» данным эксперимента и принималось - для третьего участка и - для четвертого. Показатель степени на всех участках принимался . В расчетах также было учтено, что начальное отклонение вектора напряжений не всегда соответствует углу излому траектории . На первом изломе в точке 2 после участка простого нагружения . Перед вторым и третьим изломом вектор напряжений имел отклонение, что учитывалось при численном решении в виде начального условия для участка траектории , где расчетное значение угла сближения в конце предшествующего участка перед изломом траектории. Перед вторым изломом при переходе от второго к третьему участку траектории расчетное значение при этом излом траектории в точке 3 производился в сторону отклонения вектора напряжений, поэтому в расчете принималось Перед третьим изломом (переход от третьего к четвертому участку траектории) расчетное значение а излом траектории в точке 4 производился в сторону, противоположную отклонению вектора напряжений, поэтому здесь Рис. 2. Траектория деформирования на плоскости Fig. 2. Deformation trajectory on the plane Рис. 3. Отклик по напряжениям на плоскости Fig. 3. Response of tension on the plane Рис. 4. Глобальная диаграмма деформирования Fig. 4. General chart of deformation Рис. 5. Диаграмма характеристики векторных свойств материала Fig. 5. Chart that characterize the vector material properties Рис. 6. Локальная диаграмма деформирования Fig. 6. Local chart of deformation Рис. 7. Локальная диаграмма деформирования Fig. 7. Local chart of deformation Как видно, принятые для модели данные достаточно хорошо соответствуют экспериментальным как по скалярным, так и по векторным свойствам. Заключение Сравнение результатов теоретических численных расчетов с результатами экспериментального исследования по заданной сложной многозвенной траектории показало достаточную для практических задач достоверность расчетных данных и точность построенных аппроксимаций функционалов процессов используемой математической модели теории процессов [3-5].Об авторах
В Г Зубчанинов
Тверской государственный технический университет
А А Алексеев
Тверской государственный технический университет
В И Гультяев
Тверской государственный технический университет
Список литературы
- Ильюшин А.А. Пластичность. Упругопластические деформации. - М.: Гостехиздат, 1948. - 376 с.
- Ильюшин А.А. Механика сплошной среды. - М.: Изд-во МГУ, 1971. - 310 с.
- Ильюшин А.А. Пластичность. Основы общей математической теории. - М.: Изд-во АН СССР, 1963. - 273 с.
- Ильюшин А.А. Труды (1946-1966). Т. 2. Пластичность. - М.: Физматлит, 2004. - 480 с.
- Зубчанинов В.Г. Механика процессов пластических сред. - М.: Физматлит, 2010. - 352 с.
- Зубчанинов В.Г. Устойчивость и пластичность. Т.2. Пластичность. - М.: Физматлит, 2008. - 336 с.
- Ленский В.С. Гипотеза локальной определенности в теории пластичности // Изв. АН СССР, ОТН. - 1962. - № 5. - С. 154-158.
- Дао Зуй Бик. О гипотезе локальной определенности в теории пластичности // Вестн. МГУ. Сер. «Математика, механика». - 1965. - № 2. - С. 67-75.
- Васин Р.А, Ильюшин А.А. Об одном представлении законов упругости и пластичности в плоских задачах // Изв. АН СССР. МТТ. - 1983. - № 4. - С. 114-118.
- Зубчанинов В.Г., Алексеев А.А., Гультяев В.И. Численное моделирование процессов сложного упругопластического деформирования стали по двузвенным ломаным траекториям // Проблемы прочности и пластичности: межвуз. сб. - Н. Новгород: Изд-во Нижегород. ун-та, 2014. - Вып. 76. - Ч. 1. - С. 18-25.
- Зубчанинов В.Г. Постулат изотропии и закон сложной разгрузки сплошных сред // Изв. РАН. Механика твердого тела. - 2011. - № 1. - С. 27-37. doi: 10.3103/S0025654411010043
- Ленский B.C. Экспериментальная проверка основных постулатов общей теории упругопластических деформаций // Вопросы теории пластичности. - М.: Изд-во АН СССР, 1961. - С. 58-82.
- Васин Р.A. Некоторые вопросы связи напряжений и деформаций при сложном нагружении // Упругость и неупругость. Вып. 1. - М.: Изд-во МГУ, 1971. - С. 59-126.
- Некоторые экспериментальные данные об общем законе пластичности Ильюшина / И. Охаши, М. Токуда, И. Курита, Т. Сузуки // Изв. АН СССР. МТТ. - 1981. - № 6. - С. 53-64.
- Ленский В.С., Ленский Э.В. Трехчленное соотношение общей теории пластичности // Изв. АН СССР. МТТ. - 1985. - № 4. - С. 111-115.
- Ленский В.С., Машков И.Д. Проверка законов в трехмерном пространстве девиатора деформаций // Упругость и неупругость. - М.: Изд-во МГУ, 1971. - № 2. - С. 158-166.
- Васин Р.А. Идеи Ильюшина А.А., обогатившие теорию пластичности // Упругость и неупругость: материалы междунар. науч. симпозиума по проблемам механики деформируемых тел, посвященного 105-летию со дня рождения А.А.Ильюшина. - М.: Изд-во МГУ, 2016. - С. 66-69.
- Васин Р.А. Теория упругопластических процессов и исследование структурно-механических свойств материалов // Изв. РАН. МТТ. - 2011. - № 1. - С. 19-26.
- Васин Р.А. О «памяти» материала в теории упругопластических процессов // Изв. Тул. гос. ун-та. Естественные науки. - 2013. - № 2-2. - С. 59-70.
- Муравлев А.В. Экспериментальное построение функционалов пластичности для траекторий деформаций типа двухзвенных ломаных в опытах на сплошных цилиндрических областях // Вестн. Моск. ун-та. Сер. 1: Математика. Механика. - 1996. - № 5. - С. 74-80.
- Муравлев А.В. Обобщение постулата изотропии теории упругопластических процессов А.А. Ильюшина на конечные деформации // Проблемы прочности, пластичности и устойчивости в механике деформируемого твердого тела: материалы VIII Междунар. науч. симпозиума. - Тверь: Изд-во Твер. гос. техн. ун-та, 2015. - С. 25-29.
- Муравлёв А.В., Девятов А.С. Развитие теории упругопластических процессов А.А. Ильюшина и экспериментально-теоретических методов исследования вязкопластических свойств материалов при конечных деформациях // Проблемы машиностроения и автоматизации. - 2016. - № 1. - С. 84-90.
- Жуков А.М. Некоторые особенности поведения металлов при упругопластическом деформировании // Вопросы теории пластичности. - М.: Изд-во АН СССР, 1961. - C. 30-57.
- Жуков А.М. О свойствах запаздывания в общей теории пластичности // Изв. АН. МТТ. - 1992. - № 5. - С. 110-119.
- Трусов П.В., Швейкин А.И., Янц А.Ю. Основные положения теории упругопластических процессов А.А. Ильюшина: анализ с позиций многоуровневого моделирования для случая больших градиентов перемещений // Упругость и неупругость: материалы междунар. науч. симпозиума по проблемам механики деформируемых тел, посвященного 105-летию со дня рождения А.А.Ильюшина. - М.: Изд-во МГУ, 2016. - С. 119-126.
- Дегтярев В.П. Пластичность и ползучесть в машиностроительных конструкциях. - М.: Машиностроение, 1967. - 130 с.
- Шевченко Ю.Н., Бабешко М.Е., Терехов Р.Г. Термовязкоупругопластические процессы сложного деформирования элементов конструкций. - Киев: Наук. думка, 1992. - 328 с.
- Молодцов И.Н., Бабаева Д.О. Некоторые вопросы верификации теории упругопластических процессов при сложном нагружении // Проблемы машиностроения и автоматизации. - 2016. - № 3. - С. 98-105.
- Зубчанинов В.Г., Алексеев А.А., Алексеева Е.Г. Математическое моделирование процессов пластического деформирования материалов по сложным плоским траекториям // Materials Physics and Mechanics (MPM). - 2015. - Т. 24, № 2. - С. 107-118.
- Зубчанинов В.Г., Алексеев А.А., Алексеева Е.Г. Проверка постулата изотропии и численное моделирование процессов деформирования материалов на сложных гладких траекториях // Materials Physics and Mechanics. - 2016. - Т. 29, № 2. - С. 150-157.