The Experimental and Theoretical Study of Large Deformations of Cylindrical Samples from Steel 09G2S with Stress Concentrators under Tension-Torsion Loading to Failure
- Authors: Bazhenov VG1, Kazakov DA1, Nagornykh EV1, Osetrov DL1, Ryabov AA2
- Affiliations:
- Researcher Institute of Mechanics, National Research Lobachevsky State University of Nizhni Novgorod
- FGUP RFYAC-VNIIEF
- Issue: No 4 (2018)
- Pages: 69-81
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/481
- DOI: https://doi.org/10.15593/perm.mech/2018.4.06
- Cite item
Abstract
The paper presents the experimental and numerical results of elastoplastic deformation of specimens with a cylindrical working part and specimens with a stress concentrator made of 09G2C steel under monotonic kinematic tension-torsion loadings prior to failure, taking into account large deformations and inhomogeneity of the stress-strain state (SSS). A complete system of equations describing generalized axisymmetric torsion problems is written in a cylindrical coordinate system. Kinematic relations are formulated at speeds in the metric of the current state, which makes it possible to describe large deformations. The equation of motion of a continuous medium follows from the equation of the balance of virtual capacities. Since we consider active loading processes close to proportional ones, we describe the elastoplastic properties of materials by the flow theory with nonlinear isotropic hardening. The system of equations, supplemented by kinematic boundary and initial conditions, is solved by the finite element method in a combination with an explicit integration scheme of the "cross" type. On the basis of the experimental-calculation method, true deformation diagrams for tension and torsion were constructed to supplement the plasticity model with material functions. These diagrams differ significantly with deformations of more than 15 %, which is caused by the sensitivity of the plastic properties of the test material to the form of the stressed state. To describe the tension-torsion, the dependence of the deformation diagram on the form of the stressed state is introduced. The refined deformation diagram is a linear combination of tension and torsion diagrams, which coefficients depend on the parameter of the form of the stressed state, the parameter of the stress triaxiality or the Nadai-Lode stress parameter. There is a good correspondence of the numerical results with the experiment on the integral characteristics (axial force from axial displacement and torque from the twist angle). The analysis of the SSS parameters in the neck and circular concentrator is performed. The estimated the peculiarities of the mutual influence of tension and torsion on the localization of deformations and destruction of cylindrical samples and samples with a stress concentrator.
Full Text
Введение Элементы конструкций в штатных и нештатных режимах работы могут подвергаться многокомпонентному силовому воздействию, которое приводит к сложному напряженному состоянию. В теории пластичности полагается, что закон пластического упрочнения при нагружениях, близких к простым, не зависит от вида напряженного состояния [1, 2]. Многочисленные экспериментальные проверки гипотезы «единой кривой» при простом и сложном напряженном состоянии, реализуемом при растяжении-сжатии-кручении-внутреннем или внешнем давлении металлических трубчатых и сплошных образцов как подтверждают, так и отрицают ее справедливость даже при умеренных пластических деформациях (до 15 %) [2-12]. В теоретических исследованиях при описании процессов деформирования и разрушения материалов, пластические свойства которых чувствительны к изменению вида напряженного состояния, физические соотношения (условия пластичности, критерии разрушения) дополняются зависимостями от параметров и функций вида напряженного состояния, выраженных через главные напряжения [10, 13-20]. Для описания вида напряженного состояния широко применяется параметр «жесткости» напряженного состояния (параметр трехосности напряжений), являющийся отношением среднего напряжения (первого инварианта тензора напряжений) к интенсивности касательных напряжений (квадратный корень из второго инварианта тензора напряжений) [13-16, 20-25], . (1) Показано, что предельная нагрузка существенным образом зависит от параметров, характеризующих чувствительность пластических свойств материалов к изменению вида напряженного состояния [15]. В модели накопления повреждений [16, 17] для учета зависимости величины критического значения энергии пластического разрыхления , при которой происходит разрушение в точке материала, от реализуемого вида напряженного состояния применяется параметр . (2) В [10, 19] для учета чувствительности пластических свойств от вида напряженного состояния вводится зависимость функции, описывающей поверхность нагружения с учетом микронапряжений, от параметра вида напряженного состояния, пропорционального отношению третьего инварианта девиатора активных напряжений ко второму в степени 3/2, . Процесс образования или снятия микронапряжений также зависит от вида напряженного состояния и определяется через второй и третий инварианты девиатора микронапряжений [10, 19]. Поскольку вид напряженного состояния определяется отношением главных напряжений, а этих отношений два, то и параметров вида напряженного состояния в общем случае также два [26]. В ряде работ для учета вида напряженного состояния при построении моделей пластических сред с поврежденностью рассматривается набор параметров, характеризующих вид напряженного состояния (меры напряженного состояния): параметр трехосности напряжений, параметр Надаи-Лоде по напряжениям (3) и параметр [13, 22, 27, 28]. При экспериментальном изучении поведения материалов при сложном напряженном состоянии для реализации сложного напряженного состояния в образцах применяются надрезы, выточки [29, 23]. Анализ напряженно-деформированного состояния (НДС) сплошных и полых образцов с концентраторами напряжений в стадии упругопластического деформирования связан с большими трудностями в связи с его существенной неодноосностью и неоднородностью, вызванными локализацией деформаций [29]. Перспективным является применение численного моделирования, которое позволяет изучать свойства материалов, процессы деформирования и разрушения при достаточно сложной геометрии образцов с учетом геометрической и физической нелинейности, краевых эффектов, неоднородности НДС без привлечения упрощающих предположений и априорных гипотез силового и кинематического характера. В настоящей работе выполнены экспериментальные и численные исследования процессов упругопластического деформирования сплошных цилиндрических образцов с гладкой поверхностью и с круговым концентратором напряжений из стали 09Г2С при пропорциональном кинематическом нагружении кручением и растяжением до разрушения. Изучено влияние вида напряженного состояния на деформационное упрочнение стали 09Г2С при растяжении и кручении при больших деформациях, а также проведена оценка влияния комбинированного нагружения на потерю устойчивости пластического деформирования в виде шейки. Результаты подобных исследований при больших деформациях и неоднородном НДС практически отсутствуют. 1. Методика численного решения Полная система механических уравнений, описывающих решение обобщенных двумерных задач кручения, записывается в цилиндрической системе эйлеровых координат ( - ось вращения) в текущей лагранжевой постановке [30]. Кинематические соотношения формулируются в скоростях и строятся в метрике текущего состояния, что позволяет учитывать большие формоизменения. Компоненты тензоров скоростей деформаций (симметричной части градиента скорости перемещений) и скоростей вращения при учете равенства (- угол закручивания по окружной координате ) имеют вид , , , , , , (4) , Точка над переменной означает частную производную по времени, а символ после запятой - частную производную по соответствующей координате. Уравнение движения сплошной среды следует из уравнения баланса виртуальных мощностей (5) Здесь - область, занимаемая меридиональным сечением сплошной среды; - компоненты тензоров напряжений Коши; и - компоненты векторов скорости перемещения и ускорения перемещения; - плотность; - часть поверхности, на которой задается априори известная поверхностная нагрузка; - часть поверхности, на которой задаются контактные давления, определяемые в процессе решения; и - компоненты поверхностных и контактных нагрузок . В силу осевой симметрии все искомые функции зависят от радиальной и осевой координат и не зависят от окружной. Для устранения особенности на оси вращения и повышения точности численного решения уравнений движения вводятся новые функции , , . Тогда общее уравнение динамики (5) при учете соотношений (4) преобразуется к виду (6) При моделировании процессов активного нагружения, близких к пропорциональным, упругопластические свойства материалов описываются теорией течения с нелинейным изотропным упрочнением. Связь между компонентами девиатора скоростей напряжений и упругими составляющими компонент девиатора скоростей деформаций осуществляется на основе обобщенного закона Гука в метрике текущего состояния , ,, , (7) где - производная Яуманна, описывающая поворот частицы среды как жесткого целого; , - модули объемного сжатия и сдвига; - давление; - символ Кронекера. В силу осевой симметрии в выражениях для компонент тензоров скоростей вращения (последние три равенства в (4)) полный угол поворота относительно оси не учитывается. Скорости пластических деформаций определяются ассоциированным законом течения , (8) где - зависимость интенсивности напряжений от параметра Одквиста при монотонном нагружении (истинная диаграмма деформирования материала). Система уравнений (4), (6)-(8), дополненная кинематическими граничными и начальными условиями, решается методом конечных элементов в сочетании с явной схемой интегрирования типа «крест» [8]. Для численного моделирования процессов деформирования образцов при пропорциональном комбинированном нагружении растяжением-кручением реализована зависимость диаграммы деформирования от вида напряженного состояния. Уточненная диаграмма деформирования является линейной комбинацией диаграмм растяжения и кручения, коэффициенты которой зависят от параметра вида напряженного состояния П (2), принятого в виде [18], параметра «жесткости» напряженного состояния (1) или параметра Надаи-Лоде по напряжениям (3) (9) (10) (11) Параметр вида напряженного состояния изменяется в отрезке , при этом величина соответствует всестороннему сжатию, - одноосному сжатию, - чистому сдвигу, - одноосному растяжению, - всестороннему растяжению [18]. Истинная диаграмма деформирования принимается в виде (12) Здесь - диаграмма деформирования при одноосном растяжении; - диаграмма деформирования при кручении. Параметр «жесткости» напряженного состояния изменяется в интервале , при этом величина соответствует всестороннему сжатию, - одноосному сжатию, - чистому сдвигу, - одноосному растяжению, - всестороннему растяжению [20]. Истинная диаграмма деформирования принимается в виде (13) Параметр Надаи-Лоде по напряжениям изменяется в отрезке , при этом величина соответствует одноосному растяжению, - чистому сдвигу, - одноосному сжатию. Истинная диаграмма деформирования принимается в виде (14) 2. Монотонное пропорциональное нагружение кручением и растяжением цилиндрических образцов Проводилось численное и экспериментальное исследование процессов деформирования осесимметричных образцов переменной толщины с цилиндрической рабочей частью при монотонном кинематическом нагружении растяжением-кручением с учетом больших деформаций и неоднородности НДС. Эксперименты выполнялись на модернизированном испытательном комплексе Z100 ZWICK-ROEL, позволяющем проводить испытания на растяжение и кручение цилиндрических образцов при одновременном синхронизированном по времени задании параметров: скорость изменения продольного перемещения - скорость изменения угла закручивания. Предельные значения параметров испытательной техники: по силе ±100 кН, по крутящему моменту ±1000 Нм. Измерительная аппаратура комплекса включает в себя лазерный экстензометр Laser Xtens, позволяющий одновременно измерять как продольные деформации (база измерения 20-100 мм) с разрешающей способностью ≤ 1 мкм, так и деформации при кручении, класс точности 1 по ISO 9513, датчик силы 0-250 кН, класс точности 1 по ISO 7500-1. Программное обеспечение testXpert настраивается для проведения испытаний при растяжении, кручении или комбинированном нагружении и предоставляет возможность получения необходимых экспериментальных значений параметров исследуемого процесса (осевое перемещение, угол закручивания, осевая сила, крутящий момент) в параметрическом или других видах. Геометрические параметры образца (рис. 1): радиус и длина рабочей части м, , размеры переходной части , , радиус и длина цилиндрического захвата , , общая длина образца м. Образец изготовлен из стали 09Г2С, упругие характеристики материала: модуль объемного сжатия МПа, модуль сдвига МПа, предел текучести МПа, плотность кг/м3. Рис. 1. Конструктивная схема образца с цилиндрической рабочей частью Fig. 1. The scheme of the specimen with a cylindrical working part Для оценки взаимного влияния двух видов нагружения на процесс деформирования вводятся следующие параметры: , , , , , , , Кинематический параметр характеризует соотношение кручения и растяжения, , и - условные интенсивность деформаций, осевая и сдвиговая деформации соответственно на поверхности рабочей части образца, отнесенные к начальным длине и радиусу рабочей части образца, и - осевое перемещение и угол закручивания между торцами. При чистом кручении , , при одноосном растяжении , , при комбинированном нагружении ; - осевая сила, - крутящий момент, А - суммарная работа осевой силы и крутящего момента. При выражение для работы можно преобразовать к виду ( - обобщенная сила) 73 В качестве критерия потери устойчивости пластического деформирования с образованием шейки принимается условие , характеризующее момент достижения обобщенной силой максимальных значений. Данное условие является обобщением известного условия А. Консидера, полученного при растяжении стрежней [3]. В эксперименте реализованы следующие программы нагружения: одноосное растяжение (), чистое кручение (), пропорциональное нагружение растяжением-кручением при На рис. 2 и 3 линиями с маркерами приведены экспериментальные зависимости осевой силы от условной осевой деформации и крутящего момента от условной сдвиговой деформации . Соответствующие значения параметра приведены рядом с кривыми. Точками на кривых отмечены моменты потери устойчивости пластического деформирования с образованием шейки. Рис. 2. Экспериментальные и расчетные зависимости осевой силы от условной осевой деформации при растяжении-кручении образцов с цилиндрической рабочей частью Fig. 2. Experimental and calculated dependencies of the axial force from the conditional axial strain under tension-torsion of specimens with a cylindrical working part Рис. 3. Экспериментальные и расчетные зависимости крутящего момента от условной сдвиговой деформации на поверхности рабочей части образца при растяжении-кручении образцов с цилиндрической рабочей частью Fig. 3. Experimental and calculated dependencies of the torsional moment from the conditional shear strain under tension-torsion of specimens with a cylindrical working part 74 На основе проведенных экспериментов строились истинные диаграммы деформирования для стали 09Г2С. На рис. 4 сплошной линией отмечена истинная диаграмма деформирования, построенная на основе экспериментальных данных по монотонному кручению () до разрушения образца по методике [8], пунктирной - по одноосному растяжению () с учетом образования шейки по методике [31]. На диаграмме растяжения стали 09Г2С имеется площадка текучести, предел текучести при растяжении МПа. Материальные кривые, построенные при этих двух видах нагружения, при величине параметра Одквиста практически совпадают. При величине параметра Одквиста (момент начала образования шейки при одноосном растяжении) расхождение между кривыми составляет 4,5 % при - 17 %. В работе [12] отмечается, что диаграммы деформирования при растяжении и кручении принципиально не могут совпадать, поскольку при растяжении до начала образования шейки реализуется простое нагружение, а при кручении - сложное с поворотом осей тензора напряжений. Однако на основе имеющихся у авторов экспериментальных данных [11, 12] был сделан вывод, что расхождения эти лежат в рамках погрешности эксперимента. Для ряда материалов (стали 3, ЭЯ-1Т, ШХ-15, 12Х80Н10Т, латунь [8, 11, 12]) «гипотеза единой кривой» подтверждается экспериментально при больших деформациях до 80 %. Причиной различия диаграмм растяжения и кручения стали 09Г2С, по-видимому, является крупнозернистая неоднородная структура исследуемого материала. Рис. 4. Истинные диаграммы деформирования стали 09Г2С при растяжении (пунктир) и кручении (сплошная линия) Fig. 4. The true stress-strain diagrams of steel 09G2S under tension (dotted line) and torsion (solid line) Проведено численное моделирование процессов деформирования образцов с учетом зависимости диаграммы деформирования от вида напряженного состояния. В расчетах темп нагружения задавали таким образом, чтобы вклад инерционных сил был пренебрежимо мал. На рис. 2, 3 сплошными линиями приведены результаты расчетов, в которых зависимость диаграмм деформирования от вида напряженного состояния учитывается выражениями (9), (10) с параметром вида напряженного состояния П и параметром жесткости (эти расчеты практически совпадают), пунктирными - выражением (11) с параметром Надаи-Лоде по напряжениям. При комбинированном пропорциональном нагружении расчетные кривые, учитывающие вид НДС, располагаются между кривыми, полученными с диаграммами деформирования растяжения или кручения. Закритическое поведение образцов после локализации деформаций с образованием шейки лучше описывается моделью, учитывающей вид напряженного состояния с помощью параметра П (или ). Неучет вида напряженного состояния приводит к существенным различиям в определении моментов потери устойчивости пластического деформирования в виде шейки и в описании закритического поведения. При комбинированном нагружении растяжением-кручением () в эксперименте и расчете после момента потери устойчивости в процессе образования шейки наблюдается более интенсивный рост осевой силы до разрушения, при этом величина крутящего момента уменьшается из-за уменьшения радиуса поперечного сечения шейки, что согласуется с результатами исследований для стали 12Х18Н10Т [8, 31]. Процессы деформирования до момента потери устойчивости близки к лучевым, а закритическое поведение характеризуется траекториями малой кривизны. При комбинированном нагружении () до момента потери устойчивости во всем объеме образца вид напряженного состояния медленно изменяется в сторону растяжения, так как радиус образца уменьшается. Изменение вида напряженного состояния в шейке в сторону деформации чистого сдвига после момента потери устойчивости происходит наиболее интенсивно при чистом растяжении (). До момента падения осевой силы во всем объеме образца происходит процесс активного нагружения. В месте образования шейки активный процесс продолжается до разрушения. В экспериментах шейка образуется в средней части образца, ее положение зависит от многих случайных факторов. В расчетах образование шейки происходит на плоскости симметрии. 3. Монотонное пропорциональное нагружение кручением и растяжением цилиндрических образцов с концентратором напряжений Проводилось численное и экспериментальное исследование процессов деформирования осесимметричных образцов переменной толщины с кольцевыми концентраторами напряжений на цилиндрической рабочей части из стали 09Г2С при монотонном кинематическом нагружении растяжением-кручением с учетом больших деформаций и неоднородности НДС. Геометрические параметры образца (рис. 5) и упругопластические характеристики материала такие же, как у образцов с цилиндрической рабочей частью. Радиус кольцевого концентратора . В эксперименте реализованы следующие программы нагружения: одноосное растяжение (), чистое кручение (), пропорциональное нагружение растяжением-кручением при . На рис. 6, 7 приведены зависимости осевой силы от условной осевой деформации и крутящего момента от условной сдвиговой деформации на поверхности цилиндрической части образца . Рис. 5. Расчетная схема образца с круговым концентратором на цилиндрической рабочей части Fig. 5. The computational scheme of the specimen with a circular concentrator on a cylindrical working part Рис. 6. Экспериментальные и расчетные зависимости осевой силы от условной осевой деформации при растяжении- кручении образцов с круговым концентратором Fig. 6. Experimental and calculated dependencies of the axial force from the conditional axial strain under tension-torsion of specimens with a circular concentrator on a cylindrical working part Рис. 7. Экспериментальные и расчетные зависимости крутящего момента от условной сдвиговой деформации на цилиндрической поверхности при растяжении-кручении образцов с круговым концентратором Fig. 7. Experimental and calculated dependencies of the torsional moment from the conditional shear strain under tension-torsion of specimens with a circular concentrator on a cylindrical working part 75 Обозначения такие же, как на рис. 2, 3. Соответствующие значения параметра q приведены рядом с кривыми. Расчетные зависимости хорошо согласуются с экспериментальными данными как по величине, так и по моменту достижения максимальных значений интегральных характеристик. При чистом растяжении нисходящая ветвь экспериментальной кривой расчетом описывается также удовлетворительно. При наложении кручения на растяжение и чистом кручении эта ветвь не описывается, поскольку при численном моделировании не были учтены поврежденность материала и разрушение, которое имеет многостадийный характер, включающий образование микроповреждений, слияние микропор, а затем развитие макроразрушения. Процесс макроразрушения начинается на поверхности образца и при дальнейшем нагружении (кручении) распространяется к оси вращения. Заключительный этап деформирования лабораторных образцов, характеризующийся разупрочнением и разрушением материала, данной методикой описывается только до момента начала процесса макроразрушения. 4. Анализ НДС-образцов при больших деформациях По результатам расчетов проведен анализ НДС-образцов в момент потери устойчивости пластического деформирования с образованием шейки (только для образцов с гладкой поверхностью) и в момент начала разрушения при растяжении, комбинированном нагружении (для образцов с гладкой поверхностью , для образцов с круговым концентратором напряжений ) и кручении. Получено распределение параметра Одквиста (рис. 8, 9) в меридиональном сечении образцов. На рис. 8 приведены данные для образцов с цилиндрической рабочей частью, буквами a и b обозначены распределения полей при растяжении, c и d - при комбинированном нагружении, e - при кручении, a и c - в момент потери устойчивости, b, d, e - в момент начала разрушения. На рис. 9 приведены данные для образцов с круговой выточкой в момент начала разрушения, буквой a обозначено распределение полей при растяжении, b - при комбинированном нагружении, c - при кручении. На рис. 10-19 приведены расчетные распределения параметров НДС вдоль радиуса на плоскости симметрии образцов. По оси абсцисс отложены номера ячеек конечно-элементной сетки, нумерация осуществляется от свободной поверхности к оси вращения. Сплошными линиями обозначены данные при растяжении, пунктирными - при комбинированном нагружении, штрихпунктирными - при кручении. Кривые без маркеров соответствуют моментам начала разрушения, с маркерами - моментам потери устойчивости пластического деформирования с образованием шейки (только для сплошных образцов). а б в г д Рис. 8. Распределение параметра Одквиста в меридиональном сечении образцов в момент потери устойчивости пластического деформирования с образованием шейки (a, в) и в момент начала разрушения (б, г, д) при растяжении (a, б), растяжении-кручении (в, г) и кручении (д) Fig. 8. Distribution of the Odquist parameter in the meridional section of specimens at the moment of the stability loss under plastic deformation with the neck formation (a, c) and in the beginning of fracture (b, d, e) under tension (a, b), tension-torsion (c, d) and torsion (e) а б в Рис. 9. Распределение параметра Одквиста в меридиональном сечении образцов с концентратором напряжений в момент начала разрушения Fig. 9. Distribution of the Odquist parameter in the meridional section of the specimens with a stress concentrator at the beginning of fracture 76 На рис. 10, 11 приведены распределения интенсивности напряжений и параметра Одквиста в шейке вдоль радиуса сплошного образца, на рис. 12, 13 - на плоскости симметрии образца с кольцевым концентратором. При растяжении сплошного образца до момента потери устойчивости пластического деформирования с образованием шейки наблюдается однородность НДС рабочей части как вдоль радиуса, так и вдоль оси вращения. При кручении и комбинированном нагружении до образования шейки наблюдается однородность НДС вдоль оси вращения и существенная неоднородность вдоль радиуса. При растяжении в процессе образования шейки максимальные значения интенсивности напряжений и параметра Одквиста наблюдаются в шейке на оси вращения и вдоль радиуса отличаются незначительно (см. рис. 10, 11). При кручении и кручении-растяжении максимальные значения этих функций достигаются на поверхности образца, зависимость параметра Одквиста от радиуса близка к линейной. При растяжении и/или кручении образца с круговым концентратором напряжений максимальные значения интенсивности напряжений и параметра Одквиста возникают на поверхности выточки в поперечном сечении с минимальным радиусом (см. рис. 12, 13). Рис. 10. Распределение интенсивности напряжений в шейке вдоль радиуса образца с цилиндрической рабочей частью Fig. 10. Distribution of stress intensity in the neck along the radius of the specimen with a cylindrical working part Рис. 11. Распределение параметра Одквиста в шейке вдоль радиуса образца с цилиндрической рабочей частью Fig. 11. Distribution of the Odquist parameter in the neck along the radius of the specimen with a cylindrical working part На рис. 14, 15 приведены расчетные распределения по радиусу (в шейке) образца с цилиндрической рабочей частью параметров вида напряженного состояния П и трехосности напряжений в момент потери устойчивости пластического деформирования с образованием шейки и в момент начала разрушения. При растяжении и растяжении-кручении сплошного образца в момент потери устойчивости пластического деформирования с образованием шейки значения обоих параметров практически совпадают. При растяжении их величины близки к 1, при комбинированном нагружении изменяются монотонно от 0,15 на поверхности до 1,1 на оси вращения. В процессе образования шейки их величины на поверхности образца практически не меняются, а на оси вращения становятся больше 1, что соответствует смещению вида напряженного состояния от одноосного растяжения в сторону всестороннего сжатия (см. рис. 14, 15). Рис. 12. Распределение интенсивности напряжений вдоль радиуса на плоскости симметрии образца с круговым концентратором Fig. 12. Distribution of stress intensity along the radius on the plane of symmetry of a specimen with a circular concentrator Рис. 13. Распределение параметра Одквиста вдоль радиуса на плоскости симметрии образца с круговым концентратором Fig. 13. Distribution of Odquist parameter along the radius on the plane of symmetry of a specimen with a circular concentrator На рис. 16 приведены расчетные распределения параметра Надаи-Лоде по напряжениям . В процессе образования шейки имеется немонотонная зависимость параметра Надаи-Лоде от условной интенсивности деформаций (параметра нагружения), поэтому распределения в момент потери устойчивости и в момент разрушения практически совпадают. N Рис. 14. Распределение параметра вида напряженного состояния П в шейке вдоль радиуса образца с цилиндрической рабочей частью Fig. 14. Distribution of the parameter of the form of the stressed state П in the neck along the radius of the specimen with a cylindrical working part N -χσ Рис. 15. Распределение параметра трехосности напряжений в шейке вдоль радиуса образца с цилиндрической рабочей частью Fig. 15. Distribution of the parameter of stress triaxiality in the neck along the radius of a specimen with a cylindrical working part N Рис. 16. Распределение параметра Надаи-Лоде по напряжениям в шейке вдоль радиуса образца с цилиндрической рабочей частью Fig. 16. Distribution of the Nadai-Lode parameter for stresses in the neck along the radius of the specimen with a cylindrical working part На рис. 17-19 представлены кривые на плоскости симметрии образцов с концентратором в момент начала разрушения. При растяжении и комбинированном нагружении параметры вида напряженного состояния П и трехосности напряжений на поверхности принимают значения, близкие к 1, а на оси вращения - большие 1. N ξ Рис. 17. Распределение параметра вида напряженного состояния П вдоль радиуса на плоскости симметрии образца с круговым концентратором Fig. 17. Distribution of the parameter of the form of the stressed state П along the radius in the plane of symmetry of the specimen with a circular concentrator N -χσ Рис. 18. Распределение параметра трехосности напряжений вдоль радиуса на плоскости симметрии образца с круговым концентратором Fig. 18. Distribution of the three-axis stress parameter along the radius on the plane of symmetry of the specimen with a circular concentrator N Рис. 19. Распределение параметра Надаи-Лоде по напряжениям вдоль радиуса на плоскости симметрии образца с круговым концентратором Fig. 19. Distribution of the Nadai-Lode parameter for stresses along the radius on the plane of symmetry of the specimen with a circular concentrator 5. Анализ разрушения образцов при различных видах НДС На рис. 20 и 21 приведены фотографии образцов после испытаний до разрушения при монотонном нагружении растяжением (1), кручением (2) и кручением-растяжениием (3) с для сплошного образца и с для образца с концентратором напряжений. Видно, что при кручении сплошного образца (см. рис. 20) не происходит потери устойчивости пластического деформирования, как при растяжении с образованием шейки, а разрушение происходит от сдвиговых деформаций в плоскости, перпендикулярной оси стержня. При совместном действии кручения-растяжения образуется шейка, радиус поперечного сечения которой больше, а длина меньше, чем при чистом растяжении; разрушение происходит по винтовой поверхности с максимальными сдвиговыми деформациями в зависимости от соотношения деформаций растяжения и кручения. Рис. 20. Образцы с цилиндрической рабочей частью после испытаний растяжением (1), кручением (2) и растяжением-кручением (3) Fig. 20. Specimens with a cylindrical working part after testing by tension (1), torsion (2) and tension-torsion (3) Рис. 21. Образцы с кольцевым концентратором напряжений после испытаний растяжением (1), кручением (2) и растяжением-кручением (3) Fig. 21. Specimens with a circular stress concentrator after testing by tension (1), torsion (2), and tension-torsion (3) В аналогичных экспериментах для стали 12Х18Н10Т при чистом растяжении в шейке в месте разрушения образца образуется симметричная чашка, так как сначала происходит хрупкое разрушение отрывом поперек оси вращения, а затем сдвиговое с выходом на поверхность образца [32]. В случае стали 09Г2С происходит сдвиговое разрушение под углом, близким к 45° к оси стержня. Поверхность сдвигового разрушения шероховатая, как при хрупком отрыве стали 12Х18Н10Т. Причиной различия в характере разрушения является крупнозернистая структура стали 09Г2С. При чистом кручении вид поверхности разрушения практически совпадает для обеих сталей. На образцах с кольцевым концентратором характер разрушения аналогичен характеру разрушения образцов без выточки. Заключение Построенные истинные диаграммы деформирования стали 09Г2С на основе экспериментов по растяжению и кручению цилиндрических сплошных образцов до разрушения практически совпадают до величины параметра Одквиста 15 % и существенно отличаются при больших деформациях. В результате численного моделирования процессов деформирования образцов из стали 09Г2С при монотонном кинематическом нагружении растяжением-кручением до разрушения получено хорошее совпадение с экспериментальными данными по интегральным характеристикам (усилиям и моментам) для образцов с гладкой поверхностью и образцов с кольцевым концентратором. Расчеты проведены с учетом влияния вида напряженного состояния на диаграмму деформирования. Неучет вида напряженного состояния приводит к существенным различиям в определении моментов образования шеек и в описании закритического поведения. В рассмотренных задачах растяжения-кручения до момента потери устойчивости пластического деформирования с образованием шейки выбор параметра вида напряженного состояния не оказывает существенного влияния на результаты расчетов. Закритическое поведение образцов после момента потери устойчивости пластического деформирования с образованием шейки лучше описывается моделями с параметром П или .About the authors
V G Bazhenov
Researcher Institute of Mechanics, National Research Lobachevsky State University of Nizhni Novgorod
D A Kazakov
Researcher Institute of Mechanics, National Research Lobachevsky State University of Nizhni Novgorod
E V Nagornykh
Researcher Institute of Mechanics, National Research Lobachevsky State University of Nizhni Novgorod
D L Osetrov
Researcher Institute of Mechanics, National Research Lobachevsky State University of Nizhni Novgorod
A A Ryabov
FGUP RFYAC-VNIIEF
References
- Дегтярев В.П. Деформации и разрушение в высокона-пряженных конструкциях. - М.: Машиностроение, 1987. - 105 с.
- Зубчанинов В.Г. Механика сплошных деформируемых сред. - Тверь: Изд-во ТГТУ: ЧуДо, 2000. - 703 с.
- Малинин Н.Н. Прикладная теория пластичности и пол-зучести. - М.: Машиностроение, 1975. - 400 с.
- Механические свойства материалов при статическом нагружении трубчатых образцов в условиях плоского и объ-емного напряженных состояний / А.Ю. Кузькин, Д.В. Латышев, М.Ю. Петров, В.А. Попов // Науч.-техн. ведомости C.-Петерб. гос. политехн. ун-та. - 2014. - № 2 (195). - С. 162-173.
- Шлянников В.Н., Иштыряков И.С., Яруллин Р.Р. Ха-рактеристики деформирования сплава Д16Т при совместном нагружении растяжением, сжатием, кручением и внутренним давлением // Тр. Академэнерго. - 2014. - № 3. - С. 78-90.
- Ипатова А.В., Вильдеман В.Э. Построение материаль-ных функций неупругого деформирования алюминиевого спла-ва Д16Т по результатам испытаний на растяжение и кручение // Вестн. Самар. гос. техн. ун-та. Сер. Физ.-мат. науки. - 2012. - № 4 (29). - С. 106-114. DOI: https://DOI.org/10.14498/vsgtu1106
- Высокоскоростное деформирование и разрушение стали 09Г2С / Вл.Вас. Баландин, Вл.Вл. Баландин, А.М. Брагов, Л.А. Игумнов, А.Ю. Константинов, А.К. Ломунов // Изв. РАН. МТТ. - 2014. - № 6. - С. 78-85.
- Моделирование процессов деформирования и локали-зации пластических деформаций при кручении-растяжении тел вращения / В.Г. Баженов, С.В. Зефиров, Л.Н. Крамарев, Е.В. Павленкова // ПММ. - 2008. - Т. 72. - Вып. 2. - С. 342-350.
- Бердин В.К., Кашаев Р.М. Об определении напряжен-ного состояния при растяжении с кручением сплошного ци-линдра // Проблемы прочности. - 2001. - № 1. - С. 28-37.
- Бондарь В.С., Даншин В.В. Пластичность. Непро-порциональные нагружения. - М.: Изд-во МГТУ «МАМИ», 2008. - 218 с.
- Дощинский Г.А., Максак В.И. О деформации простой сдвиг // Изв. Том. политехн. ин-та. - 1970. - Т. 173. - С. 13-17.
- Максак В.И., Дощинский Г.А. Методика и исследо-вание больших пластических деформаций при простом нагру-жении // Изв. Том. политехн. ин-та. - 1970. - Т. 173. - С. 3-9.
- Шлянников В.Н., Иштыряков И.С. Параметры функ-ций вида напряженного состояния для алюминиевого сплава Д16Т // Тр. Академэнерго. - 2014. - № 4. - С. 51-63.
- Ломакин Е.В., Мельников А.М. Задачи плоского на-пряженного состояния тел с вырезами, пластические свойства которых зависят от вида напряженного состояния // Изв. РАН. МТТ. - 2011. - № 1. - С. 77-89.
- Ломакин Е.В., Мельников А.М. Пластическое плос-кое напряженное состояние тел, свойства которых зависят от вида напряженного состояния // Вычислительная механика сплошных сред. - 2009. - Т. 2, № 2. - С. 48-64.
- Ломакин Е.В. Зависимость предельного состояния композитных и полимерных материалов от вида напряженного состояния // Мех. композит. материалов. - 1988. - № 1. - С. 3-9.
- Соотношения модели поврежденной среды для мате-риалов, подвергающихся терморадиационным воздействиям / С.А. Ка¬пустин, В.А. Горохов, О.Ю. Виленский, В.Б. Кайдалов, А.А. Руин // Проблемы прочности и пластичности. - 2012. - № 74. - С. 5-15.
- Капустин С.А., Чурилов Ю.А., Горохов В.А. Моде-лирование нелинейного деформирования и разрушения конст-рукций в условиях многофакторных воздействий на основе МКЭ: монография. - Н. Новгород: Изд-во Нижегород. гос. ун-та им. Н.И. Лобачевского, 2015. - 347 с.
- Бондарь В.С., Абашев Д.Р. Пластическое деформи-рование материалов, чувствительных к виду напряженного состояния // Вестник Пермского национального исследова-тельского политехнического университета. Механика. - 2018. - № 1. - С. 29-39. doi: 10.15593/perm.mech/2018.1.03.
- Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругопластических сред с повреждениями. - М.: Физ-матлит, 2008. - 424 с.
- Tekoğlu C, Hutchinson JW, Pardoen T. 2015 On locali-zation and void coalescence as a precursor to ductile fracture. Phil. Trans. R. Soc. A 373. doi: 10.1098/rsta.2014.0121.
- Nahshon K., Hutchinson J.W. Modification of the Gur¬son Model for shear failure // European Journal of Mechanics A/Solids. - 2008. - Vol. 27. - P. 1-17. doi: 10.1016/j.euromechsol.2007.08.002
- Kossakowski P.G. Experimental determination of the void volume fraction for s235jr steel at failure in the range of high stress triaxialities // Arch. Metall. Mater. - 2017. - Vol. 62. - No. 1. - P. 167-172. doi: 10.1515/amm-2017-0023
- Modeling of Ductile Failure in High Strength Steel / N. Her¬zig, S. Abdel-Malek, L.W. Meyera, S.J. Cimpoeru // Procedia Engineering. 2017. - Vol. 197. - Р. 285-293. doi: 10.1016/j.proeng.2017.08.106
- Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities / X Chen, Y Peng, S Peng, S Yao, C Chen, P Xu // PLoS ONE. - 2017. - Vol. 12(7): e0181983. https://DOI.org/10.1371/journal.pone.0181983
- Ломакин Е.В. Механика сред с зависящими от вида напряженного состояния свойствами // Физ. мезомеханика. - 2007. - № 10 (5). - С. 41-52.
- Calibration procedures for a computational model of ductile fracture / Z. Xue, M.G. Pontin, F.W. Zok, John W. Hutchinson // Engineering Fracture Mechanics. - 2009. - Vol. 77(3). - P. 492-509. doi: 10.1016/j.engfracmech.2009.10.007
- Nielsen K.L., John W. Hutchinson. 2012. Cohesive trac-tion-separation laws for tearing of ductile metal plates // Interna-tional Journal of Impact Engineering. - 2011. - Vol. 48. - P. 15-23. doi: 10.1016/j.ijimpeng.2011.02.009
- Писаренко Г.С., Лебедев А.А. Деформирование и проч¬ность материалов при сложном напряженном состоянии. - Киев: Наукова думка, 1976. - 416 с.
- Коробейников С.Н. Нелинейное деформирование твердых тел. - Новосибирск: Изд-во СО РАН, 2000. - 262 с.
- Экспериментально-теоретическое исследование пре-дельных состояний упругопластических стержней различного поперечного сечения при растяжении / В.Г. Баженов, А.И. Кибец, П.В. Лаптев, С.Л. Осетров // Проблемы механики: сб. ст. к 90-летию со дня рождения А.И. Ишлинского / под ред. Д.М. Климова и др. - М.: Физматлит, 2003. - С. 116-123.
- Баженов В.Г., Жегалов Д.В., Павленкова Е.В. Чис-ленное и экспериментальное исследование упругопластических процессов растяжения-кручения осесимметричных тел при больших деформациях // Изв. РАН. МТТ. - 2011. - № 2. - С. 57-66.