Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

MRI-based experimentations of fingertip flat compression: Geometrical measurements and finite element inverse simulations to investigate material property parameters

Modeling human-object interactions is a necessary step in the ergonomic assessment of products. Fingertip finite element models can help investigating these interactions, if they are built based on realistic geometrical data and material properties. The aim of this study was to investigate the fingertip geometry and its mechanical response under compression, and to identify the parameters of a hyperelastic material property associated to the fingertip soft tissues.Fingertip compression tests in an MRI device were performed on 5 subjects at either 2 or 4 N and at 15° or 50°. The MRI images allowed to document both the internal and external fingertip dimensions and to build 5 subject-specific finite element models. Simulations reproducing the fingertip compression tests were run to obtain the material property parameters of the soft tissues.Results indicated that two ellipses in the sagittal and longitudinal plane could describe the external fingertip geometry. The internal geometries indicated an averaged maximal thickness of soft tissues of 6.4 ± 0.8 mm and a 4 ± 1 mm height for the phalanx bone. The averaged deflections under loading went from 1.8 ± 0.3 mm at 2 N, 50° to 3.1 ± 0.2 mm at 4 N, 15°. Finally, the following set of parameters for a second order hyperelastic law to model the fingertip soft tissues was proposed: C₀₁ = 0.59 ± 0.09 kPa and C₂₀ = 2.65 ± 0.88 kPa.These data should facilitate further efforts on fingertip finite element modeling.     

An anatomical subject-specific FE-model for hip fracture load prediction

In order to reduce the socio-economic burden induced by osteoporotic hip fractures, finite element models have been evaluated as an additional diagnostic tool for fracture prediction. For a future clinical application, the challenge is to reach the best compromise between model relevance and computing time. Based on this consideration, the current study focused on the development and validation of a subject-specific FE-model using an original parameterised generic model and a specific personalization method. A total of 39 human femurs were tested to failure under a quasi-static compression in stance configuration. The corresponding FE- models were generated and for each specimen the numerical fracture load (FFEM) was compared with the experimental value (FEXP), resulting in a significant correlation (FEXP = 1.006 FFEM with r2 = 0.87 and SEE = 1220 N, p < 0.05) obtained with a reasonable computing time (30 mn). Further in vivo study should confirm the ability of this FE-model to improve the fracture risk prediction.     

Rib fractures under anterior–posterior dynamic loads: Experimental and finite-element study

The purpose of this study was to investigate whether using a finite-element (FE) mesh composed entirely of hexahedral elements to model cortical and trabecular bone (all-hex model) would provide more accurate simulations than those with variable thickness shell elements for cortical bone and hexahedral elements for trabecular bone (hex–shell model) in the modeling human ribs. First, quasi-static non-injurious and dynamic injurious experiments were performed using the second, fourth, and tenth human thoracic ribs to record the structural behavior and fracture tolerance of individual ribs under anterior–posterior bending loads. Then, all-hex and hex–shell FE models for the three ribs were developed using an octree-based and multi-block hex meshing approach, respectively. Material properties of cortical bone were optimized using dynamic experimental data and the hex–shell model of the fourth rib and trabecular bone properties were taken from the literature. Overall, the reaction force–displacement relationship predicted by both all-hex and hex–shell models with nodes in the offset middle-cortical surfaces compared well with those measured experimentally for all the three ribs. With the exception of fracture locations, the predictions from all-hex and offset hex–shell models of the second and fourth ribs agreed better with experimental data than those from the tenth rib models in terms of reaction force at fracture (difference <15.4%), ultimate failure displacement and time (difference <7.3%), and cortical bone strains. The hex–shell models with shell nodes in outer cortical surfaces increased static reaction forces up to 16.6%, compared to offset hex–shell models. These results indicated that both all-hex and hex–shell modeling strategies were applicable for simulating rib responses and bone fractures for the loading conditions considered, but coarse hex–shell models with constant or variable shell thickness were more computationally efficient and therefore preferred.     

Effects of idealized joint geometry on finite element predictions of cartilage contact stresses in the hip

Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error ～0.5 mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

Experimental model of tibial plateau fracture for biomechanical testing

Although adequate reduction and stable fixation have been recognized to be the prime goals in the treatment of displaced tibial plateau fractures, the optimal fixation technique remains controversial. The lack of a reliable model and a standard methodology contribute to this situation. The purpose of this study is to develop an experimental model of a tibial plateau fracture and a testing methodology that reproduces the failure mode commonly seen in the clinical setting. Using solid-foam and composite Sawbones tibiae, three different models of bi-condylar tibial plateau fracture (solid-foam, reinforced solid-foam and composite), six specimens for each model, were created and stabilized with double plating. The specimens were subjected to cyclic axial compression with increasing maximum load until failure. A femoral component of a total knee replacement of similar size and shape to the synthetic tibial surface was used as a load applicator. The experiment was repeated on six specimens of human cadaver tibiae. Among the Sawbones specimens, only the reinforced solid-foam model was found to produce a consistent failure mode (collapse in the medial plateau) comparable to that reported clinically in the literature. This mode of failure was also confirmed by the cadaver experiments. The failure load of the reinforced solid-foam model ranged from 4150 to 4260 N with a mean +/- SD of 4201 +/- 44 N and a coefficient of variance of 0.01, whereas for the cadaver model the failure load ranged from 1675 to 6096 N with a mean +/- SD of 3768 +/- 1482 N and a coefficient of variance of 0.39. We recommend the reinforced-foam model for future mechanical tests to compare different fixation methods for tibial plateau fractures.     

Apparent Young's modulus of vertebral cortico-cancellous bone specimens

Up to now, due to cortical thickness and imaging resolution, it is not possible to derive subject-specific mechanical properties on the 'vertebral shell' from imaging modalities applicable in vivo. As a first step, the goal of this study was to assess the apparent Young's modulus of vertebral cortico-cancellous bone specimens using an inverse method. A total of 22 cortico-cancellous specimens were harvested from 22 vertebral bodies. All specimens were tested in compression until failure. To compute the apparent Young's modulus of the specimen from the inverse method, the boundary conditions of the biomechanical experiments were faithfully reproduced in a finite element model (FEM), and an optimisation routine was used. The results showed a mean of the apparent Young's modulus of 374?±?208?MPa, ranging from 87 to 791?MPa. By computing an apparent Young's modulus of a cortico-cancellous medium, this study gives mechanical data for an FEM of an entire vertebra including an external shell combining both bone tissues.     

Fracture toughness and fatigue crack propagation rate of short fiber reinforced epoxy composites for analogue cortical bone

Third-generation mechanical analogue bone models and synthetic analogue cortical bone materials manufactured by Pacific Research Laboratories, Inc. (PRL) are popular tools for use in mechanical testing of various orthopedic implants and biomaterials. A major issue with these models is that the current third-generation epoxy-short fiberglass based composite used as the cortical bone substitute is prone to crack formation and failure in fatigue or repeated quasistatic loading of the model. The purpose of the present study was to compare the tensile and fracture mechanics properties of the current baseline (established PRL "third-generation" E-glass-fiber-epoxy) composite analogue for cortical bone to a new composite material formulation proposed for use as an enhanced fourth-generation cortical bone analogue material. Standard tensile, plane strain fracture toughness, and fatigue crack propagation rate tests were performed on both the third- and fourth-generation composite material formulations using standard ASTM test techniques. Injection molding techniques were used to create random fiber orientation in all test specimens. Standard dog-bone style tensile specimens were tested to obtain ultimate tensile strength and stiffness. Compact tension fracture toughness specimens were utilized to determine plane strain fracture toughness values. Reduced thickness compact tension specimens were also used to determine fatigue crack propagation rate behavior for the two material groups. Literature values for the same parameters for human cortical bone were compared to results from the third- and fourth-generation cortical analogue bone materials. Tensile properties of the fourth-generation material were closer to that of average human cortical bone than the third-generation material. Fracture toughness was significantly increased by 48% in the fourth-generation composite as compared to the third-generation analogue bone. The threshold stress intensity to propagate the crack was much higher for the fourth-generation material than for the third-generation composite. Even at the higher stress intensity threshold, the fatigue crack propagation rate was significantly decreased in the fourth-generation composite compared to the third-generation composite. These results indicate that the bone analogue models made from the fourth-generation analogue cortical bone material may exhibit better performance in fracture and longer fatigue lives than similar models made of third-generation analogue cortical bone material. Further fatigue testing of the new composite material in clinically relevant use of bone models is still required for verification of these results. Biomechanical test models using the superior fourth-generation cortical analogue material are currently in development.     

Diametral compression of non-circular diaphyseal bone sections

Many research endeavors involve strength testing of long bones, frequently using whole-bone four-point bending models. Recently, diametral compression of short sections has been used to quantify local mechanical parameters and effects of treatment, but testing of biologically derived samples entails a number of added complications, such as the non-circularity of bone sections, ambiguity of load orientation during testing, thickness variation in a section, and size and shape variation between sections in a single sample. In order to quantify the effects of these confounding factors, finite element diametral compression models of a number of bone sections were compared with simplified circular and elliptical sections. Each anatomic section was tested in all rotationally stable load configurations. A high degree of correlation was observed between the anatomic sections and their circular and elliptic analogs, indicating that meaningful comparisons may be made between bone sections of disparate geometry. The aspect ratio and shape of the bone sections did not have a significant impact on the maximum in-plane principal stresses, whereas stresses were strongly dependant on the mean thickness and spatial thickness variation. Some variation due to load orientation was observed. These results indicate that diametral ring compression testing of anatomic sections can be used effectively to measure structural and material parameters of long bones, and that anatomic variation can be successfully accommodated. The ability to use diametral compression testing should allow researchers to obtain many more samples from each specimen than whole-bone bending without the difficulty of extracting solid core or dog-bone samples.     

Subject-specific mechanical properties of vertebral cancellous bone assessed using a low-dose X-ray device

ObjectivesAccording to the inter-individual variability of bone mechanical properties, subject-specific evaluation of the cancellous bone Young's modulus is needed to build finite-element models predicting vertebral strength with accuracy. Relationships based on the density assessed by quantitative computed tomography were proposed. However, quantitative computed tomography is not always suited for the analysis of the whole spine for patients’ follow-up because of the high radiation dose. Hence, this study aims at evaluating the mechanical properties of the vertebral cancellous bone using a low-dose X-ray device.Material and methodsNineteen vertebrae were considered. Biplanar radiographs were made using the low-dose EOS^® system with a dual-energy modality to evaluate antero-posterior and lateral areal bone mineral densities. A cylindrical sample was extracted from each vertebral body and tested until failure to assess the Young's modulus and the ultimate stress of the vertebral cancellous bone.Results and discussionMechanical properties were significantly related to the EOS^® areal densities. On one hand, the relationships remained less predictive than those based on quantitative computed tomography, but on the other hand, they better predict mechanical properties than previous studies using dual X-ray absorptiometry (clinical gold standard system for density assessment).ConclusionThe study shows the feasibility to predict the Young's modulus of the vertebral cancellous bone from the whole vertebral areal bone mineral density (BMD). It gives promising prospects to build finite-element models, including both subject-specific geometry and subject-specific mechanical properties by using a low-dose X-ray device for regions where high radiation doses would limit tomography assessment possibilities.     

Age estimation from cortical index of the human clavicle in northwest Indians

The proportion of the cortical thickness to the total diameter of the bone (cortical index) was calculated in paired clavicles obtained from 128 male and 82 female medicolegal postmortem subjects who were apparently healthy prior to their accidental death. The ages of the subjects varied from 15 to 85 years. The clavicles were cut either horizontally or parasagittally, and measurements were taken at midclavicle. From the ages of 15 to 30 years, the cortical index increased. It decreased steadily thereafter, with an initial sharp decrease in the age group 31-40 years in both sexes. After the age of 40 years, this rapid decrease in the index continued in the females, but became slow and gradual in the males. Differences between left and right sides were statistically insignificant in both sexes (P greater than 0.05). However, the sexual differences were significant (P less than 0.01) in the age groups from 41 years onwards.     

Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis.

In Part I we reported the results of linear finite element models of the proximal femur generated using geometric and constitutive data collected with quantitative computed tomography. These models demonstrated excellent agreement with in vitro studies when used to predict ultimate failure loads. In Part II, we report our extension of those finite element models to include nonlinear behavior of the trabecular and cortical bone. A highly nonlinear material law, originally designed for representing concrete, was used for trabecular bone, while a bilinear material law was used for cortical bone. We found excellent agreement between the model predictions and in vitro fracture data for both the onset of bone yielding and bone fracture. For bone yielding, the model predictions were within 2 percent for a load which simulated one-legged stance and 1 percent for a load which simulated a fall. For bone fracture, the model predictions were within 1 percent and 17 percent, respectively. The models also demonstrated different fracture mechanisms for the two different loading configurations. For one-legged stance, failure within the primary compressive trabeculae at the subcapital region occurred first, leading to load transfer and, ultimately, failure of the surrounding cortical shell. However, for a fall, failure of the cortical and trabecular bone occurred simultaneously within the intertrochanteric region. These results support our previous findings that the strength of the subcapital region is primarily due to trabecular bone whereas the strength of the intertrochanteric region is primarily due to cortical bone.     

The non-linear response of a muscle in transverse compression: assessment of geometry influence using a finite element model

Most recent finite element models that represent muscles are generic or subject-specific models that use complex, constitutive laws. Identification of the parameters of such complex, constitutive laws could be an important limit for subject-specific approaches. The aim of this study was to assess the possibility of modelling muscle behaviour in compression with a parametric model and a simple, constitutive law. A quasi-static compression test was performed on the muscles of dogs. A parametric finite element model was designed using a linear, elastic, constitutive law. A multi-variate analysis was performed to assess the effects of geometry on muscle response. An inverse method was used to define Young's modulus. The non-linear response of the muscles was obtained using a subject-specific geometry and a linear elastic law. Thus, a simple muscle model can be used to have a bio-faithful, biomechanical response.     

Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.     

A new cortical thickness mapping method with application to an in vivo finite element model

Finite element modelling of musculoskeletal systems, with geometrical structures constructed from computed tomography (CT) scans, is a useful and powerful tool for biomechanical studies. The use of CT scans from living human subjects, however, is still limited. Accurate reconstruction of thin cortical bone structures from CT scans of living human subjects is especially problematic, due to low CT resolution that results from mandatory low radiation doses and/or involuntary movements of the subject. In this study, a new method for mapping cortical thickness is described. Using the method, cortical thickness measurements of a coxal (pelvis) bone obtained from CT scans of a cadaver were mapped to the coxal geometry as obtained through CT scans of a live human subject, resulting in accurate cortical thickness while maintaining geometric fidelity of the live subject. The mapping procedure includes shape-preserving parameterisation, mesh movement and interpolation of thickness using a search algorithm. The methodology is applicable to modelling of other bones where accurate cortical thickness is needed and for which such data exist.     

The Mechanical Benefit of Medial Support Screws in Locking Plating of Proximal Humerus Fractures

Background

The purpose of this study was to evaluate the biomechanical advantages of medial support screws (MSSs) in the locking proximal humeral plate for treating proximal humerus fractures.

Methods

Thirty synthetic left humeri were randomly divided into 3 subgroups to establish two-part surgical neck fracture models of proximal humerus. All fractures were fixed with a locking proximal humerus plate. Group A was fixed with medial cortical support and no MSSs; Group B was fixed with 3 MSSs but without medial cortical support; Group C was fixed with neither medial cortical support nor MSSs. Axial compression, torsional stiffness, shear stiffness, and failure tests were performed.

Results

Constructs with medial support from cortical bone showed statistically higher axial and shear stiffness than other subgroups examined (P<0.0001). When the proximal humerus was not supported by medial cortical bone, locking plating with medial support screws exhibited higher axial and torsional stiffness than locking plating without medial support screws (P≤0.0207). Specimens with medial cortical bone failed primarily by fracture of the humeral shaft or humeral head. Specimens without medial cortical bone support failed primarily by significant plate bending at the fracture site followed by humeral head collapse or humeral head fracture.

Conclusions

Anatomic reduction with medial cortical support was the stiffest construct after a simulated two-part fracture. Significant biomechanical benefits of MSSs in locking plating of proximal humerus fractures were identified. The reconstruction of the medial column support for proximal humerus fractures helps to enhance mechanical stability of the humeral head and prevent implant failure.     

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion–extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.     

Muscle parameters estimation based on biplanar radiography

The evaluation of muscle and joint forces in vivo is still a challenge. Musculo-Skeletal (musculo-skeletal) models are used to compute forces based on movement analysis. Most of them are built from a scaled-generic model based on cadaver measurements, which provides a low level of personalization, or from Magnetic Resonance Images, which provide a personalized model in lying position. This study proposed an original two steps method to access a subject-specific musculo-skeletal model in 30 min, which is based solely on biplanar X-Rays. First, the subject-specific 3D geometry of bones and skin envelopes were reconstructed from biplanar X-Rays radiography. Then, 2200 corresponding control points were identified between a reference model and the subject-specific X-Rays model. Finally, the shape of 21 lower limb muscles was estimated using a non-linear transformation between the control points in order to fit the muscle shape of the reference model to the X-Rays model. Twelfth musculo-skeletal models were reconstructed and compared to their reference. The muscle volume was not accurately estimated with a standard deviation (SD) ranging from 10 to 68%. However, this method provided an accurate estimation the muscle line of action with a SD of the length difference lower than 2% and a positioning error lower than 20 mm. The moment arm was also well estimated with SD lower than 15% for most muscle, which was significantly better than scaled-generic model for most muscle. This method open the way to a quick modeling method for gait analysis based on biplanar radiography.