Методика расчёта напряжённо-деформационного состояния волноводно-распределительных систем космических аппаратов Текст научной статьи по специальности «Строительство и архитектура»

Journal of Siberian Federal University. Engineering & Technologies 2 (2012 5) 150-161

УДК 621.372.83.001.24

Методика расчёта

напряжённо-деформационного состояния волноводно-распределительных систем космических аппаратов

П.Н. Сильченкоа*, И.В. Кудрявцев3, М.М. Михнёв6, В.Н. Наговицин6

а Сибирский федеральный университет, Россия 660041, Красноярск, пр. Свободный, 79 б ОАО «Информационные спутниковые системы

им. М.Ф. Решетнёва» Россия 662972, Красноярский край, Железногорск, ул. Ленина, 52 1

Received 4.04.2012, received in revised form 11.04.2012, accepted 18.04.2012

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

Ключевые слова: космический аппарат, волноводно-распределительная система, тонкостенные элементы, неосесимметричное поперечное сечение, прочность, напряженно-деформированное состояние, метод расчета, программное обеспечение.

Введение

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

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

* Corresponding author E-mail address: [email protected], [email protected]

1 © Siberian Federal University. All rights reserved

Функционально-эксплуатационные характеристики волноводно-распределительной системы КА в совокупности с наземными антенными системами во многом влияют на качество и надежность работы всего наземно-космического комплекса в целом.

Управление и функционально-эксплуатационными характеристиками ВРС с одновременным влиянием на увеличение срока активного существования всего КА возможно обеспечивать в некоторой степени путем достижения стабильности и точности геометрических размеров и напряжённо-деформированных состояний отдельных элементов и конструкций антенно-фидерных систем в целом. Все элементы волноводно-распределительных систем должны в процессе изготовления (сборка, монтаж, отладка, регулировка и т.д.), испытаний, наземных транспортных операций, выведения на орбиту и раскрытия трансформируемых элементов выдерживать все возникающие комбинации статических и динамических нагрузок с обеспечением глобальной и локальной прочности без появления остаточных деформаций.

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

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

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

Постановка задачи

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

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

Наиболее распространенные виды элементов представлены на рис. 2: прямые (рис. 2, а), криволинейные тонкостенные элементы, имеющие постоянный (рис. 2, б,в) или переменный (рис. 2, г) радиусы кривизны в любой из координатных плоскостей. Эти тонкостенные эле- 151 -

Рис. 1. Пример участка волнов одно-распределительной системы: а- прямой элемент; б - криволинейный элемент с переменным радиусом кривизны; в - криволинейный элемент с постоянным радиусом кривизны; г - соединительная муфта; д - фланец; е - гибкая секция; ж - промежуточная опора

Рис. 2. Основные элементы волноводно-распределительной системы: а - прямой элемент; б -криволинейный элемент с переменным радиусом кривиены; в - криволинейные элементы с постоянным радиусом кривизны; г - соединительная муфта; д - фланец

менты соединены между собой при помощи муфт и фланцев (рис. 2, д,е). Набором таких различных тонкостенных прямых и криволинейных элементов, соединительных муфт и фланцев моделируется конструкция любого (рис. 1) участка ВРС применительно к рассматриваемому КА.

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

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

Рис. 3. Варианты конструктивного исполнения опор для крепления участков ВРС: а - стержневая опора; б - пластинчатая опора

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

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

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

На этапе изготовления неизбежно возникают случайные воздействия, приводящие к появлению технологических погрешностей, которые нарушают расчетную геометрию соединения между каждым элементом участка и участками ВРС в целом. В результате, на этапе монтажа ВРС требуются дополнительные усилия для совмещения точек крепления каждого участка с опорами, что приводит к появлению в отдельных элементах внутренних монтажных напряжений, которые в сочетании с действующими нагрузками на этапе эксплуатации могут вызвать нарушение требуемых условий локальной и глобальной прочности ВРС [3]. Частично накопленные технологические погрешности сборки можно компенсировать установкой гибкой секции в конце (рис. 1, д) рассматриваемого волновода.

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

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

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

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

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

Согласно [4] для расчета напряженно-деформированного состояния (НДС) «базового» прямого элемента используем подход линейной теории тонких оболочек. Однако применение теории оболочек для муделирования тонкостенных прямо}'гольного коробчутого сечения элементов ВРС невозможно (причина указана выше). Следовательно, невозможно и представление математической модели системой дифференциальных уравнений, связанной общими граничными условиями, соответствующими реальной конструкции элемента ВРС с учетом резких переходов (складок).

Решение задачи

1

N

Рис. 4. Модель криволинейного элемента ВРС (набор А'прямых элементов под углами у)

Этот существенный недостаток не позволяет напрямую использовать линейную теорию оболочек для решения задач с тонкостенными конструктивными элементами ВРС, поэтому «базовый» прямой элемент ВРС представим в виде составной пространственной конструкции, состоящей из отдельных пластинок, как частный случай оболочки при радиусах кривизны пластинок Я,,;—» (рис. 5).

Рассматривая отдельную тонкую пластинку в ее локальной сиитеме координат а, в, z (рис. 5), воспользуемся геометрически иинейными уравнениями [4] , которые с учетом параметров Ламе А=ы=1,радиусах кривизны Я;=Я2—у и толщине пластииы к, примтут вид

И а-Д 9 ^^ М =й н о;

^ д2 о, Ж а 0.

Подс тавив в (1) вместо усилив!! Иъ Ы2, Б их выражения в функциях соответствующих на-пряже ний ф = ф(а, =) д=я каждой рассматрива.мой пластй1ны:

2 2 2 ^ —г к; ы2 —г л; 5 -л (2)

получим дифференциальное уравнение ¿4-го порядка для 1-й пластины в виде

Ек

(3)

4

4

4

2

2

2

4

2

С учетом прогиба пластины со = со(а, (3) выражения для изгибающих и крутящих моментов Ы2, М2, Н будут и меть вид

2 2 2 2 2 м1 в — — ; м2 о— — ; н (1 -, (4)

где; Гс - цилиндрическая жест ко сть пласт ины; ц - коэффициент Пуасс она.

Подставив выр ажения (4) для моментов в последние дв а уравне ния системы (1), получим дифференциальное уравнение 4-го порядка для /-й пластины в виде

4 4 4 ,2 2 2

—г —Г д Г -- Г ■ (5)

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

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

- 155 -

Рис. 5. Расчетная схема «базового» прямого элемента

в конечном итоге представляют собой! систему из 8-ми нелинейных дифференциальных уравнений 4-го порядка в частных про изводных, имею щую вид

4 4 4^2 2 2

N. -т- -— М,-^ д,-д,-дп /г-^ ;

(6)

2 2 '

2—з Нг —т —т -- —Г Ч,— Ч,— Л—

4 4 4 2 2 2 2

-¡- -Г^-Г -г Е)1 -

где i=1,2,3,4.

Полученная система (6) включает четыре подсистемы уравнений (3) и (5) для каждой из четырех отдельных пластин прямого элемента, каждая из которых описывает НДС соответствующей пластины. Все подсистемы связаны между собой зависимостями, отражающими граничные условия в местах со единения пластин между собой (рис. 5).

Например, на линиях соединения 1-й и (1+1) -й пластин граничные условия, выраженные через силовые факторы, (Зуд}' т и меть вид

'!/», Ь,

Ql

в

А ъ,

а(1 1)

!А , Ь, .

; N

'А ъ,

в,

РЧ II

!А, ь,

(7)

'I,, г>,,' ь, Ъ,;

Если выразит1> эти грани чные усло вия через функцию напряжений ф(а, в) и прогиб ю(а, в), то условия примут вид

А —г ь

с, /!

О

^II 3

I

с, , /г

(2 )-

(2

3 3 3

г г ^ /9 \_г У

2

11 11

(8)

где Б, и С,- - цилиндрические жесткости отдельных пластин в составе прямого элемента.

Полученная глобальная система дифференциальных уравнений (6) для четырёх взаимосвязанных пластин (рис. 5) совместно с граничными условиями (8) является сложной, и в имеющейся литературе не обнаружено информации по их аналитическому решению. Имеются частные случаи приближенного решения в этой системе лишь одной подсистемы, в которой отсутствует динамический член, и только для частных случаев условий закрепления и нагружения пластинки [6, 7]. Аналитическое решение рассмотренной задачи расчёта всей тонкостенной конструкции волноводно-распределительной системы космического аппарата весьма сложно, и в настоящее время мы не располагаем информацией о возможных методах её решений.

Использование численных методов решения системы (6) не всегда обеспечивает необходимую точность результатов [8]. В частности, при применении метода конечных элементов (МКЭ) результаты расчетов будут очень зависеть от типа и количества используемых конечных элементов (КЭ). В результате при расчете участков ВРС методом конечных элементов толщина стенки конструкционных элементов будет сопоставима с размером самого конечного элемента, а для получения результатов необходимой точности нужно как минимум два слоя конечных элементов. С учетом значительной протяженности тонкостенной конструкции ВРС это требует применения очень мелкой сетки разбиения, что согласно [8] может приводить к значительным вычислительным погрешностям и потере точности получаемых результатов. Использование же грубой сетки элементов приводит к неточным или же некорректным результатам расчета, особенно в динамической постановке.

Для выхода из данной ситуации и получения решения с требуемой точностью нами разработана методика [2], согласно которой любой исследуемый участок волноводно-распределительной системы (см. рис. 1) в глобальной постановке моделируется пространственной стержневой конструкцией с эквивалентными условиями нагружения и закрепления (рис. 6). Соединительные муфты и фланцы на данном этапе считаются отсутствующими, но будут учтены на следующем этапе уточнённого решения для локальной области участка ВРС. Число стержней при моделировании участка выбирается в зависимости от требуемой точности расчета и может значительно превышать число составляющих его тонкостенных элементов.

5

2

2

2

1 i 1

1 I 1

2 3

Рис. 6. Стержневая модель участка волноводно-распределительной системы

Для каждого / -го стержня строится матрицт жесткости к,- [9] , которая является частью глобальной матрицы жесткости к всей етержневой констрскции ВРС. В результате получаем систему линейных уравнений равновесия стержневой модели ВРС в виде

ку = Г

(9)

где V - вектор перемещений узлов конечно -элементное стержней; Г - вектор внешних нагрузок.

Решение системы (9) выпооняется с иипользованием метода Холецкого но основе ЬЮ-разложения с получением вектора перемещений V для всех узлов модели, применяя который определяют усилия и, следовательно, напряжения и деформации в любой точке всей конструкции ВРС. Используя по лученные резуньтаты расчета, по известным теориям прочности определяют место расположериэ опас ного сечения (рис. 6).

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

Для проведения уточнённого расчета выделенной локальной области опасного сечтния интересующий участок с опасным сечением транслируется в виде твердотельной модели в ППП Ansys (Nastran и др.) с соответствующими граничными условиями в виде векторов усилий и перемещений, полученными на этапе расчёта ВРС как пространственной стержневой конструкции.

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

Рис. 7. Пример расчета стержневой модели участка ВРС

приложению соответствующих краевых граничных условий. В программе для этого моду ля используется вну тренний язык программирования квекретного ППП. Например, для передачи выделенного участка в ППП Л№У8 разработунное программное обеспечение использует внутренний язык программирования APDL этого пакета [10].

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

Все процедуры, необходимые для задания геометрии, свойств материала, условий закрепления и нагружения, расчета участуу ВРС как стержневой конструкции, а также выделения локального участка и его трансляция в виде лвордотельной модели в ППП Лnsys (Nas8аan и др.) с соответствующими граничными условиями автоматизированы и оформлены в виде законченного программного продукта.

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

Заключение

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

Рис. 8. Расчетные модели выделенного локального участка! ВРС: а - КЭ-модель локального участка; б - результаты расчета локального участка

Созданные методики и программное обеспечение использованы в ОАО «Информационные спутниковые системы имени академика М.Ф. Решетнёва» при проектировании участков ВРС КА серии «Экспресс», «ГЛОНАСС»», «Луч-5А»с. Результаты расчетов позволили об основать уменьшение толщины стенки элементов с 1,Л дс 0,6 мм и о»еапечение условий прочности и жесткости при снижении массы ВРС на »0-30 %.

Работа выполнена по результатам НИР №20116, ГР №01200906434, г. Красноярск, »011.

Список литературы

[1] Сильченко П.Н., Корчагин А.И., Михнев М.М. и др., Пат. № 2317184 РФ // Б. И. 2008. Бюл.

№ 5.

[2] Сильченко П.Н., Михнев М.М., Кудрявцев И.В. и др. // Проблемы машиностроения и надежности машин. 2012. №1. С. 91-95.

[3] Сильченко П.Н., Кудрявцев И.В., Михнев М.М». // Технология машлностроения. 2011. № 2. С. 26-30.

[4] Новожилов В.В. Линейная теория тонких оболочек. Л.: Политехника, 1991. 656 с.

[5] Сильченко П.Н., Михнев М.М., Кудрявцев И.В. и др. // Труды Международной конференции по дифференциальным уравнениям и динамическим системам. Суздаль, 2010. С. 172-174.

[6] Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки: пер. с англ. изд.3. М.: Эдиториал УРСС, »009. 640 с.

[7]. СусотеринМ. В. . Дис. ... д-ратехн . наук. СПб., 2010. 300 а.

[8] Галанин М.П. Методы численного анализ. матоматичсских моделей. М.: МГТУ им. Н.Э. Бау мана, 2012. 591 с.

[9] Przemieniecki J.S. Theory of Matrix Structural Analysis. McGraw-Hill, New York. 2010. Р. 480.

[10] A nsys Parametric Design Language Guide (APDL) // Ansys Inc. 2009. P. а0.

[11] Сильченко П.Н., Михнев М.М., Кудрявцев И.В. и др. // Заключительный отчет о НИР 20116 ПИ СФУ; рук. П. Н. Сильченко. Красноярск, 2011. 239 с. № ГР 01.2.00304354.

- 160 -

Method of Stress-Deformation Distribution Computation for Waveguide Spacecraft Systems

Petr N. Silchenkoa, Ilya V. Kudryavceva Mikhail M. Mikhnevb and Vasiliy N. Nagovitsinb

a Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041 Russia b JSC "Information Satellite Systems" Reshetnev Company" 52 Lenin Str., Zheleznogorsk, Krasnoyarsk region, 662972 Russia

The problem question of providing deflected mode of waveguide distributive systems of communication satellite on one of operational stages is considered. The differential equations of making elements are received and the technique of their calculation on the static strength, realized in the software is developed for the analysis deflected mode of separate regions of waveguide distributive systems of communication satellite.

Keywords: communication satellite, waveguide distributive system, thin-walled structure, asymmetrical cross-section, strength, stress state, deflected mode, design procedure, software.