Vol 29, No 4 (2023)

Cover Page

Full Issue

Mechanics

Influence of the process of damage accumulation on the asymptotic behavior of stress fields under creep conditions of a sample with a central crack

Bykova Y.S., Stepanova L.V.

Abstract

The present study describes the influence of the mutual effect of damage accumulation on the stress fields at the central crack tip in the plate subjected to uniaxial tension. The objective of the study is to analyze the cracked plate experiencing uniaxial loading under creep conditions using the CAE software SIMULIA ABAQUS taking into account the damage accumulation processes near the crack tip. Computations were performed by means of the user procedure UMAT (User material), which is based on the Bailey-Norton creep power law and the Kachanov — Rabotnov damage evolution equation describing the power-law damage accumulation processes. The analysis of the obtained results showed that in the creep and elasticity zones, in computations without taking into account the damage effect, there are asymptotics of stress fields that correspond to well-known analytical solutions of fracture mechanics (Hutchinson — Rice — Rosengren asymptotics for the creep zone, the asymptotics of linear fracture mechanics corresponding to the stress dependence. inversely proportional to the square root of the distance from the tip of the crack – for the zone of elastic behavior of the material. The presence of damage in the cracked specimen affected the asymptotic behavior of the crack tip fields. It is felt that the proposed procedure could pave the way for the analytical solution of the boundary value problem and allow us to determine the structure of the asymptotic solution of the problem. The finite element analysis clearly shows that the process of damage accumulation affects the change in the asymptotic behavior of the stress field in the vicinity of the crack tip and leads to a new asymptotic distribution of the stress tensor components. The proposed procedure can pave the way to an analytical solution of the boundary value problem and allow us to determine the structure of the asymptotic solution of the problem.

Vestnik of Samara University. Natural Science Series. 2023;29(4):7-25
pages 7-25 views

Reference shape of bodies with enhanced kinematics. Part I. Geometric methods

Koifman K.G.

Abstract

The work develops differential-geometric methods for modeling finite incompatible deformations of hyperelastic solids. They are based on the representation of a body as a smooth manifold, on which a metric and a non-Euclidean connection are synthesized. The resulting geometric space is interpreted as global stress-free shape, and the physical response and material balance equations are formulated relative to it. Within the framework of the geometric approach, deformations are modeled as embeddings of a non-Euclidean shape in physical space. Measures of incompatibility are represented by invariants of the affine connection, namely, curvature, torsion and nonmetricity, and the connection itself is determined by the type of physical process.

This article is the first part of the study. The proposed geometric approach is applied to bodies whose response depends on the first deformation gradient. Compatibility conditions are obtained and their geometric interpretation is proposed.

Vestnik of Samara University. Natural Science Series. 2023;29(4):26-53
pages 26-53 views

Reference shape of bodies with enhanced kinematics. Part II. Second gradient and microstructure

Koifman K.G.

Abstract

The work develops differential-geometric methods for modeling finite incompatible deformations of hyperelastic solids with enhanced kinematics. The response of such bodies, along with the standard kinematic field represented by the deformation gradient, is characterized by additional tensor fields. As such, the paper considers: 1) the second deformation gradient and 2) the tensor field of the second rank, modeling the microstructure of the body. For each of these two cases, compatibility conditions are obtained and their geometric interpretation is proposed. Geometry is synthesized on the material manifold representing a body with enhanced kinematics. The corresponding affine connection has non-zero torsion and curvature, which can be useful for modeling a body with dislocations and disclinations.

Vestnik of Samara University. Natural Science Series. 2023;29(4):54-76
pages 54-76 views

Bending of a circular disk: from cylinder to ultrathin membrane

Lychev S.A., Digilov A.V., Pivovaroff N.A.

Abstract

The article discusses methods of mathematical modeling of the stress-strain state of a circular disc at various ratios of its thickness to radius, ranging from 1 to 103. For sufficiently thick plates, the solution of three-dimensional linear elasticity theory is used, for plates of medium thickness — the solution of linear bending equations within the Kirchhoff – Love hypotheses and nonlinear equations of Foppl – von Karman, and for ultrathin plates — the nonlinear equations of Adkins – Rivlin – Green. A comparative analysis of the solutions has been conducted, and ranges of relative thickness have been identified in which the considered solutions adequately describe the deformation process. This result enables the selection of a method for mathematical modeling of the stress-strain state of circular plates used in microelectromechanical systems that is most suitable for their relative size.

Vestnik of Samara University. Natural Science Series. 2023;29(4):77-105
pages 77-105 views

Mathematical Modelling

Experience in modeling inclined cracks in materials with cubic crystal structure

Mushankova K.A., Stepanova L.V.

Abstract

In this work, a good coincidence of atomistic and continuum stress fields at the crack tip under mixed mode loading conditions in an anisotropic medium with cubic symmetry of elastic properties is revealed. The atomic stress distributions associated with the crack tip are obtained using the molecular dynamics method. Continuum distributions are obtained from the theoretical solution of the problem of determining the stress-strain state at the crack tip, based on the methods of the elasticity theory of anisotropic media and the subsequent decomposition of complex potentials by eigenfunctions. In the framework of a molecular dynamics computational experiment, a single-crystal face-centered copper at low temperature was considered in order to isolate the elastic mode of deformation of a single crystal, and the embedded atom potential was used. A distinctive feature of the conducted molecular dynamic modeling is the consideration of a crack that makes up various angles with planes of symmetry of the crystal. In the vicinity of the crack tip, points lying in annular regions at different distances from the crack tip and of different thickness were selected, and the dependences of the stress tensor components depending on the polar angle were plotted. A comparison of the angular dependencies obtained by atomistic calculation and using a theoretical solution showed their good consistency. It is found that the similarity of the angular dependences of the stress tensor components is observed for all the studied values of two angles: the angle between the axis of symmetry of the crystal lattice (in the plane of the plate) and the direction of the crack and the angle between the action of the tensile load and the crack line.

By virtue of this property of solutions, it can be concluded that solutions of continuous fracture mechanics can serve to describe stress fields at atomistic distances from the crack tip.

Vestnik of Samara University. Natural Science Series. 2023;29(4):106-116
pages 106-116 views

Mathematical Methods in Natural Sciences

Modeling of kinematic characteristics of the line of sight of a mobile antenna of a spacecraft when it is pointed at an information reception point

Gorelov Y.N., Kurganskaya L.V.

Abstract

Within the framework of the general problem of calculating the kinematic characteristics of the mobile antenna of the spacecraft, the problem of the kinematics of the complex movement of the line of sight "Sc — IRP" ("spacecraft — information reception point") is considered. The solution of this problem is given, with the help of which the relationship of kinematic characteristics in the absolute and relative movements of the line of sight "Sc — IRP" is established in the form of theorems on the addition of its angular velocities and accelerations.

Vestnik of Samara University. Natural Science Series. 2023;29(4):117-124
pages 117-124 views

Modeling of kinematic characteristics of a mobile antenna of a spacecraft in a two-stage pivoting device

Gorelov Y.N., Kurganskaya L.V.

Abstract

The problem of calculating the kinematic characteristics of a mobile antenna device (AD) based on the specified values of kinematic characteristics of the line of sight (LS) "spacecraft — information reception point" in the associated coordinate system of the spacecraft is considered. The kinematic characteristics of the AD in the form of the current values of its rotation angles and angular velocities and accelerations (according to the corresponding control channels) are determined from the conditions for combining the LS with the axis of the directional diagram of the AD for the case of a two-stage pivoting device with mutually orthogonal axes of rotation of the AD.

Vestnik of Samara University. Natural Science Series. 2023;29(4):125-132
pages 125-132 views

Physics

Modeling of nanoheterojunction betavoltaic cell on GaN and GaP with Si and 3C-SiC/Si

Dolgopolov M.V., Chipura A.S.

Abstract

In this paper, the electrophysical properties and efficiency of energy converters — betavoltaic elements, which contain a GaN and GaP heterojunction on Si and 3C-SiC/Si substrates, are modeled. For conversion into electrical energy, external 63Ni or internal 14C radioactive sources with a test specific activity of 100 mCicm-2 are investigated in the simulation. The system of parameters and characteristics is optimized: diffusion lengths, short-circuit current, open circuit voltage, filling factor, reverse saturation current and efficiency. It was shown in simulation results, that in the device structure with junction depth of 0.1 microns, the good operation of a betavoltaic element is determined, the short-circuit current density is up to 200 nAcm-2, the open circuit voltage is up to 3.7 V, the power density is up to 700 nWcm-2, efficiency up to 25 %. The conversion efficiency reaches its maximum value when using a radioisotope source with an activity density from 25 to 100 mCicm-2. The conversion efficiency with the location of the injector source inside is estimated to be about 30 times higher than with the external location.

Vestnik of Samara University. Natural Science Series. 2023;29(4):133-142
pages 133-142 views

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies