Case

解析事例

クリーンルームの室圧制御解析

クリーンルームの室圧制御解析サムネイル画像
はじめに

 クリーンルームでは、細菌の拡散などを防ぐために周囲の部屋との圧力差を一定に保つことが重要です。簡単のため、以下の2部屋のモデルについて、目標の圧力差となる吸気量(qin)を1DCAE解析により求めた事例を紹介します。 

  • twin1
基礎式

 理想気体の状態方程式を時間微分します。(1)式に示します。

 
 
\[ \frac{dP}{dt} = \frac{RT}{V}\frac{dW}{dt} \quad\quad\quad\quad\quad (1) \]

 

 ここで、各変数は以下を示します。

  \begin{array}{ll} P & : 室内の圧力[Pa] \\ V & : 部屋の体積[m3] \\ W & : 室内の空気質量[kg] \\ R & : 空気の気体定数[N m/(kg K)] \\ T & : 室内の絶対温度[K] \\ t & : 時間[sec] \end{array}  

 

 右辺の室内空気量の変化量は(2)式で表されます。 

 

\[ \frac{dW}{dt} = \frac{\rho}{3600} \left(q_{in}-q_{out}-q_{lk}(P)\right) \quad\quad\quad\quad\quad (2) \]

 漏れ流量は(3)式となります。 

\[ q_{lk}(P)=3780\cdot ELA\cdot \left|P\right|^{0.65} \quad\quad\quad\quad\quad (3) \]

 

 ここで、各変数は以下を示します。

 \begin{array}{ll} p & : 空気密度[kg/m3] \\ q_{in} & : 給気風量[m3/h] \\ q_{out} & : 排気風量[m3/h] \\ q_{lk} & : 漏れ風量[m3/h] \\ ELA & : 有効漏れ面積[m2] \end{array} 

 また、ドアの開閉を考慮します。ドアの開閉は次式でモデル化します。 

\[ q_{12}(t)=state(t)*\alpha(P_1-P_2) \quad\quad\quad\quad\quad (4) \]

 

 ここで、state関数は以下を示します。 

\[ state(t)= \begin{cases} 0 & \text{close} \\ 1 & \text{open} \end{cases} \]
\[ \alpha : \text{係数} \]

 

 以上のことから、(5)の室圧と風量の関係式が得られます。この(5)式について1DCAEのモデルを構築し、解析を行います。 

\[ \frac{dP_1}{dt}=\frac{RT\rho}{3600V_1}\left(q_{in}-q_{out1}-q_{lk1}-q_{12}\right)\quad\quad\quad\quad\quad (5) \]
解析モデル作成

まず、基礎式を解くブロック線図モデルを作成します。ソフトウェアはAltair Twin Activateを使用しました。

  • twin2
  • room内
    room内
  • twin4
    room1,room2内
計算結果1

表に示す入力パラメータを使用し、計算を行った結果を示します。

qin[m³/h] 80
qout1[m³/h] 10
qout2[m³/h] 40
R[N m/(kg K)]] 287
T[K] 273.15+20
rho[kg/m³] 1.2
Vrm[m³] 60
ELA[m²] 0.002
alpha[-] 0.05

 

圧力差は、以下に示す結果となりました。

  • twin5
閉鎖時圧力差[Pa] 17.7
開放時圧力差[Pa] 16.7

 

制御システムの追加

 続いて、PID制御を追加し、目的の圧力差となる吸気量を求めます。 

  • twin6
計算結果2

表に示す入力パラメータを使用し、計算を行った結果を示します。目標の圧力差(dP_target)を20Paとしました。

 

dP_target[Pa] 20
P 10
I 3.5
D 5
qout1[m³/h] 10
qout2[m³/h] 40
R[N m/(kg K)]] 287
T[K] 273.15+20
rho[kg/m³] 1.2
Vrm[m³] 60
ELA[m²] 0.002
alpha[-] 0.05

 

圧力差と給気量qinの結果を以下に示します。

 

  • twin7
  • twin8

 

まとめ
  • 2部屋の簡易モデルについて1DCAEモデルを構築し、システム全体の挙動を効率的に解析しました。

  • PID制御を追加し、目的の圧力差となる給気量を求めることができました。このように、制御系モデルとの連携が容易に可能です。

  • 3次元の熱流体解析と比べ計算が高速で、わずか数秒程度で結果が得られます。このため、要求仕様を素早く検証でき、設計段階での迅速なフィードバックが可能になります。

 

 

 

 

RELATED

関連記事はこちら

Altair Twin Activateを利用したロードリミッタの事例
2026.06.04

Altair Twin Activateを利用したロードリミッタの事例

LS-DYNAとTwin Activateを組み合わせて衝突時のシートベルトのロードリミッタ制御効果を確認した事例です。
ドローンの制御事例:Twin ActivateとMotionSolveの連携
2025.09.29

ドローンの制御事例:Twin ActivateとMotionSolveの連携

テラバイトの「ドローンの制御事例:Twin ActivateとMotionSolveの連携 │ Altair Twin ...
Contact

ご相談・資料の請求はこちら

製品やサービスに関するご相談、導入に関するご質問等お気軽にお問い合わせください。
各種資料をご希望の方は、資料請求フォームよりお申し込みいただけます。