두 명의 엔지니어가 동일한 기어 쌍을 모델링하더라도 서로 다른 강성 계산 방법을 사용하면 두 자릿수 차이가 나는 예측 결과를 얻을 수 있습니다. 이러한 차이는 반올림 오차가 아니라, 모든 후속 동적 결과에 조용히 영향을 미치는 방법론적 차이에서 비롯됩니다. 대부분의 기어 엔지니어는 분석 맥락에 적합한지 여부를 고려하지 않고 소프트웨어에서 제공하는 방법이나 교수가 가르친 방법을 기본적으로 사용합니다.
일반적인 휠 바디를 사용하는 표준 산업용 스퍼 기어 및 헬리컬 기어의 경우, ISO 6336 단일 치 강성 공식으로 대부분의 정격 및 예비 동적 계산을 처리할 수 있습니다. 엔지니어들은 자체적인 오류를 유발하는 맞춤형 포텐셜 에너지 구현 방식이나, 표준 공식으로 몇 초 만에 얻을 수 있는 결과의 5% 이내의 오차를 얻기 위해 몇 주씩 설정 시간을 소모하는 유한 요소 해석(FEA) 모델에 지나치게 의존하는 경향이 있습니다. FEA는 결합된 메시, 경량 물체 또는 비표준 접촉 형상과 같은 경우에만 진정으로 필요합니다.
기어 맞물림 강성을 결정하는 요소는 무엇인가?
기어 톱니는 맞물렸다가 떼는 동작을 반복하는 고정된 블록이 아닙니다. 기어 톱니는 강성이 가변적인 스프링과 같아서, 그 유연성이 시간에 따라 변합니다. 메싱 사이클그리고 강성 변화가 기어 동역학에서 주요 여기원이라는 것입니다.
5가지 규정 준수 구성 요소
단일 치형의 강성은 직렬로 연결된 다섯 가지 구성 요소로 분해됩니다. 각 구성 요소는 서로 다른 변형 메커니즘을 나타냅니다.
- 굽힘 순응도(k_b) 치아는 횡방향 하중을 받을 때 캔틸레버 보처럼 휘어집니다. 이러한 현상은 인벌류트 프로파일의 접촉점이 많을수록 두드러지게 나타납니다.
- 전단 컴플라이언스(k_s) — 횡방향 전단 변형은 단면적의 약 1.2배에 비례합니다. 굽힘에 비해 상대적으로 작지만, 짧은 치아의 경우 무시할 수 없습니다.
- 축방향 압축 컴플라이언스(k_a) — 치축에 수직인 하중 성분은 치축을 압축합니다. 이는 높은 압력각이나 큰 나선각을 가진 헬리컬 기어에서만 중요합니다.
- 헤르츠 접촉 준수(k_h) — 맞닿는 면 사이의 접촉 영역에서의 표면 변형. 이 요소는 접촉 형상이 작용선을 따라 천천히 변화하기 때문에 메시 생성 주기 동안 비교적 일정합니다.
- 필렛 기초 준수(k_f) — 기어 본체 자체는 이뿌리가 림과 만나는 부분에서 변형됩니다. 이 변형은 맞물림 주기 초반, 이뿌리 근처에서 주로 발생하며, 접촉점이 높아질수록 인벌류트 굽힘으로 인해 변형량이 줄어듭니다.
결합된 단일 치아 강성은 직렬 컴플라이언스 관계인 1/k_t = 1/k_b + 1/k_s + 1/k_a + 1/k_h + 1/k_f를 따릅니다. 여러 쌍의 치아가 동시에 접촉하는 경우, 각 치아 쌍의 강성은 병렬로 합산됩니다.
베버와 바나셰크가 1953년에 제시한 분석 프레임워크에서 유래한 이러한 분해 방식은 단순한 역사적 유물이 아닙니다. 이는 현대 기어 분석 소프트웨어의 근간을 이루며, 전 세계 산업용 기어 등급 평가에 사용되는 ISO 6336 강성 계수의 토대가 되고 있습니다. 유한 요소 해석(FEA)을 바로 사용하는 엔지니어들은 70년 동안 검증된 분석 작업을 건너뛰는 셈입니다.
싱글 페어 vs 더블 페어 접촉 영역
시간에 따라 변하는 맞물림 강성(TVMS)은 톱니 쌍이 접촉에 들어가고 나올 때 진동합니다. 접촉비가 1에서 2 사이인 스퍼 기어의 경우, 이 주기는 단일 톱니 쌍 영역(한 톱니 쌍이 모든 하중을 지탱하고 강성이 낮음)과 이중 톱니 쌍 영역(두 톱니 쌍이 하중을 분담하고 강성이 높음) 사이를 번갈아 가며 나타납니다.
절대적인 강성 값이 아니라 강성 변화 패턴이 동적 가진을 유발합니다. 뤄양 과학기술대학의 연구진은 톱니 위상차를 조절하여 이러한 변화 패턴을 조작함으로써 기어 진동의 가속도 진폭을 35~39% 감소시킬 수 있음을 입증했습니다. 해석적 방법만으로 계산 가능한 이 강성 변화 패턴만으로도 대책을 설계하기에 충분했습니다.
메쉬 강성에 대한 폐쇄형 공식
ISO 6336 단일 치 강성
ISO 6336은 단일 치형 강성 c'를 단위 하중으로 정의합니다. 얼굴 너비 작용선을 따라 한 쌍의 치아를 1마이크로미터 변형시키는 데 필요한 조건입니다. 이 표준은 엔지니어가 각 순응 구성 요소를 개별적으로 모델링하도록 요구하는 대신 경험적으로 검증된 계수를 제공합니다.
일반적인 스퍼 기어의 경우, c'는 모듈, 톱니 수 및 프로파일 변위에 따라 일반적으로 14~20 N/(mm·μm) 범위에 속합니다. 전체 맞물림 강성 c_gamma는 접촉비를 고려합니다. 이중 쌍 접촉 시 강성은 대략 두 배가 됩니다.
ISO 방식은 평균 경사 강성(총 맞물림력을 총 맞물림 변위로 나눈 값)을 사용합니다. 이는 산업용 기어 계산 작업의 대부분을 차지하는 정적 기어 등급에 대한 정확한 공식입니다.
위치 에너지 방법
위치 에너지 방법은 각 치아를 가변 단면 캔틸레버 보로 모델링하고 하중 하에서 저장된 탄성 변형 에너지로부터 강성을 계산합니다. 총 변형 에너지는 굽힘, 전단, 축압축, 헤르츠 접촉 및 필렛-기초 항을 합산합니다.
위치 에너지로부터 얻은 단일 치아 강성:
k_i = F^2 / (2 · U_total)
여기서 F는 적용된 메쉬 힘이고 U_total은 모든 변형 에너지 성분의 합입니다.
이 방법은 ISO 6336의 표로 정리된 계수에서 직접 제공하지 않는 메싱 주기 전체에 걸친 완전한 강성 곡선을 제공합니다. 그러나 이 방법은 정밀한 기하학적 모델링을 필요로 합니다. 기준 원에서 치아를 모델링하는 기존 방식은 검증된 표준에 비해 최대 8.6%의 오차를 발생시킵니다. 뿌리 원과 기준 원 사이의 치아 프로파일을 정확하게 처리하는 개선된 방식은 이 오차를 약 2.9%로 줄입니다.
해석적 정확도와 유한요소법 정확도 비교
해석적 방법(ISO 6336 및 적절하게 구현된 위치 에너지 방법 모두)은 루트 필렛 형상을 정확하게 모델링했을 때 표준 스퍼 기어에 대한 유한 요소 결과와 최대 5.72%의 차이를 보입니다. 실제 루트 필렛을 무시한 단순화된 해석 모델은 강성을 약 14% 과대평가하는 반면, 기존의 경험적 근사치는 최대 39%까지 과소평가합니다.
5.72%의 정확도 범위는 산업용 기어 등급 평가 및 예비 동적 분석에 충분합니다. 저는 엔지니어들이 ISO 공식으로 30초 만에 얻을 수 있는 결과와 동일한 정확도를 얻기 위해 2주 동안 유한 요소 해석(FEA) 메쉬 강성 모델을 구축하는 것을 본 적이 있습니다. FEA의 계산 오버헤드는 분석 맥락에서 진정으로 요구될 때만 정당화되며, 표준 형상의 경우 그러한 경우는 드뭅니다.
기어 매개변수가 강성에 미치는 영향
모듈, 면폭 및 치아 개수
모듈과 면폭은 메쉬 강성에 거의 선형적으로 영향을 미칩니다. 즉, 면폭이 두 배가 되면 전체 메쉬 강성도 두 배가 됩니다. 톱니 개수는 톱니 형상을 통해 강성에 영향을 미칩니다. 즉, 주어진 피치 직경에서 톱니 수가 많을수록 톱니가 더 얇고 유연해집니다.
압력각은 치면을 따라 하중 방향을 변화시킵니다. 압력각을 20도에서 25도로 증가시키면 치면의 밑부분이 넓어져 강성이 증가하지만, 베어링 하중 또한 증가합니다. 이는 설계상의 절충점이며, 임의로 조절할 수 있는 변수가 아닙니다.
필렛 반경의 역설
공구 끝단 필렛 반경을 증가시키면 전체 맞물림 주기 동안 맞물림 강성이 5~6% 감소합니다. 이는 직관과 반대되는 결과입니다. 필렛이 클수록 뿌리 부분이 더 강해질 것으로 예상되기 때문입니다. 하지만 그 메커니즘은 간단합니다. 필렛이 커지면 기초 강성(k_f)이 증가하여 이뿌리가 더 유연해지기 때문입니다. 접촉비가 높은 기어는 표준 접촉비 기어보다 필렛 반경 변화에 1.2배 더 민감합니다. 치아 변형 필렛 영역을 변경하면 강성 여기 패턴이 직접적으로 변경됩니다.
제조 공차는 강성을 통해 누적됩니다.
치형 두께, 끝단 직경 및 중심 거리의 제조 오차는 맞물림 강성을 통해 안전 계수에 직접적으로 영향을 미칩니다. 공칭 성능이 유사한 두 기어 쌍은 극명하게 다른 내구성을 보였습니다. 첫 번째 기어 쌍의 최대 정적 전달 오차 표준 편차는 두 번째 기어 쌍보다 6배 더 높았습니다(0.249μm 대 0.039μm). 치아 뿌리 스트레스 제조 오류가 없을 경우 명목상 1.001이었던 안전 계수는 공차 누적을 고려했을 때 0.988로 떨어져 최소 임계값인 1.0 미만이 되었습니다.
기어 정밀도 등급은 강성 예측에 허용되는 제조 오차의 범위를 결정합니다. 공칭 형상이 동일한 6등급 기어와 10등급 기어는 실제 강성이 톱니마다 다르게 변하기 때문에 실제 동적 응답이 다릅니다.
어떤 방법을 사용해야 할까요?
ISO 6336, 위치 에너지 해석, 유한 요소 해석(FEA)은 동일한 기어 쌍에 대해 각각 다른 강성 값을 산출하며, 각 값은 서로 다른 상황에서 타당합니다. 어떤 해석 방법이 분석에 적합한지는 다음 세 가지 질문을 통해 판단할 수 있습니다. 이 강성 값은 어떤 결정을 뒷받침하는가? 그 결정에 필요한 정확도는 어느 정도인가? 모델링하는 형상은 무엇인가?
정적 등급 및 예비 설계
ISO 6336 또는 AGMA 2001에 따른 정적 기어 등급 계산 시에는 ISO 6336의 단일 치 강성 c'와 해당 경험적 계수를 사용하십시오. 이는 평균 기울기 접근법으로, 정적 하중 분포 및 강도 계산에 적합한 공식입니다. 이 공식은 즉시 계산되며 수십 년간 산업 현장에서 검증되었습니다.
대부분의 산업용 기어 엔지니어링 작업은 이 범주에 속합니다. 감속기 크기 결정, 치근 안전 계수 검사, 동적 하중 계수 K_v 추정 등 모든 작업에는 ISO 6336에서 제공하는 강성 값이 사용되며, 이 값은 해당 결정을 내리는 데 충분한 정확도를 제공합니다.
동적 분석 및 NVH
시간 영역 동적 시뮬레이션 또는 NVH 예측을 위해서는 메시 생성 주기 동안의 강성 변화, 즉 TVMS 곡선이 필요합니다. 포텐셜 에너지 방법은 이 곡선을 해석적으로 생성합니다. 동적 해석은 정상 상태 하중 주변의 작은 진동을 다루기 때문에 평균 기울기가 아닌 국부 기울기 공식(기준 작동점 주변의 증분 강성)을 사용해야 합니다.
치형 모델을 엄격하게 설계해야 합니다. 특히 톱니 수가 42개 미만인 기어의 경우, 실제 루트 필렛 프로파일과 베이스 서클과 루트 서클 사이의 간격이 크게 벌어지므로 이 점을 반드시 고려해야 합니다. 위치 에너지 구현에 소홀하면 ISO 6336을 준수하지 않음으로써 얻는 오차보다 더 큰 오차가 발생합니다.
유한요소해석이 진정으로 필요한 경우
세 가지 시나리오에 대한 유한 요소 메쉬 강성 분석을 예약합니다.
- 구멍이나 얇은 막이 있는 경량 휠 본체. 휠 본체에 경량화 구멍을 뚫으면 최대 전달 오차가 9배 증가하여 0.79μm에서 7.5μm로 커집니다. 본체의 탄성이 치아의 탄성보다 훨씬 크기 때문에 치아 수준의 분석 모델로는 이러한 변화를 포착할 수 없습니다.
- 결합된 이중 나선형 메쉬. 기어 반쪽 사이의 축 방향 힘 변동은 단일 메쉬 해석 모델로는 표현할 수 없는 결합 여기 메커니즘을 생성합니다. 이러한 형상에 대해서는 유한 요소 기반 결합 메쉬 해석만이 유일하게 실현 가능한 접근 방식입니다.
- 비표준 접촉 형상. ISO 6336 계수표가 적용되지 않는 비대칭 프로파일, 극심한 프로파일 변화 또는 비정형적인 치형.
그 외의 모든 경우, 즉 솔리드 또는 기존 림 바디를 사용하는 표준 스퍼 및 헬리컬 쌍의 경우, 해석적 방법은 FEA와 5~6% 이내의 오차 범위 내에서 훨씬 적은 비용과 시간으로 결과를 제공합니다.
메쉬 강성 계산 시 흔히 발생하는 오류
치아 뿌리 원 대신 밑면 원을 기준으로 치아를 모델링합니다. 일반적인 인벌류트 스퍼 기어의 경우, 톱니 수가 정확히 42개일 때만 베이스 서클과 루트 서클이 같아집니다. 톱니 수가 42개 미만인 경우(실제 기어 설계의 대부분을 차지함), 루트 서클과 베이스 서클 사이의 톱니 부분을 무시하면 체계적인 오차가 발생합니다. 이러한 단순화를 사용하는 기존의 위치 에너지 계산 방식은 검증된 표준과 비교했을 때 최대 8.6%의 오차를 보입니다.
루트 필렛 형상을 무시합니다. 루트 필렛을 급격한 전환으로 근사화하는 단순화된 모델은 메쉬 강성을 약 14% 과대평가합니다. 필렛-기반 컴플라이언스는 사소한 보정 항이 아니라, 메싱 주기 초기에 강성 응답을 지배하는 요소입니다.
동적 모델에서 평균 경사 강성을 사용합니다. 평균 경사 강성과 국부 경사 강성은 넓은 토크 범위와 표면 개조가 적용된 기어에서 의미 있는 차이를 보이는 예측값을 생성합니다. 동적 시뮬레이션 모델에서 ISO 6336의 평균 경사 값(정적 정격에 적합)을 사용하면 고유 진동수 및 강제 응답 예측값이 부정확해집니다. 강성 정의를 해석 유형에 맞게 조정하십시오.
유한요소해석(FEA)을 자동으로 더 정확하게 처리하는 것. 유한요소해석(FEA)의 메쉬 강성 결과는 요소 메쉬, 접촉 공식, 경계 조건의 정확도에 따라 달라집니다. 접촉 영역에서 요소 크기가 너무 큰 불량한 FEA 모델은 잘 구현된 해석 공식보다 정확도가 떨어질 수 있습니다. 즉, 방법 자체가 정확도를 보장하는 것이 아니라 구현 방식이 정확도를 좌우합니다.
메쉬 강성에 대한 체계적인 접근 방식
모든 강성 계산은 먼저 그 계산 결과가 어떤 결론을 뒷받침하는지 자문하는 것부터 시작해야 합니다. 정적 기어 등급의 경우 ISO 6336 c'를 사용하고 다음으로 넘어가십시오. 시간 영역 동역학의 경우 적절한 루트 형상을 적용하여 위치 에너지 방법을 구현하십시오. 경량 동체 또는 결합 메쉬의 경우 해석 모델이 물리적 현상을 완벽하게 포착할 수 없으므로 유한 요소 해석(FEA)에 투자해야 합니다.
치별 강성 계산에는 실제 루트 필렛 형상, 42개 미만 톱니 기어의 경우 베이스-루트 원 형상, 그리고 해석 맥락에 맞는 정확한 경사 공식이 고려되어야 합니다. 이 세 가지를 제대로 고려하면 계산 방법의 중요성은 대부분의 엔지니어가 생각하는 것보다 훨씬 낮습니다. 하지만 이 중 하나라도 잘못되면 아무리 정교한 계산 기법을 사용하더라도 결과는 변하지 않습니다.
