A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.

Crystallization conditions and controls on trace element residence in the main minerals from the Pedra Branca Syenite, Brazil: An electron microprobe and LA-ICPMS study

Major and trace-element microanalyses of the main minerals from the 610 Ma Pedra Branca Syenite, southeast Brazil, allow inferences on intensive parameters of magmatic crystallization and on the partition of trace-elements among these minerals, with important implications for the petrogenetic evolution of the pluton. Two main syenite types make up the pluton, a quartz-free syenite with tabular alkali feldspar (laminated silica-saturated syenite, LSS, with Na-rich augite + phlogopite + hematite + magnetite + titanite + apatite) and a quartz-bearing syenite (laminated silica-oversaturated syenite, LSO, with scarce corroded plagioclase plus diopside + biotite +/- hornblende + ilmenite magnetite +/- titanite + apatite). Both types share a remarkable enrichment in incompatible elements as K, Ba, Sr, P and LREE. Apatite saturation temperatures of similar to 1060-1090 degrees C are the best estimates of liquidus, whereas the pressure of emplacement, based on Al-in-hornblende barometry, is estimated as 3.3 to 4.8 khan Although both units crystallized under oxidizing conditions, oxygen fugacity was probably higher in LSS, as shown by higher mg# of the mafic minerals and higher hematite contents in Hem-Ilm(ss). In contrast with the Ca-bearing alkali-feldspar from LSO, which hosts most of the whole-rock Sr and Pb, virtually Ca-free alkali-feldspar from LSS hosts similar to 50% of whole-rock Sr and similar to 80% of Pb, the remainder of these elements being shared by apatite, pyroxene and titanite. This contrast reflects a strong crystal-chemical control, whereby a higher proportion of an element with similar ratio and charge (Ca2+) enhances the residence of Sr and Pb in the M-site of alkali feldspar. The more alkaline character of the LSS magma is inferred to have inhibited zircon saturation; Zr + Hf remained in solution until late in the crystallization, and were mostly accommodated in the structure of Ca-Na pyroxene and titanite, which are one order of magnitude richer in these elements compared to the same minerals in LSO, where most of Zr and Hf are inferred to reside in zircon. The REE, Th and U reside mostly in titanite and apatite; D(REE)Tit/Ap raises steadily from 1 to 6 from La to Tb then remains constant up to Lu in the LSO sample; these values are about half as much in the LSS sample, where lower contents of incompatible elements in titanite are attributed to its greater modal abundance and earlier crystallization. (C) 2012 Elsevier B.V. All rights reserved.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s-1 for 9.1 mu s. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.130.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

CNO and F abundances in the globular cluster M22 (NGC 6656)

Context. Recent studies have confirmed the long standing suspicion that M 22 shares a metallicity spread and complex chemical enrichment history similar to that observed in omega Cen. M 22 is among the most massive Galactic globular clusters and its color-magnitude diagram and chemical abundances reveal the existence of sub-populations. Aims. To further constrain the chemical diversity of M 22, necessary to interpret its nucleosynthetic history, we seek to measure relative abundance ratios of key elements (carbon, nitrogen, oxygen, and fluorine) best studied, or only available, using high-resolution spectra at infrared wavelengths. Methods. High-resolution (R = 50 000) and high S/N infrared spectra were acquired of nine red giant stars with Phoenix at the Gemini-South telescope. Chemical abundances were calculated through a standard 1D local thermodynamic equilibrium analysis using Kurucz model atmospheres. Results. We derive [Fe/H] = -1.87 to -1.44, confirming at infrared wavelengths that M 22 does present a [Fe/H] spread. We also find large C and N abundance spreads, which confirm previous results in the literature but based on a smaller sample. Our results show a spread in A(C+N+O) of similar to 0.7 dex. Similar to mono-metallic globular clusters, M 22 presents a strong [Na/Fe]-[O/Fe] anticorrelation as derived from Na and CO lines in the K band. For the first time we recover F abundances in M 22 and find that it exhibits a 0.6 dex variation. We find tentative evidence for a flatter A(F)-A(O) relation compared to higher metallicity globular clusters. Conclusions. Our study confirms and expands upon the chemical diversity seen in this complex stellar system. All elements studied to date show large abundance spreads which require contributions from both massive and low mass stars.

Elemental abundances and classification of carbon-enhanced metal-poor stars

We present a detailed study of carbon-enhanced metal-poor (CEMP) stars, based on high-resolution spectroscopic observations of a sample of 18 stars. The stellar spectra for this sample were obtained at the 4.2 m William Herschel Telescope in 2001 and 2002, using the Utrecht Echelle Spectrograph, at a resolving power R similar to 52 000 and S/N similar to 40, covering the wavelength range lambda lambda 3700-5700 angstrom. The atmospheric parameters determined for this sample indicate temperatures ranging from 4750 K to 7100 K, log g from 1.5 to 4.3, and metallicities -3.0 <= [Fe/H]<=-1.7. Elemental abundances for C, Na, Mg, Sc, Ti, Cr, Cu, Zn, Sr, Y, Zr, Ba, La, Ce, Nd, Sm, Eu, Gd, Dy are determined. Abundances for an additional 109 stars were taken from the literature and combined with the data of our sample. The literature sample reveals a lack of reliable abundance estimates for species that might be associated with the r-process elements for about 67% of CEMP stars, preventing a complete understanding of this class of stars, since [Ba/Eu] ratios are used to classify them. Although eight stars in our observed sample are also found in the literature sample, Eu abundances or limits are determined for four of these stars for the first time. From the observed correlations between C, Ba, and Eu, we argue that the CEMP-r/s class has the same astronomical origin as CEMP-s stars, highlighting the need for a more complete understanding of Eu production.

Inhibition of cannabinoid CB1 receptor upregulates Slc2a4 expression via nuclear factor-kappa B and sterol regulatory element-binding protein-1 in adipocytes

Evidences have suggested that the endocannabinoid system is overactive in obesity, resulting in enhanced endocannabinoid levels in both circulation and visceral adipose tissue. The blockade of cannabinoid receptor type 1 (CB1) has been proposed for the treatment of obesity. Besides loss of body weight, CB1 antagonism improves insulin sensitivity, in which the glucose transporter type 4 (GLUT4) plays a key role. The aim of this study was to investigate the modulation of GLUT4-encoded gene (Slc2a4 gene) expression by CB1 receptor. For this, 3T3-L1 adipocytes were incubated in the presence of a highly selective CB1 receptor agonist (1 mu M arachidonyl-2'-chloroethylamide) and/or a CB1 receptor antagonist/inverse agonist (0.1, 0.5, or 1 mu M AM251, 1-(2,4-dichlorophenyl)-5-(4-iodophenyl)-4-methyl-N-1-piperidinyl-1H-pyrazole-3-carboxamide). After acute (2 and 4 h) and chronic (24 h) treatments, cells were harvested to evaluate: i) Slc2a4, Cnr1 (CB1 receptor-encoded gene), and Srebf1 type a (SREBP-1a type-encoded gene) mRNAs (real-time PCR); ii) GLUT4 protein (western blotting); and iii) binding activity of nuclear factor (NF)-kappa B and sterol regulatory element-binding protein (SREBP)-1 specifically in the promoter of Slc2a4 gene (electrophoretic mobility shift assay). Results revealed that both acute and chronic CB1 receptor antagonism greatly increased (similar to 2.5-fold) Slc2a4 mRNA and protein content. Additionally, CB1-induced upregulation of Slc2a4 was accompanied by decreased binding activity of NF-kappa B at 2 and 24 h, and by increased binding activity of the SREBP-1 at 24 h. In conclusion, these findings reveal that the blockade of CB1 receptor markedly increases Slc2a4/GLUT4 expression in adipocytes, a feature that involves NF-kappa B and SREBP-1 transcriptional regulation. Journal of Molecular Endocrinology (2012) 49, 97-106

FINITE ELEMENT ANALYSIS OF MAXILLARY BONE STRESS CAUSED BY ARAMANY CLASS IV OBTURATOR PROSTHESES

Statement of problem. The retention of an Aramany Class IV removable partial dental prosthesis can be compromised by a lack of support. The biomechanics of this obturator prosthesis result in an unusual stress distribution on the residual maxillary bone. Purpose. This study evaluated the biomechanics of an Aramany Class IV obturator prosthesis with finite element analysis and a digital 3-dimensional (3-D) model developed from a computed tomography scan; bone stress was evaluated according to the load placed on the prosthesis. Material and methods. A 3-D model of an Aramany Class IV maxillary resection and prosthesis was constructed. This model was used to develop a finite element mesh. A 120 N load was applied to the occlusal and incisal platforms corresponding to the prosthetic teeth. Qualitative analysis was based on the scale of maximum principal stress; values obtained through quantitative analysis were expressed in MPa. Results. Under posterior load, tensile and compressive stresses were observed; the tensile stress was greater than the compressive stress, regardless of the bone region, and the greatest compressive stress was observed on the anterior palate near the midline. Under an anterior load, tensile stress was observed in all of the evaluated bone regions; the tensile stress was greater than the compressive stress, regardless of the bone region. Conclusions. The Aramany Class IV obturator prosthesis tended to rotate toward the surgical resection when subjected to posterior or anterior loads. The amount of tensile and compressive stress caused by the Aramany Class IV obturator prosthesis did not exceed the physiological limits of the maxillary bone tissue. (J Prosthet Dent 2012;107:336-342)

Probabilistic crack growth analyses using a boundary element model: Applications in linear elastic fracture and fatigue problems

This paper addresses the numerical solution of random crack propagation problems using the coupling boundary element method (BEM) and reliability algorithms. Crack propagation phenomenon is efficiently modelled using BEM, due to its mesh reduction features. The BEM model is based on the dual BEM formulation, in which singular and hyper-singular integral equations are adopted to construct the system of algebraic equations. Two reliability algorithms are coupled with BEM model. The first is the well known response surface method, in which local, adaptive polynomial approximations of the mechanical response are constructed in search of the design point. Different experiment designs and adaptive schemes are considered. The alternative approach direct coupling, in which the limit state function remains implicit and its gradients are calculated directly from the numerical mechanical response, is also considered. The performance of both coupling methods is compared in application to some crack propagation problems. The investigation shows that direct coupling scheme converged for all problems studied, irrespective of the problem nonlinearity. The computational cost of direct coupling has shown to be a fraction of the cost of response surface solutions, regardless of experiment design or adaptive scheme considered. (C) 2012 Elsevier Ltd. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Finite element analysis of stress distribution in anchor teeth in surgically assisted rapid palatal expansion

The treatment of a transverse maxillary deficiency in skeletally mature individuals should include surgically assisted rapid palatal expansion. This study evaluated the distribution of stresses that affect the expander's anchor teeth using finite element analysis when the osteotomy is varied. Five virtual models were built and the surgically assisted rapid palatal expansion was simulated. Results showed tension on the lingual face of the teeth and alveolar bone, and compression on the buccal side of the alveolar bone. The subtotal Le Fort I osteotomy combined with intermaxillary suture osteotomy seemed to reduce the dissipation of tensions. Therefore, subtotal Le Fort I osteotomy without a step in the zygomaticomaxillary buttress, combined with intermaxillary suture osteotomy and pterygomaxillary disjunction may be the osteotomy of choice to reduce tensions on anchor teeth, which tend to move mesiobuccally (premolar) and distobuccally (molar)

Rainwater major and trace element contents in Southeastern Brazil: an assessment of a sugar cane region in dry and wet period

The purpose of this study was to assess the composition of the rainwater in Araraquara City, Brazil, a region strongly influenced by pre-harvest burning of sugar cane crops. Chemical and mineralogical variables were measured in rainwater collected during the harvest, dry period of 2009 and the non-harvest, wet period of 2010. Ca2+ and NH4+ were responsible for 55% of cations and NO3- for 45% of anions in rainwater. Al and Fe along with K were the most abundant among trace elements in both soluble and insoluble fractions. High volume weighted mean concentration (VWM) for most of the analyzed species were observed in the harvest, dry period, mainly due to agricultural activities and meteorological conditions. The chemistry of the Araraquara rainwater and principal component analysis (PCA) quantification clearly indicate the concurrence of a diversity of sources from natural to anthropogenic especially related to agricultural activities.

Morse taper implants at different bone levels: a finite element analysis of stress distribution

AIM: To explore the biomechanical effects of the different implantation bone levels of Morse taper implants, employing a finite element analysis (FEA). METHODS: Dental implants (TitamaxCM) with 4x13 mm and 4x11 mm, and their respective abutments with 3.5 mm height, simulating a screwed premolar metal-ceramic crown, had their design performed using the software AnsysWorkbench 10.0. They were positioned in bone blocks, covered by 2.5 mm thickness of mucosa. The cortical bone was designed with 1.5 mm thickness and the trabecular bone completed the bone block. Four groups were formed: group 11CBL (11 mm implant length on cortical bone level), group 11TBL (11 mm implant length on trabecular bone level), group 13CBL (13mm implant length on cortical bone level) and group 13TBL (13 mm implant length on trabecular bone level). Oblique 200 N loads were applied. Von Mises equivalent stresses in cortical and trabecular bones were evaluated with the same design program. RESULTS: The results were shown qualitatively and quantitatively by standard scales for each type of bone. By the results obtained, it can be suggested that positioning the implant completely in trabecular bone brings harm with respect to the generated stresses. Its implantation in the cortical bone has advantages with respect to better anchoring and locking, reflecting a better dissipation of the stresses along the implant/bone interfaces. In addition, the search for anchoring the implant in its apical region in cortical bone is of great value to improve stabilization and consequently better stress distribution. CONCLUSIONS: The implant position slightly below the bone in relation to the bone crest brings advantages as the best long-term predictability with respect to the expected neck bone loss.