Failure locus of the anterior cruciate ligament: 3D finite element analysis

Anterior cruciate ligament (ACL) disruption is a common injury that is detrimental to an athlete's quality of life. Determining the mechanisms that cause ACL injury is important in order to develop proper interventions. A failure locus defined as various combinations of loadings and movements, internal/external rotation of femur and valgus and varus moments at a 25^o knee flexion angle leading to ACL failure was obtained. The results indicated that varus and valgus movements were more dominant to the ACL injury than femoral rotation. Also, Von Mises stress in the lateral tibial cartilage during the valgus ACL injury mechanism was 83% greater than that of the medial cartilage during the varus mechanism of ACL injury. The results of this study could be used to develop training programmes focused on the avoidance of the described combination of movements which may lead to ACL injury.     

3-D finite element modelling of facial soft tissue and preliminary application in orthodontics

Prediction of soft tissue aesthetics is important for achieving an optimal outcome in orthodontic treatment planning. Previously, applicable procedures were mainly restricted to 2-D profile prediction. In this study, a generic 3-D finite element (FE) model of the craniofacial soft and hard tissue was constructed, and individualisation of the generic model based on cone beam CT data and mathematical transformation was investigated. The result indicated that patient-specific 3-D facial FE model including different layers of soft tissue could be obtained through mathematical model transformation. Average deviation between the transformed model and the real reconstructed one was 0.47 ± 0.77 mm and 0.75 ± 0.84 mm in soft and hard tissue, respectively. With boundary condition defined according to treatment plan, such FE model could be used to predict the result of orthodontic treatment on facial soft tissue.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

Overrepresentation of the COL3A1 AA genotype in Polish skiers with anterior cruciate ligament injury

Although various intrinsic and extrinsic risk factors for anterior cruciate ligament (ACL) rupture have been identified, the exact aetiology of the injury is not yet fully understood. Type III collagen is an important factor in the repair of connective tissue, and certain gene polymorphisms may impair the tensile strength. The aim of this study was to examine the association of the COL3A1 rs1800255 polymorphism with ACL rupture in Polish male recreational skiers. A total of 321 male Polish recreational skiers were recruited for this study; 138 had surgically diagnosed primary ACL ruptures (ACL-injured group) and 183 were apparently healthy male skiers (control group – CON) who had no self-reported history of ligament or tendon injury. Both groups had a comparable level of exposure to ACL injury. Genomic DNA was extracted from the oral epithelial cells. All samples were genotyped on a real-time polymerase chain reaction instrument. The genotype distribution in the ACL-injured group was significantly different than in CON (respectively: AA=10.1 vs 2.2%, AG=22.5 vs 36.1, GG=67.4 vs 61.8%; p=0.0087). The AA vs AG+GG genotype of COL3A1 (odds ratio (OR)=5.05; 95% confidence interval (CI), 1.62-15.71, p=0.003) was significantly overrepresented in the ACL-injured group compared with CON. The frequency of the A allele was higher in the ACL-injured group (21.4%) compared with CON (20.2%), but the difference was not statistically significant (p=0.72). This study revealed an association between the COL3A1 rs1800255 polymorphism and ACL ruptures in Polish skiers.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

3D finite element models of shoulder muscles for computing lines of actions and moment arms

Accurate representation of musculoskeletal geometry is needed to characterise the function of shoulder muscles. Previous models of shoulder muscles have represented muscle geometry as a collection of line segments, making it difficult to account for the large attachment areas, muscle–muscle interactions and complex muscle fibre trajectories typical of shoulder muscles. To better represent shoulder muscle geometry, we developed 3D finite element models of the deltoid and rotator cuff muscles and used the models to examine muscle function. Muscle fibre paths within the muscles were approximated, and moment arms were calculated for two motions: thoracohumeral abduction and internal/external rotation. We found that muscle fibre moment arms varied substantially across each muscle. For example, supraspinatus is considered a weak external rotator, but the 3D model of supraspinatus showed that the anterior fibres provide substantial internal rotation while the posterior fibres act as external rotators. Including the effects of large attachment regions and 3D mechanical interactions of muscle fibres constrains muscle motion, generates more realistic muscle paths and allows deeper analysis of shoulder muscle function.     

Effect of stochastic resonance on proprioception and kinesthesia in anterior cruciate ligament reconstructed patients

Low amplitude mechanical noise vibration has been shown to improve somatosensory acuity in various clinical groups with comparable deficiencies through a phenomenon known as Stochastic Resonance (SR). This technology showed promising outcomes in improving somatosensory acuity in other clinical patients (e.g., Parkinson’s disease and osteoarthritis). Some degree of chronic somatosensory deficiency in the knee has been reported following anterior cruciate ligament (ACL) reconstruction surgery. In this study, the effect of the SR phenomenon on improving knee somatosensory acuity (proprioception and kinesthesia) in female ACL reconstructed (ACLR) participants (n = 19) was tested at three months post-surgery, and the results were compared to healthy controls (n = 28). Proprioception was quantified by the measure of joint position sense (JPS) and kinesthesia with the threshold to detection of passive movement (TDPM).The results based on the statistical analysis demonstrated an overall difference between the somatosensory acuity in the ACLR limb compared to healthy controls (p = 0.007). A larger TDPM was observed in the ACLR limb compared to the healthy controls (p = 0.002). However, the JPS between the ACLR and healthy limbs were not statistically significantly different (p = 0.365). SR significantly improved JPS (p = 0.006) while the effect was more pronounced in the ACLR cohort. The effect on the TDPM did not reach statistical significance (p = 0.681) in either group.In conclusion, deficient kinesthesia in the ACLR limb was observed at three months post-surgery. Also, the positive effects of SR on somatosensory acuity in the ACL reconstructed group warrant further investigation into the use of this phenomenon to improve proprioception in ACLR and healthy groups.     

Simulation of 3D tumor cell growth using nonlinear finite element method

We propose a novel parallel computing framework for a nonlinear finite element method (FEM)-based cell model and apply it to simulate avascular tumor growth. We derive computation formulas to simplify the simulation and design the basic algorithms. With the increment of the proliferation generations of tumor cells, the FEM elements may become larger and more distorted. Then, we describe a remesh and refinement processing of the distorted or over large finite elements and the parallel implementation based on Message Passing Interface to improve the accuracy and efficiency of the simulation. We demonstrate the feasibility and effectiveness of the FEM model and the parallelization methods in simulations of early tumor growth.     

An augmented Lagrangian finite element formulation for 3D contact of biphasic tissues

Biphasic contact analysis is essential to obtain a complete understanding of soft tissue biomechanics, and the importance of physiological structure on the joint biomechanics has long been recognised; however, up to date, there are no successful developments of biphasic finite element contact analysis for three-dimensional (3D) geometries of physiological joints. The aim of this study was to develop a finite element formulation for biphasic contact of 3D physiological joints. The augmented Lagrangian method was used to enforce the continuity of contact traction and fluid pressure across the contact interface. The biphasic contact method was implemented in the commercial software COMSOL Multiphysics 4.2^® (COMSOL, Inc., Burlington, MA). The accuracy of the implementation was verified using 3D biphasic contact problems, including indentation with a flat-ended indenter and contact of glenohumeral cartilage layers. The ability of the method to model multibody biphasic contact of physiological joints was proved by a 3D knee model. The 3D biphasic finite element contact method developed in this study can be used to study the biphasic behaviours of the physiological joints.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

Objective: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA).

Materials and methods: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed.

Results: Proper displacement results and strain–stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found.

Conclusion: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

Mechanical performance of endodontic restorations with prefabricated posts: sensitivity analysis of parameters with a 3D finite element model

Many studies have investigated the effect of different parameters of the endodontically restored tooth on its final strength, using in vitro tests and model simulations. However, the differences in the experimental set-up or modelling conditions and the limited number of parameters studied in each case prevent us from obtaining clear conclusions about the relative importance of each parameter. In this study, a validated 3D biomechanical model of the restored tooth was used for an exhaustive sensitivity analysis. The individual influence of 20 different parameters on the mechanical performance of an endodontic restoration with prefabricated posts was studied. The results bring up the remarkable importance of the loading angle on the final restoration strength. Flexural loads are more critical than compressive or tensile loads. Young's modulus of the post and its length and diameter are the most influential parameters for strength, whereas other parameters such as ferrule geometry or core and crown characteristics are less significant.     

On the use of AAA and AcurosXB algorithms for three different stereotactic ablative body radiotherapy (SABR) techniques: Volumetric modulated arc therapy (VMAT), intensity modulated radiation therapy (IMRT) and 3D conformal radiotherapy (3D-CRT)

AimThe purpose of this study was to investigate the dosimetric characteristics of three stereotactic ablative body radiotherapy (SABR) techniques using the anisotropic analytical algorithm (AAA) and Acuros XB algorithm. The SABR techniques include coplanar volumetric modulated arc therapy (C-VMAT), non-coplanar intensity modulated radiation therapy (NC-IMRT) and non-coplanar three-dimensional conformal radiotherapy (NC-3D CRT).BackgroundSABR is a special type of radiotherapy where a high dose of radiation is delivered over a short time. The treatment outcome and accuracy of the dose delivered to cancer patients highly depend on the dose calculation algorithm and treatment technique.Materials and methodsTwelve lung cancer patients underwent 4D CT scanning, and three different treatment plans were generated: C-VMAT, NC-IMRT, NC-3D CRT. Dose calculation was performed using the AAA and Acuros XB algorithm. The dosimetric indices, such as conformity index (CI), homogeneity index, dose fall-off index, doses received by organs at risk and planning target volume, were used to compare the plans. The accuracy of AAA and Acuros XB (AXB) algorithms for the lung was validated against measured dose on a CIRS thorax phantom.ResultsThe CIs for C-VMAT, NC-IMRT and NC-3D CRT were 1.21, 1.28 and 1.38 for the AAA, respectively, and 1.17, 1.26 and 1.36 for the Acuros XB algorithm, respectively. The overall dose computed by AcurosXB algorithm was close to the measured dose when compared to the AAA algorithm. The overall dose computed by the AcurosXB algorithm was close to the measured dose when compared to the AAA algorithm.ConclusionThis study showed that the treatment planning results obtained using the Acuros XB algorithm was better than those using the AAA algorithm in SABR lung radiotherapy.     

Stress distribution of three-unit fixed partial prostheses (conventional and pontic) supported by three or two implants: 3D finite element analysis of ductile materials

In implantology, when financial or biological feasibility limitations appear, it is necessary to use prostheses with geometries that deviate from the conventional, with a pontic in the absence of an intermediate implant. The aim of this study was analyze and understand the general differences in the stresses generated in implants, components and infrastructures according to the configuration of the prosthesis over three or two implants. Thus, this paper analyzes the von Mises equivalent stresses (VMES) of ductile materials on their external surfaces. The experimental groups: Regular Splinted Conventional Group (RCG), which had conventional infrastructures on 3 regular-length Morse taper implants (4x11?mm); Regular Splinted Pontic Group (RPG), which had infrastructures with intermediate pontics on 2 regular-length Morse taper implants (4x11?mm). The simulations of the groups were created with Ansys Workbench 10.0 software. The results revealed that the RPG presented greater areas of possible fragility due to higher stress concentrations, for example, in the cervical area of the union between the implant and component the top platform of the abutment, as well as greater coverage of the stress by the cervical implant threads. The RPG infrastructure was also more affected by stresses in the connection areas between the prostheses and on the occlusal surface. There is an advantage to using prostheses supported by a greater number of implants (RCG) because this decreases the stress in the analyzed structures and consequently improves stress dissipation to the supporting bone, which would preserve the system.     

Phenotyping of hypertensive heart disease and hypertrophic cardiomyopathy using personalized 3D modelling and cardiac cine MRI

Differential diagnosis of hypertensive heart disease (HHD) and hypertrophic cardiomyopathy (HCM) is clinically challenging but important for treatment management. This study aims to phenotype HHD and HCM in 3D + time domain by using a multiparametric motion-corrected personalized modeling algorithm and cardiac magnetic resonance (CMR). 44 CMR data, including 12 healthy, 16 HHD and 16 HCM cases, were examined. Multiple CMR phenotype data consisting of geometric and dynamic variables were extracted globally and regionally from the models over a full cardiac cycle for comparison against healthy models and clinical reports. Statistical classifications were used to identify the distinctive characteristics and disease subtypes with overlapping functional data, providing insights into the challenges for differential diagnosis of both types of disease. While HCM is characterized by localized extreme hypertrophy of the LV, wall thickening/contraction/strain was found to be normal and in sync, though it was occasionally exaggerated at normotrophic/less severely hypertrophic regions during systole to preserve the overall ejection fraction (EF) and systolic functionality. Additionally, we observed that hypertrophy in HHD could also be localized, although at less extreme conditions (i.e. more concentric). While fibrosis occurs mostly in those HCM cases with aortic obstruction, only minority of HHD patients were found affected by fibrosis. We demonstrate that subgroups of HHD (i.e. preserved and reduced EF: HHDpEF & HHDrEF) have different 3D + time CMR characteristics. While HHDpEF has cardiac functions in normal range, dilation and heart failure are indicated in HHDrEF as reflected by low LV wall thickening/contraction/strain and synchrony, as well as much reduced EF.     

RAG1-mediated ubiquitylation of histone H3 is required for chromosomal V(D)J recombination

RAG1 and RAG2 proteins are key components in V(D)J recombination. The core region of RAG1 is capable of catalyzing the recombination reaction; however, the biological function of non-core RAG1 remains largely unknown. Here, we show that in a murine-model carrying the RAG1 ring-finger conserved cysteine residue mutation (C325Y), V(D)J recombination was abrogated at the cleavage step, and this effect was accompanied by decreased mono-ubiquitylation of histone H3. Further analyses suggest that un-ubiquitylated histone H3 restrains RAG1 to the chromatin by interacting with the N-terminal 218 amino acids of RAG1. Our data provide evidence for a model in which ubiquitylation of histone H3 mediated by the ring-finger domain of RAG1 triggers the release of RAG1, thus allowing its transition into the cleavage phase. Collectively, our findings reveal that the non-core region of RAG1 facilitates chromosomal V(D)J recombination in a ubiquitylation-dependent pathway.     

Synthesis and photoluminescence characteristics of (Y,Gd)BO3:RE (RE = Eu3+, Ce3+, Dy3+ and Tb3+) phosphors for blue chip and near‐UV white LEDs

A series of Eu³⁺‐, Ce³⁺‐, Dy³⁺‐ and Tb³⁺‐doped (Y,Gd)BO₃ phosphors was synthesized by a solid‐state diffusion method. X‐Ray diffraction confirmed their hexagonal structure and the scanning electron microscopy results showed crystalline particles. The excitation spectra revealed that (Y,Gd)BO₃ phosphors doped with Eu³⁺, Ce³⁺ , Dy³⁺ and Tb³⁺ are effectively excited with near UV‐light of 395 nm/blue light, 364, 351 and 314 nm, respectively. Photoluminescence spectra of Eu³⁺‐, Ce³⁺‐ and Tb³⁺/Dy³⁺‐doped phosphor showed intense emission of reddish orange, blue and white light, respectively. The phosphor Y_0.60Gd_0.38BO₃:Ce_0.02 showed CIE 1931 color coordinates of (0.158, 0.031) and better color purity compared with commercially available blue BAM:Eu²⁺ phosphor. The phosphor (Y,Gd)BO₃ doped with Eu³⁺, Dy³⁺ and Tb³⁺ showed CIE 1931 color coordinates of (0.667, 0.332), (0.251, 0.299) and (0.333, 0.391) respectively. Significant photoluminescence characteristics of the prepared phosphors indicate that they might serve as potential candidates for blue chip and near‐UV white light‐emitting diode applications. Copyright © 2015 John Wiley & Sons, Ltd.