Heat-Driven Synchronization in Coupled Liquid Crystal Elastomer Spring Self-Oscillators

Self-oscillating coupled machines are capable of absorbing energy from the external environment to maintain their own motion and have the advantages of autonomy and portability, which also contribute to the exploration of the field of synchronization and clustering. Based on a thermally responsive liquid crystal elastomer (LCE) spring self-oscillator in a linear temperature field, this paper constructs a coupling and synchronization model of two self-oscillators connected by springs. Based on the existing dynamic LCE model, this paper theoretically reveals the self-oscillation mechanism and synchronization mechanism of two self-oscillators. The results show that adjusting the initial conditions and system parameters causes the coupled system to exhibit two synchronization modes: in-phase mode and anti-phase mode. The work conducted by the driving force compensates for the damping dissipation of the system, thus maintaining self-oscillation. The phase diagrams of different system parameters are drawn to illuminate the self-oscillation and synchronization mechanism. For weak interaction, changing the initial conditions may obtain the modes of in-phase and anti-phase. Under conditions of strong interactions, the system consistently exhibits an in-phase mode. Furthermore, an investigation is conducted on the influence of system parameters, such as the LCE elastic coefficient and spring elastic coefficient, on the amplitudes and frequencies of the two synchronization modes. This study aims to enhance the understanding of self-oscillator synchronization and its potential applications in areas such as energy harvesting, power generation, detection, soft robotics, medical devices and micro/nanodevices.


Introduction
Self-oscillation refers to the phenomenon where a system generates sustained oscillations or periodic changes without external excitation, due to internal coupling and feedback mechanisms [1][2][3][4][5][6][7]. As a result, self-oscillating systems do not require a continuous energy supply from external sources, reducing energy consumption and system complexity. These systems can be adjusted and controlled by tuning internal parameters and coupling methods. Additionally, self-oscillating systems exhibit great flexibility capable of displaying various oscillatory behaviors such as periodic oscillations [8,9] and chaotic oscillations [10][11][12]. Various feedback mechanisms have been suggested to counteract energy loss attributed to damping dissipation, including the coupling of chemical reactions and large deformations [13][14][15], as well as the self-shading mechanism [16]. Currently, self-oscillation systems are widely used in various scientific and engineering fields, such as sensor technology [10,[17][18][19][20][21][22], soft robots [23,24] and so on.
In recent years, the exploration of active materials has further expanded the possibilities of self-oscillating systems. Through ongoing research and development, scientists continue to discover new active materials with unique properties and enhanced performance, such as dielectric elastomers [25], hydrogels [26,27], ionic gels [13], thermal response polymers [28] and liquid crystal elastomer (LCE) [29][30][31]. These active substances produce different responses when stimulated by light [7], heat [9], electricity [32] and

Model and Theoretical Formulation
In the current section, a coupled self-oscillating system consisting of two LCE fibers and a linear spring under a linear temperature field is proposed. Meanwhile, the governing equations and solution methods of the system are given. Figure 1 illustrates the coupled self-oscillating system within a linear temperature field, which consists of two identical LCE spring oscillators connected by two springs. In the nonstress state, the primary length of the LCE fiber is L 1 and the primary length of the spring is L 2 , as shown in Figure 1a. According to Yakacki et al. [28], LC monomer (RM257) and cross-linking agent (PETMP), etc., are used as raw materials, and LCE fibers can be made by a two-step cross-linking reaction. First, one end of the LCE fiber is fixed, while another end is connected with a spring. The lower end of the spring is connected with another spring through a fixed pulley so that the two LCE fibers can be connected in series. To ensure that the system is force stabilized, the LCE fiber and the spring should be pre-stretched, where the prestretch amount is λ 1 , λ 2 , respectively. In the state of equilibrium, the lengths of the LCE fibers and the springs are λ 1 L 1 and λ 2 L 2 , respectively, as shown in Figure 1b. Then, the equilibrium equation of the system in the non-stress state can be obtained:

Dynamic Model of Two LCE Spring Oscillators
where F s 1 0 and F s 2 0 are the initial elastic forces of the two springs, respectively; F L 1 0 and F L 2 0 are the initial elastic forces of two LCE fibers, respectively, where F s0 = k(λ 2 L 2 − L 2 ), F L0 = K(λ 1 L 1 − L 1 ). k and K are elastic coefficients of spring and LCE fiber, respectively. In this case, we can obtain the relationship between λ 1 and λ 2 , i.e.: where F s 1 0 = F s 1 0 /mg, F s 2 0 = F s 2 0 /mg, F L 1 0 = F L 1 0 /mg, F L 2 0 = F L 2 0 /mg, k = kL 1 /mg, K = KL 1 /mg and L 2 = L 2 /L 1 . When placed in the linear temperature field, LCE fibers begin to oscillate along the vertical direction, in which the displacements of particles 1 and 2 are w 1 (t) and w 2 (t), respectively, as shown in Figure 1c. The force analysis diagram of the two particles is given in Figure 1d, where F s 1 and F s 2 are the elastic force of the two springs, respectively (referred to as spring force); F L 1 and F L 2 are the elastic force of two fibers (hereinafter referred to as the driving force); and F d ( . w 1 ) and F d ( . w 2 ) are the damping force in the process of vibration. To simplify the analysis, we make the assumption that the damping force is directly proportional to the particle's velocity and always acts in the opposite direction to the particle's motion. The dynamic governing equations of the system can then be obtained and can be applied at any time: where ..
dt , and the spring force is Since the system can vibrate continuously without divergence, only nonlinear damping is studied, and it is assumed that: where a 0 and a 1 represent the first and second damping coefficients, respectively.
where )  Figure 1. Schematic diagram of two identical LCE fibers connected by two identical springs within the linear temperature field. (a) Reference state; (b) pre-stretched state; (c) current state; (d) force analysis of the mass. Two coupled LCE oscillators can vibrate synchronously within the linear temperature field.
Since the system can vibrate continuously without divergence, only nonlinear damping is studied, and it is assumed that: where 0 a and 1 a represent the first and second damping coefficients, respectively.

Tension in the LCE Fibers
According to the non-uniform deformation of LCE fiber in a linear temperature field, the Lagrangian coordinate system 1 X , 2 X and Euler coordinate system 1 x , 2 x need to be established by taking the particle at the end of LCE fiber as the origin, as shown in Figure 1a,b. When LCE fiber vibrates, the instantaneous position and displacement of a

Tension in the LCE Fibers
According to the non-uniform deformation of LCE fiber in a linear temperature field, the Lagrangian coordinate system X 1 , X 2 and Euler coordinate system x 1 , x 2 need to be established by taking the particle at the end of LCE fiber as the origin, as shown in Figure 1a,b. When LCE fiber vibrates, the instantaneous position and displacement of a particle can be used as 2). The displacement of the particle is represented by w 1 (t) and w 2 (t), respectively.
We assume that the driving force of LCE fiber is linearly dependent on strain: where K is the elastic coefficient of LCE fiber, and one-dimensional strain ε 1 (X, t), ε 2 (X, t) is given by: We assume that the heat-induced strain ε T (X, t) is linearly related to the temperature difference T(X, t) in LCE fiber: where α represents the coefficient of thermal expansion, α < 0 represents thermal contraction, and α > 0 represents thermal expansion.
Since the driving force F L (t) is uniform and constant in the LCE fiber, it can be obtained by integrating both sides of Equation (5) from 0 to X and combining with Equations (6) and (7); the driving force at the end X = L of the LCE fiber can be written as: Since the temperature field in LCE fiber is unevenly distributed and changes with time, heat exchange occurs between the fiber and its surroundings, resulting in a temperature distribution denoted by T ext (t). For simplicity, there is an assumption that the radius R is much smaller than the length L so that the temperature field in the LCE fiber can be seen as uniform, i.e., T = T(X, t). In this case, the temperature in the fiber can be obtained: where τ = ρ c h indicates the characteristic time, ρ c is the heat capacity per unit length of the fiber, and h is the heat transfer coefficient. Assume that the steady-state temperature field in the environment is linear: where Q refers to the temperature at x = 0 and β represents the gradient of temperature. By defining the following infinitesimal constants: t = t/ L/g, F L = F L /mg, u = u/L, w = w/L, X = X/L, x = x/l, τ = τ/ L/g, K = KL/mg, α = αT L , T = T/T L , T ext = T ext /T L , β = βL/T L and Q = Q/T L (T L is the temperature at x = L).
Thus, the elastic force of LCE fiber can be obtained The solution of the temperature field is [77]: By substituting Equation (12) into Equation (11), the elastic force F L (t) of LCE fiber can be obtained:

Governing Equations
By defining F d = F d /mg, a 0 = a 0 m L g , a 1 = a 1 L m , and combing with Equations (4) and (13), Equation (3) can be rewritten as: Equation (14) is an ordinary differential equation with second-order variable coefficients, which is difficult to obtain its analytic solution. In this case, the classical fourth-order Runge-Kutta method is adopted to solve Equation (14) numerically, and the steady-state response of LCE fiber is obtained, meaning, the time-history curve of oscillation of the system.

Two Modes of Synchronization and Their Mechanisms
In the current section, two synchronization modes, namely in-phase mode and antiphase mode, are proposed according to the dynamic Equation (14), and the self-oscillation mechanism and synchronization mechanism are elaborated in detail.
To better study the synchronization behaviors of two LCE spring oscillators, it is necessary to obtain the typical values of the dimensionless system parameters. According to the existing experiments [52,54,78,79], the actual values of system parameters are summarized in Table 1, and the dimensionless system parameters are calculated in Table 2.

Two Synchronization Modes
The time histories of mass displacements can be obtained by setting system parameters K, α, β, a 0 , a 1 , τ, v 1 0 , v 2 0 . The calculation results show that there are two synchronous modes in the system, namely in-phase mode and anti-phase mode, as shown in Figure 2. In calculation, the system parameter is set to: , the two LCE fibers with the same initial velocity first vibrate in the same direction within the linear temperature field. Then, under the influence of damping, the amplitude of self-oscillation gradually decreases and finally stops on the upper side, as shown in Figure 2a,b. Although fibers convert heat into kinetic energy when heated, the converted kinetic energy does not keep them oscillating. When a 0 = 0.02, v 1 0 = 0.1, v 2 0 = 0.5, the fiber will continue to vibrate in the temperature field, and finally evolve into self-oscillation, as shown in Figure 2c,d. In this case, the energy obtained from the temperature field is greater than the damping dissipation, so the self-oscillation is guaranteed. When a 0 = 0.2, v 1 0 = 0.1, v 2 0 = −0.5, the system can maintain the static mode of anti-phase mode, as shown in Figure 2e,f. As a 0 = 0.02, Figure 2g,h plot the displacement-time diagram and phase trajectory diagram in the anti-phase mode. A similar experimental phenomenon was reported by Ghislaine et al. [69], where two liquid crystal network oscillators interacted with each other driven by light and underwent synchronized in-phase and anti-phase oscillations in the steady state.
. There exist two synchronous modes of the system, namely, in-phase mode and antiphase mode.

Self-Oscillation Mechanism
To further investigate the mechanism of self-oscillation of LCE fiber within the linear temperature field, Figure 3a,b plot the time-history curves of LCE fiber for in-phase and anti-phase modes, indicating that two LCE fibers oscillate periodically within the temperature field in in-phase and anti-phase modes. Figure 3c,e plot the curve of the tension of LCE fiber and spring changing with time in in-phase mode, indicating that the tension of LCE fiber and spring change periodically. Figure 3d,f plot the time-varying curves of the driving force and spring force in anti-phase mode, which indicate that the tension of LCE fiber and spring also maintain periodic changes in anti-phase mode. Figure 3g,i show that in the in-phase mode, the LCE fiber and spring tension, along with the displacement, form hysteresis loops, and the region surrounded by the hysteresis loops represents the work done by the LCE fiber tension and spring force. The work done by the driving force of LCE fiber represents the energy input of the system, while the work done by the spring represents the work expended by the resistance. When the energy gain There exist two synchronous modes of the system, namely, in-phase mode and anti-phase mode.

Self-Oscillation Mechanism
To further investigate the mechanism of self-oscillation of LCE fiber within the linear temperature field, Figure 3a,b plot the time-history curves of LCE fiber for in-phase and antiphase modes, indicating that two LCE fibers oscillate periodically within the temperature field in in-phase and anti-phase modes. Figure 3c,e plot the curve of the tension of LCE fiber and spring changing with time in in-phase mode, indicating that the tension of LCE fiber and spring change periodically. Figure 3d,f plot the time-varying curves of the driving force and spring force in anti-phase mode, which indicate that the tension of LCE fiber and spring also maintain periodic changes in anti-phase mode. Figure 3g,i show that in the in-phase mode, the LCE fiber and spring tension, along with the displacement, form hysteresis loops, and the region surrounded by the hysteresis loops represents the work done by the LCE fiber tension and spring force. The work done by the driving force of LCE fiber represents the energy input of the system, while the work done by the spring represents the work expended by the resistance. When the energy gain is equal to the resistance dissipation, the system will maintain self-oscillation. Figure 3h,j draw the hysteresis loops of the driving force of LCE fiber and spring force in anti-phase mode, which refers to the same energy compensation mechanism as the case of the in-phase mode.
; the parameters for the anti-phase mode are . (a,b) Time-history curves for in-phase and anti-phase modes; (c,d) change curve of driving force with time in in-phase and anti-phase modes; (e,f) spring force versus time curves for in-phase and anti-phase modes; (g,h) curves of the work done by the driving force for in-phase and anti-phase modes; (i,j) curves of the work carried out by the spring force for in-phase and anti-phase modes. The energy absorbed by the system from the external environment compensates for the damping dissipation, thus maintaining the self-oscillation of the system.

Synchronization Mechanism
To better study the mechanism of synchronization between two LCE fibers after selfoscillation in a linear temperature field, we plot some key physical quantities in the process of self-oscillation. Figure 4a,b draw the time-history curves for in-phase and antiphase modes. Figure 4c,d, respectively, draw the change curve of the phase difference between fiber 1 and fiber 2 for in-phase and anti-phase modes. Figure 4c,d show that in the in-phase mode, the phase difference gradually decreases until it reaches zero, while in  (e,f) spring force versus time curves for in-phase and anti-phase modes; (g,h) curves of the work done by the driving force for in-phase and anti-phase modes; (i,j) curves of the work carried out by the spring force for in-phase and anti-phase modes. The energy absorbed by the system from the external environment compensates for the damping dissipation, thus maintaining the self-oscillation of the system.

Synchronization Mechanism
To better study the mechanism of synchronization between two LCE fibers after selfoscillation in a linear temperature field, we plot some key physical quantities in the process of self-oscillation. Figure 4a,b draw the time-history curves for in-phase and anti-phase modes. Figure 4c,d, respectively, draw the change curve of the phase difference between fiber 1 and fiber 2 for in-phase and anti-phase modes. Figure 4c,d show that in the in-phase mode, the phase difference gradually decreases until it reaches zero, while in the anti-phase mode, the phase difference finally reaches a fixed value which is equal to half a cycle. the anti-phase mode, the phase difference finally reaches a fixed value which is equal to half a cycle.  Through careful calculation, it is found that when the initial velocity directions of two self-oscillators are the same, the system will always develop into a synchronous mode. However, when the initial velocity direction is opposite, there is a critical LCE elastic coefficient that triggers a transition between in-phase and anti-phase modes. This result is similar to the existing experiment in that the elasticity coefficient can affect the synchronization mode of the system [69]. When the elastic coefficient of LCE is K < 7100, the system can be affected by the initial velocity, where the anti-phase synchronization mode occurs. When K ≥ 7100, the system always leads to an in-phase synchronous mode. In the case of weak interaction, i.e., K < 7100, the system can be likened to being acted on by external forces, as shown in Figure 4e. The system is divided into two LCE selfexciter separately for discussion. In the in-phase mode, each LCE oscillator is equivalent to applying an additional cycle force to another oscillator, which can be expressed as F 1 = A 1 sin(ω 1 t + ϕ 1 ) and F 2 = A 2 sin(ω 2 t + ϕ 2 ). When the periodic force F is consistent with the frequency of the harmonic oscillator, the synchronization phenomenon occurs. The same is true in the anti-phase mode.
Under the weak interaction, Figure 4f-h plots the changes in the system synchronization mode by changing the velocity of L 2 when the initial velocity v 1 0 = 0.05 of L 1 is fixed. As shown in Figure 4f-h, the ring represents the phase change in the process of movement. When the velocity of L 1 is unchanged, the position of L 1 does not change. When the velocity of L 2 is in the blue region, the system can achieve anti-phase mode, because the phases of the two repel each other. It can be seen from Figure 4f-h that an increase in the LCE elasticity coefficient leads to an increase in the synchronization region until, finally, the synchronization region covers all regions.

Parametric Analysis
In Equation (14), there are seven dimensionless parameters: K, k, α, β, a 0 , α 1 and τ, which will affect the motion process of the system. This section discusses the influence of these system parameters on the amplitude and frequency of self-oscillation in in-phase and anti-phase modes.

Effect of LCE Elasticity Coefficient
In Figure 5a,b, the amplitude and frequency of the system change with the change of LCE fiber elasticity coefficient K in in-phase and anti-phase modes. Figure 5a,b show that when K ≤ 6, no matter in which mode, the system always achieves a static state, because the driving force F L is less than the initial elastic force F s0 of the spring, and the LCE oscillator cannot vibrate, i.e., the amplitude and frequency are 0. When K > 6, the amplitude and frequency of the system gradually increase with the increase of K. These results can be understood through the energy input and dissipation of damping. With the increase of K, the driving force F L gradually increases, and the energy supply from the linear temperature field gradually increases, so the amplitude and frequency of the self-oscillation increase. Figure 5c,d draw the limit cycles of different K in two modes. Figure 5c,d show that there is a limit value for the static state and self-excited state in the two synchronous modes, namely K = 6. Figure 6a,b, respectively, draw the curves of the amplitude and frequency changing with different spring elastic coefficients k in the in-phase and the anti-phase modes. Figure 6a shows that in the in-phase mode, with the increase of spring elastic coefficient k, the amplitude of the system gradually decreases and the frequency increases, because with the increase of k, the damping dissipation increases, so the amplitude decreases gradually, while the increase of k makes the spring stiffness increase, so the frequency increases gradually. Figure 6b shows that the amplitude and frequency of the system remain basically unchanged in the anti-phase mode. This is because in inverting mode, the two LCE fibers move in opposite directions and at equal distances, so the total length of the spring connected to the lower end remains the original length, and the amplitude and frequency of the system remain the same. Figure 6c,d show the different limit cycles with different spring elastic coefficients k in two modes. The results show that the system is always in a vibration state in in-phase mode or anti-phase mode, and its motion mode is independent of k.  Figure 6a,b, respectively, draw the curves of the amplitude and frequency changing with different spring elastic coefficients k in the in-phase and the anti-phase modes. Figure 6a shows that in the in-phase mode, with the increase of spring elastic coefficient k , the amplitude of the system gradually decreases and the frequency increases, because with the increase of k , the damping dissipation increases, so the amplitude decreases gradually, while the increase of k makes the spring stiffness increase, so the frequency increases gradually. Figure 6b shows that the amplitude and frequency of the system remain basically unchanged in the anti-phase mode. This is because in inverting mode, the two LCE fibers move in opposite directions and at equal distances, so the total length of the spring connected to the lower end remains the original length, and the amplitude and frequency of the system remain the same. Figure 6c,d show the different limit cycles with different spring elastic coefficients k in two modes. The results show that the system is always in a vibration state in in-phase mode or anti-phase mode, and its motion mode is independent of k .

Effect of Thermal Expansion Coefficient
In Figure 7a,b, the amplitude and frequency of the system change with different thermal expansion coefficients α in the in-phase and the anti-phase modes. It can be seen

Effect of Thermal Expansion Coefficient
In Figure 7a,b, the amplitude and frequency of the system change with different thermal expansion coefficients α in the in-phase and the anti-phase modes. It can be seen from Figure 7a,b that when |α| ≤ 0.2, the system in two modes is in a static state, and the amplitude and frequency are zero. When |α| > 0.2, the amplitude increases with the increase of the absolute value of |α|, and the frequency is unaffected, because the driving force gradually increases with the increase of |α|, so the amplitude gradually increases. Figure 7c,d draw the images of limit cycles changing with different thermal expansion coefficients α in the in-phase and the anti-phase modes. Figure 7 shows that there is a critical value for both static state and vibration state in two modes, namely α = −0.2.    Figure 8a,b plot the variation curves of amplitude and frequency of the system along with temperature gradient β in the in-phase and anti-phase modes. Figure 8a,b show that when β ≤ 0.04, the system is in a static state, where the amplitude and frequency are 0. On the contrary, when β > 0.04, the amplitude of two modes increases with the increase of β, while the frequency is unchanged, which can be understood through energy input and damping dissipation because with the increase of β, the temperature of the temperature field gradually increases, and the driving force F L gradually increases, so the amplitude gradually increases. In Figure 8c,d, the limit cycles change with different temperature gradients β in two modes. It can be seen that there is a critical value β = 0.04 for both static state and self-oscillation state in in-phase and anti-phase modes. Figure 9a,b show the variation curves of the amplitude and frequency with the change of the first damping coefficient a 0 in the in-phase and anti-phase modes. It can be seen from Figure 9a,b that when a 0 ≥ 0.05, the system is in static states, and the amplitude and frequency are 0. When a 0 < 0.05, the amplitude of both the in-phase and the anti-phase modes decreases with the increase of a 0 , while the frequency is unaffected, because with the increase of a 0 , the work carried out by damping increases and the energy dissipation of the system increases, so the amplitude gradually decreases. In Figure 9c,d, the limit cycles change with the change of the first damping coefficient a 0 in the two modes. The results show that the same limit value exists in the static state and vibration state of the in-phase and the anti-phase modes, namely, a 0 = 0.05.    Figure 10a,b show that the amplitude and frequency of the self-oscillation change with the change of the second damping coefficient a 1 in the in-phase and the anti-phase modes. It can be seen from Figure 10 that the amplitude of the system decreases gradually as the second damping coefficient a 1 increases in the in-phase and the anti-phase modes, while the frequency remains unchanged. This is because as a 1 increases, the damping dissipation increases and so the amplitude gradually decreases. Figure 10c,d draw the images of the change of limit cycles with different second damping coefficients a 1 in two modes. The results indicate that the system is always in a vibration state in in-phase and anti-phase modes, and its motion mode is independent of a 1 .

Effect of the Characteristic Time
In Figure 11a  (c,d) limit cycles of in-phase mode and anti-phase mode. In both in-phase mode and anti-phase mode, the amplitude decreases as the second damping coefficient a 1 increases, while the frequency remains constant.

Effect of the Characteristic Time
In Figure 11a,b, the amplitude and frequency change with different characteristic times τ in the in-phase and the anti-phase modes. The results show that in both modes, when τ ≤ 0.06, the system is in a static state and the amplitude and frequency are 0. At τ > 0.06, the amplitude of the self-oscillation increases as τ increases, and the frequency stays the same, because with the increase of τ, the heat transfer rate in LCE fiber increases, resulting in the gradual increase of driving force F L , so the amplitude gradually increases. Figure 11c,d draw the images of limit cycles changing with different temperature gradients τ in in-phase and anti-phase modes. It can be obtained that there is a critical value τ = 0.06 for triggering the self-oscillators in both modes.

Conclusions
The prevalence of synchronization and collective behaviors among self-excited coupled oscillators in nature necessitates investigation due to their inherent benefits, such as efficient energy harvesting, autonomous operation, and enhanced equipment portability. In this paper, based on thermally responsive LCE spring self-oscillators under a linear temperature field, the synchronous behavior of two coupled self-oscillators connected by springs is theoretically investigated. The mechanisms of self-oscillation and synchronization are theoretically investigated, integrating the well-established dynamic LCE model. According to the numerical findings, the system exhibits two synchronization modes: the in-phase mode and the anti-phase mode. Self-oscillations are sustained through a dynamic balance between damped dissipation and work carried out by the driving force.
The numerical findings indicate that the synchronization mode primarily depends on the interaction between two LCE self-oscillators. In cases of strong interaction with the elastic coefficient of LCE, the system consistently develops into the in-phase synchronous mode. However, when the interaction is weak, altering the initial conditions can lead to the in-phase and anti-phase modes. When the initial velocity direction of the two selfoscillators is the same or the initial velocity direction is opposite but the value is small, the system achieves the in-phase synchronous mode. On the contrary, as the initial velocity direction of the two self-oscillators is opposite and the relative value is large, the system evolves into an anti-phase synchronization mode.
In addition, the influences of the LCE elastic coefficient, spring elastic coefficient, thermal expansion coefficient and other system parameters on the synchronous mode, amplitude and frequency of the self-oscillations are systematically studied. The self-oscil-

Conclusions
The prevalence of synchronization and collective behaviors among self-excited coupled oscillators in nature necessitates investigation due to their inherent benefits, such as efficient energy harvesting, autonomous operation, and enhanced equipment portability. In this paper, based on thermally responsive LCE spring self-oscillators under a linear temperature field, the synchronous behavior of two coupled self-oscillators connected by springs is theoretically investigated. The mechanisms of self-oscillation and synchronization are theoretically investigated, integrating the well-established dynamic LCE model. According to the numerical findings, the system exhibits two synchronization modes: the in-phase mode and the anti-phase mode. Self-oscillations are sustained through a dynamic balance between damped dissipation and work carried out by the driving force.
The numerical findings indicate that the synchronization mode primarily depends on the interaction between two LCE self-oscillators. In cases of strong interaction with the elastic coefficient of LCE, the system consistently develops into the in-phase synchronous mode. However, when the interaction is weak, altering the initial conditions can lead to the in-phase and anti-phase modes. When the initial velocity direction of the two selfoscillators is the same or the initial velocity direction is opposite but the value is small, the system achieves the in-phase synchronous mode. On the contrary, as the initial velocity direction of the two self-oscillators is opposite and the relative value is large, the system evolves into an anti-phase synchronization mode.
In addition, the influences of the LCE elastic coefficient, spring elastic coefficient, thermal expansion coefficient and other system parameters on the synchronous mode, amplitude and frequency of the self-oscillations are systematically studied. The selfoscillation amplitude demonstrates a positive correlation with the increase of LCE elastic coefficient, thermal expansion coefficient, temperature gradient and characteristic time, while demonstrates a negative correlation with the increase of spring elastic coefficient and damping coefficient. Unlike existing work about self-oscillating synchronization systems based on active materials [7,70], this paper elucidates in detail the mechanism of the synchronization phenomenon. This study is expected to advance the comprehension of self-oscillator synchronization and provide its potential applications in diverse fields, including energy harvesting, power generation, detection, soft robotics, medical devices and micro/nanodevices. In addition, the research in this paper has the potential to be extended to large-scale synchronization systems containing a large number of coupled oscillators, which has a promising application in this field.