Design of a Noncontact Torsion Testing Device Using Magnetic Levitation Mechanism

: To facilitate torsion testing in special environments, a noncontact torsion testing device, where a noncontact torque can be transmitted to a specimen, was developed using magnetic levitation technology. A permanent magnetic gear is used to produce noncontact torque. In addition, four electromagnets and four attractive-type permanent magnetic bearings were employed to realize levitation; in more detail, the four electromagnets actively stabilized two levitation degrees of freedom (DoFs), while the four attractive-type permanent magnetic bearings passively stabilized four DoFs. Furthermore, a plant model considering the effect caused by the four attractive-type permanent magnetic bearings was built for the two levitation DoFs requiring active control. Based on the plant model, two PD-controllers were designed. Moreover, a control simulation was conducted to obtain appropriate PD-gains. Finally, experiments further validated the feasibility of the whole scheme, and it was proven that the device can apply a 0.126 N · m torque to the specimen while maintaining levitation.


Introduction
Magnetic levitation (Maglev) technology has aroused extensive attention due to its advantages of being friction-free, lubrication-free, and maintenance-free, and it has been used in many industrial applications, such as railway systems [1][2][3][4][5][6], rotating machinery [7][8][9][10][11][12], ultra-clean assembly line, and artificial heart pumps [13,14], which are the typical applications of maglev.In addition to these, there are some novel applications of maglev, such as maglev gyro [15], maglev globe [16], maglev lateral vibration attenuation system [17], permanent maglev system [18,19], and zero-power maglev system [20,21].Similarly, in our previous research, to cater to the testing requirements for the parts and materials working in special environments, we applied maglev to tension and bending testing devices [22,23].Special environments include vacuum, special gases, and special liquids, as well as environments with abnormal air pressure or temperature, which may greatly influence the mechanical performance of the parts or components; therefore, it is necessary to facilitate the mechanical test in special environments.E.g., high-purity aluminum [24], which is widely used in satellites, spacecraft, and other aerospace devices, has to work in vacuum environments because, in normal air, a layer of oxide film will form on its surface, impairing its mechanical performance, corrosion resistance, and adhesion.In addition, gas-lubricated bearings [25] need to work in helium environments to reduce air resistance and better heat dissipation during high-speed rotation because helium has a lower viscosity and higher thermal conductivity compared with ordinary air.Furthermore, silicon-based photonic chips [26], which are typically used in optical communication and sensing, need to work in nitrogen environments because nitrogen can help silicon-based photonic chips maintain optical performance and stability.Previously, if we wanted to test these parts or materials' mechanical properties in their original working environment, the entire material testing device had to be placed in a container in which the special environment was.On the contrary, with maglev technology, the specimens in the devices [22,23] can be applied with load in a noncontact way; as a result, only the specimens and a few parts need to be placed in the container where the special environment can be created, which avoids exposing sensitive components of the testing device (e.g., circuits in force sensors and actuators are sensitive to environment humidity) to the special environment, facilitating the conduct of the testing.The research [22,23] was inspired by the research [27] held by Okayama University.In the research [27], Naoya Tada and Hiroyasu Masago developed a device in which a noncontact tension force can be applied to one end of a specimen using permanent magnets.As the tension force can be transmitted to the specimen in a noncontact way, the specimen can be isolated from the fixture at the end, and the specimen can be put in a container where a liquid environment can be created, which will be helpful to liquid and humid environment's testing.However, the device in the research [27] only allows one end of the specimen to realize noncontact force; the other end still remains in contact force.Therefore, it is difficult to seal the specimen in a closed container, but a completely closed container is necessary to create a vacuum and gas environment.That is to say, this device is only helpful for testing in a liquid environment and not really helpful for testing in a gas or vacuum environment.To address this issue, in the research [22,23], we developed a completely noncontact tension testing device and a completely noncontact bending device using maglev technology.In addition, by designing various maglev mechanisms, various types of noncontact material testing devices can be developed.The Torsion test is one of the most common mechanical tests; it can directly measure the shear strength of materials, which is a very important parameter in some applications.The Torsion test is applicable to a wide variety of materials, such as metals, polymers, ceramics, and composites, and it is a short test that can usually be finished in a few minutes.Due to these advantages, it is widely used in metal processing, material research, medical equipment, and the automobile industry.Given the status and superiority of the Torsion test, in this paper, we focus on developing a noncontact torsion testing device using maglev technology.
Primarily, it is necessary to consider a control scheme for global levitation stability.Generally, there are six levitation degrees of freedom (DoF) for a rigid levitated object.If many levitation DoFs can be inherently stable, i.e., passively stable, a lot of hardware, such as displacement sensors, electromagnets, amplifiers, and controllers can be saved.Therefore, we want the designed structure to have as many passively stable levitation DoFs as possible.However, according to Earnshaw's theorem [28], it is necessary that at least one levitation DoF be actively controlled for complete non-contact.In this paper, to maximize passively stable levitation DoFs, a structure with four passively stable levitation DoFs was designed.In the structure, the passive stability is mainly realized by four attractive-type permanent magnetic bearings; each attractive-type permanent magnetic bearing consists of two axially magnetized permanent magnet rings.Furthermore, to produce noncontact torque, we employed one permanent magnetic gear, which consists of two radially magnetized permanent magnets.
Furthermore, there are still two levitation DoFs to be actively controlled.We used typical attraction-type electromagnetic levitation for these two DoFs.Generally, a PDcontroller is enough for ordinary electromagnetic levitation, and PD-gains are usually determined based on the plant model derived from the electromagnetic force model.However, in this case, since the attractive-type permanent magnetic bearings are involved in the work, the plant model will be affected by the permanent magnetic force.To obtain proper PD-gains, it is necessary to consider the effect caused by the permanent magnetic force.Therefore, before building the plant model, FEM analysis was conducted on the attractive-type permanent magnetic bearings to obtain their support characteristics.Then, taking the support characteristics into account, a plant model was built.With the plant model, two PD-controllers were designed.
This paper is organized as follows: Section 1 is the research background.Section 2 describes the structure and the working principle of the device.Section 3 is about the control of the two DoFs, including the plant model building process, the PD-controller design process, and simulation verification.Section 4 conducts an experiment to verify the feasibility of the device.Section 5 draws a summary and extension of this paper.

Structure Overview
Figure 1 gives a schematic diagram to highlight the advantages of this research.Figure 1a represents the traditional torsion testing mechanism, which applies contact torque to the specimen.As shown in Figure 1a, if a traditional torsion testing mechanism is used to test a specimen in a special environment (liquid, vacuum, or humid environment), since contact force is required, the whole device, including the torque sensor, actuator, etc., has to be put in a container where the special environment is created.In addition, as shown by the green line in Figure 1a, to provide support for the torque, both ends of the whole device have to be fixed to the inner wall of the container; that is to say, the container needs to bear the torque.However, on the one hand, the actuator and torque sensor contain circuits, which are sensitive to liquid, air humidity, or barometric pressure; on the other hand, since the container needs to bear the torque, the method is not applicable to the test in the soft body.As an improvement, the mechanism in Figure 1b not only employs noncontact torque using permanent magnetic gear, which prevents the actuator and torque sensor from entering the special environment, but also employs magnetic levitation to support the torque, which avoids the container bearing the torque.
trol of the two DoFs, including the plant model building process, the PD-controller process, and simulation verification.Section 4 conducts an experiment to verify the bility of the device.Section 5 draws a summary and extension of this paper.

Structure Overview
Figure 1 gives a schematic diagram to highlight the advantages of this researc ure 1a represents the traditional torsion testing mechanism, which applies contact to the specimen.As shown in Figure 1a, if a traditional torsion testing mechanism i to test a specimen in a special environment (liquid, vacuum, or humid environment) contact force is required, the whole device, including the torque sensor, actuator, et to be put in a container where the special environment is created.In addition, as by the green line in Figure 1a, to provide support for the torque, both ends of the device have to be fixed to the inner wall of the container; that is to say, the container to bear the torque.However, on the one hand, the actuator and torque sensor c circuits, which are sensitive to liquid, air humidity, or barometric pressure; on the hand, since the container needs to bear the torque, the method is not applicable to t in the soft body.As an improvement, the mechanism in Figure 1b not only employ contact torque using permanent magnetic gear, which prevents the actuator and sensor from entering the special environment, but also employs magnetic levitat support the torque, which avoids the container bearing the torque.Figure 2 gives two illustrations to describe the structure and the working pri As shown in Figure 2a,b, a magnetic levitation torsion testing device (MLTTD) m consists of a torque sensor, a servo motor, a permanent magnetic gear made of tw ally magnetized disc-shaped permanent magnets, four attractive-type permanen netic bearings each made of two axially magnetized annular permanent magnets, tw placement sensors, a specimen, two iron plates, a thrust bearing, and four electroma The specification of the main components is shown in Table 1.As shown in Figu there is a semi-transparent container; the contents of the container are the levitated which includes one of the disc-shaped permanent magnets, the thrust bearing, the men, the iron plates, and half of the attractive-type permanent magnetic bearing pa Figure 2 gives two illustrations to describe the structure and the working principle.As shown in Figure 2a,b, a magnetic levitation torsion testing device (MLTTD) mainly consists of a torque sensor, a servo motor, a permanent magnetic gear made of two radially magnetized disc-shaped permanent magnets, four attractive-type permanent magnetic bearings each made of two axially magnetized annular permanent magnets, two displacement sensors, a specimen, two iron plates, a thrust bearing, and four electromagnets.The specification of the main components is shown in Table 1.As shown in Figure 2b, there is a semi-transparent container; the contents of the container are the levitated object, which includes one of the disc-shaped permanent magnets, the thrust bearing, the specimen, the iron plates, and half of the attractive-type permanent magnetic bearing pairs.
In brief, the permanent magnetic gear is used to generate noncontact torque, and the torque sensor is used to measure the noncontact torque; the four electromagnets and the four attractive-type permanent magnetic bearings are used to control levitation.Since there are two levitation DoFs that require active control, the two displacement sensors are equipped for active control.In brief, the permanent magnetic gear is used to generate noncontact torque, a torque sensor is used to measure the noncontact torque; the four electromagnets an four attractive-type permanent magnetic bearings are used to control levitation.there are two levitation DoFs that require active control, the two displacement senso equipped for active control.

Principle of Levitation
Figure 3 presents an illustration to describe the six levitation DoFs, they are X, RY, Z, and RZ.Table 2 shows the control strategy for the six levitation DoFs.As sho Table 2, X and RY are actively controlled, while RX, Y, Z, and RZ are passively contr In more detail, referring to Figure 3 and Figure 2b, the two displacement sensors de and RY for feedback, the displacement data will be transmitted to a controller, an controller will adjust the currents of the four electromagnets in real-time according displacement data so that the levitated object will be kept at a constant position in the X-direction and the RY-direction.On the other hand, RX, Y, Z, and RZ are pas controlled by the four attractive-type permanent magnetic bearings; the principle i due to the edge effect, a magnetic force will be generated between two permanent ma of each pair once they are not concentric, and this magnetic force will hinder the p nent magnets from misalignment.

Principle of Levitation
Figure 3 presents an illustration to describe the six levitation DoFs, they are X, RX, Y, RY, Z, and RZ.Table 2 shows the control strategy for the six levitation DoFs.As shown in Table 2, X and RY are actively controlled, while RX, Y, Z, and RZ are passively controlled.In more detail, referring to Figures 3 and 2b, the two displacement sensors detect X and RY for feedback, the displacement data will be transmitted to a controller, and the controller will adjust the currents of the four electromagnets in real-time according to the displacement data so that the levitated object will be kept at a constant position in both the X-direction and the RY-direction.On the other hand, RX, Y, Z, and RZ are passively controlled by the four attractive-type permanent magnetic bearings; the principle is that, due to the edge effect, a magnetic force will be generated between two permanent magnets of each pair once they are not concentric, and this magnetic force will hinder the permanent magnets from misalignment.

Principle of Torque Application
Figure 4 presents an illustration to describe the torque application.As described before, torque is generated by the permanent magnetic gear.As shown in Figure 2b, the upper permanent of the permanent magnetic gear is driven by a servo motor so that it can rotate; once it rotates, a noncontact torque will be generated for the other permanent magnet.As shown in Figure 4, driven by this torque, the rotating body will rotate, and through the thrust bearing, the rotation will be transmitted to the specimen so that the specimen can be twisted.The four attractive-type permanent magnetic bearings

Principle of Torque Application
Figure 4 presents an illustration to describe the torque application.As described before, torque is generated by the permanent magnetic gear.As shown in Figure 2b, the upper permanent of the permanent magnetic gear is driven by a servo motor so that it can rotate; once it rotates, a noncontact torque will be generated for the other permanent magnet.As shown in Figure 4, driven by this torque, the rotating body will rotate, and through the thrust bearing, the rotation will be transmitted to the specimen so that the specimen can be twisted.

Principle of Torque Application
Figure 4 presents an illustration to describe the torque application.As described before, torque is generated by the permanent magnetic gear.As shown in Figure 2b, the upper permanent of the permanent magnetic gear is driven by a servo motor so that it can rotate; once it rotates, a noncontact torque will be generated for the other permanent magnet.As shown in Figure 4, driven by this torque, the rotating body will rotate, and through the thrust bearing, the rotation will be transmitted to the specimen so that the specimen can be twisted.Figure 5 presents a FEM analysis result, which was calculated using J-mag software.In the analysis, the airgap between the two permanents was 17 mm.The abscissas of both graphs denote the relative angle between the two permanents; the ordinate of the first graph denotes the generated torque; and the ordinate of the second graph denotes the attractive force between the two permanents in the vertical direction.In this paper, we want to obtain 0.1 N•m torque.As shown by the red area in Figure 5a, as long as the relative angle reaches no more than 20 • , 0.1 N•m torque can be obtained.On the other hand, as can be seen from the red area in Figure 5b, the attractive force will change from 18.8 N to about 17.7 N while the torque is increasing.The attractive force can be used to support the levitation in the Z-direction so that the burden of the attractive permanent magnetic bearings will be reduced and the attractive force will not be transmitted to the specimen because the thrust bearing will block any upward force.Although the attractive force will change a little during the torque application, it will be fine because the attractive-type permanent magnetic bearing will compensate for the change.
In the analysis, the airgap between the two permanents was 17 mm.The abscissas of both graphs denote the relative angle between the two permanents; the ordinate of the first graph denotes the generated torque; and the ordinate of the second graph denotes the attractive force between the two permanents in the vertical direction.In this paper, we want to obtain 0.1 N•m torque.As shown by the red area in Figure 5a, as long as the relative angle reaches no more than 20°, 0.1 N•m torque can be obtained.On the other hand, as can be seen from the red area in Figure 5b, the attractive force will change from 18.8 N to about 17.7 N while the torque is increasing.The attractive force can be used to support the levitation in the Z-direction so that the burden of the attractive permanent magnetic bearings will be reduced and the attractive force will not be transmitted to the specimen because the thrust bearing will block any upward force.Although the attractive force will change a little during the torque application, it will be fine because the attractive-type permanent magnetic bearing will compensate for the change.

Control of X and RY
As mentioned in Section 2.2, X and RY require active control; therefore, this Section will discuss this issue.

State Space Model Derivation
To start, a plant model is built as follows: Figure 6 shows two illustrations for the plant model.In Figure 6a, z1, z2, z3, and z4 are the four airgaps between the electromagnets and the iron plates.O is the centroid of the levitated object.M is the point on the surface of the iron plate at the same vertical height as the midpoint of z1 and z2.Obviously, . l is the vertical distance between the lower and upper attractive-type permanent magnetic bearings.h is the vertical distance between the lower and upper electromagnets.In Figure 6b, F1, F2, F3, and F4 denote the attractive forces of the four electromagnets on the iron plates, respectively.FP and TP, respectively denote the resultant force and resultant torque of the four attractive-type permanent magnetic bearings on the levitated object.Relative angle between the two permanent magnets (°) Relative angle between the two permanent magnets (°) Figure 5. FEM analysis results of the permanent magnetic gear: (a) Torque (b) Attractive force between the two permanent magnets in a vertical direction.

Control of X and RY
As mentioned in Section 2.2, X and RY require active control; therefore, this Section will discuss this issue.

State Space Model Derivation
To start, a plant model is built as follows: Figure 6 shows two illustrations for the plant model.In Figure 6a, z 1 , z 2 , z 3 , and z 4 are the four airgaps between the electromagnets and the iron plates.O is the centroid of the levitated object.M is the point on the surface of the iron plate at the same vertical height as the midpoint of z 1 and z 2 .Obviously, z M = z 1 +z 3 2 .l is the vertical distance between the lower and upper attractive-type permanent magnetic bearings.h is the vertical distance between the lower and upper electromagnets.In Figure 6b, F 1 , F 2 , F 3 , and F 4 denote the attractive forces of the four electromagnets on the iron plates, respectively.F P and T P , respectively denote the resultant force and resultant torque of the four attractive-type permanent magnetic bearings on the levitated object.The expressions of F1, F2, F3, and F4 are shown in Equation ( 1), where a and c are two coefficients depending on the shape and dimensions of the electromagnets and the iron plates, and i1, i2, i3, and i4 are the currents of the four electromagnets, respectively.As shown in Equation ( 2), F1, F2, F3, and F4 can be linearized at a working point (i0, z0), where The expressions of F 1 , F 2 , F 3 , and F 4 are shown in Equation ( 1), where a and c are two coefficients depending on the shape and dimensions of the electromagnets and the iron plates, and i 1 , i 2 , i 3 , and i 4 are the currents of the four electromagnets, respectively.As shown in Equation ( 2), F 1 , F 2 , F 3 , and F 4 can be linearized at a working point (i 0 , z 0 ), where i 0 and z 0 are a constant current and a constant airgap, respectively.In Equation ( 2), F 0 denotes the magnetic force when the current and airgap are just i 0 and z 0 , and ∆i in addition, k i and k z are the current coefficient and airgap coefficient, respectively.Their expressions are shown in Equation ( 3), which was derived by calculating the derivative of the magnetic force (Any of F 1 , F 2 , F 3, and F 4 ) with respect to the current (Any of i 1 , i 2 , i 3, and i 4 ) and the airgap (Any of z 1 , z 2 , z 3, and z 4 ) at the working point (i 0 , z 0 ).
It can be assumed that the whole levitated object is a rigid body, and point M has the same horizontal displacement as point M, so its dynamics can be described by Equation ( 4), where m is the mass of the levitated object and J is the moment of inertial around its centroid.F P and T P can be regarded as two generalized forces with stiffness, so Equation (4) can be rewritten as Equation (5), where k a and k r are the stiffnesses of F P and T P, respectively.
As mentioned before, z M equals z 1 +z 3 2 .Furthermore, if the variation range of θ is very small, it can be approximately regarded that θ equals z 1 −z 3 h .Therefore, Equation ( 5) can be rewritten as Equation (6).Then, substituting ∆i into Equation (6), Equation ( 7) can be obtained.
After that, define a state space as follows: where X is the state variable, U is the input, Y is the output, and Equation ( 13) is the state space expression.From Equation ( 9), the expressions of A, B and C can be drawn as follows: So far, the plant model has been built.

Determination of Plant Model Parameter
In Equations ( 14)-( 16), there are some parameters to be determined; they are m, J, h, l, k i , k z , k r , and k a .Among them, m, J, h, and l can be easily calculated based on the densities, dimensions, and shapes of the solids; their values are 1.731 kg, 3.3 × 10 −3 kg•m 2 , 0.143 m, and 0.0774 m, respectively.k i and k z can be calculated using Equation (3), however, the calculation requires the values of a and c.Therefore, the attractive force of a single electromagnet to an iron plate was calculated at various combinations of current and airgap by FEM analysis.After the analysis, the result data were imported into the Curve Fitting Toolbox of MATLAB to obtain a fitting surface, which was generated in the form of F = ai 2 (z+c) 2 , which is the general equation of electromagnetic force.The generated fitting surface is shown in Figure 7, where the black point is the original data.With the generation of the fitting surface, a and c were also calculated out by the toolbox; their values were 1.204 × 10 −5 Nm 2 A −2 , 2.38 × 10 −3 m, respectively.Furthermore, i 0 and z 0 were set as 2 A and 3 mm, and then k i and k z were calculated using Equation ( 3); their values were 1.6639 N/A and 618.5426N/m, respectively.l, ki, kz, kr, and ka.Among them, m, J, h, and l can be easily calculated based on the densities dimensions, and shapes of the solids; their values are 1.731 kg, 3.3 × 10 −3 kg•m 2 , 0.143 m and 0.0774 m, respectively.ki and kz can be calculated using Equation ( 3), however, the cal culation requires the values of a and c.Therefore, the attractive force of a single electromag net to an iron plate was calculated at various combinations of current and airgap by FEM analysis.After the analysis, the result data were imported into the Curve Fitting Toolbox o MATLAB to obtain a fitting surface, which was generated in the form of F= ai 2 (z+c) 2 , which is the general equation of electromagnetic force.The generated fitting surface is shown in Fig ure 7, where the black point is the original data.With the generation of the fitting surface, a and c were also calculated out by the toolbox; their values were 1.204 × 10 −5 Nm 2 A −2 , 2.38 × 10 −3 m, respectively.Furthermore, i0 and z0 were set as 2 A and 3 mm, and then ki and kz were calculated using Equation ( 3); their values were 1.6639 N/A and 618.5426N/m, respectively Furthermore, to determine values for kr and ka, another FEM analysis was conducted on the four attractive-type permanent magnetic bearings.In the FEM analysis, the four permanent magnets in the bearing holder were set to be static, and the four permanen magnets embedded in the iron plate were set to be movable; then, FP and TP are calculated under various lateral displacements (zM-z0) and rotation angles (θ) of the four permanen magnets embedded in the iron plate.The results are shown in Figure 8. Figure 8a shows Furthermore, to determine values for k r and k a , another FEM analysis was conducted on the four attractive-type permanent magnetic bearings.In the FEM analysis, the four permanent magnets in the bearing holder were set to be static, and the four permanent magnets embedded in the iron plate were set to be movable; then, F P and T P are calculated under various lateral displacements (z M -z 0 ) and rotation angles (θ) of the four permanent magnets embedded in the iron plate.The results are shown in Figure 8. Figure 8a shows the result of F P , and Figure 8b shows the result of T P .It can be drawn from the slope of the red dashed line in both graphs that the stiffnesses k a and k r are approximately +12.5 N/mm and +0.133 N•m/ • in the vicinity of the initial point (abscissa = 0).So far, all the plant model parameters have been caught, and the main parameter values are shown in Table 3.  3.

Control Design
As can be seen from Equation ( 14), the first two state variables and the last two state variables are not coupled to each other, therefore they can be controlled individually, so a feedback matrix with the following structure can be used: where P 1 , D 1 , P 2 , and D 2 can be regarded as two PD-gains of two independent PDcontrollers.If the feedback matrix is used, a closed-loop system matrix will be as follows.
According to Equation (18), two independent characteristic equations can be derived as follows: It can be seen from Equation (19) that they are two second-order systems.Generally, second-order systems have characteristic equations in the form of Equation (20), where ζ 1 and ζ 2 denote the damping ratios of the systems, ω n1 and ω n2 denote the undamped oscillation angular frequencies of the systems.Comparing Equations ( 19) and ( 20), the following equations can be derived.
A damping ratio of 0.707 can coordinate overshoot and settling time; thus, it is usually regarded as an optimal damping ratio in practice [30].Therefore, in this paper, 0.707 is taken for both ζ 1 and ζ 2 .As a result, the overshoots will approach 4.33%.In addition, 0.025 s is taken for both t s1 and t s2 , to realize that ω n1 and ω n2 should be 169.73s −1 .Then, substituting ζ 1 = ζ 2 = 0.707 and ω n1 = ω n2 = 169.73s −1 into Equations ( 23) and ( 24), PD-gains that realize the desired dynamics performance will be obtained as Table 4.

Control Simulation
To verify the feasibility of the designed PD-controllers, control simulation was conducted using Matlab-Simulink.Control diagram is shown in Figure 9.The plant model parameters and control parameters in the simulation are the same as Tables 3 and 4. To make the simulation closer to reality, the electromagnetic force model used in the simulation was the nonlinear plant model denoted by Equation (1) rather than the linearized plant model denoted by Equation (2).Moreover, to test the dynamic performance of the system, two sets of 0.1 mm step disturbances were applied to the airgaps at 0.1 s and 0.5 s, respectively; the first set was in the same direction, and the second set was in the opposite direction.The simulation result is shown in Figure 10. ( ) A damping ratio of 0.707 can coordinate overshoot and settling time; thus, it is usually regarded as an optimal damping ratio in practice [30].Therefore, in this paper, 0.707 is taken for both ζ1 and ζ2.As a result, the overshoots will approach 4.33%.In addition, 0.025 s is taken for both ts1 and ts2, to realize that ωn1 and ωn2 should be 169.73s −1 .Then, substituting ζ1=ζ2 = 0.707 and ωn1=ωn2 = 169.73s −1 into Equations ( 23) and (24), PD-gains that realize the desired dynamics performance will be obtained as Table 4.

Control Simulation
To verify the feasibility of the designed PD-controllers, control simulation was conducted using Matlab-Simulink.Control diagram is shown in Figure 9.The plant model parameters and control parameters in the simulation are the same as Tables 3 and 4. To make the simulation closer to reality, the electromagnetic force model used in the simulation was the nonlinear plant model denoted by Equation (1) rather than the linearized plant model denoted by Equation (2).Moreover, to test the dynamic performance of the system, two sets of 0.1 mm step disturbances were applied to the airgaps at 0.1 s and 0.5 s, respectively; the first set was in the same direction, and the second set was in the opposite direction.The simulation result is shown in Figure 10.As shown in Figure 10, after bearing the first step disturbance, the overshoots airgaps were about 3%, and the settling times (5% tolerance) of both airgaps wer 0.02 s.On the other hand, after bearing the second step disturbance, the overshoots airgaps were about 3%, and the settling times (5% tolerance) of both airgaps wer 0.02 s.This result implies that the dynamic performances in the simulation are different from the designed dynamic performance (overshoot: 4.33%; settling tim s); this may be due to the nonlinearity of the plant model used.Anyway, the ove (about 3%) of both responses were lower than the designed overshoot (4.33%), wh plies the stability of the system is better than the design standard.On the other ha settling times (about 0.02 s) of both responses were shorter than the designed settli (0.025 s), which implies the rapidity of the system is better than the design stand sum up, the designed PD-controllers are feasible.

Experiment
To validate the feasibility of the devices, an experiment was conducted.Fig presents a photograph of the experiment setup, and Figure 11b presents an overall structure diagram.The model of the torque sensor is FT05, made by FORSENTE Limited(Shenzhen, China); its measuring range is ±0.2 N•m.The model of the d ment sensors is GP-XC12ML, and the manufacturer is Panasonic Industry C (Osaka, Japan).The power amplifiers are ESCON made by MAXON Japan Corp. ( Japan).The control system is MicroLabBox dSPACE; it is equipped with softwa ControlDesk.In ControlDesk, we can import a Matlab-Simulink file, adjust con rameters, record, and observe experiment data in real-time.In addition, as can be Figure 11b, there are two paths, i.e., the levitation path and the torsion path; the fo used to stabilize X and RY, and the latter is used to apply noncontact torque.As men before, the noncontact torque is generated by the relative rotation between the tw shaped permanent magnets; the rotation is driven by the servo motor.In the expe the servo motor drove the upper disc-shaped permanent magnet to rotate at a c speed, so that a noncontact torque was generated on the lower disc-shaped per magnet, and through the thrust bearing, the torque was finally transmitted to th 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 As shown in Figure 10, after bearing the first step disturbance, the overshoots of both airgaps were about 3%, and the settling times (5% tolerance) of both airgaps were about 0.02 s.On the other hand, after bearing the second step disturbance, the overshoots of both airgaps were about 3%, and the settling times (5% tolerance) of both airgaps were about 0.02 s.This result implies that the dynamic performances in the simulation are slightly different from the designed dynamic performance (overshoot: 4.33%; settling time: 0.025 s); this may be due to the nonlinearity of the plant model used.Anyway, the overshoots (about 3%) of both responses were lower than the designed overshoot (4.33%), which implies the stability of the system is better than the design standard.On the other hand, the settling times (about 0.02 s) of both responses were shorter than the designed settling time (0.025 s), which implies the rapidity of the system is better than the design standard.To sum up, the designed PD-controllers are feasible.

Experiment
To validate the feasibility of the devices, an experiment was conducted.Figure 11a presents a photograph of the experiment setup, and Figure 11b presents an overall system structure diagram.The model of the torque sensor is FT05, made by FORSENTEK Co., Limited(Shenzhen, China); its measuring range is ±0.2 N•m.The model of the displacement sensors is GP-XC12ML, and the manufacturer is Panasonic Industry Co., Ltd.(Osaka, Japan).The power amplifiers are ESCON made by MAXON Japan Corp. (Tokyo, Japan).The control system is MicroLabBox dSPACE; it is equipped with software, i.e., ControlDesk.In ControlDesk, we can import a Matlab-Simulink file, adjust control parameters, record, and observe experiment data in real-time.In addition, as can be seen in Figure 11b, there are two paths, i.e., the levitation path and the torsion path; the former is used to stabilize X and RY, and the latter is used to apply noncontact torque.As mentioned before, the noncontact torque is generated by the relative rotation between the two disc-shaped permanent magnets; the rotation is driven by the servo motor.In the experiment, the servo motor drove the upper disc-shaped permanent magnet to rotate at a constant speed, so that a noncontact torque was generated on the lower disc-shaped permanent magnet, and through the thrust bearing, the torque was finally transmitted to the specimen.Figure 12 shows the experiment result, where Figure 12a shows torque data measured by the torque sensor and Figure 12b,c show airgap response and current response, respectively.As shown in Figure 12a, the torque kept increasing throughout the whole process; the maximum torque was about 0.126 N•m.As shown in Figure 12b, the two airgaps slightly fluctuated, but the airgap errors were within ±0.05 mm, which is considered fine.Figure 12 shows the experiment result, where Figure 12a shows torque data measured by the torque sensor and Figure 12b,c show airgap response and current response, respectively.As shown in Figure 12a, the torque kept increasing throughout the whole process; the maximum torque was about 0.126 N•m.As shown in Figure 12b, the two airgaps slightly fluctuated, but the airgap errors were within ±0.05 mm, which is considered fine.

Conclusions
To facilitate torsion testing in special environments, this paper proposes a noncontact torsion testing device using magnetic levitation technology, where a noncontact torque can be transmitted to a specimen.To ensure the stability of the six levitation degrees of freedom (DOF), four levitation DoFs are passively controlled by some permanent magnets, while two levitation DoFs are actively controlled by four electromagnets.Furthermore, a plant model was built for the control of the two levitation DoFs.With the plant model, two PD-controllers were designed, and the simulation result demonstrated that the designed PD-controllers make the system have good dynamic performance.Finally, an experiment was performed to test the stability of the levitation and the torque load capacity.The experiment result indicated that the device could apply a 0.126 N•m noncontact torque to the specimen while maintaining levitation.
Currently, the control mode of the torque is open-loop, which makes it impossible to accurately control the amount of torque.In the future, a torque feedback system will be

Conclusions
To facilitate torsion testing in special environments, this paper proposes a noncontact torsion testing device using magnetic levitation technology, where a noncontact torque can be transmitted to a specimen.To ensure the stability of the six levitation degrees of freedom (DOF), four levitation DoFs are passively controlled by some permanent magnets, while two levitation DoFs are actively controlled by four electromagnets.Furthermore, a plant model was built for the control of the two levitation DoFs.With the plant model, two PD-controllers were designed, and the simulation result demonstrated that the designed PD-controllers make the system have good dynamic performance.Finally, an experiment was performed to test the stability of the levitation and the torque load capacity.The experiment result indicated that the device could apply a 0.126 N•m noncontact torque to the specimen while maintaining levitation.
Currently, the control mode of the torque is open-loop, which makes it impossible to accurately control the amount of torque.In the future, a torque feedback system will be constructed using torque sensors to apply the desired amount of torque.

Figure 2 .
Figure 2. Physical structure of the noncontact torsion testing device: (a) Overall structure (b) exploded view.

Figure 2 .
Figure 2. Physical structure of the noncontact torsion testing device: (a) Overall structure (b) Partial exploded view.

Figure 3 .
Figure 3. Illustration of the six levitation DoFs.Table 2. Control strategy of the six levitation DoFs.Degree of Freedom Control Type Component That Implements Control X Active The four electromagnets Y Passive The four attractive-type permanent magnetic bearings Z Passive The four attractive-type permanent magnetic bearings RX Passive The four attractive-type permanent magnetic bearings RY Active The four electromagnets RZ Passive The four attractive-type permanent magnetic bearings

Figure 4 .
Figure 4. Illustration of the torque application.

Figure 4 .
Figure 4. Illustration of the torque application.

Figure 4 .
Figure 4. Illustration of the torque application.

Figure 5 .
Figure 5. FEM analysis results of the permanent magnetic gear: (a) Torque (b) Attractive force between the two permanent magnets in a vertical direction.

Figure 6 .
Figure 6.Illustration of the plant model: (a) Illustration of the dimensions (b) Illustration of the forces.

Figure 6 .
Figure 6.Illustration of the plant model: (a) Illustration of the dimensions (b) Illustration of the forces.

Figure 7 .
Figure 7. Fitting surface generated from the FEM analysis result data.

Figure 7 .
Figure 7. Fitting surface generated from the FEM analysis result data.

x
FOR PEER REVIEW 10 of 17 the result of FP, and Figure 8b shows the result of TP.It can be drawn from the slope of the red dashed line in both graphs that the stiffnesses ka and kr are approximately +12.5 N/mm and +0.133 N•m/° in the vicinity of the initial point (abscissa = 0).So far, all the plant model parameters have been caught, and the main parameter values are shown in Table

Figure 8 .
Figure 8. FEM analysis results of the four attractive-type permanent magnetic bearings: (a) Analysis result of FP (b) Analysis result of TP.

Figure 8 .
Figure 8. FEM analysis results of the four attractive-type permanent magnetic bearings: (a) Analysis result of F P (b) Analysis result of T P .

Table 1 .
Specification of main components.

Table 1 .
Specification of main components.

Table 2 .
Control strategy of the six levitation DoFs.

Table 3 .
Main parameters of the plant model.

Table 3 .
Main parameters of the plant model.

Table 5 .
Control parameters in the experiment.

Table 6 .
Specification of the specimen.

Table 6 .
Specification of the specimen.