CC BY
79
УДК 539.42
АВТОМОДЕЛЬНОЕ РЕШЕНИЕ ЗАДАЧИ
0 СМЕШАННОМ ДЕФОРМИРОВАНИИ ПЛАСТИНЫ С ТРЕЩИНОЙ В СРЕДЕ С ПОВРЕЖДЕННОСТЬЮ1
© 2013 Л.В. Степанова, Е.М. Адылина2
В статье получено асимптотическое решение задачи определения напряженно-деформированного состояния и поля сплошности в окрестности вершины трещины в образце, находящемся в условиях смешанного деформирования. На основании автомодельного представления решения и гипотезы о формировании области полностью дефрагментированного материала вблизи вершины трещины получено распределение напряжений, скоростей деформаций и сплошности у стационарной трещины в среде с поврежденностью в полном диапазоне смешанных форм деформирования (от чистого сдвига до нормального отрыва). Построены высшие приближения в асимптотических разложениях полей напряжений, скоростей деформаций ползучести и сплошности.
Ключевые слова: напряженно-деформированное состояние у вершины трещины, степенной определяющий закон, смешанное деформирование, метод разложения по собственным функциям, нелинейная задача на собственные значения, поврежденность, параметр поврежденности.
1. Автомодельные представления решений.
Приложение к смешанному деформированию образцов с трещинами
В нелинейной механике разрушения в последнее время сложилось представление о процессе разрушения как о многомасштабном, многоуровневом процессе [1-3]. Аналитическое моделирование многомасштабного, многоуровневого характера разрушения и изменения микроструктуры материала сопряжено с трудностями, обусловленными невозможностью континуальной механики сплошных сред различать явление на разных масштабах [2]. В [4] отмечается, что нужно понять не только явление на каждом из многих диапазонов масштабов, но и взаимодействие явлений в очень разных масштабах. Г.И. Баренблатт пишет [4], что подобные взаимодействия постоянно проявляются в таких областях, как материаловедение и биология. Здесь же можно продолжить и естественно назвать механику
1 Работа поддержана Российским фондом фундаментальных исследований (проект 12-08-00390).
2Степанова Лариса Валентиновна ([email protected]), Адылина Екатерина Михайловна ([email protected]), кафедра математического моделирования в механике Самарского государственного университета, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
разрушения, где по самому определению влияние параметров, описывающих явление в малых масштабах, распространяется по всему диапазону масштабов рассматриваемой задачи: микроструктура материала может оказывать существенное влияние на макроскопическое поведение образца. Поэтому целесообразно создать и использовать такие модели, которые бы связывали микроструктуру материала и ее влияние на макроскопическое поведение образца. Неполная автомодельность и автомодельные решения второго рода являются, по всей видимости, средством понимания того, как это влияние распространяется, и инструментом решения таких задач с помощью компьютера [4].
В механике хрупкого разрушения подобная модель может быть создана на основе использования сингулярного поведения напряжений у вершины трещины или углового выреза с различными показателями степени в асимптотическом разложении напряжений вблизи кончика трещины. В линейной механике разрушения (ЛМР) в качестве базового распределения может быть выбрано классическое решение Уильямса с корневой особенностью напряжений в непосредственной окрестности вершины трещины. Более сильная или более слабая особенности могут быть введены для моделирования степени поврежденности (степени разрушения) материала у кончика трещины на уровне мезо- и микроструктуры материала [2]. С математической точки зрения данная модель приводит к необходимости определения всего спектра собственных значений в задачах на собственные значения, следующих из проблемы нахождения напряженно-деформированного состояния у вершины трещины [5-7]. Если в ЛМР весь спектр собственных значений найден и построено полное асимптотическое представление полей напряжений и перемещений (решение М. Уильямса) [8], то в нелинейной механике разрушения остается открытым вопрос о спектре собственных значений нелинейной задачи на собственные значения, следующей из проблемы отыскания полей напряжений у вершины трещины в материалах со степенными определяющими уравнениями. Особенно актуальной проблема определения всего спектра собственных значений стала при изучении смешанного деформирования образцов с дефектами и трещинами [9-11], поскольку очевидно, что в линейной механике справедлив принцип суперпозиции, следовательно, решения, описывающие деформацию отрыва, складываются с решениями, отвечающими трещине поперечного сдвига. Для материалов с нелинейными определяющими соотношениями, очевидно, принцип суперпозиции не может быть применен, и анализ напряженно-деформированного состояния в случае смешанного нагружения элемента конструкции представляет собой самостоятельную задачу.
Вопросам анализа напряженно-деформированного состояния у вершины трещины в различных материалах в условиях пластического течения, ползучести посвящается большое количество работ [12-19]. Особенно интересными в последнее время стали вопросы смешанного нагружения элементов конструкций с дефектами [9-11]. Если для антиплоского сдвига, нормального отрыва и поперечного сдвига получен класс автомодельных решений [14], то для комбинации нескольких форм деформирования результаты полностью отсутствуют. В [16] определен весь спектр собственных значений в нелинейной задаче на собственные значения, следующей из проблемы нахождения напряженно-деформированного состояния (НДС) вблизи кончика трещины в материале со степенным определяющим законом в условиях смешанного нагружения для всех значений параметра смешанности на-гружения, характеризующего его вид. Параметр смешанности изменяется от 0 (что соответствует поперечному сдвигу) до 1 (что отвечает нормальному отрыву).
В [16] решение нелинейных задач на собственные значения получено численно, что позволило найти весь спектр собственных значений для степенного определяющего закона. Проведенное качественное исследование собственных значений в [16] позволило в настоящем исследовании использовать свойство автомодель-ности процесса накопления повреждений в твердом теле (для степенного закона накопления повреждений) и ввести новую автомодельную переменную.
В настоящей работе на основе автомодельного представления решения задачи определения полей напряжений, деформаций и сплошности у вершины трещины в твердом теле в условиях ползучести и введения автомодельной промежуточной асимптотики получен скейлинговый закон накопления повреждений на мезоуровне (автомодельный закон роста области активного накопления повреждений и области полностью поврежденного материала). Получены автомодельные решения задач отыскания НДС в образцах с трещинами для смешанных форм деформирования элементов конструкций с трещинами в среде с поврежденностью (во всем диапазоне смешанных форм деформирования). Вид смешанного нагру-жения задается параметром смешанности нагружения Ши [20-22] (и определяется как соотношение окружной и сдвиговой компонент тензора напряжений на линии продолжения трещины)
2
Мр = - аг^ п
Иш ^(г'в = 0)
>0 отв (г, в = 0)
(1.1)
Данная работа развивает подход, примененный в [14; 15], и использует автомодельное представление решения для анализа полей напряжений, скоростей деформаций и сплошности у вершины трещины в случае смешанного деформирования.
2. Автомодельное решение задачи о трещине в среде с поврежденностью в условиях смешанного деформирования (нормальный отрыв и поперечный сдвиг). Промежуточная автомодельная асимптотика
2.1. Анализ размерностей. Автомодельное представление решения
Целью настоящего анализа является определение полей напряжений, деформаций и сплошности вблизи вершины трещины, находящейся в условиях смешанного деформирования (рис. 1).
Для рассматриваемой задачи уравнения равновесия в полярной системе координат имеют вид
до-тт +1 дагв + &гт - оее =о 1 доев + дотв + 2огв =о (21)
дг г дв г ' г дв дг г
В плоской задачи разрешающим относительно искомой функции напряжений Эри выступает одно условие совместности деформаций
2 = гд£тт + гд2 (г£ее) (22)
2 дг Г дв = дв2 дг +г дг2 ■ ( )
Поведение материала описывается с помощью определяющих уравнений, построенных на основе степенного закона Бейли-Нортона, теории установившейся ползучести с применением концепции эффективного напряжения:
£гз = (3/2)Б(ае/ф)п-1в3/ф, (2.3)
где ае = \j2sij- интенсивность касательных напряжений для случая плоского деформированного состояния, определяемая равенством ае = = \f3yi(агг — адд)2 + 4о1в/2, ф - параметр сплошности, эволюционирующий в соответствии со степенным законом накопления повреждений
= —А (Она)т, (2.4)
<и\ф) у '
аеду = аае + ва1 + (1 — а — в)акк - эквивалентное напряжение.
Процесс активного накопления рассеянных повреждений происходит в непосредственной окрестности вершины трещины, где образуется область диспергированного (дефрагментированного, разрушенного) материала, а на удалении от нее параметр сплошности стремится к единице, что отвечает неповрежденному материалу. Это позволяет сформулировать асимптотическое граничное условие в бесконечно удаленной точке как условие асимптотического сближения с решением ХРР - решением для степенных определяющих уравнений (2.3) [23-25], где ф = 1 (таким образом, рассматривается непосредственная окрестность вершины трещины). В соответствии с гипотезой о маломасштабной поврежденности на больших расстояниях от вершины трещины (больших по сравнению с характерным линейным размером области полностью поврежденного материала, но все еще малых по сравнению с длиной трещины, характерным линейным размером образца) поле напряжений определяется решением ХРР (решением аналогичной задачи без учета процесса накопления повреждений ф = 1):
( С* ч1/^1)
аз (г ^ж,е,г)=1 Бцпг) (в>п)> (2.5)
где С * - инвариантный интеграл теории установившейся ползучести. Начальное
условие при г = 0 и граничное условие в бесконечно удаленной точке (2.5) совпадают, поскольку они задаются решением задачи для ф = 1.
Анализ размерности величин, входящих в уравнения (2.3) - (2.5), позволяет установить, что для определяющих соотношений (2.3), кинетического уравнения (2.4), начальных и граничных условий (2.5) существует автомодельная переменная
Я = г (Аг)-(п+1)/т Шп/С*, (2.6)
и система уравнений задачи допускает автомодельное представление решения: а^(г, в, г) = (А£)-1/т а^(Я, в), ф(г,в,г) = ф(Я,в).
Следует отметить, что граничное условие в бесконечно удаленной точке может быть сформулировано в более общей по сравнению с (2.5) форме
а^ (г ^ то, в, г) ^ Сг3а^(в,п), (2.7)
где показатель степени в подлежит определению в ходе решения задачи, С - амплитуда поля напряжений на бесконечности, определяемая геометрией реального образца и системой приложенных нагрузок. Для степенных определяющих соотношений (2.3), кинетического уравнения (2.4) и более общих граничных условий (2.7) существует автомодельная переменная
( _ \ 1/(вт)
Я = г (жст) (2.8)
(автомодельные решения, основанные на введении автомодельных переменных, степени которых должны быть определены в ходе решения задачи, относят к автомодельным решениям второго рода [4]).
После введения автомодельной переменной уравнения равновесия (2.1), определяющие уравнения (2.3), условие совместности (2.2), сохраняют свою форму, а кинетическое уравнение накопления повреждений принимает вид
( лт
Яф1П = -вт [ае/ф) . (2.9)
Далее знак обозначающий безразмерную величину, опускается.
2.2. Асимптотическое решение задачи. Метод разложения по собственным функциям
Асимптотическое решение задачи вне области диспергированного материала, формирующейся у вершины трещины (на больших расстояниях от вершины трещины Я ^ то) разыскивается в форме
^ то
х(Я,в) = £ Ях«+1Л(в), Ф(Я, в) = 1 Я^дк(в), (2.10)
к=0 к=0
где Хк < 0,7к < 0,/к(в),дк(в) подлежат определению в ходе решения задачи. В силу (2.10) асимптотическое представление компонент тензора напряжений вне области диспергированного материала имеет вид
то то
авв(Я, в) = £ Ях*-1(\к + 1)Лк/к(в), ат(Я, в) = - £ Ях"-1Лк/'к(в),
к=0 то к=0 (2.11)
аня(Е, в) = £ Ях»-1 [(Лк + 1)/к(в) + /''(в)].
к=0
Асимптотический анализ уравнений задачи позволяет установить, что справедливы равенства Ак — Ао = к(А1 — Ао) для всех к, и асимптотическое представление интенсивности касательных напряжений принимает форму
где
о е (Е,в) = ЕЛо — 1о(0\в)
1 + ]Т КХк
к=1
Л о(к)(в)
о(0) = л/[(1 — А0)/о + /0']2 + 4А§(/0)2,
о(1 = /е(1)/(о(0))2,
/(1) = [(1 — А2)/о + /0'][(1 — А2)/1 + /1'] + 4АоА1/0/1,
о(2)
/(2) = [(1—А2)/1 + /1']2+2[(1—А2)/о+/0'][(1—А2)/2+/2']+4а2/2+8А0А2/1 /2,
/(3) =2[(1 — Х2)/о+/0'][(1 —А3)/з + /З']+2[(1 — А2)/1 +/1'][(1 — А2)/2 + /2']+
+8А0А3/0/3 + 8 А1Х2 У1 ,
(2)
2
(о(0))2 (о
Г(3) = 1
(з)
(о(0)
//
)2 (о(0))4
+
(/М
\ (о (0) )2 у
о(4)
-5
(ЛЦ
\ (о (0) )2 у
+ 6
/ (1) 2 / (2)
(0))б
(о (0))
^ /(2) А
\ (о (0))2 у
2 ^(1)^(3) .(4) — 4/е /е__+4 /е
(о (0) )4
(о (0) )2
/е(4) = [(1 — А2)/2 + /2']2 + 2[(1 — А2)/о + /0'][(1 — а2)/4 + /4']+
у2\ ^ _1_ _1_ /1 \2^/\2
/2 ] + 2[(1 — А0)
+2[(1 — Х2)/1 + /1'][(1 — А3)/з + /3'] + 4А2(/2)2 + 8А0А4/' /4 + 8А1А3/1/3.
Кинетическое уравнение накопления повреждений и гипотеза о том, что порядки малости слагаемых, последовательно появляющихся в этом уравнении в правой и левой частях, позволяют связать показатели степеней в асимптотических разложениях функции напряжений Эри и параметра сплошности (2.10): 7 = (Ао — 1)т, 1к = (Ао — 1)т+к(А1 — Ао). Анализ асимптотических разложений компонент тензора скоростей деформаций ползучести дает возможность установить, что А1 — Ао = = 7 = (Ао — 1)т, откуда легко найти, что 7к = к(Ао — 1)т, и получить следующие асимптотические представления компонент тензора скоростей деформаций ползучести вне области полностью поврежденного (диспергированного) материала:
¿Ея(Е,в) = д(Ло-1)п40^ (в) + д(Ло-1)(п+т)4^ (в) + д(Ло-1)(п+2т)42^ (в)+
(п
+Д(Ло-1)(п+3т)£030(в) + Е(Ло-1)(п+4т)£(040 (в) + . ..
¿Ев(Е, в) = Е(Ло-1)п£(001(О) + Е(Ло-1)(п+т)£(11(в) + Е(Ло-1)(п+2т)£(02)в(О)+
'Ев
+Е(Ло-1)(п+3т)£(03в)(в) + Е(Ло-1)(п+4т)£(04в)(в) + ...
где угловые распределения асимптотических разложений компонент тензора деформаций ползучести определяются следующими формулами:
угловые распределения главного члена асимптотического разложения компонент тензора скоростей деформаций
¿(0) " ЕЕ(
(в) = (о(0))п-1 [(1 — А2)/о + /0'], £(Ев)(в) = —2 (о(0))"~1 Ао/0;
п— 1
коэффициенты первого слагаемого в асимптотическом разложении компонент тензора скоростей деформаций
£(1) ( ¿ЕЕ(
(в)= (о(0))п ^[(1 — Х2)/1 + /1'] + [(1 — А2)/о + /0'] [(п — 1)о(1) + п^1]}
/ п— 1 { Г //Л \
¿Е)(в) = —2 (о(0))" {А1/1 + Ао/0 [(п — 1)о((1) + пд1] } ;
3
е
е
(0^4
К
2
е
коэффициенты второго слагаемого асимптотического разложения компонент тензора скоростей деформаций
£ЯЯ(
+
г(в) = (аео})п-1 {[(1 - Л2)/2 + /2'] + [(1 - Л/ + . + [(1-Л2)/0 + /0'] {^^[(п-2)(а(1))2 +2а(2)]+п(п-1)д1а(1) + 2[(п +1)д2 + 2д2]}}.
е™(в) = -2 и0)П-иЛ2/2 + Л1/1 [(п - 1)а«1} + пдА +
/1 ] (п - 1)а(1) + пд1 п
+Л0/0 { ^(п - 2)(а(!))2 + 2а(2)] + п(п - 1)^ + 2[(п + 1)д2 + 2д2]} } ;
коэффициенты третьего слагаемого асимптотического разложения компонент тензора скоростей деформаций
(в)
аео)
1
[(1 - Л2)/з + /''] + [(1 - Л2)/2 + /2'] (п - 1)ае1} + пд1
+
+ [(1-Л2)/1 + /1'] { ^[(п-2)(а(1))2 + 2а(2)]+п(п-1)^) + |[(п + 1)д? + 2д2]} +
+ [(1-Л0)/0 + /0']
п-1
(п-2)(п-3)(а(1))3 + 6(п-2)а(1)а(2) + 6а(3)
+
+ 3! [(п + 1)(п + 2)д3 + 6(п + 1)дд + 6д^ + пд1[(п-2)(а(1))2 +2а£2)]+ п г 2 Г(п"
+ п [(п + 1)д2 +2д2^} ,
4>) = -2
а(0)
п-1
+ Л2/2 [(п - 1)а£Х) + пд1 (1)
+
+Л1 /1 { ^^[(п-2)(а(!))2 +2а£2)]+п(п-1)д1 а« + 2[(п + 1)д? + 2д2]} +
+Л0/Л п [(п + 1)(п + 2)д3 + 6(п + 1)д1д2 + 6дз] + пд1[(п-2)(а(1))2 +2а£2)] +
+
п1
3!
(п-2)(п-3)(а(1))3+6(п-2)а(1)а(2) +6а£3)] [(п + 1)д? + 2д2]
коэффициенты четвертого слагаемого асимптотического разложения компонент тензора скоростей деформаций
£}»(в) = (а(0)) П-1 ,
р(4)(в)= 2 {а(0)\П 1 р(4)
£Пв(в) = -2 {^е ) £Пв,
= (1 - Л2)/4 + /4' + [(1 - Л3)/з + /31 [(п - 1)а(1) + пд1
п1
Л2)/2 ■ г"' ■
п1
+
+ [(1-Л2)/2 + /2'] {п-^[(п-2)(а(1))2 +2а£2)]+п(п-1)^ + 2[(п+1)д2 +2д2]} +
+ [(1-Л? )/1 + /1']
3!
[(п - 2)(п - 3)(а(1))3 + 6(п - 2)а(1)а(2) + 6а£3)] +
п1
+ ^[(п - 2)(а(1))2 + 2а(2)] + п(п - ^[(п + 1)д2 + 2д2]+
+ 3[(п + 1)(п + 2)д3 + 6(п + 1)д1д2 + 6д3]} +
+ [(1-Л2)/0 + /0']{ п-1[(п-2)(п-3)(п-4)(а(1))4 + 12(п-2)(п-3)(а(1))2а£2)] +
+ О-1 [24(п - 2)(а£2))2 + 24(п - 2)а«а£3) + 24а№] +
+
4! 1
п — 1
-пд1
(п - 2)(п - 3)(а(1))3 + 6(п - 2)ае1}а(2) + 6а£3)
+
+ П^(п + ^ + ^[(п - 2)(а(1))2 + 2а(2)]+
+ -(п - 1)а(х)[(п + 1)(п + 2)д? + 6(п + 1)д^ + 6дз] +
+ П [(п + 1)(п + 2)[(п + 3)д4 + 12^2^2] + 24(п + 1)[д2 + дд] + 24^] }
^яв — л4/4 + Аз/з
(п - 1)а(1) + пд1
+
+А2/2 {п-1[(п-2)(а(1))2 +2а (2)]+п(п-1)д1а« + |[(п+1)д? +2д2]} +
п1
+А1/^ "3т[(п - 2)(п - 3)(^(1))3 + 6(п - 2)ст(1)ст(2) + 6а(3)] +
+ п^[(п - 2)(а(1))2 + 2а (2)] + 2(п - 1)а«[(п + 1)д2 + 2д2] + „з
+ 3[(п + 1)(п + 2)д3 + 6(п + 1)дд + 6дз]} +
+Ас/^ ~-Г~ [(п-2)(п-3)(п-4)(а(1))4 + 12(п-2)(п-3)(а(1))2а(2)] + + [24(п - 2)(а(2))2 + 24(п - 2)а(1Ч3) + 24а(4)]+
п1
(п - 2)(п - 3)(а(1))3 + 6(п - 2)а(1)а(2) + 6а(3)
+
+ -пр пд1
+ пп^[(п + 1)д2 + 2д2][(п - 2)(а(1))2 + 2а(2)]+ ^
+ -(п - 1)а(1)[(п + 1)(п + 2)д3 + 6(п + 1)д^ + 6дз]+ '
+ п [(п + 1)(п + 2)[(п + 3)д4 + 12д2д2] + 24(п + 1)[д2 + дд] + 24д4^ .
Из условия совместности деформаций можно получить систему обыкновенных дифференциальных уравнений для определения угловых распределений компонент тензора напряжений. Для нахождения функции /о (в) получается нелинейное ОДУ
лД0) ,2 (0)
2[(Ао - 1)п + 1]-^ — -¡Ш - (Ао - 1)п[(Ао - 1)п + 2]е%. (2.13)
Относительно функций /к (в) (к ^ 1) из условия совместности (2.2) следует система линейных обыкновенных дифференциальных уравнений
¿е(к) ¿2£{к)
, -пс, и & р т?
2[(Ао - 1)(п + кт) + 1]-
'Ев — " °дд
ив ив2
-(Ао - 1)(п + кт)[(Ао - 1)(п + кт) +2]еЕкЕЕ. (2.14)
Решение системы уравнений (2.13), (2.14) должно удовлетворять граничным условиям /к(в — ±п) — 0 и /к(в — ±п) — 0. Численное решение системы уравнений (2.13), (2.14) на отрезке [-п, п] разыскивалось с помощью описываемой ниже процедуры в пакете МаШетаИса V. 5.1.
В случае смешанного деформирования соображения симметрии и антисимметрии использованы быть не могут, и необходимо искать решение уравнения (2.13) на отрезке [-п,п]. В условиях смешанного нагружения при численном решении уравнения (2.13) отрезок интегрирования [-п, п] можно разбить на два отрезка: [-п, 0] и [0,п]. Сначала уравнение (2.13) интегрируется на отрезке [0, п], и двухточечная краевая задача сводится к задаче Коши с начальными условиями
/0(в = 0) = 1, /0(в = 0) = -(Л + 1)/tg(Mpп/2), /0'(в = 0) = А2, /0''(в = 0) = А3.
(2.15)
Значение производной первого порядка находится из заданного вида смешанного нагружения (значение параметра смешанности нагружения (1.1), задающего вид нагружения, известно). Неизвестные постоянные А2 и А3 определяются таким образом, чтобы выполнялись краевые условия на верхнем берегу трещины:
После того как подобраны постоянные А2 и А3, уравнение (2.13) интегрируется на отрезке [-п, 0], для чего двухточечная краевая задача для уравнения (2.13) с граничными условиями
/0(в = -п) = 0, /0(в = -п) = 0, /0'(-п) = В2, /0''(в = -п) = В3. (2.18)
Неизвестные постоянные В2 и В3 подбираются таким образом, чтобы выполнялись условия равновесия элемента, расположенного на луче в = 0. Уравнения равновесия данного элемента требуют непрерывности компонент тензора напряжений адд и агд на луче в = 0, что влечет за собой непрерывность функций /0(в) и /0(в) при в = 0 (и, следовательно, краевые условия (2.16)). Поэтому две неизвестные постоянные В2 и В3 определяются таким образом, чтобы решение, разыскиваемое на отрезке [-п, 0], удовлетворяло бы краевым условиям при в = 0.
2.3. Конфигурация области диспергированного (дефрагментированного) материала
После нахождения численного решения системы уравнений (2.13) и (2.14) можно определить границу области диспергированного материла в окрестности вершины трещины с помощью уравнений
ф(Я,в) = 1 - Я7д0(в) = 0
ф(Я, в) = 1 - Я^д0 (в) - Я71 д1(в) = 0,
ф(Я, в) = 1 - Яд0(в) - К'1 д1(в) - К'2д2(в) = 0,
ф(Я,в) = 1 - Я^ д0 (в) - Я^ д1(в) - Я^ д2(в) - Я^ д3(в) = 0,
ф(Я,в) = 1 - К'д0(в) - К'1 д1(в) - К'2д2(в) - К'3д3(в) - К'4д4(в) = 0,
/0(в = п) = 0, /0 (в = п) = 0.
(2.16)
(2.17)
где
Остальные коэффициенты асимптотического разложения параметра сплошности находятся подобным образом. Решение уравнения (2.13) вместе с условиями отсутствия поверхностных усилий приводит к ранее рассмотренной нелинейной задаче на собственные значения. Оказалось, что если в качестве условия в бесконечно удаленной точке принять условие асимптотического сближения с решением ХРР (Л = п/(п + 1)) и предположить, что вне области активного накопления повреждений решение выходит на асимптотику ХРР, то границы области диспергированного материала, построенные с помощью двучленного и трехчленного асимптотических разложений параметра сплошности, значительно отличаются друг от друга по своей форме и размерам, что свидетельствует о том, что разложения (2.10) не имеют асимптотической природы, ибо каждое следующее слагаемое в асимптотическом разложении должно быть малой поправкой к предыдущему. Однако, если предположить, что в окрестности вершины трещины в условиях смешанного деформирования имеется область с промежуточной асимптотикой поля напряжений (где собственные значения отличаются от значений задачи ХРР и определяются значениями, приведенными в [16]), то асимптотические разложения (2.10) приводят к границам области диспергированного материала, сходящимся к предельному контуру. Вычисленные собственные значения для всего диапазона значений параметра смешанности нагружения МР для п = 5 приведены в табл. 1. Следует отметить, что альтернативный подход решения нелинейных задач на собственные значения, базирующийся на методе малого параметра, изложен в [16; 26].
Таблица 1
Собственные значения нелинейной задачи на собственные значения для различных значений параметра смешанности нагружения (п = 5)
М Р Л /' '(0) /'' '(0) /' '(-п) /'' '(-п)
МР = 0,05 -0, 237872 3,080269 7,474145 17,168290 -16, 59101
МР = 0,1 -0, 221158 2,007918 3,081809 9, 286569 -8, 84890
МР = 0, 2 -0, 202342 1,142387 0, 988910 5, 316140 -4, 98511
МР = 0, 3 -0,194812 0, 753230 0, 302029 3, 972174 -3, 70051
МР = 0,4 -0,192320 0, 546480 -0,068735 3, 263461 -3, 03392
МР = 0, 5 -0,194290 0,472467 -0,364145 2, 830609 -2, 63612
МР = 0, 6 -0, 200055 0,463181 -0,626727 2, 539084 -2, 37534
МР = 0, 7 -0, 204759 0,430847 -0,818362 2, 269456 -2,13255
МР = 0, 8 -0, 210625 0,405497 -0,988799 2,029076 -1, 91625
МР = 0, 9 -0, 218248 0, 394459 -1,159518 1, 811569 -1, 72195
МР = 0, 95 -0, 221500 0, 394789 -1, 248002 1, 707946 -1, 62792
Поэтому асимптотическое решение (2.10) можно рассматривать как автомодельное предельное решение задачи о трещине в среде с поврежденностью в условиях смешанного нагружения для всех форм смешанного деформирования от чистого сдвига до нормального отрыва. Конфигурации областей полностью поврежденного материала для новой асимптотики показаны на рис. 2-4, где буквой к обозначена граница области диспергированного материала, построенная с помощью к + 1-членного асимптотического разложения параметра сплошности. На рис. 5-8 изображены линии уровня параметра сплошности.
Рис. 2. Геометрия области полностью поврежденного материала для Мр — 0; 0,05, 0,1; 0, 2; 0, 3 и Мр — 0.4
2.4. Амплитудный масштабный множитель С
После возвращения к исходным переменным поле напряжений вне области полностью поврежденного материала определяется с помощью равенства
а^(г, в,г) — (Агт)-1/тс [дЛ»-1а(о)(в)+дЛ1-1а|^1)(в)+дЛ2-142)(в) + ...], (2.19)
где Д — г(АгСт)1/((Л-1)т). Отличительной особенностью полученного представления является зависимость поля напряжений только от одного амплитудного
Рис. 3. Геометрия области полностью поврежденного материала для МР = = 0, 5; 0, 6; 0, 7; 0, 8;0, 9 и МР = 1
множителя С, отражающего влияние поля напряжений в окрестности вершины трещины от геометрии образца и системы приложенных нагрузок. Поэтому для нахождения поля напряжений вблизи кончика трещины достаточно определить всего один масштабный множитель (амплитудный коэффициент) С. В традиционной механике трещин используется подход, приводящий к задаче нахождения трех или более коэффициентов асимптотических разложений поля напряжений. В последнее время в нелинейной механике разрушения, одной из задач которой является определение напряженно-деформированного состояния в непосредственной
Рис. 4. Геометрия области полностью поврежденного материала и линии уровня параметра сплошности для п — 5, т — 4, Мр — 0,05; 0,1
Рис. 5. Геометрия области полностью поврежденного материала и линии уровня параметра сплошности для п — 5, т — 4, Мр — 0, 2; 0, 3
Рис. 6. Геометрия области полностью поврежденного материала и линии уровня параметра сплошности для п — 5, т — 4, Мр — 0,4; 0, 5
Рис. 8. Геометрия области полностью поврежденного материала и линии уровня параметра сплошности для п = 5, т = 4, Мр = 0, 8; 0, 9
окрестности вершины трещины, большое внимание уделяется построению высших приближений в асимптотических разложениях механических величин (полей напряжений, деформаций и перемещений) в окрестности вершины трещины в нелинейных материалах.
Для материалов, следующих степенному закону установившейся ползучести ¿щ = ЗВа^-1 вц/2, асимптотические разложения компонент тензора напряжений могут быть представлены в форме
ац(г, в) = (С*/(В1пг))1/(п+1) аЦ\в)+Азг^а^ (в) + 1А3г^а^ (в) +.... (2.20)
Первое слагаемое в (2.20) в литературе носит название решения ХРР. В (2.20) фигурирует С*-интеграл, определяемый равенством
С* = J (]¥3х2 - ацпц, (2.21)
где IV = /д ы ацсИц = —аР^е. Асимптотический анализ асимптотического выражения интенсивности напряжений позволяет установить, что показатели степеней
в (2.20) связаны между собой вз - в! — 2(в2 - в!). Детальный асимптотический анализ трехчленного асимптотического разложения компонент тензора напряжений (2.20) и интенсивности касательных напряжений показывает, что константы А2 и А3 не являются независимыми, они связаны равенством А3/А1 — (А2/А1) . В литературе имеется достаточно большое число оценок инвариантного С*-интеграла, полученных для различных образцов, наиболее часто используемых в эксперименте. Однако остаются открытыми многие вопросы, связанные с определением второго коэффициента в асимптотических разложениях, а именно А2. Естественно предположить, что А2 в общем случае является функций приложенных нагрузок, геометрических характеристик образца и показателя нелинейности материала п. Определение этого коэффициента и составляет одну из важных задач нелинейной механики разрушения, поскольку даже известные оценки, приведенные в справочной литературе, дают существенную погрешность. Например, для образца с односторонним боковым надрезом для инвариантного С*-интеграла имеется оценка С * — БаоеоШ (1 - о)о^ (Р/Ро), где Ро — 1,4555^^ (1 - а)ао, Н1 - безразмерная функция приложенной нагрузки и геометрии образца. Оказывается, что уточненные вычисления, проведенные в последнее время, показали, что значения функции могут отличаться от приведенных в справочниках более чем на сто процентов. Следовательно, с появлением вычислительных комплексов, таких как ЛЪадив, предоставляется возможность получения более точных оценок для параметров механики разрушения, что и составляет одну из задач современной механики деформируемого твердого тела. Амплитудные коэффициенты, входящие в структуру полей напряженно-деформированного состояния в области вершины трещины, как известно, являются характеристиками сопротивления материала разрушению. Эти коэффициенты могут описывать текущее состояние или принимать предельные значения в зависимости от уровня накопленных повреждений в исследуемом объекте. Для нахождения оставшейся постоянной А2 применяются различные подходы, включающие конечно-элементный расчет для образцов различной конфигурации, находящихся в условиях разных систем усилий. Тем не менее в общем случае исследование данной проблемы не является математически завершенным и требует объединения полученных разными исследователями результатов в единый справочник, подобный существующим справочникам, содержащим оценки коэффициентов интенсивности напряжений и ^интегралу.
Альтернативой описанному выше подходу является учет процессов накопления повреждений в окрестности вершины трещины, что приводит к асимптотическому представлению механических полей вблизи кончика трещины (2.19), содержащему лишь один амплитудный множитель, что значительно упрощает решение задачи. Один амплитудный фактор С может быть определен с помощью инкорпорирования скалярного параметра сплошности (поврежденности) в расчетные комплексы, осуществляющие конечно-элементный анализ.
Результаты и выводы
В работе получено новое промежуточно-асимптотическое решение задачи о стационарной трещине в среде с поврежденностью в образце, находящемся в условиях смешанного деформирования (во всем диапазоне смешанных форм деформирования от чистого сдвига до нормального отрыва). В ходе численного анализа найдена новая асимптотика дальнего поля напряжений. Асимптотику поля напряжений можно интерпретировать как промежуточную асимптотику поля напряжений
и параметра сплошности, справедливую на расстояниях, много больших характерного линейного размера области поврежденного материала, но все еще много меньших по сравнению с длиной трещины, с характерным линейным размером тела. С использованием новой асимптотики дальнего поля напряжений у вершины трещины в условиях смешанного нагружения построены конфигурации областей диспергированного материала. Выведены пятичленные асимптотические разложения компонент тензора напряжений и шестичленные асимптотические разложения скалярного параметра сплошности и показано, что шестичленные асимптотические разложения параметра сплошности приводят к контурам области полностью поврежденного материала, сходящимся к некоторому предельному контуру.
Литература
[1] Пестриков В.М., Морозов Е.М. Механика разрушения. СПб.: Профессия, 2012.
[2] Си Дж. Мезомеханика, понятие сегментации и мультискейлинговый подход: нано-микро-макро // Физическая мезомеханика. 2008. Т. 11. № 3. C. 5-18.
[3] Бьюи Х.Д. Механика разрушения: Обратные задачи и решения. М.: Физмат-лит, 2011. 412 с.
[4] Баренблатт Г.И. Автомодельные явления - анализ размерностей и скейлинг. Долгопрудный: Интеллект, 2009.
[5] Stepanova L. Eigenspectra and orders of stress singularity at a mode I crack tip for a power-law medium // Comptes Rendus Mechanique. 2008. V. 336. № 1-2. P. 232-237.
[6] Степанова Л.В. О собственных значениях в задаче о трещине антиплоского сдвига в материале со степенными определяющими уравнениями// Прикладная механика и техническая физика. 2008. Т. 49. № 1. С. 173-180.
[7] Степанова Л.В., Федина М.Е. Автомодельное решение задачи о трещине антиплоского сдвига в связанной постановке (ползучесть-поврежденность) // Прикладная механика и техническая физика. 2002. Т. 43. № 5. С. 114-123.
[8] Hello G., Tahar M.B., Roelandt J.M. Analytical determination of coefficients in crack-tip stress expansions for a finite crack in an infinite plane medium // Int. J. of Solids and Structures. 2012. V. 49. P. 556-566.
[9] Шлянников В.Н., Кислова С.Ю. Параметры смешанных форм деформирования для трещины в виде математического разреза // Известия Саратовского университета. Сер.: Математика. Механика. Информатика. 2009. Т. 9. Вып. 1. С. 77-84.
[10] Шлянников В.Н., Туманов А.В. Упругие параметры смешанных форм деформирования полуэллиптической трещины при двухосном нагружении // Известия Саратовского университета. Сер.: Математика. Механика. Информатика. 2010. Т. 10. Вып. 2. С. 73-80.
[11] Shlyannikov V.N. Elastic-Plastic Mixed-Mode Fracture Criteria and Parameters/ V.N. Shlyannikov. Berlin: Springer, 2003.
[12] Pan J., Lin P.C. Analytical solutions for crack-tip sectors in perfectly plastic Mises materials under mixed in-plane and out-of-plane shear loading conditions // Engng. Fracture Mechanics. 2006. V. 73. P. 1797-1813.
[13] Rahman M., Hancock J.W. Elastic perfectly-plastic asymptotic mixed mode crack tip fields in plane stress // Int. J. Solids and Structures. 2006. V. 43. P. 3692-3704.
[14] Степанова Л.В., Федина М.Е. О геометрии области полностью поврежденного материала у вершины трещины антиплоского сдвига в связанной постановке задаче (связка "ползучесть-поврежденность") // Вестник Самарского государственного университета. 2001. № 2. С. 87-113.
[15] Степанова Л.В., Федина М.Е. Автомодельное решение задачи о трещине отрыва в связанной постановке // Прикладная математика и механика. 2008. Т. 72. № 3. С. 516-527.
[16] Адылина Е.М., Степанова Л.В. О построении многомасштабных моделей неупругого разрушения // Вестник Самарского государственного университета. 2012. № 9. C. 70-83.
[17] Степанова Л.В. Уточненный расчет напряженно-деформированного состояния у вершины трещины в условиях циклического нагружения в среде с по-врежденностью // Вестник Самарского государственного унвиерситета. 2011. № 83. C. 105-115.
[18] Степанова Л.В. Анализ собственных значений в задаче о трещине в материале со степенным определяющим законом // Журнал вычислительной математики и математической физики. 2009. Т. 49. № 8. С. 1399-1415.
[19] Астафьев В.И., Шестериков С.А., Степанова Л.В. Асимптотика напряженно-деформированного состояния в окрестности вершины трещины в условиях ползучести // Вестник Самарского государственного университета. 1995. Спец. выпуск. С. 59-64.
[20] Shih C.F. Elastic-plastic analysis of combined mode crack problems // Ph. D. Thesis, Harvard University, Cambridge, M.A. 1973.
[21] Си Дж. Труды Американского общества инженеров-механиков // Конструирование и технология машиностроения. 1976. № 4. C. 113-120.
[22] Shih C.F. Small scale yielding analysis of mixed mode plane-strain crack problems // Fracture Analysis ASTM STP 560. 1974. P. 187-210.
[23] Hutchinson J.W. Singular behaviour at the end of tensile crack in a hardening material //J. Mech. Phys. Solids. 1968. V. 16. № 1. P. 13-31.
[24] Rice J.R., Rosengren G.F. Plane strain deformation near a crack tip in a power-law hardening material // J. Mech. Phys. Solids. 1968. V. 16. № 1. P. 1-12.
[25] Hutchinson J.W. Plastic stress and strain fields at a crack tip //J. Mech. Phys. Solids. 1968. V. 16. № 5. P. 337-347.
[26] Адылина Е.М., Игонин С.А., Степанова Л.В. О нелинейной задаче на собственные значения, следующей из анализа напряжений у вершины усталостной трещины // Вестник Самарского государственного университета. 2012. № 3/1(94). С. 83-102.
Поступила в редакцию 31/V/2013;
в окончательном варианте — 31/V/2013.
SELF-SIMILAR SOLUTION TO THE PROBLEM OF MIXED DEFORMATION OF PLATE WITH A CRACK TIP IN A DAMAGED MEDIUM
© 2013 L.V. Stepanova, E.M. Adylina3
In the article asymptotic solution of the problem of defining of stressed state of strain and field of continousness in the neighbourhood of crack tip in the sample which is in state of mixed deformation is received. On the basis of self-similar presentation of solution and hypothesis about the formation of area fully defragmentated material near the crack tip, the distribution of strains, speed of deformation and continuousness of a stationary crack in a damaged medium at full range of mixed forms deformation (from simple share to normal fracture) is obtained. Highest approximations in asymptotic resolution of fields of strain, speed of deformation of creep and continuous are obtained.
Key words: stressed state of strain near the crack tip, sedate determining law, mixed deformation, method of eigenfunction expansion, nonlinear problem on eigenvalues, nonlinear problem on eigen functions, damage, parameter of damage.
Paper received 31/V/2013. Paper accepted 31/V/2013.
3Stepanova Larisa Valentinovna ([email protected]), Adylina Ekaterina Mikhailovna ([email protected]), the Dept. of Mathematical Modelling in Mechanics, Samara State University, Samara, 443011, Russian Federation.
