Применение метода синтеза форм для расчета колебаний космического летательного аппарата
Учебное пособие по курсу «Динамика конструкций космических летательных аппаратов»
Покупка
Новинка
Тематика:
Космический транспорт
Год издания: 2009
Кол-во страниц: 16
Дополнительно
Рассмотрены вопросы расчета собственных колебаний конструкции космического летательного аппарата с применением метода синтеза форм колебаний в виде метода жестких границ по Крейгу — Бэмптону.
Для студентов, обучающихся по специальностям «Ракетостроение» и «Космические аппараты и разгонные блоки».
Скопировать запись
Фрагмент текстового слоя документа размещен для индексирующих роботов.
Для полноценной работы с документом, пожалуйста, перейдите в
ридер.
Московский государственный технический университет имени Н.Э. Баумана С.Н. Дмитриев, И.Ю. Калугин, О.Н. Тушев ПРИМЕНЕНИЕ МЕТОДА СИНТЕЗА ФОРМ ДЛЯ РАСЧЕТА КОЛЕБАНИЙ КОСМИЧЕСКОГО ЛЕТАТЕЛЬНОГО АППАРАТА Рекомендовано редсоветом МГТУ им. Н.Э. Баумана в качестве учебного пособия по курсу «Динамика конструкций космических летательных аппаратов» М о с к в а Издательство МГТУ им. Н.Э. Баумана 2 0 0 9
УДК 629.7.01(075.8) ББК 39.56 Д 532 Ре це нзе нт ы: О.В. Кузнецов, В.И. Никитенко Дмитриев С.Н., Калугин И.Ю., Тушев О.Н. Применение метода синтеза форм для расчета колебаний космического летательного аппарата: Учеб. пособие по курсу «Динамика конструкций космических летательных аппаратов». — М.: Изд-во МГТУ им. Н.Э. Баумана, 2009. — 16 с. Рассмотрены вопросы расчета собственных колебаний конструкции космического летательного аппарата с применением метода синтеза форм колебаний в виде метода жестких границ по Крейгу — Бэмптону. Для студентов, обучающихся по специальностям «Ракетостроение» и «Космические аппараты и разгонные блоки». УДК 629.7.01(075.8) ББК 39.56 Учебное издание Дмитриев Сергей Николаевич Калугин Иван Юрьевич Тушев Олег Николаевич ПРИМЕНЕНИЕ МЕТОДА СИНТЕЗА ФОРМ ДЛЯ РАСЧЕТА КОЛЕБАНИЙ КОСМИЧЕСКОГО ЛЕТАТЕЛЬНОГО АППАРАТА Редактор С.А. Серебрякова Корректор Г.С. Беляева Компьютерная верстка С.А. Серебряковой Подписано в печать 28.01.2009. Формат 60×84/16. Усл. печ. л. 0,93. Изд. № 196. Тираж 100 экз. Заказ . Издательство МГТУ им. Н.Э. Баумана Типография МГТУ им. Н.Э. Баумана 105005, Москва, 2-я Бауманская ул., 5 © МГТУ им. Н.Э. Баумана, 2009 Д 532
ВВЕДЕНИЕ При построении расчетных конечно-элементных моделей сложных динамических систем широкое применение получил метод синтеза форм колебаний. Он используется как для уменьшения объема вычислений с помощью конечно-элементных моделей большой размерности, так и для упрощения процесса обмена данными между смежными организациями. Методом синтеза форм колебаний называют метод расчета колебаний составной системы, основанный на использовании частот и форм колебаний отдельных ее частей — подконструкций. В зависимости от условий закрепления, при которых проводится расчет подконструкций, применяют различные варианты метода синтеза форм, в частности методы свободных [1, 2] и жестких границ [3]. Наибольшее распространение в настоящее время получил вариант метода жестких границ, предложенный Крейгом и Бэмптоном [3], который в терминах обмена данными часто называют форматом Крейга — Бэмптона. Уменьшения размерности в методе жестких границ достигают благодаря замене переменных. Вместо вектора узловых перемещений во всех узловых точках подконструкции вводят в рассмотрение вектор, состоящий из перемещений в избранных узловых точках, и вектор координат форм. Для узловых перемещений и координат форм используются также термины «физические координаты» и «модальные координаты». В число избранных узловых точек следует обязательно включать те, в которых данная подконструкция соединяется с другими подконструкциями. Эти узловые точки и перемещения в них называются граничными. Поскольку эти узловые перемещения учитываются при построении сокращенной динамической модели, их называют также со
храняемыми степенями свободы, а оставшиеся — исключаемыми, или внутренними. Координаты форм в векторе неизвестных сокращенной динамической модели при использовании формата Крейга — Бэмптона соответствуют формам колебаний подконструкции при закреплении всех ее сохраняемых или граничных степеней свободы, поэтому метод является разновидностью метода жестких границ. Часть высших тонов колебаний не учитывается, поэтому размерность динамической модели снижается. В работе [4] проблема построения сокращенных динамических моделей была связана с аппроксимацией амплитудно-частотных характеристик системы в некотором частотном диапазоне. Таким образом, методам конденсации было дано некоторое теоретическое обоснование применительно к решению задач о собственных колебаниях систем, состоящих из нескольких подконструкций. Метод Крейга — Бэмп- тона встроен в ряд пакетов программ, реализующих метод конечных элементов, в том числе широко используемый на практике NASTRAN. 1. СОБСТВЕННЫЕ ЧАСТОТЫ И ФОРМЫ КОЛЕБАНИЙ ОТДЕЛЬНОЙ ПОДКОНСТРУКЦИИ И ИХ СВОЙСТВА При использовании метода конечных элементов в форме метода жесткостей [5, 6] задача расчета собственных частот и форм колебаний сводится к решению обобщенной проблемы собственных значений вида 2[ ]{ } [ ]{ } 0, M K −ω ϕ + ϕ = (1) где [K] и [M] — матрицы жесткости и массы подконструкции соответственно; 2 ω = λ — собственное число; { } ϕ — собственный вектор, ω — круговая (циклическая) собственная частота. Собственные числа и собственные векторы, являющиеся решением (1), обладают свойствами, перечисленными ниже в пп. 1–5. 1. Собственные числа являются корнями характеристического уравнения 2 det( [ ] [ ]) 0. M K −ω + = (2)
Действительно, для определения вектора собственных форм необходимо решить однородную систему уравнений 2 ( [ ] [ ]){ } 0, M K −ω + ϕ = (3) которая имеет нетривиальное решение только при условии (2). Раскрыв определитель (3), получим алгебраическое уравнение относительно 2 , n ω где n — порядок матриц [K] и [M], равный числу степеней свободы конечно-элементной модели. Уравнение (3) имеет n корней. Таким образом, конечноэлементная модель подконструкции позволяет определить n частот и n форм колебаний подконструкции. Поскольку модель является некоторой приближенной моделью для описания динамического поведения подконструкции, только часть найденных частот и форм (в нижней части спектра) близки к истинным. Собственные числа можно собрать в некоторой диагональной матрице собственных чисел 2 2 2 1 2 [ ] diag( , ,..., ), n Λ = ω ω ω (4) а собственные векторы — в матрице форм, располагая их по столбцам: 1 2 [ ] { , ,..., }. n Φ = ϕ ϕ ϕ (5) 2. Собственная форма колебаний {ϕ} определяется с точностью до постоянного множителя с. Если {ϕ} — собственная форма колебаний, то { } c ϕ — тоже собственная форма колебаний. В справедливости данного свойства можно убедиться, подставив в уравнение (1) форму { } c ϕ в качестве собственного вектора. При делении правой и левой частей на константу получим верное равенство. Для определенности собственные векторы нормируют. В дальнейшем предполагается, что собственные векторы нормированы с матрицей массы, т. е. выполняется соотношение [ ]{ } 1 M ϕ ϕ = ⎢ ⎥ ⎣ ⎦ (6) (здесь символы ⎢ ⎥ ⎣ ⎦ применяются для обозначения матрицыстроки).
3. Собственные формы колебаний, соответствующие некратным собственным значениям, ортогональны с матрицами массы и жесткости, т. е. имеют место равенства [ ][ ]{ } 0; i j M ϕ ϕ = (7) [ ][ ]{ } 0, j i M ϕ ϕ = (8) где нижние индексы i и j обозначают номера тонов колебаний. Для доказательства выпишем обобщенную проблему собственных значений для i-го и j-го тонов колебаний: 2[ ]{ } [ ]{ } 0; i i i M K −ω ϕ + ϕ = (9) 2[ ]{ } [ ]{ } 0. j j j M K −ω ϕ + ϕ = (10) Умножим уравнение (9) слева на строку j ⎢ ⎥ ϕ ⎣ ⎦ , а уравнение (10) — слева на строку iϕ ⎢ ⎥ ⎣ ⎦ , затем первое из полученных соотношений транспонируем: 2 [ ]{ } [ ]{ } 0; i i j i j M K −ω ϕ ϕ + ϕ ϕ = ⎢ ⎥ ⎢ ⎥ ⎣ ⎦ ⎣ ⎦ (11) 2 [ ]{ } [ ]{ } 0. j i j i j M K −ω ϕ ϕ + ϕ ϕ = ⎢ ⎥ ⎢ ⎥ ⎣ ⎦ ⎣ ⎦ (12) Вычтем из уравнения (11) уравнение (12) и получим 2 2 ( ) [ ]{ } 0. i j i j M ω − ω ϕ ϕ = ⎢ ⎥ ⎣ ⎦ (13) Если j i ω ≠ ω , то [ ]{ } 0, i j M ϕ ϕ = ⎢ ⎥ ⎣ ⎦ что и требовалось доказать. При выполнении соотношения (13) из (11) или (12) следует, что справедливо и соотношение ортогональности с матрицей жесткости: [ ]{ } 0. j i K ⎢ ⎥ ϕ ϕ = ⎣ ⎦ (14) 4. Собственные векторы, соответствующие кратным собственным значениям, являются линейно независимыми, но необязательно ортогональными с матрицей массы. Имеет место следующее свойство: если { } iϕ и { } j ϕ — собственные векторы, соответст
вующие кратному собственному значению λ, то любая линейная комбинация этих векторов { } { } { } i i j j C C ϕ = ϕ + ϕ (15) также будет собственным вектором для того же собственного значения λ; Ci и Cj — некоторые произвольные постоянные. Проверить данное положение можно путем подстановки (15) в обобщенную проблему собственных значений (1). Собственные векторы в этом случае определяются неоднозначно. Однако подпространство, образованное ими, является единственным. Кратные собственные значения возникают, например, при наличии симметрии в механической системе. Другим типичным случаем являются нулевые частоты незакрепленной системы, число которых в общем случае свободной пространственной системы равно шести. Добиться ортогональности собственных векторов с матрицей массы можно, проведя их принудительную ортогонализацию. В дальнейшем будем считать, что все собственные векторы (включая собственные векторы при кратных собственных числах) удовлетворяют соотношению ортогональности с матрицей массы. 5. Собственные числа и собственные векторы связаны между собой соотношением Рэлея 2 [ ]{ } , [ ]{ } K M ϕ ϕ ⎢ ⎥ ⎣ ⎦ ω = ϕ ϕ ⎢ ⎥ ⎣ ⎦ (16) которое можно получить из (1), умножив его слева на ϕ ⎢ ⎥ ⎣ ⎦ . Если собственные векторы нормированы с матрицей массы (6), то из соотношения Рэлея получим 2 [ ]{ }. i i i K ω = ϕ ϕ ⎢ ⎥ ⎣ ⎦ (17) Свойства собственных чисел и векторов, рассмотренные выше (ортогональность с матрицами масс и жесткостей, условие нормирования и соотношение Рэлея), можно объединить в два матричных уравнения: т [ ] [ ][ ] [ ]; M E Φ Φ = (18) т [ ] [ ][ ] [ ], K Φ Φ = Λ (19)
которые называются соотношениями M- и K-ортогональности; [E] — единичная матрица. 2. ОСНОВНЫЕ РАСЧЕТНЫЕ ЗАВИСИМОСТИ МЕТОДА ДИНАМИЧЕСКОЙ КОНДЕНСАЦИИ КРЕЙГА — БЭМПТОНА Рассмотрим построение сокращенной динамической модели подконструкции по Крейгу — Бэмптону применительно к задаче определения собственных частот и форм колебаний составной системы. При построении сокращенной динамической модели подконструкции ее узловые точки разделяются на граничные, в которых подконструкция соединяется с соседними, и внутренние, принадлежащие только данной подконструкции. Параметры, относящиеся к граничным узлам, будем обозначать нижним индексом 1, а к внутренним — нижним индексом 2. Будучи частью конструкции, подконструкция не является изолированной. При колебаниях конструкции в целом по одной из ее собственных форм колебаний на подконструкцию со стороны других подконструкций действуют силы, изменяющиеся по гармоническому закону и приложенные в ее граничных узлах. Уравнение колебаний отдельной подконструкции является уравнением вынужденных стационарных колебаний с круговой частотой p. Обозначим амплитудные значения граничных перемещений как {U1}, а внутренних — {U2}, вектор амплитудных значений граничных сил — {P1}. Во внутренних узлах подконструкция не нагружена, следовательно, {P2} = 0. Матрицы жесткости [K] и массы [M] разобьем на блоки, соответствующие граничным и внутренним степеням свободы (предполагается, что граничные перемещения собраны в верхней части вектора перемещений). Уравнение вынужденных стационарных колебаний подконструкции с учетом введенных выше обозначений можно записать в виде 11 12 11 12 1 1 2 21 22 21 22 2 . 0 M M K K U P p M M K K U ⎡ ⎤ ⎡ ⎤⎧ ⎫ ⎧ ⎫ − + = ⎨ ⎬ ⎨ ⎬ ⎢ ⎥ ⎢ ⎥ ⎩ ⎭ ⎣ ⎦ ⎣ ⎦⎩ ⎭ (20)
Для получения сокращенной динамической модели в уравнении (20) в соответствии с методом Крейга — Бэмптона выполним частичную замену переменных. Преобразование осуществляется введением некоторой матрицы [T], связывающей исходный вектор узловых перемещений подконструкции т 1 2 { } [ , ] U U U = с вектором, составленным из вектора узловых перемещений 1 { } U и вектора координат форм {q}: 1 1 2 [ ] . U U T U q ⎧ ⎫ ⎧ ⎫ = ⎨ ⎬ ⎨ ⎬ ⎩ ⎭ ⎩ ⎭ (21) Таким образом, вектор 1 { } U сохраняется в уравнении, а вектор 2 { } U выражается через 1 { } U и {q}. Если размерности векторов в правой и левой частях уравнения (21) одинаковы, то преобразование является точным. У сокращенной модели подконструкции размерность вектора {q} меньше размерности вектора 2 { }. U Для получения матрицы [T] должны быть предварительно определены собственные частоты и формы колебаний подконструкции при закрепленных узловых степенях свободы, соответствующих вектору узловых перемещений 1 { }, U и так называемые статические формы, связывающие векторы 1 { } U и 2 { } U при статическом нагружении подконструкции в граничных узловых точках. Исследования показали, что сокращенная динамическая модель должна давать точное решение статической задачи [2 – 4]. Только в этом случае она обладает приемлемой точностью. Поэтому включение статических форм в матрицу преобразования [T] является необходимым. Обобщенная проблема собственных значений (1) для подконструкции, закрепленной в граничных узлах, имеет вид 22 22 22 22 22 [ ][ ][ ] [ ][ ] 0. M K − Λ Φ + Φ = (22) Матрицы [M22] и [K22] имеют размер , r r × где r — число внутренних степеней свободы. Диагональная матрица собственных чисел [Λ22] содержит собственные числа, а матрица [Ф22] — собственные векторы подконструкции, закрепленной в граничных
узлах. Собственные числа и векторы удовлетворяют соотношениям M- и K-ортогональности. Число найденных собственных чисел и векторов , m r ≤ поэтому матрица [Ф22] в общем случае не квадратная. Статические формы для записи матрицы [T] представляют собой формы упругих смещений внутренних узловых точек подконструкции при последовательном задании единичных смещений по всем граничным степеням свободы при условии, что перемещения во всех оставшихся граничных точках равны нулю. Обозначим число граничных степеней свободы . k n r = − Тогда матрица статических форм должна иметь размер n k × и может быть представлена в виде E S ⎡ ⎤ ⎢ ⎥ ⎣ ⎦ , (23) где [E] — -единичная k k × матрица, а [S] подлежит определению из решения задачи статики: 11 12 1 1 21 22 2 . 0 K K U P K K U ⎡ ⎤ ⎧ ⎫ ⎧ ⎫ = ⎨ ⎬ ⎨ ⎬ ⎢ ⎥ ⎩ ⎭ ⎣ ⎦ ⎩ ⎭ (24) Совокупность перемещений 1 { } U образует матрицу [E], совокупность перемещений 2 { } U — матрицу [S]. В правой части (24) набор векторов 1 { } P образует матрицу [KG] размером , k k × коэффициенты которой равны реакциям в закрепленных узлах подконструкции. Матрица [KG] также подлежит определению. Таким образом, приходим к системе уравнений 11 12 21 22 . 0 G K K K E K K S ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ = ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ (25) Решив (25) относительно [S] и [KG], получим 1 22 21 [ ] [ ] [ ]; S K K − = − (26) 1 11 12 22 21 [ ] [ ] [ ][ ] [ ]. G K K K K K − = − (27)