All issues

Statement of a problem of management of movement of a body in a viscous liquid is given. Movement bodies it is induced by moving of internal material points. On a basis the numerical decision of the equations of movement of a body and the hydrodynamic equations approximating dependencies for viscous forces are received. With application approximations the problem of optimum control of body movement dares on the set trajectory with application of hybrid genetic algorithm. Possibility of the directed movement of a body under action is established back and forth motion of an internal point. Optimum control movement direction it is carried out by motion of other internal point on circular trajectory with variable speed.

The article discusses the model of the anthropomorphic type of mechanism of the exoskeleton with links of variable length. Four models of parts of variable length are considered comprehensively: the model link of the exoskeleton of variable length with a resilient member and a rigid strong core; the model of the telescopic link; the model link with the masses in the hinge-joint between them; the link model with an arbitrary number of masses. The differential equations of motion in the form of Lagrange equations of the second kind are made. On the basis of analysis of differential equations of motion for multi-link rod of a mechanical system type, exoskeleton revealed their structure, which allowed us to represent them in vector-matrix form. The General pattern of building matrices are established for the first time and the generalization of the expressions for elements of matrices in two-dimensional case are obtained. New recursive and matrix methods of composing of differential equations of motion are given. A unified approach to constructing differential equations of motion of the exoskeleton based on the developed recursive and matrix methods write differential equations of motion of the proposed exoskeleton. Comparison of the time of writing the differential equations of motion proposed methods, in comparison with the Lagrange equations of the second kind, in the system of computer mathematics Mathematica conducted. An analytical study of the model of the exoskeleton carried out. It was found that for mechanisms with n movable links of the Cauchy problem for systems of differential equations of motion for any initial conditions there is no single and unlimited continue. Control of the exoskeleton is accomplished using the torques which are located in the hinge-joints in the joints of the links and simulating control actions. Numerical investigation of a model of the exoskeleton is made, a comparison of results of calculations for exoskeletons with various models of units is held. A numerical study of the empirical evidence about the man and his movements is used. It is established that the choice structure of the exoskeleton model with lumped masses is more preferable to a model with perfectly rigid strong core. As an exoskeleton, providing comfortable movement of people, and you should repeat the properties of the musculoskeletal system.

The work is devoted to the development of a computational algorithm of the Cartesian grid method for studying the interaction of a shock wave with moving bodies with a piecewise linear boundary. The interest in such problems is connected with direct numerical simulation of two-phase media flows. The effect of the particle shape can be important in the problem of dust layer dispersion behind a passing shock wave. Experimental data on the coefficient of aerodynamic drag of non-spherical particles are practically absent.

Mathematical model is based on the two-dimensional Euler equations, which are solved in a region with varying boundaries. The defining system of equations is integrated using an explicit scheme and the Cartesian grid method. The computational algorithm at the time integration step includes: determining the step value, calculating the dynamics of the body movement (determining the force and moment acting on the body; determining the linear and angular velocities of the body; calculating the new coordinates of the body), calculating the gas parameters. At each time step, all cells are divided into two classes – external (inside the body or intersected by its boundaries) and internal (completely filled with gas). The solution of the Euler equations is constructed only in the internal ones. The main difficulty is the calculation of the numerical flux through the edges common to the internal and external cells intersected by the moving boundaries of the bodies. To calculate this flux, we use a two-wave approximation for solving the Riemann problem and the Steger-Warming scheme. A detailed description of the numerical algorithm is presented.

The efficiency of the algorithm is demonstrated on the problem of lifting a cylinder with a base in the form of a circle, ellipse and rectangle behind a passing shock wave. A circular cylinder test was considered in many papers devoted to the immersed boundary methods development. A qualitative and quantitative analysis of the trajectory of the cylinder center mass is carried out on the basis of comparison with the results of simulations presented in eight other works. For a cylinder with a base in the form of an ellipse and a rectangle, a satisfactory agreement was obtained on the dynamics of its movement and rotation in comparison with the available few literary sources. Grid convergence of the results is investigated for the rectangle. It is shown that the relative error of mass conservation law fulfillment decreases with a linear rate.

A mathematical model is formulated for the coastal slope erosion of sandy channel, which occurs under the action of a passing flood wave. The moving boundaries of the computational domain — the bottom surface and the free surface of the hydrodynamic flow — are determined from the solution of auxiliary differential equations. A change in the hydrodynamic flow section area for a given law of change in the flow rate requires a change in time of the turbulent viscosity averaged over the section. The bottom surface movement is determined from the Exner equation solution together with the equation of the bottom material avalanche movement. The Exner equation is closed by the original analytical model of traction loads movement. The model takes into account transit, gravitational and pressure mechanisms of bottom material movement and does not contain phenomenological parameters.

Based on the finite element method, a discrete analogue of the formulated problem is obtained and an algorithm for its solution is proposed. An algorithm feature is control of the free surface movement influence of the flow and the flow rate on the process of determining the flow turbulent viscosity. Numerical calculations have been carried out, demonstrating qualitative and quantitative influence of these features on the determining process of the flow turbulent viscosity and the channel bank slope erosion.

Data comparison on bank deformations obtained as a result of numerical calculations with known flume experimental data showed their agreement.

Based on the modeling method using a mesh system, an algorithm is implemented for solving a unsteady problem with moving bodies The algorithm takes into account the movement and rotation of bodies according to a given law of motion. The algorithm is applied to analysis the flow around an infinite composed of cylinders with an elliptical cross-section, which either move across the flow or rotate with a change in the angle of attack. To simulate the flow of bodies with a sharp edge, characteristic of the profiles of gas turbine machines, an algorithm for constructing a mesh of type C with the inclusion of a certain area behind the profile is implemented. The program for modeling the flow near the profile is implemented within the framework of models of Euler equations, Navier – Stokes equations in the approximation of a thin layer with laminar viscosity and turbulent viscosity in the framework of an algebraic viscosity model. The program has also been adapted to solve the problems of internal gas dynamics of turbomachines. For this purpose, the method of setting the boundary conditions at the entrance and exit from the calculated area from the velocity to the pressure drop, as well as at the lateral boundaries from the free flow to the periodicity, was changed. This made it possible to simulate the flow of gas in the inter-blade channels of compressors and turbines of gas turbine engines. To refine the algorithm, a series of calculations of the aerodynamic parameters of several turbine cascades in various subsonic and supersonic modes and their comparison with the experiment were carried out. Calculations of turbine grating parameters were carried out within the framework of the inviscid and viscous gas model. The calculation and experiment were compared by the distribution of gas parameters near the profile, as well as by the energy losses of the flow in the cascade. Calculations have shown the applicability and correctness of the program to solve this class of problems. To test the program on the problems of external subsonic aerodynamics, calculations of the aerodynamic characteristics of an isolated airfoil in an undisturbed flow were performed. The results obtained allow us to assert the applicability of the hybrid grid method to various classes of problems of applied gas dynamics.

The paper presents results that confirm the ability to control the movement of the kink of the modified sine-Gordon equation with variable external force parameters. Three types of external influences have been considered: permanent action, periodic action with a constant frequency and a frequency-modulated periodic exposure. The dependences of the position and velocity of the kink on time for various values of the parameters of external influence were obtained using the method of McLaughlin and Scott. It is shown that by changing the settings, one can adjust the velocity and direction of movement of the kink.

This article is devoted to the study of self-propulsion of bodies in a fluid by the action of internal mechanisms, without changing the external shape of the body. The paper presents an overview of theoretical papers that justify the possibility of this displacement in ideal and viscous liquids.

A special case of self-propulsion of a rigid body along the surface of a liquid is considered due to the motion of two internal masses along the circles. The paper presents a mathematical model of the motion of a solid body with moving internal masses in a three-dimensional formulation. This model takes into account the three-dimensional vibrations of the body during motion, which arise under the action of external forces-gravity force, Archimedes force and forces acting on the body, from the side of a viscous fluid.

The body is a homogeneous elliptical cylinder with a keel located along the larger diagonal. Inside the cylinder there are two material point masses moving along the circles. The centers of the circles lie on the smallest diagonal of the ellipse at an equal distance from the center of mass.

Equations of motion of the system (a body with two material points, placed in a fluid) are represented as Kirchhoff equations with the addition of external forces and moments acting on the body. The phenomenological model of viscous friction is quadratic in velocity used to describe the forces of resistance to motion in a fluid. The coefficients of resistance to movement were determined experimentally. The forces acting on the keel were determined by numerical modeling of the keel oscillations in a viscous liquid using the Navier – Stokes equations.

In this paper, an experimental verification of the proposed mathematical model was carried out. Several series of experiments on self-propulsion of a body in a liquid by means of rotation of internal masses with different speeds of rotation are presented. The dependence of the average propagation velocity, the amplitude of the transverse oscillations as a function of the rotational speed of internal masses is investigated. The obtained experimental data are compared with the results obtained within the framework of the proposed mathematical model.

In a number of fundamental and practical problems, it is necessary to describe the dynamics of the motion of complexshaped particles in a high-speed gas flow. An example is the movement of coal particles behind the front of a strong shock wave during an explosion in a coal mine. The paper is devoted to numerical simulation of the dynamics of translational and rotational motion of a square-shaped body, as an example of a particle of a more complex shape than a round one, in a supersonic flow behind a passing shock wave. The formulation of the problem approximately corresponds to the experiments of Professor V. M. Boiko and Professor S. V. Poplavski (ITAM SB RAS).

Mathematical model is based on the two-dimensional Euler equations, which are solved in a region with varying boundaries. The defining system of equations is integrated using an explicit scheme and the Cartesian grid method which was developed and verified earlier. The computational algorithm at the time integration step includes: determining the step value, calculating the dynamics of the body movement (determining the force and moment acting on the body; determining the linear and angular velocities of the body; calculating the new coordinates of the body), calculating the gas parameters. To calculate numerical fluxes through the edges of the cell intersected by the boundaries of the body, we use a two-wave approximation for solving the Riemann problem and the Steger – Warming scheme.

The movement of a square with a side of 6 mm was initiated by the passage of a shock wave with a Mach number of 3,0 propagating in a flat channel 800 mm long and 60 mm wide. The channel was filled with air at low pressure. Different initial orientation of the square relative to the channel axis was considered. It is found that the initial position of the square with its side across the flow is less stable during its movement than the initial position with a diagonal across the flow. In this case, the calculated results qualitatively correspond to experimental observations. For the intermediate initial positions of a square, a typical mode of its motion is described, consisting of oscillations close to harmonic, turning into rotation with a constant average angular velocity. During the movement of the square, there is an average monotonous decrease in the distance between the center of mass and the center of pressure to zero.

The paper presents the physical and mathematical formulation of the problems of internal ballistics of an artillery shot for a charge consisting of a set of powder tubes and their stress-strain state. Combustion and movement of a bundle of powder tubes along the barrel channel is modeled by an equivalent tubular charge of all-round combustion. It is assumed that the equivalent tube moves along the axis of the bore. The speed of movement of an equivalent tubular charge and its current position are determined from Newton’s second law. When calculating the flow parameters, two-dimensional axisymmetric equations of gas dynamics were used, for the solution of which an axisymmetric orthogonalized difference grid is constructed, which adapts to the flow conditions. The control volume method is used to numerically solve the system of gas-dynamic equations. The gas parameters at the boundaries of the control volumes are determined using a self-similar solution to the Godunov’s problem of the decay of an arbitrary discontinuity. The stress-strain state is modeled for a separate burning powder tube located in the field of gas-dynamic parameters. The calculation of the gas-dynamic parameters of the shot is carried out without taking into account the deformed state of the powder elements. The behavior of powder elements during firing is considered under these conditions. The finite element method with the division of the calculation area into triangular elements is used to solve the problem of elasticity. In the process of powder tube burnout, the computational grid on each time layer of the dynamic problem is completely updated due to a change in the boundaries of the powder element due to combustion. The paper shows the time dependences of the parameters of the internal ballistics process and the stress-strain state of powder elements, as well as the distribution of the main parameters of the flow of combustion products at different points in time. It has been established that the tubular powder elements during the shot experience significant deformations, which must be taken into account when solving the basic problem of internal ballistics. The data obtained give an idea of the level of equivalent stresses acting at various points of the powder element. The results obtained indicate the relevance of the conjugate formulation of the problem of gas dynamics and the stress-strain state for charges consisting of tubular powders, since this allows a new approach to the design of tubular charges and opens up the possibility of determining the parameters on which the physics of the combustion process of gunpowder significantly depends, therefore, and the dynamics of the shot process.

Buckling problems of thin elastic shells have become relevant again because of the discrepancies between the standards in many countries on how to estimate loads causing buckling of shallow shells and the results of the experiments on thinwalled aviation structures made of high-strength alloys. The main contradiction is as follows: the ultimate internal stresses at shell buckling (collapsing) turn out to be lower than the ones predicted by the adopted design theory used in the USA and European standards. The current regulations are based on the static theory of shallow shells that was put forward in the 1930s: within the nonlinear theory of elasticity for thin-walled structures there are stable solutions that significantly differ from the forms of equilibrium typical to small initial loads. The minimum load (the lowest critical load) when there is an alternative form of equilibrium was used as a maximum permissible one. In the 1970s it was recognized that this approach is unacceptable for complex loadings. Such cases were not practically relevant in the past while now they occur with thinner structures used under complex conditions. Therefore, the initial theory on bearing capacity assessments needs to be revised. The recent mathematical results that proved asymptotic proximity of the estimates based on two analyses (the three-dimensional dynamic theory of elasticity and the dynamic theory of shallow convex shells) could be used as a theory basis. This paper starts with the setting of the dynamic theory of shallow shells that comes down to one resolving integrodifferential equation (once the special Green function is constructed). It is shown that the obtained nonlinear equation allows for separation of variables and has numerous time-period solutions that meet the Duffing equation with “a soft spring”. This equation has been thoroughly studied; its numerical analysis enables finding an amplitude and an oscillation period depending on the properties of the Green function. If the shell is oscillated with the trial time-harmonic load, the movement of the surface points could be measured at the maximum amplitude. The study proposes an experimental set-up where resonance oscillations are generated with the trial load normal to the surface. The experimental measurements of the shell movements, the amplitude and the oscillation period make it possible to estimate the safety factor of the structure bearing capacity with non-destructive methods under operating conditions.