Next Article in Journal
Kinetics Analysis and ADRC-Based Controller for a String-Driven Vascular Intervention Surgical Robotic System
Next Article in Special Issue
A THz Waveguide Bandpass Filter Design Using an Artificial Neural Network
Previous Article in Journal
Oscillation Characteristics of an Artificial Cochlear Sensory Epithelium Optimized for a Micrometer-Scale Curved Structure
Previous Article in Special Issue
Dispersion Engineering of Silicon Nitride Microresonators via Reconstructable SU-8 Polymer Cladding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Beating of a Spherical Liquid Crystal Elastomer Balloon under Periodic Illumination

School of Civil Engineering, Anhui Jianzhu University, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(5), 769; https://doi.org/10.3390/mi13050769
Submission received: 7 April 2022 / Revised: 6 May 2022 / Accepted: 10 May 2022 / Published: 13 May 2022
(This article belongs to the Special Issue Micro/Nano Resonators, Actuators, and Their Applications)

Abstract

:
Periodic excitation is a relatively simple and common active control mode. Owing to the advantages of direct access to environmental energy and controllability under periodic illumination, it enjoys broad prospects for application in soft robotics and opto-mechanical energy conversion systems. More new oscillating systems need to be excavated to meet the various application requirements. A spherical liquid crystal elastomer (LCE) balloon model driven by periodic illumination is proposed and its periodic beating is studied theoretically. Based on the existing dynamic LCE model and the ideal gas model, the governing equation of motion for the LCE balloon is established. The numerical calculations show that periodic illumination can cause periodic beating of the LCE balloon, and the beating period of the LCE balloon depends on the illumination period. For the maximum steady-state amplitude of the beating, there exists an optimum illumination period and illumination time rate. The optimal illumination period is proved to be equivalent to the natural period of balloon oscillation. The effect of system parameters on beating amplitude are also studied. The amplitude is mainly affected by light intensity, contraction coefficient, amount of gaseous substance, volume of LCE balloon, mass density, external pressure, and damping coefficient, but not the initial velocity. It is expected that the beating LCE balloon will be suitable for the design of light-powered machines including engines, prosthetic blood pumps, aircraft, and swimmers.

1. Introduction

Liquid Crystal Elastomer (LCE) is a responsive material synthesized from anisotropic rod-like liquid crystal molecules and long-chain stretch polymers [1,2,3,4,5,6], which can respond to external stimuli, such as thermal [7,8,9], electric [10,11], optical [12,13,14], magnetic [15] and chemical [16,17,18,19] fields. LCE generally enjoys the significant advantages of quick response to deformation, deformation recoverability and noiselessness [20], and it has broad application prospects in many fields, including artificial muscles [21,22,23], microsystems and MEMS [24,25,26,27], actuators and sensors [28,29,30,31,32,33], energy harvesters [34,35], medical devices [36,37], and telescopic optical devices [38,39].
The oscillation of plates and shells caused by electric, magnetic, optical, and thermal excitations is widely used in aerospace, military and civil engineering [40,41,42,43,44,45,46,47,48]. In recent years, light-powered oscillation of LCE has gained huge interests. It can convert light energy into mechanical energy for use in soft robots, micro factories and nanomachines without conventional motor drives [49,50,51,52]. The light-powered oscillation is caused by the large deformation of LCE due to the orientation sequence change in liquid crystal molecule [53,54,55,56,57]. For example, ultraviolet light can induce the photoisomerization of azobenzene molecules, which leads to the large deformation of LCE [58,59]. Finkelmann et al. synthesized LCE containing azobenzene molecules with a photoinduced contraction strain close to 20% [60]. Zhao et al. theoretically studied the forced oscillation of a photosensitive liquid crystal elastomer cantilever beam, and realized the control of the oscillation characteristics by adjusting the illumination period and rate [61]. Cheng et al. proposed a light-powered self-oscillating LCE balloon with self-shadowing coating, and investigated the dynamic process of the self-sustained oscillating LCE balloon based on the dynamic LCE model and ideal gas model [62].
Although many experiments and theoretical studies have been carried out to investigate the static deformation and self-sustained oscillation of LCE induced by static light sources, there are scant research works on the dynamic response of LCE under periodic illumination. In this paper, we propose a model for a spherical LCE balloon under periodic illumination, establish the governing equation of motion, and then further study the beating of the LCE balloon under periodic illumination. We can achieve the dynamic deformation and its deformation response time of the LCE balloon by controlling the light period and light time of periodic illumination. Theoretical calculations show that there exists an optimal lighting period and lighting time rate that maximizes the amplitude of the beating. We also extensively explore the effects of light intensity, contraction coefficient, amount of gaseous substance, volume of LCE balloon, mass density, external pressure, damping coefficient, and initial velocity on the amplitude of the beating. In Section 2, the theoretical model and the governing equation of motion for LCE balloon are introduced. In Section 3, the optimal illumination period and illumination time rate are analyzed, and the effects of each parameter on the beating amplitude of the system are discussed. In Section 4, the results are collated to draw the conclusions.

2. Theoretical Model and Formulation

2.1. Dynamics of the Spherical LCE Balloon

Figure 1 depicts a spherical light-powered beating LCE balloon under periodic illumination. The light illumination is on/off with amplitude I 0 and zero. We assume that the azobenzene liquid-crystal molecules are planar anchored and randomly distributed in the LCE membrane of the balloon. According to Yu et al. [63], ultraviolet light or laser with wavelength less than 400 nm can induce the trans-to-cis isomerization of LCE. When the light shines on the surface of the balloon, a driving force is generated that causes the balloon to contract. When the light is off, the driving force gradually decreases, and the balloon gradually expands. The periodic light illumination may drive the periodic vibration of the LCE balloon. The LCE balloon without stress is taken as the reference state, in which the corresponding radius, thickness and mass density of LCE balloon are denoted by r 0 , h 0 , ρ , respectively. We suppose that the material of the LCE balloon is incompressible and its volume is always V L . Then, the LCE balloon is filled with gas with amount of substance n I , inducing the radius of the balloon increased to r I , which is treated as the initial state. Afterwards, the LCE balloon is placed in the periodic light region, and its radius will vary with time under the combined action of the periodic light and the inner pressure of the balloon. LCE balloon will contract under illumination and recover under non-illumination conditions, and the instantaneous radius of the LCE balloon is described as r ( t ) . Considering the incompressibility of the material, it is exactly in this sense that the thickness of the LCE balloon is regarded to be much smaller than the radius, and thus it can be defined as h = V L / 4 π r 2 .
For simplicity, the influence of gravity is ignored. Since the sphere is completely symmetrical, we analyze the force condition per unit volume. A spherical shell volume element in the LCE balloon is selected; through force analysis, the element is subjected to the gas pressure p in inside the balloon, the tensile stress σ in the LCE membrane of the balloon, the damping force and the total external pressure p ext = p load + p am , in which p load indicates the loading pressure and p am denotes the ambient pressure. Considering Newton’s laws of dynamics, the thermodynamic versus kinetic reaction control equation for the LCE element can be depicted as
p in d s 2 p ext d s 2 4 σ h d s 1 2 d s r β d r d t d s 2 = ρ h d 2 r d t 2 d s 2
where 4 σ h d s 1 2 d s r is the force applied by the tensile stress in the LCE membrane of the LCE balloon as shown in Figure 1c, d s is the edge length of the spherical shell volume element of the LCE balloon, and β is the damping coefficient. With the assumption that gas inside the balloon is the ideal gas, the equation of state will be p in V = n I R T e , where R denotes the ideal gas constant, T e represents the thermodynamic temperature of the ideal gas, and V = 4 π r 3 / 3 is the gas volume. As stated above, the expression of the gas pressure inside the LCE balloon is obtained, which is p in = 3 n I R T e / 4 π r 3 . Equation (1) can be reduced to
p in p ext 2 σ h r β d r d t = ρ h d 2 r d t 2
To focus on the oscillation of LCE balloon under periodic illumination, for simplicity, we assume that the stiffness of the LCE is constant during the oscillation, and the stress–strain relation is given as
σ = E eff ε
with E eff = 2 E / 3 indicating the effective Young’s modulus in the equiaxial stress state, and ε = [ r r 0 ( 1 + ε L ) ] / r 0 ( 1 + ε L ) denoting the elastic strain, where ε L is the effective light-induced contraction strain. To keep things simple, the effective light-induced contraction ε L is approximated as being proportional to the number fraction of cis-isomers φ ( t ) , i.e.,
ε L = C 0 φ ( t )
with C 0 being the contraction coefficient. Section 2.2 will provide a further description of the number fraction φ ( t ) .

2.2. Dynamic LCE Model

It is clearly known that after obtaining the number fraction of cis-isomers in the LCE balloon, the light-induced contraction can be calculated, and then the dynamics of the LCE balloon can be studied. Here, we provide a brief summary of the well-established dynamic LCE model. As illustrated in the research work of Yu et al. [63], ultraviolet light or laser with a wavelength less than 400 nm can induce the trans-to-cis isomerization of LCE. With the absorption of light energy, the rod-like trans-isomers isomerize to the bent cis-isomers; meanwhile, due to the thermal effect, some of the cis isomers convert back to trans ones. Hence, the number fraction of cis-isomers is dependent on the light-driven trans-to-cis isomerization, the thermal excitation from trans to cis, and the thermal cis-to-trans relaxation. In general, the thermal excitation from trans to cis is often negligible compared to the light-driven trans-to-cis isomerization. The number fraction of cis-isomers bent in LCE can be described by the following rate equation [63]:
φ t = η 0 I 0 ( 1 φ ) φ τ
with τ representing the thermal cis-to-trans relaxation time, I 0 indicating the light intensity, and η 0 denoting the light absorption constant. By solving Equation (5), the number fraction of cis-isomers can be given as
φ ( t ) = η 0 τ I 0 η 0 τ I 0 + 1 + ( φ 0 η 0 τ I 0 η 0 τ I 0 + 1 ) exp [ t τ ( η 0 τ I 0 + 1 ) ]
where φ 0 is the number fraction of cis-isomers at the initial moment t = 0 of a process. If the initial number fraction of the cis-isomers is zero in the illuminated state, that is, φ 0 = 0 , Equation (6) can be simplified to
φ ( t ) = η 0 τ I 0 η 0 τ I 0 + 1 { 1 exp [ t τ ( 1 + η 0 τ I 0 ) ] }
In addition, for the non-illuminated state, namely I 0 = 0 , Equation (6) can be modified into
φ ( t ) = φ 0 exp ( t τ )
where the undetermined φ 0 can be set to the maximum value of φ ( t ) in Equation (7), which is φ 0 = η 0 τ I 0 / ( η 0 τ I 0 + 1 ) . Therefore, the number fraction of cis-isomers in Equation (8) can be rewritten as
φ ( t ) = η 0 τ I 0 η 0 τ I 0 + 1 exp ( t τ )
It is noted that Equations (7) and (9) provide the evolution laws of the number fraction in illuminated and non-illuminated states, and the current cis number fraction at the conversion moment between the two states does not change, as described in Section 2.4.

2.3. Nondimensionalization

To nondimensionalize the above governing equation, we take the following dimensionless quantities into account: t ¯ = t / τ , I ¯ = η 0 τ I 0 , r ¯ = r / r 0 , r ¯ I = r I / r 0 , p ¯ L = p L / E eff , p ¯ ext = p ext / E eff , n ¯ I = 3 n I R T e / 4 π E eff r 0 3 , V ¯ L = V L / 4 π r 0 3 , β ¯ = β r 0 / E eff τ , ρ ¯ = ρ r 0 2 / E eff τ 2 and φ ¯ = φ ( η 0 τ I 0 + 1 ) / η 0 τ I 0 .
In this paper, the illumination time and darkness time in one cycle are T on and T off , respectively. The light illumination period is defined as T L = T on + T off , and T on / T L denotes the illumination time rate. It is noted that T on and T off are independent of the thermal cis-to-trans relaxation time. Their corresponding dimensionless parameters are T ¯ on = T on / τ , T ¯ off = T off / τ , and T ¯ L = T ¯ L / τ . A is the amplitude of the oscillation of the balloon, and r ˙ ¯ is the dimensionless velocity.
In the illuminated state, the dimensionless expression of Equation (7) is given as
φ ¯ = 1 exp [ t ¯ ( I ¯ + 1 ) ]
Additionally, for the non-illuminated state, the nondimensionalized Equation (9) can be expressed as
φ ¯ = exp ( t ¯ )
By substituting Equations (3) and (4) into Equation (2), the dimensionless governing equation of the LCE balloon can be demonstrated as
d 2 r ¯ d t ¯ 2 = β ¯ r ¯ 2 ρ ¯ V ¯ L d r ¯ d t ¯ + n ¯ I ρ ¯ V ¯ L r ¯ 2 r ¯ 1 + C 0 φ ¯ r ¯ ρ ¯ ( 1 C 0 φ ¯ ) p ¯ ext r ¯ 2 ρ ¯ V ¯ L
where φ ¯ is determined by the evolution laws of Equations (10) and (11), and its evolution is described in detail in Section 2.4.

2.4. Solution Method

Equation (12) presents ordinary differential equations with variable coefficients for which there is no analytic solution. The ordinary differential equations are numerically solved in this paper using the classical fourth-order Runge–Kutta method using MATLAB software. The second-order ordinary differential equation with variable coefficients is transformed into two first-order ordinary differential equations with variable coefficients. Thus, the control equations can be rewritten as
{ d r ¯ ( t ¯ ) d t ¯ = r ˙ ¯ , d 2 r ¯ d t ¯ 2 = f ( t ¯ , r ¯ , r ˙ ¯ ) r ˙ ¯ ( t ¯ = 0 ) = r ˙ ¯ I , r ¯ ( t ¯ = 0 ) = r ¯ I . ,
Consequently, the final steady-state responses of the balloon can be iteratively obtained.
As Figure 2 shows, in an illumination cycle, when the LCE balloon switches between light-on and light-off states, the evolution law is correspondingly converted between Equations (10) and (11). When the LCE balloon is in the illumination state, the number of cis-isomers increases with the passage of time. When the balloon is in the non-illumination state, the number of cis-isomers decreases with the passage of time. It should be noted that at the moment of conversion between light-on and light-off states, the transient cis number fraction φ ¯ ( t ) keeps unchanged. For example, during the first illumination time T ¯ on , φ ¯ ( t ¯ ) rises from coordinate origin to B along the “light on” curve. Next, during the first darkness time T ¯ off , φ ¯ ( t ¯ ) falls from C to D along the “light off” curve. In the second illumination time T ¯ on , φ ¯ ( t ¯ ) transforms from point D to point A.

3. Results and Discussion

In this section, we present a brief discussion on the beating of LCE balloons driven by periodic illumination. Firstly, the optimal illumination period and illumination time rate are numerically determined to maximize the vibration amplitude of the system. Then, under the optimum illumination conditions, the effects of light intensity, contraction coefficient, amount of gaseous substance, volume of LCE balloon, mass density, external pressure and damping coefficient on the vibration amplitude are further investigated. The beating amplitude of the LCE balloon can be adjusted by regulating the periodic lighting, material properties and geometric parameters.

3.1. Dimensionless Parameters

To study the beating of the LCE balloon in detail, we need to estimate the typical values of the dimensionless parameters. From the accessible experiments [3,61,64,65], the material properties and geometric parameters of the system are shown in Table 1, and the dimensionless parameters are estimated as shown in Table 2. Thus, we can calculate the dimensionless light intensity, damping coefficient, amount of gaseous substances, etc. In the following, we explore the beating of balloons with different material properties under different periodic light patterns.

3.2. Optimal Illumination Period

Figure 3 shows the time course diagrams of the LCE balloon for different illumination periods, with the other parameters at T ¯ on / T ¯ L = 0.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . It can be observed from Figure 3 that after the initial non-periodic oscillation, the LCE balloon gradually stabilizes in the vibration. The beatings of the balloon are initially composed of transient vibrations related to free oscillations and steady-state vibrations relating to forced oscillations; then, the free oscillations disappear due to damping, and finally the balloon exhibits only the steady-state response of forced oscillations. It can be seen that the steady-state beating period is equal to the illumination period, and the illumination period has a great influence on the amplitude of the beating. With the different illumination periods, the balloon undergoes two different oscillation states, that is the single peak steady-state oscillation shown in Figure 3a,b, and the multimodal steady-state oscillation shown in Figure 3c,d. Multi-peak oscillation is relatively complex. Here, we only consider the case with a single-peak steady-state oscillation to keep things simple.
Figure 4 presents the variation of the beating amplitude of the balloon with the illumination period at different light time rates. The other parameters in the calculation are set to C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . As can be clearly seen from the figure, the amplitude first increases to the maximum value, and then decreases with increasing illumination period. The LCE balloon beating amplitude reaches its maximum when the illumination period is around 2.7 . It is important to note that the optimal illumination period does not vary with the light time rate.

3.3. Optimal Illumination Time Rate

The variation of the beating amplitude of the balloon with the optimal illumination time rate under the optimal illumination period is plotted in Figure 5. The other parameters are T ¯ L = 2.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing illumination time rate, the beating amplitude first increases from zero to its maximum value, and then decreases back to zero. The maximum amplitude corresponds to the illumination time rate T ¯ on / T ¯ L = 0.5 . The beating amplitude is symmetrical around the time rate of 0.5, and the beating amplitude on the axis of symmetry is the maximum. It is easy to draw the conclusion that by adjusting the illumination time rate, we can control the beating amplitude in soft robot driving and energy acquisition systems.

3.4. Effect of Light Intensity

Figure 6 shows the effect of light intensity I ¯ on the beating of the LCE balloon under the optimal illumination period and optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The time history curves and limit cycles of beating of the balloon under three different light intensities are plotted in Figure 6a,b respectively. As can be seen from the figure, with increasing light intensity, the beating amplitude increases, but the beating period remains unchanged, always being equivalent to the illumination period. The results indicate that the light intensity affects the amplitude of steady-state beating. This is because the light-induced contraction strain of the LCE balloon becomes larger with increasing light intensity, and therefore the amplitude increases.

3.5. Effect of Contraction Coefficient

Figure 7 reveals the influence of contraction coefficient C 0 on the beating of the LCE balloon under the optimal illumination period and optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . Figure 7a plots the time history of the balloon beating under three different contraction coefficients. Figure 7b shows the limit cycles of the LCE balloon under three different contraction coefficients. It can be observed that with increasing contraction coefficient, the beating amplitude shows an increasing trend, but the beating period remains at the same value as the illumination period. The results prove that the contraction coefficient has an influence on the amplitude of steady-state oscillation. This is because with larger values of contraction coefficient, the magnitude of light-induced contraction strain is larger, and thus the amplitude increases. For different contraction coefficients, the LCE balloon can complete beating under the switching between illumination state and non-illumination state.

3.6. Effect of Amount of Substance

Figure 8 depicts the effect of amount of gaseous substance n ¯ I on beating under the optimal illumination period and optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . Figure 8a plots the evolution of the balloon subjected to beatings for different amount of gaseous substance. Figure 8b plots the limit cycles of the LCE balloon for three different amounts of gaseous substance. As can be seen from the graph, the amplitude of the beating increases with increasing amount of gaseous substance; however, the period of beating remains the same as the illumination period. It is indicated that the amplitude of the steady-state oscillations is affected by the amount of gaseous substance. For different amounts of gaseous substance, the LCE balloon can accomplish beating through switching between illuminated and non-illuminated states.

3.7. Effect of LCE Volume

The influence of volume V ¯ L on the beating of the balloon at the optimal illumination period and the optimal time rate, is presented in Figure 9, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . Figure 9a shows the time history of the beating of balloons under different volumes. Figure 9b shows the limit cycles of the LCE balloon under different volumes. As shown in the figure, the LCE balloon presents a beating mode. With increasing balloon volume, the beating amplitude increases; in addition, the beating period remains equal to the illumination period. It can be concluded that the balloon volume has an effect on the amplitude of steady-state oscillation; moreover, for different volumes, the LCE balloon can complete beating under switching between illumination state and non-illumination state.

3.8. Effect of Mass Density

With the other parameters at T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , p ¯ ext = 0.18 , V ¯ L = 0.05 , β ¯ = 0.14 , and r ˙ ¯ I = 0 , Figure 10 presents the effect of mass density ρ ¯ on beating at the optimal illumination period and optimal time rate. The amplitude evolution and limit cycles of the beatings of the balloon under different mass densities are presented in Figure 10a,b, respectively. It can be seen from the figure that the LCE balloon is in the beating mode. With increasing balloon mass density, the beating amplitude increases, but the beating period always remains equivalent to the illumination period. The mass density is proved to be capable of affecting the amplitude of steady-state oscillation. For different mass densities, the LCE balloon can realize beating under the transition between illumination and non-illumination.

3.9. Effect of External Pressure

As shown in Figure 11, the influence of the external pressure p ¯ 0 on beating is illustrated under the optimal illumination period and optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . Figure 11a shows the amplitude evolution of beating of LCE balloon under different external pressures. Figure 11b shows the limit cycles of the LCE balloon under three different external pressures. It can be clearly seen from the figure that the LCE balloon experiences a beating motion. With increasing external pressure, the beating amplitude undergoes a slight decrease; however, the beating period does not change, and it remains equal to the illumination period. The results reveal that the external pressure will have an influence on the amplitude of steady-state oscillation. For different values of external pressure, the beating of the LCE balloon can be completed under switching between the illumination state and non-illumination state.

3.10. Effect of Damping Coefficient

Figure 12 provides the effect of damping coefficient β ¯ on beating under the optimal illumination period and optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and r ˙ ¯ I = 0 . Figure 12a,b show the time history diagrams and limit cycles of beating of balloon under different damping coefficients, respectively. It can be clearly observed from the figure that the LCE balloon is in the beating mode. With increasing damping coefficient, the beating amplitude decreases; additionally, the beating period remains the same, always being equivalent to the illumination period. Obviously, the damping coefficient has an effect on the amplitude of steady-state oscillation. This is due to the energy competition between energy input from illumination and damping-induced energy dissipation. With increasing damping coefficient, the energy dissipation induced by damping increases, resulting in a decrease in beating amplitude.

3.11. Effect of Initial Velocity

Figure 13 indicates the effect of initial velocity r ˙ ¯ I on the beating of the balloon under the optimal illumination period and the optimal time rate, with the other parameters set to T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and β ¯ = 0.14 . Figure 13a shows the evolution of the beating of the balloon at different initial velocities under the optimal illumination period and optimal time rate. Figure 13b shows the limit cycles of the LCE balloon at three different initial velocities. It can be clearly seen that the LCE balloon is in the beating mode, and the beating amplitude does not change with increasing initial velocity. Meanwhile, the beating period remains equal to the illumination period. The results show that the initial velocity has no influence on the amplitude of beating, which is consistent with the general characteristics of forced oscillation.

3.12. An Application Example of the Periodic Oscillation of the Balloon

The oscillating system proposed in this paper has the potential to be developed as a light-fueled micro-pump. The LCE balloon could be placed in a sealed transparent glass sphere with a check valve at each end. When the LCE balloon expands, the water in the glass sphere is squeezed out from the outlet check valve, and when the LCE balloon contracts, the water flows into the glass sphere from the inlet check valve. In practical applications, the energy/power density and energy conversion efficiency highly depend on the specific energy conversion processes. For the simple model developed in the current study, during the periodic oscillation of the LCE balloon, the light energy absorbed by the system is used to compensate for the damping energy and to do work on connected devices. For this light-driven micro-pump, the work done by the system on the connected devices can be considered to be the effective work of the pump.
For the typical values of I = 2.5   W / cm 2 , ρ = 10 3   kg / m 3 , r I = 1.164   mm , τ = 0.1   s , β = 8.02   s MPa / m , E = 1   MPa , η 0 = 0.00022   m 2 / s w , C 0 = 0.3 , and the mass density of the fiber ρ f = 10 3   kg / m 3 , the dimensionless parameters are calculated as T ¯ on / T ¯ L = 0.5 , I ¯ = 0.55 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 and r ˙ ¯ I = 0 . In this case, the maximum and minimum values of the balloon radius during beating are numerically calculated to be r ¯ 2 = 1.164 and r ¯ 1 = 1.145 , respectively, and thus the effective work on the external connected equipment during the expansion of the balloon in a period is calculated as W = E eff p ¯ ext 4 3 π r 0 3 ( r ¯ 2 3 r ¯ 1 3 ) = 3.82 × 10 5 J . In addition, the dimensionless period can be numerically calculated as T ¯ = 2.69 . By inserting T on = T ¯ on τ = 0 . 1345   s , the average power of the external work done by the pump is P = W / T on = 2.84 × 10 4   W .

4. Conclusions

The light-driven oscillation of LCE can convert light energy into mechanical energy, which is of great significance for soft robots, microfactories, and nanomachines. In this paper, the active control of spherical LCE balloon beating driven by periodic illumination is discussed theoretically. Based on dynamic LCE model and Newton’s laws of dynamics, the governing equations are derived, and the numerical solution method of the dynamic equations are given. The numerical results validate that LCE balloon can beat periodically under periodic illumination, and the beating period is related to the period of illumination. For the maximum steady-state amplitude of periodic beating, an optimal periodic illumination condition is exhibited. With increasing illumination time rate, the beating amplitude first increases, reaches its maximum, and then decreases. The amplitude of beating can be accurately controlled by adjusting the light intensity, illumination period, illumination time rate, etc. The beating LCE balloon in this study has promising application prospects in the fields of opto-mechanical energy conversion systems and light-fueled machines and equipment, such as cardiac pacemakers, advanced robotics, and so on.

Author Contributions

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

Funding

This work is supported by National Natural Science Foundation of China (Grant No. 12172001), University Natural Science Research Project of Anhui Province (Grant No. KJ2020A0449), Natural Science Foundation of Anhui Province (Grant No. 2008085QA23), and The Postgraduate Research Project of Universities in Anhui Province (Grant No. YJS20210499).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

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

References

  1. 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]
  2. Agrawal, A.; Yun, T.H.; Pesek, S.L.; Chapman, W.G.; Verduzco, R. Shape-responsive liquid crystal elastomer bilayers. Soft Matter 2014, 10, 1411–1415. [Google Scholar] [CrossRef] [PubMed]
  3. Camacho-Lopez, M.; Finkelmann, H.; Palffy-Muhoray, P.; Shelley, M. Fast liquid-crystal elastomer swims into the dark. Nat. Mater. 2004, 3, 307–310. [Google Scholar] [CrossRef] [PubMed]
  4. Jin, L.; Yan, Y.; Huo, Y. A gradient model of light-induced bending in photochromic liquid crystal elastomer and its nonlinear behaviors. Int. J. Nonlin. Mech. 2010, 45, 370–381. [Google Scholar] [CrossRef]
  5. Sawa, Y.; Urayama, K.; Takigawa, T.; DeSimone, A.; Teresi, L. Thermally driven giant bending of liquid crystal elastomer films with hybrid alignment. Macromolecules 2010, 43, 453–464. [Google Scholar] [CrossRef]
  6. Wang, Z.; Cai, S. Recent progress in programming liquid crystal elastomer through dynamic covalent chemistry. J. Mater. Chem. B. 2020, 8, 6610–6623. [Google Scholar] [CrossRef]
  7. Du, C.; Zhang, B.; Cheng, Q.; Xu, P.; Li, K. Thermally driven self-rotation of a hollow torus motor. Micromachines 2022, 13, 434. [Google Scholar] [CrossRef]
  8. Ware, T.H.; McConney, M.E.; Wie, J.J.; Tondiglia, V.P.; White, T.J. Voxelated liquid crystal elastomers. Science 2015, 347, 982–984. [Google Scholar] [CrossRef] [Green Version]
  9. Baumann, A.; Sánchez-Ferrer, A.; Jacomine, L.; Martinoty, P.; Houerou, V.L.; Ziebert, F.; Kulic, I.M. Motorizing fibres with geometric zero-energy modes. Nat. Mater. 2018, 17, 523–527. [Google Scholar] [CrossRef]
  10. Liu, J.; Gao, Y.; Wang, H.; Wang, H.; Poling-Skutvik, R.; Yang, S. Shaping and locomotion of soft robots using filament actuators made from liquid crystal elastomer–carbon nanotube composites. Adv. Intell. Syst. 2020, 2, 1900163. [Google Scholar] [CrossRef] [Green Version]
  11. 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] [PubMed] [Green Version]
  12. Vantomme, G.; Gelebart, A.H.; Broer, D.J.; Meijer, E.W. A four-blade light-driven plastic mill based on hydrazone liquid-crystal networks. Tetrahedron 2017, 73, 4963–4967. [Google Scholar] [CrossRef]
  13. Corbett, D.; Warner, M. Deformation and rotations of free nematic elastomers in response to electric fields. Soft Matter 2009, 5, 1433. [Google Scholar] [CrossRef]
  14. Hao, Z.; 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]
  15. Haberl, 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]
  16. 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]
  17. 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]
  18. Küupfer, J.; Finkelmann, H. Liquid crystal elastomer: Influence of the orientational distribution of the crosslinks on the phase behaviour and reorientation process. Macromol. Chem. Phys. 1994, 195, 1353–1367. [Google Scholar] [CrossRef]
  19. Kai, L.; Wu, P.; Cai, S. Chemomechanical oscillations in a responsive gel induced by an autocatalytic reaction. J. Appl. Phys. 2014, 116, 6379–6380. [Google Scholar]
  20. Chihyung, A.; Li, K.; Cai, S. Light or thermally powered autonomous rolling of an elastomer rod. ACS Appl. Mater. Inter. 2018, 10, 25689–25696. [Google Scholar]
  21. Li, M.H.; Keller, P. Artificial muscles based on liquid crystal elastomers philos. Trans. R. Soc. A 2006, 364, 2763–2777. [Google Scholar]
  22. 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] [PubMed]
  23. Shenoy, D.K.; Thomsen, D.L., III; Srinivasan, A.; Keller, P.; Ratna, B.R. Carbon coated liquid crystal elastomer film for artificial muscle applications. Sensor. Actuat. A-Phys. 2002, 96, 184–188. [Google Scholar] [CrossRef]
  24. Fleischmann, E.K.; Liang, H.L.; Kapernaum, N.; Giesselmann, F.; Lagerwall, J.; Zentel, R. One-piece micropumps from liquid crystalline core-shell particles. Nat. Commun. 2012, 3, 1178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. 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. [Google Scholar] [CrossRef] [PubMed]
  26. Zhao, D.; Liu, Y. Light-induced spontaneous bending of a simply supported liquid crystal elastomer rectangular plate. Phys. Rev. E 2020, 101, 042701. [Google Scholar] [CrossRef] [PubMed]
  27. Zhao, D.; Liu, Y. A prototype for light-electric harvester based on light sensitive liquid crystal elastomer cantilever. Energy 2020, 198, 117351. [Google Scholar] [CrossRef]
  28. Ohm, C.; Brehmer, M.; Zentel, R. Liquid crystalline elastomers as actuators and sensors. Adv. Mater. 2010, 22, 3366–3387. [Google Scholar] [CrossRef]
  29. Lin, G.R.; Lin, C.J.; Lin, C.T. Investigations on an integrated conducting nanoparticle–liquid crystal elastomerlayer. Nanotechnology 2007, 18, 415706–415711. [Google Scholar]
  30. Fan, H.; Li, S. Modeling microtubule cytoskeleton via an active liquid crystal elastomer model. Comp. Mater. Sci. 2015, 96, 559–566. [Google Scholar] [CrossRef]
  31. Xia, Y.; Cedillo-Servin, G.; Kamien, R.D.; Yang, S. Guided folding of nematic liquid crystal elastomer sheets into 3D via patterned 1D microchannels. Adv. Mater. 2016, 28, 9637–9643. [Google Scholar] [CrossRef] [PubMed]
  32. Dawson, N.J.; Kuzyk, M.G.; Neal, J.; Luchette, P.; Palffy-Muhoray, P. Modeling the mechanisms of the photomechanical response of a nematic liquid crystal elastomer. Josa. B 2011, 28, 2134–2141. [Google Scholar] [CrossRef] [Green Version]
  33. Lindsey, H.; Kirstin, P.; Guo, Z.L.; Metin, S. Soft actuators for small-scale robotics. Adv. Mater. 2017, 29, 1603483. [Google Scholar]
  34. Li, C.; Liu, Y.; Huang, X.; Jiang, H. Direct sun-driven artificial heliotropism for solarenergy harvesting based on a photo-thermomechanical liquid-crystal elastomer nanocomposite. Adv. Funct. Mater. 2012, 22, 5166–5174. [Google Scholar] [CrossRef]
  35. 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]
  36. Minori, A.F.; He, Q.; Glick, P.E.; Adibnazari, I.; Stopol, A.; Cai, S.; Tolley, M.T. Reversible actuation for self-folding modular machines using liquid crystal elastomer. Smart. Mater. Struct. 2020, 29, 105003. [Google Scholar] [CrossRef]
  37. Turiv, T.; Krieger, J.; Babakhanova, G.; Yu, H.; Shiyanovskii, S.V.; Wei, Q.H.; Lavrentovich, O.D. Topology control of human fibroblast cells monolayer by liquid crystal elastomer. Sci. Adv. 2020, 6, eaaz6485. [Google Scholar] [CrossRef]
  38. Dawson, N.J.; Kuzyk, M.G.; Neal, J.; Luchette, P.; Palffy-Muhoray, P. Cascading of liquid crystal elastomer photomechanical optical devices. Opt. Commun. 2011, 284, 991–993. [Google Scholar] [CrossRef] [Green Version]
  39. Wu, J.; Ye, W.; Wang, Y.; Su, C.Y. Modeling of photo-responsive liquid crystal elastomer actuators. Inform. Sci. 2021, 560, 441–455. [Google Scholar] [CrossRef]
  40. Firouzi, B.; Zamanian, M. The effect of capillary and intermolecular forces on instability of the electrostatically actuated microbeam with t-shaped paddle in the presence of fringing field. Appl. Math. Model. 2019, 71, 243–268. [Google Scholar] [CrossRef]
  41. Karimipour, I.; Beni, Y.T.; Akbarzadeh, A.H. Size-dependent nonlinear forced vibration and dynamic stability of electrically actuated micro-plates. Commun. Nonlinear. Sci. 2019, 78, 104856. [Google Scholar] [CrossRef]
  42. Montazami, R.; Spillmann, C.M.; Naciri, J.; Ratna, B.R. Enhanced thermomechanical properties of a nematic liquid crystal elastomer doped with gold nanoparticles. Sensor. Actuat. A-Phys. 2012, 178, 175–178. [Google Scholar] [CrossRef]
  43. Zhang, J.; Guo, Y.; Hu, W.; Soon, R.H.; Davidson, Z.S.; Sitti, M. Smart actuators: Liquid crystal elastomer-based magnetic composite films for reconfigurable shape-morphing soft miniature machines. Adv. Mater. 2021, 33, 2170054. [Google Scholar] [CrossRef]
  44. Dunn, M.L. Photomechanics of mono- and polydomain liquid crystal elastomer films. J. Appl. Phys. 2007, 102, 013506. [Google Scholar] [CrossRef]
  45. Bhattacharyya, A.; Smith, T.L.; Anderson, A.C. Low temperature thermal conductivity and specific heat of elastomers. J. Non-Cryst. Solids. 1979, 31, 395–400. [Google Scholar] [CrossRef]
  46. Jiang, W.; Niu, D.; Liu, H.; Wang, C.; Zhao, T.; Yin, L.; Shi, Y.; Chen, B.; Ding, Y.; Lu, B. Photoresponsive soft-robotic platform: Biomimetic fabrication and remote actuation. Adv. Funct. Mater. 2014, 24, 7598–7604. [Google Scholar] [CrossRef]
  47. Chambers, M.; Finkelmann, H.; Remškar, M.; Sánchez-Ferrer, A.; Zalar, B.; Žumer, S. Liquid crystal elastomer–nanoparticle systems for actuation. J. Mate. Chem. 2009, 19, 1524–1531. [Google Scholar] [CrossRef]
  48. Fuchi, K.; Ware, T.H.; Buskohl, P.R.; Reich, G.W.; Vaia, R.A.; White, T.J.; Joo, J.J. Topology optimization for the design of folding liquid crystal elastomer actuators. Soft Matter 2015, 11, 7288–7295. [Google Scholar] [CrossRef]
  49. Kuenstler, A.S.; Chen, Y.; Bui, P. Blueprinting photothermal shape morphing of liquid crystal elastomers. Adv. Mater. 2020, 32, 2000609. [Google Scholar] [CrossRef]
  50. Choi, H.J.; Jeong, K.U.; Chien, L.C.; Lee, M.H. Photochromic 3-dimensional actuator based on an uncrosslinked liquid crystal elastomer. J. Mater. Chem. 2009, 19, 7124. [Google Scholar] [CrossRef]
  51. Desimone, A.; Gidoni, P.; Noselli, G. Liquid crystal elastomer strips as soft crawlers. J. Mech. Phys. Solids. 2015, 84, 254–272. [Google Scholar] [CrossRef] [Green Version]
  52. Cheng, Q.; Liang, X.; Li, K. Light-powered self-excited motion of a liquid crystal elastomer rotator. Nonlinear. Dynam. 2021, 103, 2437–2449. [Google Scholar] [CrossRef]
  53. Parrany, M.A. Nonlinear light-induced vibration behavior of liquid crystal elastomer beam. Int. J. Mech. Sci. 2018, 136, 179–187. [Google Scholar] [CrossRef]
  54. Li, M.E.; Lv, S.; Zhou, J.X. Photo-thermo-mechanically actuated bending and snapping kinetics of liquid crystal elastomer cantilever. Smart. Mate. Struct. 2014, 23, 125012–125018. [Google Scholar] [CrossRef]
  55. Schüring, H.; Stannarius, R.; Tolksdorf, C.; Zentel, R. Liquid crystal elastomer balloons. Macromolecules 2001, 34, 3962–3972. [Google Scholar] [CrossRef]
  56. Wang, Y.; Liao, W.; Sun, J.; Nandi, R.; Yang, Z. Bioinspired construction of artificial cardiac muscles based on liquid crystal elastomer fibers. Adv. Mater. 2021, 7, 2100934. [Google Scholar] [CrossRef]
  57. Yusuf, Y.; Sumisaki, Y.; Kai, S. Birefringence measurement of liquid single crystal elastomer swollen with low molecular weight liquid crystal. Chem. Phys. Lett. 2003, 382, 198–202. [Google Scholar] [CrossRef]
  58. Li, K.; Cai, S. Modeling of light-driven bending vibration of a liquid crystal elastomer beam. J. Appl. Mech. 2015, 83, 031009. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, J.; Akomolafe, O.I.; Jiang, J.; Peng, C. Light-actuated liquid crystal elastomer prepared by projection display. Materials 2021, 14, 7245. [Google Scholar] [CrossRef]
  60. Dunn, M.L.; Maute, K. Photomechanics of blanket and patterned liquid crystal elastomer films. Mech. Mater. 2009, 41, 1083–1089. [Google Scholar] [CrossRef]
  61. 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]
  62. Cheng, Q.; Zhou, L.; Du, C.; Li, K. A light-fueled self-oscillating liquid crystal elastomer balloon with self-shading effect. Chaos. Soliton. Fract. 2022, 155, 111646. [Google Scholar] [CrossRef]
  63. Yu, Y.; Nakano, M.; Ikeda, T. Directed bending of a polymer film by light. Nature 2003, 425, 145. [Google Scholar] [CrossRef] [PubMed]
  64. Nägele, T.; Hoche, R.; Zinth, W.; Wachtveitl, J. Femtosecond photoisomerization of cis-azobenzene. Chem. Phys. Lett. 1997, 272, 489–495. [Google Scholar] [CrossRef]
  65. White, T.J.; Tabiryan, N.V.; Serak, S.V.; Hrozhyk, U.A.; Tondiglia, V.P.; Koerner, H.; Bunning, T. A high frequency photodriven polymer oscillator. Soft Matter 2008, 4, 1796–1798. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of a spherical light-powered beating LCE balloon. (a) Reference state in stress-free state. (b) The balloon is inflated by gas with amount of substance n I to the balanced state, namely the initial state. (c) The instant state of the balloon as a function of the radius r ( t ) under periodic lighting conditions.
Figure 1. Schematic diagram of a spherical light-powered beating LCE balloon. (a) Reference state in stress-free state. (b) The balloon is inflated by gas with amount of substance n I to the balanced state, namely the initial state. (c) The instant state of the balloon as a function of the radius r ( t ) under periodic lighting conditions.
Micromachines 13 00769 g001
Figure 2. The evolution laws of the number fraction of cis isomers φ ¯ ( t ¯ ) in the balloon with time t ¯ .
Figure 2. The evolution laws of the number fraction of cis isomers φ ¯ ( t ¯ ) in the balloon with time t ¯ .
Micromachines 13 00769 g002
Figure 3. Time course diagrams of the LCE balloon forced to beat for period illumination (a) T ¯ L = 1 , (b) T ¯ L = 3 , (c) T ¯ L = 5 , and (d) T ¯ L = 7 . The parameters are T ¯ on / T ¯ L = 0.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The illumination period affects the beating amplitude of the balloon.
Figure 3. Time course diagrams of the LCE balloon forced to beat for period illumination (a) T ¯ L = 1 , (b) T ¯ L = 3 , (c) T ¯ L = 5 , and (d) T ¯ L = 7 . The parameters are T ¯ on / T ¯ L = 0.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The illumination period affects the beating amplitude of the balloon.
Micromachines 13 00769 g003
Figure 4. The influence of illumination period on the amplitude of the beating balloon, for different illumination time rates. The parameters are C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude first increases and then decreases with increasing illumination period.
Figure 4. The influence of illumination period on the amplitude of the beating balloon, for different illumination time rates. The parameters are C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude first increases and then decreases with increasing illumination period.
Micromachines 13 00769 g004
Figure 5. Variation of the beating amplitude of the balloon with the illumination time rate under the optimal illumination period. The parameters are T ¯ L = 2.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude first increases and then decreases with increasing illumination time rate.
Figure 5. Variation of the beating amplitude of the balloon with the illumination time rate under the optimal illumination period. The parameters are T ¯ L = 2.7 , C 0 = 0.3 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude first increases and then decreases with increasing illumination time rate.
Micromachines 13 00769 g005
Figure 6. The effect of light intensity I ¯ on the light-powered beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing light intensity, the beating amplitude increases.
Figure 6. The effect of light intensity I ¯ on the light-powered beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing light intensity, the beating amplitude increases.
Micromachines 13 00769 g006
Figure 7. The influence of contraction coefficient C 0 on the beating of the LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing contraction coefficient, the beating amplitude increases.
Figure 7. The influence of contraction coefficient C 0 on the beating of the LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing contraction coefficient, the beating amplitude increases.
Micromachines 13 00769 g007
Figure 8. The effect of amount of gaseous substance n ¯ I on the beating of the LCE balloon. (a) Time histories; (b) Limt cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude of oscillations increases with increasing amount of gaseous substance.
Figure 8. The effect of amount of gaseous substance n ¯ I on the beating of the LCE balloon. (a) Time histories; (b) Limt cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The amplitude of oscillations increases with increasing amount of gaseous substance.
Micromachines 13 00769 g008
Figure 9. The effect of LCE balloon volume V ¯ L on the beating of the LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The beating amplitude increases with increasing balloon volume.
Figure 9. The effect of LCE balloon volume V ¯ L on the beating of the LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , ρ ¯ = 0.16 , p ¯ ext = 0.18 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . The beating amplitude increases with increasing balloon volume.
Micromachines 13 00769 g009
Figure 10. The effect of mass density ρ ¯ on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , p ¯ ext = 0.18 , V ¯ L = 0.05 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing balloon mass density, the beating amplitude increases.
Figure 10. The effect of mass density ρ ¯ on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , p ¯ ext = 0.18 , V ¯ L = 0.05 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing balloon mass density, the beating amplitude increases.
Micromachines 13 00769 g010
Figure 11. The influence of external pressure p ¯ 0 on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing external pressure, the beating amplitude decreases.
Figure 11. The influence of external pressure p ¯ 0 on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , β ¯ = 0.14 , and r ˙ ¯ I = 0 . With increasing external pressure, the beating amplitude decreases.
Micromachines 13 00769 g011
Figure 12. The effect of damping coefficient β ¯ on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and r ˙ ¯ I = 0 . With increasing damping coefficient, the beating amplitude decreases.
Figure 12. The effect of damping coefficient β ¯ on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and r ˙ ¯ I = 0 . With increasing damping coefficient, the beating amplitude decreases.
Micromachines 13 00769 g012
Figure 13. The effect of initial velocity r ˙ ¯ I on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and β ¯ = 0.14 . The initial velocity does not affect the beating amplitude.
Figure 13. The effect of initial velocity r ˙ ¯ I on the beating of LCE balloon. (a) Time histories; (b) Limit cycles. The parameters are T ¯ on / T ¯ L = 0.5 , T ¯ L = 2.7 , I ¯ = 0.5 , C 0 = 0.3 , n ¯ I = 0.45 ( r ¯ I = 1.164 ) , V ¯ L = 0.05 , ρ ¯ = 0.16 , p ¯ ext = 0.18 , and β ¯ = 0.14 . The initial velocity does not affect the beating amplitude.
Micromachines 13 00769 g013
Table 1. Material properties and geometric parameters.
Table 1. Material properties and geometric parameters.
ParameterDefinitionValueUnits
τ Thermal relaxation time0.1s
ρ Mass density1000~1200 kg / m 3
I 0 Light intensity15~35 kW / m 2
β Damping coefficient8~10 s · MPa / m
η 0 Light-absorption constant0.00022 m 2 / s · W
E Young’s modulus1 MPa
r I Initial radius0.001164m
C 0 Contraction coefficient0.3/
Table 2. Dimensionless parameters.
Table 2. Dimensionless parameters.
I ¯ n ¯ I V ¯ L ρ ¯ p ¯ ext β ¯ r ˙ ¯ I
0.3~0.50.4~0.50.045~0.0550.15~0.180.17~0.190.12~0.150~0.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, W.; Cheng, Q.; Du, C.; Dai, Y.; Li, K. Beating of a Spherical Liquid Crystal Elastomer Balloon under Periodic Illumination. Micromachines 2022, 13, 769. https://doi.org/10.3390/mi13050769

AMA Style

Cheng W, Cheng Q, Du C, Dai Y, Li K. Beating of a Spherical Liquid Crystal Elastomer Balloon under Periodic Illumination. Micromachines. 2022; 13(5):769. https://doi.org/10.3390/mi13050769

Chicago/Turabian Style

Cheng, Wenyan, Quanbao Cheng, Changshen Du, Yuntong Dai, and Kai Li. 2022. "Beating of a Spherical Liquid Crystal Elastomer Balloon under Periodic Illumination" Micromachines 13, no. 5: 769. https://doi.org/10.3390/mi13050769

APA Style

Cheng, W., Cheng, Q., Du, C., Dai, Y., & Li, K. (2022). Beating of a Spherical Liquid Crystal Elastomer Balloon under Periodic Illumination. Micromachines, 13(5), 769. https://doi.org/10.3390/mi13050769

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