and Computational Methods

doi: 10.18698/2309-3684-2014-1-3656

The theory of thin constructive-orthotropic plates with a two-periodic structure was suggested. Examples of such structures are honeycomb sandwich panels and backed plates. The theory is based on equations of a three-dimensional elasticity theory with the help of asymptotic expansions in terms of a small parameter being the ratio of a plate thickness and a characteristic length without introducing 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 are similar to equations of the

Dimitrienko Y., Gubareva E., Sborschikov S. Asymptotic theory of constructive-orthotropic plates with two-periodic structures. Маthematical Modeling and Coтputational Methods, 2014, №1 (1), pp. 36-56

doi: 10.18698/2309-3684-2022-2-2862

The problem of development a theory for calculating the stress-strain state of thin multilayer elastic plates, for which linearized slip conditions are specified at the interface between the layers, is considered. The solution of this problem is constructed using an asymptotic analysis of the general equations of the 3-dimensional theory of elasticity with the conditions of non-ideal contact of the layers. The asymptotic analysis is carried out with respect to a small geometric parameter representing the ratio of the plate thickness to its characteristic length. Recurrent formulations of local quasi-one-dimensional problems of the theory of elasticity with slip are obtained. Explicit analytical solutions are obtained for these problems. The derivation of the averaged equations of elastic equilibrium of multilayer plates is presented, taking into account the slippage of the layers. It is shown that due to the effect of slippage of layers, the system of averaged equations of the theory of multilayer plates has an increased - 5th order of derivatives, in contrast to the classical 4th order, which takes place in the theory of Kirchhoff-Love plates. It is shown that the asymptotic theory makes it possible to obtain an explicit analytical expression for all 6 components of the stress tensor in the layers of the plate. As a special case, the problem of calculating the stress-strain state of a 4-layer plate under uniform pressure bending with one slip coefficient is considered. A complete analytical solution of this problem is obtained, including explicit expressions for all non-zero components of the stress tensor. A numerical analysis of the solution of the averaged problem for a composite plate is carried out, in which the layers are unidirectional reinforced fibrous materials oriented at different angles. A comparative analysis of the influence of the fiber reinforcement angles and the slip coefficient of the layers on the displacement of the plate and the distribution of stresses in the layers was carried out. It is shown that the problem of bending a plate with slip admits the existence of a spectrum of critical values of the slip coefficient, when passing through which the displacements and stresses in the layers of the plate change significantly, and these critical values depend on the angle of reinforcement of the composite layers.

Димитриенко Ю.И., Губарева Е.А. Асимптотическая теория многослойных тонких упругих пластин с проскальзыванием слоев. Математическое моделирование и численные методы, 2022, № 2, с. 30–64

doi: 10.18698/2309-3684-2014-4-1836

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

doi: 10.18698/2309-3684-2019-1-326

The problem of deformation of thin two-layer plates, for which a slip condition is speci-fied at the interface between the layers, instead of the classical case of ideal contact, is considered. The method of asymptotic analysis of the general equations of the 3-dimensional theory of elasticity is used to solve this problem under the influence of transverse pressure, longitudinal and shear forces on the end surfaces. Asymptotic analysis is performed using a small geometric parameter representing the ratio of thickness to the characteristic length of the plate. Recurrent formulations of local quasi-one-dimensional problems of elasticity theory with slippage are obtained. For these problems, explicit analytical solutions are obtained. The averaged equations of elastic equilibrium of a two-layer plate with slippage of layers are derived. It is shown that, due to slippage, the order of the averaged equations of the theory of plates increases to 5 orders of magnitude, in contrast to the classical 4th order, which takes place in the theory of Kirchhoff – Love plates. Additional boundary conditions to this 5th order system are formulated and its analytical solution is obtained for the case of a rectangular plate under the influence of uniform pressure. A numerical analysis of the solution of the averaged problem is carried out. It is shown that the presence of layer slippage significantly increases the deflection of the plate in comparison with the conditions of ideal contact of the layers.

Димитриенко Ю.И., Губарева Е.А. Асимптотическая теория тонких двухслой-ных упругих пластин с проскальзыванием слоев. Математическое моделирование и численные методы. 2019. № 1. с. 3–26.

doi: 10.18698/2309-3684-2015-1-6782

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

doi: 10.18698/2309-3684-2014-2-2848

We propose a method for calculating effective viscoelastic properties of composite materials under steady-state cyclical vibrations. The method is based on asymptotic averaging of periodic structures and finite-element solution of local problems of viscoelasticity in periodicity cells of composite materials. We provide examples of numerical simulation of viscoelastic properties for composites with unidirectional reinforcement, and of calculations of complex tensors of stress concentration in a periodicity cell. The paper presents a comparative analysis of dependencies of loss tangent of complex composite elasticity

modulus on vibration frequencies obtained through FEA calculations and rough mixed formulae. We show that rough mixed formulae, often used for calculating dissipative properties of composite materials, can yield appreciable calculation errors.

Dimitrienko Y., Gubareva E., Sborschikov S. Finite element modulation of effective viscoelastic properties of unilateral composite materials. Маthematical Modeling and Coтputational Methods, 2014, №2 (2), pp. 28-48

doi: 10.18698/2309-3684-2017-1-3254

The article considers the modeling results of incompressible layered composites with finite strains deformation according to the individual layers characteristics. The article proposes an asymptotic averaging method version for layered nonlinearly elastic incompressible composites with finite deformations and periodic structure. We are using a universal representation of the defining relations for incompressible composite layers, proposed by Yu.I. Dimitrienko, which allows us to simulate simultaneously for a complex of various nonlinear elastic media models characterizedby the choice of a pair of energy tensors. We proved that if all composite layers are incompressible, the composite as a whole is also an incompressible, but anisotropic, medium. The article considers the problem of laminated plate uniaxial stretching from incompressible layers with finite deformations. Using the developed method, we calculated the effective deformation diagrams connecting the averaged Piola — Kirchhoff stress tensors components and the strain gradient, as well as the stress distribution in the composite layers.

The developed method for calculating effective deformation diagrams and stresses in composite layers can be used in the design of elastomeric composites with specified properties.

Dimitrienko Y., Gubareva E., Kolzhanova D., Karimov S. Incompressible layered composites with finite deformations on the basis of the asymptotic averaging method. Маthematical Modeling and Coтputational Methods, 2017, №1 (13), pp. 32-54

doi: 10.18698/2309-3684-2020-2-2645

The paper is devoted to the development of a method for calculating the nonlinear dielectric properties of composites with a periodic structure. Methods for predicting of the nonlinear dielectric properties of composites play an important role in the design of dielectric materials with specified properties, in particular for heterogeneous ferroelectrics, which are widely used to create various devices and electrical devices, for example, to create memory storage devices for computers. A quasi-static problem of the distribution of an electric charge in an inhomogeneous polarizable medium with a periodic structure and nonlinear dielectric properties is considered. To solve this nonlinear problem, the asymptotic homogenization method proposed by N.S. Bakhvalov, E. Sanchez-Palencia, B.E. Pobedria. As a result, local nonlinear problems of electrostatics on the periodicity cell are formulated, an algorithm for calculating effective nonlinear constitutive relations for dielectric properties, and an averaged problem for a composite with effective properties are proposed. For the case of a composite with a layered structure, the solution of local problems is obtained, and effective defining relations for the nonlinear dielectric properties of the composite are constructed. It is shown that a laminated composite is a transversely isotropic nonlinear dielectric material if it is isotropic materials. A numerical example of calculating the nonlinear properties of a 2-layer composite based on barium titanate and ferroelectric ceramic varicond VK4 is considered. A model is proposed that describes the nonlinear dependence of the dielectric constant of these materials on the vector of the electric field strength. It is shown that the nonlinear dependence of the dielectric constant tensor of the composite on the strength vector differs significantly for the direction of the field in the plane of the layers and in the transverse direction. It is shown that the developed technique can serve as a basis for designing nonlinear dielectric composite materials with anisotropic properties.

Димитриенко Ю.И., Губарева Е.А., Зубарев К.М. Моделирование нелинейных диэлектрических свойств композитов на основе метода асимптотической гомогенизации. Математическое моделирование и численные методы. 2020. № 2. с. 26–45

doi: 10.18698/2309-3684-2021-1-330

The problem of calculating the stress–strain state of a composite laminated panel during cylindrical bending under conditions of finite deformations is considered. To solve the problem, the method of asymptotic averaging of periodic nonlinear elastic structures with finite deformations was applied, which was developed in detail earlier in the previous works of the authors. A feature of this problem is the use of universal models of constitutive relations for isotropic components of the composite, as well as for the composite as a whole, which is a transversely isotropic nonlinear elastic medium. Universal models make it possible to obtain solution of problems within the framework of a single solution algorithm simultaneously for several classes of models of nonlinear elastic media corresponding to different conjugate pairs of stress tensors–deformation. An analytical solution is obtained for the problem of cylindrical bending of a composite panel. A numerical analysis of the solution is carried out using the example of a composite, the periodicity cell of which consists of two layers: polyurethane and rubber. It is shown that for thin panels the stresses, both averaged and true, practically do not depend on the class of the model of the constitutive relations. At the same time, for thicker panels, the stresses differ significantly for different classes of models of composite layers.

Димитриенко Ю.И., Губарева Е.А., Каримов С.Б., Кольжанова Д.Ю. Моделирование напряжений в композитной нелинейно упругой панели при цилиндрическом изгибе. Математическое моделирование и численные методы, 2021, № 1, с. 3–30.

doi: 10.18698/2309-3684-2022-4-330

The general asymptotic theory of thin multilayer shells developed by the authors earlier in Part 1 of this study is applied to cylindrical anisotropic thermoelastic shells. It is shown that for cylindrical shells the general theory is substantially simplified: general two-dimensional averaged thermoelasticity equations for multilayer shells are obtained. These equations are similar to the classical equations of cylindrical shells in the Kirchhoff-Love theory, but they are obtained in a completely different way: on the basis of only an asymptotic analysis of the general three-dimensional equations of the theory of thermoelasticity. No hypotheses regarding the distribution of displacements or stresses over the thickness are used in this theory, which makes it logically consistent. In addition, the developed theory makes it possible to obtain explicit analytical expressions for all 6 components of the stress tensor in cylindrical anisotropic shells. Explicit expressions are obtained for all tensor constants included in these stress formulas. An example of calculating thermal stresses in a cylindrical composite shell with axisymmetric bending due to the combined action of external pressure and one-sided non-stationary heating is given. An example of a layered-fiber 4-layer shell with different angles of helical winding of reinforcing fibers is considered. It is shown that the developed one allows one to study in detail such complex effects as the formation of significant transverse thermal stresses during heating, which significantly exceed the level of interlayer shear stresses, which are traditionally considered the most dangerous for layered composites.

Димитриенко Ю.И., Губарева Е.А., Пичугина А.Е., Белькова К.В., Борин Д.М. Моделирование термонапряжений в композитных оболочках на основе асимптотической теории. Часть 2. Расчет цилиндрических оболочек. Математическое моделирование и численные методы, 2022, № 3, с. 3–30

doi: 10.18698/2309-3684-2018-3-114132

The previously developed general asymptotic theory of thin multilayer shells is used for the case of cylindrical shells. The ratios are presented in explicit analytical form for all six components of the stress tensor in a thin multilayer elastic cylindrical shell, depending on the deformations, curvatures of the middle surface of the shell, as well as their derivatives along the longitudinal coordinates. The obtained formulas make it possible to calculate all the distributions of the components of the stress tensor over the thickness in a cylindrical shell after finding solutions to the two-dimensional problem of the theory of KirchhoffLyav shells. An example is given of the calculation of stresses in a cylindrical composite shell underaxisymmetric bending by pressure. To calculate stresses by these formulas, only a differentiation of displacements is required - a deflection and two displacements of the middle surface of the shell, for which an analytical solution is obtained.

Димитриенко Ю.И., Губарева Е.А., Пичугина А.Е. Моделирование напряжений в тонких композитных цилиндрических оболочках на основе асимптотической теории. Математическое моделирование и численные методы, 2018, № 3, с. 114–132.

doi: 10.18698/2309-3684-2021-2-1537

The aim of this work is to find the constitutive relations for a layered elastoplastic composite according to the flow theory using the method of asymptotic averaging. This goal is achieved by developing an algorithm for solving the problem of the theory of plastic flow for a layered composite material, taking into account various characteristics and properties of these layers of the material, followed by visualizing the result in the form of effective plasticity diagrams connecting the components of averaged stress tensors and components of averaged strain tensors.

Димитриенко Ю.И., Губарева Е.А., Черкасова М.С. Моделирование деформирования слоистых периодических композитов на основе теории пластического течения. Математическое моделирование и численные методы, 2021, № 2, с. 15–37.

doi: 10.18698/2309-3684-2018-4-7292

The paper considers a model of effective constitutive relations for a transversal-isotropic incompressible composite with finite strains. The model belongs to the so-called class of universal models that connect several pairs of energy stress and strain tensors simultaneously. A method is proposed for separating coupled problems of micro- and macroscopic deforming composites with finite strains that arise when the method of asymptotic homogenization (AH) of periodic structures is used. The method is based on the application of the effective constitutive relation model as an approximation dependence of the results of numerical simulation of the composite deformation curves obtained using the exact AH method. To find the elastic constants of the transversely isotropic composite model the method of minimizing the deviation of the approximation of deformation diagrams from the AH diagrams is used for a series of problems of standard deforming at finite strains. To solve minimization problems, the Nelder—Mead method is used. The results of numerical simulation by the proposed method for nonlinear elastic layered composites are presented, which showed good approximation accuracy, achieved due to application of the proposed method for the separation of coupled problems of micro- and macroscopic deforming.

Димитриенко Ю.И., Губарева Е.А., Каримов С.Б., Кольжанова Д.Ю. Модели-рование эффективных характеристик трансверсально-изотропных несжимаемых композитов с конечными деформациями. Математическое моделирование и чис-ленные методы, 2018, № 4, с. 72–92.

doi: 10.18698/2309-3684-2015-2-322

The article presents a model of microstructure of two-phase monocrystalline intermetallic alloys in the form of a periodic structure of the hexagonal type, as well as a mathematical model of elastic-plastic deformation of monocrystalline alloy, based on the method of asymptotic smoothing periodic structures. Deformation plasticity theory under loading is used for the phases with due regard for the effect of their damage level during loading. For numerical calculations of the developed model the heat-resistant monocrystalline alloy of the type VKNA-1V was used. Finite element calculations of deformation and fracture micromechanical processes in monocrystalline alloy of the type VKNA-1V were carried out. It was found that under tension maximum values of phase damagability, which determine the beginning of the alloy micro-fracture zone, are achieved in the areas adjacent to the phase interface and in areas of maximum curvature of the geometric shape of the phases. Calculations of heat-resistant alloy strain diagrams in plastic range are proved to be consistent with experimental data.

Dimitrienko Y., Gubareva E., Sborschikov S., Bazyleva O., Lutsenko A., Oreshko E. Modeling the elastic-plastic characteristics of monocrystalline intermetallic alloys based on microstructural numerical analysis. Маthematical Modeling and Coтputational Methods, 2015, №2 (6), pp. 3-22

doi: 10.18698/2309-3684-2016-2-323

The purpose of this article is to propose a model of deformation of elastic-plastic composite materials with periodic structures with an allowance for fault probability of the composite phases. The model is based on a variant of the deformation theory of plasticity with the active loading. To simulate the effective characteristics of elastic-plastic composites, we applied the method of asymptotic homogenization of periodic structures. For numerical solution of linearized problems on the periodicity cell we offered the finite elements method using SMCM software medium developed at the Scientific-Educational Center of Supercomputer Engineering Modeling and Program Software Development of the Bauman Moscow State Technical University. We provide the research with the examples of numerical computations for dispersion-reinforced metal composites (aluminum matrix filled with SiC particles). Finally, we present the results of numerical modeling of deformation processes, damage accumulation and metal-composite destruction.

Dimitrienko Y., Gubareva E., Sborschikov S. Multiscale modeling of elastic-plastic composites with an allowance for fault probability. Маthematical Modeling and Coтputational Methods, 2016, №2 (10), pp. 3-23

doi: 10.18698/2309-3684-2014-3-7488

We examined effects of typical for different climatic zones atmospheric conditions on flight program optimization for a subsonic long-haul passenger aircraft. Simulation of flight and power plant performance was based on current traditional approaches used in solving problems of this kind. The acceleration-climb flight segment has been optimized by minimizing fuel consumption at this flight segment. The cruising flight segment has been optimized considering operating limitations accepted for civil aviation. The in-built model of bypass turbojet engine was used for simulating the flight. This model allows calculating power plant performances under any flight conditions. The flight of subsonic aircraft has been examined in one vertical plane. Calculations have been performed for 6 standard air temperature variations with altitude (depending on climatic zone). Atmospheric pressure variation near Earth surface was considered and effects of atmospheric conditions on flight program optimization were estimated.

Mozzhorina T., Gubareva E. Simulating atmospheric conditions influence on flight program optimization for a subsonic passenger aircraft. Маthematical Modeling and Coтputational Methods, 2014, №3 (3), pp. 74-88

doi: 10.18698/2309-3684-2020-4-84110

An asymptotic theory of thermoelasticity of multilayer composite shells is proposed, the derivation of the basic equations of which is based on the asymptotic expansion in terms of a small geometric parameter of three-dimensional thermoelasticity equations. This method was previously developed by the authors for thin composite plates, and in this article it is applied to thin-walled shells of an arbitrary frame. According to the developed method, the original three-dimensional problem of thermoelasticity decomposes into a recurrent successor of one-dimensional local problems of thermoelasticity and an averaged two-dimensional problem of thin shells. For local problems of thermoelasticity, analytical solutions are obtained, which make it possible to close the averaged formulation of the problem of the theory of shells with respect to 5 unknown functions: longitudinal displacements, deflection, and two shear forces. It is shown that the averaged problem for multilayer shells coincides with the classical system of equations for Kirchhoff–Love shells, however, it is more substantiated, since the asymptotic theory does not contain any assumptions regarding the pattern of the distribution of permutations and stresses over thickness. In addition, the asymptotic theory makes it possible to calculate all the stresses in the shell, without solving any additional problems, but only by differentiating the averaged displacements.

Димитриенко Ю.И., Губарева Е.А., Пичугина А.Е. Моделирование термона-пряжений в композитных оболочках на основе асимптотической теории. Часть 1. Общая теория оболочек. Математическое моделирование и численные методы, 2020, № 4, с. 84–110.