Transient Thermal Stresses in FG Porous Rotating Truncated Cones Reinforced by Graphene Platelets

: The present work studies an axisymmetric rotating truncated cone made of functionally graded (FG) porous materials reinforced by graphene platelets (GPLs) under a thermal loading. The problem is tackled theoretically based on a classical linear thermoelasticity approach. The truncated cone consists of a layered material with a uniform or non-uniform dispersion of GPLs in a metal matrix with open-cell internal pores, whose effective properties are determined according to the extended rule of mixture and modiﬁed Halpin–Tsai model. A graded ﬁnite element method (FEM) based on Rayleigh–Ritz energy formulation and Crank–Nicolson algorithm is here applied to solve the problem both in time and space domain. The thermo-mechanical response is checked for different porosity distributions (uniform and functionally graded), together with different types of GPL patterns across the cone thickness. A parametric study is performed to analyze the effect of porosity coefﬁcients, weight fractions of GPL, semi-vertex angles of cone, and circular velocity, on the thermal, kinematic, and stress response of the structural member.


Introduction
Recently, multifunctional porous structures have been largely applied in many engineering applications because of their excellent properties, namely, lightweight, heat resistance, and great energy absorption [1,2]. The presence of internal pores within a metal matrix, however, can significantly reduce the structural stiffness [3][4][5], whose limitations can be overcome by adopting lightweight materials with nanofillers e.g., carbon nanotubes (CNTs) [6][7][8] or graphene nano-platelets (GPLs) [9]. In lightweight structures, such nanoparticles keep their excellent potentials, especially for even dispersions within a metal or polymer matrix [10], thus preventing any potential agglomeration of reinforcements. More specifically, GPLs feature higher mechanical properties than CNTs, a lower cost and a more special surface area and 2D geometry [11].
In such a context, FG porous structures reinforced with GPLs have increased the attention of the recent literature to improve the efficiency and capability of structural members, while controlling the density, size, and pattern of porosities within the material along with the GPL distribution [12,13]. From a numerical and analytical standpoint, several researches have investigated the influence of porosity and graphene platelets pattern on the mechanical behavior of structures subjected to different loading and boundary conditions. Chen et al. [13], for example, applied a Timoshenko beam theory and von Kármán type nonlinearity to investigate the nonlinear vibration and post-buckling behavior of FG-GPLs porous nanocomposite beams, whose problem was solved numerically Accordingly, Mohammadjani and Shariyat [51] investigated the nonlinear thermomechanical vibration of FG annular plates/disks subjected to a non-uniform thermomechanical loading condition in a dynamic sense, by applying a 3D thermoviscoelasticity theory. A hygro-thermo-mechanical nonlinear study was also proposed recently for curved thin and moderately thick shallow FG panels in [52].
Deka et al. [53] presented a semi-analytical solution for FG rotating disks by enforcing both Dirichlet and Neumann boundary conditions while implementing a homotopy perturbation method (HPM) for a thermomechanical loading condition. A similar problem was also studied by Talebitooti et al. [54] to check for the critical buckling behavior of stiffened rotating FG cylindrical shells subjected to a thermomechanical loading, as provided by a FSDT and DQM. More recently, a 1D finite element method has been proposed in [55] to investigate the 3D thermoplastic response of rotating disks with variable thickness by employing a Carrera-unified formulation (CUF). Saadatfar and Fakhri [56] studied the response of rotating FG hybrid cylindrical shells with a piezoelectric layer under hygrothermal conditions. Tornabene [57] applied a general theoretical framework to survey the dynamics of rotating doubly-curved shell structures made of FG materials. Khorsand and Tang [58] optimized the weight of a rotating FG annular disk with variable thickness subjected to a thermo-mechanical loading by using co-evolutionary particle swarm optimization and DQMs.
Based on the available literature, in most cases the investigations focus on the steady state and time dependent thermal stresses of cylindrical shells, where a little attention is given to the thermoelastic behavior of conical shells under a transient thermal loading. Starting with the preliminary work [59], the present work provides a more generalized investigation on the transient thermal stress behavior of FG rotating cones in presence of possible porosities. The multilayer cone is assumed to be reinforced with uniform and nonuniform GPLs in a metallic matrix with possible open-cell internal pores. Three different porosity distributions are assumed across the thickness direction of the structure including uniform and two symmetric FG patterns. The parametric study also includes five different dispersion patterns of GPL through the thickness direction. A 2D axisymmetric elasticity theory is herein employed to define the problem which is solved numerically in time and space, based on a classical finite element approach, Rayleigh-Ritz energy formulation and Crank-Nicolson algorithm. A large numerical investigation checks for the sensitivity of the transient response of FG porous truncated cone to different input parameters, including various rotating velocities, porosity distributions, semi-vertex angles of the cone, porosity coefficients, GPL dispersion patterns, and weight fraction of the GPL nanofillers, as useful for design purposes as well as for different mechanical and aerospace applications, such as high-speed centrifugal separators, gas turbines, or high-power aircraft jet engines. The work is organized as follows: after this introduction, the theoretical basics of the problem are defined in Section 2, whose finite element modeling is detailed in Section 3. A large numerical investigation is presented and discussed in Section 4, whose final remarks are reported in the concluding Section 5.

Geometry Description
Let us consider a FG-GPL porous truncated cone with length L, thickness h (along the x direction, 0 ≤ x ≤ h), internal and external radii of the small base a and b, respectively, and semi-vertex angle ϕ, as depicted in Figure 1. Axisymmetric cylindrical coordinates (r, z) are assumed along the radial and axial directions, respectively. Details of porosity distribution and GPL pattern of the structure are also reported in Figure 1.

Material Properties
Three different porosity distributions are assumed throughout the cone thickness ( Figure 1), i.e., two types of non-uniform symmetric distribution and a uniform one. In the first distribution (labeled as I-D1), the porosity is nonlinearly symmetric, with a higher density around the mid-plane rather than the inner and outer surfaces of the cone. In the second distribution (labeled as II-D2), a non-uniform symmetric porosity is assumed, with a higher density near the inner and outer surfaces rather than the mid-plane. For such two cases, the material properties with porosities are defined in Equations (1) and (2), respectively, whereas for a uniform distribution of porosity (labeled as D3) the material properties read as in Equation (3). At the same time, five different GPL patterns are assumed along the cone thickness, as represented in Figure 1 and defined in [43,60]. More specifically: Porosity distribution 1 (Non-uniform symmetric I-D1) Porosity distribution 2 (Non-uniform symmetric II-D2) Figure 1. Geometry of rotating truncated cone for various porosity and GPL patterns across the x axis (thickness direction).

Material Properties
Three different porosity distributions are assumed throughout the cone thickness ( Figure 1), i.e., two types of non-uniform symmetric distribution and a uniform one. In the first distribution (labeled as I-D1), the porosity is nonlinearly symmetric, with a higher density around the mid-plane rather than the inner and outer surfaces of the cone. In the second distribution (labeled as II-D2), a non-uniform symmetric porosity is assumed, with a higher density near the inner and outer surfaces rather than the mid-plane. For such two cases, the material properties with porosities are defined in Equations (1) and (2), respectively, whereas for a uniform distribution of porosity (labeled as D3) the material properties read as in Equation (3). At the same time, five different GPL patterns are assumed along the cone thickness, as represented in Figure 1 and defined in [43,60]. More specifically: Porosity distribution 1 (Non-uniform symmetric I-D1) Porosity distribution 2 (Non-uniform symmetric II-D2) Appl. Sci. 2022, 12, 3932

of 26
Uniform porosity distribution (D3) The following equations are also applied for evaluating x [61] In the previous relations, x is the normal space of points through the thickness of the cone from the inner surface; r n and z n refer to the radius and height of an arbitrary point in the domain; r in and z in refer to the radius and height of points at the inner surface changing from a to (a + L tan ϕ) and from 0 to L, respectively. Moreover, E(x), G(x), ρ(x), c(x), α(x), k(x) stand for the modulus of elasticity, shear stiffness, mass density, specific heat capacity, thermal expansion coefficient, and heat conduction coefficient of the porous nanocomposite cone, respectively. At the same time, E * , G * , ρ * , c * , α * , k * are the similar properties of GPL cone without any porosity. Furthermore, e 0 and e * 0 (0 ≤ e 0 (e * 0 ) < 1) are the coefficients of porosity belonging to D1 and D2, respectively; e m and e * m stand for the coefficients of mass density related to D1 and D2, respectively; Θ and Θ stand for the variables of a uniform porosity distribution. As the size and density of the internal cavities increase, the porosity increases, leading to a reduction in the material properties. The influence of the material properties of open-cell metal foams [62][63][64] is defined mathematically by Equation (5), from which we can derive the Relation (6) between e 0 and e m , as follows The mass of FG-GPL truncated cone is assumed to be identical, such that here employed to determine e * 0 and Θ for a known value of e 0 . It can be observed that e * 0 rises by enhancing e 0 . When e 0 reaches 0.6, e * 0 (=0.9612) is near to the upper bound. This justifies the selection of e 0 ∈ [0, 0.6], as applied hereafter.
Based on a Halpin-Tsai micromechanics model [65], the elasticity modulus for the nanocomposite without porosities E * is defined as where E GPL and E m refer to the elasticity modulus of GPLs and metallic matrix, respectively, l GPL , w GPL , and t GPL refer to the length, width, and thickness of nanofiller platelets, and V GPL is the volume content of GPLs. The rule of mixture [66] is used to compute the mass density, Poisson's ratio, heat capacitance, heat conductivity coefficient, and heat expansion coefficient of the nanocomposite, namely, where ρ GPL , ρ m , ν GPL , ν m , c GPL , c m , k GPL , k m , α GPL , α m stand for the mass density of GPLs, mass density of metal matrix, Poisson's ratio of GPLs, Poisson's ratio of metal matrix, specific heat capacity of GPLs, specific heat capacity of metal matrix, heat conduction coefficient of GPLs, heat conduction coefficient of metal matrix, thermal expansion coefficient of GPLs, thermal expansion coefficient of metal matrix, respectively. Moreover, the Poisson's ratio is kept constant for open-cell metal foams [67]. The shear modulus G * of the nanocomposite reads as follows The volume content of GPLs, V GPL , changes across the cone thickness for various dispersion patterns, as follows where t i1 , t i2 , and t i3 denote the top limit of the V GPL , and i = 1, 2, 3 corresponding to porosity distributions 1, 2, and 3, respectively. V T GPL is estimated by applying the nanofiller weight fraction ∆ GPL into Equation (19), which is used to derive t i1 , t i2 , and t i3 by the following relation

Heat Transfer Equation
The heat conduction equation can be defined in axisymmetric cylindrical coordinates as where the heat generation rate has been neglected, T is the temperature, and k r (x) and k z (x) stand for the heat conduction coefficients in the radial and axial direction, respectively. The thermal boundary conditions are considered as where T 0 refers to the temperature at t = 0, c 0 is a constant value, r out is the radius of points at the external surface that changes from b to b + L tan ϕ.

Thermo-Elasticity Equations
Based on a classical linear thermo-elasticity theory, the strain field for an axisymmetric problem reads as follows and u and v refer to the displacement components across the r and z directions, respectively. Moreover, the kinematic relations read as follows [59] [σ] = where ∆T = T − T 0 is the temperature estimated from the heat conduction equation. In addition, [ε] stands for the elastic strain and [ε T ] denotes the thermal strain field stemming from the temperature variation. The elastic stiffness matrix D d is defined as [68] The FG-GPL porous truncated cone is simply supported on its two bases, such that the following kinematic boundary conditions are enforced

Finite Element Modelling
The problem is solved numerically by using the finite element method. In such a context, we approximate the cone in the r − z plane of revolution by means of shape . Based on a Kantorovich approximation, we define the thermal distribution for each element (e) as follows Appl. Sci. 2022, 12, 3932 The indexes d and t refer to the structural and thermal solutions of the problem, respectively. By substitution of Equation (31) into Equation (26), the elastic strain matrix for the arbitrary element (e) can be defined as where , or more specifically Note that matrix [N(r, z)] (e) refers to the second-order interpolation functions in terms of its nodal values for element (e), and [T(t)] denotes the corresponding nodal temperature vector. In addition to the displacement components, the inhomogeneity of the material properties of the FG-GPL porous cone can be expressed by means of their nodal values as follows where E i , ρ i , c i , α i , k i refer to the elasticity modulus, mass density, specific heat capacity, thermal expansion coefficient, and heat conduction coefficient related to the arbitrary node i. Moreover, N,Ê,ρ,Ĉ,α, andk refer to the shape functions vector, elasticity modulus, mass density, heat capacity, heat expansion coefficient, and heat conductivity coefficient for each element, defined as By a mathematical manipulation, Equation (22) can be simplified as where V is the volume, S 1 and S 2 are the boundary surfaces of the cone, q , Q, T ∞ , and h stand for the heat flux, rate of energy generation per unit volume, surrounding temperature, and heat convection coefficient, respectively.
According to FEM, we redefine Equation (42) (43) [ where in which V(e) denotes the volume element, S 1 (e) and S 2 (e) refer to the boundary elements subjected to a heat flux and heat convection, respectively. By assembling the element matrices, the global heat transfer takes the following form An unconditionally stable scheme based on a Crank-Nicolson time stepping procedure is thus applied for numerically solving the differential Equation (49), under a suitable time step selection, see Ref. [70]. More details about the algorithm are provided in Appendix A. Once deriving the transient temperature distribution, we solve the thermoelastic problem, while using simplex linear triangular elements to approximate the kinematic field. Based on a Rayleigh-Ritz energy formulation, the structural stiffness and force element matrices owing to a varying temperature and rotational velocity are, thus, obtained as The first term in (51) is the force due to a thermal variation and the second term is that one due to a rotational velocity; moreover, Γ r defines the body force component in the r direction Γ r = ρrω 2 (53) whereas the governing equations of the problem take the following matrix form This equation must be solved at each time step with the known temperature variation as derived from Equation (49).

Verification
In this section we present the results from a transient stress analysis for a simplysupported FG-GPL porous rotating truncated cone, defined geometrically by a = 1, b = 1.5, L = 2, ϕ = 15 • , and ω = 100 rad/s. The selected structural member is made of a homogenous material with the following thermo-mechanical properties: E = 198.8 GPa, ρ = 8340 kg/m 3 , ν = 0.33, k = 512 W/(m · k), C = 876 j/(Kg · k), α = 7.74 × 10 −5 1/k, while the thermal boundary conditions are enforced for T 0 = 300 • k, T 1 = 600 • k. Our numerical results are compared to those ones from the commercial FEM software ANSYS Workbench, in terms of thermal time history, radial displacement, radial stress state at the midpoint of the cone, under the assumptions γ GPL = 0% and e 0 = 0, e m = 0 (see Figure 2). To model the problem in ANSYS Workbench, quadratic quadrilateral elements have been used for the sake of accuracy, under the assumption of 2D axisymmetric elasticity like in the present study. To get convergent results, the time step has been set equal to 1 s, and 20 × 30 elements have been considered along the radial and axial directions. As visible from all the plots in Figure 2, our results match perfectly with predictions from the ANSYS code, thus confirming the reliability and accuracy of the proposed FEM-based formulation to treat the selected thermomechanical problem. All the time histories in Figure 2 refer to the center point of the cone, whose temperature ( Figure 2a) and radial displacement ( Figure 2b) feature a smooth and monotonic increase up to a plateau value, while observing a non-monotonic variation in the compressive stress state (Figure 2c) with a rapid increase up to the maximum value in a short time-lapse, followed by a gradual decrease up to a plateau value.
The first term in (51) is the force due to a thermal variation and the second term is that one due to a rotational velocity; moreover, r Γ defines the body force component in the r direction 2 Γ ω r r ρ = (53) whereas the governing equations of the problem take the following matrix form This equation must be solved at each time step with the known temperature variation as derived from Equation (49).

Verification
In this section we present the results from a transient stress analysis for a simplysupported FG-GPL porous rotating truncated cone, defined geometrically by

Transient and Steady Responses
The study continuous with the analysis of the transient thermal response of a FG-GPL porous rotating truncated cone, with material properties in Table 1 for both GPLs and metal matrix. As plotted in Figure 3, a systematic study checks for the influence of the GPL pattern on the thermal time history (Figure 3a Figure 3, the effect of the GPL pattern is much more pronounced on the transient thermal and kinematic response compared to the stress one. More specifically, it is worth observing that the temperature response maintains almost the same for a GPLUD and GPLV distribution. At the same time, GPLUD, GPLV, and GPLO patterns yield the same effect on the radial stress response. Besides, the effect of GPLUD and GPLV patterns on the other stress components seems almost the same. The results also denote that the maximum and minimum compressive axial and radial stress is reached for a GPLA and GPLO, respectively. In addition, the maximum and minimum steady state of radial, axial, and shear stresses belongs to a GPLO and GPLA distribution, respectively, while the maximum and minimum steady state of tangential stress and radial displacement belongs to a GPLO and GPLX distribution. Figure 4 indicates the influence of porosity distribution on the time history of temperature (Figure 4a), radial displacement (Figure 4b), and stresses (Figure 4c-e) for the central point of the truncated

Transient and Steady Responses
The study continuous with the analysis of the transient thermal response of a FG-GPL porous rotating truncated cone, with material properties in Table 1 for both GPLs and metal matrix. As plotted in Figure 3, a systematic study checks for the influence of the GPL pattern on the thermal time history (Figure 3a), radial displacement (Figure 3b), and stress state (Figure 3c-e) for the centered point of the selected truncated cone. Hereafter, the stress state [σ] T = [σ rr σ θθ σ zz τ rz ] T will be labeled as [σ] T = [σ r σ t σ z τ] T for simplicity reasons. Based on the plots in Figure 3, the effect of the GPL pattern is much more pronounced on the transient thermal and kinematic response compared to the stress one. More specifically, it is worth observing that the temperature response maintains almost the same for a GPLUD and GPLV distribution. At the same time, GPLUD, GPLV, and GPLO patterns yield the same effect on the radial stress response. Besides, the effect of GPLUD and GPLV patterns on the other stress components seems almost the same. The results also denote that the maximum and minimum compressive axial and radial stress is reached for a GPLA and GPLO, respectively. In addition, the maximum and minimum steady state of radial, axial, and shear stresses belongs to a GPLO and GPLA distribution, respectively, while the maximum and minimum steady state of tangential stress and radial displacement belongs to a GPLO and GPLX distribution. Figure 4 indicates the influence of porosity distribution on the time history of temperature (Figure 4a), radial displacement (Figure 4b), and stresses (Figure 4c-e) for the central point of the truncated cone. The results show that the temperature and displacement field will be steady sooner for D1, and the maximum and minimum magnitude of steady state temperature and displacement belongs to D1 and D3, respectively. It is interesting to note that the distribution of porosity has a meaningful effect on the radial and tangential stress time histories, which, in turn, show a relaxation for a D1 porosity, since they change from a compressive to a traction state. Appl (e) (f)   (e) (f) In Figure 5 we also plot the geometrical effect of a semi-vertex angle on the time history of temperature (Figure 5a), radial displacement (Figure 5b), and stresses (Figure 5c-e) at the central point of the truncated cone. It can be seen that, by increasing a semi-vertex angle, the cone behavior will be steady sooner. Figure 6 also shows the effect of the porosity coefficient on the same time histories for the same central point of the selected structure. It is worth observing that an increased porosity coefficient yields an increased steady state radial displacement and stress field. Structures with higher levels of porosity will respond more quickly to a thermal variation and reach the equilibrium conditions in a faster way.
(e) (f) In Figure 5 we also plot the geometrical effect of a semi-vertex angle on the time history of temperature (Figure 5a), radial displacement (Figure 5b), and stresses (Figure 5ce) at the central point of the truncated cone. It can be seen that, by increasing a semi-vertex angle, the cone behavior will be steady sooner. Figure 6 also shows the effect of the porosity coefficient on the same time histories for the same central point of the selected structure. It is worth observing that an increased porosity coefficient yields an increased steady state radial displacement and stress field. Structures with higher levels of porosity will respond more quickly to a thermal variation and reach the equilibrium conditions in a faster way.  (e) (f)  Figure 7 also analyzes the sensitivity of the response for different weight fractions of nanofillers. While the weight fraction of nanofillers increases, the steady state response significantly decreases due to an increased structural stiffness. For higher weight fractions of nanofillers, indeed, the structure will be steady faster. As also plotted in Figure 8, we repeat the same systematic analysis for a varying circular velocity, where it is clearly evident that an increased circular velocity causes an increased steady state stress field (Figure  8a-d) and radial displacement (Figure 8e).  Figure 7 also analyzes the sensitivity of the response for different weight fractions of nanofillers. While the weight fraction of nanofillers increases, the steady state response significantly decreases due to an increased structural stiffness. For higher weight fractions of nanofillers, indeed, the structure will be steady faster. As also plotted in Figure 8, we repeat the same systematic analysis for a varying circular velocity, where it is clearly evident that an increased circular velocity causes an increased steady state stress field (Figure 8a  The influence of various GPL patterns on temperature, radial displacement, and stresses along the r direction (at / 2 z L = ) is depicted in Figure 9. Note that the effect of a GPLV and GPLUD pattern on the overall thermo-mechanical response is quite similar. As also plotted in Figure 10, for different porosity patterns, the D1, D3, and D2 patterns provide a higher level of temperature and radial displacement, respectively. The effect of The influence of various GPL patterns on temperature, radial displacement, and stresses along the r direction (at z = L/2) is depicted in Figure 9. Note that the effect of a GPLV and GPLUD pattern on the overall thermo-mechanical response is quite similar. As also plotted in Figure 10, for different porosity patterns, the D1, D3, and D2 patterns provide a higher level of temperature and radial displacement, respectively. The effect of D1 on the radial stress state is nonlinear, while D3 and D2 provide almost a linear distribution in the radial direction. The influence of the porosity coefficient on temperature, radial displacement, and stresses along the r-direction at z = L/2 is shown in Figure 11. As expected, by increasing the porosity coefficient, the radial displacement increases, owing the structure to achieve a lower stiffness. At the same time, the temperature increases for an increased porosity level in the structure; see also the final contour plots in Figure 12 for the radial, axial, tangential, and shear stresses in Figure 12 in the steady state situation.

Conclusions
This work focused on the transient thermal stress response of porous rotating truncated cones reinforced by GPLs based on a classical linear thermos-elasticity. A finite element approach based on a Rayleigh-Ritz formulation and Crank-Nicolson time-stepping algorithm was employed to model and solve the problem, while checking for the influence of the porosity coefficient and distributions, GPL dispersion pattern, weight fraction of nanofillers, semi vertex angle, and rotational velocity, on the transient temperature, radial displacement, and stress state of the structure. Based on a large systematic investigation, it seems that GPLUD and GPLV distributions of GPLs provide almost the same temperature response. In addition, GPLUD, GPLV, and GPLO patterns have the same effect on the radial stress response. Besides, the effect of GPLUD and GPLV patterns on the other stress components seems almost the same. Results denote that the maximum and minimum compressive axial and radial stress state is reached for a GPLA and GPLO, respectively. Moreover, the maximum and minimum steady state for the radial, axial, and shear stresses correspond to GPLO and GPLA, respectively, while the maximum and minimum steady state of tangential stress and radial displacement is reached for a GPLO and GPLX distribution. The temperature and displacement response will be steady sooner for a D1 porosity, and the maximum and minimum magnitude of steady state temperature and displacement correspond to a D1 and D3 porosity, respectively.

Conclusions
This work focused on the transient thermal stress response of porous rotating truncated cones reinforced by GPLs based on a classical linear thermos-elasticity. A finite element approach based on a Rayleigh-Ritz formulation and Crank-Nicolson time-stepping algorithm was employed to model and solve the problem, while checking for the influence of the porosity coefficient and distributions, GPL dispersion pattern, weight fraction of nanofillers, semi vertex angle, and rotational velocity, on the transient temperature, radial displacement, and stress state of the structure. Based on a large systematic investigation, it seems that GPLUD and GPLV distributions of GPLs provide almost the same temperature response. In addition, GPLUD, GPLV, and GPLO patterns have the same effect on the radial stress response. Besides, the effect of GPLUD and GPLV patterns on the other stress components seems almost the same. Results denote that the maximum and minimum compressive axial and radial stress state is reached for a GPLA and GPLO, respectively. Moreover, the maximum and minimum steady state for the radial, axial, and shear stresses correspond to GPLO and GPLA, respectively, while the maximum and minimum steady state of tangential stress and radial displacement is reached for a GPLO and GPLX distribution. The temperature and displacement response will be steady sooner for a D1 porosity, and the maximum and minimum magnitude of steady state temperature and displacement correspond to a D1 and D3 porosity, respectively.
The distribution of porosity has a meaningful effect on the stress time histories. More specifically, a D1 porosity gets a relaxation in the stress state, since the nature of radial and tangential stresses changes from a compressive to a tensional state. Moreover, the influence of a D1 porosity on the radial stress state is nonlinear, while a D3 and D2 porosity features almost a linear distribution throughout the radial direction.
For an increased semi-vertex angle, the cone behavior will be steady sooner. The steady state radial displacement and stress field increases for an increased porosity coefficient. A conic structure with an increased amount of porosity will also react more quickly to a thermal condition variation, reaching a faster equilibrium condition. For an increased weight fraction of GPLs, the steady state radial displacement significantly decreases, and the structure becomes steady faster. At the same time, conic structures with an increased circular velocity feature an increased steady state radial displacement and stress field. The present findings could serve as theoretical predictions for design purposes of high-speed centrifugal separators, gas turbines, or high-power aircraft jet engines.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
The direct integration method may be used to integrate Equation (49)  where β is a constant which may be selected by the analyst. Equation (49) is written at times t and t + ∆t, the first equation is multiplied by (1 − β) and the second equation is multiplied by β to give These two equations are added and Equation (A1) is used to eliminate the time derivatives of {θ}. The resulting equation is solved for {θ} t+∆t , which yields For β = 1 2 , we obtain the Crank-Nicolson algorithm, which is unconditionally stable, regardless of the value of ∆t.