Все выпуски

Резольвентный метод, базирующийся на преобразованиях Лежандра, применен для интегрирования уравнений баллистики в среде со степенным по скорости сопротивлением, коэффициент которого падает линейно с высотой. Во втором приближении по градиенту плотности и с учетом уменьшения с высотой ускорения свободного падения g(y) задача сведена к линейному дифференциальному уравнению. Его решением получены универсальные формулы для неоднородностной добавки к резольвентной функции f_n(b), а также к вертикальной и горизонтальной координатам δy(b), δx(b), b = tgθ - наклон траектории. Подробно рассмотрен случай квадратичного сопротивления.

The resolvent method based on Legendre transformation was applied to integrate ballistic equations of a heavy point mass in inhomogeneous medium with the drag force being power-law with respect to speed, at that the coefficient of the drag force decreases linearly with height y. General expressions were obtained for resolvent function a′′_bb(b) with a(b) being an intercept and b = tgθ, where я is inclination angle. In the second order by gradient c/m⁻¹ of perturbative approach, the universal formulas for δa′′_bb(b)-, δx(b)-, δy(b)-additions were derived. The case of Releigh resistance was considered particularly in frames of the method above and inhomogeneity influence on the motion was investigated. The falling of gravity g(y) is taken into consideration too.

Рассмотрена динамика вращения твердого тела (ротатора) вокруг неглавной оси Oz, проходящей через его центр масс, с учетом диссипативных моментов: сухого трения M_fr, возникающего в опорах из-за поперечных динамических реакций, и квадратичного по угловой скорости ω аэродинамического сопротивления M_R=-c|ω|ω. Показано, что уравнение динамики и вытекающие из него кинетики вращения тела качественно различны в общем и частном случаях инерционных и диссипативных параметров: осевого момента инерции J_zz, коэффициентов c и α=M_fr/√ε²+ω⁴ (ε - угловое ускорение). В частном случае равенства J_zz=c=α обнаружено отсутствие физически возможного решения для вращения по инерции в рамках динамики абсолютно твердого тела. Парадокс разрешается через нормализующее причинно-следственные связи введение запаздывающих величин ε(t-τ) и ω(t-τ), определяющих в согласии с принципом Даламбера поперечные реакции в опорах оси M_x,y(t-τ) и пару M_fr(t-τ). Последняя же определяла темп потери кинетического момента dK_z(t)/dt в момент времени t. Кинетика вращения при этом имеет импульсивный характер так называемого фрикционно-аэродинамического удара. Также путем численного интегрирования продемонстрирована необычная угловая кинетика φ(t) затухающих колебаний ротатора под действием упругого момента M_e=-κφ, характеризующаяся наличием двух фаз: кратковременного стартового участка, зависящего от начальных условий, затем резко переходящего в фазу почти синусоидальных колебаний с медленно убывающей амплитудой.

The article studies the rotational dynamics of a rigid body (rotator) around the central but non-principal axis Oz passing through its center of mass under the action of dry frictional torque M_fr=α√ε²+ω⁴ caused by inertia forces in the axis's supports and the drag momentum M_R=-c|ω|ω quadratic in angular speed ω. It has been shown that the dynamical equations and the equations of the kinetics of the body's rotation, which follow from the dynamical equations, are qualitatively different in general and particular cases of the inertial and dissipative parameters involved: the axial moment of inertia J_zz and the coefficients c and α=Mfr/√ε²+ω⁴ (where ε is the angular acceleration). It is found that in the particular case of the equality J_zz=c=α a physical feasible solution for the inertial rotation within the dynamics of a perfectly rigid body is absent. The paradox is resolved by the introduction of the lagged angular velocity ω(t-τ) and acceleration ε(t-τ) as factors defining due to D'Alembert principle the supports' transversal reactions M_x,y(t-τ) and hence the value of M_fr(t-τ). The last one determines the loss rate of kinetic momentum, i.e. the dK_z(t)/dt at time t. The rotational kinetics had a type of frictional-aerodynamic impact. Also, by numerical integration, there was shown the unusual angular kinetics φ(t) of the damping oscillations of the rotator under the action of the elastic torque M_e=-κφ. The kinetics was characterized by the presence of two phases: the short starting part strongly depending on initial conditions followed by the phase of almost sine wave oscillations with extremely slow damping.

Рассматриваются два подхода к решению задачи математического моделирования обтекания метаемых тел: численное решение уравнений движения сплошной среды Навье-Стокса, осредненных по Рейнольдсу (RANS - Reynolds-averaged Navier–Stokes), с использованием модели турбулентности и прямое численное моделирование (DNS - Direct Numerical Simulation). Тестирование рассматриваемых подходов проводится при решении задачи обтекания тел вращения с простой геометрией: сферы и цилиндра с конической головной частью, для которых известны значения коэффициентов сопротивления при различных числах Маха. Проведено качественное и количественное сравнение результатов обтекания рассматриваемых тел сверхзвуковым потоком, полученным по методикам RANS и DNS. Апробация методики численного моделирования проводится для метаемого тела (снаряда) характерной формы. Представлены результаты численного моделирования обтекания снаряда для широкого диапазона параметров: чисел Маха и углов нутации. Выполнено сравнение расчетных значений коэффициентов сопротивления с эмпирическими эталонными зависимостями по законам 1943 и 1958 годов.

Two approaches to the problem of numerical simulation of streamlined bodies airflow are considered. These approaches are: numerical calculation of the Reynolds-averaged Navier-Stokes equations (RANS) using the turbulence model and direct numerical simulation (DNS). Testing of the considered approaches were conducted by solving the problem of flow past bodies of revolution with simple geometries: sphere and cone cylinder, for which values of drag coefficient at different Mach numbers are known. Qualitative and quantitative comparison of the results for the supersonic flow (modelled by RANS and DNS methods) around the bodies under consideration are carried out. The numerical simulation method is tested by considering the missile body (projectile) of characteristic shape. The numerical simulation results for the flow around the projectile are presented for a wide range of parameters: Mach numbers and angles of nutation. The calculated values of the drag coefficients are compared to the empirical reference dependencies according to the laws of 1943 and 1958.