МОДЕЛИРОВАНИЕВНАУКЕОМАТЕРИАЛАХ^
УДК 539.3
А. А. Валишин, Д. В. Мищенко
ОСОБЕННОСТИ КВАЗИХРУПКОГО РАЗРУШЕНИЯ ПОЛИМЕРОВ И КОМПОЗИТОВ НА ИХ ОСНОВЕ
Рассмотрено квазихрупкое разрушение полимеров и композитов на их основе, которое наблюдается при температуре, изменяющейся от температуры хрупкости до температуры стеклования. Определены форма и размеры зоны вынужденной эластичности, возникающей в аморфных твердых полимерах перед трещиной разрушения при температуре, превышающей температуру хрупкости. Проведены исследования движения трещины с зоной вынужденной эластичности и процессов тепловыделения при формировании этой зоны.
Е-таИ:е^@та11.ги
Ключевые слова: полимеры, композиционные материалы на их основе,
трещина, зона вынужденной эластичности.
Согласно принятой в литературе прочностной классификации, твердые полимеры подразделяют на три типа: низкопрочные, высокопрочные и свехпрочные [1]. Разрушение низкопрочных материалов, о которых пойдет речь, происходит в результате роста наиболее опасной трещины вплоть или до полного разрушения изделия (конструкции), либо до достижения неприемлемого уровня повреждения.
Трещина в материале растет под действием локальных напряжений в ее вершине. При этом неважно, каким фактором обусловлены эти локальные напряжения: внешним механическим напряжением, неоднородным температурным полем или, может быть, каким-то другим фактором. Известно, что в механическом поле трещина является концентратором напряжений: локальные напряжения вблизи вершины трещины многократно превышают напряжения вдали от нее. В температурном поле трещина является еще и концентратором теплового потока: тепловой поток и градиент температурного поля вблизи трещины имеют намного большие значения, чем вдали от нее [2].
Реальная трещина представляет собой щель с асимптотически сходящимися берегами. При этом расстояние между берегами (раскрытие трещины) плавно уменьшается вплоть до межчастичных расстояний (межатомных и межмолекулярных). Отсюда следует, что в концевой части трещины есть участок, где существенны силы межатомного и межмолекулярного взаимодействия берегов трещины
[3—5]. Межчастичные силы сцепления вблизи кончика трещины обеспечивают плавное смыкание ее берегов и конечность напряжений вблизи вершины трещины.
Кончик трещины напоминает маленький клюв. Форма трещины в области так называемого клюва, его напряженное состояние и размеры полностью определяются силами межатомного и межмолекулярного взаимодействия и не зависят от внешних нагрузок и температуры. Иными словами, клюв трещины автономен по отношению ко всей трещине, и при ее росте перемещается впереди трещины, не меняя ни своих размеров, ни формы [3, 4, 6]. Форма остальной части трещины зависит от внешних напряжений механического, термического или иного происхождения. Предельное значение размеров клюва определяется формулой
где X — среднее межмолекулярное расстояние, характеризующее наибольшее раскрытие клюва; Е — модуль Юнга; и — коэффициент Пуассона; агп — удельная поверхностная энергия разрушения. Численную оценку размера клюва проведем с использованием данных для поле-метилметакрилата (ПММА) [1]: Я = 1,2-10"9 м, Е = 3,93-109 Н/м2, и = 0,25, ап = 3,9-10"2 Дж/м2 и получим й = 7,9-10"10 м. Сравним
предельное значение размеров клюва с диаметром флуктуационного объема, в котором, согласно термофлуктуационной теории разрушения [1], происходят элементарные акты разрушения. Для ПММА
флуктуационный объем Уа =(0,5...1,4)-10"28 м3 [1], отсюда диаметр флуктуационного объема й ~ ^2Уа =(4,6...6,5)-10_1° м, т. е. диаметр
флуктуационного объема и длина клюва трещины практически одинаковы. Следовательно, в области клюва — микрообъема вблизи трещины (флуктуационного объема) — и разыгрываются элементарные акты разрушения. Точнее, в этом микрообъеме происходят заключительные акты «дорывания» напряженных связей. Первичные акты предразрушения происходят перед фронтом трещины разрушения в условиях, когда вследствие больших напряжений в полимерах развивается неупругая деформация (вынужденная эластическая деформация).
Напряжение сцепления в клюве трещины определяется формулой
d
жА2 E
с помощью которой можно оценить порядок величины напряжения сцепления. Используя для ПММА численные данные из работы [1]: Е =
= 3,93-109 Н/м2, ап = 3,9-10"2 Дж/м2, получаем /шах = 7,8-108 Н/м2. Следовательно, напряжение в клюве трещины имеет конечное значение, и для ПММА его можно принять равным 7,8-108 Н/м2. Для трещины нормального отрыва оно определяется нормальным напряжением ауу. Фронт трещины поперечного сдвига (трещина второго типа)
продвигается за счет скольжения берегов один относительно другого. Поэтому напряжение сцепления в клюве этой трещины — тангенциальное напряжение аху, хотя межчастичные связи при всяком виде
макроскопических напряжений испытывают только растяжение.
При фрактографическом исследовании поверхности излома разрушенных образцов наблюдали три зоны: зеркальную, матовую и перьевую (искромсанную) [1, 7]. Зеркальная зона соответствует медленному флуктуационному подрастанию трещины. Трещина на своем пути пересекает множество микрополостей размером от нескольких до десятка микрон. Взаимодействие трещины с этими микрополостями приводит к зарождению многих вторичных микротрещин, так что трещина разрушения продвигается как бы через заранее подготовленную разрыхленную среду, насыщенную микропорами и микротрещинами. В матовой зоне, соответствующей атерми-ческой стадии движения трещины с большой скоростью, возросшие локальные напряжения становятся достаточными для активизации изолированных микрополостей и их взаимодействия. Одновременно вместе с магистральной трещиной распространяется множество вторичных микротрещин. В перьевой зоне, которая завершает разрушение, этот процесс становится еще интенсивнее.
Процессы предразрушения, т. е накопление микроповреждений (разрыхление) перед фронтом трещины, протекают параллельно с развитием неупругой вынужденно-эластической деформации. При этом процессы предразрушения и вынужденно-эластического деформирования развиваются во времени с конечной скоростью.
Наиболее известной теорией, учитывающей наличие неупругой зоны перед трещиной, является теория Леонова—Панасюка— Дагдейла (ЛПД) [5, 8]. В этой теории рассматривается трещина нормального отрыва с пластической зоной перед ней. Теория ЛПД была создана для материалов, деформирующихся пластически, в первую очередь металлов. Пластическая зона моделируется тонкой прослойкой на продолжении трещины, заполненной нелинейно пластически деформированным материалом.
В полимерах аналогом пластической деформации является вынужденно-эластическая деформация. По внешним признакам эта де-
формация напоминает пластическую, хотя по физической природе это совершенно разные виды деформации. Поэтому перед фронтом трещины возникает зона вынужденной эластичности. Развитие зоны вынужденной эластичности происходит с преодолением внутреннего трения и выделением теплоты, что приводит к нагреванию материала вблизи трещины. Это, с одной стороны, ускоряет процесс накопления повреждений, а с другой — через релаксацию напряжения тормозит их. Микрорасслоение материала в зоне вынужденной эластичности, которое начинается, когда температура зоны приближается к температуре стеклования Тс, приводит к ориентационному упрочнению и также затормаживает процессы предразрушения.
По достижении напряжения, равного пределу вынужденной эластичности о"в для данной температуры, «размораживается» сегментальная подвижность звеньев макромолекул и начинает развиваться зона вынужденной эластичности ползучести при напряжении аа, связанная с разворачиванием свернутых конформаций макромолекул и большими перемещениями материала внутри зоны вынужденной эластичности. Вне зоны вынужденной эластичности, где напряжения меньше предела вынужденной эластичности, сегментальная подвижность блокируется межмолекулярным взаимодействием, и одно только тепловое движение не в силах его преодолеть. Поэтому вне эластической зоны материал испытывает только малые упругие деформации.
Напряжение в эластической зоне конечно и равно пределу вынужденной эластичности, соответствующему данной локальной температуре зоны. Форма и размеры зоны вынужденной эластичности описываются в полярных координатах формулой
r = -K-T (4 - — sin2 в-cos^l. (1)
1 4 J W
Здесь K — коэффициент интенсивности напряжений (КИН) (K1 — для трещины нормального отрыва и K2 — для трещины поперечного сдвига); в — полярный угол; ав — предел вынужденной эластичности.
На рисунке представлена рассчитанная по формуле (1) форма зоны вынужденной эластичности. Построение выполнено в безразмерных полярных координатах с использованием безразмерного полярного радиуса р = r K- 2 . Большая часть зоны вынужденной эла-/ -ка-
стичности вытянута вдоль берегов трещины. Примерно 20 % этой зоны располагается впереди перед фронтом трещины. В этой ее ак-
тивной части, являющейся продолжением клюва трещины, процессы предразрушения способствуют дальнейшему росту трещины.
Форма зоны вынужденной эластичности вблизи трещины поперечного сдвига при плоском напряженном состоянии
Чтобы найти площадь зоны вынужденной эластичности, проинтегрируем выражение (1) по области, ограниченной контуром зоны (см. рисунок). Получаем
s=11
Ж
Г к V
у
(2)
Площадь так называемой активной части зоны
s =
акт
0,22
Г к V
Ж
\°в у
(3)
Объем активной части получается умножением на толщину образца к:
V =
0,22h
Ж
' K V4
Ч^в у
(4)
Длина активной части зоны, находящейся перед фронтом трещины,
д=075
Ж
' K V2
у
(5)
Размеры зоны вынужденной эластичности, как следует из формул (4)—(7), определяются соотношением между КИН и пределом вынужденной эластичности сгв. В определении размеров зоны вынужденной эластичности ведущую роль играет КИН. Динамический КИН для движущейся трещины, вообще говоря, отличен от статического, и зависит от скорости распространения трещины. В работе [9] для движущейся трещины нормального отрыва при постоянном одноосном растягивающем напряжении получено выражение для динамического КИН в виде
К (I ,у) = КСтат (I) КДин (у) .
Здесь Кстат (I) — статический КИН, рассчитываемый обычными методами механики разрушения (I — длина трещины); Кдин (у) — динамическая поправка (V — скорость распространения трещины):
К л. \ _ 1 - Чск
К дин ^) _-
N1 - V
где фактор g слабо зависит от скорости V и принимает значения, близкие к единице; ек — скорость поверхностных рэлеевских волн; с1 —
скорость продольных упругих волн расширения-сжатия (ек < с1). С увеличением скорости V динамическая поправка Кдин (V) уменьшается, стремясь к нулю при V ^ ек. Эксперименты показывают, что максимальная скорость трещины не достигает рэлеевской скорости и по различным данным составляет (0,2...0,7)ек [9, 10].
По мере роста трещины увеличивается статический КИН, и размеры зоны вынужденной эластичности также увеличиваются, поскольку на больших расстояниях от фронта трещины напряжения достигают предела вынужденной эластичности. При увеличении локальной температуры вблизи трещины уменьшается предел вынужденной эластичности, т. е. на больших расстояниях, напряжения достигают этого предела. Поэтому размеры зоны вынужденной эластичности увеличиваются с ростом локальной температуры. Поскольку локальная температура также зависит от длины трещины, можно предположить, что при заданных внешних условиях размеры зоны управляются размером трещины.
При локальной температуре вблизи фронта трещины, превышающей температуру хрупкости Тхр, начинает развиваться вынужден-
ная эластическая деформация и появляется зона вынужденной эластичности, размеры которой сокращаются с уменьшением локальной температуры, и при температуре Тхр зона стягивается к клюву трещины. Поэтому при температуре ниже температуры хрупкости трещина растет по механизму хрупкого разрушения без зоны вынужденной эластичности. По достижении температуры квазихрупкости Ткхр в зоне вынужденной эластичности появляется собственная
структура: материал зоны расслаивается на тяжи, разделенные между собой, и зона преобразуется в крейз (трещину серебра). Для справки приведем значения температуры Тхр, Ткхр, Тс для двух полимеров [1]:
ПММА Тхр = -29 °С, Ткхр = 52 °С, Тс =100 °С;
ПЭТФ Тхр = -20 °С, Ткхр = 45 °С, Тс = 80 °С.
хр ' кхр '
Если локальная температура выше температуры Тхр или становится таковой в процессе медленного подрастания трещины, перед ее фронтом появляется зона вынужденной эластичности, и трещина прорастает через нее, одновременно двигая ее вперед.
Для аморфных полимеров, таких как ПММА, ПС, ПЭТФ и других, ниже температуры стеклования Тс выделяют три прочностных состояния: хрупкое, квазихрупкое и нехрупкое (вязкоупругое) [1]. Они разделяются температурами хрупкости Тхр и квазихрупкости
Ткхр. В каждом состоянии распространение трещины имеет свою
специфику, что отражено в названиях трещин: хрупкая, квазихрупкая и крейз. Ниже значения Тхр реализуется хрупкая трещина без эластической зоны, в диапазоне значений Тхр ... Ткхр перед фронтом трещины появляется зона вынужденной эластичности, а выше значения Ткхр она преобразуется в крейз и трещина прорастает через него за
счет накопления в тяжах повреждений и их последующего разрыва.
В работе [10] отмечается, что для развития в пластичных материалах (в первую очередь в металлах) хрупкого или нехрупкого разрушения важную роль играет скорость трещины. Пока скорость мала, пластическая деформация успевает развиться и трещина растет через пластическую зону. Когда скорость достигает достаточно большого значения, пластическая деформация не успевает развиться за трещиной и наступает хрупкий разрыв, который затем может снова превратиться в пластический, если освобождающейся упругой энергии окажется недостаточно для обеспечения нужной скорости трещины. Этот релаксационный фактор еще в большей степени проявляется в полимерах, которые в отношении механических свойств являются
вязкоупругими средами с сильной зависимостью механических характеристик от скорости и длительности (и частоты) воздействия, а также от температуры.
Моделирование на компьютере движения трещины в условиях существования впереди нее зоны вынужденной эластичности приводит к следующим результатам. Вынужденная эластическая деформация развивается во времени с конечной скоростью, зависящей от напряжения и температуры. Скорость формирования зоны вынужденной эластичности характеризуется временем т запаздывания вынужденной эластической деформации ползучести, зависящим от интенсивности межмолекулярных взаимодействий, тормозящих процесс разворачивания свернутых конформаций макромолекул, а также от температуры и внешнего напряжения. Введем безразмерный характеристический параметр % = —, где /кр — критическая длина на флук-
1кр
туационной стадии движения. Произведение —т определяет среднее расстояние, на которое может продвинуться трещина за время развития вынужденной эластичности деформации. Если % ^ 1, за время т
трещина может пройти путь, намного меньший ее критического пути, т. е. подавляющую часть своего пути трещина будет двигаться в условиях уже развившейся вынужденной эластической деформации впереди нее. Поэтому условие % « 1 определяет наличие перед трещиной
сформировавшейся зоны вынужденной эластичности.
Развитие вынужденной эластической деформации происходит с совершением работы против сил внутреннего трения и сопровождается выделением теплоты. Можно показать, что пока скорость трещины невелика, этот поток теплоты успевает отводиться во вне, и движение трещины является изотермическим. По мере роста трещины ее скорость возрастает, напряжения вблизи фронта увеличиваются, выделяющаяся теплота уже не успевает полностью отводиться и происходит постепенный переход к адиабатическому движению трещины, зона вынужденной эластичности нагревается, увеличивается частота элементарных актов разрушения, время развития вынужденной эластической деформации уменьшается, уменьшается и предел вынужденной эластичности сгв. Все эти факторы противоречивым образом влияют на движение трещины и эволюцию зоны вынужденной эластичности. В конце медленной стадии движения трещины из-за адиабатического разогрева ее скорость резко, почти скачком, увеличивается, приближаясь к предельной скорости атер-мической стадии. В это время характеристический параметр % становится больше единицы, и вынужденная эластическая деформация
уже не успевает за трещиной, размеры зоны вначале «замораживаются», а с переходом к атермической стадии движения трещины зона вынужденной эластичности исчезает, и заключительный этап разрушения происходит по механизму хрупкого разрушения — трещина «летит» адиабатически с высокой локальной температурой в ее головной части. Следовательно, условие X » 1 определяет отсутствие
зоны вынужденной эластичности.
При температуре выше температуры хрупкости в нагруженном материале в окрестности трещины появляется зона вынужденной эластичности, заполненная нелинейно деформированным полимерным веществом. Материал зоны подвергается вынужденной эластической деформации ползучести, связанной с большими макроскопическими перемещениями. Вне зоны материал образца испытывает только упругие деформации. На микроскопическом уровне развитие вынужденной эластической ползучести, обусловленной возникновением в силовом поле вынужденной сегментальной подвижности, проявляется в изменении конформаций макромолекул от свернутых к вытянутым без их разрыва и перемещения как целого. При этом изменяется (уменьшается) конформационная энтропия макромолекул. В этом смысле вынужденная эластическая деформация имеет энтропийное происхождение. Но конформационные преобразования при вынужденной ползучести происходят с преодолением блокирующего действия межмолекулярных сил, поэтому изменяется не только энтропия, но и внутренняя энергия макромолекул. При этом характер вынужденной эластической деформации не чисто энтропийный. Тормозящее действие межмолекулярных сил проявляется как внутреннее трение, на преодоление которого затрачивается работа, рассеивающаяся в виде потока теплоты. В связи с этим при формировании зоны вынужденной эластичности она нагревается.
В зоне вынужденной эластичности деформация имеет две составляющие — упругую, которая возникает практически мгновенно, и запаздывающую — вынужденно-эластическую. Упругая составляющая обусловлена растяжением химических связей и валентных углов, в первую очередь, в цепях главной валентности полимерных макромолекул. Эластическая составляющая появляется только выше температуры хрупкости, так как только при этой температуре проявляется сегментальная подвижность как вынужденная ползучесть в стесненных условиях малого объема зоны. Ее развитие растянуто во времени, а эластическая составляющая значительно превосходит упругую. Все это позволяет записать тензор деформаций в зоне вынужденной эластичности в виде
С (М, ^ = етр (М) + ет (М, г)л(т - ТХр), К] = 1, (6)
где г][Т - Тхр ) — функция Хевисайда, зависящая от пространственных координат произвольной точки М зоны и от времени. Формула (6) действует только в пределах зоны эластической деформации, которая при понижении температуры до температуры хрупкости Тхр
стягивается к клюву трещины. При температуре ниже Тхр зона вынужденной эластичности отсутствует, и в этой формуле остается только первое слагаемое. Тензор вынужденной эластической деформации разложим на шаровую и девиаторную части:
8/эл — /л ^ з 0эл^г/',
где 0 эл = 5ккэл — след тензора эластической деформации,
8кк = 8XX + ^ + 8; 5г/ — символ КPонекеPа,
/.. =5.. -10 .
г/эл г/эл 3 эл г/
Вынужденная эластическая деформация развивается в стесненных условиях малого объема окрестности трещины, за пределами зоны вынужденной эластичности материал деформируется только упруго. Поэтому при деформации этой зоны изменение объема не происходит, и существенным является только девиатор деформации. Уточним формулу (6):
eV
(М, t) = %уупр (М) + 1т (М, t) 0 (T - гхр).
В пределах зоны вынужденной эластичности (когда она есть) устанавливается напряжение, равное пределу вынужденной эластичности. Для нормальных и тангенциальных компонент напряжений предел вынужденной эластичности, вообще говоря, различен. Как правило, для сдвига он меньше, чем для растяжения [1]. Поэтому, например, для трещины поперечного сдвига размеры зоны вынужденной эластичности определяются пределом вынужденной эластичности на сдвиг, и эластическая составляющая деформации появляется, когда тангенциальное напряжение достигает предела вынужденной эластичности для данной температуры.
В тензоре упругих деформаций 8г/упр существенны и шаровая
часть, и девиатор, так как при упругой деформации изменяются и объем, и форма. Деформация химических связей и валентных углов в макромолекулах — это всегда упругое растяжение, но оно вызывается всеми компонентами тензора напряжений, поскольку на микро-
скопическом уровне всегда имеются химические связи в макромолекулах, растягиваемые всеми компонентами макроскопического напряжения.
В процессе формирования эластической зоны возникают вязкие или диссипативные напряжения, вызванные силами внутреннего трения. Вязкие напряжения зависят от скорости вынужденной эластической деформации и исчезают, когда развитие этой деформации завершается. Работа, производимая вязкими напряжениями, определяет потери энергии при деформировании. Эта энергия рассеивается в виде теплоты, поэтому при формировании зоны вынужденной эластичности она нагревается. Потери энергии на трение в единице объема в единицу времени определяются так называемой диссипативной функцией Я (М, t), а именно выражение 2Я (М, t)АУАt определяет
потери энергии в малом объеме АУ около точки М за малое время Аt приблизительно в момент времени t. Мощность тепловыделения в единице объема
й (М, t ) = 2Я (М, t).
Количество теплоты, выделяемой во всем объеме V зоны вынужденной эластичности за промежуток времени t2],
t2
й = {{2Я(М, t)сУЖ. (7)
tl V
Диссипативная функция Я (М, t) является квадратичной функцией тензора скоростей эластической деформации ¿^. Тензор вынужденной эластической деформации представляет собой сумму девиа-тора и шаровой части:
Я = »(т )2 + 2Фэл )2.
Здесь ц — вязкость формы (на сдвиг), % — вязкость объема. Поскольку в рассматриваемом случае вынужденная эластическая деформация определяется девиатором ], имеем
Я (М, t ) = ^( 4эл )2.
Тогда мощность тепловыделения в произвольной точке М зоны вынужденной эластичности
Q(M, t) = 2ß[im (M, t)
(8)
Зависимость девиатора эластической деформации от времени определяется реологическими свойствами полимерного вещества. Для оценки этих свойств можно взять реологическую модель Кельвина (линейное стандартное тело). Эта модель основана на сочетании мгновенной упругости и запаздывающей упругости, обусловленной вязкостью материала. Мгновенная упругость характеризуется мгновенным модулем упругости Е, а запаздывающая — вязкостью ц и временем запаздывания т. В одномерном случае ползучесть материала Кельвина описывается формулой
(
-Л
1 - e
(9)
Здесь а0 — приложенное напряжения; е^ — равновесное значение запаздывающей деформации:
т — (
V
Первое слагаемое в формуле (9) определяет мгновенную упругую деформацию и соответствует первому члену формулы (6), а второе — запаздывающую вынужденную эластическую деформацию, соответствующую второму члену формулы (6). С учетом формулы (9) запишем девиатор вынужденной эластической деформации в виде
(
-Л
1 - e т
v J
V
(T - Тхр ).
//эл (М, г ) = /^{М, г)
Тогда скорость этой деформации
1 —
4эл (М, г ) = - /г^(м, г )е *л(Т - Тр).
(10)
Подставив соотношение (10) в формулу (8), получаем мощность
тепловыделения
Q (М, t ) = Т- e т [^(М, t)^v(T - Тхр ).
-t
(11)
Это означает, что в процессе формирования эластической зоны в каждой точке ее объема действует тепловой источник с такой мощ-
ностью. Количество теплоты, выделяющейся за время формирования зоны вынужденной эластичности, можно получить, если подставить (10) в формулу (7) и провести предусмотренное интегрирование:
Q = 0,86F - ГХр).
Здесь V — объем зоны вынужденной эластичности; j — средняя
равновесная эластическая деформация.
По завершении формирования зоны вынужденной эластичности в ней устанавливается однородное напряжение, равное пределу вынужденной эластичности. При этом количество выделяющейся теплоты в зоне вынужденной эластичности
Q = 0,86V ^Г - ГХр ).
Таким образом, определены размеры клюва трещины и напряжение, возникающее в нем, а также форма и размеры зоны вынужденной эластичности, возникающей в аморфных твердых полимерах перед трещиной разрушения при температуре выше температуры хрупкости.
Смоделировано движение трещины с зоной вынужденной эластичности. Исследованы процессы тепловыделения при формировании зоны вынужденной эластичности.
СПИСОК ЛИТЕРАТУРЫ
1. Бартенев Г. М. Прочность и механизм разрушения полимеров. М.: Химия. 1984.
2. Валишин А. А., Карташов Э. М., Тишаева С. Д. // Ученые записки МИТХТ, 2003. № 7. С. 40-45.
3. Баренблатт Г. И. // Проблемы механики сплошных сред. Изд-во АН СССР. М.-Л., 1961.
4. Баренблатт Г. И. // Прикладная математика и механика. 1964. № 4. С. 630-643.
5. Панасюк В. В. Предельное равновесие хрупких тел с трещинами. Киев: На-укова думка. 1968.
6. Ландау Л. Д., Лифшиц Е. М. Теория упругости. М.: Наука. 1978.
7. Партон В. З. Механика разрушения. М.: Наука. 1990.
8. Основы экспериментальной механики разрушения / И.М. Керштейн, В. Д. Клюшников, Е.В. Ломакин и др. М.: Изд-во МГУ, 1989.
9. Костров Б. В. // Прикладная математика и механика. 1974. № 3. С. 551-560.
10. Fe hieb eck D. K., O row an E. O. // Weld. Journ. Ress. Suppl. 1955. 34. Р. 570-575.
Статья поступила в редакцию 03.07.2012.