CHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS

Similar documents
Example 24 Spring-back

COMPUTER AIDED ENGINEERING. Part-1

MAE 323: Lab 7. Instructions. Pressure Vessel Alex Grishin MAE 323 Lab Instructions 1

Using MSC.Nastran for Explicit FEM Simulations

CHAPTER 5 USE OF STL FILE FOR FINITE ELEMENT ANALYSIS

Lecture 3 : General Preprocessing. Introduction to ANSYS Mechanical Release ANSYS, Inc. February 27, 2015

Using three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model

CHAPTER-10 DYNAMIC SIMULATION USING LS-DYNA

CHAPTER 1. Introduction

Coupled Analysis of FSI

CHAPTER 8 FINITE ELEMENT ANALYSIS

An Optimization Procedure for. Springback Compensation using LS-OPT

ANSYS Element. elearning. Peter Barrett October CAE Associates Inc. and ANSYS Inc. All rights reserved.

A Coupled 3D/2D Axisymmetric Method for Simulating Magnetic Metal Forming Processes in LS-DYNA

Module: 2 Finite Element Formulation Techniques Lecture 3: Finite Element Method: Displacement Approach

TryItNow! Step by Step Walkthrough: Spoiler Support

Element Order: Element order refers to the interpolation of an element s nodal results to the interior of the element. This determines how results can

Finite Element Analysis of Ellipsoidal Head Pressure Vessel

TWO-DIMENSIONAL PROBLEM OF THE THEORY OF ELASTICITY. INVESTIGATION OF STRESS CONCENTRATION FACTORS.

Module 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis

Module 3 Mesh Generation

Using the Eulerian Multiphase Model for Granular Flow

CHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING

midas NFX 2017R1 Release Note

PARALLEL ENGINEERING SIMULATIONS BASED ON FORMING SIMULATION OF A HEAT EXCHANGER PLATE

ANSYS Workbench Guide

Exercise 1. 3-Point Bending Using the GUI and the Bottom-up-Method

Design Optimization of a Weather Radar Antenna using Finite Element Analysis (FEA) and Computational Fluid Dynamics (CFD)

Types of Idealizations. Idealizations. Cylindrical Shaped Part. Cyclic Symmetry. 3D Shell Model. Axisymmetric

Offshore Platform Fluid Structure Interaction (FSI) Simulation

Elastic Analysis of a Bending Plate

Computer Life (CPL) ISSN: Fluid-structure Coupling Simulation Analysis of Wavy Lip Seals

NX Advanced FEM. fact sheet

ANSYS/LS-Dyna. Workbench Using. Steven Hale Senior Engineering Manager CAE Associates, Inc. June 13, CAE Associates

Study of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket

Some Aspects for the Simulation of a Non-Linear Problem with Plasticity and Contact

PYTHON API TO RUN CALCULIX: PYCALCULIX Justin Black Mechanical Engineer

Step Change in Design: Exploring Sixty Stent Design Variations Overnight

ANSYS 5.6 Tutorials Lecture # 2 - Static Structural Analysis

Adaptive Finite Element Simulation of Sheet Forming Operations Using Continuum Elements

2008 International ANSYS Conference

5. Finite Element Analysis of Bellows

Introduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

Virtual Die Tryout of Miniature Stamping Parts

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS

Global to Local Model Interface for Deepwater Top Tension Risers

ECE421: Electronics for Instrumentation

Recent Advances on Higher Order 27-node Hexahedral Element in LS-DYNA

CHAPTER 7 BEHAVIOUR OF THE COMBINED EFFECT OF ROOFING ELEMENTS

NX Advanced FEM. Benefits

Through Process Modelling of Self-Piercing Riveting

Rigid Dynamic Analysis in Workbench

Wall thickness= Inlet: Prescribed mass flux. All lengths in meters kg/m, E Pa, 0.3,

SOLIDWORKS Simulation

Simulation of Laminar Pipe Flows

TOPOLOGY OPTIMIZATION OF WING RIBS IN CESSNA CITATION

Simulation of Overhead Crane Wire Ropes Utilizing LS-DYNA

GEOSTUDIO Tutorials Results and Procedure Comparison

Stress Analysis of thick wall bellows using Finite Element Method

Embedded Reinforcements

Validation Report: Additional Data Mapping to Structural Analysis Packages

Example Plate with a hole

MAE 455 COMPUTER-AIDED DESIGN AND DRAFTING MIDTERM EXAM PRACTICE QUESTIONS. Name: You are allowed one sheet of notes.

Predicting the mechanical behaviour of large composite rocket motor cases

An Explanation on Computation of Fracture Mechanics Parameters in ANSYS

A Graphical User Interface for Simulating Resin-Transfer-Molding Combining LS-DYNA and OpenFOAM

Essay 5 Tutorial for a Three-Dimensional Heat Conduction Problem Using ANSYS

Ansys Mechanical APDL

Orbital forming of SKF's hub bearing units

The Effect of Element Formulation on the Prediction of Boost Effects in Numerical Tube Bending

midas FEA V320 Table of Contents

Computational Simulation Of AT595 Blast Container using LS-Dyna

Program: Advanced Certificate Program

Engineering Analysis

This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:

Newly Developed Capabilities of DYNAFORM Version 5.0

Advances in LS-DYNA for Metal Forming (I)

KEYWORDS Morphing, CAE workflow, Optimization, Automation, DOE, Regression, CFD, FEM, Python

Facundo DEL PIN / Iñaki ÇALDICHOURY / Rodrigo PAZ / / Livermore Software Technology Corporation

Creo Simulate 3.0 Tutorial

THREE-DIMENSIONAL PROBLEM OF THE THEORY OF ELASTICITY STRESS IN A THICK-WALLED PRESSURE VESSEL

Finite Element Analysis using ANSYS Mechanical APDL & ANSYS Workbench

Elfini Solver Verification

Model Set up, Analysis and Results of the Inverse Forming Tool in ANSA

Exercise 2: Mesh Resolution, Element Shapes, Basis Functions & Convergence Analyses

VALIDATE SIMULATION TECHNIQUES OF A MOBILE EXPLOSIVE CONTAINMENT VESSEL

Finite Element Buckling Analysis Of Stiffened Plates

midas NFX An insight into midas NFX

A comparison of DYNA3D, NIKE3D and LS-DYNA

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

Benchmark on the numerical simulation of a tube bundle vibration under cross flow

CFD modelling of thickened tailings Final project report

Aufgabe 1: Dreipunktbiegung mit ANSYS Workbench

Evaluation of LS-DYNA Material Models for the Analysis of Sidewall Curl in Advanced High Strength Steels

Similar Pulley Wheel Description J.E. Akin, Rice University

Tut #2 LSDyna: Impact Analysis By C. Daley

Fluid Structure Interaction - Moving Wall in Still Water


LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.

Release 10. Kent L. Lawrence. Mechanical and Aerospace Engineering University of Texas at Arlington SDC PUBLICATIONS

Transcription:

54 CHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS 4.1 INTRODUCTION In Fluid assisted deep drawing process the punch moves in the fluid chamber, the pressure is generated in the fluid. This fluid pressure depends on the process geometry and various process parameters. Estimation of fluid pressure by using flotran CFD analysis software. ANSYS - LS Dyna finite element simulation software is used to predict the radial, hoop and drawing stresses with varying fluid pressures. 4.2 DETERMINATION OF FLUID PRESSURE Ansys - Flotran CFD analysis is used to study the variation of pressure of fluid with different punch radius at constant speed of punch and also pressure variation with different punch speeds at constant punch radius using three fluids such as castor oil, olive oil and heavy machine oil. This pressure of fluid is to evaluate the blank holding pressure analization of stresses in this process. The element type is fluid 141 element from flotran CFD library is selected for meshing.the FLUID 141 element shown in fig.4.1 This figure shows FLUID 141 geometry, locations of node and coordinate system for this element. The element is defined by three nodes [triangle] or four nodes [quadrilateral] and by isotropic properties of material.

55 Fig.4.1 FLUID 141 Element geometry The fluid model is developed in Ansys preprocessing using geometric modeling approach. The resulted geometry with 2D geometric options are shown in fig.4.2 Fig.4.2 Geometric model ( axisymmetric) of the deep drawing process Using adaptive mesh, a converged mesh is shown in fig.4.3

56 Fig.4.3 Flotran CFD model of the process The total number of elements and nodes in the model are 7972 and 8364. Boundary and loading conditions: Vx = Vy = 0 on the boundary and punch velocity, Vy = 10mm/sec. The fig.4.4 shows the boundary and loading conditions of process.

57 Fig.4.4 Boundary and loading conditions of the process Table 4.1 Variation of fluid pressure for different fluids by varying punch radius Punch radius Fluid Pressure P [ N/m 2 ] r p [mm] Heavy machine Olive oil oil Castor oil 10 2.23 11.25 31.5 15 3.01 16.52 45.2 20 3.52 20.7 59.2 25 4.0 29.0 65.0 30 5.8 31.23 85.5 35 9.76 35.0 101.2 40 14.134 58.47 121.24 45 18.37 63.2 127.5 50 25.43 71.5 138.7

58 Table 4.1 represents the variation of fluid pressure for three different fluids by increasing punch radius at constant punch speed u = 10mm/sec. From this methodology the fluid pressure is increases with increase in punch radius. The fluid pressure is also depends on the viscosity. When the viscosity of fluid is more, then the fluid pressure is high. The variation of fluid pressures is observed from the table at any punch radius is P olive oil < P Heavy machine oil < P castor oil The same procedure is used for study on variation of pressure of fluid with different punch speeds at constant punch radius using three fluids such as castor oil, olive oil and heavy machine oil with consideration of the process radius of punch 40mm and radius of die opening 45mm. Table 4.2 Variation of fluid pressure for different fluids by varying punch Speed. Punch speed Fluid Pressure P [ N/m 2 ] u [mm/sec] Heavy machine Olive oil oil Castor oil 5 10.42 27.83 73.25 10 14.134 58.47 121.24 15 16.32 70.4 164.5 20 17.98 89.72 215.02 25 18.73 107.2 263.7 30 21.53 122.55 310.4 Table 4.2 represents the variation of fluid pressure for three different fluids by increasing punch speed. From this methodology the fluid pressure is increases with increase in punch speed. The fluid pressure is also depends on the viscosity. When the viscosity of fluid

59 is more, then the fluid pressure is high. The variation of fluid pressures is observed from the table at any punch speed is P olive oil < P Heavy machine oil < P castor oil. 4.3 FINITE ELEMENT ANALYSIS OF DEEP DRAWING PROCESS ANSYS - LS Dyna finite element simulation software is used to predict the radial, hoop and drawing stresses with varying blank radius, blank thickness, punch radius, fluid viscosities and fluid pressures. A 2d solid 162 ( Ansys-Ls dyna library) element is used to model the billet and punch and die is modeled with rigid elements. Coupled degrees of freedom is imposed between punch and billet and billet and die. 4.3.1 PLANE 162 Element Description PLANE162 is used for modeling 2-D solid structures in ANSYS LS- DYNA. The PLANE 162 element shown in fig.4.5.the element can be used either as a planer or as an axisymmetric element. The element is defined by four nodes having six degrees of freedom at each node: translations, velocities, and accelerations in the nodal x and y directions. A three-node triangle option is also available, but not recommended. Fig. 4.5 PLANE 162 Element geometry

60 The element is used in explicit dynamic analyses only. When using this element, the model must only contain PLANE162 elements - we cannot mix 2-D and 3-D explicit elements in the same model. Furthermore, all PLANE162 elements in the model must be the same type (plane stress, plane strain, or axisymmetric). The deep drawing process is simulated with condition axisymmetric. The Geometric model of deep drawing process along with the boundary conditions shown in fig.4.6 (a) Fig.4.6 (a) Geometric model of deep drawing process with boundary Conditions The boundary conditions with reasons are adopted for this process as follows. AB : Ux = 0 [ symmetric axis ]

61 BC, CD and DE : Ux = Uy = 0 [ fixed, i.e. Die is supported by ground ] EF : Uy = 0 [ blank holder support ] FG and EH : Ux = 0 [ normal to the blank surface ] The punch moves in y direction, hence punch velocity Vy=10mm/sec. Finite element model of the deep drawing process shown in fig.4.6 [b] Fig.4.6 [b] Finite element model of the deep drawing process The total number of elements and nodes in the model are 1160 and 1200 respectively. LS - Dyna explicit solver is used to predict the radial, hoop and drawing stresses by varying the radius and thickness of blank for different punch radius using three different fluids. The variation of these stresses are reported in chapters 5, 6 and 7.