Feasibility Study of the Electromagnetic Damper for Cable Structures Using Real-Time Hybrid Simulation

Cable structure is a major component of long-span bridges, such as cable-stayed and suspension bridges, and it transfers the main loads of bridges to the pylons. As these cable structures are exposed to continuous external loads, such as vehicle and wind loads, vibration control and continuous monitoring of the cable are required. In this study, an electromagnetic (EM) damper was designed and fabricated for vibration control and monitoring of the cable structure. EM dampers, also called regenerative dampers, consist of permanent magnets and coils. The electromagnetic force due to the relative motion between the coil and the permanent magnet can be used to control the vibration of the structure. The electrical energy can be used as a power source for the monitoring system. The effects of the design parameters of the damper were numerically analyzed and the damper was fabricated. The characteristics of the damper were analyzed with various external load changes. Finally, the vibration-control and energy-harvesting performances of the cable structure were evaluated through a hybrid simulation. The vibration-control and energy-harvesting performances for various loads were analyzed and the applicability to the cable structure of the EM damper was evaluated.


Introduction
In Korea, structural accidents, such as the collapses of the Gyeongju Maury Resort (2014) and the Busan North-South Bridge (2013), occur frequently. Recently, the cable fire accident of the Seohae Grand Bridge (2015) caused a casualty, and consumed considerable time and cost to recover the bridge. Annually, various structural failures cause casualties and financial losses. In order to prevent the structural failure, vibration control and monitoring systems are applied to various civil structures. The bridge is a representative civil engineering structure that not only plays a role in contributing to economic growth by utilizing land more efficiently, but also enables the balanced development of national industries and regions. With the continuous development of construction technology, the rate of the construction of cable bridges, such as cable-stayed and suspension bridges, is increasing. It is necessary to continuously maintain and monitor a bridge cable because the cable, which is a main component of the cable bridge, transfers various loads to the main tower by tension of the cable. Due to the characteristics of the cable structures, the bridge cables have a very low damping ratio. These low damping ratios cause fatigue loads on the cables, by continuous exposure to wind and vehicle loads. The fatigue loads acting on the cable lead to the deterioration of the lifetime of the bridge. Therefore, vibration control and monitoring of the cable are required. The electromagnetic (EM) damper, which has both vibration-control and energy-harvesting functions, has been proposed as an effective maintenance method for the cable structure. The EM damper can reduce the vibration This study aims to evaluate the feasibility of the EM damper for vibration-control and structure monitoring. In this study, a tubular-type linear damper was first fabricated as a vibration-control method of a cable structure, and characteristic and performance tests were conducted. The proposed EM damper has capabilities of vibration control and energy harvesting, which can be used as a source of the energy for the monitoring system. The characteristics of the damper were analyzed through characteristic tests. For the various external loads, the variation of the damper force was measured. Before applying the damper to the target structure, the behavior of the damper according to the external load was investigated. A real-time hybrid simulation was performed to evaluate the applicability of the damper for the cable structure. Contrary to previous studies, the performance of the EM damper was analyzed using a hybrid simulation incorporating the bridge cable properties to reflect the actual characteristics of the cable structure. The vibration-control and energy-harvesting performances of the cable structure were evaluated through a hybrid simulation that reflected the characteristics of the actual damper, rather than a numerical analysis reflecting the ideal damper model. In order to analyze the energy-harvesting performance of the EM damper, the mean wind speed acting on the actual bridge was considered.

Electromagnetic (EM) Damper Design
An EM damper is a typical regenerative damping system and has been widely used for automobile vibration control. The regenerative damping utilizes the electromagnetic force generated during the process of converting mechanical vibration energy of a structure into usable electrical energy as system damping. The regenerative damping is used in various fields as the utilization of regenerative energy generated in the vibration-control process increases. Recently, studies are being conducted to use EM dampers for vibration control of civil engineering structures.
The ideal model of the EM damper force can be expressed by the damper coefficient and the velocity of the damper. For the ideal damping force and damping coefficient, the time delay and mechanical friction of the damper are not considered. The damping force of the EM damper is proportional to the speed of the structure, like in the conventional mechanical damper. The electromagnetic damping force F d is proportional to the EM damper coefficient c d and applied velocity v: The damping coefficient of the damper must include a time delay due to the inductance of the damper coil. The damping force-velocity relation including the time delay component of the damper can be expressed as follows: where L coil is the inductance of the damper coil, R coil is the resistance of the damper coil, and K t is a motor constant of the damper. Table 1 and Figure 1 show the design parameters of the linear-type EM damper. The parameters in the table include the motor constant of the damper. The detailed derivation of the motor constant K t is given in Palomera-Arias (2005).
Using Equation (2), the damper coefficient, including the time delay, is as follows: where ω is a frequency of the input excitation. The inductance is influenced by the excitation frequency of the external load. The influence of the inductance is considerable when the external load has a high-frequency band. However, considering that the natural frequencies of most civil infrastructures are in a low-frequency band, the influence of the inductance on the damping coefficient of the EM damper is negligible.    The damping coefficient of the EM damper is influenced by the motor constant, system resistance, and inductance. The maximum damping coefficient of the EM damper, when neglecting the influence of inductance at low frequencies, is as follows: The resistance R total of the EM damper is the sum of the resistance of the solenoid coil and the resistance of the additional circuitry attached to the outside of the EM damper. Considering c d,max , only, the resistance change of the solenoid coil according to the variation of the coil wire thickness is obtained without including the external resistance. The number of turns, wire length, and coil resistance, which are the main properties of a coil, can be estimated by the following equations [17,18]: where V coil is the volume of the coil, V w is the volume of the coil wire, d w is the diameter of the coil wire, l w is the total length of the coil wire, R coil is the resistance of the coil, and ρ is a conductivity of the material of the coil. Figure 2 shows the damping coefficient and density variation of the EM damper. The damping density is obtained by considering both the damping coefficient and volume of the damper. The damping density is the value of the damping coefficient divided by the total volume of the damper, which represents the efficiency according to the total volume of the damper. The damping coefficient and density variation of the magnet radius, air gap, wire radius, stator yoke thickness, magnet length, and pole pitch, were analyzed, all of which have a major influence on the damping coefficient and damping density of the damper. The observed damping coefficient showed a tendency to increase as the air gap and pole pitch decreased, and as magnet radius, magnet length, and stator yoke thickness increased.
Unlike the damping coefficient that continuously increased, the damping density tended to converge to a constant value. For the magnet length and radius, the damping density tended not to increase significantly when a certain value was reached, and the maximum damping density was 0.02 m for the stator yoke thickness. In the case of the coil, the increase of the radius caused a decrease in the resistance per unit length, but a decrease in the total number of turns and the total length. Conversely, if the radius decreased, the magnitude of the resistance per unit length increased, but the total number of turns and length both increased. By considering all such trends, the coil radius change did not significantly affect the damping coefficient and density. Based on the variation of the damping density according to the design variables, the design parameters shown in Table 2 were selected and the damper was fabricated.

Characteristic Test of the EM Damper
The characteristic test of the damper was carried out on the change of the damping force according to the excitation frequency and displacement change. The excitation frequency was changed from 1 to 5 Hz in 1 Hz increments, and the displacement was changed from 6 to 15 mm in 3 mm increments. Table 3 and Figure 3 show the characteristic test cases and damper design with its configuration, respectively. Each column in Table 3 represents the amplitude change of the external load at the fixed frequency. The experimental cases shown in the table mean the amplitude change for the constant frequency of the external load.

Characteristic Test of the EM Damper
The characteristic test of the damper was carried out on the change of the damping force according to the excitation frequency and displacement change. The excitation frequency was changed from 1 to 5 Hz in 1 Hz increments, and the displacement was changed from 6 to 15 mm in 3 mm increments. Table 3 and Figure 3 show the characteristic test cases and damper design with its configuration, respectively. Each column in Table 3 represents the amplitude change of the external load at the fixed frequency. The experimental cases shown in the table mean the amplitude change for the constant frequency of the external load.    Figure 4 shows the damping force displacement and damping force velocity with varying excitation amplitudes. As shown in the figures, the damping force of the damper was linearly proportional to the amplitude of the applied load for a constant frequency input. The hysteresis curves of the damper showed an elliptical shape. The direction of the major axis of the ellipse was kept constant and the width of the ellipse was increased when an external load of various amplitudes was applied with a constant frequency.
The force-displacement curve showed a negative stiffness, induced by the magnetic force of the damper. A permanent magnet, which was placed inside a coil, induced an induction current in the coil, producing the Lorentz force. The Lorentz force generated in the damper acted as a damping force of the damper. The Lorentz force acted in the opposite direction to the permanent magnet. Therefore, the force-displacement hysteresis curve of the EM damper showed the negative stiffness.   Figure 4 shows the damping force displacement and damping force velocity with varying excitation amplitudes. As shown in the figures, the damping force of the damper was linearly proportional to the amplitude of the applied load for a constant frequency input. The hysteresis curves of the damper showed an elliptical shape. The direction of the major axis of the ellipse was kept constant and the width of the ellipse was increased when an external load of various amplitudes was applied with a constant frequency.
The force-displacement curve showed a negative stiffness, induced by the magnetic force of the damper. A permanent magnet, which was placed inside a coil, induced an induction current in the coil, producing the Lorentz force. The Lorentz force generated in the damper acted as a damping force of the damper. The Lorentz force acted in the opposite direction to the permanent magnet. Therefore, the force-displacement hysteresis curve of the EM damper showed the negative stiffness.
The variation of the damping force as a function of the frequency of the EM damper is shown in Figure 5. For the variations in the external frequency with constant amplitude, the damping force of the damper increased nonlinearly. The major axis of the ellipse of the hysteresis curve shifted when the excitation frequency changed for a constant amplitude. In the displacement-damper force graph, the lengths of the minor and major axes increased simultaneously. That is, as the frequency of the applied load increased, the hysteresis curve gradually became circular. As the excitation frequency increased, the slope of the hysteresis curve increased as the damping coefficient of the damper increased. In addition, as the excitation frequency increased, the area of the hysteresis curve of the displacement-damping force increased. This indicated that the amount of total energy that the damper can dissipate increased.  The variation of the damping force as a function of the frequency of the EM damper is shown in Figure 5. For the variations in the external frequency with constant amplitude, the damping force of the damper increased nonlinearly. The major axis of the ellipse of the hysteresis curve shifted when the excitation frequency changed for a constant amplitude. In the displacement-damper force graph, the lengths of the minor and major axes increased simultaneously. That is, as the frequency of the applied load increased, the hysteresis curve gradually became circular. As the excitation frequency increased, the slope of the hysteresis curve increased as the damping coefficient of the damper increased. In addition, as the excitation frequency increased, the area of the hysteresis curve of the displacement-damping force increased. This indicated that the amount of total energy that the damper can dissipate increased.  Figure 6a shows the maximum damping force of the damper for changes in the amplitude of the external load. The maximum damping force of the damper linearly increased as the displacement of the input load increased. The damper force increased linearly as the amplitude increased, even when the input load frequency was changed. Figure 6b shows the variations in the maximum damping force and its nonlinear changes for a varied excitation frequency. In the velocity-force hysteresis curve, the damping coefficient increased as the excitation frequency increased. The nonlinear changes   The variation of the damping force as a function of the frequency of the EM damper is shown in Figure 5. For the variations in the external frequency with constant amplitude, the damping force of the damper increased nonlinearly. The major axis of the ellipse of the hysteresis curve shifted when the excitation frequency changed for a constant amplitude. In the displacement-damper force graph, the lengths of the minor and major axes increased simultaneously. That is, as the frequency of the applied load increased, the hysteresis curve gradually became circular. As the excitation frequency increased, the slope of the hysteresis curve increased as the damping coefficient of the damper increased. In addition, as the excitation frequency increased, the area of the hysteresis curve of the displacement-damping force increased. This indicated that the amount of total energy that the damper can dissipate increased.  Figure 6a shows the maximum damping force of the damper for changes in the amplitude of the external load. The maximum damping force of the damper linearly increased as the displacement of the input load increased. The damper force increased linearly as the amplitude increased, even when the input load frequency was changed. Figure 6b shows the variations in the maximum damping force and its nonlinear changes for a varied excitation frequency. In the velocity-force hysteresis curve, the damping coefficient increased as the excitation frequency increased. The nonlinear changes of the hysteresis curve occurred at the input frequency change. For a constant input amplitude, the maximum damping force of the damper increased nonlinearly as the excitation frequency increased. As the magnitude of the change in the magnet flux acting on the coil increased, the damping force of the damper nonlinearly increased as the excitation frequency increased. When the excitation frequency exceeded 5 Hz, the change in the magnet flux reached saturation. The maximum damping force of the damper converged after 5 Hz.   Figure 6a shows the maximum damping force of the damper for changes in the amplitude of the external load. The maximum damping force of the damper linearly increased as the displacement of the input load increased. The damper force increased linearly as the amplitude increased, even when the input load frequency was changed. Figure 6b shows the variations in the maximum damping force and its nonlinear changes for a varied excitation frequency. In the velocity-force hysteresis curve, the damping coefficient increased as the excitation frequency increased. The nonlinear changes of the hysteresis curve occurred at the input frequency change. For a constant input amplitude, the maximum damping force of the damper increased nonlinearly as the excitation frequency increased. As the magnitude of the change in the magnet flux acting on the coil increased, the damping force of the damper nonlinearly increased as the excitation frequency increased. When the excitation frequency

Hybrid Simulation of the EM Damper
The validation of the vibration-control and energy-harvesting performance of the EM damper was performed using a hybrid simulation. The hybrid simulation is a combination of the experimental method of the physical experiment and a numerical simulation. By combining the numerical simulation and experiment, hybrid simulations can be used to verify full-scale structures or virtual structures that are difficult to verify experimentally [19][20][21]. In this study, a hybrid simulation was constructed using a numerical cable structure and a real damper model. As mentioned previously, the numerical model of the EM damper did not reflect the characteristics of a real EM damper. In this study, the hybrid simulation was constructed to experimentally reflect the behavior of a real EM damper.
In the process of the hybrid simulation, the measurement of the physical experiment and the calculation of the numerical simulation are conducted simultaneously. For the external loads, the response of the numerical cable structure is calculated in the numerical simulation. At the same time, the cable responses at the EM damper position are calculated and applied to the EM damper. Then, the damping force of the EM damper is measured and input to the numerical simulation simultaneously. All processes of calculation and measurement are performed in real-time.
Since the hybrid simulation had few constraints on the target structure, it is possible to evaluate the performance of the proposed damper by applying the damper to various structures. The physical experiment was performed using a shaking table and the EM damper. The numerical model of the cable was built on the computer. Using this hybrid simulation setup, the response of the cable structure to various external loads, the vibration-control performance of the damper, and the energyharvesting performance were evaluated. Figure 7 shows a block diagram of the hybrid simulation.
The hybrid simulation was constructed using Matlab Simulink Desktop Real-Time. In the hybrid simulation process, the numerical analysis was conducted by iterating the response of the cable at the damper location and the corresponding damper force. If the computing speed of the simulation could not reach the simulation time, a time delay occurred. However, the time delay was prevented by setting missed ticks. The hybrid simulation was stopped immediately when the number of delayed time steps exceeded the number of missed ticks. In the hybrid simulation, the number of missed ticks was ~100, which is equivalent to 0.1 s. The maximum time delay that could occur in the hybrid simulation was ~0.1 s. Therefore, the time delay due to the computing process of the hybrid simulation was negligible.

Hybrid Simulation of the EM Damper
The validation of the vibration-control and energy-harvesting performance of the EM damper was performed using a hybrid simulation. The hybrid simulation is a combination of the experimental method of the physical experiment and a numerical simulation. By combining the numerical simulation and experiment, hybrid simulations can be used to verify full-scale structures or virtual structures that are difficult to verify experimentally [19][20][21]. In this study, a hybrid simulation was constructed using a numerical cable structure and a real damper model. As mentioned previously, the numerical model of the EM damper did not reflect the characteristics of a real EM damper. In this study, the hybrid simulation was constructed to experimentally reflect the behavior of a real EM damper.
In the process of the hybrid simulation, the measurement of the physical experiment and the calculation of the numerical simulation are conducted simultaneously. For the external loads, the response of the numerical cable structure is calculated in the numerical simulation. At the same time, the cable responses at the EM damper position are calculated and applied to the EM damper. Then, the damping force of the EM damper is measured and input to the numerical simulation simultaneously. All processes of calculation and measurement are performed in real-time.
Since the hybrid simulation had few constraints on the target structure, it is possible to evaluate the performance of the proposed damper by applying the damper to various structures. The physical experiment was performed using a shaking table and the EM damper. The numerical model of the cable was built on the computer. Using this hybrid simulation setup, the response of the cable structure to various external loads, the vibration-control performance of the damper, and the energy-harvesting performance were evaluated. Figure 7 shows a block diagram of the hybrid simulation.
The hybrid simulation was constructed using Matlab Simulink Desktop Real-Time. In the hybrid simulation process, the numerical analysis was conducted by iterating the response of the cable at the damper location and the corresponding damper force. If the computing speed of the simulation could not reach the simulation time, a time delay occurred. However, the time delay was prevented by setting missed ticks. The hybrid simulation was stopped immediately when the number of delayed time steps exceeded the number of missed ticks. In the hybrid simulation, the number of missed ticks was~100, which is equivalent to 0.1 s. The maximum time delay that could occur in the hybrid simulation was~0.1 s. Therefore, the time delay due to the computing process of the hybrid simulation was negligible.
In order to construct the hybrid simulation, a numerical cable model is required. Figure 8 shows the cable structure and the loads acting on the cable. A deformation of the cable occurred, not only due to the external loads but also due to the sag from self-weight. For the cable with a total length of l, the deflection due to self-weight at an arbitrary position x is given by where m is the mass per unit length, g is the acceleration of gravity, l is the cable length, and T 0 is the tension force acting on the cable. In order to construct the hybrid simulation, a numerical cable model is required. Figure 8 shows the cable structure and the loads acting on the cable. A deformation of the cable occurred, not only due to the external loads but also due to the sag from self-weight. For the cable with a total length of l, the deflection due to self-weight at an arbitrary position x is given by where m is the mass per unit length, g is the acceleration of gravity, l is the cable length, and 0 is the tension force acting on the cable. When an external force is applied to the cable, a dynamic displacement occurs, causing In order to construct the hybrid simulation, a numerical cable model is required. Figure 8 shows the cable structure and the loads acting on the cable. A deformation of the cable occurred, not only due to the external loads but also due to the sag from self-weight. For the cable with a total length of l, the deflection due to self-weight at an arbitrary position x is given by where m is the mass per unit length, g is the acceleration of gravity, l is the cable length, and 0 is the tension force acting on the cable. When an external force is applied to the cable, a dynamic displacement occurs, causing additional tension in the cable. If the additional tension caused by the dynamic displacement is T, the equation of motion of the cable is given by  When an external force is applied to the cable, a dynamic displacement occurs, causing additional tension in the cable. If the additional tension caused by the dynamic displacement is T, the equation of motion of the cable is given by The occurrence of instantaneous deformation in a cable is due to the external force and dynamic displacement. This deformation causes additional tension inside the cable. The equation of motion, which takes into account both the dynamic displacement of the cable and the additional strain due to deformation, is expressed by where E is the young's modulus of the cable and A is the area of the cable. The transverse deflection of the cable can be approximated by the sum of the sine-shaped assumed shape functions: where n is the number of modes considered, q i (t) are the generalized displacements, and φ i (t) is a set of the shape functions. According to Johnson et al. (2000), introducing shape functions based on the deflection due to a static force at the damper can reduce the number of terms required for comparable accuracy. Conversely, hundreds of terms in the sine series are generally used as shape functions if this static deflection is not considered, resulting in considerable computation effort [22,23]. The shape function, including the deflection due to a static force, can be expressed by Using Equations (12) and (14), a non-dimensional equation of motion of the cable and its matrix form can be expressed as [24]: where M = m ij = m l 0 φ i φ j dx is the mass matrix, C = αM + βK is the damping matrix calculated form Rayleigh damping, and is the stiffness matrix. The properties of the cable damper model are shown in Table 4. The free-vibration response at the center point of the cable according to the change in the external resistance is shown in Figure 9. As the larger external resistance was connected, the smaller control force of the damper was shown. In the free-vibration test, the vibration dissipation capacity of the damper sharply decreased as the external resistance increased. In the uncontrolled case without the damper, the cable damping ratio was 0.48%. For the short circuit case with zero external resistance, the damping ratio was increased by~3 to 1.32%. As the external resistance increased, the damping ratio gradually decreased to~1%. Table 5 shows the damping ratio change of the cable due to the variations in the external resistance. The free-vibration response at the center point of the cable according to the change in the external resistance is shown in Figure 9. As the larger external resistance was connected, the smaller control force of the damper was shown. In the free-vibration test, the vibration dissipation capacity of the damper sharply decreased as the external resistance increased. In the uncontrolled case without the damper, the cable damping ratio was 0.48%. For the short circuit case with zero external resistance, the damping ratio was increased by ~3 to 1.32%. As the external resistance increased, the damping ratio gradually decreased to ~1%. Table 5 shows the damping ratio change of the cable due to the variations in the external resistance.   11 show the displacement and acceleration response of the cable with the natural frequency excitation, respectively. As shown in the figures, for a constant sinusoidal input load, the damper showed a certain level of damping performance. For each sinusoidal input load, the magnitude of the maximum response of cable was reduced to less than half in both the displacement and acceleration responses. Figure 12 shows the frequency domain response for each sinusoidal input. The magnitude of the response was reduced to less than half, even in the frequency domain with respect to each natural frequency. At least 50% of the response attenuation was shown for all the frequency input. The results of the sinusoidal input are shown in Table 6. The uncontrolled, controlled, and decay ratio of displacement, acceleration, and FFT amplitude are shown in Table 6.     Figures 10 and 11 show the displacement and acceleration response of the cable with the natural frequency excitation, respectively. As shown in the figures, for a constant sinusoidal input load, the damper showed a certain level of damping performance. For each sinusoidal input load, the magnitude of the maximum response of cable was reduced to less than half in both the displacement and acceleration responses. Figure 12 shows the frequency domain response for each sinusoidal input. The magnitude of the response was reduced to less than half, even in the frequency domain with respect to each natural frequency. At least 50% of the response attenuation was shown for all the frequency input. The results of the sinusoidal input are shown in Table 6. The uncontrolled, controlled, and decay ratio of displacement, acceleration, and FFT amplitude are shown in Table 6. Figure 13 shows the displacement and acceleration response of the cable as the resistance changed. Similar to the free-vibration test, the damper showed better damping performance as the external resistance decreased. For the short circuit case, the response decreased by more than 40% for all values of mean wind speed. Figure 14 shows the change in RMS value of the cable displacement for the variations in the mean wind speed. It represents the ratio of the RMS response for the uncontrolled case without damper. The RMS response was reduced by more than 40% for all values of mean wind speed, and the RMS response reduction was a maximum at a mean wind speed of 3 m/s. Figure 15 shows the hysteresis curve for a mean wind speed of 5 m/s. Unlike the characteristic test, the wind load acting on the cable contained various frequency components. As the cable in which the wind load was applied was excited by various frequency components, the hysteresis curve of the damper also showed characteristics including various frequency components. Both the displacement-force and velocity-force graphs did not exhibit a constant ellipse trajectory like the sinusoidal input. As the various frequency components were contained, several ellipsoidal trajectories are shown in the hysteresis curve.
Controlled 54 ohm 0.0100 Figures 10 and 11 show the displacement and acceleration response of the cable with the natural frequency excitation, respectively. As shown in the figures, for a constant sinusoidal input load, the damper showed a certain level of damping performance. For each sinusoidal input load, the magnitude of the maximum response of cable was reduced to less than half in both the displacement and acceleration responses. Figure 12 shows the frequency domain response for each sinusoidal input. The magnitude of the response was reduced to less than half, even in the frequency domain with respect to each natural frequency. At least 50% of the response attenuation was shown for all the frequency input. The results of the sinusoidal input are shown in Table 6. The uncontrolled, controlled, and decay ratio of displacement, acceleration, and FFT amplitude are shown in Table 6.        Figure 13 shows the displacement and acceleration response of the cable as the resistance changed. Similar to the free-vibration test, the damper showed better damping performance as the external resistance decreased. For the short circuit case, the response decreased by more than 40% for all values of mean wind speed. Figure 14 shows the change in RMS value of the cable displacement for the variations in the mean wind speed. It represents the ratio of the RMS response for the uncontrolled case without damper. The RMS response was reduced by more than 40% for all values of mean wind speed, and the RMS response reduction was a maximum at a mean wind speed of 3 m/s.      Figure 13 shows the displacement and acceleration response of the cable as the resistance changed. Similar to the free-vibration test, the damper showed better damping performance as the external resistance decreased. For the short circuit case, the response decreased by more than 40% for all values of mean wind speed. Figure 14 shows the change in RMS value of the cable displacement for the variations in the mean wind speed. It represents the ratio of the RMS response for the uncontrolled case without damper. The RMS response was reduced by more than 40% for all values of mean wind speed, and the RMS response reduction was a maximum at a mean wind speed of 3 m/s.     The voltage and current output of the EM damper according to the external resistance change for a 5 m/s mean wind speed are shown in Figure 16 and Table 7. For the case of zero external resistance, the output current from the damper was a maximum. However, the power output from the damper was a minimum because of the relatively small voltage output. As the magnitude of the external resistor increased, the magnitude of the current generated in the damper gradually decreased. However, the voltage generated by the damper was relatively increased. The magnitude   The voltage and current output of the EM damper according to the external resistance change for a 5 m/s mean wind speed are shown in Figure 16 and Table 7. For the case of zero external resistance, the output current from the damper was a maximum. However, the power output from the damper was a minimum because of the relatively small voltage output. As the magnitude of the external resistor increased, the magnitude of the current generated in the damper gradually decreased. However, the voltage generated by the damper was relatively increased. The magnitude   The voltage and current output of the EM damper according to the external resistance change for a 5 m/s mean wind speed are shown in Figure 16 and Table 7. For the case of zero external resistance, the output current from the damper was a maximum. However, the power output from the damper was a minimum because of the relatively small voltage output. As the magnitude of the external resistor increased, the magnitude of the current generated in the damper gradually decreased. However, the voltage generated by the damper was relatively increased. The magnitude The voltage and current output of the EM damper according to the external resistance change for a 5 m/s mean wind speed are shown in Figure 16 and Table 7. For the case of zero external resistance, the output current from the damper was a maximum. However, the power output from the damper was a minimum because of the relatively small voltage output. As the magnitude of the external resistor increased, the magnitude of the current generated in the damper gradually decreased. However, the voltage generated by the damper was relatively increased. The magnitude of the power generated by the damper gradually increased as the magnitude of the external resistance increased.
The magnitude of the power generated by the EM damper was related to the damper speed and the magnitude of the external load. The damping coefficient of the damper decreased as the external resistance increased. For a constant magnitude of the external loads, the larger the external resistance applied to the damper, the greater the magnitude of the relative velocity of the permanent magnet and solenoid coil of the damper. As a result, the magnitude of the voltage that could be obtained from the damper increased, but the magnitude of the current gradually decreased. The magnitude of the total power gradually increased, because the magnitude of the voltage rise was relatively larger than the amount of the current reduction, and it finally converged to a specific value. of the power generated by the damper gradually increased as the magnitude of the external resistance increased.
The magnitude of the power generated by the EM damper was related to the damper speed and the magnitude of the external load. The damping coefficient of the damper decreased as the external resistance increased. For a constant magnitude of the external loads, the larger the external resistance applied to the damper, the greater the magnitude of the relative velocity of the permanent magnet and solenoid coil of the damper. As a result, the magnitude of the voltage that could be obtained from the damper increased, but the magnitude of the current gradually decreased. The magnitude of the total power gradually increased, because the magnitude of the voltage rise was relatively larger than the amount of the current reduction, and it finally converged to a specific value.   According to the Korea Meteorological Administration (KMA), the mean wind speed in the past 5 years from 2010 was~5.3 m/s at the 2nd Jindo bridge. The proposed EM damper can acquire power of more than 100 mW when using an external resistance of 4 ohm or more, at a mean wind speed of 5 m/s. According to the Park et al. (2010),~42 mW of power is required when two measurements are performed on the Imote2 wireless sensor [25]. The proposed system shows a sufficient energy-harvesting performance to drive the wireless monitoring sensor, even considering the power loss of the energy-harvesting circuit. However, as shown in Figure 17, if the size of the external resistance is increased to improve the energy-harvesting performance, the vibration-control performance of the cable is remarkably decreased. There is a limit to simultaneously satisfy the energy-harvesting and vibration-control performances. Therefore, it is necessary to develop a proper system operation plan according to changes in the external loads. The proposed EM damper is operated as an energy-harvesting device for the ambient wind loads. Considering the mean wind speed acting on the bridge, a sufficient power can be acquired to drive the monitoring sensors. In addition, the battery cell can be used for storing the remaining power, which can be used when the operation of the energy-harvesting system is difficult. The EM damper is used as the vibration-control system for external loads that cause excessive vibration of cables, such as typhoons and earthquakes. If excessive vibration occurs in the cable, the system is operated in the direction of maximizing the vibration-control performance of the EM damper by minimizing the energy-harvesting performance.  According to the Korea Meteorological Administration (KMA), the mean wind speed in the past 5 years from 2010 was ~5.3 m/s at the 2nd Jindo bridge. The proposed EM damper can acquire power of more than 100 mW when using an external resistance of 4 ohm or more, at a mean wind speed of 5 m/s. According to the Park et al. (2010), ~42 mW of power is required when two measurements are performed on the Imote2 wireless sensor [25]. The proposed system shows a sufficient energyharvesting performance to drive the wireless monitoring sensor, even considering the power loss of the energy-harvesting circuit. However, as shown in Figure 17, if the size of the external resistance is increased to improve the energy-harvesting performance, the vibration-control performance of the cable is remarkably decreased. There is a limit to simultaneously satisfy the energy-harvesting and vibration-control performances. Therefore, it is necessary to develop a proper system operation plan according to changes in the external loads. The proposed EM damper is operated as an energyharvesting device for the ambient wind loads. Considering the mean wind speed acting on the bridge, a sufficient power can be acquired to drive the monitoring sensors. In addition, the battery cell can be used for storing the remaining power, which can be used when the operation of the energyharvesting system is difficult. The EM damper is used as the vibration-control system for external loads that cause excessive vibration of cables, such as typhoons and earthquakes. If excessive vibration occurs in the cable, the system is operated in the direction of maximizing the vibrationcontrol performance of the EM damper by minimizing the energy-harvesting performance.

Conclusions
In this study, the EM damper was first fabricated as a vibration-control method of a cable structure, and characteristic and performance tests were conducted. Contrary to previous studies, the performance of the EM damper was analyzed using a hybrid simulation incorporating the bridge cable properties to reflect the actual characteristics of the cable structure. Also, the energy-harvesting performance was analyzed by using the mean wind speed acting on the actual bridge. In order to

Conclusions
In this study, the EM damper was first fabricated as a vibration-control method of a cable structure, and characteristic and performance tests were conducted. Contrary to previous studies, the performance of the EM damper was analyzed using a hybrid simulation incorporating the bridge cable properties to reflect the actual characteristics of the cable structure. Also, the energy-harvesting performance was analyzed by using the mean wind speed acting on the actual bridge. In order to evaluate the feasibility of the damper, the performance of the EM damper was evaluated by reflecting actual wind load.
In order to evaluate the applicability of the EM damper to the cable structure, a characteristic test of the damper and a performance evaluation using a hybrid simulation were carried out. As a result of the characteristic test, a linear damping force change was confirmed in the amplitude change of the external load, and nonlinear damping force change was observed in the frequency variation.
The hybrid simulation results showed that the EM damper had a certain level of damping performance for the cable structure. In the free-vibration test, the damper showed a maximum damping ratio of 0.0132 at zero external resistance. As the external resistance increased, the damping ratio gradually decreased. It was confirmed that the attenuation performance was more than 50% in both the frequency and time domains for the sinusoidal input of the first to third natural frequencies.
The RMS displacement reduction of the cable was 42-44% in the case of wind load. As the magnitude of the external resistance increased, the cable vibration-control performance decreased, and the magnitude of the power generated by the damper gradually increased.
The applicability of the EM damper to a cable structure needs to be approached, not only from the vibration-control performance, but also from the maintenance aspect using the generated power. The vibration-control performance of the damper was insufficient compared to conventional damper systems. However, as excessive vibration of a cable can be reduced to an appropriate level, it could be applied to a cable structure in which vibration is continuously generated. In addition, as the damping coefficient of the damper can be adjusted using an external resistance, it could be utilized as a semi-active system, and better vibration-control performance is expected.
The energy-harvesting performance of the damper gradually increased as the external resistance increased. However, as the external resistance increased, the vibration-control performance of the damper decreased; therefore, it is necessary to maintain a proper vibration-control performance level.
According to the KMA, the mean wind speed was~4.8 m/s at the 2nd Jindo bridge, and 80% of the one-minute wind speed was more than 3 m/s. This can be expected to generate sufficient power to operate a wireless sensor for monitoring. The wind speed above 10 m/s was~4.6%. Most wind speeds are expected to be used for energy harvesting. As the power consumption of low-power sensors used for civil structure monitoring is several tens of mW, the power obtained from the damper is expected to be used as a power source for monitoring sensors.