JP2019049283A - Gas charging device and gas charging method - Google Patents

Gas charging device and gas charging method Download PDF

Info

Publication number
JP2019049283A
JP2019049283A JP2017173016A JP2017173016A JP2019049283A JP 2019049283 A JP2019049283 A JP 2019049283A JP 2017173016 A JP2017173016 A JP 2017173016A JP 2017173016 A JP2017173016 A JP 2017173016A JP 2019049283 A JP2019049283 A JP 2019049283A
Authority
JP
Japan
Prior art keywords
pressure
gas
filling
flow rate
pdot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2017173016A
Other languages
Japanese (ja)
Inventor
大 牛島
Dai Ushijima
大 牛島
ジョナタン・シュバリエ
Chevalier Jonathan
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Azbil Corp
Original Assignee
Azbil Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Azbil Corp filed Critical Azbil Corp
Priority to JP2017173016A priority Critical patent/JP2019049283A/en
Publication of JP2019049283A publication Critical patent/JP2019049283A/en
Pending legal-status Critical Current

Links

Landscapes

  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

To surely charge a gas into a pressure vessel at a high speed.SOLUTION: A control device 9 is provided with a pressure fluctuation rate calculation portion 9-1, a mass flow rate calculation portion 9-2, and a volume flow rate calculation portion 9-3. The pressure fluctuation rate calculation portion 9-1 acquires a charging pressure P (Pmes(measured pressure)) of a gas in a pressure vessel 7 at a flow rate update interval Δt, and calculates a pressure fluctuation rate Pdot of the charging pressure P to make a change curve of the charging pressure P of the gas approximate an ideal pressure curve as a control target in controlling the charging amount of the gas in each acquirement of the charging pressure P of the gas. The mass flow rate calculation portion 9-2 determines a mass flow rate mdot of the gas from the pressure fluctuation rate Pdot calculated by the pressure fluctuation calculation portion 9-1. The volume flow rate calculation portion 9-3 determines a volume flow rate Q of the gas from the mass flow rate mdot determined by the mass flow rate calculation portion 9-2. The control device 9 controls an operation of the charging device 11 so that the gas of the volume flow rate Q calculated by the volume flow rate calculation portion 9-3 is delivered to the pressure vessel 7.SELECTED DRAWING: Figure 2

Description

本発明は、圧力容器にガスを充填するガス充填装置およびガス充填方法に関する。   The present invention relates to a gas filling apparatus and a gas filling method for filling a pressure vessel with gas.

従来より、圧力容器にガスを充填するガス充填装置として、例えば特許文献1に示されるような圧力発生充填装置がある。   Conventionally, as a gas filling apparatus for filling a pressure vessel with gas, there is a pressure generating and filling apparatus as shown in, for example, Patent Document 1.

この圧力発生充填装置では、流量コントローラを3台並列に接続し、この流量コントローラを切り替えることによって、3段階に分けて圧力容器へのガスの充填を行うようにしている。   In this pressure generating and filling apparatus, three flow controllers are connected in parallel, and by switching the flow controllers, the pressure vessel is filled with gas in three stages.

各流量コントローラは、例えば、第1の流量コントローラを流量大、第2の流量コントローラを流量中、第3の流量コントローラを流量小というように、その流量の制御範囲が異なっている。各流量コントローラは、制御装置によって、そこを流れるガスの流量が一定となるようにコントロールされる。   Each flow rate controller has different control ranges of the flow rate, for example, the first flow rate controller has a large flow rate, the second flow rate controller has a flow rate, and the third flow rate controller has a small flow rate. Each flow rate controller is controlled by the controller so that the flow rate of gas flowing therethrough becomes constant.

図15に、特許文献1に示された圧力容器へのガスの充填圧力の変化曲線を示す。図15において、Pは圧力容器内のガスの充填圧力、P0は圧力容器内のガスの充填圧力の初期値(初期圧力)である。P1,P2,P3は1段目,2段目,3段目の圧力設定値である。t1,t2,t3は、ガスの充填圧力Pが圧力設定値P1,P2,P3に到達した時刻である。圧力設定値P3は最終目標圧力Pfに相当する。   The change curve of the filling pressure of the gas to the pressure vessel shown in FIG. 15 by patent document 1 is shown. In FIG. 15, P is the gas filling pressure in the pressure vessel, and P0 is the initial value (initial pressure) of the gas filling pressure in the pressure vessel. P1, P2 and P3 are pressure setting values for the first, second and third stages. t1, t2 and t3 are times when the gas filling pressure P reaches the pressure set values P1, P2 and P3. The pressure set value P3 corresponds to the final target pressure Pf.

この例では、圧力容器内のガスの充填圧力PがP0からP1に達するまでは、第1の流量コントローラを稼働して圧力容器へのガスの充填を行う。ガスの充填圧力PがP1に達すると、第2の流量コントローラを稼働して圧力容器へのガスの充填を行う。ガスの充填圧力PがP2に達すると、第3の流量コントローラを稼働して圧力容器へのガスの充填を行う。そして、ガスの充填圧力Pが最終目標圧力Pfに達した後は、この最終目標圧力Pfの状態を保つようにする。すなわち、ガスの充填圧力Pを最終目標圧力Pfに整定させる。   In this example, the first flow rate controller is operated to fill the pressure vessel with gas until the filling pressure P of the gas in the pressure vessel reaches P0 to P1. When the gas filling pressure P reaches P1, the second flow rate controller is operated to fill the pressure vessel with gas. When the gas filling pressure P reaches P2, the third flow rate controller is operated to fill the pressure vessel with gas. Then, after the gas filling pressure P reaches the final target pressure Pf, the final target pressure Pf is maintained. That is, the gas filling pressure P is settled to the final target pressure Pf.

特開2015−113938号公報JP, 2015-113938, A

しかしながら、この特許文献1に示された圧力発生充填装置では、図16に示すように、流量コントローラの切替時にガスの充填圧力Pがオーバーシュートし、ガスの充填圧力Pが最終目標圧力Pfに整定されるまでの時間が長引くという問題があった。すなわち、圧力容器へのガスの充填を開始してから、ガスの充填が完了(整定)するまでの充填時間が長引くという問題があった。   However, in the pressure generating and filling apparatus shown in this patent document 1, as shown in FIG. 16, the gas filling pressure P overshoots when the flow rate controller is switched, and the gas filling pressure P settles to the final target pressure Pf. There was a problem that time until it was done was prolonged. That is, there has been a problem that the filling time from the start of gas filling to the pressure vessel to the completion (setting) of the gas filling is prolonged.

本発明は、このような課題を解決するためになされたもので、その目的とするところは、高速かつ正確に圧力容器へのガスの充填を行うことが可能なガス充填装置およびガス充填方法を提供することにある。   The present invention has been made to solve such problems, and its object is to provide a gas filling apparatus and a gas filling method capable of filling a pressure vessel with gas at high speed and accurately. It is to provide.

このような目的を達成するために本発明は、圧力容器(7)にガスを充填するように構成されたガス充填装置(100)において、圧力容器内のガスの充填圧力Pを計測するように構成された圧力計測部(8)と、圧力計測部が計測するガスの充填圧力Pを流量更新間隔Δtで取得し、その取得したガスの充填圧力Pに基づいて圧力容器へのガスの充填量を制御し、ガスの充填圧力Pを最終目標圧力Pfに整定させるように構成された制御装置(9)とを備え、制御装置(9)は、ガスの充填圧力Pを取得する毎に、ガスの充填量を制御する際の制御目標として、ガスの充填圧力Pの変化曲線を下記(A)式で規定される理想圧力曲線に近似させるような充填圧力Pの圧力変動率Pdotを算出するように構成された圧力変動率算出部(9−1)を備えることを特徴とする。   In order to achieve such an object, the present invention measures the filling pressure P of the gas in the pressure vessel in a gas filling apparatus (100) configured to fill the pressure vessel (7) with gas. The pressure measurement unit (8) configured and the filling pressure P of the gas measured by the pressure measurement unit are acquired at the flow rate update interval Δt, and the filling amount of the gas into the pressure vessel based on the acquired filling pressure P of the gas And a controller (9) configured to control the filling pressure P of the gas to the final target pressure Pf, and the control device (9) controls the gas filling pressure P every time the gas filling pressure P is obtained. As a control target at the time of controlling the filling amount of P, calculate the pressure fluctuation rate Pdot of the filling pressure P that approximates the change curve of the filling pressure P of the gas to the ideal pressure curve defined by the following equation (A) Pressure fluctuation rate calculation unit (9-1) Characterized in that it comprises.

この発明において、制御装置は、ガスの充填圧力Pを取得する毎に、ガスの充填圧力Pの変化曲線を(A)式で規定される理想圧力曲線に近似させるような圧力変動率Pdotを算出し、この算出した圧力変動率Pdotを制御目標として圧力容器へのガスの充填量を制御する。これにより、ガスの充填圧力Pを取得する毎に、ガスの充填圧力Pが(A)式で規定される理想圧力曲線に近似して変化するものとなり、オーバーシュートを生じさせないようにして、高速かつ正確に圧力容器へのガスの充填を行うことが可能となる。   In the present invention, the control device calculates the pressure fluctuation rate Pdot which approximates the change curve of the gas filling pressure P to the ideal pressure curve defined by the equation (A) every time the gas filling pressure P is obtained. The amount of gas filling the pressure vessel is controlled with the calculated pressure fluctuation rate Pdot as a control target. As a result, whenever the gas filling pressure P is obtained, the gas filling pressure P changes in approximation to the ideal pressure curve defined by the equation (A), so that overshoot does not occur, and high speed is achieved. And it becomes possible to carry out the filling of the gas to a pressure vessel correctly.

本発明では、ガスの充填圧力Pを取得する毎に圧力変動率Pdotを算出するが、この圧力変動率Pdotの算出式として下記(B)式や(C)式、(D)式などを用いることが考えられる。
本発明では、(C)式を用いて圧力変動率Pdotを算出する手法を技法1(局所微分)と呼び、(D)式を用いて圧力変動率Pdotを算出する手法を技法2(中間点微分)と呼び、(B)式を用いて圧力変動率Pdotを算出する手法を技法3(線形補完)と呼ぶものとする。技法1,2には後述するように短所があり、技法1,2,3のうち技法3が最も適した技法であると言える。アプリケーションによっては技法1や2を用いてもよく、技法3に技法1や技法2を組み合わせるなどしてもよい。
In the present invention, the pressure fluctuation rate Pdot is calculated each time the filling pressure P of the gas is obtained, but the following equations (B), (C), (D), etc. are used as the calculation formula of the pressure fluctuation rate Pdot. It is conceivable.
In the present invention, a method of calculating the pressure fluctuation rate Pdot using the equation (C) is called technique 1 (local differentiation), and a method of calculating the pressure fluctuation rate Pdot using the equation (D) is technique 2 (midpoint A method of calculating the pressure fluctuation rate Pdot by using the equation (B) is referred to as a technique 3 (linear interpolation). Techniques 1 and 2 have disadvantages as will be described later, and it can be said that technique 3 is the most suitable of techniques 1, 2 and 3. Depending on the application, techniques 1 and 2 may be used, and technique 3 may be combined with techniques 1 and 2.

なお、上記説明では、一例として、発明の構成要素に対応する図面上の構成要素を、括弧を付した参照符号によって示している。   In the above description, as an example, constituent elements on the drawing corresponding to constituent elements of the invention are indicated by reference numerals in parentheses.

以上説明したことにより、本発明によれば、ガスの充填圧力Pを取得する毎に、ガスの充填圧力Pの変化曲線を(A)式で規定される理想圧力曲線に近似させるような充填圧力Pの圧力変動率Pdotを算出し、この算出した圧力変動率Pdotを制御目標として圧力容器へのガスの充填量を制御するようにしたので、ガスの充填圧力Pを取得する毎に、ガスの充填圧力Pが(A)式で規定される理想圧力曲線に近似して変化するものとなり、オーバーシュートを生じさせないようにして、高速かつ正確に圧力容器へのガスの充填を行うことが可能となる。   As described above, according to the present invention, the filling pressure is such that the change curve of the filling pressure P of the gas is approximated to the ideal pressure curve defined by the equation (A) every time the filling pressure P of the gas is obtained. Since the pressure fluctuation rate Pdot of P is calculated and the calculated pressure fluctuation rate Pdot is used as a control target to control the filling amount of gas into the pressure vessel, the gas filling pressure P is obtained each time the gas filling pressure P is acquired. The filling pressure P changes close to the ideal pressure curve defined by the equation (A), and it is possible to fill the pressure vessel with gas at high speed and accurately without causing overshoot. Become.

図1は、本発明の実施の形態に係るガス充填装置の要部の構成を示す図である。FIG. 1 is a view showing the configuration of the main part of a gas filling apparatus according to an embodiment of the present invention. 図2は、このガス充填装置における制御装置の機能ブロック図である。FIG. 2 is a functional block diagram of a control device in this gas filling device. 図3は、理想的な充填プロセスにおける充填圧力の変化を示す図である。FIG. 3 is a diagram showing changes in filling pressure in an ideal filling process. 図4は、初期圧力をP0とした場合の最終目標圧力Pfまでの理想圧力曲線を示す図である。FIG. 4 is a diagram showing an ideal pressure curve up to the final target pressure Pf when the initial pressure is P0. 図5は、時定数τが異なるガスの充填圧力の変化曲線を例示する図である。FIG. 5 is a diagram illustrating a change curve of the filling pressure of the gas having different time constants τ. 図6は、技法1(局所微分)を説明する図である。FIG. 6 is a diagram for explaining the technique 1 (local differentiation). 図7は、理想圧力曲線と対比して、技法1を用いた場合のガスの充填圧力の変化曲線を示す図である。FIG. 7 is a diagram showing a change curve of the gas filling pressure when the technique 1 is used, in contrast to the ideal pressure curve. 図8は、技法2(中間点微分)を説明する図である。FIG. 8 is a diagram for explaining the technique 2 (mid-point differentiation). 図9は、理想圧力曲線と対比して、技法2を用いた場合のガスの充填圧力の変化曲線を示す図である。FIG. 9 is a diagram showing a change curve of the gas filling pressure in the case of using the technique 2 in contrast to the ideal pressure curve. 図10は、技法3(線形補間)を説明する図である。FIG. 10 is a diagram for explaining the technique 3 (linear interpolation). 図11は、理想圧力曲線と対比して、技法3を用いた場合のガスの充填圧力の変化曲線を示す図である。FIG. 11 is a diagram showing a change curve of the gas filling pressure when the technique 3 is used, in contrast to the ideal pressure curve. 図12は、目標圧力に対してある割合まで到達した時点で流量の更新を止めて望む圧力まで同じ流量を保つようにした例を説明する図である。FIG. 12 is a diagram for explaining an example in which updating of the flow rate is stopped when a certain percentage of the target pressure is reached, and the same flow rate is maintained to the desired pressure. 図13は、連続する圧力変動率Pdotの比を用いて時定数τの下限を見つける方法を説明する図である。FIG. 13 is a diagram for explaining a method of finding the lower limit of the time constant τ by using the ratio of continuous pressure fluctuation rates Pdot. 図14は、理想圧力曲線と対比して、技法1,2,3を用いた場合のガスの充填圧力の変化曲線を示す図である。FIG. 14 is a diagram showing a change curve of the gas filling pressure in the case of using the techniques 1, 2 and 3 in contrast to the ideal pressure curve. 図15は、特許文献1に示されたガスの充填圧力の変化曲線を示す図である。FIG. 15 is a view showing a change curve of the gas filling pressure shown in Patent Document 1. As shown in FIG. 図16は、特許文献1に示されたガスの充填圧力の変化曲線において流量コントローラの切替時にオーバーシュートが生じる様子を示す図である。FIG. 16 is a diagram showing how overshoot occurs when the flow rate controller is switched in the change curve of the gas filling pressure shown in Patent Document 1. As shown in FIG.

〔発明の概要〕
先ず、本発明の実施の形態の説明に入る前に、本発明の概要について説明する。
最初に、静圧(流れがない圧力および安定後の圧力)は、測定された圧力の関数と仮定し、ガス流量、ガス流量の変化、ガス温度によって変化すると考える。これらの因子がすべてゼロである場合、静圧は測定圧力に等しいので、これらの因子のすべてが小さい値では、静圧と測定圧力の差がゼロになると推測できる。
[Summary of the invention]
First, prior to describing the embodiments of the present invention, an outline of the present invention will be described.
First, static pressure (pressure without flow and pressure after stabilization) is assumed to be a function of measured pressure, and it is considered to change with gas flow rate, change of gas flow rate, gas temperature. If all these factors are zero, the static pressure is equal to the measured pressure, so it can be inferred that the difference between the static pressure and the measured pressure is zero at small values of all these factors.

したがって、1次コントローラ(後述)の流量更新を高頻度に行うことで(変化量を少なくする事で)、オーバーシュートおよび安定化時間の問題を解決する事ができる。実際、その指数関数的性質のために、1次コントローラ出力は、圧力差が小さい充填プロセスの終了時にガスに小さな変化を生じさせる。   Therefore, the problems of overshoot and stabilization time can be solved by frequently updating the flow rate of the primary controller (described later) (by reducing the amount of change). In fact, due to its exponential nature, the primary controller output causes small changes in the gas at the end of the filling process with small pressure differences.

本発明では、物理的コントローラの実施を避けるために、最終目標圧力、実際の測定圧力、ガス温度、アプリケーションに依存するいくつかの定数(経験的にも容易に見つけることができる)と同様なコントローラの速度を表す時定数τだけを、1次コントローラに数学的(仮想的な)実装する事を提案する。   In the present invention, a controller similar to the final target pressure, the actual measured pressure, the gas temperature, some application dependent constants (which can be easily found empirically) in order to avoid the implementation of the physical controller We propose to implement mathematically (virtually) in the first-order controller only for the time constant τ that represents the velocity of.

これらすべては、どのプログラミング言語でも迅速にプログラミングできる。最大の利点は、時定数τを最適化する事で充填プロセスを望むように高速かつ正確にする事が可能となることである。   All these can be programmed quickly in any programming language. The biggest advantage is that optimizing the time constant τ allows the filling process to be as fast and accurate as desired.

〔実施の形態〕
以下、本発明の実施の形態を図面に基づいて詳細に説明する。図1は本発明の実施の形態に係るガス充填装置100の要部の構成を示す図である。
Embodiment
Hereinafter, embodiments of the present invention will be described in detail based on the drawings. FIG. 1 is a view showing the configuration of the main part of a gas filling apparatus 100 according to an embodiment of the present invention.

このガス充填装置100は、流量コントローラ1−1〜1−3と、流量コントローラ用遮断弁2−1〜2−3と、遮断弁3,4と、ゲート弁5−1,5−2と、真空ポンプ6と、圧力容器7と、圧力計8と、制御装置9と、温度センサ10とを備えている。   The gas filling apparatus 100 includes flow controllers 1-1 to 1-3, shutoff valves 2-1 to 2-3 for a flow controller, shutoff valves 3, 4 and gate valves 5-1, 5-2. A vacuum pump 6, a pressure vessel 7, a pressure gauge 8, a control device 9, and a temperature sensor 10 are provided.

また、ガス充填装置100は、図示してはいないが、圧力容器7へ充填するガス(G)を生成する圧力発生装置を備えている。この圧力発生装置としては、一定圧のガスを生成するガス圧縮ポンプなどを想定している。充填するガスが、空気以外のガスで、例えば窒素のときは、圧力発生装置のガスの取り入れ口に窒素ガスボンベを繋ぐ。   Further, although not shown, the gas filling apparatus 100 is provided with a pressure generator that generates a gas (G) to be filled into the pressure vessel 7. As the pressure generator, a gas compression pump or the like that generates a gas at a constant pressure is assumed. When the gas to be charged is a gas other than air, for example, nitrogen, a nitrogen gas cylinder is connected to the gas inlet of the pressure generator.

このガス充填装置100において、流量コントローラ1−1〜1−3は、圧力容器7へのガスの供給路L中に並列に接続されている。流量コントローラ1−1〜1−3は、例えば、流量コントローラ1−1が流量大、流量コントローラ1−2が流量中、流量コントローラ1−3が流量小というように、その流量の制御範囲が異なっている。この流量コントローラ1−1〜1−3に対して流量コントローラ用遮断弁2−1〜2−3が設けられている。なお、流量を確実に遮断する機能を持つ流量コントローラ1(1−1,1−2,1−3)を用いる場合は、流量コントローラ用遮断弁2(2−1,2−2,2−3)は使用しなくてもよい。   In the gas filling device 100, the flow rate controllers 1-1 to 1-3 are connected in parallel in the gas supply path L to the pressure vessel 7. In the flow controllers 1-1 to 1-3, for example, the control range of the flow is different such that the flow controller 1-1 has a large flow, the flow controller 1-2 has a flow, and the flow controller 1-3 has a small flow. ing. Flow control valve blocking valves 2-1 to 2-3 are provided for the flow controllers 1-1 to 1-3. When using the flow controller 1 (1-1, 1-2, 1-3) having a function to shut off the flow reliably, the flow control valve 2 (2-1, 2-2, 2-3) ) May not be used.

遮断弁3は、圧力容器7へのガスの供給路L中に設けられており、この遮断弁3の前後にゲート弁5−1,5−2を介して真空ポンプ6が接続されている。圧力計8は、圧力容器7内のガスの充填圧力Pを計測し、その計測した充填圧力Pを制御装置9へ送る。   The shutoff valve 3 is provided in the gas supply passage L to the pressure vessel 7, and the vacuum pump 6 is connected to the front and back of the shutoff valve 3 via the gate valves 5-1 and 5-2. The pressure gauge 8 measures the filling pressure P of the gas in the pressure vessel 7, and sends the measured filling pressure P to the controller 9.

圧力計8と圧力容器7との間には、圧力計8を保護するための安全装置として、高圧時に自動的に流路を遮断する遮断弁4が設けられている。温度センサ10は、圧力容器7へのガスの温度Tを検出し、その検出したガスの温度Tを制御装置9へ送る。   Between the pressure gauge 8 and the pressure vessel 7, as a safety device for protecting the pressure gauge 8, a shutoff valve 4 is provided which automatically shuts off the flow path at high pressure. The temperature sensor 10 detects the temperature T of the gas to the pressure vessel 7, and sends the detected temperature T of the gas to the controller 9.

制御装置9は、流量コントローラ1(1−1〜1−3)、流量コントローラ用遮断弁2(2−1〜2−3)、遮断弁3、ゲート弁5(5−1,5−2)、真空ポンプ6などの動作を制御することによって、圧力容器7へのガスの充填を行う。   The controller 9 includes a flow controller 1 (1-1 to 1-3), a shutoff valve 2 (2-1 to 2-3) for the flow controller, a shutoff valve 3, and a gate valve 5 (5-1 and 5-2). The pressure container 7 is filled with gas by controlling the operation of the vacuum pump 6 and the like.

この制御装置9は、プロセッサや記憶装置からなるハードウェアと、これらのハードウェアと協働して各種機能を実現させるプログラムとによって実現される。なお、図1において、制御装置9によってその動作が制御される構成(1点鎖線で囲んだ構成)を充填装置11と呼ぶ。   The control device 9 is realized by hardware including a processor and a storage device, and a program that realizes various functions in cooperation with the hardware. Note that, in FIG. 1, a configuration (a configuration surrounded by an alternate long and short dash line) whose operation is controlled by the control device 9 is referred to as a filling device 11.

図2に制御装置9の機能ブロック図を示す。なお、図2には、制御装置9の機能ブロックと合わせて、この機能ブロックが果たす役割を説明する必要上、充填装置11や圧力容器7も示している。   FIG. 2 shows a functional block diagram of the control device 9. Note that FIG. 2 also shows the filling device 11 and the pressure vessel 7 in order to explain the role played by this functional block, together with the functional block of the control device 9.

制御装置9は、本実施の形態特有の機能ブロックとして、圧力変動率算出部9−1と、質量流量算出部9−2と、体積流量算出部9−3とを備えている。   The control device 9 includes a pressure fluctuation rate calculation unit 9-1, a mass flow rate calculation unit 9-2, and a volumetric flow rate calculation unit 9-3 as functional blocks unique to the present embodiment.

この制御装置9において、圧力変動率算出部9−1は、圧力計8が計測する圧力容器7内のガスの充填圧力P(Pmes(計測圧))を流量更新間隔Δtで取得し、ガスの充填圧力Pを取得する毎に、ガスの充填量を制御する際の制御目標として、ガスの充填圧力Pの変化曲線を理想圧力曲線(後述する(1)式で規定される理想圧力曲線)に近似させるような充填圧力Pの圧力変動率Pdotを算出する。   In the control device 9, the pressure fluctuation rate calculation unit 9-1 acquires the filling pressure P (Pmes (measurement pressure)) of the gas in the pressure vessel 7 measured by the pressure gauge 8 at the flow rate update interval Δt. Every time when the filling pressure P is acquired, a change curve of the filling pressure P of the gas is set as an ideal pressure curve (an ideal pressure curve defined by equation (1) described later) as a control target when controlling the filling amount of gas. The pressure variation rate Pdot of the filling pressure P is calculated so as to approximate.

この圧力変動率算出部9−1が前述した1次コントローラに相当する。また、圧力変動率算出部9−1において制御目標として圧力変動率Pdotを算出する周期が、前述した1次コントローラの流量更新の頻度に相当する。   The pressure fluctuation rate calculation unit 9-1 corresponds to the above-described primary controller. Further, a cycle in which the pressure fluctuation rate Pdot is calculated as a control target in the pressure fluctuation rate calculation unit 9-1 corresponds to the frequency of the flow rate update of the primary controller described above.

質量流量算出部9−2は、圧力変動率算出部9−1によって算出された圧力変動率Pdotからガスの質量流量mdotを求める。体積流量算出部9−3は、質量流量算出部9−2によって求められた質量流量mdotからガスの体積流量Qを求める。   The mass flow rate calculation unit 9-2 obtains the mass flow rate mdot of gas from the pressure variation rate Pdot calculated by the pressure variation rate calculation unit 9-1. The volumetric flow rate calculation unit 9-3 obtains the volumetric flow rate Q of the gas from the mass flow rate mdot obtained by the mass flow rate calculation unit 9-2.

なお、質量流量算出部9−2での質量流量mdotの算出や体積流量算出部9−3での体積流量Qの算出では、ガスの温度Tなども用いられる。ガスの質量流量mdotや体積流量Qの算出に用いられる式などについては後述する。   In the calculation of the mass flow mdot in the mass flow calculation unit 9-2 and the calculation of the volumetric flow Q in the volumetric flow calculation unit 9-3, the temperature T of the gas or the like is also used. The equations used to calculate the mass flow rate mdot of the gas and the volumetric flow rate Q will be described later.

制御装置9は、体積流量算出部9−3で算出されたガスの体積流量Qが圧力容器7に送られるように、充填装置11の動作を制御する。すなわち、圧力計8が計測するガスの充填圧力Pを流量更新間隔Δtで取得し、その取得したガスの充填圧力Pに基づいて圧力容器7へのガスの充填量(流量)を制御し、ガスの充填圧力Pを最終目標圧力(以下、単に目標圧力と呼ぶ場合もある。)Pfに整定させるようにする。   The control device 9 controls the operation of the filling device 11 so that the volumetric flow rate Q of the gas calculated by the volumetric flow rate calculation unit 9-3 is sent to the pressure vessel 7. That is, the filling pressure P of the gas measured by the pressure gauge 8 is acquired at the flow rate update interval Δt, and the filling amount (flow rate) of the gas into the pressure vessel 7 is controlled based on the acquired filling pressure P of the gas. The filling pressure P is set to a final target pressure (hereinafter sometimes referred to simply as the target pressure) Pf.

なお、この例において、制御装置9は、流量コントローラ1−1,1−2,1−3を選択的に稼働し、体積流量算出部9−3で算出された体積流量Qのガスが圧力容器7に送られるように、流量コントローラ1−1,1−2,1−3を流れるガスの流量をコントロールするが、流量コントローラ1−1,1−2,1−3を組み合わせて稼働するようにしてもよい。また、この例では、流量コントローラ1を3台としているが、3台に限られるものではなく、1台としてもよい。   In this example, the control device 9 selectively operates the flow controllers 1-1, 1-2, and 1-3, and the gas having the volumetric flow rate Q calculated by the volumetric flow rate calculating unit 9-3 is a pressure vessel. Control the flow rate of the gas flowing through the flow rate controllers 1-1, 1-2, and 1-3 so as to be sent to 7, but operate in combination with the flow rate controllers 1-1, 1-2, and 1-3. May be Moreover, although three flow controllers 1 are provided in this example, the number of flow controllers 1 is not limited to three, and may be one.

制御装置9は、圧力容器7へのガスの充填を開始する場合、ゲート弁5−2を開き、真空ポンプ6を動かして、圧力容器7の真空排気を行う。次に、ゲート弁5−1を開き、弁2と弁3との間の残圧を真空排気する。これにより、圧力容器7の初期圧力P0を真空圧とする。なお、真空ポンプ6を使わず残圧充填ガスを排出する機構は別なものを採用してもかまわない。   When starting the filling of the pressure vessel 7 with gas, the controller 9 opens the gate valve 5-2 and moves the vacuum pump 6 to evacuate the pressure vessel 7. Next, the gate valve 5-1 is opened, and the residual pressure between the valve 2 and the valve 3 is evacuated. Thereby, the initial pressure P0 of the pressure vessel 7 is set to a vacuum pressure. The mechanism for discharging the residual pressure filling gas without using the vacuum pump 6 may adopt another one.

〔理想圧力曲線について〕
理想的な充填プロセスは、時間とともに圧力の指数関数的な変化になる。すなわち、図3に示すように、充填圧力Pは時間tの経過とともに指数関数的に大きくなり、目標圧力Pfから遠く離れたところでは充填率が高くなり、目標圧力Pfに近づくと充填率は小さくなる。
[On the ideal pressure curve]
An ideal filling process results in an exponential change in pressure over time. That is, as shown in FIG. 3, the filling pressure P increases exponentially with the passage of time t, the filling rate becomes high at a point far away from the target pressure Pf, and the filling rate becomes small when approaching the target pressure Pf. Become.

この理想的な充填プロセスでは、1次制御の時間応答(圧力と時間の関数)として、理想圧力曲線の式が求められる。図4に、初期圧力をP0とした場合の目標圧力Pfまでの理想圧力曲線を示す。下記に、理想圧力曲線の式の算出過程を示す。   In this ideal filling process, the equation of the ideal pressure curve is determined as the time response (function of pressure and time) of primary control. FIG. 4 shows an ideal pressure curve up to the target pressure Pf when the initial pressure is P0. Below, the calculation process of the formula of an ideal pressure curve is shown.

この理想圧力曲線の式(e)において、τは時定数であり、この時定数τが圧力の上昇スピードを決め、時定数τが小さいほどPf−Pで与えられる流量は大きくなる。すなわち、図5に示すように、τが小さいほど圧力の上昇スピードが速くなり、τが大きいほど圧力の上昇スピードが緩くなる。図5において、τ=5の曲線C1は圧力の上昇スピードが速く、その分流量の変動が激しい。τ=45の曲線C2は圧力の上昇スピードが遅く、整定するまでに時間がかかる。τ=20の曲線C3はC1とC2の中間に位置している。例えば、この曲線C3を理想圧力曲線とする。   In equation (e) of this ideal pressure curve, τ is a time constant, and this time constant τ determines the speed at which the pressure rises, and the smaller the time constant τ, the larger the flow rate given by Pf−P. That is, as shown in FIG. 5, the speed of pressure rise is faster as τ is smaller, and the speed of pressure rise is slower as τ is larger. In FIG. 5, the curve C1 with τ = 5 has a fast rise speed of pressure, and the fluctuation of the flow rate is intense accordingly. The curve C2 with τ = 45 has a slow rise in pressure and takes time to settle. A curve C3 of τ = 20 is located between C1 and C2. For example, this curve C3 is an ideal pressure curve.

理想圧力曲線の式(e)において、tiを流量更新間隔Δt毎のガスの充填圧力Pの取得タイミングとし、Piを今回取得した圧力Pとすると、この式(e)は下記(1)式として示される。この(1)式が本発明でいう(A)式に相当する。   Assuming that in the equation (e) of the ideal pressure curve, ti is the acquisition timing of the gas filling pressure P at every flow rate update interval Δt, and Pi is the pressure P acquired this time, this equation (e) is as the following equation (1) Indicated. The equation (1) corresponds to the equation (A) in the present invention.

〔圧力変動率について〕
充填圧力Pの変動率(圧力変動率)は、理想気体に対応する標準体積流量に直線的に関連することができる。したがって、理想圧力曲線上の圧力変動率(微分値)を求めれば、その圧力変動率を制御流量に関連することができる。上記の(1)式で示される理想圧力曲線上の圧力変動率をPdot(t)とすると、この圧力変動率Pdot(t)は下記(2)式で表される。
[About pressure fluctuation rate]
The variation rate of the filling pressure P (pressure variation rate) can be linearly related to the standard volumetric flow rate corresponding to the ideal gas. Therefore, if the pressure fluctuation rate (differential value) on the ideal pressure curve is determined, the pressure fluctuation rate can be related to the control flow rate. Assuming that the pressure fluctuation rate on the ideal pressure curve represented by the above equation (1) is Pdot (t), this pressure fluctuation rate Pdot (t) is represented by the following equation (2).

本実施の形態において、圧力変動率算出部9−1は、圧力計8が計測する圧力容器7内のガスの充填圧力Pを流量更新間隔Δtで取得し、ガスの充填圧力Pを取得する毎に、ガスの充填量を制御する際の制御目標として、ガスの充填圧力Pの変化曲線を上記の(1)式で規定される理想圧力曲線に近似させるような充填圧力Pの圧力変動率Pdotを算出する。この場合、圧力変動率Pdot(制御目標)の算出方法として、次のような技法1、2、3が考えられる。   In the present embodiment, the pressure fluctuation rate calculation unit 9-1 acquires the filling pressure P of the gas in the pressure vessel 7 measured by the pressure gauge 8 at the flow rate update interval Δt, and acquires the filling pressure P of the gas. In addition, as a control target when controlling the gas filling amount, the pressure fluctuation rate Pdot of the filling pressure P that approximates the change curve of the filling pressure P of the gas to the ideal pressure curve defined by the above equation (1) Calculate In this case, the following techniques 1, 2 and 3 can be considered as a method of calculating the pressure fluctuation rate Pdot (control target).

〔技法1:局所微分(図6参照)〕
技法1では、流量更新間隔Δtが短く、上記の(2)式において「t≒ti」であると考える。すなわち、t−ti≒0と考える。そして、t−ti≒0として得られる圧力変動率Pdotを、すなわち下記(3)式で示される圧力変動率Pdotを、今回の圧力の取得ポイント(今回の流量更新タイミング)での制御目標とする。
Technique 1: Local differentiation (see FIG. 6)
In the first technique, it is considered that the flow rate update interval Δt is short and “ttti” in the above equation (2). That is, it is considered that t-ti ≒ 0. Then, the pressure fluctuation rate Pdot obtained as t-ti00, that is, the pressure fluctuation rate Pdot represented by the following equation (3) is set as a control target at the current pressure acquisition point (current flow rate update timing) .

図7に、理想圧力曲線Iと対比して、技法1を用いた場合のガスの充填圧力Pの変化曲線S1を示す。技法1では、どの瞬間のtについても実際の圧力Pを測定し、理想圧力曲線に適合させるための圧力変動率Pdotを求めるようにしている。しかし、この技法1では、上記の(2)式において「t−ti≒0」とするため、理想圧力曲線よりも圧力が高めになり、目標圧力Pfを超過する虞がある。   In contrast to the ideal pressure curve I, FIG. 7 shows a change curve S1 of the gas filling pressure P when the technique 1 is used. In the first technique, the actual pressure P is measured at any instant t, and the pressure variation rate Pdot for finding an ideal pressure curve is obtained. However, in this technique 1, in order to set it as "t-ti に お い て 0" in said (2) Formula, a pressure becomes high rather than an ideal pressure curve, and there exists a possibility that the target pressure Pf may be exceeded.

〔技法2:中間点微分(図8参照)〕
技法2では、上記の(2)式において「t−ti=Δt/2」であるとする。そして、t−ti=Δt/2として得られる圧力変動率Pdotを、すなわち下記(4)式で示される圧力変動率Pdotを、今回の圧力の取得ポイント(今回の流量更新タイミング)での制御目標とする。
Technique 2: Midpoint Derivation (See FIG. 8)
In Technique 2, it is assumed that “t−ti = Δt / 2” in the above equation (2). Then, the pressure fluctuation rate Pdot obtained as t-ti = Δt / 2, that is, the pressure fluctuation rate Pdot represented by the following equation (4), is the control target at the current pressure acquisition point (current flow update timing) I assume.

図9に、理想圧力曲線Iと対比して、技法2を用いた場合のガスの充填圧力Pの変化曲線S2を示す。断続的な更新パターンを要求するアプリケーションでは、流量を全体的に過大にならないように、流量を更新するポイントの中間点から予測する式が必要となる。但し、技法2では、オーバーシュートは確実に回避することができるが、充填圧力Pが目標圧力Pfに近づくまでに時間がかかり、目標圧力Pfに達しない虞がある。   In contrast to the ideal pressure curve I, FIG. 9 shows a change curve S2 of the gas filling pressure P when the technique 2 is used. In an application requiring an intermittent update pattern, a formula that predicts from the midpoint of the points at which the flow rate is updated is required so that the flow rate is not excessively excessive overall. However, in the technique 2, although overshoot can be reliably avoided, it takes time for the filling pressure P to approach the target pressure Pf, and there is a possibility that the target pressure Pf is not reached.

〔技法3:線形補間(図10参照)〕
技法3では、次回の圧力の取得ポイントでの充填圧力PをPi+1とし、今回の圧力の取得ポイントでの制御目標とされる圧力変動率Pdotが下記(5)式で表されるものとする。そして、上記(1)式において「t−ti=Δt」とすることにより、次回の圧力の取得ポイントで予想される充填圧力Pi+1を求め(下記(6)式)、この求めたPi+1を上記(5)式に代入することによって、下記(7)式で示される圧力変動率Pdotを得る。この(7)式で示される圧力変動率Pdotを今回の圧力の取得ポイントでの制御目標とする。
[Method 3: linear interpolation (see FIG. 10)]
In technique 3, assuming that the filling pressure P at the next acquisition point of pressure is P i + 1 , the pressure fluctuation rate Pdot to be controlled at the acquisition point of the current pressure is represented by the following equation (5) I assume. Then, by setting “t−ti = Δt” in the above equation (1), the filling pressure P i + 1 expected at the next acquisition point of pressure is determined (the following equation (6)), and this determined P By substituting i + 1 into the equation (5), a pressure fluctuation rate Pdot represented by the following equation (7) is obtained. The pressure fluctuation rate Pdot represented by the equation (7) is set as a control target at the current pressure acquisition point.

図11に、理想圧力曲線Iと対比して、技法3を用いた場合のガスの充填圧力Pの変化曲線S3を示す。この技法3では、実測した圧力から次に更新される圧力までの圧力変動率を線形補間して求めるため、理想圧力曲線によくフィットし、オーバーシュートしにくく、高速かつ正確に圧力容器7へのガスの充填を行うができる。また、最終的な充填時間を最も正確に見積もることができる。   In contrast to the ideal pressure curve I, FIG. 11 shows a change curve S3 of the gas filling pressure P when the technique 3 is used. In this technique 3, since the pressure fluctuation rate from the measured pressure to the pressure to be updated next is obtained by linear interpolation, it fits well to the ideal pressure curve, is unlikely to overshoot, and can be applied to the pressure vessel 7 quickly and accurately. It is possible to do gas filling. Also, the final fill time can be estimated most accurately.

〔充填時間の見積もり〕
aを精度(最終目標値±α%)、Pf-を許容圧力下限値とすると、Pf-は下記(8)式により示される。また、tfを初期圧力P0から許容圧力下限値Pf-に到達するまでの時間とすると、上記(1)式より許容圧力下限値Pf-は下記(9)式として表される。この(9)式から、許容圧力下限値Pf-に達するまでの時間tfは、下記(10)式として表される。この(10)式に、許容圧力下限値Pf-、初期圧力P0、目標圧力Pfを代入することにより、許容圧力下限値Pf-に達するまでの時間(充填時間)tfを見積もることができる。
[Estimate of filling time]
Assuming that a is an accuracy (final target value ± α%) and Pf is an allowable pressure lower limit value, Pf is expressed by the following equation (8). Further, the tf from the initial pressure P0 allowable pressure lower limit value Pf - When time to reach the permissible from the above equation (1) pressure lower limit value Pf - is expressed as the following equation (9). From this equation (9), the time tf until reaching the allowable pressure lower limit Pf is expressed as the following equation (10). By substituting the allowable pressure lower limit Pf , the initial pressure P 0, and the target pressure Pf into the equation (10), the time (filling time) tf until reaching the allowable pressure lower limit Pf can be estimated.

〔充填プロセスの終わり〕
本実施の形態では、ガスの充填圧力Pを流量更新間隔Δtで取得し、その取得したガスの充填圧力Pに基づいて圧力容器7へのガスの充填量(流量)を制御し、ガスの充填圧力Pを最終目標圧力Pfに整定させるようにする。
[End of filling process]
In the present embodiment, the filling pressure P of the gas is acquired at the flow rate update interval Δt, the filling amount (flow rate) of the gas into the pressure vessel 7 is controlled based on the acquired filling pressure P of the gas, and the gas filling is performed. The pressure P is allowed to settle to the final target pressure Pf.

しかし、充填圧力Pの上昇は指数関数的であるために、(理想的な条件で)充填圧力Pが目標圧力Pfに達するまで待つことは得策ではない(指数関数的曲線が実際に最終値に達することはない)。   However, it is not a good idea to wait for the filling pressure P to reach the target pressure Pf (under ideal conditions), as the increase in the filling pressure P is exponential (exponential curve is actually at its final value) Never reach).

また、アプリケーションにより、許容誤差の下限で、充填プロセスを一度止めたいことがあるかもしれない。しかし、ほとんどの場合、(許容差内の)圧力不足が発生することになる。   Also, depending on the application, it may be desirable to stop the filling process once, with a lower tolerance limit. However, in most cases an underpressure (within tolerance) will occur.

これに対しては、例えば、アプリケーションが目標圧力の下限側(または上限側)を要求している場合、目標圧力に対してある割合(代表的には97%)まで到達した時点で流量の更新を止めて(図12参照)、望む圧力まで同じ流量を保つようにするとよい。   For this, for example, if the application requires the lower limit side (or upper limit side) of the target pressure, the flow rate is updated when a certain percentage (typically 97%) of the target pressure is reached. (See FIG. 12) to maintain the same flow rate to the desired pressure.

〔圧力変動率Pdotからのガスの質量流量mdotの算出〕
質量流量算出部9−2は、圧力変動率算出部9−1によって算出された圧力変動率Pdotからガスの質量流量mdotを求める。この場合、質量流量算出部9−2は、ボイル・シャルルの法則から導かれる下記の(12)式を用いて、ガスの質量流量mdotを求める。
[Calculation of gas mass flow rate mdot from pressure fluctuation rate Pdot]
The mass flow rate calculation unit 9-2 obtains the mass flow rate mdot of gas from the pressure variation rate Pdot calculated by the pressure variation rate calculation unit 9-1. In this case, the mass flow rate calculating unit 9-2 obtains the mass flow rate mdot of the gas using the following equation (12) derived from Boyle-Charles' law.

〔ガスの質量流量mdotからのガスの体積流量Qの算出〕
体積流量算出部9−3は、質量流量算出部9−2によって求められた質量流量mdotからガスの体積流量Qを求める。この場合、体積流量算出部9−3は、下記の(13)式を用いてガスの標準状態の体積流量Qstdを求め、この標準状態の体積流量Qstdを単位換算して((14)式)、ガスの体積流量Qを求める。
[Calculation of volumetric flow rate Q of gas from mass flow rate mdot of gas]
The volumetric flow rate calculation unit 9-3 obtains the volumetric flow rate Q of the gas from the mass flow rate mdot obtained by the mass flow rate calculation unit 9-2. In this case, the volumetric flow rate calculation unit 9-3 obtains the volumetric flow rate Qstd of the standard state of the gas using the following equation (13), converts the volumetric flow rate Qstd of this standard state into a unit (equation (14)) , Determine the volumetric flow rate Q of the gas.

なお、上記の(13)式は、理想気体の状態方程式から次のようにして導かれる   The above equation (13) is derived as follows from the equation of state of the ideal gas

〔システムの特性による定数〕
技法1,2,3の方程式((3)式、(4)式、(7)式)から圧力変動率Pdotが分かる。圧力変動率Pdotは圧力差(Pf−P)、つまり誤差に対して比例関係にあると見ることができる。下記の(15)式にこれらの方程式の1つを代入すると、流量は(Pf−P)/Tに対して線形の関係を持つことができる。そのため、システムの特性による定数Kを見つけることができる。
[Constant according to system characteristics]
The pressure fluctuation rate Pdot is known from the equations (3), (4) and (7) of techniques 1, 2 and 3. The pressure fluctuation rate Pdot can be seen to be proportional to the pressure difference (Pf-P), that is, the error. Substituting one of these equations into the following equation (15), the flow rate can have a linear relationship with (Pf−P) / T. Therefore, it is possible to find a constant K depending on the characteristics of the system.

〔チューニング〕
本発明では、たった1つのパラメータ(時定数τ)を実験で求めるだけで、ガスの充填圧力Pの変化曲線を理想圧力曲線にフィットさせるようにすることが可能となる。この発明において、流量更新間隔Δtとしては、機器で許容される最小の値を推奨するが、必要に応じてその値を拡大させてもよい。
〔tuning〕
According to the present invention, it is possible to fit the change curve of the gas filling pressure P to the ideal pressure curve only by experimentally determining only one parameter (time constant τ). In this invention, although the minimum value permitted by the device is recommended as the flow rate update interval Δt, the value may be expanded as necessary.

本発明の主とする目的の1つは、望ましくないオーバーシュートおよび長い安定化期間の原因となる、充填プロセスの突然の変化を避けることである。本発明において、曲線の滑らかさは、圧力変動率Pdotの比を用いて定量化することができる。これは、τの関数で求めることができる。   One of the main objects of the present invention is to avoid sudden changes in the filling process that cause undesirable overshoot and long stabilization periods. In the present invention, the smoothness of the curve can be quantified using the ratio of pressure regulation Pdot. This can be obtained as a function of τ.

また、本発明において、例えば技法3を用いた場合、連続する圧力変動率Pdotの比を用いて(図13参照)、下記に示す(23)式より、時定数τの下限を見つけることができる。   In the present invention, for example, when technique 3 is used, the lower limit of the time constant τ can be found from the following equation (23) using the ratio of continuous pressure fluctuation rates Pdot (see FIG. 13). .

また、充填の97%が完了する前に、更新回数に下限を強制することによって、下記に示す(26)式により、時定数τの下限を見つけることができる(技法1〜3共通)。   In addition, the lower limit of the time constant τ can be found by the following equation (26) by forcing the lower limit to the number of updates before 97% of the filling is completed (common to techniques 1 to 3).

図14に、理想圧力曲線Iと対比して、技法1,2,3を用いた場合のガスの充填圧力Pの変化曲線S1,S2,S3を示す。S1,S2,S3を比較すると、技法3を用いた場合の充填圧力Pの変化曲線S3が最も理想圧力曲線Iにフィットしており、最適であることが分かる。   FIG. 14 shows, in contrast to the ideal pressure curve I, change curves S1, S2 and S3 of the gas filling pressure P when the techniques 1, 2 and 3 are used. When S1, S2 and S3 are compared, it can be seen that the change curve S3 of the filling pressure P when the technique 3 is used fits the most ideal pressure curve I and is optimum.

技法1では、常に理想圧力曲線Iよりも高くなるので、圧力変動率Pdotは常に過剰となり、目標圧力Pfに速く到達するが、整定するのに時間がかかる。但し、P0とPfとの差が大きい場合に、スタート時の所定時間だけ技法1で制御し、その後、技法3に切り替えて最後まで制御するという制御方法はあり得る。   In the technique 1, the pressure fluctuation rate Pdot is always excessive because it is always higher than the ideal pressure curve I, and the target pressure Pf is quickly reached, but it takes time to settle. However, when the difference between P0 and Pf is large, there may be a control method of controlling with the technique 1 for a predetermined time at the start and then switching to the technique 3 and controlling until the end.

技法2では、オーバーシュートは確実に避けられるが、圧力変動率Pdotの上昇が緩いので、目標圧力Pfに整定するまでに時間がかかる。但し、技法3で制御を開始し、目標圧力Pfに整定する直前で技法2に切り替えて、オーバーシュートを確実に回避するという制御方法はあり得る。   In the technique 2, overshooting is surely avoided, but since the rise of the pressure fluctuation rate Pdot is slow, it takes time to settle to the target pressure Pf. However, there may be a control method of starting control by technique 3 and switching to technique 2 immediately before settling to the target pressure Pf to reliably avoid overshoot.

〔実施の形態の拡張〕
以上、実施の形態を参照して本発明を説明したが、本発明は上記の実施の形態に限定されるものではない。本発明の構成や詳細には、本発明の技術思想の範囲内で当業者が理解し得る様々な変更をすることができる。
[Extension of the embodiment]
Although the present invention has been described above with reference to the embodiments, the present invention is not limited to the above embodiments. Various changes that can be understood by those skilled in the art can be made to the configuration and details of the present invention within the scope of the technical idea of the present invention.

1(1−1〜1−3)…流量コントローラ、2(2−1〜2−3)…流量コントローラ用遮断弁、3,4…遮断弁、5(5−1,5−2)…ゲート弁、6…真空ポンプ、7…圧力容器、8…圧力計、9…制御装置、9−1…圧力変動率算出部、9−2…質量流量算出部9−2、9−3…体積流量算出部、10…温度センサ、11…充填装置、100…ガス充填装置。   1 (1-1 to 1-3) ... flow controller, 2 (2-1 to 2-3) ... shut off valve for flow controller, 3, 4 ... shut off valve, 5 (5-1 to 5-2) ... gate Valve 6: Vacuum pump 7: Pressure vessel 8. Pressure gauge 9. Controller: 9-1: pressure fluctuation rate calculation unit 9-2: mass flow rate calculation unit 9-2, 9-3: volumetric flow rate Calculation unit, 10: temperature sensor, 11: filling device, 100: gas filling device.

Claims (8)

圧力容器にガスを充填するように構成されたガス充填装置において、
前記圧力容器内の前記ガスの充填圧力Pを計測するように構成された圧力計測部と、
前記圧力計測部が計測する前記ガスの充填圧力Pを流量更新間隔Δtで取得し、その取得したガスの充填圧力Pに基づいて前記圧力容器への前記ガスの充填量を制御し、前記ガスの充填圧力Pを最終目標圧力Pfに整定させるように構成された制御装置とを備え、
前記制御装置は、
前記ガスの充填圧力Pを取得する毎に、前記ガスの充填量を制御する際の制御目標として、前記ガスの充填圧力Pの変化曲線を下記(A)式で規定される理想圧力曲線に近似させるような前記充填圧力Pの圧力変動率Pdotを算出するように構成された圧力変動率算出部
を備えることを特徴とするガス充填装置。
In a gas filling apparatus configured to fill the pressure vessel with gas,
A pressure measurement unit configured to measure a filling pressure P of the gas in the pressure vessel;
The filling pressure P of the gas measured by the pressure measurement unit is acquired at a flow rate update interval Δt, and the filling amount of the gas into the pressure vessel is controlled based on the acquired filling pressure P of the gas, A controller configured to settle the filling pressure P to a final target pressure Pf;
The controller is
Every time the filling pressure P of the gas is obtained, the change curve of the filling pressure P of the gas is approximated to the ideal pressure curve defined by the following equation (A) as a control target when controlling the filling amount of the gas A gas filling apparatus, comprising: a pressure fluctuation rate calculation unit configured to calculate a pressure fluctuation rate Pdot of the filling pressure P that causes the pressure fluctuation to occur.
請求項1に記載されたガス充填装置において、
前記圧力変動率算出部は、
前記圧力変動率Pdotを下記(B)式に基づいて算出する
ことを特徴とするガス充填装置。
In the gas filling device according to claim 1,
The pressure fluctuation rate calculation unit
A gas charging apparatus, wherein the pressure fluctuation rate Pdot is calculated based on the following equation (B).
請求項1に記載されたガス充填装置において、
前記圧力変動率算出部は、
前記圧力変動率Pdotを下記(C)式に基づいて算出する
ことを特徴とするガス充填装置。
In the gas filling device according to claim 1,
The pressure fluctuation rate calculation unit
A gas charging apparatus, wherein the pressure fluctuation rate Pdot is calculated based on the following equation (C).
請求項1に記載されたガス充填装置において、
前記圧力変動率算出部は、
前記圧力変動率Pdotを下記(D)式に基づいて算出する
ことを特徴とするガス充填装置。
In the gas filling device according to claim 1,
The pressure fluctuation rate calculation unit
A gas filling apparatus, which calculates the pressure fluctuation rate Pdot based on the following equation (D).
請求項2に記載されたガス充填装置において、
前記圧力変動率算出部は、
前記制御装置の制御開始から所定期間のみ一時的に下記(E)式に基づいて前記圧力変動率Pdotを算出し、その後、前記(B)式に基づいて前記圧力変動率Pdotを算出する
ことを特徴とするガス充填装置。
In the gas filling apparatus according to claim 2,
The pressure fluctuation rate calculation unit
The pressure variation rate Pdot is temporarily calculated based on the following equation (E) only for a predetermined period from the start of control of the control device, and thereafter the pressure variation rate Pdot is calculated based on the equation (B). Characteristic gas filling device.
請求項2に記載されたガス充填装置において、
前記圧力変動率算出部は、
前記制御装置の制御開始から前記(B)式に基づいて前記圧力変動率Pdotを算出し、制御終了直前で下記(F)式に基づいて前記圧力変動率Pdotを算出する
ことを特徴とするガス充填装置。
In the gas filling apparatus according to claim 2,
The pressure fluctuation rate calculation unit
The pressure variation rate Pdot is calculated based on the equation (B) from the start of control of the control device, and the pressure variation rate Pdot is calculated based on the following equation (F) immediately before the end of control. Filling device.
請求項1〜6の何れか1項に記載されたガス充填装置において、
前記制御装置は、
さらに、
前記圧力変動率算出部によって算出された前記圧力変動率Pdotから前記ガスの質量流量を求めるように構成された質量流量算出部と、
前記質量流量算出部によって求められた前記質量流量から前記ガスの体積流量を求めるように構成された体積流量算出部と
を備えることを特徴とするガス充填装置。
In the gas charging apparatus according to any one of claims 1 to 6,
The controller is
further,
A mass flow rate calculation unit configured to obtain a mass flow rate of the gas from the pressure fluctuation rate Pdot calculated by the pressure fluctuation rate calculation unit;
And a volumetric flow rate calculating unit configured to determine a volumetric flow rate of the gas from the mass flow rate determined by the mass flow rate calculating unit.
圧力容器にガスを充填するガス充填方法において、
前記圧力容器内の前記ガスの充填圧力Pを計測する圧力計測ステップと、
前記圧力計測ステップによって計測される前記ガスの充填圧力Pを流量更新間隔Δtで取得し、その取得したガスの充填圧力Pに基づいて前記圧力容器への前記ガスの充填量を制御し、前記ガスの充填圧力Pを最終目標圧力Pfに整定させる制御ステップとを備え、
前記制御ステップは、
前記ガスの充填圧力Pを取得する毎に、前記ガスの充填量を制御する際の制御目標として、前記ガスの充填圧力Pの変化曲線を下記(G)式で規定される理想圧力曲線に近似させるような前記充填圧力Pの圧力変動率Pdotを算出する圧力変動率算出ステップ
を備えることを特徴とするガス充填方法。
In a gas filling method for filling a pressure vessel with gas,
A pressure measurement step of measuring a filling pressure P of the gas in the pressure vessel;
The filling pressure P of the gas measured in the pressure measurement step is acquired at a flow rate update interval Δt, and the filling amount of the gas into the pressure vessel is controlled based on the acquired filling pressure P of the gas, the gas Control step of settling the filling pressure P at the final target pressure Pf,
The control step is
Every time the filling pressure P of the gas is acquired, the change curve of the filling pressure P of the gas is approximated to the ideal pressure curve defined by the following equation (G) as a control target when controlling the filling amount of the gas A pressure fluctuation rate calculating step of calculating a pressure fluctuation rate Pdot of the filling pressure P which causes the gas filling method to have
JP2017173016A 2017-09-08 2017-09-08 Gas charging device and gas charging method Pending JP2019049283A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017173016A JP2019049283A (en) 2017-09-08 2017-09-08 Gas charging device and gas charging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017173016A JP2019049283A (en) 2017-09-08 2017-09-08 Gas charging device and gas charging method

Publications (1)

Publication Number Publication Date
JP2019049283A true JP2019049283A (en) 2019-03-28

Family

ID=65904955

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017173016A Pending JP2019049283A (en) 2017-09-08 2017-09-08 Gas charging device and gas charging method

Country Status (1)

Country Link
JP (1) JP2019049283A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112556774A (en) * 2020-11-27 2021-03-26 国网安徽省电力有限公司电力科学研究院 Gradient inflatable SF based on weighing method6Method for measuring volume of gas chamber
JP7464080B2 (en) 2022-06-10 2024-04-09 トヨタ自動車株式会社 Gas Charging System

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112556774A (en) * 2020-11-27 2021-03-26 国网安徽省电力有限公司电力科学研究院 Gradient inflatable SF based on weighing method6Method for measuring volume of gas chamber
JP7464080B2 (en) 2022-06-10 2024-04-09 トヨタ自動車株式会社 Gas Charging System

Similar Documents

Publication Publication Date Title
EP1321836A1 (en) Controller, temperature controller and heat processor using same
US7971604B2 (en) Flow controller delivery of a specified-quantity of a fluid
JP5342489B2 (en) Chemical supply system
EP1996988B1 (en) Pressure control system with optimized performance
JP6778476B2 (en) How to fill a container with a proportional valve
US9975177B2 (en) Method for a pouring control and a storage medium for storing programs for causing a computer to carry out a process for controlling pouring
JPH10318496A (en) Charging method of gas cylinder
JP7270988B2 (en) Flow control device and flow control method
RU2017135262A (en) DEVICE AND METHOD FOR FILLING A CONTAINER WITH A PRODUCT
JP2019049283A (en) Gas charging device and gas charging method
US11378073B2 (en) Estimation device and valve control device
US7270141B2 (en) Methods and systems for controlling viscosity in real time
US10316835B2 (en) Method of determining output flow rate of gas output by flow rate controller of substrate processing apparatus
JP7273596B2 (en) Flow rate calculation device, flow rate calculation system, and flow rate calculation device program
JP2011202553A (en) Feedback control device
US10248137B2 (en) Method for controlling flow rate of fluid, mass flow rate control device for executing method, and mass flow rate control system utilizing mass flow rate control device
JP2018127902A (en) Automatic pressure regulating valve and vacuum pumping system
US11326914B2 (en) Flow rate measurement apparatus and method for more accurately measuring gas flow to a substrate processing system
JPS6233156B2 (en)
JP2020009111A (en) Estimation device and pressure control valve
RU2587526C1 (en) Air pressure regulator in antechamber of wind tunnel
KR20170137640A (en) A closed-loop control device for controlling at least one control value of at least one tempering circle
JP3580353B2 (en) Pressure generator
JP2020124696A (en) Flow quantity control system, flow quantity control method and flow quantity control program
US20190094892A1 (en) Pressure control method for process chamber and pressure control device for process chamber