Simultaneous measurements of thickness and temperature profile of the lubricant film at chip-tool interface during machining process using luminescent sensors

Simultaneous measurements of thickness and temperature profile of the lubricant film at chip-tool interface during machining have been studied in this experimental programme. Conventional techniques such as thermography can only provide temperature measurement under controlled environment in a laboratory and without the addition of lubricant. The present study builds on the capabilities of luminescent sensors in addition to direct image based observations of the chip-tool interface. A suite of experiments conducted using different types of sensors are reported in this paper, especially noteworthy are concomitant measures of thickness and temperature of the lubricant. (C) 2014 Elsevier Ltd.

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.

Stability of the flow in a soft tube deformed due to an applied pressure gradient

A linear stability analysis is carried out for the flow through a tube with a soft wall in order to resolve the discrepancy of a factor of 10 for the transition Reynolds number between theoretical predictions in a cylindrical tube and the experiments of Verma and Kumaran J. Fluid Mech. 705, 322 (2012)]. Here the effect of tube deformation (due to the applied pressure difference) on the mean velocity profile and pressure gradient is incorporated in the stability analysis. The tube geometry and dimensions are reconstructed from experimental images, where it is found that there is an expansion and then a contraction of the tube in the streamwise direction. The mean velocity profiles at different downstream locations and the pressure gradient, determined using computational fluid dynamics, are found to be substantially modified by the tube deformation. The velocity profiles are then used in a linear stability analysis, where the growth rates of perturbations are calculated for the flow through a tube with the wall modeled as a neo-Hookean elastic solid. The linear stability analysis is carried out for the mean velocity profiles at different downstream locations using the parallel flow approximation. The analysis indicates that the flow first becomes unstable in the downstream converging section of the tube where the flow profile is more pluglike when compared to the parabolic flow in a cylindrical tube. The flow is stable in the upstream diverging section where the deformation is maximum. The prediction for the transition Reynolds number is in good agreement with experiments, indicating that the downstream tube convergence and the consequent modification in the mean velocity profile and pressure gradient could reduce the transition Reynolds number by an order of magnitude.

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.

COMPUTATIONAL MODELING OF A MULTI-LAYERED PIEZO-COMPOSITE BEAM MADE UP OF MFC

The Variational Asymptotic Method (VAM) is used for modeling a coupled non-linear electromechanical problem finding applications in aircrafts and Micro Aerial Vehicle (MAV) development. VAM coupled with geometrically exact kinematics forms a powerful tool for analyzing a complex nonlinear phenomena as shown previously by many in the literature 3 - 7] for various challenging problems like modeling of an initially twisted helicopter rotor blades, matrix crack propagation in a composite, modeling of hyper elastic plates and various multi-physics problems. The problem consists of design and analysis of a piezocomposite laminate applied with electrical voltage(s) which can induce direct and planar distributed shear stresses and strains in the structure. The deformations are large and conventional beam theories are inappropriate for the analysis. The behavior of an elastic body is completely understood by its energy. This energy must be integrated over the cross-sectional area to obtain the 1-D behavior as is typical in a beam analysis. VAM can be used efficiently to approximate 3-D strain energy as closely as possible. To perform this simplification, VAM makes use of thickness to width, width to length, width multiplied by initial twist and strain as small parameters embedded in the problem definition and provides a way to approach the exact solution asymptotically. In this work, above mentioned electromechanical problem is modeled using VAM which breaks down the 3-D elasticity problem into two parts, namely a 2-D non-linear cross-sectional analysis and a 1-D non-linear analysis, along the reference curve. The recovery relations obtained as a by-product in the cross-sectional analysis earlier are used to obtain 3-D stresses, displacements and velocity contours. The piezo-composite laminate which is chosen for an initial phase of computational modeling is made up of commercially available Macro Fiber Composites (MFCs) stacked together in an arbitrary lay-up and applied with electrical voltages for actuation. The expressions of sectional forces and moments as obtained from cross-sectional analysis in closed-form show the electro-mechanical coupling and relative contribution of electric field in individual layers of the piezo-composite laminate. The spatial and temporal constitutive law as obtained from the cross-sectional analysis are substituted into 1-D fully intrinsic, geometrically exact equilibrium equations of motion and 1-D intrinsic kinematical equations to solve for all 1-D generalized variables as function of time and an along the reference curve co-ordinate, x(1).

A thermodynamic framework for the evolution of damage in concrete under fatigue

A closed-form expression for the dual of dissipation potential is derived within the framework of irreversible thermodynamics using the principles of dimensional analysis and self-similarity. Through this potential, a damage evolution law is proposed for concrete under fatigue loading using the concepts of damage mechanics in conjunction with fracture mechanics. The proposed law is used to compute damage in a volume element when a member is subjected to fatigue loading. The evolution of damage from microcracking to macrocracking of the entire member is captured through a series of volume elements failing one after the other. The number of loading cycles to failure of the member is obtained as the summation of number of cycles to failure for each individual volume element. A parametric study is conducted to determine the effect of the size of the volume element on the model's prediction of fatigue life. A global damage index is also defined, and the residual moment carrying capacity of damaged beams is evaluated. Through a deterministic sensitivity analysis, it is found that the load range and maximum aggregate size are the most influencing parameters on the fatigue life of a plain concrete beam.

VAM Applied to Dimensional Reduction of Nonlinear Multifunctional Film-fabric Laminates

This work aims at asymptotically accurate dimensional reduction of non-linear multi-functional film-fabric laminates having specific application in design of envelopes for High Altitude Airships (HAA). The film-fabric laminate for airship envelope consists of a woven fabric core coated with thin films on each face. These films provide UV protection and Helium leakage prevention, while the core provides required structural strength. This problem is both geometrically and materially non-linear. To incorporate the geometric non-linearity, generalized warping functions are used and finite deformations are allowed. The material non-linearity is handled by using hyper-elastic material models for each layer. The development begins with three-dimensional (3-D) nonlinear elasticity and mathematically splits the analysis into a one-dimensional through-the-thickness analysis and a two-dimensional (2-D) plate analysis. The through-the-thickness analysis provides the 2-D constitutive law which is then given as an input to the 2-D reference surface analysis. The dimensional reduction is carried out using Variational Asymptotic Method (VAM) for moderate strains and very small thickness-to-wavelength ratio. It features the identification and utilization of additional small parameters such as ratio of thicknesses and stiffness coefficients of core and films. Closed form analytical expressions for warping functions and 2-D constitutive law of the film-fabric laminate are obtained.

Global response sensitivity analysis using probability distance measures and generalization of Sobol's analysis

The study introduces two new alternatives for global response sensitivity analysis based on the application of the L-2-norm and Hellinger's metric for measuring distance between two probabilistic models. Both the procedures are shown to be capable of treating dependent non-Gaussian random variable models for the input variables. The sensitivity indices obtained based on the L2-norm involve second order moments of the response, and, when applied for the case of independent and identically distributed sequence of input random variables, it is shown to be related to the classical Sobol's response sensitivity indices. The analysis based on Hellinger's metric addresses variability across entire range or segments of the response probability density function. The measure is shown to be conceptually a more satisfying alternative to the Kullback-Leibler divergence based analysis which has been reported in the existing literature. Other issues addressed in the study cover Monte Carlo simulation based methods for computing the sensitivity indices and sensitivity analysis with respect to grouped variables. Illustrative examples consist of studies on global sensitivity analysis of natural frequencies of a random multi-degree of freedom system, response of a nonlinear frame, and safety margin associated with a nonlinear performance function. (C) 2015 Elsevier Ltd. All rights reserved.

Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

Stability analysis of thick piezoelectric metal based FGM plate using first order and higher order shear deformation theory

This paper presents the stability analysis of functionally graded plate integrated with piezoelectric actuator and sensor at the top and bottom face, subjected to electrical and mechanical loading. The finite element formulation is based on first order and higher order shear deformation theory, degenerated shell element, von-Karman hypothesis and piezoelectric effect. The equation for static analysis is derived by using the minimum energy principle and solutions for critical buckling load is obtained by solving eigenvalue problem. The material properties of the functionally graded plate are assumed to be graded along the thickness direction according to simple power law function. Two types of boundary conditions are used, such as SSSS (simply supported) and CSCS (simply supported along two opposite side perpendicular to the direction of compression and clamped along the other two sides). Sensor voltage is calculated using present analysis for various power law indices and FG (functionally graded) material gradations. The stability analysis of piezoelectric FG plate is carried out to present the effects of power law index, material variations, applied mechanical pressure and piezo effect on buckling and stability characteristics of FG plate.

Handling large variations in mechanics: Some applications

There is a need to use probability distributions with power-law decaying tails to describe the large variations exhibited by some of the physical phenomena. The Weierstrass Random Walk (WRW) shows promise for modeling such phenomena. The theory of anomalous diffusion is now well established. It has found number of applications in Physics, Chemistry and Biology. However, its applications are limited in structural mechanics in general, and structural engineering in particular. The aim of this paper is to present some mathematical preliminaries related to WRW that would help in possible applications. In the limiting case, it represents a diffusion process whose evolution is governed by a fractional partial differential equation. Three applications of superdiffusion processes in mechanics, illustrating their effectiveness in handling large variations, are presented.

Lifshitz transition and modulation of electronic and transport properties of bilayer graphene by sliding and applied normal compressive strain

Using density functional theory (DFT) we investigate the changes in electronic and transport properties of graphene bilayer caused by sliding one of the layers. Change in stacking pattern breaks the lattice symmetry, which results in Lifshitz transition together with the modulation of the electronic structure. Going from AA to AB stacking by sliding along armchair direction leads to a drastic transition in electronic structure from linear to parabolic dispersion. Our transport calculations show a significant change in the overall transmission value for large sliding distances along zigzag direction. The increase in interlayer coupling with normal compressive strain increases the overlapping of conduction and valence band, which leads to further shift in the Dirac points and an enhancement in the Lifshitz transition. The ability to tune the topology of band structure by sliding and/or applying normal compressive strain will open doors for controlled tuning of many physical phenomenon such as Landau levels and quantum Hall effect in graphene. (C) 2015 Elsevier Ltd. All rights reserved.

Effect of applied pressure on densification of monolithic ZrCx ceramic by reactive hot pressing

The effect of applied pressure on reactive hot pressing (RHP) of zirconium (Zr):graphite (C) in molar ratios of 1:0.5, 1:0.67, 1:0.8, and 1:1 was studied at 1200 degrees C for 60 min. The relative density achievable increased with increasing pressure and ranged from 99% at 4 MPa for ZrC0.5 to 93% for stoichiometric ZrC at 100 MPa. The diminishing influence of pressure on the final density with increasing stoichiometry is attributed to two causes: the decreasing initial volume fraction of the plastically deforming Zr metal which leads to the earlier formation of a contiguous, stress shielding carbide skeleton and the larger molar volume shrinkage during reaction which leads to pore formation in the final stages. A numerical model of the creep densification of a dynamically evolving microstructure predicts densities that are consistent with observations and confirm that the availability of a soft metal is primarily responsible for the achievement of such elevated densification during RHP. The ability to densify nonstoichiometric compositions like ZrC0.5 at pressures as low as 4 MPa offers an alternate route to fabricating dense nonstoichiometric carbides.