- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 709178, 9 pages

http://dx.doi.org/10.1155/2012/709178

## Design of an Annular Disc Subject to Thermomechanical Loading

^{1}A.Yu. Ishlinskii Institute for Problems in Mechanics, Russian Academy of Sciences, 101-1 Prospect Vernadskogo, Moscow 119526, Russia^{2}Department 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.

#### 1. Introduction

Thin annular discs subject to various loading conditions are a class of commonly used structures in mechanical engineering. Therefore, there is a vast amount of literature on this topic. These studies can be conveniently divided into two groups, namely, analysis and design. The present paper deals with the design of thin discs. Therefore, previous works solely related to the analysis of discs are not considered here. Reviews of some results on the design of thin discs are provided in [1, 2]. An analytical solution to design an annular disc of variable thickness under internal and external pressures has been proposed in [1]. It has been assumed that the initiation of plastic yielding is controlled by Tresca yield criterion. The goal of the design has been to find conditions under which the yield criterion is simultaneously satisfied at all points of the disc. A disc rotating in a temperature field has been considered in [2]. The objective function in the problem is the disc weight. A numerical method has been used to solve the problem. The present paper concerns with an annular disc inserted into a rigid container and subject to thermomechanical loading. One of the loading parameters is uniform temperature and the other is pressure over the inner radius of the disc. Both vary with the time. Two design criteria are adopted. The requirement of one of the design criteria is that the distribution of stresses is uniform over the entire disc at the final stage. Using this criterion, a relation between optimal loading parameters is obtained. The solution is given in analytical form. The additional requirement of the other design criterion is that the state of strain is uniform. The same requirement is adopted to design Michell structures (see [3]). It is shown that in the case of the disc under consideration there is the unique solution for the second design problem. In particular, the thickness of the disc is constant.

#### 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 [4–6]. 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.

#### 7. Conclusions

New solutions for design of a thin annular disc subject to thermomechanical loading have been proposed. Two design criteria have been adopted. One of the criteria requires that the state of stress is uniform at the final stage of loading. This criterion leads to a relation between optimal values of the loading parameters for each specific shape of the disc. A more restrictive criterion additionally requires, by analogy to Michell structures, that the state of strain is uniform at the final stage. Application of this criterion has shown that the state of stress and strain required appears in the disc of constant thickness at certain values of the loading parameters. This design also satisfies the requirements formulated in [1]. Possible developments of the approach proposed include plastic anisotropy, pressure-dependency of the yield criterion, and variation of some material properties along the radius.

#### Acknowledgment

The research described in this paper has been supported by the Grants NSC-99-2218-E-194-003-MY3, RFBR-11-01-93001, and GK No. 11.519.11.3015. A part of this work was done while the first author was with National Chung Cheng University (Taiwan) as a Research Scholar under the recruitment program supported by the National Science Council of Taiwan (contact 99-2811-E-194-009).

#### References

- K. Hein and M. Heinloo, “The design of nonhomogeneous equi-strength annular discs of variable thickness under internal and external pressures,”
*International Journal of Solids and Structures*, vol. 26, no. 5-6, pp. 617–630, 1990. View at Google Scholar · View at Scopus - B. Farshi, H. Jahed, and A. Mehrabian, “Optimum design of inhomogeneous non-uniform rotating discs,”
*Computers and Structures*, vol. 82, no. 9-10, pp. 773–779, 2004. View at Publisher · View at Google Scholar · View at Scopus - P. Dewhurst, “Analytical solutions and numerical procedures for minimum-weight Michell structures,”
*Journal of the Mechanics and Physics of Solids*, vol. 49, no. 3, pp. 445–467, 2001. View at Publisher · View at Google Scholar · View at Scopus - J. Chakrabarty, “Elastic-plastic expansion of a circular hole in an infinite plate of variable thickness,”
*International Journal of Mechanical Sciences*, vol. 13, no. 5, pp. 439–444, 1971. View at Google Scholar · View at Scopus - J. W. Phillips and M. Schrock, “Note on shear stresses in accelerating disks of variable thickness,”
*International Journal of Mechanical Sciences*, vol. 13, no. 5, pp. 445–449, 1971. View at Google Scholar · View at Scopus - A. N. Sherbourne and D. N. S. Murthy, “Stresses in disks with variable profile,”
*International Journal of Mechanical Sciences*, vol. 16, no. 7, pp. 449–459, 1974. View at Google Scholar · View at Scopus