# In Situ Measurement of Seeking Speed and Seeking Induced Head-Disk Interface Instability in Hard Disk Drives

Authors: Yu Wang; Xiongfei Wei; Yanyang Zi; Kwok-Leung Tsui
Journal: Shock and Vibration
 (2015)

---

## Abstract

This paper investigated the instability of head-disk interface caused by the voice coil motor (VCM) end crashing the crash stop during the seeking of magnetic head. To make the whole process of that clear, an in situ measurement method based on maximum likelihood estimation and extended Kalman filter for seeking speed at component level was developed first and was then calibrated by a high speed camera. Given a crash between VCM end and crash stop that may be a consequence of the continuous increasing seeking speed, the seeking speed was carefully controlled by using our developed method to find a critical value that may induce vigorous head-disk interface instability. Acoustic emission sensor and laser Doppler vibrometer were used to capture the transient dynamic behaviors of magnetic head when the crash is happening. Damage analysis and mode identification were carried out to reveal the relationship between the damage patterns on disk surface and head dynamics. The results of this study are helpful to optimize the track seeking profile during the HDD operation, as well as the design of components such as head and head arm.


---

## Body

## 1. Introduction

The areal density of the hard disk drive (HDD) is expected to reach multiseveral Tb/in2 in next few years [1, 2]. At such high density, the data transfer needs to further speed up due to the competition pressure from the solid-state drive (SSD) and the requirement of users. The access time or response time of a rotating drive that measures the time of HDD components’ response before the drive transferring data is mainly related to the mechanical nature of the rotating disks and moving heads. However, due to the limitations of mechanical components in HDDs, low-speed motion has become one of the chief obstacles for further increasing the speed of data transfer [3, 4]. Track following and seeking is a key step in the loop of data transfer in HDDs, and thus increasing seeking speed is raised as a critical task in HDDs.Generally, head seeking is implemented by inputting a digital current into the voice coil motor (VCM) to drive the head arm assembly move. To provide a specific speed, an available method is to add a constant current to VCM for a prescribed amount of time. However, this method may not achieve the precise speed control due to the disturbances, such as friction of the pivot bearings and the resistance force of data flex cable [5, 6]. An alternative method taking advantage of back electromotive force (EMF) was proposed to create a negative acceleration in order to produce the VCM velocity which is widely adopted in current HDDs [7–9]. The mechanism of this method is that a back EMF representing the velocity of the VCM will be generated when the power is cut off. The seeking speed can be attained by adjusting the back EMF.For a safety consideration, two crash stops were designed to limit the rotation of the VCM to prevent the head from moving off of the disk or hitting the spindle during head seeking [10]. No matter which manner used for VCM control, when the seeking speed is too high, it has a risk that the end of VCM may crash the crash stop. Such a crash may induce the HDD instability such as head slap on disk and head/disk damage. Although this topic is essential, there is rarely report on that in literature. It is therefore necessary to have a better understanding for the whole process of that, especially at component level, which is generally followed by system integration and thus is highly related to the component design and the yielding of the final products.Given the above considerations, this paper investigated the instability of HDI caused by the VCM crashing the crash stop during head seeking at the component level. To investigate the effect of seeking speed on anti-impact performance of head-disk interface, a maximum likelihood estimation (MLE) and extended Kalman filter (EKF), MLE-EKF, method was developed to estimated and control the speed flexibly. The EKF method is an improved method of the Kalman filter that is incapable of handling the nonlinear problems [11, 12]. Then, different speeds were carefully controlled to find a critical value that may induce HDI instability. Acoustic emission (AE) sensor and laser Doppler vibrometer (LDV) were used to capture the transient dynamic behaviors of head when the crash is happening. Damage analysis and mode identification were carried out to reveal the relationship between the damage patterns on disk surface and head dynamics.


## 2. Design and Calibration of Seeking Speed

### 2.1. Mathematic Model

Before implementing the experiment, the seeking speed of head should be well controlled and calibrated. To this end, a drive-based tester was used to implement the seeking test. Head seeking controlling method by the back EMF-based control technology applies the back EMF that is generated by coil movement in the magnetic field to drive the VCM to move. A simplified electronic circuit is shown in Figure1, and the electrical dynamics of the VCM can be described as [13](1)Lcdictdt+Rcict+emt=ect,where Lc is the inductance of coil, Rc is the resistance of coil, ec is applied voltage, and em is the induced back EMF in the coil due to the permanent magnets. According to Lenz’s law, the running of VCM generates a back EMF proportional to rotational velocity. Therefore, the back EMF can be calculated by em=Keω, and Ke is the back EMF constant.Figure 1
Simplified electronic circuit of a typical VCM.According to the mechanical dynamic model of the VCM, the current can be transferred to the rotation speed by [14, 15] (2)Ktict-Tf-Kc·x=Jω˙t,where Kt is the torque constant, ic is the coil current, J is the inertia of head arm, Tf is the torque generated by pivot friction force, Kc is the stiffness of flex cable, x is the displacement, and ω is the angular velocity. The values of Kt and J used in this experiment are 0.026 Nm/A and 0.1309 × 10−5 kgm2, respectively. By integrating ω˙, the rotation speed of arm can be obtained. Given the arm length l, the head seeking speed can be calculated. It is therefore necessary to know the implemented current and estimate the key parameters including Kt, Tf, and J before identifying the seeking speed.A programmable power source was used to control the current of VCM so as to control the seeking speed. Figure2 shows the designed waveform of VCM current and the measured current profile by an oscilloscope. The measured current profile consists of three phases. Phase I provided the power for VCM to disengage the latch. Phase II generated the seeking speed. Phase III was produced by the back EMF. It can be observed at phase III that a sudden drop occurred at 0.018 s, implying that a crash happened. A more detailed explanation was given in Section 3. When the current turned off, the rotational speed of VCM can be directly determined by back EMF. Using this property, the unknown parameters Tf and Kc can be estimated accordingly.Figure 2
The measured VCM current profiles during seeking.


### 2.2. Parameter Estimation

To estimate the parameters in (2), it is necessary to reconsider the mechanics of VCM. When the VCM current is cut off, (2) becomes (3)Tf+Kc·x=-Jω˙t,assuming that the displacement of flex cable has a linear relationship with the duration(4)x=a+bt.By integratingω˙(t), the angular speed of VCM can be calculated by (5)ω^=ω0-1JTf+Kcat-12Kcbt2.For the sake of simplicity, (5) can be simplified by combining the parameters, that is, letting C1=(1/J)Tf+Kca, C2=Kcb/J, and (5) become(6)ω^=ω0-C1t-12C2t2.Due to the disturbance, the measurement always companied with measurement error. The measurement error can be regarded as an independent and identically distributed normal random variable. Given the measurement errorε~N(0,σ2), the angular speed can be written as(7)ω=ω^+ε.The angular speed followsω~N(ω^,σ2). To estimate the parameters C1, C2, and σ, the maximum likelihood estimation method was used, which then turns to the logarithmic likelihood function of the angular speed:(8)log⁡L=∑i=1mni2-ln⁡2π-2ln⁡σ-∑k=1niωi-ω0+C1t+1/2C2t222σ2.By maximizing the log-likelihood function, the parametersC1, C2, and σ can be obtained.Three VCM current waveforms were used to build the likelihood function, as illustrated in Figure3. By maximizing the log-likelihood function (8), the estimated parameters are C1=245.9, C2=0, and σ=0.54. The seeking speed can be calculated by integrating ω˙ (when VCM current was supplied) or by em/Ke (when VCM current was cut off).Figure 3
Three VCM current waveforms to generate different seeking speed.


### 2.3. Real-Time Speed Estimation

We treated these estimated parameters as the initial parameters for an EKF; a real-time VCM speed can be obtained through the prediction and update steps in the EKF. We assume the state of the EKF isx=C1C2ω0. According to the mechanical dynamic model of the VCM, the state and measurement function of the EKF can be expressed as follows.State model is(9)xk=fxk-1∣k-1,k-1+qk-1,fxk-1∣k-1,k-1=100010001C1,k-1C2,k-1ω0,k-1,where qk-1=vC1vC2vω0, vC1~N(0,σC1), vC2~N(0,σC2), and vω0~N(0,σω0).Measurement model is(10)ωk=hxk,k+v,hxk,k=ω0,k-1-C1,kt-12C2,kt2,where v~N0,σy.Based on this state-space function, the first order EKF can be used to smooth the raw data. The steps for the first order EKF are as follows [12].Predict(11)xk∣k-1=fxk-1∣k-1,k-1,Pk∣k-1=Fxxk-1∣k-1,k-1Pk-1∣k-1FxTxk-1∣k-1,k-1+Qk-1,where Fx is defined to be ∂f/∂x.Update(12)vk=ωk-hxk,k,Sk=Hxxk-1∣k-1,k-1Pk-1∣k-1HxTxk-1∣k-1,k-1+Rk,Kk=Pk∣k-1HxTxk-1∣k-1,k-1Sk-1,xk∣k=xk∣k-1+Kkvk,Pk∣k=Pk∣k-1-KkSkKkT,where Hx is defined to be ∂h/∂x.To illustrate the effectiveness of the proposed method, the EMF induced speed was estimated by EKF. As shown in Figure4, the estimated speed follows the measured speed well as the maximum relative error is 8%. This indicates that our proposed speed estimation method is valid for estimating the head speed.Figure 4
Real-time estimation of the head speed by EKF.


### 2.4. Speed Calibration by High Speed Camera

To verify the aforementioned calculated seeking speed, a high speed camera taking the snapshots for head arm with a speed of 4500 frames per second was used. The head arm pattern at any location, as shown in Figure5, was continuously captured by the camera.Figure 5
The pattern captured by high speed camera.Using a National Instruments image identification software to measure the angular displacementΔθ of the featured pattern of head arm at a small time interval Δt during seeking process, the speed can be calculated by Δθ/Δt.A correlation between our estimated speed and the measured speed by camera was conducted, as illustrated by Figure6. As observed from this figure, a good correlation exists between the estimated and the measured results as R2 is approximately 0.95. This means the seeking speed can be well quantified by our estimated speed.Figure 6
The seeking speed calibrated by high speed camera.


## 2.1. Mathematic Model

Before implementing the experiment, the seeking speed of head should be well controlled and calibrated. To this end, a drive-based tester was used to implement the seeking test. Head seeking controlling method by the back EMF-based control technology applies the back EMF that is generated by coil movement in the magnetic field to drive the VCM to move. A simplified electronic circuit is shown in Figure1, and the electrical dynamics of the VCM can be described as [13](1)Lcdictdt+Rcict+emt=ect,where Lc is the inductance of coil, Rc is the resistance of coil, ec is applied voltage, and em is the induced back EMF in the coil due to the permanent magnets. According to Lenz’s law, the running of VCM generates a back EMF proportional to rotational velocity. Therefore, the back EMF can be calculated by em=Keω, and Ke is the back EMF constant.Figure 1
Simplified electronic circuit of a typical VCM.According to the mechanical dynamic model of the VCM, the current can be transferred to the rotation speed by [14, 15] (2)Ktict-Tf-Kc·x=Jω˙t,where Kt is the torque constant, ic is the coil current, J is the inertia of head arm, Tf is the torque generated by pivot friction force, Kc is the stiffness of flex cable, x is the displacement, and ω is the angular velocity. The values of Kt and J used in this experiment are 0.026 Nm/A and 0.1309 × 10−5 kgm2, respectively. By integrating ω˙, the rotation speed of arm can be obtained. Given the arm length l, the head seeking speed can be calculated. It is therefore necessary to know the implemented current and estimate the key parameters including Kt, Tf, and J before identifying the seeking speed.A programmable power source was used to control the current of VCM so as to control the seeking speed. Figure2 shows the designed waveform of VCM current and the measured current profile by an oscilloscope. The measured current profile consists of three phases. Phase I provided the power for VCM to disengage the latch. Phase II generated the seeking speed. Phase III was produced by the back EMF. It can be observed at phase III that a sudden drop occurred at 0.018 s, implying that a crash happened. A more detailed explanation was given in Section 3. When the current turned off, the rotational speed of VCM can be directly determined by back EMF. Using this property, the unknown parameters Tf and Kc can be estimated accordingly.Figure 2
The measured VCM current profiles during seeking.


## 2.2. Parameter Estimation

To estimate the parameters in (2), it is necessary to reconsider the mechanics of VCM. When the VCM current is cut off, (2) becomes (3)Tf+Kc·x=-Jω˙t,assuming that the displacement of flex cable has a linear relationship with the duration(4)x=a+bt.By integratingω˙(t), the angular speed of VCM can be calculated by (5)ω^=ω0-1JTf+Kcat-12Kcbt2.For the sake of simplicity, (5) can be simplified by combining the parameters, that is, letting C1=(1/J)Tf+Kca, C2=Kcb/J, and (5) become(6)ω^=ω0-C1t-12C2t2.Due to the disturbance, the measurement always companied with measurement error. The measurement error can be regarded as an independent and identically distributed normal random variable. Given the measurement errorε~N(0,σ2), the angular speed can be written as(7)ω=ω^+ε.The angular speed followsω~N(ω^,σ2). To estimate the parameters C1, C2, and σ, the maximum likelihood estimation method was used, which then turns to the logarithmic likelihood function of the angular speed:(8)log⁡L=∑i=1mni2-ln⁡2π-2ln⁡σ-∑k=1niωi-ω0+C1t+1/2C2t222σ2.By maximizing the log-likelihood function, the parametersC1, C2, and σ can be obtained.Three VCM current waveforms were used to build the likelihood function, as illustrated in Figure3. By maximizing the log-likelihood function (8), the estimated parameters are C1=245.9, C2=0, and σ=0.54. The seeking speed can be calculated by integrating ω˙ (when VCM current was supplied) or by em/Ke (when VCM current was cut off).Figure 3
Three VCM current waveforms to generate different seeking speed.


## 2.3. Real-Time Speed Estimation

We treated these estimated parameters as the initial parameters for an EKF; a real-time VCM speed can be obtained through the prediction and update steps in the EKF. We assume the state of the EKF isx=C1C2ω0. According to the mechanical dynamic model of the VCM, the state and measurement function of the EKF can be expressed as follows.State model is(9)xk=fxk-1∣k-1,k-1+qk-1,fxk-1∣k-1,k-1=100010001C1,k-1C2,k-1ω0,k-1,where qk-1=vC1vC2vω0, vC1~N(0,σC1), vC2~N(0,σC2), and vω0~N(0,σω0).Measurement model is(10)ωk=hxk,k+v,hxk,k=ω0,k-1-C1,kt-12C2,kt2,where v~N0,σy.Based on this state-space function, the first order EKF can be used to smooth the raw data. The steps for the first order EKF are as follows [12].Predict(11)xk∣k-1=fxk-1∣k-1,k-1,Pk∣k-1=Fxxk-1∣k-1,k-1Pk-1∣k-1FxTxk-1∣k-1,k-1+Qk-1,where Fx is defined to be ∂f/∂x.Update(12)vk=ωk-hxk,k,Sk=Hxxk-1∣k-1,k-1Pk-1∣k-1HxTxk-1∣k-1,k-1+Rk,Kk=Pk∣k-1HxTxk-1∣k-1,k-1Sk-1,xk∣k=xk∣k-1+Kkvk,Pk∣k=Pk∣k-1-KkSkKkT,where Hx is defined to be ∂h/∂x.To illustrate the effectiveness of the proposed method, the EMF induced speed was estimated by EKF. As shown in Figure4, the estimated speed follows the measured speed well as the maximum relative error is 8%. This indicates that our proposed speed estimation method is valid for estimating the head speed.Figure 4
Real-time estimation of the head speed by EKF.


## 2.4. Speed Calibration by High Speed Camera

To verify the aforementioned calculated seeking speed, a high speed camera taking the snapshots for head arm with a speed of 4500 frames per second was used. The head arm pattern at any location, as shown in Figure5, was continuously captured by the camera.Figure 5
The pattern captured by high speed camera.Using a National Instruments image identification software to measure the angular displacementΔθ of the featured pattern of head arm at a small time interval Δt during seeking process, the speed can be calculated by Δθ/Δt.A correlation between our estimated speed and the measured speed by camera was conducted, as illustrated by Figure6. As observed from this figure, a good correlation exists between the estimated and the measured results as R2 is approximately 0.95. This means the seeking speed can be well quantified by our estimated speed.Figure 6
The seeking speed calibrated by high speed camera.


## 3. Experimental Setup for Crashing Characterization

Based on back EMF control technology, a component-level evaluation system for seeking was developed. The schematic of the experimental configuration is shown in Figure7. An acoustic emission sensor was mounted on the top head arm to monitor the arm dynamics. A laser Doppler vibrometer (LDV) was used to detect the head-disk contact and monitor the transient head dynamics when a crash is happening. In this experiment, a typical Pemto slider (1.235 × 0.7 × 0.23 mm3) covered by a 2 nm diamond-like carbon (DLC) film was tested. The nominal flying attitude was 10 nm. The tested disk was a carbon coated 3.5-inch disk platter with a rotation speed of 5400 rpm. Signals were collected by an oscilloscope and a spectrum analyzer. The disk damage was inspected by an optical surface analyzer (OSA).Figure 7
The schematic of the experiment setup.Figure8 shows the signals monitored during the whole seeking process from inner diameter (ID) to outer diameter (OD) and then to ID. As observed from this figure, at phase III that the VCM current experienced, a vigorous vibration that occurred after the maximum back EMF was just generated. This vigorous vibration indicates a crash happening. Due to the crash, the seeking speed has a significant loss, leading to a sudden drop of the back EMF. By calculating the displacement that the head moved, this vigorous vibration can be identified as the crashing point between VCM end and OD crash stop. A very small response was observed at ID crash stop, indicating that the seeking speed has a significant loss after the crashing at OD crash stop.Figure 8
The monitored signals during seeking from inner diameter (ID) to outer diameter (OD), and then to ID. CS is the abbreviation of crash stop.


## 4. The Effect of Seeking Speed on HDD Instability

### 4.1. Seeking Speed and Crashing Energy

The crash between VCM end and crash stop may induce a vigorous vibration of head and arm and further result in HDD instability such as head slap on disk and head/disk damage depending on the crashing energy. Seeking speed is the main cause of the crashing energy. To reveal the relationship between speed and crashing energy, crashing energy characterized by AE signal under different seeking speeds was investigated, as illustrated by Figure9. It can be observed that, with the increase of seeking speed during crashing, the crash-induced vibration increases quickly.Figure 9
AE signals under different crashing speeds. ips is the abbreviation of inch per second.To describe the crashing energy, the root mean squares (RMS) of AE signals were calculated under different crashing speeds. A clear increasing tread can be observed in Figure10. From around 160 ips, the energy has begun to increase sharply, which indicated that HDD may become instable after this critical value. Further increasing the speed may lead to head/disk damage.Figure 10
The energy generated by crashing under different seeking speeds.


### 4.2. Head Dynamics

To monitor the head transient dynamics when crashing just happened, the laser spot of LDV was fixed on the position of OD crash stop to capture the vibration of head trailing edge center when head passed OD position. During the experiment, the laser will scan the trailing edge of head back when the head passed. Figure11 shows the laser scanned area on the head back when head passed OD.Figure 11
Laser scanned area on head when it passed to OD position.A spectrum analyzer was used to acquire and analyze the laser signals with the cutoff of the frequency range of 10 kHz–200 kHz. Three different crashing speeds were conducted in order to reveal the relationship between the head dynamics and speed, as illustrated in Figure12. As observed from this figure, when the crashing speed reached 162 ips, a large head vibration happened, while at 66 ips and 106 ips, there is nearly no significant response. This suggests that, at 162 ips, the air bearing may collapse and leads to HDI instability. The frequency spectrum of head response at 162 ips shows that the dominant modes are at 29 kHz and 63 kHz. The responses of these two modes at 162 ips are about 10 times bigger than that of the other two speeds.Figure 12
The frequency spectrum of LDV signal during crashing at OD crash stop.To reveal the dynamic nature of the head during crashing, a finite element model of whole head arm was built to identify the mode shapes of the two modes. The air bearing was modeled by four springs whose stiffness was identified by a parameter identification code [16]. The finite element model and the constraints are shown in Figure 13. Both modal and harmonic analysis were carried out using the finite element model. The mode shape of head under a lateral excitation (a unit force on coil end) showed that a torsion mode of head gimbal and the rolling mode of the air bearing dominated the head side response. The frequencies of these two modes (29.8 kHz and 62.3 kHz), as illustrated in Figure 14, are in accordance with the measured frequencies by LDV, respectively.Figure 13
The finite element model of head arm.Figure 14
Head side response with a lateral excitation.


### 4.3. Damage Inspection of Disk Surface

After the experiment, the tested disks were inspected by OSA. Figure15 shows the tested disk at crashing speed of 162 ips. As seen from this figure, some special patterns (damage) were observed, providing the evidence of air bearing collapse. The width of the patterns is about 0.7 mm, which is equivalent to the head width, indicating that the head rolling motion dominates the head vibration when the arm end is crashing crash stop. This result is in good accordance with the mode shape that was identified by finite element model.Figure 15
The OSA inspection result of disk after crashing test (crashing speed = 162 ips).


## 4.1. Seeking Speed and Crashing Energy

The crash between VCM end and crash stop may induce a vigorous vibration of head and arm and further result in HDD instability such as head slap on disk and head/disk damage depending on the crashing energy. Seeking speed is the main cause of the crashing energy. To reveal the relationship between speed and crashing energy, crashing energy characterized by AE signal under different seeking speeds was investigated, as illustrated by Figure9. It can be observed that, with the increase of seeking speed during crashing, the crash-induced vibration increases quickly.Figure 9
AE signals under different crashing speeds. ips is the abbreviation of inch per second.To describe the crashing energy, the root mean squares (RMS) of AE signals were calculated under different crashing speeds. A clear increasing tread can be observed in Figure10. From around 160 ips, the energy has begun to increase sharply, which indicated that HDD may become instable after this critical value. Further increasing the speed may lead to head/disk damage.Figure 10
The energy generated by crashing under different seeking speeds.


## 4.2. Head Dynamics

To monitor the head transient dynamics when crashing just happened, the laser spot of LDV was fixed on the position of OD crash stop to capture the vibration of head trailing edge center when head passed OD position. During the experiment, the laser will scan the trailing edge of head back when the head passed. Figure11 shows the laser scanned area on the head back when head passed OD.Figure 11
Laser scanned area on head when it passed to OD position.A spectrum analyzer was used to acquire and analyze the laser signals with the cutoff of the frequency range of 10 kHz–200 kHz. Three different crashing speeds were conducted in order to reveal the relationship between the head dynamics and speed, as illustrated in Figure12. As observed from this figure, when the crashing speed reached 162 ips, a large head vibration happened, while at 66 ips and 106 ips, there is nearly no significant response. This suggests that, at 162 ips, the air bearing may collapse and leads to HDI instability. The frequency spectrum of head response at 162 ips shows that the dominant modes are at 29 kHz and 63 kHz. The responses of these two modes at 162 ips are about 10 times bigger than that of the other two speeds.Figure 12
The frequency spectrum of LDV signal during crashing at OD crash stop.To reveal the dynamic nature of the head during crashing, a finite element model of whole head arm was built to identify the mode shapes of the two modes. The air bearing was modeled by four springs whose stiffness was identified by a parameter identification code [16]. The finite element model and the constraints are shown in Figure 13. Both modal and harmonic analysis were carried out using the finite element model. The mode shape of head under a lateral excitation (a unit force on coil end) showed that a torsion mode of head gimbal and the rolling mode of the air bearing dominated the head side response. The frequencies of these two modes (29.8 kHz and 62.3 kHz), as illustrated in Figure 14, are in accordance with the measured frequencies by LDV, respectively.Figure 13
The finite element model of head arm.Figure 14
Head side response with a lateral excitation.


## 4.3. Damage Inspection of Disk Surface

After the experiment, the tested disks were inspected by OSA. Figure15 shows the tested disk at crashing speed of 162 ips. As seen from this figure, some special patterns (damage) were observed, providing the evidence of air bearing collapse. The width of the patterns is about 0.7 mm, which is equivalent to the head width, indicating that the head rolling motion dominates the head vibration when the arm end is crashing crash stop. This result is in good accordance with the mode shape that was identified by finite element model.Figure 15
The OSA inspection result of disk after crashing test (crashing speed = 162 ips).


## 5. Conclusion

This paper investigated the effect of seeking speed on HDD instability. The head seeking was controlled by a back EMF-based technology. A mathematical method based on MLE-EKF for deriving the seeking speed was developed. AE sensor and LDV were used to characterize the crashing energy and capture the head transient dynamics at the different speeds, respectively. With the help of finite element model and OSA inspection, the rolling mode of air bearing was found to dominate the head dynamics. This mode will lead to a serious damage of disk when a high seeking speed was adopted. The results revealed in this paper will help to optimize the drive-level track seeking profile, as well as the design of components such as head and head arm.



---

*Source: 989854-2015-11-15.xml*