Вариант теории термопластичности
- Авторы: Бондарь ВС1, Даншин ВВ1, Кондратенко АА1
- Учреждения:
- Московский государственный машиностроительный университет (МАМИ)
- Выпуск: № 2 (2015)
- Страницы: 21-35
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/245
- DOI: https://doi.org/10.15593/perm.mech/2015.2.02
- Цитировать
Аннотация
Рассматриваются основные положения и уравнения теории термопластичности, относящейся к классу теорий пластического течения при комбинированном упрочнении. Тензор скоростей деформаций представляется в виде суммы тензоров скоростей упругой и пластической деформаций. Упругая деформация следует обобщенному закону Гука, распространенному на неизотермическое нагружение. Вводится поверхность нагружения, которая изотропно расширяется или сужается и смещается в процессе нагружения. Для радиуса поверхности нагружения формулируется эволюционное уравнение, учитывающее дополнительное изотропное упрочнение при непропорциональном нагружении, а также обобщенное на неизотермическое нагружение. В качестве параметра, характеризующего меру сложности процесса нагружения, принимается параметр Кадашевича-Мосолова, соответствующий углу между векторами скоростей деформаций и напряжений. Смещение поверхности нагружения описывается на основе модели Новожилова-Шабоши, подразумевающей, что полное смещение есть сумма смещений, для каждого из которых имеет место свое эволюционное уравнение. Проведенный ранее анализ петли пластического гистерезиса позволил выделить три типа микронапряжений (смещений) и сформулировать три типа эволюционных уравнений. Здесь эти эволюционные уравнения обобщены на неизотермическое нагружение. Для определения тензора скоростей пластической деформации используется ассоциированный (градиентальный) закон течения. Для жестких и мягких режимов нагружения получены выражения для определения скорости накопленной пластической деформации. Сформулированы условия упругого и упругопластического состояний. Для описания нелинейных процессов накопления повреждений вводятся кинетические уравнения накопления повреждений, где в качестве энергий, расходуемых на создание повреждений в материале, принимаются энергии, равные работам микронапряжений первого и второго типов на поле пластических деформаций. Здесь эти уравнения обобщены на неизотермическое нагружение. Выделяются материальные функции, замыкающие вариант теории, формулируется базовый эксперимент и метод идентификации материальных функций. Приводится описание верификации варианта теории термопластичности на широком спектре конструкционных сталей и сплавов и программ экспериментальных исследований. Новыми результатами работы являются адекватные описания в рамках одной теории следующих явлений: - посадка петли пластического гистерезиса при несимметричных жестких циклических нагружениях; - вышагивание (ratcheting) петли пластического гистерезиса при несимметричных мягких циклических нагружениях; - закономерности сложного нагружения как по плоским, так и пространственным траекториям; - эффекты дополнительного изотропного упрочнения при непропорциональных (сложных) циклических нагружениях; - эффекты нелинейного суммирования повреждений для произвольных процессов нагружения; - закономерности неизотермического нагружения.
Полный текст
Введение Математическое моделирование процессов деформирования и накопления повреждений при произвольных процессах сложного неизотермического нагружения в условиях повторности термосиловых воздействий в основном строится на вариантах теорий пластического течения при комбинированном упрочнении, обзор и анализ которых содержится в работах [1-18]. Основной проблемой построения этих вариантов является формулировка достаточно адекватных эволюционных уравнений для радиуса поверхности нагружения (изотропное упрочнение), для смещения центра поверхности нагружения (анизотропное упрочнение), а также формулировка кинетических уравнений накопления повреждений для произвольных процессов нагружения. Актуальной является и проблема обобщения теории на неизотермическое нагружение. Для описания изменения радиуса поверхности нагружения с учетом дополнительного изотропного упрочнения и неизотермического нагружения принимается эволюционное уравнение, предложенное в работах [4-6], где в качестве меры сложности процесса непропорционального нагружения принимается параметр Кадашевича-Мосолова [19], равный квадрату синуса угла между векторами скоростей напряжений и деформаций. Обоснование выбора этого параметра с учетом экспериментальных эффектов дополнительного изотропного упрочнения и разупрочнения приведено в работе [5]. Для описания смещения поверхности нагружения используется модель Новожилова-Шабоши [20, 21], подразумевающая, что полное смещение есть сумма смещений, для каждого из которых имеет место свое эволюционное уравнение. В качестве таких уравнений в настоящей работе принимаются уравнения Ишлинского-Прагера [22, 23], Амстронга-Фредерика-Кадашевича [24, 25 ], а также аналог уравнений Оно-Ванга [6, 26]. Для описания явлений посадки петли пластического гистерезиса и вышагивания (ratcheting) при несимметричных жестких и мягких циклических нагружениях на основе принципа симметрии циклических свойств [27] параметр, входящий в эволюционное уравнение для смещения (микронапряжения) первого типа, принимается зависящим [6, 27] не только от температуры, но и от накопленной пластической деформации. Для описания нелинейного процесса накопления повреждений формулируются [28] кинетические уравнения накопления повреждений, где в качестве энергий, расходуемых на создание повреждений в материале, принимаются работы микронапряжений первого и второго типов на поле пластических деформаций. Ответственность микронапряжений за процесс накопления повреждений следует из гипотезы Новожилова-Рыбакиной [29] о пропорциональности скорости накопления повреждений интенсивности микронапряжений. Обоснование применимости кинетических уравнений накопления повреждений на основе критериев работы микронапряжений первого и второго типов содержится в работах [3, 5, 9, 28, 30, 31]. Для определения материальных функций, замыкающих вариант теории термопластичности, формулируется базовый эксперимент и метод идентификации материальных функций. Приводится описание верификации [3, 5, 9, 27, 28, 30, 31] варианта теории термопластичности на широком спектре конструкционных сталей и сплавов и программ экспериментальных исследований. 1. Основные положения и уравнения теории Материал однороден и начально изотропен. Рассматриваются только поликристаллические конструкционные стали и сплавы. В процессе упругопластического деформирования в материале может возникать только пластическая деформационная анизотропия. Рассматриваются малые деформации при температурах, когда нет фазовых превращений, и скоростях деформаций, когда динамическими эффектами можно пренебречь. Случаи больших градиентов температуры не рассматриваются. Тензор скоростей деформаций представляется в виде суммы тензоров скоростей упругой и пластической деформаций: . (1) Упругие деформации следуют обобщенному закону Гука (2) где - соответственно модуль Юнга, коэффициент Пуассона, коэффициент температурного расширения, являющиеся функциями температуры Т; - тензор напряжений; - среднее напряжение; - символ Кронекера ( при при ). Полагается, что в пространстве составляющих тензора напряжений существует поверхность нагружения, разделяющая области упругого и упругопластического состояний. Поверхность нагружения изотропно расширяется или сужается и смещается в процессе нагружения. Уравнение поверхности нагружения принимается в следующем виде: (3) Здесь - девиатор активных напряжений [2]; - девиатор напряжений; - девиатор смещения (микронапряжений, добавочных напряжений, остаточных микронапряжений [1]); - размер (радиус) поверхности нагружения. Тензор характеризует анизотропное (направленное) упрочнение, а скаляр - изотропное упрочнение. Тензор и скаляр являются функционалами процесса нагружения. Для радиуса поверхности нагружения принимается следующее эволюционное уравнение: (4) Здесь - накопленная пластическая деформация (длина дуги траектории пластической деформации, параметр Одквиста); - интенсивность скоростей пластической деформации (скорость накопленной пластической деформации); - определяющие функции, которые выражаются через материальные. В случае отсутствия дополнительного изотропного упрочнения [3-6] радиус поверхности нагружения является функцией накопленной пластической деформации и температуры, (5) и, следовательно, . (6) Для описания дополнительного изотропного упрочнения при непропорциональном нагружении определяющая функция в уравнении (4) принимает [4-6] следующий вид: , (7) где - радиус поверхности нагружения при пропорциональном (простом) нагружении; - описывает дополнительное изотропное упрочнение при непропорциональном (сложном) нагружении. Для принимается [4-6, 32] следующее выражение: . (8) Здесь характеризует интенсивность (скорость) дополнительного упрочнения или разупрочнения, а - величину дополнительного упрочнения или разупрочнения; - параметр (мера) непропорциональности нагружения. Параметр при упрочнении и разупрочнении принимает различные значения, т.е. (9) Значение параметра при разупрочнении, т.е. , зависит [4-6, 32] от максимальной величины радиуса , достигнутой при упрочнении. Поэтому для принимается следующее выражение: . (10) Для параметра , характеризующего дополнительное упрочнение, можно также принять зависимость в долях от , т.е. (11) Здесь - модули дополнительного упрочнения и разупрочнения. В качестве параметра непропорциональности принимается параметр Кадашевича-Мосолова [19], соответствующий квадрату синуса угла между векторами скоростей деформаций и напряжений, обоснование выбора которого проведено в работе [5], , (12) Здесь - девиаторы деформаций и напряжений. Для максимального радиуса вводится [4, 6] следующее эволюционное уравнение: (13) При получении связи между определяющими функциями и материальными функциями при изменении температуры принимаются [1, 4, 6, 7, 8] неизменными относительные положения величин радиуса и максимального радиуса между их предельными значениями, т.е. (14) Тогда для и имеют место следующие выражения: (15) Итак, для описания изотропного упрочнения необходимо экспериментально определить: - функция изотропного упрочнения при пропорциональном (простом) нагружении; - модули дополнительного изотропного упрочнения и разупрочнения при непропорциональном (сложном) нагружении. Смещение поверхности нагружения описывается на основе модели Новожилова-Шабоши [20, 21], подразумевающей, что полное смещение есть сумма смещений, для каждого из которых имеет место свое эволюционное уравнение, (16) Проведенный [28] анализ петли пластического гистерезиса позволил выделить три типа микронапряжений (смещений) и сформулировать три типа эволюционных уравнений. В качестве первого эволюционного уравнения для микронапряжений первого типа принимается уравнение Ишлинского-Прагера [22, 23], обобщенное на неизотермическое нагружение, (17) Здесь определяющие функции и выражаются через экспериментально определяемые материальные функции. Для определяющей функции с учетом явления вышагивания (ratcheting) петли пластического гистерезиса при мягком несимметричном циклическом нагружении и явления посадки петли пластического гистерезиса при жестком несимметричном циклическом нагружении и на основании принципа симметрии циклических свойств [27] принимается следующее выражение: (18) Здесь - модуль анизотропного упрочнения; - модули вышагивания. При изменении температуры принимаются неизменными относительные положения величин микронапряжений первого типа между их предельными значениями, т.е. (19) Тогда для имеет место следующее выражение: (20) В качестве второго эволюционного уравнения для микронапряжений второго типа принимается уравнение Амстронга-Фредерика-Кадашевича [24, 25], обобщенное на неизотермическое нагружение, (21) Определяющие функции и выражаются через материальные следующим образом [3-6]: (22) Здесь - модули анизотропного упрочнения, определяемые экспериментально. При изменении температуры принимается неизменными относительные положения величин микронапряжений второго типа между их предельными значениями, т.е. (23) Тогда для имеет место следующее выражение: (24) Последующие эволюционные уравнения для микронапряжений третьего типа соответствуют простейшему аналогу [6, 28, 30] уравнений Оно-Ванга [26], обобщенному на неизотермическое нагружение, (25) Определяющие функции выражаются через материальные следующим образом [6, 28, 30]: (26) При изменении температуры принимаются неизменными относительные положения величин микронапряжений третьего типа между их предельными значениями, т.е. (27) Тогда для имеет место следующее выражение: (28) Окончательно уравнение для смещения поверхности нагружения с учетом (16), (17), (21), (25) будет иметь вид (29) (30) Пластические деформации определяются на основе ассоциированного с поверхностью (3) закона течения следующим образом: (31) Дифференцируя уравнение (3) по времени, подставляя в полученное выражение (29), (31) и далее разрешая относительно скорости накопленной пластической деформации , можно получить следующее уравнение для скорости накопленной пластической деформации при мягком нагружении, т.е. при заданных напряжениях: (32) Для получения уравнения при жестком нагружении (при заданных деформациях) на основании (1), (2) следует выражение (33) Далее подставляя (33) в (32) и разрешая относительно , можно получить следующее уравнение для скорости накопленной пластической деформации при жестком нагружении: (34) Для смешанных режимов нагружения уравнения для приводятся в работах [3-6, 9]. Условия упругого и упругопластического состояний имеют вид (35) Здесь скорость накопленной пластической деформации задается выражениями (32) или (34) или любым другим выражением, связывающим скорость накопленной пластической деформации и скорости напряжений и деформаций (смешанные режимы нагружения). Для описания нелинейных процессов накопления повреждений вводятся кинетические уравнения накопления повреждений, базирующиеся на энергетическом принципе, где в качестве энергий, расходуемых на создание повреждений в материале, принимаются энергии, равные работам микронапряжений первого и второго типов на поле пластических деформаций. Эти уравнения аналогичны [28], но здесь они обобщены на неизотермическое нагружение. Тогда кинетические уравнения накопления повреждений будут иметь следующий вид: (36) (37) Здесь - мера повреждения; - энергии разрушения соответственно при вышагивании петли гистерезиса, при пропорциональном и непропорциональном нагружении; и - функция и параметр нелинейности процесса накопления повреждений; - интенсивность микронапряжений второго типа. Критерием разрушения будет достижение повреждением предельного значения, обычно принимаемого равным единице. Таким образом, данный вариант теории термопластичности замыкают следующие материальные функции, подлежащие экспериментальному определению: - упругие параметры; - модули анизотропного упрочнения; - модули вышагивания; - модули анизотропного упрочнения, соответствующие аналогу модели Оно-Ванга; - функция изотропного упрочнения; - модули дополнительного изотропного упрочнения; - энергия разрушения при пропорциональном нагружении ; - энергия разрушения при непропорциональном нагружении ; - параметр нелинейности процесса накопления повреждений; - энергия разрушения при вышагивании. 2. Базовый эксперимент и метод идентификации материальных функций Материальные функции определяются по результатам испытаний в условиях упругопластического одноосного пропорционального и двухосного непропорционального напряженных состояний при различных уровнях температуры. Базовый эксперимент включает в себя следующий набор данных: - упругие параметры, которые определяются традиционными методами; - диаграмма пластического деформирования при растяжении до деформации ; - циклические диаграммы при симметричном растяжении-сжатии при постоянной амплитуде деформации ; - циклические диаграммы при несимметричном растяжении-сжатии при постоянной амплитуде деформации и средней деформации цикла ; - данные по малоцикловой усталости при одноблочном и двухблочном жестком симметричном циклическом нагружении; - данные по малоцикловой усталости при одноблочном мягком несимметричном циклическом нагружении; - диаграмма максимальных значений интенсивности напряжений на цикле от накопленной пластической деформации при непропорциональном циклическом нагружении по траектории деформаций в виде окружности радиуса до стабилизации дополнительного упрочнения и последующем пропорциональном циклическом нагружении до стабилизации разупрочнения; - данные по усталостному разрушению при непропорциональном циклическом нагружении по траекториям деформаций в виде окружностей с разными радиусами. Метод идентификации материальных функций подробно изложен в работах [3-6, 27, 28, 30] и применяется здесь для каждого из рассматриваемых уровней температуры. 3. Верификация варианта теории термопластичности Далее рассматривается верификация на основе сопоставления расчетных и экспериментальных результатов для различных процессов пропорционального (простого) и непропорционального (сложного), изотермического и неизотермического режимов нагружения. В работах [3, 5, 9, 27, 28, 30, 31] приводится верификация варианта теории при пропорциональных стационарных и нестационарных, симметричных и несимметричных, жестких и мягких режимах циклического изотермического нагружения. Иллюстрируется адекватное описание эффекта посадки петли пластического гистерезиса при жестком несимметричном циклическом нагружении, а также явление вышагивания (ratcheting) при мягком несимметричном циклическом нагружении. В этих же работах анализируются процессы нелинейного накопления повреждений при одноблочных и многоблочных режимах жестких и мягких циклических нагружений. Рассматриваются процессы от малоцикловой до многоцикловой усталости (циклов). Неизотермические процессы при пропорциональном (простом) циклическом нагружении рассматриваются в работах [3, 9]. Показывается адекватное описание теорией экспериментальных циклических диаграмм и разрушения в условиях малоцикловой усталости. Процессы непропорционального (сложного) нагружения достаточно полно анализируются в монографии [5], где рассматриваются плоские и пространственные траектории деформаций в широком диапазоне кривизн и круток от малых до больших. В этой же монографии показывается адекватное описание теорией эффекта дополнительного изотропного упрочнения при непропорциональных (сложных) циклических нагружениях. Там же рассматривается разрушение при циклических пропорциональных и непропорциональных режимах нагружения от малоцикловой до многоцикловой усталости. Снижение долговечности в условиях непропорционального циклического нагружения по сравнению с пропорциональным при одинаковых размахах деформаций достигает почти порядка, это показывает и эксперимент, и расчет. Соответствие расчетных и экспериментальных результатов для широкого спектра конструкционных сталей и сплавов и режимов нагружения говорит о достаточной работоспособности предложенного варианта теории термопластичности. Заключение Сформулированы основные положения и уравнения варианта теории термопластичности, адекватно описывающие кинетику напряженно-деформированного состояния и нелинейные процессы накопления повреждений при произвольном сложном неизотермическом нагружении. Выделены материальные функции, замыкающие теорию, сформулированы базовый эксперимент и метод идентификации материальных функций. Приводится описание верификации варианта теории термопластичности на широком спектре конструкционных сталей и сплавов и программ экспериментальных исследований. Следует отметить, что здесь в рамках одной теории описываются, экспериментально полученные в последнее время закономерности сложного нагружения как по плоским, так и по пространственным траекториям деформаций, эффекты посадки и вышагивания (ratcheting) петли пластического гистерезиса при несимметричных циклических нагружениях, закономерности неизотермического нагружения и нелинейного суммирования повреждений, а также эффекты дополнительного изотропного упрочнения при непропорциональных циклических нагружениях.Об авторах
В С Бондарь
Московский государственный машиностроительный университет (МАМИ)
В В Даншин
Московский государственный машиностроительный университет (МАМИ)
А А Кондратенко
Московский государственный машиностроительный университет (МАМИ)
Список литературы
- Термопрочность деталей машин: справочник / под ред. И.А. Биргера, Б.Ф. Шорра. - М.: Машиностроение, 1975. - 455 с.
- Новожилов В.В., Кадашевич Ю.И. Микронапряжения в конструкционных материалах. - Л.: Машиностроение. - М.: Изд-во МГУ, 1990. - 224 с.
- Бондарь В.С. Неупругое поведение и разрушение материалов и конструкции при сложном неизотермическом нагружении: дис. … д-ра физ.-мат. наук. - М.: Изд-во МАМИ, 1990. - 314 с.
- Бондарь В.С. Неупругость. Варианты теории.- М.: Физматлит, 2004. - 144 с.
- Бондарь В.С., Даншин В.В. Пластичность. Пропорциональные и непропорциональные нагружения. - М.: Физматлит, 2008. - 176 с.
- Bondar V.S. Inelasticity. Variants of the theory. - New York: Begell House, 2013. - 194 p.
- Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругопластических сред с повреждениями. - М.: Физматлит, 2008. - 424 с.
- Нелинейная механика материалов / Ж. Бессон, Ж. Каето, Ж.-Л. Шабоши, Т.С. Форест. - Санкт-Петербург: Изд-во Политехн. ун-та, 2010. - 397.
- Макаров Д.А. Математическое моделирование процессов неизотермического неупругого деформирования и накопления повреждений в конструкционных материалах: дис. … канд. физ.-мат.наук. - М.: Изд-во МАМИ, 2005. - 108 с.
- Bari S., Hassan T. An advancement in cyclic plasticity modeling for multiaxial ratcheting simulation // International Journal of Plasticity. - 2002. - Vol. 18. - P. 873-894.
- Uniaxial ratcheting and fatigue failure of tempered 42CrMo steel: Damage evolution and damage-coupled viscoplastic constitutive model / G. Kang, Y. Liu, J. Ding, Q. Gao // Int. J. of Plasticity. - 2009. - Vol. 25. - P. 838-860.
- Kan Q., Kang G. Constitutive model for uniaxial transformation ratcheting of super-elastic NiTi shape memory alloy at room temperature // Int. J. of Plasticity. - 2009. doi: 10.1016/j.ijplas.2009.08.005
- Chaboche J.-L. A review of some plasticity and viscoplasticity constitutive theories // Int. J. of Plasticity. - 2008. - Vol. 24. - P. 1642-1692.
- Rahman S.M., Hassan T., Corona E. Evaluation of cyclic plasticity models in ratcheting simulation of straight pipes under cyclic bending and steady internal pressure // Int. J. of Plasticity. - 2008. - Vol. 24. - P. 1756-1791.
- Abdel-Karim M. Modified kinematic hardening rules for simulations of ratcheting // Int. J. of Plasticity. - 2009. - Vol. 25. - P. 1560-1587.
- Abdel-Karim M. An evaluation for several kinematic hardening rules on prediction of multiaxial stress-controlled ratchetting // Int. J. of Plasticity. - 2010. - Vol. 26. - P. 711-730.
- Dafalias Y.F., Feigenbaum H.P. Biaxial ratchetting with novel variations of kinematic hardening // Int. J. of Plasticity. - 2011. - Vol. 27. - P. 479-491.
- Chaboche J.-L., Kanouté P., Azzouz F. Cyclic inelastic constitutive equations and their impact on the fatigue life predictions // Int. J. of Plasticity. - 2012. - Vol. 35. - P. 44-66.
- Кадашевич Ю.И., Мосолов А.Б. О соотношениях эндохронной теории пластичности с «новой» мерой внутреннего времени при сложном циклическом нагружении // Технология легких сплавов. - 1990. - № 3. - С. 32-36.
- Новожилов В.В. О сложном нагружении и перспективах феноменологического подхода к исследованию микронапряжений // ПММ. - 1964. - Т. 28. - Вып. 3. - С. 393-400.
- Chaboche J.-L., Dang-Van K., Cordier G. Modelization of the strain memory effect on the cyclic hardening of 316 stainless steel // Proceedings of the 5th International Conference on SMiRT. Div L, Berlin, 1979. - P. No. L. 11/3.
- Ишлинский А.Ю. Общая теория пластичности с линейным упрочнением // Укр. матем. журн. - 1954. - Т. 6. - Вып. 3. - С. 314-324.
- Prager W. A new method of analyzing stresses and strains in work hardening plastic colids // ASME J. Appl. Mech. - 1956. - Vol. 23. - P. 493-496.
- Amstrong P.J., Frederick C.O. A mathematical represention of the multiaxial bauscinger effect // CEGB Report No. RD/B/N/ 731. - 1966.
- Кадашевич Ю.И. О различных тензорно-линейных соотношениях в теории пластичности // Исследования по упругости и пластичности. - Л.: Изд-во ЛГУ, 1967. - Вып. 6. - С. 39-45.
- Ohno N., Wang J.-D. Kinematic hardening rules with critical state of dynamic recovery, part 1: formulations and basic features for ratcheting behavior // International Journal of Plasticity. - 1993. - Vol. 9. - P. 375-390.
- Бондарь В.С. Некоторые новые результаты исследования пластичности материалов при сложном нагружении // Упругость и неупругость. - М: ЛЕНАНД, 2006. - С. 94-109.
- Бондарь В.С., Даншин В.В., Макаров Д.А. Математическое моделирование процессов деформирования и накопления повреждений при циклических нагружениях // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 2. - С. 125-152.
- Новожилов В.В., Рыбакина О.Г. О перспективах построения критерия прочности при сложном нагружении // Прочность при малом числе нагружения. - М.: Наука, 1969. - С. 71-80.
- Бондарь В.С., Бурчаков С.В., Даншин В.В. Математическое моделирование процессов упругопластического деформирования и разрушения материалов при циклических нагружениях // Проблемы прочности и пластичности: межвуз. сб. Вып. 72. - Нижний Новгород: Изд-во Нижегород. гос. ун-та, 2010. - С. 18-27.
- Бондарь В.С., Даншин В.В., Семенов П.В. Нелинейные процессы накопления повреждений при нестационарных циклических нагружениях // Проблемы прочности и пластичности. - 2012. - Вып. 75. - Ч. 2. - С. 96-104.
- Benallal A., Marquis D. Constitutive Equations for No proportional. Cyclic Elasto-Viscoplasticity // Journal of Engineering Materials and Technology. - 1987. - Vol. 109. - P. 326-337.