Free vibration and wave propagation analysis of uniform and tapered rotating beams using spectrally formulated finite elements

This paper presents a formulation of an approximate spectral element for uniform and tapered rotating Euler-Bernoulli beams. The formulation takes into account the varying centrifugal force, mass and bending stiffness. The dynamic stiffness matrix is constructed using the weak form of the governing differential equation in the frequency domain, where two different interpolating functions for the transverse displacement are used for the element formulation. Both free vibration and wave propagation analysis is performed using the formulated elements. The studies show that the formulated element predicts results, that compare well with the solution available in the literature, at a fraction of the computational effort. In addition, for wave propagation analysis, the element shows superior convergence. (C) 2007 Elsevier Ltd. All rights reserved.

Fractional damping: Stochastic origin and finite approximations

Fractional-order derivatives appear in various engineering applications including models for viscoelastic damping. Damping behavior of materials, if modeled using linear, constant coefficient differential equations, cannot include the long memory that fractional-order derivatives require. However, sufficiently great rnicrostructural disorder can lead, statistically, to macroscopic behavior well approximated by fractional order derivatives. The idea has appeared in the physics literature, but may interest an engineering audience. This idea in turn leads to an infinite-dimensional system without memory; a routine Galerkin projection on that infinite-dimensional system leads to a finite dimensional system of ordinary differential equations (ODEs) (integer order) that matches the fractional-order behavior over user-specifiable, but finite, frequency ranges. For extreme frequencies (small or large), the approximation is poor. This is unavoidable, and users interested in such extremes or in the fundamental aspects of true fractional derivatives must take note of it. However, mismatch in extreme frequencies outside the range of interest for a particular model of a real material may have little engineering impact.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

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.

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.

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.

The Line Spectral Frequency Model of a Finite-Length Sequence

The line spectral frequency (LSF) of a causal finite length sequence is a frequency at which the spectrum of the sequence annihilates or the magnitude spectrum has a spectral null. A causal finite-length sequencewith (L + 1) samples having exactly L-LSFs, is referred as an Annihilating (AH) sequence. Using some spectral properties of finite-length sequences, and some model parameters, we develop spectral decomposition structures, which are used to translate any finite-length sequence to an equivalent set of AH-sequences defined by LSFs and some complex constants. This alternate representation format of any finite-length sequence is referred as its LSF-Model. For a finite-length sequence, one can obtain multiple LSF-Models by varying the model parameters. The LSF-Model, in time domain can be used to synthesize any arbitrary causal finite-length sequence in terms of its characteristic AH-sequences. In the frequency domain, the LSF-Model can be used to obtain the spectral samples of the sequence as a linear combination of spectra of its characteristic AH-sequences. We also summarize the utility of the LSF-Model in practical discrete signal processing systems.

Comparison of elasticity and shell-theory solutions for finite circular cylindrical shells

A finite circular cylindrical shell subjected to a band of uniform pressure on its outer rim was investigated, using three-dimensional elasticity theory and the classical shell theories of Timoshenko (or Donnell) and Flügge. Detailed comparison of the resulting stresses and displacements was carried out for shells with ratios of inner to outer shell radii equal to 0.80, 0.85, 0.90 and 0.93 and for ratios of outer shell diameter to length of the shell equal to 0.5, 1 and 2. The ratio of band width to length of the shell was 0.2 and Poisson's ratio used was equal to 0.3. An Elliot 803 digital computer was used for numerical computations.

Influence of finite size effect on magnetic and magnetotransport properties of La0.5Sr0.5CoO3 thin films

We present a comprehensive study of the thickness dependent structural, magnetic and magnetotransport properties of oriented La0.5Sr0.5CoO3 thin films grown on LaAlO3 by Pulsed Laser Deposition. We observe that these films undergo a reduction in Curie temperature (T-c) with a decrease in film thickness, and it is found to be primarily caused by the finite size effect since the finite scaling law [T-c(infinity) T-c(t)/T-c(infinity) = (c/t)lambda holds good over the studied thickness range. We rule out the contribution from the strain induced suppression of Curie temperature with decreasing film thickness since all the films exhibit a constant out of plane tensile strain (0.5%) irrespective of their varying thickness. However, we observe that the coercivity of the films is an order of magnitude higher than that of the bulk due to the tensile strain. In addition, we also observe an increase in the magneto resistance peak and a decrease in coercivity and electrical resistivity with an increase in film thickness. (C) 2010 Elsevier Ltd. All rights reserved.

Three-Dimensional simulation of deformation fields underneath vickers indenter:Effects of power-law Plasticity

The effects of power-law plasticity (yield strength and strain hardening exponent) on the plastic strain distribution underneath a Vickers indenter was systematically investigated by recourse to three-dimensional finite element analysis, motivated by the experimental macro-and micro-indentation on heat-treated Al-Zn-Mg alloy. For meaningful comparison between simulated and experimental results, the experimental heat treatment was carefully designed such that Al alloy achieve similar yield strength with different strain hardening exponent, and vice versa. On the other hand, full 3D simulation of Vickers indentation was conducted to capture subsurface strain distribution. Subtle differences and similarities were discussed based on the strain field shape, size and magnitude for the isolated effect of yield strength and strain hardening exponent.