GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 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.

Buckling and post-buckling of extensible rods revisited: A multiple-scale solution

An exact non-linear formulation of the equilibrium of elastic prismatic rods subjected to compression and planar bending is presented, electing as primary displacement variable the cross-section rotations and taking into account the axis extensibility. Such a formulation proves to be sufficiently general to encompass any boundary condition. The evaluation of critical loads for the five classical Euler buckling cases is pursued, allowing for the assessment of the axis extensibility effect. From the quantitative viewpoint, it is seen that such an influence is negligible for very slender bars, but it dramatically increases as the slenderness ratio decreases. From the qualitative viewpoint, its effect is that there are not infinite critical loads, as foreseen by the classical inextensible theory. The method of multiple (spatial) scales is used to survey the post-buckling regime for the five classical Euler buckling cases, with remarkable success, since very small deviations were observed with respect to results obtained via numerical integration of the exact equation of equilibrium, even when loads much higher than the critical ones were considered. Although known beforehand that such classical Euler buckling cases are imperfection insensitive, the effect of load offsets were also looked at, thus showing that the formulation is sufficiently general to accommodate this sort of analysis. (c) 2008 Elsevier Ltd. All rights reserved.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

The natural force density method for the shape finding of taut structures

This work presents, with the aid of the natural approach, an extension of the force density method for the initial shape finding of cable and membrane structures, which leads to the solution of a system of linear equations. This method, here called the natural force density method, preserves the linearity which characterizes the original force density method. At the same time, it overcomes the difficulties that the original procedure presents to cope with irregular triangular finite element meshes. Furthermore, if this method is applied iteratively in the lines prescribed herewith, it leads to a viable initial configuration with a uniform, isotropic plane Cauchy stress state. This means that a minimal surface for the membrane can be achieved through a succession of equilibrated configurations. Several numerical examples illustrate the simplicity and robustness of the method. (C) 2008 Elsevier B.V. All rights reserved.

A contribution to the exact modal solution of in-plane beam structures

The exact vibration modes and natural frequencies of planar structures and mechanisms, comprised Euler-Bernoulli beams, are obtained by solving a transcendental. nonlinear, eigenvalue problem stated by the dynamic stiffness matrix (DSM). To solve this kind of problem, the most employed technique is the Wittrick-Williams algorithm, developed in the early seventies. By formulating a new type of eigenvalue problem, which preserves the internal degrees-of-freedom for all members in the model, the present study offers an alternative to the use of this algorithm. The new proposed eigenvalue problem presents no poles, so the roots of the problem can be found by any suitable iterative numerical method. By avoiding a standard formulation for the DSM, the local mode shapes are directly calculated and any extension to the beam theory can be easily incorporated. It is shown that the method here adopted leads to exact solutions, as confirmed by various examples. Extensions of the formulation are also given, where rotary inertia, end release, skewed edges and rigid offsets are all included. (C) 2008 Elsevier Ltd. All rights reserved.

Evidences of the evolution from solid solution to surface segregation in Ni-doped SnO(2) nanoparticles using Raman spectroscopy

Ni-doped SnO(2) nanoparticles, promising for gas-sensing applications, have been synthesized by a polymer precursor method. X-ray diffraction (XRD) and transmission electron microscopy (TEM) data analyses indicate the exclusive formation of nanosized particles with rutile-type phase (tetragonal SnO(2)) for Ni contents below 10 mol%. The mean crystallite size shows a progressive reduction with the Ni content. Room-temperature Raman spectra of Ni-doped SnO(2) nanoparticles show the presence of Raman active modes and modes activated by size effects. From the evolution of the A(1g) mode with the Ni content, a solubility limit at similar to 2 mol% was estimated. Below that content, Raman results are consistent with the occurrence of solid solution (ss) and surface segregation (seg.) of Ni ions. Above similar to 2 mol% Ni, the redshift of A(1g) mode suggests that the surface segregation of Ni ions takes place. Disorder-activated bands were determined and their integrated intensity evolution with the Ni content suggest that the solid-solution regime favors the increase of disorder; meanwhile, that disorder becomes weaker as the Ni content is increased. Copyright (C) 2010 John Wiley & Sons, Ltd.

Galerkin Method and Weighting Functions Applied to Nonlinear H(infinity) Control with Output Feedback

This paper considers two aspects of the nonlinear H(infinity) control problem: the use of weighting functions for performance and robustness improvement, as in the linear case, and the development of a successive Galerkin approximation method for the solution of the Hamilton-Jacobi-Isaacs equation that arises in the output-feedback case. Design of nonlinear H(infinity) controllers obtained by the well-established Taylor approximation and by the proposed Galerkin approximation method applied to a magnetic levitation system are presented for comparison purposes.

Extension of a Recent Method for Obtaining Exact Solutions of the Bruce and Klute Equation

A method based on a specific power-law relationship between the hydraulic head and the Boltzmann variable was recently presented. We generalized this relationship to a range of powers and extended the solution to include the saturated zone. As a result, the new solution satisfies the Bruce and Klute equation exactly.

Characterization of tree-rings of Pinus caribaea var. hondurensis Barr. et Golf. trees by X ray densitometry

This work aimed to determining the anatomical structure of wood, through methodology of histology and X-ray densitometry, of resin-tapped and not resin-tapped Pinus caribaea var. hondurensis trees samples, of three diameter classes. Pine trees, in forest plantation established in 1969, in the Ecological Experimental Station of Itirapina, from the Forestry Institute of Sao Paulo State, were measured and stratified into three classes of trunk diameter. The pine trees were resin-tapped since 2004, with the opening of two simultaneous and opposing panels. Sixty samples of pine wood trees were extracted from the tree trunk through a non-destructive method and in the laboratory. Tree rings were determined in the laboratory and wood apparent density by X-ray densitometry. The test results showed that: (i) false tree rings occur in the early wood and late wood of the tree rings due to climate change; (ii) the X-ray densitometry allowed the demarcation of the tree rings limits; (iii) the wood apparent density average was significantly different between the trees in high class diameter and in the medium-low class; (iv) the wood characteristics from the resin-tapped and non resin-tapped faces did not show significant differences.

Variation of wood density and relationship with the tree-ring width of mahogany trees, Swietenia macrophylla, from Amazonian tropical forest of Peru

Mahogany trees, Swietenia macrophylla, occur in open rainforest, semi deciduous and deciduous and dense rainforest of Peruvian Amazonian tropical forest. They occur, preferentially, in areas with a defined dry season, with typical phenology and seasonal variation activity, forming distinct tree-rings. The present work had as aim to determine the wood density radial variation of 14 mahogany trees, of two populations of the Peruvian Amazonian tropical forest, through the X-ray densitometry and to evaluate their application as methodology, compared to the classic method of measurement table, for the determination of the treering width. The radial wood apparent density of the trees profiles rendered it possible to delimit the areas of juvenile-adult wood and of the heartwood-sapwood, relative to the anatomical structure and chemical composition differences, due to the extractives and the vessels obstruction by tyloses. The mean, minimum and maximum wood apparent density of the mahogany trees for the Populations A and B were of 0.70; 0.29; 1.01 g.cm(-3) and 0.81; 0.29; 1.19 g.cm(-3), respectively. The analysis of the variance and mean test indicate differences of mean wood density among the mahogany trees of each population, probably due to the age of the trees. There was no correlation between mean wood density of mahogany trees among the two populations, as well as, between the tree-ring width and the respective mean density. The X-ray densitometry technique is an important tool in the evaluation of the radial variation of wood apparent density and the delimitation of tree-ring boundaries, with correlations of 0.94 and 0.93 in relation to measurement table, for each sampled population.

Definition of priority areas for forest conservation through the ordered weighted averaging method

The general objective of this study was to evaluate the ordered weighted averaging (OWA) method, integrated to a geographic information systems (GIS), in the definition of priority areas for forest conservation in a Brazilian river basin, aiming at to increase the regional biodiversity. We demonstrated how one could obtain a range of alternatives by applying OWA, including the one obtained by the weighted linear combination method and, also the use of the analytic hierarchy process (AHP) to structure the decision problem and to assign the importance to each criterion. The criteria considered important to this study were: proximity to forest patches; proximity among forest patches with larger core area; proximity to surface water; distance from roads: distance from urban areas; and vulnerability to erosion. OWA requires two sets of criteria weights: the weights of relative criterion importance and the order weights. Thus, Participatory Technique was used to define the criteria set and the criterion importance (based in AHP). In order to obtain the second set of weights we considered the influence of each criterion, as well as the importance of each one, on this decision-making process. The sensitivity analysis indicated coherence among the criterion importance weights, the order weights, and the solution. According to this analysis, only the proximity to surface water criterion is not important to identify priority areas for forest conservation. Finally, we can highlight that the OWA method is flexible, easy to be implemented and, mainly, it facilitates a better understanding of the alternative land-use suitability patterns. (C) 2008 Elsevier B.V. All rights reserved.

`Swingle` citrumelo propagation by cuttings for citrus nursery tree production or inarching

`Swingle` citrumelo [Citrus paradisi MacFaden x Poncirus trifoliata (L.) Raf.] has been extensively used as a rootstock in several citrus growing regions of the World, including Southern Brazil where `Rangpur` lime (Citrus limonia Osbeck) is still the predominant variety despite being affected by several important pathogens. in this case, `Swingle` citrumelo is used to produce nursery trees to establish new orchards or to be inarched in adult and healthy groves in order to change the rootstock. We report herein a system to produce trees on `Swingle` citrumelo more rapidly by budding onto non-rooted cuttings, as well as assessing potential to rapidly multiply `Swingle` through rooting of non-budded cuttings. Therefore, two potential products are described: budded trees and rooted rootstock cuttings. `Valencia` sweet orange [Citrus sinensis (L.) Osbeck] was budded at different heights on cuttings derived from eight-month old rootstocks. Grafted and additional non-budded cuttings were then treated with indole-3-butyric acid (500 mg L(-1)) or left untreated before rooting. Three types of cuttings were evaluated: softwood, semi-hardwood and hardwood. The use of nursery trees derived from pre-budded hardwood cuttings of `Swingle` citrumelo is an alternative grafting method on this cultivar. Softwood cuttings with one leaf pair were considered the most adequate material for rapid multiplication of `Swingle` citrumelo by cutting. This could be particularly useful for inarching production or conventional budding after transplant of cutting-derived rootstocks. (C) 2008 Elsevier B.V. All rights reserved.

Growth-active gibberellins overcome the very slow shoot growth of Hancornia speciosa, an important fruit tree from the Brazilian ""Cerrado""

Shoot elongation of Hancornia speciosa, an endangered tree from the Brazilian savannah ""Cerrado"", is very slow, thus limiting nursery production of plants. Gibberellins (GAs) A(1), A(3), and A(5), and two inhibitors of GA biosynthesis, trinexapac-ethyl and ancymidol were applied to shoots of Hancornia seedlings. GA(1) and GA(3) significantly stimulated shoot elongation, while GA(5) had no significant effect. Trinexapac-ethyl and ancymidol, both at 100 A mu g per seedling, inhibited shoot elongation up to 45 days after treatment, though the effect was statistically significant only for ancymidol. Somewhat surprisingly, exogenous GA(3) more effectively stimulated shoot elongation in SD-grown plants, than in LD-grown plants. The results from exogenous application of GAs and inhibitors of GA biosynthesis imply that Hancornia shoot growth is controlled by GAs, and that level of endogenous growth-active GAs is likely to be the limiting factor for shoot elongation in Hancornia. Application of GAs thus offer a practical method for nursery production of Hancornia seedlings for outplanting into the field.