서론: 리니어 압축기 설계와 진동 해석의 중요성
본 포스팅에서는 파이썬을 이용한 수치 적분(Numerical Integration) 방식으로 시스템의 동적 거동을 예측하는 로직을 구현해 보겠습니다. 이는 초기 파라미터 튜닝 단계에서 설계 사이클 타임을 획기적으로 줄여주는 강력한 도구가 될 것입니다.
본론 1: 1자유도 시스템의 지배 방정식
질량 \( m \), 스프링 상수 \( k \), 감쇠 계수 \( c \)로 구성된 1자유도 감쇠 진동 시스템의 운동 방정식은 뉴턴의 제2법칙에 의해 다음과 같이 정의됩니다.
$$ m\ddot{x} + c\dot{x} + kx = F(t) $$
여기서 \( x \)는 변위, \( \dot{x} \)는 속도, \( \ddot{x} \)는 가속도를 나타냅니다. 우리는 이 2차 미분 방정식을 1차 미분 방정식들의 시스템으로 변환하여, 파이썬의 scipy.integrate.odeint를 통해 수치적으로 해를 구할 수 있습니다.
상태 변수를 \( y_1 = x \), \( y_2 = \dot{x} \)라고 정의하면 다음의 시스템을 얻습니다.
$$ \dot{y}_1 = y_2 $$
$$ \dot{y}_2 = \frac{1}{m} (F(t) – c y_2 – k y_1) $$
본론 2: 파이썬 기반 시뮬레이션 구현
아래 코드는 외부 구동력 \( F(t) = A \sin(\omega t) \)가 인가되는 상황을 가정하여 시스템의 시간 응답을 계산하고 시각화하는 예제입니다.
import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint # 시스템 파라미터 정의 m = 1.5 # 피스톤 등가 질량 [kg] k = 20000 # 스프링 강성 [N/m] c = 5.0 # 감쇠 계수 [N·s/m] # 외부 구동력 함수: F(t) = A * sin(omega * t) def force(t): return 10.0 * np.sin(2 * np.pi * 5.0 * t) # 시스템 미분 방정식 (상태 공간 모델) def system_dynamics(y, t, m, c, k): x, v = y dxdt = v dvdt = (force(t) - c * v - k * x) / m return [dxdt, dvdt] # 시뮬레이션 타임 스텝 설정 t = np.linspace(0, 2, 1000) y0 = [0.0, 0.0] # 초기 변위와 초기 속도 # 수치 적분 수행 solution = odeint(system_dynamics, y0, t, args=(m, c, k)) # 결과 시각화 plt.figure(figsize=(10, 5)) plt.plot(t, solution[:, 0], label='Displacement (x)') plt.title('1-DOF Linear Compressor Model Simulation') plt.xlabel('Time [s]') plt.ylabel('Amplitude [m]') plt.grid(True) plt.legend() plt.show()
본론 3: 해석 결과의 공학적 해석
위 코드를 실행하면 시스템이 과도 응답(Transient Response) 단계를 거쳐 정상 상태(Steady State)에 도달하는 과정을 확인할 수 있습니다. 여기서 중요한 점은 다음 두 가지입니다.
- 위상 지연(Phase Lag): 댐핑 계수 \( c \)가 커질수록 구동력과 변위 사이의 위상차는 90도에 가까워집니다. 이는 압축기 밸브 개폐 타이밍과 직결되는 핵심 설계 요소입니다.
- 진폭 제어: 고유 진동수 \( \omega_n = \sqrt{k/m} \) 부근에서 외부 구동 주파수와 일치시키면 진폭이 극대화됩니다. 리니어 압축기는 이 공진점을 활용하여 최소 에너지로 최대 행정을 확보합니다.
결론: 실무 엔지니어를 위한 제언
위와 같은 1D 시뮬레이션은 설계의 방향성을 결정하는 ‘나침반’과 같습니다. 복잡한 3D FEA로 넘어가기 전, 이러한 간소화된 물리 모델을 통해 파라미터 민감도 분석(Sensitivity Analysis)을 수행하십시오. 파이썬은 단순한 언어를 넘어, 여러분의 공학적 직관을 수치적 근거로 뒷받침해 줄 강력한 엔지니어링 동료가 될 것입니다.