Biomechanical analysis of INFINITY rehabilitation method for treatment of low back pain.

[Purpose] Low back pain is a pervasive problem in modern societies. Physical rehabilitation in treatment of low back pain should reduce pain, muscle tension and restore spine stability and balance. The INFINITY® rehabilitation method that is based on a figure of eight movement pattern was proved to be effective in low back pain treatment. The aim of the paper is to estimate the effect of a figure of eight motion on the L5/S1 load and lumbar spine muscle activation in comparison to other motion patterns.
[Subjects and Methods] Three-dimensional model of lumbar spine musculoskeletal system is used to simulate effect of various load motion pattern induced by displacement of the center of gravity of the upper body. Four motion patterns were examined: lateral and oblique pendulum-like motion, elliptical motion and figure of eight motion.
[Results] The simple pendulum-like and elliptical-like patterns induce harmonic muscle activation and harmonic spinal load. The figure of eight motion pattern creates high-frequency spinal loading that activates remodeling of bones and tendons. The figure of eight pattern also requires muscle activity that differs from harmonic frequency and is more demanding on muscle control and could also improve muscle coordination.
[Conclusion] The results of the study indicate that complex motion pattern during INFINITY® rehabilitation might enhance the spine stability by influencing its passive, active and neural components.

INTRODUCTION

Statistical reports from many developed countries indicate that low back pain is very common with almost half the adult population reporting low back pain lasting for at least 24 hours at some time in the year. Moreover, low back pain is a very frequent cause of long term disability in middle age (35 to 55 years of age)1). It is proposed that the most efficient prevention against low back pain is keeping movement activity supported with a suitable rehabilitation regime that does not overload the motoric system but strengthens and stabilizes the spinal system (for review see Guzmán et al.2); Laird et al.3)). Usually, the physical therapy involves supervised exercise therapy, spinal manipulative therapy or home exercise and advice4).

A rehabilitation program should relieve pain and muscle tension; restore motion and balance; develop strength, endurance and power; as well as maintain cardiovascular fitness5). Traditional physical therapy or fitness training has emphasized strengthening of individual muscle groups. It has been proposed that the skill or the quality of the movement pattern that is used could be more important than achievement of muscle strength6). Stability and not strength can be considered as the objective in rehabilitation of low back pain7, 8). Clinical studies also indicate improvements in disability degree and/or pain level using spine stabilization approach9). Various guidelines for treating low back pain have been proposed in different countries and regions10). There is no therapy providing full recovery in all patients with low back pain yet and new therapeutic methods are being tested10, 11).

An original rehabilitation method named INFINITY® has been developed at the Rehabilitation Clinic Brandys nad Orlici, in the Czech Republic for low back pain therapy. The method has been verified in a clinical study12). Its name comes from the English word infinity because it utilizes movement in a figure of eight as a key segment of the exercises. The INFINITY® method is based on micro-movements in the range of millimeters and proved to be effective in patients with small spinal ROM (range of motion). The INFINITY® method can even be applied in the acute phase when a patient may be suffering from intense pain and it offers extended variability of autotherapeutic exercises. It was hypothesized that a figure of eight motion contributes to the stabilization and the strengthening of trunk muscles including the deep stabilization system. However, there seems to be a lack of science-based knowledge whether and how the figure of eight motion pattern during INFINITY® exercise influences spinal and muscular load.

Advances in mathematical modeling and technology are creating new opportunities that may quantify the effect of physiotherapy on lumbar spine13, 14). The overall aim of the study is to compare the effect of a figure of eight movement on deep muscle activation and spinal load with simple movements during rehabilitation. The spinal muscle and L5/S1 load is assessed using mathematical modeling.

SUBJECTS AND METHODS

Biomechanical modelling was adopted to estimate muscle and joint forces that are difficult to study through experimentation. The input into the biomechanical model is the anatomical model coupled with the physics of the musculoskeletal system. Considering the body static and dynamic equilibrium, active and passive internal forces (muscle and joint forces, respectively) acting within anatomical structure could be determined. Upright standing was taken as reference body posture within the present study and the spinal load is estimated considering equilibrium of gravitational, intra-abdominal, muscle and spinal reaction forces and torques. The L5/S1 junction was chosen as a point where spinal load is being estimated. The human body could be divided into two segments with respect to L5/S1: the part above L5/S1 junction, and the part below L5/S1 junction, i.e. the upper body and the lower body, respectively (Fig. 1). Center of gravity of the upper body is located in front of the lumbar spine in upright standing and activity of postural muscles is needed to maintain equilibrium.

Fig. 1.

Coordinate system and musculoskeletal model. Abdominal muscles are not shown. Origin of coordinate system corresponds to projection of center of gravity of upper body to L5/S1 transversal plane.

The musculoskeletal model of the lumbar spine described in Daniel13) is adopted (Fig. 1). Muscles attached over a large area were divided into individual musculoskeletal units according to Stokes and Gardner-Morse15). Muscles with curved line of action were defined point-wise using wrapping points. The model contains 284 individual muscle units15).

The number of muscles included in the model exceeds number of equilibrium equations considerably. In principle, there are an infinite number of different combinations of muscle and joint forces that can be applied to maintain a particular posture or produce a given movement16). Therefore, the muscle activity required to maintain equilibrium in a given position of body is computed using the method of inverse dynamic optimization16). This method is based on assumption that loads are shared by muscles according to the minimization of some performance criterion. Minimization of the sum of muscle stresses cubed was used in this study as a performance criterion. This optimization criterion is based on the idea that muscles are activated in a way that minimizes the energy expenditure and maximizes their endurance time17). The force of each muscle is constrained not to exceed the maximum isometric force18). After computation of the muscle forces, the components of the L4/S5 loading force are determined from the force equilibrium equations for the L5/S1 segment of spine. Magnitude of spine loading forces is used to express loading of the L5/S1 segment.

Rehabilitation motion pattern is represented by the motion of the center of gravity of the upper body with respect to its neutral position. Projection of the center of gravity in a neutral position of the body to a transversal plane passing the center of L5/S1 intervertebral disc is chosen as the center of the coordinate system (Fig. 1). Orientation of Cartesian coordinates is defined according to Stokes19): the x-axis points anteriorly, the y-axis points laterally and the z-axis points superiorly, defining four segments: anterior-left (AL), anterior-right (AR), posterior-left (PL), and posterior-right (PR) segment (Fig. 1).

Although the motion of the center of gravity of the upper body may follow a three-dimensional curve, the spine torque is influenced only by displacement of the center of gravity in the transversal (xy) plane. The periodical motion of the center of gravity of the upper body in the xy plane during the exercise can be described by kinematic equation.

(1)

where ω=2π/T is the angular frequency with the period T, f is function ranging from 0 to 1 and ax and ay are amplitudes of motion in x and y direction, respectively. As proposed by Schneck20), the motion during the physical exercises can be described using simple harmonic functions (Table 1).

Table 1.

Kinematic equations of the center of the upper body in the plane of section

Four types of motion were modeled (Fig. 2, Table 1) and were chosen on the basis of our clinical experience. In upright standing a patient actively moves the center of gravity of their upper body along a specified curve. The linear L-R motion (Fig. 2A) and linear oblique motion AL-PR (Fig. 2B) represents simple pendulum movement in frontal and oblique planes, respectively. The circular motion is represented by ellipse trajectory (Fig. 2C) starting at anterior position. The most complex trajectory has figure of eight (∞) shape (Fig. 2D), denoted as lemniscate in algebraic geometry. A figure of eight is consistent with the INFINITY® method. A modified equation of Bernoulli lemniscate is used (Table 1).

Fig. 2.

Trajectories of upper body center of gravity and its projection into the transversal plane. (A) L-R pendulum motion, (B) AL-PR pendulum motion, (C) ellipse motion, (D) ∞ shape motion. Displacements not in scale.

There exist two principal types of INFINITY motion: the macro-movement (amplitudes in a range of centimeters) and the micro-movement (amplitudes in a range of millimeters)12). The period of movement depends on the movement range: the former has the period T from 7 to 14 s, while the periord T is up to 5 s in the latter movement. Considering anteroposterior and mediolateral displacement a and a to be equal to 2 cm and 5 mm for macro and micro/movement respectively, the speed is lower than 0.05 m/s and therefore inertial components are considered to be negligible, i.e. slow motion quasi-static analysis is performed13). It is further assumed that small displacements of the center of gravity of the upper body do not significantly change musculoskeletal geometry of the lower lumbar spine, the motion occurs mainly in upper spine segments with larger mobility. Based on the INFINITY® method guidelines, the maximum anteroposterior (AP) displacement (a) and left-right (LR) displacement (a) was taken to be 2 and 4 cm, respectively12). Calculations were performed for reference patient 65 kilograms weight and 170 cm height.

RESULTS

Figure 3 presents dependence of L5/S1 spinal load on general displacement of the center of gravity of the upper body. The larger the anterior displacement (x>0), the higher the spinal load. The larger the posterior displacement (x<0), the lower the spinal load. Spinal load is less sensitive to LR displacement () than to AP displacement. Within the studied range of motion, the spinal load is within 10% deviation of load in neutral position (shaded area in Fig.3).

Fig. 3.

Dependence between the L5/S1 spinal load and displacement of the upper body center of gravity in the transversal plane. Shaded area depicts zone of 10% deviation of spinal load from neutral position.

The values of spinal load for the four types of motion (Fig. 2, Table 1) are shown in Fig. 3 and Table 2. LR pendulum motion does not influence L5/S1 load (Fig. 4A). The range of spinal load is approximately the same in all studied motions except LR pendulum motion. The period of the spinal load cycle of AL-PR pendulum motion (Fig. 3B) and ellipse motion (Fig. 3C) is the same as the period of driving kinematics (Table 1) and the loading curve follows simple harmonic function. The ∞ shape displacements used in INFINITY® method exhibit two maxima and two minima of spinal force during one cycle (Fig. 3D), i.e. the frequency of spinal load cycle is twice the frequency of motion cycle.

Table 2.

Maximum and minumum values of L5/S1 spinal load and musculus iliocostalis (ILC) force for various types of upper body center of gravity motion pattern

Motion		(A)	(B)	(C)	(D)
		Linear L-R	Linear AL-PR	Ellipse	∞ shape
Spine load	Max [N]	980.1	1,003.7	1,003.6	1,003.7
	Min [N]	980	956.5	956.4	956.4
ILC load	Max [N]	172.2	178.8	173.9	172.2
	Min [N]	148.9	142.3	147.2	148.9

Fig. 4.

L5/S1 spinal load for various types of upper body center of gravity motion pattern. Inset shows trajectories of center of gravity of upper body in the transversal plane for (A) L-R pendulum motion, (B) AL-PR pendulum motion, (C) ellipse motion, (D) ∞ shape motion.

The changes of muscle force in musculus iliocostalis is shown in Fig. 5. The magnitude of muscle force is almost the same in all studied motions. All studied motions except AL-PR pendulum motion load the left and right segments symmetrically, while AL-PR motion considerably overloads the right muscle segment. The frequency of muscle load follows the harmonic frequency of driving motion except the ∞ motion.

DISCUSSION

The effect of physiotherapy could be measured by means of spinal kinematics21) or electromyography22). However, measurement of muscle force of deep stabilizing muscles and spinal load is technically complex and invasive electromyography techniques are required23). The aim of this study was to estimate the effect of various motion patterns on load of musculoskeletal system of spine using a mathematical model. The benefit of biomechanical modelling is that it can test various motion pattern in controlled environment. The INFINITY method® has been proven as a suitable method in low back pain rehabilitation improving plantar functions, postural control in the upright stance and pain12).

Values of the spinal load proposed within this model are below the threshold proposed by NIOSH (1981) as low back injury risk in the workplace that is 3400 N24). Keeping the force low in rehabilitation of low back pain is required so as not to cause harm to already damaged tissues. In the patients with low back pain, the spinal force could be further elevated by increased muscle coactivation25). We have found, that during micro-motions proposed by the INFINITY® method, the spinal load is not considerably increased. This provides rationale for application of the method in patients with a spine with low load bearing capacity.

A common belief is that muscle activation and muscle tension in physical therapy changes in direct response to specific interventions and enhances spinal stability26). This statement has been confirmed in the presented study: the more complex the motion pattern, the more complex the spine and the muscle loading cycle (Figs. 4, 5).

Fig. 5.

Force in the left and right segments of musculus iliocostalis (ILC) for various types of upper body center of gravity motion pattern. (A) L-R pendulum motion, (B) AL-PR pendulum motion, (C) ellipse motion, (D) ∞ shape motion.

The study further indicates influence of complex motion on passive structures of muscles. ∞ motion pattern (Fig. 2D) activates passive structure by higher frequency loading (Fig. 4D) that is known to be associated with bone remodeling and maintaining bone health27). Higher frequency load could also decrease the stiffness of viscoelastic tendons28) and enhance spine mobility.

The symmetry of the loading cycle is important as well. Acute and chronic low back pain might be associated with ipsilateral atrophy of the multifidus muscles29). Symmetric motion pattern conditions symmetric activation of the muscles (Fig. 5ACD) and hence strengthens the active structures symmetrically. Unsymmetrical motion pattern, e.g. oblique pendulum motion AL-PR, could overload right muscle segments (Fig. 5B). Asymmetric muscle loading predicted by the model in oblique AL-PR pendulum motion is in accordance with EMG measurements for asymmetric load30). An interesting finding is that small LR pendulum motion could effectively exercise deep spine muscles without overloading the spine (Fig. 4A). This conclusion should be further verified in an EMG study.

Studied motion patterns were chosen on the basis of our clinical experience as they required active patient approach to muscle activation31). Clinical study32) indicates that active patient approach is significantly better than chiropractic spinal manipulation and home exercise in terms of patient satisfaction with treatment and trunk muscle endurance. Using an active patient approach offers extensive possibilities of autotherapy of motoric system and simultaneously brain training26). While simple pendulum or ellipse motion pattern results in simple harmonic muscle activation pattern (Fig. 5), the more complex ∞ motion pattern requires a qualitatively more complex muscle activation. Complex muscle activation should be accompanied with active training of the central neural system, that could further increase spinal dynamic stability33).

The mathematical model gives values of spinal load and pattern of EMG activation of back muscles that agree with the previous experimental measurements13). However, a mathematical modeling approach presents several limitations. The model is limited to slow subtle motions that do not considerably change the lumbar musculoskeletal geometry. The spinal load could be increased considerably in case of large displacements due to lateral bending34). The presented results are valid for subtle motions proposed in the INFINITY® method. Including the extensive motion in decimeters would require improvement of the muscle model, mainly by considering muscle wrapping of curved muscle35). The mathematical model further assumes that muscles behave as the ideal force generators and underestimates the effect of nonlinear passive structures. Considering the muscle tendon dynamics and passive properties of the muscle-tendon unit would provide a more accurate force estimation, especially for fast motions36). To model the muscles force dynamics more accurately, a dynamics optimization approach which takes into account dynamic properties of neuromusculoskeletal systems should be used37). The quasi-static model used in this study is in principle not able to fully consider the effects of motion, acceleration or stability goals on the muscle recruitment patterns of the trunk38).

The model was designed to study the INFINITY® method that is primarily based on small micro-motions12). Improvement in the mathematical model would allow its application in fast motion exercise, but would not influence the general conclusions of the study which shows the positive effect of a complex motion pattern. According to Panjabi theory39), there are three components contributing to spinal stability: the passive system constituting bones, joints and ligaments; the active system constituting spine muscles and tendons; and the neural system comprising the nerves and central nervous system which direct and control the active system in providing dynamic stability. The study indicates that more complex motion provides higher frequency load to passive structure and requires complex activation of muscles enhancing neural system involvement. Results support hypothesis that the quality of the movement pattern that is used in rehabilitation is important to achieve spine stability6). General conclusions of the study are in agreement with clinical findings showing positive effect of complex Tai-Chi exercises on low back pain22).

We conclude that the complexity of motion pattern directly influences the spinal load and muscle activation pattern. The simple pendulum-like or circular-like motions induce harmonic muscle activation and spinal load. The figure of eight pattern (∞) proposed in the INFINITY® method doubles the frequency of spinal loading while conditions complex muscle activation pattern that could induce remodeling in passive structures, strengthen active structures, and train the central neural system.