SYNCHRONIZATION OF A CLASS OF SECOND-ORDER NONLINEAR SYSTEMS

The asymptotic behavior of a class of coupled second-order nonlinear dynamical systems is studied in this paper. Using very mild assumptions on the vector-field, conditions on the coupling parameters that guarantee synchronization are provided. The proposed result does not require solutions to be ultimately bounded in order to prove synchronization, therefore it can be used to study coupled systems that do not globally synchronize, including synchronization of unbounded solutions. In this case, estimates of the synchronization region are obtained. Synchronization of two-coupled nonlinear pendulums and two-coupled Duffing systems are studied to illustrate the application of the proposed theory.

Experimental results on the nonlinear H(infinity) control via quasi-LPV representation and game theory for wheeled mobile robots

In this paper, nonlinear dynamic equations of a wheeled mobile robot are described in the state-space form where the parameters are part of the state (angular velocities of the wheels). This representation, known as quasi-linear parameter varying, is useful for control designs based on nonlinear H(infinity) approaches. Two nonlinear H(infinity) controllers that guarantee induced L(2)-norm, between input (disturbances) and output signals, bounded by an attenuation level gamma, are used to control a wheeled mobile robot. These controllers are solved via linear matrix inequalities and algebraic Riccati equation. Experimental results are presented, with a comparative study among these robust control strategies and the standard computed torque, plus proportional-derivative, controller.

Application of lumped dissipation model in nonlinear analysis of reinforced concrete structures

This paper deals with the application of the lumped dissipation model in the analysis of reinforced concrete structures, emphasizing the nonlinear behaviour of the materials The presented model is based on the original models developed by Cipollina and Florez-Lopez (1995) [12]. Florez-Lopez (1995) [13] and Picon and Florez-Lopez (2000) [14] However, some modifications were introduced in the functions that control the damage evolution in order to improve the results obtained. The efficiency of the new approach is evaluated by means of a comparison with experimental results on reinforced concrete structures such as simply supported beams, plane frames and beam-to-column connections Finally, the adequacy of the numerical model representing the global behaviour of framed structures is investigated and the limits of the analysis are discussed (C) 2009 Elsevier Ltd All rights reserved

A simple way to introduce fibers into FEM models

This communication proposes a simple way to introduce fibers into finite element modelling. This is a promising formulation to deal with fiber-reinforced composites by the finite element method (FEM), as it 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 into the pre-existent finite element numerical system to consider any distribution 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 is the reduced work required by the user to introduce fibers, avoiding `rebar` elements, node-by-node geometrical definitions or even complex mesh generation. An additional characteristic of the technique is the possibility of representing unbounded stresses at the end of fibers using a finite number of degrees of freedom. Further studies are required for non-linear applications in which localization may occur. Along the text the linear formulation is presented and the bounded connection between fibers and continuum is considered. Four examples are presented, including non-linear analysis, to validate and show the capabilities of the formulation. Copyright (c) 2007 John Wiley & Sons, Ltd.

A numerical model for the description of the nonlinear behaviour of multi-leaf masonry walls

A nonlinear finite element model was developed to simulate the nonlinear response of three-leaf masonry specimens, which were subjected to laboratory tests with the aim of investigating the mechanical behaviour of multiple-leaf stone masonry walls up to failure. The specimens consisted of two external leaves made of stone bricks and mortar joints, and an internal leaf in mortar and stone aggregate. Different loading conditions, typologies of the collar joints, and stone types were taken into account. The constitutive law implemented in the model is characterized by a damage tensor, which allows the damage-induced anisotropy accompanying the cracking process to be described. To follow the post-peak behaviour of the specimens with sufficient accuracy it was necessary to make the damage model non-local, to avoid mesh-dependency effects related to the strain-softening behaviour of the material. Comparisons between the predicted and measured failure loads are quite satisfactory in most of the studied cases. (c) 2007 Elsevier Ltd. All rights reserved.

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.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

Characterization and kinetics of sulfide-oxidizing autotrophic denitrification in batch reactors containing suspended and immobilized cells

Sulfide-oxidizing autotrophic denitrification is an advantageous alternative over heterotrophic denitrification, and may have potential for nitrogen removal of low-strength wastewaters, such as anaerobically pre-treated domestic sewage. This study evaluated the fundamentals and kinetics of this process in batch reactors containing suspended and immobilized cells. Batch tests were performed for different NO(x)(-)/S(2-) ratios and using nitrate and nitrite as electron acceptors. Autotrophic denitrification was observed for both electron acceptors, and NO(x)(-)/S(2-) ratios defined whether sulfide oxidation was complete or not. Kinetic parameter values obtained for nitrate were higher than for nitrite as electron acceptor. Zero-order models were better adjusted to profiles obtained for suspended cell reactors, whereas first-order models were more adequate for immobilized cell reactors. However, in the latter, mass transfer physical phenomena had a significant effect on kinetics based on biochemical reactions. Results showed that sulfide-oxidizing autotrophic denitrification can be successfully established for low-strength wastewaters and have potential for nitrogen removal from anaerobically pre-treated domestic sewage.

Analytical solutions for geometrically nonlinear trusses

This paper presents an analytical method for analyzing trusses with severe geometrically nonlinear behavior. The main objective is to find analytical solutions for trusses with different axial forces in the bars. The methodology is based on truss kinematics, elastic constitutive laws and equilibrium of nodal forces. The proposed formulation can be applied to hyper elastic materials, such as rubber and elastic foams. A Von Mises truss with two bars made by different materials is analyzed to show the accuracy of this methodology.

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

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

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.

Considerations on the kinetics of froth flotation of ultrafine coal contained in tailings

This article presents a kinetic evaluation of froth flotation of ultrafine coal contained in the tailings from a Colombian coal preparation plant. The plant utilizes a dense-medium cyclones and spirals circuit. The tailings contained material that was 63% finer than 14 mu m. Flotation tests were performed with and without coal ""promoters"" (diesel oil or kerosene) to evaluate the kinetics of flotation of coal. It was found that flotation rates were higher when no promoter was added. Different kinetic models were evaluated for the flotation of the coal from the tailings, and it was found that the best fitted model was the classical first-order model.

Importância da hidrodinâmica na cinética de flotação de partículas grossas

The processes that govern the rate of particle recovery in a flotation cell include the following sub-processes: collision, attachment, and stability of the aggregate formed by particles and bubbles. Collision is controlled by bulk hydrodynamics inside the flotation cell, while attachment is largely dominated by variables that belong to the domain of surface chemistry (contact angle, induction time). As for the stability of the particle/bubble aggregate, its efficiency depends on both hydrodynamics plus surface chemistry variables of the system. The flotation recovery of coarse particles of apatite and glass spheres was measured by micro-flotation and batch flotation tests in which hydrodynamic parameters were evaluated, such as impeller rotational speed, diameter, and geometry, as well as particle size and density. Results revealed that a proper impeller rotational speed yielded turbulence levels, which enabled to keep particles fully suspended, this way optimizing the collision efficiency between particles and bubbles, without jeopardizing the stability of the particle-bubble aggregates.

Kinetics and Modeling of Fatty Alcohol Ethoxylation in an Industrial Spray Loop Reactor

The kinetics of the ethoxylation of fatty alcohols catalyzed by potassium hydroxide was studied to obtain the rate constants for modeling of the industrial process. Experimental data obtained in a lab-scale semibatch autoclave reactor were used to evaluate kinetic and equilibrium parameters. The kinetic model was employed to model the performance of an industrial-scale spray tower reactor for fatty alcohol ethoxylation. The reactor model considers that mass transfer and reaction occur independently in two distinct zones of the reactor. Good agreement between the model predictions and real data was found. These findings confirm the reliability of the kinetic and reactor model for simulating fatty alcohol ethoxylation processes under industrial conditions.

Hydrolysis of acetic anhydride: Non-adiabatic calorimetric determination of kinetics and heat exchange

A simple calorimetric method to estimate both kinetics and heat transfer coefficients using temperature-versus-time data under non-adiabatic conditions is described for the reaction of hydrolysis of acetic anhydride. The methodology is applied to three simple laboratory-scale reactors in a very simple experimental setup that can be easily implemented. The quality of the experimental results was verified by comparing them with literature values and with predicted values obtained by energy balance. The comparison shows that the experimental kinetic parameters do not agree exactly with those reported in the literature, but provide a good agreement between predicted and experimental data of temperature and conversion. The differences observed between the activation energy obtained and the values reported in the literature can be ascribed to differences in anhydride-to-water ratios (anhydride concentrations). (C) 2010 Elsevier Ltd. All rights reserved.