Geometric nonlinear dynamic analysis of plates and shells using eight-node hexahedral finite elements with reduced integration

This work presents a geometric nonlinear dynamic analysis of plates and shells using eight-node hexahedral isoparametric elements. The main features of the present formulation are: (a) the element matrices are obtained using reduced integrations with hourglass control; (b) an explicit Taylor-Galerkin scheme is used to carry out the dynamic analysis, solving the corresponding equations of motion in terms of velocity components; (c) the Truesdell stress rate tensor is used; (d) the vector processor facilities existing in modern supercomputers were used. The results obtained are comparable with previous solutions in terms of accuracy and computational performance.

Pulse response based control for positive definite systems

Pulse Response Based Control (PRBC) is a recently developed minimum time control method for flexible structures. The flexible behavior of the structure is represented through a set of discrete time sequences, which are the responses of the structure due to rectangular force pulses. The rectangular force pulses are given by the actuators that control the structure. The set of pulse responses, desired outputs, and force bounds form a numerical optimization problem. The solution of the optimization problem is a minimum time piecewise constant control sequence for driving the system to a desired final state. The method was developed for driving positive semi-definite systems. In case the system is positive definite, some final states of the system may not be reachable. Necessary conditions for reachability of the final states are derived for systems with a finite number of degrees of freedom. Numerical results are presented that confirm the derived analytical conditions. Numerical simulations of maneuvers of distributed parameter systems have shown a relationship between the error in the estimated minimum control time and sampling interval

Non-linear analysis of laminated metal matrix composites by an integrated micro/macro-mechanical model

The present paper describes an integrated micro/macro mechanical study of the elastic-viscoplastic behavior of unidirectional metal matrix composites (MMC). The micromechanical analysis of the elastic moduli is based on the Composites Cylinder Assemblage model (CCA) with comparisons also draw with a Representative Unit Cell (RUC) technique. These "homogenization" techniques are later incorporated into the Vanishing Fiber Diameter (VFD) model and a new formulation is proposed. The concept of a smeared element procedure is employed in conjunction with two different versions of the Bodner and Partom elastic-viscoplastic constitutive model for the associated macroscopic analysis. The formulations developed are also compared against experimental and analytical results available in the literature.

A temperature predicting model for manufacturing processes requiring coiling

A model for predicting temperature evolution for automatic controling systems in manufacturing processes requiring the coiling of bars in the transfer table is presented. Although the method is of a general nature, the presentation in this work refers to the manufacturing of steel plates in hot rolling mills. The predicting strategy is based on a mathematical model of the evolution of temperature in a coiling and uncoiling bar and is presented in the form of a parabolic partial differential equation for a shape changing domain. The mathematical model is solved numerically by a space discretization via geometrically adaptive finite elements which accomodate the change in shape of the domain, using a computationally novel treatment of the resulting thermal contact problem due to coiling. Time is discretized according to a Crank-Nicolson scheme. Since the actual physical process takes less time than the time required by the process controlling computer to solve the full mathematical model, a special predictive device was developed, in the form of a set of least squares polynomials, based on the off-line numerical solution of the mathematical model.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

The snail hosts of schistosomiasis: some evolutionary and ecological perspectives in relation to control

Despite opportunities for radiation provided by spatio-temporal isolation, the basic morphological plan of pulmonate snails has remained conservative. In consequence of the resulting dearth of morphological characters and their plasticity, there is a case for using biochemical characters such as exogenous chemicals released by the snails (e.g. amino acids) and their chemoreception niche as taxonomic aids to classify snails of medical importance. As these same chemicals are used by snails to distinguish conspecifics they could also be used as "environmental antibodies" in controlled release formulations (CRF's) designed to remove target snails in a specific, cost-effective and ecologically acceptable manner. The snails, surface-living bacteria, algae and macrophytic plants are considered as co-evolved, interactive modular systems with strong mutualistic elements. Recently, anthropogenic perturbations such as deforestation, and damming of flowing waters, have benefited these modules whereas others such as river canalization, acid deposition, accumulation of pesticide residues and eutrophication have harmed them. Research is needed to elucidate the factors which limit the growth of snails in primitive habitats, uninfluenced by man, as well as in those subject to harmful anthropogenic factors. The understanding thus gained could be applied to develop cost-effective primary health care strategies to reduce or prevent transmission of schistosomiasis and other water related diseases.

Modelagem e simulação da desidratação osmótica em pedaços de abacaxi utilizando o método de elementos finitos

Water loss and sugar gain were modelling during the osmotic dehydration process of pieces of pineaplle. The transfer of solute to the fruit and the water to the solution was based on Fick's 2nd law. The three dimensional model was solved by the finite element method with the usage of the software COMSOL Multiphysics 3.2. The main and cross diffusion coefficients and the Biot number were determined on the simulation and the deviation between the experimental and the simulated data were 4,28% to sucrose and 1,66 to the water.

ALLOMETRIC MODELS FOR ESTIMATING ABOVEGROUND BIOMASS AND BIOMASS ALLOCATION OF CAPIXINGUI TREES (Croton floribundus Spreng.) IN AN AGRISILVICULTURAL SYSTEM

ABSTRACT The objective of this study was to select allometric models to estimate total and pooled aboveground biomass of 4.5-year-old capixingui trees established in an agrisilvicultural system. Aboveground biomass distribution of capixingui was also evaluated. Single- (diameter at breast height [DBH] or crown diameter or stem diameter as the independent variable) and double-entry (DBH or crown diameter or stem diameter and total height as independent variables) models were studied. The estimated total biomass was 17.3 t.ha-1, corresponding to 86.6 kg per tree. All models showed a good fit to the data (R2ad > 0.85) for bole, branches, and total biomass. DBH-based models presented the best residual distribution. Model lnW = b0 + b1* lnDBH can be recommended for aboveground biomass estimation. Lower coefficients were obtained for leaves (R2ad > 82%). Biomass distribution followed the order: bole>branches>leaves. Bole biomass percentage decreased with increasing DBH of the trees, whereas branch biomass increased.

Evaluation of the heating of broilers in shed by means of gas hoods and underfloor heating

This study was conducted to evaluate and compare the behavior of the heating inside poultry shed through gas hood and underfloor heating. The experiment was conducted in poultry shed belonging to the Federal Institute of Education, Science and Technology of Triângulo Mineiro Region, Uberlândia city – state of Minas Gerais (MG), Brazil. The dimensions of the shed are 24 meters long and 9.6 meters wide and with a ceiling height of 3.2 meters. The temperature was measured with an optical thermometer of Minipa brand, MT 350 model. It was used, to the analysis of temperature behavior, the public domain software FEMM 4.2, which uses finite elements techniques, with data collected from two lots. Underfloor heating is made using hot water flowing through a serpentine type system, which is installed below the bed; this hot water is from solar heaters. An energetic and economic assessment of the warming shed for raising chickens was realized. From the results obtained with the simulations, it may observe that the heating through the floor provides a more homogeneous distribution of temperature when compared with the hood heating. The flow of heat is upwards supplying, thus, the greatest need of heating of the bird, which is the pectoral part.

On the numerical simulation of machining processes

Numerical simulation of machining processes can be traced back to the early seventies when finite element models for continuous chip formation were proposed. The advent of fast computers and development of new techniques to model large plastic deformations have favoured machining simulation. Relevant aspects of finite element simulation of machining processes are discussed in this paper, such as solution methods, material models, thermo-mechanical coupling, friction models, chip separation and breakage strategies and meshing/re-meshing strategies.

Numerical study of Petrov-Galerkin formulations for the shallow water wave equations

The behavior of Petrov-Galerkin formulations for shallow water wave equations is evaluated numerically considering typical one-dimensional propagation problems. The formulations considered here use stabilizing operators to improve classical Galerkin approaches. Their advantages and disadvantages are pointed out according to the intrinsic time scale (free parameter) which has a particular importance in this kind of problem. The influence of the Courant number and the performance of the formulation in dealing with spurious oscillations are adressed.

Contact with friction using the augmented Lagrangian Method: a conditional constrained minimization problem

This work presents a formulation of the contact with friction between elastic bodies. This is a non linear problem due to unilateral constraints (inter-penetration of bodies) and friction. The solution of this problem can be found using optimization concepts, modelling the problem as a constrained minimization problem. The Finite Element Method is used to construct approximation spaces. The minimization problem has the total potential energy of the elastic bodies as the objective function, the non-inter-penetration conditions are represented by inequality constraints, and equality constraints are used to deal with the friction. Due to the presence of two friction conditions (stick and slip), specific equality constraints are present or not according to the current condition. Since the Coulomb friction condition depends on the normal and tangential contact stresses related to the constraints of the problem, it is devised a conditional dependent constrained minimization problem. An Augmented Lagrangian Method for constrained minimization is employed to solve this problem. This method, when applied to a contact problem, presents Lagrange Multipliers which have the physical meaning of contact forces. This fact allows to check the friction condition at each iteration. These concepts make possible to devise a computational scheme which lead to good numerical results.

Response of Grilage-Like Structures by the Mobility Approach

This work presents an application of the Mobility Approach to the analysis of the power flow through grillage-like structures. Such structures are usually found in offshore platforms, supporting large and heavy machines. Different wave kinds (longitudinal, flexural and torsional) were initially considered in the power flow analysis between two beams joined in L. Beams excited by an in-plane point force showed strong coupling between longitudinal-flexural waves, while that for out-of-plane point force excitation, flexural-torsional waves coupling represents the most important mechanism of energy transmission. The response determination of grillages by the mobility approach requires the structure to be separated into simple beam-like structural components. Equations for rotations and displacements at the joints of all beams are written for as mobility functions, and moments and forces acting at the joints. A system of equations relating all such internal forces and moments is obtained. This approach was applied to simple grillages. Response results showed good agreement when compared to those provided by Finite Elements.

Simulation of Steady-State Nonlinear Heat Transfer Problems Through the Minimization of Quadratic Functionals

In this work it is presented a systematic procedure for constructing the solution of a large class of nonlinear conduction heat transfer problems through the minimization of quadratic functionals like the ones usually employed for linear descriptions. The proposed procedure gives rise to an efficient and easy way for carrying out numerical simulations of nonlinear heat transfer problems by means of finite elements. To illustrate the procedure a particular problem is simulated by means of a finite element approximation.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).