Другие журналы
|
научное издание МГТУ им. Н.Э. БауманаНАУКА и ОБРАЗОВАНИЕИздатель ФГБОУ ВПО "МГТУ им. Н.Э. Баумана". Эл № ФС 77 - 48211. ISSN 1994-0408![]()
77-30569/348787 Анализ сверхзвукового флаттера осесимметричной оболочки с меридианом в виде сплайна в среде SolidWorks на базе API
# 04, апрель 2012
Файл статьи:
![]() УДК 004.4 Мордовский государственный университет им. Н. П. Огарева Разрушение авиационных и ракетных конструкций может происходить не только из-за развития в них чрезмерных напряжений или из-за потери статической устойчивости, но и вследствие динамической неустойчивости в сверхзвуковом потоке газа, т.е. флаттера. Явления флаттера, связанные с элементами обшивки летательных аппаратов, принято относить к категории панельного флаттера (или флаттера пластин и оболочек) [1]. Задача о флаттере включает в себя несколько аспектов (рис. 1): · аэродинамический, связанный с определением параметров потока на поверхности обтекаемого тела; · упругий, связанный с анализом напряженно-деформированного состояния конструкции под действием внешних нагрузок; · инерционный, связанный с необходимостью учета сил инерции. Рис. 1. Основные аспекты решения задачи о флаттере
В теории панельного флаттера известны два вида флаттера: связанный и одномодовый. Первый возникает при слиянии частот колебаний оболочки в комплексной плоскости при больших числах Маха. Сравнение границ флаттера связанного типа с экспериментами при M>1.7 удовлетворительное [2]. Одномодовый флаттер возникает при 1<M<2 и изучен в настоящее время слабо, поскольку требует для своего исследования привлечения уточненной аэродинамической теории, а это сопряжено со значительными математическими трудностями [1]. Цель данной работы состоит: · в разработке методики, алгоритма и программной реализации для решения задачи о связанном флаттере для осесимметричной оболочки с меридианом в виде Безье-сплайна, при этом толщина оболочки вдоль меридиана также предполагается заданной Безье-функцией; · в создании программной среды для дальнейшего исследования оболочек рассматриваемого класса в отношении флаттера одномодового типа. Программное обеспечение, разрабатываемое нами в этих целях, назовем условно FlutterWorks. В такой постановке для решения задачи лучше всего подходит метод конечных элементов (МКЭ), а задача в целом дополняется еще двумя аспектами: · геометрическим, предполагающим создание геометрической модели оболочки, удобной для ее последующей дискретизации; · алгебраическим, связанным с решением несимметричной квадратичной обобщенной задачи на собственные значения высокого порядка и с матрицей блочно-ленточной структуры. Одной из особенностей современного этапа развития САПР в целом и, в частности, CAE-систем [3], является открытость их архитектуры, что дает возможность, используя средства API, разрабатывать специализированные библиотеки, расширяющие и дополняющие базовый функционал САПР. В качестве такой системы мы выбрали SolidWorks (SolidWorksFlowWorksSimulation). В этом смысле данная работа является продолжением работы [4]. Тогда рассматриваемая задача реализуется в виде следующей последовательности этапов: 1. Формирование геометрической модели оболочки, т.е. геометрический аспект, включающий в себя процедуру считывания параметров геометрии модели SolidWorks из дерева конструирования, в частности, координат ключевых точек сплайнов (рис. 2 а). Эта функция реализуется специальным программным модулем, разработанным на базе APISolidWorks. 2. Решение аэродинамической задачи средствами FlowWorksSimulation с последующим импортом результатов решения в FlutterWorks (аэродинамический аспект). Эту функцию также реализует специальный программный модуль, разработанный на базе API FlowWorksSimulation (рис. 2 б). 3. Формирование конечноэлементной модели для решения задачи о флаттере осесимметричной оболочки рассматриваемого типа. 4. Решение обобщенной квадратичной несимметричной задачи на собственные значения.
Геометрический аспект задачи и конечноэлементная дискретизация. После того как координаты
где
Густота сетки задается при этом величиной шага для параметра
где
К сожалению, базовый функционал SolidWorks не предусматривает возможности формирования оболочки, толщина которой была бы задана Безье-функцией, поэтому координаты ключевых точек для толщины задаются в виде уравнений, а визуализация оболочки в целом осуществляется в отдельном окне. Все связанные с названными обстоятельствами особенности построения конечноэлементной сетки иллюстрирует собой рис. 3. Рис. 3. Конечноэлементная дискретизация
Конечноэлементная модель. Как уже указывалось, для решения рассматриваемой задачи в качестве конечного элемента используется усеченный конус. На рис. 4 показаны компоненты перемещений и внутренних силовых факторов, имеющих место в поперечных и меридиональных сечениях оболочки. Конечный элемент построен на основе теории С.П.Тимошенко с учетом деформации поперечного сдвига в меридиональном направлении [6]. Таким образом, для неосесимметричной деформации имеем:
Здесь Рисунок 4. Конечный элемент оболочки
Тогда уравнение колебаний конечного элемента оболочки можно получить на основе принципа Гамильтона и в соответствии со стандартной процедурой МКЭ имеем [7]:
где
Здесь Объединяя вклад каждого конечного элемента на основе принципа прямой жесткости, и осуществляя преобразование перемещений к общей системе координат, получим для ансамбля элементов следующее уравнение:
где Тестовая задача 1. С целью тестирования разработанного программного обеспечения была решена задача по определению собственных частот колебаний жестко защемленной на торцах усеченной конической оболочки со следующими параметрами: меньший торцевой радиус На рис. 5 и 6 показано сравнение результатов решения этой задачи для различных гармоник и форм колебаний по меридиану с известным экспериментальным и аналитическим решением этой задачи [8] и результатами SolidWorksSimulation. Следует отметить совпадение этих решений с относительной погрешностью не более 2 %. Рис. 5. Результаты решения задачи анализа колебаний конической оболочки
Флаттер оболочки. Будем использовать для определения аэродинамической нагрузки поршневую теорию:
Здесь
В этом уравнении Aи D- матрицы ансамбля конечных элементов, соответственно аэродинамической жесткости, и аэродинамического демпфирования. A и
При этом следует иметь в виду, что параметры ( Порядок уравнения (4) высок, но матрицы имеют блочно-ленточную структуру. Тогда, следуя В.Г.Бунькову [10], есть смысл перейти к обобщенным координатам, в качестве которых использовать собственные векторы задачи (2), т.е. формы колебаний оболочки в вакууме. Учитывая ортогональность собственных форм и, нормируя их по кинетической энергии
где Дальнейшее решение сводится к анализу изменения собственных значений задачи (6) при возрастании параметра (плотности) потока. При этом используется дифференциальный прогноз, т.е. для каждого следующего собственного значения Для определения производной сформируем сопряженную систему:
где
в этих формулах
Тогда производная от собственного значения по плотности потока определяется в виде: Cоответственно, очередное собственное значение аппроксимируется как здесь
Здесь
При возрастании плотности потока собственные значения сближаются друг с другом и при некотором критическом значении становятся кратными (близкими к кратным), такая ситуация называется бифуркационной. Затем одно из собственных значений выходит на положительную полуплоскость, чему соответствуют колебания с возрастающей амплитудой, т.е. флаттер. На рис. 7 показаны годографы нескольких низших комплексных частот (рис. 7а), а также бифуркационная ситуация для частот, движущихся навстречу друг другу в более крупном масштабе (рис. 7 б). Тестовая задача 2. С целью тестирования разработанного программного обеспечения была решена задача по определению границы флаттера жестко защемленной на торцах усеченной конической оболочки со следующими параметрами: меньший торцевой радиус
Рис. 8. Возмущенный около оболочки поток. Распределение плотности потока.
Рис. 9. Критическое значение динамического давления, соответствующее наступлению флаттера
Пример. Задача 3. Рассмотрим задачу определения границы флаттера для почти конической оболочки с такими же параметрами, как и в Задаче 2, но являющейся элементом составной конструкции (рис. 10) и заданной сплайном. От конической оболочки Задачи 2 эта оболочка отличается тем, что в точках сопряжения сплайна, являющимся ее меридианом, заданы условия касательности с дугой окружности на левом торце и с прямой – на правом торце.
Рис. 10. Составная оболочечная конструкция На рис. 11 показан характер распределения плотности потока при обтекании конструкции (получено средствами FlowWorksSimulation), а граница флаттера, как показывают наши расчеты, изменилась по сравнению с Задачей 2 существенным образом: динамическое давление флаттера на бесконечности, соответствующее сочетанию первой и второй формы колебаний, составило Рис. 11. Качественный характер распределения плотности потока около обтекаемой оболочки (получено средствами FlowWorksSimulation)
Программная реализация. Структура программного обеспечения в целом приведена на рис. 12. AddInмодуль создает COM-интерфейс и через соответствующую структуру данных передает основной указатель m_iSldWorks = CComQIPtr<ISldWorks, &__uuidof(ISldWorks)>(ThisSW); на API SolidWorks Application. Этот же модуль обеспечивает создание в меню SolidWorksпункта MechanicsOptimum– Flutter, при обращении к которому производится загрузка основного dll-модуля FLUTTER.DLL. Функции API приложения FlowWorksSimulationреализованы в динамической библиотеке NIKCommonApi.dll, а доступ к ним организуется посредством созданного в нашей программной среде пространства имен: #import "NIKCommonApi.dll" namespace API = NIKCommonApiLib; Рис. 12. Структура программного обеспечения
Интерфейс пользователя (рис. 6) полностью встроен в окно TaskPaneSolidWorks. Работа выполнена в соответствии с базовым соглашением о сотрудничестве с компанией SolidWorksRussiaи партнерской программой SolidWorksCorp. Partner (Research). Автор выражает признательность сотрудникам компании SolidWorksRussiaД.Б. Собянину и И.В. Щекину за квалифицированную техническую поддержку, оказанную в процессе выполнения данной работы. Литература: 1. Алгазин С.Д., Кийко И.А. Флаттер пластин и оболочек. Ин-т проблем механики РАН. – М.: Наука, 2006. 247 с. 2. Веденеев В.В. Численное исследование сверхзвукового флаттера пластины с использованием точной аэродинамической теории //Механика жидкости и газа. 2009, № 2. С.169-178. 3. Норенков И. П. Основы автоматизированного проектирования: учеб. для вузов. — 4-е изд., перераб. и доп.— М.: Изд-во МГТУ им. Н. Э. Баумана, 2009. 430 с. 4. Чугунов М. В., Небайкина Ю. А. Программный модуль для решения задач оптимального проектирования в среде SolidWorks на базе API //Наука и образование: электронное научно-техническое издание, 2011, № 9. (<http://technomag.edu.ru/doc/206217.html>). 5. Шикин Е.В., Плис А.И. Кривые и поверхности на экране компьютера. – М.: «Диалог», МИФИ, 1996. 280 с. 6. Королев В.И. Слоистые анизотропные пластины и оболочки из армированных пластмасс. – М.:Машиностроение, 1965. 272 с. 7. Бате К., Вилсон Е. Численные методы анализа и метод конечных элементов: Пер. с англ. М.: Стройиздат, 1982. 447с. 8. Кольман Э.Р. Экспериментальное исследование собственных частот колебаний стальных усеченных конических оболочек вращения // В кн.: Расчеты на прочностью – М.: Машиностроение, 1968. – Вып. 8. 9. Вольмир А.С. Оболочки в потоке жидкости и газа. Задачи гидроупругости. – М.: Наука, 1979. 387 с. 10. Буньков В.Г. Проблема собственных значений матриц в задаче о флаттере // Труды ЦАГИ. 1978.Т. 17, № 2. С. 93-99. 11. Sunder P.J., Ramakrishnan C.V., Sengupta S. Optimum cone angles in aeroelastic flutter // Comp. and Struct. 1983. Vol.17, N 1. P. 25-29. 12. Ueda T., Kobayashi S., Kihira M. Supersonic Flutter of Truncated Conical Shells. // Trans. Jap. Soc. Aer. Sp. Sci. 1977. V.20, N 47. P.13 -30. Публикации с ключевыми словами: метод конечных элементов, САПР, SolidWorks, API, флаттер осесимметричной оболочки, Безье-сплайны, аэродинамическое давление, поршневая теория, обобщенная несимметричная задача на собственные значения, COM-интерфейс Публикации со словами: метод конечных элементов, САПР, SolidWorks, API, флаттер осесимметричной оболочки, Безье-сплайны, аэродинамическое давление, поршневая теория, обобщенная несимметричная задача на собственные значения, COM-интерфейс Смотри также: Тематические рубрики: Поделиться:
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|