WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

загрузка...
   Добро пожаловать!

Pages:     || 2 |

На правах рукописи

Данг Нгок Ань

РАЗРАБОТКА МЕТОДА РАСЧЁТА ПАРАМЕТРОВ ГЕОМЕТРИЧЕСКОЙ СТАБИЛЬНОСТИ КОРПУСОВ АНТЕННЫХ РЕФЛЕКТОРОВ

Специальность 01.02.06 – Динамика, прочность машин и аппаратуры

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук

Москва – 2008

Работа выполнена в Московском государственном техническом университете им. Н.Э. Баумана

Научный руководитель: доктор технических наук, профессор

Попов Борис Глебович

Официальные оппоненты: доктор технических наук, профессор

Сарбаев Борис Сафиулович

кандидат технических наук

Никитенко Вячеслав Иванович

Ведущая организация РКК «Энергия»

Защита состоится « 16 » октября 2008 г. в 14 ч. 30 мин. на заседании диссертационного совета Д 212.141.03 при МГТУ им. Н.Э. Баумана по адресу: 107005, г. Москва, 2-ая Бауманская ул., д.5.

С диссертацией можно ознакомиться в библиотеке Московского государственного технического университета им. Н.Э. Баумана.

Автореферат разослан «16» сентября 2008 г.

Ученый секретарь

диссертационного совета Карпачев А. Ю.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

- геометрической точности;

- гладкости поверхности;

- размеростабильности при нагреве;

- высокому качеству отражающей поверхности

- устойчивости к ультрафиолетовому излучению, к атомарному кислороду и другим факторам космического пространства.

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

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

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

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

Для достижения этой цели были поставлены и решены следующие задачи:

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

2. Разработать алгоритм и программу расчёта, позволяющие дать заключение о соответствии конструкции по заданным геометрическим параметрам.

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

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

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

Практическая ценность работы. Разработанный новый метод расчёта позволяет выбрать структуры несущих композитных слоёв, обеспечивающие наименьшие отклонения от требуемой геометрической формы рефлектора. Результаты диссертационной работы использованы и внедрены на предприятии РКК «Энергия».

Апробация работы. Основные положения диссертации и полученные результаты докладывались на Международной научно-практической конференции «Участие молодых учёных, инженеров и педагогов в разработке и реализации инновационных технологий» (Москва, 20-24 ноября 2006), III Международной конференции «Ракетно- космическая техника: фундаментальные и прикладные проблемы» (г. Москва, 19 - 23 ноября 2007 г.), научном семинаре кафедры «Прикладная механика» МГТУ им Н.Э. Баумана (2008).

Публикации. По результатам исследований опубликованы 2 печатные работы и тезисы, список публикаций приведен в конце автореферата.

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

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

Для математического описания деформированного состояния применён подход Лагранжа (материальный), когда в качестве независимых переменных использованы материальные координаты частицы тела. Учет геометрической нелинейности в квадратичном приближении выполнен на основе гипотез Кирхгофа-Лява. Деформированное состояние описано компонентами тензора деформаций Грина-Лагранжа.

Рассмотрена оболочка вращения, представляющая собой поверх­ность, полученную вращением образующей вокруг вертикальной оси (рис. 1).

Pиc.1. Схема оболочки вращения

Рис. 2. Кинематическая схема перемещений точек срединной поверхности оболочки

В качестве перемещений в локальной системе координат (ЛСК) принимают меридианальные перемещения u, окружные v и нормальные перемещения w (рис. 1, 2). В качестве перемещений в глобальной системе координат (ГСК) принимают радиальные перемещения ur, осевые перемещения uz и окружные v. Перемещения u, w связаны с ur, uz следующим образом:

или (1)

Деформационные соотношения, связывающие обобщённые деформации с перемещениями, соответствуют гипотезам Кирхгофа-Лява, изложенным в работах Бидермана В. Л., Феодосьева В. И., Усюкина В. И..

В случае линейности деформаций деформационные соотношения выражаются формулами:

В случае учета нелинейности эти соотношения имеют вид:

где

- матрица, связывающая деформации и изменения кривизны с перемещениями.

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

Рис. 3. Система координат и напряжения в отдельном слое

Согласно закону Гука напряжения связаны с деформациями соотношениями упругости

, (4)

По определению погонные усилия [Н/м]

,,, (5)

погонные моменты [Н]

,,. (6)

Рис. 4. К определению

внутренних силовых

факторов

Рис. 5. К Выбору координатной поверхности

в многослойном пакете:

а) Координатная поверхность совпадает

со срединной поверхностью

б) Координатная поверхность совпадает с

нижней лицевой поверхностью

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

или в матричной форме

, (7)

где,

- вектор-столбец обобщенных деформаций,

, - вектор-столбцы начальных и температурных внутренних силовых факторов.

В этих соотношениях коэффициенты Bij называют коэффициентами мембранных жесткостей:

; (8)

коэффициенты Cij называют коэффициентами смешанных жесткостей:

; (9)

коэффициенты Dij называют коэффициентами изгибных жесткостей:

; (10)

n - число слоев в многослойном пакете; - коэффициент матрицы Eij соотношений упругости для k-го однонаправленного композиционного слоя;, - нормальные координаты нижней и верхней поверхностей k-го слоя.

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

где - длина образующей конической оболочки;

=[]Т- вектор-столбец обобщенных возможных деформаций; =[N1, N2, N12, M1, M2, M12, ]T- вектор-столбец внутренних силовых факторов; U=[u, w,v]T – вектор-столбец возможных перемещений; p=[pu, pw, pv]T- вектор-столбец распределённых поверхностных касательных и нормальных сил.

Вторая глава посвящена построению конечного элемента (КЭ) многослойной композитной оболочки вращения. Оболочки вращения набирали КЭ конических оболочек (рис. 6). Отдельный КЭ конической оболочки будет определяться нормальными круговыми сечениями 1 и 2 (рис. 7). Для текущего сечения с аргументом x радиус параллели будет равен:

, (12)

где s1 – расстояние вдоль образующей от начальной параллели до параллели первого сечения КЭ.

Рис. 6. Конечные элементы конических оболочек

Рис. 7. Координаты отсчета конечного элемента

Для равновесного состояния отдельного КЭ запишем формулировку ПВП:

, (13)

где: q=[ur1, uz1, v1, 1, ur2, uz2, v2, 1]Т- глобальные узловые возможные степени свободы; R=[ r1Nr1, r1Nz1, r1N12, r1M1, r2Nr2, r2Nz2, r2N12, r2M2]T- вектор-столбец обобщенных узловых реакций; - угол поворота нормального сечения в плоскости меридиана.

Для приближенного описания свойств КЭ задана аппроксимация полей перемещений:

или

U = Ф a, (14)

где а = [ a1, a2, a3, a4, a5, a6, a7, a8]T - коэффициенты аппроксимации.

Ф = - матрица функции формы;

U =- вектор узловых перемещений

Согласно деформационным соотношениям (2) для аппроксимаций перемещений (14) получим

= L1U = L1Ф a = B a, (15)

где B=L1Ф=

=

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

U = Фa, (16)

где a вектор-столбец (8х1) произвольных коэффициентов.

Подстановка (14), (16) в вариационное уравнение (13) с учетом соотношений упругости (7) дает

ФТp rdx - qTR = 0,

или (17)

где

. (18)

Для того, чтобы перейти от коэффициентов аппроксимации а = [ a1, a2, a3, a4, a5, a6, a7, a8]T к степеням свободы в ГСК q=[ur1, uz1, v1, 1, ur2, uz2, v2, 2]Т, cначала перейдем к степеням свободы в ЛСК q л = [ u1, w1, v1, 1, u2, w2, v2,2 ] T.

Дадим определение степеням свободы в ЛСК

Эти уравнения рассмотрим как систему уравнений для определения коэффициентов

T a = qл, (19)

где

Решив (19), выразим коэффициенты аппроксимации а через q л :

a = T -1q л.

Аналогичным образом определим связь

a = T -1q л .

Тогда формулировка равновесного состояния для КЭ (17) примет вид

или (20)

где. (21)

Теперь перейдем к степеням свободы в ГСК. Для этого выразим степени свободы в ЛСК q л через степени свободы в ГСК q :

q л = А q, (22)

или

Аналогично q л = А q.

Тогда формулировку равновесного состояния для КЭ (20) запишем через степени свободы в ГСК:

(23)

где (24)

Из уравнений (23) в силу произвольности коэффициентов q следует искомое уравнение

, (25)

Pages:     || 2 |






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»