A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

A meshless adaptive multiscale method for fracture

The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111) plane of an FCC crystal in the fine scale region. The Lennard-Jones potential is used to model the atom-atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. (C) 2014 Elsevier B. V. All rights reserved.

Beyond classical dynamic structural plasticity using mesh-free modelling techniques

The problem of modelling the transient response of an elastic-perfectly-plastic cantilever beam, carrying an impulsively loaded tip mass, is,often referred to as the Parkes cantilever problem 25]; The permanent deformation of a cantilever struck transversely at its tip, Proc. R. Soc. A., 288, pp. 462). This paradigm for classical modelling of projectile impact on structures is re-visited and updated using the mesh-free method, smoothed particle hydrodynamics (SPH). The purpose of this study is to investigate further the behaviour of cantilever beams subjected to projectile impact at its tip, by considering especially physically real effects such as plastic shearing close to the projectile, shear deformation, and the variation of the shear strain along the length and across the thickness of the beam. Finally, going beyond macroscopic structural plasticity, a strategy to incorporate physical discontinuity (due to crack formation) in SPH discretization is discussed and explored in the context of tip-severance of the cantilever beam. Consequently, the proposed scheme illustrates the potency for a more refined treatment of penetration mechanics, paramount in the exploration of structural response under ballistic loading. The objective is to contribute to formulating a computational modelling framework within which transient dynamic plasticity and even penetration/failure phenomena for a range of materials, structures and impact conditions can be explored ab initio, this being essential for arriving at suitable tools for the design of armour systems. (C) 2014 Elsevier Ltd. All rights reserved.

Novel supersonic nozzles for mixing enhancement in supersonic ejectors

Two novel supersonic nozzles Tip Ring Supersonic Nozzle and Elliptic Sharp Tipped Shallow (ESTS) Lobed Nozzle have been developed to enhance mixing at high speeds which is beneficial to supersonic ejectors. A circular ring protruding at the exit of a conical nozzle forms the tip ring nozzle. The innovative ESTS lobed nozzle comprising of four elliptic lobes with sharp tips that do not protrude deep into the core supersonic flow is produced by a novel yet simple methodology. A comparative experimental study is conducted between a conical nozzle, an ESTS lobed nozzle and a tip ring nozzle with exit Mach number of 2.3. For the first time, the three dimensional flow structure from ESTS lobed nozzle and tip ring nozzle is revealed from laser scattering flow visualization experiments on the free jet. A doubling of jet spreading rate is observed in the ESTS lobed nozzle. When applied to a supersonic ejector, both nozzles achieve a 30% increase in entrainment of secondary flow. The loss of compression ratio is 15% for the ESTS lobed nozzle while it is 50% for the tip ring nozzle. Further, the behavior of wall static pressure profile corroborates mixing enhancement. (C) 2014 Elsevier Ltd. All rights reserved.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

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.

Fracture in metallic glasses: mechanics and mechanisms

Significant progress in understanding the mechanical behavior of metallic glasses (MGs) was made over the past decade, particularly on mechanisms of plastic deformation. However, recent research thrust has been on exploring the mechanics and physics of fracture. MGs can be very brittle with K-Ic values similar to silicate glasses and ceramics or very tough with K-Ic akin to high toughness crystalline metals. Even the tough MGs can become brittle with structural relaxation following annealing at temperatures close to glass transition temperature (T-g). Detailed experimental studies coupled with complementary numerical simulations of the recent past have provided insights on the micromechanisms of failure as well as nature of crack tip fields, and established the governing fracture criteria for ductile and brittle glasses. In this paper, the above advances are reviewed and outstanding issues in the context of fracture of amorphous alloys that need to be resolved are identified.

Development of cup shaped microneedle array for transdermal drug delivery

Microneedle technology is one of the attractive methods in transdermal drug delivery. However, the clinical applications of this method are limited owing to: complexity in the preparation of multiple coating solutions, drug leakage while inserting the microneedles into the skin and the outer walls of the solid microneedle can hold limited quantity of drug. Here, the authors present the fabrication of an array of rectangular cup shaped silicon microneedles, which provide for reduced drug leakage resulting in improvement of efficiency of drug delivery and possibility of introducing multiple drugs. The fabricated solid microneedles with rectangular cup shaped tip have a total height of 200 mu m. These cup shaped tips have dimensions: 60 x 60 mu m (length x breadth) with a depth of 60 mu m. The cups are filled with drug using a novel in-house built drop coating system. Successful drug dissolution was observed when the coated microneedle was used on mice. Also, using the above method, it is possible to fill the cups selectively with different drugs, which enables simultaneous multiple drug delivery. (C) 2015 American Vacuum Society.

Multiobjective Optimization of Piezoelectric Bimorph Actuator with Rigid Extension

Production of high tip deflection in a piezoelectric bimorph laminar actuator by applying high voltage is limited by many physical constraints. Therefore, piezoelectric bimorph actuator with a rigid extension of non-piezoelectric material at its tip is used to increase the tip deflection of such an actuator. Research on this type of piezoelectric bending actuator is either limited to first order constitutive relations, which do not include non-linear behavior of piezoelectric element at high electric field, or limited to curve fitting techniques. Therefore, this paper considers high electric field, and analytically models tapered piezoelectric bimorph actuator with a rigid extension of non-piezoelectric material at its tip. The stiffness, capacitance, effective tip deflection, block force, output strain energy, output energy density, input electrical energy and energy efficiency of the actuator are calculated analytically. The paper also discusses the multi-objective optimization of this type of actuator subjected to the mechanical and electrical constraints.

Design and Evaluation of Torsional Probes for Multifrequency Atomic Force Microscopy

Multifrequency atomic force microscopy is a powerful nanoscale imaging and characterization technique that involves excitation of the atomic force microscope (AFM) probe and measurement of its response at multiple frequencies. This paper reports the design, fabrication, and evaluation of AFM probes with a specified set of torsional eigen-frequencies that facilitate enhancement of sensitivity in multifrequency AFM. A general approach is proposed to design the probes, which includes the design of their generic geometry, adoption of a simple lumped-parameter model, guidelines for determination of the initial dimensions, and an iterative scheme to obtain a probe with the specified eigen-frequencies. The proposed approach is employed to design a harmonic probe wherein the second and the third eigen-frequencies are the corresponding harmonics of the first eigen-frequency. The probe is subsequently fabricated and evaluated. The experimentally evaluated eigen-frequencies and associated mode shapes are shown to closely match the theoretical results. Finally, a simulation study is performed to demonstrate significant improvements in sensitivity to the second-and the third-harmonic spectral components of the tip-sample interaction force with the harmonic probe compared to that of a conventional probe.

Direct Measurement of Three-Dimensional Forces in Atomic Force Microscopy

Direct measurement of three-dimensional (3-D) forces between an atomic force microscope (AFM) probe and the sample benefits diverse applications of AFM, including force spectroscopy, nanometrology, and manipulation. This paper presents the design and evaluation of a measurement system, wherein the deflection of the AFM probe is obtained at two points to enable direct measurement of all the three components of 3-D tip-sample forces in real time. The optimal locations for measurement of deflection on the probe are derived for a conventional AFM probe. Further, a new optimal geometry is proposed for the probe that enables measurement of 3-D forces with identical sensitivity and nearly identical resolution along all three axes. Subsequently, the designed measurement system and the optimized AFM probe are both fabricated and evaluated. The evaluation demonstrates accurate measurement of tip-sample forces with minimal cross-sensitivities. Finally, the real-time measurement system is employed as part of a feedback control system to regulate the normal component of the interaction force, and to perform force-controlled scribing of a groove on the surface of polymethyl methacrylate.

NANOLOCALIZED SINGLE CELL MEMBRANE NANOELECTROPORATION

Today single cell research is a great interest to analyze cell to cell or cell to environment behavior with their intracellular compounds, where bulk measurement can provide average value. To deliver biomolecules precise and localized way into single living cell with high transfection rate and high cell viability is a challenging and promisible task for biological and therapeutic research. In this report, we present a nano-localized single cell nano-electroporation technique, where electroporation take place in a very precise and localized area on a single cell membrane to achieve high efficient delivery with high cell viability. We fabricated 60nm gap with 40 nm triangular Indium Tin Oxide (ITO) based nano-eletcrode tip, which can intense electric field in a nano-localized area of a single cell to permeabilize cell membrane and deliver exogenous biomolecules from outside to inside of the cell. This device successfully deliver dyes, proteins into single cell with high cell viability (98%). The process not only control precise delivery mechanism into single cell with membrane reversibility, but also it can provide special, temporal and qualitative dosage control, which might be beneficial for therapeutic and biological cell studies.

EXPERIMENTAL INVESTIGATION OF UNSTEADY FLOW IN A TRANSONIC UNI-STAGE AXIAL COMPRESSOR

Detailed steady and unsteady experimental measurements and analysis were performed on a Single stage Transonic Axial Compressor with asymmetric rotor tip clearance for studying the compressor stall phenomena. The installed compressor had asymmetric tip clearance around the rotor casing varying from about 0.65mm to 1.25mm. A calibrated 5-hole aerodynamic probe was traversed radially at exit of rotor and showed the characteristics of increased flow angle at lower mass flow rates for all the speeds. Mach number distribution and boundary layer effects were also clearly captured. Unsteady measurements for velocity were carried out to study the stall cell behavior using a single component calibrated hotwire probe oriented in axial and tangential directions for choke/free flow and near stall conditions. The hotwire probe was traversed radially across the annulus at inlet to the compressor and showed that the velocity fluctuations were dissimilar when probe was aligned axial and tangential to the flow. Averaged velocities across the annulus showed the reduction in velocity as stall was approached. Axial mean flow velocity decreased across the annulus for all the speeds investigated. Tangential velocity at free flow condition was higher at the tip region due to larger radius. At stall condition, the tangential velocity showed decreased velocities at the tip and slightly increased velocities at the hub section indicating that the flow has breakdown at the tip region of the blade and fluid is accelerated below the blockage zone. The averaged turbulent intensity in axial and tangential flow directions increased from free flow to stall condition for all compressor rated speeds. Fast Fourier Transform (FFT) of the raw signals at stall flow condition showed stall cell and its corresponding frequency of occurrence. The stalling frequency of about half of rotational speed of the rotor along with large tip clearance suggests that modal type stall inception was occurring.

Lower-Bound Axisymmetric Formulation for Geomechanics Problems Using Nonlinear Optimization

A lower-bound limit analysis formulation, by using two-dimensional finite elements, the three-dimensional Mohr-Coulomb yield criterion, and nonlinear optimization, has been given to deal with an axisymmetric geomechanics stability problem. The optimization was performed using an interior point method based on the logarithmic barrier function. The yield surface was smoothened (1) by removing the tip singularity at the apex of the pyramid in the meridian plane and (2) by eliminating the stress discontinuities at the corners of the yield hexagon in the pi-plane. The circumferential stress (sigma(theta)) need not be assumed. With the proposed methodology, for a circular footing, the bearing-capacity factors N-c, N-q, and N-gamma for different values of phi have been computed. For phi = 0, the variation of N-c with changes in the factor m, which accounts for a linear increase of cohesion with depth, has been evaluated. Failure patterns for a few cases have also been drawn. The results from the formulation provide a good match with the solutions available from the literature. (C) 2014 American Society of Civil Engineers.

An automatic algorithm for mixed mode crack growth rate based on drop potential method

A new automatic algorithm for the assessment of mixed mode crack growth rate characteristics is presented based on the concept of an equivalent crack. The residual ligament size approach is introduced to implementation this algorithm for identifying the crack tip position on a curved path with respect to the drop potential signal. The automatic algorithm accounting for the curvilinear crack trajectory and employing an electrical potential difference was calibrated with respect to the optical measurements for the growing crack under cyclic mixed mode loading conditions. The effectiveness of the proposed algorithm is confirmed by fatigue tests performed on ST3 steel compact tension-shear specimens in the full range of mode mixities from pure mode Ito pure mode II. (C) 2015 Elsevier Ltd. All rights reserved.