Next Article in Journal
Equations of Disturbed Motion of the Moving Part of the Gyroscope Suspension
Next Article in Special Issue
Trunk Posture from Randomly Oriented Accelerometers
Previous Article in Journal
Detecting Pest-Infested Forest Damage through Multispectral Satellite Imagery and Improved UNet++
Previous Article in Special Issue
Accounting for Viscoelasticity When Interpreting Nano-Composite High-Deflection Strain Gauges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Planar Strain Sensor Design for Capturing 3-Dimensional Fingertip Forces from Patients Affected by Hand Paralysis

1
Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
2
Department of Neurology, Johns Hopkins Medicine, Baltimore, MD 21287, USA
3
Department of Kinesiology, University of Georgia, Athens, GA 30602, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2022, 22(19), 7441; https://doi.org/10.3390/s22197441
Submission received: 23 August 2022 / Revised: 21 September 2022 / Accepted: 27 September 2022 / Published: 30 September 2022
(This article belongs to the Special Issue Biomedical Sensing for Human Motion Monitoring)

Abstract

:
Assessment and therapy for individuals who have hand paresis requires force sensing approaches that can measure a wide range of finger forces in multiple dimensions. Here we present a novel strain-gauge force sensor with 3 degrees of freedom (DOF) designed for use in a hand assessment and rehabilitation device. The sensor features a fiberglass printed circuit board substrate to which eight strain gauges are bonded. All circuity for the sensor is routed directly through the board, which is secured to a larger rehabilitative device via an aluminum frame. After design, the sensing package was characterized for weight, capacity, and resolution requirements. Furthermore, a test sensor was calibrated in a three-axis configuration and validated in the larger spherical workspace to understand how accurate and precise the sensor is, while the sensor has slight shortcomings with validation error, it does satisfy the precision, calibration accuracy, and fine sensing requirements in orthogonal loading, and all structural specifications are met. The sensor is therefore a great candidate for sensing technology in rehabilitation devices that assess dexterity in patients with impaired hand function.

1. Introduction

Neurological conditions such as stroke and spinal cord injury can cause significant hand impairment—both in terms of dexterity and strength [1,2,3]. Therefore, there is a need to accurately quantify these two components of hand control.
It has been demonstrated that increasing instructed force in a finger correlates with increasing involuntary force in non-instructed fingers [4]—a phenomenon called enslaving. The ability to suppress enslaving is known as finger individuation [4]. By quantifying the forces exerted by the instructed and non-instructed fingers, an index for assessing finger individuation independently from strength can be calculated. Additionally, strength and precision control of the fingers have been shown to recover via partially separable control systems [4]. This potentially allows for finger dexterity training even in the presence of significant weakness.
Finger individuation is usually assessed in the context of finger flexion [4]. Knowledge of finger individuation in three dimensions is desirable given that dexterity is characterized by flexion/extension/abduction/adduction movements of the finger. Previously developed research devices such as the device by Diedrichsen and colleagues [4] and the HandCARE Device [5] offer force sensing capabilities for the hand, but still cannot capture the rich 3-dimensional workspace of the hand. Devices with multi-DOF sensing such as the Gifu Exoskeleton [6] incorporate more degrees of freedom per finger, but opt to sense finger torque at each joint rather than fingertip force. Thus, fingertip forces can only be estimated, which can result in inaccuracies in kinematic analyses. Additionally, the transmission forces between joints have to be accounted for in a good kinematic model of the hand or else further inaccuracies can develop. Therefore, in order to accurately measure three-dimensional finger individuation, a novel force sensor is needed that is capable of measuring three-dimensional forces at the fingertip regardless of patient strength. To detect very subtle movements in near-plegic hands, this sensor design should also be very sensitive to small forces.
High quality sensing of 3D forces at the fingertip requires fine resolution, low hysteresis, high precision, high linearity, resistance to fatigue effects, low package mass, and low cost. There are plenty of research devices that capture fine forces, but they only detect 1 DOF [7,8,9,10,11,12] and/or do not report detailed sensor performance [13,14,15,16,17,18,19]. There have been several attempts in the academic literature to capture fine 3D force signals from the fingertips with clear performance measures. Artificial fingers have been embedded with small force-sensing packages based on piezo-resistive MEMS technology [20]. Capacitive surfaces have also been attached to the pads of human fingertips to capture pressure and shear contact forces [21]. Additionally, 3D contact forces can be estimated from a photoplethysmograph of blood volume underneath the fingernail [22]. Despite 3-DOF capability and compact sizing, these sensors have limited directional range and resolve down to only hundredths of a Newton. The most suitable sensor used commercially and in academia is the ATI Nano17, which can capture the 3D force at milli-Newton resolution [23]. However, current market price for a single sensor is over $5000, making it fiscally challenging for a five-finger, low-cost home rehabilitation device. Identifying a current gap of implementation in the hand rehabilitation space, we propose a low-cost sensor technology that captures milli-Newton-level forces in three axes (six-directions) with high capacity for strength characterization, and accuracy suitable for integration into rehabilitation therapy devices. Our sensor is based on a strain-gauge approach, given the technology’s prior success in measuring insect and vertebrate biomechanics [24,25] and surgical instrument forces [26]. In the following sections, we describe the design and characterization of this sensor.

2. Materials and Methods

2.1. Design Requirements

A selective summary of sensors used in upper extremity rehabilitation is detailed in Table 1. In most cases, the sensing technology implemented involved strain gauges or piezoresitive elements.
As described in Table 1, devices in academia and on the market have established sensing capabilities, which we use to inform the specification criteria of our novel sensor. The criteria are as follows:
  • The sensor must fit the fixture geometry indicated in Figure 1. This geometry has a circular area that approximates the cross-sectional area of a fingertip, which is designed to avoid structural collisions with other adjustable finger interfaces (for a multi-finger interface).
  • The sensor, in combination with appropriate amplification and digitization, must sense a resolution of 1 mN or better. For finger force control tasks, oscillatory hand tremor at 10 Hz exhibits typical force power of approximately 0.001 N [29], so this resolution is needed to capture sufficient fine kinematics.
  • The sensor should weigh at most 10 grams so that five would not weigh down their portable chassis in a mutli-finger setup. This is also consistent with the ATI Nano17.
  • The sensor should sense a full-load capacity of 20 N or greater in all spherical directions. Previous literature has determined that fingertip forces can reach up to 60 N for the fingers and up to 110 N for the thumb in dedicated 1D pinch tasks [30], but general 3D finger tasks in a palm-down hand posture report an average finger strength of 20.35 N [4].
  • The maximum accuracy, hysteresis, repeatability, and nonlinearity errors should be at most 5% of the full-scale output. This strikes an appropriate tradeoff from the maximum 3% errors achieved by 1-DOF sensors and is reasonable given that most force sensors for fingertips do not specify their accuracy errors.
  • The sensor should have a maximum displacement of 5 mm at the finger interface, which approximates isometric force sensing.
  • The cost of the sensor from fabrication to assembly should be at most $1000 after economics of scale.

2.2. Design Overview

For sensing forces, a planar strain gauge based architecture was opted for due to the technology’s low cost, ease of inspection, and ease of measurement.
FR-4 fiberglass resin was selected as the planar substrate for the sensor board. This substrate choice allows for simplified routing of shorter and smaller circuit traces, as compared to a sheet metal surface, where thicker insulation for wiring is needed. FR-4 also supports power and ground interior layers, making the power layers accessible anywhere on the board, while FR-4 has a much smaller tensile strength and is more expensive than spring metal [31,32], the simplified wiring and resulting smaller assembly volume compared to other implementations outweighs the material’s limitations.
To generate bending loads, the sensor board had to be fixed to mechanical ground. To properly mechanically ground the board, the aluminum fixture, shown in Figure 2, was machined. Holes were drilled in the sensor board to accommodate fastening to the fixture. To transmit the force from a user’s finger to the sensor board, an aluminum frame with a finger-cup interface was secured to the top of the sensor board, while an identical aluminum frame was secured to the bottom of the sensor board (see Figure 3). The finger cup interface consists of a silicone cup attached to a plastic base. A central hole and alignment holes were used to precisely align the frame to the sensor board (using aluminum dowel pins).

2.3. Force Sensing Principle

To reduce assembly complexity, all sensing is done on the surfaces of a disk of FR-4 material, which for the purposes of numerical analysis can be modeled as a plate. According to the theory of bending plates, induced bending of a thin beam can be approximated as deforming a 2D plate at the beam’s midplane. The midplane is orthogonal to the plate’s thickness, and the deformation acts parallel to the thickness dimension. The resulting deformation generates strain experienced throughout the plate [33]. Essentially if the strain at select locations on the surface are identified, forces from the fingertip, based on angular vector of deflection, can then be back-derived in 3D space. Strain gauges placed at these unique locations identify the relevant strains by stretching or contracting with the deformation, changing the effective electrical resistance of the gauge. The voltages probed from these resistances can then be mapped to corresponding forces. Since strain gauges work under bending rather than shearing forces, any force component parallel to the plate will not induce strain on the gauges. Therefore, all forces from the fingertip should be applied at a consistent location above the plate via the finger-cup interface. With this design, any parallel force component will cause a bending moment about the center of the sensor, which the strain gauges can detect. To capture three degrees of freedom, at least three strain gauges of unique location and orientation are needed. A fourth gauge was also added to provide information about error by introducing redundancy.
At each planar placement on the disk surfaces, a pair of strain gauges are surface-bonded, one gauge on the top of the sensor board and one gauge underneath, as shown in Figure 4. By connecting the top and bottom strain gauges into a voltage divider circuit as illustrated in Figure 5, the output voltage between the gauges can be determined from a coupled strain-dependent pair of resistance changes. This coupled pair of gauges improves the measurement sensitivity, and increases the resolution of any signal acquired from the sensor. The equation for the voltage divider output is expressed as
V d i v = R G , g n d + Δ R g n d R G , + + Δ R + + R G , + Δ R
where R G , + , R G , , and R G , g n d are the gauge resistances of the top surface gauge, the bottom gauge, and the gauge closest to electrical ground, respectively; Δ R + , Δ R , and Δ R g n d are the resistance changes of the top, bottom, and ground-adjacent gauge, respectively, from applied strain; and V d i v is the output voltage of the voltage divider.
To prevent voltage comparisons from canceling out during computation of relative Z-axis signal later, gauge locations 2 and 4 had voltage dividers with excitation on top and ground on bottom shown in Figure 5, whereas locations 1 and 3 had voltage dividers with excitation on bottom and ground on top in a flipped configuration. As a result, the former configuration has the bottom gauge designated as a ground-adjacent gauge, while the latter has the top gauge as ground-adjacent.
Any change in resistance Δ R caused by an external strain on a single gauge is calculated as
Δ R = G F × ε S × R G ,
where G F is the gauge factor, ε S is the strain on surface S, and R G is the nominal resistance of the gauge package. Therefore, Equation (1) expands to
V d i v = R G , g n d + ( R G , g n d × G F g n d × ε g n d ) R G , + + ( R G , + × G F + × ε + ) + R G , + ( R G , × G F × ε )
where G F + , G F , and G F g n d are the gauge factors of the top, bottom, and ground-adjacent gauges; and ε + , ε , and ε g n d are the strains applied at the top, bottom, and ground-adjacent gauges, respectively. Any strain gauge on the top surface, ε + , is assumed to experience equal but opposite strain to the bottom surface, ε , such that ε + = ε . Additionally, all gauges are assumed to have the same gauge factor and nominal resistance. Therefore, Equation (3) simplifies to
V d i v = 1 + G F g n d × ε g n d 2
There are four planar placements where pairs of strain gauges are bonded, specified in Figure 6. Voltage data is interpreted by making comparisons between two different locations. Out of the six unique comparisons, only four are needed to get full voltage information. For example, the comparison signal Δ V 1 2 between two given numerical locations 1 and 2, is computed as:
Δ V 1 2 = V e x c × G F 2 Δ ε 1 2 × G + V e x c 2 ,
where V e x c is the excitation voltage, G F is gauge factor, G is any gain applied to magnify the difference in divider outputs, and Δ ε 1 2 = ε 1 ε 2 is the difference in strain between two locations. The comparison signal voltage should be half the excitation voltage at equal (or zero) strains to minimize saturation of the comparison signal under load. Again, the equations assumes a constant nominal gauge factor, nominal zero-strain resistances across strain gauges, and no thickness change upon deflection.
These comparisons compensate for changes resulting from temperature effects and perform common-mode rejection to eliminate unnecessary direct-current (DC) bias. Therefore, a virtual Wheatstone bridge is created when comparing different placements. The relative voltage signals for each force direction can then be described by the following relationships:
V x = Δ V 4 2 cos ϕ Δ V 3 1 sin ϕ V e x c 2 ( cos ϕ sin ϕ ) V y = Δ V 4 2 sin ϕ + Δ V 3 1 cos ϕ V e x c 2 ( sin ϕ + cos ϕ ) V z = Δ V 3 4 + Δ V 1 2 V e x c e = Δ V 4 2 Δ V 3 1 + Δ V 3 4 Δ V 1 2
where ϕ is the rotational offset of the quadrilateral pattern of gauges from the -Y axis going counter-clockwise on top surface (or clockwise on bottom surface), V e x c is the excitation voltage, and e is the signed gauge board error. Ideally, gauge board error should always be equal to zero from redundancy, and any non-zero value indicates a board defect. In this application, ϕ is 45°. Therefore, the relative signals in Equation (6) simplify to
V x = Δ V 4 2 Δ V 3 1 2 V y = Δ V 4 2 + Δ V 3 1 V e x c 2 V z = Δ V 3 4 + Δ V 1 2 V e x c e = Δ V 4 2 Δ V 3 1 + Δ V 3 4 Δ V 1 2

2.4. FEA Simulation

To model how strain propagates on the sensor, a static bending finite element analysis (FEA) was conducted using the Simulation package in Solidworks 2020 (Dassault Systèmes; Vélizy-Villacoublay, France). For the bending plate simulation, it is assumed that:
  • The strain gauges and substrate operate within their elastic regions.
  • The gauges are symmetrically and perfectly bonded.
  • The axes of bending deformation align with the central point of the disk.
  • The plate is transversely isotropic, which can later be verified by entries of the material’s stiffness matrix and the calibration matrix. The glass fibers that are cross-woven in FR-4 have dimensions on the same order of magnitude, so the bulk material model is approximately isotropic in any planar direction [34].
To break the simulation problem into simpler sub-problems, a representative sub-unit of the sensor (roughly a fourth of the total disk) was modeled that contained a strain gauge, a pair of fastening holes, and a pair of alignment holes. FEA simulation results are shown in Figure 7. A 2.5 N vertical force (directed down into the plate) was chosen for simulation loading so that this loading scenario could be easily replicated during physical characterization; this force amount is a fourth of 10 N, which is expected to be a typical magnitude of exerted force. The strain field shows compression within the strain gauge area when the vertical force downwards is applied. This compression shortens the gauge length, changing its resistance. Therefore, our FEA modeling is reasonable for strain gauge simulation.
After verifying the sub-unit model, the full disk that approximates four sub-units in a diamond configuration was modeled. Geometry of the disk is shown in Figure 8. Boundary conditions were set by the circumference of the disk as well as the interior holes. Here, 10.0 N was applied down the Z-axis so that each sub-unit would still experience 2.5 N of downwards force. An illustration of reaction forces, applied forces, and strain gauge directions is shown in Figure 9. The resulting Von Mises stress field is shown in Figure 10a, with FEA simulation results of strain summarized in Figure 10b,c for each surface coordinate defined for strain analysis, respectively. Since the analysis force is parallel to the Z-axis, the vertical distance of the finger generating the force from the board does not matter. However, simulating sideways forces from the finger through the interface to the board does need careful consideration of the vertical distance since the interface would be imparting a bending moment on the board. For simulations of X-Y forces that act parallel to the board’s X-Y plane, a moment arm of 10.0 N is applied 15.5 mm from the sensor plane. These simulation results are in Appendix A.

2.5. Data Acquisition

Four unique voltage pairs (described in Section 2.3) are each compared and amplified by a corresponding programmable amplifier (Texas Instruments PGA309). For four total comparisons on a single sensor, four amplifiers are needed. The amplifier model has a 100 Hz low-pass filter to remove high frequency noise [35]. The amplifier transforms an input voltage difference into a single output voltage according to the transfer function
V o u t = ( ( Δ V i n + V c ) G i + V f ) G f × G o
where V o u t is the amplified output voltage, Δ V i n is the difference between probing voltages of two voltage dividers, V c is the coarse offset voltage, V f is the fine offset voltage, G i is an input gain, G f is a fine gain, and G o is an output gain. The corresponding inverse function is
Δ V i n = V o u t G f × G o V f G i V c
Table 2 shows the default settings for all parameters. V c is set equal to the zero-load V i n via the standard zeroing procedure detailed in Algorithm 1, and V f is set such that a zero-load causes an output voltage of V e x c / 2 .
Algorithm 1:Procedure for tareing sensor readings via writing amplifier settings
  • Set tolerances 1,2
  • for each finger do
  •     for each amplifier do
  •         Set gains, offsets to default                                                                          ▹ See Table 2
  •         Set counters 1,2 to zero
  •         while counter 1 is less than 10 do
  •            Obtain 1000 readings
  •            Compute average output reading
  •            Compute input estimate from average output with Equation (9)
  •            Round input estimate to nearest discrete coarse offset ▹ Refer to PGA309 User’s Manual Table 6–10 [36]
  •            Set coarse offset to rounded estimate
  •            Write coarse offset to amplifier
  •            Compute error between offset and original input estimate
  •            if tolerance 1 exceeds error then
  •                Break while loop
  •            else
  •                Increment counter 1
  •            end if
  •         end while
  •         while counter 2 is less than 10 do
  •            Obtain 1000 readings
  •            Compute average output reading
  •            Compute difference between half scale output and average reading
  •            if tolerance 2 exceeds difference then
  •                Break while loop
  •            else
  •                Compute step as negative scaled difference
  •                Compute sum of step and current fine offset
  •                Set fine offset to sum
  •                Write new fine offset to amplifier
  •                Increment counter 2
  •            end if
  •         end while
  •     end for
  • end for
After amplification and offset, the resulting differential voltages are digitized through analog-digital conversion (ADC) and organized into data packets inside a Teensy 3.5 microcontroller. The microcontroller also conditions the data with a 60 Hz notch filter to remove electromagnetic interference (EMI) caused by electric noise generated from wall power. Upon packet construction, the data with four differential voltages are transmitted to a personal computer (PC), where the data is used to compute relative force signals described in Section 2.3. Afterwards, the relative signals are bandpass filtered between 0.5 Hz and 80 Hz with a 2nd-order Butterworth filter to capture normal human motion as well as stroke-related tremor [29,37]. For summary statistics, mean averaging was utilized.

2.6. Sensor Assembly

The sensor assembly is shown in Figure 11 with the finger interface, fixture, and aluminum frames attached to the sensor. (It should be noted for convention’s sake that the upper (top) aluminum frame is connected against the bottom surface of the sensor on the -Z side, whereas the lower (bottom) aluminum frame is fitted against the top surface of the sensor on the +Z side.) The associated free body diagram in Figure 12 demonstrates how applying a fingertip force in the center of the interface generates Z forces and X/Y torques at the sensor surface, which in turn generate bending strain on the sensor surface. Realized assembly photos are displayed in Figure 13. Silicone used for the interface cup was Smooth-On Dragon Skin 10, the plastic cup base was 3D printed from PA 11 nylon, and all structural metal machined was 6061 aluminum alloy. The specific type of FR4 material used is 185HR laminate by Isola Group [31]. The FR4 substrate and copper tracing were provided by Advanced Circuits [38]. The strain sensing packages are S5024 5000 Ω linear gauges provided and installed by Micro-Measurements [39].

2.7. Experiments

2.7.1. Characterization

After full assembly, the full-scale output ( F S O ), sensitivity, and signal-to-noise ratio ( S N R ) of the sensor were characterized. For this test, a sensor is affixed into a 2-stage assembly of perpendicular rotary tables (see Figure 14 and Figure 15). A collection of calibration weights is incrementally suspended from the sensor until at least one axial voltage reports over 10% non-linearity from a dynamic best fit line. This non-linearity represents breakage of the sensor. The maximum loading weight is the capacity, and the resulting voltage reading is the F S O . The sensitivity of the sensor is then the ratio of voltage to weight (measured using a gram-weight scale) at full capacity. A new sensor is then loaded at 50% capacity and the post-digitization signals are measured for fluctuation. The maximum and minimum levels of the fluctuations are then compared to the F S O to calculate the S N R as follows:
S N R = 10 × l o g 10 F S O n o i s e h i g h n o i s e l o w 2
where n o i s e h i g h and n o i s e l o w are the upper and lower voltage limits of the noise fluctuation, respectively.
The resolution of the sensor post-digitization is characterized as
r e s = V F S O V e x c × V A D C b i t s A D C b i t s n o i s e , h i g h b i t s n o i s e , l o w
where V F S O is the full-scale output voltage of the sensor, V e x c is the excitation voltage, V A D C is the operating voltage of the ADC, b i t s n o i s e , h i g h and b i t s n o i s e , l o w are the digital levels of the high and low limits of the noise, respectively, and b i t s A D C are the digital levels of the ADC.

2.7.2. Calibration

To map the relative axial voltages to real-world force units in all spherical directions, a calibration routine for six orthogonal directions is performed. Orthogonal loading directions were chosen since three bases are needed to fully define a 3D Euclidean workspace. Bidirectional loading was also chosen to validate the sign-flipping symmetry of the workspace. For this test, the same 2-stage assembly of perpendicular rotary tables is utilized (see Figure 15). For each direction, the rotary stages are rotated to their corresponding positions. Due to the specific mounting and markings of the stages, the zero-position of the assembly is actually [270°, 344°], which is the home position for the first and second stages, respectively. To establish orthogonality, each loading direction was a multiple of 90 degrees from the zero-position or “home direction”. The four perpendicular loading directions on the X-Y plane had a 45 degree offset from the canonical X-Y axes to verify that interdependence between the X and Y bases was not present. Furthermore, the two opposite loading directions orthogonal to the X-Y plane were established for Z-axis loading. Table 3 shows the corresponding pairs between loading directions and the angular positions of stage 1 and stage 2. After the sensor is installed in the calibration setup, the sensor is progressively weighed in three cycles of increasing and decreasing loads. A summary of the loading procedure is shown in Table 4.
The x, y, and z components of a weighted loading direction for each trial i are calculated as
W R i = W i × cos ( θ 1 , i θ 1 , h o m e ) × cos ( θ 2 , i θ 2 , h o m e ) W i × sin ( θ 1 , i θ 1 , h o m e ) W i × cos ( θ 1 , i θ 1 , h o m e ) × sin ( θ 2 , i θ 2 , h o m e )
where W i is the hanging weight for trial i, θ k , i is the angle of a given stage k for trial i, and the θ k , h o m e is [270°, 344°] for stage 1 and stage 2, respectively. After collecting a set of four readings (each from a corresponding amplifier and strain gauge pair) and appending to an array of voltage readings A measured during the hanging of weights, the least-squares fit calibration matrix C is computed via the Penrose-Moore pseudoinverse
C = W R T ( A + ) T
where
A + = ( A T A ) 1 A T
and W R is the vector of x, y, and z loading components.
This calibration matrix assumes that the columns of A are linearly independent. Given the large number of readings with measurement noise, it is assumed this is true. However, a rank check on A is performed during analysis to verify this assumption.

Nonlinearity Error ( N L )

Nonlinearity error is the sensor’s inability to output a straight line of readings over the full range of inputs, and is computed as
N L = m a x | A T C + W R T | F S O
where A, C, and W R are described in Equation (13), and F S O is described in Equation (10).

Hysteresis Error ( H Y S )

Hysteresis error is the sensor’s inability to output the same reading at a given input depending on whether the input falls or rises to that point, and is calculated as:
H Y S = m a x | A , k ( W i ) A , k ( W i ) | F S O ; k = 2 , 3
where i represents trial index (which also indexes the weight used), k represent cycle index, A , k represents a reading on the increasing load curve, and A , k represents a reading on the decreasing load curve.

Repeatability Error ( R P T )

Repeatability error is the sensor’s inability to output the same reading upon repeated application of the same input, and is calculated as:
R P T = m a x | A , k ( W i ) A , j ( W i ) | F S O ; j , k = 2 , 3 ; j k
where j represents a cycle index to be compared. i, k, and A , k are described in Equation (16).
The cycle index k indicates the current repetition of the cyclic loading, while the trial index i indicates the running count of weights over the whole session. Table 4 details how these indices increment in a loading direction. The reason the first cycle is not considered for error metrics is due to a drifting shift at tare input, known as zero-drift or tare-shift, present from the plasticity of the circuit board.

Accuracy ( A C C )

Accuracy is the root mean squared of the stated error components, and is calculated as:
A C C = N L 2 + H Y S 2 + R P T 2 3
Inspecting the relative-force to real-force (3-by-3) calibration matrix reveals the degree of constitutive anisotropy. The more similar the diagonal entries are and the smaller the off-diagonal entries are, the more isotropic the calibration is. However, the more dissimilar the diagonal entries are, the more orthotropic the material is. Furthermore, the larger the off-diagonal entries are, the more anisotropic the sensor is due to strain in one direction generating signals that represent another direction.

2.7.3. Validation

After obtaining an estimation matrix from calibration, it is checked against non-calibration directions to test its ability to predict loads in the spherical workspace from observed voltages. To systematically test the workspace validation, a random set of 30 rotary stage orientations with associated random weights are generated (see example listing of orientations and loads in Table 5). After the sensor is oriented in a given angular configuration, a zero reading is taken followed by the random load reading. The validation procedure then repeats for the next orientation.
The component-wise residuals R, which compare the predicted weight components from the actual load components, are given by
R = C A T W R T
where W R is the vector of actual load components, C is the calibration matrix, and A is the matrix of voltage readings from the sensor-ADC system. In terms of percent error, Equation (19) becomes
[ % R ] u = 100 % * [ C A T ] u [ W R T ] u [ W R T ] u
where [ · ] u indicates an x, y, or z component of a vector. It should be noted in Equation (20) that if there is zero anticipated load from regression despite a non-zero actual load, then percent error becomes −100%. If the two load types have inverted signs, then percent error is more negative than −100%.
The vector of system-wise residuals R m , which compares overall weight predicted from calibration with actual overall validation weight, are similarly calculated as:
R m = | | C A m T | | 2 W m
where A m and W m are the voltage readings and applied load at configuration m, respectively, and R m is the residual of that configuration m. In terms of percent error, Equation (21) becomes
% R m = 100 % * | | C A m T | | 2 W m W m
Finally, the aggregate root-square of all residuals is computed as:
R R S = | | R m | | 2
where R R S is the 2-norm of the system-wide residuals R m .
When plotting the overall residuals and component errors relative to the angular configuration space, surface fits can predict how non-uniform or even non-linear residual errors are depending on the direction of the force vector. Simple inspection of higher-order polynomial coefficients beyond the fitted intercept can reveal predictive non-uniformity or non-linearity, while R-squared coefficients of determination can indicate how well the fits predict the errors.
All fitting curves and surfaces for validation results are generated with the Statistics and Machine Learning Toolbox in MATLAB 2021a (MathWorks; Natick, MA, USA). More generally, all matrix calculations mentioned above are programmatically computed in a MATLAB 2021a script. The script also generates all the characterization, calibration, and validation plotting results below.

3. Results

A summary of the sensor’s specifications is shown in Table 6.

3.1. Characterization

After fabrication, the sensor unit had a mass of 2.60 g. At breakage, the first test sensor exhibited a capacity of 14.71 N or 1500.0 gf. Effectively, the sensor can read 576.9 times its body weight, a strength-to-weight ratio comparable to a TAS501 S-Type load cell [40]. The sensor also had a maximum full-scale output of 0.289 V. Accordingly, the sensitivity of the device is estimated at a linear 19.64 mV/N or 0.192 mV/gf. The envelope fitting for noise computation is shown in Figure 16. The maximum noise voltage difference of the X, Y, and Z signals was 3.80 mV, which at an ADC reference of 3.3 V and ADC resolution of 16 bits corresponds to 75 digital levels or 6.22 bits. The signal-to-noise ratio was 36.81 dB. The maximum resolution of the sensor-ADC system is 2.564 mN or 0.262 gf.

3.2. Calibration

A plot of the calibration voltage response over the loading routine in the +X, −Y loading direction of the test sensor is shown in Figure 17. (Voltage responses for the other five loading directions are shown in Appendix B.) The full voltage response indicates strong repeatability between cycles and strong linearity overall. A selection of three cycles from the first configuration is shown in Figure 18, which identify mild hysteresis in the first cycle and the lack of hysteresis in the second and third cycles. The associated sensitivity matrix for mapping loads to voltages is depicted in Table 7. For this test sensor, the hysteresis distance between loading curves on the same cycle is 2.29% of the full-scale voltage. The maximum deviation from the best fit curve is 2.94% FSO, and the maximum variation of repeated readings is 0.61% FSO. Aggregating these metrics using a root-mean-sum-square formula yields an average accuracy of 2.18% FSO.

3.3. Validation

Figure 19 details the loading validation by comparing the predicted axial-component of loading (from the calibration matrix) with the actual axial-components of loading, which are elements of the real number set within plus or minus capacity. Overall composite residuals as a function of overall loading are depicted in Figure 20.
Component-wise, the coefficients of determination from fitting an ideal 1:1 weight prediction curve to the X, Y, and Z components are 88.46%, 96.05%, and 99.12% respectively. System-wise, the coefficient of determination from fitting an ideal weight prediction curve is 96.98%, whereas the coefficient of determination from a best-fit curve is 97.06%. The aggregate root-square of all residuals for the tested sensor was 6.21%. The maximum residual error during validation was 211.4 g from the predicted value of 930.0 g.
Component and overall residual errors as a percentage of applied load are shown within an angular configuration space (or C-space) in Figure 21. The coefficients of the mixed-quadratic polynomial surface fits for each residual error type, as well as associated coefficients of determination, are shown in Table 8. Summary statistics of error by load and a histogram of residual error are provided in Appendix C.

4. Discussion

The sensor was successfully manufactured to fit within the fixture geometry shown in Figure 1. The sensor’s mass of 2.6 g after manufacture exceeds the design goal of 10 g, which helps any larger system the sensor would be integrated into stay lightweight and portable. While the final resolution of 2.564 mN does not meet the design resolution of 1 mN, the sensor system is still capable of capturing fine hand tremor signals given that certain tasks only need 20 mN or less of precision [29]. The design result and goal are in the same order of magnitude, and the resolution can be improved with better filtering and higher-bit ADC units. Additionally, while the sensor’s capacity of 14.71 N may have not met the design goal for maximum voluntary contraction measurements, it can capture forces of moderate exertion without any issue. Furthermore, the repeatability error of 0.61% FSO, the hysteresis error of 2.29% FSO, and the nonlinearity error of 2.94% FSO, all do not exceed the design goal of below 5%. This demonstrates that the sensor transduces loads consistently and linearly within and between applications. Combining these errors into an accuracy metric of 2.18% meets the design goal and demonstrates good general sensor accuracy. Due in part to the rigidity of the metal components used to assemble the system, the sensor is able to collect signals without flexing more than 5 mm from neutral state. This was verified with calipers as the sensor was being loaded. Finally, the sensor cost estimate of $350 is well below the stated goal of $1000, making this sensor solution an ideal alternative compared to other systems in the market. Comparisons between the design sensor and current sensors used in finger rehabilitation are given in Table 9. The design sensor is not intended to completely supersede all metrics of any of these sensors and to replace them in the field, but highlight that the proposed sensor design addresses a missing implementation gap in 3D finger sensing. Overall, the sensor’s ability to record 3D forces of small to relatively large magnitude in a light package highlights a great design success.
Despite the general success of characterization and calibration, there are some caveats that warrant discussion. Higher fluctuations in noise and repeatability will affect the system’s ability to distinguish finer motion signals. Stated accuracy metrics from calibration were achieved only after voltage readings from the first cycle of loading were excluded from analysis. Plastic deformation of the tare reading when a configuration is first loaded could explain why the loading curves’ zeroes shift to another configuration. Although hysteresis error within a configuration can meet design requirements, that does not imply the zero-load output will not shift significantly between configurations. The plastic zero-shifting could be resolved by improving the bonding technique and the substrate material, although verification is yet to be done. Additionally, validation of the sensor in the entire 3D workspace of possible applied forces is a more significant challenge. The root-mean-sum-square aggregate of the overall residuals does not meet the target of 5% FSO in the specification criteria. Unusually high residuals may partially be due to the aforementioned plastic tare shifting. Since each of the 30 random readings is its own configuration, a single cycle of loading may be insufficient to fix nonlinearity issues. Another plausible reason is that relative errors from validation loading over the entire C-space is not as uniform as previously understood. This reason is supported by the presence of high error outliers shown in Figure 21. For example, there is a cluster of −100% relative X errors from validating around angular configuration [270°, 254°]. Furthermore, applying surface polynomial fits to the C-space residuals do not exhibit a uniform “horizontal” plane. Instead of just having zero-order affine coefficients, all 2nd-order surface fits in Table 8 exhibit non-zero linear, mixed-interaction, and squared terms. The associated R-squared metrics imply correlations that are weak (component-wise) to moderate (overall). Generating higher-order surface fits that predict residual errors from angular position better would further demonstrate that the C-space is not just non-uniform in error, but also incredibly non-linear. Despite low nonlinearity error being demonstrated across levels of calibration weights, it is still possible for nonlinear residuals to appear across non-calibration loading angles during validation. The introduction of holes in the substrate or the transversely isotropic structure may mechanically explain these non-uniform and non-linear predicted load error trends across the C-space. Furthermore, cross-sensitivity can exist on this sensor due to the lack of mechanical isolation between the strain gauges on the sensor substrate, which may explain the non-zero values in non-diagonal entries of the calibration matrix.
In the future, the dynamics of the sensor should be explored by collecting metrics such as rise/fall time after a step input, vibrational parameters, and fatigue cycles before breakage. Given that human subjects will operate any peripheral device with the sensor in a constantly dynamic manner, understanding mechanics beyond steady-state conditions is important. Furthermore, accuracy errors and residual errors could be improved by either switching the substrate to a very elastic element (sheet metal) with good insulation for signal routing, or by switching to a more robust layout of measuring bending deflection (e.g., capacitive sensing with a stiff dielectric as demonstrated by Jeong et al. [41]). By moving away from a brittle material like fiberglass resin, the sensor should exhibit less variation in plastic tare-shift between configurations and in maximum input capacity. Additionally, we can investigate the effects of cross-sensitivity between gauges with FEA modeling with forces in all axes, and with a semi-supervised method of regressing the calibration matrix rather than a simple linear regression.
Despite these limitations, the presented sensor offers an inexpensive and relatively accurate sensing package for capturing the fine forces generated at the fingertips in all directions. By implementing the sensors in a larger rehabilitative device, dexterous finger control, such as finger individuation and object manipulation, can be assessed and potentially trained independent of hand weakness [4,42]. Precision control can be tracked and improved without needing a separate hand model to estimate task-space forces from joint torques. The possibility of training hand dexterity despite marked weakness increases the chances for patients to regain control needed for activities of daily living (e.g., manipulating and tapping an egg, single-tapping a touchscreen, twisting a dial on a car console, etc.).

Author Contributions

Conceptualization, K.O., J.X. and J.C.; Methodology, K.O. and J.C.; Software, J.C.; Formal Analysis, J.C.; Investigation, J.C.; Resources, K.O., J.D.B., J.W.K. and J.X.; Data Curation, J.C.; Writing—Original Draft Preparation, J.C.; Writing—Review & Editing, J.D.B., J.W.K. and J.X.; Visualization, J.C.; Project Administration, K.O., J.W.K. and J.D.B.; Funding Acquisition, J.W.K. All authors have read and agreed to the published version of the manuscript.

Funding

Early work on this research was funded by the James S. McDonnell Foundation in 2010 and 2013.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://github.com/Jake-Carducci/fingersensorpaper (27 September 2022).

Acknowledgments

Thanks to Omar Ahmad, Promit Roy, Alex Forrence, and Timothy Ma for assisting with hardware and firmware development, and thanks to Aiden Devaney, Sapna Kumar, and Teresa George for assisting in sensor calibration and validation.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DOFDegrees of freedom
FEAFinite element analysis
dimDimensions
FSOFull scale output
DCDirect-current
EMIElectromagnetic interference
PCPersonal computer
SNRSignal-to-noise ratio
ADCAnalog-to-digital conversion
CalibCalibration
C-spaceConfiguration space

Appendix A

Figure A1. FEA simulation of a 10.0 N force in the X-direction applied 15.5 mm from the center of the 3-DOF sensor. Applied forces are represented by pink arrows with the moment arm in dashed pink lines, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Figure A1. FEA simulation of a 10.0 N force in the X-direction applied 15.5 mm from the center of the 3-DOF sensor. Applied forces are represented by pink arrows with the moment arm in dashed pink lines, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Sensors 22 07441 g0a1
Figure A2. FEA simulation of a 10.0 N force in the Y-direction applied 15.5 mm from the center of the 3-DOF sensor. Applied forces are represented by pink arrows with the moment arm in dashed pink lines, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Figure A2. FEA simulation of a 10.0 N force in the Y-direction applied 15.5 mm from the center of the 3-DOF sensor. Applied forces are represented by pink arrows with the moment arm in dashed pink lines, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Sensors 22 07441 g0a2

Appendix B

Figure A3. Voltage response to 3D calibration loading in different configurations.
Figure A3. Voltage response to 3D calibration loading in different configurations.
Sensors 22 07441 g0a3
Figure A4. Voltage response to 3D calibration loading in different configurations.
Figure A4. Voltage response to 3D calibration loading in different configurations.
Sensors 22 07441 g0a4
Figure A5. Voltage response to 3D calibration loading in different configurations.
Figure A5. Voltage response to 3D calibration loading in different configurations.
Sensors 22 07441 g0a5

Appendix C

Figure A6. Central measures and spread of component-wise residuals as a function of the overall actual load applied. For this sensor, each load has 3–4 associated configurations each with a residual error. (a) Means of multiple residuals tied to each overall load are represented by the centerline, while standard deviation is represented by error bars, extending above and below by 1 s.d. (b) Medians of residuals tied to each overall load are represented as a centerline in the envelope plot, 25-percent and 75-percent quartiles are represented by a darker region, and upper/lower extremes are represented by lighter regions.
Figure A6. Central measures and spread of component-wise residuals as a function of the overall actual load applied. For this sensor, each load has 3–4 associated configurations each with a residual error. (a) Means of multiple residuals tied to each overall load are represented by the centerline, while standard deviation is represented by error bars, extending above and below by 1 s.d. (b) Medians of residuals tied to each overall load are represented as a centerline in the envelope plot, 25-percent and 75-percent quartiles are represented by a darker region, and upper/lower extremes are represented by lighter regions.
Sensors 22 07441 g0a6
Figure A7. Histogram of overall residual error. The bin size of 11 g.
Figure A7. Histogram of overall residual error. The bin size of 11 g.
Sensors 22 07441 g0a7

References

  1. Jankovic, J. Parkinson’s disease: Clinical features and diagnosis. J. Neurol. Neurosurg. Psychiatry 2008, 79, 368–376. [Google Scholar] [CrossRef] [PubMed]
  2. Francisco, G.E.; McGuire, J.R. Poststroke Spasticity Management. Stroke 2012, 43, 3132–3136. [Google Scholar] [CrossRef] [PubMed]
  3. Charles, S.K.; Okamura, A.M.; Bastian, A.J. Does a basic deficit in force control underlie cerebellar ataxia? J. Neurophysiol. 2013, 109, 1107–1116. [Google Scholar] [CrossRef] [PubMed]
  4. Jing Xu, X.; Naveed Ejaz, X.; Hertler, B.; Branscheidt, M.; Widmer, M.; Faria, A.V.; Harran, M.D.; Cortes, J.C.; Kim, N.; Celnik, P.A.; et al. Separable systems for recovery of finger strength and control after stroke Even healthy people show in-voluntary force production (enslaving) on the uninstructed. J. Neurophysiol. 2017, 118, 1151–1163. [Google Scholar] [CrossRef]
  5. Dovat, L.; Lambercy, O.; Gassert, R.; Maeder, T.; Milner, T.; Leong, T.C.; Burdet, E. HandCARE: A cable-actuated rehabilitation system to train hand function after stroke. IEEE Trans. Neural Syst. Rehabil. Eng. 2008, 16, 582–591. [Google Scholar] [CrossRef] [PubMed]
  6. Ito, S.; Kawasaki, H.; Ishigure, Y.; Natsume, M.; Mouri, T.; Nishimoto, Y. A design of fine motion assist equipment for disabled hand in robotic rehabilitation system. J. Frankl. Inst. 2011, 348, 79–89. [Google Scholar] [CrossRef]
  7. Heo, J.S.; Chung, J.H.; Lee, J.J. Tactile sensor arrays using fiber Bragg grating sensors. Sens. Actuators A Phys. 2006, 126, 312–327. [Google Scholar] [CrossRef]
  8. Da Silva, J.; de Carvalho, A.; da Silva, D. A strain gauge tactile sensor for finger-mounted applications. IEEE Trans. Instrum. Meas. 2002, 51, 18–22. [Google Scholar] [CrossRef]
  9. Castellini, C.; Ravindra, V. A wearable low-cost device based upon Force-Sensing Resistors to detect single-finger forces. In Proceedings of the 5th IEEE RAS/EMBS International Conference on Biomedical Robotics and Biomechatronics, Sao Paulo, Brazil, 12–15 August 2014; pp. 199–203. [Google Scholar] [CrossRef]
  10. Nakatani, M.; Shiojima, K.; Kinoshita, S.; Kawasoe, T.; Koketsu, K.; Wada, J. Wearable contact force sensor system based on fingerpad deformation. In Proceedings of the 2011 IEEE World Haptics Conference, Istanbul, Turkey, 21–24 June 2011; pp. 323–328. [Google Scholar] [CrossRef]
  11. Park, J.; Heo, P.; Kim, J.; Na, Y. A Finger Grip Force Sensor with an Open-Pad Structure for Glove-Type Assistive Devices. Sensors 2019, 20, 4. [Google Scholar] [CrossRef]
  12. Ponce Wong, R.D.; Posner, J.D.; Santos, V.J. Flexible microfluidic normal force sensor skin for tactile feedback. Sens. Actuators A Phys. 2012, 179, 62–69. [Google Scholar] [CrossRef]
  13. Troncossi, M.; Mozaffari-Foumashi, M.; Parenti-Castelli, V. An Original Classification of Rehabilitation Hand Exoskeletons; Technical Report 4; Verizona Publisher Limited: London, UK, 2016. [Google Scholar]
  14. Sarakoglou, I.; Tsagarakis, N.G.; Caldwell, D.G. Occupational and physical therapy using a hand exoskeleton based exerciser. In Proceedings of the 2004 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Sendai, Japan, 28 September–2 October 2004; Volume 3, pp. 2973–2978. [Google Scholar] [CrossRef]
  15. Ben-Tzvi, P.; Ma, Z. Sensing and Force-Feedback Exoskeleton (SAFE) Robotic Glove. IEEE Trans. Neural Syst. Rehabil. Eng. 2015, 23, 992–1002. [Google Scholar] [CrossRef] [PubMed]
  16. Brokaw, E.B.; Black, I.; Holley, R.J.; Lum, P.S. Hand Spring Operated Movement Enhancer (HandSOME): A portable, passive hand Exoskeleton for stroke rehabilitation. IEEE Trans. Neural Syst. Rehabil. Eng. 2011, 19, 391–399. [Google Scholar] [CrossRef] [PubMed]
  17. Winter, S.H.; Bouzit, M. Testing and usability evaluation of the MRAGES force feedback glove. In Proceedings of the Fifth International Workshop on Virtual Rehabilitation, IWVR 2006, New York, NY, USA, 29–30 August 2006; pp. 82–87. [Google Scholar] [CrossRef]
  18. Wege, A.; Hommel, G. Development and control of a hand exoskeleton for rehabilitation of hand injuries. In Proceedings of the 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems, IROS, Edmonton, AB, Canada, 2–6 August 2005; pp. 3046–3051. [Google Scholar] [CrossRef]
  19. Lambercy, O.; Dovat, L.; Gassert, R.; Burdet, E.; Teo, C.L.; Milner, T. A haptic knob for rehabilitation of hand function. IEEE Trans. Neural Syst. Rehabil. Eng. 2007, 15, 356–366. [Google Scholar] [CrossRef] [PubMed]
  20. Shirakawa, K.; Tanaka, Y.; Hashimoto, M.; Watarai, E.; Igarashi, T. Wearable Artificial Fingers with Skin Vibration and Multi-Axis Force Sensors. IEEE Trans. Haptics 2021, 14, 242–247. [Google Scholar] [CrossRef] [PubMed]
  21. Hale, N.; Valero, M.; Tang, J.; Moser, D.; Jiang, L. A preliminary study on characterisation of finger interface kinetics using a pressure and shear sensor system. Prosthetics Orthot. Int. 2018, 42, 60–65. [Google Scholar] [CrossRef] [Green Version]
  22. Mascaro, S.; Asada, H. Measurement of Finger Posture and Three-Axis Fingertip Touch Force Using Fingernail Sensors. IEEE Trans. Robot. Autom. 2004, 20, 26–35. [Google Scholar] [CrossRef]
  23. ATI Industrial Automation. F/T Sensor: Nano17. Available online: https://www.ati-ia.com/products/ft/ft_models.aspx?id=Nano17 (accessed on 13 January 2022).
  24. Dai, Z.; Wang, Z.; Ji, A. Dynamics of gecko locomotion: A force-measuring array to measure 3D reaction forces. J. Exp. Biol. 2011, 214, 703–708. [Google Scholar] [CrossRef]
  25. Reinhardt, L.; Blickhan, R. Ultra miniature force plate for measuring triaxial forces in the micro newton range. J. Exp. Biol. 2013. [Google Scholar] [CrossRef]
  26. Kim, S.; Kim, C.; Park, S.; Lee, D.Y. A 3-DOF sensor to estimate the force applied to the tip of a surgical instrument. In Proceedings of the 2017 18th International Conference on Advanced Robotics (ICAR), Hong Kong, China, 10–12 July 2017; pp. 143–148. [Google Scholar] [CrossRef]
  27. Sun, Y. Fingertip Force and Contact Position and Orientation Sensor. In Proceedings of the IEEE International Conference on Robotics and Automation, Shanghai, China, 9–13 May 2013; pp. 1114–1119. [Google Scholar]
  28. Brookhuis, R.A.; Wiegerink, R.J.; Lammerink, T.S.; Ma, K.; Krijnen, G.J. Large range multi-axis fingertip force sensor. In Proceedings of the Transducers and Eurosensors XXVII: The 17th International Conference on Solid-State Sensors, Actuators and Microsystems, Transducers and Eurosensors 2013, Barcelona, Spain, 16–20 June 2013; pp. 2737–2740. [Google Scholar] [CrossRef]
  29. Endo, H.; Kawahara, K. Relationship between hand stability and the 10-Hz physiological tremor during various manual tasks. Ergonomics 2010, 53, 491–501. [Google Scholar] [CrossRef]
  30. Bretz, K.J.; Jobbágy, Á.; Bretz, K. Force measurement of hand and fingers. Biomech. Hung. 2010, 3, 61–66. [Google Scholar] [CrossRef]
  31. Isola Group. 185HR. Available online: https://www.isola-group.com/pcb-laminates-prepreg/185hr-laminate-and-prepreg/ (accessed on 13 January 2022).
  32. C.D. Wälzholz GmbH & Co. Spring Steel with High Yield Strength. Available online: https://www.waelzholz.com/en/steel-materials/cold-rolled-steel-strip/spring-steel.html (accessed on 15 June 2022).
  33. Mindlin, R.D. Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates. In The Collected Papers of Raymond D. Mindlin Volume I; Springer New York: New York, NY, USA, 1989; pp. 225–232. [Google Scholar] [CrossRef]
  34. Isola. Understanding Glass Fabric; Technical Report; Isola: Chandler, AZ, USA, 2021. [Google Scholar]
  35. Texas Instruments. PGA309 Noise Filtering; Technical Report; Texas Instruments Inc.: Dallas, TX, USA, 2005. [Google Scholar]
  36. Texas Instruments. PGA309 Voltage Output Programmable Sensor Conditioner User’s Guide; Technical Report; Texas Instruments Inc.: Dallas, TX, USA, 2011. [Google Scholar]
  37. Mcauley, J.H.; Marsden, C.D. Physiological and pathological tremors and rhythmic central motor control. Brain 2000, 123, 1545–1567. [Google Scholar] [CrossRef] [PubMed]
  38. Advanced Circuits. Available PCB Materials. Available online: https://www.4pcb.com/laminate-material-vendors.html (accessed on 21 September 2022).
  39. Micro-Measurements. S5024. Available online: https://micro-measurements.com/pca/detail/s5024 (accessed on 21 September 2022).
  40. HT Sensor Technology, Co., Ltd. TAS501 S Type Load Cell. Available online: http://www.htc-sensor.com/products/105.html (accessed on 13 January 2022).
  41. Jeong, H.; Choi, K.; Park, S.J.; Park, C.H.; Choi, H.R.; Kim, U. Rugged and compact three-axis force/torque sensor for wearable robots. Sensors 2021, 21, 2770. [Google Scholar] [CrossRef] [PubMed]
  42. Mawase, F.; Cherry-Allen, K.; Xu, J.; Anaya, M.; Uehara, S.; Celnik, P. Pushing the Rehabilitation Boundaries: Hand Motor Impairment Can Be Reduced in Chronic Stroke. Neurorehabilit. Neural Repair 2020, 34, 733–745. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Bounding geometry of the sensor fixture used to fit the sensor to the hand rehabilitation device.
Figure 1. Bounding geometry of the sensor fixture used to fit the sensor to the hand rehabilitation device.
Sensors 22 07441 g001
Figure 2. The aluminum fixture that fits and grounds the 3-DOF sensor. The sensor is secured to the fastening holes from within the sensor cavity. The long fixture body allows large-scale adjustment of the sensor to accommodate the user’s unique neutral hand posture.
Figure 2. The aluminum fixture that fits and grounds the 3-DOF sensor. The sensor is secured to the fastening holes from within the sensor cavity. The long fixture body allows large-scale adjustment of the sensor to accommodate the user’s unique neutral hand posture.
Sensors 22 07441 g002
Figure 3. The finger-cup interface attached to a two-part aluminum frame assembly. The 3-DOF sensor fits between the two halves of the frame assembly, which secure together to create a rigid structure on the sensor. A plastic base attaches to the top of this structure to provide a rigid base for a silicone cup. Once secure (using adhesive), the silicone cup provides a comfortable grip onto a user’s finger. The finger forces (and torques) can then be transmitted in total to the body of the sensor.
Figure 3. The finger-cup interface attached to a two-part aluminum frame assembly. The 3-DOF sensor fits between the two halves of the frame assembly, which secure together to create a rigid structure on the sensor. A plastic base attaches to the top of this structure to provide a rigid base for a silicone cup. Once secure (using adhesive), the silicone cup provides a comfortable grip onto a user’s finger. The finger forces (and torques) can then be transmitted in total to the body of the sensor.
Sensors 22 07441 g003
Figure 4. Exploded view illustration of strain gauge units bonding to both sides of sensor board. Stain gauges are placed from above and underneath the sensor body to a location on the surface indicated by the arrows. Holes, for purposes of fastening to the aluminum fixture and for fitting alignment pins between the frame assembly, are rendered in the sensor body for reference.
Figure 4. Exploded view illustration of strain gauge units bonding to both sides of sensor board. Stain gauges are placed from above and underneath the sensor body to a location on the surface indicated by the arrows. Holes, for purposes of fastening to the aluminum fixture and for fitting alignment pins between the frame assembly, are rendered in the sensor body for reference.
Sensors 22 07441 g004
Figure 5. Schematic drawing representing a pair of strain gauges connected in a voltage divider circuit on the 3-DOF sensor. Ideally, the two gauges experience equal and opposite resistance (as shown) resulting in the voltage reading at the probe location having a “doubling” effect compared to only a single-gauge configuration.
Figure 5. Schematic drawing representing a pair of strain gauges connected in a voltage divider circuit on the 3-DOF sensor. Ideally, the two gauges experience equal and opposite resistance (as shown) resulting in the voltage reading at the probe location having a “doubling” effect compared to only a single-gauge configuration.
Sensors 22 07441 g005
Figure 6. The numerical references of gauge locations 1 to 4 on the (a) top surface and (b) bottom surface of the sensor. A pair of gauges on top and bottom will share a reference number, but will appear at “mirrored” locations after flipping the body of the sensor. These reference numbers are useful in describing the relationship between Euclidean-space force components and changes in probe voltages between gauge pairs in the voltage divider configurations.
Figure 6. The numerical references of gauge locations 1 to 4 on the (a) top surface and (b) bottom surface of the sensor. A pair of gauges on top and bottom will share a reference number, but will appear at “mirrored” locations after flipping the body of the sensor. These reference numbers are useful in describing the relationship between Euclidean-space force components and changes in probe voltages between gauge pairs in the voltage divider configurations.
Sensors 22 07441 g006
Figure 7. Sub-unit of the 3-DOF sensor with one strain gauge (indicated by the white arrow), two alignment pins, and two fastening holes. A 2.5 N downwards force was simulated, represented by pink arrows at the edge of the sub-unit. Fixed surfaces on the holes are represented by brown arrows. Resulting strain of the sub-unit surface, with respect to the white axis in the u direction, can be negative and compressive (which is colored more blue) or can be positive and tensile (which is colored more red).
Figure 7. Sub-unit of the 3-DOF sensor with one strain gauge (indicated by the white arrow), two alignment pins, and two fastening holes. A 2.5 N downwards force was simulated, represented by pink arrows at the edge of the sub-unit. Fixed surfaces on the holes are represented by brown arrows. Resulting strain of the sub-unit surface, with respect to the white axis in the u direction, can be negative and compressive (which is colored more blue) or can be positive and tensile (which is colored more red).
Sensors 22 07441 g007
Figure 8. Geometry of disk approximated from four strain sub-units arranged in diamond-shaped configuration. Each sub-unit is distinguished by a diagonal striping pattern; subunit areas are represented by unique rectangular patterns. Pink diamonds indicate strain gauge placements. Alignment and fastening holes are rendered for reference.
Figure 8. Geometry of disk approximated from four strain sub-units arranged in diamond-shaped configuration. Each sub-unit is distinguished by a diagonal striping pattern; subunit areas are represented by unique rectangular patterns. Pink diamonds indicate strain gauge placements. Alignment and fastening holes are rendered for reference.
Sensors 22 07441 g008
Figure 9. Illustration of applied forces and reaction forces from fixed supports to sensor board, as well as strain-sensing directions. Applied forces in this scenario are all into the page, designated by purple crossed circles. Fixed supports are all out of the page, designated by lime-green dotted circles.
Figure 9. Illustration of applied forces and reaction forces from fixed supports to sensor board, as well as strain-sensing directions. Applied forces in this scenario are all into the page, designated by purple crossed circles. Fixed supports are all out of the page, designated by lime-green dotted circles.
Sensors 22 07441 g009
Figure 10. FEA simulation of a vertical 10.0 N downwards force distributed around the strain gauges on the 3-DOF sensor. Applied forces are represented by pink arrows pointed down, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Figure 10. FEA simulation of a vertical 10.0 N downwards force distributed around the strain gauges on the 3-DOF sensor. Applied forces are represented by pink arrows pointed down, while support fixtures are represented by brown arrows. Surface coordinates u and vs. are represented by a grey frame. All results are viewed from top surface. Results include (a) von Mises stress field, (b) u-direction strain, and (c) v-direction strain fields. The axis being referenced for strain measurement is indicated by a white bidirectional arrow on the surface of the sensor at the gauge locations. The two strain axes are perpendicular, but are offset from established X-Y axes by 45 degrees.
Sensors 22 07441 g010
Figure 11. Side view of final assembled sensor-fixture-interface system. The 3-DOF sensor fits and aligns between the aluminum frames, and is fastened inside the beam structure. The plastic base holds the silicone cup for human interfacing, and is held against the top aluminum frame by a retaining spring.
Figure 11. Side view of final assembled sensor-fixture-interface system. The 3-DOF sensor fits and aligns between the aluminum frames, and is fastened inside the beam structure. The plastic base holds the silicone cup for human interfacing, and is held against the top aluminum frame by a retaining spring.
Sensors 22 07441 g011
Figure 12. Free-body diagram of the human finger and the 3-DOF sensor within the overall assembly, viewed from the side (a) and front (b). Total force exerted by the finger onto the assembly is broken down into axial components, designated by green arrows. The resultant finger force is represented by a red arrow. The axial directions are indicated by black coordinate frames to the side. Reaction force and torque are also indicated at the sensor by blue arrows. The reaction components contribute to bending of the sensor, which are sensed by the gauges.
Figure 12. Free-body diagram of the human finger and the 3-DOF sensor within the overall assembly, viewed from the side (a) and front (b). Total force exerted by the finger onto the assembly is broken down into axial components, designated by green arrows. The resultant finger force is represented by a red arrow. The axial directions are indicated by black coordinate frames to the side. Reaction force and torque are also indicated at the sensor by blue arrows. The reaction components contribute to bending of the sensor, which are sensed by the gauges.
Sensors 22 07441 g012
Figure 13. Realized implementation of overall finger assembly viewed from top and underneath.
Figure 13. Realized implementation of overall finger assembly viewed from top and underneath.
Sensors 22 07441 g013
Figure 14. The loading materials needed to characterize, calibrate, and validate a 3-DOF sensor: (a) 2-stage rotary table setup for configuring 3D angular orientation of sensor; (b) set of loading weights (ranging from 1 g to 2 kg) and loading tube for holding weights and providing line of tension from the sensor.
Figure 14. The loading materials needed to characterize, calibrate, and validate a 3-DOF sensor: (a) 2-stage rotary table setup for configuring 3D angular orientation of sensor; (b) set of loading weights (ranging from 1 g to 2 kg) and loading tube for holding weights and providing line of tension from the sensor.
Sensors 22 07441 g014
Figure 15. Assembled two-stage loading setup with (a) overall view and (b) zoomed view of sensor installed in second stage. The second stage is mounted perpendicular to the first stage. The finger assembly is mounted to the second stage such that the sensor’s Z-axis is pointed perpendicular to the normal axis of the second stage. The tension line of the loading tube is secured to the top of the sensor using a custom plastic fastener.
Figure 15. Assembled two-stage loading setup with (a) overall view and (b) zoomed view of sensor installed in second stage. The second stage is mounted perpendicular to the first stage. The finger assembly is mounted to the second stage such that the sensor’s Z-axis is pointed perpendicular to the normal axis of the second stage. The tension line of the loading tube is secured to the top of the sensor using a custom plastic fastener.
Sensors 22 07441 g015
Figure 16. Voltage response of one voltage divider unit during constant 50% loading capacity. A fitted envelope has a green upper bound at the highest noise level and a blue lower bound at the lowest noise level. The height of the envelope represents the magnitude of noise in the signal, which here is 5.79 mV. Note that the maximum noise of 3.8 mV when calculated from X, Y, and Z signals is lower than noise calculated from raw divider signals due to how X, Y, and Z voltages are computed from raw voltages.
Figure 16. Voltage response of one voltage divider unit during constant 50% loading capacity. A fitted envelope has a green upper bound at the highest noise level and a blue lower bound at the lowest noise level. The height of the envelope represents the magnitude of noise in the signal, which here is 5.79 mV. Note that the maximum noise of 3.8 mV when calculated from X, Y, and Z signals is lower than noise calculated from raw divider signals due to how X, Y, and Z voltages are computed from raw voltages.
Sensors 22 07441 g016
Figure 17. Voltage response to 3D calibration loading in the +X −Y direction. All three forward-backward cycles are plotted in dashed lines, and corresponding lines of best fit are solid.
Figure 17. Voltage response to 3D calibration loading in the +X −Y direction. All three forward-backward cycles are plotted in dashed lines, and corresponding lines of best fit are solid.
Sensors 22 07441 g017
Figure 18. Voltage response to 3D calibration loading in 3 separate cycles in the +X −Y loading direction.
Figure 18. Voltage response to 3D calibration loading in 3 separate cycles in the +X −Y loading direction.
Sensors 22 07441 g018
Figure 19. Axial components of actual loading compared with loads predicted by the calibration matrix. The black line indicates ideal estimation matching with actual load.
Figure 19. Axial components of actual loading compared with loads predicted by the calibration matrix. The black line indicates ideal estimation matching with actual load.
Sensors 22 07441 g019
Figure 20. Overall loading residuals as a function of the overall loads applied. The discrete levels of loading are highlighted as orange-yellow vertical lines. The dashed grey line indicates zero error and that actual overall mass matches what is predicted by calibration. Meanwhile, the red diagonal line indicates the line of best fit.
Figure 20. Overall loading residuals as a function of the overall loads applied. The discrete levels of loading are highlighted as orange-yellow vertical lines. The dashed grey line indicates zero error and that actual overall mass matches what is predicted by calibration. Meanwhile, the red diagonal line indicates the line of best fit.
Sensors 22 07441 g020
Figure 21. Relative percents of overall loading residuals and axial-component errors as a function of C-space angular orientation, normalized to either the total applied load or respective applied load component (x, y, or z). Magenta crosses indicate the positions at which the sensor was calibrated (with mirrored counterparts), while the circles indicates at what angular positions the errors were sampled and how large they were (indicated by color) for the component or system. Relative component error scaling is limited to 150% in magnitude for easier interpretation of color differences at smaller errors. Component errors more negative than −100% indicate actual and anticipated loads used to calculate error have inverted signs.
Figure 21. Relative percents of overall loading residuals and axial-component errors as a function of C-space angular orientation, normalized to either the total applied load or respective applied load component (x, y, or z). Magenta crosses indicate the positions at which the sensor was calibrated (with mirrored counterparts), while the circles indicates at what angular positions the errors were sampled and how large they were (indicated by color) for the component or system. Relative component error scaling is limited to 150% in magnitude for easier interpretation of color differences at smaller errors. Component errors more negative than −100% indicate actual and anticipated loads used to calculate error have inverted signs.
Sensors 22 07441 g021
Table 1. Tabular summary of a selection of sensors used in upper-limb rehabilitation.
Table 1. Tabular summary of a selection of sensors used in upper-limb rehabilitation.
MetricUnitDiedrichsen [4]HANDCare [5]FFPO [27]MESA+ [28]
Sensor Honeywell FSG15N1AAngst-Pfister MilliNewtonATI Nano17Wafer ring
DOF 1165
Massg3N/A9.075.93
SpanN150.4 to 21260
Accuracy% FSO *2.251N/AN/A
Hysteresis% FSO *3N/AN/AN/A
Nonlinearity% FSO *1.5N/AN/A1
Repeatability% FSO *21N/AN/A
Cost$118N/A5000N/A
Planar dimmm12.7 × 18.225.4179
Depthmm912.714.51
N/A: Specifications are not publicly known. * FSO: Full scale output.
Table 2. Default parameter settings for the PGA309 amplifier transfer function (see Equation (8)) converting a pair of probing voltages from a gauge bridge into a larger output voltage.
Table 2. Default parameter settings for the PGA309 amplifier transfer function (see Equation (8)) converting a pair of probing voltages from a gauge bridge into a larger output voltage.
SettingValue
G i , 0 32
G f , 0 1
G o , 0 9
V f , 0 3.33 / 18
V c , 0 0
Table 3. Mapping between cardinal loading direction and angle locations on rotary stages. θ 1 and θ 2 refer to the first and second stage angles, respectively.
Table 3. Mapping between cardinal loading direction and angle locations on rotary stages. θ 1 and θ 2 refer to the first and second stage angles, respectively.
Loading DirectionStage Angle (°)
XYZ θ 1 θ 2
++N/A135344
+N/A45344
N/A315344
+N/A225344
N/AN/A+270254
N/AN/A27074
+: positive loading direction, −: negative loading direction, N/A: no loading direction.
Table 4. Loading cycle trials for a single loading direction. The mass of the tray was 230 grams.
Table 4. Loading cycle trials for a single loading direction. The mass of the tray was 230 grams.
TrialCycleWeight + Tray (g)
110
21530
311030
41530
520
62530
721030
82530
930
103530
1131030
123530
1330
Table 5. Example loading cycle routines for a validation protocol. Each cycle has two rotary stage angles and a loading weight. Stated weights do not include 230 g loading tray.
Table 5. Example loading cycle routines for a validation protocol. Each cycle has two rotary stage angles and a loading weight. Stated weights do not include 230 g loading tray.
θ 1 (°) θ 2 (°)W (g)
259194400
28348700
83195100
71309800
33771100
25656300
2522600
259238500
527400
300105200
3233510
90275200
3788600
297246300
8650100
68227700
54309800
19324200
239125600
1091750
293245500
177253500
1041660
17131800
58101300
23727700
111160200
33360400
15144400
120331500
Table 6. Experimentally-obtained specifications for a representative 3-DOF fingertip sensor, with comparison to original requirements.
Table 6. Experimentally-obtained specifications for a representative 3-DOF fingertip sensor, with comparison to original requirements.
SpecificationUnitsExperimental
Value
Requirement
Sensor massg2.65
Range/capacityN14.7120
Full scale outputV0.289N/A
SNRdB36.8N/A
Noisebits6.22N/A
Final resolution *mN2.5641
gf0.2620.101
Repeatability **% FSO0.615
Hysteresis **% FSO2.295
Nonlinearity **% FSO2.945
Accuracy **% FSO2.185
Root-square residual% range6.215
Dimensionsmm23.00 × 16.18 × 10.24N/A
Estimated unit cost$3501000
* Assumes noise fluctuations are negligible after averaging filter. ** First cycle removed to avoid initial plastic zero drift effects. N/A.
Table 7. 3-by-3 sensitivity matrix of test sensor. For any given value, its column represents the axial component of the voltage reading being mapped from, and its row represents the axial component of the estimated weight being mapped to. For example, −180.02 in row 1, column 2 is the scaling factor applied to the Y-component of the voltage reading to get a least-squares regression estimate of the X-component of the load.
Table 7. 3-by-3 sensitivity matrix of test sensor. For any given value, its column represents the axial component of the voltage reading being mapped from, and its row represents the axial component of the estimated weight being mapped to. For example, −180.02 in row 1, column 2 is the scaling factor applied to the Y-component of the voltage reading to get a least-squares regression estimate of the X-component of the load.
X VoltageY VoltageZ Voltage
X load−2903.2−367.68−77.003
Y load201.952918.978.264
Z load−61.526−3.72583943.6
Table 8. Coefficients and R-squared metrics of mixed-quadratic surface fits for each relative residual type.
Table 8. Coefficients and R-squared metrics of mixed-quadratic surface fits for each relative residual type.
Coefficient TypeTotal RelativeX RelativeY RelativeZ Relative
Intercept1.93 × 10 1 6.80 × 10 2 5.32 × 10 15 1.39 × 10 1
q 1 6.52 × 10 2 1.26 × 10 1 7.45 × 10 13 5.59 × 10 1
q 2 1.81 × 10 1 1.73 3.31 × 10 13 9.93 × 10 2
q 1 2 6.35 × 10 5 3.85 × 10 2 1.96 × 10 11 1.75 × 10 3
q 1 × q 2 4.69 × 10 4 7.12 × 10 3 1.21 × 10 11 7.09 × 10 4
q 2 2 2.05 × 10 4 4.00 × 10 3 2.60 × 10 10 1.66 × 10 5
R 2 32.16%13.51%7.14%9.93%
Table 9. Comparative performance of 3-DOF strain sensor to existing sensors used in upper-limb rehabilitation.
Table 9. Comparative performance of 3-DOF strain sensor to existing sensors used in upper-limb rehabilitation.
MetricUnitTest
Sensor
[4][5][27][28]
Sensor CustomFSG15N1AMilliNewtonNano17Custom
DOF 31165
Massg2.63N/A9.075.93
SpanN14.71150.4 to 21260
Accuracy% FSO *2.182.251N/AN/A
Hysteresis% FSO *2.293N/AN/AN/A
Nonlinearity% FSO *2.941.5N/AN/A1
Repeatability% FSO *0.6121N/AN/A
Cost$350118N/A5000N/A
Planar dimmm23 × 16.212.7 × 18.225.4179
Depthmm10.24912.714.51
N/A: Specifications are not publicly known. * FSO: Full scale output.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Carducci, J.; Olds, K.; Krakauer, J.W.; Xu, J.; Brown, J.D. Novel Planar Strain Sensor Design for Capturing 3-Dimensional Fingertip Forces from Patients Affected by Hand Paralysis. Sensors 2022, 22, 7441. https://doi.org/10.3390/s22197441

AMA Style

Carducci J, Olds K, Krakauer JW, Xu J, Brown JD. Novel Planar Strain Sensor Design for Capturing 3-Dimensional Fingertip Forces from Patients Affected by Hand Paralysis. Sensors. 2022; 22(19):7441. https://doi.org/10.3390/s22197441

Chicago/Turabian Style

Carducci, Jacob, Kevin Olds, John W. Krakauer, Jing Xu, and Jeremy D. Brown. 2022. "Novel Planar Strain Sensor Design for Capturing 3-Dimensional Fingertip Forces from Patients Affected by Hand Paralysis" Sensors 22, no. 19: 7441. https://doi.org/10.3390/s22197441

APA Style

Carducci, J., Olds, K., Krakauer, J. W., Xu, J., & Brown, J. D. (2022). Novel Planar Strain Sensor Design for Capturing 3-Dimensional Fingertip Forces from Patients Affected by Hand Paralysis. Sensors, 22(19), 7441. https://doi.org/10.3390/s22197441

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop