Thermally Driven Continuous Rolling of a Thick-Walled Cylindrical Rod

Self-sustained motion can take advantage of direct energy extraction from a steady external environment to maintain its own motion, and has potential applications in energy harvesting, robotic motion, and transportation. Recent experiments have found that a thermally responsive rod can perform self-sustained rolling on a flat hot plate with an angular velocity determined by the competition between the thermal driving moment and the friction moment. A rod with a hollow cross section tends to greatly reduce the frictional resistance, while promising improvements in thermal conversion efficiency. In this paper, through deriving the equilibrium equations for steady-state self-sustained rolling of the thick-walled cylindrical rod, estimating the temperature field on the rod cross-section, and solving the analytical solution of the thermally induced driving moment, the dynamic behavior of the thermally driven self-sustained rolling of the thick-walled cylindrical rod is theoretically investigated. In addition, we investigate in detail the effects of radius ratio, heat transfer coefficient, heat flux, contact angle, thermal expansion coefficient, and sliding friction coefficient on the angular velocity of the self-sustained rolling of the thick-walled cylindrical rod to obtain the optimal ratio of internal and external radius. The results are instructive for the application of thick-walled cylindrical rods in the fields of waste heat harvesters and soft robotics.

Recently, the study of self-sustained motion for thermally responsive polymeric material has gained plenty of attention [4,17,55].It has been observed experimentally that a polymer rod can roll or rotate steadily on a flat hot plate with uniform temperature.During its steady motion, the driving moment originates from the inhomogeneous expansion of the thermally responsive polymeric material.The self-sustained motion device can continuously obtain energy from the steady environment to maintain its motion, converting thermal energy into mechanical energy.Li et al. developed a thermodynamic model of a rod rolling on a hot plate and calculated the thermally induced deformation and stress fields of the rod during steady-state rolling, which also successfully predicted the bistability of the rod [55].In addition, Du et al. developed a theoretical model of thermally driven self-sustained rotation of a hollow torus, and found that for a given heat flux, there exists an optimal radius ratio maximizing the energy efficiency [17].Similar to other self-sustained motions [56][57][58], self-sustained rolling can harvest energy directly from the environment to maintain its own motion under the action of steady external stimuli [59][60][61], which has considerable application prospects in the field of soft robotics [13,15,62].
Similar to the overturning characteristics of the ring, the central part of the torus contributes less to the driving moment, but instead increases its friction moment and reduces the overturning angular velocity of the torus.Other things being equal, a rod with hollow section is expected to greatly reduce the friction moment and obtain a larger rolling angular velocity than a solid rod.It is possible for the rod to have an optimum radius ratio that maximizes the angular velocity and efficiency of the rod rolling, so it is of great interest to study thick cylindrical rod rolling on hot plate.In contrast to the hollow torus with constant curvature [17], the curvature of the steadily rolling thickwalled cylindrical rod varies with parameters such as heat flux, heat transfer coefficient, etc. Constrained by complex experimental conditions, this paper focuses on theoretical modelling to theoretically investigate the effects of various system parameters, including the radius ratio, on the angular velocity of a thermally responsive rod rolling itself on a flat plate; the effect of the radius ratio on the thermal conversion efficiency of a hollow self-sustaining rolling rod; and ultimately to predict the optimum internal and external radius ratio.
According to the existing literature, there are many studies on the self-rolling of the solid rod, but that of the thick-walled cylindrical rod has not been explored.The objective of this paper is to theoretically analyze the influence of a hollow section on a rod's selfrolling, with the aim of playing a guiding role in improving the energy efficiency.The layout of this paper is as follows.In Section 2, the theoretical model of thermally driven self-sustained rolling of a thick-walled cylindrical rod on a hot plate is established, the temperature field on the rod cross-section is studied, and the analytical solution of the thermally induced driving moment of the rod is derived.In Section 3, the equilibrium equation for the steady-state self-sustained rolling of a thick-walled cylindrical rod is given.Meanwhile, the effects of radius ratio, heat transfer coefficient, heat flux, contact angle, thermal expansion coefficient, and sliding friction coefficient on the rolling angular velocity of the thick-walled cylindrical rod are investigated in detail.In addition, the effects of radius ratio and dimensionless heat flux on the energy efficiency are also investigated.In Section 4, a summary of this study is presented.

Temperature Field in the Steadily Rolling Rod
Our model is sketched in Figure 1.As shown in Figure 1a,b, a thick-walled cylindrical rod with internal radius a, external radius b, and length L was placed on a horizontal hot plate.For thermal expansion materials, the rod rolled in the concave direction (Figure 1a), and for thermal shrinkage materials, the rod rolled in the convex direction as shown in Figure 1b and Video S1.In Video S1, we use the nylon 6 material as an example.Figure 1c shows the driving moment applied to the thick-walled cylindrical rod and the friction moment acting on the rod. Figure 1d shows the stresses developed in the section of the thick-walled cylindrical rod and the moments in all directions.In Figure 1e, q indicates the heat flux, T e is the external environment temperature, and θ 0 is the half contact angle.
Micromachines 2022, 13, x FOR PEER REVIEW Figure 1b and Video S1.In Video S1, we use the nylon 6 material as an example.Fig shows the driving moment applied to the thick-walled cylindrical rod and the fr moment acting on the rod. Figure 1d shows the stresses developed in the section thick-walled cylindrical rod and the moments in all directions.In Figure 1e, q ind the heat flux, e T is the external environment temperature, and 0  is the half c angle.In this study, our proposed model adopts the direct heating method.Of cours entirely possible to use indirect heating methods such as light or electromagnetic he as in the study of Ref. [14], the self-sustained rolling of solid LCE rod was achieved light.Therefore, the light-powered self-sustained rolling of the hollow LCE rod wi very worthwhile study.It is worth mentioning that, in order to simplify the proble assume that a thermal insulation layer is coated on the inner wall of the thick-wall lindrical rod, and the inner wall of the thick-walled cylindrical rod was seen as an batic boundary.
The coordinate systems in this paper are established with respect to plane, as s in Figure 1.The Poisson effect of the thick-walled cylindrical rod and the rod end are neglected in the following analysis of this paper, so as to simplify the problem plane problem. Figure 1c shows the driving moment applied to the thick-walled cy cal rod and the friction moment acting on the rod.In addition, since the deformat this study is small, we assume that the contact angle is constant in the calculation rolling angular velocity of the thick-walled cylindrical rod at moment t is repres by ( ) t  .The heat transfer coefficient, the heat conduction coefficient, the mass de The stresses on the cross-section of the thick-walled cylindrical rod and the moments in all directions.(e) The cross-section of the thick-walled cylindrical rod.Considering that the length of the rod is much greater than its thickness, we assume that the temperature distribution in the rod cross-section remains the same throughout the length of the rod.An inhomogeneous temperature field in the rod cross-section leads to uniform thermal expansion/shrinkage, which generates a driving moment M drive that equilibrates with the friction moment M f and causes the thick-walled cylindrical rod to self-roll.
In this study, our proposed model adopts the direct heating method.Of course, it is entirely possible to use indirect heating methods such as light or electromagnetic heating, as in the study of Ref. [14], the self-sustained rolling of solid LCE rod was achieved using light.Therefore, the light-powered self-sustained rolling of the hollow LCE rod will be a very worthwhile study.It is worth mentioning that, in order to simplify the problem, we assume that a thermal insulation layer is coated on the inner wall of the thick-walled cylindrical rod, and the inner wall of the thick-walled cylindrical rod was seen as an adiabatic boundary.
The coordinate systems in this paper are established with respect to plane, as shown in Figure 1.The Poisson effect of the thick-walled cylindrical rod and the rod end effect are neglected in the following analysis of this paper, so as to simplify the problem as a plane problem. Figure 1c shows the driving moment applied to the thick-walled cylindrical rod and the friction moment acting on the rod.In addition, since the deformation in this study is small, we assume that the contact angle is constant in the calculation.The rolling angular velocity of the thick-walled cylindrical rod at moment t is represented by ω(t).The heat transfer coefficient, the heat conduction coefficient, the mass density, and the specific heat of the rod are denoted by h, k, ρ and c, respectively.During the steady-state rolling of the thick-walled cylindrical rod, we have T(r, θ, t) = T(r, θ) and ω(t) = ω.
For the problem studied in this paper, we assume that except for the small area in contact with the hot plate, the remaining external surface of the thick-walled cylindrical rod is cooled by convection, while the internal surface of the rod is adiabatic.Due to the complexity of the thermal boundary of the rod, the thermal radiation effect has been neglected in order to simplify the analysis.To facilitate the calculation and analysis, we introduce the dimensionless parameters: radius ratio λ = a/b, r = r/b, T = T/T e , q = qb/kT e , h = hb/k, β m = β m b, and ω = ωb 2 /ψ.Here, T is the dimensionless temperature filed, q is the dimensionless heat flux, h is the dimensionless heat transfer coefficient, β m is the dimensionless air damping, and ω is the dimensionless rolling angular velocity.In the steady-state rolling, the temperature field distribution in the rod crosssection can be expressed by following the previous work as [17,63,64] where where B e λ, h and B i λ, h are given by with β m λ, h being the positive root of the following characteristic equation [64] with J n being the first-class Bessel functions, and Y n the second-class Bessel functions.
From the accessible experiments [14,42], the typical values of material properties and geometrical parameters available for the experiments are presented in Table 1, and the estimates of dimensionless parameters are listed in Table 2.The deformation of the thickwalled cylindrical rod for the given parameters in this paper is small.Figure 2a-d shows the temperature fields on the cross section of a thick-walled cylindrical rod in steady-state rolling for different combinations of rolling angular velocity ω and heat flux q.We set h = 0.3, λ = 0.3, and θ 0 = 0.2.For a given heat flux q, the temperature difference in the rod cross-section increases and then decreases with the increase in rolling angular velocity ω, and the temperature field tends to be uniform.In contrast, for a given rolling angular velocity ω, the temperature difference in the rod cross-section increases with the increase in heat flux q, and the non-uniformity of the temperature field increases.
Figure 2e-h depicts the temperature fields on the cross section of a thick-walled cylindrical rod in steady-state rolling for different combinations of λ and h.In the computation, we set ω = 4, q = 35, and θ 0 = 0.2.For a given heat transfer coefficient h, along with the increase in radius ratio λ of the thick-walled cylindrical rod, we witness a slight increase in the temperature difference on the rod cross-section, and the temperature field tends to be inhomogeneous.This is because the smaller the radius ratio λ, the more unfavorable the heat transfer inside the thick-walled cylindrical rod, leading to a decrease in temperature in the low-temperature area and an increase in temperature in the high-temperature area.Additionally, for a given radius ratio λ, the temperature difference of the thick-walled cylindrical rod cross-section decreases with the increase in heat transfer coefficient h, and the temperature field tends to be homogeneous.This is mainly owing to the fact that a larger heat transfer coefficient is conducive to heat transfer, resulting in a decrease in temperature in the high temperature region of the cross section and an increase in temperature in the low temperature region, so that the temperature field tends to be homogeneous.2e-h.In this case, the rolling velocity of the rod is one of the input parameters.When the thick-walled cylindrical rod rolls steadily, the non-uniformity of the temperature field on the cross section decreases with the increase in rolling angular velocity  and heat transfer coefficient h , and increases with the increase in heat flux q and radius ratio  .Temperature fields on the cross section of a thick-walled cylindrical rod in steady-state rolling for different combinations of rolling angular velocity ω and heat flux q: (a) ω = 4, q = 20, (b) ω = 8, q = 20, (c) ω = 4, q = 35, (d) ω = 8, q = 35, and different combinations of radius ratio λ and heat transfer coefficient h: (e) λ = 0.1, h = 0.3, (f) λ = 0.3, h = 0.3, (g) λ = 0.1, h = 0.6, (h) λ = 0.3, h = 0.6.We set h = 0.3, λ = 0.3, and θ 0 = 0.2 in a-d, and ω = 4, q = 35 and θ 0 = 0.2 in e-h.In this case, the rolling velocity of the rod is one of the input parameters.When the thick-walled cylindrical rod rolls steadily, the non-uniformity of the temperature field on the cross section decreases with the increase in rolling angular velocity ω and heat transfer coefficient h, and increases with the increase in heat flux q and radius ratio λ.

Driving Moment for the Rolling of the Thick-Walled Cylindrical Rod
In this study, we assume that the thermally induced strain ε T is proportional to the temperature change on the rod cross-section for the purpose of analysis.The thermally induced strain is given by ε where C T is the thermal expansion coefficient of the material.For thermal expansion material, C T is positive, while for thermal shrinkage material, C T is negative.
Through applying the linear thermoelastic model, we can calculate the thermal stress σ(r, θ) on the rod cross-section along its normal direction as where E is the elastic modulus of the material.The thermal bending moment of the rod around z-axis on the cross section can be described as Considering that the rod can bend freely in the lateral plane, the thermally driven lateral curvature can be calculated as During the steady-state rolling of the thick-walled cylindrical rod, the total strain on the rod cross-section can be expressed as Thus, the axial stress on the rod cross-section can be obtained as, and the direction of the axial stress on the cross-section of the rod is shown in Figure 1d.The total bending moment of the rod around x-axis on the cross section can be described as Therefore, the total net moment applied on the thick-walled cylindrical rod can be expressed as [55] M For the convenience of computational analysis, we introduce 2 , and κ z = κ z b.By combining Equations ( 1) and ( 11)-( 18), the driving moment of a thick-walled cylindrical rod in a steady rolling process can be rewritten as The details of the derivation of Equation ( 19) are given in the Appendix A. In Equation (19), P λ, h can be expressed as with MeijerG {s 1 . . .s n }, s n+1 . . .s p , {t 1 . . .t m }, s m+1 . . .s q , z being defined as the MeijerG function.
Figure 3 reflects the influence of ω and λ of the thick-walled cylindrical rod on the lateral curvature κ z .In the calculation, we set h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.It is apparent from the diagram that for a given radius ratio λ, the lateral curvature κ z increases first and then decreases with the increasing rolling angular velocity ω.This is because the inhomogeneity of the steady-state temperature field in the rod cross-section increases and then decreases with the increase in rolling velocity (as shown in Figure 2).Meanwhile, the thermal bending moment and the thermally driven lateral curvature are positively correlated with the inhomogeneity of the temperature field (which can be derived from Equations ( 13) and ( 14)).For a given rolling angular velocity ω, in the process of enlarging the radius ratio λ, the lateral curvature κ z shows a trend of rising first and then falling.This is because the inhomogeneity of the temperature field in the rod cross-section increases slightly with the increase in λ (as shown in Figure 2).It is easily concluded from the figure that under a given rolling angular velocity ω, there is an optimal radius ratio for the thick-walled cylindrical rod that maximizes the lateral curvature κ z .In addition, it is worth mentioning that the dimensionless heat transfer coefficient h can influence the lateral curvature κ z of the thick-walled cylindrical rod by affecting the temperature field of the rod cross-section.increases and then decreases with the increase in rolling angular velocity  .This is because the inhomogeneity In the computation, we set h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.In this case, the rolling velocity of the rod is one of the output parameters.For a given angular velocity ω, there exists a radius ratio λ corresponding to an optimal lateral curvature κ z .
Figure 4 illustrates the effects of the rolling angular velocity ω and the radius ratio λ of the thick-walled cylindrical rod on the driving moment M drive .In the computation, we set: h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.As can be seen from the figure, for a given radius ratio λ, the driving moment M drive increases and then decreases with the increase in rolling angular velocity ω.This is because the inhomogeneity of the temperature field of the rod cross-section increases and then decreases as the angular velocity increases, as shown in Figure 2, which ultimately leads to an increase and then a decrease in the driving moment M drive .For a given rolling angular velocity ω, the driving moment M drive first increases and then decreases with the increase in radius ratio λ.This is because the temperature field in the rod cross-section becomes more nonuniform as the radius ratio λ grows larger.At a smaller radius ratio, the increased driving moment M drive due to the increase in temperature field inhomogeneity is large enough to compensate for the loss of driving moment M drive due to the reduction in the rod cross-section.For a given rolling angular velocity ω, there exists an optimal radius ratio of the thick-walled cylindrical rod corresponding to a maximum driving moment.However, we have to note that the actual maximum rolling angular velocity is dependent on the competition between the driving moment and the friction moment.

Self-Rolling of the Thick-Walled Cylindrical Rod on a Hot Plate
Considering the driving moment in Equation ( 19), the equilibrium equation during steady rolling of the thick-walled cylindrical rod is further derived in this section.Then, the effects of radius ratio, thermal expansion coefficient, heat flux, contact angle, heat transfer coefficient and sliding friction coefficient on the rolling angular velocity of the thick-walled cylindrical rod are studied in detail, and the critical value for triggering the self-sustained rolling is found.Furthermore, the dependence of the energy efficiency on the radius ratio is also investigated.

Equilibrium Equations
During the steady rolling, the surface of the thick-walled cylindrical rod is also subjected to sliding friction between the rod and the hot plate.The magnitude of the friction force on the thick-walled cylindrical rod can be expressed as  .The effects of rolling angular velocity ω and radius ratio of the thick-walled cylindrical rod λ on the driving moment M drive .In the computation, we set: h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.In this case, the rolling velocity of the rod is one of the output parameters.For a given angular velocity ω, there exists an optimal radius ratio λ corresponding to a maximum driving moment M drive .

Self-Rolling of the Thick-Walled Cylindrical Rod on a Hot Plate
Considering the driving moment in Equation ( 19), the equilibrium equation during steady rolling of the thick-walled cylindrical rod is further derived in this section.Then, the effects of radius ratio, thermal expansion coefficient, heat flux, contact angle, heat transfer coefficient and sliding friction coefficient on the rolling angular velocity of the thick-walled cylindrical rod are studied in detail, and the critical value for triggering the self-sustained rolling is found.Furthermore, the dependence of the energy efficiency on the radius ratio is also investigated.

Equilibrium Equations
During the steady rolling, the surface of the thick-walled cylindrical rod is also subjected to sliding friction between the rod and the hot plate.The magnitude of the friction force on the thick-walled cylindrical rod can be expressed as where C f denotes the sliding friction coefficient between the rod and the hot plate, and g is the gravitational acceleration.Moreover, the friction moment on the thick-walled cylindrical rod can be expressed as During the steady-state rolling, the driving moment is equal to the friction moment.By combining Equations ( 19) and ( 22) and separating the fixed parameters, we can derive the equilibrium equation for the steady-state rolling as where It can be observed from Equation ( 23) that although there are many parameters related to the angular velocity ω, including C f , ρ, b, E, q, C T , θ 0 , and h, we only need to analyze the effects of three parameters, namely radius ratio λ, heat-transfer coefficient h and dimensionless parameter Ω = C f ρgr 2 /E qC T sin θ 0 2 , on the rolling angular velocity ω.

Angular Velocity of the Self-Rolling of Thick-Walled Cylindrical Rod
Figure 5 shows the effect of dimensionless parameter Ω on the rolling angular velocity ω of the thick-walled cylindrical rod for different radius ratios λ.In the computation, we set h = 0.3.As can be seen from the diagram, there is a critical Ω crit for maintaining the steady-state rolling of the thick-walled cylindrical rod, and Ω crit increases with the increase in radius ratio λ.Obviously, the static state of the thick-walled cylindrical rod is always a steady state.When Ω increases to a critical the rod has another high-velocity steady state.Between static state and high velocity state, the rod has an intermediate velocity state, which is an unstable state because a perturbation of increasing (decreasing) the speed of the rod increases (decreases) its driving moment, leading to further increase or decrease in the rolling speed of the rod [55].From Equations ( 21)- (23) and Ω = C f ρgr 2 /E qC T sin θ 0 2 , it can been seen that the dimensionless parameter Ω increases as the friction moment M f increases or as the driving moment M drive decreases, so that the rolling angular velocity decreases.
Taking Ω = C f ρgr 2 /E qC T sin θ 0 2 into consideration, we can also conclude that the rolling angular velocity ω will increase with the increases in heat flux q, contact angle θ 0 , thermal expansion coefficient C T , external radius b, and elastic modulus of the material E, and that conversely, it will decrease with increases in sliding friction coefficient C f and mass density ρ.The effects of these parameters on the rolling angular velocity of the thick-walled cylindrical rod are consistent with that of the solid rod on a hot plate [7,55].Figure 6 presents the effect of heat transfer coefficient h on the rolling angular velocity ω of the thick-walled cylindrical rod for different radius ratios λ.In the computation, we set Ω = 4.As observed from the plots, there is a critical heat transfer coefficient h crit for keeping the steady rolling of the thick-walled cylindrical rod, and h crit increases with the increase in radius ratio λ.It is obvious that the static state of the thick-walled cylindrical rod is always a stable state.When the heat transfer coefficient h approaches the critical value, the rod presents another high-velocity steady state.Between the static state and high velocity, the rod has an intermediate rolling velocity, behaving as an unstable state.When the thick-walled cylindrical rod is rolling at a high velocity, under the given radius ratio λ, the increase in heat transfer coefficient h causes the rolling angular velocity ω to increase first and then decrease.It is well understood from a physical point of view that for a smaller heat transfer coefficient h, the inhomogeneity of the steady-state temperature field in the rod cross-section increases with the increase in heat transfer coefficient h.an intermediate velocity state, which is an unstable state because a perturbation of increasing (decreasing) the speed of the rod increases (decreases) its driving moment, leading to further increase or decrease in the rolling speed of the rod [55].From Equations ( 21)-( 23) and ( )    ) into consideration, we can also conclude that the rolling angular velocity  will increase with the increases in heat flux q , contact angle There is a critical Ω crit for maintaining the steady-state rolling of the thick-walled cylindrical rod, and Ω crit increases with the increase in radius ratio λ. Figure 6 presents the effect of heat transfer coefficient h on the rolling angular ve- locity  of the thick-walled cylindrical rod for different radius ratios  .In the compu- tation, we set 4 Ω = .As observed the plots, there is a critical heat transfer coefficient crit h for keeping the steady rolling of the thick-walled cylindrical rod, and crit h increases with the increase in radius ratio  .It is obvious that the static state of the thick-walled cylindrical rod is always a stable state.When the heat transfer coefficient h approaches the critical value, the rod presents another high-velocity steady state.Between the static state and high velocity, the rod has an intermediate rolling velocity, behaving as an unstable state.When the thick-walled cylindrical rod is rolling at a high velocity, under the given radius ratio  , the increase in heat transfer coefficient h causes the rolling an- gular velocity  to increase first and then decrease.It is well understood from a physical point of view that for a smaller heat transfer coefficient h , the inhomogeneity of the steady-state temperature field in the rod cross-section increases with the increase in heat transfer coefficient h .When the heat transfer coefficient h is large, the inhomogeneity of the steady-state temperature field decreases with the increase in heat transfer coefficient h .For a given heat transfer coefficient h , the rolling angular velocity  increases with the increase in radius ratio  .This is because when the radius ratio  increases, the heat transfer inside the thick-walled cylindrical rod is more favorable, and in turn the rolling angular When the heat transfer coefficient h is large, the inhomogeneity of the steady-state temperature field decreases with the increase in heat transfer coefficient h.For a given heat transfer coefficient h, the rolling angular velocity ω increases with the increase in radius ratio λ.This is because when the radius ratio λ increases, the heat transfer inside the thick-walled cylindrical rod is more favorable, and in turn the rolling angular velocity ω is greater.For a larger heat transfer coefficient h, the temperature difference on the rod cross-section decreases with the increase in heat transfer coefficient h, and the temperature field tends to be homogeneous.
Figure 7 shows the effect of radius ratios λ on the rolling angular velocity ω of the thick-walled rod for different heat fluxes q.In the computation, we set h = 0.3.As can be seen from the diagram, a critical radius ratio λ crit exists to maintain the steady rolling of the thick-walled cylindrical rod, and λ crit decreases as the heat flux q increases.It can be clearly seen that static state of the rod is always a stable state.When the radius ratio λ is increased towards a critical value, there is another stable state of the rod with high velocity.Between static and high velocity, there is an intermediate rolling velocity of the rod, which is an unstable state.For the thick-walled cylindrical rod rolling at high velocity, under a given heat flux q, the rolling angular velocity ω increases with the increase in radius ratios λ.This is due to the fact that the larger the radius ratio λ, the more conducive it is to the heat transfer inside the rod, the larger the rolling angular velocity ω.
ence on the rod cross-section decreases with the increase in heat transfer coefficient h , and the temperature field tends to be homogeneous.
Figure 7 shows the effect of radius ratios  on the rolling angular velocity  of the thick-walled rod for different heat fluxes q .In the computation, we set 3 .0 = h .As can be seen from the diagram, a critical radius ratio crit  exists to maintain the steady rolling of the thick-walled cylindrical rod, and crit  decreases as the heat flux q in- creases.It can be clearly seen that static state of the rod is always a stable state.When the radius ratio  is increased towards a critical value, there is another stable state of the rod with high velocity.Between static and high velocity, there is an intermediate rolling velocity of the rod, which is an unstable state.For the thick-walled cylindrical rod rolling at high velocity, under a given heat flux q , the rolling angular velocity  increases with the increase in radius ratios  .This is due to the fact that the larger radius ratio  , the more conducive it is to the heat transfer inside the rod, the larger the rolling angu- lar velocity  .For a given radius ratio  , the rolling angular velocity  increases with the in- crease in heat flux q .The reason for this phenomenon is that with the increase in heat flux q , the temperature difference on the rod cross-section increases, and the non-uni- formity of temperature field on the cross section increases, eventually resulting in the increases in driving moment drive M and rolling angular velocity  .

Energy Efficiency of the Self-Rolling Thick-Walled Cylindrical Rod
The self-rolling system studied in this paper has the potential to be applied as a thermally driven motor or energy harvester.In the above theoretical model, the thick-walled cylindrical rod compensates the energy dissipated by damping through absorbing the thermal energy from the environment.If it is used as a thermally driven motor or energy harvester, we can regard the work performed by the damping force as the effective work For a given radius ratio λ, the rolling angular velocity ω increases with the increase in heat flux q.The reason for this phenomenon is that with the increase in heat flux q, the temperature difference on the rod cross-section increases, and the non-uniformity of temperature field on the cross section increases, eventually resulting in the increases in driving moment M drive and rolling angular velocity ω.

Energy Efficiency of the Self-Rolling Thick-Walled Cylindrical Rod
The self-rolling system studied in this paper has the potential to be applied as a thermally driven motor or energy harvester.In the above theoretical model, the thickwalled cylindrical rod compensates the energy dissipated by damping through absorbing the thermal energy from the environment.If it is used as a thermally driven motor or energy harvester, we can regard the work performed by the damping force as the effective work output of the system.In the steady rolling of the thick-walled cylindrical rod, the input thermal power can be expressed in terms of heat flux density as P in = 2qθ 0 Lb and the effective power of the rod as P e = M f ω.Combining Equations ( 21) and ( 22), the energy efficiency of a self-rolling thick-walled cylindrical rod is obtained from the following formula: It can be clearly seen from Equation ( 25) that the energy efficiency η is related to many parameters, including C f , ρ, b, a, q, and θ 0 .In the following, we take the combination of radius ratio λ and dimensionless heat flux q as an example to study the energy efficiency of the self-rolling thick-walled cylindrical rod.The typical values of parame- Figure 8 illustrates the dependence of energy efficiency η on the radius ratio λ and dimensionless heat flux q.The parameters are set: h = 0.3, κ z = 0.025, C T = 5 × 10 −3 , C f = 0.45 and θ 0 = 0.2.The results show that for a stationary thick-walled cylindrical rod, the energy efficiency of the rod is zero.When the thick-walled cylindrical rod undergoes a steady-state rolling, for a given heat flux q, the energy efficiency increases first and then gradually decreases along with the increase in radius ratio.There exists an optimal radius ratio maximizing the energy efficiency of the rod, and the optimal radius ratio presents a certain decrease as the heat flux increases.This results from the competition between the kinetic energy and the heat dissipation, which both increase with the increase in the radius ratio.Therefore, the radius ratio of the thick-walled cylindrical rod is not as large as possible.In practical applications, it is necessary to seek a balance between the rolling angular velocity and the energy efficiency.
It can be clearly seen from Equation ( 25) that the energy efficiency  is related to many parameters, including f C ,  , b , a , q , and 0  .In the following, we take the combination of radius ratio  and dimensionless heat flux q as an example to study the energy efficiency of the self-rolling thick-walled cylindrical rod.The typical values of


. The results show that for a stationary thickwalled cylindrical rod, the energy efficiency of the rod is zero.When the thick-walled cylindrical rod undergoes a steady-state rolling, for a given heat flux q , the energy effi- ciency increases first and then gradually decreases along with the increase in radius ratio.There exists an optimal radius ratio maximizing the energy efficiency of the rod, and the optimal radius ratio presents a certain decrease as the heat flux increases.This results from the competition between the kinetic energy and the heat dissipation, which both increase with the increase in the radius ratio.Therefore, the radius ratio of the thick-walled cylindrical rod is not as large as possible.In practical applications, it is necessary to seek a balance between the rolling angular velocity and the energy efficiency.There exists an optimal radius ratio that maximizes the energy efficiency of the thick-walled cylindrical rod, and the optimal radius ratio decreases as the heat flux increases.The parameters are h = 0.3, κ z = 0.025, C T = 5 × 10 −3 , C f = 0.45, and θ 0 = 0.2.There exists an optimal radius ratio that maximizes the energy efficiency of the thick-walled cylindrical rod, and the optimal radius ratio decreases as the heat flux increases.
Table 3 lists the optimal radius ratio, the maximum energy efficiency of the thickwalled cylindrical rod, the energy efficiency of the solid rod and the energy efficiency improvement of the thick-walled cylindrical rod for three different heat fluxes q in Figure 8.As shown in Table 3, q = 30 corresponds to the maximum energy efficiency improvement, i.e., the increase from 0% to 16.95%.Overall, the results of this study show that, in practical applications, we can improve energy efficiency by adjusting the radius ratio of thick-walled cylindrical rod.Similarly, by substituting Equations ( 1), ( 11) and ( 16) into Equation ( 17) and simplifying, the analytical solution for the total bending moment of the rod around x-axis can be obtained as Equation ( 19) can be obtained by substituting Equations (A3) and (A4) Equation ( 18) then integrating.

Figure 1 .
Figure 1.Schematic of a thermally driven rolling thick-walled cylindrical rod on a hot plate.thermal expansion materials, the thermal expansion coefficient

Figure 1 .
Figure 1.Schematic of a thermally driven rolling thick-walled cylindrical rod on a hot plate.(a) For thermal expansion materials, the thermal expansion coefficient C T is positive, and the rod rolls in a concave direction.(b) For thermal shrinkage materials, the thermal expansion coefficient C T is negative, and the rod rolls in a convex direction.(c) The driving moment applied to the thick-walled cylindrical rod and the friction moment acting on the rod.(d)The stresses on the cross-section of the thick-walled cylindrical rod and the moments in all directions.(e) The cross-section of the thick-walled cylindrical rod.Considering that the length of the rod is much greater than its thickness, we assume that the temperature distribution in the rod cross-section remains the same throughout the length of the rod.An inhomogeneous temperature field in the rod cross-section leads to uniform thermal expansion/shrinkage, which generates a driving moment M drive that equilibrates with the friction moment M f and causes the thick-walled cylindrical rod to self-roll.

Figure 2 .
Figure 2. Temperature fields on the cross section of a thick-walled cylindrical rod in steady-state rolling for different combinations of rolling angular velocity  and heat flux q : (a)

Figure 3 ..
Figure 3.The effects of rolling angular velocity  and radius ratio  of the thick-walled cylin- drical rod on the lateral curvature z  .In the computation, we set 3 .0 = h , 10 = L , 3 T 10 5

Figure 4 .
Figure 4 illustrates the effects of the rolling angular velocity  and the radius ratio  of the thick-walled cylindrical rod on the driving moment drive M. In the computation,

Figure 3 .
Figure 3.The effects of rolling angular velocity ω and radius ratio λ of the thick-walled cylindrical rod on the lateral curvature κ z .In the computation, we set h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.In this case, the rolling velocity of the rod is one of the output parameters.For a given angular velocity ω, there exists a radius ratio λ corresponding to an optimal lateral curvature κ z .

Figure 4 ..
Figure 4.The effects of rolling angular velocity  and radius ratio of the thick-walled cylindrical rod  on the driving moment

C
denotes the sliding friction coefficient between the rod and the hot plate, and g is the gravitational acceleration.Moreover, the friction moment on the thick-walled cylindrical rod can be expressed as

Figure 4
Figure 4.The effects of rolling angular velocity ω and radius ratio of the thick-walled cylindrical rod λ on the driving moment M drive .In the computation, we set: h = 0.3, L = 10, C T = 5 × 10 −3 , q = 30 and θ 0 = 0.2.In this case, the rolling velocity of the rod is one of the output parameters.For a given angular velocity ω, there exists an optimal radius ratio λ corresponding to a maximum driving moment M drive .
can been seen that the dimensionless parameter Ω in- creases as the friction moment f M increases or as the driving moment drive M decreases, so that the rolling angular velocity decreases.

Figure 5 .
Figure 5.The effect of dimensionless parameter Ω on the rolling angular velocity  of the thick- walled cylindrical rod for

Figure 5 .
Figure 5.The effect of dimensionless parameter Ω on the rolling angular velocity ω of the thickwalled cylindrical rod for h = 0.3.The rolling angular velocity ω decreases with the increase in Ω.There is a critical Ω crit for maintaining the steady-state rolling of the thick-walled cylindrical rod, and Ω crit increases with the increase in radius ratio λ.

Micromachines 2022 , 18 0
13,  x FOR PEER REVIEW 12 of , thermal expansion coefficient T C , external radius b , and elastic modulus of the material E , and that conversely, it will decrease with increases in sliding friction coeffi- cient f C and mass density  .The effects of these parameters on the rolling angular velocity of the thick-walled cylindrical rod are consistent with that of the solid rod on a hot plate[7,55].

Figure 6 .
Figure 6.The effect of heat transfer coefficient h on the rolling angular velocity  of the thick- walled cylindrical rod for 4 Ω= .The rolling angular velocity  increases first and then decreases with the increase in heat transfer coefficient h .There is a critical transfer coefficient crit h for keep-

Figure 6 .
Figure 6.The effect of heat transfer coefficient h on the rolling angular velocity ω of the thick-walled cylindrical rod for Ω = 4.The rolling angular velocity ω increases first and then decreases with the increase in heat transfer coefficient h.There is a critical transfer coefficient h crit for keeping the steady-state rolling of the thick-walled cylindrical rod, and h crit increases with the increase in radius ratio λ.

Figure 7 ..
Figure 7.The effect of radius ratio  on the rolling angular velocity  of the thick-walled cylin- drical rod for 3 .0 = h .The rolling angular velocity  increases with the increase in radius ratio  .There is a critical radius ratio crit  for keeping the steady rolling of the thick-walled cylindrical

Figure 7 .
Figure 7.The effect of radius ratio λ on the rolling angular velocity ω of the thick-walled cylindrical rod for h = 0.3.The rolling angular velocity ω increases with the increase in radius ratio λ.There is a critical radius ratio λ crit for keeping the steady rolling of the thick-walled cylindrical rod, and λ crit decreases with the increase in heat flux q.

Figure 8
Figure 8  illustrates the dependence of energy efficiency  on the radius ratio  and dimensionless heat flux q .The parameters are set:

Figure 8 .
Figure 8.The dependence of energy efficiency  on the radius ratio  and dimensionless heat flux q .The parameters are

Figure 8 .
Figure 8.The dependence of energy efficiency η on the radius ratio λ and dimensionless heat flux q.The parameters are h = 0.3, κ z = 0.025, C T = 5 × 10 −3 , C f = 0.45, and θ 0 = 0.2.There exists an optimal radius ratio that maximizes the energy efficiency of the thick-walled cylindrical rod, and the optimal radius ratio decreases as the heat flux increases.

Table 1 .
Material properties and geometric parameters.

Table 3 .
Energy efficiency improvement of the thick-walled cylindrical rod.