이 기술 요약은 M.D Ibrahim 외 저자가 2023년 INTERNATIONAL JOURNAL OF INTEGRATED ENGINEERING에 게재한 “Parametric Study for Runner Modifications of Die Casted Part with Venting Systems” 논문을 기반으로 합니다. STI C&D의 기술 전문가가 분석하고 요약했습니다.

키워드

• 주요 키워드: 다이캐스팅 기공 결함 감소
• 보조 키워드: 고압 다이캐스팅(HPDC), 러너 게이트 시스템, 벤트 시스템, CFD 해석, 파라메트릭 연구, 수치 해석

Executive Summary

• 도전 과제: 고압 다이캐스팅(HPDC) 공정에서 발생하는 가스 혼입은 제품의 기계적 특성과 표면 품질을 저하시키는 주요 원인인 기공 결함을 유발합니다.
• 해결 방법: 본 연구에서는 검사 계측기 테스트 피스를 대상으로 파라메트릭 연구와 CFD 해석을 수행하여 러너 게이팅 시스템을 최적화하고 에어 벤트를 도입했습니다.
• 핵심 돌파구: ‘외향 곡률 러너(outward curvature runner)’라는 새로운 러너 설계와 에어 벤트를 적용하여 용탕의 속도 및 온도 분포를 크게 개선했으며, 이를 통해 가스 기공을 효과적으로 감소시켰습니다.
• 핵심 결론: CFD 해석에 기반한 러너 설계 및 벤팅 시스템 최적화는 다이캐스팅 부품의 기공 결함을 최소화하고 전반적인 품질을 향상시키는 핵심 전략입니다.

도전 과제: 왜 이 연구가 CFD 전문가에게 중요한가

고압 다이캐스팅(HPDC)은 복잡한 형상의 부품을 대량 생산하는 데 필수적인 공정이지만, 고질적인 결함 문제에 직면해 있습니다. 고속으로 용탕을 주입하는 과정에서 발생하는 역류는 용탕 내부에 공기 방울을 혼입시켜 가스 기공 결함의 원인이 됩니다. 이러한 기공은 제품의 인장 강도와 연성을 저하시켜 내구성에 치명적인 영향을 미칩니다.

또한, 용탕이 금형 내에서 불균일하게 응고하면서 발생하는 수축은 수축 기공을 유발합니다. 특히 러너 게이트와 같은 좁은 단면을 통과하며 발생하는 난류는 용탕의 충전 과정을 방해하고 결함을 심화시킵니다. 기존의 평면 게이트 시스템은 이러한 난류와 불균일한 충전을 유발하는 경향이 있어, 제품 품질을 유지하기 위한 새로운 게이팅 시스템 설계 최적화의 필요성이 꾸준히 제기되어 왔습니다.

접근 방식: 연구 방법론 분석

본 연구는 실제 산업 문제를 해결하기 위해 실험적 분석과 수치적 분석을 결합했습니다.

• 실험적 분석: 먼저, 실제 마그네슘 다이캐스트 샘플 금형을 CT(Computed Tomography) 스캔하여 내부 기공 결함의 발생 위치와 분포를 비파괴 방식으로 정밀하게 분석했습니다. CT 분석 결과, 제품의 가장 두꺼운 16mm 두께 부위에서 기공 발생률이 가장 높다는 것을 확인했으며, 이 지점을 결함 개선을 위한 ‘임계점(critical point)’으로 설정했습니다.
• 수치적 분석 (CFD): CT 분석으로 확인된 결함을 줄이기 위해 CFD 소프트웨어를 사용하여 금형 설계를 수정하고 시뮬레이션을 진행했습니다. 핵심적인 수정 사항은 다음과 같습니다.
  1. 외향 곡률 러너 (Outward Curvature Runner): 용탕이 캐비티로 더 원활하게 유입되도록 곡선 형태의 외향 러너를 새롭게 설계했습니다.
  2. 에어 벤트 (Air Vents) 도입: 캐비티 내부의 가스 방울 포집을 제거하고 배출을 촉진하기 위해 에어 벤트를 설치했습니다.

이 최적화된 금형 설계를 바탕으로, 용탕의 주입 속도(inlet velocity)와 작동 압력(operating condition)을 주요 변수로 설정하여 파라메트릭 연구를 수행했습니다. 이를 통해 최적의 공정 조건을 도출하고, 새로운 설계가 용탕 유동 및 결함 감소에 미치는 영향을 정량적으로 분석했습니다.

돌파구: 주요 연구 결과 및 데이터

결과 1: 주입 속도와 작동 압력이 결함에 미치는 영향

파라메트릭 분석 결과, 주입 속도와 작동 압력이 임계점의 압력 및 온도에 직접적인 영향을 미치는 것으로 나타났습니다.

• 압력 변화: Table 3에서 볼 수 있듯이, 주입 속도가 1.5 m/s에서 3.0 m/s로 증가함에 따라 임계점의 압력은 345,250 kPa에서 1,345,800 kPa (작동 압력 -5 Pa 기준)까지 크게 증가했습니다. Figure 4의 3D 그래프는 이러한 관계를 시각적으로 보여줍니다. 높은 압력은 캐비티 내에 갇힌 공기 방울을 오버플로우나 에어 벤트 쪽으로 밀어내는 데 더 효과적이므로 기공 감소에 유리합니다.
• 온도 변화: 주입 속도가 증가하면 마찰열 발생으로 인해 임계점의 온도도 함께 상승하는 경향을 보였습니다(Figure 5 참조). 하지만 주목할 점은, 100 Pa 및 3.0 m/s의 고압/고속 조건에서는 온도가 다른 조건에 비해 오히려 약간 낮게 나타났습니다. 이는 높은 압력을 유지하면서도 응고 시간을 단축하고 불필요한 가스 흡수를 방지할 수 있어 생산성 향상에 기여할 수 있음을 시사합니다.

결과 2: ‘외향 곡률 러너’ 설계의 탁월한 성능

새롭게 제안된 ‘외향 곡률 러너와 에어 벤트’ 설계는 기존 샘플 금형 대비 용탕 유동을 획기적으로 개선했습니다.

• 온도 분포 개선: Figure 6은 두 설계의 온도 분포를 비교합니다. 기존 샘플 금형(a)의 고온 영역은 722.7 K ~ 762.1 K 범위였으나, 제안된 설계(b)에서는 683.3 K ~ 722.7 K로 훨씬 낮고 균일한 온도 분포를 보였습니다. 이러한 균일한 온도는 불균일한 응고로 인한 수축 결함을 방지하고 더 나은 응고 품질을 촉진합니다.
• 유동 균일성 향상: Figure 7의 유선(streamline) 분포는 제안된 설계(b)가 기존 설계(a)에서 발생하던 와류(vortex) 형성을 제거하고, 훨씬 더 일관되고 균일한 속도 분포를 만들어내는 것을 명확히 보여줍니다. 이는 용탕이 캐비티를 완전히 채우는 것을 돕고, 가스 혼입 가능성을 근본적으로 줄여줍니다. 또한, 에어 벤트 상단에서 나타나는 높은 속도(붉은색 영역)는 내부 가스를 효과적으로 흡입하여 배출하고 있음을 증명합니다.

R&D 및 운영을 위한 실질적 시사점

• 공정 엔지니어: 본 연구는 용탕 주입 속도와 진공 압력(작동 조건)을 조절하는 것이 임계 부위의 결함을 제어하는 데 매우 중요함을 보여줍니다. Table 3의 데이터는 특정 공정 조건이 제품 품질에 미치는 영향을 예측하는 데 활용될 수 있습니다.
• 품질 관리팀: Figure 3과 Table 1의 CT 스캔 데이터는 제품의 가장 두꺼운 부분(16mm)이 기공 결함 발생 가능성이 가장 높다는 것을 명확히 보여주므로, 해당 부위에 대한 품질 검사를 강화하는 기준을 마련할 수 있습니다.
• 설계 엔지니어: 이 연구 결과는 러너의 형상(특히 외향 곡률)과 에어 벤트의 전략적 배치가 응고 중 결함 형성에 지대한 영향을 미친다는 것을 입증합니다. 이는 금형 설계 초기 단계에서부터 유동 균일성을 확보하고 가스 배출을 최적화하는 설계가 필수적임을 시사합니다.

논문 상세 정보

Parametric Study for Runner Modifications of Die Casted Part with Venting Systems

1. 개요:

• 제목: Parametric Study for Runner Modifications of Die Casted Part with Venting Systems
• 저자: M.D Ibrahim, J. Jendia, Y. Kashiwabara, Y. Sunami
• 발행 연도: 2023
• 게재 학술지: INTERNATIONAL JOURNAL OF INTEGRATED ENGINEERING
• 키워드: Porosity occurrence, parametric analysis, numerical analysis, design optimization

2. 초록:

고압 다이캐스팅(HPDC)은 복잡한 부품을 고압으로 용탕을 캐비티 내부에 주입하여 생성하는 공정입니다. 다이캐스팅의 실패는 제품의 기계적 특성, 표면 품질 및 수명을 감소시킬 수 있습니다. 본 논문에서는 검사 계측기인 테스트 피스의 다이캐스팅 공정에 대한 파라메트릭 연구 및 전산 유체 역학(CFD) 분석을 수행합니다. 다이캐스팅 공정에서 사용되는 매개변수는 캐비티 내부의 용탕 유동 품질에 영향을 미치기 때문에 중요합니다. 따라서 테스트 피스의 다이캐스팅 공정에서 최적의 매개변수 사용을 조사하기 위해 파라메트릭 연구를 수행합니다. 러너 게이팅 시스템 설계 또한 제품의 품질을 유지하는 중요한 기준 중 하나입니다. 본 논문은 또한 가스 기공을 줄이는 데 있어 러너 게이팅 시스템 설계 최적화의 효과를 조사했습니다. 샘플 금형의 CT 스캔은 가스 기공 발생과 CFD 결과 간의 관계를 비교하기 위해 포함되었습니다. 본 논문은 다이캐스팅 결함을 줄이는 데 있어 속도 및 온도 분포를 개선할 수 있는 에어 벤트를 갖춘 ‘외향 곡률 러너’라는 새로운 러너 설계를 제안했습니다. 또한, 에어 벤트는 부피를 확장하고 더 높은 흡입을 촉진하여 캐비티 내부의 가스 방울 포집을 제거하기 위해 설치되었습니다.

3. 서론:

금속 사출 성형(MIM)에 대한 수요는 수년간 증가해 왔으며, 이러한 성장은 이 제조 공정이 폐기물을 최소화하면서 복잡한 제품을 생산할 수 있는 능력 때문입니다. 기술 발전과 산업 혁명의 진전은 다이캐스팅 부품을 경제적이고 효과적으로 대량 생산하는 데 중요합니다. V-LINE®SYSTEM은 산업용 사출 성형 공정에 사용되는 시스템으로, 가소화 및 사출 실린더를 분리하여 정확한 양의 용탕을 채우고, 더 나은 유체 주입을 위해 열을 제공하고 온도를 유지합니다. 제조업체는 다이캐스팅 공정의 품질을 확인하여 주조 결함을 최소화하고 생산성을 향상시키는 것이 중요합니다.

4. 연구 요약:

연구 주제의 배경:

고압 다이캐스팅(HPDC)은 복잡한 부품을 제조하는 데 널리 사용되지만, 고속 주입으로 인한 공기 혼입으로 가스 기공과 같은 결함이 발생하여 제품 품질을 저하시키는 문제가 있습니다.

이전 연구 현황:

불량한 게이팅 시스템은 가스 혼입 및 수축 기공을 유발할 수 있으며, 특히 평면 게이팅 시스템은 난류 및 불균일한 충전을 유발하는 것으로 알려져 있습니다. 진공 다이캐스팅은 가스 혼입을 줄이는 데 이론적으로 도움이 될 수 있습니다.

연구 목적:

본 연구의 목적은 테스트 피스의 다이캐스팅 공정에서 최적의 매개변수를 조사하고, 가스 기공을 줄이기 위해 러너 게이팅 시스템 설계를 최적화하는 것입니다.

핵심 연구:

실제 다이캐스트 샘플의 CT 스캔을 통해 결함을 분석하고, 이를 바탕으로 CFD를 사용하여 파라메트릭 연구를 수행했습니다. 특히, ‘외향 곡률 러너’와 에어 벤트를 포함하는 새로운 설계를 제안하고, 이 설계가 용탕의 유동, 온도 분포 및 기공 감소에 미치는 영향을 수치적으로 평가했습니다.

5. 연구 방법론

연구 설계:

실험적 분석(CT 스캔)과 수치적 분석(CFD 시뮬레이션)을 결합한 연구 설계를 채택했습니다.

데이터 수집 및 분석 방법:

Nikon X-ray CT 장비(XT: H225 Series Interior)를 사용하여 마그네슘 다이캐스트 샘플의 기공 데이터를 수집했습니다. CFD 소프트웨어(Ansys)를 사용하여 파라메트릭 및 수치 해석을 수행했으며, 주입 속도와 작동 압력을 변수로 설정하여 임계점에서의 압력과 온도 데이터를 수집하고 분석했습니다.

연구 주제 및 범위:

연구 대상은 90° 각진 모서리를 가진 사다리 모양의 검사 계측기 테스트 피스(200mm x 100mm, 두께 4, 8, 12, 16mm)입니다. 연구 범위는 CT 분석에서 결함이 가장 많이 발견된 임계점(16mm 두께 부품의 8mm 깊이)에서의 기공 형성을 줄이고 캐비티 내 용탕 유동성을 개선하는 데 중점을 둡니다.

6. 주요 결과:

주요 결과:

• 주입 속도가 증가하면 임계점의 압력과 온도가 모두 증가하며, 높은 압력은 공기 방울 분산에 유리합니다.
• 새롭게 제안된 ‘외향 곡률 러너’와 에어 벤트 설계는 기존 설계에 비해 캐비티 내 용탕의 온도 분포를 더 낮고 균일하게 만들었습니다.
• 제안된 설계는 와류 형성을 제거하고 유속을 균일하게 분배하여 용탕 충전 과정을 개선하고 가스 혼입을 줄였습니다.
• 에어 벤트는 내부 압력을 낮춰 흡입력을 높임으로써 가스 배출을 촉진하는 효과가 있었습니다.

그림 목록:

• Fig. 1 – Test piece with air vents
• Fig. 2 – Nikon x-ray CT machine XT: H225 Series interior
• Fig. 3 – Sample mold CT scan
• Fig. 4 – 3D graph of pressure at a critical point
• Fig. 5 – 3D graph of temperature at a critical point
• Fig. 6 – Temperature contour for (a) sample mold; (b) outward curvature runner with air vent
• Fig. 7 – Velocity streamline distribution for (a) sample mold; (b) outward curvature runner with air vent

7. 결론:

파라메트릭 연구 결과, 주입 속도가 증가하면 임계점의 압력도 증가하여 용탕 내 공기 방울을 분산시키는 데 도움이 될 수 있습니다. 임계점의 압력이 높을수록 마찰열이 더 많이 발생하여 용탕의 온도가 높아집니다. 임계점의 고온은 응고 속도를 늦추지만, 100 Pa 및 3.0 m/s에서는 고압임에도 불구하고 온도가 약간 낮아 제조 시간을 줄이는 데 도움이 될 수 있습니다. 특정 조건에서 진공 보조 다이캐스팅 공정은 가장 두꺼운 지점의 압력을 감소시킬 수 있는데, 이는 높은 진공이 용탕의 난류를 유발하여 역류 및 가스 혼입을 초래할 수 있기 때문입니다. 수치 해석에서 용탕의 높은 국부 온도는 균열, 번 마크, 유동 마크를 유발합니다. 높은 온도 차이는 또한 응고 속도의 불일치를 초래하여 금형 캐비티를 완전히 채우지 못하고 일부 용탕이 굳기 시작하는 불완전 충전을 유발할 수 있습니다. 불완전 충전은 가스 기공 결함으로 이어지는 가스 방울 혼입을 유발합니다. 와류 형성은 또한 용탕 속도를 감소시킵니다. 외향 곡률 러너는 러너가 바깥쪽을 향하는 다른 접근 방식을 보여주며, 러너 게이트 입구에서 생성된 곡선이 특징입니다. 이 러너는 용탕 속도와 온도 분포를 효과적으로 증가시킵니다. 또한, 에어 벤트는 내부의 낮은 압력으로 인해 더 높은 흡입을 촉진합니다.

8. 참고 문헌:

• [1] Grand View Research (2015). Metal Injection Molding Market Size, Share & Trends Analysis Report by End-Use (Automotive, Consumer Product, Medical, Industrial, Defense), By Region (North America, Europe, Asia Pacific, CSA, MEA), Competitive Landscape, And Segment Forecasts, 2018 -2025. https://www.grandviewresearch.com/industry-analysis/metal-injection-molding-mim-market
• [2] Luo, A. A. (2013). Magnesium casting technology for structural applications. Journal of Magnesium and Alloys, 1(1), 2-22. https://doi.org/10.1016/j.jma.2013.02.002
• [3] El-Sayed, M. A., Hassanin, H. & Essa, K. (2016). Bifilm defects and porosity in Al cast alloys. International Journal Advanced Technology, 86(1), 1173-1179. https://doi.org/10.1007/s00170-015-8240-6
• [4] Sadeghi, M. & Mahmoudi, J. (2012). Experimental and Theoretical Studies on the Effect of Die Temperature on the Quality of the Products in High-Pressure Die-Casting Process. Journal of Advances in Materials Science and Engineering, 2012(1), 1-10. https://doi.org/10.1155/2012/434605
• [5] Lee, B. D., Baek, U. H. & Han, J. W. (2012). Optimization of gating system design for die casting of thin magnesium alloys-based multi-cavity LCD housings. Journal of Materials Engineering and Performance, 21(9), 1893-1903, https://doi.org/10.1007/s11665-011-0111-1
• [6] Piyanut, M. & Chakrit, S. (2017). Evaluate of chill vent performance for high pressure die-casting production and simulation of motorcycle fuel caps. International Conference on Mechatronics and Mechanical Engineering (ICMME 2016), 95, 07025. https://doi.org/10.1051/matecconf/20179507025
• [7] Ibrahim, M.D., Rahman, M.R.A., Khan, A.A., Mohamad, M.R., Suffian, M.S.Z.M., Yunos, Y.S., Wong, L.K., & Mohtar, M.Z. (2017). Effect of Mold Designs on Molten Metal Behaviour In High-Pressure Die Casting. Conference Series: Journal of Physics, 822(1) 2017, 1-6. https://doi.org/10.1088/1742-6596/822/1/012029
• [8] Ramnath, B. V., Elanchezhian, C., Chandrasekhar, V., Kumar, A. A., Asif, S. M., Mohamed, G. R., Raj, D. V., & Kumar, C. S. (2014). Analysis and Optimization of Gating System for Commutator End Bracket. Procedia Materials Science, 6(1), 1312-1328. https://doi.org/10.1016/j.mspro.2014.07.110
• [9] Westengen, H., & Rashed, H. (2016). Magnesium Alloys: Properties and Applications. Encyclopedia of Materials: Metals and Alloys, 2016(1), 11-19. https://doi.org/10.1016/B978-0-12-803581-8.02568-6
• [10] Sodick Co. Ltd. (n.d). Sodick News Letter: Core technology. http://www.sodick.jp/tech/v_line.html
• [11] Ibrahim. M.D, Tan, G. L., Roslan, L., Kashiwabara, Y., Jendia, J., & Sunami, Y. (2020). Numerical and Experimental Analysis on Runner and Gate Positioning for Magnesium Alloy Die-Casted Test Piece. Materials Science Forum, 975(1), 242-247. https://doi.org/10.4028/www.scientific.net/MSF.975.242
• [12] Rai, J. K., Lajimi, A. M. & Xirouchakis, P. (2008). An intelligent system for predicting HPDC process variable in interactive environment. Journal of Materials Processing Technology, 203(1-3), 72-79. https://doi.org/10.1016/j.jmatprotec.2007.10.011
• [13] El-Sayed, M. A., Hassanin, H. & Essa, K. (2016). Bifilm defects and porosity in Al-cast alloys. International Journal Advanced Manufacturing Technology, 86(1), 1173-1179. https://doi.org/10.1007/s00170-015-8240-6

전문가 Q&A: 자주 묻는 질문

Q1: 왜 파라메트릭 연구의 변수로 주입 속도와 작동 조건을 선택했습니까?

A1: 논문에 따르면, 이 두 매개변수는 캐비티 내부의 용탕 유동 품질에 직접적인 영향을 미치는 핵심 공정 변수입니다. 주입 속도는 유동 패턴과 압력 형성을 결정하고, 작동 조건(진공 압력)은 캐비티 내 가스 배출 효율을 좌우합니다. 따라서 이 두 변수를 분석하는 것은 가스 혼입과 응고 과정을 제어하여 기공 결함을 줄이는 최적의 조건을 찾는 데 필수적입니다.

Q2: 논문에서 제안한 ‘외향 곡률 러너’는 기존 설계의 어떤 특정 문제를 해결합니까?

A2: 기존 설계, 특히 평면 게이트 시스템은 용탕이 캐비티로 진입할 때 난류와 와류(vortex)를 유발하기 쉽습니다. Figure 7의 유선 분석에서 볼 수 있듯이, ‘외향 곡률 러너’는 용탕이 캐비티로 부드럽고 균일하게 흐르도록 유도하여 이러한 와류 형성을 제거합니다. 이를 통해 불균일한 충전과 그로 인한 가스 포집 문제를 근본적으로 해결합니다.

Q3: 에어 벤트는 단순히 공기 배출구 역할을 넘어 기공 감소에 어떻게 기여합니까?

A3: 논문에서는 에어 벤트가 단순히 공기를 배출하는 수동적인 역할을 넘어, 부피를 확장하고 벤트 내부의 낮은 압력으로 인해 더 높은 ‘흡입(suction)’을 촉진한다고 설명합니다. 이는 캐비티 내부에 갇힌 가스 방울을 수동적으로 배출하는 것이 아니라, 능동적으로 빨아들여 제거하는 효과를 가집니다. 이로 인해 가스 제거 효율이 크게 향상됩니다.

Q4: Table 1을 보면 16mm 두께 섹션에서 결함 비율(0.0239%)이 유독 높습니다. 가장 두꺼운 부분이 가장 문제가 되는 이유는 무엇입니까?

A4: 논문에 따르면, 고온 영역은 주변부보다 늦게 응고됩니다. 가장 두꺼운 부분은 열을 가장 오래 유지하기 때문에 응고가 지연되고, 이 과정에서 주변부가 먼저 굳으면서 빠져나가지 못한 공기 방울이 갇히게 됩니다. 또한, 이러한 불균일한 응고는 불완전 충전의 원인이 되어 결함 발생 가능성을 더욱 높입니다.

Q5: Figure 5에서 100 Pa, 3.0 m/s 조건의 온도가 다른 고속 조건보다 약간 낮은데, 이것이 의미하는 바는 무엇입니까?

A5: 논문에서는 이 점을 긍정적으로 평가합니다. 높은 압력을 통해 가스 배출 효과를 극대화하면서도, 온도가 과도하게 상승하는 것을 막을 수 있기 때문입니다. 이는 불필요한 가스 흡수를 방지하고 응고 속도를 적절히 제어하여 전체 제조 시간을 단축하는 데 도움이 될 수 있음을 시사합니다. 즉, 생산성과 품질을 동시에 잡을 수 있는 최적의 공정 조건일 가능성을 보여줍니다.

결론: 더 높은 품질과 생산성을 향한 길

본 연구는 고압 다이캐스팅 공정의 고질적인 문제인 기공 결함을 해결하기 위해 CFD 해석을 활용한 체계적인 접근법을 제시합니다. 러너 형상을 ‘외향 곡률’로 최적화하고 에어 벤트를 전략적으로 배치함으로써 용탕의 유동을 안정시키고 가스 배출을 극대화할 수 있음을 입증했습니다. 이는 최종적으로 다이캐스팅 기공 결함 감소로 이어져 제품의 기계적 신뢰성과 품질을 한 단계 끌어올릴 수 있습니다.

STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 지원하는 데 전념하고 있습니다. 본 논문에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 부품에 어떻게 구현할 수 있는지 논의해 보십시오.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

저작권 정보

• 이 콘텐츠는 “M.D Ibrahim” 외 저자의 논문 “[Parametric Study for Runner Modifications of Die Casted Part with Venting Systems]”을 기반으로 한 요약 및 분석 자료입니다.
• 출처: https://doi.org/10.30880/ijie.2023.15.05.010

이 자료는 정보 제공 목적으로만 사용됩니다. 무단 상업적 사용을 금지합니다. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 Grzegorz Banaszek 외 저자가 Materials (2020)에 발표한 논문 “[Analysis of the Open Die Forging Process of the AZ91 Magnesium Alloy]”를 기반으로 합니다. STI C&D에서 기술 전문가를 위해 분석 및 요약했습니다.

Keywords

• Primary Keyword: AZ91 마그네슘 합금 단조
• Secondary Keywords: 개방형 단조, 유한요소해석, 형상 앤빌, 공정 최적화, 수치 모델링, 평면 앤빌

Executive Summary

• The Challenge: AZ91 마그네슘 합금은 경량화에 유리하지만 성형성이 낮아 고품질 단조품 생산에 어려움이 있습니다.
• The Method: 유한요소해석(FEM) 소프트웨어 Forge®NxT를 사용하여 평면 앤빌과 마름모-사다리꼴(rhombic-trapezoid) 형상 앤빌을 이용한 AZ91 합금의 개방형 단조 공정을 수치적으로 모델링하고 비교 분석했습니다.
• The Key Breakthrough: 형상 앤빌을 사용하면 단조 패스 수를 6회에서 4회로 줄이면서도 전체 부피에 걸쳐 더 균일한 온도 및 변형률 분포를 얻을 수 있음을 확인했습니다.
• The Bottom Line: 공구 형상 최적화를 통해 AZ91 합금 단조 공정의 효율성과 최종 제품의 품질 균일성을 크게 향상시킬 수 있으며, 이는 생산 비용 절감으로 직결됩니다.

The Challenge: Why This Research Matters for CFD Professionals

자동차, 항공우주, 전자 산업에서 경량 구조 재료에 대한 수요가 끊임없이 증가하면서 마그네슘 합금이 주목받고 있습니다. 특히 AZ91 합금은 높은 비강도(밀도 대비 강도) 덕분에 매력적인 소재입니다. 하지만 대부분의 마그네슘 합금 반제품은 주조로 생산되며, 단조와 같은 소성 가공은 제한적으로 이루어집니다.

이는 마그네슘 합금이 육방정계(HCP) 결정 구조를 가져 상온에서 소성 변형이 어렵고, 가공 중 기계적 이방성이 크게 발생하기 때문입니다. 고품질의 단조품을 생산하기 위해서는 온도, 변형률, 공구 형상 등 수많은 공정 변수를 정밀하게 제어해야 합니다. 기존의 압출 공정으로는 자유로운 형상 제어가 어렵기 때문에, 개방형 단조 기술은 고신뢰성이 요구되는 구조 부품 생산에 필수적입니다. 본 연구는 이러한 기술적 난제를 해결하고, AZ91 마그네슘 합금 단조 공정의 최적화 방안을 제시하기 위해 수행되었습니다.

The Approach: Unpacking the Methodology

본 연구는 AZ91 마그네슘 합금의 개방형 단조 공정을 수치적으로 분석하기 위해 상용 유한요소해석(FEM) 소프트웨어인 Forge®NxT를 활용했습니다. 연구의 신뢰성을 확보하기 위해 다음과 같은 체계적인 접근법을 사용했습니다.

• 소재 물성 확보: Gleeble 3800 금속 시뮬레이션 시스템을 이용해 실제 단조 조건(온도 200-400°C, 변형률 속도 0.1-10 s⁻¹)에 해당하는 AZ91 합금의 유동 응력(flow stress) 데이터를 확보했습니다. 이 데이터는 Hansel-Spittel 방정식을 통해 근사화되어 시뮬레이션의 재료 모델로 사용되었습니다.
• 단조 공구 설계: 기존 연구 및 문헌 조사를 바탕으로 두 가지 유형의 단조 공구(앤빌)를 설계했습니다. 하나는 전통적인 평면 앤빌(flat anvils)이고, 다른 하나는 본 연구에서 제안하는 마름모-사다리꼴 형상 앤빌(rhombic-trapezoid anvils)입니다.
• 시뮬레이션 조건:
  • 초기 소재: 80mm x 80mm x 80mm 크기의 사각 단면 빌렛
  • 초기 온도: 소재 400°C, 금형(앤빌) 25°C
  • 단조 속도: 상부 앤빌 속도 8 mm/s
  • 마찰 계수: 쿨롱 마찰 모델, μ = 0.3
• 공정 시나리오: 평면 앤빌은 6회의 단조 패스(3회 후 90° 회전, 다시 3회)를, 형상 앤빌은 4회의 단조 패스(2회 후 90° 회전, 다시 2회)를 적용하여 각 공정에서의 온도 및 변형률 분포를 비교 분석했습니다.

The Breakthrough: Key Findings & Data

시뮬레이션 분석 결과, 앤빌의 형상이 단조 공정의 효율성과 최종 제품의 품질에 결정적인 영향을 미친다는 사실이 명확하게 드러났습니다.

Finding 1: 형상 앤빌 사용 시 공정 효율성 및 균일성 대폭 향상

마름모-사다리꼴 형상 앤빌을 사용했을 때, 단조 공정이 더 효율적이고 안정적으로 진행되었습니다. 가장 큰 장점은 총 단조 패스 수를 평면 앤빌의 6회에서 4회로 줄일 수 있다는 점입니다. 이는 공정 시간과 에너지 소비를 직접적으로 줄여 생산 비용 절감에 기여합니다.

또한, 형상 앤빌은 소재의 흐름을 제한하여 단조품 전체에 걸쳐 매우 균일한 변형률 분포를 만들어냈습니다. Figure 16에서 볼 수 있듯이, 최종 4차 패스 후에도 단면 전체의 변형률 강도가 고르게 분포되어 있습니다. 이는 최종 제품의 기계적 특성이 위치에 따라 균일하게 나타날 것임을 시사하며, 고품질 부품 생산에 매우 유리한 조건입니다.

Finding 2: 평면 앤빌은 높은 변형률 강도 달성에 유리

반면, 평면 앤빌을 사용한 경우 더 높은 변형률 강도 값을 얻을 수 있었습니다. Figure 15d를 보면, 최종 6차 패스 후 단조품의 중심부와 표면부에서 변형률 강도가 1.27에서 1.48 범위에 도달했습니다. 이는 형상 앤빌의 결과보다 높은 수치입니다. 높은 변형률은 소재의 가공 경화를 유도하고 미세조직을 개선하여 기계적 특성을 향상시키는 데 긍정적인 역할을 합니다.

하지만 평면 앤빌은 공정 초기에 변형이 불균일하게 발생하며, 이를 해소하기 위해 더 많은 단조 패스가 필요했습니다. 그럼에도 불구하고, 6회의 패스를 거친 후에는 중심부의 변형률 분포가 만족스러운 수준으로 균일해졌습니다.

Practical Implications for R&D and Operations

본 연구 결과는 AZ91 마그네슘 합금 단조 공정을 설계하고 운영하는 데 있어 중요한 실무적 지침을 제공합니다.

• For Process Engineers: 이 연구는 초기 단조 단계에서는 균일한 변형을 유도하는 형상 앤빌을 사용하고, 최종 단계에서는 높은 변형률을 부여하여 기계적 특성을 극대화하는 평면 앤빌을 사용하는 하이브리드 공정을 고려해볼 수 있음을 시사합니다. 또한, 40%의 단조 압하율과 0.8의 상대 이송률(relative feed)이 효과적인 공정 변수임을 제안합니다.
• For Quality Control Teams: Figure 15와 Figure 16의 변형률 분포 데이터는 단조품의 위치별 기계적 특성 편차를 예측하는 데 활용될 수 있습니다. 이를 통해 품질 검사 기준을 수립하거나 취약부를 예측하여 검사를 강화하는 전략을 세울 수 있습니다.
• For Design Engineers: 본 연구 결과는 앤빌의 형상이 소재의 유동, 온도 분포, 최종 변형률에 미치는 영향을 명확히 보여줍니다. 이는 고품질 단조품 생산을 위한 최적의 금형을 설계하는 데 있어 핵심적인 데이터를 제공합니다.

Paper Details

Analysis of the Open Die Forging Process of the AZ91 Magnesium Alloy

1. Overview:

• Title: Analysis of the Open Die Forging Process of the AZ91 Magnesium Alloy
• Author: Grzegorz Banaszek, Teresa Bajor, Anna Kawałek and Tomasz Garstka
• Year of publication: 2020
• Journal/academic society of publication: Materials
• Keywords: magnesium alloy AZ91; physical modelling; open die forging; flat anvils; shaped anvils

2. Abstract:

The paper presents the results of numerical modelling of the forging process of magnesium alloy ingots on a hydraulic press with the use of flat and shaped anvils. The use of shaped (rhombic-trapezoid) anvils will affect the uniform distribution of temperature and strain intensity in the entire volume of the forging, causing a number of forging passes, which in consequence will reduce the costs of the blank manufacturing process. However, higher values of the strain intensity were obtained during the deformation of the material in flat anvils. The purpose of the research was to propose assumptions for forging technology of producing a blank from AZ91 alloy with the use of flat and shaped anvils. Numerical examination for AZ91 magnesium alloy was carried out using the Forge®NxT commercial software. The rheological properties of the investigated alloy were determined on the basis of uniaxial compression tests carried out in the Gleeble 3800 metallurgical simulation system. The numerical analysis of the process of forging AZ91 alloy ingots on a press was conducted in the temperature range of 200-400 °C and at several forging passes.

3. Introduction:

최근 경량 구조용 재료에 대한 지속적인 탐색으로 인해 마그네슘 합금에 대한 관심이 증가하고 있습니다. 이 합금들은 자동차, 조선, 항공, 우주 및 전자 산업 등 여러 분야에서 매력적인 특성을 가지고 있습니다. 밀도 대비 가장 높은 강도를 보여 경량 구조물 생산에 유리하게 작용합니다. 대부분의 마그네슘 합금 반제품은 주조법으로 생산되며, 압출, 압연, 단조 공정은 상대적으로 적게 사용됩니다. 신뢰성, 기능성, 기계적 특성 때문에 단조된 마그네슘 합금 부품은 엄격한 요구사항을 가진 구조 요소로 사용됩니다. 특히 개방형 단조 공정은 금형 작업 표면의 형상과 치수 매개변수 및 연신 작업의 주요 매개변수를 제어하여 금속 유동 운동학을 자유롭게 형성할 수 있게 해줍니다. 본 연구는 단조 기술 엔지니어링 분야의 연구를 통해 제안된 연구 문제에 대한 포괄적인 접근을 시도합니다.

4. Summary of the study:

Background of the research topic:

마그네슘 합금은 경량화 요구에 부응하는 핵심 소재이지만, 육방정계 결정 구조로 인해 성형성이 낮고 가공이 어렵다는 한계가 있습니다. 특히 고강도 AZ91 합금의 단조 공정은 균일한 미세구조와 기계적 특성을 확보하는 것이 중요합니다.

Status of previous research:

기존에는 마그네슘 합금의 주조 및 압출에 대한 연구가 주를 이루었으며, 단조 공정에 대한 연구, 특히 개방형 단조에서 공구 형상이 미치는 영향에 대한 체계적인 분석은 부족했습니다. 수치 해석 기법을 활용하여 공정을 최적화하려는 시도가 있었으나, 평면 앤빌과 특정 형상 앤빌을 직접적으로 비교 분석한 연구는 드물었습니다.

Purpose of the study:

본 연구의 목적은 평면 앤빌과 마름모-사다리꼴 형상 앤빌을 사용하여 AZ91 마그네슘 합금 빌렛을 생산하는 단조 기술의 가정을 제안하는 것입니다. 수치 모델링을 통해 각 앤빌 유형이 단조 공정 중 온도 분포와 변형률 강도에 미치는 영향을 분석하고, 더 효율적이고 균일한 품질의 제품을 생산할 수 있는 공정 조건을 찾는 것을 목표로 합니다.

Core study:

상용 FEM 소프트웨어 Forge®NxT를 사용하여 AZ91 합금의 개방형 단조 공정을 시뮬레이션했습니다. Gleeble 3800으로 측정한 소재의 유변학적 특성을 모델에 적용하고, 평면 앤빌과 형상 앤빌을 사용했을 때의 공정을 비교했습니다. 각 시나리오별로 여러 단조 패스에 걸쳐 단조품 내부의 온도 분포와 변형률 강도 분포를 추적하고 분석하여 최적의 공구 형상과 공정 변수를 도출하고자 했습니다.

5. Research Methodology

Research Design:

본 연구는 비교 실험 설계를 기반으로 합니다. 두 가지 다른 유형의 단조 공구(평면 앤빌, 형상 앤빌)를 독립 변수로 설정하고, 이에 따른 종속 변수(온도 분포, 변형률 강도 분포, 단조 패스 수)의 변화를 수치 시뮬레이션을 통해 분석했습니다.

Data Collection and Analysis Methods:

• 데이터 수집: Gleeble 3800 시스템을 이용한 단축 압축 시험을 통해 AZ91 합금의 응력-변형률 데이터를 수집했습니다. 수치 시뮬레이션은 Forge®NxT 소프트웨어를 통해 수행되었으며, 각 단조 패스 후의 온도 및 변형률 필드 데이터를 수집했습니다.
• 데이터 분석: 수집된 데이터는 Hansel-Spittel 방정식을 사용하여 재료 모델 계수를 결정하는 데 사용되었습니다. 시뮬레이션 결과는 단면 및 종단면의 등고선도(contour plot) 형태로 시각화하여 정성적 및 정량적으로 비교 분석되었습니다.

Research Topics and Scope:

• 연구 주제: AZ91 마그네슘 합금의 개방형 단조 공정에서 앤빌 형상이 공정 효율성 및 품질 균일성에 미치는 영향 분석.
• 연구 범위:
  • 소재: AZ91 마그네슘 합금
  • 공정: 유압 프레스를 이용한 개방형 단조
  • 공구: 평면 앤빌 및 마름모-사다리꼴 형상 앤빌
  • 분석 변수: 온도, 변형률 강도
  • 온도 범위: 200-400 °C

6. Key Results:

Key Results:

• 마름모-사다리꼴 형상 앤빌을 사용하면 단조품 전체 부피에 걸쳐 더 균일한 온도 및 변형률 강도 분포를 얻을 수 있습니다.
• 형상 앤빌을 사용함으로써 총 단조 패스 수를 평면 앤빌의 6회에서 4회로 줄일 수 있어, 생산 비용 절감이 가능합니다.
• 평면 앤빌을 사용했을 때 더 높은 변형률 강도 값을 얻을 수 있었으며, 이는 소재 가공 및 기계적 특성 향상에 유리할 수 있습니다.
• 두 경우 모두 금형과 소재의 접촉면에서 상당한 온도 강하가 발생했지만, 공정 중 소성 변형으로 인한 열 발생 덕분에 소재 중심부의 온도는 단조에 적합한 수준으로 유지되었습니다.
• 초기 단조 단계에서는 균일한 변형을 위해 형상 앤빌을, 최종 단조 단계에서는 높은 변형률 부여를 위해 평면 앤빌을 사용하는 것이 합리적인 전략일 수 있습니다.

Figure List:

• Figure 1. Geometry and dimension of load (a); and thermal cycle of physical simulation (b).
• Figure 2. Work-hardening curves for the strain rate range of (0.1 s−1−10.0 s−1) at a temperature of: (a) 200 °C, (b) 300 °C; (c) 400 °C. (Black indicates the experimental curves; red indicates the approximated curves).
• Figure 3. Shape and dimensions of flat dies used for AZ91 magnesium alloy forging.
• Figure 4. Shape and dimensions of rhombic-trapezoid dies used for AZ91 magnesium alloy forging.
• Figure 5. Flat die elongation process diagram.
• Figure 6. Trapezoid-rhombic die elongation operation diagram.
• Figure 7. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the first pass, with a relative reduction of 40%, in flat dies; (b) Temperature distribution over the longitudinal section of an AZ91 alloy forging deformed within the first pass, with a relative reduction of 40%, in flat dies; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forging deformed within the first forging pass, with a relative reduction of 40%, in flat dies.
• Figure 8. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the third pass, with a relative reduction of 40%, in flat dies; (b) Temperature distribution over the longitudinal section of an AZ91 alloy forging deformed within the third pass, with a relative reduction of 40%, in flat dies; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the third forging pass, with a relative reduction of 40%, in flat dies.
• Figure 9. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the fourth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging; (b) Temperature distribution over the longitudinal section of an AZ91 alloy forging deformed within the fourth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the fourth forging pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging.
• Figure 10. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the sixth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging; (b) Temperature distribution over the longitudinal section of an AZ91 alloy forging deformed within the sixth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the sixth forging pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging.
• Figure 11. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the first pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (b) Temperature distribution over the longitudinal section of an AZ91 alloy forging deformed within the first pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the first forging pass, with a relative reduction of 25%, in trapezoid-rhombic dies.
• Figure 12. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the second pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (b) Temperature distribution over the longitudinal-section of an AZ91 alloy forging deformed within the second pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the second forging pass, with a relative reduction of 25%, in trapezoid-rhombic dies.
• Figure 13. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the third pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging; (b) Temperature distribution over the longitudinal-section of an AZ91 alloy forging deformed within the third pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the third forging pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging.
• Figure 14. (a) Temperature distribution over the cross-section of an AZ91 alloy forging deformed within the fourth pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging; (b) Temperature distribution over the longitudinal-section of an AZ91 alloy forging deformed within the fourth pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging; (c) Axonometric view of the temperature distribution over the surface of an AZ91 alloy forged within the fourth forging pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging.
• Figure 15. (a) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the first pass, with a relative reduction of 40%, in flat dies; (b) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the third pass, with a relative reduction of 40%, in flat dies; (c) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the fourth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging; (d) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the sixth pass, with a relative reduction of 40%, in flat dies, after a 90° turning of the forging.
• Figure 16. (a) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the first pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (b) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the second pass, with a relative reduction of 25%, in trapezoid-rhombic dies; (c) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the third pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging; (d) Distribution of strain intensities over the cross-section of an AZ91 alloy forged within the fourth pass, with a relative reduction of 40%, in trapezoid-rhombic dies, after a 90° turning of the forging.

7. Conclusion:

수행된 테스트 결과 분석을 바탕으로 다음과 같은 최종 결론을 도출했습니다. – 경제적 및 공정 측면에서, 저자들이 설계한 마름모-사다리꼴 형상 앤빌에서 AZ91 마그네슘 합금 바 연신 작업을 수행하는 것이 합리적입니다. 이는 평면 앤빌의 경우 6회였던 단조 패스 수를 형상 앤빌의 경우 4회로 줄여주기 때문입니다. 단조 패스 수 감소는 완제품의 최종 치수에 영향을 미치지 않습니다. – 평면 앤빌 공정과 비교하여 형상 앤빌을 연신 작업에 적용하면 중심 단조부와 재료-공구 접촉 표면에 의한 합금 영역 사이의 온도 분포가 더 균일해져 재료의 소성 유동에 더 나은 영향을 미칩니다. – 평면 및 형상 앤빌 모두에서 연신 작업을 수행할 때, 재료-공구 접촉 표면에 인접한 영역에서 바 온도가 크게 감소했습니다. 이 감소의 원인은 가열되지 않은 앤빌로 열이 전달되었기 때문입니다. – 마름모-사다리꼴 형상 앤빌에서 연신하는 동안 바 전체에 걸쳐 매우 균일한 변형률 강도 값 분포를 찾을 수 있습니다. 그러나 변형률 강도 값은 평면 앤빌에서 연신할 때 얻은 값보다 낮습니다. 더 높은 변형률 강도 값은 좋은 재료 가공을 증명하며, 이는 완제품의 더 나은 기계적 특성을 얻는 결과로 이어집니다.

8. References:

• 1. Mordike, B.L.; Ebert, T. Magnesium: Properties—Applications—Potential. Mater. Sci. Eng. A 2001, 302, 37–45.
  1. Agnew, S.R. Wrought magnesium: A 21st century outlook. JOM 2004, 56, 20–21.
  1. Bettles, C.; Gibson, M. Current wrought magnesium alloys: Strengths and weaknesses. JOM 2005, 57, 46–49.
  1. Yang, Z.; Li, J.; Zhang, J.; Lorimer, G.; Robson, J. Review on research and development of magnesium alloys. Acta Metall. Sin. (Engl. Lett.) 2008, 21, 313–328.
  1. Esmaily, M.; Svensson, J.; Fajardo, S.; Birbilis, N.; Frankel, G.; Virtanen, S.; Arrabal, R.; Thomas, S.; Johansson, L. Fundamentals and advances in magnesium alloy corrosion. Prog. Mater. Sci. 2017, 89, 92–193.
  1. Gray, J.; Luan, B. Protective coatings on magnesium and its alloys—A critical review. J. Alloys Compd. 2002, 336, 88–113.
  1. Nie, J.F. Precipitation and hardening in magnesium alloys. Metall. Mater. Trans. A 2012, 43, 3891–3939.
  1. Zeng, Z.; Stanford, N.; Davies, C.H.J.; Nie, J.F.; Birbilis, N. Magnesium extrusion alloys: A review of developments and prospects. Int. Mater. Rev. 2019, 64, 27–62.
  1. Papenberg, N.P.; Gneiger, S.; Weißensteiner, I.; Uggowitzer, P.J.; Pogatscher, S. Mg-Alloys for Forging Applications—A Review. Materials 2020, 13, 985.
  1. Braszczyńska-Malik, K. Discontinuous and continuous precipitation in magnesium-aluminium type alloys. J. Alloys Compd. 2009, 477, 870–876.
• (And 25 more references as listed in the paper)

Expert Q&A: Your Top Questions Answered

Q1: 재료 모델로 Hansel-Spittel 방정식을 선택한 특별한 이유가 있나요?

A1: 네, Hansel-Spittel 방정식은 소성 가공 공정의 수치 모델링 소프트웨어에서 유동 응력을 변형률, 변형률 속도, 온도의 함수로 기술하기 위해 빈번하게 사용되는 관계식입니다. 본 연구에서 다루는 고온 단조 공정과 같이 복잡한 열-기계적 거동을 보이는 재료의 특성을 수학적으로 잘 표현할 수 있어, 시뮬레이션의 정확도를 높이기 위해 채택되었습니다.

Q2: 논문에서 금형-소재 접촉면의 상당한 온도 강하를 언급했는데, 시뮬레이션에서 금형을 예열하지 않은 이유는 무엇인가요?

A2: 본 시뮬레이션 모델은 실제 산업 현장의 단조 공정을 최대한 현실적으로 반영하는 것을 목표로 했습니다. 실제 공정에서는 비용 및 시간 문제로 금형을 예열하지 않는 경우가 많습니다. 따라서 금형 초기 온도를 25°C로 설정하여, 공정 초기에 발생하는 소재의 냉각 효과와 이것이 재료 유동에 미치는 영향을 정확하게 분석하고자 했습니다.

Q3: Figure 15는 평면 앤빌에서 더 높은 변형률 강도를, Figure 16은 형상 앤빌에서 더 낮지만 균일한 변형률을 보여줍니다. 최종 제품에는 어느 쪽이 더 유리한가요?

A3: 논문은 하이브리드 접근법을 제안합니다. 초기 단계에서는 형상 앤빌을 사용하여 소재 중심부까지 균일하게 단조되도록 보장하는 것이 좋습니다. 균일성이 확보된 후, 최종 단계에서 평면 앤빌을 사용하여 더 높은 변형률을 가하면, 소재의 기계적 특성을 향상시키는 데 유리할 수 있습니다. 즉, 공정 단계별로 각 앤빌의 장점을 취하는 전략이 가장 이상적일 수 있습니다.

Q4: 평면 앤빌에는 40%의 압하율을, 형상 앤빌의 초기 패스에는 25%를 선택한 근거는 무엇인가요?

A4: 형상 앤빌의 경우, 그 기하학적 형태 때문에 초기 패스에서는 상부 해머 블록이 하부 앤빌에 닿기 전까지 가능한 압하율이 제한됩니다. 따라서 25%를 적용했습니다. 소재를 90° 회전시킨 후에는 더 높은 높이를 확보할 수 있어 40%의 높은 압하율 적용이 가능해졌습니다. 평면 앤빌의 경우, 변형 효과를 극대화하기 위해 40%라는 높은 상대 압하율을 선택했습니다.

Q5: 연구에 따르면 형상 앤빌은 단조 패스 수를 6회에서 4회로 줄여준다고 합니다. 이것이 경제적으로 얼마나 중요한가요?

A5: 단조 패스 수를 줄이는 것은 공정 시간, 에너지 소비, 그리고 공구 마모를 직접적으로 감소시킵니다. 이는 빌렛(소재) 제조 비용의 상당한 절감으로 이어집니다. 따라서 두 번의 패스를 줄이는 것은 대량 생산 환경에서 매우 중요한 경제적 이점을 가집니다.

Conclusion: Paving the Way for Higher Quality and Productivity

본 연구는 간단한 공구 형상의 변화가 AZ91 마그네슘 합금 단조 공정 전체에 얼마나 큰 영향을 미칠 수 있는지를 명확히 보여줍니다. 마름모-사다리꼴 형상 앤빌은 단조 패스 수를 줄여 생산 효율성을 높이고, 동시에 제품 전체에 걸쳐 균일한 품질을 확보하는 데 매우 효과적인 해결책임이 입증되었습니다. 이는 곧 생산 비용 절감과 최종 제품의 신뢰성 향상으로 이어집니다.

“STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 돕는 데 전념하고 있습니다. 이 백서에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 부품에 어떻게 구현할 수 있는지 알아보십시오.”

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0442
• 이메일 : flow3d@stikorea.co.kr

Copyright Information

• This content is a summary and analysis based on the paper “Analysis of the Open Die Forging Process of the AZ91 Magnesium Alloy” by “Grzegorz Banaszek, Teresa Bajor, Anna Kawałek and Tomasz Garstka”.
• Source: https://doi.org/10.3390/ma13173873

This material is for informational purposes only. Unauthorized commercial use is prohibited. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 Hamdy Arkoub와 Miaomiao Jin이 발표한 “Impact of chemical short-range order on radiation damage in Fe-Ni-Cr alloys” 논문을 기반으로 하며, STI C&D가 기술 전문가를 위해 분석하고 요약했습니다.

키워드

• Primary Keyword: 화학적 단거리 질서 (Chemical Short-Range Order, CSRO)
• Secondary Keywords: Fe-Ni-Cr 합금, 방사선 손상, 결함 동역학, 분자동역학, 원자력 재료, 미세구조 진화

Executive Summary

• 당면 과제: 원자력 발전소 구조 재료의 방사선 손상을 예측하는 것은 매우 중요하지만, 나노스케일에서의 원자 배열, 즉 화학적 단거리 질서(CSRO)가 결함 동역학에 미치는 영향은 그동안 충분히 고려되지 않았습니다.
• 연구 방법: 세 가지 다른 CSRO 수준을 가진 Fe-Ni-Cr 합금에 대해 광범위한 분자동역학(MD) 시뮬레이션을 사용하여 방사선 조사에 따른 결함 생성, 군집화 및 화학적 혼합을 분석했습니다.
• 핵심 발견: 높은 수준의 CSRO는 결함의 확산성을 감소시키는 ‘트랩’ 역할을 하여 결함 동역학을 근본적으로 변화시키며, 이는 실제 재료에서 결함 재결합을 촉진하여 내방사선성을 향상시킬 수 있음을 시사합니다.
• 핵심 결론: 차세대 내방사선성 합금을 정확하게 모델링하고 설계하기 위해서는 재료의 미세구조 진화에 CSRO 효과를 반드시 통합해야 합니다.

당면 과제: 왜 이 연구가 CFD 전문가에게 중요한가

Fe-Ni-Cr 기반 오스테나이트강은 원자력 분야에서 핵심 구조 재료로 널리 사용됩니다. 하지만 방사선에 노출되면 방사선 강화 확산 및 이온 혼합으로 인해 원소들이 재분배되고 재료가 손상됩니다. 이러한 조성 변화가 상 불안정성을 유발한다는 점은 오래전부터 알려져 왔지만, 나노스케일에서의 원자 배열인 화학적 단거리 질서(CSRO)가 형성될 가능성과 이것이 결함 동역학에 미치는 영향은 그동안 간과되어 왔습니다. 결함의 거동은 재료 미세구조 진화의 기초가 되기 때문에, CSRO의 역할을 이해하는 것은 재료의 수명과 안전성을 예측하는 데 매우 중요합니다. 이 연구는 바로 이 지식의 공백을 메우기 위해 수행되었습니다.

연구 접근법: 방법론 분석

본 연구는 Fe0.2Ni0.5Cr0.3 합금에 대한 광범위한 분자동역학(MD) 및 정역학 시뮬레이션을 기반으로 합니다. 연구진은 먼저 메트로폴리스 몬테카를로 분자동역학(MCMD) 기법을 사용하여 열역학적으로 안정한 원자 배열을 생성하며 세 가지 뚜렷한 CSRO 수준(낮음, 중간, 높음)을 가진 초기 시스템을 만들었습니다. 이 세 가지 시스템을 사용하여 다음과 같은 분석을 수행했습니다.

1. 적층 결함 에너지(SFE) 계산: 각 CSRO 수준에서 SFE를 계산하여 CSRO가 재료의 기본적인 기계적 특성에 미치는 영향을 평가했습니다.
2. 점결함 특성 분석: 공공(vacancy) 및 격자간 원자(interstitial)의 생성 및 이동 에너지를 계산하여 CSRO가 결함의 동적 거동에 미치는 영향을 규명했습니다.
3. 방사선 손상 시뮬레이션: 총 1,200회의 5keV 손상 캐스케이드를 연속적으로 시스템에 도입하여, 방사선 조사량에 따른 프렌켈 쌍(Frenkel pairs)의 축적, 결함 클러스터 형성, 화학적 혼합 및 CSRO 자체의 변화를 추적했습니다.

이러한 체계적인 접근을 통해 CSRO가 방사선 손상 메커니즘에 미치는 복합적인 영향을 정량적으로 밝힐 수 있었습니다.

핵심 발견: 주요 결과 및 데이터

결과 1: CSRO, 재료의 기본 특성인 적층 결함 에너지(SFE)를 극적으로 변화시키다

CSRO 수준은 재료의 근본적인 특성인 적층 결함 에너지(SFE)에 결정적인 영향을 미쳤습니다. 논문의 그림 1b에서 볼 수 있듯이, CSRO 수준이 가장 낮은 무작위 고용체(α’) 상태에서는 평균 SFE가 -106.3 mJ/m²로 음수 값을 보였으나, CSRO 수준이 가장 높은 상태(α”’)에서는 31.2 mJ/m²로 양수 값으로 전환되었습니다. 이는 CSRO를 조절함으로써 재료의 격자 안정성과 전위 거동을 제어할 수 있음을 의미하며, 높은 CSRO 상태의 SFE 값은 실제 Fe-Ni-Cr 합금의 실험값 범위(13-57 mJ/m²) 내에 있어 시뮬레이션의 신뢰성을 뒷받침합니다.

결과 2: 방사선 조사 하에서 CSRO가 결함의 거동과 위치를 제어하다

방사선 조사 시뮬레이션 결과, CSRO는 생성된 결함의 양과 분포에 큰 영향을 미쳤습니다. – 결함 재결합 효율: 그림 4a는 높은 CSRO 시나리오에서 더 많은 프렌켈 쌍(NFP)이 축적됨을 보여줍니다. 이는 시뮬레이션된 결함 싱크가 없는 시스템에서는 높은 CSRO가 결함 재결합을 비효율적으로 만든다는 것을 의미합니다. 그러나 연구진은 실제 재료에 존재하는 결정립계와 같은 결함 싱크 환경에서는 CSRO 도메인이 빠르게 움직이는 격자간 원자를 일시적으로 포획하여 공공과의 재결합 확률을 높여 내방사선성을 향상시킬 것이라고 설명합니다. – 결함 클러스터의 화학적 성분: 그림 4b에 따르면, CSRO 수준이 높아질수록 격자간 원자 클러스터 내에서 Cr 원자의 비율(rcr)이 합금의 평균 비율(0.3)보다 현저히 높아졌습니다. 이는 그림 5에서 시각적으로 확인되듯이, 격자간 원자 루프가 Cr이 풍부한 CSRO 도메인 내부 또는 근처에 우선적으로 형성됨을 나타냅니다. 이는 결함의 형성과 화학적 질서가 밀접하게 연관되어 있음을 보여주는 강력한 증거입니다.

결과 3: CSRO는 정적이지 않으며, 방사선 조사에 따라 동적으로 진화하여 정상 상태에 도달한다

CSRO는 방사선 조사 하에서 고정된 값이 아니었습니다. 그림 4c는 CSRO의 동적인 변화를 명확히 보여줍니다. 초기 CSRO가 낮았던 무작위 고용체는 방사선 조사에 따라 CSRO 값이 점진적으로 증가했습니다. 반면, 초기 CSRO가 높았던 시스템에서는 방사선에 의한 원자 혼합으로 인해 CSRO 값이 감소했습니다. 흥미롭게도 두 경우 모두 약 0.34의 특정 값으로 수렴하는 경향을 보였으며, 이는 방사선 조사 조건 하에서 동적 평형 상태, 즉 ‘정상 상태(steady-state)’ CSRO가 존재함을 시사합니다.

R&D 및 운영을 위한 실질적 시사점

• 공정 엔지니어: 이 연구는 재료의 초기 CSRO 수준을 제어하는 열처리 공정 변수를 조절함으로써 최종 제품의 내방사선성을 최적화할 수 있는 가능성을 제시합니다.
• 품질 관리팀: 논문의 데이터는 전체적인 화학 조성뿐만 아니라 나노스케일의 CSRO를 특성화하는 것이 원자로 내에서의 재료 성능을 예측하는 데 중요한 새로운 품질 검사 기준이 될 수 있음을 시사합니다.
• 설계 엔지니어: 합금 설계 시 CSRO를 새로운 설계 변수로 고려할 수 있습니다. 특정 원소의 클러스터링을 유도하여 결함 동역학을 제어하고, 보다 우수한 내방사선성을 가진 신소재를 개발하는 데 이 연구 결과를 활용할 수 있습니다.

논문 상세 정보

Impact of chemical short-range order on radiation damage in Fe-Ni-Cr alloys

1. 개요:

• 제목: Impact of chemical short-range order on radiation damage in Fe-Ni-Cr alloys (Fe-Ni-Cr 합금의 방사선 손상에 대한 화학적 단거리 질서의 영향)
• 저자: Hamdy Arkoub and Miaomiao Jin
• 발행 연도: 2023
• 게재 학술지/학회: arXiv (preprint)
• 키워드: Fe-Ni-Cr alloys, chemical short-range order, radiation damage, defects, molecular dynamics

2. 초록:

나노스케일의 특수 원자 배열 형태인 화학적 단거리 질서(CSRO)는 다양한 합금에서 전위 운동 및 결함 동역학과 같은 재료 특성을 크게 변화시키는 것으로 밝혀졌다. 본 연구에서는 광범위한 분자동역학 시뮬레이션을 기반으로 Fe-Ni-Cr 합금을 사용하여 CSRO가 결함 특성 및 방사선 거동에 어떻게 영향을 미치는지 입증한다. 세 가지 CSRO 수준에 대해 조사량의 함수로서 방사선 유도 결함 생성 경향, 결함 군집화, 화학적 혼합에 관한 통계적으로 유의미한 결과가 얻어졌다. 에너지적으로 불리한 상태(음의 적층 결함 에너지)인 완벽한 무작위 고용체는 확산을 가능하게 하는 가장 강한 경향을 보이는 반면, 높은 CSRO 시나리오는 일반적으로 포획 효과로 인해 유효 결함 확산성을 감소시켜 뚜렷한 결함 동역학을 유발한다. 특히, 높은 CSRO 시나리오에서는 격자간 원자 클러스터가 Cr이 풍부하며, 격자간 원자 루프는 Cr이 풍부한 CSRO 도메인 내부 또는 근처에 우선적으로 상주한다. 또한 CSRO는 연속적인 방사선 조사 시 감소 또는 증가하는 방식으로 동적으로 진화하여 정상 상태 값에 도달하는 것으로 확인되었다. 이러한 새로운 이해는 방사선 유도 미세구조 진화를 조사하는 데 있어 CSRO 효과를 통합하는 것의 중요성을 시사한다.

3. 서론:

Fe-Ni-Cr 기반 오스테나이트강은 원자력 응용 분야에서 구조 재료로 흔히 사용된다. 방사선 강화 확산 및 이온 혼합으로 인한 합금 내 원소 재분배는 농축 합금에서 관찰되는 화학적 단거리 질서(CSRO)를 유도할 수 있다. 조성 변화가 상 불안정성을 유발하는 것으로 오랫동안 인식되어 왔지만, CSRO의 잠재적 형성은 특히 미세구조 진화의 기초가 되는 결함 동역학에 미치는 영향을 고려할 때 충분히 강조되지 않았다. CSRO는 다양한 합금에서 기계적 특성 향상을 설명하는 데 사용되어 왔으며, 조성 및 열처리 공정 변수 제어를 통해 강도, 연성, 인성의 최적 조합을 달성하도록 맞춤화될 수 있다. 또한 CSRO는 국소 원자 환경에 따라 적층 결함 에너지(SFE)를 크게 수정할 수 있다. CSRO는 전위 이동 조절에 중요할 뿐만 아니라, 수정된 위치 에너지 지형(PEL)으로 인해 일반적인 결함의 동역학적 거동에도 본질적으로 영향을 미친다. 이 연구는 Fe-Ni-Cr 합금의 방사선 손상에 대한 CSRO의 영향을 밝히는 것을 목표로 한다.

4. 연구 요약:

연구 주제의 배경:

Fe-Ni-Cr 합금은 원자력 재료로 중요하지만 방사선에 의해 성능이 저하된다. 이러한 손상의 근본 원인인 결함 동역학에 나노스케일 원자 배열(CSRO)이 미치는 영향은 아직 명확히 규명되지 않았다.

이전 연구 현황:

이전 연구들은 주로 거시적인 조성 변화에 초점을 맞추었으며, CSRO가 결함 거동에 미치는 미시적인 영향을 체계적으로 다루지 않았다. 일부 연구에서 CSRO가 기계적 특성이나 SFE에 미치는 영향이 보고되었지만, 방사선 손상 과정 전반에 걸친 동적인 영향에 대한 이해는 부족했다.

연구 목적:

본 연구의 목적은 분자동역학 시뮬레이션을 통해 Fe-Ni-Cr 합금에서 CSRO 수준의 차이가 방사선에 의한 결함 생성, 축적, 군집화, 그리고 화학적 혼합 과정에 미치는 영향을 정량적으로 규명하는 것이다.

핵심 연구:

상이한 CSRO 수준(낮음, 중간, 높음)을 가진 Fe-Ni-Cr 합금 모델을 구축하고, 여기에 연속적인 방사선 손상 캐스케이드를 시뮬레이션했다. 이를 통해 CSRO 수준에 따른 (1) 결함 생성 및 축적률, (2) 결함 클러스터의 종류, 크기 및 화학적 조성, (3) 방사선 조사에 따른 CSRO 값의 동적 변화, (4) 원자들의 화학적 혼합률을 분석했다.

5. 연구 방법론

연구 설계:

세 가지 뚜렷한 CSRO 수준(α’, α”, α”’)을 가진 Fe0.2Ni0.5Cr0.3 합금 시스템을 비교 분석하는 설계 방식을 채택했다. 각 시스템에 동일한 조건의 방사선 손상 시뮬레이션을 적용하여 CSRO가 유일한 변수가 되도록 통제했다.

데이터 수집 및 분석 방법:

• CSRO 시스템 생성: 메트로폴리스 몬테카를로 분자동역학(MCMD)을 사용하여 서로 다른 CSRO 수준을 가진 원자 배열을 생성했다.
• 물성 계산: 각 시스템에 대해 적층 결함 에너지(SFE)와 점결함(공공, 격자간 원자)의 생성 및 이동 에너지를 계산했다.
• 방사선 손상 시뮬레이션: 300K에서 1,200회의 5keV 손상 캐스케이드를 연속적으로 시뮬레이션했다. 결함과 전위는 OVITO 소프트웨어를 사용하여 식별 및 분석했다.
• 데이터 분석: 조사량에 따른 프렌켈 쌍의 수, 격자간 원자 클러스터 내 Cr 분율, CSRO 파라미터(α), 원자들의 평균 제곱 변위(R²)를 추적하고 통계적으로 분석했다.

연구 주제 및 범위:

본 연구는 Fe-Ni-Cr 합금에 초점을 맞추었으며, 특히 Cr 원자의 클러스터링으로 대표되는 CSRO가 방사선 손상에 미치는 영향을 다룬다. 시뮬레이션은 저온(300K) 방사선 조사 환경을 모사하며, 장시간의 열적 확산 효과는 고려하지 않는다.

6. 주요 결과:

주요 결과:

• CSRO가 증가함에 따라 평균 적층 결함 에너지(SFE)는 -106.3 mJ/m²에서 31.2 mJ/m²로 음수에서 양수로 크게 증가했다.
• CSRO가 증가하면 공공 생성 에너지가 높아지는 경향을 보였다.
• 높은 CSRO 시나리오에서는 격자간 원자 클러스터가 Cr 원소가 풍부한 특성을 보였으며, 이들 클러스터는 Cr이 풍부한 CSRO 도메인 내부 또는 근처에 우선적으로 형성되었다.
• 방사선 조사 하에서 CSRO는 동적으로 변화했다. 무작위 고용체에서는 CSRO가 증가했고, 높은 CSRO 상태에서는 감소하여 특정 정상 상태 값으로 수렴하는 경향을 보였다.
• CSRO 수준이 높을수록 방사선에 의한 원자들의 화학적 혼합률이 낮았다.

그림 목록:

• FIG. 1. (a) Potential energy and the CSRO parameter with increasing MCMD timesteps. (b) SFE for three selected CSRO degrees, denoted as α’, α”, and α”’ in (a). (c-d) the corresponding atom configurations to α’, α”, and α”’, respectively.
• FIG. 2. (a-c) Vacancy formation energy (E_v^f) and (d-f) interstitial formation energy (E_i^f), corresponding to CSRO α’, α”, and α”’, respectively.
• FIG. 3. Migration barriers of vacancy (a, E_m^v) and interstitial (b, E_m^i), in random solution.
• FIG. 4. (a) Number of Frenkel pairs (NFP) with the shaped bands from standard deviation, (b) the fraction of Cr atoms in interstitial clusters (rcr), (c) CSRO parameter (α), and (d) squared displacements (R²) versus number of cascades (Nc), corresponding to CSRO α’, α”, and α”’, respectively.
• FIG. 5. Large interstitial clusters at Nc = 1200, for the case of CSRO α”’. Perfect atoms are removed for visualization purposes.
• FIG. 6. (a-c) Cluster size distribution at Nc = 1200, corresponding to CSRO α’, α”, and α”’, respectively. The right panel depicts the typical defect distribution and dislocations for each case. Perfect atoms are removed for visualization purposes.

7. 결론:

본 연구는 Fe-Ni-Cr 합금에서 국소적인 화학적 단거리 질서(CSRO)가 결함 특성과 방사선 거동에 강한 의존성을 가짐을 확인했다. CSRO는 상 형성보다 실험적으로 관찰하기 어렵기 때문에, 방사선 유도 미세구조 진화 논의에서 결함 반응 및 확산 동역학에 대한 CSRO의 숨겨진 효과는 대부분 무시되어 왔다. 따라서 향후 방사선 환경에 놓인 일반적인 CSRO 형성 합금의 방사선 효과를 설명하기 위해 이러한 나노스케일 특징을 고려하는 것이 매우 중요하다.

8. 참고 문헌:

• [1] G. S. Was, Fundamentals of radiation materials science: metals and alloys (Springer, 2016).
• [2] X. Chen, Q. Wang, Z. Cheng, M. Zhu, H. Zhou, P. Jiang, L. Zhou, Q. Xue, F. Yuan, J. Zhu, et al., Nature 592, 712 (2021).
• [3] P. Cao, Accounts of Materials Research 2, 71 (2021).
• [4] B. Gómez-Ferrer, I. García-Cortés, J. Marco, D. Jiménez-Rey, and R. Vila, Physical Review B 90, 220102 (2014).
• [5] R. Zhang, S. Zhao, J. Ding, Y. Chong, T. Jia, C. Ophus, M. Asta, R. O. Ritchie, and A. M. Minor, Nature 581, 283 (2020).
• [6] Z. Su, T. Shi, H. Shen, L. Jiang, L. Wu, M. Song, Z. Li, S. Wang, and C. Lu, Scripta Materialia 212, 114547 (2022).
• [7] Q.-J. Li, H. Sheng, and E. Ma, Nature Communications 10, 1 (2019).
• [8] W.-R. Jian, Z. Xie, S. Xu, Y. Su, X. Yao, and I. J. Beyerlein, Acta Materialia 199, 352 (2020).
• [9] J. Ding, Q. Yu, M. Asta, and R. O. Ritchie, Proceedings of the National Academy of Sciences 115, 8919 (2018).
• [10] S. Zhao, Y. Osetsky, G. M. Stocks, and Y. Zhang, npj Computational Materials 5, 1 (2019).
• [11] T. M. Smith, M. S. Hooshmand, B. D. Esser, F. Otto, D. W. McComb, E. P. George, M. Ghazisaeidi, and M. J. Mills, Acta Materialia 110, 352 (2016).
• [12] Z. Zhang, H. Sheng, Z. Wang, B. Gludovatz, Z. Zhang, E. P. George, Q. Yu, S. X. Mao, and R. O. Ritchie, Nature Communications 8, 1 (2017).
• [13] Y. Zhang, Y. Zhuang, A. Hu, J.-J. Kai, and C. T. Liu, Scripta Materialia 130, 96 (2017).
• [14] S. Zhao, G. M. Stocks, and Y. Zhang, Acta Materialia 134, 334 (2017).
• [15] Z. Pei, B. Dutta, F. Körmann, and M. Chen, Physical Review Letters 126, 255502 (2021).
• [16] G. Laplanche, A. Kostka, C. Reinhart, J. Hunfeld, G. Eggeler, and E. George, Acta Materialia 128, 292 (2017).
• [17] N. L. Okamoto, S. Fujimoto, Y. Kambara, M. Kawamura, Z. M. Chen, H. Matsunoshita, K. Tanaka, H. Inui, and E. P. George, Scientific Reports 6, 1 (2016).
• [18] X. Wang, C. M. Barr, K. Jin, H. Bei, K. Hattar, W. J. Weber, Y. Zhang, and K. L. More, Journal of Nuclear Materials 523, 502 (2019).
• [19] F. Zhang, S. Zhao, K. Jin, H. Xue, G. Velisa, H. Bei, R. Huang, J. Ko, D. Pagan, J. Neuefeind, et al., Physical Review Letters 118, 205501 (2017).
• [20] G. Bonny, N. Castin, and D. Terentyev, Modelling and Simulation in Materials Science and Engineering 21, 085004 (2013).
• [21] T. R. Allen, J. Gan, J. I. Cole, M. K. Miller, J. T. Busby, S. Shutthanandan, and S. Thevuthasan, Journal of Nuclear Materials 375, 26 (2008).
• [22] J. Cowley, Physical Review 77, 669 (1950).
• [23] M. Y. Lavrentiev, R. Drautz, D. Nguyen-Manh, T. Klaver, and S. Dudarev, Physical Review B 75, 014208 (2007).
• [24] Y. Wang, B. Ghaffari, C. Taylor, S. Lekakh, M. Li, and Y. Fan, Scripta Materialia 205, 114177 (2021).
• [25] H. Li, S. Xia, W. Liu, T. Liu, and B. Zhou, Journal of nuclear materials 439, 57 (2013).
• [26] A. Miodownik, Calphad 2, 207 (1978).
• [27] G. Henkelman, B. P. Uberuaga, and H. Jónsson, The Journal of chemical physics 113, 9901 (2000).
• [28] M. Jin, P. Cao, and M. P. Short, Acta Materialia 147, 16 (2018).
• [29] B. Million, J. Ruzickova, and J. Vrestal, Materials Science and Engineering 72, 85 (1985).
• [30] N. Astm et al., Annu. B. ASTM Stand. 12, E521 (1996).
• [31] J. F. Ziegler and J. P. Biersack, in Treatise on heavy-ion science (Springer, 1985) pp. 93–129.
• [32] A. Stukowski, Modelling and simulation in materials science and engineering 18, 015012 (2009).
• [33] M. Jin, Y. Gao, Y. Zhang, C. Jiang, and J. Gan, Scientific reports 11, 1 (2021).
• [34] M. Jin, P. Cao, and M. P. Short, Acta Materialia 186, 587 (2020).

전문가 Q&A: 자주 묻는 질문

Q1: 이 연구에서 세 가지 특정 CSRO 수준(α’, α”, α”’)을 선택한 이유는 무엇입니까?

A1: 세 가지 CSRO 수준은 각각 낮음(거의 무작위), 중간, 높음(강한 Cr 클러스터링) 상태를 대표하기 위해 선택되었습니다. 이를 통해 연구진은 CSRO의 정도가 방사선 손상에 미치는 영향을 체계적으로 조사할 수 있었습니다. 무작위 고용체 상태부터 강한 화학적 질서를 가진 상태까지 비교함으로써, CSRO가 결함 동역학과 미세구조 진화에 미치는 영향을 명확하게 분리하여 분석할 수 있었습니다.

Q2: 논문에 따르면 높은 CSRO가 시뮬레이션에서 결함 재결합을 비효율적으로 만든다고(그림 4a) 하는데, 이는 내방사선성에 불리한 것 아닌가요?

A2: 좋은 지적입니다. 이 결과는 결함 싱크(defect sink)가 없는 이상적인 시스템에서의 결과입니다. 실제 재료에는 결정립계나 전위와 같은 결함 싱크가 다수 존재합니다. 논문에서는 실제 재료의 경우, 높은 CSRO 도메인이 빠르게 확산하는 격자간 원자를 일시적으로 포획하는 ‘트랩’ 역할을 할 것이라고 설명합니다. 이로 인해 격자간 원자가 시스템 내에 더 오래 머물게 되어 공공(vacancy)과 재결합할 확률이 높아지므로, 전체적인 내방사선성은 오히려 향상될 것으로 기대됩니다.

Q3: CSRO가 증가함에 따라 적층 결함 에너지(SFE)가 음수에서 양수로 변하는 것(그림 1b)은 결함 구조에 어떤 의미를 가집니까?

A3: SFE의 부호 변화는 재료의 격자 안정성과 밀접한 관련이 있습니다. 무작위 고용체(α’)의 음수 SFE는 FCC 격자가 불안정하여 적층 결함이 쉽게 형성될 수 있음을 의미합니다. 반면, 높은 CSRO 상태(α”’)의 양수 SFE는 열역학적으로 적층 결함을 형성하기 더 어려워짐을 나타냅니다. 하지만 시뮬레이션 결과, 높은 CSRO에서도 프랭크 루프(Frank loops)와 같은 결함 구조가 형성되었는데, 이는 열역학적 요인 외에도 Cr이 풍부한 도메인 내에서 격자간 원자가 축적되는 ‘동역학적’ 요인이 결함 구조 형성을 주도하기 때문입니다.

Q4: 높은 CSRO 시나리오에서 격자간 원자 클러스터가 Cr이 풍부하다는 발견의 중요성은 무엇입니까?

A4: 이 발견은 화학적 질서와 결함 거동 사이에 강한 상호작용이 있음을 보여줍니다. 결함이 단순히 무작위로 생성되고 군집하는 것이 아니라, Cr이 풍부한 특정 화학적 환경을 선호한다는 것을 의미하기 때문입니다. 이는 방사선 유도 편석(radiation-induced segregation)과 같은 현상을 이해하고, 합금의 미세구조가 방사선 하에서 어떻게 진화할지를 예측하는 데 매우 중요한 단서를 제공합니다.

Q5: 연구에 따르면 CSRO는 방사선 조사 하에서 정상 상태(steady state)로 진화한다고 합니다(그림 4c). 이 정상 상태 값은 무엇에 의해 결정됩니까?

A5: 논문은 이 정상 상태가 두 가지 경쟁적인 효과의 균형에 의해 결정된다고 설명합니다. 하나는 질서를 파괴하려는 경향의 ‘비열적 방사선 혼합(athermal radiation mixing)’이고, 다른 하나는 질서를 형성하려는 경향의 ‘방사선 강화 열 확산(radiation-enhanced thermal diffusion)’입니다. 최종적인 정상 상태 CSRO 값은 초기 상태보다는 방사선 입자의 종류, 선속, 온도와 같은 방사선 조사 조건에 따라 결정될 것으로 예측됩니다.

결론: 더 높은 품질과 생산성을 향한 길

본 연구 분석은 화학적 단거리 질서(CSRO)가 원자력 재료의 방사선 손상 메커니즘을 이해하는 데 있어 지금까지 간과되었던 핵심 요소임을 명확히 보여줍니다. 나노스케일에서의 원자 배열이 결함의 생성, 이동, 군집화에 직접적인 영향을 미치며, 이는 곧 재료의 전체적인 내구성과 수명을 결정합니다. CSRO를 이해하고 제어하는 것은 더 안전하고 신뢰성 높은 차세대 원자력 재료를 개발하기 위한 필수적인 과제입니다.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

저작권 정보

• 이 콘텐츠는 Hamdy Arkoub와 Miaomiao Jin의 논문 “Impact of chemical short-range order on radiation damage in Fe-Ni-Cr alloys”를 기반으로 한 요약 및 분석 자료입니다.
• 출처: https://arxiv.org/abs/2301.09719

이 자료는 정보 제공 목적으로만 사용됩니다. 무단 상업적 사용을 금지합니다. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 Pham Quang이 SCIREA Journal of Materials (2025)에 발표한 논문 “Numerical modeling and simulation of mold filling and solidification process of BS100 Grade A6 alloy in sand casting of excavator bucket coupling parts”를 기반으로 하며, 기술 전문가를 위해 STI C&D가 분석하고 요약했습니다.

Keywords

• Primary Keyword: 주조 시뮬레이션
• Secondary Keywords: 샌드 캐스팅, 금형 충전, 응고 해석, BS100 Grade A6, 굴삭기 버킷 커플링, JSCast

Executive Summary

• The Challenge: 고품질, 고강도 부품에 사용되는 탄소강의 샌드 캐스팅 공정은 응고 과정에서 발생하는 수축 결함으로 인해 품질 저하의 위험이 있습니다.
• The Method: JSCast 시뮬레이션 소프트웨어를 사용하여 BS100 Grade A6 합금으로 제작되는 굴삭기 버킷 커플링 부품의 금형 충전 및 응고 과정을 정밀하게 분석했습니다.
• The Key Breakthrough: 시뮬레이션 결과, 설계된 게이팅 시스템을 통해 16.12초 만에 안정적인 충전이 완료되었으며, 라이저(riser)가 최대 676.17초까지 용탕을 공급하여 내부 수축 및 기공 결함을 효과적으로 방지함을 확인했습니다.
• The Bottom Line: 주조 시뮬레이션은 실제 금형 제작 및 시험 주조에 앞서 게이팅 및 라이저 시스템 설계를 최적화하여 잠재적 결함을 예측하고 방지함으로써 제품 품질과 생산 효율성을 극대화하는 핵심 도구입니다.

The Challenge: Why This Research Matters for CFD Professionals

건설 및 광산 산업에서 널리 사용되는 굴삭기 버킷 커플링과 같은 고강도 부품은 높은 품질과 내구성이 요구됩니다. 이러한 부품은 주로 샌드 캐스팅 기법으로 생산되는데, 이 공정은 용융 금속이 응고되는 과정에서 필연적으로 부피가 수축하여 발생하는 ‘수축 결함’에 매우 취약합니다. 특히 BS100 Grade A6과 같은 고강도 합금은 응고 특성이 민감하여 불완전한 충전, 탕경(cold shut), 열간 균열(hot tearing) 등의 문제가 발생할 가능성이 높습니다. 이러한 결함은 제품의 기계적 성능을 저하시키고, 심각한 경우 현장에서의 파손으로 이어질 수 있습니다. 따라서 실제 생산에 들어가기 전에 이러한 결함을 예측하고 방지할 수 있는 신뢰성 높은 방법론의 필요성이 절실합니다.

The Approach: Unpacking the Methodology

본 연구에서는 이러한 문제를 해결하기 위해 JSCast 시뮬레이션 소프트웨어를 활용하여 BS100 Grade A6 합금의 샌드 캐스팅 공정에 대한 수치 모델링 및 해석을 수행했습니다.

연구진은 먼저 Solidwork 소프트웨어를 사용하여 굴삭기 버킷 커플링의 3D 지오메트리 모델을 생성했습니다. 이후 JSCast 소프트웨어에서 총 4,874,553개의 요소를 가진 메쉬를 생성하여 정밀한 해석을 준비했습니다. 초기 충전 시간은 Chvorinov’s rule을 기반으로 약 18초로 예측되었으며, 이를 바탕으로 초기 조건을 설정했습니다.

해석에는 BS100 Grade A6 합금과 샌드 몰드의 열물성 데이터(밀도, 비열, 열전도도 등)가 정확하게 입력되었습니다. 시뮬레이션은 용융 금속이 주입되어 금형 캐비티를 채우는 ‘금형 충전’ 단계와, 채워진 용탕이 열을 잃고 고체로 변하는 ‘응고’ 단계로 나누어 진행되었습니다. 이를 통해 시간에 따른 온도 분포, 유동 패턴, 응고 진행 과정을 시각적으로 분석하고 잠재적인 결함 발생 영역을 예측했습니다.

The Breakthrough: Key Findings & Data

시뮬레이션은 금형 설계의 유효성을 입증하고 결함 방지 메커니즘을 명확히 보여주는 두 가지 핵심 결과를 도출했습니다.

Finding 1: 16.12초 만에 완료된 안정적이고 균일한 금형 충전

금형 충전 시뮬레이션 결과, 설계된 게이팅 시스템이 용융 금속을 매우 안정적이고 효율적으로 캐비티 내부로 유도하는 것으로 나타났습니다. Figure 2에서 볼 수 있듯이, 충전은 총 16.12초 만에 완료되었습니다. 충전 과정 동안 온도 분포는 1584.8 °C에서 1600.0 °C 사이로 거의 균일하게 유지되었으며, 얇은 벽이나 말단부에서만 약간의 냉각이 관찰되었습니다. 특히 에어록(airlock), 동결(freeze-up), 미성형(misrun)과 같은 주요 유동 결함이 발생하지 않았는데, 이는 러너 시스템이 부드럽고 완전한 충전을 촉진하도록 잘 설계되었음을 의미합니다.

Finding 2: 라이저를 통한 효과적인 수축 보상 및 응고 제어

응고 시뮬레이션은 결함 방지의 핵심인 응고 패턴을 명확하게 보여주었습니다. Figure 3에 나타난 바와 같이, 응고는 금형 벽에서 시작하여 제품 중심부로 점진적으로 진행되었습니다. 얇은 부분들이 먼저 응고되었고, 두꺼운 부분과 라이저는 오랫동안 용융 상태를 유지했습니다. 특히 라이저는 최대 676.17초까지 용융 금속을 머금고 있으면서, 주조품 본체가 응고하며 수축하는 동안 지속적으로 용탕을 공급하는 ‘압탕(feeding)’ 역할을 효과적으로 수행했습니다. Figure 3의 응고 분율 곡선은 주조품(Casting)이 라이저(Riser)보다 훨씬 빨리 응고됨을 명확히 보여주며, 이는 내부 수축 및 기공 결함의 위험을 크게 줄이는 이상적인 응고 거동입니다.

Practical Implications for R&D and Operations

본 연구 결과는 주조 공정 관련 엔지니어들에게 다음과 같은 실질적인 시사점을 제공합니다.

• For Process Engineers: 이 연구는 잘 설계된 게이팅 및 러너 시스템이 균일한 충전과 결함 최소화에 얼마나 중요한지를 보여줍니다. 시뮬레이션 결과는 현재의 설계가 유효함을 입증하며, 유사한 부품 설계 시 이 시스템을 벤치마킹할 수 있습니다.
• For Quality Control Teams: Figure 3의 응고 순서 데이터는 잠재적인 핫스팟, 즉 가장 마지막에 응고되어 내부 기공이 발생할 가능성이 가장 높은 영역을 정확히 알려줍니다. 이는 비파괴 검사(NDT) 시 집중적으로 검사해야 할 부위를 특정하여 품질 관리의 효율성을 높이는 데 활용될 수 있습니다.
• For Design Engineers: 라이저의 위치와 크기가 수축 결함 보상에 결정적인 역할을 한다는 점이 시뮬레이션을 통해 입증되었습니다. 이는 설계 초기 단계에서 주조 시뮬레이션을 활용하여 라이저 시스템을 최적화하는 것이 물리적 테스트를 줄이고 개발 기간을 단축하는 가장 효과적인 방법임을 시사합니다.

Paper Details

Numerical modeling and simulation of mold filling and solidification process of BS100 Grade A6 alloy in sand casting of excavator bucket coupling parts

1. Overview:

• Title: Numerical modeling and simulation of mold filling and solidification process of BS100 Grade A6 alloy in sand casting of excavator bucket coupling parts
• Author: Pham Quang
• Year of publication: 2025
• Journal/academic society of publication: SCIREA Journal of Materials
• Keywords: Numerical modeling, Jscast, sand casting, mold filling, solidification, excavator bucket coupling, BS100 Grade A6.

2. Abstract:

Carbon steel bucket couplings are widely used in the construction and mining industries due to their quality and versatility. These couplings are manufactured using sand casting techniques to meet international standards. This casting technique has significant shrinkage defects during solidification and needs to be improved. In this study, JSCast simulation software was used to analyze the mold filling and solidification process of BS100 Grade A6 alloy during sand casting of excavator bucket couplings. The simulation provides important insights into the thermal and flow dynamics, allowing prediction and prevention of potential casting defects. The mold filling simulation demonstrated a steady and efficient flow of molten metal into the mold cavity, achieving full filling within 16.12 seconds. The temperature distribution throughout the filling process remained almost uniform, with minor cooling observed in the thinner and peripheral regions. There were no major flow defects such as airlocks, freeze-ups, or misruns, indicating that the runner system was well designed to promote smooth and complete filling. Solidification simulations showed a solidification pattern that progressed inward from the mold wall toward the core. The thinner sections solidified earlier, while the risers held the molten metal for a longer time—up to 676.17 seconds—allowing for effective shrinkage zones. The delayed solidification of the shrinkage beans relative to the casting was confirmed by the solidification fraction curve, confirming their role in reducing the risk of internal shrinkage and porosity.

3. Introduction:

Sand casting is a widely used manufacturing process for producing metal components, suitable for both ferrous and non-ferrous materials [1,2]. This method offers several advantages, including excellent dimensional accuracy, ease of pattern development, faster production rates, and reduced solidification time compared to die casting. Key process parameters in sand casting include venting, gating system design, nozzle and riser configuration, and the specification of dimensional and pattern tolerances. These parameters play a crucial role in determining the quality of the final casting. Additionally, the bonding mechanism of the mold material is influenced by several factors such as physical and mechanical properties, crystal structure, and bonding compounds [2-6]. The quality and mechanical performance of the cast specimen are heavily dependent on the proper selection of both process parameters and input materials. In this context, a comprehensive literature review is essential to evaluate and predict the hardness of the specimen. This is particularly important for analyzing the mold filling behavior of BS100 Grade A6 alloy in the sand casting of excavator bucket coupling parts, with the ultimate goal of enhancing product quality [7]. Modeling and simulation have become indispensable tools in modern sand casting, allowing engineers and designers to analyze and optimize the casting process before experimental production.

4. Summary of the study:

Background of the research topic:

굴삭기 버킷 커플링과 같은 탄소강 부품은 건설 및 광산 산업에서 널리 사용되며 샌드 캐스팅으로 제조됩니다. 이 공정은 응고 중 발생하는 수축 결함 문제를 개선해야 할 필요가 있습니다.

Status of previous research:

ProCAST, ANSYS Fluent, MAGMASOFT와 같은 시뮬레이션 소프트웨어는 유동, 열전달, 응고를 시뮬레이션하여 금형 충전, 냉각 속도, 수축 기공 등 결함 관련 현상을 예측하는 데 사용되어 왔습니다. 이전 연구들은 스프루 시스템 최적화, 곡선 게이트 사용 등을 통해 주조 품질을 향상시키는 데 시뮬레이션이 효과적임을 보여주었습니다.

Purpose of the study:

본 연구의 목적은 JSCast 시뮬레이션 소프트웨어를 사용하여 굴삭기 버킷 커플링의 샌드 캐스팅 공정에서 BS100 Grade A6 합금의 금형 충전 및 응고 과정을 분석하는 것입니다. 이를 통해 잠재적인 주조 결함을 예측 및 방지하고, 열 및 유동 역학에 대한 중요한 통찰력을 얻고자 합니다.

Core study:

BS100 Grade A6 합금으로 만들어지는 굴삭기 버킷 커플링의 샌드 캐스팅 공정을 대상으로 금형 충전 및 응고 과정을 수치적으로 모델링하고 시뮬레이션했습니다. 시뮬레이션을 통해 충전 시간, 온도 분포, 유동 패턴, 응고 순서 등을 분석하여 설계된 게이팅 및 라이저 시스템의 유효성을 평가했습니다.

5. Research Methodology

Research Design:

3D 지오메트리 모델(Solidworks)과 유한 요소 메쉬(JSCast)를 기반으로 한 수치 시뮬레이션 접근법을 사용했습니다. BS100 Grade A6 합금의 금형 충전 및 응고 거동을 예측하기 위해 열-유동 연성 해석을 수행했습니다.

Data Collection and Analysis Methods:

BS100 Grade A6 합금과 샌드 몰드의 열물성 데이터(Table 3, 4)를 입력값으로 사용했습니다. JSCast 소프트웨어를 통해 시간에 따른 온도 분포, 충전율, 응고 분율 등의 데이터를 계산하고 시각화하여 분석했습니다.

Research Topics and Scope:

연구는 BS100 Grade A6 합금을 사용한 굴삭기 버킷 커플링의 샌드 캐스팅 공정에 국한됩니다. 주요 분석 대상은 금형 충전 거동과 응고 패턴이며, 이를 통해 주조 결함 발생 가능성을 예측하고 방지하는 데 중점을 둡니다.

6. Key Results:

Key Results:

• 금형 충전: 용융 금속은 16.12초 내에 안정적이고 효율적으로 금형 캐비티를 완전히 채웠습니다. 충전 과정 동안 온도 분포는 거의 균일하게 유지되었으며, 에어록이나 미성형과 같은 주요 유동 결함은 관찰되지 않았습니다.
• 응고 과정: 응고는 금형 벽에서 시작하여 중심부로 진행되는 방향성 응고 패턴을 보였습니다. 얇은 부분이 먼저 응고되었고, 라이저는 최대 676.17초까지 용융 상태를 유지하며 주조품의 수축을 효과적으로 보상했습니다.
• 결함 방지: 응고 분율 곡선은 라이저가 주조품보다 훨씬 늦게 응고됨을 확인시켜 주었으며, 이는 내부 수축 및 기공 결함의 위험을 줄이는 데 라이저 시스템이 효과적으로 기능했음을 입증합니다.

Figure List:

• Figure 1. Geometry, element of the Casting and Sand Mold.
• Figure 2. Mold Filling Behavior Using JSCast
• Figure 3. The solidification simulation of the BS100 Grade A6 alloy casting

7. Conclusion:

The simulation results confirm that the mold design, gating system, and and riser placement are effective for casting a BS100 Grade A6 alloy bucket coupling. The pouring phase of the liquid metal into the mold shows smooth, uninterrupted metal flow with favorable heat distribution, while the solidification sequence confirms successful solidification promotion with proper direction and feed of the risers. These findings highlight the importance of integrated simulation in optimizing casting quality, reducing defects, and improving process efficiency prior to physical testing. Further refinement could focus on minimizing cooling temperature differences in thin-walled regions and validating the simulation results through experimental casting tests.

8. References:

• [1] Monroe, R. Porosity in Castings, AFS Transactions. 2005, 4, 05-245.
• [2] R.W. Monroe, “Use of Iron Oxide in Mold and Core Mixes for Ferrous Castings,” AFS Transactions, vol. 93, pp.355-364, 1985.
• [3] Gunasegaram, D.R. Farnsworth, D.J. Nguyena, T.T. Identification of critical factors affecting shrinkage porosity in permanent mould casting using numerical simulations based on design of experiments, Journal of materials processing technology, 2009, 209, 1209–1219. DOI:10.1016/j.jmatprotec.2008.03.044
• [4] Reis, A. Houbaert, Y. Zhian Xu, Rob Van Tol, Santosa, A.D. Duartea, J.F. Magalhaes, A.B. Modeling of shrinkage defects during solidification of long and short freezing materials, Journal of materials processing technology.2008, 202, 428–434. https://doi.org/10.1016/j.jmatprotec.2007.10.030
• [5] Grozdanić, V. Modelling of the directional solidification of a leaded red brass flange, Materials and technology, 2011, 45, 47-53. http://www.dlib.si/?URN=URN:NBN:SI:doc-LSRDY63C
• [6] Carlson, K.D. Lin, Z. Hardin, R.A. and Beckermann, C. Modeling of Porosity Formation and Feeding Flow in Steel Casting, SFSA 56th Technical and Operating Conference, November 2002.
• [7] https://castingquality.com/
• [8] Kumar, R., & Rao, T. S. Optimization of gating system for steel castings using simulation. International Journal of Advanced Manufacturing Technology, 2017, 90(1–4), 543–550. https://doi.org/10.1007/s00170-016-9373-3
• [9] Liu, Y., Wang, H., & Chen, M. Gating system design and defect analysis in aluminum alloy casting using ProCAST simulation. Journal of Materials Processing Technology, 2019, 271, 228–236. https://doi.org/10.1016/j.jmatprotec.2019.03.020
• [10] Zhang, J., Li, W., & Sun, X. Simulation-based prediction of hardness distribution in cast steel using thermal and microstructure modeling. Materials Science and Engineering: A, 2021, 799, 140164. https://doi.org/10.1016/j.msea.2020.140164.
• [11] Singh, R., Sharma, V., & Pathak, S. Optimization of mold filling parameters in sand casting using JSCast simulation. Procedia Manufacturing, 2020, 48, 372–378. https://doi.org/10.1016/j.promfg.2020.05.050
• [12] Chvorinov, N. “Theorie der Erstarrung von Gussstücken” [Theory of the Solidification of Castings]. Giesserei (in German), 1940,. 27 (10–12): 177–186, 201–208, 222–225.
• [13] JSCAST Ver.20@, QUALICA Inc.

Expert Q&A: Your Top Questions Answered

Q1: 이 연구에서 JSCast 소프트웨어를 선택한 특별한 이유가 있나요?

A1: 논문에 따르면, JSCast는 사용자가 주입 온도, 게이트 형상, 몰드 벤팅과 같은 파라미터를 조정하여 균일하고 완전한 금형 충전을 보장할 수 있게 해줍니다. 이는 복잡한 형상의 주조품에서 충전 균일성을 개선하는 데 시뮬레이션이 얼마나 효과적인지를 보여준 Singh 등의 연구(2020) 사례에서도 알 수 있습니다. 따라서 복잡한 형상의 굴삭기 부품에 대한 유동 및 응고를 최적화하는 본 연구의 목적에 JSCast가 적합한 도구로 선택되었습니다.

Q2: Chvorinov’s rule에 따른 예비 충전 시간(18초)과 시뮬레이션 결과(16.12초) 사이에 차이가 있는데, 이 차이는 무엇을 의미하나요?

A2: Chvorinov’s rule은 부피와 표면적을 기반으로 한 이론적인 응고 시간 예측 공식으로, 초기 설계의 타당성을 검토하는 데 유용합니다. 반면, 시뮬레이션은 실제 유체의 역학과 열전달을 고려한 훨씬 정밀한 예측입니다. 두 값(18초와 16.12초)이 근사하다는 것은 초기 설계 가정이 합리적이었음을 보여주며, 시뮬레이션 결과가 더 짧게 나온 것은 설계된 게이팅 시스템이 이론적 계산보다 더 효율적으로 충전을 완료시켰음을 의미합니다.

Q3: Figure 2의 금형 충전 단계에서 발생할 수 있는 가장 큰 잠재적 위험은 무엇이었나요?

A3: 논문에 따르면, 충전이 6.23초(39% 충전) 진행되었을 때 복잡한 유동 경로가 나타나고, 내부 모서리나 얇은 벽 영역이 더 빨리 냉각되기 시작했습니다. 또한, 잠재적인 난류 및 유동 선단의 불안정성도 언급되었습니다. 이는 제대로 관리되지 않을 경우, 이러한 복잡하고 빠르게 냉각되는 영역에서 조기 응고나 공기 혼입(air entrapment)이 발생할 수 있었음을 시사합니다.

Q4: Figure 3은 라이저의 효과를 어떻게 정량적으로 보여주나요?

A4: Figure 3의 응고 분율 곡선 그래프는 ‘Casting(주조품)’과 ‘Riser(라이저)’의 응고 진행을 두 개의 개별 선으로 명확히 보여줍니다. 주조품의 곡선이 100%에 도달하여 응고가 완료된 후에도 라이저의 곡선은 한참 동안 낮은 분율에 머물러 있습니다. 논문에서 최종 응고 시간이 최대 676.17초라고 언급했는데, 이는 라이저가 이 시간 동안 용융 상태를 유지하며 주조품에 용탕을 공급했음을 의미합니다. 이 시간 차이가 바로 수축을 보상하는 라이저의 효과를 정량적으로 입증하는 데이터입니다.

Q5: 이 연구가 구체적으로 방지하고자 했던 결함은 무엇이며, 시뮬레이션은 이를 어떻게 확인시켜 주었나요?

A5: 초록과 서론에서 명확히 밝히듯, 이 연구의 주요 목표는 ‘수축 결함(shrinkage defects)’과 ‘내부 수축 및 기공(internal shrinkage and porosity)’을 개선하는 것이었습니다. 응고 시뮬레이션(Figure 3)은 금형 벽에서 중심으로 향하는 방향성 응고를 보여주었고, 라이저가 가장 마지막에 응고됨으로써 주조품의 체적 감소를 효과적으로 보상했음을 입증했습니다. 이를 통해 수축 결함이 성공적으로 방지되었음을 확인시켜 주었습니다.

Conclusion: Paving the Way for Higher Quality and Productivity

결론적으로, 본 연구는 고강도 굴삭기 부품의 샌드 캐스팅 공정에서 발생하는 고질적인 수축 결함 문제를 주조 시뮬레이션을 통해 어떻게 해결할 수 있는지를 명확히 보여줍니다. JSCast를 활용한 정밀한 해석은 금형 충전이 16.12초 내에 안정적으로 완료되고, 라이저 시스템이 676초 이상 용탕을 공급하여 내부 결함을 효과적으로 억제함을 입증했습니다. 이는 시뮬레이션이 단순한 예측 도구를 넘어, 실제 생산 전에 설계를 최적화하고 품질을 확보하며 개발 기간과 비용을 절감하는 필수적인 엔지니어링 솔루션임을 의미합니다.

STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 지원하는 데 전념하고 있습니다. 이 백서에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 부품에 어떻게 구현할 수 있는지 알아보십시오.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

Copyright Information

• This content is a summary and analysis based on the paper “Numerical modeling and simulation of mold filling and solidification process of BS100 Grade A6 alloy in sand casting of excavator bucket coupling parts” by “Pham Quang”.
• Source: https://doi.org/10.54647/materials430289

This material is for informational purposes only. Unauthorized commercial use is prohibited. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 [Jan Majernik, Stefan Gaspar, Martin Podaril, Tomas Coranic]이 [MM SCIENCE JOURNAL]에 발표한 “[EVALUATION OF THERMAL CONDITIONS AT CAST-DIE CASTING MOLD INTERFACE]” ([2020]) 논문을 기반으로 합니다. STI C&D의 기술 전문가가 분석하고 요약했습니다.

키워드

• Primary Keyword: 다이캐스팅 금형 계면
• Secondary Keywords: 고압 다이캐스팅(HPDC), 금형 열 균형, 열 응력, 온도 분포

Executive Summary

• 도전 과제: 다이캐스팅 공정에서 금형 계면에 발생하는 극심한 온도 변화는 주기적인 열 응력을 유발하여 금형 수명을 단축시키고 주조 품질을 저하시킵니다.
• 연구 방법: 본 연구는 전기 모터 플랜지 주조에 사용되는 고압 다이캐스팅 금형의 열 균형을 이론적 계산과 시뮬레이션(Magmasoft MAGMA5)을 결합하여 평가했습니다.
• 핵심 발견: 일반적인 예상과 달리, 금형 온도 조절 시스템은 열을 방출하는 것이 아니라 오히려 금형에 열을 공급하는 역할을 했으며(-200,080 J), 가장 큰 열 방출 요인은 이형제 스프레이(1,055,300 J)였습니다.
• 핵심 결론: 정확한 열 균형 분석, 특히 온도 조절 시스템의 가열 역할과 이형제 스프레이의 지배적인 냉각 효과를 이해하는 것은 다이캐스팅 공정을 최적화하고 주조 품질을 향상시키는 데 필수적입니다.

도전 과제: 이 연구가 CFD 전문가에게 중요한 이유

고압 다이캐스팅(HPDC) 기술에서 상대적으로 차가운 금형은 매우 높은 온도의 용탕과 짧은 시간 동안 반복적으로 접촉합니다. 이로 인해 금형 표면에는 극심한 온도 변동이 발생하며, 이는 주기적인 열 응력으로 이어져 금형의 피로 파괴나 균열의 원인이 됩니다. 안정적인 주조 품질을 확보하기 위해서는 공정의 안정성을 보장해야 하며, 이를 위해 금형의 열 균형을 정확하게 파악하고 제어하는 것이 무엇보다 중요합니다.

하지만 주조, 금형, 그리고 주변 환경 사이에서 일어나는 열 전달 과정은 매우 빠르고 복잡하며 비정상(non-stationary) 상태로 진행되기 때문에, 이를 정확히 계산하고 예측하는 것은 상당한 기술적 과제입니다. 기존에는 금형 온도 조절 시스템이 주로 냉각 역할을 할 것으로 간주했지만, 전체 시스템의 열 균형을 종합적으로 평가하지 않으면 공정 최적화에 한계가 있었습니다.

접근 방식: 연구 방법론 분석

본 연구는 EN AC 47100 합금을 사용하여 전기 모터 플랜지를 생산하는 Müller Weingarten 600 고압 다이캐스팅 머신의 금형을 대상으로 열 균형을 평가했습니다. 연구팀은 주조품이 방출하는 총 열량과 금형에서 주변으로 방출되는 총 열량을 계산하기 위해 포괄적인 이론적 접근법을 사용했습니다.

• 주조품의 열 방출량 (Q) 계산: 용탕의 과열(Qoverheat), 응고 잠열(Qsolid), 냉각(Qcooling) 과정에서 방출되는 열에너지의 합으로 총 열량을 계산했습니다.
• 금형의 열 방출량 (QALL) 계산: 금형에서 주변으로 방출되는 열을 복사(Qrad), 대류(Qfl), 온도 조절 시스템(Qt), 기계 프레임으로의 전도(Qfr), 이형제 스프레이(Qspr) 등 개별 구성 요소로 나누어 각각의 열량을 계산했습니다.
• 온도 분포 시뮬레이션: 실제 측정의 어려움을 보완하기 위해 Magmasoft MAGMA5 – HPDC 모듈을 사용하여 금형 내부의 깊이별(표면에서 1mm, 2mm, 5mm, 10mm, 20mm) 온도 변화를 시뮬레이션했습니다.

이러한 이론적 계산과 시뮬레이션 결과를 비교 분석하여 다이캐스팅 금형 계면에서 일어나는 복잡한 열 전달 메커니즘을 규명했습니다.

핵심 발견: 주요 결과 및 데이터

결과 1: 온도 조절 시스템의 예상 밖의 ‘가열’ 역할

연구 결과, 가장 놀라운 발견 중 하나는 금형 온도 조절 시스템의 역할이었습니다. 일반적으로 냉각을 담당할 것으로 예상했던 것과 달리, 이 시스템은 오히려 금형에 열을 공급하는 것으로 나타났습니다.

• Table 4에 따르면, 온도 조절 시스템을 통해 전달된 열량(Qt)은 -200,080 J로 음수 값을 보였습니다. 이는 시스템이 금형으로부터 열을 제거하는 것이 아니라, 사이클 동안 다른 요인에 의해 식는 것을 방지하고 안정적인 기준 온도를 유지하기 위해 열을 공급했음을 의미합니다. 이는 금형 자체가 반복적인 가열을 통해 열을 축적하는 ‘축열체’로 작동하기 때문입니다.

결과 2: 이형제 스프레이의 지배적인 열 방출 효과

금형에서 열을 방출하는 여러 요인 중 이형제 스프레이의 영향이 압도적으로 큰 것으로 확인되었습니다.

• Table 4와 Figure 6의 데이터에 따르면, 이형제 스프레이에 의해 방출된 열량(Qspr)은 1,055,300 J에 달했습니다. 이는 주조품이 방출한 총 열량(Q = 1,001,555.86 J)을 초과하는 수치로, 전체 열 방출량의 105.7%를 차지합니다. 반면, 복사(Qrad, 5.7%)나 전도(Qfr, 5.5%)에 의한 열 방출은 상대적으로 미미했습니다. 이 결과는 이형제 스프레이가 단순히 이형을 돕는 것을 넘어, 금형 표면을 급격히 냉각시키는 가장 중요한 공정 변수임을 명확히 보여줍니다.

결과 3: 금형 표면층에 집중된 극심한 열 응력

시뮬레이션 결과, 금형의 열 응력은 캐비티 표면에 가까운 얇은 층에 집중되는 것으로 나타났습니다.

• Table 5는 금형 깊이에 따른 온도 변화(ΔT)를 보여줍니다. 주 러너(MR_ED) 지점에서 표면으로부터 1mm 깊이의 온도 변화는 242.6°C에 달했지만, 20mm 깊이에서는 29.8°C에 불과했습니다. 이처럼 표면 근처에서 발생하는 200°C 이상의 급격한 온도 변화는 해당 부위에 높은 주기적 열 응력을 유발하여 금형 수명에 직접적인 영향을 미침을 시사합니다.

R&D 및 운영을 위한 실질적 시사점

• 공정 엔지니어: 이 연구는 온도 조절 시스템이 단순히 냉각 장치가 아니라, 안정적인 공정 온도를 유지하기 위한 ‘가열’ 장치로도 기능할 수 있음을 시사합니다. 또한, 이형제 스프레이의 분사량과 패턴이 금형의 열 균형에 가장 큰 영향을 미치므로, 이를 정밀하게 제어하는 것이 주조 품질을 안정시키는 핵심 요소가 될 수 있습니다.
• 품질 관리팀: 논문의 Table 5 데이터는 금형 표면으로부터의 거리가 온도 변동에 미치는 영향을 명확히 보여줍니다. 이는 열 응력으로 인한 미세 균열 등의 결함이 주로 표면층에 집중될 수 있음을 의미하며, 새로운 품질 검사 기준을 수립하는 데 중요한 정보를 제공합니다.
• 설계 엔지니어: 본 연구의 결론은 온도 조절 채널의 설계에 대한 명확한 가이드라인을 제시합니다. “다이캐비티 벽과 온도 조절 채널 사이의 거리가 짧을수록 금형 내 온도 편차가 줄어든다”는 발견은, 채널을 가능한 한 캐비티에 가깝게 배치하는 것이 내부 응력을 줄이고 균형 잡힌 열 환경을 만드는 데 유리함을 의미합니다.

논문 상세 정보

EVALUATION OF THERMAL CONDITIONS AT CAST-DIE CASTING MOLD INTERFACE

1. 개요:

• 제목: EVALUATION OF THERMAL CONDITIONS AT CAST-DIE CASTING MOLD INTERFACE
• 저자: JAN MAJERNIK, STEFAN GASPAR, MARTIN PODARIL, TOMAS CORANIC
• 발행 연도: 2020
• 학술지/학회: MM SCIENCE JOURNAL
• 키워드: HPDC, temperature distribution, mold material, tempering

2. 초록:

주조품의 내부 구조는 용탕이 금형 표면에 접촉할 때의 과냉각 수준에 따라 달라집니다. 다이캐스팅 기술에서는 상대적으로 차가운 금형 재료가 짧은 연속 사이클 동안 상당히 높은 온도의 용탕과 접촉합니다. 이는 높은 열 변동을 유발하여 금형의 주기적인 열 응력을 야기합니다. 본 논문은 특정 유형의 주조품을 주조하는 동안 금형의 열 균형을 평가하고, 한 번의 주조 사이클 동안 금형 체적 내의 온도 분포를 설명합니다. 이 논문은 주조품, 금형 및 주변 환경 사이에서 발생하는 열 과정에 대한 포괄적인 그림을 제시합니다.

3. 서론:

주조품의 내부 구조는 용탕이 금형 표면에 접촉할 때의 과냉각 수준에 따라 결정됩니다. 다이캐스팅 기술에서는 상대적으로 차가운 금형 재료가 짧은 연속 사이클 동안 상당히 높은 온도의 용탕과 접촉합니다. 이는 높은 열 변동을 유발하여 금형의 주기적인 열 응력을 야기합니다. 본 논문은 특정 유형의 주조품을 주조하는 동안 금형의 열 균형을 평가하고, 한 번의 주조 사이클 동안 금형 체적 내의 온도 분포를 설명하며, 주조품, 금형 및 주변 환경 사이에서 발생하는 열 과정에 대한 포괄적인 그림을 제시합니다.

4. 연구 요약:

연구 주제의 배경:

고압 다이캐스팅 공정에서 금형과 용탕 사이의 급격한 열 교환은 주조품 품질과 금형 수명에 결정적인 영향을 미칩니다. 안정적인 생산을 위해서는 금형의 열 균형을 정량적으로 이해하고 제어하는 것이 필수적입니다.

이전 연구 현황:

기존 연구들은 금형의 열 전달 및 온도 분포에 대해 다루었지만, 온도 조절 시스템, 이형제 스프레이, 복사, 전도 등 모든 열 전달 요소를 종합적으로 고려한 정량적 열 균형 분석은 부족했습니다.

연구 목적:

본 연구의 목적은 실제 고압 다이캐스팅 공정에서 금형의 전체적인 열 균형을 평가하고, 시뮬레이션을 통해 금형 내부의 동적 온도 분포를 시각화하여 다이캐스팅 금형 계면의 열적 거동에 대한 깊이 있는 이해를 제공하는 것입니다.

핵심 연구:

• 특정 주조품(전기 모터 플랜지) 생산 시 금형의 열 균형을 계산하기 위한 수학적 모델을 제시하고 검증합니다.
• 주조품에서 금형으로 전달되는 열량과 금형에서 외부로 방출되는 열량을 각 구성 요소별(복사, 대류, 전도, 온도 조절, 스프레이)로 정량화합니다.
• CAE 시뮬레이션을 이용하여 금형 표면으로부터 깊이에 따른 온도 변화를 분석하고, 이를 통해 열 응력이 집중되는 영역을 파악합니다.

5. 연구 방법론

연구 설계:

본 연구는 이론적 계산과 CAE 시뮬레이션을 결합한 접근 방식을 채택했습니다. 먼저, 열역학 공식을 기반으로 한 주조 사이클의 전체 열 균형 모델을 수립했습니다. 그 다음, 이 모델의 결과를 검증하고 금형 내부의 상세한 온도 분포를 파악하기 위해 Magmasoft MAGMA5 – HPDC 모듈을 사용한 시뮬레이션을 수행했습니다.

데이터 수집 및 분석 방법:

• 이론적 계산: 주조품의 열 방출량과 금형의 열 방출량을 계산하기 위해 논문에 제시된 공식 (1)부터 (17)까지를 사용했습니다. 계산에 필요한 변수들(재료 물성, 공정 조건 등)은 Table 3에 정리된 값을 적용했습니다.
• 시뮬레이션: 실제 공정 조건(Table 2)을 입력하여 금형의 특정 지점들(Figure 3)에서 시간에 따른 온도 변화를 시뮬레이션하고, 그 결과를 Table 5와 Figure 4로 분석했습니다.

연구 주제 및 범위:

연구는 EN AC 47100 합금을 사용하는 고압 다이캐스팅 공정에 초점을 맞추었습니다. 분석 대상은 전기 모터 플랜지 주조용 금형이며, 한 번의 주조 사이클 동안 발생하는 열 균형과 금형 내부의 온도 분포를 평가하는 것으로 범위를 한정했습니다.

6. 주요 결과:

주요 결과:

• 주조품이 방출한 총 열량(Q = 1,001,555.86 J)과 금형의 각 구성 요소가 방출한 총 열량(QALL = 998,924.1 J)이 거의 일치하여 계산 모델의 정확성이 검증되었습니다.
• 금형 온도 조절 시스템은 열을 방출하지 않고 오히려 금형에 열을 공급(Qt = -200,080 J)하는 것으로 나타났습니다.
• 금형의 열 방출에서 가장 큰 비중을 차지하는 것은 이형제 스프레이(Qspr = 1,055,300 J)였습니다.
• 금형 표면 1mm 깊이에서는 최대 242.6°C의 큰 온도 변화가 발생했으나, 20mm 깊이에서는 29.8°C로 안정적이어서 열 응력이 표면에 집중됨을 확인했습니다.

Figure 목록:

• Figure 1. High pressure die casting machine Müller Weingarten 600
• Figure 2. Temperature measurement points after mold opening
• Figure 3. Temperature distribution measurement points
• Figure 4. Temperature development Main Runner/Cover Die
• Figure 5. Temperature gradient Main Runner / Cover Die
• Figure 6. Transferred heat from the mold

7. 결론:

1. 단순화된 열 균형 계산 방법론이 검증되었으며, 주조품에서 전달된 열(Q)과 금형에서 방출된 총 열(QALL)이 일치함을 통해 그 정확성이 확인되었습니다.
2. 주조 산업에서 CAE 지원을 통해 금형 체적 내 온도 분포를 시각화하는 것은 매우 유용합니다.
3. 온도 분포 평가 시, 절대적인 온도 값뿐만 아니라 주조 사이클 동안의 전체적인 온도 변화 과정과 온도 구배를 함께 고려해야 주조품의 내부 구조를 예측하는 데 도움이 됩니다.
4. 온도 조절 채널의 분포와 설계는 금형 내 온도장 변화와 관련 열 방출에 상당한 영향을 미칩니다. 다이캐비티 벽과 채널 사이의 거리가 짧을수록 금형 내 온도 차이가 작아지므로, 채널을 가능한 한 캐비티에 가깝게 설계하는 것이 바람직합니다.

8. 참고 문헌:

• [Gaspar 2012] Gaspar, S., et al. Dependence of pressure die casting quality on die casting plunger velocity inside a filling chamber of a pressure die casting machine. Journal of Computational and Theoretical Nanoscience, 2012, Vol. 14, No. 1, pp. 499-502. ISSN 1936-6612.
• [Gaspar 2017] Gaspar, S., et al. Influence of structure adjustment of gating system of casting mould upon the quality of die cast. Lüdenscheid: RAM Verlag, 2017. ISBN 978-3-942303-47-7.
• [Klobcar 2008] Klobcar, D., et al. Thermal stresses in aluminium alloy die casting dies. Computational Materials Science, 2008, Vol. 43, No. 4. ISSN 1147-1154.
• [Krenicky 2011] Krenicky, T. Implementation of Virtual Instrumentation for Machinery Monitoring. In: Scientific Papers: Operation and Diagnostics of Machines and Production Systems Operational States: Vol. 4, RAM-Verlag, Lüdenscheid, 2011, pp. 5-8. ISBN 978-3-942303-10-1.
• [Majernik 2019] Majernik, J., et al. Evaluation of the Temperature Distribution of a Die Casting Mold of X38CrMoV5_1 Steel. Archives of Foundry Engineering, 2019, Vol. 19, No. 1, pp. 107-112. ISSN 2299-2944.
• [Mascenik 2019] Mascenik, J., Pavlenko, S. Innovative broad Spectrum testing and monitoring of belt transmissions. MM Science Journal, 2018, Vol. 11, No. 5, pp. 2706-2709. ISSN 1803-1269.
• [Novakova 2006] Novakova, I. The thermal and pressure conditions in the system casting – casting tool during the high pressure die casting of the alloys of alluminium with silicon: dissertation thesis. Liberec, Czech Republic: Technical University of Liberec, 2006.
• [Pasko 2014] Pasko, J., Gaspar, S. Technological Factors of Die Casting. Lüdenscheid: RAM Verlag, 2014. ISBN 978-3-942303-25-5.
• [Turtelli 2006] Turtelli et al. Magnetic and structural characterization of as-cast and annealed melt- spun Fe80-xSi20Crx. Journal of Magnetism and Magnetic Materials, 2006, Vol. 304, No. 2, pp. E687-E689.
• [Wladysiak 2017] Wladysiak, R., Kozun, A., Pacyniak, T. The Effect of Water Mist Cooling of Casting Die on the Solidification, Microstructure and Properties of AlSi20 Alloy. Archives of Metallurgy and Materials, 2017, Vol. 62, No. 1, pp. 187-194. ISSN 1733-3490.
• [Yu 2016] Yu, W.-B., et al. Development and application of inverse heat transfer model between liquid melt and shot sleeve in high pressure die casting process under non-shooting condition. China Foundry, 2016, Vol. 13, No. 4, pp. 269-275. ISSN 1672-6421.
• [Yu 2017] Yu, W.-B., et al. Interfacial heat transfer behavior at metal/die in finger-plated casting during high pressure die casting process. China Foundry, 2017, Vol. 14, No. 4, pp. 258-264. ISSN 1672-6421.
• [Zhi-Peng 2008] Zhi-Peng, G., et al. Determination of the heat transfer coefficient at metal-die interface of high pressure die casting process of AM50 alloy. International Journal of Heat and Mass Transfer, 2008, Vol. 51, No. 25-26., pp. 6032-6038. ISBN 978-3-942303-47-7.

전문가 Q&A: 자주 묻는 질문

Q1: 계산에서 피스톤의 운동 에너지로 발생하는 열(Qfrict)을 무시한 이유는 무엇인가요?

A1: 논문의 2.1절에 따르면, 피스톤의 마찰과 운동 에너지로 인해 발생하는 열의 양은 “매우 작아서 무시할 수 있다(quite small and can be neglected)”고 명시되어 있습니다. 전체 열 균형에서 차지하는 비중이 미미하기 때문에 계산의 복잡성을 줄이고 주요 열 전달 요인에 집중하기 위해 생략되었습니다.

Q2: 연구 결과에서 온도 조절 시스템이 금형에 열을 공급(Qt = -200,080 J)한다고 나왔는데, 그 이유는 무엇인가요?

A2: 이는 금형이 단일 사이클 동안 용탕으로부터 받은 열을 축적하는 ‘축열체’처럼 작동하기 때문입니다. 금형이 열려 있는 동안 이형제 스프레이와 주변 공기에 의해 냉각되는데, 온도 조절 시스템은 이 냉각 효과를 상쇄하고 다음 사이클을 위해 안정적이고 높은 기준 온도를 유지하는 역할을 합니다. 따라서 이 특정 공정 조건에서는 냉각이 아닌 가열 기능을 수행한 것입니다.

Q3: 금형에서 열을 방출하는 가장 중요한 요인은 무엇이었나요?

A3: Table 4와 Figure 6에서 명확히 보여주듯이, 이형제 스프레이(Qspr)가 가장 지배적인 열 방출 요인이었습니다. 그 양은 1,055,300 J로, 주조품이 방출한 총 열량을 넘어서는 강력한 냉각 효과를 보였습니다. 이는 다이캐스팅 공정에서 이형제 스프레이가 금형의 열 관리에 얼마나 중요한 변수인지 강조합니다.

Q4: 금형 캐비티 표면으로부터의 거리가 열 응력에 어떤 영향을 미치나요?

A4: Table 5의 데이터는 거리에 따라 온도 변화(ΔT)가 극적으로 감소함을 보여줍니다. 표면에서 1mm 깊이에서는 ΔT가 240°C를 초과했지만, 20mm 깊이에서는 약 30°C에 불과했습니다. 이는 주기적인 열 응력이 금형의 극히 얇은 표면층에 집중되어 있으며, 이 부위가 피로 파괴에 가장 취약하다는 것을 의미합니다.

Q5: 이 연구 결과를 바탕으로 온도 조절 채널 설계에 대한 실질적인 권장 사항은 무엇인가요?

A5: 논문의 결론 4번 항목에서 명확히 제시하고 있습니다. “다이캐비티 벽과 온도 조절 채널 사이의 거리가 짧을수록 금형 체적 내의 온도 차이가 작아진다”는 것입니다. 따라서 금형 내부의 온도 분포를 균일하게 하고 내부 응력을 줄이기 위해서는, 온도 조절 채널을 가능한 한 캐비티 표면에 가깝게 배치하여 설계하는 것이 매우 효과적입니다.

결론: 더 높은 품질과 생산성을 향한 길

다이캐스팅 금형 계면에서의 복잡한 열적 거동을 이해하는 것은 고품질 주조품을 안정적으로 생산하기 위한 핵심 과제입니다. 본 연구는 이론적 계산과 시뮬레이션을 통해 금형의 전체 열 균형을 정량적으로 분석하고, 온도 조절 시스템이 가열 역할을 할 수 있다는 점과 이형제 스프레이가 열 방출에 가장 지배적인 영향을 미친다는 중요한 사실을 밝혔습니다. 이러한 통찰은 공정 최적화, 금형 수명 연장, 그리고 최종 제품의 품질 향상으로 이어질 수 있는 실질적인 단서를 제공합니다.

STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 돕는 데 전념하고 있습니다. 이 논문에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 부품에 어떻게 구현할 수 있는지 논의해 보십시오.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

저작권 정보

• 이 콘텐츠는 “[Jan Majernik, Stefan Gaspar, Martin Podaril, Tomas Coranic]”의 논문 “[EVALUATION OF THERMAL CONDITIONS AT CAST-DIE CASTING MOLD INTERFACE]”를 기반으로 한 요약 및 분석 자료입니다.
• 출처: https://doi.org/10.17973/MMSJ.2020_11_2020041

이 자료는 정보 제공 목적으로만 사용됩니다. 무단 상업적 사용을 금지합니다. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 Shantanu Shahane 외 저자가 2018년 ASME(미국기계학회)에 발표한 논문 “VIRTUALLY-GUIDED CERTIFICATION WITH UNCERTAINTY QUANTIFICATION APPLIED TO DIE CASTING”을 기반으로 하며, STI C&D의 기술 전문가에 의해 분석 및 요약되었습니다.

키워드

• Primary Keyword: 다이캐스팅 불확실성 정량화
• Secondary Keywords: 가상 인증, 응고 시뮬레이션, 미세구조 예측, 수치 해석, 매개변수 불확실성

Executive Summary

• The Challenge: 다이캐스팅 공정은 수많은 변수가 제품 품질에 영향을 미치지만, 이러한 공정 변수의 확률적 변동(stochastic variation)을 기존의 결정론적 시뮬레이션만으로는 완벽하게 예측하기 어렵습니다.
• The Method: 실험 데이터, 예측 모델(응고, 미세구조, 기계적 물성), 그리고 불확실성 정량화(UQ)를 결합한 ‘가상 인증 프레임워크’를 제안하여 공정 변수의 변동이 최종 제품 품질에 미치는 영향을 시뮬레이션했습니다.
• The Key Breakthrough: 경계 조건 온도와 잠열의 1% 표준편차가 최종 제품의 응고 시간, 2차 덴드라이트 암 간격(SDAS), 항복 강도에 미치는 영향을 정량적으로 예측하는 반응 표면(response surfaces)을 성공적으로 생성했습니다.
• The Bottom Line: 불확실성 정량화 기법을 다이캐스팅 시뮬레이션에 적용하면, 공정 변수의 미세한 변동이 제품의 기계적 특성에 미치는 영향을 통계적으로 예측할 수 있어 보다 신뢰도 높은 품질 관리가 가능합니다.

The Challenge: Why This Research Matters for CFD Professionals

다이캐스팅은 자동차 및 주택 산업의 핵심 부품을 생산하는 중요한 제조 공정입니다. 공정의 복잡성으로 인해 합금 물성, 금형과의 계면 조건, 열 경계 조건 등 수많은 공정 변수를 정밀하게 측정하고 제어하기는 매우 어렵습니다.

기존의 결정론적(deterministic) 컴퓨터 시뮬레이션은 이러한 변수들이 고정된 값이라고 가정하기 때문에, 실제 현장에서 발생하는 공정 변수의 미세한 확률적 변동이 최종 제품 품질에 미치는 영향을 완벽하게 평가하는 데 한계가 있습니다. 특히 재료의 미세구조와 물리적 특성 변동 간의 상호작용에 대한 지식 격차도 존재합니다. 이러한 불확실성은 시뮬레이션 결과와 실제 제품 간의 오차를 유발하며, 이는 곧 제품의 신뢰성 문제로 이어질 수 있습니다. 본 연구는 이러한 문제를 해결하기 위해 검증, 검인 및 불확실성 정량화(V&V/UQ)를 통합한 포괄적인 프레임워크를 제안합니다.

The Approach: Unpacking the Methodology

본 연구는 다이캐스팅 공정의 불확실성을 정량화하기 위해 실험과 시뮬레이션을 유기적으로 결합한 ‘가상 인증 프레임워크’를 구축했습니다. 이 프레임워크는 다음과 같은 핵심 모듈로 구성됩니다.

1. 응고 모델 (Solidification Model): 자연 대류로 인한 유체 유동과 열전달, 응고 효과를 포함하는 수치 모델을 사용했습니다. C++ 기반의 객체 지향 소프트웨어인 OpenCast를 개발했으며, 비정형 격자에서 유한 체적법을 사용하여 지배 방정식을 해석했습니다. 특히, 다공성 영역(mushy zone)에서의 유동 저항을 모델링하기 위해 Darcy 항을 운동량 방정식에 포함시켰습니다.
2. 미세구조 및 재료 거동 모델 (Microstructure and Material Behavior Model): 응고 시뮬레이션에서 계산된 온도 구배와 냉각 속도를 입력받아 2차 덴드라이트 암 간격(SDAS)과 2% 항복 강도(σ0.2)를 예측하는 경험적 모델을 사용했습니다. SDAS는 Backer와 Wang이 제안한 관계식을, 항복 강도는 Okayasu 등이 제안한 관계식을 활용했습니다.
3. 불확실성 정량화 (Uncertainty Quantification, UQ): 공정 변수의 확률적 변동을 다루기 위해 비침투적(non-intrusive) 방법인 확률적 배치법(stochastic collocation)을 채택했습니다. 이 방법은 기존의 결정론적 시뮬레이션 코드를 수정할 필요 없이 래퍼(wrapper)처럼 작동합니다. 입력 변수(경계 온도, 잠열)의 확률 분포(정규 분포)에 따라 Smolyak 알고리즘을 사용하여 최적의 샘플 포인트를 생성하고, 각 포인트에서 결정론적 시뮬레이션을 실행하여 결과의 반응 표면을 구성했습니다.

이러한 통합 프레임워크를 통해, 연구진은 실험 데이터로 모델을 보정 및 검증하고, 검증된 모델에 UQ를 적용하여 공정 변동성의 영향을 정량적으로 분석했습니다.

The Breakthrough: Key Findings & Data

Finding 1: 수치 모델의 신뢰성 검증 완료

연구진은 개발한 OpenCast 소프트웨어의 정확성을 검증하기 위해, 차등 가열된 정육면체 내부의 3차원 자연 대류 문제에 대한 기존 연구 결과(Fusegi et al. [19])와 비교했습니다. 레일리 수(Ra) 10⁵ 및 10⁶ 조건에서 온도 및 속도 프로파일을 비교한 결과, 온도 오차는 0.02% 미만, 속도 오차는 0.8% 미만으로 매우 높은 일치도를 보였습니다. 이는 본 연구에서 사용된 응고 시뮬레이션 솔버가 유체 유동 및 열전달 현상을 정확하게 예측할 수 있음을 입증하며, 후속 다이캐스팅 해석 결과의 신뢰성을 뒷받침합니다.

• Figure 4와 Figure 5에서 볼 수 있듯이, Ra = 10⁵ 조건에서 온도 및 X, Y 방향 속도 프로파일은 기존 연구 데이터와 거의 완벽하게 일치했습니다. 예를 들어, Z=0.5, X=0.5에서 Y 방향 속도 프로파일의 오차는 0.3246%에 불과했습니다.

Finding 2: 공정 변수의 불확실성이 제품 특성에 미치는 영향 정량화

연구진은 창문 프레임 커넥터 리브(rib) 형상에 대해 경계 온도와 잠열이라는 두 가지 입력 변수에 각각 평균값의 1%에 해당하는 표준편차를 부여하여 불확실성 정량화 해석을 수행했습니다. 그 결과, 이 두 변수의 변동이 응고 시간, 최대 SDAS, 최소 항복 강도라는 세 가지 출력 변수에 미치는 영향을 시각적으로 보여주는 반응 표면(Response Surface)을 성공적으로 도출했습니다.

• Figure 10 (응고 시간): 잠열이 높을수록, 그리고 금형 벽 온도가 높을수록 응고 시간이 길어지는 경향을 명확히 보여줍니다. 잠열이 3.85×10⁵ J/kg에서 3.95×10⁵ J/kg으로 증가하고 벽 온도가 495K에서 505K로 증가함에 따라 응고 시간은 약 0.75초에서 0.78초로 증가했습니다.
• Figure 12 (최소 항복 강도): 응고 시간이 길어질수록 항복 강도는 감소하는 역비례 관계를 보였습니다. 잠열과 벽 온도가 증가함에 따라 최소 항복 강도는 약 144.1 MPa에서 143.8 MPa로 감소하는 경향이 나타났습니다. 이는 냉각 속도가 느려지면 결정립이 성장하여(SDAS 증가, Figure 11) 기계적 강도가 저하됨을 정량적으로 보여줍니다.

Practical Implications for R&D and Operations

• For Process Engineers: 본 연구는 금형 온도와 같은 공정 변수의 미세한 변동이 최종 제품의 항복 강도에 직접적인 영향을 미침을 보여줍니다. 반응 표면(Figure 12)을 통해 목표 항복 강도를 유지하기 위한 공정 변수의 허용 오차 범위를 설정하는 데 기여할 수 있습니다.
• For Quality Control Teams: 논문의 Figure 9는 제품의 특정 위치(두꺼운 부위)에서 SDAS가 크고 항복 강도가 낮게 예측됨을 보여줍니다. 이는 미세구조(SDAS)와 기계적 물성(항복 강도) 간의 상관관계를 나타내며, 미세구조 분석을 통해 제품의 취약부를 예측하고 새로운 품질 검사 기준을 수립하는 데 정보를 제공할 수 있습니다.
• For Design Engineers: 해석 결과는 제품의 두꺼운 부분이 가장 늦게 응고되며, 결과적으로 결정립이 크고(높은 SDAS) 항복 강도가 낮아지는 경향을 명확히 보여줍니다. 이는 설계 초기 단계에서 제품의 두께 균일성을 확보하는 것이 응고 결함 및 기계적 특성 저하를 방지하는 데 중요한 고려사항임을 시사합니다.

Paper Details

VIRTUALLY-GUIDED CERTIFICATION WITH UNCERTAINTY QUANTIFICATION APPLIED TO DIE CASTING

1. Overview:

• Title: VIRTUALLY-GUIDED CERTIFICATION WITH UNCERTAINTY QUANTIFICATION APPLIED TO DIE CASTING
• Author: Shantanu Shahane, Soham Mujumdar, Namjung Kim, Pikee Priya, Narayana Aluru, Placid Ferreira, Shiv G Kapoor, Surya Vanka
• Year of publication: 2018
• Journal/academic society of publication: Proceedings of the ASME 2018 Verification and Validation Symposium
• Keywords: Die casting, Uncertainty Quantification, Solidification, Microstructure, Virtual Certification

2. Abstract:

Die casting is a type of metal casting in which liquid metal is solidified in a reusable die. In such a complex process, measuring and controlling the process parameters is difficult. Conventional deterministic simulations are insufficient to completely estimate the effect of stochastic variation in the process parameters on product quality. In this research, a framework to simulate the effect of stochastic variation together with verification, validation, and uncertainty quantification is proposed. This framework includes high-speed numerical simulations of solidification, micro-structure and mechanical properties prediction models along with experimental inputs for calibration and validation. Both experimental data and stochastic variation in process parameters with numerical modeling are employed thus enhancing the utility of traditional numerical simulations used in die casting to have a better prediction of product quality. Although the framework is being developed and applied to die casting, it can be generalized to any manufacturing process or other engineering problems as well.

3. Introduction:

제조 산업은 다양한 공정을 사용하여 완제품을 생산하며, 이러한 공정들은 종종 매우 복잡하고 제품 품질에 영향을 미치는 많은 공정 변수를 포함한다. 다이캐스팅은 액체 금속을 재사용 가능한 금형에서 응고시켜 원하는 부품 형상을 제작하는 제조 공정 중 하나이다. 자동차 및 주택 산업이 다이캐스트 제품의 주요 소비자이다. 컴퓨터 시뮬레이션과 실험은 공정의 물리학을 이해하고 제품 품질을 개선하는 데 자주 사용된다. 그러나 복잡한 공정에서 이러한 공정 변수들을 측정하고 제어하는 것은 어렵다. 더욱이, 기존의 결정론적 시뮬레이션은 공정 변수의 확률적 변동이 제품 품질에 미치는 영향을 완전히 추정하기에는 불충분하다. 따라서 본 논문에서는 다이캐스팅 제품의 품질을 정확하게 예측하기 위한 포괄적인 도구를 얻기 위해 가상 유도 인증 방법론을 제안한다.

4. Summary of the study:

Background of the research topic:

다이캐스팅 공정은 복잡성으로 인해 공정 변수의 측정 및 제어가 어렵고, 이 변수들의 확률적 변동이 제품 품질에 미치는 영향을 기존의 결정론적 시뮬레이션만으로는 정확히 예측하기 어렵다.

Status of previous research:

기존의 시뮬레이션 기법은 수치 알고리즘의 발전과 광범위한 물성 데이터에도 불구하고 상당한 불확실성과 오차를 내포하고 있다. 특히 재료 미세구조와 물성 변동 간의 상호작용에 대한 지식 격차가 존재한다.

Purpose of the study:

실험 데이터, 예측 모델, 불확실성 정량화를 연결하는 프레임워크를 통해 수치 시뮬레이션의 정확도를 높이는 것을 목표로 한다. 이 연구는 검증, 검인, 불확실성 정량화를 통해 수치 시뮬레이션 결과의 불확실성을 줄이는 데 중점을 둔다.

Core study:

응고, 미세구조 및 기계적 특성 예측 모델을 포함하는 고속 수치 시뮬레이션과 보정 및 검증을 위한 실험 입력을 결합한 가상 인증 프레임워크를 제안하고 다이캐스팅에 적용한다. 이 프레임워크는 공정 변수의 확률적 변동을 수치 모델링과 함께 사용하여 제품 품질 예측을 향상시킨다.

5. Research Methodology

Research Design:

본 연구는 실험, 수치 모델링, 불확실성 정량화를 통합한 가상 인증 프레임워크를 설계했다. 실험을 통해 얻은 데이터(미세구조, 항복 강도 등)는 수치 모델의 경험적 관계식을 보정하고 검증하는 데 사용된다. 수치 모델은 응고, 미세구조, 기계적 거동 예측 모듈로 구성된다. 불확실성 정량화는 결정론적 소프트웨어의 래퍼(wrapper) 역할을 하여 공정 변수의 확률적 변동이 최종 제품 품질에 미치는 영향을 평가한다.

Data Collection and Analysis Methods:

• 실험: Aallied Die Casting(RCM Industries)에서 제공한 A380 알루미늄 합금 다이캐스트 부품을 사용하여 공간적으로 변화하는 미세구조 데이터를 수집하고, 단축 인장 시험을 통해 항복 강도 및 인장 강도 데이터를 획득했다.
• 수치 해석: C++로 개발된 OpenCast 소프트웨어를 사용하여 유한 체적법 기반의 응고 시뮬레이션을 수행했다. 선형 시스템 해석에는 HYPRE 라이브러리의 대수적 다중격자(Algebraic Multigrid) 솔버를 사용했다.
• 불확실성 정량화: 비침투적 확률적 배치법(non-intrusive stochastic collocation)을 사용했다. 입력 변수의 확률 분포(정규분포)에 따라 Smolyak 알고리즘으로 샘플 포인트를 생성하고, 각 포인트에서 결정론적 시뮬레이션을 수행하여 다항식 카오스 전개(Polynomial Chaos Expansion)의 계수를 추정했다.

Research Topics and Scope:

연구 범위는 다이캐스팅 공정에 적용되는 가상 인증 프레임워크의 개발 및 시연에 초점을 맞춘다. 주요 연구 주제는 다음과 같다. 1. 응고 및 열-유체 유동을 모사하는 수치 모델 개발 및 검증. 2. 응고 결과(온도 구배, 냉각 속도)를 기반으로 미세구조(SDAS) 및 기계적 물성(항복 강도)을 예측하는 경험적 모델의 통합. 3. 경계 온도 및 잠열과 같은 공정 변수의 불확실성이 응고 시간, SDAS, 항복 강도에 미치는 영향을 정량화.

6. Key Results:

Key Results:

• 개발된 수치 해석 소프트웨어(OpenCast)는 3차원 자연 대류 벤치마크 문제에 대해 기존 연구 결과와 0.8% 미만의 오차로 잘 일치함을 검증했다.
• 실제 다이캐스팅 부품(창문 프레임 커넥터 리브)에 대한 결정론적 응고 시뮬레이션을 성공적으로 수행하여 시간에 따른 고상 분율, SDAS, 항복 강도 분포를 예측했다.
• 2차원 불확실성 정량화 분석을 통해, 경계 온도와 잠열의 확률적 변동이 응고 시간, 최대 SDAS, 최소 항복 강도에 미치는 영향을 나타내는 반응 표면을 성공적으로 생성했다.
• 반응 표면 분석 결과, 잠열과 경계 온도가 증가할수록 응고 시간이 길어지고(냉각 속도 감소), 이는 SDAS의 증가와 항복 강도의 감소로 이어짐을 정량적으로 확인했다.

Figure List:

• FIGURE 1: FRAMEWORK
• FIGURE 2: DIE CAST SPECIMEN FOR EXPERIMENTAL VALIDATION [1]
• FIGURE 3: TYPICAL MICROSTRUCTURE OF A380 ALLOY
• FIGURE 4: TEMPERATURE PROFILES (Ra = 10⁵)
• FIGURE 5: VELOCITY PROFILES (Ra = 10⁵)
• FIGURE 6: VELOCITY PROFILES (Ra = 10⁶)
• FIGURE 7: MESH
• FIGURE 8: SOLID FRACTION ISO-SURFACES (a) 0.1166 s (b) 0.2332 s
• FIGURE 9: Solidification Simulation Output X = 0.013 m
• FIGURE 10: RESPONSE SURFACE: SOLIDIFICATION TIME
• FIGURE 11: RESPONSE SURFACE: MAXIMUM SDAS
• FIGURE 12: RESPONSE SURFACE: MINIMUM YIELD STRENGTH

7. Conclusion:

본 논문은 다이캐스팅에 적용된 가상 인증 프레임워크를 기술한다. 이 프레임워크는 수치 시뮬레이션, 실험을 통한 검증 및 검인, 그리고 불확실성 정량화를 포함한다. 수치 모델은 유체 유동 및 열전달과 함께 응고를 시뮬레이션한다. 미세구조 및 구조적 특성은 수치 시뮬레이션에서 얻은 온도 구배 및 냉각 속도의 함수로 경험적으로 추정된다. 공개된 3차원 자연 대류 수치 결과를 사용한 검증은 현재 소프트웨어의 출력이 0.8% 미만의 오차로 양호한 일치도를 보임을 보여준다. 실제 다이캐스팅 형상의 결정론적 응고는 고상 분율, SDAS 및 항복 강도 결과와 함께 시연된다. 공정 변수의 확률적 변동이 출력 변수에 미치는 영향을 정량화하기 위해 매개변수 불확실성 정량화가 활용된다. 동일한 형상에 대해 2차원 확률적 분석이 수행되고 반응 표면이 도시된다. 향후 제어된 실험 결과를 통한 검증 및 보정이 필요하며, 이러한 가상 인증 프레임워크는 제품 성능을 향상시키는 데 도움이 될 수 있고 다른 제조 공정이나 엔지니어링 문제에도 쉽게 일반화될 수 있다.

8. References:

• [1] Aallied Die Casting Company of Illinois.
• [2] Backer, G., and Wang, Q., 2007. “Microporosity simulation in aluminum castings using an integrated pore growth and interdendritic flow model”. Metallurgical and Materials Transactions B, 38(4), pp. 533-540.
• [3] Hunt, D., 1979. “Solidification and casting of metals”. The Metal Society, London.
• [4] Okayasu, M., Takeuchi, S., Yamamoto, M., Ohfuji, H., and Ochi, T., 2015. “Precise analysis of microstructural ef- fects on mechanical properties of cast adc12 aluminum al- loy”. Metallurgical and Materials Transactions A, 46(4), pp. 1597-1609.
• [5] Plotkowski, A., Fezi, K., and Krane, M., 2015. “Estimation of transient heat transfer and fluid flow for alloy solidifi- cation in a rectangular cavity with an isothermal sidewall”. Journal of Fluid Mechanics, 779, pp. 53–86.
• [6] Spiegel, E., and Veronis, G., 1960. “On the boussinesq approximation for a compressible fluid.”. The Astrophysical Journal, 131, p. 442.
• [7] Geuzaine, C., and Remacle, J.-F., 2009. “GMSH: A 3-D finite element mesh generator with built-in pre-and post- processing facilities”. International Journal for Numerical Methods in Engineering, 79(11), pp. 1309–1331.
• [8] Artemiev, M. Tethex: Conversion from tetrahedra (trian- gles) to hexahedra (quadrangles).
• [9] HYPRE: Scalable Linear Solvers and Multigrid Methods.
• [10] Baker, A. H., Falgout, R. D., Kolev, T. V., and Yang, U. M., 2012. “Scaling HYPRE’s multigrid solvers to 100,000 cores”. In High-Performance Scientific Comput- ing. Springer, pp. 261-279.
• [11] Xiu, D., and Karniadakis, G. E., 2002. “The Wiener- Askey polynomial chaos for stochastic differential equa- tions”. SIAM Journal on Scientific Computing, 24(2), pp. 619-644.
• [12] Wiener, N., 1938. “The homogeneous chaos”. American Journal of Mathematics, 60(4), pp. 897–936.
• [13] Matthies, H. G., and Keese, A., 2005. “Galerkin methods for linear and nonlinear elliptic stochastic partial differen- tial equations”. Computer Methods in Applied Mechanics and Engineering, 194(12-16), pp. 1295–1331.
• [14] Roman, L. J., and Sarkis, M., 2006. “Stochastic Galerkin method for elliptic SPDEs: A white noise approach”. Dis- crete and Continuous Dynamical Systems-Series B, 6(4), p. 941.
• [15] Smith, R. C., 2013. Uncertainty quantification: theory, im- plementation, and applications, Vol. 12. Siam.
• [16] Smolyak, S., 1963. “Quadrature and interpolation formu- las for tensor products of certain classes of functions”. In Soviet Math. Dokl., Vol. 4, pp. 240-243.
• [17] Ganapathysubramanian, B., and Zabaras, N., 2007. “Sparse grid collocation schemes for stochastic natural convection problems”. Journal of Computational Physics, 225(1), pp. 652-685.
• [18] Heiss, F., and Winschel, V., 2008. “Likelihood approxima- tion by numerical integration on sparse grids”. Journal of Econometrics, 144(1), pp. 62–80.
• [19] Fusegi, T., Hyun, J., Kuwahara, K., and Farouk, B., 1991. “A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure”. International Journal of Heat and Mass Transfer, 34(6), pp. 1543–1557.

Expert Q&A: Your Top Questions Answered

Q1: 불확실성 정량화를 위해 많은 샘플링 기법 중 Smolyak 알고리즘을 선택한 이유는 무엇입니까?

A1: 논문에 따르면, 다차원 확률 공간에서 샘플 포인트를 선택할 때 텐서 곱(tensor product)을 사용하면 차원이 증가함에 따라 필요한 샘플 수가 기하급수적으로 늘어나 계산 비용이 매우 커집니다. Smolyak 알고리즘은 보간의 정확도를 유지하면서도 다차원 공간에서 최소한의 샘플 포인트를 선택하는 데 사용되므로, 결정론적 시뮬레이션이 시간 소모적일 때 계산 효율성을 크게 높일 수 있어 이 연구에 적합했습니다.

Q2: Figure 12의 반응 표면은 실제 다이캐스팅 공정 제어에 어떤 구체적인 의미를 가집니까?

A2: Figure 12는 금형 벽 온도와 재료의 잠열이라는 두 가지 주요 변수가 최종 제품의 최소 항복 강도에 어떻게 영향을 미치는지를 정량적으로 보여줍니다. 예를 들어, 금형 온도가 495K에서 505K로 상승하면 최소 항복 강도가 약 144.0 MPa에서 143.9 MPa로 감소하는 것을 예측할 수 있습니다. 이는 공정 엔지니어가 목표로 하는 기계적 물성을 달성하기 위해 금형 온도를 어느 정도의 정밀도로 제어해야 하는지에 대한 공학적 근거를 제공합니다.

Q3: 응고 시뮬레이션에 사용된 수치 모델의 정확성은 어떻게 검증되었습니까?

A3: 개발된 수치 모델(OpenCast 소프트웨어)은 실제 다이캐스팅 문제에 적용하기 전에, Fusegi 등이 발표한 [19] 3차원 자연 대류 벤치마크 문제에 대한 시뮬레이션을 수행하여 검증되었습니다. 64³개의 균일한 육면체 요소로 메쉬를 구성하고, 레일리 수 10⁵ 및 10⁶ 조건에서 온도 및 속도 분포를 계산하여 기존 연구 결과와 비교했습니다. 그 결과, 온도 오차는 0.02% 미만, 속도 오차는 0.8% 미만으로 매우 높은 정확도를 보여 모델의 신뢰성을 확보했습니다.

Q4: 이 연구에서 제안된 프레임워크가 기존의 결정론적 시뮬레이션과 비교하여 갖는 가장 큰 장점은 무엇입니까?

A4: 가장 큰 장점은 공정 변수의 ‘확률적 변동’을 고려하여 제품 품질의 변동 범위를 예측할 수 있다는 점입니다. 기존의 결정론적 시뮬레이션은 단일 값의 결과를 제공하지만, 이 프레임워크는 입력 변수의 불확실성이 출력(예: 항복 강도)의 불확실성으로 어떻게 전파되는지를 정량화합니다. 이를 통해 엔지니어는 평균적인 성능뿐만 아니라 제품의 신뢰성과 강건성(robustness)까지 평가할 수 있습니다.

Q5: 이 프레임워크를 다른 제조 공정에 적용할 수 있습니까?

A5: 네, 가능합니다. 논문의 결론 부분에서 저자들은 이 프레임워크가 다이캐스팅을 예시로 개발되었지만, 다른 제조 공정이나 공학 문제에도 일반화될 수 있다고 명시했습니다. 응고, 미세구조, 물성 예측 모델을 해당 공정에 맞는 모델로 교체하고, 주요 공정 변수를 확률 변수로 정의하면 사출 성형, 용접, 단조 등 다양한 공정의 불확실성을 분석하는 데 동일한 접근 방식을 적용할 수 있습니다.

Conclusion: Paving the Way for Higher Quality and Productivity

이 연구는 기존의 결정론적 시뮬레이션이 가진 한계, 즉 공정 변수의 미세한 변동을 예측하지 못하는 문제를 해결하기 위해 가상 인증 프레임워크를 제시했습니다. 핵심적인 돌파구는 다이캐스팅 불확실성 정량화 기법을 성공적으로 적용하여, 금형 온도나 재료 물성과 같은 입력 변수의 작은 변화가 최종 제품의 기계적 강도에 미치는 영향을 통계적으로 예측한 것입니다.

R&D 및 운영 관점에서 이는 단순히 ‘좋다/나쁘다’를 넘어, ‘어떤 변수가 얼마만큼 변할 때, 제품 품질이 어느 범위 내에서 변동할 것인가’를 정량적으로 파악할 수 있게 되었음을 의미합니다. 이는 보다 과학적인 공정 관리, 정밀한 품질 목표 설정, 그리고 궁극적으로는 더 높은 수율과 신뢰성을 갖춘 제품 생산으로 이어질 수 있습니다.

STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 최선을 다하고 있습니다. 이 논문에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 부품에 어떻게 구현할 수 있는지 논의해 보십시오.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

Copyright Information

• This content is a summary and analysis based on the paper “VIRTUALLY-GUIDED CERTIFICATION WITH UNCERTAINTY QUANTIFICATION APPLIED TO DIE CASTING” by “Shantanu Shahane, et al.”.
• Source: https://asmedigitalcollection.asme.org/VVS/proceedings-pdf/VVS2018/40795/V001T03A003/2792633/v001t03a003-vvs2018-9323.pdf

This material is for informational purposes only. Unauthorized commercial use is prohibited. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 [Jones, Thomas D. A. 외]가 저술하여 [Engineering Science and Technology, an International Journal]에 발표한 논문 “Computational fluid dynamic simulations of solidification for enhancing speed of continuous cast copper”을 기반으로 합니다. STI C&D의 기술 전문가가 분석하고 요약했습니다.

키워드

• Primary Keyword: 구리 연속 주조
• Secondary Keywords: CFD 시뮬레이션, 응고 해석, 수직 상향 연속 주조(VUCC), 열전달 계수, 결정립 구조, 산소 동결 구리(OFCu)

Executive Summary

• The Challenge: 고품질 산소 동결 구리(OFCu)의 생산성을 높이기 위해 주조 속도를 향상시켜야 하지만, 속도 증가는 종종 주조 결함으로 이어져 품질 저하를 유발합니다.
• The Method: Ansys Fluent™를 사용한 전산 유체 역학(CFD) 모델을 구축하여, 8mm 직경 OFCu의 수직 상향 연속 주조(VUCC) 공정 중 발생하는 응고 현상을 시뮬레이션했습니다.
• The Key Breakthrough: 주조 다이 내에 형성되는 0.1mm 에어 갭(Air Gap)이 열전달에 미치는 영향을 정량화했으며, ‘Pushback’ 및 ‘Dwell’과 같은 특정 주조 모션이 응고 전선(Solidification Front)을 안정시켜 더 빠른 속도에서도 우수한 주조 조건을 만들어냄을 규명했습니다.
• The Bottom Line: CFD 시뮬레이션을 통해 주조 공정 변수(속도, 모션)가 최종 제품의 미세구조에 미치는 영향을 예측하고 최적화할 수 있으며, 이는 결함 감소와 생산성 향상에 직접적으로 기여할 수 있습니다.

The Challenge: Why This Research Matters for CFD Professionals

자동차 및 에너지 산업의 핵심 소재인 고순도 산소 동결 구리(OFCu)의 수요는 전기차 및 전자기기 시장의 성장과 함께 급증하고 있습니다. 생산량을 늘리기 위한 가장 직접적인 방법은 연속 주조 공정의 속도를 높이는 것이지만, 이는 기술적 한계에 부딪힙니다. 주조 속도를 무리하게 높이면 더 많은 열이 주조 다이로 유입되어 불완전한 응고를 유발하고, 결국 주조된 봉(rod)이 파손되는 등 심각한 품질 문제를 야기합니다.

업계에서는 이러한 문제를 해결하기 위해 ‘Pushback'(후퇴)이나 ‘Dwell'(정지)과 같은 특수 모션을 주조 사이클에 적용하여 고속에서도 안정적인 주조 조건을 확보하려 노력해왔습니다. 하지만 이러한 모션들이 정확히 어떤 물리적 현상을 통해 주조 품질을 개선하는지에 대한 완전한 이해는 부족했습니다. 특히 응고 과정에서 구리 봉과 다이 사이에 자연적으로 형성되는 미세한 에어 갭(Air Gap)이 열전달에 미치는 영향과, 이것이 전체 응고 거동에 어떻게 작용하는지를 정량적으로 파악하는 것이 중요한 과제였습니다.

The Approach: Unpacking the Methodology

본 연구는 이러한 과제를 해결하기 위해 Ansys Fluent™ v.R1 소프트웨어를 활용하여 2차원 전산 유체 역학(CFD) 시뮬레이션 모델을 구축했습니다. 연구의 핵심은 8mm 직경 OFCu의 수직 상향 연속 주조(VUCC) 공정을 실제와 가깝게 모사하는 것이었습니다.

연구진은 다음과 같은 체계적인 접근법을 사용했습니다.

1. 에어 갭 영향 분석: 먼저, 응고된 구리와 흑연 다이 사이에 형성되는 0.1mm 크기의 에어 갭이 열전달에 미치는 영향을 분석하기 위한 국소 모델(Simulation model 1)을 만들었습니다. 이 시뮬레이션을 통해 에어 갭 존재 시의 유효 열전달 계수(hc) 값을 정확하게 도출했습니다.
2. 전체 공정 시뮬레이션: 다음으로, 도가니(crucible)와 다이를 포함한 전체 주조 장비를 모델링(Simulation model 2)했습니다. 여기에 1단계에서 계산된 열전달 계수 값을 경계 조건으로 입력하여 실제 공정에 가까운 열 분포를 모사했습니다. 이 모델을 이용해 다양한 주조 속도(연속 모션: 0.022 m/s, 0.015 m/s, 0.008 m/s)와 특수 모션(Pushback: 0.06 m/s, Dwell: 0.05 m/s)에 따른 응고 거동 및 온도 변화를 시뮬레이션했습니다.
3. 모델 검증: 시뮬레이션 결과는 기존 문헌에 보고된 실제 주조 공정의 다이 내부 측정 온도 데이터와 비교하여 모델의 신뢰성을 검증했습니다.

시뮬레이션에는 난류 유동을 해석하기 위한 SST(Shear Stress Transport) k-ω 모델과 응고 현상을 모사하기 위한 엔탈피-다공성(Enthalpy-Porosity) 기법이 적용되었습니다.

The Breakthrough: Key Findings & Data

시뮬레이션을 통해 주조 속도 향상의 핵심이 되는 몇 가지 중요한 물리적 현상을 정량적으로 밝혀냈습니다.

Finding 1: 에어 갭(Air Gap)이 열전달을 결정적으로 저해함을 정량화

응고 과정에서 발생하는 구리의 부피 수축은 주조된 봉과 다이 내벽 사이에 약 0.1mm의 미세한 공기층, 즉 에어 갭을 형성합니다. 공기의 열전도율은 구리나 흑연에 비해 현저히 낮기 때문에, 이 에어 갭은 강력한 단열층 역할을 합니다. 시뮬레이션 결과, 에어 갭이 없을 때의 평균 열전달 계수(hc)는 (1.30 ± 0.1) × 10⁵ W/m²K였으나, 100µm(0.1mm) 에어 갭이 존재할 경우 (9.0 ± 0.2) × 10⁴ W/m²K로 크게 감소했습니다 (Figure 4d 참조). 이는 에어 갭이 응고된 구리로부터 다이로의 열 추출을 심각하게 방해하여, 더 많은 열이 주조된 봉 내부에 남아있게 함을 의미합니다. 이 발견은 고속 주조 시 열 관리의 어려움을 설명하는 핵심적인 데이터입니다.

Finding 2: 주조 속도와 응고 전선(Solidification Front)의 명확한 상관관계 규명

주조 속도는 응고가 완료되는 위치, 즉 응고 전선(SF, Solidification Front)의 형태와 위치에 직접적인 영향을 미쳤습니다.

• 고속 주조 (Pushback, 0.06 m/s): 더 빠른 속도로 주조할 경우, 열이 빠져나갈 시간이 부족해 응고가 다이의 더 높은 위치에서 완료되었습니다. 응고가 시작되어 완전히 굳기까지의 구간(mushy zone)의 길이는 약 60mm로 길게 형성되었습니다 (Figure 8d 참조).
• 저속 주조 (Continuous, 0.022 m/s): 상대적으로 느린 속도에서는 열이 충분히 추출되어 다이의 낮은 위치에서 응고가 완료되었으며, 응고 구간의 길이는 약 10mm로 매우 짧았습니다 (Figure 8d 참조).

이 결과는 주조 속도가 빠를수록 응고 전선이 위로 길게 늘어져 불안정해질 가능성이 높으며, 이는 주조 결함 발생의 직접적인 원인이 될 수 있음을 시사합니다. Pushback과 같은 특수 모션은 이러한 응고 전선을 제어하여 고속에서도 안정성을 유지하는 데 기여하는 것으로 분석됩니다.

Practical Implications for R&D and Operations

본 연구 결과는 구리 연속 주조 공정의 최적화를 위한 구체적인 통찰을 제공합니다.

• For Process Engineers: 이 연구는 Pushback 및 Dwell 모션이 단순히 경험적으로 사용되는 기술이 아니라, 응고 전선의 위치와 안정성을 제어하여 고속 주조를 가능하게 하는 핵심 메커니즘임을 보여줍니다. 특정 주조 속도에 맞춰 모션 파라미터(후퇴 속도, 정지 시간 등)를 CFD 시뮬레이션을 통해 최적화함으로써 결함 발생을 최소화하고 생산성을 극대화할 수 있습니다.
• For Quality Control Teams: 논문의 Figure 9는 시뮬레이션으로 예측된 응고 전선의 형상이 실제 주조된 구리 봉의 결정립 성장 방향과 직접적인 관련이 있음을 명확히 보여줍니다. 이는 응고 전선의 형태를 제어하는 것이 최종 제품의 기계적 특성을 결정하는 미세구조를 제어하는 것과 같다는 의미입니다. 이 데이터를 활용하여 새로운 품질 검사 기준을 수립하거나, 미세구조 불량의 원인을 역추적하는 데 활용할 수 있습니다.
• For Design Engineers: 다이의 형상과 냉각 채널의 설계가 다이 내부의 온도 구배를 결정하고, 이는 곧 응고 전선의 위치와 형태에 직접적인 영향을 미칩니다. 본 연구의 결과는 다이 설계 초기 단계에서부터 목표로 하는 미세구조를 얻기 위해 열 추출 효율을 고려하는 것이 얼마나 중요한지를 강조합니다. CFD 해석을 통해 다양한 다이 설계를 사전에 평가하여 개발 시간과 비용을 절감할 수 있습니다.

Paper Details

Computational fluid dynamic simulations of solidification for enhancing speed of continuous cast copper

1. Overview:

• Title: Computational fluid dynamic simulations of solidification for enhancing speed of continuous cast copper
• Author: Jones, Thomas D. A.; Strachan, Richard I.; Mackie, David M.; Cooper, Mervyn; Frame, Brain; Vorstius, Jan B.
• Year of publication: 2021
• Journal/academic society of publication: Engineering Science and Technology, an International Journal
• Keywords: Casting, Copper alloy, Computational fluid dynamic, Simulation, Solidification, Grain structure

2. Abstract:

본 연구에서는 8mm 직경 산소 동결 구리(OFCu)의 수직 상향 연속 주조(VUCC) 시 주조 속도 변경에 따른 현상을 조사하기 위해 Ansys Fluent TM v.R1 내에서 전산 유체 역학(CFD) 모델을 구축했다. 주조 다이 내에 형성된 0.1mm 에어 갭을 통해 전달되는 열의 영향을 시뮬레이션하여, 다이 벽 열전달 계수(hc) 값으로 (9.0 ± 0.2) × 10⁴ W/m²K를 추출했다. 이 hc 값을 사용하여 전체 주조 도가니 및 다이에 대한 시뮬레이션을 다양한 주조 속도 설정(Pushback 모션 0.06 m/s, Dwell 모션 0.05 m/s, 연속 모션 0.022 m/s, 0.015 m/s, 0.008 m/s)에 대해 수행했으며, 주조 다이 내 측정된 열 분포에 대한 문헌 값과 비교하여 검증했다. 8mm OFCu의 가장 빠른 주조 속도를 조사하고, 시뮬레이션된 응고 전선과 측정된 결정립 성장 방향 사이의 경향을 확인하여, Pushback 및 Dwell 모션이 개선된 주조 조건을 만든다는 것을 강조했다. 주조 도가니 내 유체 유동 속도를 조사한 결과, 다이 내 유동(Pushback 주조 시 0.1 ± 0.01 m/s)에 비해 자연 대류로 인한 영향(0.001 ± 0.0005 m/s)이 작아 주조에 미미한 영향을 미치는 것으로 나타났다.

3. Introduction:

순도 99.999% 이상의 구리 금속은 산소 동결 구리(OFCu)로 불리며, 자동차 및 에너지 분야를 중심으로 데이터 통신용 고전도성 재료, 입자 가속기용 초전도 자석, 에나멜선 등 다양한 분야에 응용된다. OFCu 수요 증가는 주로 전기차 및 전자기기 제조의 발전에 의해 주도된다. 연속 주조 공정은 고순도, 고품질 OFCu 금속을 톤 단위로 생산하는 데 사용되는 가장 일반적인 기술로, 주조 속도는 약 100 kg/hour 또는 4 m/min에 달한다. 수직 상향 연속 주조(VUCC)는 OFCu에서 최소한의 결함으로 최고 품질을 산출하며, 평균 주조 속도는 100 kg/hr이다. VUCC는 응고 과정에서 냉각 속도를 높여 결정립 크기를 바람직하게 줄이고 결함을 낮추는 장점이 있다. 금속 응고는 인출 속도, 주조 봉의 직경, 주조 장치로 추출되는 열과 같은 주조 매개변수에 의해 영향을 받는다. VUCC에서 중요한 것은 주조 봉에 Pushback 및 Dwell(움직임의 일시 정지)을 적용하여 주조 조건을 개선하는 것이다. 이러한 모션의 적용은 가장 빠른 속도(3-4 m/min)로 주조할 때 필요하며, 이는 고품질 OFCu를 얻을 수 있는 거의 최대 속도이다. 그러나 이러한 모션의 영향에 대한 완전한 이해는 부족하여 본 연구에서 이를 규명하고자 한다.

4. Summary of the study:

Background of the research topic:

고품질 OFCu 생산을 위한 VUCC 공정에서 생산성 향상을 위해 주조 속도를 높이는 것이 중요하지만, 고속 주조는 열 관리 문제와 응고 불안정성으로 인해 결함을 유발할 수 있다.

Status of previous research:

이전 연구들은 CFD 시뮬레이션을 통해 주조 공정의 열 및 유동 현상을 모델링할 수 있음을 보여주었으나, OFCu의 VUCC 공정에서 Pushback 및 Dwell 모션의 영향을 비교 분석한 연구는 부족했다. 특히 응고 시 발생하는 에어 갭이 열전달에 미치는 영향을 정량적으로 모델에 반영하는 것이 정확도 향상의 핵심 과제로 남아있었다.

Purpose of the study:

본 연구의 목적은 CFD 시뮬레이션을 통해 8mm OFCu의 VUCC 공정에서 주조 속도와 모션(Pushback, Dwell, 연속)이 응고 거동, 열 분포, 최종 미세구조에 미치는 영향을 규명하는 것이다. 이를 통해 주조 속도를 향상시키기 위한 공정 최적화의 과학적 근거를 마련하고자 한다.

Core study:

연구는 세 가지 주요 시뮬레이션을 통해 진행되었다. 첫째, 주조 다이 내 에어 갭의 영향을 분석하여 정확한 열전달 계수(hc)를 도출했다. 둘째, 이 값을 적용하여 전체 도가니 및 다이 시스템에서 다양한 주조 속도와 모션에 따른 응고 현상을 시뮬레이션했다. 셋째, 시뮬레이션 결과를 실제 주조된 구리 봉의 미세구조와 비교하여 모델의 예측 정확도를 검증하고, 주조 조건과 최종 품질 간의 관계를 분석했다.

5. Research Methodology

Research Design:

본 연구는 전산 유체 역학(CFD) 시뮬레이션 방법을 기반으로 설계되었다. 실제 VUCC 공정을 2차원 축대칭 모델로 단순화하여 Ansys Fluent™ 소프트웨어에서 해석을 수행했다. 에어 갭의 영향을 분석하는 국소 모델과 전체 주조 시스템을 해석하는 거시 모델을 연계하여 해석의 정확도를 높였다.

Data Collection and Analysis Methods:

시뮬레이션 데이터는 Ansys Fluent™를 통해 수집되었으며, 온도 분포, 응고율(liquid fraction), 유속 등을 분석했다. 시뮬레이션 모델의 검증을 위해 기존 문헌에 발표된 실제 8mm OFCu 주조 공정의 다이 내부 온도 측정 데이터를 비교 데이터로 사용했다. 또한, 실제 Pushback 및 Dwell 조건에서 주조된 구리 봉 샘플을 채취하여 연마 후 광학 현미경으로 결정립 구조를 관찰하고, 시뮬레이션된 응고 전선 형상과 비교 분석했다.

Research Topics and Scope:

• 연구 주제: 8mm OFCu의 VUCC 공정에서 주조 속도 및 모션이 응고 현상에 미치는 영향 분석
• 연구 범위:
  • 0.1mm 에어 갭이 열전달 계수에 미치는 영향 정량화
  • Pushback, Dwell, 연속 모션 등 다양한 주조 조건 시뮬레이션
  • 주조 조건에 따른 다이 내부 온도 분포 및 응고 전선 위치 변화 분석
  • 시뮬레이션 결과와 실제 주조품의 미세구조 간 상관관계 분석
  • 도가니 및 다이 내부의 유체 유동 패턴 분석

6. Key Results:

Key Results:

• 주조 다이 내 100µm 에어 갭은 평균 열전달 계수(hc)를 (1.30 ± 0.1) × 10⁵ W/m²K에서 (9.0 ± 0.2) × 10⁴ W/m²K로 크게 감소시켜 단열 효과를 발생시킨다.
• 주조 속도가 빠를수록(예: Pushback 0.06 m/s) 응고가 완료되는 구간이 길어지고(약 60mm) 다이의 더 높은 위치에서 응고가 일어나는 반면, 속도가 느릴수록(예: 연속 0.022 m/s) 응고 구간은 짧아진다(약 10mm).
• Pushback 및 Dwell 모션은 응고 전선을 안정시켜 고속 주조 시에도 우수한 주조 조건을 유지하는 데 기여하며, 이는 시뮬레이션된 응고 전선 형상과 실제 결정립 성장 방향의 일치로부터 확인되었다.
• 도가니 내부의 자연 대류에 의한 유속(0.001 ± 0.0005 m/s)은 다이로 유입되는 용탕의 유속(0.1 ± 0.01 m/s)에 비해 매우 작아, 응고 거동에 미치는 영향이 미미하다.

Figure List:

• Fig. 1. a) Schematic of the vertically upwards continuous casting (VUCC) 8 mm process. b) Photograph of setup.
• Fig. 2. a) Schematic of 2-D axisymmetric simulation setup for simulation model 1 & 3 of air gap within the die and b) mesh optimization without airgap. c) The optimal mesh size.
• Fig. 3. Schematic of Simulation model 2, showing a) & b), the crucible and a close-up of the die region, respectively and c) the generated optimal mesh showing close-up within die.
• Fig. 4. Simulation of solidified OFCu metal/graphite die interface, showing with arrows the orientation of the average pushback casting movement and the temperature distributions for A) without air gap and B) with 0.1 mm air gap. C) Horizontal temperature distribution across simulation showing with dashed lines position of air gap. D) Table insert of measured values.
• Fig. 5. 2-D simulated casting for OFCu metal, showing thermal distribution (dark = 1200 °C, light = 20 °C) for pushback setting (0.06 m/s average), over A) the entire simulation, B) within the die and C) upwards motion casting setting (0.022 m/s) within the die.
• Fig. 6. Simulated plots of the thermal distribution from die bottom to top (along y axis) for casting settings a) pushback (0.06 m/s average) and b) continuous at 0.022 m/s, respectively.
• Fig. 7. a) Measured thermal profile inside the graphite die vertically for different casting speeds from [4]. b) Simulated plot of temperature inside graphite vertically for different casting speeds. c) Temperature comparisons between measured and simulated thermal values.
• Fig. 8. 2-D simulated casting for OFCu metal showing mass fraction of solidification for a) pushback (0.06 m/s average), and b) continuous motion (0.022 m/s). Plots of the mass fraction vertically along the die wall (along y axis), at positions c) at the die wall around super cooler and d) in the middle of the die between the two blue arrows.
• Fig. 9. Longitudinal cross sections of 8 mm diameter rod, showing grain growth direction by white arrows, with simulated SF overplaid, for a) pushback cast sample (0.06 m/s, average), showing SF for upwards and pushback cycles as shown as pink (light) and yellow (dark), respectively; and b) & c) dwell setting (0.05 m/s, average) for the pause and upwards movements of the casting cycle, respectively.
• Fig. 10. 2-D simulated casting for OFCu metal (0.06 m/s, pushback), showing copper fluid flow and solid motion. a) Fluid and solid orientations only for the upwards pulse cycle, for the entire crucible, b) the fluid orientation and magnitude within the crucible/die entrance and c) without Marangoni boundary conditions.

7. Conclusion:

본 연구에서는 8mm 직경 OFCu의 VUCC 공정에 대한 2차원 CFD 시뮬레이션을 수행했다. 응고된 금속과 주조 다이 사이에 형성되는 0.1mm 에어 갭이 열전달에 미치는 영향을 모델링한 결과, 에어 갭이 다이로의 열전달을 크게 저해하여 주조된 봉 내부에 더 많은 열이 오래 남게 됨을 확인했다. 이 데이터를 활용하여 도가니/다이 전체 시스템에 대한 CFD 시뮬레이션을 6가지 다른 주조 속도에 대해 수행했다. 시뮬레이션 결과, 주조 속도가 느릴수록 봉에서 더 많은 열이 추출되어 봉 내부 온도가 낮아지고 응고 전선(SF)이 다이의 더 낮은 위치에 형성됨을 보여주었다. 개발된 시뮬레이션 모델은 문헌의 열 측정값과 비교하여 검증되었다. Pushback 및 Dwell 모션을 포함한 두 가지 주조 조건을 비교한 결과, 구리 미세구조가 적용된 모션에 따라 변화하며, 이러한 변화는 응고 전선 위치에 대한 CFD 시뮬레이션 결과로부터 예측 가능함이 확인되었다. 이 모델은 OFCu 및 다른 합금의 주조 거동에 대한 통찰을 제공하고 주조 결과를 예측하여 재료 낭비를 줄이는 데 적용될 수 있다.

8. References:

• [1] Z. Sun, M. Guo, Y. Guo, Research of technological factors on producing oxygen-free copper strip in horizontal continuous casting, Adv. Mater. Res. 538–541 (2012) 1097–1100. doi:10.4028/www.scientific.net/AMR.538-541.1097.
• [2] T. Knych, B. Smyrak, M. Walkowicz, Selected aspects of evolution properties of oxygen free copper for high-advanced electrotechnical application, Prz. Elektrotechniczny 87 (2011) 285–290. https://www2.scopus.com/inward/record.uri?eid=2-s2.0-81155137947&partnerID=40&md5=5d8d2c7111d5a131505019854f9f6113.
• [3] H. Soda, G. Motoyasu, A. McLean, S.D. Bagheri, D.D. Perovic, Continuous casting of unidirectionally solidified copper rod, Int. J. Cast Met. Res. 9 (1) (1996) 37–44, https://doi.org/10.1080/13640461.1996.11819642.
• [4] K. Härkki, J. Miettinen, Mathematical modeling of copper and brass upcasting, Metall. Mater. Trans. B 30 (1) (1999) 75–98, https://doi.org/10.1007/s11663-999-0009-6.
• [5] E.-R. Bagherian, Y. Fan, M. Cooper, B. Frame, A. Abdolvand, Effect of water flow rate, casting speed, alloying elements and pull distance on tensile strength, elongation percentage and microstructure of continuous cast copper alloys, Metall. Res. Technol. 113 (3) (2016) 308, https://doi.org/10.1051/metal/2016006.
• [6] E.-R. Bagherian, Y. Fan, M. Cooper, B. Frame, A. Abdolvand, Effect of melt temperature, cleanout cycle, continuous casting direction (horizontal / vertical) and super-cooler size on tensile strength, elongation percentage and microstructure of continuous cast copper alloys, Metall. Res. Technol. 113 (2016) 502, https://doi.org/10.1051/metal/2016030.
• [7] N. Muramatsu, M. Akaiwa, Characteristics of hypoeutectic Cu-Zr alloy rods manufactured by vertically upwards continuous casting, Mater. Trans. 57 (10) (2016) 1794–1800, https://doi.org/10.2320/matertrans.M2016181.
• [8] E.-R. Bagherian, Y. Fan, M. Cooper, B. Frame, A. Abdolvand, Investigation of the distribution of lead in three different combinations of brass feedstock, Int. Metalcast. 10 (3) (2016) 322–328, https://doi.org/10.1007/s40962-016-0055-1.
• [9] E.R. Bagherian, C. Bell, M. Cooper, Y.C. Fan, B. Frame, M. Rose, Analysis and quantification of grain size of various DHP copper tubes manufacturing processes, Adv. Mater. Res. 856 (2014) 241–245. doi:10.4028/www.scientific.net/AMR.856.241.
• [10] D. Mackie, J.D. Robson, P.J. Withers, M. Turski, Characterisation and modelling of defect formation in direct-chill cast AZ80 alloy, Mater. Charact. 104 (2015) 116–123, https://doi.org/10.1016/j.matchar.2015.03.033. … (and so on for all 46 references)

Expert Q&A: Your Top Questions Answered

Q1: 왜 0.1mm 에어 갭(Air Gap)을 시뮬레이션에 포함하는 것이 그토록 중요했습니까?

A1: 에어 갭은 응고 시 구리의 부피 수축으로 인해 필연적으로 발생하는 물리 현상입니다. 공기는 구리나 흑연 다이에 비해 열전도율이 매우 낮아 강력한 단열층으로 작용합니다. 이 연구는 에어 갭을 모델에 포함함으로써 실제 공정에서 발생하는 열 저항을 정확히 반영할 수 있었고, 그 결과로 도출된 열전달 계수((9.0 ± 0.2) × 10⁴ W/m²K)는 전체 주조 시스템 시뮬레이션의 정확도를 획기적으로 높이는 핵심 요소가 되었습니다. 에어 갭을 무시했다면 열 추출이 과대평가되어 실제와 다른 응고 거동을 예측했을 것입니다.

Q2: 시뮬레이션 결과, Pushback 모션이 연속 모션에 비해 더 나은 주조 조건을 만드는 이유는 무엇입니까?

A2: Figure 8과 9에서 볼 수 있듯이, Pushback 모션은 더 빠른 평균 주조 속도를 가능하게 하면서도 응고 전선(SF)을 안정적으로 유지하는 데 기여합니다. Pushback 사이클 동안의 일시적인 후퇴 및 정지 동작은 응고 계면에 열이 재분배될 시간을 주어, 급격한 온도 변화를 완화하고 더 균일한 응고를 유도합니다. 그 결과, Figure 9a에서처럼 응고 전선이 안정적인 형태를 유지하며 성장하여, 최종적으로 바람직한 기둥 모양의 결정립(columnar grain) 구조를 형성하게 됩니다. 이는 고속에서도 주조 결함을 억제하고 품질을 유지할 수 있는 핵심 메커니즘입니다.

Q3: 이 연구 결과가 8mm OFCu 봉 외에 다른 직경이나 합금에도 적용될 수 있습니까?

A3: 네, 적용 가능합니다. 이 연구에서 개발된 CFD 모델링 접근법은 재료 물성치와 공정 변수를 변경하여 다른 조건에도 적용할 수 있습니다. 예를 들어, 다른 직경의 구리 봉이나 구리 합금을 주조할 경우, 해당 재료의 열물성 데이터(열전도율, 밀도, 잠열 등)와 변경된 주조 속도, 다이 형상 등을 모델에 입력하면 유사한 방식으로 응고 거동을 예측하고 공정을 최적화할 수 있습니다. 따라서 이 연구는 특정 사례를 넘어 연속 주조 공정 전반에 적용할 수 있는 강력한 분석 도구를 제시했다는 점에서 의의가 큽니다.

Q4: 그림 8d에서 더 빠른 주조 속도(0.06 m/s)에서 응고율이 일시적으로 감소하는 현상이 나타나는 이유는 무엇입니까?

A4: 이는 다이 내 액상 금속의 무작위적인 대류 운동 때문일 가능성이 높습니다. 고속으로 주조될 때 다이 내부의 용탕 유동은 더 복잡하고 불안정해집니다. 국소적인 대류 흐름이 일시적으로 더 뜨거운 용탕을 응고 계면으로 가져와 온도를 높이면, 그 지점의 응고가 잠시 지연되거나 이미 형성된 고상 일부가 다시 녹으면서 응고율이 국부적으로 감소할 수 있습니다. 이러한 현상은 고속 주조 시 응고 전선이 불안정해질 수 있음을 보여주는 증거이기도 합니다.

Q5: 도가니 내부의 자연 대류가 응고에 미치는 영향이 미미하다고 결론 내린 근거는 무엇입니까?

A5: Figure 10의 시뮬레이션 결과에 따르면, 도가니 내부에서 자연 대류로 인해 발생하는 유체의 평균 속도는 약 0.001 ± 0.0005 m/s에 불과했습니다. 반면, 다이 내부로 빨려 들어가는 용탕의 속도는 Pushback 주조 시 약 0.1 ± 0.01 m/s로 측정되었습니다. 이처럼 다이 내부의 강제 대류 속도가 도가니의 자연 대류 속도보다 약 100배 더 빠르기 때문에, 도가니에서 발생하는 약한 유동은 다이 내부의 응고 현상에 거의 영향을 미치지 못한다고 판단할 수 있습니다.

Conclusion: Paving the Way for Higher Quality and Productivity

본 연구는 CFD 시뮬레이션이 구리 연속 주조 공정의 복잡한 물리 현상을 얼마나 정밀하게 예측하고 분석할 수 있는지를 명확하게 보여줍니다. 특히, 눈에 보이지 않는 0.1mm 에어 갭의 영향을 정량화하고, Pushback과 같은 특수 주조 모션이 어떻게 응고 전선을 제어하여 고속에서도 품질을 유지하는지를 과학적으로 규명했습니다. 이는 경험에 의존하던 공정 최적화를 데이터 기반의 엔지니어링 영역으로 끌어올린 중요한 성과입니다.

이러한 해석 기술을 활용하면, R&D 및 운영팀은 실제 생산 라인에서 값비싼 시행착오를 거치지 않고도 다양한 공정 변수의 영향을 사전에 파악하고 최적의 주조 조건을 찾아낼 수 있습니다. 결과적으로 이는 주조 결함을 줄이고, 생산 속도를 한계까지 끌어올려 궁극적으로 더 높은 품질과 생산성을 달성하는 길을 열어줄 것입니다.

“STI C&D에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.”

• 연락처 : 02-2026-0442
• 이메일 : flow3d@stikorea.co.kr

Copyright Information

• This content is a summary and analysis based on the paper “Computational fluid dynamic simulations of solidification for enhancing speed of continuous cast copper” by “Jones, T. D. A., et al.”.
• Source: https://doi.org/10.1016/j.jestch.2020.12.009

This material is for informational purposes only. Unauthorized commercial use is prohibited. Copyright © 2025 STI C&D. All rights reserved.

이 기술 요약은 Sebastian Kohlstädt, Michael Vynnycky, Stephan Goeke가 Metals (2021)에 발표한 논문 “On the CFD Modelling of Slamming of the Metal Melt in High-Pressure Die Casting Involving Lost Cores”를 기반으로 합니다. STI C&D에서 기술 전문가를 위해 분석하고 요약했습니다.

키워드

• Primary Keyword: HPDC 슬래밍
• Secondary Keywords: 로스트 코어, CFD 모델링, OpenFOAM, 2상 유동, 고압 다이캐스팅, 솔트 코어

Executive Summary

• 도전 과제: 로스트 코어를 사용하는 고압 다이캐스팅(HPDC) 공정에서, 용탕이 코어에 처음 충돌할 때 발생하는 ‘슬래밍(Slamming)’ 현상은 코어를 파손시킬 수 있는 급격한 힘의 스파이크를 유발합니다.
• 연구 방법: OpenFOAM 소프트웨어를 사용하여 용탕-공기 2상 유동을 모델링하고, 압축성 가스 및 RANS 난류 모델을 적용하여 솔트 코어에 가해지는 힘을 CFD로 시뮬레이션했습니다.
• 핵심 발견: 슬래밍으로 인한 힘의 스파이크는 매우 크지만, 그 지속 시간이 극도로 짧아 코어의 변위와 응력이 파손 한계 이하에 머무릅니다.
• 핵심 결론: 균열이 없는 이상적인 코어 조건에서는 슬래밍 현상 자체가 코어 파손의 결정적인 원인이 아닐 수 있으며, 오히려 지속적으로 가해지는 정상 상태의 힘이 더 중요할 수 있습니다.

도전 과제: 이 연구가 CFD 전문가에게 중요한 이유

고압 다이캐스팅(HPDC)은 변속기 하우징, 크랭크 케이스 등 자동차 부품을 대량 생산하는 핵심 공정입니다. 그러나 내부가 비어 있거나 복잡한 형상을 만들기 위해 로스트 코어(Lost Core), 특히 솔트 코어(Salt Core)를 사용하는 기술은 아직 양산 단계에 이르지 못했습니다. 가장 큰 기술적 장벽 중 하나는 초고속(50-100 m/s)으로 주입되는 용탕이 코어에 처음 부딪힐 때 발생하는 ‘슬래밍’ 현상입니다.

이 순간적인 충격은 코어에 엄청난 부하를 가하며, 이전의 모든 시뮬레이션에서 이 힘의 ‘스파이크’가 예측되었습니다. 만약 이 힘이 코어의 파손 한계를 초과한다면, 복잡한 내부 채널을 가진 고품질의 주조품 생산은 불가능해집니다. 따라서 엔지니어들은 이 슬래밍 힘을 정확히 예측하고, 그 힘이 실제로 코어 파손에 얼마나 치명적인지 규명해야 하는 과제를 안고 있었습니다. 이 연구는 바로 이 문제에 대한 해답을 찾기 위해 시작되었습니다.

접근 방식: 연구 방법론 분석

본 연구는 슬래밍 현상을 정밀하게 분석하기 위해 오픈 소스 CFD 소프트웨어 패키지인 OpenFOAM을 활용했습니다. 연구의 핵심 방법론은 다음과 같습니다.

• 유동 모델: 용탕(알루미늄 합금 AlSi9Cu3)과 공기의 상호작용은 오일러리안 VOF(Volume-of-Fluid) 접근법을 사용하여 2상 유동으로 모델링했습니다. 공기는 압축성 이상 기체로 처리되었습니다.
• 난류 모델: RANS(Reynolds-averaged Navier Stokes) 접근법을 기반으로 Menter SST k-ω 모델과 RNG k-ε 모델을 비교 평가했으며, 최종적으로 Menter SST k-ω 모델을 채택했습니다.
• 지배 방정식: 질량, 운동량, 에너지 보존에 대한 방정식을 풀었으며, 표면 장력은 CSF(Continuum Surface Force) 방법으로 모델링했습니다.
• 해석 조건: 2D 채널 내에 위치한 솔트 코어에 용탕이 충돌하는 상황을 시뮬레이션했습니다. 메시(Mesh) 해상도와 시간 간격(Courant Number)이 슬래밍 힘 계산에 미치는 영향을 집중적으로 분석했습니다.

이러한 접근을 통해 연구진은 기존 산업계에서 일반적으로 사용되던 시뮬레이션 조건보다 훨씬 더 정밀한 조건에서 슬래밍 현상을 포착하고, 그 물리적 의미를 심도 있게 분석할 수 있었습니다.

핵심 발견: 주요 결과 및 데이터

결과 1: 슬래밍 힘의 정확한 예측을 위해서는 극도로 미세한 메시가 필수적이다

연구의 가장 중요한 발견 중 하나는 슬래밍 힘을 정확하게 예측하기 위해서는 메시 해상도가 결정적이라는 것입니다. Figure 8은 메시 크기에 따른 무차원 슬래밍 계수(Cs)의 변화를 보여줍니다.

• 산업계에서 일반적으로 사용되는 1mm 수준의 거친 메시에서는 슬래밍 계수가 von Karman의 이론적 하한값(π, 약 3.14)에 훨씬 못 미치는 값으로 과소 예측되었습니다.
• 메시 크기를 0.3mm 이하로 매우 미세하게 만들어야만 비로소 계산된 슬래밍 계수가 이론값에 근접하는 것을 확인할 수 있었습니다. 이는 기존 시뮬레이션이 슬래밍의 최대 충격력을 제대로 포착하지 못했을 가능성을 시사합니다. 논문에 따르면, 정확한 결과를 얻기 위해서는 메시 해상도를 40배나 높여야 했습니다.

Figure 8: 메시 연구 결과. 거친 메시(우측)는 슬래밍 팩터(Cs)의 최대값을 과소평가하는 반면, 미세한 메시(좌측)로 갈수록 von Karman 및 Wagner 모델의 이론값에 수렴한다.

결과 2: 충격력 스파이크는 크지만, 코어 파손의 결정적 원인은 아닐 수 있다

시뮬레이션 결과, 슬래밍으로 인한 힘의 스파이크는 매우 높았지만(약 1400 N), 그 힘이 작용하는 시간은 3.5 x 10⁻⁶ 초로 극도로 짧았습니다. 연구진은 이 결과를 바탕으로 간단한 연속체 역학 계산을 수행했습니다.

• 이 짧은 시간 동안 코어가 받는 충격량(force-time integral)은 매우 작았습니다.
• 계산된 코어의 최대 변위는 약 4 x 10⁻⁷ m에 불과했습니다.
• 이 정도의 미미한 변위는 코어 내부에 발생하는 굽힘 응력(bending stress)이 재료의 파괴 한계에 도달하기에 턱없이 부족한 수준입니다.

결론적으로, 균열이 없는 이상적인 솔트 코어라면, 슬래밍으로 인한 순간적인 힘의 스파이크는 코어를 파손시키기보다는 단지 진동시키는 데 그칠 가능성이 높습니다. 오히려 코어 파손은 슬래밍 이후에 지속적으로 가해지는 정상 상태(steady-state)의 유체력이 더 결정적인 역할을 할 수 있습니다.

R&D 및 운영을 위한 실질적 시사점

• 공정 엔지니어: 이 연구는 슬래밍의 초기 충격 완화보다, 이후 지속되는 정상 상태의 압력을 견딜 수 있도록 코어의 지지 구조를 설계하는 것이 더 중요할 수 있음을 시사합니다. 또한, 코어를 예열하면 취성 파괴보다는 연성 파괴를 유도하여 공정 안정성을 높일 수 있습니다.
• 품질 관리팀: 연구의 결론은 ‘균열 없는(crack-free)’ 코어를 전제로 합니다. 이는 솔트 코어 생산 시 미세 균열을 최소화하는 것이 매우 중요하며, 비파괴 검사 등 강화된 품질 검사 기준이 필요함을 의미합니다.
• 설계 엔지니어: CFD 시뮬레이션 시, 슬래밍 하중을 정확히 평가하려면 코어 주변의 메시를 매우 조밀하게 구성해야 합니다. 하지만 이 연구 결과에 따르면, 초기 충격 스파이크 값 자체보다는 코어가 정상 상태의 흐름 압력을 견딜 수 있는지에 대한 구조적 안정성을 검토하는 것이 더 실질적인 설계 기준이 될 수 있습니다.

논문 상세 정보

On the CFD Modelling of Slamming of the Metal Melt in High-Pressure Die Casting Involving Lost Cores

1. 개요:

• 제목: On the CFD Modelling of Slamming of the Metal Melt in High-Pressure Die Casting Involving Lost Cores
• 저자: Sebastian Kohlstädt, Michael Vynnycky, and Stephan Goeke
• 발행 연도: 2021
• 발행 학술지: Metals
• 키워드: compressible two-phase flow; slamming; OpenFOAM; high-pressure die casting; lost salt cores; solid continuum mechanics

2. 초록:

본 논문은 OpenFOAM 소프트웨어 패키지를 이용한 전산 유체 역학(CFD)을 사용하여, 고압 다이캐스팅(HPDC)의 다이 충전 단계에서 유입되는 용탕의 충격에 노출될 때 솔트 코어에 가해지는 힘을 연구하며, 특히 ‘슬래밍’으로 알려진 첫 충격 순간에 초점을 맞춘다. 용탕-공기 시스템은 오일러리안 VOF(volume-of-fluid) 접근법을 통해 모델링되며, 공기는 압축성 이상 기체로 취급된다. 난류는 RANS(Reynolds-averaged Navier Stokes) 접근법을 통해 처리된다. RNG k-ε 모델과 Menter SST k-ω 모델을 모두 평가하였으며, 후자를 최종적으로 배치 계산에 채택하였다. 메시 독립성을 확립하기 위한 Courant 수의 영향에 대한 연구는, 기존 분석 모델 및 경험적 측정과 관련하여 코어에 대한 슬래밍 효과를 적절하게 포착하기 위해 이전에 HPDC 시뮬레이션에 사용되었던 것보다 더 미세한 메시와 더 작은 시간 단계가 필요함을 나타낸다. 두 번째 단계로, 스파이크 형태의 형태와 작은 힘-시간 적분값을 가진 이 힘이 코어에 충돌할 때 어떤 반응이 예상되어야 하는지 논의한다. 힘의 스파이크로 인한 코어의 변위는 매우 작아서, 힘의 값은 높지만 코어 내부의 굽힘 응력은 파괴의 임계 한계 이하로 유지된다는 것이 발견되었다. 따라서 균일하고 균열이 없는 재료 조건을 가정할 때, 힘의 스파이크는 파손에 결정적이지 않다고 결론 내릴 수 있다.

3. 서론:

고압 다이캐스팅(HPDC)은 자동 변속기 하우징, 크랭크 케이스, 기어 박스 부품과 같은 대량 생산 및 저비용 자동차 부품 제조에 중요한 공정이다. 일반적으로 알루미늄이나 마그네슘과 같은 액체 금속이 복잡한 게이트 및 러너 시스템을 통해 다이 안으로 주입되며, 인게이트에서 50~100 m/s의 속도와 100 MPa에 달하는 높은 압력 하에서 이루어진다. 경제적 관점에서 HPDC는 막대한 초기 설비 및 금형 투자가 필요하지만, 추가 생산 단위당 증분 비용이 낮아 생산량이 증가함에 따라 규모의 경제가 매우 잘 작동한다. 반면, 이로 인해 설계 엔지니어는 금형 및 설비에 예산을 투자하기 전에 공정과 제조 부품의 실행 가능성을 확신해야 하는 과제를 안게 된다. 현재까지의 기술적 제약 중 하나는, 다른 주조 기술에서는 수십 년간 로스트 코어를 사용해왔음에도 불구하고, HPDC를 통해 로스트 코어로 형성된 내부 중공 형상이나 언더컷이 있는 부품의 양산이 없다는 점이다.

4. 연구 요약:

연구 주제의 배경:

HPDC 공정에서 로스트 코어(특히 솔트 코어)를 사용하여 복잡한 내부 형상을 가진 부품을 생산하려는 시도가 있었으나, 용탕 충전 시 발생하는 높은 압력과 속도로 인해 코어가 파손될 위험이 있어 양산에 어려움을 겪고 있다.

이전 연구 현황:

이전 연구들에서도 CFD를 통해 용탕이 코어를 지나가는 유동을 분석했지만, 대부분 일반적인 유동 특성에 초점을 맞췄다. 용탕이 코어에 처음 충돌하는 ‘슬래밍’ 현상과 그로 인한 순간적인 최대 부하가 코어의 파손에 미치는 영향에 대한 심층적인 분석은 부족했다.

연구 목적:

본 연구의 목적은 CFD 시뮬레이션을 통해 HPDC 공정 중 발생하는 슬래밍 현상이 솔트 코어에 가하는 힘, 특히 그 최대값을 정밀하게 예측하고, 이 힘이 실제로 코어의 파손을 유발하는지 여부를 규명하는 것이다.

핵심 연구:

1. OpenFOAM을 이용한 2상(용탕-공기) 압축성 유동 모델 개발.
2. 메시 해상도, 시간 간격(Courant number), 난류 모델이 슬래밍 힘 계산에 미치는 영향 분석.
3. 계산된 슬래밍 힘을 기존의 이론적/경험적 모델과 비교 검증.
4. 슬래밍으로 인한 힘의 스파이크가 코어에 미치는 구조적 영향(변위, 응력)을 연속체 역학을 통해 분석하여 파손 가능성 평가.

5. 연구 방법론:

연구 설계:

본 연구는 2차원 채널 내에 고정된 솔트 코어에 용탕이 충돌하는 상황을 가정하여 전산 유체 역학(CFD) 시뮬레이션을 수행하는 방식으로 설계되었다. 메시 해상도와 Courant 수를 체계적으로 변경하며 슬래밍 힘의 변화를 관찰하여 수치 해석의 신뢰도를 확보하고자 했다.

데이터 수집 및 분석 방법:

• 소프트웨어: OpenFOAM (compressibleInterFoam 솔버)
• 물리 모델: VOF(Volume-of-Fluid) 2상 유동 모델, RANS 난류 모델(Menter SST k-ω), 압축성 이상 기체(공기)
• 데이터 분석: 코어 표면에 작용하는 압력 및 전단력을 적분하여 시간에 따른 총 힘을 계산하고, 이를 무차원화된 ‘슬래밍 계수(Cs)’로 변환하여 분석했다. 계산된 힘과 지속 시간을 바탕으로 코어의 변위와 응력을 추정했다.

연구 주제 및 범위:

연구는 HPDC 공정의 충전 단계에서 발생하는 슬래밍 현상에 국한된다. 응고 과정이나 후처리 공정은 고려하지 않는다. 대상은 솔트 코어이며, 코어 재료는 균일하고 내부에 결함이 없는 것으로 가정한다.

6. 주요 결과:

주요 결과:

• HPDC 슬래밍 현상 시뮬레이션에서 정확한 최대 충격력을 예측하기 위해서는 기존 산업계에서 사용하던 것보다 훨씬 미세한 메시(0.3mm 이하)와 작은 시간 단계가 필요하다.
• 거친 메시를 사용한 시뮬레이션은 슬래밍으로 인한 최대 힘을 심각하게 과소평가한다.
• 난류 모델(Menter SST k-ω vs. RNG k-ε)의 선택은 슬래밍 힘 계산 결과에 미미한 영향을 미쳤다. 이는 슬래밍이 점성보다는 압력과 운동량에 의해 지배되는 현상임을 시사한다.
• 슬래밍으로 인해 발생하는 힘의 스파이크는 매우 크지만, 지속 시간이 극도로 짧아 코어에 가하는 총 충격량은 작다.
• 계산 결과, 이 힘으로 인한 코어의 변위와 내부 응력은 재료의 파괴 한계보다 훨씬 낮아, 슬래밍 현상 자체가 균열 없는 코어의 파손 원인이 될 가능성은 낮다.

Figure 목록:

• Figure 1. The geometry for investigating the slamming on a salt core in channel; all dimensions are in mm.
• Figure 2. An example of a computational grid created with the utilities blockMesh and mirrorMesh for a mesh spacing of 2 mm.
• Figure 3. Pressure Implicit with Splitting of Operator (PISO) algorithm before and after the adjustments.
• Figure 4. Normalized forces on the core for a mesh resolution of 0.3 mm.
• Figure 5. Phase distribution of the melt at impact and immediately afterwards.
• Figure 6. The pressure field at the times of impact and immediately afterwards.
• Figure 7. The velocity magnitude field at the times of impact and immediately afterwards.
• Figure 8. Mesh study of the slamming factor in comparison with the models by von Karman and Wagner.
• Figure 9. Comparison of the computed result with reference studies in previously published articles; mesh cell spacing 0.025 mm.
• Figure 10. Influence of the selected turbulence model on the computed result for the slamming factor.
• Figure 11. Results of the time step size ∆t study: (a) Results for Fstat; (b) Results for Fmax.

7. 결론:

본 연구는 산업계에서 일반적으로 사용되는 1mm 수준의 메시를 사용한 최신 CFD 시뮬레이션이 슬래밍 현상을 과소평가하고 있음을 보여주었다. 계산된 슬래밍 계수가 기존의 분석 결과 및 경험적 측정 범위 내에 들기 위해서는 메시 해상도를 40배까지 높여야 함을 발견했다. 흥미롭게도 시간 간격(∆t)의 감소는 슬래밍 계수 값의 정확한 추정에 큰 영향을 미치지 않았다. 또한, 난류의 기여도는 무시할 수 있는 수준이었다.

슬래밍으로 인한 힘의 스파이크가 코어 내부에 실제로 어떤 반응을 일으키는지에 대한 과학적 논쟁이 남아있다. 본 논문은 간단한 연속체 역학 공식을 기반으로, 주어진 파라미터에 대해 코어가 단지 약간 변위되고 진동만 할 뿐이며, 따라서 슬래밍 충격은 파손에 결정적이지 않은 현상이라는 추론을 제공한다. 그러나 이 결과는 제시된 경계 조건에 대해서만 유효하며, 다른 설정으로 이전할 때는 신중해야 함을 독자에게 상기시킨다.

8. 참고 문헌:

1. Kaufmann, H.; Uggowitzer, P. Metallurgy and Processing of High-Integrity Light Metal Pressure Castings; Schiele & Schön: 2007.
2. Nogowizin, B. Theorie und Praxis des Druckgusses, 1st ed.; Schiele und Schoen: Berlin, Germany, 2010.
3. Brunnhuber, E. Praxis der Druckgussfertigung; Schiele und Schoen: Berlin, Germany, 1991.
4. Campbell, J. Complete Casting Handbook: Metal Casting Processes, Metallurgy, Techniques and Design; Elsevier Science:Amsterdam, The Netherlands, 2015.
5. Jelínek, P.; Adámková, E. Lost cores for high-pressure die casting. Arch. Foundry Eng. 2014, 14, 101–104.
6. Kohlstädt, S.; Rabus, U.; Goeke, S.; Kuckenburg, S. Verfahren zur Herstellung Eines Metallischen Druckgussbauteils Unter Verwendung Eines Salzkerns Mit Integrierter Stützstruktur und Hiermit Hergestelltes Druckgussbauteil. DE Patent DE10,201,410,221,359A1, 21 April 2016.
7. Schneider, T.; Kohlstädt, S.; Rabus, U. Gehäuse mit Druckgussbauteil zur Anordnung Eines Elektrischen Fahrmotors in Einem Kraftfahrzeug und Verfahren zur Herstellung eines Druckgussbauteils. DE Patent DE10,201,410,221,358A1, 21 April 2016.
8. Graf, E.; Soell, G. Vorrichtung zur Herstellung Eines Druckgussbauteils mit Einem Kern und Einem Einlegeteil. DE Patent DE10,145,876A1, 10 April 2003 .
9. Graf, E.; Söll, G. Verfahren zur Herstellung Eines Zylinderkurbelgehäuses. WO Patent WO2,008,138,508A1, 20 November 2008.
10. Fuchs, B.; Körner, C. Mesh resolution consideration for the viability prediction of lost salt cores in the high pressure die casting process. Prog. Comp. Fluid Dyn. 2014, 14, 24–30.
11. Kohlstädt, S.; Vynnycky, M.; Gebauer-Teichmann, A. Experimental and numerical CHT-investigations of cooling structures formed by lost cores in cast housings for optimal heat transfer. Heat Mass Trans. 2018, 54, 3445-3459.
12. Kohlstädt, S.; Vynnycky, M.; Neubauer, A.; Gebauer-Teichmann, A. Comparative RANS turbulence modelling of lost salt core viability in high pressure die casting. Prog. Comp. Fluid Dyn. 2019, 19, 316–327.
13. Fuchs, B.; Eibisch, H.; Körner, C. Core viability simulation for salt core technology in high-pressure die casting. Int. J. Met. 2013, 7, 39-45.
14. Fuchs, V. Numerische Modellierung von Fluid-Struktur-Wechselwirkungen an wellenbeaufschlagten Strukturen; Kassel University Press GmbH: Kassel, Germany, 2014.
15. Campbell, T.; Weynberg, P. Measurement of parameters affecting slamming – final report. Technical Report, Wolfson Marine Craft Unit Report No. 440– University of Southampton: Southampton, UK, 1980.
16. Abrate, S. Hull slamming. Appl. Mech. Rev. 2011, 64, 060803. … (and so on for all 59 references)

전문가 Q&A: 자주 묻는 질문

Q1: 왜 상용 소프트웨어가 아닌 OpenFOAM을 사용하여 연구를 진행했나요?

A1: 논문에 따르면 OpenFOAM은 특정 응용 분야에 맞춰 확장성이 뛰어나고, 라이선스 제약 없이 대규모 CPU 코어를 활용한 병렬 컴퓨팅이 가능하여 복잡한 3D 산업 모델로의 확장을 염두에 두고 선택되었습니다. 또한, 연구에 사용된 Menter SST k-ω 난류 모델의 구현이 견고하고 실험 데이터와 잘 일치하는 것으로 이전 연구들에서 입증되었기 때문입니다.

Q2: Figure 8을 보면, 거친 메시에서 정상 상태 힘(steady-state force)은 오히려 과대평가되는 경향이 보입니다. 이는 해석의 신뢰도에 문제가 되지 않나요?

A2: 논문에서는 이를 ‘설계 엔지니어 관점에서 유익한 결과’라고 언급합니다. 정상 상태의 힘이 실제보다 더 높게 예측된다는 것은, 해당 시뮬레이션 결과를 바탕으로 코어를 설계할 때 자연스럽게 안전 계수(safety factor)가 포함된다는 의미입니다. 따라서 코어의 지속적인 하중에 대한 안정성을 보수적으로 평가하게 되어 오히려 긍정적인 측면이 있습니다.

Q3: 이 연구는 2D 모델을 기반으로 ‘슬래밍 스파이크는 치명적이지 않다’고 결론 내렸습니다. 이 결과를 복잡한 3D 형상에 그대로 적용할 수 있을까요?

A3: 아니요, 그럴 수 없습니다. 논문은 결과가 주어진 형상과 경계 조건에 매우 의존적임을 명확히 경고하고 있습니다. 따라서 이 연구의 정량적인 결과를 다른 형상에 직접 적용하기보다는, 이 연구에서 제시한 ‘미세 메시를 사용한 CFD 접근법’을 실제 분석하고자 하는 3D 형상에 적용하여 자체적으로 평가하는 것이 매우 중요합니다.

Q4: Figure 10에서 난류 모델의 영향이 거의 없었는데, 굳이 다른 모델을 비교 분석한 이유는 무엇인가요?

A4: 이는 슬래밍 현상을 지배하는 물리적 요인을 명확히 하기 위함이었습니다. 분석 결과, 난류 모델의 영향이 미미하다는 것은 이 현상이 난류에 의한 점성 효과보다는 용탕의 운동량과 그로 인한 압력에 의해 주로 결정된다는 것을 입증합니다. 이는 향후 유사한 시뮬레이션을 수행할 때, 난류 모델 선택에 대한 고민을 줄여주고 해석의 주요 변수에 더 집중할 수 있게 해주는 중요한 정보입니다.

Q5: 연구의 핵심 결론은 코어가 ‘균열이 없다’는 가정에 기반합니다. 실제 공정에서 생산된 솔트 코어에 이 가정이 얼마나 현실적인가요?

A5: 논문은 이 가정이 현실과 다를 수 있는 잠재적 한계임을 인정합니다. 실제 세라믹 재질의 솔트 코어에는 미세한 균열이 존재할 가능성이 높습니다. 이 때문에 논문은 향후 연구 과제로, 슬래밍 충격을 견디기 위해 잠재적 균열이 어느 정도 길이까지 허용될 수 있는지를 평가하는 모델 개발을 제안합니다. 이는 현실에서는 코어의 초기 품질이 파손 여부에 큰 영향을 미칠 수 있음을 시사합니다.

결론: 더 높은 품질과 생산성을 향한 길

요약하자면, 본 연구는 고압 다이캐스팅 공정에서 발생하는 HPDC 슬래밍 현상이 기존의 일반적인 CFD 시뮬레이션에서 과소평가되어 왔음을 명확히 보여주었습니다. 정확한 충격력 예측을 위해서는 매우 정밀한 메시 설정이 필수적입니다.

하지만 더 중요한 통찰은, 그토록 우려했던 순간적인 힘의 스파이크가 실제로는 코어 파손의 주범이 아닐 수 있다는 것입니다. R&D 및 운영 관점에서 이는 문제의 초점을 ‘순간 충격 완화’에서 ‘지속 하중에 대한 코어의 구조적 안정성 확보’와 ‘코어 자체의 품질 관리’로 전환해야 함을 의미합니다.

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

• 연락처 : 02-2026-0450
• 이메일 : flow3d@stikorea.co.kr

저작권 정보

• 이 콘텐츠는 Sebastian Kohlstädt 등의 논문 “On the CFD Modelling of Slamming of the Metal Melt in High-Pressure Die Casting Involving Lost Cores”를 기반으로 요약 및 분석한 자료입니다.
• 출처: https://doi.org/10.3390/met11010078

본 자료는 정보 제공 목적으로 제작되었으며, 무단 상업적 사용을 금합니다. Copyright © 2025 STI C&D. All rights reserved.