Rubric: "01.02.00 Mechanics"
The suggested thermocreep theory for thin multilayer plates is based on analysis of general three dimensional nonlinear theory of thermalcreep by constructing asymptotic expansions in terms of a small parameter being the ratio of a plate thickness and a characteristic length. Here we do not introduce any hypotheses on a distribution character for displacements and stresses through the thickness. Local problems were formulated for finding stresses in all structural elements of a plate. It was shown that the global (averaged by the certain rules) equations of the plate theory were similar to equations of the Kirchhoff–Love plate theory, but they differed by a presence of the three-order derivatives of longitudinal displacements. The method developed allows to calculate all six components of the stress tensor including transverse normal stresses and stresses of interlayer shear. For this purposes one needs to solve global equations of thermal creep theory for plates, and the rest calculations are reduced to analytical formulae use.
Dimitrienko Y., Gubareva E., Yurin Y. Asymptotic theory of thermocreep for multilayer thin plates. Маthematical Modeling and Coтputational Methods, 2014, №4 (4), pp. 18-36
On the basis of mathematical model of thermal interaction between inclusion and the matrix we estimated influence of inclusions deviations from spherical shape on the effective thermal conductivity coefficient of the composite and associated with such deviation a possible occurrence of the anisotropy of the composite with respect to the property of thermal conductivity. Using the dual variational formulation of the stationary problem of heat conduction in an inhomogeneous body we built bilateral estimates of effective thermal conductivity.
Zarubin V., Kuvyrkin G., Savelyeva I. Effective thermal conductivity of a composite in case of inclusions shape deviations from spherical ones. Маthematical Modeling and Coтputational Methods, 2014, №4 (4), pp. 3-17
A model for calculation of a rock stress-strain state considering creep is suggested. The algorithm for finite element solving the three-dimensional creep problem using finite-difference scheme of Euler's method with respect to time is presented. The specialized software is developed allowing the computer to build 3D-models of rock areas based on the initial series of 2D images, obtained with the seismic data, and to perform finite element calculation of variations in rock strain-stress state with time. Numerical simulation of rock stress-strain state was conducted on the example of a zone of the Astrakhan oil and gas field. It was found that there occurs rock mass rising in some points, and in the other points it can slope down with time. The creep rate of different layers is not the same — the highest values of the creep rate are realized in the layers of clay and sand, filled with fluid, which have the most notable creep properties. The developed algorithm and software for numerical simulation proved to be quite effective and can be applied to the study of rock stress-strain state.
Dimitrienko Y., Yurin Y. Finite element simulation of the rock stress-strain state under creep. Маthematical Modeling and Coтputational Methods, 2015, №3 (7), pp. 101-118
The article describes a class of promising anisogrid structures representing mesh shell of unidirectional carbon. A brief analysis of existing approaches to modeling deformation of grid structures is presented. New mathematical and numerical models are proposed for reliable description of complex behavior of anisogrid structures under different kinds of loads. A high degree of accuracy and stability of the numerical model based on the expansions of unknown functions in Chebyshev polynomials and Fourier series is caused by the lack of saturation of such approximations. Efficiency of the proposed models and techniques is demonstrated on the example of solving test boundary-value problems and a problem of axial compression of anisogrid cylindrical shell.
Golushko S., Semisalov B. Numerical modeling of anisogrid structures deformation using schemes of high accuracy without saturation. Маthematical Modeling and Coтputational Methods, 2015, №2 (6), pp. 23-45
In the paper we propose an algorithm of parameters (time constants of the turbine) iden-tification using the gradient method with an adaptive model. The adaptive mathematical model has the same structure as the identification object. The identification criterion is based on the loss function, which is the misalignment between the left and right sides of the equation, which describes the adaptive model. Thus it is avoided the need of finding the solution of a nonlinear equation for the adaptive model in an explicit form. In the model the signal observed at the output of the identified object is used instead of the output signal. Since mathematical models are nonlinear, the Newton – Kantorovich linearization and the matrix operator apparatus are applied to solve the problem. The features of gradient vector computation and features of the identification algorithm and its organization are considered. The results of the two time constants identification for the mathematical model of the turbine PT-12/15-35/10M are presented.
Kornyushin Y., Egupov N., Kornyushin P. Identification of parameters of regulator actuators for steam power turbines using matrix operator apparatus. Маthematical Modeling and Coтputational Methods, 2015, №2 (6), pp. 73-86
The article describes performed simulation of force action on streamlined horizontal elements of engineering structures in the upper layer of sharply stratified flow associated with the generation of waves at the interface between the liquid layers. We obtained an integral representation of the wave drag and lift, made numerical calculations for a real marine environment. The conditions under which there is a significant increase in the hydrodynamic reactions on streamlined structural elements were revealed.
Vladimirov I., Korchagin N., Savin A. Simulation of wave action on horizontal structure elements in the upper layer of stratified flow. Маthematical Modeling and Coтputational Methods, 2014, №4 (4), pp. 74-87
The article presents a suggested method of numerical finite-element solving the ‘hole ovalization’ problem. This method can be applied for experimental development of advanced aviation materials with the aim of determining structure element resistance against deforming with stress concentrators, mainly, connectors. The method is based on three-dimensional finite element solution of the problem of lastoplastic deformation of plates with a hole under crushing. It is appropriate for reduction of xperimental studies and replacing them by the numerical experiments. The Ilyushin model of small lastoplastic deformations has been used. The results of numerical simulation of a threedimensional stress-strain state of elastoplastic plates under crushing are presented as well as results of experimental nvestigations of deforming plates of Al-alloy 163. It is shown that the results of numerical and experimental modeling for deforming plates under crushing agree quite well.
Dimitrienko Y., Gubareva E., Sborschikov S., Erasov V., Yakovlev N. Computational modeling and experimental investigation of elastic-plastic plates deforming under crushing. Маthematical Modeling and Coтputational Methods, 2015, №1 (5), pp. 67-82
The article describes a method for calculating the dynamic stability of cylindrical shell under axial compressive time-varying load. The case of linearly varying load was con-sidered as an example.
Dubrovin V., Butina T. Modeling of the dynamic stability of a cylindrical shell under the axial compressive load. Маthematical Modeling and Coтputational Methods, 2015, №2 (6), pp. 46-57
The article presents an analytical expression for calculating pressure on the surface of blunted cones in a supersonic gas flow, taking into account the curvature discontinuity along the generatrix. We used a genetic algorithm and multi-stage functional optimisation methods for the least-squares method to determine free parameters of the expression. We compare the results obtained to the rigorous numerical solution to the inviscid problem. The comparison shows that it is possible to use the analytical expression for pressure distribution over a surface in a wide Mach number range for various cone halfangles. The expression proposed accounts for the curvature discontinuity along the generatrix at the point where the sphere is tangent to the conical surface, unlike the expressions found in previously published works.
Bulgakov V.N., Kotenev V.P., Sapozhnikov D.A. Modeling supersonic flow around blunted cones, taking into account the curvature discontinuity along the generatrix of the solid. Маthematical Modeling and Coтputational Methods, 2017, №2 (14), pp. 81-93