A generalized finite difference method for modeling cardiac electrical activation on arbitrary, irregular computational meshes

A generalized finite difference (GFD) method is presented that can be used to solve the bi-domain equations modeling cardiac electrical activity. Classical finite difference methods have been applied by many researchers to the bi-domain equations. However, these methods suffer from the limitation of requiring computational meshes that are structured and orthogonal. Finite element or finite volume methods enable the bi-domain equations to be solved on unstructured meshes, although implementations of such methods do not always cater for meshes with varying element topology. The GFD method solves the bi-domain equations on arbitrary and irregular computational meshes without any need to specify element basis functions. The method is useful as it can be easily applied to activation problems using existing meshes that have originally been created for use by finite element or finite difference methods. In addition, the GFD method employs an innovative approach to enforcing nodal and non-nodal boundary conditions. The GFD method performs effectively for a range of two and three-dimensional test problems and when computing bi-domain electrical activation moving through a fully anisotropic three-dimensional model of canine ventricles.     

Effects of friction on the unconfined compressive response of articular cartilage: a finite element analysis

A finite element analysis is used to study a previously unresolved issue of the effects of platen-specimen friction on the response of the unconfined compression test; effects of platen permeability are also determined. The finite element formulation is based on the linear KLM biphasic model for articular cartilage and other hydrated soft tissues. A Galerkin weighted residual method is applied to both the solid phase and the fluid phase, and the continuity equation for the intrinsically incompressible binary mixture is introduced via a penalty method. The solid phase displacements and fluid phase velocities are interpolated for each element in terms of unknown nodal values, producing a system of first order differential equations which are solved using a standard numerical finite difference technique. An axisymmetric element of quadrilateral cross-section is developed and applied to the mechanical test problem of a cylindrical specimen of soft tissue in unconfined compression. These studies show that interfacial friction plays a major role in the unconfined compression response of articular cartilage specimens with small thickness to diameter ratios.     

Patient-specific initial wall stress in abdominal aortic aneurysms with a backward incremental method

Patient-specific wall stress simulations on abdominal aortic aneurysms may provide a better criterion for surgical intervention than the currently used maximum transverse diameter. In these simulations, it is common practice to compute the peak wall stress by applying the full systolic pressure directly on the aneurysm geometry as it appears in medical images. Since this approach does not account for the fact that the measured geometry is already experiencing a substantial load, it may lead to an incorrect systolic aneurysm shape. We have developed an approach to compute the wall stress on the true diastolic geometry at a given pressure with a backward incremental method. The method has been evaluated with a neo-Hookean material law for several simple test problems. The results show that the method can predict an unloaded configuration if the loaded geometry and the load applied are known. The effect of incorporating the initial diastolic stress has been assessed by using three patient-specific geometries acquired with cardiac triggered MR. The comparison shows that the commonly used approach leads to an unrealistically smooth systolic geometry and therefore provides an underestimation for the peak wall stress. Our backward incremental modelling approach overcomes these issues and provides a more plausible estimate for the systolic aneurysm volume and a significantly different estimate for the peak wall stress. When the approach is applied with a more complex material law which has been proposed specifically for abdominal aortic aneurysm similar effects are observed and the same conclusion can be drawn.     

Methods for quasi-linear viscoelastic modeling of soft tissue: application to incremental stress-relaxation experiments

The quasi-linear viscoelastic (QLV) model was applied to incremental stress-relaxation tests and an expression for the stress was derived for each step. This expression was used to compare two methods for normalizing stress data prior to estimating QLV parameters. The first and commonly used normalization method was shown to be strain-dependent. Thus, a second normalization method was proposed and shown to be strain-independent and more sensitive to QLV time constants. These analytical results agreed with representative tendon data. Therefore, this method for normalizing stress data was proposed for future studies of incremental stress-relaxation, or whenever comparing stress-relaxation at different strains.     

Application of Markov Chain Model to Calculate the Average Number of Times of Use of a Material in Society. An Allocation Methodology for Open-Loop Recycling. Part 2: Case Study for Steel (6 pp)

- Preamble. In this series of two papers, a methodology to calculate the average number of times a material is used in a society from cradle to grave is presented and applied to allocation of environmental impact of virgin material. Part 1 focused on methodology development and showed how the methodology works with hypothetical examples of material flows. Part 2 presents case studies for steel recycling in Japan, in which the methodology is applied and allocation of environmental impact of virgin steel is conducted. - Abstract Goal, Scope and Background. The life cycle of steel begins with the mining of iron ore from the earth. Steel is produced in steel works and used in various products. Some of the steels are recycled at the products' end of life and used as a resource for the production of new steel in electric furnaces, while the remaining steel is used just once in products before being discarded (landfilled). In this paper, case studies were conducted to analyze the average number of times the element of iron is used and its residence time in society, in which the methodology developed in Part 1 of the paper was applied. CO2 emissions caused by steel productions and recycling were allocated by the number of times the element of steel is used in a society. Results and Discussion On the basis of the material flows of steel in Japan in 2000, it was calculated that at least 70% of the BF crude iron produced in Japan in 2000 was ultimately exported. On the assumption that steel is used in other countries in the same way as it is in Japan, the average number of times of use and the residence time of elemental iron in society are 2.67 and 62.9 years, respectively. Both of these values depend significantly on the recycling ratios of steel from construction and automobiles. Our model indicated that if the recycling ratio of steel from civil engineering and construction increased from 50% to 60%, the average number of times used would increase to 3.17 and the residence time of elemental iron in society would increase to 75.8 years. If CO2 emissions caused by steel productions and recycling are allocated by the number of times the element of steel is used in a society, it was calculated that steel use of one time generates in average an environmental burden of 1.03 t-CO2/t. Conclusion A method was developed to calculate the average number of times a material is used in a society from cradle to grave. Our methodology is based on Markov chain model using matrix-based numerical analysis, and has been successfully applied to steel. The results obtained by this methodology, i.e. the average number of times the element of iron is used in society, could be used for allocation of environmental burdens of virgin material as well as an indicator for assessing the state of material use in a certain year, based on material flow of material in that year. Recommendation and Perspective It is recognized that further researches must be conducted to gather data on steel production, use, and recovery in other countries and incorporate them into the transition probability matrix to obtain more precise results. Although this paper deals only with steel, this method can also be applied to other materials.     

An efficient method of modeling material properties using a thermal diffusion analogy: an example based on craniofacial bone

The ability to incorporate detailed geometry into finite element models has allowed researchers to investigate the influence of morphology on performance aspects of skeletal components. This advance has also allowed researchers to explore the effect of different material models, ranging from simple (e.g., isotropic) to complex (e.g., orthotropic), on the response of bone. However, bone's complicated geometry makes it difficult to incorporate complex material models into finite element models of bone. This difficulty is due to variation in the spatial orientation of material properties throughout bone. Our analysis addresses this problem by taking full advantage of a finite element program's ability to solve thermal-structural problems. Using a linear relationship between temperature and modulus, we seeded specific nodes of the finite element model with temperatures. We then used thermal diffusion to propagate the modulus throughout the finite element model. Finally, we solved for the mechanical response of the finite element model to the applied loads and constraints. We found that using the thermal diffusion analogy to control the modulus of bone throughout its structure provides a simple and effective method of spatially varying modulus. Results compare favorably against both experimental data and results from an FE model that incorporated a complex (orthotropic) material model. This method presented will allow researchers the ability to easily incorporate more material property data into their finite element models in an effort to improve the model's accuracy.     

Numerical modelling of shape regulation and growth stresses in trees

This paper is the first of a series focusing on the biomechanical analysis of live trees. The finite element method (fem) is the most common method used for the analysis of complex mechanical structures. Several fem industrial codes exist, but they need to be adapted to calculate the mechanical behaviour of growing trees. A general incremental model has been developed for this specific application. In this model, time was discretised and for any developmental stage, a new equilibrium was written considering the increment of weight due to the mass of new wood layers and new vegetative elements being added. Maturation strains of new-formed cells were also considered for the simulation of the shoot reorientation process. This model was intended for use at the whole plant level. A multi-layer beam finite element is presented, which is well adapted to discretise tree limbs. The shape evolution of the structure was represented at each time step by the nodal displacement vector. The mechanical stresses induced as a result of growth were determined within the stem using a cumulative process taking into account the past history of each growth ring. The first basic results of growth stresses and shape evolution were compared with already published results at the branch level.     

Proximal Femoral Density Patterns are Consistent with Bicentric Joint Loads

We developed an alternate method for density-based load estimation and applied it to estimate hip joint load distributions for two femora. Two-dimensional finite element models were constructed from single energy quantitative computed tomography (QCT) data. Load estimation was performed using five loading regions on the femoral head. Within each loading region, individual nodal loads, normal to the local surface, were supplied as input to the load estimation. An optimization procedure independently adjusted individual nodal load magnitudes in each region, and the magnitudes of muscle forces on the greater trochanter, such that the applied tissue stimulus approached the reference stimulus throughout the model. Dominant estimated load resultant directions were generally consistent with published experimental data for loads during gait. The estimated loads also suggested that loads near the extremes of the articulating surface may be important (even required) for development and maintenance of normal bone architecture. Estimated load distributions within nearly all regions predicted bicentric loading patterns, which are consistent with observations of hip joint incongruity. Remodeling simulations with the estimated loads predicted density distributions with features qualitatively similar to the QCT data sets. This study illustrates how applications of density-based bone load estimation can improve understanding of dominant loading patterns in other bones and joints. The prediction of bicentric loading suggests a very fine level of local adaptation to details of joint loading.     

Proximal Femoral Density Patterns are Consistent with Bicentric Joint Loads

We developed an alternate method for density-based load estimation and applied it to estimate hip joint load distributions for two femora. Two-dimensional finite element models were constructed from single energy quantitative computed tomography (QCT) data. Load estimation was performed using five loading regions on the femoral head. Within each loading region, individual nodal loads, normal to the local surface, were supplied as input to the load estimation. An optimization procedure independently adjusted individual nodal load magnitudes in each region, and the magnitudes of muscle forces on the greater trochanter, such that the applied tissue stimulus approached the reference stimulus throughout the model. Dominant estimated load resultant directions were generally consistent with published experimental data for loads during gait. The estimated loads also suggested that loads near the extremes of the articulating surface may be important (even required) for development and maintenance of normal bone architecture. Estimated load distributions within nearly all regions predicted bicentric loading patterns, which are consistent with observations of hip joint incongruity. Remodeling simulations with the estimated loads predicted density distributions with features qualitatively similar to the QCT data sets. This study illustrates how applications of density-based bone load estimation can improve understanding of dominant loading patterns in other bones and joints. The prediction of bicentric loading suggests a very fine level of local adaptation to details of joint loading.     

Viscoelastic material properties of the myocardium and cardiac jelly in the looping chick heart

Accurate material properties of developing embryonic tissues are a crucial factor in studies of the mechanics of morphogenesis. In the present work, we characterize the viscoelastic material properties of the looping heart tube in the chick embryo through nonlinear finite element modeling and microindentation experiments. Both hysteresis and ramp-hold experiments were performed on the intact heart and isolated cardiac jelly (extracellular matrix). An inverse computational method was used to determine the constitutive relations for the myocardium and cardiac jelly. With both layers assumed to be quasilinear viscoelastic, material coefficients for an Ogden type strain-energy density function combined with Prony series of two terms or less were determined by fitting numerical results from a simplified model of a heart segment to experimental data. The experimental and modeling techniques can be applied generally for determining viscoelastic material properties of embryonic tissues.     

A finite-element model for the mechanical analysis of skeletal muscles

In the present paper, a finite-element model for simulating muscle mechanics is described. Based on nonlinear continuum mechanics an algorithm is proposed that includes the contractile active and passive properties of skeletal muscle. Stress in the muscle is assumed to result from the superposition of a passive and an active part. The passive properties are described by a hyperelastic constitutive material law whereas the active part depends on the fibre length, shortening velocity and an activation function. The constraint of approximate incompressibility of the muscle element is satisfied as a property of the constitutive equations. Because of the nonlinear behaviour of the material and the highly dynamical performance an incremental procedure including iterative methods is used. The advantage of the model over previous formulations is the possibility to integrate the element into an engineering standard finite-element programme ANSYS using advanced numerical tools. The model allows simulations of muscle recruitment, calculations of stress and strain distributions and predictions of muscle shape. Other possible applications are studies of the muscle architecture, the effect of inertia and impacts. First, simple examples are presented.     

Method for the design and analysis of a non-linear anisotropic lower limb prosthetic socket

A simplified computer-based structural analysis procedure has been developed for evaluating the stresses in a lower limb prosthetic socket when subjected to external loads. This technique has been tested by studying a simple linear axisymmetric shell subjected to uniform loads. It was found that the stress and strain distributions obtained are close to those determined by a finite element technique. The method that has been proposed involves an incremental procedure, which can be used to analyse the behaviour of a prosthetic socket of irregular shape made from non-linear anisotropic materials. By employing the proposed procedure to study an example, in which a simplified lower limb prosthetic socket is subject to internal pressure and frictional forces, the effects of the material non-linearity have been found to be significant.     

Biaxial incremental homeostatic elastic moduli of coronary artery: two-layer model

The detailed mechanical properties of various layers of the coronary artery are important for understanding the function of the vessel. The present article is focused on the determination of the incremental modulus in different layers and directions in the neighborhood of the in vivo state. The incremental modulus can be defined for any material subjected to a large deformation if small perturbations in strain lead to small perturbations of stresses in a linear fashion. This analysis was applied to the porcine coronary artery, which was treated as a two-layered structure consisting of an inner intima-media layer and an outer adventitia layer. We adopted a theory based on small-perturbation experiments at homeostatic conditions for determination of incremental moduli in circumferential, axial, and cross directions in the two layers. The experiments were based on inflation and axial stretch. We demonstrate that under homeostatic conditions the incremental moduli are layer- and direction dependent. The incremental modulus is highest in the circumferential direction. Furthermore, in the circumferential direction, the media is stiffer than the whole wall, which is stiffer than the adventitia. In the axial direction, the adventitia is stiffer than the intact wall, which is stiffer than the media. Hence, the coronary artery must be treated as a composite, nonisotropic body. The data acquire physiological relevance in relation to coronary artery health and disease.     

Identifying heterogeneous anisotropic properties in cerebral aneurysms: a pointwise approach

The traditional approaches of estimating heterogeneous properties in a soft tissue structure using optimization-based inverse methods often face difficulties because of the large number of unknowns to be simultaneously determined. This article proposes a new method for identifying the heterogeneous anisotropic nonlinear elastic properties in cerebral aneurysms. In this method, the local properties are determined directly from the pointwise stress–strain data, thus avoiding the need for simultaneously optimizing for the property values at all points/regions in the aneurysm. The stress distributions needed for a pointwise identification are computed using an inverse elastostatic method without invoking the material properties in question. This paradigm is tested numerically through simulated inflation tests on an image-based cerebral aneurysm sac. The wall tissue is modeled as an eight-ply laminate whose constitutive behavior is described by an anisotropic hyperelastic strain energy function containing four parameters. The parameters are assumed to vary continuously in the sac. Deformed configurations generated from forward finite element analysis are taken as input to inversely establish the parameter distributions. The delineated and the assigned distributions are in excellent agreement. A forward verification is conducted by comparing the displacement solutions obtained from the delineated and the assigned material parameters at a different pressure. The deviations in nodal displacements are found to be within 0.2% in most part of the sac. The study highlights some distinct features of the proposed method, and demonstrates the feasibility of organ level identification of the distributive anisotropic nonlinear properties in cerebral aneurysms.     

Leukocyte deformability: finite element modeling of large viscoelastic deformation.

An axisymmetric deformation of a viscoelastic sphere bounded by a prestressed elastic thin shell in response to external pressure is studied by a finite element method. The research is motivated by the need for understanding the passive behavior of human leukocytes (white blood cells) and interpreting extensive experimental data in terms of the mechanical properties. The cell at rest is modeled as a sphere consisting of a cortical prestressed shell with incompressible Maxwell fluid interior. A large-strain deformation theory is developed based on the proposed model. General non-linear, large strain constitutive relations for the cortical shell are derived by neglecting the bending stiffness. A representation of the constitutive equations in the form of an integral of strain history for the incompressible Maxwell interior is used in the formulation of numerical scheme. A finite element program is developed, in which a sliding boundary condition is imposed on all contact surfaces. The mathematical model developed is applied to evaluate experimental data of pipette tests and observations of blood flow.     

Connecting incremental shear modulus and Poisson's ratio of lung tissue with morphology and rheology of microstructure

The width and curvature of the collagen and elastin fiber bundles in the human pulmonary interalveolar septa and alveolar mouths are measured. The data, together with the known mechanical properties of collagen and elastin fibers, are used to derive the incremental elastic moduli of the lung tissue. The constitutive equation for small incremental stress and strain superposed on a homeostatic inflated lung is linear and isotropic, and characterized by two material constants.     

A three-dimensional inverse finite element analysis of the heel pad

Quantification of plantar tissue behavior of the heel pad is essential in developing computational models for predictive analysis of preventive treatment options such as footwear for patients with diabetes. Simulation based studies in the past have generally adopted heel pad properties from the literature, in return using heel-specific geometry with material properties of a different heel. In exceptional cases, patient-specific material characterization was performed with simplified two-dimensional models, without further evaluation of a heel-specific response under different loading conditions. The aim of this study was to conduct an inverse finite element analysis of the heel in order to calculate heel-specific material properties in situ. Multidimensional experimental data available from a previous cadaver study by Erdemir et al. ("An Elaborate Data Set Characterizing the Mechanical Response of the Foot," ASME J. Biomech. Eng., 131(9), pp. 094502) was used for model development, optimization, and evaluation of material properties. A specimen-specific three-dimensional finite element representation was developed. Heel pad material properties were determined using inverse finite element analysis by fitting the model behavior to the experimental data. Compression dominant loading, applied using a spherical indenter, was used for optimization of the material properties. The optimized material properties were evaluated through simulations representative of a combined loading scenario (compression and anterior-posterior shear) with a spherical indenter and also of a compression dominant loading applied using an elevated platform. Optimized heel pad material coefficients were 0.001084 MPa (μ), 9.780 (α) (with an effective Poisson's ratio (ν) of 0.475), for a first-order nearly incompressible Ogden material model. The model predicted structural response of the heel pad was in good agreement for both the optimization (<1.05% maximum tool force, 0.9% maximum tool displacement) and validation cases (6.5% maximum tool force, 15% maximum tool displacement). The inverse analysis successfully predicted the material properties for the given specimen-specific heel pad using the experimental data for the specimen. The modeling framework and results can be used for accurate predictions of the three-dimensional interaction of the heel pad with its surroundings.     

Probabilistic finite element analysis of the uncemented hip replacement—effect of femur characteristics and implant design geometry

In the present study, a probabilistic finite element tool was assessed using an uncemented total hip replacement model. Fully bonded and frictional interfaces were investigated for combinations of three proximal femurs and two implant designs, the Proxima short stem and the IPS hip stem prostheses. The Monte Carlo method was used with two performance indicators: the percentage of bone volume that exceeded specified strain limits and the maximum nodal micromotion. The six degrees of freedom of bone-implant relative position, magnitude of the hip contact force (L), and spatial direction of L were the random variables. The distal portion of the proximal femurs was completely constrained and some of the main muscle forces acting in the hip were applied. The coefficients of the linear approximation between the random variables and the output were used as the sensitivity values. In all cases, bone-implant position related parameters were the most sensitive parameters. The results varied depending on the femur, the implant design and the interface conditions. Values of maximum nodal micromotion agreed with results from previous studies, confirming the robustness of the implemented computational tool. It was demonstrated that results from a single model study should not be generalised to the entire population of femurs and that bone variability is an important factor that should be investigated in such analyses.     

Multivariate analysis of metastasis risk in laryngeal carcinoma. II. Immune response

The presence of inflammatory reaction, plasma cells, and eosinophils in peritumoral connective tissue and in neoplastic stroma was evaluated with morphometrical method in 181 patients affected by laryngeal carcinoma. A logistic multiple regression model was applied making it with the use of an independent variable represented by the "infiltrating" or "expansive" types of tumor growth, in order to evaluate the probability of nodal metastatsis of each parameter. The results suggest an inverse correlationship between plasma cells and inflammatory infiltration and incidence of nodal metastatsis only in the comparison of the extreme conditions: those with scarce infiltration versus the ones with large infiltration. Inflammatory or plasmacellular infiltration may represent both a defense mechanism against cancer and an aspecific or allergic reaction. The eosinophilic infiltration shows no value in the prevention of nodal involvement.     

Bone Load Estimation for the Proximal Femur Using Single Energy Quantitative CT Data

A density-based load estimation method was applied to determine femoral load patterns. Two-dimensional finite element models were constructed using single energy quantitative computed tomography (QCT) data from two femora. Basic load cases included parabolic pressure joint loads and constant tractions on the greater trochanter. An optimization procedure adjusted magnitudes of the basic load cases, such that the applied mechanical stimulus approached the ideal stimulus throughout each model. Dominant estimated load directions were generally consistent with published experimental data for gait. Other estimated loads suggested that loads at extreme joint orientations may be important to maintenance of bone structure. Remodeling simulations with the estimated loads produced density distributions qualitatively similar to the QCT data sets. Average nodal density errors between QCT data and predictions were 0.24 g/cm(3) and 0.28 g/cm(3). The results indicate that density-based load estimation could improve understanding of loading patterns on bones.