Abstract

With the evolution of space exploration, large flexible appendages have been developed in space structures. External perturbation or attitude maneuvering stimulates the vibration (low frequencies) of the aforementioned structures. In this research, a novel joint mechanism was developed in conjunction with active control to inhibit the low-frequency vibration of a large flexible appendage. A compact active joint, based on the idea of electromagnetic direct drive, was designed. The dynamic equations for the large flexible appendage system and the active joint were derived using the Lagrange function with the assumed-modes approach. Single- and multi-frequency excitations were simulated by two noncontact strategies for periodic vibration stimulation along the direction of rotation. The research results revealed that the interference signal had a primary frequency bandwidth of 0.07–0.63 Hz, and the vibration attenuation was prominent between 5.95 and 32.41 dB within the valid bandwidth. Effective inhibition of both the larger and the smaller amplitude vibrations at frequencies lower than 1 Hz could be realized using the proposed active joint without attachment of intelligent materials onto the flexible appendage surface.

1. Introduction

Due to low damping, lightweight, long span, and low stiffness of space flexible structures, the structures inherently exhibit low-frequency nonlinear vibrations. Because temperature shocks or space particles affect the spacecraft during its entry/departure from the Earth shadow or during attitude maneuvering, the vibrations of flexible appendage inevitably exist in structural terms. Typically, the vibration frequencies of such large structures are lower than 1 Hz, and achieving natural damping in an elevated atmosphere and under micro-gravity conditions is difficult. The spacecraft pointing accuracy becomes compromised in case vibrations are not suppressed efficiently. Spacecraft is considerably affected by such low-frequency vibrations. Therefore, the investigation of low-frequency vibrations for flexible spacecraft is critical [16].

The following two major types of control strategies have been reported in literature: the passive control methods based on viscoelastic damping materials and the active control methods based on piezoelectric materials. The implementation of passive control is easy, and its structure is simple [7, 8]. Vibrations can be isolated through appropriate arrangement of damper materials at the appropriate positions of mechanisms/structures. However, passive control cannot satisfactorily suppress abrupt contextual changes, which are influenced by the weights of structure and materials. For various design conditions, passive vibration controllers can become ineffective rapidly. Active vibration isolation techniques differ considerably from passive techniques. In control and modeling literature, the flexible beams driven by motors [911] have been investigated broadly, and because of multipurpose use including accurate tracking control of the beam tip end-effector and suppression of related vibrations, vibration control becomes challenging [1214]. As an intelligent material, piezoelectric is typically used for additional stiffening and damping activities to inhibit unwanted vibration. To improve control performance, numerous active control algorithms have been proposed (e.g., neural network control [15], sliding mode control [16], fuzzy control [1719], robust control [20, 21], adaptive control [22, 23], nonlinear control algorithm [17, 24], and data-driven control [25]). Despite excellent micro-amplitude drive, intelligent materials exhibit a few shortcomings related to large stroke drive and control [2631].

Accordingly, in this study, a novel active joint mechanism is proposed to suppress both the larger and the smaller amplitude vibrations at low frequencies below 1 Hz just by active joint, without attachment of intelligent material onto the flexible appendage surface. The proposed joint mechanism shows a certain novelty for realizing compact structure, micro-amplitude drive, integration of structure and function, and real-time control. Our major innovations are as follows:(a)Electromagnetic Direct Drive. The proposed joint mechanism, based upon the idea of direct driving, exhibits a compact construction and operates in conjunction with an electro-permanent-magnet driving mechanism. Because the spindle is directly linked to the load, the clearance and friction generated by intermediary transmission devices, such as gears and belts, are minimized.(b)Integration of Structural and Functional Characteristics. The joint is placed between adjacent appendages, which simultaneously serve as a functional part and a structural connector. The joint mechanism performs as a structural element, establishing firm connections and supporting neighboring appendages. It also behaves as a functional portion, allowing appendages to deploy, retract, and vibrate.(c)Integration of Deployment/Retraction and Vibration Control. The described joint system has the potential to be employed in flexible appendages, such as solar panels, to create appendages with self-deploying/folding functionalities and vibration control capacity.When alternating current (AC) power is used, deployment and retraction can be actively conducted. The appendage can be spread freely with no impact, and particularly, they can be recycled. Each 60° rotational movement can be achieved using the driving unit with a single impulse-current excitation, thereby facilitating the control easier and making the power supply highly efficient. The panel can be directly driven to deploy/fold for 180° through only three impulse signals in the coils. Moreover, the damping effects can be naturally generated through attraction forces between the rotational permanent magnet and yoke. This damping provides protection from the overshooting of deploying/folding action or vibration.When direct current (DC) power is supplied, the joint mechanism shows varying stiffness characteristics. When the external stimulation is near the system’s resonance frequencies, the joint torsional stiffness can be altered by varying the electric current, resulting in the system’s frequency-shift phenomenon. Consequently, the isolation effect can be obtained. In addition, the active control algorithm is used to actively suppress the low-frequency vibration of the flexible appendage. When the appendage is subjected to periodic excitation disturbance, the joint provides the output torque of same frequency but in reverse direction to reduce the appendage vibration. The real-time angular displacement feedback is employed to control the input current. The output torque of the joint mechanism is proportional to the input current.

2. System Model

2.1. Design of the Active Joint

Figure 1 displays the schematic of the joint mechanism, which mainly comprises a double-magnetic-turntable rotor and stator. In the stator, six coils, which are grouped into EM1, EM2, and EM3, each of which is outfitted using a magnetic yoke. The stator is fixed in part to the housing. The rotor comprises double magnet discs with embedded type layout, which are arranged on the bilateral coil sides symmetrically. For each disc, four permanent magnets are assigned sequentially, and a magnet-plate yoke is used to outfit each disc to mitigate magnetic leakage. The adjoining permanent magnets exhibit magnetic poles of opposite directions. The sector permanent magnet forms a small air gap (only 1.7 mm) with the yoke, which allows the establishment of a uniform and high magnetic flux density at the air gap between two magnetic discs. The joint mechanism actuation is possible using the axial magnetic field through the permanent magnet-electromagnet interaction.

2.2. Modeling of the Joint Driving Torque

In this research, three sets of electromagnetic coils are assumed to be static, whereas the permanent magnet turntables are assumed to be mobile. Accordingly, following Newton’s law of motion, the joint driving torque τ can be formulated as follows [3235]:where I refers to the input current of the coil; Nac and Nturn denote the number and turns of the coil, respectively; indicates the equivalent radius of the magnet plate, where , with R1 and R2 representing the inner and outer magnet radii, respectively; le denotes the equivalent length of the coil; and Bz is the air gap flux density. For the equation of the joint driving torque, the structural parameters of the joint mechanism are illustrated in Figure 2. The specific deduction process of equation (1) is presented in detail in reference [34].

2.3. Modeling of the Large Flexible Appendage System

Consider a large flexible appendage system comprising an active joint and a flexible appendage as displayed in Figure 3. The flexible appendage is assumed to be slender, and the rotary mass moment of inertia of cross section as well as shear deformation was not considered. The deflection of the flexible appendage is deemed small, and the gravity action is ignored. Here, τ is the external torque and θ is the angular displacement of the joint.

In accordance with Lagrange’s method, the following differential equations should be satisfied using a dynamic system that is completely positioned using generalized coordinate qi:where L is the so-called Lagrangian expressed as follows:where T and V denote the kinetic and potential energies of the modeled system, respectively. Furthermore, Qi refers to the generalized external force imposed on the relevant coordinate qi.

Figure 3 illustrates the deformation of P, an arbitrary point (P), on the appendage. The constant inertial frame is represented by the coordinates Oξη [36]. In the local system of coordinates Oxy, which is fixed to the appendage, the position of P0 is x away from the origin of the frame in the Ox direction. The movement of P0 to the P location is noted following the displacements. The location vector of P can be expressed as follows:where denotes P0’s location vector in the Oxy frame and denotes the P0 to P location vector in the same frame, with and u(x, t) representing the axial and transverse deformations of the appendage, respectively. Typically, is considerably smaller than u(x,t). Thus, is assumed to be negligible. In this case, can be considered. The conversion matrix of Oxy to Oξη is expressed as follows:

The velocity of P can be expressed as follows:

Thus, the kinetic energy of the system can be expressed as follows:where ρ denotes the mass density of the plate and J denotes the inertia moment of the joint.

The computational formula for the system’s potential energy is expressed as follows:where E denotes Young’s modulus of the appendage and I denotes its geometrical inertia moment.

By using the assumed modes method, the appendage’s elastic displacement u(x,t) can be described as the linear combination of mode shapes ϕ(x) multiplied by the generalized time-related coordinates q(t) as follows:where ϕ(x) and q(t) denote the 1 × N and N × 1 vectors, respectively, which are expressed as follows:where ϕ(x) is the mode shape function. The boundary conditions of the flexible appendage are as follows.

At x = 0,

At x = l,

For the dynamic equations of the coupling system of the flexible appendage and joint mechanism, both the kinetic energy T presented in equation (7) and the potential energy V given in equation (8) can be expressed using mode shapes together with the generalized coordinates provided in equation (9). Then, by using Lagrange’s equation (equation (2)), the dynamic equations of the coupling system can be obtained from equations (7)–(12):where is the displacement vector; refers to the generalized force; τ(t) represents the output torque of the joint mechanism; and M, K, and h are the mass matrix, stiffening matrix, and nonlinear force, respectively.where the variables of the mass matrix M can be expressed as follows:

3. Vibration Control Simulations

3.1. Simulations of the Joint Driving Torque

Table 1 presents the geometric variables and material traits of the aluminum plate, where the appendage structural damping is disregarded. Table 2 presents the structural variables of the joint mechanism. The simulation was conducted according to the joint driving torque model in Section 2.2. The starting position and the rotation direction of the articulated stator and rotor are shown in Figure 4. The driving torque characteristics of the joint are shown in Figure 5.

The simulation results findings revealed the following: (a)In 0°–40°, when the coils were supplied with DC-6 A power, the maximum output torque was 0.29 Nm. Within this rotation angle range, the joint driving torque remained nearly constant, and the torque was the largest. In addition, the driving torque grew linearly with the increase in the coil input current. The output torque of the joint mechanism was proportional to the input current. Therefore, in 0°–40°, the active control method could be used to actively suppress system vibrations. When the system was subjected to periodic excitation disturbance, the joint was used to output the torque of the same frequency but in the opposite direction to reduce the vibrations of space structures. The real-time angular displacement feedback was used to control the input current.(b)At the angle of 60°, the driving torque was zero, and the joint output torque showed a nearly linear relationship with the rotation angle of 50°–70°, suggesting that the joint torsional stiffness was almost constant. Therefore, when DC supply was provided, the joint showed the characteristic of variable stiffness at this equilibrium position (60°). When the frequency of external stimulation approached the system’s resonance frequency, the joint torsional stiffness could be changed by modulating the electric current, resulting in the frequency-shift phenomenon. Consequently, the isolation effect was realized. Therefore, the joint mechanism offered the potential for semi-active vibration control based on varying stiffness characteristics. However, in this study, the focus was not on semi-active control but on the active vibration control method.

3.2. Simulations of the System Vibration Control

The vibration of space flexible appendage mainly exhibits frequencies lower than 1 Hz. Therefore, in this study, the simulations were conducted to theoretically verify the low-frequency vibration suppression effect of the flexible appendage at less than 1 Hz and to further guide the structural design of the active joint.

To detail the vibration inhibitory efficiency of the active joint, a sinusoidal disturbance torque was utilized. The external disturbance amplitude was set to 0.01 Nm, and the disturbance frequencies were 0.2, 0.4, 0.6, and 0.8 Hz.

To suppress the low-frequency vibration, just simple algorithm for P control is adopted. In the case of the appendage disturbance by sinusoidal stimulation, the input current is controlled based on the real-time feedback of the angular displacement. The expression for the input current of the joint can be I = kp(θd − θ), where kp is the proportional gain of the P controller, and the desired position of the joint is θd = 0. The output torque of the joint mechanism is proportional to the input current, and the relationship between them can be obtained using equation (1).

The comparison of various proportional gains is presented with numerical simulation results. According to the simulation results in Figure 6, the flexible structure vibration can be reduced to small amplitudes. The larger the excitation frequency is, the smaller the amplitude is, and the proposed active joint exhibits a suitable and an efficient performance to inhibit the appendage’s low-frequency vibration just by using the simple P control algorithm. The larger the proportional gain is, the superior the vibration inhibitory efficiency is and the more the energy consumption of the joint mechanism is.

4. Vibration Control Experiments

Spacecraft in space is typically affected by various complex disturbances, resulting in low-frequency vibration of space flexible appendages. Although simulating diverse interference sources during the ground experiment is difficult, equivalent external perturbation in space can be obtained by introducing some typical periodic excitations to study low-frequency vibration isolation performance.

4.1. Vibration Control Experiments under Single-Frequency Excitation

For the functionality validation of the designed active joint mechanism, as displayed in Figure 7, an experimental platform based on rotary flexible appendage was constructed and a self-fabricated stimulator for electromagnetic noncontact periodic excitations was developed along the direction of rotation. One permanent magnet (PM1) was attached to the output shaft of the servo motor and bonded the other (PM2) to the connecting member of the joint mechanism. The two PMs were kept almost 1 mm apart regarding air gap. The motor output angular velocity and displacement based on the PM1–PM2 interaction were regulated, which enabled production of sinusoidal stimulation for experimentation within a certain frequency band, such as varying amplitudes and frequencies.

Figure 8 depicts the block diagram of the experimental system. During the experiments, the low-frequency signal of sinusoidal perturbation, which is stimulated by electromagnetic noncontact periodic excitation, is exerted on the large appendage. The desired angle position θd is 0. The applied current is proportional to the error signal. The input current is controlled using the real-time feedback of angular displacement. If the appendage is interfered with periodic stimulation, the joint inhibits the appendage vibration by outputting identical frequency and reverse directional torque. The sampling period is 1 ms. Two sets of experiments were conducted, namely Experiments I and II.

For the excitations of same frequency but distinct amplitudes in Experiment I, the time domain and frequency domain responses are displayed in Figure 9. Considering Figure 9(a) as an example, the perturbation signal has a primary frequency of 0.26 Hz. For the proportional gain kp of 0.02, the amplitude with control in time domain decreased from –2.40°–2.39° to –1.44°–1.34°. The amplitude with control in the frequency domain was suppressed from 2.26° to 1.05°, with 53.54% isolating effect and 6.66 dB attenuation. For control gain kp of 0.03, the results revealed a decline of the amplitude with control in the temporal domain to –0.27°–0.31° from –2.40°–2.39°. The amplitude with control in the frequency domain is attenuated from 2.26° to 0.16°, which indicates a decrease of 92.92%, corresponding to 23 dB suppression. Subsequently, based on the real-time feedback of angular displacement, the input current was adjusted to achieve active vibration isolation. The three sets of coils had the same input current, and the root mean square of the input current in each coil was 0.431 A. Next, in Figures 9(b) and 9(c), the same analytical procedures were performed. The experiment results are listed in Table 3. The minimum control scope of amplitude exhibits a decrease to –0.27°–0.31° (peak-to-peak 0.58°) from –2.40°–2.39° (peak-to-peak 4.79°), whereas the maximum control scope exhibits repression to –0.67°–1.06° (peak-to-peak 1.73°) from –14.33°–15.28° (peak-to-peak 29.61°). The frequency of perturbation almost remained unchanged (0.26 Hz), and the amplitude varied, with the vibration attenuation of 23–32.41 dB. As indicated by the experimental findings, the vibration prominently diminished merely by the simple P control algorithm, particularly in the case of large-amplitude vibrations at low frequencies. The larger the proportional gain was, the better the isolation effect was, and the larger the current required. The best vibration isolation effect up to 97.60% was achieved in Experiment I (b).

For the excitations of various frequencies and amplitudes in Experiment II, the time- and frequency-domain response are displayed in Figure 10. In Figure 10(a), the perturbation signal has a primary frequency of 0.46 Hz. For proportional gain kp of 0.10, the amplitude with control in time domain decreased from –4.30°–4.33° to –2.78°–2.82°; the amplitude with control in the frequency domain was suppressed from 3.90° to 2.53°, with 35.13% isolating effect and 3.76 dB attenuation. For kp of 0.30, the results revealed a decline in the amplitude with control in the temporal domain from –4.30°–4.33° to –1.32°–0.93°. The amplitude with control in the frequency domain was attenuated from 3.9° to 0.69°, with a decrease of 82.31%, which corresponded to 15.04 dB suppression. The root mean square of the joint input current was 1.479 A. Next, in Figures 10(b) and 10(c), the same analytical procedures were performed. The experiment results are listed in Table 4, the control scope of amplitude exhibits a decrease to –0.30°–0.29° (peak-to-peak 0.59°) from –0.98°–0.90° (peak-to-peak 1.88°), whereas the maximum control scope exhibits the decrease to –1.32°–0.93° (peak-to-peak 2.25°) from –4.30°–4.33° (peak-to-peak 8.63°). The perturbation signal exhibited a primary frequency bandwidth of 0.46–0.63 Hz, whereas, the vibration attenuation ranged between 12.26 and 15.04 dB within the valid bandwidth. The larger the excitation frequency is, the smaller the amplitude is. As detailed by the foregoing findings, the inhibition of vibrations with smaller amplitudes is possible, and the control strategy is simple and easy to implement.

In Section 4.1, the rotational vibration of flexible appendage can be frequently excited by the external noncontact disturbances with single-frequency periodic excitation. The joint input current was adjusted to actively maintain the real-time feedback of angular displacement to suppress the vibration of the appendage. The experimental results revealed that for the single-frequency excitation, even with the simple P control algorithm, the proposed joint mechanism effectively suppressed both the larger and the smaller amplitude vibrations at low frequency. The perturbation signal exhibited a primary frequency bandwidth of 0.26–0.63 Hz, where the vibration attenuation ranged between 12.26 and 32.41 dB in effective bandwidth.

4.2. Vibration Control Experiments under Multi-frequency Stimulations

Two regular shaking head fans with adjustable 2 to 3 power levels are applied for simulating the periodic noncontact stimulations at multi-frequencies, as displayed in Figure 11. During the experiments, the fans were located on the bilateral sides of the flexible appendage. Through differing in modes of combination, varying perturbation traits can be derived. Four cases were provided to illustrate the vibration isolation efficiency of the proposed actuator under two fans airflow periodic excitation, as presented in Table 5. Although this excitation method is simple, the appendages random perturbation was presented accurately.

In Experiment I, as illustrated in Figure 12(a), for kp of 0.05, the amplitude under active control decreased from –11.97°–9.99° to –5.78°–5.84°, the primary frequencies of perturbation signal were 0.08 and 0.17 Hz, respectively, the vibration amplitude was weakened under active control from 3.22° to 2.64°, and 3.17° to 1.84°, with 18.01% (–1.73 dB) and 53.31% (–6.62 dB) isolation effect, respectively. For current kp of 0.20, the amplitude declined from –11.97°–9.99° to –3.13°–2.34° under active control. The vibration amplitude was weakened under active control from 3.22° to 1.07° at 0.08 Hz and 3.17° to 0.48° at 0.17 Hz, with 66.67% (–9.57 dB) and 84.86% (–14.60 dB) isolation effect, respectively. Furthermore, the higher was the proportional gain, the more obvious was the active control effect. In Experiment II, the isolation effect is discussed following the same analytical procedures. Figure 12(b) displays the perturbation signal and control efficiency, and Table 6 details the two groups of experimental results.

According to the single-frequency stimulation strategy, in the process of mutual driving between the two PMs, passive damping effect is produced and vibration isolation may be achieved. Therefore, to prove that the actual isolation effect was achieved by the proposed joint mechanism rather than passive damping, Experiments III and IV were performed. As displayed in Figure 13, the motor output shaft was locked, and PM1 remained stationary.

In Experiment III, as illustrated in Figure 14(a), for kp of 0.05, the amplitude declined under active control from –5.83°–4.76° to –4.14°–3.15°. The vibration amplitude weakened under active control from 1.37° to 1.23° at 0.08 Hz, and 1.06° to 0.96° at 0.17 Hz, showing a decrease of 10.22% (–0.94 dB) and 9.4% (–0.86 dB), respectively. For kp of 0.20, the amplitude under active control decreased from –5.83°–4.76° to –2.18°–1.63°. The vibration amplitude weakened under active control from 1.37° to 0.67° at 0.08 Hz and 1.06° to 0.43°at 0.17 Hz, showing a decrease of 51.09% (–6.21 dB) and 59.43% (–7.84 dB), respectively. The higher was the proportional gain, the more obvious was the active control effect. In Experiment IV, the active control effect was studied using the same analytical procedure. Figure 14(b) displays the perturbation signal and control efficiency, and Table 7 details relevant summarizations. As suggested by the results of fans airflow stimulation-based experimentation, the passive damper was efficient under severe vibrations. Thus, the experimental comparison results revealed the advantages and the robustness of the proposed method.

In Section 4.2, for multi-frequency excitations, four groups of experiments were conducted. The perturbation signal has a primary frequency bandwidth of 0.07–0.17 Hz, where the vibration attenuation is prominent between 5.95 and 27.47 dB within the valid bandwidth. Higher proportional gain indicates efficient vibration inhibition. The optimal efficiency of the vibration inhibition can be up to 95.77% at 0.13 Hz in Experiment II. Moreover, the novel active joint probably allows repression of vibrations with large amplitudes at low frequencies to some extent. Throughout the aforementioned experiments, the simple P control proved to be a promising and versatile control technique based on the proposed active joint mechanism.

For the external disturbances/noises introduced by the bilateral fan, it is difficult to use only a single fan to simulate low-frequency periodic reciprocating disturbances along the rotation direction. Therefore, two fans are attached onto both the sides of the flexible appendage; one is a shaking-head-type fan, and the other fan does not shake its head. The noncontact airflow force generated by the two fans acts on the surface of the flexible appendage, and then, the periodic interference is obtained. Furthermore, different interference characteristics can be obtained by varying fan power combinations.

For these effects, the low-frequency rigid body angular vibration generated by the airflow excitation presents a multi-frequency phenomenon, mainly because the low-frequency rigid body angular vibration of the system is coupled with the elastic vibration generated by the flexible appendage. The first dominant frequency (approximately 0.08 Hz), that is, the frequency of the shaking head fan, represents the low-frequency rigid body angular vibration caused by this fan. However, other dominant frequency components are the noise interferences caused by the elastic vibration generated by the flexible appendage, and the elastic vibration generated by the flexible appendage results from the angular motion of the rigid body of the system. The airflow excitation of the noncontact fans has a certain guiding importance for the experimental method of low- or ultra-low-frequency interference source simulation of space flexible appendages.

5. Conclusions

This study focused on the suppression of low-frequency vibration of a large flexible appendage. A novel active joint based on direct electromagnetic drive was proposed. The joint was mounted at the root of the flexible appendage and provides control actions to suppress both the larger and the smaller amplitude vibrations at low frequency.

In this study, modeling, simulation, and experimental validation of vibration control, based on a novel joint mechanism, were performed. For the proposed system, its motion equations are acquired based on the assumed-modes approach. Numerical predictions were verified through extensive comparisons. Two noncontact strategies were developed for periodic vibration stimulation along the rotation direction. A self-fabricated electromagnetic vibration stimulator was applied to the appendage root for single-frequency periodic stimulation. As displayed by the results, the perturbation signal exhibits a primary frequency bandwidth of 0.07–0.63 Hz, where the vibration attenuation was prominent between 5.95 and 27.47 dB within the valid bandwidth. Vibrations at low frequencies can be suppressed using the active joint, without additional intelligent material onto the appendage surface. Furthermore, the implementation of the P control scheme considerably attenuated vibration over the uncontrolled case. The numerical results and experimental findings revealed that the proposed joint mechanism exhibits robust vibration suppression at low-frequency ranges for large flexible appendage systems.

The prominent features of such a P-like controller include simple control, high reliability, high robustness, and easy implementation for the proposed joint mechanism, which is particularly suitable for aerospace applications.

However, this paper does not focus on the control algorithm but on the complete solution. The main contribution of this study is the method of low-frequency vibration control for space large flexible appendages with active joint mechanisms. This method has the advantage of electromagnetic direct drive. The spindle is directly linked to the load, and thus, the clearance and friction caused by intermediate transmissions are minimized. Furthermore, it presents the benefits of integrated structure and function design. It is not only a structural member but also a functional member. It acts as a structural component, providing support to neighboring appendages. It also behaves as a functional piece, which can achieve appendage vibration control. Furthermore, an active control algorithm was proposed, and the experimental results indicated that the control algorithm is feasible and easy to implement.

Currently, considerable efforts have been devoted to the structural design and ultimate prototyping of the joint mechanism; the P-like control algorithm provides only a preliminary functional verification. In the future research, we will focus on the intelligent control algorithms for the proposed joint mechanism to realize active vibration control of space flexible appendages highly effectively.

Data Availability

The raw data are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was sponsored by Shanghai Sailing Program (No. 20YF1417100) and China Postdoctoral Science Foundation (No. 2021M692066).