토크 입력 완료 0 / 90 케이스
형식: No, M_inner, M_cage, M_element, M_drag — 탭/쉼표 구분, 한 줄에 한 케이스. 엑셀에서 5개 열을 복사해 붙여넣으면 됩니다.
| No | Case | Type | Z |
ωi [rpm] | ωc [rpm] | ωR [rpm] |
Minner [N·m] | Mcage [N·m] | Melement [N·m] | Mdrag [N·m] |
PCFD [W] | PISO [W] | PSKF [W] | 입력자 |
PCFD = Σ|Pi| (수식·부호 규약은 Methods §4).
Melement는 전동체 1개의 자전축 기준 토크.
입력은 팀 공유 DB에 저장되어 실시간 동기화.
빨간색 = PCFD가 [0.1×min(PISO,PSKF), 10×max] 범위 밖 → 해석 확인 필요.
1 · Kinematics — 속도 경계조건
Harris & Kotzalas, Rolling Bearing Analysis (5th ed.), Eq. 10.14–10.15. 내륜 회전, 외륜 고정, 순수 구름(pure rolling) 가정.
ωc = ωi2(1 − γ)
ωR = dm ωi2Dw(1 − γ2)
γ = Dw cos αeffdm
ωR: 케이지 동행좌표계(MRF) 기준 전동체 자전 속도, 방향은 내륜과 반대. αeff — DGBB: 0°, 4PCBB: 30° (nominal contact angle), TRB: 롤러 축각 β = αcup − ε/2. TRB의 경우 공통 꼭지점 순수구름 해와 동치:
ωc = ωi sin αisin αi + sin αo
αi = αcup − ε, αo = αcup
전동체 회전축 단위벡터 = cos β · x̂ + sin β · r̂ (BRG 경계조건 파일의 축벡터와 전수 일치 확인).
2 · Oil volume fraction XSNU
단상 CFD의 윤활유–공기 등가물성에 사용하는 oil volume fraction. Parker(1984) oil-injection 상관식의 oil bath 재피팅.
XSNU = 64.7 × 107ω0.26 dm4.80
ω = ωi [rpm], dm [mm]
3 · Lubricant & air properties — 온도 의존 물성
윤활유(ATF, OSU 측정 기반 상관식)와 공기의 온도 의존 물성. TK = T + 273.15.
ρoil = 1024.4 − 0.6388 TK [kg/m³]
μoil = 1012 exp(−2.303 · 15.91 T + 96.3621.3 + T + 96.36) [Pa·s]
α1 = 3.5952 × 10−8 · T−0.22133 [Pa⁻¹] (압력–점도 계수, EHL용)
ρair = pR TK, p = 101325 Pa, R = 287.05
μair = 1.716 × 10−5 (TK273.15)1.5 383.55TK + 110.4 [Pa·s] (Sutherland)
윤활유 (OSU)
| T [°C] | ρ [kg/m³] | μ [Pa·s] |
공기
| T [°C] | ρ [kg/m³] | μ [Pa·s] |
케이스별 표기 온도 T는 점도식 역산값 (샘플 μ가 OSU 점도식과 일치, T = 20–96.6 °C). 샘플 ρ는 μ와 선형 상관으로 생성됨.
4 · Power loss — CFD 토크 → 손실, 참조식 비교
P = M · ω9.5493 [W]
Pchurning = Z · Melement · ωR9.5493
PCFD = Σ |Pi|
(M [N·m], ω [rpm])
- Inner ring: ωi · Cage: ωc · Element churning: Melement(전동체 1개) × ωR(<0) × Z · Drag: ωc
참조식 — oil bath, 유면 = 최하단 전동체 중심, 오일 실물성 ν = μ/ρ. 무부하(회전) 손실 항만 비교. 원문 대조 완료 (ISO/TR 14179-1:2001 Eq. 27–29, Table 5, Fig. 3 · SKF General Catalogue drag losses, Table 4, Diagram 4).
ISO 14179-1: M0 = 10−10 f0 (ν ωi)2/3 dm3 [N·m], ν ωi ≥ 2000
f0 = f0,min + (f0,max − f0,min) HDOR, H = D − dm2, DOR = D
- Table 5 (단열, oil bath): (f0,min, f0,max) = DGBB (2, 4) · 4PCBB (6, 12) · TRB (4, 8). H: 베어링 최하단(외륜 외경)→유면 (Fig. 3)
- 참조식의 dm = 0.5(d+D) (ISO Eq. 10 · SKF 정의) — 운동학·XSNU는 실제 피치원 Dpw 사용
SKF: Mdrag = 0.4 VM Kball dm5 ωi2 + M2 (ball)
Mdrag = 4 VM Kroll Cw B dm4 ωi2 + M2 (roller)
M2 = 1.093×10−7 ωi2 dm3 (ωi dm2 ftν)−1.379 Rs
- H(SKF): 외륜 궤도–전동체 최저 접촉점→유면. 접촉점 직경 추정 = ball류 0.5(D+D1), TRB는 D → Hball = (D+D1)/4 − dm/2, HTRB = (D−dm)/2
- VM(x = H/dm): Diagram 4 근사 — ball 0.0006(1−e−5.04x1.5), roller 0.0014(1−e−2.84x1.5). SKF Bearing Select 온라인 계산기 실행값(32011 X·6011, ν=9.83 cSt, H=8.75 mm, 4000 rpm → Mdrag 76 / 19.8 N·mm)으로 보정·검증
- Kball = irwKZ(d+D)/(D−d)×10−12, Kroll = KLKZ(d+D)/(D−d)×10−12 · Table 4: KZ = 3.1 (DGBB, four-point) / 6.0 (TRB), KL = 0.7 (TRB) · Cw = f(lD), lD = 5KLB/dm · Rs = 0.36dm2(t−sint)fA, t = 2cos−1((0.6dm−H)/0.6dm)
- 속도 의존성: PISO ∝ ω5/3, PSKF ∝ ω3 (drag 1항) — 고속·고유면에서 SKF가 역전. 본 유면(최하단 전동체 중심)에서는 VM이 작아 대부분 케이스에서 PSKF < PISO
판정 기준: PCFD ∉ [0.1×min(PISO, PSKF), 10×max(PISO, PSKF)] → 빨간색 경고. 표준식의 실측 대비 산포(±2–3×)를 감안한 범위.
5 · Data QA — 원본 데이터 검증에서 발견·처리한 사항