COMPARATIVE CALCULATION OF THE LOW FREQUENCY COMPONENT OF DYNAMIC OVERLOADS FOR ONBOARD SPACE VEHICLE EQUIPMENT

Full Text

ВВЕДЕНИЕ В состав современного космического аппарата (КА) входит множество систем служебного и целевого назначения, обеспечивающих функционирование изделия и выполнение им полётного задания. Системы представляют собой навесное оборудование (НО), состоящее из приборов, агрегатов, датчиков. НО крепится к корпусу изделия с помощью крепёжных элементов (кронштейнов, посадочных плоскостей, рам) как с внешней стороны, так и внутри соответствующих отсеков. В процессе эксплуатации изделия, приборы и агрегаты испытывают различные виды нагрузок, изменяющихся с течением времени (вибрации, удары). Эти процессы включают в себя статическую и динамическую составляющую. Определение последней является сложным и ответственным этапом расчётов. В частности, динамическое воздействие обычно выступает определяющим для бортовой аппаратуры. Оно включает в себя свободные и вынужденные колебания. Цель данной статьи - проведение сравнительного анализа определения динамического поведения приборов и агрегатов космического аппарата, при воздействии на неё внешних возмущающих факторов. ВЫСОКОЧАСТОТНЫЕ И НИЗКОЧАСТОТНЫЕ КОЛЕБАНИЯ Вибрационная устойчивость - важная задача при проектировании КА. Действие вибраций (высокочастотных процессов) на этапе полёта в составе ракеты космического назначения (РКН) негативно сказывается на вибропрочности бортовой аппаратуры, закреплённой на корпусе изделия при помощи кронштейнов и посадочных плоскостей [1]. Источником вибраций преимущественно выступают: двигательная установка в процессе работы, колебания топлива в магистралях, работа силовых гироскопов системы управления КА. Вибрации и удары приводят к образованию усталостных трещин в кронштейнах, ослаблению болтовых соединений и образованию других дефектов, что может привести к разрушению. Помимо описанных вибраций на аппаратуру КА воздействуют низкочастотные корпусные процессы, совместно они показаны на рис. 1. Данный вид нагрузки даёт наибольший отклик конструкции в моменты старта, максимального скоростного напора, разделения ступеней и др. Каждая из перечисленных составляющих оказывает значительное воздействие на поведение конструкции КА на этапе выведения. Смоделировать общую динамику КА, максимально приближенную к реальной - главная задача инженера-конструктора. ДИНАМИЧЕСКИЙ РАСЧЕТ КОНСТРУКЦИИ Для выполнения расчёта строится конечно-элементная модель (КЭМ). Она должна учитывать многие нюансы конструкции объекта исследования, так как от этого напрямую зависит точность полученных данных. Для расчёта нагрузок на установки агрегатов, приборов используются следующие исходные данные: а) массовые, центровочные и инерционные характеристики установок; б) узлы крепления основных установок на силовом корпусе изделия (тип закрепления); в) жесткостные характеристики установок агрегатов, приборов в виде перемещения центра масс приборов, агрегатов под действием собственного веса в любом направлении. В качестве объекта исследования в данной статье выступает прибор системы управления движением (СУД) малого космического аппарата (МКА) «Аист-2Д» [2]. Задача состоит в определении низкочастотной динамической перегрузки, действующей в центре масс прибора при расчетном случае «максимальный скоростной напор» на этапе полета в составе РКН. Предварительно проведен анализ амплитудно-частотных характеристик (АЧХ) конструкции МКА «Аист-2Д» (рис. 2) [3]. С помощью программного комплекса Femap with NX Nastran вычислены максимальные перегрузки в центре масс прибора (рис. 3). Использовался анализ переходных процессов [4]. Результат отражен в табл. 1. МЕТОД ОЦЕНКИ НАГРУЗОК НА БОРТОВУЮ АППАРАТУРУ ПРОГРАММОЙ «AGENA 4» Систему прибор-кронштейн называют установкой и, для упрощения математических расчётов, представляют как «массу на пружине» (рис. 4). Каждая из таких систем имеет свои динамические характеристики и определённым образом оказывает воздействие на всё изделие. В данном методе, для расчёта нагрузок на установки приборов и агрегатов, используется программа, написанная на языке Fortran. Она позволяет определить коэффициент динамичности, связывающий значения ускорения в центре масс установки и ускорения на корпусе, в месте крепления кронштейна. Для её использования необходимо знать, помимо частоты и перегрузки, которые реализуются на корпусе (в месте крепления кронштейна), динамические характеристики рассматриваемой установки. Собственные частоты колебаний рассматриваемой установки, полученные по результатам анализа АЧХ, приведены в табл. 2. Система координат совпадает с системой координат изделия. Работа программы основана на следующих принципах. Для расчётной схемы «корпус - прибор» (рис. 4) коэффициент динамичности определяется из решения дифференциального уравнения, описывающего колебания одномассовой упругозакрепленной системы с колеблющимся основанием, как отношение максимального ускорения на массе к максимальному ускорению основания системы: (1) где ym- перемещение массы, wm- круговая частота собственных колебаний системы, wm=2pλуст, dm - логарифмический декремент собственных колебаний системы, dm = dуст i, yосн(t)- перемещение основания системы, определяемое заданным ускорением основания системы: , (2) λс, δуст - частота и логарифмический декремент собственных колебаний установки прибора, агрегата; λki , nki(x)- частота колебаний корпуса изделия и амплитуда перегрузки (ускорения) в месте крепления прибора, агрегата по длине корпуса изделия на i-ой частоте колебаний корпуса изделия; ,(3) где dki - демпфирующий коэффициент, для учёта нестационарности колебаний корпуса изделия [5]. Графическое представление колебаний системы изображено на рис. 5. В структуре программы приняты следующие обозначения: Nstrok - количество обрабатываемых строк; wk - частота колебаний корпуса (Гц); Dk - логарифмический декремент корпуса; Pk - ускорение на корпусе (м/с2); wm - частота колебаний установки прибора (Гц); Dm - логарифмический декремент установки прибора; Nzikl - количество циклов корпусных колебаний; ndt - количество шагов за период колебаний; dPkk - изменение ускорения корпуса в квадрате. Исходные данные и результаты решения отражены в табл. 3 - 4. Каждая из осей соответствует системе координат изделия (X - вертикальное направление, Y - поперечное, Z - дополняет систему до правой тройки векторов). ОБСУЖДЕНИЕ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ Мы получили низкочастотную составляющую (λki на рис. 5) максимального динамического отклика установки прибора. Для получения суммарной динамической нагрузки необходимо вычислить вибрационную составляющую (λc на рис. 5). В табл. 5 приведено сравнение результатов решения двумя методами. Значения показывают величину перегрузки, реализующейся в центре масс прибора в каждом из направлений. Видно, что показатели имеют пренебрежимо малое расхождение по величине. ВЫВОДЫ Проведен расчёт низкочастотных нагрузок на прибор ориентации по Земле МКА «Аист-2Д» с помощью программы конечно-элементного моделирования и с использованием программы расчёта динамики бортовой аппаратуры. Результаты показали высокую сходимость, что подтверждает возможность использования программы. Изложенная программа может быть использована при расчётах нагрузок на бортовую аппаратуру не только ракетно-космической техники, но и авиационной техники, автомобилестроения. Она позволяет с достаточной степенью точности спрогнозировать динамическое поведение конструкции в различных расчётных случаях. Рис. 1. Высокочастотные и низкочастотные процессы Рис. 2. Прибор ориентации по Земле Рис. 3. Максимальный отклик конструкции в направлении оси Y Таблица 1. Максимальная низкочастотная составляющая динамической перегрузки в центре масс прибора (использование программного обеспечения Femap/NX Nastran) Рис. 4. Расчётная схема установки прибора (агрегата) Таблица 2. Собственные частоты колебаний установки Рис. 5. Колебания установки прибора (агрегата) по времени Таблица 3. Входные данные Таблица 4. Выходные данные Таблица 5. Сравнение результатов

About the authors

Aleksandr Alekseevich Igolkin

Samara National Research University named after Academician S.P. Korolyov

Email: igolkin97@gmail.com
Doctor of Science (Engineering), Associate Professor, Professor, Department of Power Plant Automatic Systems

Evgeny Vladimirovich Shakhmatov

Samara National Research University named after Academician S.P. Korolyov

Email: shakhm@ssau.ru
Corresponding Member of the Russian Academy of Sciences, Doctor of Science (Engineering), Professor, Research Supervisor of the University

Artem Andreevich Popkov

Samara National Research University named after Academician S.P. Korolyov

Email: artiem.popkov.95@mail.ru
Postgraduate Student, Department of Power Plants Automatic Systems

Aleksandr Gennadievich Filipov

Samara National Research University named after Academician S.P. Korolyov

Email: iskander-filipov@yandex.ru
Postgraduate Student, Department of Power Plants Automatic Systems

References

Supplementary files