Abstract

The thumb prosthesis mechanism is optimally designed by using five performance criteria including the following: least square structural error, mechanical manufacturing imprecision error, driving optimal torque, mechanical strength reliability, and production cost of the thumb mechanism. This paper was devoted to the optimization of the thumb prosthesis’s mechanism by taking into consideration the manufacturing cost model based on machining cost theory which took into detail the shape of the workpieces and the strength reliability of all the parts composing the entire mechanism. Every optimization problem displays a particular set of an independent vector of optimal parameters, showing the impact of each objective function on the configuration of the prosthetic device. The multiobjective optimization showed that the mechanical reliability and the production cost included in any combination of the simultaneous optimization enabled the achievement of the same optimum variables design, though with some exceptions. With the inclusion of the labor charges, the depreciation rate of the equipment, and production assets in the mathematical’s manufacturing cost model, the optimal manufacturing cost generated from the numerical simulation was 501.0021 USD. Therefore, the global manufacturing cost and the mechanical strength reliability of the whole prosthesis mechanism have a real impact on the customization of the structure, due to the stochastic nature of the trajectory of the cutting tools during the manufacturing processes.

1. Introduction

The invention of new prosthetic devices over the past decades has observed advances as a consequence of recent technological development, advancing in the aspect of greater dexterous hand devices. The use of a mechanical prosthesis to replace a missing limb is one of the most fascinating challenges faced by prosthetists in the prosthetic rehabilitation protocol [14]. The level of autonomy and the capabilities of performing social activities, trade, and daily tasks negatively influenced amputees. However, the state-of-the-art prosthesis failed to coalesce smooth functioning, reliability, suitable personalized aesthetic, adequate maintenance, affordability, and sometimes adequate grasping [5, 6]. Most prostheses fail to closely match the user-specific requirement and have poor design, resulting in the abandonment of the artificial device by the amputees [7]. One of the most popular works on prosthesis optimization in the literature has been done by Vinet et al. [8]. They obtained all the bending angles of the hand morphology using 3D graphical software CATIA on four planar bar mechanisms.

The thumb plays a paramount role in the various grasping modes. Liken to unfavorable conditions such as sickness, accident, or congenital malformation, the rehabilitation process of the thumb is greatly affected especially by its contribution to six grasping modes retained by many prosthetists [9, 10]. The design of a thumb prosthetic device is critical since forty percent of the entire functionality of the human depends on the hand. In housekeeping and laundry, the grasping mode is dominated at sixty percent by the contribution of the thumb. Many prosthetists faced the challenge to choose a simple, accurate, and less mobility or degree of freedom, which allows equal force distribution. The least square structural error on the second phalanx of the thumb and the trajectory of the model is determined by the end reflector which is the fingertip [11, 12]. To obtain the safe, smooth, accurate, and affordable functionality of the thumb, one of the most interesting solutions is the use of the planar articulated crossed four-bar mechanism to design the prosthetic mechanism as shown in Figure 1. The synthesis followed by the design methodology of mechanisms has been the subject of several investigations [14, 8, 13].

The main objective of this paper is to introduce into upper limb optimization processes the mechanical reliability of each element of the basic mechanism, used to design the artificial device. Their applications are numerous and their synthesis is generally done by reference to input positions and output positions of the considered parameters [14]. The four-link articulated planar mechanism proposed for synthesis finds application in a hand prosthesis that can perform various types of grasping daily objects. Its optimal synthesis will be made by reference to positions of the coupling and output members. For this purpose, two performance criteria are used: the least square structural error on the angle of flexion of the second phalanx of a thumb and the maximum motor torque necessary to balance a given force at the tip of this finger. This paper will also integrate the global manufacturing cost model taking into account the stochastic nature of the trajectory of the cutting tools during the machining process on the numerical control machine center. This cost will also integrate the social charges, the labor charges, and the depreciation rate of the equipment and productivity asset [15]. All these parameters were adjustable according to the economical standard of the area of the production. In order to study the virtual optimal mechanism, the mathematical model of the objective function relative to the various performance criteria is established [16]. Assuming that in the biomechanics of artificial limb, all the fingers contribute actively to accomplish if necessary each of the six grasping object modes retained by Pelletier and Vinet, Ruben, and Carrozza et al., taking into account the fact that the humanoid thumb mobility in particular among other fingers, the optimal design will be carried out on the base of minimization or maximization of single or a set of objectives function submitted to the various constraints [1719]. Finally, the entire optimization process will be necessary to appreciate the influence of the said criteria on the prosthesis’s human hand design. After the analysis of the mechanism consisting of the development of all the mathematical relations necessary to be summited to the optimization processes, for its synthesis, we will make the optimal synthesis of the said mechanism by reference to seven positions. Using the data relating to the thumb, the results obtained will be compared and discussed according to the state of art [20].

To mechanically simulate the different grasping mode of the object by an artificial hand, planar articulated mechanisms with one simple degree of freedom are generally recommended in hand prosthetic applications [2125]. Based on this configuration, 1DOF four-crossed bar mechanism is considered in this study to design the thumb which is usually optimized by the means of transmission angle and energy consumed which have a significant effect on the mechanism design. To improve the use of these devices for the patients in need, new additional criteria will be modeled and integrated with the optimization process as objective functions. The first is the manufacturing cost, which would impact the affordability of the prosthesis by the amputees, and the second is the reliability of the setting mechanism and its effect on the prosthesis design. In the first approach, each criterion will be optimized individually and the three first criteria will be optimized simultaneously. Strength reliability and manufacturing cost will be considered individually, and finally, we will carry out the optimization of all the five criteria simultaneously [20].

2. Material and Methods

The thumb’s mechanism description, its performance criteria, and its objective functions are presented in this section.

2.1. Thumb Finger Description

The thumb finger associated with the four-bar mechanism fit out with its driving system is depicted in Figure 1.

Figure 1 is made of a cable connected both to the driving link at the point and to the driving pulley at the point whose connection angle is . The structure of thumb mechanism is added in the revised manuscript as Figure 2.

Table 1 defines the different parameters shown in Figure 1.

Freudenstein’s relation [2325, 28] between and at any mechanism position is given by where , , , , , , , and .

The values of some parameters in Equation (1) are similar to those used in [24] and are given in Table 2.

The design parameter value bounds of the thumb’s mechanism are taken in [14] and are given in Table 3.

The vector of the independent variables of the design is defined as since the length of is fixed and .

In this study, the weight of the thumb mechanism was assumed to be negligible due to the dimensionless bars constituting the four-crossed bar mechanism used to model this artificial thumb [20].

2.2. Performance Criteria and Objective Functions of the Thumb Prosthesis Mechanism

The thumb prosthesis design is optimized by using the following five performance criteria: (i)Least square structural error on the bending angle of the second phalanx should be minimized to ensure that the prosthetic finger should align at the prescribed positions(ii)Mechanical error on the bending angle of the second phalanx, due to dimensional tolerances and clearances on the articulations, should be minimized(iii)Maximum shaft driving torque should be minimized and able to counterbalance a grasping force applied at the end of the thumb finger(iv)Strength reliability of the mechanism should be maximized to ensure the mechanical strength during grasping, holding, or pinching operations(v)Manufacturing cost of the mechanism should be minimized to evaluate its effect on the design variables

The objective function associated with the least square structural error on bending angle is defined as where is the desired angle and is the computed angle. It is also called mean square error, which describes the error, given in the form of least square function, describing the error observed between the computed value and the given value of the bending angle between the palmer plan and the second phalanx TP2 of the thumb [14].

The objective function associated with the mechanical error on bending angle is defined as

By assuming that the frictional forces at joints are negligible, the objective function associated with the maximum shaft driving torque is defined as where is the traction angle between the flexible cable on the driving system and the driving bar, is the driving torque applied to the shaft, and is torque reduced to point of the driving bar. The objective function associated with the strength reliability of the mechanism is defined as

where , , and are the reliability of driven, junction, and driving bars, respectively [28, 29]. Since the prostheses are manufactured in very small sets, the manufacturing computer-aided design is appropriate for this purpose [30]. The manufacturing unitary cost of the th manufacturing operation on the th bar of the mechanism is where , , , and are the sum of the expenses related directly to the preparation, cutting duration, out of cutting time, and the cutting tool, respectively; is the cost cutting tools; is the global rate of exploitation of the preparation section, general expenses, and labor; is hourly machine cost; , , , , , and are technical amortization, financial expenses, maintenance expenses and repair, expenses on local or on clutter, expenses on energy, and wage costs and social, respectively; , , , and are linear speed, mass production set, feed rate in turn/tooth, and Taylor constant function of material of the cutting tool, respectively; are the cutting tools used relatively at th bar and th manufacturing operation (milling cutter 2 cuttings ∅20, cutter 2 cuttings ∅6, drill ∅2, drill ∅4, piloting milling cutter to counter boring and machining tapping ); are the tool trajectories relative to the facing milling, lateral milling, drilling, counter boring, and tapping of the th bar of the mechanism and the th manufacturing operation, respectively; , , and are the holes 1, 2, and 3, respectively, and is execution duration allowed to the preparation of a set of workpiece. The manufacturing cost of the junction bar is where is the facing milling, is the lateral milling, is the drilling, and is the counterboring. The manufacturing cost of the driving bar is where is the facing milling, is the lateral milling, is the drilling, and is the counterboring. The manufacturing cost of the driven bar is where is the facing milling, is the lateral milling, is the drilling, is the counterboring, and is the tapping. The global manufacturing cost of the mechanism is with where is in USD and the dimensions of all mechanism pieces are in millimeter [26]. The objective functions related to individual and multiobjective optimization are

Each of these multiobjective functions consists of a set of objective functions which are simultaneously optimized. The above functions expressed in terms of design variables are equality constraints and inequality constraints according to the MATLAB optimization toolbox. The convergence criteria on the design variables , the objective functions , and the constraints are, respectively, equal to 10−7, 10−9, and 10−13.

3. Results and Discussions

Three arbitrary points are chosen to conduct the optimization process to validate the results generated during the numerical simulations. Each objective function (known as least square structural error, mechanical error and strength reliability are dimensionless, maximum driving torque is expressed in N.mm, and manufacturing cost is expressed in USD) which is subjected to three equalities and eighteen inequality constraints is optimized by taking into consideration the design variable bounds in Table 3. Table 3 also gives the unit relative to each design variable. By using the optimum absolute value of every single function, the multiobjective optimizations , , , and subjected to the same constraints are conducted. The optimization problem solved in this paper has a complex nature. We have five nonlinear objective functions subjected to all type of constraints and must be minimize simultaneously. There is no unique solution to the problem of multiobjective optimization of several objective functions expressed by mathematical expressions which are performance criteria. The selection of optimal solutions is based on Pareto-optimality conditions [31, 32]. A solution is called “pareto-optimal” if there is no other solution for which an improvement of an objective function does not lead to a degradation of at least one of the other functions. It is therefore important, for a problem of simultaneous minimization of several criteria, to generate the “pareto-optimal” solutions. It is what we achieve with “fgoalattain,” which solve multiobjective goal attainment problems. This method is implemented in “MATLAB” software using the “fgoalattain” minimization function. We have chosen the set of vector of the weighting coefficients where absolute (goal) is the absolute value of the vector of the objectives set from individual optimization [33]. This condition guarantees the same percentage of achievement of the objectives set for all active objective functions [20]. The optimal design variables of the thumb mechanism are presented in Table 4.

Table 4 contains all the generated values of the design variables resulting from the optimization of each objective function and a set of coupled functions according to the optimization strategy adopted. The vectors of transversal, longitudinal, or angular dimensions of the virtual mechanisms are proposed at each optimum obtained. Individual optimization is represented by to , and multiobjective optimization is represented by to . The units of the various objective function are indicated in Table 5 as follows: and are in degree Celsius, do not have unit, is expressed in N.mm, and is given in USD. The set of optimal design variable validated at each optimization problem corresponds to the situation where an objective function reaches a feasible solution on the computing solver, where all the constraints are satisfied with the lower and upper bounds. Figures 35 are plotted with the unit of , 4, and 5 in degree Celsius and with , 2, 6, 7, 8, and 9 in millimeters on -axis. The optimal design variables versus individual and multiobjective optimization of thumb are presented in Figure 3.

Figure 3 helps to interpret the results relating to the optimal design variable, and it is associated with the optimization problem . The -axis shows the numerical values of the optimal design variable , and the abscissa axis is associated with the independent design variable. Table 5 presents the optimal values taken by the objective functions at the pareto-optimal solution for multiobjective optimization.

The other nonoptimal values calculated earlier when obtaining the optimal values of the functions to be optimized are removed from Table 5 for the reasons of clarity. Figure 4 is plotted by using the partial data from Table 4 from the singular optimization.

Figure 4 facilitates the interpretation of the results relating to the optimal design variables relative to individual optimization. In essence, these variables are expressed as a function of the objective functions relating to the various optimization problems to be solved. Figures 35 are associated with a function , the -axis shows the numerical values of the optimal design variable with the unit of , 4, and 5 in degree Celsius and with , 2, 6, 7, 8, and 9 in millimeters as indicated in Table 3. Figure 6 is associated with optimization problem and the abscissa axis is associated with the functional transversal and longitudinal dimensions of the mechanism (). The units of the various objective functions are indicated in Table 5 as follows: and are in degree Celsius, and do not have unit, while is expressed in N.mm and is given in USD. The optimal design variables versus multiobjective optimization of the thumb are presented in Figure 5.

Figure 5 is associated with the objective functions fi optimization problem , the -axis shows the numerical values of the optimal design variable , and the abscissa axis is associated with the variable design parameters.

Figure 5 is plotted based on the partial data from Table 4. The ordinate axis represents the optimal numerical values given by the design variables at the end of the process of simultaneous optimization of a set of coupled objective functions according to the strategy adopted. Moreover, Figure 3 helps to identify the behavior of individual design variables with respect to each simultaneous optimization problem. Figure 4 which presents the behavior of each objective function versus the optimization problem is plotted by using Table 5.

Figure 6 demonstrates a proper understanding of the behavior of all the objective functions relating to the optimization problems during the optimization process of each of the objective functions relating to the different problems to be solved. The ordinate axis gives the values of the function considered, and the abscissa axis shows all the functions for all the resulting optimization strategies. Figure 6 also gives a panoramic idea of the behavior of each objective function according to the optimization problem. The -axis represents the optimal values taken by objective functions when an optimum is obtained at the end of the optimization process of the different optimization problems. Figures 46 show simultaneous optimization of all the criteria confer on the thumb, with the same length of the fixed bar , the same thickness of all the bars, and about the same diameter of the driving pulley. Two observations can be highlighted from the beginning of the manufacturing cost and reliability relating to the optimization process: the first is related to the reliability regarding all criteria being appreciably equal to the unit and the second observation is related to the minimum manufacturing cost, from which the largest design variable values are obtained. The manufacturing constraints such as geometrical and dimensional tolerances, accuracy (0.1 μm), production set (), cutting tool life span, and the shape of the workpiece guide the designer to choose the suitable manufacturing model cost parameters. The optimization process was conducted by varying three times the value of the stopping criteria as indicated at the end of the last paragraph of Section 2, to ensure the accuracy provided by the CNC machining process. The final results displayed at the end of optimization process are the optimal output, respecting all the constraints including manufacturing ones. The optimal objective function values from individual optimization were used as the weighting absolute goal to be achieved by the multiobjective “fgoalattain” solver. The steps where the system stop with the “nonfeasible” answer of the algorithm are cancelled and not validated for this study. The accuracy configures in the computing solver for design variable, objective function, and constraints are indicated in Section 2.2. From there, the precision observed on the results of Table 4 is justified and imposed to the manufacturer the used of computing-aided manufacturing to produce the part of this artificial mechanism.

The value of the optimal function obtained during the entire optimization process relative to the reliability of the whole mechanism is 0.9999, compared to 1 obtained by Ngale et al. [28], when only the driving bar was taken in consideration during the optimization processes. This difference of 0.1111 justifies the contribution of each bar in the reliability mathematical model. Furthermore, the multiobjective optimization showed that the mechanical reliability and the production cost included in any combination of the simultaneous optimization enabled the achievement of the same optimum variables design, though with some exceptions.

Earlier works from literature realized by Ngale et al. in 2003 and 2016, Lim et al. in 2018, and Tsamo et al. in 2020 on the prosthesis mechanism’s optimization show that optimal design variables obtained from one single optimization to multiobjective optimization must be different to justify the contribution of the said performance criteria on the improvement of the design [24, 26, 27]. In this paper, the global manufacturing cost and the strength reliability of the whole mechanism are included in all combinations of the simultaneous optimization necessary to obtain the random optimum design variables. The optimal manufacturing cost generated from the numerical simulation is 501.0021 USD, for the thumb versus 961 USD for the thumb proposed by Ventimiglia [34] and 553 USD proposed by Choi et al., for one compliant finger, obtained by 3D printing [16]. The difference is 459.9979 USD with Ventimiglia price and 51.9979 USD with Choi et al. price. This difference is due to the inclusion of the power system, actuation, and miscellaneous component cost by Ventimiglia. The particularity of the cost generate on this study is that the design of the mechanism is based directly on the objective function modelled. The gap observed at the level of the price can be affected to numerous factors such as human factor error while evaluating empirical quantities, machine imprecision, labor charges, the depreciation rate of the equipment, and productivity assets which are included in the mathematical’s manufacturing cost model of this study [3537]. Otherwise, the simultaneous optimization of all the criteria made it possible to reduce the standard deviation between all design variables. One notices that the virtual mechanisms of the fingers and the thumb are sensitive to each of the new performance criteria such as manufacturing cost and the strength reliability of the whole mechanism. The value of the optimal function obtained during the entire optimization process relative to the reliability is quite different with some exceptions of those obtained by Ngale and Vinet and Ngale who focused their study only on the reliability of the driving bar of the four-crossed bar mechanism [14, 20].

4. Conclusion

The paramount goal of this study was to optimize the thumb prosthesis’s mechanism by taking into consideration the manufacturing cost model based on machining cost theory which took into detail the shape of the workpieces and the strength reliability of all the parts composing the entire mechanism. These criteria were added to those choosing from the state-of-the-art based on their popular contribution relative to the structure’s design optimization. The virtual mechanisms of the thumb are sensitive to each of the new performance criteria employed. The value of the optimal function obtained during the entire optimization process relative to the reliability of the whole mechanism is 0.9999, compared to 1 obtained when only the driving bar was taken in consideration during the optimization processes. This difference of 0.1111 justifies the contribution of each bar in the reliability mathematical model. Furthermore, the multiobjective optimization showed that the mechanical reliability and the production cost included in any combination of the simultaneous optimization enabled the achievement of the same optimum variables design, though with some exceptions. It was found that the mechanical error has little influence on the design variables resulting from the simultaneous optimization. For the simultaneous optimization of all the criteria where it was integrated, its influence was reduced considerably because all the optimal values of the objective functions are lower than the value resulting from the individual optimization. The optimal manufacturing cost generated from the numerical simulation was 501.0021 USD, for the thumb versus 961 USD, and 553 USD for one compliant finger, obtained by 3D printing. The difference is 459.9979 USD and 51.9979 USD more than the price generated on this study. This difference is due to the inclusion of the power system, actuation, and miscellaneous component cost by Ventimiglia. The particularity of the cost generate on this study is that the design of the mechanism is based directly to the objective function modelled. The gap observed at the level of the price can be affected to numerous factors such as human factor error while evaluating empirical quantities, machine imprecision, labor charges, the depreciation rate of the equipment, and productivity assets which are included in the mathematical’s manufacturing cost model of this study. Therefore, further studies in the biomechanics field could integrate the weight and the tribological phenomenon in the optimization processes of the mechanism to achieve the real design of these artificial devices.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.