251
|
Fatigue responses of the human cervical spine intervertebral discs. J Mech Behav Biomed Mater 2016; 69:30-38. [PMID: 28033533 DOI: 10.1016/j.jmbbm.2016.11.026] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 11/22/2016] [Accepted: 11/30/2016] [Indexed: 11/20/2022]
Abstract
Numerous studies have been conducted since more than fifty years to understand the behavior of the human lumbar spine under fatigue loading. Applications have been largely driven by low back pain and human body vibration problems. The human neck also sustains fatigue loading in certain type of civilian occupational and military operational activities, and research is very limited in this area. Being a visco-elastic structure, it is important to determine the stress-relaxation properties of the human cervical spine intervertebral discs to enable accurate simulations of these structures in stress-analysis models. While finite element models have the ability to incorporate viscoelastic material definitions, data specific to the cervical spine are limited. The present study was conducted to determine these properties and understand the responses of the human lower cervical spine discs under large number of cyclic loads in the axial compression mode. Eight disc segments consisting of the adjacent vertebral bodies along with the longitudinal ligaments were subjected to compression, followed by 10,000 cycles of loading at 2 or 4Hz frequency by limiting the axial load to approximately 150 N, and subsequent to resting period, subjected to compression to extract the stress-relaxation properties using the quasi-linear viscoelastic (QLV) material model. The coefficients of the model and disc displacements as a function of cycles and loading frequency are presented. The disc responses demonstrated a plateauing effect after the first 2000 to 4000 cycles, which were highly nonlinear. The paper compares these responses with the "work hardening" phenomenon proposed in clinical literature for the lumbar spine to explain the fatigue behavior of the discs. The quantitative results in terms of QLV coefficients can serve as inputs to complex finite element models of the cervical spine to delineate the local and internal load-sharing responses of the disc segment.
Collapse
|
252
|
Nogueira W, Litvak LM, Landsberger DM, Büchner A. Loudness and pitch perception using Dynamically Compensated Virtual Channels. Hear Res 2016; 344:223-234. [PMID: 27939418 DOI: 10.1016/j.heares.2016.11.017] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Revised: 11/26/2016] [Accepted: 11/29/2016] [Indexed: 10/20/2022]
Abstract
Reducing power consumption is important for the development of smaller cochlear implant (CI) speech processors. Simultaneous electrode stimulation may improve power efficiency by minimizing the required current applied to a given electrode. Simultaneous in-phase stimulation on adjacent electrodes (i.e. virtual channels) can be used to elicit pitch percepts intermediate to the ones provided by each of the physical electrodes in isolation. Virtual channels are typically implemented in monopolar stimulation mode, producing broad excitation patterns. Focused stimulation may reduce the excitation patterns, but is inefficient in terms of power consumption. To create a more power efficient virtual channel, we developed the Dynamically Compensated Virtual Channel (DC-VC) using four adjacent electrodes. The two central electrodes are current steered using the coefficient α (0<α<1 ) whereas the two flanking electrodes are used to focus/unfocus the stimulation with the coefficient σ (-1<σ<1). With increasing values of σ, power can be saved at the potential expense of generating broader electric fields. Additionally, reshaping the electric fields might also alter place pitch coding. The goal of the present study is to investigate the tradeoff between place pitch encoding and power savings using simultaneous electrode stimulation in the DC-VC configuration. A computational model and psychophysical experiments in CI users have been used for that purpose. Results from 10 adult Advanced Bionics CI users have been collected. Results show that the required current to produce comfortable levels is significantly reduced with increasing σ as predicted by the computational model. Moreover, no significant differences in the estimated number of discriminable steps were detected for the different values of σ. From these results, we conclude that DC-VCs can reduce power consumption without decreasing the number of discriminable place pitch steps.
Collapse
|
253
|
El-Zayat BF, Heyse TJ, Fanciullacci N, Labey L, Fuchs-Winkelmann S, Innocenti B. Fixation techniques and stem dimensions in hinged total knee arthroplasty: a finite element study. Arch Orthop Trauma Surg 2016; 136:1741-1752. [PMID: 27704204 DOI: 10.1007/s00402-016-2571-0] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/05/2016] [Indexed: 11/27/2022]
Abstract
INTRODUCTION No evidence-based guidelines are available to determine the appropriate stem length, and whether or not to cement stems in revision total knee arthroplasty (TKA). Therefore, the objective of this study was to compare stresses and relative movement of cemented and uncemented stems of different lengths using a finite element analysis. MATERIALS AND METHODS A finite element model was created for a synthetic tibia. Two stem lengths (95 and 160 mm) and two types of fixation (cemented or press fit) of a hinged TKA were examined. The average compressive stress distribution in different regions of interest, as well as implant micromotions, was determined and compared during lunge and squat motor tasks. RESULTS Both long and short stems in revision TKA lead to high stresses, primarily in the region around the stem tip. The presence of cement reduces the stresses in the bone in every region along the stem. Short stem configurations are less affected by the presence of cement than the long stem configuration. Press-fit stems showed higher micromotions compared to cemented stems. CONCLUSIONS Lowest stresses and micromotion were found for long cemented stems. Cementless stems showed more micromotion and increased stress levels especially at the level of the stem tip, which may explain the clinical phenomenon of stem-end pain following revision knee arthroplasty. These findings will help the surgeon with optimal individual implant choice.
Collapse
|
254
|
Bukač M, Alber M. Multi-component model of intramural hematoma. J Biomech 2016; 50:42-49. [PMID: 27876369 DOI: 10.1016/j.jbiomech.2016.11.040] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2016] [Accepted: 11/02/2016] [Indexed: 01/22/2023]
Abstract
A novel multi-component model is introduced for studying interaction between blood flow and deforming aortic wall with intramural hematoma (IMH). The aortic wall is simulated by a composite structure submodel representing material properties of the three main wall layers. The IMH is described by a poroelasticity submodel which takes into account both the pressure inside hematoma and its deformation. The submodel of the hematoma is fully coupled with the aortic submodel as well as with the submodel of the pulsatile blood flow. Model simulations are used to investigate the relation between the peak wall stress, hematoma thickness and permeability in patients of different age. The results indicate that an increase in hematoma thickness leads to larger wall stress, which is in agreement with clinical data. Further simulations demonstrate that a hematoma with smaller permeability results in larger wall stress, suggesting that blood coagulation in hematoma might increase its mechanical stability. This is in agreement with previous experimental observations of coagulation having a beneficial effect on the condition of a patient with the IMH.
Collapse
|
255
|
Diciotti S, Nobis A, Ciulli S, Landini N, Mascalchi M, Sverzellati N, Innocenti B. Development of digital phantoms based on a finite element model to simulate low-attenuation areas in CT imaging for pulmonary emphysema quantification. Int J Comput Assist Radiol Surg 2016; 12:1561-1570. [PMID: 27838881 DOI: 10.1007/s11548-016-1500-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2016] [Accepted: 10/31/2016] [Indexed: 11/28/2022]
Abstract
PURPOSE To develop an innovative finite element (FE) model of lung parenchyma which simulates pulmonary emphysema on CT imaging. The model is aimed to generate a set of digital phantoms of low-attenuation areas (LAA) images with different grades of emphysema severity. METHODS Four individual parameter configurations simulating different grades of emphysema severity were utilized to generate 40 FE models using ten randomizations for each setting. We compared two measures of emphysema severity (relative area (RA) and the exponent D of the cumulative distribution function of LAA clusters size) between the simulated LAA images and those computed directly on the models output (considered as reference). RESULTS The LAA images obtained from our model output can simulate CT-LAA images in subjects with different grades of emphysema severity. Both RA and D computed on simulated LAA images were underestimated as compared to those calculated on the models output, suggesting that measurements in CT imaging may not be accurate in the assessment of real emphysema extent. CONCLUSIONS Our model is able to mimic the cluster size distribution of LAA on CT imaging of subjects with pulmonary emphysema. The model could be useful to generate standard test images and to design physical phantoms of LAA images for the assessment of the accuracy of indexes for the radiologic quantitation of emphysema.
Collapse
|
256
|
Fradet L, Wang X, Lenke LG, Aubin CE. Biomechanical analysis of proximal junctional failure following adult spinal instrumentation using a comprehensive hybrid modeling approach. Clin Biomech (Bristol, Avon) 2016; 39:122-128. [PMID: 27750079 DOI: 10.1016/j.clinbiomech.2016.10.008] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Revised: 10/04/2016] [Accepted: 10/10/2016] [Indexed: 02/07/2023]
Abstract
BACKGROUND Proximal junctional failure is a severe proximal junctional complication following adult spinal instrumentation and involving acute proximal junctional kyphotic deformity, mechanical failure at the upper instrumented vertebra or just above, and/or proximal junctional osseoligamentous disruption. Clinical studies have identified potential risk factors, but knowledge on their biomechanics is still lacking for addressing the proximal junctional failure issues. The objective of this study was to develop comprehensive computational modeling and simulation techniques to investigate proximal junctional failure. METHODS A 3D multibody biomechanical model based on a 47year old lumbar scoliosis surgical case that subsequently had traumatic proximal junctional failure was first developed to simulate patient-specific spinal instrumentation (from T11 to S1), compute the postoperative geometry of the instrumented spine, simulate different physiological loads and movements. Then, a highly detailed finite element model of the proximal junctional spinal segment was created using as input the geometry and displacements from the multibody model. It enabled to perform detailed stress and failure analysis across the anatomical structures. FINDINGS The simulated postoperative correction and traumatic failure (wedge fracture at upper instrumented vertebra) agreed well with the clinical report (within 2° difference). Simulated stresses around the screw threads (up to 4.7MPa) generated during the instrumentation and the buckling effect of post-operative functional loads on the proximal junctional spinal segment, were identified as potential mechanical proximal junctional failure risk factors. INTERPRETATION Overall, we demonstrated the feasibility of the developed hybrid modeling technique, which realistically allowed the simulation of the spinal instrumentation and postoperative loads, which constitutes an effective tool to further investigate proximal junctional failure pathomechanisms.
Collapse
|
257
|
Areias B, Santos C, Natal Jorge RM, Gentil F, Parente MP. Finite element modelling of sound transmission from outer to inner ear. Proc Inst Mech Eng H 2016; 230:999-1007. [PMID: 27591576 DOI: 10.1177/0954411916666109] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
The ear is one of the most complex organs in the human body. Sound is a sequence of pressure waves, which propagates through a compressible media such as air. The pinna concentrates the sound waves into the external auditory meatus. In this canal, the sound is conducted to the tympanic membrane. The tympanic membrane transforms the pressure variations into mechanical displacements, which are then transmitted to the ossicles. The vibration of the stapes footplate creates pressure waves in the fluid inside the cochlea; these pressure waves stimulate the hair cells, generating electrical signals which are sent to the brain through the cochlear nerve, where they are decoded. In this work, a three-dimensional finite element model of the human ear is developed. The model incorporates the tympanic membrane, ossicular bones, part of temporal bone (external auditory meatus and tympanic cavity), middle ear ligaments and tendons, cochlear fluid, skin, ear cartilage, jaw and the air in external auditory meatus and tympanic cavity. Using the finite element method, the magnitude and the phase angle of the umbo and stapes footplate displacement are calculated. Two slightly different models are used: one model takes into consideration the presence of air in the external auditory meatus while the other does not. The middle ear sound transfer function is determined for a stimulus of 60 dB SPL, applied to the outer surface of the air in the external auditory meatus. The obtained results are compared with previously published data in the literature. This study highlights the importance of external auditory meatus in the sound transmission. The pressure gain is calculated for the external auditory meatus.
Collapse
|
258
|
Tagliabue S, Rossi E, Baino F, Vitale-Brovarone C, Gastaldi D, Vena P. Micro-CT based finite element models for elastic properties of glass-ceramic scaffolds. J Mech Behav Biomed Mater 2016; 65:248-255. [PMID: 27592293 DOI: 10.1016/j.jmbbm.2016.08.020] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Revised: 08/12/2016] [Accepted: 08/13/2016] [Indexed: 11/30/2022]
Abstract
In this study, the mechanical properties of porous glass-ceramic scaffolds are investigated by means of three-dimensional finite element models based on micro-computed tomography (micro-CT) scan data. In particular, the quantitative relationship between the morpho-architectural features of the obtained scaffolds, such as macroscopic porosity and strut thickness, and elastic properties, is sought. The macroscopic elastic properties of the scaffolds have been obtained through numerical homogenization approaches using the mechanical characteristics of the solid walls of the scaffolds (assessed through nanoindentation) as input parameters for the numerical simulations. Anisotropic mechanical properties of the produced scaffolds have also been investigated by defining a suitable anisotropy index. A comparison with morphological data obtained through the micro-CT scans is also presented. The proposed study shows that the produced glass-ceramic scaffolds exhibited a macroscopic porosity ranging between 29% and 97% which corresponds to an average stiffness ranging between 42.4GPa and 36MPa. A quantitative estimation of the isotropy of the macroscopic elastic properties has been performed showing that the samples with higher solid fractions were those closest to an isotropic material.
Collapse
|
259
|
Song JW, Lim JK, Lee KJ, Sung SJ, Chun YS, Mo SS. Finite element analysis of maxillary incisor displacement during en-masse retraction according to orthodontic mini-implant position. Korean J Orthod 2016; 46:242-52. [PMID: 27478801 PMCID: PMC4965595 DOI: 10.4041/kjod.2016.46.4.242] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2015] [Revised: 12/22/2015] [Accepted: 01/08/2016] [Indexed: 11/25/2022] Open
Abstract
Objective Orthodontic mini-implants (OMI) generate various horizontal and vertical force vectors and moments according to their insertion positions. This study aimed to help select ideal biomechanics during maxillary incisor retraction by varying the length in the anterior retraction hook (ARH) and OMI position. Methods Two extraction models were constructed to analyze the three-dimentional finite element: a first premolar extraction model (Model 1, M1) and a residual 1-mm space post-extraction model (Model 2, M2). The OMI position was set at a height of 8 mm from the arch wire between the second maxillary premolar and the first molar (low OMI traction) or at a 12-mm height in the mesial second maxillary premolar (high OMI traction). Retraction force vectors of 200 g from the ARH (-1, +1, +3, and +6 mm) at low or high OMI traction were resolved into X-, Y-, and Z-axis components. Results In M1 (low and high OMI traction) and M2 (low OMI traction), the maxillary incisor tip was extruded, but the apex was intruded, and the occlusal plane was rotated clockwise. Significant intrusion and counter-clockwise rotation in the occlusal plane were observed under high OMI traction and -1 mm ARH in M2. Conclusions This study observed orthodontic tooth movement according to the OMI position and ARH height, and M2 under high OMI traction with short ARH showed retraction with maxillary incisor intrusion.
Collapse
|
260
|
Saniei E, Setayeshi S, Akbari ME, Navid M. Parameter estimation of breast tumour using dynamic neural network from thermal pattern. J Adv Res 2016; 7:1045-1055. [PMID: 27857851 PMCID: PMC5106462 DOI: 10.1016/j.jare.2016.05.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2016] [Revised: 05/27/2016] [Accepted: 05/29/2016] [Indexed: 11/19/2022] Open
Abstract
This article presents a new approach for estimating the depth, size, and metabolic heat generation rate of a tumour. For this purpose, the surface temperature distribution of a breast thermal image and the dynamic neural network was used. The research consisted of two steps: forward and inverse. For the forward section, a finite element model was created. The Pennes bio-heat equation was solved to find surface and depth temperature distributions. Data from the analysis, then, were used to train the dynamic neural network model (DNN). Results from the DNN training/testing confirmed those of the finite element model. For the inverse section, the trained neural network was applied to estimate the depth temperature distribution (tumour position) from the surface temperature profile, extracted from the thermal image. Finally, tumour parameters were obtained from the depth temperature distribution. Experimental findings (20 patients) were promising in terms of the model's potential for retrieving tumour parameters.
Collapse
|
261
|
Yang J, Song H, Fu L, Wang M, Li W. A new solution method for wheel/rail rolling contact. SPRINGERPLUS 2016; 5:471. [PMID: 27217986 PMCID: PMC4835415 DOI: 10.1186/s40064-016-2128-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Accepted: 04/07/2016] [Indexed: 11/24/2022]
Abstract
To solve the problem of wheel/rail rolling contact of nonlinear steady-state curving, a three-dimensional transient finite element (FE) model is developed by the explicit software ANSYS/LS-DYNA. To improve the solving speed and efficiency, an explicit–explicit order solution method is put forward based on analysis of the features of implicit and explicit algorithm. The solution method was first applied to calculate the pre-loading of wheel/rail rolling contact with explicit algorithm, and then the results became the initial conditions in solving the dynamic process of wheel/rail rolling contact with explicit algorithm as well. Simultaneously, the common implicit–explicit order solution method is used to solve the FE model. Results show that the explicit–explicit order solution method has faster operation speed and higher efficiency than the implicit–explicit order solution method while the solution accuracy is almost the same. Hence, the explicit–explicit order solution method is more suitable for the wheel/rail rolling contact model with large scale and high nonlinearity.
Collapse
|
262
|
Nishida N, Kanchiku T, Imajo Y, Suzuki H, Yoshida Y, Kato Y, Nakashima D, Taguchi T. Stress analysis of the cervical spinal cord: Impact of the morphology of spinal cord segments on stress. J Spinal Cord Med 2016; 39:327-34. [PMID: 25832134 PMCID: PMC5073768 DOI: 10.1179/2045772315y.0000000012] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/31/2022] Open
Abstract
OBJECTIVE Although there are several classifications for cervical myelopathy, these do not take differences between spinal cord segments into account. Moreover, there has been no report of stress analyses for individual segments to date. METHODS By using the finite element method, we constructed 3-dimensional spinal cord models comprised of gray matter, white matter, and pia mater of the second to eighth cervical vertebrae (C2-C8). We placed compression components (disc and yellow ligament) at the front and back of these models, and applied compression to the posterior section covering 10%, 20%, 30%, or 40% of the anteroposterior diameter of each cervical spinal cord segment. RESULTS Our results revealed that, under compression applied to an area covering 10%, 20%, or 30% of the anteroposterior diameter of the cervical spinal cord segment, sites of increased stress varied depending on the morphology of each cervical spinal cord segment. Under 40% compression, stress was increased in the gray matter, lateral funiculus, and posterior funiculus of all spinal cord segments, and stress differences between the segments were smaller. CONCLUSION These results indicate that, under moderate compression, sites of increased stress vary depending on the morphology of each spinal cord segment or the shape of compression components, and also that the variability of symptoms may depend on the direction of compression. However, under severe compression, the differences among the cervical spinal segments are smaller, which may facilitate diagnosis.
Collapse
|
263
|
Wu JZ, Herzog W, Federico S. Finite element modeling of finite deformable, biphasic biological tissues with transversely isotropic statistically distributed fibers: toward a practical solution. ZEITSCHRIFT FUER ANGEWANDTE MATHEMATIK UND PHYSIK 2016; 67:26. [PMID: 27330228 PMCID: PMC4908457 DOI: 10.1007/s00033-015-0598-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The distribution of collagen fibers across articular cartilage layers is statistical in nature. Based on the concepts proposed in previous models, we developed a methodology to include the statistically distributed fibers across the cartilage thickness in the commercial FE software COMSOL which avoids extensive routine programming. The model includes many properties that are observed in real cartilage: finite hyperelastic deformation, depth-dependent collagen fiber concentration, depth- and deformation-dependent permeability, and statistically distributed collagen fiber orientation distribution across the cartilage thickness. Numerical tests were performed using confined and unconfined compressions. The model predictions on the depth-dependent strain distributions across the cartilage layer are consistent with the experimental data in the literature.
Collapse
|
264
|
3D finite element model of the chinchilla ear for characterizing middle ear functions. Biomech Model Mechanobiol 2016; 15:1263-77. [PMID: 26785845 DOI: 10.1007/s10237-016-0758-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Accepted: 01/05/2016] [Indexed: 12/16/2022]
Abstract
Chinchilla is a commonly used animal model for research of sound transmission through the ear. Experimental measurements of the middle ear transfer function in chinchillas have shown that the middle ear cavity greatly affects the tympanic membrane (TM) and stapes footplate (FP) displacements. However, there is no finite element (FE) model of the chinchilla ear available in the literature to characterize the middle ear functions with the anatomical features of the chinchilla ear. This paper reports a recently completed 3D FE model of the chinchilla ear based on X-ray micro-computed tomography images of a chinchilla bulla. The model consisted of the ear canal, TM, middle ear ossicles and suspensory ligaments, and the middle ear cavity. Two boundary conditions of the middle ear cavity wall were simulated in the model as the rigid structure and the partially flexible surface, and the acoustic-mechanical coupled analysis was conducted with these two conditions to characterize the middle ear function. The model results were compared with experimental measurements reported in the literature including the TM and FP displacements and the middle ear input admittance in chinchilla ear. An application of this model was presented to identify the acoustic role of the middle ear septa-a unique feature of chinchilla middle ear cavity. This study provides the first 3D FE model of the chinchilla ear for characterizing the middle ear functions through the acoustic-mechanical coupled FE analysis.
Collapse
|
265
|
Miller LE, Urban JE, Stitzel JD. Development and validation of an atlas-based finite element brain model. Biomech Model Mechanobiol 2016; 15:1201-14. [PMID: 26762217 DOI: 10.1007/s10237-015-0754-1] [Citation(s) in RCA: 68] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2015] [Accepted: 12/21/2015] [Indexed: 11/30/2022]
Abstract
Traumatic brain injury is a leading cause of disability and injury-related death. To enhance our ability to prevent such injuries, brain response can be studied using validated finite element (FE) models. In the current study, a high-resolution, anatomically accurate FE model was developed from the International Consortium for Brain Mapping brain atlas. Due to wide variation in published brain material parameters, optimal brain properties were identified using a technique called Latin hypercube sampling, which optimized material properties against three experimental cadaver tests to achieve ideal biomechanics. Additionally, falx pretension and thickness were varied in a lateral impact variation. The atlas-based brain model (ABM) was subjected to the boundary conditions from three high-rate experimental cadaver tests with different material parameter combinations. Local displacements, determined experimentally using neutral density targets, were compared to displacements predicted by the ABM at the same locations. Error between the observed and predicted displacements was quantified using CORrelation and Analysis (CORA), an objective signal rating method that evaluates the correlation of two curves. An average CORA score was computed for each variation and maximized to identify the optimal combination of parameters. The strongest relationships between CORA and material parameters were observed for the shear parameters. Using properties obtained through the described multiobjective optimization, the ABM was validated in three impact configurations and shows good agreement with experimental data. The final model developed in this study consists of optimized brain material properties and was validated in three cadaver impacts against local brain displacement data.
Collapse
|
266
|
Mobility and stress analysis of different surgical simulations during a sacral colpopexy, using a finite element model of the pelvic system. Int Urogynecol J 2016; 27:951-7. [PMID: 26755057 DOI: 10.1007/s00192-015-2917-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Accepted: 12/02/2015] [Indexed: 10/22/2022]
Abstract
INTRODUCTION AND HYPOTHESIS We aim to analyze the combined influence of the size of the mesh, the number of sutures, the combined use of an anterior and posterior mesh, and the tension applied to the promontory, on the mobility of the pelvic organs and on the sutures, using a Finite Element (FE) model of the female pelvic system during abdominal sacral colpopexy. METHODS We used a FE model of the female pelvic system, which allowed us to simulate the mobility of the pelvic system and to evaluate problems related to female prolapse. The meshes were added to the geometrical model and then transferred to computing software. This analysis allowed us to compare the stress and mobility during a thrust effort in different situations. RESULTS The bigger the mesh, the less mobility of both anterior and posterior organs there would be. This is accompanied by an increase in stress at the suture level. The combination of a posterior mesh with an anterior one decreases mobility and stress at the suture level. There is a particularly relevant stressing zone on the suture at the cervix. The increase in the number of sutures induces a decrease in the tension applied at each suture zone and has no impact on organ mobility. CONCLUSION Our model enables us to simulate and analyze an infinite number of surgical hypotheses. Even if these results are not validated at a clinical level, we can observe the importance of the association of both an anterior and a posterior mesh or the number of sutures.
Collapse
|
267
|
Albro MB, Nims RJ, Durney KM, Cigan AD, Shim JJ, Vunjak-Novakovic G, Hung CT, Ateshian GA. Heterogeneous engineered cartilage growth results from gradients of media-supplemented active TGF-β and is ameliorated by the alternative supplementation of latent TGF-β. Biomaterials 2015; 77:173-185. [PMID: 26599624 DOI: 10.1016/j.biomaterials.2015.10.018] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2015] [Revised: 10/05/2015] [Accepted: 10/06/2015] [Indexed: 02/07/2023]
Abstract
Transforming growth factor beta (TGF-β) has become one of the most widely utilized mediators of engineered cartilage growth. It is typically exogenously supplemented in the culture medium in its active form, with the expectation that it will readily transport into tissue constructs through passive diffusion and influence cellular biosynthesis uniformly. The results of this investigation advance three novel concepts regarding the role of TGF-β in cartilage tissue engineering that have important implications for tissue development. First, through the experimental and computational analysis of TGF-β concentration distributions, we demonstrate that, contrary to conventional expectations, media-supplemented exogenous active TGF-β exhibits a pronounced concentration gradient in tissue constructs, resulting from a combination of high-affinity binding interactions and a high cellular internalization rate. These gradients are sustained throughout the entire culture duration, leading to highly heterogeneous tissue growth; biochemical and histological measurements support that while biochemical content is enhanced up to 4-fold at the construct periphery, enhancements are entirely absent beyond 1 mm from the construct surface. Second, construct-encapsulated chondrocytes continuously secrete large amounts of endogenous TGF-β in its latent form, a portion of which undergoes cell-mediated activation and enhances biosynthesis uniformly throughout the tissue. Finally, motivated by these prior insights, we demonstrate that the alternative supplementation of additional exogenous latent TGF-β enhances biosynthesis uniformly throughout tissue constructs, leading to enhanced but homogeneous tissue growth. This novel demonstration suggests that latent TGF-β supplementation may be utilized as an important tool for the translational engineering of large cartilage constructs that will be required to repair the large osteoarthritic defects observed clinically.
Collapse
|
268
|
Bijar A, Rohan PY, Perrier P, Payan Y. Atlas-Based Automatic Generation of Subject-Specific Finite Element Tongue Meshes. Ann Biomed Eng 2015; 44:16-34. [PMID: 26577253 DOI: 10.1007/s10439-015-1497-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2015] [Accepted: 10/27/2015] [Indexed: 10/22/2022]
Abstract
Generation of subject-specific 3D finite element (FE) models requires the processing of numerous medical images in order to precisely extract geometrical information about subject-specific anatomy. This processing remains extremely challenging. To overcome this difficulty, we present an automatic atlas-based method that generates subject-specific FE meshes via a 3D registration guided by Magnetic Resonance images. The method extracts a 3D transformation by registering the atlas' volume image to the subject's one, and establishes a one-to-one correspondence between the two volumes. The 3D transformation field deforms the atlas' mesh to generate the subject-specific FE mesh. To preserve the quality of the subject-specific mesh, a diffeomorphic non-rigid registration based on B-spline free-form deformations is used, which guarantees a non-folding and one-to-one transformation. Two evaluations of the method are provided. First, a publicly available CT-database is used to assess the capability to accurately capture the complexity of each subject-specific Lung's geometry. Second, FE tongue meshes are generated for two healthy volunteers and two patients suffering from tongue cancer using MR images. It is shown that the method generates an appropriate representation of the subject-specific geometry while preserving the quality of the FE meshes for subsequent FE analysis. To demonstrate the importance of our method in a clinical context, a subject-specific mesh is used to simulate tongue's biomechanical response to the activation of an important tongue muscle, before and after cancer surgery.
Collapse
|
269
|
Pasha S, Aubin CE, Labelle H, Parent S, Mac-Thiong JM. The biomechanical effects of spinal fusion on the sacral loading in adolescent idiopathic scoliosis. Clin Biomech (Bristol, Avon) 2015; 30:981-7. [PMID: 26190098 DOI: 10.1016/j.clinbiomech.2015.06.019] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Revised: 06/23/2015] [Accepted: 06/25/2015] [Indexed: 02/07/2023]
Abstract
BACKGROUND Posterior spinal surgical correction is performed to correct spinal deformities in adolescent idiopathic scoliosis. Although the relative spino-pelvic alignment changes after spinal surgery, pelvis remains unfused in idiopathic scoliosis surgery. The impact of the spinal fusion on the transferred load to the pelvis via sacrum is not documented in the scoliotic subgroups. METHOD Bi-planar radiographs of 9 scoliotic subjects before and in average 16 months after spinal instrumentation surgery, and 12 controls were selected retrospectively. Patient-specific 3D reconstruction and finite element models of the spine, ribcage, and pelvis were developed. Spinal parameters (Cobb angles, kyphosis, lordosis), sacro-pelvic parameters (pelvic incidence, pelvic tilt, sacral slope), frontal and sagittal balances, the position of the trunk center of mass, and the centroid of the stress distribution on the sacrum superior endplate were measured and computed before operation and in the last follow-up. FINDINGS The position of the stress distribution centroid on the sacrum superior endplate with respect to the central hip vertical axis was significantly different between pre-operative and post-operative patients p<0.05. The distance between the anterior-posterior position of the trunk center of mass and the center of pressure on the superior sacral endplate significantly decreased after the spinal surgery p<0.05. INTERPRETATION The impact of the scoliosis spinal fusion on the transferred load between the spine and pelvis was evaluated. The biomechanical loading of the sacrum endplate was related to the post-operative postural balance and compensatory changes in the spino-pelvic alignment after scoliosis surgery.
Collapse
|
270
|
Cutrì E, Serrani M, Bagnoli P, Fumero R, Costantino ML. The cardiac torsion as a sensitive index of heart pathology: A model study. J Mech Behav Biomed Mater 2015; 55:104-119. [PMID: 26580023 DOI: 10.1016/j.jmbbm.2015.10.009] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2015] [Revised: 10/13/2015] [Accepted: 10/22/2015] [Indexed: 01/29/2023]
Abstract
The torsional behaviour of the heart (i.e. the mutual rotation of the cardiac base and apex) was proved to be sensitive to alterations of some cardiovascular parameters, i.e. preload, afterload and contractility. Moreover, pathologies which affect the fibers architecture and cardiac geometry were proved to alter the cardiac torsion pattern. For these reasons, cardiac torsion represents a sensitive index of ventricular performance. The aim of this work is to provide further insight into physiological and pathological alterations of the cardiac torsion by means of computational analyses, combining a structural model of the two ventricles with simple lumped parameter models of both the systemic and the pulmonary circulations. Starting from diagnostic images, a 3D anatomy based geometry of the two ventricles was reconstructed. The myocytes orientation in the ventricles was assigned according to literature data and the myocardium was modelled as an anisotropic hyperelastic material. Both the active and the passive phases of the cardiac cycle were modelled, and different clinical conditions were simulated. The results in terms of alterations of the cardiac torsion in the presence of pathologies are in agreement with experimental literature data. The use of a computational approach allowed the investigation of the stresses and strains in the ventricular wall as well as of the global hemodynamic parameters in the presence of the considered pathologies. Furthermore, the model outcomes highlight how for specific pathological conditions, an altered torsional pattern of the ventricles can be present, encouraging the use of the ventricular torsion in the clinical practice.
Collapse
|
271
|
Urinary catheter capable of repeated on-demand removal of infectious biofilms via active deformation. Biomaterials 2015; 77:77-86. [PMID: 26584348 DOI: 10.1016/j.biomaterials.2015.10.070] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2015] [Revised: 10/17/2015] [Accepted: 10/26/2015] [Indexed: 12/25/2022]
Abstract
Biofilm removal from biomaterials is of fundamental importance, and is especially relevant when considering the problematic and deleterious impact of biofilm infections on the inner surfaces of urinary catheters. Catheter-associated urinary tract infections are the most common cause of hospital-acquired infections and there are over 30 million Foley urinary catheters used annually in the USA. In this paper, we present the design and optimization of urinary catheter prototypes capable of on-demand removal of biofilms from the inner luminal surface of catheters. The urinary catheters utilize 4 intra-wall inflation lumens that are pressure-actuated to generate region-selective strains in the elastomeric urine lumen, and thereby remove overlying biofilms. A combination of finite-element modeling and prototype fabrication was used to optimize the catheter design to generate greater than 30% strain in the majority of the luminal surface when subjected to pressure. The catheter prototypes are able to remove greater than 80% of a mixed community biofilm of Proteus mirabilis and Escherichia coli on-demand, and furthermore are able to remove the biofilm repeatedly. Additionally, experiments with the prototypes demonstrate that biofilm debonding can be achieved upon application of both tensile and compressive strains in the inner surface of the catheter. The fouling-release catheter offers the potential for a non-biologic, non-antibiotic method to remove biofilms and thereby for impacting the thus far intractable problem of catheter-associated infections.
Collapse
|
272
|
Goff MG, Lambers FM, Sorna RM, Keaveny TM, Hernandez CJ. Finite element models predict the location of microdamage in cancellous bone following uniaxial loading. J Biomech 2015; 48:4142-4148. [PMID: 26522622 DOI: 10.1016/j.jbiomech.2015.10.023] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2015] [Revised: 10/10/2015] [Accepted: 10/18/2015] [Indexed: 10/22/2022]
Abstract
High-resolution finite element models derived from micro-computed tomography images are often used to study the effects of trabecular microarchitecture and loading mode on tissue stress, but the degree to which existing finite element methods correctly predict the location of tissue failure is not well characterized. In the current study, we determined the relationship between the location of highly strained tissue, as determined from high-resolution finite element models, and the location of tissue microdamage, as determined from three-dimensional fluoroscopy imaging, which was performed after the microdamage was generated in-vitro by mechanical testing. Fourteen specimens of human vertebral cancellous bone were assessed (8 male donors, 2 female donors, 47-78 years of age). Regions of stained microdamage, were 50-75% more likely to form in highly strained tissue (principal strains exceeding 0.4%) than elsewhere, and generally the locations of the regions of microdamage were significantly correlated (p<0.05) with the locations of highly strained tissue. This spatial correlation was stronger for the largest regions of microdamage (≥1,000,000μm(3) in volume); 87% of large regions of microdamage were located near highly strained tissue. Together, these findings demonstrate that there is a strong correlation between regions of microdamage and regions of high strain in human cancellous bone, particularly for the biomechanically more important large instances of microdamage.
Collapse
|
273
|
Franklyn M, Oakes B. Aetiology and mechanisms of injury in medial tibial stress syndrome: Current and future developments. World J Orthop 2015; 6:577-589. [PMID: 26396934 PMCID: PMC4573502 DOI: 10.5312/wjo.v6.i8.577] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/31/2015] [Revised: 05/31/2015] [Accepted: 07/27/2015] [Indexed: 02/06/2023] Open
Abstract
Medial tibial stress syndrome (MTSS) is a debilitating overuse injury of the tibia sustained by individuals who perform recurrent impact exercise such as athletes and military recruits. Characterised by diffuse tibial anteromedial or posteromedial surface subcutaneous periostitis, in most cases it is also an injury involving underlying cortical bone microtrauma, although it is not clear if the soft tissue or cortical bone reaction occurs first. Nuclear bone scans and magnetic resonance imaging (MRI) can both be used for the diagnosis of MTSS, but the patient’s history and clinical symptoms need to be considered in conjunction with the imaging findings for a correct interpretation of the results, as both imaging modalities have demonstrated positive findings in the absence of injury. However, MRI is rapidly becoming the preferred imaging modality for the diagnosis of bone stress injuries. It can also be used for the early diagnosis of MTSS, as the developing periosteal oedema can be identified. Retrospective studies have demonstrated that MTSS patients have lower bone mineral density (BMD) at the injury site than exercising controls, and preliminary data indicates the BMD is lower in MTSS subjects than tibial stress fracture (TSF) subjects. The values of a number of tibial geometric parameters such as cross-sectional area and section modulus are also lower in MTSS subjects than exercising controls, but not as low as the values in TSF subjects. Thus, the balance between BMD and cortical bone geometry may predict an individual's likelihood of developing MTSS. However, prospective longitudinal studies are needed to determine how these factors alter during the development of the injury and to find the detailed structural cause, which is still unknown. Finite element analysis has recently been used to examine the mechanisms involved in tibial stress injuries and offer a promising future tool to understand the mechanisms involved in MTSS. Contemporary accurate diagnosis of either MTSS or a TSF includes a thorough clinical examination to identify signs of bone stress injury and to exclude other pathologies. This should be followed by an MRI study of the whole tibia. The cause of the injury should be established and addressed in order to facilitate healing and prevent future re-occurrence.
Collapse
|
274
|
Begum MS, Dinesh MR, Tan KFH, Jairaj V, Md Khalid K, Singh VP. Construction of a three-dimensional finite element model of maxillary first molar and it's supporting structures. JOURNAL OF PHARMACY AND BIOALLIED SCIENCES 2015; 7:S443-50. [PMID: 26538895 PMCID: PMC4606637 DOI: 10.4103/0975-7406.163496] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2015] [Revised: 04/28/2015] [Accepted: 05/22/2015] [Indexed: 12/05/2022] Open
Abstract
The finite element method (FEM) is a powerful computational tool for solving stress-strain problems; its ability to handle material inhomogeneity and complex shapes makes the FEM, the most suitable method for the analysis of internal stress levels in the tooth, periodontium, and alveolar bone. This article intends to explain the steps involved in the generation of a three-dimensional finite element model of tooth, periodontal ligament (PDL) and alveolar bone, as the procedure of modeling is most important because the result is based on the nature of the modeling systems. Finite element analysis offers a means of determining strain-stress levels in the tooth, ligament, and bone structures for a broad range of orthodontic loading scenarios without producing tissue damage.
Collapse
|
275
|
Miles B, Kolos E, Walter WL, Appleyard R, Li Q, Chen Y, Ruys AJ. Subject-specific finite element model with an optical tracking system in total hip replacement surgery. Proc Inst Mech Eng H 2015; 229:280-90. [PMID: 25934257 DOI: 10.1177/0954411915578688] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Intra-operative peri-prosthetic femoral fractures are a significant concern in total hip arthroplasty and can occur at any time during surgery, with the highest incidence during implant insertion. This study combines subject-specific finite element analysis modeling with an optical tracking system to characterize the resultant strain in the bone and results of impaction during total hip replacement surgery. The use of ABG II femoral stem (Stryker Orthopaedics, Mahwah, NJ, USA) in the model yielded the following results. Hammer velocity was measured experimentally using a three-dimensional optical tracking system and these data were input into the finite element analysis model so that intra-operative loading scenario could be simulated. A quasi-static explicit simulation and a dynamic loading step using two implant-bone interface friction (0.1 and 0.4 friction coefficients) states were simulated. The maximum swing velocity of a mallet was experimentally measured at 1.5 m/s and occurred just before impaction of the hammer with implant introducer. Two friction states resulted in different results with the lower friction coefficient generating higher strains in the anterior regions of the model and higher displacement of the implant with respect to the femur when compared to the high friction state.
Collapse
|