Numerical simulation for development of methodology of stress-strain state control of composite bulkhead for aviation application with the usage of FBG sensors
- Authors: Anoshkin AN1, Zuiko VY.1, Pelenev KA1, Pisarev PV1, Shipunov GS1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 4 (2018)
- Pages: 47-57
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/479
- DOI: https://doi.org/10.15593/perm.mech/2018.4.04
- Cite item
Abstract
Fiber-optic sensors are used to monitor the state of structures in various industries. The best solution for monitoring deformations in the products made of polymer composite materials are fiber Bragg gratings (FBGs). In this paper, we consider the possibility of studying the stress-strain state (SSS), tracking structural changes in the construction with the help of FBGs, using the example of a "U" -shaped composite frame of an aircraft engine. Based on the revealed regularities and peculiarities of the FBGs operation, the main control points have been determined. In the course of the research, we created a three-dimensional computer model of an aircraft propulsion system composite frame which was used for the analysis of the stress-strain state under various loading conditions. The developed model makes it possible to perform a layer-by-layer analysis of the composite frame structure to estimate the normal and shear interlayer stresses, which determine the destruction of the structure. Numerical modeling of this problem was carried out by the finite element method (FEM) in the software package ANSYS Workbench. Numerical results were compared with the experimental data obtained in laboratory tests of a composite frame segment equipped with fiber-optic sensors. The tests were carried out on the equipment of the unique scientific installation "Unique scientific and technological complex of automated computation" using the Zwick servo-hydraulic machine and the Astro X327 interrogator. The results obtained will be used to determine the size of the effective "sensitive zone" for embedding fiber-optic sensors (FOS) based on Bragg gratings into the construction of a composite frame. The main goal of this work is the development of a methodology for detecting and identifying various types of damage and the maximum load during the working lifespan of the structure.
Full Text
Введение В последнее время полимерные композиционные материалы (ПКМ) начинают широко применяться в авиационных двигательных установках. Опыт использования таких материалов в подобных конструкциях, подверженных воздействию сложных циклических нагрузок в течение длительного периода времени (более 15 лет), невелик. Поэтому большой интерес вызывает разработка систем мониторинга состояния таких деталей для использования как во время стендовых и летных испытаний, так и в процессе эксплуатации. Целью настоящей работы является разработка концептуальной схемы мониторинга и контроля НДС композитного шпангоута авиационного двигателя с использованием оптоволоконных датчиков. Композитный силовой шпангоут представляет собой кольцо П-образного профиля (рис. 1) и предназначен для крепежа решеток реверсивного устройства авиационного двигателя. Шпангоут изготавливается из текстильного углепластика на эпоксидном связующем по технологии пропитки под давлением (RTM). Конструктивно в шпангоуте можно выделить сплошную кольцевую зону П-образной полки и два фланца, которые на ряде участков имеют фигурные вырезы и крепежные отверстия. a b Рис. 1. Общий вид модели (a) и сегмент натурного шпангоута (б) реверсивного устройства Fig. 1. General view of the model (a) and the segment of the bulkhead (b) of the reversing device К одному из фланцев шпангоута крепятся решетки реверсивного устройства, передающие нагрузки в осевом и радиальном направлении от газодинамических сил. Нагрузки распределены неравномерно как по радиальной, так и окружной координате, что приводит к изгибу шпангоута и возникновению не только высоких напряжений в плоскости армирования слоев, но и высоких межслоевых напряжений (рис. 2). Подробно задача расчета НДС и оценки прочности шпангоута рассмотрена ранее в работах [1-2]. Эти напряжения могут вызвать появление межслоевых трещин, их развитие может привести к снижению жесткости шпангоута, что является опасным для нормального функционирования реверсивного устройства двигателя. Поэтому мониторинг и контроль напряженно-деформированного состояния шпангоута является актуальной задачей. Рис. 2. Результаты численного моделирования шпангоута реверсивного устройства при воздействии эксплуатационных нагрузок Fig. 2. Simulation results of the bulkhead of the reversing device under operational loads Для решения задач мониторинга состояния авиационных конструкций в последнее время активно используются оптоволоконные датчики различного типа [3-14]. Наиболее распространенными из них, обеспечивающими хорошую точность измерений с приемлемой стоимостью как самих датчиков, так и регистрирующей аппаратуры, являются датчики на основе решеток Брэгга (ВБР) (рис. 3) [15-19]. Для ВБР-датчиков разработаны компактные модели интеррогаторов с малым весом и высокой надежностью, что позволяет размещать их на борту воздушного судна [20-24] и создавать, таким образом, систему мониторинга состояния конструкции в процессе эксплуатации. Рис. 3. Схема работы волоконной Брэгговской решетки Fig. 3. Operation scheme of a fiber Bragg grating При создании такой системы одной из основных задач является определение схемы расположения датчиков в конструкции и разработка методики оценки напряженно-деформированного состояния по показаниям датчиков. Датчики должны позволять оценивать напряженно-деформированное состояние, в том числе в наиболее нагруженных зонах, и обеспечивать оценку запаса прочности по регистрируемым данным. Для большинства оболочечных изделий, изготовленных из металлических сплавов, линейные и угловые деформации, регистрируемые на поверхности в различных точках, позволяют оценить максимальные напряжения и запас прочности по интегральным критериям. Для мониторинга состояния неоднородных слоистых композитных конструкций существует необходимость оценки линейных и сдвиговых межслойных деформаций и напряжений, перпендикулярных плоскости слоя. Несмотря на то, что значения указанных деформаций и напряжений невелики, предел прочности на межслойный отрыв и сдвиг слоистых композиционных материалов также в десятки раз меньше пределов прочности в направлении армирования. Особенно опасными данные напряжения становятся в местах перегибов слоев, складок, а также при изгибе. По типу расслоения от действия межслойных напряжений часто разрушаются фланцевые узлы композитных оболочек, многослойные сотовые конструкции [25-30]. Оценка межслойных напряжений является актуальной задачей для системы мониторинга рассматриваемого композитного силового шпангоута авиационного двигателя. Следует отметить, что стандартные ВБР-датчики позволяют регистрировать деформации только вдоль оптического волокна. Поэтому при обычном размещении на поверхности или внутри многослойной композитной конструкции с помощью датчиков можно получить только деформацию в плоскости армирующих слоев. Межслоевые деформации и напряжения можно при этом оценить косвенно на основе расчета напряженно-деформированного состояния конструкции. Считая деформирование шпангоута упругим, а также учитывая единственность упругого решения задачи, можно найти зависимость значения межслойных напряжений в наиболее нагруженных зонах шпангоута от других компонентов напряжений или деформаций, регистрируемых оптоволоконными датчиками по результатам анализа полей напряжений и деформаций при тестовой нагрузке. Эта зависимость будет справедливой при любом пропорциональном изменении нагружающих усилий относительно тестовых значений для шпангоута с заданной схемой армирования до появления расслоений. Таким образом, можно приблизительно определять значения межслойных напряжений при эксплуатационной нагрузке, если тестовая нагрузка ей пропорциональна, и оценивать запас прочности конструкции до расслоения. При возникновении и развитии межслойных трещин поля напряжений и деформаций в слоях композитной конструкции будут отличаться от тех, которые были получены в результате решения упругой задачи, и установленная зависимость будет нарушена. Жесткость конструкции на изгиб при этом будет изменяться, что приведет к изменению деформаций в плоскости слоев и показаний ВБР-датчиков в некоторой окрестности зоны расслоения. Таким образом, по показаниям ВБР-датчиков, расположенных в плоскости слоя, можно фиксировать возникновение и развитие расслоений в конструкции и оценивать изменение её интегральной жесткости. При этом датчики должны быть расположены в зоне, имеющей наибольшие межслойные напряжения, и ориентированы преимущественно в направлении наибольших линейных деформаций. В настоящей работе представлены результаты расчетных и экспериментальных исследований для оценки возможности мониторинга НДС сегмента композитного шпангоута с помощью ВБР-датчиков при лабораторных испытаниях. Проведена расчетно-экспериментальная оценка возможности регистрации возникновения межслоевых трещин при лабораторных механических испытаниях сегмента шпангоута с помощью ВБР-датчиков, расположенных в плоскости слоев. 1. Механические испытания сегмента композитного шпангоута с ВБР-датчиками Сегмент шпангоута - конструктивно подобный элемент натурной детали (рис. 4) был изготовлен из текстильного углепластика со схемой армирования [0°, 90°] по технологии RTM, аналогично используемой при изготовлении натурной детали. Использовалась равнопрочная ткань, основа ткани ориентировалась в осевом либо кольцевом направлении сегмента. Механическое нагружение сегмента проводилось на испытательной машине Zwick Z100 Pro Line, входящей в состав «Уникального научно-технологического комплекса автоматизированной выкладки» ПНИПУ, с помощью специально изготовленной оснастки, обеспечивающей жесткое закрепление нижней полки сегмента и приложение вертикальной растягивающей нагрузки к верхней полке. Для контроля деформированного состояния сегмента на наружную поверхность П-образной полки шпангоута были приклеены две оптоволоконные линии. Первая линия, ориентированная в кольцевом направлении, содержала четыре датчика, позволяющих контролировать деформацию сегмента шпангоута в кольцевом направлении. Вторая линия включала участок, ориентированный в осевом направлении, возвратную петлю и участок, ориентированный под 45° к осевому направлению. На каждом участке второй линии содержалось по четыре датчика, что позволяло контролировать осевую деформацию и деформацию под 45° соответственно. В лабораторных испытаниях желательно обеспечивать фиксацию оптоволоконных линий на поверхности испытываемой детали с возможностью последующего снятия. Поэтому на данном этапе исследования для этого использовалась двухсторонняя липкая лента авиационного назначения SK2AT200-1. Рис. 4. Схема испытаний и расположения датчиков ВБР на наружной поверхности П-образной полки сегмента шпангоута Fig. 4. Tests (a) and location of FBG sensors on the outer surface of the U-shaped flange of the bulkhead segment Волоконно-оптические датчики на основе Брэгговских решеток, предоставленные для проведения исследований индустриальным партнёром ООО «Инверсия-Сенсор» и изготовленные методом фазовой маски на оптическом волокне марки SMF-28 эксимерным лазером, имели следующие характеристики: коэффициент отражения (70±3)%; ширина пика на полувысоте (0,2±0,01) нм; уровень боковых пиков (9,5±0,3) дБ.; длина решетки 6,5 мм; резонансные длины волн опрашиваемых датчиков представлены в табл. 1. Для регистрации деформаций с помощью датчиков использовался интеррогатор ASTRO327 с точностью определения деформаций (±4·10-5). Испытания сегмента проводились при комнатной температуре +22°C с постоянной скоростью нагружения 5 мм/мин. В процессе испытаний регистрировалось усилие нагружения, перемещение захватов испытательной машины, и снимались показания с оптоволоконных датчиков. На рис. 5 приведены диаграммы нагружения сегмента и синхронизированные с ней осевые деформации на поверхности П-образной полки, регистрируемые датчиками в осевом, кольцевом и направлении под углом 45° к осевому направлению. Диаграмма нагружения имеет начальный участок с изменяющимся наклоном, который обусловлен смещением верхнего фланца сегмента шпангоута в испытательной оснастке при нагружении. Оценка этого смещения приведена в работе [2]. Кроме того, на диаграмме имеются скачки, свидетельствующие о появлении расслоений в композиционном материале и об изменении жесткости конструкции. Следует отметить, что в ряде сигналов с ВБР, регистрируемых в процессе испытаний, присутствовали значительные шумы, что, вероятно, вызвано нарушением контакта между волокном и поверхностью образца. В последующих испытаниях композитных образцов другого типа [31] данная проблема была решена за счёт замены клеевой ленты на цианоакрилатный клей с каучуковым наполнением, что позволяло получать сигналы с ВБР-датчиков, наклеенных на поверхность композитных образцов без шумов, и легко снимать оптоволоконные линии после проведения испытаний, используя растворитель. Таблица 1 Резонансные длины волн решеток Брэгга ВБР-датчиков Table 1 Resonant wavelengths of Bragg gratings FBG sensors Обозначение датчика B009 B010 B011 B012 B103 Резонансная длина волны Брэгга, нм 1530,842 1524,023 1517,008 1510,026 1571,88 Обозначение датчика B104 B105 B111 B112 - Резонансная длина волны Брэгга, нм 1565,514 1558,667 1516,784 1509,875 - a б в г Рис. 5. Диаграммы нагружения сегмента шпангоута и деформации от датчиков ВБР в осевом направлении (а); в кольцевом направлении (б); под углом 45° к осевому направлению (в); (г) характер разрушения шпангоута при механических испытаниях при нагрузке 1,7 кН [2]. Характерные значения нагружающего усилия: 1,27 кН (А); 1,39 кН (В); 1,7 кН (С) Fig. 5. Loading diagrams of the bulkhead segment and deformations from the FBG sensors in the axial direction (a); in the circumferential direction (b); the direction at an angle of 45° to the axial direction (c) and the type of destruction of the bulkhead during mechanical tests under 1.7 kN (d) [2]; representative values of the loading force: 1.27 kN (A); 1.39 kN (B); 1.7 kN (C) Появление первого макроскопического расслоения в конструкции с выходом межслоевой трещины на боковую часть образца фиксируется при нагрузке 1,27 кН. При этом на диаграмме нагружения (см. рис. 5) регистрируется первое скачкообразное падение усилия, отражающее мгновенное перераспределение поля деформаций. При дальнейшем увеличении нагрузки до 1,39 кН на диаграмме нагружения наблюдается второе скачкообразное падение, при этом на боковой поверхности конструкции отмечается разрастание образовавшегося ранее расслоения и образование новых менее явных межслоевых трещин. При достижении нагрузки в 1,70 кН (см. рис. 5) регистрируется последнее скачкообразное падение усилия, при этом размер расслоения достигает критического значения, при котором происходит раскрытие шпангоута и потеря его несущей способности. Далее сегмент шпангоута начинает разгибаться при усилии приблизительно 1,2 кН. Следует отметить, что после каждого скачкообразного падения усилия меняется наклон следующего участка диаграммы нагружения, что свидетельствует об изменении жесткости конструкции вследствие расслоений. Указанные расслоения регистрировались ранее при проведении рентгеновского контроля in-situ в процессе проведения аналогичных механических испытаний сегмента шпангоута [2]. На графиках деформаций, регистрируемых датчиками ВБР, в соответствующих точках также видны скачки (см. рис. 5), выделяемые из общего шума. После этого наклон графиков изменялся, что также отражает изменение жесткости конструкции. Особенно отчетливо указанные скачки видны на графиках деформаций в осевом направлении (вдоль нагружающего усилия) и в направлении под 45° к оси. Результаты эксперимента сравнивались с результатами математического моделирования лабораторных испытаний сегмента шпангоута, представленного в следующем разделе. 2. Математическое моделирование сегмента композитного шпангоута, сравнение с результатами механических испытаний Моделирование лабораторных испытаний сегмента шпангоута проводилось в результате решения упругой краевой задачи, расчетная схема которой представлена на рис. 6. Рассматривался сегмент композитного шпангоута с явным описанием армирующих слоев. Это позволяет проводить оценку напряжений как в плоскости слоя, так и межслоевых отрывных и сдвиговых напряжений, которые определяют начало разрушения слоистой конструкции и потерю её жесткости. В модели использовалось две системы координат: глобальная система координат конструкции OXYZ (см. рис. 6) и локальная система координат слоя, где ось OX1 совпадает с направлением нитей основы, ось OX2 с направлением нитей утка, ось OX3 перпендикулярна плоскости слоя. Кроме композитного сегмента шпангоута в расчетную схему был включен участок стальной оснастки, контактирующий с верхним фланцем шпангоута. В качестве граничных условий краевой задачи принималась равномерно распределенная нагрузка вдоль стальной оснастки и жесткое закрепление нижнего фланца шпангоута по площадкам контакта в отверстиях болтового соединения по границе 1 (рис. 6, б). Эквивалентное значение суммарного нагружающего усилия принималось равным 1 кН для удобства последующего анализа результатов решения упругой задачи. Для моделирования смещения верхнего фланца шпангоута в испытательной оснастке в процессе механических испытаний на площадке контакта верхнего фланца шпангоута (по границе 2) и испытательной оснастки задавались контактные граничные условия с трением, коэффициент трения принимался равным 0,3. а б Рис. 6. Схема испытаний (а) и расчетная схема (б) для сегмента шпангоута: 1 - жесткое закрепление в области болтовых соединений; 2 - нагрузка 1 кН Fig. 6. Testing scheme (a) and calculation scheme (b) for the bulkhead segment. Scheme of boundary conditions: 1 - bolting; 2 - 1 kN load Решение упругой задачи осуществлялось методом конечных элементов с использованием программного комплекса ANSYS Workbench. Упругие и прочностные свойства анизотропных слоев углепластика, использованные в расчетах, приведены табл. 2, 3. Свойства слоя текстильного углепластика в плоскости армирования 0X1X2 принимались по результатам испытаний стандартных образцов, свойства перпендикулярно плоскости армирования принимались по оценочным данным из работ [33]. Таблица 2 Упругие свойства композиционного материала Table 2 Elastic properties of the composite material E11, ГПа E22, ГПа E33, ГПа G12, ГПа G13, ГПа G23, ГПа n21 n13 n32 Равнопрочный углепластик 63,9 63,9 20,0 19,5 2,7 2,7 0,04 0,3 0,3 Таблица 3 Прочностные свойства композиционного материала Table 3 Strenght properties of composite material , МПа , МПа , МПа , МПа , МПа , МПа , МПа , МПа , МПа Равнопрочный углепластик 809 804 809 804 44 128 150 77 77 Примечание: E11 - модуль упругости вдоль основы; E22 - модуль упругости поперек основы (вдоль утка); E33 - модуль поперек плоскости слоя; G12, G13, G23 - модули сдвига;n21, n13, n32 - коэффициенты Пуассона. Sii+, Sii-, Sij - пределы статической прочности материала на растяжение, сжатие и сдвиг в локальной системе координат слоя. При создании конечно-элементной модели шпангоута использовались твердотельные трехмерные восьмиузловые элементы Solid185. Исследование сходимости численного решения показало, что максимальные напряжения в слоях изменяются не более чем на 5% при использовании конечных элементов размером не менее половины толщины слоя h/2 (0,125 мм). Общее количество конечных элементов, используемых для дискретизации объекта, составило порядка 2 миллионов. На рис. 7 представлен общий вид конечно-элементной сетки, используемой в расчетах. а б Рис. 7. Общий вид (а) и фрагмент (б) конечно-элементной сетки сегмента шпангоута Fig. 7. General view (а) and a fragment (b) of the finite element mesh of the bulkhead segment Математическая постановка задачи расчета НДС рассматриваемого образца, включающая соотношения теории упругости неоднородного анизотропного тела, приводилась авторами при расчетах подобных конструкций в работах [32, 33]. В вариационной формулировке эта постановка для области V c границей S заключается в поиске минимума функционала Лагранжа, вариация которого имеет вид , (1) где εij и δεij - тензор и вариация тензора деформаций; Cijkl - тензор упругих модулей; δui - вариация вектора перемещений; Fi - вектор внешних сил, действующих на участке границы St. На рис. 8 показано поле осевых деформаций εхх на поверхности сегмента шпангоута при нагрузке 1 кН. Сравнение полей деформаций, полученных в результате численного расчета, с данными, которые зарегистрированы датчиками при механических испытаниях, приведены в табл. 4. Рис. 8. Поле осевых деформаций ε11 на поверхности сегмента шпангоута при нагрузке 1 кН Fig. 8. The axial strain field ε11 on the surface of the bulkhead segment at a load of 1 kN Анализируя данные табл. 4, можно отметить следующее: наибольшее расхождение (свыше 90%) с результатами расчетов показывают датчики В009, В010, В011, В012, регистрирующие кольцевые деформации, значения которых при выбранной схеме испытаний меньше погрешности используемого интеррогатора (±4·10-5). Таким образом, можно сказать, что показания датчиков В009, В010, В011, В012 при нагрузке 1 кН недостоверны. Показания датчиков В103, В104 и В112, регистрирующих деформацию под углом 45° и осевую деформацию, удовлетворительно соответствуют расчетным значениям (отклонение от 1,5 до 15%). Деформации, регистрируемые датчиком В105, находятся на границе погрешности интеррогатора, что может объяснять их отличие от расчетных значений на 47%. Датчик В111 характеризуется высокой зашумленностью сигнала, что объясняется плохой приклейкой оптического волокна на поверхность шпангоута. Поэтому показания датчика В111 расходятся с расчётными значениями на 23%. Следует отметить, что диаграммы деформации, регистрируемые всеми указанными датчиками независимо от отклонений, расчетных значений, фиксируют изменение жесткости конструкции при возникновении расслоений. 3. Оценка прочности композитного шпангоута Появление первой макроскопической трещины в сегменте шпангоута при испытаниях зафиксировано при нагрузке 1,27 кН, при этом на диаграмме нагружения и диаграммах, регистрируемых ВБР, наблюдаются скачки (см. рис. 5). Значения деформаций при этой нагрузке, как зарегистрированных ВБР, так и рассчитанных в результате математического моделирования, приведены в табл. 5. Таблица 4 Деформации по показаниям ВБР при нагрузке 1 кН и расчетные значения деформаций в окрестностях датчиков, направленные по оси чувствительности ВБР Table 4 Deformations according to the readings of the FBG sensors under 1 kN and the calculated values of deformations in the vicinity of the sensors directed along the sensitivity axis of the FBG sensors Параметр Кольцевое направление Под углом 45° к осевому направлению Осевое направление B009 B010 B011 B012 B103 B104 B105 B111 B112 Деформации по показаниям ВБР·10-3 -0,15 0,055 0,055 0,13 -0,73 -0,81 -0,13 -0,90 -0,69 Расчетные значения деформаций·10-3 0,001 0,0005 0,0028 0,011 -0,70 -0,69 -0,07 -0,69 -0,70 Отклонение, % 99,33 99,05 94,91 91,54 4,11 14,81 46,92 23,33 1,45 Таблица 5 Деформации, регистрируемые ВБР в момент разрушения при нагрузке 1,27кН Table 5 Deformations according to the readings of the FBG sensors under 1,27 kN Параметр Кольцевое направление B009 B010 B011 B012 - Деформации по показаниям ВБР*10-3 -0,18 0,049 0,057 0,12 - Расчетные значения деформаций*10-3 0,001 283 0,00066 0,003 556 0,01397 - Отклонение, % 99,29 98,65 93,76 88,36 - Параметр Под углом 45° к осевому направлению Осевое направление B103 B104 B105 B111 B112 Деформации по показаниям ВБР*10-3 -0,81 -0,89 -0,15 -0,92 -0,78 Расчетные значения деформаций*10-3 -0,889 -0,8763 -0,08763 -0,8763 -0,889 Отклонение, % 8,89 1,54 41,58 4,75 12,26 Как видно из табл. 5, деформации, определяемые по показаниям датчиков В009, В010, В011, В012, В105, остались недостоверны и значительно отличаются от расчетных, поскольку деформация в кольцевом направлении при данном виде испытаний меньше погрешности интеррогатора. Показания датчиков В103, В104, В111 и В112, регистрирующих деформацию под углом 45° и осевую деформацию, удовлетворительно соответствуют расчетным значениями. Датчик В111 в данном случае показывает хорошее совпадение с расчётом; скачки в показаниях этого датчика связаны с высокой зашумленностью сигнала, вызванной плохой приклейкой. Однако регистрируемые значения деформаций и соответствующие рассчитываемые значения упругих напряжений в плоскости армирования слоев в сегменте шпангоута далеки от предельных. Для определения места разрушения и выявления критических компонент тензора напряжений использовались поля напряжений в слоях и критерий максимальных напряжений. Предельные значения напряжений приведены в табл. 3. Анализ показал, что наибольшие напряжения в сегменте шпангоута возникают в области изгиба слоев при переходе от П-образной полки на фланцы. При этом в плоскости армирования наибольшие значения растягивающих напряжений вдоль основы σ11 и вдоль утка σ22 возникают в первом слое на внутренней поверхности образца. При нагрузке 1,27 кН они составляют 210,9 и 53,4 МПа соответственно. Максимальные сжимающие напряжения вдоль основы σ11 и вдоль утка σ22 возникают в последнем слое на наружной поверхности образца и составляют -145,4 и -41,5 МПа. Наибольшие касательные напряжения в плоскости армирующих слоев τ12 не превышают по модулю 29,7 МПа. Сравнивая полученные значения напряжений с предельными значениями (см. табл. 3), можно отметить, что при нагрузке 1,27 КН, которая вызывает появление первой трещины, запас прочности на разрыв, сжатие и сдвиг в плоскости армирующих слоев составляет приблизительно 3,8, 5,5 и 5 соответственно. Анализ межслойных напряжений показал, что межслойные касательные напряжения при нагрузке 1,27 кН невелики, их максимальные значения τ13 = 16 МПа и τ23 = 4,2 МПа. При этом запас прочности сегмента шпангоута по критерию максимальных напряжений можно оценить как 4,8 и 18. Наиболее опасными для сегмента шпангоута при данной схеме нагружения являются межслойные отрывные напряжения σ33. Максимальные значения этих напряжений наблюдаются в центральных слоях, в отмеченных зонах скругления профиля и при расчетной нагрузке 1,27 кН составляют порядка 46,3 МПа (рис. 9). Рис. 9. Поле межслойных нормальных напряжений s33 [MПa] при нагрузке 1,27 кН Fig. 9. The field of interlayer normal stresses s33 [MPa] under 1,27 кN Именно эти напряжения определяют начало расслоения сегмента шпангоута в процессе механических испытаний, что приводит к скачкам на диаграмме деформирования. Следует отметить, что прочность на межслоевой сдвиг и отрыв слоистого композита в значительной степени зависят не только от материала связующего, но и от технологии изготовления и, в частности, давления прессования. При недостаточном или избыточном давлении прессования можно получить значительные отклонения этой характеристики. Такой же эффект может возникнуть и в зонах искривления или перегиба слоев в конструкции. Предельные значения межслойных отрывных напряжений имеют значительный разброс, указанный в табл. 3. Сравнивая значение межслоевых отрывных напряжений 46 МПа при критической нагрузке 1,27 кН, полученное в результате расчета, со средним значением соответствующего диапазона прочности слоистого углепластика (см. табл. 3) на поперечное растяжение 50 МПа, можно отметить, что отличие составляет приблизительно 8%. Таким образом, математическое моделирование лабораторных испытаний позволяет спрогнозировать характер разрушения и оценивать прочность сегмента шпангоута с достаточной точностью по критерию максимальных напряжений, используя известные значения предельных характеристик для слоистого композита. Разрушение сегмента композитного шпангоута происходит в виде расслоения, вызванного высокими межслоевыми отрывными напряжениями в зонах перегиба слоев. Показания оптоволоконных датчиков регистрирующих при этом напряжения в плоскости слоев в отдельных точках шпангоута могут являться индикаторами уровня критичности напряженного состояния. Используя значения этих датчиков при мониторинге состояния шпангоута в процессе испытаний и результаты расчета НДС шпангоута при схеме нагружения, можно косвенно оценивать значения межслойных отрывных напряжений и контролировать достижение критического состояния конструкции, при котором появятся межслоевые трещины. Кроме того, скачкообразное изменение показаний датчиков будет свидетельствовать о появлении межслоевых трещин в композитной конструкции. Заключение Проведенные исследования показали, что оптоволоконные датчики, расположенные в плоскости слоев композитной конструкции, позволяют регистрировать момент возникновения и развития межслойных разрушений за счет скачкообразного изменения деформации в плоскости слоя, связанного с локальными изменениями жесткости композита в окрестности расслоения. Сопоставляя показания оптоволоконных датчиков на основе ВБР и результаты математического моделирования, можно оценивать запас прочности конструкции при эксплуатационной нагрузке как по напряжениям в плоскости слоя, так и по межслойным напряжениям. Результаты исследования демонстрируют возможность косвенной оценки межслоевых напряжений и деформаций в композитной конструкции при механических испытаниях по показаниям оптоволоконных датчиков, расположенных в плоскости слоя, используя результаты расчета напряженно-деформированного состояния конструкции в упругой постановке с явным описанием схемы армирования. Представляет интерес проведение дальнейших исследований для оценки схемы расположения датчиков в слоистых композитных конструкциях, позволяющих более точно оценивать межслойные напряжения и регистрировать изменение продольной жесткости слоя при появлении расслоений. Полученные результаты будут использованы для разработки системы мониторинга и контроля НДС композитного силового шпангоута авиационного двигателя для стендовых и летных испытаний.About the authors
A N Anoshkin
Perm National Research Polytechnic University
V Yu Zuiko
Perm National Research Polytechnic University
K A Pelenev
Perm National Research Polytechnic University
P V Pisarev
Perm National Research Polytechnic University
G S Shipunov
Perm National Research Polytechnic University
References
- Anoshkin A.N., Zuiko V.Yu., Artemev V.V. Design and manufacture of aircraft engine frame from polymer composites depending on technological features // Conference on Mechanics of Composite Materials. - Riga, 2014.
- Numerical and experimental study of composite bulkhead partition strength with in-situ x-ray monitoring / N.N. Potrahov, A.N. Anoshkin, V.Yu. Zuiko, V.M. Osokin, P.V. Pisarev, K.A. Pele¬nev // PNRPU Mechanics Bulletin. - 2017. - No. 1. - Р. 118-133. doi: 10.15593/perm.mech/2017.1.08
- Staszewski W., Tomlinson G., Boller C. Health monitoring of aerospace structures smart sensor technologies and signal pro-cessing. - 2004.
- Farrar C.R., Worden K. An introduction to structural health monitoring // Philos. Trans. A. Math. Phys. Eng. Sci. - 2007. - Vol. 365. - No. 1851. - Р. 303-315, doi: 10.3390/s131216551
- Diamanti K., Soutis C. Structural health monitoring tech-niques for aircraft composite structures // Prog. Aerosp. Sci. - 2010. - Vol. 46. - No. 8. - Р. 342-352.
- Advanced design for lightweight structures: Review and prospects / D.F.O. Braga, S.M.O. Tavares, L.F.M. Da Silva, P.M.G.P. Mo¬reira, P.M.S.T. De Castro // Prog. Aerosp. Sci. 2014. - Vol. 69. - Р. 29-39. doi: 10.1016/j.paerosci.2014.03.003
- SMART Layer and SMART Suitcase for structural health monitoring applications / M. Lin, X. Qing, A. Kumar, S.J. Beard // Proc. SPIE. - 2001. - Vol. 4332. - No. 408. - Р. 98-106. doi: 10.1117/12.429646
- Fan Y., Kahrizi M. Characterization of a FBG strain gage array embedded in composite structure // Sensors Actuators, A Phys. - 2005. - Vol. 121. - No. 2. - Р. 297-305. doi: 10.1016/j.sna.2005.01.021
- Application of optical fiber as strain gauges in polymer composite materials / E.N. Kablov [et al.] // Polym. Sci. Ser. D. - 2011. - Vol. 4. - No. 3. - P. 246-251. doi: 10.1134/S1995421211030063
- Structural health monitoring of composite aircraft struc-tures using fiber bragg grating sensors / N. Gupta [et al.] // J. Indian Inst. Sci. - 2013. - Vol. 93. - P. 735-750.
- A survey of health monitoring systems for wind turbines / M.L. Wymore, J.E. Van Dam, H. Ceylan, D. Qiao // Renew. Sus-tain. Energy Rev. - 2015. - Vol. 52. - No. 1069283. - P. 976-990. doi: 10.1016/j.rser.2015.07.110
- Sierra-Pérez J., Torres-Arredondo M.A., Güemes A. Damage and nonlinearities detection in wind turbine blades based on strain field pattern recognition. FBGs, OBR and strain gauges comparison // Compos. Struct. - 2016. - Vol. 135. - P. 156-166. doi: 10.1016/j.compstruct.2015.08.137
- Application of FBG sensors for geotechnical health monitoring, a review of sensor design, implementation methods and packaging techniques / C.Y. Hong, Y.F. Zhang, M.X. Zhang, L.M.G. Leung, L.Q. Liu // Sensors and Actuators, A: Physical. - 2016. - Vol. 244. - P. 184-197. doi: 10.1016/j.sna.2016.04.033
- Hafizi Z.M., Epaarachchi J., Lau K.T. Impact location determination on thin laminated composite plates using an NIR-FBG sensor system // Meas. J. Int. Meas. Confed. - 2015. - Vol. 61. - P. 51-57. doi: 10.1088/1757-899X/100/1/012056
- Beukema R.P. Embedding technologies of FBG sensors in composites: Technologies, applications and practical use // Proc. 6th Eur. Work. - Struct. Heal. Monit. 2012, EWSHM-2012. - 2012. - Vol. 1. - P. 341-348.
- A novel time-division multiplexing fiber Bragg grating sensor interrogator for structural health monitoring / Y. Dai, Y. Liu, J. Leng, G. Deng, A. Asundi // Opt. Lasers Eng. - 2009. - Vol. 47. - No. 10. - P. 1028-1033. doi: 10.1016/j.optlaseng.2009.05.012
- Use of FBG sensors for SHM in aerospace structures / G.C. Kahandawa, J. Epaarachchi, H. Wang, K.T. Lau // Photonic Sensors. - 2012. - Vol. 2. - No. 3. - P. 203-214. doi: 10.1007/s13320-012-0065-4
- On the improvement of strain measurements with FBG sensors embedded in unidirectional composites / G. Pereira, C. Frias, H. Faria, O. Frazão, A.T. Marques // Polym. Test. - 2013. - Vol. 32. - No. 1. - P. 99-105.
- Fibre Bragg gratings in structural health monitoring - Present status and applications / M. Majumder, T.K. Gango-padhyay, A.K. Chakraborty, K. Dasgupta, D.K. Bhattacharya // Sensors and Actuators a-Physical. - 2008. - Vol. 147. - No. 1. - P. 150-164. doi: 10.1016/j.sna.2008.04.008
- For an issue of embedded optical fibre in CFRP and strain measurement with fibre Bragg gratings sensors / M. Makhsi¬dov, V. Fedotov, M. Shiyonok, A. Zuev // J. Compos. Mech. Des. - 2014. - Vol. 20. - No. 4. - P. 568-574.
- Güemes A., Menendez J.M. Fiber-Optic Sensors // Structural Health Monitoring. - 2010. - P. 225-285. doi: 10.3390/s150715494
- Применение оптического волокна в качестве дат¬чиков деформаций в полимерных композиционных материа¬лах / E. Каб¬лов [и др.] // Все материалы. Энциклопедический справочник. - 2010. - № 3. - С. 10-15.
- Медведков О., Королев И., Васильев С. Запись воло-конных Бреговских решёток в схеме с интерферометром Лой¬да и моделирование их спектральных свойств / НЦВО РАН. - М., 2004. - № 6. - С. 1-46.
- Применение оптоволоконных технологий при создании встроенных систем самодиагностики авиационных кон¬струк¬ций / А. Сережнов, A. Кузнецов, А. Лукьянов, А. Бра-гин // Научный Вестник НГТУ. - 2016. - Т. 64, № 3. - С. 95-105.
- Технологии и задачи механики композиционных материалов для создания лопатки спрямляющего аппарата авиационного двигателя / А.Н. Аношкин, В.Ю. Зуйко, Г.С. Ши-пунов, А.А. Третьяков // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 4. - С. 5-44. doi: 10.15593/perm.mech/2014.4.01
- 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.
- 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
- Березин А.В., Козинкина А.И. Особенности диагнос-тики повреждений и оценки прочности композита // Механика композиционных материалов и конструкций. - 1999. - № 1. - С. 99-122.
- Аношкин А.Н., Ташкинов А.А. Прогнозирование не-су¬щей способности композитных фланцев корпусных деталей авиадвигателей / Перм. гос. техн. ун-т. - Пермь, 1998. - С. 101.
- Лукина Н.Ф. Клеевые препреги на основе тканей Porcher - перспективные материалы для деталей и агрегатов из ПКМ [Электронный ресурс] // Труды ВИАМ. - 2014. - № 6. - URL: http://viam-works.ru/ru/articles?art_id=677 (дата обращения: 03.11.16). doi: 10.18577/2307-6046-2014-0-6-10-10
- Исследование НДС и оценка прочности композитной лопатки спрямляющего аппарата авиационного двигателя / М.А. Гринев, А.Н. Аношкин, П.В. Писарев, В.Ю. Зуйко, Г.С. Ши¬пунов // Вестник Пермского национального исследо-вательского политехнического университета. Механика. - 2015. - № 4. - С. 293-307. doi: 10.15593/perm.mech/2015.4.17
- Компьютерное моделирование механического пове¬де-ния композитной лопатки спрямляющего аппарата авиа¬ци-онного двигателя / М.А. Гринев, А.Н. Аношкин, П.В. Писарев, В.Ю. Зуйко, Г.С. Шипунов // Вестник Пермского националь-ного исследовательского политехнического университета. Механика. - 2015. - № 3. - С. 38-51. doi: 10.15593/perm.mech/2015.3.04
- Experimental study of non-uniform strains in composite structures with embedded fiber Bragg grating sensors / A.N. Anosh¬kin, V.P. Matveenko [et al.] // The Conference on Smart Structures and Materials SMART 2017, 6th International Conference on Smart Materials and Nanotechnology in Engineering SMN2017, 5-8 June 2017. - Madrid, Spain, 2017.