FEM Analysis of Punching-Process in Consideration of Micro Die Wear

Similar documents
Coming Out Mechanism of Steel Shaft from Ceramic Sleeve

Analysis of axisymmetric cup forming of metal foil and micro hydroforming process

Analysis of forming- Slipline Field Method

Modelling of ductile failure in metal forming

Available online at ScienceDirect. Procedia Engineering 81 (2014 )

End forming of thin-walled tubes

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

Optimization of blank dimensions to reduce springback in the flexforming process

Generation Mechanism of Driving Out Force of the Shaft from the Shrink Fitted Ceramic Roll by Introducing Newly Designed Stopper

The plastic behaviour of silicon subjected to micro-indentation

Inverse identification of plastic material behavior using. multi-scale virtual experiments

Finite Element Simulation with Coupled Thermo-Mechanical Analysis of Superplastic Dieless Tube Drawing Considering Strain Rate Sensitivity*

Study on the Energy Absorption and the Local Buckling of the Expansion Tubes

Thermal Stress Analysis of BOF Accounting for Friction Force Experimentally Measured

FRACTURE BEHAVIOR OF RIVETED LAP JOINTS DUE TO PROJECTILE IMPACTS

Significance of the local sheet curvature in the prediction of sheet metal forming limits by necking instabilities and cracks

Abstract. 1 Introduction

Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture

Simulation of the effect of DIE Radius on Deep Drawing Process

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.

Advanced Structural Analysis EGF Cylinders Under Pressure

Failure from static loading

Neutral Radius Value Determination by Numerical Simulation Method at Ring Upsetting Test

Indentation cracking of glass II: Relationship between crack resistance and equibiaxial flexural strength

Arch. Metall. Mater. 62 (2017), 1, 85-90

An Analytical Model for Long Tube Hydroforming in a Square Cross-Section Die Considering Anisotropic Effects of the Material

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

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

INFLUENCE OF PROCESSING PARAMETERS ON FORMING QUALITY OF NON-CIRCULAR SPINNING

Influence of impact velocity on transition time for V-notched Charpy specimen*

Evaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis

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

Module 5: Theories of Failure

INTERFACIAL STRENGTH EVALUATION IN A GLASS FIBER REINFORCED COMPOSITE USING CRUCIFORM SPECIMEN METHOD

FCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering

Advanced Friction Modeling in Sheet Metal Forming

Analysis of shear cutting of dual phase steel by application of an advanced damage model

ULTRASONIC CHARACTERIZATION OF RESIDUAL STRESS AND TEXTURE IN CAST STEEL RAILROAD WHEELS

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

MMJ1133 FATIGUE AND FRACTURE MECHANICS A - INTRODUCTION INTRODUCTION

University of Sheffield The development of finite elements for 3D structural analysis in fire

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

INVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING ANISOTROPY

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

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

A slip-line solution to metal machining using a cutting tool with a step-type chip-breaker

The Nonlinear Dynamic Behavior Analysis and Research for Sheet Steel Stamping Forming Problems Based on Finite Element Method

INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM

MAE 322 Machine Design Lecture 2. Dr. Hodge Jenkins Mercer University

Finite element simulations of fretting contact systems

VALIDATING DYNAMIC TENSILE MECHANICAL PROPERTIES OF SHEET STEELS FOR AUTOMOTIVE CRASH APPLICATIONS

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

Wear 265 (2008) Contents lists available at ScienceDirect. Wear. journal homepage:

Lecture #8: Ductile Fracture (Theory & Experiments)

Discrete Element Modelling of a Reinforced Concrete Structure

Engineering Science OUTCOME 1 - TUTORIAL 4 COLUMNS

Influence of Back Pressure on Slab Edge Deformation Behavior during Width Reduction Pressing *1

IMPACT STRENGTH AND RESPONSE BEHAVIOR OF CFRP GUARDER BELT FOR SIDE COLLISION OF AUTOMOBILES

A NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES

Influence of the collision speed and angle of a bullet: experimental reconstruction of bullet configuration and FE-analysis

5. STRESS CONCENTRATIONS. and strains in shafts apply only to solid and hollow circular shafts while they are in the

Special edition paper

#695 Numerical and experimental study on determination of the contact temperature in an industrial ironing operation

Mesoscopic Simulation of Failure of Mortar and Concrete by 3D RBSM

Effects of the semi die/plug angles on cold tube drawing with a fixed plug by FEM for AISI 1010 steel tube

Stress Analysis of Mandrel Used In Pilger Technology

LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL

Boundary and Mixed Lubrication Friction Modeling under Forming Process Conditions

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

Effect of intermediate principal stresses on compressive strength of Phra Wihan sandstone

Modelling the behaviour of plastics for design under impact

ME 2570 MECHANICS OF MATERIALS

Three-Dimensional Finite Element Analysis of Material Nonlinearity

Design and Analysis of Progressive Tool

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

Plane Strain Test for Metal Sheet Characterization

Compression Molding of Green Composite made of Wood Shavings

Effect of Mandrel, Its Clearance and Pressure Die on Tube Bending Process via Rotary Draw Bending Method

When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress?

Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

Investigation of the deformation stability in the incremental sheet forming process

MAE 322 Machine Design. Dr. Hodge Jenkins Mercer University

On Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models

Analysis of a Lap Joint Including Fastener Hole Residual Stress Effects

Predicting the shear strength of round clinched joint

Effect of Deformation Mode of Cylindrical Tubes with Corrugated Surface Dimensional Errors

Tool edge radius effect on cutting temperature in micro-end-milling process

Friction Properties of Surface with Circular Micro-patterns

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

Stress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy

Influence of fibre proportion and position on the machinability of GFRP composites- An FEA model

Simulation of the cutting action of a single PDC cutter using DEM

THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT

NUMERICAL AND EXPERIMENTAL STUDY OF FAILURE IN STEEL BEAMS UNDER IMPACT CONDITIONS

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

Stress Intensity Factors Analyses of a Three-Dimensional Interface Crack Between Anisotropic Dissimilar Materials Under Thermal Stress

Bulk Metal Forming II

A study of forming pressure in the tube-hydroforming process

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

Transcription:

FEM Analysis of Punching-Process in Consideration of Micro Die Wear Takashi Ueda 1,a, Takashi Iizuka 1 and Shinichi Enoki 2 1 Department of Mechanical and System Engineering, Kyoto Institute of Technology, Matsugasaki Goshokaido-cho, Sakyo-ku, 606-8585 Kyoto, Japan 2 National Institute of Technology, Nara college, Yata-cho, Yamatokooriyama-shi, 639-1080 Nara, Japan Abstract. Detailed punching-process analysis is difficult because the punch diameter, sheet thickness, clearance, and die wear of the punching-process are approximately 10 mm, 1 mm, 100 m, and 10 m, respectively. Furthermore, such analyses have rarely considered die wear. Therefore, in the present study, we set the mesh size in the neighborhood of the die corner radius and the punch to 1/10 m. Moreover, we analyze the punching process considering die wear. We investigated the stress, incremental strain, and equivalent strain for different levels of die wear. 1 Introduction Burr formation in the punching process has adverse effects on post-processing, product accuracy, and worker safety. Thus, a burr-free punching method is desired. Therefore, burr-free blanking using servo press machine on a few steps punching method was reported[1, 2]. However, preventing burr formation during the one step punching process is difficult due to the nature of the formation mechanism of the punching process. Reducing burrs by a simpler method may be possible if the mechanism of burr formation can be clarified. A number of experimental and analytical investigations dealing with burr shape formation have been conducted [3]. The results of these studies revealed that the clearance [4] and die wear [5] primarily affected burr shape formation. However, the formation mechanism of burr shape in the punching process has not yet been clarified. Therefore, Iizuka [6, 7] and Sakamoto [8] focused on die wear that affects the formation of burrs and investigated the burr formation mechanism in the punching process. The burr formation mechanism in the punching process has been investigated experimentally while focusing on die wear. However, in order to support the experimental results, a detailed numerical analysis of the burr formation mechanism in the punching process is needed. Finite element method (FEM) analysis has been reported as being effective in plastic formation. A number of FEM analyses have been conducted in order to clarify the formation mechanism [9, 10] and crack growth [11, 12] of the punching process. However, detailed analysis of the punching process is more difficult than the analysis of other formation mechanisms. Since orders of magnitude on the punch diameter, the sheet thickness, the clearance, and the die wear of the punching-process are a Corresponding author: m5623104@edu.kit.ac.jp approximately 10 mm, 1 mm, 100 m, and 10 m, respectively, they are very difference order of magnitude. So, there are reports that the mesh size in the neighborhood of the die corner radius and punch was set 10 m [4, 10]. But, Conventional analyses of the punching process rarely consider micro die wear. Therefore, we set the mesh size in the neighborhood of the die corner radius and punch to 1/10 m. Moreover, we analyze the punching process while considering die wear. We varied the die wear amount as 10 m, 20 m, 30 m, and 40 m. We investigated the hydrostatic pressure distribution, the equivalent plastic strain distribution, and the plastic incremental strain distribution for different levels of die wear in the neighborhood of the die corner radius. 2 Analysis method The punching process analysis of the present study was conducted using Simufact.forming12.0.3 (Simufact Engineering System). In the present study, an axisymmetric two-dimensional model was used as an analysis model, as shown in Figure 1. At this time,r d, R p and R h are die corner radius, punch corner radius and die holder shoulder round. We simulated the case in which a blank was punched by a metal punch with a diameter of 4 mm. Cold-rolled steel sheet (SPCC) with a thickness of 1.6 mm was used as blank. The punch and die holder are rigid bodies. The coefficient of friction between the punch and die and the blank was 0.2. The blank was subjected to a forced displacement of 1.4 mm with a processing velocity of 500 mm/s. The base mesh size of the blank in this analysis was set to 37.5 m. The mesh size in the neighborhood of the die corner radius and the

MATEC Web of Conferences punch was set to 4.7 m. The clearance was 4.69%. In the analysis, the die corner radius R d was varied as 10 m, 20 m, 30 m, and 40 m, and the punch corner radius R p and die holder shoulder round R h were set to 20 m. and the die side in the case of R d = 20 m. However, the plastic incremental strain is distributed on the die corner radius only in case of R d = 40 m. Next, we confirm the plastic incremental strain distribution for a stroke of 1.2 mm. In this case, a gap between the die Figure 1. Analysis model. 3 Analysis result 3.1 Plastic incremental strain For a stroke of 1 mm, the plastic incremental strain distribution is shown in Figure 2 for die corner radius radii of 10 m, 20 m, 30 m, and 40 m. Figure 3. Plastic incremental strain distribution on neighborhood of die edge with every stroke. Figure 2. Plastic incremental strain distribution with a stroke of 1mm. As shown in Figure 2, thin layer shear is thought to occur in the domain linking the punch corner radius and the die corner radius because high plastic incremental strain is confirmed in this region. Next, we confirm plastic incremental strain in the neighborhood of the die corner radius for stroke of 1 mm, 1.2 mm, or 1.28 mm. The plastic incremental strain distributions are shown for die corner radius radii of 20 m and 40 m in Figure 3. As shown in Figure 3, we confirm that the plastic incremental strain that occurs in the neighborhood of the die corner radius is distributed on the die corner radius and the blank occurs at the die corner radius in the cases of R d = 20 m and 40 m. The gap is bigger in the case of R d = 40 m than in case of R d = 20 m. For a stroke of 1.28 mm, we can confirm that the gap is distributed on the die corner radius and the die side in the case of R d = 20 m. However, for the case in which R d = 40 m, the gap is larger than that in the case of a stroke of 1.2 mm. We cannot confirm the existence of a gap on the die side. Next, the plastic incremental strain distributions for a stroke of 1.2 mm in the cases of R d = 10 m, 20 m, 30 m, and 40 m are shown in Figure 4. As shown in Figures 4(a) and 4(b), in the cases of R d = 10 m and 20 m, the gap is distributed on the die corner radius and the die side. In case of R d = 10 m, the gap is focused on the die side. However, in the case of R d = 20 m, the gap is focused on the die corner radius. As shown in Figures 4(c) and 4(d), in the cases of R d = 30 m and 40 m, the gap is distributed on only on the die corner radius. We cannot confirm the existence of a gap

NUMIFORM 2016 on the die side. Therefore, when the die corner radius R d increases, the domain that the gap is focused moves from the die corner radius to the die side. The gap is thought to satisfy the fracture condition, because the contact pressure is not focused in the gap. As shown in Figure 3, fracture is thought to occur from the end part of the gap for all die conditions, because, for a stroke of 1.28 mm, high plastic incremental strain occurs at the end of the gap. The hydrostatic pressure distribution, the maximum principal stress distribution, and the minimum principal stress distribution with a stroke of 1 mm in the case of R d = 20 m are shown in Figure 5, Figure 6 and Figure 7. Figure 7. Minimum principal stress distribution with a stroke of 1mm in case of R d = 20μm Figure 4. Plastic incremental strain distribution with a stroke of 1.28mm on ever die abrasion status. 3.2 Hydrostatic pressure Figure 5. Hydrostatic pressure distribution with a stroke of 1mm in case of R d = 20μm Figure 8. Hydrostatic pressure distribution with every stroke in case of R d = 20μm Figure 6. Maximum principal stress distribution with a stroke of 1mm in case of R d = 20μm In the present study, the punching process exhibits a circumferential plane strain field because, as shown in Figure 5, Figure 6 and Figure 7, the mean of the maximum principal stress distribution and the minimum principal stress distribution is the same as the hydrostatic pressure distribution. Therefore, the mechanism of the punching process can be explained by the hydrostatic

MATEC Web of Conferences pressure. As shown in Figure 5, tensile stress occurs in the neighborhood of the die corner radius and the punch corner radius. The tensile stress in the neighborhood of the die edge is shown as a tensile stress field in Figure 6 and Figure 7. Compression stress occurred at the in bottom face of the punch, at the top face of the die, at the die corner radius is shown as a tensile stress field the punch side, and at the die side. Therefore, the blank on the die top face is expected to flow toward the punch side. As for the blank on punch top face is toward die side. Next, the hydrostatic pressure distributions for strokes of 1 mm, 1.2 mm, and 1.28 mm in the case of R d = 20 m are shown in Figure 8. than elsewhere. In Step 2, the tensile stress appears in the pure shear domain in Step 1. In Step 3, the tensile stress in this domain is approximately the same as that in the neighborhood of the die corner radius. Next, the hydrostatic pressure distributions in the neighborhood of the die corner radius for R d = 20 m and 40 m are shown in Figure 9. As shown in Figure 9, in Step 1, the tensile stress is distributed only on the end part of the die corner radius. Compression stress occurs at locations other than the die corner radius. The tensile stress in Step 2 is higher than that in Step 1. Moreover, a gap appears between the die corner radius and the blank. The high tensile stress domain in Step 3 is larger than that in Step 2. The tensile stress domain for R d = 40 m is larger than that for R d = 20 m. Therefore, it's expected that gap form varies depending on the level of die corner radius. 3.3 Equivalent plastic strain The equivalent plastic strain distributions for a stroke of 1 mm in the cases of R d = 10 m, 20 m, 30 m, and 40 m are shown in Figure 10. As shown in Figure 10, high equivalent plastic strain is distributed over the range of the die side from the die corner radius. Therefore, the blank over the range from the die corner radius to the die side is thought to be brittle. We confirmed that die wear has little effect on equivalent plastic strain. Figure 9. Hydrostatic pressure distribution on neighborhood of die edge. As shown in Figure 8, we confirmed that the stress distribution for the punching process is divided into three steps. In Step 1, the hydrostatic pressure in the domain linking the punch corner radius and the die corner radius is approximately 0 MPa. This domain is thought to be a pure shear domain. Next, the hydrostatic pressure distributions in the neighborhood of the die corner radius for strokes of 1 mm, 1.2 mm, and 1.28 mm in the case of R d = 20 m and 40 m are shown in Figure 9. The tensile stress in the neighborhood of the domain is shown. The tensile stress in the neighborhood of the die corner radius and the punch corner radius is higher Figure 10. Equivalent plastic strain distribution with a stroke of 1mm. 4 Discussion The plastic incremental strain, hydrostatic pressure, and equivalent plastic strain are high for a stroke of 1 mm, as shown in Figure 11. As shown in Figure 11, the high plastic incremental strain occurred at the die corner radius. And, the high equivalent plastic strain is distributed over the range of the die side from the die corner radius. However, high hydrostatic pressure

NUMIFORM 2016 occurred at the neighborhood of the die corner radius. The domain that high equivalent plastic strain was shown is thought to be brittle. The plastic incremental strain shows the deformation volume. Therefore, the location at which the equivalent plastic strain overlaps the plastic incremental strain is likely to exhibit crack formation. Cracks are thought to progresses as a result of the hydrostatic pressure of the tensile field. Therefore, in case of Rd = 20 μm, crack is thought to occurred over the range of the die side from the die corner radius. In case of Rd = 40 μm, crack is thought to occurred on die corner radius. The origin of crack formation differs according to the level of die wear. (1) In the plastic incremental strain distribution, the thin layer shear is thought to occur in the domain linking the punch corner radius and the die corner radius. (2) The gap between the die and the blank occurs at the die corner radius for a stroke of 1.2 mm. Moreover, the gap differs depending of the level of die corner radius. (3) In this punching process, the stress distribution of the thin layer shear is divided into three steps. (4) We confirmed that the level of die wear has little effect on the equivalent plastic strain. (5) The origin of crack formation differs according to the level of die wear. Reference Figure 11. High value of plastic incremental strain, hydrostatic pressure and equivalent plastic strain with a stroke of 1mm. 5 Conclusion We analyze the punching process while considering die wear. We varied the die wear amount as 10 m, 20 m, 30 m, and 40 m. We investigated the hydrostatic pressure distribution, the equivalent plastic strain distribution, and the plastic incremental strain distribution for different levels of die wear in the neighborhood of the die corner radius. The results obtained are summarized as follows. 1. K. Junlapen, M. Kubo, K. Mizuochi, N. Koga, P. Keawtatip, ICTP, 158-163, (2008) 2. S. Thipprakmas, M. Jin, M. Murakawa, ICTP, Oct. 28-31, 889-894, (2002) 3. T. Maeda, I. Aoki, J. Fac. Eng. Univ. Tokyo (B), 34-3,537-577 (1978) 4. N. Hatanaka, K. Yamaguchi, N. Takakura, T. Iizuka, J. Mat. Pro. Tec., 140, 628-634(2003) 5. H. Jaafar, K. Mori, Y. Abe, K. Nakanishi, Journal of Materials Processing Technology, 227, 190-199, (2016) 6. T. Iizuka, N. Hatanaka, T. Ishihara, N. Takakura, ICTP, 4-1, 2330-2335, (2008) 7. T. Iizuka, Wei. Yu, N. Hatanaka, N. Takakura, J. Jpn. Soc. Technol., (2nd report), 337, (2009) 8. D. Sakamoto, T. Iizuka, N. Hatanaka, N. Takakura, Metal Forming, 537, (2012) 9. S. H. Jeong, J. J. Kang, S. I. Oh, J. ICTP, 631-634, (1996) 10. N. Hatanaka, K. Yamaguchi, N. Takakura, Journal of Materials Technology, 139, 64-70, (2003) 11. Y. Yoshida, S. Sumikawa, N. Yukawa, T. Ishikawa, Ad. Tec. Pla., 2 883-888, (2002) 12. A. Ahmed, H. Muhamad, IIUM, Vol. 13, 153-159, (2012)