Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

A positional FEM Formulation for geometrical non-linear analysis of shells

This work presents a fully non-linear finite element formulation for shell analysis comprising linear strain variation along the thickness of the shell and geometrically exact description for curved triangular elements. The developed formulation assumes positions and generalized unconstrained vectors as the variables of the problem, not displacements and finite rotations. The full 3D Saint-Venant-Kirchhoff constitutive relation is adopted and, to avoid locking, the rate of thickness variation enhancement is introduced. As a consequence, the second Piola-Kirchhoff stress tensor and the Green strain measure are employed to derive the specific strain energy potential. Curved triangular elements with cubic approximation are adopted using simple notation. Selected numerical simulations illustrate and confirm the objectivity, accuracy, path independence and applicability of the proposed technique.

Anisotropic Confinement, Electronic Coupling and Strain Induced Effects Detected by Valence-Band Anisotropy in Self-Assembled Quantum Dots

A method to determine the effects of the geometry and lateral ordering on the electronic properties of an array of one-dimensional self-assembled quantum dots is discussed. A model that takes into account the valence-band anisotropic effective masses and strain effects must be used to describe the behavior of the photoluminescence emission, proposed as a clean tool for the characterization of dot anisotropy and/or inter-dot coupling. Under special growth conditions, such as substrate temperature and Arsenic background, 1D chains of In(0.4)Ga(0.6) As quantum dots were grown by molecular beam epitaxy. Grazing-incidence X-ray diffraction measurements directly evidence the strong strain anisotropy due to the formation of quantum dot chains, probed by polarization-resolved low-temperature photoluminescence. The results are in fair good agreement with the proposed model.

Dynamical Casimir effect for a massless scalar field between two concentric spherical shells with mixed boundary conditions

In this work we analyze the dynamical Casimir effect for a massless scalar field confined between two concentric spherical shells considering mixed boundary conditions. We thus generalize a previous result in literature [Phys. Rev. A 78, 032521 (2008)], where the same problem is approached for the field constrained to the Dirichlet-Dirichlet boundary conditions. A general expression for the average number of particle creation is deduced considering an arbitrary law of radial motion of the spherical shells. This expression is then applied to harmonic oscillations of the shells, and the number of particle production is analyzed and compared with the results previously obtained under Dirichlet-Dirichlet boundary conditions.

Dynamical Casimir effect for a massless scalar field between two concentric spherical shells

In this work we consider the dynamical Casimir effect for a massless scalar field-under Dirichlet boundary conditions-between two concentric spherical shells. We obtain a general expression for the average number of particle creation, for an arbitrary law of radial motion of the spherical shells, using two distinct methods: by computing the density operator of the system and by calculating the Bogoliubov coefficients. We apply our general expression to breathing modes: when only one of the shells oscillates and when both shells oscillate in or out of phase. Since our results were obtained in the framework of the perturbation theory, under resonant breathing modes they are restricted to a short-time approximation. We also analyze the number of particle production and compare it with the results for the case of plane geometry.

Ultra-high-energy cosmic rays from centaurus A: jet interaction with gaseous shells

Ultra-high-energy cosmic rays (UHECRs), with energies above similar to 6 x 10(19) eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model, many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 2: shells

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.

Time flies with music whatever its emotional valence

The present study used a temporal bisection task to investigate whether music affects time estimation differently from a matched auditory neutral stimulus, and whether the emotional valence of the musical stimuli (i.e., sad vs. happy music) modulates this effect. The results showed that, compared to sine wave control music, music presented in a major (happy) or a minor (sad) key shifted the bisection function toward the right, thus increasing the bisection point value (point of subjective equality). This indicates that the duration of a melody is judged shorter than that of a non-melodic control stimulus, thus confirming that ""time flies"" when we listen to music. Nevertheless, sensitivity to time was similar for all the auditory stimuli. Furthermore, the temporal bisection functions did not differ as a function of musical mode. (C) 2010 Elsevier B.V. All rights reserved.

On an Electrode Modified by a Supramolecular Ruthenium Mixed Valence (Ru(II)/Ru(III)) Diphosphine-Porphyrin Assembly

Three novel polymetallic ruthenium (III) meso-tetra(4-pyridyl)porphyrins containing peripheral ""RuCl(3)(dppb)"" moieties have been prepared and characterized. The X-ray structure of the tetraruthenated {NiTPyP[RuCl(3)(dppb)](4)} porphyrin complex crystallizes in the triclinic space group FT. This structure is discussed and compared with the crystal data for the mer-[RuCl(3)(dppb)(py)]. The {TPyP[RuCl(3)(dppb)](4)} and {CoTPyP[RuCl(3)(dppb)](4)} porphyrins were used to obtain electrogenerated films on ITO and glass carbon electrode surfaces, respectively. Such tetraruthenated porphyrins form films of a mixed-valence species {TPyP[Ru(dppb)](4)(mu Cl(3))(2)}(2)(4n2+) and {CoTPyP[Ru(dppb)](4)(mu Cl(3))(2)}(2n)(4n2+) on the electrode surface. The modified electrode with {CoTPyP[RuCl(3)(dppb)](4)} is very stable and can be used to detect organic substrates such as catechol.

Study on the observation of Eu(2+) and Eu(3+) valence states in low silica calcium aluminosilicate glasses

The optical, magnetic and structural properties of Eu doped low silica calcium aluminosilicate glasses were investigated. The optical absorption coefficient presented two bands at 39 246 and 29 416 cm(-1), which were assigned respectively to the 4f(7) ((8)S(7/2)) -> 4f(6) (4F(J)) 5d (T(2g)), and 4f(7) ((8)S(7/2)) -> 4f(6) (4F(J)) 5d (E(g)) transitions of Eu(2+). The fluorescence measured at 300 K on a sample doped with 0.5 wt% of Eu(2)O(3) exhibited a broad band centered at 17 350 cm(-1), which is attributed to the 4f(6)5d -> 4f(7) transition of Eu(2+), whereas the additional peaks are due to the (5)D(0) -> (7)F(J) (J = 1, 2, 4) transitions of Eu(3+). From magnetization and XANES data it was possible to evaluate the fractions of Eu(2+) and Eu(3+) for the sample doped with 0.5 and 5.0 wt% of Eu(2)O(3), the values of which were approximately 30 and 70%, respectively.

Linear dimensional changes in plaster die models using different elastomeric materials

Dental impression is an important step in the preparation of prostheses since it provides the reproduction of anatomic and surface details of teeth and adjacent structures. The objective of this study was to evaluate the linear dimensional alterations in gypsum dies obtained with different elastomeric materials, using a resin coping impression technique with individual shells. A master cast made of stainless steel with fixed prosthesis characteristics with two prepared abutment teeth was used to obtain the impressions. References points (A, B, C, D, E and F) were recorded on the occlusal and buccal surfaces of abutments to register the distances. The impressions were obtained using the following materials: polyether, mercaptan-polysulfide, addition silicone, and condensation silicone. The transfer impressions were made with custom trays and an irreversible hydrocolloid material and were poured with type IV gypsum. The distances between identified points in gypsum dies were measured using an optical microscope and the results were statistically analyzed by ANOVA (p < 0.05) and Tukey's test. The mean of the distances were registered as follows: addition silicone (AB = 13.6 µm, CD=15.0 µm, EF = 14.6 µm, GH=15.2 µm), mercaptan-polysulfide (AB = 36.0 µm, CD = 36.0 µm, EF = 39.6 µm, GH = 40.6 µm), polyether (AB = 35.2 µm, CD = 35.6 µm, EF = 39.4 µm, GH = 41.4 µm) and condensation silicone (AB = 69.2 µm, CD = 71.0 µm, EF = 80.6 µm, GH = 81.2 µm). All of the measurements found in gypsum dies were compared to those of a master cast. The results demonstrated that the addition silicone provides the best stability of the compounds tested, followed by polyether, polysulfide and condensation silicone. No statistical differences were obtained between polyether and mercaptan-polysulfide materials.

Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

Obtenção e caracterização de α-quitina e quitosanas de cascas de Macrobrachium rosembergii

The shells of Macrobrachium rosenbergii were submitted to deproteinization (Dp) and demineralization (Dm) aiming the extraction of α-chitin. The different parts of the shells were processed independently by carrying out sequence 1 (Dp/Dm) and sequence 2 (Dm/Dp). Both sequences allowed the extraction of chitins with low contents of calcium and magnesium, regardless of the part being processed. The sequence 1 lead to higher extraction yields while sequence 2 resulted in lower contents of inorganic compounds. Extensively deacetylated chitosans (GA<10%) of medium molecular weight (0,9 x 10(5) < Mv < 2 x 10(5) g/mol) resulted from the deacetylation of chitin.

Two dimensional analysis of inflatable structures by the positional FEM

This paper presents a, simple two dimensional frame formulation to deal with structures undergoing large motions due to dynamic actions including very thin inflatable structures, balloons. The proposed methodology is based on the minimum potential energy theorem written regarding nodal positions. Velocity, acceleration and strain are achieved directly from positions, not. displacements, characterizing the novelty of the proposed technique. A non-dimensional space is created and the deformation function (change of configuration) is written following two independent mappings from which the strain energy function is written. The classical New-mark equations are used to integrate time. Dumping and non-conservative forces are introduced into the mechanical system by a rheonomic energy function. The final formulation has the advantage of being simple and easy to teach, when compared to classical Counterparts. The behavior of a bench-mark problem (spin-up maneuver) is solved to prove the formulation regarding high circumferential speed applications. Other examples are dedicated to inflatable and very thin structures, in order to test the formulation for further analysis of three dimensional balloons.