Fig. 7. Temperature Distribution of Specimen in the vicinity of Failure site by Computer Solidification Simulation.

이 기술 요약은 한국주조공학회지(2014)에 게재된 김헌주 저자의 “AC2B 알루미늄 주조합금의 정량적 응고균열 강도 평가” 논문을 바탕으로, (주)에스티아이씨앤디의 기술 전문가에 의해 분석 및 요약되었습니다.

Keywords

  • Primary Keyword: 정량적 응고균열 강도 평가
  • Secondary Keywords: 응고균열, 열간균열, AC2B, 알루미늄 주조, CAE, 응고해석, 주조성 평가

Executive Summary

  • The Challenge: CAE를 이용한 주조품의 응고균열 예측 시뮬레이션의 정확도를 높이기 위해서는 고액공존상태(mushy zone)에 있는 합금의 정량적 응력-변형률 데이터가 필수적이지만, 기존의 상대적 평가 방법으로는 이 데이터를 얻기 어렵습니다.
  • The Method: AC2B 알루미늄 합금이 응고하는 과정에서 발생하는 응력-변형률 관계를 정량적으로 측정할 수 있는 인장형 응고균열 평가 장치를 자체 개발하여 물리적 시험을 수행했습니다.
  • The Key Breakthrough: 물리적 시험과 컴퓨터 시뮬레이션(파단 위치의 온도 예측) 및 열역학 계산(고상율 계산)을 결합하여, 특정 조건에서 합금의 정확한 응고균열 강도 데이터를 확보하는 표준화된 평가 절차를 확립했습니다.
  • The Bottom Line: 본 연구는 주조품의 응고균열을 예측하고 방지하기 위한 CAE 시뮬레이션의 신뢰도를 획기적으로 향상시키는 데 필요한 핵심 물성 데이터를 확보하는 강력한 방법론을 제공합니다.

The Challenge: Why This Research Matters for CFD Professionals

자동차 산업을 중심으로 경량화 요구가 증가하면서 알루미늄 합금의 사용이 확대되고 있습니다. 특히 복잡한 형상의 부품을 고압금형주조(HPDC)법으로 생산할 때, 응고 과정에서 발생하는 수축 응력으로 인한 응고균열(hot tearing)은 고질적인 문제입니다.

엔지니어들은 제품 개발 초기 단계에서 CAE 열응력 해석을 통해 이러한 문제를 예측하고 해결하고자 합니다. 하지만 시뮬레이션의 정확도는 입력되는 재료 물성 데이터의 신뢰성에 크게 좌우됩니다. 기존의 응고균열 평가법(예: CRC 몰드)은 여러 합금의 균열 민감도를 ‘상대적으로’ 비교할 수는 있었지만, 고액공존상태에서 균열이 발생하는 한계 응력과 변형률을 ‘정량적’ 수치로 제공하지 못했습니다. 이 데이터의 부재는 CAE 해석 결과와 실제 현상 간의 괴리를 유발하여, 결국 많은 시간과 비용이 소요되는 시행착오를 야기하는 원인이었습니다.

The Approach: Unpacking the Methodology

이 연구는 AC2B 알루미늄 주조합금의 응고균열 특성을 정량적으로 평가하기 위해 새로운 접근법을 채택했습니다. 연구의 핵심은 자체 제작한 인장형 응고균열 평가 장치입니다(그림 1).

  • 시험 장비 및 시편: H형상의 시험편(그림 2)이 응고되는 동안, 장비의 한쪽 끝에 장착된 볼트를 0.75kW 모터로 일정한 변형 속도(60 mm/min)로 인장시킵니다. 이때 발생하는 하중과 변위는 로드셀(Load Cell)과 LVDT를 통해 실시간으로 측정되어 응력-변형률 곡선을 얻습니다.
  • 재료 및 공정 조건: 상용 AC2B 알루미늄 합금(Table 1)을 750°C에서 용해하고, Al-5wt%Ti-1wt%B로 결정립 미세화 처리 및 질소 가스로 탈가스 처리를 수행했습니다. 용탕은 710°C의 온도로 150°C로 예열된 금형에 주입되었습니다(Table 3).
  • 데이터 확보를 위한 하이브리드 접근법: 이 연구의 독창성은 물리적 실험과 컴퓨터 시뮬레이션을 결합한 데 있습니다.
    1. 물리적 시험: 위에서 설명한 장치를 이용해 응고 중인 시편에 인장력을 가하여 파단 시점의 최대 응력과 변형률을 직접 측정합니다.
    2. 컴퓨터 응고해석 (AnyCasting®): 실제 파단이 일어나는 위치에 온도계를 설치하면 응력 집중점으로 작용하여 측정의 정확성을 해치므로, 컴퓨터 시뮬레이션을 통해 비침습적으로 파단 부위의 정확한 온도 변화를 예측했습니다(Table 2).
    3. 열역학 계산 (ThermoCalc®): 시뮬레이션으로 예측된 파단 시점의 온도(600°C)를 바탕으로, Scheil 응고 모델을 적용하여 해당 온도의 정확한 고상율(Solid Fraction)을 계산했습니다.
Fig. 7. Temperature Distribution of Specimen in the vicinity of Failure site by Computer Solidification Simulation.
Fig. 7. Temperature Distribution of Specimen in the vicinity of Failure site by Computer Solidification Simulation.

이러한 통합적 접근을 통해 특정 공정 조건 하에서 응고균열이 발생하는 순간의 물리적 상태(온도, 고상율, 응력, 변형률)를 정확하게 규명할 수 있었습니다.

The Breakthrough: Key Findings & Data

본 연구를 통해 도출된 핵심 성과는 정량적 평가 절차의 확립과 이를 통해 얻어진 구체적인 데이터입니다.

Finding 1: 정량적 응고균열 강도 평가 절차의 표준화

가장 중요한 성과는 신뢰성 있는 데이터를 얻기 위한 일련의 평가 절차를 확립한 것입니다(그림 6). 이 절차는 [시편 준비 → 열간균열 시험 수행 → 파단 위치 측정 → 컴퓨터 시뮬레이션을 통한 파단점 온도 예측 → 열역학 계산을 통한 고상율 산출 → 합금 특성과의 상관관계 분석]으로 구성됩니다. 이는 기존의 경험적, 상대적 평가 방식에서 벗어나, CAE 해석에 직접 활용 가능한 물리적 데이터를 생성하는 과학적이고 체계적인 방법론을 제시합니다.

Finding 2: 특정 조건 하 AC2B 합금의 구체적인 응고균열 강도 데이터 확보

제안된 절차를 통해 AC2B 합금에 대한 구체적인 데이터를 확보했습니다. – 금형 예열온도 150°C 조건에서, 시편의 파단은 시뮬레이션을 통해 600°C에서 발생하는 것으로 예측되었습니다(그림 7). – ThermoCalc® 계산 결과, 이 온도에서 고상율(Solid Fraction)은 29.3%였습니다(Table 4). – 이 시점에서 측정된 최대 응력(응고균열 강도)은 12.9 kgf/cm²였으며, 파단 변형률은 2.5%로 나타났습니다(Table 5).

이러한 정량적 데이터는 특정 합금이 응고 과정의 어느 단계에서, 얼마만큼의 기계적 강도를 가지는지 명확히 보여주며, 이는 CAE 시뮬레이션의 정확도를 비약적으로 향상시키는 핵심 정보가 됩니다.

Practical Implications for R&D and Operations

본 연구 결과는 주조 공정, 품질 관리, 제품 설계 등 다양한 분야의 전문가들에게 실질적인 시사점을 제공합니다.

  • For Process Engineers: 이 연구는 금형 예열 온도와 같은 공정 변수가 파단 시점의 고상율에 직접적인 영향을 미치며, 이는 곧 응고균열 강도와 직결됨을 보여줍니다. 즉, 금형 온도를 조절하여 균열에 취약한 고상율 구간을 제어함으로써 균열 발생을 억제하는 공정 최적화가 가능합니다.
  • For Quality Control Teams: 논문의 Table 5에 제시된 데이터(최대 응력 12.9 kgf/cm², 파단 변형률 2.5%)는 AC2B 합금이 응고 중 가장 취약한 상태일 때의 기계적 물성에 대한 정량적 기준을 제공합니다. 이는 보다 정교한 품질 검사 기준을 수립하는 데 활용될 수 있습니다.
  • For Design Engineers: 이 연구에서 제시된 방법론은 CAE 열응력 해석의 정확도를 높이는 데 필요한 핵심 응력-변형률 데이터를 제공합니다. 이를 통해 제품 설계 초기 단계에서 응고균열 발생 가능성을 훨씬 더 정확하게 예측할 수 있으며, 이는 값비싼 시제품 제작과 금형 수정을 최소화하여 개발 기간과 비용을 단축시키는 효과를 가져옵니다.

Paper Details

AC2B 알루미늄 주조합금의 정량적 응고균열 강도 평가

1. Overview:

  • Title: AC2B 알루미늄 주조합금의 정량적 응고균열 강도 평가 (Quantitative Evaluation of Solidification Crack Strength of AC2B Aluminum Casting Alloy)
  • Author: 김헌주 (Heon-Joo Kim)
  • Year of publication: 2014
  • Journal/academic society of publication: 한국주조공학회지 (Journal of Korea Foundry Society)
  • Keywords: Solidification crack test, Solidification crack strength, Solidification crack test procedure

2. Abstract:

CAE 해석을 통한 열응력의 수치 해석은 주조품의 개발 단계에서 응고균열 문제를 예측하고 처리하는 효과적인 방법이 될 수 있다. 이 경우 정량적인 응력-변형률 데이터가 필요하다. 본 연구에서는 응력-변형률 관계를 정량적으로 측정할 수 있는 인장형 응고균열 시험 장치를 개발하고 시험 절차를 확립하였다. 다음 시험 절차로부터 얻어진 응고균열 강도는 고체의 결정립 크기, 결정립 형태, 고체 결정립의 분포 등과 같은 응고 합금의 열-소성 특성에 대한 영향 인자 측면에서 평가하는 데 활용될 수 있다. 제안된 시험 절차는 다음과 같다: 응고 시뮬레이션을 통한 응고균열 시편의 파단 부위 온도 예측, Scheil 조건 하의 응고 열역학적 해로부터 파단 부위의 고상 분율 계산.

3. Introduction:

자동차 산업에서 연비 향상 및 환경 규제 대응을 위해 차체 경량화가 중요해지면서, 철계 부품을 알루미늄 경량 소재로 대체하려는 연구가 활발히 진행 중이다. 주조법으로 부품을 생산할 때에는 합금의 주조성(수축성, 유동성, 응고균열 민감성) 평가가 필수적이다. 특히 고압금형주조(HPDC) 시 Al-Mg계, Al-Cu-Si계 합금에서 응고균열 문제가 빈번히 발생한다. 기존의 응고균열 평가법(열간균열 시험)은 주로 CRC 몰드와 같이 상대적인 민감도를 비교하는 방식이며, 고액공존상태에서 균열이 발생하는 한계 응력에 대한 정량적 데이터를 제공하지 못한다. 실제 제품 개발에서 CAE를 이용한 열응력 해석의 정확도를 높이려면 정량적인 응력-변형률 곡선 자료가 필요하다. 따라서 본 연구는 주조용 알루미늄 합금의 응고균열 민감성을 정량적으로 평가하기 위한 장치를 개발하고 실험 방법론을 정립하는 것을 목표로 한다.

4. Summary of the study:

Background of the research topic:

자동차 경량화 추세에 따라 알루미늄 주조 부품의 수요가 증가하고 있으며, 이와 함께 생산성 높은 주조 공법에서 발생하는 응고균열 문제가 중요한 기술적 과제로 부상했다.

Status of previous research:

기존의 응고균열 평가는 주로 여러 합금의 균열 민감도를 상대적으로 비교하는 정성적 방법에 머물러 있었다. 이는 CAE 해석에 필요한 정량적인 기계적 물성 데이터를 제공하는 데 한계가 있었다.

Purpose of the study:

주조용 AC2B 알루미늄 합금을 대상으로, 응고 과정에서 발생하는 균열 강도를 정량적으로 평가할 수 있는 시험 장치를 개발하고, 신뢰성 있는 데이터 확보를 위한 표준화된 시험 절차를 확립하고자 한다.

Core study:

자체 개발한 인장형 응고균열 시험 장치를 이용해 AC2B 합금의 응고 중 기계적 거동을 측정하고, 컴퓨터 응고해석 및 열역학 계산을 결합하여 파단 시점의 온도와 고상율에 따른 응고균열 강도를 정량적으로 규명하였다.

5. Research Methodology

Research Design:

실제 주조 공정과 유사한 조건에서 응고 중인 시편에 직접 인장력을 가하여 파단 시점의 기계적 특성을 측정하는 실험적 연구 설계를 채택했다. 실험 데이터의 한계를 보완하기 위해 컴퓨터 시뮬레이션과 열역학 계산을 병행하는 하이브리드 접근법을 사용했다.

Data Collection and Analysis Methods:

  • 데이터 수집: 로드셀과 LVDT를 사용하여 응고 중인 시편의 하중 및 변위 데이터를 실시간으로 수집했다. K-type 열전대를 이용해 금형 내 특정 위치의 냉각 곡선을 측정했다.
  • 데이터 분석: 수집된 하중-변위 데이터로 응력-변형률 곡선을 작성하여 최대 응력(균열 강도)과 파단 변형률을 도출했다. 주조해석 S/W (AnyCasting®)를 사용해 파단 위치의 온도를 예측하고, 열역학 S/W (ThermoCalc®)를 이용해 해당 온도의 고상율을 계산하여 실험 결과와 연관지어 분석했다.

Research Topics and Scope:

본 연구는 상용 AC2B 알루미늄 주조합금을 대상으로 하며, 금형 예열온도 150°C, 주입온도 710°C의 특정 주조 조건 하에서 인장형 응고균열 강도를 정량적으로 평가하는 방법론을 개발하고 검증하는 데 초점을 맞추었다.

6. Key Results:

Key Results:

  • 응고 중인 합금의 응고균열 강도를 정량적으로 측정할 수 있는 H자형 시험편 기반의 인장형 시험 장치 및 평가 방법론을 성공적으로 개발 및 확립하였다.
  • 금형 예열온도 150°C 조건에서 AC2B 합금의 파단은 600°C에서 발생하며, 이때 고상율은 29.3%로 계산되었다.
  • 해당 조건에서 측정된 AC2B 합금의 최대 응고균열 강도는 12.9 kgf/cm², 파단 변형률은 2.5%로 정량화되었다.
  • 컴퓨터 응고해석을 통해 예측한 냉각 속도와 실제 측정한 냉각 속도가 잘 일치함을 확인하여(그림 9), 시뮬레이션을 통한 파단점 온도 예측의 유효성을 검증하였다.
Fig. 3. Photograph of Solidification Cracking Test Apparatus.
Fig. 3. Photograph of Solidification Cracking Test Apparatus.

Figure List:

  • Fig. 1. Constituents of Solidification Cracking Test Apparatus.
  • Fig. 2. Schematic Mold Drawing of Solidification Cracking Test.
  • Fig. 3. Photograph of Solidification Cracking Test Apparatus.
  • Fig. 4. Comparison of Cooling Curves measured from different location of Mold Preheated to 150°C.
  • Fig. 5. Typical Stress-Strain Curve and Cooling Curve of Solidification Cracking Test.
  • Fig. 6. Test Procedure for Solidification Cracking Test.
  • Fig. 7. Temperature Distribution of Specimen in the vicinity of Failure site by Computer Solidification Simulation.
  • Fig. 8. Comparison of Simulated Cooling Curves to Measured Cooling Curves.
  • Fig. 9. Comparison of Measured Cooling Rates to Simulated Cooling Rates.

7. Conclusion:

  1. 본 연구에서 개발한 H자형 시험편을 사용한 인장형 응고균열 강도 실험 방법으로 응고 중인 합금의 응고균열 강도를 정량적으로 측정할 수 있었다.
  2. 자체 개발한 인장형 응고균열 장치는 마찰 문제를 보완한 H자형 시험편 금형, 인장응력을 가하는 구동부(모터), 응력/변형을 측정하는 데이터 획득부(로드셀, LVDT) 및 데이터 출력부로 구성되어 있다.
  3. 응고해석을 통해 시험편 균열 지점의 온도를 예측하고, 열역학 계산 프로그램을 통해 해당 온도의 고상율을 환산하는 과정을 통해, 실제 실험으로 측정된 응고균열 강도값에 해당하는 응고 합금의 고상율을 유추 가능하였다.
  4. H자형 시험편으로 인장형 응고균열 강도를 평가한 본 실험 방법은 응고과정에서 형성된 고상의 결정립 크기, 형상, 분포 등 금속 조직 특성 변화가 응고균열 강도에 미치는 영향을 평가할 수 있다.

8. References:

  1. Makhlouf M. and Makhlouf D. Apelian, ACRC, Paper No. DEFC07-99ID13716, “Casting characteristics of Aluminum die casting alloys”, (2002) 1-46.
  2. Jorstad J. L, Die Casting Engineer, “Future technology in Aluminum die casting”, 9 (2006) 18-25.
  3. Koch H, Hielsher. U, Sternau H. and Franke A. J., TMS, “Silafont-36, the new low iron high-pressure die-casting alloy”, (1995) 1011-1018.
  4. Watanabe T, Kimura R, Nakazawa T, Chiba H, Tanaka S, Ueki T, Toriyama T. and Yoshida M, J. of Japan Institute of Light Metals, “Development of the theories of solidification cracking and measuring methods for mechanical properties of solid-liquid coexistence alloys”, 58 (2008) 395-405.
  5. Clyne T. W. and Davies G. J, Proc. Conf. on Solidification and Casting of Metals, “Solidification and casting of metals”, (1979) 275.
  6. Jorstad J, Wayne L. and Zalensas D. L, American Foundrymen’s Society Inc., “Aluminum Casting Technology”, 2nd Edi., (1997) 21-26.
  7. Viano D, StJohn D, Grandfield J. and Caceres C, Light Metals 2005, “Hot Tearing in Aluminium-Copper Alloys” (2005) 1069-1073.
  8. Bishop H. F. and Pellini W. S, AFS Trans., “Investigation of Metallurgical and Mechanical Effects in Development of Hot Tearing”, 65 (1957) 247-258.
  9. Novikov I. I, Goryachelomkost tsvetnykh metallov I splavov, Moscow, Nauka “Investigation of Metallurgical and Mechanical Effects in Development of Hot Tearing”, (1966) 299.
  10. Novikov I. I, Matveeva KT Zavod Lab, “Hot Shortness of Non-Ferrous Metals and Alloys”, 11 (1957) 1362.
  11. Davidson C, Viano D, Lu I. and StJohn D, International Journal of Cast Material Research, “Observation of crack initiation during hot tearing”, 19 (2006) 59-65.
  12. Mitchell J. B, Cockcroft S. I, Viano D, Davidson C. and StJohn D, Metallurgical and Materials Transaction A, “Determination of Strain during Hot Tearing by Image Correlation”, 38A (2007) 2503-2511.
  13. Instone S, Ph. D. Thesis, “The Effect of Alloy Composition and Microstructure on Hot Cracking of Vertical Direct Chill Cast Aluminum Alloy Billet”, University of Queensland, (1999).
  14. Instone S, StJohn D. and Grandfield J, International Journal of Cast Material Research, “New apparatus for characterising tensile strength development and hot cracking in the mushy zone”, 12 (2000) 441-456.
  15. Dahle A. K, Instone S. and Sumitomo T, Metallurgical and Materials Transaction A, “Relationship between Tensile and Shear Strengths of the Mushy Zone in Solidifying Aluminum Alloys”, 34A (2003) 105-113.

Expert Q&A: Your Top Questions Answered

Q1: 파단 위치의 온도를 측정하기 위해 열전대 대신 컴퓨터 시뮬레이션을 사용한 이유는 무엇인가요?

A1: 실제 시험편의 파단이 예상되는 지점에 열전대를 직접 설치할 경우, 열전대 자체가 응력 집중점으로 작용하여 자연스러운 파단 위치와 거동에 영향을 미칠 수 있습니다. 이는 측정 결과의 신뢰도를 떨어뜨릴 수 있으므로, 본 연구에서는 비침습적인 방법인 컴퓨터 응고해석을 통해 실제 파단 위치의 온도를 정확하게 예측함으로써 시험의 객관성과 정확성을 확보했습니다.

Q2: 논문에서 파단 시점의 고상율을 29.3%로 계산했는데, 이 값은 사용된 계산 모델(Scheil)에 얼마나 민감한가요?

A2: 본 연구에서는 고상 내 확산은 없고 액상 내 완전 혼합을 가정하는 Scheil 모델을 사용했습니다. 이는 실제 주조 공정의 빠른 냉각 속도를 잘 모사하기 때문입니다. 만약 고상과 액상 모두에서 완전한 확산을 가정하는 평형(Equilibrium) 모델을 사용했다면, 같은 온도(600°C)에서 고상율은 더 낮게 계산되었을 것입니다. 따라서 정확한 고상율을 얻기 위해서는 실제 공정 조건을 가장 잘 반영하는 열역학 모델을 선택하는 것이 매우 중요합니다.

Q3: 응고해석 모델(AnyCasting®)의 온도 예측 신뢰성은 어떻게 검증되었나요?

A3: 시뮬레이션 모델의 신뢰성을 확보하기 위해, 연구진은 별도의 검증용 시험편 내 4개의 특정 위치에 열전대를 설치하여 실제 냉각 곡선을 측정했습니다. 그 후, 시뮬레이션에서 동일한 위치의 온도 변화를 예측하여 두 결과를 비교했습니다. 그림 8과 9에서 볼 수 있듯이, 측정된 냉각 곡선 및 냉각 속도와 시뮬레이션 예측치가 매우 잘 일치하여 모델의 정확성이 검증되었습니다.

Q4: Table 5에서 최대 응력이 12.9 kgf/cm²로 나타났습니다. 만약 금형 예열온도를 더 낮추면 이 값은 어떻게 변할 것으로 예상할 수 있나요?

A4: 논문에 따르면 응고균열 강도는 고상율에 비례하는 경향을 보입니다. 금형 예열온도를 낮추면 냉각 속도가 빨라져 동일한 시간 내에 더 많은 고상이 형성되므로, 인장력이 가해지는 시점의 고상율이 더 높아질 것입니다. 따라서 금형 예열온도를 150°C보다 낮게 설정하면 측정되는 최대 응고균열 강도는 12.9 kgf/cm²보다 더 높게 나타날 것으로 예상할 수 있습니다.

Q5: 결론에서 이 방법론이 미세조직의 영향을 평가할 수 있다고 언급했는데, 구체적으로 어떤 방식인가요?

A5: 예를 들어, 결정립 미세화제의 첨가량을 조절하여 결정립 크기를 변화시키거나, 냉각 속도를 제어하여 덴드라이트 간격을 바꾸는 등 미세조직 인자를 체계적으로 변경할 수 있습니다. 이렇게 미세조직이 다른 시편들을 대상으로 본 연구에서 제안한 정량적 응고균열 강도 평가를 각각 수행하면, 특정 미세조직 특성(예: 미세한 등축정)과 측정된 균열 강도 값의 변화를 직접적으로 연관시킬 수 있습니다. 이를 통해 합금 설계 및 공정 개발을 보다 과학적으로 수행할 수 있습니다.

Conclusion: Paving the Way for Higher Quality and Productivity

응고균열은 복잡한 알루미늄 주조품의 품질과 생산성을 저해하는 오랜 난제였습니다. 본 연구는 AC2B 합금의 정량적 응고균열 강도 평가를 위한 혁신적인 방법론을 제시함으로써 이 문제에 대한 중요한 돌파구를 마련했습니다. 실험과 시뮬레이션을 결합하여 확보한 신뢰도 높은 물성 데이터는 CAE 해석의 정확도를 극대화하여, 개발 초기 단계에서부터 균열 발생을 예측하고 최적의 설계를 가능하게 합니다. 이는 결국 시제품 제작 비용 절감과 개발 기간 단축으로 이어져 기업의 경쟁력을 강화할 것입니다.

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

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

Copyright Information

  • This content is a summary and analysis based on the paper “AC2B 알루미늄 주조합금의 정량적 응고균열 강도 평가” by “김헌주”.
  • Source: http://dx.doi.org/10.7777/jkfs.2014.34.4.136

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

PDF View