Содержание

Стр.
|ВВЕДЕНИЕ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | |
|. . . . . . . . . . . . . . . | |
|1. ОПИСАНИЕ КОНСТРУКЦИИ МОТОГОНДОЛЫ . . . . . . . . . . . . . . . | |
|2. СИЛОВОЙ РАСЧЕТ ВОЗДУХОЗАБОРНИКА . . . . . . . . . . . . . . . . . | |
|. | |
| 2.1. |Исходные данные для силового расчета . . . . . . . . | |
| |. . . . | |
| 2.2 |Распределение расчетных аэродинамических нагрузок по | |
| |длине воздухозаборника . . . . . . . . . . . . . . . | |
| |. . . . . . | |
| 2.3. |Распределение нагрузок по длине и по сечениям | |
| |воздухозаборника . . . . . . . . . . . . . . . . . . .| |
| |. . . . . . . . . . | |
| 2.4. |Распределение аэродинамических нагрузок по внутренней | |
| |поверхности воздухозаборника . . . . . . . . . . | |
| 2.5. |Определение равнодействующей по сечениям | |
| |воздухозаборника от внешних и внутренних | |
| |аэродинамических нагрузок . . . . . . . . . . . . . . | |
| |. . . . . . . | |
| 2.6. |Нагрузки на болты крепления воздухозаборника к | |
| |проставке . . . . . . . . . . . . . . . . . . . . . . | |
| |. . . . . . . . . . . . | |
| 2.7. |Проверка прочности воздухозаборника самолета . . . . .| |
| |. | |
| 2.8. |Автоматизация расчета аэродинамических нагрузок | |
| |воздухозаборника . . . . . . . . . . . . . . . . . . .| |
| |. . . . . . . . . . | |
|3. Технологический процесс изготовления воздухозаборника канала | |
|сотовой звукопоглощающей конструкции . . . . . . . . . . . . . . . . | |
|. . . . | |
|3.1. Технологичность конструкции воздухозаборника . . . . . . . . . .| |
|. . . . . | |
|3.2. Применяемые материалы и оборудование . . . . . . . . . . . . . .| |
|. . . . . . | |
|3.3. Технологический процесс сборки обшивок и элементов каркаса | |
|3.4. Использование в конструкции воздухозаборника композиционных | |
|материалов . . . . . . . . . . . . . . . . . . . . . . . . . . . . .| |
|. . . . . . . . . . . . . . . . | |
| 3.4.1|Методы получения ПКМ . . . . . . . . . . . . . . . . .| |
| |. . . . . . | |
|4. ОХРАНА ТРУДА И ОКРУЖАЮЩЕЙ СРЕДЫ . . . . . . . . . . . . . . . . . | |
|5. ЭКОНОМИКА И ОРГАНИЗАЦИЯ ПРОИЗВОДСТВА . . . . . . . . . . . . | |
|ЛИТЕРАТУРА | |
|ПРИЛОЖЕНИЕ | |
| | |

ВВЕДЕНИЕ

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

Они служат для торможения потока воздуха перед поступлением его в двигатель, а основными требованиями, предъявляемыми к входным устройствам, являются:

– обеспечение высоких значений коэффициента сохранения полного давления;

– создание равномерного потока на входе в двигатель или желаемой
(допустимой) неравномерности;

– минимальное аэродинамическое сопротивление;

– обеспечение устойчивой и эффективной работы во всем требуемом диапазоне режимов полета и режимов работы двигателя.

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

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

Положение мотоустановок на самолете

Рис. 1

– работы гондол двигателей и их пилонов как дополнительного горизонтального оперения;

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

– двигатели удалены от пассажирской кабины и от топливных баков; ж) повысить эксплуатационные характеристики силовой установки и всего самолета в целом за счет:

– обеспечения возможности замены целиком всей гондолы вместе с двигателем;

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

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

1) дозвуковые – для дозвуковых летательных аппаратов;

2) сверхзвуковые – для сверхзвуковых летательных аппаратов.

К дозвуковому диффузору ТРД относится не только сам внутренний канал, по которому воздух поступает к двигателю, но и примыкающая к нему входная часть – заборник воздуха. Заборник должен иметь плавное очертание входных кромок, что необходимо для предотвращения срыва потока на входе.

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

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

Существует три основных способа достижения этой цели: применение малошумных двигателей, более совершенные приемы эксплуатации самолетов и двигателей и рациональная установка двигателей на самолете.

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

Наибольшее влияние на уровень шума оказывает скорость истечение газа, поэтому действенным способом снижения шума является переход в пассажирской авиации от ТРД к двухконтурным двигателям, шум реактивной струи которых меньше из-за существенно меньшей ее скорости. Однако главным источником шума у ДТРД стал вентилятор. В настоящее время разработаны следующие основные способы снижения шума одноступенчатого вентилятора: отказ от ВНА вентилятора, пониженная окружная скорость рабочего колеса, оптимальное соотношение чисел лопаток выходного направляющего аппарата и рабочего колеса, увеличенное расстояние между этими рядами лопаток. Следует отметить, что, хотя применение турбовентиляторов с высокой частотой вращения позволяет снизить массу двигателя, требование по уровню шума заставляет ограничивать частоту вращения значениями, соответствующими окружным скоростям вентиляторов 400–450 м/с. Кроме того, рассматриваются другие предложения по снижению шума вентилятора одним из которых является способ снижения шума в процессе распространения его из воздухозаборника и выходного устройства. Этот способ включает облицовку стенок проточной части звукопоглощающими конструкциями (ЗПК). Пример применения таких конструкции в мотогондоле двигателя RB.211 для самолета L-1011 показан на рис. 2.
Применение ЗПК важно и тем, что при этом в конструкцию двигателя никаких изменений не вносится.

Акустически обработанная мотогондола двигателя пассажирского самолета

[pic] а – мотогондола с ЗПК; б – многослойная звукопоглощающая конструкция;

1 – перфорированная обечайка; 2 – сотовый заполнитель; 3 – опорная поверхность.

Рис. 2

1. ОПИСАНИЕ КОНСТРУКЦИИ МОТОГОНДОЛЫ

На самолете установлены мотогондолы с использованием в конструкции композиционных материалов (звукопоглощающие панели воздухозаборника).

Мотогондола (рис. 3) состоит из:

– передней части воздухозаборника;

– задней части (створки мотогондолы);

– панелей крепления створок мотогондолы.

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

Канал – трехслойная оболочка. Внутренняя обшивка (перфорированная) выполнена из алюминиевого сплава Д19чАТВ толщиной 1,8 мм, нагруженная обшивка – из сплава Д19чАТ = 1,2 мм.

Заполнитель: ТССП-Ф-10П, сотовый, с шестигранной ячейкой а = 10 мм.

Толщина панели – 20 мм.

Внешняя поверхность воздухозаборника – обечайка представляет собой клепанную оболочку с обшивкой из материала Д16-АТВ (травленая) с толщиной обшивки 1,8 мм, под двумя подкрепляющими до толщины равной 1,2 мм между ними.

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

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

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

В конструкцию воздухозаборника встроена противообледенительная система
(ПОС) с отбором горячего воздуха от третьей ступени компрессора высокого давления двигателя.

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

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

Носок воздухозаборника состоит из обшивки, поперечной диафрагмы, на которой крепится коллектор с частью трубы ПОС и шпангоута № 1. Шпангоут № 1 сборной конструкции имеет кольцевую форму и состоит из стенки, усиленной поясами и диафрагмами.

Коллектор входит в конструкцию противообледенительной системы воздухозаборника (ПОС). Звукопоглощающая канальная панель (ЗПК) конструктивно выполнена в виде двух дюралюминиевых обшивок, между которыми вклеен сотовый заполнитель. Со стороны проточной части обшивка перфорирована. ПО торцам панели приклеены профили для стыковки с носком по шпангоуту № 1 и со шпангоутом № 4 воздухозаборника.

2. СИЛОВОЙ РАСЧЕТ ВОЗДУХОЗАБОРНИКА

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

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

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

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

2.1. Исходные данные для силового расчета

Аэродинамические нагрузки на мотогондолу приведены в табл. 1,
(xy и xz даны в долях длины мотогондолы. В носке мотогондолы х = 0).

Таблица 1

Характеристика расчетных случаев А' и Д' для установок под двигатели
| | |
|Расчетные |Значения характеристик |
|случаи |nyэ |(, |(, |(звнутр|q, |yэмг , |xy |zэмг, |xz |
| | |град|град|, град |кг/м3 |кг | |кг | |
|А' |2,5 |10 |0 |0 |2000 |1600/ |0,16( |(190 |0,16(|
| | | | | | |1100 |0,83 | | |
| | | | | | | | | |0,55 |
|Д' |-1,0|-4 |0 |0 |2000 |-2210/ |0,16( |(160 |0,16(|
| | | | | | |-1810 |90,55 | | |
| | | | | | | | | |0,55 |

Нагрузки распределяются по внешней поверхности следующим образом:
– избыточное давление по поверхности определяется по формуле (1.1)

(Pэ = pq ,

(1.1)

где (Pэ – избыточное давление на поверхности; q – скоростной напор; p – рассчитывается по формуле:

p = p1+ py + pz .

(1.2)

Величина p1 определяется по графику на рис. 4

Величина py для случая Д' дается на прилагаемом графике (рис. 5). Для других режимов величина py пересчитывается пропорционально Yмг.

Значение pz определяется по формуле:

pz = pz( + pz( .

(1.3)

Распределение pz( по контуру и длине воздухозаборника дается на графике (рис. 6). При этом pz( определяется по выражению:

pz( = (z(()мг/q)Kz( .

(1.4)

В случаях А' и Д' z(()мг = zмг, в других расчетных случаях следует принимать z(()мг = (180 кг. Kz( определяется по графику на рис. 6.

Распределение pz( по контуру принимается таким же как и для pz(. При этом:

pz( = ((zмг – 180)/q)Kz( .

(1.5)

где zмг – берется из таблиц;

Kz( – определяется по графику на рис. 7.

2.2. Распределение расчетных аэродинамических нагрузок по длине воздухозаборника

Нагрузки на внутреннюю поверхность воздухозаборника представлены в таблицах 2 и 3.

Таблица 2

Расчетные значения нагрузок в случае А'
|х |0° |60° |120° |180° |240° |300° |
|0 |-1105 |-545 |564 |+1105 |545 |-564 |
| | |-804 |305 | |804 |-305 |
|0,05 |-940 |-464 |476 |+940 |464 |-476 |
| | |-679 |261 | |679 |-261 |
|0,1 |-774 |-383 |391 |+774 |383 |-391 |
| | |-553 |221 | |553 |-221 |
|0,153 |-597 |-296 |302 |+597 |296 |-302 |
| | |-431 |167 | |431 |-167 |

Таблица 3

Расчетные значения нагрузок в случае Д'
|х |0° |60° |120° |180° |240° |300° |
|0 |+442 |207 |-235 |-442 |-207 |235 |
| | |-12 |-454 | |12 |454 |
|0,05 |+376 |177 |-199 |-376 |-177 |199 |
| | |-3 |-379 | |3 |379 |
|0,1 |+310 |146 |-164 |-310 |-146 |164 |
| | |2 |-308 | |-2 |308 |
|0,153 |+239 |113 |-127 |-239 |-113 |127 |
| | |-1 |-241 | |1 |241 |

2.3. Распределение нагрузок по длине и по сечениям воздухозаборника

2.3.1. Несимметричное распределение нагрузки

Характер несимметричного распределения максимальных нагрузок по длине воздухозаборника в случае А' показан на рис. 8, а по сечению воздухозаборника на рис. 9

Распределение нагрузок по длине воздухозаборника

Рис. 8

Изменение максимальных нагрузок по сечению воздухозаборника

[pic]

Рис. 9
Расчетные нагрузки в случае А' и Д' определяются по формуле:

p = f·q?·(z/q)· Kz(

(1.6)

Нагрузки по длине мотогондолы определим, подставляя значения для случая А':

p = 2·2000?·(±190/2000)· Kz( = ±380Kz( .

В случае Д':

p = 2·2000?·(±160/2000)·Kz( = ±320Kz( .

Нагрузки по контуру мотогондолы определим, подставляя значения для случая А':

p = ((±190 – 180)/2000)·2·2000·Kz( = (20;-740)Kz( .

В случае Д':

p = ((±160 – 180)/2000)·2·2000·Kz( = (-40;-680)Kz( .

Суммарные нагрузки:

В случае А':

p = ±380 Kz( Kz(·(+20;–740) .

В случае Д':

p = ±320 Kz( Kz(·(-40;–680) .

2.3.2. Равномерное распределение нагрузки

Характер распределения нагрузки p1 по сечениям воздухозаборника приведен на рис. 10

Характер распределения нагрузки p1 по сечениям воздухозаборника

[pic]

Рис. 10

Таблица 4

| |Угол |Для всех углов |
| | |Расчетный случай |
| | |А' |Д' | | |
| | |скоростной напор – q, кг/м2 |
|х |Д' |2000 |2000 |680 |2000 |
|0 |1,66 |-6640 |-6140 | | |
|0,05 |1,02 |-4080 |-4080 | | |
|0,1 |0,86 |-3440 |-3440 | | |
|0,153 |0,76 |-3040 |-3040 | | |

2.3.3. Распределение py по воздухозаборнику

Характер распределения нагрузки py приведен на рис. 11.
Величина нагрузки py по воздухозаборнику:

py = (1600/2210)·2·2000 = 2895,93py* .

Распределение py по воздухозаборнику

[pic]

Рис. 11

Значения py* приведены в табл. 5.

Таблица 5

Значение нагрузки py*
|Сечение | |( |
|х |py* |0 | | | | | |
|0 |0,435 |-1259 |-630 |630 |1260 |630 |-630 |
|0,05 |0,370 |-1072 |-536 |536 |1072 |536 |-536 |
|0,1 |0,305 |-883 |-442 |442 |883 |442 |-442 |
| | |-883 |-883 |883 |883 | | |
|0,153 |0,235 |-681 |-681 |681 |681 |341 |-341 |
|0,1716 |0,210 |-608 |-608 |608 |608 |304 |-304 |

Коэффициент пересчета для случая Д':

Л = -1,3812 и py = -4000py*

2.3.4. Распределение нагрузки по воздухозаборнику от силы pz

Для случая А'

pz = ±380 Kz( ·(+20;-740) Kz(

Таблица 5

Распределение нагрузки по длине и по контуру от силы pz
| | | |( |
|х |Kz( |Kz( |0° |60° |120° |180° |240° |300° |
|0 |0,55 |0,395 |0 |-174 |-174 |0 |174 |174 |
| | | | |-72 |-72 | |72 |72 |
|0,05 |0,51 |0,325 |0 |-162 |-162 |0 |162 |162 |
| | | | |-40 |-40 | |40 |40 |
|0,1 |-0,42 |0,260 |0 |-134 |-134 |0 |134 |134 |
| | | | |-28 |-28 | |28 |28 |
|0,153 |-0,27 |0,205 |0 |-85 |-85 |0 |85 |85 |
| | | | |-42 |-42 | |42 |42 |

Суммарные аэродинамические нагрузки на воздухозаборник приведены в табл. 6, 7, 8 и 9

Таблица 6
Суммарные аэродинамические нагрузки на воздухозаборник в случае А' и

L = 3,8 м (Рр, кг/м2)
| |(, град |
|х |0° |60° |120° |180° |240° |300° |
|0 |-7900 |-7444 |-6184 |-5380 |-5836 |-7096 |
| | |-7342 |-6082 | |-5938 |-7198 |
|0,05 |-5752 |-4778 |-3706 |-3008 |-3382 |-4454 |
| | |-4656 |-3584 | |-3504 |-4576 |
|0,1 |-4323 |-4016 |-3132 |-2557 | | |
| | |-4457 |-2691 | |-2864 |-3748 |
| | |-3910 |-3026 | | | |
| | |-4351 |-2585 | |-2970 |-3854 |
|0,153 |-3721 |-3806 |-2444 |-2353 |-2614 |-3296 |
| | |-3763 |-2401 | |-2657 |-3339 |
|0,1716 |-3528 |-3581 |-2315 |-2312 |-2563 |-3171 |
| | |-3591 |-2375 | |-2553 |-3161 |

Таблица 7
Суммарные аэродинамические нагрузки на воздухозаборник в случае Д' (р =
±3200, и Kz(·(-40;-680) Kz(

| | | |( |
|х |Kz( |Kz( |0° |60° |120° |180° |240° |300° |

Продолжение табл. 7

|0 |-0,55 |0,395 |-6640* |-166 |-166 |0 |166 |166 |
| | | |0 |-80,2 |-80,2 | |80,2 |80,2 |
|0,05 |-0,51 |0,325 |-4080* |-152,5 |-152,5 |0 |153 |153 |
| | | |0 |-50 |-50 | |50 |50 |
|0,1 |-0,42 |0,260 |-3440* |-1254 |-1254 |0 |125,4 |125,4 |
| | | |0 |-36,7 |-36,7 | |36,7 |36,7 |
|0,153 |-0,27 |0,205 |-2920* |-82 |-82 |0 |82 |82 |
| | | |0 |-46 |-46 | |46 |46 |
|0,1716 |-0,17 |0,185 |-2560* |-54 |-54 | |54 |54 |
| | | |0 |-62 |-62 | |62 |62 |

*) Указаны значения равномерного распределения р1 по сечениям и по длине воздухозаборника

Таблица 8
Суммарные аэродинамические нагрузки на воздухозаборник в случае Д'

(К = -1,3812, py = -4000·py* (кг/м2)
| |( |
|х |0° |60° |120° |180° |240° |300° |
|0 |1740 |870 |-870 |-1740 |-870 |870 |
|0,05 |1486 |740 |-740 |-1486 |-740 |740 |
|0,1 |1220 |610 |-610 |-1220 |-610 |610 |
| | |1220 |-1220 | | | |
|0,153 |941 |941 |-941 |-941 |-471 |471 |
|0,1716 |840 |840 |-840 |-840 |-420 |420 |

Таблица 9
Суммарные расчетные аэродинамические нагрузки на воздухозаборник в случае
Д'
| |( |
|х |0° |60° |120° |180° |240° |300° |
|0 |-4900 |-5936 |-7676 |-8380 |-7344 |-5604 |
| | |-5850 |-7590 | |-7430 |-5690 |
|0,05 |-2600 |-3493 |-4973 |-5560 |-4667 |-3187 |
| | |-3390 |-4870 | |-4770 |-3290 |
|0,1 |-2220 |-2955 |-4175 |-4660 | | |
| | |-2345 |-4785 | |-3925 |-2705 |
| | |-2867 |-4087 | | | |
| | |-2257 |-4697 | |-4013 |-2793 |
|0,153 |-2100 |-2181 |-4063 |-3980 |-3429 |-2487 |
| | |-2145 |-4027 | |-3465 |-2523 |
|0,1716 |-2080 |-2134 |-3814 |-3760 |-3286 |-2446 |
| | |-2142 |-3822 | |-3278 |-2438 |

2.4. Распределение аэродинамических нагрузок на внутренней поверхности воздухозаборника

Нагрузки в канале от py в случае А':

q = 2000 кг/м2, Dвх = 1,6 м, f = 2,0, ( = -10°;

Sвх = (r2 = 2,01 м2, (' = 0,1745;

Y = Sвх·q·( = 2,01·2·2000·0,1745 = 1403 кг .

Нагрузки в канале от py в случае Д':

q = 2000 кг/м2, Dвх = 1,6 м, f = 2,0, ( = -4°;

Sвх = (r2 = 2,01 м2, (' = 0,0698;

Y = Sвх·q·( = -2,01·2·2000·0,0698 = -561 кг .

В случае А':

pz = (20;-740)Кz( ; py = (1403/2210)·2·2000·py* = 2539,3py* (кг/м2)

В случае Д':

pz = (-40;-680)Кz( ; py = (-561/2210)·2·2000·py* = -1015py* (кг/м2)

Таблица 10

Значения нагрузок в случае А' и Д' при ( = 0°
| | |Расчетный случай |
| | |А' |Д' |
|х |Кz( |pz = (20;-740), |pz = (-40;-680), |
| | |кг/м2 |кг/м2 |
|0 |0,395 |8 |-16 |
| | |-292 |-269 |
|0,05 |0,325 |7 |-13 |
| | |-241 |-221 |
|0,1 |0,260 |5 |-10 |
| | |-192 |-177 |
|0,153 |0,260 |4 |-8 |
| | |-152 |-140 |

Таблица 11

Значения нагрузок в случае А’ и Д’ при ( = 90°
| | |Расчетный случай |
| | |А’ |Д’ |
|х |py* |py = 2539,2 |py = -1015, |
| | |кг/м2 |кг/м2 |
|0 |-0,435 |1105 |-442 |
|0,05 |-0,370 |940 |-376 |
|0,1 |-0,307 |774 |-310 |
|0,153 |-0,235 |594 |-239 |

Таблица 12

Нагрузки на внутреннюю поверхность воздухозаборника в случае А’
| |(p = pycos( + pzsin( |
| |( |
|x |0° |60° |120° |180° |240° |300° |
|0 |-1105 |-552 |557 |1105 |557 |-552 |
| | |7 |7 |0 |-7 |-7 |
| | |-252 |-252 | |252 |252 |
|( |-1105 |-545 |564 |1105 |545 |-564 |
| | |-804 |305 | |804 |-305 |
|0,05 |-940 |-470 |470 |940 |470 |-470 |
| | |6 |6 |0 |-6 |-6 |
| | |-209 |-209 | |209 |209 |
|( |-940 |-464 |476 |940 |464 |-476 |
| | |-679 |261 | |678 |-261 |
|0,1 |-774 |-387 |387 |774 |387 |-387 |
| | |6 |4 |0 |-4 |-4 |
| | |-166 |-166 | |166 |166 |
|( |-774 |-383 |391 |774 |383 |-391 |
| | |-553 |221 | |553 |-221 |
|0,153 |-597 |-299 |299 |597 |299 |-299 |
| | |3 |3 |0 |-3 |-3 |
| | |-132 |-132 | |132 |122 |
|( |-597 |-296 |302 |597 |296 |-302 |
| | |-431 |167 | |431 |-167 |

Таблица 13

Нагрузки на внутреннюю поверхность воздухозаборника в случае Д’
| | |(p = pycos( + pzsin( |
| | |( |
| | |0° |60° |120° |180° |240° |300° |
|х |pycos( |442 |221 |-221 |-442 |-221 |221 |
|0 |pzsin( | |-14 |-14 | |14 |14 |
| | | |-233 |-233 | |233 |233 |
| |( |442 |207 |-235 |-442 |-207 |235 |
| | | |-12 |-454 | |-12 |454 |
| |pycos( |376 |188 |-188 |-376 |-188 |188 |

Продолжение табл. 13
|0,05 |pzsin( | |-11 |-11 |0 |11 |11 |
| | | |-191 |-191 | |191 |191 |
| |( |376 |177 |-199 |-376 |-177 |199 |
| | | |-3 |-379 | |3 |379 |
| |pycos( |310 |155 |-155 |-310 |-155 |155 |
|0,1 |pzsin( | |-9 |-9 | |9 |9 |
| | | |-153 |-153 | |153 |153 |
| |( |310 |146 |-164 |-310 |-146 |164 |
| | | |2 |-308 | |-2 |308 |
| |pycos( |239 |120 |-120 |-239 |-120 |120 |
|0,153 |pzsin( | |-7 |-7 | |7 |7 |
| | | |-121 |-121 | |121 |121 |
| |( |239 |113 |-127 |-239 |-113 |127 |
| | | |-1 |-241 | |1 |241 |

2.5. Определение равнодействующих нагрузок по сечениям воздухозаборника от внешних и внутренних аэродинамических нагрузок

Суммарное распределение нагрузки в поперечном сечении воздухозаборника

[pic]

Рис. 16

Расчет нагрузок от внешних аэродинамических сил (для нижних значений pz производится по формулам:

[pic] , (1.7)

[pic] . (1.8)

Принимаем значение (r = 2,826 м;

Рассчитанные значения нагрузок по формулам (1.7) и (1.8) представлены в табл. 14, 15

Таблица 14

Суммарные значения нагрузок в случае А'
|х |pycos( |pzsin( |qy |qz |q(, кг/м |(, град |
|0 |-1260cos( |-83sin( |-3561 |-235 |-3569 |3,8 |
|0,05 |-1072cos( |-47sin( |-3029 |-133 |-3032 |2,5 |
|0,1 |-883cos( |-33sin( |-2495 |-93 |-2497 |2,13 |
| | | |-2838 | |-2840 |1,88 |
|0,153 |-681cos( |-49sin( |-1925 |-138 |-2143 |3,7 |
| | | |-2138 | | | |

x = 0,1; -1589,4 – 2495·0,5 = -2838 кг/м; х = 0,153; -1226 – 1925·0,5 = -2139 кг/м.

Таблица 15

Суммарные значения нагрузок в случае Д'
|х |pycos( |pzsin( |qy |qz |q(, кг/м |(, град |
|0 |1740cos( |-93sin( |4917 |-263 |4924 |3,06 |
|0,05 |1481cos( |-58sin( |4185 |-164 |4188 |2,25 |
|0,1 |1220cos( |-42sin( |3448 |-119 |3450 |1,98 |
| | | |3893 | |3895 |1,75 |
|0,153 |941cos( |-53sin( |2659 |-150 |3028 |2,84 |
| | | |3024 | | | |

x = 0,1; 0,5·3448 + 2169 = 3893 кг/м; х = 0,153; 0,5·2659 + 1694 = 3024 кг/м.

2.6. Нагрузки на болты крепления воздухозаборника к проставке

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

2.6.1. Определение нагрузок на болты крепления в случае А'

Для определения нагрузок на болты крепления воздухозаборника к проставке примем:

– число болтов n = 12;

– Dокр. болтов = 1440 мм;

Распределение суммарной погонной нагрузки в точках Д, С, В, А определяем как:

qД = 3835 + 3777 = 7607 кг/м; qС = 3157 + 3046 = 6203 кг/м; qВ = 2526 + 2425 = 4951 кг/м; qА = 2000 + 1977 = 3977 кг/м.

Распределение суммарной погонной нагрузки по длине представлено на рис. 17

Распределение суммарной погонной нагрузки по длине воздухозаборника

[pic]

Рис. 17

Величина суммарной приведенной нагрузки R в центре давления определяется как:

R = ((7607 + 6703)/2 + (6203 + 4956)/2)·0,19 + ((4956 + 3977)/2)·0,202 =

= 3274 (кг).

Для определения координаты центра давления определим суммарный изгибающий момент МА:

МА = 6203·0,19·0,487 + 4956·0,19·0,297 + 3977·0,202·0,101 +

1404·0,19·0,5·0,518 + 1247·0,19·0.5·0,329 + 979·0,202·0,5·0,135 = 1056 кг·м.

Координата центра давления хц.д.= 1056/3274 = 0,3225 м .

Расчетные нагрузки на болты определяем по формулам [6]:

Рmax = 4M/nDокр.б. ,

(1.9)

Рmax = (4·0,3235·3274)/(12·1,44) = 245 кг .
Срезающая нагрузка буртика (зуба) проставки:

Рсрр = 3274 кг .

Вес воздухозаборника Gв-ка = 93 кг, хц.т. = 350 мм вперед от плоскости крепления к проставке.

Нагрузки на болты крепления воздухозаборника от инерционных нагрузок представлены на рис. 18.

[pic]

Рис. 18

Задаемся коэффициентом перегрузки n = 1,5, тогда

Ринрц = Gв-ка·n = 93·1,5 = 140 (кг) .

М = 0,35·140 = 49 (кг·м) .

Рб = (4·43)/(12·1,44) = 11,34 (кг) .

Суммарный Мизг = [pic]= 1059,271 (кг·м) .

Максимальная растягивающая нагрузка на болт Рболт = 245,2 кг.

Схема расположения крепежных болтов по контуру и суммарные действующие нагрузки приведены на рис. 19

Схема расположения крепежных болтов по контуру и суммарные действующие нагрузки

[pic]

Рис. 19

2.6.2. Определение нагрузок на болты крепления в случае Д'

Расчетные нагрузки на воздухозаборник по сечениям и по длине, точки приложения равнодействующих, рассчитанные значения моментов и перерезывающих сил приведены на рис. 20.
Максимальное растягивающее усилие на болт:

Рmaxр = 4М/4d = (4·0,392·2400)/(12·1,440) = 218 кг.
Срезающая нагрузка воспринимается буртиком проставки – Рсрр = 2400 кг.

2.7. Проверка прочности воздухозаборника самолета

2.7.1. Исходные данные для расчета

Внутренняя обшивка: ( = 1,8 мм, материал: сплав Д19, перфорация – диаметром 2 мм.

Расчетные нагрузки на воздухозаборник в случае Д'

Рис. 20

Заполнитель: ТССП-Ф-10П (ТУ-596-258-87), удельный весь заполнителя – (
= 35±5 кг/м3; (сж = 15 кг/см2. Параметры заполнителя и перфорированной обшивки приведены на рис. 21.

Параметры заполнителя и перфорированной обшивки

[pic]

Рис. 21

Внешняя обшивка: ( = 1,2 мм, материал: сплав Д19.

Обечайка изготовлена из Д16Т, ( = 1,8 мм, травленная с ( = 1,8 мм до (
= 1,2 мм. Максимальный размер клетки 101 на 120 мм. Характерные размеры и сечения представлены на рис. 22

Типовое сечение обечайки
[pic]

Рис. 22

2.7.2. Расчет сечения в районе проставки в расчетном случае А'

Размеры рассчитываемого сечения приведены на рис. 23

Рис. 23

Находим момент инерции сечения:

I = ((0,4D3() = 0,4·142,53·0,12 + 0,4·1923·0,12 + 0,4·138,53·(пр =

638037,84 см4 .

Приведенная толщина внутренней общивки:

(пр. внутр. обш. = [(((·138,5)/12,0208)·0,2·0,18 – (·138,5·0,18]/((·138,5)

.

Нормальные напряжения от изгиба воздухозаборника:

( = (М·d)/J·2 = (3948·192·38,2)/(638037,84·2) = 22,69 (кг/см2),

Избыток прочности

( = 2750/22,69 – 1 >> 1.

М = Р·l ; Р = 3948 кг; l = 38,2 см. q = 22,69·0,12 = 2,72 кг/см

Проверяем ячейку травления на устойчивость от q = 2,72 кг/см. Схема нагружения ячейки приведена на рис. 24. Принимаем, что длинные края ячейки обшивки оперты

Схема нагружения ячейки обшивки

[pic]

Рис. 24

Величина a/b = 101/120 = 0,841; К = 3,6.

(кр = 2750 кг/см2,

( = 2750/355 – 1 = 6,746 ,

( >> 1

2.7.3. Проверка прочности внутреннего канала на осевое сжатие

Проверку прочности внутреннего канала на осевое сжатие проведем по методике изложенной в [6]:

Тдейств. = [P·l·(d + d1)((в + (н)((d + d1)]/2J = [3948·38,22(138,5 +

142,5)2(0,12 + + 0,15)·3,14]/(638037,94·4) = 3958 (кг)

Действующая сжимающая нагрузка от qp равна 2000·1,5 = 3000 (кг/м2).

Т = ((/4)(1922 – 1382)·0,3 = 4198,74 (кг).

Суммарная нагрузка: (Т = 8157 кг.
Заполнитель маложесткий. Расчетные формулы для трехслойных панелей (6):

(зап < 1,21qEпр ,

Li = E1H/E1B = 1,

( = (C + 1/[2,6(1 + 50)],

[pic]

b = 1,21qEпр/Gзап

С1 = D1рас/D1

D1 = 4(z0 – h – (н)3 + 4(H – z0)3 + 4li[z03 – (z0 – (н)3] ,

z0 = [(в2 + 2(в((н + h) + li(н2]/[2((в + li(н)] .

Расчет по приведенным выше формулам дает:

z0 = [0,152 + 2·0,15(0,12 + 2,8) + 1·0,122]/[2(0,15 + 1·0,12] = 1,246 ,

D1 = 4(1,246 – 2,8 – 0,12)3 + 4(2,27 – 1,246)3 + 4·1(1,2463 – (1,246 –

1,123) = = 3,652 ,

B2 = 0,15 + 0,12 = 0,27 (мм) .

Епр = 6,8·105 кг/см2. Принимаем для маложесткого заполнителя К = 0,2.

D1рас = (в3 + li(н3 = 0,153 + 0,123 = 0,00513 .

C1 = 0,005103/3,652 = 0,001397 . q = [2,0(1 – 0,001397)·(0,27·3,652]/[70,25(2,0 + 2,27)2] = 0,001548 .

Приведенный модуль сдвига:

Gзап = Gxz = 1,5·((c/t)·Gм ,

Gзап = Gyz = ((c/t)·Gм ,

Модуль сдвига заполнителя: Gм = Ем/[2(1 + ()],

Gм = 6000/[2(1 + 0,25)] = 2400 (кг/см2) .

Gзап = Gxz = 1,5·(0,025/1,732)·2400 = 52 (кг/см2) ,

Gзап = Gyz = (0,025/1,732)·2400 = 35 (м/см2),

Gзап = (52·35 = 42,7 (кг/см2)

42,7 < 1,21·0,001548·6,8·105;

42,7 < 1273,7 т.е. заполнитель маложесткий.

b = 1273,7/42,7 = 29,83

( = (0,001397 + [1/2·29,83·(1 + 5 – 0,001397)] = 0,054025

Критическая осевая сила Ткр:

Ткр = 2·(·К·Епр(ВzD1 ·( = 2·(·0,2·6,8·105(0,27·3,652 ·0,054025 = 45842 кг.

Избыток прочности

( = 45842/8157 – 1 = 4,62 .

Расчет напряжений во внутренних и наружных слоях трехслойной панели проведем в соответствии с [6]. Схема нагружения представлена на рис. 25.

Схема нагружения трехслойной панели

[pic]

Рис. 25

qв = q(1/(1 + ()); qн = q(q/(1 + ());

( = l·((н/(в) , l = E1н/Е1в = 6,8·105/6,8·105 = 1 , q = 8157/((·140,5) = 18,48 (кг/см2),

( = 1·(1,2/1,5) = 0,8 , qв = 18,48(1/(1 + 0,8)) = 10,27 (кг/см)

[(02] = 27,5 (кг/мм2). qн = 18,48(0,8/1,8) = 8,21 (кг/см),

(вн = 1027/0,15 = 68,5 (кг/см2),

(н = 8,21/0,12 = 68,42 (кг/см2) .

Избыток прочности:

( = 27,5/0,685 - 1 = >> 39,14 .

2.7.4. Проверка прочности внутреннего канала на внешнее давление

Расчетные нагрузки:

1. Установившийся режим Н = 0; М = 0;
Разрежение на входе в заборник распространяется на всю длину канала:

(p( = -0,645 кг/см2 ; (ст.соты = 0,04 ;

(м = 2400 кг/см2 ; (xz = 83 кг/см2 ;

Gyz = 55,42 кг/см2 ;

Gзап = (35,4·83 = 67,8 кг/см2 .

Определяем Ркр для несимметричной трехслойной оболочки с мягким средним слоем (рис. 26)

[pic]

Рис. 26

[pic]

li = 1 = Eн/Ев ; К = 0,8 .
Заполнитель маложесткий: g = h·B1/l·R0,5 ,

[pic]

a = 5g·Eпр/Gзап , с1 = D2рас/D2 ,

D2рас = (в3 + li(н3 = 0,153 + 0,123 = 0,00513 (см3) .

D2 = 4(z0 – h – (н)3 + 4(H – z1)3 + 4li[z03 – (z0 – (н)3] , z0 = [(в2 + 2(в((н + h) + li(н2]/[2((в + li(н)] . z0 = [0,152 + 2·0,15(0,12 + 2,8) + 0,122]/[2(0,15 + 0,12] = 1,2461 .

D2 = 4(1,246 – 2,8 – 0,12)3 + 4(2,27 – 1,246)3 + 4·1(1,2463 – (1,246 –

1,123) = = 3,6515 , с1 = 0,005103/3,6515 = 0,0013975 .

[pic] = 6,627·10-4 .

Gзап = 67,8 кг/см2 . а = 5·6,627·10-4·6,8·105/67,8 = 33,22.

Заполнитель маложесткий:

Gзап 40 |200 |

При конструировании и изготовлении защитного устройства должны быть учтены особенности конкретного пресса и условия работы на нем.

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

Винтовые прессы с балансирами должны иметь ограждения пути, проходимого этими балансирами.

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

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

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

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

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

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

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

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

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

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

4.2. Обеспечение чистоты производственного процесса.

Работа на прессе не сопровождается: а) загрязнением воздуха значительными вредными выделениями

(ядовитыми газами ,парами); б) значительными вибрациями ; в) воздействием на рабочего электромагнитными полями;

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

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

5. Расчет технико - экономической эффективности изготовления воздухозаборника из композиционных материалов.

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

В процессе эксплуатации конструкций из КМ были выявлены основные преимущества:

- Малая масса по сравнению с традиционными типами подкрепленных пластин и оболочек.

- Экономичность по сравнению с традиционными конструкциями.

- Хорошие теплоизолирующие свойства.

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

Себестоимость канала воздухозаборника определяем по формуле :

С = М + ПФ + Зо + Зд + Зсс + НРу, (1.10)

где М -стоимость материалов;

ПФ -стоимость полуфабрикат;

Зо -зарплата основная;

Зд -зарплата дополнительная;

Зсс-отчисление на соцстрахование;

НРу-накладн. цеховых расходов.

Веса деталей канала воздухозаборника приведнны в табл. 32

ЛИТЕРАТУРА

1. C.М. Егер ‘Проектирование самолетов’ 1983г.

2. A.Н. Глаголев ‘Конструкция самолетов’ 1975г.

3. C.И. Зоншайн ‘Аэродинамика и конструкция летательных аппаратов’ 1966г.

4. И.А. Максимов, В.А. Секистов ‘Двигатели самолетов и вертолетов’

1977г.

5. Сборник трудов ‘Теория и практика проектирования пассажирских самолетов‘ 1976г.

6. В.Т. Лизин ,В.А. Пяткин ‘Проектирование тонкостенных конструкций’.

Приложение

Листинг программы расчета аэродинамических нагрузок на мотогондолу

real*8 fi,Po(14)/1.65,1.69,1.67,0.98,0.88,0.82,0.78,0.56,

* .35,.23,.17,.15,.14,.15/ real*8 Py(14)/3.4,3.4,3.4,3.3,2.5,1.1,

* .7,.55,.55,.55,.56,.45,-2.,-1./,

* Pza(14)/.26,.26,.26,.26,.26,.26,.26,.27,.43,.67,

* .91,.9,.6,0./,

* Pzb(14)/3.95,3.95,3.9,3.75,2.7,1.6,1.2,.62,.4,.3,

* .2,.15,.12,0./,

* X(14)/0.,.0125,.025,.05,.1103,.15,.181,.3,.4304,

* .55,.65,.774,.9,1./ real*8 gamm(14)/44.,24.,16.,7.,5.,2.,1.5,-2.,-8.,-10.,-
11.,

* -12.,-12.,-12./ real*8
R(14)/1.,1.06,1.11,1.145,1.2,1.21,1.225,1.24,1.16,

* 1.15,1.07,.95,.815,.72/ real*8 alf,bett,q,t,Ln,Lk,ref,fin,sh,dfi,r1ef,r2ef,fif,

*
Myx(5),Mzx(5),L,Pyc(5),Pxc(5),Pzc(5),P(5),Xc(5),Xcc(5),xt,pi write(*,*) 'Введите начальное значение угла Fi' read(6,*) fin write(*,*) 'Введите конечное значение угла Fi' read(6,*) fi write(6,*) 'Введите начальную Xотн' read(6,*) Ln write(6,*) 'Введите конечную Xотн' read(6,*) Lk write(*,*) 'Введите значение угла BETTA' read(*,*) bett write(*,*) 'Введите значение угла ALFA' read(*,*) alf write(*,*) 'Введите значение q' read(*,*) q pi=4*datan(1) fi=fi*pi/180. fin=fin*pi/180. alf=alf*pi/180. bett=bett*pi/180.

L=5.6 sh=.003 dfi=2.*pi/1257. xt=Ln do 4 i=1,13 if(xt.lt.x(i)) goto 5

4 continue

5 i=i-1 do 103 ik=1,5

Pyc(ik)=0.

Pzc(ik)=0.

Pxc(ik)=0.

Myx(ik)=0.

Mzx(ik)=0.

Xc(ik)=0.

103 Xcc(ik)=0. write(*,*) 'Номер участка',i gamm(i)=gamm(i)*pi/180. do 1 t=Ln,Lk,sh if(xt.gt.x(i+1)) then i=i+1 write(*,*) 'Номер участка',i gamm(i)=gamm(i)*pi/180. endif do 2 fif=fin,fi,dfi c write(*,*) fif*180./pi,xt,r(i) c Учет пилона **********************************

C if((xt.gt.x(5).and.xt.lt.x(6)).and.

C * (fif.lt.1.449.or.fif.gt.1.693)) goto 6

C if((xt.gt.x(6).and.xt.lt.x(7)).and.

C * (fif.lt.1.344.or.fif.gt.1.798)) goto 6

C if(xt.gt.0.1103.and.(fif.gt.1.2915.and.fif.lt.1.85))

C * goto 3

6 continue r1ef=dcos(fif+dfi/2.) r2ef=dcos(fif+dfi/2.) if(xt.gt.0.1103.and.(fif.gt.0..and.fif.le.pi/2.))then r1ef=1. r2ef=1. c write(*,*) '1 *********',fif c pause endif if(xt.gt.0.1103.and.(fif.gt.pi/2.and.fif.le.pi))then r1ef=-1. r2ef=-1. c write(*,*) '2 *********',fif c pause endif c write(*,*) fif*180./pi,xt,r1ef,r2ef c PAUSE' '

P(1)=(Po(i)+Py(i)*alf*dsin(fif+dfi/2.)+Pza(i)*alf*r1ef

* -Pzb(i)*bett*r2ef)*q*L*dfi*R(i)*sh c write(*,*) P/(sh*L*R(i)*dfi),xt,fif c pause

P(2)=Po(i)*q*L*dfi*R(i)*sh

P(3)=Py(i)*alf*dsin(fif+dfi/2.)*q*L*dfi*R(i)*sh

P(4)=Pza(i)*alf*r1ef*q*L*dfi*R(i)*sh

P(5)=-Pzb(i)*bett*r2ef*q*L*dfi*R(i)*sh do 1000 ik=1,5

Pyc(ik)=Pyc(ik)+P(ik)*dsin(fif+dfi/2.)

Pzc(ik)=Pzc(ik)+P(ik)*dcos(fif+dfi/2.)

Myx(ik)=Myx(ik)+P(ik)*dsin(fif+dfi/2.)*xt

Mzx(ik)=Mzx(ik)+P(ik)*dcos(fif+dfi/2.)*xt

Pxc(ik)=Pxc(ik)+P(ik)*dtan(gamm(i))

1000 continue

3 continue

2 continue xt=xt+sh

R(i)=R(i)+sh*L*dtan(gamm(i))

1 continue do 100 ik=1,5 if(dabs(Pyc(ik)).gt.1.d-10) then

Xc(ik)=Myx(ik)/Pyc(ik) else write(*,*)ik,'Myx=',Myx(ik) endif if(dabs(Pzc(ik)).gt.1.d-10) then

Xcc(ik)=Mzx(ik)/Pzc(ik) else write(*,*)ik,'Mzx=',Mzx(ik) endif write(*,*) 'Значения аэр. сил в Ц.Д.'

100 write(*,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L

12 format(1x,'Pxc=',f10.2,/,1x,'Pyc=',f10.2/,1x,'Pzc=',f10.2/,

* 1x,'Xyc=',f15.7,/,1x,'Xzc=',f15.7) open(1,file='aer.res') write(1,*) 'Начальное значение угла Fi' write(1,*) fin*180./pi write(1,*) 'Конечное значение угла Fi' write(1,*) fi*180./pi write(1,*) 'Начальная Xнач' write(1,*) Ln*5.6 write(1,*) 'Конечная Xкон' write(1,*) Lk*5.6 write(1,*) 'Значение угла BETTA' write(1,*) bett*180./pi write(1,*) 'Значение угла ALFA' write(1,*) alf*180./pi write(1,*) 'Значение q' write(1,*) q write(1,*) 'Значения аэр. сил в Ц.Д.' do 102 ik=1,5

102 write(1,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L close(1) stop' ' end