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.

Multiple random crack propagation using a boundary element formulation

This paper proposes a boundary element method (BEM) model that is used for the analysis of multiple random crack growth by considering linear elastic fracture mechanics problems and structures subjected to fatigue. The formulation presented in this paper is based on the dual boundary element method, in which singular and hyper-singular integral equations are used. This technique avoids singularities of the resulting algebraic system of equations, despite the fact that the collocation points coincide for the two opposite crack faces. In fracture mechanics analyses, the displacement correlation technique is applied to evaluate stress intensity factors. The maximum circumferential stress theory is used to evaluate the propagation angle and the effective stress intensity factor. The fatigue model uses Paris` law to predict structural life. Examples of simple and multi-fractured structures loaded until rupture are considered. These analyses demonstrate the robustness of the proposed model. In addition, the results indicate that this formulation is accurate and can model localisation and coalescence phenomena. (C) 2010 Elsevier Ltd. All rights reserved.

On the analysis of notched concrete beams: From measurement with digital image correlation to identification with boundary element method of a cohesive model

A way of coupling digital image correlation (to measure displacement fields) and boundary element method (to compute displacements and tractions along a crack surface) is presented herein. It allows for the identification of Young`s modulus and fracture parameters associated with a cohesive model. This procedure is illustrated to analyze the latter for an ordinary concrete in a three-point bend test on a notched beam. In view of measurement uncertainties, the results are deemed trustworthy thanks to the fact that numerous measurement points are accessible and used as entries to the identification procedure. (C) 2010 Elsevier Ltd. 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.

Modeling and Analysis of Piezoelectric Energy Harvesting From Aeroelastic Vibrations Using the Doublet-Lattice Method

Multifunctional structures are pointed out as an important technology for the design of aircraft with volume, mass, and energy source limitations such as unmanned air vehicles (UAVs) and micro air vehicles (MAVs). In addition to its primary function of bearing aerodynamic loads, the wing/spar structure of an UAV or a MAV with embedded piezoceramics can provide an extra electrical energy source based on the concept of vibration energy harvesting to power small and wireless electronic components. Aeroelastic vibrations of a lifting surface can be converted into electricity using piezoelectric transduction. In this paper, frequency-domain piezoaeroelastic modeling and analysis of a canti-levered platelike wing with embedded piezoceramics is presented for energy harvesting. The electromechanical finite-element plate model is based on the thin-plate (Kirchhoff) assumptions while the unsteady aerodynamic model uses the doublet-lattice method. The electromechanical and aerodynamic models are combined to obtain the piezoaeroelastic equations, which are solved using a p-k scheme that accounts for the electromechanical coupling. The evolution of the aerodynamic damping and the frequency of each mode are obtained with changing airflow speed for a given electrical circuit. Expressions for piezoaeroelastically coupled frequency response functions (voltage, current, and electrical power as well the vibratory motion) are also defined by combining flow excitation with harmonic base excitation. Hence, piezoaeroelastic evolution can be investigated in frequency domain for different airflow speeds and electrical boundary conditions. [DOI:10.1115/1.4002785]

A 2-D Delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement.

Post-processing and visualization techniques in 2D boundary element analysis

Most post-processors for boundary element (BE) analysis use an auxiliary domain mesh to display domain results, working against the profitable modelling process of a pure boundary discretization. This paper introduces a novel visualization technique which preserves the basic properties of the boundary element methods. The proposed algorithm does not require any domain discretization and is based on the direct and automatic identification of isolines. Another critical aspect of the visualization of domain results in BE analysis is the effort required to evaluate results in interior points. In order to tackle this issue, the present article also provides a comparison between the performance of two different BE formulations (conventional and hybrid). In addition, this paper presents an overview of the most common post-processing and visualization techniques in BE analysis, such as the classical algorithms of scan line and the interpolation over a domain discretization. The results presented herein show that the proposed algorithm offers a very high performance compared with other visualization procedures.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Structure and properties of low temperature plasma carburized austenitic stainless steels

Austenitic stainless steels cannot be conventionally surface treated at temperatures close to 550 degrees C due to intense precipitation of nitrides or carbides. Plasma carburizing allows introducing carbon in the steel at temperatures below 500 degrees C without carbide precipitation. Plasma carburizing of AISI 316L was carried out at 480 degrees C and 400 degrees C, during 20 h, using CH(4) as carbon carrier gas. The results show that carbon expanded austenite (gamma(c)), 20 mu m in depth, was formed on the surface after the 480 degrees C treatment. Carbon expanded austenite (gamma(c)), 8 mu m in depth, was formed on the surface after the 400 degrees C treatment. DRX results showed that the austenitic FCC lattice parameter increases from 0.358 nm to 0.363 nm for the 400 degrees C treatment and to 0.369 nm for the 480 degrees C treatment, giving an estimation of circa 10 at.% carbon content for the latter. Lattice distortion, resulting from the expansion and the associated compressive residual stresses increases the surface hardness to 1040 HV(0.025). Micro-scale tensile tests were conducted on specimens prepared with the conditions selected above, which has indicated that the damage imposed to the expanded austenite layer was more easily related to each separated grain than to the overall macro-scale stresses imposed by the tensile test. (C) 2009 Elsevier B.V. All rights reserved.

Topology Optimized Design, Microfabrication and Characterization of Electro-Thermally Driven Microgripper

This article presents a systematic and logical study of the topology optimized design, microfabrication, and static/dynamic performance characterization of an electro-thermo-mechanical microgripper. The microgripper is designed using a topology optimization algorithm based on a spatial filtering technique and considering different penalization coefficients for different material properties during the optimization cycle. The microgripper design has a symmetric monolithic 2D structure which consists of a complex combination of rigid links integrating both the actuating and gripping mechanisms. The numerical simulation is performed by studying the effects of convective heat transfer, thermal boundary conditions at the fixed anchors, and microgripper performance considering temperature-dependent and independent material properties. The microgripper is fabricated from a 25 mm thick nickel foil using laser microfabrication technology and its static/dynamic performance is experimentally evaluated. The static and dynamic electro-mechanical characteristics are analyzed as step response functions with respect to tweezing/actuating displacements, applied current/power, and actual electric resistance. A microgripper prototype having overall dimensions of 1mm (L) X 2.5mm (W) is able to deliver the maximum tweezing and actuating displacements of 25.5 mm and 33.2 mm along X and Y axes, respectively, under an applied power of 2.32 W. Experimental performance is compared with finite element modeling simulation results.

Modelling grain boundary migration during geometric dynamic recrystallization

High-angle grain boundary migration is predicted during geometric dynamic recrystallization (GDRX) by two types of mathematical models. Both models consider the driving pressure due to curvature and a sinusoidal driving pressure owing to subgrain walls connected to the grain boundary. One model is based on the finite difference solution of a kinetic equation, and the other, on a numerical technique in which the boundary is subdivided into linear segments. The models show that an initially flat boundary becomes serrated, with the peak and valley migrating into both adjacent grains, as observed during GDRX. When the sinusoidal driving pressure amplitude is smaller than 2 pi, the boundary stops migrating, reaching an equilibrium shape. Otherwise, when the amplitude is larger than 2 pi, equilibrium is never reached and the boundary migrates indefinitely, which would cause the protrusions of two serrated parallel boundaries to impinge on each other, creating smaller equiaxed grains.

Hardness and structure characterization of Ti(6)Al(4)V films produced by reactive magnetron sputtering on a conventional austenitic stainless steel

Ti(6)Al(4)V thin films were grown by magnetron sputtering on a conventional austenitic stainless steel. Five deposition conditions varying both the deposition chamber pressure and the plasma power were studied. Highly textured thin films were obtained, their crystallite size (C) 2008 Elsevier Ltd. All rights reserved.

A stable MPC with zone control

Several MPC applications implement a control strategy in which some of the system outputs are controlled within specified ranges or zones, rather than at fixed set points [J.M. Maciejowski, Predictive Control with Constraints, Prentice Hall, New Jersey, 2002]. This means that these outputs will be treated as controlled variables only when the predicted future values lie outside the boundary of their corresponding zones. The zone control is usually implemented by selecting an appropriate weighting matrix for the output error in the control cost function. When an output prediction is inside its zone, the corresponding weight is zeroed, so that the controller ignores this output. When the output prediction lies outside the zone, the error weight is made equal to a specified value and the distance between the output prediction and the boundary of the zone is minimized. The main problem of this approach, as long as stability of the closed loop is concerned, is that each time an output is switched from the status of non-controlled to the status of controlled, or vice versa, a different linear controller is activated. Thus, throughout the continuous operation of the process, the control system keeps switching from one controller to another. Even if a stabilizing control law is developed for each of the control configurations, switching among stable controllers not necessarily produces a stable closed loop system. Here, a stable M PC is developed for the zone control of open-loop stable systems. Focusing on the practical application of the proposed controller, it is assumed that in the control structure of the process system there is an upper optimization layer that defines optimal targets to the system inputs. The performance of the proposed strategy is illustrated by simulation of a subsystem of an industrial FCC system. (C) 2008 Elsevier Ltd. All rights reserved.

Porous silicon optical cavity structure applied to high sensitivity organic solvent sensor

The present work reports the thermal annealing process, the number of layer and electrochemical process effect in the optical response quality of Bragg and microcavity devices that were applied as organic solvent sensors. These devices have been obtained by using porous silicon (PS) technology. The optical characterization of the Bragg reflector, before annealing, showed a broad photonic band-gap structure with blue shifted and narrowed after annealing process. The electrochemical process used to obtain the PS-based device imposes the limit in the number of layers because of the chemical dissolution effect. The interface roughness minimizations in the devices have been achieved by using the double electrochemical cell setup. The microcavity devices showed to have a good sensibility for organic solvent detection. The thermal annealed device showed better sensibility feature and this result was attributed to passivation of the surface devices. (c) 2007 Elsevier Ltd. All rights reserved.

Optimized network structure for full-synchronization

A network of Kuramoto oscillators with different natural frequencies is optimized for enhanced synchronizability. All node inputs are normalized by the node connectivity and some important properties of the network Structure are determined in this case: (i) optimized networks present a strong anti-correlation between natural frequencies of adjacent nodes: (ii) this anti-correlation should be as high as possible since the average path length between nodes is maintained as small as in random networks: and (iii) high anti-correlation is obtained without any relation between nodes natural frequencies and the degree of connectivity. We also propose a network construction model with which it is shown that high anti-correlation and small average paths may be achieved by randomly rewiring a fraction of the links of a totally anti-correlated network, and that these networks present optimal synchronization properties. (C) 2008 Elsevier B.V. All rights reserved.