^{1}

^{*}

^{1}

^{2}

Many researches point out that intervertebral pressure and transformation are key parameters for evaluating intervertebral disc degeneration. Aiming at avoiding the damage caused by direct and indirect measuring methods, this research proposes a cylindroid hypothesis and measuring method, which can monitor the strain condition of the intervertebral disc in vivo and real-time without being damaged.

According to incomplete statistical analysis, around 7% - 10% of the 1.3 billion people in China have some kind of cervical disease. It shows a growing trend. After some kind of treatment, the motion function of the fused segments is lost forever and the adjacent level becomes degenerate easily [

Based on the physiology structure and properties of the intervertebral disc, this study raised a solution that could monitor the intervertebral strain in vivo without being damaged. It can solve the intervertebral disc strain measure dilemma that direct measure methods can only be executed on a few healthy volunteers and other simulation or indirect measures are too hard to achieve an agreement.

Aiming at getting the strain condition of the intervertebral disc, a mathematic model and a measurement system based on strain gauges was built. In our hypothesis, the complete strain condition and the dimension of the disc can be calculated by only fixing 6 strain gauges on the front side of an intervertebral disc.

In some imaging measurement about intervertebral disc, we can easily find out that the transverse section of the disc is an axialsymmetrypolygonwith a smooth border. At first, the study assumes the geometric shape of the intervertebral disc as a cylindroid, sets the origin at the geometric center and establishes a coordinate in the vertical view. As

According to the result of actual measurement, the length in the coronal plane is longer than the length in the sagittal plane, the expression of the side face of the disc in Cartesian coordinates is:

A is the long axis of the ellipse, B is the short axis of the ellipse, and A > B. As the dimension of the ellipse can’t be measured in vivo, A and B are both unknown. With the help of the transfer equation between Cartesian coordinates and polar coordinates:

We can rewrite (1) as:

Link a certain point on the ellipse’s angle with the origin, the angle between the X axis and it is, is the distance between from the point to the origin. Any point on the ellipse can be described as.

In order to realize the nondestructive measurement, only a bit of arc symmetric with Y axis is available.

As Figures 2 and 3 show, 6 strain gauges were packaged in flexible PCB and fixed on the front surface of the intervertebral disc, symmetrical with YZ plane. ab = de and bc = cd.

a, b, c, d, e are tensile strain gauges and sensitive for the tension in Z direction. f is a shear strain gauge and sensitive for the in torsion Z direction. The measuring result of f is independent to others’ equation, we focus on the geometric equation and strain equation of a, b, c, d, e.

After fixing the flexible PCB right in the front of the intervertebral disc, there are only 4 independent parameters:

The location of a, b, c, d, e are respective:,

, , ,.

We can get 2 independent equations:

Add the arc length formula of plane curve:

We can get:

In the equation above: and are known quantity, are unknown quantity.

In the normal stress analysis theory, the entire deformation of the intervertebral disc can be resolved into 4 kinds:

1) Bending deformation in X direction;

2) Bending deformation in Y direction;

3) Bending deformation in Z direction;

4) Torsional deformation in Z direction.

According to the equation of strain compatibility, the strain at a, b, c, d, e can be described as below:

In summary, the model including 7 known quantities:, details are shown in

9 unknown quantities:, details are shown in

9 independent equations:

After solving the equation system above, we can bring B and into the shear strain equation:

Then we can get the parameter describing torsional strain in Z direction:.

For convenience, the equation system can be described as below:

Because the elliptic integral exists in Equations (18) and (19), the equation system is hard to be solved directly. So we assume the long axis radius, short axis radius and angle of the measure point at first and then select the fitting answer. The specific steps were presented in

The model above is theoretical and hypothesis based. Some modification and improvement will make the model more effective and accurate.

Firstly, the error the cylindroid simply introduced into the model needs to be calculated and modificated. Some imaging research tells that the front height of the intervertebral disc is longer than the back height, so an angle about 3˚ - 6˚ exists between the top and bottom surface of the intervertebral disc naturally.

Secondly, the strain gauge should stick on the surface directly. But confined by the biocompatibility and operation time, the strain gauge had to be packed in a biocompatibility flexible board and sticked at the same time. The

board’s properties especially the thickness will bring some error to the result.

Thirdly, it is hard to modify the model in vivo because of the unavoidable specimen specificity, installation error, model simplification error and follow-up circuit error. It will also lose the real-time advantage. A more effective method is to practice the experiment on cadaveric cervical with both the ellipsoid hypothesis method and the direct pressure measure method. The data of ellipsoid hypothesis method will be adjusted by taking the data from pressure measurement method into consideration. It will enable this new model more convincing when it was taken in vivo.

Last but not least, comparison study with other intervertebral pressure or strain measurement method will make our model more persuasive. Before the whole model and hypothesis can be taken into practice, a good number of experimental data on cadaveric cervical should be studied and parameter in the model should be adjusted.

Aiming at getting the strain condition of the intervertebral disc, a mathematic model and a measurement system based on strain gauges were built. In our hypothesis, by only fixing 6 strain gauges on the front side of anintervertebral disc, the complete strain condition and the dimension of the disc can be calculated.

Compared with the former method, the model in this paper is without being damaged and can be operated on the adjacent disc incidentally when a treatment is executed on a patient. It can solve the intervertebral disc strain measure dilemma that direct measure methods can only be executed on a few healthy volunteers and other simulation or indirect measure is too hard to achieve an agreement.

So far, the hypothesis is in exploring stage. We will move to experiment stage and further modify and adjust the model accordingly.