Finite element simulation of wedge indentation

A finite element simulation of frictionless wedge indentation of a copper strip has been carried out under plane strain conditions. The problem was first modelled using an one-pass contact algorithm. The difficulties associated with using this method to model wedge indentation problems are explained. An alternative procedure which alleviates some of the problems associated with the one-pass contact algorithm is proposed for modelling frictionless wedge indentation. Also, a re-meshing procedure which has to be carried out when the distortion of the elements around the indenter becomes significant, is discussed. A sample problem involving indentation of a 4 mm copper strip by a rigid wedge indenter has been modelled and the results are compared with experimental and theoretical results.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.

Finite element simulation of BAW propagation in inhomogeneous plate due to piezoelectric actuation

A set of finite elements (FEs) is formulated to analyze wave propagation through inhomogeneous material when subjected to mechanical, thermal loading or piezo-electric actuation. Elastic, thermal and electrical properties of the materials axe allowed to vary in length and thickness direction. The elements can act both as sensors and actuators. These elements are used to model wave propagation in functionally graded materials (FGM) and the effect of inhomogeneity in the wave is demonstrated. Further, a surface acoustic wave (SAW) device is modeled and wave propagation due to piezo-electric actuation from interdigital transducers (IDTs) is studied.

Partial delamination modeling in composite beams using a finite element method

A new method of modeling partial delamination in composite beams is proposed and implemented using the finite element method. Homogenized cross-sectional stiffness of the delaminated beam is obtained by the proposed analytical technique, including extension-bending, extension-twist and torsion-bending coupling terms, and hence can be used with an existing finite element method. A two noded C1 type Timoshenko beam element with 4 degrees of freedom per node for dynamic analysis of beams is implemented. The results for different delamination scenarios and beams subjected to different boundary conditions are validated with available experimental results in the literature and/or with the 3D finite element simulation using COMSOL. Results of the first torsional mode frequency for the partially delaminated beam are validated with the COMSOL results. The key point of the proposed model is that partial delamination in beams can be analyzed using a beam model, rather than using 3D or plate models. (c) 2013 Elsevier B.V. All rights reserved.

Vave propagation and bandgaps in a parametrically modulated composite laminate

Wave propagation and its frequency bandgaps in a parametrically modulated composite laminate are reported in this paper. The modulated properties under considerations are due to periodic microstructure, for example honeycomb core sandwich composite, which can be parameterized and homogenized in a suitable scale. Wave equations are derived by assuming a third-order shear deformation theory. Homogenization of the wave equations is carried out in the scale of wavelength. In-plane wave and flexural-shear wave dispersions are obtained for a range of values of a stiffness modulation coefficient (alpha). A clear pattern of stop-bands is observed for alpha >= 4. To validate the band-gap phenomena, we take recourse to time domain response obtained from finite element simulation. As predicted by the proposed analytical technique, a distinct correlation between the chosen frequency band and the simulated wave arrival time and amplitude reduction is found. This promises practical applications of the proposed analytical technique to designing parametrically modulated composite laminate for wave suppression. (C) 2009 Elsevier B.V. All rights reserved.

Effect of Dilatation on the Elasto-Plastic Response of Bulk Metallic Glasses Under Indentation

Elasto-plastic response of bulk metallic glasses (BMGs) follows closely the response of granular materials through pressure dependent (or normal stress) yield locus and shear stress induced material dilatation. On a micro-structural level, material dilatation is responsible for stress softening and formation of localized shear band, however its influence on the macro-scale flow and deformation is largely unknown. In this work, we systematically analyze the effect of material dilatation on the gross indentation response of Zr-based BMG via finite element simulation. The strengthening/softening effect on the load-depth response and corresponding stress-strain profiles are presented in light of differences in elastic-plastic regimes under common indenters. Through comparison with existing experimental results, we draw conclusions regarding selection of suitable dilatation parameters for accurately predicting the gross response of BMGs

Interaction between a notch and cylindrical voids in aluminum single crystals: Experimental observations and numerical simulations

A combined 3D finite element simulation and experimental study of interaction between a notch and cylindrical voids ahead of it in single edge notch (tension) aluminum single crystal specimens is undertaken in this work. Two lattice orientations are considered in which the notch front is parallel to the crystallographic 10 (1) over bar] direction. The flat surface of the notch coincides with the (010) plane in one orientation and with the (1 (1) over bar1) plane in the other. Three equally spaced cylindrical voids are placed directly ahead of the notch tip. The predicted load-displacement curves, slip traces, lattice rotation and void growth from the finite element analysis are found to be in good agreement with the experimental observations for both the orientations. Finite element results show considerable through-thickness variation in both hydrostatic stress and equivalent plastic slip which, however, depends additionally on the lattice orientation. The through-thickness variation in the above quantities affects the void growth rate and causes it to differ from the center-plane to the free surface of the specimen. (c) 2012 Elsevier Ltd. All rights reserved.

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3-8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal-India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs. Crown Copyright (C) 2013 Published by Elsevier Ltd. All rights reserved.

Shape Memory Alloy based Active Chiral Composite Cells

Wing morphing is one of the emerging methodology towards improving aerodynamic efficiency of flight vehicle structures. In this paper a morphing structural element is designed and studied which has its origin in the well known chiral structures. The new aspect of design and functionality explored in this paper is that the chiral cell is actuated using thermal Shape Memory Alloy (SMA) actuator wires to provide directional motion. Such structure utilizes the potential of different actuations concepts based on actuator embedded in the chiral structure skin. This paper describes a new class of chiral cell structure with integrated SMA wire for actuation. Chiral topological constructs are obtained by considering passive and active load path decoupling and sub-optimal shape changes. Single cell of chiral honeycomb with actuators are analyzed using finite element simulation results and experiments. To this end, a multi-cell plan-form is characterized showing interesting possibilities in structural morphing applications. The applicability of the developed chiral cell to flexible wing skin, variable stiffness based design and controlling longitudinal-to-transverse stiffness ratio are discussed.

Dynamic response of a piezoelectric flapping wing

Piezo-composite membranes have advantages over motorized flapping where frequencies are high and certain coupling between bending and twisting is useful to generate lift and forward flight. We draw examples of fruit fly and bumble bee. Wings with Piezo ceramic PZT coating are realized. The passive mechanical response of the wing is characterized experimentally and validated using finite element simulation. Piezoelectric actuation with uniform electrode coating is characterized and optimal frequencies for flapping are identified. The experimental data are used in an empirical model and advanced ratio for a flapping insect like condition for various angular orientations is estimated.

A low-power, low-cost soil-moisture sensor using dual-probe heat-pulse technique

This paper presents the development and testing of an integrated low-power and low-cost dual-probe heat-pulse (DPHP) soil-moisture sensor in view of the electrical power consumed and affordability in developing countries. A DPHP sensor has two probes: a heater and a temperature sensor probe spaced 3 mm apart from the heater probe. Supply voltage of 3.3V is given to the heater-coil having resistance of 33 Omega power consumption of 330 mW, which is among the lowest in this category of sensors. The heater probe is 40 mm long with 2 mm diameter and hence is stiff enough to be inserted into the soil. The parametric finite element simulation study was performed to ensure that the maximum temperature rise is between 1 degrees C and 5 degrees C for wet and dry soils, respectively. The discrepancy between the simulation and experiment is less than 3.2%. The sensor was validated with white clay and tested with red soil samples to detect volumetric water-content ranging from 0% to 30%. The sensor element is integrated with low-power electronics for amplifying the output from thermocouple sensor and TelosB mote for wireless communication. A 3.7V lithium ion battery with capacity of 1150 mAh is used to power the system. The battery is charged by a 6V and 300 mA solar cell array. Readings were taken in 30 min intervals. The life-time of DPHP sensor node is around 3.6 days. The sensor, encased in 30 mm x 20 mm x 10 mm sized box, and integrated with electronics was tested independently in two separate laboratories for validating as well as investigating the dependence of the measurement of soil-moisture on the density of the soil. The difference in the readings while repeating the experiments was found out to be less than 0.01%. Furthermore, the effect of ambient temperature on the measurement of soil-moisture is studied experimentally and computationally. (C) 2015 Elsevier B.V. All rights reserved.

Quantitative In Situ Analysis of Deformation in Sliding Metals: Effect of Initial Strain State

Using in situ, high-speed imaging of a hard wedge sliding against pure aluminum, and image analysis by particle image velocimetry, the deformation field in sliding is mapped at high resolution. This model system is representative of asperity contacts on engineered surfaces and die-workpiece contacts in deformation and machining processes. It is shown that large, uniform plastic strains of 1-5 can be imposed at the Al surface, up to depths of 500 mu m, under suitable sliding conditions. The spatial strain and strain rate distributions are significantly influenced by the initial deformation state of the Al, e.g., extent of work hardening, and sliding incidence angle. Uniform straining occurs only under conditions of steady laminar flow in the metal. Large pre-strains and higher sliding angles promote breakdown in laminar flow due to surface fold formation or flow localization in the form of shear bands, thus imposing limits on uniform straining by sliding. Avoidance of unsteady sliding conditions, and selection of parameters like sliding angle, thus provides a way to control the deformation field. Key characteristics of the sliding deformation such as strain and strain rate, laminar flow, folding and prow formation are well predicted by finite element simulation. The deformation field provides a quantitative basis for interpreting wear particle formation. Implications for engineering functionally graded surfaces, sliding wear and ductile failure in metals are discussed.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Real-time Computer Simulation of Three Dimensional Elastostatics using the Finite Point Method

Real-time simulation of deformable solids is essential for some applications such as biological organ simulations for surgical simulators. In this work, deformable solids are approximated to be linear elastic, and an easy and straight forward numerical technique, the Finite Point Method (FPM), is used to model three dimensional linear elastostatics. Graphics Processing Unit (GPU) is used to accelerate computations. Results show that the Finite Point Method, together with GPU, can compute three dimensional linear elastostatic responses of solids at rates suitable for real-time graphics, for solids represented by reasonable number of points.