Next Article in Journal
Synapse-Neuron-Aware Training Scheme of Defect-Tolerant Neural Networks with Defective Memristor Crossbars
Previous Article in Journal
An Experiment-Based Study of Formation Damage Using a Microetching Model Displacement Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Self-Sustained Collective Motion of Two Joint Liquid Crystal Elastomer Spring Oscillator Powered by Steady Illumination

Department of Civil Engineering, Anhui Jianzhu University, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(2), 271; https://doi.org/10.3390/mi13020271
Submission received: 26 December 2021 / Revised: 25 January 2022 / Accepted: 31 January 2022 / Published: 8 February 2022

Abstract

:
For complex micro-active machines or micro-robotics, it is crucial to clarify the coupling and collective motion of their multiple self-oscillators. In this article, we construct two joint liquid crystal elastomer (LCE) spring oscillators connected by a spring and theoretically investigate their collective motion based on a well-established dynamic LCE model. The numerical calculations show that the coupled system has three steady synchronization modes: in-phase mode, anti-phase mode, and non-phase-locked mode, and the in-phase mode is more easily achieved than the anti-phase mode and the non-phase-locked mode. Meanwhile, the self-excited oscillation mechanism is elucidated by the competition between network that is achieved by the driving force and the damping dissipation. Furthermore, the phase diagram of three steady synchronization modes under different coupling stiffness and different initial states is given. The effects of several key physical quantities on the amplitude and frequency of the three synchronization modes are studied in detail, and the equivalent systems of in-phase mode and anti-phase mode are proposed. The study of the coupled LCE spring oscillators will deepen people’s understanding of collective motion and has potential applications in the fields of micro-active machines and micro-robots with multiple coupled self-oscillators.

1. Introduction

Self-excited oscillation is a kind of periodic motion that is maintained by constant external excitations [1,2,3,4]. Similar to biological active feeding, it can directly harvest energy from a constant environment to maintain its periodic motion [5,6]. In addition, the period and amplitude of the self-oscillation generally depend on the intrinsic parameters of the system and are independent of the initial conditions, which makes the system robust [7,8,9]. Due to the unique advantages of self-excited oscillation systems, they have broad application prospects in the fields of energy acquisition, sensing with electronic skins [10], soft robotics [11,12,13], medical instruments [14,15,16,17], and motors [18]. The self-excited oscillations are generally based on responsive materials, including liquid crystal elastomers (LCEs) [19,20], dielectric elastomers [21], hydrogels [22,23,24], and ion gels [25,26]. Based on different stimuli-responsive materials and structures, different feedback mechanisms are proposed to realize energy compensation, such as the coupling of chemical reactions and large deformation [25,26], the self-shadowing effect [14,27], the coupling of liquid volatilization, and membrane deformation [28].
Based on the self-oscillation systems that have been previously reported [29], the coupling and synchronization phenomena of two or more self-excited oscillation systems and their collective motion have attracted extensive attention [30,31,32,33,34]. Synchronization and collective motion are ubiquitous in nature, such as the overall movement of a school of fish or a flock of wild geese or a group of fireflies flickering together [35,36,37]. As early as 1673, C. Huygens studied the synchronization phenomenon of simple pendulum coupling by observing two identical clocks oscillating synchronously with two pendulums swinging in opposite directions [38]. Recent research has confirmed that the coupling between two pendulums is caused by tiny mechanical oscillations that propagate through the wooden structure in which the clocks are mounted. Similar experiments have enabled a large number of metronomes to swing synchronously on a freely moving base [39,40]. Recently, based on optically responsive LCE, Ghislaine et al. experimentally studied the synchronized oscillations of thin plastic actuators fueled by light and found two kinds of in-phase and anti-phase synchronous oscillation phenomena in the steady-state [30,31]. Their numerical simulations qualitatively explained the origin of synchronized motion and found that motion can be regulated by the mechanical properties of coupling.
LCEs are advanced multifunctional materials that combine the flexibility of polymeric networks with the nematic structure of liquid crystals [41,42], which have the advantages of having a fast response, recoverable deformation, and low noise [10,43,44,45]. This special composition and structure enable LCEs to respond to external light [46,47,48,49], heat [50], electric fields [51,52], magnetic fields [53] and chemical substances [54]. Based on LCE materials, several self-exciting motion modes have been constructed, such as rolling [20], vibration [17], swinging [10,55], stretching and shrinking [56], rotation [57], eversion or inversion [9,58], torsion [59], jumping [60], and buckling [61] modes. These self-exciting motion modes provide good ideas for studying the coupling of multiple systems and their collective motion.
Studying the self-excited oscillation coupling of two or more systems and their collective motion is beneficial to the construction of richer and more complex types of motion, allowing them to demonstrate more versatile functions in the micro-robots. In view of this, we used two identical LCE fibers connected by a spring to construct a new optically responsive LCE spring oscillation coupling system. The self-oscillation mechanism and the possible in-phase mode, anti-phase mode, and non-phase-locked mode are discussed, and the influence of some key physical quantities on its synchronization modes, amplitude, and period are analyzed. The layout of this paper is as follows. In Section 2, based on the LCE dynamic model first proposed by Finkelmann et al. [62,63], the governing equation is derived, and the difference schemes of the dynamics equation and the solution method are given. In Section 3, the different synchronization modes of the self-excited oscillation are discussed, and the detailed mechanisms are revealed. In Section 4, the effects of various system parameters on the self-excited oscillation in the three synchronization modes are studied in detail. In Section 5, the equivalent systems in the in-phase mode and anti-phase mode are provided. Finally, the concluding remarks are shared in Section 6.

2. Model and Formulation

2.1. Dynamic Model of the Two LCE Spring Oscillators

Figure 1 sketches a dynamic model of the self-oscillation coupling system under uniform and constant illumination, which is composed of two identical LCE spring oscillators and a spring. The illuminated zone is represented by the shaded area. The original length of the LCE fiber in a stress-free state is L f , and the original length of the spring is L s , as shown in Figure 1a. Next, fix one end of the LCE fiber and tie the mass on the other end, and connect the two fibers with a spring. For simplicity, the gravity of the mass blocks is ignored, and an initial strain λ p 1 is assigned to the LCE fibers to ensure that the system hang horizontally, as shown in Figure 1b. The displacement of the two mass blocks is denoted by u 1 t and u 2 t , respectively. The current lengths of LCE fiber 1 and LCE fiber 2 are l f 1 t = L f λ p + u 1 t and l f 2 t = L f λ p u 2 t , respectively, as shown in Figure 1c. F f 1 t and F f 2 t are the spring forces of LCE fiber 1 and LCE fiber 2, respectively, which are also called driving forces hereinafter. F s t is the spring force of the spring (abbreviated as spring force), and F d 1 t and F d 2 t are the damping forces during oscillation. For simplicity, it is assumed that the damping forces are proportional to the velocity of the mass, and the direction is always opposite to the velocity of the mass.
In order to analyze the inhomogeneous deformation of the two LCE fibers, the Lagrangian coordinate systems X 1 and X 2 are fixed and established for the initial configuration of LCE fiber 1 and LCE fiber 2, and the Eulerian coordinate systems x 1 and x 2 in the current configuration are also established. The instantaneous position of a material point X 1 ( X 2 ) of LCE fiber 1 (LCE fiber 2) can be represented as x 1 = x 1 X 1 , t ( x 2 = x 2 X 2 , t ) during the oscillation. Since the pulling forces of the LCE fibers are much greater than the gravity of the mass blocks and the LCE fibers, and we ignore the gravity of the mass blocks and the LCE fibers for simplicity. According to Newtonian mechanics, the following governing equation holds at any moment during mass oscillation.
m u ¨ 1 = F f 1 t + F s t c u ˙ 1 m u ¨ 2 = F f 2 t F s t c u ˙ 2 ,
where c is the damping coefficient, u ˙ and u ¨ indicate the velocity d u t d t and acceleration d 2 u ( t ) d t 2 of the mass, respectively.
For simplicity, the force of the LCE fiber F f t is assumed to be proportional to the elastic strain ε e t , i.e.,
F f 1 t = K f L f ε e 1 t F f 2 t = K f L f ε e 2 t ,  
where K f is the spring constant of the LCE fiber, and ε e t = ε tot t ε t with ε tot t are the total strain, and ε tot t is the light-driven contraction. For simplicity, the total strain is defined as ε tot t = λ X , t 1 [64], where λ ( X , t ) are written as
λ 1 ( X 1 , t ) = d x 1 X 1 , t d X 1 ,   λ 2 ( X 2 , t ) = d x 2 X 2 , t d X 2 .  
Then, Equation (2) can be rewritten as
F f 1 t = K f L f λ 1 X 1 , t 1 ε 1 X 1 , t F f 2 t = K f L f λ 2 X 2 , t 1 ε 2 X 2 , t ,
where the light-driven contraction strain ε ( X , t ) is assumed to be proportional to the volume fraction of the isomers in the cis state φ ( X , t ) in the LCE fiber, which can be written as
ε 1 X 1 , t = C 0 φ 1 X 1 , t ε 2 X 2 , t = C 0 φ 2 X 2 , t ,
where C 0 is the contraction coefficient.
In order to find the instantaneous position x 1 ( x 2 ) of the LCE fiber 1 (LCE fiber 2) at any time, we first rewrite F f 1 ( t ) ( F f 2 ( t ) ) as u 1 t ( u 2 t ) and ε 1 ( X 1 , t ) ( ε 2 ( X 1 , t ) ). Considering that the LCE fiber is in uniaxial tensile state, the axial force is homogeneous although the contraction is inhomogeneous. Noting that F f 1 ( t ) and F f 2 ( t ) are homogeneous and constant in the LCE fibers, by integrating Equation (4) from 0 to L f on both sides, we can obtain
F f 1 t = K f L f λ p 1 + u 1 ( t ) 0 L f ε 1 X 1 , t d X 1 F f 2 t = K f L f λ p 1 u 2 ( t ) 0 L f ε 2 X 2 , t d X 2 ,
where λ p is the pre-stretch of the LCE fiber.
Then, from Equation (2), λ 1 X 1 , t , λ 2 X 2 , t can be expressed by F f 1 ( t ) , F f 2 ( t ) as
λ 1 X 1 , t = F f 1 t K f L f + 1 + ε 1 X 1 , t λ 2 X 2 , t = F f 2 t K f L f + 1 + ε 2 X 2 , t .
By combining Equations (3) and (7), we can obtain
d x 1 X 1 , t = F f 1 t K f L f + 1 + ε 1 X 1 , t d X 1 d x 2 X 2 , t = F f 1 t K f L f + 1 + ε 1 X 1 , t d X 2 .  
Combining Equations (6) and (8), we achieve
d x 1 X 1 , t = ε 1 X 1 , t + λ p + u 1 ( t ) 0 L f ε 1 X 1 , t d X 1 L f d X 1 d x 2 X 2 , t = ε 1 X 1 , t + λ p u 2 ( t ) + 0 L f ε 2 X 2 , t d X 2 L f d X 2 .  
By integrating both sides of the above formula from 0 to X , we achieve
x 1 X 1 , t = λ p X 1 + 0 X 1 ε 1 X 1 , t d X 1 + X 1 L f u 1 t 0 L f ε 1 X 1 , t d X 1 x 2 X 2 , t = λ p X 2 + 0 X 2 ε 2 X 2 , t d X 2 X 2 L f u 2 t + 0 L f ε 2 X 2 , t d X 2 .
The calculated x 1 and x 2 from Equation (10) can be compared to λ p L s to determine whether the LCE fiber matter point is in the illuminated or non-illuminated area.
The two LCE fibers are coupled by the spring in the middle. Therefore, it is necessary for us to express the spring force F s t . During the oscillation process, the elongation Δ L s ( t ) of the spring at any moment can be expressed by the displacement of the two mass blocks u 1 ( t ) and u 2 ( t ) as
Δ L s ( t ) = u 2 ( t ) u 1 ( t ) + K f L f λ p 1 K s .  
According to the deformation Δ L s ( t ) of the spring and spring constant (coupling stiffness) K s of the spring, the spring force can be expressed as
F s ( t ) = K s u 2 ( t ) u 1 ( t ) + K f L f λ p 1 .  

2.2. Evolution Law of Number Fraction in the Two LCE Fibers

In order to calculate the light-driven contraction strain, we need to obtain the number fractions in the LCE fibers. According to the research of Yu et al., the trans-to-cis isomerization of LCE can be induced by UV or laser with a wavelength of less than 400 nm [62]. Generally, under UV light excitation, light-driven cis-to-trans isomerization can be neglected [65], and the number fraction of cis-isomers ϕ ( t ) depends on thermal excitation from trans to cis, thermally driven relaxation from cis to trans, and light-driven trans-to-cis isomerization. Considering that the thermal excitation from trans to cis is often negligible compared to the light-driven excitation, we define the number fraction of cis isomers in LCE fibers and use the following governing equations to describe the evolution of the number fraction of the cis isomers [66].
ϕ 1 X 1 , t t = η 0 I 0 1 ϕ 1 X 1 , t T 0 1 ϕ 1 X 1 , t ϕ 2 X 2 , t t = η 0 I 0 1 ϕ 2 X 2 , t T 0 1 ϕ 2 X 2 , t ,
where T 0 is the thermal relaxation time responding to the cis state to trans state, I 0 is the light intensity, and η 0 is a light-absorption constant.

2.3. Nondimensionalization

By defining the following dimensionless parameters: F ˜ f 1 t = F f 1 t T 0 2 / m L f , F ˜ f 2 t = F f 2 t T 0 2 / m L f , c ˜ = c T 0 / m , F ˜ s ( t ) = F s t T 0 2 / m L f , u ˜ 1 ( t ) = u 1 ( t ) / L f , u ˜ 2 ( t ) = u 2 ( t ) / L f , t ˜ = t / T 0 , K ˜ f = K f T 0 2 / m , and K ˜ s = K s T 0 2 / m , Equation (10) can be rewritten as
x ˜ 1 X ˜ 1 , t ˜ = λ p + u ˜ 1 t ˜ X ˜ 1 X ˜ 1 0 1 ε 1 X ˜ 1 , t ˜ d X ˜ 1 + 0 X ˜ 1 ε 1 X ˜ 1 , t ˜ d X ˜ 1 x ˜ 2 X ˜ 2 , t ˜ = λ p u ˜ 2 t ˜ X ˜ 2 X ˜ 2 0 1 ε 2 X ˜ 2 , t ˜ d X ˜ 2 + 0 X ˜ 2 ε 2 X ˜ 2 , t ˜ d X ˜ 2 .  
Combining Equations (1), (6), and (12), we achieve
u ¨ ˜ 1 = K ˜ s u ˜ 2 ( t ) u ˜ 1 ( t ) K ˜ f u ˜ 1 ( t ˜ ) 0 1 ε 1 X ˜ 1 , t ˜ d X ˜ 1 c u ˙ ˜ 1 u ¨ ˜ 2 = K ˜ f u ˜ 2 ( t ˜ ) + 0 1 ε 1 X ˜ 1 , t ˜ d X ˜ 1 K ˜ s u ˜ 2 ( t ˜ ) u ˜ 1 ( t ) c u ˙ ˜ 2 ,
where u ˙ ˜ and u ¨ ˜ indicate the velocity d u ˜ t ˜ d t ˜ and acceleration d 2 u ˜ t ˜ d t 2 of the mass, respectively.
By defining the parameter I ˜ = T 0 η 0 I 0 , Equation (13) can be rewritten as
ϕ 1 X ˜ 1 , t ˜ t ˜ = I ˜ 1 + I ˜ ϕ 1 X ˜ 1 , t ˜ ϕ 2 X ˜ 2 , t ˜ t ˜ = I ˜ 1 + I ˜ ϕ 2 X ˜ 2 , t ˜ .  

2.4. Solution Method

In this paper, the LCE fiber is discretized into J material points (let J = 500 for both LCE fibers in this paper) in the numerical calculation process. The position vector of each material point of LCE fiber 1 and LCE fiber 2 in the Lagrangian coordinate system can be expressed as X 1 = X 1 1 , X 1 2 X 1 J , X 2 = X 2 1 , X 2 2 X 2 J , and is called the material coordinate. The position vector of each material point of LCE fiber 1 and LCE fiber 2 in the Eulerian coordinate system can be expressed as x 1 = x 1 1 , x 1 2 x 1 J , x 2 = x 2 1 , x 2 2 x 2 J , which is called the spatial coordinate. By using the different methods to solve the dimensionless equation in Equation (16), the expression for the cis number fraction at any moment in the LCE fiber can be written as
ϕ 1 n + 1 = ϕ 1 n + I ˜ ( 1 + I ˜ ) ϕ 1 n Δ t ˜ ϕ 2 n + 1 = ϕ 2 n + I ˜ ( 1 + I ˜ ) ϕ 2 n Δ t ˜ .  
Substituting Equation (17) into Equation (4), we are able to obtain the light-driven contraction strain ε 11 , ε 21 for LCE fiber 1 and LCE fiber 2 located in the illuminated area
ε 11 = C 0 ϕ 1 n + I ˜ 1 + I ˜ ϕ 1 n Δ t ˜ ε 21 = C 0 ϕ 2 n + I ˜ 1 + I ˜ ϕ 2 n Δ t ˜ .  
Similarly, when the LCE fibers are located in a non-illuminated area, the shrinkage strain ε 12 , ε 22 can be written as
ε 12 = C 0 ϕ 1 n ϕ 1 n Δ t ˜ ε 22 = C 0 ϕ 2 n ϕ 2 n Δ t ˜ .  
Equation (15) is an ordinary differential equation with variable coefficients of second order, so no analytical solution can be obtained. Herein, we used the classical fourth-order Runge–Kutta method to solve the differential equation using Matlab software. By iterating Equations (15) and (17)–(19) we are able to obtain the final steady-state response of the LCE spring oscillator, i.e., the relationship between displacement and velocity with time histories.

3. Three Synchronization Modes and Their Mechanisms

3.1. Three Synchronization Modes

To investigate the collective motion of the two spring oscillators, we first need to estimate the typical values of the dimensionless parameters in the model. From the accessible experiments [19,49,66,67], the typical values of the material properties and geometric parameters are listed in Table 1.
Figure 2 shows three steady synchronization modes of the self-oscillation coupling system. In the computation, we fix K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0.2 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . For K ˜ s = 0.05 ; the result shows that the two LCE oscillators oscillate the anti-phase, and the time–history curve and domain of attraction of the anti-phase mode are given in Figure 2b,c, respectively. For K ˜ s = 0.15 , the results show that the amplitude of the two LCE oscillators changes periodically, which means that the self-oscillation coupling system is in non-phase-locked mode, as shown in Figure 2d,e. Figure 2f is the domain of attraction in the non-phase-locked mode. For K ˜ s = 0.5 , the calculation shows that the two curves of the two LCE oscillators coincide, which means that the two LCE oscillators are in in-phase mode, as show in Figure 2g,h. Figure 2i is the domain of attraction of u ˜ 1 and u ˜ 2 in in-phase mode. When other physical parameters remain unchanged, the conversion of different synchronization modes can be realized by changing one of the parameters, such as by changing the coupling stiffness K ˜ s .

3.2. The Mechanism of Self-Excited Oscillation

To investigate the mechanism of the self-excited oscillation of the LCE oscillator under uniform and constant illumination, Figure 3 plots the mechanism of the self-excited oscillation in anti-phase mode for K ˜ f = 5 , K ˜ s = 0.1 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.4 , and u ˙ ˜ 2 0 = 0.4 . During self-excited oscillation, the number fraction ϕ 1 ϕ 2 of cis-isomers, the shrinkage strain ε 1 ε 2 , the driving force F ˜ f 1 F ˜ f 2 , the spring force F ˜ s , and the displacement u ˜ 1 u ˜ 2 of mass block change periodically with the different time histories. In Figure 3d, the driving force and the displacement of the mass show a closed-loop relationship, and the area enclosed by the closed loop represents the network created by the driving force, which compensates the energy loss of the system and maintains the periodic self-excited oscillation of the system. Figure 3f represents the dependence of the displacement on time. It can be seen from this figure that the two mass blocks have a phase difference of half a cycle.
Numerical calculations show that the mechanisms of the self-excited oscillation of in-phase mode and non-phase-locked mode are similar to that of anti-phase mode. The number fraction of cis-isomers, the shrinkage strain, the driving force, the spring force, and the displacement of two mass blocks change periodically with time histories too. Driving force and mass displacement are also in a closed-loop relationship, and the area enclosed by the closed loop represents the network created by the driving force. This is because a part of the LCE fiber continuously enters and exits the illuminated area, causing the periodic contraction and relaxation of the LCE fiber, which causes the self-excited oscillation of the two mass blocks.

3.3. Triggering Conditions for Three Synchronization Modes

Numerical calculations show that the synchronization mode depends on the combination of various parameters. We took the combination of coupling stiffness and initial state as an example to study the effects of the system parameters on synchronous mode. In the computation, we fixed K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 and studied the effects of the phase difference between the two mass blocks and the coupling stiffness on the three synchronization modes at the initial moment by changing u ˙ ˜ 2 0 and K ˜ s . By analyzing a large amount of calculation data, the attractive domains of the three steady synchronization modes can be drawn, as shown in Figure 4.
Figure 4 shows the phase diagram of three steady synchronization modes under different coupling stiffnesses and different initial states. It can be seen from the figure that the smaller the coupling stiffness and the greater the phase difference between the two mass blocks of the initial state are, the easier it is to obtain the anti-phase mode. Conversely, it is conducive to the realization of the in-phase mode. In general, the in-phase mode is more easily achieved than the anti-phase mode, indicating a difference in robustness and the corresponding basins of attraction between the three synchronization modes. When the phase difference between the two LCE spring oscillators corresponding to the initial condition is large enough and the coupling stiffness K ˜ s is of moderate size, then there is a stable non-phase-locked mode between the in-phase mode and anti-phase mode. It is worth mentioning that when the two mass blocks are in anti-phase at the initial state (i.e., u ˙ ˜ 1 0 = 0.3 ), then the phase diagram of the anti-phase synchronization mode is the largest. The system is in anti-phase mode when K ˜ s < 2.9 , and when K ˜ s 2.9 , it is in in-phase mode.
After sorting out a large amount of calculation data, we briefly summarize the influence of other parameters on the three steady synchronization modes. When other parameters remain unchanged and the smaller the light intensity and contraction coefficient are, the easier it is to obtain the anti-phase mode, and the larger the damping coefficient and the spring coefficient of the LCE fiber are, the easier it is to obtain the anti-phase mode. Conversely, it is conducive to the realization of the in-phase mode. When the size of other parameters is moderate, there is a stable non-phase-locked mode between the in-phase mode and anti-phase mode. Studying the phase diagram of three steady synchronization modes under different parameters is extremely important in the design of automatic machines and soft robots.

4. Parametric Study

4.1. Effect of Coupling Stiffness

Figure 5a,b plot the effect of three different coupling stiffnesses on the self-oscillation model in the in-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.2 . Under the action of different coupling stiffness, the domain of attraction u ˜ 1 t and u ˜ 2 t overlap, and the corresponding limit cycles also overlap, indicating that the coupling stiffness K ˜ s has no effect on self-oscillation in in-phase mode. This is because when the LCE oscillators are in in-phase mode, the distance between the two mass blocks remains unchanged, and the spring force is a constant force. Figure 5c,d plot the effect of three different coupling stiffnesses on the self-oscillation in the anti-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.65 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . The results show that the larger the coupling stiffness, the smaller the limit cycle and the stronger the suppression of self-excited oscillation. The reason for this phenomenon is because the light energy absorbed from the environment decreases as the coupling stiffness K ˜ s increases in anti-phase mode. Figure 5e–h plot the effect of coupling stiffness on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.25 . The results are similar to the anti-phase mode; that is, the larger the coupling stiffness, the smaller the limit cycle and the stronger the suppression of self-excited oscillation. The change of the amplitude within a doubling period increases as the coupling stiffness K ˜ s increases within a doubling period. More calculations show that the effects of the coupling stiffness are the same in non-phase-locked mode. By calculating the data, it was determined that the self-oscillation frequencies in the anti-phase mode and non-phase-locked mode increase as the coupling stiffness K ˜ s increases.

4.2. Effect of Light Intensity

Figure 6a,b plot the effect of three different light intensities on self-oscillation in in-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . Figure 6c,d plot the effect of three different light intensities on self-oscillation in anti-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.1 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . The light intensities have the same effect on the in-phase mode and the anti-phase mode. The results show that the smaller the light intensity, the smaller the limit cycle, and the stronger the suppression of self-excited oscillation, as shown in Figure 6b,d. This is because the light-driven contraction of the LCE fiber increases as the light intensity increases, and the light energy absorbed by the LCE fiber from the environment increases. This is consistent with the self-excited oscillation of a single oscillator [55]. Figure 6e–h plot the effect of light intensity on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , C 0 = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . It can be seen that the greater the light intensity, the greater the amplitude and the corresponding limit cycle and the greater the change of the amplitude within a doubling period, until it evolves into in-phase mode. The results show that the light intensity has no effect on the frequency of the three synchronization modes.

4.3. Effect of Contraction Coefficient

Figure 7a,b plot the effect of three different contraction coefficients on the self-oscillation in the in-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . Figure 7c,d plot the effect of three different contraction coefficients on the self-oscillation in the anti-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.1 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . The contraction coefficients have the same effect on the in-phase mode and the anti-phase mode. It can be seen form Figure 7a,c that the amplitudes increase as the contraction coefficient increases. The smaller the contraction coefficient and the smaller the corresponding limit cycle are shown in Figure 7b,d. The reason for this phenomenon is that the light-driven contraction of the LCE fiber increases as the contraction coefficient C 0 increases and as the light energy that is absorbed from the environment increases. This is consistent with the self-excited oscillation of a single oscillator [55]. Figure 7e–h plot the effect of the contraction coefficient on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . It can be seen from the figures that the greater the contraction coefficient, the greater the amplitude and the corresponding limit cycle, and the greater the change of the amplitude within a doubling period, until it evolves into in-phase mode. More calculations show that the effect of the contraction coefficient is the same during non-phase-locked mode and that it has no effect on the frequency of the three synchronization modes.

4.4. Effect of Damping Coefficient

Figure 8a,b plot the effect of three different damping coefficients on self-oscillation in the in-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . Figure 8c,d plot the effect of three different damping coefficients on self-oscillation in the anti-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . The damping coefficients have the same effect on the in-phase mode and the anti-phase mode. The results show that the greater the damping coefficient, the smaller the corresponding limit cycle, as shown in Figure 8b,d. This is because that as the damping coefficient increases, the energy consumed by the system increases, and the net energy input during the oscillating decreases. This is consistent with physical tuition and other self-oscillating systems [57]. Figure 8e–h plot the effect of the damping coefficient on self-oscillation in non-phase-locked mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . It can be seen from the figure that the greater the damping coefficient, the smaller the amplitude, limit cycle, and the change of the amplitude, and it then evolves into anti-phase mode until it is in static mode. More calculations show that the frequency of the three synchronization modes does not change with the changes in the damping coefficient.

4.5. Effect of Spring Constant of LCE Fiber

Figure 9a,b plot the effect of three different spring constants of the LCE fiber on the self-oscillation in in-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . Figure 9c,d plot the effect of three different spring constants on self-oscillation in anti-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . The spring constants of LCE fiber K ˜ f have the same effect on in-phase mode and anti-phase mode. The results show that the greater the spring constant of the LCE fiber, the greater the corresponding limit cycle is, as shown in Figure 9b,d. The reason for this phenomenon is that the light energy that is absorbed from the environment increases as the spring constant of the LCE fiber K ˜ f increases. This is consistent with the self-excited oscillation of a single oscillator [56]. Figure 9e–h plot the effect of the spring constant of the LCE fiber on self-oscillation in non-phase-locked mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . It can be seen from the figure that the greater the spring constant of the LCE fiber, the greater the amplitude, the limit cycle, and the change in the amplitude within a doubling period, and it then evolves into anti-phase mode. This is because as the spring constant of the LCE fiber K ˜ f increases, the ratio of the spring force F ˜ s to the driving force F ˜ f 1 F ˜ f 2 decreases, which is equivalent to reducing the coupling stiffness K ˜ s and is more conducive to the realization of anti-phase mode, and this result corresponds to Figure 4. More calculations show that the self-oscillation frequencies of the three synchronization modes increase as the spring constant of the LCE fiber K ˜ f increases.

4.6. Effect of Initial Condition

Figure 10a,b plot the effect of the initial condition on self-oscillation in in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . Figure 10c,d plot the effect of the initial condition on self-oscillation in anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.05 , I ˜ = 0.65 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . Under the action of the different initial conditions, the domains of attraction u ˜ 1 t and u ˜ 2 t overlap, and the corresponding limit cycles also overlap, indicating that the initial conditions have no effect on the self-excited oscillation in in-phase mode and anti-phase mode. This is because the initial condition does not affect the elastic strain energy input to the system. More calculations also show that the influence of initial condition is the same, and this is consistent with the self-excited oscillation of a single oscillator [55,60]. Figure 10e–h plot the effect of initial condition on the self-oscillation in non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . It can be seen from the figures that the smaller the phase difference between the two mass blocks of the initial state, the greater the amplitude of the amplitude change until it evolves into in-phase mode. More calculations show that the initial conditions have no effect on the frequency of the three synchronization modes.

5. Equivalent Systems

Figure 11a,b plot the equivalent systems of the in-phase mode and anti-phase mode, respectively. In in-phase mode, the length of the spring does not change with time, and thus, the system in in-phase mode is equivalent to a single self-excited oscillator with a constant spring force (Figure 11a). Furthermore, it can be easily predicted that the amplitude of the equivalent system increases as the light intensity, contraction coefficient, and spring coefficient of the LCE fiber, and decreases as the damping coefficient increases. The frequency of the self-excited oscillation increases as the spring coefficient of the LCE fiber increases.
In the anti-phase mode, the midpoint of the spring does not move during the oscillation, and thus, the system is equivalent to single self-excited oscillator constrained by a fixed spring with half of its original length, as shown in Figure 11b. Therefore, it can also be predicted that the amplitude of the equivalent system increases as the light intensity, contraction coefficient, and spring coefficient of the LCE fiber increases and decreases as the damping coefficient and spring coefficient of the spring increase. The frequency of self-excited oscillation increases as the spring coefficient of the LCE fiber and spring coefficient of the spring increase.

6. Conclusions

The study of the self-excited oscillation coupling of two or more systems and their collective motion is beneficial to the construction of richer and more complex motions, and the more abundant motion modes are constructed, the richer the functions of micro-robots will be. In this article, we constructed two joint LCE spring oscillators that are connected by a spring and theoretically investigated their synchronization phenomenon based on the well-established dynamic LCE model. The numerical calculations show that the self-oscillation coupling system has three steady synchronization modes: in-phase mode, anti-phase mode, and non-phase-locked mode. The synchronization mode depends on the combination of various parameters, and we provided phase diagrams for three steady synchronization modes under different coupling stiffnesses and different initial states. In addition, the results show that the amplitude of three synchronization modes and the changes in the amplitude within a doubling period in the non-phase-locked mode increased as the light intensity, contraction coefficient, and spring coefficient of the LCE fiber increased and decreased as the damping coefficient decreased. The frequencies of the anti-phase mode and non-phase-locked mode increased as the coupling stiffness and spring coefficients of the LCE fiber increased, and the frequency of the in-phase mode increased as the spring coefficients of the LCE fiber increased. The coupling stiffness and initial conditions had no effect on the in-phase mode. These results are consistent with the self-excited oscillation of a single oscillator, for which we also proposed equivalent systems in the in-phase mode and the anti-phase mode. The coupled self-oscillation system proposed in this paper may be used to design complex systems and have broad application prospects in the fields of signal monitoring, medical equipment, micro-active machines, and micro-robots.

Author Contributions

Conceptualization, K.L., C.D., Y.Y. and Q.C.; Methodology, K.L., C.D. and Y.Y.; Software, C.D.; Validation, C.D., K.L. and Y.Y.; Formal Analysis, C.D. and Y.Y.; Investigation, K.L., Y.Y., C.D. and Q.C.; Resources, Y.Y. and K.L.; Data Curation, C.D., K.L. and Y.Y.; Writing—Original Draft Preparation, C.D. and K.L.; Writing—Review & Editing, C.D., K.L. and Y.Y.; Visualization, C.D. and K.L.; Supervision, Y.Y. and K.L.; Project Administration, Y.Y. and K.L.; Funding Acquisition, Y.Y. and K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by University Natural Science Research Project of Anhui Province under Grant Nos. KJ2020A0449 and KJ2020A0453, National Natural Science Foundation of China under Grant No. 12172001, and The University Synergy Innovation Program of Anhui Province under Grant No. GXXT-2021-019.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Ding, W.J. Self-Excited Vibration; Tsing-Hua University Press: Beijing, China, 2009. [Google Scholar]
  2. Jenkins, A. Self-oscillation. Phys. Rep. 2013, 525, 167–222. [Google Scholar] [CrossRef] [Green Version]
  3. Nocentini, S.; Parmeggiani, C.; Martella, D.; Wiersma, D.S. Optically driven soft micro robotics. Adv. Opt. Mater. 2018, 6, 1800207. [Google Scholar] [CrossRef]
  4. Ge, F.; Yang, R.; Tong, X.; Camerel, F.; Zhao, Y. A multifunctional dyedoped liquid crystal polymer actuator: Light-guided transportation, turning in locomotion, and autonomous motion. Angew. Chem. Int. Ed. 2018, 57, 11758. [Google Scholar] [CrossRef]
  5. Lindsey, H.; Kirstin, P.; Guo, Z.L.; Metin, S. Soft Actuators for Small-Scale Robotics. Adv. Mater. 2017, 29, 1603483. [Google Scholar]
  6. Wang, Z.; Xu, W.; Wang, Z.; Lyu, D.; Mu, Y.; Duan, W.; Wang, Y. Polyhedral Micromotors of Metal-Organic Frameworks: Symmetry Breaking and Propulsion. J. Am. Chem. Soc. 2021, 143, 19881–19892. [Google Scholar] [CrossRef] [PubMed]
  7. Papangelo, A.; Putignano, C.; Hoffmann, N. Self-excited vibrations due to viscoelastic interactions. Mech. Syst. Signal Process. 2020, 144, 106894. [Google Scholar] [CrossRef]
  8. Huang, H.; Aida, T. Towards molecular motors in unison. Nat. Nanotechnol. 2019, 14, 407. [Google Scholar] [CrossRef]
  9. Chatterjee, S. Self-excited oscillation under nonlinear feedback with time-delay. J. Sound Vib. 2011, 330, 1860–1876. [Google Scholar] [CrossRef]
  10. Rothemund, P.; Ainla, A.; Belding, L.; Preston, D.J.; Kurihara, S.; Suo, Z.; Whitesides, G.M. A soft bistable valve for autonomous control of soft actuators. Sci. Robot. 2018, 3, eaar7986. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, L.; Chang, L.; Hu, Y.; Huang, M.; Ji, Q.; Lu, P.; Liu, J.; Chen, W.; Wu, Y. An autonomous soft actuator with light-driven self-sustained wavelike oscillation for phototactic self-locomotion and power generation. Adv. Funct. Mater. 2020, 30, 1908842. [Google Scholar] [CrossRef]
  12. Byun, J.; Lee, Y.; Yoon, J.; Lee, B.; Oh, E.; Chung, S.; Lee, T.; Cho, K.; Kim, J.; Hong, Y. Electronic skins for soft, compact, reversible assembly of wirelessly activated fully soft robots. Sci. Robot. 2018, 3, 9020. [Google Scholar] [CrossRef] [Green Version]
  13. Yang, J.C.; Mun, J.; Kwon, S.Y.; Park, S.; Bao, Z.; Park, S. Electronic skin: Recent progress and future prospects for skin-attachable devices for health monitoring, robotics, and prosthetics. Adv. Mater. 2019, 31, 1904765. [Google Scholar] [CrossRef] [Green Version]
  14. Rihani, R.T.; Kim, H.; Black, B.J.; Atmaramani, R.; Saed, M.O.; Pancrazio, J.J.; Ware, T.H. Liquid crystal elastomer-based microelectrode array for in vitro neuronal recordings. Micromachines 2018, 9, 416. [Google Scholar] [CrossRef] [Green Version]
  15. Lu, H.; Zou, Z.; Wu, X.; Shi, C.; Liu, Y.; Xiao, J. Biomimetic Prosthetic Hand Enabled by Liquid Crystal Elastomer Tendons. Micromachines 2021, 12, 736. [Google Scholar] [CrossRef]
  16. Ahn, C.; Li, K.; Cai, S. Light or Thermally Powered Autonomous Rolling of an Elastomer Rod. ACS. Appl. Mater. Interfaces 2018, 10, 25689–25696. [Google Scholar] [CrossRef]
  17. Serak, S.; Tabiryan, N.; Vergara, R.; White, T.J.; Vaia, R.A.; Bunning, T.J. Liquid crystalline polymer cantilever oscillators fueled by light. Soft Matter 2010, 6, 779–783. [Google Scholar] [CrossRef]
  18. Yamada, M.; Kondo, M.; Mamiya, J.; Yu, Y.; Kinoshita, M.; Barrett, C.J.; Ikeda, T. Photomobile Polymer Materials: Towards Light-Driven Plastic Motors. Angew. Chem. Int. Ed. 2008, 47, 4986–4988. [Google Scholar] [CrossRef]
  19. He, Q.G.; Wang, Z.J.; Wang, Y.; Wang, Z.J.; Li, C.H.; Annapooranan, R.; Zeng, J.; Chen, R.K.; Cai, S. Electrospun liquid crystal elastomer microfiber actuator. Sci. Robot. 2021, 6, eabi9704. [Google Scholar] [CrossRef]
  20. Li, K.; Su, X.; Cai, S. Self-sustained rolling of a thermally responsive rod on a hot surface. Extreme Mech. Lett. 2021, 42, 101116. [Google Scholar] [CrossRef]
  21. Wu, J.; Yao, S.; Zhang, H.; Man, W.; Bai, Z.; Zhang, F.; Wang, X.; Fang, D.; Zhang, Y. Liquid Crystal Elastomer Metamaterials with Giant Biaxial Thermal Shrinkage for Enhancing Skin Regeneration. Adv. Mater. 2021, 33, 2106175. [Google Scholar] [CrossRef]
  22. Sun, B.; Jia, R.; Yang, H.; Chen, X.; Tan, K.; Deng, Q.; Tang, J. Magnetic Arthropod Millirobots Fabricated by 3D-Printed Hydrogels. Adv. Intell. Syst. 2021, 4, 2100139. [Google Scholar] [CrossRef]
  23. He, X.; Aizenberg, M.; Kuksenok, O.; Zarzar, L.D.; Shastri, A.; Balazs, A.C.; Aizenberg, J. Synthetic homeostatic materials with chemo-mechano-chemical self-regulation. Nature 2012, 487, 214–218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Shin, A.; Ha, J.; Lee, M.; Park, K.; Park, G.H.; Choi, T.H.; Cho, K.J.; Kim, H.Y. Hygrobot: A self-locomotive ratcheted actuator powered by environmental humidity. Sci. Robot. 2018, 3, 2629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Su, H.; Yan, H.; Zhong, Z. Deep neural networks for large deformation of photo-thermo-pH responsive cationic gels. Appl. Math. Model. 2021, 100, 549–563. [Google Scholar] [CrossRef]
  26. Li, K.; Wu, P.Y.; Cai, S. Chemomechanical oscillations in a responsive gel induced by an autocatalytic reaction. J. Appl. Phys. 2014, 116, 043523. [Google Scholar] [CrossRef]
  27. Cheng, Y.; Lu, H.; Lee, X.; Zeng, H.; Priimagi, A. Kirigami-based light-induced shape-morphing and locomotion. Adv. Mater. 2019, 32, 1906233. [Google Scholar] [CrossRef] [Green Version]
  28. Chakrabarti, A.; Choi, G.P.T.; Mahadevan, L. Self-excited motions of volatile drops on swellable sheets. Phys. Rev. Lett. 2020, 124, 258002. [Google Scholar] [CrossRef]
  29. Cheng, Q.; Zhou, L.; Du, C.; Li, K. A light-fueled self-oscillating liquid crystal elastomer balloon with self-shading effect. Chaos Solitons Fractals 2021, 155, 111646. [Google Scholar] [CrossRef]
  30. Vantomme, G.; Elands, L.C.; Gelebart, A.H.; Meijer, E.W.; Pogromsky, A.Y.; Nijmeijer, H.; Broer, D.J. Coupled liquid crystalline oscillators in Huygens’ synchrony. Nat. Mater. 2021, 18, 1702–1706. [Google Scholar] [CrossRef]
  31. Dudkowski, D.; Czołczyński, K.; Kapitaniak, T. Multistability and synchronization: The co-existence of synchronous patterns in coupled pendula. Mech. Syst. Signal Process. 2022, 166, 108446. [Google Scholar] [CrossRef]
  32. O’Keeffe, K.P.; Hong, H.; Strogatz, S.H. Oscillators that sync and swarm. Nat. Commun. 2017, 8, 1504. [Google Scholar] [CrossRef]
  33. Dudkowski, D.; Czołczyński, K.; Kapitaniak, T. Synchronization of two self-excited pendula: Influence of coupling structure’s parameters. Mech. Syst. Signal Process. 2018, 112, 1–9. [Google Scholar] [CrossRef]
  34. Bennett, M.; Schatz, M.F.; Rockwood, H.; Wiesenfeld, K. Huygens’s clocks. Proc. R. Soc. A Math. Phys. Eng. Sci. 2019, 458, 563–579. [Google Scholar] [CrossRef]
  35. Pikovsky, A.; Rosenblum, M.; Kurths, J.; Hilborn, R.C. Synchronization: A Universal Concept in Nonlinear Science. Am. J. Phys. 2002, 70, 655–658. [Google Scholar] [CrossRef]
  36. Vicsek, T.; Zafeiris, A. Collective motion. Physiol. Rep. 2012, 517, 71–140. [Google Scholar] [CrossRef] [Green Version]
  37. Yu, Y. Synchronized dancing under light. Nat. Mater. 2021, 20, 1594–1595. [Google Scholar] [CrossRef]
  38. Huygens, C. Oeuvres Complètes de Christiaan Huygens; Martinus Nijhoff: Hague, The Netherlands, 1899. [Google Scholar]
  39. Ramirez, J.P.; Olvera, L.A.; Nijmeijer, H.; Alvarez, J. The sympathy of two pendulum clocks: Beyond Huygens’ observations. Sci. Rep. 2016, 6, 23580. [Google Scholar] [CrossRef] [Green Version]
  40. Ramirez, J.P.; Fey, R.H.B.; Aihara, K.; Nijmeijer, H. An improved model for the classical Huygens’ experiment on synchronization of pendulum clocks. J. Sound Vib. 2014, 333, 7248–7266. [Google Scholar] [CrossRef]
  41. Zhang, M.; Wiederhecker, G.S.; Manipatruni, S.; Barnard, A.; McEuen, P.; Lipson, M. Synchronization of micromechanical oscillators using light. Phys. Rev. Lett. 2012, 109, 233906. [Google Scholar] [CrossRef] [Green Version]
  42. Dey, S.; Agra-Kooijman, D.M.; Ren, W.; McMullan, P.J.; Griffin, A.C.; Kumar, S. Soft elasticity in main chain liquid crystal elastomers. Crystals 2013, 3, 363–390. [Google Scholar] [CrossRef]
  43. Lan, R.; Sun, J.; Shen, C.; Huang, R.; Zhang, Z.; Zhang, L.; Wang, L.; Yang, H. Near-infrared photodriven self-sustained oscillation of liquid-crystalline network film with predesignated polydopamine coating. Adv. Mater. 2020, 32, 1906319. [Google Scholar] [CrossRef] [PubMed]
  44. Lan, R.; Sun, J.; Shen, C.; Huang, R.; Zhang, Z.; Ma, C.; Bao, J.; Zhang, L.; Wang, L.; Yang, D.; et al. Light-Driven Liquid Crystalline Networks and Soft Actuators with Degree-of-Freedom-Controlled Molecular Motors. Adv. Funct. Mater. 2020, 30, 2000252. [Google Scholar] [CrossRef]
  45. Lu, X.; Guo, S.; Tong, X.; Xia, H.; Zhao, Y. Tunable photo controlled motions using stored strain energy in malleable azobenzene liquid crystalline polymer actuators. Adv. Mater. 2017, 29, 1606467. [Google Scholar] [CrossRef]
  46. Zeng, H.; Wani, O.M.; Wasylczyk, P.; Kaczmarek, R.; Priimagi, A. Self-regulating iris based on light-actuated liquid crystal elastomer. Adv. Mater. 2017, 29, 1701814. [Google Scholar] [CrossRef] [PubMed]
  47. Kumar, K.; Knie, C.; Bléger, D.; Peletier, M.A.; Friedrich, H.; Hecht, S.; Broer, D.J.; Debije, M.G.; Schenning, A.P.H.J.S. Chaotic self-oscillating sunlight-driven polymer actuator. Nat. Commun. 2016, 7, 11975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kuenstler, A.S.; Chen, Y.; Bui, P.; Kim, H.; DeSimone, A.; Jin, L.; Hayward, R.C. Blueprinting photothermal shape-morphing of liquid crystal elastomers. Adv. Mater. 2020, 32, 2000609. [Google Scholar] [CrossRef] [PubMed]
  49. Zhao, J.; Xu, P.; Yu, Y.; Li, K. Controllable vibration of liquid crystal elastomer beams under periodic illumination. Int. J. Mech. Sci. 2020, 170, 105366. [Google Scholar] [CrossRef]
  50. Li, K.; Du, C.; He, Q.; Cai, S. Thermally driven self-oscillation of an elastomer fiber with a hanging weight. Extreme Mech. Lett. 2021, 50, 101547. [Google Scholar] [CrossRef]
  51. Na, Y.H.; Aburaya, Y.; Orihara, H.; Hiraoka, K. Measurement of electrically induced shear strain in a chiral smectic liquid-crystal elastomer. Phys. Rev. E 2011, 83, 061709. [Google Scholar] [CrossRef] [Green Version]
  52. Corbett, D.; Warner, M. Deformation and rotations of free nematic elastomers in response to electric fields. Soft Matter 2009, 5, 1433–1439. [Google Scholar] [CrossRef]
  53. Haber, J.M.; Sanchez-Ferrer, A.; Mihut, A.M.; Dietsch, H.; Hirt, A.M.; Mezzenga, R. Liquid-crystalline elastomer-nanoparticle hybrids with reversible switch of magnetic memory. Adv. Mater. 2013, 25, 1787–1791. [Google Scholar] [CrossRef]
  54. Harris, K.D.; Bastiaansen, C.W.; Lub, J.; Broer, D.J. Self-assembled polymer films for controlled agent-driven motion. Nano Lett. 2005, 5, 1857–1860. [Google Scholar] [CrossRef]
  55. Liang, X.; Chen, Z.; Zhu, L.; Li, K. Light-powered self-excited oscillation of a liquid crystal elastomer pendulum. Mech. Syst. Signal Process. 2022, 163, 108140. [Google Scholar] [CrossRef]
  56. Zeng, H.; Lahikainen, M.; Liu, L.; Ahmed, Z.; Wani, O.M.; Wang, M.; Yang, H.; Priimagi, A. Light-fuelled freestyle self-oscillators. Nat. Commun. 2019, 10, 5057. [Google Scholar] [CrossRef] [Green Version]
  57. Cheng, Q.; Liang, X.; Li, K. Light-powered self-excited motion of a liquid crystal elastomer rotator. Nonlinear Dyn. 2021, 103, 2437–2449. [Google Scholar] [CrossRef]
  58. Li, K.; Chen, Z.; Wang, Z.; Cai, S. Self-sustained eversion or inversion of a thermally responsive torus. Phys. Rev. E 2021, 103, 033004. [Google Scholar] [CrossRef]
  59. Hu, Z.; Li, Y.; Lv, J. Phototunable self-oscillating system driven by a self-winding fiber actuator. Nat. Commun. 2021, 12, 3211. [Google Scholar] [CrossRef]
  60. Xu, P.; Jin, J.; Li, K. Light-powered self-excited bouncing of a liquid crystal elastomer ball. Int. J. Mech. Sci. 2021, 208, 106686. [Google Scholar] [CrossRef]
  61. Gelebart, A.H.; Mulder, D.J.; Varga, M.; Konya, A.; Ghislaine, V.; Meijer, E.W.; Selinger, R.L.B.; Broer, D.J. Making waves in a photoactive polymer film. Nature 2017, 546, 632–636. [Google Scholar] [CrossRef] [Green Version]
  62. Yu, Y.; Nakano, M.; Ikeda, T. Photomechanics: Directed bending of a polymer film by light. Nature 2003, 425, 145. [Google Scholar] [CrossRef]
  63. Finkelmann, H.; Nishikawa, E.; Pereira, G.G.; Warner, M. A new opto-mechanical effect in solids. Phys. Rev. Lett. 2001, 87, 015501. [Google Scholar] [CrossRef]
  64. Brannon, R. Kinematics: The Mathematics of Deformation; Course Notes, ME EN; University of Utah: Salt Lake City, UT, USA, 2008; 39p. [Google Scholar]
  65. Warner, M.; Terentjev, E.M. Liquid Crystal Elastomers; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  66. Hogan, P.M.; Tajbakhsh, A.R.; Terentjev, E.M. UV Manipulation of order and macroscopic shape in nematic elastomers. Phys. Rev. E 2002, 65, 041720. [Google Scholar] [CrossRef] [Green Version]
  67. Torras, N.; Zinoviev, K.E.; Marshall, J.E.; Terentjev, E.M.; Esteve, J. Bending kinetics of a photo-actuating nematic elastomer cantilever. Appl. Phys. Lett. 2011, 99, 254102. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The dynamic model of the self-oscillation coupling system is composed of two identical LCE fibers and a spring. (a) Reference configuration. (b) Pre-stretched state. (c) Current configuration. Under uniform and constant illumination, under uniform and constant illumination, the self-oscillation of the spring oscillators can be triggered by the coupling between the light-driven contraction of the fibers and their movement.
Figure 1. The dynamic model of the self-oscillation coupling system is composed of two identical LCE fibers and a spring. (a) Reference configuration. (b) Pre-stretched state. (c) Current configuration. Under uniform and constant illumination, under uniform and constant illumination, the self-oscillation of the spring oscillators can be triggered by the coupling between the light-driven contraction of the fibers and their movement.
Micromachines 13 00271 g001
Figure 2. Three steady synchronization modes of the self-oscillation coupling system. The other parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0.2 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . (ac) Anti-phase synchronization mode ( K ˜ s = 0.05 ). (df) Non-phase-locked mode ( K ˜ s = 0.15 ). (gi) In-phase synchronization mode ( K ˜ s = 0.5 ).
Figure 2. Three steady synchronization modes of the self-oscillation coupling system. The other parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0.2 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . (ac) Anti-phase synchronization mode ( K ˜ s = 0.05 ). (df) Non-phase-locked mode ( K ˜ s = 0.15 ). (gi) In-phase synchronization mode ( K ˜ s = 0.5 ).
Micromachines 13 00271 g002
Figure 3. Mechanism of the self-excited oscillation in anti-phase mode. The parameters are K ˜ f = 5 , K ˜ s = 0.1 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . The illuminated zone is represented by the shaded area. (a) Time histories of the number fractions of cis-isomers in the two LCE fibers. (b) Time histories of contraction strains. (c) The spring forces of the two LCE fibers (driving forces). (d) The dependence of the spring force of LCE fiber on the mass displacement. (e) The spring force of the spring (spring force); (f) The displacements of the two mass blocks.
Figure 3. Mechanism of the self-excited oscillation in anti-phase mode. The parameters are K ˜ f = 5 , K ˜ s = 0.1 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . The illuminated zone is represented by the shaded area. (a) Time histories of the number fractions of cis-isomers in the two LCE fibers. (b) Time histories of contraction strains. (c) The spring forces of the two LCE fibers (driving forces). (d) The dependence of the spring force of LCE fiber on the mass displacement. (e) The spring force of the spring (spring force); (f) The displacements of the two mass blocks.
Micromachines 13 00271 g003
Figure 4. Dependence of the synchronization mode on the coupling stiffness and initial state. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . Where the purple area is the phase diagram of the in-phase mode, the yellow area is the phase diagram of the anti-phase mode, and the green area is the phase diagram of the anti-phase mode.
Figure 4. Dependence of the synchronization mode on the coupling stiffness and initial state. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . Where the purple area is the phase diagram of the in-phase mode, the yellow area is the phase diagram of the anti-phase mode, and the green area is the phase diagram of the anti-phase mode.
Micromachines 13 00271 g004
Figure 5. Effect of coupling stiffness on the self-oscillation for a constant synchronous mode. (a,b) Effect of coupling stiffness on the in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.2 . (c,d) Effect of coupling stiffness on the anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.65 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . (eh) Effect of coupling stiffness on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.25 . The coupling stiffness has no effect on the self-oscillation in in-phase mode and has an effect on the period and amplitude of the self-oscillation in anti-phase mode and non-phase-locked mode.
Figure 5. Effect of coupling stiffness on the self-oscillation for a constant synchronous mode. (a,b) Effect of coupling stiffness on the in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.2 . (c,d) Effect of coupling stiffness on the anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.65 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.3 . (eh) Effect of coupling stiffness on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 2 0 = 0.25 . The coupling stiffness has no effect on the self-oscillation in in-phase mode and has an effect on the period and amplitude of the self-oscillation in anti-phase mode and non-phase-locked mode.
Micromachines 13 00271 g005
Figure 6. Effect of light intensity on self-oscillation for a constant synchronous mode. (a,b) Effect of light intensity on in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of light intensity on the anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.1 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of light intensity on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , C 0 = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . Light intensity only affects the amplitudes of the three synchronization modes.
Figure 6. Effect of light intensity on self-oscillation for a constant synchronous mode. (a,b) Effect of light intensity on in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of light intensity on the anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.1 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of light intensity on the non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , C 0 = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . Light intensity only affects the amplitudes of the three synchronization modes.
Micromachines 13 00271 g006
Figure 7. Effect of contraction coefficient on the self-oscillation for a constant synchronous mode. (a,b) Effect of three different contraction coefficients on self-oscillation in in-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of three different contraction coefficients on self-oscillation in anti-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of the contraction coefficient on the self-oscillation in the non-phase-locked mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . The contraction coefficient only affects the amplitude of the three synchronization modes.
Figure 7. Effect of contraction coefficient on the self-oscillation for a constant synchronous mode. (a,b) Effect of three different contraction coefficients on self-oscillation in in-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of three different contraction coefficients on self-oscillation in anti-phase mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of the contraction coefficient on the self-oscillation in the non-phase-locked mode. The parameters are K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . The contraction coefficient only affects the amplitude of the three synchronization modes.
Micromachines 13 00271 g007
Figure 8. Effect of the damping coefficient on self-oscillation for a constant synchronous mode. (a,b) Effect of three different damping coefficients on self-oscillation in in-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of three different damping coefficients on the self-oscillation in anti-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of the damping coefficient on the self-oscillation in the non-phase-locked mode for K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . The damping coefficient only affects the amplitude of the three synchronization modes.
Figure 8. Effect of the damping coefficient on self-oscillation for a constant synchronous mode. (a,b) Effect of three different damping coefficients on self-oscillation in in-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . (c,d) Effect of three different damping coefficients on the self-oscillation in anti-phase mode. The parameters are K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of the damping coefficient on the self-oscillation in the non-phase-locked mode for K ˜ f = 5 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.8 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . The damping coefficient only affects the amplitude of the three synchronization modes.
Micromachines 13 00271 g008
Figure 9. Effect of spring constant of LCE fiber on self-oscillation for a constant synchronous mode. (a,b) Effect of three different spring constants on self-oscillation in in-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . (c,d) Effect of three different spring constants on self-oscillation in anti-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of spring constant of the LCE fiber on the self-oscillation in non-phase-locked mode for c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . The spring constants of the LCE fiber have effects on the frequency and amplitude of the three synchronization modes.
Figure 9. Effect of spring constant of LCE fiber on self-oscillation for a constant synchronous mode. (a,b) Effect of three different spring constants on self-oscillation in in-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.1 . (c,d) Effect of three different spring constants on self-oscillation in anti-phase mode. The parameters are c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.6 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of spring constant of the LCE fiber on the self-oscillation in non-phase-locked mode for c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.3 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , u ˙ ˜ 1 0 = 0.3 , and u ˙ ˜ 1 0 = 0.2 . The spring constants of the LCE fiber have effects on the frequency and amplitude of the three synchronization modes.
Micromachines 13 00271 g009
Figure 10. Effect of the initial condition on self-oscillation for a constant synchronous mode. (a,b) Effect of initial condition on self-oscillation in in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . (c,d) Effect of initial condition on the self-oscillation in anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.05 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of initial condition on the self-oscillation in non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . The initial conditions affect the amplitude of the non-phase-locked mode and the change of the amplitude within a doubling period.
Figure 10. Effect of the initial condition on self-oscillation for a constant synchronous mode. (a,b) Effect of initial condition on self-oscillation in in-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . (c,d) Effect of initial condition on the self-oscillation in anti-phase mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.05 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . (eh) Effect of initial condition on the self-oscillation in non-phase-locked mode for K ˜ f = 5 , c ˜ = 0.1 , λ p = 1.15 , K ˜ s = 0.2 , I ˜ = 0.5 , C 0 = 0.7 , u ˜ 1 0 = 0 , u ˜ 2 0 = 0 , and u ˙ ˜ 1 0 = 0.3 . The initial conditions affect the amplitude of the non-phase-locked mode and the change of the amplitude within a doubling period.
Micromachines 13 00271 g010
Figure 11. Equivalent systems of (a) in-phase mode and (b) anti-phase mode. In in-phase mode, the system is equivalent to a single oscillator with a constant force. In anti-phase mode, the system is equivalent to a single oscillator constrained by a fixed spring with half original length.
Figure 11. Equivalent systems of (a) in-phase mode and (b) anti-phase mode. In in-phase mode, the system is equivalent to a single oscillator with a constant force. In anti-phase mode, the system is equivalent to a single oscillator constrained by a fixed spring with half original length.
Micromachines 13 00271 g011
Table 1. Material properties and geometric parameters.
Table 1. Material properties and geometric parameters.
ParameterDefinitionValueUnits
L f Original length of the LCE fiber 0.1 m
K f Internal radius 20 ~ 50 N / m
K s Coupling stiffness 0 ~ 5 N / m
m Mass 0.5 g
c Damping coefficient 5 ~ 8 × 10 3 kg / s
T 0 Thermal relaxation time 10 2 s
η 0 Light-absorption constant 3 × 10 3 1 / s
I Light intensity 1 ~ 3 W / cm 2
C 0 Contraction coefficient 0.4 ~ 0.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Du, C.; Cheng, Q.; Li, K.; Yu, Y. Self-Sustained Collective Motion of Two Joint Liquid Crystal Elastomer Spring Oscillator Powered by Steady Illumination. Micromachines 2022, 13, 271. https://doi.org/10.3390/mi13020271

AMA Style

Du C, Cheng Q, Li K, Yu Y. Self-Sustained Collective Motion of Two Joint Liquid Crystal Elastomer Spring Oscillator Powered by Steady Illumination. Micromachines. 2022; 13(2):271. https://doi.org/10.3390/mi13020271

Chicago/Turabian Style

Du, Changshen, Quanbao Cheng, Kai Li, and Yong Yu. 2022. "Self-Sustained Collective Motion of Two Joint Liquid Crystal Elastomer Spring Oscillator Powered by Steady Illumination" Micromachines 13, no. 2: 271. https://doi.org/10.3390/mi13020271

APA Style

Du, C., Cheng, Q., Li, K., & Yu, Y. (2022). Self-Sustained Collective Motion of Two Joint Liquid Crystal Elastomer Spring Oscillator Powered by Steady Illumination. Micromachines, 13(2), 271. https://doi.org/10.3390/mi13020271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop