Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.

Transitory second-order reciprocal connection for two surfaces in point contact

This paper investigates the instantaneous spatial higher pair to lower pair substitute-connection which is kinematically equivalent up to acceleration analysis for two smooth surfaces in point contact. The existing first-order equivalent substitute-connection consisting of a Hooke's joint (U-joint) and a spherical joint (S-joint) connected by an additional link is extended up to second-order. A two step procedure is chalked out for achieving this equivalence. First, the existing method is employed for velocity equivalence. In the second step, the two centers of substitution are obtained as a conjugate relationship involving the principal normal curvatures of the surfaces at the contact point and the screw coordinates of the instantaneous screw axis (ISA) of the first-order relative motion. Unlike the classical planar replacement, this particular substitution cannot be done by merely examining the profiles of the contacting surfaces. An illustrative example of a three-link direct-contact mechanism is presented. (C) 2014 Elsevier Ltd. All rights reserved.

Sphere to ring morphological transformation in drying nanofluid droplets in a contact-free environment

Understanding the transients of buckling in drying colloidal suspensions is pivotal for producing new functional microstructures with tunable morphologies. Here, we report first observations and elucidate the buckling instability induced morphological transition (sphere to ring structure) in an acoustically levitated, heated nanosuspension droplet using dynamic energy balance. Droplet deformation featuring the formation of symmetric cavities is initiated by capillary pressure that is two to three orders of magnitude greater than the acoustic radiation pressure, thus indicating that the standing pressure field has no influence on the buckling front kinetics. With an increase in heat flux, the growth rate of surface cavities and their post-buckled volume increase while the buckling time period reduces, thereby altering the buckling pathway and resulting in distinct precipitate structures. However, irrespective of the heating rate, the volumetric droplet deformation exhibits a linear time dependence and the droplet vaporization is observed to deviate from the classical D-2-law.

Characterization of friction at three contact pairs by photoelastic Isotropic point (IP)

Friction coefficient between a circular-disk periphery and V-block surface was determined by introducing the concept of isotropic point (IP) in isochromatic field of the disk under three-point symmetric loading. IP position on the symmetry axis depends on active coefficient of friction during experiment. We extend this work to asymmetric loading of circular disk in which case two frictional contact pairs out of three loading contacts, independently control the unconstrained IP location. Photoelastic experiment is conducted on particular case of asymmetric three-point loading of circular disk. Basics of digital image processing are used to extract few essential parameters from experimental image, particularly IP location. Analytical solution by Flamant for half plane with a concentrated load, is utilized to derive stress components for required loading configurations of the disk. IP is observed, in analytical simulations of three-point asymmetric normal loading, to move from vertical axis to the boundary along an ellipse-like curve. When friction is included in the analysis, IP approaches the center with increase in loading friction and it goes away with increase in support friction. With all these insights, using experimental IP information, friction angles at three contact pairs of circular disk under asymmetric loading, are determined.

Theoretical Insights to Niobium-Doped Monolayer MoS2-Gold Contact

We report a first principles study of the electronic properties for a contact formed between Nb-doped monolayer MoS2 and gold for different doping concentrations. We first focus on the shift of energy levels in band structure and the density of states with respect to the Fermi level for a geometrically optimized 5 x 5 MoS2 supercell for both pristine and Nb-doped structures. The doping is achieved by substituting Mo atoms with Nb atoms at random positions. It is observed that for an experimentally reported sheet hole doping concentration of (rho(2D)) 1.8 x 10(14) cm(-2), the pristine MoS2 converts to degenerate p-type semiconductor. Next, we interface this supercell with six layers of < 111 > cleaved surface of gold to investigate the contact nature of MoS2-Au system. By careful examination of projected band structure, projected density of states, effective potential and charge density difference, we demonstrate that the Schottky barrier nature observed for pure MoS2-Au contact can be converted from n-type to p-type by efficient Nb doping.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Physical Insights on the Ambiguous Metal-Graphene Interface and Proposal for Improved Contact Resistance

The ambiguous behavior of metal-graphene interface has been addressed in this paper using density functional theory and nonequilibrium Green's function formalism. For the first time, the fundamental chemistry of metal-graphene interface, in particular role of sp-hybridized and sp(2)-hybridized carbon atoms, has been emphasized and discussed in detail in this paper. It was discovered that the sp-hybridized sites at the edge of a graphene monolayer contribute to 40% of current conduction when compared with sp(2)-hybridized atom sites in the graphene-metal overlap region. Moreover, we highlighted the insignificance of an additional metal layer, i.e., sandwiched contact, due to lacking sp-hybridized carbon sites. A fundamental way of defining the contact resistance, while keeping chemical bonding in mind, has been proposed. The bonding insight has been further used to propose the novel ways of interfacing metal with graphene, which results in a 40% reduction in contact resistance.

Air microjet system for non-contact force application and the actuation of micro-structures

We demonstrate a non-contact technique to apply calibrated and localized forces in the micro-Newton to milli-Newton range using an air microjet. An electromagnetically actuated diaphragm controlled by a signal generator is used to generate the air microjet. With a nozzle diameter of 150 mu m, the microjet diameter was maintained to a maximum of 1 mm at a distance of 5 mm from the nozzle. The force generated by the microjet was measured using a commercial force sensor to determine the velocity profile of the jet. Axial flow velocities of up to 25 m s(-1) were obtained at distances as long as 6 mm. The microjet exerted a force up to 1 mu N on a poly dimethyl siloxane (PDMS) micropillar (50 mu m in diameter, 157 mu m in height) and 415 mu N on a PDMS membrane (3 mm in diameter, 28 mu m thick). We also demonstrate that from a distance of 6 mm our microjet can exert a peak pressure of 187 Pa with a total force of about 84 mu N on a flat surface with 8 V operating voltage. Out of the cleanroom fabrication and robust design make this system cost effective and durable.

Intrinsic Limit for Contact Resistance in Exfoliated Multilayered MoS2 FET

A new method for the separation of contact resistance (R-contact) into Schottky barrier resistance (R-SB) and interlayer resistance (R-IL) is proposed for multilayered MoS2 FETs. While R-SB varies exponentially with Schottky barrier height (Phi(bn)), R-IL essentially remains unchanged. An empirical model utilizing this dependence of R-contact versus Phi(bn) is proposed and fits to the experimental data. The results, on comparison with the existing reports of lowest R-contact, suggest that the extracted R-IL (1.53 k Omega.mu m) for an unaltered channel would determine the lower limit of intrinsic R-contact even for barrierless contacts for multilayered exfoliated MoS2 FETs.

General Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters of Arbitrary Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expression is well known for indentation in elastic and in elastic-plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. Furthermore, the same expression holds true for both fast loading and unloading. These results should provide a sound basis for using the relationship for determining properties of viscoelastic solids using indentation techniques.

Effects of 'sinking in' and 'piling up' on estimating the contact area under load in indentation

The phenomena of the 'piling up' and 'sinking-in' of surface profiles in conical indentation in elastic-plastic solids with work hardening are studied using dimensional and finite-element analysis. The degree of sinking in and piling up is shown to depend on the ratio of the initial yield strength Y to Young's modulus E and on the work-hardening exponent n. The widely used procedure proposed by Oliver and Pharr for estimating contact depth is then evaluated systematically. By comparing the contact depth obtained directly from finite-element calculations with that obtained from the initial unloading slope using the Oliver-Pharr procedure, the applicability of the procedure is discussed.

Adhesion Elastic Contact and Hysteresis Effect

In this paper, we study the relationship between the pull-off force and the transition parameter (or Tabor number) as well as the variation of the pull-off radius with the transition parameter in the adhesion elastic contact. Hysteresis models are presented to describe the contact radius as a function of external loads in loading and unloading processes. Among these models, we verified the hysteresis model from Johnson{Kendall{Roberts theory, based on which the calculated results are in good agreement with experimental ones.

Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expres

Comparison of Various Adhesion Contact Theories and the Influence of Dimensionless Load Parameter

Three models, JKR (Johnson, Kendall and Roberts), DMT (Derjaguin, Muller, and Toporov) andMD (Maugis-Dugdale),are compared with the Hertz model in dealing with nano-contact problems. It has been shown that both the dimensionless load parameter, P D P=.1/4