Manufacturing defects in composite flanges: modelling and analysis of influence on static strength
- Authors: Anoshkin AN1, Zuiko VY.1, Osokin VM1, Tretyakov AA1, Pisarev PV1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 2 (2016)
- Pages: 5-21
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/179
- DOI: https://doi.org/10.15593/perm.mech/2016.2.01
- Cite item
Abstract
This work presents an approximate methodology based on the numerical simulation aimed at assessing how defects affect static strength of polymer composite flanges in load bearing elements used in aeronautical engineering. Various defects may occur in the process of flange forming in the areas of bending, e.g. curvature of the layers, resin pockets, voids, delamination and others. Thus, there is an actual problem of assessing how technological defects affect the strength of this part. The main defects in composite laminates are defined and a review of literature related to strength problems for such structures are presented in this work. The numerical stress-strain analysis of composite flange with main types of defects under force loading condition was carried out with ANSYS software. A two - dimensional axisymmetric finite element model was used. For the development of the structural model a parametric modelling approach was applied, including defect size, configuration and location parameters. The problem was solved in a general statement for an anisotropic elastic body. In order to keep the original part profile, we reduced the overall thickness of the layers in the local area near the defect by the size of the defect, so the overstated value of structural strength was obtained in simulation. The safety factor of the flange was estimated by stress components using the maximum stress criterion. The stresses in the material principal directions and interlaminar shear and normal stresses were determined for analysis. A comparison of the obtained data with the results for a defect-free flange allows estimating the impact of each defect or its combination on the static strength of flange.
Full Text
Введение Проблема расчетной оценки конструкционной прочности многослойных армированных пластиков на этапе проектирования приобрела в последние годы большую актуальность в связи с ростом объемов использования этих материалов в различных областях современной техники. Ее решение связано с целым комплексом задач и отличается большой сложностью, обусловленной, в частности, многообразием механизмов разрушения армированных пластиков, анизотропией материалов, значительным числом варьируемых переменных, а также возможным присутствием дефектов различного характера [1]. Виды дефектов в полимерных композиционных материалах (ПКМ) можно разделить на производственно-технологические и эксплуатационные. Производственно-технологические дефекты, в свою очередь, можно разделить на геометрические (отклонения в углах укладки, толщинах слоев, складки, смоляные карманы), химические (погрешности химического состава) и механические, вызванные погрешностями при механической обработке. Эксплуатационные дефекты могут возникать под воздействием внешних природных факторов (град, удары молнии), человеческого фактора (падение рабочего инструмента при обслуживании техники), а также являться следствием циклически повторяемых нагрузок (характерных, в частности, для авиационной техники), которые приводят к появлению усталостных трещин. В композитных конструкциях выделяют следующие основные типы дефектов: сколы, трещины, расслоение, неравномерная пропитка армирующего материала, посторонние включения, поры, вздутие на поверхности материала, смоляные карманы, усадочные раковины, складки и искривления в слоях, подрезка слоев. Наиболее часто встречаемыми дефектами в конструкции деталей авиационной техники являются смоляные карманы, поры, складки (или волны в одном или нескольких слоях) и расслоения. Последние являются наиболее распространенным видом дефекта и часто считаются определяющим фактором при решении вопроса об использовании композитов в конструкциях. Отслоения могут возникнуть в местах концентрации напряжений, в районе скачкообразного изменения толщины материала в результате несовершенств технологии производства или действия эксплуатационных нагрузок. Одной из особенностей процесса расслоения является локальная потеря устойчивости и выпучивание отслоившегося слоя с последующим ростом зоны дефекта. Это может стать причиной общего разрушения элемента конструкции при нагрузке значительно меньше расчетного уровня. Вопрос влияния различных дефектов на механические характеристики композитных конструкций весьма детально изложен в работах [2, 3]: проведен всесторонний обзор существующих математических и экспериментальных методов, отмечены заслуги отечественных ученых В. В. Болотина, Л.М. Качанова, Ю.М. Тарнопольского, зарубежных исследователей Г. Чея, Ч.Д. Бэбкока, В. Боттего, У. Йина, Р.К. Кананья, Д.Э. Кардоматиса, Г. Симитсеса, С. Салама и др. [4-21]. Тем не менее, несмотря на большое количество исследований, оценка прочности композитов с типичными повреждениями на этапе проектирования является по-прежнему весьма сложной задачей ввиду недостаточного развития теоретической базы, позволяющей обоснованно свести реальный эксплуатационный дефект к его расчетному аналогу [1]. В настоящей работе предложена приближенная методика оценки влияния технологических дефектов на статическую прочность фланцев из ПКМ корпусных деталей авиационной техники, основанная на схематичном описании дефектов и использовании современных методов численного моделирования для расчета конструкции. Проиллюстрировано применение методики для сравнения полей напряжений и оценок прочности дефектного и бездефектного фланца [24]. 1. Постановка задачи и подходы к моделированию дефектов На первом этапе решения задачи об оценке влияния дефектов на прочность композитной детали необходимо оценить характерные размеры, конфигурацию и расположение локальных дефектов по толщине и площади элемента конструкции. В качестве объекта исследования рассматривается типичный отгибной фланец кожуха авиадвигателя, изготовленный из стеклопластика. Фланец является наиболее нагруженным элементом, определяющим прочность кожуха в целом [22]. На рис. 1 показаны основные типы дефектов, возникающих в двух конструктивных вариантах композитного фланца. Так, в зоне перехода цилиндрической части фланца к его отгибному участку происходит образование складок в слоях материала, что приводит к появлению поры или смоляного кармана. Кроме того, возможно образование и продольного смоляного кармана вследствие недостаточного или неравномерно распределенного давления при формовании детали. Отметим, что аналогичные дефекты могут наблюдаться и в других ответственных деталях типа шпангоутов, балок, узлов крепления лопаток и т.п., где имеются подобные участки с перегибами слоев. Дефекты в зоне перегибов слоев, особенно типа складок, плохо диагностируются всеми известными методами неразрушающего контроля (шерография, Рис. 1. Видимые дефекты в конструкции фланцев из ПКМ: а - расслоение; б - смоляные карманы; в - поры; г - складки термография, ультразвуковой и рентгеновский контроль). Опыт изготовления образцов - сегментов конструкций с фланцевыми соединениями показывает, что возникновение таких дефектов, как правило, происходит равномерно вдоль зоны перегиба по всей ширине образца. Кроме того, эти дефекты видны на срезах препарированных конструкций после испытаний. Вопрос, насколько такой дефект будет протяженным по окружной координате при изготовлении конструкции, остается открытым. Можно предположить, что смоляные карманы и расслоения будут иметь ограниченную протяженность в окружном направлении, а складкообразование может быть распределено по всей окружности, например, в случае неправильной установки кольцевых упорных элементов формообразующей оснастки. Для оценки плотности материала в зоне дефекта было проведено исследование содержания связующего. Из дефектной конструкции (см. рис. 1) были взяты два типа проб материала: пробы первого типа использовались для оценки содержания связующего на масштабе 3-4 толщин слоев, пробы второго типа - для оценки содержания связующего по всей толщине конструкции. Пробы первого типа были получены в четырех точках (100 ± 0,2) мг в виде стружки после засверливания отверстий диметром 1,5 мм на глубину 30 мм: две пробы - непосредственно из зоны дефекта, и две пробы - на значительном удалении от этой зоны. Пробы второго типа были получены выпиливанием трех образцов по толщине конструкции размером приблизительно 25×10×10 мм, взятых из полок фланца и зоны перегиба с расслоением. Объемная доля связующего определялась с помощью термогравиметрического анализа на установке TGA/DSC 1/1100 LF (Mettler Toledo) в Институте технической химии УрО РАН. В результате анализа было установлено, что в окрестности дефекта массовое содержание связующего в слоях стеклопластика выше. Пробы первого типа в окрестности дефекта показали значения 37,1 и 41,0 %, проба второго типа - 35,4 %. В остальной части конструкции пробы первого типа показали значение 30,8 и 25,5 %, пробы второго типа - 30,0 и 26,0 % соответственно. Это отличие объясняется малым давлением формования в зоне перегиба слоев конструкции, что и приводит к образованию поры или смоляного кармана. В случае возникновения касательных или продольных напряжений в неотвержденном композиционном материале на этапе изготовления конструкции в этой зоне может произойти искривление слоев. Ниже изложены подходы к моделированию типовых дефектов типа искривления слоев, пор и смоляных карманов в композитных конструкциях, реализованные в работе. Решение задачи в двумерной постановке позволяет получить оценки НДС в кольцевом дефекте либо в средней зоне дефекта, достаточно протяженного по окружной координате, что представляет практический интерес. В дальнейших исследованиях может быть рассмотрено решение трехмерной задачи о фрагменте конструкции с локальным дефектом для уточнения полученного двумерного решения. В рассматриваемой конструкции выделим участок, включающий три ячейки звукопоглощающего контура (ЗПК) коробчатого типа и отгибной фланец (рис. 2, а). Построение геометрической модели и последующий анализ напряженно-деформированного состояния (НДС) такого образца проводились в пакете конечно-элементного анализа ANSYS в осесимметричной постановке (рис. 2, б). В состав детали входят слои стеклопластика с различной ориентацией основы тканого наполнителя - как вдоль образующей, так и поперек. Между основными слоями детали и ячейками ЗПК введен слой смолы. Крепежные отверстия, расположенные по окружности фланца, в расчетной схеме не рассматривались. Рис. 2. Геометрическая модель многослойного образца фланца: а - трехмерная; б - осесимметричная Складки (или волны в одном или нескольких слоях), как правило, образуются в процессе формования стыковочного фланца на цилиндрической части детали и в области перехода на отгибной (торцевой) участок [22]. Моделирование искривления слоев на базе исходной осесимметричной модели выполняется с помощью двухзвенной кривой средствами геометрического моделирования ANSYS, как, например, показано на рис. 3. При этом исходная геометрия и толщина конструкции не менялись, а уменьшалась толщина слоев в окрестности дефекта (см. рис. 2, 3) без изменения механических свойств материала. Таким образом, рассматривалась конструкция с меньшей суммарной толщиной слоев армирующего материала в окрестности дефекта. Целью исследования при этом являлась оценка влияния изгиба слоев в зоне дефекта на межслоевые напряжения, которые определяют прочность конструкции к растрескиванию и снижению жесткости фланца на изгиб. При этом оценка прочности конструкции при упрощенном описании дефекта с уменьшением толщины слоев в его окрестности получается с заведомо завышенным коэффициентом запаса. В настоящей работе дефект искривления слоев описывается двумя размерами, зависящими для удобства от одного параметра - толщины слоя h. Три слоя в данном примере на рис. 3, а (1-3) имеют постоянную толщину, остальные (4-7) поджимаются на половину толщины слоя в зоне максимального искривления (0,5h). Для правильного описания ортотропных свойств слоев в зоне искривления необходимо задать локальные системы координат (СК) для каждого из слоев. Таким образом, каждый из слоев 1а-7а и 1b-7b имеет свою локальную систему координат. Рис. 3. Искривление слоев: стыковочного фланца (а) и в зоне перехода цилиндрической части к фланцу (б) Более точно геометрически дефект в виде искривления слоев описывается с помощью многозвенной кривой. В этом случае увеличивается количество локальных СК, необходимых для правильного описания свойств материала. Таким же способом моделируется искривление слоев на внешней поверхности зоны перехода цилиндрической части детали к фланцу (рис. 3, б). В данном случае требуемое искривление достигается смещением узловых точек нескольких слоев по линии раздела слоев в зоне перегиба. Для каждого вновь образованного участка в области искривления также необходимо задать свою локальную СК. В обоих рассмотренных вариантах в зоне искривления слоев образуются пустоты - «поры». В реальных конструкциях они могут полностью или частично заполняться смолой и называются в этом случае «смоляными карманами». С точки зрения компьютерного моделирования области со смолой и поры одинаковы по геометрическим параметрам и отличаются лишь физико-механическим свойствами. На рис. 1 видно, что на цилиндрической части рассмотренных образцов присутствует продольный смоляной карман - область, заполненная смолой. Наличие такой прослойки может привести к утолщению цилиндрической части с некоторым искривлением внешних слоев. Рис. 4. Модели фланца с искривлением слоев и смоляными карманами в зонах перегибов (а) и в цилиндрической части фланца (б) На рис. 4 показаны модели, иллюстрирующие искривление слоев в наиболее характерных областях фланца с образованием пор или смоляных карманов. Ниже приведены результаты численного анализа НДС и оценки прочности фланца с использованием модели, показанной на рис. 4, а. При этом было рассмотрено два варианта: искривление слоев с образованием поры или смоляного кармана. 2. Математическая постановка задачи расчета НДС композитного фланца с дефектами Для расчета НДС и оценки запасов прочности рассматриваемой детали использовалась математическая постановка осесимметричной задачи теории упругости анизотропного неоднородного тела, система уравнений которой для подобных конструкций приведена авторами в работах [22, 25]. При задании граничных условий учитывалось, что в процессе эксплуатации на корпусные детали авиационного двигателя действует отрывающая осевая сила от реактивной струи сопла и изгибающий момент от присоединенной массы [22]. В предлагаемой методике в соответствии с общепринятым подходом для расчета фланцевых соединений принимается, что система внешних нагрузок сводится к суммарной растягивающей силе Р. Этой силе эквивалентна нагрузка р, равномерно распределенная по поверхности торцевых силовых слоев конструкции. Значение распределенной нагрузки p принималось равным 1 Н/м2. Так как в упругой постановке справедлива линейная зависимость значений напряжений в слоях конструкции от прикладываемых перемещений или усилий, то, анализируя результаты расчетов при единичной нагрузке, можно оценить запас прочности конструкции относительно любой величины эксплуатационной нагрузки. На противоположном конце фланца задаются площадки свободного опирания размером с диаметр шайбы болтового соединения. Схема граничных условий показана на рис. 5. Механические свойства стеклопластика, использованные в расчетах, были получены экспериментально в рамках предыдущих исследований авторов [26, 27] и представлены в табл. 1. Рис. 5. Граничные условия Таблица 1 Свойства стеклопластика, использованные при расчете Материал E11, ГПа E22, ГПа E33, ГПа G12, ГПа G13, ГПа G23, ГПа n12 n13 n23 Стеклопластик 24,6 18,6 6,0 4,0 3,0 3,0 0,15 0,42 0,18 Примечание: E11 - модуль упругости вдоль основы; E22 - модуль упругости вдоль утка; E33 - модуль поперек плоскости слоя; G12, G13, G23 - модули сдвига; n21, n13, n32 - коэффициенты Пуассона. Для областей, заполненных смолой (слой в районе ячеек ЗПК и дефекты в виде смоляных карманов), задавались свойства эпоксидной смолы с модулем упругости Е = 2,9 ГПа и коэффициентом Пуассона n = 0,356. В случае дефектов в виде пор вместо смолы задавался материал с существенно меньшим модулем упругости - 0,1 МПа. Такой подход принят в силу особенностей параметрического описания разработанных моделей, делающих нежелательной процедуру непосредственного удаления геометрических областей, имитирующих данные дефекты. Для оценки запаса прочности ортотропного материала каждого из слоев использовался критерий максимальных напряжений. Соответствующие пределы прочности приведены в табл. 2. Таблица 2 Пределы статической прочности материалов (МПа) Материал Стеклопластик 600 450 300 295 40 - - 45 45 Примечание: , , - пределы статической прочности материала на растяжение, сжатие и сдвиг соответственно в локальной системе координат слоя. При создании конечно-элементной модели использовался конечный элемент PLANE183 в четырехугольной модификации с включенной опцией осевой симметрии. Было проведено исследование сходимости результатов численного расчета при использовании трех типов конечно-элементных сеток с характерными средними размерами конечного элемента, равными h, 0,5h и 0,1h (h - толщина слоя), аналогично работе [25]. Оценка сходимости проводилась по эпюрам средних по элементу межслойных отрывающих напряжений, определяющих начало расслоения конструкции в наиболее опасных сечениях. Было показано, что переход от сетки с размерами конечного элемента, равными 0,5h, к сетке с размерами 0,1h приводит к отличию максимальных значений рассчитанных межслойных отрывающих напряжений не более чем на 1,5 %. Таким образом, для настоящих расчетов была выбрана сетка с характерным размером, равным половине толщины слоя (0,5h). Примерное количество конечных элементов в разработанных моделях составило порядка 70 тысяч. 3. Анализ НДС и оценка запаса статической прочности композитного фланца с дефектами В результате решения осесимметричной задачи теории упругости при силовом нагружении фланца получены поля напряжений и деформаций. На рис.6 приведена картина распределения осевых перемещений Uzz в глобальной системе координат конструкции. Характер деформирования фланца показывает, что заданные силовые условия нагружения создают изгибающий момент в верхних силовых слоях кожуха, что приводит к изгибу как цилиндрического, так и торцевого участков фланцевого соединения. Данный изгибающий момент и определяет наиболее опасные зоны в конструкции, которые показаны на рис. 6, б. В данной работе анализ НДС и оценку прочности многослойной композитной конструкции предлагается проводить по осредненным напряжениям, полученным в конечных элементах на сетке с характерным размером, равным половине толщины слоя, на основе анализа полей напряжений и графиков изменения межслойных напряжений (эпюр) вблизи границ между слоями. Запас прочности для межслойного характера разрушения слоистого пакета определялся с использованием пределов прочности при межслойном отрыве или сдвиге , . Запас прочности слоев на разрыв или сжатие вдоль волокон основы или утка определялся в результате сравнения напряжений с соответствующими предельными значениями для композита в направлении армирования , , , . Рис. 6. Поле осевых перемещений Uzz (мм) и наиболее нагруженные зоны композитного фланца Сопоставление результатов численного расчета НДС для моделей фланца с дефектами с соответствующими результатами для неповрежденной конструкции, полученными ранее в работе [25], позволяет оценить степень влияния различных типов дефектов на статическую прочность данного фланца. На рис. 7 показано изменение картины распределения полей осевых напряжений σzz (в глобальной СК) на примере зоны № 1 (см. рис. 6, б) при наличии дефектов. Видно, что с появлением поры или смоляного кармана происходит существенное перераспределение напряжений, при этом область локализации максимума в случае поры смещается к ее вершине. Рис. 7. Поля напряжений σzz (МПа) в глобальной СК в зоне №1 фланца: а - без дефекта, б - с дефектом в виде поры, в - с дефектом в виде смоляного кармана На рис. 8 показаны сечения, использованные для анализа межслойных напряжений в опасных зонах № 1-3 фланца (см. рис. 6, б). Заметим, что для первой зоны исследовались не только отмеченные сечения, но и более удаленные от поверхности. Однако в силу заметного затухания анализируемых напряжений в этих сечениях данные результаты в работе не приводятся. Для третьей зоны был проведен анализ межслойных напряжений для наружных наиболее нагруженных слоев, переходящих с цилиндрической части фланца на наружную оболочку кожуха и ориентированных под углом 55º к образующей оболочке. Напряжения в сечениях, направленных под углом к осям глобальной декартовой СК, необходимо рассматривать в локальных СК соответствующих слоев. Сечения № 18-20 (см. рис. 8) используются для оценки напряжений в слоях вдоль направления армирования (по основе или утку). Рис. 8. Сечения для построения эпюр межслойных напряжений: а - зона № 1, б - зона № 2, в - зона № 3 На рис. 9 приведены эпюры межслойных нормальных и касательных напряжений для выбранных сечений в зонах № 1 и № 3 для фланца с дефектом в виде поры. Аналогичные напряжения в области искривления слоев в зоне № 2 значительно ниже (в 10 и более раз), поэтому из рассмотрения их можно исключить. Анализ результатов расчета НДС показал следующее. Эпюры межслойных напряжений в различных слоях исследуемой области подобны - максимум имеет место, приблизительно, на расстоянии толщины слоя от линии перегиба слоев. Исключение составляют эпюры отрывающих и касательных напряжений между первым и вторым слоем (рис. 9, а, б, эпюра № 1 и рис. 9, в, г, эпюра № 9) в окрестности расположения поры. Для данных эпюр максимальные напряжения возникают в месте соединения слоев после «огибания» поры (см. рис. 8, а). Межслойные напряжения достигают наибольшего значения на торцевом участке фланца между вторым и третьим слоем, их значения составляют: для нормальных отрывающих - 3,21 МПа, для касательных - 1,35 МПа (см. рис. 9, а, б, эпюра № 2). Далее во внутренних слоях межслойные напряжения уменьшаются. Межслойные отрывающие напряжения в месте соединения первого и второго слоя после «огибания» поры достигают значения 2,56 МПа как на торцевой (рис. 9, а, эпюра № 1), так и на цилиндрической части фланца (рис. 9, в, эпюра № 9). Эти напряжения несколько ниже максимальных, использованных в работе для оценки нагрузки, определяющей начало расслоения конструкции. Сравнение с аналогичными напряжениями в исходной бездефектной конструкции показывает, что произошло изменение локализации области максимальных межслойных напряжений в торцевой части фланца и их незначительное увеличение: отрывных напряжений - на 6 % (см. рис. 9, а), касательных напряжений - на 2 % (см. рис. 9, б). Для цилиндрического участка фланца изменение напряженного состояния оказалось существеннее - повышение межслойных нормальных напряжений до 1,5 раз (см. рис. 9, в), касательных напряжений - до 1,3 раз (см. рис. 9, г). При таком перераспределении напряжений произошла разгрузка участка перехода с цилиндрической части фланца на наружную оболочку ЗПК (зона № 3), характеризуемого эпюрами № 21-22 (см. рис. 9, д-е). Снижение межслойных нормальных напряжений составило порядка 7 %, касательных напряжений - 10 %. Рис. 9. Эпюры межслойных нормальных (а, в, д) и касательных (б, г, е) напряжений в локальной СК слоя для фланца с дефектом в виде поры. Цифра у кривых соответствует номеру сечения на рис. 8. Пунктиром обозначена эпюра напряжений в сечении соответствующего номера для бездефектного фланца [25]. По оси абсцисс откладывается расстояние вдоль слоя от линии перегиба слоев (см. рис. 8) Таким образом, наименьший запас прочности для рассматриваемой конструкции с дефектом в виде поры наблюдается в области искривления слоев в вершине поры и определяется величиной межслойных нормальных растягивающих напряжений, возникающих между вторым и третьим слоем торцевого участка фланца (эпюра № 2). Принимая в качестве предельного значения отрывающих межслойных напряжений = 40 МПа, получим соответствующее значение предельной распределенной нагрузки p, вызывающей расслоение на данном участке, равной 12,5 МПа. Запас прочности конструкции по межслойным сдвиговым напряжениям более чем в два раза выше запаса прочности по нормальным межслойным напряжениям. Как было показано ранее [25], в исходном образце фланца предполагаемое разрушение по механизму расслоения ожидается в зоне № 3 при нагрузке порядка 12 МПа. Таким образом, рассматриваемый дефект в виде поры приводит к перераспределению напряжений в конструкции со сменой области прогнозируемого начала разрушения, при незначительном (4 %) повышении общей несущей способности. Следует отметить, что анализируемые поля напряжения являются осредненными по половинной толщине слоя. Дополнительные расчеты исследуемого фланца на сетках со сгущением в окрестности соединения слоев после «огибания» поры показали наличие высокой концентрации напряжений. Таким образом, полученные осредненные напряжения требуют последующего уточнения для оценки возможности появления сингулярности в данной зоне. Кроме того, представляет интерес последующая постановка и решение задачи о межслоевой трещине в окрестности поры в зоне концентрации напряжений. Это обстоятельство может привести к снижению расчетного значения нагрузки, определяющей начало расслоения исследуемой конструкции с дефектом. С другой стороны, следует учитывать, что частичное расслоение приводит, как правило, к некоторому снижению жесткости конструкции, которая может оставаться работоспособной и обладать еще значительным запасом прочность по напряжениям в плоскости армирующих слоев. Для фланца с дефектом в виде смоляного кармана отмечается равномерное падение межслойных напряжений во всех рассмотренных зонах по сравнению с бездефектной конструкцией [25]. Причем для зоны № 1 снижение отрывных напряжений составило 40-50 %, касательных напряжений - до 70 %. В зоне № 3 уменьшение рассматриваемых напряжений составило порядка 4-6 %. Эта зона и является наиболее опасным участком конструкции, а ее предельное состояние будет достигаться при величине распределенной нагрузки p, равной 12,5 МПа, как и в случае с дефектом в виде поры. Можно отметить, что в рассматриваемом случае в окрестности соединения слоев после «огибания» смоляного кармана концентрация напряжений отсутствует. Это объясняется тем, что модули Юнга смолы и стеклопластика в направлении, перпендикулярном слою E33, отличаются только в два раза. На рис. 10 представлены эпюры осевых напряжений, действующих в поперечных сечениях фланца № 18-20 в плоскости армирующих слоев. Данные напряжения отражают тот факт, что цилиндрический участок фланца подвержен осесимметричному изгибу вследствие изгибающего момента. Сравнение результатов, полученных для двух рассмотренных расчетных случаев, с исходными данными для фланца без дефектов (рис. 10, а) показывает, что наличие пор в зонах № 1 и № 2 приводит к небольшому, порядка 8 %, снижению максимальных растягивающих напряжений в поверхностных слоях (рис. 10, б). Интересно отметить, что в случае образования смоляных карманов в зоне № 3 вблизи перехода с цилиндрической части фланца на наружную оболочку ЗПК происходит резкое (на 50 %) увеличение растягивающих напряжений в наружных слоях - эпюра № 20 Рис. 10. Эпюры продольных напряжений в глобальной СК в слоях фланца: а - без дефектов; б - с дефектом в виде поры; в - с дефектом в виде смоляного кармана. Цифра у кривых соответствует номеру сечения на рис. 9 (рис. 10, в). Однако, как и для исходной конструкции фланца без дефектов, запас прочности по напряжению в плоскости слоев на данном участке минимум в 5 раз выше аналогичных оценок для межслойного характера разрушения. Заключение Проведено исследование влияния технологических дефектов: искривления слоев, пор и смоляных карманов на НДС и статическую прочность типового композитного фланца, являющегося составной частью звукопоглощающего кожуха авиадвигателя. Рассмотренные дефекты могут возникать в подобных конструкциях вследствие отклонения параметров технологического процесса, при этом в настоящее время не существует надежной методики их выявления средствами неразрушающего контроля. Целью работы являлась оценка влияния дефекта подобного типа на НДС и прочность конструкции, определяемую величиной поперечных и сдвиговых межслойных напряжений. Из предыдущих исследований [22-25] известно, что именно эти напряжения обусловливают начало растрескивания и расслоения в зоне перегиба слоев фланца, где и возникают технологические дефекты. При моделировании дефектов учитывалось изменение геометрии слоя в окрестности дефекта с сохранением его эффективных механических свойств. Для сохранения исходного контура конструкции общая толщина армирующего пакета слоев в локальной зоне дефекта уменьшалась на величину дефекта. Таким образом, оценка прочности конструкции с дефектом осуществлялась с запасом. Задача решалась на двумерной осесимметричной модели, что соответствует равномерному распределению дефекта в конструкции по окружной координате. Решение задачи осуществлялось численно в среде ANSYS с учетом анизотропии механических свойств материала в упругой постановке, оценка статической прочности проводилась по критерию максимальных напряжений. Результаты расчета сравнивались с аналогичными данными для фланца без дефектов, полученными в работе [27]. Проведенные расчеты показали, что даже при задании уменьшенной толщины слоя изменение геометрии в зоне образования поры или смоляного кармана не приводит к повышению межслойных напряжений по сравнению с бездефектной конструкцией. Таким образом, рассмотренные дефекты незначительно снижают статическую прочность подобных конструкций. В последующих исследованиях представляет интерес оценка влияния разброса механических свойств материала слоев в окрестности дефекта на уровень напряжений в слоях, оценка сингулярности напряжений в окрестности смыкания слоев вокруг поры, а также оценка стойкости материала к распространению трещин. Это обстоятельство может привести к снижению расчетного значения нагрузки, определяющей начало расслоения рассматриваемой конструкции с дефектом в виде поры. При этом отметим, что частичное расслоение конструкции, как правило, приводит только к снижению жесткости конструкции, которая может оставаться работоспособной и обладать еще значительным запасом прочности по напряжениям в плоскости армирующих слоев.About the authors
A N Anoshkin
Perm National Research Polytechnic University
V Yu Zuiko
Perm National Research Polytechnic University
V M Osokin
Perm National Research Polytechnic University
A A Tretyakov
Perm National Research Polytechnic University
P V Pisarev
Perm National Research Polytechnic University
References
- Сапожников C.Б. Дефекты и прочность армированных пластиков / Челяб. гос. техн. ун-т. - Челябинск, 1994. - C. 162.
- Defects in composite structures: Its effects and prediction methods - A comprehensive review / K. Senthil, A. Arockiarajan, R. Palaninathan, B. Santhosh, K.M. Usha // Composite Structures. - 2013. - Vol. 106. - P. 139-149. doi: 10.1016/j.compstruct.2013.06.008
- Бохоева Л.А. Особенности расчёта на прочность элементов конструкций из изотропных и композиционных материалов с допустимыми дефектами. - Улан-Удэ: Изд-во Вост.-Сиб. технол. ун-та, 2007. - С. 192.
- Болотин В.В. Дефекты типа расслоений в конструкциях из композитных материалов // Механика композитных материалов. - 1984. - № 2. - С. 239-255.
- Болотин В.В. Межслойное разрушение композитов при комбинированном нагружении // Механика композитных материалов. - 1988. - № 3. - С. 410-418.
- Болотин В.В. О динамическом распространении трещин // Прикладная математика и механика. - 1992. - Т. 56, № 1- C. 150-162.
- Трошин В.П. Влияние продольного расслоения в слоистой цилиндрической оболочке на величину критического внешнего давления // Механика композитных материалов. - 1982. - № 5. - C. 838-842.
- Трошин В.П. К устойчивости цилиндрических оболочек с расслоениями // Механика композитных материалов. -1981. - № 4. - C. 729-731.
- Тарнопольский Ю.М. Расслоение сжимаемых стержней из композитов // Разрушение композитных материалов.- 1979. - № 2. - С. 160-166.
- Доан Чак Луат, Лурье С.А., Дудченко А.А. Моделирование деградации свойств композита при растрескивании и расслоении при статическом и циклическом нагружении // Механика композиционных материалов и конструкций. - 2008. - T 14, № 4. - С. 623-637.
- Bottega W.J., Maewal A. Delamination buckling and growth in laminates // Journal Applied Mechanics. - 1983. - Vol. 50. - No 1. - P. 184-189. doi: 10.1115/1.3166988
- Chai H., Babcock C.D., Knauss W.G. One dimensional modeling of failure in laminated plates by delamination buckling // International Journal of Solids and Structures. - 1981. - Vol. 27. - No. 11. - P. 1069-1083. doi: 10.1016/0020-7683(81)90014-7
- Chai H., Babcock C.D. Two-dimensional modeling of compressive failure in delaminated laminates // Journal of Composite materials. - 1985. - Vol. 19. - No. 1. - P. 67-91. DOI: 10.1177/ 002199838501900105
- Bolotin V.V. Delaminations in composite structures: its origin, buckling, growth and stability // Composites Part B: Engineering. - 1996. - Vol. 27. - No. 2. - P. 129-145. doi: 10.1016/1359-8368(95)00035-6
- Simitses G.J., Sallam S., Yin W.H. Effect of delamination of axially loaded homogeneous laminated plates // AIAA Journal. - 1985. - Vol. 23. - No. 9. - P. 1437-1444. doi: 10.2514/3.9104
- Muc A, Stawiarski A. Identification of damages in composite multilayered cylindrical panels with delaminations // Composite Structures. - 2012. - Vol. 94. - No. 5. - P. 1871-1879. DOI:10.1016/ j.compstruct.2011.11.026
- Prediction of flange debonding in composite stiffened panels using an analytical crack tip element-based methodology / Z. Mikulik, D.W. Kelly, B.G. Prusty, R.S. Thomson // Composite Structures. - 2008. - Vol. 85. - No. 3. - P. 233-244. doi: 10.1016/j.compstruct.2007.10.027
- Jian Li. Flange delamination prediction in composite structures with ply waviness // AIAA Journal. - 2000. - Vol. 38. - No. 5. - P. 893-897. doi: 10.2514/2.1044
- Camanho P.P., Davila C.G., Pinho S.T. Fracture analysis of composite co-cured structural joints using decohesion elements // Fatigue & Fracture of Engineering Materials & Structures. - 2004. - Vol. 27. - No. 9. - P. 745-757. doi: 10.1111/j.1460-2695.2004.00695.x
- Steeves C.A., Fleck N. A. Compressive strength of composite laminates with terminated internal plies // Composites Part A: Applied Science and Manufacturing. - 2005. - Vol. 36. - No. 6. - P. 798-805. doi: 10.1016/j.compositesa.2004.10.024
- McElroy M., Leone F., Ratcliffe J. Simulation of delamination-migration and core crushing in a CFRP sandwich structure // Composites Part A: Applied Science and Manufacturing. - 2015. - Vol. 79. - P. 192-202. doi: 10.1016/j.compositesa.2015.08.026
- Аношкин А.Н., Ташкинов А.А. Прогнозирование несущей способности композитных фланцев корпусных деталей авиадвигателей / Перм. гос. техн. ун-т. - Пермь, 1998. - С. 101.
- Аношкин А.Н., Ташкинов А.А. Нестационарные процессы накопления повреждений композитных фланцев при циклических нагрузках // Механика композитных материалов. - 1997. - Т. 33, № 5. - С. 449-454.
- Аношкин А.Н., Ташкинов А.А., Грицевич А.М. Прогнозирование несущей способности композитных фланцев корпусных деталей авиадвигателей // Механика композитных материалов. - 1997. - Т. 33, № 3. - С. 255-262. doi: 10.1007/s00158-010-0617-4
- Расчет НДС и оценка прочности композитного фланца стеклопластикового кожуха авиационного газотурбинного двигателя / А.Н. Аношкин, М.В. Рудаков, И.С. Страумит, Е.Н. Шустова // Вестн. Уфим. гос. авиац. техн. ун-та. - 2011. - Т. 15, № 1 (41). - С. 67-75.
- Repair of damage in aircraft composite sound-absorbing panels / A.N. Anoshkin, V.Y. Zuiko, M.A. Tashkinov, V.V. Silberschmidt // Composite Structures. - 2015. - Vol. 120. - P. 153-166. doi: 10.1016/j.compstruct.2014.10.001
- Experimental-theoretical research of mechanical properties of perforated composite sandwich panels / A.N. Anoshkin, V. Yu. Zuiko, A.V. Tchugaynova, E.N. Shustova // Solid State Phenomena. - 2016. - Vol. 243. - P. 1-10. doi: 10.4028/www.scientific.net/SSP.243.1