Cost function criteria using muscle synergies: Exploring the potential of muscle synergy hypothesis.

BACKGROUND AND OBJECTIVE: Solving the redundant optimization problem for human muscles depends on the cost function. Choosing the appropriate cost function helps to address a specific problem. Muscle synergies are currently limited to those obtained by electromyography. Furthermore, debate continues regarding whether muscle synergy is derived or real. This study proposes new cost functions based on the muscle synergy hypothesis for solving the optimal muscle force output problem through musculoskeletal modeling. METHODS: We propose two new computational cost functions involving muscle synergies, which are extracted from muscle activations predicted by musculoskeletal modelling rather than electromyography. In this study, we constructed a musculoskeletal model for simulation using the "Grand Challenge Competition to Predict In Vivo Knee Loads" dataset. Muscle synergies were obtained using non-negative matrix factorization. Two cost functions with muscle synergies were constructed by integrating the polynomial and min/max criterion. Two new functions were verified and validated in normal, smooth, and bouncy gaits. RESULTS: The muscle synergies based on normal gaits were classified into four modules. The cosine similarities of the first three modules were all >0.9. In the normal and smooth gaits, the forces in most muscles predicted using the two new functions were within three standard deviations of the root mean square error for electromyographic comparisons. Predicted muscle force curves using the four methods as well as characteristic points (i.e., time points in the gait cycle when the significant difference was observed between normal and bouncy gaits) were obtained to validate their predictive capabilities. CONCLUSIONS: This study constructed two new cost functions involving muscle synergies, verified and validated the ability, and explored the potential of muscle synergy hypothesis.

Muscle dynamics analysis by clustered categories during jogging in patients with anterior cruciate ligament deficiency.

BACKGROUND: Patients with anterior cruciate ligament (ACL) deficiency (ACLD) tend to have altered lower extremity dynamics. Little is known about the changes in dynamic function and activation during jogging in patients with ACLD. METHODS: Twenty patients with an injured ACL before ACL reconstruction (ACLD group) and nine healthy male volunteers (control group) were recruited. Each volunteer repeated the jogging experiment five times. Based on the experimental data measured, a musculoskeletal multibody dynamics model was employed to simulate the tibiofemoral joint dynamics during jogging. Eighteen muscles were used for analysis. The obtained dynamics data were used for clustering and curve difference analysis. RESULTS: The 18 muscles studied were divided into 3 categories. All the quadriceps, the soleus, the gastrocnemius, and the popliteus were classified as label 1. All the hamstrings were classified as label 2, and the sartorius muscles were classified as label 3. Among them, the classification of the short head of the biceps femoris was significantly different between the two groups (P < 0.001). The force curves of all 18 muscles and the between-group differences were studied according to clustered categories. Most muscle force in label 1 was approaching zero in the terminal stance phase, which was significantly lower than that in the control group (P < 0.05). The muscle force in label 2 had areas with significant differences in the stance phase. Muscle force in label 3 was significantly lower than that in the control group in the pre-swing phase. CONCLUSIONS: This study showed that there are various changes of muscle function and activation in patients with ACLD. Through clustering and curve analysis, the joint reactions and changes of different muscle forces in the gait cycle between the ACLD and control groups could be further clarified.

An Analysis of Trabecular Bone Structure Based on Principal Stress Trajectory.

To understand the mechanism of Wolff's law, a finite element analysis was performed for a human proximal femur, and the principal stress trajectories of the femur were extracted using the principal stress visualization method. The mechanism of Wolff's law was evaluated theoretically based on the distribution of the principal stress trajectories. Due to the dynamics of the loads, there was no one-to-one correspondence between the stress trajectories of the fixed load and the trabeculae in the cancellous architecture of the real bone. The trabeculae in the cancellous bone were influenced by the magnitude of the principal stress trajectory. Equivalent principal stress trajectories suitable for different load changes were proposed through the change in load cycle and compared with the anatomical structure of the femur. In addition, the three-dimensional distribution of the femoral principal stress trajectory was established, and the adaptability potential of each load was discussed. The principal stress visualization method could also be applied to bionic structure design.

Orthodontic craniofacial pattern diagnosis: cephalometric geometry and machine learning.

Efficient and reliable diagnosis of craniofacial patterns is critical to orthodontic treatment. Although machine learning (ML) is time-saving and high-precision, prior knowledge should validate its reliability. This study proposed a craniofacial ML diagnostic workflow base on a cephalometric geometric model through clinical verification. A cephalometric geometric model was established to determine the landmark location by analyzing 408 X-ray lateral cephalograms. Through geometric information and feature engineering, nine supervised ML algorithms were conducted for sagittal and vertical skeleton patterns. After dimension reduction, plane decision boundary and landmark contribution contours were depicted to demonstrate the diagnostic consistency and the consistency with clinical norms. As a result, multi-layer perceptron achieved 97.56% accuracy for sagittal, while linear support vector machine reached 90.24% for the vertical. Sagittal diagnoses showed average superiority (91.60 ± 5.43)% over the vertical (82.25 ± 6.37)%, where discriminative algorithms exhibited more steady performance (93.20 ± 3.29)% than the generative (85.98 ± 9.48)%. Further, the Kruskal-Wallis H test was carried out to explore statistical differences in diagnoses. Though sagittal patterns had no statistical difference in diagnostic accuracy, the vertical showed significance. All aspects of the tests indicated that the proposed craniofacial ML workflow was highly consistent with clinical norms and could supplement practical diagnosis.

Effect of remaining pericervical dentin on biomechanical behavior of endocrown-restored molars with different materials: Three-dimensional finite element and Weibull analyses.

To evaluate the effect of remaining pericervical dentin (PCD) on the biomechanical behavior of endocrown-restored molars with different materials, six three-dimensional finite element (FE) models were reconstructed with different thicknesses and heights of pulp-chamber lateral dentinal wall (PCLDW). IPS Empress 2, In-Ceram Zirconia, and Lava Ultimate were selected as the materials. Compared with the Lava Ultimate FE models, the maximum tensile stress in the FE models using ceramics was higher in the endocrown and lower in the PCD surrounding it, and the overall failure probabilities with different PCLDW thicknesses and heights were similar, ranging from 9.8% to 12.9% under the normal lateral masticatory force, which were lower than the FE models using Lava Ultimate (ranging from 13.4% to 15.1%). Considering the bonding properties of ceramics, endocrown-restored molars using etchable lithium disilicate-reinforced glass ceramic exhibit superior longevity due to the stress shielding effect, regardless of the thickness and height of PCLDW.

Leveraging Multivariable Linear Regression Analysis to Identify Patients with Anterior Cruciate Ligament Deficiency Using a Composite Index of the Knee Flexion and Muscle Force.

Patients with anterior cruciate ligament (ACL) deficiency (ACLD) tend to have altered lower extremity kinematics and dynamics. Clinical diagnosis of ACLD requires more objective and convenient evaluation criteria. Twenty-five patients with ACLD before ACL reconstruction and nine healthy volunteers were recruited. Five experimental jogging data sets of each participant were collected and calculated using a musculoskeletal model. The resulting knee flexion and muscle force data were analyzed using a t-test for characteristic points, which were the time points in the gait cycle when the most significant difference between the two groups was observed. The data of the characteristic points were processed with principal component analysis to generate a composite index for multivariable linear regression. The accuracy rate of the regression model in diagnosing patients with ACLD was 81.4%. This study demonstrates that the multivariable linear regression model and composite index can be used to diagnose patients with ACLD. The composite index and characteristic points can be clinically objective and can be used to extract effective information quickly and conveniently.

Finite element analysis of stress in oral mucosa and titanium mesh interface.

BACKGROUND: The stiffness of titanium mesh is a double-blade sword to repair larger alveolar ridges defect with excellent space maintenance ability, while invade the surrounding soft tissue and lead to higher mesh exposure rates. Understanding the mechanical of oral mucosa/titanium mesh/bone interface is clinically meaningful. In this study, the above relationship was analyzed by finite elements and verified by setting different keratinized tissue width in oral mucosa. METHODS: Two three-dimensional finite element models were constructed with 5 mm keratinized tissue in labial mucosa (KM cases) and 0 mm keratinized tissue in labial mucosa (LM cases). Each model was composed of titanium mesh, titanium screws, graft materials, bone, teeth and oral mucosa. After that, a vertical (30 N) loadings were applied from both alveolar ridges direction and labial mucosa direction to stimulate the force from masticatory system. The displacements and von Mises stress of each element at the interfaces were analyzed. RESULTS: Little displacements were found for titanium mesh, titanium screws, graft materials, bone and teeth in both LM and KM cases under different loading conditions. The maximum von Mises stress was found around the lingual titanium screw insertion place for those elements in all cases. The keratinized tissue decreased the displacement of oral mucosa, decreased the maximum von Mises stress generated by an alveolar ridges direction load, while increased those stress from labial mucosa direction load. Only the von Mises stress of the KM cases was all lower than the tensile strength of the oral mucosa. CONCLUSION: The mucosa was vulnerable under the increasing stress generated by the force from masticatory system. The adequate buccal keratinized mucosa width are critical factors in reducing the stress beyond the titanium mesh, which might reduce the titanium exposure rate.

Finite element analysis on pushing the molar backwards using invisible aligner with different migration displacement.

PURPOSE: This paper examines the biomechanical mechanism behind the effect of the invisible aligner technique on tooth movement processes. METHODS: To compare the effects of different target positions on tooth movement and the periodontal ligament (PDL), two kinds of aligners were designed to provide displacements of 0.2 mm (Model A) and 0.3 mm (Model B). Different displacements of the maxillary second molar were simulated using the finite element (FE) method. RESULTS: The results of numerical simulations showed that the maximum stress was in the PDL of the distal surface and the palatal surface. The stress of the PDL in Model B was larger than Model A, with the displacement of the second molar 0.027 mm in Model A, by 44.9% lesser than that in Model B. CONCLUSIONS: The aligner that provided a displacement of 0.2 mm was more suitable for pushing the second molar backward in the initial stage. During the tooth movement processes, the displacement of the crown was larger than that of the root and the displacement decreased gradually from the crown to the root. In addition, the displacement and rotation of teeth during orthodontic treatment were measured and analysed.

Effect of partial restorative treatment on stress distributions in non-carious cervical lesions: a three-dimensional finite element analysis.

BACKGROUND: Partial restoration combined with periodontal root coverage surgery can be applied to the treatment of non-carious cervical lesions (NCCLs) accompanied with gingival recessions in clinical practice. However, the feasibility of NCCL partial restorative treatment from a biomechanical perspective remains unclear. This study aimed to investigate the effect of partial restorations on stress distributions in the NCCLs of mandibular first premolars via three-dimensional finite element analysis. METHODS: Three-dimensional finite element models of buccal wedge-shaped NCCLs in various locations of a defected zenith (0 mm, 1 mm, and 2 mm) were constructed and divided into three groups (A, B, and C). Three partially restored NCCL models with different locations of the lower restoration border (1 mm, 1.5 mm, and 2 mm), and one completely restored NCCL model were further constructed for each group. The following restorative materials were used in all restoration models: composite resin (CR), glass-ionomer cement (GIC), and mineral trioxide aggregate (MTA). The first principal stress distributions under buccal oblique loads of 100 N were analyzed. Restoration bond failures were also evaluated based on stress distributions at dentin-restoration interfaces. RESULTS: When the partial restoration fully covered the defected zenith, the first principal stress around the zenith decreased and the maximum tensile stress was concentrated at the lower restoration border. When the partial restoration did not cover the defected zenith, the first principal stress distribution patterns were similar to those in unrestored models, with the maximum tensile stress remaining concentrated at the zenith. As the elastic modulus of the restorative material was altered, the stress distributions at the interface were not obviously changed. Restoration bond failures were not observed in CR, but occurred in GIC and MTA in most models. CONCLUSIONS: Partial restorations that fully covered defected zeniths improved the stress distributions in NCCLs, while the stress distributions were unchanged or worsened under other circumstances. CR was the optimal material for partial restorations compared to GIC and MTA.

Identification of Kinetic Abnormalities in Male Patients after Anterior Cruciate Ligament Deficiency Combined with Meniscal Injury: A Musculoskeletal Model Study of Lower Limbs during Jogging.

There is little known about kinetic changes in anterior cruciate ligament deficiency (ACLD) combined with meniscal tears during jogging. Therefore, 29 male patients with injured ACLs and 15 healthy male volunteers were recruited for this study to investigate kinetic abnormalities in male patients after ACL deficiency combined with a meniscal injury during jogging. Based on experimental data measured by an optical tracking system, a subject-specific musculoskeletal model was employed to estimate the tibiofemoral joint kinetics during jogging. Between-limb and interpatient differences were compared by the analysis of variance. The results showed that decreased knee joint forces and moments of both legs in ACLD patients were detected during the stance phase compared to the control group. Meanwhile, compared with ACLD knees, significantly fewer contact forces and flexion moments in ACLD combined with lateral and medial meniscal injury groups were found at the mid-stance, and ACLD with medial meniscal injury group showed a lower axial moment in the loading response (p < 0.05). In conclusion, ACLD knees exhibit reduced tibiofemoral joint forces and moments during jogging when compared with control knees. A combination of meniscus injuries in the ACLD-affected side exhibited abnormal kinetic alterations at the loading response and mid-stance phase.

Relationship between high intra-abdominal pressure and compliance of the pelvic floor support system in women without pelvic organ prolapse: A finite element analysis.

Previous studies mainly focused on the relationship between the size of the prolapse and injury to the supporting tissues, but the strain and stress distributions of the supporting tissues as well as high-risk areas of injury are still unknown. To further investigate the effect of supporting tissues on organs and the interactions between organs, this study focused on the relationship between high intra-abdominal pressure and the compliance of the pelvic floor support system in a normal woman without pelvic organ prolapse (POP), using a finite element model of the whole pelvic support system. A healthy female volunteer (55 years old) was scanned using magnetic resonance imaging (MRI) during rest and Valsalva maneuver. According to the pelvic structure contours traced by a gynecologist and anatomic details measured from dynamic MRI, a finite element model of the whole pelvic support system was established, including the uterus, vagina with cavity, cardinal and uterosacral ligaments, levator ani muscle, rectum, bladder, perineal body, pelvis, and obturator internus and coccygeal muscles. This model was imported into ANSYS software, and an implicit iterative method was employed to simulate the biomechanical response with increasing intra-abdominal pressure. Stress and strain distributions of the vaginal wall showed that the posterior wall was more stable than the anterior wall under high intra-abdominal pressure. Displacement at the top of the vagina was larger than that at the bottom, especially in the anterior-posterior direction. These results imply potential injury areas with high intra-abdominal pressure in non-prolapsed women, and provide insight into clinical managements for the prevention and surgical repair plans of POP.

Three-dimensional radiological anatomy of condyle trabecular bone based on a Volume-of-Interest analysis.

OBJECTIVES: Three-dimensional radiological anatomic characteristics of condyle trabeculae was obtained quantitatively based on a volume-of-interest (VOI) analysis. METHODS: Nine human mandibular condyle specimens were scanned by micro-computed tomography (micro-CT). A total of 34 VOIs were selected from each condyle specimen, which were divided into six layers and four parts to analyze the morphological characteristics of trabeculae based on cylindrical VOIs with a diameter and height of 2 mm. One-way analysis of variance was used to compare the regional differences of morphological parameters among each layer and part. RESULTS: Values for bone mineral density, bone volume/total volume, trabecular thickness, and trabecular bone number were greater in the anterior part compared with the posterior part; and the lateral part was larger than the medial part in the first, second, and third layers, while the medial part was larger in the fourth and fifth layers; these values in the first and sixth layers were much larger, while those in the third and fourth layers were smaller. Bone surface area/bone volume, trabecular spacing, and trabecular bone pattern factor were larger in the posterior part than in the anterior part; and the lateral part was larger than the medial part in the fourth and fifth layers, while the medial part was larger in the first and second layers. CONCLUSIONS: The morphological distribution of VOIs was anisotropic within trabecular bone of human mandibular condyles. The upper and lower ends of trabecular bone were much more compact, with higher bone density, trabecular thickness, and trabecular number than in the middle layers.

Influence of different diameter reductions in the labial neck region on the stress distribution around custom-made root-analogue implants.

This study was designed to investigate the influence of diameter reductions on the stress distribution around root-analogue implants via 3D finite element analysis. Four root-analogue implant models with different diameter reductions (0, 1, 2, or 3 mm), a traditional threaded implant and congruent bone models were created through reverse engineering. A 100-N force was applied parallel with and in a 45° angle to the implant axis, respectively. The stress concentration in the labial neck area around implants with 1-2 mm diameter reduction was lower than seen with no reduction. When the implant diameter was reduced by 3 mm, there were obvious stress concentrations in both implant and bone (the maximum stress was 206  and 111 MPa, respectively). In other groups, the maximum stress was 65.1 MPa in the bone and 108 MPa in the implant. Additionally, the stress concentration in the bone around the root-analogue implant when the implant diameter was reduced by 0-2 mm (maximum stress of 65.1 MPa) was obviously smaller than that around the traditional implant (maximum stress 130.4 MPa). Reducing the diameter of maxillary central incisor root-analogue implants by up to 2 mm next to the labial cortical bone could help disperse stress.

A Finite Element Analysis Model is Suitable for Biomechanical Analysis of Orbital Development.

ABSTRACT: The authors investigated orbital bone development in congenital microphthalmia (CM) using a three-dimensional finite element analysis model, after the orbital dimension deficiency was improved with a self-inflating hydrogel expander implant.Data of a 2-year-old male CM patient were examined. The orbital structure, eyeball, eye muscles, and self-inflating hydrogel expander were constructed according to computed tomography examination data. The effects of self-expanding spherical hydrogel at various locations in the muscle cone on orbital bone development were examined using 3-mL self-expanding expanders placed at shallow (model 1: 2 mm depth) and deep (model 2: 8 mm depth) muscle cone positions. This model simulated the hydrogel expansion process; the orbital bone biomechanics and radial displacement nephograms were obtained when the hydrogel volume increased 3, 5, 7, and 9 times and analyzed.The orbital bone biomechanics were concentrated at the medial orbital wall center, gradually spreading to the whole orbital wall. Biomechanics and radial displacement of the inferior temporal and lateral distal orbital wall region were the most significant, and were more significant in model 1 than in model 2.Finite element analysis is suitable for the biomechanical analysis of orbital development in CM. The shallow position inside CM patients' muscle cone is the optimal site for hydrogel implantation.

The Biomechanical Response of the Cornea in Orthokeratology.

Orthokeratology has been widely used to control myopia, but the mechanism is still unknown. To further investigate the underlying mechanism of corneal reshaping using orthokeratology lenses via the finite element method, numerical models with different corneal curvatures, corneal thicknesses, and myopia reduction degrees had been developed and validated to simulate the corneal response and quantify the changes in maximum stress in the central and peripheral corneal areas during orthokeratology. The influence of the factors on corneal response had been analyzed by using median quantile regression. A partial eta squared value in analysis of variance models was established to compare the effect size of these factors. The results showed central and peripheral corneal stress responses changed significantly with increased myopia reduction, corneal curvature, and corneal thickness. The target myopia reduction had the greatest effect on the central corneal stress value (partial eta square = 0.9382), followed by corneal curvature (partial eta square = 0.5650) and corneal thickness (partial eta square = 0.1975). The corneal curvature had the greatest effect on the peripheral corneal stress value (partial eta square = 0.5220), followed by myopia reduction (partial eta square = 0.2375) and corneal thickness (partial eta square = 0.1972). In summary, the biomechanical response of the cornea varies significantly with the change in corneal conditions and lens designs. Therefore, the orthokeratology lens design and the lens fitting process should be taken into consideration in clinical practice, especially for patients with high myopia and steep corneas.

Influence of sagittal root positions on the stress distribution around custom-made root-analogue implants: a three-dimensional finite element analysis.

BACKGROUND: Stress concentration may cause bone resorption even lead to the failure of implantation. This study was designed to investigate whether a certain sagittal root position could cause stress concentration around maxillary anterior custom-made root-analogue implants via three-dimensional finite element analysis. METHODS: The von Mises stresses in the bone around implants in different groups were compared by finite element analysis. Six models were constructed and divided into two groups through Geomagic Studio 2012 software. The smooth group included models of unthreaded custom-made implants in Class I, II or III sagittal root positions. The threaded group included models of reverse buttress-threaded implants in the three positions. The von Mises stress distributions and the range of the stresses under vertical and oblique loads of 100 N were analyzed through ANSYS 16.0 software. RESULTS: Stress concentrations around the labial lamella area were more prominent in the Class I position than in the Class II and Class III positions under oblique loading. Under vertical loading, the most obvious stress concentration areas were the labial lamella and palatal apical areas in the Class I and Class III positions, respectively. Stress was relatively distributed in the labial and palatal lamellae in the Class II position. The maximum von Mises stresses in the bone around the custom-made root-analogue implants in this study were lower than around traditional implants reported in the literature. The maximum von Mises stresses in this study were all less than 25 MPa in cortical bone and less than 6 MPa in cancellous bone. Additionally, compared to the smooth group, the threaded group showed lower von Mises stress concentration in the bone around the implants. CONCLUSIONS: The sagittal root position affected the von Mises stress distribution around custom-made root-analogue implants. There was no certain sagittal root position that could cause excessive stress concentration around the custom-made root-analogue implants. Among the three sagittal root positions, the Class II position would be the most appropriate site for custom-made root-analogue implants.

Finite Element Analysis of Two- and Three-Dimensional Fixation in Treating Mandibular Symphyseal Fracture Combined With Bilateral Condylar Intracapsular Fractures.

ABSTRACT: The aim of this study was to compare through finite element analysis two- and three-dimensional (2D and 3D) fixation in the treatment of mandibular symphyseal fracture combined with bilateral condylar intracapsular fractures. The authors created 2 fixation models for the above fracture, and analyzed the stress and displacement in the mandible and fixation materials under 3 loading conditions. The von Mises stress of the mandible and plates peaked during lateral occlusion, and was lowest during central occlusion. In all conditions, stresses in the fixation materials did not exceed the yield stress of titanium. The inferior border of the symphyseal fracture segments showed opposing displacements, and the mandible tended to widen in the 2D fixation model. However, the fracture displacement did not exceed 150 µm for either fixation method. The results suggested that after well reduction and fixation of condylar intracapsular fractures, either 2D or 3D fixation for symphyseal fracture can provide adequately strong fixation. Compared with 2D fixation, 3D fixation has more advantages in controlling the mandibular width and preventing the fixation materials from enduring excessive stress.

Use of the normalcy index for the assessment of abnormal gait in the anterior cruciate ligament deficiency combined with meniscus injury.

The normalcy index (NI) has been implemented by several studies as a simple index for quantitatively analyzing diffident gait abnormalities, such as children with cerebral palsy and idiopathic toe-walkers. However, whether the NI can be used in anterior cruciate ligament (ACL) deficiency with different types of meniscus injuries or not, has not been reported yet. In this study, 25 patients who combined different types of ACL and meniscus injuries were evaluated by the NI analysis. 12 healthy subjects were used to define the normal range of NI. The result showed that NI values of patients were significantly larger than the control group (P < 0.05). Meanwhile, the tendency of increasing NI values associated with increasing pathology were significant with only 5 subjects in the smallest group (Jonkheere-Terpsta test: P < 0.001). These results indicated that the NI was a concise yet effective tool to evaluate combined ACL and meniscus injury patients. Increasing severity degree of meniscus tears in ACL rupture patients is corresponded to increasing NI values. It also demonstrates that the proposed NI can be applied as a robustness factor to detect the discrepancy between healthy and patient subjects clinically, and has the potential in the quantitative evaluation of pre- or post-surgery and rehabilitation.

Effect of central retainer shape and abduction angle during preparation of teeth on dentin and cement layer stress distributions in endocrown-restored mandibular molars.

To study the effect of central retainer shape and abduction angle during tooth preparation on stress distribution in endocrown-restored molars via finite element (FE) analysis, we constructed five FE models with different central retainer shapes and abduction angles. Under an oblique load, the distributions of maximum tensile stress in cervical dentin around the endocrown and on the cement layer, as well as maximum shear stress on the cement layer, were more balanced in the FE model in which the central retainer shape was generated based on the anatomical form of the pulp chamber. Moreover, there were no differences in stress distributions among FE models with different abduction angles. Therefore, the shape of the central retainer should be designed on the basis of the anatomical form of the pulp chamber; abduction angle during tooth preparation does not influence the repair effect of endocrown-restored mandibular molars.

Biomechanical comparison between the custom-made mandibular condyle prosthesis and total temporomandibular joint prosthesis in finite element analysis.

PURPOSE: The purpose of this study was to compare the biomechanical behavior of the custom-made mandibular condyle prosthesis and total TMJ prosthesis. METHODS: Three models of one beagle dog, the condyle prosthesis (Model 1, replacing the right condyle only), the TMJ prosthesis (model 2, replacing the whole right TMJ) and the intact TMJ (model 3) were established, and the mechanical responses under muscle forces loading were analyzed using finite element method. RESULTS: Models 1 and 3 had the similar stress distribution on the right disc, which suggested that the condyle prosthesis did not change the disc stress so much when the muscle forces were applied. The stress of the right TMJ prosthesis in Model 2 was larger than both Models 1 and 3, and the stress of the contralateral TMJ reduced by 12% in Model 2. The anterior border of the condyle seemed to be a stress concentration region, not only for the intact condyle, but also for the condyle prosthesis and the total TMJ prosthesis. CONCLUSIONS: The total TMJ prosthesis changed the biomechanical balance of the bilateral TMJ. When the condyle prosthesis iss applied, the custom-made profile is recommended.