Stress–strain behavior of mitral valve leaflets in the beating ovine heart

Excised anterior mitral leaflets exhibit anisotropic, non-linear material behavior with pre-transitional stiffness ranging from 0.06 to 0.09 N/mm² and post-transitional stiffness from 2 to 9 N/mm². We used inverse finite element (FE) analysis to test, for the first time, whether the anterior mitral leaflet (AML), in vivo, exhibits similar non-linear behavior during isovolumic relaxation (IVR). Miniature radiopaque markers were sewn to the mitral annulus, AML, and papillary muscles in 8 sheep. Four-dimensional marker coordinates were obtained using biplane videofluoroscopic imaging during three consecutive cardiac cycles. A FE model of the AML was developed using marker coordinates at the end of isovolumic relaxation (when pressure difference across the valve is approximately zero), as the reference state. AML displacements were simulated during IVR using measured left ventricular and atrial pressures. AML elastic moduli in the radial and circumferential directions were obtained for each heartbeat by inverse FEA, minimizing the difference between simulated and measured displacements. Stress–strain curves for each beat were obtained from the FE model at incrementally increasing transmitral pressure intervals during IVR. Linear regression of 24 individual stress–strain curves (8 hearts, 3 beats each) yielded a mean (±SD) linear correlation coefficient (r²) of 0.994±0.003 for the circumferential direction and 0.995±0.003 for the radial direction. Thus, unlike isolated leaflets, the AML, in vivo, operates linearly over a physiologic range of pressures in the closed mitral valve.     

Fitted hyperelastic parameters for Human brain tissue from reported tension,compression, and shear tests

The mechanical properties of human brain tissue are the subject of interest because of their use in understanding brain trauma and in developing therapeutic treatments and procedures. To represent the behavior of the tissue, we have developed hyperelastic mechanical models whose parameters are fitted in accordance with experimental test results. However, most studies available in the literature have fitted parameters with data of a single type of loading, such as tension, compression, or shear. Recently, Jin et al. (Journal of Biomechanics 46:2795−2801, 2013) reported data from ex vivo tests of human brain tissue under tension, compression, and shear loading using four strain rates and four different brain regions. However, they do not report parameters of energy functions that can be readily used in finite element simulations. To represent the tissue behavior for the quasi-static loading conditions, we aimed to determine the best fit of the hyperelastic parameters of the hyperfoam, Ogden, and polynomial strain energy functions available in ABAQUS for the low strain rate data, while simultaneously considering all three loading modes. We used an optimization process conducted in MATLAB, calling iteratively three finite element models developed in ABAQUS that represent the three loadings. Results showed a relatively good fit to experimental data in all loading modes using two terms in the energy functions. Values for the shear modulus obtained in this analysis (897−1653 Pa) are in the range of those presented in other studies. These energy-function parameters can be used in brain tissue simulations using finite element models.     

Modelling and simulation of porcine liver tissue indentation using finite element method and uniaxial stress–strain data

We hypothesize that both compression and elongation stress–strain data should be considered for modeling and simulation of soft tissue indentation. Uniaxial stress–strain data were obtained from in vitro loading experiments of porcine liver tissue. An axisymmetric finite element model was used to simulate liver tissue indentation with tissue material represented by hyperelastic models. The material parameters were derived from uniaxial stress–strain data of compressions, elongations, and combined compression and elongation of porcine liver samples. in vitro indentation tests were used to validate the finite element simulation. Stress–strain data from the simulation with material parameters derived from the combined compression and elongation data match the experimental data best. This is due to its better ability in modeling 3D deformation since the behavior of biological soft tissue under indentation is affected by both its compressive and tensile characteristics. The combined logarithmic and polynomial model is somewhat better than the 5-constant Mooney–Rivlin model as the constitutive model for this indentation simulation.     

Mechanics of the right whale mandible: Full scale testing and finite element analysis

In an effort to better understand the mechanics of ship-whale collision and to reduce the associated mortality of the critically endangered North Atlantic right whale, a comprehensive biomechanical study has been conducted by the Woods Hole Oceanographic Institution and the University of New Hampshire. The goal of the study is to develop a numerical modeling tool to predict the forces and stresses during impact and thereby the resulting mortality risk to whales from ship strikes.Based on post-mortem examinations, jaw fracture was chosen as a fatal endpoint for the whales hit by a vessel. In this paper we investigate the overall mechanical behavior of a right whale mandible under transverse loading and develop a finite element analysis model of the bone. The equivalent elastic modulus of the cortical component of right whale mandible is found by comparing full-scale bending tests with the results of numerical modeling. The finite element model of the mandible can be used in conjunction with a vessel-whale collision event model to predict bone fracture for various ship strike scenarios.     

Development of structural and material clavicle response corridors under axial compression and three point bending loading for clavicle finite element model validation

Clavicle injuries were frequently observed in automotive side and frontal crashes. Finite element (FE) models have been developed to understand the injury mechanism, although no clavicle loading response corridors yet exist in the literature to ensure the model response biofidelity. Moreover, the typically developed structural level (e.g., force–deflection) response corridors were shown to be insufficient for verifying the injury prediction capacity of FE model, which usually is based on strain related injury criteria. Therefore, the purpose of this study is to develop both the structural (force vs deflection) and material level (strain vs force) clavicle response corridors for validating FE models for injury risk modeling. 20 Clavicles were loaded to failure under loading conditions representative of side and frontal crashes respectively, half of which in axial compression, and the other half in three point bending. Both structural and material response corridors were developed for each loading condition. FE model that can accurately predict structural response and strain level provides a more useful tool in injury risk modeling and prediction. The corridor development method in this study could also be extended to develop corridors for other components of the human body.     

Finite element analysis of dental implant loading on atrophic and non-atrophic cancellous and cortical mandibular bone – a feasibility study

The first aim of this study was to assess displacements and micro-strain induced on different grades of atrophic cortical and trabecular mandibular bone by axially loaded dental implants using finite element analysis (FEA). The second aim was to assess the micro-strain induced by different implant geometries and the levels of bone-to-implant contact (BIC) on the surrounding bone. Six mandibular bone segments demonstrating different grades of mandibular bone atrophy and various bone volume fractions (from 0.149 to 0.471) were imaged using a micro-CT device. The acquired bone STL models and implant (Brånemark, Straumann, Ankylos) were merged into a three-dimensional finite elements structure. The mean displacement value for all implants was 3.1±1.2 µm. Displacements were lower in the group with a strong BIC. The results indicated that the maximum strain values of cortical and cancellous bone increased with lower bone density. Strain distribution is the first and foremost dependent on the shape of bone and architecture of cancellous bone. The geometry of the implant, thread patterns, grade of bone atrophy and BIC all affect the displacement and micro-strain on the mandible bone. Preoperative finite element analysis could offer improved predictability in the long-term outlook of dental implant restorations.     

Biomechanics of changes in ACL and PCL material properties or prestrains in flexion under muscle force-implications in ligament reconstruction

The effects of changes in cruciate ligament material and prestrain on knee joint biomechanics following ligament reconstruction surgery by a tendon are not adequately known. A 3D nonlinear finite element model of the entire knee joint was used to investigate the joint response at different flexion angles under a quadriceps force while varying ACL and PCL initial strains or material properties. The ACL and PCL forces as well as tibiofemoral contact forces/areas substantially increased with greater ACL or PCL initial strains or stiffness. The patellofemoral contact force slightly increased whereas the tibial extensor moment slightly decreased with tenser or stiffer ACL. Reverse trends were predicted with slacker ACL. Results confirm the hypotheses that changes in the prestrain of one cruciate ligament substantially influence the force in the other cruciate ligament and the entire joint and that the use of the patellar tendon (PT) as a replacement for cruciate ligaments markedly alters the joint biomechanics with trends similar to those predicted when increasing prestrains. Forces in both ACL and PCL ligaments increased as one of them became tenser or stiffer and diminished as it became slacker. These results have important consequences in joint biomechanics following ligament injuries or replacement and tend to recommend the use of grafts with smaller prestrains (i.e. slacker than intact) when using the PT as the replacement material with stiffness greater than that of replaced ligament itself.     

人体骨骼肌肉层对心电传导和体表电位标测的影响

应用电磁场数值计算方法求解三维心电图正问题,重点考察了人体骨骼肌肉层的各向异性导电性对心电传导的影响。推导了描述人类心电场的有限元模型和边界元模型,以及二者结合的形式,并在一个包含异性导电肌肉层的三维人体模下进行了模拟计算。     

Achilles tendon stress is more sensitive to subject-specific geometry than subject-specific material properties: A finite element analysis

This study used subject-specific measures of three-dimensional (3D) free Achilles tendon geometry in conjunction with a finite element method to investigate the effect of variation in subject-specific geometry and subject-specific material properties on tendon stress during submaximal isometric loading. Achilles tendons of eight participants (Aged 25–35 years) were scanned with freehand 3D ultrasound at rest and during a 70% maximum voluntary isometric contraction. Ultrasound images were segmented, volume rendered and transformed into subject-specific 3D finite element meshes. The mean (±SD) lengths, volumes and cross-sectional areas of the tendons at rest were 62 ± 13 mm, 3617 ± 984 mm³ and 58 ± 11 mm² respectively. The measured tendon strain at 70% MVIC was 5.9 ± 1.3%. Subject-specific material properties were obtained using an optimisation approach that minimised the difference between measured and modelled longitudinal free tendon strain. Generic geometry was represented by the average mesh and generic material properties were taken from the literature. Local stresses were subsequently computed for combinations of subject-specific and generic geometry and material properties. For a given geometry, changing from generic to subject-specific material properties had little effect on the stress distribution in the tendon. In contrast, changing from generic to subject-specific geometry had a 26-fold greater effect on tendon stress distribution. Overall, these findings indicate that the stress distribution experienced by the living free Achilles tendon of a young and healthy population during voluntary loading are more sensitive to variation in tendon geometry than variation in tendon material properties.     

Modeling capsule tissue growth around disk-shaped implants: a numerical and in vivo study

We propose a new mathematical model that describes the growth of fibrous tissue around rigid, disk-shaped implants. A solution methodology based on an efficient regularized iterative method is presented to calibrate the model from some measurements of the capsule tissue concentration. Numerical results obtained with synthetic data are presented to demonstrate the ability of the proposed solution methodology to determine the model parameters corresponding to a given implant. In addition, numerical results obtained with experimental data are presented to illustrate the validity of the proposed model.     

人体骨骼肌肉层对心电传导和体表电位标测的影响

应用电磁场数值计算方法求解三维心电图正问题,重点考察了人体骨骼肌肉层的各向异性导电性对心电传导的影响。推导了描述人类心电场的有限元模型和边界元模型,以及二者结合的形式,并在一个包含各向异性导电性肌肉层的三维人体模型下进行了模拟计算。讨论了骨骼肌肉层在不同的肌肉纤维方向组合下的胸腔模型,以及纤维的各向异性导电性对体表电位图的影响     

一种处理人体骨骼肌肉层各向异性导电性的数值方法

重点讨论了应用电磁场数值计算方法求解心电图正问题和逆问题中人体骨骼肌肉层的各向异性导电性的处理方法．文中应用有限元和边界元结合的方法,构造了一个包含不同的纤维方向组合的骨骼肌肉层的三维胸腔模型,并在此模型下引入了局部坐标内的各向异性导电率向整体坐标转换的方法．据此进行的模拟计算结果以图像的方式清晰地展示了各向异性导电性对体表电位图的影响．     

有限元法计算各向异性介质球头模型电位分布

在脑电正问题研究中,脑神经元所产生的电活动可用电流偶极子来模拟.本文提出把大脑看作各向异性介质球,即同时考虑电容效应、电导效应对脑内电流偶极子产生的电位的影响,并用有限元法推导出偶极子在各向异性介质球模型中的电位分布计算公式.结果表明介质的电容效应对电位分布是有影响的,反映了大脑活组织电特性,对外来不同频率的信号刺激有不同的响应.同时有限元法对大脑内某一区域内电位分布求解表明,测量较深层组织的电特性变化敏感的特点,可获得更多的测量信息.     

Finite element analysis of customized reconstruction plates for mandibular continuity defect therapy

Large mandibular continuity defects pose a significant challenge in oral maxillofacial surgery. One solution to this problem is to use computer-guided surgical planning and additive manufacturing technology to produce patient-specific reconstruction plates. However, when designing customized plates, it is important to assess potential biomechanical responses that may vary substantially depending on the size and geometry of the defect.     

Compressive mechanical properties of atherosclerotic plaques—Indentation test to characterise the local anisotropic behaviour

Accurate material models and associated parameters of atherosclerotic plaques are crucial for reliable biomechanical plaque prediction models. These biomechanical models have the potential to increase our understanding of plaque progression and failure, possibly improving risk assessment of plaque rupture, which is the main cause of ischaemic strokes and myocardial infarction. However, experimental biomechanical data on atherosclerotic plaque tissue is scarce and shows a high variability. In addition, most of the biomechanical models assume isotropic behaviour of plaque tissue, which is a general over-simplification. This review discusses the past and the current literature that focus on mechanical properties of plaque derived from compression experiments, using unconfined compression, micro-indentation or nano-indentation. Results will be discussed and the techniques will be mutually compared. Thereafter, an in-house developed indentation method combined with an inverse finite element method is introduced, allowing analysis of the local anisotropic mechanical properties of atherosclerotic plaques. The advantages and limitations of this method will be evaluated and compared to other methods reported in literature.