УДК 629.539.3
А.И. Уваров
ОПРЕДЕЛЕНИЕ ДЛИНЫ ШЕЙКИ В РАСТЯГИВАЕМОМ КРУГЛОМ МЕТАЛЛИЧЕСКОМ ОБРАЗЦЕ ДЛЯ РАСЧЕТНОЙ ОЦЕНКИ УДЕЛЬНОЙ ЭНЕРГОЕМКОСТИ МАТЕРИАЛА
Нижегородский государственный технический университет им. Р.Е. Алексеева
С целью расчета энергоемкости материалов, используемых в элементах конструкции автомобиля и обеспечивающих его пассивную безопасность, разработана аналитическая модель начала процесса локального деформирования (образования шейки) в растягиваемом цилиндрическом пластическом теле при оси симметричной деформации. Использован кинематический метод теории предельного равновесия. Получено приближенное аналитическое решение, определяющее локальное поле скоростей перемещений, соответствующее такой же величине растягивающей нагрузки, как и при равномерном деформировании. При этом размер района локального деформирования вдоль оси симметрии пропорционален радиусу цилиндра. Произведено качественное (по форме недеформируемых областей) и количественное (по минимальному диаметру шейки) сопоставление полученного аналитического решения с результатами численного моделирования. Получено удовлетворительное совпадение. Результаты могут быть использованы для расчетной оценки удельной энергоемкости материалов энергопоглощающих элементов конструкций пассивной безопасности.
Ключевые слова: пассивная безопасность, шейка в круглом образце, удельная энергоемкость.
Введение
В современных транспортных системах в качестве элементов систем пассивной безопасности применяются конструкции, поглощающие энергию соударений движущихся тел за счет собственных деформаций [1-3]. В автомобилестроении специальными элементами, выполняющими функцию снижения ударных нагрузок на пассажиров и груз за счет собственного деформирования при ударах, являются бамперы и ударные амортизаторы [4, 5]. Принцип поглощения энергии движения используется также при создании конструктивной защиты атомных энергетических установок от ударов движущихся объектов [6]. При проектировании энергопоглощающего элемента конструкции необходимо выбрать оптимальный материал. Критерием эффективности материала, используемого для изготовления энергопо-глощающих элементов, может служить энергия, поглощаемая единицей объема материала, которая равна площади диаграммы напряжения - относительное удлинение, получаемое при механических испытаниях образцов на растяжение [7].
Разрушению растягиваемого образца из металлических сплавов предшествует образование района локального деформирования (шейки). Энергию, поглощаемую материалом образца, можно разделить на энергию равномерного деформирования и энергию деформации в шейке. Доля энергии, поглощаемой в шейке, зависит от соотношения размеров образца, формы образца, и при определении энергопоглощающей способности материалов ее целесообразно отделить и учитывать площадь диаграммы при равномерном деформировании до образования шейки, когда деформации одинаковы во всем объеме материала (заштрихованная площадь на рис. 1). Площадь диаграммы напряжений, ограниченная величиной относительного равномерного удлинения öp, не зависит от длины образца и равна av - энергии, поглощаемой единицей объема материала деформируемого образца до образования шейки (штриховые линии на рис. 1). Приближенно величина av может быть вычислена (без учета упругих деформаций), как площадь трапеции:
aV = 1(Gr В ) '5 р (1)
© Уваров А.И., 2019.
ав
ат
О
Рис. 1. Схема диаграммы напряжений при растяжении металлического образца:
а - нормальное напряжение в поперечном сечении образца, 8 - относительное удлинение, аТ - предел текучести, аВ - временное сопротивление
8
Величина равномерного остаточного удлинения ёр не всегда содержится в данных сертификата на материал, который, как правило, содержит величины остаточного удлинения образца вместе с шейкой 5о и относительного поперечного сужения образца при разрыве у. При известной форме шейки величина ёр может быть вычислена. Основными геометрическими параметрами шейки являются ее длина и наименьший диаметр. Эти величины определяются путем испытаний образца на растяжение с соответствующими измерениями. Разработке метода расчетного определения этих размеров посвящена настоящая работа.
Экспериментальное исследование возможности прогнозирования места разрушения при пластическом растяжении проведено в работе [8]. В работе [9] исследовано влияние структуры кристаллической решетки на локализацию пластической деформации. Известно аналитическое решение задачи о напряженном состоянии в шейке круглого образца [10]. Исходными данными для расчетов при этом служат геометрические параметры: диаметр шейки и радиус кривизны ее контура, получаемые в эксперименте, то есть расчету должны предшествовать испытания образцов. В работе [11] разработана аналитическая модель определения начального размера района локального деформирования при растяжении широкой металлической пластины при плоской деформации.
Целью настоящей работы является расчетное определение размера района локальных деформаций (шейки) в начальный период ее возникновения в образце цилиндрической формы при осесимметричной деформации. В качестве метода исследования выбрана разработка аналитической модели локализации деформаций и сопоставление аналитической модели с результатами численного моделирования.
В первом разделе рассмотрена осесимметричная деформация при растяжении круглого образца. Получено аналитическое решение, аналогичное [11], определяющее длину шейки в начале ее образования. Выполнен анализ полученного приближенного решения. Во втором разделе содержится обоснование локализации деформаций. В третьем - приведены сопоставления, полученного аналитического решения с результатами численного моделирования и экспериментальными фактами.
Построение аналитической модели локального деформирования
Рассмотрим тело в форме прямого кругового цилиндра радиуса R (диаметра Do) длины L из идеально-жесткопластического материала с пределом текучести от, нагруженное растягивающей силой величины Р, направленной вдоль оси цилиндра. Найдем значение
нагрузки Р, соответствующее началу пластической деформации (предельную нагрузку). Используем кинематический метод теории предельного равновесия [12]. Ось Z цилиндрической системы координат г2в совпадает с осью цилиндра, и, м> - перемещения точек материала в направлениях осей г, Z соответственно.
Будем полагать, что имеет место осесимметричная деформация. При этом считаем, что перемещения точек материала тела, имевших одинаковые координаты г, Z до деформации, получают одинаковые перемещения и, ж, а перемещения в направлении перпендикулярном осям г, Z (в окружном направлении) отсутствуют.
На приведенных ниже рисунках показаны сечения плоскостью, проходящей через ось цилиндра. Вначале рассмотрим простейший случай равномерного деформирования, с которым сопоставим разработанную далее модель локального деформирования.
Равномерное деформирование
Пусть все элементы тела испытывают одинаковое удлинение в направлении действия нагрузки (вдоль оси Z )во всем объеме.
Пусть 5 - бесконечно малое относительное перемещение краев тела при растяжении силой Р (рис. 2).
р 1 "г р
-----------------t
< 1 ->
< Ь
рч ^
Рис. 2. Схема равномерного деформирования
Относительное удлинение в направлении оси Z во всех точках тела равно: ег = е = — • Из условия сохранения объема следует: 8 = ей = -— •
г В 2
Предположим, что все элементы тела находятся в линейном напряженном состоянии. При этом условию текучести Треска-Сен-Венана будут соответствовать напряжения = °г = °е = О-
Диссипация энергии в единице объема тела:
, §
а = е ■ ст,, = — стт т ь т
Диссипация энергии в объеме тела:
с
Б = а ■ (5 • I) = - • стг • (5 • X) = 5 • стг -5,
где 5 = %Я2 - площадь поперечного сечения тела. Работа внешней нагрузки:
А = Р-5
Дифференцируя по параметру 5 и приравнивая мощность внешней нагрузки мощности диссипации, находим величину предельной нагрузки для равномерного деформирования:
К* = Рь = (2)
Поскольку удовлетворяется условие текучести, и все элементы находятся в равновесии, то полученное решение можно считать точным.
Локальное деформирование. Общее описание модели
Предположим, что деформации происходят только в части тела длины I, которую в дальнейшем будем называть шейкой. Материал вне шейки остается недеформируемым (рис. 3).
р Р
1 <-> L
\ у
Рис. 3. Расположение локального района деформирования
Предположим, что при локализации деформаций (образовании шейки) деформации распределяются неравномерно по деформируемому объему, уменьшаясь от оси симметрии к наружным поверхностям и от середины шейки к ее краям. Данное предположение основано на результатах конечно элементного моделирования и известных экспериментальных данных об образовании трещины внутри шейки растягиваемого образца [12]. Ему соответствует рассмотренный далее однопараметрический кинематический механизм, который определяет различную деформацию в разных областях объема шейки (единственным параметром, определяющим поле перемещений, является длина деформируемого района I). Деформируемые области отделены друг от друга и от недеформируемого района коническими поверхностями.
Пусть локальный объем деформирования состоит из трех областей: одной внутренней (средней) и двух одинаковых внешних (крайних) областей, примыкающих к недеформируе-мым внешним областям тела. Внутренняя (средняя) область деформирования испытывает вдвое большее постоянное удлинение в направлении векторов нагрузки (кусочно-постоянная аппроксимация), чем внешние области деформирования (рис. 4).
Внутренняя область деформирования испытывает линейную деформацию величины 8 в направлении действия нагрузки. Данная часть имеет форму двух конусов радиуса R и высоты 1/4 . Ее объем (двух конусов) равен:
V =1 %Д2 - 2 =1 %Д 21. е 3 4 6
Две внешние области деформирования испытывает линейную деформацию вдвое меньшей величины 8/2. Суммарный объем двух внешних областей (части четырех цилиндров за вычетом объемов конусов) равен:
Ve/2
1 1 А 2
„Д 2 - _ 1 „д 2 - 4 = 2 „Д 21 .
. 4 3 4) "
3
К внешним областям деформирования примыкают недеформируемые части тела в форме конусов, в которых деформации отсутствуют 8=0 (рис. 4).
Рис. 4. Расположение областей деформирования с разным удлинением
Пусть недеформируемые части тела под действием нагрузки Р получили малое взаимное перемещение величины 5. Положение границ деформированных областей при этом показано на рис. 5.
Сплошными линиями показаны положения образующих (следов границ областей) областей в деформированном состоянии.
Все элементы внутренней области деформирования испытают удлинение в продольном направлении величины
26 (3)
8=8 =
/
а элементы внешних областей деформирования испытают удлинение величины 82 = ^ = у ■> а
также сдвиговую деформацию.
Суммарный угол сдвига во внешних областях (см. рис.5): у = уг + у2 ,
Я
1_ 2
гЯ I
ёи dz
8 1_ 24 Я
_8/_ 8Д
Ск —1 28 Г я —1
< 1 8Я, ~ / , 1 8Я ,
гЯ г/
Тогда у = уг + уг = -— — = г / о К
Предположим, что действующие в направлениях, перпендикулярных растягивающему усилию кольцевые и радиальные нормальные напряжения пренебрежимо малы, а осевые
<зт
нормальные и касательные напряжения равны максимально возможным: а2 = аГ ;тг2 = — .
Рис. 5. Положение деформированных границ деформированных областей и сдвиговая деформация в наружных областях деформирования
Тогда диссипация энергии в деформируемом объеме будет равна:
И = Уе<5Т8 + V, 12<3Т | + V,/2 ^у =
о 28 ( I I ^ = -пЯ 21ст — + -пЯ 21от- + -%Я2/——
1
6
28 2
8 2
I 3
/ 3
2 /
Я %Я
= кЯ ат8
1 + -
\Л
Я _ _/ / 8Д))
Работа внешней нагрузки: А=Рё. Приравнивая работу внешней нагрузки мощности диссипации и дифференцируя обе части равенства по параметру 5, найдем:
Р = жЯ а,
1 + -
3
Я I
\\
I Ш))
Выражение во внутренних скобках равно нулю ( у = 0 ) при
/ = 2у[2Я,
(4)
тогда
Р = Рпред = = Яе? ,
что равно предельной нагрузке для равномерного деформирования (2). При этом все элементы деформируемой области испытывают одноосное растяжение.
Локальное деформирование. Детальное рассмотрение найденного поля перемещений
г
Рис. 6. Расположение системы координат для описания поля перемещений в левой половине района локального деформирования (шейки)
Рассмотрим полученное приближенное решение для локального деформирования более подробно. Начало отсчета цилиндрической системы координат поместим в вершине конической поверхности, ограничивающей левую недеформируемую область (рис. 6). Получим зависимости для левой половины деформируемого района (шейки). Поле (скоростей) перемещений деформируемого района (для - < г < при данном расположении системы координат) зададим в следующем виде.
Уравнение прямой (образующей конической поверхности), отделяющей область = 0 от области £2 = £ / 2:
ЛЯ г = -4— г . /
Уравнение прямой (образующей конической поверхности), отделяющей область = £ / 2 от области е„ = £ :
„Я
г = 4—г . /
Для внутренней области £2 = е.
Перемещения и вдоль оси г и перемещения вдоль оси 2 :
и = —е г/2, (5)
ж = гг.
Данному полю (скоростей) перемещений соответствует следующее поле (скоростей) деформаций для области £г = е:
ди /о
8г = — = -е/2,
дг
и /о 8е = — = —8/2,
г
дw
8 2 = "& =8'
ди дм>
У гг =— + — = гг дг дг
Для внешних областей 82 =8/2 .
Перемещения и вдоль оси г и перемещения м> вдоль оси Z :
8, 4 R
и=—4(г+~Г2 8 1 I
w = — (--г + 2) .
2 4 Л
Данному полю (скоростей) перемещений соответствует следующее поле (скоростей) деформаций:
ди 8 8 г = дг ~~ 4, и 8.. йг.
8е----4(1 + 4--),
г 4 I г
дw 8 8 2 = &" = 2 , ди дw _ 8Й 81
У гг = &+"дТ = _ Т+■
При этом у гг = 0 при Для I = 2л/2й :
I = 242й.
8„ 42 гч 8 8 42 г
8е =--(1 +--) --------.
е 4 2 г 4 4 2 г
ема максимально и равно 8 +8й + 8 =-----= Н-----¡^ = + —, знак «+» относится
Р г е г 42 г 4 2 42г 8
8 42 г
Изменение объема в точках материала равно 8 г + 8е+8 г = ——- ——.
й й _
На поверхностях раздела областей при г = ±4—г = ±4—-¿=— г = ±42г изменении объ-
I 2л/ 2К
8 42 г _ 8 42 г _ 8
4' 2 г ~
к граничной поверхности недеформируемой области.
8 42 г
Условие сохранения объема: 8г +8е+8г = —-— = 0 выполняется только при
г = 0 и интегрально для области в целом: для любой точки с осевой координатой + г имеется ответная точка с осевой координатой — г, для которой сумма линейных деформаций имеет такую же величину и противоположный знак. Поэтому суммарное изменение объема внешних областей равно нулю.
Для оценки погрешности, вносимой изменением объема, рассмотрим следующие рас-
8
суждения. Среднее изменение объема по модулю во внешних областях равно —. Поскольку
16
продольная линейная деформация здесь равна S, то лишний объем может быть получен добавкой лишней деформации величины —. Тогда погрешность можно оценить, взяв отноше-
16
ние этих величин, которое равно 2/16=0,125 или 12,5 %. Данные рассуждения не являются строгими и полученные результаты нуждаются в дополнительных проверках, которые содержатся в последующих разделах.
Перемещения на границе с внутренней областью sz = s :
s 4R s 4R 7 s
u =--(r H--z) =--(r H---r) =--r ,
4V I ' 2V l 4R 2
s 1 I Л s 1 I R
w = — (--r + z) = — (--4—z + z) = sz,
2 4 R 2 4 R I
что совпадает с (5).
Перемещения на границе с недеформируемой областью s z = 0:
r 4R s 4R I ч Л
u =--(r H--z) =--(r---r) = 0 ,
4 I 4 I 4 R
s Л I ч s Л I I . n
w = — (--r + z) = — (--r--r) = 0 .
2 4 R 2 4 R 4R
Для правой половины деформируемого района (шейки) аналогичные соотношения
могут быть получены путем расположения системы координат симметрично относительно
середины шейки.
Пусть Za - угол наклона образующих конических поверхностей, разграничивающих области деформирования, к оси симметрии (рис.6). Величина этого угла однозначно связана
с длиной района локального деформирования. Для I = 2yÍ2R согласно рис. 6 следует:
R К
tga = ^2^/2^/4 = 2 , при этом a = 54,736°. Относительное удлинение вдоль образующих
разграничивающих конических поверхностей вычислим по известной формуле [13]:
2 • 2 sa=sz cos a + sr sin a + yzr smacosa.
- 2 tg2a 1 . 2 tg2a 2
С учетом cos a =-= — ;sin a =-= — и yxz = 0 находим
1 + tg a 3 1 + tg a 3
s 1 s 2 ^ ^
sa =-----+ 0 = 0.
a 23 4 3
Выводы по результатам детального рассмотрения найденного поля скоростей перемещений:
• при I = 2V2R минимум энергии, поглощенной пластическими деформациями, обеспечивается отсутствием сдвигов в плоскостях, проходящих через ось симметрии;
• поле (скоростей) перемещений непрерывно - на границах областей перемещения одинаковы;
• плоским остается только поперечное сечение внутренней области деформирования плоскостью симметрии;
• все элементы тела находятся в состоянии растяжения напряжениями az = aT аналогично случаю равномерного деформирования и при I = 2V2R предельная нагрузка равна P = nR 2aT = SaT, выполняются условия равновесия, как внутренних элементов, так и нагруженных внешней нагрузкой;
• условие постоянства объема выполняется для внутренней области, но не выполняется для точек внешних областей, однако выполняется для внешних областей в целом, интегрально;
• относительное удлинение вдоль образующих разграничивающих конических поверхностей равно нулю.
Условие локализации деформаций
В работе [11] показано, что условием локализации деформаций при растяжении образца из материала с упрочнением является выполнение неравенства
йа
(--а) < 0
йг
(6)
при равенстве усилий равномерного и локального деформирования. Если это условие выполняется, то при дальнейшем растяжении образца нагрузка уменьшается, причем скорость уменьшения нагрузки для локального деформирования больше, чем для равномерного.
В рассматриваемой задаче одинаковой предельной нагрузке Pпред=Saт соответствуют два различных поля (скоростей) перемещений: равномерное и локальное, для которого поле (скоростей) перемещений определяется указанными зависимостями. Найденное локальное поле скоростей не принадлежит точному решению, однако, поскольку локализация деформаций есть экспериментальный факт, следовательно, имеется локальное поле скоростей, для которого предельная нагрузка равна нагрузке, соответствующей равномерному полю скоростей, а погрешность полученного решения можно отнести к неточности задаваемого поля скоростей. Для проверки способности разработанной аналитической модели отражать существенные особенности процесса локального деформирования и оценки количественной погрешности расчетов по полученным формулам, сопоставим ее с численным моделированием и экспериментальными фактами.
Качественное сравнение.
Форма шейки на испытанных образцах разного диаметра
На рис. 7 показаны фотографии шейки круглых образцов из пластичной стали, диаметры которых отличаются в 2,5 раза.
Рис. 7. Фото образцов разного диаметра с развитой шейкой:
левое фото Do=6мм, правое фото Do=15мм
Сходство формы позволяет заключить, что шейки геометрически подобны и длина района локальных деформаций пропорциональна радиусу образца. Полученное соотношение
/ = 2^ (то есть длина шейки пропорциональна радиусу образца) качественно совпадает с данным наблюдением, следовательно, погрешность в данном соотношении может содержаться только в числовом коэффициенте перед радиусом цилиндра.
Проведенные сопоставления свидетельствуют о качественной правильности полученного аналитического решения.
Результаты конечно элементного моделирования.
Форма недеформируемой части образца при локализации деформаций
На рис.8 приведены результаты конечно элементного моделирования. Показано распределение интенсивности деформаций в верхней половине продольного сечения цилиндрического тела, растягиваемого в горизонтальном направлении. В численном расчете учтены геометрическая и физическая нелинейности. Приведенная на рисунке форма соответствует общему относительному удлинению 4,2375 %. Соотношение размеров цилиндра L/Do=5, длина цилиндра L=0,4 м. Материал - пластический с билинейной зависимостью напряжений от деформаций. Модуль Юнга равен 200000МПа, предел текучести - 200Мпа, модуль упрочнения - 20МПа. Сетка имеет одинаковый шаг по вертикали и горизонтали и содержит по радиусу тела 11 узлов. Нагружение производилось сообщением перемещения левой границы относительно правой.
Отметим, что в численном расчете принят низкий модуль упрочнения с тем, чтобы равномерная деформация была минимальной (см. критерий (6)), и расчету подвергалась в основном деформация в районе шейки (в численном расчете шейка образуется практически без равномерной деформации). Это также приближает постановку численного решения к постановке аналитического, где материал полагается жестко-идеально-пластическим.
Сначала сопоставим результаты данного численного расчета (осесимметричный случай) с результатами аналогичного расчета для случая плоской деформации (рис. 9). Размеры продольного сечения, свойства материала, разбивка на элементы и нагрузка идентичны таковым для случая осевой симметрии.
Рис. 8. Осесимметричная деформация wl = 0,01695м, 80 =0.042375
кк
Рис. 9. Плоская деформация Wl = 0,01695м, 8 0 =0.042375
По результатам численных расчетов определение длины шейки в момент начала ее образования затруднено: начало локализации захватывает область с нечеткой границей, которая затем отчетливо локализуется уже при значительных деформациях. При дальнейшем деформировании хорошо заметны сходство и отличие в форме недеформируемых внешних по отношению к шейке областей. Сходство в том, что в обоих случаях (как осесимметрич-ной, так и плоской деформации) разграничивающая поверхность не является плоским нормальным поперечным сечением, и длина района деформаций возрастает от оси или плоскости симметрии к наружной поверхности. Также в обоих случаях имеет место увеличение деформаций от поверхности внутрь образца в среднем сечении. Отличие заключается в том, что для осесимметричного случая линия, ограничивающая недеформируемый объем (в продольном сечении), заметно меньше отклонена от нормали к оси (плоскости) симметрии.
Для плоской деформации след границы недеформируемого объема может быть адекватно аппроксимирован отрезком прямой под углом около 45 ° к оси симметрии. Для осе-симметричного случая образующая граничной поверхности недеформируемой области искривлена. Однако очевидно, что при аппроксимации отрезком прямой угол наклона к оси симметрии будет больше, чем для случая плоской деформации. Отметим, что данные численные расчеты легко воспроизводимы при определенном опыте и навыках численных расчетов, содержащих физическую и геометрическую нелинейности.
Далее сопоставим численные и аналитические решения. В аналитических моделях след границы недеформируемой области есть отрезок прямой, проведенный под углом к оси симметрии сечения. Полученные в [11] для плоской деформации и в первом разделе настоящей работы для осесимметричной деформации решения показывают различные значения данных углов: 45° и 54,7° соответственно. В обоих случаях величины углов наклона к оси симметрии найдены из условия отсутствия сдвигов в плоскостях продольных сечений. Форма следа разграничивающей поверхности в численных расчетах имеет аналогичную тенденцию (рис. 8, 9).
Отмеченное выше сходство формы граничной поверхности недеформируемой области в аналитических моделях и в численных расчетах, полученных различными способами, подтверждает правильность основных принципов построения аналитических моделей.
Проведенные сопоставления показывают адекватность отражения в рассмотренных расчетных моделях основных закономерностей явления локализации деформации при растяжении цилиндрических тел и позволяют перейти к количественному сравнению аналитических моделей с результатами численных расчетов.
Количественное сравнение расчетов по полученным зависимостям.
Наименьший диаметр шейки
Применим кинематический механизм, использованный для определения предельной нагрузки при локальном деформировании, для расчета конечных деформаций в районе шейки. Получим зависимости для расчета наименьшего диаметра шейки и максимального продольного относительного удлинения. Пусть I - первоначальная длина шейки. Обозначим Wl относительное перемещение недеформируемых областей после образования шейки (рис. 10).
Рис. 10. Схема размеров шейки для расчета конечных деформаций
Используем геометрические соотношения (см. первый раздел данной статьи), заменив 8 (бесконечно малое перемещение) на И7. Тогда, полагая в (3) д=Wl, а в (5) г=Я, найдем уменьшение диаметра в середине шейки:
ДО = 2■ кД)| = 2= 2• = ^ = ^ .
1 1 2 2 I I
Обозначим: От - минимальный диаметр шейки, Оо=2Я - первоначальный диаметр образца. Тогда
А» = А-АО = А
^ = АО - - )■
I 0 I
При этом в соответствии с (4):
Ал = ДО -
w,
л/2Д
> = < О.)•
Относительное поперечное сужение образца:
2 2 лД —кИ т у = —--;-= 1-
г \2
д.
лД
о
V А
= 1-
Щ
Дл/2
Максимальное продольное относительное удлинение в шейке (см. (2) и (3)):
2-и^
8тах Ъп42'
(?)
(8)
(9)
Сопоставим количественные результаты численных расчетов, описанных в третьей части настоящей работы для случая осесимметричной деформации и аналитических расчетов по приведенным в четвертой части формулам для трех значений И7. Распределение интен-
сивностей деформаций и форма продольного сечения для тех же значений Wl приведены на рис. 8, 11, 12. Результаты расчетов существенных параметров сведены в табл. 1. Общее относительное удлинение вычислено как отношение относительного перемещения концов тела
w!
к первоначальной длине: 80 = — .
L
Таблица 1
Результаты расчетов параметров деформирования шейки
wl ,м 80 у по численному расчету у по ф. (8) Расхождение по величине у, % Расхождение по величине ^, % 8 тах числ./ аналитич. Расхождение по величине 8 , %
0,01695 0.042375 0,2141 0,2772 29,45 4,27 0,2818/0,2996 6,33
0,03295 0.082375 0,4735 0,4977 5,1 2,38 0,7336/0,5825 25,94
0,04835 0,122675 0,7244 0,6781 6,8 8,1 1,329/0,8653 53,5
ЭТЕР=1 5ИВ =50 Т1КЕ=.411873 ЕРЮ1НТ (АТО) ОМХ =.,03295 ЗЫТ =. 429Е-03 5МХ =1.09934
Рис. 11. Осесимметричная деформация Wl = 0,03295м, 8 о =0.082375
5ТЕР=1 5ИВ =70 Т1ИЕ=.511873 ЕРТ01НТ (АУБ) ОМХ =.04895 5МН =.2Ё2Е-03 5МХ =1.99456
Рис. 12. Осесимметричная деформация wl = 0,04835м, 80 =0,122675
Значительное расхождение в величине у для общего удлинения 0.042375 связано с тем, что при малом отличии Dm от Do малое различие в величине Dm приводит к значитель-
ному различию в величине у. В остальных сопоставлениях различия в величинах Dm и у в численном и аналитическом решении достаточно малы.
Значительное расхождение в величине 8max для общего удлинения 0.082375 связано с тем, что при больших величинах wi, как следует из результатов численного расчета, деформации в шейке распределяются неравномерно.
За исключением двух указанных случаев, имеет место удовлетворительное совпадение результатов численных и аналитических расчетов. Это показывает, что аналитическая модель может быть использована для расчетной оценки геометрических параметров шейки круглого образца при разрыве. Если известно поперечное сужение у и остаточное удлинение образца вместе с шейкой So, то с использованием соотношения (8) может быть вычислено относительное перемещение недеформируемых областей wi, а, следовательно, равномерное удлинение, что может быть использовано при расчете удельной энергоемкости металлических материалов. Для ориентировочной оценки величины продольной деформации внутри шейки может быть использована зависимость (9).
Заключение
Получено приближенное аналитическое решение задачи о предельном равновесии растягиваемого круглого цилиндрического тела для локального поля скоростей, основанное на предположении об уменьшении длины деформированного района от поверхности к оси симметрии. Решение определяет длину шейки в начале ее образования в долях радиуса цилиндра и коническую форму поверхностей, разграничивающих шейку и недеформируемые внешние области.
Аналитическое решение сопоставлено с результатами численного решения. Обнаружено удовлетворительное совпадение формы поверхностей недеформируемых областей и величины наименьшего диаметра шейки.
Сопоставление формы шейки реальных образцов различного диаметра подтверждает вывод аналитического решения о том, что длина шейки пропорциональна диаметру образца.
Вычисленные по полученным зависимостям размеры шейки могут быть использованы для расчета величины равномерного удлинения при определении удельной энергоемкости металлических материалов по формуле (1).
Информация об удельной энергоемкости материала может быть использована при проектировании элементов конструкции систем пассивной безопасности автомобиля, а также различных видов конструктивной защиты.
Библиографический список
1. Ефремов, А. К. Системы защиты от ударных воздействий // Наука и образование. - 2015. - № 11. - С. 344-369.
2. Пат. 2432295 Российская Федерация, МПК В 63 В 3/24, В 63 В 3/14. Корпус судна с защитой отсека при столкновении с другим судном / Петухов В.Н. - № 2010108777/11; заявл. 09.03.10; опубл. 27.10.11, Бюл. № 30.
3. Правила классификации и постройки атомных судов и плавучих сооружений. НД № 2-020101069. - СПб., 2012. - 151 с.
4. Гончаров, Р.Б. Проблемы поиска оптимальных конструктивных параметров бампера автомобиля при ударном воздействии с позиций пассивной безопасности / Р.Б. Гончаров, В.Н. Зузов // Труды НГТУ им. Р.Е. Алексеева. - 2018. - № 3. - С. 130-136.
5. Фентон, Дж. Несущий каркас кузова автомобиля и его расчет / Дж. Фентон; пер. с англ.; под ред. Э.И. Григолюка. - М.: Машиностроение, 1984. - 200 с.
6. Пат. 2186916 Российская Федерация, МПК Е 04 Н 9/04, Е 04 Н 9/14. Устройство для защиты объектов с атомной энергетической установкой от внешних ударных воздействий / В.Е. Спиро, А.В. Агафонов, А.И. Дульнев, В.И. Лебедев, В.А Никитин., О.М Палий., В.М Пашин., Э.Л Петров, К Б. Проскуряков. - № 2000127925/03; заявл. 08.11.00; опубл. 10.08.02.
7. Малков, В.П. Энергоемкость механических систем / В.П. Малков. - Н. Новгород: Изд-во Нижегородского университета, 1995. - 258 с.
8. Данилов, В.И. О локализации пластической деформации на стадии предразрушения и возможности прогнозирования места и времени вязкого разрыва / В.И. Данилов, Д.В. Орлова, Л.Б. Зуев, И.О. Болотина // Журнал технической физики. - 2011. - Т. 81, вып. 2. - С. 51-57.
9. Данилов, В.И. Типы локализации пластической деформации и стадии диаграмм нагружения металлических материалов с различной кристаллической структурой / В.И. Данилов, Л.Б. Зуев, Е.В. Летахова, Д.В. Орлова, И.А. Охрименко // Прикладная механика и техническая физика. - 2006 -№ 2. - С. 176-184.
10. Давиденков, Н.Н. Анализ напряженного состояния в шейке растянутого образца / Н.Н. Дави-денков, Н.И. Спиридонова // Заводская лаборатория. - 1945. - Т. 11, № 6. - С. 583-593.
11. Уваров, А.И. Аналитическое определение длины района локальной пластической деформации (шейки) при растяжении широкой металлической пластины / А.И. Уваров // Наука и образование: научное издание МГТУ им. Н.Э. Баумана. - 2016. - № 9. - С. 138-154.
12. Малинин, Н.Н. Прикладная теория пластичности и ползучести / Н.Н. Малинин. - М.: Машиностроение, 1975. - 400 с.
13. Биргер, И.А. Сопротивление материалов: Учебное пособие / И.А. Биргер, Р.Р. Мавлютов. - М.: Наука, 1986. - 560 с.
14. Уваров, А.И. Пластические деформации оболочечных элементов конструктивной защиты / А.И. Уваров // Судостроительная промышленность. Серия: Проектирование судов. - 1986. -Вып. 3. - С. 43-48.
Дата поступления в редакцию: 06.02.2019
A.I. Uvarov
DETERMINATION OF THE NECK LENGTH IN A TENSILE ROUND METAL SPECIMEN FOR THE CALCULATED ASSESSMENT OF THE SPECIFIC ENERGY
ABSORPTION OF THE MATERIAL
Nizhny Novgorod state technical university n.a. R.E. Alekseev
Purpose: The aim of this work was to determine the size of the area of the local plastic deformation (neck) tensile cylindrical sample to determine the energy intensity of the material in the calculation of the passive safety of the car. Design/methodology/approach: The process of the beginning of localization of strain (necking) under tensile cylindrical body of plastic material under conditions of axisymmetric strain state.
Findings: Used the kinematic method of limit equilibrium theory. Adopted by the local velocity field of moving consists of two regions of deformation, internal and external, and elongation in the elements of the inner region is twice that in the exterior region. Region of deformation separated from one another and from non-deforming regions inclined planes. From the condition of absence of shear deformations is determined by the length of the area of local deformation, equal to the product of the cylinder diameter to the square root of two. The value of the ultimate tensile load at the local strain equal to the ultimate tensile load with uniform deformation.
Detailed analysis found the velocity field showed that in this case the conditions of equilibrium, the displacement field has no breaks, and separates the surface no deformation. The condition of constancy of the volume as a whole integrally, is observed in the presence of compensating each other of deviations with different signs. Produced qualitative (in the form of non-deforming regions) and quantitative (minimum diameter of the neck) comparison of analytical solutions with numerical simulations. Satisfactory coincidence.
Originality/value: The obtained results can be used to assess the effectiveness of the metal alloys used for the manufacture of structural elements, providing passive safety of the car.
Key words: passive safety. neck in a round sample, specific power capacity.