Energy Loss in the Impact of Elastic Spheres on a Rigid Half-Space in Presence of Adhesion

Adhesion can cause energy losses in asperities or particles coming into dynamic contact resulting in frictional dissipation, even if the deformation occurring is purely elastic. Such losses are of special significance in impact of nanoparticles and friction between surfaces under low contact pressure to hardness ratio. The objective of this work is to study the effect of adhesion during the normal impact of elastic spheres on a rigid half-space, with an emphasis on understanding the mechanism of energy loss. We use finite element method for modeling the impact phenomenon, with the adhesion due to van der Waals force and the short-range repulsion included as body forces distributed over the volume of the sphere. This approach, in contrast with commonly used surface force approximation, helps to model the interactions in a more precise way. We find that the energy loss in impact of elastic spheres is negligible unless there are adhesion-induced instabilities. Significant energy loss through elastic stress waves occurs due to jump-to-contact and jump-out-of-contact instabilities and can even result in capture of the elastic sphere on the half-space.

Note on modelling directionality in piezoresistivity

When computing the change in electrical resistivity of a piezoresistive cubic material embedded in a deforming structure, the piezoresistive and the stress tensors should be in the same coordinate system. While the stress tensor is usually calculated in a coordinate system aligned with the principal axes of a regular structure, the specified piezoresistive coefficients may not be in that coordinate system. For instance, piezoresistive coefficients are usually given in an orthogonal cartesian coordinate system aligned with the <100> crystallographic directions and designers sometimes deliberately orient a crystallographic direction other than <100> along the principal directions of the structure to increase the gauge factor. In such structures, it is advantageous to calculate the piezoresistivity tensor in the coordinate system along which the stress tensors are known rather than the other way around. This is because the transformation of stress will have to be done at every point in the structure but piezoresistivity tensor needs to be transformed only once. Here, using tensor transformation relations, we show how to calculate the piezoresistive tensor along any arbitrary Cartesian coordinate system from the piezoresistive coefficients for the <100> coordinate system. Some of the software packages that simulate the piezoresistive effect do not have interfaces for calculation of the entire piezoresistive tensor for arbitrary directions. This warrants additional work for the user because not considering the complete piezoresisitive tensor can lead to large errors. This is illustrated with an example where the error is as high as 33%. Additionally, for elastic analysis, we used hybrid finite element formulation that estimates stresses more accurately than displacement-based formulation. Therefore, as shown in an example where the change in resistance can be calculated analytically, the percentage error of our piezoresistive program is an order of magnitude lower relative to displacement-based finite element method.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Crack stability in edge-notched clamped beam specimens: modeling and experiments

Stability of a fracture toughness testing geometry is important to determine the crack trajectory and R-curve behavior of the specimen. Few configurations provide for inherent geometric stability, especially when the specimen being tested is brittle. We propose a new geometrical construction called the single edge notched clamped bend specimen (SENCB), a modified form of three point bending, yielding stable cracking under load control. It is shown to be particularly suitable for small-scale structures which cannot be made free-standing, (e.g., thin films, coatings). The SENCB is elastically clamped at the two ends to its parent material. A notch is inserted at the bottom center and loaded in bending, to fracture. Numerical simulations are carried out through extended finite element method to derive the geometrical factor f(a/W) and for different beam dimensions. Experimental corroborations of the FEM results are carried out on both micro-scale and macro-scale brittle specimens. A plot of vs a/W, is shown to rise initially and fall off, beyond a critical a/W ratio. The difference between conventional SENB and SENCB is highlighted in terms of and FEM simulated stress contours across the beam cross-section. The `s of bulk NiAl and Si determined experimentally are shown to match closely with literature values. Crack stability and R-curve effect is demonstrated in a PtNiAl bond coat sample and compared with predicted crack trajectories from the simulations. The stability of SENCB is shown for a critical range of a/W ratios, proving that it can be used to get controlled crack growth even in brittle samples under load control.

Effect of Process Deviations on Performance of Piezoresistive Pressure Sensors

Diaphragm thickness and the corresponding piezoresistor locations change due to over or under etching in bulk micromachined piezoresistive pressure sensor which intern influences the device performance. In the present work, variation of sensitivity and nonlinearity of a micro electro mechanical system low pressure sensor is investigated. The sensor is modeled using finite element method to analyze the variation of sensitivity and nonlinearity with diaphragm thickness. To verify the simulated results, the sensors with different diaphragm thicknesses are fabricated. The models are verified by comparing the calculated results with experimental data. This study is potentially useful for the researchers as most of the times the diaphragm is either over-etched or under-etched due to inherent variation in wafer thickness and involving manual operations.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

Role of interface curvature on stress distribution under indentation for ZrN/Zr multilayer coating

Contact damage in curved interface nano-layeredmetal/nitride (150 (ZrN)/10 (Zr) nm) multilayer is investigated in order to understand the role of interface morphology on contact damage under indentation. A finite element method (FEM) model was formulated with different wavelengths of 1000 nm, 500 nm, 250 nm and common height of 50 nm, which gives insight on the effect of different curvature on stress field generated under indentation. Elastic-plastic properties were assigned to the metal layer and substrate while the nitride layer was assigned perfectly elastic properties. Curved interface multilayers show delamination along the metal/nitride interface and vertical cracks emanating from the ends of the delamination. FEM revealed the presence of tensile stress normal to the interface even under the contact, along with tensile radial stresses, both present at the valley part of the curve, which leads to vertical cracks associated with interfacial delamination. Stress enhancement was seen to be relatively insensitive to curvature. (C) 2014 Elsevier B.V. All rights reserved.

Matrix Crack Detection in Composite Plate with Spatially Random Material Properties using Fractal Dimension

Fractal dimension based damage detection method is investigated for a composite plate with random material properties. Composite material shows spatially varying random material properties because of complex manufacturing processes. Matrix cracks are considered as damage in the composite plate. Such cracks are often seen as the initial damage mechanism in composites under fatigue loading and also occur due to low velocity impact. Static deflection of the cantilevered composite plate with uniform loading is calculated using the finite element method. Damage detection is carried out based on sliding window fractal dimension operator using the static deflection. Two dimensional homogeneous Gaussian random field is generated using Karhunen-Loeve (KL) expansion to represent the spatial variation of composite material property. The robustness of fractal dimension based damage detection method is demonstrated considering the composite material properties as a two dimensional random field.

Analysis of Entropy Generation During Mixed Convective Heat Transfer of Nanofluids Past a Rotating Circular Cylinder

The entropy generation due to mixed convective heat transfer of nanofluids past a rotating circular cylinder placed in a uniform cross stream is investigated via streamline upwind Petrov-Galerkin based finite element method. Nanosized copper (Cu) particles suspended in water are used with Prandtl number (Pr)=6.9. The computations are carried out at a representative Reynolds number (Re) of 100. The dimensionless cylinder rotation rate, a, is varied between 0 and 2. The range of nanoparticle volume fractions (phi) considered is 0 <= phi <= 5%. Effect of aiding buoyancy is brought about by considering two fixed values of the Richardson number (Ri) as 0.5 and 1.0. A new model for predicting the effective viscosity and thermal conductivity of dilute suspensions of nanoscale colloidal particles is presented. The model addresses the details of the agglomeration-deagglomeration in tune with the pertinent variations in the effective particulate dimensions, volume fractions, as well as the aggregate structure of the particulate system. The total entropy generation is found to decrease sharply with cylinder rotation rates and nanoparticle volume fractions. Increase in nanoparticle agglomeration shows decrease in heat transfer irreversibility. The Bejan number falls sharply with increase in alpha and phi.

Matrix crack detection in spatially random composite structures using fractal dimension

Fractal dimension based damage detection method is studied for a composite structure with random material properties. A composite plate with localized matrix crack is considered. Matrix cracks are often seen as the initial damage mechanism in composites. Fractal dimension based method is applied to the static deformation curve of the structure to detect localized damage. Static deflection of a cantilevered composite plate under uniform loading is calculated using the finite element method. Composite material shows spatially varying random material properties because of complex manufacturing processes. Spatial variation of material property is represented as a two dimensional homogeneous Gaussian random field. Karhunen-Loeve (KL) expansion is used to generate a random field. The robustness of fractal dimension based damage detection methods is studied considering the composite plate with spatial variation in material properties.

Compressive flow behavior of Al-Si based alloy: Role of heat treatment

The flow characteristics of a near eutectic Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of the flow behavior on heat treatment is studied by testing the alloy in non-heat treated (NHT) and heat treated (HT) conditions. The heat treatment has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. It is observed that the strength of the alloy increases with increase in strain rate and it increases more rapidly above the strain rate of 10(-1) s(-1) in HT condition at all the temperatures, and at 100 degrees C and 200 degrees C in NHT condition. The thermally dependent process taking place in the HT matrix is responsible for the observed greater SRS in HT condition. The alloy in HT condition exhibits a larger work hardening rate than in NHT condition during initial stages of straining. However, the hardening rate decreases more sharply at higher strains in HT condition due to precipitate shearing and higher rate of Si particle fracture. Thermal hardening is observed at 200 degrees C in NHT condition due to precipitate formation, which results in increased SRS at higher temperatures. Thermal softening is observed in HT condition at 200 C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by a finite element method support the experimentally observed particle and matrix fracture behavior. A negative SRS and serrated flow are observed in the lower strain rate regime (3 x 10(-4)-10(-2) s(-1)) at RT and 100 degrees C, in both NHT and HT conditions. The observations show that both dynamic strain aging (DSA) and precipitate shearing play a role in serrated flow. (C) 2015 Elsevier B.V. All rights reserved.

Mixed Mode Analysis of Obliquely Oriented Edge Cracked Semicircular Disk Under Compression

A new mixed-mode compression fracture specimen, obliquely oriented edge cracked semicircular disk (OECSD) is analyzed by extending pure opening mode configuration of edge cracked semicircular disk (ECSD) under Hertzian compression. Photoelastic experiments are conducted on two different specimens of OECSD of same size and different crack lengths and inclinations. Finite element method (FEM) is used to solve a number of cases of the problem varying crack length and crack inclination. FE results show a good match with experiments. Inclination of edge crack in OECSD can be so made as to obtain any mode-mixity ratio between zero and one and beyond for any crack length. The new specimen can be used for fracture testing under compression more conveniently than the existing ones in several ways.

A Study of Pressure-Dependent Squeeze Film Stiffness as a Resonance Modulator Using Static and Dynamic Measurements

We report on the resonant frequency modulation of inertial microelectromechanical systems (MEMS) structures due to squeeze film stiffness over a range of working pressures. Squeeze film effects have been studied extensively, but mostly in the context of damping and Q-factor determination of dynamic MEMS structures, typically suspended over a fixed substrate with a very thin air gap. Here, we show with experimental measurements and analytical calculations how the pressure-dependent air springs (squeeze film stiffness) change the resonant frequency of an inertial MEMS structure by as much as five times. For capturing the isolated effect of the squeeze film stiffness, we first determine the static stiffness of our structure with atomic force microscope probing and then study the effect of the air spring by measuring the dynamic response of the structure, thus finding the resonant frequencies while varying the air pressure from 1 to 905 mbar. We also verify our results by analytical and Finite Element Method calculations. Our findings show that the pressure-dependent squeeze film stiffness can affect a rather huge range of frequency modulation (>400%) and, therefore, can be used as a design parameter for exploiting this effect in MEMS devices. 2014-0310]

Role of Si modification on the compressive flow behavior of Al-Si based alloy

The flow characteristics of a near-eutectic heat-treated Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of flow behavior on modification is examined by testing the alloy in both the unmodified and modified conditions. Modification has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. The strength of the alloy is found to increase with increase in strain rate for both the conditions. The increase is more rapid above the strain rate of 10(-1) s(-1) for the unmodified alloy at all the temperatures. This rapid increase is observed at 1 s(-1) at RT and 100 degrees C, and at 10(-2) s(-1) at 200 degrees C for the modified alloy. The thermally dependent process of the Al matrix is rate controlling in the unmodified alloy. On the other hand, the thermally dependent process of both Al matrix and Si particles are rate controlling, which is responsible for the higher strain rate sensitivity (SRS) in the modified alloy. The unmodified alloy exhibits a larger work hardening rate than the modified alloy during the initial stages of straining due to fiber loading of unmodified Si particles. However, the hardening rate decreases sharply at higher strains for the unmodified alloy due to a higher rate of Si particle fracture. Thermal softening is observed for both alloys at 200 degrees C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by microstructure based finite element method support the experimentally observed particle and matrix fracture behavior. Negative SRS and serrated flow are observed at lower strain rate regime (3 x 10(-4) to 10(-2) s(-1)) at RT and 100 degrees C, in both alloys. The critical onset strain is found to be lower and the magnitude of serration is found to be higher for the modified alloy, which suggests that, in addition to dynamic strain aging, Si particle size and morphology also play a role in serrated flow. (C) 2015 Elsevier Inc All rights reserved.