Design and Analysis of Novel Actuation Mechanism with Controllable Stiffness

Actuators intended for human–machine interaction systems are usually designed to be mechanically compliant. Conventional actuators are not suitable for this purpose due to typically high stiffness. Advanced powered prosthetic and orthotic devices can vary their stiffness during a motion cycle and are power-efficient. This paper proposes a novel actuator design that modulates stiffness by means of a flexible beam. A motorized drive system varies the active length of the cantilever beam, thus achieving stiffness modulation. New large deflection formulation for cantilever beams with rolling contact constraints is used to determine the moment produced by the actuator. To validate the proposed solution method, an experiment was performed to measure large static deformations of a cantilever beam with the same boundary conditions as in the actuator design. The experiments indicate excellent agreement between measured and calculated contact forces between beam and roller, from which the actuator moment is determined.


Introduction
The overall performance of motion devices found in industry, entertainment, rehabilitation, prosthetics, and general robotics is dictated by the capabilities and efficiency of the driving actuator.Due to the increasing number of ways in which actuators can work alongside humans, safety and energy considerations are pivotal to their development [1].Variable-stiffness actuators (VSA) were introduced to offer a compliant mechanical interface and opportunities for energy optimization through control.Many promising VSA designs have been proposed, with low stiffness and small energy losses due to friction.VSA concepts have also enabled the development of lightweight actuator designs that present advantages for powered human-assist devices, such as orthoses, prostheses, and exoskeletons.Comprehensive reviews of variable-stiffness actuation through the use of compliant elements may be found, for instance, in References [2,3].
VSA designs, and variable impedance actuators more generally, may be classified in a number of ways.In terms of the means to achieve impedance variations, a classification was made in Reference [4] according to the presence of physical springs, dampers, or inertial elements.Active impedance systems do not rely on such physical elements, obtaining variable impedance effects solely by means of motor feedback control.On the other hand, inherent compliance, damping, or inertial designs include the corresponding physical elements fitted with a mechanism to modulate their respective parameters.As in this paper, designs may exist that cut across categories.
In terms of actuator power requirements, those reported in the literature are often classified as active or passive actuators.An active device relies on an external source of power (such as a battery) to operate a motor providing force or torque.On the other hand, a passive device produces forces from the action of spring or damping elements.When a motor or other powered device is used to adjust the damping or stiffness of a physical element, the design is referred to as semiactive or quasi-active.In this case, the required power by the adjustment motor is negligible in comparison with that corresponding to the actuator output.
Any given variable impedance actuator draws characteristics from the above classes.The Mechanical Impedance Adjuster (MIA) of Reference [5] is of the inherent compliance and damping type.The prototype consists of three separate actuators to adjust compliance, position, and damping, and it requires four sensors.In terms of power classification, the MIA is active since an externally powered motor delivers the main actuation forces.Successful applications of the MIA design in single and multiple degrees of freedom (DOF) have been reported [6,7].
An important subclass of inherent compliance designs is given by Series Elastic Actuators (SEA), designed to include an elastic element between the joint and an externally powered actuator.These designs offer fast impact-absorption response and are therefore safe for operation with humans.The spring operates as a mechanical low-pass filter, whereby elevated force peaks are rapidly absorbed by the actuator and not transmitted to the support.Physical damping elements can also be added [4], or active control techniques may be implemented to achieve damping-like behavior [8].SEA design may be optimized with the goal of decreasing energy consumption by proper selection of the compliant element's stiffness [9].Furthermore, optimization techniques may be used to select the motor and the mechanism parameters [10][11][12].
Active impedance control greatly reduces reliance on physical elements and the corresponding complexity [3].This approach has been used and accepted for several decades in industrial robotics [13,14].For instance, in Reference [15], the compliant motion of an industrial robot manipulator is mimicked measuring joint torques and employing them as feedback to manipulate control gains.Recent studies by groups such as the German Aerospace Center (DLR) and KUKA Roboter also feature active compliance designs.The Lightweight Robot (LWR3) [16,17] and other achievements in physical human-robot interaction [18,19] are further examples.
Active compliance, however, presents disadvantages when compared to physical compliance mechanisms.First, the creation of a compliance effect (virtual spring) from an intrinsically noncompliant actuator, such as an electric motor, involves power losses that are not present when energy is stored and released from a physical spring.This is because a positive Joule dissipation power is associated with any current, regardless of its direction.System operation also relies on sensors, which are susceptible to failure, and there is additional complexity associated with the implementation of a feedback control system.Active compliance actuators have been primarily aimed at ensuring safety for humans [20] rather than improving energy efficiency.
Variable impedance designs have entered the field of powered prostheses and human-assist devices.They have presented a variety of advantages for users, such as comfort and energy-consumption reduction.Historically, the majority of prostheses (transtibial and transfemoral) presented fixed or limited stiffness variations during gait motion.Recently, researchers [21] introduced a nonbiologically inspired compliant prosthetic limb to modulate stiffness.A stack of leaf springs was introduced to provide the required compliance.Actuator stiffness is modulated with feedback control, allowing high-speed motion with small energy costs.A similar approach was taken in a human-robot collaboration to perform an isometric bicep curl with load [22].In Reference [23], a novel variable-impedance controller modulated by force sensing was introduced to operate a powered leg prosthesis with energy regeneration.Surplus power associated with portions of the gait cycle is regenerated and stored in supercapacitors.Near-normal gait at net-zero power consumption was demonstrated in level walking by an amputee.The above designs clearly cut across the categories proposed in Reference [4].
Other mechanisms were found in the literature that are capable of stiffness modulation.In References [24,25], an ankle-foot prosthesis with stiffness modulation was proposed.To achieve different ankle stiffnesses during ambulation gait, the prosthetic foot uses cam-based transmission and a variable leaf spring to regulate the torque curve.From the rotation of the ankle, a deflection of the leaf spring is created.The active length of the leaf spring is controlled by a motorized lead screw, which is also responsible for adjusting the stiffness of the spring.This design is semiactive, since a small amount of power is used to adjust stiffness, and the required resisting torque at the ankle is generated without an actuating motor.Simultaneously, the design combines the features of active impedance control and inherent compliance designs.
The VSA concept of this paper can be regarded as a combination of active-impedance and inherent-compliance design principles.Specifically, a physical compliant element is used to ensure the primary requirements of the actuator, such as intrinsic compliance, load bearing and efficient energy storage, and release its capabilities.At the same time, the actuator has controllable features to modulate stiffness and damping.As in the MIA, the new device modulates stiffness by controlling the active length of a leaf spring.However, unlike the MIA, the mechanism is far simpler, backdrivable, and does not rely on friction to regulate damping.The proposed variable-stiffness mechanism is intended for use with a regenerative electromechanical drive system containing an electric storage element such as supercapacitor.Its design allows energy to be stored in the beam and in the supercapacitor and it is based on a flexible beam capable of large deflections, where stiffness is modulated by varying the position of a contact roller that determines the active length of the beam.When integrated with the electromechanical drive, the proposed VSA corresponds to the semiactive category.
Large beam deflections have been studied for at least 70 years [26].Since then, a variety of new of methods have become available to solve this problem with different degrees of complexity and accuracy.Finite-element analysis (FEA) is extremely accurate and applicable to complex loading conditions.However, its computational complexity is the highest [27] compared with other methods, such as the circle-arc method [28] and the elliptic integral method [29].Another approach to solving this problem is to use a Pseudorigid-Body Model (PRBM) [30].The approach of Reference [31] involves the numerical solution of elliptic integrals.The same initial problem formulation is used in this paper.
The remainder of the paper is organized as follows: Section 2 describes the design concept for the proposed actuator; Section 3 presents the static solution for the geometry and large beam deflections; Section 4 presents an experimental validation of the solution method; and Section 5 concludes and presents future work discussion.

Design Concept
The proposed VSA in this paper was designed for use between any two rigid links connected by a rotational DOF.The design is aimed at storing and releasing energy with high efficiency while providing controllability and keeping the weight and complexity small.The actuator provides a controllable joint torque by means of the variable stiffness accomplished with a flexible cantilever beam.A roller is mounted on a pulley-driven shuttle used to vary the point of application of the force on the beam.This changes beam stiffness and consequently controls the joint moment.Due to beam elasticity, the mechanism can store and return potential energy.
A small electric motor and a servo drive with regenerative electronics are used to shift the contact roller.A supercapacitor is used to deliver or recover the electric energy that is converted by the DC motor/generator from the elastic potential energy of the beam.Supercapacitors are better suited as power-storage elements than batteries due to their high power density and ability to accept large charging currents [23,32].The control system for this actuator, and its energy storage and release capabilities, are described in Reference [32].
Figure 1 shows how the proposed VSA could be used in transfemoral prosthesis.A pair of twin swivelling rollers are used to deflect the cantilever beam, and roller position can be shifted by the motor-driven pulley system.As detailed below, a controllable knee moment is obtained by modulating roller position.With proper control, the design offers the opportunity to minimize energy consumption while maintaining normal gait patterns.Since rollers are not preloaded against the beam, we assume that contact occurs between the beam and only one of the rollers at a time.The corresponding reaction force can be decomposed into components perpendicular and parallel to the length of Link 2 (represented as the leg shank in Figure 1).The normal component creates a moment in the joint.The parallel component is a thrust force controlled with the motor-driven belt.When the beam releases energy, negative work associated with the thrust force is regenerated by the servo amplifiers and stored in batteries or supercapacitors.Conversely, power may be drawn from the electric storage to drive the belt, deflect the beam, and store elastic energy Such bidirectional exchange of electric and elastic energy is expected to significantly contribute to overall system efficiency.

Geometry and Beam Deflection Statics
A desired knee-moment profile is obtained by suitably varying the point of application of force on the beam, which is, in turn, determined by knee angle and roller position.In order to study the use of the beam as a compliant mechanism, a cantilever beam, subject to large deflections, is examined.As shown in Figure 2, a co-ordinate frame (ub, vb) was established to analyze deflections.Contact between beam and rollers occurs at single point T. Therefore, normal force F must be perpendicular to both the beam and the roller in contact.Two component forces, nP and P, are defined along the ub and vb axes, respectively.The beam is expected to undergo deflections exceeding the assumptions of elementary bending theory.That is, a given point on the beam is displaced in a plane, not merely in the vb direction.Euler-Bernoulli beam theory with geometric nonlinearity is suitable to analyze such cases whenever beam thickness is small relative to its width and length, so that shear deformations may be neglected.
Modeling large deflections in beams was extensively studied with different levels of complexity and accuracy [31].However, the case under study introduces difficulty in that neither the active length of the beam nor the magnitude and orientation of the applied force are known a priori.Rather, these quantities become unknowns that appear in constraining equations associated with tangency with the roller.Therefore, only roller position xr and the angle between the links are known.Note that the cantilever beam is subjected to single-point force.For this study, this point is considered the end forces used in Reference [31], and it defines the active length of the beam.The following sections describe the proposed solution method that returns the contact-force components, the active beam length, and the end deflections and other variables of interest, such as the joint moment and elastic energy stored in the beam.

Large Beam Deflection Formulation
The initial formulation of Reference [31] is used here, followed by a new approach to handle the unknown active length and applied force and orientation.The force applied at contact point T, as shown in Figure 2, is decomposed along the (ub,uv) frame.These forces produce a moment at the origin of the frame that is related to beam geometry by the Euler-Bernoulli moment-curvature relationship:

dt) M(xb,yb) = El-= P(a-x)+ nP(b -y)
(1) where 6 and s are the beam slope and arc length, E is the elastic module, and I is the moment of inertia of the cross section.Additional geometric constraints due to the cantilever configuration determine a boundary-value problem: located at angle £ measured from vr, as shown in Figure 3.The roller is located at constant offset dr from the us axis and at a variable distance xr from the origin along the us axis.The joint angle is denoted as 0; (angle between Links 1 and 2, assumed positive in flexion).The fixed end of the beam is at distances d0 and d1 from the joint, at a slant mounting angle αw, as shown in Figure 1.
Note that the slope of the beam at T relative to ub is given by tan(θ0).Letting φ be the angle between the axis and the resultant force F, a straightforward derivation shows that the rolling contact constraint can be written as: (p -Oo = n/1 (6) tan(</>) = -1/rz (7) Moreover, it is straightforward to show that A = 0, and the following equation applies: To relate end deflections a and b appearing in Equations ( 4) and ( 5) to the above quantities, a set of homogeneous coordinate transformations are defined following the notational conventions of robotics [33].The transformation between frames b and s can be written as: where Tr contains the coordinates of T in the roller frame, namely, (rsin(z),rcos(z)).A symbolic computation of the above transformations yields the pertinent expressions for Tb = (a, b): where Equation ( 8) has been used to eliminate £.The above equations contain the unknown θ0, which must be determined from the beam-deflection relationships presented earlier.

Numerical Solution Method
Given xr and θj, the governing equations can be obtained by Algorithm 1. First, observe that the ratio b/a from Equations ( 4) and ( 5) is a function of θ0 alone.That is: where n = n(θ0) from the rolling contact constraints (Equations ( 6) and ( 7)).Therefore, using this relationship and Equations ( 10) and ( 11), a single nonlinear equation in θ0 is obtained, which can be solved using a tool such as Matlab's fsolve or a custom line search.Additionally, parameters a and b can be found by this process.Beam active length L is solved for by evaluating the integrals presented in Equation ( 4) or ( 5) (one may be used for verification) using the values of a and b.Nondimensional parameter α0 and force P is determined by evaluating the integral in Equation ( 3) with the value found of θ0.
This method yields two functions of interest for actuator control: Moment at the joint due the deflection of the beam Mb(θj, xr) and thrust force Fb(θj,xr) on the actuation mechanism mounted on link 2. Given a set of physical parameters for the beam and the links, both functions can be readily evaluated by a computer program and used in control system development.

Experimental Validation
To validate the modeling approach and solution method used in this paper, a test rig was built to measure large beam deflections under the same boundary conditions as in the actuator design.As shown in Figures 4 and 5 the test rig consisted of two links connected by a joint and a cantilever beam that was fixed at Link 2. The roller can be moved to any desired position and locked.The beam is attached as a cantilever to an L-shaped bracket that freely rotates about the joint.A cable is attached to the top arm of the L-shaped bracket and routed through a pulley.A deadweight is then applied to the end of the cable, which rotates the bracket and produces a deflection against the roller.

Experiment Methodology
A camera and image-processing software were used to capture the shape of the beam under a known load.The experiment was performed with four distinct rectangular beams with different dimensions and properties as shown in Table 1.Euler-Bernoulli theory remains accurate for sufficiently slender beams, that is, beams whose thickness is much smaller than their width and length [34].Although a universal threshold for this ratio is not available for all situations, a lower limit of 5.4 for the length-to-thickness ratio has been mentioned [35] for cantilever beams.Table 2 shows the ratios for the active lengths involved in the experiments, largely meeting this threshold.The validity of the Euler-Bernoulli assumptions is best supported by the accuracy of the experiment predictions.
The elastic modulus of each beam was estimated using free-oscillation experiments [36] and direct measurements of cross-section properties and beam density.A set of known weights were applied to the cable, with the roller position locked.As a result, the beam underwent deflection, resulting in a reaction force from roller F. Dimensions a, b, and e were extracted by image processing.Beam shape was estimated by fitting a third-degree polynomial function to image points.Other geometric properties of the beam were estimated by linear fitting and by direct measurements from the test rig.From Figure 6, contact force P can be written as the following equation: p = mge + fa + famgsinfa + (yd + rfjmgcosfa + i/j + 0 + tm The above equation contains the mass of Link 1, intended to account for the gravitational moment due to the weight of Link 1.However, this influence was subsequently determined to be insignificant.Aiming to validate Equation ( 13), the roller structure was replaced by a load cell (Omega Engineering, Norwalk, CT, USA).The difference between experiment results and actual force P was small across all tests.Specifically, an overall root mean square (RMS) of 1.956 N was found using four beams and three materials according to the following formula: where N is the total number of trials and Pn is the difference between forces Pexp extracted experimentally and forces calculated by the proposed algorithm.A total of N = 56 trials were performed, with contact forces varying between 2 and 20 N. The accuracy of the proposed approach can also be evaluated using nondimensional parameter α0 from Equation (15).For this, theoretical reaction P was found for each α0 calculated with experimental data.These theoretical forces were compared with predictions using the algorithm and forces extracted from the experiment using image processing and static analysis, as described earlier.
The results of these comparisons are presented in Figures 7 and 8 for L = 0.750 m and L = 0.430 m, respectively.In each figure, solid lines are calculated with the proposed algorithm.RMS errors between measurements and calculations for each beam were 1.1525, 2.9869, 0.92602, and 1.8255 N, indicating excellent agreement.Beam dimensions were selected to cover a wide range of nondimensional parameter α0.This ensures that the validity of the proposed calculation method is not restricted to the relatively small dimensional range of the specimens used in the experiment.It is expected that slender beams with arbitrary dimensions under loading magnitudes matching the range of α0 result in similarly accurate predictions.

Conclusions and Future Work
A new design for an actuation mechanism with controllable stiffness was presented.A formulation of large beam deflections with a rolling contact constraint was developed, and a numerical solution algorithm was established.Experimental procedures using a custom test rig validated the formulation and solution algorithm.
A marked advantage of the proposed design over previous VSA concepts is that the used mechanism to vary stiffness does not incur significant power losses, but rather is designed to convert and store the energy released by the beam in electric form by means of a regenerative drive system containing a storage element such as a supercapacitor.Moreover, damping effects can be introduced by means of control without the intervention of lossy physical damping elements.
The proposed mechanism uses a straight beam that may reduce the applicability of the concept in prosthetic legs due to a range of motion limitations.The cantilever beam may exhibit high root stresses and protrude inconveniently when large link rotations are attempted.One way to mitigate this limitation is to consider curved or variable-section beams that can be optimally designed.The method of this paper would require substantial modifications to allow for such variable geometries.The proposed VSA may be more conveniently used as an orthotic device around human joints having smaller ranges of motion, for instance, for ankle or wrist dorsiflexion.
The proposed calculation method requires the solution of a nonlinear equation in which numerical integrations are involved.However, model-based real-time control algorithms may require faster computations.A possible solution is to use the proposed method offline to generate moment data for a range of roller positions and link angles.The control algorithm may then access the solutions in the form of a look-up table.
Further work related to the VSA of this paper may involve optimization of the shape and material of the beam, including the possibility of using functionally graded designs.Direct nonlinear control techniques could be explored that achieve impedance regulation in combination with energy-optimization methods.

Figure 1 .
Figure 1.(a) General schematic of the proposed variable-stiffness actuator (VSA) design in a prosthetic leg application.(b) Geometry used for mechanism analysis.

Figure 2 .
Figure 2. Geometry of beam deflection.Active length L is the arc length between fixed end and point T.
an,Transx_dlTransy _doRotZ/e.and between frames b and frames r: Hs r(u) = Transx,XrTransy/dr The explicit form of the above transformation matrices is found in Appendix A.l.The coordinates of T expressed in the b-frame can now be found as Tb = HbTs = HbHf.Tr (9)

Figure 4 .
Figure 4. Schematic of test-rig structure.Cable is used to apply forces using deadweights.

Figure 5 .
Figure 5.An instance of beam-deflection measurements with the test rig.

Figure 6 .
Figure 6.Geometry and loading diagram of the test rig used to derive contact force between roller and beam.

Figure 7 .
Figure 7. Contacted force P versus nondimensional parameter for beam length of 0.75 m.

Figure 8 .
Figure 8. Contacted force P versus nondimensional parameter for beam length of 0.43 m.

Table 2 .
Width-to-length and thickness-to-length ratios for the two active lengths used in the experiments (l1 = 0.75 m and l2 = 0.43 m).