Next Article in Journal
A Novel Method of Calibrating Micro-Scale Parameters of PFC Model and Experimental Validation
Previous Article in Journal
Autonomous Navigation Framework for Intelligent Robots Based on a Semantic Environment Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solution Behavior in the Vicinity of Characteristic Envelopes for the Double Slip and Rotation Model

1
National Key Laboratory for Precision Hot Processing of Metals, Harbin Institute of Technology, Harbin 150001, China
2
School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, China
3
Ishlinsky Institute for Problems in Mechanics RAS, Moscow 119526, Russia
4
Federal State Autonomous Educational Institution of Higher Education, South Ural State University (national research university), Chelyabinsk 454080, Russia
5
Division of Computational Mathematics and Engineering, Institute for Computational Science, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
6
Faculty of Civil Engineering, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(9), 3220; https://doi.org/10.3390/app10093220
Submission received: 28 January 2020 / Revised: 11 March 2020 / Accepted: 20 March 2020 / Published: 6 May 2020
(This article belongs to the Section Mechanical Engineering)

Abstract

:
The boundary conditions significantly affect solution behavior near rough interfaces. This paper presents general asymptotic analysis of solutions for the rigid plastic double slip and rotation model in the vicinity of an envelope of characteristics under plane strain and axially symmetric conditions. This model is used in the mechanics of granular materials. The analysis has important implications for solving boundary value problems because the envelope of characteristics is a natural boundary of the analytic solution. Moreover, an envelope of characteristics often coincides with frictional interfaces. In this case, the regime of sticking is not possible independently of the friction law chosen. It is shown that the solution is singular in the vicinity of envelopes. In particular, the profile of the velocity component tangential to the envelope is described by the sum of the constant and square root functions of the normal distance to the envelope in its vicinity. As a result, some components of the strain rate tensor approach infinity. This finding might help to develop an efficient numerical method for solving boundary value problems and provide the basis for the interpretation of some experimental results.

1. Introduction

Models of pressure-dependent plasticity based on the Mohr–Coulomb yield criterion are widely used for describing deformation of granular materials and soil [1,2]. Boundary conditions for the flow of granular materials are attracting considerable interest due to their effect on solution behavior [3,4,5,6,7,8,9,10]. In the case of hyperbolic models, the envelope of characteristics is a natural boundary of the analytic solution. The envelope of characteristics often coincides with frictional interfaces. A distinguished feature of this boundary condition is that the regime of sliding occurs independently of the friction law chosen. In the case of pressure-dependent plasticity, such solutions have been found in [11,12,13]. In these works, the double shearing model proposed in [14] has been adopted. The solutions given in [11,12,13] show that the velocity field is singular (some derivatives approach infinity in the vicinity of envelopes).
The present paper deals with the rigid plastic double slip and rotation model under plane strain and axial symmetry conditions [15]. The solutions found in [11,12] are also solutions for this model if the intrinsic spin involved in the double slip and rotation model vanishes. Other singular solutions for the rigid plastic double slip and rotation model have been derived in [16,17].
General asymptotic analyses of the singular behavior of solutions help to develop efficient numerical methods for solving boundary value problems and provide the basis for the interpretation of some experimental results. In the case of rigid plastic solids, such analyses have been carried out in [18] for isotropic perfectly plastic material, [19] for anisotropic perfectly plastic material, and [20] for a model that is used in the mechanics of polymers. These findings show that the exact asymptotic representation of solutions near envelopes of characteristics is affected by the shape of the yield surface and other constitutive equations. This paper supplies the exact asymptotic representation of solutions near envelopes of characteristics for the double slip and rotation model.
The applied aspect of the present paper is twofold. Firstly, the behavior of soil–structure interfaces is very important for predicting the overall response of soil–structure systems. It is known that a very narrow layer of soil adjacent to the structure controls the shearing deformation near the interface. Several models have been proposed for describing the behavior of material within this layer. One group of models introduces special constitutive equations in the vicinity of soil–structure interfaces (for example, [21,22,23]). The other group treats the interface layer as a geometric surface subject to certain physical conditions (for example, [24,25,26]). The exact asymptotic representation of solutions found in the present paper can be used in conjunction with the approach adopted in models of the first group. Secondly, suitable interface elements should be used to accurately simulate the behavior of material near soil–structure interfaces via the finite element method [27,28]. The present paper demonstrates that the solution may be singular. Therefore, traditional linear approximations with finite elements, for example [28], are not capable of representing the real solution. The exact asymptotic representation of solutions found can be used in conjunction with the extended finite element method [29].

2. Basic Equations

The double slip and rotation model is a model of pressure-dependent plasticity [15]. The version of the model considered in the present paper is rigid/plastic (i.e., the elastic portion of strain is neglected). The constitutive equations of the model are a pressure-dependent yield criterion and a flow rule. The present paper deals with the Mohr–Coulomb yield criterion. The flow rule has been formulated in [15].

2.1. Plane Strain Deformation

Let ( x 1 , x 2 ) be an orthogonal coordinate system in planes of flow. The x 2 coordinate curves are straight (Figure 1). Then, it is possible to assume with no loss of generality that the scale factor of these coordinate curves is unity. The scale factor of the x 1 coordinate lines may be represented as
H 1 = 1 + x 2 R ( x 1 ) .
Here R ( x 1 ) is the radius of curvature of the x 1 coordinate curve determined by the equation x 2 = 0 (S–line in Figure 1). The equilibrium equations referred to the coordinate system chosen are [30]
σ 11 x 1 + [ 1 + x 2 R ( x 1 ) ] σ 12 x 2 + 2 σ 12 R ( x 1 ) = 0 , [ 1 + x 2 R ( x 1 ) ] σ 22 x 2 + σ 12 x 1 + ( σ 22 σ 11 ) R ( x 1 ) = 0 .
Here, Equation (1) has been used. Furthermore, σ 11 , σ 22 and σ 12 are the stress components referred to the ( x 1 , x 2 ) coordinate system.
As usual, it is convenient to introduce the following stress variables:
p = 1 2 ( σ 11 + σ 22 ) , q = 1 2 ( σ 11 σ 22 ) 2 + 4 σ 12 2 > 0 .
Then, the Mohr–Coulomb yield criterion reads
q p sin ϕ = k cos ϕ .
Here k is the cohesion and ϕ is the angle of internal friction. Both are material constants. Let ψ be the orientation of the major principal stress σ 1 relative to the x 1 direction (Figure 2). Then, using (3), the stress components referred to the ( x 1 , x 2 ) coordinate system are expressed as
σ 11 = p + q cos 2 ψ , σ 22 = p q cos 2 ψ , σ 12 = q sin 2 ψ .
Substituting (5) into (2) leads to
p x 1 + cos 2 ψ q x 1 2 q sin 2 ψ ψ x 1 + [ 1 + x 2 R ( x 1 ) ] sin 2 ψ q x 2 + 2 q cos 2 ψ [ 1 + x 2 R ( x 1 ) ] ψ x 2 + 2 q sin 2 ψ R ( x 1 ) = 0 , [ 1 + x 2 R ( x 1 ) ] ( p x 2 + cos 2 ψ q x 2 ) + 2 q [ 1 + x 2 R ( x 1 ) ] sin 2 ψ ψ x 2 + sin 2 ψ q x 1 + 2 q cos 2 ψ ψ x 1 2 q cos 2 ψ R ( x 1 ) = 0 .
One can eliminate p in these equations using (4). As a result,
( cos 2 ψ 1 sin ϕ ) q x 1 2 q sin 2 ψ ψ x 1 + [ 1 + x 2 R ( x 1 ) ] sin 2 ψ q x 2 + 2 q cos 2 ψ [ 1 + x 2 R ( x 1 ) ] ψ x 2 + 2 q sin 2 ψ R ( x 1 ) = 0 , [ 1 + x 2 R ( x 1 ) ] ( cos 2 ψ + 1 sin ϕ ) q x 2 + 2 q [ 1 + x 2 R ( x 1 ) ] sin 2 ψ ψ x 2 + sin 2 ψ q x 1 + 2 q cos 2 ψ ψ x 1 2 q cos 2 ψ R ( x 1 ) = 0 .
The flow rule is [15]
ξ 11 + ξ 22 = 0 , ξ 12 cos 2 ψ 1 2 ( ξ 11 ξ 22 ) sin 2 ψ + ( ω 12 + Ω ) sin ϕ = 0 .
Here ξ 11 , ξ 22 and ξ 12 are the strain rate components referred to the ( x 1 , x 2 ) coordinate system. Furthermore, ω 12 is the only non-zero spin component and Ω is the intrinsic spin due to grain rotation. In general, the latter is an unknown variable which is governed by the equation of rotational motion. For the objective of the present paper it is sufficient to assume that Ω may be represented by a Taylor series in the vicinity of x 2 = 0 . Using the first equation in (8) the second equation can be rewritten as
ξ 12 cos 2 ψ ξ 11 sin 2 ψ + ( ω 12 + Ω ) sin ϕ = 0 .
The strain rate and spin components are expressed as
ξ 11 = ( u 1 x 1 + u 2 R ( x 1 ) ) [ 1 + x 2 R ( x 1 ) ] 1 , ξ 22 = u 2 x 2 , 2 ξ 12 = u 1 x 2 + ( u 2 x 1 u 1 R ( x 1 ) ) [ 1 + x 2 R ( x 1 ) ] 1 , 2 ω 12 = u 1 x 2 ( u 2 x 1 u 1 R ( x 1 ) ) [ 1 + x 2 R ( x 1 ) ] 1 .
Here u 1 and u 2 are the velocity components referred to the ( x 1 , x 2 ) coordinate system. Substituting (10) into (8) and (9) yields
u 1 x 1 + [ 1 + x 2 R ( x 1 ) ] u 2 x 2 = u 2 R ( x 1 ) , ( cos 2 ψ sin ϕ 1 ) u 2 x 1 + ( cos 2 ψ sin ϕ + 1 ) [ 1 + x 2 R ( x 1 ) ] u 1 x 2 2 sin 2 ψ sin ϕ u 1 x 1 = 2 u 2 sin 2 ψ R ( x 1 ) sin ϕ ( sin ϕ cos 2 ψ ) u 1 R ( x 1 ) sin ϕ 2 Ω [ 1 + x 2 R ( x 1 ) ] = 0 .

2.2. Axisymmetric Deformation

Let ( x 1 , x 2 ) be an orthogonal coordinate system in the generic plane of a cylindrical coordinate system ( r , θ , z ) and γ is the angle between the r-axis and the x 1 direction measured from the r-axis anticlockwise (Figure 3). As in the case of plane strain deformation, the x 2 lines are straight. Then, the scale factor of these coordinate lines are unity and the scale factor of the x 1 coordinate lines is given by (1). The solution is independent of θ . Moreover, the circumferential velocity and the shear stresses σ 1 θ and σ 2 θ vanish everywhere. The non-trivial equilibrium equations are [30]
σ 11 x 1 + [ 1 + x 2 R ( x 1 ) ] σ 12 x 2 + ( σ 11 σ θ θ ) r r x 1 + σ 12 [ ( 1 + x 2 R ( x 1 ) ) r r x 2 + 2 R ( x 1 ) ] = 0 , σ 22 x 2 + [ 1 + x 2 R ( x 1 ) ] 1 σ 12 x 1 + ( σ 22 σ θ θ ) r r x 2 + ( σ 22 σ 11 ) R ( x 1 ) + x 2 + ( 1 + x 2 R ( x 1 ) ) 1 σ 12 r r x 1 = 0 .
It follows from the geometry of Figure 3 that r / x 1 = ( 1 + x 2 / R ( x 1 ) ) cos γ and r / x 2 = sin γ . Then, the equations in (12) become
σ 11 x 1 + [ 1 + x 2 R ( x 1 ) ] σ 12 x 2 + ( σ 11 σ θ θ ) r [ 1 + x 2 R ( x 1 ) ] cos γ + σ 12 [ 2 R ( x 1 ) ( 1 + x 2 R ( x 1 ) ) sin γ r ] = 0 , σ 22 x 2 + [ 1 + x 2 R ( x 1 ) ] 1 σ 12 x 1 ( σ 22 σ θ θ ) sin γ r + ( σ 22 σ 11 ) R ( x 1 ) + x 2 + σ 12 cos γ r = 0 .
There are several flow regimes associated with difference faces and edges of the Mohr–Coulomb yield criterion (Figure 4). One of the principal stresses is the circumferential stress σ θ θ . The other principal stresses are denoted as σ 1 and σ 2 . It is assumed that σ 1 > σ 2 . In this case, there are two regimes in which the hypothesis of Haar and von Karmen, that σ 1 = σ θ θ or σ 2 = σ θ θ , is satisfied. These regimes correspond to point A and F in Figure 4. Using (3) the yield criterion can be represented as
σ θ θ = p m q , q = p sin ϕ + k cos ϕ
where m = + 1 for regime A, and m = 1 for regime F. The equations in (5) are valid. Substituting (5) and the first equation in (14) into (13) yields
p x 1 + cos 2 ψ q x 1 2 q sin 2 ψ ψ x 1 + [ 1 + x 2 R ( x 1 ) ] ( sin 2 ψ q x 2 + 2 q cos 2 ψ ψ x 2 ) + [ 1 + x 2 R ( x 1 ) ] ( cos 2 ψ + m ) q cos γ r + [ 2 R ( x 1 ) ( 1 + x 2 R ( x 1 ) ) sin γ r ] q sin 2 ψ = 0 , p x 2 cos 2 ψ q x 2 + 2 q sin 2 ψ ψ x 2 + [ 1 + x 2 R ( x 1 ) ] 1 ( sin 2 ψ q x 1 + 2 q cos 2 ψ ψ x 1 ) + ( cos 2 ψ m ) q sin γ r 2 q cos 2 ψ R ( x 1 ) + x 2 + q sin 2 ψ cos γ r = 0 .
One can eliminate p here using the second equation in (14). Then, the equations in (15) become
( cos 2 ψ 1 sin ϕ ) q x 1 2 q sin 2 ψ ψ x 1 + [ 1 + x 2 R ( x 1 ) ] ( sin 2 ψ q x 2 + 2 q cos 2 ψ ψ x 2 ) + [ 1 + x 2 R ( x 1 ) ] ( cos 2 ψ + m ) q cos γ r + [ 2 R ( x 1 ) ( 1 + x 2 R ( x 1 ) ) sin γ r ] q sin 2 ψ = 0 , ( cos 2 ψ + 1 sin ϕ ) q x 2 + 2 q sin 2 ψ ψ x 2 + [ 1 + x 2 R ( x 1 ) ] 1 ( sin 2 ψ q x 1 + 2 q cos 2 ψ ψ x 1 ) + ( cos 2 ψ m ) q sin γ r 2 q cos 2 ψ R ( x 1 ) + x 2 + q sin 2 ψ cos γ r = 0 .
The flow rule is
ξ 11 + ξ 22 + ξ θ θ = 0 , ξ 12 cos 2 ψ 1 2 ( ξ 11 ξ 22 ) sin 2 ψ + ( ω 12 + Ω ) sin ϕ = 0 .
Here ξ θ θ is the circumferential strain rate. The equations in (10) are valid. Moreover, ξ θ θ = u r / r where u r is the radial velocity. It follows from the geometry of Figure 3 that u r = u 1 cos γ u 2 sin γ . Using this equation and (10) it is possible to transform (17) into
u 1 x 1 + [ 1 + x 2 R ( x 1 ) ] u 2 x 2 = ( u 2 sin γ u 1 cos γ ) r [ 1 + x 2 R ( x 1 ) ] u 2 R ( x 1 ) , ( cos 2 ψ sin ϕ 1 ) u 2 x 1 + ( cos 2 ψ sin ϕ + 1 ) [ 1 + x 2 R ( x 1 ) ] u 1 x 2 s i n 2 ψ sin ϕ u 1 x 1 + [ 1 + x 2 R ( x 1 ) ] sin 2 ψ sin ϕ u 2 x 2 = ( cos 2 ψ sin ϕ ) u 1 R ( x 1 ) sin ϕ + sin 2 ψ u 2 R ( x 1 ) sin ϕ 2 Ω [ 1 + x 2 R ( x 1 ) ] .

3. Characteristic Analysis

3.1. Plane Strain Deformation

Multiplying the first equation in (7) by cos 2 ψ , the second by sin 2 ψ and summing gives
( 1 cos 2 ψ sin ϕ ) q x 1 sin 2 ψ sin ϕ [ 1 + x 2 R ( x 1 ) ] q x 2 + 2 q [ 1 + x 2 R ( x 1 ) ] ψ x 2 = 0 .
Multiplying the first equation in (7) by sin 2 ψ , the second by cos 2 ψ and summing gives
sin 2 ψ sin ϕ q x 1 ( 1 + cos 2 ψ sin ϕ ) [ 1 + x 2 R ( x 1 ) ] q x 2 + 2 q ψ x 1 = 2 q R ( x 1 ) .
Assume that S – curve in Figure 1 is a characteristic curve. Then, ψ = ± ( π / 4 + ϕ / 2 ) on this curve [14]. In particular, ψ / x 1 = 0 . Equation (20) supplies the characteristic relation in the form
± cot ϕ q x 1 2 q R ( x 1 ) = 0 .
In the proceeding, it is assumed that S – curve in Figure 1 is an envelope of characteristics. Therefore, Equation (21) is not satisfied even though ψ = ± ( π / 4 + ϕ / 2 ) on this curve. It is possible only if
| q x 2 |
as x 2 0 . In this case, the second term in (20) reduces to the expression 0 at x 2 = 0 .

3.2. Axisymmetric Deformation

Multiplying the first equation in (16) by cos ψ , the second by sin ψ and summing gives
( 1 cos 2 ψ sin ϕ ) q x 1 sin 2 ψ sin ϕ [ 1 + x 2 R ( x 1 ) ] q x 2 + 2 q [ 1 + x 2 R ( x 1 ) ] ψ x 2 = F 1 .
Multiplying the first equation in (16) by sin ψ , the second by cos ψ and summing gives
sin 2 ψ sin ϕ q x 1 ( 1 + cos 2 ψ sin ϕ ) [ 1 + x 2 R ( x 1 ) ] q x 2 + 2 q ψ x 1 = F 2 .
The terms F 1 and F 2 include no derivatives.
Assume that S – curve in Figure 3 is a characteristic curve. Then, ψ = ± ( π / 4 + ϕ / 2 ) on this curve [31]. In particular, ψ / x 1 = 0 . Equation (24) supplies the characteristic relation in the form
± cot ϕ q x 1 = F 2 .
In the proceeding, it is assumed that S – curve in Figure 3 is an envelope of characteristics. Therefore, Equation (25) is not satisfied even though ψ = ± ( π / 4 + ϕ / 2 ) on this curve. It is possible only if
| q x 2 |
as x 2 0 . In this case, the second term in (24) reduces to the expression 0 at x 2 = 0 .

4. Asymptotic Analysis

The asymptotic analysis below is based on the assumptions that all stress and velocity components are bounded everywhere and all derivatives with respect to x 1 are bounded everywhere, and the solution is represented by power series in x 2 near the curve x 2 = 0 (Figure 1 and Figure 3).

4.1. Plane Strain Deformation

Let ψ f be the value of ψ at x 2 = 0 . For definiteness, it is assumed that ψ f = π / 4 + ϕ / 2 . The case ψ f = ( π / 4 + ϕ / 2 ) can be treated in a similar manner. The condition (22) requires that
q = q 0 + q 1 x 2 α + o ( x 2 α )
as x 2 0 and
0 < α < 1 .
Here q 0 and q 1 are independent of x 2 . The first term in (19) is bounded. Therefore, Equations (27) and (19) combine to give
ψ = ψ f + ψ 1 x 2 α + o ( x 2 α )
as x 2 0 . Turning to Equation (20), the first and third terms of this equation are bounded. The term on its right-hand side is also bounded. Therefore, the second term must be bounded as x 2 0 . Expanding cos 2 ψ in the Taylor series in the vicinity of ψ = ψ f one gets
1 + cos 2 ψ sin ϕ = 2 cot ϕ ( ψ ψ f ) + o [ ( ψ ψ f ) ]
as ψ ψ f . Using (27), (29) and (30) the second term in (20) is represented as
( 1 + cos 2 ψ sin ϕ ) [ 1 + x 2 R ( x 1 ) ] q x 2 = 2 α q 1 ψ 1 cot ϕ x 2 2 α 1 + o ( x 2 2 α 1 )
as x 2 0 . Substituting (27), (29) and (31) into (20) shows that x 2 2 α 1 = O ( x 2 α ) or x 2 2 α 1 = O ( 1 ) as x 2 0 . The former yields α = 1 , which contradicts (28). The latter results in
α = 1 2 .
Then, (27) and (29) become
q = q 0 + q 1 x 2 + o ( x 2 )   and   ψ = ψ f + ψ 1 x 2 + o ( x 2 )
as x 2 0 , respectively. The distribution of p neat the curve is readily found from (4) and (33) as
p = q 0 k cos ϕ sin ϕ + q 1 sin ϕ x 2 + o ( x 2 )
as x 2 0 .
It is now necessary to demonstrate that (33) is compatible with the plastic flow rule. The velocity components are represented as
u 1 = U 0 + U 1 x 2 β + o ( x 2 β )   and   u 2 = V 0 + V 1 x 2 λ + o ( x 2 λ )
as x 2 0 . Here U 0 , U 1 , V 0 and V 1 are independent of x 2 . Moreover, β > 0 and λ > 0 . The first equation in (11) shows that u 2 / x 2 = O ( 1 ) as x 2 0 . Therefore λ = 1 , and the second equation in (34) becomes
u 2 = V 0 + V 1 x 2 + o ( x 2 )
as x 2 0 . Substituting (30), (33), (35) and (36) into the second equation in (11) results in
2 ψ 1 V 1 β cot ϕ x 2 x 2 β 1 = O ( x 2 )
as x 2 0 . It follows from this equation that β = 1 / 2 and the first equation in (35) becomes
u 1 = U 0 + U 1 x 2 + o ( x 2 )
as x 2 0 .
To the best of the authors’ knowledge, two semi-analytical solutions for the double slip and rotation are available [16,17]. The asymptotic analysis of these solutions shows that they satisfy (33), (34), (36) and (38). Several semi-analytic solutions for the double-shearing model, which is another widely used model of pressure-dependent plasticity, have been derived in [11,12,13] and [16]. The asymptotic analysis of these solutions also shows that they satisfy (33), (34), (36), and (38).

4.2. Axisymmetric Deformation.

It is seen from (19), (20), (23) and (24) that the left-hand sides of the equations under plane strain and axial symmetry coincide. The right-hand sides of these equations have no effect on asymptotic analysis. Therefore, the equations in (33) and (34) are valid for axisymmetric deformation. The circumferential stress is readily determined from the first equation in (14).
Turning to the velocity equations in (18), the derivative u 2 / x 2 in the second equation can be eliminated using the first equation. As a result,
u 1 x 1 + [ 1 + x 2 R ( x 1 ) ] u 2 x 2 = ( u 2 sin γ u 1 cos γ ) r [ 1 + x 2 R ( x 1 ) ] u 2 R ( x 1 ) , ( cos 2 ψ sin ϕ 1 ) u 2 x 1 + ( cos 2 ψ sin ϕ + 1 ) [ 1 + x 2 R ( x 1 ) ] u 1 x 2 2 sin 2 ψ sin ϕ u 1 x 1 + = [ 1 + x 2 R ( x 1 ) ] ( u 1 cos γ u 2 sin γ ) sin 2 ψ r sin ϕ + ( cos 2 ψ sin ϕ ) u 1 R ( x 1 ) sin ϕ + 2 u 2 sin 2 ψ R ( x 1 ) sin ϕ 2 Ω [ 1 + x 2 R ( x 1 ) ] .
The left-hand sides of these equations coincide with the corresponding left-hand sides of the equations in (11). Therefore, Equations (36) and (38) are also valid.

5. Conclusions

The exact asymptotic representation of solutions in the vicinity of envelopes of characteristics has been found for the double slip and rotation model under plane strain and axially symmetric conditions. This representation coincides with the behavior of specific solutions found in [11,12,13,16,17]. This representation also coincides with that for isotropic rigid plastic solids obeying an arbitrary smooth yield criterion [18], and for some anisotropic yield criteria [19]. It is seen from Equations (33), (34) and (38) that the stress and velocity fields are singular in the vicinity of the envelope of characteristics, in the sense that some derivatives of the stress and velocity components approach infinity near the envelope.
The applied aspect of the asymptotic representation of solutions derived in the present paper is twofold. Firstly, as known from other models of rigid plasticity, numerical solutions using traditional finite elements do not converge if the exact solution is described by equations similar to (33), (34) and (38) [32,33]. The asymptotic representation given in (33), (34) and (38) can be used in conjunction with the extended finite element method, outlined in [29], for developing an accurate numerical approach for solving corresponding boundary value problems.
Secondly, the asymptotic representation given in (33), (34) and (38) can be used in conjunction with models for soil–structure interfaces [21,22,23]. In particular, an envelope of characteristics usually coincides with the interface between a structure and material, for example between a structure and soil. In this case, experiment confirms some mathematical features of the general theoretical solution. If an envelope of characteristics coincides with the interface then the regime of sliding occurs independently of the friction law chosen. This means that the friction stress is independent of the roughness of the contact surface. It is in general known that the surface roughness affects the material behavior near the interface [34,35]. However, this effect is usually significant only if the surface roughness is small enough. For example, in the case of the interface between Yongdinghe sand and low carbon steel the friction stress is practically independent of the surface roughness if R n > R c r = 0.1 where R n is the relative interface roughness and R c r is its critical value [35]. In [35], the soil–structure interfaces are classified as smooth or rough depending on the value of R n . When R n < R c r the interface is smooth. If R n > R c r then the interface is rough. In the latter case, strain localization occurs near the interface. This feature of material behavior is in agreement with (38) as this equation predicts that the shear strain rate approaches infinity as x 2 0 . Moreover, it has been found in [35] that the residual friction angle is close to the internal friction angle of sand if R n > R c r . This feature is in agreement with the second equation in (33). Therefore, it is reasonable to assume that the condition that an envelope of characteristics coincides with the soil–structure interface is adequate when R n > R c r . In this case, it is natural to expect the similarity between the behavior of soil and interface. This similarity has been investigated in [23].

Author Contributions

Conceptualization, S.A.; writing, Y.W.; formal analysis, E.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Hebei Province of China with Grant No. E2019202224 and the RFBR (Project No.18-51-76001).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cox, G.M.; Thamwattana, N.; McCue, S.W.; Hill, J.M. Coulomb-Mohr granular materials: Quasi-static flows and the highly frictional limit. Appl. Mech. Rev. 2008, 61, 060802. [Google Scholar] [CrossRef]
  2. Goddard, J.D. Continuum modeling of granular media. Appl. Mech. Rev. 2014, 66, 050801. [Google Scholar] [CrossRef] [Green Version]
  3. Hui, K.; Haff, P.K.; Ungar, J.E.; Jackson, R. Boundary conditions for high-shear grain flows. J. Fluid Mech. 1984, 145, 223–233. [Google Scholar] [CrossRef] [Green Version]
  4. Gutt, G.M.; Haff, P.K. Boundary conditions on continuum theories of granular flow. Int. J. Multiph. Flow 1991, 17, 621–634. [Google Scholar] [CrossRef]
  5. Jenkins, J.T. Boundary conditions for rapid granular flow: Flat, frictional walls. J. Appl. Mech. 1992, 59, 120–127. [Google Scholar] [CrossRef]
  6. Savage, S.B.; Dai, R. Studies of granular shear flows: Wall slip velocities, “layering” and self-diffusion. Mech. Mater. 1993, 16, 225–238. [Google Scholar] [CrossRef]
  7. Zheng, X.M.; Hill, J.M. Boundary effects for Couette flow of granular materials: Dynamical modelling. Appl. Math. Model. 1996, 20, 82–92. [Google Scholar] [CrossRef] [Green Version]
  8. Ehsan, B.; Ali, N.; Shima, Z. A state boundary surface model for improving the dilatancy simulation of granular material in reinforced anchors. Arab. J. Geosci. 2017, 10, 281–293. [Google Scholar]
  9. Yang, L.; Padding, J.T.; Kuipers, J.A.M. Partial slip boundary conditions for collisional granular flows at flat frictional walls. AIChE J. 2017, 63, 1853–1871. [Google Scholar] [CrossRef] [Green Version]
  10. Sarno, L.; Carleo, L.; Papa, M.N.; Villani, P. Experimental investigation on the effects of the fixed boundaries in channelized dry granular flows. Rock Mech. Rock Eng. 2018, 51, 203–225. [Google Scholar] [CrossRef]
  11. Pemberton, C.S. Flow of imponderable granular materials in wedge-shaped channels. J. Mech. Phys. Solids 1965, 13, 351–360. [Google Scholar] [CrossRef]
  12. Marshall, E.A. The compression of a slab of ideal soil between rough plates. Acta Mech. 1967, 3, 82–92. [Google Scholar] [CrossRef]
  13. Spencer, A.J.M. Compression and shear of a layer of granular material. J. Eng. Math. 2005, 52, 251–264. [Google Scholar] [CrossRef]
  14. Spencer, A.J.M. A theory of the kinematics of ideal soils under plane strain conditions. J. Mech. Phys. Solids 1964, 2, 337–351. [Google Scholar] [CrossRef]
  15. Harris, D. A hyperbolic augmented elasto-plastic model for pressure-dependent yield. Acta Mech. 2014, 225, 2277–2299. [Google Scholar] [CrossRef] [Green Version]
  16. Alexandrov, S.; Harris, D. Comparison of solution behaviour for three models of pressure-dependent plasticity: A simple analytical example. Int. J. Mech. Sci. 2006, 48, 750–762. [Google Scholar] [CrossRef] [Green Version]
  17. Alexandrov, S.; Harris, D. An exact solution for a model of pressure-dependent plasticity in an un-steady plane strain process. Eur. J. Mech. A Solids 2014, 29, 966–975. [Google Scholar] [CrossRef]
  18. Alexandrov, S.; Richmond, O. Singular plastic flow fields near surfaces of maximum friction stress. Int. J. Non-Linear Mech. 2001, 36, 1–11. [Google Scholar] [CrossRef]
  19. Alexandrov, S.; Jeng, Y.R. Singular rigid/plastic solutions in anisotropic plasticity under plane strain conditions. Cont. Mech. 2013, 25, 685–689. [Google Scholar] [CrossRef]
  20. Alexandrov, S.; Lang, L.; Lyamina, E.; Date, P.P. Solution behavior near envelopes of characteristics for certain constitutive equations used in the mechanics of polymers. Materials 2019, 12, 2725. [Google Scholar] [CrossRef] [Green Version]
  21. Desai, C.S.; Zaman, M.M. Thin-layer element for interfaces and joints. Int. J. Numer. Anal. Methods 1984, 8, 19–43. [Google Scholar] [CrossRef]
  22. De Gennaro, V.; Frank, R. Elasto-plastic analysis of the interface behavior between granular media and structure. Comput. Geotech. 2002, 29, 547–572. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, H.; Song, E.; Ling, H.I. Constitutive modeling of soil-structure interface through the concept of critical state soil mechanics. Mech. Res. Commun. 2006, 33, 515–531. [Google Scholar] [CrossRef]
  24. Gens, A.; Carol, I.; Alonso, E.E. An interface element formulation for the analysis of soil reinforcement interaction. Comput. Geotech. 1989, 7, 133–151. [Google Scholar] [CrossRef]
  25. Boulo, M. Basic features of soli structure interface behavior. Comput. Geotech. 1989, 7, 115–131. [Google Scholar] [CrossRef]
  26. Ghionna, V.H.; Mortara, G. An elastoplastic model for sand-structure interface behavior. Geotechnique 2002, 52, 41–50. [Google Scholar] [CrossRef]
  27. Kaliakin, V.N. Insight into deficiencies associated with commonly used zero-thickness interface elements. Comput. Geotech 1995, 17, 225–232. [Google Scholar] [CrossRef]
  28. Villard, P. Modelling of interface problems by the finite element method with considerable displacements. Comput. Geotech. 1996, 19, 23–45. [Google Scholar] [CrossRef]
  29. Fries, T.P.; Belytschko, T. The extended/generalized finite element method: An overview of the method and its applications. Int. J. Numer. Meth. Eng. 2010, 84, 253–304. [Google Scholar] [CrossRef]
  30. Malvern, L.E. Introduction to the Mechanics of a Continuous Medium; Prentice Hall: Englewood Cliffs, NJ, USA, 1969. [Google Scholar]
  31. Spencer, A.J.M. Deformation of ideal granular materials. In Mechanics of Solids, the Rodney Hill 60th Anniversary Volume; Hopkins, H.G., Sewell, M.J., Eds.; Pergamon Press: Oxford, UK, 1982; pp. 607–652. [Google Scholar]
  32. Chen, J.C.; Pan, C.; Rogue, C.M.O.L.; Wang, H.P. A Lagrangian reproducing kernel particle method for metal forming analysis. Comp. Mech. 1998, 22, 289–307. [Google Scholar] [CrossRef]
  33. Facchinetti, M.; Miszuris, W. Analysis of the maximum friction condition for green body forming in an ANSYS environment. J. Eur. Ceram. Soc. 2016, 36, 2295–2302. [Google Scholar] [CrossRef] [Green Version]
  34. Uesugi, M.; Kishida, H. Frictional resistance at yield between dry sand and mild steel. Soils Found. 1986, 26, 139–149. [Google Scholar] [CrossRef] [Green Version]
  35. Hu, L.; Pu, J.L. Testing and modeling of soil-structure interface. J. Geotech. Geoenviron. 2003, 120, 851–860. [Google Scholar] [CrossRef]
Figure 1. Orthogonal coordinate system ( x 1 , x 2 ) near a given curve (S–line) under plane strain conditions.
Figure 1. Orthogonal coordinate system ( x 1 , x 2 ) near a given curve (S–line) under plane strain conditions.
Applsci 10 03220 g001
Figure 2. Orientation of the major principal stress.
Figure 2. Orientation of the major principal stress.
Applsci 10 03220 g002
Figure 3. Orthogonal coordinate system ( x 1 , x 2 ) near a given curve (S–line) under axially symmetric conditions.
Figure 3. Orthogonal coordinate system ( x 1 , x 2 ) near a given curve (S–line) under axially symmetric conditions.
Applsci 10 03220 g003
Figure 4. Mohr–Coulomb yield locus under axially symmetric conditions.
Figure 4. Mohr–Coulomb yield locus under axially symmetric conditions.
Applsci 10 03220 g004

Share and Cite

MDPI and ACS Style

Wang, Y.; Alexandrov, S.; Lyamina, E. Solution Behavior in the Vicinity of Characteristic Envelopes for the Double Slip and Rotation Model. Appl. Sci. 2020, 10, 3220. https://doi.org/10.3390/app10093220

AMA Style

Wang Y, Alexandrov S, Lyamina E. Solution Behavior in the Vicinity of Characteristic Envelopes for the Double Slip and Rotation Model. Applied Sciences. 2020; 10(9):3220. https://doi.org/10.3390/app10093220

Chicago/Turabian Style

Wang, Yao, Sergei Alexandrov, and Elena Lyamina. 2020. "Solution Behavior in the Vicinity of Characteristic Envelopes for the Double Slip and Rotation Model" Applied Sciences 10, no. 9: 3220. https://doi.org/10.3390/app10093220

APA Style

Wang, Y., Alexandrov, S., & Lyamina, E. (2020). Solution Behavior in the Vicinity of Characteristic Envelopes for the Double Slip and Rotation Model. Applied Sciences, 10(9), 3220. https://doi.org/10.3390/app10093220

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