2D calculation method based on composite beam theory for the determination of local homogenised stiffnesses of long bones

A calculation method using the finite element technique is presented. Its main objective was to determine strains, stresses and more particularly stiffnesses in any cross section of a tibia, thus enabling the localisation of tibial torsion in vivo. Each tibial cross section was considered to be a non-uniform cross section of a composite beam with arbitrary orientation of fibres. The determination of stresses, strains and stiffnesses within a composite beam cross section has been defined by solving a variational problem. The validation of this method was performed on a tibial diaphysis of which each cross section was assumed to be the cross section of a composite beam made of orthotropic materials with orthotropic axes of any orientation with respect to the principal axis of the bone. The comparison of the results, from our model and that of a three-dimensional one, was performed on each nodal value (strains, stresses) of the meshed cross section as it was impossible to obtain local stiffnesses by experimentation. The good agreement between the results has validated our finite element program. Actually, this method has enabled to treat directly 2D geometric reconstructions from CT scan images with a good accuracy to determine locally the homogenised mechanical characteristics of human tibia in vivo, and particularly to quantify torsional tibial abnormalities of children without approximation of the shape of the cross section and by calculating the real moment of inertia J. The importance of the fibre orientation with regards to the stiffness values has been emphasised. This 2D method has also allowed to reduce CPU time of the 3D modelling and calculation.