Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device
Identifieur interne : 001B88 ( Pmc/Checkpoint ); précédent : 001B87; suivant : 001B89Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device
Auteurs : Jakob Oblak [Slovénie] ; Zlatko Matja I [Slovénie]Source :
- Journal of NeuroEngineering and Rehabilitation [ 1743-0003 ] ; 2011.
Abstract
Variable structure parallel mechanisms, actuated with low-cost motors with serially added elasticity (series elastic actuator - SEA), has considerable potential in rehabilitation robotics. However, reflected masses of a SEA and variable structure parallel mechanism linked with a compliant actuator result in a potentially unstable coupled mechanical oscillator, which has not been addressed in previous studies.
The aim of this paper was to investigate through simulation, experimentation and theoretical analysis the necessary conditions that guarantee stability and passivity of a haptic device (based on a variable structure parallel mechanism driven by SEA actuators) when in contact with a human. We have analyzed an equivalent mechanical system where a dissipative element, a mechanical damper was placed in parallel to a spring in SEA.
The theoretical analysis yielded necessary conditions relating the damping coefficient, spring stiffness, both reflected masses, controller's gain and desired virtual impedance that needs to be fulfilled in order to obtain stable and passive behavior of the device when in contact with a human. The validity of the derived passivity conditions were confirmed in simulations and experimentally.
These results show that by properly designing variable structure parallel mechanisms actuated with SEA, versatile and affordable rehabilitation robotic devices can be conceived, which may facilitate their wide spread use in clinical and home environments.
Url:
DOI: 10.1186/1743-0003-8-3
PubMed: 21251299
PubMed Central: 3063208
Affiliations:
Links toward previous steps (curation, corpus...)
Links to Exploration step
PMC:3063208Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device</title><author><name sortKey="Oblak, Jakob" sort="Oblak, Jakob" uniqKey="Oblak J" first="Jakob" last="Oblak">Jakob Oblak</name><affiliation wicri:level="1"><nlm:aff id="I1">University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana, Slovenia</nlm:aff><country xml:lang="fr">Slovénie</country><wicri:regionArea>University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana</wicri:regionArea><wicri:noRegion>1000 Ljubljana</wicri:noRegion></affiliation></author><author><name sortKey="Matja I, Zlatko" sort="Matja I, Zlatko" uniqKey="Matja I Z" first="Zlatko" last="Matja I">Zlatko Matja I</name><affiliation wicri:level="1"><nlm:aff id="I1">University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana, Slovenia</nlm:aff><country xml:lang="fr">Slovénie</country><wicri:regionArea>University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana</wicri:regionArea><wicri:noRegion>1000 Ljubljana</wicri:noRegion></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">21251299</idno><idno type="pmc">3063208</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3063208</idno><idno type="RBID">PMC:3063208</idno><idno type="doi">10.1186/1743-0003-8-3</idno><date when="2011">2011</date><idno type="wicri:Area/Pmc/Corpus">002043</idno><idno type="wicri:Area/Pmc/Curation">002043</idno><idno type="wicri:Area/Pmc/Checkpoint">001B88</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device</title><author><name sortKey="Oblak, Jakob" sort="Oblak, Jakob" uniqKey="Oblak J" first="Jakob" last="Oblak">Jakob Oblak</name><affiliation wicri:level="1"><nlm:aff id="I1">University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana, Slovenia</nlm:aff><country xml:lang="fr">Slovénie</country><wicri:regionArea>University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana</wicri:regionArea><wicri:noRegion>1000 Ljubljana</wicri:noRegion></affiliation></author><author><name sortKey="Matja I, Zlatko" sort="Matja I, Zlatko" uniqKey="Matja I Z" first="Zlatko" last="Matja I">Zlatko Matja I</name><affiliation wicri:level="1"><nlm:aff id="I1">University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana, Slovenia</nlm:aff><country xml:lang="fr">Slovénie</country><wicri:regionArea>University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana</wicri:regionArea><wicri:noRegion>1000 Ljubljana</wicri:noRegion></affiliation></author></analytic><series><title level="j">Journal of NeuroEngineering and Rehabilitation</title><idno type="eISSN">1743-0003</idno><imprint><date when="2011">2011</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p>Variable structure parallel mechanisms, actuated with low-cost motors with serially added elasticity (series elastic actuator - SEA), has considerable potential in rehabilitation robotics. However, reflected masses of a SEA and variable structure parallel mechanism linked with a compliant actuator result in a potentially unstable coupled mechanical oscillator, which has not been addressed in previous studies.</p></sec><sec><title>Methods</title><p>The aim of this paper was to investigate through simulation, experimentation and theoretical analysis the necessary conditions that guarantee stability and passivity of a haptic device (based on a variable structure parallel mechanism driven by SEA actuators) when in contact with a human. We have analyzed an equivalent mechanical system where a dissipative element, a mechanical damper was placed in parallel to a spring in SEA.</p></sec><sec><title>Results</title><p>The theoretical analysis yielded necessary conditions relating the damping coefficient, spring stiffness, both reflected masses, controller's gain and desired virtual impedance that needs to be fulfilled in order to obtain stable and passive behavior of the device when in contact with a human. The validity of the derived passivity conditions were confirmed in simulations and experimentally.</p></sec><sec><title>Conclusions</title><p>These results show that by properly designing variable structure parallel mechanisms actuated with SEA, versatile and affordable rehabilitation robotic devices can be conceived, which may facilitate their wide spread use in clinical and home environments.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Harwin, Ws" uniqKey="Harwin W">WS Harwin</name></author><author><name sortKey="Patton, J" uniqKey="Patton J">J Patton</name></author><author><name sortKey="Edgerton, Vr" uniqKey="Edgerton V">VR Edgerton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hidler, Jm" uniqKey="Hidler J">JM Hidler</name></author><author><name sortKey="Nichols, D" uniqKey="Nichols D">D Nichols</name></author><author><name sortKey="Pelliccio, M" uniqKey="Pelliccio M">M Pelliccio</name></author><author><name sortKey="Brady, K" uniqKey="Brady K">K Brady</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Krebs, Hi" uniqKey="Krebs H">HI Krebs</name></author><author><name sortKey="Palazzolo, Jj" uniqKey="Palazzolo J">JJ Palazzolo</name></author><author><name sortKey="Dipietro, L" uniqKey="Dipietro L">L Dipietro</name></author><author><name sortKey="Ferraro, M" uniqKey="Ferraro M">M Ferraro</name></author><author><name sortKey="Krol, J" uniqKey="Krol J">J Krol</name></author><author><name sortKey="Rannekleiv, K" uniqKey="Rannekleiv K">K Rannekleiv</name></author><author><name sortKey="Volpe, Bt" uniqKey="Volpe B">BT Volpe</name></author><author><name sortKey="Hogan, N" uniqKey="Hogan N">N Hogan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pignolo, L" uniqKey="Pignolo L">L Pignolo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Krebs, Hi" uniqKey="Krebs H">HI Krebs</name></author><author><name sortKey="Hogan, N" uniqKey="Hogan N">N Hogan</name></author><author><name sortKey="Volpe, Bt" uniqKey="Volpe B">BT Volpe</name></author><author><name sortKey="Aisen, Ml" uniqKey="Aisen M">ML Aisen</name></author><author><name sortKey="Edelstein, L" uniqKey="Edelstein L">L Edelstein</name></author><author><name sortKey="Diels, C" uniqKey="Diels C">C Diels</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Coote, S" uniqKey="Coote S">S Coote</name></author><author><name sortKey="Stokes, E" uniqKey="Stokes E">E Stokes</name></author><author><name sortKey="Murphy, B" uniqKey="Murphy B">B Murphy</name></author><author><name sortKey="Harwin, W" uniqKey="Harwin W">W Harwin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Der Linde, Rq" uniqKey="Van Der Linde R">RQ Van der Linde</name></author><author><name sortKey="Lammertse, P" uniqKey="Lammertse P">P Lammertse</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Burgar, Cg" uniqKey="Burgar C">CG Burgar</name></author><author><name sortKey="Lum, Ps" uniqKey="Lum P">PS Lum</name></author><author><name sortKey="Shor, Pc" uniqKey="Shor P">PC Shor</name></author><author><name sortKey="Van Der, Loos" uniqKey="Van Der L">Loos Van der</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hesse, S" uniqKey="Hesse S">S Hesse</name></author><author><name sortKey="Schulte Tigges, G" uniqKey="Schulte Tigges G">G Schulte-Tigges</name></author><author><name sortKey="Konrad, M" uniqKey="Konrad M">M Konrad</name></author><author><name sortKey="Bardeleben, A" uniqKey="Bardeleben A">A Bardeleben</name></author><author><name sortKey="Werner, C" uniqKey="Werner C">C Werner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nef, T" uniqKey="Nef T">T Nef</name></author><author><name sortKey="Mihelj, M" uniqKey="Mihelj M">M Mihelj</name></author><author><name sortKey="Kiefer, G" uniqKey="Kiefer G">G Kiefer</name></author><author><name sortKey="Perndl, C" uniqKey="Perndl C">C Perndl</name></author><author><name sortKey="Muller, R" uniqKey="Muller R">R Müller</name></author><author><name sortKey="Riener, R" uniqKey="Riener R">R Riener</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sanchez, Rj" uniqKey="Sanchez R">RJ Sanchez</name></author><author><name sortKey="Wolbrecht, E" uniqKey="Wolbrecht E">E Wolbrecht</name></author><author><name sortKey="Smith, R" uniqKey="Smith R">R Smith</name></author><author><name sortKey="Liu, J" uniqKey="Liu J">J Liu</name></author><author><name sortKey="Cramer, S" uniqKey="Cramer S">S Cramer</name></author><author><name sortKey="Rahman, T" uniqKey="Rahman T">T Rahman</name></author><author><name sortKey="Bobrow, Je" uniqKey="Bobrow J">JE Bobrow</name></author><author><name sortKey="Reinkensmeyer, Dj" uniqKey="Reinkensmeyer D">DJ Reinkensmeyer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stienen, A" uniqKey="Stienen A">A Stienen</name></author><author><name sortKey="Hekman, E" uniqKey="Hekman E">E Hekman</name></author><author><name sortKey="Van Der Helm, Fct" uniqKey="Van Der Helm F">FCT Van der Helm</name></author><author><name sortKey="Prange, Gb" uniqKey="Prange G">GB Prange</name></author><author><name sortKey="Jannink, Mja" uniqKey="Jannink M">MJA Jannink</name></author><author><name sortKey="Aalsma, Amm" uniqKey="Aalsma A">AMM Aalsma</name></author><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H Van der Kooij</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pratt, G" uniqKey="Pratt G">G Pratt</name></author><author><name sortKey="Williamson, Mm" uniqKey="Williamson M">MM Williamson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robinson, W" uniqKey="Robinson W">W Robinson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wyeth, G" uniqKey="Wyeth G">G Wyeth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Oblak, J" uniqKey="Oblak J">J Oblak</name></author><author><name sortKey="Cikajlo, I" uniqKey="Cikajlo I">I Cikajlo</name></author><author><name sortKey="Matja, Z" uniqKey="Matja Z">Z Matjačć</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Oblak, J" uniqKey="Oblak J">J Oblak</name></author><author><name sortKey="Perry, Jc" uniqKey="Perry J">JC Perry</name></author><author><name sortKey="Jung, Jh" uniqKey="Jung J">JH Jung</name></author><author><name sortKey="Cikajlo, I" uniqKey="Cikajlo I">I Cikajlo</name></author><author><name sortKey="Keller, T" uniqKey="Keller T">T Keller</name></author><author><name sortKey="Matja I, Z" uniqKey="Matja I Z">Z Matjaćić</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Massie, Th" uniqKey="Massie T">TH Massie</name></author><author><name sortKey="Salisbury, Jk" uniqKey="Salisbury J">JK Salisbury</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Williamson, Mm" uniqKey="Williamson M">MM Williamson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vallery, H" uniqKey="Vallery H">H Vallery</name></author><author><name sortKey="Veneman, J" uniqKey="Veneman J">J Veneman</name></author><author><name sortKey="Van Asseldonk, E" uniqKey="Van Asseldonk E">E Van Asseldonk</name></author><author><name sortKey="Ekkelenkamp, R" uniqKey="Ekkelenkamp R">R Ekkelenkamp</name></author><author><name sortKey="Buss, M" uniqKey="Buss M">M Buss</name></author><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H Van der Kooij</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Colgate, Je" uniqKey="Colgate J">JE Colgate</name></author><author><name sortKey="Schenkel, G" uniqKey="Schenkel G">G Schenkel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Colgate, Je" uniqKey="Colgate J">JE Colgate</name></author><author><name sortKey="Grafing, Pe" uniqKey="Grafing P">PE Grafing</name></author><author><name sortKey="Stanley, Mc" uniqKey="Stanley M">MC Stanley</name></author><author><name sortKey="Schenkel, G" uniqKey="Schenkel G">G Schenkel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Colgate, E" uniqKey="Colgate E">E Colgate</name></author><author><name sortKey="Hogan, N" uniqKey="Hogan N">N Hogan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hogan, N" uniqKey="Hogan N">N Hogan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gil, Jj" uniqKey="Gil J">JJ Gil</name></author><author><name sortKey="Avello, A" uniqKey="Avello A">A Avello</name></author><author><name sortKey="Rubio, A" uniqKey="Rubio A">Á Rubio</name></author><author><name sortKey="Fl Rez, J" uniqKey="Fl Rez J">J Flórez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Diolaiti, N" uniqKey="Diolaiti N">N Diolaiti</name></author><author><name sortKey="Niemeyer, G" uniqKey="Niemeyer G">G Niemeyer</name></author><author><name sortKey="Barbagli, F" uniqKey="Barbagli F">F Barbagli</name></author><author><name sortKey="Salisbury, Jk" uniqKey="Salisbury J">JK Salisbury</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">J Neuroeng Rehabil</journal-id><journal-title-group><journal-title>Journal of NeuroEngineering and Rehabilitation</journal-title></journal-title-group><issn pub-type="epub">1743-0003</issn><publisher><publisher-name>BioMed Central</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">21251299</article-id><article-id pub-id-type="pmc">3063208</article-id><article-id pub-id-type="publisher-id">1743-0003-8-3</article-id><article-id pub-id-type="doi">10.1186/1743-0003-8-3</article-id><article-categories><subj-group subj-group-type="heading"><subject>Methodology</subject></subj-group></article-categories><title-group><article-title>Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes" id="A1"><name><surname>Oblak</surname><given-names>Jakob</given-names></name><xref ref-type="aff" rid="I1">1</xref><email>jakob.oblak@ir-rs.si</email></contrib><contrib contrib-type="author" id="A2"><name><surname>Matjačić</surname><given-names>Zlatko</given-names></name><xref ref-type="aff" rid="I1">1</xref><email>zlatko.matjacic@ir-rs.si</email></contrib></contrib-group><aff id="I1"><label>1</label>University Rehabilitation Institute, Republic of Slovenia, Linhartova 51, 1000 Ljubljana, Slovenia</aff><pub-date pub-type="collection"><year>2011</year></pub-date><pub-date pub-type="epub"><day>20</day><month>1</month><year>2011</year></pub-date><volume>8</volume><fpage>3</fpage><lpage>3</lpage><history><date date-type="received"><day>14</day><month>5</month><year>2010</year></date><date date-type="accepted"><day>20</day><month>1</month><year>2011</year></date></history><permissions><copyright-statement>Copyright ©2011 Oblak and Matjačć; licensee BioMed Central Ltd.</copyright-statement><copyright-year>2011</copyright-year><copyright-holder>Oblak and Matjačć; licensee BioMed Central Ltd.</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0">http://creativecommons.org/licenses/by/2.0</ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p></license></permissions><self-uri xlink:href="http://www.jneuroengrehab.com/content/8/1/3"></self-uri><abstract><sec><title>Background</title><p>Variable structure parallel mechanisms, actuated with low-cost motors with serially added elasticity (series elastic actuator - SEA), has considerable potential in rehabilitation robotics. However, reflected masses of a SEA and variable structure parallel mechanism linked with a compliant actuator result in a potentially unstable coupled mechanical oscillator, which has not been addressed in previous studies.</p></sec><sec><title>Methods</title><p>The aim of this paper was to investigate through simulation, experimentation and theoretical analysis the necessary conditions that guarantee stability and passivity of a haptic device (based on a variable structure parallel mechanism driven by SEA actuators) when in contact with a human. We have analyzed an equivalent mechanical system where a dissipative element, a mechanical damper was placed in parallel to a spring in SEA.</p></sec><sec><title>Results</title><p>The theoretical analysis yielded necessary conditions relating the damping coefficient, spring stiffness, both reflected masses, controller's gain and desired virtual impedance that needs to be fulfilled in order to obtain stable and passive behavior of the device when in contact with a human. The validity of the derived passivity conditions were confirmed in simulations and experimentally.</p></sec><sec><title>Conclusions</title><p>These results show that by properly designing variable structure parallel mechanisms actuated with SEA, versatile and affordable rehabilitation robotic devices can be conceived, which may facilitate their wide spread use in clinical and home environments.</p></sec></abstract></article-meta></front><body><sec><title>Background</title><p>Rehabilitation robotics is a rapidly evolving field [<xref ref-type="bibr" rid="B1">1</xref>-<xref ref-type="bibr" rid="B4">4</xref>]. Numerous haptic robots were developed for movement training of upper extremities following neurological disorder. According to works published to date, robots for upper extremity motor rehabilitation are usually serial linkage mechanisms that can be in general divided in two groups. The first consists of serial linkage mechanisms with only 1 to 3 degrees of freedom (DOF), where the end-effector of the robot is in contact with the user's hand, making it suitable for only one activity of upper extremity movement (either arm reaching movement or wrist movements). Clinical use of such low-DOF serial mechanisms [<xref ref-type="bibr" rid="B5">5</xref>-<xref ref-type="bibr" rid="B10">10</xref>] necessitates the use of two or more devices in order to provide comprehensive upper extremity movement training. This is neither convenient from a practical nor cost-effective point of view. On the other hand, exoskeleton mechanisms may have up to 7 DOF [<xref ref-type="bibr" rid="B11">11</xref>-<xref ref-type="bibr" rid="B13">13</xref>] and can provide comprehensive upper extremity movement rehabilitation. However, such mechanisms require high quality back drivable actuators for each DOF, which necessitates complex and thus expensive design.</p><p>Few rehabilitation robots have implemented a parallel kinematic structure. Parallel mechanisms usually have mechanical linkages with many DOF that greatly exceed the resulting DOF of the whole mechanism. This property allows for a design where some of the joints may be easily locked or unlocked, thus resulting in different workspace configurations suitable for different aspects of arm or wrist movements training. Another characteristic of parallel robots is that the actuators are located at the robot's base. This feature allows the implementation of series elastic actuators (SEA) [<xref ref-type="bibr" rid="B14">14</xref>-<xref ref-type="bibr" rid="B16">16</xref>] that utilize standard off-the-shelf DC motors with suitable planetary gearheads and suitable springs, providing similar overall haptic performance as their high quality back-drivable counterparts. Universal Haptic Drive (UHD) [<xref ref-type="bibr" rid="B17">17</xref>] and Variable Structure Pantograph (VSP) [<xref ref-type="bibr" rid="B18">18</xref>] are the two devices in which mechanisms with lockable joints and SEA actuation were successfully implemented and tested in clinical practice.</p><p>However, from a control point of view, both beneficial aspects; parallel kinematic structure (such as VSP) and SEA based drive, result in a mechanical system where the reflected masses of the SEA and the parallel kinematic structure (serially connected with a spring) become comparable, resulting in a coupled mechanical oscillator. Suitable control of such a rehabilitation robot, which should provide stable haptic interaction when in contact with a human, may present a considerable challenge, not addressed in previous studies.</p><p>The aim of this paper was to investigate theoretically, through simulations and experimentally the necessary conditions that guarantee stability and passivity of a haptic device, based on a variable structure parallel mechanism driven by SEA actuators, when in contact with a human. We have analyzed an equivalent mechanical system where a dissipative element, a mechanical damper, was placed in parallel with a spring in SEA. The goal was to derive conditions that must be met in order to allow use of a SEA driven variable structure parallel mechanism as a stable haptic interface in upper extremity rehabilitation.</p></sec><sec sec-type="methods"><title>Methods</title><sec><title>Variable structure pantograph: mechanical linkage</title><p>Figure <xref ref-type="fig" rid="F1">1</xref> shows a VSP haptic device that is composed of a variable structure parallel linkage and two SEA actuators. A brief description on the device is provided here, while more detailed information can be found in [<xref ref-type="bibr" rid="B17">17</xref>,<xref ref-type="bibr" rid="B18">18</xref>].</p><fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>The photograph of variable structure pantograph (VSP)</bold>. The essence of the VSP is variable structure parallel mechanism, which is driven by a visco-elastic actuator. The VSP promises high suitability for training of upper extremity tasks involving hand positioning and orientation.</p></caption><graphic xlink:href="1743-0003-8-3-1"></graphic></fig><p>The main parts of the parallel mechanism are the three joints that can be either locked or released (Figure <xref ref-type="fig" rid="F2">2(A)</xref>). By locking or releasing DOFs in these joints (I, II, III), the mechanical configuration of VSP can be changed, enabling use of the device in several operational modes shown in Figure <xref ref-type="fig" rid="F2">2</xref>. The three possible modes ("ARM", "WRIST" and "REACH") are briefly described below.</p><fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>A variable structure parallel mechanism</bold>. A variable structure parallel mechanism enables using the device in different operational modes. Switching between modes can be easily achieved by locking or releasing joints I, II, III (A). Workspace in the "ARM" mode (B), "WRIST" mode (C) and "REACH" mode (D) are presented.</p></caption><graphic xlink:href="1743-0003-8-3-2"></graphic></fig><p>"ARM" mode: locking joint I and releasing joints II and III, results in 2 DOF quasi-planar movements in Forward/Backward/Left/Right directions, as shown in Figure <xref ref-type="fig" rid="F2">2(B)</xref>. The movement prescribed by the workspace "ARM" mode is similar to required for reaching and/or moving objects on a table, desk, or countertop.</p><p>"WRIST" mode: the mechanical configuration, termed as "WRIST" mode, is achieved by releasing joint I and locking joints II and III. A subject holding on the handle bar can perform movements in wrist as shown in Figure <xref ref-type="fig" rid="F2">2(C)</xref>. By setting the offset orientation of the handle bar in the horizontal or vertical position, movement of all 3-DOFs in wrist (Flexion/Extension, Radial/Ulnar deviation and Pronation/Supination as shown in Figure <xref ref-type="fig" rid="F2">2(C)</xref>) can be achieved. The resulting movement of the user's wrist is similar to what is required for performing wrist-orienting motions in the following activities: pouring from a bottle, brushing teeth, or stirring a pot.</p><p>"REACH" mode: locking joints I and III and releasing joint II results in a mechanical configuration, which allows training of Forward and Up/Lateral reach movements. These motions are therefore similar to activities such as reaching for a high drawer or cupboard, or moving objects from one side of the cupboard to the other.</p></sec><sec><title>Variable structure pantograph: series visco-elastic actuation</title><p>The variable structure parallel linkage of VSP is actuated by a SEA based drive as shown in Figure <xref ref-type="fig" rid="F1">1</xref>. The implemented drive consists of two sets of DC motors (Maxon, RE40, 150 W) with gearheads (GP 52 C, 81:1) and incremental encoders. Torques from both gearheads' shafts impose force on the actuated bar through serially added mechanical springs and string wires, see Figure <xref ref-type="fig" rid="F3">3(A, B)</xref>. The string wires are connected to the actuated bar perpendicular one to another, which enables actuation of the VSP in 2DOFs.</p><fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Series visco-elastic actuator and parallel mechanism</bold>. (A) Actuation of the VSP consists of: 1-two sets of DC motors with gears and encoders, 2-elastic springs, 3- mechanical dampers, 4-pulleys and 5-an actuated bar. Schematic presentation of the series visco-elastic actuator and parallel mechanism (only 1 DOF is shown for clarity), with characteristic values of mechanical component parameters used in the VSP. Impedance felt at the arm is 16 times smaller than at the bottom of the actuated bar, because the actuated bar is divided by a spherical joint in ratio 4:1.</p></caption><graphic xlink:href="1743-0003-8-3-3"></graphic></fig><p>Introduction of an elastic element (mechanical spring) in series with the motor provides many benefits, including: more accurate and stable force control, attenuation of both backlash and friction nonlinear effects and the actuators' own impedance as well as providing greater shock tolerance (important for safety concerns). On the other hand, introducing SEA in haptic drive leads to reduction of the achievable force bandwidth. Since relatively slow movements can be expected during rehabilitation training, reduced force bandwidth does not present a significant problem.</p><p>Utilization of a variable structure parallel mechanism is essential in designing a versatile rehabilitation device. However, using a parallel mechanism where considerable endpoint mass is in series with both SEA's spring and motor mass, results in a coupled oscillator needing appropriate damping. Adequate dissipation of mechanical energy is needed, to achieve a stable haptic interaction when the device is in contact with a human. A convenient location for a mechanical damper may be in parallel with the SEA spring. Figure <xref ref-type="fig" rid="F3">3(A, B)</xref> presents a schematic illustration of an implemented parallel mechanism driven by a series visco-elastic actuator.</p></sec><sec><title>Variable structure pantograph: linearized model</title><p>In Figure <xref ref-type="fig" rid="F4">4</xref>, the open loop model is illustrated, where M and m denotes masses, X<sub>I </sub>and X<sub>o </sub>positions, and F<sub>I </sub>and F<sub>O </sub>forces on the motor and the actuated bar, respectively. Attaching parallel mechanism on the actuated bar, significantly increases endpoint mass. The motor is coupled to the parallel mechanism via a mechanical spring K and damper b. The equivalent viscous friction in the motor and planetary gearhead is marked with B [<xref ref-type="bibr" rid="B17">17</xref>].</p><fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>A linearized model</bold>. The motor with a gearhead (reflected mass M) is connected to the variable structure parallel mechanism (reflected mass m) via a spring (compliance K) and a viscous damper (damping B).</p></caption><graphic xlink:href="1743-0003-8-3-4"></graphic></fig><p>The relation between motor mass and the mass of the parallel mechanism can be given by two differential equations:</p><p><disp-formula id="bmcM1"><label>(1)</label><mml:math id="M1" name="1743-0003-8-3-i1" overflow="scroll"><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>I</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>M</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>¨</mml:mo></mml:mover><mml:mtext>I</mml:mtext></mml:msub><mml:mo>−</mml:mo><mml:mtext>B</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>I</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>K</mml:mtext><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>I</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>b</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>I</mml:mtext></mml:msub><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math></disp-formula></p><p><disp-formula id="bmcM2"><label>(2)</label><mml:math id="M2" name="1743-0003-8-3-i2" overflow="scroll"><mml:mrow><mml:mtext>K</mml:mtext><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>I</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>K</mml:mtext><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>I</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>b</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:mtext>m</mml:mtext><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>¨</mml:mo></mml:mover><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>−</mml:mo><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext> </mml:mtext><mml:mo>=</mml:mo><mml:mn>0.</mml:mn></mml:mrow></mml:math></disp-formula></p><p>By taking Laplace transforms and relating equations, an expression relating F<sub>O</sub>, X<sub>O </sub>and F<sub>I </sub>can be found:</p><p><disp-formula id="bmcM3"><label>(3)</label><mml:math id="M3" name="1743-0003-8-3-i3" overflow="scroll"><mml:mtable><mml:mtr><mml:mtd><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>m</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>+</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>I</mml:mtext></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p><p>This equation is important, for it determines the motor torque/force F<sub>I</sub>, needed to achieve a given output force F<sub>O</sub>, when the handle bar of the VSP is moving. If the parallel mechanism of the device is assumed to be clamped <inline-formula><mml:math id="M4" name="1743-0003-8-3-i4" overflow="scroll"><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>O</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula>, then the transfer function between output force and motor force is:</p><p><disp-formula id="bmcM4"><label>(4)</label><mml:math id="M5" name="1743-0003-8-3-i5" overflow="scroll"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>I</mml:mtext></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext></mml:mrow><mml:mrow><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula></p><p>On the other hand, the position to force transfer function of the uncontrolled plant (F<sub>I </sub>= 0) is defined as:</p><p><disp-formula id="bmcM5"><label>(5)</label><mml:math id="M6" name="1743-0003-8-3-i6" overflow="scroll"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>m</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula></p><p>The negative sign before X<sub>O </sub>comes from the definition of the directions of F<sub>O </sub>and X<sub>O</sub>. These two equations define the model of the plant to be controlled. The motion of the handle bar (X<sub>O</sub>) is modelled as a disturbance on the output force (F<sub>O</sub>), see shaded block in Figure <xref ref-type="fig" rid="F5">5</xref>.</p><fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Impedance based force control</bold>. The inner force loop compares desired the virtual force F<sub>V </sub>with the force feedback F<sub>O</sub>. The output impedance loop calculates the desired virtual force F<sub>V </sub>from the position of the handle bar X<sub>O</sub>.V is the desired virtual stiffness and P the proportional gain of the controller.</p></caption><graphic xlink:href="1743-0003-8-3-5"></graphic></fig><p>The principal purpose of haptic devices is to allow human operators to touch, feel and manipulate objects in a virtual environment. For this reason, the impedance felt at the handle bar should be as close as possible to the desired virtual impedance (V, Figure <xref ref-type="fig" rid="F5">5</xref>).</p><p>Usually, three criteria are employed when designing haptic devices [<xref ref-type="bibr" rid="B19">19</xref>]: (1) movement in free space (LOW impedance) should be opposed with minimal possible force, (2) solid virtual objects (HIGH impedance) must feel stiff, and (3) virtual constraints must not be easily saturated, which requires a suitable impedance- based force control. In series elastic actuators, a variety of control strategies are possible. Williamson [<xref ref-type="bibr" rid="B20">20</xref>] proposed a control strategy for SEA with a feed-forward model and PID controller. Vallery [<xref ref-type="bibr" rid="B21">21</xref>] chose the concept of cascade force control with proportional-integral controller. We decided to implement the simplest approach, which is a proportional controller, in order to have a clearer picture on the influence of various mechanical components on passivity of haptic interaction [<xref ref-type="bibr" rid="B22">22</xref>,<xref ref-type="bibr" rid="B23">23</xref>].</p><p>Impedance based force control (Figure <xref ref-type="fig" rid="F5">5</xref>), was implemented in MATLAB (Simulink). In computer simulation, the VSP's haptic performance was investigated by simulating sinusoidal movements of the handle bar (X<sub>O</sub>). In simulation, LOW and HIGH impedance virtual force F<sub>V </sub>was compared to calculated force on the handle bar F<sub>O</sub>, for different values of damper b and proportional gain P. In order to investigate the option of using low-cost motors with potentially redundant backlash, different values of backlash were considered in the simulation model.</p></sec><sec><title>Conditions regarding stability and passivity of the system</title><p>For a haptic system, the output impedance is usually defined as the transfer function from the velocity of the gripper, in our case handle bar <inline-formula><mml:math id="M7" name="1743-0003-8-3-i7" overflow="scroll"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mover accent="true"><mml:mtext>X</mml:mtext><mml:mo>˙</mml:mo></mml:mover><mml:mtext>O</mml:mtext></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>, to the opposing force (F<sub>O</sub>):</p><p><disp-formula id="bmcM6"><label>(6)</label><mml:math id="M8" name="1743-0003-8-3-i8" overflow="scroll"><mml:mrow><mml:mtext>Z</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext>s</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext>s</mml:mtext></mml:mrow></mml:mfrac><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p><p>Colgate [<xref ref-type="bibr" rid="B24">24</xref>] and Hogan [<xref ref-type="bibr" rid="B25">25</xref>] have proven that a system will be stable while in contact with changing environments if and only if the output impedance Z(s), obeys the following rules:</p><p>1. Z(s) has no poles in the right half plane</p><p>2. Re(Z(jw)) is nonnegative for all frequencies w</p><p>If these conditions are met, the impedance is a stable function of frequency and the system exhibits passivity.</p><p>By means of Equation 6 and simplified system model in Figure <xref ref-type="fig" rid="F5">5</xref>, output impedance of our controlled system is given by:</p><p><disp-formula id="bmcM7"><label>(7)</label><mml:math id="M9" name="1743-0003-8-3-i9" overflow="scroll"><mml:mtable><mml:mtr><mml:mtd><mml:mtext>Z</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext>s</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext>s</mml:mtext></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>m</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext>V</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>P</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>s</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>P</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p><p>First, we will check the condition of asymptotic stability. The characteristic equation for output impedance can be written as:</p><p><disp-formula id="bmcM8"><label>(8)</label><mml:math id="M10" name="1743-0003-8-3-i10" overflow="scroll"><mml:mrow><mml:msup><mml:mrow><mml:mtext>Ms</mml:mtext></mml:mrow><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b P</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>K</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K P</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mtext>s</mml:mtext><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math></disp-formula></p><p>It is obvious that Hurwitz determinants of the characteristic equation for Z(s) are nonnegative for all technically realizable values of mechanical components. Henceforth, Z(s) has no poles in right half plane and the first rule is met.</p><p>By replacing the complex variable "s" in Equation 7 by "jw", the frequency response Z(jw) can be obtained. The real part of the impedance frequency response is given by:</p><p><disp-formula id="bmcM9"><label>(9)</label><mml:math id="M11" name="1743-0003-8-3-i11" overflow="scroll"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mtext>Re</mml:mtext><mml:mo></mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>Z</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext>j</mml:mtext><mml:mtext>w</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mtext>Re</mml:mtext><mml:mo></mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>O</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mtext>X</mml:mtext><mml:mtext>O</mml:mtext></mml:msub><mml:mtext>j</mml:mtext><mml:mtext>w</mml:mtext></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mtext>z</mml:mtext><mml:mn>0</mml:mn><mml:mo>+</mml:mo><mml:mtext>z</mml:mtext><mml:mn>2</mml:mn><mml:msup><mml:mtext>w</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext>z</mml:mtext><mml:mn>4</mml:mn><mml:msup><mml:mtext>w</mml:mtext><mml:mn>4</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>b</mml:mtext><mml:mo>+</mml:mo><mml:mtext>B</mml:mtext><mml:mo>+</mml:mo><mml:mtext>b</mml:mtext><mml:mtext>P</mml:mtext><mml:mo stretchy="false">)</mml:mo><mml:mtext>w</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mtext>K</mml:mtext><mml:mo>+</mml:mo><mml:mtext>K</mml:mtext><mml:mtext>P</mml:mtext><mml:mo>−</mml:mo><mml:mtext>M</mml:mtext><mml:msup><mml:mtext>w</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p><p>Re(Z(jw)) is nonnegative for all frequencies w, if all values z0, z2 and z4 are nonnegative (see Figure <xref ref-type="fig" rid="F6">6</xref>). This gives three conservative conditions for the passivity of the system that have to be considered.</p><fig id="F6" position="float"><label>Figure 6</label><caption><p><bold>Conservative conditions regarding passivity of the system</bold>.</p></caption><graphic xlink:href="1743-0003-8-3-6"></graphic></fig><p>First, the virtual stiffness V is limited by the stiffness of the mechanical spring K and controller's proportional gain P. However, if we look at the third condition, P is limited by the reflected motor mass M and reflected mass of the parallel mechanism m. It is straightforward, that value P can be increased by reducing m. The second condition demands that there must be sufficient damping between the relative position of motor X<sub>I </sub>and parallel mechanism X<sub>O</sub>. Usually, in technical realization of the system, damping is always present due to viscous friction in mechanical components but is not sufficient. For this reason, an additional damping element should be inserted in parallel with the spring to satisfy the damping condition for passivity of the system. The derived conservative conditions for passivity are general. In the following section, these conditions will be applied to characteristic values of the mechanical components used in technical realization of the VSP (listed in Figure <xref ref-type="fig" rid="F3">3(B)</xref>).</p></sec></sec><sec><title>Results</title><sec><title>Variable structure pantograph: application of derived passivity conditions</title><p>In the "ARM" and "REACH" mode of the VSP, estimated reflected mass of the parallel mechanism is relatively high (m = 18 kg). From the third condition on passivity, achievable controller's proportional gain is relatively small <inline-formula><mml:math id="M12" name="1743-0003-8-3-i12" overflow="scroll"><mml:mrow><mml:mo mathvariant="bold" stretchy="false">(</mml:mo><mml:mi mathvariant="bold">P</mml:mi><mml:mo mathvariant="bold">≤</mml:mo><mml:mfrac><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="bold">m</mml:mi></mml:mfrac><mml:mo mathvariant="bold">≅</mml:mo><mml:mn mathvariant="bold">1</mml:mn><mml:mo mathvariant="bold">.</mml:mo><mml:mn mathvariant="bold">9</mml:mn><mml:mo mathvariant="bold" stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>. This is due to the high reflected mass of the parallel mechanism m. From the first condition, maximal virtual stiffness V can be determined as <inline-formula><mml:math id="M13" name="1743-0003-8-3-i13" overflow="scroll"><mml:mrow><mml:mi mathvariant="bold">V</mml:mi><mml:mo mathvariant="bold">≤</mml:mo><mml:mi mathvariant="bold">K</mml:mi><mml:mfrac><mml:mrow><mml:mi mathvariant="bold">P</mml:mi><mml:mo mathvariant="bold">+</mml:mo><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mi mathvariant="bold">P</mml:mi></mml:mfrac><mml:mo mathvariant="bold">≅</mml:mo><mml:mn mathvariant="bold">1</mml:mn><mml:mn mathvariant="bold">2</mml:mn><mml:mn mathvariant="bold">0</mml:mn><mml:mn mathvariant="bold">0</mml:mn><mml:mn mathvariant="bold">0</mml:mn><mml:mi mathvariant="bold"> </mml:mi><mml:mi mathvariant="bold"> </mml:mi><mml:mi mathvariant="bold">N</mml:mi><mml:mo mathvariant="bold">/</mml:mo><mml:mi mathvariant="bold">m</mml:mi></mml:mrow></mml:math></inline-formula>. Finally, from the second condition, damping of <inline-formula><mml:math id="M14" name="1743-0003-8-3-i14" overflow="scroll"><mml:mrow><mml:mtext mathvariant="bold">b</mml:mtext><mml:mo>≥</mml:mo><mml:mfrac><mml:mrow><mml:mo>−</mml:mo><mml:msup><mml:mtext mathvariant="bold">B</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">M</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">V</mml:mtext><mml:mo>+</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msup><mml:mtext mathvariant="bold">B</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>−</mml:mo><mml:mtext mathvariant="bold">M</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">V</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>−</mml:mo><mml:mn mathvariant="bold">4</mml:mn><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">K</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">m</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msqrt></mml:mrow><mml:mrow><mml:mn mathvariant="bold">2</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac><mml:mo>≅</mml:mo><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mn mathvariant="bold">8</mml:mn><mml:mn mathvariant="bold">0</mml:mn><mml:mn mathvariant="bold">0</mml:mn><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">N</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">s</mml:mtext><mml:mo>/</mml:mo><mml:mtext mathvariant="bold">m</mml:mtext></mml:mrow></mml:math></inline-formula> is needed, if we set virtual stiffness as <bold>V = 12000 N/m </bold>(solid virtual objects). On the other hand, if we want to emulate free space, where <bold>V = 0 N/m </bold>and all the other parameters remain the same, a much smaller damping of <bold>b ≅ 250 Ns/m </bold>satisfies the second condition for passivity.</p><p>In the "WRIST" mode, parameters needed to meet conditions of passivity are different. Given that almost all mass of the parallel mechanism is supported by joints II and III and therefore does not move (see Figure <xref ref-type="fig" rid="F2">2</xref>) and that only the segment with the handle bar is moving, the reflected mass is much smaller (m = 10 kg). For this reason values that satisfy conditions of passivity in the "WRIST" mode are different; see Table <xref ref-type="table" rid="T1">1</xref>.</p><table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Values of P, V and b that satisfy conditions regarding passivity.</p></caption><table frame="hsides" rules="groups"><thead><tr><th></th><th align="center" colspan="2">"ARM" and "REACH" mode [m = 18 kg]</th><th align="center" colspan="2">"WRIST" mode [m = 10 kg]</th></tr></thead><tbody><tr><td align="center">P</td><td align="center" colspan="2">1,9</td><td align="center" colspan="2">3,5</td></tr><tr><td colspan="5"><hr></hr></td></tr><tr><td align="center">V [N/m]</td><td align="center" colspan="2">0 ≤ V ≤ 12000</td><td align="center" colspan="2">0 ≤ V ≤ 10000</td></tr><tr><td colspan="5"><hr></hr></td></tr><tr><td></td><td align="center">V = 0 N/m</td><td align="center">V = 12000 N/m</td><td align="center">V = 0 N/m</td><td align="center">V = 10000 N/m</td></tr><tr><td></td><td colspan="4"><hr></hr></td></tr><tr><td align="center">b [Ns/m]</td><td align="center">250</td><td align="center">800</td><td align="center">210</td><td align="center">800</td></tr></tbody></table></table-wrap><p>It is obvious from Figure <xref ref-type="fig" rid="F6">6</xref> that it is more demanding to meet conditions of passivity in the case where the end point mass (reflected mass of parallel mechanism m) is higher. For this reason, in further analysis higher mass of parallel mechanism (m = 18 kg in the "ARM" and "REACH" modes) was considered.</p><p>Results listed in Table <xref ref-type="table" rid="T1">1</xref> present conservative conditions for passivity of the VSP, where all values z0, z2, and z4 are nonnegative. However, Re(Z(jw)) can be nonnegative for all frequencies w and desired V, also with suitable selection of P and b (see Figure <xref ref-type="fig" rid="F7">7</xref>). The maximal achievable virtual stiffness V in a HIGH impedance environment is 12000 N/m, which is sufficient for rehabilitation purposes. As can be seen from Figure <xref ref-type="fig" rid="F7">7</xref>, passive VSP behavior in any mode can be achieved (Re(Z(jw)) ≥ 0), if P ≤ 19 and parallel damping is at least b ≥ 780 Ns/m.</p><fig id="F7" position="float"><label>Figure 7</label><caption><p><bold>System can exhibit passivity with suitable selection of P and b, for desired V</bold>. If proportional gain of the controller P is higher than 1.9, Re(Z(jw)) becomes negative and therefore the system does not exhibit passivity. When emulating a LOW impedance environment (V = 0 N/m), damping of at least b ≥ 190 Ns/m is needed, while for a HIGH impedance environment, damping of at least b ≥ 780 Ns/m is needed.</p></caption><graphic xlink:href="1743-0003-8-3-7"></graphic></fig><p>When considering technical realization, a parallel damper with damping coefficient of b ≥ 780 Ns/m would be a rather heavy duty mechanical element. For this reason, it was decided to set b ≈ 200 Ns/m, which can be technically realized, however at the expense of reduced controller's proportional gain P. By demanding VSP' passivity in a HIGH impedance virtual environment (V = 12000 N/m) where parallel damping is b = 200 Nm/s, P should not exceed a value of 0.95 (see Figure <xref ref-type="fig" rid="F8">8(A)</xref>). By reducing the virtual stiffness V, proportional gain P can be increased (Figure <xref ref-type="fig" rid="F8">8(B)</xref>, while maintaining VSP passivity.</p><fig id="F8" position="float"><label>Figure 8</label><caption><p><bold>By reducing virtual stiffness V, proportional gain P can be increased</bold>. By reducing the controller' proportional gain to P = 0.95, parallel damping of b = 200 Ns/m is sufficient for the VSP' passivity when emulating a HIGH impedance environment (V = 12000 N/m) (a). However, P can be increased if a lower virtual stiffness V is desired (b).</p></caption><graphic xlink:href="1743-0003-8-3-8"></graphic></fig></sec><sec><title>Variable structure pantograph: Simulation evaluation</title><p>Based on the results obtained in the previous subsection, simulation evaluation of VSP's haptic performance was undertaken (MATLAB, Simulink). In simulation model, a damper with b = 200 Ns/m was added parallel to the spring as depicted in Figure <xref ref-type="fig" rid="F3">3(B)</xref> and Figure <xref ref-type="fig" rid="F4">4</xref>. In terms of system passivity, the proportional gain of the controller P was varied from 1.9 in a LOW impedance to P = 0.95 in a HIGH impedance simulated environment. Haptic performance was investigated by simulating sinusoidal movements of the handle bar for ± 8 cm at frequencies of 1.0 Hz, 0.5 Hz and 0.1 Hz. It is important to point out that due to the design of the VSP (see Figure <xref ref-type="fig" rid="F3">3(B)</xref>), the displacement on the bottom of the actuated bar (X<sub>O</sub>) is 4 times smaller than the movement of the handle bar. Similarly, the force on the handle bar is 4 times smaller than the force on the bottom of the actuated bar where the cable wire is attached. For this reason, impedance felt by the subject holding the handle bar is 16 times smaller than on the bottom of the actuated bar. Therefore, the impedance felt by the user on the handle bar in a HIGH impedance simulated environment should be approximately 750 N/m (12000 N/m: 16) and the maximal force while repeating sinusoidal movements with amplitude ± 8 cm should be approximately 60 N (750 N/m * 8 cm). Desired force felt by the user in a LOW impedance simulated environment (0 N/m) should be 0 N. Additionally, the influence of backlash (1 mm and 4 mm), which is typically present in DC motors with planetary gears, was investigated.</p><p>Results of the VSP's haptic performance simulation with a parallel added damper (b = 200 Ns/m) are presented in Figure <xref ref-type="fig" rid="F9">9</xref>. The values of the forces presented in Figure <xref ref-type="fig" rid="F8">8</xref> are interaction forces between the user and the handle bar and are approximately 4 times smaller than on the bottom of the actuated bar.</p><fig id="F9" position="float"><label>Figure 9</label><caption><p><bold>Simulation evaluation</bold>. Simulation of VSP performance in the ARM mode with a damper added parallel to the spring (b = 200 Ns/m) in LOW and HIGH impedance environments and for two values of simulated backlash (1 mm and 4 mm). Haptic performance was investigated by simulating sinusoidal movements (1.0 Hz, 0.5 Hz and 0.1 Hz respectively) of the handle bar X<sub>O</sub>, where the desired force (F<sub>V</sub>), was compared to the actual force (F<sub>O</sub>) on the handle bar.</p></caption><graphic xlink:href="1743-0003-8-3-9"></graphic></fig></sec><sec><title>Variable structure pantograph: experimental evaluation</title><p>To verify results of the theoretical analysis and the simulation evaluation in previous subsections, the VSP's haptic performance was also experimentally examined on the recently developed VSP haptic robot [<xref ref-type="bibr" rid="B18">18</xref>]. Parallel to the spring, a damper with b = 200 Ns/m was added. The damper was technically realized by means of a pneumatic cylinder (FESTO), where damping was adjusted by setting appropriate air flow on the input and the output of the cylinder. Similar to simulation, the subject holding the handle bar imposed quasi-sinusoidal movements in the forward and backward direction for ± 8 cm at a frequency of approximately 1 Hz. Haptic performance was examined in the ARM mode by simulating LOW and HIGH impedance environments (see Figure <xref ref-type="fig" rid="F10">10</xref>). In a LOW impedance environment (V = 0 N/m), the actual interaction force between the user and the handle bar did not exceed 5 N, which is comparable to results from computer simulation. On the other hand, when simulating a HIGH impedance environment (V = 12000 N/m), the actual force corresponded to the desired values more precisely than predicted by simulation.</p><fig id="F10" position="float"><label>Figure 10</label><caption><p><bold>Experimental evaluation of VSP performance in ARM mode with damper added parallel to spring (b = 200 Ns/m)</bold>. Desired force is a product of the current position of the handle bar and desired virtual stiffness (V), divided by 16 due to the corresponding leverage implemented in the VSP design.</p></caption><graphic xlink:href="1743-0003-8-3-10"></graphic></fig><p>The passivity of the VSP was also tested experimentally. In order to destabilize haptic interaction between the human and the VSP, the handle bar was exposed to fast movements, without and with an added damper parallel to the spring. Haptic interaction in both cases is presented in Figure <xref ref-type="fig" rid="F11">11</xref>. It is obvious from Figure <xref ref-type="fig" rid="F11">11</xref>, that in the case when there is no damper added parallel to the spring, response of the VSP to fast movements (which normally do not occur in the rehabilitation movement) is much more oscillatory than in the case where a damper is added. In the oscillating interval, the subject is still holding the handle bar and thereby prevents the VSP from experiencing an unstable response. If the subjects were to release the handle bar during the oscillating interval, the VSP's response would become unstable, which could potentially lead to mechanical destruction of the device. This demonstrates that the VSP has stable behavior, if parallel to the spring, sufficient damping is presented. That was not the case when damper was omitted. Quantitative comparisons between simulations and experiments would not give meaningful results since the movement in the experimental evaluation is human-driven and therefore highly variable.</p><fig id="F11" position="float"><label>Figure 11</label><caption><p><bold>VSP movement and interaction forces induced by fast movements of the handle bar</bold>. Adding a damping element (b = 200 Ns/m) parallel to the mechanical spring significantly stabilizes haptic performance.</p></caption><graphic xlink:href="1743-0003-8-3-11"></graphic></fig></sec></sec><sec><title>Discussion</title><p>Actuators with series elasticity have been extensively studied in the field of robotics [<xref ref-type="bibr" rid="B14">14</xref>-<xref ref-type="bibr" rid="B16">16</xref>,<xref ref-type="bibr" rid="B20">20</xref>], where they were predominantly used in actuation of walking machines. Use of these actuators in haptic devices was limited to cases where the endpoint mass of devices are negligible as compared to the reflected inertia of the actuator [<xref ref-type="bibr" rid="B21">21</xref>], (included references reflect only a portion of the relevant literature). In case of the Variable Structure Pantograph haptic device, endpoint mass is not negligible due to a variable structure parallel mechanism. The main contribution of this paper is the derivation of passivity conditions that need to be fulfilled for a rehabilitation robot with a mechanism mass comparable to the reflected mass of the geared actuator. The results show that appropriate damping must be provided parallel to the SEA spring in order to obtain stable and passive behavior of the device when it is in contact with a human.</p><p>Three necessary conditions in terms of passivity were obtained from the theoretical analysis, which should be used when designing haptic devices with SEA actuation:</p><p>1.) The maximum proportional gain P of the controller must be limited by the ratio of the actuator and parallel mechanism reflected masses: <inline-formula><mml:math id="M15" name="1743-0003-8-3-i15" overflow="scroll"><mml:mrow><mml:mi mathvariant="bold">P</mml:mi><mml:mo mathvariant="bold">≤</mml:mo><mml:mfrac><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="bold">m</mml:mi></mml:mfrac></mml:mrow></mml:math></inline-formula>. Hence, better force control can be achieved (higher P) either by use of a motor with higher reflected inertia or by use of lighter parallel mechanism.</p><p>2.) The maximum achievable virtual stiffness V must be limited by controller's proportional gain P and the stiffness of mechanical spring K added in series to the motor: <inline-formula><mml:math id="M16" name="1743-0003-8-3-i16" overflow="scroll"><mml:mrow><mml:mi mathvariant="bold">V</mml:mi><mml:mo mathvariant="bold">≤</mml:mo><mml:mi mathvariant="bold">K</mml:mi><mml:mfrac><mml:mrow><mml:mi mathvariant="bold">P</mml:mi><mml:mo mathvariant="bold">+</mml:mo><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mi mathvariant="bold">P</mml:mi></mml:mfrac></mml:mrow></mml:math></inline-formula>. This condition is similar to results obtained by Vallery [<xref ref-type="bibr" rid="B21">21</xref>], where it was shown that the SEA cannot display higher pure stiffness than the spring stiffness when passivity is desired.</p><p>3.) Third, to achieve haptic device passivity, sufficient damping b should be presented in parallel to the SEA spring: <inline-formula><mml:math id="M17" name="1743-0003-8-3-i17" overflow="scroll"><mml:mrow><mml:mtext mathvariant="bold">b</mml:mtext><mml:mo>≥</mml:mo><mml:mfrac><mml:mrow><mml:mo>−</mml:mo><mml:msup><mml:mtext mathvariant="bold">B</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">M</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">V</mml:mtext><mml:mo>+</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msup><mml:mtext mathvariant="bold">B</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>−</mml:mo><mml:mtext mathvariant="bold">M</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">V</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn mathvariant="bold">2</mml:mn></mml:msup><mml:mo>−</mml:mo><mml:mn mathvariant="bold">4</mml:mn><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">K</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">m</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msqrt></mml:mrow><mml:mrow><mml:mn mathvariant="bold">2</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mo>+</mml:mo><mml:mtext mathvariant="bold">B</mml:mtext><mml:mtext> </mml:mtext><mml:mtext mathvariant="bold">P</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:math></inline-formula>. Necessity of appropriate damping was also derived by Colgate and Schenkel [<xref ref-type="bibr" rid="B22">22</xref>], where the passivity of systems comprising a continuous time plant and discrete time controller was considered. This means that a damper inserted parallel to the spring ensures required dissipation of mechanical energy.</p><p>To verify results of theoretical analysis, simulation evaluation was undertaken. Simulation results predicted stable haptic performance for both HIGH and LOW impedance simulated environments. Simulation results revealed that haptic performance is also adequate in the case when higher values of backlash are presented in the system. This means that high-cost precision motors and gearheads that are currently used in VSP haptic device may be replaced by low-cost motors with greater backlash. Experimental evaluation has confirmed the simulation results and has shown that when appropriate damping and controller's proportional gain are used, stable interaction between machine and human are achieved in LOW and HIGH impedance environments, which was not the case when the damper was omitted. Generation of a HIGH impedance environment is limited to a virtual stiffness of 750 N/m, because the impedance felt at the arm is 16 times smaller than what the actuator can provide (12000 N/m:16). However, this does not present a notable limitation to rehabilitation where more compliant and thus gentle guiding in performance of training tasks is necessary. Also, experimental evaluation revealed that the achievable impedance range is sufficient [<xref ref-type="bibr" rid="B17">17</xref>].</p><p>In most cases described in literature, discrete linear models are used when dealing with general purpose sampled haptic environments, which are characterized with high Z-bandwidth [<xref ref-type="bibr" rid="B26">26</xref>,<xref ref-type="bibr" rid="B27">27</xref>]. It is important to point out that in this paper we utilized a continuous linear model of the studied haptic robot. In the particular case of the rehabilitation robot actuated with SEA, the typical Z-bandwidth is much lower (in our case the virtual stiffness is limited to 750 N/m). Also, the force bandwidth of the SEA as well as movement during rehabilitation are limited to app. 1 Hz [<xref ref-type="bibr" rid="B17">17</xref>], while the sampling rate is relatively high (1 kHz). Furthermore, it has been demonstrated that the effects of digitalization in conjunction with a usually high Z-bandwidth, (that a haptic interface should be able to render) can cause instabilities at frequencies of several hundred Hz, while at frequencies below 10 Hz, the effects of A/D and D/A devices placed within the control loop are negligible [<xref ref-type="bibr" rid="B27">27</xref>]. This enabled the use of a continuous linear model, which is much more intuitive to comprehend. The decision to model the parallel mechanism with a simple mass is related to the fact that the range of motion of the VSP is rather limited and relatively slow, meaning that the predominant dynamics will be dominated by the mass properties of the mechanism. Finally, the use of a linear model to mimic the dynamics of a geared DC motor has been experimentally validated in our previous paper describing the UHD robot [<xref ref-type="bibr" rid="B17">17</xref>].</p></sec><sec><title>Conclusions</title><p>In conclusion the results of our study have shown that by properly designing rehabilitation device that uses a parallel mechanism and actuators with series elasticity, stability and passivity of haptic performance can be obtained. Because such a haptic system may be composed entirely of off-the-shelf mechanical components, versatile and affordable rehabilitation robotic devices can be produced, which may facilitate their wide spread use in clinical and home environments.</p></sec><sec><title>Competing interests</title><p>University Rehabilitation Institute, Republic of Slovenia, received financial support for taking part in the development of the Variable Structure Pantograph device from FATRONIK Tecnalia, Spain. Both authors are co-authors of a patent application describing essential features of the Variable Structure Pantograph device while FATRONIK Tecnalia is the assignee of the same patent application.</p></sec><sec><title>Authors' contributions</title><p>Both authors significantly contributed to the conception, theoretical analysis, simulation and experimental evaluation and writing of the manuscript. Both authors revised and approved the final manuscript.</p></sec></body><back><sec><title>Acknowledgements</title><p>The authors acknowledge the financial support from FATRONIK Tecnalia and Slovenian research agency (grant no. P2-0228).</p></sec><ref-list><ref id="B1"><mixed-citation publication-type="journal"><name><surname>Harwin</surname><given-names>WS</given-names></name><name><surname>Patton</surname><given-names>J</given-names></name><name><surname>Edgerton</surname><given-names>VR</given-names></name><article-title>Challenges and opportunities for robot mediated neurorehabilitation</article-title><source>Proceedings of the IEEE Special issue on medical robotics</source><year>2006</year><volume>94</volume><issue>9</issue><fpage>1717</fpage><lpage>1726</lpage></mixed-citation></ref><ref id="B2"><mixed-citation publication-type="journal"><name><surname>Hidler</surname><given-names>JM</given-names></name><name><surname>Nichols</surname><given-names>D</given-names></name><name><surname>Pelliccio</surname><given-names>M</given-names></name><name><surname>Brady</surname><given-names>K</given-names></name><article-title>Advances in the understanding and treatment of stroke impairment using robotic devices</article-title><source>Top Stroke Rehabil</source><year>2005</year><volume>12</volume><issue>2</issue><fpage>22</fpage><lpage>35</lpage><pub-id pub-id-type="doi">10.1310/RYT5-62N4-CTVX-8JTE</pub-id><pub-id pub-id-type="pmid">15940582</pub-id></mixed-citation></ref><ref id="B3"><mixed-citation publication-type="journal"><name><surname>Krebs</surname><given-names>HI</given-names></name><name><surname>Palazzolo</surname><given-names>JJ</given-names></name><name><surname>Dipietro</surname><given-names>L</given-names></name><name><surname>Ferraro</surname><given-names>M</given-names></name><name><surname>Krol</surname><given-names>J</given-names></name><name><surname>Rannekleiv</surname><given-names>K</given-names></name><name><surname>Volpe</surname><given-names>BT</given-names></name><name><surname>Hogan</surname><given-names>N</given-names></name><article-title>Rehabilitation robotics: Performance-based progressive robot-assisted therapy</article-title><source>Autonomous Robots</source><year>2003</year><volume>15</volume><issue>1</issue><fpage>7</fpage><lpage>20</lpage><pub-id pub-id-type="doi">10.1023/A:1024494031121</pub-id></mixed-citation></ref><ref id="B4"><mixed-citation publication-type="journal"><name><surname>Pignolo</surname><given-names>L</given-names></name><article-title>Robotics in neuro-rehabilitation</article-title><source>J Rehab Med</source><year>2009</year><volume>41</volume><fpage>955</fpage><lpage>960</lpage><pub-id pub-id-type="doi">10.2340/16501977-0434</pub-id></mixed-citation></ref><ref id="B5"><mixed-citation publication-type="journal"><name><surname>Krebs</surname><given-names>HI</given-names></name><name><surname>Hogan</surname><given-names>N</given-names></name><name><surname>Volpe</surname><given-names>BT</given-names></name><name><surname>Aisen</surname><given-names>ML</given-names></name><name><surname>Edelstein</surname><given-names>L</given-names></name><name><surname>Diels</surname><given-names>C</given-names></name><article-title>Overview of clinical trials with MIT-MANUS: a robot-aided neuro-rehabilitation facility</article-title><source>Technology and health care: official journal of the European Society for Engineering and Medicine</source><year>1999</year><volume>7</volume><issue>6</issue><fpage>419</fpage><lpage>423</lpage><pub-id pub-id-type="pmid">10665675</pub-id></mixed-citation></ref><ref id="B6"><mixed-citation publication-type="book"><name><surname>Coote</surname><given-names>S</given-names></name><name><surname>Stokes</surname><given-names>E</given-names></name><name><surname>Murphy</surname><given-names>B</given-names></name><name><surname>Harwin</surname><given-names>W</given-names></name><article-title>The effect of GENTLE/s robot-mediated therapy on upper extremity dysfunction post stroke</article-title><source>Proceedings of the 8th International Conference on Rehabilitation Robotics 23-25 April 2003</source><publisher-name>Daejon, Korea</publisher-name><fpage>59</fpage><lpage>61</lpage></mixed-citation></ref><ref id="B7"><mixed-citation publication-type="journal"><name><surname>Van der Linde</surname><given-names>RQ</given-names></name><name><surname>Lammertse</surname><given-names>P</given-names></name><article-title>HapticMaster-a generic force controlled robot for human interaction</article-title><source>Industrial Robot: An International Journal</source><year>2003</year><volume>30</volume><issue>6</issue><fpage>515</fpage><lpage>524</lpage><pub-id pub-id-type="doi">10.1108/01439910310506783</pub-id></mixed-citation></ref><ref id="B8"><mixed-citation publication-type="other"><article-title>REO Therapy</article-title><ext-link ext-link-type="uri" xlink:href="http://www.motorika.com">http://www.motorika.com</ext-link></mixed-citation></ref><ref id="B9"><mixed-citation publication-type="journal"><name><surname>Burgar</surname><given-names>CG</given-names></name><name><surname>Lum</surname><given-names>PS</given-names></name><name><surname>Shor</surname><given-names>PC</given-names></name><name><surname>Van der</surname><given-names>Loos</given-names></name><article-title>Development of robots for rehabilitation therapy: The Palo Alto VA/Stanford experience</article-title><source>Journal Rehabil Res Dev</source><year>2000</year><volume>37</volume><issue>6</issue><fpage>663</fpage><lpage>673</lpage></mixed-citation></ref><ref id="B10"><mixed-citation publication-type="journal"><name><surname>Hesse</surname><given-names>S</given-names></name><name><surname>Schulte-Tigges</surname><given-names>G</given-names></name><name><surname>Konrad</surname><given-names>M</given-names></name><name><surname>Bardeleben</surname><given-names>A</given-names></name><name><surname>Werner</surname><given-names>C</given-names></name><article-title>Robot-assisted arm trainer for the passive and active practice of bilateral forearm and wrist movements in hemiparetic subjects</article-title><source>Arch Phys Med Rehabil</source><year>2003</year><volume>84</volume><issue>6</issue><fpage>915</fpage><lpage>920</lpage><pub-id pub-id-type="doi">10.1016/S0003-9993(02)04954-7</pub-id><pub-id pub-id-type="pmid">12808550</pub-id></mixed-citation></ref><ref id="B11"><mixed-citation publication-type="other"><name><surname>Nef</surname><given-names>T</given-names></name><name><surname>Mihelj</surname><given-names>M</given-names></name><name><surname>Kiefer</surname><given-names>G</given-names></name><name><surname>Perndl</surname><given-names>C</given-names></name><name><surname>Müller</surname><given-names>R</given-names></name><name><surname>Riener</surname><given-names>R</given-names></name><article-title>ARMin - Exoskeleton for Arm Therapy in Stroke Patients</article-title><source>Proceedings of the 2007 IEEE 10th International Conference on Rehabilitation Robotics</source><year>2007</year><fpage>68</fpage><lpage>74</lpage></mixed-citation></ref><ref id="B12"><mixed-citation publication-type="other"><name><surname>Sanchez</surname><given-names>RJ</given-names><suffix>Jr</suffix></name><name><surname>Wolbrecht</surname><given-names>E</given-names></name><name><surname>Smith</surname><given-names>R</given-names></name><name><surname>Liu</surname><given-names>J</given-names></name><name><surname>Cramer</surname><given-names>S</given-names></name><name><surname>Rahman</surname><given-names>T</given-names></name><name><surname>Bobrow</surname><given-names>JE</given-names></name><name><surname>Reinkensmeyer</surname><given-names>DJ</given-names></name><article-title>A pneumatic robot for re-training arm movement after stroke: Rationale and mechanical design</article-title><source>Proceedings of the 2005 IEEE 9th International Conference on Rehabilitation Robotics</source><year>2005</year><fpage>500</fpage><lpage>504</lpage><comment>full_text</comment></mixed-citation></ref><ref id="B13"><mixed-citation publication-type="other"><name><surname>Stienen</surname><given-names>A</given-names></name><name><surname>Hekman</surname><given-names>E</given-names></name><name><surname>Van der Helm</surname><given-names>FCT</given-names></name><name><surname>Prange</surname><given-names>GB</given-names></name><name><surname>Jannink</surname><given-names>MJA</given-names></name><name><surname>Aalsma</surname><given-names>AMM</given-names></name><name><surname>Van der Kooij</surname><given-names>H</given-names></name><article-title>Dampace: dynamic force-coordination trainer for the upper extremities</article-title><source>Proceedings of the 10th IEEE International Conference on Rehabilitation Robotics 12-15 June 2007; Netherlands</source><fpage>820</fpage><lpage>826</lpage></mixed-citation></ref><ref id="B14"><mixed-citation publication-type="other"><name><surname>Pratt</surname><given-names>G</given-names></name><name><surname>Williamson</surname><given-names>MM</given-names></name><article-title>Series Elastic Actuators</article-title><source>Proceedings of the IEEE International Conference on Intellegent Robots and Systems 1</source><year>1995</year><fpage>399</fpage><lpage>406</lpage></mixed-citation></ref><ref id="B15"><mixed-citation publication-type="book"><name><surname>Robinson</surname><given-names>W</given-names></name><article-title>Design and analysis of series elasticity in closed-loop actuator force control</article-title><source>Ph.D, dissertation</source><year>2000</year><publisher-name>Department of Mechanical Engineering, MIT</publisher-name></mixed-citation></ref><ref id="B16"><mixed-citation publication-type="book"><name><surname>Wyeth</surname><given-names>G</given-names></name><article-title>Control issues for velocity sourced series elastic actuators</article-title><source>Proceedings of the 2006 Australian Conference on Robotics and Automation</source><year>2006</year><publisher-name>ACRA</publisher-name></mixed-citation></ref><ref id="B17"><mixed-citation publication-type="journal"><name><surname>Oblak</surname><given-names>J</given-names></name><name><surname>Cikajlo</surname><given-names>I</given-names></name><name><surname>Matjačć</surname><given-names>Z</given-names></name><article-title>Universal haptic drive: A robot for arm and wrist rehabilitation</article-title><source>IEEE Trans Neural Syst Rehabil Eng</source><year>2010</year><volume>18</volume><issue>3</issue><fpage>293</fpage><lpage>302</lpage><pub-id pub-id-type="doi">10.1109/TNSRE.2009.2034162</pub-id><pub-id pub-id-type="pmid">19846386</pub-id></mixed-citation></ref><ref id="B18"><mixed-citation publication-type="book"><name><surname>Oblak</surname><given-names>J</given-names></name><name><surname>Perry</surname><given-names>JC</given-names></name><name><surname>Jung</surname><given-names>JH</given-names></name><name><surname>Cikajlo</surname><given-names>I</given-names></name><name><surname>Keller</surname><given-names>T</given-names></name><name><surname>Matjaćić</surname><given-names>Z</given-names></name><article-title>A variable structure pantograph mechanism for comprehensive upper extremity haptic movement training</article-title><source>Engineering in Medicine and Biology Society (EMBC), 2010 Annual International Conference of the IEEE, Aug. 31 2010-Sept. 4 2010</source><publisher-name>Buenos Aires</publisher-name><fpage>5859</fpage><lpage>5862</lpage></mixed-citation></ref><ref id="B19"><mixed-citation publication-type="journal"><name><surname>Massie</surname><given-names>TH</given-names></name><name><surname>Salisbury</surname><given-names>JK</given-names></name><article-title>The phantom haptic interface: a device for probing virtual objects</article-title><source>Proceedings of the ASME Winter Annual Meeting: Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systems</source><year>1994</year><volume>55</volume><issue>1</issue><fpage>295</fpage><lpage>302</lpage></mixed-citation></ref><ref id="B20"><mixed-citation publication-type="book"><name><surname>Williamson</surname><given-names>MM</given-names></name><article-title>Series Elastic Actuators</article-title><source>Master's dissertation</source><year>1995</year><publisher-name>Department of Electrical Engineering and Computer Science, MIT</publisher-name></mixed-citation></ref><ref id="B21"><mixed-citation publication-type="journal"><name><surname>Vallery</surname><given-names>H</given-names></name><name><surname>Veneman</surname><given-names>J</given-names></name><name><surname>Van Asseldonk</surname><given-names>E</given-names></name><name><surname>Ekkelenkamp</surname><given-names>R</given-names></name><name><surname>Buss</surname><given-names>M</given-names></name><name><surname>Van der Kooij</surname><given-names>H</given-names></name><article-title>Compliant Actuation of Rehabilitation Robots - Benefits and Limitations of Series Elastic Actuators</article-title><source>IEEE Robotics and Automation Magazine</source><year>2008</year><volume>15</volume><issue>3</issue><fpage>60</fpage><lpage>69</lpage><pub-id pub-id-type="doi">10.1109/MRA.2008.927689</pub-id></mixed-citation></ref><ref id="B22"><mixed-citation publication-type="journal"><name><surname>Colgate</surname><given-names>JE</given-names></name><name><surname>Schenkel</surname><given-names>G</given-names></name><article-title>Passivity of a class of sampled-data systems: application to haptic interfaces</article-title><source>Journal of Robotic Systems</source><year>1997</year><volume>14</volume><issue>1</issue><fpage>37</fpage><lpage>47</lpage><pub-id pub-id-type="doi">10.1002/(SICI)1097-4563(199701)14:1<37::AID-ROB4>3.0.CO;2-V</pub-id></mixed-citation></ref><ref id="B23"><mixed-citation publication-type="other"><name><surname>Colgate</surname><given-names>JE</given-names></name><name><surname>Grafing</surname><given-names>PE</given-names></name><name><surname>Stanley</surname><given-names>MC</given-names></name><name><surname>Schenkel</surname><given-names>G</given-names></name><article-title>Implementation of stiff virtual walls in force-reflecting interfaces</article-title><source>Proc IEEE Virtual Reality Annu Int Symp</source><year>1993</year><fpage>202</fpage><lpage>207</lpage><comment>full_text</comment></mixed-citation></ref><ref id="B24"><mixed-citation publication-type="other"><name><surname>Colgate</surname><given-names>E</given-names></name><name><surname>Hogan</surname><given-names>N</given-names></name><article-title>An analysis of Contact Instability in Terms of Passive Equivalents</article-title><source>IEEE Intl Conf on Robotics and Automation</source><year>1989</year><fpage>404</fpage><lpage>409</lpage></mixed-citation></ref><ref id="B25"><mixed-citation publication-type="journal"><name><surname>Hogan</surname><given-names>N</given-names></name><article-title>On the Stability of Manipulators Performing Contact Tasks</article-title><source>IEEE Journal of Robotics and Automation</source><year>1988</year><volume>4</volume><issue>6</issue><pub-id pub-id-type="doi">10.1109/56.9305</pub-id></mixed-citation></ref><ref id="B26"><mixed-citation publication-type="journal"><name><surname>Gil</surname><given-names>JJ</given-names></name><name><surname>Avello</surname><given-names>A</given-names></name><name><surname>Rubio</surname><given-names>Á</given-names></name><name><surname>Flórez</surname><given-names>J</given-names></name><article-title>Stability Analysis of a 1 DOF Haptic Interface Using the Routh-Hurwitz Criterion</article-title><source>IEEE Trans Con-trol Syst Techno</source><year>2004</year><volume>12</volume><issue>4</issue><fpage>583</fpage><lpage>588</lpage><pub-id pub-id-type="doi">10.1109/TCST.2004.825134</pub-id></mixed-citation></ref><ref id="B27"><mixed-citation publication-type="journal"><name><surname>Diolaiti</surname><given-names>N</given-names></name><name><surname>Niemeyer</surname><given-names>G</given-names></name><name><surname>Barbagli</surname><given-names>F</given-names></name><name><surname>Salisbury</surname><given-names>JK</given-names></name><article-title>Stability of Haptic Rendering: Discretization, Quantization, Time-Delay and Coulomb Effects</article-title><source>IEEE Trans Robot</source><year>2006</year><volume>22</volume><issue>2</issue><fpage>256</fpage><lpage>268</lpage><pub-id pub-id-type="doi">10.1109/TRO.2005.862487</pub-id></mixed-citation></ref></ref-list></back></pmc><affiliations><list><country><li>Slovénie</li></country></list><tree><country name="Slovénie"><noRegion><name sortKey="Oblak, Jakob" sort="Oblak, Jakob" uniqKey="Oblak J" first="Jakob" last="Oblak">Jakob Oblak</name></noRegion><name sortKey="Matja I, Zlatko" sort="Matja I, Zlatko" uniqKey="Matja I Z" first="Zlatko" last="Matja I">Zlatko Matja I</name></country></tree></affiliations></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Pmc/Checkpoint
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001B88 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd -nk 001B88 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Ticri/CIDE
|area= HapticV1
|flux= Pmc
|étape= Checkpoint
|type= RBID
|clé= PMC:3063208
|texte= Design of a series visco-elastic actuator for multi-purpose rehabilitation haptic device
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/RBID.i -Sk "pubmed:21251299" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd \
| NlmPubMed2Wicri -a HapticV1
| This area was generated with Dilib version V0.6.23. |  |