Ukljucite javascript
Register Log in

Archive > Year 2008, Number 1

Computer Simulation And Modeling Of Cartilage Deformation During Athletes' Land Training: One Case Study


The knee is the anatomical spot most often treated by orthopaedic surgeons [2]. Knee problems result in millions of visits to doctors’ offices annually. One of the most common knee injuries in sport is tearing. Intensive research has been carried out over the last 30 years to understand the conditions necessary for maintaining healthy cartilage and the mechanical and biochemical environment that leads to disease [2]. Much has been learned about the morphology, biochemistry and mechanics of cartilage, but many of the important questions remain unanswered. For example, what are the dynamic loading conditions to which cartilage is exposed during complex training as well as players' game activity, and how does this environment influence the health of the tissue? The answers are fundamental for sports training as well as diagnosing and treating joint disease, since dynamic loading affects the movement of tissue growth factors, which must be transported into the cartilage layer to keep it healthy. One key difficulty in understanding these issues has been the inability to measure knee movements accurately in vivo. Conventional motion capture technologies employ video-based systems to track markers attached to the skin. While these systems are non-invasive, the markers affixed to the skin shift relative to the underlying bone [3]. Skin motion effect can be avoided by using medical imaging techniques such as fluoroscopy. Newer technologies, such as cine magnetic resonance imaging (MRI) and single- and biplane fluoroscopy, are capable of accurately recording 3-D joint motion non-invasively. Cine MRI is limited by the restriction it imposes on the movements that the subject can perform, and it cannot be used to investigate weight-bearing activities such as gait [4].

In order to understand knee injuries we need first to have an appreciation of the anatomy and structures of the knee joint. Basically the knee is a hinge joint and classified as a synovial joint, which means it is surrounded by a capsule filled with fluid. This fluid helps nourish and lubricate the structures within [2].

To diagnose a knee problem, doctors often need to observe knee behavior during the gait cycle. Knee motion simulation provides valuable information for physicians to make correct diagnosis as well as for sports trainers to make better training conditions. Specifically, knee motion simulation should describe the kinematic behavior of tibia, femur, patella, and other elements involved in various knee positions. Three dimensional motion tracking systems capture and measure the motion data by attaching markers to various positions of the knee. This information, along with the knowledge of knee kinematics, can help create subject-specific knee motion simulation [3, 4, 5]. However, it is usually difficult to simulate the motion due to the lack of non-invasive methods to measure the knee motion parameters. The major challenge is the definition of femur and tibia rotation axes.



The average test results are expressed in the stages as it is presented in Table 1. The participant is a male soccer player (Age – 18.0 yrs; BH – 1.84 m; BM 71.0 kg; BMI 21.01 kg/m2; VO2max – 68.25 ml/kg/min; Anaerobic threshold – 173 bmp). The starting speed is 8 km/h and each minute the speed is increased by decreasing the interval between the beeps.

Testing Procedure

The Shuttle Run determines the maximum aerobic endurance of a player. The test was developed in 1982 by Leger [6]. There are many names now for the Shuttle Run Test: the Beep Test, the Bleep Test, or the YoYo Test. With this test it is possible to run across the gym between two lines twenty meters apart until maximum endurance of the participant, i.e. until testee quitting. The test is performed indoors and the speed of the participants is determined by the interval between two beeps played on a cassette.

During the testing, the pulse is followed through the pulse-meter (HR Polar). The software note from the pulse-meter is entered into the computer and then put to further use for additional diagnostic purposes.

Indirect calculation of the oxygen consumption, expressed in ml per min per kg, is achieved through a pattern: VO2max= (5.87 x S)-19.458, where S is running time [11]. The values are directly proportional to the aerobic capacity of the testee.

Table 1. Shuttle Run Test protocol for the maximum aerobic endurance

Knee Motion Capture

We used our own motion capture system to collect the specific motion data [7]. The system includes a computer, a high-speed camera and a number of colored markers. For unilateral gait analysis, five markers are attached to the subject at various positions, where three were on the subject's dominant leg (at knee, hip and ankle). Each marker was identified by a different color. The computerized camera system captured the exact motion of these markers, which identified the movement of the body while the subject ran in a straight line. The camera was connected /Basler A601fc camera, Graftek Imaging, Inc, USA/ to a computer that collected gait cycle kinematical data. The 3D tracking software [7] of our own design was then used to identify and convert 2D camera information into 3D motion data at different time frames (Figure 1). The result of motion tracking was a series of 3D coordinates for each numbered marker.

Figure 1. 3D motion capture system

The knee motion involves a series of three rotations (flexion/extension, abduction/adduction, and internal/external rotations) and three translations (anterior/posterior, superior/inferior, and medial/lateral translations). The rotation angles and the translations can be calculated according to the 3D motion data of the stick model recognized from a DVD video movie.

Finite Element Method

We consider large displacements of a porous deformable body filled with fluid, occupying the whole pore volume. The physical quantities at the material point are: the displacement of solid u, relative fluid velocity with respect to the solid (Darcy’s velocity) q, fluid pressure p, swelling pressure pc, and electrical potential f.

The governing equations for the coupled problem are described above. First, consider the solid equilibrium equation,

where ss is stress in the solid phase, n is porosity, k is the permeability matrix, rs is density of solid, b is body force per unit mass, q is relative velocity of fluid, and ü is acceleration of the solid material. The operator LT is

The Equation (1) and the others that follow, correspond to the current configuration tB and are essential for the calculation of the change of material properties during deformation (e.g. porosity). In this section we omit every index to indicate configuration, or time. However, it will be necessary to include them in the later development of an incremental-iterative scheme for the solution of the nonlinear finite element equations.

The equilibrium equation of the fluid phase (no electrokinetic coupling) is

where p is pore fluid pressure, rf is fluid density and n is fluid acceleration. This equation is also known as the generalized Darcy's law. Both equilibrium equations are written per unit volume of the mixture. Combining Equations (1) and (3) we obtain

where s is the total stress which can be expressed in terms of ss and p, as

and s=(1-n)ss+nrf is the mixture density. Here m is a constant vector defined as mT={1 1 1 0 0 0} to indicate that the pressure component contributes to the normal stresses only. We have also taken into account that the pressure has positive sign in compression, while tensional stresses and strains are considered positive. In the following analysis we employ the effective stress, s', defined as

which is relevant for the constitutive relations of the solid. Using the definition of relative velocity q as the volume of the fluid passing in a unit time through a unit area of the mixture (Darcy's velocity),

we transform (3) into

The final continuity equation using the elastic constitutive law and fluid incompressibility is given in the form [1]

The resulting FE system of equations is solved incrementally [1], with time step Δt. We impose the condition that the balance equations are satisfied at the end of each time step (t+Δt). Hence, we have the following system of equations

Results and Discussion

Here our goal is to calculate cartilage deformation at the knee joint based on the subject's knee model and classical biomechanical analysis of the shuttle run test, and then we can visualize the change of cartilage deformation during knee motion. The overall forces on the knee joint are presented in Figure 2.

The geometry of the problem is shown in Figure 3. We analyze the knee joint as electrokinetic transduction in charged, homogenous, isotropic, hydrated material. In the analysis we first use the material constants by Frank and Grodzinsky [8, 9]. After that we fit values of the material constants to experimental data by employing our numerical solutions.

Following the approach of Frank and Grodzinsky [9], the fluid and the solid are assumed incompressible, and the tissue compression is due to the prescribed force of constant amplitude. The initial result for menisci deformation is presented in Figure 4 [10].

Figure 4. Menisci deformation for a peak force. The units presented in the pallete color for deformation are “cm”[10].

Practical Application

Our main goal is to integrate the various steps described above to produce a practical and economical diagnostic tool for non-invasively assessing musculoskeletal function on a subject-specific basis. In summary, virtual biomechanical simulation has phenomenal potential in improving medical science and health care as well as sports trainings. Due to the complexity of medical science and the computer technologies involved in building such simulation systems, there are still a lot of technical challenges in front of us. This work provides solutions to a part of this big problem. Indeed, the 3D reconstruction, motion simulation, and biomechanical visualization techniques presented can be seen as some basic elements of the next generation virtual biomechanical simulation system, a system where all the components of visualization, automated model generation, surgery simulation and surgery assistance come together.


Support of Serbian Ministry of Science and Environmental Protection, Projects ТR – 6209 and OI– 144028, is gratefully acknowledged.