Although alterations in knee joint loading resulting from injury have been shown to influence the development of osteoarthritis, actual in vivo loading conditions of the joint remain unknown. A method for determining in vivo ligament loads by reproducing joint specific in vivo kinematics using a robotic testing apparatus is described. The in vivo kinematics of the ovine stifle joint during walking were measured with 3D optical motion analysis using markers rigidly affixed to the tibia and femur. An additional independent single degree of freedom measuring device was also used to record a measure of motion. Following sacrifice, the joint was mounted in a robotic/universal force sensor test apparatus and referenced using a coordinate measuring machine. A parallel robot configuration was chosen over the conventional serial manipulator because of its greater accuracy and stiffness. Median normal gait kinematics were applied to the joint and the resulting accuracy compared. The mean error in reproduction as determined by the motion analysis system varied between 0.06 mm and 0.67 mm and 0.07 deg and 0.74 deg for the two individual tests. The mean error measured by the independent device was found to be 0.07 mm and 0.83 mm for the two experiments, respectively. This study demonstrates the ability of this system to reproduce in vivo kinematics of the ovine stifle joint in vitro. The importance of system stiffness is discussed to ensure accurate reproduction of joint motion.

Disorders of the first ray of the foot (defined as the hard and soft tissues of the first metatarsal, the sesamoids, and the phalanges of the great toe) are common, and therapeutic interventions to address these problems range from alterations in footwear to orthopedic surgery. Experimental verification of these procedures is often lacking, and thus, a computational modeling approach could provide a means to explore different interventional strategies. A three-dimensional finite element model of the first ray was developed for this purpose. A hexahedral mesh was constructed from magnetic resonance images of the right foot of a male subject. The soft tissue was assumed to be incompressible and hyperelastic, and the bones were modeled as rigid. Contact with friction between the foot and the floor or footwear was defined, and forces were applied to the base of the first metatarsal. Vertical force was extracted from experimental data, and a posterior force of 0.18 times the vertical force was assumed to represent loading at peak forefoot force in the late-stance phase of walking. The orientation of the model and joint configuration at that instant were obtained by minimizing the difference between model predicted and experimentally measured barefoot plantar pressures. The model were then oriented in a series of postures representative of push-off, and forces and joint moments were decreased to zero simultaneously. The pressure distribution underneath the first ray was obtained for each posture to illustrate changes under three case studies representing hallux limitus, surgical arthrodesis of the first ray, and a footwear intervention. Hallux limitus simulations showed that restriction of metatarsophalangeal joint dorsiflexion was directly related to increase and early occurrence of hallux pressures with severe immobility increasing the hallux pressures by as much as 223%. Modeling arthrodesis illustrated elevated hallux pressures when compared to barefoot and was dependent on fixation angles. One degree change in dorsiflexion and valgus fixation angles introduced approximate changes in peak hallux pressure by 95 and 22 kPa, respectively. Footwear simulations using flat insoles showed that using the given set of materials, reductions of at least 18% and 43% under metatarsal head and hallux, respectively, were possible.

Electromyographic (EMG) activity is associated with several tasks prior to landing in walking and running including positioning the leg, developing joint stiffness and possibly control of soft tissue compartment vibrations. The concept of muscle tuning suggests one reason for changes in muscle activity pattern in response to small changes in impact conditions, if the frequency content of the impact is close to the natural frequency of the soft tissue compartments, is to minimize the magnitude of soft tissue compartment vibrations. The mechanical properties of the soft tissue compartments depend in part on muscle activations and thus it was hypothesized that changes in the muscle activation pattern associated with different impact conditions would result in a change in the acceleration transmissibility to the soft tissue compartments. A pendulum apparatus was used to systematically administer impacts to the heel of shod male participants. Wall reaction forces, EMG of selected leg muscles, soft tissue compartment and shoe heel cup accelerations were quantified for two different impact conditions. The transmissibility of the impact acceleration to the soft tissue compartments was determined for each subject/soft tissue compartment/shoe combination. For this controlled impact situation it was shown that changes in the damping properties of the soft tissue compartments were related to changes in the EMG intensity and/or mean frequency of related muscles in response to a change in the impact interface conditions. These results provide support for the muscle tuning idea—that one reason for the changes in muscle activity in response to small changes in the impact conditions may be to minimize vibrations of the soft tissue compartments that are initiated at heel-strike.

A 15 degrees of freedom lumped parameter vibratory model of human body is developed, for vertical mode vibrations, using anthropometric data of the 50th percentile US male. The mass and stiffness of various segments are determined from the elastic modulii of bones and tissues and from the anthropometric data available, assuming the shape of all the segments is ellipsoidal. The damping ratio of each segment is estimated on the basis of the physical structure of the body in a particular posture. Damping constants of various segments are calculated from these damping ratios. The human body is modeled as a linear spring-mass-damper system. The optimal values of the damping ratios of the body segments are estimated, for the 15 degrees of freedom model of the 50th percentile US male, by comparing the response of the model with the experimental response. Formulating a similar vibratory model of the 50th percentile Indian male and comparing the frequency response of the model with the experimental response of the same group of subjects validate the modeling procedure. A range of damping ratios has been considered to develop a vibratory model, which can predict the vertical harmonic response of the human body.

Currently there is no commonly accepted way to define, much less quantify, locomotor stability. In engineering, “orbital stability” is defined using Floquet multipliers that quantify how purely periodic systems respond to perturbations discretely from one cycle to the next . For aperiodic systems, “local stability” is defined by local divergence exponents that quantify how the system responds to very small perturbations continuously in real time . Triaxial trunk accelerations and lower extremity sagittal plane joint angles were recorded from ten young healthy subjects as they walked for 10 min over level ground and on a motorized treadmill at the same speed. Maximum Floquet multipliers (Max FM) were computed at each percent of the gait cycle (from 0% to 100%) for each time series to quantify the orbital stability of these movements. Analyses of variance comparing Max FM values between walking conditions and correlations between Max FM values and previously published local divergence exponent results were computed. All subjects exhibited orbitally stable walking kinematics (i.e., magnitudes of Max FM < 1.0 ), even though these same kinematics were previously found to be locally unstable. Variations in orbital stability across the gait cycle were generally small and exhibited no systematic patterns. Walking on the treadmill led to small, but statistically significant improvements in the orbital stability of mediolateral ( p = 0.040 ) and vertical ( p = 0.038 ) trunk accelerations and ankle joint kinematics ( p = 0.002 ) . However, these improvements were not exhibited by all subjects ( p ⩽ 0.012 for subject × condition interaction effects). Correlations between Max FM values and previously published local divergence exponents were inconsistent and 11 of the 12 comparisons made were not statistically significant ( r 2 ⩽ 19.8 % ; p ⩾ 0.049 ). Thus, the variability inherent in human walking, which manifests itself as local instability, does not substantially adversely affect the orbital stability of walking. The results of this study will allow future efforts to gain a better understanding of where the boundaries lie between locally unstable movements that remain orbitally stable and those that lead to global instability (i.e., falling).

A three-dimensional nonlinear finite element model (FEM) was developed for a parametric study that examined the effect of synthetic augmentation on nonfractured vertebrae. The objective was to isolate those parameters primarily responsible for the effectiveness of the procedure; bone cement volume and bone density were expected to be highly important. Injection of bone cement was simulated in the FEM of a vertebral body that included a cellular model for the trabecular core. The addition of 10% and 20% cement by volume resulted in an increase in failure load, and the larger volume resulted in an increase in stiffness for the vertebral body. Placement of cement within the vertebral body was not as critical a parameter as cement amount. Simulated models of very poor bone quality saw the best therapeutic benefits.

We explored how hip joint actuation can be used to control locomotive bifurcations and chaos in a passive dynamic walking model that negotiated a slightly sloped surface ( γ < 0.019 rad ) . With no hip actuation, our passive dynamic walking model was capable of producing a chaotic locomotive pattern when the ramp angle was 0.01839 rad < γ < 0.0190 rad . Systematic alterations in hip actuation resulted in rapid transition to any locomotive pattern available in the chaotic attractor and induced stability at ramp angles that were previously considered unstable. Our results detail how chaos can be used as a control scheme for locomotion.

The knee joint is partially stabilized by the interaction of multiple ligament structures. This study tested the interdependent functions of the anterior cruciate ligament (ACL) and the medial collateral ligament (MCL) by evaluating the effects of ACL deficiency on local MCL strain while simultaneously measuring joint kinematics under specific loading scenarios. A structural testing machine applied anterior translation and valgus rotation (limits 100 N and 10 N m , respectively) to the tibia of ten human cadaveric knees with the ACL intact or severed. A three-dimensional motion analysis system measured joint kinematics and MCL tissue strain in 18 regions of the superficial MCL. ACL deficiency significantly increased MCL strains by 1.8% ( p < 0.05 ) during anterior translation, bringing ligament fibers to strain levels characteristic of microtrauma. In contrast, ACL transection had no effect on MCL strains during valgus rotation (increase of only 0.1%). Therefore, isolated valgus rotation in the ACL-deficient knee was nondetrimental to the MCL. The ACL was also found to promote internal tibial rotation during anterior translation, which in turn decreased strains near the femoral insertion of the MCL. These data advance the basic structure-function understanding of the MCL, and may benefit the treatment of ACL injuries by improving the knowledge of ACL function and clarifying motions that are potentially harmful to secondary stabilizers.

Total replacement of the glenohumeral joint provides an effective means for treating a variety of pathologies of the shoulder. However, several studies indicate that the procedure has not yet been entirely optimized. Loosening of the glenoid component remains the most likely cause of implant failure, and generally this is believed to stem from either mechanical failure of the fixation in response to high tensile stresses, or through osteolysis of the surrounding bone stock in response to particulate wear debris. Many computational studies have considered the potential for the former, although only few have attempted to tackle the latter. Using finite-element analysis an investigation, taking into account contact pressures as well as glenohumeral kinematics, has thus been conducted, to assess the potential for polyethylene wear within the artificial shoulder. The relationships between three different aspects of glenohumeral design and the potential for wear have been considered, these being conformity, polyethylene thickness, and fixation type. The results of the current study indicate that the use of conforming designs are likely to produce slightly elevated amounts of wear debris particles when compared with less conforming joints, but that the latter would be more likely to cause material failure of the polyethylene. The volume of wear debris predicted was highly influenced by the rate of loading, however qualitatively it was found that wear predictions were not influenced by the use of different polyethylene thicknesses nor fixation type while the depth of wearing was. With the thinnest polyethylene designs ( 2 mm ) the maximum depth of the wear scar was seen to be upwards of 20% higher with a metal-backed fixation as opposed to a cemented design. In all-polyethylene designs peak polymethyl methacrylate tensile stresses were seen to reduce with increasing polyethylene thickness. Irrespective of the rate of loading of the shoulder joint, the current study indicates that it is possible to optimize glenoid component design against abrasive wear through the use of high conformity designs, possessing a polyethylene thickness of at least 6 mm .

Because fall experiments with volunteers can be both challenging and risky, especially with older volunteers, we wished to develop computer simulations of falls to provide a theoretical framework for understanding and extending experimental results. To perform a preliminary validation of the articulated total body (ATB) model for passive falls, we compared the model predictions of fall direction, impact location, and impact velocity as a function of disturbance type (faint, slip, step down, trip) and gait speed (fast, normal, slow) to experimental results with young adult volunteers. The three-dimensional ATB model had 17 segments and 16 joints. Its physical characteristics, environment definitions, contact functions, and initial conditions were representative of our experiment. For each combination of disturbance and gait speed, the ATB model was left to fall passively under gravity once disturbed, i.e., no joint torques were applied, until impact with the floor occurred. Finally, we also determined the sensitivity of the model predictions to changes in the model’s parameters. Our model predictions of fall angles and impact angles were qualitatively in agreement with those observed experimentally for ten and seven of the 12 original simulations, respectively. Quantitatively, the model predictions of fall angles, impact angles, and impact velocities were within one experimental standard deviation for seven, three, and nine of the 12 original simulations, respectively, and within two experimental standard deviations for ten, nine, and 11 of the 12 original simulations, respectively. Finally, the fall angle and impact angle region did not change for 92% and 95% of the 74 input variation simulations, respectively, and the impact velocities were within the experimental standard deviations for 78% of the 74 input variation simulations. Based on our simulations and a sensitivity analysis, we conclude that our preliminary validation of the ATB model for passive falls was successful. In fact, these ATB model simulations represent a significant step forward in fall simulations. We believe that with additional work, the ATB model could be used to accurately simulate a variety of human falls and may be useful in further understanding the etiology and mechanisms of fall injuries such as hip fractures.

Two areas not well researched in the field of seating mechanics are the distribution of normal and shear forces, and how those forces change with seat position. The availability of these data would be beneficial for the design and development of office, automotive and medical seats. To increase our knowledge in the area of seating mechanics, this study sought to measure the normal and shear loads applied to segmental supports in 12 seated positions, utilizing three inclination angles and four levels of seat back articulation that were associated with automotive driving positions. Force data from six regions, including the thorax, sacral region, buttocks, thighs, feet, and hand support were gathered using multi-axis load cells. The sample contained 23 midsized subjects with an average weight of 76.7 kg and a standard deviation of 4.2 kg , and an average height of 1745 mm with a standard deviation of 19 mm . Results were examined in terms of seat back inclination and in terms of torso articulation for relationships between seat positions and support forces. Using a repeated measures analysis, significant differences ( p < 0.05 ) were identified for normal forces relative to all inclination angles except for forces occurring at the hand support. Other significant differences were observed between normal forces behind the buttocks, pelvis, and feet for torso articulations. Significant differences in the shear forces occurred under the buttocks and posterior pelvis during changes in seat back inclination. Significant differences in shear forces were also identified for torso articulations. These data suggest that as seat back inclination or torso articulation change, significant shifts in force distribution occur.

Increasingly complex models of the neck neuromusculature need detailed muscle and kinematic data for proper validation. The goal of this study was to measure the electromyographic activity of superficial and deep neck muscles during tasks involving isometric, voluntary, and reflexively evoked contractions of the neck muscles. Three male subjects ( 28 - 41 years ) had electromyographic (EMG) fine wires inserted into the left sternocleidomastoid, levator scapulae, trapezius, splenius capitis, semispinalis capitis, semispinalis cervicis, and multifidus muscles. Surface electrodes were placed over the left sternohyoid muscle. Subjects then performed: (i) maximal voluntary contractions (MVCs) in the eight directions ( 45 deg intervals) from the neutral posture; (ii) 50 N isometric contractions with a slow sweep of the force direction through 720 deg ; (iii) voluntary oscillatory head movements in flexion and extension; and (iv) initially relaxed reflex muscle activations to a forward acceleration while seated on a sled. Isometric contractions were performed against an overhead load cell and movement dynamics were measured using six-axis accelerometry on the head and torso. In all three subjects, the two anterior neck muscles had similar preferred activation directions and acted synergistically in both dynamic tasks. With the exception of splenius capitis, the posterior and posterolateral neck muscles also showed consistent activation directions and acted synergistically during the voluntary motions, but not during the sled perturbations. These findings suggest that the common numerical-modeling assumption that all anterior muscles act synergistically as flexors is reasonable, but that the related assumption that all posterior muscles act synergistically as extensors is not. Despite the small number of subjects, the data presented here can be used to inform and validate a neck model at three levels of increasing neuromuscular–kinematic complexity: muscles generating forces with no movement, muscles generating forces and causing movement, and muscles generating forces in response to induced movement. These increasingly complex data sets will allow researchers to incrementally tune their neck models’ muscle geometry, physiology, and feedforward/feedback neuromechanics.

Joint injuries during sporting activities might be reduced by understanding the extent of the dynamic motion of joints prone to injury during maneuvers performed in the field. Because instrumented spatial linkages (ISLs) have been widely used to measure joint motion, it would be useful to extend the functionality of an ISL to measure joint motion in a dynamic environment. The objectives of the work reported by this paper were to (i) design and construct an ISL that will measure dynamic joint motion in a field environment, (ii) calibrate the ISL and quantify its static measurement error, (iii) quantify dynamic measurement error due to external acceleration, and (iv) measure ankle joint complex rotation during snowboarding maneuvers performed on a snow slope. An “elbow-type” ISL was designed to measure ankle joint complex rotation throughout its range ( ± 30 deg for flexion/extension, ± 15 deg for internal/external rotation, and ± 15 deg for inversion/eversion). The ISL was calibrated with a custom six degree-of-freedom calibration device generally useful for calibrating ISLs, and static measurement errors of the ISL also were evaluated. Root-mean-squared errors (RMSEs) were 0.59 deg for orientation (1.7% full scale) and 1.00 mm for position (1.7% full scale). A custom dynamic fixture allowed external accelerations ( 5 g , 0 - 50 Hz ) to be applied to the ISL in each of three linear directions. Maximum measurement deviations due to external acceleration were 0.05 deg in orientation and 0.10 mm in position, which were negligible in comparison to the static errors. The full functionality of the ISL for measuring joint motion in a field environment was demonstrated by measuring rotations of the ankle joint complex during snowboarding maneuvers performed on a snow slope.

Skin and garment constitute a dynamic contact system for human body comfort and protection. Although dermatological injuries due to fabric actions during human body movement are common, there is still no general guidance or standard for measuring or evaluating skin/garment contact interactions, especially, during intense sports. A three-dimensional explicit finite element (EFE) model combined with Augmented Lagrange algorithm (ALA) is developed to simulate interactions between skin and fabric during rotation of the arm. Normalized effective shear stresses at the interface between skin and the sleeve during the arm rotation are provided to reflect the severity of the interactions. The effects due to changes in fabric properties, fabric-skin gap, and arm rotation rate are also illustrated. It has been demonstrated from our predictions that factors such as elastic modulus, friction coefficients, density of fabric, and the initial gap between skin and fabric influence significantly the shear stress and thus the discomfort and even injury potential to skin during intensive body movement such as sports and military. Thus this study for the first time confirms quantitatively that poorly chosen fabric with inappropriate garment design renders adverse actions on human skin.

Background: The management of soft tissue balance during surgery is essential for the success of total knee arthroplasty (TKA) but remains difficult, leaving it much to the surgeon’s feel. Previous assessments for soft tissue balance have been performed under unphysiological joint conditions, with patellar eversion and without the prosthesis only at extension and 90 deg of flexion. We therefore developed a new tensor for TKA procedures, enabling soft tissue balance assessment throughout the range of motion while reproducing postoperative joint alignment with the patellofemoral (PF) joint reduced and the tibiofemoral joint aligned. Our purpose in the present study was to clarify joint gap kinematics using the tensor with the CT-free computer assisted navigation system. Method of Approach: Joint gap kinematics, defined as joint gap change during knee motion, was evaluated during 30 consecutive, primary posterior-stabilized (PS) TKA with the navigation system in 30 osteoarthritic patients. Measurements were performed using a newly developed tensor, which enabled the measurement of the joint gap throughout the range of motion, including the joint conditions relevant after TKA with PF joint reduced and trial femoral component in place. Joint gap was assessed by the tensor at full extension, 5 deg, 10 deg, 15 deg, 30 deg, 45 deg, 60 deg, 90 deg, and 135 deg of flexion with the patella both everted and reduced. The navigation system was used to obtain the accuracy of implantations and to measure an accurate flexion angle of the knee during the intraoperative joint gap measurement. Results: Results showed that the joint gap varied depending on the knee flexion angle. Joint gap showed an accelerated decrease during full knee extension. With the PF joint everted, the joint gap increased throughout knee flexion. In contrast, the joint gap with the PF joint reduced increased with knee flexion but decreased after 60 deg of flexion. Conclusions: We clarified the characteristics of joint gap kinematics in PS TKA under physiological and reproducible joint conditions. Our findings can provide useful information for prosthetic design and selection and allow evaluation of surgical technique throughout the range of knee motion that may lead to consistent clinical outcomes after TKA.

Shoulder motion is complex and significant research efforts have focused on measuring glenohumeral joint motion. Unfortunately, conventional motion measurement techniques are unable to measure glenohumeral joint kinematics during dynamic shoulder motion to clinically significant levels of accuracy. The purpose of this study was to validate the accuracy of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics. We have developed a model-based tracking technique for accurately measuring in vivo joint motion from biplane radiographic images that tracks the position of bones based on their three-dimensional shape and texture. To validate this technique, we implanted tantalum beads into the humerus and scapula of both shoulders from three cadaver specimens and then recorded biplane radiographic images of the shoulder while manually moving each specimen’s arm. The position of the humerus and scapula were measured using the model-based tracking system and with a previously validated dynamic radiostereometric analysis (RSA) technique. Accuracy was reported in terms of measurement bias, measurement precision, and overall dynamic accuracy by comparing the model-based tracking results to the dynamic RSA results. The model-based tracking technique produced results that were in excellent agreement with the RSA technique. Measurement bias ranged from − 0.126 to 0.199 mm for the scapula and ranged from − 0.022 to 0.079 mm for the humerus. Dynamic measurement precision was better than 0.130 mm for the scapula and 0.095 mm for the humerus. Overall dynamic accuracy indicated that rms errors in any one direction were less than 0.385 mm for the scapula and less than 0.374 mm for the humerus. These errors correspond to rotational inaccuracies of approximately 0.25 deg for the scapula and 0.47 deg for the humerus. This new model-based tracking approach represents a non-invasive technique for accurately measuring dynamic glenohumeral joint motion under in vivo conditions. The model-based technique achieves accuracy levels that far surpass all previously reported non-invasive techniques for measuring in vivo glenohumeral joint motion. This technique is supported by a rigorous validation study that provides a realistic simulation of in vivo conditions and we fully expect to achieve these levels of accuracy with in vivo human testing. Future research will use this technique to analyze shoulder motion under a variety of testing conditions and to investigate the effects of conservative and surgical treatment of rotator cuff tears on dynamic joint stability.

The locations of the joint axes of the ankle complex vary considerably between subjects, yet no noninvasive method with demonstrated accuracy exists for locating these axes. The moments of muscle and ground reaction forces about the joint axes are dependent on axis locations, making knowledge of these locations critical to accurate musculoskeletal modeling of the foot and ankle. The accuracy of a computational optimization method that fits a two-revolute model to measured motion was assessed using computer-generated data, a two-revolute mechanical linkage, and three lower-leg cadaver specimens. Motions were applied to cadaver specimens under axial load while bone-mounted markers attached to the tibia, talus, and calcaneus were tracked using a video-based motion analysis system. Estimates of the talocrural and subtalar axis locations were computed from motions of the calcaneus relative to the tibia using the optimization method. These axes were compared to mean helical axes computed directly from tibia, talus, and calcaneus motions. The optimization method performed well when the motions were computer-generated or measured in the mechanical linkage, with angular differences between optimization and mean helical axes ranging from 1 deg to 5 deg . In the cadaver specimens, however, these differences exceeded 20 deg . Optimization methods that locate the anatomical joint axes of the ankle complex by fitting two revolute joints to measured tibia-calcaneus motions may be limited because of problems arising from non-revolute behavior.

This study presents an optimized matching algorithm for a dual-orthogonal fluoroscopic image system used to determine six degrees-of-freedom total knee arthroplasty (TKA) kinematics in-vivo. The algorithm was evaluated using controlled conditions and standard geometries. Results of the validation demonstrate the algorithm’s robustness and capability of realizing a pose from a variety of initial poses. Under idealized conditions, poses of a TKA system were recreated to within 0.02±0.01 mm and 0.02±0.03 deg for the femoral component and 0.07±0.09 mm and 0.16±0.18 deg for the tibial component. By employing a standardized geometry with spheres, the translational accuracy and repeatability under actual conditions was found to be 0.01±0.06 mm. Application of the optimized matching algorithm to a TKA patient showed that the pose of in-vivo TKA components can be repeatedly located, with standard deviations less than ±0.12 mm and ±0.12 deg for the femoral component and ±0.29 mm and ±0.25 deg for the tibial component. This methodology presents a useful tool that can be readily applied to the investigation of in-vivo motion of TKA kinematics.

The determination of biomechanical force systems of implanted femurs to obtain adequate strain measurements has been neglected in many published studies. Due to geometric alterations induced by surgery and those inherent to the design of the prosthesis, the loading system changes because the lever arms are modified. This paper discusses the determination of adequate loading of the implanted femur based on the intact femur-loading configuration. Four reconstructions with Lubinus SPII, Charnley Roundback, Müller Straight and Stanmore prostheses were used in the study. Pseudophysiologic and nonphysiologic implanted system forces were generated and assessed with finite element analysis. Using an equilibrium system of forces composed by the Fx (medially direction) component of the hip contact force and the bending moments Mx (median plane) and My (coronal plane) allowed adequate, pseudo-physiological loading of the implanted femur. We suggest that at least the bending moment at the coronal plane must be restored in the implanted femur-loading configuration.

This study investigated the role of the material properties assumed for articular cartilage, meniscus and meniscal attachments on the fit of a finite element model (FEM) to experimental data for meniscal motion and deformation due to an anterior tibial loading of 45 N in the anterior cruciate ligament-deficient knee. Taguchi style L18 orthogonal arrays were used to identify the most significant factors for further examination. A central composite design was then employed to develop a mathematical model for predicting the fit of the FEM to the experimental data as a function of the material properties and to identify the material property selections that optimize the fit. The cartilage was modeled as isotropic elastic material, the meniscus was modeled as transversely isotropic elastic material, and meniscal horn and the peripheral attachments were modeled as noncompressive and nonlinear in tension spring elements. The ability of the FEM to reproduce the experimentally measured meniscal motion and deformation was most strongly dependent on the initial strain of the meniscal horn attachments ( ε 1 H ) , the linear modulus of the meniscal peripheral attachments ( E P ) and the ratio of meniscal moduli in the circumferential and transverse directions ( E θ ∕ E R ) . Our study also successfully identified values for these critical material properties ( ε 1 H = − 5 % , E P = 5.6 MPa , E θ ∕ E R = 20 ) to minimize the error in the FEM analysis of experimental results. This study illustrates the most important material properties for future experimental studies, and suggests that modeling work of meniscus, while retaining transverse isotropy, should also focus on the potential influence of nonlinear properties and inhomogeneity.