Adhesion of micro-cantilevers subjected to mechanical point loading: modeling and experiments

This paper presents experimental and theoretical results that characterize the adhesion of MEMS cantilevers by means of mechanical actuation. Micro-cantilever beams are loaded at various locations along the freestanding portion of the beam using a nanoindenter. Transitions between three equilibrium configurations (freestanding, arc-shaped, and s-shaped beams) and the response to cyclic loading are studied experimentally. The resulting mechanical response is used to estimate the interface adhesion energy (using theoretical models), and to quantify the energy dissipated during cyclic loading. The experiments reveal interesting behaviors related to adhesion: (i) path dependence during mechanical loading of adhered beams, (ii) history dependence of interfacial adhesion energy during repeated loading, and (iii) energy dissipation during cyclic loading, which scales roughly with estimated cyclic changes in the size of the adhered regions. The experimental results are interpreted in the context of elementary fracture-based adhesion and contact models, and briefly discussed in terms of their implications regarding the nature of adhesion and future modeling to establish adhesion mechanisms.     

Nonlinear convection of a viscoelastic fluid with inclined temperature gradient

The energy method is used to analyze the viscoelastic fluid convection problem in a thin horizontal layer, subjected to an applied inclined temperature gradient. The boundaries are considered to be rigid and perfectly conducting. Both linear and nonlinear stability analyses are carried out. The eigenvalue problem is solved by the Chebyshev Tau-QZ method and comparisons are reported between the results of the linear theory and energy stability theory.Received: 12 March 2004, Accepted: 19 April 2004, Published online: 17 September 2004PACS: 47.20 Ky, 47-27 Te, 83.60 Wc Correspondence to: P.N. Kaloni     

基于黏弹性的微注射成型流动模拟

应用有限元方法研究了微注射成型中瞬态、可压缩、非牛顿熔体流动的黏弹性对流动前沿及流动平衡的影响。基于Phan-Thien-Tanner模型建立了熔体流动的本构方程,利用Hele-Shaw假设和简化建立了瞬态、可压缩、非牛顿熔体流动的连续性方程、动量方程、能量方程;为了有效地描述微注射成型的尺寸效应,采用了边界滑移和表面张力边界条件。通过分部积分和待定系数法导出了带有边界信息的变分方程和求解应力分量的半解析公式,构造了有限元离散求解及超松驰迭代算法。模拟结果表明:熔体的黏弹性对浇口附近的压力和后续的熔体流动前沿有重要影响;与黏性模型相比,黏弹性模型可以控制模拟压力的快速增长,减少不同型腔之间的充填差异,与短射实验结果也更吻合。     

Regulation of Cell Adhesion Free Energy by External Sliding Forces

Cells experience a variety of forces and external velocities in vivo. While a number of continuum based models have addressed cell substrate interactions in these dynamic environments, a molecular picture of adhesion under external sliding forces and velocities is lacking. Using a molecular thermodynamic model, that incorporates entropic and steric repulsions, molecular conformations and constraints, penetrable substrates and explicit binding interactions, we study the effect of external sliding velocities on cell adhesion. We map the free energy landscapes under a broad range of external forces, binding energies and receptor surface coverages. Our calculations predict the regimes where free energy landscapes become resistant to external forces. Our model shows good agreement with experimental studies and lays out a scalable framework for analyzing and quantifying a broad spectrum of in vivo and in vitro adhesion studies.     

Two-dimensional planar flow of a viscoelastic plastic medium

The present study is concerned with finite element simulation of the planar entry flow of a viscoelastic plastic medium exhibiting yield stress. The numerical scheme is based on the Galerkin formulation. Flow experiments are carried out on a carbon black filled rubber compound. Steady-state pressure drops are measured on two sets of contraction or expansion dies having different lengths and a constant contraction or expansion ratio of 4:1 with entrance angles of 90, 45 and 15 degrees. The predicted and measured pressure drops are compared. The predicted results indicate that expansion flow has always a higher pressure drop than contraction flow. This prediction is in agreement with experimental data only at low flow rates, but not at high flow rates. The latter disagreement is possibly an indication that the assumption of fully-developed flow in the upstream and downstream regions is not realistic at high flow rates, even for the large length-to-thickness ratio channels employed. The evolution of the velocity, shear stress, and normal stress fields in the contraction or expansion flow and the location of pseudo-yield surfaces are also calculated.     

Adhesion and friction of an elastic half-space in contact with a slightly wavy rigid surface

The paper deals with the estimation of the pressure distribution, the shape of contact and the friction force at the interface of a flat soft elastic solid moving on a rigid half-space with a slightly wavy surface. In this case an unsymmetrical contact is considered and justified with the adhesion hysteresis. For soft solids as rubber and polymers the friction originates mainly from two different contributions: the internal friction due to the viscoelastic properties of the bulk and the adhesive processes at the interface of the two solids. In the paper the authors focus on the latter contribution to friction. It is known, indeed, that for soft solids, as rubber, the adhesion hysteresis is, at least qualitatively, related to friction: the larger the adhesion hysteresis the larger the friction. Several mechanisms may govern the adhesion hysteresis, such as the interdigitation process between the polymer chains, the local small-scale viscoelasticity or the local elastic instabilities. In the paper the authors propose a model to link, from the continuum mechanics point of view, the friction to the adhesion hysteresis. A simple one-length scale roughness model is considered having a sinusoidal profile. For partial contact conditions the detached zone is taken to be a mode I propagating crack. Due to the adhesion hysteresis, the crack is affected by two different values of the strain energy release rate at the advancing and receding edges respectively. As a result, an unsymmetrical contact and a friction force arise. Additionally, the stability of the equilibrium configurations is discussed and the adherence force for jumping out of contact and the critical load for snapping into full contact are estimated.     

Adhesion between elastic cylinders based on the double-Hertz model

A cohesive zone model for two-dimensional adhesive contact between elastic cylinders is developed by extending the double-Hertz model of Greenwood and Johnson (1998). In this model, the adhesive force within the cohesive zone is described by the difference between two Hertzian pressure distributions of different contact widths. Closed-form analytical solutions are obtained for the interfacial traction, deformation field and the equilibrium relation among applied load, contact half-width and the size of cohesive zone. Based on these results, a complete transition between the JKR and the Hertz type contact models is captured by defining a dimensionless transition parameter μ, which governs the range of applicability of different models. The proposed model and the corresponding analytical results can serve as an alternative cohesive zone solution to the two-dimensional adhesive cylindrical contact.     

Linear viscoelastic model for bending and torsional modes in fluid membranes

This paper develops and applies a linear viscoelastic model for bending and torsional modes of fluid membranes, based on the nonlinear Cosserat surface fluid model. The linearized fluid membrane model in spherical and cylindrical geometries is shown to decouple bending and torsional viscoelastic modes. It is found that solutions of the membrane viscoelastic model to small-amplitude oscillatory bending and torsion allows for the measurement of the bending and the torsion viscosity. The model and its potential in characterizing the bending and torsion viscoelasticity of membranes complements the on-going efforts to establish the role of curvature in dissipative process of biological membranes.     

Stability of inhomogeneously aging viscoelastic bodies

Moscow Automechanical Institute. Translated from Prikladnaya Mekhanika, Vol. 25, No. 8, pp. 39–44, August, 1989.     

Run-up and spin-up in a viscoelastic fluid. IV

Three problems are discussed — run-up when the fluid is contained between infinite parallel plates and longitudinal run-up and spin-up when it is contained in an infinitely long circular cylinder. The procedure adopted for solving these problems differs from that employed in Part I of this series, where these three problems were previously discussed, and yields results for the velocity fields in quite different forms. It is similar to that used in Part III in the context of the problem of run-up between parallel plates.     

Run-up and spin-up in a viscoelastic fluid. III

The problem is discussed of run-up in an incompressible viscoelastic fluid contained between infinite parallel rigid plates which are simultaneously given equal parallel velocities. The problem is analyzed in terms of the disturbances which spread from the boundaries into the fluid and are reflected back and forth at the boundaries.     

Pressure wave propagation in liquid-filled tubes of viscoelastic material

The propagation of small-amplitude waves in a thick-walled long viscoelastic tube of variable cross-section, filled with a viscous incompressible fluid, is considered with account for wave reflection at the tube end in application to arterial pulse wave propagation. A solution is obtained in the form of expansions in a small parameter. The effect of the coefficient of wave reflection at the tube end and the wall material parameters on the fluid volume flow-rate and the tube wall displacement is investigated. It is shown that the volume flow-rate phase spectrum characteristics depend only slightly on the wall properties and can be used in clinical diagnostics for finding the reflection coefficient from pressure and flow-rate records.     

Computation of viscoelastic cable coating flows

A viscoelastic analysis is presented for model tube tooling, draw-down and combined geometry flows encountered in the cable coating industries. The work investigates the development of stress fields and studies the effect of varying entry flow stress boundary conditions. The analysis takes into account tube tooling and draw-down flow sections individually, and in combination. The flow behaviour of cable-coating grade low density polyethylene is studied assuming a viscoelastic, isothermal flow, and employing a Taylor–Petrov–Galerkin finite element scheme with an exponential Phan-Thien–Tanner constitutive model. © 1998 John Wiley & Sons, Ltd.     

Flow of viscoelastic fluid between confocal elliptic cylinders

Summary Using elliptic coordinates, the flow pattern of a fluid of grade four between two elliptic tubes is determined. A comparison between the position of the maximum of the axial velocity in the present case and in the case of two concentric circular tubes shows a basic difference. In the elliptic case the maximum is shifted towards the external wall, while in the case of concentric circular tubes the shift is in the direction of the internal wall. The secondary flow shows dissymmetry with reference to the intermediate line , which itself lies nearer to the external wall.

---

Zusammenfassung Unter Benutzung elliptischer Koordinaten wird die Strömung zwischen zwei elliptischen Rohren bestimmt. Ein Vergleich zwischen der Lage des axialen Geschwindigkeitsmaximums im vorliegenden Fall und im Fall zweier konzentrischer Kreisrohre ergibt einen grundsätzlichen Unterschied: Das Maximum ist im elliptischen Fall zur äußeren Wand hin verschoben, während die Verschiebung im Fall der konzentrischen Kreisrohre zur inneren Wand hin erfolgt. Die Sekundärströmung ist unsymmetrisch relativ zur mittleren Stromlinie , die selbst näher zur äußeren Wand liegt.

---

A planar domain representing the annular region - vector inx ₁,x ₂-plane - x _i rectangular coordinates - rectangular unit vectors - , elliptic coordinates - ₁, ₂ ellipses representing respectively the internal and external tubes - = ₂ – ₁ annular widthy = ( – ₁)/ - µ 1^st grade material constant - _i 2^nd grade material constants - _i 3^rd grade material constants - _i 4^th grade material constants - I unit tensor - T _E extra stress (T + pI) - V potential of body forces - material density = (p/) + V = –ax ₃ + () - a specific driving force - arbitrary scalar function - A _k Rivlin-Eriksen tensors - S stress scalar defined onA - t stress vector defined onA - P stress tensor defined onA - v axial velocity - v _i i ^th term in the approximation ofv - u velocity vector perpendicular to the axis 4( ₃ + ₄ + ₅ + 1/2₆) –2/µ(2 ₁ + ₂)( ₂ + ₃) - T stress tensor - p arbitrary hydrostatic pressure - u _i i ^th term in the approximation ofu - stream function definingu - _i i ^th term in the approximation of With 8 figures and 1 table     

An empirical constitutive equation of integral type for viscoelastic liquids

A 5-constant constitutive equation is proposed. The analytical form for the relaxation modulus as a function of flow conditions was chosen based on experimental data for stress-relaxation in solid polymers. The resulting formulae for the material functions in simple and oscillatory shear flow fulfil the empirical Cox-Merz rule as well as other phenomenological relations formulated by Coleman and Markowitz. The theoretical results are compared with experimental data obtained by Han for various polymer melts. Good agreement between theory and experiment is found.     

Solitary waves on the surface of a viscoelastic fluid running down an inclined plane

In this paper, we study the existence and the role of solitary waves in the finite amplitude instability of a layer of a second-order fluid flowing down an inclined plane. The layer becomes unstable for disturbances of large wavelength for a critical value of Reynolds number which decreases with increase in the viscoelastic parameter M. The long-term evolution of a disturbance with an initial cosinusoidal profile as a result of this instability reveals the existence of a train of solitary waves propagating on the free surface. A novel result of this study is that the number of solitary waves decreases with in crease in M. When surface tension is large, we use dynamical system theory to describe solitary waves in a moving frame by homoclinic trajectories of an associated ordinary differential equation.     

The probability stability of viscoelastic plates

The almost sure stability of homogeneous viscoelastic plates subjected to a random wide-band stationary in-plane load is investigated. The viscoelastic behavior of the plate is described in terms of the Boltzmann superposition principle, the relaxation kemels of which are represented by the sums of exponents. On the assumption that the in-plane load is random wide-band stationary process, sufficient conditions for almost sure stability of viscoelastic plates are obtained by the applications of Lyapunov's direct method. Project supported by the National Natural Science Foundation of China (No. 59635140) and the National Postdoctoral Foundation.     

Nonlinear dynamic buckling of a viscoelastic model under impact loading

A simple nonlinear buckling analysis is applied to a one-degree-of-freedom arch under impact loading in which viscous damping may also be included. Such a loading consists of a falling body striking centrally the joint mass of the arch in such a way that a completely plastic impact can be postulated. When there is no damping the exact dynamic buckling load for such a kind of loading-associated with an unbounded motion can be established by using a static criterion (approach). More specifically, it was shown that the dynamic buckling load corresponds to that unstable equilibrium state where the total potential energy of the system is zero. Furthermore, it was proved that the second variation of the total potential energy at the foregoing unstable equilibrium state is negative definite. This implies that the curve loading versus displacement resulting by the vanishing of the total potential energy has always a maximum on the afore mentioned unstable state. It was also found that the system may become sensitive to initial conditions. If damping is included the foregoing static criterion yields lower bound buckling estimates. These findings were verified by employing a highly efficient approximate technique as well as the numerical scheme of Runge-Kutta for solving any nonlinear initial-value problem.