#### Abstract

The paper presents an exact analytic solution for a class of elastic-plastic models with damage evolution. The boundary value problem consists of a planar deformation comprising the simultaneous shearing and expansion of a hollow cylindrical specimen of material and involves a bimaterial interface at which the materials stick to each other. With no loss of generality for understanding the qualitative behaviour of the solution near the bimaterial interface, an extreme case when the hard material is rigid is considered. The solution is reduced to a transcendental equation for the value of the equivalent plastic strain at the bimaterial interface. This equation predicts that the equivalent plastic strain attains a maximum under certain conditions. The existence of the solution of the boundary value problem depends on the value of the damage parameter at fracture, which is a material constant. In particular, if this value is larger than the value of the damage parameter at the bimaterial interface corresponding to the maximum possible value of the equivalent strain at this interface, then no solution exists. Experimental data available in the literature are used to assess whether Lemaitre’s model is applicable.

#### 1. Introduction

The ductile damage evolution equations of continuum mechanics can be conveniently divided into three groups: (i) uncoupled models, (ii) partly coupled models, and (iii) coupled models. In the case of the models of group (i) the damage evolution equation is solved having found the solution for the system of equations of plasticity theory. A recent review of models of this group is provided in [1]. The models of group (iii) include the fully coupled system of equations of plasticity theory and the damage evolution equation. In particular, the material is plastically compressible. A recent review of models of this group is provided in [2]. In the present paper, a wide class of models of group (ii) for elastic-plastic materials is adopted. These models containing one scalar measure of damage are based on the pressure-independent yield criterion and the associated flow rule for plastic strain rates. The evolution of the damage parameter is controlled by an additional equation. Such models are presented, for example, in [3–7]. The yield criterion and elastic moduli depend on the damage parameter. Many structures and machine parts are composed of two or more different materials and, therefore, contain bimaterial interfaces. Solution behaviour in the vicinity of such interfaces may reveal some qualitative features that provide an assessment of the model adopted to find the solution. For example, using a rigid perfectly plastic material model, it has been shown in [8] that the velocity field is singular in the vicinity of some bimaterial interfaces. This singularity in the velocity field is not compatible with the class of damage evolution models reviewed in [9]. In order to assess the mathematical model, it is possible to study any boundary value problem, even if its practical application is not feasible. It is more important to find the exact solution of the problem without relaxing the boundary conditions and field equations. Such a boundary value problem has been proposed and solved for a rigid plastic linear hardening material model in [10]. Then, this boundary value problem has been successfully used in [11, 12] to study various qualitative features of solutions for several material models. However, all of these solutions are for rigid plastic materials with no damage evolution. In the present paper, the same boundary value problem is combined with a class of elastic-plastic models with damage evolution at small strains. The solution is reduced to a transcendental equation for the value of the equivalent plastic strain at the bimaterial interface. The solution to this equation predicts that this value attains a maximum under certain conditions. The existence of the solution of the boundary value problem depends on the value of the damage parameter at fracture, which is supposed to be a constitutive parameter. In particular, if the value of the damage parameter at the bimaterial interface corresponding to the maximum possible value of the equivalent plastic strain at this interface is smaller than the damage parameter at fracture, then no solution exists. The model proposed by Lemaitre [3] is studied in some detail. In particular, experimental data provided in [13] are used to assess whether the model for this specific set of constitutive parameters is applicable (in the sense that the solution exists).

#### 2. Statement of the Problem

Consider an infinite elastic-plastic circular cylinder of internal radius and external radius subject to the action of tangential stresses uniformly distributed over its external radius. It is assumed that a rigid rod of radius is inserted into the hole of the cylinder. The rod is motionless and fixed to the cylinder. Due to this system of loading and constraints, the cylinder is twisted. It is convenient to introduce a cylindrical polar coordinate system with its *z*-axis coinciding with the axis of symmetry of the cylinder. The boundary conditions are
for and
for . Here is the circumferential displacement and is the angular displacement of the external surface of the cylinder.

A typical dependence of the shear yield stress, , on the equivalent plastic strain, , and damage parameter, , is [3, 4, 7] Here is the initial shear yield stress when and . The function satisfies the conditions and for any . The equivalent plastic strain is defined by the equation Here is the time and is the convected derivative. The equivalent plastic strain rate, , is defined by where are the plastic strain rate components. The yield criterion is where are the deviatoric components of the stress tensor. The flow rule associated with (6) is where is a scalar factor of proportionality. The elastic strains, , obey Hooke’s law. The total strains are In the case of isotropic models with a single scalar measure of damage typical damage evolution equations are [3, 4, 7] Here is the hydrostatic stress, is the threshold equivalent plastic strain for the initiation of damage, is constant, and is a given function of its arguments. The fracture criterion is where is a material constant whose value for metallic materials lies in the range [7]. Equations (3) to (11) along with Hooke’s law represent a typical isotropic elastic-plastic material model with damage evolution at small strains. These equations should be supplemented with the equilibrium equations.

#### 3. Solution

In this section, an exact solution of the equations of the previous section is derived in the cylindrical coordinates assuming that the solution is independent of and *z*. It is then demonstrated that the solution is compatible with the boundary conditions (1) and (2). Suppose that all components of the stress tensor vanish except for the shear component . The yield criterion (6) is satisfied if
The only nontrivial equilibrium equation is
Denote the radius of the elastic-plastic boundary by , . The material just on the elastic side of the boundary must be on the point of yielding. Therefore, at . Integrating (13) with the use of this condition gives
It follows from Hooke’s law that
where is the modulus of shear.

Consider the elastic region, . In this region, and , where is the shear modulus of undamaged material. Therefore, substituting (14) into (15) yields

In the plastic region, , combining (12) and (14) gives A typical assumption is that [3, 4, 7]. Then, substituting (12) into (15) results in Since the only nonzero stress component is , it follows from (7) that the only nonzero plastic strain rate component is . Then, (5) reduces to It has been taken into account here that , as follows from (12). Since strains are small, and . Therefore, integrating (19) with respect to gives The total shear strain in the plastic region is found from (8), (18), and (20) as Since , (9) becomes where . Assuming that at the initial instant, (22) can be integrated to give Here where is a dummy variable of integration. It follows from (3), (17), and (23) that where is the function inverse to . It has been assumed here that . Substituting (3) into (21) yields For the problem under consideration, the total shear strain is defined by Combining (26) and (27) gives It is convenient to introduce new independent variables as Then, where . Using (29), (30), (25) and (28) transforms to respectively. Eliminating between (31) and (32) leads to where . Let be the value of at . Since at , the solution to (33) is The value of should be found using the solution in the elastic region. Substituting (27) into (16) and using (29) and (30) give The solution of this equation satisfying the boundary condition (2) is It follows from (29) that if . Therefore, the value of is determined from (36) as Substituting (37) into (34) gives Let be the value of at . Then, combining (38) and the boundary condition (1) results in It follows from (29) that at , and, therefore, at . Then, (31) shows that . Using this equation to eliminate in (39) yields It follows from (23) that the value of at is given by It has been assumed here that .

#### 4. Fracture Initiation and Nonexistence of the Solution

Equation (40) determines as a function of . It follows from this equation that The right hand side of this equation vanishes if or Using (31) and (32), the left hand side of (43) can be represented in the form where . Consider the coefficient of in (45). It follows from the definition for the function that is equal to the value of at . Therefore, and the denominator of the coefficient of in (45) is smaller than . On the other hand, it is evident that its numerator is larger than . Therefore, the coefficient is larger than 1. However, the coefficient of in (43) is smaller than 1. Hence (43) has no solution. Let be the solution to (44). The corresponding values of and are determined from (40) and (41) at . Expanding the left hand side of (42) in a series in the vicinity of gives to leading order. Here is constant. Integrating this equation with the use of the condition at results in to leading order. This solution can be rewritten in the form . It is evident from this expression that cannot exceed the value of . Thus if (43) has a solution, then the solution given in the previous section does not exist for . On the other hand, according to the fracture criterion (11), the initiation of fracture at the bimaterial interface occurs if . The corresponding value of is determined from (41) as . Hence the solution to the boundary value problem of the previous section, including the fracture criterion, exists if and only if In particular, the solution predicts the initiation of fracture at the bimaterial interface at . It will be shown in the next section that the inequality (46) is not satisfied for some specific damage evolution laws and constitutive parameters used in applications.

#### 5. Illustrative Example

In the case under consideration the damage evolution model [3, 4] is obtained if where is constant. Therefore, evaluating the integral in (24) gives Then, it follows from (23) that The equation of the strain hardening curve for the steel tested in [13], is Using experimental data provided in [13] it is possible to find that Substituting (48)–(50) into (31) gives . Hence, (44) becomes . Its solution is . Substituting this value of into (41) with the use of (48) gives . This means that the solution to the boundary value problem without using the fracture criterion breaks down when . On the other hand, the fracture criterion (11) is not satisfied as follows from (50). Therefore, it is necessary to conclude that the damage evolution model adopted is not compatible with the constitutive parameters recommended in [13].

#### 6. Conclusions

An exact analytic solution for an elastic-plastic model coupled with a damage evolution equation has been found. Special attention has been devoted to the qualitative behavior of the solution in the vicinity of the bimaterial interface. Even though the boundary value problem solved is not feasible for practical applications, the main result obtained is useful for assessing damage evolution models by mathematical methods. In particular, it has been shown that the boundary value problem may have no solution because the equation relating the equivalent plastic strain at the bimaterial interface and one of the boundary conditions has no solution. This conclusion is not valid if the initiation of fracture occurs at a lower value of this strain. Therefore, having specific damage evolution equation and constitutive parameters, it is possible to verify that the model is compatible with conditions at bimaterial interfaces where shear straining prevails. As an example, the constitutive equations proposed in [13] have been used, and it has been demonstrated that the aforementioned equation has no solution. This means that some boundary value problems may have no solution for this model.

Even though (9) has been considered explicitly, it is evident that (10) can be treated in a similar manner. It is just necessary to replace (23) with where , , and .

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The research described in this paper has been supported by the Grants RFBR-11-01-00987 and NSH-3842.2012.1.