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.

Numerical studies on channel formation and growth during solidification: Effect of process parameters

In the present work, solidification of a hyper-eutectic ammonium chloride solution in a bottom-cooled cavity (i.e. with stable thermal gradient) is numerically studied. A Rayleigh number based criterion is developed, which determines the conditions favorable for freckles formation. This criterion, when expressed in terms of physical properties and process parameters, yields the condition for plume formation as a function of concentration, liquid fraction, permeability, growth rate of a mushy layer and thermophysical properties. Subsequently, numerical simulations are performed for cases with initial and boundary conditions favoring freckle formation. The effects of parameters, such as cooling rate and initial concentration, on the formation and growth of freckles are investigated. It was found that a high cooling rate produced larger and more defined channels which are retained for a longer durations. Similarly, a lower initial concentration of solute resulted in fewer but more pronounced channels. The number and size of channels are also found to be related to the mushy zone thickness. The trends predicted with regard to the variation of number of channels with time under different process conditions are in accordance with the experimental observations reported in the literature.

Self-interrupted regenerative metal cutting in turning

A new approach is used to study the global dynamics of regenerative metal cutting in turning. The cut surface is modeled using a partial differential equation (PDE) coupled, via boundary conditions, to an ordinary differential equation (ODE) modeling the dynamics of the cutting tool. This approach automatically incorporates the multiple-regenerative effects accompanying self-interrupted cutting. Taylor's 3/4 power law model for the cutting force is adopted. Lower dimensional ODE approximations are obtained for the combined tool–workpiece model using Galerkin projections, and a bifurcation diagram computed. The unstable solution branch off the subcritical Hopf bifurcation meets the stable branch involving self-interrupted dynamics in a turning point bifurcation. The tool displacement at that turning point is estimated, which helps identify cutting parameter ranges where loss of stability leads to much larger self-interrupted motions than in some other ranges. Numerical bounds are also obtained on the parameter values which guarantee global stability of steady-state cutting, i.e., parameter values for which there exist neither unstable periodic motions nor self-interrupted motions about the stable equilibrium.

Secondary flows induced by the rotation of a sphere or coaxial cones in micropolar fluids

The micropolar fluids like Newtonian and Non-Newtonian fluids cannot sustain a simple shearing motion, wherein only one component of velocity is present. They exhibit both primary and secondary motions when the boundaries are subject to slow rotations. The primary motion, as in Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., resembles that of Newtonian fluid for slow steady rotation. We further notice that the micro-rotation becomes identically equal to the vorticity present in the fluid and the condition b) of "Wall vorticity" can alone be satisfied at the boundaries. As regards, the secondary motion, we notice that it can be determined by the above procedure for a special class of fluids, namely that for which j0(n2-n3)=4 n3/l2. Moreover for this class of fluids, the micro-rotation is identical with the vorticity of the fluid everywhere. Also the stream function for the secondary flow is identical with that for the Newtonian fluid with a suitable definition of the Reynolds number. In contrast with the Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., this class of micropolar fluids does not show separation. This is in conformity with the statement of Condiff and Dahler (3) that in any steady flow, internal spin matches the vorticity everywhere provided that (i) spin boundary conditions are satisfied, (ii) body torques and non-conservative body forces are absent, and (iii) inertial and spin-inertial terms are either negligible or vanish identically.

Four generations: SUSY and SUSY breaking

We revisit four generations within the context of supersymmetry. Wecompute the perturbativity limits for the fourth generation Yukawa couplings and show that if the masses of the fourth generation lie within reasonable limits of their present experimental lower bounds, it is possible to have perturbativity only up to scales around 1000 TeV. Such low scales are ideally suited to incorporate gauge mediated supersymmetry breaking, where the mediation scale can be as low as 10-20 TeV. The minimal messenger model, however, is highly constrained. While lack of electroweak symmetry breaking rules out a large part of the parameter space, a small region exists, where the fourth generation stau is tachyonic. General gauge mediation with its broader set of boundary conditions is better suited to accommodate the fourth generation.

Large deflections of rectangular plates

Approximate solutions for the non-linear bending of thin rectangular plates are presented considering large deflections for various boundary conditions. In the case of stress-free edges, solutions are given for von Kármán's equations in terms of the stress function and the deflection of the plate. In the case of immovable edges, equations are constructed in terms of the three displacements and these are solved. The solution is given by using double series consisting of the appropriate Beam Functions which satisfy the boundary conditions. The differential equations are satisfied by using the orthogonality properties of the series. Numerical results for square plates with uniform lateral load indicate good convergence of the series solution presented here.

Thermal stresses in rectangular plates

A method of determining the thermal stresses in a flat rectangular isotropic plate of constant thickness with arbitrary temperature distribution in the plane of the plate and with no variation in temperature through the thickness is presented. The thermal stress have been obtained in terms of Fourier series and integrals that satisfy the differential equation and the boundary conditions. Several examples have been presented to show the application of the method.

Vibration of rectangular orthotropic plates

An approximate analytical procedure has been given to solve the problem of a vibrating rectangular orthotropic plate, with various combinations of simply supported and clamped boundary conditions. Numerical results have been given for the case of a clamped square plate. Nomenclature 2a, 2b sides of the rectangular plate h plate thickness Eprime x , Eprime y , EPrime, G elastic constants of te orthotropic material D x Eprime x h 3/12 D y Eprime y h 3/12 H xy EPrimeh 3/12+Gh 3/6 D x , D y and H xy are rigidity constants of the orthotropic platergr mass of the plate per unit area ngr Poisson's ratio W deflection of the plate p circular frequency gamma b/a ratio X m , Y characteristic functions of the vibrating beam problem -lambda rgrp 2 a 2 b 2/H xy the frequency parameter.

Flow of non-newtonian fluid between torsionally oscillating discs

The flow of an incompressible non-Newtonian viscous fluid contained between two torsionally oscillating infinite parallel discs is investigated. The two specific cases studied are (i) one disc only oscillates while the other is at rest and (ii) both discs oscillate with the same frequency and amplitude but in opposite directions. Assuming that the amplitude of oscillation,Ω/n, is small and neglecting the squares and higher powers ofΩ/n, the equations of motion have been solved exactly for velocity and pressure satisfying all the boundary conditions. The effect of both positive and negative coefficients of cross-viscosity on the steady components of the flow has been represented graphically.

Effect of matrix cracking and material uncertainty on composite plates

A laminated composite plate model based on first order shear deformation theory is implemented using the finite element method.Matrix cracks are introduced into the finite element model by considering changes in the A, B and D matrices of composites. The effects of different boundary conditions, laminate types and ply angles on the behavior of composite plates with matrix cracks are studied.Finally, the effect of material property uncertainty, which is important for composite material on the composite plate, is investigated using Monte Carlo simulations. Probabilistic estimates of damage detection reliability in composite plates are made for static and dynamic measurements. It is found that the effect of uncertainty must be considered for accurate damage detection in composite structures. The estimates of variance obtained for observable system properties due to uncertainty can be used for developing more robust damage detection algorithms. (C) 2010 Elsevier Ltd. All rights reserved.

The transient response of a two-terminal solion to a step-current input

The transition time associated with the time-variation of the voltage across a two-terminal diaphragm-less solion in response to a step-current stimulus has been studied experimentally. A theoretical analysis has also been made by solving the diffusion problem under the appropriate initial and boundary conditions. The behaviour of the theoretically predicted transition times is in agreement with the observed behaviour. The systems under study have been shown to be different from those used hitherto in thin-layer chronopotentiometry.

The transient response of a two-terminal solion to a stepcrrent input

The transition time associated with the time-variation of the voltage across a two-terminal diaphragm-less solion in response to a step-current stimulus has been studied experimentally. A theoretical analysis has also been made by solving the diffusion problem under the appropriate initial and boundary conditions. The behaviour of the theoretically predicted transition times is in agreement with the observed behaviour. The systems under study have been shown to be different from those used hitherto in thin-layer chronopotentiometry.

Thermodynamic consistency of the ternary functions

A ternary thermodynamic function has been developed based on statistico-thermodynamic considerations, with a particular emphasis on the higher-order terms indicating the effects of truncation at the various stages of the treatment. Although the truncation of a series involved in the equation introduces inconsistency, the latter may be removed by imposing various thermodynamic boundary conditions. These conditions are discussed in the paper. The present equation with higher-order terms shows that the α function of a component reduces to a quadratic function of composition at constant compositional paths involving the other two components in the system. The form of the function has been found to be representative of various experimental observations.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Elastic property estimation using Ultrasound Assisted Optical Elastography through Remote Palpation - A simulation study - art. no. 61640Q

We propose an effective elastography technique in which an acoustic radiation force is used for remote palpation to generate localized tissue displacements, which are directly correlated to localized variations of tissue stiffness and are measured using a light probe in the same direction of ultrasound propagation. The experimental geometry has provision to input light beam along the ultrasound propagation direction, and hence it can be prealigned to ensure proper interception of the focal region by the light beam. Tissue-mimicking phantoms with homogeneous and isotropic mechanical properties of normal and malignant breast tissue are considered for the study. Each phantom is insonified by a focusing ultrasound transducer (1 MHz). The focal volume of the transducer and the ultrasound radiation force in the region are estimated through solving acoustic wave propagation through medium assuming average acoustic properties. The forward elastography problem is solved for the region of insonification assuming the Lame's parameters and Poisson's ratio, under Dirichlet boundary conditions which gives a distribution of displacement vectors. The direction of displacement, though presented spatial variation, is predominantly towards the ultrasound propagation direction. Using Monte Carlo (MC) simulation we have traced the photons through the phantom and collected the photons arriving at the detector on the boundary of the object in the direction of ultrasound. The intensity correlations are then computed from detected photons. The intensity correlation function computed through MC simulation showed a modulation whose strength is found to be proportional to the amplitude of displacement and inversely related to the storage (elastic) modulus. It is observed that when the storage modulus in the focal region is increased the computed displacement magnitude, as indicated by the depth of modulation in the intensity autocorrelation, decreased and the trend is approximately exponential.