skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Joint velocity dependence of fatigue in isokinetic tasks
The ability to predict the decline in muscle strength over the course of an activity (i.e., fatigue) can be a crucial aid to task design, injury prevention, and rehabilitation efforts. Current models of muscle fatigue have been hitherto validated only for isometric contractions, but most real-world tasks are dynamic in nature, involving continuously varying joint velocities. It has previously been proposed that a three-compartment-controller (3CC) model might be used to predict fatigue for such tasks by using it in conjunction with joint- and direction-specific torque-velocity-angle (TVA) surfaces. This allows for the calculation of a time-varying target load parameter that can be used by the 3CC model, but it increases model complexity and has not been validated by experimental data. An alternative approach is proposed where the effect of joint velocity is modeled by a velocity parameter and integrated into the fatigue model equations, removing the dependence on external TVA surfaces. The predictions using both methods are contrasted against experimental data collected from 20 subjects in a series of isokinetic tests involving the knee and shoulder joints, covering a range of velocities encountered in day-to-day tasks. A much lower degree of fatigue is observed for moderate velocities compared to that for very low or very high velocities. Predictions using the integrated velocity parameter are computationally less expensive than using TVA surfaces and are also closer to experimentally obtained values. The modified fatigue model can therefore be applied to dynamic tasks with varying velocities when the task is discretized into several isokinetic tasks.  more » « less
Award ID(s):
2014281 1849279
PAR ID:
10428809
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Proceedings of the 7th International Digital Human Modeling Symposium
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract In this study, a hybrid predictive model is used to predict 3D asymmetric lifting motion and assess potential musculoskeletal lower back injuries for asymmetric lifting tasks. The hybrid model has two modules: a skeletal module and an OpenSim musculoskeletal module. The skeletal module consists of a dynamic joint strength based 40 degrees of freedom spatial skeletal model. The skeletal module can predict the lifting motion, ground reaction forces (GRFs), and center of pressure (COP) trajectory using an inverse dynamics based optimization method. The equations of motion are built by recursive Lagrangian dynamics. The musculoskeletal module consists of a 324-muscle-actuated full-body lumbar spine model. Based on the generated kinematics, GRFs and COP data from the skeletal module, the musculoskeletal module estimates muscle activations using static optimization and joint reaction forces through the joint reaction analysis tool. Muscle activation results between simulated and experimental EMG are compared to validate the model. Finally, potential lower back injuries are evaluated for a specific-weight asymmetric lifting task. The shear and compression spine loads are compared to NIOSH recommended limits. At the beginning of the dynamic lifting process, the simulated compressive spine load beyond the NIOSH action limit but less than the permissible limit. This is due to the fatigue factors considered in NIOSH lifting equation. 
    more » « less
  2. ABSTRACT Recent studies have demonstrated that muscle force is not determined solely by activation under dynamic conditions, and that length history has an important role in determining dynamic muscle force. Yet, the mechanisms for how muscle force is produced under dynamic conditions remain unclear. To explore this, we investigated the effects of muscle stiffness, activation and length perturbations on muscle force. First, submaximal isometric contraction was established for whole soleus muscles. Next, the muscles were actively shortened at three velocities. During active shortening, we measured muscle stiffness at optimal muscle length (L0) and the force response to time-varying activation and length perturbations. We found that muscle stiffness increased with activation but decreased as shortening velocity increased. The slope of the relationship between maximum force and activation amplitude differed significantly among shortening velocities. Also, the intercept and slope of the relationship between length perturbation amplitude and maximum force decreased with shortening velocity. As shortening velocities were related to muscle stiffness, the results suggest that length history determines muscle stiffness and the history-dependent muscle stiffness influences the contribution of activation and length perturbations to muscle force. A two-parameter viscoelastic model including a linear spring and a linear damper in parallel with measured stiffness predicted history-dependent muscle force with high accuracy. The results and simulations support the hypothesis that muscle force under dynamic conditions can be accurately predicted as the force response of a history-dependent viscoelastic material to length perturbations. 
    more » « less
  3. In this study, a 3D asymmetric lifting motion is predicted by using a hybrid predictive model to prevent potential musculoskeletal lower back injuries for asymmetric lifting tasks. The hybrid model has two modules: a skeletal module and an OpenSim musculoskeletal module. The skeletal module consists of a dynamic joint strength based 40 degrees of freedom spatial skeletal model. The skeletal module can predict the lifting motion, ground reaction forces (GRFs), and center of pressure (COP) trajectory using an inverse dynamics-based motion optimization method. The musculoskeletal module consists of a 324-muscle-actuated full-body lumbar spine model. Based on the predicted kinematics, GRFs and COP data from the skeletal module, the musculoskeletal module estimates muscle activations using static optimization and joint reaction forces through the joint reaction analysis tool in OpenSim. The predicted asymmetric motion and GRFs are validated with experimental data. Muscle activation results between the simulated and experimental EMG are also compared to validate the model. Finally, the shear and compression spine loads are compared to NIOSH recommended limits. The differences between asymmetric and symmetric liftings are also compared. 
    more » « less
  4. Physics-based simulations of Arctic sea ice are highly complex, involving transport between different phases, length scales, and time scales. Resultantly, numerical simulations of sea ice dynamics have a high computational cost and model uncertainty. We employ data-driven machine learning (ML) to make predictions of sea ice motion. The ML models are built to predict present-day sea ice velocity given present-day wind velocity and previous-day sea ice concentration and velocity. Models are trained using reanalysis winds and satellite-derived sea ice properties. We compare the predictions of three different models: persistence (PS), linear regression (LR), and a convolutional neural network (CNN). We quantify the spatiotemporal variability of the correlation between observations and the statistical model predictions. Additionally, we analyze model performance in comparison to variability in properties related to ice motion (wind velocity, ice velocity, ice concentration, distance from coast, bathymetric depth) to understand the processes related to decreases in model performance. Results indicate that a CNN makes skillful predictions of daily sea ice velocity with a correlation up to 0.81 between predicted and observed sea ice velocity, while the LR and PS implementations exhibit correlations of 0.78 and 0.69, respectively. The correlation varies spatially and seasonally: lower values occur in shallow coastal regions and during times of minimum sea ice extent. LR parameter analysis indicates that wind velocity plays the largest role in predicting sea ice velocity on 1-day time scales, particularly in the central Arctic. Regions where wind velocity has the largest LR parameter are regions where the CNN has higher predictive skill than the LR. 
    more » « less
  5. Abstract In this study, the fatigue progression and optimal motion trajectory during repetitive lifting task is predicted by using a 10 degrees of freedom (DOFs) two-dimensional (2D) digital human model and a three-compartment controller (3CC) fatigue model. The numerical analysis is further validated by conducting an experiment under similar conditions. The human is modeled using Denavit-Hartenberg (DH) representation. The task is mathematically formulated as a nonlinear optimization problem where the dynamic effort of the joints is minimized subjected to physical and task specific constraints. A sequential quadratic programming method is used for the optimization process. The design variables include control points of (1) quartic B-splines of the joint angle profiles; and (2) the three compartment sizes profiles for the six physical joints of interest — spine, shoulder, elbow, hip, knee, and ankle. Both numerical and experimental liftings are performed with a 15.2 kg box as external load. The simulation reports the human joint torque profiles and the progression of joint fatigue. The joint torque profiles show periodic trends. A maximum of 17 cycles are predicted before the repetitive lifting task fails, which also reasonably agrees with that of the experimental results (16 cycles). This formulation is also a generalized one, hence it can be used for other repetitive motion studies as well. 
    more » « less