**Towards a Realistic and Self-Contained Biomechanical Model of the Hand**

Joaquín L. Sancho-Bru et al.\* *Universitat Jaume I Spain* 

### **1. Introduction**

210 Theoretical Biomechanics

RANKL product data sheet. Product data sheet alx-522-131: Fc (human):rankl, soluble

http://www.enzolifesciences.com/fileadmin/reports/Fc\_humanRANKL\_Soluble\_

Rice J. C., Cowin S. C., Bowman J. A. On the dependence of the elasticity and strenght of

Robling AG, Castillo AB, and Turner CH. Biomechanical and molecular regulation of bone remodeling. Annual Review of Biomedical Engineering, 8:455–498, 2006. Rubin J, Murphy T, Nanes MS et al. Mechanical strain inhibits expression of osteoclast

Rubin J, Murphy TC, Zhu L et al. Mechanical strain differentially regulates endothelial

via erk1/2 mapk. Journal of Biological Chemistry, 278:34018–34025, 2004. Rubin J, Rubin C, and Jacobs RC. Molecular pathways mediating mechanical signaling in

Ruimerman R, Hilbers PAJ, van Rietbergen B et al. A theoretical framework for strain-

Saquib N., von Muhlen D., Garland CF et al. Serum 25-hydroxyvitamin D, parathyroid

Sikavitsas VI., Temenoff JS. and Mikos AG. Biomaterials and bone mechanotransduction,

Tanaka S., Takahashi N., Udagawa N et al. Macrophage colony-stimulating factor is

Turner CH, Anne V, and Pidaparti RMV. A uniform strain criterion for trabecular bone

Valle F, Sandal M, and Samori B. The interplay between chemistry and mechanics in the

van't Hof RJ., Ralston SH. Cytokine-induced nitric oxide inhibits bone resorption by

Wimalawansa SJ., Chapa MT., Yallampalli C et al. Prevention of corticosteroid-induced bone loss with nitric oxide donor nitroglycerin in male rats. Bone 21 (3), 275-280, 1997 Wimalawansa SJ. Nitroglycerin therapy is as efficacious as standard estrogen replacement

Journal of Bone and Mineral Research 12 (11), 1797-1804, 1997. van't Hof RJ., Ralston SH. Nitric oxide and bone. Immunology 103, 255-261, 2001.

cancellous bone on apparent density. Journal of Biomechanics 21 (2), 155-168. 1988.

differentiation factor by murine stromal cells. Am J Physiol-Cell Ph, 278:C1126–32,

nitric-oxide synthase and receptor activator of nuclear kappa b ligand expression

related trabecular bone maintenance and adaptation. J Biomech, (38):931–941, 2005.

hormone, and bone mineral density in men: the rancho bernardo study.

indispensable for both proliferation and differentiation of osteoclast progenitors.

adaptation: Do continuum-level strain gradients drive adaptation? J Biomech,

transduction of a mechanical signal into a biochemical function. Physics of Life

inducing apoptosis of osteoclast progenitors and supressing osteoclast activity.

therapy (premarin) in prevention of oophorectomy-induced bone loss: A human pilot clinical study. Journal of Bone and Mineral Research 15 (11), 2240-2244, 2000. Wimalawansa SJ. Rationale for using nitric oxide donor therapy for prevention of bone loss and treatment of osteoporosis in humans. Ann NY Acad Sci 1117, 283-297, 2007. Zhang, N. et al. Relationship between age-related serum concentrations of tgf-β1 and tgf- β 2

and those of osteoprotegerin and leptin in native chinese women. Clinical Chimica

(mouse) (rec.), enzo life sciences. 2008 [online]

Osteoporosis International 17, 1734-1741, 2006

Journal of Clinical Investigation 91, 257-263, 1993

Acta. Doi:10.1016/j.cca.2009.01.021, 2009.

mouse\_rec\_rep\_Xcq55b.pdf

bone. Gene, 367:1–16, 2006.

Biomaterials 22: 2581-2593, 2001

30(6):555–63, 1997.

Reviews, 4:157–188, 2007.

2000.

Most of human mechanical interactions with the surrounding world are performed by the hands. They allow us to perform very different tasks; from exerting high forces (e.g. using a hammer) to executing very precise movements (e.g. cutting with a surgical tool). This versatility is possible because of a very complex constitution: a great number of bones connected through different joints, a complicated musculature and a dense nervous system. This complexity is already evident from the kinematics point of view, with more than 20 degrees of freedom (DOF) controlled by muscles, tendons and ligaments.

Mathematical representations are used in order to perform qualitative or quantitative analyses on this complex reality. These representations are known as biomechanical models of the hand. In biomechanics, their use allows studying problems that cannot be analysed directly on humans or that have an experimental cost that is too high; e.g., the study of new alternatives for restoring hand pathologies. Biomechanical models are a description of the hand as a mechanical device: the different elements of the hand are defined in terms of rigid bodies, joints and actuators, and the mechanical laws are applied. As they are simplified mathematical models of the reality, their use and validity depends on the simplifications considered.

The first biomechanical models of the hand were developed to explain and clarify the functionality of different anatomical elements. In this regard, we can find many works that studied the function of the intrinsic muscles (Leijnse & Kalker, 1995; Spoor, 1983; Spoor & Landsmeer, 1976; Storace & Wolf, 1979, 1982; Thomas et al., 1968) and many others that tried to give an insight into the movement coordination of the interphalangeal joints (Buchner et al., 1988; Lee & Rim, 1990). Models for studying the causes and effects of different pathologies of the hand also appeared early on, such as the swan neck and boutonnière deformities or the rupture of the triangular ligament or the volar displacement of the extensor tendon (Smith et al., 1964; Storace & Wolf, 1979, 1982). All these models were, though, very limited, twodimensional models allowing only the study of flexion-extension movements, they modelled only one finger, and they included important simplifications. By the year 2000, few threedimensional models had been developed (Biryukova & Yourovskaya, 1994; Casolo & Lorenzi,

<sup>\*</sup> Antonio Pérez-González, Marta C. Mora, Beatriz E. León, Margarita Vergara, José L. Iserte,Pablo J. Rodríguez-Cervantes and Antonio Morales *Universitat Jaume I, Spain* 

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 213

Storace & Wolf, 1979). Others were developed to help in medical planning and surgery for patients; they are usually dynamic models and are used to perform quantitative analyses, such as the study of the tendon excursions in the medical planning of tendon transfers (Giurintano & Hollister, 1991) or to study the nervous stimulation required to restore the grasping ability in muscular dysfunction patients (Esteki & Mansour, 1997). Yet others studied the hand while performing specific tasks with different aims, so as to have approximate values for the articular forces for testing prosthetic designs (Weightman &

Recent models do not differ much from the ones developed before 2000 (Fok & Chou, 2010; Kamper et al., 2006; Kurita et al., 2009; Lee et al., 2008a, 2008b; Qiu et al., 2009; Roloff et al., 2006; Sancho-Bru et al., 2001, 2003a, 2003b, 2008; Valero-Cuevas, 2000; Valero-Cuevas et al., 2000, 2005; Vigouroux et al., 2006, 2008; Wu et al., 2010). All models present a similar configuration. The kinematics are modelled without considering the restraining structure, just the resultant physiological articular movement. The concept of instantaneous centre of rotation has been used to define an axis of rotation in joints with a single predominant DOF. Much effort has also been spent on finding the rotation axes of joints with two DOF (Brand & Hollister, 1992), through the consideration of a virtual link connecting the axes (Giurintano et al., 1995). Thus, all works use fixed axes of rotation; depending on the joint, one or two axes of rotation are considered. This approximation has been found to be good enough for most of the cases, particularly if there is no interest in analysing the role of the

All works in the literature consider the ideal case of a non-friction belt around a pulley to model the tendons on a joint. Therefore, the tensional force on a tendon is the same along its pathway if no split or connection to other tendon exists. Two different approaches have been used to model tendon action on the joints. The first one considers the tendon freely running when crossing the joint between two points attached one to the proximal segment of the joint and the other to the distal segment. This approach is the basis of the first serious attempt to develop a 3D normative model of the hand (An et al., 1979), in which the position of the tendons with respect the bone segments were obtained from the measurement on 10 fresh cadaveric specimens. The second approach comes from the application of the virtual work principle and considers the moment arm created by the tendon as the first derivative of the tendon excursion with respect to the rotated angle about the rotation axis under study (Storace & Wolf, 1979). This second approach is not strictly correct (Casolo & Lorenzi, 1994), as it does not take into account the work due to the deformation of the sheaths and other structures that constrain the tendon's trajectory along its pathway. Although both approaches present advantages and disadvantages, the second one is difficult to implement in 3D modelling, mainly because of the complexity in the tendon excursion calculation at

Most of the works in the literature use Hill's model to account for the muscles' mathematical modelling. This simple model allows the consideration of the three main parameters, i.e., muscle activation level and variation of the maximum deliverable muscle force with muscle

Finally, the dynamic equilibrium equations lead to an indeterminate system of equations, with more unknowns (muscle forces) than available equations. Inequality constraints taking into account the maximal forces that may be delivered by each muscle and that tendons cannot support compressive forces have to be considered as well. The problem is usually solved by minimising some cost function. Different functions have been investigated, most

Amis, 1982). These too are quantitative analyses performed on dynamic models.

articular soft tissue or the articular stresses (Youm et al., 1978).

joints with more than one DOF.

length and muscle contraction velocity.

1994; Chao et al., 1976; Chao & An, 1978; Esteki & Mansour, 1997; Mansour et al., 1994; Valero-Cuevas et al., 1998), and none of them modelled the complete hand.

Since 2000, many three-dimensional biomechanical models can be found in literature, having been developed for very different purposes (Fok & Chou, 2010; Kamper et al., 2006; Kurita et al., 2009; Lee et al., 2008a, 2008b; Qiu et al., 2009; Roloff et al., 2006; Sancho-Bru et al., 2001, 2003a, 2003b, 2008; Valero-Cuevas, 2000; Valero-Cuevas et al., 2000, 2005; Vigouroux et al., 2006, 2008; Wu et al., 2010): to understand the role of the different anatomical elements, to understand the causes and effects of pathologies, to simulate neuromuscular abnormalities, to plan rehabilitation, to simulate tendon transfer and joint replacement surgeries, to analyse the energetics of human movement and athletic performance, to design prosthetics and biomedical implants, to design functional electric stimulation controllers, to name a few. These models, however, do not differ much from the ones developed before 2000, and many limitations are still evident. For example, contact forces and zones need to be measured experimentally and input to the model.

In contrast, much research has been carried out on animation techniques over the past years, mainly for use in developing computer games. Lately, these advances have been cleverly used by some ergonomics researchers to develop improved graphical and kinematics hand models for evaluating the use of products (Endo et al., 2007, 2009; Goussous, 2007; Kawaguchi, 2009), with good results.

On the other hand, robot hand grasps have been extensively studied for years. Although until 2000 little attention was paid to human hand grasping, this too has become a hot topic in robotics. The experience in modelling the robot hand grasps has been used to reach a better understanding of human grasping (Miller & Allen, 2004; Peña-Pitarch, 2007). The hand is considered as the end-effector for humans. These models, however, are not appropriate for studying many of the above-mentioned objectives, as their interest is different. The focus in robotics research is on planning the grasp and finding an optimum grasp, and quality grasp measures that have been developed for robots are used. These models do not include muscles and tendons in the formulation.

The latest developments in ergonomic hand models and human hand grasp models can be used to improve the existing biomechanical models of the hand and extend their functionality. A promising research area lies ahead with scientist, aiming to obtain a more comprehensive model of the hand, integrating knowledge and developments from the fields of biomechanics, ergonomics, robotics, and computer animation.

In this chapter, a review of the literature regarding biomechanical models of the hand, ergonomics hand models and human hand grasp models is presented. The three approaches are used to draw out the rules for developing an improved biomechanical model, able to tackle any of the above-mentioned objectives in a virtual environment, without external experimental data.

### **2. Literature review**

### **2.1 Biomechanical models of the hand**

Over the years, biomechanical models of the hand have been developed for different purposes. Some of them tried to study the functionality of different anatomical elements with the aim of gaining a deeper understanding of the causes and effects of many hand pathologies. These are usually very simplified (mostly two-dimensional) kinematic models (sometimes dynamic) that are used to perform qualitative analyses (Leijnse et al., 1992;

1994; Chao et al., 1976; Chao & An, 1978; Esteki & Mansour, 1997; Mansour et al., 1994; Valero-

Since 2000, many three-dimensional biomechanical models can be found in literature, having been developed for very different purposes (Fok & Chou, 2010; Kamper et al., 2006; Kurita et al., 2009; Lee et al., 2008a, 2008b; Qiu et al., 2009; Roloff et al., 2006; Sancho-Bru et al., 2001, 2003a, 2003b, 2008; Valero-Cuevas, 2000; Valero-Cuevas et al., 2000, 2005; Vigouroux et al., 2006, 2008; Wu et al., 2010): to understand the role of the different anatomical elements, to understand the causes and effects of pathologies, to simulate neuromuscular abnormalities, to plan rehabilitation, to simulate tendon transfer and joint replacement surgeries, to analyse the energetics of human movement and athletic performance, to design prosthetics and biomedical implants, to design functional electric stimulation controllers, to name a few. These models, however, do not differ much from the ones developed before 2000, and many limitations are still evident. For example, contact

In contrast, much research has been carried out on animation techniques over the past years, mainly for use in developing computer games. Lately, these advances have been cleverly used by some ergonomics researchers to develop improved graphical and kinematics hand models for evaluating the use of products (Endo et al., 2007, 2009; Goussous, 2007;

On the other hand, robot hand grasps have been extensively studied for years. Although until 2000 little attention was paid to human hand grasping, this too has become a hot topic in robotics. The experience in modelling the robot hand grasps has been used to reach a better understanding of human grasping (Miller & Allen, 2004; Peña-Pitarch, 2007). The hand is considered as the end-effector for humans. These models, however, are not appropriate for studying many of the above-mentioned objectives, as their interest is different. The focus in robotics research is on planning the grasp and finding an optimum grasp, and quality grasp measures that have been developed for robots are used. These

The latest developments in ergonomic hand models and human hand grasp models can be used to improve the existing biomechanical models of the hand and extend their functionality. A promising research area lies ahead with scientist, aiming to obtain a more comprehensive model of the hand, integrating knowledge and developments from the fields

In this chapter, a review of the literature regarding biomechanical models of the hand, ergonomics hand models and human hand grasp models is presented. The three approaches are used to draw out the rules for developing an improved biomechanical model, able to tackle any of the above-mentioned objectives in a virtual environment, without external

Over the years, biomechanical models of the hand have been developed for different purposes. Some of them tried to study the functionality of different anatomical elements with the aim of gaining a deeper understanding of the causes and effects of many hand pathologies. These are usually very simplified (mostly two-dimensional) kinematic models (sometimes dynamic) that are used to perform qualitative analyses (Leijnse et al., 1992;

Cuevas et al., 1998), and none of them modelled the complete hand.

forces and zones need to be measured experimentally and input to the model.

models do not include muscles and tendons in the formulation.

of biomechanics, ergonomics, robotics, and computer animation.

Kawaguchi, 2009), with good results.

experimental data.

**2. Literature review** 

**2.1 Biomechanical models of the hand** 

Storace & Wolf, 1979). Others were developed to help in medical planning and surgery for patients; they are usually dynamic models and are used to perform quantitative analyses, such as the study of the tendon excursions in the medical planning of tendon transfers (Giurintano & Hollister, 1991) or to study the nervous stimulation required to restore the grasping ability in muscular dysfunction patients (Esteki & Mansour, 1997). Yet others studied the hand while performing specific tasks with different aims, so as to have approximate values for the articular forces for testing prosthetic designs (Weightman & Amis, 1982). These too are quantitative analyses performed on dynamic models.

Recent models do not differ much from the ones developed before 2000 (Fok & Chou, 2010; Kamper et al., 2006; Kurita et al., 2009; Lee et al., 2008a, 2008b; Qiu et al., 2009; Roloff et al., 2006; Sancho-Bru et al., 2001, 2003a, 2003b, 2008; Valero-Cuevas, 2000; Valero-Cuevas et al., 2000, 2005; Vigouroux et al., 2006, 2008; Wu et al., 2010). All models present a similar configuration. The kinematics are modelled without considering the restraining structure, just the resultant physiological articular movement. The concept of instantaneous centre of rotation has been used to define an axis of rotation in joints with a single predominant DOF. Much effort has also been spent on finding the rotation axes of joints with two DOF (Brand & Hollister, 1992), through the consideration of a virtual link connecting the axes (Giurintano et al., 1995). Thus, all works use fixed axes of rotation; depending on the joint, one or two axes of rotation are considered. This approximation has been found to be good enough for most of the cases, particularly if there is no interest in analysing the role of the articular soft tissue or the articular stresses (Youm et al., 1978).

All works in the literature consider the ideal case of a non-friction belt around a pulley to model the tendons on a joint. Therefore, the tensional force on a tendon is the same along its pathway if no split or connection to other tendon exists. Two different approaches have been used to model tendon action on the joints. The first one considers the tendon freely running when crossing the joint between two points attached one to the proximal segment of the joint and the other to the distal segment. This approach is the basis of the first serious attempt to develop a 3D normative model of the hand (An et al., 1979), in which the position of the tendons with respect the bone segments were obtained from the measurement on 10 fresh cadaveric specimens. The second approach comes from the application of the virtual work principle and considers the moment arm created by the tendon as the first derivative of the tendon excursion with respect to the rotated angle about the rotation axis under study (Storace & Wolf, 1979). This second approach is not strictly correct (Casolo & Lorenzi, 1994), as it does not take into account the work due to the deformation of the sheaths and other structures that constrain the tendon's trajectory along its pathway. Although both approaches present advantages and disadvantages, the second one is difficult to implement in 3D modelling, mainly because of the complexity in the tendon excursion calculation at joints with more than one DOF.

Most of the works in the literature use Hill's model to account for the muscles' mathematical modelling. This simple model allows the consideration of the three main parameters, i.e., muscle activation level and variation of the maximum deliverable muscle force with muscle length and muscle contraction velocity.

Finally, the dynamic equilibrium equations lead to an indeterminate system of equations, with more unknowns (muscle forces) than available equations. Inequality constraints taking into account the maximal forces that may be delivered by each muscle and that tendons cannot support compressive forces have to be considered as well. The problem is usually solved by minimising some cost function. Different functions have been investigated, most

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 215

direct kinematics, but do not allow for example inverse kinematics for the joints of the hand, even when it is incorporated for the other joints of the body. In recent developments some attempts to improve the hand model incorporated into some programs have been done

Early models of the hand (Davidoff & Freivalds, 1993) were actually kinematic models that simulated roughly the external geometry of the hand and its movements. The geometry of the hand has been modelled mainly by jointed cylinders (Fig. 1) and cones (Armstrong, 2009; Sancho-Bru et al., 2003a, 2003b). However, if the geometry of the hand model is not very accurate, the algorithms for inverse kinematics are not precise enough. Recently, some efforts have been made in accurately modelling the surface of real hands to be incorporated into 3D hand models. Rhee et al. (2006) presented an automated method to make a specific human hand model from an image of the palm of the hand. Different algorithms were used in the process: principal creases are extracted, joint locations are estimated from them and the skin geometry of a generic hand model deformed based on hand contours. Rogers et al (2008) made a scalable 3-D geometric model of the hand based on 66 landmarks of the palm surface from 100 subjects in four functional postures. The purpose was to analyse the deformation of the palm surface during the grasp of an object. Recent models incorporate the surface of the hand as a mesh object with more or less realism, obtained from the location of a number of landmarks of the hand or from digital 3D-scanning of the hand (Endo et al., 2007; Peña-Pitarch, 2007, van Nierop et al., 2008). The mesh is linked to a skeleton whose movement controls the deformation of the mesh with different types of

Fig. 1. Different views of the geometric model used in Sancho-Bru et al. (2003a) simulating a

hand gripping two cylinders of different diameters.

(Peña-Pitarch, 2007, Yang et al., 2007).

algorithms.

of them without any physiological basis. The most often used criterion is the minimisation of the sum of the squared muscle stresses, which has been related to the maximisation of fatigue resistance (Crowninshield & Brand, 1981).

All the effort in biomechanics has been focused on appropriately modelling the different hand components (kinematics, muscles, tendons, etc.). Little effort has been spent on the formulation of the grasping problem when using a biomechanical model. In this sense, many limitations persist. Current models do not allow the estimation of the contact information required to use biomechanical models for simulating the grasping of objects. Forces and zones of contact still need to be measured experimentally and input to the model.

### **2.2 Hand models in ergonomics**

Ergonomics, according to the International Ergonomics Association, is 'the scientific discipline concerned with the understanding of interactions among humans and other elements of a system, and the profession that applies theory, principles, data and methods to design in order to optimise human well-being and overall system performance'. Hand models in Ergonomics are used to simulate postures adopted while grasping objects with different purposes. One of the main goals of physical ergonomics is the study of the size and shape of objects according to the anthropometry of the different people that have to interact with them. Thus, the main feature of a model for Ergonomics is that it has to allow representing different populations and percentiles. People having hands of different sizes and proportions will adopt different postures in grasping the same object for the same functions. For example, pressing a button of a phone with the thumb while holding it with the same hand can be easily achieved for a specific hand size while keeping the grasp. However, other people with different size of hand will need to change the grasping posture to achieve pressing the button. This is a typical problem of reach that needs to be solved in ergonomic assessment.

In recent years, virtual humans have been incorporated into the design process for ergonomic assessment of different types of products, mainly in the aerospace and automotive industry but also in others like product design, tasks simulation, personnel training or simulation of other worker environments (Colombo & Cugini, 2005; Yang et al., 2007). Several commercial software programs such as Jack, RAMSIS, HumanCAD, Safework and SantosHuman are available and other studies have been conducted on digital human models such as SAMMIE (Case et al. 1990) or the Boeing Human Modeling System for the same purposes. A virtual human in these packages is defined as a kinematic chain composed of a number of rigid links connected by joints. These joints have the DOF and allowable motion limits corresponding to the anatomical joint of the human being. Direct and inverse kinematics is incorporated into the models so they can replicate human body movements and also evaluate forces acting in joints. Moreover, different population and percentiles may be selected for the size of the model, usually from known anthropometric databases. With these capabilities the problems of reach and clearance usual in ergonomics may be solved easily. Other useful capabilities of these models are the simulation of the sense of sight with virtual cameras located in the eyes or the possibility to change any particular data of the model, like dimensions of limbs or motion limits of some joints, in order to simulate a particular person or disability. However, the majority of these models focuses on the whole body and does not pay attention to the accuracy of the hand model. Most of them just incorporate a list of hand postures (grasping or others) to be chosen, i.e.

of them without any physiological basis. The most often used criterion is the minimisation of the sum of the squared muscle stresses, which has been related to the maximisation of

All the effort in biomechanics has been focused on appropriately modelling the different hand components (kinematics, muscles, tendons, etc.). Little effort has been spent on the formulation of the grasping problem when using a biomechanical model. In this sense, many limitations persist. Current models do not allow the estimation of the contact information required to use biomechanical models for simulating the grasping of objects. Forces and zones of contact still need to be measured experimentally and input to the

Ergonomics, according to the International Ergonomics Association, is 'the scientific discipline concerned with the understanding of interactions among humans and other elements of a system, and the profession that applies theory, principles, data and methods to design in order to optimise human well-being and overall system performance'. Hand models in Ergonomics are used to simulate postures adopted while grasping objects with different purposes. One of the main goals of physical ergonomics is the study of the size and shape of objects according to the anthropometry of the different people that have to interact with them. Thus, the main feature of a model for Ergonomics is that it has to allow representing different populations and percentiles. People having hands of different sizes and proportions will adopt different postures in grasping the same object for the same functions. For example, pressing a button of a phone with the thumb while holding it with the same hand can be easily achieved for a specific hand size while keeping the grasp. However, other people with different size of hand will need to change the grasping posture to achieve pressing the button. This is a typical problem of reach that needs to be solved in

In recent years, virtual humans have been incorporated into the design process for ergonomic assessment of different types of products, mainly in the aerospace and automotive industry but also in others like product design, tasks simulation, personnel training or simulation of other worker environments (Colombo & Cugini, 2005; Yang et al., 2007). Several commercial software programs such as Jack, RAMSIS, HumanCAD, Safework and SantosHuman are available and other studies have been conducted on digital human models such as SAMMIE (Case et al. 1990) or the Boeing Human Modeling System for the same purposes. A virtual human in these packages is defined as a kinematic chain composed of a number of rigid links connected by joints. These joints have the DOF and allowable motion limits corresponding to the anatomical joint of the human being. Direct and inverse kinematics is incorporated into the models so they can replicate human body movements and also evaluate forces acting in joints. Moreover, different population and percentiles may be selected for the size of the model, usually from known anthropometric databases. With these capabilities the problems of reach and clearance usual in ergonomics may be solved easily. Other useful capabilities of these models are the simulation of the sense of sight with virtual cameras located in the eyes or the possibility to change any particular data of the model, like dimensions of limbs or motion limits of some joints, in order to simulate a particular person or disability. However, the majority of these models focuses on the whole body and does not pay attention to the accuracy of the hand model. Most of them just incorporate a list of hand postures (grasping or others) to be chosen, i.e.

fatigue resistance (Crowninshield & Brand, 1981).

**2.2 Hand models in ergonomics** 

ergonomic assessment.

model.

direct kinematics, but do not allow for example inverse kinematics for the joints of the hand, even when it is incorporated for the other joints of the body. In recent developments some attempts to improve the hand model incorporated into some programs have been done (Peña-Pitarch, 2007, Yang et al., 2007).

Early models of the hand (Davidoff & Freivalds, 1993) were actually kinematic models that simulated roughly the external geometry of the hand and its movements. The geometry of the hand has been modelled mainly by jointed cylinders (Fig. 1) and cones (Armstrong, 2009; Sancho-Bru et al., 2003a, 2003b). However, if the geometry of the hand model is not very accurate, the algorithms for inverse kinematics are not precise enough. Recently, some efforts have been made in accurately modelling the surface of real hands to be incorporated into 3D hand models. Rhee et al. (2006) presented an automated method to make a specific human hand model from an image of the palm of the hand. Different algorithms were used in the process: principal creases are extracted, joint locations are estimated from them and the skin geometry of a generic hand model deformed based on hand contours. Rogers et al (2008) made a scalable 3-D geometric model of the hand based on 66 landmarks of the palm surface from 100 subjects in four functional postures. The purpose was to analyse the deformation of the palm surface during the grasp of an object. Recent models incorporate the surface of the hand as a mesh object with more or less realism, obtained from the location of a number of landmarks of the hand or from digital 3D-scanning of the hand (Endo et al., 2007; Peña-Pitarch, 2007, van Nierop et al., 2008). The mesh is linked to a skeleton whose movement controls the deformation of the mesh with different types of algorithms.

Fig. 1. Different views of the geometric model used in Sancho-Bru et al. (2003a) simulating a hand gripping two cylinders of different diameters.

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 217

grasp synthesis is the problem of finding a suitable set of contacts given an object and some

In the following sections, a detailed description of the contact models and the most common

A contact can be defined as a joint between the finger and the object. Their shape, stiffness and frictional characteristics define the nature of this joint (Mason, 2001). The force applied by a finger at a contact point generates a *wrench* on the object with force and torque components. The contact model maps the wrench at some reference point of the object, usually the centre of mass. Salisbury (1982) proposed a taxonomy of eight contact models. Among these, the more common contact models used in robotic grasping (Fig. 2) are the point contacts with and without friction and the soft-finger contacts (Roa Garzón, 2009). Point contact models, also named rigid-body contact models, assume rigid-body models for the hand and the grasped object while the soft-finger contact models, also called compliant or regularised models, assume that the hand is a deformable element grasping a rigid body (Kao et al., 2008). The former models assume the collision to be an instantaneous and discontinuous phenomenon (discrete event) and the equations of motion are derived by balancing the system's momenta before and after the impact. In contrast, compliant models

A *point contact without friction* can only transmit forces along the normal to the object surface at the contact point. No deformations are allowed at the points of contact between the two bodies and, instead, contact forces arise from the constraint of incompressibility and impenetrability between the rigid bodies. These models do not represent the real contact situations that appear in robotic manufacturing operations (Cutkosky, 1989; Lin et al., 2000) and, when used, the machine accuracy is negatively affected. Moreover, they are not capable of predicting the individual contact forces of a multiple-contact fixture (Bicchi, 1994; Harada

A *point contact with friction* can also transmit forces in the tangential directions to the surface at the contact point. If Coulomb's friction model is used, all the forces that lie within the friction cone with an angle *atan(µ)* can be applied, where *µ* is the friction coefficient of the contacting materials. Here, contact forces arise from two sources: the rigid-body model assumption for both the hand and the object, and the frictional forces. The use of this contact model in the manipulation planning problem has led to some interesting conclusions. There may be multiple solutions to a particular problem (ambiguity) or there may be no solutions

Finally, the soft contact model is used to model the contact between a soft finger and a rigid object allowing the finger to apply an additional torsional moment with respect to the normal at the contact point (Ciocarlie et al., 2005, 2007; Howe et al., 1988; Howe & Cutkosky, 1996; Kao & Cutkosky, 1992; Kao & Yang, 2004). A typical contact between a soft finger and a contact surface can be modelled by the Hertzian contact model (Hertz, 1882; Johnson, 1985). However, robotic fingertips are made of nonlinear elastic materials. For that reason, the Hertzian contact model does not accurately represent this contact. In Xydas & Kao (1999) and Xydas et al. (2000) a power-law theory is presented for modelling nonlinear elastic contacts present in robotic fingers. It subsumes the Hertzian contact theory. More realistic, and complicated, models have been developed in the last few years that better represent the

constrains on the allowable contacts.

**2.3.1 Grasp contact models** 

et al., 2000).

(inconsistency) (Erdmann, 1994).

approaches for grasp analysis and synthesis is presented.

describe the normal and tangential compliance relations over time.

Other important aspect of hand models for ergonomics is associated with the study of musculoskeletal disorders. Early epidemiological studies (Mital and Kilbom, 1992) showed that the use of hand tools with an improper design for the worker or the task could lead to a high risk of developing cumulative hand trauma disorders (CHTD). The factors influencing the development of CHTD have been reported in different works (Keyserling, 2000; Kong et al., 2006; Muggleton et al., 1999; Schoenmarklin et al., 1994; Spielholz et al., 2001) and different methods have been used in these studies: epidemiological studies, physiological measurements (electromyography activity, pressure in tissues, posture of hand and wrist, tactile sensitivity), biomechanical models of hand and wrist structures and psychophysical assessments. These studies report that CHTD are associated with repetitive tasks, high forces, extreme or awkward postures of hand and wrist, velocity and acceleration of wrist motions and exposure time, among others. Different theories of injury development have been proposed (Kumar, 2001). All of them assume that CHTD and other musculoskeletal disorders are of biomechanical nature.

Therefore, biomechanical hand models able to predict movements, postures and internal forces of hand and wrist structures can be used to assess the risk of developing CHTD. Tendon excursions or maximum gripping strength have been used as index in different works to assess gripping posture for health (Armstrong et al., 2009, Sancho-Bru et al., 2003b).

None of the reported biomechanical models of the hand for ergonomics accounts for all the above-mentioned requirements, although some attempts have been made. Armstrong et al (2009) have developed a scalable kinematic model of the hand with simple geometry (cones and cylinders). The model includes a posture prediction algorithm for fingers that reproduces in a high percentage the observed postures and is able to compute tendon excursions and wrist movements. The model is used to assess how much space is required for hands in an assembly task and to calculate the risk of CHTD from tendon forces and hand strength. Other group of researchers (Endo et al 2007, 2009; Kawaguchi et al., 2009) have developed a scalable digital hand model with an accurate shape of the hand that includes a semiautomatic grasp planning function with robotics indexes of quality (see next section). The model incorporates a 'comfort database' obtained from experimental measurements to assess comfort of postures and is used in the assessment of physical interaction with electronic appliances.

### **2.3 Grasping in robotics**

For many years the robotics community has been studying the autonomous handling of objects by robots. A robot should be able to locate the object and then grasp it, and possibly transport it to a specified destination. The purpose of a grasp is to constrain the potential movements of the object in the event of external disturbances. For a specific robotic hand, different grasp types are planned and analysed in order to decide which one to execute.

A grasp is commonly defined as a set of contacts on the surface of the object. A contact model should be defined to determine the forces or torques that the robot manipulator must exert on the contact areas. Most of the work in robotics assumes point contacts, and larger areas of contact are usually discretised to follow this assumption (Bicchi & Kumar, 2000).

Two main problems can be distinguished in robotic grasping: analysis and synthesis (Mason, 2001). Grasp analysis consists on finding whether the grasp is stable using common closure properties, given an object and a set of contacts. Then, a quality measure can be evaluated in order to enable the robot to select the best grasp to execute. On the other hand, grasp synthesis is the problem of finding a suitable set of contacts given an object and some constrains on the allowable contacts.

In the following sections, a detailed description of the contact models and the most common approaches for grasp analysis and synthesis is presented.

### **2.3.1 Grasp contact models**

216 Theoretical Biomechanics

Other important aspect of hand models for ergonomics is associated with the study of musculoskeletal disorders. Early epidemiological studies (Mital and Kilbom, 1992) showed that the use of hand tools with an improper design for the worker or the task could lead to a high risk of developing cumulative hand trauma disorders (CHTD). The factors influencing the development of CHTD have been reported in different works (Keyserling, 2000; Kong et al., 2006; Muggleton et al., 1999; Schoenmarklin et al., 1994; Spielholz et al., 2001) and different methods have been used in these studies: epidemiological studies, physiological measurements (electromyography activity, pressure in tissues, posture of hand and wrist, tactile sensitivity), biomechanical models of hand and wrist structures and psychophysical assessments. These studies report that CHTD are associated with repetitive tasks, high forces, extreme or awkward postures of hand and wrist, velocity and acceleration of wrist motions and exposure time, among others. Different theories of injury development have been proposed (Kumar, 2001). All of them assume that CHTD and other musculoskeletal

Therefore, biomechanical hand models able to predict movements, postures and internal forces of hand and wrist structures can be used to assess the risk of developing CHTD. Tendon excursions or maximum gripping strength have been used as index in different works to assess gripping posture for health (Armstrong et al., 2009, Sancho-Bru et al.,

None of the reported biomechanical models of the hand for ergonomics accounts for all the above-mentioned requirements, although some attempts have been made. Armstrong et al (2009) have developed a scalable kinematic model of the hand with simple geometry (cones and cylinders). The model includes a posture prediction algorithm for fingers that reproduces in a high percentage the observed postures and is able to compute tendon excursions and wrist movements. The model is used to assess how much space is required for hands in an assembly task and to calculate the risk of CHTD from tendon forces and hand strength. Other group of researchers (Endo et al 2007, 2009; Kawaguchi et al., 2009) have developed a scalable digital hand model with an accurate shape of the hand that includes a semiautomatic grasp planning function with robotics indexes of quality (see next section). The model incorporates a 'comfort database' obtained from experimental measurements to assess comfort of postures and is used in the assessment of physical

For many years the robotics community has been studying the autonomous handling of objects by robots. A robot should be able to locate the object and then grasp it, and possibly transport it to a specified destination. The purpose of a grasp is to constrain the potential movements of the object in the event of external disturbances. For a specific robotic hand, different grasp types are planned and analysed in order to decide which one to execute. A grasp is commonly defined as a set of contacts on the surface of the object. A contact model should be defined to determine the forces or torques that the robot manipulator must exert on the contact areas. Most of the work in robotics assumes point contacts, and larger areas of contact are usually discretised to follow this assumption (Bicchi & Kumar, 2000). Two main problems can be distinguished in robotic grasping: analysis and synthesis (Mason, 2001). Grasp analysis consists on finding whether the grasp is stable using common closure properties, given an object and a set of contacts. Then, a quality measure can be evaluated in order to enable the robot to select the best grasp to execute. On the other hand,

disorders are of biomechanical nature.

interaction with electronic appliances.

**2.3 Grasping in robotics** 

2003b).

A contact can be defined as a joint between the finger and the object. Their shape, stiffness and frictional characteristics define the nature of this joint (Mason, 2001). The force applied by a finger at a contact point generates a *wrench* on the object with force and torque components. The contact model maps the wrench at some reference point of the object, usually the centre of mass. Salisbury (1982) proposed a taxonomy of eight contact models. Among these, the more common contact models used in robotic grasping (Fig. 2) are the point contacts with and without friction and the soft-finger contacts (Roa Garzón, 2009). Point contact models, also named rigid-body contact models, assume rigid-body models for the hand and the grasped object while the soft-finger contact models, also called compliant or regularised models, assume that the hand is a deformable element grasping a rigid body (Kao et al., 2008). The former models assume the collision to be an instantaneous and discontinuous phenomenon (discrete event) and the equations of motion are derived by balancing the system's momenta before and after the impact. In contrast, compliant models describe the normal and tangential compliance relations over time.

A *point contact without friction* can only transmit forces along the normal to the object surface at the contact point. No deformations are allowed at the points of contact between the two bodies and, instead, contact forces arise from the constraint of incompressibility and impenetrability between the rigid bodies. These models do not represent the real contact situations that appear in robotic manufacturing operations (Cutkosky, 1989; Lin et al., 2000) and, when used, the machine accuracy is negatively affected. Moreover, they are not capable of predicting the individual contact forces of a multiple-contact fixture (Bicchi, 1994; Harada et al., 2000).

A *point contact with friction* can also transmit forces in the tangential directions to the surface at the contact point. If Coulomb's friction model is used, all the forces that lie within the friction cone with an angle *atan(µ)* can be applied, where *µ* is the friction coefficient of the contacting materials. Here, contact forces arise from two sources: the rigid-body model assumption for both the hand and the object, and the frictional forces. The use of this contact model in the manipulation planning problem has led to some interesting conclusions. There may be multiple solutions to a particular problem (ambiguity) or there may be no solutions (inconsistency) (Erdmann, 1994).

Finally, the soft contact model is used to model the contact between a soft finger and a rigid object allowing the finger to apply an additional torsional moment with respect to the normal at the contact point (Ciocarlie et al., 2005, 2007; Howe et al., 1988; Howe & Cutkosky, 1996; Kao & Cutkosky, 1992; Kao & Yang, 2004). A typical contact between a soft finger and a contact surface can be modelled by the Hertzian contact model (Hertz, 1882; Johnson, 1985). However, robotic fingertips are made of nonlinear elastic materials. For that reason, the Hertzian contact model does not accurately represent this contact. In Xydas & Kao (1999) and Xydas et al. (2000) a power-law theory is presented for modelling nonlinear elastic contacts present in robotic fingers. It subsumes the Hertzian contact theory. More realistic, and complicated, models have been developed in the last few years that better represent the

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 219

by Ferrari & Canny (1992) being the most widely-used. They proposed to calculate the radius of the largest ball inscribed in the convex hull centred in the origin. Force-closure

Many approaches have been proposed to measure the quality of a grasp. Some of the measures focus on evaluating the ability to resist external disturbances, others on evaluating the dexterity. These measures can be classified into two groups depending on whether they consider the location of the contact points on the object or the configuration of the endeffector. There are also some that are a combination of these two approaches (see Roa

Measures from the first group take into account the geometric properties of the objects, their materials and closure properties to evaluate the grasp. For example, Li & Sastry (1998) proposed to calculate the smallest singular value of the grasp matrix, which indicates how far the grasp configuration is from losing the capability of withstanding external wrenches. Others have proposed to favour the grasps whose contact points are distributed in a uniform way on the object surface, which improves their stability (Mirtich & Canny, 1994; Park & Starr, 1992). This can be done by measuring either the angles or the area of the polygon whose vertices are the contact points. The centroid of this polygon is also used to calculate its distance to the object centre of mass (Ding et al., 2001; Ponce et al., 1997). The smaller this distance the better the grasp can resist the effect of external forces. In addition, some other measures take into account the uncertainty in the position of the fingers; therefore instead of contact points they calculate contact regions in which force closure grasps are assured (Nguyen, 1988; Roa Garzón, 2009). The quality of the grasp is measured

The previous approaches do not consider any limitation on the finger forces, so that in some cases the fingers have to apply very large forces to resist small perturbations. Other measures do consider limitations on the magnitudes of the finger forces. They can limit the force on each finger or the sum of forces applied by all fingers. Ferrari & Canny (1992) used the largest ball not only to evaluate grasp closure but also to measure the grasp quality. This is a geometric representation of the smallest perturbation wrench that breaks the grasp, independently of its direction. It has been widely used by the robotics community (Borst et al., 2003; Miller & Allen, 1999; Roa Garzón, 2009). The volume of the ball is also considered as a quality measure with the advantage that it remains constant independently of the used

When a task is specified to be performed after the object is grasped, the quality of the grasp can be measured with its ability to counteract the expected disturbances during the task execution. The set of all wrenches that are expected to be applied on the object defines the *task wrench space (TWS)* and can be approximated as an ellipsoid (Li & Sastry, 1988) or as a convex polytope (Haschke et al., 2005; Zhu at al., 2001). The problem with these approaches is that modelling the TWS can be quite complicated (Borst et al., 2004). Pollard (2004) introduced the concept of an *object wrench space (OWS)*, which is the set of wrenches generated by applying a distribution of disturbance forces on the surface of the object. Borst et al. (2004) proposed the use of the largest factor by which the OWS can fit the GWS as the

On the other hand, there are measures that consider the configuration of the end-effector, requiring the hand-object Jacobian for their calculation (Roa Garzón, 2009). An example of

grasps are the ones where the sphere's radius is larger than zero.

**Measures of grasp performance** 

Garzón (2009) for a thorough review).

by the size of these regions.

measure of the grasp quality.

torque reference system (Miller & Allen, 1999).

contact mechanics for soft fingers (Ciocarlie et al., 2005, 2007; Gonthier, 2007). However, it is the hard finger contacts with friction that are used more often in robotics.

Fig. 2. Contact models commonly used in robotics: a) Point contact without friction; b) Point contact with friction; c) Soft-finger contact

### **2.3.2 Grasp analysis**

After establishing the contact model, the set of contacts defining each grasp can be analysed in order to test its ability to resist disturbances and its dexterity properties. As it is presented afterwards, a grasp can resist disturbances in any direction if it fulfils one of the two closure conditions*.* However, there is usually more than one grasp that fulfils them. That is why many metrics and approaches have been proposed to evaluate the dexterity of the selected grasps and determine which one is the best to be executed.

### **Disturbance resistance**

The first test for evaluating a grasp consists of determining its ability to constrain the motions of the manipulated object and to apply arbitrary contact forces on the object without violating friction constraints at the contacts (Bicchi, 1995). Two commonly used properties have been proposed to ensure this condition: force and form closure. A grasp is in *force-closure* if the fingers can apply, through the set of contacts, arbitrary wrenches on the object, which means that any motion of the object is resisted by the contact forces (Nguyen, 1988). On the other hand, a grasp is in *form-closure* if the location of the contact points on the object ensures its immobility (Bicchi, 1995).

Form closure is a stronger condition than force closure and it is mostly used when executing power grasps (Siciliano & Khatib, 2008). Force closure is possible with fewer contacts, making it suitable for executing precision grasps, but it requires the ability to control internal forces.

In order to verify the form or force closure property of a grasp, many tests have been proposed (see Liu et al. (2004b) and Roa Garzón (2009) for a review). Most of them dene conditions to be satisfied by the grasp wrenches in the wrench space. A *grasp wrench space (GWS)* is the space of wrenches that can be applied to the object at each contact point. The boundary of the wrench space can be calculated as a convex hull. Force-closure then can be determined verifying if the origin of the wrench space lies inside this convex hull (Mishra et al., 1987). Several tests have been proposed to verify this condition, with the one developed by Ferrari & Canny (1992) being the most widely-used. They proposed to calculate the radius of the largest ball inscribed in the convex hull centred in the origin. Force-closure grasps are the ones where the sphere's radius is larger than zero.

### **Measures of grasp performance**

218 Theoretical Biomechanics

contact mechanics for soft fingers (Ciocarlie et al., 2005, 2007; Gonthier, 2007). However, it is

Fig. 2. Contact models commonly used in robotics: a) Point contact without friction; b) Point

After establishing the contact model, the set of contacts defining each grasp can be analysed in order to test its ability to resist disturbances and its dexterity properties. As it is presented afterwards, a grasp can resist disturbances in any direction if it fulfils one of the two closure conditions*.* However, there is usually more than one grasp that fulfils them. That is why many metrics and approaches have been proposed to evaluate the dexterity of the selected

The first test for evaluating a grasp consists of determining its ability to constrain the motions of the manipulated object and to apply arbitrary contact forces on the object without violating friction constraints at the contacts (Bicchi, 1995). Two commonly used properties have been proposed to ensure this condition: force and form closure. A grasp is in *force-closure* if the fingers can apply, through the set of contacts, arbitrary wrenches on the object, which means that any motion of the object is resisted by the contact forces (Nguyen, 1988). On the other hand, a grasp is in *form-closure* if the location of the contact points on the

Form closure is a stronger condition than force closure and it is mostly used when executing power grasps (Siciliano & Khatib, 2008). Force closure is possible with fewer contacts, making it suitable for executing precision grasps, but it requires the ability to control

In order to verify the form or force closure property of a grasp, many tests have been proposed (see Liu et al. (2004b) and Roa Garzón (2009) for a review). Most of them dene conditions to be satisfied by the grasp wrenches in the wrench space. A *grasp wrench space (GWS)* is the space of wrenches that can be applied to the object at each contact point. The boundary of the wrench space can be calculated as a convex hull. Force-closure then can be determined verifying if the origin of the wrench space lies inside this convex hull (Mishra et al., 1987). Several tests have been proposed to verify this condition, with the one developed

the hard finger contacts with friction that are used more often in robotics.

contact with friction; c) Soft-finger contact

object ensures its immobility (Bicchi, 1995).

grasps and determine which one is the best to be executed.

**2.3.2 Grasp analysis** 

**Disturbance resistance** 

internal forces.

Many approaches have been proposed to measure the quality of a grasp. Some of the measures focus on evaluating the ability to resist external disturbances, others on evaluating the dexterity. These measures can be classified into two groups depending on whether they consider the location of the contact points on the object or the configuration of the endeffector. There are also some that are a combination of these two approaches (see Roa Garzón (2009) for a thorough review).

Measures from the first group take into account the geometric properties of the objects, their materials and closure properties to evaluate the grasp. For example, Li & Sastry (1998) proposed to calculate the smallest singular value of the grasp matrix, which indicates how far the grasp configuration is from losing the capability of withstanding external wrenches. Others have proposed to favour the grasps whose contact points are distributed in a uniform way on the object surface, which improves their stability (Mirtich & Canny, 1994; Park & Starr, 1992). This can be done by measuring either the angles or the area of the polygon whose vertices are the contact points. The centroid of this polygon is also used to calculate its distance to the object centre of mass (Ding et al., 2001; Ponce et al., 1997). The smaller this distance the better the grasp can resist the effect of external forces. In addition, some other measures take into account the uncertainty in the position of the fingers; therefore instead of contact points they calculate contact regions in which force closure grasps are assured (Nguyen, 1988; Roa Garzón, 2009). The quality of the grasp is measured by the size of these regions.

The previous approaches do not consider any limitation on the finger forces, so that in some cases the fingers have to apply very large forces to resist small perturbations. Other measures do consider limitations on the magnitudes of the finger forces. They can limit the force on each finger or the sum of forces applied by all fingers. Ferrari & Canny (1992) used the largest ball not only to evaluate grasp closure but also to measure the grasp quality. This is a geometric representation of the smallest perturbation wrench that breaks the grasp, independently of its direction. It has been widely used by the robotics community (Borst et al., 2003; Miller & Allen, 1999; Roa Garzón, 2009). The volume of the ball is also considered as a quality measure with the advantage that it remains constant independently of the used torque reference system (Miller & Allen, 1999).

When a task is specified to be performed after the object is grasped, the quality of the grasp can be measured with its ability to counteract the expected disturbances during the task execution. The set of all wrenches that are expected to be applied on the object defines the *task wrench space (TWS)* and can be approximated as an ellipsoid (Li & Sastry, 1988) or as a convex polytope (Haschke et al., 2005; Zhu at al., 2001). The problem with these approaches is that modelling the TWS can be quite complicated (Borst et al., 2004). Pollard (2004) introduced the concept of an *object wrench space (OWS)*, which is the set of wrenches generated by applying a distribution of disturbance forces on the surface of the object. Borst et al. (2004) proposed the use of the largest factor by which the OWS can fit the GWS as the measure of the grasp quality.

On the other hand, there are measures that consider the configuration of the end-effector, requiring the hand-object Jacobian for their calculation (Roa Garzón, 2009). An example of

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 221

a) b)

c) d)

Arm System (courtesy of DLR Robotics and Mechatronics Center)

fast) during the object manipulation.

Fig. 3. Anthropomorphic robot hands: a) Barrett Hand (courtesy of the UJI Robotics Intelligent Lab); b) ARMAR III Hand (courtesy of the Institute for Anthropomatics at KIT); c) Shadow Hand C5 (courtesy of Shadow Robot Company); d) Anthropomorphic DLR Hand

different robotic grasping postures and could be adapted to human grasping.

neuromuscular control. The features that we require in order to create a model are:

models lack realism for assessing the use of handheld products from an ergonomics point of view. Hand models in ergonomics have reached a high level of realism but do not allow for mechanical analyses. On the other hand, biomechanical models are not self-contained, as they need contact information to be input to the model. Current models do not allow predicting grasping postures nor evaluating contact forces and zones, much less predicting the movements while grasp planning. Quality grasp measures in robotics allow comparing

A detailed proposal for modelling the different components of the hand is provided below: joints-kinematics, muscles, ligaments and passive tissues, skin, contact with objects and

• The model has to simulate the complete hand in order to allow the study of any grasp. • The model has to be scalable to allow the simulation of different population groups. • The model has to simulate and show the grasping of an object in a realistic way.

• The model has to estimate the muscular patterns required to perform a movement while counteracting the system of external forces that define the object manipulation.

Furthermore, the model has to estimate the articular forces at the hand joints. • The model has to be dynamic in order to allow the study of any grasping task (slow or

this group is a measure that favours a grasp that, given certain velocities in the finger joints, produces the largest velocities on the grasped object, calculated with the volume of the manipulability ellipsoid (Yoshikawa, 1985). There is another measure that penalises the joints of the hand being in their maximum limits, calculating the deviation of the joint angles from their centres (Liegeois, 1977). Additionally, there are other measures in this group that also consider the task, giving higher quality indexes to the grasps which ensure the maximum transformation ratio along the direction wrenches more likely to be applied on the object when executing it (Chiu, 1988).

### **2.3.3 Grasp synthesis**

Given an object, grasp synthesis algorithms should provide a suitable set of contacts on the object surface and determine an appropriate hand configuration. Usually they take the geometry of the object as an input to select optimal force-closure contact locations. These contacts are the starting point for grasp analysis and dexterous manipulation methods.

Some approaches give only information about the finger contact locations on the object without considering the hand constraints. They can result in stable grasps that are not reachable in practice by the robot hand. Moreover, even if they are reachable, it is difficult to position the fingers precisely on the contact points because there will be always unavoidable errors locating the end-effector (Morales et al., 2006).

Alternative approaches, called knowledge-based approaches, have considered the configuration of the hand by generating the grasp with a predefined set of hand postures. The idea of hand preshapes started with studies of the human prehension capabilities (Napier, 1956) that introduced the distinction between power and precision grasps. Following this work, Cutkosky (1989) created a taxonomy in which details of the task and the object geometry are taken into account. Since then, several papers have adopted this approach for grasping (Morales et al., 2006; Stansfield, 1991; Wren, 1995). Miller et al. (2003) used a simulator called GraspIt! to test the set of hand preshapes on a 3D model of the object. Using a simulator has many advantages, including the ability to plan grasps in complex environments involving obstacles and also to check the reachability constraints of the robot arm. More recently OpenRAVE, a planning architecture that has a more exible design, has been proposed to automate this process (Diankov, 2010).

Despite many years of research and all the advances we have reviewed, the robotics community is still not able to build a manipulator with similar capabilities to the human hand. The robot hands constructed until now are only simplifications (Fig. 3), given the complexities not only at the sensor and actuator level, but also at the control level. They vary from the easiest to control, such as 2-jaw grippers, to more anthropomorphic hands like the Salisbury Hand, the Utah-MIT Hand, the Barrett Hand, the ARMAR III Hand or the DLR Hand II (see Biagiotti et al. (2002) and Parada et al. (2008) for a review).

### **3. Hand biomechanical model proposal**

In this section, the current knowledge on biomechanical, ergonomics and robotics hand models is used to draw out the rules for developing a realistic and self-contained biomechanical model of the hand.

Based on the literature review, current hand biomechanical models allow estimating the muscular patterns required to perform a movement while counteracting a system of external forces. But their use for studying object grasping is limited. On the one hand, biomechanical

this group is a measure that favours a grasp that, given certain velocities in the finger joints, produces the largest velocities on the grasped object, calculated with the volume of the manipulability ellipsoid (Yoshikawa, 1985). There is another measure that penalises the joints of the hand being in their maximum limits, calculating the deviation of the joint angles from their centres (Liegeois, 1977). Additionally, there are other measures in this group that also consider the task, giving higher quality indexes to the grasps which ensure the maximum transformation ratio along the direction wrenches more likely to be applied on

Given an object, grasp synthesis algorithms should provide a suitable set of contacts on the object surface and determine an appropriate hand configuration. Usually they take the geometry of the object as an input to select optimal force-closure contact locations. These contacts are the starting point for grasp analysis and dexterous manipulation methods. Some approaches give only information about the finger contact locations on the object without considering the hand constraints. They can result in stable grasps that are not reachable in practice by the robot hand. Moreover, even if they are reachable, it is difficult to position the fingers precisely on the contact points because there will be always unavoidable

Alternative approaches, called knowledge-based approaches, have considered the configuration of the hand by generating the grasp with a predefined set of hand postures. The idea of hand preshapes started with studies of the human prehension capabilities (Napier, 1956) that introduced the distinction between power and precision grasps. Following this work, Cutkosky (1989) created a taxonomy in which details of the task and the object geometry are taken into account. Since then, several papers have adopted this approach for grasping (Morales et al., 2006; Stansfield, 1991; Wren, 1995). Miller et al. (2003) used a simulator called GraspIt! to test the set of hand preshapes on a 3D model of the object. Using a simulator has many advantages, including the ability to plan grasps in complex environments involving obstacles and also to check the reachability constraints of the robot arm. More recently OpenRAVE, a planning architecture that has a more exible

Despite many years of research and all the advances we have reviewed, the robotics community is still not able to build a manipulator with similar capabilities to the human hand. The robot hands constructed until now are only simplifications (Fig. 3), given the complexities not only at the sensor and actuator level, but also at the control level. They vary from the easiest to control, such as 2-jaw grippers, to more anthropomorphic hands like the Salisbury Hand, the Utah-MIT Hand, the Barrett Hand, the ARMAR III Hand or the DLR

In this section, the current knowledge on biomechanical, ergonomics and robotics hand models is used to draw out the rules for developing a realistic and self-contained

Based on the literature review, current hand biomechanical models allow estimating the muscular patterns required to perform a movement while counteracting a system of external forces. But their use for studying object grasping is limited. On the one hand, biomechanical

the object when executing it (Chiu, 1988).

errors locating the end-effector (Morales et al., 2006).

design, has been proposed to automate this process (Diankov, 2010).

Hand II (see Biagiotti et al. (2002) and Parada et al. (2008) for a review).

**3. Hand biomechanical model proposal** 

biomechanical model of the hand.

**2.3.3 Grasp synthesis** 

Fig. 3. Anthropomorphic robot hands: a) Barrett Hand (courtesy of the UJI Robotics Intelligent Lab); b) ARMAR III Hand (courtesy of the Institute for Anthropomatics at KIT); c) Shadow Hand C5 (courtesy of Shadow Robot Company); d) Anthropomorphic DLR Hand Arm System (courtesy of DLR Robotics and Mechatronics Center)

models lack realism for assessing the use of handheld products from an ergonomics point of view. Hand models in ergonomics have reached a high level of realism but do not allow for mechanical analyses. On the other hand, biomechanical models are not self-contained, as they need contact information to be input to the model. Current models do not allow predicting grasping postures nor evaluating contact forces and zones, much less predicting the movements while grasp planning. Quality grasp measures in robotics allow comparing different robotic grasping postures and could be adapted to human grasping.

A detailed proposal for modelling the different components of the hand is provided below: joints-kinematics, muscles, ligaments and passive tissues, skin, contact with objects and neuromuscular control. The features that we require in order to create a model are:


Towards a Realistic and Self-Contained Biomechanical Model of the Hand 223

defining a flexion/extension axis of rotation connecting the carpus to each metacarpal. The orientation of the axes is defined oblique in order to appropriately simulate the arching of the palm (Kapandji, 1998). Due to the important role that the shape of the palm plays in grasping, this model is considered more suitable for grasping simulation than others in the

The data for the location and orientation of the rotation axes comes from An et al. (1979), Buchholz et al. (1992) and Hollister et al. (1995). Axes data and link lengths are fully scaled with respect to the hand length and hand breadth (Sancho-Bru, 2000). Limits for the joints

In order to study the forward and inverse kinematics of the hand, the Denavit- Hartenberg method from the robotics field (Denavit and Hartenberg, 1955) was adapted to define the

Muscles and tendons control the movement of the skeletal chains. Muscles have been considered using a simple Hill three-component model (Hill, 1938) that takes into account the muscle activation level (α) and the force-length and force-velocity relationships, as well as the different index of architecture of muscles. The model considers a contractile element (CE), which is the basic component that generates force, a parallel elastic element (PEE), which is responsible for the passive force generated by the muscle when it is stretched, and a series elastic element (SEE), the muscle tendon unit, which has been considered to be

SEE

= ⋅ *F PCSA S* max max , (1)

α

PEE

conditions is proportional to its physiological cross-sectional area (PCSA):

inextensible, so that the force the muscle exerts (*F*) can be written as:

The force a muscle can exert depends on the actual muscle length and contraction velocity. It is widely accepted (An et al., 1991) that the maximum force a muscle can exert in optimal

where *Smax* is the maximum stress the muscle can bear, which has been considered the same

The strain of tendons is insignificant for the magnitude of forces developed by the muscles (Goldstein et al., 1987). Under this consideration, the SEE has been considered to

CE

Fig. 5. Hill's three-component model for the muscles

for each muscle (An et al., 1991).

have been obtained from Tubiana (1981) and Tubiana et al. (1996).

literature.

position of any segment point.

**3.2 Musculo-tendon action** 

inextensible (Fig. 5).


The model proposed in this section has been developed in a scalable way, choosing two very well known anthropometric parameters of the hand that are easy to measure and representative of the hand size. The parameters are the hand length (HL) and hand breadth (HB) and are shown in Fig. 4.

Fig. 4. Parameters used to scale the model: HL (hand length) and HB (hand breadth)

### **3.1 Kinematics**

In order to achieve realistic grasping postures, care has to be taken when selecting the appropriate DOF among the different hand bones. The DOF have to allow the hand model to reach the hand posture for any grasping task. In this sense, it is important that the model considers not only the thumb and finger movements but also the palm arching.

The hand has been considered as five skeletal open chains of rigid bodies connected to the carpus through different joints which characterise the kinematic behaviour of the chains.

Distal and proximal interphalangeal (DIP and PIP) joints of the fingers as well as the interphalangeal (IP) joint of the thumb are trochlear joints, capable only of flexion/extension movements (Brand & Hollister, 1992). These joints are modelled as one DOF joints by means of defining a rotation axis connecting the adjacent phalanxes (hinge joint).

Thumb and fingers metacarpophalangeal (MCP) joints are condylar joints, capable of flexion/extension and abduction/adduction movements (Brand & Hollister, 1992). The thumb carpometacarpal (CMC) joint is a saddle joint, capable also of flexion/extension and abduction/adduction movements (Brand & Hollister, 1992). All these joints are modelled as two DOF joints by defining two axes of rotation connecting the adjacent segments. The axes are neither intersecting nor orthogonal (Brand & Hollister, 1992), so that a virtual link is used to connect both axes (Giurintano et al., 1995).

Finally, the hand model allows the arching of the palm by modelling the CMC joints of the little and ring fingers. These joints are arthrodial joints, with a very limited range of movement (Kapandji, 1998). They have been modelled as one DOF joints by means of defining a flexion/extension axis of rotation connecting the carpus to each metacarpal. The orientation of the axes is defined oblique in order to appropriately simulate the arching of the palm (Kapandji, 1998). Due to the important role that the shape of the palm plays in grasping, this model is considered more suitable for grasping simulation than others in the literature.

The data for the location and orientation of the rotation axes comes from An et al. (1979), Buchholz et al. (1992) and Hollister et al. (1995). Axes data and link lengths are fully scaled with respect to the hand length and hand breadth (Sancho-Bru, 2000). Limits for the joints have been obtained from Tubiana (1981) and Tubiana et al. (1996).

In order to study the forward and inverse kinematics of the hand, the Denavit- Hartenberg method from the robotics field (Denavit and Hartenberg, 1955) was adapted to define the position of any segment point.

### **3.2 Musculo-tendon action**

222 Theoretical Biomechanics

• The model has to predict feasible grasping postures for a given object and provide the

• The model has to incorporate quality grasping measures for evaluating the grasp. The model proposed in this section has been developed in a scalable way, choosing two very well known anthropometric parameters of the hand that are easy to measure and representative of the hand size. The parameters are the hand length (HL) and hand breadth

Fig. 4. Parameters used to scale the model: HL (hand length) and HB (hand breadth)

considers not only the thumb and finger movements but also the palm arching.

of defining a rotation axis connecting the adjacent phalanxes (hinge joint).

used to connect both axes (Giurintano et al., 1995).

In order to achieve realistic grasping postures, care has to be taken when selecting the appropriate DOF among the different hand bones. The DOF have to allow the hand model to reach the hand posture for any grasping task. In this sense, it is important that the model

The hand has been considered as five skeletal open chains of rigid bodies connected to the carpus through different joints which characterise the kinematic behaviour of the chains. Distal and proximal interphalangeal (DIP and PIP) joints of the fingers as well as the interphalangeal (IP) joint of the thumb are trochlear joints, capable only of flexion/extension movements (Brand & Hollister, 1992). These joints are modelled as one DOF joints by means

Thumb and fingers metacarpophalangeal (MCP) joints are condylar joints, capable of flexion/extension and abduction/adduction movements (Brand & Hollister, 1992). The thumb carpometacarpal (CMC) joint is a saddle joint, capable also of flexion/extension and abduction/adduction movements (Brand & Hollister, 1992). All these joints are modelled as two DOF joints by defining two axes of rotation connecting the adjacent segments. The axes are neither intersecting nor orthogonal (Brand & Hollister, 1992), so that a virtual link is

Finally, the hand model allows the arching of the palm by modelling the CMC joints of the little and ring fingers. These joints are arthrodial joints, with a very limited range of movement (Kapandji, 1998). They have been modelled as one DOF joints by means of

contact information required for evaluating the grasp

(HB) and are shown in Fig. 4.

**3.1 Kinematics** 

Muscles and tendons control the movement of the skeletal chains. Muscles have been considered using a simple Hill three-component model (Hill, 1938) that takes into account the muscle activation level (α) and the force-length and force-velocity relationships, as well as the different index of architecture of muscles. The model considers a contractile element (CE), which is the basic component that generates force, a parallel elastic element (PEE), which is responsible for the passive force generated by the muscle when it is stretched, and a series elastic element (SEE), the muscle tendon unit, which has been considered to be inextensible (Fig. 5).

Fig. 5. Hill's three-component model for the muscles

The force a muscle can exert depends on the actual muscle length and contraction velocity. It is widely accepted (An et al., 1991) that the maximum force a muscle can exert in optimal conditions is proportional to its physiological cross-sectional area (PCSA):

$$F\_{\text{max}} = \text{PCSA} \cdot \text{S}\_{\text{max}} \,\prime \tag{1}$$

where *Smax* is the maximum stress the muscle can bear, which has been considered the same for each muscle (An et al., 1991).

The strain of tendons is insignificant for the magnitude of forces developed by the muscles (Goldstein et al., 1987). Under this consideration, the SEE has been considered to inextensible, so that the force the muscle exerts (*F*) can be written as:

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 225

finger muscles have been taken from Lee & Rim (1990); data for the remaining extrinsic muscles have been obtained from Lemay & Crago (1996) and for the remaining intrinsic muscles from Jacobson et al. (1992). The muscle maximal lengthening velocity ( max

> Index Medial Ring Little Thumb 1st FP 2nd FP 3rd FP 4th FP APB 1st FS 2nd FS 3rd FS 4th FS FPB 1st EDC+EI 2nd EDC 3rd EDC EDQ OPP 1st LU 2nd LU 3rd LU 4th LU ADD 1st DI 2nd DI 4th DI 3rd VI 1st DI 1st VI 3rd DI 2nd VI FDQ APL ADQ EPB FPL EPL

Table 1. Muscles modelled on each skeletal chain (acronyms in the nomenclature section)

points defining the tendon actions have been obtained from An et al. (1979).

amount of fibres of each branch (Eyler & Markee, 1954) has been considered.

circular path has been considered.

Most of the muscles do not act directly on the bones, but transmit the force to the tendons, which finally insert into the bones. To model the tendon action crossing the joints, straight lines connecting 2 points have been considered, one fixed with respect to the proximal bone and the other one with respect to the distal bone (Fig. 6a). This approximation has been found to be close enough to the behaviour of all tendons with the exception of extensors (An et al*.*, 1979), for which Landsmeer's model I has been considered (Fig. 6b). The data for the

Fig. 6. Models for the tendons crossing the joints: a) Straight lines; b) Landsmeer's model I The extensor hood mechanisms of the fingers are modelled as a tendon net. The net allows for the connection and division of the tendon paths. The insertions and connection points considered for the tendon nets on each skeletal chain are shown in Fig. 7. Appropriate force balances have been considered in the connecting points of this deformable tendon net. Second DI, fourth DI and ADQ tendons do present a double insertion into the proximal phalanxes and into the extensor aponeuroses. A force distribution proportional to the

The muscle force-length and force-velocity relationships presented above require the calculation of the lengthening of the muscles from *lo* as a function of time. Having considered the tendons inextensible, the muscle lengthening coincides with the tendon excursion. To calculate the length of the tendon path crossing each joint (*li*), straight lines connecting the points have been considered, except for the extensor tendons, for which a

b)

been taken to be 2.5 s-1 (Kaufman et al*.*, 1991).

a)

ε) has

$$F = F\_{\text{max}} \left( F\_{\text{CE}} + F\_{\text{PEE}} \right) \tag{2}$$

where *FCE* and *FPEE* are the normalised forces delivered by the CE and PEE, respectively. The force exerted by the muscle can be decomposed into an active force and a passive force corresponding to the forces delivered by the CE and PEE, respectively. The force delivered by the CE is related to the muscle architecture and is a function of the muscle length *lCE*, the contraction velocity *vCE*, and the muscle activation level α (from 0 to 1), which is controlled by the central nervous system (Kaufman et al., 1991):

$$F\_{\rm CE} = \mathcal{Q} \cdot F\_l(l\_{\rm CE}) \cdot F\_v(v\_{\rm CE}) \tag{3}$$

where *Fl* and *Fv* are the non-dimensional force-length and force-velocity relationships. A characteristic bell-shaped curve exists between force and length of the muscle. To model this dependence, the expression proposed by Kaufman et al. (1991) has been used:

$$F\_l(\mathcal{E}, \dot{i}\_a) = e^{-\left[\frac{\left(\mathcal{E} \ast \mathbf{1}\right)^{0.8537} \left(1 - \frac{1}{i\_a}\right) - 1.0}{0.35327 \left(1 - i\_a\right)}\right]^2} \quad \text{for } \dot{i}\_a < \mathbf{1} \tag{4}$$

$$F\_l(\mathcal{E}, \mathbf{i}\_a) = e^{-[2.727277 \ln(\mathbf{c} \times \mathbf{i})]^2} \quad \text{for } \mathbf{i}\_a = \mathbf{1} \tag{5}$$

where *ia* is the muscle architecture index, defined as the ratio between the muscle fibre length and the muscle belly length, and ε is the muscle strain due to its lengthening from *lo*, the muscle length for the optimal conditions.

The force a muscle can exert decreases when the contraction velocity of the muscle fibres increases. To model this dependence the expression proposed by Hatze (1981) has been used

$$F\_v(\dot{\eta}) = \frac{0.1433}{0.1074 + e^{-1.409 \sinh(3.2 \cdot \eta + 1.6)}} \tag{6}$$

where η is the normalised contractile element velocity, given by the ratio between the lengthening velocity of the muscle ( ε ), and its maximal value ( max ε).

The force generated by the PEE is a function only of its length. An exponential relationship has been considered in this case (Lee & Rim, 1990; Kaufman et al*.*, 1991), with *b1* and *b2* muscle dependent constants:

2 1 1 *<sup>b</sup> F be b PEE* ⋅ε=⋅ − , (7)

The scalability of the muscular action is achieved by scaling the PCSA of the muscles with respect to the product of hand length and hand breadth parameters (Sancho-Bru et al., 2008) from its value for *HL* = 18.22 cm and *HB* = 8.00 cm.

$$\frac{PCSA(HL, HB)}{PCSA(\overline{HL}, \overline{HB})} = 1 + 0.01333 \cdot \text{(HB} \cdot HL - \overline{HB} \cdot \overline{HL}) \tag{8}$$

The muscles considered on each skeletal chain are listed in Table 1. PCSA data for index finger muscles have been taken from Valero-Cuevas et al*.* (1998); data for the remaining muscles have been obtained from Brand & Hollister (1992). Muscle stress limit (*Smax*) has been obtained from Zajac (1989). Fibre and muscle lengths and the constants *b1*, *b2* for index

where *FCE* and *FPEE* are the normalised forces delivered by the CE and PEE, respectively. The force exerted by the muscle can be decomposed into an active force and a passive force corresponding to the forces delivered by the CE and PEE, respectively. The force delivered by the CE is related to the muscle architecture and is a function of the muscle length *lCE*, the contraction velocity *vCE*, and the muscle activation level α (from 0 to 1), which is controlled

> =⋅ ⋅ α

where *Fl* and *Fv* are the non-dimensional force-length and force-velocity relationships. A characteristic bell-shaped curve exists between force and length of the muscle. To model

this dependence, the expression proposed by Kaufman et al. (1991) has been used:

( )

ε

0.1433 ( ) 0.1074 *Fv <sup>e</sup>*

ε

*Fi e l a ai*

[ ] <sup>2</sup> 2.727277 ln( 1) (, ) for 1 *Fi e l a ai* ε

where *ia* is the muscle architecture index, defined as the ratio between the muscle fibre

The force a muscle can exert decreases when the contraction velocity of the muscle fibres increases. To model this dependence the expression proposed by Hatze (1981) has been used

The force generated by the PEE is a function only of its length. An exponential relationship has been considered in this case (Lee & Rim, 1990; Kaufman et al*.*, 1991), with *b1* and *b2*

> 2 1 1 *<sup>b</sup> F be b PEE* ⋅ε

=+ ⋅ ⋅ − ⋅ (,) 1 0.01333 ( ) (,) *PCSA HL HB HB HL HB HL*

*PCSA HL HB* (8)

The scalability of the muscular action is achieved by scaling the PCSA of the muscles with respect to the product of hand length and hand breadth parameters (Sancho-Bru et al., 2008)

The muscles considered on each skeletal chain are listed in Table 1. PCSA data for index finger muscles have been taken from Valero-Cuevas et al*.* (1998); data for the remaining muscles have been obtained from Brand & Hollister (1992). Muscle stress limit (*Smax*) has been obtained from Zajac (1989). Fibre and muscle lengths and the constants *b1*, *b2* for index

1.409 sinh(3.2 1.6)

is the normalised contractile element velocity, given by the ratio between the

), and its maximal value ( max

η

<sup>2</sup> <sup>1</sup> 0.96343 1 1 1.0 0.35327 (1 ) (, ) for 1 *ia ai*

 <sup>⋅</sup> <sup>−</sup> + − <sup>−</sup> ⋅ −

by the central nervous system (Kaufman et al., 1991):

ε

ε

η

ε

length and the muscle belly length, and

lengthening velocity of the muscle (

muscle dependent constants:

where

η

the muscle length for the optimal conditions.

from its value for *HL* = 18.22 cm and *HB* = 8.00 cm.

= + max ( ) *FF F F CE PEE* , (2)

() ( ) *F Fl F v CE l CE v CE* (3)

= < (4)

is the muscle strain due to its lengthening from *lo*,

− ⋅+ = = (5)

− ⋅ ⋅+ <sup>=</sup> <sup>+</sup> (6)

ε).

=⋅ − , (7)

finger muscles have been taken from Lee & Rim (1990); data for the remaining extrinsic muscles have been obtained from Lemay & Crago (1996) and for the remaining intrinsic muscles from Jacobson et al. (1992). The muscle maximal lengthening velocity ( max ε ) has been taken to be 2.5 s-1 (Kaufman et al*.*, 1991).


Table 1. Muscles modelled on each skeletal chain (acronyms in the nomenclature section)

Most of the muscles do not act directly on the bones, but transmit the force to the tendons, which finally insert into the bones. To model the tendon action crossing the joints, straight lines connecting 2 points have been considered, one fixed with respect to the proximal bone and the other one with respect to the distal bone (Fig. 6a). This approximation has been found to be close enough to the behaviour of all tendons with the exception of extensors (An et al*.*, 1979), for which Landsmeer's model I has been considered (Fig. 6b). The data for the points defining the tendon actions have been obtained from An et al. (1979).

Fig. 6. Models for the tendons crossing the joints: a) Straight lines; b) Landsmeer's model I

The extensor hood mechanisms of the fingers are modelled as a tendon net. The net allows for the connection and division of the tendon paths. The insertions and connection points considered for the tendon nets on each skeletal chain are shown in Fig. 7. Appropriate force balances have been considered in the connecting points of this deformable tendon net. Second DI, fourth DI and ADQ tendons do present a double insertion into the proximal phalanxes and into the extensor aponeuroses. A force distribution proportional to the amount of fibres of each branch (Eyler & Markee, 1954) has been considered.

The muscle force-length and force-velocity relationships presented above require the calculation of the lengthening of the muscles from *lo* as a function of time. Having considered the tendons inextensible, the muscle lengthening coincides with the tendon excursion. To calculate the length of the tendon path crossing each joint (*li*), straight lines connecting the points have been considered, except for the extensor tendons, for which a circular path has been considered.

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 227

The data for the location of the points defining the tendon paths comes from An et al. (1979) and Buchholz et al. (1992), and are also scaled with respect to the hand length and hand

In previous work, we showed the importance of modelling the effect of ligaments for studying free finger movements. In the case of grasping, their consideration is not so relevant. Their effect can be neglected for studying power grasps, but they can play an important role in the case of some precision grasps, particularly those involving fast

In the case of DIP and PIP joints of fingers and thumb, the insertion of the collateral ligaments on the proximal segment of the joint corresponds to the flexion-extension axis (Dubousset, 1981). Therefore, they do not develop any flexion-extension moment over the joint and they do not need to be modelled. In the case of MCP joints, the proximal insertion of the lateral ligament on the metacarpal head remains dorsal to the center of the articular curvature (Fig. 8), so that collateral ligaments are lax in extension, but they become taut in flexion, decreasing significantly the range of lateral movement (Craig, 1992; Dubousset, 1981; Kapandji, 1998). Tension on the radial and ulnar ligaments increases with adduction and abduction of the MCP joint, respectively. Furthermore, the line of action of the ligaments remains dorsal to the flexion- extension axis of the joint (Craig, 1992), developing

an extension moment over the joint, in addition to the abduction-adduction moment.

Both ulnar and radial ligaments over MCP joints have been considered. A unique fibre for each ligament has been considered, joining two points representing the insertions into the bones. One point is fixed with respect to the metacarpal, and the other one with respect to the proximal phalanx. No interaction between bone and ligament has been considered; therefore the ligament path is a straight line between the insertion points. Its non-linear behaviour has been taken into account considering a quadratic relationship between the

( )<sup>2</sup>

where *K* is the characteristic constant of the ligament, *Llig* the length of the fibre representing

*F KL L lig* =⋅ − *lig lig o*, , (10)

force developed by the ligament (*Flig*) and its elongation (Mommersteeg et al., 1996)

Fig. 8. Collateral ligament over MCP joints becomes taut with flexion.

the ligament, and *Llig,o* the unstrained length of the ligament.

Flexion

breadth (Sancho-Bru, 2000).

**3.3 Ligaments** 

movements.

Fig. 7. Sketch of the extensor mechanisms of the fingers and thumb (dorsal view) showing the insertions into the bones () and the connections and splittings considered (): a) little finger; b) ring finger; c) medial finger; d) index finger; e) thumb.

The data for the location of the points defining the tendon paths comes from An et al. (1979) and Buchholz et al. (1992), and are also scaled with respect to the hand length and hand breadth (Sancho-Bru, 2000).

### **3.3 Ligaments**

226 Theoretical Biomechanics

**Distal Proximal**

**EDQ**

**4th LU**

**3rd UI**

**4th DI**

**3rd EDC**

**3rd LU 2ond UI**

**3rd DI**

**2ond EDC**

**2ond LU 2ond DI**

**1st UI**

**1st LU**

**ADD**

**EPL**

**APB**

**1st EDC+EI**

**ADQ**

**DIP PIP MCP**

**DIP PIP MCP**

**DIP PIP MCP**

**DIP PIP MCP**

**IP MCP**

Fig. 7. Sketch of the extensor mechanisms of the fingers and thumb (dorsal view) showing the insertions into the bones () and the connections and splittings considered (): a) little

finger; b) ring finger; c) medial finger; d) index finger; e) thumb.

a)

b)

c)

d)

e)

In previous work, we showed the importance of modelling the effect of ligaments for studying free finger movements. In the case of grasping, their consideration is not so relevant. Their effect can be neglected for studying power grasps, but they can play an important role in the case of some precision grasps, particularly those involving fast movements.

In the case of DIP and PIP joints of fingers and thumb, the insertion of the collateral ligaments on the proximal segment of the joint corresponds to the flexion-extension axis (Dubousset, 1981). Therefore, they do not develop any flexion-extension moment over the joint and they do not need to be modelled. In the case of MCP joints, the proximal insertion of the lateral ligament on the metacarpal head remains dorsal to the center of the articular curvature (Fig. 8), so that collateral ligaments are lax in extension, but they become taut in flexion, decreasing significantly the range of lateral movement (Craig, 1992; Dubousset, 1981; Kapandji, 1998). Tension on the radial and ulnar ligaments increases with adduction and abduction of the MCP joint, respectively. Furthermore, the line of action of the ligaments remains dorsal to the flexion- extension axis of the joint (Craig, 1992), developing an extension moment over the joint, in addition to the abduction-adduction moment.

Fig. 8. Collateral ligament over MCP joints becomes taut with flexion.

Both ulnar and radial ligaments over MCP joints have been considered. A unique fibre for each ligament has been considered, joining two points representing the insertions into the bones. One point is fixed with respect to the metacarpal, and the other one with respect to the proximal phalanx. No interaction between bone and ligament has been considered; therefore the ligament path is a straight line between the insertion points. Its non-linear behaviour has been taken into account considering a quadratic relationship between the force developed by the ligament (*Flig*) and its elongation (Mommersteeg et al., 1996)

$$F\_{\rm lig} = K \cdot \left(L\_{\rm lig} - L\_{\rm lig,o}\right)^2 \tag{10}$$

where *K* is the characteristic constant of the ligament, *Llig* the length of the fibre representing the ligament, and *Llig,o* the unstrained length of the ligament.

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 229

joints at observation-based rates. To select the appropriate rotation rate we propose the use of neural networks, similar to those used in Kyota et al. (2005) and Rezzoug & Gorce (2008). This will require intensive experimental work beforehand to record the postures for grasping objects (of different shapes, sizes, weights, etc.) when performing different tasks (power and precision). The experimental data have to be analysed in order to characterise the human grasp and find the parameters affecting the grasping posture. These parameters will be used as input to the neural network to estimate a tentative clenched posture. The rotation rate is defined by the difference between the angles of the fully opened state and

In order to generate the grasp, a contact model is required. We need to check whether the surface skin model makes contact with the surface of the object model. In reality, the surface of a hand is deformed when making contact with the object. Generally, this deformation has a non-linearly elastic property, and it could be simulated using finite element analyses. But this would need a long execution time. This is unacceptable for our model where a large number of different grasp postures have to be generated and tested within a practical time. Therefore, we propose to consider a simple geometric collision-detection algorithm based on the one used by Endo et al. (2007). The algorithm allows the penetration of the surface skin model and the object model. This penetration is limited by a tolerance that relates to the

The distances between the points on the skin surface and the object are calculated while the joint angles of each joint rotate according to the specific joint rotation algorithm. When the maximum penetration distance between the skin surface points and the object reaches the given tolerance, the contact is achieved and the joint rotation ends. When distal segments of

As we have mentioned previously, the model has to provide the contact information required for evaluating the grasp. If a classical robotics quality measure of the grasp is to be performed, the only data needed are the contact points and associated normal vectors at

When trying to estimate the muscular pattern associated with a grasp, the model needs more contact information. The contact forces between the object and the hand have to be considered in this case. Unlike what happens with robots, real human fingers conform to the grasped object shape. As the contact finger surface is deformable the contact does not occur at just one point but over some finite area that increases as the normal forces increase. Due to this effect, in addition to the normal force and tangential force due to friction, human finger contact may support frictional torsional moments with respect to the normal at the contact point. This clearly shows that the consideration of rigid contact commonly used in robotics is not appropriate for use in studying the human grasp, and a soft contact has to be modelled. Most objects manipulated by human hands are much stiffer than human hands, and it is reasonable to consider the objects to be grasped as rigid bodies, and the hand as a

Different soft contact models have been investigated and proposed in order to better account for this deformation effect in the context of soft finger contacts (Ciocarlie et al., 2005, 2007; Gonthier, 2007). In Ciocarlie et al. (2007) friction constraints are derived based on general expressions for non-planar contacts of elastic bodies, taking into account the local geometry and structure of the objects in contact. The following approximation can be used

> τ*n*):

to express the constraint relating the magnitudes of frictional force (*ft*) and moment (

all four fingers make contact with the object, the simulation terminates.

these points, which are easily obtained from the proposed contact model.

the tentative clenched one.

hand stiffness of each contact region.

deformable body.

The data for the ligament insertion points have been obtained from the geometric model presented in Youm et al. (1978), and the stiffness constant has been estimated to be 750 N/cm2 from Minami et al. (1985).

### **3.4 Skin and contact with objects**

One of the applications of the biomechanical model is its use in assessing the use of handheld products from an ergonomics point of view. To accomplish that goal, the model has to incorporate a realistic model of the skin from the visual point of view. The advances in computer animation have made possible the development of a number of convincing surface skin models.

We propose to use a surface skin model similar to that of Endo et al. (2007) or Goussous (2007). The surface skin model is a 3-dimensional polygonal mesh for the hand surface generated from CT images. The geometry of the skin model is defined at only one opened posture. A surface skin deformation algorithm defines the deformed geometry of the surface skin model when the posture of the kinematic model is changed (Fig. 9). The algorithm assigns each bone a capsule-shaped *envelope*. Vertices of the modified skin within these envelopes move with the bones. Where envelopes overlap, vertex motion is a blend between the envelopes. The influence of each bone for vertices within the intersection of two bones' envelopes is controlled by assigning weight values. The ratio of a vertex's weight values, which always total 1.0, determine the relative extent to which each bone's motion affects the vertex. Furthermore, the model gets scaled when the kinematic model is scaled.

As stated before, the model has to simulate and show the grasping of an object in a realistic way. To satisfy this requirement, it is not enough to have a visually realistic model of the surface skin. The model must also be able to predict feasible grasping postures.

### Fig. 9. Surface skin model

In order to generate grasp postures automatically, we propose to use a grasping algorithm based on that of Choi (2008). This algorithm uses a function to automatically generate a natural grasping motion path of the hand model from a fully opened state to a clenched one. The goal is to find contacts between the surface hand skin and the object surface while rotating the joint angles of the fingers. Care has to be taken to properly choose the rotation rate of the finger joints, as it affects the final posture prediction. Based on the results from Choi (2008), we propose to use a variable rotation algorithm, by describing rotations of all

The data for the ligament insertion points have been obtained from the geometric model presented in Youm et al. (1978), and the stiffness constant has been estimated to be 750

One of the applications of the biomechanical model is its use in assessing the use of handheld products from an ergonomics point of view. To accomplish that goal, the model has to incorporate a realistic model of the skin from the visual point of view. The advances in computer animation have made possible the development of a number of convincing

We propose to use a surface skin model similar to that of Endo et al. (2007) or Goussous (2007). The surface skin model is a 3-dimensional polygonal mesh for the hand surface generated from CT images. The geometry of the skin model is defined at only one opened posture. A surface skin deformation algorithm defines the deformed geometry of the surface skin model when the posture of the kinematic model is changed (Fig. 9). The algorithm assigns each bone a capsule-shaped *envelope*. Vertices of the modified skin within these envelopes move with the bones. Where envelopes overlap, vertex motion is a blend between the envelopes. The influence of each bone for vertices within the intersection of two bones' envelopes is controlled by assigning weight values. The ratio of a vertex's weight values, which always total 1.0, determine the relative extent to which each bone's motion affects the

As stated before, the model has to simulate and show the grasping of an object in a realistic way. To satisfy this requirement, it is not enough to have a visually realistic model of the

In order to generate grasp postures automatically, we propose to use a grasping algorithm based on that of Choi (2008). This algorithm uses a function to automatically generate a natural grasping motion path of the hand model from a fully opened state to a clenched one. The goal is to find contacts between the surface hand skin and the object surface while rotating the joint angles of the fingers. Care has to be taken to properly choose the rotation rate of the finger joints, as it affects the final posture prediction. Based on the results from Choi (2008), we propose to use a variable rotation algorithm, by describing rotations of all

vertex. Furthermore, the model gets scaled when the kinematic model is scaled.

surface skin. The model must also be able to predict feasible grasping postures.

N/cm2 from Minami et al. (1985).

**3.4 Skin and contact with objects** 

surface skin models.

Fig. 9. Surface skin model

joints at observation-based rates. To select the appropriate rotation rate we propose the use of neural networks, similar to those used in Kyota et al. (2005) and Rezzoug & Gorce (2008). This will require intensive experimental work beforehand to record the postures for grasping objects (of different shapes, sizes, weights, etc.) when performing different tasks (power and precision). The experimental data have to be analysed in order to characterise the human grasp and find the parameters affecting the grasping posture. These parameters will be used as input to the neural network to estimate a tentative clenched posture. The rotation rate is defined by the difference between the angles of the fully opened state and the tentative clenched one.

In order to generate the grasp, a contact model is required. We need to check whether the surface skin model makes contact with the surface of the object model. In reality, the surface of a hand is deformed when making contact with the object. Generally, this deformation has a non-linearly elastic property, and it could be simulated using finite element analyses. But this would need a long execution time. This is unacceptable for our model where a large number of different grasp postures have to be generated and tested within a practical time.

Therefore, we propose to consider a simple geometric collision-detection algorithm based on the one used by Endo et al. (2007). The algorithm allows the penetration of the surface skin model and the object model. This penetration is limited by a tolerance that relates to the hand stiffness of each contact region.

The distances between the points on the skin surface and the object are calculated while the joint angles of each joint rotate according to the specific joint rotation algorithm. When the maximum penetration distance between the skin surface points and the object reaches the given tolerance, the contact is achieved and the joint rotation ends. When distal segments of all four fingers make contact with the object, the simulation terminates.

As we have mentioned previously, the model has to provide the contact information required for evaluating the grasp. If a classical robotics quality measure of the grasp is to be performed, the only data needed are the contact points and associated normal vectors at these points, which are easily obtained from the proposed contact model.

When trying to estimate the muscular pattern associated with a grasp, the model needs more contact information. The contact forces between the object and the hand have to be considered in this case. Unlike what happens with robots, real human fingers conform to the grasped object shape. As the contact finger surface is deformable the contact does not occur at just one point but over some finite area that increases as the normal forces increase. Due to this effect, in addition to the normal force and tangential force due to friction, human finger contact may support frictional torsional moments with respect to the normal at the contact point. This clearly shows that the consideration of rigid contact commonly used in robotics is not appropriate for use in studying the human grasp, and a soft contact has to be modelled. Most objects manipulated by human hands are much stiffer than human hands, and it is reasonable to consider the objects to be grasped as rigid bodies, and the hand as a deformable body.

Different soft contact models have been investigated and proposed in order to better account for this deformation effect in the context of soft finger contacts (Ciocarlie et al., 2005, 2007; Gonthier, 2007). In Ciocarlie et al. (2007) friction constraints are derived based on general expressions for non-planar contacts of elastic bodies, taking into account the local geometry and structure of the objects in contact. The following approximation can be used to express the constraint relating the magnitudes of frictional force (*ft*) and moment (τ*n*):

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 231

 <sup>=</sup>

with *n* between 2.0 and 4.0, and where *Fi* represents the force exerted by muscle *i*, and *PCSAi* its physiological cross-sectional area. In this case, *n* = 2 will be used. This function is minimised when subjected to Eq. 13 together with the force balances of the tendon nets. Additional constraints are that tendon forces must be non-negative, and the limits of muscle

A global grasp evaluation can be performed through the use of the proposed model, merging the knowledge from ergonomics, robotics and biomechanics. The classical ergonomics evaluation of grasp posture and reachability is possible, for different percentiles of the population represented by the corresponding anthropometric parameters. Furthermore, CHTD evaluation can be performed by using the predicted postures and

It is advisable to use force closure from robotics as a part of the proposed model; once a grasping posture is estimated by the grasping algorithm, force closure should be assured before consuming time in determining the contact forces. Any of the robotics quality measures could be used for evaluating the grasp. Depending on the task to be performed, it would be better to use a grasp quality measure to evaluate the disturbance resistance or a

But the most relevant contribution to grasp evaluation has to come from biomechanics analysis. Grasp measures related to the muscle and articular forces have to be investigated. Just to provide insight into this sense, and to ensure coherence with our model formulation, we propose the use of Eq. 15 as a quality measure related to fatigue that we can call *fatigue index*: the smaller the fatigue index the better will be the grasp. For power grasps, an alternative measure can be the difference between the maximal force the hand can exert on the grasped object for the posture being analysed and the real contact forces estimated by the contact model; this alternative measure can be seen as a safety margin for the muscle forces, that we can call *muscle safety margin index*. Additional measures can be investigated,

A realistic and self-contained biomechanical model of the hand has been proposed by merging the current knowledge of biomechanics, ergonomics and robotics. The model simulates the complete hand and can easily be scaled to study different percentiles of populations. It has a realistic representation that allows the ergonomic evaluation of products. The model is dynamic and can be used to study the muscular patterns associated with a specific grasp. It allows predicting feasible grasping postures and provides the contact information required for evaluating the grasp. Finally, the model incorporates original quality grasping measures such as the fatigue index and the muscle safety margin index, in addition to the usual robotics and ergonomics metrics and evaluations. All the

grasp quality measure to evaluate the manipulability.

such as the maximal contact pressure, etc.

forces obtained from Eqs. 2 and 3 varying the muscle activation level from 0 to 1

**3.6 Grasp evaluation** 

muscle forces.

**4. Conclusion** 

 

*i i PCSA*

*n*

*<sup>F</sup> OBJ* , (15)

max max *PEE* ( *vl PEE* )⋅+⋅≤≤⋅ *FFFFFFF* . (16)

$$\left\|f\right\|^2 + \frac{\left\|\mathbf{r}\right\|^2}{\left\|e\_n\right\|^2} \le \mu^2 \cdot \mathbf{P}^2 \tag{11}$$

where *P* is the total load applied in the direction of the contact normal, μ is the friction coefficient and *en* is called the eccentricity parameter (height of the ellipsoid described by Eq. 11). Considering a Winkler elastic foundation (Johnson, 1985) of depth *h* and elastic modulus *K*, the eccentricity parameter is given by:

$$
\omega\_n = \frac{8}{15} \sqrt{a \cdot b} \quad , \tag{12}
$$

where *a* and *b* can be calculated from the relative radii of curvature *R'* and *R''* of the objects in contact and the compression δof the elastic layer:

$$a = \sqrt{2 \cdot \delta \cdot R'}; b = \sqrt{2 \cdot \delta \cdot R'}; \delta = \sqrt{\frac{P \cdot h}{K \cdot \pi \cdot \left(R' \cdot R'\right)^{1/2}}}\tag{13}$$

The actual grasping forces for a given posture will be obtained by considering that they have to satisfy the dynamic equilibrium of the grasped object. There is not a unique set of forces that ensures the equilibrium but we have to take into account the biomechanical limitations (maximal muscle forces) and the control performed by the central nervous system. In an effort to minimise the computational cost, we propose to uncouple the computation of the contact forces from the neuromuscular control model. This can be done by considering that the central nervous system is trying to attempt performing the grasp with minimal contact forces, as implemented for robots in the work of Liu et al. (2004a).

### **3.5 Neuromuscular control**

The movement of the skeletal chains, together with the contact forces and the corresponding application points are input to the model. The problem to be solved is the derivation of the muscle activation levels required to produce the given motion under the external loads. It is, therefore, an inverse dynamics problem.

The dynamics equations of the open chain of rigid bodies have been derived using the Lagrange method (García de Jalón & Bayo, 1994). For a system with *m* generalised coordinates *qk*, this equation is expressed as:

$$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}\_k} - \frac{\partial L}{\partial q\_k} = Q\_k^{uc} \quad k = 1, \dots, m \tag{14}$$

where *L* is the Lagrangian function and *nc Qk* are the generalised non-conservative forces. The generalised coordinates have been considered coincident with the system DOF (*m=23*).

Eqs. 14 together with the force balances of the tendon nets lead to an indeterminate problem. For example, in the case of the index finger, there are 12 equations (four corresponding to the DOF considered and eight to force balances in the tendon net) and 18 unknowns (six muscle forces and 12 branch forces of the tendon net). There is not a unique combination of muscular efforts that satisfy the dynamic equilibrium constraints. To solve the problem, a criterion chosen by the central nervous system to determine the muscle action control must be introduced. Our proposal is to maximise the endurance. According to Crowninshield and Brand (1981), this is achieved by minimising the non-linear objective function

$$COB = \sum \left(\frac{F\_i}{PCSA\_i}\right)^{\nu} \tag{15}$$

with *n* between 2.0 and 4.0, and where *Fi* represents the force exerted by muscle *i*, and *PCSAi* its physiological cross-sectional area. In this case, *n* = 2 will be used. This function is minimised when subjected to Eq. 13 together with the force balances of the tendon nets. Additional constraints are that tendon forces must be non-negative, and the limits of muscle forces obtained from Eqs. 2 and 3 varying the muscle activation level from 0 to 1

$$F\_{\rm PE} \cdot F\_{\rm max} \le F \le \{F\_l \cdot F\_v + F\_{\rm PE}\} \cdot F\_{\rm max} \,. \tag{16}$$

### **3.6 Grasp evaluation**

230 Theoretical Biomechanics

2 2 <sup>2</sup> ·*P e*

coefficient and *en* is called the eccentricity parameter (height of the ellipsoid described by Eq. 11). Considering a Winkler elastic foundation (Johnson, 1985) of depth *h* and elastic

*baen* ·· <sup>15</sup>

where *a* and *b* can be calculated from the relative radii of curvature *R'* and *R''* of the objects

· ;''··2;'··2 *RRK hP RbRa*

δδ

The actual grasping forces for a given posture will be obtained by considering that they have to satisfy the dynamic equilibrium of the grasped object. There is not a unique set of forces that ensures the equilibrium but we have to take into account the biomechanical limitations (maximal muscle forces) and the control performed by the central nervous system. In an effort to minimise the computational cost, we propose to uncouple the computation of the contact forces from the neuromuscular control model. This can be done by considering that the central nervous system is trying to attempt performing the grasp with minimal contact forces, as implemented for robots in the work of Liu et al. (2004a).

The movement of the skeletal chains, together with the contact forces and the corresponding application points are input to the model. The problem to be solved is the derivation of the muscle activation levels required to produce the given motion under the external loads. It is,

The dynamics equations of the open chain of rigid bodies have been derived using the Lagrange method (García de Jalón & Bayo, 1994). For a system with *m* generalised co-

*k*

where *L* is the Lagrangian function and *nc Qk* are the generalised non-conservative forces. The generalised coordinates have been considered coincident with the system DOF (*m=23*). Eqs. 14 together with the force balances of the tendon nets lead to an indeterminate problem. For example, in the case of the index finger, there are 12 equations (four corresponding to the DOF considered and eight to force balances in the tendon net) and 18 unknowns (six muscle forces and 12 branch forces of the tendon net). There is not a unique combination of muscular efforts that satisfy the dynamic equilibrium constraints. To solve the problem, a criterion chosen by the central nervous system to determine the muscle action control must be introduced. Our proposal is to maximise the endurance. According to Crowninshield and

*dL L Qk m dt q q*

*k k*

Brand (1981), this is achieved by minimising the non-linear objective function

1, , *nc*

∂ ∂ −= = ∂ ∂ , (14)

of the elastic layer:

<sup>=</sup>

*n n*

τ

*f*

where *P* is the total load applied in the direction of the contact normal,

modulus *K*, the eccentricity parameter is given by:

δ

δ=

in contact and the compression

**3.5 Neuromuscular control** 

therefore, an inverse dynamics problem.

ordinates *qk*, this equation is expressed as:

*t*

22

≤+ , (11)

<sup>8</sup> <sup>=</sup> , (12)

= , (13)

( ) 2/1 '''···

π

μ

is the friction

μ

A global grasp evaluation can be performed through the use of the proposed model, merging the knowledge from ergonomics, robotics and biomechanics. The classical ergonomics evaluation of grasp posture and reachability is possible, for different percentiles of the population represented by the corresponding anthropometric parameters. Furthermore, CHTD evaluation can be performed by using the predicted postures and muscle forces.

It is advisable to use force closure from robotics as a part of the proposed model; once a grasping posture is estimated by the grasping algorithm, force closure should be assured before consuming time in determining the contact forces. Any of the robotics quality measures could be used for evaluating the grasp. Depending on the task to be performed, it would be better to use a grasp quality measure to evaluate the disturbance resistance or a grasp quality measure to evaluate the manipulability.

But the most relevant contribution to grasp evaluation has to come from biomechanics analysis. Grasp measures related to the muscle and articular forces have to be investigated. Just to provide insight into this sense, and to ensure coherence with our model formulation, we propose the use of Eq. 15 as a quality measure related to fatigue that we can call *fatigue index*: the smaller the fatigue index the better will be the grasp. For power grasps, an alternative measure can be the difference between the maximal force the hand can exert on the grasped object for the posture being analysed and the real contact forces estimated by the contact model; this alternative measure can be seen as a safety margin for the muscle forces, that we can call *muscle safety margin index*. Additional measures can be investigated, such as the maximal contact pressure, etc.

### **4. Conclusion**

A realistic and self-contained biomechanical model of the hand has been proposed by merging the current knowledge of biomechanics, ergonomics and robotics. The model simulates the complete hand and can easily be scaled to study different percentiles of populations. It has a realistic representation that allows the ergonomic evaluation of products. The model is dynamic and can be used to study the muscular patterns associated with a specific grasp. It allows predicting feasible grasping postures and provides the contact information required for evaluating the grasp. Finally, the model incorporates original quality grasping measures such as the fatigue index and the muscle safety margin index, in addition to the usual robotics and ergonomics metrics and evaluations. All the

Towards a Realistic and Self-Contained Biomechanical Model of the Hand 233

Bernabeu, J. & Tornero, J. (2002). Hough Transform for Distance Computation and Collision

Biagiotti, L., Lotti, F., Melchiorri, C. & Vassura, G. (2002). Design Aspects for Advanced

Bicchi, A. (1994). On the Problem of Decomposing Grasp and Manipulation Forces in

Bicchi, A. (1995). On the Closure Properties of Robotic Grasping. *International Journal of Robotics Research*, 14, 4, (August 1995), pp. (319–334), ISSN 0278-3649 Bicchi, A. & Kumar, V. (2000). Robotic grasping and contact: a review, *Robotics and* 

Biryukova, E.V., Yourovskaya, V.Z. (1994). A model of human hand dynamics. *Advances in* 

Borst, C., Fischer, M. & Hirzinger, G. (2003). Grasping the dice by dicing the grasp, *Proc.* 

Brand, P.W. & Hollister, A.M. (1992). *Clinical Mechanics of the Hand*, Elsevier Science Health

Buchholz, B., T.J. Armstrong, &S.A. Goldstein. Anthropometric data for describing the

Buchner, H.J., Hines, M.J., Hemami, H (1988). A dynamic model for finger interphalangeal coordination. *Journal of Biomechanics*, 21, 6,pp. (459-468), ISSN 0021-9290 Case, K., Porter, J.M. and Bonney, M.C., SAMMIE: a man and workplace modelling system.

Casolo, F., Lorenzi, V. (1994). Finger mathematical modelling and rehabilitation. En

Chao, E.Y., Opgrande, J.D., Axmear, F.E. (1976). Three dimensional force analysis of finger joints in selected isometric hand functions. *Journal of Biomechanics*, 9: 387-395. Chao, E.Y., An, K.N. (1978). Internal forces in human hand. *Journal of the Engineering* 

Chiu, S.L. (1988). Task compatibility of manipulator pustures. *Int. J. Robotics Research*, 7, 5)

Ciocarlie, M., Miller, A. & Allen, P. (2005). Grasp Analysis Using Deformable Fingers, 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2005) Ciocarlie, M., Lackner, C. & Allen, P. (2007). Soft Finger Model with Adaptive Contact

Geometry for Grasping and Manipulation Tasks, *Second Joint EuroHaptics Conference* 

*Computer Aided Ergonomics*, 1990, 31–56 (Taylor & Francis: London).

Science div, third edition, ISBN13 9780815127864, St. Louis

Towards Intelligent Robotic Manipulation., pp. (1–16)

13, 2, (July 1994), pp. (127-147), ISSN 0921-8890

(393-398), ISSN 1042-296X

ISBN: 0-306-44580-8, New York.

348–353.

*May 2004*

0014-0139

pp. (13-21)

New York. Plenum Press.

*Mechanics Division*, ASCE, 104 (EM1): 255-272.

Choi, J. (2008). Doctoral dissertation, University of Michigan.

Avoidance, *IEEE Transactions on Robotics and Automation*, 18, 3, (June 2002), pp.

Robot Hands, IEEE/RSJ *International Conference on Intelligent Robots and Systems*.

Multiple Whole-limb Manipulation, *Int. Journal of Robotics and Autonomous Systems,*

*Automation, 2000. Proceedings. ICRA '00. IEEE International Conference on*, IEEE, pp.

*the Biomechanics of the Hand and Wrist*, Eds. F. Schuind, pp. (107-122), Plenum Press,

*IEEE/RSJ IROS, ISBN 0-7803-7860-1, Las Vegas, NV, USA, Octuber 2003*Borst, C., Fischer, M. & Hirzinger, G. (2004). Grasp planning: how to choose a suitable task wrench space, *Proc. IEEE ICRA, ISBN 0-7803-8232-3, New Orleans, LA, USA, April-*

kinematics of the human hand. *Ergonomics* 35, 3, (March 1992), pp.(261-273), ISSN

Advances in the Biomechanics of the Hand and Wrist, Eds. F. Schuind, pp. 197-223,

abovementioned features are performed in a virtual environment, without external experimental data.

## **5. Acknowledgment**

We are grateful for the financial support of the Fundació Caixa-Castelló and the Universitat Jaume I throughout the project P1-1B2009-40; of the Spanish Research and Innovation Ministry, and the EU (FEDER funds) throughout the project DPI2010-18177; and of European Community's Seventh Framework Programme through the project GRASP, IST-FP7-IP-215821. This research is related to these three projects.

### **6. Nomenclature**


### **7. References**


abovementioned features are performed in a virtual environment, without external

We are grateful for the financial support of the Fundació Caixa-Castelló and the Universitat Jaume I throughout the project P1-1B2009-40; of the Spanish Research and Innovation Ministry, and the EU (FEDER funds) throughout the project DPI2010-18177; and of European Community's Seventh Framework Programme through the project GRASP, IST-

FP7-IP-215821. This research is related to these three projects.

Cumulative hand trauma

EDQ Extensor digiti quinti

**7. References** 

0021-9290

Ltd., New York, pp. 1-50.

3D Three-dimensional FP Flexor profundus ADD Adductor pollicis FPB Flexor pollicis brevis ADQ Abductor digiti quinti FPL Flexor pollicis longus APB Abductor pollicis brevis FS Flexor superficialis APL Abductor pollicis longus GWS Grasp wrench space CE Contractile element HB Hand breadth

disorders HL Hand length CMC Carpometacarpal IP Interphalangeal DI Dorsal interosseous LU Lumbrical

DIP Distal interphalangeal MCP Metacarpophalangeal DOF Degrees of freedom OPP Opponens pollicis EDC Extensor digitorum communis OWS Object wrench space

EFM Elastic foundation model PEE Parallel elastic element EI Extensor indicis PIP Proximal interphalangeal EPB Extensor pollicis brevis SEE Series elastic element EPL Extensor pollicis longus TWS Task wrench space FDQ Flexor digiti quinti VI Volar interosseous

An, K.N., Chao, E.Y., Cooney, W.P. & Linscheid, R.L. (1979). Normative model of human

An, K.N., Chao, E.Y.S. and Kaufman, K.R., 1991. Analysis of muscle and joint loads. In:

Armstrong, T.J., Best, C., Bae, S., Choi, J., Grieshaber, D.C., Park, D., Woolley, C. & Zhou, W.

hand for biomechanical analysis. J*ournal of Biomechanics*, 12, 10, pp. (775-788), ISSN

Mow, V.C and Hayes, W.C. (Eds.), *Basic Orthopaedics Biomechanics*. Raven Press.

Development of a Kinematic Hand Model for Study and Design of Hose Installation. *Digital Human Modeling*, HCII 2009, LNCS 5620, pp. 85–94, 2009.

 PCS A

area

Physiological cross-sectional

experimental data.

**6. Nomenclature** 

CHT D

**5. Acknowledgment** 


Towards a Realistic and Self-Contained Biomechanical Model of the Hand 235

Goldstein, S.A., Armstrong, T.J., Chaffin, D.B. & Matthews, L.S. (1987). Analysis of

Gonthier, Y. (2007). Contact Dynamics Modelling for Robotic Task Simulation. Ph. d. thesis.

Goussous, F.A. Grasp planning for digital humans. PhD dissertation. Iowa University. 2007 Harada, K., Kaneko, M. & Tsuji, T. (2000). Rolling Based Manipulation for Multiple Objects. *Proc. IEEE Int. Conf. on Robotics and Automation (ICRA),* pp. 3887-3894. Haschke, R. Steil, J.J., Steuwer, I. & Ritter, H. (2005). Task-oriented quality measures for

Hatze, H., 1981. Myocibernetic control models of skeletal Muscle. University of South

Hertz, H. (1882). On the Contact of Rigid Elastic Solids and on Hardness, Ch 6: Assorted

Hill, A.V., 1938, "The of shortening and dynamic constants of muscle", Proc R Soc London

Hollister, A.M., Giurintano, D.J., Buford, W.L., Myers, L.M. & Novich, A. (1995). The axes of

Howe, R. Kao, I. & Cutkosky, M. (1988). Sliding of robot fingers under combined tornsion

Howe, R. & Cutkosky, M. (1996). Practical force-motion models for sliding manipulation *Int.* 

Jacobson, M.D., Raab, R., Faxelo. B.M., Abrams, R.A., Botte, M.J. & Lieber, R.L. (1996).

Kamper DG, Fischer HC, Cruz EG. (2006). Impact of finger posture on mapping from muscle activation to joint torque. Clinical Biomechanics. 21(4):361-9. Kao, I. & Cutkosky, M (1992). Dexterous manipulation with compliance and sliding. *Int.* 

Kao, I. & Yang, F. (2004). Stiffness and Contact Mechanics for Soft Fingers in Grasping and Manipulation. *IEEE Transactions on Robotics and Automation*, 20(1), pp. 132-135. Kao, I., Linch, K. & Burdick, J. W. (2008). Contact Modeling and Manipulation. *Handbook of* 

Kaufman, K.R., An, K.N., Litchy, W.J. and Chao, E.Y.S., 1991. Physiological prediction of muscle forces I. Theoretical formulation. J. Neuroscience 40 (3), 781-792. Kawaguchi, K., Endo, Y., Kanai, S. (2009). Database-Driven Grasp Synthesis and Ergonomic

Keyserling, W.M. 2000; Workplace risk factors and occupational musculoskeletal disorders,

Assessment for Handheld Product Design. Lecture Notes in Computer Science,

part 2: A review of biomechanical and Psychophysical research on risk factors

Kapandji, A.I., 1998. *Fisiologie articulaire. Membre Supérieur*. Editions Maloine, Paris.

rotation of the thumb interphalangeal and metacarpophalangeal joints. *Clinical* 

and shear loading. Proceedings *IEEE Int. Conf. on Robotics and Automation*, pp. 103-

Architectutural design of the human intrinsic hand muscles. *The Journal of Hand* 

1-6.

105.

University of Waterloo.

*and Automation*, pp. 689-694.

Papers by H. Hertz . MacMillan, New York, 1882.

*Journal of Robotics Research* 15(6), pp. 555-572.

*Journal of Robotics Research* 11(1), pp. 20-40.

Johnson, K. L., (1985). Contact Mechanics. Cambridge University Press.

Africa, Pretoria.

B, Vol 126, pp. 136-195.

*Orthopaedics*, (320), pp. 188-193.

*Surgery*, 17A (5), pp. 804-811.

*Robotics, chapter 27*, pp. 647-668.

Volume 5620/2009, 642-652.

cumulative strain in tendons and tendon sheaths. *Journal of Biomechanics*, 20 (1), pp.

dextrous grasping. In: *Proc. 6th IEEE Int. Conf. Computational Intelligence in Robotics* 

*and Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systems (WHC'07),* ISBN 978-0-7695-2738-3, Tsukuba, Japan, March 2007


Colombo and Cugini, 2005; Virtual humans and prototypes to evaluate ergonomics and

Craig, S.M., 1992. Anatomy of the joints of the fingers. *Hand Clinics,* 8, 4, (November 1992),

Crowninshield, R.D. & Brand, R.A. (1981). A physiologically based criterion of muscle force

Cutkosky, M. (1989). On grasp choice, grasp models, and the design of hands for

Davidoff and Freivalds, 1993) Davidoff N.A., Freivalds A. A graphic model of the human

Denavit, J. and Hartenberg, R.S. (1955) 'A kinematic notation for lower-pair mechanisms based on matrices', *Journal of Applied Mechanics*, Vol. 77, pp.215–221. Diankov, R. (2010). Automated Construction of Robotic Manipulation Programs, PhD thesis,

Ding, D., Lee, Y.-H. & Wang, S. (2001). Computation of 3-d form-closure grasps, IEEE

Dubousset, J.F., 1981. The digital joints. In: Tubiana, R. (Ed.), The Hand. Saunders

Endo, Y., Kanai, S., Kishinami, T., Miyata, N., Kouchi, M., Mochimaru, M. (2007). Virtual

Endo, Y., Kanai, S., Miyata, N., Kouchi, M., Mochimaru, M., Konno, J., Ogasawara, M.,

Erdmann, M.A. (1994). On the representation of firction in configuration space. *Int. Journal of* 

Esteki, A., Mansour, J.M. (1997). A dynamic model of the hand with application in

Eyler, D.L. and Markee, J.E., 1954, "The anatomy and function of the intrinsic musculature of the fingers", *The Journal of Bone and Joint Surgery,* Vol. 36-A, pp. 1-18. Ferrari, C. & Canny, J. (1992). Planning optimal grasps, *Proceedings 1992 IEEE International* 

Fok KS, Chou SM. Development of a finger biomechanical model and its considerations.

García de Jalón, J. and Bayo, E., 1994. Kinematic and dynamic simulation of multibody

Giurintano, D.J. & Hollister, A.M. (1991). Force analysis of the thumb for a five-link system.

Giurintano, D.J., Hollister, A.M., Buford, W.L., Thompson, D.E. & Myers, L.M. (1995). A

virtual 5-link model of the thumb. *Medical Engineering & Physics,* 17, 4, (June 1995),

Ergonomics Conference: Celebrating the Past - Shaping the Future.

grasping assessment using 3D digital hand model. 10th Annual Applied

Shimokawa, M. (2009). Optimization-Based Grasp Posture Generation Method of Digital Hand for Virtual Ergonomics Assessment. SAE Int. J. Passeng. Cars -

functional neuromuscular stimulation. *Annals of Biomedical Enginnering*, 25: 440-451.

*(WHC'07),* ISBN 978-0-7695-2738-3, Tsukuba, Japan, March 2007

prediction in locomotion. *Journal of Biomechanics*. 14: 793-801.

Carnegie Mellon University, Robotics Institute.

Electron. Electr. Syst. Volume 1, Issue 1, 590-598.

*Conference on Robotics and Automation* pp. 2290–2295.

*Journal of Biomechanics*. 2010 Mar 3;43(4):701-13.

ASME Biomechanics Symposium, 120:213-217.

systems. Springer-Verlag, New York.

pp. (297-303), ISSN 1350-4533

Company, Philadelphia, pp. 191-201.

*Robotics Research,* 13(3),pp. 240-271.

Transactions on Robotics and Automation **17**(4): 515 –522.

4828

pp. (693-700), ISSN 0749-0712

*and Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systems* 

safety. *Journal of Engineering Design,* 16, 2, (April 2005), pp. (195–207), ISSN 0954-

manufacturing tasks, IEEE Transactions on Robotics and Automation **5**(3): 269–279.

hand using CATIA. *International Journal of Industrial Ergonomics*, 1993, 12 (4) 255-264


Towards a Realistic and Self-Contained Biomechanical Model of the Hand 237

Mason, M. T. (2001). *Mechanics of robotic manipulation*, The MIT Press, Cambridge, MA, USA. Miller, A. T. & Allen, P. K. (1999). Examples of 3D grasp quality computations, *Robotics and* 

Miller, A. T., Knoop, S., Christensen, H. & Allen, P. K. (2003). Automatic grasp planning

Miller, A., Allen, P.K. (2004). "Graspit!: A Versatile Simulator for Robotic Grasping". *IEEE* 

Minami, A., An, K.N., Cooney, W.P. and Linscheid, R.L., 1985. Ligament stability of the MCP joint: A biomechanical study. *J. Hand Surgery* 10A (2), 255-260. Mirthich, B. & Canny, J. (1994). Easily computable optimum grasps in 2D and 3D. In: *Proc.* 

Mishra, B., Schwartz, J. T. & Sharir, M. (1987). On the existence and synthesis of multinger

Mital, A. and Kilbom, A., 1992, Design, selection and use of hand tools to alleviate trauma of

Mommersteeg, T.J.A., Blankevoort, L., Huiskes, R., Kooloos, J.G.M. and Kauer, J.M.G., 1996.

Morales, A., Sanz, P., Del Poblil, A. & Fagg, A. (2006). Vision-based three-finger grasp

Muggleton, J.M., Allen, R. & Chappell, P.H. Hand and arm injuries associated with

Parada, J. E., Nava, N. E. & Ceccarelli, M. (2008). A Methodology for the Design of Robotic

Park, Y.C. & Starr, G.P. (1992). Grasp synthesis of polygonal objects using a three-fingered robotic hand. *International Journal of Robotics Research,* 11 (3), pp. 163-184. Peña-Pitarch, E. Virtual human hand. PhD dissertation. Universitat Politècnica de

Pollard, S.S. (2004). Closure and quality equivalence for efficient synthesis of grasps from

Ponce, J., Sullivan, S., Sudsang, A., Boissonat, J.D. & Merlet, J.P. (1997). On computing for-

Qiu D, Fischer HC, Kamper DG. Muscle activation patterns during force generation of the

Rezzoug, N. & Gorce, P. (2008). Prediction of fingers posture using artificial neural

index finger. Conf Proc IEEE Eng Med Biol Soc. 2009;2009:3987-90.

examples. *Int. J. Robotics Research*, 23 (6), pp. 595-614.

networks. *Journal of Biomechanics*. 41, pp. 2743-2749.

*Journal of Robotics Research*, 16 (1), pp. 11-35.

the upper extremities: Part II- The scientific basis (knowledge base) for the guide,

Characterization of the mechanical behavior of human knee ligaments: a

synthesis constrained by hand geometry. *Robotics and Autonomous Systems,* 54, 6,

repetitive manual work in industry: a review of disorders, risk factors and

Hands with Multiple Fingers, *International Journal of Advanced Robotic Systems* **5**(2):

finger equilibrium and force-closure grasps of polyhedral objects. *International* 

IEEE International Conference on, Vol. 2, IEEE, pp. 1824–1829.

*Robotics and Automation Magazine*, V. 11, No.4, pp. 110-122.

*International Journal of Industrial Ergonomics*, Vol. 10, pp. 7-21

numerical-experimental approach. J. Biomechanics 29 (2), 151-160.

preventive measures. *Ergonomics*, 1999, VOL. 42, NO. 5, 714 - 739 Napier, J. (1956). The prehensile movements of the human hand, Surger **38**(4): 902–913. Nguyen, V.-D. (1988). *Constructing force-closure grasps*, Institute of Electrical and Electronics

1240–1246.

Engineers.

177–184.

Catalunya. 2007.

*IEEE ICRA 1994*, pp. 739-747.

positive grips, *Algorithmica* **2**(1-4): 541–558.

(June 2006), pp. (496–512), ISSN 0921-8890

*Automation, 1999. Proceedings. 1999 IEEE International Conference on*, Vol. 2, IEEE, pp.

using shape primitives, Robotics and Automation, 2003. Proceedings. ICRA'03.

associated with upper extremity disorders. AIHA Journal; Mar/Apr 2000; 61 (2) pg. 231-243


Kong, Y., Jang, H. & Freivalds, A. (2006) Wrist and Tendon Dynamics as Contributory Risk

Kumar S. 2001. Theories of musculoskeletal injury causation. Ergonomics, 2001, VOL. 44,

Kurita Y, Onoue T, Ikeda A, Ogasawara T. Biomechanical analysis of subjective pinching

Kyota, F., Watabe, T., Saito, S. & Nakajima, M. (2005). Detection and Evaluation of Grasping

Lee, J.W., Rim, K. (1990). Maximum finger force prediction using a planar simulation of the

Lee SW, Chen H, Towles JD, Kamper DG. (2008a) Effect of finger posture on the tendon

Lee SW, Chen H, Towles JD, Kamper DG. (2008b). Estimation of the effective static moment

Leijnse, J.N.A.L., Bonte, J.E., Landsmeer, J.M.F., Kalker, J.J., Van Der Meulen, J.C., Snijders,

Leijnse, J.N.A.L., Kalker, J.J. (1995). A two-dimensional kinematic model of the lumbrical in

Lemay, M.A. & Crago, P.E. (1996). A dynamic model for simulating movements of the elbow, forearm and wrist. *Journal of Biomechanics*, 29 (19, pp. 1319-1330. Li, Z & Sastry S. (1988). Task-oriented optimal grasping by multifingered robotic hands.

Li Z., Chang C.-C., Dempsey P. G., Ouyang L., Duan J.. Validation of a three-dimensional

Liegeois, A. (1977). Automatic supervisory control for the configuration and behavior of

Lin, Q., Burdick, J. M. & Rimon, E. (2000). A stiffness-based quality measure for compliant

Liu, G., Xu, J., & Li, Z. (2004a). On geometric algorithms for real-time grasping force optimization. *IEEE Transactions on Control Systems* Technology, 12 (6): 843–859. Liu, G., Xu, J., Wang, X. & Li, Z. (2004b). On Quality Functions for Grasp Synthesis, Fixture

Mansour, J.M., Rouvas, C., Sarangapani, J. (1994). Quantitative functional anatomy of finger

*and Wrist*, Eds. F. Schuind, pp. 177-188, New York. Plenum Press.

hand scanning and dimension extraction method with dimension data.

multibody mechanisms. *IEEE Trans. System, Man and Cybernetics*, 7 (12), pp. 842-

grasps and fixtures. *IEEE Transactions on Robotics and Automation*, 16(6), pp. 675-688.

Planning, and Coordinated Manipulation, *IEEE Transactions on Automation Science* 

muscles: Appliction to controlled grasp. *En Advances in the Biomechanics of the Hand* 

*in Manufacturing*, 16, 1, pp. (83–105), ISSN 1090-8471

C.J. (1992). *Journal of Biomechanics*, 25 (11):1 253-1264.

the human finger. *Journal of Biomechanics,* 28 (3): 237-249.

*IEEE Journal of Robotics and Automation*, 4 (1), pp. 32-44.

Ergonomics, Vol. 51, No. 11, November 2008, 1672–1692

231-243

NO. 1, 17 ± 47

2009;2009:5231-4.

204: 169-178.

41(7):1567-73.

868.

*Engineering*. 130(5):051014.

*and Engineering* **1**(2): 146–162.

Conference on Cyberworlds (CW'05).

associated with upper extremity disorders. AIHA Journal; Mar/Apr 2000; 61 (2) pg.

Factors in Work-Related Musculoskeletal Disorders. *Human Factors and Ergonomics* 

effort based on tendon-skeletal model. Conf Proc IEEE Eng Med Biol Soc.

Positions for Autonomous Agents. Proceedings of the 2005 International

middle finger. Proc. Instn. Mech. Engnrs. Part H: *Journal of Engineering in Medicine*,

force distribution within the finger extensor mechanism. *Journal of Biomechanical* 

arms of the tendons in the index finger extensor mechanism. *Journal of Biomechanics*.

Mason, M. T. (2001). *Mechanics of robotic manipulation*, The MIT Press, Cambridge, MA, USA.


Towards a Realistic and Self-Contained Biomechanical Model of the Hand 239

Storace, A., Wolf, B. (1982). Kinematic analysis of the role of the finger tendons. *Journal of* 

Thomas, D.H., Long, C., Landsmeer, J.M.F. (1968). Biomechanical consideration of lumbricalis behaviour in the human finger. *Journal of Biomechanics*, 1: 107-115.

Tubiana, R., Thomine, J. and Mackin, E. (1996), 'Examination of the hand and wrist', Second

Valero-Cuevas, F.J., Zajac, F.E., Burgar, C.G. (1998). Large index-fingertip forces are

Valero-Cuevas, F.J. (2000). Predictive modulation of muscle coordinatin pattern magnitude

Valero-Cuevas, F.J. Towles, JD, Hentz, VR. (2000). Quantification of fingertip force reduction

Valero-Cuevas, F.J. (2005). An integrative approach to the biomechanical function and neuromuscular control of the fingers. *Journal of Biomechanics*, 38(4): p.673-84. van Nierop O.A., van der Helm A., Overbeeke K.J., Djajadiningrat T.J.P. A natural human

Vigouroux L, Quaine F, Labarre-Vila A, Moutet F. (2006) Estimation of finger muscle tendon

Vigouroux L, Ferry M, Colloud F, Paclet F, Cahouet V, Quaine F. Is the principle of

Weightman, B. & Amis, A.A. (1982). Finger joint force predictions related to design of joint

Wren, D. & Fisher, R. (1995). Dextrous hand grasping strategies using preshapes and digit

Wu JZ, An KN, Cutlip RG, Dong RG. A practical biomechanical model of the index finger

Xydas, N. & Kao, I. (1999). Modeling of contact mechanics and friction limit surface for soft

Xydas, N., Bhagavat, M. & Kao, I. (2000). Study of soft-nger contact mechanics using nite

Yang, J., Kim, J.H., Abdel-Malek, K., Marler, T., Beck, S. & Kopp, G.R. A new digital human

Yoshikawa, T. (1985). Manipulability of robotics mechanisms. *Int. J. Robotics Research*, 4, 2,

produced by subject-independent patterns of muscle excitation. *Journal of* 

scales fingertip force magnitude over the voluntary range. *Journal of* 

in the forefinger following simulated paralysis of extensor and intrinsic muscles.

tensions and pulley forces during specific sport-climbing grip techniques. *Journal of* 

minimization of secondary moments validated during various fingertip force

trajectories, IEEE International Conference on Systems, Man and Cybernetics, Vol.

simulating the kinematics of the muscle/tendon excursions. *Biomed Mater Eng.*

fingers with experimental results, *Int. Journal of Robotics Research*, 18(9), pp. 947-950.

elements analysis and experiments, *Proceedings 2000 ICRA. Millennium Conference. IEEE International Conference on Robotics and Automation. Symposia Proceedings (Cat.* 

environment and assessment of vehicle interior design. Computer-Aided Design

*Biomechanics*, 15, 5, pp. (391-393), ISSN 0021-9290

Tubiana, R (1981), 'The hand', Vol. 1, W.B. Sanders Company, 1981.

Edition, Martin Dunitz, 1996.

*Neurophysiology*, 83 (3): 1469-1479.

*Journal of Biomechanics*, 33(12):1601-1609.

hand model. Visual Comput (2008) 24: 31–44

*Biomechanics*, 39, 14, pp. (2583-2592), ISSN 0021-9290

production conditions? *Hum Mov Sci*. 2008 Jun;27(3):396-407.

replacements. *Journal of Biomedical Engineering.* 4:197-205.

*Biomechanics*, 31: 693:703.

1, pp. 910 –915 vol.1.

*No.00CH37065)* (April): 2179–2184

2010;20(2):89-97.

39 (2007) 548–558

pp. (735-747)


Rhee T., Neumann U., Lewis J.P. 2006, Human Hand Modeling from Surface Anatomy. I3D

Roa Garzón, M. (2009). *Grasp Planning Methodology for 3D Arbitrary Shaped Objects*, Ph. d.

Rogers M. S., Barr A. B., Kasemsontitum B.; Rempel D. M.. A three-dimensional

Roloff I, Schöffl VR, Vigouroux L, Quaine F. Biomechanical model for the determination of

Salisbury, J. (1982). *Kinematics and Force analysis of Articulated Hands*, Ph.d. thesis, Standford

Sancho-Bru, J.L., 2000, "Model Biomecànic de la mà orientat al disseny d'eines manuals",

Sancho-Bru, J.L., Perez-Gonzalez, A., Vergara-Monedero, M., Giurintano,D.J., (2001). A 3-D

Sancho-Bru, J.L., Giurintano, D.J., Pérez-González, A., and Vergara, M., (2003b), "Optimum tool handle diameter for a cylinder grip", *J. Hand. Ther*., 16(4), pp. 337-342. Sancho-Bru, J.L, Vergara, M., Rodríguez-Cervantes, P.J., Giurintano, D., Pérez-González, A.

Schoenmarklin R.W., Marras W.S., Leurgans S.E. industrial wrist motions and incidence of hand/wrist cumulative trauma disorders. *Ergonomics*, 1994; 37: 1449-1459 Siciliano, B. & Khatib, O. (eds) (2008). *Springer Handbook of Robotics*, Vol. 15, Springer, Berlin,

Smith, E.M., Juvinall, R.C., Bender, L.F. Pearson, J.R. (1964). Role of the finger flexors in

Spielholz et al., 2001; P. Spielholz, B. Silverstein, M. Morgan, H. Checkoway And J.

Spoor, C.W., Landsmeer, J.M.F. (1976). Analysis of the zigzag movement of the human

Spoor, C.W. (1983). Balancing a force on the fingertip of a two-dimensional finger model

Stansfield, S. (1991). Robotic grasping of unknown objects: A knowledge-based approach,

Storace, A., Wolf, B. (1979). Functional role of the finger tendons. *Journal of Biomechanics,* 12:

without intrinsic muscles. *Journal of Biomechanics*, 16: 497-504.

*The International Journal of Robotics Research* **10**(4): 314–326.

*Biomechanics*, 34, 11, (November 2001), pp. (1491–1500), ISSN 0021-9290 Sancho-Bru, J.L., Perez-Gonzalez, A., Vergara, M., Giurintan, D.J., (2003a). A 3D

ACM New York

23.

University.

0090-6964

Heidelberg.

575-578.

*Engineering* 125, 78–83.

thesis, Universidad Politécnica de Cataluña.

*Ergonomics*, Vol. 51, No. 4, April 2008, 511–526

Ph.D. thesis, Universitat Jaume I, Castelló.

Ergonomics, 2001, VOL. 44, NO. 6, 588 - 613

tendon. *Journal of Biomechanics*, 9: 561-566.

'06 Proceedings of the 2006 symposium on Interactive 3D graphics and games,

anthropometric solid model of the hand based on landmark measurements.

the forces acting on the finger pulley system. *Journal of Biomechanics.* 2006;39(5):915-

dynamic model of human finger for studying free movements. *Journal of* 

biomechanical model of the hand for power grip. *Journal of Biomechanical* 

(2008). Scalability of the Muscular Action in a Parametric 3D Model of the Index Finger. *Annals of Biomedical Engineering,* 36, 1, (January 2008, pp. (102-107), ISSN

rheumatoid deformities of the metacarpophalangeal joints. *Arth. Rheum*., 7: 467-480.

Kaufman. Comparison of self-report, video observation and direct measurement methods for upper extremity musculoskeletal disorder physical risk factors.

finger under influence of the extensor digitorum tendon and the deep flexor


**Part 3** 

**Locomotion and Joint Biomechanics** 

