Научная статья на тему 'Моделирование трещинообразования в полосе переменной толщины'

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

CC BY
153
61
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОЛОСА ПЕРЕМЕННОЙ ТОЛЩИНЫ / ЗОНА ПРЕДРАЗРУШЕНИЯ / ЗАРОЖДЕНИЕ ТРЕЩИНЫ / STRIP OF VARIABLE THICKNESS / PREFRACTURE ZONE / NUCLEATION OF CRACKS

Аннотация научной статьи по физике, автор научной работы — Мирсалимов М. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Мирсалимов М. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Modelling of Cracking in the Strip of Variable Thickness

The mathematical description of the model crack nucleation in a strip of variable thickness was conducted. Definition of the unknown parameters that characterize the nucleus crack is reduced to solving a system of singular integral equations. The condition that determines the critical value of external load at which cracking occurs was obtained.

Текст научной работы на тему «Моделирование трещинообразования в полосе переменной толщины»

9. Додд Р., Эйлбек Дж., Гиббон Дж., Моррис Х. Со-литоны и нелинейные волновые уравнения. М. : Мир, 1988. 694 с.

10. Морозов А. Д., Драгун Т. Н. Визуализация и анализ инвариантных множеств динамических систем. М.; Ижевск : ИКИ, 2003. 304 с.

11. Goirgilli A., Lazutkin V. F., Simo C. Visualization of УДК 539.375

МОДЕЛИРОВАНИЕ ТРЕЩИНООБРАЗОВАНИЯ В ПОЛОСЕ ПЕРЕМЕННОЙ ТОЛЩИНЫ

М. В. Мирсалимов

Азербайджанский технический университет, Баку E-mail: [email protected]

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

Ключевые слова: полоса переменной толщины, зона предраз-рушения, зарождение трещины.

a hyperbolic structure in area preserving maps // Reg. & Chaot. Dyn. 1997. Vol. 2, № 3/4. P. 47-61.

12. Асланов В. С., Ледков А. С. Особенности вращательного движения космического аппарата при спуске в атмосфере Марса // Космические исследования. 2007. Т. 45, № 4. С. 351-357.

Modelling of Cracking in the Strip of Variable Thickness M. V. Mirsalimov

The mathematical description of the model crack nucleation in a strip of variable thickness was conducted. Definition of the unknown parameters that characterize the nucleus crack is reduced to solving a system of singular integral equations. The condition that determines the critical value of external load at which cracking occurs was obtained.

Key words: strip of variable thickness, prefracture zone, nucleation of cracks.

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

ПОСТАНОВКА ЗАДАЧИ

Рассматривается однородная изотропная полоса переменной толщины. Обозначим через 2Ь и 2Н ширину и толщину полосы соответственно (рис. 1). Срединная плоскость (х, у) является плоскостью симметрии. Полоса переменной толщины находится в обобщенном плоском напряженном состоянии. Пусть полоса подвергается силовому нагружению. Считается, что заданы напряжения на границе у = ±Ь полосы (условие первой основной задачи). Принято, что в материале полосы после многократного нагружения появляется концентратор напряжений (область ослабленных межчастичных связей материала). При нагружении в прослойке перенапряженного материала образуется зона пластического течения. Пусть для определенности силовая нагрузка изменяется так, что в области ослабленных межчастичных связей материала осуществляется пластическое деформирование. Через некоторое число нагружений возможность пластического деформирования в области ослабленных межчастичных связей материала исчерпывается, и раскрытие берегов зоны пластического течения резко увеличивается. Если раскрытие берегов зоны предразрушения в точке максимальной концентрации достигает предельного значения 5с для данного материала полосы, то в этой точке зарождается трещина. Важно отметить, что в предлагаемой расчетной модели удлинение зародышевой трещины происходит непрерывно от нулевой длины.

Рис. 1. Расчетная схема задачи о зарождении трещины в полосе переменной толщины

© Мирсалимов М. В., 2012

95

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

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

Твердое тело (полоса), деформируемое за пределом упругости, представляется как тело, деформируемое упруго всюду, кроме некоторых поверхностей зон предразрушения. Напомним, что упругая деформация представляет собой изменение в расстояниях между элементарными частицами без изменения порядка их расположения, в то время как пластическая деформация сопровождается изменением относительного расположения атомов за счет сдвигов одного слоя атомов относительно другого (элементарные пластические деформации). Поскольку указанные зоны (прослойки перенапряженного материала) малы по сравнению с остальной упругой частью полосы, их можно мысленно заменить разрезами, поверхности которых взаимодействуют между собой по некоторому закону, соответствующему действию удаленного материала. Таким образом, задача, о деформации полосы переменной толщины за пределами упругости сводится к задаче о деформации некоторого упругого тела, имеющего условные разрезы. Силы, действующие по поверхностям этих зон, называются силами ослабленных связей, а зоны, в которых они появляются, — зонами предразрушения, или областями ослабленных связей. Считается, что толщина полосы 2Н(х,у) удовлетворяет условиям 0 < Н1 < Н(х,у) < Н2, где Н1, Н2 — наименьшее и наибольшее значение толщины полосы соответственно. Функция толщины может быть представлена в виде

Н(х, у) = Но [1 + еН(х, у)] ,

где Н0 = (Н1 + Н2)/2, е = (Н2 — Н1 )/(Н2 + Н1), —1 < Н(х,у) < 1 — некоторая известная безразмерная непрерывная функция.

Зона предразрушения ориентирована в направлении максимальных растягивающих напряжений. В центре зоны предразрушения разместим начала локальной системы координат х101 у1, ось х1 которой совпадает с линией расположения зоны предразрушения и образует угол а1 с осью х (см. рис. 1). Взаимодействие берегов зоны предразрушения (связи между берегами) сдерживает зарождение трещины. При действии силовой нагрузки на границе полосы у = ±Ь в связях, соединяющих берега зоны предразрушения, возникают нормальные аУ1 = а3 и касательные тХ1У1 = т3 напряжения, где а3 — предел текучести материала на растяжение, т3 — предел текучести материала полосы на сдвиг.

Таким образом, к берегам зоны ослабленных межчастичных связей материала приложены нормальные а3 и касательные т3 напряжения.

Запишем общие уравнения статического деформирования полосы:

- уравнения равновесия

дЫх + ту = 0 + Шху = 0

дх ду ду дх

закон Гука

N 2ЕН (ди ду\ ^ 2ЕН (дь ди\ м ЕН (ди ду

х 1 — V2 удх + ^ ду ) ' у 1 — V2 у ду + ^ дх) ' ху 1 + V \ду + дх

где Ых, Ыу, Ыху — нормальные и сдвигающие усилия, приходящиеся на единицу длины соответственно; и, у — компоненты вектора смещений; E — модуль упругости материала полосы; V — коэффициент Пуассона материала.

Краевые условия задачи имеют вид

при у = = Ь ау - *тху Р1(х), — 8 < х < го.

при у = = —Ь ау - *тх у = Р2(х). — 8 < х < го.

при у1 : = 0 Су1 — *тх 1 у1 = — ■ , — < х1 <

(1)

где Су, тху — компоненты тензора напряжений; р (х) и р2(х) — заданные функции, убывающие на бесконечности и удовлетворяющие условиям статики.

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

|-и+ — V- — г(и+ — и-)| = 5С, (2)

где <5С — характеристика сопротивления материала полосы трещинообразованию. Это дополнительное условие позволяет найти параметры полосы, при которых появляется трещина.

МЕТОД РЕШЕНИЯ

Решение системы уравнений статического деформирования полосы ищем в виде

N = №> + ежХ1} +..., N = №0> + е^у1) +..., Жху = №> + еж(у> +...,

и = и0 + еи1 + . . . , V = V0 + еv1 + . . . , 11 = I0 + е*1 + ...

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

дй

дй

дх ' ху ду

у = ж(0> — + №> —

1 у ду + ху дх.

(3)

Аналогично определяются составляющие Х2, У2 объемной силы для второго и последующих приближений.

Граничные условия задачи (1) принимают соответственно вид:

- для нулевого приближения

при

при

при

У = Ь

V = —Ь

VI =0

г(0>

у

г(0>

у

(0)

(0)

— гтх у = Р1 (х),

(0)

— 2тх у = Р2 (х),

— 2Т-

у

(0) = х1 у1 = а

—го < х < го. —го < х < го. —I0 < х1 < I0;

(4)

для первого приближения

при II Ж* у — ^у

при II — Ж* у — ^у

при 0 = у1 Ж* у — ^у

где р — гр2 = 2Лю(р1(х) — «Р2(х)), С1 — «т1 =

—го < х < го. —го < х < го. т1), —11 < х1

Здесь при выводе уравнений первого приближения были приняты следующие обозначения:

N = — й(х, у)Ж(0>, № = №1> — Л,(х, у)№0>,

^у = №> — Л(х,у)№>.

(5)

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

*

*

приложены усилия, равные по величине и обратные по знаку, определяемые первым напряженным состоянием при У1 = 0.

Рассмотрим нулевое приближение. Первое состояние для нулевого приближения (для полосы постоянной толщины) при отсутствии зоны предразрушения известно [2, 3].

Граничные условия задачи, определяющей второе состояние в нулевом приближении, принимает

вид

—го < Х < го,

—го < X < го. (6)

гтЦ^ = ст* — гт* — (р(0) (Х1) — гд(0) (Х1)), —I0 < XI < I0.

:»(0) (х1)и д; ' (х1) — нормальные и касательные напряжения, возникающие в фектной полосе на оси х1 (у1 = 0) от действия напряжений на границе у = ±Ьполосы в нулевом приближении.

Граничные условия задачи (6) с помощью формул Колосова-Мусхелишвили [4] запишем в виде краевой задачи для отыскания двух аналитических функций Ф0(г) и Ф0(г):

при у=Ь ст(0) ст° = 0,

при у II — ст(0) ст° = 0,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

при у1 = 0 ст(0) ст°1 • (0) гТХ 1°1

(х1)и д(0)(х1) — нормальные

т(°) = 0 т X ° — 0.

т(°) = 0 т ° = 0,

при у = Ь Ф0 (г) + Ф0 (г) + гФ0(г ) + Ф0(г)=0,

при у = —Ь Ф0 (г) + Ф0 (г) + гФ0(г ) + Ф0(г)=0,

при у1 = 0 Ф0 (Х1) + Ф0 (Х1) + Х1Ф0 (Х1) + Ф0 (Х1) = / 0(Х1),

где /0(Х1) = ст* — гт* — (р(0) — гд(0)), |Х11 < I0.

Комплексные потенциалы Ф0(г) и Ф0(г) ищем в виде [3]

ф0 м = 2л £/

^=0 , _/к

ф0(г) = ^Ее"2'ак /

&=0 , —1к

#0 (Ь) т* в1

г0°(Ь)

(7)

А,

(8)

где = Ьвг“к + ^0, = в-г“к (г — г^), а0 = а2 = 0, = гЬ, = —гЬ, 10 = го, 12 = го.

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

#о(Ь) Ь — Х

#2 (Ь)

+ #2 (Ь)К0,2 (Ь — Х) + 0° (Ь)Ь0,2 (Ь — х)

+ 0о(£)К2,0 (Ь — Х) + #°(Ь)Ь2,0(Ь — Х)

#0(Ь) Ь — Х

А = — J ^0(Ь)К0>1(Ь, х) + #0(Ь)Ь0,1 (Ь, Х^ А, |х| < го,

-г?

(9)

г?

А = —J [#0(Ь)К2>1 (Ь,х) + #0(Ь)Ь2>1 (Ь,х^ ЙЬ, |х| < го,

_/0

г1

+ #1(Ь)К1,1 (Ь,х) + #0(Ь)Ь1,1 (Ь,х)

ЙЬ + / #о(Ь)К1,0(Ь,х)+ (Ь)Ь1,0(Ь,х)

+ / [#2(Ь)К1,2(Ь,х)+ 0°(Ь)Ь1,2(Ь,х)] А = п/0(х), |х| < I0.

(10)

Ввиду некоторой громоздкости величины КП)£, ЬП)& (&,п = 0,1, 2) не приводятся. Из системы трех интегральных уравнений исключаются две неизвестные функции д°(Ь) и д° (£)• Подставляя в (10) функции д°(х)и д° (х), найденные из решения интегральных уравнений (9), после некоторых преобразований получим одно комплексное сингулярное интегральное уравнение относительно неизвестной функции д°(х):

#0(Ь) Ь — х

+ ^ (Ь)Яц (Ь, х) + #0 (Ь)$11 (Ь, х)

А = п/0(х), |х| < I0.

(11)

к

0

1

ОС

ОС

О

_/0

г1

О

_/0

г1

Выражения для функций Я 11(Ь, х) и 511(Ь, х) из-за некоторой громоздкости величины не приводятся (они имеют вид, подобный (У.41) из [3]).

К сингулярному интегральному уравнению (11) для внутренней зоны предразрушения добавляется дополнительное условие:

J д0 (Ь) ^ = о, (12)

-г?

обеспечивающее однозначность смещений при обходе контура зоны предразрушения в нулевом приближении. Комплексное сингулярное интегральное уравнение (11) при дополнительном условии (12) сводится к системе M алгебраических уравнений относительно приближенных значений искомой функции д0(х 1) в узловых точках. Для получения системы алгебраических уравнений сначала приведем в интегральном уравнении (11) и в условии (12) все интервалы интегрирования к одному отрезку [—1,1] с помощью замены переменных Ь = I0т, х = 10^ (|Ь| < I0, |х| < I0). Решение сингулярного интегрального уравнения ищем в виде

.0^_ 91 (О

g0«) =

где ^і(С) — регулярная в интервале [—1,1] функция.

Используя квадратурные формулы типа Гаусса [3, 5], сингулярное интегральное уравнение (11) при дополнительном условии (12) сведем к системе М алгебраических уравнений для определения М неизвестных д|(£т) (т = 1, 2,..., М):

І M

l 0 [g*(tm)Rl 1 (l 0 tm ,l 0Xr )+ gl(tm )ST i (l 0 tm ,l 0 Xr )] = /° (Xr ), m=1 (13)

M

Yh gl (tm) = 0 (r = І, 2,..., M — І; m = І, 2,..., M),

m=

2m — І nr

где tm = cos ——tt-n, xr = cos —.

2M M

Так как решение сингулярного интегрального уравнения (11) ищется в классе всюду ограниченных функций (напряжений), то к системе (13) необходимо добавить условия ограниченности напряжений на концах зоны предразрушения xr = ±l 0 в нулевом приближении. Эти условия имеют вид

M 2m І M 2m І

^(_ ^M+mgl(tm)tg n = 0, J](— ^g^m) ctg ш n = 0. (14)

m= m=

Полученная алгебраическая система (13)-(14) из-за неизвестного размера зоны предразрушения нелинейная. Система уравнений (13)—(14) относительно д*(Ьт) и I0 позволяет при заданной внешней нагрузке найти напряженно-деформированное состояние полосы при наличии зоны предразрушения в нулевом приближении. Решение системы уравнений (13)—(14) позволяет найти приближенные значения функций г>0(Ьт), и0(Ьт)и размеры зоны предразрушения. Для ее решения используется метод последовательных приближений, суть которого состоит в следующем. Решаем систему (13) при некотором определенном значении размера зоны предразрушения относительно остальных неизвестных. Остальные неизвестные входят в систему линейным образом. Значения размера зоны предразрушения и соответствующие значения остальных неизвестных не будут, вообще говоря, удовлетворять уравнениям (14). Подбирая значения параметров длины зоны предразрушения, будем многократно повторять вычисления, пока уравнения (14) не будут удовлетворяться с достаточной точностью. В каждом приближении алгебраическая система решалась методом Гаусса с выбором главного элемента.

Перейдем к построению решения задачи в первом приближении. Зная решение задачи в нулевом приближении, находим компоненты усилий Ж(0), ^У0), По формулам (3) находим компоненты

объемной силы ^ = X + гУ1 в первом приближении. Для определения первого состояния в первом приближении имеем уравнения плоской теории упругости при наличии объемных сил (3). При

наличии объемных сил решение плоской задачи для первого состояния ищем в виде N * = N(1) + N(1) N * = N(1) + N(1) N * = N(1) + N(1)

х х* 1 х? ’ у у* У1 5 х У х У* х У1 '

где N(*), жУ1), ^у* — частное решение уравнений плоской теории упругости при наличии объемной силы (3); , ^^, N^1 — общее решение уравнений плоской теории упругости при отсутствии

объемных сил.

Для усилий первого состояния в первом приближении при наличии объемных сил имеем общие представления

N + N _ / ^ 1 др

х у = 4 Ие ( Ф1 0 (г) - 1

2Л0 у 2(1 + к) дг

_ _ (15)

N - ^ +2*^*у 2 М.-ГМ, 1 д(к^1 - а^

——— -2 ^ ®1о(г)+ф1о(г) + 2(Г+^) —дг—; •

В эти соотношения входят две аналитические функции Ф10(г) и Г10(г) комплексной переменной г = х+гу и две функции р (г, г), ^1 (г, г), представляющие собой любые частные решения уравнений

д2 Р (г, г) = д 2^1 (г, г) = р

дгдг , дг2

Граничные условия для нахождения комплексных потенциалов Ф10 (г) и Г10 (г) на основании первых двух краевых условий (5) и общих представлений (15) имеют вид

при у = Ь Ф10(-) + Фю(-)+ -Ф'ю(г) + Фю(-) = Л(г) - Лр*(х),

при у = -Ь Ф10(-) + Фю(-)+ -Ф'ю(г) + Фю(-) = Д(г) - Лр*(х),

где

(16)

/ ( ) 1 И дР 1 / дР дС1 1 (17)

Л (2)-TT^ Ие эТ - 2(ГГ^У ^‘дТ - ~аГ^ (17)

Здесь функции р (г, г) и ^1(г, г) определяются по формулам

Р (г, г) = у йг у Р(-,-)й-, ^(г, г) = у йг у Р(г,-)й— (18)

Решение граничной задачи (16) известно [2, 3].

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Согласно принципу суперпозиции граничные условия задачи, определяющие второе напряженно-деформированное состояние в первом приближении, получим в виде

при У = Ь N = 0, ^у = o,

при у = -Ь N = 0, Х*у = 0, (19)

при у1 = 0, |х11 < 1} Ху? = -Ла* - р*(1), Я*?у? = -Л/г* - д* (1).

Здесь р*(1)(х1) = 2Л0р(1) и д* (1)(х1) = 2Л0д(1) — нормальные и касательные напряжения, возникающие в сплошной бездефектной полосе без зоны предразрушения по оси х1 (у1=0) в первом состоянии первого приближения соответственно.

Для определения комплексных потенциалов Ф1(г) и Г1(г) второго состояния первого приближения имеем граничную задачу

при у = Ь Ф1 (-) + Ф1(-)+ -Ф1 (г) + Г1 (г) = 0,

при у = -Ь Ф1 (г) + Ф1(г)+ -Ф1 (г) + Г1 (г) = 0, (20)

при у1 =0, |х11< 1} Ф1 (х1) + Ф1(х1)+ х1 Ф1(х1 ) + Г1(х1) = / (х1),

где / 1(х1) = -Л(а5 - гт*) - (р(1) - гд(1))

Комплексные потенциалы Ф^г) и Ф1 (г) ищем в виде, подобном (8). Поступая аналогично вышеизложенному, получим комплексное сингулярное интегральное уравнение относительно искомой функции ^1(х1) с дополнительным условием типа (12).

Как и в нулевом приближении, интегральное уравнение с дополнительным условием типа (12) сводим к системе комплексных алгебраических уравнений относительно M неизвестных д*(1) (£ш) (т = 1, 2,..., М)

1 м _________

М ^ I1 [^1 (1) (*ш)Я*1 (I1 *ш ,1^ )+ 0* (1)(^ш £ш ,11 Хг )] = /1 (жг ),

т=1 м (21)

Х>*(1 (^ш) = 0 (г = 1,2,..., М - 1).

т=1

Для замкнутости алгебраической системы (21) не хватает двух уравнений, определяющих размеры зоны предразрушения в первом приближении. Записывая условия конечности напряжений в вершинах зоны предразрушения х1 = ±1}, получим недостающие два уравнения

м 2т 1 м 2т 1

]Г(-1)шд*(1)^ш^ га п = 0, ]Г(-1)М+шд*(1)^ш^ Т, п = 0. (22)

Для определения раскрытия берегов зоны предразрушения и размеров зоны предразрушения задается закон изменения толщины полосы. Функцию Л(ж, у) разложим в ряд Тейлора в окрестности начала координат, ограничиваясь для простоты несколькими первыми членами разложения. Проинтегрировав (18), находим функции ^1 (г, г) и ^1 (г, г). По найденным функциям ^1 (г, г) и ^1(г,г) находим функцию /1(г). О способе решения алгебраической системы (21), (22) сказано выше. После решения алгебраической системы находим

01 (Ж1) = д0 (х1) + е^ (Х1), 11 = 10 + е^.

Для определения предельного значения внешней нагрузки, при которой появится трещина, используем соотношение (2).

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

Х1

1 | к С

д1(х1)^х1 = V*(х1, 0) — -ш *(х1? 0).

2д J

-її

где V* (ж, 0) = -и+ — V-, и * (ж, 0) = и+ — и-, д — модуль сдвига материала.

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

X 0

1 + К [

----2— ді(жі) ^жі = V*(жо, 0) — іи *(жо, 0).

-її

Применяя замену переменных и заменяя интеграл суммой, находим

1 . / М!

----ді(*т) = V*(жо, 0) — іи*(жо, 0), (23)

т=і

где Мі — число узловых точек, содержащихся в отрезке (—7і,ж0).

Из (23) находим

1 + 7 Мї 1 + 7 Мї

V* (ж0, 0) =-----2^- ^ (^^ и * (ж0, 0) =---2ГГ иі (^т).

т=і т=і

Тогда модуль вектора смещений на берегах зоны предразрушения при жі = ж0 будет

. + к 7г/і 2д М

V:. = Vи2 + «2 = 1+-К а2 + в2, (24)

Мї Мї

где А = £ -иі(£т), в = £ V!(£т).

т=і т=і

Таким образом, условием, определяющим предельное значение внешней нагрузки, будет

1 + к nl1 2д M

л/А2 + B2 = 5С.

(25)

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

На рис. 2 представлена зависимость длины зоны предразрушения ^/2Ь от безразмерной нагрузки ро/о^ для полосы, толщина которой меняется по линейному закону. В расчетах были приняты следующие значения параметров: V = 0.3; Е =1.8 ■ 105М Па, а1 = п/4, г° =0, М = 30.

р°л/£

На рис. 3 приведена зависимость безразмерной предельной нагрузки pc =

от относитель-

ной толщины полосы Л,1 /Л0.

Если зона предразрушения одним концом выходит на наружную поверхность полосы (у = ±£), необходимость в дополнительном условии (12) в каждом приближении отпадает.

Рис. 2. Зависимость длины зоны предразрушения li/2L от безразмерной нагрузки po/ Us

ЗАКЛЮЧЕНИЕ

yoVL \/Eas Sc

от относительной толщины полосы

hi

ho

Анализ модели зарождения трещины в полосе переменной толщины при силовом нагружении сводится к параметрическому совместному исследованию разрешающей алгебраической системы задачи и критерия появления трещины (25) при различных значениях свободных параметров полосы. Это различные механические характеристики материала и геометрические данные полосы. Построенная расчетная модель позволяет при изменении параметров а1 и г° исследовать влияние расположения зоны ослабленных связей материала в полосе на трещинообразование.

В заключение отметим, что построенная модель, описывающая начало разрушения материала полосы, основана на дислокационных механизмах зарождения трещины. Дислокационные механизмы зарождения трещины получили экспериментальное подтверждение [1].

Библиографический список

1. Панасюк В. В. Механика квазихрупкого разрушения нах и оболочках. Киев : Наук. думка, 1976. 443 с.

материалов. Киев : Наук. думка, 1991. 416 с. 4. Мусхелишвили Н. И. Некоторые основные задачи

2. Уфлянд Я. С. Интегральные преобразования в зада- математической теории упругости. М. : Наука, 1966.

чах теории упругости. Л. : Наука, 1967. 420 с. 707 с.

3. Панасюк В. В., Саврук М. П., Дацышин А. П. 5. Мирсалимов В. М. Неодномерные упругопластиче-Распределение напряжений около трещин в пласти- ские задачи. М. : Наука, 1987. 256 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.