유한요소 전자기

편집 역사 토론
전자기해석 전자기학 마지막 수정: 2026-07-07 09:14:32
유한요소 전자기
Finite Element Electromagnetics
분야전자기학 × 수치해석 × 전기기기 설계
지배 방정식맥스웰 방정식 (저주파 근사)
핵심 기법에지 요소, 자기 벡터 퍼텐셜 $\mathbf{A}$
주 응용모터·변압기·와전류·자기장 해석
대표 소프트웨어Ansys Maxwell, JMAG, Altair Flux, COMSOL

1. 개요[편집]

고주파는 FDTD가 잡고, 저주파는 FEM이 잡는다. 그리고 모터 설계자의 밥줄은 저주파에 있다.

유한요소 전자기(Finite Element Electromagnetics)는 맥스웰 방정식에서 변위전류를 무시하거나 근사할 수 있는 저주파·준정적(magneto-quasistatic) 영역의 자기장·와전류·힘 문제를 유한요소법으로 이산화해 푸는 해석 분야다. 파장이 소자 크기보다 훨씬 큰 60 Hz 전력기기부터 수 kHz 인버터 구동 모터까지, 전파(radiation)가 아니라 자속(flux)이 주인공인 세계를 다룬다.

안테나 해석처럼 전파가 날아다니는 고주파 문제에서는 시간영역 FDTD나 모멘트법이 강세지만, 모터·변압기·솔레노이드·자기부상처럼 “장이 갇혀 있는” 문제에서는 복잡한 형상과 비선형 자성체를 자연스럽게 담아내는 유한요소법이 압도적으로 유리하다.1 회전 기계의 토크 리플 한 방울까지 잡아내야 하는 전기차 모터 시대에 이 분야의 몸값이 유난히 높아진 이유다.

2. 지배 방정식과 저주파 근사[편집]

완전한 맥스웰 방정식에서 앙페르 법칙은 다음과 같다.

∇×H=J+∂D∂t\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t}

저주파에서는 변위전류 ∂D/∂t\partial \mathbf{D}/\partial t가 전도전류 J\mathbf{J}에 비해 무시할 만큼 작다. 이 근사가 “준정적(quasi-static)“의 핵심이다. 자속밀도를 자기 벡터 퍼텐셜 A\mathbf{A}로 표현(B=∇×A\mathbf{B} = \nabla \times \mathbf{A})하면 와전류를 포함한 지배 방정식은 다음 형태로 정리된다.

∇×(1μ∇×A)+σ∂A∂t=Js\nabla \times \left( \frac{1}{\mu} \nabla \times \mathbf{A} \right) + \sigma \frac{\partial \mathbf{A}}{\partial t} = \mathbf{J}_s

여기서 μ\mu는 투자율, σ\sigma는 도전율, Js\mathbf{J}_s는 외부 전원 전류밀도다. 강자성체에서는 μ\mu가 B\mathbf{B}에 의존하는 비선형 함수이므로, 이 방정식은 매 시간 스텝마다 뉴턴-랩슨법으로 반복해 풀어야 한다. 자기 포화(saturation)라는 성가신 비선형성이 바로 여기서 튀어나온다.

3. 에지 요소: 절점 요소가 실패하는 이유[편집]

전자기 FEM이 다른 물리 분야와 결정적으로 다른 점이 하나 있다. 절점(node)에 스칼라 자유도를 놓는 전통적 갤러킨 방법이 벡터 장 A\mathbf{A}에서는 처참하게 실패한다는 것이다.

절점 기반 벡터 요소는 서로 다른 투자율의 경계에서 접선 성분 연속 조건을 강제로 만족시키려다 물리적으로 존재하지 않는 가짜 해(spurious solution)를 만들어 낸다. 이를 해결한 것이 에지 요소(edge element), 즉 자유도를 절점이 아니라 요소의 변(edge)에 놓는 벡터 요소다. 이들은 수학적으로 H(curl)H(\text{curl}) 공간에 속하며, 다음 성질을 자연스럽게 보장한다.

  • 접선 성분은 요소 경계를 가로질러 연속 (물리적으로 옳음)
  • 법선 성분은 불연속 허용 (매질 경계에서 실제로 그렇다)
  • ∇×\nabla \times 연산이 이산 공간에서 정확히 표현됨

에지 요소의 등장은 전자기 FEM을 실용 궤도에 올린 사건이었다. Whitney 요소라고도 불리는 이 저차 벡터 요소는 오늘날 거의 모든 상용 코드의 기본값이다.2

4. 와전류와 스킨 효과[편집]

교류 자기장이 도체에 침투하면 패러데이 법칙에 의해 유도 전류, 즉 와전류(eddy current)가 흐른다. 와전류는 도체 표면으로 밀려나(스킨 효과) 손실을 만들고 자속을 밀어낸다. 침투 깊이는 다음 표피 깊이(skin depth)로 특징지어진다.

δ=2ωμσ\delta = \sqrt{\frac{2}{\omega \mu \sigma}}

문제는 이 δ\delta가 고주파에서 매우 얇아진다는 점이다. 표피 깊이를 제대로 해상하려면 도체 표면에 요소를 극도로 촘촘히 깔아야 하고, 이는 적응 격자 세분화나 경계층 요소가 필수가 되는 이유다. 변압기 철심의 규소강판을 얇게 적층하는 것도, 모터 설계에서 이 와전류 손실을 억제하기 위한 물리적 대응이다.

응용관심 주파수핵심 물리
전력 변압기50/60 Hz자속 분포, 철손·동손
유도 전동기수십 Hz~kHz회전자 와전류, 토크
인버터 구동 모터~수 kHz고조파 손실, 온도
무선 충전 코일수십~수백 kHz스킨·근접 효과

5. 회전 기계와 토크 계산[편집]

모터 해석의 백미는 회전자를 실제로 돌리는 것이다. 회전자와 고정자 사이 공극(air gap)에 미끄럼 경계(sliding surface)나 이동 대역(moving band) 기법을 두어, 매 스텝마다 격자를 재접합하며 상대 회전을 구현한다. 토크는 공극을 감싸는 폐곡면에서 맥스웰 응력 텐서(Maxwell stress tensor)를 적분해 얻는다.

T=∮Sr×(1μ0(B⋅n^)B−B22μ0n^)dST = \oint_{S} \mathbf{r} \times \left( \frac{1}{\mu_0} (\mathbf{B} \cdot \hat{\mathbf{n}})\mathbf{B} - \frac{B^2}{2\mu_0}\hat{\mathbf{n}} \right) dS

토크 리플, 코깅 토크, 역기전력 파형 같은 지표가 여기서 나오며, 이들은 전기차의 정숙성과 효율을 좌우한다. 실무에서는 이 전자기 해석 결과를 열전달 해석이나 구조해석과 엮어 다물리 연성해석으로 확장하는 경우가 많다. 자기력이 계산되면 그것이 곧 열원이자 진동원이기 때문이다.

6. 현업에서의 현실[편집]

  • 모터 설계자의 90%는 3D 대신 2D 단면 해석으로 산다. 축방향 대칭 덕분에 계산량이 극적으로 줄기 때문. 3D는 엔드 와인딩 누설이 궁금할 때만 꺼낸다.
  • 비선형 BB-HH 곡선 데이터를 강판 제조사에서 못 받으면 해석 정확도의 절반이 날아간다. 재료 데이터가 곧 실력이다.
  • 영구자석 감자(demagnetization)를 무시하고 설계했다가 고온에서 자석이 죽어버리는 사고가 은근히 흔하다.
  • “시뮬 토크는 나오는데 실측이 안 나와요”의 상당수는 적층 철심의 등가 투자율을 잘못 넣은 탓이다.

7. 관련 문서[편집]

8. Footnotes[편집]

  1. FDTD로 60 Hz 모터를 풀겠다고 나서면 파장 5000 km를 해상하려고 우주적 크기의 격자를 짜야 한다. 도구는 주파수에 맞춰 골라야 하는 법이다. ↩

  2. 에지 요소가 없던 시절엔 벡터 FEM 결과에서 “고윳값이 0 근처에 잔뜩 몰린 유령 모드”가 나와 연구자들을 괴롭혔다. Whitney(1957)의 미분형식 이론이 뒤늦게 구원투수로 등판했다. ↩