Трещинообразование в круговом диске под действием объемных сил Текст научной статьи по специальности «Физика»

Теория упругости

ТРЕЩИНООБРАЗОВАНИЕ В КРУГОВОМ ДИСКЕ ПОД ДЕЙСТВИЕМ

ОБЪЕМНЫХ СИЛ

Н.М. КАЛАНТАРЛЫ, канд. физ.-мат. наук, доцент Институт математики и механики НАН Азербайджана, AZ1141, ул. Б. Вахабзаде, 9, e-mail: nailyak1975@gmail.com

Рассматривается круговой металлический диск под действием объемных сил. Дается решение задачи механики разрушения о зарождении трещин в круговом диске. Используется модель зоны ослабленных межчастичных связей в состоянии пластического течения при постоянном напряжении.

КЛЮЧЕВЫЕ СЛОВА: круговой диск, зоны ослабленных межчастичных связей материала, трещинооброзование, объемные силы.

Круговые диски широко применяются в современных машинах, в паровых и газовых турбинах, в компрессорах, вентиляторах и машинах химической промышленности. Диски подвергаются нагрузкам, вызывающим их растяжение. Задача о трещинообразовании в конструкциях и изделиях является важной проблемой теории прочности [1-8]. Разработка расчетных моделей исследования трещинообразования в дисках представляет собой актуальную проблему механики материалов.

Рассмотрим металлический круговой диск в процессе работы. Считается, что на частицы среды действуют объемные силы F = X + iY , где X(x, y) и

Y (x, y) заданные функции; i = . В процессе деформации в некоторых точках диска могут появляться зоны, в которых закон Гука не выполняется, т.е. в этих областях напряжения превосходят предел упругости. Указанные зоны (прослойки перенапряженного материала) малы по сравнению с остальной упругой частью диска, и их можно мысленно удалить, заменив разрезами, поверхности которых взаимодействуют между собой по некоторому закону, соответствующему действию удаленного материала. По мере нагружения металлического кругового диска в нем будут возникать зоны предразрушения (области ослабленных межчастичных связей материала), где материал диска деформируется за пределом упругости. В исследуемом случае возникновение зародышевой трещины представляет собой процесс перехода зоны предразрушения в область разорванных связей между поверхностями диска. Размер зоны ослабленных межчастичных связей материала при этом является неизвестным и должен быть определен в процессе решения задачи механики разрушения. Зона предразру-шения ориентирована в направлении максимальных растягивающих напряжений в диске. Исследования [9-11] возникновения областей с нарушенной структурой материала показывают, что в начальной стадии зоны предразрушения представляют собой узкий вытянутый слой, затем с ростом нагрузки внезапно появляется вторичная система зон, содержащих материал с частично нарушенными связями.

Будем считать, что в диске имеется N прямолинейных зон предразрушения длиной 2lk (k = 1,2,...N) (рис. 1). В центрах зон предразрушения разместим начала локальных систем координат х^О^уъ ось xk которой совпадают с линией скольжения и образуют углы ak с осью х (в = 0). Полагаем, что в зоне предразрушения имеет место пластическое течение при постоянном напряжении.

Отнесем диск к полярной системе координат гв , выбрав начало координат в центре окружности L с радиусом R (рис. 1). Граничные условия рассматриваемой задачи механики разрушения о зарождении трещины в круговом диске имеют вид:

«г = 0; тгв = 0 при r=R (1)

на берегах зоны предразрушения:

\хк\< £k, к=1,2,...^ (2)

«Ук = « ;

ТхкУк Т:

где о8 - предел текучести материала диска на растяжение; т. - предел текучести материала диска на сдвиг; аг, ад, тгд - компоненты тензора напряжений в полярных координатах; 1к - полудлины зон предразрушения.

,, Для определения значений внешней

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

1 :__х

{/ 1 0 к/ \ А \

[2]:

(ик~ик)" Ф+- Ч)

= 8С , к = 1, 2, ..., N

Рис. 1. Расчетная схема задачи о

где (ид! - ик. ) - нормальная составляющая раскрытия берегов ¿-той зоны предразрушения; (и+ - и-) - касательная состав-заРождении трещин в круг°в°м даже ляющая раскрытия берегов к-той зоны предразрушения; 8С - характеристика сопротивления материала диска трещи-

нообразованию, определяемая опытным путем [10].

При наличии объемных сил решение задачи представим в виде суммы:

(3)

« = « + «1;

0 , 1 «У = «У + «у:

Т = Г0 +Г1

ху ху ху

где , ст°, тХу - любое частное решение уравнений плоской теории упругости

при наличии объемных сил; «I, «1у , Г - решение уравнений плоской теории

упругости при отсутствии объемных сил. Для напряжений а0, «0 , тХу имеем формулы [12]:

«х + «у =-

-^Яе «,

1 + к дz

(4)

1

а - а + 2го =

у х ^ 1 + к дг

д(кё - ),

в которые входят две функции Q(z, г) и ¿1 (г, г), представляющие собой любые частные решения уравнений

д ^

¿ = ¿(г, г), ^ = ЖЩ.

дгдг ^ ' дг2 V ' Здесь к - постоянная Мусхелишвили для плоского напряженного состояния. Используя формулы Колосова-Мусхелишвили [13]:

«1 + «1 = 4ЯеФ(г), «1 -«1 + 21Т1ху = 2[гФ'(г) + ^(г)],

граничные условия (1) - (2) с учетом (3) запишем в виде краевой задачи для отыскания комплексных потенциалов ф(г ) и %(z ) в виде

Ф(т) + Ф(Т)-e2'в[|ф'(т) + ^(г)] = -(тг0 -iT%), (6)

Ф^) + Ф^ + hФ'(гк) + %(tk) = ^ + 1ts - (^ + ixltyt ), k = 1, 2, ..., N.

Здесь т = R exp(ie) ; tk - аффикс точек берегов зоны предразрушения. Комплексные потенциалы, дающие решение краевой задачи (6), ищем в виде:

Ф^ ) = Ф 0 (z ) + Ф! (z ) ; %(z ) = % (z ) + % (z ), (7)

где потенциалы Ф1 (z ) и % (z ) ищем в виде

! N j k

Mz)=-L § J

gk (t)dt

2m -1 t - zk

jk

(8)

N jk LU k=1 -j.

gk (t) t - z,

Tke

(t - zk ):

-gk (t )

dt ,

Здесь Tk = teliXk + z0; zk = e iak (z - z0); gk (t) - искомые функции, характеризующие раскрытие берегов зон предразрушения:

gk(xk) =

2/и ô i(1 + k) ÔXk

[u+ (Xk ,0) - u+ (Xk ,0) + i(v+ (Xk ,0) - v- (Xk ,0))].

Неизвестные функции gk (t) и комплексные потенциалы Ф 0 (z) и (z) должны быть определены из граничных условий на контуре L (r = R) и берегах зон предразрушения. Используя (7) - (8) для отыскания комплексных потенциалов Ф 0 (z ) и (z ), краевое условие (6) на границе диска представим в виде

Ф 0 (т)+Ф 0 (т) -

-2ie

тФ'

ÖD+= -.0 + iT0e )- (f1 (e)+if2 (e)), (9)

где ^ (0) + if2 0)=Ф1 (г) + Ф1 (г) - е [г Ф^ + ^Дг)]. Для решения краевой задачи (9) использовали метод Н.И. Мусхелишвили [13]:

Ф0 (z) = 2- JN* (т)Г^- - 2- IdT

2m J It- z 2т

% (z )= Л Ф 0 (z ) + Л Ф 0 f1V1Ф0 (z )

z z v z J z

+1

N

■(t)=-(t0 + iT0e )-(f1 (e) + if2 (e)).

После интегрирования и некоторых преобразований получаем:

Ф 0>-2m § I {[ жл+2 )т"е^-(t >+

e fdt + Ф 00 (z ),

Tk z2 Tk - 2z + Tk

( zTk -1)2

1 N j

% >=2m § j

k=1 -j k

. —3

eiak T ___

-"—T-y gk (t) + ( z2 t2 + 4 - 3zTk + ( zTk-1)

+

-2 —Л Tke + zT^Tk -3TkTk |-=

-гак

( zTk -1)

-gk

It)

+ ЛФ00 (z) + Лф00( - I- -ф0с (z);

z ) z

ф*« = ¿J^ + - z 2r

1 -± Vr.

Удовлетворяя функциями (7) - (8), граничным условиям (6) на берегах зон предразрушения, после некоторых преобразований, получим систему N комплексных сингулярных интегральных уравнений относительно неизвестных

функций gk (xk )

N e k р 1

Z {R* (t, Xk )gk (t ) + S„k (t, Xk )gjt)\dt = n[as - irs - (a°yi - ir0^ )+ /0 (Xk )J. (10)

k=1 -e k

Здесь:

/0 (Xk ) = -[ф00 (Xk ) + Ф00 (Xk ) + Xk ф00(Xk ) + %0 (Xk )] ,

^00 (z ) = \ ф 00 (z ) + -1 ф 00 f1 ] -1 ф 00 (z ) ; z z ^ z ) z

Rnk Xk )= e

z

,iak (

V

e -2iak ---+ =——

Tk - Xn T - X n

— 2 X„T k

-= +

1 - XnTk

+

X2nTk - 2Xn + Tk + e-2ian 2Xn (t»T» -1)+Tk(Xn+ Tk\XnTk - з)+ 4Tk (1 -TkXn)2 (1 -TkXn)3

Snk (t, Xk ) =

1

T - X

Âk лп e-2ia„

-iak

Tk - Xn (Tk - Xn )

Tk2+ Xn2Tk - 2Xn + Tk jk2 (Xn - Tk )

n k + k \ n k>

1 - TkXn (1 - XnTk )2

(1 - TkXn)

Xn = X„eia" + z0.

n n n

Переменные Xn, z, zn, e n - безразмерные, отнесенные к радиусу диска R.

Для внутренних зон предразрушения к системе сингулярных интегральных уравнений необходимо добавить дополнительные равенства

k

{ gk (t )V

t = 0,

(11)

-e »

обеспечивающие однозначность перемещений при обходе контуров зон пред-разрушения.

Система комплексных сингулярных интегральных уравнений (10) при дополнительных условиях (11) с помощью процедуры алгебраизации [9, 14] сводится к системе N*M алгебраических уравнений для определения N*M неизвестных gn(tm) (k = 1, 2, ..., N; m = 1, 2, ...,M):

1 M N

MM ZZ e k

M m =1 k=1

gk (tm )Rnk (e ktm,e nXr ) + gk (tm )Snk (e ktm,e nXr )

(12)

= ^ -irs - (< (Xr) -r Vt (Xr)) + /0(Xr) (n = 1, 2, ..., N; r = 1,2,.,M-1),

XtÄ v r '

z

z

-ia

e

2

e

e

2

2

E g°n(tm)=о.

m=1

„ 2m -1 r

Здесь: tm = cos-ж ; xr = cos — я.

m 2M M

Перейдя в (12) к комплексно сопряженным значениям, получим еще N*M алгебраических уравнений. Решение каждого интегрального уравнения ищется в классе всюду ограниченных функций (напряжений), так как напряжения в диске всюду ограничены. Полученные алгебраические системы уравнений не являются пока замкнутыми, так как для их замкнутости не хватает 2N уравнений, выражающих условия разрешимости интегрального уравнения. Для замыкания алгебраических систем необходимо добавить условия конечности напряжений у вершин зон предразрушения (условия разрешимости краевой задачи в классе всюду ограниченных функций). В случае, когда все N зон предразруше-ния являются внутренними, эти условия можно представить в виде: м 2 m — 1

E (— l)m g0 (tm )ctg m-1 ж = о (n = 1, 2, ..., N), (13)

m=1

E (— 1)M '"SO ^m )tg1M Ж = 0.

m =1

Совокупность 2N уравнений позволяет определить размеры зон предраз-рушения. Из-за неизвестных заранее размеров зон предразрушения алгебраическая система (12), (13) оказалась нелинейной. Для ее решения сначала отделяем действительные и мнимые части. В результате получаем две действительные системы, состоящая каждая из N*M+2N уравнений.

Для решения нелинейных алгебраических систем использовали метод последовательных приближений. Используя решение задачи о зонах пластического течения в круговом диске, находим смещения на берегах зон предразруше-ния

1 + к xс

— I gk(xk)dxk = Vk(xk,0) — iuk(xk,0), —■5 k

Здесь ц - модуль сдвига материала диска;

vk(xk,0) = V+(xk,0) — V (xk,0); uk(xk,0)=uk+(xk,0)—uk(xk,0).

Смещения на берегах зоны предразрушения при xk = x0 будут

1 + Kx0C

— — I gk(xk) dxk = Vk(x0k ,0) — iuk(x0k ,0).

2М —5 k

Используя замену переменной интегрирования и заменяя интеграл с помощью квадратурной формулы Гаусса суммой, найдем

1 + К Ж I k ^ .

— —--— E gk (tm ) = Vk (x0k ,0) — iuk (x0k ,0) ,

2и M m=1 k k

где mч - число узловых точек в отрезке (—5 k, x0 ).

Из последнего соотношения находим:

Л Mu ли Mu

1 + к ж5 k ^ 1 + к ж5 k . .

Vk (x0k ,0) = — ----k E Vk (tm ) , uk (x0k ,0) = — ----k E uk (tm )

k м m=1 k м m=

и для модуля вектора перемещений на берегах зон предразрушения при

xk = хок получаем:

v0t Чv++-v-)2 + (u¡-uk)2 = ij+ B¡

к 2 u M

где Ак (^); Вк ик (^) (к = 1, 2, ..., Ы).

т=1 т=1

Таким образом, условием, определяющим критическое значение объемных сил, при котором происходит появление трещин, будет

1 + к л1 к

лДТ+BÍ = . (14)

2j M

Совместное решение нелинейной системы уравнений, состоящей из алгебраических систем, заменяющих интегральные уравнения с условиями разрешимости задачи и критерия (14) позволяет определить размер зоны предразрушения, значения искомых функций vk (хк) и uk (хк) в узловых точках, величину внешней нагрузки, при которой происходит появление трещины в диске.

Для упрощения вычислений функции X(х, y) и Y(х, y) разлагались в ряды Тейлора, причем ограничивались несколькими первыми членами разложения. С помощью интегрирования уравнений (5) находим

z z

Q(z, z) = JdzJF(z, z)dz ;

zz

F1( z, z) = J dz J F (z, z) dz

По найденным функциям Q(z, z) и F1 (z, z) согласно (4) находим

_0 o

^r , ?re, ^, тчУ1 .

Заключение. Предложена эффективная схема расчета трещинообразования в круговом диске внешними объемными силами.

Л и т е р а т у р а

1. Болотин В.В. Механика зарождения трещин и начального развития усталостных трещин // ФХММ, 1986, т. 22, №1, с. 18-23.

2. Мирсалимов В.М. Зарождение дефекта типа трещины во втулке контактной пары // Математическое моделирование, 2005, том 17, №2, с. 35-45.

3. Ахмедова М.В. Зарождение трещин в тонкой пластине, ослабленной периодической системой криволинейных отверстий // Вестник ЧПГУ им. И.Я. Яковлева, серия: Механика предельного состояния, 2013, №4 (18), с. 3-14.

4. Искендеров Р.А. Зарождение трещины при поперечном изгибе изотропной пластины, ослабленной периодической системой круговых отверстий // Строительная механика инженерных конструкций и сооружений, 2013, №3, с. 18-28.

5. Гасанов Ф.Ф. Зарождение трещин в изотропной среде с периодической системой круговых отверстий, заполненных жесткими включениями, при продольном сдвиге // Строительная механика инженерных конструкций и сооружений, 2014, №3, с. 44-50.

6. Мир-Салим-заде М.В. Зарождение трещин в перфорированной подкрепленной пластине // Прикл. мех. и техн. физика, 2008, том 49, № 6, с. 170-180.

7. Zolgharnein E., Mirsalimov V.M. Nucleation of a Crack under Inner Compression of Cylindrical Bodies // ActaPolytechnicaHungarica, 2012, Vol. 9, № 2, p. 169-183.

8. Mirsalimov V.M., Hasanov Sh.G. Modeling of crack nucleation in covering on an elastic base // International Journal of Damage Mechanics, 2014, Vol. 23(3), p. 430-450.

9. Мирсалимов В.М. Неодномерные упругопластические задачи. - М.: Наука, 1987. 256 с.

10. Панасюк В.В. Механика квазихрупкого разрушения материалов. - Киев: Наук. думка, 1991, 416 с.

11. Rusinko A., Rusinko K. Plasticity and Creep of Metals. - Springer, Verlag Berlin Heidelberg, 2011, 436 с.

12. Mirsalimov M. V., Rustamov B. E. Simulation of Partial Closure of a Crack-Like Cavity with Cohesion between the Faces in an Isotropic Medium // Journal of Applied Mechanics and Technical Physics, 2013, № 6, p. 1021-1029.

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

14. Panasyuk V.V., Savruk M.P. and Datsyshyn A.P. A general method of solution of two-dimensional problems in the theory of cracks. // Eng. Fract. Mech, Vol. 9, №2, 1977, p. 481-497.

R e f e r e n c e s

1. Bolotin, V.V. (1986). Mechanics of the initiation and initial development of fatigue cracks. Soviet materials science, Vol. 22, Issue 1, p. 14-19.

2. Mirsalimov, V.M. (2005). Initiation of defects such as a crack in the bush of contact pair. Matem. Mod. Vol. 17, № 2, p. 35-45.

3. Akhmedova, M.V. (2013). Cracks nucleation in thin plate, weakened by the periodic system of the curvilinear holes. Vestnik of I. Yakovlev Chuvash State Pedagogical University. Line: Mechanics of definable state. № 4 (18), p. 3-14.

4. Iskenderov ,R.A. (2013). The crack nucleation in the isotropic plate, weakened by a periodical system of circular holes under transverse bending. Structural Mechanics of Engineering Constructions and Buildings, №3, P. 18-28.

5. Hasanov, F.F. (2014). Nucleation of cracks in isotropic medium with periodic system of the circular holes filled with rigid inclusions, at longitudinal shear. Structural Mechanics of Engineering Constructions and Buildings, №3, p. 44-50.

6. Mir-Salim-zade, M.V. (2008). Generation of cracks in a perforated reinforced plate. Journal of Applied Mechanics and Technical Physics, Vol. 49, Issue 6, p. 1030-1039.

7. Zolgharnein, E., Mirsalimov, V.M. (2012). Nucleation of a Crack under Inner Compression of Cylindrical Bodies. ActaPolytechnicaHungarica, Vol. 9, №2, p. 169-183.

8. Mirsalimov, V.M., Hasanov, Sh.G. (2014). Modeling of crack nucleation in covering on an elastic base. International Journal of Damage Mechanics, Vol. 23(3), p. 430-450.

9. Mirsalimov, V.M. (1987). Non-one-dimensional elastoplasticproblems. Moscow: Nauka. 256 p.

10. Panasyuk, V.V. (1991). Mechanics of quasibrittle fracture of material. Naukova Dumka, Kiev. 416 p.

11. Rusinko, A., Rusinko, K. (2011). Plasticity and Creep of Metals. Springer, Verlag Berlin Heidelberg, 436 p.

12. Mirsalimov, M.V., Rustamov, B.E. (2013). Simulation of Partial Closure of a Crack-Like Cavity with Cohesion between the Faces in an Isotropic Medium. Journal of Applied Mechanics and Technical Physics, №6, p. 1021-1029.

13. Muskhelishvili, N.I. (1977). Some Basic Problems of Mathematical Theory of Elasticity. Amsterdam: Kluwer, 707 p.

14. Panasyuk, V. V., Savruk, M.P. and Datsyshyn, A.P. (1977). A general method of solution of two-dimensional problems in the theory of cracks. Eng. Fract. Mech., Vol. 9, № 2, p. 481-497.

CRACKING IN CIRCULAR DISK UNDER THE INFLUENCE OF VOLUME FORCES

N.M. Kalantarly

Institute of Mathematics and Mechanics of NAS of Azerbaijan

A circular metal disk under the influence of volume forces is considered. The solution of fracture mechanics problem of crack nucleation in the circular disk under action of volume forces is given. The model of weakened interparticle bonds zone on condition of plastic flow under constant stress is used.

KEY WORDS: circular disk, zone of weakened interparticle bonds of material, cracking, volume forces.