Refine
Year of publication
Document Type
- conference proceeding (article) (68) (remove)
Has Fulltext
- no (68)
Is part of the Bibliography
- no (68)
Keywords
- Biomechanik (3)
- Simulation (2)
- 3D-Druck (1)
- AnyBody (1)
- Arbeitsbedingungen (1)
- Beckenbruch (1)
- Bewegungsapparat (1)
- Biomechanische Analyse (1)
- Computersimulation (1)
- Computertomographie (1)
Institute
Begutachtungsstatus
- peer-reviewed (4)
High loads can appear in the individual joints of the human foot while the driver uses the pedals, in particular, during breaking. Measuring these internal forces is very difficult or almost impossible; therefore, advanced models are necessary to perform musculoskeletal simulations. The objective of this investigation was to see what loads are acting in the individual foot joints from the phalanges to calcaneus and talus during different brake scenarios. The Glasgow-Maastricht AnyBody Foot Model with 26 separate segments, connected by joints, ligaments and muscles was used inside the AnyBody Modeling System to compute individual mid foot joint loads. The amount, the direction of the force and additionally also the load insertion point was varied for several simulations. Figure 1: Seated musculoskeletal body model with applied brake force and forces for the lateral, intermediate and medial cuneiform-navicular joint for two different brake forces. The simulation showed that for the different brake scenarios, different muscles will be activated in the human and therefore different loads are apply in the fore-and mid-foot, respectively. The torso of the subject was assumed to be fixed in the seat. Further studies are ongoing to simulate the seat as an elastic element that allows different H-point locations according to the different loadings in the foot from the brake pedal using a new inverse dynamics analysis method called force-dependent kinematics.
Soccer is one of the most popular sports all around the world. It is an injurious type of sport with a focus on lower extremities and high psychological pressure during matches. The stressor is linked with injuries and an increased musculoskeletal loading. This study investigates the influence of cognitive stress on the load profile of the knee joint. Twelve professional youth soccer players performed highly dynamic runs with and without additional cognitive stress. The runs were analysed with a musculoskeletal simulation software. The data analysis shows no difference in knee joint reaction loading under additional mental stress compared to the baseline. Yet running times are significantly lower in the baseline. While there is no increase in the joint loads, the running times indicate an altered movement behaviour when the subjects are exposed to additional mental demand.
BACKGROUND:
Validation and verification of multibody musculoskeletal models sEMG is a difficult process because of the reliability of sEMG data and the complex relationship of muscle force and sEMG.
OBJECTIVE:
This work aims at comparing experimentally recorded and simulated muscle activities considering a numerical model for crosstalk. METHODS:
For providing an experimentally derived reference data set, subjects were performing elevations of the arm, where the activities of the contemplated muscle groups were measured by sEMG sensors. Computed muscle activities were further processed and transformed into an artificial electromyographical signal, which includes a numerical crosstalk model. In order to determine whether the crosstalk model provides a better agreement with the measured muscle activities, the Pearson correlation coefficient has been computed as a qualitative way of assessing the curve progression of the data sets.
RESULTS:
The results show an improvement in the correlation coefficient between the experimental data and the simulated muscle activities when taking crosstalk into account.
CONCLUSIONS:
Although the correlation coefficient increased when the crosstalk model was utilized, it is questionable if the discretization of both, the crosstalk and the musculoskeletal model, is accurate enough.
The gold standard for the analysis of human kinematics and kinetics is a camera-based motion capture system in combination with force measurement platforms. Alternatively, inertial measurement units can be utilized to obtain human kinematics, while ground reaction forces are computed from full body dynamics. This setup represents a system independent from the spatial confinement of a gait laboratory. The aim of this study is the comparison of the two methods by the investigation of lower limb kinematics and the resulting joint reaction forces within the ankle-, knee- and hip joints. For this purpose, human motion during gait was captured simultaneously by both measurement techniques. 13 trials from 8 different test subjects were evaluated in total. IMU data was processed with a quaternion based Kalman Filter. The data sets were implemented into a musculoskeletal simulation program in order to drive a virtual human body model. Each sensor was aligned to the gravitational and magnetic field vectors of the earth. The angles of flexions, extensions and rotations were analyzed to determine kinematic differences. Joint reaction forces defined kinetic dissimilarities. The overall kinematic differences of both models yielded root mean square errors of 7.62°, 6.02°, 4.95°, 2.79°, 2.38° and 3.56° for ankle flexion, subtalar eversion, knee flexion, hip external rotation, hip abduction and hip flexion, respectively. The proximo-distal differences in force peaks between the models yielded overall for the ankle, 57.33 %Bodyweight(BW) ± 46.86 %BW (16.66 %(Maximum peak to peak) ± 13.62 %) for the knee 37.09 %BW ± 29.33 %BW (17.65 % ± 15.44 %) and 32.03 %BW ± 24.33 %BW (15.6 % ± 12.54 %) for the hip. The overall outcome of this work investigated an approach independent of the common setup of the gait laboratory, thus enabling a cheaper and more flexible technology as an alternative. However, kinematic and thus kinetic differences remain rather large. Future work aims to improve the contact criterion for the calculation of the ground reaction forces and the implementation of a full-body calibration algorithm for the IMU system in order to counteract magnetic field disturbances.
Biomechanik des Alterns
(2016)
Einführung in die Mechanik
(2013)
The development of biomimetic robots has gained research interest in the last years as it may both help under-standing processes of motion execution in biological systems as well as developping a novel generation of intelligent and energy efficient robots. However, exact model generation that builds up on observations and robot design is very time intensive. In this paper we present a novel pipeline for co-development of biomimetic hardware and simulation models based on biological Computer Tomography (CT) data. For this purpose we exploit State of the Art rapid prototyping technologies such as 3D Printing and the Neurorobotics Platform for musculoskeletal simulations in virtual environments. The co-development integrates both advantages of virtual and physical experimental models and is expected to increase development speed of controllers that can be tested on the simulated counterpart before application to a printed robot model. We demonstrate the pipeline by generating a one year old infant model as a musculoskeletal simulation model and a print-in-place 3D printed skeleton as a single movable part. Even though we hereonly introduce the initial body generation and only a first testsetup for a modular sensory and control framework, we can clearly spot advantages in terms of rapid model generation and highly biological related models. Engineering costs are reducedand models can be provided to a wide research community for controller testing in an early development phase.
Work places and conditions strains the human body, both psychologically and biomechanically. In order to analyse working conditions and in the following to improve them, detailed knowledge about the effect of the different stressors on the body is needed. This manuscript discusses methods on how to evaluate biomechanical and mental loading and its effect on the musculoskeletal system. A possible workflow for the analysis is presented.
Fragestellung
Es soll in dieser Computersimulationsstudie untersucht werden, wie der Osteosyntheseverbund Platte mit Schrauben im Verbund mit einer im mittleren Drittel gebrochenen Clavicula durch das Bewegungsausmaß
in vivo belastet ist. Was sind die grundlegenden Kräfte die auf Clavicula und Implantat wirken und welchen Einfl uss hat die Bruchform.
Methodik
Die Muskel- und Gelenkkräfte sowie die Belastung des Implantatverbundes wurden mit einer muskuloskelletalen
Simulationssoftware (AnyBody Technology, V.4) berechnet. Hierfür wurden mit einem komplexen Model des menschlichen Körpers folgende Bewegungen analysiert: eine Flexion von 160° und Abduktion 160° mit
einem Gewicht von 2 kg in der Hand. Aus CT-Patientendaten wurden zwei dreidimensionale Modelle des Clavicula-Implantat Verbundes gebildet, die sich in der Frakturform unterscheiden (Querfraktur und vertikale
Fraktur). In beiden Modellen wurde eine Claviculaosteosynthese in superiorer Position mit einer 6 Loch LCP mit 2 Schrauben pro Hauptfragment verwendet. Die Materialeigenschaften wurden aus der Dichte des Materials sowie aus Literaturdaten verwendet. Die Muskel- und Gelenkkräfte aus der muskuloskelletalen Berechnung wurden auf das Finite Elemente Modell übertragen und die Spannungen und Dehnungen des Implantat-Knochenverbundes
wurden berechnet.
Ergebnisse
Es zeigte sich, dass die simulierte in vivo Belastung stark abhängig vom Flexionswinkel ist. Das Implantat ist in der superioren Lage auf Biegung belastet, welche maximale Werte im Überschulterniveau erreicht. Die Bruchform mit anatomischer Reposition und Kontakt der Hauptfragmente zueinander führt zu einer deutlichen Entlastung des Osteosyntheseverbundes im Vergleich zu einer Bruchform mit vertikaler Fraktur.
Schlussfolgerung
Aus den Analysen ist eine Positionierung der Plattenosteosynthese für die im mittleren Drittel frakturierte Clavicula in anterior-superiorer Lage wünschenswert. Die anatomische Reposition entlastet den Osteosyntheseverbund und sollte möglichst erreicht werden. Die Nachbehandlung sollte ein Bewegungsausmaß für den Arm für 4 Wochen für einfache Bruchformen auf 70° Flexion und Abduktion limitieren und für komplexe Bruchformen diese Limitierung ausgedehnt werden.
INTRODUCTION
Several clinical studies demonstrated long-term adjacent-level effects after implantation of spinal fusion devices[1]. These effects have been reported as adjacent joint degeneration and the development of new symptoms correlating with adjacent segment degeneration[2] and the trend has therefore gone to motion preservation devices; however, these effects have not been understood very well and have not been investigated thoroughly[3].
The aim of this study is to investigate the effect of varying the stiffness of spinal fusion devices on the adjacent vertebral levels. Disc forces, moments and facet joint forces were analyzed.
METHODS
The AnyBody Modeling System was used to compute the in-vivo muscle and joint reaction forces of a musculoskeletal model. The full body model used in this study consists of 188 muscle fascicles in the lumbar spine and more than 1000 individual muscle branches in total. The model has been proposed by de Zee et al.[3], validated by Rasmussen et al.[4] and by Galibarov et al.[5]. The new model[5] determines the individual motions between vertebrae based on the equilibrium between forces acting on the vertebrae from muscles and joints and the passive stiffness in disks and ligaments, figure 1a. An adult of 1.75 m and 75 kg with a spinal implant in L4L5 was modeled. This model was subjected to a flexion-extension motion using different elastic moduli to analyze and compare to a non-implanted scenario. The analyzed variables were vertebral motion, the disc reaction forces and moments, as well as facet joint forces in the treated and the adjacent levels: L2L3, L3L4, L4L5 and L5-Sacrum.
RESULTS
When introducing a spinal fusion device in the L4L5 joint the reaction forces and moments decreased in this joint with stiffer devices leading to lower joint loads. However, in the adjacent joints, L3L4 and L5Sacrum, an increase was observed when implanting stiffer devices. Similar trends could be found for the L2L3 joint. The loads in the facet joints showed the same trends. While introducing a spinal fusion device reduced the facet joint forces in the treated joint, the loads in the adjacent facet joints were increased according to the stiffness of the implanted device, figure 1b.
DISCUSSION
While the treated disc joint showed reduced motion and loads, the adjacent levels demonstrated a significant increase. In particular, the increased facet joint forces in the adjacent levels can lead to adjacent level facet pain or accelerated facet joint degeneration. Introducing a device resulted in preventing facet contact and therefore facet joint loads, even using the device with the lowest stiffness.
CONCLUSION
The presented model shows that clinical complications such as facet joint degeneration in adjacent levels after implantation of spinal fusion device are consistent with the change in the mechanical-stimulus distribution in the system.
On modelling spine curvature dependent on muscular and external forces in multibody dynamics system
(2011)
This paper presents a computational approach for investigating effect of muscular and external forces on curvature of the lumbar spine. Multibody dynamics system is used to compute the lumbar spine curvature using a force-dependent kinematics facility, e.g. this method allows releasing some degrees of freedom in order to be computed based on the current load configuration.
In recent years, musculoskeletal computation has become a widely used tool to investigate joint and muscle forces within the human body. However, the issue of muscle fatigue is not considered adequately in most models and is a challenging task. One aspect that needs to be examined is the interaction of muscles during an exhausting task. Therefore, an experimental study was designed to analyze the changes of back muscle recruitment pattern during such exercises.
In this study 38 subjects (27 male, 11 female, height = 177±8.5 cm, weight = 74.0±13.6 kg) participated. Each subject had to perform three static and three dynamic exhausting exercises where the back muscles were loaded with subject specific forces using a dynamometer adapter especially designed for the trunk muscles. To collect the muscle activity, twelve surface electromyography sensors were applied on the back, and four on the abdominal muscles. Muscle activity and fatigue were analyzed by calculating the maximum voluntary contraction normalized signal and the median frequency. At first the fatigue of m. erector spinae and m. multifidi was analyzed, since these muscles carry the main load during the exercises. Subsequently the activity of the m. trapezius, m. rectus abdominis and m. obliquus externus were investigated to determine recruitment patterns. To gain more detailed information of these patterns a numerical model was built using the AnyBody Modeling System™. Analyzing the measurements, we can observe an increasing muscle activity during isokinetic exercises while the force is constant. Since the activity in the simulation is defined as the current force output divided by the strength of the muscle, the strength parameter was scaled down based on the measured data, assuming a linear force – activity correlation, and using a numerical algorithm considering the influence of cross talk.
The results show, that changes in recruitment pattern can be divided into three major subgroups. Prior to total exhaustion, some of the subjects show additional activation of muscles in the trapezius region, while other subjects show an additional activation of abdominal muscles, increasing the intra-abdominal pressure which supports the spine. In the third group an activation in both regions can be observed.
The numerical simulations show an increasing activity of abdominal muscles as well as muscles in the upper back. Especially the m. latissimus dorsi shows a significantly higher activity.
The results lead to the conclusion that prior to total exhaustion, additional muscles are recruited to support the main muscles. It was shown that abdominal muscles are activated to support back muscles by pressurizing the trunk cavity to delay total exhaustion as long as possible. In conclusion, the results show that changes in muscle recruitment pattern need to be considered when introducing muscle fatigue to musculoskeletal models.
Thoracolumbar spine model with articulated rigcage for the prediction of dynamic spinal loading
(2015)
BACKGROUND: Validation and verification of multibody musculoskeletal models sEMG is a difficult process because of the reliability of sEMG data and the complex relationship of muscle force and sEMG. OBJECTIVE: This work aims at comparing experimentally recorded and simulated muscle activities considering a numerical model for crosstalk. METHODS: For providing an experimentally derived reference data set, subjects were performing elevations of the arm, where the activities of the contemplated muscle groups were measured by sEMG sensors. Computed muscle activities were further processed and transformed into an artificial electromyographical signal, which includes a numerical crosstalk model. In order to determine whether the crosstalk model provides a better agreement with the measured muscle activities, the Pearson correlation coefficient has been computed as a qualitative way of assessing the curve progression of the data sets. RESULTS: The results show an improvement in the correlation coefficient between the experimental data and the simulated muscle activities when taking crosstalk into account. CONCLUSIONS: Although the correlation coefficient increased when the crosstalk model was utilized, it is questionable if the discretization of both, the crosstalk and the musculoskeletal model, is accurate enough.
Common Schanz screw systems can be used to stabilize pelvic ring fractures. In order to accommodate for different patient’s requirements, implants can be placed in cranio-caudal direction into the os ilium (T1), or into the supraacetabular bone canal, and thus, in dorso-ventral direction (T2). Whereas both techniques are currently used, no data of the biomechanical behavior is available up to this date. The aim of this study is to analyze, whether T2 shows biomechanical advantages with respect to tissue and implant stresses due to the enlarged bone-implant interface. Forces acting on the pelvis were analyzed using motion capture data of a gait cycle obtained by the utilization of a musculoskeletal simulation program. A three dimensional finite element (FE) model of the pelvis with grayscale-based material properties was generated. The muscle and joint reaction forces at toe-off were applied to the FE model and instable pelvis fractures were implemented. The osteosynthesis systems were positioned within the model in order to enable the comparison between the two different surgical techniques. Stresses and displacements were analyzed for bone tissue, fracture zone and implant.
T2 lead to approx. 30% larger displacements in the fracture zone. Von-Mises stresses were larger for T2 in the implant (80 MPa vs. 227 MPa), whereas T1 leads to larger stresses in the bone tissue (200 MPa vs. 140 MPa).
Both implantation techniques showed a good biomechanical behavior. Differences could be found with respect to tissue strains and deformations in the fracture zone. If bone quality or fracture healing are of concern, T2 or T1 should be used, respectively. However, both techniques seem to be applicable for cases with no special requirements. Further analyses aim to investigate the behavior under cyclic loading.
The AnyBody™ Modeling System (AMS) [1], is an universally used musculoskeletal simulation software using inverse dynamics. Until now, no complete human hand model is known in the AMS. Also considering other musculoskeletal software platforms, just one detailed entire hand model is recently published [2] but is only based on one subject. The aim of this work is to implement a full detailed hand model for the AMS including all extrinsic and intrinsic muscles using data by the UWB gained through an anatomical study of ten cadaver hands.
Subject-specific Musculoskeletal Simulation of Hip Dislocation Risk in Activities of Daily Living
(2011)
Ergonomic workplaces lead to fewer work-related musculoskeletal disorders and thus fewer sick days. There are various guidelines to help avoid harmful situations. However, these recommendations are often rather crude and often neglect the complex interaction of biomechanical loading and psychological stress. This study investigates whether machine learning algorithms can be used to predict mechanical and stress-related muscle activity for a standardized motion. For this purpose, experimental data were collected for trunk movement with and without additional psychological stress. Two different algorithms (XGBoost and TensorFlow) were used to model the experimental data. XGBoost in particular predicted the results very well. By combining it with musculoskeletal models, the method shown here can be used for workplace analysis but also for the development of real-time feedback systems in real workplace environments.
Measuring functional outcome after total hip replacement with subject-specific hip joint loading
(2012)
Total hip replacement is an often-performed orthopedic surgical procedure; the amount of procedures undertaken will increase since our life expectancy is growing. In order to optimize function, hip biomechanics should be restored to as near normal as possible. The goal of this pilot study was to determine whether or not it is feasible to compute the vectorial hip reaction force pathways on the head of the prosthesis and the force angles relative to the cup of the prosthesis that occur during gait in total hip replacement patients, serving as an objective measurement of the functional outcome following hip replacement. A three-dimensional gait analysis, measuring ground reaction forces and kinematics, was performed. The data retrieved from the gait analysis was used as the input for the musculoskeletal model to compute vectorial joint reaction forces for data processing. To evaluate the position and orientation of the joint reaction forces, the force path, as well as the force angles for the operated and non-operated joint, has been calculated during the stance phase of the specific leg. The force path for subject 2 on the non-operated side is only located in the posterior-lateral quarter, as is the force path for subject 1. In contrast to this subject, the force path for subject 2 at the operated hip joint can be found only within the anterior quarter of the head of the implant, where it is nearly equally distributed in the medio-lateral half of the prosthesis head. The force-inclination angles on the cup of subject 1, with respect to the plane of the socket face, indicates that the force vector is mainly positioned in the same quadrant when compared with subject 2 (in a cup-fixed coordinate system). The force-anteversion angle behaves similarly to the force-inclination angle, even when the effects are not as pronounced. The proposed methods in this article are aiming to define two functional outcomes of total hip replacement that are related to wear and rim loading. It is accepted that wear is not only a function of time, but a function of use. Owing to the methods listed in this article, we are able to determine a) the applied force and b) the sliding distance (force pathway) in a subject-specific manner. The computed hip-reaction force angles and the distance to the rim cup are a measurement for cup or rim loading, and occurs in the so-called safe-zones. This method may well give us insight into the biomechanical situation during gait, after receiving total hip replacement, that we need to fully understand the mechanisms acting on a hip joint and to prove a possible increase of functional outcome after receiving total hip replacement.
Navigated Femur First Total Hip Arthroplasty leads to improved Biomechanical Outcome after surgery
(2015)