A Computationally Efficient Muscle Model

Author(s):  
Matthew Millard ◽  
Scott Delp

Modern biomechanical models often include many degrees of freedom, complex joint geometry, and hundreds of muscles [1]. Such complex models may provide realistic representations of musculoskeletal dynamics, but they require a great deal of time to simulate. The state equations of the muscles are frequently the stiffest numerical element in a biomechanical model and determine the simulation time. Here we introduce a muscle model that generates similar forces as a conventional muscle model, but requires fewer explicit integration steps.

2018 ◽  
Vol 63 (1) ◽  
pp. 7-15 ◽  
Author(s):  
László Bencsik ◽  
Ambrus Zelei

Biomechanical models of different complexity are used to understand the dynamics of human running. Low degrees-of-freedom models are appropriate for the prediction of the effect of certain parameter changes. We present a minimally complex biomechanical model which characterizes the effects of foot strike pattern and shank angle on the ground-foot impact intensity, which influences the risk of injuries and energy efficiency.A three segment leg model (thigh, shank and foot) is proposed combined with the mass of the rest of the body parts concentrated in the hip. The ground-foot impact intensity and the absorbed kinetic energy are analyzed using multibody dynamics tools. The impact intensity was discovered in the parameter space of the angle of the thigh, the angle of the shank, the foot strike pattern and the running speed.The results regarding the effect of strike pattern are in coincidence with the literature: forefoot strike implies lower impact intensity and energy absorption than rearfoot strike. However, in contrast of the previous result of a two segment foot model from the related literature, the calculations indicated that the shank angle highly affects the impact intensity: the impact intensity can be reduced by foot touchdown under the hip. We showed that foot and shank cannot be analyzed in itself without considering the thigh and the total body weight, and we also confirmed that the horizontal velocity cannot be neglected when foot impact is analyzed.


2021 ◽  
Vol 82 (5) ◽  
Author(s):  
Hannah J. Pybus ◽  
Amanda L. Tatler ◽  
Lowell T. Edgar ◽  
Reuben D. O’Dea ◽  
Bindi S. Brook

AbstractPrecision-cut lung-slices (PCLS), in which viable airways embedded within lung parenchyma are stretched or induced to contract, are a widely used ex vivo assay to investigate bronchoconstriction and, more recently, mechanical activation of pro-remodelling cytokines in asthmatic airways. We develop a nonlinear fibre-reinforced biomechanical model accounting for smooth muscle contraction and extracellular matrix strain-stiffening. Through numerical simulation, we describe the stresses and contractile responses of an airway within a PCLS of finite thickness, exposing the importance of smooth muscle contraction on the local stress state within the airway. We then consider two simplifying limits of the model (a membrane representation and an asymptotic reduction in the thin-PCLS-limit), that permit analytical progress. Comparison against numerical solution of the full problem shows that the asymptotic reduction successfully captures the key elements of the full model behaviour. The more tractable reduced model that we develop is suitable to be employed in investigations to elucidate the time-dependent feedback mechanisms linking airway mechanics and cytokine activation in asthma.


Author(s):  
Zachary Merrill ◽  
April Chambers ◽  
Rakié Cham

Body segment parameters (BSPs) such as segment mass and center of mass are used as inputs in ergonomic design and biomechanical models to predict the risk of musculoskeletal injuries. These models have been shown to be sensitive to the BSP values used as inputs, demonstrating the necessity of using accurate and representative parameters. This study aims to provide accurate BSPs by quantifying the impact of age and body mass index on torso and thigh mass and center of mass in working adults using whole body dual energy x-ray absorptiometry (DXA) scan data. The results showed significant effects of gender, age, and body mass index (BMI) on torso and thigh mass and center of mass, as well as significant effects of age and BMI within genders, indicating that age, gender, and BMI need to be taken into account when predicting BSPs in order to calculate representative ergonomic and biomechanical model outputs.


Author(s):  
Derek Lura ◽  
Rajiv Dubey ◽  
Stephanie L. Carey ◽  
M. Jason Highsmith

The prostheses used by the majority of persons with hand/arm amputations today have a very limited range of motion. Transradial (below the elbow) amputees lose the three degrees of freedom provided by the wrist and forearm. Some myoeletric prostheses currently allow for forearm pronation and supination (rotation about an axis parallel to the forearm) and the operation of a powered prosthetic hand. Older body-powered prostheses, incorporating hooks and other cable driven terminal devices, have even fewer degrees of freedom. In order to perform activities of daily living (ADL), a person with amputation(s) must use a greater than normal range of movement from other body joints to compensate for the loss of movement caused by the amputation. By studying the compensatory motion of prosthetic users we can understand the mechanics of how they adapt to the loss of range of motion in a given limb for select tasks. The purpose of this study is to create a biomechanical model that can predict the compensatory motion using given subject data. The simulation can then be used to select the best prosthesis for a given user, or to design prostheses that are more effective at selected tasks, once enough data has been analyzed. Joint locations necessary to accomplish the task with a given configuration are calculated by the simulation for a set of prostheses and tasks. The simulation contains a set of prosthetic configurations that are represented by parameters that consist of the degrees of freedom provided by the selected prosthesis. The simulation also contains a set of task information that includes joint constraints, and trajectories which the hand or prosthesis follows to perform the task. The simulation allows for movement in the wrist and forearm, which is dependent on the prosthetic configuration, elbow flexion, three degrees of rotation at the shoulder joint, movement of the shoulder joint about the sternoclavicular joint, and translation and rotation of the torso. All joints have definable restrictions determined by the prosthesis, and task.


Author(s):  
Mortadha Graa ◽  
Mohamed Nejlaoui ◽  
Ajmi Houidi ◽  
Zouhaier Affi ◽  
Lotfi Romdhane

In this paper, an analytical reduced dynamic model of a rail vehicle system is developed. This model considers only 38 degrees of freedom of the rail vehicle system. This reduced model can predict the dynamic behaviour of the rail vehicle while being simpler than existing dynamic models. The developed model is validated using experimental results found in the bibliography and its results are compared with existing more complex models from the literature. The developed model is used for the passenger comfort evaluation, which is based on the value of the weighted root mean square acceleration according to the ISO 2631 standard. Several parameters of the system, i.e., passenger position, loading of the railway vehicle and its speed, and their effect on the passenger comfort are investigated. It was shown that the level of comfort is mostly affected by the speed of the railway vehicle and the position of the seat. The load, however, did not have a significant effect on the level of comfort of the passenger.


Author(s):  
Raj Desai ◽  
Anirban Guha ◽  
Pasumarthy Seshu

Long duration automobile-induced vibration is the cause of many ailments to humans. Predicting and mitigating these vibrations through seat requires a good model of seated human body. A good model is the one that strikes the right balance between modelling difficulty and simulation results accuracy. Increasing the number of body parts which have been separately modelled and increasing the number of ways these parts are connected to each other increase the number of degrees of freedom of the entire model. A number of such models have been reported in the literature. These range from simple lumped parameter models with limited accuracy to advanced models with high computational cost. However, a systematic comparison of these models has not been reported till date. This work creates eight such models ranging from 8 to 26 degrees of freedom and tries to identify the model which strikes the right balance between modelling complexity and results accuracy. A comparison of the models’ prediction with experimental data published in the literature allows the identification of a 12 degree of freedom backrest supported model as optimum for modelling complexity and prediction accuracy.


2009 ◽  
Vol 12 (01) ◽  
pp. 31-43 ◽  
Author(s):  
Rositsa T. Raikova

Less attention is paid to joint reactions when optimization tasks are solved aiming to predict individual muscle forces driving a biomechanical model. The reactions are important, however, for joint stability and for prevention from injuries, especially for fast motions and submaximal loading. The purpose of the paper is to investigate the influence of the joint reaction as a criterion in an objective function and to study the possibilities for prediction of antagonistic co-contraction. Planar elbow flexions in the sagittal plane with duration from 0.4 to 2 s are simulated, and muscle forces and elbow joint reaction are calculated solving numerically optimization tasks formulated for models with one (elbow moment equation only) and two (elbow and shoulder moment equations) degrees of freedom (DOF). The objective function is a weighted sum of muscle forces and joint reaction raised to different powers. The following conclusions can be made: (1) if the joint reaction is included in the objective function, antagonistic co-contraction can be predicted even for 1 DOF model; in some situations the use of such objective function can destroy the synergistic muscles' action; (2) the prediction of antagonistic muscles' co-contraction for 2 DOF model depends on the way the biarticular muscles are modeled, and this is valid for both dynamic and quasistatic conditions; if there are no biarticular muscles, antagonistic co-contraction cannot be predicted in one joint using popular objective functions, like minimum of sum of muscle forces or muscle stresses raised to a power.


Author(s):  
Fernando Geremias Toni ◽  
Clóvis de Arruda Martins

Abstract Due to the number of layers and their respective geometrical complexities, finite element analyzes of flexible pipes usually require large-scale schemes, with a high number of elements and degrees-of-freedom. If proper precautions are not taken, such as suitable algorithms and numerical methods, the computational costs of these analyzes may become unfeasible to the current computational standards. Finite macroelements are finite elements formulated for the solution of a specific problem considering and taking advantage of its particularities, such as geometry patterns, in order to obtain computational advantages, as reduced number of degrees-of-freedom and ease of problem description. The element-by-element method (EBE) also fits very well in this context, since it is characterized by the elimination of the global stiffness matrix and its memory consumption grows linearly with the number of elements, besides being highly parallelizable. Over the last decades, several works regarding the EBE method were published in the literature, but none of them directly applied to flexible pipes. Due to the contact elements between the layers, problems with flexible pipes are usually characterized by very large matrix-bandwidth, making the implementation of EBE method more challenging, so that its efficiency and scalability are not compromised. Therefore, this work presents a parallelized implementation of an element-by-element architecture for structural analysis of flexible pipes using finite macroelements, consisting of an extension of a previous work from the same authors. New synchronization algorithms were developed, with scalability improvements, the methodology was extended to other finite macroelements and comparisons were made with a well-stablished FEM software, with significant gains in simulation time and memory consumption.


2020 ◽  
Vol 39 (10-11) ◽  
pp. 1239-1258
Author(s):  
Shameek Ganguly ◽  
Oussama Khatib

Multi-surface interactions occur frequently in articulated-rigid-body systems such as robotic manipulators. Real-time prediction of contact-interaction forces is challenging for systems with many degrees of freedom (DOFs) because joint and contact constraints must be enforced simultaneously. While several contact models exist for systems of free rigid bodies, fewer models are available for articulated-body systems. In this paper, we extend the method of Ruspini and Khatib and develop the contact-space resolution (CSR) model by applying the operational space theory of robot manipulation. Through a proper choice of contact-space coordinates, the projected dynamics of the system in the contact space is obtained. We show that the projection into the dynamically consistent null space preserves linear and angular momentum in a subspace of the system dynamics complementary to the joint and contact constraints. Furthermore, we illustrate that a simultaneous collision event between two articulated bodies can be resolved as an equivalent simultaneous collision between two non-articulated rigid bodies through the projected contact-space dynamics. Solving this reduced-dimensional problem is computationally efficient, but determining its accuracy requires physical experimentation. To gain further insights into the theoretical model predictions, we devised an apparatus consisting of colliding 1-, 2-, and 3-DOF articulated bodies where joint motion is recorded with high precision. Results validate that the CSR model accurately predicts the post-collision system state. Moreover, for the first time, we show that the projection of system dynamics into the mutually complementary contact space and null space is a physically verifiable phenomenon in articulated-rigid-body systems.


2020 ◽  
Vol 17 (170) ◽  
pp. 20200598 ◽  
Author(s):  
Mohammad S. Razavi ◽  
J. Brandon Dixon ◽  
Rudolph L. Gleason

The lymphatic system transports lymph from the interstitial space back to the great veins via a series of orchestrated contractions of chains of lymphangions. Biomechanical models of lymph transport, validated with ex vivo or in vivo experimental results, have proved useful in revealing novel insight into lymphatic pumping; however, a need remains to characterize the contributions of vasoregulatory compounds in these modelling tools. Nitric oxide (NO) is a key mediator of lymphatic pumping. We quantified the active contractile and passive biaxial biomechanical response of rat tail collecting lymphatics and changes in the contractile response to the exogenous NO administration and integrated these findings into a biomechanical model. The passive mechanical response was characterized with a three-fibre family model. Nonlinear regression and non-parametric bootstrapping were used to identify best-fit material parameters to passive cylindrical biaxial mechanical data, assessing uniqueness and parameter confidence intervals; this model yielded a good fit ( R 2 = 0.90). Exogenous delivery of NO via sodium nitroprusside (SNP) elicited a dose-dependent suppression of contractions; the amplitude of contractions decreased by 30% and the contraction frequency decreased by 70%. Contractile function was characterized with a modified Rachev–Hayashi model, introducing a parameter that is related to SNP concentration; the model provided a good fit ( R 2 = 0.89) to changes in contractile responses to varying concentrations of SNP. These results demonstrated the significant role of NO in lymphatic pumping and provide a predictive biomechanical model to integrate the combined effect of mechanical loading and NO on lymphatic contractility and mechanical response.


Sign in / Sign up

Export Citation Format

Share Document