A study on the pretreatment of a sugarcane bagasse sample with dilute sulfuric acid

Experiments based on a 2(3) central composite full factorial design were carried out in 200-ml stainless-steel containers to study the pretreatment, with dilute sulfuric acid, of a sugarcane bagasse sample obtained from a local sugar-alcohol mill. The independent variables selected for study were temperature, varied from 112.5A degrees C to 157.5A degrees C, residence time, varied from 5.0 to 35.0 min, and sulfuric acid concentration, varied from 0.0% to 3.0% (w/v). Bagasse loading of 15% (w/w) was used in all experiments. Statistical analysis of the experimental results showed that all three independent variables significantly influenced the response variables, namely the bagasse solubilization, efficiency of xylose recovery in the hemicellulosic hydrolysate, efficiency of cellulose enzymatic saccharification, and percentages of cellulose, hemicellulose, and lignin in the pretreated solids. Temperature was the factor that influenced the response variables the most, followed by acid concentration and residence time, in that order. Although harsher pretreatment conditions promoted almost complete removal of the hemicellulosic fraction, the amount of xylose recovered in the hemicellulosic hydrolysate did not exceed 61.8% of the maximum theoretical value. Cellulose enzymatic saccharification was favored by more efficient removal of hemicellulose during the pretreatment. However, detoxification of the hemicellulosic hydrolysate was necessary for better bioconversion of the sugars to ethanol.

Characterization of commercial cellulases and their use in the saccharification of a sugarcane bagasse sample pretreated with dilute sulfuric acid

This study aimed to correlate the efficiency of enzymatic hydrolysis of the cellulose contained in a sugarcane bagasse sample pretreated with dilute H(2)SO(4) with the levels of independent variables such as initial content of solids and loadings of enzymes and surfactant (Tween 20), for two cellulolytic commercial preparations. The preparations, designated cellulase I and cellulase II, were characterized regarding the activities of total cellulases, endoglucanase, cellobiohydrolase, cellobiase, beta-glucosidase, xylanase, and phenoloxidases (laccase, manganese and lignin peroxidases), as well as protein contents. Both extracts showed complete cellulolytic complexes and considerable activities of xylanases, without activities of phenoloxidases. For the enzymatic hydrolyses, two 2(3) central composite full factorial designs were employed to evaluate the effects caused by the initial content of solids (1.19-4.81%, w/w) and loadings of enzymes (1.9-38.1 FPU/g bagasse) and Tween 20 (0.0-0.1 g/g bagasse) on the cellulose digestibility. Within 24 h of enzymatic hydrolysis, all three independent variables influenced the conversion of cellulose by cellulase I. Using cellulase II, only enzyme and surfactant loadings showed significant effects on cellulose conversion. An additional experiment demonstrated the possibility of increasing the initial content of solids to values much higher than 4.81% (w/w) without compromising the efficiency of cellulose conversion, consequently improving the glucose concentration in the hydrolysate.

Aerodynamic study of three-dimensional larynx models using finite element methods

The airflow velocities and pressures are calculated from a three-dimensional model of the human larynx by using the finite element method. The laryngeal airflow is assumed to be incompressible, isothermal, steady, and created by fixed pressure drops. The influence of different laryngeal profiles (convergent, parallel, and divergent), glottal area, and dimensions of false vocal folds in the airflow are investigated. The results indicate that vertical and horizontal phase differences in the laryngeal tissue movements are influenced by the nonlinear pressure distribution across the glottal channel, and the glottal entrance shape influences the air pressure distribution inside the glottis. Additionally, the false vocal folds increase the glottal duct pressure drop by creating a new constricted channel in the larynx, and alter the airflow vortexes formed after the true vocal folds. (C) 2007 Elsevier Ltd. All rights reserved.

Finite element homogenization technique for the characterization of d(15) shear piezoelectric macro-fibre composites

A finite element homogenization method for a shear actuated d(15) macro-fibre composite (MFC) made of seven layers (Kapton, acrylic, electrode, piezoceramic fibre and epoxy composite, electrode, acrylic, Kapton) is proposed and used for the characterization of its effective material properties. The methodology is first validated for the MFC active layer only, made of piezoceramic fibre and epoxy, through comparison with previously published analytical results. Then, the methodology is applied to the seven-layer MFC. It is shown that the packaging reduces significantly the shear stiffness of the piezoceramic material and, thus, leads to significantly smaller effective electromechanical coupling coefficient k(15) and piezoelectric stress constant e(15) when compared to the piezoceramic fibre properties. However, it is found that the piezoelectric charge constant d(15) is less affected by the softer layers required by the MFC packaging.

Improved finite element for 3D laminate frame analysis including warping for any cross-section

The most ordinary finite element formulations for 3D frame analysis do not consider the warping of cross-sections as part of their kinematics. So the stiffness, regarding torsion, should be directly introduced by the user into the computational software and the bar is treated as it is working under no warping hypothesis. This approach does not give good results for general structural elements applied in engineering. Both displacement and stress calculation reveal sensible deficiencies for both linear and non-linear applications. For linear analysis, displacements can be corrected by assuming a stiffness that results in acceptable global displacements of the analyzed structure. However, the stress calculation will be far from reality. For nonlinear analysis the deficiencies are even worse. In the past forty years, some special structural matrix analysis and finite element formulations have been proposed in literature to include warping and the bending-torsion effects for 3D general frame analysis considering both linear and non-linear situations. In this work, using a kinematics improvement technique, the degree of freedom ""warping intensity"" is introduced following a new approach for 3D frame elements. This degree of freedom is associated with the warping basic mode, a geometric characteristic of the cross-section, It does not have a direct relation with the rate of twist rotation along the longitudinal axis, as in existent formulations. Moreover, a linear strain variation mode is provided for the geometric non-linear approach, for which complete 3D constitutive relation (Saint-Venant Kirchhoff) is adopted. The proposed technique allows the consideration of inhomogeneous cross-sections with any geometry. Various examples are shown to demonstrate the accuracy and applicability of the proposed formulation. (C) 2009 Elsevier Inc. 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 layered finite element for reinforced concrete beams with bond-slip effects

The behaviour of reinforced concrete members is affected by the slipping of steel bars inserted in the concrete matrix. A tension-stiffening effect and crack evolution occur from the beginning of slipping; thus, the assessment of those phenomena requires the introduction of a bond-slip interaction model. This work presents a beam-layered model, including the constitutive relationships of materials and their interaction, according to the CEB-FIP Model Code 1990. To eliminate the finite element sub-division procedure, a continuous slip function is imposed into the element domain. The results are continuous descriptions of bond stress in the steel-concrete interface, as well as concrete and steel stresses along the element. (C) 2007 Elsevier Ltd. 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.

An electromechanical finite element model for piezoelectric energy harvester plates

Vibration-based energy harvesting has been investigated by several researchers over the last decade. The goal in this research field is to power small electronic components by converting the waste vibration energy available in their environment into electrical energy. Recent literature shows that piezoelectric transduction has received the most attention for vibration-to-electricity conversion. In practice, cantilevered beams and plates with piezoceramic layers are employed as piezoelectric energy harvesters. The existing piezoelectric energy harvester models are beam-type lumped parameter, approximate distributed parameter and analytical distributed parameter solutions. However, aspect ratios of piezoelectric energy harvesters in several cases are plate-like and predicting the power output to general (symmetric and asymmetric) excitations requires a plate-type formulation which has not been covered in the energy harvesting literature. In this paper. an electromechanically coupled finite element (FE) plate model is presented for predicting the electrical power output of piezoelectric energy harvester plates. Generalized Hamilton`s principle for electroelastic bodies is reviewed and the FE model is derived based on the Kirchhoff plate assumptions as typical piezoelectric energy harvesters are thin structures. Presence of conductive electrodes is taken into account in the FE model. The predictions of the FE model are verified against the analytical solution for a unimorph cantilever and then against the experimental and analytical results of a bimorph cantilever with a tip mass reported in the literature. Finally, an optimization problem is solved where the aluminum wing spar of an unmanned air vehicle (UAV) is modified to obtain a generator spar by embedding piezoceramics for the maximum electrical power without exceeding a prescribed mass addition limit. (C) 2009 Elsevier Ltd. All rights reserved.

Collapse behavior of soil in a Brazilian region affected by a rising water table

Collapsible soils are usually nonsaturated, low density, and metastable-structured soils that are known to exhibit a volume reduction following an episode of moisture increase or suction reduction. This paper describes the collapsible behavior of clayey sand based on controlled soil suction tests carried out on undisturbed samples from the city of Pereira Barreto, in the State of Sao Paulo, Brazil. Foundation settlements due to soil collapse are common in this region and occurred during the filling of the reservoir of the Tres Irmaos Dam, which induced the elevation of the groundwater table in different parts of Pereira Barreto. This paper shows that collapse strains depend on the stress and soil suction acting in the sample and that saturation is not necessary for a collapse to occur. The influence of soil suction, gradual wetting, and the wetting and drying cycle on the collapsible behavior of the soil is also shown and discussed.

Efficient Modeling of Thin Wires in a Lossy Medium by Finite Elements Applied to Grounding Systems

A procedure is proposed to accurately model thin wires in lossy media by finite element analysis. It is based on the determination of a suitable element width in the vicinity of the wire, which strongly depends on the wire radius to yield accurate results. The approach is well adapted to the analysis of grounding systems. The numerical results of the application of finite element analysis with the suitably chosen element width are compared with both analytical results and those computed by a commercial package for the analysis of grounding systems, showing very good agreement.

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

Finite element computational modeling of externally bonded CFRP composites flexural behavior in RC beams

This paper focuses on the flexural behavior of RC beams externally strengthened with Carbon Fiber Reinforced Polymers (CFRP) fabric. A non-linear finite element (FE) analysis strategy is proposed to support the beam flexural behavior experimental analysis. A development system (QUEBRA2D/FEMOOP programs) has been used to accomplish the numerical simulation. Appropriate constitutive models for concrete, rebars, CFRP and bond-slip interfaces have been implemented and adjusted to represent the composite system behavior. Interface and truss finite elements have been implemented (discrete and embedded approaches) for the numerical representation of rebars, interfaces and composites.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

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.