Table of Contents Author Guidelines Submit a Manuscript
Advances in Acoustics and Vibration
Volume 2015 (2015), Article ID 707043, 8 pages
Research Article

Modeling of Laterally Sliding Motion of a Magnetic Clamp

Département de Génie Mécanique, École Polytechnique, 2900 Boulevard Édouard-Montpetit, Montréal, QC, Canada H3T 1J4

Received 5 October 2015; Revised 26 November 2015; Accepted 29 November 2015

Academic Editor: Kim M. Liew

Copyright © 2015 A. Mahmud 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.


A sliding magnetic clamp is used to hold a thin aluminum panel during a milling operation. The design includes a permanent magnet group follower (slave module) which slides laterally over the panel attracted by another permanent magnet group (master module) attached to the industrial robot end effector from the machined side of the panel. The lateral sliding motion of the slave module in response to the master module motion is studied using a transfer function based motion model established considering the lateral magnetic stiffness. The model is validated experimentally.

1. Introduction

This paper proposes a model for the lateral sliding motion of the master-to-slave module by magnetic attraction forces. The grasping end effector is intended to mill pockets in thin aircraft fuselage skin panels [1] with the support of magnetic clamping. The system is developed for holding a panel during milling operation. The clamping system includes two magnetic modules which are placed on either side of the panel. Figure 1 shows a general schematic of the grasping machining end effector [2].

Figure 1: General schematic diagram of the grasping machining end effector concept.

The master module has three cylindrical permanent magnet sets, each set comprising three magnets (Figure 2). The milling cutter is at the center of the triangle formed by these three magnet sets. The slave module has a similar magnet arrangement. The master and slave work together as a clamp.

Figure 2: Master with three extended legs (a). Slave with three sets of permanent magnets (b) [7].

Elles and Lemarquand [3] describe radial stability characteristics in magnetic couplings used to transmit a torque through an air gap separation wall without any mechanical contact between the two rotors. The paper proposes two cylindrical-airgap magnetic couplings structures that transmit a torque but no radial displacement or vibration from the inner rotor to the outer rotor.

A cylindrical permanent magnet can generate both axial and radial forces. Lateral sliding is the result of radial forces. Yonnet [4] defined lateral magnetic stiffness while discussing different magnetic properties for permanent magnetic bearings and couplings.

The slave module laterally follows the master module motion during machining operation. This slave motion caused by the master-slave magnetic attraction force has not yet been modeled. Vokoun et al. [5] calculated the attraction force between two cylindrical permanent magnets with a common axis considering the impact of lateral displacement on the axial attraction force. Later Mahmud et al. [6] modeled the lateral attraction force between cylindrical permanent magnets at different axial gap considering lateral displacement. A model is needed to describe the slave motion resulting from lateral magnetic attraction force between master and slave module. Without a model it is not possible to predict whether the slave would accurately follow the master or experience sluggish or stick-slip motion.

This paper develops the transfer function between the master module motion and slave module motion. First a simplified model is prepared for the mechanical and magnetic components. Then classical motion equations are used to describe the slave motion where lateral magnetic stiffness has been incorporated. Experiments are then conducted to validate the transfer function model. Finally conclusions are drawn.

2. Transfer Function Modeling

A simplified 2D schematic of the system is presented in Figure 3. On top, the master module is fixed to the milling spindle outer body. The spindle is moved by a robot or a CNC machine tool. Each magnet of the slave module is subjected to a magnetic pull force , out-of-plane, in the axial direction and a lateral force , in-plane, in the radial direction [8]. Both the master and slave modules are subjected to frictional forces and , respectively.

Figure 3: Master-slave magnetic forces in a simplified 2D view.

This model is a two-degree-of-freedom (in-plane and out-of-plane) mechanical system. The master module is the driver of the system. The slave module follows the master module due to the radial magnetic attraction force .

The machine tool or robotic arm moves the master and causes a lateral displacement . An displacement response is imparted to the slave by the radial magnetic force where is the lateral magnetic stiffness as defined by Yonnet [4].

A free body diagram for each of the system components (master, panel, and slave) is shown in Figure 4. The master is held by the robot by axial force . Master motion is dependent on the applied lateral force . In this analysis the robot is assumed to be sufficiently strong and stiff to generate the necessary force to generate the desired input displacement . For the slave module the axial contact force is and for the master it is which is equal to the applied axial magnetic force . The cutter applies a milling thrust force and radial force to the thin skin panel.

Figure 4: Free body diagram for panel-master-slave.

First let us consider a simplified model without friction. If the master and slave are just pulled apart by holding the slave by a string arrangement (Figures 7 and 8) without placing any plate between them there is no friction. All vertical force components can be ignored just for the time being since they do not have any contribution to lateral motion. Milling forces are ignored in the absence of milling operation. After elimination of all friction forces ( and ) and vertical forces (, , and ), only the radial magnetic force remains. The free body diagram is also simplified as in Figure 5.

Figure 5: Modified free body diagram for master and slave.

Assuming the slave has a mass and damping coefficient (air viscosity, string attachments) and the magnetic lateral stiffness an equivalent single degree of freedom mass-spring-damper mechanical system is used to model the system as shown in Figure 6.

Figure 6: Master-slave equivalent mechanical system without the panel and associated friction and milling force.
Figure 7: Lateral stiffness measurement by manual lateral force measurement with spring balance.
Figure 8: Test setup for system parameter identification.

A single degree of freedom damped spring mass system with a moving support or base (base excitation) has been described by Hundal [9]. Here the master motion is the system input so the differential equation of motion without friction is Applying Laplace transform to this equation gives (assuming zero initial conditions) and the system transfer function becomes

We now include friction to complete the transfer function modeling. This requires considering the skin panel between the master and slave.

Sliding friction or dry friction which is considered as coulomb friction is modeled as a viscously damped system with an equivalent viscous damping coefficient [10]:where is the approximate amplitude of the steady state motion assuming both the slave and panel surface are well lubricated ensuring a viscous medium. Consider

Here is the frequency of motion, is the natural frequency of the slave, and is the driving force applied on the master causing the motion.

Equations (4) and (5) are based on the applied force applied on the master. Applied driving force input is equal to the master displacement multiplied by stiffness (). So (5) can be rewritten considering master displacement as input instead of the applied force:

The slave applies an axial magnetic force on the skin panel. For the slave the equivalent viscous damping coefficient is . The master-to-slave transfer function (3) is rewritten with equivalent viscous damping coefficient to consider the fiction interfering during the milling operation. Consider

3. System Parameter Identification

The slave mass viscous damping coefficient and lateral magnetic stiffness need to be identified. The slave mass is 0.67 kg as measured on a weighing scale.

Direct force measurement with a spring balance provides stiffness value and subsequent calculation yields the undamped natural frequency. Figure 7 shows the preliminary test setup where a gap of 5 mm is maintained by string arrangement between the master and slave which represent the thickness of a panel to machine. For a 0.003 m lateral (radial) displacement 43.51 N is needed yielding stiffness of 14500 N/m. The undamped natural frequency is calculated as

Using the same string arrangement setup (Figure 8) a series of impact hammer tests were conducted to identify the damping and lateral magnetic stiffness for different gap values.

While making impact test on the slave module it was ensured that the hammer impact direction remained aligned with the center line as shown in Figure 9. This additional precaution ensures that the impact force did not create any rotational torque between the master and slave modules.

Figure 9: Impact force direction on the slave module during impact hammer test.

One miniature single axis accelerometer (model number 352C22, sensitivity 10 mV/g) and an impact hammer (model number 086D05, sensitivity 0.23 mV/N) both from PCB PIEZOTRONICS were connected to MATLAB data acquisition program via NI data acquisition card DAQCard-A1-16E-4. As soon as the hammer impacts the slave, vibrations are recorded. Acquired acceleration data were integrated twice to obtain the displacement data with respect to time. The calculated displacement profiles have been analyzed for slave damping and stiffness value.

Figure 10 shows the displacement pattern for a 5 mm gap between the master and slave. Immediately after the impulse force the slave continuously vibrates at its damped natural frequency.

Figure 10: Displacement graph for a 5 mm gap between the master and slave.

Logarithmic decrement is applied for identifying the damping ratio [11]:where is the displacement of the th peak at the time and is the number of complete oscillations considered. Brincker et al. [12] convert the logarithmic decrement value to the damping ratio :

The damped natural frequency is calculated using and the oscillation period :

Then the required lateral magnetic stiffness coefficient value is calculated as follows:

Finally, the viscous damping coefficient is

A stiffness value of  N/m (Table 1) associated with a viscous damping value of 0.9090 N-sec/m has been calculated for the slave module for a gap of 5 mm between the master and slave (Figure 11). This is in good agreement with the value of  N/m obtained from the static load test.

Table 1: System identification experiments results summary.
Figure 11: A 5 mm gap between the master and slave is maintained by a string arrangement.

The measured mass and stiffness values yielded a calculated natural frequency of 147.13 rad/sec for a gap of 5 mm where the accelerometer based measurement and mathematical analysis results in a natural frequency of 149.46 rad/sec.

The procedure was repeated for five gap values. Table 1 shows the results. The table also includes the natural frequency values and corresponding damping coefficients.

Damping coefficient generally decreased with an increased gap. In case of stiffness a decreasing pattern is prominent. Figure 12 shows that the stiffness value is decreasing as the gap increases.

Figure 12: Stiffness reduces as the master-slave gap increases (thicker plate).

During a real milling operation the master and slave will hold the aluminum panel and friction will occur during the sliding motion. So, an experimental setup has been realized to test the friction. The plate is clamped cantilevered. Both the manual (spring balance) and accelerometer based process are used for this setup to get the friction parameter.

A thin aluminum plate of thickness 2.54 mm with 1.25 mm slide pad on both sides makes a gap of 5.04 mm. The impact hammer test procedure is completed. Since the captured accelerometer data did not show any regular sinusoidal vibration pattern as shown in Figure 14, logarithmic decrement based equation was not feasible for identifying the system parameter.

Figure 13: Plate between master and slave placed where accelerometer is positioned on the slave.
Figure 14: Slave displacement result for the impact hammer test of Figure 13.

Another test has been conducted with a spring balance to measure the lateral and axial magnetic force as shown in Figure 15. At this condition spring balance direct lateral force reading is 52.92 N, and axial magnetic force, , is 66.64 N. Polyoxymethylene pad over aluminum panel has a friction coefficient of 0.17 [13]. The frictional force is 11.32 N which is used for the calculation of equivalent viscous damping coefficient. Table 2 lists the axial magnetic force between master and slave for plates of different thickness.

Table 2: Friction force at different plate thickness.
Figure 15: Test setup for axial force measurement by a spring balance.

4. Experiments and Results

The derived transfer function model does not include milling force; only magnetic attraction forces are required. So, the transfer function validation experiment is designed without milling operation.

The master was moved by manual push and pull over a skin panel with the slave following the master. Both the master and slave experienced friction.

Two identical accelerometers are attached to the master module and the slave module (Figure 16). The first test is performed with a thin aluminum plate (2.54 mm, pad included total master-slave gap of 5.04 mm).

Figure 16: Simultaneous acceleration data acquisition for both master and slave.

Table 1 was obtained for tests without a plate between master and slave. Still, the same stiffness value would be accepted even if a plate is placed between master and salve because the stiffness is a function of distance between master and slave. This is because the stiffness is expected to come from the magnetic force and a nonmagnetic plate does not influence the magnetic field. Inserting a plate would only change the damping because of friction. So an equivalent viscous damping coefficient has been calculated based on (4).

Data acquisition sampling frequency value is 10,000 Hz. Frictional force value from Table 2 has been used for calculating steady state slave motion amplitude, . Equivalent viscous damping coefficient has been calculated for each steady state slave displacement by (6). For a 2.54 mm thick skin panel the transfer function has been defined as follows based on the transfer function model (7). Consider

Figure 17 is obtained after plotting the master and slave motion response from accelerometer reading. On the top master displacement is presented. On the bottom slave displacement is presented in a solid line. An additional dotted line is also plotted to show the simulated slave displacement based on the derived master-to-slave transfer function.

Figure 17: Comparison of transfer function based simulated motion to the actual slave motion (master-slave gap 5.04 mm).

-square value [14] is calculated for the simulated slave displacement with respect to the measured slave displacement for a gap of 5.04 mm. An -square value of 0.89 has been calculated showing the acceptability of the proposed slave motion model.

Magnetic stiffness was used as the prime constituent of the proposed master-slave motion transfer function. Since magnetic stiffness is highly dependent on the gap between the master and slave, one more experiment has been accomplished with a thicker skin panel of 5.05 mm (total gap 7.55 mm) for a reduced stiffness value of  N/m. In this case the axial magnetic force is 50.96 N. The transfer function is still valid to predict the slave motion to an acceptable level, associated with an -square value of 0.8340 as shown in Figure 18.

Figure 18: Comparison of transfer function based simulated motion to the actual slave motion (master-slave gap 7.55 mm).

5. Conclusion

During the skin panel milling operation the robot drive master module will follow the instructed tool path based on the required pocket dimension, but the slave motion is dependent on the magnetic attraction force against its inertia and the frictional force. A master-to-slave motion transfer function is proposed to calculate the slave response. The proposed transfer function based motion model was able to predict the slave position based on the input master motion. Though this paper only considers the lateral sliding motion in case of real pocket milling operation there will be rotational motion too. An in-depth rotational motion analysis would be relevant as future work to predict the complete slave module motion response.

Spring balance direct measurement provided the lateral stiffness values later verified by logarithmic decrement formula based calculation. But the most critical part of the system parameter identification is damping values. While sliding over the skin panel the slave confronts coulomb friction. Equivalent viscous damping coefficient has been used for the transfer function.

The master-to-slave transfer function is based on three parameters, slave mass, lateral magnetic stiffness, and damping coefficient. This specific project had a slave module of 0.67 kg with an approximate lateral stiffness of 14500 N/m. Two independent experiments with different thickness skin panel have been tested to verify the transfer function model acceptability. Both of the test results satisfactorily validated the model suitability.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was supported by the project of “Machining of Skin Panels (Al or Al-Li),” CRIAQ MANU-412–NSERC RDCPJ 411911-10.


  1. S. Castagne, R. Curran, and P. Collopy, “Implementation of value-driven optimisation for the design of aircraft fuselage panels,” International Journal of Production Economics, vol. 117, no. 2, pp. 381–388, 2009. View at Publisher · View at Google Scholar · View at Scopus
  2. A. Mahmud, J. R. R. Mayer, and L. Baron, “Determining the minimum clamping force by cutting force simulation in aerospace fuselage pocket machining,” The International Journal of Advanced Manufacturing Technology, vol. 80, no. 9–12, pp. 1751–1758, 2015. View at Publisher · View at Google Scholar
  3. P. Elles and G. Lemarquand, “Analytical study of radial stability of permanent-magnet synchronous couplings,” IEEE Transactions on Magnetics, vol. 35, no. 4, pp. 2133–2136, 1999. View at Publisher · View at Google Scholar · View at Scopus
  4. J. P. Yonnet, “Permanent magnet bearings and couplings,” IEEE Transactions on Magnetics, vol. 17, no. 1, pp. 1169–1173, 1981. View at Publisher · View at Google Scholar · View at Scopus
  5. D. Vokoun, M. Beleggia, L. Heller, and P. Šittner, “Magnetostatic interactions and forces between cylindrical permanent magnets,” Journal of Magnetism and Magnetic Materials, vol. 321, no. 22, pp. 3758–3763, 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. A. Mahmud, J. Mayer, and L. Baron, “Magnetic attraction forces between permanent magnet group arrays in a mobile magnetic clamp for pocket machining,” CIRP Journal of Manufacturing Science and Technology, vol. 11, pp. 82–88, 2015. View at Publisher · View at Google Scholar
  7. A. Chouinard, Remplacement de l'usinage Chimique des Tooles Aeronautiques Minces en Aluminium par de l'usinage Mecanique, École Polytechnique de Montréal, Montréal, Canada, 2011.
  8. A. K. Ghosh and S. S. Hecker, “Stretching limits in sheet metals: in-plane versus out-of-plane deformation,” Metallurgical Transactions, vol. 5, no. 10, pp. 2161–2164, 1974. View at Publisher · View at Google Scholar · View at Scopus
  9. M. S. Hundal, “Response of a base excited system with Coulomb and viscous friction,” Journal of Sound and Vibration, vol. 64, no. 3, pp. 371–378, 1979. View at Publisher · View at Google Scholar · View at Scopus
  10. D. J. Inman and R. C. Singh, Engineering Vibration, vol. 3, Prentice-Hall, Upper Saddle River, NJ, USA, 2001.
  11. C.-H. Lamarque, S. Pernot, and A. Cuer, “Damping identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement—part 1: theory,” Journal of Sound and Vibration, vol. 235, no. 3, pp. 361–374, 2000. View at Publisher · View at Google Scholar · View at Scopus
  12. R. Brincker, C. Ventura, and P. Andersen, “Damping estimation by frequency domain decomposition,” in Proceedings of the 19th International Modal Analysis Conference (IMAC '01), Kissimmee, Fla, USA, February 2001.
  13. H. Unal and A. Mimaroglu, “Friction and wear behaviour of unfilled engineering thermoplastics,” Materials & Design, vol. 24, no. 3, pp. 183–187, 2003. View at Publisher · View at Google Scholar · View at Scopus
  14. A. C. Cameron and F. A. G. Windmeijer, “R-squared measures for count data regression models with applications to health-care utilization,” Journal of Business & Economic Statistics, vol. 14, no. 2, pp. 209–220, 1996. View at Google Scholar · View at Scopus