Elasto-plastic contact stress analysis of joints subjected to cyclic loading

Rigorous elastic-plastic finite element analysis of joints subjected to cyclic loading is carried out. An incremental-iterative algorithm is developed in a modular form combining elasto-plastic material behaviour and contact stress analysis. For the case of the interference fit, the analysis sequentially carries out insertion of the pin and application of the load on the joint, covering possible initiation of separation (and/or yielding) and progressively the receding/advancing contact at the pin-plate interface. Deformations of both the plate and the pin are considered in the analysis. Numerical examples are presented for the case of an interference fit pin in a large plate under remote cyclic tension, and for an interference fit pin lug joint subjected to cyclic loading. A detailed study is carried out for the latter problem considering the effect of change in contact/separation at the pin-plate interface on local stresses, strains and redistribution of these stresses with the spread of a plastic zone. The results of the study are a useful input for the estimation of the fatigue life of joints. Copyright (C) 1996 Elsevier Science Ltd

Interaction analysis of strip footing resting on a non-homogeneous elastic medium

In this paper, a plane stress solution for the interaction analysis of strip footing resting on (i) a non-homogeneous elastic half-plane and (ii) a non-homogeneous elastic layer resting on a rigid stratum has been presented. The analysis has been done using a combined analytical and FEM method in which the discretization of the half-plane is not required and thereby minimizes the computational efforts considerably. The contact pressure distribution and the settlement profile for the selected cases of varying modulus half-plane, which has more relevance to foundation engineering, have been given. Experimental verification through a photoelastic method of stress analysis has been carried out for the case of footing on Gibson elastic half-plane, and the contact pressure distribution thus obtained has been compared with the theoretical results. Copyright (C) 1996 Elsevier Science Ltd

Role of elasto-plastic analysis under cyclic loading in fatigue crack growth studies

Linear Elastic Fracture Mechanics (LEFM) has been widely used in the past for fatigue crack growth studies, but this is acceptable only in situations which are within small scale yielding (SSY). In many practical structural components, conditions of SSY could be violated and one has to look for fracture criteria based on elasto-plastic analysis. Crack closure phenomenon, one of the most striking discoveries based on inelastic deformations during crack growth, has significant effect on fatigue crack growth rate. Numerical simulation of this phenomenon is computationally intensive and involved but has been successfully implemented. Stress intensity factors and strain energy release rates lose their meaning, J-integral (or its incremental) values are applicable only in specific situations, whereas alternate path independent integrals have been proposed in the literature for use with elasto-plastic fracture mechanics (EPFM) based criteria. This paper presents certain salient features of two independent finite element (numerical) studies of relevance to fatigue crack growth, where elasto-plastic analysis becomes significant. These problems can only be handled in the current day computational environment, and would have been only a dream just a few years ago.

Diffusion and Growth of the mu Phase (Ni6Nb7) in the Ni-Nb System

Incremental diffusion couple experiments are conducted to determine the average interdiffusion coefficient and the intrinsic diffusion coefficients of the species in the Ni6Nb7 (mu phase) in the Ni-Nb system. Further, the tracer diffusion coefficients are calculated from the knowledge of thermodynamic parameters. The diffusion rate of Ni is found to be higher than that of Nb, which indicates higher defect concentration in the Ni sublattice.

Aerodynamics of delta wings with flaps at hypersonic speeds

Control surface effectiveness is an important parameter for any aeroplane. For a hypersonic aircraft, though the power required to operate the flaps is determined by low speed flying conditions, it is imperative to know the effect of flaps at hypersonic speeds. Hence, studies have been done on this topic by aerodynamicists for over 40 years. In spite of this, only a limited data is available in the literature on this subject. This paper discusses the experimental study of the effect of sweep on the aerodynamic characteristics of thin slab delta wings with flaps at hypersonic speeds. For the purpose of this investigation, a novel special thin six-component balance, which has a thickness of 4mm and can be housed inside wings with 8mm thickness, has been designed. The wings had a sweep of 76degrees, 70degrees and 65degrees, t/c of 0.053 and flaps with 12% of wing area and 12% of wing chord. Testing were done at Mach 8.2, Re number of 2.13 x 10(6) (based on chord), from alpha = -12degrees to 12degrees and flap angle of 20degrees, 30degrees and 40degrees. Separation lengths, measured from Schlieren pictures, clearly show that there is 'no appreciable' effect of sweep on them. Also, using a simple local flow field calculation, the separation has been identified to be transitional in nature. These features of separation reflect in the force data. Because of the small separation length, the flaps (inspite of their small size) were very effective in generating additional C-N, C-M and C-l, which increased with increase in flap angle. In general, the C-N, C-M and X-CP were unaffected by sweep for symmetric flap deflection at positive incidences and asymmetric flap case, For symmetric flap case at negative incidences, only C-N was not influenced by the sweep but C-M decreased and X-CP moved upstream as the sweep is decreased, The wing with lower sweep produces higher CA and lower (L/D)(max) for both symmetric and asymmetric flaps. The rolling moment and adverse yaw increased with decrease in sweep for asymmetric flap deflection. Newtonian theory is shown to be incapable of predicting the effect of sweep on C-l, C-n and on the incremental values of C-N, C-M and C-A. In conclusion, it can be said that a small flap is generally adequate for hypersonic aeroplanes provided they operate at altitudes where transitional and turbulent separation can be expected to occur. This would make the flaps effective and thus enable ample control authority.

A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Nonlinear Finite-Element Modeling of Batter Piles under Lateral Load

In this paper, a finite-element model is developed in which the nonlinear soil behavior is represented by a hyperbolic relation for static load condition and modified hyperbolic relation, which includes both degradation and gap for a cyclic load condition. Although batter piles are subjected to lateral load, the soil resistance is also governed by axial load, which is incorporated by considering the P-Δ moment and geometric stiffness matrix. By adopting the developed numerical model, static and cyclic load analyses are performed adopting an incremental-iterative procedure where the pile is idealized as beam elements and the soil as elastoplastic spring elements. The proposed numerical model is validated with published laboratory and field pile test results under both static and cyclic load conditions. This paper highlights the importance of the degradation factor and its influence on the soil resistance-displacement (p-y) curve, number of cycles of loading, and cyclic load response.

A nonlinear spectral finite element model for analysis of wave propagation in solid with internal friction and dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.

Dynamic subspace clustering for very large high-dimensional databases

Emerging high-dimensional data mining applications needs to find interesting clusters embeded in arbitrarily aligned subspaces of lower dimensionality. It is difficult to cluster high-dimensional data objects, when they are sparse and skewed. Updations are quite common in dynamic databases and they are usually processed in batch mode. In very large dynamic databases, it is necessary to perform incremental cluster analysis only to the updations. We present a incremental clustering algorithm for subspace clustering in very high dimensions, which handles both insertion and deletions of datapoints to the backend databases.

Synthesis and characterization of novel cationic lipid and cholesterol-coated gold nanoparticles and their interactions with dipalmitoylphosphatidylcholine membranes

Novel gold nanoparticles bearing cationic single-chain, double-chain, and cholesterol based amphiphilic units have been synthesized. These nanoparticles represent size-stable entities in which various cationic lipids have been immobilized through their thiol group onto the gold nanoparticle core. The resulting colloids have been characterized by UV-vis, (1)H NMR, FT-IR spectroscopy, and transmission electron microscopy. The average size of the resultant nanoparticles could be controlled by the relative bulkiness of the capping agent. Thus, the average diameters of the nanoparticles formed from the cationic single-chain, double-chain, and cholesterol based thiolate-coated materials were 5.9,2.9, and 2.04 nm, respectively. We also examined the interaction of these cationic gold nanoparticles with vesicular membranes generated from dipalmitoylphosphatidylcholine (DPPC) lipid suspensions. Nanoparticle doped DPPC vesicular suspensions displayed a characteristic surface plasmon band in their UV-vis spectra. Inclusion of nanoparticles in vesicular suspensions led to increases in the aggregate diameters, as evidenced from dynamic light scattering. Differential scanning calorimetric examination indicated that incorporation of single-chain, double-chain, and cholesteryl-linked cationic nanoparticles exert variable effects on the DPPC melting transitions. While increased doping of single-chain nanoparticles in DPPC resulted in the phases that melt at higher temperatures, inclusion of an incremental amount of double-chain nanoparticles caused the lowering of the melting temperature of DPPC. On the other hand, the cationic cholesteryl nanoparticle interacted with DPPC in membranes in a manner somewhat analogous to that of cholesterol itself and caused broadening of the DPPC melting transition.

Severe wear of a near eutectic aluminium-silicon alloy

The severe wear of a near eutectic aluminium silicon alloy is explored using a range of electron microscopic, spectroscopic and diffraction techniques to identify the residually strained and unstrained regions, microcracks and oxidized regions in the subsurface. In severe wear the contact pressure exceeds the elastic shakedown limit. Under this condition the primary and eutectic silicon particles fragment drastically. The fragments are transported by the matrix as it undergoes incremental straining with each cyclic contact at the asperity level. The grains are refined from similar to 2000 nm in the bulk to 30 nm in the near surface region. A large reduction in the interparticle distance compared with that for a milder stage of wear gives rise to high strain gradients which contribute to an enhancement of the dislocation density. The resulting regions of very high strain in the boundaries of the recrystallized grains as well as within the subgrains lead to the formation of microvoidskracks. This is accompanied by the formation of brittle oxides at these subsurface interfaces due to enhanced diffusion of oxygen. We believe that the abundance of such microcracks in the near surface region, primed by severe plastic deformation, is what distinguishes a severe wear regime from mild wear. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Micro stepping of Shape Memory Alloy based Poly Phase motor

In literature we find broadly two types of shape memory alloy based motors namely limited rotation motor and unlimited rotation motor. The unlimited rotation type SMA based motor reported in literature uses SMA springs for actuation. An attempt has been made in this paper to develop an unlimited rotation type balanced poly phase motor based on SMA wire in series with a spring in each phase. By isolating SMA actuation and spring action we are able achieve a constant force by the SMA wire through out its range of operation. The Poly phase motor can be used in stepping mode for generating incremental motion and servo mode for generating continuous motion. A method of achieving servo motion by micro stepping is presented. Micro stepping consists of controlling single-phase temperature with a position feedback. The motor has been modeled with a new approach to the SMA wire Hysterysis model. Motor is simulated for different responses and the results are compared with the experimental data.

Modeling the Field Emission Current Fluctuation in Carbon Nanotube Thin Films

Owing to their distinct properties, carbon nanotubes (CNTs) have emerged as promising candidate for field emission devices. It has been found experimentally that the results related to the field emission performance show variability. The design of an efficient field emitting device requires the analysis of the variabilities with a systematic and multiphysics based modeling approach. In this paper, we develop a model of randomly oriented CNTs in a thin film by coupling the field emission phenomena, the electron-phonon transport and the mechanics of single isolated CNT. A computational scheme is developed by which the states of CNTs are updated in time incremental manner. The device current is calculated by using Fowler-Nordheim equation for field emission to study the performance at the device scale.

Influence of Spatial Scales on the Performance of Saturated Subsurface Flow and Transport Models

Numerical modeling of saturated subsurface flow and transport has been widely used in the past using different numerical schemes such as finite difference and finite element methods. Such modeling often involves discretization of the problem in spatial and temporal scales. The choice of the spatial and temporal scales for a modeling scenario is often not straightforward. For example, a basin-scale saturated flow and transport analysis demands larger spatial and temporal scales than a meso-scale study, which in turn has larger scales compared to a pore-scale study. The choice of spatial-scale is often dictated by the computational capabilities of the modeler as well as the availability of fine-scale data. In this study, we analyze the impact of different spatial scales and scaling procedures on saturated subsurface flow and transport simulations.

Sensitivity of terrestrial water and energy budgets to CO(2)-physiological forcing: an investigation using an offline land model

Increasing concentrations of atmospheric carbon dioxide (CO(2)) influence climate by suppressing canopy transpiration in addition to its well- known greenhouse gas effect. The decrease in plant transpiration is due to changes in plant physiology (reduced opening of plant stomata). Here, we quantify such changes in water flux for various levels of CO(2) concentrations using the National Center for Atmospheric Research's (NCAR) Community Land Model. We find that photosynthesis saturates after 800 ppmv (parts per million, by volume) in this model. However, unlike photosynthesis, canopy transpiration continues to decline at about 5.1% per 100 ppmv increase in CO(2) levels. We also find that the associated reduction in latent heat flux is primarily compensated by increased sensible heat flux. The continued decline in canopy transpiration and subsequent increase in sensible heat flux at elevated CO(2) levels implies that incremental warming associated with the physiological effect of CO(2) will not abate at higher CO(2) concentrations, indicating important consequences for the global water and carbon cycles from anthropogenic CO(2) emissions.