Gaussian curvature analysis allows for automatic block placement in multi-block hexahedral meshing

Musculoskeletal finite element analysis (FEA) has been essential to research in orthopaedic biomechanics. The generation of a volumetric mesh is often the most challenging step in a FEA. Hexahedral meshing tools that are based on a multi-block approach rely on the manual placement of building blocks for their mesh generation scheme. We hypothesise that Gaussian curvature analysis could be used to automatically develop a building block structure for multi-block hexahedral mesh generation. The Automated Building Block Algorithm incorporates principles from differential geometry, combinatorics, statistical analysis and computer science to automatically generate a building block structure to represent a given surface without prior information. We have applied this algorithm to 29 bones of varying geometries and successfully generated a usable mesh in all cases. This work represents a significant advancement in automating the definition of building blocks.     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

Hexahedral meshing of subject-specific anatomic structures using mapped building blocks

To extend the use of computational techniques like finite element analysis to clinical settings, it would be beneficial to have the ability to generate a unique model for every subject quickly and efficiently. This work is an extension of two previously developed mapped meshing tools that utilised force and displacement control to map a template mesh to a subject-specific surface. The objective of this study was to map a template block structure, common to multiblock meshing techniques, to a subject-specific surface. The rationale is that the blocks are considerably less refined and may be readily edited after mapping, thereby yielding a mesh of high quality in less time than mapping the mesh itself. In this paper, the versatility and robustness of the method was verified by processing four data-sets. The method was found to be robust enough to cope with the variability of bony surface size, spatial position and geometry, producing building block structures (BBSs) that generated meshes comparable to those produced using BBSs that were created manually.     

Patient-specific computational haemodynamics: generation of structured and conformal hexahedral meshes from triangulated surfaces of vascular bifurcations

Measuring the blood flow is still limited by current imaging technologies and is generally overcome using computational fluid dynamics (CFD) which, because of the complex geometry of blood vessels, has widely relied on tetrahedral meshes. Hexahedral meshes offer more accurate results with lower-density meshes and faster computation as compared to tetrahedral meshes, but their use is limited by the far more complex mesh generation. We present a robust methodology for conformal and structured hexahedral mesh generation – applicable to complex arterial geometries as bifurcating vessels – starting from triangulated surfaces. Cutting planes are used to slice the lumen surface and to construct longitudinal Bezier splines. Afterwards, an isoparametric transformation is used to map a parametrically defined quadrilateral surface mesh into the vessel volume, resulting in stacks of sections which can then be used for sweeping. Being robust and open source based, this methodology may improve the current standard in patient-specific mesh generation and enhance the reliability of CFD to patient-specific haemodynamics.     

Biomechanical analysis of foot with different foot arch heights: a finite element analysis

Clinically, different foot arch heights are associated with different tissue injuries to the foot. To investigate the possible factors contributing to the difference in foot arch heights, previous studies have mostly measured foot pressure in either low-arched or high-arched feet. However, little information exists on stress variation inside the foot with different arch heights. Therefore, this study aimed to implement the finite element (FE) method to analyse the influence of different foot arches. This study established a 3D foot FE model using software ANSYS 11.0. After validating the FE model, this study created low-arched, high-arched and normal-arched foot FE models. The FE analysis found that both the stress and strain on the plantar fascia and metatarsal were higher in the high-arched foot, whereas the stress and strain on the calcaneous, navicular and cuboid were higher in low-arched foot. Additionally, forefoot pressure was increased with an increase in arch height.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Unstructured hexahedral mesh generation of complex vascular trees using a multi-block grid-based approach

The trend towards realistic numerical models of (pathologic) patient-specific vascular structures brings along larger computational domains and more complex geometries, increasing both the computation time and the operator time. Hexahedral grids effectively lower the computational run time and the required computational infrastructure, but at high cost in terms of operator time and minimal cell quality, especially when the computational analyses are targeting complex geometries such as aneurysm necks, severe stenoses and bifurcations. Moreover, such grids generally do not allow local refinements. As an attempt to overcome these limitations, a novel approach to hexahedral meshing is proposed in this paper, which combines the automated generation of multi-block structures with a grid-based method. The robustness of the novel approach is tested on common complex geometries, such as tree-like structures (including trifurcations), stenoses, and aneurysms. Additionally, the performance of the generated grid is assessed using two numerical examples. In the first example, a grid sensitivity analysis is performed for blood flow simulated in an abdominal mouse aorta and compared to tetrahedral grids with a prismatic boundary layer. In the second example, the fluid–structure interaction in a model of an aorta with aortic coarctation is simulated and the effect of local grid refinement is analyzed.     

Parametric study of orthopedic insole of valgus foot on partial foot amputation

Orthopedic insole was important for partial foot amputation (PFA) to achieve foot balance and avoid foot deformity. The inapposite insole orthosis was thought to be one of the risk factors of reamputation for foot valgus patient, but biomechanical effects of internal tissues on valgus foot had not been clearly addressed. In this study, plantar pressure on heel and metatarsal regions of PFA was measured using F-Scan. The three-dimensional finite element (FE) model of partial foot evaluated different medial wedge angles (MWAs) (0.0°–10.0°) of orthopedic insole on valgus foot. The effect of orthopedic insole on the internal bone stress, the medial ligament tension of ankle, plantar fascia tension, and plantar pressure was investigated. Plantar pressure on medial heel region was about 2.5 times higher than that of lateral region based on the F-Scan measurements. FE-predicted results showed that the tension of medial ankle ligaments was the lowest, and the plantar pressure was redistributed around the heel, the first metatarsal, and the lateral longitudinal arch regions when MWA of orthopedic insole ranged from 7.5° to 8.0°. The plantar fascias maintained about 3.5% of the total load bearing on foot. However, the internal stresses from foot bones increased. The simulation in this study would provide the suggestion of guiding optimal design of orthopedic insole and therapeutic planning to pedorthist.     

Implication of obesity on motion,posture and internal stress of the foot: an experimental and finite element analysis

Obesity causes increased loading on the foot which can damage the soft tissue and bone ultimately leading to foot problems. Experimental and computational methods were used to analyse the chain of biomechanical changes in the lower limb due to obesity. The experimental study shows some changes in foot posture and gait where obese subjects were more likely to have pronated feet, smaller joint angles in the sagittal and frontal planes, smaller cadence, and smaller stride length. Anatomically correct finite element models generated on obese subjects showed increased and altered internal and plantar stress. Altered foot posture was identified as a key indicator of increased internal stress indicating the importance of foot posture correction.     

Development of a patient-specific anatomical foot model from structured light scan data

The use of anatomically accurate finite element (FE) models of the human foot in research studies has increased rapidly in recent years. Uses for FE foot models include advancing knowledge of orthotic design, shoe design, ankle–foot orthoses, pathomechanics, locomotion, plantar pressure, tissue mechanics, plantar fasciitis, joint stress and surgical interventions. Similar applications but for clinical use on a per-patient basis would also be on the rise if it were not for the high costs associated with developing patient-specific anatomical foot models. High costs arise primarily from the expense and challenges of acquiring anatomical data via magnetic resonance imaging (MRI) or computed tomography (CT) and reconstructing the three-dimensional models. The proposed solution morphs detailed anatomy from skin surface geometry and anatomical landmarks of a generic foot model (developed from CT or MRI) to surface geometry and anatomical landmarks acquired from an inexpensive structured light scan of a foot. The method yields a patient-specific anatomical foot model at a fraction of the cost of standard methods. Average error for bone surfaces was 2.53 mm for the six experiments completed. Highest accuracy occurred in the mid-foot and lowest in the forefoot due to the small, irregular bones of the toes. The method must be validated in the intended application to determine if the resulting errors are acceptable.     

Protocol for constructing subject-specific biomechanical models of knee joint

A robust protocol for building subject-specific biomechanical models of the human knee joint is proposed which uses magnetic resonance imaging, motion analysis and force platform data in conjunction with detailed 3D finite element models. The proposed protocol can be used for determining stress and strain distributions and contact kinetics in different knee elements at different body postures during various physical activities. Several examples are provided to highlight the capabilities and potential applications of the proposed protocol. This includes preliminary results on the role of body weight on the stresses and strains induced in the knee articular cartilages and meniscus during single-leg stance and calculations of the induced stresses and ligament forces during the gait cycle.     

Finite elements/Taguchi method based procedure for the identification of the geometrical parameters significantly affecting the biomechanical behavior of a lumbar disc

The aim of this work is to show a quick and simple procedure able to identify the geometrical parameters of the intervertebral disc that strongly affect the behavior of the FEM model. First, we allocated a selection criterion for the minimum number of geometrical parameters that describe, with a good degree of approximation, a healthy human vertebra. Next, we carried out a sensitivity analysis using the ‘Taguchi orthogonal array’ to arrive at a quick identification of the parameters that strongly affect the behavior of the Fem model.     

Intramyocardial Injections to De-Stiffen the Heart: A Subject-Specific in Silico Approach

We hypothesized that minimally invasive injections of a softening agent at strategic locations in stiff myocardium could de-stiffen the left ventricle (LV) globally. Physics-based finite element models of the LV were created from LV echocardiography images and pressures recorded during experiments in four swine. Results confirmed animal models of LV softening by systemic agents. Regional de-stiffening of myocardium led to global de-stiffening of LV. The mathematical set up was used to design LV global de-stiffening by regional softening of myocardium. At an end diastolic pressure of 23 mmHg, when 8 ml of the free wall was covered by intramyocardial injections, end diastolic volume (EDV) increased by 15.0%, whereas an increase up to 11 ml due to intramyocardial injections in the septum and free wall led to a 26.0% increase in EDV. Although the endocardial intramyocardial injections occupied a lower LV wall volume, they led to an EDV (44 ml) that was equal compared to intramyocardial injections in the mid-wall (44 ml) and larger compared to intramyocardial injections in the epicardium (41 ml). Using an in silico set up, sites of regional myocardium de-stiffening could be planned in order to globally soften overly stiff LV in heart failure with preserved ejection fraction. This novel treatment is built on subject-specific data. Hypothesis-testing of these simulation findings in animal models is warranted.     

Factors affecting subject-specific finite element models of implant-fitted rat bone specimens: critical analysis of a technical protocol

The authors propose a protocol to derive finite element (FE) models from micro computer tomography scans of implanted rat bone. A semi-automatic procedure allows segmenting the images using specimen-specific bone mineral density (BMD) thresholds. An open-source FE model generator processes the segmented images to a quality tetrahedral mesh. The material properties assigned to each element are integrated from the BMD field. Piecewise, threshold-dependent density–elasticity relationships are implemented to limit the effects of metal artefacts. A detailed sensitivity study highlights the coherence of the generated models and quantifies the influence of the modelling parameters on the results. Two applications of the protocol are proposed. The stiffness of bare and implanted rat tibiae specimens is predicted by simulating three-point bending and inter-implant displacement, respectively. Results are compared with experimental tests. The mean value and the variability between the specimens are well captured in both tests.     

Fluid-induced osteolysis: modelling and experiments

A model to calculate bone resorption driven by fluid flow at the bone–soft tissue interface is developed and used as a basis for computer calculations, which are compared to experiments where bone is subjected to fluid flow in a rat model. Previous models for bone remodelling calculations have been based on the state of stress, strain or energy density of the bone tissue as the stimulus for remodelling. We believe that there is experimental support for an additional pathway where an increase in the amount of the cells directly involved in bone removal, the osteoclasts, is caused by fluid pressure, flow velocity or other parameters related to fluid flow at the bone–soft tissue interface, resulting in bone resorption.     

Patient-specific simulation of a stentless aortic valve implant: the impact of fibres on leaflet performance

In some cases of aortic valve leaflet disease, the implant of a stentless biological prosthesis represents an excellent option for aortic valve replacement (AVR). In particular, if compared with the implant of mechanical valves, it provides a more physiological haemodynamic performance and a reduced thrombogeneticity, avoiding the use of anticoagulants. The clinical outcomes of AVR are strongly dependent on an appropriate choice of both prosthesis size and replacement technique, which is, at present, strictly related to surgeon's experience and skill. This represents the motivation for patient-specific finite element analysis able to virtually reproduce stentless valve implantation. With the aim of performing reliable patient-specific simulations, we remark that, on the one hand, it is not well established in the literature whether bioprosthetic leaflet tissue is isotropic or anisotropic; on the other hand, it is of fundamental importance to incorporate an accurate material model to realistically predict post-operative performance. Within this framework, using a novel computational methodology to simulate stentless valve implantation, we test the impact of using different material models on both the stress pattern and post-operative coaptation parameters (i.e. coaptation area, length and height). As expected, the simulation results suggest that the material properties of the valve leaflets affect significantly the post-operative prosthesis performance.