서론: 파라미터 최적화가 제자리걸음하는 이유
금속 레이저 분말 베드 용융(LPBF) 및 고분자 소결 공정의 현장 셋업에서 가장 흔한 상황은 다음과 같다. 엔지니어가 레이저 출력을 250 W에서 280 W로 조정했더니 상대 밀도가 0.3% 상승하여 계속 출력을 높인다. 300 W에서 키홀 현상이 발생하면 다시 출력을 낮춘다. 한 사이클 동안 10개 이상의 시편을 만들지만, 결과로 얻는 것은 단편적인 경험 곡선일 뿐이다. 장비를 바꾸거나 분말 배치를 교체하면 다시 처음부터 시작해야 한다. 문제는 엔지니어의 노력 부족이 아니라, 단일 인자 순차 조정법 자체가 상호작용을 처리할 수 없다는 데 있다. 출력과 스캔 속도는 함께 체적 에너지 밀도를 결정하고, 분말층 두께는 다시 용융풀의 방열 조건을 바꾼다. 세 인자가 서로 결합되어 있으므로 한 번에 하나의 변수만 변경해서는 전체 양상을 결코 볼 수 없다.
실험계획법(Design of Experiments, DoE)은 이러한 다인자 결합 문제를 해결하기 위한 엔지니어링 경로를 제공한다. 최소한의 시편 수로 인자의 주효과, 상호작용 효과 및 반응 곡면을 한 번에 추정하고, 신뢰구간을 포함한 허용 공정 윈도우를 제시할 수 있다. 본문에서는 3D 프린팅 공정 개발을 대상으로 직교 스크리닝, 반응표면 모델링 및 윈도우 검증까지의 완전하고 실행 가능한 방법을 소개한다.
1. 먼저 세 가지 실험을 구분해야 한다: 스크리닝, 모델링 및 검증
공정 개발은 한 번의 실험으로 완료되지 않으며, 일반적으로 세 단계로 진행된다. 각 단계는 목적, 인자 수 및 시편 수가 다르다.
첫 번째 단계는 스크리닝 실험으로, 6~10개의 후보 인자 중 실제로 유의한 2~4개를 선별하는 것이 목적이다. 대표적인 후보 인자는 레이저 출력, 스캔 속도, 해치 간격, 층 두께, 예열 온도, 스캔 전략(스트라이프 폭), 층간 회전각 및 보호가스 유량이다. 스크리닝 단계에서는 부분 요인 설계 또는 Plackett-Burman 설계를 사용하며, 8~12회 시험으로 완료할 수 있다. 응답은 상대 밀도와 표면 조도라는 두 가지 거시 지표만 평가한다.
두 번째 단계는 모델링 실험으로, 선별된 2~3개의 핵심 인자에 대해 반응표면 모델을 구축하는 것이 목적이다. 일반적으로 중심합성설계(CCD) 또는 Box-Behnken 설계를 사용한다. 3인자 CCD는 축점 6개, 큐브점 8개 및 중심점 반복 6개로 구성되어 총 20회 시험이 필요하다. 응답은 상대 밀도, 인장강도, 치수 편차 및 표면 조도 Ra의 네 항목으로 확대된다.
세 번째 단계는 검증 실험으로, 모델이 예측한 최적 영역에서 3~5개의 확인점을 선정하여 반복 출력하고, 모델 예측값과 실측값의 편차가 ±5% 이내인지 확인한다. 검증을 통과한 후 윈도우 경계를 10% 축소하여 양산 공정 카드의 파라미터 범위로 설정한다.
2. 직교 실험: 최소 시편으로 최대 공간을 포괄하기
인자 수가 많지 않고 각 인자의 수준 수가 같을 때 직교표는 가장 경제적인 선택이다. 316L 스테인리스강 LPBF를 예로 들어 레이저 출력 P, 스캔 속도 v, 해치 간격 h의 세 인자를 각각 3수준으로 검토한다고 하자. 완전 요인 설계에는 27회 시험이 필요하지만, L9(3^4) 직교표를 사용하면 9회만 필요하다.
구체적인 설계는 다음과 같다. P는 200/240/280 W, v는 800/1000/1200 mm/s, h는 0.09/0.11/0.13 mm로 설정한다. L9에 따라 시험을 배치한 뒤 각 수준 조합에서 10×10×10 mm 큐브 시편 한 세트를 출력하고, 아르키메데스법으로 상대 밀도를 측정하면서 표면 조도 Ra도 기록한다.
데이터 처리는 범위 분석으로 빠르게 순위를 정할 수 있다. 각 인자의 수준별 평균 응답을 계산하고, 범위 R = max(평균 응답) − min(평균 응답)으로 구한다. R이 클수록 해당 인자의 영향이 더 유의하다. R(P)=2.1%, R(v)=1.4%, R(h)=0.6%가 얻어졌다면 중요도 순서는 출력 > 속도 > 간격이다. 추가로 분산분석(ANOVA)을 수행하여 F 검정으로 유의성을 판단한다. F 값이 임계값보다 크고 p 값이 0.05 미만인 인자가 다음 모델링 단계에 포함된다.
직교 실험의 한계는 인자 간 상호작용이 없다고 가정한다는 점이다. 출력과 속도 간 강한 결합이 의심되는 경우, 즉 체적 에너지 밀도 E = P/(v·h·t) 자체가 곱셈 관계인 경우에는 직교표를 계속 적용하기보다 모델링 단계에서 반응표면 설계로 전환해야 한다.
3. 반응표면법: 이산 점에서 연속 곡면으로
반응표면법(RSM)의 핵심은 2차 다항식으로 응답과 인자의 관계를 적합하는 것이다.
Y = β0 + Σβi·xi + Σβii·xi² + Σβij·xi·xj + ε
여기서 교차항 βij는 상호작용 효과를 정량적으로 표현한다. CCD 설계를 예로 들면, 3인자 설계의 핵심 파라미터는 축 방향 거리 α이다. 회전 가능 설계에서는 α = 2^(k/4) = 2^(3/4) ≈ 1.682를 사용하며, 여기서 k는 인자 수이다. 축점은 코딩 공간에서 ±1.682에 위치하고, 큐브점은 ±1에 위치하며, 중심점은 순수 오차 추정을 위해 6회 반복한다.
적합 완료 후에는 반드시 세 가지 검증을 수행해야 한다. 첫째, 적합 부족(lack-of-fit) 검정에서 p 값이 0.05보다 크면 모델 형식과 데이터 사이에 유의한 불일치가 없음을 의미한다. 둘째, 결정계수 R²와 수정 R²를 확인하며, 두 값의 차이는 0.2 미만이어야 한다. 그렇지 않으면 과적합 가능성이 있다. 셋째, 잔차의 정규성과 등분산성을 점검한다. 실제 측정에서는 R²가 0.95에 달하지만 수정 R²가 0.86에 불과한 경우가 흔하다. 이때는 유의하지 않은 3차항을 제거한 뒤 다시 적합해야 한다.
다중 응답 최적화는 3D 프린팅에서 특히 어려운 부분이다. 상대 밀도는 최대화해야 하고, 조도는 최소화해야 하며, 치수 편차는 0에 가까워야 한다. 이 세 목표는 대개 서로 충돌한다. 엔지니어링에서는 만족도 함수(desirability function)를 사용하여 각 응답을 0~1 구간으로 매핑한 다음, 기하평균으로 종합 만족도 D를 계산하고 D가 최대가 되는 지점에서 파라미터 조합을 결정한다. 가중치는 고객의 인수 우선순위에 따라 배정한다. 예를 들어 항공 부품은 상대 밀도 가중치를 0.5, 조도를 0.3, 치수 편차를 0.2로 설정할 수 있다.
4. 공정 윈도우 설정: 최적점을 허용 구간으로 전환하기
최적점을 찾는 것은 첫 단계일 뿐이며, 양산에는 하나의 구간이 필요하다. 실용적인 방법은 중첩 등고선도(overlay contour plot)를 사용하는 것이다. 각 응답의 합격 임계값을 등고선으로 표시하고(상대 밀도 ≥ 99.5%, Ra ≤ 12 μm, 치수 편차 ≤ ±0.1 mm), 모든 등고선이 겹치는 영역을 공정 윈도우로 설정한다.
윈도우가 결정된 후에는 반드시 안전 여유를 확보해야 한다. 권장 방법은 윈도우 경계를 안쪽으로 10% 축소하는 것이다. 예를 들어 모델이 출력의 허용 구간을 235~265 W로 제시한다면, 공정 카드에는 238~262 W로 기재한다. 분말 배치 변동, 보호가스 순도 변화 및 광학 시스템 출력 감쇠가 실제 공정점을 서서히 이동시킬 수 있기 때문이다. 여유를 확보하면 공정능력지수 Cpk를 1.0에서 1.33 이상으로 높일 수 있다.
윈도우 검증은 3개 배치 확인 방식으로 수행한다. 배치당 5개 시편을 제작하고, 각각 윈도우 중심, 상한 경계 안쪽 및 하한 경계 안쪽에서 한 세트씩 출력한다. 연속 3개 배치가 모두 합격하고 배치 간 평균값 이동이 응답 표준편차의 1/3 미만이면 윈도우가 안정적이라고 판정할 수 있다. 검증 데이터는 이후 추적성의 기준선으로 활용할 수 있도록 공정 이력대장에도 기록한다.
5. 일반적인 오류와 점검표
오류 1: 시편 수만을 비용으로 간주하는 것. 실제로 9회의 직교 실험과 20회의 CCD 실험에 드는 총 시편 비용은 30회의 무계획 단일 인자 순차 조정보다 훨씬 낮으며, 재사용 가능한 모델도 얻을 수 있다.
오류 2: 블록 효과를 무시하는 것. 20회의 CCD 실험을 3일에 걸쳐 수행하면 장비 상태, 분말 배치, 환경 온도 및 습도의 변화가 잔차에 혼입될 수 있다. 올바른 방법은 날짜를 블록 인자로 설계에 포함하거나, 설계 내에 블록별 중심점 반복을 배치하여 드리프트를 모니터링하는 것이다.
오류 3: 상대 밀도만 확인하는 것. 상대 밀도가 99.5%에 도달한 공정점이라도 기계적 특성이 합격하지 않을 수 있으며, 미융합 결함은 밀도 검사에서 과소평가될 수 있다. 최소한 인장강도와 연신율을 동시에 측정하고, 필요 시 미세조직 금속학 검사와 산업용 CT 샘플 검사를 추가하는 것이 권장된다.
실행 점검표: 인자 수준이 장비가 안정적으로 작동할 수 있는 범위를 포괄하는가? 순수 오차 추정을 위한 중심점 반복을 배치했는가? 적합 부족 검정을 수행했는가? 다중 응답 가중치가 고객 인수 기준과 일치하는가? 윈도우에 10% 여유를 두었는가? 검증 배치가 윈도우 경계를 포함하는가?
결론
실험계획법은 3D 프린팅 공정 개발을 경험적 시행착오에서 정량적이고 재현 가능한 엔지니어링 활동으로 전환한다. 직교 실험은 9~12회 시험으로 인자 스크리닝을 완료하고, 반응표면법은 약 20회 시험으로 상호작용을 포함하는 2차 모델을 구축하며, 만족도 함수는 다목적 충돌을 해결한다. 중첩 등고선은 안전 여유가 포함된 허용 공정 윈도우를 제공한다. Blueprint 3D와 같이 설계 도면부터 실물 납품까지 전 수명주기를 관리해야 하는 플랫폼에서 DoE의 가치는 단순한 파라미터 조정에 그치지 않는다. 공정 지식을 검색, 재사용 및 검증 가능한 파라미터 모델로 축적하여, 모든 공정 개발이 다시 원점에서 시작하는 일이 아니라 다음 개발의 출발점이 되게 한다.