Numerical Analysis of An Anastomotic Device

The explicit dynamic finite element method was utilized to investigate the deformation behaviour of a woven wire mesh tubular device that is used in a side-to-side anastomotic procedure for achieving gastrointestinal anastomosis. The numerical model was initially verified by comparison to experimental results that were obtained using a specialized testing mechanism. Once validated, the finite element model (FEM) was parameterized to ascertain the influence of several device parameters on its deformation behaviour. The importance of these parameters, as related to its optimal design for use in minimally invasive surgery (MIS), was subsequently ascertained and discussed.     

An approach to the mechanical constitutive modelling of arterial tissue based on homogenization and optimization

This paper is concerned with characterizing the quasistatic mechanical behaviour of arterial tissue undergoing finite deformation through hyperelastic constitutive functions. Commonly the parameters of constitutive functions are established by a process of optimization based on experimental data. Instead we construct a finite element model of a representative volume element of the material and compute its homogenized response to a range of deformations. These data are then used to provide objective functions for optimizing the parameters of two analytical models from the literature.     

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.     

Vacuum-assisted closure: microdeformations of wounds and cell proliferation

The mechanism of action of the Vacuum Assisted Closure Therapy (VAC; KCI, San Antonio, Texas), a recent novel innovation in the care of wounds, remains unknown. In vitro studies have revealed that cells allowed to stretch tend to divide and proliferate in the presence of soluble mitogens, whereas retracted cells remain quiescent. The authors hypothesize that application of micromechanical forces to wounds in vivo can promote wound healing through this cell shape-dependent, mechanical control mechanism. The authors created a computer model (finite element) of a wound and simulated VAC application. Finite element modeling is commonly used to engineer complex systems by breaking them down into simple discrete elements. In this model, the authors altered the pressure, pore diameter, and pore volume fraction to study the effects of vacuum-induced material deformations. The authors compared the morphology of deformation of this wound model with histologic sections of wounds treated with the VAC. The finite element model showed that most elements stretched by VAC application experienced deformations of 5 to 20 percent strain, which are similar to in vitro strain levels shown to promote cellular proliferation. Importantly, the deformation predicted by the model also was similar in morphology to the surface undulations observed in histologic cross-sections of the wounds. The authors hypothesize that this tissue deformation stretches individual cells, thereby promoting proliferation in the wound microenvironment. The application of micromechanical forces may be a useful method with which to stimulate wound healing through promotion of cell division, angiogenesis, and local elaboration of growth factors. Finite element modeling of the VAC device is consistent with this mechanism of action.     

A numerical study of the flow-induced vibration characteristics of a voice-producing element for laryngectomized patients

A computational model for exploring the design of a voice-producing voice prosthesis, or voice-producing element (VPE), is presented. The VPE is intended for use by laryngectomized patients who cannot benefit from current speech rehabilitation techniques. Previous experiments have focused on the design of a double-membrane voice generator as a VPE. For optimization studies, a numerical model has been developed. The numerical model introduced incorporates the finite element (FE) method to solve for the flow-induced vibrations of the VPE system, including airflow coupled with a mass-loaded membrane. The FE model includes distinct but coupled fluid and solid domains. The flow solver is governed by the incompressible, laminar, unsteady Navier–Stokes equations. The solid solver allows for large deformation, large strain, and collision. It is first shown that the model satisfactorily represents previously published experimental results in terms of frequency and flow rate, enabling the model for use as a design tool. The model is then used to study the influence of geometric scaling, membrane thickness, membrane stiffness, and slightly convergent or divergent channel geometry on the model response. It is shown that physiological allowable changes in the latter three device parameters alone will not be sufficient to generate the desired reduction in fundamental frequency. However, their effects are quantified and it is shown that membrane stiffness and included angle should be considered in future designs.     

Ventricular mechanics in diastole: material parameter sensitivity

Models of ventricular mechanics have been developed over the last 20 years to include finite deformation theory, anisotropic and inhomogeneous material properties and an accurate representation of ventricular geometry. As computer performance and the computational efficiency of the models improve, clinical application of these heart mechanics models is becoming feasible. One such application is to estimate myocardial material properties by adjusting the constitutive parameters to match wall deformation from MRI or ultrasound measurements, together with a measurement (or estimate) of ventricular pressure. Pigs are now the principal large animal model for these studies and in this paper we present the development of a new three-dimensional finite element model of the heart based on measurements of the geometry and the fibre and sheet orientations of pig hearts. The end-diastolic deformation of the model is computed using the "pole-zero" constitutive law which we have previously used to model the mechanics of passive myocardial tissue specimens. The sensitivities of end-diastolic fibre-sheet material strains and heart shape to changes in the material parameters are computed for the parameters of the pole-zero law in order to assess the utility of the models for inverse material property determination.     

Strain measurement in coronary arteries using intravascular ultrasound and deformable images

Atherosclerotic plaque rupture is responsible for the majority of myocardial infarctions and acute coronary syndromes. Rupture is initiated by mechanical failure of the plaque cap, and thus study of the deformation of the plaque in the artery can elucidate the events that lead to myocardial infarction. Intravascular ultrasound (IVUS) provides high resolution in vitro and in vivo cross-sectional images of blood vessels. To extract the deformation field from sequences of IVUS images, a registration process must be performed to correlate material points between image pairs. The objective of this study was to determine the efficacy of an image registration technique termed Warping to determine strains in plaques and coronary arteries from paired IVUS images representing two different states of deformation. The Warping technique uses pointwise differences in pixel intensities between image pairs to generate a distributed body force that acts to deform a finite element model. The strain distribution estimated by image-based Warping showed excellent agreement with a known forward finite element solution, representing the gold standard, from which the displaced image was created. The Warping technique had a low sensitivity to changes in material parameters or material model and had a low dependency on the noise present in the images. The Warping analysis was also able to produce accurate strain distributions when the constitutive model used for the Warping analysis and the forward analysis was different. The results of this study demonstrate that Warping in conjunction with in vivo IVUS imaging will determine the change in the strain distribution resulting from physiological loading and may be useful as a diagnostic tool for predicting the likelihood of plaque rupture through the determination of the relative stiffness of the plaque constituents.     

The influence and biomechanical role of cartilage split line pattern on tibiofemoral cartilage stress distribution during the stance phase of gait

This study presents an evaluation of the role that cartilage fibre ‘split line’ orientation plays in informing femoral cartilage stress patterns. A two-stage model is presented consisting of a whole knee joint coupled to a tissue-level cartilage model for computational efficiency. The whole joint model may be easily customised to any MRI or CT geometry using free-form deformation. Three ‘split line’ patterns (medial–lateral, anterior–posterior and random) were implemented in a finite element model with constitutive properties referring to this ‘split line’ orientation as a finite element fibre field. The medial–lateral orientation was similar to anatomy and was derived from imaging studies. Model predictions showed that ‘split lines’ are formed along the line of maximum principal strains and may have a biomechanical role of protecting the cartilage by limiting the cartilage deformation to the area of higher cartilage thickness.     

A finite element formulation and program to study transient swelling and load-carriage in healthy and degenerate articular cartilage

The theory of poroelasticity is extended to include physico-chemical swelling and used to predict the transient responses of normal and degenerate articular cartilage to both chemical and mechanical loading; with emphasis on isolating the influence of the major parameters which govern its deformation. Using a new hybrid element, our mathematical relationships were implemented in a purpose-built poroelastic finite element analysis algorithm (u-pi-c program) which was used to resolve the nature of the coupling between the mechanical and chemical responses of cartilage when subjected to ionic transport across its membranous skeleton. Our results demonstrate that one of the roles of the strain-dependent matrix permeability is to limit the rate of transmission of stresses from the fluid to the collagen-proteoglycan solid skeleton in the incipient stages of loading, and that the major contribution of the swelling pressure is that of preventing any excessive deformation of the matrix.     

基于心脏模型推断心肌病理状态的仿真研究

采用有限元方法建立了人体心脏力学模型, 并基于该仿真模型提出了一种无创推断心肌病理状态的新方案,并通过心肌梗塞定位试验加以验证。仿真研究表明这种基于模型参数优化解的方法能将心肌梗塞定位到模型基本单元的空间范围,可用于提取心肌疾病信息。     

A transversely isotropic viscoelastic constitutive equation for brainstem undergoing finite deformation

The objective of this study was to define the constitutive response of brainstem undergoing finite shear deformation. Brainstem was characterized as a transversely isotropic viscoelastic material and the material model was formulated for numerical implementation. Model parameters were fit to shear data obtained in porcine brainstem specimens undergoing finite shear deformation in three directions: parallel, perpendicular, and cross sectional to axonal fiber orientation and determined using a combined approach of finite element analysis (FEA) and a genetic algorithm (GA) optimizing method. The average initial shear modulus of brainstem matrix of 4-week old pigs was 12.7 Pa, and therefore the brainstem offers little resistance to large shear deformations in the parallel or perpendicular directions, due to the dominant contribution of the matrix in these directions. The fiber reinforcement stiffness was 121.2 Pa, indicating that brainstem is anisotropic and that axonal fibers have an important role in the cross-sectional direction. The first two leading relative shear relaxation moduli were 0.8973 and 0.0741, respectively, with corresponding characteristic times of 0.0047 and 1.4538 s, respectively, implying rapid relaxation of shear stresses. The developed material model and parameter estimation technique are likely to find broad applications in neural and orthopaedic tissues.     

Biomechanical properties of human articular cartilage under compressive loads

The function of articular cartilage is to support and distribute loads and to provide lubrication in the diarthrodial joints. Cartilage function is described by proper mechanical and rheological properties, strain and depth-dependent, which are not completely assessed. Unconfined and confined compression are commonly used to evaluate the Young's modulus (E) and the aggregate modulus (H(A)), respectively. The Poisson's ratio (nu) can be calculated indirectly from the equilibrium compression data, or using the biphasic indentation technique; it has recently been optically evaluated by using video microscopy during unconfined compression. The transient response of articular cartilage during confined compression depends on its permeability k; a constant value of k can be easily identified by a simple analytical model of confined compression tests, whereas more complex models or direct measurements (permeation tests) are needed to study the permeability dependence on deformation. A poroelastic finite element model of articular cartilage was developed for this purpose. The elastic parameters (E,nu) of the model were evaluated performing unconfined compression creep tests on human articular cartilage disks, whereas k was identified from the confined test response. Our combined experimental and computational method can be used to identify the parameters that define the permeability dependence on deformation, as a function of depth from articular surface.     

Evaluation of accuracy of non-linear finite element computations for surgical simulation: study using brain phantom

In this paper, the accuracy of non-linear finite element computations in application to surgical simulation was evaluated by comparing the experiment and modelling of indentation of the human brain phantom. The evaluation was realised by comparing forces acting on the indenter and the deformation of the brain phantom. The deformation of the brain phantom was measured by tracking 3D motions of X-ray opaque markers, placed within the brain phantom using a custom-built bi-plane X-ray image intensifier system. The model was implemented using the ABAQUS(TM) finite element solver. Realistic geometry obtained from magnetic resonance images and specific constitutive properties determined through compression tests were used in the model. The model accurately predicted the indentation force-displacement relations and marker displacements. Good agreement between modelling and experimental results verifies the reliability of the finite element modelling techniques used in this study and confirms the predictive power of these techniques in surgical simulation.     

A nonlinear biphasic fiber-reinforced porohyperviscoelastic model of articular cartilage incorporating fiber reorientation and dispersion

A nonlinear biphasic fiber-reinforced porohyperviscoelastic (BFPHVE) model of articular cartilage incorporating fiber reorientation effects during applied load was used to predict the response of ovine articular cartilage at relatively high strains (20%). The constitutive material parameters were determined using a coupled finite element-optimization algorithm that utilized stress relaxation indentation tests at relatively high strains. The proposed model incorporates the strain-hardening, tension-compression, permeability, and finite deformation nonlinearities that inherently exist in cartilage, and accounts for effects associated with fiber dispersion and reorientation and intrinsic viscoelasticity at relatively high strains. A new optimization cost function was used to overcome problems associated with large peak-to-peak differences between the predicted finite element and experimental loads that were due to the large strain levels utilized in the experiments. The optimized material parameters were found to be insensitive to the initial guesses. Using experimental data from the literature, the model was also able to predict both the lateral displacement and reaction force in unconfined compression, and the reaction force in an indentation test with a single set of material parameters. Finally, it was demonstrated that neglecting the effects of fiber reorientation and dispersion resulted in poorer agreement with experiments than when they were considered. There was an indication that the proposed BFPHVE model, which includes the intrinsic viscoelasticity of the nonfibrillar matrix (proteoglycan), might be used to model the behavior of cartilage up to relatively high strains (20%). The maximum percentage error between the indentation force predicted by the FE model using the optimized material parameters and that measured experimentally was 3%.     

Constitutive model of brain tissue suitable for finite element analysis of surgical procedures.

Realistic finite element modelling and simulation of neurosurgical procedures present a formidable challenge. Appropriate, finite deformation, constitutive model of brain tissue is a prerequisite for such development. In this paper, a large deformation, linear, viscoelastic model, suitable for direct use with commercially available finite element software packages such as ABAQUS is constructed. The proposed constitutive equation is of polynomial form with time-dependent coefficients. The model requires four material constants to be identified. The material constants were evaluated based on unconfined compression experiment results. The analytical as well as numerical solutions to the unconfined compression problem are presented. The agreement between the proposed theoretical model and the experiment is good for compression levels reaching 30% and for loading velocities varying over five orders of magnitude. The numerical solution using the finite element method matched the analytical solution very closely.     

Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability

Knowledge of the wall stresses in an abdominal aortic aneurysm (AAA) may be helpful in evaluating the need for surgical intervention to avoid rupture. This must be preceded by the development of a more suitable finite strain constitutive model for AAA, as none currently exists. Additionally, reliable stress analysis of in vivo AAA for the purposes of clinical diagnostics requires patient-specific values of the material parameters, which are difficult to determine noninvasively. The purpose of this work, therefore, was three-fold: (1) to develop a finite strain constitutive model for AAA; (2) to estimate the variation of model parameters within a sample population; and (3) to evaluate the sensitivity of computed stress distribution in AAA due to this biologic variation. We propose here a two parameter, hyperelastic, isotropic, incompressible material model and utilize experimental data from 69 freshly excised AAA specimens to both develop the functional form of the model and estimate its material parameters. Parametric analyses were performed via repeated finite element computations to determine the effect of varying each of the two model parameters on the stress distribution in a three-dimensional AAA model. The agreement between experimental data and the proposed functional form of the constitutive law was very good (R2 > 0.9). Our finite element simulations showed that the computed AAA wall stresses changed by only 4% or less when both the parameters were varied within the 95% confidence intervals for the patient population studied. This observation indicates that in lieu of the patient-specific material parameters, which are difficult to determine the use of population mean values is sufficiently accurate for the model to be reasonably employed in a clinical setting. We believe that this is an important advancement toward the development of a computational tool for the estimation of rupture potential for individual AAA, for which there is great clinical need.     

Static and Dynamic Mechanics Analysis on Artificial Hip Joints with Different Interface Designs by the Finite Element Method

Four different structural models of artificial joints were developed and the finite element method （FEM） was employed to investigate their mechanical characteristics under static and dynamic conditions. The materials used in the FEM calculation were ultra-high molecular weight polyethylene （UHMWPE）, 316L stainless steel, CoCrMo alloy and Ti6A14V alloy. The stress distribution, strain, and elastic deformation under static and dynamic conditions were obtained. Analysis and comparison of the ~alculation results of different models were conducted. It is shown that with the same parameters the model of a metallic femur head covered with an artificial cartilage layer is more similar to the structure of the natural human joint and its mechanical characteristics are the best of the four models.     

Statistical finite element method for real-time tissue mechanics analysis

The finite element (FE) method can accurately calculate tissue deformation. However, its low speed renders it ineffective for many biomedical applications involving real-time data processing. To accelerate FE analysis, we introduce a novel tissue mechanics simulation technique. This technique is suitable for real-time estimation of tissue deformation of specific organs, which is required in computer-aided diagnostic or therapeutic procedures. In this method, principal component analysis is used to describe each organ shape and its corresponding FE field for a pool of patients by a small number of weight factors. A mapping function is developed to relate the parameters of organ shape to their FE field counterpart. We show that irrespective of the complexity of the tissue's constitutive law or its loading conditions, the proposed technique is highly accurate and fast in estimating the FE field. Average deformation errors of less than 2% demonstrate the accuracy of the proposed technique.     

Analysis of transformation plasticity in steel using a finite element method coupled with a phase field model

An implicit finite element model was developed to analyze the deformation behavior of low carbon steel during phase transformation. The finite element model was coupled hierarchically with a phase field model that could simulate the kinetics and micro-structural evolution during the austenite-to-ferrite transformation of low carbon steel. Thermo-elastic-plastic constitutive equations for each phase were adopted to confirm the transformation plasticity due to the weaker phase yielding that was proposed by Greenwood and Johnson. From the simulations under various possible plastic properties of each phase, a more quantitative understanding of the origin of transformation plasticity was attempted by a comparison with the experimental observation.     

A Finite Element Formulation and Program to Study Transient Swelling and Load-carriage in Healthy and Degenerate Articular Cartilage

The theory of poroelasticity is extended to include physico-chemical swelling and used to predict the transient responses of normal and degenerate articular cartilage to both chemical and mechanical loading; with emphasis on isolating the influence of the major parameters which govern its deformation. Using a new hybrid element, our mathematical relationships were implemented in a purpose-built poroelastic finite element analysis algorithm (u–π–c fortran program) which was used to resolve the nature of the coupling between the mechanical and chemical responses of cartilage when subjected to ionic transport across its membranous skeleton. Our results demonstrate that one of the roles of the strain-dependent matrix permeability is to limit the rate of transmission of stresses from the fluid to the collagen-proteoglycan solid skeleton in the incipient stages of loading, and that the major contribution of the swelling pressure is that of preventing any excessive deformation of the matrix.