A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Karman`s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.

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.

Stepped and smooth spillways: resistance effects on stilling basin lengths

The Darcy-Weisbach equation was used in the analysis of flow over spillways, furnishing theoretical tools to design stilling basins. Predictions for the length of hydraulic jump stilling basins downstream of stepped and smooth spillways are presented, together with ranges of values for the Darcy-Weisbach friction factor of both spillways. The experimental data were compared with results of the theoretical solution of the gradually varied flow equation. All comparisons were made in non-dimensional form. The values of the Darcy-Weisbach friction factor were roughly five times smaller for smooth spillways than for stepped spillways. The theoretical predictions and the experimental data allow to present approximate equations for a preliminary evaluation of the length and the bed level of hydraulic jump stilling basins. In the same way, approximate equations were presented for the evaluation of the friction factor in smooth and stepped spillways, as a function of the Froude number at the downstream cross-section.

An upflow fixed-bed anaerobic-aerobic reactor for removal of organic matter and nitrogen from L-lysine plant wastewater

This paper reports on the design of a new reactor configuration - an upflow fixed-bed combined anaerobic-aerobic reactor - can operate as a single treatment unit for the removal of nitrogen (approximate to 150 mg N/L) and organic matter (approximate to 1300 mg COD/L) from Lysine plant wastewater. L-Lysine, an essential amino acid for animal nutrition, is produced by fermentation from natural raw materials of agricultural origin, thus generating wastewater with high contents of organic matter and nitrogen. The best operational condition of the reactor was obtained with a hydraulic retention time of 35 h (21 h in the anaerobic zone and 14 h in the aerobic zone) and a recycling ratio (R) of 3.5. In this condition, the COD, total Kjeldahl nitrogen (TKN), and total nitrogen (TN) removal efficiencies were 97%, 96%, and 77%, respectively, with average effluent concentrations of 10 +/- 36 mg COD/L, 2 +/- 1 mg NH(4)(+)-N/L, 8 +/- 3 mg Org-N/L, 1 +/- 1 mg NH(2)(-)-N/L, and 26 +/- 23 mg NH(3)(-)-N/L.

Vegetable Oil Quenchants: Calculation and Comparison of The Cooling Properties of a Series of Vegetable Oils

The compositions of canola, soybean, corn, cottonseed and sunflower oils suggest that they exhibit substantially different propensity for oxidation following the order of Canola < corn < cottonseed < sunflower approximate to soybean. These data suggest that any of the vegetable oils evaluated could be blended with minimal impact on viscosity although compositional differences would surely affect oxidative stability. Cooling curve analysis showed that similar cooling profiles were obtained for different vegetable oils. Interestingly, no film boiling or transition nucleate boiling was observed with any of the vegetable oils and heat transfer occurs only by pure nucleate boiling and convection. High-temperature cooling properties of vegetable oils are considerable faster than those observed for petroleum oil-based quenchants. (C)2010 Journal of Mechanical Engineering. All rights reserved.

VEGETABLE OIL STRUCTURE AND ANTIOXIDANTS

Five vegetable oils: canola, soybean, corn, cottonseed and sunflower oils were characterized with respect to their composition by gas chromatography and viscosity. The compositions of the vegetable oils suggest that they exhibit substantially different propensity for oxidation following the order of: canola < corn < cottonseed < sunflower approximate to soybean. Viscosities at 40 degrees C and 100 degrees C and the viscosity index (VI) values were determined for the vegetable oils and two petroleum oil quenchants: Microtemp 157 (a conventional slow oil) and Microtemp 153B (an accelerated or fast oil). The kinematic viscosities of the different vegetable and petroleum oils at 40 degrees C were similar. The VI values for the different vegetable oils were very close and varied between 209-220 and were all much higher than the VI values obtained for Microtemp 157 (96) and Microtemp 153B (121). These data indicate that the viscosity variations of these vegetable oils are substantially less sensitive to temperature variation than are the parafinic oil based Microtemp 157 and Microtemp 153B. Although these data suggest that any of the vegetable oils evaluated could be blended with minimal impact on viscosity, the oxidative stability would surely be substantially impacted. Cooling curve analysis was performed on these vegetable oils at 60 degrees C under non-agitated conditions. These results were compared with cooling curves obtained for Microtemp 157, a conventional, unaccelerated petroleum oil, and Microtemp 153B, an accelerated petroleum oil under the same conditions. The results showed that cooling profiles of the different vegetable oils were similar as expected from the VI values. However, no boiling was observed wit any of the vegetable oils and heat transfer occurs only by convection since there is no full-film boiling and nucleate boiling process as typically observed for petroleum oil quenchants, including those of this study. Therefore, high-temperature cooling is considerable faster for vegetable oils as a class. The cooling properties obtained suggest that vegetable oils would be especially suitable fur quenching low-hardenability steels such as carbon steels.

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.

A transient three-dimensional heat transfer model of the human body

The objective of this work is to develop an improved model of the human thermal system. The features included are important to solve real problems: 3D heat conduction, the use of elliptical cylinders to adequately approximate body geometry, the careful representation of tissues and important organs, and the flexibility of the computational implementation. Focus is on the passive system, which is composed by 15 cylindrical elements and it includes heat transfer between large arteries and veins. The results of thermal neutrality and transient simulations are in excellent agreement with experimental data, indicating that the model represents adequately the behavior of the human thermal system. (C) 2009 Elsevier Ltd. All rights reserved.

Sphericity of apatite particles determined by gas permeability through packed beds

Because shape is an assessment of the three-dimensional form of a particle, it may be described in terms of sphericity (Psi), which is a measure of how closely a particle approaches a spherical configuration. In this study, Darcy`s law and the Kozeny-Carman model for fluid flow through porous media were applied to packed beds to determine the sphericity (Psi) of apatite particles. The beds were composed of glass spheres or particles of apatite (igneous from Brazil and sedimentary from the United States) of three classes of size (Class 1: -297 +210 mu m; Class 2: -210 +149 mu m; Class 3: -149 +105 mu m). Glass spheres were used to validate the model because of its known sphericity (Psi = 1.00). Apatite particles, either igneous or sedimentary, showed very close values for particle sphericity (Psi approximate to 0.6). Observations on particle images conducted by scanning electron microscopy illustrated that igneous (Psi = 0.623) and sedimentary (Psi = 0.644) particles of apatite of Class 2 predominantly exhibit elongated shape. The close value of particle sphericity (Psi approximate to 0.6) showed by either igneous or sedimentary apatite may be justified by the similarity in particle shape.

Surface tension of flotation solution and its influence on the selectivity of the separation between apatite and gangue minerals

This work presents and discusses the influence of the surface tension (gamma(LV)) of methanol/water mixtures on the flotation response of apatite versus gangue minerals conditioned with flotation reagents (75 g/t cornstarch and 100 g/t Berol 867) at pH 10.6. Berol 867 is a collector composed of sodium alkyl sarcosinate plus nonionic surfactant. The highest Schulz efficiency of separation (recovery of apatite minus recovery of gangue) was achieved at approximate to 51.0 mN/m. The critical surface tension of wettability (gamma(C)) of apatite was found to occur at 34.7 mN/m when determined by means of gamma flotation experiments, , and it was 33.9 mN/m when determined by Zisman`s approach.

Concentration inequalities and laws of large numbers under epistemic and regular irrelevance

This paper presents concentration inequalities and laws of large numbers under weak assumptions of irrelevance that are expressed using lower and upper expectations. The results build upon De Cooman and Miranda`s recent inequalities and laws of large numbers. The proofs indicate connections between the theory of martingales and concepts of epistemic and regular irrelevance. (C) 2010 Elsevier Inc. All rights reserved.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Probabilistic logic with independence

This paper investigates probabilistic logics endowed with independence relations. We review propositional probabilistic languages without and with independence. We then consider graph-theoretic representations for propositional probabilistic logic with independence; complexity is analyzed, algorithms are derived, and examples are discussed. Finally, we examine a restricted first-order probabilistic logic that generalizes relational Bayesian networks. (c) 2007 Elsevier Inc. All rights reserved.

Impact of plastic deformation on magnetoacoustic properties of Fe-2%Si alloy

The paper presents the results of a complementary study including magnetic hysteresis loops B(H), magnetic Barkhausen noise (MBN) and magnetoacoustic emission (MAE) signals measurements for plastically deformed Fe-2%Si samples. The investigated samples had been plastically deformed with plastic strain level (epsilon(p)) up to 8%. The properties of B(H) loops are quantified using the coercivity H(C) and maximum differential permeability mu(rmax) as parameters. The MBN and MAE voltage signals were analysed by means of rms-like voltage (Ub and Ua, respectively) envelopes, plotted as a function of applied field strength. Integrals of the Ub and Ua voltages over half of a period of magnetization were then calculated. It has been found that He and integrals of Ub increase, while mu(rmax) decreases monotonically with increasing epsilon(p). The MAE (Ua) peak voltage at first decreases, then peaks at epsilon(p) approximate to 1.5% and finally decreases again. The integral of the Ua voltage at first increases for low epsilon(p) and then decreases for epsilon(p) > 1.5%. All those various dependence types suggest the possibility of detection of various stages of microstructure change. The above-mentioned results are discussed qualitatively in the paper. Some modelling of the discussed dependency is also presented. (C) 2008 Elsevier Ltd. All rights reserved.

Anomalous compressibility of ferropericlase throughout the iron spin cross-over

The thermoelastic properties of ferropericlase Mg(1-x)Fe(x)O (x = 0.1875) throughout the iron high-to-low spin cross-over have been investigated by first principles at Earth`s lower mantle conditions. This cross-over has important consequences for elasticity such as an anomalous bulk modulus (K(S)) reduction. At room temperature the anomaly is somewhat sharp in pressure but broadens with increasing temperature. Along a typical geotherm it occurs across most of the lower mantle with a more significant K(S) reduction at approximate to 1,400-1,600 km depth. This anomaly might also cause a reduction in the effective activation energy for diffusion creep and lead to a viscosity minimum in the mid-lower mantle, in apparent agreement with results from inversion of data related with mantle convection and postglacial rebound.