Coupled polynomial field approach for elimination of flexure and torsion locking phenomena in the Timoshenko and Euler-Bernoulli curved beam elements

The curvature related locking phenomena in the out-of-plane deformation of Timoshenko and Euler-Bernoulli curved beam elements are demonstrated and a novel approach is proposed to circumvent them. Both flexure and Torsion locking phenomena are noticed in Timoshenko beam and torsion locking phenomenon alone in Euler-Bernoulli beam. Two locking-free curved beam finite element models are developed using coupled polynomial displacement field interpolations to eliminate these locking effects. The coupled polynomial interpolation fields are derived independently for Timoshenko and Euler-Bernoulli beam elements using the governing equations. The presented of penalty terms in the couple displacement fields incorporates the flexure-torsion coupling and flexure-shear coupling effects in an approximate manner and produce no spurious constraints in the extreme geometric limits of flexure, torsion and shear stiffness. the proposed couple polynomial finite element models, as special cases, reduce to the conventional Timoshenko beam element and Euler-Bernoulli beam element, respectively. These models are shown to perform consistently over a wide range of flexure-to-shear (EI/GA) and flexure-to-torsion (EI/GJ) stiffness ratios and are inherently devoid of flexure, torsion and shear locking phenomena. The efficacy, accuracy and reliability of the proposed models to straight and curved beam applications are demonstrated through numerical examples. (C) 2012 Elsevier B.V. All rights reserved.

An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Determination of Isospectral Nonuniform Rotating Beams

In this paper we look for nonuniform rotating beams that are isospectral to a given uniform nonrotating beam. A rotating nonuniform beam is isospectral to the given uniform nonrotating beam if both the beams have the same spectral properties, i.e., both the beams have the same set of natural frequencies under a given boundary condition. The Barcilon-Gottlieb type transformation is proposed that converts the governing equation of a rotating beam to that of a uniform nonrotating beam. We show that there exist rotating beams isospectral to a given uniform nonrotating beam under some special conditions. The boundary conditions we consider are clamped-free and hinged-free with an elastic hinge spring. An upper bound on the rotation speed for which isospectral beams exist is proposed. The mass and stiffness distributions for these nonuniform rotating beams which are isospectral to the given uniform nonrotating beam are obtained. We use these mass and stiffness distributions in a finite element analysis to show that the obtained beams are isospectral to the given uniform nonrotating beam. A numerical example of a beam having a rectangular cross section is presented to show the application of our analysis. DOI: 10.1115/1.4006460]

Wave propagation in a porous composite beam: Porosity determination, location and quantification

A wave-based method is developed to quantify the defect due to porosity and also to locate the porous regions, in a composite beam-type structure. Wave propagation problem for a porous laminated composite beam is modeled using spectral finite element method (SFEM), based on the modified rule of mixture approach, which is used to include the effect of porosity on the stiffness and density of the composite beam structure. The material properties are obtained from the modified rule of mixture model, which are used in a conventional SFEM to develop a new model for solving wave propagation problems in porous laminated composite beam. The influence of the porosity content on the group speed and also the effect of variation in theses parameters on the time responses are studied first, in the forward problem. The change in the time responses with the change in the porosity of the structure is used as a parameter to find the porosity content in a composite beam. The actual measured response from a structure and the numerically obtained time responses are used for the estimation of porosity, by solving a nonlinear optimization problem. The effect of the length of the porous region (in the propagation direction), on the time responses, is studied. The damage force indicator technique is used to locate the porous region in a beam and also to find its length, using the measured wave propagation responses. (C) 2012 Elsevier Ltd. All rights reserved.

MACHINING OF MMCs: A REVIEW

Over the last few decades, Metal Matrix Composites (MMCs) have emerged as a material system offering tremendous potential for future applications. The primary advantages offered by these materials are their improved mechanical properties, particularly in the areas of wear, strength and stiffness. Of the MMCs, Aluminum matrix composites have grown in prominence due to their low density, low melting point and low cost. However, machining these materials remains a challenging task mainly due to the high abrasiveness of the reinforcing phases. Conventional machining processes such as turning, milling or drilling are adopted for machining MMCs. In this article, the existing and ongoing developments in machining MMCs vis-a-vis tool life, tool wear, machinability and understanding chip formation mechanism have been highlighted. Most of the studies discussed in this review will focus on Aluminum matrix composites. Certain areas of machining studies which have hitherto not been investigated have also been detailed.

Analogy Between Rotating Euler-Bernoulli and Timoshenko Beams and Stiff Strings

The governing differential equation of a rotating beam becomes the stiff-string equation if we assume uniform tension. We find the tension in the stiff string which yields the same frequency as a rotating cantilever beam with a prescribed rotating speed and identical uniform mass and stiffness. This tension varies for different modes and are found by solving a transcendental equation using bisection method. We also find the location along the rotating beam where equivalent constant tension for the stiff string acts for a given mode. Both Euler-Bernoulli and Timoshenko beams are considered for numerical results. The results provide physical insight into relation between rotating beams and stiff string which are useful for creating basis functions for approximate methods in vibration analysis of rotating beams.

Non-rotating beams isospectral to a given rotating uniform beam

In this paper, we seek to find non-rotating beams with continuous mass and flexural stiffness distributions, that are isospectral to a given uniform rotating beam. The Barcilon-Gottlieb transformation is used to convert the fourth order governing equation of a non-rotating beam, to a canonical fourth order eigenvalue problem. If the coefficients in this canonical equation match with the coefficients of the uniform rotating beam equation, then the non-rotating beam is isospectral to the given rotating beam. The conditions on matching the coefficients leads to a pair of coupled differential equations. We solve these coupled differential equations for a particular case, and thereby obtain a class of non-rotating beams that are isospectral to a uniform rotating beam. However, to obtain isospectral beams, the transformation must leave the boundary conditions invariant. We show that the clamped end boundary condition is always invariant, and for the free end boundary condition to be invariant, we impose certain conditions on the beam characteristics. We also verify numerically that the frequencies of the non-rotating beam obtained using the finite element method (FEM) are the exact frequencies of the uniform rotating beam. Finally, the example of beams having a rectangular cross-section is presented to show the application of our analysis. Since experimental determination of rotating beam frequencies is a difficult task, experiments can be easily conducted on these rectangular non-rotating beams, to calculate the frequencies of the rotating beam. (c) 2012 Elsevier Ltd. All rights reserved.

Rotating beams isospectral to axially loaded nonrotating uniform beams

In this paper, the stiffness and mass per unit length distributions of a rotating beam, which is isospectral to a given uniform axially loaded nonrotating beam, are determined analytically. The Barcilon-Gottlieb transformation is extended so that it transforms the governing equation of a rotating beam into the governing equation of a uniform, axially loaded nonrotating beam. Analysis is limited to a certain class of Euler-Bernoulli cantilever beams, where the product between the stiffness and the cube of mass per unit length is a constant. The derived mass and stiffness distributions of the rotating beam are used in a finite element analysis to confirm the frequency equivalence of the given and derived beams. Examples of physically realizable beams that have a rectangular cross section are shown as a practical application of the analysis.

Modeling the dynamic response of wrap-faced reinforced soil retaining walls

This paper describes the development of a numerical model for simulating the shaking table tests on wrap-faced reinforced soil retaining walls. Some of the physical model tests carried out on reinforced soil retaining walls subjected to dynamic excitation through uniaxial shaking tests are briefly discussed. Models of retaining walls are constructed in a perspex box with geotextile reinforcement using the wraparound technique with dry sand backfill and instrumented with displacement sensors, accelerometers, and soil pressure sensors. Results showed that the displacements decrease with the increase in number of reinforcement layers, whereas acceleration amplifications were not affected significantly. Numerical modeling of these shaking table tests is carried out using the Fast Lagrangian Analysis of Continua program. The numerical model is validated by comparing the results with experiments on physical models. Responses of wrap-faced walls with varying numbers of reinforcement layers are compared. Sensitivity analysis performed on the numerical models showed that the friction and dilation angle of backfill material and stiffness properties of the geotextile-soil interface are the most affecting parameters for the model response.

Some studies on evaluation of degradation in composite adhesive joints using ultrasonic techniques

Experimental and theoretical studies on degradation of composite-epoxy adhesive joints were carried out on samples having different interfacial and cohesive properties. Oblique incidence ultrasonic inspection of bonded joints revealed that degradation in the adhesive can be measured by significant variation in reflection amplitude as also by a shift in the minima of reflection spectrum. It was observed that severe degradation of the adhesive leads to failure dominated by interfacial mode. Through this investigation it is demonstrated that a correlation exists between the bond strength and a frequency shift in reflection minimum. The experimental data was validated using analytical models. Though both bulk adhesive degradation and interfacial degradation influences the shift in spectrum minimum, the contribution of the latter was found to be significant. An inversion algorithm was used to determine the interfacial transverse stiffness using the experimental oblique reflection spectrum. The spectrum shift was found to depend on the value of interfacial transverse stiffness using which a qualitative assessment can be made on the integrity of the joint.

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

Topology optimization-based design of a compliant aircraft compliant aircraft wing for morphing leading and trailing edges

In this paper, we present a methodology for designing a compliant aircraft wing, which can morph from a given airfoil shape to another given shape under the actuation of internal forces and can offer sufficient stiffness in both configurations under the respective aerodynamic loads. The least square error in displacements, Fourier descriptors, geometric moments, and moment invariants are studied to compare candidate shapes and to pose the optimization problem. Their relative merits and demerits are discussed in this paper. The `frame finite element ground structure' approach is used for topology optimization and the resulting solutions are converted to continuum solutions. The introduction of a notch-like feature is the key to the success of the design. It not only gives a good match for the target morphed shape for the leading and trailing edges but also minimizes the extension of the flexible skin that is to be put on the airfoil frame. Even though linear small-displacement elastic analysis is used in optimization, the obtained designs are analysed for large displacement behavior. The methodology developed here is not restricted to aircraft wings; it can be used to solve any shape-morphing requirement in flexible structures and compliant mechanisms.

An approach for assessing the accuracy of accelerometer responses during axial impact tests

The primary objective of the present study is to show that for the most common configuration of an impactor system, the accelerometer cannot exactly reproduce the dynamic response of a specimen subjected to impact loading. An equivalent Lumped Parameter Model (LPM) of the given impactor set-up has been formulated for assessing the accuracy of an accelerometer mounted in a drop-weight impactor set-up for an axially loaded specimen. A specimen under the impact loading is represented by a non-linear spring of varying stiffness, while the accelerometer is assumed to behave in a linear manner due to its high stiffness. Specimens made of steel, aluminium and fibre-reinforced composite (FRC) are used in the present study. Assuming the force-displacement response obtained in an actual impact test to be the true behaviour of the test specimen, a suitable numerical approach has been used to solve the governing non-linear differential equations of a three degrees-of-freedom (DOF) system in a piece-wise linear manner. The numerical solution of the governing differential equations following an explicit time integration scheme yields an excellent reproduction of the mechanical behaviour of the specimen, consequently confirming the accuracy of the numerical approach. However, the spring representing the accelerometer predicts a response that qualitatively matches the assumed force-displacement response of the test specimen with a perceptibly lower magnitude of load.

Substrate conductivity dependent modulation of cell proliferation and differentiation in vitro

In designing and developing various biomaterials, the influence of substrate properties, like surface topography, stiffness and wettability on the cell functionality has been investigated widely. However, such study to probe into the influence of the substrate conductivity on cell fate processes is rather limited. In order to address this issue, spark plasma sintered HA-CaTiO3 (Hydroxyapatite-Calcium titanate) has been used as a model material system to showcase the effect of varying conductivity on cell functionality. Being electroactive in nature, mouse myoblast cells (C2C12) were selected as a model cell line in this study. It was inferred that myoblast adhesion/growth systematically increases with substrate conductivity due to CaTiO3 addition to HA. Importantly, parallel arrangement of myoblast cells on higher CaTiO3 containing substrates indicate that self-adjustable cell patterning can be achieved on conductive biomaterials. Furthermore, enhanced myoblast assembly and myotube formation were recorded after 5 days of serum starvation. Overall, the present study conclusively establishes the positive impact of the substrate conductivity towards cell proliferation and differentiation as well as confirms the efficacy of HA-CaTiO3 biocomposites as conductive platforms to facilitate the growth, orientation and fusion of myoblasts, even when cultured in the absence of external electric field.

Rotating beams and non-rotating beams with shared eigenpair for pinned-free boundary condition

In this paper we look for a rotating beam, with pinned-free boundary conditions, whose eigenpair (frequency and mode-shape) is same as that of a uniform non-rotating beam for a particular mode. It is seen that for any given mode, there exists a flexural stiffness function (FSF) for which the ith mode eigenpair of a rotating beam with uniform mass distribution, is identical to that of a corresponding non-rotating beam with same length and mass distribution. Inserting these derived FSF's in a finite element code for a rotating pinned-free beam, the frequencies and mode shapes of a non-rotating pinned-free beam are obtained. For the first mode, a physically realistic equivalent rotating beam is possible, but for higher modes, the FSF has internal singularities. Strategies for addressing these singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test functions for rotating beam codes and also for targeted destiffening of rotating beams.