M. Lengsfeld

Comparison of geometry-based and CT voxel-based finite element modelling and experimental validation.

The objectives of the paper presented here is the comparison of a geometry-based and voxel-based finite element (FE) method preprocessor of the human femur. The models were experimentally validated by strain gauge measurements (principal stress). The correlation coefficients (r) between the three methods (geometry-based FEM, voxel-based FEM, strain gauge measurements) were found to be in the range 0.91-0.94 (r2: 0.84-0.88). The relationships between the samples are highly significant (P = 0.001), where the strain gauge results are the independent variables. These results suggest that the validity with respect to the principal stress of a voxel-based modelling is similar to the validity of geometry-based modelling. In summary, therefore, we conclude that voxel-based meshing allows a straightforward interfacing with computerized tomography (CT) scans and might contribute to a clinically applicable FEM technology.

Lumbar spine curvature during office chair sitting

Prolonged sitting is generally accepted as a high risk factor in low back pain and it is frequently suggested that a lordotic posture of the lumbar spine should be maintained during sitting. We asked whether the sagittal curvature of the lumbar spine during sitting is affected by the seat tilt, backrest and the direction of the synchronised mechanism of the back and seat tilt (synchro tilt). Two office chairs were tested by multibody analysis interfacing a human model with a chair model. Results indicate that a synchronised mechanism of an office chair representing a posterior tilt of the seat while the backrest is reclined maintains an evenly distributed lumbar lordosis. The segmental angles are between 3.1 and 3.6 degrees at the lumbar vertebrae 1/2-4/5 (L1/2-L4/5). These lumbar spine segmental angles are not sensitive to the backrest height. In contrast, a synchro tilt concept with a reduction of the seats posterior tilt while the backrest is reclined causes a strong reduction of the lumbar lordosis in backrest recline with a maximum reduction from 11.7 to 2.8 degrees in L4/5. As a consequence of these results, a synchro tilt concept with a posterior tilt of the seat while the backrest is reclined is preferable from the lumbar spine kinematics point of view.

Mechanical effects of continuous passive motion on the lumbar spine in seating

The aim of this study was to develop a model which describes the mechanical spinal response to small alternating pelvic stimulation induced by an active rotational movement of a normal chair. The rotary continuous passive motion (RCPM) of the seat about a vertical axis of only 0.6 degrees resulted in an increased in spinal length as opposed to the normal daily shrinkage, and back patients experienced pain relief. Passive and active exercies have been broadly applied for treating and healing spinal disorders. A rigid body package (ADAMS Android) was used to translate the stimulation of the ischial tubersoity in caudo-cranial handing-over visualisation. The parameters of the model were set so that the values of the global stiffness and geometry of the intervertebral discs could be changed. In vivo validation of the model was based on force and moment measurements using an internal AO fixator. The predicitons of the model concerning natural frequency (4.5Hz) in vertical direction and the axial torsion response on small pelvic torsion are comparable with experimental data.

Validation data for periprosthetic bone remodelling theories

Periprosthetic adaptive bone remodelling after total hip arthroplasty (THA) has been frequently simulated in computer models, combining bone remodelling theory with finite element analysis. Unfortunately, there still subsist a lack of clinical data, which are necessary for validation of these simulation results. Therefore, the objective of the current project is to collect prospective volumetric bone density data with a clinical computerized tomography study in seven patients after THA. A retrospective study 12 years after implantation in 11 patients was added. A data set of about 100000 bone voxels for each femur was collected. In all prospective cases, the predominant change is seen during the first year. The average density reduction in the horizontal slices was between 50 and 150 Hounsfield units (HU) (approx. 10%; p<0.001) after 2 years. Loss of density is particularly strong distal of the minor trochanter and decreases from proximal to distal. For the 12 years retrospective study, the contralateral femur provided the control. Similar trends comparable to the prospective 2-year follow-up CT density values were seen in most cases with density reductions of up to 400 HU (30%). However, in one of these cases there was no difference between the operated and the control density. As far as we are aware, this is the first collection of fully prospective 3D validation data in vivo for periprosthetic adaptive bone remodelling theories. The data are also unique as they are suitable for direct patient-specific 3D finite element meshing and individual weight-related loading.

Size and direction of hip joint forces associated with various positions of the acetabulum.

When total hip replacement is performed, the position of the acetabular component may affect wear and component survival time. We considered the questions: In what way does displacement of the hip joint center alter (1) the magnitude and (2) the direction of the resultant force? Biomechanical tests were carried out on a human multibody model. After displacement of the hip joint center, the resultant forces were calculated for the single leg stance. With the flexed single leg stance, maximum hip joint forces were observed with lateral, cranial, posterior displacement. The peak forces were affected by the modeling of a gluteus maximus wrapping point at the ischial tuberosity and were overestimated when this was removed. With the straight single leg stance, posterior displacement decreases the total load on the hip joint because of the increased leverage of the rectus femoris. With regard to the direction of the resultant force, medial displacement increases the angles in both planes, cranial displacement increases it in the sagittal plane (cranial, posterior-caudal, anterior), and anterior displacement decreases the angle in the sagittal plane and increases it in the frontal plane (medial, cranial-lateral, caudal). The direction of the force is relatively insensitive to displacement of the hip joint center. The results presented here indicate a marked increase in the force after lateral, cranial, posterior displacement of the center in the flexed single leg stance. To avoid extreme joint loading and to reduce the wear after total hip arthroplasty, the cranial and posterior regions of the acetabulum should be fully reconstructed. A high hip joint center has an adverse effect on the magnitude of the force, although the directions are hardly affected by it.

Functions of hip joint muscles

The relative torque components (percentages of total hip torque) of 26 hip muscles about three joint axes were calculated using a human multibody model to determine muscle actions as a function of hip joint angle. It was demonstrated that the muscle actions obtained from this straight line model of hip musculature were in good agreement with published functions; the data on muscle function changes during hip joint motions exceed the data available in the literature and could extend the anatomical knowledge which is available today.

Development of a hybrid finite element model for individual simulation of intertrochanteric osteotomies

Intertrochanteric osteotomies of the proximal femur are used to improve the anatomy and function of the hip joint in a number of orthopaedic diseases. To investigate the geometrical and biomechanical aspects of pre-operative planning we created a set of programs to automatically perform a simulation of intertrochanteric osteotomies on a three-dimensional finite element model of the human proximal femur based on computed tomography (CT) data and using uniform brick-shaped elements. To eliminate artefacts resulting from the rough surface of the brick elements, the femoral head was represented by a tetrahedron-based head that included a cartilage layer and a subchondral cortical zone. Applicability of the procedure was tested by performing a parametric study using a model created from CT scans taken in vivo, by applying individually calculated force conditions for the one-leg stance situation. We found a large influence of osteotomy angle on the observed stress in the femoral head cartilage, especially in a situation with insufficient containment of the femoral head. The model presented here is a biomechanical tool to simulate intertrochanteric osteotomies patient-specifically for a better understanding of the effects of such operations in the individual case. The open design of the described programs allows future interfacing with surgical navigation and robot systems.

Lengths and lever arms of hip joint muscles: geometrical analyses using a human multibody model

Abstract Multibody analysis was applied in order to construct a model of the human body including the large joints and 31 muscles crossing the hip and knee. The model represents the 50th percentile rank of a male adult. The model was used to study the influence of hip joint angles on the distances between muscular insertions and origins and on muscular lever arms. Maximum and minimum muscle lengths were determined during inplane and combined hip joint motions. Maximum and minimum lever arms and their flexion/extension, abduction/adduction and rotation components during inplane hip motions were computed. A good qualitative agreement was obtained with the literature data; being suitably compared with the results presented here. It is concluded that a large geometrical data base on hip joint muscle lengths and lever arms is provided which was not previously available and which may be of significant value in morphometry, musculoskeletal research and biomechanical modelling.

The long-term (8-12 years) results of valgus and lengthening osteotomy of the femoral neck.

Abstract Shortening of the femoral neck and proximal displacement of the greater trochanter are the principal complications following avascular necrosis of the capital epiphysis head in early childhood. We report here the long-term follow-up of a series of osteotomies performed to lengthen the femoral neck and thus to restore the normal anatomy and function of the hip joint. Out of a sample of 24 patients, 15 (62.5%) were reviewed at the end of a mean follow-up of 10 years and 2 months. Fourteen of them (93.3%) had originally complained of pain on walking which varied in severity. After the operation, five were free from pain at the end of the follow-up, 9 still experienced pain on walking, and 1 also reported pain at rest. Before the operation, a positive Trendelenburg’s sign was found in 8 of the patients. This was still present in 2 at follow-up. The average perpendicular distance from the center of the femoral head to a horizontal plane passing through the tip of the greater trochanter (the centro-trochanteric distance) was reduced from 33 mm to 7 mm. Apart from correcting the anatomical deformity, it is suggested that this operation may well provide lasting relief from pain and increase the power of the abductor muscles. It may also delay the onset of osteoarthritic change and in this way postpone the necessity for an endoprosthesis.

Automatisierte Generierung von 3-D Finite Elemente Codes des menschlichen Femurs - Automatic Preprocessing of 3-D Finite Element Codes of the Human Femur

Schlüsselwörter: FEM, Modellgenerierung, humaner Femur Stark unregelmäßige geometrische Vorgaben des Knochens erschweren eine dreidimensionale Aufbereitung und regelmäßige Vernetzung im Rahmen der Methode der Finiten Elemente (FEM) eine auch in der Biomechanik zunehmend an Bedeutung gewinnende Simulationsmethode für Spannungsvorhersagen und Streßkompatibilitäts-Tests von Endoprothesen. Es wurde daher für das Beispiel Femur ein Präprozessor entwickelt, der, gestützt auf Konturdaten nachbearbeiteter Computertomogramme, eine schnelle, flexible und automatisierte Generierung von berechnungsfähigen 3-D-FEM-Modellen des Femurs gestattet. Besonderes Kennzeichen ist eine regelmäßige und standardisierte Vernetzung mit hinsichtlich Rechenzeit, -genauigkeit und Speicherplatzbedarf sehr günstigen Quaderelementen. Neben einer wesentlichen Vereinfachung der Modellgenerierung leistet das entwickelte Softwarepaket auch einen Beitrag zur besseren Untersuchbarkeit individueller Variationen und schafft die Voraussetzung für eine Etablierung der Methode in der Klinik.