# Thermoelastic Characteristics of Functionally Graded Circular Disk Models under the Loading of Contact Forces

Authors: Jaegwi Go
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

A rotating functionally graded circular disk undergoing a contact load is taken into account to investigate the thermoelastic characteristics. By considering contact force, a pair of partial differential equations is induced as the governing equations based on Hooke’s law. The behavior of circular disk modes is described with the variations of contact force and homogeneous thickness. A finite volume model is introduced to obtain approximate solutions for the governing equations because of the complexity of the equations. Contact force is highly influential in the radial direction compared to the circumferential direction in the displacement distribution, while a large radial stress appears near the area of the contact point. In the strain distribution, the magnitude increases as the angle grows near part of the outer boundary in the circular domain. The radial distribution profiles are susceptible to the growth of contact force in nearby area of the outer boundary, whereas the influences on the circumferential direction profiles are trivial. The increase of homogeneous thickness dwindles the radial magnitude of displacement, stress, and strain distribution profiles over nearby area of the outer boundary of the circular domain. As a result, numerical approach demonstrates that contact force and homogeneous thickness are indispensable parameters and provide deep influence on the thermoelastic movements of a rotating circular disk. Thus, the results obtained may be useful to design an appropriate FGM circular disk model for the industrial area by controlling the above parameters.


---

## Body

## 1. Introduction

The analysis of thermoelastic characteristics of a rotating circular disk has been a challenging work due to complex processes generating during operations. Friction, wear, heat generation, and temperature deformation are representative processes developed during the load of a contact force. By the operation of a rotating circular disk, all parameters are connected mutually and complicated interactions progress. The thermoelastic characteristics of a circular disk are under complicated interaction influences and the variation of each parameter may not be ignorable to determine the thermoelastic movements. Thus, the influences of temperature, contact pressure, and blushed wear on the rotating circular disk are interesting and valuable for research.Ceramic is a typical material used at one surface to resist severe environmental effects such as high temperature, wear, and friction. The other surface is composed of a metal to ensure higher toughness and thermal conductivity. The complementary functions of each material invoke the introduction of functionally graded materials (FGMs). FGMs are nonhomogeneous composites composed of two or more different material phases and the volume fractions continuously vary with space variables. Thus, FGMs are structures involving successively varying mechanical and thermal properties from one surface to another surface. FGMs were introduced to prosper thermal barrier coatings for propulsion system of space planes to endure high temperature and ensure high thermal conductivity at the beginning [1]. The graded material property provides improved bonding strength in various applications undergoing a high temperature gradient field.On the pursuit superior properties such as strength, impact toughness, and fiber resistance, fiber metal laminates (FMLs) have drawn attention as promising hybrid composite materials in last few decades. FMLs are composed of metallic and composite materials. The combination of fiber-reinforced composite layers and composite materials creates FMLs, which exhibit excellent durability. Thus, FMLs are applicable in various industries and are used in aircraft, wind turbine, and automobile. A laminated composite made with thin sheet/grass was investigated by BnvsGaneshGupta et al. [2] for the flexural properties. They evaluated flexural strength and failure analysis using stereo microscope. Zimmermann and Wang [3] analyzed failure modes and characteristics of adhesively and mechanically fastened joints in aircraft composite materials.Meanwhile, the shape of circular disk is commonly employed for practical use in various industrial applications, and the physical behavior of circular disk models is a matter of interest to many engineering researchers. Obata and Noda [4] and Liew et al. [5] analyzed the thermal stresses of FGM circular cylinder and hollow sphere. Afsar et al. [6] investigated the brittle fracture characteristics of thick-walled FGM cylinders by considering incompatible eigenstrain, which is a cause of thermal residual stress. The characterization that oscillatory movements reduce cutting forces in operating process of a circular disk cutter was presented by Kovalyshen [7], and Huang and Li [8] studied the radial component effects on the bending of thin circular plates using Kirchhoff plate theory. The angular speed of a rotating circular disk varies during the start and stop process in machine operation. Dai et al. [9] investigated the displacement and stress fields of a FGM hollow circular disk subjecting an angular acceleration due to a changing temperature. Static behavior of functionally graded auxetic-porous structures was presented by Rad [10] considering torsional interaction and horizontal friction force. The static response of nonuniform bidirectional functionally graded auxetic-porous material circular plates is carried out based on poroelasticity theory. Friction induced a vibration in rotary contact system, and the vibration causes a misalignment problem. Tadokoro et al. [11] described the stabilizing effect arising from parallel misalignment in circular sliding contact. Infinitesimal theory of plane elasticity and Complementary Functions Method was employed to investigate hydrogen-induced stresses in functionally graded axisymmetric spheres, cylinders, and disks by Yildirim [12]. Sondhi et al. [13] used element-based gradation of material properties over the discretized domain to study linear elastic stress analysis of functionally graded rotating disks.However, the comprehensive understanding and quantifying for the thermoelastic field of a rotating FGM circular disk may be useful in design FGM cutter or a grinding disk. In general, FGM cutter or a grinding disk experiences the thermal load and inertia force due to rotation of the disk. The previous study [14] demonstrated that proper control of some parameters such as temperature distribution, angular speed, and radial thickness improves thermoselastic characteristics in the design of a circular or grinding disk and reduces failure mechanism such as crack and brittle, but most of previous works progressed their studies without considering contact load even though contact load is a crucial parameter in determination of thermoelastic characteristics. The present study focuses on the analysis of thermoelastic characteristics of a rotating Al2O3/Al FGM disk subjected to a contact force. Young’s modulus, CTE, and density of the FGM disk are assumed to vary exponentially only in the radial direction due to symmetry with respect to the axis of the disk, whereas Poisson’s ratio is assumed to be constant because of insignificant effect on the thermoelastic characteristics of the disk. The research addresses the following: (i) the contact force and homogeneous thickness are crucial factors in the determination of the thermoelastic characteristics of FGM circular models, (ii) the radial direction is more susceptible to the variation of the magnitude of contact force and homogeneous thickness, and (iii) the variation of homogeneous thickness is highly influential to thermoelastic characteristics over the area of 0.2<r−a/b−a<0.8 in FGM circular domain.The process for the mathematical approach is presented in Section2. A pair of partial differential equations is derived based on Hooke’s law in Section 2.1. Due to the complexity of the governing equation, a finite volume method is introduced in Section 2.2 to obtain numerical solutions. In Section 3, numerical solutions are displayed and the elastic characteristics are discussed according to the results obtained. The significances of the results are addressed in the conclusions in Section 4.


## 2. Mathematical Modeling

A rotating FGM circular disk with a concentric circular hole is considered (see Figure1). The origin of the polar coordinate system r−θ is assumed to be located at the center of the disk and hole. A and B, denoted by the dark and the white colors, represent constituent materials of FGM circular disk, as shown in the figure. The distribution of each material varies continuously along the radial direction only. The radii of the hole and outer surface of the disk are designated by a and b. The angular velocity ω can be determined from the relation ω=2πN/60, N represents the revolutions per minute (rpm).Since the material distributions vary along the radial direction only, the FGM disk can be reduced to an axisymmetric problem. Thus, all properties of the disk can be treated as a function of r only. Young’s modulus, coefficient of thermal expansion, and the density of the disk are denoted by E, α, and ρ, respectively, and each property is assumed to vary exponentially as follows (for various and detailed distribution profiles, see [15]):(1a)E=E0eβr,(1b)α=α0eξr,(1c)ρ=ρ0eμr.Figure 1
(a) Schematic diagram of functionally graded circular disk models. (b) Notations of finite control volumes.
(a)(b)By the assumption that the disk is composed of 100% materialA at the surface of the hole (r=a) and 100% material B at the outer surface (r=b), the constants in equations (1a)–(1c) can be determined as follows:(2a)E0=EAe−βa,(2b)α0=αAe−ξa,(2c)ρ0=ρAe−μa,(3a)β=1a−blnEAEB,(3b)ξ=1a−blnαAαB,(3c)μ=1a−blnρAρB.The subscriptsA and B represent the properties of the constituent materials A and B, respectively, but the nonsubscripted variables are used to denote the properties of FGM composed of materials A and B.


### 2.1. Thermoelastic Formulation

In the polar coordinate, the relations between the strain components and deformation components are as follows:(4)εr=∂u∂r,εθ=1r∂v∂θ+ur,εz=0,γrθ=1r∂u∂θ+∂v∂r−vr,γzθ=0,γrz=0.The plane stress exposed with thermal expansion, by Hooke’s law, produces the strain-stress relations undergoing thermal expansion ([16], p8):(5)εr=1Eσr−υσθ+αT,εθ=1Eσθ−υσr+αT,τrθ=E21+νγrθ,τθz=0,τrz=0.Tr is the function of temperature variation with radius variable r. The temperature function can be expressed as(6)Tr=c1lnr+c2,under the assumption that the circular disk is subjected by a loading of symmetric temperature to the radial direction only [17]. c1and c2 will be determined based on the boundary conditions. The equilibrium equations in polar coordinates are(7a)∂σr∂r+1r∂τrθ∂θ+σr−σθr+ρω2r=0,(7b)∂σθ∂θ+r∂τrθ∂r+2τrθ=0,and the algebra among equations (4), (5), (7a), and (7b) yields the following governing equations:(8a)∂∂rr∂u∂r+1−υ2r∂2u∂θ2−3−υ2r∂v∂θ+1+υ2∂2v∂r∂θ−ur+1−υ2Eρω2r2=υ+1αrdTrdr,(8b)1−υ2∂∂rr∂v∂r+1r∂2v∂θ2+3−υ2r∂u∂θ+1+υ2∂2u∂r∂θ−1−υ2vr=0.The governing equations (8a) and (8b) are a pair of partial differential equations in the polar coordinate system expressed with the radial and circumferential displacement distribution profiles under the loading of temperature. The thermoelastic characteristics of a functionally graded rotating circular disk subjected to a contact force are investigated according to the boundary conditions:(9)σra,θ=0,σrb,θ−0=0,σrb,0=P,at a contact point.The boundary condition implies that the circular model is under the loading of thermal expansion over the entire domain and of radial pressureP at the outer boundary contact point b,0.


### 2.2. Finite Volume Formulation

Due to the complexity of the governing equations, a numerical technique is required, and a finite volume method is applied for approximation. The domain is divided up into control volume and integrates the field equations over each control volume (see Figure1(b)). Equations (8a) and (8b) thus are discretized as follows [18]:(10)A11ui+1,j+A12ui,j+1+A13ui,j+A14ui,j−1+A15ui−1,j+B11vi,j+1+B12vi,j+B13vi,j−1+B14vi−1,j+1+B15vi−1,j+B16vi−1,j−1+B17vi−2,j+1+B18vi−2,j+B19vi−2,j−1=fi,j,A21ui,j+1+A22ui,j+A23ui,j−1+A24ui−1,j+1+A25ui−1,j+A26ui−1,j−1+A27ui−2,j+1+A28ui−2,j+A29ui−2,j−1+B21vi+1,j+B22vi,j+1+B23vi,j+B24vi,j−1+B25vi−1,j=0,(11)A11=1−υ2ΔrΔθ1ri,j,A12=ΔθΔrri,j+1/2,A13=−ΔθΔrr+i,j+1/2ri,j−1/2−1−υΔrΔθ1ri,j−ΔrΔθ1ri,j,A14=ΔθΔrri,j−1/2,A15=1−υ2ΔrΔθ1ri,j,B11=381+υ,B12=−343−υΔr1ri,j,B13=−7161+υ,B14=−121+υ,B15=523−υ1ri,j,B16=9161+υ,B17=181+υ,B18=−343−υΔr1ri,j,B19=−181+υ,fi,j=1+υαri,jΔθTi,j+1/2−Ti,j−1/2−ΔrΔθρω21−υ2Eri,j2,A21=381+υ,A22=343−υΔr1ri,j,A23=−7161+υ,A24=−121+υ,A25=−3−υΔr1ri,j,A26=9161+υ,A27=181+υ,A28=3−υΔr1ri,j,A29=−181+υ,B21=ΔrΔθ1ri,j,B22=ΔθΔrri,j+1/2,B23=−ΔθΔrr+i,j+1/2ri,j−1/2−1−υ2ΔrΔθ1ri,j−2ΔrΔθ1ri,j,B24=ΔθΔrri,j−1/2,B25=ΔrΔθ1ri,j.The validation and convergence of the above finite volume method (FVM) were presented in [17], and Mathematica (version 5.0) is used to obtain numerical solutions.


## 2.1. Thermoelastic Formulation

In the polar coordinate, the relations between the strain components and deformation components are as follows:(4)εr=∂u∂r,εθ=1r∂v∂θ+ur,εz=0,γrθ=1r∂u∂θ+∂v∂r−vr,γzθ=0,γrz=0.The plane stress exposed with thermal expansion, by Hooke’s law, produces the strain-stress relations undergoing thermal expansion ([16], p8):(5)εr=1Eσr−υσθ+αT,εθ=1Eσθ−υσr+αT,τrθ=E21+νγrθ,τθz=0,τrz=0.Tr is the function of temperature variation with radius variable r. The temperature function can be expressed as(6)Tr=c1lnr+c2,under the assumption that the circular disk is subjected by a loading of symmetric temperature to the radial direction only [17]. c1and c2 will be determined based on the boundary conditions. The equilibrium equations in polar coordinates are(7a)∂σr∂r+1r∂τrθ∂θ+σr−σθr+ρω2r=0,(7b)∂σθ∂θ+r∂τrθ∂r+2τrθ=0,and the algebra among equations (4), (5), (7a), and (7b) yields the following governing equations:(8a)∂∂rr∂u∂r+1−υ2r∂2u∂θ2−3−υ2r∂v∂θ+1+υ2∂2v∂r∂θ−ur+1−υ2Eρω2r2=υ+1αrdTrdr,(8b)1−υ2∂∂rr∂v∂r+1r∂2v∂θ2+3−υ2r∂u∂θ+1+υ2∂2u∂r∂θ−1−υ2vr=0.The governing equations (8a) and (8b) are a pair of partial differential equations in the polar coordinate system expressed with the radial and circumferential displacement distribution profiles under the loading of temperature. The thermoelastic characteristics of a functionally graded rotating circular disk subjected to a contact force are investigated according to the boundary conditions:(9)σra,θ=0,σrb,θ−0=0,σrb,0=P,at a contact point.The boundary condition implies that the circular model is under the loading of thermal expansion over the entire domain and of radial pressureP at the outer boundary contact point b,0.


## 2.2. Finite Volume Formulation

Due to the complexity of the governing equations, a numerical technique is required, and a finite volume method is applied for approximation. The domain is divided up into control volume and integrates the field equations over each control volume (see Figure1(b)). Equations (8a) and (8b) thus are discretized as follows [18]:(10)A11ui+1,j+A12ui,j+1+A13ui,j+A14ui,j−1+A15ui−1,j+B11vi,j+1+B12vi,j+B13vi,j−1+B14vi−1,j+1+B15vi−1,j+B16vi−1,j−1+B17vi−2,j+1+B18vi−2,j+B19vi−2,j−1=fi,j,A21ui,j+1+A22ui,j+A23ui,j−1+A24ui−1,j+1+A25ui−1,j+A26ui−1,j−1+A27ui−2,j+1+A28ui−2,j+A29ui−2,j−1+B21vi+1,j+B22vi,j+1+B23vi,j+B24vi,j−1+B25vi−1,j=0,(11)A11=1−υ2ΔrΔθ1ri,j,A12=ΔθΔrri,j+1/2,A13=−ΔθΔrr+i,j+1/2ri,j−1/2−1−υΔrΔθ1ri,j−ΔrΔθ1ri,j,A14=ΔθΔrri,j−1/2,A15=1−υ2ΔrΔθ1ri,j,B11=381+υ,B12=−343−υΔr1ri,j,B13=−7161+υ,B14=−121+υ,B15=523−υ1ri,j,B16=9161+υ,B17=181+υ,B18=−343−υΔr1ri,j,B19=−181+υ,fi,j=1+υαri,jΔθTi,j+1/2−Ti,j−1/2−ΔrΔθρω21−υ2Eri,j2,A21=381+υ,A22=343−υΔr1ri,j,A23=−7161+υ,A24=−121+υ,A25=−3−υΔr1ri,j,A26=9161+υ,A27=181+υ,A28=3−υΔr1ri,j,A29=−181+υ,B21=ΔrΔθ1ri,j,B22=ΔθΔrri,j+1/2,B23=−ΔθΔrr+i,j+1/2ri,j−1/2−1−υ2ΔrΔθ1ri,j−2ΔrΔθ1ri,j,B24=ΔθΔrri,j−1/2,B25=ΔrΔθ1ri,j.The validation and convergence of the above finite volume method (FVM) were presented in [17], and Mathematica (version 5.0) is used to obtain numerical solutions.


## 3. Numerical Results and Discussion

In this section, the differential equation induced in Section2.1 and the finite volume formula developed in Section 2.2 are used to evaluate the approximated results of different components of displacement, stress, and strain for an Al2O3/Al FGM circular disk. The mechanical and thermal properties of these ingredient materials are shown in Table 1.Table 1
Mechanical and thermal properties used for analyzing thermoelastic characteristics of rotating CM circular disks.
Material/propertyElastic modulus (MPa)Thermal expansion coefficient (10−6/°C)Density (g/cm3)Substrate (Al)7123.12.70Bond coating (Al/Al2O3)164.313.61.61Top (Al2O3)3808.00.96Stress distribution profiles are presented in Figure2. As shown in Figure 2(a), the circular disk is undergoing tensile radial stress. All areas of the disk except near outer boundary experience trivial influence in radial stress distribution, whereas the contact point suffers from the largest radial stress. But different phase develops in circumferential stress. At normalized values r−a/b−a=0.1 and 0.5, the magnitude of circumferential stresses is ignorable in comparison with circumferential stress at r−a/b−a=0.9 (see Figure 2(b)). In addition, the entire area of the disk is under the loading of tensile circumferential stress at r−a/b−a=0.9 and the circumferential stress distribution exhibits exponential decay with the increase of angle. The strain distribution profiles are described in Figure 3. Near area of outer boundary in circular domain presents that the magnitude increases with the increase of angle value in the radial strain distribution (see Figure 3(a)). Positive radial strain distribution appears over the disk, and the largest radial strain value occurring at the contact point generates. However, the circumferential strain distribution shows different pattern (see Figure 3(b)). At r−a/b−a=0.9, the circumferential strain deploys to the positive direction, while negative circumferential strain distributions are obtained at r−a/b−a=0.1 and 0.5. The changing shape of the circumferential strain is a parabola as the angle increases at r−a/b−a=0.1 and 0.5, whereas the circumferential strain distribution profile fluctuates at the beginning and logarithmic decay appears after θ=2 at r−a/b−a=0.9. The results inform that (i) contact force is highly influential to near area of contact point in radial direction, (ii) the circumference stress is more sensitive to contact force than the radial stress distribution, and (iii) various phenomena appear in circumferential strain distribution depending on the value of r−a/b−a.Figure 2
Stress distribution profiles in the (a) radial direction and (b) circumferential direction.
(a)(b)Figure 3
Strain distribution profiles in the (a) radial direction and (b) circumferential direction.
(a)(b)Figure4 presents the influences of contact force on the stress distribution by choosing the representative angles θ=0 and 180 and the normalized radius r−a/b−a=0.1 and 0.9. The radial stress distribution is susceptible to the change of contact force, especially near the area of contact point (see Figures 4(a) and 4(b)). Entire domain of the circular disk suffers from a tensile stress and the magnitude of the radial stress is getting larger as the contact force increases. As shown in Figure 4(c), the area near the inner boundary of the disk sits loose to the variation of contact force, but the area near the outer boundary reacts sensitively to the change of contact force, and noticeable impact appears when θ=0 (see Figure 4(d)). By the increase of contact force, the radial and circumferential stresses exhibit exponential growth near area of contact point, even though minor effects appear over entire circular domain.Figure 4
Effects of contact force on the stress in the radial direction (a) atθ=0 and (b) at θ=180 and in the circumferential direction (c) at r−a/b−a=0.1 and (d) at r−a/b−a=0.9.
(a)(b)(c)(d)The effects of contact force on the strain distributions are described in Figure5. The reactions of the radial strain distribution to the change of contact force are similar to those of radial stress distribution (see Figures 5(a) and 5(b)). Nearby area of the outer boundary is particularly susceptible to the increase of contact force and the magnitude of the radial strain grows with the increment of contact force. Trivial reaction appears in the circumferential strain distribution to the variation of contact force for r−a/b−a=0.1 (see Figure 5(c)), whereas nearby part of the outer boundary is sensitive to the change of contact force (see Figure 5(d)). Pronounced alterations in the circumferential strain distribution appear on area around θ=0. The numerical solutions imply that (i) the magnitude of contact force is a critical factor to determine the thermoelastic characteristics of circular FGM models and (ii) the radial direction is more susceptible to the variation of the magnitude of contact force.Figure 5
Effects of contact force on the strain in the radial direction (a) atθ=0 and (b) at θ=180 and in the circumferential direction (c) at r−a/b−a=0.1 and (d) at r−a/b−a=0.9.
(a)(b)(c)(d)The influences of homogeneous thickness on thermoelastic characteristics, denoted byH, are explained through Figures 6–8. H=1 implies that the circular disk is composed of homogeneous material. The representative values H=0.25, 0.5, and 0.75 are applied. Figure 6 expresses the effects of H on the displacement distributions.Figure 6
Effects of homogeneous thickness on the displacement in the radial direction (a) atθ=0 and (b) at θ=180 and in the circumferential direction (c) at r−a/b−a=0.1 and (d) at r−a/b−a=0.9.
(a)(b)(c)(d)Figure 7
Effects of homogeneous thickness on the stress in the radial direction (a) atθ=0 and (b) at θ=180 and in the circumferential direction (c) at r−a/b−a=0.1 and (d) at r−a/b−a=0.9.
(a)(b)(c)(d)Figure 8
Effects of homogeneous thickness on the strain in the radial direction (a) atθ=0 and (b) at θ=180 and in the circumferential direction (c) at r−a/b−a=0.1 and (d) at r−a/b−a=0.9.
(a)(b)(c)(d)The area between normalized valuesr−a/b−a=0.2 and r−a/b−a=0.85 reveals noticeable reactions to homogeneous thickness variation. The magnitude of the radial displacement increases with the growth of homogeneous part (see Figures 6(a) and 6(b)), but nearby area of the outer boundary displays that the magnitude to the positive direction decreases as the value of H increases. The results imply that the growth of functionally graded area in the circular disk yields more stability in the radial displacement distribution profiles.As shown in Figures6(c) and 6(d)), minor influence appears in the circumferential displacement except nearby area of contact force. The circumferential displacement increases with growth of the value of H over the area around contact force.Figure7 presents the effects of H on the stress distributions. The radial distribution exhibits smaller tensile stress over the nearby area of outer boundary with increase of the value of H when θ=0 (see Figure 7(a)). But the circular disk undergoes larger compressive radial stress over area between normalized values r−a/b−a=0.1 and r−a/b−a=0.8 when θ=180 (see Figure 7(b)). According to the radial stress distributions, the circular disk with larger functionally graded area may not guarantee higher stability over nearby area of outer boundary. Different phase develops in circumferential stress distributions. Over nearby area of the inner boundary, entire domain of the disk experiences compressive circumferential stress for the values of H=0.5 and 0.75, whereas some part of the disk with H=0.25 undergoes tensile circumferential stress (see Figure 7(c)). Over nearby area of the outer boundary of the disk, the circumferential stress with H=0.75 is trivial in comparison with the value of H=0.25 and 0.5 (see Figure 7(d)). The magnitude of the circumferential stress is higher near contact point as the homogeneous thickness decreases. Based on the circumferential stress distributions, greater functionally graded area affords higher stability in the circular disk. The influences of homogeneous thickness H on the strain distributions are described in Figure 8. Nearby area of the outer boundary of the disk is sensitive to the change of homogeneous thickness in the radial strain distributions. The magnitude of strain decreases with the growth of the value of H (see Figures 8(a) and 8(b)), which signifies that the circular disk with larger homogeneous area is more stable over nearby area of the outer boundary. Similar movements occur in the circumferential strain. As shown in Figures 8(c) and 8(d), the magnitude of the circumferential strain is getting larger according to the decline of the value of H. The outcomes describe that the circular disk stability with greater functionally graded area is greater in the circumferential strain distribution. The numerical approaches make point of the following: (i) the variation of homogeneous thickness is highly influential to thermoelastic characteristics over the area of 0.2<r−a/b−a<0.8 and (ii) the radial direction is more susceptible to the change of homogeneous thickness.


## 4. Conclusions

A rotating FGM circular disk with a concentric circular hole is taken into account to investigate the thermoelastic characteristics. The circular disk is under the loading of a contact force. The effects of the variation of contact force and homogeneous thickness on the components of displacement, stress, and strain have been described. Over nearby area of the contact point, the thermoselastic behavior reacts sensitively to the change of contact force and exhibits exponential growth in the magnitude, whereas the influence of contact force over the other area is greater in the circumferential stress distributions. The outcomes explain that the outer boundary area is exposed to become brittle with the growth of a contact force. In the radial displacement distribution, the growth of functionally graded area grants more stability in the circular disk, while, according to the radial stress distributions, the circular disk with larger functionally graded area may not guarantee higher stability over nearby area of outer boundary. Based on the circumferential distribution profiles, the circular disk with greater functionally area evinces higher stability. Overall, the results demonstrate that contact force and homogeneous thickness are crucial parameters to determine the thermoelastic characteristics of FGM circular disks, and the movements of an FGM circular disk can be controlled by controlling of these parameters. Therefore, a FGM circular cutter or grinding disk suffering from the loading of a contact force can be designed applying the finite volume analysis used in this investigation to promote proper and reliable thermoelastic characteristics in service.



---

*Source: 8845028-2021-01-06.xml*