Version of the theory of thermoplasticity
- Authors: Bondar VS1, Danshin VV1, Kondratenko AA1
- Affiliations:
- Moscow State University of Mechanical Engineering (MAMI)
- Issue: No 2 (2015)
- Pages: 21-35
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/245
- DOI: https://doi.org/10.15593/perm.mech/2015.2.02
- Cite item
Abstract
The paper explains the basic concepts and equations of the theory of thermoplasticity belonging to the class of theories of plastic flow in the combined hardening. The tensor strain rate is represented as a sum of tensors of the velocities of elastic and plastic deformations. The elastic deformation follows the generalized Hooke's law distributed to non-isothermal loading. The surface loading is introduced which isotropically widens or narrows and shifts in the process of loading. For the radius of the surface loading the authors formulated evolutionary equation taking into account the additional isotropic hardening under non-proportional loading, also generalized to non-isothermal loading. We have taken the parameter of Kadashevich-Mosolova corresponding to the angle between the velocity vectors of strain and stress as the parameter characterizing the degree of complexity of the process of loading. Surface displacement loading is described based on the model of Novozhilov-Chaboche implying that the total displacement is the sum of the displacements and each displacement has its own evolutionary equation. An earlier analysis of a loop of plastic hysteresis allowed distinguishing three types of micro-strains (displacements) and formulating three types of evolution equations. Here these evolutionary equations are generalized to the non-isothermal loading. To determine the axial velocity of plastic deformation we used associate (gradient) flow law. It became possible to determine expressions for speed of the accumulated plastic strain for hard and soft modes of loading. The authors have formulated conditions of elastic and elastic-plastic conditions. The kinetic equation of damage accumulation is introduced for the description of nonlinear processes of damage accumulation. Here, energies equal to the work of microstresses of the first and second types to the field of plastic deformations are accepted as the energy spent on creating damage in the material. These equations are generalized to the non-isothermal loading. We have highlighted material options completing theory option, formulated the basic experiment and the method of identifying material functions. The paper describes the verification of thermoplasticity theory on a wide range of structural steels and alloys and programs of experimental studies. New results have adequate descriptions within one theory of the following phenomena: - landing loop of plastic hysteresis in nonsymmetrical rigid cyclic loading; - ratcheting of loop plastic hysteresis in nonsymmetrical soft cyclic loading; - the regularities of complex loading as on planar or spatial trajectories; - the effects of additional isotropic hardening under disproportionate (complex) cyclic loading; - the effects of non-linear summation of damage to arbitrary loading processes; - the patterns of non-isothermal loading.
Full Text
Введение Математическое моделирование процессов деформирования и накопления повреждений при произвольных процессах сложного неизотермического нагружения в условиях повторности термосиловых воздействий в основном строится на вариантах теорий пластического течения при комбинированном упрочнении, обзор и анализ которых содержится в работах [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) петли пластического гистерезиса при несимметричных циклических нагружениях, закономерности неизотермического нагружения и нелинейного суммирования повреждений, а также эффекты дополнительного изотропного упрочнения при непропорциональных циклических нагружениях.About the authors
V S Bondar
Moscow State University of Mechanical Engineering (MAMI)
V V Danshin
Moscow State University of Mechanical Engineering (MAMI)
A A Kondratenko
Moscow State University of Mechanical Engineering (MAMI)
References
- Термопрочность деталей машин: справочник / под ред. И.А. Биргера, Б.Ф. Шорра. - М.: Машиностроение, 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.