`Mathematical Problems in EngineeringVolume 2012 (2012), Article ID 709178, 9 pageshttp://dx.doi.org/10.1155/2012/709178`
Research Article

## Design of an Annular Disc Subject to Thermomechanical Loading

1A.Yu. Ishlinskii Institute for Problems in Mechanics, Russian Academy of Sciences, 101-1 Prospect Vernadskogo, Moscow 119526, Russia
2Department of Mechanical Engineering and Advanced Institute for Manufacturing with High-tech Innovations, National Chung Cheng University, 168 University Road, Chia-Yi 62102, Taiwan

Received 27 July 2012; Accepted 3 September 2012

Academic Editor: Valery Yakhno

Copyright © 2012 Sergei Alexandrov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Two solutions to design a thin annular disc of variable thickness subject to thermomechanical loading are proposed. It is assumed that the thickness of the disc is everywhere sufficiently small for the stresses to be averaged through the thickness. The state of stress is plane. The initiation of plastic yielding is controlled by Mises yield criterion. The design criterion for one of the solutions proposed requires that the distribution of stresses is uniform over the entire disc. In this case there is a relation between optimal values of the loading parameters at the final stage. The specific shape of the disc corresponds to each pair of such parameters. The other solution is obtained under the additional requirement that the distribution of strains is uniform. This solution exists for the disc of constant thickness at specific values of the loading parameters.

#### 2. Statement of the Problem

Consider a thin annular disc of outer radius b and inner radius a inserted into a rigid container of radius a. It is convenient to introduce a cylindrical coordinate system (r, , z) with its z-axis coinciding with the axis of symmetry of the disc. The initial thickness of the disc, h, is a function of r. The disc is subject to thermal loading by a uniform temperature field varying with the time. The disc has no stress at the initial temperature. Uniform pressure varying with the time is applied over the inner radius of the disc. The outer radius is fixed to the container. It is evident that the problem is axisymmetric. In particular, the solution is independent of . Moreover, the normal stresses in the cylindrical coordinates, , , and are the principal stresses. It is also assumed that the state of stress is two-dimensional, . The pressure applied, thermal expansion caused by a rise of temperature, and the constraints imposed on the disc affect the initial zero-stress state. It is assumed that the rise of temperature above the reference state, T, and the pressure over the inner radius are monotonically nondecreasing functions of the time, t. The boundary conditions are at and at . Here u is the radial displacement, is a function of the time, and is a constant introduced for further convenience. The circumferential displacement vanishes everywhere.

It is assumed that the thickness of the disc is everywhere sufficiently small for the stresses to be averaged through the thickness. In this case the only nontrivial equilibrium equation becomes The total radial, , and circumferential, , strains are defined by where the superscript denotes the thermal portions of the total strains, the superscript e the elastic portions of the total strains, and the superscript p the plastic portions of the total strains. It follows from Hooke’s law that where E is Young’s modulus and is Poisson’s ratio. The thermal portions of the total strains are given by where is the thermal coefficient of linear expansion. In the plastic range, Mises yield criterion is adopted. For the problem under consideration this criterion reduces to where is the yield stress in tension, a material constant for perfectly plastic materials. This quantity is also involved in (2.1). The associated flow rule is written in terms of the strain rate components. A consequence of this rule is where and are the plastic portions of the total radial and circumferential strain rates. Another essential equation following from the associated flow rule expresses plastic incompressibility, , where is the plastic portion of the total axial strain rate. This equation serves to determine and is not important for the present solution. At small strains,

According to the design criterion proposed by Michell (see [3]), all of the structural elements must be strained by exactly the same strain magnitude in either simple tension or pure compression. This criterion can be too restrictive for the structure under consideration. Therefore, in the present paper two design criteria are adopted. First, it is required that an equistressed state occurs in the entire disc. Then, the possibility to obtain a uniform distribution of strains is explored.

#### 3. Restriction on Thickness Variation

The same magnitude of the elastic portion of strains can be obtained if and only if the distribution of stress components is uniform. Then, it follows from (2.3) that where is constant. Let be the thickness of the disc at . Then, the solution of (3.1) satisfying this condition is Note that this function is often adopted in studies devoted to analysis of thin discs, for example [46]. The uniform distribution of stresses is only required in the final stage of loading. Using (3.2) the equation of equilibrium (2.3) for intermediate stages becomes

#### 4. Thermoelastic Solution

At the beginning of the process of loading the entire disc is elastic. At this stage, Eliminating u between these two equations, using (2.5) and (2.6), and taking into account that T is independent of r yield Eliminating the stress in (4.2) by means of (3.3) gives It is convenient to introduce the dimensionless radius by . Then, the general solution of (4.3) is Where and are constants of integration and Substituting (4.4) into (2.5) determines . Then, using this expression for and (2.6) the radial displacement can be found from the equation . As a result, where and . Substituting the boundary conditions (2.1) and (2.2) into (4.4) and (4.6) leads to where .

#### 5. Thermoelastic-Plastic Solution for Design

The yield criterion (2.7) is satisfied by the following substitution: where is a function of and . Substituting (5.1) into (3.3) results in The zone where the yield criterion is satisfied should occupy the entire disc at the final stage. As it has been mentioned before, the design criterion chosen is satisfied if and only if the distribution of stress is uniform at this stage. Therefore, it should be uniform over the domain where (5.2) is valid. It follows from (5.1) that the condition that the distribution of the stresses in the plastic zone is uniform is equivalent to the condition that is independent of . It is evident that the general solution of (5.2) does not satisfy this requirement. However, this equation has a special solution in the form , where It is seen from (2.1), (5.1) and (5.3) that this special solution takes place if and only if

Let R be the dimensionless radius of the elastic/plastic boundary. The general solution (4.4) and (4.6) is valid in the elastic zone. However, A and B are not given by (4.7). The stresses and are continuous across the elastic/plastic boundary. Therefore, it follows from (4.4) and (5.1) that The boundary condition (2.2) combined with (4.6) gives Solving (5.5) for and results in Thus, the quantities and are independent of . Eliminating A and B in (5.6) by means of (5.7) leads to

It is convenient to introduce the following quantities and . Then, (2.8) becomes Since the stresses are constant in the plastic zone, the elastic strain rates vanish. Moreover, it follows from (2.6) that and . Therefore, the total strain rates in the equation of compatibility can be replaced with their plastic portions. Then, using the definition for and this equation is reduced to Substituting (5.1) at into (5.9) and then eliminating in (5.10) yield Here, (5.3) has been used to eliminate . Equation (5.11) can be immediately integrated to give where is a function of integration. Introduce the notation . Note that w is proportional to the radial velocity. Using (5.12) the value of w on the plastic side of the elastic/plastic boundary is determined as Differentiating (4.6) with respect to yields Thus, the value of w on the elastic side of the elastic/plastic boundary is Since , it follows from (5.13) and (5.15) that This equation can be rewritten in the following equivalent form: Eliminating here the derivatives and by means of (5.7) gives Eliminating here by means of (5.3) yields Substituting (5.19) into (5.12) gives Integrating with respect to determines the circumferential plastic strain as where is an arbitrary function of . This function should be found using the condition at the elastic/plastic boundary. Then, it follows from (5.21) that Substituting (5.1) at into (5.9), eliminating by means of (5.3), and integrating with the respect to using the condition that when give

#### 6. Design of the Disc

The solution found can be used to search for two kinds of optimal conditions. In particular, it is possible to search for a uniform distribution of stresses at the final stage of loading. This kind of design requires that the plastic zone occupies the entire disc. The stresses at any point of the disc are given by (5.1), where should be replaced with . Putting in (5.8) determines the value of at which the entire disc becomes plastic Using (4.5), (5.3), (5.4), (5.7), and (6.1) it is possible to find a relation between the two optimal loading parameters at the final stage, and , at a given value of n. The distribution of the elastic and thermal portions of the strain tensor is uniform at these values of the parameters. However, the plastic portion of the strain tensor varies with the radius according to (5.22) and (5.23). Therefore, not all of the requirements of Michell structures are satisfied. The following solution enables the total strain distribution to be uniform. It is evident from (5.12) that it is possible if and only if . Then, it follows from (5.19) that or . These conditions along with (4.5) provide two equations for n. The equation corresponding to has no solution. The other equation gives . In this case . Thus, the distribution of strains in the disc of constant thickness is uniform if , as follows from (5.3) and (5.4). Moreover, it is seen from (5.3) and (5.7) that at and . Therefore, (5.8) does not provide any relation between R and . The physical meaning of this feature of the solution is that the plastic zone simultaneously occupies the entire disc of any size. Thus, this design satisfies the criterion adopted in [1]. The corresponding value of can be found from the thermoelastic solution. In particular, since is given, and in (4.7) are solely dependent of . Therefore, replacing A and B in (4.4) with and , respectively, and putting determine the stresses as functions of . Finally, substituting these functions into the yield criterion (2.7) gives the equation for the value of at which the entire disc becomes plastic.