CC BY
116
2. Белан, Н. В. Стационарные плазменные двигатели : учеб. пособие для студентов вузов / Н. В. Белан, В. П. Ким, А. И. Оранский, В. Б. Тихонов ; Харьков. авиац. ин-т. Харьков, 1989.
3. Гришин, С. Д. Электрические ракетные двигатели космических аппаратов / С. Д. Гришин, Л. В. Лесков, Н. П. Козлов. М. : Машиностроение, 1975.
4. Пат. РФ № 2110131. МГД-способ преобразования тепловой энергии в электрическую замкнутого цикла / В. В. Данилов, В. С. Славин. Б. И. 1998. N° 2.
5. Пат. РФ № 2162958. Способ ускорения потока рабочего тела в канале ракетного двигателя / В. С. Славин,
В. В. Данилов, М. В. Краев. Б. И. 2001. № 4.
6. Славин, В. С. Космическая транспортная система на основе МГД-генератора замкнутого цикла и ускорителя газоплазменного потока / В. С. Славин, В. В. Данилов, М. В. Краев // Перспективы МГД и плазменных технологий в аэрокосмических приложениях : сб. трудов. М. : И ВТ РАН, 1999. С. 32-36.
V M. Kraev, V S. Slavin, M. V Kraev
ACCELERATION OF THE GAS STREAM IN THE CHANNEL OF THE ELECTROJET ENGINE
The opportunity of the further perfection of electrojet engines is examined due to gas stream acceleration in its channel due to electrodynamic force of Lorentz in non-uniform gas-plasma stream. The calculation results of the engine key parameters with draft 1 400 Н are given.
Keywords: the electro jet engine, a gas stream, acceleration.
УДК629.7.192:656.735.3
Н. А. Махутов, В. В. Москвичев, А. С. Серегин, М. Н. Степнов
ОЦЕНКА ПРЕДЕЛЬНЫХ СОСТОЯНИЙ ЛОПАСТЕЙ ВОЗДУШНЫХ ВИНТОВ САМОЛЕТОВ
Рассмотрены основные типы предельных состояний для лопастей воздушных винтов самолетов. Выполнены расчетно-экспериментальные оценки остаточной прочности и ресурса лопастей при статическом и циклическом нагружениях. Установлены критические и допускаемые размеры трещин, получены расчетные кривые усталости, позволившие сделать выводы об изменении коэффициентов запаса по долговечности при наличии исходных дефектов.
Ключевые слова: воздушные винты, усталостная прочность, конструкции самолетов.
Лопасти воздушных винтов относятся к наиболее нагруженным элементам конструкций самолетов. Основные причины, приводящие к отказу лопастей, связаны с нестационарным циклическим нагружением от совместного действия центробежной силы и изгибающего момента в процессе эксплуатации, что в сочетании с естественной статистической природой характеристик усталостной прочности приводит к разрушению лопастей при числах циклов N = 105.. ,106 [1]. В этих условиях заданный срок службы лопасти определяется по результатам натурных испытаний лопастей на основе статистических методов расчета. Однако в силу значительного разброса характеристик усталостной прочности наблюдаются случаи отказа лопастей до полной выработки установленного ресурса [2], тогда как большее количество лопастей снимается с эксплуатации, не исчерпав своей несущей способности.
В процессе циклического нагружения происходит формирование магистральной усталостной трещины от технологических и эксплуатационных дефектов, конфигурация которой определяется схемой приложения действую-
щих нагрузок и геометрией лопасти в месте расположения дефектов. Наиболее вероятным дефектом являются поверхностные полуэллиптические трещины в области прикомлевых сечений, возникающие под действием центробежной силы и изгибающего момента. Нельзя также исключать возможности развития боковых трещин от передней и задней кромок, а также кольцевых трещин в местах крепления лопасти от центробежных усилий. При наличии усталостных трещин наиболее опасными становятся динамические перегрузки, возникающие в результате скоростных маневров самолета и турбулентных порывов ветра. В большинстве случаев именно с ними связано окончательное разрушение лопасти при длинах трещин меньше рассчитанных по результатам циклических испытаний.
Достижение предельного состояния в лопасти может осуществляться за счет:
- образования и развития трещин многоцикловой усталости от вибрационных нагрузок с большим числом циклов нагружения (N = 106);
- образования и развития трещин малоцикловой усталости от повторных нагрузок при взлетах и посадках;
- образование разрушения типа хрупкого или квази-хрупкого на определенной стадии развития усталостных трещин в результате динамических перегрузок.
Исходя из этого, на основе экспериментальных и расчетных методов механики разрушения, была разработана программа исследования предельных состояний, возникающих в лопастях воздушных винтов самолетов. Программа включала в себя следующее:
- анализ эксплуатационной нагруженности лопасти;
- проведение испытаний образцов, вырезанных из различных участков лопастей с целью определения характеристик трещиностойкости при статическом, циклическом и ударном нагружениях алюминиевого сплава Д1, из которого изготовляется лопасть;
- обоснование расчетных схем и методов расчета несущей способности тел с трещинами, применительно к лопасти воздушного винта;
- определение критических и допускаемых размеров трещин и чисел циклов нагружения по результатам испытаний образцов и их сопоставление с данными натурных испытаний лопастей.
Схема разрезки лопастей предполагала изготовление различных типов образцов с трещинами, которые наиболее полно отражали схему нагружения данного участка лопасти в соответствии с действующими усилиями. Так из комлевой части лопасти, нагруженной осевыми инерционными усилиями, изготавливали цилиндрические образцы с кольцевой трещиной для испытания на осевое растяжение. Были предусмотрены испытания на осевое растяжение и изгиб образцов с боковыми и полуэллип-тическими трещинами при статическом и циклическом нагружениях для оценки влияния растягивающих напряжений от инерционных усилий и изгибающих моментов. Совместное влияние этих нагрузок на распространение усталостных трещин исследовалось на образцах внецент-ренного растяжения.
Особое внимание было уделено оценке предельного состояния комлевой части лопасти, в связи с имевшимися случаями отрыва лопастей во время набора высоты [2]. Комлевая часть имеет осесимметричную форму, по поверхности которой для крепления лопасти наносится круглая резьба, при этом по дну канавок возможно возникновение и развитие кольцевой трещины в процессе эксплуатации.
При проведении расчетов были приняты средние значения механических характеристик сплава Д1 - а02 = 300 МПа и а в = 500 МПа. Эксплуатационное максимальное растягивающие напряжение от действующих усилий в расчетном сечении составляет стэ = 80 МПа. С учетом коэффициента безопасности Е = 1,5 и коэффициента максимальной перегрузки пшах = 3,0 расчетное напряжение составит ар = Е • пшах -стэ = 360 МПа. Тогда коэффициент запаса прочности по пределу прочности будет п = ав/ ар = 1,4. Примем значение п = 1,5.
Оценка остаточной прочности лопасти в случае возникновения кольцевых трещин может быть выполнена на основе концепции предела трещиностойкости [3] с использованием результатов экспериментальных исследований характеристик трещиностойкости сплава Д1, полученных на цилиндрических образцах с кольцевой тре-
щиной в широком диапазоне длин трещин I и диаметров образцов б [4; 5]. Результаты таких исследований представлены на рис. 1 в виде зависимостей критических значений коэффициентов интенсивности напряжений К*, определенных по максимальной нагрузке на диаграмме «нагрузка-смещение берегов трещины» от длины трещины I для значений Б = 90; 70; 35; 17 мм.
Расчет предельного состояния согласно концепции предела трещиностойкости ведется по критерию [3]
К* = 1С, (1)
где 1С - предел трещиностойкости, представляющий собой непрерывную совокупность значений предельных коэффициентов интенсивности напряжений КС* для всего диапазона длин трещин (включая I = 0), и записанного в виде функции от обратного коэффициента запаса прочности по пределу прочности:
Vп = стС Ь в, (2)
где стС - разрушающее напряжение в брутто-сечении.
Тогда зависимость К*с -<з6с/ств, согласно определения, есть предел трещиностойкости и она может быть записана в виде 1С -1/п . Расчет значений К* проводили по формуле [6]:
К =°С— Е ), (3)
где ё - диаметр нетто-сечения; Е (ё/Б) - поправочная
0 5 10 15 20 25 30 I, мм
Рис. 1. Зависимость предельных значений коэффициентов интенсивности напряжений от длины трещины для цилиндрических образцов из сплава Д1: о - Б = 90 мм; □ - 70 мм; Д - 35 мм; ◊ - 17 мм
Формулу (3) с учетом соотношений (1) и (2) можно записать через величину ав и п
= * = ^2%Б 11 - ё/Б Е ( , .
!с = Кс =ав 4п ^ё/о)3 Е (ё/В). (4)
Из выражения (4) следует, что коэффициент интенсивности напряжений зависит от коэффициента запаса прочности п и имеет вид
К* = К** (п) = !с (п) . (5)
Варьируя величину п можно получить серию кривых К* - I (I = О(1 -ё/О)/2) для заданных значений диаметра О (рис. 2).
Для определения допускаемой длины трещины в выражениях (1) и (5) (рис. 2, а) необходимо ввести коэффициент запаса прочности по пределу трещиностойкости т, показывающий, во сколько раз необходимо увеличить значения К* за счет увеличения длины трещины при неизменных нагрузках и запасах прочности п , чтобы они достигли предельных величин, т. е.
1С (п)/т=1с (п )/т=кс* (п). (6)
По этому выражению определяются значения Кс (п) для разных запасов т, используя экспериментальные данные, представленные на рис. 1. Соответствующие результаты расчета по формулам (4) и (6) для коэффициентов т = 1,0; 1,5; 2,0; 2,5 приведены на рис. 2, б. Критическая длина трещины 1с при конкретном значении коэффициента п определяется графически (рис. 2) при т = 1. В нашем случае п = 1,5, и тогда критическая длина трещины по результатам эксперимента при т = 1 составит
1с = 6,15 мм. Коэффициент запаса по длине трещины рекомендуется принимать т = 2,0... 4,0. Зная, что т = 2,5, получим допускаемую длину трещины [I ] = 2,6 мм.
Значения предела трещиностойкости могут быть определены расчетным путем. Если предположить, что независимо от исходной длины трещины разрушающее напряжение по нетто-сечению аН = а в (согласно I теории прочности) и имея в виду, что ст^1 =<56с (ё/О)2 , получим связь между относительной длиной трещины и коэффициентом запаса п:
1=# <7'
О V п
Справедливость данного предположения подтверждается экспериментальными результатами испытаний гладких образцов и образцов с трещинами при различных значениях ё/О (рис. 3). Подставляя формулу (7) в (4) получаем зависимость
к;,
Н/мм3/2
2000
1600
1200
800
400
0
п = 1 1 /1,5/ У 2 // / -2,5
г/ т
1/
\ 1
| 1 4-1 1 1
Ш У ТЧ ''
(8)
Рис. 3. Зависимость обратного коэффициента запаса прочности 1/ п = абС /а в от относительной длины трещины: о - О = 90 мм; □ - 70 мм; А - 35 мм; ◊ - 17 мм
При варьировании п могут быть построены расчетные кривые 1С /т -1/ п (рис. 2, б). При т = 1 расчетная кривая с достаточной степенью точности соответствует экспериментальным точкам, полученным при испытаниях образцов с О = 70 мм. Критическая длина трещины при определении по расчетной зависимости 1с равна 1с = 6,54 мм.
Значения К* зависят от набольшего диаметра образцов (рис. 1), и это определяет существование аналогичной зависимости для значений 1с (рис. 4). Учет этого факта возможен путем введения относительного предела трещиностойкости в форме (рис. 4, а):
1С
4п
1(7^)‘
(9)
1с/т
т=1
-~а /У V- \ |Я\ I
/ / □ / -- -'■и*- . Л5 ■О-
// и Р ;,0 С а *. "™П. ч \
—-НИ-/ & \
2,5
[I] и
10
а
15 20 /,
мм 0 0,2 0,4 0,6 а С/ав=1/п
б
Рис. 2. Зависимость расчетных значений Кс (уравнение 4) при различных коэффициентах п от длины трещины I (а) и предела трещиностойкости 1с (б) при различных коэффициентах т для серии образцов О = 70 мм: □ - экспериментальные точки; I - расчет по уравнению (8) при т = 1
Расчетная кривая по уравнению (9) удовлетворительно описывает (рис. 4, б) результаты испытаний образцов с В = 90; 70 мм и дает заниженные значения по сравнению с результатами испытаний образцов с диаметром В = 35; 17 мм. Несоответствие этих результатов можно объяснить тем, что образцы малого диаметра разрушаются при развитых пластических деформациях, при этом стН >ств, т. е. нарушается принятое при выводе формулы (9) условие стН =ств.
Из рис. 4, б следует, что критические длины трещин могут определяться с удовлетворительной точностью по уравнению (8), начиная с диаметра В = 70 мм и выше, т. е. отсюда следует вывод, что размеры критических 1С и допускаемых [/] длин трещин, определенные на образцах с В =70 мм (см. рис. 2), будут соответствовать критическим и допускаемым длинам трещин, возникающих в натурном сечении концевой части лопасти.
Для оценки предельных состояний лопасти на стадии развития дефекта при циклическом нагружении в качестве исходной информации необходимо иметь диаграммы циклического разрушения, устанавливающие связь между скоростью роста трещины и величиной коэффициента инстенсивности напряжений. Такие диаграммы были получены по результатам усталостных испытаний образцов, имеющих боковые и полуэллиптические трещины и вырезанных из различных участков лопасти, на трехточечный изгиб, осевое и внецентренное растяжение. Диаграммы циклического разрушения, представленные для отдельных испытаний на рис. 5, удовлетворительно описываются степенным уравнением Пэриса
— = с (ак )п
(10)
где С, п - постоянные материала и условий испытания; АК - размах коэффициента интенсивности напряжений.
Возможность разрушения лопасти на определенной стадии развития усталостной трещины при наличии динамических перегрузок требует определения характеристик трещиностойкости материала при динамическом
нагружении. Были проведены испытания образцов с боковой трещиной двух типоразмеров по ГОСТ 9454-78 на трехточечный изгиб в широком диапазоне скоростей нагружения, характеризующихся скоростью изменения коэффициента интенсивности напряжений К1 = 10-1...105 Н-мм-3/2/мин. Кроме того, были проведены ударные испытания (К @ 105 Н-мм-32/мин) при пониженных температурах (до -180 ° С). Как и следовало ожидать, учитывая малую чувствительность алюминиевых сплавов к скорости деформирования, полученные значения Кдс динамического коэффициента интенсивности напряжений практически не зависят от скорости нагружения. Совместное действие динамического нагружения и пониженных температур, как одних из основных охрупчивающих факторов, также мало сказывается на сопротивлении сплава Д1 развитию трещин. В процессе эксплуатации для лопасти следует ожидать изменение величины К1 от 103 до 105 Н-мм-3/2/мин. Тогда, полученные в этом диапазоне значения Кдс , средняя величина которых составила Кдс = 500 Н/мм32, могут быть использованы в качестве ориентировочных расчетных характеристик для консервативной оценки остаточной прочности лопасти при наличии дефектов в условиях возможных ударных перегрузок.
Результаты циклических и динамических испытаний образцов были использованы при рассмотрении кинетики роста полуэллиптической трещины, являющейся наиболее вероятным дефектом для лопасти. Исходные размеры трещины назначались с учетом возможностей дефектоскопического контроля, позволяющего выявлять поверхностные дефекты длиной 1,5 .6,0 мм, в зависимости от типа искательных головок дефектоскопа [7]. Начальная глубина трещины была принята 10 = 1,0 и 3,0 мм при соотношении полуосей 10/2а = 0,2. Расчет числа циклов Ы, на основе выражения (10), при увеличении начального дефекта с размером 10 производился по формуле
2 ( . \
N = -
(п - 2) СМ”12 Аст'
1
1
<п-2)/2 ] (п - 2)/2
/
(11)
а б
Рис. 4. Зависимость предела трещиностойкости 1с для образцов различных диаметров (а) и относительного предела трещиностойкости 1^ (стг 43) (б): о - 3> = 90 мм; □ - 70 мм; Д - 35 мм; ◊ - 17 мм
I
С
где ^ - текущая длина усталостной трещины; М- параметр формы трещины в выражении, определяющем коэффициент интенсивности напряжений [8].
ДК = ДаТМ7. (12)
50
100
150
1
(п - 2)
V (п-2)
N• С• Мп2 -Аст'
при которых возможно разрушение лопасти при динамических перегрузках. Рассматривая связь между напряжениями и критическими числами циклов можно построить кривую усталости (рис. 7), соответствующую предельному состоянию разрушения лопасти, имеющей исходный дефект 10 = 3 мм (кривая 1). Разрушение лопасти только под действием циклического нагружения без учета динамических перегрузок описывается кривой усталости 2, а для лопасти, имеющей исходный дефект 10 = 1 мм
- кривой 3. Здесь же представлена кривая усталости лопастей воздушных винтов без исходной дефектности (кривая 4), построенная по стадии окончательного разрушения в результате проведения натурных стендовых испытаний [1].
250 д К, кгс/мм32
Рис. 5. Диаграммы циклического разрушения образцов, вырезанных из наиболее нагруженных участков лопасти: о • - трехточечный изгиб образцов с боковой трещиной;
□ - внецентренное растяжение. о - горбушка;
• - рабочая поверхность лопасти
На рис. 6, а для разных величин размаха напряжений Да = 60; 80; 120 МПа при М = 2,9 показаны, рассчитанные по уравнению (11), кривые роста трещины ^ в зависимости от числа циклов N. При расчете использованы значения коэффициентов С = 3, 110-14 и п = 5,0, полученные при обработке результатов циклических испытаний образцов (см. рис. 5). Если из (11) выразить длину трещины и подставить в (12), как это сделано в работе [9], то получим связь между коэффициентом интенсивности напряжений К1 и числом циклов нагружения (рис. 6):
Да М12
(13)
Рис. 6. Зависимости длины усталостной трещины (а) и коэффициента интенсивности напряжений (б) от числа циклов нагружения для заданных значений Даа
Анализ кривых усталости позволяет сделать определенные выводы. При наличии в лопасти повышенной исходной макродефектности (поверхностная трещина глубиной до 3 мм) по данным расчетно-экспериментальной оценки наблюдается снижение долговечности в 6.8 раз. Снижение исходной дефектности от 10 = 3 мм до 10 = 1 мм при том же уровне циклических напряжений приводит к увеличению долговечности в 3.6 раз. При допущенном уровне повреждения лопасти трещинами (/0 = 3 мм) снижение запаса по отношению к стадии разрушения лопасти, не содержащей исходной дефектности (рис. 7), по напряжениям составляет 1,5___1,7 раза, а по числу циклов
нагружения с учетом динамических перегрузок снижение запаса увеличилось до 10 раз.
/(п-2)/2 2
С использованием этого выражения для исходной длины трещины /0 = 3,0 мм рассчитаны зависимости К1 - N (рис. 6, б) для тех же самых значений Да . Данные зависимости отражают изменение коэффициента интенсивности напряжений за счет увеличения размера дефекта с ростом числа циклов нагружения. Учитывая возможность динамических перегрузок лопасти уровень значений К1 (рис. 6, б) должен быть ограничен минимальной величиной Кдс = 500 Н/мм32. Исходя из этого можно определить для различных амплитуд циклических напряжений критические размеры трещин и числа циклов нагружения,
а ,
МПа
130
90
50
% Ч X \ >
N \ N V 2\ ч г 4
С учетом *■> , перегрузки 1
10б 2
6 10' 2
б 10
К,цикл
Рис. 7. Кривые усталости лопастей: 1, 2, 3 - расчет; 4 - эксперимент; 2 - /0 = 3,0 мм; 3 - /0 = 1,0 мм
Таким образом, расчетно-экспериментальный анализ результатов статических, циклических и динамических испытаний позволил рассмотреть вопрос о допускаемых размерах дефектов в лопасти и критических числах циклов нагружения. Обращает на себя внимание резкое увеличение коэффициентов интенсивности напряжений и повышение скоростей развития усталостных трещин на определенной стадии циклического нагружения, что указывает на необходимость более тщательного дефектоскопическо -го контроля трещин в лопастях не по количеству часов налета, а по числу взлетов и посадок. Назначение сроков периодического контроля, обеспечивающих безопасность полетов при эксплуатации лопастей по состоянию, требует дополнительных статистических данных об эксплуатационной нагруженности, законах развития усталостных трещин и характеристиках трещиностойкости материала.
Библиографический список
1. Степнов, М. Н. Сопротивление лопастей воздушных винтов самолетов Проблемы прочности / М. Н. Степнов, А. С. Серегин, О. В. Леонова, Ю. Л. Сухоросов, Е. Н. Куликов, Э. В. Дунаев. 1977. N° 5.
2. Савкина, И. А. Прочность и долговечность авиационных конструкций. Обзор по материалам иностран-
ной печати / И. А. Савкина, В. В. Кирин. М. : ЦНТИ ГА, 1976.
3. Партон, В. З. Механика упругопластического разрушения / В. З. Партон, Е. М. Морозов. М. : Наука, 1974.
4. Махутов, Н. А. Характеристики механики разрушения в связи с размерами трещин / Н. А. Махутов,
B. В. Москвичев // Машиноведение. 1983. № 6. С. 67-71.
5. Трещиностойкость и механические свойства конструкционных материалов технических систем / В. В. Москвичев, Н. А. Махутов, А. П. Черняев и др. Новосибирск : Наука, 2002.
6. Андрейкив, А. Е. К теории предельного равновесия цилиндрического образца с внешними кольцевыми трещинами / А. Е. Андрейкив, В. В. Панасюк, Н. Н. Панько // ФХММ. 1974. № 3. С. 29-39.
7. Дефектоскопия деталей при эксплуатации авиационной техники / под ред. П. И. Беды. М. : Воениздат, 1978.
8. Браун, У. Испытания высокопрочных металлических материалов на вязкость разрушения при плоской деформации / У Браун, Дж. Сроули ; пер. с англ. под ред. Б. А. Дроздовского. М. : Мир, 1972.
9. Васильченко, Г. С. Расчет на прочность массивных конструкций, содержащих дефекты / Г. С. Васильченко, Е. М. Морозов // Вестник машиностроения. 1977. № 3.
C. 72-77.
N. A. Makhutov, V. V. Moskvichev, A. S. Seregin, M. N. Stepanov ESTIMATION OF PLANE’S BLADE MARGINAL STATES
The main types of plane’s blade marginal states are considered. Experimental calculated estimation of the retained strength and blade life time at static and cyclic loading are calculated. The ultimate and permissible dimension of cracks are determined. The calculated exhaust curves are received. It let make a conclusion about safety life time factor with the initial malfunction availability.
Keywords: exhaust safety, blades, constructions, planes.
УДК 531:519.6
Р. Н. Хамитов
УГЛОВЫЕ КОЛЕБАНИЯ АМОРТИЗИРУЕМОГО ОБЪЕКТА С ПАССИВНОЙ СИСТЕМОЙ ДЕМПФИРОВАНИЯ
Исследуются угловые движения амортизируемого объекта с упругим подвесом и пассивным демпфированием. Рассмотрены приближенные решения для нестационарного сейсмовоздействия и приводятся результаты численного моделирования переходных процессов с помощью прикладной программы Matlab.
Ключевые слова: угловые движения, пассивная система демпфирования.
Особое место при анализе динамики амортизируе- ния и при кинематическом гармоническом возмуще-мых объектов (АО) занимают угловые движения. Круп- нии со стороны основания цель защиты объекта может ногабаритный объект можно рассматривать в первом заключаться в уменьшении амплитуды абсолютного приближении как абсолютно твердое тело на упругом ускорения (перегрузки) объекта. Степень реализации подвесе, а возмущение, действующее на объект со сто- цели характеризуется безразмерным коэффициентом роны основания, как кинематическое возмущение. В эффективности К:
этом случае, как это следует из справочника [1], для простейшего одномерного случая прямолинейного движе-
К = a / (w2 х )
