Tooling design for ecm

Tooling design for ecm V.K. Jain* and P.C. Pandey** This paper critically reviews existing analytical models, ie classical ecm theory, the cos 0 method, the analogue method, the finite difference technique, the complex variable technique, a perturbation method and nomo, graphic and empirical approaches to metal removal in ecm. A review of the progress to date made in the fields of application of the finite elements technique to several electrochemical operations, ie ecm, ec drilling, boring, wire cutting etc is also presented. An introduction to the application of three dimensional finite element analysis to ecm problems is included. The paper suggests areas where future work could improve the accuracy of exis'ting models pertaining to ecru tool design

Electrochemical machining (ecm), described 1'2's as the process of accelerated, controlled anodic corrosion, is achieved by a high velocity electrolyte flowing between the anode and cathode (Fig 1). The electrolyte serves multifarious purposes, including removing the heat generated during machining 4, allowing a high rate of metal dissolution, removing the reaction products etc. Rate and direction of anodic dissolution in the process depend on the applied voltage, electrolyte conductivity, current density, anode material, electrolyte flow velocity, presence of anodic film'etc. Over the years, ecm has found applications in a number of practical machining operations such as turning 6-9'116, trepanning, broaching l°'12s, grinding HA2'I21 fine hole drilling ~3, die sinking ~, cavity sinking ls'~6 , piercing 17'is and deburring tg'sg, and is being used widely in the aeronautics, nuclear 3°, biomedical a and space industries 22. Some typicalexamples of ecru applications are: machining of turbine blades 3z made of high temperature and high strength alloys 21'2a-2s, copying of internal and external surfaces, cutting of curvilinear slots ~4'26,machining of intricate patterns s'2°ns, production of long curved profiles 27, machining of gears2s'29, production of integrally bladed nozzle-rings for use in diesel locomotives a~ and machining of thin (80/Jm) large diameter diaphragms. The process has several advantages to offer, eg machinability of the work material is independent 32'33'37-40 of its mechanical properties, production of stress and burr free surfaces, comparatively good surface finish 34-36, better corrosion resistance of the parts produced, absence of swarf or chip disposal problems, low overall machining times, and working accuracy independent of the operator's skill 36. Many research data have been published, but the potential of ecm as a production process has not been explored fully due to practical difficulties faced in the design of tools, the complex nature 71'72 of the process of metal removal and the difficulties encountered in prediction of the anode shape. This paper attempts to highlight some of the problems of tooling design in ecm and the prediction of anode shape. *Mechanical Engineering Department, MNR Engineering College, Allahabad 211004, India and **Professor, Department of Mechanical and Industrial Engineering, University of Roorkee, Roorkee-247 672, India

P R E C I S I O N E N G I N EE R I N G

Tooling design in ecm Tool design in ecru is mainly concerned with the prediction of anode (work) shape obtainable from a tool while operating under specified conditions of machining. Conversely 64'6s'7°, it also deals with the computation of the tool shape which under specified conditions would produce a work piece having a prescribed shape and accu racy 99'z=s. C o m p u t a t i o n o f anode shape

The workshape obtainable from a given tool in ecm under ideal conditions can be computed from the corresponding equilibrium gap. However, these results are influenced by a large number of parameters such as the presence of an anodic film a6'42-4a, electrolyte flow rate 3a, work material microstructure, intergranular attack s°'s~, change in valency 33 of the work material during cutting, type of electrolye s2-e°A26 and the role of additives 4a'49, current density, electrolyte throwing =2s power, electrolyte pressure distribution, stray current attack and passivity s2'sa's9. In some cases, the exact influence of these parameters is not clear. Fig 1 shows the ec hole sinking operation with the interelectode gap (lEG) divided into four different regions. Most of the work 67-s7 available in the literature pertains to the front and side gaps (zones 2 and 4) only, whereas scanty or no information is available about the mode of material removal in the transient 7L~z'l°° and stagnation regions 87'ss (zones 3 and 1 ).

LElectrolyte

t

/

/

Sol i;,

aQ ~ , // / ~"/ /~_// / ~)~.-~_~ /

Transition

I-Stagnation

Fig I Schematic diagram of ecd with outward mode of electrolyte flow. Zone 1 (stagnation), zone 2 (front) zone 3 (transition) and zone 4 (side)

0141-6359/80/040195-12 $02.00 © 1980 IPC Business Press

195

For the case of machining with plane parallel electrodes, the local metal removal s9 and the equilibrium gap =° are given by Equations (1) and (2) respectively:

In Equation (4), Hopenfield and Cole2° assume a value of n of 1.5 for the case of uniform void distribution; for nonuniform distribution, especially when bubbles are concentrated near the cathode, a value of n of 2.0 has been suggested~6'63'69'sl. Experimentally, the void fraction (~v can be evaluated4'20'Ts'Sz'93 from:

El c t rn =

(1)

F EvEK~

Ye

F F F Pm cos 0

(2)

C~v- 1 + ~

Equation (2) is valid in zones 2 and 4 only and when nonpassivating electrolytes are used. For a passivating electrolyte, E v is as follows Ev = A v - A Ev

where = Qg/QI and

(3) Og -

The overpotential 4,/X Ev, in Equation (3) consists of reversible, concentration and chemical overpotentials and can sometimes be high 1° (up to about 5V). Equation (2) above is based on the assumption that electrolyte conductivity K remains constant during the operation• However, K has been known to be a function of electrolyte temperature, void fraction, distribution of voids within the l E G ~6'63'69'sl , electrolyte concentration, efficiency of electrolyte filtrationS6 etc. For the computation of K, Tipton ~ proposed the following equation which accounts for the effects of electrolyte temperature and hydrogen liberation: K = K o (1 + a AT)(1 - OZv)n

(5)

/ c Rg (273+T+A T) ZFP

The temperature rise4 z~Tof the electrolyte in ecm is due to chemical reaction in the bulk of the electrolyte, the reversible heat of reaction and the viscous heat generated within the electrolyte• However, the major factor that contributes to temperature rise is the Joule heating: dT dx

E2v K - V Pe C e y2

Ac K

which, when solved, yields 1 A T = ( T - T o ) = -d

[Exp {(A c Ko o~x) ( 1-O~v)n} -1]

(6)

(7)

The electrolyte temperature at any point within the l E G depends upon the length of the flow path, x. Equations (2),

(4)

Nomenclature A

Cross sectional area Applied voltage t Interelectrode gaps (See Fig 1) Bg,a ,as ,a s Constant Bare length of the electrode bb bc,Bl -Bs Constants Specific heat of the electrolyte Ce Machine tool-rate Cm Cost of drilling tool Ct Unitary matrix Dm E Electrochemical equivalent Effective voltage Ev F Faraday's constant Front feed rate FF Current Ic ! Function lEG I nterelectrode gap J Current density Stiffness matrix Km I( Electrolyte conductivity K,j Coefficients of stiffness matrix k Wave number L Length of electric resistance m Metal removal Total number of sparks before the end of toolms life mr/" Metal removal rate n Exponent n r Length of normal to anode surface Number of sparks per unit length of drilling ns Mean gap between electrodes P P Electrolyte pressure within l E G Q Electrolyte flow rate R~p Gap resistance f Radius A v

196

T t tc V v x W Y Z AE v Po p 0

o~v v,v' e ¢ X 77

Temperature Machining time Time to change the tool Electrolyte flow velocity Volume Distance along electrolyte flow direction Width Interelectrode gap Valency Overpotential Volumetric concentric concentration of hydrogen Density Angle of inclination between feed direction and normal to the tool surface Temperature coefficient of electrolyte conductivity Void fraction Correction factors Amplitude of surface irregularities Electric field potential Wavelength Machining efficiency

Subscripts

c d e F g i,j

I m

s t V

corner drill electrolyte (Ce ), element (H e ), equilibrium front gas node number liquid work material sparks, spike, side condition at time t voltage

PRECISION ENGINEERING

(4) and (7) predict that the lEG would become nonparallel after a time t even for the simplest machining configuration. Theoretically, the ec reaction would attain an equilibrium only after t approaches infinity and the corresponding lEG can be computed from Equations (8) and (9) while machining with zero and finite feed rates, respectively. Y = (Y2o + 2 C A t) ½

shapes, the electric field gets distorted and in such cases Ohm's Law is not obeyed 27 . Under such conditions J can be computed from 37-39 : J=-~ (2K o

Ippolito 1~7 has also mentioned that the conductivity of the meters and hence application of Ohm's Law to ecm is not justified. According to him

(8)

Rgap

Ev E K Ye z~t=-~F [Yo-Yt + Ye In ( -o "e) ] 1

Y

-

Re =

(9)

Yt- Ye

Y= Yo+C' At

(10)

Ye = Yo + ( C ' - F F) A t

(11)

where

C' = T?J/F Pm and z& t is very small Equations (10) and (11) yield values of Y which are close to those obtained from (8) and (9) and agree well with the experimental9~']°~ data. Assuming that the electrolyte behaves like a pure ohmic resistance, the current density J can be computed from Equation (12).

y

1 Va QbpC

; Radd = ~ K / c ~c

(15)

(16)

For a simple shaped lEG, Equation (16) would give satisfactory results; however, for higher accuracy and when a complex shaped lEG is involved, the conformal mapping technique is recommended. O p t i m i s a t i o n o f process variables

(12)

Rgap

A~

Vl W Y1 = V2 W Y2

Ev -

(14)

and the values of the constants have been evaluated experimentally nT as a = -0.32, b = -0.12, c = -0.11, ~ = 0.151. Current density distribution within the lEG has been measured by Kawafune ~s who found it to be lower than the analytically calculated distribution. However, in some cases~ , experimentally measured current density has been reported to be higher than that obtained analytically (see also Fig 2). The conflicting nature of the results can be explained on the basis of actual anodic dissolution efficiency, r/, achieved and a change in anode material valency with time. lEG varies along the electrolyte flow direction, with the result that the flow velocity also varies and can be evaluated from the equation of continuity for flow of electrolyte, ie

Equation (9) is implicit and can be solved iteratively. Equation (8) predicts an infinite value of the gap as t approaches infinity. However, in practice, as the lEG increases, the current density decreases so that the mrr (metal removal rate) gradually diminishes and for a large gap the process would come to a standstill. Equations (8) and (9) have been re-derived ]13 with the lEG expressed in terms of the current density, J:

Ev K

= R e + Rad d

where

C=11 F p m

-

(13)

lEG is a function of the machining current and other para-

where

J

)

During ecm, at certain values of the limiting tool feed rate, the process 16'8° can change from boiling to non-boiling and choking to non-choking. Feed rate greater than the

Equation (12) is valid within a certain range of experimental conditions only 2° . While machining complicated anode

/ 10 ~--

x = 84mm

L

f

Machining time

x = 57mm Current 2"3

Current 2"4

5"

o

0

~"

£~

T

o

o

0

t=0S

~3

•

o

t = 140S

o

n

0

t : 320s

o [I

I I I f II

I

c

I 0 I--

L

3

-t r t f t

I I I I I I I I

x = 2.Smm

--

Current Z I

--

it

I

I

III

x = 30mm Current Z 2

b

x

5

Oo

o

000

°0

0

o

o

9

O

t = 40s

o

o

O

o

t =60s

_,', ~, z~ ~ ,', ,,, ~, ~ ,

I I o

II

IOO

200

3oo

o

Machining t i m e , s

I I

1111 IOO

I

II

I I

200

Machining t i m e , s

Zero feed rote Experimental x, c,,o, (3 F E T - II - -

V = 21.0m/s x = Distance

I

I

--

~z~ ~, z~ t : 3 2 0 s

Machining accuracy

II 5

I

I II

,c E

2.02

2.00

3OO

I IIII

% 2.04

1.98

Ill

10

E ~

14

Distance, x , m x I0 - z

Fig 2 Comparison of analytical and experimental results

PRECISION ENGINEERING

197

value corresponding to the intersection point m' of choke limit and boil limit (Fig 3) implies that the flow in the lEG is either choked or boiling. Maximum mrr is defined to correspond to the point m'. Optimisation analysis carried out by Bhattacharyya62 with mrr (or feed rate for constant area of machining) as the objective function and the following three constraints Vo i> 1.65 F F (Electrolyte boiling)

(17)

Vo/> 1.44 F}: (Hydrodynamic instability)

(18)

Vo ~< 2.66 F F (Passivity)

(19)

predicted an optimum feed rate, in ecm, of 1.825 mm/min. This analysis does not include the onset of sparking as a constraint. This is important from a practical viewPoint as beyond a certain frequency of sparking (or feed rate) the tool damage due to sparking could become significant and make the process uneconomical 73'a8. Higher feed rates in excess of the optimum lead to lower tool life. Larsson and coworkers 7a suggest that the feed rate given by Equation (20) would result in minimum machining cost/part. ms o.5 F F ={b c

(tc+~,t/cm)}

(20)

For the condition of zero sparking rate, Larsson and co-workers ~a's7 have derived an optimum feed rate of 2.28 mm/min which means the optimum feed rate62 of 1.825 mm/min takes care of the sparking constraint automatically. However in general, high mrr in ecm can be realised by increasing the feed rate, supply voltage, electrolyte pressure, flow velocity and the electrolyte inlet temperature. To achieve higher feed the quantity J r/should be high, but increase in J leads to a decrease in r/. On the other hand, if Y is decreased beyond a certain value, increase in potential gradient would result in arcing and the process would come to a stop 76 It has also been reported ~27 that, in some cases*,the sludge produced during the process could coat the electrodes resulting in a complete halt of the anodic dissolution and short circuiting of the lEG. From the discussion above, it is evident that the relationship between the'parameters governing the ecm process is complex and hence difficult to model analytically. Some progress made in this direction is discussed below.

have been used extensively. An alternative method for tool profile correction is to calculate the working gap and account for it during tool manufacture. With the availability of high speed computers it is now possible to determine the tool shape, that would yield a workpiece of specified geometry, by successive approximations. Since the computing time is only a fraction of that required for carrying out the experiments, applying tool correction analytically has proved to be time Saving and economical. A number of models for anode shape prediction and tooling design in ecm have been developed. Since the ability to predict the variation in lEG for any given operating conditions is a pre-requisite of proper design of ecm tools, these models are discussed in terms of the equilibrium gap. Cos 0 m e t h o d

This method is based on the computation of equilibrium gap for the given conditions but excludes the consideration of the mode of e!ectrolyte flow, overpotential, variation in electrolyte conductivity, heat transferred to the environment etc. In this method, equilibrium work shape is computed corresponding to the tool whose profile has to be approximated by a large number of planar sections inclined at different angles. This method was first proposed by Tipton 6s and is described in detail in reference 92, but its scope is limited by the following Regions with sharp corners 92 cannot be analysed. Generally applicable only for 0 < 45 ° It is not possible to account for the effects of the mode of electrolyte flow, void fraction, change in electrolyte temperature and conductivity, overpotential, heat conducted away to the tool and work etc. Researchers have also expressed conflicting opinions about the reference surface for the measurement of the angle 0. Tsuei 7° et al suggest that 0 should be taken as the angle between the tool-feed direction and the normal to the anode surface,and not the cathode surface, as suggested by others 6s . However, in view of the approximations involved, its use is not recommended, especially when complex shaped workpieces are to be analysed.

A n a l y t i c a l m o d e l s f o r t o o l design in e c m

Finite difference technique (fdt) :

The classical theory of ecm is valid under ideal conditions only. Trial and error methods for tool profile corrections

The cos 0 method is based on the assumption that the lines of electric potential are straight and normal to the electrode surfaces, but this is not true when electrodes with small radii are used. In such cases, it is necessary to know the electric potential distribution in order to determine the current density, mrr, transient and equilibrium anode profile etc. fdt has been employed by several investigators4'4~'64'7°'83-8s'~19 for ecm tooling design using a nonpassivating electrolyte with constant conductivity and temperature. In such cases, the electric field flow lines are governed by Laplace's Equation (21) and the boundary conditions given by Equation (22)

4

-

Theoretical

,,

.c_

_~ g

3

fn I

I

o Run 2 z~ R u n I

I I o 0.06 o.12 0.18 Feedrote.in/min Fig 3 The ecm characteristic curve'S°.m' corresponds to the point o f maximum metal removal rate

198

a2~ a2~ ~2~ . +--+--=0 ~)x2 ~3y2 ~)z2

(21)

~b= 0 at the cathode and ~=Ev=Av-AEv at the anode

(22)

Equations (21) and (22) would yield a set of simultaneous equations that could be solved for ~ by the backward differences, forward differences, or central differences9s'96

PRECISION ENGINEERING

technique. Corresponding to the potential distributibh obtained, the instantaneous current density s3 at the anode surface can be evaluated from ~ r l l l ] l l l l l l l l l

a ~b) J = K ( an---;

(23)

i i i J ~ l l ] l

........

yJ Fig 4 shows the tool-wo'rk surfaces and lEG drawn in square mesh. The initial potentials at the grid points within the lEG region are set by linear interpolation along the vertical gridlines betweeen the tool and work boundaries. For a point 0 located in a mesh of spacing h (i,j) we can write s3

@ii= (@i+1,i + @i-1 ,j+ @i,i+1 + @i,i-1 )/4

I

,. ,. ,. ,.

i l r l l l

iJl

'~'

Ii ~ iI i iI I Ji Ii I[ ii l l i i

. . . . .

0 = B l x 2 + B2Y2 + B3x + B4y+Bs

(25)

Analogue m e t i l o d Laplace Equation (21) was first solved by the conducting paper analogue s3'9~'97. In this case, equipotential surfaces representing the anode and cathode, to an approximate scale, have to be drawn on a conducting paper. The work-" boundary is segmented to evaluate local current density and then it is moved to the position they would occupy at the end of the time interval At. This could be obtained by finding point to point movement by the vectorial addition of feed rate vector (F F /k t cos ~) and cut vector

PRECISION ENGINEERING

J

;4

::''"

j-li_ I

il

i+l

~x

a

b 3 0

4

I

4

/"z I

After each computational interval of Ats, work and tool boundaries are moved to new locations, according to the cut and feed vector. For simplicity, Tipton sl has only accounted for the vertical component of the cut vector which is proportional to the vertical potential gradient at the work boundary. However, for precision in the results, movenlent of the cut vector in both x and y directions should be considered. It is thus evident that fdt would yield approximate results. Further, in case of a complex shaped lEG, the tool work boundaries cannot be matched accurately using square mesheswhich introduce still greater approximations. On account of the approximation involved in the use of this technique, the authors recommended 92'u° the use of the fet for tool design in ecm.

/

j2,

~ i l ~- -.......

(24)

Equation (24) shows that the potential at any point 0 (i,j) on the square mesh is equal to the mean of the potentials at the four nearest adjacent points (ie the points 1,2,3 and 4 in Fig 4). The procedure of computation in this case is to consider each grid point (i,j) within the field in turn and adjust its potential value to the mean of those at the four points around it and to repeat the process until the potential values attained are correct within the prescribed tolerances. It should be noted that the work and tool boundaries are at fixed known potentials and, therefore, their values are not adjusted during the relaxation process. In some cases, all the points on the tool and work boundary may not lie on grid points (Fig 4) and the regular stars may not be formed; such problems have been attempted by the use of the over-relaxation technique s3. To solve such problems, Nanayakkara and Larsson n9 have suggested the use of irregular grids along with regular grids. In this case, a polynomial Equation (25) was used instead of a linear interpolation Equation (24) to describe the potential distribution in the area around a nodal point including its near neighbours.

h1

[llll

"3 I

d

c

e

I

Fig 4 Finite difference analysis model s3 of ecru process (a) I EG-tool-work system discretised in square meshes (b) A square mesh of spacing h having four equiplaced adjacent points 1,2,3,4 which form a regular star. (c) Boundary forms irregular star with one short arm-03. (d) Boundary forms irregular star with two short arms 02,03 and (e) Boundary forms irregular star with three short arms 02,03,04 ( K i ~ , ) A t). The method is based on the assumption that feed rate velocity and cut velocities are constant over the time A t. This transforms point A to point B to make a new surface on the work and this process is repeated until all points A are transformed to points B which are on the worksurface. Before repeating the process, the applied voltage must be adjusted so that the same current values as used earlier are obtained. The process has to be repeated until the workshape does not change appreciably between two successive steps and the final equilibrium shape obtained s3. The analogue method is approximate and cumbersome, its accuracy depends on the skill of the operator and is not advisable for use when a high degree of precision is desired.

Empirical formulae The exact path of the electric current flow lines, within the lEG, is difficult to determine analytically. Therefore, normally, the chordal distance between two stations27 is taken as the length of current flow line and this assumption in the majority of cases is responsible for the discrepancy between the analytical and experimental results. It has also been found that the conformity of the surface radii of the tool and the anode cavity decreases as the angle 0 increases. Therefore, attempts 22,3~,71,1oo have been made to derive empirical equations for the evaluation of the

lEG. Based on experimental data, empirical Equation (26) has been suggested by K6nig and Pah122.

ao = r 0"35 0.35 [(10e*) ye] 0.5 c

(26)

for 0.15 ~< Ye ~ 0.6 mm and 0.5 ~
199

K6nig and Degenhardt 7~ also suggest Equations (27-28) when the bare length of the tool (b b) varies between 1 mm ~< r c ~< 5 mm. ao = (0.1 + Ye)(0.314 r c + 1.17) for b b = 0

(27)

ao,= 2Ye + 0.1 [6.283(rc-1)] 0"5 for bb/> 1

(28)

Equations (27-28) have been found 91'127 to yield erroneous results at low feed rates (F F <~ 0.006 mm/s) or when the equilibrium gap is large ( Ye ~> 1.0 mm). Based on the regression analysis 127 of experimental data plotted in Fig 5 (mild steel anode, brass cathode, electrolyte 10% w/v NaCI and feed rate 0.0057 mm/s, outward electrolyte flow, coefficient of correlation 0.8) the authors suggest Equation (29) for the computation of a n, mm. ao = 0.19536 Ye + 0.5779

(29)

The magnitude ofa's (Fig 1) can be evaluated 22 from Equation (30), as t

,

2,0

5

as = (2 b b Y e t aol •

(30}

!

Analytical evaluation of as, which is affected by stray current, is difficult. It can, however, be evaluated =z from the empirical equation

gaps, Ye can be used as a basic parameter. In some of the cases the side gap has been demonstrated to be independent of the machining time which, however, is not true in practice. It is to be noted that in the majority of cases the effect of electrolyte flow mode (inward, outward etc) has been neglected. Empirical equations are normally valid under specified working conditions only, which limits their use.

Nomographic approach For the evaluation of equilibrium anode shape, a nomographic approach has also been used22:1'1°1. K6nig 22 has prepared a nomogram for the evaluation of side gap for the known equilibrium gap, tool-radius and the bare tool length. Fig 6 shows such a nomogram for bb=0 and bb/>0 for a given tool. Heitmann a4 has used the nomograms for evaluation of electrolyte temperature rise (A T) while working under specified conditions. Use of nomograms has proved to be advantageous in planning an ecm operation. However, the nomograms are very often based on a number of simplified assumptions 8,9 and hence cannot be recommended for general use; they would be more useful if available in a generalised form.

a s = [2 b b Ye + rc0.7 0.123(10e*) Ye] 0.5 + 0.65 Ye

(31)

Test results show that the overcut in ecm is a function of the machining parameters but does not depend on the tool di mensions22 . The authors' experimental data 9~ , however, point to the contrary. Further, it should be noted that Equations (26-31) are valid only for the case of non-passivating electrolytes. I ppolito and Fassolio l°° checked the validity of Equation (28) and report that for rc=0 it gives the upper bound whereas Equation (32) gives the lower bound of the side gap; Equation (33) was found to cross the experimental data. I

as = (2 b b Ye +/)2 y2}0.5 -e"

(32)

t

as = (2 b b Ye + 2.9 Ye2)0"5

(33)

The lateral gap at the end of the non-coated zone of the electrode can be computed from: a's = (2Yeb b + (1; Ye*)2 - c ) 0"5

(34)

where

Anode shape prec~iction in ecm has been attempted by employing the complex variables approach for the case of shaped workpieces and using simplified assumptions. Collett et al 1°3 have determined lEG for both completely side insulated and bare and straight sided tools and arrived at Equation (36) Overcut at corner = 1.159 Machine gap

=

1.16 mm/min

From the results given in reference 100 it is evident that Equation (34) is valid only for a specified tool-work combination and specified machining conditions and the discrepancy between experimental and analytical results increases with an increase in K Ev/F F. Equation (35) also predicts variation in overcut with the conical taper and it was found l°° to be practically independent of the uncoated length b b. ,

2B'Y

as=as+ ~

e

,

larctan

2x+B' ~r

- arctan (BI/D')]

(35)

where

B'= v'E K E v ; D' = (4F~= a's2 -B'2 )0.5 FPm v' is a correction factor of conductivity and takes into account the fact that the lines of electrical potential are not straight and the work-surface will be lightly passivated by the presence of the metallic oxides and hydroxides. It is thus evident that for the evaluation of side and front

(36)

However, PERA 118 suggests a value of 1.7 for this ratio. This is because Collett et al 1e3 assumed that the electrolyte conductivity and void fraction remain constant throughout. Hewson-Browne 1°4 extended the work of Collett et al and used the conformal mapping technique for the analysis of two dimensional machining using straight sided tools with a finite land width b b. For an insulated tool it was assumed that the sides do not participate in metal removal; it was shown for such a case that 0.731 ~ao/Y e ~ 1.159

¢$--

Ye-l.29Ye, c = 0.29 and F F

200

Complex variables approach

(37)

where the lower limit applies to the case of a side insulated tool and the higher limit is for a bare tool. However, the authors have found experimentally that for the case of machining with bare tools the ratio ao/Ye is a function of machining conditions and tool-work combinations. Nilson and Tsuei41'86'94 have solved the problem of anode shape predicted by the inverted method, in which the spatial coordinates were selected as dependent variables on the plane of complex potential and the transformed boundary conditions are known explicitly on the free boundary. Known asymptotic solutions for the side gap region reduce the size of the field which must be determined. Use of mathematical models has not so far been generalised because workpieces have complex shapes and there are many parameters to consider.

A perturbation method This analysis has been developed by McGeough ani:l co-workers 77-79 and can be used for the analysis of ec

PRECISION ENGINEERING

operations like deburring, anode smoothing or shaping when the amplitude of micro-irregularities on the cathode and.anode is small compared to lEG. This analysis has limited applications and is based on the following simplified assumptionsS9 : Electrolyte conductivity is constant Effects of Joule heating and hydrogen gas bubbles are suppressed by sufficient agitation o.f the electrolyte Machining efficiency is 100% Perturbation boundaries change so slowly that their motion may be ignored at any instant. Consider a case in which cathode and anode shapes (Fig 7) -are defined mathematically by Equations (38)a-b respectively

e=eoEXp[_MEvk~ c o t h k s ds] (41} Po M Ev-FFS dp provided d-~ =/=0. Here, M = E El(pro F) and s is the coordinate

Y=0

(38)a

Y =p + e sin k x where e < < p

(38)b

that specifies the average surface. Equation (41) has been applieds9 to analyse irregularities on the anode shape of the type of short wavelength, long wavelength, arbitrarily shaped and even and arbitrarily shaped. This analysis is based on the solution of Laplace's Equation (21) and therefore accounts for the field concentration effects. This explains the more rapid smoothing expected with short wavelength irregularities which have a much greater field concentration than larger wavelength irregularities. This analysis uses the Fourier series expansion and has been extended to analyse a more general class of electrode shaping by Fourier Transforms. An equation similar to (40) has also been derived by Tipton 6s in terms of dimensionless quantities for the case of surface smoothing to the required tolerance.

Laplace's Equation (21) can be solved with the boundary condition given below. ~b=0on Y = O a n d ~ = E v o n

Y=p+esinkx

It has been shown s9 that the potential at any point Y can be obtained from

= E v Y Eve s i n k x s i n h k y P

p

Finite elements technique (fet) In view of the shortcomings in the analytical models discussed so far, the authors recommend l°s the use of the fet where the choice of the shape and size of the elements is convenient ~°6, it is easy to incorporate different boundary conditions 1°~ and it is possible to analyse non-homogeneous situations 9s . In this case, the effects of simultaneous variations in different parameters could easily be accounted for. The authors have so far developed one, two and three dimensional models for the anode shape obtainable in ecm.

(39)

sinhkp

Using Equation (39), a set of equations has been derived which describe separately the behaviours of the anode surface and irregularities respectively: 1

A t = FFF [P°-P+Pe in

Pe-Po]

(40)

Pe-P

18 -Unidimensional x

14 --

A

ro ,o

A

E ~o

o

O

S

~ ~"

--a%~

~

~

0 0

6

0

I

I

I

I

0.5

I

1.0

1 1.5

I

I

I

2.0

1

l

I

2.5

:3.0

I

analysis

(Model

FET-11)

For the analysis of the unidimensional electrochemical wire cutting process (ecwcp) l°s'z°9 and the ecru process92'zl°'ln the lEG has been discretised into a number of elements (Fig 8). Electrolyte temperature at any nodal point in the continuum, ie lEG is evaluated from (6) which could be written as I

:3.5

I

51=dT 6Tfj -

4.0

dx

xi

~j dT ~xi

(6

d_TT) dx dx

(42}

= 0

~, mm

2.5

r, mm :3.03.54.0

yC x

~mm~

Cathode

/ . z / i / / / / / / / / / / / / / / / / / / / / / / .

"~

p(,I

2.0

Electrolyte

I'--

)' - ' t

1.5 1.0 I

----I=

0.5

Fig 5 (Top left) Relationship between equilibrium gap Ye and overcut, ao

-

b,.~"~ I

0.7

I

0.5

Ill

[

410.:3

o.s

Fig 6 (Left) Nomogram 22 for the determination of the side gap, as

t

0.1 0

0.6 0.4 0.2 Og, mrn

PRECISION ENGINEERING

0

0.2

III I I ~1~).4 0.6

Fig 7 (Above) Configuration of electrodes s9

201

On simplification, this reduces to X = D H e.= K m Te

(43)

where D =~

i

1

; He =

Pe y; Ce V and x•

Km = f l

BTBdx = [

1 -1

1

--1

xi

] for an element

(44)

X2

=

2

X3

-

-1

(45)

T~

Two dimensional analysis

fet using triangular elements can be used for the analysis of two dimensional machining problems assuming that the temperature rise over an element follows a linear law. But in the case of non-parallel lEG the electric current flow lines would not be normal 84 to the surfaces and the electric field potential in such cases would be governed by Laplace Equation (21). The authors H3'H4 have developed

k

a

Work

~,

I

Work

'- -

k

V

2 X

C

Fig 8 (a) One, (b) two and (c) three dimensional discretisation of interelectrode gap

202

I= fffv 1 {g~T

(a0)2] d x d y d z oz

[ D m ] { g } dv

(46)

(47)

09 09 o30] wheregT= [~x Oy Oz

(48) !

I(e)

/ = / ( 1 ) +/(2) + / ( 3 ) + . . . + / ( E ' ) = ~

(49)

e=l

where I (e) is the contribution of a single element to I. The minimisation of / occurs when

01

-

63{0]

0

E'

E'

Z; / e = : S 0 [ 0 ] e=l e=l

0 I (e)

-0

(50)

0[0]

T o evaluate derivatives a-~-}' 63 / the integral should be w r i t t e n

in terms of nodal values, ie

0 (e) = [N re) ] (0]

(51)

Therefore,

-aN l (o) ax g(e) = a N z (e) ay ON t (°) az

aN 2 (e) ax

aN¢(eF """

01

ax

0 N 2 (e) ay "

a N¢ (e) ay

92

63N2 (e) o~z "'"

0 N e (e) az

0e

=[B] (e)[o] (52)

Here B contains the information related to derivatives of the shape functions. After simplification 66 it can be shown that 8/ a9 - [ K e ] {0}

(53)

[K e] = fffv[Bre)] T [D(e)] [B(e)] dv

(54)

where

a0

X

b

(~-~t +

or

0/ -- = E ~,

Tool I

ox.

The final system of equations is obtained by substituting (53) into (50)

-F%.-~----~---~;~I j'

1 [(a0)2 + 09,2

I=fffv 2

T3

The set of simultaneous equations can be solved by Gaussian elimination technique 122 for unknown nodal temperatures using a boundary condition of A T = 0 at entry. Thereafter. current density, mrr, conductivity etc can be computed. The analytical results have been compared with experimental data and a good correlation between them has been shown. As an illustration, a comparison between temperature computed from model FET-11 and obtained experimentally9° is given in Fig 9. The authors have also conducted a parametric study H2 of the ecru process in attempting to form some guide lines and to evolve a standard procedure for ecm tool designers. This model has also been extended 113 (model SGFET-11) for the anode shape prediction during ec drilling which has been treated as a problem with moving coordinates. It has also been shown that the overcut on the side (a's) is basically governed by the computed value of ao . A comparison between the computed and experimentally measured anode shapes shows 91 a good correlation between the two. However, the degree of correlation could further be improved if more realistic values of the machining efficiency, T/, anode material valency and heat transfer data were available.

i

Three dimensional analysis

Equation (21) can be used for the evaluation of three dimensional potential distribution within lEG. The problem is equivalent to finding the 0 distribution that satisfies the boundary condition and minimises the function /

Assembly of the equations thus obtained for different elements, say three, would result in a symmetric, sparse and banded matrix of the form given below:

IXlt fi 1 !l

a computational procedure (Fig 10) for anode shape prediction (model SGFET-22), evaluation of current density and other related parameters by solving the Laplace's Equation (21). Fig 11 shows a comparison between the theoretical and experimental 9° electrolyte temperature distribution within lEG. Further, improvements in FET-22 and SGFET-22 are possible by the use of isoparametric elements 9s'123 instead of triangular elements.

[K(e !] [0]

(55)

e=l

Evaluation of stiffness matrix K - - m(e) involves the solution of the integral Equation (54) and its value is a function of the type of element being used. Fig 8 illustrates the use of tetrahedral elements for the analysis. Assuming linear interpolation functions, field variable O(x,y,z) can be uniquely and continuously represented in the solution domain by Equation (56)

~e (x,y,z) = [N (e) ] {0 e] = NiOi + NjOj + NkOk + NI~ I

(56)

PRECISION ENGINEERING

Experiment . . . .

FET-II

ko = 0.01412 ~-I mm-S Yo : 0.Smm

To = 30oC A t : 30S t : 500s

_ Ev=2O.OV

--// I

45

O = 23.33 = I03 mmS/s F"F = O.OI5mmIs

E v:12.0 v

I

Readand print inputdata

I

Computeno.of computationcycles

J

/

40

I GenerateK,V and Y for different nodes I I Computefeed rote at everynode I

35

I

Determinebandwidth Initiolise parameters

.= o

30

J

J I

Generateand assemble elementmatrices

J

&

E

Formglobalmatrix 45 [

O : 31.67=103mmS/s FF = O.OIOmm/s / Ev = 16.0V / / /

40 --

35

-

/

/

~

/ / /

/

/

O : 25,0=103mm3/s F"F : O.OIImm/s Ev = 12.0V

Apply boundaryconditions j

_

f

J

Solve set of eqns.to determine

I

I

I Compufe I

//

-

J

I ComputeT I :50 O

,

,

8

16

F

Y" 24

0

Distance, mm

I

I

8

16

~

Distance, mm

Fig 9 Comparison of analytical (FET- 11) and experimental results

x o •

Experiment. . . . FF, mm/s x I0 0.0011 0.0015 0.0020

FET-22 V, m/s x 4.74 o 6.02 • 2.32

Computeheat transferto air

I

Computedecreasein electrolytetemperature

1

Jco~p~e

Fig 10 (Right) Computational scheme used in mode/SGFET-22. KiN is an index which counts the number of computational cycles. NC is the total number of machining cuts (or number of cycles) for which the computation has been carried out V, m/s &528 2.:55 2.74

I 24

effectiveelectrolytetemberatureJ

I. ComputeK,Y,mrrL,mrrvondmr I

FF, mm/s x I0 0.0015 0.0011 0.0020

fes coordinates Outward flow

40

Yes I

Reassignx-y coordinates ~ Side flow

~0 Yes Reassign x-y coordinates inward flow

30

//

I ComputeV for new value of Y

E ¢p :

20

/

E /

I0

/

/

/

/

/

/e

NN= numberof nodes Electrolyte flow for 0

I0 Distance. mm

20

25

O

I0 Distance, mm

20

25

IOWO outwards

Fig 11 Comparison of analytical (FET-22) and experimental results

PRECS IO I NENGN I EERN IG

203

Differentiation of functional / and simplification 66 leads to Equation (57)

K['Kji i Kki

Kij Kji Kkj

K"/[ Kkk

Kkl[

= .(0}-

In the opinion of the authors, it is better to go for numerical methods such as the f d t and fet. The fet is the stronger of the t w o and should develop into a standard procedure for ecm tool design.

1571

Ck

References

KIk K I I j ¢1 LKli KIj where Kii, Kjk . . . . . KII are coefficients of stiffness matrix K (e) that can be evaluated as follows:

1.

Bannard J.E. Effect of Density on the Electrical Conductance of Aqueous Sodium Chloride Solutions. J. App. Electrochem, 1975, 5, 43-53

K(e)_SSSv(a__~N- O~N" ~/Vi ~__N/N- ONi ~NIN'~dxdydz

2.

Idem Electrochemical Machining. ibid, 1977, 7, 1-29

3.

Idem The Use of Electrochemical Machining in Biomedical Engineering. ibid, 1975, 5, 89--90

4.

Hopenfield J. et al. Prediction of the One Dimensional Equilibrium Cutting Gap in Electrochemical Machining. Trans. ASME, 1969, 755-.765

5.

LaBoda M.A. et al. Intricate Pattern ECM on Ferrous Alloys. J. Electrochem. Soc. 1975. 122, 1489--14.91

6.

BergsmaF. Electrochemical Machining o1" Metals. Annals CIRP, 1968, 93-99

7.

Hofstede A. et al. Some Remarks on Electrochemical Turning. ibid, 1970, 18, 93--106

8.

Mukherjee S. et al. On Electrochemical Turning. J.Inst.Eng. (India) 1977, 57,210-211

9.

M,lkherjee S.K. et al. Analysis of Electrochemical Turning. Froc. 8th AIMTDR'Conf, 1978, 554-557

10.

DeBarr A.E. and Oliver D.A. Electrochemical Machining. MacJonalo and Co. Ltd, 1975

11.

Kuppuswamy G. Electrochemical Grinding. Proc. 7th AIMTDR Conf, 1976, 337-340

12.

RanganathanV. Electrochemical Grinding of Titanium. ibid, 165--168

13.

Bannard J. Fine Hole Drilling Using Electrochemical Machining. Proc. 19th Int. MTDR Conf, 1978, 503-510

14.

Chikamori K. et al. Electrochemical Machining in Stagnant Electrolyte. Int. Conf. Prod. Engg. held at New Delhi, 1977, viii, 31

15.

Kawafune K. etal. Accuracy in Cavity Sinking by ECM. Annals CIRP, 1967, 15,443-455 Thorpe J.F. et al. A Theoretical Analysis of the Equilibrium Sinking of Shallow, Axially Symmetric Cavities. Fundamentals of Electrochemical Machining lEd. C.L. Faust). Published by Electrochemical Society Inc, Princeton NJ, USA, 1-39

ij -

\Ox

(58)

+ ov ov +Tz Oz I

Using the definition of interpolation function, Equation (58} can be shown to be reduced as 1

Ki j =

36v2 (bibj

+ c i c j + did/)

(59)

where a,b,c and d are constants for an element and can be computed w~ This analysis w o u l d give three dimensional potential distribution w i t h i n the lEG f r o m which current density, mrr, temperature distribution etc could be computed as done in earlier models.

Discussion and conclusions The mechanism of anodic passivation at high current densities in ecru is not well 36'42-47 understood; it has been found, however, t o reduce the anodic dissolution efficiency. Therefore, attempts have been made to destroy the passivating f i l m during the process by the use of reducing a~jents 39'46 as additives t o the electrolyte. Quantitative data regarding the effects of grain size, grain boundaries, density and orientation, on mrr and current density are not available and hence their effect in modelling cannot be incorporated. Also, governing equations do not exist that can relate the influence of electrolyte mixtures on surface finish °~, dimensional control 43-4s, pH of the electrolyte 124 and mrr, especially w i t h comoosite materials. Further, most of the cases do not account for the change in valency 33 w i t h time. This can often seriously affect the values of machining efficiencies s9 . Presence of lines or striations on the machined surfaces have also been observed experimentally but their basic cause has not been resolved 4° However, remedies, set by experience, have been suggested, eg injection 32 of gas into the electrolyte and use of multiholed electrodes 37-4°. It is evident that the problems concerning ecm tool design are highly complex and involve a number of interrelated parameters, such as temperature, electrolyte conductivity, current density. Further, a number of parameters (effect of microstructure, grain boundary, passivation f i l m etc) whose effects are not well k n o w n , make the process more d i f f i c u l t to analyse. It can thus be concluded that most of the models of the ecru process make use of simplified assumptions, such as constant c o n d u c t i v i t y and neglecting the void fraction which give rise to discrepancies between the analytical and experimental results. In practice, ecm tool designers often resort to an empirical approach which, however, requires extensive research data, including the effect of passivating electrolytes on metal removal. A nomographic approach can also be useful for tool designers and can be developed to a satisfactory level.

204

16.

17.

Mashinostroeniya V. Electrochemical Piercing of Holes in Mono - and Polycrystalline Alloys. REJ, 1974, 54(4), 75- 76

18.

Petrov Yu. N. et al. Electrochemical Piercing of Holes in VK Type Cemented Carbides. REJ, 1972, 52(10), 54-57

19.

Rolsten R.F. and Manning E.J. Electrochemical Deburring. The Tool and Manufacturing Engineer, 1969, 14-- 16

20.

Hopenfield J. et al. Electrochemical Machining -- Prediction and Correlation of Process Variables. Trans. ASME, 1966, 88B, 1-7

21.

Hoare J.P. et al. Electrochemical Machining of High Temperature Alloys in NaCI0~ Solutions. J. Electrochem. Soc, 1973, 120(8), 1072-1077

22.

K6nig W. and Pahl D. Accuracy and Optimal Working Conditions in ECM. Annals CIRP, 1970, 18, 223-230

23.

Idem Electrochemical Machining of Heat Resistant Alloys. ibid, 1972, 20, 43--44

24.

IsaeveA.J. Electrochemical Shaping of High Strength Bodies of Rotation. REJ, 1972, 52(5), 40--42

25.

Meleka A.H. ]-he Application of Electrochemical Machining. Proc. Conf. on Machinability, Iron and Steel Institute, 1965, 221-224

PRECISION ENGINEERING

26.

27.

Teterev A.G. et al. Assessing Errors in the Electrochemical Machining of Curvilinear Slots. Machines and Tooli,ng, 1970, 41 (7), 4 3 - 4 6 Moir P.J. et al. Electrochemical Machining of Curved Long Flow Path Profiles. Proc. 16th Int. MTDR Conf, 1976,

283-289 28.

Hoare J.P. et al. An Investigation into the Differences Between NaCI and NaCI03 as Electrolytes in Electrochemical Machining. ibid, 1969, 116(2), 199-203

54.

Hoare J.P. et al. Anion Effects on the Dissolution o f Steel in ECM Binary Electrolyte System. ibid, 1973, 120(11),

55.

Hoare J.P. et al. ECM Behaviour of Steel in NaCI04 Electrolytes. Corrosion, 1971, 27(5), 2 1 1 - 2 1 5

56.

Hoare J.P. et al. Current Efficiency During the Electrochemical Machining of Iron and Nickel. Corrosion Science, 1975, 1 5 , 4 3 5 - 4 4 0

142-145

1452-1456

Sorkhel S. et al. Electrochemical Machining. Proc. 5th

AIMTDR Conf, 1972, 6 5 1 - 6 5 6 31.

53.

Capello G. et al. A New Approach to Electrochemical Finishing of Hardened Cylindrical Gear Tooth Face. Annals

CIRP, 1979, 28/1, 103-107 30.

Hoare J.P. Oxide Film Studies on Iron in Electrochemical Machining Electrolytes. J. Electrochem. Soc., 1970, 117(1),

Olenchich V.P. et al. Electrochemical Machining of Gears.

Machines and Tooling, 1975, 46(7), 5 9 - 6 0 29.

ASME, 1976, 98B, 360-368 52.

Sorkhel S.K. et al. Passivity Reduction in ECM Process.

7th ibid, 1976, 385-388 32.

Cole R.R. Basic Research in Electrochemical Machining -Present Status and Future Directions. Int.J. Prod. Res.,

1965, 4(2), 7 5 - 8 5 33.

35.

Heitmann H. Surface Roughness Formation and Determination of Optimum Working Conditions in Electrochemical Die Sinking. Proc. 2nd AIMTDR Conf, 1968, 4 1 - 4 7

58.

59.

Conf, 1976, 341-344 Ramani A.M. et al. Electrochemical Machining -- Developments in N.A.L. ibid, 3 5 7 - 3 6 0

37.

Kawafune K. Studies of Electrochemical Machining (1st Report). Bull. JSME, 1968, 11 (45), 554-564

Hoare J.P. et al. Electrochemical Machining. Research Publication GMR-2427, 11 May 1977, Research Laboratories Genera/Motors Corporation, Warren, Michigan, USA Hoare J.P. et al. Electrochemical Machining. Scientific

American, 1974, 230(I), 3 0 - 3 7 60.

Hoare J.P. et al. Some Effects of Solution Flow Rate on the Stability of Anodic Films on Steel. Corrosion, National

Association o{ Corrosion Engineers, USA, Paper No. 108, 1976, 1 - 7

Mahalingam V. et al. Measurement of Anodic Potential in ECM and its Relation to Surface Finish. Proc. 7th AIMTDR,

36.

Ito S. et al. Study of Electrolytic Machining Process(XI I).

J.Mech. Lab. Japan, 1967, 13(I), 15-24

Moir P.J. et al. Some Effects of Electrolyte Flow Rates on the Electrochemical Machining Process. Proc. 16th Int.

MTDR Conf, 1976, 275-281 34.

57.

61.

Evans J.M. et al. Effect of Surface Finish Produced During ECM Upon the Fatigue Life of Nimonic 80A. Proc. 12th

Int. MTDR Conf, 1972, 2 7 1 - 2 7 7 62.

Bhattacharyya A. et al. Analysis for Optimum Parametric Combination in Electrochemical Machining. Annals ClRP,

ibid, 565-569

63.

39.

Kawafune K. et al. Effect of the Working Parameters on the Working Gap During Electrochemical Machining. Annals CIRP, 1970, 1 8 , 3 0 5 - 3 1 7

Flueranbrack F. et al. Verification of a One-dimensional Two Phase Flow Model of the Frontal Gap in Electrochemical Machining. Trans. ASME, 1976, 98B(2), 4 3 1 - 4 3 7

64.

40.

Idem Study on Shape Processed by ECM. ibid, 1968,

Tipton H. The Determination of Tool Shape for Electrochemical Machining. Machinery and Production Engg., 14

38.

1973, 22/1, 5 9 - 6 0

Idem Studies of Electrochemical Machining (2nd Report).

Feb 1968, 3 2 5 - 3 2 8

345-352 65.

Idem The Dynamics of Electrochemical Machining. Proc.

41.

Nilson R.H. et al. Free Boundary Problem of E~M by Alternating Field Technique on Inverted Plane. Corn.Math.

42.

Mashinostroeniya V. Electrochemical Machining of Shaped Surfaces on Round Components by Tubular Shaped Cathode Tools. RE J, 1975, 55(11), 6 5 - 6 6

67.

Chin D.T. et al. Anodic Current Efficiency and Dimensional Control in Electrochemical Machining. J.Electrochem.Soc. 1973, 120(I 1), 1487-1493

68.

Cook N.H. et alo Experimental Studies in Electromachining. Trans. ASME, 1973, 95B, 9 4 5 - 9 5 0

69.

Fluerenbrack F. et al. Compressibility Effects in Electrochemical Machining. ibid, 1976, 98B(2), 4 2 3 - 4 3 0

70.

Tsuei Y.G. at al. Theoretical and Experimental Study of Workpiece Geometry in Electrochemical Machining. ibid,

71.

K~nig W. et al. The Influence of Process Parameters and Tool-Electrode Geometry on the Development of the Overcut in ECM with High Current Densities. Principles

4pp.Mech. and Engg, 1975, 6, 2 6 5 - 2 8 2

43.

44.

Chin D.T: Anodic Films and ECM Dimensional Control. A Study of Steel Electrodes in Solutions Containing Na 2 SO 4 and NaCI03 . ibid, 1972, 119(8), 1043--1048

45.

Idem Anodic Mechanism of Electrochemical Machining. Study of Current Transient on a Rotating Electrode. ibid,

46.

Datta M. et al. Film Breakdown on Nickel under Transpassive Dissolution Conditions in Sodium Nitrate Solutions.

5th Int. MTDR Conf, 1964, 509-522 66.

1977, 76-WA/Prod, 1 - 5

ibid, 1977, 124(4), 4 8 3 - 4 8 9

of Electrochemical Machining (Ed. C.L. Faust), 1971, 63-88. Published by Electrochemical Society Inc, PO Box 2071, Princeton N J, USA

Landolt D. et al. High Rate Anodic Dissolution o f Copper.

ibid, 1969, 116(10), 1384-1390 48.

Hoare J.P. Anodic Polarization Studies on Soft Iron in Some Concentrated Salt Solutions. Nature, 1968, 219(5158),

49.

Bannard J. The Effect of a Surface-Active Additive on the Dissolution Efficiency and Surface Finish in ECM. J.App.

72.

K~nig W. et al. Mathematical Model for the Calculation of the Contour of the Anode in Electrochemical Machining. Annals CIRP, 1977, 25, 8 3 - 8 7

73.

Larsson C.N. et al. Electrochemical Drilling: Tool Life and Feed Rate. Proc. lEE Conf. on Electrical Methods of

1034-1035

Electrochem, 1977, 7, 189--195 50.

Kops L. Effect of Pattern of Grain Boundary Network on Metal Removal-Rate in Electrochemical Machining. Annals

Machining, Forming and Coating, 1975, 161-168 74.

Loutrel S.P. et al. Sodium Chloride Electrolyte Data at High Temperatures and Pressures. Trans. ASME, 1973, 95B, 997-1002

75.

Idem A Theoretical Model for High Rate Electrochemical Machining. ibid, 1003-- 1008

ClRP, 1976, 25/1, 125-130 51.

Kops L. et al. Investigation on the Role of Workpiece Grain Size in Electrochemical Machining Process. Trans.

PRECISION ENGINEERING

Collett D.E. et al. A Complex Variable Approach to Electrochemical Machining Problems. J.Eng.Math, 1970,

4(I), 2 9 - 3 7

1971, 118(1), 174-179

47.

Segerlind L.J. Applied Finite Elements Analysis. John Wiley

and Sons Inc, New York, 1976

205

76.

Idem High Rate Electrochemical Machining. ibid, 992--996

77.

McGaough J.A. and Rasmuseen H. Theoretical Analysis of Electrochemical Machining Process. ibid, 69-WA/Prod-10

78.

Idem A Perturbation Analysis of the Electroforming Process. J.Mech. Engg. Sci, 1976, 18(6), 2 7 1 - 2 7 8

79.

chemical Machining Problems. J.Engg.Maths, 1970, 4(1), 29--37 104.

Hewson-Browne R.C. Further Applica'dons of Complex Variable Methods to Electrochemical Machining Problems. ibid, 1971, 5(4), 2 3 3 - 2 4 0

Idem On the Derivation of the Quasi Steady Model in Electrochemical Machining. J.Inst. Maths. Applcs, 1974, 13, 13--21

105.

Jain V.K. and Pandey P.C. On Some Aspects of Tooling Design in ECM. Mechanical Engg. Bulletin, 1977, 8(3), 66- 72

80.

Thorpe J.F. Operating Characteristics of Electrochemical Machining. Fertigung, 1975, 6(2), 4 7 - 5 0

106.

Zienkiewicz O.C. et ah Finite Elements in the Solution of Field Problems. The Engineer, Sept. 1965

81.

Thorpe J.F. et al. Analytical Determination of the Equilibrium Electrode Gap in Electrochemical Machining. Int.J. MTDR, 1969, 9, 131-144

107.

Desai C.S. et al. Introduction to the Finite Element Method. Affiliated East-West Press Pvt. Ltd., New Delhi, 1977

108.

82.

Tipton H. Metal Machining Improved by Electrical Methods. Electrical Review, 12 March 1965, 4 0 4 - 4 0 6

Jain V.K. and Pandey P.C. An Analysis of ECWCP Using Finite Elements Technique. Proc. 20th Int.MTDR Conf, 1979, 6 3 1 - 5 3 6

83.

Tipton H. The Determination of Shape of Tools for Use in Electrochemical Machining. Res. Report No. 40, 1971, MTIRA, UK

109.

Idem Application of Finite Elements Technique for the Analysis of Electrochemical Wire Cutting Process. (To be published in J.Institution o f Engrs (India))

84.

Idem The Calculation of Tool Shapes for Electrochemical Machining. Fundamentals o f Electrochemical Machining (Ed.C.L. Faust), 1971, 87-102. Published by Electrochemical Society Inc, Princeton N J, USA

110.

Idem Design of ECM Tooling Using Finite Elements Technique. Proc. 8th AIMTDR Conf., 1978, 5 6 6 - 5 7 0

111.

Idem Application of Finite Elements Techniques to ECM. Int.J.Engg.Prod. fin the press)

85.

Idem Tool Design for Electrochemical Machining. Machinery and Production Engineering, 1968, 869--871

112.

86.

Nilson R.H. et al. Inverted Cauchy Problems for the Laplace Equation in Engineering Design. J.Engg.Math. 1974, 8(4), 329-337

Idem An Analysis and Parametric Study of ECM Process. To be published in International Journal o f Engg. Production, New Delhi, 1979, 3(3)

113.

Idem Anode Shape Prediction in Deep Hole Drilling During ECM. Proc. 4th Biennial Seminar on Gas Turbine Organised by G TRE (India)

114.

Idem Finite Element Approach to the Two Dimensional Analysis of Electrochemical Machining. Precision Engineering, Jan 1980, 2(1), 2 3 - 2 8

87.

Larsson C.N. et al. Tool Damage by Sparking in ECM. Proc. 18th Int. MTDR Conf, 1977, 4 9 9 - 5 0 5

88.

Larsson C.N. at ah Electrochemical Effects on Shape Reproduction in ECM. 19th ibid, 1978, 5 3 3 - 5 4 0

89.

McGeough J.A. Principles of Electrochemical Machining. Chapman and Hall, London, 1974

115.

Arantsev M.A. et al. New Electrolyte Delivery Method in Electrochemical Machining. REJ, 1972, 52(17), 4 8 - 4 9

90.

Bhatia S.M. Effects of Electrolyte Conductivity on ECM Process. M. Tech. Dissertation, H T Bombay, India, 1971

116.

91.

Jain V.K. and Pandey P.C. Anode Profile in ECD. (To be published)

Dietz H. et al. Electrochemical Turning Considerations on Machining Rates which can be Attained. Annals CIRP, 1979, 28/I, 9 3 - 9 7

117.

92.

Idem Design and Analysis of ECM Process. Precision Engineering, October 1979, 1 (4), 199-206

Ippolito R. Effects of Fluidodynamic Conditions on the Gap Resistance in ECM. Annals CIRP, 1975, 24/I, 115-119

118.

Electrochemical Machining. PERA Report No. 145, 1965

93.

Chikamori K. Cell Voltage vs Current Density Curves During Electrolysis of NaCI Solutions. Denki Kagaku, 1973, 41(4), 275-280

119.

Nanayakkara M.B. et al. Computation and Verification of Workpiece Shape in Electrochemical Machining. 20th Int. MTDR Conf. 1979, 817-624

94.

Nilson R.H. et al. Free Boundary Problem for the Laplace Equation with Application to ECM Tool Design. Trans. ASME, 1976, 98(I), 5 4 - 5 8

120.

Rasch F.O. et al. Summary Report of Cooperate Research. Annals CIRP, 1978, 27/2, 561-564

95.

McCormick J.M. et al. Numerical Methods in Fortran. Prentice Hall o f India Pvt. Ltd., New Delhi, 1971

121.

Kuppuswamy G. et al. Magnetic Field on Electrolytic Grinding. Proc 8th AIMTDR Conf, 1978, 5 4 6 - 5 4 9

96.

Vitkovitch D. (Ed) Field Analysis, Experimental and Computational Methods. Van Nostrand, 1966

122.

Bhirud L.L. Matrix Operations on the Computers. Oxford and IBH Publishing Co., New Delhi, 1975

97.

Bewley L.V,,Two Dimensional Fields in Electrical Engg. Dover, 1963

123.

Bathe K.J. and Wilson E.L. Numerical Methods in Finite Element Analysis. Prentice Hall o f India Pvt. Ltd., New Delhi, 1978

98.

Martin H.C. at al. Introduction to Finite Element Analysis.

124.

Chin D.T. Dissolution of Mild Steel with a Rotating Ring Hemi-spherical Electrode. J. Electrochem. Soc, 1974, 121(12) 1592-1595

125.

Wilson J.F. Practice and Theory of Electrochemical Machining. Wiley Interscience, New. York, 1971

126.

Pandey P.C. and Shan H.S. Modern Machining Process. Tara McGraw Hill Publishing Co. Ltd., New Delhi, 1"980

127.

Jain M.K. An Analysis of ECM Process for Anode Shape Prediction. Ph.D. Thesis, University of Roorkee, Roorkee, India, 1980

128.

Landolt D. Throwing Power Measurements During High Rate Nickel Dissolution Under Active and Transpassive Conditions. J. Electrochem. Soc, 1972, 119(6), 708--712

Tara McGraw Hill Publishing Co. Ltd., New Delhi, 1975 99.

Dietz H. et al. Reproduction Accuracy with Electrochemical Machining: Determination of Side Gap. Annals CIRP, 1973, 22/1, 6 1 - 6 2

100.

Ippolito R. et al. The Workpiece Shape in the Electrochemical Machining of External Profiles Using Coated Electrodes. Int.J.MTDR, 1976, 16, 129-136

101.

Ganesh S. et al. Nomographic Approach for Planning in ECM. J. Inst.Eng. (India), 1979, 2 0 0 - 2 0 4

102.

Bhatia S.M. at al. Effect of Some Variables on ECM Accuracy. ibid, 1977, 58, 8 4 - 8 7

103.

Collett D.E. et al. A Complex Variable Approach to Electro-

206

PRECISION ENGINEERING