Trajectory Control in Discrete-Time Nonlinear Coupling Dynamics of a Soft Exo-Digit and a Human Finger Using Input–Output Feedback Linearization

: This paper presents a quasi-static model-based control algorithm for controlling the motion of a soft robotic exo-digit with three independent actuation joints physically interacting with the human ﬁnger. A quasi-static analytical model of physical interaction between the soft exo-digit and a human ﬁnger model was developed. Then, the model was presented as a nonlinear discrete-time multiple-input multiple-output (MIMO) state-space representation for the control system design. Input–output feedback linearization was utilized and a control input was designed to linearize the input–output, where the input is the actuation pressure of an individual soft actuator, and the output is the pose of the human ﬁngertip. The asymptotic stability of the nonlinear discrete-time system for trajectory tracking control is discussed. A soft robotic exoskeleton digit (exo-digit) and a 3D-printed human-ﬁnger model integrated with IMU sensors were used for the experimental test setup. An Arduino-based electro-pneumatic control hardware was developed to control the actuation pressure of the soft exo-digit. The effectiveness of the controller was examined through simulation studies and experimental testing for following different pose trajectories corresponding to the human ﬁnger pose during the activities of daily living. The model-based controller was able to follow the desired trajectories with a very low average root-mean-square error of 2.27 mm in the x-direction, 2.75 mm in the y-direction, and 3.90 degrees in the orientation of the human ﬁnger distal link about the z-axis.


Introduction
Every year, 800,000 people suffer a new stroke in the United States (this number globally is 13,000,000) [1]. Stroke remains the leading cause of long-term disability in the human upper extremities (more than 80%), and only 10% can regain partial mobility [2]. Clinical studies have shown that early access to physical and occupational therapy interventions with intensive/repetitive movements, such as continuous passive motion (CPM), help restore functionalities in the impaired hand of post-stroke individuals [3]. However, the lack of patient compliance and limited availability of resources and therapists hinders the efficacy and effectiveness of these approaches. On the other hand, robotic systems are very effective in performing repetitive and labor-intensive tasks. Several robotic devices have been developed and studied for upper limb rehabilitation [4]. Particularly, wearable robots for upper-body rehabilitation and physical assistance have benefited from intrinsic mechanical compliance inherited in soft robotic approaches. This feature makes them suitable for safe physical interaction with the human body [5,6]. Additionally, soft wearable robots for upper limb rehabilitation address the limitations of their rigid counterparts (i.e., conventional rigid robots), including complex mechanisms, heavy weight, lack of safety, and cost [7][8][9][10][11]. More than 50 soft exoskeletons have been developed over the past 15 years for hand rehabilitation [12,13]; these robots still provide limited capabilities for motion sensors, such as inertial measurement units (IMUs) [10,14], electromagnetic tracking [9], and infrared distance sensors [31]; internal pressure and tactile force sensors; or a combination of these sensors [15] have been used. Though these approaches provide intuitive gross hand and finger motion control, the aforementioned algorithms are mainly for a pure sensor-based position or grip force control, which lacks the required hybrid force-position-based dexterous manipulation to perform fine motor skills. Additionally, extensive reviews of literature on soft robotic hand exoskeletons [52][53][54]63] have shown that current control algorithms do not consider multi-joint control and the multi-contactpoint interaction of the human fingers and the interaction of soft wearable robots, thereby impeding the achievement of complex fine motor functional tasks. Thus, a new control paradigm from a human-robot interaction perspective (i.e., a bilateral approach to facilitate performing collaborative tasks) is needed to take into account the physical interaction between the soft robot and the human hand.
The effective operation of soft robotic hand exoskeletons for physical therapy and assistive motion requires the utilization of a multi-joint control scheme that takes into account the physical interactions between the robot and the human hand [34,64,65]. The highly nonlinear nature of soft robots' dynamics and control, along with the complexity of pHRI and users' biomechanical variation, have impeded the research related to soft wearable robots' control [8,64,65]. Some early works in this field have developed dynamic models and simple control algorithms for soft pneumatic actuators in the form of a continuous fluidic elastomer actuator (FSA) [66,67]. Physical interactions between an FSA and a human finger model were experimentally studied [66] using tactile sensors embedded between the human finger model and the soft robot. In another work, the interactions of a fabric-based pneumatic exoskeleton with a human finger were studied through modeling and experimental validation [35]. None of these works has considered multi-contact-point interaction between the soft robotic digit and the human finger. Our previous work includes extensive studies of design, modeling, and kinematic-based control algorithms on a hybrid soft-and-rigid actuator-based hand rehabilitation robot (REHAB Glove) [10,11,34]; modeling and characterization of physical human-robot interaction (pHRI) in wearable soft hand robots [68,69]; a kinematic-based adaptive control of a bilateral rehabilitation robot [14]; and simulation-based quasi-static force and position control [70]. In our recent work, single-input single-output (SISO) quasi-static model-based adaptive position control of soft exo-digit and the human finger physical interaction was developed and experimentally verified for a step input and trajectory tracking cases [63].
This work presents a multi-joint quasi-static model-based control algorithm for controlling the motion of a soft robotic exoskeleton digit (exo-digit) with three individually controlled joints while physically interacting with the human finger. The contributions of this work are listed as follows, • Analytical modeling formulation of quasi-static physical interaction between the human finger model and a soft robotic exo-digit with individually controlled joints was derived. • A feedback linearization control algorithm was derived for a nonlinear discrete-time multiple-input multiple-output (MIMO) state-space representation. • The feedback linearizability and stability conditions of the nonlinear discrete-time state-space model of soft exo-digit and human finger physical interaction were studied. • Experimental testing of the soft exo-digit interacting with a human finger model for tracking a desired trajectory was carried out.
The structure of this paper is set as follows: The quasi-static analytical model is derived in the Section 2 with the forward and inverse kinematics formulations of the human finger model and the Jacobian of each model. Section 3 presents the nonlinear discrete-time state-space model of a soft robot and the human finger physical interaction. Additionally, a model-based control law was derived based on the input-output feedback linearization for trajectory tracking of the human fingertip pose while assisted by the soft exo-digit. Section 4 describes the implementation of the electro-pneumatic control system hardware for testing the trajectory tracking cases, and finally, Section 5 discusses the results obtained from simulations and experimental testing, which validated the presented model-based controller.

Quasi-Static Model of Physical Human-Soft Robot Interaction
The physical interaction between a human finger and a soft-wearable robot is studied using a model of soft robotic exo-digit and a corresponding human finger model, as shown in Figure 1a. Continuous passive motion (CPM) therapy is the main application of this human-robot interaction. In CPM, the robot actively applies slow, repetitive motions to the user's hand while the user is passive. Thus, the physical interaction can be modeled as a quasi-static process. The soft robotic exo-digit consists of three soft continuous actuating sections (i.e., half-bellow shape hollow structures) representing the joints of the soft-bodied robot and four semi-rigid blocks linked through those pseudo joints. The human finger is modeled as an articulated rigid multi-body with three revolute joints representing the metacarpophalangeal (MCP), proximal interphalangeal (PIP), and distal interphalangeal (DIP) joints of the human finger. Rotational (torsional) springs are integrated at each joint of the human finger model to replicate the biomechanical stiffness of the human joints. The torsional springs are selected from the vendor with a stiffness of k 1h = 5.9 × 10 −4 N·m/deg, k 2h = 3.83 × 10 −4 N·m/deg, and k 3h = 1.07 × 10 −4 N·m/deg for MCP, PIP, and DIP joints, respectively. The semi-rigid links of the soft robotic exo-digit are attached to the links of the human finger model (as shown in Figure 1a), where the physical interaction occurs; see Figure 1b. Frames are attached to each soft-bodied link to describe the kinematic motion of each segment. The quasi-static interaction of the soft robotic digit and the human finger model were decoupled to derive the equations of motion (balance) for each body where the physical interaction (contact forces) is distributed along the length of the two bodies. It is assumed that the reaction force between the human finger and the soft robotic digit happens at the center of each link and is dominant in the direction normal to the soft-bodied links. Figure 2 shows the workspace of the human finger model and the soft exo-digit. In both cases, the joint angles range from 0-90 • for MCP, 0-100 • for PIP, and 0-70 • for DIP, and the link length is measured by the distance from one joint to another. The biomechanical range of motion of a human finger served as the basis for setting these ranges. From Figure 2 it can be seen that the shape of the workspace for both the human finger model and the soft exo-digit is the same.

Kinematics of the Human Finger Model
The human finger was modeled as a serial kinematic chain of three moving rigid links and a fixed-base link with a single-degree-of-freedom hinge joint connecting adjacent links as shown in Figure 1. The motion of the joints is described by q H = [q H,1 , q H,2 , q H,3 ] T . The forward kinematics of the human finger model, L(q H ), can be described similar to a classical planar three-link robotic arm where the pose of the tip of the human finger, (x, y, θ), with respect to the reference frame is given as follows, where the following notations are used in the entire paper for sine and cosine functions, s i = sin(q i ), c i = cos(q i ), s ij = sin(q i + q j ), c ij = cos(q i + q j ), s ijk = sin(q i + q j + q k ), and c ijk = cos(q i + q j + q k ). Additionally, the Jacobian matrices, J H,i , which map the joint velocitiesq H to the linear υ H,i and angular ω H,i velocities of each link can be obtained aŝ whereυ H,i is a vector that combines the planar linear and angular velocity components, i is the length of ith link. Note these Jacobian matrices are derived for the center point of each link where the interaction force applies.

Coupled Human-Robot Interaction Quasi-Static Model
The multi-contact points quasi-static models of the human finger and the soft robotic exo-digit, developed separately, are presented here. For simplifying the form of equations, the contact forces, expressed in the reference frame, are written in terms of body-fixed frames as follows: which combines the bodyfixed frame forces, the quasi-static equations can be written in the following compact form, where the quasi-static equations for the human finger model are defined as follows: where the special Jacobian matrix (J * H ) T and the joint stiffness matrix K H are defined as follows: Similarly, the quasi-static equations for the soft robotic exo-digit are defined as and p = [p 1 , p 2 , p 3 ] T is the actuation pressure vector (control inputs), and B R is a diagonal matrix of the first moment of the cross-section area of each soft segment where the internal pressure actuation p i acts on them, It is obvious that both upper triangular special Jacobian matrices (J * H ) T and (J * R ) T are invertible by assuming q = 0. Now, Equations (8) and (9) are combined to derive a collective model for pHRI. The combined model will be used for the control algorithm development. Eliminating the force vectors, f * RH and f * HR between Equations (8) and (9) which represents the coupled quasi-statics of the human-robot model. The torque term is given as The coefficients of the polynomials of the rational functions in Equation (11) are given in the Appendix file of the Supplementary Materials folder. Note that these coefficients are functions of material properties (i.e., the Yeah third-order hyperelastic model [11]) and the geometry of the soft actuator segments).

Nonlinear Discrete-Time State-Space Representation
Control problem: a square multiple-input multiple-output (MIMO) tracking control problem is considered here such that the pose of the human fingertip (i.e., (x, y, θ)) follows the reference inputs using the soft robotic exo-digit. For this purpose, the coupled quasistatic equation, shown in Equation (10), can be written into a nonlinear discrete-time state-space form with affine control as follows for position control, where x k = q(k) is the state vector at time index k, p j k is the jth control input at time index k in p k = [p 1 k , p 2 k , p 3 k ] T , and the targeted output vector y k = h(x k ) is defined as the position vector of the human fingertip [x(k), y(k)] T = L(x k ), given in Equations (1) and (2), and the bending of the human finger distal end (i.e., θ(k) = x 1 k + x 2 k + x 3 k ). f(x k ) and g j (x k ) are given as follows: and We assume that there is an equilibrium point (x e k , p e k ) that satisfies x k+1 = 0 in Equation (12).

Input-Output Feedback Linearization
Using the notion of input-output linearization [71,72], the discrete derivative of the output y k in Equation (12) yields where are Lie derivatives of h(x k ) with respect to the field functions f(x k ) and g j (x k ), respectively. A remarkable result from the work by Jayaraman and Chizeck [73] has shown that a nonlinear discrete-time system like Equation (12) with dimension n is feedback linearizable if it has vector relative degree {r 1 , . . . , r m }, where m is the dimension of the control input p k , such that r 1 + r 2 + · · · + r m = n. Note that the MIMO system has a uniform relative degree of one given that L g j h(x k ) = 0 in Equation (12) for j = 1, 2, 3 in the neighborhood of x e k . The details of analysis on the value of L g j h(x k ) and its variation with respect to the configuration space (q 1 , q 2 , q 3 ) of the soft exo-digit are given in Appendix A.
The goal is to find a coordinate transformation, z k = T(x k ) and a feedback control p k = γ(x k , υ k ) to linearize Equation (12) for the mapping between the inputs and outputs for the trajectory tracking problem. Thus, z k = y k − y d k , i.e., the trajectory tracking error, is selected where y d k is the vector of the reference inputs (i.e., the desired outputs) which the human finger controlled by the soft exo-digit must follow. Given the discrete derivative of z k , and substituting Equation (15) into Equation (16) yields the new state space model in the new coordinates expressed by a set of scalar equations for i = 1, 2, 3 given as follows: exists, then the control input vector (p k ) (i.e., the actuation pressure of the soft segments) can be derived as a nonlinear state feedback control [73], to linearize the input-output map which transforms Equation (12) to its equivalent linear system as follows: Invertibility of Equation (18) is shown in Appendix A. The additional control input υ k will be defined by a linear full-state feedback control as follows: where K υ ∈ R 3×3 is a diagonal feedback gains matrix. Choosing K υ to be positivedefinite guarantees exponentially asymptotic stability of the input-output linear system, Equation (20), and consequently, guarantees the stability of the original nonlinear statespace system, Equation (12), to follow the reference signal y d k . The control block diagram is shown in Figure 3.

Control System Implementation
An electro-pneumatic control system was developed to examine the model-based control algorithms derived in Section 3. The control system is based on modulating the actuation pressure, where the control action will be determined based on the error between the desired trajectories (pose of the finger distal end) and the actual feedback (i.e., determined through IMU sensors reading and the forward kinematics of the human finger) as shown in Figure  , and the current internal pressure (p c ) measured by in-line pressure sensors, p = p c + ∆p, where the pressure change (∆p) is associated with the tracking error and quasi-static model of the human-robot interaction and will be determined by Equation (19). A PWM control scheme was used, as shown in Figure 4b, through an Arduino to control the opening and closing cycles of three-way solenoid valves proportional to the input control value for achieving the pressure change. A two-way valve is integrated into the exhaust port of each three-way valve to control the pressure relief process. A PWM signal was sent to the two-way valves for closing and opening. Note that the positive change in pressure leads to pressurizing the soft actuator while the negative change in pressure causes vacuuming (or air bleeding). Each soft actuator segment (soft exo-digit module) is individually controlled, while their collective controlled actuation leads to tracking the desired movement of the human finger by the soft exo-digit. Table 1 provides the description and justification of each component used in the experiment. Figure 4c shows the actual electro-pneumatic board, the soft-robotic exo-digit, and the human finger model experimental setup with each component labeled according to the corresponding number in Table 1.  Table 1.

Analysis of the Control Law
The control law obtained for the actuation pressure of the soft robotic exo-digit expressed in Equation (19) was examined for the configuration space of the robot where y d k+1 = 0. This examination provides insights into how the controlled pressure of each soft actuator segment varies with respect to the variation of joint angles q = [q 1 , q 2 , q 3 ] T . Note that the actual movement of the soft robotic exo-digit along with the human finger (i.e., particularly the pose of the tip of the finger) in the Cartesian space will be corresponding to a subset of this full configuration space of the soft robot. Understanding the overall behavior of required control pressure helps to better determine the constraints on the physical human-robot interaction and the design of the control system for experimental testing. As mentioned before, the range of each joint angle is set based on the biomechanical range of motion of the human finger (MCP: 0-90 • , PIP: 0-100 • , and DIP: 0-70 • ). The plots are provided for each actuation pressure in Figures 5-7. The color variation in these plots corresponds to changes in pressure, with blue indicating the lowest pressure and yellow indicating the highest. The exact pressure values can be determined from the vertical line on the right side of each figure. It is shown that the pressure actuation for the first joint varies with respect to the change of all three joint angles, although this behavior is different for the other two joints where they are independent of the variation of q 1 (see Figures 6 and 7). Additionally, for the third joint actuation pressure p 3 k , it seems its value does not depend on the variation of the second joint angle q 3 . These trends are consistent with the coupling in the quasi-static model of the soft robotic exo-digit and the human finger interaction, where most of the matrices are upper-triangles that make the third joint to be decoupled from the other two joints, and consequently, the second joint is only coupled with the third joint. In contrast, the first joint is coupled with both the second and third joints as it is expected. Additionally, Figure 5 shows that the actuation pressure for the first joint becomes constant (i.e., it is independent of the variation of the other two joints) after a certain value of the second joint (q 2 ). The threshold for the start of the constant pressure region decreases (i.e., it occurs at a lower q 2 ) with an increase in the first joint angle (q 1 ). The reason for this behavior could be due to the counter-effect of the second and third joint motions and their coupled quasi-static torque reactions after reaching certain joint angles that decouples the effect of the last two joint motion from the actuation pressure of the first soft actuator.

Trajectory Tracking of the Desired Fingertip Pose in Simulation
In this part, we theoretically examined the performance of the controller by following the desired fingertip pose for two different trajectories. Figure 8a shows the desired path (x d and y d ) of the fingertip in Cartesian space, and the corresponding joint angles (motions) are depicted in Figure 8b. The path is arbitrary, but we considered a smooth curve function. Based on the latter, we can determine the desired orientation of the fingertip (i.e., θ d = q 1 + q 2 + q 3 ). The desired continuous trajectories (x d , y d , and θ d ) are discretized into a series of discrete points, and the collective set of these values was used in Equation (19) to calculate the required actuation pressure at each soft actuator joint in order for the soft robotic exo-digit to move the human finger along the desired trajectories. The results are shown in Figure 9 for the varying values of the three joint angles. The color of each point in the plots represents the values of the pressure magnitude according to the provided color map. The results show that the control input solutions exist for tracking these specific trajectories in the configuration space of the soft robotic exo-digit with a trend of increasing pressure. The required theoretical pressure inputs are in the range of 0-30 kPa, 0-35kPa, and 0-80 kPa for MCP, PIP, and DIP joints, respectively, which are achievable by the control hardware system (Figure 4c).
In the second case, we used the experimental data of an actual human finger movement (flexion) to generate the desired trajectories, as shown in Figure 10a. The corresponding joint angles are determined based on the inverse kinematics of the human finger model, which will be used to calculate the desired orientation of the fingertip (i.e., θ d = q 1 + q 2 + q 3 ) for this case. Similar to the first case, these continuous functions were discretized into a series of way-points where the actuation pressure of each soft actuator was determined, as shown in Figure 11. The result shows the existence of the control solutions for following the second desired trajectory, while the magnitude of the required actuation pressure can be obtained from the provided color map. Note that the magnitude of the actuation pressures is higher compared to the first case as the desired trajectory requires more joint variation close to the full range of human finger range of motion. The maximum pressures (160 kPa) are beyond the range of our control hardware. Additionally, the rate of change of joint motion (joint velocities) is selected to be slow in order to be consistent with the quasi-static motion assumption. This is required by the therapeutic interventions for stroke patients as the fast movement of their fingers leads to spasticity, a condition in which there is an abnormal increase in muscle tone or stiffness of muscle [14].

Trajectory Tracking of the Desired Fingertip Pose
A desired trajectory of a human fingertip is provided in the form of equations x = 0.002002t 3 − 0.1278t 2 − 1.848t + 99.72, y = 0.0002511x 3 − 0.0407x 2 + 2.728x − 150.5, θ = −4t − 36. With the given pose (x, y, θ) of the fingertip, where each parameter varies with time, t, the three joint angles associated with the finger are calculated using inverse kinematics. The experimental setup includes an Arduino Mega 2560, four IMUs MPU6050 to measure the angular position of each joint in a human finger model, three MPRLS pressure sensors to measure the pressure buildup in three joints, three FA0520E three-way air valves to increase the pressure in each joint, three FA0520D two-way valves to release pressure from the joints, and six HM MOS modules for opening and closing the valves. An I2C Multiplexer (TCA9548A) is used to connect the pressure sensors and IMUs to the Arduino. The system is controlled by an Arduino program that adjusts the pressure applied to the joints to match the desired joint angles. The code starts by defining the output pins of the valves, reading the IMUs to obtain the joint angles, calculating the desired joint angles using inverse kinematics and comparing the actual and desired angles. If the desired angle is greater than the actual angle, the three-way valve of the joint opens to increase the pressure. If the desired angle is less than the actual angle, the two-way valve releases air from the joint, causing the red point in the fingertip, as shown in Figure 12, to move along the desired trajectory. Figure 13 depicts the joint angles for three joints (MCP, PIP, and DIP) and compares the actual angles with the desired angles for experiment-1. The outcomes reveal that the finger was able to accurately trace the desired joint angle for all three joints, with a root-mean-square error (RMSE) of 1.32, 1.46, and 0.86 degrees for each joint, respectively. The actual joint motions show an oscillating behavior around the desired joint trajectories. This behavior is due to the use of solenoid valves and a pulse width modulation (pwm) scheme for modulating the pressure of the soft actuators based on the presented control law. Replacing the solenoid valves with proportional valves, which operate with continuous opening and closing rather than pulsing, could help eliminate the chattering responses of the soft actuators. Additionally, the range of motion of the trajectory tracking can be increased by using a pressure sensor that measures in a higher range compared to the current sensors; the current pressure sensors in the control system measure in about the 0-70 kPa range, which led to a limited range of motion available for each soft actuator. Figure 14 demonstrates the expected and actual path of the fingertip as measured by both the IMU-based tracking and the camera. The IMU-based tracking uses the measured joint angles of the human finger and the forward kinematics equations to determine the fingertip pose. The picture confirms that the fingertip followed the desired path closely, with an RMSE of 7.49 mm in the x direction and 12.35 mm in the y direction for the camera measurement, 1.98 mm in the x direction, and 2.79 mm in the y direction for the IMU measurement. The deviation between the camera and IMU-based tracking is mainly due to the error of the camera software to accurately and consistently detect the exact location of the center of the marker attached to the fingertip.
The experiment was repeated for the second time, and the results for the joint angle and the tip trajectory are shown in Figure 15 and Figure 16, respectively. The obtained results confirm the repeatability of the control algorithm and the hardware. Once again, the actual joint angles followed the desired joint angles with a low RMSE of 1.63, 1.60, and 0.92 degrees for each joint, respectively. The RMSE for the camera measurement was 6.77 mm in the x direction and 11.15 mm in the y direction, while for the IMU measurement, it was 2.55 mm in the x direction and 2.71 mm in the y direction. Note that the resulting RMSE for IMU-based tracking and the camera-based tracking are very close and consistent between the two cases.
The Supplementary Materials folder contains the videos for both experiment-1 and experiment-2.

Conclusions
An algorithm for regulating the motion of a soft robotic exo-digit with three separate actuation joints that physically interact with a human finger is presented in this research using a quasi-static model-based approach. The motion of the soft exo-digit was controlled by regulating the position of the tip of the human finger model. The physical interaction between the soft exo-digit and a model of a human finger was developed using a quasi-static analytical model. The model was provided as a nonlinear discrete-time MIMO state-space representation for the control system design. The actuation pressure of a single soft actuator was used as the input, and a control input was created to linearize the input-output, with the output being the position of a human fingertip. The nonlinear discrete-time system for trajectory tracking control's asymptotic stability is examined. The experimental test setup included a soft robotic exo-digit and a 3D-printed model of a human finger that was integrated with IMU sensors. For regulating the actuation pressure of the soft exo-digit, Arduino-based electro-pneumatic control hardware was developed. The controller's performance was evaluated through simulated research and experimental testing for following various position trajectories similar to the human finger pose throughout daily tasks. The experiment was repeated twice, verifying the control algorithm's repeatability. With a low average root-mean-square error of 2.27 mm in the x-direction, 2.75 mm in the y-direction, and 3.90 degrees in the orientation of the human finger distal link about the z-axis, the model-based controller was able to follow the desired trajectories.
In this study, we considered a model of human-robot interaction with simplifying assumptions such as linear torsional springs with a constant stiffness at the joint of the human finger model and the soft exo-digit's soft actuation segments with a constant length and curvature. Moreover, model parameters uncertainties involved in both the soft robotic exo-digit and human finger model as well as variation of biomechanical parameters of the human finger model (i.e., joint stiffness and damping, rigid links' inertia, and overall size), have not been considered. Apart from that, the human finger model was passive, so it only contributed a resistive force to the motions of the soft exo-digit. In contrast, the active participation of the human hand in daily activities is inevitable. These factors are ascribed to the study's limitations. Considering all these limitations, future research will adopt a more realistic model of pHRI, develop adaptive robust hybrid force-position control laws to accommodate multi-task aspects (i.e., required simultaneous force and position controls) of the human hand fine motor motion, adaptability to variations in the biomechanics of the users, and robustness against uncertainties of the models.  Data Availability Statement: Access to data will be provided upon request.

Acknowledgments:
The authors would like to thank Joshua Kirkland for helping with 3D-printing parts for the experimental test setup used in this study.

Conflicts of Interest:
The authors declare no conflicts of interest.
L g 3 h(x k ) = A p3 r p3 x 2 3 k 3 k 2 x 2 σ 9 σ 8 σ 5   −l 2 3 σ 11 k 2 σ 8 σ 5 − (σ 1 + l 3 σ 11 )σ 4 k 3 σ 5 − (σ 1 + 2l 1 s 1 + l 3 σ 11 )σ 3 k 3 k 2 σ 9 σ 8 l 2 3 σ 10 k 2 σ 8 σ 5 x 2 + (σ 2 + l 3 σ 10 )σ 4 k 3 σ 5 x 2 + (σ 2 + 2l 1 c 1 + l 3 σ 10 )σ 3 k 3 k 2 x 2 σ 9 σ 8 l 3 k 2 x 2 σ 8 σ 5 + σ 4 k 3 σ 5 + σ 3 k 3 k 2 x 2 σ 9 σ 8   (A3) where σ 1 ,...,σ 8 are given in the appendix file of the Supplementary Materials folder. The determinant of L G h(x k ) is provided here, and the values of the determinant are obtained for the configuration space of the soft robotic exo-digit as illustrated in Figure A1 where the pseudo joint bending angles varied in the range of motion of the human finger joints. The results are presented for varying values of q 1 , while the surface demonstrates the variation of the determinant for variation q 2 and q 3 . Excluding the q 2 = 0, which is the singularity of L G h(x k ), the determinant stayed finite and there is a trend of increasing values towards the upper limit of q 2 and q 3 for the range of 0 < q 1 ≤ 50 • . As shown in Figure A1, after q 1 larger than 60 • , there is a discontinuity in the determinant manifold towards the upper limit q 2 and q 3 . The discontinuity is displayed on the 2D projection of the manifold in Figure A2. Despite the discontinuity, 4A p1 A p2 A p3 r p1 r p2 r p3 l 2 1 l 2 2 l 3 x 4 1 x 2 x 3 3 s 1 k 1 k 2 k 3 σ 4 σ 5 x 2 1 (x 3 c 1l2 σ 9 + x 2l3 σ 7 ) + x 2 x 3 (l 1 σ 6 + x 1 (l 1 σ 8 + x 1 σ 3 )) where σ 1 ,...,σ 9 are given in the appendix file of the Supplementary Materials folder. Figure A1. Variation of L G h(x k ) of with respect to q 2 and q 3 for different values of q 1 . Figure A2. 2D plot of surface plots in Figure A1, the red arrow shows the discontinuity in L G h(x k ) for that region of the configuration space.