Kinematic plasticity during flight in fruit bats: individual variability in response to loading

All bats experience daily and seasonal fluctuation in body mass. An increase in mass requires changes in flight kinematics to produce the extra lift necessary to compensate for increased weight. How bats modify their kinematics to increase lift, however, is not well understood. In this study, we investigated the effect of a 20% increase in mass on flight kinematics for Cynopterus brachyotis, the lesser dog-faced fruit bat. We reconstructed the 3D wing kinematics and how they changed with the additional mass. Bats showed a marked change in wing kinematics in response to loading, but changes varied among individuals. Each bat adjusted a different combination of kinematic parameters to increase lift, indicating that aerodynamic force generation can be modulated in multiple ways. Two main kinematic strategies were distinguished: bats either changed the motion of the wings by primarily increasing wingbeat frequency, or changed the configuration of the wings by increasing wing area and camber. The complex, individual-dependent response to increased loading in our bats points to an underappreciated aspect of locomotor control, in which the inherent complexity of the biomechanical system allows for kinematic plasticity. The kinematic plasticity and functional redundancy observed in bat flight can have evolutionary consequences, such as an increase potential for morphological and kinematic diversification due to weakened locomotor trade-offs.     

Direct measurements of the kinematics and dynamics of bat flight

Experimental measurements and analysis of the flight of bats are presented, including kinematic analysis of high-speed stereo videography of straight and turning flight, and measurements of the wake velocity field behind the bat. The kinematic data reveal that, at relatively slow flight speeds, wing motion is quite complex, including a sharp retraction of the wing during the upstroke and a broad sweep of the partially extended wing during the downstroke. The data also indicate that the flight speed and elevation are not constant, but oscillate in synchrony with both the horizontal and vertical movements of the wing. PIV measurements in the transverse (Trefftz) plane of the wake indicate a complex 'wake vortex' structure dominated by a strong wing tip vortex shed from the wing tip during the downstroke and either the wing tip or a more proximal joint during the upstroke. Data synthesis of several discrete realizations suggests a 'cartoon' of the wake structure during the entire wing beat cycle. Considerable work remains to be done to confirm and amplify these results.     

Parameter optimization model of learning in stepping motion

In this study we combine the representation of motion by a finite number of hardwired functions with parameter optimization to model learning during a stepping motion. Representation of experimental kinematic data by a finite number of predetermined functions and undetermined coefficients was analyzed. Least squares approximation was used to represent experimental data of stepping motions over obstacles of different heights. Functional relationships between coefficients and obstacles heights were also obtained. Learning of stepping over an obstacle was then formulated as a finite dimensional optimization problem. The pattern of foot path, and joint angles trajectories obtained by this learning model, were then compared to the experimental data. The results of the data fitting analysis and of the optimization process as a model for motion learning, indicate that motion can be adequately represented by a set of hardwired functions, and a finite number of task dependent coefficients.     

In vitro analysis of kinematics and elastostatics of the human rib cage during thoracic spinal movement for the validation of numerical models

Neither kinematic nor stiffness properties of the rib cage during thoracic spinal motion were investigated in previous studies, while being essential for the accurate validation of numerical models of the whole thorax. The aim of this in vitro study therefore was to quantify the kinematics and elastostatics of the human rib cage under defined boundary conditions. Eight fresh frozen human thoracic spine specimens (C7-L1, median age 55 years, ranging from 40 to 60 years) including entire rib cages were loaded quasi-statically in flexion/extension, lateral bending, and axial rotation using pure moments of 5 Nm. Relative motions of ribs, thoracic vertebrae, and sternal structures as well as strains on the ribs were measured using optical motion tracking of 150 reflective markers per specimen, while specimens were loaded displacement-controlled with a constant rate of 1°/s for 3.5 cycles. The third full cycle was used to determine relative angles and strains at full loading of the spine for all motion directions. Largest relative angles were found in the main loading directions with only small motions at the mid-thoracic levels. Highest strains of the intercostal spaces were detected in the anterior section of the lowest fourth of the rib cage, showing compressions and elongations of more than 10% in all spinal motion planes. Elastostatic rib deformation was generally less than 1%. Rib-sternum relative motions exhibited complex motion patterns, overall showing relative angles below 2°. The results indicate that rib cage structures are not macroscopically deformed during spinal motion, but exhibit characteristic reproducible kinematics patterns.     

Minimum toe clearance adaptations to floor surface irregularity and gait speed

Toe speed during gait generally nears its maximum while its height reaches a local minima approximately halfway through swing phase. Trips are thought to frequently occur at these local minima (minimum toe clearance or MTC events) and trip risk has been quantified using the minimum distance between the toe and ground here (MTC). This study investigated MTC on floor surfaces with and without multiple small obstacles. After shoes and floor surfaces were digitized, 14 unimpaired subjects (half women) each traversed a 4.88 m walkway 4 times at slow, preferred, and fast speeds across surfaces with no obstacles, visible obstacles, and hidden obstacles. Both surfaces with obstacles had the same random obstacle configuration. Shoe and body segment motions were tracked using passive markers and MTC and joint kinematics calculated. All MTC and kinematic variables tested significantly increased with faster instructed gait speed except the likelihood of MTC event occurrence (local minima in minimum toe clearance trajectory when foot is in upper quartile of speed). MTC events were less frequent for swing phases on surfaces with obstacles (80% vs. 98% for no obstacles). MTC values, when present, were doubled by the presence of visible obstacles (22.2 ± 7.3mm vs. 11.1 ± 5.7 mm) and further increased to 26.8 ± 7.1mm when these obstacles were hidden from view (all comparisons p ≤ 0.0003). These substantial floor surface-related changes in MTC event occurrences and values resulted from alterations in toe- and heel-clearance trajectories caused by subtle but significant changes in joint kinematics that did not exceed 10% each joint's swing phase range of motion.     

Kinematics of lower limbs during walking are emulated by springy walking model with a compliantly connected,off-centered curvy foot

The dynamics of the center of mass (CoM) during walking and running at various gait conditions are well described by the mechanics of a simple passive spring loaded inverted pendulum (SLIP). Due to its simplicity, however, the current form of the SLIP model is limited at providing any further information about multi-segmental lower limbs that generate oscillatory CoM behaviors and their corresponding ground reaction forces. Considering that the dynamics of the CoM are simply achieved by mass-spring mechanics, we wondered whether any of the multi-joint motions could be demonstrated by simple mechanics. In this study, we expand a SLIP model of human locomotion with an off-centered curvy foot connected to the leg by a springy segment that emulates the asymmetric kinematics and kinetics of the ankle joint. The passive dynamics of the proposed expansion of the SLIP model demonstrated the empirical data of ground reaction forces, center of mass trajectories, ankle joint kinematics and corresponding ankle joint torque at various gait speeds. From the mechanically simulated trajectories of the ankle joint and CoM, the motion of lower-limb segments, such as thigh and shank angles, could be estimated from inverse kinematics. The estimation of lower limb kinematics showed a qualitative match with empirical data of walking at various speeds. The representability of passive compliant mechanics for the kinetics of the CoM and ankle joint and lower limb joint kinematics implies that the coordination of multi-joint lower limbs during gait can be understood with a mechanical framework.     

Validation of a model-based inverse kinematics approach based on wearable inertial sensors

Abstract

Wearable inertial measurement units (IMUs) are a promising solution to human motion estimation. Using IMUs 3D orientations, a model-driven inverse kinematics methodology to estimate joint angles is presented. Estimated joint angles were validated against encoder-measured kinematics (robot) and against marker-based kinematics (passive mechanism). Results are promising, with RMS angular errors respectively lower than 3 and 6?deg over a minimum range of motion of 50?deg (robot) and 160?deg (passive mechanism). Moreover, a noise robustness analysis revealed that the model-driven approach reduces the effects of experimental noises, making the proposed technique particularly suitable for application in human motion analysis.     

Use of dual Euler angles to quantify the three-dimensional joint motion and its application to the ankle joint complex

This paper presents a modified Euler angles method, dual Euler angles approach, to describe general spatial human joint motions. In dual Euler angles approach, the three-dimensional joint motion is considered as three successive screw motions with respect to the axes of the moving segment coordinate system; accordingly, the screw motion displacements are represented by dual Euler angles. The algorithm for calculating dual Euler angles from coordinates of markers on the moving segment is also provided in this study. As an example, the proposed method is applied to describe motions of ankle joint complex during dorsiflexion–plantarflexion. A Flock of Birds electromagnetic tracking device (FOB) was used to measure joint motion in vivo. Preliminary accuracy tests on a gimbal structure demonstrate that the mean errors of dual Euler angles evaluated by using source data from FOB are less than 1° for rotations and 1 mm for translations, respectively. Based on the pilot study, FOB is feasible for quantifying human joint motions using dual Euler angles approach.     

In situ comparison of A-mode ultrasound tracking system and skin-mounted markers for measuring kinematics of the lower extremity

Skin-mounted marker based motion capture systems are widely used in measuring the movement of human joints. Kinematic measurements associated with skin-mounted markers are subject to soft tissue artifacts (STA), since the markers follow skin movement, thus generating errors when used to represent motions of underlying bone segments. We present a novel ultrasound tracking system that is capable of directly measuring tibial and femoral bone surfaces during dynamic motions, and subsequently measuring six-degree-of-freedom (6-DOF) tibiofemoral kinematics. The aim of this study is to quantitatively compare the accuracy of tibiofemoral kinematics estimated by the ultrasound tracking system and by a conventional skin-mounted marker based motion capture system in a cadaveric experimental scenario. Two typical tibiofemoral joint models (spherical and hinge models) were used to derive relevant kinematic outcomes. Intra-cortical bone pins equipped with optical markers were inserted in the tibial and femoral bones to serve as a reference to provide ground truth kinematics. The ultrasound tracking system resulted in lower kinematic errors than the skin-mounted markers (the ultrasound tracking system: maximum root-mean-square (RMS) error 3.44° for rotations and 4.88 mm for translations, skin-mounted markers with the spherical joint model: 6.32° and 6.26 mm, the hinge model: 6.38° and 6.52 mm). Our proposed ultrasound tracking system has the potential of measuring direct bone kinematics, thereby mitigating the influence and propagation of STA. Consequently, this technique could be considered as an alternative method for measuring 6-DOF tibiofemoral kinematics, which may be adopted in gait analysis and clinical practice.     

Mandibular kinematics represented by a non-orthogonal floating axis joint coordinate system

There are many methods used to represent joint kinematics (e.g., roll, pitch, and yaw angles; instantaneous center of rotation; kinematic center; helical axis). Often in biomechanics internal landmarks are inferred from external landmarks. This study represents mandibular kinematics using a non-orthogonal floating axis joint coordinate system based on 3-D geometric models with parameters that are "clinician friendly" and mathematically rigorous. Kinematics data for two controls were acquired from passive fiducial markers attached to a custom dental clutch. The geometric models were constructed from MRI data. The superior point along the arc of the long axis of the condyle was used to define the coordinate axes. The kinematic data and geometric models were registered through fiducial markers visible during both protocols. The mean absolute maxima across the subjects for sagittal rotation, coronal rotation, axial rotation, medial-lateral translation, anterior-posterior translation, and inferior-superior translation were 34.10 degrees, 1.82 degrees, 1.14 degrees, 2.31, 21.07, and 6.95 mm, respectively. All the parameters, except for one subject's axial rotation, were reproducible across two motion recording sessions. There was a linear correlation between sagittal rotation and translation, the dominant motion plane, with approximately 1.5 degrees of rotation per millimeter of translation. The novel approach of combining the floating axis system with geometric models succinctly described mandibular kinematics with reproducible and clinician friendly parameters.     

Feasibility of using a video-based motion analysis system for measuring thumb kinematics

While several different methods have been used to measure hand kinematics, fluoroscopy is generally considered to be the most accurate. Recently, video-based motion analysis has been developed for the measurement of joint kinematics. This method is versatile, easy to use, and can measure motions dynamically. Surface markers are most commonly used in the video-based motion systems. However, whether the surface markers placed on the thumb accurately represent the true kinematics of the underlying bony segment is questionable.In this study, the feasibility of surface markers to represent thumb kinematics was investigated by fluoroscopy. Both the positions of surface markers and bony landmarks were simultaneous recorded and then digitized. The Ra(2) values comparing the angular changes of the thumb interphalangeal, metacarpal and carpometacarpal joints derived using the surface markers or bony landmarks were 0.9986, 0.9730 and 0.9186 in the flexion/extension plane respectively, 0.8837, 0.9697 and 0.8775 in the abduction/adduction plane; and 0.9884, 0.9643 and 0.9431 in the opposition plane. The ranges, mean and standard deviation of the absolute differences between calculated angles of different marker sets were also compared. These data revealed that the similarities of the two different marker techniques throughout the motion cycle were high. The differences between the two methods were also within clinically allowable range of +/-5 degrees. It is concluded that the application of the video-based motion analysis system with surface markers to thumb kinematics is warranted.     

Nonlinear smooth orthogonal decomposition of kinematic features of sawing reconstructs muscle fatigue evolution as indicated by electromyography

Tracking or predicting physiological fatigue is important for developing more robust training protocols and better energy supplements and/or reducing muscle injuries. Current methodologies are usually impractical and/or invasive and may not be realizable outside of laboratory settings. It was recently demonstrated that smooth orthogonal decomposition (SOD) of phase space warping (PSW) features of motion kinematics can identify fatigue in individual muscle groups. We hypothesize that a nonlinear extension of SOD will identify more optimal fatigue coordinates and provide a lower-dimensional reconstruction of local fatigue dynamics than the linear SOD. Both linear and nonlinear SODs were applied to PSW features estimated from measured kinematics to reconstruct muscle fatigue dynamics in subjects performing a sawing motion. Ten healthy young right-handed subjects pushed a weighted handle back and forth until voluntary exhaustion. Three sets of joint kinematic angles were measured from the right upper extremity in addition to surface electromyography (EMG) recordings. The SOD coordinates of kinematic PSW features were compared against independently measured fatigue markers (i.e., mean and median EMG spectrum frequencies of individual muscle groups). This comparison was based on a least-squares linear fit of a fixed number of the dominant SOD coordinates to the appropriate local fatigue markers. Between subject variability showed that at most four to five nonlinear SOD coordinates were needed to reconstruct fatigue in local muscle groups, while on average 15 coordinates were needed for the linear SOD. Thus, the nonlinear coordinates provided a one-order-of-magnitude improvement over the linear ones.     

Segmental motion of forefoot and hindfoot as a diagnostic tool

Segmental motions derived from non-invasive motion analysis are being used to investigate the intrinsic functional behavior of the foot and ankle in health and disease. The goal of this research was to examine the ability of a generic segmented model of the foot to capture and differentiate changes in internal skeletal kinematics due to neuromuscular disease and/or trauma. A robotic apparatus that reproduces the kinematics and kinetics of gait in cadaver lower extremities was employed to produce motion under normal and aberrant neuromuscular activation patterns of tibialis posterior and/or tibialis anterior. Stance phase simulations were conducted on 10 donor limbs while recording three-dimensional kinematic trajectories of (1) skin-mounted markers used clinically to construct segmented foot models, and (2) bone-mounted marker clusters to capture actual internal bone motion as the gold standard for comparison. The models constructed from external marker data were able to differentiate the kinematic behaviors elicited by different neuromuscular conditions in a manner similar to that using the bone-derived data. Measurable differences between internal and externally measured kinematics were small, variable and random across the three axes of rotation and neuromuscular conditions, with a tendency toward more differences noted during early and late stance. Albeit slightly different, three-dimensional motion profiles of the hindfoot and forefoot segments correlated well with internal skeletal motion under all neuromuscular conditions, thereby confirming the utility of measuring segmental motions as a valid means of clinical assessment.     

Hindlimb Motion during Steady Flight of the Lesser Dog-Faced Fruit Bat,Cynopterus brachyotis

In bats, the wing membrane is anchored not only to the body and forelimb, but also to the hindlimb. This attachment configuration gives bats the potential to modulate wing shape by moving the hindlimb, such as by joint movement at the hip or knee. Such movements could modulate lift, drag, or the pitching moment. In this study we address: 1) how the ankle translates through space during the wingbeat cycle; 2) whether amplitude of ankle motion is dependent upon flight speed; 3) how tension in the wing membrane pulls the ankle; and 4) whether wing membrane tension is responsible for driving ankle motion. We flew five individuals of the lesser dog-faced fruit bat, Cynopterus brachyotis (Family: Pteropodidae), in a wind tunnel and documented kinematics of the forelimb, hip, ankle, and trailing edge of the wing membrane. Based on kinematic analysis of hindlimb and forelimb movements, we found that: 1) during downstroke, the ankle moved ventrally and during upstroke the ankle moved dorsally; 2) there was considerable variation in amplitude of ankle motion, but amplitude did not correlate significantly with flight speed; 3) during downstroke, tension generated by the wing membrane acted to pull the ankle dorsally, and during upstroke, the wing membrane pulled laterally when taut and dorsally when relatively slack; and 4) wing membrane tension generally opposed dorsoventral ankle motion. We conclude that during forward flight in C. brachyotis, wing membrane tension does not power hindlimb motion; instead, we propose that hindlimb movements arise from muscle activity and/or inertial effects.     

A mechanical supination sprain simulator for studying ankle supination sprain kinematics

This study presents a free-fall mechanical supination sprain simulator for evaluating the ankle joint kinematics during a simulated ankle supination sprain injury. The device allows the foot to be in an anatomical position before the sudden motion, and also allows different degrees of supination, or a combination of inversion and plantarflexion. Five subjects performed simulated supination sprain trials in five different supination angles. Ankle motion was captured by a motion analysis system, and the ankle kinematics were reported in plantarflexion/dorsiflexion, inversion/eversion and internal/external rotation planes. Results showed that all sprain motions were not pure single-plane motions but were accompanied by motion in other two planes, therefore, different degrees of supination were achieved. The presented sprain simulator allows a more comprehensive study of the kinematics of ankle sprain when compared with some previous laboratory research designs.     

The efficacy of a video-based marker-less tracking system for gait analysis

Most clinical gait analyses are conducted using motion capture systems which track retro-reflective markers that are placed on key landmarks of the participants. An alternative to a three-dimensional (3D) motion capture, marker-based, optical camera system may be a marker-less video-based tracking system. The aim of our study was to investigate the efficacy of the use of a marker-less tracking system in the calculation of 3D joint angles for possible use in clinical gait analysis. Ten participants walked and jogged on a treadmill and their kinematic data were captured with a marker and marker-less tracking system simultaneously. The hip, knee and ankle angles in the frontal, sagittal and transverse planes were computed. Root Mean Square differences (RMS_diff) between corresponding angles for each participant’s support phase were calculated and averaged to derive the mean within-subject RMS_diff. These within-subject means were averaged to obtain the mean between-subject RMS_diff for the relevant joint angles in the two gait conditions (walking and jogging). The RMS_diff between the two tracking systems was less than 1° for all rotations of the three joint angles of the hip and knee. However, there were slightly larger differences in the ankle joint angles. The results of this study suggest a potential application in gait analysis in clinical settings where observations of anatomical motions may provide meaningful feedback.     

The effects of anatomical errors on shoulder kinematics computed using multi-body models

Joint motion calculated using multi-body models and inverse kinematics presents many advantages over direct marker-based calculations. However, the sensitivity of the computed kinematics is known to be partly caused by the model and could also be influenced by the participants’ anthropometry and sex. This study aimed to compare kinematics computed from an anatomical shoulder model based on medical images against a scaled-generic model and quantify the effects of anatomical errors and participants’ anthropometry on the calculated joint angles. Twelve participants have had planar shoulder movements experimentally captured in a motion lab, and their shoulder anatomy imaged using an MRI scanner. A shoulder multi-body dynamics model was developed for each participant, using both an image-based approach and a scaled-generic approach. Inverse kinematics have been performed using the two different modelling procedures and the three different experimental motions. Results have been compared using Bland–Altman analysis of agreement and further analysed using multi-linear regressions. Kinematics computed via an anatomical and a scaled-generic shoulder models differed in average from 3.2 to 5.4 degrees depending on the task. The MRI-based model presented smaller limits of agreement to direct kinematics than the scaled-generic model. Finally, the regression model predictors, including anatomical errors, sex, and BMI of the participant, explained from 41 to 80% of the kinematic variability between model types with respect to the task. This study highlighted the consequences of modelling precision, quantified the effects of anatomical errors on the shoulder kinematics, and showed that participants' anthropometry and sex could indirectly affect kinematic outcomes.

     

Comparative assessment of bone pose estimation using Point Cluster Technique and OpenSim

Estimating the position of the bones from optical motion capture data is a challenge associated with human movement analysis. Bone pose estimation techniques such as the Point Cluster Technique (PCT) and simulations of movement through software packages such as OpenSim are used to minimize soft tissue artifact and estimate skeletal position; however, using different methods for analysis may produce differing kinematic results which could lead to differences in clinical interpretation such as a misclassification of normal or pathological gait. This study evaluated the differences present in knee joint kinematics as a result of calculating joint angles using various techniques. We calculated knee joint kinematics from experimental gait data using the standard PCT, the least squares approach in OpenSim applied to experimental marker data, and the least squares approach in OpenSim applied to the results of the PCT algorithm. Maximum and resultant RMS differences in knee angles were calculated between all techniques. We observed differences in flexion/extension, varus/valgus, and internal/external rotation angles between all approaches. The largest differences were between the PCT results and all results calculated using OpenSim. The RMS differences averaged nearly 5° for flexion/extension angles with maximum differences exceeding 15°. Average RMS differences were relatively small (< 1.08°) between results calculated within OpenSim, suggesting that the choice of marker weighting is not critical to the results of the least squares inverse kinematics calculations. The largest difference between techniques appeared to be a constant offset between the PCT and all OpenSim results, which may be due to differences in the definition of anatomical reference frames, scaling of musculoskeletal models, and/or placement of virtual markers within OpenSim. Different methods for data analysis can produce largely different kinematic results, which could lead to the misclassification of normal or pathological gait. Improved techniques to allow non-uniform scaling of generic models to more accurately reflect subject-specific bone geometries and anatomical reference frames may reduce differences between bone pose estimation techniques and allow for comparison across gait analysis platforms.     

Kinematic performance of a six degree-of-freedom hand model (6DHand) for use in occupational biomechanics

Upper extremity musculoskeletal disorders represent an important health issue across all industry sectors; as such, the need exists to develop models of the hand that provide comprehensive biomechanics during occupational tasks. Previous optical motion capture studies used a single marker on the dorsal aspect of finger joints, allowing calculation of one and two degree-of-freedom (DOF) joint angles; additional algorithms were needed to define joint centers and the palmar surface of fingers. We developed a 6DOF model (6DHand) to obtain unconstrained kinematics of finger segments, modeled as frusta of right circular cones that approximate the palmar surface. To evaluate kinematic performance, twenty subjects gripped a cylindrical handle as a surrogate for a powered hand tool. We hypothesized that accessory motions (metacarpophalangeal pronation/supination; proximal and distal interphalangeal radial/ulnar deviation and pronation/supination; all joint translations) would be small (less than 5° rotations, less than 2mm translations) if segment anatomical reference frames were aligned correctly, and skin movement artifacts were negligible. For the gripping task, 93 of 112 accessory motions were small by our definition, suggesting this 6DOF approach appropriately models joints of the fingers. Metacarpophalangeal supination was larger than expected (approximately 10°), and may be adjusted through local reference frame optimization procedures previously developed for knee kinematics in gait analysis. Proximal translations at the metacarpophalangeal joints (approximately 10mm) were explained by skin movement across the metacarpals, but would not corrupt inverse dynamics calculated for the phalanges. We assessed performance in this study; a more rigorous validation would likely require medical imaging.     

Kinematic Analysis and Experimental Verification on the Locomotion of Gecko

This paper presents a kinematic analysis of the locomotion of a gecko,and experimental verification of the kinematicmodel.Kinematic analysis is important for parameter design,dynamic analysis,and optimization in biomimetic robot research.The proposed kinematic analysis can simulate,without iteration,the locomotion of gecko satisfying the constraint conditionsthat maintain the position of the contacted feet on the surface.So the method has an advantage for analyzing the climbing motionof the quadruped mechanism in a real time application.The kinematic model of a gecko consists of four legs based on 7-degreesof freedom spherical-revolute-spherical joints and two revolute joints in the waist.The motion of the kinematic model issimulated based on measurement data of each joint.The motion of the kinematic model simulates the investigated real gecko’smotion by using the experimental results.The analysis solves the forward kinematics by considering the model as a combinationof closed and open serial mechanisms under the condition that maintains the contact positions of the attached feet on the ground.The motions of each joint are validated by comparing with the experimental results.In addition to the measured gait,three othergaits are simulated based on the kinematic model.The maximum strides of each gait are calculated by workspace analysis.Theresult can be used in biomimetic robot design and motion planning.