We study the features of propagation of a longitudinal wave in an acoustic (mechanical) metamaterial, modeled as a one-dimensional chain, containing equal masses, connected by elastic elements (springs), and having the same rigidity. Each mass contains within itself a series connection of another mass and viscous element (damper). The mass-to-mass model is free from the drawbacks of a number of other mechanical models of metamaterials: i.e. it eliminates the need to have the property of a deformable body to possess a negative mass, density, and (or) a negative elastic modulus. It is shown that the model under consideration makes it possible to describe the dispersion and frequency-dependent attenuation of a longitudinal wave, the character of which essentially depends on the ratio of the external and internal mass of the metamaterial. The behavior of the phase and group velocities of the wave is studied, as well as the evolution of its profile, both in the low-frequency and high-frequency ranges. The mass ratios were found at which the phase velocity exceeds the group velocity (normal dispersion) in magnitude and those at which the group velocity exceeds the phase velocity (anomalous dispersion) in a wide frequency range. Having the same asymptotic values when the frequency tends to infinity, the phase and group velocities have significant differences in behavior, namely, that the phase velocity is a monotonic function of frequency, and the group velocity has a maximum. In addition, in the region of normal dispersion, the group velocity may be negative, i.e. the so-called “reverse wave” effect is true, when, despite the fact that the phase velocity is directed in the positive direction of the spatial axis, the energy in such a wave is transferred in the negative direction.

The weakening of the material begins reaching a critical level of stress state, is characterized by a decrease in the level of stress during growing deformations and can develop with an equilibrium accumulation of structural damage. The equilibrium accumulation of damage is possible if the given displacements of the boundary points are provided (that is, with “hard” loading) and if the rigidity of the loading system is sufficient. The design becomes unable to withstand the load only when zones with weakened connections are developed enough. Therefore, taking into account the full deformation diagram in the calculations allows to more accurately determine the load bearing capacity of the design. This paper gives an analytical solution for the problem of a homogeneous cylindrical solid torsion with a circular cross section with its hard loading taking into account the material weakening. Piecewise linear approximations of elastic and elastoplastic medium with a linear weakening at the supercritical deformation stage are considered. The diagrams are plotted regarding stress distribution over the cross section are given; the graphs of the maximum torque value and the extreme value of the relative angle of rotation on the parameters of the deformation diagram. The dependences of the torque on the relative angle of rotation of the sections for the stage of initial supercritical deformation, as well as the stage of supercritical deformation and fracture are determined. The graphs of the dependence of torque on the angle of rotation of the section are given. Reserves of the load bearing capacity of the design are identified. It is noted that taking into account the weakening of the material is expedient in strength calculations and in determination of the system’s safety factor.

The loading of a crack-like defect in mode II is considered. In contrast to the classical representation of a crack in the form of a mathematical cut, the proposed model defines a crack in the form of a physical cut with a characteristic linear size. The mental continuation of a physical cut in a solid forms an interaction layer. It is significant that the stress-strain state of the layer does not introduce a singularity to the crack model. The product of the increment of the specific free energy in the face square element of the layer by the linear size determines its energy product. The object of the study is a double-cantilever sample, and the subject of study is the energy product in the face element of the interaction layer. The external load of the cantilevers leads to their horizontal antisymmetric displacements, which form uniform shear deformations in the interaction layer. From the equilibrium conditions of the cantilevers in the variation form, taking into account the hypothesis of axial deformation homogeneity and their reduction, a system of differential equations is obtained, which relates the stress state in the layer and the cantilevers. The solution of the characteristic equation of the system is investigated for various ratios of layer thickness and cantilevers. It is shown that when the relationship is less than a certain value, depending on the Poisson's ratio, real roots take place. In the framework of the real roots of the characteristic equation, an analytical solution of the problem is obtained. Subject to the neglect of compression cantilevers found a simplified solution. The deformations in the layer are determined taking into account the compression of the consoles and without it. The analysis of the dependence of the energy product on the relationship of the thickness of the layer and cantilevers. It is shown that with a thickness ratio of 10-6 or less, the energy product practically does not change its value. Accounting for the compression of cantilevers gives a difference in the values of the energy product of the order of 20 % in relation to the simplified solution of the problem.

Currently, modern CAE programs like ANSYS, LS-DYNA, ABAQUS, Simcenter 3D etc. are widely used to develop each engineering product. CAE reduces development costs and speeds up product launches. Computer modelling allows eliminating unsuccessful design options at early design stages and minimising or eliminating critical design changes at the prototype testing stage. Computer simulations are especially needed for the development of high-performance perfection designs like jet engines. The use of composites in the engines requires the availability of appropriate design tools to assess the mechanical behaviour of structural elements made of such materials under operational loads and in an emergency. Predictive modelling of the deformation and fracture of the composite fan blades and fan case subjected to dynamic loading cause considerable interest of engineers. As in the case of simulation of shock loading of metals, a model of a composite material must be verified. However, the problem is that the fibrous composite material itself is already a structure and can be modelled in various ways: without considering the layered structure (homogeneous approach) taking the mesostructure into account (at the ply-level or at the yarn-level) taking the microstructure into account (at the filament level). The requirements for the initial data, the number of parameters determined during the verification process, the complexity of creating geometric models will differ in each case. This paper briefly describes the main approaches to the numerical simulation of composite elements under high-velocity impact loading. The main advantages and disadvantages of these approaches are also considered. On the example of the meso-scale approach, the main parameters of the computational models that affect the results of calculations are shown. Based on the obtained data, the main recommendations were formulated on the validation of meso-scale models of composites.

Depletion of traditional hydrocarbon reserves leads to the development of extracting methods for heavy crude oil and bitumen characterized by extremely high viscosity. The most effective technology is the steam-assisted gravity drainage. The aim of this method is to decrease oil viscosity by injection of hot steam into the reservoir. Increase of temperature, pore pressure and change of stress-strain state during this process significantly affect porosity which is the key storage parameter of the reservoir. This work is devoted to the analysis of models for porosity evolution during the steam-assisted gravity drainage process. The authors have developed an original model to describe steam-assisted gravity drainage which includes the mass balance equation for a three-phase flow, the energy balance equation involving latent heat due to vaporization/condensation of water/steam and Darcy’s law for fluid filtration. Numerical implementation of the proposed equations was based on the pressure-saturation algorithm. The results have shown a substantial qualitative and quantitative disagreement between the considered models. Coupling of porosity with volumetric strain leads to the rise of its magnitude. Models relating porosity to pore pressure show simultaneous existence of high-porous (near the injection well) and low-porous (near the production well) areas. In case when porosity is dependent on effective stress a circular area of a compacted soil is formed. Therefore, to obtain a correct estimation of the oil production rate in an arbitrary reservoir it is necessary to define the prevailing mechanism of porosity evolution (volumetric strain, pore pressure or effective stress).

Since 90’s carbonic nanotubes are broadly used in nanophysics, nanobiology and nanomechanics in nanofluidic devices, nanocontainers for gas storage and nanopipes conveying fluid. They have a perfect hollow cylindrical geometry and superior mechanical strength. The flowing fluid can be water, oil, dynamic flow of methane, ethane and ethylene molecules. The problem of the fluid-structure interaction could be considered in the case of nanoscale. However, the experiments at the nanoscale are difficult and expensive. That is why the continuum elastic models have been used to study the fluid-structure interaction. The carbon nanotubes are considered with Euler- and Timoshenko-beam models. In this paper the dynamic stability of a single-walled carbon nanotube is investigated on the basis of the Euler-beam model and with the employment of the Generalized Differential Quadrature Method. The tube under investigation is assumed hinged at its both ends and is embedded in a polymer matrix. To study the influence of the surrounding elastic medium (for example, a polymer) on the stability of the pipe, an elastic base of Pasternak is introduced. A differential equation is presented that describes the transverse vibrations of a nanotube embedded in a polymer matrix. Dimensionless parameters are introduced. The scheme of Chebycheva-Gauss-Lobato is used for sampling. The coefficients are calculated using Lagrange interpolation functions. A system of homogeneous equations is written in the matrix form. The obtained numerical results are for flowing fluids with different densities. In order to study the effect of the surrounding elastic medium (such as polymer) on the stability of the pipe the Pasternak elastic foundation is introduced. The critical velocities of each type of fluid are determined for different stiffnesses of this matrix. A decrease in the critical speed with the increasing mass ratio has been established.

A number of problems in the theory of elasticity, the theory of heterogeneous media, the theory of thin shells and plates is reduced to solving boundary value problems for systems of inhomogeneous polyharmonic equations. The paper proposes a numerical algorithm for solving systems of polyharmonic equations of the form in single-connected and multi-connected areas with a piecewise smooth contour with specified boundary conditions. Two cases are considered when the function is a known polyharmonic function and when the function is also the desired polyharmonic function. The boundary conditions can have the form similar to Dirichlet conditions, Neiman conditions, and can have a mixed form when on one part of the boundary conditions of the Dirichlet type are given, and on the other - hand, the conditions of the Neiman type. On the basis of multiple applications of the Laplace operator and the boundary element method, which is based on the green integral identity, the given system is reduced to a system of integral identities. After approximating the boundary by an inscribed n-gon and discretizing the system of integral identities, the latter is reduced to a system of linear algebraic equations, which is conveniently represented as a system of matrix equations. The existence and uniqueness of the solution follows from the existence of a unique solution of a system of linear algebraic equations. Special attention is paid to the application of the algorithm to the solution of problems on the bending of thin plates, and the bending load can be a known function, and can be an unknown polyharmonic function of an arbitrary order with given boundary conditions. The problem of bending a thin plate of elliptic shape with a known load on the surface is solved, as well as the problem of bending a thin square plate with an unknown load, which is the solution of a harmonic equation with given boundary conditions. The level lines are constructed and the forms of curved plates are given.

The effect of a discrete aggregate on the bending stability of thin-walled cylindrical shells was studied. The samples for the study were made of aluminum alloy 3004 (temper H19). The samples were cantilevered; a vertical concentrated cross force was applied to its free end. Ten empty samples and ten samples 90% filled with iron powder were tested. The samples were loaded in stages by 10N, and when approaching the moment of buckling the samples were loaded by 1N or less. The force, deflection of the free end of the sample, and axial strain were recorded at each stage of loading. Displacement of the free end of the sample versus cross force diagrams is plotted. The buckling force was determined by the inflection point of the diagrams. The samples lost the stability in elasticity. The diagrams for the empty and filled samples for each series almost coincide until the loss of stability of the empty sample. This means that when the empty sample loses its stability, the buckle directed towards the center line is formed; and for the filled samples, the formation of buckles is prevented by the discrete material, therefore the critical force increased by 18.8%. The influence of discrete aggregate is considered on the basis of approximation [24] for tank trucks manufactured at the enterprise. The critical stress is calculated using the superposition principle, since the stability is lost in the elasticity. The critical stress for the filled tank is determined by the sum of the critical stress for an empty shell with the stresses created by the weight load and the hydrostatic pressure of the discrete filler. The calculation of critical stresses showed that for the pattern samples the influence of discrete filler is 8.3%, and for the full-sized tanks the influence of discrete filler on the value of the critical voltage is significant and amounts to 62%. We studied the effect of various discrete fillers on stability, such as river sand, iron and copper powders at different filling ratios of the samples. Due to an increase in powder density and filling volume, the value of the critical force increased. For the samples filled to 90% with the river sand had the critical force increase by 13.3%, while those filled with the iron powder had their critical force increase by 40.5%, the copper powder samples had an increase of the critical force by 43.1%.

A new method is developed for solving boundary value problems of elastoplastic deformation of polycrystalline materials based on the field-theoretical approach. The boundary value problem for inhomogeneous global strain fields in differential form is transformed into a system of integral equations for mesostrain tensors in grains. In this approach, strain at any point of any grain represented as a superposition of homogeneous macrostrain and contributions of interactions with strain in given grain and all another grains of polycrystalline body . It is shown that the effects of the interaction of strains in polycrystal grains can be described using tensors of the fourth rank. This tensor has 36 independent components. The interactions are additive in nature, that drastically simplifies the solution of some problems, for example, search for extreme microstructures of a polycrystal in which critical localized phenomena arise, such as nucleation of the first plastic slips occur. Constitutive equations of deformation type are used for whole body and separate grains. A model of elastoplastic deformation of single crystals of grains is constructed. The physical mechanisms of plastic deformation are shifts in slip systems of crystals. General expressions are obtained for calculating the secant modules of single crystals for any multiaxial deformation. To solve systems of integral equations for mesostrains in grains the perturbation theory upon intergrains interaction used. Nonlinear systems of equations for plastic strains are solved by the iteration method. The features of the elastoplastic interaction of grains are theoretically investigated. The intensity of the elastoplastic interaction depends on the deformed state of the grains. For two identical grains, the elastoplastic interaction of the pair is several times more intense than the elastic one. In this case, the effect of plastically deformed grain on elastically deformed grain is much higher than the inverse effect. An increase in the intensity of interaction with the development of plastic deformations leads to the effect of homogenization of mesodeformations. Computational experiments were performed using polycrystalline titanium as an example.

18 декабря 2019 года ушел из жизни прекрасный человек и замечательный ученый – профессор Рудольф Алексеевич Васин (род. 11 августа 1937 года). Вся его сознательная жизнь связана с Московским государственным университетом, выпускником которого он стал в 1959 году. После завершения аспирантуры с 1962 года Рудольф Алексеевич пришел в НИИ механики при МГУ, где и работал до своего последнего дня, пройдя путь от младшего да главного научного сотрудника и заведующего лабораторией, от кандидата до доктора физико-математических наук и профессора. Очень быстро Рудольф Алексеевич превратился в ведущего специалиста страны в области механики деформируемого твердого тела, особенно – в экспериментальной механике, где он обладал непререкаемым авторитетом. В связи с этим огромное время ему приходилось тратить на деятельность, отвлекающую от основной научной работы, – на оппонирование огромного числа докторских и кандидатских диссертаций, рецензирование статей в ведущих журналах (вначале – Советского Союза, позднее – России). И к этой работе Рудольф Алексеевич относился чрезвычайно ответственно, можно сказать – истово. Он всегда считал своим долгом «поднимать» молодых исследователей, особенно не имеющих административной поддержки в столичных городах. С большим удовольствием Рудольф Алексеевич приезжал в Пермь, где, по его мнению, сложилась серьезная школа механиков. Несмотря на свою занятость, он без малейших колебаний вошел в состав редколлегии вначале сборника ППИ, а позднее – журнала «Вестник ПНИПУ. Механика», где активно работал более 20 лет.

Особенно следует отметить человеческие качества Рудольфа Алексеевича. В общении его всегда отличала какая-то внутренняя глубокая интеллигентность, деликатность. Будучи человеком чрезвычайно принципиальным (особенно – в вопросах науки и этики) и страстным, он умел так выстроить любой разговор, что даже его идейные противники относились к нему с глубоким уважением.

Память о Рудольфе Алексеевиче навсегда останется с нами. Он всегда будет для нас примером Человека
и Ученого.

 

Редколлегия журнала «Вестник ПНИПУ. Механика»