Большинство серийно выпускаемых в мире вертолетов оснащены шарнирным несущим винтом (НВ) с малым разносом горизонтальных шарниров. Такая конструктивная схема вполне соответствовала требованиям своего времени. При этом общим конструктивным недостатком всех типов винтов, использующих различные шарниры, является необходимость устранения проблем износных явлений при эксплуатации вертолета. Кроме этого, любой подшипник предполагает наличие определенных люфтов, что в свою очередь отрицательно сказывается на уровне вибраций вертолета, а также неизбежно ведет к увеличению массы самой втулки из-за необходимости плавной передачи нагрузок с подшипника на упругие части втулки.
Выполнение задач сегодняшнего дня требует от вертолета улучшения его потребительских качеств: уменьшения уровня шума, повышения маневренности и удешевления стоимости его эксплуатации. Это возможно при условии совершенствования несущей системы вертолета совместно с широким применением композиционных материалов, увеличивающих ресурс. Поэтому в последние десятилетия на новых перспективных вертолетах стали применяться так называемые «бесшарнирные» несущие винты, роль шарниров в которых выполняют специальные упругие элементы. С одной стороны применение элементов такого типа обеспечивает существенное упрощение конструкции несущего винта, и в значительной степени повышает эксплуатационные характеристики и маневренность вертолета.
Как показывает опыт эксплуатации, такие винты имеют больший ресурс по сравнению с обычными несущими винтами с горизонтальными и вертикальными шарнирами. Тем самым удешевляется эксплуатация вертолета, но при этом увеличиваются затраты на проектирование и изготовление таких конструкций. Поэтому точность прогнозирования нагружения и, соответственно, оценки ресурса несущей системы вертолета является на сегодняшний день одной из ключевых задач вертолетостроения.
На ОАО «Казанский вертолетный завод» в серийное производство вышел легкий многоцелевой вертолет «Ансат» (рис. 1).
Рис. 1. Вертолет «Ансат».
Рис. 2. Бесшарнирная втулка вертолета Ансат.
Несущий винт этого вертолета оснащен бесшарнирной втулкой, в которой функции горизонтального, вертикального и осевого шарниров выполняет упругий элемент протяженного типа — торсион.
Торсион обеспечивает восприятие центробежных растягивающих усилий, и передачу на втулку сил и моментов, приходящих с лопасти. Передача управляющих усилий происходит через кожух торсиона жестко прикрепленного к комлевой части лопасти и шарнирно опертого в комлевой части торсиона. Торсион состоит из трех основных участка: комлевого, рабочего и концевого (рис. 3). Комлевый участок, предназначен для передачи сил и изгибающих моментов на вал главного редуктора. Рабочий упруго-деформируемый участок разделен продольными прорезями на ручьи для снижения крутильной жесткости балки. Изгибом и закруткой этого участка торсиона обеспечиваются маховое движение и качание лопасти, а также изменение угла установки лопасти. Концевой участок предназначен для соединения с переходником, обеспечивающего крепление лопасти.
Новый конструктивный элемент несущей системы привлек внимание многих исследователей. Вопросам расчета напряженно-деформированного состояния универсального торсиона НВ вертолета «Ансат», представляющего собой стержневую многослойную композитную конструкцию, посвящено достаточно много работ [1, 2, 3, 4]. В которых в рамках сдвиговой модели С. П. Тимошенко выводится система нелинейных дифференциальных уравнений упругого деформирования торсиона и граничных условий. В предположении отсутствия распределенной нагрузки по длине получено аналитическое решение для отдельного стержня и определены границы этого решения. Отмечено, что касательные напряжения в модели типа С. П. Тимошенко в общем случае не удовлетворяют статическим граничным условиям на боковой поверхности стержня и, следовательно, требуются уточнения.
В исследованиях А. Ю. Лисса [5] уточнен подход, предложенный в работе [3] и представлена методика расчета торсиона балочного типа на изгиб в двух плоскостях и кручение.
Исследования напряженно-деформированного и предельного состояния многослойных композиционных торсионов сложной формы с использованием конечно-элементных методов представлены в работах [6, 7, 8,9,10,11,12].
Эти работы позволили получить хорошее представление о прочностных и деформационных характеристиках изолированного торсиона.
С другой стороны, торсион является конструктивным элементом втулки и не эксплуатируется вне несущего винта. При этом нагрузки, действующие на торсион, определяются силами и моментами, приходящими со стороны лопасти, которые в свою очередь во многом зависят от упруго-махового движения, определяемого условиями закрепления, т. е. деформационными свойствами упругого элемента, а также режимом полета вертолета. Кроме этого, отсутствие горизонтального и вертикального шарниров приводит к передаче изгибающих моментов, действующих на лопасти и значительной мере определяемых ее упругостью, на вал НВ и, следовательно, оказывают влияние на баланс нагрузок вертолета в целом. Следовательно, для того, чтобы определить нагружение торсиона в полете необходимо решить задачи балансировки и динамики полета вертолета в комплексе с аэроупругим несущим винтом.
Математическое моделирование нагружения вертолета можно выполнить двумя способами. Первый — вычислением в процессе решения уравнения движения. Второй — определением сил и моментов заранее и вводом их в вычислительные машины в виде таблиц, графиков или неких апроксимационных зависимостей. Каждый из способов имеет свои преимущества и недостатки, и применяется в зависимости от решаемой задачи. Наиболее трудоемким при решении уравнений равновесия вертолета является вычисление сил и моментов, действующих на вертолет со стороны несущего и рулевого винтов. Поэтому в большинстве существующих моделей расчета балансировки и динамики полета вертолета принято использовать второй способ, хотя при этом приходится прибегать к упрощающим допущениям, так как возникают сложности моделирования зависимостей нагружения от всех параметров, влияющих на них. Кроме этого изменения параметров втулки и лопастей требует повторного пересчета характеристик винта, но возможность проведения расчетов в режиме реального времени, конечно, окупает и эти сложности. Этот способ действительно хорошо применим, но только в задачах моделирования процесса пилотирования летчиком. Следует отметить, что этот способ тоже исследовался в данной работе.
В зарубежной науке выбрано первое направление, но инструментарий, который при этом выбран, требует очень больших вычислительных ресурсов и затрат времени. При наличии высокопроизводительных кластерных вычислительных систем это вполне оправдано, особенно для решения сложных вопросов аэродинамики. В задачах же аэроупругости, особенно вращающихся систем особую роль играют проблемы механики, и прочности, решение большинства которых на сегодняшний день основано, на конечно-элементных программных пакетах. Это опять же требует больших затрат вычислительных ресурсов. Поэтому на этапах летных испытаний и сертификации нового образца вертолетной техники применение таких сложных и дорогостоящих систем не всегда приемлемо.
Поэтому для комплексного решения задач динамики полета, балансировки вертолета и прочности бесшарнирного несущего винта необходимы другие подходы. Пусть менее универсальные методы и алгоритмы, но более эффективные и направленные на решение конкретных задач.
Таким образом, в данной работе решается научная проблема — развития теории и методов расчета нагружения вертолета оснащенного бесшарнирным несущим винтом с упругим элементом протяженного типа в произвольном полете.
В первой главе представлена математическая модель аэроупругого бесшарнирного винта механика, которой основана на стержневой теории. Первые работы по аэродинамическому и прочностному расчету винтов появились еще в конце девятнадцатого в начале двадцатого веков. Это в основном работы Н. Е. Жуковского, С. К. Джевецкого, Г. Х. Сабинина, Б. Н. Юрьева [13], В. П. Ветчинкина. Особенно много работ в исследовании аэродинамики и аэроупругости несущих винтов появилось с началом активного применения вертолетов в 50-х годах нашего столетия [14, 15, 16, 17, 18,19, 20,21].
Основой аэродинамики несущего винта в горизонтальном полете являются работы Глауэрта и Локка, содержащие значительное количество упрощающих допущений, которые пришлось ввести ввиду математической сложности уравнений, описывающих поведение несущего винта на этом режиме полета. По мере развития теории и создания летающих винтокрылых аппаратов некоторые из допущений, содержащихся в работах Глауэрта, изучались и заменялись более точными положениями, которые совместно с теорией элемента лопасти позволяют получать надежные результаты! Более точное представление воздушного потока в зоне работы несущего винта вертолета дает вихревая теория [23], но широкое применение этой теории для текущих расчетов затруднено в виду недостаточности быстродействия даже самых современных вычислительных средств.
Достаточно хорошо освещены в литературе также вопросы расчета лопастей на статическую и динамическую прочность [22, 23, 24]. В традиционных методах при расчете лопастей несущих винтов в качестве расчетной схемы принимается тонкий, естественно закрученный стержень с прямолинейной осью жесткости. Упругие перемещения такого стержня под действием нагрузки полагаются малыми, что позволяет исключить нелинейные члены в записи уравнений равновесия. В настоящее время теория тонких стержней, известная как теория Кирхгофа-Клебша, имеет широкое применение при исследовании прочности и устойчивости стрежневых систем. *.
В общем случае, без предположения малости упругих перемещений, задача расчета деформаций изгиба лопасти в двух плоскостях и кручения является геометрически нелинейной, так как силы в срединной плоскости вызывают не только изгиб, но и кручение. Поэтому при решении задач аэроупругости конечность перемещений скажется на составляющих аэродинамической и инерционной нагрузки. Применение теории больших перемещении до сих пор имеет ограниченный характер, хотя довольно широко используется теория упругой линии двоякой кривизны [25, 26, 27]. Много работ посвящено также исследованию задач статики, динамики и устойчивости авиационных конструкций, базирующихся на стрежневой расчетной схеме. Как правило, в них учитывается конечность перемещений, но в разрешающих уравнениях накладываются ограничения на величину этих перемещений или линеаризуются сами уравнения.
С появлением и быстрым развитием средств вычислительной техники новый подъем происходит и в развитии методов расчета. В середине 60-х годов нашла широкое применение методика расчета деформаций лопасти несущего винта, разработанная А. В. Некрасовым [28, 29]. На основе этой методики была разработана программа для расчета изгиба лопастей в плоскости взмаха и кручения. Деформации лопасти разлагаются по формам собственных изолированных колебаний.
В начале 70-х годов наиболее существенный вклад в развитие методов расчета деформаций лопастей несущих винтов внесли работы А. Ю. Лисса [30, 31, 32]. А. Ю. Лиссом в разложении деформаций применены формы связанных собственных колебаний лопасти с учетом изгиба в двух плоскостях и кручения.
Можно также отметить работы Т. Д. Смоляниновой [33, 34, 35] в которых предложены уравнения изгибно-изгибно-крутильных колебаний лопасти с произвольной формой упругой оси до нагружения. Уравнения колебаний и выражения для погонных сил и моментов получены в матричной форме и в проекциях на оси различных систем координат.
Решение задачи расчета деформаций без предположения малости упругих перемещений стало возможным благодаря развитию эффективных численных методов решения задач строительной механики. Они позволяют заменить дифференциальные уравнения системой нелинейных алгебраических уравнений. С появлением этих методик теория больших перемещений тонких стержней получила дальнейшее развитие в работах В. А. Павлова и его учеников.
В этих и других ранее опубликованных работах при выводе уравнений колебания лопасти несущий винт рассматривался изолированно. В такой постановке силы и моменты, возникающие на упругих лопастях за счет движения вертолета в пространстве, не учитывались. Этот подход оправдан при рассмотрении установившихся режимов полета, и конечно, при отсутствии вращения вертолета в пространстве. Однако при моделировании неустановившихся режимов полета пренебрежение этими составляющими может привести значимым погрешностям.
В статье [36] представлены наиболее полные и математически строгие на сегодняшний день уравнения колебания упругой лопасти в матрично-векторной форме при произвольном пространственном движении вертолета. В этой работе вывод уравнений, определяющих массово-инерционную нагрузку, выполнен в следующем порядке. Вначале радиус-вектор сечения лопасти путем ряда последовательных обратных переходов проецируется в земную систему координат. Затем в ней выполняется дифференцирование. Далее путем прямых последовательных переходов полученные ускорения проецируются на оси связанной с сечением лопасти системы координат. Результат представлен в компактной матрично-векторной форме. Но следует отметить, что практически все составляющие матриц и векторов зависят от времени, и соответственно при двукратном дифференцировании суммы произведений обратного и прямого перехода количество компонент значительно увеличиться. Очевидно, что при преобразовании этих уравнений в алгебраический вид они будут достаточно громоздкими, и это приведет к значительному росту затрат машинного времени при вычислениях.
В данной работе предлагается применить другой способ, основанный на разделении движения на переносное и относительное. При этом дифференцирование будет производиться во вращающейся системе координат, что позволит минимизировать число переходов, и тем самым получить более рациональные уравнения упруго-махового движения лопасти винта вертолета, построенные на основе нелинейной теории больших перемещений тонких стержней.
Применение принципа Даламбера позволяет представить инерционную нагрузку в виде внешних сил и моментов и, таким образом, проинтегрировать уравнения равновесия по времени с использованием одного из численных подходов.
Аэродинамическая нагрузка в этой работе на лопастях несущего винта определяется по элементно-импульсной теории. При этом применение теории несущей линии не вполне оправдано вблизи концов крыла. Если в концевом сечении хорда лопасти конечна, то теория элемента лопасти дает ненулевую подъемную силу. Однако в действительности нагрузка на конце лопасти уменьшается до нуля, причем спад происходит довольно быстро. Поэтому для этой проблемы применим известный метод приближенного расчета концевых потерь.
Для вычисления неравномерного распределения индуктивных скоростей используются формулы, которые основаны на результатах классической вихревой теории несущего винта [22]. В этой работе учитывается первая гармоника неравномерности поля индуктивных скоростей. Считается, что для решения задач балансировки и динамики полета вертолета на большинстве режимов полета эта теория достаточна. Это удалось подтвердить и в численных исследованиях автора представленных в последующих главах.
Но основным недостатком методов численной реализации нелинейной теории пространственно-деформированных стержневых конструкций применительно к решению задачи аэроупругого несущего винта являются большие затраты машинного времени на поиск численного решения. Это связано в основном с условной сходимостью численного решения задачи упругих колебаний лопастей по времени (азимуту). Для решения таких задач нужны более эффективные методы и алгоритмы.
Во второй главе представлены эффективные методы и алгоритмы аэроупругого расчета винтов вертолета. Применение нелинейной теории пространственно-деформированных стержней крыльевого профиля позволяет получить систему интегро-дифференциальных уравнений аэроупругих колебаний лопастей несущего винта. Эту систему невозможно решить аналитически в том виде, в каком она была получена. Поэтому для приведения к матричной алгебраической форме применяются интегрирующие матрицы на основе интерполяции «напряженными» сплайнами (сплайн с растяжением).
Полет вертолета от взлета до посадки можно разделить на ряд чередующихся неустановившихся (маневр) и установившихся режимов полета. Интегрирование уравнений колебаний лопастей несущего винта на установившихся режимах полета вертолета по времени проводится при помощи методики, основанной на использовании разложения изгибных и крутильных колебаний лопасти в тригонометрический ряд Фурье. Идея применения коэффициентов разложения ряда Фурье рассмотрена в работе [37], как один из возможных путей решения уравнений вынужденных колебаний. Но проблема практического применения не была тогда решена. Главная особенность и преимущество данной методики состоит в наперед известной зависимости между прогибами, скоростями и ускорениями расчетных точек на лопасти. Эта особенность позволяет принципиально изменить путь поиска решения по азимуту и избавиться от итераций, связанных с нахождением скоростей и ускорений. Что позволяет на установившихся режимах полета существенно сократить время получения нагружения несущей системы вертолета в целом, не снижая при этом точности расчета. Данная задача была решена в рамках диссертационной работы автора на соискание степени кандидата наук. Дальнейшее развитие этой методики позволило включить в единую матрицу неизвестных и коэффициенты тяги НВ, по которым определяется поле индуктивных скоростей, а также решить задачу моделирования карданного винта по одной лопасти.
При моделировании неустановившегося движения лопастей НВ в качестве основы предлагается использовать обратный способ интегрирования по времени на основе кубической сплайн-интерполяции. Вывод зависимостей и алгоритмы решения уравнений неустановившегося движения также показаны в этой главе. Показаны результаты сравнительного анализа, и сходимости применяемых методик интегрирования по времени.
При исследовании задач динамики полета вертолета выполняются многократные расчеты для определения величин управляющих параметров, обеспечивающих безопасное управление вертолетом. Этот итерационный процесс может длиться недели. Поэтому проблема сокращения временных затрат при исследовании динамики полета вертолета достаточно актуальна. В данной главе также показано решение задачи имитационного моделирования нагружения несущего винта посредством искусственных нейронных сетей. f I Ij, *, 111 '.
Результаты исследований показали возможность применения такой модели в пилотажном тренажере вертолета.
Третья глава содержит результаты моделирования деформационных свойств торсиона вертолета «Ансат» производства ОАО «Казанский вертолетный завод». На основании проведенного анализа работ ряда исследователей можно сделать вывод, что на сегодняшний день конечно-элементные модели наиболее корректно обеспечивают полноценное моделирование упругого деформирования торсиона. Эти модели достаточно точно описывают НДС торсиона, особенно в части оценки зон концентрации напряжений. При этом они состоят из многих десятков или сотен тысяч элементов. Это требует значительных временных затрат, а при вычислении только перемещений концевой части торсиона, необходимых для решения задачи аэроупругих колебаний лопасти, такие модели избыточны. '.
В данной главе показаны способ и результаты построения матриц податливостей, позволяющих моделировать перемещения конца торсиона в зависимости от внешних сил, приложенных к нему. Построение матриц проведено при помощи упрощенной конечно-элементной модели разработанной в MSC. Nastran for Windows v4.0. Достоверность полученных результатов доказывает посредством сравнения с экспериментальными данными и результатами расчетов по модели других исследователей. Так же представлена методика и результаты и оценки влияния температуры внешней среды на деформационные свойства торсиона.
Полученные матрицы податливостей торсиона позволили провести сравнительное исследование уровня упругих деформаций шарнирного и бесшарнирного несущих винтов.
Важно заметить, что построение матриц податливости торсиона это долговременная и сложная научно-исследовательская работа. Для выполнения такого исследования необходимо конечно-элементная модель торсиона и комплекс экспериментальных работ по определению перемещений конца торсиона при различном наборе нагружения.
Моделирование шарнира таких сложностей не требует. Поэтому было проведено исследование возможности подобрать эквивалентный шарнирный винт, и решить вопрос можно ли обойтись без длительных расчетных и экспериментальных исследований.
Четвертая глава посвящена решению задачи пространственной балансировки и вибробалансировки. Проблема решения балансировки вертолета сравнительно слабо отражена в открытой печати. Помимо работ опубликованных в открытой печати, в промышленности были разработаны более эффективные методики расчета без традиционного разделения на продольную и поперечную балансировку, учитывающие конструктивные особенности несущих винтов с шарнирным креплением лопастей (наличие демпферов вертикальных шарниров, деформации лопастей и т. д.). Здесь следует отметить диссертацию А. Ю. Лисса, защищенную в 1974 году, где разработан метод корректировки балансировочных характеристик с учетом упругости лопастей.
В современных прогрессивных конструкциях несущего винта доля, А упругих деформаций в маховых движениях лопастей существенно повышается. Высокие маневренные характеристики вертолета с бесшарнирным винтом подразумевают более жесткую связь между управляющими параметрами и нагрузкой на несущем винте. То есть изменение управляющих параметров винта за очень короткий промежуток времени вызывает изменение по величине и направлению переменных нагрузок. С появлением таких конструкций несущей системы допущения, принятые для шарнирного винта неоправданны. Требуются новые подходы к решению задачи балансировки, аэроупругих колебаний лопастей и, соответственно, в прогнозировании нагрузок, действующих на агрегаты несущей системы вертолета. Решение же задачи определения переменных при оценке ресурса лопастей требует более точного определения нагрузок. Поэтому в прочностном расчете проводится уточнение нагружения несущей системы вертолета с учетом упругости лопастей и при этом корректировка балансировочных характеристик обычно не проводится. Таким образом, какой бы сложной не была математическая модель несущего винта при определении нагружения для оценки усталостной прочности, допущения, традиционно вводимые в математическую модель балансировки, не позволят получить результаты с желаемой точностью.
В первой части данной главы показана математическая модель пространственной балансировки одновинтового вертолета, предложенная в рамках кандидатской диссертации автора. Показаны новые результаты сравнения с летными данными, а так же результаты исследования влияния первой гармоники неравномерности поля индуктивных скоростей на балансировочные характеристики и нагружение лопастей несущего винта. Результаты сравнения доказывают возможность применения элементно-импульсной теории для поставленных перед программным комплексом задач. Исследованы вопросы применимости эквивалентного шарнирного винта взамен бесшарнирному при решении задачи балансировки вертолета. Которые доказывают эту возможность для решения задачи балансировки, но при этом не представляется возможным определить нагружение упругого элемента протяженного типа. Предложена методика решения задачи приведения численной модели к натурным характеристикам вертолета с дополнительным оборудованием. Показаны результаты сравнения с летными данными. Исследованы вопросы нагружения торсиона в полете в зависимости от различных факторов. Показано развитие данной модели вертолета в пространственную балансировку нового СВВП — «дисколета Павловых».
Во второй части главы предложена нелинейная модель пространственной балансировки с учетом периодических коэффициентов ряда Фурье, так как в реальности, на вал передаются гармонические составляющие нагружения кратные числу лопастей. Гармоники, которые при суммировании сил отдельных лопастей не уравновешиваются и формируют нагружение на втулке, называются проходными. Эти проходные гармоники вызывают вибрации вертолета. Однако это справедливо, в том случае, когда все лопасти идентичны и совершают одно и тоже периодическое движение. Источники этих нагрузок — след винта, также эффекты срыва и сжимаемости на больших скоростях полета, и для классического шарнирного винта они практически не зависят от деформаций изгиба лопасти [22].
С другой стороны с появлением в серийном производстве вертолетов с бесшарнирным несущим винтом возникли другие вопросы. В отличие от шарнирного у такого винта изгибающий момент передается на вал, следовательно, упругость лопасти и возможные отличия в деформационных свойствах одного из рукавов бесшарнирной втулки должны оказывать влияние на уровень вибраций такого вертолета. При этом упругость лопасти будет влиять только на проходные гармоники, а вот отличие в деформационных свойствах одного из элементов бесшарнирной «втулки вызовет появление низших гармоник. Но самое важное на данном этапе: — необходимость найти точное положение каждой лопасти в любой точке азимута с учетом динамики установившихся колебаний вертолета в пространстве, вызванных проходными гармониками нагружения. Это необходимо для моделирования динамики полета вертолета и движения лопастей на неустановившихся режимах полета. Поэтому в первую очередь для решения этих задач разработана математическая модель балансировки вертолета с учетом периодических коэффициентов.
В последней части главы представлена качественная оценка тех критических параметров и условий, которые оказывают значительное влияние на точность расчетов, в том числе и влияние учета упругости лопастей, правда, на больших скоростях полета, где влияние следа винта нет. Следует отметить, что для количественной оценки уровня вибраций следует иметь более совершенную модель аэродинамики несущего винта вертолета.
Пятая глава посвящена решению задачи динамики полета одновинтового вертолета с бесшарнирным несущим винтом. Моделированию динамики полета вертолета оснащенного классическим шарнирным винтом посвящено достаточно много работ [15, 24, 89]. Новый бесшарнирный винт поставил много новых задач перед исследователями.
В данной работе предложена численная модель движения вертолета в трехмерном Евклидовом пространстве и кинематические соотношения, в 1 которых отсутствует известный недостаток — составляющая -, что при собО п «угле тангажа — приведет к делению на ноль. В качестве результатов доказывающих достоверность моделирования неустановившегося движения представлен ряд исследований, на всем протяжении которых осуществляется сравнение с летными данными.
Исследуются способы учета нагрузок, создаваемых несущим винтом, так как в большинстве случаев для расчета динамики полета вертолета применяют допущение, что в каждый момент времени упруго-маховое движение лопастей, силы и моменты несущего винта соответствуют их мгновенным значениям. В этом случае переходные процессы движения лопастей винта не учитываются, а вводится некий коэффициент торможения по управлению. Этот коэффициент обычно находится посредством сравнения с результатами летных испытаний уже созданного вертолета. А на этапе проектирования вертолета он обычно принимается на основании опыта разработчиков или данных аналогичного вертолета, которые обычно засекречены. Особенно этот вопрос интересен для нового вертолета с бесшарнирным несущим винтом с высокими моментными характеристиками. Проверка результатов подбора эквивалентного шарнирного винта показала, что возможно использовать эквивалентный бесшарнирному шарнирный винт при решении задач динамики полета вертолета на неустановившихся режимах полета.
Сертификация, то есть получение сертификата типа летательного аппарата является необходимым условием выхода на серийное производство и эксплуатацию вертолета. На этом этапе проверяются конструктивные решения, прочностные расчеты и т. д. заложенные на этапах проектирования вертолета. При этом следует заметить, что в некоторых, но очень важных пунктах норм летной годности не заложена методическая база по решению поставленной задачи. Поэтому в ходе научно-исследовательских задач были решены также и задачи по методическому обеспечению решения задач сертификации:
• определение максимально и минимально возможных перегрузок;
• возможности обеспечения безопасной авторотационной посадки;
• и ряда других не представленных в этой работе.
При проведении расчетов динамики полета вертолет с бесшарнирным винтом возникли трудности с пилотированием модели. Вследствие достаточно небольших значений моментов инерции вертолета подбор I управления по циклическому шагу занимал достаточно длительное время. Ведь даже небольшое отличие параметров управления от необходимых значений приводило с течением времени к значимому отклонению от заданной траектории. Это приводило к неоправданным затратам времени, что противоречило основным целям разработанного программного комплекса. Поэтому был разработан алгоритм автоматизации рабочего места оператора программного комплекса, т. е. реализована возможность полета по заданной траектории. Также исследован вопрос влияния второй гармоники управления на динамику колебательного движения вертолета в пространстве.
Практическая значимость. Программный комплекс балансировки и динамики полета вертолета с бесшарнирным аэроупругим несущим винтом нашел свое применение при проектировании, испытаниях и сертификации вертолета «Ансат», а полученные результаты одобрены специалистами сертификационных центров РФ. В ходе этих работ автором был проведен ряд научно-исследовательских работ. При этом оформлено около 40 научно-технических отчетов.
Основные результаты диссертации изложены в научных статьях опубликованных: в журналах реферируемых ВАК [38, 39, 40, 41, 42, 43, 44, 45], в зарубежных научных изданиях [46, 47, 48] и докладывались на форумах Европейского [49, 50, 51] и Российского вертолетных обществ [52, 53, 54, 55], а также на других научно-технических конференциях [56, 57, 58, 59, 60,61,62, 63,64, 65].
Общие положения.
Системы координат.
Земная О0×0у0г0 — оси и начало этой системы координат связаны с Землей и выбираются из условий задачи.
Нормальная земная — отличается от земной тем, что начало координат центре масс вертолета.
Связанная Охуг — с началом в центре масс вертолета. Ось Ох направлена вперед параллельно так называемой строительной горизонтали вертолета. Ось Оу перпендикулярна Ох и направлена вверх, а Ог образует правую систему координат.
Для перехода из Ох ув связанную систему координат используется матрица перехода Эйлера вида ь] =.
В матрице (0.1.1) углы ц/, 0, у — углы последовательных переходов, определяющие ориентацию вертолета в пространстве.
Угол рыскания ц/ - угол между осью Ох£ нормальной систем координат и проекцией продольной оси Ох на горизонтальную плоскость 0Уг.
Угол тангажа 0 — угол между продольной осью и горизонтальной плоскостью Ох у нормальной системы координат. сое Особц/.
— соз^вт^.
— втдсоэ^созу + со8ц/зту + сое в соэ у 8И11|/8ту + 8Н1 $ 8Н11|/С08у втОсовувту + +эту сое у.
— созЗвту соеусову—ятОвцнувту.
0.1.1).
Угол крена у — угол между поперечной осью и осью Огё, смещенной в положение, соответствующее нулевому углу рыскания.
Г'.
Рис. 1. Системы координат вертолета: — вид сбоку.
Невращающаяся Оътхпунгн с началом в центре втулки винта. Обычно вал несущего винта наклонен относительно фюзеляжа на угол заклинения фзак, тогда матрица перехода.
С08фзак фзак О.
— 8Шфзак СОБ фзак О.
О 0 1.
0.1.2).
Вращающаяся Опхвувгв — оси этой системы координат получены путем поворота системы координат Овтхиуиги на угол азимутального положения лопасти ц/н. Положительное направление азимутального вращения по часовой стрелке при виде сверху.
С08Ц/" 0 8Н1)/н ч]=.
0 1 0 -8И1Ц/Н 0 С081]/н.
Угол азимута НВ отсчитывается от хвостовой балки, т. е. начало отсчета угла сдвинуто на 90 градусов. В этом случае необходим еще один поворот на 90 градусов.
0.1.3).
О 0 1 О 1 о -10 0.
0.1.4) у=180.
90°.
Рис. 2. Системы координат вертолета: — вид сверху.
Маховая 0Kxbybzb — это система координат с началом Ок, находящимся в центре жесткости комлевого сечения с координатами хк, ук, zK. Оси повернуты на углы взмаха р, отставания г) и осевого поворота С, комля лопасти относительно вращающейся с.к. sin^cos (3+ sin^sinp COS^SInriSInPCOSI^SInr|COSp cos^cospCOS (^SInp +.
— sin P sin rising +sin^sinrjcosP sint] -sinpcosri cos P cost].
K = cos^cosrjsinpcosri.
0.1.5).
Связанная с деформированным элементом лопасти 0 $х5у$гсистема координат с началом 0Е в центре жесткости каждого сечения деформированной лопасти, и ее оси направлены по главным центральным осям сечения и повернуты относительно Окхьуьгь на углы последовательных поворотов фь Ф2, фзвт^сов^ + сое (р^ сое (р2 -$т (ръсо%(рг вШ (ръ БШ (рх -+сое (ръ вт (р2 вт <�рх — сое (ръ БШ (рг сое (рх сов (ръ сов фх — совз вт (рх +.
— вт срх вт <�рг вт (ръ ч-вт (ръ вт (рг сое (рх.
— ъъ (рхо.оъ<�рг сов^сов^.
0.1.6).
Связанная с центром масс сечения лопасти Отхтут2тсистема координат с началом От в центре масс каждого сечения деформированной лопасти, а ее оси параллельны 05×5у5г5.
Рис. 3. Лопастные системы координат.
Основные обозначения.
Физические параметры: т, 0 — масса и сила тяжести вертолетатя — погонная масса элемента лопасти- § - ускорение свободного паденияМ.
X /гУ ^Х2.
— 3 3 -3 ух у уг.
— 3 -3 3 гх гу г тензор инерции вертолета;
Л" л] - тензор инерции сечения лопасти. Параметры пространственного движения вертолета: о{К>Уу>К)>{Кг>Уу8>^) -скорость поступательного движения центра масс вертолета и ее проекции на оси связанной и нормальной земной системы координат;
Ух, Уу, У2 — локальная производная по времени от скорости движения центра масс в осях связанной (неинерциальной) системы координат;
П0-угловая скорость вращения связанной системы координат относительно нормальной земной системы координат;
— ускорение центра масс вертолета и ее проекция на оси нормальной земной (инерциальной) системы координат;
Рх, Ру, Рг и Мх, Му, Мг — проекции главного вектора внешних сил Р и момента М на оси связанной системы координат;
Хт9Гш, г" и Мхпл, Му ш, Мгш — проекции аэродинамических сил, создаваемых планером вертолета, на оси связанной системы координат;
ХНД, 2Н и Мт, Му1{, Мт — проекции сил и моментов, создаваемых несущим винтом вертолета;
Хри М^М^Мцпроекции сил и моментов, создаваемых рулевым винтом вертолета.
Управляющие параметры:
Фн — общий шаг несущего винта;
0р02 — продольный и поперечный угол циклического шага несущего винта;
Фр — шаг рулевого винта.
Разные параметры:
В — коэффициент концевых потерьпнл, прл — число лопастей несущего и рулевого винтов;
Ки, Кр — количество гармоник разложения в ряд Фурье по азимуту несущего и рулевого винтовсон, сор — угловые скорости вращения несущего и рулевого винтов.
Сокращения: НВ — несущий винтРВ — рулевой винт;
СУ — система управления общим и циклическим шагом;
АП — автомат перекоса;
ИНС — искусственная нейронная сеть.
Индексы: н — несущий винтр — рулевой винтпл — планер.
Заключение
.
1. На основе геометрически нелинейной теории пространственно-деформируемых стержней крыльевого профиля выведен вариант уравнений упруго-маховых колебаний лопасти бесшарнирного несущего винта с учетом произвольного движения вертолета в пространстве. В данных уравнениях математически и физически корректно учтены граничные условия, определяемые бесшарнирной втулкой винта.
2. Разработаны эффективные методики и алгоритмы численного решения уравнений упруго-маховых колебаний лопасти. В этих способах, в отличие от известных, предложено:
— интегрирование по длине проводить при помощи интегрирующих матриц на основе сплайна с растяжением;
— интегрирование по времени построено на сочетании методики разложения в ряд Фурье и обратного способа на основе кубического сплайна.
Проведен сравнительный анализ и исследована сходимость применяемых методик. Показана их высокая эффективность.
Предложен способ имитационного моделирования нагружения бесшарнирного несущего винта, основанный на алгоритмах искусственных нейронных сетей. Исследована задача влияния топологии нейронной сети на точность моделирования нагружения. Полученные результаты исследований показывают возможность применения такой модели в пилотажном тренажере вертолета.
3. Разработана упрощенная конечно-элементная модель упругого элемента (торсион) бесшарнирной втулки МБС. Каз^ап. Посредством сравнения с экспериментальными данными и результатами расчетов другой КЭ модели показана достоверность моделирования деформационных свойств торсиона. Проведены расчетно-экспериментальные исследования деформационных свойств торсиона. Выявлены зависимости перемещений и углов наклона конца торсиона от нагружения и температуры окружающей среды. На основе этих результатов построена высокоскоростная имитационная модель торсиона, позволяющая выполнить прямое включение в модель аэроупругого несущего винта.
4. Предложена модель пространственной балансировки одновинтового вертолета с бесшарнирным аэроупругим несущим винтом. Путем сравнения с летными данными показана достоверность предложенной модели и корректность применения принятой модели поля индуктивных скоростей, создаваемых несущим винтом вертолета. Доказана необходимость учета упругости лопастей при расчете нагружения упругого элемента протяженного типа в полете. Показана возможность расчета балансировочных характеристик и динамики неустановившихся режимов полета с применением эквивалентного шарнирного винта (увеличением разноса ГШ). Показано, что в этом случае не представляется возможным достаточно точно вычислить нагружение упругого элемента протяженного типа. Оценено влияние составляющей деформации сдвига торсиона на балансировку вертолета. Разработана модель балансировки одновинтового вертолета с периодическими коэффициентами. Исследована сходимость по числу учитываемых гармоник. Оценено влияние упругости лопастей на величину первой проходной гармоники момента тангажа, создаваемого бесшарнирным несущим винтом вертолета.
5. Предложен способ решения уравнений динамики полета вертолета. В этих кинематических соотношениях отсутствует критическая точка. Проведено исследование способа моделирования динамики полета вертолета, когда в каждый момент времени движение лопастей, силы и моменты НВ соответствуют их мгновенным значениям, или когда моделируется движение каждой лопасти по времени. Показано влияние переходных процессов. Представлены результаты решения задачи определения максимальных и минимальных перегрузок при сопровождении сертификационных испытаний легкого многоцелевого вертолета. Выполнено исследование возможности безопасной посадки на авторотации.