Анализ неоднозначности угла атаки для модели аэродинамического маятника в потоке среды
Authors: Беляков Д.В. 
Published: December 2024

Работа посвящена построению и исследованию математической модели автоколебаний аэродинамического маятника в потоке среды. В качестве модели воздействия среды на тело принята модель квазистатического обтекания пластинки средой. Согласно этой гипотезе, аэродинамические силы, действующие на тело, прикладываются в центре давления. В рассматриваемой задаче центр давления является подвижным относительно пластинки. Получены уравнения движения для рассматриваемого тела. Проведен переход к новым безразмерным переменным. Показано нарушение единственности при определении угла атаки. Проведен параметрический анализ областей неоднозначности. Показано, что в наиболее характерном положении равновесия, соответствующем состоянию покоя областей неоднозначности нет. Проведено исследование устойчивости различных положения равновесия, в которых реализован критерий Гурвица и изображены области устойчивости. Показано, что силы аэродинамического воздействия для тел с одними формами могут способствовать развитию автоколебаний, а для других затуханию. В математическом пакете MATLAB предложен комплекс программ, позволяющих проводить численные исследования, реализующий численное интегрирование уравнений описывающих колебания пластинки с неподвижным центром давления. Такая модель возможна при условии того, что длина стержня намного больше ширины пластинки. При запуске программы строится область устойчивости и на ней вводятся геометрические параметры: жесткость пружины и длина стержня. Далее вводится вектор начальных условий. При поиске численного решения используется процедура ode45, реализующая методы Рунге-Кутта четвертого и пятого порядка с переменным шагом. В файле-функции на каждом шаге метода Рунге-Кутта решается нелинейное уравнение из которого определяется угол атаки. При поиске численного решения экспериментальные аэродинамические функции интерполируются кубическим сплайном. Полученное путем интегрирования решение изображается на графике в виде фигур Лиссажу. Таким образом, разработана математическая модель колебаний пластинки, проведен параметрический анализ устойчивости, с помощью комплекса программ на базе специализированной системы компьютерной математики есть возможность подтвердить полученные аналитические результаты.
The work is devoted to the construction and study of a mathematical model of self-oscillations of an aerodynamic pendulum in a medium flow. The model of quasi-static flow of the plate by the medium was adopted as a model of the impact of the medium on the body. According to this hypothesis, the aerodynamic forces acting on the body are applied at the center of pressure. In the problem under consideration, the center of pressure is movable relative to the plate. A transition to new dimensionless variables is carried out. The violation of uniqueness in determining the angle of attack is shown. A parametric analysis of ambiguity domains is carried out. It is shown that in the most characteristic position of equilibrium, corresponding to the state of rest, there are no ambiguity regions. The stability of the equilibrium position corresponding to the state of rest was studied, in which the Hurwitz criterion was implemented, and the stability region was depicted. It is shown that aerodynamic forces for bodies with some forms can contribute to the development of self-oscillations, and for others, attenuation. In the MATLAB mathematical package, a set of programs is proposed allowing to carry out numerical research, realizing the numerical integration of equations describing the vibrations of a plate with a fixed pressure center. Such a model is possible, provided that the length of the rod is much greater than the width of the plate. When starting the program, a stability area is built and geometric parameters are entered on it: spring stiffness and rod length. The initial conditions vector is then entered. To find a numerical solution, the ode45 procedure is used, which implements the Runge-Kutt methods of the fourth and fifth orders with a variable step. In the function-file, at each step of the Runge-Kutta method, a nonlinear equation is solved from which the angle of attack is determined. When searching for a numerical solution, the experimental aerodynamic functions are interpolated by a cubic spline. The solution obtained by integration is shown on the graph in the form of Lissajous figures. Thus, a mathematical model of plate oscillations has been developed, a parametric analysis of stability has been carried out, and with the help of a set of programs based on a specialized system of computer mathematics, it is possible to confirm the obtained analytical results.