New rational interpolation functions for finite element analysis of rotating beams

A rotating beam finite element in which the interpolating shape functions are obtained by satisfying the governing static homogenous differential equation of Euler–Bernoulli rotating beams is developed in this work. The shape functions turn out to be rational functions which also depend on rotation speed and element position along the beam and account for the centrifugal stiffening effect. These rational functions yield the Hermite cubic when rotation speed becomes zero. The new element is applied for static and dynamic analysis of rotating beams. In the static case, a cantilever beam having a tip load is considered, with a radially varying axial force. It is found that this new element gives a very good approximation of the tip deflection to the analytical series solution value, as compared to the classical finite element given by the Hermite cubic shape functions. In the dynamic analysis, the new element is applied for uniform, and tapered rotating beams with cantilever and hinged boundary conditions to determine the natural frequencies, and the results compare very well with the published results given in the literature.

The effects of manganese content and mould size on abrasion and slurry erosion behaviour of chromium-manganese iron systems investigated by positron lifetime spectroscopy

Abrasion and slurry erosion behaviour of chromium-manganese iron samples with chromium (Cr) in the range similar to 16-19% and manganese (Mn) at 5 and 10% levels have been characterized for hardness followed by microstructural examination using optical and scanning electron microscopy. Positron lifetime studies have been conducted to understand the defects/microporosity influence on the microstructure. The samples were heat treated and characterized to understand the structural transformations in the matrix. The data reveals that hardness decreased with increase in Mn content from 5 to 10% in the first instance and then increase in the section size in the other case, irrespective of the sample conditions. The abrasion and slurry erosion losses show increase with increase in the section size as well as with increase in Mn content. The positron results show that as hardness increases from as-cast to heat treated sample, the positron trapping rate and hence defect concentration showed opposite trend as expected. So a good correlation between defects concentration and the hardness has been observed. These findings also corroborate well with the microstructural features obtained from optical and scanning electron microscopy. (C) 2009 Elsevier B. V. All rights reserved.

Fracture process zone size and true fracture energy of concrete using acoustic emission

Acoustic emission (AE) energy, instead of amplitude, associated with each of the event is used to estimate the fracture process zone (FPZ) size. A steep increase in the cumulative AE energy of the events with respect to time is correlated with the formation of FPZ. Based on the AE energy released during these events and the locations of the events, FPZ size is obtained. The size-independent fracture energy is computed using the expressions given in the boundary effect model by least squares method since over-determined system of equations are obtained when data from several specimens are used. Instead of least squares method a different method is suggested in which the transition ligament length, measured from the plot of histograms of AE events plotted over the un-cracked ligament, is used directly to obtain size-independent fracture energy. The fracture energy thus calculated seems to be size-independent.

Parametric finite element analyses of geocellsupported embankments

This paper presents the results from parametric finite element analyses of geocell-supported embankments constructed on weak foundation soils. A composite model is used to numerically simulate the improvement in the strength and stiffness of the soil as a result of geocell confinement. The shear strength of the geocell-encased soil is obtained as a function of the additional confining pressure due to the geocell encasement considering it as a thin cylinder subjected to internal pressure. The stiffness of the geocell-encased soil is obtained from the stiffness of the unreinforced soil and the tensile modulus of the geocell material using an empirical equation. The validity of the model is verified by simulating the laboratory experiments on model geocell-supported embankments. Parametric finite element analyses of the geocell-supported embankments are carried out by varying the dimensions of the geocell layer, the tensile strength of the material used for fabricating the geocell layer, the properties of the infill soil, and the depth of the foundation layer. Some important guidelines for selecting the geocell reinforcement to support embankments on weak foundation soils are established through these numerical studies.

Parametrized actor-critic algorithms for finite-horizon MDPs

Due to their non-stationarity, finite-horizon Markov decision processes (FH-MDPs) have one probability transition matrix per stage. Thus the curse of dimensionality affects FH-MDPs more severely than infinite-horizon MDPs. We propose two parametrized 'actor-critic' algorithms to compute optimal policies for FH-MDPs. Both algorithms use the two-timescale stochastic approximation technique, thus simultaneously performing gradient search in the parametrized policy space (the 'actor') on a slower timescale and learning the policy gradient (the 'critic') via a faster recursion. This is in contrast to methods where critic recursions learn the cost-to-go proper. We show w.p 1 convergence to a set with the necessary condition for constrained optima. The proposed parameterization is for FHMDPs with compact action sets, although certain exceptions can be handled. Further, a third algorithm for stochastic control of stopping time processes is presented. We explain why current policy evaluation methods do not work as critic to the proposed actor recursion. Simulation results from flow-control in communication networks attest to the performance advantages of all three algorithms.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Context-dependency of a complex fruit-frugivore mutualism: temporal variation in crop size and neighborhood effects

The quantity of fruit consumed by dispersers is highly variable among individuals within plant populations. The outcome Of Such selection operated by firugivores has been examined mostly with respect to changing spatial contexts. The influence of varying temporal contexts on frugivore choice, and their possible demographic and evolutionary consequences is poorly understood. We examined if temporal variation in fruit availability across a hierarchy of nested temporal levels (interannual, intraseasonal, 120 h, 24 h) altered frugivore choice for a complex seed dispersal system in dry tropical forests of southern India. The interactions between Phyllanthus emblica and its primary disperser (ruminants) was mediated by another frugivore (a primate),which made large quantities of fruit available on the ground to ruminants. The direction and strength of crop size and neighborhood effects on this interaction varied with changing temporal contexts.Fruit availability was higher in the first of the two study years, and at the start of the season in both years. Fruit persistence on trees,determined by primate foraging, was influenced by crop size andconspecific neighborhood densities only in the high fruit availability year. Fruit removal by ruminants was influenced by crop size in both years and neighborhood densities only in the high availability year. In both years, these effects were stronger at the start of the season.Intraseasonal reduction in fruit availability diminished inequalities in fruit removal by ruminants and the influence of crop size and fruiting neighborhoods. All trees were not equally attractive to frugivores in a P. emblica population at all points of time. Temporal asymmetry in frugivore-mediated selection could reduce potential for co-evolution between firugivores and plants by diluting selective pressures. Inter-dependencies; formed between disparate animal consumers can add additional levels of complexity to plant-frugivore mutualistic networks and have potential reproductive consequences for specific individuals within populations.

Size-dependent magnetic properties of iron carbide nanoparticles embedded in a carbon matrix

Here we report on the magnetic properties of iron carbide nanoparticles embedded in a carbon matrix. Granular distributions of nanoparticles in an inert matrix, of potential use in various applications, were prepared by pyrolysis of organic precursors using the thermally assisted chemical vapour deposition method. By varying the precursor concentration and preparation temperature, compositions with varying iron concentration and nanoparticle sizes were made. Powder x-ray diffraction, transmission electron microscopy and Mossbauer spectroscopy studies revealed the nanocrystalline iron carbide (Fe3C) presence in the partially graphitized matrix. The dependence of the magnetic properties on the particle size and temperature (10 K < T < 300 K) were studied using superconducting quantum interference device magnetometry. Based on the affect of surrounding carbon spins, the observed magnetic behaviour of the nanoparticle compositions, such as the temperature dependence of magnetization and coercivity, can be explained.

Investigations of the scaling criteria for a mild combustion burner

In this paper, a new strategy for scaling burners based on "mild combustion" is evolved and adopted to scaling a burner from 3 to a 150 kW burner at a high heat release Late of 5 MW/m(3) Existing scaling methods (constant velocity, constant residence time, and Cole's procedure [Proc. Combust. Inst. 28 (2000) 1297]) are found to be inadequate for mild combustion burners. Constant velocity approach leads to reduced heat release rates at large sizes and constant residence time approach in unacceptable levels of pressure drop across the system. To achieve mild combustion at high heat release rates at all scales, a modified approach with high recirculation is adopted in the present studies. Major geometrical dimensions are scaled as D similar to Q(1/3) with an air injection velocity of similar to 100 m/s (Delta p similar to 600 mm water gauge). Using CFD support, the position of air injection holes is selected to enhance the recirculation rates. The precise role of secondary air is to increase the recirculation rates and burn LIP the residual CO in the downstream. Measurements of temperature and oxidizer concentrations inside 3 kW, 150 kW burner and a jet flame are used to distinguish the combustion process in these burners. The burner can be used for a wide range of fuels from LPG to producer gas as extremes. Up to 8 dB of noise level reduction is observed in comparison to the conventional combustion mode. Exhaust NO emissions below 26 and 3 ppm and temperatures 1710 and 1520 K were measured for LPG and producer gas when the burner is operated at stoichiometry. (c) 2004 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Study of gas cavity size hysteresis in a packed bed using DEM

When a high velocity gas jet is introduced into a packed bed a cavity is formed. The size of the cavity shows hysteresis on increasing and decreasing gas flow rates. This hysteresis leads to different cavity sizes at same gas flow rate depending on the bed history. The size of cavity affects the gas flow profiles in the packed bed. In this study the cavity size hysteresis phenomenon has been modeled using discrete element method along with turbulent gas flow. A reasonable agreement has been found between computed and experimental results on cavity size ysteresis. The effect of various parameters, such as nozzle height from the bed bottom and packing height, on the cavity size hysteresis has been studied. It is found that inter-particle interaction forces along with gas drag and bed porosity play an important role in describing the cavity size hysteresis. The injection of gas flow allows the particles to go to an unconstrained state than they were previously in, and their ability to remain in that state, even under decreased gas drag force, leads to the phenomenon of cavity size hysteresis. (c) 2007 Elsevier Ltd. All rights reserved.

The dynamics of solvation of an electron in the image potential state by a layer of polar adsorbates

Recently, ultrafast two-photon photoemission has been used to study electron solvation at a two-dimensional metal/polar adsorbate interfaces [A. Miller , Science 297, 1163 (2002)]. The electron is bound to the surface by the image interaction. Earlier we have suggested a theoretical description of the states of the electron interacting with a two-dimensional layer of the polar adsorbate [K. L. Sebastian , J. Chem. Phys. 119, 10350 (2003)]. In this paper we have analyzed the dynamics of electron solvation, assuming a trial wave function for the electron and the solvent polarization and then using the Dirac-Frenkel variational method to determine it. The electron is initially photoexcited to a delocalized state, which has a finite but large size, and causes the polar molecules to reorient. This reorientation acts back on the electron and causes its wave function to shrink, which will cause further reorientation of the polar molecules, and the process continues until the electron gets self-trapped. For reasonable values for the parameters, we are able to obtain fair agreement with the experimental observations. (c) 2005 American Institute of Physics.

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Unusual dielectric response in B-site size-disordered hexagonal transition metal oxides

We discover that hexagonal holmium copper titanate (Ho2CuTiO6), has a unique and highly desirable combination of high dielectric constant, low losses, very small temperature coefficient, and low frequency dependence. Our first-principles calculations indicate that these exceptional properties result from a size-difference at the Cu/Ti B-site that suppresses the expected ferroelectric transition, combined with the dominance of intermediate-frequency polar vibrational modes in the dielectric response. Our results suggest that the use of such B-site disorder in alloys of hexagonal transition-metal oxides should generally result in similar robust dielectrics.

Thermal stresses in a finite hollow cylinder due to an axisymmetric temperature field at the end surface

A three-dimensional exact solution for determining the thermal stresses in a finite hollow cylinder subject to a steady state axisymmetric temperature field over one of its end surfaces has been given. Numerical results for a hollow cylinder, having lenght to outer diameter ratio equal to one and inner to outer diameter ratio equal to 0.75, subjected to a symmetric temperature variation over the end surfaces of the cylinder have been given.

Thermal stresses in a finite solid cylinder due to an axisymmetric temperature field at the end surface.

A three-dimensional rigorous solution for determining thermal stresses in a finite solid cylinder due to a steady state axisymmetric temperature field over one of its end surfaces is given. Numerical results for a solid cylinder having a length to diameter ratio equal to one and subjected to a symmetric temperature variation over half the radius of the cylinder at the end surfaces are included. These results have been compared with the results of the approximate solution given by W. Nowacki.