영구자석 동기 직선 전동기의 시스템 수준 모델링, 비선형 시뮬레이션 및 다목적 최적화

1. 시스템 수준 기초 모델링: 회전형 전동기 기반의 등가 직선 모델

반도체 노광 장비의 웨이퍼 스테이지나 고속 정밀 플로터와 같은 초정밀 위치 결정 시스템에서는 영구자석 동기 직선 전동기(PMLSM)가 핵심 구동기로 사용됩니다. 이러한 시스템의 고정밀 제어를 위해서는 전자기적 특성과 기계적 동역학을 정확히 반영하는 수학적 모델이 필수적입니다. 초기 설계 단계에서는 3차원 유한요소해석(FEA) 대신, 회전형 영구자석 동기 전동기(PMSM)의 수학적 모델을 직선 운동으로 매핑한 집중 정수 모델을 활용하여 제어 루프의 타당성을 빠르게 검증하는 것이 효율적입니다.

1.1 수학적 모델 매핑 및 파라미터 변환

PMLSM의 수학적 모델은 PMSM의 모델을 직선 좌표계로 전개하여 유도할 수 있습니다. 주요 파라미터의 매핑 관계는 다음과 같습니다.

속도 매핑: 회전형의 각속도 $\omega$는 직선형의 선속도 $v$로 변환됩니다. $v = \frac{2\pi n \tau}{60}$ (여기서 $\tau$는 극피치, $n$은 rpm).
추력 방정식: 전자기 추력 $F_{em}$은 회전형의 토크 방정식을 반경 대신 극피치 기반으로 변환한 것입니다. $F_{em} = \frac{3\pi}{p\tau} [\psi_f i_q + (L_d - L_q)i_d i_q]$ ($p$는 극쌍수, $\psi_f$는 영구자석 쇄교 자속, $L_d, L_q$는 d-q축 인덕턴스).
기계적 파라미터: 회전 관성 모멘트($J$)는 가동자의 질량($M$)으로, 회전 점성 마찰 계수는 직선 점성 마찰 계수($B_l$)로 대체됩니다.

1.2 Simulink/Simscape 기반 등가 모델 구현

MATLAB/Simulink 환경에서 Simscape Electrical 라이브러리를 활용하면 물리 기반의 등가 모델을 신속하게 구축할 수 있습니다. 아래 스크립트는 PMSM 블록을 직선 운동 인터페이스와 연결하기 위한 핵심 파라미터 초기화 과정을 보여줍니다.

% PMLSM 기초 등가 모델 초기화 및 파라미터 매핑 스크립트
modelName = 'PMLSM_Equiv_System';
new_system(modelName);
open_system(modelName);

% 전자기 및 기계적 파라미터 정의
motorParams.poles = 4;             % 극 수 (Pole pairs = 2)
motorParams.tau = 0.025;           % 극피치 (m)
motorParams.flux_pm = 0.085;       % 영구자석 쇄교 자속 (Wb)
motorParams.Ld = 1.2e-3;           % d축 인덕턴스 (H)
motorParams.Lq = 1.2e-3;           % q축 인덕턴스 (H)
motorParams.Rs = 0.65;             % 고정자 저항 (Ohm)
motorParams.mass = 0.15;           % 가동자 질량 (kg)

% 회전-직선 운동 변환 비율 계산
% 직선 속도 v = 각속도 omega * (극피치 * 극쌍수)
motorParams.lead_screw_pitch = motorParams.tau * (motorParams.poles / 2); 

% Simscape 블록 추가 및 파라미터 할당
add_block('simscape/Electrical/Electromechanical/Electromechanical/Permanent Magnet Synchronous Machine', [modelName '/PMSM_Core']);
set_param([modelName '/PMSM_Core'], 'p', num2str(motorParams.poles/2));
set_param([modelName '/PMSM_Core'], 'FluxPM', num2str(motorParams.flux_pm));
set_param([modelName '/PMSM_Core'], 'Ld', num2str(motorParams.Ld));
set_param([modelName '/PMSM_Core'], 'Lq', num2str(motorParams.Lq));
set_param([modelName '/PMSM_Core'], 'Rs', num2str(motorParams.Rs));

% 기계적 인터페이스 및 부하 블록 구성
add_block('simscape/Mechanical/Translational/Translational Mass', [modelName '/Moving_Mass']);
set_param([modelName '/Moving_Mass'], 'mass', num2str(motorParams.mass));

2. 고급 비선형 모델링: 외란 및 동역학 요인 통합

기초 선형 모델은 제어기 설계의 초기 검증에는 유용하지만, 실제 PMLSM에서 발생하는 비선형 외란을 반영하지 못합니다. 고정밀 서보 제어를 위해서는 치슬력(Cogging force)과 비선형 마찰력을 모델에 통합하여 시뮬레이션의 충실도를 높여야 합니다.

2.1 공간 고조파를 고려한 치슬력 모델링

치슬력은 고정자 철심의 슬롯과 영구자석 간의 자속 변화로 인해 발생하는 위치 의존적 추력 맥동입니다. 단일 정현파 모델보다는 공간 고조파를 포함한 모델이 실제 현상을 더 정확히 모사합니다.

function F_cogging = calculate_cogging_force(position)
    % 다중 공간 고조파를 고려한 치슬력(Cogging Force) 계산
    pole_pairs = 2;
    pole_pitch = 0.025;
    
    % 고조파 진폭 (N) - FEA 결과 기반 피팅 값
    A1 = 0.45; % 기본파 진폭
    A3 = 0.12; % 3차 고조파 진폭
    A6 = 0.05; % 6차 고조파 진폭
    
    % 공간 전기각 계산
    theta_e = (2 * pi * pole_pairs * position) / pole_pitch;
    
    % 치슬력 합성 (기본파 및 고조파 중첩)
    F_cogging = A1 * sin(theta_e) + A3 * sin(3 * theta_e) + A6 * sin(6 * theta_e);
end

2.2 비선형 마찰 동역학 통합

직선 안내면에서 발생하는 마찰력은 저속 구간에서 스트리벡(Stribeck) 효과와 쿨롱 마찰을 유발하여 위치 결정 오차의 원인이 됩니다. Simscape의 Translational Friction 블록을 활용하여 정지 마찰, 쿨롱 마찰, 그리고 점성 마찰을 포함하는 복합 마찰 모델을 기계적 부하단에 직렬로 연결합니다. 이를 통해 가동자의 시작 및 정지 구간에서 발생하는 추력 피크와 위치 지체 현상을 시뮬레이션에서 관찰할 수 있습니다.

3. 다목적 최적화: 제어 알고리즘 및 전자기 구조 설계

비선형 모델이 구축되면, 시스템의 성능을 극대화하기 위해 제어 파라미터와 전자기 구조 파라미터에 대한 최적화를 수행합니다.

3.1 유전 알고리즘을 활용한 슬라이딩 모드 제어기 파라미터 최적화

치슬력과 마찰 외란에 강인한 슬라이딩 모드 제어(SMC)를 적용할 때, 슬라이딩 평면 계수와 지수 도달칙 파라미터의 수동 튜닝은 한계가 있습니다. 유전 알고리즘(GA)을 사용하여 위치 추적 오차의 적분 절댓값(IAE)을 최소화하는 방향으로 파라미터를 자동 최적화할 수 있습니다.

% 유전 알고리즘(GA)를 활용한 슬라이딩 모드 제어기(SMC) 파라미터 최적화
objective_func = @(vars) evaluate_smc_performance(vars);
num_vars = 4; % 최적화 변수: [c1, c2, epsilon, k]

% 탐색 공간 제약 조건 설정
lower_bounds = [5, 5, 0.05, 0.5];
upper_bounds = [150, 150, 5.0, 8.0];

% GA 옵션 구성
ga_opts = optimoptions('ga', ...
    'PopulationSize', 40, ...
    'MaxGenerations', 100, ...
    'CrossoverFraction', 0.8, ...
    'Display', 'final', ...
    'PlotFcn', @gaplotbestf);

% 최적화 실행
[opt_vars, min_iae] = ga(objective_func, num_vars, [], [], [], [], lower_bounds, upper_bounds, [], ga_opts);

% 최적화된 파라미터를 Simulink 워크스페이스로 전달
c1_opt = opt_vars(1);
c2_opt = opt_vars(2);
epsilon_opt = opt_vars(3);
k_opt = opt_vars(4);

3.2 전자기 구조 파라미터 최적화 (Maxwell 및 Isight 연동)

제어 알고리즘의 한계를 넘어 물리적인 추력 맥동 자체를 줄이기 위해서는 전자기 구조 최적화가 필요합니다. Ansys Maxwell 2D 과도 해석 모델과 Isight 최적화 프레임워크를 연동하여 극피치($\tau$)와 영구자석 두께($h$)를 설계 변수로 설정합니다. 목적 함수는 치슬력 피크($F_{cog0}$) 최소화 및 rated 추력 유지로 설정하며, 다목적 최적화 알고리즘(예: NSGA-II)을 적용합니다. 이를 통해 고정자 슬롯 폭과의 반정수 슬롯 조합을 유도하여 치슬력을 획기적으로 저감하면서도 필요한 정격 추력을 만족하는 최적의 전자기 설계를 도출할 수 있습니다.

태그: PMLSM Matlab Simulink Simscape Ansys Maxwell

6월 1일 02:01에 게시됨

괴물 클럽