14 mm. Wrap Angle= 2.2. x=0. Tens. x=l GAP. 23 um. 2.5 mm. 3.2 mm. x=x x=x TE. Leading Edge Trailing Edge

Size: px
Start display at page:

Download "14 mm. Wrap Angle= 2.2. x=0. Tens. x=l GAP. 23 um. 2.5 mm. 3.2 mm. x=x x=x TE. Leading Edge Trailing Edge"

Transcription

1 1 Contact Tape Recording with a Flat Head Contour Hans F. Hinteregger and Sinan Muftu Massachusetts Institute of Technology Haystack Observatory Westford, MA 01886

2 Tape Recording with a Flat Head 2 ABSTRACT A row-bar of commercial, IBM thin-lm disk heads, with 3:75m wide, 2m high, magnetoresistive read sensors, was evaluated on linear (Metrum 96) tape drive. TiC The tape was wrapped 2 o over the hard Al 2 O 3 leading corner of the bar and not wrapped at all over the relatively soft TF trailing corner. After more than 2000 hours of running in contact at 8m=s under 90N=m tension we can place an upper limit of 4nm on head wear. Output, hence spacing, was stable and independent of the tape speed in the range of 0:5 8m=s which is proof of contact with the substrate. 3 A model of the head tape interface showed a very good match to the experimentally observed results. Air pressure becomes subambient in the tape-rowbar interface. 3 The tape is \sucked" down into contact with the row-bar surface in the indicated speed range. 3 The highest contact pressure region is near the wrapped (leading) edge. 3 Contact pressure on most of the remaining at part of the tape bearing surface is only 10% of the ambient pressure.

3 Tape Recording with a Flat Head 3 Wrap Geometry Wrap Angle= 2.2 x=0 Tens. o V x GAP 14 mm 23 um x=l 3.2 mm 2.5 mm x=x x=x TE LE Leading Edge Trailing Edge Figure 1: Row-bar wrap geometry. EXPERIMENTAL RESULTS 1. WVHS and SVHS tapes were shuttled in contact for 120 and 2000 hours, respectively. 2. For the 2000 hours case shuttling was at V = 8m=s and tape tension was at T =2:3N. 3. During the same hours of operation a ferrite head wore 10m. 4. No speed dependence was observed in the MR head output in the range 0 8m=s for the SVHS tape. 5. No change in the MR head resistance was observed during this time and no loss of output was seen either (SVHS tape recorded with 2:25f c=m to within 0:5dB).

4 Tape Recording with a Flat Head 4 WVHS tape at 4.5 fc=m ( 114kfci) with 3:75m MR V (m/s) S/EN (db) S/TN (db) :5 43:5 1 51:5 45:5 2 51:5 45:5 4 51:5 46:5 8 49:5 46:5 Table 1: Output variation with tape speed. The MR sensor outperformed a 38m ferrite head. Output signal (S) to tape noise (TN) is high enough to project the feasibility of 20; 000tpi recording at 18dB SNR with SVHS- and WVHS-like tape at 3 and 6f c=m, respectively.

5 Tape Recording with a Flat Head 5

6 Tape Recording with a Flat Head 6 Model of the Head-Tape Interface Tape Eqn.: r t = D d4 w +( dx 4 a V 2 x T x ) d2 w dx 2 (p P a ) P c =0 (a) Reyn. Eqn.: r p = d dp [ph3 dx dx (1+6 a )] h 6V d(ph) x =0 (b) dx Contact P.: P c = P max(h t 2 t ) 2 H( (h t )) (c) Spacing: h = w + (d) Disp. BC: w = w o ; dw dx = m at x =0 (e) w = 0; dw dx =0atx=L x (f) Pressure BC: P = P a at x = L LE ;L TE (g) (1) x Coordinate axis P c Contact pressure w Tape displacement P max =10MPa Contact pressure at h =0 h Head-tape spacing is t =24;48nm Asperity engagement height D(= Ec3 ) 12(1 2 ) Bending stiness H Heaviside step function T x Tape tension E(= 4GP a) Modulus of elasticity V x Tape speed (= 0:3) Poisson's ratio p Air pressure c(= 15m) Tape thickness P a (= 101:3kP a) Ambient pressure w o Tape disp., x =0 a (= 63:5nm) Molecular mean-free path m = w L x LE Tape slope, x =0 (= 18:5N sm 2 ) Air viscosity r t ;r p Residuals

7 Tape Recording with a Flat Head 7 The Wear Relation According to Archard's wear law the wear volume is proportional to the applied load; Vwear = kl W H (2) where, V wear Wear volume k Wear Coecient l Sliding distance (56; 000km) W Normal load H Hardness of the material (450GP a for Al 2 O 3 TiC) Using the apparent area of contact, A a, this equation can be transformed into the following form; wear(x) =CPc(x) (3) where, C(= ) Wear coecient chosen arbitrarily wear Wear depth Apparent contact pressure P c

8 Tape Recording with a Flat Head 8 The Numerical Solution The governing dierential equations (1.a,b) are discretized with O(x 2 ) accurate nite dierence approximations. The residuals are linearized using the Newton-Raphson approach. Equations are solved iteratively; The tape equation [J t ]fwg I t+1 = fr t g I t (4) The Reynolds equation [Jp]fpg I p+1 = frpg I p (5) where, I t ; I p Iteration counters for the tape and pressure solutions fr t g I t ; fr p g I p The residuals dened by (1.a,b) [J t ]; [J p ] The non-linear Jacobian matrices for the tape and Reynolds eqns. fwg I t+1 ; fpgi p+1 The correction vectors Coupling is provided by a relaxation coecient, C relax. The tape equation is updated at the end of each iteration step; fwg (I t+1) = fwg (I t) + C relax fwg (I t+1) (6) The relaxation coecient is set initially to C relax =0:8 1 and it is cut by 1=2every time the Euclidean norm, l 2, of the tape residual grows; l (I t) 2 (w) = 0 B@ NN X k=1 r 2 t k 1 CA 1=2 (7)

9 Tape Recording with a Flat Head 9 The Solution Algorithm 0. Start with\initial conditions" 1. Start wear iterations I w = I w +1 I t =0 2. Start tape iterations, I t = I t +1 (a) If h 0 in the contact zone then, Solve Reynolds eqn. using NR-method I p =0 i. Start pressure iterations, I p = I P +1 ii. Solve [J p ]fpg I p+1 = fr p g I p iii. Calculate the new residual, r p iv. If r p < goto step 2.(b), otherwise goto step 2.(a)i. (b) Calculate the contact pressure (c) Solve [J t ]fwg I t+1 = fr t g I t (d) Calculate tape eqn. residual, r t (e) Update fwg I t+1, fhg I t+1 and check C relax (f) If r t < goto step 3 otherwise goto step 2 3. Calculate wear amount, wear 4. Calculate the new head shape, I w+1 = I w + wear 5. Calculate the new spacing, h 6. If I w limit is reached then stop otherwise goto step 1

10 Tape Recording with a Flat Head 10

11 Tape Recording with a Flat Head 11

12 Tape Recording with a Flat Head 12

13 Tape Recording with a Flat Head 13 The Edge Wear and Corresponding Pressures 0 V=8 m/s, Asperity Height=48 nm Wear depth (nm) Wear Iter.=100 Wear Iter.=500 Wear It.= Distance Along the Row-Bar (m) Air Pressure [(p-pa)/pa] V=8 m/s, Asperity Height=48 nm Wear Iter.=100 Wear Iter.=500 Wear Iter.= Distance Along the Row-Bar (m) Contact Pressure [Pc/Pa] Wear Iter.=100 Wear Iter.=500 Wear Iter.= Distance Along the Row-Bar (m)

14 Tape Recording with a Flat Head 14 Spacing on the Gap as a Function of Tape Speed T x =2:2N,Leading Edge Worn to 300 Its. 49 Asperity Size=48 nm, Wrap Ang.= 2.2 deg. Spacing, h, at the gap (nm) Tape Speed, V, (m/s) Total Force/Width (N/m) Asperity Size=48 nm, Wrap Ang.= 2.2 deg. Total Contact Force Total Air Bearing Force Total Force Due to Bending Tape Speed, V, (m/s)

15 Tape Recording with a Flat Head 15 The Eect of Tape Tension at V x =8m=s Leading Edge Worn to 300 Its. (Exp. Value) Head-tape spacing at the gap (nm) V=8 m/s, Asperity Size=48 nm, Wrap Ang.= 2.2 deg Tape Tension, Tx (N/m) Total Force/Width (N/m) V=8 m/s, Asperity Size=48 nm, Wrap Ang.= 2.2 deg. Total Contact Force Total Air Bearing Force Total Force due to Bending Tape Tension, Tx (N/m)

16 Tape Recording with a Flat Head 16 SUMMARY AND CONCLUSIONS As shown here by experiment and theory, contact between a tape and a at bearing surface can be maintained at high speeds, due the sub-ambient air pressure in the interface. Wear of the wrapped corner does not lead to ying for a wide range of tape speed and tension provided that the leading edge is worn to equilibrium. The necessity of incorporating the Reynolds equation into the wear calculations was shown. We have seen no experimental evidence of signicant wear on the at part of the Al 2 O 3 TiC substrate, nor of the softer layers near the trailing corner. The latter area which contains the TF heads is protected from wear because of the low contact pressure due to the combination of no-wrap and as-lapped TF recession. FUTURE WORK Test and model symmetrically-wrapped bi-directional contact congurations. These can be produced by capping a row-bar with substrate material or perhaps by just reducing the wrap angle to 0:5 o.

Sinan Müftü Associate Professor Department of Mechanical Engineering Northeastern University, 334 SN Boston, MA

Sinan Müftü Associate Professor Department of Mechanical Engineering Northeastern University, 334 SN Boston, MA TAPE MECHANICS OVER A FLAT RECORDING HEAD UNDER UNIFORM PULL-DOWN PRESSURE Sinan Müftü Associate Professor Department of Mechanical Engineering Northeastern University, 334 SN Boston, MA 115-5 Submitted

More information

The Mechanics of CMP and Post-CMP Cleaning

The Mechanics of CMP and Post-CMP Cleaning The Mechanics of CMP and Post-CMP Cleaning Sinan Müftü Ahmed Busnaina George Adams Department of Mechanical, Industrial and Manuf. Engineering Northeastern University Boston, MA 02115 Introduction Objective

More information

Traction Between a Web and a Smooth Roller. Sinan Müftü Associate Professor, Member ASME. John J. Jagodnik, Graduate Student

Traction Between a Web and a Smooth Roller. Sinan Müftü Associate Professor, Member ASME. John J. Jagodnik, Graduate Student Sinan Müftü Associate Professor, Member ASME John J. Jagodnik, Graduate Student Northeastern University Department of Mechanical Engineering Boston, MA 2115 Submitted for review in Journal of Tribology,

More information

JRConsulting Atwood CO Phone: (970) FAX: (775)

JRConsulting Atwood CO Phone: (970) FAX: (775) Bi-Directional Magnetic Read/Write Recording Head Surface Contour with Plurality of Bernoulli Pocket Cavities for Generating Very Low Media-to-Head Separations J. Robert (Bob) Cope JRConsulting Atwood

More information

As an example, the two-bar truss structure, shown in the figure below will be modelled, loaded and analyzed.

As an example, the two-bar truss structure, shown in the figure below will be modelled, loaded and analyzed. Program tr2d 1 FE program tr2d The Matlab program tr2d allows to model and analyze two-dimensional truss structures, where trusses are homogeneous and can behave nonlinear. Deformation and rotations can

More information

Engineering Mechanics. Friction in Action

Engineering Mechanics. Friction in Action Engineering Mechanics Friction in Action What is friction? Friction is a retarding force that opposes motion. Friction types: Static friction Kinetic friction Fluid friction Sources of dry friction Dry

More information

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS 6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement

More information

Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2

More information

Computational Modelling of the Surface Roughness Effects on the Thermal-elastohydrodynamic Lubrication Problem

Computational Modelling of the Surface Roughness Effects on the Thermal-elastohydrodynamic Lubrication Problem Proceedings of the International Conference on Heat Transfer and Fluid Flow Prague, Czech Republic, August 11-12, 2014 Paper No. 192 Computational Modelling of the Surface Roughness Effects on the Thermal-elastohydrodynamic

More information

U U B U P x

U U B U P x Smooth Quantum Hydrodynamic Model Simulation of the Resonant Tunneling Diode Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-184 Abstract Smooth

More information

CRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th

CRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA

More information

Mechanical Properties of Materials

Mechanical Properties of Materials Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of

More information

Heat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is

Heat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi

More information

197 1st Avenue, Suite 120, Needham MA Tel Fax

197 1st Avenue, Suite 120, Needham MA Tel Fax 197 1st Avenue, Suite 120, Needham MA 02494 Tel 781-444-2250 Fax 781-444-2251 USinfo@csm-instruments.com www.csm-instruments.com //// T 09 113 Wear and Friction Analysis of Thin Coatings An in-depth study

More information

A CONTACT-MECHANICS BASED MODEL FOR DISHING AND EROSION IN

A CONTACT-MECHANICS BASED MODEL FOR DISHING AND EROSION IN Mat. Res. Soc. Symp. Proc. Vol. 671 001 Materials Research Society A CONTACT-MECHANICS BASED MODEL FOR DISHING AND EROSION IN CHEMICAL-MECHANICAL POLISHING Joost J. Vlassak Division of Engineering and

More information

Modelling mechanical systems with nite elements

Modelling mechanical systems with nite elements Modelling mechanical systems with nite elements Bastian Kanning February 26, 2010 Outline 1 Introduction 2 Virtual displacements 3 Getting started The ansatz functions 4 Mass and stiness elementwise 5

More information

General elastic beam with an elastic foundation

General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

1 Nonlinear deformation

1 Nonlinear deformation NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected

More information

MECHANICAL PROPERTIES OF HYDROGEL USING NANOINDENTATION

MECHANICAL PROPERTIES OF HYDROGEL USING NANOINDENTATION MECHANICAL PROPERTIES OF HYDROGEL USING NANOINDENTATION Prepared by Duanjie Li, PhD & Jorge Ramirez 6 Morgan, Ste156, Irvine CA 9618 P: 949.461.99 F: 949.461.93 nanovea.com Today's standard for tomorrow's

More information

Sensitivity and Reliability Analysis of Nonlinear Frame Structures

Sensitivity and Reliability Analysis of Nonlinear Frame Structures Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,

More information

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,

More information

Notes on Rubber Friction

Notes on Rubber Friction Notes on Rubber Friction 2011 A G Plint Laws of Friction: In dry sliding between a given pair of materials under steady conditions, the coefficient of friction may be almost constant. This is the basis

More information

Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN

Transactions on Engineering Sciences vol 14, 1997 WIT Press,  ISSN On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,

More information

This chapter introduces the description of the surface interaction mechanism based on the friction, wear and excavation laws.

This chapter introduces the description of the surface interaction mechanism based on the friction, wear and excavation laws. Chapter 5 Surface interaction 5.1 Introduction This chapter introduces the description of the surface interaction mechanism based on the friction, wear and excavation laws. Whenever two solids touch each

More information

Analysis of contact deformation between a coated flat plate and a sphere and its practical application

Analysis of contact deformation between a coated flat plate and a sphere and its practical application Computer Methods and Experimental Measurements for Surface Effects and Contact Mechanics VII 307 Analysis of contact deformation between a coated flat plate and a sphere and its practical application T.

More information

Measurement and modelling of contact stiffness

Measurement and modelling of contact stiffness Measurement and modelling of contact stiffness D. Nowell, University of Oxford, UK Joints workshop, Chicago, 16/8/12 Difficulties in modelling contacts In general, the normal and tangential stiffnesses

More information

1. Tasks of designing

1. Tasks of designing 1 Lecture #18(14) Designing calculation of cross section of a highly aspect ratio wing Plan: 1 Tass of designing Distribution of shear force between wing spars Computation of the elastic center 4 Distribution

More information

[5] Stress and Strain

[5] Stress and Strain [5] Stress and Strain Page 1 of 34 [5] Stress and Strain [5.1] Internal Stress of Solids [5.2] Design of Simple Connections (will not be covered in class) [5.3] Deformation and Strain [5.4] Hooke s Law

More information

PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation

PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In

More information

E1. Column subjected to bending and axial action, Non-dissipative material

E1. Column subjected to bending and axial action, Non-dissipative material E1. Column subjected to bending and axial action, Non-dissipative material Problem definition In this exercise, calculate roughly the maximum load, due to the material behavior, of a column structure made

More information

For an imposed stress history consisting of a rapidly applied step-function jump in

For an imposed stress history consisting of a rapidly applied step-function jump in Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For

More information

Centro de Processamento de Dados, Universidade Federal do Rio Grande do Sul,

Centro de Processamento de Dados, Universidade Federal do Rio Grande do Sul, A COMPARISON OF ACCELERATION TECHNIQUES APPLIED TO THE METHOD RUDNEI DIAS DA CUNHA Computing Laboratory, University of Kent at Canterbury, U.K. Centro de Processamento de Dados, Universidade Federal do

More information

INFLUENCE OF NORMAL FORCE AND HUMIDITY ON FRICTION AND WEAR OF UNLUBRICATED STEEL/ STEEL COUPLES

INFLUENCE OF NORMAL FORCE AND HUMIDITY ON FRICTION AND WEAR OF UNLUBRICATED STEEL/ STEEL COUPLES INFLUENCE OF NORMAL FORCE AND HUMIDITY ON FRICTION AND WEAR OF UNLUBRICATED STEEL/ STEEL COUPLES D. KLAFFKE Federal Institute for Materials Research and Testing (BAM), Lab. VIII.2, Unter den Eichen 87,

More information

MODELING THE EFFECTIVE ELASTIC MODULUS OF RC BEAMS EXPOSED TO FIRE

MODELING THE EFFECTIVE ELASTIC MODULUS OF RC BEAMS EXPOSED TO FIRE Journal of Marine Science and Technology, Vol., No., pp. -8 () MODELING THE EFFECTIVE ELASTIC MODULUS OF RC BEAMS EXPOSED TO FIRE Jui-Hsiang Hsu*, ***, Cherng-Shing Lin**, and Chang-Bin Huang*** Key words:

More information

Transition from Boundary Lubrication to Hydrodynamic Lubrication of Slider Bearings

Transition from Boundary Lubrication to Hydrodynamic Lubrication of Slider Bearings R. C. Tseng F. E. Talke Transition from Boundary Lubrication to Hydrodynamic Lubrication of Slider Bearings Abstract: The transition from boundary lubrication to fully hydrodynamic lubrication is investigated

More information

quantum semiconductor devices. The model is valid to all orders of h and to rst order in the classical potential energy. The smooth QHD equations have

quantum semiconductor devices. The model is valid to all orders of h and to rst order in the classical potential energy. The smooth QHD equations have Numerical Simulation of the Smooth Quantum Hydrodynamic Model for Semiconductor Devices Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-84 Abstract

More information

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free

More information

Table 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad

Table 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad Chapter 1 Introduction to Boundary element Method - 1D Example For reference: Hong-Ki Hong and Jeng-Tzong Chen, Boundary Element Method, Chapter 1 Introduction to Boundary element Method - 1D Example,

More information

Optimization Methods. Lecture 23: Semidenite Optimization

Optimization Methods. Lecture 23: Semidenite Optimization 5.93 Optimization Methods Lecture 23: Semidenite Optimization Outline. Preliminaries Slide 2. SDO 3. Duality 4. SDO Modeling Power 5. Barrier lgorithm for SDO 2 Preliminaries symmetric matrix is positive

More information

Characterisation Programme Polymer Multi-scale Properties Industrial Advisory Group 22 nd April 2008

Characterisation Programme Polymer Multi-scale Properties Industrial Advisory Group 22 nd April 2008 Characterisation Programme 6-9 Polymer Multi-scale Properties Industrial Advisory Group nd April 8 SE: Improved Design and Manufacture of Polymeric Coatings Through the Provision of Dynamic Nano-indentation

More information

Sinan Muftu. Submitted in Partial Fulllment. of the. Requirements for the Degree. Doctor of Philosophy. University ofrochester. Rochester, New York

Sinan Muftu. Submitted in Partial Fulllment. of the. Requirements for the Degree. Doctor of Philosophy. University ofrochester. Rochester, New York The Transient Foil Bearing Problem in Magnetic Recording by Sinan Muftu Submitted in Partial Fulllment of the Requirements for the Degree Doctor of Philosophy Supervised by Professor Richard C. Benson

More information

Design of a hydrostatic symmetric-pad bearing with the membrane-type restrictor

Design of a hydrostatic symmetric-pad bearing with the membrane-type restrictor Design of a hydrostatic symmetric-pad bearing with the membrane-type restrictor Professor: Shih-Chieh Lin Manufacturing and Production System Lab Dept. of Power Mechanical Engineering, National Tsing Hua

More information

Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method

Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July 11-15, 2016 ICCFD9-0113 Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume

More information

On the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method

On the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846

More information

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,

More information

Biaxial Analysis of General Shaped Base Plates

Biaxial Analysis of General Shaped Base Plates Biaxial Analysis of General Shaped Base Plates R. GONZALO ORELLANA 1 Summary: A linear model is used for the contact stresses calculation between a steel base plate and a concrete foundation. It is also

More information

Linear stochastic approximation driven by slowly varying Markov chains

Linear stochastic approximation driven by slowly varying Markov chains Available online at www.sciencedirect.com Systems & Control Letters 50 2003 95 102 www.elsevier.com/locate/sysconle Linear stochastic approximation driven by slowly varying Marov chains Viay R. Konda,

More information

Performance of Low Density Parity Check Codes. as a Function of Actual and Assumed Noise Levels. David J.C. MacKay & Christopher P.

Performance of Low Density Parity Check Codes. as a Function of Actual and Assumed Noise Levels. David J.C. MacKay & Christopher P. Performance of Low Density Parity Check Codes as a Function of Actual and Assumed Noise Levels David J.C. MacKay & Christopher P. Hesketh Cavendish Laboratory, Cambridge, CB3 HE, United Kingdom. mackay@mrao.cam.ac.uk

More information

Boundary and Mixed Lubrication Friction Modeling under Forming Process Conditions

Boundary and Mixed Lubrication Friction Modeling under Forming Process Conditions Boundary and Mixed Lubrication Friction Modeling under Forming Process Conditions V.T. Meinders, J. Hol, and A.H. van den Boogaard University of Twente, Faculty of Engineering Technology, P.O.Box 217,

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

Quiescent Steady State (DC) Analysis The Newton-Raphson Method

Quiescent Steady State (DC) Analysis The Newton-Raphson Method Quiescent Steady State (DC) Analysis The Newton-Raphson Method J. Roychowdhury, University of California at Berkeley Slide 1 Solving the System's DAEs DAEs: many types of solutions useful DC steady state:

More information

Figure 43. Some common mechanical systems involving contact.

Figure 43. Some common mechanical systems involving contact. 33 Demonstration: experimental surface measurement ADE PhaseShift Whitelight Interferometer Surface measurement Surface characterization - Probability density function - Statistical analyses - Autocorrelation

More information

Numerical Limit Analysis of Rigid Plastic Structures

Numerical Limit Analysis of Rigid Plastic Structures Numerical Limit Analysis of Rigid Plastic Structures Hector Andres Tinoco Navarro Institute of Technology and Innovation University of Southern Denmark Supervisors: Prof. Ph.D. Linh Cao Hoang University

More information

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION Supplementary Information for Manuscript: Nanoscale wear as a stress-assisted chemical reaction Supplementary Methods For each wear increment, the diamond indenter was slid laterally relative to the silicon

More information

Lecture 7. Pile Analysis

Lecture 7. Pile Analysis Lecture 7 14.5 Release Pile Analysis 2012 ANSYS, Inc. February 9, 2013 1 Release 14.5 Pile definition in Mechanical - There are a number of methods that can be used to analyze piled foundations in ANSYS

More information

ENGN 2290: Plasticity Computational plasticity in Abaqus

ENGN 2290: Plasticity Computational plasticity in Abaqus ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element

More information

quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner

quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size

More information

Advanced Mechanical Principles

Advanced Mechanical Principles Unit 36: Unit code Advanced Mechanical Principles R/615/1504 Unit level 5 Credit value 15 Introduction A mechanical engineer is required to have an advanced knowledge of most of the machinery used within

More information

Elements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004

Elements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004 Elements of Polymer Structure and Viscoelasticity David M. Parks Mechanics and Materials II 2.002 February 18, 2004 Outline Elements of polymer structure Linear vs. branched; Vinyl polymers and substitutions

More information

PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.

PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9. PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional

More information

PILE SOIL INTERACTION MOMENT AREA METHOD

PILE SOIL INTERACTION MOMENT AREA METHOD Pile IGC Soil 2009, Interaction Moment Guntur, INDIA Area Method PILE SOIL INTERACTION MOMENT AREA METHOD D.M. Dewaikar Professor, Department of Civil Engineering, IIT Bombay, Mumbai 400 076, India. E-mail:

More information

Simulating Wear in Disc Brakes

Simulating Wear in Disc Brakes Simulating Wear in Disc Brakes Nagi H. Elabbasi, Matthew J. Hancock, and Stuart B. Brown, Massachusetts 47A Kearney Road Needham, Massachusetts 02494 phone 781.433.0433 multiphysics@veryst.com www.veryst.com

More information

BUTT SPLICE HINGING. KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated

BUTT SPLICE HINGING. KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated BUTT SPLICE HINGING BY KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated Introduction Splicing is a process used to join the tail of an expiring roll to the start

More information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how

More information

Strength of Materials Prof. S.K.Bhattacharya Dept. of Civil Engineering, I.I.T., Kharagpur Lecture No.26 Stresses in Beams-I

Strength of Materials Prof. S.K.Bhattacharya Dept. of Civil Engineering, I.I.T., Kharagpur Lecture No.26 Stresses in Beams-I Strength of Materials Prof. S.K.Bhattacharya Dept. of Civil Engineering, I.I.T., Kharagpur Lecture No.26 Stresses in Beams-I Welcome to the first lesson of the 6th module which is on Stresses in Beams

More information

Measuring State Parameters of the Atmosphere

Measuring State Parameters of the Atmosphere Measuring State Parameters of the Atmosphere Some Applications of Atmospheric Thermodynamics Earth Observing Laboratory, NCAR IDEAS-4 Tutorial Introduction Goals of This Presentation Present two complementary

More information

Introductory guide to measuring the mechanical properties of nanoobjects/particles

Introductory guide to measuring the mechanical properties of nanoobjects/particles Jeremias Seppä MIKES Metrology, VTT Technical Research Centre of Finland Ltd P.O. Box 1000, FI-02044 VTT, Finland Contents: AFM Cantilever calibration F-d curves and cantilever bending Hitting the particles

More information

Strain Gages. Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, Shear Modulus, (S) N/m 2

Strain Gages. Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, Shear Modulus, (S) N/m 2 When you bend a piece of metal, the Strain Gages Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, 1979 Material Young's Modulus, (E) 10 11 N/m 2 Shear Modulus,

More information

Creep Compliance Analysis Technique for the Flattened Indirect Tension Test of Asphalt Concrete

Creep Compliance Analysis Technique for the Flattened Indirect Tension Test of Asphalt Concrete Creep Compliance Analysis Technique for the Flattened Indirect Tension Test of Asphalt Concrete Eshan V. Dave Andrew F. Braham Prof. William G. Buttlar Prof. Glaucio H. Paulino CONCREEP8, Ise-Shima, Japan

More information

How to measure the shear viscosity properly?

How to measure the shear viscosity properly? testxpo Fachmesse für Prüftechnik 10.-13.10.2016 How to measure the shear viscosity properly? M p v Rotation Capillary Torsten Remmler, Malvern Instruments Outline How is the Shear Viscosity defined? Principle

More information

Finite element analysis of membrane wrinkling

Finite element analysis of membrane wrinkling INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 50:1017 1038 Finite element analysis of membrane wrinkling K. Lu, M. Accorsi ; and J. Leonard Department of Civil

More information

Numerical and Experimental Studies on Thermoforming Process. Sogang University

Numerical and Experimental Studies on Thermoforming Process. Sogang University Numerical and Experimental Studies on Thermoforming Process Thermoforming Process Hot plate Atmosphere Seal Mold Air on Air on Vacuum or atmosphere Introduction Thermoforming Process Advantage Low forming

More information

INTRODUCTION TO STRAIN

INTRODUCTION TO STRAIN SIMPLE STRAIN INTRODUCTION TO STRAIN In general terms, Strain is a geometric quantity that measures the deformation of a body. There are two types of strain: normal strain: characterizes dimensional changes,

More information

Advanced Friction Modeling in Sheet Metal Forming

Advanced Friction Modeling in Sheet Metal Forming Advanced Friction Modeling in Sheet Metal Forming J.Hol 1,a, M.V. Cid Alfaro 2, T. Meinders 3, J. Huétink 3 1 Materials innovation institute (M2i), P.O. box 58, 26 GA Delft, The Netherlands 2 Tata Steel

More information

Analytical solution for polish-rate decay in chemical mechanical polishing

Analytical solution for polish-rate decay in chemical mechanical polishing J Eng Math DOI 10.1007/s10665-010-9369-9 LETTER TO THE EDITOR Analytical solution for polish-rate decay in chemical mechanical polishing Hong Shi Terry A. Ring Received: 17 August 2009 / Accepted: 15 March

More information

A new nite-element formulation for electromechanical boundary value problems

A new nite-element formulation for electromechanical boundary value problems INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2002; 55:613 628 (DOI: 10.1002/nme.518) A new nite-element formulation for electromechanical boundary value problems

More information

Supplementary Figures

Supplementary Figures Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image

More information

1 Outline Part I: Linear Programming (LP) Interior-Point Approach 1. Simplex Approach Comparison Part II: Semidenite Programming (SDP) Concludin

1 Outline Part I: Linear Programming (LP) Interior-Point Approach 1. Simplex Approach Comparison Part II: Semidenite Programming (SDP) Concludin Sensitivity Analysis in LP and SDP Using Interior-Point Methods E. Alper Yldrm School of Operations Research and Industrial Engineering Cornell University Ithaca, NY joint with Michael J. Todd INFORMS

More information

Lab 01: Harmonic Motion I. Theory: Three experiments. The first we measure the oscillation of a spring, the second of a rubber band (non-linear).

Lab 01: Harmonic Motion I. Theory: Three experiments. The first we measure the oscillation of a spring, the second of a rubber band (non-linear). Dr. W. Pezzaglia Physics 8C Lab, Spring 04 Page Las Positas College Lab # Harmonic Motion 04Jan3 Lab 0: Harmonic Motion I. Theory: Three experiments. The first we measure the oscillation of a spring, the

More information

Geometric Stiffness Effects in 2D and 3D Frames

Geometric Stiffness Effects in 2D and 3D Frames Geometric Stiffness Effects in D and 3D Frames CEE 41. Matrix Structural Analsis Department of Civil and Environmental Engineering Duke Universit Henri Gavin Fall, 1 In situations in which deformations

More information

Hakwoon Kim Gunhee Jang Sanghoon Lee. 1 Introduction

Hakwoon Kim Gunhee Jang Sanghoon Lee. 1 Introduction Microsyst Technol (2011) 17:749 759 DOI 10.1007/s00542-010-1188-4 TECHNICAL PAPER Complete determination of the dynamic coefficients of coupled journal and thrust bearings considering five degrees of freedom

More information

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Influential Factors on Adhesion between Wheel and Rail under Wet Conditions

Influential Factors on Adhesion between Wheel and Rail under Wet Conditions Influential Factors on Adhesion between Wheel and Rail under Wet Conditions H. Chen, M. Ishida, 2 T. Nakahara Railway Technical Research Institute, Tokyo, Japan ; Tokyo Institute of Technology, Tokyo,

More information

The Rotating Spring. ω Rotating spring. Spring at rest. x L. x' L' Jacob Gade Koefoed

The Rotating Spring. ω Rotating spring. Spring at rest. x L. x' L' Jacob Gade Koefoed The Rotating Spring A spring of mass m, which is uniformly distributed along the spring, spring constant k and unstretched length L is rotated about an ais perpendicular to the spring at a constant angular

More information

Mechanical Principles

Mechanical Principles Unit 4: Mechanical Principles Unit code: F/601/1450 QCF level: 5 Credit value: 15 Aim This unit aims to develop learners understanding of an extended range of mechanical principles that underpin the design

More information

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield

More information

3.5 Reinforced Concrete Section Properties

3.5 Reinforced Concrete Section Properties CHAPER 3: Reinforced Concrete Slabs and Beams 3.5 Reinforced Concrete Section Properties Description his application calculates gross section moment of inertia neglecting reinforcement, moment of inertia

More information

Supplementary Material

Supplementary Material Mangili et al. Supplementary Material 2 A. Evaluation of substrate Young modulus from AFM measurements 3 4 5 6 7 8 Using the experimental correlations between force and deformation from AFM measurements,

More information

Reglerteknik, TNG028. Lecture 1. Anna Lombardi

Reglerteknik, TNG028. Lecture 1. Anna Lombardi Reglerteknik, TNG028 Lecture 1 Anna Lombardi Today lecture We will try to answer the following questions: What is automatic control? Where can we nd automatic control? Why do we need automatic control?

More information

Sliding Bearings. Fig.(1) (a) Full-journal bearing and (b) partial-journal bearing

Sliding Bearings. Fig.(1) (a) Full-journal bearing and (b) partial-journal bearing Sliding Bearings The goal of a bearing is to provide relative positioning and rotational freedom while transmitting a load between two parts, commonly a shaft and its housing. The object of lubrication

More information

2 marks Questions and Answers

2 marks Questions and Answers 1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and

More information

Chapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship

Chapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction

More information

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST

More information

A new computational method for threaded connection stiffness

A new computational method for threaded connection stiffness Research Article A new computational method for threaded connection stiffness Advances in Mechanical Engineering 2016, Vol. 8(12) 1 9 Ó The Author(s) 2016 DOI: 10.1177/1687814016682653 aime.sagepub.com

More information

Modeling of a Circular Plate with Piezoelectric Actuators of Arbitrary Shape

Modeling of a Circular Plate with Piezoelectric Actuators of Arbitrary Shape Vol. 13 013 ACTA PHYSICA POLONICA A No. 6 Acoustic Biomedical Engineering Modeling of a Circular Plate with Piezoelectric Actuators of Arbitrary Shape M. Wiciak a, R. Trojanowski b a Cracow University

More information

1 Introduction IPICSE-2016

1 Introduction IPICSE-2016 (06) DOI: 0.05/ matecconf/06860006 IPICSE-06 Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method

More information

CE 102: Engineering Mechanics. Minimum Potential Energy

CE 102: Engineering Mechanics. Minimum Potential Energy CE 10: Engineering Mechanics Minimum Potential Energy Work of a Force During a Finite Displacement Work of a force corresponding to an infinitesimal displacement, Work of a force corresponding to a finite

More information