Research Article | Open Access
Vibration Antiresonance Design for a Spacecraft Multifunctional Structure
Spacecraft must withstand rigorous mechanical environment experiences such as acceleration, noise, vibration, and shock during the process of launching, satellite-vehicle separation, and so on. In this paper, a new spacecraft multifunctional structure concept designed by us is introduced. The multifunctional structure has the functions of not only load bearing, but also vibration reduction, energy source, thermal control, and so on, and we adopt a series of viscoelastic parts as connections between substructures. Especially in this paper, a vibration antiresonance design method is proposed to realize the vibration reduction. The complex zero-point equations of the vibration system are firstly established, and then the vibration antiresonance design for the system is achieved. For solving the difficulties due to viscoelastic characteristics of the connecting parts, we present the determining formulas to obtain the structural parameters, so that the complex zero-point equations can be satisfied. Numerical simulation and ground experiment demonstrate the correctness and effectiveness of the proposed method. This method can solve the structural vibration control problem under the function constraints of load bearing and energy supplying and will expand the performance of spacecraft functional modules.
Before getting into orbit and carrying out the normal work, spacecraft should undergo rigorous prelaunch environment. As we know, large level mechanical vibration connected with noise field and structural vibration will be generated during the launch process of rocket, and strong acceleration and severe noise will be accompanied during the ascending period of rocket. Furthermore, intensive mechanical shock will be caused by the satellite-vehicle separation and the operation of initiating explosive device and release mechanism, such as the deployment of solar arrays. Such vibration and shock loads will directly transfer to the spacecraft through the structure and hence bring a potential risk of damage of critical components or structures on spacecraft and accuracy reduction and functional degeneration of spacecraft system or payloads [1, 2]. Therefore, adopting efficient structural vibration control skills to make sure the spacecraft can withstand the rigorous mechanical environment experiences, such as acceleration, noise, vibration, and shock in the boost phase of carrier rocket, is always an inevitable and close concerned problem in the field of aerospace engineering.
In recent two decades, with the development of multichip module technology, high performance composite material technology, flexible circuit connection technology, and advanced thermal control technology, the spacecraft multifunctional structure technology breaks away from cocoon and gets rapid development [3–6]. Such technology organically combines the functions of data processing, radiation protection, thermal control, electric power storage, and so on as a whole and achieves an integrated design of structure, function, and material, which can eliminate the redundant mass and volume of traditional satellite, improve the available free space and payload mass ratio of satellite, shorten the time period of satellite design and fabrication substantially, and extend the life of satellite effectively. Up to date, several researchers have developed such multifunctional structure technology, prompted a series of fundamental researches, and conducted spaceflight test verification on some key techniques. For example, NASA, Lockheed Martin Space Systems, and the U.S. Air Force Research Laboratory have performed a large number of studies on the multifunctional structure, which highly integrates structure, electron, thermal control, and radiation protection, and have produced such multifunctional structure.
In addition, some structural vibration control technologies such as particle damper, impact damper, particle tuned mass damper, and viscoelastic damper [7–10] are used in the aerospace field at present. A reaction wheel isolator assembly  made up of a set of six machined spring elements with bonded viscoelastic material was designed to reduce microvibration transmitted to the optical payload of Chandra Observatory. A 1 Hz isolator  based on constrained layer viscoelastic damping technology is used between the spacecraft bus and the optical payload for the James Webb Space Telescope. An active vibration isolation system , which consists of a piezoelectric stack actuator, is used to isolate sensitive spacecraft payloads from the spacecraft in satellite ultra-quiet isolation technology experiment. Particle impact dampers (PIDs) [14–16] are enclosures/containers or structural voids partially filled with one or more particles (e.g., small balls of metals, sand, ceramics, and glass), which work on by dissipating energy through collisions and frictions between the particles or between the particles and the wall of the enclosure. Such PIDs also have effective vibration absorption effect on the spacecraft.
In summary, the state of the art of spacecraft multifunctional structure technology has taken the functional-integrations such as structure/battery, structure/electron, structure/cable, structure/thermal control, and structure/protection into consideration. Most often the existing vibration control strategies are designed and applied to spacecraft in an independent way. However, the function integration of structure and vibration reduction and the functional-integration design on structure have not been well considered as yet. In our research, we designed a new spacecraft multifunctional structure (MFS) concept, which incorporates the functions of load bearing, vibration reduction, energy source, thermal control, and so on. This paper briefly introduces such a multifunctional structure concept. In the design of such structure, we adopt a series of viscoelastic parts as connections between substructures. Considering the significance of structural vibration reduction for spacecraft or their payloads, a vibration antiresonance design method is especially proposed to realize the vibration reduction in the paper. This proposed method can give a new way to solve the structural vibration control problem under the function constraints of load bearing and energy supplying and will expand the performance of spacecraft functional modules.
The paper is organized as follows. Firstly, the structural composition and functions of the new-concept spacecraft multifunctional structure designed by us are described in brief, and its structural equivalent simplified model is established. Secondly, the complex zero-point equation of the vibration system is derived, and the principle of the vibration antiresonance design method is explained. Thirdly, for solving the difficulties due to viscoelastic characteristics of the connecting parts, the determining formulas of unknown structural parameters are given based on the obtained complex zero-point equation, where the stiffness and damping parameters of the viscoelastic connecting parts are obtained by using the nonlinear mechanical constitutive relation model of the rubber materials. Finally, the effectiveness of the vibration control method proposed in this paper is verified by both the numerical simulation and ground experiment.
2. System and Its Equivalent Simplified Model
In order to improve the payload ratio of spacecraft and increase the available space in the spacecraft, we have proposed a new spacecraft MFS design concept. Such structure mainly consists of a frame structure, four power modules, and multiple viscoelastic connecting parts, as shown in Figure 1.
The frame structure, made of carbon fiber composite materials, is acting as the primary force bearing structure and has the function of load bearing. The power modules use all-solid-state polymer Li-ion batteries and have the functions of electric power supplying and storage. In addition, each power module can behave as an oscillator in a vibration control system to consume the system vibration energy. The viscoelastic connecting parts form a ring-shaped discontinuous cluster as shown in Figure 1. Each viscoelastic connecting parts member can provide an elastic structural support for the power modules and meanwhile provide a good vibration protection for the power modules by generating elastic deformation to dissipate part of the vibration energy. Furthermore, the viscoelastic connecting parts have the capacities of insulation and heat-sinking, since they are made from a kind of high polymer material that is known as silicone rubber.
From Figure 1, we can see that the global configuration of the frame structure is taken on a crisscross shape and thus the frame structure forms four squared cavities. In the center of each squared cavity, there exists a power module. The whole multifunctional structure is symmetric about the x-, y-, and -axis. Here, coordinate system o-xyz denotes the rectangular inertial coordinate system of the multifunctional structure, where the -axis and -axis are perpendicular to each other and are in a plane parallel to the power module surface; the -axis is perpendicular to xoy plane and determined by the right-hand rule. The viscoelastic connecting parts are used to connect the frame structure and the power modules at a precompressed way, and this connection is just surface-contact, which means the viscoelastic connecting parts are not fixed or limited in determined positions.
By embedding the power modules into the interior of structure at an elastic restraint way, the spacecraft multifunctional structure achieves a highly combination of multifunction such as structural load bearing, electric power supplying, vibration reduction, and thermal control. At the same time, by applying a variety of lightweight materials such as carbon fiber composite materials in the fabrication of such multifunctional structure, the total structural mass is reduced drastically. Since the original space used for mounting the electric power systems is totally saved, the available space of spacecraft payloads is improved largely. Consequently, the multifunctional structure technology used in such structure has improved the payload mass ratio, the payload volume ratio, and the function structure ratio of the system platform significantly.
For convenience of subsequent discussions, we just choose a quarter of the structure model shown in Figure 1 as the analysis object according to the axial symmetry properties of the structure. The final-obtained simplified physical model is shown in Figure 2.
In Figure 2, because the stiffness of the power module is far greater than that of the viscoelastic connecting parts and thus the power module almost cause no elastic deformation in the vibration process, the power module could be equivalent as a seismic mass (denoted by ), which just has a mass property influence on the structural dynamic responses. The frame structure is the primary load bearing structure with high stiffness and high strength; thus, it can also be taken as a seismic mass (denoted by ). Since the viscoelastic connecting parts members are made from a high elasticity material and can be acting as the energy-sapping components that perform the function of vibration suppressing, they are equivalent as an assembly of a spring with certain stiffness and a damper with certain damping. The action effects of multiple viscoelastic connecting parts used for supporting the power module are equal to that of connecting multiple springs and dampers in parallel. In this way, the equivalent model of the spacecraft multifunctional structure can be viewed as a constrained system with damping, as shown in Figure 3.
In Figure 3, and denote the mass of the frame structure and the power module, respectively; subscripts and denote the component of one physical variable in the direction and direction, respectively; and denote the complex stiffness and the damping of the viscoelastic connecting parts, respectively; denotes the structural stiffness of the spacecraft multifunctional structure; denotes the external acceleration excitation, while and denote the distance between the edge of the power module and the frame structure in the and direction, respectively. According to the equivalent vibration reduction model, the frame structure, the power module, and the viscoelastic connecting parts constitute a vibration reduction system, which is the same as a spring-oscillator model, where the viscoelastic connecting parts are equal to elastic supports with certain stiffness and damping, and the power module is equal to an oscillator with certain mass. This vibration reduction system can partly absorb the vibration energy by the damping properties of viscoelastic connecting parts on the one hand and can partly consume the vibration energy by transferring the structural vibration energy to the kinetic energy of the oscillator on the other hand. The result of integration of the two factors above leads to a significant attenuation of vibration transmission in the structure and makes an efficient management for the mechanical environment of the power module.
3. Complex Zero-Point Equation of System
In this section, we discuss the detailed analysis of the equivalent simplified vibration reduction model shown in Figure 3. Taking the axial symmetry of the model into consideration, and meanwhile supposing that the stiffness and damping provided by the spring and damper in the direction can be merged into those in the direction, we can thus obtain a more simplified system model in which the mechanical behaviors in the direction and in the direction are decoupled, as shown in Figure 4.
In Figure 4, and denote the equivalent complex stiffness and damping in the direction, respectively. For this simplified model shown in Figure 4, the equations of motion of the vibration system can be written aswhere is the external acceleration exciting force acting in the direction.
It can be further written in matrix form aswhere is the mass matrix, is the damping matrix, is the stiffness matrix, denotes the input, and denotes the response.
Make the Fourier transformation for (2); we obtainwhere .
In (3), the term of has the property of stiffness, which hinders the movement of the system and can be called dynamic stiffness. Under the same excitation, the larger the dynamic stiffness is, the smaller the response will be. Hence, it is also known as impedance. For simplicity, we omit the subscript in the subsequent analysis: namely, we use , , and to substitute , , and , respectively. Thus, the impedance matrix of the system can be expressed as
The frequency response function of the system is the inverse matrix of the impedance matrix; namely,
It further can be written as
Therefore, the system equation can be expressed as follows:or
In (6), each element in the frequency response function matrix has a specific physical meaning and plays a different role in the system transfer characteristics. Generally, is referred to as the frequency response function of point on the input of point . When l = p, it is also known as the origin frequency response function; when l ≠ p, it is called the cross-point frequency response function.
It can be obtained from (8) thatwhere the origin frequency response function is
The origin frequency response function is a complex number. Its real part and imaginary part can be, respectively, written as below
We have analyzed the amplitude-frequency characteristic of the origin frequency response function , to get its logarithmic amplitude-frequency characteristic curve, which is shown in Figure 5.
The following can be seen from Figure 5:
() When the denominator of is equal to zero, namely, , the system reaches a resonance state.
() When the numerator of is equal to zero, gets the minimum value . It can be known from (9) that the vibration response (i.e., displacement ) of mass in the system becomes a minimum value, so the system is in the antiresonance state. The so-called antiresonance refers to the fact that, under the action of harmonic excitations with some specific frequencies, the elastic system causes the situation that the harmonic response of some parts of system is approximate to zero. In other words, the antiresonance situation means that the dynamic compliance of some parts of the system is approximate to zero at certain frequencies [17–20]. Here, we called the equationsoras the basic complex zero-point equations of two-degree-of-freedom constrained system with damping, where superscripts and denote for the real part and imaginary part, respectively.
It is easy to derive that the conditions for making (14) hold are
Antiresonance is an important parameter reflecting the local characteristics of vibration system (or the dynamic characteristics of excitation points), which is related to the physical parameters of some substructures in the system. The principle of vibration antiresonance design is that the structural parameters in the system should be reasonably designed to ensure that the complex zero-point equation of the vibration system will be satisfied: namely, the transfer function is satisfied to the minimum value condition; thus, the vibration of the concerned position in the system achieves a minimum level.
4. Vibration Antiresonance Design
For the spacecraft multifunctional structure described in this paper, the result that the frame structure will almost generate no vibration under external excitations is desired. According to the principle of antiresonance explained in the above section, we can choose appropriate parameters to make (13) or (14) hold, so that the vibration of the frame structure of the multifunctional structure remains at a minimum level under external exciting force and thus the purpose of vibration reduction is achieved. Supposing that the structural stiffness of the multifunctional structure and the mass of the frame structure are already known preliminarily, the main unknown parameters existing in the complex zero-point equation of the system mainly include the mass of the supported structure , the stiffness of the viscoelastic connecting parts , and damping of the viscoelastic connecting parts . The following context will give a detailed discussion on the determining method of these parameters.
4.1. Stiffness and Damping of the Viscoelastic Connecting Parts
The viscoelastic connecting parts in the spacecraft multifunctional structure are an elastic member made of a rubber-like material, which has both the characteristics of an elastic solid and a viscous liquid under dynamic loadings. The elasticity demonstrates the stiffness characteristic of the member, while the viscous property demonstrates the damping characteristic of the member. Because the viscoelastic connecting parts have the anelasticity phenomenon, its stress-strain relationship is not linear. Therefore, it needs accurate viscoelastic constitutive model to predict its dynamic mechanical response. In this section, the Havriliak-Negami model (referred to as H-N model) is used to model the constitutive properties of the viscoelastic connecting parts in spacecraft multifunctional structures.
By combining the complex modulus of the viscoelastic material , with the low-frequency rubbery state modulus and the high-frequency glass state modulus , the viscoelastic constitutive model based on the H-N model is obtained as below :where is the angular frequency, is the relaxation time related to the temperature, and is the imaginary unit; , , , and are model parameters independent of temperature, and , is the storage modulus of viscoelastic material, is the loss modulus of viscoelastic material. When , ; when , .
The relationship curve between and is commonly known as the Cole-Cole curve, as shown in Figure 6, where denotes the storage modulus axis and denotes the loss modulus axis.
The parameters and in the viscoelastic constitutive model based on H-N model can be determined by an iterative method given in . Thus, for the viscoelastic connecting parts structure with known shape, the computing formula of the stiffness coefficient and the damping coefficient can be given. In the paper, the viscoelastic connecting parts structure shown in Figure 2 is a solid cubic structure with length L, width , and height , respectively; thus, we havewhere is the number of viscoelastic connecting parts at one side of the power module.
It can be seen from (17) that, when the parameters of the viscoelastic constitutive model are determined based on the H-N model, the material storage modulus and loss modulus of viscoelastic connecting parts can be obtained, and thus the stiffness coefficient and damping coefficient of viscoelastic connecting parts can be obtained.
4.2. Mass of the Power Module
According to the analysis result of Section 4.1, we think the stiffness and damping parameters of the viscoelastic connecting parts are already obtained by using the nonlinear mechanical constitutive relation model of the rubber materials. Based on the complex zero-point equation of the vibration system, the mass of the power module must be satisfied to the following equation:where is the vibration angular frequency of the system; we can use the antiresonance frequency need to be concerned in the calculation, which is the frequency under the antiresonance state.
By substituting (17) into (14), we can get the final complex zero-point equations of the spacecraft multifunctional structure with the viscoelastic material properties as follows:where the physical dimensions including length L, width , and height and the number of the viscoelastic connecting parts are the structural parameters need to be determined in the design of the spacecraft multifunctional structure.
5. Verification of the Effectiveness of the Method
In this section, both the numerical simulation and ground experiment are used to verify the effectiveness of the vibration antiresonance design method proposed in this paper.
5.1. Numerical Simulation Analysis
It is assumed that the frame structure of the spacecraft multifunctional structure is made of M40J composite material, whose Equivalent Young’s modulus is 76.88 GPa, Poisson ratio is 0.319, and mass density is 1510 kg/m3. We suppose that the spacecraft multifunctional structure is under the action of a harmonic acceleration excitation. Here we analyze the structural response under harmonic excitation in the condition of . The results of acceleration response of the power module and the frame structure in the direction are shown in Figure 7, respectively.
We can know from Figure 7 that the response amplitude of the power module is less than that of the exciting loads, and the larger the frequency is, the smaller the response amplitude of the power module will be. Moreover, there is an obvious phase difference between the response amplitude of the power module and the exciting loads, and the higher the exciting frequency, the more obvious the nonlinearity of the acceleration response of the power module.
Make a further analysis on the energy dissipation of the spacecraft multifunctional structure under harmonic excitation. Let the harmonic acceleration excitation be , where (1 g = 9800 mm/s2) and . We conduct a numerical simulation of the viscous dissipation energy and the dissipation rate on the conditions of Hz, Hz, Hz, and Hz, and the time duration used in the simulation is set to be 0.02 s. The viscous dissipation energy and the dissipation rate of the spacecraft multifunctional structure under harmonic excitation in the x direction are shown in Figures 8(a) and 8(b), respectively.
(a) The change of viscous dissipative energy with time
(b) The change of dissipation rate with time
We can see from Figure 8 that the energy dissipation rate on the condition of (close to the resonance frequency of system) is much larger than that on the conditions of Hz, Hz, and Hz. When the exciting frequency is Hz, the viscous dissipation energy increases monotonically with time and reaches a maximum of 0.0001238 J (in the x direction) and 0.0001376 J (in the y direction) at the time of 0.02 s. As to the viscous dissipation rate, it tends to change periodically with the increasing of time.
5.2. Ground Experimental Verification
The satellite is located in the fairing of rocket and will be subjected to the vibration load transmitted by the transition section during rocket launch. In order to analyze the dynamic characteristics of the spacecraft multifunctional structure under sinusoidal vibration loads, we conduct a sinusoidal vibration experiment research. The sinusoidal vibration experimental conditions are shown in Table 1, and the sinusoidal vibration experimental photo is shown in Figure 9.
|: “” denotes single peak value.|
In the dynamic test of the spacecraft multifunctional structure, the distribution of the acceleration pick-up sensors A, B, C, and D is shown in Figure 10. The type of the acceleration pick-up sensors used in the experiment is BK4517 (Denmark Bruel & Kjaer Ltd.), their physical sizes are 11.48 mm × 6.35 mm × 3.8 mm, the sensitivity is 10 mV/g, the frequency range is from 1 Hz to 20000 Hz, the resonance frequency is 55 kHz, the working temperature range is from −55°C to 121°C, and the mass is 0.6 × 10−3 kg.
The vibration response of the spacecraft multifunctional structure in the direction as measured by sensors A and B is shown in Figure 11.
From Figure 11, we can see that the acceleration output amplitude of measuring point B is obviously lower than that of measuring point A. Hence, we can conclude that the external vibration is attenuated significantly through the transmission of the spacecraft multifunctional structure.
For the new spacecraft multifunctional structure design concept mainly composed of the frame structure, the power modules, and the viscoelastic connecting parts, a vibration antiresonance design method is proposed in this paper to reduce the vibration transmission through the structure. The basic complex zero-point equation of the vibration system is established. As to the special difficulties caused by the viscoelastic characteristics of the connecting parts, we deduce the determining formulas to obtain the structural parameters, such as the stiffness and damping of the viscoelastic connecting parts, which will make the complex zero-point equations hold. Furthermore, we get the final complex zero-point equations of the spacecraft multifunctional structure, which has taken the viscoelastic material properties into account. Therefore, the vibration of the frame structure decreases under the action of external exciting forces and thus the wanted result of structural vibration reduction is achieved. Numerical simulation and ground experiment results show that the vibration reduction function of such spacecraft multifunctional structure is achieved successfully by using this vibration antiresonance design. The method proposed in this paper solves the structural vibration control problem under the function constraints of load bearing and energy supplying and expands the performance of spacecraft functional modules.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
- P. Fortescue, G. Swinerd, and J. Stark, Spacecraft Systems Engineering, John Wiley & Sons, 2011.
- T. Dunham, Aerospace Engineering: Design, Concepts and Applications, Clanrye International, 2016.
- S. C. Roberts and G. S. Aglietti, “Structural performance of a multifunctional spacecraft structure based on plastic lithium-ion batteries,” Acta Astronautica, vol. 67, no. 3-4, pp. 424–439, 2010.
- E. Harris N and R. Morgenthaler D, “Design & testing of multifunctional structure concept for spacecraft,” in Proceedings of the 41st AIAA/ASME/ASCE/AHS/ASC Structure, Structure Dynamics, and Materials Conference and Exhibit, Atlanta, Ga, USA, 2000.
- D. M. Barnett, S. Rawal, and K. Rummel, “Multifunctional structures for advanced spacecraft,” Journal of Spacecraft and Rockets, vol. 38, no. 2, pp. 226–230, 2001.
- S. E. Hahn, R. Usami, and T. Ozaki, “Multifunctional structure spacecraft bus technology,” in Proceedings of the A Collection of the 22nd AIAA International Communications Satellite Systems Conference and Exhibit, pp. 226–239, Monterey, Calif, USA, May 2004.
- Z. Xie, W. S. Shepard, and K. A. Woodbury, “Design optimization for vibration reduction of viscoelastic damped structures using genetic algorithms,” Shock and Vibration, vol. 16, no. 5, pp. 455–466, 2009.
- X. P. Do, K. Shah, and S.-B. Choi, “Damping force tracking control of MR damper system using a new direct adaptive fuzzy controller,” Shock and Vibration, vol. 2015, Article ID 947937, 2015.
- Z. Lu, D. Wang, and P. Li, “Comparison study of vibration control effects between suspended tuned mass damper and particle damper,” Shock and Vibration, vol. 2014, Article ID 903780, 2014.
- Y. Du, S. Wang, Y. Zhu, L. Li, and G. Han, “Performance of a new fine particle impact damper,” Advances in Acoustics and Vibration, vol. 2008, pp. 1–6, 2008.
- P. Y. Bely, K. J. Pendergast, C. J. Schauwecker, and J. B. Breckinridge, “Use of a passive reaction wheel jitter isolation system to meet the Advanced X-Ray Astrophysics Facility imaging performance requirements,” in Proceedings of the Astronomical Telescopes & Instrumentation, p. 1078, Kona, Hawaii, USA.
- A. J. Bronowicki, “Vibration isolator for large space telescopes,” Journal of Spacecraft and Rockets, vol. 43, no. 1, pp. 45–53, 2006.
- V. Babuška, R. S. Erwin, and L. A. Sullivan, “System identification of the suite isolation platform: Comparison of ground and flight experiments,” in Proceedings of the 44th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 2003, Norfolk, Va, USA, April 2003.
- V. L. Shinde and A. K. Pathak, “Review on particle damping technique for vibration suppression [J],” International Journal of Innovative Research in Science, Engineering and Technology, vol. 5, no. 3, pp. 2890–2895, 2016.
- V. U. Sakhare, A. S. Aradhye, and B. S. Gandhare, “A review on design and validation of performance of particle impact damper,” International Journal of Current Engineering and Technology, vol. 6, no. 6, pp. 2038–2041, 2016.
- P. S. Kachare and Bimleshkumar, “Effective location and packing ratio of a particle impact damper for a mild steel cantilever beam under transient and forced vibration,” International Review of Mechanical Engineering, vol. 5, no. 7, pp. 1226–1234, 2011.
- E. Kausel, Advanced Structural Dynamics, Cambridge University Press, 2017.
- M. Berenguel and E. F. Camacho, “Frequency-based adaptive control of systems with antiresonance modes,” Control Engineering Practice, vol. 4, no. 5, pp. 677–684, 1996.
- K. Jones and J. Turcotte, “Finite element model updating using antiresonant frequencies,” Journal of Sound and Vibration, vol. 252, no. 4, pp. 717–727, 2002.
- V. Arora, “Use of resonance and antiresonance frequencies for better matching of frequency response function,” International Journal of Structural Engineering, vol. 5, no. 1, pp. 13–23, 2014.
- D. Hao and D. Li, “Determination of dynamic mechanical properties of carbon black filled rubbers at wide frequency range using Havriliak-Negami model,” European Journal of Mechanics, A/Solids, vol. 53, pp. 303–310, 2015.
Copyright © 2017 Dong-Xu Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.