Dynamical Behaviors of a Translating Liquid Crystal Elastomer Fiber in a Linear Temperature Field

Liquid crystal elastomer (LCE) fiber with a fixed end in an inhomogeneous temperature field is capable of self-oscillating because of coupling between heat transfer and deformation, and the dynamics of a translating LCE fiber in an inhomogeneous temperature field are worth investigating to widen its applications. In this paper, we propose a theoretic constitutive model and the asymptotic relationship of a LCE fiber translating in a linear temperature field and investigate the dynamical behaviors of a corresponding fiber-mass system. In the three cases of the frame at rest, uniform, and accelerating translation, the fiber-mass system can still self-oscillate, which is determined by the combination of the heat-transfer characteristic time, the temperature gradient, and the thermal expansion coefficient. The self-oscillation is maintained by the energy input from the ambient linear temperature field to compensate for damping dissipation. Meanwhile, the amplitude and frequency of the self-oscillation are not affected by the translating frame for the three cases. Compared with the cases of the frame at rest, the translating frame can change the equilibrium position of the self-oscillation. The results are expected to provide some useful recommendations for the design and motion control in the fields of micro-robots, energy harvesters, and clinical surgical scenarios.


Introduction
The phenomenon of continuous periodic motion of a system under the influence of a steady external environment stimulation is known as self-excited motion [1][2][3][4]. To compensate for the energy dissipation consumed by the damping, the periodic motion of self-excited oscillation is maintained by collecting the energy directly from the external environment [5][6][7]. This property minimizes the complexity of self-excited oscillation systems, making complicated control system design easier to implement and allowing for high load capacity [8,9]. Self-excited motion possesses the autonomous characteristics [10][11][12][13][14], which aids in the study of non-equilibrium thermodynamic processes [15,16], and has several applications in the fields of active machinery [17][18][19][20], mobile robots [21,22], energy acquisition [23][24][25][26], and motors [27].
One of the prevalent external stimuli for triggering the actuation of responsive materials is heat [61][62][63]. In a steady-temperature field, self-excited oscillation of a liquid crystal elastomer fiber with a hanging weight was recently reported [46]. In the experiment, the anisotropic rod-like liquid crystalline molecules and stretchy long-chain polymers are used to synthesize and prepare the LCE material. The LCE microfiber is hung from a rigid flat plate in the temperature field and the mass is placed on the bottom free end of the fiber. The fiber system was discovered to be able to vibrate continuously in an inhomogeneous temperature field. Theoretical research demonstrates that the self-oscillation is caused by the combination between the vibration of the fiber and heat exchange [47]. When the fiber vibrates, it can collect heat energy from the environment, to compensate for the damping dissipation and keep its own motion [47]. This complex, thermodynamic-coupled nonlinear dynamics problem is further studied in [61]. It abstractly presented a fiber engine module and established its constitutive model and asymptotic relationship, which is similar to the Kelvin-Voigt viscoelastic model consisting of a dashpot and a spring. It is found that the one-end-fixed fiber in a linear temperature field has three kinds of behaviors: damper, spring, and engine.
Although considerable research has been carried out on the self-oscillation motion of a thermal-responsive LCE fiber with one end fixed [46,47,61], the dynamics of a translating thermally responsive fiber with a moving frame need to be explored for further possible applications. To investigate the dynamical behaviors of the LCE fiber translating in a temperature field varying linearly with position, we propose a theoretic constitutive model of LCE fiber and its asymptotic relationship and study the fiber-mass system in three typical cases of the frame, at rest, uniform, and accelerating translation. The objective of this work is to establish a constitutive model, obtain its dynamical behavior, and provide some useful recommendations for the design and motion control of micro-robots of various structures and applications.
The article is organized as follows. In Section 2, a constitutive model of translating LCE fiber in a linear temperature field is proposed, and its asymptotic relationship is derived for the small characteristic time. In Section 3, we establish a thermally responsive fiber-mass system with zero characteristic time, and formulate its corresponding governing equations based on the constitutive model, to investigate its dynamical behaviors for three typical cases of the frame, at rest, uniform, and accelerating translation. In Section 4, by considering the small characteristic time, we study the dynamics of the fiber-mass system through the proposed constitutive model and asymptotic relationship. In Section 5, we further investigate the dynamics of the fiber-mass system with a finite characteristic time based on the constitutive model. Finally, a short summary is provided in Section 6.

Constitutive Model of a Translating Thermally Responsive Fiber
In this section, we propose a constitutive model of a LCE fiber connected with a translating frame in a linear temperature field. For a small characteristic time, its asymptotic relationship is further derived. When the displacements and velocities of the two ends of the LCE fiber are given, the tensional force in the fiber can be calculated.

Constitutive Model
The model for a LCE fiber connected with a translating frame in a linear temperature field is shown in Figure 1. We define a reference state of the fiber, which represents the free-standing fiber without thermal expansion at reference temperature T r , and the original length is denoted by L. w 1 (t), . w 1 (t), w 2 (t), and . w 2 (t) are the displacements and velocities of the translating frame and the free end of the LCE fiber, which is regarded as a current state. To analyze the deformation of a translating LCE fiber, a Lagrangian coordinate system, X, is fixed in the reference configuration of the fiber, and the Eulerian coordinate system, x, is also built in the current state to describe the spatial coordinate. The origins of coordinates are fixed at the O point. Meanwhile, x = x(X, t) is the instantaneous position of a material at point X of the translating fiber, and u(X, t) is its instantaneous displacement.
velocities of the translating frame and the free end of the LCE fiber, which is regarded as a current state. To analyze the deformation of a translating LCE fiber, a Lagrangian coordinate system, X , is fixed in the reference configuration of the fiber, and the Eulerian coordinate system, x , is also built in the current state to describe the spatial coordinate. The origins of coordinates are fixed at the O point. Meanwhile, is the instantaneous position of a material at point X of the translating fiber, and ( ) t X u , is its instantaneous displacement.
are the displacements and velocities of the translating frame and the free end of the LCE fiber. When the displacements and velocities of the two ends of the LCE fiber are given, the tension force of the fiber can be determined.
LCE material can generally generate a rapid deformation response and recoverable deformation, which is synthesized and prepared by liquid crystalline molecules and flexible long-chain polymers [46]. When a tensile force is applied onto the LCE fiber, it is stretched in the monodomain state. When heated up in the temperature field, the stretched fiber contracts along the axial direction because of the transition from the monodomain phase to the isotropic phase. When cooled down, the fiber can completely relax to its original length. Hence, the nematic-isotropic phase transition and the reorientation of the liquid crystal mesogens occur during the contraction and expansion of the fiber. For simplicity, we assume that the tension force, ( ) t F L , of the translating fiber is proportional to its elastic strain, as: is the one-dimensional strain, namely, ε . (2) is the thermally induced strain, and we also assume that it is proportional to the temperature difference, ( ) t X T , , in the fiber between the actual temperature of the fiber and the reference temperature r T , i.e., w 1 (t), w 2 (t), and . w 2 (t) are the displacements and velocities of the translating frame and the free end of the LCE fiber. When the displacements and velocities of the two ends of the LCE fiber are given, the tension force of the fiber can be determined. LCE material can generally generate a rapid deformation response and recoverable deformation, which is synthesized and prepared by liquid crystalline molecules and flexible long-chain polymers [46]. When a tensile force is applied onto the LCE fiber, it is stretched in the monodomain state. When heated up in the temperature field, the stretched fiber contracts along the axial direction because of the transition from the monodomain phase to the isotropic phase. When cooled down, the fiber can completely relax to its original length. Hence, the nematic-isotropic phase transition and the reorientation of the liquid crystal mesogens occur during the contraction and expansion of the fiber. For simplicity, we assume that the tension force, F L (t), of the translating fiber is proportional to its elastic strain, as: in which K is the spring constant. ε(X, t) is the one-dimensional strain, namely, ε T (X, t) is the thermally induced strain, and we also assume that it is proportional to the temperature difference, T(X, t), in the fiber between the actual temperature of the fiber and the reference temperature T r , i.e., with α being the thermal expansion coefficient. α < 0 indicates thermal contraction, while α > 0 indicates thermal expansion.
Considering that in the fiber, F L (t) is uniform and constant, Equation (1) is integrated on both sides from 0 to X. By combining Equations (2) and (3), Equation (1) can further be expressed as: When X = L, the tensional force at the end of the translating fiber can be reduced as: It is noted that the temperature field in the translating fiber varies with time and is inhomogeneous. Hence, there exits heat exchange between the translating fiber and the surrounding environment, and we denote the temperature distribution by T ext (x). We also assume that the thin fiber's radius, R, is significantly smaller than its length, L. As a result, the temperature field in the translating fiber is essentially homogeneous over its radius. The temperature field, T = T(X, t), in the fiber is governed by: where τ is the characteristic time, which is related to the heat capacity of the fiber per unit length, ρ c , and the heat exchange coefficient, h. The characteristic time can be expressed by τ = ρ c /h. The actual external temperature field is generally complex and nonlinear and may be produced by a light on-off or hot plate. In this study, for simplicity, we assume that the external steady temperature field is linear with the spatial coordinate, and the distribution of the linear temperature field can be described by: where β is the temperature gradient. It is noted that the origin of spatial coordinates is fixed at O point (x = 0), and the temperature difference, T ext (0), is zero, i.e., the actual external temperature is equal to the reference temperature T r . It should be noted that the spatial coordinate of the material point, x(X, t), can be calculated by: x(X, t) = u(X, t) + X, where u(X, t) can be derived from Equations (4) and (5) as: By introducing the following dimensionless parameters: t = t/ L/g, F L = F L /mg, u = u/L, w 1 = w 1 /L, w 2 = w 2 /L, X = X/L, x = x/L, τ = τ/ L/g, K = KL/mg, α = αT L , T = T/T L , T ext = T ext /T L , and β = βL/T L (T L is the environmental temperature at x = L), we combine Equations (5)- (9) and derive the formulas as follows: Equation (10) shows that the temperature field inside the fiber can be determined for a given displacement of the two ends of the fiber. Then, the tension force of the fiber can be calculated using Equation (11). In this study, the constitutive model of the translating fiber is considered Equations (10) and (11). The four dimensionless parameters, such as the spring constant K, the temperature gradient β, the characteristic time τ, and the thermal expansion coefficient α, are important influencing factors of the constitutive model. Conveniently, the constitutive model can be used to analyze the translating fiber-mass system in the following sections.

Asymptotic Relationship
In the following, we develop the asymptotic relationship of the translating LCE fiber with a small characteristic time, τ, which means the situations of a high heat transfer coefficient or a low heat capacity. For τ << 1, the temperature field can be written as: Then, we insert Equation (12) into Equation (10) and compare the coefficients of the same order of τ. The temperature field can be further given as: By analytically solving Equation (13), we can obtain the zero-order part of the temperature field as follows: Then, we combine Equation (14) with Equation (15), and the first-order part of the temperature field can be calculated as: Therefore, the temperature field can be solved as: We combine Equation (11) with Equation (17) and obtain the tensional force, F L (t), of the translating fiber as: Equation (18) is the asymptotic relationship of the constitutive model of translating LCE fiber described by Equations (10) and (11). Next, we further establish a dynamic fiber-mass system model and utilize the constitutive model to obtain the corresponding governing equation. For three typical cases of the frame at rest, uniform, and accelerating translation, the asymptotic relationship is utilized in Section 3 to study the dynamical behaviors of the fiber-mass system with zero characteristic time, and both the constitutive model and approximate solution are used in Section 4 to study the case of a small characteristic time. The constitutive model is also utilized to investigate the case of finite characteristic time, τ, in Section 5. The typical values of material properties and geometric parameters from accessible experiments are listed in Table 1 [46,47], and the dimensionless parameters are estimated in Table 2.
ρ c heat capacity of the fiber per unit length 0.01 J/m 2 / • C Table 2. Dimensionless parameters.

Dynamics of the Fiber-Mass System with Zero Characteristic Time
In this section, we establish a thermally responsive fiber-mass system with zero characteristic time, and formulate its corresponding governing equations based on the constitutive model, to investigate its dynamical behaviors for three typical cases of the frame at rest, uniform, and accelerating translation. By obtaining the dynamical behavior of this translating thermally responsive fiber-mass system, we can provide some useful recommendations for the design and motion control of micro-robots of various structures and applications.

Governing Equation
We first study the dynamics of the fiber-mass system translating in a temperature field varying linearly with position for zero characteristic time, τ = 0, which is sketched in Figure 2. One end of the LCE fiber is attached to a mass block with a mass of m, and the other end is connected to a frame that ignores gravity. Therefore, the governing equation for the dynamical model of the mass block can be built as below: where F d . w 2 is the damping force of the mass block due to the damping, mg is the gravitational force of the mass block, and F L (t) is the tensional force from the fiber. We assume that the damping force is linearly proportional to its velocity, which can be written as: where a 0 is the damping coefficient. By defining F d = F d /mg and a 0 = a 0 m L g , and submitting Equations (18) and (20) into Equation (19), the dimensionless form of the governing equation can be obtained as: ..
, the governing equation can be simplified as: It can be seen from Equation (22) that the motion of the fiber-mass system is determined by the displacement, velocity, and acceleration of the mass connected to the end of the LCE fiber, as well as the displacement of the frame attached to the other end of the LCE fiber. With the help of Matlab software, we adopt the four-order Runge-Kutta method described in detail in [47,61] to solve the differential Equation (22) with variable By It can be seen from Equation (22) that the motion of the fiber-mass system is determined by the displacement, velocity, and acceleration of the mass connected to the end of the LCE fiber, as well as the displacement of the frame attached to the other end of the LCE fiber. With the help of Matlab software, we adopt the four-order Runge-Kutta method described in detail in [47,61] to solve the differential Equation (22) with variable coefficients. In this paper, we investigate the vibration behavior of the mass block for three typical cases of the frame at rest, uniform, and accelerating translation.

Static Frame with
.
In the case of a static frame with . w 1 t = 0, Equation (22) can be simplified as: ..
Equation (23) is the standard form of the damped free vibration differential equation. By solving the equation, the solution can be analytically derived as: where the natural angular frequency ω 0 = Kαβ e αβ −1 , δ = a 0 2 , and the parameters are related to the system's initial condition setting, especially,ŵ 2 (0) = x 0 andˆ. w 2 (0) = v 0 . The variation of the displacement of the mass with time for the initial condition of the system setting as w 2 (0) = 0 and . w 2 (0) = 0 is plotted in Figure 3 by using Equation (24). We set the system variables K = 20, α = −0.3, β = 1, and a 0 = 0.01 for the numerical calculations. It is demonstrated that the mass block vibrates up and down when subjected to its own gravity and the tensional force of the fiber. Due to the damping consuming the energy, the amplitude of the vibration reduces dramatically over time and rapidly decays to zero.
where the natural angular frequency x v x δ δ ω θ are related to the system's initial condition setting, especially, . The variation of the displacement of the mass with time for the initial condition of the system setting as is plotted in Figure 3 by using Equation (24). We set the system variables 20 for the numerical calculations. It is demonstrated that the mass block vibrates up and down when subjected to its own gravity and the tensional force of the fiber. Due to the damping consuming the energy, the amplitude of the vibration reduces dramatically over time and rapidly decays to zero.  , the fiber-mass system also translates with the translating frame, as

Uniform Translational Frame with
.
We further investigate the dynamics of a system consisting of a fiber and a mass in the case of a uniform translational frame with . w 1 t = C. By setting the dimensionless parameters K = 10, α = −0.38, β = 1, and a 0 = 0.45, Figure 4 plots the time histories of the displacements and phase trajectories in the cases of . w 1 t = 0 and . w 1 t = 0.01. For . w 1 t = 0, due to the dissipation of damping, the vibration amplitude of the mass block gradually decreases and the mass finally rests at the static mode, as shown in Figure 4a,b. For . w 1 t = 0.01, the fiber-mass system also translates with the translating frame, as shown in Figure 4c,d. However, the mass first vibrates and finally stops vibration due to the damping dissipation. It is noted that the equilibrium position of the mass and the length of the LCE fiber vary with time because of the variation of contraction of the LCE fiber translating in the steady linear temperature field.

Uniformly Accelerated Translational Frame with
.
Then, we analyze the dynamics of the fiber-mass system in the case of the frame at uniformly accelerated translation with . w 1 t = a m t. By setting the dimensionless parameters K = 10, α = −0.38, β = 1, a 0 = 0.45, and the accelerated velocity a m = 0.0001, Figure 5 plots the time history and phase trajectory of the displacement of the fiber-mass system with uniformly accelerated translational frame in a linear temperature field. A similar conclusion is found in the cases of the frame at uniform translation as that shown in Figure 4c,d. The mass first vibrates and finally stops vibration due to the damping dissipation. Meanwhile, whether the frame is at uniform translation or uniformly accelerated translation, the equilibrium position of the mass and the length of the LCE fiber vary with time because of the variation of contraction of the LCE fiber translating in the steady linear temperature field. shown in Figure 4c,d. However, the mass first vibrates and finally stops vibration due to the damping dissipation. It is noted that the equilibrium position of the mass and the length of the LCE fiber vary with time because of the variation of contraction of the LCE fiber translating in the steady linear temperature field.

Dynamics of the Fiber-Mass System with Small Characteristic Time
In this section, we utilize the constitutive model and asymptotic relationship to study the dynamical behaviors of the translating fiber-mass system with small characteristic time τ << 1 in the typical cases of the frame at rest, uniform, and accelerating translation.

Governing Equations
The translating fiber-mass system in a temperature field varying linearly with position with a small characteristic time is represented by Figure 2. It is noted that the system only has two modes for linear damping: converging or diverging [47]. In contrast, the system can continuously vibrate without diverging in practice. [47]. Only the nonlinear damping would be considered in the current study, and it is assumed that: where a 0 and a 1 are the corresponding damping coefficients. With the definition of the dimensionless parameters F d = F d /mg, a 0 = a 0 m L g , and a 1 = a 1 L m , we substitute the Equations (18) and (25) into Equation (19), and the governing equation for the mass block can be rewritten as: ..
It can be seen from Equation (27) that when the damping coefficient is positive, it dissipates energy during movement. When the damping coefficient is negative, it contributes to the motion of the system because energy input from the environment compensates for the damping dissipation [47,61]. Next, we study the dynamics of the fiber-mass system for three typical cases of the frame at rest, uniform, and accelerating translation through the proposed constitutive model and asymptotic relationship. In order to obtain the exact solution, the constitutive model of Equations (11) and (12) is coupled with the governing Equation (19), and then solved by utilizing the fourth-order Runge-Kutta method. In the calculation of the asymptotic solution, the asymptotic relationship of Equation (27) can be numerically calculated to obtain the relationship between the displacement and time of the mass block and to analyze the translating fiber-mass system.

Static Frame with
. w 1 t = 0 Figure 6 plots the variation of the displacement of mass with time in the fiber-mass system with a static frame in a linear temperature field. In the numerical calculations, we set the parameters K = 20, τ = 0.01, α = −0.3, β = 1, a 0 = 0.01, a 1 = 0.05, and the initial state of w 2 (0) = 0 and . w 2 (0) = 0. The blue solid line in Figure 6 shows the asymptotic solution with Equation (27) of the dynamics of the system. Meanwhile, the red dotted line in Figure 6 represents its exact solution of the translating thermally responsive fiber-mass system, which is solved by combining the constitutive model with Equations (11) and (12) and the dynamics governing Equation (19). The comparation between the asymptotic results and the exact numerical results is also presented. It is shown that the two solutions are consistent with each other. As shown in Figure 6, the amplitude of the displacement of the mass block is steady and the mass vibrates periodically in the linear temperature field. This is because during the vibration, the LCE fiber contracts in high temperature and relaxes in low temperature, which leads to continuous stretch or shrink of the LCE fiber and eventually achieves a periodic self-sustained oscillation. During the process of the self-oscillation, the fiber can absorb heat energy from its surrounding environment to sustain its motion due to the compensation for damping dissipation. of the mass block is steady and the mass vibrates periodically in the linear temperature field. This is because during the vibration, the LCE fiber contracts in high temperature and relaxes in low temperature, which leads to continuous stretch or shrink of the LCE fiber and eventually achieves a periodic self-sustained oscillation. During the process o the self-oscillation, the fiber can absorb heat energy from its surrounding environment to sustain its motion due to the compensation for damping dissipation. Here, we take temperature gradient 1 = β to describe the linear temperature field distribution, which indicates that the temperature increases from top to bottom. During the movement of the fiber, the LCE fiber contracts in high temperature and relaxes in low temperature, which leads to energy absorption and compensates for damping dissipation Consequently, the mass block is in the periodic self-excited motion. To further analyze the influence of the temperature gradient, we set the temperature gradient 1 β = − , which indicates that the temperature decreases from top to bottom. The relationship of the dis placement and the time of the mass in the fiber-mass system is shown in Figure 7. It i found that the mass block no longer moves periodically and becomes stationary due to energy dissipation, which is much different from the self-oscillation for the case of  (27), the damp ing coefficient is negative, and the fiber engine module can absorb energy from the envi ronment and compensate for the damping dissipation like an engine during the move ment, which contributes to the self-oscillation of the system. Here, we take temperature gradient β = 1 to describe the linear temperature field distribution, which indicates that the temperature increases from top to bottom. During the movement of the fiber, the LCE fiber contracts in high temperature and relaxes in low temperature, which leads to energy absorption and compensates for damping dissipation. Consequently, the mass block is in the periodic self-excited motion. To further analyze the influence of the temperature gradient, we set the temperature gradient β = −1, which indicates that the temperature decreases from top to bottom. The relationship of the displacement and the time of the mass in the fiber-mass system is shown in Figure 7. It is found that the mass block no longer moves periodically and becomes stationary due to energy dissipation, which is much different from the self-oscillation for the case of β = 1. The different behaviors of the dynamic fiber can be understood by the proposed theoretical model. For the case of β = −1, from Equation (27), the damping coefficient is positive, and the fiber engine module represents the damper that dissipates energy during movement and leads to a stationary mass. For the case of β = 1, from Equation (27), the damping coefficient is negative, and the fiber engine module can absorb energy from the environment and compensate for the damping dissipation like an engine during the movement, which contributes to the self-oscillation of the system. . Figure 8a,b show that the LCE fiber-mass system first vibrates and eventually develops into periodic self-sustained oscillation for the frame at rest. For the translating frame with , the fiber-mass system also translates with the frame, as shown in Figure 8c,d. However, the mass block can also self-oscillate due to the energy compensation between energy input and damping dissipation. Similarly, the equilibrium position of the mass and the length of the LCE fiber also vary with time because of the variation of contraction of the LCE fiber translating in the steady linear Figure 7. Variation of the displacement of the mass block with time in a linear temperature field with temperature gradient β = −1. It is found that the mass block develops into stationary due to energy dissipation, which is much different from the self-oscillation for the case of β = 1.  Figure 8a,b show that the LCE fiber-mass system first vibrates and eventually develops into periodic self-sustained oscillation for the frame at rest. For the translating frame with . w 1 t = 0.005, the fiber-mass system also translates with the frame, as shown in Figure 8c,d. However, the mass block can also self-oscillate due to the energy compensation between energy input and damping dissipation. Similarly, the equilibrium position of the mass and the length of the LCE fiber also vary with time because of the variation of contraction of the LCE fiber translating in the steady linear temperature field.

Uniformly Translational Frame with
Next, we consider the dynamics of the translating LCE fiber in the case of a uniform translational frame with . By setting the dimensionless parameters . Figure 8a,b show that the LCE fiber-mass system first vibrates and eventually develops into periodic self-sustained oscillation for the frame at rest. For the translating frame with , the fiber-mass system also translates with the frame, as shown in Figure 8c,d. However, the mass block can also self-oscillate due to the energy compensation between energy input and damping dissipation. Similarly, the equilibrium position of the mass and the length of the LCE fiber also vary with time because of the variation of contraction of the LCE fiber translating in the steady linear temperature field.

Uniformly Accelerated Translational Frame with
. w 1 t = a m t Next, we investigate the case of the translating LCE fiber-mass system in the case of a frame at uniformly accelerated translation with . w 1 t = a m t. By setting the dimensionless parameters K = 20, τ = 0.01, α = −0.3, β = 1, a 0 = 0.02, a 1 = 0.05, and the accelerated velocity a m = 0.0001, the time histories of displacement and phase trajectories of self-oscillation are plotted in Figure 9. For the case of the frame at uniformly accelerated translation, the mass block can also self-oscillate periodically due to the energy compensation between energy input and damping dissipation. Meanwhile, the equilibrium position of the mass and the length of the LCE fiber both vary with time. This is because the LCE fiber translates in the steady linear temperature field and its contraction varies with time. These results are similar to the cases of the frame at uniform translation, as shown in Figure 8c,d.  Figure 9. For the case of the frame at uniformly accelerated translation, the mass block can also self-oscillate periodically due to the energy compensation between energy input and damping dissipation. Meanwhile, the equilibrium position of the mass and the length of the LCE fiber both vary with time. This is because the LCE fiber translates in the steady linear temperature field and its contraction varies with time. These results are similar to the cases of the frame at uniform translation, as shown in Figure 8c,d.

Dynamics of the Fiber-Mass System with Finite Characteristic Time
The asymptotic relationship of Equation (27) and the constitutive model of Equations (11) and (12) have been used to analyze the dynamical behaviors of the translating LCE

Dynamics of the Fiber-Mass System with Finite Characteristic Time
The asymptotic relationship of Equation (27) and the constitutive model of Equations (11) and (12) have been used to analyze the dynamical behaviors of the translating LCE fiber-mass system with zero τ or small τ in a temperature field varying linearly with position in the above study. In this section, by considering the different finite characteristic times, we study the dynamical behaviors of the translating LCE fiber-mass system in the typical cases of the frame at rest, uniform, and accelerating translation. Since the characteristic time τ is not small, the asymptotic relationship is no longer applicable, and we will only adopt the constitutive model for the numerical calculations. Figures 10-12 plot the time histories of the displacement of mass in the fiber-mass system translating in a linear temperature field for different finite characteristic times τ. Three typical cases of the frame at rest, uniform, and accelerating translation are considered in the numerical calculations, by setting the dimensionless parameters K = 20, α = −0.3, β = 1, and a 0 = 0.02, as well as a 1 = 0.05. For the fiber-mass system with a static frame in a linear temperature field, when the characteristic time τ increases, the equilibrium position and frequency are unaffected, while the amplitude first increases and then decreases obviously, as shown in Figure 10. Hence, there exists an optimal characteristic time for the maximum of amplitude, which can provide convenience and guidance for engineering applications. Similarly, the conclusion can also be found in the case of the frame at uniform translation, as shown in Figure 11, and the case of the frame at uniform accelerating translation, as shown in Figure 12.

Conclusions
Self-oscillation of a LCE fiber translating in an inhomogeneous temperature field is worth investigating to widen its applications. In this paper, we proposed a theoretic con- For a given finite characteristic time, by comparing Figure 10a, Figure 11a, and Figure 12a of the three cases of the frame at rest, uniform, and accelerating translation, the fiber-mass system remains in the period of self-excited oscillation by the energy input from the ambient linear temperature field to compensate for damping dissipation. Meanwhile, for the frame at uniform translation or uniformly accelerated translation, the equilibrium position of the mass varies with time, which is different from the case of the frame at rest. A similar conclusion is also found in Figure 10b, Figure 11b, and Figures 12b and 10c,  Figure 11c, and Figure 12c for different finite characteristic times. These results mean that the finite characteristic time and the frame translation do not affect the self-oscillation.

Conclusions
Self-oscillation of a LCE fiber translating in an inhomogeneous temperature field is worth investigating to widen its applications. In this paper, we proposed a theoretic constitutive model and the asymptotic relationship of a LCE fiber translating in a temperature field varying linearly with position and investigated the dynamical behaviors of a fiber-mass system translating in a linear temperature field. The main conclusions are summarized as follows: (i) For zero characteristic time, the LCE fiber-mass system in the linear temperature field vibrates freely and cannot develop into self-oscillation in the three typical cases of the frame at rest, uniform, and accelerating translation. (ii) For a small characteristic time, the LCE fiber-mass system in the linear temperature field can self-oscillate by absorbing energy from the surrounding temperature field to compensate for the damping dissipation. (iii) For a finite characteristic time, the amplitude of the selfoscillation of the LCE fiber-mass system increases first, then decreases as the characteristic time increases, and at an optimal characteristic time, the amplitude can attain the maximum value. (iv) In the three cases of the frame at rest, uniform, and accelerating translation, the equilibrium position of the self-oscillation is different, while its amplitude and frequency are the same. The results are expected to provide some useful recommendations for the design and motion control of micro-robots, motors, and active machines with various structures and applications.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that support the findings of this study are available upon reasonable request from the authors.

Conflicts of Interest:
The authors declare no conflict of interest. Nomenclature a 0 first damping coefficient of the mass block a 1 second damping coefficient of the mass block a m accelerated velocity of the translational frame A initial amplitude of the free vibration of the mass block C uniform velocity of the translational frame F d . w 2 damping force of the mass block F L (t) tensional force of the translating LCE fiber g gravitational acceleration h heat exchange coefficient between the LCE fiber and the environment K spring constant of the LCE fiber L original length of the LCE fiber m mass of the mass block R radius of the thin LCE fiber t time T(X, t) temperature field of the LCE fiber T ext (x) environment temperature field T L environmental temperature at x = L T r reference temperature in reference state u(X, t) instantaneous displacement of a material point w 1 (t) displacement of the translating frame . w 1 (t) velocity of the translating frame w 2 (t) displacement of the free end of the LCE fiber . w 2 (t) velocity of the free end of the LCE fiber X Lagrangian coordinate of the LCE fiber