iDocSlide.Com

Free Online Documents. Like!

All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.

Share

Description

Three-dimensional nonlinear finite element model of lumbar intervertebral disc

Tags

Transcript

Acta of Bioengineering and Biomechanics Vol. 7, No. 2, 2005
Three-dimensional nonlinear finite element model of lumbar intervertebral disc
T
OMASZ
Ł
ODYGOWSKI
,
W
ITOLD
K
Ą
KOL
,
M
ARCIN
W
IERSZYCKI
Institute of Structural Engineering, Pozna
University of Technology, Piotrowo 5, 60-965 Pozna
, Poland. Tomasz.Lodygowski@put.poznan.pl, Witold.Kakol@put.poznan.pl, Marcin.Wierszycki@put.poznan.pl
M
A
Ł
GORZATA
B.
O
GURKOWSKA
Department of Biomechanics, University School of Physical Education, Królowej Jadwigi 27/35, 61-871 Pozna
, Poland.
ogurkmal@man.poznan.pl The main objective of this study is to model a three-dimensional nonlinear finite-element vertebra disc, which can be used in further simulation of human lumbar spinal segments in surgery, analyses of spinal equilibrium and stability. Because of complexity of modelling it is proposed to carry out analyses on the simplified model built in such a way that the nonlinear response of the disc is replaced by a spring-type behaviour whose characteristics are obtained by computer simulations of an isolated disc. The geometry data of a human lumbar spinal segment, including a disc, is acquired from the computer tomography or magnetic nuclear resonance measurements, and a CAD model is designed and imported into FEA program. The simplified model was validated for loading schemes, including axial compression, bending and torsion acting on the spinal L4-L5 segment. Two models of intervertebra disc are shown and theirs advantages and accuracy are discussed.
Key words: biomechanics, nonlinear FEA, lumbar spine
1. Introduction
The complexity of carrying out numerical simulations of human lumbar spinal segment is rather obvious [2]–[4]. It arises from complexity of its geometry, difficulties in defining material properties (bones, ligaments etc.) and simplifications which have to be done according to recognition of boundary conditions. The
T.
Ł
ODYGOWSKI
et al. 30
experimental verification made on the living system is not easy either; it has rather qualitative and integrated character. The intervertebral disc and ligaments play the essential role in the spine kinematics and contribute to general stiffness of spinal segment [1]. It exhibits pronounced time-dependent deformations when subjected to load variations. These deformations are caused by fluid that flows to and from the disc and by viscoelastic deformation of annulus fibrous. The fluid flow is due to the differences between mechanical and osmotic pressure. The intervertebral disc is composed of three different parts. The essential component of the intervertebral disc is the annulus fibrous. It consists of 10–20 sheets of collagen called a lamellar collagen. Superior and inferior vertebral endplates are the other components of the intervertebral disc. These are the plates of cartilage that covers the superior and inferior aspects of the disc and binds the disc to their respective vertebral bodies. Besides enabling bending movements between vertebral bodies, the intervertebral disc allows twisting and small sliding movements. Their amplitudes depend on elasticity and tensile stiffness of the annulus.
2. The aim of the study
The aim of the study was to develop and validate the simplified model of the intervertebral disc and further the whole lumbar spinal motion segment. This simplification is based on replacement of the disc by a connector-type element, whose elasticity behaviour is spring-like in available components of relative motion. If nonlinear behaviour is required, a force or moment versus relative displacement or rotation can be specified. Alternatively 12 by 12 elasticity matrix (as a function of field variables) can be given. It was applied in our studies. To follow this approach and to reach the goal that leads to this crucial simplification of the spine model we have to examine first the disc (figure 2) which is subjected to six loading schemes: axial compression, shearing in two directions, bending in two planes and torsion. By applying the unit loads (axial, shear and moments) and recording displacements at the appropriate points, where relative motion is estimated, a compliance matrix and then the stiffness matrix of the disc can be built. In the range of our interest, i.e. small rotations and small displacements, a tangent stiffness matrix can be used because of the assumed linearity of load–displacement relations. In order to design a computational model of the spinal motion segment, we need some precise geometrical data of the real object. Besides topology, the additional data such as volume density, surface texture, etc. can be of interest. Different methods of acquisition of geometrical data can be used. They may be obtained with contact scanners and non-contact scanners. These methods are described in detail in [13], [16], [19]. The geometry data of a human lumbar spinal segment, including a disc, was collected using the computer tomography and magnetic nuclear resonance. The CAD
Three-dimensional nonlinear finite element model of disc
31
model was developed and imported into FEA commercial code ABAQUS and meshed. A final result of meshing is shown in figure 1.
3. FEA modelling
The structure and properties of the spinal motion segment analyzed are complex. The first step taken to allow simulation of the behaviour of biomechanical system is a proper definition of the model geometry. In the kinematics of the spinal motion segment, intervertebral disc and ligaments play the most essential role. In this study, we focused our attention on a proper modelling of disc in numerical simulation of L4– L5 segment. The FEA model of motion segment is shown in figure 1.
Fig. 1. FE model of human lumbar spinal segment: VB – vertebral bodies, ID – intervertebral disc, L – ligaments, PE – posterior elements
The annulus fibrous of the intervertebral disc is a highly structured material made up of alternating tissue layers with collagen fibers oriented at +30
°
and –30
°
with respect to the circumferential axis. The disc can be modelled as anisotropic material of hyperelastic properties [4] or by structural elements which induce this anisotropy. The ground matrix of the disc annulus is modelled by 3D solid elements. The collagen fibers can be modelled by truss elements carrying only tensile stresses [6] and by rebar type elements embedded in 3D solid elements. When using the first approach each truss element in each layer of the annulus fibrous has to be connected with 3D solid
ID PE
VB
L L
VB
T.
Ł
ODYGOWSKI
et al. 32
element nodes in such a way as to keep up the orientation mentioned above (figure 2a). Development of such a model is a quite tedious and time-consuming occupation. Additionally, this approach is mesh-dependent because any remeshing requires re-defining the data concerning truss elements (their orientation and properties). In the present study, the second method was employed which has some advantages over the previous one being mentioned above. The rebar elements are uniaxial reinforcements in solid elements which can be defined as surface layers (figure 2b) with uniformly spaced reinforcing bars. Such layers are treated as smeared layers of a constant thickness equal to the area of each reinforcing bar divided by the reinforcing bar spacing. This approach allows ease altering the number of layers, section properties and section orientation and is independent of element re-meshing. The calibration of the 3D embedded elements by means of published data [9] was performed in our study by selecting an appropriate number of fibers (rebars) with their cross-sectional areas. The nucleus pulposus, assumed here as an incompressible body, was modelled as fluid-filled cavity using hydrostatic fluid elements with initial pressure equal to 2 MPa [9]. Hydrostatic fluid elements cover the boundaries of the nucleus pulposus. They share the nodes at the cavity boundary with the standard elements of annulus fibrous.
annulus fibrosusnucleus pulposuscollagen fibers
a) b) Fig. 2. FE models of intervertebral disc: a) with truss elements, b) with surfaces or rebar layers
Selected constitutive data used in the model of a disc for further computations are summarized in table 1.
annulus fibrous
Three-dimensional nonlinear finite element model of disc
33 Table 1. Material properties of disc components Anatomic part Material Properties Endplates isotropic
E
= 23.8 MPa
ν
= 0.4 Annulus fibrous – ground matrix isotropic
E
= 4.0 MPa
ν
= 0.4 Annulus fibrous – fibers isotropic, no compression
E
= 45 MPa – Nucleus pulposus – incompressible body incompressible fluid
ρ
= 1.0E–6 kg/mm
3
– Ligaments isotropic
E
= 6–12 Mpa –
4. FE
analysis
A basic idea of the simplification of a motion segment modelling is to replace the complex structure of the intervertebral disc by one connector-type element of complex properties. The elastic behaviour of the element of this type can be described as equivalent stiffness matrix that has in general case the form of
F
=
KD
, where
F
is the vector (12 components) of generalized forces that act on the segment,
D
is the vector (12 components) of mutual displacements between bones, and
K
(12
×
12 matrix) is the stiffness operator of the segment. The values of these matrix components were obtained based on FEA simulations. For the behaviour of the segment in the range of our interest, i.e. small rotations and small displacements, a tangent stiffness matrix was used because of assumed linearity of load displacement relations. In order to follow this approach, the isolated disc, which was subjected to twelve loading schemes, was examind in the first place. By means of applying the unit loads (axial, shear and moments) and recording displacements at the appropriate points, where relative motion was estimated, a compliance matrix and then the stiffness matrix were built for the model. Four concepts of connector elements were tested. The first three were based on a two-node, twelve dof (six transverse and six rotational) element whose behaviour was described by means of 12 by 12 stiffness matrix. Three definitions of this matrix were prepared: the full matrix with 144 non-zero components, the reduced matrix with 92 significant non-zero components and the symmetrized matrix with 52 non-zero components. The last one was based on a special type of connector element. It was a two-node element with 6 by 6 stiffness matrix. This matrix described relative movements and rotations of these two nodes (table 2). The last stage was the comparison of the numerical models of motion segments containing one element of intervertebral disc with a model containing a multi-element, complex disc definition. All four concepts of the equivalent element were studied. The six load cases discussed above were taken into consideration (figure 3). The recorded relative displacements and rotations allowed validation of different concepts of intervertebral disc simplifications. The loads were applied to the reference point which coupled all points lying on the upper surface of the disc, while the bottom surface of the disc was protected against

We Need Your Support

Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

SAVE OUR EARTH

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...Sign Now!

We are very appreciated for your Prompt Action!

x