Research Article | Open Access

Anarbay Kudaykulov, Erkin Arinov, Nurlybek Ispulov, Abdul Qadir, Kalamkas Begaliyeva, "Numerical Study of a Thermally Stressed State of a Rod", *Advances in Mathematical Physics*, vol. 2019, Article ID 8986010, 9 pages, 2019. https://doi.org/10.1155/2019/8986010

# Numerical Study of a Thermally Stressed State of a Rod

**Academic Editor:**Zine El Abiddine Fellah

#### Abstract

This paper considers a new method based on the law of energy conservation for the study of thermo-stress-strain state of a rod of limited length with simultaneous presence of local heat fluxes, heat exchanges, and thermal insulation. The method allows determining the field of temperature distribution and the three components of deformations and stresses, as well as the magnitude of the rod elongation and the resulting axial force with an accuracy of satisfying the energy conservation laws. For specific initial data, all the sought-for ones are determined numerically with high accuracy. We found that all solutions satisfy the laws of energy conservation.

#### 1. Introduction

This work is devoted to the study of the thermo-stress-strain state of a rod of constant cross section and of limited length. In this case, a heat flux of constant intensity is supplied to a closed local surface. The rest of the side surface of the rod is fully thermally insulated. Through the cross-sectional area of the two ends of the rod, convective heat exchange with the environment occurs. In this case, the coefficients of heat transfer and the temperature of the surrounding medium at the two ends of the rod are different. To determine the temperature field, the energy conservation law is used in the form of the total heat energy functional taking into account the existing dissimilar types of heat sources, physical and mechanical properties of the rod material, and its geometrical dimensions. Using this, the law of temperature distribution along the length of the rod under study is constructed. They determine the magnitude of the elongation of the rod, in the case of pinching one end. If both ends of the rod are clamped, then the magnitude of the resulting axial force is determined. The laws of distribution of thermoelastic, temperature, and elastic components of deformations and stresses are also determined depending on the presence of local heat flux, thermal insulation of heat exchanges, the geometry of the rod, and the physic-mechanical properties of the material of the rod. To determine the displacement field, the energy conservation law is used in the form of the potential energy functional of elastic deformation taking into account the presence of a temperature field. Further, the displacement field is determined taking into account the actual operating conditions. The developed program allows varying the values of the source data.

In [1], the fundamental laws of the theory of thermoelasticity for deformable solids are given. In [2, 3], the results of a numerical study of a thermo-stress-deformed state of a rod under the action of laser beams are presented. In this case, the finite element method was used. The application of the finite element method is given in [4]. In [5], the exact solution of a two-dimensional definition problem in an ideally elastic-plastic cylindrical rod with a given uniform internal temperature is given. The dependence of voltage on temperature in a rod of limited length is given in [6]. In this work, the temporary factor is also taken into account. In [7], different statements of problems of thermal elasticity are given for structural elements under temperature effects. In this case, the theoretical basis of the method is focused on the Maysel formula of uncoupled thermoelasticity. In [8], on the basis of the small parameter method, the problem of determining the stress-strain state of an elastic-plastic pipe in the presence of temperature is considered. It uses the terms of Mises. In [9], the fundamentals of the theory of thermoelasticity and methods for their implementation in solving specific applied problems are presented. This takes into account the power and temperature factors. In the works [10, 11] methods and computational algorithms for numerical solution of the class of applied problems of mechanics are presented.

In contrast to the above works, this paper uses the fundamental energy conservation laws in combination with the constructed quadratic spline functions to solve a particular applied problem.

#### 2. Formulation of the Problem and Methods

We consider a horizontal bar of limited length L [cm], and a constant cross section . The horizontal axis* ox* we will direct from the left to right coincides with the axis of the rod. The side surfaces of the sections () and the core are fully thermally insulated. A heat flux of constant intensity is supplied on the lateral surface of the rod section . Through the cross-sectional area of the left and right ends of the rod, heat exchange with the environment takes place.

Heat transfer coefficients at* x = 0* and at* x = L *. The temperature of the environments of these areas are and [], respectively. The design scheme of the problem is shown in Figure 1.

#### 3. The Solution to the Problem Using the Energy Conservation Law

The rod under consideration is discretized by elements of length L [cm]. Within the length of one discrete element, we approximate the temperature field by a second-order polynomialwhere a, b, c are constants, the values of which are still unknown. The law of temperature distribution within the length of one discrete element is shown in Figure 2.

In the local coordinate system , we fix three nodes with coordinates .

The temperature values in these sections will be denoted, respectively, by the following.

Then substituting (2) into (1), we get the following system.

Solving the system, we define the following.

Substituting (4) into (1), we have the following.

Here we introduce the following notations.

These functions are called quadratic spline functions in the local coordinate system. They have the following properties.

The temperature gradient within the length of one discrete element in the local coordinate system has the following form.

Here the following should be noted.

To construct resolving systems of equations taking into account the natural boundary conditions, we discretize the rod under study by three elemental lengths. For the first discrete element, the total energy functional with the thermal insulation of the side surface is as follows.

Here the first integral over the cross-sectional area of the left end takes place only for points of this surface. In the local coordinate system, taking into account the global node indexing, (11) can be rewritten in the following form:where is the heat transfer coefficient.

Here it must be said that the lateral surface of the first discrete element is fully thermally insulated. In expression (12), it should be noted that the sum of the coefficients before the nodal values of temperature will be zero. For example, in the first bracket (1-1) = 0, and in the second bracket 7-16 + 2-16 + 16 + 7 = 0. Now in Figure 1 we consider the second discrete element. This item is internal. But on the side surface of this element summed thermal constant intensity , the total heat energy functional will have the following form.

Here is the volume of the second discrete element;* P = 2 πk* is the perimeter of the cross section. Finally go to the last third discrete element. The side surface of this element is fully insulated, but through the cross-sectional area of the right end there is a heat exchange with its environment. In this case, the heat transfer coefficient is , and the ambient temperature is . The length of this item is . The total heat energy functional for the third discrete element will be as follows:where is the volume of the third discrete element and

*F*is the cross-sectional area of the rod; then the total heat energy functional for the rod under study has the following form.

To construct resolving systems of linear algebraic equations for the nodal values of temperatures, the functional is minimized by* T1, T2, …, T7*.

Solving this system, the node values of temperatures are calculated. They are used to construct the law of temperature distribution along the length of three discrete sections of the rod.

If the coefficient of thermal expansion of the material of the rod is a constant value, then the value of thermal elongation of the rod in case of pinching with one end is determined in accordance with the theory of thermal physics [11].

If both ends of the rod are clamped, then it cannot lengthen, but there is an axial compressive force R [kG], which is determined by the conditions of compatibility deformation.

Then substituting (19) in (18), we get the following.If both ends of the rod are clamped, then it cannot lengthen, but there is an axial compressive force R [kG], which is determined by the conditions of compatibility deformation. The essence of this approach is as follows. First, consider the horizontal rod clamped by the left end (Figure 3).

This rod is under the influence of the compressive force R which is applied at the right free end. Then, according to Hooke’s law, it is reduced by the value of where is the total length of the test rod, modulus of elasticity of the material of the rod, and rod area.

If both ends of the rod under study are rigidly clamped, then naturally it does not lengthen or shorten; i.e.,

Then, taking into account (21), we will determine the value of the axial compressive force R that arises in the clamped rod by the two ends of the rod under study.

In this case, the distribution field of the thermoelastic component of stress also arises, which is determined in accordance with Hooke’s law [1].

In this case, the distribution field of the thermoelastic component of stress , which is determined in accordance with Hooke’s law [1], also appears.

Then, according to Hooke’s law, the field distribution of the thermoelastic component of deformation has the following form.

The field distribution of the temperature component of the strain and stress is determined on the basis of the general theory of thermoelasticity.

From these relations, in accordance with the theory of thermoelasticity, the field distribution of elastic components of deformations and stresses is determined.The displacement field along the length of one discrete element of length l [cm] is approximated by a complete second-order polynomial.

The gradient of displacement is determined from herewhere . The potential energy functional of elastic deformation in the presence of a temperature field has the following form [1].where is the elastic component of deformation; is the elastic component of stress.

Given these relations, the expression of the potential energy of other strains for the rod under study can be rewritten as follows.

Due to the pinch of both ends of the rod, .

By minimizing this functional on the nodal displacements , a resolving system of linear algebraic equations is constructed for the desired quantities, taking into account the simultaneous presence of heterogeneous heat sources and thermal insulation

By solving the latter system, the values are determined. According to them, the displacement field is constructed within three discrete elements.

#### 4. Analysis of the Obtained Results

For the practical application of the above method and algorithm, we take the following initial data:

with these initial data, the resulting solution is shown in Figures 4–7.

For the numerical solution of the problem under consideration, the rod under study was discretized by n = 24—discrete elements of the same length . The side surface of the first and last 8 elements is fully insulated. The heat flux is supplied to the lateral surface of the middle 8 elements.

From Figure 4, it is seen that, due to the symmetry of the problem under consideration, at the ends of the rod the temperature values will be T (x = 0) = T (x = L) = 540^{0} K. At that time, the highest temperature occurs in the middle of the rod, It can also be seen from Figure 4 that the temperature distribution field along the length of the rod under study has a parabolic character.

The laws of distribution of the three component deformations are shown in Figure 5. Here, the distribution field of the temperature component of deformation has a compressive character along the entire length of the rod under study. It has a parabolic shape with a bulge down. Its value at the ends of the investigated rod will be At the time in the middle of the rod, it has This shows that in the middle of the rod the temperature component of the deformation will be 2.157 times greater than that at the end of the rod. At that time, the value of the thermoelastic component of deformations will be constant along the length of the rod and The behavior of the elastic component of deformation will be of alternating sign. For example, in the middle part of the rod, where the heat flux is supplied, the elastic component of deformation has a tensile character. Amplitude At that time, the first and last 1/3 of the rod will experience compressive . It should be noted that at the ends of the rod In general, the distribution field has a parabolic character, convex at the top. In the relevant generalized Hooke’s law, the behavior of the stress components will be appropriate for the corresponding component strain (Figure 6).

The values of the thermoelastic component of stress along the entire length of the rod will be constant However, the law of distribution of the temperature component of voltage will have a parabolic character, bulge down, and it will have a compressive character along the entire length of the rod. The values of at the ends of the rod will be . At that time in the middle of the rod . This shows that in the middle of the bar the value of will be 2.157 times more than that at its ends. The elastic component of stress in the first and last 1/3 of the rod behaves as compressive, and in the middle 1/3 as tensile. With this

The distribution law of the displacement of the rod sections is not shown in Figure 7. From this figure, it is clear that all sections of the left half of the rod are moved against the direction of the axis* ox*, and sections of the right half are shifted in the direction of the axis of the oh. This process is due to the presence on the surface of the middle 1/3 of the investigated rod summed heat flux. At that time, the amplitude of the displacement of the section whose coordinate x =7.5 cm will be equal to U=(x=7.5)= - 0.01627 cm. Similarly, we have that U=(x=22.5)=0.01627 cm. This shows that the process under study is strictly symmetric about the middle of the rod.

#### 5. Conclusion

The developed method based on the fundamental laws of energy conservation, a computational algorithm, and an application program in Python allows automating the construction of resolving systems of equations taking into account the natural boundary conditions for rods of limited length under the influence of dissimilar types of heat sources. The developed system also allows determining the laws of temperature distribution, all components of deformations and stresses, and displacements. In this regard, it can be stated that the developed method, algorithm, and Python programs are relatively universal in the sense of studying the steady thermo-stress-strain state of the supporting core elements of strategic equipment under the influence of various local heat sources. In this case, the obtained numerical results will differ in high accuracy, since these results satisfy the fundamental laws of energy conservation.

#### Data Availability

The data used to support the findings of this paper are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### References

- S. Timoshenko and J. N. Goodyear,
*Theory of Elasticity*, McGRAW-Hill. Book. Company. Inc., 1987. - K. S. Shibib, M. A. Minshid, and N. E. Alattar, “Thermal and stress analysis in Nd:YAG laser rod with different double end pumping methods,”
*Thermal Science*, vol. 15, suppl. 2, pp. S399–S407, 2011. View at: Publisher Site | Google Scholar - K. S. Shibib, M. M. Tahir, and H. I. Qatta, “Analytical model of transient temperature and thermal stress in continuous wave double-end-pumped laser rod: thermal stress minimization study,”
*Pramana—Journal of Physics*, vol. 79, no. 2, pp. 287–297, 2012. View at: Publisher Site | Google Scholar - D. L. Logan,
*A First Course in the Finite Element Method*, CENGAGE Learning, 2012. - Y. Orçan, “Thermal stresses in a heat generating elastic-plastic cylinder with free ends,”
*International Journal of Engineering Science*, vol. 32, no. 6, pp. 883–898, 1994. View at: Publisher Site | Google Scholar - V. I. Andreev and R. A. Turusov, “Nonlinear modeling of the kinetics of thermal stresses in polymer rods,” in
*Advanced Materials and Structural Engineering*, Hu, Ed., Tailor & Francis Group, London, UK, 2016. View at: Google Scholar - Y. I. Nyashin, V. Y. Kiryukhin, and F. Ziegler, “Control of thermal stress and strain,”
*Journal of Thermal Stresses*, vol. 23, no. 4, pp. 309–326, 2000. View at: Publisher Site | Google Scholar - K. K. Gornostaev, A. V. Kovalev, and Y. V. Malygina, “Stress–strain state in an elastoplastic pipe taking into account the temperature and compressibility of the material,”
*Journal of Physics: Conference Series*, vol. 973, 2018. View at: Google Scholar - B. F. Shorr, “Thermal integrity in mechanics and engineering,” in
*Foundations of Thermoelasticity*, pp. 33–55, Springer, Berlin, Germany, 2015. View at: Google Scholar - O. C. Zienkiewicz and K. Morgan,
*Finite Elements & Approximation*, John Wilei & Sons, Inc., New York, NY, USA, 1983. - O. C. Zienkiewicz and R. L. Taylor,
*The Finite Element Method*, Butterworth-Heinemann, Oxford, UK, 5th edition, 2000.

#### Copyright

Copyright © 2019 Anarbay Kudaykulov 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.