A finite element analysis on the indication for extracting partially impacted mandibular third molars considering mandibular trauma.

BACKGROUND: Patients presenting with partially impacted lower third molars (M3) have a higher likelihood of experiencing angle fractures while simultaneously decreasing the risk of condylar fractures. However, the specific biomechanical mechanism responsible for this occurrence remains unclear. Moreover, there is an ongoing debate regarding whether the removal of M3s might actually increase the risk of condylar fractures. This study aimed to evaluate how the presence of M3s influences mandibular fractures resulting from blows to the symphysis and lateral mandibular body, and to determine the indication for extracting M3s in such cases. METHODS: Models of the mandible with a partially M3-impacted model (M3I), M3-extracted model (M3E), and M3-absent model (M3A) were generated using a computer. A traumatic blown force of 2000 N was applied to the symphysis and the right body of the mandible. Von Mises and principal stresses were analyzed, and failure indexes were determined. Two cases of mandibular linear fractures were chosen for model verification and interpretation. RESULTS: When force was applied to the symphysis, the condylar region exhibited the highest stress levels, while stress in the mandibular angle region was much less regardless of the M3 state. On applying the force to the right mandibular body, stress in the condylar region decreased while stress in the mandibular body increased, especially in the blown regions. Impacted tooth or cavity formation post-M3 extraction led to uneven stress distribution on the blown side of the mandible, increasing the risk of mandibular angle fractures. In cases where M3 was absent or the extraction socket had healed, stress from lateral traumatic blown force was evenly distributed along both the inner and outer oblique lines of the mandible, thereby reducing the risk of mandibular fractures. CONCLUSIONS: The reduced risk of condylar fractures in patients with partially impacted lower M3s and mandibular angle fractures is mainly due to lateral blows on the mandible, which generate less stress in the condylar region than blows on the mandibular symphysis, rather than being caused by the M3 itself. Extraction of the lower M3 can decrease the risk of mandibular fractures, with a minor influence on condylar fractures.

Influence of rotator cuff preload on fracture configuration in proximal humerus fractures: a proof of concept for fracture simulation.

INTRODUCTION: In regard of surgical training, the reproducible simulation of life-like proximal humerus fractures in human cadaveric specimens is desirable. The aim of the present study was to develop a technique that allows simulation of realistic proximal humerus fractures and to analyse the influence of rotator cuff preload on the generated lesions in regards of fracture configuration. MATERIALS AND METHODS: Ten cadaveric specimens (6 left, 4 right) were fractured using a custom-made drop-test bench, in two groups. Five specimens were fractured without rotator cuff preload, while the other five were fractured with the tendons of the rotator cuff preloaded with 2 kg each. The humeral shaft and the shortened scapula were potted. The humerus was positioned at 90° of abduction and 10° of internal rotation to simulate a fall on the elevated arm. In two specimens of each group, the emergence of the fractures was documented with high-speed video imaging. Pre-fracture radiographs were taken to evaluate the deltoid-tuberosity index as a measure of bone density. Post-fracture X-rays and CT scans were performed to define the exact fracture configurations. Neer's classification was used to analyse the fractures. RESULTS: In all ten cadaveric specimens life-like proximal humerus fractures were achieved. Two III-part and three IV-part fractures resulted in each group. The preloading of the rotator cuff muscles had no further influence on the fracture configuration. High-speed videos of the fracture simulation revealed identical fracture mechanisms for both groups. We observed a two-step fracture mechanism, with initial impaction of the head segment against the glenoid followed by fracturing of the head and the tuberosities and then with further impaction of the shaft against the acromion, which lead to separation of the tuberosities. CONCLUSION: A high energetic axial impulse can reliably induce realistic proximal humerus fractures in cadaveric specimens. The preload of the rotator cuff muscles had no influence on initial fracture configuration. Therefore, fracture simulation in the proximal humerus is less elaborate. Using the presented technique, pre-fractured specimens are available for real-life surgical education. LEVEL OF EVIDENCE: III.

Applying Machine Learning to Finger Movements Using Electromyography and Visualization in Opensim.

Electromyographic signals have been used with low-degree-of-freedom prostheses, and recently with multifunctional prostheses. Currently, they are also being used as inputs in the human-computer interface that controls interaction through hand gestures. Although there is a gap between academic publications on the control of an upper-limb prosthesis developed in laboratories and its service in the natural environment, there are attempts to achieve easier control using multiple muscle signals. This work contributes to this, using a database and biomechanical simulation software, both open access, to seek simplicity in the classifiers, anticipating their implementation in microcontrollers and their execution in real time. Fifteen predefined finger movements of the hand were identified using classic classifiers such as Bayes, linear and quadratic discriminant analysis. The idealized movements of the database were modeled with Opensim for visualization. Combinations of two preprocessing methods-the forward sequential selection method and the feature normalization method-were evaluated to increase the efficiency of these classifiers. The statistical methods of cross-validation, analysis of variance (ANOVA) and Duncan were used to validate the results. Furthermore, the classifier with the best recognition result was redesigned into a new feature space using the sparse matrix algorithm to improve it, and to determine which features can be eliminated without degrading the classification. The classifiers yielded promising results-the quadratic discriminant being the best, achieving an average recognition rate for each individual considered of 96.16%, and with 78.36% for the total sample group of the eight subjects, in an independent test dataset. The study ends with the visual analysis under Opensim of the classified movements, in which the usefulness of this simulation tool is appreciated by revealing the muscular participation, which can be useful during the design of a multifunctional prosthesis.

Finite element method for estimation of applanation force and to study the influence of intraocular pressure of eye on tonometry.

PURPOSE: Discover the associations of force of applanation on the eye with the plunging depth of the cornea and quantify them. The results will be utilized as the feedback parameter in the new prototype development of eye care instruments as additional force may damage the internal structure of the eye or may result in erroneous output. METHOD: A finite element-based eye model is designed utilizing the actual dimensions of the human eye. A standardized tonometer is designed and the simulation is carried out at predetermined deformation of the cornea to find the force of applanation on the cornea during tonometry. Adding on, the influence of IOP during tonometry is analyzed for a range of plunging depths of the cornea. RESULTS: The graphical results inferred the linear relation between the force of applanation with the deformation of the cornea and the results are quantified. The resulting deformation and stress plot of FEM based simulation approach is analyzed and observations regarding deformations and stress are made. CONCLUSION: The human eye is successfully developed and also computed force on the cornea during tonometry is validated. The inference drawn from the deformation plot and stress plot is that the junction of cornea-sclera along with cornea-tonometer periphery undergo maximum deformation and experiences the highest stress compared to other areas of the eye while during tonometry.

Musculoskeletal biomechanics of patients with or without adjacent segment degeneration after spinal fusion.

STUDY DESIGN: A retrospective, single center, case-control study was performed. OBJECTIVE: The present study employed patient-specific biomechanical modeling to find potential biomechanical differences after spinal fusion at L4/5 in patients with and without subsequent development of adjacent segment disease (ASD). METHODS: The study population comprised patients who underwent primary spinal fusion at L4/5 and were either asymptomatic during > 4 years of follow-up (CTRL; n = 18) or underwent revision surgery for ASD at L3/4 (n = 20). Landmarks were annotated on preoperative and follow-up lateral radiographs, and specific musculoskeletal models were created using a custom-built modeling pipeline. Simulated spinal muscle activation and lumbar intervertebral shear loads in unfused segments were analyzed in upright standing and forward flexion. Differences between the pre- and postoperative conditions were computed for each patient. RESULTS: The average postoperative muscle activity in the upright standing posture was 88.4% of the preoperative activity in the CTRL group (p <  0.0001), but did not significantly change from pre- to postoperatively in the ASD group (98.0%). The average shear load magnitude at the epifusional joint L3/4 during upright standing increased from pre- to postoperatively in the ASD group (+ 3.9 N, +/- 17.4 (n = 18)), but decreased in the CTRL group (- 4.6 N, +/- 23.3 (n = 20); p <  0.001). CONCLUSION: Patient-specific biomechanical simulation revealed that spinal fusion surgery resulted in greater shear load magnitude and muscle activation and therefore greater forces at the epifusional segment in those with ASD compared with those without ASD. This is a first report of patient-specific disc load and muscle force calculation with predictive merits for ASD.

The Biomechanical Mechanism of Upper Airway Collapse in OSAHS Patients Using Clinical Monitoring Data during Natural Sleep.

Obstructive sleep apnea hypopnea syndrome (OSAHS) is a common sleep disorder characterized by repeated pharyngeal collapse with partial or complete obstruction of the upper airway. This study investigates the biomechanics of upper airway collapse of OSASH patients during natural sleep. Computerized tomography (CT) scans and data obtained from a device installed on OSASH patients, which is comprised of micro pressure sensors and temperature sensors, are used to develop a pseudo three-dimensional (3D) finite element (FE) model of the upper airway. With consideration of the gravity effect on the soft palate while patients are in a supine position, a fluid-solid coupling analysis is performed using the FE model for the two respiratory modes, eupnea and apnea. The results of this study show that the FE simulations can provide a satisfactory representation of a patient's actual respiratory physiological processes during natural sleep. The one-way valve effect of the soft palate is one of the important mechanical factors causing upper airway collapse. The monitoring data and FE simulation results obtained in this study provide a comprehensive understanding of the occurrence of OSAHS and a theoretical basis for the individualized treatment of patients. The study demonstrates that biomechanical simulation is a powerful supplementation to clinical monitoring and evaluation.

Quantal biomechanical effects in speech postures of the lips.

The unique biomechanical and functional constraints on human speech make it a promising area for research investigating modular control of movement. The present article illustrates how a modular control approach to speech can provide insights relevant to understanding both motor control and observed variation across languages. We specifically explore the robust typological finding that languages produce different degrees of labial constriction using distinct muscle groupings and concomitantly distinct lip postures. Research has suggested that these lip postures exploit biomechanical regions of nonlinearity between neural activation and movement, also known as quantal regions, to allow movement goals to be realized despite variable activation signals. We present two sets of computer simulations showing that these labial postures can be generated under the assumption of modular control and that the corresponding modules are biomechanically robust: first to variation in the activation levels of participating muscles, and second to interference from surrounding muscles. These results provide support for the hypothesis that biomechanical robustness is an important factor in selecting the muscle groupings used for speech movements and provide insight into the neurological control of speech movements and how biomechanical and functional constraints govern the emergence of speech motor modules. We anticipate that future experimental work guided by biomechanical simulation results will provide new insights into the neural organization of speech movements.NEW & NOTEWORTHY This article provides additional evidence that speech motor control is organized in a modular fashion and that biomechanics constrain the kinds of motor modules that may emerge. It also suggests that speech can be a fruitful domain for the study of modularity and that a better understanding of speech motor modules will be useful for speech research. Finally, it suggests that biomechanical modeling can serve as a useful complement to experimental work when studying modularity.

Biomechanical Optimization-Based Planning of Periacetabular Osteotomy.

Modern computerized planning tools for periacetabular osteotomy (PAO) use either morphology-based or biomechanics-based methods. The latter rely on estimation of peak contact pressures and contact areas using either patient-specific or constant thickness cartilage models. We performed a finite element analysis investigating the optimal reorientation of the acetabulum in PAO surgery based on simulated joint contact pressures and contact areas using patient-specific cartilage model. Furthermore we investigated the influences of using patient-specific cartilage model or constant thickness cartilage model on the biomechanical simulation results. Ten specimens with hip dysplasia were used in this study. Image data were available from CT arthrography studies. Bone models were reconstructed. Mesh models for the patient-specific cartilage were defined and subsequently loaded under previously reported boundary and loading conditions. Peak contact pressures and contact areas were estimated in the original position. Afterward we used validated preoperative planning software to change the acetabular inclination by an increment of 5° and measured the lateral center-edge angle (LCE) at each reorientation position. The position with the largest contact area and the lowest peak contact pressure was defined as the optimal position. In order to investigate the influence of using patient-specific cartilage model or constant thickness cartilage model on the biomechanical simulation results, the same procedure was repeated with the same bone models but with a cartilage mesh of constant thickness. Comparison of the peak contact pressures and the contact areas between these two different cartilage models showed that good correlation between these two cartilage models for peak contact pressures (r = 0.634 ∈[0.6, 0.8], p < 0.001) and contact areas (r = 0.872 > 0.8, p < 0.001). For both cartilage models, the largest contact areas and the lowest peak pressures were found at the same position. Our study is the first study comparing peak contact pressures and contact areas between patient-specific and constant thickness cartilage models during PAO planning. Good correlation for these two models was detected. Computer-assisted planning with FE modeling using constant thickness cartilage models might be a promising PAO planning tool when a conventional CT is available.

Real-time prediction of postoperative spinal shape with machine learning models trained on finite element biomechanical simulations.

PURPOSE: Adolescent idiopathic scoliosis is a chronic disease that may require correction surgery. The finite element method (FEM) is a popular option to plan the outcome of surgery on a patient-based model. However, it requires considerable computing power and time, which may discourage its use. Machine learning (ML) models can be a helpful surrogate to the FEM, providing accurate real-time responses. This work implements ML algorithms to estimate post-operative spinal shapes. METHODS: The algorithms are trained using features from 6400 simulations generated using the FEM from spine geometries of 64 patients. The features are selected using an autoencoder and principal component analysis. The accuracy of the results is evaluated by calculating the root-mean-squared error and the angle between the reference and predicted position of each vertebra. The processing times are also reported. RESULTS: A combination of principal component analysis for dimensionality reduction, followed by the linear regression model, generated accurate results in real-time, with an average position error of 3.75 mm and orientation angle error below 2.74 degrees in all main 3D axes, within 3 ms. The prediction time is considerably faster than simulations based on the FEM alone, which require seconds to minutes. CONCLUSION: It is possible to predict post-operative spinal shapes of patients with AIS in real-time by using ML algorithms as a surrogate to the FEM. Clinicians can compare the response of the initial spine shape of a patient with AIS to various target shapes, which can be modified interactively. These benefits can encourage clinicians to use software tools for surgical planning of scoliosis.

Evaluation of women's aging influence on sperm passage inside the fallopian tube using 3D dynamic mechanical modeling.

The fallopian tubes play an important role in human fertility by facilitating the spermatozoa passage to the oocyte as well as later actively facilitating the fertilized oocyte transportation to the uterus cavity. The fallopian tubes undergo changes involving biological, physical, and morphological processes due to women aging, which may impair fertility. Here, we have modelled fallopian tubes of women at different ages and evaluated the chances of normal and pathological sperm cells reaching the fertilization site, the ampulla. By utilizing a unique combination of simulative tools, we implemented dynamic three-dimensional (3D) detailed geometrical models of many normal and pathological sperm cells swimming together in 3D geometrical models of three fallopian tubes associated with different women's age groups. By tracking the sperm cell swim, we found that for all age groups, the number of normal sperm cells in the ampulla is the largest, compared with the pathological sperm cells. On the other hand, the number of normal sperm cells in the fertilization site decreases due to the morphological and mechanical changes that occur in the fallopian tube with age. Moreover, in older ages, the normal sperm cells swim with lower velocities and for shorter distances inside the ampulla toward the ovary. Thus, the changes that the human fallopian tube undergoes due to women's aging have a significant influence on the human sperm cell motility. Our model of sperm cell motility through the fallopian tube in relation to the woman's age morphological changes provides a new scope for the investigation and treatment of diseases and infertility cases associated with aging, as well as a potential personalized medicine tool for evaluating the chances of a natural fertilization per specific features of a man's sperm and a woman's reproductive system.

Experimental and virtual testing of bone-implant systems equipped with the AO Fracture Monitor with regard to interfragmentary movement.

Introduction: The management of fractured bones is a key domain within orthopedic trauma surgery, with the prevention of delayed healing and non-unions forming a core challenge. This study evaluates the efficacy of the AO Fracture Monitor in conjunction with biomechanical simulations to better understand the local mechanics of fracture gaps, which is crucial for comprehending mechanotransduction, a key factor in bone healing. Through a series of experiments and corresponding simulations, the study tests four hypotheses to determine the relationship between physical measurements and the predictive power of biomechanical models. Methods: Employing the AO Fracture Monitor and Digital Image Correlation techniques, the study demonstrates a significant correlation between the surface strain of implants and interfragmentary movements. This provides a foundation for utilizing one-dimensional AO Fracture Monitor measurements to predict three-dimensional fracture behavior, thereby linking mechanical loading with fracture gap dynamics. Moreover, the research establishes that finite element simulations of bone-implant systems can be effectively validated using experimental data, underpinning the accuracy of simulations in replicating physical behaviors. Results and Discussion: The findings endorse the combined use of monitoring technologies and simulations to infer the local mechanical conditions at the fracture site, offering a potential leap in personalized therapy for bone healing. Clinically, this approach can enhance treatment outcomes by refining the assessment precision in trauma trials, fostering the early detection of healing disturbances, and guiding improvements in future implant design. Ultimately, this study paves the way for more sophisticated patient monitoring and tailored interventions, promising to elevate the standard of care in orthopedic trauma surgery.

The effect of the 2-UPS/RR ankle rehabilitation robot with coupling biomechanical model on muscle behaviors.

With the popularization of biomechanical simulation technology, aiming at the rehabilitation of ankle joint injury, we imported simplified model of proposed 2-UPS/RR (two identical unconstraint kinematic branches with a universal-prismatic-spherical (UPS) structure and two rotating pair (R)) ankle rehabilitation robot into AnyBody Modeling System. Therefore, a human-machine model was established using the HILL-type muscle model and muscle recruitment criteria. This paper investigated the effects of rehabilitation trajectories on biomechanical response during rehabilitation. Additionally, three main lower limb muscles (soleus, peroneal brevis, and extensor digitorum longus) were examined under different rehabilitation trajectories (plantar dorsiflexion, varus or valgus, and compound movement) in the present study. Based on the biomechanical response of lower limbs, the results showed that different muscles had different sensitivities to the change of rehabilitation trajectories. The correlation coefficient between joint force and plantar dorsiflexion angle reached 0.99 (P < 0.01), indicating that the change of joint force was mainly dominated by plantar dorsiflexion/plantar flexion, but also affected by varus or valgus. Safe rehabilitation training can be achieved by controlling the designed 2-UPS/RR rehabilitation robot. The behavior of muscle force and joint force under different rehabilitation trajectories can meet the needs of rehabilitation and treatment of joint diseases, and provide more reasonable suggestions for early rehabilitation.

Prediction of post-operative adding-on or compensatory lumbar curve correction after anterior vertebral body tethering.

PURPOSE: Anterior Vertebral Body Tethering (AVBT), a fusionless surgical technique based on growth modulation, aims to correct pediatric scoliosis over time. However, medium-term curvature changes of the non-instrumented distal lumbar curve remains difficult to predict. The objective was to biomechanically analyze the level below the LIV to evaluate whether adding-on or compensatory lumbar curve after AVBT can be predicted by intervertebral disc (ID) wedging and force asymmetry. METHODS: 33 retrospective scoliotic cases instrumented with AVBT were used to computationally simulate their surgery and 2-year post-operative growth modulation using a finite element model. The cohort was divided into two subgroups according to the lumbar curvature evolution over 2 years: (1) correction > 10° (C); (2) maintaining ± 10° (M). The lumbar Cobb angle and residual ID wedging angle under LIV were measured. Simulated pressures and moments at the superior endplate of LIV + 1 were post-processed. These parameters were correlated at 2 years postoperatively. FINDINGS: On average, the LIV + 1 simulated moment was 538 Nmm for subgroup C, 155 Nmm for subgroup M with lumbar Cobb angle > 20° and 34 Nmm for angle < 20° whereas the ID angle was 1° for C and 0° for M. INTERPRETATION: On average, a positive moment on the LIV + 1 superior growth plate led to correction of the lumbar curvature, whereas a null moment kept it stable, and a parallel immediate postoperative ID under LIV contributed to its correction or preservation. Nevertheless, the significant interindividual variability suggested that other parameters are involved in the distal non-instrumented curvature evolution. LEVEL OF EVIDENCE: IV.

Case study: The influence of Achilles tendon rupture on knee joint stress during counter-movement jump - Combining musculoskeletal modeling and finite element analysis.

Background: Presently, the current research concerning Achilles tendon rupture repair (ATR) is predominantly centered on the ankle joint, with a paucity of evidence regarding its impact on the knee joint. ATR has the potential to significantly impede athletic performance and increase tibiofemoral contact forces in athletes. The purpose of this study was to prognosticate the distribution of stress within the knee joint during a countermovement jump through the use of a simulation method that amalgamated a musculoskeletal model of a patient who underwent Achilles tendon rupture repair with a finite element model of the knee joint. Methods: A male elite badminton player who had suffered an acute Achilles tendon rupture in his right leg one year prior was selected as our study subject. In order to analyze his biomechanical data, we employed both the OpenSim musculoskeletal model and finite element model to compute various parameters such as joint angles, joint moments, joint contact forces, and the distribution of knee joint stress. Results: During the jumping phase, a significantly lower knee extension angle (p < 0.001), ankle dorsiflexion angle (p = 0.002), peak vertical ground reaction force (p < 0.001), and peak tibiofemoral contact force (p = 0.009) were observed on the injured side than on the uninjured side. During the landing phase, the ankle range of motion (ROM) was significantly lower on the injured side than on the uninjured side (p = 0.009), and higher peak vertical ground reaction forces were observed (p = 0.012). Additionally, it is logical that an injured person will put higher load on the uninjured limb, but the finite element analysis indicated that the stresses on the injured side of medial meniscus and medial cartilage were significantly greater than the uninjured side. Conclusions: An Achilles tendon rupture can limit ankle range of motion and lead to greater joint stress on the affected area during countermovement jumps, especially during the landing phase. This increased joint stress may also transfer more stress to the soft tissues of the medial knee, thereby increasing the risk of knee injury. It is worth noting that this study only involves the average knee flexion angle and load after ATR in one athlete. Caution should be exercised when applying the conclusions, and in the future, more participants should be recruited to establish personalized knee finite element models to validate the results.

Estimation of the biomechanical mammographic deformation of the breast using machine learning models.

BACKGROUND: A typical problem in the registration of MRI and X-ray mammography is the nonlinear deformation applied to the breast during mammography. We have developed a method for virtual deformation of the breast using a biomechanical model automatically constructed from MRI. The virtual deformation is applied in two steps: unloaded state estimation and compression simulation. The finite element method is used to solve the deformation process. However, the extensive computational cost prevents its usage in clinical routine. METHODS: We propose three machine learning models to overcome this problem: an extremely randomized tree (first model), extreme gradient boosting (second model), and deep learning-based bidirectional long short-term memory with an attention layer (third model) to predict the deformation of a biomechanical model. We evaluated our methods with 516 breasts with realistic compression ratios up to 76%. FINDINGS: We first applied one-fold validation, in which the second and third models performed better than the first model. We then applied ten-fold validation. For the unloaded state estimation, the median RMSE for the second and third models is 0.8 mm and 1.2 mm, respectively. For the compression, the median RMSE is 3.4 mm for both models. We evaluated correlations between model accuracy and characteristics of the clinical datasets such as compression ratio, breast volume, and tissue types. INTERPRETATION: Using the proposed models, we achieved accurate results comparable to the finite element model, with a speedup of factor 240 using the extreme gradient boosting model. These proposed models can replace the finite element model simulation, enabling clinically relevant real-time application.

Automated tooth crown design with optimized shape and biomechanics properties.

Despite the large demand for dental restoration each year, the design of crown restorations is mainly performed via manual software operation, which is tedious and subjective. Moreover, the current design process lacks biomechanics optimization, leading to localized stress concentration and reduced working life. To tackle these challenges, we develop a fully automated algorithm for crown restoration based on deformable model fitting and biomechanical optimization. From a library of dental oral scans, a conditional shape model (CSM) is constructed to represent the inter-teeth shape correlation. By matching the CSM to the patient's oral scan, the optimal crown shape is estimated to coincide with the surrounding teeth. Next, the crown is seamlessly integrated into the finish line of preparation via a surface warping step. Finally, porous internal supporting structures of the crown are generated to avoid excessive localized stresses. This algorithm is validated on clinical oral scan data and achieved less than 2 mm mean surface distance as compared to the manual designs of experienced human operators. The mechanical simulation was conducted to prove that the internal supporting structures lead to uniform stress distribution all over the model.

Dynamic 3D Modeling for Human Sperm Motility through the Female Cervical Canal and Uterine Cavity to Predict Sperm Chance of Reaching the Oocyte.

Sperm motility in the female genital tract is a key factor in the natural selection of competent cells that will produce a healthy offspring. We created a dynamic three-dimensional (3D) mechanical model of human sperm cells swimming inside cervical canal and uterine cavity dynamic 3D models, all generated based on experimental studies. Using these simulations, we described the sperm cells' behaviors during swimming inside the 3D tract model as a function of 3D displacement and time. We evaluated normal- and abnormal-morphology sperm cells according to their chances of reaching the oocyte site. As expected, we verified that the number of normal sperm cells that succeeded in reaching the fallopian tube sites is greater than the number of abnormal sperm cells. However, interestingly, after inspecting various abnormal sperm cells, we found out that their scores changed compared to swimming in an infinite medium, as is the case with in vitro fertilization. Thus, the interactions of abnormal sperm cells and the complicated geometry and dynamics of the uterus are significant factors in the filtering of abnormal sperm cells until they reach the oocyte site. Our study provides an advanced tool for sperm analysis and selection criteria for fertility treatments.

Personalisation of Plantarflexor Musculotendon Model Parameters in Children with Cerebral Palsy.

Neuromusculoskeletal models can be used to evaluate aberrant muscle function in cerebral palsy (CP), for example by estimating muscle and joint contact forces during gait. However, to be accurate, models should include representative musculotendon parameters. We aimed to estimate personalised parameters that capture the mechanical behaviour of the plantarflexors in children with CP and typically developing (TD) children. Ankle angle (using motion capture), torque (using a load-cell), and medial gastrocnemius fascicle lengths (using ultrasound) were measured during slow passive ankle dorsiflexion rotation for thirteen children with spastic CP and thirteen TD children. Per subject, the measured rotation was input to a scaled OpenSim model to simulate the torque and fascicle length output. Musculotendon model parameters were personalised by the best match between simulated and experimental torque-angle and fascicle length-angle curves according to a least-squares fit. Personalised tendon slack lengths were significantly longer and optimal fibre lengths significantly shorter in CP than model defaults and than in TD. Personalised tendon compliance was substantially higher in both groups compared to the model default. The presented method to personalise musculotendon parameters will likely yield more accurate simulations of subject-specific muscle mechanics, to help us understand the effects of altered musculotendon properties in CP.

Capturing contact in mitral valve dynamic closure with fluid-structure interaction simulation.

PURPOSE: Realistic fluid-structure interaction (FSI) simulation of the mitral valve opens the way toward planning for surgical repair. In the literature, blood leakage is identified by measuring the flow rate, but detailed information about closure efficiency is missing. We present in this paper an FSI model that improves the detection of blood leakage by building a map of contact. METHODS: Our model is based on the immersed boundary method that captures a map of contact and perfect closure of the mitral valve, without the presence of orifice holes, which often appear with existing methods. We also identified important factors influencing convergence issues. RESULTS: The method is demonstrated in three typical clinical situations: mitral valve with leakage, bulging, and healthy. In addition to the classical ways of evaluating MV closure, such as stress distribution and flow rate, the contact map provides easy detection of leakage with identification of the sources of leakage and a quality assessment of the closure. CONCLUSIONS: Our method significantly improves the quality of the simulation and allows the identification of regurgitation as well as a spatial evaluation of the quality of valve closure. Comparably fast simulation, ability to simulate large deformation, and capturing detailed contact are the main aspects of the study.

Consideration of specific key points improves outcome of decompression treatment in patients with endocrine orbitopathy: pre-/post-OP comparison and biomechanical simulation.

Endocrine orbitopathy is typically treated by resecting orbital walls. This procedure reduces intraorbital pressure by releasing intraorbital tissue, effectively alleviating the symptoms. However, selection of an appropriate surgical plan for treatment of endocrine orbitopathy requires careful consideration because predicting the effects of one-, two-, or three-wall resections on the release of orbital tissues is difficult. Here, based on our experience, we describe two specific orbital sites ('key points') that may significantly improve decompression results. Methodological framework of this work is mainly based on comparative analysis pre- and post-surgery tomographic images as well as image- and physics-based simulation of soft tissue outcome using the finite element modelling of mechanical soft tissue behaviour. Thereby, the optimal set of unknown modelling parameters was obtained iteratively from the minimum difference between model predictions and post-surgery ground truth data. This report presents a pre-/post-surgery study indicating a crucial role of these particular key points in improving the post-surgery outcome of decompression treatment of endocrine orbitopathy which was also supported by 3D biomechanical simulation of alternative two-wall resection plans. In particular, our experimental results show a nearly linear relationship between the resection area and amount of tissue released in the extraorbital space. However, a disproportionately higher volume of orbital outflow could be achieved under consideration of the two special key points. Our study demonstrates the importance of considering natural biomechanical obstacles to improved outcomes in two-wall resection treatment of endocrine orbitopathy. Further investigations of alternative surgery scenarios and post-surgery data are required to generalize the insights of this feasibility study.