Next Article in Journal
Modeling 3D–1D Junction via Very-Weak Formulation
Next Article in Special Issue
Displacements and Stress Functions of Straight Dislocations and Line Forces in Anisotropic Elasticity: A New Derivation and Its Relation to the Integral Formalism
Previous Article in Journal
Interval Valued T-Spherical Fuzzy Soft Average Aggregation Operators and Their Applications in Multiple-Criteria Decision Making
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spherically Symmetric Tensor Fields and Their Application in Nonlinear Theory of Dislocations

by
Evgeniya V. Goloveshkina
* and
Leonid M. Zubov
Institute of Mathematics, Mechanics, and Computer Science of Southern Federal University, Milchakova Str. 8a, 344090 Rostov-on-Don, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(5), 830; https://doi.org/10.3390/sym13050830
Submission received: 12 April 2021 / Revised: 28 April 2021 / Accepted: 6 May 2021 / Published: 9 May 2021
(This article belongs to the Special Issue Recent Advances in the Study of Symmetry and Continuum Mechanics II)

Abstract

:
The concept of a spherically symmetric second-rank tensor field is formulated. A general representation of such a tensor field is derived. Results related to tensor analysis of spherically symmetric fields and their geometric properties are presented. Using these results, a formulation of the spherically symmetric problem of the nonlinear theory of dislocations is given. For an isotropic nonlinear elastic material with an arbitrary spherically symmetric distribution of dislocations, this problem is reduced to a nonlinear boundary value problem for a system of ordinary differential equations. In the case of an incompressible isotropic material and a spherically symmetric distribution of screw dislocations in the radial direction, an exact analytical solution is found for the equilibrium of a hollow sphere loaded from the outside and from the inside by hydrostatic pressures. This solution is suitable for any models of an isotropic incompressible body, i.e., universal in the specified class of materials. Based on the obtained solution, numerical calculations on the effect of dislocations on the stress state of an elastic hollow sphere at large deformations are carried out.

1. Introduction

In this work, a spherically symmetric state in which the dislocation density tensor has a spherically symmetric representation is determined. Herewith, the system of equations consisting of constitutive equations, equilibrium equations, and an incompatibility equation admits a spherically symmetric solution. Despite the simplicity of representation, such solutions help to describe complex physical and mechanical phenomena. In gas dynamics, for the problem of the propagation of a strong blast, spherical symmetry enables one to consider the charge as a point charge and is used to find the adiabatic unsteady motions of a perfect gas [1]. In the paper [2], a static spherically symmetric solution for the gravitational field equations obtained in a spherically symmetric Kalb–Ramond background allows us to conclude that in string gravity, spacetime can have such features as a wormhole or a naked singularity. The work [3] is devoted to static spherically symmetric spacetime with nonlinear electrodynamic sources, where conditions under which the metric of static spherically symmetric solutions of the Einstein equations associated with nonlinear electrodynamics is regular (without singularities) at the center (origin) are established. In [4], it is investigated which representations of the Higgs fields of the gauge group can be spherically symmetric or at least contribute to spherical symmetry of stress-energy tensor, and a simple criterion for the existence of nontrivial spherically symmetric Higgs fields is obtained. Two classes of spherically symmetric static vacuum solutions of the Poincaré gauge field theory, obtained using the double duality ansatz for the corresponding Riemann–Cartan curvature tensors are considered in [5]. With the help of the Petrov classification scheme and the Einstein’s field equations using a radial coordinate dependent symmetric tensor of the second-rank, invariant with respect to inflation in the radial direction, in [6], it is investigated what the algebraic structure of the cosmological term can be in the spherically symmetric case. In the paper [7], for an arbitrary number of dimensions with arbitrary torsion it is shown that in the most general theory of gravity, only static spherically symmetric asymptotically flat solutions of the variational equations of motion contain gravitational singularities.
Microstructure defects in the form of dislocations began to be investigated at the beginning of the twentieth century. Some works [8,9,10,11,12,13,14,15,16,17,18,19] devoted to this topic are based on the continuum theory of dislocations, which assumes a transition from a discrete set of dislocations to their continuous distribution. It is advisable to carry out such a transition in the case when the number of dislocations in a limited volume is large enough.
This paper solves the problem for an elastic sphere containing continuously distributed dislocations. Problems about a sphere with dislocations were previously considered in the papers [20,21,22,23]. The main result of this work is a new exact solution to the nonlinear theory of dislocations, which is of spherical symmetry. This solution is universal in the class of isotropic incompressible materials. Universal solutions are such solutions of the equilibrium equations that are valid for any constitutive equations from a certain class of materials. These solutions are used to experimentally determine the constitutive equations and to check numerical results. For incompressible isotropic nonlinear elastic bodies, five families of nonuniform universal deformations are known [24,25,26,27,28,29,30,31,32]. Universal deformations of solids with constraints different from the incompressibility condition are studied in [33,34,35]. Universal solutions for generalized models of continuous media with complex physical and mechanical properties [36,37,38,39] are built up in [40,41].
The obtained universal solution supplements the small list of known exact solutions of the nonlinear dislocation theory [16,17,18,19,42] and describes the influence of distributed screw dislocations in the radial direction on large spherically symmetric deformations of a hollow elastic sphere made of an incompressible material in the presence of an external or internal hydrostatic pressure. The effect of dislocations on the stress state of thin and thick spherical shells is analyzed.

2. Input Relations of the Theory of Dislocations

Let R = X k i k and r = x s i s be the radius vectors of the point of the elastic body in the reference and deformed configurations, respectively, X k and x s ( k , s = 1 , 2 , 3 ) are the Cartesian coordinates of the reference and final states of the body, i k are the coordinate unit vectors. In what follows, we will use the definitions of the gradient, divergence, and curl operators in the coordinates of the reference configuration
Grad Ψ = R N Ψ Q N , Div Ψ = R N · Ψ Q N , Curl Ψ = R N × Ψ Q N , R N = i k Q N X k .
Here Ψ is an arbitrary differentiable tensor field, Q N = Q N X 1 , X 2 , X 3 , N = 1 , 2 , 3 are some Lagrangian curvilinear coordinates, ⊗ denotes tensor multiplication, × and · denote vector multiplication and scalar multiplication, respectively.
The system of equilibrium equations for a nonlinear elastic medium in the absence of mass forces consists of the equilibrium equations
Div D = 0 ,
constitutive equations
D = d W d F , W = W ( G ) ,
and geometric relationships
F = Grad r , G = F · F T .
Here D is the asymmetric Piola stress tensor, also called the first Piola - Kirchhoff stress tensor, F is the deformation gradient, G is the metric tensor, W is the specific strain energy.
Let V be the region occupied by an elastic body in the reference configuration. To introduce the concept of dislocation density in an elastic medium, consider the problem of determining the displacement field u ( R ) = r ( R ) R from the tensor field F ( R ) given in the multiply connected domain V, which is assumed to be differentiable and single-valued in the domain V and satisfies the compatibility equation Curl F = 0 . Taking into account that
Grad u = F I ,
where I is the unit tensor, we see that in the case of a multiply connected domain the vector field u ( R ) , and hence the vector field r ( R ) , generally speaking, are not uniquely determined. This means the presence of isolated translational dislocations in the body [15], each of which, due to (5), is characterized by the Burgers vector b n :
b n = Γ n d R · F I = Γ n d R · F , n = 1 , 2 , . . . n 0 .
Here, Γ n is a simple closed contour encompassing the line of only one n-th dislocation. The total Burgers vector of a discrete set of n 0 dislocations according to (6) is determined by the formula
B = n = 1 n 0 b n = n = 1 n 0 d R · F .
Due to the well-known properties of curvilinear integrals, in (7) it can be replaced by a single integral over the closed contour Γ 0 covering the lines of all n 0 dislocations
B = Γ 0 d R · F .
Following [43,44,45], we pass to the limit from a discrete set of dislocations to their continuous distribution and transform the contour integral (8) into a surface integral by the Stokes formula
B = Σ 0 N · Curl F d Σ .
Here, Σ 0 is the surface spanning the contour Γ 0 , N is the unit normal to Σ 0 . The relation (9) allows us to introduce [43,45,46] the density of continuously distributed dislocations as a second-rank tensor α = Curl F , the flow of which through any surface is equal to the total Burgers vector of dislocations crossing this surface. In what follows, the dislocation density tensor is assumed to be a given function of the Lagrangian coordinates Q N , which must satisfy the solenoidality requirement
Div α = 0 .
If dislocations with a given tensor density are distributed in a body, then the displacement field and the vector field r ( R ) do not exist. In this case, the relation (4) is replaced by the incompatibility equation
Curl F = α ,
and the tensor F is called the distortion tensor.
The full system of equilibrium equations for an elastic body with distributed dislocations contains the components of the tensor distortion field F as unknown functions and consists of the equilibrium Equation (2), the incompatibility Equation (11), and the constitutive Equation (3).
If a given hydrostatic pressure q acts on a part of the body surface with normal N , then the boundary conditions on this part of the boundary have the form [24]
N · D = q ( det F ) F 1 · N .
If an elastic material is isotropic, then the specific energy W is the function of invariants of the metric tensor G
W ( G ) = W ( I 1 , I 2 , I 3 ) ,
I 1 = tr G , I 2 = 1 2 tr 2 G tr G 2 , I 3 = det G = ( det F ) 2 .
In the following, we will consider an isotropic incompressible material. In this case, the specific energy is a function of only two arguments: W = W ( I 1 , I 2 ) , since the third invariant I 3 is equal to unity due to the incompressibility condition
det F = 1 .
The constitutive equation of an isotropic incompressible elastic body has the form [24,31,32]
D = D * p F T , D * = ( τ 1 + I 1 τ 2 ) F τ 2 G · F ,
τ 1 = 2 W ( I 1 , I 2 ) I 1 , τ 2 = 2 W ( I 1 , I 2 ) I 2 .
Here, p is a pressure in an incompressible material, not expressed in terms of deformation, τ β ( I 1 , I 2 ) ( β = 1 , 2 ) are the scalar response functions of an incompressible elastic medium.
Sometimes the elastic properties of an isotropic incompressible material are specified using the function of the first and second invariants of the tensor U [24,32]:
W ( G ) = W ( J 1 , J 2 ) ,
J 1 = tr U , J 2 = 1 2 tr 2 U tr U 2 , U = G 1 / 2 .
The tensor D * in the constitutive Equation (16) is expressed through the specific energy (18) as follows [24]
D * = ( η 1 + J 1 η 2 ) A η 2 F , η β = W ( J 1 , J 2 ) J β , β = 1 , 2 .
Here, A is the rotation tensor. The symmetric positive definite tensor U and the properly orthogonal tensor A form the polar decomposition of the distortion tensor
F = U · A .

3. Spherically Symmetric Tensor Fields

In the reference configuration of the elastic body, let us introduce the spherical coordinates Q 1 = Φ (longitude), Q 2 = Θ (latitude), Q 3 = R (radial distance) by the formulas
X 1 = R cos Φ cos Θ , X 2 = R sin Φ cos Θ , X 3 = R sin Θ .
In what follows, we will use the orthonormal vector basis e Φ , e Θ , and e R consisting of unit vectors tangent to the coordinate lines:
e Φ = i 1 sin Φ + i 2 cos Φ , e Θ = i 1 cos Φ + i 2 sin Φ sin Θ + i 3 cos Θ , e R = i 1 cos Φ + i 2 sin Φ cos Θ + i 2 sin Θ .
Consider a second-rank tensor field P and represent it in the form of the decomposition
P = P m n e m e n , m , n = 1 , 2 , 3 , e 1 = e Φ , e 2 = e Θ , e 3 = e R .
A tensor field P will be called spherically symmetric if the components P m n depend only on the radial coordinate R and it is invariant under any rotations around the radial axis, i.e., around the ort e R . The latter means [47,48] that for any real values of χ ( R ) , the following equality holds:
K · P · K T = P , K = g cos χ ( R ) + d sin χ ( R ) + e R e R ,
g = e Φ e Φ + e Θ e Θ = I e R e R , d = e Φ e Θ e Θ e Φ = I × e R .
Let us find a general representation of the second-rank spherically symmetric tensor field. In (24), set χ π . Then K = 2 e R e R I and the relation (24) gives the equality
( 2 e R e R I ) · P = P · ( 2 e R e R I ) .
From (26) it follows that e R · P · e Φ = e Φ · P · e R = e Θ · P · e R = e R · P · e Θ = 0 . Setting χ π / 2 in (24), we obtain the equalities e Φ · P · e Φ = e Θ · P · e Θ , e Φ · P · e Θ = e Θ · P · e Φ . Thus, the tensor P has the form
P = P 1 ( R ) g + P 2 ( R ) d + P 3 ( R ) e R e R .
With the help of (25) it is easy to check that the tensors g , d and e R e R satisfy the relation (24) for any values of χ . Therefore, the representation (27) is necessary and sufficient for spherical symmetry of the tensor field P .
Let us point out some elementary but useful facts related to second-rank spherically symmetric tensors.
1. Two spherically symmetric tensors P and L commute: P · L = L · P .
2. If P 1 2 + P 2 2 P 3 0 , then there is an inverse tensor P -1 which is expressed like this:
P -1 = P 1 P 1 2 + P 2 2 g P 2 P 1 2 + P 2 2 d + 1 P 3 e R e R .
3. The polar decomposition, i. e., the representation in the form P = H · S , where H is a symmetric positive definite tensor, and S is an orthogonal tensor, for a spherically symmetric tensor field is written in a simple and explicit form:
H = P 1 2 + P 2 2 g + P 3 e R e R , S = P 1 P 1 2 + P 2 2 g + P 2 P 1 2 + P 2 2 d + sgn P 3 e R e R .
Arbitrary rank tensor field B = B k l p q e k e l e p e q is called spherically symmetric if its components B k l p q depend only on the radial coordinate and for any values of χ the following equality holds:
B k l p q ( R ) ( e k · K ) ( e l · K ) ( e p · K ) ( e q · K ) = B .
4. An arbitrary spherically symmetric vector field has the form f ( R ) e R .
5. Using the concept of an isotropic tensor function [47,48] one can prove that the value of an isotropic function of any number of spherically symmetric tensor arguments of arbitrary rank is itself a spherically symmetric tensor. In particular, the product of two second-rank spherically symmetric tensors is also a second-rank spherically symmetric tensor.
6. With the help of Formulas (1), (22), and (27), we can derive the expressions for the gradient, curl, and divergence of a second-rank spherically symmetric tensor field P
Grad P = e R d P 1 d R g + d P 2 d R d + d P 3 d R e R e R + P 3 P 1 R e Φ e Φ e R + e R e Φ + e Θ e Θ e R + e R e Θ + P 2 R e Φ e Θ e R e R e Θ + e Θ e R e Φ e Φ e R ,
Curl P = 1 R d d R ( R P 2 ) g + P 3 P 1 R d P 1 d R d + 2 P 2 R e R e R ,
Div P = d P 3 d R + 2 R ( P 3 P 1 ) e R .
Based on the definition (30), using (31) one can check that the third rank tensor field Grad P is spherically symmetric.

4. Spherically Symmetric State of Elastic Hollow Sphere with Distributed Dislocations in Radial Direction

Consider an elastic body in the form of a hollow sphere with an outer radius R 0 and an inner radius R 1 . Let us assume that the dislocation density tensor is of spherical symmetry, i.e., has the form
α = α 1 ( R ) g + α 2 ( R ) d + α 3 ( R ) e R e R .
The solenoidality condition (10) due to (33) and (34) leads to the following constraint on the scalar dislocation densities α k ( R ) ( k = 1 , 2 , 3 )
2 α 3 + R d α 3 d R = 2 α 1 .
We will seek the solution of the equilibrium Equation (2) and the incompatibility Equation (11) in the form of a spherically symmetric distortion tensor field
F = F 1 ( R ) g + F 2 ( R ) d + F 3 ( R ) e R e R .
Based on (32) and (36), the tensor Equation (11) is transformed to three scalar ordinary differential equations
2 F 2 = R α 3 , d d R R F 2 = R α 1 ,
F 3 d d R ( R F 1 ) = R α 2 .
Since the dislocation densities α k ( R ) , k = 1 , 2 , 3 , are assumed to be given functions, the first equation of (37) determines the function F 2 ( R ) , the second equation is a consequence of the first one and the relation (35). It remains to satisfy Equation (38).
From the constitutive Equations (3) and (13) and the expression (36) it follows that in the case of an isotropic material, the Piola stress tensor will be spherically symmetric, i.e., has the form
D = D 1 ( R ) g + D 2 ( R ) d + D 3 ( R ) e R e R .
Due to (32) and (39), the vector equilibrium Equation (2) is reduced to one scalar equation
d D 3 d R + 2 ( D 3 D 1 ) R = 0 .
The stresses D 1 and D 3 are expressed in terms of the distortion components F 1 , F 2 , and F 3 with the help of the constitutive relations of the material. Therefore, with the help of the relation (38), the Equation (40) is transformed to a second order nonlinear differential equation with an unknown function F 1 ( R ) . Thus, the problem of large deformations of an elastic sphere with a spherically symmetric distribution of dislocations is reduced to a nonlinear boundary value problem for an ordinary differential equation.
Some problems for elastic sphere with dislocations, the distribution of which is of spherical symmetry, are considered in the works [20,21,22,23,49]. The paper [20] describes the stress state of a hollow sphere made of a compressible semi-linear (harmonic) material containing edge dislocations and subject to external or internal hydrostatic pressures. Thanks to the selected material model and dislocation distribution, an exact solution was built. In [21], a universal solution for the class of isotropic incompressible elastic bodies was found, and eigenstresses in a solid sphere (in a space with a spherical cavity) with a special distribution of edge and screw dislocations, as well as stresses in the presence of hydrostatic pressure applied to the surface of the sphere (cavity), were determined. In [22], for compressible and incompressible materials, the stresses in a hollow sphere, caused by edge and screw dislocations, are numerically determined. The paper [49] contains an exact formulation and solution of the stability problem for a hollow sphere made of a semi-linear material with edge dislocations in the framework of the three-dimensional theory. In the work [23], in the framework of the model of a micropolar medium and in the classical theory of elasticity without couple stresses, an exact solution about the eigenstresses caused by the spherically symmetric distribution of dislocations and disclinations in a hollow linearly elastic sphere was found.
In this work, we investigate the equilibrium problem for an elastic sphere at large deformations taking into account distributed screw dislocations in the radial direction, i.e., in the case when α 1 = α 2 = 0 . In this case, the solenoidality requirement (35) completely defines the character of the dependence of the scalar screw dislocation density on the radial coordinate
α 3 ( R ) = 2 a R 2 ,
where a is a real constant. Denoting h ( R ) = R F 1 ( R ) , based on (36)–(38) and (41), we obtain
F = h ( R ) R g + a R d + d h ( R ) d R e R e R .
For comparison, consider the problem of inflation or compression of a hollow sphere in the absence of dislocations. In this case, the displacement field exists and deformation, as is known [24,26,32], has the form
r = r ( R ) , φ = Φ , θ = Θ ,
where φ , θ , and r are the spherical coordinates of the points of the sphere in the deformed state. The distortion tensor corresponding to the deformation (43) has the form
F = r ( R ) R g + d r ( R ) d R e R e R .
Comparing (42) and (44), we see that in the absence of dislocations, the function h ( R ) R represents the radial displacement of the points of the sphere.
Considering the material incompressible, from the incompressibility condition (15) and the expression (42), we obtain a differential equation for determining the function h ( R ) :
d h d R = R 2 a 2 + h 2 .
The solution to Equation (45) is written as follows:
h 3 h 0 3 + 3 a 2 ( h h 0 ) = R 3 R 0 3 , h 0 = h ( R 0 ) ,
where h 0 is an unknown constant. The cubic Equation (46) with respect to h ( R ) has only one real root
h ( R ) = ε + ε 2 + a 6 1 / 3 a 2 ε + ε 2 + a 6 1 / 3 , ε = 1 2 R 3 R 0 3 + 1 2 h 0 2 + 3 a 2 h 0 .
It can be shown [20] that in the absence of dislocations ( a = 0 ), negative values of the function h ( R ) = R 3 R 0 3 + h 0 3 3 correspond to the sphere eversion deformation given by the deformation tensor
F = d r d R e R e R r R e Φ e Φ r R e Θ e Θ .
Taking into account the incompressibility condition (45), we have the expression for the distortion tensor (42) in the form
F = h R g + a R d + R 2 h 2 + a 2 e R e R .
With the help of (48), the metric tensor (4) and its invariants (14) are calculated:
G = h 2 + a 2 R 2 g + R 4 h 2 + a 2 2 e R e R ,
I 1 = 2 h 2 + a 2 R 2 + R 4 h 2 + a 2 2 , I 2 = 2 R 2 h 2 + a 2 + h 2 + a 2 2 R 4 .
Using Formulas (28) and (29), we find the tensor F -T , the rotation tensor A , the stretch tensor U , and its invariants (18)
F -T = R h 2 + a 2 h g + a d + h 2 + a 2 R 2 e R e R ,
A = h h 2 + a 2 g + a h 2 + a 2 d + e R e R ,
U = h 2 + a 2 R g + R 2 h 2 + a 2 e R e R ,
J 1 = 2 h 2 + a 2 R + R 2 h 2 + a 2 , J 2 = 2 R h 2 + a 2 + h 2 + a 2 R 2 .
Although the displacement field does not exist in an elastic body with distributed dislocations, it is possible to determine the elongation of the material fibers and the change of the material surface area during deformation. Applying the well-known [47] formula for transforming an elementary oriented area under deformation, we calculate the area s of the surface into which the surface R = R 0 turns after deformation. Using (49), we obtain
s = R 0 2 e R · G -1 · e R 0 2 π d Φ π 2 π 2 cos Θ d Θ = 4 π h 0 2 + a 2 .
Formula (55) shows that for a 0 the value h 0 cannot be considered the outer radius of the deformed sphere.
Please note that distributed screw dislocations in the radial direction create a nonlinear effect in an elastic sphere. The influence of such dislocations cannot be detected within the framework of the linear theory of elasticity. Indeed, as is known [44,50], the incompatibility equation in the linear theory of elasticity has the form
Curl ( Curl e ) T = 1 2 Curl α T + 1 2 Curl α T T , e = Δ 1 2 F I + 1 2 F T I .
Here, e is the linear strain tensor. With the help of (32), one can check that if α = α 3 ( R ) e R e R , then the right-hand side of Equation (56) vanishes.
In addition to the incompatibility equation, the equilibrium Equation (40) and the boundary conditions must be fulfilled. We will assume that the outer surface of the sphere is loaded by a uniform hydrostatic pressure q 0 , and the inner one is loaded by a pressure q 1 . Then by virtue of (12), (39), and (51), we have the following boundary conditions
D 3 ( R 0 ) = h 0 2 + a 2 R 0 2 q 0 , D 3 ( R 1 ) = h 2 ( R 1 ) + a 2 R 1 2 q 1 .
Based on (16), (17), (39), and (51), the constitutive equations of an isotropic incompressible material in a spherically symmetric problem take the form
D 1 = D 1 * R h h 2 + a 2 p , D 2 = D 2 * a R h 2 + a 2 p , D 3 = D 3 * h 2 + a 2 R 2 p ,
D 1 * = h R τ 1 + I 1 τ 2 h 2 + a 2 R 2 τ 2 , D 2 * = a h D 1 * , D 3 * = R 2 h 2 + a 2 τ 1 + I 1 τ 2 R 4 h 2 + a 2 2 τ 2 .
If the specific energy is specified as a function of the tensor U , then, in accordance with (20), (52)–(54), the components of the tensor D * are written as
D 1 * = h η 1 + J 1 η 2 h 2 + a 2 η 2 R , D 2 * = a h D 1 * , D 3 * = η 1 + J 1 η 2 R 2 η 2 h 2 + a 2 .
The components of the tensor D * are known functions of the coordinate R, since they are expressed in terms of the function h ( R ) which is defined by Formula (47). The pressure function p ( R ) remains unknown. Using the relations (58), we transform the equilibrium Equation (40) into a differential equation for the function p ( R ) :
h 2 + a 2 R 2 d p d R = d D 3 * d R + 2 ( D 3 * D 1 * ) R .
To solve the boundary value problem (57) and (61), we introduce a new unknown function σ ( R ) :
σ ( R ) = R 2 h 2 + a 2 D 3 * p ( R ) .
With the help of the incompressibility condition (45), we find
d σ d R = 2 R D 1 * h 2 + a 2 2 h R 4 D 3 * h 2 + a 2 3 .
The boundary conditions (57) are written through the function σ ( R ) as follows:
σ ( R 1 ) = q 1 , σ ( R 0 ) = q 0 .
The solution of Equation (63), satisfying the first boundary condition in (64) has the form
σ ( R ) = 2 R 1 R R h 2 + a 2 D 1 * h R 3 h 2 + a 2 2 D 3 * d R q 1 .
The second boundary condition in (64) is used to determine the constant h 0 . The function p ( R ) and hence all stresses are now determined with the help of (62).
The solution (65) yields an analytical expression for the loading diagram of a sphere with dislocations, i.e., the dependence of the difference between the internal and external pressures on the deformation parameter h 0
q 1 q 0 = 2 R 1 R 0 R h 2 + a 2 D 1 * h R 3 h 2 + a 2 2 D 3 * d R .
The exact solution obtained here for the problem of large spherically symmetric deformations of an elastic hollow sphere with radial screw dislocations is universal in the class of isotropic incompressible materials.
Next, we will consider separately two cases of sphere loading: hydrostatic compression and inflation.
Let us write the solution of the differential Equation (61) in the form
p ( R ) = R 1 R f ( R ) d R + K ,
where K is the constant of integration determined from the boundary condition (57) on an unloaded surface.
For arbitrary incompressible material, we have
f ( R ) = R 2 h 2 + a 2 R 2 τ 1 h 2 + a 2 + 2 τ 2 + 2 R 3 h 3 a 2 h R 2 h 2 + a 2 + 2 τ 2 R 2 h h 2 + a 2 R 3 h R 3 h 2 + a 2 2 ,
and for Bartenev–Khazanovich material [24] (this model is also known as the Varga model), with η 1 = const = 2 and η 2 = 0 in (20), this function has the form
f ( R ) = 2 η 1 R h 2 + a 2 h h 2 + a 2 h 2 + a 2 2 .
In (68), ‘ ’ means the derivative of the function with respect to the argument R.

5. Numerical Results

For Section 5.2 and Section 5.3, the numerical results are given for models of two incompressible materials: neo-Hookean and Bartenev–Khazanovich. For each of them, the cases of thin ( R 1 = 0.95 R 0 ) and thick ( R 1 = 0.5 R 0 ) shells are considered.
Table 1 shows the pressure values for each specified case of the problem.

5.1. Problem for Solid Sphere

Since in the absence of dislocations the function h ( R ) defined by the expression (47) means a radial coordinate, then for a solid sphere the boundary condition
h ( 0 ) = 0
means that the sphere maintains the continuity property. Therefore, this condition is set in a natural way.
Suppose that the boundary condition (70) is also valid for a sphere with dislocations. Let us check this hypothesis in the case of hydrostatic compression of a hollow sphere by pressure q 0 , assuming the inner radius to be small. For specificity, we will assume R 1 = 10 7 . Please note that the radial coordinate is referred to the outer radius R 0 , which is equivalent to R 0 = 1 .
Consider a model of the neo-Hookean material [24], for which τ 1 = const = 1 (due to dimensionless stresses), τ 2 = 0 . Table 2 compares the values of h ( 0 ) for a solid sphere and h ( 10 7 ) for a hollow sphere for different dislocation densities. For the corresponding dislocation density, the parameter h 0 is determined from the problem of a solid sphere, for which the pressure q 0 is calculated from the problem of a hollow sphere. In this case, the more accurate the parameter h 0 is determined, the more the values of the functions h ( 0 ) and h ( 10 7 ) coincide.
Thus, the numerical analysis showed the validity of the boundary condition (70) for a solid sphere with dislocations, for which the solution taking into account this condition takes the form
h ( R ) = 1 2 R 3 + 1 2 R 6 + 4 a 6 1 / 3 a 2 1 2 R 3 + 1 2 R 6 + 4 a 6 1 / 3 .
Such a function is positive R due to the relation R 3 = h 3 + 3 a 2 h obtained from (45).

5.2. Inflation of Sphere

The inflation of the sphere is defined by the boundary condition (57) in the form
D 3 ( R 1 ) = q 1 F 1 2 ( R 1 ) + 1 4 α 3 2 ( R 1 ) R 1 2
which allows calculating h 0 in (46), and by the condition that the outer surface is not loaded
D 3 ( R 0 ) = 0 ,
with the help of which the constant K is found from (67).
For arbitrary incompressible material, the constant has the form
K = F 3 ( R 0 ) F 3 ( R 0 ) τ 1 + I 1 ( R 0 ) τ 2 τ 2 F 3 3 ( R 0 ) R 1 R 0 f ( R ) d R ,
and for Bartenev–Khazanovich material, we have
K = η 1 F 3 ( R 0 ) R 1 R 0 f ( R ) d R ,
where the function f ( R ) is determined by Formulas (68) and (69), respectively.
For inflation of the sphere made of the neo-Hookean material, the numerical results are presented in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 in the case of a thin shell and in Figure 6, Figure 7, Figure 8 and Figure 9 in the case of a thick shell. In these figures and below, T 3 and T 1 are the radial and circumferential Cauchy stresses, which are components of the symmetric Cauchy stress tensor T = ( det F ) 1 F T · D [24,31,32]. In this paper, the stresses are determined by the formulas: T 3 = F 3 D 3 , T 1 = F 1 D 1 + F 2 D 2 .
Since, in the absence of dislocations, h 0 is the outer radius of the sphere after deformation, and h 1 is the inner radius, the curves at a = 0 can be considered to be the dependences of the applied pressure on the outer and inner radii, respectively. Figure 1, Figure 2, Figure 6, Figure 7 show that with an increase in dislocation density, the stability loss occurs earlier. In addition, dislocations reduce the sphere’s resistance to inflation. In both figures, the graphs differ slightly.
As shown in Section 4, the values h 0 < 0 describe the eversion of the sphere. For h 0 > 0 , the right-hand branch of the curve with a maximum describes an unstable equilibrium position. To calculate the inflation stresses, we will use the left branch ( h 0 > 1 for a = 0 ) corresponding to a stable equilibrium position.
For a thin shell ( R 1 = 0.95 R 0 ) made of the neo-Hookean material, pressure q 1 = 0.02 corresponds to the following pairs ( a , h 0 ) used in Figure 3 (the scaled diagram is in Figure 4) and Figure 5: ( 0 , 1.034 ) , ( 0.5 , 0.911 ) , ( 0.9 , 0.585 ) , ( 1.1 , 0.428 ) .
For a thin shell ( R 1 = 0.95 R 0 ) made of the Bartenev–Khazanovich material, pressure q 1 = 0.02 corresponds to the following pairs ( a , h 0 ) used in Figure 10 (the scaled diagram is in Figure 11) and Figure 12: ( 0 , 1.035 ) , ( 0.5 , 0.91 ) , ( 0.9 , 0.579 ) , ( 1.1 , 0.41 ) .
For a thick shell ( R 1 = 0.5 R 0 ) made of the neo-Hookean material, pressure q 1 = 0.2 corresponds to the following pairs ( a , h 0 ) used in Figure 8 and Figure 9: ( 0 , 1.009 ) , ( 0.5 , 0.838 ) , ( 0.9 , 0.485 ) , ( 1.1 , 0.436 ) .
For a thick shell ( R 1 = 0.5 R 0 ) made of the Bartenev–Khazanovich material, pressure q 1 = 0.2 corresponds to the following pairs ( a , h 0 ) used in Figure 13 and Figure 14: ( 0 , 1.008 ) , ( 0.5 , 0.836 ) , ( 0.9 , 0.479 ) , ( 1.1 , 0.436 ) .
Dislocations do not significantly affect the radial stress in a thin shell (Figure 3 and Figure 4) and significantly increase this stress in a thick shell (Figure 8). In a thin shell, the circumferential stress (Figure 5) is distributed more uniformly than in a thick one (Figure 9). In both cases, dislocations increase the stresses in absolute value.
For the Bartenev–Khazanovich material, the solution to the sphere inflation problem is presented in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. The dependence of the applied pressure on the parameter h 0 , which is the outer radius after deformation in the absence of dislocations, is shown in Figure 15 (thin shell) and Figure 16 (thick shell). For h 0 > 0 , the points on the part of the curve with an extremum correspond to a stable equilibrium position up to the maximum point and an unstable one — after it. The dependence q 1 ( h 0 ) , as well as the stresses T 3 and T 1 for the Bartenev–Khazanovich material, are of the same nature as for the neo-Hookean material. For a thin shell, the radial stress is shown in Figure 10 and, to scale, in Figure 11, and the circumferential stress is in Figure 12. With this, the stresses in the Bartenev–Khazanovich material differ insignificantly from the stresses in the neo-Hookean material. Additionally, the stresses in a thick shell (Figure 13 and Figure 14) are close in magnitude to the corresponding stresses of the neo-Hookean material (Figure 8 and Figure 9).

5.3. Compression of Sphere by Hydrostatic Pressure

For hydrostatic compression of the sphere by the pressure q 0 , the boundary condition (57) is given by the relation
D 3 ( R 0 ) = q 0 F 1 2 ( R 0 ) + 1 4 α 3 2 ( R 0 ) R 0 2
and allows calculating the constant h 0 from (46). The no load condition on the inner surface
D 3 ( R 1 ) = 0
determines the unknown K from (67).
For an arbitrary incompressible material, the constant has the form
K = τ 1 + I 1 ( R 1 ) τ 2 F 3 2 ( R 1 ) τ 2 F 3 4 ( R 1 ) ,
and for Bartenev–Khazanovich material, we have
K = η 1 F 3 ( R 1 ) .
In the case when the hydrostatic pressure q 0 acts on the outer surface of the sphere R = R 0 , in Figure 19, Figure 20, Figure 21 and Figure 22 the stable branch for h 0 > 0 is the right one. For the convenience of distinguishing closely spaced curves in Figure 22, the scaled diagram is shown in Figure 23. The results obtained for the dependence q 0 ( h 0 ) agree with the tension diagram from [13].
For a thin shell, the radial stress (Figure 24) and the circumferential stress (Figure 25) in the case of a neo-Hookean material differ little from the corresponding stresses of the Bartenev–Khazanovich material (Figure 26 and Figure 27). The stresses in a thick shell (Figure 28, Figure 29, Figure 30 and Figure 31) show that dislocations increase the nonuniformity of dislocation distribution. In this case, the lower density of dislocations in some parts of the sphere causes higher stresses in absolute value.
For a thin shell ( R 1 = 0.95 R 0 ) made of the neo-Hookean material, pressure q 0 = 0.005 corresponds to the following pairs ( a , h 0 ) used in Figure 24 and Figure 25: ( 0 , 0.993 ) , ( 0.5 , 0.852 ) , ( 0.9 , 0.342 ) , ( 1.1 , 0.095 ) .
For a thin shell ( R 1 = 0.95 R 0 ) made of the Bartenev–Khazanovich material, pressure q 0 = 0.005 corresponds to the following pairs ( a , h 0 ) used in Figure 26 and Figure 27: ( 0 , 0.993 ) , ( 0.5 , 0.852 ) , ( 0.9 , 0.341 ) , ( 1.1 , 0.092 ) .
For a thick shell ( R 1 = 0.5 R 0 ) made of the neo-Hookean material, pressure q 0 = 0.2 corresponds to the following pairs ( a , h 0 ) used in Figure 28 and Figure 29: ( 0 , 0.994 ) , ( 0.5 , 0.793 ) , ( 0.9 , 0.052 ) , ( 1.1 , 0.155 ) .
For a thick shell ( R 1 = 0.5 R 0 ) made of the Bartenev–Khazanovich material, pressure q 0 = 0.2 corresponds to the following pairs ( a , h 0 ) used in Figure 30 and Figure 31: ( 0 , 0.994 ) , ( 0.1 , 0.987 ) , ( 0.2 , 0.965 ) , ( 0.5 , 0.792 ) .
As can be seen from the numerical analysis, the results for the neo-Hookean material and the Bartenev–Khazanovich material differ insignificantly. However, as seen in Figure 7 and Figure 18, the q 1 pressure values that correspond to the stable branch for the neo-Hookean material, for the Bartenev–Khazanovich material, may already correspond to the unstable branch.
Both these models are significant because they are widely used in research, experimentally confirmed and mathematically correct since they satisfy the condition of ellipticity at any deformation [24].

5.4. Eigenstresses in Sphere

The case of eigenstresses arising due to the presence of dislocations is investigated for a shell made of the neo-Hookean material (Figure 32 and Figure 33). In the absence of dislocations, the stresses are identically zero.
For a thin shell ( R 1 = 0.95 R 0 ) made of the neo-Hookean material, pressures q 0 = q 1 = 0 correspond to the following pairs ( a , h 0 ) used in Figure 32 and Figure 33: ( 0 , 1 ) , ( 0.5 , 0.862 ) , ( 0.9 , 0.401 ) , ( 1.1 , 0.021 ) .
Figure 34 shows the dependence on the dislocation parameter a 2 of the quantity h 0 for a thick shell R 1 = 0.5 R 0 in the absence of an external load.

6. Conclusions

In this work, we have solved the problem of large deformations of a hollow elastic sphere with distributed screw dislocations in the radial direction under loading by distributed pressure. The stress state of thin and thick shells was investigated for inflation and hydrostatic compression.
When solving the problem, we used the spherical symmetry of the distortion and stress tensors. The tensor equilibrium equation is reduced to a single scalar one which is transformed to a differential equation with respect to the pressure function in an incompressible body, not expressed in terms of deformation. The equation for one of the distortions, obtained from the incompressibility condition is reduced by introducing an auxiliary function to a cubic equation for this function. Together with the constitutive equations, the system of governing equations is closed by three scalar incompatibility equations, of which only one is independent.
The boundary conditions are the conditions of the presence of pressure on the inner surface of the sphere for inflation and on the outer surface — during hydrostatic compression. It is found that the boundary condition on the inner surface of a sphere of infinitesimal radius can be set with a small error in the same way as for a solid sphere. For numerical analysis, we used models of two incompressible materials: neo-Hookean and Bartenev–Khazanovich. A comparison of the stress state in the presence of dislocations and in their absence is carried out. It is observed that the dislocations have a nonlinear effect on the stresses.

Author Contributions

Conceptualization, L.M.Z.; methodology, L.M.Z.; software, E.V.G.; validation, L.M.Z.; formal analysis, L.M.Z. and E.V.G.; investigation, E.V.G.; data curation, E.V.G.; writing—original draft preparation, L.M.Z. and E.V.G.; writing—review and editing, L.M.Z.; visualization, E.V.G.; supervision, L.M.Z.; project administration, L.M.Z. Both authors have read and agreed to the published version of the manuscript.

Funding

The reported study was funded by the Russian Foundation of Basic Research, project number 19-31-90045.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Sedov, L.I. Propagation of strong shock waves. J. Appl. Math. Mech. 1946, 10, 241–250. [Google Scholar]
  2. SenGupta, S.; Sur, S. Spherically symmetric solutions of gravitational field equations in Kalb-Ramond background. Phys. Lett. B 2001, 521, 350–356. [Google Scholar] [CrossRef] [Green Version]
  3. Garcia-Diaz, A.A.; Gutierrez-Cano, G. Regularity conditions for spherically symmetric solutions of Einstein-nonlinear electrodynamics equations. Ann. Phys. 2020, 422, 168323. [Google Scholar] [CrossRef]
  4. Künzle, H.P.; Oliynyk, T.A. Spherically symmetric Einstein–Yang–Mills–Higgs fields for general compact gauge groups. Nonlinear Anal. 2005, 63, 473–480. [Google Scholar] [CrossRef]
  5. Baekler, P. Spherically symmetric solutions of the Poincaré gauge field theory. Phys. Lett. 1983, 96A, 279–282. [Google Scholar] [CrossRef]
  6. Dymnikova, I.G. The algebraic structure of a cosmological term in spherically symmetric solutions. Phys. Lett. B 2000, 472, 33–38. [Google Scholar] [CrossRef] [Green Version]
  7. Wheeler, J.T. Symmetric solutions to the maximally Gauss – Bonnet extended Einstein equations. Nuclear Phys. 1986, B273, 732–748. [Google Scholar] [CrossRef] [Green Version]
  8. Kondo, K. On the geometrical and physical foundations in the theory of yielding. In Proceedings of the 2nd Japan National Congress of Applied Mechanics, Istanbul, Turkey, 20–28 August 1952; pp. 41–47. [Google Scholar]
  9. Bilby, B.A.; Bullough, R.; Smith, E. Continuous distributions of dislocations: A new application of the methods of non-Riemannian geometry. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 1955, A231, 263–273. [Google Scholar]
  10. Kröner, E. Allgemeine Kontinuumstheorie der Versetzungen und Eigenspannungen. Arch. Ration. Mech. Anal. 1960, 4, 273–334. [Google Scholar] [CrossRef]
  11. Zorski, H. Statistical theory of dislocations. Int. J. Solids Struct. 1968, 4, 959–974. [Google Scholar] [CrossRef]
  12. Mossakowska, Z. Displacement description of dislocation lines. I. Cyclic functions. Arch. Mech. 1979, 31, 533–545. [Google Scholar]
  13. Zubov, L.M. Nonlinear Theory of Dislocations and Disclinations in Elastic Bodies; Springer: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  14. Derezin, S.V.; Zubov, L.M. Equations of a nonlinear elastic medium with continuously distributed dislocations and disclinations. Dokl. Phys. 1999, 44, 391–394. [Google Scholar]
  15. Derezin, S.V.; Zubov, L.M. Disclinations in nonlinear elasticity. Ztsch. Angew. Math. Mech. 2011, 91, 433–442. [Google Scholar] [CrossRef]
  16. Yavary, A.; Goriely, A. Riemann—Cartan geometry of nonlinear dislocation mechanics. Arch. Ration. Mech. Anal. 2012, 205, 59–118. [Google Scholar] [CrossRef]
  17. Zelenina, A.A.; Zubov, L.M. Bending and twisting of nonlinear elastic bodies with continuously distributed dislocations. Vestn. Yuzhn. Nauchn. Tsentr. RAN 2009, 3, 15–22. [Google Scholar]
  18. Zubov, L.M. Spherically symmetric solutions in the nonlinear theory of dislocations. Dokl. Phys. 2014, 59, 419–422. [Google Scholar] [CrossRef]
  19. Zelenina, A.A.; Zubov, L.M. Nonlinear effects during the tension, bend, and torsion of elastic bodies with distributed dislocations. Dokl. Phys. 2013, 58, 354–357. [Google Scholar] [CrossRef]
  20. Zhbanova, E.V.; Zubov, L.M. The Influence of Distributed Dislocations on Large Deformations of an Elastic Sphere. In Advanced Methods of Continuum Mechanics for Materials and Structures; Advanced Structured Materials; Naumenko, K., Aßmus, M., Eds.; Springer: Singapore, 2016; Volume 60, pp. 61–76. [Google Scholar]
  21. Goloveshkina, E.V.; Zubov, L.M. Universal spherically symmetric solution of nonlinear dislocation theory for incompressible isotropic elastic medium. Arch. Appl. Mech. 2019, 89, 409–424. [Google Scholar] [CrossRef]
  22. Goloveshkina, E.V.; Zubov, L.M. Eigenstresses in a Nonlinearly Elastic Sphere with Distributed Dislocations. In New Achievements in Continuum Mechanics and Thermodynamics; Advanced Structured Materials; Abali, B., Altenbach, H., dell’Isola, F., Eremeyev, V., Öchsner, A., Eds.; Springer International Publishing: Cham, Switzerland, 2019; Volume 108, pp. 137–155. [Google Scholar]
  23. Zelenina, A.A.; Zubov, L.M. Spherically Symmetric Deformations of Micropolar Elastic Medium with Distributed Dislocations and Disclinations. In Advances in Mechanics of Microstructured Media and Structures; dell’Isola, F., Eremeyev, V.A., Porubov, A., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 357–369. [Google Scholar]
  24. Lurie, A.I. Nonlinear Theory of Elasticity; North-Holland: Amsterdam, The Netherlands, 1990. [Google Scholar]
  25. Rivlin, R. Large Elastic Deformations of Isotropic Materials. IV. Further Developments of the General Theory. Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci. 1948, 241, 379–397. [Google Scholar]
  26. Truesdell, C.; Noll, W. The Non-Linear Field Theories of Mechanics. In Handbuch der Physik; Flügge, S., Ed.; Springer: Berlin, Germany, 1965; Volume 3, pp. 1–602. [Google Scholar]
  27. Rivlin, R. Large elastic deformations of isotropic materials. V. The problem of flexure. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 1949, 195, 463–473. [Google Scholar]
  28. Rivlin, R. Large elastic deformations of isotropic materials. VI. Further results in the theory of torsion, shear and flexure. Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci. 1949, 242, 173–195. [Google Scholar]
  29. Ericksen, J.L. Deformations possible in every isotropic, incompressible, perfectly elastic body. Z. Angew. Math. Phys. ZAMP 1954, 5, 466–489. [Google Scholar] [CrossRef]
  30. Klingbeil, W.W.; Shield, R.T. On a class of solutions in plane finite elasticity. Z. Angew. Math. Phys. ZAMP 1966, 17, 489–511. [Google Scholar] [CrossRef]
  31. Truesdell, C. A First Course in Rational Continuum Mechanics; Academic Press: New York, NY, USA, 1977. [Google Scholar]
  32. Ogden, R.W. Non-Linear Elastic Deformations; Dover: New York, NY, USA, 1997. [Google Scholar]
  33. Pucci, E.; Saccomandi, G. Universal motions for constrained simple materials. Int. J. Non-Linear Mech. 1998, 34, 469–484. [Google Scholar] [CrossRef]
  34. Saccomandi, G. Universal Results in Finite Elasticity. Non-linear Elasticity: Theory and Applications. Lond. Math. Soc. Lect. Note Ser. 2001, 283, 97–134. [Google Scholar]
  35. Saccomandi, G.; Beatty, M.F. Universal Relations for Fiber-Reinforced Elastic Materials. Math. Mech. Solids 2002, 7, 95–110. [Google Scholar] [CrossRef]
  36. Eremeyev, V.; Lebedev, L.P.; Altenbach, H. Foundations of Micropolar Mechanics; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  37. Lurie, K.A. The problem of effective parameters of a mixture of two isotropic dielectrics distributed in space-time and the conservation law for wave impedance in one-dimensional wave propagation. Proc. R. Soc. A Math., Phys. Eng. Sci. 1998, A454, 1767–1779. [Google Scholar] [CrossRef]
  38. Lurie, K.A. An Introduction to the Mathematical Theory of Dynamic Materials. In Advances in Mechanics and Mathematics; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2007; Volume 15. [Google Scholar]
  39. Maugin, G.A. Continuum Mechanics Through the Twentieth Century; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  40. Dorfmann, A.; Ogden, R.W.; Saccomandi, G. Universal relations for non-linear magnetoelastic solids. Int. J. Non-Linear Mech. 2004, 39, 1699–1708. [Google Scholar] [CrossRef]
  41. Zubov, L.M. Universal deformations of micropolar isotropic elastic solids. Math. Mech. Solids 2016, 21, 152–167. [Google Scholar] [CrossRef]
  42. Zubov, L.M. Universal Solutions of Nonlinear Dislocation Theory for Elastic Cylinder. Mech. Solids 2020, 55, 701–709. [Google Scholar] [CrossRef]
  43. Nye, J.F. Some geometrical relations in dislocated crystals. Acta Met. 1953, 1, 153–162. [Google Scholar] [CrossRef]
  44. Eshelby, J.D. The continuum theory of Lattice Defects. In Solid State Physics; Seitz, F., Turnbul, D., Eds.; Academic Press: New York, NY, USA, 1956; Volume 3, pp. 79–144. [Google Scholar]
  45. Landau, L.D.; Lifshitz, E.M. Theory of Elasticity; Theoretical Physics; Pergamon: Oxford, UK, 1975; Volume 7. [Google Scholar]
  46. Vakulenko, A.A. The relationship of micro- and macroproperties in elastic-plastic media. Itogi Nauk. Tekh. Ser. Mekh. Deform. Tverd. Tela 1991, 22, 3–54. [Google Scholar]
  47. Zubov, L.M.; Karyakin, M.I. Tensor Calculus; Vuzovskaya kniga: Moscow, Russia, 2006. (In Russian) [Google Scholar]
  48. Eremeyev, V.A.; Cloud, M.J.; Lebedev, L.P. Applications of Tensor Analysis in Continuum Mechanics; World Scientific: Hackensack, NJ, USA, 2018. [Google Scholar]
  49. Goloveshkina, E.V.; Zubov, L.M. Equilibrium stability of nonlinear elastic sphere with distributed dislocations. Contin. Mech. Therm. 2020, 32, 1713–1725. [Google Scholar] [CrossRef]
  50. Teodosiu, C. Elastic Models of Crystal Defects; Springer: Berlin/Heidelberg, Germany, 1982. [Google Scholar]
Figure 1. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 1. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g001
Figure 2. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 2. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g002
Figure 3. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , radial Koshi stress T 3 .
Figure 3. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , radial Koshi stress T 3 .
Symmetry 13 00830 g003
Figure 4. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , scaled radial Koshi stress T 3 ; R 0 = 1 , R 1 = 0.95 .
Figure 4. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , scaled radial Koshi stress T 3 ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g004
Figure 5. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , circumferential Koshi stress T 1 .
Figure 5. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g005
Figure 6. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 6. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g006
Figure 7. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 7. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g007
Figure 8. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , radial Koshi stress T 3 .
Figure 8. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , radial Koshi stress T 3 .
Symmetry 13 00830 g008
Figure 9. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , circumferential Koshi stress T 1 .
Figure 9. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g009
Figure 10. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , radial Koshi stress T 3 .
Figure 10. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , radial Koshi stress T 3 .
Symmetry 13 00830 g010
Figure 11. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , scaled radial Koshi stress T 3 .
Figure 11. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , scaled radial Koshi stress T 3 .
Symmetry 13 00830 g011
Figure 12. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , circumferential Koshi stress T 1 .
Figure 12. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 = 0.02 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g012
Figure 13. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , radial Koshi stress T 3 .
Figure 13. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , radial Koshi stress T 3 .
Symmetry 13 00830 g013
Figure 14. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , circumferential Koshi stress T 1 .
Figure 14. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 = 0.2 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g014
Figure 15. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 15. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g015
Figure 16. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 16. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g016
Figure 17. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 17. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thin shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g017
Figure 18. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 18. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), inflation of thick shell by pressure q 1 , dependence q 1 ( h 1 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g018
Figure 19. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 19. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g019
Figure 20. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 20. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g020
Figure 21. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Figure 21. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.95 .
Symmetry 13 00830 g021
Figure 22. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 22. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell, dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g022
Figure 23. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell, scaled dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Figure 23. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell, scaled dependence q 0 ( h 0 ) ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g023
Figure 24. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , radial Koshi stress T 3 .
Figure 24. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , radial Koshi stress T 3 .
Symmetry 13 00830 g024
Figure 25. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell by pressure q 1 = 0.005 , circumferential Koshi stress T 1 .
Figure 25. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thin shell by pressure q 1 = 0.005 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g025
Figure 26. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , radial Koshi stress T 3 .
Figure 26. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , radial Koshi stress T 3 .
Symmetry 13 00830 g026
Figure 27. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , circumferential Koshi stress T 1 .
Figure 27. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thin shell by pressure q 0 = 0.005 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g027
Figure 28. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , radial Koshi stress T 3 .
Figure 28. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , radial Koshi stress T 3 .
Symmetry 13 00830 g028
Figure 29. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , circumferential Koshi stress T 1 .
Figure 29. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g029
Figure 30. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , radial Koshi stress T 3 .
Figure 30. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , radial Koshi stress T 3 .
Symmetry 13 00830 g030
Figure 31. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , circumferential Koshi stress T 1 .
Figure 31. Bartenev–Khazanovich material ( η 1 = 2 , η 2 = 0 ), hydrostatic compression of thick shell by pressure q 0 = 0.2 , circumferential Koshi stress T 1 .
Symmetry 13 00830 g031
Figure 32. Thin shell made of neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), radial Koshi eigenstress T 3 .
Figure 32. Thin shell made of neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), radial Koshi eigenstress T 3 .
Symmetry 13 00830 g032
Figure 33. Thin shell made of neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), circumferential Koshi eigenstress T 1 .
Figure 33. Thin shell made of neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), circumferential Koshi eigenstress T 1 .
Symmetry 13 00830 g033
Figure 34. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), eigenstresses in sphere, dependence of quantity h 0 on dislocation parameter a 2 ; R 0 = 1 , R 1 = 0.5 .
Figure 34. Neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ), eigenstresses in sphere, dependence of quantity h 0 on dislocation parameter a 2 ; R 0 = 1 , R 1 = 0.5 .
Symmetry 13 00830 g034
Table 1. Considered cases of applied loads for material models: external pressure q 0 and internal pressure q 1 .
Table 1. Considered cases of applied loads for material models: external pressure q 0 and internal pressure q 1 .
MaterialThin ShellThick Shell
neo-Hookean
Bartenev–Khazanovich
q 1 = 0.02 , q 0 = 0.005 q 1 = 0.2 , q 0 = 0.2
Table 2. Comparison of the values of h ( 0 ) in the case of a solid sphere and h ( 10 7 ) in the case of a hollow sphere for neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ).
Table 2. Comparison of the values of h ( 0 ) in the case of a solid sphere and h ( 10 7 ) in the case of a hollow sphere for neo-Hookean material ( τ 1 = 1 , τ 2 = 0 ).
a h 0 q 0 for a Hollow Sphere h ( 0 ) h ( 10 7 )
01000
0.5 0.7564 0.582 0 0.00009
1 0.3221 0.107 0 0.00009
1.5 0.1477 0.032 0 0.0003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Goloveshkina, E.V.; Zubov, L.M. Spherically Symmetric Tensor Fields and Their Application in Nonlinear Theory of Dislocations. Symmetry 2021, 13, 830. https://doi.org/10.3390/sym13050830

AMA Style

Goloveshkina EV, Zubov LM. Spherically Symmetric Tensor Fields and Their Application in Nonlinear Theory of Dislocations. Symmetry. 2021; 13(5):830. https://doi.org/10.3390/sym13050830

Chicago/Turabian Style

Goloveshkina, Evgeniya V., and Leonid M. Zubov. 2021. "Spherically Symmetric Tensor Fields and Their Application in Nonlinear Theory of Dislocations" Symmetry 13, no. 5: 830. https://doi.org/10.3390/sym13050830

APA Style

Goloveshkina, E. V., & Zubov, L. M. (2021). Spherically Symmetric Tensor Fields and Their Application in Nonlinear Theory of Dislocations. Symmetry, 13(5), 830. https://doi.org/10.3390/sym13050830

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop