# Self-Sustained Oscillation of a Photothermal-Responsive Pendulum under Steady Illumination

Authors: Dali Ge; Peibao Xu; Kai Li
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Self-sustained oscillation has the advantages of harvesting energy from the environment and self-control, and thus, the development of new self-oscillating systems can greatly expand its applications in active machines. In this paper, based on conventional photothermal shrinkable material or photothermal expansive material, a simple pendulum is proposed. The light-powered self-sustained oscillation of the simple pendulum is theoretically studied by establishing a dynamic model of the photothermal-responsive pendulum. The results show that there are two motion modes of the simple pendulum, which are the static mode and the oscillation mode. Based on the photothermal-responsive model, this paper elucidates the mechanism of the self-excited oscillation. The condition for triggering self-excited oscillation is further studied. In addition, the influence of the system parameters on the amplitude and frequency is also obtained. This study may have potential applications in energy harvesting, signal monitoring, and soft machines.


---

## Body

## 1. Introduction

Photothermal-responsive materials exhibit photothermal effects under light illumination, which results in some macroscopic changes, such as shape change and color change. Due to the advantages of light stimulus, such as environmentally friendly, remote controllability, and instantaneity, photothermal-responsive materials have attracted many researchers' attention. Suzuki and Tanaka [1] introduced the photosensitive component chlorophyll in the gel and designed a light-sensitive gel. The temperature of the illuminated gel increases and the gel shrinks. Hu et al. [2] found that PNIPAM gel will reduce the transmittance of visible laser light when infrared radiation makes the gel locally heated up. Kubo et al. [3] filled silicon dioxide inverse opal with liquid crystal molecules and made a device that changes the color of reflected light caused by the light response.Photothermal-responsive materials with reasonable design can produce photo-induced deformation and motion for a variety of flexible intelligent actuators, which have wide application prospects in the fields of artificial muscles [4], retractable optical devices [5], light-driven microbionic propellers [6–8], light-powered microfluidic pump [9], light-fueled microgenerator [10], and so on. Ikeda et al. [11] used alternating ultraviolet and visible light to achieve cis-trans isomer photomechanical actuation of liquid crystal elastomers containing azobenzene. Carbon nanotubes, graphene, and graphene oxide have good photothermal effects, which can convert near-infrared light and visible light into heat [12–15]. They can be used as nanoheaters to generate local heat and cause conformational modification or phase transition of elastomer matrix [15–18] to realize light-powered motions.Self-excited oscillations have been observed in many fields of science and engineering [19–23] and can be utilized to produce periodic motions which are often required in numerous engineering applications [24–28]. Different from the forced vibration with damping dissipation, which often applies periodic stimulus to cause periodic vibration, the self-excited oscillation utilizes the nonlinear feedback within the system to result in and maintain steady-state periodic motions, such as the flapping of leaves under the breeze, the peristalsis of the digestive tract, the periodic beating of the heart, and the beautiful music produced by the violin bow. Under the steady environment, the self-excited oscillation system takes its own motion state as the regulator to balance the input energy and the dissipated energy, so as to realize the continuous periodic motion. The characteristics of the self-excited oscillation system, such as frequency and amplitude, are often determined by the physical and geometrical parameters of the system and have nothing to do with the initial conditions. The particularity of the self-excited oscillation arouses the interest of researchers and has been applied in robots and other fields. Xu et al. [24] have carried out dynamic simulation for a woodpecker robot running under friction and collision. Ono et al. [25] studied the mechanism of achieving steady walking motion of a self-oscillating bipedal knee robot. Morazzani et al. [26] proposed a new type of self-excited three-legged walking robot. According to the principle of bionics, Li et al. [27] proposed a pipeline crawling microrobot based on the self-excited oscillation. Yamano et al. [28] constructed a multilink adaptive control swimming robot model, which can generate different self-excited oscillations according to the change of environment, so as to change the motion mode.Recently, the self-excited oscillation of optomechanical oscillators has been extensively explored and applied in macrodevices [29], microdevices [30], and nanodevices [31]. There exist various self-excited periodic motions in systems dominated by photothermal effects [32, 33]. With the in-depth study of photothermal-responsive materials, more progress has been made in the study of self-excited oscillation based on photothermal-responsive materials. Hollander et al. [34] studied a nonlinear thermoviscoelastic cantilever coated with a light-absorbing film. The self-excited oscillation of the cantilever is triggered by the photothermal effect. Dirk J. Gelebart’s research group [35] added light-absorbing dyes to liquid crystal network polymers and prepared a continuous oscillating composite film by the photothermal effect. Houri et al. [36] irradiate the graphene barrel with a red continuous wave laser beam, and it oscillates autonomously by photothermal feedback. Buks and Martin [37] studied a fully on-fiber optomechanical cavity which was formed by patterning a suspended metallic mirror near the tip of an optical fiber and by introducing a static reflector inside the fiber and characterized its performance as a sensor. Optically induced self-excited oscillation is observed above a threshold value of the injected laser power.Although some self-excited oscillations of photothermal-responsive materials and their applications have been studied, the modes of self-excited oscillation are not enough, which limits the potential applications of self-excited oscillation. In this paper, a simple pendulum is constructed by using conventional photothermal expansive material (PEM) or photothermal shrinkable material (PSM) as a cycloid, and the governing equation of variable pendulum length is established based on the photothermal effect. Through numerical calculation, it is verified that the photothermal-responsive pendulum can realize self-excited oscillation under steady light illumination. The mechanism of self-excited oscillation is elucidated, the conditions for realizing self-excited oscillation are discussed, and the dependence of amplitude and frequency on system parameters is extensively studied.


## 2. Theoretical Model and Formulation

### 2.1. Governing Equations for Dynamics of the Simple Pendulum

The dynamic model and governing equations of a photothermal-responsive pendulum under steady illumination are established in this section. As shown in Figure1, a simple pendulum made of photothermal-responsive materials is established and illuminated steadily. The end of the cycloid made of photothermal-responsive material with negligible mass is tied to a fixed point O. A pendulum ball with mass m is suspended by the cycloid. By establishing a polar coordinate system, the position of a simple pendulum can be expressed by pendulum length r and swing angle θ. During the motion of the pendulum, the length r and swing angle θ of the pendulum both vary with the time. In the polar coordinate system, two mutually perpendicular unit vectors er and eθ are along the directions of increasing pendulum length and pendulum angle, respectively, and they also change during the swing of the pendulum.Figure 1
Schematics of the light-driven swing of a photothermal-responsive pendulum based on (a) PSM cycloid and (b) PEM cycloid.
(a)(b)To obtain the self-excited motion, a PSM pendulum should be illuminated in the middle zone as shown in Figure1(a), while a PSE pendulum should be illuminated on the two sides as shown in Figure 1(b), where the angle φ denotes the illuminated zone or dark zone in the middle. For the two pendulum systems in Figure 1, the cycloid gradually shrinks when it is in the middle zone, while gradually expands when it swings out of the middle zone. Eventually, the pendulum could swing continuously as discussed as follows.In the polar coordinate system, the velocity and acceleration of the pendulum along the directionseθ and er are as follows:(1)vθ=rθ˙,aθ=rθ¨+2r˙θ˙,vr=r˙,ar=r¨−rθ˙2,where r˙ and r¨ are the first derivative dr/dt and the second derivative d2r/dt2, respectively, and θ˙ and θ¨ are the first derivative dθ/dt and the second derivative d2θ/dt2, respectively.As shown in Figure1, mg is the gravity of the pendulum, F is the tension of the cycloid, and Fθ and Fr are the components of the air damping force in eθ and er directions, respectively. Since the fiber is generally much thinner than the pendulum ball, the air damping force of the fiber is ignored for simplicity. It is assumed that the air damping force of the pendulum is proportional to the velocity of the ball, and the direction is opposite to the velocity of the ball. The governing equation for the dynamic of the pendulum can be given as(2)mrθ¨+2r˙θ˙=−mgsinθ−βθ˙r,mr¨−rθ˙2=mgcosθ−βr˙−F,where β is the linear air damping coefficient. To study the light-powered motion of a simple pendulum with the photothermal response and its mechanism, equation (2) is rewritten as(3)θ¨=−2r˙r+βmθ˙−grsinθ,F=mgcosθ−βr˙−mr¨−rθ˙2.


### 2.2. Model of the Photothermal-Responsive Pendulum

This section mainly describes the dynamics of temperature and length of the photothermal-responsive pendulum under stable light illumination and in dark. Since the simple pendulums made with PSM and PEM in Figure1 have the same photothermal effect, in the following, the simple pendulum of PSM in Figure 1(a) is chosen as the representative of the photothermal-responsive pendulum.The cycloid of the pendulum has a very small radius, i.e., the well-known Biot number is very small. Therefore, heat exchange in the pendulum is assumed to be very fast and the temperature in the pendulum is homogeneous. Due to the photothermal effect, the pendulum can convert light illumination into heat. The heat converted from light per second is denoted byq. The pendulum can also exchange heat with the environment, and the heat flux is assumed to be linear to the temperature difference T between the pendulum and the environment. Under light illumination, the temperature difference T is determined by(4)T˙=q−kTρc,where ρc is the specific heat capacity and k is the heat transfer coefficient. By solving equation (4), in the illumination zone, the temperature difference T is(5)T=T01−e−t/τ,while in the dark zone, i.e., q=0, the temperature difference T is(6)T=T0e−t/τ,where T0=q/k represents the limit temperature difference of photothermal-responsive fiber under longtime illumination, τ=ρc/k reflects the characteristic time for heat exchange between photothermal-responsive fiber and environment, and the larger τ indicates the longer time required for attaining the limit temperature difference T0 of the photothermal-responsive fiber.Deformation of the photothermal-responsive pendulum can be driven by the temperature change, and the thermal strain of the fiber is assumed to be linear to the temperature change. The elastic modulus is assumed to be very large, and the elastic strain is ignored. Therefore, the length of the PSM fiber can be calculated as, in the illumination zone:(7)r=r01−αT01−e−t/τ,in the dark zone: (8)r=r01−αT0e−t/τ,where r0 is the origin length of the PSM fiber without thermal strain and α is the linear thermal contraction coefficient of the PSM fiber.


### 2.3. Nondimensionalization

To nondimensionalize the governing equations above, we introduce the following dimensionless parameters:t¯=t/τ, T¯=T/T0, r¯=r/r0, ε0=αT0, β¯=βτ/m, g¯=gτ2/r0, and F¯=Fτ2/mr0. Here, ε0 is the contraction strain of the photothermal-responsive cycloid at the limit temperature difference T0, and β¯ reflects the air damping. g¯ can be rewritten as g¯=τ/τ02 by the natural period τ0=r0/g of the single pendulum, so it represents the speed of heat transfer relative to natural oscillation. The larger g¯ is, the faster the heat transfers. The dimensionless parameter F¯ characterizes the tension of the photothermal-responsive cycloid. Then, equations (5)–(8) can also be rewritten as, in the illumination zone, i.e., −φ≤θ≤φ.(9)T¯=1−e−t¯,r¯=1−ε01−e−t¯.in the dark zone, i.e., θ<−φ or θ>φ:(10)T¯=e−t¯,r¯=1−ε0e−t¯.The following differential equation is introduced:(11)r˙=drdt=1τdrdt¯=r¯˙τ,r¨=d2rdt2=r¯¨τ2.By inserting equations (9)–(11), equation (3) can be derived as, in the illumination zone, i.e., −φ≤θ≤φ:(12)d2θdt¯2=2ε0et¯1−ε01−et¯−β¯dθdt¯−g¯1−ε01−et¯sinθ,F¯=g¯cosθ+β¯ε0e−t¯−ε0e−t¯+1−ε01−e−t¯dθdt¯2,in the dark zone, i.e., θ<−φ or θ>φ:(13)d2θdt¯2=−2ε0et¯1−ε0et¯−β¯dθdt¯−g¯1−ε0et¯sinθ,F¯=g¯cosθ−β¯ε0e−t¯+ε0e−t¯+1−ε0e−t¯dθdt¯2.


### 2.4. Solution Method

The second-order differential equations (12) and (13) provide a complete description of the dynamics of the photothermal-responsive pendulum. The swing angle θ can be obtained by the first equation in equation (13), and then, the pendulum force F¯ can be further obtained by the second equation in equation (13). It is very challenging to obtain an analytical solution to equation (13). In this paper, we will numerically solve the swing of the pendulum based on the well-known fourth-order Runge–Kutta method. By letting ψ=θ˙, the second-order ordinary differential equation (13) is transformed into a first-order differential equation system as follows:(14)dθdt¯=ψ,in the illumination zone, i.e., −φ≤θ≤φ:(15)dψdt¯=2ε0et¯1−ε01−et¯−β¯ψ−g¯1−ε01−et¯sinθ,in the dark zone, i.e., θ<−φ or θ>φ:(16)dψdt¯=−2ε0et¯1−ε0et¯−β¯ψ−g¯1−ε0et¯sinθ.The initial conditions areθ|t¯=0=θ0 and ψ|t¯=0=θ˙0. Given the dimensionless parameters φ, β¯, ε0, and g¯, the solution to equations (14)–(16) can be obtained numerically by programming in software Matlab. When the angle θt¯ and the angular velocity θ˙t¯ of the simple pendulum are obtained, T¯t¯, r¯t¯, and F¯t¯ can be further calculated.


## 2.1. Governing Equations for Dynamics of the Simple Pendulum

The dynamic model and governing equations of a photothermal-responsive pendulum under steady illumination are established in this section. As shown in Figure1, a simple pendulum made of photothermal-responsive materials is established and illuminated steadily. The end of the cycloid made of photothermal-responsive material with negligible mass is tied to a fixed point O. A pendulum ball with mass m is suspended by the cycloid. By establishing a polar coordinate system, the position of a simple pendulum can be expressed by pendulum length r and swing angle θ. During the motion of the pendulum, the length r and swing angle θ of the pendulum both vary with the time. In the polar coordinate system, two mutually perpendicular unit vectors er and eθ are along the directions of increasing pendulum length and pendulum angle, respectively, and they also change during the swing of the pendulum.Figure 1
Schematics of the light-driven swing of a photothermal-responsive pendulum based on (a) PSM cycloid and (b) PEM cycloid.
(a)(b)To obtain the self-excited motion, a PSM pendulum should be illuminated in the middle zone as shown in Figure1(a), while a PSE pendulum should be illuminated on the two sides as shown in Figure 1(b), where the angle φ denotes the illuminated zone or dark zone in the middle. For the two pendulum systems in Figure 1, the cycloid gradually shrinks when it is in the middle zone, while gradually expands when it swings out of the middle zone. Eventually, the pendulum could swing continuously as discussed as follows.In the polar coordinate system, the velocity and acceleration of the pendulum along the directionseθ and er are as follows:(1)vθ=rθ˙,aθ=rθ¨+2r˙θ˙,vr=r˙,ar=r¨−rθ˙2,where r˙ and r¨ are the first derivative dr/dt and the second derivative d2r/dt2, respectively, and θ˙ and θ¨ are the first derivative dθ/dt and the second derivative d2θ/dt2, respectively.As shown in Figure1, mg is the gravity of the pendulum, F is the tension of the cycloid, and Fθ and Fr are the components of the air damping force in eθ and er directions, respectively. Since the fiber is generally much thinner than the pendulum ball, the air damping force of the fiber is ignored for simplicity. It is assumed that the air damping force of the pendulum is proportional to the velocity of the ball, and the direction is opposite to the velocity of the ball. The governing equation for the dynamic of the pendulum can be given as(2)mrθ¨+2r˙θ˙=−mgsinθ−βθ˙r,mr¨−rθ˙2=mgcosθ−βr˙−F,where β is the linear air damping coefficient. To study the light-powered motion of a simple pendulum with the photothermal response and its mechanism, equation (2) is rewritten as(3)θ¨=−2r˙r+βmθ˙−grsinθ,F=mgcosθ−βr˙−mr¨−rθ˙2.


## 2.2. Model of the Photothermal-Responsive Pendulum

This section mainly describes the dynamics of temperature and length of the photothermal-responsive pendulum under stable light illumination and in dark. Since the simple pendulums made with PSM and PEM in Figure1 have the same photothermal effect, in the following, the simple pendulum of PSM in Figure 1(a) is chosen as the representative of the photothermal-responsive pendulum.The cycloid of the pendulum has a very small radius, i.e., the well-known Biot number is very small. Therefore, heat exchange in the pendulum is assumed to be very fast and the temperature in the pendulum is homogeneous. Due to the photothermal effect, the pendulum can convert light illumination into heat. The heat converted from light per second is denoted byq. The pendulum can also exchange heat with the environment, and the heat flux is assumed to be linear to the temperature difference T between the pendulum and the environment. Under light illumination, the temperature difference T is determined by(4)T˙=q−kTρc,where ρc is the specific heat capacity and k is the heat transfer coefficient. By solving equation (4), in the illumination zone, the temperature difference T is(5)T=T01−e−t/τ,while in the dark zone, i.e., q=0, the temperature difference T is(6)T=T0e−t/τ,where T0=q/k represents the limit temperature difference of photothermal-responsive fiber under longtime illumination, τ=ρc/k reflects the characteristic time for heat exchange between photothermal-responsive fiber and environment, and the larger τ indicates the longer time required for attaining the limit temperature difference T0 of the photothermal-responsive fiber.Deformation of the photothermal-responsive pendulum can be driven by the temperature change, and the thermal strain of the fiber is assumed to be linear to the temperature change. The elastic modulus is assumed to be very large, and the elastic strain is ignored. Therefore, the length of the PSM fiber can be calculated as, in the illumination zone:(7)r=r01−αT01−e−t/τ,in the dark zone: (8)r=r01−αT0e−t/τ,where r0 is the origin length of the PSM fiber without thermal strain and α is the linear thermal contraction coefficient of the PSM fiber.


## 2.3. Nondimensionalization

To nondimensionalize the governing equations above, we introduce the following dimensionless parameters:t¯=t/τ, T¯=T/T0, r¯=r/r0, ε0=αT0, β¯=βτ/m, g¯=gτ2/r0, and F¯=Fτ2/mr0. Here, ε0 is the contraction strain of the photothermal-responsive cycloid at the limit temperature difference T0, and β¯ reflects the air damping. g¯ can be rewritten as g¯=τ/τ02 by the natural period τ0=r0/g of the single pendulum, so it represents the speed of heat transfer relative to natural oscillation. The larger g¯ is, the faster the heat transfers. The dimensionless parameter F¯ characterizes the tension of the photothermal-responsive cycloid. Then, equations (5)–(8) can also be rewritten as, in the illumination zone, i.e., −φ≤θ≤φ.(9)T¯=1−e−t¯,r¯=1−ε01−e−t¯.in the dark zone, i.e., θ<−φ or θ>φ:(10)T¯=e−t¯,r¯=1−ε0e−t¯.The following differential equation is introduced:(11)r˙=drdt=1τdrdt¯=r¯˙τ,r¨=d2rdt2=r¯¨τ2.By inserting equations (9)–(11), equation (3) can be derived as, in the illumination zone, i.e., −φ≤θ≤φ:(12)d2θdt¯2=2ε0et¯1−ε01−et¯−β¯dθdt¯−g¯1−ε01−et¯sinθ,F¯=g¯cosθ+β¯ε0e−t¯−ε0e−t¯+1−ε01−e−t¯dθdt¯2,in the dark zone, i.e., θ<−φ or θ>φ:(13)d2θdt¯2=−2ε0et¯1−ε0et¯−β¯dθdt¯−g¯1−ε0et¯sinθ,F¯=g¯cosθ−β¯ε0e−t¯+ε0e−t¯+1−ε0e−t¯dθdt¯2.


## 2.4. Solution Method

The second-order differential equations (12) and (13) provide a complete description of the dynamics of the photothermal-responsive pendulum. The swing angle θ can be obtained by the first equation in equation (13), and then, the pendulum force F¯ can be further obtained by the second equation in equation (13). It is very challenging to obtain an analytical solution to equation (13). In this paper, we will numerically solve the swing of the pendulum based on the well-known fourth-order Runge–Kutta method. By letting ψ=θ˙, the second-order ordinary differential equation (13) is transformed into a first-order differential equation system as follows:(14)dθdt¯=ψ,in the illumination zone, i.e., −φ≤θ≤φ:(15)dψdt¯=2ε0et¯1−ε01−et¯−β¯ψ−g¯1−ε01−et¯sinθ,in the dark zone, i.e., θ<−φ or θ>φ:(16)dψdt¯=−2ε0et¯1−ε0et¯−β¯ψ−g¯1−ε0et¯sinθ.The initial conditions areθ|t¯=0=θ0 and ψ|t¯=0=θ˙0. Given the dimensionless parameters φ, β¯, ε0, and g¯, the solution to equations (14)–(16) can be obtained numerically by programming in software Matlab. When the angle θt¯ and the angular velocity θ˙t¯ of the simple pendulum are obtained, T¯t¯, r¯t¯, and F¯t¯ can be further calculated.


## 3. Two Motion Modes and Mechanism of the Self-Excited Oscillation

### 3.1. Two Motion Modes

By setting the parametersθ0, θ˙0, φ, β¯, ε0, and g¯, the time histories and phase diagrams of the light-powered swing of the pendulum can be obtained. The numerical calculation shows that there exist two motion modes of the simple pendulum: static mode and oscillation mode, as shown in Figure 2. In the computation, we set θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. For ε0=0, the pendulum with initial angular velocity first swings, and then, the amplitude of swing gradually decreases under the air damping, as shown in Figures 2(a) and 2(b). Eventually, the pendulum comes to cease at the middle, which is named as static mode. Although the pendulum is illuminated and the cycloid temperature rises, the pendulum length for ε0=0 will not change and cannot convert light energy into kinetic energy. For ε0=0.3, the pendulum under light illumination can swing continuously and finally develops into self-sustained oscillation, which is named as oscillation mode (Figures 2(c) and 2(d)). Although there is damping in the system, the energy compensation from the light maintains the oscillation. In the following section, we further explore the mechanism of self-excited oscillation and in detail elucidate how the light energy converts into kinetic energy.Figure 2
Time histories and phase diagrams of two motion modes of the pendulum: (a, b) the static mode forε0=0; (c, d) the oscillation mode for ε0=0.3.
(a)(b)(c)(d)


### 3.2. Mechanism of the Self-Excited Oscillation

To explore the mechanism of self-excited oscillation, the evolution laws of the temperature and length of the cycloid are first studied. Based on equations (9) and (10), Figure 3 plots the evolution law of T¯ and r¯ during swing in the illumination zone and the dark zone. As shown in Figure 3(a), the temperature of cycloid increases under continuous illumination, and the dimensionless temperature difference parameter T¯ can gradually increase from 0 to nearly 1. In the dark zone, T¯ can gradually decrease from the maximum value of 1 to 0. When the pendulum passes through the interface of the illumination and dark zones, the temperature of the cycloid will switch to obeying the other evolution law, i.e., from point A to point B or from point C to point D as shown in Figure 3(a). Similarly, the evolution law of the dimensionless parameter r¯ for ε0=0.3 in Figure 3(b) shows that the pendulum length r¯ contracts into 1−ε0 in the illumination zone, while recovers gradually to 1 in the dark zone. The pendulum length also switches when the pendulum passes through the interface between the illumination and dark zones.Figure 3
Evolution laws of (a)T¯ and (b) r¯ of the PSM pendulum.
(a)(b)Based on the evolution law in Figure3, the mechanism of the self-excited oscillation in Figure 2(b) can be elucidated by investigating the time histories of T¯ and r¯ in detail. Figure 4(a) plots the dependence of the dimensionless temperature difference parameter T¯ on the swing angle θ. For −0.1≤θ≤0.1, the cycloid is in the illumination zone and T¯ increases. When the pendulum swings into the dark zone, T¯ decreases. Similarly, Figure 4(b) plots the dependence of r¯ on θ. The length r¯ decreases in the illumination zone, while it increases in the dark zone. Further, Figure 4(c) plots the dependence of F¯ on θ, where F¯ is discontinuous at some points of θ. The tension F¯ sharply increases at the moment of swinging into the dark region (point ②) from the illumination region (point ①). This is because the temperature of the pendulum swinging into the dark zone decreases and in turn the cycloid length increases. As the pendulum swings in the dark area, F¯ continuously varies from point ② to ③. Then, the tension F¯ sharply decreases from point ③ to point ④ when the pendulum reenters into the illumination zone. Eventually, tension F¯ and swing angle θ form a symmetrical closed loop according to the trajectory cycle of ①②③④⑤⑥⑦⑧① in a swing period.Figure 4
Mechanism of self-excited oscillation of the PSM pendulum, forθ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. (a) T¯ vs. θ. (b) r¯ vs. θ. (c) F¯ vs. θ. (d) F¯ vs. r¯.
(a)(b)(c)(d)Based on Figures4(a)–4(c), Figure 4(d) plots the dependence of F¯ with r¯, where the shadow area enclosed by curve of F¯ with r¯ represents the magnitude of the work done by the tension F¯ in a half swing period of the simple pendulum. In an oscillation period, the work done by F¯ is twice shadow area in Figure 4(d), which is the source of self-excited oscillation. When the net work is enough for compensating the dissipation by air damping, the pendulum can oscillate continuously.


## 3.1. Two Motion Modes

By setting the parametersθ0, θ˙0, φ, β¯, ε0, and g¯, the time histories and phase diagrams of the light-powered swing of the pendulum can be obtained. The numerical calculation shows that there exist two motion modes of the simple pendulum: static mode and oscillation mode, as shown in Figure 2. In the computation, we set θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. For ε0=0, the pendulum with initial angular velocity first swings, and then, the amplitude of swing gradually decreases under the air damping, as shown in Figures 2(a) and 2(b). Eventually, the pendulum comes to cease at the middle, which is named as static mode. Although the pendulum is illuminated and the cycloid temperature rises, the pendulum length for ε0=0 will not change and cannot convert light energy into kinetic energy. For ε0=0.3, the pendulum under light illumination can swing continuously and finally develops into self-sustained oscillation, which is named as oscillation mode (Figures 2(c) and 2(d)). Although there is damping in the system, the energy compensation from the light maintains the oscillation. In the following section, we further explore the mechanism of self-excited oscillation and in detail elucidate how the light energy converts into kinetic energy.Figure 2
Time histories and phase diagrams of two motion modes of the pendulum: (a, b) the static mode forε0=0; (c, d) the oscillation mode for ε0=0.3.
(a)(b)(c)(d)


## 3.2. Mechanism of the Self-Excited Oscillation

To explore the mechanism of self-excited oscillation, the evolution laws of the temperature and length of the cycloid are first studied. Based on equations (9) and (10), Figure 3 plots the evolution law of T¯ and r¯ during swing in the illumination zone and the dark zone. As shown in Figure 3(a), the temperature of cycloid increases under continuous illumination, and the dimensionless temperature difference parameter T¯ can gradually increase from 0 to nearly 1. In the dark zone, T¯ can gradually decrease from the maximum value of 1 to 0. When the pendulum passes through the interface of the illumination and dark zones, the temperature of the cycloid will switch to obeying the other evolution law, i.e., from point A to point B or from point C to point D as shown in Figure 3(a). Similarly, the evolution law of the dimensionless parameter r¯ for ε0=0.3 in Figure 3(b) shows that the pendulum length r¯ contracts into 1−ε0 in the illumination zone, while recovers gradually to 1 in the dark zone. The pendulum length also switches when the pendulum passes through the interface between the illumination and dark zones.Figure 3
Evolution laws of (a)T¯ and (b) r¯ of the PSM pendulum.
(a)(b)Based on the evolution law in Figure3, the mechanism of the self-excited oscillation in Figure 2(b) can be elucidated by investigating the time histories of T¯ and r¯ in detail. Figure 4(a) plots the dependence of the dimensionless temperature difference parameter T¯ on the swing angle θ. For −0.1≤θ≤0.1, the cycloid is in the illumination zone and T¯ increases. When the pendulum swings into the dark zone, T¯ decreases. Similarly, Figure 4(b) plots the dependence of r¯ on θ. The length r¯ decreases in the illumination zone, while it increases in the dark zone. Further, Figure 4(c) plots the dependence of F¯ on θ, where F¯ is discontinuous at some points of θ. The tension F¯ sharply increases at the moment of swinging into the dark region (point ②) from the illumination region (point ①). This is because the temperature of the pendulum swinging into the dark zone decreases and in turn the cycloid length increases. As the pendulum swings in the dark area, F¯ continuously varies from point ② to ③. Then, the tension F¯ sharply decreases from point ③ to point ④ when the pendulum reenters into the illumination zone. Eventually, tension F¯ and swing angle θ form a symmetrical closed loop according to the trajectory cycle of ①②③④⑤⑥⑦⑧① in a swing period.Figure 4
Mechanism of self-excited oscillation of the PSM pendulum, forθ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. (a) T¯ vs. θ. (b) r¯ vs. θ. (c) F¯ vs. θ. (d) F¯ vs. r¯.
(a)(b)(c)(d)Based on Figures4(a)–4(c), Figure 4(d) plots the dependence of F¯ with r¯, where the shadow area enclosed by curve of F¯ with r¯ represents the magnitude of the work done by the tension F¯ in a half swing period of the simple pendulum. In an oscillation period, the work done by F¯ is twice shadow area in Figure 4(d), which is the source of self-excited oscillation. When the net work is enough for compensating the dissipation by air damping, the pendulum can oscillate continuously.


## 4. Influence of Systematic Parameters on the Self-Excited Oscillation

In equations (15) and (16), there are six dimensionless parameters: φ, β¯, ε0, g¯, θ0, and θ˙0. The swing of the PSM pendulum in Figure 1(a) is controlled by the above six parameters. In this section, the influence of the six parameters on the onset, frequency, and amplitude of the self-excited oscillation is extensively investigated. Considering the parameter θ0 can be transformed into the corresponding value θ˙0 through the energy transformation between gravitational energy and kinetic energy, θ0 is set to be zero in the calculation. The dimensionless oscillation frequency f¯ of the simple pendulum is defined as the number of oscillations per unit time t¯, and the angular amplitude of the self-excited oscillation of the simple pendulum is denoted by A.


### 4.1. Effect of Initial Angular Velocityθ˙0

Figure5 shows the influence of the initial angular velocity θ˙0 on the swing of the simple pendulum, for θ0=0, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. Figure 5(a) plots the phase diagram of the swing angle for θ˙0=0.12. The result shows that the simple pendulum swings with gradually attenuating amplitude and eventually comes to rest, which is static mode. The maximum swing angle θmax of simple pendulum is less than φ in the initial stage of vibration, i.e., the pendulum is always in the illumination region. Therefore, the pendulum length r¯ decreases continuously and approaches the value of 1−ε0. Figure 5(b) plots the phase diagram of the swing angle for θ˙0=0.13. The result shows that the amplitude of the single pendulum first increases and then asymptotically approaches a fixed value, which is oscillation mode. As explained previously, the single pendulum absorbs light autonomously during alternately swinging between the illumination and dark zones, which results in self-excited oscillation of the simple pendulum. Figure 5(c) plots the limit cycles, in which there exists a critical θ˙0=0.12 for transition between static mode and oscillation mode.Figure 5
The effect ofθ˙0 on the swing of the simple pendulum, for θ0=0, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. (a) Phase diagram of the static mode for θ˙0=0.12. (b) Phase diagram of the oscillation mode for θ˙0=0.13. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)Figure5(d) plots the frequency and amplitude of self-excited oscillation for different θ˙0, respectively. For θ˙0>0.12, the frequency and amplitude are both zero because the pendulum is in static mode. For θ˙0<0.12, the frequency and amplitude of light-powered self-excited oscillation are the same and independent on θ˙0. These results mean that the initial state affects the onset of the oscillation mode, while does not affect the frequency and amplitude of the oscillation. This is consistent with the general characteristic of self-excited motions [38].


### 4.2. Effect of Illumination Zone

Figure6 shows the influence of the angle φ of the illumination zone on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, ε0=0.3, β¯=0.12, and g¯=2. Figure 6(a) plots the phase diagram of the swing angle of the pendulum, for φ=0.154. Initially, the maximum swing angle of the simple pendulum is larger than φ, and then, the simple pendulum swings with increasing amplitude alternately in the illuminated zone and the dark zone. The result shows that the oscillation mode of the pendulum is triggered. Figure 6(b) plots the phase diagram of the swing angle of the pendulum, for φ=0.155. The amplitude of the simple pendulum decreases continuously and the pendulum eventually comes to cease. This is because the pendulum with φ larger than the maximum swing angle is always in the illumination zone and gradually develops into static mode. Figure 6(c) plots the limit cycles, in which there also exists a critical φ for the phase transition between static mode and oscillation mode, which is numerically determined to be 0.154.Figure 6
The effect ofφ on the swing of the pendulum, for θ0=0, θ˙0=0.2, β¯=0.12, ε0=0.3, and g¯=2. (a) Phase diagram of the oscillation mode for φ=0.154. (b) Phase diagram of the static mode for φ=0.155. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)Figure6(d) plots the frequency and amplitude of self-excited oscillation for different φ values, respectively. For φ>0.154, the frequency and amplitude are both zero because the pendulum is in static mode. For φ<0.154, the frequency of light-powered self-excited oscillation is the same and independent on φ, while the amplitude is proportional to φ. These results mean that the illuminated zone affects the amplitude of the oscillation mode, while it does not affect the frequency of the oscillation. The result can be understood by the energy competition between energy input and damping dissipation. The increasing illuminated area amplifies the energy input and in turn enlarges the amplitude of the oscillation.


### 4.3. Effect of Damping Coefficientβ¯

Figure7 shows the influence of β¯ on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, ε0=0.3, and g¯=2. Figure 7(a) plots the phase diagram of the swing angle of the pendulum, for β¯=0.15. The result shows that the self-excited oscillation of the pendulum is triggered. For β¯=0.17, the pendulum develops into static mode, as shown in Figure 7(b). These results can be understood by the energy competition between energy input and damping dissipation. With the increase of β¯, the heat exchange between pendulum and environment becomes slower, and the net energy input during the swing decreases. Once the net energy input is not enough for compensating the damping dissipation, the pendulum will come to cease. Figure 7(c) plots the limit cycles, in which there also exists a critical β¯ for the transition between static mode and oscillation mode, which is numerically determined to be 0.17.Figure 7
The effect ofβ¯ on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, ε0=0.3, and g¯=2. (a) Phase diagram of the oscillation mode for β¯=0.15. (b) Phase diagram of the static mode for β¯=0.17. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)For smallβ¯, the amplitude calculated may be larger than π/2, i.e., the simple pendulum may rotate or the tension of the cycloid is negative, which is beyond the scope of this paper. Here, β¯ is chosen as β¯≤0.02. Figure 7(d) plots the frequency and amplitude of self-excited oscillation for different β¯, respectively. For β¯>0.17, the frequency and amplitude are both zero because the pendulum is in static mode. For β¯<0.17, the frequency of light-powered self-excited oscillation increases slightly, while the amplitude decreases obviously and gradually approaches the fixed value, with the increase of β¯. The result can also be understood by the energy competition between energy input and damping dissipation. The increasing β¯ decreases the energy input and in turn attenuates the amplitude of the oscillation. The result implies that, in order to improve the absorption of light energy by the pendulum, β¯ can be appropriately reduced for large energy harvesting in engineering.


### 4.4. Effect of Contraction Coefficientε0

Figure8 shows the influence of ε0 on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. Figure 8(a) plots the phase diagram of the swing angle of the pendulum, for ε0=0.22. The pendulum is in static mode because the small strain of the cycloid cannot provide enough net energy input to compensate the energy dissipation by air damping. Figure 8(b) plots the phase diagram of the swing angle of the pendulum, for ε0=0.25. The result shows that the oscillation mode of the pendulum is triggered. In this case, the light-driven strain can provide enough net energy input for the oscillation. Figure 8(c) plots the limit cycles, in which there also exists a critical ε0 for the phase transition between static mode and oscillation mode, which is numerically determined to be 0.25.Figure 8
The effect ofε0 on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. (a) Phase diagram of the static mode at ε0=0.22. (b) Phase diagram of the oscillation mode at ε0=0.25. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)For largerε0, the pendulum can also oscillate continuously. However, according to the physical meaning of contraction strain ε0, the maximum value should be smaller than 1. Here, ε0 is chosen as 0.25≤ε0≤0.9. Figure 8(d) plots the frequency and amplitude of self-excited oscillation for different ε0, respectively. For ε0<0.25, the frequency and amplitude are both zero because the pendulum is in static mode. For ε0>0.25, the increase of ε0 has no obvious effect on the self-excited oscillation frequency, while the amplitude increases dramatically with the increase of ε0. Therefore, in order to improve the engineering application of photothermal-responsive pendulum, light-driven contractile or expansive strain of the material should be amplified for the efficient conversion of light energy to mechanical energy.


### 4.5. Effect of Gravitational Accelerationg¯

Figure9 shows the influence of g¯ on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and ε0=0.3. Figure 9(a) plots the phase diagram of the swing angle of the pendulum for g¯=0.4, which is static mode. Figures 9(b) and 9(c) plot the phase diagrams of the swing angle of the pendulum, for g¯=0.5 and g¯=4.7, respectively. The result shows that the oscillation mode of the pendulum is triggered. Careful calculations show that, for 0.4≤g¯≤4.7, the simple pendulum is in oscillation mode. However, for g¯=4.8, the pendulum swings in static mode again, as shown in Figure 9(d). Figure 9(e) also plots the limit cycles of the light-powered swing of the pendulum. The results can be understood by the energy compensation between net energy input and damping dissipation. For small g¯, i.e., the heat exchange is too fast, the cycloid rapidly deforms just only when passing through the interface between the illumination and dark zones. For large g¯, i.e., the heat exchange is too slow, the cycloid almost does not deform during the swing. Therefore, the net work done of single pendulum by the light is too small and cannot compensate the energy dissipated by air damping to maintain the oscillation.Figure 9
The effect ofg¯ on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and ε0=0.3. (a) Phase diagram of the static mode at g¯=0.4. (b) Phase diagram of the oscillation mode at g¯=0.5. (c) Phase diagram of the oscillation mode at g¯=4.7. (d) Phase diagram of the static mode at g¯=4.8. (e) Limit cycles. (f) Amplitude and frequency.
(a)(b)(c)(d)(e)(f)Figure9(f) plots the frequency and amplitude of self-excited oscillation of the simple pendulum for 0.4≤g¯≤4.7, respectively. It can be seen that the oscillation frequency and amplitude of the pendulum can be significantly increased by increasing g¯. According to the physical meaning of g¯, reducing the heat transfer rate within a certain range can amplify the oscillation frequency and amplitude and improve the efficiency of light energy to mechanical energy in engineering application.


## 4.1. Effect of Initial Angular Velocityθ˙0

Figure5 shows the influence of the initial angular velocity θ˙0 on the swing of the simple pendulum, for θ0=0, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. Figure 5(a) plots the phase diagram of the swing angle for θ˙0=0.12. The result shows that the simple pendulum swings with gradually attenuating amplitude and eventually comes to rest, which is static mode. The maximum swing angle θmax of simple pendulum is less than φ in the initial stage of vibration, i.e., the pendulum is always in the illumination region. Therefore, the pendulum length r¯ decreases continuously and approaches the value of 1−ε0. Figure 5(b) plots the phase diagram of the swing angle for θ˙0=0.13. The result shows that the amplitude of the single pendulum first increases and then asymptotically approaches a fixed value, which is oscillation mode. As explained previously, the single pendulum absorbs light autonomously during alternately swinging between the illumination and dark zones, which results in self-excited oscillation of the simple pendulum. Figure 5(c) plots the limit cycles, in which there exists a critical θ˙0=0.12 for transition between static mode and oscillation mode.Figure 5
The effect ofθ˙0 on the swing of the simple pendulum, for θ0=0, φ=0.1, β¯=0.12, ε0=0.3, and g¯=2. (a) Phase diagram of the static mode for θ˙0=0.12. (b) Phase diagram of the oscillation mode for θ˙0=0.13. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)Figure5(d) plots the frequency and amplitude of self-excited oscillation for different θ˙0, respectively. For θ˙0>0.12, the frequency and amplitude are both zero because the pendulum is in static mode. For θ˙0<0.12, the frequency and amplitude of light-powered self-excited oscillation are the same and independent on θ˙0. These results mean that the initial state affects the onset of the oscillation mode, while does not affect the frequency and amplitude of the oscillation. This is consistent with the general characteristic of self-excited motions [38].


## 4.2. Effect of Illumination Zone

Figure6 shows the influence of the angle φ of the illumination zone on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, ε0=0.3, β¯=0.12, and g¯=2. Figure 6(a) plots the phase diagram of the swing angle of the pendulum, for φ=0.154. Initially, the maximum swing angle of the simple pendulum is larger than φ, and then, the simple pendulum swings with increasing amplitude alternately in the illuminated zone and the dark zone. The result shows that the oscillation mode of the pendulum is triggered. Figure 6(b) plots the phase diagram of the swing angle of the pendulum, for φ=0.155. The amplitude of the simple pendulum decreases continuously and the pendulum eventually comes to cease. This is because the pendulum with φ larger than the maximum swing angle is always in the illumination zone and gradually develops into static mode. Figure 6(c) plots the limit cycles, in which there also exists a critical φ for the phase transition between static mode and oscillation mode, which is numerically determined to be 0.154.Figure 6
The effect ofφ on the swing of the pendulum, for θ0=0, θ˙0=0.2, β¯=0.12, ε0=0.3, and g¯=2. (a) Phase diagram of the oscillation mode for φ=0.154. (b) Phase diagram of the static mode for φ=0.155. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)Figure6(d) plots the frequency and amplitude of self-excited oscillation for different φ values, respectively. For φ>0.154, the frequency and amplitude are both zero because the pendulum is in static mode. For φ<0.154, the frequency of light-powered self-excited oscillation is the same and independent on φ, while the amplitude is proportional to φ. These results mean that the illuminated zone affects the amplitude of the oscillation mode, while it does not affect the frequency of the oscillation. The result can be understood by the energy competition between energy input and damping dissipation. The increasing illuminated area amplifies the energy input and in turn enlarges the amplitude of the oscillation.


## 4.3. Effect of Damping Coefficientβ¯

Figure7 shows the influence of β¯ on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, ε0=0.3, and g¯=2. Figure 7(a) plots the phase diagram of the swing angle of the pendulum, for β¯=0.15. The result shows that the self-excited oscillation of the pendulum is triggered. For β¯=0.17, the pendulum develops into static mode, as shown in Figure 7(b). These results can be understood by the energy competition between energy input and damping dissipation. With the increase of β¯, the heat exchange between pendulum and environment becomes slower, and the net energy input during the swing decreases. Once the net energy input is not enough for compensating the damping dissipation, the pendulum will come to cease. Figure 7(c) plots the limit cycles, in which there also exists a critical β¯ for the transition between static mode and oscillation mode, which is numerically determined to be 0.17.Figure 7
The effect ofβ¯ on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, ε0=0.3, and g¯=2. (a) Phase diagram of the oscillation mode for β¯=0.15. (b) Phase diagram of the static mode for β¯=0.17. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)For smallβ¯, the amplitude calculated may be larger than π/2, i.e., the simple pendulum may rotate or the tension of the cycloid is negative, which is beyond the scope of this paper. Here, β¯ is chosen as β¯≤0.02. Figure 7(d) plots the frequency and amplitude of self-excited oscillation for different β¯, respectively. For β¯>0.17, the frequency and amplitude are both zero because the pendulum is in static mode. For β¯<0.17, the frequency of light-powered self-excited oscillation increases slightly, while the amplitude decreases obviously and gradually approaches the fixed value, with the increase of β¯. The result can also be understood by the energy competition between energy input and damping dissipation. The increasing β¯ decreases the energy input and in turn attenuates the amplitude of the oscillation. The result implies that, in order to improve the absorption of light energy by the pendulum, β¯ can be appropriately reduced for large energy harvesting in engineering.


## 4.4. Effect of Contraction Coefficientε0

Figure8 shows the influence of ε0 on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. Figure 8(a) plots the phase diagram of the swing angle of the pendulum, for ε0=0.22. The pendulum is in static mode because the small strain of the cycloid cannot provide enough net energy input to compensate the energy dissipation by air damping. Figure 8(b) plots the phase diagram of the swing angle of the pendulum, for ε0=0.25. The result shows that the oscillation mode of the pendulum is triggered. In this case, the light-driven strain can provide enough net energy input for the oscillation. Figure 8(c) plots the limit cycles, in which there also exists a critical ε0 for the phase transition between static mode and oscillation mode, which is numerically determined to be 0.25.Figure 8
The effect ofε0 on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and g¯=2. (a) Phase diagram of the static mode at ε0=0.22. (b) Phase diagram of the oscillation mode at ε0=0.25. (c) Limit cycles. (d) Amplitude and frequency.
(a)(b)(c)(d)For largerε0, the pendulum can also oscillate continuously. However, according to the physical meaning of contraction strain ε0, the maximum value should be smaller than 1. Here, ε0 is chosen as 0.25≤ε0≤0.9. Figure 8(d) plots the frequency and amplitude of self-excited oscillation for different ε0, respectively. For ε0<0.25, the frequency and amplitude are both zero because the pendulum is in static mode. For ε0>0.25, the increase of ε0 has no obvious effect on the self-excited oscillation frequency, while the amplitude increases dramatically with the increase of ε0. Therefore, in order to improve the engineering application of photothermal-responsive pendulum, light-driven contractile or expansive strain of the material should be amplified for the efficient conversion of light energy to mechanical energy.


## 4.5. Effect of Gravitational Accelerationg¯

Figure9 shows the influence of g¯ on the swing of the simple pendulum, for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and ε0=0.3. Figure 9(a) plots the phase diagram of the swing angle of the pendulum for g¯=0.4, which is static mode. Figures 9(b) and 9(c) plot the phase diagrams of the swing angle of the pendulum, for g¯=0.5 and g¯=4.7, respectively. The result shows that the oscillation mode of the pendulum is triggered. Careful calculations show that, for 0.4≤g¯≤4.7, the simple pendulum is in oscillation mode. However, for g¯=4.8, the pendulum swings in static mode again, as shown in Figure 9(d). Figure 9(e) also plots the limit cycles of the light-powered swing of the pendulum. The results can be understood by the energy compensation between net energy input and damping dissipation. For small g¯, i.e., the heat exchange is too fast, the cycloid rapidly deforms just only when passing through the interface between the illumination and dark zones. For large g¯, i.e., the heat exchange is too slow, the cycloid almost does not deform during the swing. Therefore, the net work done of single pendulum by the light is too small and cannot compensate the energy dissipated by air damping to maintain the oscillation.Figure 9
The effect ofg¯ on the swing of the pendulum for θ0=0, θ˙0=0.2, φ=0.1, β¯=0.12, and ε0=0.3. (a) Phase diagram of the static mode at g¯=0.4. (b) Phase diagram of the oscillation mode at g¯=0.5. (c) Phase diagram of the oscillation mode at g¯=4.7. (d) Phase diagram of the static mode at g¯=4.8. (e) Limit cycles. (f) Amplitude and frequency.
(a)(b)(c)(d)(e)(f)Figure9(f) plots the frequency and amplitude of self-excited oscillation of the simple pendulum for 0.4≤g¯≤4.7, respectively. It can be seen that the oscillation frequency and amplitude of the pendulum can be significantly increased by increasing g¯. According to the physical meaning of g¯, reducing the heat transfer rate within a certain range can amplify the oscillation frequency and amplitude and improve the efficiency of light energy to mechanical energy in engineering application.


## 5. Conclusions

Photothermal-responsive materials can convert light energy into mechanical energy for driving active machines. In this work, a simple pendulum capable of self-excited oscillation is innovatively constructed by using ordinary PSM or PEM as the cycloid, and the dynamic model is established for investigating the light-powered swing of the photothermal-responsive pendulum. By theoretical formulation and numerical calculation, the self-excited oscillation behaviors of the simple pendulum are extensively studied. When the PSM pendulum is irradiated in the central zone of the swing, or conversely the PEM pendulum is illuminated on the two sides of the swing, the simple pendulum could oscillate autonomously. There exist two motion modes of the pendulum, which are static mode and oscillation mode. The numerical calculation provides time histories of the cycloid length and cycloid tension in detail and in turn elucidates the mechanism of self-excited oscillation that results from the coupling of heat exchange, photothermal deformation, and swing. When the net input energy compensates the damping dissipation, the pendulum keeps swinging. This paper further investigates the influence of physical parameters on the self-excited oscillation. The results show that the amplitude of the self-excited oscillation can be effectively increased by increasing the illumination zone, reducing the air damping dissipation, increasing the limit contraction strain of the material, or decreasing the heat transfer rate of the photothermal material, while the frequency of the self-excited oscillation is almost only affected by the heat transfer coefficient of the photothermal-responsive material. Unfortunately, as far as we know, there are no experimental studies in the literature to verify our findings. Next, we do plan to carry out corresponding experiments to further confirm our results. We hope the light-powered self-excited oscillation of photothermal-responsive pendulums can not only deepen the understanding of self-excited oscillations but also have potential applications in harvesting energy, signal sensing, soft robots, and so on.



---

*Source: 5531530-2021-02-16.xml*