Free vibration of a stochastic beam-column using stochastic fem

The free vibrational characteristics of a beam-column, which is having randomly varying Young's modulus and mass density and subjected to randomly distributed axial loading is analysed. The material property fluctuations and axial loadings are considered to constitute independent one-dimensional, uni-variate, homogeneous real, spatially distributed stochastic fields. Hamilton's principle is used to formulate the problem using stochastic FEM. Vibration frequencies and mode shapes are analysed for their statistical descriptions. A numerical example is shown.

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Design, modelling and simulation of vibratory micromachined gyroscopes

Among various MEMS sensors, a rate gyroscope is one of the most complex sensors from the design point of view. The gyro normally consists of a proof mass suspended by an elaborate assembly of beams that allow the system to vibrate in two transverse modes. The structure is normally analysed and designed using commercial FEM packages such as ANSYS or MEMS specific commercial tools such as Coventor or Intellisuite. In either case, the complexity in analysis rises manyfolds when one considers the etch hole topography and the associated fluid flow calculation for damping. In most cases, the FEM analysis becomes prohibitive and one resorts to equivalent electrical circuit simulations using tools like SABER in Coventor. Here, we present a simplified lumped parameter model of the tuning fork gyro and show how easily it can be implemented using a generic tool like SIMULINK. The results obtained are compared with those obtained from more elaborate and intense simulations in Coventor. The comparison shows that lumped parameter SIMULINK model gives equally good results with fractional effort in modelling and computation. Next, the performance of a symmetric and decoupled vibratory gyroscope structure is also evaluated using this approach and a few modifications are made in this design to enhance the sensitivity of the device.

Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

Wrinkled and slack membranes: nonlinear 3D elasticity solutions via smooth DMS-FEM and experiment

The smooth DMS-FEM, recently proposed by the authors, is extended and applied to the geometrically nonlinear and ill-posed problem of a deformed and wrinkled/slack membrane. A key feature of this work is that three-dimensional nonlinear elasticity equations corresponding to linear momentum balance, without any dimensional reduction and the associated approximations, directly serve as the membrane governing equations. Domain discretization is performed with triangular prism elements and the higher order (C1 or more) interelement continuity of the shape functions ensures that the errors arising from possible jumps in the first derivatives of the conventional C0 shape functions do not propagate because the ill-conditioned tangent stiffness matrices are iteratively inverted. The present scheme employs no regularization and exhibits little sensitivity to h-refinement. Although the numerically computed deformed membrane profiles do show some sensitivity to initial imperfections (nonplanarity) in the membrane profile needed to initiate transverse deformations, the overall patterns of the wrinkles and the deformed shapes appear to be less so. Finally, the deformed profiles, computed through the DMS FEM-based weak formulation, are compared with those obtained through an experiment on an ultrathin Kapton membrane, wherein wrinkles form because of the applied boundary displacement conditions. Comparisons with a reported experiment on a rectangular membrane are also provided. These exercises lend credence to the feasibility of the DMS FEM-based numerical route to computing post-wrinkled membrane shapes. Copyright (c) 2012 John Wiley & Sons, Ltd.

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.

FEM Simulations and Optimization about Residual Stresses in Coating Structures with Functionally Graded Materials Layer

An elasto-plastic finite element method is developed to predict the residual stresses of thermal spraying coatings with functionally graded material layer. In numerical simulations, temperature sensitivity of various material constants is included and mix

Application of laser cladding technology to improve the interface of thermal barrier coatings

The present study is focused on improvement of the adhesion properties of the interface between plasma-sprayed coatings and substrates by laser cladding technology (LCT), Within the laser-clad layer there is a gradient distribution in chemical composition and mechanical properties that has been confirmed by SEM observation and microhardness measurement. The residual stress due to mismatches in thermal and mechanical properties between coatings and substrates can be markedly reduced and smoothed out. To examine the changes of microstructure and crack propagation in the coating and interface during loading, the three-point bending test has been carried out in SEM with a loading device. Analysis of the distribution of shear stress near the interface under loading has been made using the FEM code ANSYS, The experimental results show clearly that the interface adhesion can be improved with LCT pretreatment, and the capability of the interface to withstand the shear stress as well as to resist microcracking has been enhanced.

Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Characterising Of Internal Stresses In Duplex Coating By Fem

The internal stresses in a duplex coating involving a prequenched layer are believed to change if it is exposed to thermal loading. To characterise the internal stresses in such a duplex coating, a gradient model of finite element method is set up. The initial stress within the substrate developed in as quenching and the internal stresses due to the tempering of the prequenched layer ( TPQL) in such a duplex coating are calculated. The synthetical internal stresses in coating can be estimated by superposing uniform initial stresses developed during plating. The results indicate that the residual tensile stresses due to fabrication in coating will be decreased greatly, or even synthetical compressive internal stresses may arise in the coating.

Dynamic analysis of arrest of buckle propagation on a beam on a nonlinear elastic foundation by FEM

Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.