Определение напряженного состояния и параметров разрушения тонкостенных клееных и клееклепаных элементов авиационных конструкций с трещинами тема диссертации и автореферата по ВАК РФ 05.07.03, кандидат технических наук Тягний, Анатолий Владимирович

• Тягний, Анатолий Владимирович
• кандидат технических наук
• 2004

Оглавление диссертации кандидат технических наук Тягний, Анатолий Владимирович

ВВЕДЕНИЕ.

1. РАСЧЕТ ОСТАТОЧНОЙ ПРОЧНОСТИ И ДОЛГОВЕЧНОСТИ

КЛЕЕНЫХ И КЛЕЕКЛЕПАНЫХ ДВУХСЛОЙНЫХ ПАНЕЛЕЙ С

ТРЕЩИНАМИ МЕТОДОМ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ.

1.1. Основные соотношения теории упругости и основные решения анизотропных пластин.

1.2. Эффективные упругие постоянные пластин из слоистого композитного материала.

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

1.4. Моделирование комбинированного клееклепаного соединения.

1.5. Коэффициенты интенсивности напряжений.

1.6. Дискретизация областей и численная реализация уравнений.

1.7. Влияние изгиба клееклепаной панели на коэффициент интенсивности напряжений.

1.8. Деформирование двухслойной панели с трещиной.

1.9. Численные и экспериментальные исследования.

1.1 (^Прогнозирование развития усталостных трещин.

2. РАСЧЕТНО-ЭКСПЕРИМЕНТАЛЬНЫЙ МЕТОД ОПРЕДЕЛЕНИЯ

ПАРАМЕТРОВ РАЗРУШЕНИЯ КОНСТРУКЦИИ С ТРЕЩИНОЙ ПО

РАСКРЫТИЮ ЕЕ БЕРЕГОВ.

2.1. Метод расчетной оценки коэффициентов интенсивности напряжений.

2.2. Определение напряжений на месте трещины, действовавших до ее появления.

2.3. Численная реализация при конечном числе точек измерений раскрытия трещины.

2.4. Определение положения вершин трещины.

2.5. Численные исследования.

2.6. Экспериментальные исследования.

2.7. Прогнозирование длительности роста трещины с использованием данных о раскрытии ее берегов.

Введение диссертации (часть автореферата) на тему «Определение напряженного состояния и параметров разрушения тонкостенных клееных и клееклепаных элементов авиационных конструкций с трещинами»

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

Развитие методов линейной механики разрушения и их практическое использование при проектировании создали предпосылки появления надежных, безопасных и экономичных конструкций летательных аппаратов (JIA). Стремление использовать экономическую отдачу каждого самолета до его полного изнашивания, т.е. эксплуатировать до появления трещин в силовых элементах, и при этом гарантировать безопасность полетов, привело к новому подходу в определении срока службы конструкции - принципу безопасных повреждений [6, 146, 149]. Этот принцип основан на предположении, что во время эксплуатации в конструкции присутствуют трещины, размер которых меньше или равен минимально обнаруживаемому средствами неразрушающего контроля. Задача конструктора - на стадии проектирования JIA предусмотреть, чтобы повреждение, существующее в элементе конструкции, независимо от его характера (случайное, коррозионное, усталостное и т.д.) не привело к катастрофическому разрушению от действия однократной высокой нагрузки (остаточная прочность) или низких, но часто повторяющихся нагрузок (остаточная долговечность). При этом сочетание свойств материалов, конструктивных особенностей и уровней допускаемых напряжений должно обеспечить достаточно медленный рост и достаточно большой предельный размер повреждений, а периодичность, методы и качество осмотров должны способствовать надежному обнаружению повреждения, прежде чем оно достигнет опасного размера и распространится на жизненно важные элементы конструкции. Этот принцип позволяет гарантировать безопасность даже при появлении ранних случайных трещин неусталостного происхождения.

Задача определения остаточной прочности и долговечности клееных и кле-еклепаных элементов конструкций возникла в связи с широким использованием клеевых и клеемеханических1 соединений в авиакосмической технике. При производстве ДА клеи используются для соединения металлических и неметаллических частей конструкций - обшивки со стрингерами и другими элементами каркаса крыла, фюзеляжа, хвостового оперения, в конструкциях герметических кабин, составных клееных лонжеронов, других силовых конструкциях [19, 25, 58, 61, 67, 68, 101, 145]. Для предотвращения хрупкого разрушения конструкций в период эксплуатации или для его локализации, а также для ремонта поврежденных элементов в авиационной технике в последние годы активно используются приклеенные или клееклепаные подкрепляющие элементы - узкие стопоры трещин и широкие ремонтные накладки, в том числе и из КМ [108, 109, 111,128].

Повышенный интерес к клеевым соединениям объясняется рядом преимуществ, которыми они обладают по сравнению с заклепочными и сварными. При замене заклепочных соединений клеевыми значительно снижается вес конструкции (по некоторым оценкам до 15% [85]), что дает немалый экономический эффект, особенно для больших самолетов. Склеивание снижает стоимость производства (до 25% [99]), улучшает аэродинамические характеристики поверхностей, повышает технические свойства конструкций, такие, как усталостная прочность [61], стойкость к коррозионному растрескиванию, что в свою оче

1 Под клеемеханическими соединениями будем понимать клееклепаные, клее-болтовые, клеесварные соединения. В дальнейшем часто будем пользоваться термином «клееклепаные соединения» для обозначения всех указанных типов клеемеханических соединений. редь снижает расходы по техническому обслуживанию и увеличивает срок эксплуатации.

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

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

При создании JIA и доработке его конструкции важную роль играют натурные испытания, одна из основных целей которых - выявление трещино-опасных зон. В процессе испытаний зоны с повреждениями требуют особого внимания исследователя, контроля за распространением трещин. Оперативный контроль (измерение) параметров, определяющих степень опасности состояния поврежденной конструкции при нагружении, дает возможность прогнозировать поведение испытываемого объекта и принимать обоснованное решение о своевременной остановке испытаний и о дальнейшем ремонте или доработке конструкции. Прогнозирование остаточной прочности и роста трещины основано на понятии коэффициента интенсивности напряжений (КИН) в вершине трещины. Для его достоверного определения на основе данных, измеренных во время проведения эксперимента (испытаний), необходима разработка соответствующих методов.

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

Анализ распределения напряжений и деформаций в составном элементе конструкции с повреждением является начальным этапом для расчета его остаточной прочности и долговечности. Поэтому разработка эффективных расчетных и расчетно-экспериментальных методов определения НДС плоских многослойных клееных и клееклепаных элементов из изотропных и анизотропных материалов — весьма актуальная проблема как для задач проектирования, так при испытаниях конструкций JIA.

Разрушение композитного материала - одна из наиболее сложных областей механики деформируемого твердого тела. Применение методов линейной механики разрушения и теории упругости к этим материалам усложнено прежде всего из-за анизотропии и неоднородности структуры КМ. Однако существует ряд подходов, в которых разрушение КМ с некоторыми допущениями удается согласовать с подходами линейной механики разрушения. В работах [139, 140, 153] представлен критерий разрушения для трещин, базирующийся на коэффициенте плотности энергии деформации, который определяется на основе КИН. С. Поу [130] предложил деформационный критерий разрушения волокон перед вершиной трещины в слоях, несущих основную нагрузку. Использование этого критерия позволяет определять критический КИН для трещины в композитном материале любой толщины и укладки по характеристикам разрушения для однонаправленного материала, за исключением случаев расслоения и расщепления материала. В работе [51] развит предложенный в [47, 49] силовой критерий для определения критической нагрузки и направления развития трещины в ор-тотропном материале с использованием известного асимптотического поведения упругих напряжений около вершины трещины.

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

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

Среди возможных моделей клеевого соединения простейшей являются жестко соединенные упругие мембраны [95, 124]. В этой модели предполагают, что толщина клеевого слоя равна нулю, а распределением напряжений по толщине слоев можно пренебречь. Неучет толщины и податливости клеевого слоя приводит к незначительным контактным напряжениям внутри области склеивания и высокой концентрации касательных напряжений вдоль границы [124], а значит и к неоправданно заниженной расчетной предельной нагрузке.

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

В большинстве работ, посвященных клеевым соединениям, используются модели, занимающие промежуточное положение между двумя указанными моделями и отличающиеся различной степенью упрощений и приближений. Одна из таких моделей, заключающаяся в представлении элементов клеевого соединения в виде мембран, соединенных клеевым слоем конечной толщины, работающим только на сдвиг (см., например, [116]), получила широкое распространение. Ряд приложений этой модели к исследованию тонких подкрепляющих элементов, приклеенных к упругой пластине с трещиной, двухслойных склеенных панелей с трещиной в одном из слоев, а также ступенчатых и клиновидных соединений появились впервые в работах [79, 94, 96].

В работах К.Арина [77, 78, 79], Т. Свифта [143, 144] на основе решения задачи о действии сосредоточенной силы в бесконечной изотропной и ортотроп-ной пластине с прямолинейным разрезом и с использованием модели контакта по линии проведены обширные расчетные и расчетно-экспериментальные исследования для пластин с трещинами и приклеенными узкими ребрами. Позднее этот подход развивался в работе [87] для изотропных материалов в предположении, что касательные усилия по толщине клееной конструкции изменяются линейно, а в клеевом слое постоянны.

В работах В.Н. Максименко, В.Н. Павшока [32-35] методом интегральных уравнений решены задачи о деформировании бесконечной (полубесконечной, конечной) пластины с криволинейными трещинами, отверстиями и произвольно расположенными приклеенными и клееклепаными стопорами трещин. Исследовано влияние на остаточную прочность поврежденной панели: жесткост-ных параметров элементов, размеров и взаимного расположения трещины, подкрепления, отверстия; упруго-пластического поведения материалов клея, заклепок, подкреплений; частичного повреждения клеевого слоя, заклепок, разрыва подкрепляющего элемента. Ряд результатов подтвержден экспериментально [23].

Что касается анализа НДС кпеемеханических соединений, количество работ, посвященных им, относительно невелико [2, 7, 15, 17, 34, 44, 102]. Отметим работу [102], в которой задача расчета НДС ьслееболтового ступенчатого соединения сводится к решению линейного дифференциального уравнения четвертого порядка, и работы [2, 15], в которых методом конечного элемента исследуется НДС нахлесточных клеемеханических соединений, а также конструкций типа пластины с трещиной, подкрепленной ремонтной накладкой. Анализу разрушения и оценке долговечности клееных и клееклепаных соединений посвящены работы [ 55, 56, 69, 70,114].

Для торможения образовавшихся в тонкостенной конструкции трещин наряду с одномерными подкрепляющими элементами используются широкие приклеенные или клееклепаные накладки из высокопрочных конструкционных материалов, которые, как правило, устанавливаются во время ремонта конструкции [80, 81, 86-89, 92, 108-112, 119, 122, 125, 126, 132, 133, 135, 136, 151, 152]. Такие ремонтные мероприятия представляют собой эффективное средство блокирования трещин и усиления листовых элементов конструкций с концентраторами напряжений. Обширные исследования в этом направлении были проведены Р.Джонсом и Р. Кэллинаном [108-111], которые использовали для расчетов метод конечных элементов и подход, развитый ранее Р. Митчеллом, Р. Вули, Д. Чивирутом [122], позволяющий учитывать линейное изменение касательных напряжений по толщине клееной конструкции. Аналогичные задачи методом конечного элемента для случая кругового выреза или прямолинейной трещины, металлической или композитной накладки, клеевого или клеемеха-нического способа присоединения решали Н.С. Галкина, В.И. Гришин и Т.К. Бегеев [12, 15], Д. Андерсон, С. Джу, У. Макги [75], М. Ратвани [132, 133]. Р. Чандра с соавторами [86-89] рассмотрел проблему численного и экспериментального (фотоупругость) определения КИН в вершинах трещины, подкрепленной накладкой, а также исследовал закономерности роста трещины при нагружении, имитирующем спектр полетных нагрузок. Аналитическим методом Д. Картрайт, А. Юнг, Д. Рук, Г. Доурик [92, 151, 152] определяли КИН в вершинах прямолинейной трещины в изотропной пластине, подкрепленной эллиптической или прямоугольной (представляемой в виде набора одномерных элементов жесткости) изотропной накладкой.

Широкое применение в авиационной технике находят многослойные клееные элементы конструкций. Это объясняется более высокой вязкостью разрушения по сравнению с монолитными элементами той же толщины [106, 107, 132]. В полной мере это относится к появившимся в 80-е годы материалам типа AJlOP (металлоорганопластики - многослойные алюминиевые пластины, армированным органическими арамидными волокнами) [100, 120, 138, 150], к появившимся позже материалам типа СИАЛ (Стеклопластик И Алюминий) [57, 62], которые обладают рядом уникальных свойств. Для них применимы технологические процессы обработки традиционных металлических листовых материалов, но в то же время, как и любые другие КМ, они обладают очень хорошими статическими и усталостными характеристиками, особенно в направлении волокон. Сильное замедление или прекращение роста трещин объясняется тем, что позади вершины трещины, движущейся в тонких алюминиевых листах, армирующие волокна остаются целыми и сдерживают раскрытие трещины, а также снимают часть нагрузки с листов.

Начало исследования многослойных клееных конструкций аналитическим методом было положено в работе Ф. Эрдогана и К. Арина [94]. Авторы рассматривали слои как мембраны, соединенные склеивающим слоем конечной толщины, работающим только на сдвиг. В работе дано аналитическое решение задачи для случая двухслойной панели с одной прямолинейной трещиной в металлическом слое и неповрежденным композитным слоем. Позднее этот подход был развит в работах [76, 84, 98, 103, 104, 113, 115, 132, 133]. М. Джесит и Ф. Эрдоган [98] рассмотрели влияние толщины и упругих свойств клеевых слоев на разрушение многослойных конструкций, ввели две модели клеевого слоя - приближенную в виде комбинации пружин, работающих на растяжение и сдвиг, и точную, континуальную, не использующую какие-либо упрощающие предположения. М. Ратвани [132, 133] предложил метод учета влияния изгибной жесткости слоев в зоне повреждения, провел параметрические исследования роста усталостной трещины в двухслойной клееной металлической панели. В.Ко [115] обобщил метод решения на случай ортотропных слоев при сложном нагружении, а С. Хонг, X. Ро и К. Джеонг [103, 104] исследовали влияние направления анизотропии и типа нагружения на КИН в вершинах трещины. К. Бигелоу [84] решила задачу для случая ортотропной пластины с прямолинейной трещиной, подкрепленной ортотропной полуплоскостью в предположении упруго-пластического деформирования клеевого слоя. В. Энг привел решение задачи о трещине в среднем слое трехслойной анизотропной панели [76].

На основании теоретических и экспериментальных исследований процесса роста усталостных трещин, проведенных, в основном, на плоских образцах с надрезом и инициирующей трещиной достаточно большой длины, установлено, что определяющими параметром скорости роста трещины при циклическом нагружении в большинстве случаев является соответствующее циклическое изменение коэффициентов интенсивности напряжений, в общем случае плоского напряженного состояния - К\ и Ki [46, 50, 64].

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

Для определения статистических характеристик значений случайных нагрузок используют кривые их повторяемости, которые показывают вероятное число повторений нагрузки, равной или большей данной, за час наработки (полет, год эксплуатации и т.д.). Такие кривые получают в результате соответствующих методик обработки действительных процессов измерения во времени нагрузок, замеренных на конструкциях [11, 20].

При известном законе изменения во времени эксплуатационных нагрузок необходимо использовать какой-либо метод подсчета циклов. К ним относятся методы максимумов, экстремумов, размахов, полных циклов с различными видоизменениями. Наиболее широко для подсчета числа циклов применяются метод полных циклов и метод дождя [11, 20].

Результирующие кривые повторяемости нагрузок, методы схематизации нерегулярных спектров нагружения не дают никакой информации о последовательности действия нагрузок, поэтому с их использованием при расчете длительности роста усталостных трещин возникают определенные трудности, связанные с взаимным влиянием нагрузок. Это влияние проявляется, в основном, в торможении роста трещины после действия высокой растягивающей нагрузки из-за остаточных напряжений сжатия в пластической зоне у вершины трещины, а также ослаблением эффекта торможения при действии сжимающих внешних нагрузок после больших растягивающих. Для учета эффектов торможения предложено много расчетных схем [6] (модель замедления Уилера, обобщенная модель замедления Уиленборга и др.).

Распространение усталостной трещины в элементах конструкций под действием нерегулярных нагрузок прогнозируется интегрированием зависимостей, связывающих скорость роста трещины и параметры цикла изменения КИН (размах КИН, коэффициент асимметрии), которые представлены кинетической диаграммой усталостного разрушения материала. Предложено более ста способов аналитического представления (аппроксимации) кинетических диаграмм, среди них наиболее известны формулы П. Пэриса, Р. Формана, Г.П. Черепанова и др. [6, 46, 64]. Они отличаются способом построения (теоретические, эмпирические), количеством учитываемых физико-механических констант материала и вводимых эмпирических параметров и имеют различные области применимости. Интегрирование производится либо по каждому циклу нагрузки с использованием различных моделей замедления, либо приведением внешних нагрузок к эквивалентному пульсирующему циклу по повреждаемости [52] и последующим интегрированием. Первый путь требует значительных вычислительных затрат, разработки высокоэффективных алгоритмов и программ расчета КИН, но более достоверно описывает процесс разрушения. Второй путь дает существенно большие погрешности в прогнозировании роста трещины, но значительно снижает трудоемкость вычислительных работ.

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

Экспериментальные методы определения КИН охватывают весьма большой и разнообразный круг вопросов. К ним относятся методики и средства измерения, основанные на различных физических принципах, методы расчета КИН, связанные с различными математическими моделями представления полей НДС в окрестности трещины, а такие анализ и обобщение результатов экспериментальных исследований.

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

Основными методами экспериментального исследования НДС являются: тензометрический, поляризационно-оптический, рентгенографический, методы хрупких покрытий, делительных сеток, метод голографического муара.

Для экспериментального определения КИН можно использовать любую методику, позволяющую измерить деформации или перемещения. Наиболее простой способ заключается в том, что по замеренным деформациям или перемещениям в окрестности вершины трещины, используя асимптотическое разложение НДС около вершины трещины, строят график зависимости КИН в зависимости от расстояния до вершины трещины и экстраполируют его до точки г — 0. Очевидно, что результаты зависят от способа экстраполяции. С использованием данной методики в работе [142] методом фотоупругости определен КИН в образце с краевой трещиной. В работе [83] получены весьма точные результаты определения КИН в панели с центральной трещиной, измеренного с помощью тензорезисторов. В работе [18] для измерения деформаций использовалась прецизионная координатная сетка с шагом 0,1 мм и специальные цепочки миниатюрных фольговых тензорезисторов с базой 0,8 мм. Переход от измеренных деформаций к напряжениям в пластической области осуществляется в рамках деформационной теории Ильюшина с учетом переменного коэффициента поперечной деформации.

Джеймсом и Андерсеном [105] впервые была применена иная экспериментальная методика оценки КИН. В этой методике используется связь скорости распространения усталостной трещины с КИН aN

Функцию /(Л/Г) можно определить в испытании на распространение усталостной трещины в образце, для которого величина АК известна. Определив скорость распространения трещины в образце или конструкции со сложной геометрией, с помощью приведенного уравнения можно определить КИН.

Некоторые исследователи используют многопараметрическое представление полей деформаций и напряжений [16, 82, 91, 121 ], или перемещений [14] в окрестности вершины трещины в анизотропной пластине. Далее неизвестные параметры определяются методом наименьших квадратов по нескольким измерениям деформаций или перемещений в произвольных точках около вершины трещины.

С.В. Шкараев [72-74], используя принцип суперпозиции, приводит исходную задачу определения КИН для трещины в нагруженной пластине произвольной конфигурации к упрощенной схеме: ненагруженная бесконечная или полубесконечная изотропная пластина с внутренним или краевым разрезом, на берега которого действуют компенсирующие усилия. Неизвестные параметры аппроксимации функции компенсирующих усилий восстанавливались по данным скачков смещений берегов трещины с использованием известных фундаментальных решений для указанных областей.

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

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

Различные расчетно-экспериментальные методы используются для определения номинальных напряжений [8, 66, 134, 137, 147, 148] (вырезание слоев, сверление отверстий, нанесение надрезов) и КИН по полям деформаций и перемещений в области перед вершиной трещины [8, 90,93,97,121, 123,127,129, 147, 148] (тензометрия, регистрация раскрытия вблизи вершины трещины с помощью датчиков или компьютерной обработки визуальных изображений, лазерные спекл-методы, фотоупругие и голографические методы, пьезоэлектрические, волоконно-оптические датчики и др.).

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

Представляется актуальным разработать метод расчетно-эксперимен-тального определения КИН в тонкостенных элементах конструкций JIA для общего случая криволинейной трещины в анизотропной пластине при наличии концентраторов напряжений (отверстий) и подкрепляющих элементов (приклеенных, клееклепаных стопоров, накладок и т.п.).

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

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

- метод расчета НДС, остаточной прочности, остаточной долговечности металлических и композитных двухслойных клееных и клееклепаных панелей конструкций ДА с трещинами, однослойных поврежденных панелей, с приклеенными и клееклепаными подкрепляющими элементами в виде широких ремонтных накладок;

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

- результаты параметрических исследований ряда новых задач для составных конструктивных элементов с трещинами и клеемеханическим креплением;

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

- общие аналитические решения, позволяющие на основе данных о распределении разности перемещений (скачков смещений) берегов вдоль линии трещины получать КИН и номинальные напряжения на месте трещины, действовавшие до ее появления;

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

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

В работе отражены исследования автора, выполненные в СибНИА им. С.А. Чаплыгина и НГТУ.

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

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

Предложенные методики и пакеты программ дают возможность проводить анализ остаточной прочности и остаточной долговечности элементов авиационных конструкций, устанавливать величину допустимых повреждений, критических нагрузок, обосновывать пути повышения живучести конструкций, выбирать оптимальные конструктивные параметры на стадии проектирования с учетом принципа безопасных повреждений. Разработанные методики и пакеты прикладных программ внедрены в практику СибНИА им. С.А. Чаплыгина.

Апробация работы. Основные положения диссертационной работы докладывались на Всесоюзном симпозиуме «Проблемы автоматизации в прочностном эксперименте», (Новосибирск, 1988), на Всесоюзном симпозиуме по механике разрушения (Житомир, 1990), на Российско-Китайской научной конференции по проблемам прочности авиаконструкций (Новосибирск, 1995), на Всероссийской школе-семинаре по современным проблемам механики деформируемого твердого тела (Новосибирск, 2003), на международном симпозиуме КОРУС (Томск, 2004).

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

Основные результаты работы изложены в 11 научных публикациях [36, 38, 39,41-43, 45 , 59, 60, 117, 118].

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

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

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

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

Достоверность предлагаемых моделей клеевых, клеемеханических соединений и расчетной методики иллюстрируется сравнением с некоторыми известными результатами, полученными методом конечных элементов или экспериментально.

Решен ряд задач по оценке характеристик живучести клееных и клееклепаных панелей. Предлагаются способы учета:

- изгибного эффекта панели;

- отслоения в клеевом слое;

- частичного разрушения подкрепляющего элемента;

- симметрии задачи.

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

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

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

В заключении работы сформулированы основные научные результаты.

В приложении представлен документ об использовании результатов работы.

Заключение диссертации по теме «Прочность и тепловые режимы летательных аппаратов», Тягний, Анатолий Владимирович

ЗАКЛЮЧЕНИЕ

Основные научные результаты работы заключаются в следующем:

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

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

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

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

5. Разработан расчетно-экспериментальный метод определения КИН и номинальных напряжений, метод определения положения вершин трещины по экспериментально полученным значениям скачков смещений берегов трещины в конечном числе дискретных точек.

6. Созданные расчетный и расчетно-экспериментальный методы позволяют:

- устанавливать критические и допускаемые размеры повреждений,

- делать прогнозы развития трещин и несущей способности поврежденных конструкций,

- выбирать рациональные конструктивные параметры на стадии проектирования конструкции с учетом повышенной живучести,

- осуществлять оперативное определение критических нагрузок и запаса остаточной прочности при статическом нагружении, а также их изменение при периодическом нагружении,

- устанавливать периодичность осмотров элементов авиационных конструкций во время эксплуатации и при проведении испытаний.

Список литературы диссертационного исследования кандидат технических наук Тягний, Анатолий Владимирович, 2004 год

1. Белоус А.А., Хватан A.M. Расчет клеемеханического соединения внахлестку // Проектирование, расчет и испытания конструкций из композитных материалов.- 1979.- Вып.7- С. 49-55.

2. Белоцерковский С.М., Лифанов И.К. Численные методы в сингулярных интегральных уравнениях. М.: Наука, 1985.

3. Березин И.С., Жидков Н.П. Методы вычислений. Т. 1. М.: Наука, 1966.

4. Берт Ч. Расчет пластин // Композиционные материалы. Т. 5. Анализ и проектирование конструкций / Ред. Л. Браутман и Р. Крох М.: Машиностроение, 1978.-С. 154-209.

5. Броек Д. Основы механики разрушения.- М.: Высшая школа, 1980- 368 с.

6. Брушковский А.Л., Брушковская И.Д. Унификация приближенных методов расчета заклепочных и клееклепаных соединений // Вопр. мех. деформируем. тверд.тела 1983-№4 -С. 5-10.

7. Вайншток В.А., КравецП.Я. Расчет КИН и номинальных напряжений в плоскости трещины по раскрытию ее берегов // Пробл. прочности.- 1990-№6.-С. 27-34.

8. Вапник В.Н. Алгоритмы и программы восстановления зависимостей. М.: Наука, 1984.

9. Васильев Ф.П. Численные методы решения экстремальных задач. М.: Наука, 1980,-520 с.

10. Воробьев А.З., Олькин Б.И., Стебенев В.Н., Родченко Т.С. Сопротивление усталости элементов конструкций. М.: Машиностроение, 1990. - 240 с.

11. Галкина Н.С., Гришин В.И. О концентрации напряжений около отверстий в пластинах с приклеенными накладками // Учен, записки ЦАГИ.-1978 Т. 9, № 1.-С. 137-141.

12. Гахов Ф.Д. Краевые задачи. М.: Физматгиз, 1963 -582с.

13. Городниченко В.И., Дементьев А.Д. Определение коэффициентов интенсивности напряжений в вершине сквозной трещины по полям перемещений //Учен. зап. ЦАГИ.- 1988.- Т. 19, № 6.- С. 82-93.

14. Гришин В.И., Бегеев Т.К. Коэффициенты интенсивности напряжений в пластине с центральной поперечной трещиной, усиленной накладками из композитного материала // Механ. композ. матер 1986 — № 4 - С. 696-700.

15. Дементьев А.Д. Расчет коэффициентов интенсивности напряжении в вершине сквозной трещины по данным тензометрии //Учен. зап. ЦАГИ-1987.- Т. 18, № 5.- С. 83-88.

16. Заривняк И.С., Заривняк Г.Р. Расчет клееклепаных и клеесварных соединений с учетом концентрации напряжений и начальных неправильностей // Вопр. проектир. самолет, конструкций 1983.- № 4.- С. 135-138.

17. Злочевский А.Б., Бондарович А.А., Шувалов А.Н. Определение коэффициента интенсивности напряжений тензометрическими методами //Проблемы прочности 1979 - № 6 - С. 44-47.

18. Кейгл Ч. Клеевые соединения.- М.: Мир, 1971,- 295 с.

19. Когаев В.П. Расчеты на прочность при напряжениях переменных во времени. М.: Машиностроение, 1977. - 324 с.

20. Когаев В.П., Махутов Н.А., Гусенков А.П. Расчеты деталей машин и конструкций на прочность и долговечность: Справочник. М.: Машиностроение, 1985.-488 с.

21. Космодамианский А.С. Напряженное состояние анизотропных сред с отверстиями или полостями. Донецк, Киев: Вища школа, 1976 - 200 с.

22. Краснов Л.А., Макснменко В.Н., Павшок В.Н., Тырнн В.П. Определение коэффициентов интенсивности напряжений для трещин в панелях, подкрепленных приклеенными накладками // Пробл. прочности 1989 - №5.- С.73-76.

23. Крылов В.И. Приближенное вычисление интегралов. М.: Физматгиз — 1959. - 328 с.

24. Крысин В.Н. Слоистые клееные конструкции в самолетостроении- М.: Машиностроение.- 1980.- 228 с.

25. Левин М.Г. Программное обеспечение для решения задач численного анализа на СМ ЭВМ. Кишинев: Щтиинца, 1991. - 240 с.

26. Лехницкий С.Г. Анизотропные пластинки. М.: ГИТТЛ, 1957 - 464 с.

27. Максименко В.Н. Задача о трещине в анизотропной полуплоскости, подкрепленной упругими накладками // Динамика сплошной среды: Сб. науч. тр. / АН СССР, Сиб. отд-ние. Ин-т гидродинамики. 1990. Вып. 99. С. 41-46.

28. Максименко В.Н. К контактной задаче для анизотропной пластины, подкрепленной ребром жесткости // Изв. АН СССР, МТТ.- 1981.-№1- С. 159-165.

29. Максименко В.Н. Павшок В.Н. К расчету анизотропной пластины с трещинами, усиленной клееклепаными ребрами // Журнал прикл. механ. и техн. физ. 1992. - № 1. - С. 133-140.

30. Максименко В.Н. Предельное равновесие анизотропной пластины, ослабленной эллиптическим отверстием и системой трещин сложной формы // Учен. зап. ЦАГИ. 1987. - Т. 18, № 3. - С. 24-29.

31. Максименко В.Н., Павшок В.Н. Влияние приклеенных накладок на интенсивность напряжений в вершинах трещин в анизотропной пластине // Прикл. механ 1989.- Т.25, №5.- С. 69-75.

32. Максименко В.Н., Павшок В.Н. Влияние приклеенных ребер жесткости на характеристики остаточной прочности панелей из анизотропных и изотропных материалов // Учен, записки ЦАГИ 1987- Т. 18, №2 - С. 84-92.

33. Максименко B.H., Павшок В.Н. К расчету анизотропной пластины с трещинами, усиленной клееклепаными ребрами // Журнал прикл. механ. и техн. физ.- 1992,- №1.- С. 133-140.

34. Максименко В.Н., Павшок В.Н. Расчет подкрепленной пластины с трещиной в случае нелинейной работы накладок и склеивающего слоя // Прочность и аэроупругость авиац. конструкций.- Казань: КАИ, 1988 С. 27-33.

35. Максименко В.Н., Судаков Н.С. Взаимодействие приклеенного ребра жесткости с пластиной из композитного материала // Механ. композ. матер-1983.-№3.-С. 493-498.

36. Максименко В.Н., Тягний А.В. Расчет напряженного состояния клееклепаных слоистых пластин с трещиной // Ученые записки ЦАГИ. 1990.- Т.21, №5.-С. 92-101.

37. Максименко В.Н. Влияние приклепанных ребер жесткости на развитие трещин возле отверстия Журнал прикл. механ. и техн. физ. 1988, № 2-С. 133-140.

38. Максименко В.Н., Тягний А.В. Расчетно-экспериментальный метод определения параметров разрушения конструкций с трещинами //ПМТФ — 2004. Т. 45, №4.-С. 168-175.

39. Матвиенко В.А. Напряженно-деформированное состояние клееклепаных соединений // Авиационная промышленность 1987.- №4 — С. 5-7.

40. Механика разрушения и прочность материалов: Справочное пособие: В 4 т. / Под общей редакцией Панасюка В.В. Киев: Наукова думка, 1988.

41. Морозов Н.Ф., ПаукпггоМ.В. Дискретные и гибридные модели механики разрушения. СПб.: Изд-во СПбГУ, 1995. 157 с.

42. Мусхелишвили Н.И. Сингулярные интегральные уравнения. М.: Физматгиз, 1968.-512 с.

43. Новожилов В.В. К основам теории равновесных трещин в упругих телах // ПММ. 1969. - Т. 33. Вып. 5. - С. 797-812.

44. Партон В.З., Морозов Е.М. Механика упругопластического разрушения. -М.: Наука, 1985. 502 с.

45. Петров Ю.В., Поникаров Н.В. О направлении роста трещины в ортотропном материале // МТТ. 1998. - № 4. - С. 180-184.

46. Прочность конструкций при малоцикловом нагружении / Под ред. Н.А. Махутова, А.Н. Романова. М.: Наука, 1983. - 272 с.

47. Ратвани М.М. Исследование напряжений в клееных слоистых конструкциях, ослабленных трещинами // Ракетная техника и космонавтика. 1979. — Т. 17, №9.-С. 77-85.

48. Савин Г.Н. Распределение напряжений около отверстий. Киев: Наукова думка, 1968 - 888 с.

49. Семин М.И. Расчетные модели по анализу разрушения и оценке долговечности клеемеханических соединений // Динам, и мех. поврежд. авиац. конструкций.- М 1982.- С. 19-23.

50. Семин М.И., Турчинский JI.A. Методика расчета долговечности клеевых соединений при динамическом нагружении // Пробл. прочности.- 1985-№9.-С. 92-97.

51. Тризно М.С., Москалев Е.В. Клеи и склеивание Л.: Химия, 1980.-174с.

52. Тягний А.В. Метод определения параметров разрушения конструкций с трещинами // Сб. научн. трудов НГТУ 2002. - № 2(28).- С. 89-94.

53. Фрейдин А.С. Прочность и долговечность клеевых соединений. -М.: Химия, 1981.-270 с.

54. Химмельблау Д. Прикладное нелинейное программирование. Пер. с англ. яз. М.: Мир,1975 536 с.

55. Черепанов Г.П. Механика хрупкого разрушения- М.: Наука, 1974.- 640 с.

56. Черепанов Г.П. Механика разрушения композиционных материалов. М.: Наука, 1983.-296 с.

57. Чжен В., Финни И. Метод измерения осесимметричных продольных остаточных напряжений в тонкостенных цилиндрах, сваренных кольцевым швом//Теорет. основыинж.расчетов 1985-№3 -С. 52-57.

58. Шавырин В.Н., Андреев Н.Х., Ицкович А.А. Клеемеханические соединения в технике-М.: Машиностроение, 1968.-231 с.

59. Шавырин В.Н., Рязанцев В.И. Клеесварные конструкции- М.: Машиностроение, 1981.-168 с.

60. Швечков Е.И. Закономерности разрушения и циклическая долговечность клееклепаных соединений // Машиноведение 1984.- №2 - С. 71-76.

61. Швечков Е.И. Методика определения долговечности клееклепаных соединений //Учен, записки ЦАГИ.- 1983.- Т.14, №4.- С. 137- 141.

62. Шерман Д.И. К решению плоской задачи теории упругости для анизотропной среды // Прикл. матем. и механ.- 1942.- Т.6, Вып.6 С. 509-514.

63. Шкараев С.В. Метод определения параметров разрушения элементов конструкций с краевыми трещинами // Физико-химическая механика материалов. -1992.-№6.-С. 29-35.

64. Шкараев С.В. Теоретико-экспериментальный метод определения коэффициентов интенсивности напряжений // Физико-химическая механика материалов. 1989. - № 4. - С. 97-101.

65. Шкараев С.В. Метод определения коэффициентов интенсивности напряжений в подкрепленных панелях // Физико-химическая механика материалов. 1993. -№ 1. - С. 84-89.

66. Andersson J.M., Chu C.S., McGee W.M. Growth characteristics of a fatigue crack approaching and growing beneath an adhesively bonded doubler // Trans. ASME, ser.D.- 1978.-V.100, №1.-P. 52-56.

67. Ang W.T. An anisotropic layered material with a crack // Acta mech.- 1988 — V.72, №3-4.- P. 297-308.

68. Arm K. A note on effect of lateral bending stiffness of stringers attached to a plate with a crack // Eng. Fract. Mech 1975 - V.7, №1- P. 173-179.

69. Arm K. Several intact or broken stringers attached to an orthotropic sheet with a crack// Eng. Fract. Mech.- 1979.- V.ll, №1.-P. 1-8.

70. Arin K. A plate with a crack, stiffened by a partially debonded stringer // Eng. Fract. Mech.- 1974.-V.6, №1.-P. 133-140.

71. Baker A.A. A summary of work on application of advanced fibre composites at ARL Australia // Composites.- 1978 V.9, №1P. 11-16.

72. Barut A., Hanauska J., Madenci E., Ambur D.R. Analysis method for bonded patch repair of a skin with a cutout // Composite Structures 2002 - V.55, №3.-P. 277-294.

73. Berger J.R., Dally J.W. An overdeterministic approach for measuring K\ using strain gages // Exp. Mech. 1988 № 6. - P. 142-145.

74. Bhandari S.K. Etude experimentale du facteur d'intensite des containtes au voisinage de la pointe d'une fissure de fatigue centrale dans une tole mince au moyen des mesures extensometriques. These, Ecole Nat. Superieuere de L'Aeronautique, Paris, 1969.

75. Bigelow C.A. Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet // Eng. Fract. Mech 1986.- V. 25, № l.-P. 79-90.

76. Briens G. Les applications aeronautiques du collage // Mater, et techn- 1984.-V.72, №6-7.- P. 239-246.

77. Chandra R., Guruprasad K. Numerical estimation of stress intensity factors in patched cracked plates // Eng. Fract. Mech.- 1987 V. 27, № 5 - P. 559-569.

78. Chandra R., Murthy M.V.V., Ramamurthy T.S., Rao A.K. Analytical estimation of stress intensity factors in patched cracked plates //Eng. Fract. Mech- 1985 — V. 21, № 3 —P. 479—494.

79. Chandra R., Subramanian A. Stress intensity factors in plates with a partially patched central crack // Exper. Mechanics 1989- V. 29, № 1- P. 1-5.

80. Chandra R., Sunder R. Estimation of fatique crack growth in patched cracked panels // Adv. Fract. Res. / Proc. 7th Int. Conf.Fract/Houston, Tex., 1989 Oxford, 1989.- V. 2- P. 989-997.

81. Chiang F.P., Haveesh T.V. Three dimensional crack tip deformation measured by laser speckles //Proc. SEM Conf. Exp. Mech. (Las Vegas, June, 1985- Brook field, 1985-P. 41-44.

82. Dally I.W., Sanford R.I. Strain-Gage Methods for Measuring the Opening-Mode Stress-Intensity Factor, Ki // Exp. Mech. 1987.- 27, №4. - P. 381-388.

83. Dowrick G., Cartwright D.J., Rooke D.P. The effects of repair patches on the stress distribution in a cracked sheet // Numer. Meth. Fract. Mech. / Proc. 2nd Int. Conf. / Swansea, 1980 Swansea, 1980 - P. 763-775.

84. Durig B.R., McNeill S.R., Sutton M.A., Chao Y.Y., Peters W.H. Determination of mixed mode stress intensity factors using digital image correlation // 6-th Congr. Exp. Mech., Portland, 1988: Proc.-London, 1988. V.l-P. 17-20.

85. Erdogan F., Arin K. A sandwich plate with a part-through and a debonding crack // Eng. Fract. Mech.- 1972.-V.4, №3.- P. 449-458.

86. Erdogan F., Civelek M.B. Contact problem for an elastic reinforcement bonded to an elastic plate // Trans. ASME.- 1974.- V.41, №4.- P. 1014-1018.

87. Erdogan F., Ratwani M.M. Stress distribution in bonded joints // J. Compos. Mater.- 1971.- V.5, №6.- P. 278-393.

88. Fujimoto Y., Shintaku E., Pirker G., Gang Liu. Piezoelectric sensor for stress intensity factor measurement of two dimensional cracks // Eng. Fract. Mech.-2003.-V. 70, № 9.-P. 1203-1218.

89. Gecit M.R., Erdogan F. The effect of adhesive layers on the fracture of laminated structures // Trans. ASME, ser.D.- 1978.- V.100, №1.- P. 1-10.

90. Gluing aluminum // Aerospace America 1986 - V.24, №3 — P. 30.

91. Gunnink J.W. Design studies of primary aircraft structures in ARALL laminates // J. Aircraft.- 1988.- V.25, №11.- p. 1023-1032.

92. Hart-Smith L.J. Adhesive bonding of aircraft primary structures //SAE Techn. Pap. Ser.- 1980.-№801209.- 15 p.

93. Hart-Smith L.J. Bonded-bolted composite joints // J. Aircraft- 1985.- V.22, №11.-P. 993-1000.

94. Hong C.S., Jeong K.Y. Stress intensity factors in anisotropic sandwich plate with a part-through crack under mixed mode deformation // Eng. Fract. Mech-1985.-V. 21, № 11.-P. 285-292.

95. Hong C.S., Ro H.S. Adhesive bonded orthotropic structures with a part-through crack // AIAA J.- 1984.- V. 22, № 3.- P. 443-445.

96. James L.A., Anderson W.E. A simple experimental procedure for stress intensity factor calibration //Eng. Fract. Mech. 1969. - V. 1. - P. 565-568.

97. Johnson W.S. Damage tolerance evaluation of adhesively laminated titanium // Trans. ASME, ser.D.- 1983.-V. 105, № 3.-P. 183-187.

98. Johnson W.S., Rister W.C., Spamer T. Spectrum crack growth in adhesively bonded structure // Trans. ASME, ser.D.- 1978 V.100, №1.- P. 57-64.

99. Jones R., Callinan R.J. A design study in crack patching// Fibre Sci. and Tech-nol.— 1981.-V. 14, № 2.-P. 99-111.

100. Jones R., Callinan R.J. Bonded repairs to surface flaws // Adv.Fract.Res./Proc. 6th Int. Conf. Fract./ New Delhi,1984.- Oxford, 1984.- V.5.- P. 3517-3522.

101. Jones R., Callinan RJ. Finite element analysis of patched cracks // J. Struct. Mech.- 1979.- V.7, №2.- P. 107-130.

102. Jones R., Callinan R.J., Aggarwal K.S. Analysis of bonded repairs to damaged fibre composite structures // Eng. Fract. Mech 1983- V.17, №1.-P. 37-46.

103. Kan H.P., Ratwani M.M. Composite patch repair of cracked thick aluminum structures // AIAA J.- 1984.- V.22, №4.- P. 421^27.

104. Kan H.P., Ratwani M.M. Nonlinear behavior effects in cracked metal-to-composite bonded structures // Eng. Fract. Mech.- 1980 V.15, №1- P. 123-130.

105. Keer L.M., Lin C.T., Mura T. Fracture analysis of adhesively bonded sheets // Trans. ASME, ser. E.- 1976.-V.43, №4.-P. 651-656.

106. Ко W.L. An orthotropic sandwich plate containing a part-through crack under mixed mode deformation // Eng. Fract. Mech.- 1978 V. 10, №1.- P. 15-23.

107. Lubkin J.L., Lewis L.C. Adhesive shear flow for an axially-loaded finite stringer bonded to an infinite sheet // Quart. J. Mech. Appl. Math.- 1970 V.33, №4.-P. 521-533.

108. Mall S., Naboulsi S. Modeling of a cracked metallic structure with bonded composite patch using the three layer technique // Composite Structures- 1996-V. 35, №3.-P. 295-308.

109. Marissen R. Flight simulation behaviour of aramid reinforced aluminium laminates (ARALL) // Eng. Fract. Mech.- 1984.- V. 19, № 2.- P. 261-277.

110. Marur P.R. and Tippur H.V. A strain gage method for determination of fracture parameters in bimaterial systems // Eng. Fract. Mech- 1999- V. 64, № 1,-P. 87-104.

111. Mitchell R.A., Wooly R.M., Chivirut D.J. Analysis of composite-reiforced cutouts and cracks // AIAA J.- 1975.- V.13, №6.- P. 744-749.

112. Miura N., Sakai S., Okamura H. A determination of mode I stress intensity factor from photoelastic isochromatic fringe patterns assisted by image processing technique //9th Int. Conf. Exp. Mech., Copenhagen, 1990: Proc- Copenhagen, 1990.- P. 27-34.

113. Muki R., Sternberg E. On the stress analysis of overlapping bonded elastic sheets // J. Solids and Struct.- 1968 V.4, №1.- P. 75-94.

114. Naboulsi S., Mall S. Characterization of fatigue crack growth in aluminium panels with a bonded composite patch // Composite Structures.- 1997.- V.37, № 3-4.-P. 321-334.

115. Naboulsi S., Mall S. Fatigue crack growth analysis of adhesively repaired panel using perfectly and imperfectly composite patches // Theoretical and Applied Fracture Mechanics 1997.- V. 28, № 1.- P. 13-28.

116. Narendran N., Shukla A., Letcher S. Application of fiber optic sensors to fracture mechanics problem // Eng. Fract. Mech.- 1991- V. 38, № 6.- P. 68-76.

117. Norris M. The use of bonded crack-stoppers in pressure cabins // Aeronaut. J-1981.-V. 85, № 841.-P. 57-62.

118. Picon R., Paris F., Caiias J., Marin J. A complete field method for the photoelastic determination of K\ and K\\ in general mixed-mode fracture // Eng. Fract. Mech.- 1995.- V. 51, № 3.- P. 505-511.

119. Рое C.C., Jr. A unifying strain criterion for fracture of fibrous composite lami-nates//Eng. Fract.Mech.-1983.-V. 17,№2.-P. 153-171.

120. Рое C.C., Jr. Tensile strength of composite sheets with unidirectional stringers and crack-like damage //ACEE Compos. Struct. Technol., NASA CP 2321.-1984.-P. 97-111.

121. Ratwani M.M. A parametric study of fatique crack growth in adhesively bonded structures // Trans. ASME, ser.D.- 1978 V.100, №1- P. 46-51.

122. Ratwani M.M. Analysis of cracked adhesive by bonded laminated structures // AIAA J.- 1979.- V. 17, № 9. P. 988-994.

123. Rendler N.J., Vigness I. Hole-drilling strain-gage method of measuring residual stress //Ibid.- 1966.- V. 6, № 12.- P. 577-586.

124. Roderick G.L. Prediction of cyclic growth of cracks and debonds in aluminum sheets reinforced with boron/epoxy // Fibrous Composites in Structural Design-New York, London.- 1980.- P. 467-481.

125. Rose L.R.F. A cracked plate repaired by bonded reinforcements // Int. J. Fracture.- 1982, II.- V. 18, № 2.- P. 135-144.

126. Rubicki E.F., Shadley J.R., Shealy W.S. An consistent-splitting model for experimental residual stress analysis //Exp. Mech- 1983- V. 23, №4 P. 102109.

127. Schijve J. Crack stoppers and ARALL laminates //Eng. Fract. Mech.—1990.— V. 37, №2.-P. 405—421.

128. Sih G.C. Strain energy density factor applied to mixed mode crack problems // Int. J. Fracture.- 1974.- V.10, №3.- P. 135-144.

129. Sih G.C., Chen E.P. Mechanics of fracture. V.6. Cracks in composite materials-Hague/Boston/London: Martinus Nijhoff, 1981- 546 p.

130. Sih G.C., Paris P.G., Irwin G.R. On cracks in rectangular anisotropic bodies//Int. J. Fract. Mech.- 1965.- V.l, №3.- P. 189-203.

131. Smith D.C., Smith C.W. Photoelastic determination of mixed mode stress intensity factors Eng. Fract. Mech. - 1972. - V.4. - P. 357-366.

132. Swift T. Fracture analysis of adhesively bonded cracked panels // Trans. ASME, ser.D- 1978.-V. 100,№ l.-P. 10-15.

133. Swift T. The effects of fastener flexibility and stiffener geometry on the stress intensity in stiffened cracked sheets // Prosp. Fract. Mech-1974 -P. 419-436.

134. Thrall E.W. Primary adhesively bonded structure technology (PABST) // J. Aircraft.- 1977.- V. 14, № 6.- P. 588-594.

135. Toor P.M. On damage tolerance design of fuselage structure-circumferential cracks // Eng. Fract. Mech.- 1987.-V. 26, № 5.- P. 771-782.

136. Torn Т., Houda K., Fujibayashi Т., Hamano T. A method of evaluating crack opening stress distributions and stress intensity factors based on opening displacement along a crack // JASME Intern. J. Ser I. -1990. V. 32, № 2. - P. 47-54.

137. Vainstok V.A., Kravets P.Y. Estimation of stress intensity factors and nominal stresses from discrete COD values //Eng. Fract. Mech- 1991- V. 38, № 4,5 — P. 89-94.

138. Vlieger H. Damage tolerance of stiffened-skin structures: prediction and experimental verification // Fract. Mech./ Proc. 19th Nat. Symp./ San-Antonio, Tex., 1986.-Philadelphia (Pa), 1988.-P. 169-219.

139. Yeh J.R. Fracture mechanics of delamination in ARALL laminates // Eng. Fract. Mech.- 1988.- V. 30, № 6.- P. 827-837.

140. Young A., Cartwright D.J., Rooke D.P. The boundary element method for analysing repair patches on cracked finite sheets // Aeronaut. J- 1988- V. 92, №920.-P. 416-420.

141. Young A., Rooke D.P., Cartwright DJ. Numerical study of balanced patch repairs to cracked sheets // Aeronaut. J.- 1989 V. 93, № 929 - P. 327-334.

142. Zhang S.Q., Jang B.Z., Valaire B.T., Suhling J.C. A new criterion for composite material mixed mode fracture analysis // Eng. Fract. Mech 1989 - V. 34, № 3-P. 749-769.