Modeling of Linear Displacement and Design of Adaptive Robust Controller

Twisted-string actuation method has attracted much more attention due to its compact size and considerable output force. Many studies have been proposed to estimate its linear displacement of actuation based on the angular measurement of the rotary device. However, the estimated results are not always accurate, which requires an additional process to make the estimation accuracy. In this paper, an ARC control scheme with force feedback has been developed for a high-accuracy twisted-string actuation system.


Introduction
As robotic applications have started to gain much more attention recently, it is desired to have actuators that can generate large torque while keeping its size tiny. However, with the addition of the required power source, most available solutions of actuators are still bulky. ese actuators might not be too heavy for the applications of wearable robotic devices. Among all different actuators, twisted-string actuation has the advantages of lightweight by moving bulky components to other locations while providing enough force or torque to drive moving parts. It has attracted many research groups to adopt this twisted-string actuator in various exoskeleton robotic applications [1][2][3][4].
A twisted-string actuation system comprises of two or more parallel strands that convert rotational movement to linear displacement. e configuration is shown in Figure 1. Once the string starts to be twisted, the length is shortened. With carefully selected material and an accurate model, such an actuation method can serve as an actuator that delivers linear movement. It also has the advantage of keeping the rotary driving component, which is typically bulky, at a different location to keep the weight at the location of output torque, or force, being light. Such actuators have attracted much attention since Suzuki et al. adopted it to control the stiffness of the links in the application of robotic hands [5].
Würtz et al. have presented a simple kinetostatic model of a twisted-string-based transmission system [6].
However, the simple model does not describe specific nonlinear phenomena in the physical equations. us, more complex models have been proposed to better approximate the contraction and extension for the twisted-string actuation [7,8]. ese models discuss more on the string with dual-strand, which might limit the applicable range of such an actuation device. Shisheie et al. have investigated the performance of various combinations of strands [9]. eir result has illustrated that the trajectories of hysteresis are different when the number of integrated strands is changed. A string with more strands can yield larger friction. However, the model has still not been investigated yet. A model that can accurately estimate the actuated length of both contraction and expansion is particularly important since the changed length typically cannot be measured directly. Available sensors can only measure other physical properties and estimate the real length change indirectly. ough the modified model proposed by Gaponov et al. [8], it does not perfectly estimate the internal reaction among strands. It is also challenging to determine the output force of such kind of actuator. In this study, a nonlinear approach is proposed to enhance the accuracy of the estimation of length estimation and estimate the output force of twisted-string actuation.
Instead of using an expensive optical linear encoder, this study adopted an adaptive robust control algorithm with the force feedback to compensate for nonlinear factors of the twisted-string actuation.

Experimental Setup
To validate different mathematical models and the corresponding reliabilities, an experimental platform was built to characterize the relationship between the revolutions of the driving motor and the amount of contraction or expansion of the twisted-string using different string types and configurations.
To ensure the contraction and expansion to be linear movements, two 400 mm Richelieu slides were used to guide the movements. Figure 2 illustrates the schematic configuration of the experimental setup. In this setup, a DC motor was placed on the upper side of the testbed with a string holder attached to the shaft. e string holder has 7 holes with one hole located at the center and 6 holes evenly on a circular path [9]. Depending on the types of tested string, as shown in Figure 3, different types of strings with the required number of strands were used in individual tests. e same holder was fixed on a slider, which transforms the rotational motion into linear movement. e original length of a single strand is designed to be 500 mm or 1,000 mm. e DC motor, a high-torque Pololu motor with a 19 : 1 gear ratio and a 64 CPR encoder, used to twist the string is driven by a PWM mode L298 H-Bridge. In addition to the rotary encoder, several sensors were equipped as well. A WXY linear encoder is used to measure the linear displacements of the twisted-string; the resolution of linear encoder is 0.04 mm. A high precision sliding potentiometer was also used to detect the linear displacement. An AC712 current sensor is used to measure the current of the armature of the DC motor. A SHIMPO FGV-XY20 is used to measure the tensile force generated by the twisted-string, the resolution of the force gauge 0.1 N. To acquire the movement and force data, a desktop PC with a Pentium IV 2.4 GHz and 1 GB memory is used for data acquisition. Two data acquisition I/O cards (PCI-6221 and PCI-6601) manufactured by the National Instrument is used to generate command signals to motor drivers and to acquire signals from position sensors, current sensor, and force sensor. e data acquisition program was developed using MATLAB Simulink with xPC Toolbox. e sampling rate is 100 Hz.
To guarantee smooth movement during each test, a mass was added to the lower end of the twisted-string. e friction of the guided movement in both moving directions was investigated ahead. With the proposed experimental setup, the friction of upward motion is around 0.52 N, and it is around 0.14 N for the downward movement.

Models of Twisted-String Actuation
e basic concept of the twisting actuation method is simple. With one end fixed on an actuator shaft (electric motor) and the other end fixed on the load, two or more parallel strands can be twisted by the rotational movement of the attached actuator. As the twist starts, the string connected to the motor shaft starts to reduce the distance between both ends and generates linear motion.
A mathematical model for such a transmission mechanism has been proposed by Würtz et al. [6]. is model assumes the stiffness of adopted strands being infinite, which can only be correct if the exerted force on the strands is small relative to the elastic modulus of the string. ere are eight configurations of strand combinations used by twistedstring actuation, distinguished by different strands of actuation and with/without a neutral strand, as shown in Figure 3. One configuration includes a neutral strand, which increases the effective radius about as the surrounding strands wrap the neutral one. e other configuration does not include a neutral strand. e distance r between the effective strands and center of rotation becomes r � r s + r c , where r s is the radius of the central neutral strand, and r c is the radius of the effective strands.   e kinematic relationship between the motor rotation θ and the length of the transmission system p is derived according to Pythagoras' theorem, which is shown in Figure 4. e actuation length p is a function of the number of rotations of the actuator: (1) where L is the total length of the string. Starting from completely untwisted-strings, the change in the actuation length Δp will be To maintain the static equilibrium at different positions, sufficient torque τ m must be generated by the actuator. In other words, to stay at equilibrium, the static compound force (nF i ) generated by F z as the tangential force (load), and F τ as the axial force corresponding to τ m must be aligned with the physical direction of the unwrapped helix. To provide dynamic actuation, the provided torque should be more or less than the required static torque. e amount of actuator force depends directly on the amount of load and helix slope of the strand, which is where n and α are the number of strands and helix slope, respectively. By combining equations of (3), the torque τ m becomes In cases of high loads or accelerated motion, the assumption of infinite stiffness of the string is no longer valid and the elongation of the strands due to the axial force must be taken into account. In this case, the total length of a string L is a function of the tension F i , the strand stiffness k, and its unloaded length L 0 . us, (1) can then be modified to With (5), the equation to determine the length of pretwist actuation can be explained as where θ 0 is the number of revolutions of the actuator associated with the pretwist. Assuming that this is the initial length of the cable, p(θ) as a function of the revolution of the motor can be derived as follows: where P(θ, F i ) is the amount of displacement from the pretwist position. According to the geometric relationship, there exists a maximum angle α where all strands are wrapped tightly on each other. After this point, strands begin entangling around themselves. With the maximum angle α, the pitch parameter q of the helix is minimized: On the other hand, the pitch parameter of a helix as a function of motor rotation angle θ and string length p will be 2πp � qθ. (9) Assuming no load (F z � 0) was applied to the transmission and combining (5) with the geometry of the strands, (8) and (9), the maximum achievable θ angle (θ max ) and minimum contraction (p min ) are, respectively,

Verification of Theoretical Estimation
To verify the correctness of the theoretical estimation and select the optimal configuration and material of the strands for the twisted-string actuation which requires low hysteresis, quick response, and better durability, three types of material that were tested: a monofilament fishing line, a utility cable, and a tactical cable. In these tests, the behavior of the actuation within the low-hysteresis area was investigated. Since the actuations were twisted to the maximum limit and in some cases beyond, the low-hysteresis area was determined for each configuration by testing firstly. en, the actuation was twisted to that limit, and its behavior was recorded. e response of the actuation beyond this limit demonstrates hysteresis, and the higher the number of revolutions becomes, the bigger the hysteresis will be. e monofilament fishing line was the first test, which was tested with two different lengths: 500 mm and 1,000 mm. To acquire the result without losing generality, the results were obtained by performing each test for 50 cycles and selecting one cycle as sampled data after each consecutive 10 cycles. In this test, the experiment was including a 500 g loading. e utility cable and tactical cables were tested for 8 different configurations and under 4 different loadings: 500 g, 1,500 g, 2,500 g, and 3500 g. ese tests were done 5 times for each configuration under each load in order to ensure the consistency of the response and to demonstrate the hysteresis in the response if it exists. e actuation length of the tests is set to be 500 mm and 1,000 mm anymore.

Monofilament Fishing Line.
For the monofilament fishing lines (TEBCO Outcast), the rated maximum working load is 88N and the diameter of each strand of the fishing line is 0.46 mm. e shape of the filament was confirmed to be perfectly round by measuring the diameter along the length before testing. However, after several testing, not only the diameter of the string but also the overall shape of the string may be affected based on how much it was twisted or whether it was deformed beyond its flexibility range. e configuration of the string for this experiment is a twostrand string. Before it is failed due to the repetitive movements, the variation of the diameters is 0.5 ± 0.05 mm. Figure 5 shows the relationship between the contracted length and motor revolution based on experimental results and theoretical calculations. It is clear that the basic model cannot precisely predict the length due to the elastic behavior caused by the additional force created by the attached weight or other factors. From the experiment, the adopted line starts to suffer from permanent distortion after 410 revolutions. e inaccuracy can be patched by the modified model with a calibrated stiffness coefficient.
To ensure the consistency of the twisted-string mechanism, the strands must be tested under repetitive motion/ force. Figure 6 demonstrates the lengths of the string of repetitive motions of both contraction and expansion for different string lengths (500 mm and 1,000 mm). From the test, the variations of contraction do not match with the trajectories of expansion, which implies that the contraction/ expansion process is nonlinear and some hysteresis exists in this type of actuation method.  Figure 7, hysteresis increase as the actuator revolution increases, and after reaching a peak at half of the motor revolution, it decreases. e same trend was observed for both contraction and expansion. e maximum hysteresis for 1,000 mm long string is 1 cm and for 500 mm long string 0.4 cm. If the strands were twisted more than the maximum allowed revolutions, the line begins kinking, which can yield unexpected nonlinear behavior. Kinking also distorts strands permanently and is irreversible, as shown in Figure 8. In addition, though modified model can be used to predict the contraction/expansion length, the coefficient is not a constant and can be changed when the loading is different. e friction between the two strands was not taken into consideration, either. As the revolution of the shaft increases, the friction increases as well. is part has not been included in the modified model.

Utility Line.
e major drawbacks of the fishing line are its high elasticity and slow contraction. us, a more rigid string with higher capacity of applicable loading is desired. In this section, utility lines with 0.9 mm diameter and 133 N maximum working load was selected to work as the strands of the string. e responses of the 2-strand actuation with and without the neutral strand are shown in Figure 9(a). Four different loading were applied to the string. ey are 500 g, 1,500 g, 2,500 g, and 3,500 g masses. For each load, four repetitive experiments were performed and recorded. e twistedstring actuation has different possible motion trajectory, depending on the contraction or expansion of the strands. Plots of the motion trajectory of the twisted-string actuation often form a loop or hysteresis curve, where there are different trajectories depending on the contraction or expansion of the strands, as shown in Figures 7(a) and 7(b). In this study, the area in which the average difference between the contraction and the expansion trajectories was less than 0.5 mm was defined as a low-hysteresis area. e maximal allowable motor revolutions to keep the string within the low-hysteresis area for the actuation with or without neutral strand were 47 or 60 revolutions, subsequently. From the experimental result, adding a neutral strand can reduce the maximum average number of revolutions without going to the high-hysteresis area from almost 60 revolutions to 47 revolutions, which demonstrates almost 20% improvement in the actuation response. is reduction is because of the 100% increase in the effective radius (r) shown in Figure 3. Moreover, the maximum contraction, within the low-hysteresis area, of the average of four tests decreased from 150 mm for the actuation without a neutral strand to 135 mm with a neutral strand, which demonstrates a 10% reduction in full contraction. From Figure 9(a), it is clear that changed length was affected greatly due to the flexibility of the string, which implies that the modified model can be a better solution as long as the stiffness coefficient can be derived. With a neutral strand existing in the string can improve the precision of actuation. e actuation without a neutral strand shows a 30 mm variation under different loadings for the maximum contraction, while the same configuration with a neutral strand is a 22 mm variation, which is a 25% reduction.
For 3-strand string actuation with and without a neutral strand, the experimental results are shown in Figure 9(b). e maximum motor revolutions within low-hysteresis area with or without neutral strand are 38 or 55 revolutions, subsequently. Similar to the previous actuation configuration, adding a neutral strand can reduce the maximum average number of revolutions (without going to the highhysteresis area) from almost 55 revolutions to 38 revolutions. is reduction is created by the 73% increase of effective radius (r). e maximum contraction of the average of the four loadings decreases from ∼150 mm without a neutral strand to ∼110 mm with a neutral strand. e actuation without a neutral strand shows a 15 mm variation under different loadings for the maximum contraction, while the same configuration with a neutral strand is a 13 mm variation, which is a 13% reduction. e responses of the actuation for 4-, 5-, and 6-strand with a neutral strand are shown in Figure 9(c). e maximum number of motor revolutions to keep actuation within low-hysteresis area is 35, 31, and 31.5 revolutions, respectively. As expected, adding effective strands can reduce maximum revolutions without yielding high-hysteresis. e     e maximum number of motor revolutions to keep the string within lowhysteresis area with and without a neutral strand is 31.5 and 31 revolutions, respectively. e maximum contraction for these two strings is 140 mm and 115 mm for the actuation with and without the neutral strand. Unlike other configurations, adding a neutral strand does not significantly change the behavior of actuation. is is mainly due to the multiple effective strand in the configuration, which greatly reduce the effectiveness of the neutral strand. However, the deviation caused by the flexibility of adopted string is not improved with the added neutral strand. From the experimental results, the deviations are both 13 mm under different loadings at maximum contractions. e responses of actuations for all configurations under a 3,500 g loading are shown in Figure 10. It demonstrates the change of behavior with additional strands. Table 1 summarizes the characteristics of strings with different combination of effectiveness and neutral strands. From the results, it is clear that more effective strands can make the motion less sensitive to the applied loading while reducing the maximum contraction, which can also respond faster to the changing of motor speed.

Tactical Cable.
Although utility line can have fast contracting response, it still does not provide accurate precision, rigidity, and resistance against load. us, a more rigid material with higher loading limit is desired. In this section, a tactical cable with 2.032 mm thickness and 200 N (as stated in catalog) maximum working load was adopted as the strands.
is strand is made of an abrasion-resistant material (Technora), which has low stretch and high strength. Moreover, the PTFE (Teflon) coating increases the abrasion resistance and adds UV protection. e responses of 2-strand actuation with and without a neutral strand are shown in Figure 11(a). e maximum number of motor revolutions within the low-hysteresis area for the actuation with or without neutral strand is 50 or 63 revolutions, subsequently. Adding a neutral strand can yield a 20% improvement in the actuation response for full contraction. e maximum contraction of the average of four loading tests is decreased from 165 mm without a neutral strand to 155 mm with a neutral strand (6% reduction). e actuation without a neutral strand demonstrated 26 mm of variation under different loadings for maximum contraction. It is 16 mm with a neutral strand, which is a 38% reduction. e responses of the 3-strand actuation with and without the neutral strand are shown as Figure 11(b). e maximum motor revolutions within the low-hysteresis area with or without the neutral strand were 42 or 51 revolutions, respectively. Adding a neutral strand can improve the actuation response of full contraction by 19%. e maximum contraction of the average of four tests is decreased from 165 mm without a neutral strand to ∼130 mm with a neutral strand (16% reduction). e actuation without a neutral strand shows a 11 mm deviation under different loadings for its maximum contraction, while the same configuration with a neutral strand is only 9 mm, which is an 18% reduction. e responses of the actuation for 4 , 5 , and 6 effective strands and a neutral strand are shown in Figure 11(c). e maximum motor revolutions within the low-hysteresis area for the 4-, 5-, and 6-strand configurations are 36, 33, and 28 revolutions, respectively. As expected, adding effective strands can reduce the maximum average number of revolutions without going to high-hysteresis area. e maximum contraction of the average of four tests is 125, 125, and 110 mm for the 4-, 5-, and 6-strand strings, respectively. e maximum variation under different loadings at maximum contraction was 15 mm for 4-strand string actuation, 13 mm for 5-strand string actuation, and 9 mm for 6-strand string actuation. As to 6 effective strands with and without a neutral strand, the maximum motor revolutions within lowhysteresis area are 28 and 32 revolutions, respectively. e responses of actuations for all configurations under a 3,500 g loading are shown in Figure 12. e maximum contraction of the average of four loading tests is both 110 mm for both actuations with or without the neutral strand. Table 2 summarizes the characteristics of different configurations using tactile cable as the strands.
As shown in Table 2, adding a neutral strand decreases the maximum number of revolutions of the DC motor, which also increases the contraction speed. To demonstrate the trend of change in the behavior of the actuation system by adding more strings, the responses under a 3,500 g mass loading are shown in Figure 12

Durability of Strands.
One additional critical issue of twisted-string actuation is the durability of adopted strands and wearing resistance. Figure 13 illustrates the strands of different materials after 200 working cycles. In this study, after 200 working cycles, the fishline conditions are not shown in Figure 13 because it is hard to see the difference between twisted and untwisted fishline in the picture due to the tiny size of the fishing.
As shown in Figure 13, though strands have shown some defeats after 200 working cycles, the distortions do not    permanently damage the strands to affect the behavior of actuation. However, some permanent distortion can be observed. From Figure 13, it is clear that the tactile cable, which is abrasive resistive, shows much less distortion after intensive working cycles. us, it is suggested to check the status of the strands after certain number of working cycles. It is also important to know the lifecycle of adopted strands.
In addition to available lifecycle, some inappropriate operations can damage the strands as well. One most common operation is to rotate the twisted-string beyond the allowable range of adopted strands. While being operated beyond the allowable range, the string becomes kinking because it is over twisted. Once the string becomes kinking, the string can suffer from irreversible damage. However, the nonlinear area covers the response of the actuation after going to the unpredictable nonlinear range. e distribution of low-hysteresis, high-hysteresis, and nonlinear areas is shown in Figure 14.

ARC Control Synthesis of the Twisted String Actuation
Twisted-string actuation has the advantages of compact size, flexible mechanism, and lightweight at locations where actuations are required. However, though twisted-string actuation has many advantages, its unique properties of nonlinearity including hysteresis, friction, operational dead-zone, and poor repeatability can deteriorate the desired control performance. Improper control methods without consideration of nonlinearities can lead to limit cycles or instability. To deal with both nonlinear factors, adaptive robust control (ARC) was adopted for such a motion system. Although this kinetostatic model described in the last section can reflect the relationship between axial forces, torque   S is the number of effective strands. N is the number of neutral strand.   Figure 14: e response of a 2-stand string actuation in lowhysteresis, high-hysteresis, and nonlinear areas. generated by motor, and string length, there are still obvious differences between mathematical modeling and actual performance of a twisted string. e main reason for the differences can be caused by not considering of friction among strands, the gap among different strands, elasticity of single strands, and deformation of individual strands. erefore, the twisted-string actuation system falls into the class of the real-world processes that are not easy to obtain the mathematical representation. In real-time control system, there is a large variety of plants, similar to twisted-string actuator, whose system model is highly nonlinear and cannot be physically derived. Figure 14 demonstrates the relationship between motor revolutions and force generated by the twisted-string actuator. Since the force transmission of the twisted-string actuator plays the key role in controller design, the constant current source is used to drive the DC motor, which can ensure the torque generated by DC motor stays unchanged during the whole process of twisting string. e data, as shown in Figure 15, are obtained when the constant current source provides 0.7 Amp current. From Figure 15(a), it is obvious that the force transmission keeps increasing along with motor revolutions when the torque generated by the DC motor is constant, which means that the force transmission efficiency can be increased by pretwisting string to specify operation range. e two stages of transition are shown in Figures 15(a) and 15(b), and the two phases of the force transmission are initial phase, in which the twisted-string actuator does not generate the linear displacement, and normal phase. It is easy to make a conclusion that the force transmission will increase to infinity based on simply observing Figure 15(a); however, with the increasing of motor revolutions, the parts of torque generated by the DC motor used to overcome the friction among strands will also be increased.

Utility line a er 200 cycles
Based upon above experimental results and above analysis, it is clear that if we pretwist string before using the twisted-string actuator to drive the external load, the hysteresis effect, time delay, and the irregular force transmission in the initial phase can be avoided, which is very important for designing a high-accuracy controller of motion systems. Five experiments are performed with the goal of verifying this assumption. In these five experiments, the string is pretwisted five turns before the experiments starting. Figure 16 illustrates the relationship among motor revolutions, force, and linear displacements generated by the twistedstring actuator after pretwisting string. Both the force and linear displacement curves express better continuity than no pretwisted curves, as shown in Figure 16. erefore, pretwisting the string is a method, efficient and easy to achieve, to reduce the complexity of controller design. Although this kinetostatic model described in the last section can reflect the relationship between axial forces, the torque generated by motor, and string length, there are still obvious differences between mathematical modeling and actual performance of a twisted string. For these nonlinear plants, the experiments can be performed to extract approximate model by using sophisticated curve fitting approaches. In this study, the approximation of the static characteristic of the twist string actuator is made from Figure 17 using the Matlab toolbox called Curve Fitting Tool. e armature current of DC motor which is proportional to torque is the input of the motion system in this study. erefore, the function of force transmission and linear displacements are required to design the controller. e approximate function of force transmission and linear displacement are found by the curve fitting approach.
en, the force generated by the twisted-string actuator can be predicted by this approximate function. e approximate function is exponentiation function: e free-body diagram for the load is shown in Figure 2. e equation of motion of the mass m is where m is the mass of the load, F is the tensile force generated by the twisted-string actuator, x is the displacement of the load, F L is the external force, which can be measured by the force sensor, f is the friction among strands, and Δ is unstructured uncertainties of the model. Based on equation (11), we can rewrite equation of motion (12) as follows: where K τ is the relationship between the displacement of the load and the tensile force generated by twisted-string actuator and τ L is the torque generated by the DC motor, which is equal: where K t is the torque constant of the DC motor and I is the armature current of the DC motor. We can build the dynamic model of the twisted-string actuation system in a state-space form as follows: x] is the system state. e traditional velocity and position controller design method focuses on the mechanical dynamics with an inertial load but ignoring the electrical dynamics. With the aim of high tracking performance, it is required to consider the electrical system and the connection between the electrical system and mechanical systems. Form Figure 2, it is obvious that the control input voltage does not directly influence the displacement of the load. To cope with this, the concept of backstepping can be employed to deal with the first equation in (15) not including the control input. We can rewrite the state-space model by integrating the electrical subsystem as follows:

Mathematical Problems in Engineering
where L is the constant rotor inductance, R is the constant rotor resistance, K b is the constant back EMF coefficient, ω is the speed of DC motor, and u is the input control voltage. Define a switch function, which is original from slide mode control, as where e is error between desired trajectory and actual trajectory of external load's displacement, and k 1 positive feedback gain. e transfer function form Z to e is Since the Z(s) is a linear system and k1 is positive, the e(t) will be converging zero if the Z(t) converges to zero. en, the E(s) will be converging to zero if the Z(t) and e(t) are converging to zero. View the state variable I as an independent virtual "input." e first step of the controller design is to ensure this virtual "input" can both stabilizes the second equation in (15) and let the Z(t) eventually goes to zero, which guarantee that the tacking error e(t) converges to zero. Differentiating Z with respect to time leads to By combing the second equation in (11), we have Rewrite (19) as Substitute (11) into (21): Rewrite (22) as where e adaptive robust control mechanism compensate the unknown parameters and structured uncertainties by adjusting parameters online. However, even a small disturbance may lead the parameters adaptive process to instability and large transient tracking error. By integrating the essence of the deterministic robust control (DRC), the adaptive robust controller can attenuate the effect of unstructured uncertainties and improve the transient performance. erefore, the desired virtual input I can be given by I � I a + I s , where I a is adjustable feedforward controller used to compensate the model by parameter adaption, I s is the robust controller in which I s1 is linear feedback to stabilize the nominal model of the twisted-string actuation and I s2 is used to attenuate the parameter estimation error and nonlinear term of model. In equation (23), the vector θ are all unknown parameters that need to be estimated. erefore, the vector θ will be replaced by the θ, the estimated value of the vector θ, in the controller. Since the K τ (x 1 ) and K t are positive function and constant, the Lyapunov candidate function can be defined as follows: where θ is the difference between the real value of the vector θ and θ, the estimated value of the vector. Taking the derivative of (25) gives Since the _ θ � _ θ − _ θ and θ is constant, we have Substituting (23) Since the I a is adjustable feedforward controller used to compensate the model by parameter adaption, the I a is Letting I s1 � KZ, we get where K is a positive feedback gain. In order to ensure the stability of the system, the derivative of V should be less or equal to zero by letting Solving equation (31) gives where erefore, the adaptation mechanisms for adjusting parameters are Assuming the error between the actual virtual input and the desired virtual input is where the I is the actual virtual input and I d is the desired virtual input. Substituting (34) into (23) gives Substituting (29) into (35) gives Rewrite (36) as For the deterministic robust control, the robust term I s2 should be any continuous function which satisfies the following conditions: where ε is a positive design parameter that can be set arbitrarily small. Furthermore, the Lyapunov function can be defined as follows: e derivative of the V s becomes One example of I s2 satisfying equation (38) is where θ m � θ max − θ min , θ max is the upper limit of the θ, and θ min is the lower limit of the θ. Substituting (6) and (10) into (9) gives From the definition of the θ m , the derivative of the V s satisfies In the first case, assuming e 2 � 0, the inequality (44) then becomes For the deterministic robust control, the robust term u s2 should be any continuous function which satisfies the following conditions: (1)e 2 W 1 θ 1 + u s2 ≤ ε 1 , where ε 1 is a positive design parameter that can be set arbitrarily small. Furthermore, the Lyapunov function can be defined as follows: e derivative of the V s1 becomes Substituting (65) into (68) gives One example of u s2 satisfying the (66) is where θ 1m � θ 1max − θ 1min , θ 1max is the upper limit of the θ 1 , and θ 1min is the lower limit of the θ 1 . Substituting (70) into (69) gives From the definition of the θ 1m ,, the derivative of the V s1 satisfies Case 1: if |e 2 ||W 1 θ 1m | > ε 1 , then Case 2: if |e 2 ||W 1 θ 1m | < ε 1 , then In the first case, the derivative of the V s1 satisfies e following inequality then can be derived by integrating (75). at is, where t is time variant and ‖L‖ max is the upper bound of L. According to Bellman-Gronwall, the positive definite V s1 is bounded by erefore, e 2 is bounded by In the second case, the derivative of the V s1 satisfies e following inequality then can be derived by integrating (79). at is, where t is time variant and ‖L‖ max is the upper bound of L. According to Bellman-Gronwall, the positive definite V s1 is bounded by erefore, e 2 is bounded by From the above inequalities, the tracking error of the twisted-string system exponentially approaches to a ball with converging rate K 1 . Within a finite duration, the asymptotic output tracking is achieved. Furthermore, the control input u can be derived as follows: In this section, we present simulation studies to illustrate the performance of the adaptive robust controller based on the force. e following parameters are used: the mass m of the load is 0.5 kg, the stiffness of the spring generating the external force is 4000 N/m, the armature resistant of R of the DC motor is 9.7 ohm, the inductance of the DC motor is 0.0003H, the torque constant Kt of DC motor is 0.1186, and the back EMF constant Kb of DC motor is 0.0611. e controller parameters are designed by selecting k 1 � 10,    K � 3, and K 1 � 7. According to the system properties, the mass m of the load is 0.5 kg, the stiffness of the spring generating the external force is 4000 N/m, and the bounded unknown parameters are specified in Table 3.

Conclusion
In this study, we have proposed a method to estimate the contracted/expanded length of twisted-string actuation. In this study, two types of control structures will be investigated to reduce the effect of the unknown disturbances and the nonlinearity of the twisted-string actuation. ey are (1) LQR controller with the force feedback and (2) ARC controller with the force feedback. e goal of the controller is to compensate for the tracking performances of the twistedstring actuators. Figures 18 and 19 demonstrate the experimental results of tracking performance of the ARC controller for both sinusoidal and trapezoidal motions. For such a sinusoidal motion, the amplitude of the periodic displacement is 60 mm and the duration of each cycle is 12 seconds. e ARC control scheme with force feedback demonstrates excellent results for tracking performance. It is clear that the ARC controller can improve both the transient and steadystate tracking performance over that of the LQR controller from Figure 20. From both Figures 18 and 19, the maximal tracking error of the proposed controller is ∼2 mm. It is clear that either tracking errors of desired trajectory or the difference between the desired virtual input and actual virtual input (armature current of the DC motor) asymptotically approaches to zero, which proves the feasibility of the proposed controller, as shown in Figure 21. Nonlinearity and unknown external disturbance were big issues of twistedstring actuation, which can make the controller design much more difficult. With the proposed methods, the twisted string can demonstrate effective controlled outcome. However, some other methods still need to be developed. An eventtrigger based adaptive control scheme should be developed for the twisted-string actuation systems with unknown external disturbance. Under the event-trigger based adaptive control scheme [10], the effects caused by external disturbance and unknown modeling error can be compensated by constructing the estimator of its unknown constant upper bound.

Data Availability
e data used to support the findings of this study were supplied by the Nature Science Funding of China under rules and so cannot be made freely available. Requests for access to these data should be made to Lei Jiang, leijiang@ cumtb.edu.cn].

Conflicts of Interest
e author declares that they have no conflicts of interest.