Kiran H. Shivanna

Cartilage contact pressure elevations in dysplastic hips: a chronic overload model

BackgroundDevelopmental dysplasia of the hip (DDH) is a condition in which bone growth irregularities subject articular cartilage to higher mechanical stresses, increase susceptibility to subluxation, and elevate the risk of early osteoarthritis. Study objectives were to calculate three-dimensional cartilage contact stresses and to examine increases of accumulated pressure exposure over a gait cycle that may initiate the osteoarthritic process in the human hip, in the absence of trauma or surgical intervention.MethodsPatient-specific, non-linear, contact finite element models, constructed from computed tomography arthrograms using a custom-built meshing program, were subjected to normal gait cycle loads.ResultsPeak contact pressures for dysplastic and asymptomatic hips ranged from 3.56 – 9.88 MPa. Spatially discriminatory cumulative contact pressures ranged from 2.45 – 6.62 MPa per gait cycle. Chronic over-pressure doses, for 2 million cycles per year over 20 years, ranged from 0.463 – 5.85 MPa-years using a 2-MPa damage threshold.ConclusionThere were significant differences between the normal control and the asymptomatic hips, and a trend towards significance between the asymptomatic and symptomatic hips of patients afflicted with developmental dysplasia of the hip. The magnitudes of peak cumulative contact pressure differed between apposed articular surfaces. Bone irregularities caused localized pressure elevations and an upward trend between chronic over-pressure exposure and increasing Severin classification.

IA-FEMesh: An open-source, interactive, multiblock approach to anatomic finite element model development

Finite element (FE) analysis is a valuable tool in musculoskeletal research. The demands associated with mesh development, however, often prove daunting. In an effort to facilitate anatomic FE model development we have developed an open-source software toolkit (IA-FEMesh). IA-FEMesh employs a multiblock meshing scheme aimed at hexahedral mesh generation. An emphasis has been placed on making the tools interactive, in an effort to create a user friendly environment. The goal is to provide an efficient and reliable method for model development, visualization, and mesh quality evaluation. While these tools have been developed, initially, in the context of skeletal structures they can be applied to countless applications.

Validation of a C2–C7 cervical spine finite element model using specimen-specific flexibility data

This study presents a specimen-specific C2-C7 cervical spine finite element model that was developed using multiblock meshing techniques. The model was validated using in-house experimental flexibility data obtained from the cadaveric specimen used for mesh development. The C2-C7 specimen was subjected to pure continuous moments up to +/-1.0 N m in flexion, extension, lateral bending, and axial rotation, and the motions at each level were obtained. Additionally, the specimen was divided into C2-C3, C4-C5, and C6-C7 functional spinal units (FSUs) which were tested in the intact state as well as after sequential removal of the interspinous, ligamentum flavum, and capsular ligaments. The finite element model was initially assigned baseline material properties based on the literature, but was calibrated using the experimental motion data which was obtained in-house, while utlizing the ranges of material property values as reported in the literature. The calibrated model provided good agreement with the nonlinear experimental loading curves, and can be used to further study the response of the cervical spine to various biomechanical investigations.

An interactive multiblock approach to meshing the spine

Finite element (FE) analysis is a useful tool to study spine biomechanics as a complement to laboratory-driven experimental studies. Although individualized models have the potential to yield clinically relevant results, the demands associated with modeling the geometric complexity of the spine often limit its utility. Existing spine FE models share similar characteristics and are often based on similar assumptions, but vary in geometric fidelity due to the mesh generation techniques that were used. Using existing multiblock techniques, we propose mesh generation methods that ease the effort and reduce the time required to create subject-specific allhexahedral finite element models of the spine. We have demonstrated the meshing techniques by creating a C4-C5 functional spinal unit and validated it by comparing the resultant motions and vertebral strains with data reported in the literature.

An Analytical Framework for Quadrilateral Surface Mesh Improvement with an Underlying Triangulated Surface Definition

Surface mesh quality plays a very important role in the solution accuracy and in the quality of the ensuing volumetric mesh. Amongst the techniques available, optimization based methods are commonly used for mesh quality improvement. Optimization methods for volumetric and planar surface mesh quality improvement are very well researched. However, this is not true for non-planar meshes. In this manuscript, we focus on quadrilateral non-planar surface meshes obtained during hexahedral mesh generation of anatomic structures. A modified untangling function based on node normals for quadrilateral elements is proposed. A parameterization-based method available is enhanced by giving it an analytical framework. A new projection-based method is proposed and its performance is comparable to the parametric method. The results of the enhanced/proposed methods are superior to the results obtained from Laplacian smoothing.

Feature-based multiblock finite element mesh generation

Hexahedral finite element mesh development for anatomic structures and biomedical implants can be cumbersome. Moreover, using traditional meshing techniques, detailed features may be inadequately captured. In this paper, we describe methodologies to handle multi-feature datasets (i.e., feature edges and surfaces). Coupling multi-feature information with multiblock meshing techniques has enabled anatomic structures, as well as orthopaedic implants, to be readily meshed. Moreover, the projection process, node and element set creation are automated, thus reducing the user interaction during model development. To improve the mesh quality, Laplacian- and optimization-based mesh improvement algorithms have been adapted to the multi-feature datasets.

Toward the Development of Virtual Surgical Tools to Aid Orthopaedic FE Analyses

Computational models of joint anatomy and function provide a means for biomechanists, physicians, and physical therapists to understand the effects of repetitive motion, acute injury, and degenerative diseases. Finite element models, for example, may be used to predict the outcome of a surgical intervention or to improve the design of prosthetic implants. Countless models have been developed over the years to address a myriad of orthopaedic procedures. Unfortunately, few studies have incorporated patient-specific models. Historically, baseline anatomic models have been used due to the demands associated with model development. Moreover, surgical simulations impose additional modeling challenges. Current meshing practices do not readily accommodate the inclusion of implants. Our goal is to develop a suite of tools (virtual instruments and guides) which enable surgical procedures to be readily simulated and to facilitate the development of all-hexahedral finite element mesh definitions.

Surgically oriented measurements for three-dimensional characterization of tunnel placement in anterior cruciate ligament reconstruction

Objective: To develop and evaluate the feasibility and reliability of an alternative three-dimensional (3D) measurement system capable of characterizing tunnel position and orientation in ACL reconstructed knees. Methods: We developed a surgically oriented 3D measurement system for characterizing femoral and tibial drill tunnels from ACL reconstructions. This is accomplished by simulating the positioning of the drill bit originally used to create the tunnels within the bone, which allows for angular and spatial descriptions along defined axes that are established with respect to previously described anatomic landmarks and radiographic views. Computer-generated digital phantoms composed of simplified geometries were used to verify proper calculation of angular and spatial measurements. We also evaluated the inter-observer reliability of the measurements using 10 surfaces generated from cadaveric knees in which ACL tunnels were drilled. The reliability of the measurements was evaluated by intraclass correlation coefficients. Results: The digital phantom evaluation verified the measurement methods by computing angular and spatial values that matched the known values in all cases. The intraclass correlation coefficient was calculated for four users and was found to range from 0.95 to 0.99 for the femoral and tibial measurements, demonstrating near-perfect agreement. Conclusions: The characterization of ACL tunnels has historically concentrated on two-dimensional (2D) measurements; however, it can be difficult to define ACL tunnel placement using 2D methods. We have presented novel techniques for defining graft tunnel placement from 3D surface representations of the ACL reconstructed knee. These measurements provide exact tunnel location spatially and along axes that offer the potential to comparatively analyze ACL reconstructions post-operatively using advanced imaging. These methods are reliable, and have been demonstrated to be applicable to multiple single-bundle techniques for ACL reconstruction.

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 surfaces 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.