Tissue-level remodeling simulations of cancellous bone capture effects of in vivo loading in a rabbit model

The adaptation of cancellous bone to mechanical stimuli occurs throughout normal skeletal growth and aging, as well as in response to surgery, disease and device implantation. Previously we developed an in vivo cancellous loading model in the distal lateral femur of the rabbit. In response to daily in vivo loading for four weeks, bone mass increased, trabeculae thickened and the apparent modulus of the underlying cancellous bone increased. Here, we simulated our prior in vivo rabbit loading experiment using a cell-based tissue remodeling algorithm (Mullender et al., 1994) and compared the results to the in vivo experimental data published previously. Cancellous bone tissue was added or removed from the surface of trabeculae in regions of high and low mechanical stimulus, respectively. To examine the effect of material properties on mechanically regulated adaptation, we implemented both a homogeneous material model and a model where the relative density of tissue was lower for new and surface bone tissue compared to interior tissue. The simulations captured the changes in histomorphometric parameters and mechanical properties measured in the in vivo experiment illustrating the ability of computational simulations to predict the effect of mechanically regulated adaptation on cancellous bone histomorphometry and apparent modulus.

A tissue adaptation model based on strain-dependent collagen degradation and contact-guided cell traction

Soft biological tissues adapt their collagen network to the mechanical environment. Collagen remodeling and cell traction are both involved in this process. The present study presents a collagen adaptation model which includes strain-dependent collagen degradation and contact-guided cell traction. Cell traction is determined by the prevailing collagen structure and is assumed to strive for tensional homeostasis. In addition, collagen is assumed to mechanically fail if it is over-strained. Care is taken to use principally measurable and physiologically meaningful relationships. This model is implemented in a fibril-reinforced biphasic finite element model for soft hydrated tissues. The versatility and limitations of the model are demonstrated by corroborating the predicted transient and equilibrium collagen adaptation under distinct mechanical constraints against experimental observations from the literature. These experiments include overloading of pericardium explants until failure, static uniaxial and biaxial loading of cell-seeded gels in vitro and shortening of periosteum explants. In addition, remodeling under hypothetical conditions is explored to demonstrate how collagen might adapt to small differences in constraints. Typical aspects of all essentially different experimental conditions are captured quantitatively or qualitatively. Differences between predictions and experiments as well as new insights that emerge from the present simulations are discussed. This model is anticipated to evolve into a mechanistic description of collagen adaptation, which may assist in developing load-regimes for functional tissue engineered constructs, or may be employed to improve our understanding of the mechanisms behind physiological and pathological collagen remodeling.

Functional and Biomechanical Effects of the Edge-to-Edge Repair in the Setting of Mitral Regurgitation: Consolidated Knowledge and Novel Tools to Gain Insight into Its Percutaneous Implementation

Mitral regurgitation is the most prevalent heart valve disease in the western population. When severe, it requires surgical treatment, repair being the preferred option. The edge-to-edge repair technique treats mitral regurgitation by suturing the leaflets together and creating a double-orifice valve. Due to its relative simplicity and versatility, it has become progressively more widespread. Recently, its percutaneous version has become feasible, and has raised interest thanks to the positive results of the Mitraclip(®) device. Edge-to-edge features and evolution have stimulated debate and multidisciplinary research by both clinicians and engineers. After providing an overview of representative studies in the field, here we propose a novel computational approach to the most recent percutaneous evolution of the edge-to-edge technique. Image-based structural finite element models of three mitral valves affected by posterior prolapse were derived from cine-cardiac magnetic resonance imaging. The models accounted for the patient-specific 3D geometry of the valve, including leaflet compound curvature pattern, patient-specific motion of annulus and papillary muscles, and hyperelastic and anisotropic mechanical properties of tissues. The biomechanics of the three valves throughout the entire cardiac cycle was simulated before and after Mitraclip(®) implantation, assessing the biomechanical impact of the procedure. For all three simulated MVs, Mitraclip(®) implantation significantly improved systolic leaflets coaptation, without inducing major alterations in systolic peak stresses. Diastolic orifice area was decreased, by up to 58.9%, and leaflets diastolic stresses became comparable, although lower, to systolic ones. Despite established knowledge on the edge-to-edge surgical repair, latest technological advances make its percutanoues implementation a challenging field of research. The modeling approach herein proposed may be expanded to analyze clinical scenarios that are currently critical for Mitraclip(®) implantation, helping the search for possible solutions.

Injury risk assessment of non-lethal projectile head impacts

Kinetic energy non-lethal projectiles are used to impart sufficient effect onto a person in order to deter uncivil or hazardous behavior with a low probability of permanent injury. Since their first use, real cases indicate that the injuries inflicted by such projectiles may be irreversible and sometimes lead to death, especially for the head impacts. Given the high velocities and the low masses involved in such impacts, the assessment approaches proposed in automotive crash tests and sports may not be appropriate. Therefore, there is a need of a specific approach to assess the lethality of these projectiles. In this framework, some recent research data referred in this article as “force wall approach” suggest the use of three lesional thresholds (unconsciousness, meningeal damages and bone damages) that depend on the intracranial pressure. Three corresponding critical impact forces are determined for a reference projectile. Based on the principle that equal rigid wall maximal impact forces will produce equal damage on the head, these limits can be determined for any other projectile. In order to validate the consistence of this innovative method, it is necessary to compare the results with other existing assessment methods. This paper proposes a comparison between the “force wall approach” and two different head models. The first one is a numerical model (Strasbourg University Finite Element Head Model-SUFEHM) from Strasbourg University; the second one is a mechanical surrogate (Ballistics Load Sensing Headform-BLSH) from Biokinetics.

The in vivo plantar soft tissue mechanical property under the metatarsal head: implications of tissues׳ joint-angle dependent response in foot finite element modeling

Material properties of the plantar soft tissue have not been well quantified in vivo (i.e., from life subjects) nor for areas other than the heel pad. This study explored an in vivo investigation of the plantar soft tissue material behavior under the metatarsal head (MTH). We used a novel device collecting indentation data at controlled metatarsophalangeal joint angles. Combined with inverse analysis, tissues׳ joint-angle dependent material properties were identified. The results showed that the soft tissue under MTH exhibited joint-angle dependent material responses, and the computed parameters using the Ogden material model were 51.3% and 30.9% larger in the dorsiflexed than in the neutral positions, respectively. Using derived parameters in subject-specific foot finite element models revealed only those models that used tissues׳ joint-dependent responses could reproduce the known plantar pressure pattern under the MTH. It is suggested that, to further improve specificity of the personalized foot finite element models, quantitative mechanical properties of the tissue inclusive of the effects of metatarsophalangeal joint dorsiflexion are needed.

Biomechanical analysis and clinical effects of bridge combined fixation system for femoral fractures

This work aimed to compare the stress distribution and mechanical properties of our bridge combined fixation system and commonly used metal locking plate screw system by finite element analysis and by using the Zwick/Z100 testing machine. In addition, we also investigated the clinical outcome of our bridge combined fixation system for femoral fractures in 59 patients from June 2005 to January 2013. As a result, the stress distribution in the bone plate and screws of metal locking plate screw system during walking and climbing stairs was significantly lower than that of metal locking plate screw system. No significant difference in the displacement was observed between two systems. The equivalent bending stiffness of bridge combined fixation system was significantly lower than that of metal locking plate screw system. There were no significant differences in the bending strength, yield load, and maximum force between two systems. All the cases were followed up for 12-24 months (average 18 months). The X-ray showed bone callus was formed in most patients after 3 months, and the fracture line was faint and disappeared at 6-9 months postoperatively. No serious complications, such as implant breakage and wound infection, occurred postoperatively. According to self-developed standard for bone healing, clinical outcomes were rated as excellent or good in 55 out of 59 patients (success rate: 93.2%). Therefore, our findings suggest that our bridge combined fixation system may be a promising approach for treatment of long-bone fractures.

Finite element based model predictive control for active vibration suppression of a one-link flexible manipulator

This paper presents a unique approach for active vibration control of a one-link flexible manipulator. The method combines a finite element model of the manipulator and an advanced model predictive controller to suppress vibration at its tip. This hybrid methodology improves significantly over the standard application of a predictive controller for vibration control. The finite element model used in place of standard modelling in the control algorithm provides a more accurate prediction of dynamic behavior, resulting in enhanced control. Closed loop control experiments were performed using the flexible manipulator, instrumented with strain gauges and piezoelectric actuators. In all instances, experimental and simulation results demonstrate that the finite element based predictive controller provides improved active vibration suppression in comparison with using a standard predictive control strategy.

Molecular mechanisms underlying skeletal growth arrest by cutaneous scarring

In pediatric surgeries, cutaneous scarring is frequently accompanied by an arrest in skeletal growth. The molecular mechanisms responsible for this effect are not understood. Here, we investigated the relationship between scar contracture and osteogenesis. An excisional cutaneous wound was made on the tail of neonatal mice. Finite element (FE) modeling of the wound site was used to predict the distribution and magnitude of contractile forces within soft and hard tissues. Morphogenesis of the bony vertebrae was monitored by micro-CT analyses, and vertebral growth plates were interrogated throughout the healing period using assays for cell proliferation, death, differentiation, as well as matrix deposition and remodeling. Wound contracture was grossly evident on post-injury day 7 and accompanying it was a significant shortening in the tail. FE modeling indicated high compressive strains localized to the dorsal portions of the vertebral growth plates and intervertebral disks. These predicted strain distributions corresponded to sites of increased cell death, a cessation in cell proliferation, and a loss in mineralization within the growth plates and IVD. Although cutaneous contracture resolved and skeletal growth rates returned to normal, vertebrae under the cutaneous wound remained significantly shorter than controls. Thus, localized contractile forces generated by scarring led to spatial alterations in cell proliferation, death, and differentiation that inhibited bone growth in a location-dependent manner. Resolution of cutaneous scarring was not accompanied by compensatory bone growth, which left the bony elements permanently truncated. Therefore, targeting early scar reduction is critical to preserving pediatric bone growth after surgery.

Determination of the elastic properties of rabbit vocal fold tissue using uniaxial tensile testing and a tailored finite element model

The aim of the present study was to quantify the effects of the specimen shape on the accuracy of mechanical properties determined from a shape-specific model generation strategy. Digital images of five rabbit vocal folds (VFs) in their initial undeformed conditions were used to build corresponding specific solid models. The displacement field of the VFs under uniaxial tensile test was then measured over the visible portion of the surface using digital image correlation. A three-dimensional finite element model was built, using ABAQUS, for each solid model, while imposing measured boundary conditions. An inverse-problem method was used, assuming a homogeneous isotropic linear elastic constitutive model. Unknown elastic properties were identified iteratively through an error minimization technique between simulated and measured force-time data. The longitudinal elastic moduli of the five rabbit VFs were calculated and compared to values from a simple analytical method and those obtained by approximating the cross-section as elliptical. The use of shape-specific models significantly reduced the standard deviation of the Young׳s moduli of the tested specimens. However, a non-parametric statistical analysis test, i.e., the Friedman test, yielded no statistically significant differences between the shape-specific method and the elliptic cylindrical finite element model. Considering the required procedures to reconstruct the shape-specific finite element model for each tissue specimen, it might be expedient to use the simpler method when large numbers of tissue specimens are to be compared regarding their Young׳s moduli.

Design, analysis and verification of a knee joint oncological prosthesis finite element model

BACKGROUND

The aim of this paper was to design a finite element model for a hinged PROSPON oncological knee endoprosthesis and to verify the model by comparison with ankle flexion angle using knee-bending experimental data obtained previously.

METHOD

Visible Human Project CT scans were used to create a general lower extremity bones model and to compose a 3D CAD knee joint model to which muscles and ligaments were added. Into the assembly the designed finite element PROSPON prosthesis model was integrated and an analysis focused on the PEEK-OPTIMA hinge pin bushing stress state was carried out. To confirm the stress state analysis results, contact pressure was investigated. The analysis was performed in the knee-bending position within 15.4-69.4° hip joint flexion range.

RESULTS

The results showed that the maximum stress achieved during the analysis (46.6 MPa) did not exceed the yield strength of the material (90 MPa); the condition of plastic stability was therefore met. The stress state analysis results were confirmed by the distribution of contact pressure during knee-bending.

CONCLUSION

The applicability of our designed finite element model for the real implant behaviour prediction was proven on the basis of good correlation of the analytical and experimental ankle flexion angle data.

The effects of misalignment during in vivo loading of bone: techniques to detect the proximity of objects in three-dimensional models

Theories of mechanical adaptation of bone suggest that mechanical loading causes bone formation at discrete locations within bone microstructure experiencing the greatest mechanical stress/strain. Experimental testing of such theories requires in vivo loading experiments and high-resolution finite element models to determine the distribution of mechanical stresses. Finite element models of in vivo loading experiments typically assume idealized boundary conditions with applied load perfectly oriented on the bone, however small misalignments in load orientation during an in vivo experiment are unavoidable, and potentially confound the ability of finite element models to predict locations of bone formation at the scale of micrometers. Here we demonstrate two different three-dimensional spatial correlation methods to determine the effects of misalignment in load orientation on the locations of high mechanical stress/strain in the rodent tail loading model. We find that, in cancellous bone, the locations of tissue with high stress are maintained under reasonable misalignments in load orientation (p<0.01). In cortical bone, however, angular misalignments in the dorsal direction can alter the locations of high mechanical stress, but the locations of tissue with high stress are maintained under other misalignments (p<0.01). We conclude that, when using finite element models of the rodent tail loading model, small misalignments in loading orientation do not affect the predicted locations of high mechanical stress within cancellous bone.

Varying degrees of nonlinear mechanical behavior arising from geometric differences of urogynecological meshes

Synthetic polypropylene meshes were designed to restore pelvic organ support for women suffering from pelvic organ prolapse; however, the FDA released two notifications regarding potential complications associated with mesh implantation. Our aim was to characterize the structural properties of Restorelle and UltraPro subjected to uniaxial tension along perpendicular directions, and then model the tensile behavior of these meshes utilizing a co-rotational finite element model, with an imbedded linear or fiber-recruitment local stress-strain relationship. Both meshes exhibited a highly nonlinear stress-strain behavior; Restorelle had no significant differences between the two perpendicular directions, while UltraPro had a 93% difference in the low (initial) stiffness (p=0.009) between loading directions. Our model predicted that early alignment of the mesh segments in the loading direction and subsequent stretching could explain the observed nonlinear tensile behavior. However, a nonlinear stress-strain response in the stretching regime, that may be inherent to the mesh segment, was required to better capture experimental results. Utilizing a nonlinear fiber recruitment model with two parameters A and B, we observed improved agreement between the simulations and the experimental results. An inverse analysis found A=120 MPa and B=1.75 for Restorelle (RMSE=0.36). This approach yielded A=30 MPa and B=3.5 for UltraPro along one direction (RMSE=0.652), while the perpendicular orientation resulted in A=130 MPa and B=4.75 (RMSE=4.36). From the uniaxial protocol, Restorelle was found to have little variance in structural properties along these two perpendicular directions; however, UltraPro was found to behave anisotropically.

Relationship between medial plica and medial femoral condyle--a three-dimensional dynamic finite element model

BACKGROUND

Many researches reported that the pathologic medial plica impinges on the facing medial femoral condyle during knee motion and leads to erosive changes of the articular cartilage. The purpose of this study was to construct a simplified three-dimensional dynamic finite element human knee model to evaluate the dynamics behaviour between different types of medial plicae with the facing medial femoral condyles during knee motion.

METHODS

A three-dimensional dynamic finite element model composed of femur, tibia, covering cartilage and medial plica was developed. The kinematics of this simulation model was verified by previous findings during arthroscopic examination. The validated model was used to investigate and compare the magnitudes of the cyclic pressures acting on the cartilage of the medial femoral condyles by three different types of medial plicae with various Young's moduli.

FINDINGS

All types of plicae remained in contact with the medial femoral condyles and shifted medially when the knees moved from extension to flexion. The contact pressures were positively correlated with the Young's moduli of the medial plicae. During the whole range of motion, the maximum contact pressures of all simulation scenarios occurred when the knees moved beyond 50° of flexion. When the Young's moduli of medial plicae were set greater than 60 MPa, all types of medial plicae would elicit contact pressures greater than 10 MPa on the medial femoral condyles.

INTERPRETATION

The close relationship and possible high contact pressure between fibrotic medial plica and medial femoral condyle during knee motion might be a cause of cartilage damage on the medial femoral condyle and warrants further investigation.

High strength steels, stiffness of vehicle front-end structure, and risk of injury to rear seat occupants

Previous research has shown that rear seat occupant protection has decreased over model years, and front-end stiffness is a possible factor causing this trend. In this research, the effects of a change in stiffness on protection of rear seat occupants in frontal crashes were investigated. The stiffness was adjusted by using higher strength steels (DP and TRIP), or thicker metal sheets. Finite element simulations were performed, using an LS Dyna vehicle model coupled with a MADYMO dummy. Simulation results showed that an increase in stiffness, to the extent it happened in recent model years, can increase the risk of AIS3+ head injuries from 4.8% in the original model (with a stiffness of 1,000 N/mm) to 24.2% in a modified model (with a stiffness of 2,356 N/mm). The simulations also showed an increased risk of chest injury from 9.1% in the original model to 11.8% in the modified model. Distribution of injuries from real world accident data confirms the findings of the simulations.

Analysis of the effects of surface stiffness on the contact interaction between a finger and a cylindrical handle using a three-dimensional hybrid model

Contact interactions between the hand and handle, such as the contact surface softness and contact surface curvature, will affect both physical effort and musculoskeletal fatigue, thereby the comfort and safety of power tool operations. Previous models of hand gripping can be categorized into two groups: multi-body dynamic models and finite element (FE) models. The goal of the current study is to develop a hybrid FE hand gripping model, which combines the features of conventional FE models and multi-body dynamic models. The proposed model is applied to simulate hand-gripping on a cylindrical handle with covering materials of different softness levels. The model included three finger segments (distal, middle, and proximal phalanxes), three finger joints (the distal interphalangeal (DIP), proximal interphalangeal (PIP), and metacarpophalangeal (MCP) joint), and major anatomical substructures. The model was driven by joint moments, which are the net effects of all passive and active muscular forces acting about the joints. The finger model was first calibrated by using experimental data of human subject tests, and then applied to investigate the effects of surface softness on contact interactions between a finger and a cylindrical handle. Our results show that the maximal compressive stress and strain in the soft tissues of the fingers can be effectively reduced by reducing the stiffness of the covering material.

Comparison of eight published static finite element models of the intact lumbar spine: predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centers around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

3-D finite element analysis of the influence of synovial condition in sacroiliac joint on the load transmission in human pelvic system

The anterior part of the sacroiliac joint (SIJ) is a synovial joint, with little gliding and rotary movement between the contact surfaces of SIJ during locomotion. Due to its complex structure, especially when considering the surrounding ligaments, it is difficult to construct an accurate three-dimensional (3-D) finite element model for the human pelvis. Most of the pelvic models in the previous studies were simplified with either SIJ fusing together or without the sacral bone. However, the influence of those simplifications on the load transmission in human pelvis has not been studied, so the reliability of those studies remains unclear. In this study, two 3-D pelvic models were constructed: an SIJ fusing model and an SIJ contacting model. In the SIJ fusing model, the SIJ interfaces were fused together. In the SIJ contacting model, the SIJ interfaces were just in contact with each other without fusion. Compared with the SIJ contacting model, the SIJ fusing model have smaller movements in the SIJ. The stress distribution area in the SIJ fusing model on sacroiliac cartilages was also different. Those differences contributed to the decline of tensile force in the SIJ surrounding ligaments and the re-distribution of stress in the pelvic bones. In addition, the SIJ fusing model was far less sensitive to the increase in modulus of the sacroiliac cartilages, and decrease in stiffness of the ligaments surrounding the SIJ. The presence of synovia in the SIJ had greater influence on the load transmission in the human pelvic system. Therefore, the effect of the presence of synovia should not be neglected when the biomechanical behavior of human pelvis is being studied, especially for those studies related to clinical applications.

Computational approaches for analyzing the mechanics of atherosclerotic plaques: a review

Vulnerable and stable atherosclerotic plaques are heterogeneous living materials with peculiar mechanical behaviors depending on geometry, composition, loading and boundary conditions. Computational approaches have the potential to characterize the three-dimensional stress/strain distributions in patient-specific diseased arteries of different types and sclerotic morphologies and to estimate the risk of plaque rupture which is the main trigger of acute cardiovascular events. This review article attempts to summarize a few finite element (FE) studies for different vessel types, and how these studies were performed focusing on the used stress measure, inclusion of residual stress, used imaging modality and material model. In addition to histology the most used imaging modalities are described, the most common nonlinear material models and the limited number of models for plaque rupture used for such studies are provided in more detail. A critical discussion on stress measures and threshold stress values for plaque rupture used within the FE studies emphasizes the need to develop a more location and tissue-specific threshold value, and a more appropriate failure criterion. With this addition future FE studies should also consider more advanced strain-energy functions which then fit better to location and tissue-specific experimental data.

Numerical simulation of arterial dissection during balloon angioplasty of atherosclerotic coronary arteries

Balloon angioplasty is a standard clinical treatment for symptomatic coronary artery disease. In this procedure, controlled damage is applied intraluminally to the wall of a stenotic artery. Dissection of the coronary artery is a commonly observed clinical complication of angioplasty; however, not all dissections can be detected angioscopically. This work focuses on studying the dissection mechanisms triggered during the early stages of angioplasty in an atherosclerotic coronary artery, addressing the problem by means of a parametric study based on a simplified finite element model and cohesive interface modeling. Our results emphasize the presence of several damage mechanisms, at different locations, that are triggered near the very beginning of the process and evolve competitively, depending on both geometry and material properties of the atherosclerotic vessel. Small-scale damage was evidenced, which would not be detectable by angiography or intravascular ultrasound, but could potentially be sufficient to stimulate smooth muscle cell activation, promoting late-onset complications such as restenosis.

Straightening of curved pattern of collagen fibers under load controls aortic valve shape

The network of collagen fibers in the aortic valve leaflet is believed to play an important role in the strength and durability of the valve. However, in addition to its stress-bearing role, such a fiber network has the potential to produce functionally important shape changes in the closed valve under pressure load. We measured the average pattern of the collagen network in porcine aortic valve leaflets after staining for collagen. We then used finite element simulation to explore how this collagen pattern influences the shape of the closed valve. We observed a curved or bent pattern, with collagen fibers angled downward from the commissures toward the center of the leaflet to form a pattern that is concave toward the leaflet free edge. Simulations showed that these curved fiber trajectories straighten under pressure load, leading to functionally important changes in closed valve shape. Relative to a pattern of straight collagen fibers running parallel to the leaflet free edge, the concave pattern of curved fibers produces a closed valve with a 40% increase in central leaflet coaptation height and with decreased leaflet billow, resulting in a more physiological closed valve shape. Furthermore, simulations show that these changes in loaded leaflet shape reflect changes in leaflet curvature due to modulation of in-plane membrane stress resulting from straightening of the curved fibers. This effect appears to play an important role in normal valve function and may have important implications for the design of prosthetic and tissue engineered replacement valves.

ADSEAT--Adaptive Seat to Reduce Neck Injuries for Female and Male Occupants

Neck injuries sustained in low severity vehicle crashes are of worldwide concern and the risk is higher for females than for males. The objective of the study was to provide guidance on how to evaluate protective performance of vehicle seat designs aiming to reduce the incidence of neck injuries for female and male occupants. The objective was achieved by reviewing injury risk, establishing anthropometric data of an average female, performing dynamic volunteer tests comprising females and males, and developing a finite element model, EvaRID, of an average female. With respect to injury criteria, it was concluded based on the tests that using NIC (with a lower threshold value) and Nkm (with reduced intercept values) for females would be a suitable starting point. Virtual impact simulations with seats showed that differences were found in the response of the BioRID II and EvaRID models, for certain seats.

Experimental measurement of tympanic membrane response for finite element model validation of a human middle ear

The middle ear consists of a tympanic membrane, ligaments, tendons, and three ossicles. An important function of the tympanic membrane is to deliver exterior sound stimulus to the ossicles and inner ear. In this study, the responses of the tympanic membrane in a human ear were measured and compared with those of a finite element model of the middle ear. A laser Doppler vibrometer (LDV) was used to measure the dynamic responses of the tympanic membrane, which had the measurement point on the cone of light of the tympanic membrane. The measured subjects were five Korean male adults and a cadaver. The tympanic membranes were stimulated using pure-tone sine waves at 18 center frequencies of one-third octave band over a frequency range of 200 Hz ~10 kHz with 60 and 80 dB sound pressure levels. The measured responses were converted into the umbo displacement transfer function (UDTF) with a linearity assumption. The measured UDTFs were compared with the calculated UDTFs using a finite element model for the Korean human middle ear. The finite element model of the middle ear consists of three ossicles, a tympanic membrane, ligaments, and tendons. In the finite element model, the umbo displacements were calculated under a unit sound pressure on the tympanic membrane. The UDTF of the finite element model exhibited good agreement with that of the experimental one in low frequency range, whereas in higher frequency band, the two response functions deviated from each other, which demonstrates that the finite element model should be updated with more accurate material properties and/or a frequency dependent material model.

Experimental and finite element analysis of tibial stress fractures using a rabbit model

AIM: To determine if rabbit models can be used to quantify the mechanical behaviour involved in tibial stress fracture (TSF) development.

METHODS: Fresh rabbit tibiae were loaded under compression using a specifically-designed test apparatus. Weights were incrementally added up to a load of 30 kg and the mechanical behaviour of the tibia was analysed using tests for buckling, bone strain and hysteresis. Structural mechanics equations were subsequently employed to verify that the results were within the range of values predicted by theory. A finite element (FE) model was developed using cross-sectional computer tomography (CT) images scanned from one of the rabbit bones, and a static load of 6 kg (1.5 times the rabbit's body weight) was applied to represent running. The model was validated using the experimental strain gauge data, then geometric and elemental convergence tests were performed in order to find the minimum number of cross-sectional scans and elements respectively required for convergence. The analysis was then performed using both the model and the experimental results to investigate the mechanical behaviour of the rabbit tibia under compressive load and to examine crack initiation.

RESULTS: The experimental tests showed that under a compressive load of up to 12 kg, the rabbit tibia demonstrates linear behaviour with little hysteresis. Up to 30 kg, the bone does not fail by elastic buckling; however, there are low levels of tensile stress which predominately occur at and adjacent to the anterior border of the tibial midshaft: this suggests that fatigue failure occurs in these regions, since bone under cyclic loading initially fails in tension. The FE model predictions were consistent with both mechanics theory and the strain gauge results. The model was highly sensitive to small changes in the position of the applied load due to the high slenderness ratio of the rabbit’s tibia. The modelling technique used in the current study could have applications in the development of human FE models of bone, where, unlike rabbit tibia, the model would be relatively insensitive to very small changes in load position. However, the rabbit model itself is less beneficial as a tool to understand the mechanical behaviour of TSFs in humans due to the small size of the rabbit bone and the limitations of human-scale CT scanning equipment.

CONCLUSION: The current modelling technique could be used to develop human FE models. However, the rabbit model itself has significant limitations in understanding human TSF mechanics.

The biomechanical influence of the facet joint orientation and the facet tropism in the lumbar spine

BACKGROUND CONTEXT

Facet joint orientation and facet tropism (FT) are presented as the potential anatomical predisposing factors for lumbar degenerative changes that may lead in turn to early degeneration and herniation of the corresponding disc or degenerative spondylolisthesis. However, no biomechanical study of this concept has been reported.

PURPOSE

To investigate the biomechanical influence of the facet orientation and FT on stress on the corresponding segment.

STUDY DESIGN

Finite element analysis.

METHODS

Three models, F50, F55, and F60 were simulated with different facet joint orientations (50°, 55°, and 60° relative to coronal plane) at both L2-L3 facet joints. A FT model was also simulated to represent a 50° facet joint angle at the right side and a 60° facet joint angle at the left side in the L2-L3 segment. In each model, the intradiscal pressures were investigated under four pure moments and anterior shear force. Facet contact forces at the L2-L3 segment were also analyzed under extension and torsion moments and anterior shear force. This study was supported by 5000 CHF grant of 2011 AO Spine Research Korea fund. The authors of this study have no topic-specific potential conflicts of interest related to this study.

RESULTS

The F50, F55, and F60 models did not differ in the intradiscal pressures generated under four pure moments: but under anterior shear force, the F60 and FT models showed increases of intradiscal pressure. The F50 model under extension and the F60 model under torsion each generated an increase in facet contact force. In all conditions tested, the FT model yielded the greatest increase of intradiscal pressure and facet contact force of all the models.

CONCLUSIONS

The facet orientation per se did not increase disc stress or facet joint stress prominently at the corresponding level under four pure moments, but FT could make the corresponding segment more vulnerable to external moments or anterior shear force.

Finite element modeling of embolic coil deployment: multifactor characterization of treatment effects on cerebral aneurysm hemodynamics

Endovascular coiling is the most common treatment for cerebral aneurysms. During the treatment, a sequence of embolic coils with different stiffness, shapes, sizes, and lengths is deployed to fill the aneurysmal sac. Although coil packing density has been clinically correlated with treatment success, many studies have also reported success at low packing densities, as well as recurrence at high packing densities. Such reports indicate that other factors may influence treatment success. In this study, we used a novel finite element approach and computational fluid dynamics (CFD) to investigate the effects of packing density, coil shape, aneurysmal neck size, and parent vessel flow rate on aneurysmal hemodynamics. The study examines a testbed of 80 unique CFD simulations of post-treatment flows in idealized basilar tip aneurysm models. Simulated coil deployments were validated against in vitro and in vivo deployments. Among the investigated factors, packing density had the largest effect on intra-aneurysmal velocities. However, multifactor analysis of variance showed that coil shape can also have considerable effects, depending on packing density and neck size. Further, linear regression analysis showed an inverse relationship between mean void diameter in the aneurysm and mean intra-aneurysmal velocities, which underscores the importance of coil distribution and thus coil shape. Our study suggests that while packing density plays a key role in determining post-treatment hemodynamics, other factors such as coil shape, aneurysmal geometry, and parent vessel flow may also be very important.