Friction Properties of Surface with Circular Micro-patterns

Similar documents
EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY

A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius

Measurement and modelling of contact stiffness

Figure 43. Some common mechanical systems involving contact.

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

INTRODUCTION TO STRAIN

A Finite Element Study of the Residual Stress and Deformation in Hemispherical Contacts

Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation

A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

Influential Factors on Adhesion between Wheel and Rail under Wet Conditions

PREDICTION OF OUT-OF-PLANE FAILURE MODES IN CFRP

Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact

N = Shear stress / Shear strain

Experimental Investigation of Fully Plastic Contact of a Sphere Against a Hard Flat

Pullout Tests of Geogrids Embedded in Non-cohesive Soil

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

Introduction to Engineering Materials ENGR2000. Dr. Coates

Mechanical Properties of Materials

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

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method

TRACTION AND WEAR MECHANISMS DURING ROLL-SLIP CONTACT

[5] Stress and Strain

CONTACT MODEL FOR A ROUGH SURFACE

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

BE Semester- I ( ) Question Bank (MECHANICS OF SOLIDS)

A novel technique of friction and material property measurement by tip test in cold forging

Determination of Friction Coefficient of a Press-Fit Pin in Thin Plating

Numerical simulation the bottom structures. grounding test by LS-DYNA

Experiment Two (2) Torsional testing of Circular Shafts

Design Optimization of Stress Relief Grooves in Lever Guide of Pressure Vessel for Food Processing

Stress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole

Chapter Two: Mechanical Properties of materials

Design of a fastener based on negative Poisson's ratio foam adapted from

A multiscale framework for lubrication analysis of bearings with textured surface

MECHANICS OF MATERIALS

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion

MEASUREMENTS OF THE COLLISION PROPERTIES OF SMALL SPHERES. Adam Lorenz and Michel Y. Louge. Sibley School of Mechanical and Aerospace Engineering

Discrete Element Modelling of a Reinforced Concrete Structure

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

Notes on Rubber Friction

Bending Load & Calibration Module

STRESS ANALYSIS AND STRENGTH EVALUATION OF SCARF ADHESIVE JOINTS SUBJECTED TO STATIC TENSILE LOADINGS. Graduate School of Mechanical Engineering

Finite-Element Analysis of Stress Concentration in ASTM D 638 Tension Specimens

1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor.

Temperature analysis of a pin-on-disc tribology test using experimental and numerical approaches

Contact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

ISHIK UNIVERSITY DEPARTMENT OF MECHATRONICS ENGINEERING

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS

Unloading of an elastic plastic loaded spherical contact

Finite element simulations of fretting contact systems

Engineering Solid Mechanics

Optimization of blank dimensions to reduce springback in the flexforming process

Process Modeling and Thermal/Mechanical Behavior of ACA/ACF Type Flip-Chip Packages

ScienceDirect. Bauschinger effect during unloading of cold-rolled copper alloy sheet and its influence on springback deformation after U-bending

A SOFTWARE SOLUTION FOR ADVANCED FRICTION MODELING APPLIED TO SHEET METAL FORMING

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases

The University of Melbourne Engineering Mechanics

UNIT I SIMPLE STRESSES AND STRAINS

Supplementary Material

Failure analysis of serial pinned joints in composite materials

THE EFFECTS OF LOCAL BUCKLING ON THE CRASH ENERGY ABSORPTION OF THIN-WALLED EXPANSION TUBES

TE 75R RESEARCH RUBBER FRICTION TEST MACHINE

Two Tier projects for students in ME 160 class

A fatigue limit diagram for plastic rail clips

Chapter 6: Mechanical Properties of Metals. Dr. Feras Fraige

Peak Strain and Displacement Sensors for Structural Health Monitoring

Practice Final Examination. Please initial the statement below to show that you have read it

Objectives: After completion of this module, you should be able to:

FEA MODELING OF A TRIBOMETER S PIN AND DISK INTERACTION

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.

SN QUESTION YEAR MARK 1. State and prove the relationship between shearing stress and rate of change of bending moment at a section in a loaded beam.

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Agricultural Science 1B Principles & Processes in Agriculture. Mike Wheatland

Modeling of the Bending Stiffness of a Bimaterial Beam by the Approximation of One-Dimensional of Laminated Theory

EE C247B / ME C218 INTRODUCTION TO MEMS DESIGN SPRING 2016 C. NGUYEN PROBLEM SET #4

Evaluation and Description of Friction between an Electro-Deposited Coating and a Ceramic Ball under Fretting Condition

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

A study of action point correction factor for L type flanges of wind turbine towers

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM

DEFORMATION PATTERN AND FAILURE CRITERIA OF WOVEN COMPOSITE PREFORM IN GENERAL BIAS EXTENSION

ME 243. Mechanics of Solids

Size Effects In the Crushing of Honeycomb Structures

Prediction of geometric dimensions for cold forgings using the finite element method

Structural Analysis Laboratory. Michael Storaker, Sam Davey and Rhys Witt. JEE 332 Structural Analysis. 4 June 2012.

Model tests and FE-modelling of dynamic soil-structure interaction

Effects of Ball Groupings on Ball Passage Vibrations of a Linear Guideway Type Ball Bearing Pitching and Yawing Ball Passage Vibrations

Chapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

GB/T / ISO 527-1:1993

9 MECHANICAL PROPERTIES OF SOLIDS

Effects of Forming Conditions of Roll Offset Method on Sectional Shape at the Corner of Square Steel Pipe +

Behavior of Tubular Hat Structure Under Three Point Bending

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 20, 2011 Professor A. Dolovich

Tensile stress strain curves for different materials. Shows in figure below

Numerical modeling of sliding contact

SHAKING TABLE TEST FOR FRICTIONAL ISOLATED BRIDGES AND TRIBOLOGICAL NUMERICAL MODEL OF FRICTIONAL ISOLATOR

Chapter 2 Surface Acoustic Wave Motor Modeling and Motion Control

Transcription:

Friction Properties of Surface with Circular Micro-patterns Hideo Koguchi Mechanical Engineering, 603- Kamitomioka, Nagaoka Univ. of Tech., Nagaoka, Niigata, Japan Email: koguchi@mech.nagaokaut.ac.jp Takayoshi Tanaka Graduate of Nagaoka Univ. of Tech., 603- Kamitomioka, Nagaoka Univ. of Tech., Nagaoka, Niigata, Japan Email: s064 @stn.nagaokaut.ac.jp predict pressure distribution in elasto plastic contacts. Lin et al. [5] presented a method for solving the twodimensional (2-D) isothermal rough surface contact problem of general anisotropic materials with friction. Furthermore, a lot of studies on contact analysis have been done until now [6]-[5]. It is known that the surfaces with regular patterns have a specified friction property and wettability. Recently, it is required that the flow of materials into a mold is managed by the suppression part of wrinkles in a press forming. Now, the flow of material during the press process may be controlled using the press die surface with a fine micro-pattern. Arbitrary shapes of micro fine patterns can be produced using a recent advanced process technology. Then, there are too many combinations of the shape, the sie and the arrangement of patterns. Furthermore, eperiments on the investigation of friction property for the patterned surface have been conducted. However, all combinations of the shape and the number of pattern cannot be investigated in eperiment. So, in the present study, a method that the friction property of the surface with a fine pattern can be estimated from a numerical simulation of three-dimensional contact analysis will be presented, and the evaluated values of friction coefficient are compared with those obtained from eperiment. Abstract From the previous study, it is found that the spring-back of plate in a press forming can be reduced by changing the friction property of the plate and die in a press process. In the present study, the friction property of the die surface is varied by machining a micro-pattern on the surface. However, it is very hard to conduct the eperiment for all combinations of different shapes of pattern and the number of the pattern per area. Hence, it is desired to estimate the friction property using a numerical simulation. In the present paper, a three-dimensional contact analysis is carried out for a circular micro-pattern and various contact conditions, an apparent friction coefficient of the plate with the pattern will be estimated through the comparison of eperiment and numerical simulation. It will be shown that the friction coefficient can be evaluated using the numerical simulation. Inde Terms Micro-pattern, contact simulation I. INTRODUCTION Until now, a lot of eperiment on friction control for improving the performance of machines and the quality of industrial products has been carried out. There is several ways for controlling the friction property of surface. One solution is to manufacture micro-patterns on a surface. Gong et al. [] conducted an elasto-plastic finite element analysis of a sphere in normal and sliding contact with a layered medium possessing a patterned surface with regularly spaced rectangular pads. They investigated the effect of pattern geometry on the contact pressure distribution and subsurface stress-strain field. Tianiang et al. [2] developed an element-free Galerkinfinite element (EFG-FE) coupling method, combined with the linear mathematical programming technique, to solve two-dimensional elasto-plastic contact problems. Dini et al. [3] investigated the contact pressure and shear traction distributions for a sphere pressed onto an elastically similar half-space whose surface is populated by a uniform array of spherical asperities. They showed that the rough contact absorbs less energy than a smooth one subject to the same loading history under light shear loads. Loan et al. [4] developed a simplified algorithm to II. A. Eperiment Setup Lever Grip in a tensile testing machine Dead load Jig fiing plate with a micro-pattern Pulley thin plate Figure. A schematic view of friction test apparatus. Eperimental apparatus is developed for evaluating the friction property of the plate with a pattern. A flat thin Manuscript received December 24, 206; revised June 0, 207. doi: 0.878/ijmerr.6.4.285-289 EXPERIMENTS 285

region of flat thin plate and the patterned plate is 0 20 = 200mm2, and the number of circular pattern in the contact area is 42. The pitch and sie of the circular pattern are shown in Fig. 4 and Table I. Material in the patterned plate is SKD. For the evaluation of the friction coefficient, the normal load applied to the flat thin plate and the tensile load are recorded using computer and AD converter during the test. Eperiment for several test conditions is prepared. Conditions in eperiment are shown in Table II. For instance, conditions: ) The flat thin plate are pressed by smooth flat plates, 2) The flat thin plate is pressed by the patterned plates, 3) The patterned plates processed by is used, 4) and 5) One side surface in the flat thin plate is pressed by the plate with the pattern and the other is pressed by the smooth flat plate. These conditions are referred to as conditions, 2, 3, 4 and 5, respectively. plate is set between two plates with the micro-pattern and is a constant load in the normal direction is applied to the lat thin plate. Fig. shows a schematic view of the friction test apparatus. In this apparatus, a constant load is applied to the plate with the pattern using the principle of leverage. A flat thin plate is pressed by the patterned plates and is pulled at a constant tensile speed, mm/min, using a tensile testing machine as shown in Fig. 2. TABLE I. PARAMETERS IN EXPERIMENTS AND SIMULATIONS thin plate Parameters Figure 2. Tensile testing machine (Tolerable load: 50kN) and friction test apparatus. Micro-patterns B. Eperimental Condition A precise figure for contact status is presented in Fig. 3. Plate with micro-pattern Normal load P0 Young s modulus, Poisson s ratio Normal load P0 Tensile speed Diameter d Pitch Pt Circle Height Angle Young s modulus, Poisson s ratio steel Aluminu m 20 GPa 68.3 GPa 0.33 0.34.6kN 520N.0 mm/min 00 mm 724 mm 30.0 mm 60.0º 20 GPa 0.33 TABLE II. CONDITIONS IN EXPERIMENT Condition 2 3 4 5 Plate with micro-pattern thin plate Upper surface Coating III. Lower surface Coating NUMERICAL SIMULATION Tensile load F Figure 3. Precise contact status of flat thin plate with patterned surface. A B pt d q D C Figure 5. Model for contact analysis using MARC. The present study aims to investigate how to evaluate the friction coefficient between the flat thin plate and the patterned plate using a numerical simulation. In the present simulation, a general-purpose finite element program called as MARC is used for a three-dimensional contact analysis. All materials used in the simulation are isotropic and elastic. Contact status in eperiment is Figure 4. Configuration of micro-pattern. In the present study, a circular pattern with sies as shown in Table I is processed using a chemical etching technique. After the process, hard chromium electroplating or physical vapor deposition () of chromium is carried out on the patterned surface. Contact 286

modelled as shown in Fig. 5. The total number of element is 370,004 and total number of node is 365,09. the flat plate, and then the plate with the micro-patterns slip on the flat thin plate. Then, several scratches of micro-pattern are observed. [mm] 0.724 2.6 0. 0.36.44 0.7 0.03 y 0.72 Figure 8. Photograph of surface with fine circular pattern. Figure 6. The sie of model for analysis. Displacement u.0mm (a) -y plane y P0 Figure 9. Photograph of scratched pattern in the steel flat plate. mma=0.368 600 Load F, N (b) - cross section The minimum sie of element is 7m. The sies of model for analysis are 2.6mm in length, 0.724mm in width and 0.7mm in thickness. Micro-pattern shape is circle in the analysis as shown in 6. The micro-pattern is 30 m in height. One circular pattern is placed at the center of the upper and lower surfaces and four one-half patterns are set on the both surfaces as shown in Fig. 5. In modeling of contact status, the symmetry of contact status with respect to the -ais is considered. Boundary condition for the analysis is shown in Fig. 7. Fig. 7(a) and Fig. 7(b) represent the boundary conditions on the -y plane and - cross section in the model. At the first stage, the flat thin plate is pressed by the patterned plate, at the second stage, the left side surface in the flat thin plate is pulled in the -direction at a constant speed. IV. mma=0.245 m min=0.28 400 =0.236 300 200 : surface : surface 00 0 0.0 0.5.0.5 2.0 2.5 Displacement of grip, mm (a) Steel (Normal load:.6kn) mma=0.368 600 =0.36 m min=0.353 500 Load F, N Figure 7. Boundary conditions in the simulation. =0.36 m min=0.353 500 mma=0.245 m min=0.28 400 =0.236 300 200 : surface : surface 00 0 0.0 0.5.0.5 2.0 2.5 Displacement of grip, mm RESULTS IN EXPERIMENT AND NUMERICAL ANALYSIS (b) Aluminum (Normal load: 520N) Figure 0. Load against the displacement of grip. Fig. 8 represents a photograph of surface in the patterned surface. Very fine circular patterns are observed. Fig. 9 demonstrates a photograph of scratches formed in the flat thin plate of steel after eperiment. At first, several fine traces of the circular micro-pattern remain on 287 Fig. 0 demonstrates a relationship between the tensile load and the displacement of grip in the tensile testing machine. Fig. 0(a) represents the result of steel flat plate and 0(b) that of aluminum flat plate. In the steel flat

displacement is integrated for the cross section of thin plate. This force balances with the total of shear stress. The friction coefficient, a, can be evaluated using Eq. (). () ma = ò s da P0 A where A represents the area of cross section of the thin plate and P0 is the normal load. The friction coefficient of surface with circular micropatterns in condition 2 for the steel flat plate was estimated as ma = 0.38. This value in simulation is a little bit larger than that in eperiment. In the present analysis, it is assumed that all materials are elastic and isotropic. In the net step, at least the flat thin plate should be set to be elasto-plastic materials. Though there are several issues in the present simulation, the calculated value is fairly well. So, it is supposed that the proposed method can be adequately estimated the friction coefficient. plate, the grip do not move apparently until the tensile load of 24N, after that, the load increases gradually. The steel plate without patterns slips about 400N under the condition and 600N for the plate with pattern under the condition 2. In the aluminum plate, the grip does not move until 22.5N of the tensile load. The tensile load increases with the displacement of grip, the aluminum plate slips around 80N for smooth flat surface of condition and 5N for the patterned plate of condition 2. Then, an apparent friction coefficient is evaluated from a ratio of the normal load to the tensile load at slip. It is found that the values of friction coefficient in the patterned surface are about.2-.5 times larger than those without the pattern. The values of friction coefficient for different conditions in eperiment are shown in Fig.. The results in the steel flat plate are shown in Fig. (a) and those in the aluminum flat plate are in Fig. (b). It is found that the value in condition 2 is the largest one than that in every other condition. The friction coefficients in condition 4 and 5, in which patterned surface is used in one side, yield the middle value of conditions and 2. (a) Steel flat plate Figure 2. Deformation of whole model. Upper surface P0 Micro-pattern δ Thin plate Lower surface (b) Aluminum flat plate Figure 3. Mapping plane of shear stress. Figure. A comparison of friction coefficients for several contact conditions. Then, the friction coefficient will be controlled by combining these contact conditions. Net, the results in the simulation are shown. In the simulation, the friction coefficient is set to be 0.24, which is obtained from eperiment. Fig. 2 shows the deformation of the thin plate in the model for analysis. The flat thin plate firstly deforms by the normal load. The distribution map of shear stress on the upper surface indicated in Fig. 3 is shown in Fig. 4. It is found that shear stress concentrates at the edge of each micropattern in the tensile direction. Net, an idea for evaluating a friction coefficient will be eplained. It is very hard to integrate the shear stress for the whole contact areas due to the compleity of geometry of contact area. So, tensile stress,, in the tensile direction at the side surface applying the forced 288 Figure 4. Contour map of shear stress, y V. CONCLUSIONS In the present paper, friction coefficients were evaluated by eperiment and the numerical simulation for improving the controllability of metal forming in the press working. In particular, a method calculating the friction coefficient using the result obtained from the

simulation was presented. Then, the stress distribution around the pattern was also investigated. The results in the present study are summaried as follows:. Newly developed friction testing machine yielded repeatedly the similar value of the friction coefficient for the same condition in eperiment. The friction coefficient in the aluminum plates was larger than that in the steel plates. The friction coefficient in the patterned surface plates was larger than that in the smooth surface. 2. The distribution of shear stress on the patterned plate concentrated at the edge of the circular micropattern. 3. Evaluated friction coefficients were almost similar values to the values in eperiment. This result indicates that the calculation method of the friction coefficient in the simulation is adequate. [8] [9] [0] [] [2] [3] [4] REFERENCES [] [2] [3] [4] [5] [6] [7] Z. Q. Gong and K Komvopoulos, Mechanical and thermomechanical elastic-plastic contact analysis of layered media with patterned surfaces, Journal of Tribology, vol. 26, pp. 9-7, January 2004. L. Tianiang, L. Geng, and Q. J. Wang, An element-free Galerkin-finite coupling method for elasto-plastic contact problems, Journal of Tribology, vol. 28, pp. -9, January 2006. D. Dini and A. D. Hills, Frictional energy dissipation in a rough hertian contact, Journal of Tribology, vol. 3, pp. 0240--8, April 2009. U. Loan and S. Sergiu, A simplified model for pressure distribution in elastic perfectly plastic contacts, Fascicle of Management and Technological Engineering, vol. 9, pp. 07-4, 2002. Y. Lin and T. C. Ovaert, A rough surface contact model for general anisotropic materials, Journal of Tribology, vol. 26, pp. 4-49, January 2004. A. Beheshti and M. M. Khonsari, Asperity micro-contact models as applied to the deformation of rough line contact, Tribology International, vol. 52, pp. 6-74, February 202. S. Kucharski and G. Starynski, Study of contact of rough surfaces: Modeling and eperiment, Wear, vol. 3, pp. 67-79, January 204. [5] K. Poulios and P. Klit, Implementation and applications of a finite-element model for the contact between rough surfaces, Wear, vol. 303, pp. -8, February 203. W. Wayne and Q. J. Wang, A numerical static friction model for spherical contacts of rough surfaces, influence of load, material, and roughness, Journal of Tribology, vol. 3, pp. 8, April 2009. J. M. Garcia and A. Martini, Measured and predicted static friction for real rough surfaces in point contact, Journal of Tribology, vol. 34, pp. 0350-8, July 202. H. Sojoudi and M. M. Khonsari, On the behavior of friction in lubricated point contact with provision for surface roughness, Journal of Tribology, vol. 32, pp. 0202--8, January 200. R. M. N. Fleury, D. A. Hills, and J. R. Barber, A corrective solution for finding the effects of edge-rounding on complete between elastically similar bodies. Part I: Contact law and normal contact considerations, International Journal of Solids and Structures, vol. 85-86, pp. 89-96, February 206. D. M. Mulvihill, M. E. Kartal, D. Nowell, and D. A. Hills, An elastic-plastic asperity interaction model for sliding friction, Tribology International, vol. 44, pp. 679-694, June 20. L. Xiao, S. Bjorklund, and B. G. Rosen, The influence of surface roughness and the contact pressure distribution on friction in rolling/sliding contacts, Tribology International, vol. 40, pp. 694698, November 2005. M. Satoru, I. Fumihiro, and N. Takashi, Effect of normal load on friction coefficient for sliding contact between rough rubber surface and rigid smooth plane, Tribology International, vol. 92, pp. 335-343, July 205. Hideo Koguchi has completed Postgraduate Doctor Program in Mechanical Engineering at Tokyo University In March 98, his research is on strength of materials, such as engineering polymers. He participated in various international research collaborations with Washington University (989) and Pennsylvania University (989). His current job is as Full Professor in the field of Material Engineering at Department of Mechanical Engineering, Engineering Faculty, Nagaoka University of Technology, Nagaoka, Niigata, Japan. Prof. Koguchi is a member of American Society of Mechanical Engineers (ASME), The Society of Materials Science of Japan(JSMS) and The Japan Society of Mechanical Engineers(JSME). Takayoshi Tanaka is Graduate of Nagaoka University of Technology. A member of The Japan Society of Mechanical Engineers (JSME). 289