草地利用方式对温性典型草原优势种植物功能性状的影响

当外界环境发生变化后植物能够改变自身功能性状及时调整适应策略, 因此植物功能性状能够有效地反映植物对草地利用变化的响应, 然而在内蒙古草原从植物功能性状角度开展草地利用方式影响的研究略少。该研究以内蒙古典型草原大针茅(Stipa grandis)、羊草(Leymus chinensis)、糙隐子草(Cleistogenes squarrosa)和冷蒿(Artemisia frigida) 4种主要优势种为研究对象, 探讨在长期自由放牧、割草、短期围封和长期无干扰的影响下优势种植物功能性状的差异, 以期从功能性状视角, 揭示植物在受到外界干扰后的适应策略, 旨为天然草地的可持续管理提供基础数据支持和科学依据。结果表明: 1)除糙隐子草外, 在长期放牧后内蒙古典型草原优势植物植株高度、根长和植物碳氮含量降低, 这些性状的变化能够使植物个体小型化, 适口性降低, 表明植物通过逃避放牧的策略适应长期自由放牧的干扰; 在割草管理方式下, 优势种的高度和比叶面积有增加的趋势, 其中冷蒿的氮含量对割草响应最敏感, 其根、茎、叶中的氮含量均在割草样地最低; 围封和长期无干扰处理下植物的碳氮含量增加, 表明在干扰强度降低后, 植物通过功能性状的改变从资源获取策略向资源储藏策略转变。2)对优势种功能性状集合分析表明, 糙隐子草具有较低的植株高度和较高的比叶面积, 冷蒿具有较高的木质素含量和氮含量, 这些性状能够使两种植物被家畜采食量减少, 并保证其具有较强的再生能力, 这可能是糙隐子草和冷蒿耐牧的原因; 大针茅具有最高的植株高度、最大的叶片干物质含量, 以及最高的茎、叶纤维素含量, 说明大针茅是非常典型的竞争物种, 在干扰较低的条件下, 大针茅采取竞争策略对其他物种产生较大的竞争压力可能是其占优势的重要原因。     

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³ and 0·28 g/cm³. The results indicate that density-based load estimation could improve understanding of loading patterns on bones.     

Finite Element Analysis of Tibial Implants — Effect of Fixation Design and Friction Model

Abstract

A three dimensional nonlinear finite element model was developed to investigate tibial fixation designs and friction models (Coulomb's vs nonlinear) in total knee arthroplasty in the immediate postoperative period with no biological attachment. Bi-directional measurement-based nonlinear friction constitutive equations were used for the bone-porous coated implant interface. Friction properties between the polyethylene and femoral components were measured for this study. Linear elastic isotropic but heterogeneous mechanical properties taken from literature were considered for the bone. The Tensile behaviour of polyethylene was measured and subsequently modeled by an elasto-plastic model. Based on the earlier finite element and experimental pull-out studies, pegs and screws were also realistically modeled. The geometry of every component was obtained through measurement. The PCA tibial baseplate with three different configurations was considered; one with three screws, one with one screw and two short inclined porous-coated pegs, and a third one with no fixation for the sake of comparison. The axial load of 2000N was applied through the femoral component on the medial plateau of articular insert. It was found that Coulomb's friction significantly underestimates the relative micromotion at the bone-implant interface. The lowest micromotion and lift-off were found for the design with screws. Relative micromotion and stress transfer at the bone-implant interface depended significantly on the friction model and on the baseplate anchorage configuration. Cortical and cancellous bones carried, respectively, 10–13% and 65–86% of the axial load depending on the fixation configuration used. The remaining portion was transmitted as shear force by screws and pegs. Normal and Mises stresses as well as contact area in the polyethylene insert were nearly independent of the baseplate fixation design. The Maximum Mises stress in the polyethylene exceeded yield and was found 1–2 mm below the contact surface for all designs.     

Adaptive Finite Element Analysis of the Anisotropic Biphasic Theory of Tissue-Equivalent Mechanics

Abstract

The nonlinear partial differential equations of the anisotropic biphasic theory of tissue-equivalent mechanics are solved with axial symmetry by an adaptive finite element system. The adaptive procedure operates within a method-of-lines framework using finite elements in space and backward difference software in time. Spatial meshes are automatically refined, coarsened, and relocated in response to error indications and material deformation. Problems with arbitrarily complex two-dimensional regions may be addressed. With meshes graded in high-error regions, the adaptive solutions have fewer degrees of freedom than solutions with comparable accuracy obtained on fixed quasi-uniform meshes. The adaptive software is used to address problems involving an isometric cell traction assay, where a cylindrical tissue equivalent is adhered at its end to fixed circular platens; a prototypical bioartificial artery; and a novel configuration that is intended as an initial step in a study to determine bioartificial arteries having optimal collagen and cell concentrations.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Patient-specific modelling of the foot: automated hexahedral meshing of the bones

Subject-specific finite element modelling is a powerful tool for carrying out controlled investigations of the effects of geometric and material property differences on performance and injury risk. Unfortunately, the creation of suitable meshes for these models is a challenging and time-intensive task. This paper presents an automated method of generating fully hexahedral meshes of the bones of the feet which requires only surface representations as inputs. The method is outlined and example meshes, using two human feet and the foot of a Japanese macaque, are given to demonstrate its flexibility. Mesh quality is also evaluated for the calcaneus, first metatarsal, navicular and talus. Streamlining the generation of finite element meshes of the foot will ease investigations into the patient-specific biomechanics of injury.     

Three-Dimensional Trabecular Alignment Model

The Shear-slip Mesh Update Method (SSMUM) is being used in flow simulations involving large but regular displacements of one or more boundaries of the computational domain. We follow up the earlier discussion of the method with notes on practical implementation aspects. In order to establish a benchmark problem for this class of flow problems, we define and report results from a two-dimensional viscous flow around a rotating stirrer in a square chamber. The application potential of the method is demonstrated in the context of biomedical design problem, as we perform an analysis of blood flow in a centrifugal left ventricular assist device, or blood pump, which involves a rotating impeller in a non-axisymmetric housing.     

The influence of uncemented femoral stem length and design on its primary stability: a finite element analysis

One of the crucial factors for short- and long-term clinical success of total hip arthroplasty cementless implants is primary stability. Indeed, motion at the bone–implant interface above 40 μm leads to partial bone ingrowth, while motion exceeding 150 μm completely inhibits bone ingrowth. The aim of this study was to investigate the effect of two cementless femoral stem designs with different lengths on the primary stability. A finite element model of a composite Sawbones^® fourth generation, implanted with five lengths of the straight prosthesis design and four lengths of the curved prosthesis design, was loaded with hip joint and abductor forces representing two physiological activities: fast walking and stair climbing. We found that reducing the straight stem length from 146 to 54 mm increased the average micromotion from 17 to 52 μm during fast walking, while the peak value increased from 42 to 104 μm. With the curved stem, reducing length from 105 to 54 mm increased the average micromotion from 10 to 29 μm, while the peak value increased from 37 to 101 μm. Similar findings are obtained for stair climbing for both stems. Although the present study showed that femoral stem length as well as stem design directly influences its primary stability, for the two femoral stems tested, length could be reduced substantially without compromising the primary stability. With the aim of minimising surgical invasiveness, newer femoral stem design and currently well performing stems might be used with a reduced length without compromising primary stability and hence, long-term survivorship.     

Numerical investigations of the mechanical properties of a braided non-vascular stent design using finite element method

This paper discusses various issues relating to the mechanical properties of a braided non-vascular stent made of a Ni–Ti alloy. The design of the stent is a major factor which determines its reliability after implantation into a stenosed non-vascular cavity. This paper presents the effect of the main structural parameters on the mechanical properties of braided stents. A parametric analysis of a commercial stent model is developed using the commercial finite element code ANSYS. As a consequence of the analytical results that the pitch of wire has a greater effect than other structural parameters, a new design of a variable pitch stent is presented to improve mechanical properties of these braided stents. The effect of structural parameters on mechanical properties is compared for both stent models: constant and variable pitches. When the pitches of the left and right quarters of the stent are 50% larger and 100% larger than that of the central portion, respectively, the radial stiffness in the central portion increases by 10% and 38.8%, while the radial stiffness at the end portions decreases by 128% and 164.7%, the axial elongation by 25.6% and 56.6% and the bending deflection by 3.96% and 10.15%. It has been demonstrated by finite element analysis that the variable pitch stent can better meet the clinical requirements.