최근 급변하는 사회와 시장에 유연하게 대응하기 위해서 제조업은 제품공급 중심에서 다양한 수요와 서비스 중심으로 빠르게 변화하고 있다. 이에 따라, 제조업 경쟁력을 강화하기 위한 제조산업 디지털 전환이 적극적으로 추진되고 있으며 스마트공장을 그 핵심적인 예로 들 수 있다. 스마트공장을 구현하기 위해서는 가공시스템과 가공공정의 정밀한 물리 모델을 포함한 실제 공정과 동기화된 디지털트윈의 개발이 필수적으로 요구된다. 이를 기반으로 장비와 공정의 상태 진단 및 가공물 품질 예측이 가능하며 나아가 장비와 공정의 최적화 및 제품의 생산 관리가 가능해지기 때문이다. 절삭가공 중 발생하는 절삭력은 장비 구조물의 휨이나 진동을 발생시켜 가공물의 형상정밀도와 표면거칠기를 떨어뜨리며 제어 특성에도 영향을 미치기 때문에 절삭가공 공정의 예측과 최적화, 장비구조 및 제어시스템의 고강성 설계를 위해서 가장 우선 고려하는 요소이다. 따라서, 가공시 절삭력을 실시간으로 정확하게 예측하는 모델은 공작기계 디지털트윈의 핵심요소라고 할 수 있다.

밀링가공에서의 절삭력을 예측하기 위해 많은 연구가 진행되었다. 역학적 절삭 모델(Mechanistic Cutting Force Model) 연구에서부터 단일 경로 또는 코너[1] 등 국소영역에서의 절삭 해석에 대한 연구들이 수행되었다. 반면, 최근 활발히 연구 및 개발되고 있는 가상공작기계(Virtual Machine Tool)나 디지털트윈 기술은 가공품질 예측과 장비 및 공정 최적화를 위해 활용되는데[2,3], 이를 위해서는 NC 밀링가공 시뮬레이션이 요구된다. NC 밀링가공 시뮬레이션은 G코드로 입력되는 복잡한 공구경로에 의해 가공이 진행됨에 따라 변화하는 가공물의 형상과 각 방향으로의 절삭력 등을 예측하는 것을 의미한다. 이러한 시뮬레이션은 다양하게 변하는 공구와 가공물 사이의 맞물림(Cutter-Workpiece Engagement, CWE)을 효율적으로 연산할 수 있어야 하며, 전체 공구경로에 대한 시뮬레이션이 가능해야 하므로 연산 속도가 빨라야 한다. Devore는 NC 밀링가공 시뮬레이션을 위해 공구 날의 작은 미세 요소에 의해 발생하는 힘들을 합하여 밀링가공의 절삭력을 계산하는 방법을 제시하였다[4,5]. 이후 이러한 방법을 기반으로 보다 정확하고 효율적으로 절삭력을 예측하는 수많은 연구가 진행되었다. 하지만, 현재까지 NC 가공에서의 절삭력 시뮬레이션에 있어 빠른 연산속도를 기반으로 한 실시간 시뮬레이션을 강조하는 연구는 찾아보기 힘들다. 실가공 동기화된 실시간 디지털트윈에서의 절삭력 시뮬레이션은 CNC에서 짧은 시간 간격으로 이송축 위치 등의 정보를 입력 받아 다음 정보 입력 시까지 시뮬레이션을 수행하는 것을 의미하며, 이는 실제 공정보다 빠른 절삭 시뮬레이션이 필수적이므로, 향후 고속 절삭력 시뮬레이션에 대한 연구 및 기대가 증가할 것으로 생각된다.

이 연구에서는 NC 밀링가공에서의 공구경로에 따른 재료 제거 및 절삭력 고속 시뮬레이션 방법을 제시한다. 3차원 기하 형상을 표현하는 가장 간단한 방법 중 하나이며 전산화에 적합한 복셀(Voxel)을 기반으로 가공물의 재료 제거를 시뮬레이션하고, 공구의 미소 요소에 의해 제거되는 복셀 요소 정보를 이용하여 기존의 복잡한 CWE 및 칩 두께 계산 없이 절삭력을 빠르게 연산한다. 제안한 시뮬레이션을 검증하기 위하여 가공 정밀도 평가용 시편(KS B ISO 10791-7 A160)의 가공을 수행하여 측정된 절삭력 및 실제 가공 시간과 시뮬레이션된 절삭력 및 연산 시간을 비교하였다.

개발된 시뮬레이션 모델의 예측정확도 및 연산성능을 검증하기 위하여, 밀링가공 정밀도 평가용 시편(KS B ISO 10791-7 A160)의 가공을 수행하여 측정 절삭력과 시뮬레이션 결과를 비교하였다. Fig. 4는 가공 실험 모습을 보여준다. 가공 실험은 수직형 3축 장비(SM400DH, SMEC)를 사용하여 수행하였으며, 알루미늄 6061 소재의 시편과 직경 및 날 수가 각각 10 mm 3날, 12 mm 3날, 14 mm 2날인 세 개의 공구를 사용하였다. 다축 공구동력계(9257B, Kistler)를 사용하여 가공 중 X, Y, Z 방향의 절삭력을 측정하였다. 약 25 Hz 빈도로 CNC 제어기(M800, Mitsubishi)로부터 시간, 이송축 위치, 이송속도, 주축 회전속도, 공구번호를 입력받아 그 시간 간격 동안 측정된 평균 절삭력을 동기화하여 함께 기록하였다.

Fig. 4 
Experimental setup

위와 같이 저장된 시간, 이송축 위치, 주축 회전속도, 공구 정보를 입력으로 시뮬레이션을 수행하였다. 시뮬레이션 프로그램은 C언어로 작성되었으며, 병렬 연산은 수행하지 않았다. 시뮬레이션에 사용된 연산 프로세서는 i7-6700K (Intel)이다. Table 1은 시뮬레이션 조건을 나타낸다.

Fig. 5는 실제 가공된 가공물과 복셀 정보를 이용하여 가시화한 시뮬레이션된 가공물 형상을 나타낸다. Fig. 6은 CNC 정보 입력마다 계산된 평균 절삭력과 측정된 평균 절삭력을 비교한 결과를 보여준다. Figs. 6(a)와 6(b)는 각각 X방향과 Y방향의 측정 및 시뮬레이션 절삭력을 보여주며, Figs. 6(c)와 6(d)는 가공시간 1,701초부터 1,714초 사이의 절삭력을 별도로 나타낸다. 평엔드밀 절삭 모델이 정확하게 적용되지 않는 드릴링과 유사한 구멍 가공 구간과 이송계 급속이송으로 인해 측정된 힘을 제외하면 약 25 N 이내의 오차 범위에서 실제 절삭력과 유사하게 예측이 가능함을 확인하였다. 또한, 0.09 mm의 비교적 작은 메시 크기를 적용하고 최신 연산 프로세서를 사용하지 않았음에도 불구하고, 시뮬레이션 시간은 638초로 실제 가공 시간 2,273초보다 약 3.5배 빠른 속도로 시뮬레이션이 가능함을 확인하였다.

Fig. 5 
Comparison between machined and simulated workpiece

Fig. 6 
Comparison between measured and simulated cutting force

이 연구에서는 NC 밀링가공에서의 공구경로에 따른 재료 제거 및 절삭력 고속 시뮬레이션 모델을 제시하였다. NC 밀링가공에서는 공구경로에 의해 공구와 가공물과의 기하 관계가 복잡하고 다양하게 변하게 되는데, 고속 시뮬레이션을 위해서는 이에 연산 복잡도 및 속도가 영향을 받지 않는 단순하고 효과적인 시뮬레이션 모델이 필요하다. 이를 위해 전산화가 용이한 복셀 기반의 가공물 모델과 공구중심점과 날의 미소 요소 위치로 정의된 공구 모델을 개발하였다. 제시된 모델은 공구 각 미소 요소가 지나간 영역 내의 가공물 복셀 요소들을 제거하여 재료 제거를 시뮬레이션하고, 제거된 복셀 개수만을 이용하여 간단하게 미변형 칩 두께 및 절삭력을 계산한다. 기존의 복잡한 CWE 및 칩 두께 계산 없이 절삭력을 빠르게 연산할 수 있으며, 이산화된 순간과 공간에 대한 절삭력 계산이 아닌 아주 복잡한 형상 가공에 대응할 수 있는 세밀하고 연속적인 시뮬레이션이 가능하다.

개발된 시뮬레이션 모델을 검증하기 위하여 가공 정밀도 평가용 시편(KS B ISO 10791-7 A160)의 가공을 수행하여 측정된 절삭력 및 실제 가공 시간과 시뮬레이션된 절삭력 및 연산 시간을 비교하였다. 전체 가공을 대상으로 25 N 이내의 오차로 실제 절삭력을 유사하게 예측하였고, 시뮬레이션 시간은 실제 가공시간의 28%에 불과해 실시간 시뮬레이션이 가능함을 검증하였다. 제안된 시뮬레이션은 가공품질 예측, 공정·장비 진단 및 최적화를 위한 공작기계 실시간 디지털트윈의 핵심기술로 활용될 것으로 기대된다.