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.

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e. g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3577999]

Intercellular transport of epidermis-expressed MADS domain transcription factors and their effect on plant morphology and floral transition

P>During the lifetime of an angiosperm plant various important processes such as floral transition, specification of floral organ identity and floral determinacy, are controlled by members of the MADS domain transcription factor family. To investigate the possible non-cell-autonomous function of MADS domain proteins, we expressed GFP-tagged clones of AGAMOUS (AG), APETALA3 (AP3), PISTILLATA (PI) and SEPALLATA3 (SEP3) under the control of the MERISTEMLAYER1 promoter in Arabidopsis thaliana plants. Morphological analyses revealed that epidermal overexpression was sufficient for homeotic changes in floral organs, but that it did not result in early flowering or terminal flower phenotypes that are associated with constitutive overexpression of these proteins. Localisations of the tagged proteins in these plants were analysed with confocal laser scanning microscopy in leaf tissue, inflorescence meristems and floral meristems. We demonstrated that only AG is able to move via secondary plasmodesmata from the epidermal cell layer to the subepidermal cell layer in the floral meristem and to a lesser extent in the inflorescence meristem. To study the homeotic effects in more detail, the capacity of trafficking AG to complement the ag mutant phenotype was compared with the capacity of the non-inwards-moving AP3 protein to complement the ap3 mutant phenotype. While epidermal expression of AG gave full complementation, AP3 appeared not to be able to drive all homeotic functions from the epidermis, perhaps reflecting the difference in mobility of these proteins.

Multifractal regime transition in a modified minority game model

The search for more realistic modeling of financial time series reveals several stylized facts of real markets. In this work we focus on the multifractal properties found in price and index signals. Although the usual minority game (MG) models do not exhibit multifractality, we study here one of its variants that does. We show that the nonsynchronous MG models in the nonergodic phase is multifractal and in this sense, together with other stylized facts, constitute a better modeling tool. Using the structure function (SF) approach we detected the stationary and the scaling range of the time series generated by the MG model and, from the linear (non-linear) behavior of the SF we identified the fractal (multifractal) regimes. Finally, using the wavelet transform modulus maxima (WTMM) technique we obtained its multifractal spectrum width for different dynamical regimes. (C) 2009 Elsevier Ltd. All rights reserved.

Non-linear boundary element formulation applied to contact analysis using tangent operator

This work presents a non-linear boundary element formulation applied to analysis of contact problems. The boundary element method (BEM) is known as a robust and accurate numerical technique to handle this type of problem, because the contact among the solids occurs along their boundaries. The proposed non-linear formulation is based on the use of singular or hyper-singular integral equations by BEM, for multi-region contact. When the contact occurs between crack surfaces, the formulation adopted is the dual version of BEM, in which singular and hyper-singular integral equations are defined along the opposite sides of the contact boundaries. The structural non-linear behaviour on the contact is considered using Coulomb`s friction law. The non-linear formulation is based on the tangent operator in which one uses the derivate of the set of algebraic equations to construct the corrections for the non-linear process. This implicit formulation has shown accurate as the classical approach, however, it is faster to compute the solution. Examples of simple and multi-region contact problems are shown to illustrate the applicability of the proposed scheme. (C) 2011 Elsevier Ltd. All rights reserved.

A simple method for non-linear analysis of steel fiber reinforced concrete

This paper proposes a physical non-linear formulation to deal with steel fiber reinforced concrete by the finite element method. The proposed formulation allows the consideration of short or long fibers placed arbitrarily inside a continuum domain (matrix). The most important feature of the formulation is that no additional degree of freedom is introduced in the pre-existent finite element numerical system to consider any distribution or quantity of fiber inclusions. In other words, the size of the system of equations used to solve a non-reinforced medium is the same as the one used to solve the reinforced counterpart. Another important characteristic of the formulation is the reduced work required by the user to introduce reinforcements, avoiding ""rebar"" elements, node by node geometrical definitions or even complex mesh generation. Bounded connection between long fibers and continuum is considered, for short fibers a simplified approach is proposed to consider splitting. Non-associative plasticity is adopted for the continuum and one dimensional plasticity is adopted to model fibers. Examples are presented in order to show the capabilities of the formulation.

Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.

A solid-like FEM for geometrically non-linear 3D frames

This study presents a solid-like finite element formulation to solve geometric non-linear three-dimensional inhomogeneous frames. To achieve the desired representation, unconstrained vectors are used instead of the classic rigid director triad; as a consequence, the resulting formulation does not use finite rotation schemes. High order curved elements with any cross section are developed using a full three-dimensional constitutive elastic relation. Warping and variable thickness strain modes are introduced to avoid locking. The warping mode is solved numerically in FEM pre-processing computational code, which is coupled to the main program. The extra calculations are relatively small when the number of finite elements. with the same cross section, increases. The warping mode is based on a 2D free torsion (Saint-Venant) problem that considers inhomogeneous material. A scheme that automatically generates shape functions and its derivatives allow the use of any degree of approximation for the developed frame element. General examples are solved to check the objectivity, path independence, locking free behavior, generality and accuracy of the proposed formulation. (C) 2009 Elsevier B.V. All rights reserved.

A FEM procedure based on positions and unconstrained vectors applied to non-linear dynamic of 3D frames

This study presents an alternative three-dimensional geometric non-linear frame formulation based on generalized unconstrained vector and positions to solve structures and mechanisms subjected to dynamic loading. The formulation is classified as total Lagrangian with exact kinematics description. The resulting element presents warping and non-constant transverse strain modes, which guarantees locking-free behavior for the adopted three-dimensional constitutive relation, Saint-Venant-Kirchhoff, for instance. The application of generalized vectors is an alternative to the use of finite rotations and rigid triad`s formulae. Spherical and revolute joints are considered and selected dynamic and static examples are presented to demonstrate the accuracy and generality of the proposed technique. (C) 2010 Elsevier B.V. All rights reserved.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Electrochemical study of modified non-functional bis-silane layers on Al alloy 2024-T3

In the last few years great efforts have been made in order to find and to develop environmentally friendly substitutes for Cr6+ pre-treatments applied on aluminium alloys used in the aircraft industry. Among the potential substitutes, silane layers have attracted considerable interest from researchers and from the industry. The present work investigates the anti-corrosion behaviour of (bis-1, 2-(triethoxysilyl) ethane (BTSE)) silane layers modified with Ce ions and/or silica nanoparticles applied on Al alloy 2024-T3 substrates. The corrosion behaviour was investigated in 0.1 M NaCl solution via d.c. polarization and electrochemical impedance spectroscopy (EIS). Contact angle measurements and XPS were used to assess information on the chemistry of the silane pre-treated surfaces. The results have shown that the introduction of additives improves the corrosion protection properties of the silane layer. (c) 2008 Elsevier Ltd. All rights reserved.

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes. (C) 2008 Elsevier Ltd. All rights reserved.

Inorganic elements in the fat bodies of Diatraea saccharalis (Lepidoptera: Crambidae) larvae parasitized by Cotesia flavipes (Hymenoptera: Braconidae)

Koinobiont parasitoids use several strategies to regulate the host`s physiological processes during parasitism. Although many aspects of host-parasitoid interactions have been explored, studies that attempted to assess the effects of parasitism on the availability of inorganic elements in the host are virtually nonexistent. Therefore, we aimed to evaluate the effects of parasitism on the concentrations of inorganic elements in the fat bodies of larvae of Diatraea saccharalis (Lepidoptera: Crambidae) during the development of the parasitoid Cotesia flavipes (Hymenoptera: Braconidae), by using total reflection X-ray fluorescence (TXRF). TXRF analysis allowed comparisons of the changes in the availability of the elements P. S. K, Ca, Cr, Fe, Ni, Cu, and Zn in the fat body tissues of D. saccharalis larvae parasitized by C. flavipes. Overall, the concentration of inorganic elements was higher early in parasitoid development (1 and 3 days after parasitism) compared to non-parasitized larvae, but much lower towards the end of parasitoid development (7 and 9 days after parasitism). Ca, K, and S were reduced after the fifth day of parasitism, which affected the total abundance of inorganic elements observed in the fat bodies of the parasitized hosts. The regulatory mechanisms or pathological effects related to the observed variation of the host inorganic elements induced by the parasitoid remain unknown, but there might be a strategy to make these elements available to the parasitoid larvae at the end of their development, when higher metabolic activity of the host fat body is required to sustain parasitoid growth. The observed variation of the host`s inorganic elements could also be related to the known effects of parasitism on the host`s immune response. (C) 2010 Elsevier Inc. All rights reserved.

Sudden onset of log-periodicity and superdiffusion in non-Markovian random walks with amnestically induced persistence: exact results

Random walks can undergo transitions from normal diffusion to anomalous diffusion as some relevant parameter varies, for instance the L,vy index in L,vy flights. Here we derive the Fokker-Planck equation for a two-parameter family of non-Markovian random walks with amnestically induced persistence. We investigate two distinct transitions: one order parameter quantifies log-periodicity and discrete scale invariance in the first moment of the propagator, whereas the second order parameter, known as the Hurst exponent, describes the growth of the second moment. We report numerical and analytical results for six critical exponents, which together completely characterize the properties of the transitions. We find that the critical exponents related to the diffusion-superdiffusion transition are identical in the positive feedback and negative feedback branches of the critical line, even though the former leads to classical superdiffusion whereas the latter gives rise to log-periodic superdiffusion.

The attentional demands of preferred and non-preferred gait patterns

The purpose of this study was to determine the attentional demands of natural and imposed gait, as well as the attentional costs of transitions between the walking and running co-ordination patterns. Seven healthy young men and four healthy young women undertook an auditory probe reaction time task concurrently with self-selected gait (Experiment 1) and imposed walking and running (Experiment 2) at different speeds on a motor-driven treadmill. In Experiment 1, where participants were free to choose their own movement pattern to match the speed of travel of the treadmill, normal gait control was shown to have a significant attentional cost, and hence not be automatic in the classical sense. However, this attentional cost did not differ between the two gait modes or at the transition point. In Experiment 2, where participants were required to maintain specific gait modes regardless of the treadmill speed, the maintenance of walking at speeds normally associated with running was found to have an attentional cost whereas this was not the case for running at normal walking speeds. Collectively the findings support a model of gait control in which the normal switching between gait modes is determined with minimal attention demand and in which it is possible to sustain non-preferred gait modes although, in the case of walking, only at a significant attentional/cognitive cost. © 2002 Elsevier Science B.V. All rights reserved.