CC BY
50
В.В. Стружанов, В.И. Миронов, К.А. Тарташник
ОБ ОДНОМ ПОДХОДЕ К РАСЧЕТУ ДОЛГОВЕЧНОСТИ ПРИ ЦИКЛИЧЕСКОМ НАГРУЖЕНИИ
Предлагается новый подход к расчету долговечности элементов конструкций при циклическом нагружении, основанный на взаимосвязи разрушения и потери устойчивости процесса деформирования. Для этого введены в рассмотрение неустойчивые состояния материала, характеризуемые падающими ветвями полных диаграмм деформирования, деградирующих в ходе переменного нагружения. Основные положения методики проиллюстрированы на примерах усталостного разрушения однонаправленного композита (пучка волокон) при одноосном циклическом деформировании.
Введение. Элементы машин, как правило, работают в условиях циклического нагружения. Достоверная оценка их долговечности является одним из основных факторов, влияющих на создание высоконадежной техники. В настоящее время в расчетах на долговечность наиболее распространено использование критериев, которые утверждают, что разрушение произойдет тогда, когда в какой-либо точке тела в результате циклического нагружения некоторая скалярная величина достигнет своего критического значения [1]. Однако возможен и другой подход, заключающийся в том, что разрушение связывается с моментом потери устойчивости процесса деформирования конструктивного элемента. Для статического нагружения это утверждение обоснованно в работах [2,3] и опирается на введении в рассмотрение неустойчивых состояний материала (стадии разупрочнения), характеризуемых падающими ветвями полных диаграмм деформирования. В данной работе сформулированный в [2,3] подход распространяется на задачу усталостного разрушения. В основу положено вырождение полной диаграммы квазистати-ческого деформирования в ходе циклического нагружения, которое в расчетах усталостной долговечности, по-видимому впервые, было использовано в [4].
Критерии разрушения образца при квазистатическом растяжении в испытательной
мещение правого конца образца, скрепленного со стержнем 2, обозначим V. Полагаем, что свойства образца определяет полная диаграмма деформирования в координатах растягивающее усилие д - удлинение V. Основной особенностью данной диаграммы является наличие ниспадающей до нуля ветви, т. е. после достижения предела прочности дв, деформирование образца происходит в неустойчивом режиме, когда усилие падает, а удлинение растет. Если удается сохранить равновесие, то разделение образца на части (разрушение) происходит в конце падающего участка. Отметим, что данные диаграммы неоднократно были построены на установках различного типа [5,6]. Для дальнейшего анализа конкретный вид падающего участка полной диаграммы не имеет принципиального значения. Существенен лишь факт его существования.
Следуя [2,7], запишем выражение для полной энергии системы образец-машина для случая мягкого нагружения. В результате имеем
Р и с. 1. Схематическое представление системы образец (1) - машина (2)
машине. Используя аппарат, приведенный в [2], и опуская необязательные для данного случая математические подробности, определим критерии разрушения образца, растягиваемого в некоторой испытательной машине. Систему образец-машина схематически представим в виде последовательного соединения образца 1 и передающего усилие упругого стержня 2 с жесткостью с, равной жесткости машины (рис.1). Нагружение осуществляем, прилагая к свободному концу стержня 2 силу р (мягкое нагружение), либо задавая перемещение и (жесткое нагружение). Пере-
V
и
Здесь первое слагаемое - энергия деформации образца, второе - потенциальная энергия упруго деформированных деталей машины (стержня 2), третье - работа внешней силы р, взятая со знаком минус. Задаваемые величины с, р играют роль параметров управления системой, переменные V, и - параметров состояния. Так как при активном нагружении система является консервативной, то для всевозможных значений параметров с и р все положения равновесия системы определяются следующими уравнениями:
йП = д - с(и - v)= 0, сйр = с(и _ V)_ р = о. (1)
^ йи
Совокупность критических точек (V, и, с, р) функции П (решений уравнений (1)) образует так называемое многообразие катастрофы [7].
В случае жесткого нагружения параметрами управления являются и и с, а параметром состояния - переменная V. Критические точки находятся из решения первого уравнения системы
(1). Их отображение ^,и, с) ® (V, и, с-1) также определяет многообразие катастрофы, которое для многих технических приложений с двумя параметрами управления и одним параметром состояния представляется катастрофой сборки Уитни [7].
Особую роль в данном анализе играют вырожденные критические точки, где происходит смена устойчивых положений равновесия системы на неустойчивые и наоборот [7]. Вырожденные критические точки обращают в нуль детерминант матрицы Г ессе функции П [7]. Используя это свойство для мягкого нагружения получаем уравнение
:>2^ ^2^ ( л2р ^
= 0.
д 2П д2п
йы йу
йу2 йи2
Отсюда в вырожденных критических точках выполняется равенство
йд
В случае жесткого нагружения гессиан функции энергии системы имеет одну компоненту
= 0.
(2)
д Пду2 • Приравнивая ее к нулю, получаем условие
йд йу
(3)
Р и с. 2. Бифуркационные множества при мягком нагружении
выделяющее из множества критических точек вырожденные точки. Решая совместно систему
(1) и уравнение (2) и первое уравнение в системе (1) и уравнение (3), получаем в пространствах уравнений бифуркационные множества В\ и В2 для мягкого и жесткого нагружений (рис.2,3). Известно [7], что потеря устойчивости процесса деформирования системы происходит тогда, когда путь нагружения в пространстве управления пересекает кривую В2 (путь 1 на рис.2 и 3). В случае мягкого нагружения после пересечения прямой В2 у системы исчезают устойчивые положения равновесия. Следовательно, образец в данный момент разрушается динамически на пределе прочности. В случае жесткого нагружения пересечение кривой
В2 приводит к скачкообразному переходу из одного устойчивого положения равновесия в другое. В зависимости от величины скачка образец может либо сохранить ограниченную несущую способность при существенных повреждениях микротрещинами, либо разрушиться полностью.
йд.
Если величина с достаточно велика (всегда выполняется неравенство
/> _ с), то путь нагружения может не пересекать кривую В 2 (путь 2 на рис.3). Тогда разрушение образца происходит равновесно в конце падающего участка диаграммы деформирования.
Очевидно, что в общем случае бифуркационные множества, обладающие свойствами множества В2, играют определяющую роль в нахождении момента потери устойчивости деформирования механических систем, так как здесь впервые происходит смена устойчивого положения равновесия системы на неустойчивое. Вне зависимости от числа бифуркационных множеств, присущих конкретной системе, будем обозначать такие множества символом В2. Для
95
Р и с. 3. Бифуркационные множества при жестком нагружении
определения аналитического выражения кривой В2 в случае мягкого нагружения из уравнений (1), (2) исключаем параметры состояния V, и. Получаем уравнение прямой в пространстве управления (множество В2)
Р = Рв = )> (4)
где V5 - наименьший корень уравнения
(2), т. е. д^5 )= дВ . При жестком нагружении исключаем из первого уравнения системы (1) и уравнения (3) параметр состояния V. Находим уравнение кривой в
пространстве управления (множество В2)
иВ = и(с) = с-1д (ук)+ Vй , где ук = ук (с) - значения наименьших корней уравнения (3) в зависимости от параметра с.
Отметим, что функция и(с) определена, если с е [0,с*], с* =
ёд/ёу , так как только при
данных значениях с уравнение (3) имеет решение.
Из приведенного анализа вытекает, что разрушение образца связано с потерей устойчивости процесса деформирования всей системы образец-машина. Поэтому критериями разрушения могут служить критерии потери устойчивости. Таким образом, при мягком нагружении критерием разрушения выступает равенство (4). Критерием разрушения для жесткого нагружения является равенство
и = иВ (5)
при конкретном значении параметра с.
Циклическое нагружение образца. Зададим отнулевой цикл нагружения системы с максимальным усилием р = ра, величина которого такова, что материал образца (по крайне мере на начальных циклах) работает в области упругости. Полагаем, что процесс циклирования состоит из последовательности п статических растяжений. В пространстве управления в каждом цикле данное нагружение происходит по пути 1 до точки А (рис.2). На некотором цикле растяжения ( п=^р ) образец разрушается. Как следует из рассуждений, приведенных выше, при разрушении путь нагружения должен пересечь бифуркационное множество В2. Так как ра=сот1, то такое пересечение возможно лишь в том случае, когда прямая В2 займет В2', т. е. опустится до точки А (рис. 2). Положение бифуркационных множеств в пространстве управления зависит от характеристик полной диаграммы деформирования образца. Следовательно, в ходе циклиро-вания эта диаграмма изменяет свой вид (деградирует). Данное изменение описывается некоторой двухпараметрической функцией ц=ц(у, ра, п). Тогда из условия (4) следует критерий усталостного разрушения
Ра = $(У , Ра , п^ (6)
отражающий снижение предела прочности образца до уровня максимального усилия цикла. Здесь V5 = V5 (ра, п) - наименьшее решение уравнения (2), в котором д = д(у, ра, п). Разрешая уравнение (6) относительно п получим число циклов до разрушения Ир, принадлежащее, очевидно, диаграмме Велера.
Рассмотрим теперь стационарное циклирование с контролем максимального перемещения иа. Очевидно, что и при этом варианте циклирования происходит вырождение диаграммы д^),
определяемое теперь двухпараметрической функцией д = д(V, иа, п). Пусть путь нагружения в пространстве управления для каждого цикла задает прямая, лежащая на пути 1 до точки К (рис.3). Образец разрушится, если бифуркационное множество В2, изменяя свое положение от цикла к циклу, при некотором п=Ии пройдет через точку К (кривая В2" рис.3). Используя условие (5), получаем критерий усталостного разрушения в виде уравнения
= иВ (с, иа, п) = с 1д(уh, иа, п)+ V"
(7)
и
а
где Vй = Vй (иа, п) - наименьшее решение уравнения (3) при д = д(у, иа, п), которое существу-
ет для с е [0, с*], с * = тГ dд(v, иа, n)/dv
. Решение уравнения (7) относительно п при фик-
сированном параметре с дает искомую долговечность Ли.
При данной схеме циклирования возможен еще вариант, когда уравнение (7) не имеет решения, т. е. в процессе деградации диаграммы д^) падающая ветвь становится слишком пологой (уравнение (3) не имеет решения для заданного с). Тогда образец разрушается равновесно с падением несущей способности и, следовательно, максимального усилия цикла до нуля. В этом случае бифуркационное множество В2 уходит от пути нагружения.
Таким образом, в основу проведенного анализа, как в случае квазистатического растяжения, так и в случае циклического нагружения, положен единый подход, основанный на определении разрушения образца как катастрофы (потери устойчивости процесса деформирования системы образец - машина). Отличие заключается лишь в том, что при статическом нагружении система приближается к катастрофе, а при циклическом нагружении катастрофа надвигается на систему, причем иногда это приближение не заканчивается катастрофическим разрушением (катастрофа «проходит стороной»).
Вырождение диаграммы деформирования. В качестве примера рассмотрим пучок одинаковых упруго-хрупких волокон, расположенных параллельно друг другу. При продольном растяжении деформация волокон равна деформации пучка е. Прочность волокон равномерно распределена на замкнутом интервале напряжений [у, 52]. Поперечным взаимодействием пренебрегаем.
Сначала найдем связь между условным напряжением пучка и деформацией е при квазиста-тическом растяжении. Имеем [8]:
о(е) = Е е(1 -ю(е)),
где Ю = {0, 0 <е< е1; (Ее- 51 )/(у2 - у), е1 <е< е2} - поврежденность, равная объемной
2Х-,. Качествен/ Е
доле разрушенных волокон; Е - модуль упругости волокон; е1 = , е2
ный характер диаграммы о(е) показан на рис.4 (кривая 1) для случая, когда у2 > 2 у1.
Подвергнем пучок пульсирующему циклическому нагружению. Полагаем, что характеристики упругого деформирования волокон не изменяются, их долговечности не зависят от уровня напряжений. Для любого множества волокон, пределы прочности которого расположены в конечных отрезках | \ . л'; | сг | л,, л\ ]
(у,. > 51, sJ < у2, sJ > у.) долговечность равномерно распределена на замкнутом интервале числа циклов нагружения Ль Л2].
Кроме того, разрушенные в ходе циклирования волокна не изменяют интервала [у, у2]. Данные допущения весьма существенны. Однако они позволяют существенно упростить задачу, не искажая качественной картины.
Пусть в ходе циклирования поддерживается постоянное значение амплитуды деформации пучка е а (е а < е1). Тогда зависимость между условным напряжением и деформацией пучка приобретает вид
а (е,п) = Е (1 - к(п)) (1 - (0(8))
Р и с. 4. Вырождение статической диаграммы пучка при постоянном уровне амплитудной деформации
где
к(п) = {0, 0 < п < Л1; (п - Л1 )/(Л2 - Л1), Л1 < п < Л2} -
циклическая поврежденность,
равная объемной доле усталостно разрушенных волокон. Характер изменения статической диаграммы растяжения пучка показан на рис.4. Очевидно, что полное разрушение пучка произойдет при нулевом напряжении, когда диаграмма выродится в прямую линию.
Если поддерживать постоянным амплитудное условное напряжение Оа (оа < 51), то характер изменения статической диаграммы будет более сложным. Вместе с падением модуля пучка
из-за усталостного разрушения волокон будет монотонно расти деформация є a. В этом случае, после того, как деформация є a достигнет значения e1, в рамках цикла неизбежен выход на нелинейную часть диаграммы деформирования. Тогда
О(є, n)= En є(І — Юе), (8)
где En = E(l — к(п))(і — « (є)) - модуль упругости пучка после n-го цикла нагружения;
Wa ={ 0 Єa (n) < Єі; (E є a (n) — S1 V(s2 — S1 X Єa (n) > Єі } - поврежденность от растяжения в
рамках цикла; ЮЕ = Ю , если Єa < e1, и we ={0,e <ea;(Ee — Eea(n))/(s2 — Eea(n)),e >ea} при Єa > e1 - поврежденность, возникающая при послецикловом квазистатическом растяжении.
<та
Р и с. 5. Вырождение статической диаграммы пучка при постоянном уровне амплитудного условного напряжения
Характер изменения статической диаграммы растяжения пучка показан на рис.5. В ходе вырождения статической диаграммы происходит падение предела прочности пучка. Разрушение произойдет динамически тогда, когда Оa = ОB (n) и Є a (n) = ЄB , где
B Є /
Є = у2 - деформация, отвечающая пределу прочности пучка при квазистатическом растяжении.
Реализация в течении длительного времени приведенных схем нагружения для элемента конструкции, находящегося в составе механической системы, практически невозможна. Поэтому рассмотрим еще циклическое нагружение при монотонно возрастающей амплитуде деформации є a = є a (n) . Вырождение диаграммы и в этом случае определяет формула (8) с той разницей, что Є a является известной функцией. На рис.б на качественном уровне показано изменение диаграммы пучка О(є) (кривая 1) в зависимости от возрастающего числа циклов нагружения и возрастающей функции £и(п). Кривые 2, 3 (n3 > n2) при £а < ех, кривые 4, 5 (п4 > пъ) для £а > ех.
Р и с. 6. Вырождение статической диаграммы пучка В ходе вырождения статической диа-при монотонном возрастании амплитудной дефор- граммы происходит падение предела проч-
мации £а ности пучка. Когда £,а {п) < £/;. где
ЄB = e2 /2 - деформация, отвечающая пределу прочности пучка при квазистатическом растяжении, падение прочности пучка происходит в соответствии с законом
ОB (n) = 0,25 s2 (l — K(n))/(s2 — s1). (9)
Если же Єa (n) > єB, то ОB (n) = En єa (n). Отметим, что полученные соотношения не зависят
от длины волокон и площади поперечного сечения пучка.
Циклическое растяжение пучка в испытательной машине. Поместим теперь пучок единичной длины и единичной площади поперечного сечения вместо образца в испытательную установку (рис.1) и осуществим пульсирующее циклическое нагружение с постоянным максимальным усилием цикла pa < s1. При этом на пучок передается усилие qa = pa и среднее напряжение О a по величине равно pa . Пучок получает удлинение va , по величине совпадающее с Є a. В ходе циклирования в связи с падением модуля пучка параметр £a постепенно возрастает
по закону
е а (п)= Ра1 Еп-1 . Разрушение произойдет тогда, когда величина еа достигнет значе-
ния e (в данный момент выполняется равенство (6)). Используя выражение (9), из формулы (6) находим число циклов до разрушения N = N2 - 4pa (s2 - s1) (N2 - N1)/S). Заметим, что в момент разрушения к < 1. В пучке еще сохранились целые волокна, которые при n = N разрушаются динамически.
Исследуем теперь циклическое нагружение с контролем максимального перемещения ua=const, величина которого удовлетворяет неравенству
ua < (e1 + Ec le1) l, (10)
гарантирующему, по крайне мере на начальных циклах, деформирование волокон в области упругости (1=1 мм).
Сначала найдем значение c = inf do/de , где о задано выражением (8). Получаем
e
c* = Es2 (N2 — n)/[(N2 — N1 )(s2 — s1)] (касательный модуль полной диаграммы в данной задаче всегда имеет наименьшее значение при e = e2). Исходя из выражения для c , зададим жесткость упругого элемента установки выражением c = Es2/[k(s2 — s1)], k > 1. Тогда наименьшее решение уравнения (3) (оно и единственное) равно
e h = s^_ N2(k +1) - kn - N1 2E k(N2 - n) .
Отметим, что это решение существует только для n < [N2 (k -1) + Nj ]/k. Далее, считая e2 = a e1 (a > 2), из неравенства (10) имеем предельное значение
1-k
1
/-1
Подставляя выражение для е и данное Ыа в формулу (7), получаем уравнение
1 + k
1
1
a
_ а [к (ы2 - п)+(^2 - N )]2 у / 2 2к(N2 - п)( - N) ,
которое не имеет решения, так как правая часть всегда больше левой.
Таким образом, в условиях данной задачи при циклировании катастрофы не происходит. Пучок разрушается равновесно и его долговечность равна N2 (катастрофа уходит с пути нагружения).
Долговечность защитного покрытия трубы. В заключении исследуем на долговечность тонкое (единичной толщины) защитное покрытие толстостенной трубы, выполненное из волокнистого композитного материала. Волокна расположены перпендикулярно к образующей трубы. Влиянием связующего пренебрегаем, т. е. полагаем, что свойства покрытия идентичны свойствам пучка, рассмотренного выше. Труба нагружается пульсирующим внутренним давлением с максимальной величиной давления в цикле рт , при котором напряжение в композите меньше 51. Считаем, что циклические свойства материала трубы значительно выше, чем у композита, то есть циклическое разрушение покрытия происходит раньше, чем в трубе появляются повреждения, приводящие к изменению механических свойств.
В случае плоской деформации, имеющем место в данной задаче, поведение системы покрытие-труба характеризует полная энергия, определенная только для одного сечения. Полагая площадь поперечного сечения покрытия равной единице, имеем
П с _П к + П Т - 2 я а | р dw ,
0
я Е (1 -к) V 2Ь-1,
ь0
где Пk = 2pb Jo0(e0,n)de0
v <be;
be1 < v <be2;
[ + g 2v 2 + g3 ](1 -к)
[1 (be2 )3 + g 2 (be2 )2 + g3 ](1 -KX V >be2 -
энергия деформации покрытия;
ua = e1
о
П т = 2 р | ( е'г + 0 е'0 )г йг = ш1
V2 -ш2уи + ш3и2 -
потенциальная энергия деформации трубы. Здесь а, Ь - внутренний и внешний радиусы трубы; 00 (е0, п) - напряжение в покрытии, заданное зависимостью (8); е0 = ^ ~ деформация покрытия; и, V - радиальные перемещения точек внутренней и внешней поверхностей трубы; 0Г, 00 , е'г, е0 - радиальные и тангенциальные напряжения и деформации в трубе, которые получаем из известного решения задачи Ляме при заданных на границах перемещениях и и V; р -квазистатически возрастающее давление;
g1 = -2рЕ2 [зЬ2 (^2 - у )]^, g2 = р&2 [Ь(^2 - у )]-1, g3 = —2рЬ^13 [бЕ(^2 - ^ )]-1, ш1 = п а [(1 - 2п') а2 + Ь2 J(Ь2 - а2) , ш2 = 4аЬ (1 -п') па (Ь2 - а2) ,
ш3 = ра[а2 + (1 - 2у')Ь2]ь2 - а2)1, а = Е'[(1 + у')(1 - 2у')]-1;
Е', V - модуль Юнга и коэффициент Пуассона материала трубы. Параметрами управления являются р и й = Ь — а, параметрами состояния - и, V.
Отметим, что в ходе циклирования деформация покрытия возрастает от цикла к циклу по закону
е а = 2 р аш2 рш/ [4 р 2а 2Ь(1 - 2у') + 4 рЕш3(1 -к)] и может превысить значение е1. Например, для
Рш =еШ [4р2а2Ь (1 - 2у') + 4рЕш3 ]2раш2, (11)
где еш = Ь е1 (Ь £ 1), это произойдет после числа циклов
п = Н2 - |р - (1 - Р)р а 2Ь(1 - 2/)/Еш3 ](Н2 - Их).
Тогда энергию деформации покрытия при статическом нагружении будет определять выражение
[Е (1 -К(П))(1 - иа (еа ))рУ "ЬЛ У £ Ь е0а (е0а =еа );
П * =
Іg1 (Ь е0а )3 + g2(Ь е0а )2 + g3 ] - К(П)), У > Ье2,
где g1 = g2 = g2, g3 = -рЬЕ"еа [3(*2 - *1 )]-1.
Бифуркационные множества в пространстве управления определяются решениями системы уравнений
ЭП „ 0 ЭП^ = 0 Э2П „ Э2П
Эу Эи Эу Эи
Э2 П с Эу Эи
0, (12)
последнее из которых получается приравниванием к нулю детерминанта матрицы Гессе функции П с.
Построим множества В2 для конкретных числовых значений: у = 610 МПа, *2 = 1250 МПа, Е = 2,9-105 МПа (бериллиевые волокна [9]), Е' = 0,7-105 МПа (латунь), V' = 0,3, Ь = 50 мм, И1 = 6-105, М2 = 10б. Вычисления проводим, изменяя внутренний радиус трубы от а = 10 мм до а = 49 мм. Подставляя выражение для П с в формулы (12) и исключая параметры состояния и и V,
находим критические величины р для заданных значений циклов нагружения п,. По ним для различных а строим бифуркационные множества В2, где происходит переход от устойчивых положений равновесия к неустойчивым. Эти множества изображены на рис.7, где кривая 1 построена для п1 = М, кривая 2 для п2 = 7-105, кривая 3 для п3 = 8-105, кривая 4 для п4 = 9-105. Далее по формуле (11), где полагаем егпа = е1, рассчитаем для различных внутренних радиусов
предельные значения рш , обеспечивающие на начальных циклах работу покрытия в области упругости. Получаем рш = 32,39; 27,7; 23,32; 19,17; 15,37; 11,77 МПа соответственно для а = 43, ... , 48 мм. Пустьрш = 10 МПа. Точки пересечения соответствующей прямой с бифирку-
ционными кривыми (рис.7) определяют толщину трубы и число циклов нагружения, при которых реализуется катастрофическое разрушение покрытия. Например, если а = 48,57 мм, то катастрофа произойдет при п4 = 9-105. Данные результаты кроме того показывают, что всегда существуют такие размеры трубы, которые обеспечивают равновесное (не катастрофическое) разрушение покрытия. Для рш = 9 МПа равновесное разрушение покрытия обеспечивает величина а < 48,57 мм. Если же а £ 43,5 мм, то равновесное разрушение происходит при любом усилии в цикле.
Р, МПа
Р и с. 7. Бифуркационные множества для покрытия при циклическом нагружении трубы
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Коллинз Дж. Повреждение материалов в конструкциях. Анализ, предсказание, предотвращение. М.: Мир, 1984. 624 с.
2. Стружанов В.В. О разрушении диска с ослабленной центральной зоной, //Изв. АН СССР. МТТ. 1986. №1. С. 135-141.
3. Стружанов В.В. Об одном подходе к исследованию разрушения механических систем // Проблемы прочности. 1987. №6. С. 57-63.
4. Волков С.Д., Миронов В.И. Устойчивость процесса усталостного разрушения конструкционных материалов // Уральский полит. ин-т. Свердловск, 1981. 72с.: ил. библ. 12 назв. Деп. в ВИНИТИ 02.12.81. № 5459-81.
5. Лебедев А.А., Чаусов Н.Г. Установка для испытания материалов с построением полностью равновесных диаграмм деформирования // Проблемы прочности. 1981. №12. С. 104-105.
6. Миронов В.И., Микушин В.И., Владимиров А.П. и др. Установка для определения механических свойств материала на стадии разупрочнения // Завод. лаборатория. 2001. Т.67. №3. С. 48-51.
7. Постон Т., Стюарт И. Теория катастроф и ее приложения. М.: Мир, 1980. 608 с.
8. РаботновЮ.Н. Механика деформируемого твердого тела. М.: Наука, 1988. 752 с.
9. Анциферов В.Н., Соколкин Ю.В., Ташкинов А.А. и др. Волокнистые композиционные материалы на основе титана. М.: Наука, 1990. 136 с.
Поступила 21.01.2004 г.
