계단 교란은 지속적인 변화로 모델링되는 반면, 임펄스 교란은 집중된 단기 입력으로 모델링됩니다. 크기 (A)의 계단은 (u_d(t)=A,\mathcal{H}(t))로 표기되며, 여기서 (\mathcal{H}(t))는 단위 계단 함수이고 라플라스 변환은 (U_d(s)=A/s)입니다. 임펄스는 (u_d(t)=B,\delta(t))로 표기되며, 여기서 (B)는 면적 또는 총 주입량이고 라플라스 변환은 (U_d(s)=B)입니다. 단위 임펄스는 (B=1)입니다.

계단 입력은 파일럿 플랜트가 한 운전 조건에서 다른 조건으로 어떻게 이동하는지 테스트하며, 임펄스 입력은 시간에 따른 플랜트의 동적 응답 분포를 보여줍니다. 공정 모델 (G(s))이 있을 때, 대응하는 출력은 (Y(s)=G(s)U_d(s))에서 구할 수 있습니다.

파일럿 플랜트의 교란 모델링

계단 입력은 지속적인 운전 변화를 나타냄

크기 (A)의 계단 교란은 다음과 같이 모델링됩니다.

[ d(t)=A,\mathcal{H}(t) ]

여기서

[ \mathcal{H}(t)= \begin{cases} 0, & t<0\ 1, & t\geq 0 \end{cases} ]

라플라스 변환은 다음과 같습니다.

[ D(s)=\mathcal{L}{A,\mathcal{H}(t)}=\frac{A}{s}. ]

이는 새로운 밸브 위치, 공급 유량, 가열 부하 또는 제어기 설정값과 같이 효과가 지속되는 갑작스러운 변화를 나타냅니다.

임펄스 입력은 집중된 사건을 나타냄

임펄스 교란은 다음과 같이 모델링됩니다.

[ d(t)=B,\delta(t), ]

여기서 (\delta(t))는 디락 델타 함수(Dirac delta function)이고 (B)는 임펄스 면적입니다.

[ \int_{-\infty}^{\infty} B,\delta(t),dt=B. ]

라플라스 변환은 다음과 같습니다.

[ D(s)=\mathcal{L}{B,\delta(t)}=B. ]

단위 임펄스의 경우 (B=1)이므로 변환은 단순히 다음과 같습니다.

[ D(s)=1. ]

파일럿 플랜트에서 이상적인 임펄스는 공정 응답 시간에 비해 지속 시간이 무시할 수 있는 매우 빠른 추적자, 촉매 또는 반응물 주입을 근사합니다.

과도 응답 계산

공정 모델은 입력을 출력과 연결

전달 함수를 가진 선형화된 공정 모델

[ G(s)=\frac{Y(s)}{D(s)}, ]

에서 출력 변환은 다음과 같습니다.

[ Y(s)=G(s)D(s). ]

시간 영역의 과도 응답은 (Y(s))의 역라플라스 변환을 통해 얻습니다.

이 프레임워크를 사용하면 기본 공정 동역학을 변경하지 않고도 동일한 플랜트 모델을 다양한 교란 유형에 대해 테스트할 수 있습니다.

계단 교란에 대한 응답

크기 (A)의 계단에 대해,

[ Y_{\text{step}}(s)=G(s)\frac{A}{s}. ]

(1/s) 인자는 교란이 변화 후에도 새로운 값으로 유지된다는 사실을 설명합니다.

1차 공정의 경우,

[ G(s)=\frac{K}{\tau s+1}, ]

여기서 (K)는 정상 상태 이득이고 (\tau)는 시상수입니다. 계단 응답은 다음과 같습니다.

[ Y_{\text{step}}(s)

\frac{KA}{s(\tau s+1)}. ]

시간 영역 형태는 다음과 같습니다.

[ y_{\text{step}}(t)=KA\left(1-e^{-t/\tau}\right), \qquad t\geq 0. ]

최종 출력 변화는 (KA)이며, 시상수는 플랜트가 새로운 정상 상태에 얼마나 빨리 도달하는지를 나타냅니다.

임펄스 교란에 대한 응답

면적 (B)의 임펄스에 대해,

[ Y_{\text{impulse}}(s)=G(s)B. ]

단위 임펄스의 경우, 이는 다음과 같이 축소됩니다.

[ Y_{\text{impulse}}(s)=G(s). ]

따라서 임펄스 응답은 공정 전달 함수 자체의 역라플라스 변환입니다.

동일한 1차 공정에 대해,

[ Y_{\text{impulse}}(s)=\frac{KB}{\tau s+1}, ]

따라서

[ y_{\text{impulse}}(t)=\frac{KB}{\tau}e^{-t/\tau}, \qquad t\geq 0. ]

계단 응답과 달리 임펄스 응답은 일반적으로 교란이 지나간 후 0으로 돌아갑니다.

각 테스트가 밝혀내는 것

계단 테스트는 이득과 시상수를 추정

계단 테스트는 다음을 식별하는 데 특히 유용합니다.

정상 상태 이득: 최종 출력 변화량을 입력 변화량으로 나눈 값.
시상수: 과도 응답의 특성 속도.
공정 차수: 응답이 1차, 2차 또는 더 복잡한 시스템처럼 거동하는지 여부.
안정성: 출력이 정착하는지, 진동하는지, 발산하는지 여부.

1차 응답의 경우, 출력은 한 시상수 후에 최종 변화량의 약 (63.2%)에 도달합니다.

임펄스 테스트는 체류 시간 거동을 밝힘

흐름 시스템에서 임펄스 추적자 테스트는 하류 지점에서 농도 응답을 생성합니다. 정규화된 응답은 일반적으로 체류 시간 분포(RTD)로 해석됩니다.

(C(t))가 측정된 추적자 농도인 경우, 정규화된 RTD는 종종 다음과 같이 표기됩니다.

[ E(t)=\frac{C(t)}{\int_0^\infty C(t),dt}. ]

함수 (E(t))는 각 체류 시간에서 시스템을 떠나는 물질의 분율을 설명합니다. 그 적분값은 1입니다.

[ \int_0^\infty E(t),dt=1. ]

해당 누적 RTD는 다음과 같습니다.

[ F(t)=\int_0^t E(\lambda),d\lambda. ]

따라서 임펄스 실험은 분포 (E(t))를 제공하고, 계단 추적자 실험은 그 누적 형태 (F(t))를 제공합니다.

계단 및 임펄스 응답은 수학적으로 관련됨

인과적 선형 시스템의 경우, 단위 계단 응답은 임펄스 응답의 적분입니다.

[ g_{\text{step}}(t)=\int_0^t g_{\text{impulse}}(\lambda),d\lambda. ]

반대로, 임펄스 응답은 계단 응답의 시간 미분입니다.

[ g_{\text{impulse}}(t)=\frac{d}{dt}g_{\text{step}}(t). ]

이 관계는 계단 테스트가 공정이 어떻게 새로운 조건으로 축적되는지를 식별하는 반면, 임펄스 테스트는 공정의 내부 동적 응답의 타이밍과 형태를 노출하는 이유를 설명합니다.

사시간(Dead Time) 및 지연 표현

운송 지연은 응답을 이동시킴

파일럿 플랜트에는 일반적으로 배관, 유체 체류 시간, 센서 배치 및 분석기 응답으로 인한 운송 지연이 포함됩니다. (\theta) 단위의 순수 지연은 라플라스 영역에서 다음과 같이 표현됩니다.

[ e^{-\theta s}. ]

따라서 1차 사시간 모델(first-order-plus-dead-time model)은 다음과 같이 작성할 수 있습니다.

[ G(s)=\frac{K e^{-\theta s}}{\tau s+1}. ]

지연된 계단 입력은 다음 변환을 가집니다.

[ D_{\text{step}}(s)=\frac{A e^{-\theta s}}{s}, ]

그리고 지연된 임펄스 응답은 다음 변환을 가집니다.

[ D_{\text{impulse}}(s)=B e^{-\theta s}. ]

시간 영역에서 응답은 (t=\theta) 이전에는 0이며, 그 이후에 공정 과도가 시작됩니다.

지연은 제어 평가에 영향을 미침

사시간은 공정 출력이 변화를 반영하기 전에 제어기가 대응할 수 있는 능력을 감소시킵니다. 이를 무시하면 지나치게 공격적인 PID 설정이 생성될 수 있으며, 시뮬레이션된 제어 루프가 실제 파일럿 플랜트보다 더 안정적인 것처럼 보이게 만들 수 있습니다.

(e^{-\theta s})를 포함하면 수학적 모델이 실제 측정 및 운송 거동을 더 가깝게 표현하게 됩니다.

트레이드오프 이해

이상적인 입력은 근사치임

물리적 밸브는 무한히 빠르게 위치를 변경할 수 없으며, 추적자 주입은 항상 유한한 지속 시간을 가집니다. 따라서 실제 테스트는 빠른 밸브 움직임이나 짧은 주입 펄스를 사용하여 이상적인 계단 및 임펄스를 근사합니다.

이 근사는 입력 지속 시간이 지배적인 공정 시간 척도에 비해 작을 때만 유용합니다.

계단 테스트는 운전 지점을 교란할 수 있음

큰 계단은 온도, 압력, 조성 또는 유량을 안전하거나 선형적인 운전 범위를 벗어나게 할 수 있습니다. 또한 실험 중에 공정 이득이나 시상수를 변경할 수 있습니다.

테스트 크기는 측정 노이즈를 상회할 만큼 충분히 커야 하지만, 선형 거동을 대략적으로 유지할 만큼 충분히 작아야 합니다.

임펄스 테스트는 측정 품질에 민감함

주입량이 너무 적거나, 센서가 노이즈가 많거나, 주입이 잘 특성화되지 않은 경우 임펄스 응답을 해석하기 어려울 수 있습니다. 유한 폭 펄스는 이상적인 델타 함수가 아닌 실제 입력으로 처리해야 할 수도 있습니다.

RTD 작업의 경우, 응답의 적분이 감지된 총 추적자 양을 나타내므로 측정된 농도도 일관되게 정규화되어야 합니다.

선형 모델은 범위가 제한됨

관계식

[ Y(s)=G(s)D(s) ]

은 일반적으로 선택된 운전 지점 주변에서 플랜트를 선형화하여 얻은 선형 시불변 모델을 가정합니다. 강한 비선형성, 상 변화, 포화, 오염 및 조성 의존적 특성은 비선형 동적 시뮬레이션이나 여러 운전 조건에서의 테스트가 필요할 수 있습니다.

프로젝트에 적용하는 방법

질문하는 물리적 문제와 일치하는 교란 모델을 사용하십시오.

주요 초점이 정상 상태 이득과 공정 속도인 경우: 작은 계단 교란을 적용하고, 전이 곡선을 피팅하며, 이득, 시상수, 공정 차수 및 가능한 사시간을 추정하십시오.
주요 초점이 체류 시간 분포인 경우: 짧은 추적자 펄스를 적용하고, 하류 농도 이력을 측정하며, 정규화하여 (E(t))를 얻으십시오.
주요 초점이 제어기 튜닝인 경우: PID 매개변수를 선택하기 전에 공정 이득, 시상수 및 운송 지연을 포함하는 계단 응답 모델을 사용하십시오.
주요 초점이 빠른 주입 거동인 경우: 단위 임펄스를 자동으로 가정하는 대신 측정된 주입량 또는 펄스 면적인 (B)를 사용하여 교란을 (B\delta(t))로 모델링하십시오.
주요 초점이 모델 검증인 경우: 입력 지속 시간, 센서 지연 및 운전 범위가 선형 근사를 정당화하는지 확인하면서 예측된 곡선과 측정된 과도 곡선을 비교하십시오.

계단 입력은 운전 상태 간의 플랜트 전이를 특성화하고, 임펄스 입력은 플랜트가 시간에 따라 교란을 어떻게 분산하고 소산시키는지 보여줍니다.

요약 표:

교란	수학적 모델	평가에 가장 적합	주요 결과
계단	`d(t) = A H(t)`; `D(s) = A/s`	이득, 시상수, 공정 차수, 안정성 및 제어기 튜닝	새로운 운전 조건으로의 전이
임펄스	`d(t) = B delta(t)`; `D(s) = B`	체류 시간 분포 및 동적 응답 형태	짧은 사건의 타이밍과 소산
1차 계단 응답	`y(t) = KA(1 - e^(-t/tau))`	정상 상태 변화 및 공정 속도	최종 변화는 `KA`; `t = tau`에서 63.2% 도달
1차 임펄스 응답	`y(t) = (KB/tau)e^(-t/tau)`	단기 주입 거동	응답이 0으로 감쇠
사시간 모델	`G(s) = K e^(-theta s)/(tau s + 1)`	운송 지연 및 현실적인 제어 평가	지연 `theta` 후 응답 시작
추적자 임펄스 테스트	`E(t) = C(t)/integral C(t)dt`	흐름 시스템의 체류 시간 분포	물질 체류 시간 프로필