Научная статья на тему 'Асимптотические методы нелинейной механики разрушения: результаты, современное состояние и перспективы'

Асимптотические методы нелинейной механики разрушения: результаты, современное состояние и перспективы Текст научной статьи по специальности «Физика»

CC BY
441
124
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АСИМПТОТИЧЕСКИЕ МЕТОДЫ / НЕЛИНЕЙНАЯ МЕХАНИКА РАЗРУШЕНИЯ / ПОЛЕ НАПРЯЖЕНИЙ ВБЛИЗИ КОНЧИКА ТРЕЩИНЫ / НЕЛИНЕЙНЫЕ ЗАДАЧИ НА СОБСТВЕННЫЕ ЗНАЧЕНИЯ / ASYMPTOTIC METHODS / NONLINEAR FRACTURE MECHANICS / CRACK TIP STRESS FIELD / NONLINEAR EIGENVALUE PROBLEMS

Аннотация научной статьи по физике, автор научной работы — Степанова Лариса Валентиновна, Адылина Екатерина Михайловна

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

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

Asymptotic methods of nonlinear fracture mechanics: results, contemporary state and perspectives

In the paper the brief review of the important results of nonlinear fracture mechanics recently obtained by the asymptotic methods and perturbation techniques is given. The asymptotic solutions of the whole class of crack problems are discussed. Applications of the asymptotic analysis for crack problems of nonlinear fracture mechanics such that crack tip fields in damaged medium, fatigue crack growth problems in damaged medium, nonlinear eigenvalue problems arising from fracture mechanics analysis, nonlinear eigenvalue analysis including higher order terms are elucidated.

Текст научной работы на тему «Асимптотические методы нелинейной механики разрушения: результаты, современное состояние и перспективы»

УДК 539.42

АСИМПТОТИЧЕСКИЕ МЕТОДЫ НЕЛИНЕЙНОЙ МЕХАНИКИ РАЗРУШЕНИЯ: РЕЗУЛЬТАТЫ, СОВРЕМЕННОЕ СОСТОЯНИЕ И ПЕРСПЕКТИВЫ

Л. В. Степанова, Е. М. Адылина

Самарский государственный университет,

Россия, 443011, Самара, ул. Академика Павлова, 1.

E-mails: lst@ssu. samara.ru, kateadulina@mail .ru

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

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

1. Асимптотический анализ в линейной механике разрушения. Асимптотический анализ распределений напряжений, деформаций и перемещений вблизи вершины трещины является одной из фундаментальных задач механики трещин [1]. Вершина трещины может находиться в упругом материале, пластической зоне или зоне, занятой деформациями ползучести; кончик трещины может принадлежать линии раздела нескольких изотропных или анизотропных материалов [2]. Все перечисленные факторы усложняют определение механических полей вблизи устья трещины, и до настоящего времени многие вопросы, связанные с нахождением напряжённо-деформированного состояния в нелинейной механике разрушения, остаются открытыми. Так, в настоящее время в механике трещин и в целом в механике разрушения сложилось понимание процесса разрушения как процесса многомасштабного и многоуровневого, для описания основных закономерностей которого следует прибегать к многомасштабным моделям [3-5]. В рамках многоуровневого подхода процесс разрушения моделируется различными соотношениями на различных расстояниях от вершины трещины с помощью введения многих масштабов. При построении решения задачи в целом представления, работающие на разных расстояниях от кончика трещины, сращиваются в зонах, где справедливы асимптотики соседних областей. Для понимания многоуровневого или многомасштабного подхода целесообразно обратиться к исходным математическим моделям, используемым как в теории механики трещин, так

Лариса Валентиновна Степанова (д.ф.-м.н., доц.), профессор, каф. математического моделирования в механике. Екатерина Михайловна Адылина, аспирант, каф. математического моделирования в механике.

и в инженерной практике. Одним из наиболее распространенных математических методов построения распределений напряжений, деформаций и перемещений у кончика трещин является асимптотический анализ и подходы, развитые в асимптотической теории [6,7]. В рамках линейной теории упругости с помощью методов теории возмущений построено большое количество решений, ставших классическими в механике хрупкого разрушения. Поэтому, обращаясь к современным методам оценки напряжённо-деформированно-го состояния у вершины трещины в материалах с нелинейными определяющими уравнениями, естественно сначала обратиться к результатам асимптотического анализа линейной теории упругости, ибо результаты, полученные посредством асимптотических методов в линейной теории упругости, служат теоретическим основанием инженерной практики.

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

По всей видимости, первый пример сингулярности поля напряжений в линейной теории упругости дает задача о сосредоточенной силе, действующей на прямолинейную границу бесконечно большой линейной упругой пластины (задача Фламана) [8,9]. Существует фундаментальное решение задачи Фла-мана, называемое простым радиальным распределением напряжений. Любой элемент С, расположенный на расстоянии г от точки приложения сосредоточенной силы, подвергается простому сжатию в радиальном направлении. Компоненты тензора напряжений и поле перемещений определяются при этом выражениями

(jrr = — 2Pcos9/(irr), ивв = <тГ0 = 0, ur(9 = 0) = (2Р/(ттЕ))\пг + В. (1)

Важной особенностью решения задачи Фламана является сингулярность как поля напряжений, так и поля деформаций в точке приложения нагрузки. Более того, эта особенность поля напряжений и перемещений ведет к неин-тегрируемой особенности энергии деформации:

( Гж/2 fR a2 \ (Р2 R\

W = lim / / -^-rdr =lim( — In—].

e-S-0\J-7T/2 Je 2E J s^0\ttE eJ

Поле напряжений (1) обладает особенностью вида 1 /г при г —> 0. Математически решение (1) краевой задачи линейной теории упругости о сосредоточенной силе, действующей на прямолинейной границе, приводит к парадоксу, поскольку это решение удовлетворяет всем фундаментальным уравнениям линейной теории упругости, однако в то же время нарушает упрощающие предположения, принимаемые в линейной теории упругости: напряжения не

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

Другая особенность задачи Фламана заключается в том, что ее решение послужило основой для дальнейшего развития асимптотического анализа полей напряжений в линейно упругих телах. В 1907 г. К. Вигхардт впервые обратился к задаче определения поля напряжений вблизи углового выреза, нагруженного двумя сосредоточенными силами. Именно К. Вигхардт, решая обобщенную задачу Фламана, обнаружил факторизацию поля напряжений: поле напряжений может быть представлено в виде функции, представляющей собой произведение двух функций, одна из которых зависит от г, тогда как другая является функцией полярного угла 9. Для случая трещины К. Вигхардт обнаружил корневую особенность поля напряжений и начал дискуссию, касающуюся обоснования сингулярности поля напряжений у вершины трещины в линейно упругом материале.

В 1933 г. Дж. Братц опубликовал статью [10], в которой найдено поле напряжений у вершины острого выреза. В то время, когда работа К. Вигхарда была полностью забыта, Дж. Братц воспользовался идеей фон Кармана о введении функции напряжения Эри. Дж. Братц продемонстрировал возможность определения поля напряжений в упругом полупространстве (в двумерной постановке задачи) с общим распределением нормальных и касательных усилий, заданных на границе упругой области. Дж. Братц обобщил подход фон Кармана, показав, что метод может быть применим к любой клиновидной области, нагруженной нормальными и касательными усилиями на прямолинейных границах. В 1933 г. Вестергаард опубликовал работу, посвященную анализу напряжений у вершины трещины в бетонных элементах конструкций в случае изгиба. В его статье, являющейся, по всей видимости, одной из первых попыток использования механики разрушения для описания процесса разрушения квазихрупких материалов, было предложено следующее выражение для функции напряжений Эри:

Несмотря на упомянутые пионерские работы в области применения асимптотических методов в линейной механике разрушения, в полной мере возможность разделения переменных в решении задач для клиновидных областей была раскрыта М. Уильямсом, предложившим метод разложения по собственным функциям функции напряжений Эри, компонент тензора напряжений и вектора перемещений. Именно в его работах [11, 12] функция напряжений Эри была представлена в форме произведения двух функций, одна из которых зависела от радиальной переменной г (г — расстояние от вершины клина до рассматриваемой точки), а вторая — от полярного угла 9:

где и /]{9) - собственные значения и соответствующие им собственные

I) (Ш_ \) й1п(т + \)в~ (т + I) 8[п(т~ \)

9 .

(2)

функции, определяемые из решения задач на собственные значения. Сравнительная простота и легкость использования представления (2) при решении краевых задач механики трещин были важным аспектом, который обусловил широкое распространение этого подхода [7]. Следует отметить, что, с другой стороны, природа разделения переменных в решении задач, сводящихся к бигармоническому уравнению, была уже установлена в механике линейновязкой жидкости при изучении двумерного течения Стокса — течения вязкой несжимаемой жидкости, которое может быть описано посредством функции тока Ф, удовлетворяющей бигармоническому уравнению. В механике жидкости и газа применительно к бигармоническому уравнению относительно функции тока метод разложения по собственным функциям был предложен Дином и Монтагоном [13]. Они приложили метод разложения по собственным функциям к задачам о движении линейно-вязкой жидкости вблизи угловых областей, что с математической точки зрения аналогично решению задач определения напряжённо-деформированного состояния вблизи угловых вырезов в линейно-упругих материалах. Упомянутые работы, посвященные построению асимптотик механических полей близи трещин и угловых вырезов, обеспечили основу для ответа на фундаментальные вопросы механики трещин: при каких условиях трещина получает возможность распространяться и продолжает расти? К. Вигхардт в своей работе сформулировал следующие вопросы: при каком значении сосредоточенной силы Р для материала с заданными упругими и прочностными свойствами начнется разрушение материала, в каком месте и в каком направлении будет развиваться разрушение?

В современной нелинейной механике разрушения одним из наиболее распространенных методов анализа напряжённо-деформированного состояния вблизи вершины трещины как в линейно-упругих материалах, так и в материалах с нелинейными определяющими уравнениями, является метод разложения по собственным функциям, восходящий к работам М. Уильямса [11,12], в которых впервые было использовано разложение функции напряжений Эри в ряд по степеням расстояния от кончика трещины или вершины углового выреза в линейно-упругом материале (разложение по собственным функциям). В ставших классическими работах Дж. Райса, Дж. Розенгрена и Дж. Хатчинсона [14-16] асимптотическое представление решения задачи определения напряжённо-деформированного состояния в окрестности вершины трещины в материале с определяющими уравнениями Рамберга—Осгуда е = a/E-\-aan разыскивалось в форме

оф, в) = Кг-1/^агз{в), <7е(г, в) = Кг~1/^ае{в),

е?-(г, в) = аКпг~п^п+1Ч^{е), щ{г, в) = аКпг1^п+1^йг{9). { }

С тех пор асимптотический анализ распределения напряжений у кончика трещины стал неотъемлемой частью исследований, проводимых в механике разрушения [2,17]. Только в самое последнее время метод разложения по собственным функциям полей напряжений и перемещений у вершины трещины или углового выреза был применен в целом ряде исследований [18,19].

Асимптотическое представление решения в форме (3) широко используется в нелинейной механике разрушения для анализа полей напряжений, деформаций и перемещений в условиях смешанного нагружения [20-22]. Первое обращение к исследованию смешанных форм деформирования, по всей види-

мости, было сделано в работах Си [23,24], где рассмотрены поля напряжений и деформаций вблизи вершины трещины, находящейся под действием растягивающей и сдвиговой нагрузок (т.о. рассматривалось смешанное нагружение, отвечающее трещинам типа I и типа II). Си впервые ввёл коэффициент смешанности нагружения:

2

Мр = — arctg

»-—»о are(r, 9 = 0)

(4)

Параметр смешанности нагружения Мр принимает значение Мр = 0 для чистого поперечного сдвига и значение Мр = 1 для нормального отрыва;

О < Мр < 1 для всех промежуточных типов приложенной нагрузки.

В [23] приведены распределения напряжений, деформаций и перемещений у вершины трещины в материале, который подчиняется степенному закону деформационной теории пластичности для отдельных значений параметра смешанности нагружения в условиях плоской деформации. Там же могут быть найдены контуры областей пластического течения у вершины трещины для различных значений показателя нелинейности материала и параметра смешанности нагружения. Подробный анализ напряжённо-деформированно-го состояния у вершины наклонной трещины для степенного определяющего закона выполнен в [25] (подход, используемый авторами, также изложен в [26]). В этой работе авторами предлагается процедура решения краевой задачи для нелинейного дифференциального уравнения, получающейся из проблемы нахождения полей напряжений и деформаций у вершины трещины в условиях смешанного нагружения. В [27] собраны существенные результаты, полученные для смешанных форм деформирования к 2003 году.

Расчёты угловых распределений напряжений и деформаций у вершины трещины в материале со степенным определяющим законом для различных значений параметра смешанности нагружения для плоского деформированного состояния приведены в [22] в полном диапазоне смешанных форм деформирования — от нормального отрыва до чистого сдвига.

2. Смешанное нагружение элемента конструкции с трещиной (нормальный отрыв и поперечный сдвиг) в условиях плоского напряжённого состояния. Рассмотрим стационарную трещину в материале со степенным законом упрочнения. Введём прямоугольную декартову систему координат с центром в вершине трещины и полярную систему координат (г, 9) с полюсом в вершине трещины. Уравнения равновесия представляются в форме

Г(тгг,г + (тг0,$ + (arr — иве) = 0, (7вв,в + таг$,г + 2ar$ = 0.

Условие совместности деформаций формулируется как

2 (Г£гв,в),г = £гг,вв - ге„,г + Г (гевв);ГГ .

Определяющие уравнения e%j = З-Всг”-1 Sij/2 для случая плоского напряжённого состояния принимают форму

е„ = Во-™~1 (2(jrr - (гое)/2, еое = Ba™~1 (2aee - агг)/2, (5)

£rQ — ЗВае сггд/2, ае — crrr -|- Oqq огг<т$в 3<тг^.

Граничные условия задачи есть условия отсутствия поверхностных усилий на берегах трещины a00 (г, в = ±7г) = 0, (тгв(г, в = ±7г) = 0.

3. Асимптотическое представление компонент тензора напряжений и деформаций вблизи вершины трещины. Будем искать асимптотическое представление функции напряжений Эри в окрестности вершины трещины в форме F (г, в) = rx+2f(d). Тогда компоненты тензора напряжений принимают вид

<тоо(г, в) = г\А + 2)(Л + 1 )/(0), (Т„{г, в) = гх [(Л + 2)f(6) + f"(6)\ ,

агв(г,в) = -гх(Х + 1)Г(в).

Асимптотическое представление интенсивности напряжений в окрестности вершины трещины может быть представлено как

Veir, в) = rxfe{0), /g = [(Л + 2)/ + f"}2 + [(Л + 2)(Л + I)/]2 —

- [(Л + 2)/ + /"](Л + 2)(Л + 1)/ + 3(А + I)2/'2.

В силу соотношений (5) компоненты тензора деформаций в окрестности вершины трещины определяются выражениями

£„ = ВгХп£гг(в), £09 = BrXn£ee{6), £ге = ВгХп£гв{в), (6)

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

ётг = /г1 [(А+2)(1 —А)/+2/"] /2, ëoe = frl [(Л+2)(2А+1 )/-/"] /2,

ёгв(в) = -г%-1(\ + 1)Пв)/2.

Выражения (6) для компонент тензора деформаций позволяют получить нелинейное обыкновенное дифференциальное уравнение для функции f{9):

2(А n + 1 )ёгв,в = £гг,вв \nërr + (А п + 1)А пёвв

или

fIVfe {(« - !) [2f + s(3 - s)/]2/2 + 2/е } + 6(А + l)(An + 1 ){п- 1)/е2/'х X {{Sf + f")(sf + П + s2(s - 1 ?ff - S{s - 1 ){sf + П//2--s(s - 1 )(sf + /")/'/2 + 3(s - l)2/7"} +

+ 6(A + 1)(An + 1 )(n - l)/e4/" + (n - 1 )(n - 3)[2(sf + /") - s(s - 1 )/]x x {(sf + f”)(sf + П + s2(s - 1 )2ff - s(s - 1 )(sf + /"')//2--s(s - 1 )(sf + f")f/2 + 3(s - l)2/'/"}2 +

+ (n - 1 )/2 [2/" + s(3 - s)/] {(s/' + f'")2 + (s/ + /">/"+

+S2(S _ l)2(//2 + fn _ s2(s _ 1)/Г/2 _ ф _ 1)(s// + nf_

-s(s - l)(sf + Л/72 + 3(s - 1)2(///2 + f'f"')} +

+ 2(n - l)f2[2f”' + s(3 - s)/'] {(sf + /'Ois/' + n + s2(s - 1 )2//--ф - l)(s/' + Г0//2 - s(s - l)(sf + f")f'/2 + 3(s - l)2/7"} +

+ /e(3 - s)s/" - An/e4[2/" + s(3 - s)/] +

+ (An + l)An/g [s(2s - 3)/ - /"] =0. (7)

Решение уравнения (7) должно удовлетворять условиям, следующим из требований отсутствия поверхностных усилий на берегах трещины: /(0 = = ±7Г) = 0, /'(0 = ±7Г) = 0.

Вместе с этим необходимо учесть условия, задающие вид смешанного нагружения, которые формулируются на линии продолжения трещины (0 = 0). Из определения параметра смешанности нагружения (4) следует

/'(0 = 0) = -(А + 2)/(0 = О)^е(тгМ*/2),

что позволяет сформулировать начальную задачу для отрезка [0,7г]. Необходимо найти решение уравнения (7), удовлетворяющее условиям

/(0) = 1, //(0) = -(А + 2)^(тгМ72), /"(0 )=А2, Г(0)=А3,

где неизвестные на данном этапе построения решения константы А2 и А% подбираются таким образом, чтобы выполнялись краевые условия на берегу трещины: /(тг) = 0, /'(к) = 0.

Константы А2 и А% были найдены для всех значений показателя упрочнения материала п от 2 до 13 и для всех значений параметра смешанности нагружения 0 ^ Мр ^ 1. Затем осуществлялось построение решения уравнения (7) на отрезке [—7Г,0]. С целью построения решения на этом отрезке формулировалась начальная задача с условиями при 0 = —тт:

/(-тг) = 0, /'(-тг)=0, /"(-тг ) = В2, /"'(-*) = Яз-

Постоянные В2 и В% находятся из условий равновесия элемента, расположенного на луче 0 = 0, которые в терминах функции /(0) имеют вид

/-(0 = 0) = /+(0 = 0), (/')“ (0 = 0) = (/') + (0 = 0).

Константы В2 и В% были найдены для всех значений показателя упрочнения материала п от 2 до 13 и для всех значений параметра смешанности нагружения 0 < Мр < 1.

Угловые распределения компонент тензоров напряжений и деформаций в окрестности вершины трещины в условиях плоского напряжённого состояния для различных значений параметра смешанности нагружения для п = 5 показаны на рис. 1,2.

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

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

Рис. 1. Угловые распределения компонент тензора напряжений агг,агв для п = 5

л МР= 0.1 п—5

%в/ / 7 ^ V.

\ \ /

\&тв \

\ /\Дт

4

Рис. 2. Угловые распределения компонент тензора напряжений овв для п = 5

в материалах с нелинейными определяющими уравнениями. Для материалов, СЛеДуЮЩИХ СТепеННОМу Закону Теории установившейся ПОЛЗучеСТИ ¿ij = = З-Всг”_1Sij/2, асимптотические разложения компонент тензора напряжений в непосредственной окрестности вершины трещины (г —>■ 0) могут быть представлены в форме

4!>W + Лгг'Ч2’«') + A3r»aü>(ff) + .... (8)

где С* - интеграл, определяемый равенством

С* = J (^Wdx 2 - aijiij^ds^ , W = j aijdiij = ^¿e-

В литературе имеется достаточно большое число оценок инвариантного С*-интеграла, полученных для различных образцов, наиболее часто используемых в эксперименте. Однако остаются открытыми многие вопросы, связанные с определением второго коэффициента в асимптотических разложениях (8), а именно А2. Естественно предположить, что А2 в общем случае является функцией приложенных нагрузок, геометрических характеристик образца и показателя нелинейности материала п. Определение этого коэффициента и составляет задачу теории, поскольку даже известные оценки, полученные для первого амплитудного коэффициента и приведенные в справочной литературе [28], дают существенную погрешность. Например, для образца с односторонним боковым надрезом в [28] для инвариантного J-интеграла приведена формула

J = B(jo£oW(l — a)ah\ (P/Po),

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

где h\ — безразмерная функция приложенной нагрузки и геометрии образца. Оказывается, что уточненные вычисления, проведенные в последнее время, показали, что значения функции h\ могут отличаться от приведенных в [28] более чем на сто процентов. Следовательно, с появлением вычислительных комплексов, таких как SIMULIA Abaqus, предоставляется возможность получения более точных оценок для параметров механики разрушения.

5. Нелинейные задачи на собственные значения. Метод разложения по собственным функциям механических величин в механике разрушения приводит к нелинейным задачам на собственные значения, которые решаются, как правило, численно. Метод возмущений позволяет найти аналитическое решение нелинейных задач на собственные значения, возникающих в нелинейной механике разрушения. Метод возмущений для построения аналитического решения задачи определения напряжённо-деформированного состояния и поля сплошности у растущей в условиях циклического нагружения трещины был использован в [29], где показывается, что метод разложения по собственным функциям, когда функция напряжений Эри и параметр сплошности разыскиваются в форме F(r,6) = rxf(6), *ф(г,в) = rßg(6), редуцирует задачу к нелинейной задаче на собственные значения. Обычно в механике разрушения для построения решения нелинейных задач на собственные значения прибегают к численному интегрированию уравнений задачи, что может привести к

ошибочным результатам. Для преодоления указанных сложностей при решении нелинейной задачи на собственные значения в ходе применения метода возмущений построены четырехчленные асимптотические разложения искомых функций. В [29] получены аккуратные выражения для искомых функций вплоть до слагаемых порядка 0(е4), где е = /л — /хо — разность между собственными значениями, отвечающими нелинейной и линейной задачам:

П = 1 + ЕП\ + е2П2 + £3Пз + е4Щ + . . . , m = 1 + emi + e2m2 + е3тз + е4пг4 + ...,

/(0) = fo(0) + + £2/2(0) + e3fs(0) + ...,

0(0) = 9o(0) + £5-1(6») + e252(ö) + е35з(0) + - - - -

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

А = и = -----------.

1 + n — m

6. Направления дальнейшего развития асимптотических методов в нелинейной механике разрушения Одними из важных проблем являются разработка и развитие математических многомасштабных моделей накопления повреждений на макро-, мезо- и микроуровнях у вершины трещины в элементах конструкций, находящихся в условиях пластического течения, высокотемпературной ползучести, усталостного нагружения, а также формулировка критериев разрушения неупругих материалов, учитывающих процессы повреждаемости материала в процессе деформирования [3-5]. Многомасштабность или многоуровневость процесса разрушения можно отразить, вводя в рассмотрение поля напряжений на различных расстояниях от кончика трещины с разной особенностью (с различной степенью сингулярности). Таким образом, сильная и слабая особенности механических полей у вершины трещины отражают многомасштабный процесс разрушения. В силу этих особенностей моделирования процесса разрушения представляется важным знание всех возможных решений задачи определения напряжённо-деформированного состояния у вершины трещины, а не только тех слагаемых в асимптотических разложениях полей напряжений и перемещений, которые традиционно удерживаются в классических решениях. Таким образом, многоуровневые модели разрушения требуют исследования всего спектра собственных значений в задаче определения полей напряжений и перемещений в окрестности вершины трещины в материалах с нелинейными определяющими уравнениями. Если для упругих материалов в линейной механике разрушения построено полное разложение (решение Уильямса [11,12]), то для материалов со степенными определяющими соотношениями многие вопросы остаются открытыми.

Работа выполнена при поддержке РФФИ (проект № 12—08—00390—а).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. H. D. Bui, Fracture Mechanics: Inverse Problems and Solutions / Solid Mechanics and Its

Applications. Vol. 139. Dordrecht: Springer, 2006. xxiii+398 pp.; русск. пер.: X. Д. Бьюи,

Механика разрушения: Обратные задачи и решения / пер. с англ. J1. В. Степановой.

Москва: Физматлит, 2011. 412 с.

2. J. Ы, N. Recho, Méthodes asymptotiques en mécanique de la rupture. Paris: Hermès-

Lavoisier, 2002. 262 pp.

3. G. С. Sih, “Crack tip mechanics based on progressive damage of arrow: Hierarchy of singularities and multiscale segment” // J. Theoret. Appl. Fract. Mech., 2009. Vol. 51, no. 1. Pp. 11-32.

4. G. C. Sih, X. S. Tang, “Simultaneity of multiscaling for macro-meso-micro damage model represented by strong singularities” // J. Theoret. Appl. Fract. Mech., 2004. Vol. 42, no. 3. Pp. 199-225.

5. G. C. Sih, X. S. Tang, “Weak and strong singularities reflecting multiscale damage: microboundary conditions for free-free, fixed-fixed and free-fixed constraints” // J. Theoret. Appl. Fract. Mech., 2005. T. 43, №1. C. 5-62.

6. И. И. Аргатов, Введение в асимптотическое моделирование в механике. СПб.: Политехника, 2004. 302 с. [I. I. Argatov, Introduction to Asymptotic Modelling in Mechanics. St. Petersburg: Polytechnics, 2004. 302 pp.]

7. Л. В. Степанова, Математические методы механики разрушения. М.: Физматлит, 2009. 336 с. [L. V. Stepanova, Mathematical methods of fracture mechanics. Moscow: Fizmatlit, 2009. 336 pp.]

8. A. Flamant, “Sur la répartition des pressions dans un solide rectangulaire chargé transversalement” // Compte. Rendu. Acad. Sci. Paris, 1892. Vol. 114. Pp. 1465-1468.

9. S. P. Timoshenko, J. N. Goodier, Theory of elasticity. New York: McGraw-Hill, 1970. 567 pp.; русск. пер.: С. П. Тимошенко, Дж. Гудьер, Теория упругости. М.: Наука, 1979. 560 с.

10. J. Н. A. Brahtz, “Stress distribution in a reentrant corner” // Trans. Amer. Soc. Mech. Eng., 1933. Vol. 55. Pp. 31-71.

11. M. L. Williams, “Stress singularities resulting from various boundary conditions in angular corners of plates in extension” // J. Appl. Mech,., ASME, 1952. Vol. 74. Pp. 526-528.

12. M. L. Williams, “On the stress distribution at the base of a stationary crack” // J. Appl. Mech,., ASME, 1957. Vol. 24. Pp. 109-114.

13. W. R. Dean, P. E. Montagnon, “On the steady motion of viscous liquid in a corner” // Proc. Cambridge Philos. Soc., 1948. no. 45. Pp. 389-395.

14. J. M. Hutchinson, “Singular behaviour at the end of tensile crack in a hardening material” // J. Mech,. Phys. Solids, 1968. Vol. 16, no. 1. Pp. 13-31.

15. J. M. Hutchinson, “Plastic stress and strain fields at a crack tip” // J. Mech. Phys. Solids, 1968. Vol. 16, no. 5. Pp. 337-347.

16. J. R. Rice, G. F. Rosengren, “Plane strain deformation near a crack tip in a power-law hardening material” // J. Mech. Phys. Solids, 1968. Vol. 16, no. 1. Pp. 1-12.

17. A. Carpinteri, M. Paggi, “Asymptotic analysis in Linear Elasticity: From the pioneering studies by Wieghardt and Irwin until today” // Engn. Fract. Mech., 2009. Vol. 76, no. 12. Pp. 1771-1784.

18. Z. Niu, C. Cheng, N. Recho, “A new boundary element approach of modelling singular stress fields of plane V-notch problems” // Int. J. Solids Struct., 2009. Vol. 46, no. 16. Pp. 2999-3008.

19. G. Hello, М. B. Taha, J. M. Roelandt, “Analytical determination of coefficients in crack-tip stress expansions for a finite crack in an infinite plane medium” // Int. J. Solids Struct., 2012. Vol. 49, no. 3-4. Pp. 556-566.

20. В. H. Шляпников, “Смешанные моды развития трещин при сложном напряженном состоянии (обзор)”// Завод, лаб., 1990. Т. 56. С. 77-90. [V. N. Shlyannikov, “Mixed-modes crack growth under complex stress strain state (Review)”// Zavod. Lab., 1990. Vol. 56. Pp. 77-90].

21. В. H. Шляпников, Вычислительная механика деформирования и разрушения. Казань: Казан, гос. энерг. ун-т, 2002. 228 с. [V. N. Shlyannikov, Computational Mechanics of Deformation and Fracture. Kazan: Kazan State Power Engn. Univ., 2002. 228 pp.]

22. В. H. Шляпников, С. Ю. Кислова, “Параметры смешанных форм деформирования для трещины в виде математического разреза” // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика, 2009. Т. 9, №1. С. 77-84. [V. N. Shlyannikov,

S. Yu. Kislova, “Mode mixity parameters for mathematical crack type” // Izv. Saratov. Univ. Mat. Mekh. Inform,., 2009. Vol. 9, no. f. Pp. 77-84].

23. C. F. Shih, Elastic-plastic analysis of combined mode crack problems: Ph. D. Thesis. Cambridge, M.A.: Harvard University, 1973.

24. C. F. Shih, “Small scale yielding analysis of mixed mode plane-strain crack problems” / In: Fracture Analysis'. Proceedings of the National Symposium on Fracture Mechanics. Part 2 (College Park, Md., August 27-29, 1973). Philadelphia, Pa.: American Society for Testing and Materials, 1974. 187-210 pp.

25. В. И. Астафьев, A. H. Крутов, “Распределение напряжений вблизи вершины наклонной трещины в нелинейной механике разрушения” // Изв. РАН. МТТ, 2001. № 5. С. 125-133; англ. пер.: V. I. Astaf’ev, А. N. Krutov, “Stress distribution near the tip of an inclined crack in the nonlinear fracture mechanics” // Mech. Solids, 2001. Vol. 36, no. 5. Pp. I0I-I08.

26. Астафьев В. И., Радаев Ю.Н., Степанова Л. В., Нелинейная механика разрушения. Самара: Самарский университет, 2001. 632 с. [V. I. Astaf’ev, Yu. N. Radayev, L. V. Stepanova, Nonlinear Fracture Mechanics. Samara: Samara State Univ., 2001. 632 pp.]

27. V. N. Shlyannikov, Elastic-Plastic Mixed-Mode Fracture Criteria and Parameters / Lecture Notes in Applied and Computational Mechanics. Vol. 7. Berlin: Springer, 2003. 246 pp.

28. V. Kumar; M. D. German, C. F. Shih, An engineering approach for elastic-plastic fracture analysis: General Electric Company, NP-I93I, Research Project I237-I, Topical Report. Schenectady, New York, 1981. 71 pp.

29. E. М. Адылина, С. А. Игонин, Л. В. Степанова, “О нелинейной задаче на собственные значения, следующей из анализа напряжений у вершины усталостной трещины” // Вестн. СамГУ, 2012. №3/1(94). С. 83-102. [Е. М. Adylina, S. A. Igonin, L. V. Stepanova, “About a non-linear task on eigenvalues incurring from the analysis of tensions at the fatigue crack tip” // Vestnik SamGJJ, 2012. no. 3/1(94). Pp. 83-102].

Поступила в редакцию 10/XI/2012; в окончательном варианте — 15/XII/2012.

MSC: 74Е35, 74К20

ASYMPTOTIC METHODS OF NONLINEAR FRACTURE MECHANICS: RESULTS, CONTEMPORARY STATE AND PERSPECTIVES

L. V. Stepanova, E. M. Adylina

Samara State University,

1, Academician Pavlov St., Samara, 443011, Russia.

E-mails: lst@ssu. samara.ru, kateadulina@mail .ru

In the paper the brief review of the important results of nonlinear fracture mechanics recently obtained, by the asymptotic methods and, perturbation techniques is given.

The asymptotic solutions of the whole class of crack problems are discussed. Applications of the asymptotic analysis for crack problems of nonlinear fracture mechanics such that crack tip fields in damaged medium, fatigue crack growth problems in damaged medium, nonlinear eigenvalue problems arising from fracture mechanics analysis, nonlinear eigenvalue analysis including higher order terms are elucidated,.

Key words: asymptotic methods, nonlinear fracture mechanics, crack tip stress field, nonlinear eigenvalue problems.

Original article submitted 10/XI/2012; revision submitted 15/XII/2012.

Larisa V. Stepanova (Dr. Sci. (Phys. & Math.)), Professor, Dept, of Mathematical Modelling in Mechanics. Ekaterina M. Adylina, Postgraduate Student, Dept, of Mathematical Modelling in Mechanics.

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