Finite Element Analysis of Stress Evolution in Al-Si Alloy

A 2D multi-particle model is carried out to understand the effect of microstructural variations and loading conditions on the stress evolution in Al-Si alloy under compression. A total of six parameters are varied to create 26 idealized microstructures: particle size, shape, orientation, matrix temper, strain rate, and temperature. The effect of these parameters is investigated to understand the fracture of Si particles and the yielding of Al matrix. The Si particles are modeled as a linear elastic solid and the Al matrix is modeled as an elasto-plastic solid. The results of the study demonstrate that the increase in particle size decreases the yield strength of the alloy. The particles with high aspect ratio and oriented at 0A degrees and 90A degrees to the loading axis show higher stress values. This implies that the particle shape and orientation are dominant factors in controlling particle fracture. The heat treatment of the alloy is found to increase the stress levels of both particles and matrix. Stress calculations also show that higher particle fracture and matrix yielding is expected at higher strain rate deformation. Particle fracture decreases with increase in temperature and the Al matrix plays an important role in controlling the properties of the alloy at higher temperatures. Further, this strain rate and temperature dependence is more pronounced in the heat-treated microstructure. These predictions are consistent with the experimentally observed Si particle fracture in real microstructure.

Stress Engineering using Si3N4 for stiction free release of SOI beams

We report on the effect of thin silicon nitride (Si3N4) induced tensile stress on the structural release of 200nm thick SOI beam, in the surface micro-machining process. A thin (20nm / 100nm) LPCVD grown Si3N4 is shown to significantly enhance the yield of released beam in wet release technique. This is especially prominent with increase in beam length, where the beams have higher tendency for stiction. We attribute this yield enhancement to the nitride induced tensile stress, as verified by buckling tendency and resonance frequency data obtained from optical profilometry and laser doppler vibrometry.

Investigation on two magnon scattering processes in pulsed laser deposited epitaxial nickel zinc ferrite thin film

Ferromagnetic resonance (FMR) measurements are employed to evaluate the presence of the two magnon scattering contribution in the magnetic relaxation processes of the epitaxial nickel zinc ferrite thin films deposited using pulsed laser deposition (PLD) on the (0 0 1) MgAl2O4 substrate. Furthermore, the reciprocal space mapping reveals the presence of microstructural defects which acts as an origin for the two magnon scattering process in this thin film. The relevance of this scattering process is further discussed for understanding the higher FMR linewidth in the in-plane configuration compared to the out-of-plane configuration. FMR measurements also reveal the presence of competing uniaxial and cubic anisotropy in the studied films.

A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

Non-uniform beams and stiff strings isospectral to axially loaded uniform beams and piano strings

In this paper, we derive analytical expressions for mass and stiffness functions of transversely vibrating clamped-clamped non-uniform beams under no axial loads, which are isospectral to a given uniform axially loaded beam. Examples of such axially loaded beams are beam columns (compressive axial load) and piano strings (tensile axial load). The Barcilon-Gottlieb transformation is invoked to transform the non-uniform beam equation into the axially loaded uniform beam equation. The coupled ODEs involved in this transformation are solved for two specific cases (pq (z) = k (0) and q = q (0)), and analytical solutions for mass and stiffness are obtained. Examples of beams having a rectangular cross section are shown as a practical application of the analysis. Some non-uniform beams are found whose frequencies are known exactly since uniform axially loaded beams with clamped ends have closed-form solutions. In addition, we show that the tension required in a stiff piano string with hinged ends can be adjusted by changing the mass and stiffness functions of a stiff string, retaining its natural frequencies.

A faster algorithm for testing polynomial representability of functions over finite integer rings

Given a function from Z(n) to itself one can determine its polynomial representability by using Kempner function. In this paper we present an alternative characterization of polynomial functions over Z(n) by constructing a generating set for the Z(n)-module of polynomial functions. This characterization results in an algorithm that is faster on average in deciding polynomial representability. We also extend the characterization to functions in several variables. (C) 2015 Elsevier B.V. All rights reserved.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Energy Sharing for Multiple Sensor Nodes With Finite Buffers

We consider the problem of finding optimal energy sharing policies that maximize the network performance of a system comprising of multiple sensor nodes and a single energy harvesting (EH) source. Sensor nodes periodically sense the random field and generate data, which is stored in the corresponding data queues. The EH source harnesses energy from ambient energy sources and the generated energy is stored in an energy buffer. Sensor nodes receive energy for data transmission from the EH source. The EH source has to efficiently share the stored energy among the nodes to minimize the long-run average delay in data transmission. We formulate the problem of energy sharing between the nodes in the framework of average cost infinite-horizon Markov decision processes (MDPs). We develop efficient energy sharing algorithms, namely Q-learning algorithm with exploration mechanisms based on the epsilon-greedy method as well as upper confidence bound (UCB). We extend these algorithms by incorporating state and action space aggregation to tackle state-action space explosion in the MDP. We also develop a cross entropy based method that incorporates policy parameterization to find near optimal energy sharing policies. Through simulations, we show that our algorithms yield energy sharing policies that outperform the heuristic greedy method.

Analysis of Transverse Shear Strains in Pre-Twisted Thick Beams using Variational Asymptotic Method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam reference curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

A finite variant of the Toom Model

We present results for a finite variant of the one-dimensional Toom model with closed boundaries. We show that the steady state distribution is not of product form, but is nonetheless simple. In particular, we give explicit formulas for the densities and some nearest neighbour correlation functions. We also give exact results for eigenvalues and multiplicities of the transition matrix using the theory of R-trivial monoids in joint work with A. Schilling, B. Steinberg and N. M. Thiery.

Differential viability response of prokaryotes and eukaryotes to high strength pulsed magnetic stimuli

The present study examines the efficacy of a high strength pulsed magnetic field (PMF) towards bacterial inactivation in vitro, without compromising eukaryotic cell viability. The differential response of prokaryotes Staphylococcus aureus (MESA), Staphylococcus epidermidis, and Escherichia coli], and eukaryotes C2C12 mouse myoblasts and human mesenchymal stem cells, hMSCs] upon exposure to varying PMF stimuli (1-4 T, 30 pulses, 40 ms pulse duration) is investigated. Among the prokaryotes, similar to 60% and similar to 70% reduction was recorded in the survival of staphylococcal species and E. coli, respectively at 4 T PMF as evaluated by colony forming unit (CPU) analysis and flow cytometry. A 2-5 fold increase in intracellular ROS (reactive oxygen species) levels suggests oxidative stress as the key mediator in PMF induced bacterial death/injury. The 4 T PMF treated staphylococci also exhibited longer doubling times. Both TEM and fluorescence microscopy revealed compromised membranes of PMF exposed bacteria. Under similar PMF exposure conditions, no immediate cytotoxicity was recorded in C2C12 mouse myoblasts and hMSCs, which can be attributed to the robust resistance towards oxidative stress. The ion interference of iron containing bacterial proteins is invoked to analytically explain the PMF induced ROS accumulation in prokaryotes. Overall, this study establishes the potential of PMF as a bactericidal method without affecting eukaryotic viability. This non-invasive stimulation protocol coupled with antimicrobial agents can be integrated as a potential methodology for the localized treatment of prosthetic infections. (C) 2015 Elsevier B.V. All rights reserved.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

NON-DIMENSIONAL KINETOELASTOSTATIC MAPS FOR COMPLIANT MECHANISMS

How do we assess the capability of a compliant mechanism of given topology and shape? The kinetoelastostatic maps proposed in this paper help answer this question. These maps are drawn in 2D using two non-dimensional quantities, one capturing the nonlinear static response and the other the geometry, material, and applied forces. Geometrically nonlinear finite element analysis is used to create the maps for compliant mechanisms consisting of slender beams. In addition to the topology and shape, the overall proportions and the proportions of the cross-sections of the beam segments are kept fixed for a map. The finite region of the map is parameterized using a non-dimensional quantity defined as the slenderness ratio. The shape and size of the map and the parameterized curves inside it indicate the complete kinetoelastostatic capability of the corresponding compliant mechanism of given topology, shape, and fixed proportions. Static responses considered in this paper include input/output displacement, geometric amplification, mechanical advantage, maximum stress, etc. The maps can be used to compare mechanisms, to choose a suitable mechanism for an application, or re-design as may be needed. The usefulness of the non-dimensional maps is presented with multiple applications of different variety. Non-dimensional portrayal of snap-through mechanisms is one such example. The effect of the shape of the cross-section of the beam segments and the role of different segments in the mechanism as well as extension to 3D compliant mechanisms, the cases of multiple inputs and outputs, and moment loads are also explained. The effects of disproportionate changes on the maps are also analyzed.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.