Next Article in Journal
Reproducible Molecularly Imprinted Piezoelectric Sensor for Accurate and Sensitive Detection of Ractopamine in Swine and Feed Products
Previous Article in Journal
Optical Tracking Data Validation and Orbit Estimation for Sparse Observations of Satellites by the OWL-Net
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Vibration Suppression Control of a Modular Elastic Joint

State Key Laboratory of Robotics and System, Harbin Institute of Technology, West Dazhi Street, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Submission received: 16 April 2018 / Revised: 30 May 2018 / Accepted: 4 June 2018 / Published: 7 June 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, a novel mechatronic design philosophy is introduced to develop a compact modular rotary elastic joint for a humanoid manipulator. The designed elastic joint is mainly composed of a brushless direct current (DC) motor, harmonic reducer, customized torsional spring, and fail-safe brake. The customized spring considerably reduces the volume of the elastic joint and facilitates the construction of a humanoid manipulator which employs this joint. The large central hole along the joint axis brings convenience for cabling and the fail-safe brake can guarantee safety when the power is off. In order to reduce the computational burden on the central controller and simplify system maintenance, an expandable electrical system, which has a double-layer control structure, is introduced. Furthermore, a robust position controller for the elastic joint is proposed and interpreted in detail. Vibration of the elastic joint is suppressed by means of resonance ratio control (RRC). In this method, the ratio between the resonant and anti-resonant frequency can be arbitrarily designated according to the feedback of the nominal spring torsion. Instead of using an expensive torque sensor, the spring torque can be obtained by calculating the product of spring stiffness and deformation, due to the high linearity of the customized spring. In addition, to improve the system robustness, a motor-side disturbance observer (DOb) and an arm-side DOb are employed to estimate and compensate for external disturbances and system uncertainties, such as model variation, friction, and unknown external load. Validity of the DOb-based RRC is demonstrated in the simulation results. Experimental results show the performance of the modular elastic joint and the viability of the proposed controller further.

1. Introduction

Robots, especially those working in industrial environments, have been generally designed according to the principle “rigidity by design, safety by sensors and control”. The safety to interact with rigid environments is usually guaranteed by means of active compliance control, i.e., position/force or impedance/admittance control [1,2,3]. However, the delays existing in the software system make the rigid actuator unsafe during sudden impacts and is incapable of handling contact transients. To overcome these difficulties, series elastic actuators (SEA) are proposed as a simple and effective solution to realize inherent compliance [4]. Unlike a rigid actuator, a SEA contains an elastic element in series, which introduces low mechanical output impedance, passive mechanical energy storage, and increases peak power output. Moreover, the inherent compliance brought by the elastic element can filter impact, protect drivetrains, and provide additional time for active controllers to regulate the impedance of the actuator.
SEAs have been successfully applied in numerous scenarios for more than 20 years. A prismatic SEA, using a ball-screw mechanism for speed reduction, was designed in reference [5]. This mechanism is highly efficient and has a high tolerance for impact loads. However, the main drawback of using a ball-screw mechanism in a SEA is that an extra output mechanism is needed for converting translational motion to rotation. In reference [6], a new compact soft actuation unit intended to be used in multi-DOF and small scale robotic systems was introduced. Although the designed actuator shows high integration density and wide passive compliance, the torsional spring designed based on an arrangement of linear springs is still too thick to increase the size of the actuator and suffers from nonlinear stiffness. Cummings et al. [7] designed a compact modular rotary SEA, which can meet the requirements of a wide range of applications. A customized torsional spring designed for the NASA Robonaut 2 platform was used in this SEA. However, the employment of planetary gearhead made the SEA unsuitable for a manipulator, despite a novel sensor geometry being presented for mitigating the negative impact of backlash. Moreover, SEAs are also widely employed in bipedal walking robots [8], robotic devices for rehabilitation [9,10], and wearable robotics for assistance [11]. Those SEAs also suffer from large size, non-modularity, and mechanical design complexity. Recently, variable stiffness actuators (VSAs) [12,13] have been introduced to maximize benefits of SEA with respect to safety during interaction, as well as to improve the energy efficiency of the robotic system. VSAs typically employ two motors, which respectively control the compliance and the equilibrium position of the actuated joint. The introduction of an extra motor to achieve the functionality of the stiffness adjustment may lead to the increase of complexity, size, and the total weight of the actuation system.
Another problem is that introduction of the elastic element reduces the actuator bandwidth and makes position control difficult [14]. Moreover, it is prone to vibrate, which deteriorates system stability and requires extra power for position control [15]. Several approaches have been proposed to control the position of elastic-joint robots. Tomei [16] proved that a simple proportional-differential (PD) controller can globally stabilize about a reference position. Kim et al. [17] proposed a robust PD control scheme for flexible-joint robots based on a DOb which was only applied to the motor-side dynamics. Jin et al. [18] designed an adaptive tracking controller using a time-delay estimation technique. While his approach was to achieve accurate motor-side position before the elastic element and let the arm side show passive compliant behaviors. Additionally, many vibration suppression approaches for elastic joints or elastic-joint robots have been proposed and explicit analyses have been presented. As an effective vibration suppression control method, joint damping controllers have been reported in the literature. Well-known approaches such as feedback linearization [19], model-based state-feedback controllers [20], learning control schemes [21], or linear-quadratic regulators [22], have been developed. These schemes are highly effective and show great performance. But, a limitation of these approaches is their dependence on an accurate mathematical model. In order to solve this problem, Petit et al. [23] proposed a model-free damping control approach, which utilized the joint elasticities to convert kinetic energy into elastic energy to achieve vibration suppression. However, this control approach was employed mainly to deal with the vibration introduced by disturbances. Besides, a number of feedforward control techniques have been applied to robots with elastic joints. In reference [24], a command shaping method has been used for flexible joint to reduce vibration. The same goal has been achieved for a two-inertia system [25], in which a robust control scheme combining backstepping technique and partial DOb has been proposed. This controller also included a reference generator, which provided a reference state and a feedforward input signal. Although it plays an important role in vibration suppression, feedforward control reaches a compromise between rapidity and shaper robustness for positive input shapers and smoothly shaped reference commands [26]. Furthermore, controllers of a prismatic SEA for force and position tracking based on combinations of PID, model-based, and DOb structures were presented in reference [5]. It is noted that actuator resonance is suppressed not only by the DOb but also by the joint structure in which the spring is placed between the motor housing and the chassis ground.
In this paper, a novel compact, modular rotary elastic joint is designed in detail. The designed elastic joint is evolved from the traditional rigid joint in which the expensive torque sensor is replaced by a torsional spring. This design philosophy greatly reduces the design cycle and cost. As to the mechanical system, the joint adopts a new customized torsional spring as an elastic element. The customized torsional spring is optimized using Finite Element Method (FEM) simulations to satisfy the specific performance requirements. A large central hole along the joint axis is introduced to ease cabling. In order to maintain position when an outage happens, a fail-safe brake is designed. In the electrical system, a double-layer control structure, comprising a central controller and a joint controller, is employed to reduce the computational burden and simplify system maintenance. Besides, the double-layer control structure of electrical system is suitable for constructing a humanoid manipulator in future. The central controller is mainly responsible for the implementation of a high-level control algorithm, and joint controller is responsible for sensor information processing, motor and fail-safe brake control. Data communication between central controller and joint controller is achieved by a Point-to-Point High Speed Serial Communication (PPSeCo) protocol. All the configurations above make the designed elastic joint light weight, small volume, and suitable as a module for a humanoid manipulator. In addition, a high-performance robust position controller is proposed for the elastic joint. Vibration is suppressed by using RRC, of which the control gains are determined theoretically. In order to improve the system robustness, a motor-side DOb and arm-side DOb are employed to deal with external disturbances and system uncertainties, such as model variation, friction, and load. Simulation and experimental results are given to show the validity of the DOb-based RRC.
This paper highlights research in design and vibration suppression control of a modular elastic joint. Our contributions include: (1) A mechanical design of the elastic joint which is more compact and lightweight than previous SEA designs; (2) a design of the non-uniform spiral spring employing only one profile variable parameter, which greatly reduces the complexity of optimization process; and (3) precise point-to-point position control with which vibration is effectively suppressed. This paper is organized as follows. In Section 2, the mechanical and electrical designs of the elastic joint are elaborated. Section 3 proposed a DOb-based RRC with a detailed explanation and analysis. In Section 4, position control simulations and experimental results are given to validate the efficacy of the proposed control scheme. Section 5 presents the conclusions.

2. Elastic Joint Design

For a humanoid manipulator, a modular design can strengthen interchangeability of the parts and eliminate structural differences between the joints. In terms of the configuration and workload, a humanoid manipulator can be established by two types of modular joints, namely shoulder/elbow joint and wrist joint. They are identical in the mechanical structure and the electrical system, except for the selection of some standard components. In this paper, an elastic joint in shoulder/elbow is designed. Design indexes and actual values are listed in Table 1.

2.1. Mechanical Design

The ideal objective of the elastic joint design is to achieve a compact and modular mechanical architecture. In order to reduce design circle and cost, the designed elastic joint is evolved from the traditional rigid joint, which replaces the expensive torque sensor with the torsional spring. The design diagram of the elastic joint is shown in Figure 1b. The final assembly of the elastic joint consists of a brushless DC motor, harmonic reducer, customized torsional spring, bearings, magnetic encoders, fail-safe brake, etc. The motor stator and circular spline of the harmonic reducer are fixed on the joint base. The motor rotor drives the wave generator through the motor shaft and the flexspline of the harmonic reducer drives the joint output through the torsional spring. The rotation disk and shell of the fail-safe brake are fixed on the motor shaft and the shell of the joint respectively. All housing parts are made of aluminum alloy 7A09 for the purpose of weight reduction.
  • Motor components: A brushless DC motor 72-L designed by Harbin Institute of Technology (HIT) is selected. The rotor has large center hole and the maximum diameter of the stator is 72 mm. Its rated torque and maximum torque are respectively 0.71 Nm and 0.85 Nm respectively.
  • Reducer components: A harmonic reducer CSD-25 of Harmonic Drive Company is selected. Its single-stage reduction ratio is 160 and the corresponding peak torque is 123 Nm.
  • Customized torsional spring: To minimize the volume of the joint, a customized torsional spring is designed. Compared to the existing torsional springs, this spring has the characteristics of small volume and high linearity between applied torque and deformation. See Section 2.3 for more details.
  • Bearings: A deep groove ball bearing is used to guarantee that the two rings of the torsional spring rotate coaxially. Two deep groove ball bearings are used to support the motor shaft. By replacing a pair of traditional angular contact ball bearings, a cross roller bearing THK-RA8008C, which supports the joint output, can make the joint more compact.
  • Magnetic encoders: All position information is measured by three 16-bit absolute magnetic rotary encoders designed by HIT. Motor position, motor position after reduction, and spring deformation β can be measured respectively. The motor-side position θ , used in the control system is calculated by means of sensor data fusion [27] between motor position and motor position after reduction. The arm-side position q is calculated by the difference between the motor-side position and the spring deformation, i.e., q = θ β . The magnetic encoder consists of two inductors, i.e., a code disk and a signal processing part. It possesses small-size, light, robust, and easy-to-integrate properties.
  • Fail-safe brake: In order to maintain current position when an outage happens, the elastic joint integrates with the fail-safe brake, which is mainly composed of a rotation disk, wear-resisting friction disks, a pressing disk, an electromagnet, springs, and a shell. Eight wear-resisting friction disks are evenly distributed on the rotation disk. The force generated by the electromagnet separates the pressing disk against the spring force from the friction disks when the power is on. When the power is off, the spring force presses the pressing disk on the friction disks, which makes the motor shaft stop. The fail-safe brake designed in this paper has the advantages of low energy consumption and fast braking.
Employments of standard (CSD-25, THK-RA8008C) and customized components (brushless DC motor, magnetic encoders, torsional spring) make the designed elastic joint more compact. The large central hole diameter of 17 mm along the joint axis brings convenience for cabling. The integrated fail-safe brake can stop the joint and guarantee safety when the power is off. The designed elastic joint only uses three 16-bits absolute magnetic rotary encoders for sensor information, rather than expensive sensors such as a torque sensor. Compared to the existing elastic joints, the designed elastic joint has the characteristics of light-weight, small volume and high energy, which is suitable for a module for a humanoid manipulator. Further, the designed joint is able to achieve full rotation and has little hysteresis.

2.2. Electrical System Architecture

In this section, an expandable electrical system architecture is proposed as shown in Figure 2. The dashed part represents other joints, which can together constitute a humanoid manipulator. The electrical system adopts double-layer control structure, which mainly consists of a central controller and a joint controller. Some key factors of this scheme are taken into account. First of all, the double-layer control structure can reduce the computational burden and simplify system maintenance compared to a single controller. Another reason for the choice of double-layer control structure is that it facilitates building a humanoid manipulator in the future.
The main tasks performed by the central controller are motion planning of elastic joint or manipulator, implementation of control algorithm, and providing interface with host computer. The joint controller is responsible for sensor information processing, motor and fail-safe brake control, etc. In order to increase the transmission reliability and speed of data communication, reduce cabling and noise in sensor signals, a PPSeCo protocol proposed by our lab is employed here for the data communication between central controller and joint controller [28]. All communication programs are written in VHDL and run in FPGAs. The frequency of system is 500 Hz determined by the software interrupt.
The hardware of the central controller is a Compact PCI-based DSP/FPGA board. The main processor of the central controller is a TI floating-point DSP TMS320C6713 with a maximum 1,350 MFLOPS. The DSP is an excellent choice to easily realize complex control algorithm and fast computation because of its high performance floating-point capabilities. Furthermore, the DSP provides a series of peripherals, which make it easy for user to access and extend the hardware resources. Serial signals from the joint controller are converted to parallel signals and transmitted to the DSP via the parallel interface by the FPGA, and vice versa. The DSP/FPGA board and host computer are packed into a cabinet based on a Compact PCI structure, as well as a DC/DC converter which provides power for joint motor, sensors, and fail-safe brake.
Joint controller is composed of a driver board, a control board, and a connecting board. The control board adopts the Cyclone IV series FPGA produced by Altera Company (San Jose, CA, USA) as the microprocessor to realize the control of the motor and the fail-safe brake, sensor information processing, and high-speed serial communications with the central controller. The drive board receives the commands from the control board to drive the motor and the fail-safe brake. In order to reduce the modular joint’s weight and volume, the joint controller is placed outside the joint and is intended to be placed inside the link in the future.

2.3. Spring Design

Being a fundamental component of the SEA, the elastic element has great influence on the size and mechanical properties of the joint. Thus, the design of the torsional spring is the key for the rotary SEA. Generally, there are two forms of torsional spring design in the literature. The first form [6,29] is usually based on an arrangement of linear springs, which leads to a larger volume. The second form adopts a customized design to minimize the volume and improve performance of the torsional spring. A systematic approach to the design of a monolithic torsional elastic module has been presented in reference [30]. The designed elastic module shows a good linearity characteristic between torque and angle. By arranging a number of identical springs in series or in parallel, it is possible to render different torque vs. angle characteristics, in order to match the special application requirements. Also employing the customized form, various torsional springs [7,10,31] have been designed and used in SEAs. Those torsional springs have the advantages of light-weight and small volume. In the design of the customized torsional spring, the first step is choosing a set of topologies. Commonly there are three topologies of the customized the torsional spring:
  • Groups of symmetrical flexible beams are placed between the inner and outer ring of the spring. Due to the limited space, flexible beams in such springs are usually thin and thus can only withstand relatively low torque [31].
  • Uniform spiral beams. The typical representation is an Archimedean spiral torsional spring. Although the spring employing this topology has a good ability to withstand high torque, hysteresis and backlash in its structure may deteriorate the dynamic performance [10].
  • Non-uniform spiral beams. The spring employing this topology can withstand high torque. However, the profile of the non-uniform spiral beams, which is depicted by a multi-segment curve is difficult to construct [7]. Furthermore, this kind of spring has more profile parameters, which leads to greater complexity of the optimization.
In order to realize small volume and high-performance index of the output torque, a customized torsional spring which has non-uniform spiral beams is designed in this paper. Different from the previous torsional spring, the spring designed in this paper employs eccentric-arcs to build non-uniform spiral beams, which requires only one profile variable parameter. Figure 3a shows a schematic perspective view of the customized torsional spring. It is composed of two rings interconnected by a pair of non-uniform spiral beams, which are outlined by eccentric-arcs R5 and R6 with different radii. R5 and R6 are equal to 50 mm and 58 mm, respectively. The arc lengths of R5 and R6 are about 3π/2. The values of R1, R2, R3, and R4 are constant (70 mm, 83 mm, 32 mm, and 17 mm) and constrained by the elastic joint dimensions. To simplify the design, the values of L1 and L2 are set to be constant. The corresponding design variables are, respectively, vertical eccentric distance L and spring thickness T. The thickness of the flexible beam changes with different L.
Flexible beams generate elastic deformation to transfer the torque and absorb the impact when two rings rotate relatively. Spring steel 50CrVA is selected as the spring material of which the Young’s modulus and yield strength are 205 GPa and 1,320 MPa, respectively. The stress distribution and deformation of the customized torsional spring was analyzed by FEM to ensure that the maximum stress is less than the yield strength of the material. The theoretical spring stiffness is computed by the ratio between the applied torque and the corresponding angular deformation obtained in the simulation.
The spring is responsible for almost all the passive compliance of the elastic joint and affects the dynamics of the elastic joint, the torque sensing range and resolution. Thus, proper choice of the spring stiffness is critical. However, choosing the optimal spring stiffness K s d e s i r e for an elastic joint of humanoid manipulator is still under investigation and controversial discussion [32]. Some indications lead to a value of about 350 Nm/rad in the elbow [33]. Then, the maximum angular deformation of the joint β max can be calculated as follows.
β max = τ max K s d e s i r e = 60   Nm 350   Nm / rad = 0.171   rad ,
where τ max is the max elastic output torque of the elastic joint. The methodology to determine the geometry parameters of the spring is presented as follows. The value of distance L is varied from 3.8 mm to 5.2 mm with a step of 0.1 mm for each thickness T from 9.5 mm to 12 mm with a step of 0.5 mm. Figure 4 shows the results of the adopted methodology. It can be concluded that spring becomes more stiff and robust when T is thicker. From Figure 4a, it can be observed that Von Mises stress roughly decreases as the variable L increases. As shown in Figure 4b, the relationship between variable L and angular deformation at different thicknesses of the spring can be characterized by a linear approximation.
When the geometry parameters L and T are, respectively, 4.5 mm and 10.5 mm, the spring stiffness obtained in the simulation is approximately 350 Nm/rad and the torsional spring is subjected to minimal stress. Correspondingly, the maximum value of Von Mises stress obtained in the simulation is 1261.9 MPa as shown in Figure 5. The customized torsional spring, shown in Figure 3b, is manufactured by using the Wire Electrical Discharge Machining (WEDM) process.
Actual stiffness of the designed spring should be calibrated. In this paper, it can be achieved by measuring joint stiffness by virtue of loading experiments. Loading and unloading are implemented on the flange which is mounted on the joint output. The torque exerted on the joint output can be calculated. Figure 6 shows that the calculated torque and the corresponding angular deformation of the elastic joint measured by the magnetic encoder 1. As shown in Figure 6, there is little hysteresis with loading and unloading on the joint. A linear regression was performed in both directions and the corresponding value of the joint stiffness can be calculated as 464.2 Nm/rad. The value of the joint stiffness, experimentally determined, is approximately 32.6% larger than the spring stiffness obtained by finite element analysis. This discrepancy probably arises in the imperfection of model used in the analysis and the different properties between the actual material after thermal treatment and the nominal counterparts. Besides, measurement errors of the sensor, the load and the arm may have a great influence on the experimental results.

3. Position Controller Design

3.1. Dynamic Model

As shown in Figure 1, the motor side is connected to the arm side with the torsional spring of which its elasticity easily leads to resonance in the system. The dynamic equation of the elastic joint can be described by:
M θ ¨ + f m + τ s = τ u ;
J q ¨ + f l + g ( q ) = τ s ;
τ s = K s ( θ q ) ,
where M , J are the inertias of motor side and arm side. f m and f l are viscous and coulomb frictions of motor side and arm side respectively. τ u is the motor-side input torque and τ s is the spring torque due to the joint deformation specified in Equation (4). g ( q ) is the gravitational torque and K s denotes the joint stiffness.
In practice, real model parameters cannot be achieved, i.e., we can only obtain nominal model parameters. So, the dynamic model of the elastic joint can be reformulated as follows:
M n θ ¨ + D m + τ s n = τ u ;
J n q ¨ + D l = τ s n ,
τ s n = K s n ( θ q ) ,
where M n and J n are nominal inertia of motor side and arm side respectively. D m and D l denote lumped terms including all nonlinearities, uncertainties, and external disturbances, which can be defined as:
D m = ( M M n ) θ ¨ + f m + ( τ s τ s n ) ;
D l = ( J J n ) q ¨ + f l + g ( q ) ( τ s τ s n ) ,
when D m and D l are neglected, joint dynamic equation can be simplified into two-mass resonant system and the block diagram is shown in Figure 7.
From Figure 7, the transfer function from the motor-side torque τ u to the motor-side position θ and the transfer function from the motor-side torque τ u to the arm-side position q can be expressed as follows:
θ τ u = 1 M n ( s 2 + K s n J n ) s 4 + ( K s n J n + K s n M n ) s 2 ;
q τ u = K s n M n J n s 4 + ( K s n J n + K s n M n ) s 2 .
Then, the anti-resonant frequency ω a r in the arm side, which represents the natural frequency of the elastic joint and the resonant frequency, ω r , in the motor side can be described as follows:
ω a r = K s n J n ;
ω r = K s n J n + K s n M n .
The characteristics of the two-mass resonant system can be determined if model parameters K s n , M n and J n are known. Joint stiffness ( K s n = 464.2   Nm / rad ) can be obtained by the experiment stated in Section 2.3. The inertia of motor side M n can be calculated by three-dimensional model of the elastic joint established in PRO/ENGINEER. The inertia of motor side is mainly determined by the rotor, motor shaft. and wave generator. Multiplying the sum of these inertia values by square of the reduction ratio N, the inertia of motor side can be obtained, namely, M n = M total N 2 = 1.567   Kg m 2 . According to Equation (12), the motor-side inertia J n is equal to K s n / ω a r 2 . Due to the low arm-side friction, K s n is close to K s . In addition, the anti-resonant frequency ω a r can be obtained by a free oscillation experiment, in which the motor is controlled to rotate to a certain position and then locked by the brake. Figure 8 shows the arm-side position in free oscillation experiment, from which we can know that the oscillation period T n is about 0.214 s. So, the inertia of motor side J n can be calculated as 0.528 Kg·m2.

3.2. Disturbance Observer

Robustness is essential to design a high-performance motion control system in practice. A DOb. which estimates external disturbances and system uncertainties, such as external load, friction, and inertia variation is widely used to achieve the robustness of motion control systems due to its simplicity and efficiency [34,35]. In this paper, a motor-side DOb and an arm-side DOb are both employed, as shown in Figure 9. Here, the disturbance torques on the motor side and arm side are defined in Equations (8) and (9), respectively. The outputs of the DObs, which can be described as the following two equations, are the estimation values of external disturbance:
D ^ m = g m s + g m ( D m + τ s n ) ;
D ^ l = g l s + g l D l ,
where g m and g l are the bandwidths of motor-side and arm-side DOb respectively. The bandwidth of DOb is directly related to the robustness of the motion control system. The higher the bandwidth of DOb is, the more robustness can be expected. It should be noted that the actual bandwidth of DOb is limited by practical constraints such as noise, sampling time, etc. [36].

3.3. Resonance Ratio Control

As a novel control method, the resonance ratio control first proposed by Yuki et al. [37] can suppress the vibration for the position control in two-mass resonant system. If the spring torque multiplying K Γ M n is feedback to the two-mass resonant system, the system block diagram can be transformed from Figure 7 into Figure 10. When K Γ = 1 / M n , the system shown in Figure 10 has equivalent performance of the uncontrolled plant shown in Figure 7.
Then, the anti-resonant frequency ω a r and the resonant frequency ω a of the new system can be described as follows:
ω a r = K s n J n ;
ω r = K s n J n ( 1 + K Γ J n ) = K ω a r ,
where K is called “Resonance Ratio” and defined as the following equation:
K = ( 1 + K Γ J n ) .
It is apparent from the definition in Equation (18) that the resonance ratio is a monotone increasing function of the spring torque feedback gain K Γ . This means that the resonance frequency and the resonance ratio of the system can be tuned to an arbitrary value. Intuitively, the role of K Γ is to depict to what extent the motor side is affected by the spring torque and represent the sensitivity of the motor side against the spring torque, i.e., the arm side. A PD control in the motor side is constructed as shown in Figure 11. This controller consists of a proportional gain K p , a differential gain K υ , and a spring torque feedback gain K Γ .
From Figure 11, the transfer functions from q d e s to θ and q are given as follows, respectively:
θ = K p ( s 2 + ω a r 2 ) s 4 + K υ s 3 + ( K p + ω r 2 ) s 2 + K υ ω a r 2 s + K p ω a r 2 q d e s ;
q = K p ω a r 2 s 4 + K υ s 3 + ( K p + ω r 2 ) s 2 + K υ ω a r 2 s + K p ω a r 2 q d e s .
In the position RRC, it is shown that the vibration can be suppressed by controlling the ratio between resonant and anti-resonant frequencies via spring torque feedback gain K Γ . Expressions of each control gain are given as follows [37]:
K = ω r ω a r = 5 ;
K p = ω a r 2 ,
K υ = 4 ω a r ;
K Γ = 4 / J n .
In order to improve the robustness of the position control system, a motor-side DOb and an arm-side DOb explained in Section 3.2 are employed here, as shown in Figure 12.
Robust position control of the motor side is achieved by feedback to the inner-loop DOb. The arm-side disturbance, estimated by the arm-side DOb, is compensated by feedback through the inverse system. The inverse system can be calculated by the transfer function from q d e s to q as shown in Equation (20) and from D l to q as follows:
q D l = ( s 2 + K υ s + K p + K s n K Γ ) / J n s 4 + K υ s 3 + ( K p + ω r 2 ) s 2 + K υ ω a r 2 s + K p ω a r 2 q d e s .
Then, the inverse system can be represented by the following equation:
q d e s D l = s 2 + K υ s + K p + K s n K Γ K p K s n .
For the sake of simplicity, the following approximate transfer function is generally used with an arm-side DOb:
q d e s D l = s 2 + K υ s + K p + K s n K Γ K p K s n K p + K s n K Γ K p K s n .
It should be noted that, only constant disturbances can be suppressed if the higher order derivatives are neglected in Equation (27).

4. Simulation and Experimental Results

Firstly, numerical simulation is given to confirm the validity of the proposed control approach. In the simulation, the sampling period of simulation system was 2 ms and the real model parameters were set to 20% larger than the nominal counterparts. The arm-side inertia varies a lot with the load in practice, so the real arm-side inertia was also set to double for comparison. The load and the arm length were set to 4 Kg and 0.3 m respectively. Coulomb friction and viscous friction coefficients of the motor side as well as the arm side were 0.05 Nm and 0.1 Nm·s/rad respectively. Controller gains, nominal model parameters of the elastic joint obtained in Section 3.1 and the bandwidths of DObs are listed in Table 2.
Figure 13a–d shows the simulation results of position control when a step reference input ( q d e s = 60 ° ) was applied. If the arm-side disturbance was neglected, the position RRC without motor-side DOb is shown in Figure 13a. Due to large K p and small motor-side disturbance, there was almost no discrepancy between q and q d e s . It is noteworthy that arm-side position responded rapidly, but there was a small overshoot. If external torque of 20 Nm was deliberately applied to the motor side, the steady-state error appeared. As it is shown in Figure 13b, the steady-state error caused by motor-side disturbance can be eliminated by virtue of the proposed motor-side DOb. It is clear from the figure that the motor-side DOb improved not only the steady-state performance but also the transient response of the system. Figure 13c shows the result of position RRC with DObs. An arm-side external torque of 10 Nm was applied at 3 s. After a short adjustment period, the arm-side position then kept following the step reference input quickly and smoothly. The robustness of the system was improved by using the arm-side DOb. Moreover, the proposed DOb-based RRC was robust against much larger disturbance on the nominal model, as shown in Figure 13c. The outputs of the arm-side DOb and motor-side DOb are shown in Figure 13d. Besides system uncertainties such as model variation, friction, and load, the estimated value of the arm-side disturbance contained external torque. Figure 13a–c clearly shows that the vibration can be suppressed by using the proposed RRC method. Simulation results show the validity of the DOb-based RRC.
Secondly, to validate the performance of the modular elastic joint and the DOb-based RRC, an experimental platform was established as shown in Figure 14. This platform mainly consisted of six parts, i.e., base frame, elastic joint with load, joint controller, Compact PCI-based cabinet, DC power and monitor. An input link of the elastic joint was fixed on the base frame and an output one held the load. All the actual parameters of the elastic joint were obtained by corresponding calibration experiments and listed in Table 1. From Table 1, it can be concluded that the designed elastic joint basically meets the design requirements.
The position control experiments of the elastic joint were carried out with the same control parameters as in simulation and the results are shown in Figure 15. Figure 15a shows the step responses when conventional PD controller was employed. There was a steady-state error mainly due to arm-side load when PD controller was used. The steady-state error can be decreased if the gravity component of load is compensated. But, the steady-state error cannot be eliminated completely because the gravity component of load cannot be accurately estimated. Besides, uncertainties of friction, joint stiffness and external torque will result in the steady-state error. As it is shown in the figure, another way to reduce the steady-state error was to increase K p . However, the stability deteriorated significantly. Figure 15c shows the comparison between PD control with gravity compensation and DOb-based RRC. When the position error was large, increasing rates of arm-side position in both control strategies were identical and basically equal to the maximum angular velocity of the elastic joint. Arm-side angular velocity fell rapidly in DOb-based RRC when the position error was small and arm-side position followed the reference position without vibration. Control torques of conventional PD controller and the proposed DOb-based RRC are shown in Figure 15d. The proposed DOb-based RRC avoided the oscillation of the control input. Figure 15e,f shows the trajectory tracking performances of the designed joint for sinusoidal references of different frequencies. It is clear from the Figure 15e that the proposed DOb-based RRC showed relatively high tracking performance. As shown in Figure 15f, the trajectory tracking performance deteriorates and the delay grew as the frequency increased, especially when the reference frequency was more than 2 Hz. It is easy to understand that the elasticity of the joint mainly answers for the delay and the DOb-based RRC for achieving vibration suppression will further pull down the system bandwidth. Experimental results show that the vibration can be suppressed by means of RRC and the excellent performance of the DOb-based RRC in point-to-point position control.

5. Conclusions

This paper presents the design of a compact modular rotary elastic joint for a humanoid manipulator. A customized torsional spring is designed based on FEM, which meets the performance requirements of admissible peak load, low stiffness, and compactness. The actual value of the spring stiffness, which was calibrated by means of experiment, was significantly higher than obtained by simulation. The difference between the actual and nominal properties of the material, the imperfections in the model used in the simulation and measurement errors probably led to the discrepancy. Fortunately, this discrepancy only makes the stiffness value of the elastic joint different from the expectation and has little effect on the performance of the elastic joint. The large central hole and fail-safe brake were respectively designed for cabling and power-off protection of the elastic joint. The designed elastic joint only used magnetic encoders for sensor information, which greatly reduces the cost. Employment of double-layer control structure in electrical system can reduce the burden on the controller and simplify system maintenance. Besides, the designed electrical system can be easily expanded to build a humanoid manipulator in future. The PPSeCo protocol employed in this paper for data communication between central controller and joint controller can ensure the transmission reliability, speed of data communication, and reduce noise in sensor signals. The specifications listed in Table 1 demonstrate the elastic joint suitable as a module for a humanoid manipulator.
This paper also proposed a high-performance robust point-to-point position control scheme for the designed elastic joint. Vibration of the elastic joint was suppressed by virtue of RRC, which can arbitrarily designate the ratio between the resonant and anti-resonant frequency according to the feedback of the nominal spring torsion. Moreover, robustness of the control system was achieved by using motor-side DOb and arm-side DOb.
The validity of the DOb-based RRC was demonstrated in the simulation results. Experimental results further showed the performances of the modular elastic joint and the efficacy of the controller. Design and control of a multi-DOF humanoid manipulator, which was equipped with the designed elastic joint, are the focus of our future work.

Author Contributions

S.C. and Y.L. designed the elastic joint and built the experimental platform. Y.R. designed the DOb-based RCC and carried out experiments of the elastic joint. Y.S. built the electrical system. H.L. wrote the paper, guided the experiments and put forward suggestions as the head of the research group that conducted this study.

Acknowledgments

This work is supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China under grant No. 51521003, and the Major Research Plan of the National Natural Science Foundation of China under grant No. 61690210, and the National Key R&D Program of China under grant No. 2017YFB1300400.

Conflicts of Interest

The authors declare no conflicts of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript.
PPSeCoPoint-to-Point High Speed Serial Communication
VHDLVery-High-Speed Integrated Circuit Hardware Description Language
FPGAField Programmable Gate Array
PCIPeripheral Component Interconnect
DSPDigital Signal Processor
MFLOPSMillion Floating-point Operations per Second

References

  1. Ren, Y.; Chen, Z.; Liu, Y.; Gu, Y.; Minghe, J.; Liu, H. Adaptive hybrid position/force control of dual-arm cooperative manipulators with uncertain dynamics and closed-chain kinematics. J. Frankl. Inst. 2017, 354, 7767–7793. [Google Scholar] [CrossRef] [Green Version]
  2. Hogan, N. Impedance control—An approach to manipulation. I—Theory. II—Implementation. III—Applications. J. Dyn. Syst. Meas. Control 1985, 107, 1–24. [Google Scholar] [CrossRef]
  3. Ott, C.; Mukherjee, R.; Nakamura, Y. Unified impedance and admittance control. In Proceedings of the IEEE International Conference on Robotics and Automation, Anchorage, AK, USA, 3–7 May 2010; pp. 554–561. [Google Scholar]
  4. Pratt, G.A.; Williamson, M.M. Series elastic actuators. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Pittsburgh, PA, USA, 5–9 August 1995; pp. 399–406. [Google Scholar]
  5. Paine, N.; Oh, S.; Sentis, L. Design and Control Considerations for High-Performance Series Elastic Actuators. IEEE/ASME Transact. Mechatron. 2014, 19, 1080–1091. [Google Scholar] [CrossRef] [Green Version]
  6. Tsagarakis, N.G.; Laffranchi, M.; Vanderborght, B.; Caldwell, D.G. A compact soft actuator unit for small scale human friendly robots. In Proceedings of the IEEE International Conference on Robotics and Automation, Kobe, Japan, 12–17 May 2009; pp. 4356–4362. [Google Scholar]
  7. Cummings, J.P.; Ruiken, D.; Wilkinson, E.L.; Lanighan, M.W.; Grupen, R.A.; Sup, F.C., IV. A Compact, Modular Series Elastic Actuator. J. Mech. Robot. 2016, 8. [Google Scholar] [CrossRef]
  8. Deng, K.; Zhao, M.; Xu, W. Level-ground walking for a bipedal robot with a torso via hip series elastic actuators and its gait bifurcation control. Robot. Autonom. Syst. 2016, 79, 58–71. [Google Scholar] [CrossRef]
  9. Sergi, F.; Lee, M.M.; Malley, M.K.O. Design of a series elastic actuator for a compliant parallel wrist rehabilitation robot. In Proceedings of the 13th IEEE International Conference on Rehabilitation Robotics, Seattle, WA, USA, 24–26 June 2013; pp. 1–6. [Google Scholar]
  10. Lagoda, C.; Schouten, A.C.; Stienen, A.H.A.; Hekman, E.E.G.; van der Kooij, H. Design of an electric series elastic actuated joint for robotic gait rehabilitation training. In Proceedings of the International Conference on Biomedical Robotics and Biomechatronics, Tokyo, Japan, 26–29 September 2010; pp. 21–26. [Google Scholar]
  11. Kong, K.; Bae, J.; Tomizuka, M. A compact rotary series elastic actuator for knee joint assistive system. In Proceedings of the IEEE International Conference on Robotics and Automation, Anchorage, AK, USA, 3–7 May 2010; pp. 2940–2945. [Google Scholar]
  12. Tsagarakis, N.G.; Sardellitti, I.; Caldwell, D.G. A new variable stiffness actuator (CompAct-VSA): Design and modelling. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, San Francisco, CA, USA, 25–30 September 2011; pp. 378–383. [Google Scholar]
  13. Cui, S.; Liu, Y.; Sun, Y.; Liu, H. Design of a Series Variable Stiffness Joint Based on Antagonistic Principle. In Proceedings of the International Conference on Intelligent Robotics and Applications, Wuhan, China, 16–18 August 2017; pp. 238–249. [Google Scholar]
  14. Robinson, D.W. Design and Analysis of Series Elasticity in Closed-Loop Actuator Force Control. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, June 2000. [Google Scholar]
  15. Vallery, H.; Veneman, J.; van Asseldonk, E.; Ekkelenkamp, R.; Buss, M.; Kooij, H. Compliant Actuation of Rehabilitation Robots—Benefits and Limitations of Series Elastic Actuators. IEEE Robot. Autom. Mag. 2008, 15, 60–69. [Google Scholar] [CrossRef]
  16. Tomei, P. A simple PD controller for robots with elastic joints. IEEE Trans. Autom. Control 1991, 36, 1208–1213. [Google Scholar] [CrossRef]
  17. Kim, M.J.; Chung, W.K. Disturbance-Observer-Based PD Control of Flexible Joint Robots for Asymptotic Convergence. IEEE Trans. Robot. 2015, 31, 1508–1516. [Google Scholar] [CrossRef]
  18. Jin, M.; Lee, J.; Tsagarakis, N.G. Model-Free Robust Adaptive Control of Humanoid Robots With Flexible Joints. IEEE Trans. Ind. Electron. 2017, 64, 1706–1715. [Google Scholar] [CrossRef]
  19. Palli, G.; Melchiorri, C.; Luca, A.D. On the feedback linearization of robots with variable joint stiffness. In Proceedings of the IEEE International Conference on Robotics and Automation, Pasadena, CA, USA, 19–23 May 2008; pp. 1753–1759. [Google Scholar]
  20. Petit, F.; Albu-Schäffer, A. State feedback damping control for a multi DOF variable stiffness robot arm. In Proceedings of the 2011 IEEE International Conference on Robotics and Automation, Shanghai, China, 9–13 May 2011; pp. 5561–5567. [Google Scholar]
  21. Miyazaki, F.; Kawamura, S.; Matsumori, M.; Arimoto, S. Learning control scheme for a class of robot systems with elasticity. In Proceedings of the 25th IEEE Conference on Decision and Control, Athens, Greece, 10–12 December 1986; pp. 74–79. [Google Scholar]
  22. Sardellitti, I.; Medrano-Cerda, G.A.; Tsagarakis, N.; Jafari, A.; Caldwell, D.G. Gain Scheduling Control for a Class of Variable Stiffness Actuators Based on Lever Mechanisms. IEEE Trans. Robot. 2013, 29, 791–798. [Google Scholar] [CrossRef]
  23. Petit, F.; Ott, C.; Albu-Schäffer, A. A model-free approach to vibration suppression for intrinsically elastic robots. In Proceedings of the IEEE International Conference on Robotics and Automation, Hong Kong, China, 31 May–7 June 2014; pp. 2176–2182. [Google Scholar]
  24. Kapucu, S.; Baysec, S.; Alici, G. Residual Vibration Suppression of a Flexible Joint System Using a Systematic Command Shaping Technique. Arab. J. Sci. Eng. 2006, 31, 139–152. [Google Scholar]
  25. Bang, J.S.; Shim, H.; Park, S.K.; Seo, J.H. Robust Tracking and Vibration Suppression for a Two-Inertia System by Combining Backstepping Approach With Disturbance Observer. IEEE Trans. Ind. Electron. 2010, 57, 3197–3206. [Google Scholar] [CrossRef]
  26. Conker, Ç.; Yavuz, H.; Bilgiç, H. A review of command shaping techniques for elimination of residual vibrations in flexible-joint manipulators. J. Vibroeng. 2016, 18, 2947–2958. [Google Scholar] [CrossRef]
  27. Dang, J.; Ni, F.; Gu, Y.; Jin, M.; Liu, H. A highly integrated and flexible joint test system based on DSP/FPGA-FPGA. In Proceedings of the IEEE International Conference on Robotics and Biomimetics, Guilin, China, 19–23 December 2009; pp. 1877–1882. [Google Scholar]
  28. Liu, H.; Meusel, P.; Seitz, N.; Willberg, B.; Hirzinger, G.; Jin, M.H.; Liu, Y.W.; Wei, R.; Xie, Z.W. The modular multisensory DLR-HIT-Hand. Mech. Mach. Theory 2007, 42, 612–625. [Google Scholar] [CrossRef]
  29. Yoon, S.-S.; Kang, S.; Kim, S.-J.; Kim, Y.-H.; Kim, M.; Lee, C.-W. Safe arm with MR-based passive compliant joints and visco-elastic covering for service robot applications. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Las Vegas, NV, USA, 27–31 October 2003; pp. 2191–2196. [Google Scholar]
  30. Carpino, G.; Accoto, D.; Sergi, F.; Tagliamonte, N.; Guglielmelli, E. A Novel Compact Torsional Spring for Series Elastic Actuators for Assistive Wearable Robots. J. Mech. Des. 2012, 134, 1–10. [Google Scholar] [CrossRef]
  31. Dos Santos, W.M.; Caurin, G.A.P.; Siqueira, A.A.G. Design and control of an active knee orthosis driven by a rotary Series Elastic Actuator. Control Eng. Pract. 2017, 58, 307–318. [Google Scholar] [CrossRef]
  32. Tee, K.P.; Burdet, E.; Chew, C.M.; Milner, T.E. A model of force and impedance in human arm movements. Biol. Cybern. 2004, 90, 368–375. [Google Scholar] [CrossRef] [PubMed]
  33. Wolf, S.; Eiberger, O.; Hirzinger, G. The DLR FSJ: Energy based design of a variable stiffness joint. In Proceedings of the IEEE International Conference on Robotics and Automation, Shanghai, China, 9–13 May 2011; pp. 5082–5089. [Google Scholar]
  34. Sariyildiz, E.; Yu, H.; Yu, K.; Ohnishi, K. A nonlinear stability analysis for the robust position control problem of robot manipulators via disturbance observer. In Proceedings of the IEEE International Conference on Mechatronics (ICM), Nagoya, Japan, 6–8 March 2015; pp. 28–33. [Google Scholar]
  35. Sariyildiz, E.; Chen, G.; Yu, H. An Acceleration-Based Robust Motion Controller Design for a Novel Series Elastic Actuator. IEEE Trans. Ind. Electron. 2016, 63, 1900–1910. [Google Scholar] [CrossRef]
  36. Sariyildiz, E.; Ohnishi, K. Stability and Robustness of Disturbance-Observer-Based Motion Control Systems. IEEE Trans. Ind. Electron. 2015, 62, 414–422. [Google Scholar] [CrossRef]
  37. Yuki, K.; Murakami, T.; Ohnishi, K. Vibration control of 2 mass resonant system by resonance ratio control. In Proceedings of the International Conference on Industrial Electronics, Maui, HI, USA, 15–19 November 1993; pp. 2009–2014. [Google Scholar]
Figure 1. The design process of the elastic joint. (a) Traditional rigid joint with torque sensor; (b) Cross section view of the prototype elastic joint; (c) Physical prototype of the elastic joint.
Figure 1. The design process of the elastic joint. (a) Traditional rigid joint with torque sensor; (b) Cross section view of the prototype elastic joint; (c) Physical prototype of the elastic joint.
Sensors 18 01869 g001
Figure 2. Electrical system diagram of the elastic joint.
Figure 2. Electrical system diagram of the elastic joint.
Sensors 18 01869 g002
Figure 3. The designed torsional spring. (a) Schematic perspective view of the customized torsional spring; (b) Physical model of the customized torsional spring.
Figure 3. The designed torsional spring. (a) Schematic perspective view of the customized torsional spring; (b) Physical model of the customized torsional spring.
Sensors 18 01869 g003
Figure 4. Results of the methodology adopted. The arrows indicate the direction of increasing spring thickness. (a) Geometry parameter L vs. Von Mises stress; (b) Geometry parameter L vs. Angular deformation of the spring.
Figure 4. Results of the methodology adopted. The arrows indicate the direction of increasing spring thickness. (a) Geometry parameter L vs. Von Mises stress; (b) Geometry parameter L vs. Angular deformation of the spring.
Sensors 18 01869 g004
Figure 5. Analysis results of spring by FEM. (a) Static simulation for stress distribution; (b) Static simulation for angular deformation.
Figure 5. Analysis results of spring by FEM. (a) Static simulation for stress distribution; (b) Static simulation for angular deformation.
Sensors 18 01869 g005
Figure 6. Characteristics of the customized torsion spring.
Figure 6. Characteristics of the customized torsion spring.
Sensors 18 01869 g006
Figure 7. Model of two-mass resonant system.
Figure 7. Model of two-mass resonant system.
Sensors 18 01869 g007
Figure 8. Results of the arm-side free oscillation experiment. (a) Arm-side position; (b) Magnified view of (a).
Figure 8. Results of the arm-side free oscillation experiment. (a) Arm-side position; (b) Magnified view of (a).
Sensors 18 01869 g008
Figure 9. Block diagram of motor-side DOb and arm-side DOb.
Figure 9. Block diagram of motor-side DOb and arm-side DOb.
Sensors 18 01869 g009
Figure 10. Spring torque feedback multiplying KΓMn.
Figure 10. Spring torque feedback multiplying KΓMn.
Sensors 18 01869 g010
Figure 11. PD control of the two-mass resonant system.
Figure 11. PD control of the two-mass resonant system.
Sensors 18 01869 g011
Figure 12. The position resonance ratio control based on arm-side and motor-side DOb.
Figure 12. The position resonance ratio control based on arm-side and motor-side DOb.
Sensors 18 01869 g012
Figure 13. Simulation results of position RRC. (a) Position RRC without DObs; (b) Position RRC with motor-side Dob; (c) Position RRC with DObs; (d) Outputs of DObs.
Figure 13. Simulation results of position RRC. (a) Position RRC without DObs; (b) Position RRC with motor-side Dob; (c) Position RRC with DObs; (d) Outputs of DObs.
Sensors 18 01869 g013
Figure 14. Experimental platform of the designed elastic joint.
Figure 14. Experimental platform of the designed elastic joint.
Sensors 18 01869 g014
Figure 15. Experimental results of position control. (a) PD position control; (b) Magnified view of (a); (c) Comparison between PD control and RRC control; (d) Joint control torque; (e) Sinusoidal trajectory tracking under 0.2 Hz; (f) Sinusoidal trajectory tracking under different frequencies.
Figure 15. Experimental results of position control. (a) PD position control; (b) Magnified view of (a); (c) Comparison between PD control and RRC control; (d) Joint control torque; (e) Sinusoidal trajectory tracking under 0.2 Hz; (f) Sinusoidal trajectory tracking under different frequencies.
Sensors 18 01869 g015
Table 1. Design indexes and actual values of the elastic joint.
Table 1. Design indexes and actual values of the elastic joint.
ParametersDesign IndexesActual Values
Maximum elastic output torque≥60 Nm67 Nm
Maximum angular velocity≥0.6 rad/s0.72 rad/s
Maximum diameter≤120 mm111 mm
Maximum length≤150 mm136 mm
Mass≤2.5 Kg2.4 Kg
Joint stiffness≈350 Nm/rad464.2 Nm/rad
Maximum angular range360°360°
Table 2. Control parameters of the simulation and experiment.
Table 2. Control parameters of the simulation and experiment.
ParametersDescriptionsValues
J n Nominal arm-side inertia0.528 Kg·m2
M n Nominal motor-side inertia1.567 Kg·m2
K s n Nominal joint stiffness464.2 Nm/rad
K p Proportional gain862
K υ Differential gain117.44
K Γ Feedback gain7.58
g m Motor-side bandwidth500 rad/s
g l Arm-side bandwidth500 rad/s

Share and Cite

MDPI and ACS Style

Liu, H.; Cui, S.; Liu, Y.; Ren, Y.; Sun, Y. Design and Vibration Suppression Control of a Modular Elastic Joint. Sensors 2018, 18, 1869. https://doi.org/10.3390/s18061869

AMA Style

Liu H, Cui S, Liu Y, Ren Y, Sun Y. Design and Vibration Suppression Control of a Modular Elastic Joint. Sensors. 2018; 18(6):1869. https://doi.org/10.3390/s18061869

Chicago/Turabian Style

Liu, Hong, Shipeng Cui, Yiwei Liu, Yi Ren, and Yongjun Sun. 2018. "Design and Vibration Suppression Control of a Modular Elastic Joint" Sensors 18, no. 6: 1869. https://doi.org/10.3390/s18061869

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop