Напряженно-деформированное состояние и плотность энергии деформации в вершине предельно узких U-вырезов Текст научной статьи по специальности «Физика»

Ovsyannikov Victor Evgenievich, candidate of technical science, docent, vik9800@yandex. ru, Russia, Kurgan, Kurgan State University,

Tereshenko Vladimir Yurievich, postgraduate, tereshenko@,list. ru, Russia, Kurgan, Kurgan State University

УДК 539.375

НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ И ПЛОТНОСТЬ ЭНЕРГИИ ДЕФОРМАЦИИ В ВЕРШИНЕ ПРЕДЕЛЬНО УЗКИХ U-ВЫРЕЗОВ

Ю. Н. Овчаренко

С учетом сингулярных решений в линейной механике разрушения исследовано упругое напряженно-деформированное состояние вершины предельно узкого U - выреза., Получена соответствующая аналитическая формула для определения характеристики плотности энергии деформации в качестве оценки этого состояния.

Ключевые слова: линейная механика разрушения, U-вырез, mode I, mode II, сингулярные решения, плотность энергии деформации.

Работ, посвященных U-вырезам с радиусом в вершине р Ф 0, в настоящее время имеется достаточно много, как теоретических, так и экспериментальных. Однако в представляемой работе интересы лежали в области линейной механики разрушения, что существенно сузило круг изучаемых литературных источников.

В 1966 году в магистерской диссертации 23-летний американский стажер частного Lehigh University аналитически получил [1] формулы для описания напряженного состояния у вершины предельно узкого U-выреза («blunt-трещины») для задач mode I и mode II линейной механики разрушения. Покажем формулы:

- " 1 f

s x

1 4.

s y - KI r-2 V2P 1 f 41 1 f

_T xy _ _ 4 .

В5 3cos — + cos—В 2 2

< В 5 В л 5cos — cos — В

2 2 у

• В . 5 В л - sin — + sin — В 22

P 3 В

- — cos—В

2r 2

p3 + — cos—В

2r 2

-—sin—В

2r 2

О

У

t

xy

К

II

V2P

4 1( 4

1 (_ . 0 . 5ПЛ

7 sin — + sin—0 22

V 1 ( 4

• 0 • 5 0Л

sin — + sin—0 22

J

05 3cos — + cos — 0

v 2 2 y

P • 3 0

- —sin—0

2r 2

P • 3 0

—sin—0 2r 2

P 3 0

■—cos—0 2r 2

(2)

Привлекает в (1) и (2) простота и непосредственная связь с известными формулами для трещин Вестергарда-Ирвина (р = 0). Формулы предлагается применять, когда a priori известны коэффициенты интенсивности напряжений Ki и Кц для того же тела, но с виртуальной трещиной, аналогичной узкому U-вырезу по плоскостному расположению и размеру.

Рис. 1. Рассматриваемые в [1]узкие: а - эллиптическое отверстие; б - гиперболический вырез

Схемы ^вырезов, с необходимой информацией, показаны на рис.

1. Реально такой вид может иметь обычная трещина, подвергшаяся коррозионному воздействию среды, или трещиноподобный дефект в сварном соединении (непровар, подрез), или специально выполненная узкая прорезь в какой-либо детали. В теоретическом плане интерес к узкому ^вырезу заключается в том, что реальная трещина обычно затупляется до какого-то р прежде, чем стартует, и это по возможности следует учитывать. Есть также теоретический подход при хрупком разрушении, когда предполагается существование весьма малой цилиндрической центральной зоны радиуса р у вершины трещины, где материал нельзя считать сплошной средой [2].

Обращает на себя внимание, что начало полярной системы координат находится на расстоянии Р от дна ^выреза. Именно такое расположение полярной системы позволило стажеру Greager получить вышеуказанную форму для представления напряжений.

Для получения формул (1) и (2) Creager применил комплексные

функции j(z) и c(z) для эллиптического отверстия и комплексные функции j'(z) и c'(z) для гиперболического выреза из монографии [3]. Для эллиптического отверстия это:

4j(z) = Sc e 2Xo cos 2 g cosh Z + (l - e 2X0 +2ig )sinh Z 4C(z ) = -Sc 2

(cosh 2X0 - cos 2g)Z +1 e2X0 cosh 2(Z - Xo - ig)

(З)

Здесь С = Х + *Л, где X и Л - эллиптические координаты; с и X о -параметры эллиптического отверстия; £ - внешняя нагрузка с углом приложения у.

Для гиперболического выреза:

Ф (* ) =

2(p - 2h o + sin 2h o) a sinh Z

,/ч iF|z + (l -2cos2 ho )cothZ| С (z ) = —^ „--------■ J ч—

iz)

2(p- 2ho + sin 2%)

где с и ho - параметры гиперболического выреза; F - внешняя нагрузка вдоль вертикальной оси у.

Creager рассматривал предельно узкие как эллиптическое отверстие, так и гиперболический вырез, где р очень мало. Он применял при рассмотрении вышеуказанных комплексных функций различные предельные переходы с отбрасыванием членов высокой степени малости, а также некоторые упрощающие условия. В результате получил тождественный результат, как для эллиптического отверстия, так и для гиперболического выреза. А именно.

Для mode I:

1

^ о Kj -2 е оx + оy = 2 г— r 2 cos

V2P

о

Kj

2

1,

y

оx + 2itxy '-r 2

p

i sin е + —

Л -Де

(4)

Для mode II:

Kjj

о x + о v = 2r

x y V2P

l

2 ■ е

2 sin—

о y о x + 2it xy

= KJJ V2P

r

2

-1 г 2

,p

i2cos е-i — - sin е

V r у

\ -Де 2

(5)

е

Исходя из формул (4) и (5), были получены вышеприведенные формулы (1) и (2). Они использованы в дальнейших работах Creager и дру-

2

е

r

so

гих авторов [4, 5] для прочностных исследований «blunt-трещин» в коррозионных средах и при некоторых видах усталостных нагружений.

В представляемой работе на основе сингулярных решений линейной механики разрушения проведено дальнейшее исследование упругого напряженно-деформированного состояния (НДС) около вершины предельно узкого U-выреза. Получена соответствующая аналитическая формула для определения характеристики плотности энергии деформации в качестве оценки этого состояния.

Теоретические исследования. Первоначальное знакомство с формулами (1) и (2) происходило по косвенным работам [4, 5], так как сама магистерская диссертация Creager была недоступна. Анализ конструкции вышеуказанных формул показал, что они могут быть получены из следующего представления комплексных функций j( z) и y(z) [6].

Для mode I :

где C1

= Kl .

V2P

Для mode II :

j(z ) = C1z 2 ¥(z ) = C1

1 2

— z 2 2

P z

О 2

У

(6)

где C2 =

K

ll

j(z ) = -C2 z 2

y(z ) = iC2

1

3 2 ■

— z 2 + P z 2

о 2

У

(7)

Выражения (6) и (7) при р = 0 - это соответствующее формулам (1) и (2) (тоже при р = 0) комплексное описание трещины в линейной механике разрушения. Покажем.

Общее представление функции напряжений Э. Гурса через две комплексные функции ф(г) и у(г) имеет вид:

и (г) = 2 [ф) + 1ф(1) + %(г) + %(г)], (8)

где %'(2 ) = у(7 )• Используя результаты работы [7], можно получить описание этой комплексной функции конкретно для трещины:

1

1

1

U (z ) =l 2

azz2 + azz2 +

2b

з

з

z

(9)

где С = 1 (а + a), С2 = —(а - a), Cз = 1 (ь + Ь ), C4 = — (ь - Ь ). Здесь а и Ь

- комплексные постоянные коэффициенты. Из простого сравнения (8) и (9) следует, что

1

• 2

1’ (10)

j(z) = az ■

y(z ) = bz 2.

l

Для mode I имеем [б]: а = Cl, b = — Cl. Для mode II:a = -iC2,

з

b = 1—C2. Что и требовалось показать.

Если продифференцировать (10), то имеем известные комплексные функции Вестергарда для описания напряженного состояния тела с трещи-

1

ной. Дополнительные слагаемые + рz 2 в формулах (6) и (7) были легко подобраны, исходя из внешнего вида (1) и (2).

В дальнейшем представилась возможность познакомиться с самой магистерской диссертацией Creager [1], в частности с формулами (4) и (5). Выяснилось, что предложенные автором статьи формулы (6) и (7) непосредственно следуют из (4) и (5). Ошибки не случилось.

Комплексные функции (6) и (7) позволяют без особого труда добавить к (1) и (2) формулы для перемещений, исходя из известного [8] комплексного выражения:

2ц (u + iv) = k j(z) - zjz)- y(z).

Для mode I в прямоугольной системе координат имеем :

у 1Л

2m u 2m v

= Cl r

k —

V 2 У

(/ 11 k +-----

V 2 У

0 l з p 0

cos------cos—0 + — cos—

2 2 2 r 2

0 l з p 0

sin------sin — 0 + — sin —

2 2 2 r 2

(11)

где m - модуль упругости II рода; к

з-v 1+v

для плоского напряженного

состояния; к = 3 - 4у - для плоской деформации; V - коэффициент Пуассона.

l

l

з

Для mode II:

2m u 2m v

= C2r:

Ґ Зл к + З

v 2 у

Ґ Зл - к + З 2

е l З p е

sin— + — sin—е — sin—

2 2 2 r 2

е l З p е

cos--cos—е + cos—

2 2 2 r 2

(12)

Тем самым завершен полный набор асимптотических формул для описания НДС в вершине предельно узкого и-выреза.

Плотность энергии деформации. Ясно, что применять напрямую коэффициенты интенсивности напряжений С\ и С2 для оценки НДС у вершины узкого И-выреза не корректно. Наличие радиуса р у вершины, хотя и очень малого, вносит существенное изменение. По сравнению с острой трещиной сингулярность напряжений пропадает, они становятся конечными и с другими законами развития. Для оценки соответствующего НДС можно использовать любые подходы, которые имеются в теории упругости. Однако повторимся, интересы автора лежали в области механики разрушения. Поэтому в работе предлагается оценивать НДС у вершины рассматриваемого и-выреза через плотность энергии упругой деформации. Эта характеристика вписана не только в арсенал теории упругости, но и линейной механики разрушения [9, 10].

В полярной системе координат плотность энергии деформации можно представить как

w = 1

2

Э ur Э r

+ Ое

ur + 1 Э u е r r Эе

+ Т гЄ

1 Э u

r + Э u е

r Эе Э r

ие

r

. (1З)

Выпишем формулы для напряжений и перемещений в полярной системе координат, используя основополагающие комплексные выражения [8]: _

О к +°е= 2[ф'(ї )+ф'(ї)] ,

00-0 г + 2/т ге = 2е 2гЄ [ї ф* (г) + у * (г)] ,

2т(иг + /ие) = е _г0 [ф (ї) - їф* (ї) - у(ї)].

Суммарно для mode I и mode II имеем:

« 0 З е Л

5cos—cos—е

2 2 у

^ е З е Л

3cos— + cos—е

2 2 у

_ _ ' 1 f

О r

1 2 v

=1 Clr 2 2l 1 r 2 v 1 f

Ое

_х г0 _ _ 2 v

p е

— — cos— r2

p е

+ — cos— r2

е • З 0л

sin— + sin—е

2 2 у

p . е

+ — sin— r2

+

1

8З

1

1

+ _ Со г 2 2 2

1 ( * • 0 о • 3 ^

5 БІП — + 3 Бій — 0

2 2 у

2

V

3

р . 0

— Бій —

г 2

2

03 бій — + бій—0 22

1 а і 2 і 1 - С1Г2

2ц и 0

1 ( 0 3

— СОБ— + 3С0Б — 0

2 V 2 2

(* -1л

р0

+ — бій — г2

Р 0

- — СОБ — г2

V

2

0 1 3. р 1.

СОБ--------СОБ — 0 + — СОБ — 0

2 2 2 г 2

ґ 1 ^ . 0 1 . 3Л р . 0

бій —+ — бій — 0 - — бій —

2 2 2 г 2

к +

V 2

у

+

+ ^2 г-

( 1 ^

к —

V 2 у

( 1 ^

— к +

V 2 у

0 3 . 3 л р . 0

бій — + — бій — 0 + — бій —

2 2 2 г 2

0 3 3 р 0

СОБ— + — СОБ — 0 + — СОБ — 2 2 2 г 2

Эти формулы для напряжений и перемещений подставляем в выражение для плотности энергии деформации (11). После операций дифференцирования и большого числа элементарных преобразований имеем:

Ж - г (а11С12 + 2а12С1С2 + а22С2

(14)

где

"а11" " ( ^ 2 (1 + СОБ 0)(к _ СОБ 0)+ Р Vг у

а12 = 1 8ц [2 соб 0 — (к — 1)] бій 0 — — бій 0 г

_а22 _ [к(1 _ СОБ 0) + (1 + 3 СОБ 0) СОБ 0] _ — СОБ 0 + _г 2 (Р 1 V г у

При р = 0 выражение (14) переходит в известное выражение [9, 10] для плотности энергии деформации в механике острых трещин.

Логично считать, что разрушение в точке с координатами 0с, гс на, контуре вершины и-выреза происходит тогда, когда плотность энергии деформации Ж в этом месте равна или выше ее критической величины Жс (гипотеза Бельтрами):

Ж > Жс. (15)

Следует обратить внимание на следующее. В эллиптических координатах плотность энергии деформации Ж на поверхности выреза пред-

ставляется как Ж -1 , где =°Г + °0 и еЛ =

1 -V'

Е

(плоская де-

1

формация). Фактически имеет место критерий максимального напряжения . Это отметил БШ О.С. [11].

Некоторые расчеты и графические построения [12]. Из выражения для о у в формулах (1) можно найти связь между коэффициентом интенсивности напряжений С1 и упругим коэффициентом концентрации напряжений ао в вершине узкого выреза с некоторым малым радиусом р.

Пусть 0 = 0, г = Р, о у = аО оном, где оном - номинальные напряжения. 2

Имеем

С1 = 2^2°номlimаал/Р, (16)

р ® 0.

Отсюда может быть предложен способ нахождения С1 через коэффициенты концентрации напряжений аО. Он заключается в том, что подсчитываются несколько значений правой части формулы (16) при изменении р и, соответственно, аО. Эти значения ложатся, как правило, на некоторую прямую на графике зависимости «правой части формулы» от параметра р . Эта прямая затем экстраполируется на значение радиуса р = 0 . В результате графическим путем находится искомое значение С.

Рис. 2. Графическое определение коэффициента С

На рис. 2 показан пример графического определения С1 для плоского симметричного образца с узким ^образным вырезом при чистом изгибе. Коэффициенты концентрации аО взяты из монографии Петерсона

[13]. Точность полученных таким способом результатов для оценки может быть довольно высокой. Однако данный метод можно рассматри-

вать скорее как принципиальный в плане связи С и аа, нежели практический.

Рис. 3. Эпюры напряжений оx, оу и в сечении 0 = 0

На рис. 3 показаны характерные графики для напряжений о x, о у и

txy. Они построены для 0 = 0, и условных Q = С2 = 1, р = 1. На рис. 3, а

показана эпюра напряжений о у - mode I. Именно эти напряжения будут

играть основополагающую роль в разрушении. Несмотря на их резкое возрастание при приближении к кромке U-трещины, они конечны. Пунктиром также показано как выглядела бы эпюра указанных напряжений, если бы предельно узкий U-вырез выродился в острую трещину. Пути резкого возрастания напряжений заметно различаются. Из рис. 3, б следует, что напряжения оx - mode I и txy - mode II имеют одинаковую эпюру (ввиду

тождественности аналитического описания). Она уходит в ноль у кромки U-выреза, что и следовало ожидать. Пунктиром показано, как бы выглядела эпюра указанных напряжений, если бы узкий U-вырез выродился в острую трещину с вершиной в начале координат О . Разница обнаруживается существенная в связи с отсутствием или наличием сингулярности.

На рис. 4 показаны эпюры плотности энергии деформации на кромках вершины предельно узкого U-выреза для mode I и mode II. Они построены по формулам (14) для эллиптического контура.

а б

Рис.4. Эпюры плотности энергии деформации: а - mode I; б - mode II

Эпюры построены для задачи плоской деформации при следующих

исходных данных. Модуль сдвига m =——E—т, где модуль упругости

2 (1 + v)

E = 1, коэффициент Пуассона v = 0,25. Коэффициенты интенсивности напряжений Q = С2 = 1, радиус в вершине U-выреза р = 1. Для нормального разрыва (mode I ) максимальная плотность W наблюдается (рис. 4, а) для точки контура в направлении 0 = 0 . Это направление возможного разрушения не вызывает сомнений ни с точки зрения теории, ни с точки зрения практики.

Для поперечного сдвига (mode II) расчеты показывают, что максимальная плотность W наблюдается при 0 = ±9О,2 °. Положительные значения 0 соответствуют углу, где окружные напряжения о0 отрицательны, а отрицательные значения 0 - углу, где о0 положительны. Из здравого смысла следует, что разрушение можно ожидать только в точке, где окружные напряжения о0 положительны, т.е. растягивающие, что обеспечивает раскрытие возможной микротрещины. В данном примере это имеет

место при 0 = -90,2 °.

Поговорим о точности полученных асимптотических формул (14) для оценки плотности энергии деформации W кромке вершины предельно узкого U-выреза. В качестве объекта исследования выбрано mode I и направление 0 = 0 . Сравнивались точное решение, полученное по формуле

W = 2(ох ex +оу еу) (17)

и асимптотическое (14), представленное в виде:

W = — S 2 -r 8m 2

2(k — 1) +

2

Р

V r J

(18)

Для (17) точные значения напряжений подсчитывались по комплексным выражениям (3), которые после необходимых преобразований имели вид:

о х = - (А - В - С),

о, = | (А + В + С),

где А = 1 + е2X0) ^п1'------е2Хо ;

' соэЬ 2Х -1

В = ^ со$Ь^ -е2^0 + со$Ь2Х0 -е-Хо (2э1пЬ2X- 1)$тЬ2с;0 2м1Л3 X

C = e2^0 cosh 2X0

Точные значения деформаций ех и еу определялись из обобщенного закона Гука для плоской деформации.

На рис. 5 показаны две графической зависимости относительных погрешностей вычисления плотности энергии деформации Ж по асимптотической формуле (18) в сравнении с точными решениями по формуле (17). Одна графическая зависимость соответствует р = 1, другая р = 0,1. В обоих вариантах принималось: £ = 1, т = 1, а = 100, у = 0,3. Из рисунка следует, что для обоих вариантов погрешность вблизи вершины и-выреза (х Р102) находится в пределах ± 10%, что для практических расчетов вполне допустимо. Вариант с р = 0,1 дает четкую картину уменьшения погрешности при приближении к контуру вершины выреза вплоть до «ми-нус»3%. Т. е. с уменьшением радиуса р решения у контура вершины для

Ж по формуле (14) становятся все более точными. Как следствие, это характеризует также точность представленных в работе асимптотических формул для напряжений и перемещений.

‘Ошн. погр., %

п - П 1

о~1 О

1021——-

1оо — 104 106 108

Рис. 5. Относительная погрешность расчета W по асимптотической формуле (14)

Выводы

1. Представлены более приемлемые для исследований НДС у вершины предельно узкого U-выреза выражения комплексных функций j(z) и y(z). Получено выражение для плотности энергии деформации через коэффициенты интенсивности напряжений для mode I и mode II, которое предлагается использовать в качестве оценочного показателя НДС у вершины.

2. Предлагается считать, что трещинообразование в вершине происходит, когда плотность энергии деформации на ее кромке при некоторой координате 0 достигает критической величины Wc.

3. С уменьшением радиуса р асимптотические решения для плотности энергии деформации W на контуре вершины U-выреза по формуле

(14) становятся все более точными. Это характеризует также точность представленных в работе асимптотических формул для напряжений и перемещений.

Работа представлена на второй Международной Интернет-конференции по металлургии и металлообработке, проведенной в ТулГУ

1 мая - 30 июня 2013 г.

Список литературы

1. Creager Matthew. The Elastic Stress Field Near the Tip of a Blunt Crack: Master’s Thesis. Lehigh University, 1966. 40 p.

2. Sih G.C. Strain energy density factor applied to mixed mode crack problems // International Journal of Fracture. 1974. V. 10, No.3. P. 305-321.

3. Тимошенко С.П., Гудьер Дж. Теория упругости: пер. англ. / под ред. Г.С. Шапиро. 2-е изд. М.: Наука, 1979. 560 с.

4. Creager M., Paris P. C. Elastic Field Equations for Blunt Cracks with Reference to Stress Corrosion Cracking // Int. J. Fract. Mech. 1967. Vol. 3, № 4. P. 247-252.

5. Heckel K., Wagner R. The tensile fatigue behavior of CT-specimens with small notch root radius. Int. J. Fract., 1975, Vol. 11, № 1. P. 135-140.

6. Овчаренко Ю.Н. Теория и практика V-образных вырезов в механике разрушения: Монография / под ред. Г.И. Макарова. Тула: ТулГУ, 2003. 168 с.

7. Каландия А.И. Замечания об особенности упругих решений вблизи углов // Прикладная математика и механика. 1969. Т. 33. Вып. 1. С. 132-135.

8. Мусхелишвили Н.И. некоторые основные задачи математической теории упругости. М.: Наука, 1966. 707 с.

9. Партон В.З., Морозов Е.М. Механика упругопластического разрушения. - 2-е изд., перераб. И доп. М.: Наука. Главная редакция физикоматематической литературы, 1985. 504 с.

10. Плювинаж Г. Механика упругопластического разрушения: Пер. с франц. М.: Мир, 1993. 450 с.

11. Sih G.C. Strain energy density and surface layer energy for blunt cracks or notches. Mechanics of fracture. Edited by G.E. Sih. The Netherlands, 1978, Vol. 5, p. XIII-XXXVI.

12. Овчаренко Ю.Н. Упругое напряженно-деформированное состояние и плотность энергии деформации у вершины предельно узких U-вырезов // Изв. ТулГУ. Естественные науки. 2010. Вып. 2. С.97-108.

13. Петерсон Р. Коэффициенты концентрации напряжений. М.: Мир, 1977. 302 с.

Овчаренко Юрий Николаевич, канд. техн. наук, доц., ovcharenkosarambler.ru, Россия, Тула, Тульский государственный университет

THE ELASTIC STRESS-DEFORMED CONDITIONS AND DENSITY OF ENERGY OF DEFORMA TION A T TOP EXTREMELY

NARROW U-NOTCHES

J.N. Ovcharenko

In view of singular decisions in the linear mechanics of destruction the elastic stress-deformed condition about top extremely narrow U-выреза. is investigated. The corresponding analytical formula for definition of the characteristic of density of energy of deformation as an estimation of this condition is received.

Key words: linear mechanics of destruction, U-notch, mode I, mode II, singular decisions, density of deformation energy.

Ovcharenko Jury Nikolaevich, candidate of technical sciences, docent, ovcharenkosarambler.ru, Russia, Tula, Tula State University

УДК 681.513.2

АВТОМАТИЗАЦИЯ ПОДГОТОВКИ УПРАВЛЯЮЩИХ ПРОГРАММ ДЛЯ СТАНКОВ С ЧПУ НА ОСНОВЕ МЕТОДОЛОГИИ РАСПОЗНАВАНИЯ ОБРАЗОВ

А.Б. Орлов, И.А. Антамонов

В настоящей статье предложена новая методика автоматизации подготовки управляющих программ для станков с ЧПУ, основанная на базе теории распознавания образов. Рассмотрены преимущества и недостатки существующих способов нанесения гравировки, а так же произведен анализ некоторых алгоритмов теории компьютерного зрения.

Ключевые слова: гравировка, компьютерное зрение, распознавание образов, станки с ЧПУ, Graph Cut, детектор границ Кенни.

В настоящее время существует значительное количество различных информационных систем, позволяющих автоматизировать разработку управляющих программ для станков с ЧПУ. Особенно эффективными эти системы могут быть в случае обработки сложных поверхностей, требующих многопроходной обработки с различными траекториями и возможной сменой инструмента. Одним из примеров подобных сложных поверхностей являются художественные гравировки на различных изделиях. Под гравировкой понимается процесс, связанный с нанесением изображения (текстов, узоров, геометрических фигур, логотипов и так далее) путем неполного погружения гравировального инструмента в какой-либо материал. Такими материалами могут быть различные виды металлов и сплавов, ак-

90