View metadata, citation and similar papers at core.ac.uk brought to you by CORE

provided by Cleveland-Marshall College of Law

Cleveland State University EngagedScholarship@CSU

Mechanical Engineering Faculty Publications Mechanical Engineering Department

10-2000 The nflueI nce of Orthotic Devices and Vastus Medialis Strength and Timing on Patellofemoral Loads During Running R. R. Neptune University of Calgary, [email protected]

I. C. Wright University of Calgary

Antonie J. van den Bogert Cleveland State University, [email protected] Follow this and additional works at: https://engagedscholarship.csuohio.edu/enme_facpub Part of the Biomechanical Engineering Commons How does access to this work benefit oy u? Let us know! Publisher's Statement NOTICE: this is the author’s version of a work that was accepted for publication in Clinical Biomechanics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Clinical Biomechanics, 15, 8, (10-01-2000); 10.1016/S0268-0033(00)00028-0

Original Citation Neptune, R. R., Wright, I. C., and van den Bogert, A. J., 2000, "The nflueI nce of Orthotic Devices and Vastus Medialis Strength and Timing on Patellofemoral Loads during Running," Clinical Biomechanics, 15(8) pp. 611-618.

This Article is brought to you for free and open access by the Mechanical Engineering Department at EngagedScholarship@CSU. It has been accepted for inclusion in Mechanical Engineering Faculty Publications by an authorized administrator of EngagedScholarship@CSU. For more information, please contact [email protected]. The influence of orthotic devices and vastus medialis strength and timing on patellofemoral loads during running

R.R. Neptune "', I.e. Wright " A.l. van den Bogert b • Hllmon Performance u'/wra/ory, The Unh'ersily of Calgory. Calgary. A 8, Conar/a T2N I N4 b Departllll!/!/ of Biomedical Engineering. The CIl't'ellmd Clinic FOll/lt/alion, CIeL'eflllld, 01144195, USA

1. Introduction (32%) and older (41%) age populations [I ). Depending on the study citcd, the overall yearly incidence ra le for Running is a popular form of recreation and fitness running injuries is between 37% and 56% [2] with 42% of and is cen tral to many sports. However, running is also these inju ries related to the and 26% of knee in­ associated with a high injury rate for both younger juries associated with patellofemoral pain (PFP) [3]. PFP includes all d isorders associated wi th d iscomfort on the anterior side of the knee joint. Two of the more • Corresponding author. Present address: Rehabilitation R&D common, conservative PFP treatments are exercise Center (\53). VA Palo Alto He3lth Care System. 380\ Miranda Avenue. Palo Alto. CA 94304. USA. programs and the prescription of orthoses [4]. Although £-m(li/ (lddress: [email protected] (R.R. Neptune). conservative treatments have been reported to be 612 successful, the efectiveness is often based on subjective between excessive pronation and PFP and found a or empirical data [5]. Therefore, there is considerable signifcant association. The mechanism behind the pain disagreement within the scientifc community on the is speculated to be a kinematic coupling between ex­ mechanisms causing PFP and consequently, it is not cessive foot pronation and prolonged tibial rotation [18- clear which treatment is the most efective for specifc 20]. Tiberio [20] used a theoretical model to examine the patients [5-7]. Further, running mechanics and muscu­ relationship between excessive pronation and lateral loskeletal properties vary among individuals and may patellofemoral joint compression. Tiberio [20] suggested play an important role in the efectiveness of individual that during gait, the prolonged tibial rotation prevents treatments. the knee from extending freely. To compensate for this, Studies have suggested that diferences in relative the rotates internally relative to the tibia causing muscle forces exerted on the between the vastus changes to the patella contact force and pressure pat­ medialis (VM) and vastus lateralis (VL) contribute to terns. This biomechanical explanation was termed lateral patellar tracking and malalignment, which have compensatory internal rotation of the femur, although been linked to PFP [4-7]. The lateral tracking and other studies have provided alternative explanations of malalignment occurs when the lateral pull of VL is not the relationship between excessive foot pronation and adequately balanced by pull of the vastus medialis lon­ PFP [5]. These studies suggest that further research is gus (VML) and vastus medialis oblique (VMO). This needed to understand the biomechanical relationship imbalance can lead to altered patellofemoral contact between excessive pronation and PFP before the efec­ forces and pressures that may lead to PFP. The VMO is tiveness of foot orthoses verses exercise programs can be deemed the primary medial patellar stabilizer due to the evaluated. anatomical diferences between the VML and VMO [8]. Running is a highly dynamic and complex movement Therefore, much efort has been directed towards iden­ that makes it difcult to predict the efect of these tifying techniques to selectively strengthen the VMO to treatments on musculoskeletal loading during move­ improve the force balance on the patella [4,5,7], al­ ment. Changes in the movement caused by these me­ though the efectiveness and feasibility of selective VMO chanical treatments (e.g., orthoses, [21]) result in strengthening have been widely questioned [6,7,9,10]. changes in the muscle kinematics and therefore the Much of the controversy regarding the efectiveness of muscle forces through the intrinsic muscle force, length VMO strengthening is caused by the difculty of iso­ and velocity relationships. The dynamic interaction lating the treatment efect in patients. Electrical stimu­ within the musculoskeletal system makes such changes lation may be the only way to strengthen the VMO difcult to predict and interpret and are often counter­ selectively [11], but should not be recommended unless intuitive [22]. Thus, the mechanisms behind the changes its efect on patellofemoral joint mechanics is shown to in patellofemoral joint loading are difcult to identify. be signifcant. Therefore, a theoretical analysis is needed Forward dynamic simulations of running have been to quantify the efect of VMO strengthening on patel­ developed that allow for the precise identifcation of lofemoral joint loads to justify further studies seeking to muscle forces and joint loading and have contributed to identify techniques to selectively strengthen the VMO. the understanding of mechanisms involved in loading From a neuromechanical perspective, selective VMO the leg during running (e.g., [23-25]). These models activation can also reduce the lateral dominance of VL were, however, limited to the initial impact phase (0- by either increasing the intensity of VMO relative to VL, 50 ms). Since the etiology and mechanics of PFP and the or by initiating VMO activity prior to VL [6]. In patients efectiveness of conservative treatments are not well with PFP, it is speculated that VL is activated earlier understood, a theoretical model and simulation ap­ than VMO [6,12] and at a higher relative intensity [13]. proach may contribute to our understanding and pro­ But studies examining this dynamic imbalance theory vide the information necessary to design efective have provided conficting results [7,12-16], and the rehabilitation protocols. Therefore, the goal of this functional signifcance of the identifed timing difer­ study was to use a three­dimensional forward dynamic ences, on the order of 5 ms, has been questioned [6,7]. simulation of running to investigate the efectiveness of Methodological diferences between studies and the in­ selective VMO strength and activation and the pre­ herent difculty identifying muscle excitation onset and scription of orthoses on patellofemoral loads. ofset make reconciling diferences difcult, but these studies suggest a theoretical investigation into the in­ fuence of vasti timing and intensity is warranted. 2. Methods Another conservative treatment for PFP has been the prescription of foot orthoses [4]. Eng and Pierrynowski 2.1. Musculoskeletal model [17] found that foot orthoses combined with an exercise program can be an efective method to reduce PFP in A forward dynamic musculoskeletal model was used young females. Bahlsen [18] examined the relationship to produce simulations of heel-toe running. The model 613 was previously described in detail [26] and will be de­ scribed briefy here. The musculoskeletal model was developed using DADS 8.5 software (CADSI, Coral­ ville, IA, USA) and consisted of rigid segments repre­ senting the foot, talus, shank, patella, of the support leg, pelvis, and a rest­of­body segment. The model was dimensioned to represent a male subject with a height of 180 cm and a mass of 75 kg. Musculoskeletal geometry was based on the work of Delp [27] and seg­ ment masses and inertial properties were determined using regression equations [28,29]. The joint models were based on existing literature [27] and further developed to allow three­dimensional movements at the knee. The hip joint was modeled as a spherical joint while rotations about the subtalar and talocrural joints were allowed at the ankle [30]. The tibiofemoral joint was modeled with three degrees of freedom with a moving center of rotation for fexion/ Fig. 1. Musculoskeletal simulation model. The muscles included in the model were the (GMAX), (PSOAS), ad­ extension and appropriate passive stifness for ab­ad­ ductor magnus (ADDMAG), (GMED), duction and internal/external rotation. The patella was (HAMS), rectus femoris (RF), vasti lateralis, intermedius, medialis constrained to move along a prescribed trajectory rela­ longus and medialis obliques (VAS), gastrocnemius (GAS), peroneus tive to the femur [27]. The rotation and anterior/poste­ longus (PER), soleus (SOL), tibialis posterior (TIBPOST), tibialis rior translation were functions of the superior/inferior anterior (TA), extensor digitorum longus (EXTDIG) and fexor digi­ torum longus (FLEXDIG). The four vasti muscles all received the translation relative to the femur. The patella was at­ same excitation signal. tached to the tibia by an inextensible patellar tendon. The foot consisted of three segments: rear­foot, mid­foot and toes. Flexion/extension and internal/external rota­ tions were allowed between the rear-mid and mid­foot, and fexion/extension was allowed between the mid­foot and toes. Passive stifness torques were applied at these joints so that realistic displacements were achieved during the mid­stance. Only translation, controlled by a passive spring and damper, was allowed between the pelvis and the rest­of­body segment. The model had a total of 20 degrees of freedom. Fig. 2. Array of ground contact elements on the rear­foot, mid­foot Fourteen functionally independent muscle groups and toes representing a soft running shoe. The solid lines crossing the foot segment indicate the intersections between the rear­foot, mid­foot were used to drive the model (Fig. 1). The vasti muscle and toes. The point of contact is the center of each rectangular area. group consisted of four individual muscles all receiving the same muscle excitation; the VL, VML, VMO and vastus intermedius with geometry based on Javadpour sue. The anterior-posterior, medial-lateral, and vertical et al. [31]. The force-length-velocity characteristics of force calculations as well as the determination of shoe the muscles were represented by a Hill­based lumped specifc parameters were presented in detail [26]. parameter model (e.g., [32]). The complete set of mus­ culoskeletal model parameters can be found in Neptune 2.3. Subject specifc simulations et al. [26]. Subject specifc simulations were produced by 2.2. Ground contact model searching for the muscle excitation patterns that mini­ mized the diference between simulated and experimen­ The contact between the foot and the ground was tally measured segment orientations and external modeled by 66 discrete independent visco­elastic ele­ ground reaction force profles (Eq. (1)) for each subject ments, each attached to one of the three foot segments in (n = 9) using a simulated annealing optimization algo­ locations that describe the three­dimension exterior rithm [33]. The objective function to be minimized was surface of a shoe when the foot joints are in a neutral defned as position (Fig. 2). Each element permitted deformation ( )2 perpendicular to the foor and represented the mechan­ ( )= (1) 2 ical properties of the shoe sole and underlying soft tis­ =1 =1 614 where is the measurement of variable j at time step i, toe­of was indicated when the vertical ground reaction 2 the simulation data corresponding to , the force fell below 20 N. All trials were normalized to the average inter­trial variability of variable j,and duration from heel strike to toe­of, and the forces and =( 1 ... 42) are the parameters describing the neu­ joint angles were sampled at 100ths of this stance du­ ral excitation patterns. ration. The forces were then normalized to body weight. A similar form of the objective function (Eq. (1)) was Neither the kinematic nor the ground reaction force used in a previous pedaling simulation study [34] and data were fltered, but only averaged across trials. efectively produced simulations that reproduced ex­ perimental data. The specifc quantities evaluated in Eq. 2.5. Simulation experiments (1) were the four Euler parameters describing the global attitude matrix of the rear­foot, shank and thigh, and Experiments were performed with the nine subject the three (x y z) components of the ground reaction specifc simulations to examine the infuence of vastus force, resulting in a total of m = 15 variables included in medialis oblique strength (VMOS), vastus medialis ob­ the objective function (Eq. (1)). The optimization was lique excitation timing (VMOT) and orthoses (ORTH) terminated when the cost function did not decrease the on lateral patella-femoral loads. VMOS was increased cost function by 1% within 500 function calls. This op­ by 10% from the normal strength and VMOT was both timization procedure was used to obtain nine diferent delayed and advanced 5 ms relative to the other vasti simulations of running stance, based on movement data muscles. To examine the infuence of ORTH, a medial from nine diferent subjects. arch orthosis obtained from a sports rehabilitation clinic was modeled in the shoe by increasing the stifness of the 2.4. Movement data individual shoe elements in the medial arch area in the pattern of the orthosis (Fig. 2). Once the stifness con­ Experimental data were collected from a group of tour replicating the orthosis was generated, the stifness nine healthy male subjects (height = 177.0 cm (S.D., values were uniformly scaled until changes in the ankle 10.4); weight = 73.3 kg (S.D., 12.0), age = 22.2 years joint kinematics were similar to those experimentally (S.D., 2.1)) during heel-toe running to provide initial measured [21]. The lateral patella-femoral joint load conditions for the simulation (positions and velocities of was defned as the medio­lateral component of the joint the body segments at heel strike) and tracking data (} ) constraint force vector that keeps the patella in the for the optimization algorithm. The subjects volunteered femoral groove and extracted from the simulation data to participate in the study and informed consent was using the DADS 8.5 post­processor. obtained before the data collection. To assess the efect of these treatments on lateral The subjects performed ten trials of heel-toe running patellofemoral loads, a two­way mixed factor analysis of at 4.0 m/s, while kinematic and ground reaction force variance was performed on the peak and average lateral data were collected. A high speed video system (Motion joint contact forces (P < 0.05). When signifcant treat­ Analysis Corp., Santa Rosa, CA, USA) was used to ment efects were detected, a Bonferroni pairwise com­ record the three­dimensional motion of the body seg­ parison was performed to identify which treatments ments at 240 Hz. Three retro­refective markers were were signifcantly diferent from the normal condition attached to each subject's right shoe (lateral head of the (P < 0.05). ffth metatarsal, posterior heel, superior lateral aspect of the navicular), shank (head of fbula, anterior mid­shaft of tibia, and distal fbula just proximal to the lateral 3. Results maleolus), thigh (greater trochanter, anterior mid­thigh, lateral femoral epicondyle) and pelvis (left and right The optimization was able to fnd the muscle controls anterior superior illiac spines, and right posterior supe­ (excitation onset, ofset and magnitude) to reproduce the rior illiac spine). The marker data were used to recon­ salient features of the experimentally collected data for struct the position and orientation of each body segment the nine subjects. Typical simulation results have been [35]. From the limb segment orientations, hip and knee previously reported [26]. The major joint angles and angles were determined using the joint­coordinate sys­ ground reaction forces were almost always within 2 S.D. tem of Grood and Suntay [36] and the subtalar and of the subject's data, and the muscle stimulation pat­ talocrural joint angles were determined using the joint­ terns compared well with the collected EMG data. coordinate system of Inman [30]. Typical medio­lateral patellofemoral loading patterns Ground reaction force data were collected simulta­ are shown in Fig. 3. In this subject, there was an impact neously with the kinematic data at 2400 Hz using a force peak on the lateral side, followed by a phase of medial platform (Kistler Instumente AG, Winterthur, Switzer­ loading and a second phase of lateral loading during the land). The time of touch­down was determined when the active push­of. The impact loading was on the lateral vertical ground reaction force frst exceeded 20 N and side of the joint for all subjects, except for two subjects 615

Table 1 Group average (N = 9, S.D. in parentheses) of the peak and average lateral constraint force for the normal condition and the four treat­ ments: VMOS, VMOT+, VMOT- and ORTHa

Condition Peak lateral contact Average lateral force (N) contact force (N)

Normal 230 (206) -5 (236) ORTH 221 (198) -14 (235)' VMOS 220 (200)' -17 (237)' VMOT+ 220 (204) -12 (236)' VMOT- 255 (208)' 1 (234)'

a Note: a negative value for average lateral contact force indicates that contact occurred more on the medial than on the lateral side, when Fig. 3. Typical simulated mediolateral loading patterns for the patel­ averaged over the stance phase. lofemoral joint during the stance phase of running. The following * Signifcantly diferent from the normal condition (P < 0.05). treatments were applied: ORTH (orthosis), VMOS (strengthening of the vastus medialis oblique muscle), and VMOT+ (time advance of activation of the vastus medialis oblique muscle). vastus medialis oblique excitation timing advance (VMOT+) produced either a reduction in lateral joint loading or no efect at all. The VMO time delay pro­ duced an increase in lateral joint loading, which was about twice as large as the decrease produced by the VMO time advance. The statistical analysis revealed that VMOS was the only treatment to signifcantly re­ duce the peak lateral constraint force (Table 1). Al­ though ORTH and VMOT+ decreased the constraint force, the efect was not statistically signifcant. Vastus medialis oblique excitation timing delay (VMOT-) sig­ nifcantly increased the peak force across all subjects. For the average lateral constraint force, all treatments had a signifcant efect with ORTH, VMOS and VMOT+ decreasing the lateral force while VMOT- in­ creased the load (Table 1). Table 1 also illustrates the large variation in joint loading patterns between the modeled subjects.

4. Discussion

The goal of this work was to examine the infuence of two commonly prescribed treatments (vastus medialis strengthening and orthoses) for PFP and the functional Fig. 4. Average and peak lateral loads in the patellofemoral joint. signifcance of timing defcits between VMO and VL on Values represent the diference in the joint load relative to the normal lateral patellofemoral joint forces. The forward dynamic simulation. Positive and negative values indicate an increase and de­ simulation model used in the present study had been crease in the constraint force, respectively. shown to closely replicate experimentally collected ki­ netic, kinematic and EMG data during the same run­ who had no impact peak. A lateral active peak occurred ning conditions as the present study [26]. The close in four of the nine subjects, and in three of those this agreement between the subject and model kinematic, peak was higher than the impact peak (in the control ground reaction force and muscle stimulation data in­ condition). The impact peak was never afected by the dicated that the model's performance was indistin­ orthosis, and the active peak was not or only slightly guishable from the corresponding human subject. In afected by VMO activation timing. addition, the model's behavior is based on a combina­ The efect of the four treatments on the peak and tion of well­established knowledge: musculoskeletal average lateral patellofemoral load for the group is anatomy, Newtonian mechanics, and Hill's three­com­ presented in Fig. 4. The treatments ORTH, VMOS, and ponent muscle model. These considerations provided 616 confdence in the model for experiments examining the efectiveness of common treatments for PFPs that can­ not be performed on human subjects. The simulation approach also allowed for the explicit examination of the orthosis efect without the confounding efect of muscle coordination adaptation. A potential limitation of our simulation approach is that the musculoskeletal model was not customized to the anatomy of individual subjects. It is well known that joint anatomy may predispose individuals to overuse injuries. However, having nine models with identical anatomy but diferent running mechanics allowed us to unequivocally attribute diferent responses to treatments Fig. 5. Sensitivity of the peak lateral patellofemoral load to increases in to diferences in running mechanics rather than anato­ VMOS for a typical subject. VMOS was systematically increased from -50% to +50% relative to the normal muscle strength that showed a my. Such insight is clinically important since running linear association between VMOS and the magnitude of the lateral mechanics is more easily altered than anatomy. A po­ joint load. tential limitation of this approach is that all joint con­ tact force results were reported using the same femoral joint coordinate system. Specifcally, the lateral force identifying efective techniques to selectively strengthen was taken as the force perpendicular to the femoral the VMO. sagittal plane, which contained the patellofemoral The relative intensity and timing of VMO activity groove. This plane may have had a diferent orientation relative to VL activity has also been suggested to play an in diferent subjects. Therefore, we cannot have conf­ important role in reducing the lateral loads in the knee dence in the absolute magnitudes of the lateral con­ joint [6]. Although we did not explicitly assess the efect straint force. However, the objective of this investigation of increased VMO intensity, this was indirectly ac­ was to examine the efectiveness of the treatments and counted for in the VMOS experiment. Since muscle identify mechanisms associated with the reduction in force generation is scaled by the magnitude of muscle joint loading. To this end, comparisons were performed activation, increasing the muscle force can be achieved with the normal simulation and the emphasis was placed two ways, either by increasing the maximum isometric on the diferences (Fig. 4). force or increasing the level of activation. Therefore, the The results showed that VMOS had the greatest efect VMOS results suggest that increased neural VMO ex­ in reducing the peak lateral constraint force and all citation may play an important role in reducing the treatments were efective in decreasing the average lat­ lateral joint loading, although the ability of the nervous eral constraint force except VMOT+ (Table 1). Al­ system to selectively activate muscles within a synergistic though not presented in the results, the treatments were group is highly debated [6]. also applied in combination which showed that the joint Previous work has identifed a correlation between loading response was linear (i.e., the reduction or in­ subjects with timing defcits between the VMO, VL and crease in joint loading was additive). Therefore, there PFP [6,12], although the functional signifcance of the did not appear to be any non­linear coupling associated timing defcit (on the order of 5 ms) has been questioned with the combinations of treatments. [6,7]. The results of the present study suggest that a 5 ms Increasing VMO strength was shown to be the most timing advance or delay is indeed important, although efective treatment for reducing both the peak and this theoretical result may depend on the time constants average lateral constraint forces. A post­hoc sensitivity used in the muscle activation model. The timing advance analysis was performed on the degree of strengthening reduced the peak and average force by 9.8% and 6.6%, from -50% to +50% that showed the reduction in the respectively. These reductions in the joint loading were lateral constraint force decreased linearly with in­ similar in magnitude to the orthosis. An intriguing result creased strength (Fig. 5). These results would seem to was the efect of the timing delay. The peak and average support the treatment as a potentially efective method lateral constraint force increased by 25.5% and 5.9%, in reducing lateral patellofemoral loading. While clini­ respectively. Although the timing delay does not have a cal studies have shown that methods including selective great impact on the average force throughout the stance electrical stimulation and biofeedback can produce phase, it does signifcantly afect the peak force that appreciable results [11,37-39], selective VMO corresponds with the impact peak (Fig. 3(b)). The im­ strengthening through isometric and isokinetic exer­ pact peak force is suggested to be a leading cause of cises has been widely questioned [6,7,9,10]. Neverthe­ overuse injuries during repetitive movement tasks [40]. less, the theoretical results of the present study suggest An interesting result was the relatively small and less future research might be usefully directed towards consistent orthosis efect, especially for the peak lateral 617 contact force. Proponents of orthoses have advocated References orthoses as an important component in the treatment of PFP [17]. But the simulations showed that the [1] Matheson GO, Macintyre JG, Taunton JE, Clement DB, orthosis examined in the present study was not as Lloyd­Smith R. Musculoskeletal injuries associated with physi­ cal activity in older adults. Med Sci Sports Exerc 1989;21:379- consistent as VMOS in reducing the joint loading. The 85. lack of a stronger orthosis efect could be related to the [2] van Mechelen W. Running injuries. A review of the epidemio­ specifc orthosis used in the study, which did not logical literature. Sports Med 1992;14:320-35. change the shape of the foot-ground interface, but [3] Clement DB, Taunton JE, Smart GW, McNicol KL. A survey of only increased stifness in the arch area. This type of overuse running injuries. Phys Sports Med 1981;9:47-58. [4] Cutbill JW, Ladly KO, Bray RC, Thorne P, Verhoef M. Anterior orthosis was chosen for this study because it is less knee pain: a review. Clin J Sports Med 1997;7:40-5. likely to require adaptations in muscle activation pat­ [5] Callaghan MJ, Baltzopoulos V. Anterior knee pain: the need for terns. Muscle coordination adaptation to the orthosis objective measurement. Clin Biomech 1992;7:67-74. may be important and produce favorable or adverse [6] Grabiner MD, Koh TJ, Draganich LF. Neuromechanics of the efects. The advantage of the simulation model is that it patellofemoral joint. Med Sci Sports Exerc 1994;26:10-21. [7] Powers CM, Landel R, Perry J. Timing and intensity of vastus allowed us to assess the efectiveness of the orthosis muscle activity during functional activities in subjects with and with all the parameters in the system preserved (e.g., without patellofemoral pain. Phys Ther 1996;76:946-55. muscle coordination). Further, the subjects used in the [8] Lieb FJ, Perry J. Quadriceps function. An anatomical and study were normal, healthy subjects without a history mechanical study using amputated limbs. J Bone Joint Surg Am of PFP and most likely would not be candidates for the 1968;50:1535-48. [9] Garrick JG. Anterior knee pain (chondromalacia patella). Phys orthosis treatment. Specifcally, the lack of orthosis Sports Med 1989;17:75-89. efect in certain subjects may have been due to the fact [10] Mirzabeigi E, Jordan C, Gronley JK, Rockowitz NL, Perry J. that those subjects had no pressure under the medial Isolation of the vastus medialis oblique muscle during exercise. arch during running. It is possible that a simulation Am J Sports Med 1999;27:50-3. control group of subjects with chronic PFP and can­ [11] Werner S, Arvidsson H, Arvidsson I, Eriksson E. Electrical stimulation of vastus medialis and stretching of lateral thigh didates for orthoses would produce greater changes in muscles in patients with patello­femoral symptoms. Knee Surg the joint loading. The specifc orthosis examined in the Sports Traumatol Arthrosc 1993;1:85-92. present study is also one of many types used to treat [12] Voight ML, Wieder DL. Comparative refex response times of PFP. Future studies should be directed towards ex­ vastus medialis obliquus and vastus lateralis in normal subjects amining diferent confgurations and their efectiveness and subjects with extensor mechanism dysfunction. An electr­ omyographic study. Am J Sports Med 1991;19:131-7. towards reducing patellofemoral loads. [13] Mariani PP, Caruso I. An electromyographic investigation of These fndings showed that the orthosis treatment subluxation of the patella. J Bone Joint Surg Br 1979;61­B:169- had a benefcial efect in some subjects and no efect in 71. other subjects. An adverse efect was never observed. [14] Karst GM, Willett GM. Onset timing of electromyographic The efect of selective VMO strengthening was more activity in the vastus medialis oblique and vastus lateralis muscles in subjects with and without patellofemoral pain syndrome. Phys consistent across subjects, but the feasibility of this Ther 1995;75:813-23. treatment is debatable. Activation timing of the VMO [15] MacIntyre DL, Robertson DG. Quadriceps muscle activity in muscle had a consistent efect on lateral patellofemoral women runners with and without patellofemoral pain syndrome. loading, especially the peak force developed during the Arch Phys Med Rehabil 1992;73:10-4. impact phase. The advantage of both the orthosis and [16] Souza DR, Gross MT. Comparison of vastus medialis obliquus: integrated electromyographic ratios be­ strengthening treatments is that these do not require tween healthy subjects and patients with patellofemoral pain. changing the complex neuromuscular coordination Phys Ther 1991;71:310-16. patterns of the patient. Future studies should be di­ [17] Eng JJ, Pierrynowski MR. Evaluation of soft foot orthotics in the rected towards identifying techniques that selectively treatment of patellofemoral pain syndrome. Phys Ther increase VMO strength as a potential treatment for 1993;73:62-8. [18] Bahlsen A. The etiology of running injuries: a longitudinal, PFP, and towards predicting the efcacy of orthotic prospective study. Calgary, AB: University of Calgary, 1988. treatment in patients based on their individual running [19] Nigg BM, Cole GK, Nachbauer W. Efects of arch height of the mechanics. foot on angular motion of the lower extremities in running. J Biomech 1993;26:909-16. [20] Tiberio D. The efect of excessive subtalar joint pronation on patellofemoral mechanics: a theoretical model. J Orthop Sports Acknowledgements Phy Ther 1987;9:160-5. [21] Eng JJ, Pierrynowski MR. The efect of soft foot orthotics on This research was supported by NSERC of Canada, three­dimensional lower­limb kinematics during walking and running. Phys Ther 1994;74:836-44. The Whitaker Foundation and Adidas America. We are [22] Zajac FE, Gordon ME. Determining muscle's force and action in thankful to Trisha Bellchamber and Phuong Ngo for multi­articular movement. Exerc Sport Sci Rev 1989;17:187- collecting the experimental data. 230. 618

[23] Cole GK, Nigg BM, van den Bogert AJ, Gerritsen KG. Lower [32] Soest AJv, Bobbert MF. The contribution of muscle properties in extremity joint loading during impact in running. Clin Biomech the control of explosive movements. Biol Cybern 1993;69: 1996;11:181-93. 195-204. [24] Gerritsen KG, van den Bogert AJ, Nigg BM. Direct dynamics [33] Gofe WL, Ferrier GD, Rogers J. Global optimization of simulation of the impact phase in heel-toe running. J Biomech statistical functions with simulated annealing. J Econometrics 1995;28:661-8. 1994;60:65-99. [25] Wright IC, Neptune RR, van den Bogert AJ, Nigg BM. Passive [34] Neptune RR, Hull ML. Evaluation of performance criteria for regulation of impact forces in heel-toe running. Clin Biomech simulation of submaximal steady­state cycling using a forward 1998;13:521-31. dynamic model. J Biomech Eng 1998;120:334-41. [26] Neptune RR, Wright IC, van den Bogert AJ. A method for [35] Soderkvist I, Wedin PA. Determining the movements of the numerical simulation of single limb ground contact events: skeleton using well­confgured markers. J Biomech 1993;26: application to heel-toe running. Comput Methods Biomech 1473-7. Biomed Eng 1999 (in press). [36] Grood ES, Suntay WJ. A joint coordinate system for the clinical [27] Delp SL, Surgery simulation: a computer graphics system to description of three­dimensional motions: Application to the analyze and design musculoskeletal reconstructions of the lower knee. J Biomech Eng 1983;105:136-44. limb. Ph.D. Dissertation, Department of Mechanical Engineer­ [37] Wise HH, Fiebert IM, Kates JL. EMG biofeedback as treatment ing. CA (Stanford): Stanford University, 1990. for patellofemoral pain syndrome. J Orthop Sports Phys Ther [28] Chandler RF, Clauser CE, McConville JT, Reynolds HM, Young 1984;6:95-103. JW. Investigation of inertial properties of the human body. Wright [38] LeVeau BF, Rogers C. Selective training of the vastus medialis Patterson Air Force Base, Ohio (AMRL­TR­75­137), 1975. muscle using EMG biofeedback. Phys Ther 1980;60:1410-5. [29] Clauser CE, McConville JT, Young JW. Weight, volume and [39] King AC, Ahles TA, Martin JE, White R. EMG biofeedback­ center of mass of segments of the human body. Wright Patterson controlled exercise in chronic arthritic knee pain. Arch Phys Med Air Force Base, Ohio (AMRL­TR­69­70), 1969. Rehabil 1984;65:341-3. [30] Inman VT. The joints of the ankle. Baltimore: Williams & [40] Radin EL, Yang KH, Riegger C, Kish VL, O'Connor JJ. Wilkins, 1976, p. 1-117. Relationship between lower limb dynamics and knee joint pain. [31] Javadpour SM, Finegan PJ, O'Brien M. The anatomy of the J Orthop Res 1991;9:398-405. extensor mechanism and its clinical relevance. Clin J Sport Med 1991;1:229-35.