Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Similar documents
Prediction of Chatter Stability for Milling Process Using an Improved Semi-discretization Method

Delay equations with engineering applications Gábor Stépán Department of Applied Mechanics Budapest University of Technology and Economics

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

Lecture 8 Modal Analysis

2.29 Numerical Fluid Mechanics

Adaptive sliding mode reliable excitation control design for power systems

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

Portfolios with Trading Constraints and Payout Restrictions

Analysis of Dynamic Cross Response between Spindles in a Dual Spindle Type Multi-Functional Turning Machine

Dynamic Analysis Based On ANSYS of Turning and Grinding Compound Machine Spindle Box Zanhui Shu and Qiushi Han

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

MODAL ANALYSIS AND TESTING OF ROTATING MACHINES FOR PREDICTIVE MAINTENANCE: EFFECT OF GYROSCOPIC FORCES

Stability of Serrated Milling Cutters

Technical Report TR05

Research. Mechanical Engineering Article

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

MEEM 3700 Mechanical Vibrations

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Modeling of Dynamic Systems

UIC University of Illinois at Chicago

Global Sensitivity. Tuesday 20 th February, 2018

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Constitutive Modelling of Superplastic AA-5083

Field computation with finite element method applied for diagnosis eccentricity fault in induction machine

Computer Based Porosity Design by Multi Phase Topology Optimization

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

PROPERTIES I. INTRODUCTION. Finite element (FE) models are widely used to predict the dynamic characteristics of aerospace

Multi degree of freedom measurement of machine tool movements. W. Knapp, S.Weikert Institute for Machine Tools and Manufacturing (IWF), Swiss Federal

FUZZY FINITE ELEMENT METHOD

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.

Note 10. Modeling and Simulation of Dynamic Systems

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

Evaluation of the accuracy of the Multiple Support Response Spectrum (MSRS) method

Chapter 11 Angular Momentum

Numerical Nonlinear Analysis with the Boundary Element Method

Survey of applications of discrete vortex method in civil engineering

Appendix B. The Finite Difference Scheme

Advanced Mechanical Elements

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Dynamic Contact Problem for Slide Hinge

VIBRATION OF MULTI-STAGE GEAR DRIVES INFLUENCED BY NONLINEAR COUPLINGS

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Probability, Statistics, and Reliability for Engineers and Scientists SIMULATION

Efficient, General Point Cloud Registration with Kernel Feature Maps

IDENTIFICATION OF DYNAMIC RIGIDITY FOR HIGH SPEED SPINDLES SUPPORTED ON BALL BEARINGS

DYNAMIC MODELING AND PATH PLANNING OF FLEXIBLE-LINK MANIPULATORS

DYNAMIC ANALYSIS OF A CONSTRAINED ROTATING FLEXIBLE EXTENSIBLE LINK WITH SEMI-PERIODIC IMPACT

ES 240 Solid Mechanics Z. Suo. Vibration

Tokyo Institute of Technology Periodic Sequencing Control over Multi Communication Channels with Packet Losses

MODEL OF HYDROPNEUMATIC THREE POINT HITCH

Study on Nonlinear Dynamic Characteristics of 2K-H Planetary Gear System

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Spin-rotation coupling of the angularly accelerated rigid body

HELICOPTER GROUND RESONANCE ANALYSIS USING MULTIBODY DYNAMICS

Effects of System Parameters and Controlled Torque on the Dynamics of Rigid-Flexible Robotic Manipulator

DUE: WEDS FEB 21ST 2018

Rigid body simulation

PLATE BENDING ELEMENTS

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

NON-LINEAR CONVOLUTION: A NEW APPROACH FOR THE AURALIZATION OF DISTORTING SYSTEMS

«ENERGETIC MACROSCOPIC REPRESENTATION (EMR)»

PRATICAL STATIC CALCULATION METHOD FOR ESTIMATING ELASTO-PLASTIC DYNAMIC RESPONSES OF SPACE FRAMES

Nonlinear Effects in Stiffness Modeling of Robotic Manipulators

High Efficient Dynamics Calculation Approach for Computed-Force Control of Robots with Parallel Structures

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures

Dynamic Similarity Design of Geared Rotor System in Five-Shaft Integrally Centrifugal Compressor

829. An adaptive method for inertia force identification in cantilever under moving mass

Chapter 3: Kinematics and Dynamics of Multibody Systems

Invariant deformation parameters from GPS permanent networks using stochastic interpolation

Finite Element Modelling of truss/cable structures

A Computational Viewpoint on Classical Density Functional Theory

Neuro-Adaptive Design - I:

FEM NUMERICAL ANALYSIS OF TOOL HOLDER OF TURNING LATHE. Aleksandras Stulginskis University

Rotor Noise Modeling Kenneth S. Brentner Penn State University

Modeling acoustic transducer surface waves by Transmission Line Matrix method

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Smoothing and Adaptive Stretching of Discontinuous Fields in Finite Element Applications

Computing Nonequilibrium Conformational Dynamics of Structured Nucleic Acid Assemblies

Static and Modal Analysis and Optimization of Cross- Sectional Area of the Manipulator Arm

risk and uncertainty assessment

3D Estimates of Analysis and Short-Range Forecast Error Variances

DESIGN ANALYSIS OF A RUBBER MOUNT SYSTEM FOR A PUSH- TYPE CENTRIFUGE Carl Nel

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.

Indeterminate pin-jointed frames (trusses)

Chapter - 2. Distribution System Power Flow Analysis

Computational Modelling of the Unbalanced Magnetic Pull by Finite Element Method

4DVAR, according to the name, is a four-dimensional variational method.

Transcription:

Smpack User Meetng 27 Modelng and Smulaton of a Heapod Machne Tool for the Dynamc Stablty Analyss of Mllng Processes C. Hennnger, P. Eberhard Insttute of Engneerng project funded by the DFG wthn the framework of the SPP 199 Parallel Knematc Machne Tools

Outlne motvaton: PKM machne tools and dynamc stablty of cuttng processes modelng of the machne structure modelng of the mllng process analyss of dynamc stablty technques for effcent computaton of stablty dagrams pose-dependent stablty dagrams for the heapod machne tool Paral Insttute of Engneerng

PKM Machne Tool PARALIX spatal assembly of lnks hgh statc stffness no bendng moments n lnks, only normal forces lght-weght desgn of struts IfW, Unversty of Stuttgart equal components n dentcal lnks reducton of producton costs 6-DOF heapod machne tool struts wth constant length sledges drven by lnear drect drves 5-DOF hgh-speed mllng of metal Insttute of Engneerng nonlnear knematcs and dynamcs nonlnear poston control necessary pose-dependency of stffness force/velocty transmsson and dynamc propertes knematc sngulartes loss of control of DOF, reducton of workspace

characterstcs of PKM: lght-weght structures pose-dependent dynamcs Dynamc Stablty of Cuttng Processes for PKM? eternal ectaton unbalance ectaton nerta forces process forces machne structure process relablty dynamc stablty surface qualty dsplacements mechansms of ectaton eternal and self-ected vbratons Insttute of Engneerng process process parameters spndle speed, feed, depth of cut,...

machne structure jont elastctes Elastc Multbody System of PKM Paral rheonomc jonts tool + toolholder elastc struts k m c y y z k y c y sledges, platform as rgd bodes fnte element modelng of struts (beam elements) elastc Cardan jonts and lnear jonts modal reducton of strut model (9 dof) Insttute of Engneerng modal reducton of complete model for each pose n workspace

Mode 1 y-z Mode 2 y-z Egenfrequences of Machne Structure n Workspace Mode 3 y-z Mode 1 α-β Mode 2 α-β Mode 3 α-β Insttute of Engneerng

measured tool-tp frequency repsonse Role of Machne Structure n Tool-Tp Dynamcs bmodal system representng machne structure and tool-toolholder-spndle Insttute of Engneerng

nfluence on egenfrequences and tool-tp frequency response Pose-dependency of Dynamc Behavour Insttute of Engneerng

two-dmensonal mllng process Modelng of Mllng Process Insttute of Engneerng lnear cuttng force model F = F k bh F = F k bh t t + ct r r + dynamc depth of cut wth regeneratve effect cr [( ξ(t) ξ(t τ) ) sn( ) ( μ( t) μ(t τ) ) cos( ϕ) ] h = h( ϕ) + g( ϕ) ϕ dynamc process force equaton ( y(t) (t )) F = F ( t) + bkc(t) y τ overall system consderng machne structure and process dynamcs M && y( t) + D y& (t) + ( K + Kc ) y(t) Kc y(t τ) = F (t) + delay-dfferental equaton (DDE) wth perodc coeffcents dynamc stablty?

Insperger & Stépán, 23: dscretzaton of system hstory Analyss of Dynamc Stablty by Sem-Dscretzaton Method homogeneous perodc tme-varant DDE & = A( t) (t) + Q(t) (t τ) appromatng ODE on [ ] & t t + 1 ( w b w ) = A ( t) + Q m + a m+ 1 appromatons on (t τ) w b m + (t w + Δt / 2 τ) a [, ] t t + 1 m+ 1 A ( t) A(t = A + Δt / 2) : + Δt / 2) : Q ( t) Q(t = Q + 1 = P + war m+ 1 + wbr m wth P = ep( AΔt ) Insttute of Engneerng R = (ep( AΔt) I) A 1 Q

Insttute of Engneerng Stablty Analyss transton matr for t Δ = + + + m 1 m 2 1 T b T a 1 m 2 1 1 ~ ~ ~ ~ w w ~ ~ ~ ~ I I I C R C R C P M L M M O M M M L L L L M transton matr for system perod T 1 2 k 1 k Φ Φ Φ Φ Φ = L Floquet theory for stablty of perodc tme-varant lnear systems ( ) < => 1 11 ) eg( ma Φ unstable stablty boundary stable ~, C = vector of delay-actve states, d n C n d <

Stablty Boundary Dagram dynamc stablty dependng on process parameters spndle speed N and aal depth of cut a p up-mllng process of alumnum alloy wth sngle-bladed cutter, a/d=.5 f [Hz] 722 D [-].87 m [kg].2 unstable k ct [MPa] 644 k cr [MPa] 238 Insttute of Engneerng stable

S mplctly defned stablty boundary ( N,a ) ma eg( (N,a )) { ( ) 1} = Φ p p = Effcent Computaton of Stablty Dagrams full dscretzaton equdstant dscretzaton of complete doman computatonally epensve curve trackng local search method adaptve dscretzaton hgh computatonal effcency etra treatment of cusps and near-branch zones necessary dscretzaton of dsconnected regons not possble Sptzkehre Engstelle Insttute of Engneerng

Insttute of Engneerng Curve Trackng

Comparson: Curve Trackng vs. Full Dscretzaton full dscretzaton equdstant grd 5555 ponts 325 functon evaluatons curve trackng (2952 functon eval.) full dscretzaton (325 functon eval.) curve trackng 27 curve ponts 2952 functon evaluatons equvalent to full dscretzaton wth 2525 = 625 ponts Insttute of Engneerng

dea: ntal dscretzaton wth coarse mesh successve refnement of subspaces beng crossed by curve Subspace Dvson advantages: robust algorthm dscretzaton of dsconnected curves possble eample: base dscretzaton: 57 ponts 7 refnement levels 666 functon evaluatons equvalent grd: 257385 ponts (98945 ponts) Insttute of Engneerng

nfluence of rotaton on dynamc stablty Pose-dependent Stablty Charts Insttute of Engneerng

nfluence of translaton on dynamc stablty Pose-dependent Stablty Charts Insttute of Engneerng

Conclusons strong pose-dependency of dynamc behavor at tool-tp for parallel knematc machne tool Paral multbody system of machne structure can be used for analyss of dynamc stablty of cuttng processes mass dstrbuton between machne structure and tool-toolholder-spndle s responsble for rato of resonance ampltudes methods for effcent computaton of stablty charts have been developed Insttute of Engneerng dynamc stablty of process s nfluenced by pose of tool platform n workspace