74 연구논문 유한요소시뮬레이션을통한저온분사코팅의계면접합에대한연구 배규열 * 강기철 * 윤상훈 * 이창희 * * 한양대학교신소재공학부 Finite Element Simulation of Interface Bonding in Kinetic Sprayed Coatings *Gyuyeol Bae, *Kicheol Kang, *Sanghoon Yoon and *Changhee Lee *Kinetic Spray Coating Laboratory (NRL), Division of Materials Science & Engineering, Hanyang University, Seoul 33-79, Korea *Corresponding author : chlee@hanyang.ac.kr (Received September 7, 28 ; Revised October 5, 28 ; Accepted October 22, 28) Abstract A finite element modeling approach has been described for the simulation and analysis of the micron-scaled solid particle impact behavior in kinetic spraying process, using an explicit code (ABAQUS 6.7-2). High-strain-rate plastic deformation and interface bonding features of the copper, nickel, aluminum, and titanium were investigated via FEM in conjunction with the Johnson-Cook plasticity model. Different aspects of adiabatic shear instabilities of the materials were characterized as a concept of thermal boost-up zone (TBZ), and also discussed based upon energy balance concept with respect to relative recovery energy (RRE) for the purpose of optimizing the bonding process. Key Words : Kinetic spraying, Finite element modeling (FEM), High-strain-rate impact behavior, Interface bonding, Coating. 서론 저온분사 (Kinetic spraying 또는 cold gas dynamic spraying) 는기존의열용사 (thermal spraying) 와는달리, 코팅대상면과의접촉전분말의온도가융점이하를유지하기때문에, 고온의대기환경에서발생할수있는산화및상변태의영향을최소화할수있는장점을지닌, 최근각광받고있는신코팅기술이다 ). 전형적인금속코팅이외에도벌크비정질 (bulk metallic glass) 코팅 2), 다양한복합재료코팅 3-5), 손상금속재의보수 (repair) 코팅 6) 등이본기술의다양한산업분야에의적용잠재가능성을보여주고있다. 수렴-발산 (convergent-divergent) de Laval 형태의노즐 (nozzle) 을통해형성된초음속가스 ( 주로질소또는헬륨 ) 는운동량전달을통해수십마이크론 ( 5 μm ) 크기의금속입자를초고속 ( 일반적으로 3-2 ms - ) 으로가속시킨다 7). 매우짧은시간 ( -8 - -7 s) 동안의충돌을통해야기 된입자와코팅대상재의고변형률소성변형 (high-strain-rate plastic deformation) 은충돌계면의국부적인단열온도상승 (localized adiabatic heating) 을유발하여고상접합 (solid- state bonding 또는 adhesion) 에의한재료의적층을가능하게한다 8). 이러한금속간접합은압연공정 (rolling process, Czichos 969) 을통해확인된저온용접 (cold welding) 에의한동종혹은이종금속간결합 (metallic bonding) 으로이해될수있다 9). 하지만, 수십나노초의극히짧은시간에발생하는입자충돌의고변형률소성변형현상에의한금속간접합과정은좀더복잡하며실험적접근및규명이어려워, 그접합기구가아직까지정확하게밝혀지지않았다. 이와관련하여, 재료고유의적층결함에너지 (stacking fault energy) 에따른충돌접합계면의동적회복 (dynamic recovery) 및동적재결정 (dynamic recrystallization) 현상이보고되었고 ), 최근알루미늄의동적비정질화 (dynamic amorphization) 와같은국부적상변태현상이보고되었다 ). 이러한재료의 68 Journal of KWJS, Vol. 26, No. 6, December, 28
유한요소시뮬레이션을통한저온분사코팅의계면접합에대한연구 75 미세구조변화는충돌시계면의단열온도상승에기인한단열소성불안정 (adiabatic shear instability) 현상 8) 과관련이있으며, 이에따른열적활성화반응 (thermally activated process) 이계면접합 (interface bonding) 을형성하는접합에너지 (adhesion energy) 로작용하여, 이에너지가반발에너지 (rebound energy) 보다클때, 성공적인접합이가능해진다 2). 한편, 이러한접합현상을현상학적으로규명하기위해, 여러연구자들이유한요소해석 (FEA) 방법을이용한컴퓨터시뮬레이션을수행한결과, 각재료의소성불안정현상이나타나는속도와실험적으로측정된접합임계속도 (critical velocity) 가유사함을밝혀냈고이를기준 (criterion) 으로하여, 본기술의공정최적화를위해, 다양한재료의임계속도를산출해냈다 8, 3-6). 특히, 이단열소성불안정에대해서는구리 (copper) 를통한단하나의정형화된현상으로소개되고설명되어왔으나, 이것이모든재료에적용되기는어렵다. 본질적으로이현상은각재료의물리적, 기계적특성과더불어고변형률변형조건하에서의가공경화 (work hardening) 및열적연화 (thermal softening) 특성에따라다르게나타나며, 이것이결국이소성불안정의개시 (onset) 에영향을미치게된다. 본연구에서는유한요소시뮬레이션을통해, 재료마다다르게나타나는단열소성불안정현상을근간으로계면접합의일반적인형태와원리를재료의물성및고변형률변형특성에따른열적상승지역 (thermal boost-up zone, TBZ) 으로특징지어설명하고, 에너지균형 (energy balance) 원리에입각하여정의된상대적회복에너지 (relative recovery energy, RRE) 와의상관관계를규명함으로써, 각재료의임계속도에따른접합형태및코팅특성을고찰하고예측하고자하였다. 2. 유한요소모델링 2. 유한요소해석방법 마이크론크기를갖는입자의고속충돌비선형과도 (nonlinear transient) 현상을모델링하기위해, 상업용유한요소해석 (FEA) 패키지인 ABAQUS 6.7-2를이용하였다. 상대적으로큰변형과짧은반응시간을갖는입자접합의동적과도현상을모사하기위해, 외연시간적분법 (explicit time integration) 알고리즘 7) 을적용하였으며, 축대칭 (axisymmetric) 모델을이용한완전히결합된열-응력 (fully coupled thermal- stress) 해석이수행되었다. Fig. 은입자충돌모델의기하학적형태와요 Fig. Geometry and mesh design of finite element model for particle/substrate impact simulation 소분할 (mesh) 을보여주고있다. 충돌계면의보다정확한해석을위해조밀한요소를사용하였고, 요소의극심한변형및일그러짐 (distortion) 으로인한수학적오차를최소화하기위해, Arbitrary Lagrangian Eulerian(ALE) 요소재분할 (remeshing) 방법 7) 을적용하였다. 2.2 소성변형해석모델 고속입자충돌에의한고변형률소성변형을해석하기위해, 아래와같이재료의다양한온도와변형률에따른변형경화 (strain hardening), 변형률경화 (strain-rate hardening), 열적연화현상을묘사할수있는 Johnson-Cook 소성모델 8) 을이용하였다. & ε p T T n σ = [ A + Bε ] + p C ln & ε Tm T ref ref m () 여기서, 는등가유동응력 (equivalent flow stress), 는등가소성변형, 는등가소성변형률, 는노멀라이징 ( 일반적으로 s - ) 된기준변형률 (normalizing reference strain rate) 을각각나타낸다. A, B, C, n, 그리고 m 값은재료고유의특성치로서, A는준안정 (quasi-static) 상태의단순인장또는압축실험을통해얻어진항복응력, B는가공경화계수, C는변형률가공경화계수, n은가공경화지수, m은열적연화지수를각각나타낸다. 또한, T m 은재료의융점이고, T ref 는기준온도로서, 일반적으로상온 (298 K) 의값을갖는다. 본모델에서는매우짧은입자충돌시간 ( 나노초 ) 으로인하여, 입자와모재 (substrate) 의변형시, 열의이동이없는단열 (adiabatic) 상태가가정되었다. 이것은 大韓熔接 接合學會誌第 26 卷第 6 號, 28 年 2 月 69
76 배규열 강기철 윤상훈 이창희 접촉시간 (contact time) 과재료의열확산율 (thermal diffusivity) 에의해결정된열확산거리 (thermal diffusion distance) 가시스템의크기 (characteristic system dimension) 에미치지못하기때문이다 8). 따라서, 열전도가고려되지않은소성변형에너지와내부에너지의균형은다음의식으로표현될수있다. ρ υ c dt = β σ dε p (2) 여기서, 는밀도, 는열용량을나타내며, (.9) 는 Taylor-Quinney 상수로서소성변형이열로전환된비율과같다. 2.3 입자충돌시뮬레이션 본연구에서는물성과열-기계적 (thermo-mechanical) 특성이다소다른 4가지재료 ( 구리, 니켈, 알루미늄, 티타늄 ) 를선택하여, 입자와모재의동종재료간고속충돌에의한고변형률변형현상을시뮬레이션하였다. 입자의충돌속도는 2- ms - 범위내에서 5 ms - 간격으로증가시켰으며, 입자의크기는 25 μm로고정하였다. 또한, 입자와모재의초기온도는상온 (298 K) 으로가정하였다. Table 의재료물성치와고변형률 Johnson-Cook 상수들은재료물성자료기지 (database) 와문헌을참고하였다 8-2). (d) Fig. 2 Simulated deformation features of a copper particle on a copper substrate, for the impact velocities of 4, 45, 5, and (d) 55 ms -. The arrows represent the velocity vectors of particle and substrate 변형이완료된시점에서의각충돌속도별속도벡터를통해확인할수있듯이, 상대적으로속도가낮은경우 (4 ms - ) 에는변형에의한에너지소모보다반발에너지의양이더큰것을예측할수있다. 충돌속도가증가함에따라, 입자와모재의변형이점차극심해지고, 변형에의한에너지소모가커짐으로인해반발되는정도가현저하게줄어듦을알수있다. 특히, 속도가 5 ms - 이상이되면충돌계면의분출 (jet) 현상이시작되어, 55 ms - 에서그현상이두드러지게관찰된다. 이두속도구간에서변형거동을좀더자세히살펴보면, Fig. 3에서볼수있듯이, 입자충돌계면의등가소성 3. 유한요소해석결과및고찰 Fig. 2는충돌속도에따른구리입자와모재의동종간변형거동해석결과를보여주고있다. 충돌에의한 Table Material properties used in FE models Parameter / Material Cu Al Ni Ti Density (Kg m -3 ) 896 27 889 45 Young s modulus (GPa) 24 68.9 27 6. Poison s ratio.34.33.3.34 Heat capacity (J Kg - K - ) 383 94 456 528 Melting temperature (K) 356 96 726 923 A (MPa) 9. 48.4 63. 86.57 B (MPa) 292. 345.5 648. 48.6 n.3.83.33.39 C m Reference temperature (K) Reference strain rate (s - ).25.9 298..895 293.6.44 298.94.655 298 Fig. 3 Deformation and equivalent plastic strain profiles at the side of copper particle after impact onto copper substrate, for the impact velocities of 5 and 55 ms - 62 Journal of KWJS, Vol. 26, No. 6, December, 28
유한요소시뮬레이션을통한저온분사코팅의계면접합에대한연구 77 변형 (equivalent plastic strain) 의최대크기가두속도구간에서현저하게증가하며, 최대변형값을갖는계면영역또한확장된다. ( 이것은 2차원축대칭해석결과를축을중심으로 8도회전하여얻은결과이다.) 이때, 계면의변형구배 (gradient) 에따른유동응력 (flow stress) 의전개양상을살펴보면, Fig. 4와같이충돌계면의위치에따라그양상이다소다르게나타남을알 Stress(MPa) Stress(MPa) 6 5 4 3 2..2.4.6.8..2.4 Equivalent plastic strain 6 5 4 3 2 수있다. 각위치에서의응력-변형곡선은입자충돌시접촉시간 (62.7 ns) 동안 ns 간격으로얻어진결과이다. 고변형률소성변형에기인한유동응력의변화는식 () 에의해변형및변형률에따른가공경화와단열온도상승에의한열적연화와의경쟁으로이해될수있다. 지역 a에서는가공경화우세에의한유동응력의증가가두드러지며, 이것은충돌후탄성회복에의한반발력으로작용하게된다. 지역 b에서는어느정도열적연화가작용하지만, 가공경화와의치열한경쟁으로인해유동응력의오르내림 (fluctuation) 현상으로이어진다. 이와반면, 계면분출현상이두드러지는지역 c에서는극심한변형에기인한단열온도상승이열적연화의우세를초래하여, 유동응력이붕괴 (collapse) 되며그값은거의 에가까워진다. 이러한충돌계면의응력붕괴현상은계면의온도상승과밀접한관련이있는데, Fig. 5를통해확인할수있듯이, 임계속도 (55 ms - ) 에서의급격한온도상승전개로설명될수있다. 이것은앞서소개한단열소성불안정이론으로뒷받침되며, 열적상승지역 (TBZ) 6) 으로정의되어재료에따른계면접합을특징짓고설명할수있는이론적도구가된다. 저온분사공정에서입자의초기운동에너지는대부분충돌후입자와모재의변형에너지로변환되며, 일부는점성유동과접촉면의마찰에의한에너지로소진된다. 한편, 초기운동에너지의일정량이회복가능한탄성변형에너지 (recoverable elastic strain energy) 로재료내에저장된다. 이러한에너지변환관계는열역학제 법칙에따른에너지보존원리에의해다음과같은식으로간단히표현될수있다. 2 3 4 5 6 Stress(MPa) 6 5 4 3 2 Equivalent plastic strain 4 2 8 6 4 55 m/s 5 m/s 45 m/s 4 m/s Fig. 4 5 5 2 Equivalent plastic strain Stress-equivalent plastic strain curves for center, intermediate, and edge region of impacting interface, during Cu particle impact onto Cu substrate (at υ = 55 ms -, for t = 62.7 ns) 2 3 4 5 6 Time( -9 sec) Fig. 5 Temporal development of temperature in the most severely deformed region designated as thermal boost-up zone (TBZ) 6) during impact for different velocities from 4 to 55 ms -, indicating adiabatic shear instability 大韓熔接 接合學會誌第 26 卷第 6 號, 28 年 2 月 62
78 배규열 강기철 윤상훈 이창희 E + K = EI + EV EF (3) 여기서, E K 는입자의초기운동에너지, E I 는변환된재료의내부에너지, E V 는재료의점성유동에의해소진된에너지, E F 는재료간접촉면의마찰에너지를각각나타낸다. 유한요소시뮬레이션을통해산출된 E V 와 E F 의양은 E I 양에비해상대적으로매우작다 ( 대부분 % 이하 ). 또한, 재료의내부에너지 (E I ) 는다음과같이탄성회복에너지 (E R ) 와소성변형에너지 (E P ) 의합으로표현될수있다. E I t c el = E + E = ( σ : ε& dv)dτ R P V t ( c pl + σ : ε& dv)dτ V (4) 여기서, σ c 는 Johnson-Cook 소성모델식에따른유동응력, 은탄성변형률, 는소성변형률, V는단위부피, t는단위시간을각각나타낸다. 유한요소해석에의한외연시간적분계산을통해, 응력-변형관계에따른탄성회복에너지 (E R) 와소성변형에너지 (E P) 의산출이가능하다. 에너지균형 (energy balance) 원리에입각할때, 이소성변형에너지와점성및마찰에너지는충돌계면의단열온도상승을유발하여, 접합에너지 2) 를형성하며, 나머지재료내에저장된탄성회복에너지가접합을방해하는반발에너지 2) 로작용하게된다. Fig. 6은구리입자와모재의동종간충돌시속도에따른 E R 과 E P 의변화를보여주고있다. E R 은속도에따라증가하는반면, E R/E P(relative recovery energy, RRE) 의값은감소하는것을알수있다. 이것은앞서설명했듯이, 충돌속도 ( 즉, 운동에너지 ) 가증가함에따 라, 소성변형에너지가상대적으로증가하므로, 충돌계면의단열온도상승에의한열적연화가촉진되어, 유동응력의감소가우세해지기때문이다. 또한, 앞서 Fig. 5에서보인바와같이단열소성불안정현상이나타나는임계속도 (55 ms - ) 에서그값이최소가됨을알수있다. 이와같이, 저온분사코팅의성공적인계면접합은각재료의임계속도이상에서열적상승지역의형성에의한상대회복에너지의감소로이해될수있다. 여기서, 재료의열적상승지역은입자와모재의고변형률변형시수반되는단열소성불안정현상과직접적인관련이있으며, 재료의변형능 (deformability) 및열- 기계적특성의영향을받는다. Fig. 7은서로다른 4가지재료의각각의산출된임계속도에서의변형형태및온도구배해석결과를보여주고있다. 앞서보인바와같이, 구리가상당한변형과함께가장낮은임계속도를나타내고, 니켈은상대적으로높은임계속도에서구리에상응하는변형형태를보인다. 이에반해, 알루미늄과티타늄은구리와니켈에비해월등히높은임계속도를갖는데, 이것은이두재료의낮은변형에기인한것이다. 여기서, 각재료의임계속도는충분히견고하게접합된상태를고려하여, 단열소성불안정개시후, 계면의최대상승온도가포화한계 (saturation limit) 에도달할때의속도값으로취하였다. 이때, 각재료는서로다른변형능과밀도및열용량을갖고, 또한고변형률변형시의열적연화및가공경화효과가각각다르기때문에, Fig. 8에서와같이각임계속도에서단열소성불안정의양상 (aspect) 이다르게나타 25 36 E R ( -9 J) 2 5 E R / E P E R 34 32 3 28 E R / E P ( -3 ) 5 2 3 4 5 Particle velocity (m/s) Fig. 6 Plots of ratio between recoverable strain energy and plastic dissipation energy (E R/E P, termed as relative recovery energy, i.e., RRE), in case of Cu particle impact onto Cu substrate (υ = -55 ms - ) 26 24 (d) Fig. 7 Deformation features and temperature contours of particle impact onto substrate: Cu particle/cu substrate (at υ = 55 ms -, t = 62.7 ns), Ni particle/ni substrate (at υ = 7 ms -, t = 48.8 ns), Al particle/al substrate (at υ = 775 ms -, t = 32.5 ns), and (d) Ti particle/ti substrate (at υ = 865 ms -, t = 23.3 ns) 622 Journal of KWJS, Vol. 26, No. 6, December, 28
유한요소시뮬레이션을통한저온분사코팅의계면접합에대한연구 79 8 6 4 4 2 2 8 6 TBZ( - 2 ) 8 6 8 6 4 RRE( -3 ) 4 4 2 2 3 4 5 2 2 9 8 7 Time( -9 sec) Ti Al Cu Ni Fig. 9 TBZ columns versus RRE plots for the cases of Ti/Ti, Al/Al, Cu/Cu, and Ni/Ni (particle/substrate) at critical velocities of 865 ms -, 775 ms -, 55 ms -, and 7 ms -, respectively 6 5 4 3 2 8 6 4 2 8 6 4 5 5 2 25 3 35 Time( -9 sec) Critical velocity of 25 μm particle [m/s] 9 8 7 6 5 4 3 2 Cu Ni Al Ti Fig. Critical velocities (columns) and maximum temperature (plots) at interface for the cases of Cu/Cu, Ni/Ni, Al/Al, and Ti/Ti (particle/substrate) 8 Max. temp. at interface(% of Tm) 2 5 5 2 25 Time( -9 sec) Fig. 8 Temporal evolution of temperature in TBZ: Ni particle/ni substrate (υ= 7 ms - ), Al particle/al substrate (υ = 775 ms - ), and Ti particle/ti substrate (υ = 865 ms - ) 나며, 이에따라산술적으로정의된열적상승지역의크기 ( 노멀라이징값 ) 6) 를계산할수있다. 그결과, Fig. 9에도시된바와같이각임계속도에서티타늄, 알루미늄, 구리, 니켈순으로열적상승지역의크기가증가하며, 흥미롭게도, 앞서논의한에너지균형원리에의해산출된상대회복에너지의값이이와반비례관계를가짐을알수있다. 한편, Fig. 에서와같이각임계속도에서계면의최대온도는구리의경우융점의 99.7% 에도달하고, 그다음으로니켈, 알루미늄, 티타늄순으로임계속도값에반비례하여그값이감소한다. 접촉온도 (contact temperature) 와더불어, 접합인자 (adhesion factor) 2) 인접촉면적 (contact area) 과접촉시간 (contact time) 또한 Fig. 에서와같이앞서나열한재료순으로동일한양상을보인다. 따라서, 재료의각임계속도에서나타나는다른양상의단열소성불안정현상은열적상승지역의크기뿐아니라, 에너지균형원리에따른상대회복에너지와더불어, 접촉면적과접촉시간과같은접합인자와상호연관성을가짐을알수있다. 4. 결론 본연구에서는유한요소시뮬레이션을통해, 물리 / 기계적특성이다른구리, 니켈, 알루미늄, 티타늄의입자- 모재간고속충돌에따른고변형률열-기계적거동을분석및고찰하였다. 고변형률변형시, 충돌계면의열적연화효과가가공경화효과보다우세할때나타 大韓熔接 接合學會誌第 26 卷第 6 號, 28 年 2 月 623
8 배규열 강기철 윤상훈 이창희 Contact area( -2 m 2 ) 25 2 5 Fig. 5 Cu Ni Al Ti 8 7 6 5 4 3 2 Contact time( -9 sec) Contact area (dark gray columns) and contact time (light gray columns) for the cases of Cu/Cu, Ni/Ni, Al/Al, and Ti/Ti (particle/substrate) at critical velocity 나는단열소성불안정현상은재료마다다른양상을보이고, 이를기준으로예측된각임계속도에서의열적상승지역의크기는에너지균형원리에따라산출된상대회복에너지와반비례관계를갖는다. 또한, 각재료별임계속도는구리 (55 ms - ), 니켈 (7 ms - ), 알루미늄 (775 ms - ), 티타늄 (865 ms - ) 순으로높게예측되고, 접촉온도와더불어접합인자로작용하는접촉면적과접촉시간은이와반비례하는상관관계를보인다. 따라서, 재료의단열소성불안정양상에따른열적상승지역이저온분사코팅의계면접합에영향을미치며, 결과적으로코팅의형태및특성을결정한다. 향후높은접합효율과접합강도및낮은기공도를갖는다양한재료의최적화된코팅공정을위하여, 본연구에서논의한접합인자들뿐아니라입자의크기및온도등을함께고려한범용화된모델의연구및개발이필요하다. 후 기 본연구는 26년도정부 ( 과학기술부 ) 의재원으로한국과학재단의지원을받아수행하였습니다 (No. 26-2289). 참고문헌. T.H. Van Steenkiste, J.R. Smith, R.E. Teets, J.J. Moleski, D.W. Gorkiewicz, R.P. Tison, D.R. Marantz, K.A. Kowalsky, W.L. Riggs, P.H. Zajchowski, B. Pilsner, R.C. McCune, K.J. Barnett, Surf Coat Technol, (999), 62-7 2. S.H. Yoon, H.J. Kim, C. Lee, Surf Coat Technol, 2(26), 622-629 3. S.H. Yoon, H.J. Kim, C. Lee, Surf Coat Technol, 2(27), 9524-9532 4. H.J. Kim, C. Lee, S.Y. Hwang, Mater Sci Eng A, 39(25), 243-248 5. S.M. Shin, Y. Xiong, Y.G. Ji, H.J. Kim, C. Lee, Appl Surf Sci, 254(28), 2269-2275 6. J.C. Lee, H.J. Kang, W.S. Chu, S.H. Ahn, CIRP Ann Manufact Technol, 56(27), 577-58 7. J.W. Wu, H.Y. Fang, S.H. Yoon, H.J. Kim, C. Lee, Appl Surf Sci, 252(25), 368-377 8. H. Assadi, F. Gärtner, T. Stoltenhoff, H. Kreye, Acta Mater, 5(23), 4379-4394 9. H. Czichos, J Phys D: Appl Phys, 5(972), 89-897. C. Borchers, F. Gärtner, T. Stoltenhoff, H. Assadi, H. Kreye, J Appl Phys, 93(23), 64-7. Y. Xiong, K. Kang, G. Bae, S.H. Yoon, C. Lee, Appl Phys Lett, 92(28), 94 2. J.W. Wu, H.Y. Fang, S.H. Yoon, H.J. Kim, C. Lee, Scripta Mater, 54(26), 665-669 3. M. Grujicic, J.R. Saylor, D.E. Beasley, W.S. DeRosset, D. Helfritch, Appl Surf Sci, 29(23), 2-227 4. M. Grujicic, C.L. Zhao, W.S. DeRosset, D. Helfritch, Mater Design, 25(24), 68-688 5. T. Schmidt, F. Gärtner, H. Assadi, H. Kreye, Acta Mater, 54(26), 729-742 6. G. Bae, Y. Xiong, S. Kumar, K. Kang, C. Lee, Acta Mater, 56(28), 4858-4868 7. ABAQUS TM 6.7-2 user manual, Hibbitt, Karlsson and Soerensen, Pawtucket, RI, USA, 27 8. G.R. Johnson, W.H. Cook, Proc of 7 th Int Symp Ball, Hague, The Netherlands (983), 54-547 9. N.K. Gupta, M.A. Iqbal, G.S. Sekhon, Int J Imp Eng, 32(26), 92-944 2. S.M. Lee, S.W. Seo, K.J. Park, O.K. Min, The Korean Society of Mechanical Eng A, 27(7)(23), 52-58 2. MataWeb, available online at www.matweb.com 624 Journal of KWJS, Vol. 26, No. 6, December, 28