JP5912547B2 - Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber - Google Patents

Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber Download PDF

Info

Publication number
JP5912547B2
JP5912547B2 JP2012003052A JP2012003052A JP5912547B2 JP 5912547 B2 JP5912547 B2 JP 5912547B2 JP 2012003052 A JP2012003052 A JP 2012003052A JP 2012003052 A JP2012003052 A JP 2012003052A JP 5912547 B2 JP5912547 B2 JP 5912547B2
Authority
JP
Japan
Prior art keywords
air chamber
pressure
gas
volume
furnace
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.)
Expired - Fee Related
Application number
JP2012003052A
Other languages
Japanese (ja)
Other versions
JP2013141691A (en
Inventor
大西 亮
亮 大西
伊藤 俊一郎
俊一郎 伊藤
正木 健
健 正木
荘太 鳥居
荘太 鳥居
憲一 増田
憲一 増田
亮平 堀井
亮平 堀井
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.)
Honda Motor Co Ltd
Original Assignee
Honda Motor Co Ltd
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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2012003052A priority Critical patent/JP5912547B2/en
Publication of JP2013141691A publication Critical patent/JP2013141691A/en
Application granted granted Critical
Publication of JP5912547B2 publication Critical patent/JP5912547B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Casting Support Devices, Ladles, And Melt Control Thereby (AREA)

Description

本発明は、鋳造に用いる溶湯を保持する炉における空気室の体積を得る技術、およびそれを利用した技術に関する。   The present invention relates to a technique for obtaining a volume of an air chamber in a furnace for holding a molten metal used for casting, and a technique using the technique.

低圧鋳造装置として、炉(坩堝)内の湯面の上の空間である空気室に高圧空気を供給し、溶湯に浸漬したストークを介して溶湯を金型に鋳造する構造のものが知られている。このような低圧鋳造装置では、高圧空気の供給量を制御して金型内での溶湯の挙動を適正にする必要がある。ところで、上記のような低圧鋳造装置では、圧縮性のある空気によって溶湯という慣性流体を加圧するため、フィードバック制御では指示した空気圧力に対して実際の圧力上昇に遅れが生じる。そして、その遅れがフィードバック制御に反映されて次は指示した圧力よりも実際の圧力が高くなるオーバーシュートが生じる。このため、空気圧力が大きく変動し、金型に向かって輸送される溶湯に脈動が発生する。その結果、金型内で湯面が波立ったり溶湯が逆流したりし、湯境や未充填、あるいは空気の巻き込みなどが生じて鋳造不良が発生する確率が高くなる。   As a low-pressure casting apparatus, a structure is known in which high-pressure air is supplied to an air chamber, which is a space above a molten metal surface in a furnace (crucible), and the molten metal is cast into a mold through stalk immersed in the molten metal. Yes. In such a low-pressure casting apparatus, it is necessary to control the amount of high-pressure air supplied so that the behavior of the molten metal in the mold is appropriate. By the way, in the low-pressure casting apparatus as described above, an inertial fluid called a molten metal is pressurized by compressible air, and therefore the actual pressure rise is delayed with respect to the commanded air pressure in the feedback control. Then, the delay is reflected in the feedback control, and an overshoot occurs in which the actual pressure becomes higher than the instructed pressure. For this reason, the air pressure fluctuates greatly, and pulsation occurs in the molten metal transported toward the mold. As a result, the surface of the molten metal undulates or the molten metal flows backward in the mold, and there is a high probability that a casting defect will occur due to a hot water boundary, unfilled, or air entrainment.

この指示した値に対して遅れやオーバーシュートを生じることなく湯面の位置を制御することができ、金型内での溶湯の挙動を適正にすることができる低圧鋳造の制御技術として、特許文献1に記載の内容が提案されている。   As a low pressure casting control technique that can control the position of the molten metal surface without causing delay or overshoot with respect to the indicated value, and can make the behavior of the molten metal in the mold appropriate, Patent Literature The content of 1 is proposed.

特開2010−227974号公報JP 2010-227974 A

特許文献1に記載の技術では、ある時点から所定の時間が経過した段階におけるストーク内の湯面の高さを演算により予測し、この予測された湯面の高さが指定する高さとなるように溶湯に加える圧力を制御する。この技術では、算出式に基づいて、溶湯を加圧するガスの圧力が算出されるのであるが、この算出において、炉内の空気室(空気が占めている空間)の体積の値が必要となる。   In the technique described in Patent Document 1, the height of the molten metal surface in the stalk is predicted by calculation at a stage where a predetermined time has elapsed from a certain point in time, and the predicted molten metal surface height becomes the specified height. To control the pressure applied to the melt. In this technique, the pressure of the gas for pressurizing the molten metal is calculated based on the calculation formula. In this calculation, the value of the volume of the air chamber (the space occupied by air) in the furnace is required. .

ところで、溶湯を保持する炉の内壁は、耐火性の材料で構成されているのであるが、以下の点が明らかになっている。
(1)炉の内壁は、耐火材料のキャスタブルを手作業で塗ったものであり、設計値とはその内部寸法が異なり、またその値も炉毎に異なっている。
(2)炉の内壁を構成する上記の耐火材料は、多孔質構造を有し、そこにガスが浸透する。
これらの理由により、上記の演算で用いる空気室の体積には、誤差が含まれており、それが湯面制御の不確定要素となる。また、炉は、密閉構造とされているが、内部を正圧とした場合、僅かなガスのリークがあり、しかもその値に個体バラツキや経時変化があり、これも湯面制御の不確定要素となる。
By the way, although the inner wall of the furnace holding the molten metal is made of a refractory material, the following points have been clarified.
(1) The inner wall of the furnace is made by hand-casting a refractory material, and its internal dimensions are different from the design values, and the values are also different for each furnace.
(2) The above-mentioned refractory material constituting the inner wall of the furnace has a porous structure, and gas permeates there.
For these reasons, an error is included in the volume of the air chamber used in the above calculation, which becomes an uncertain element of the hot water surface control. In addition, the furnace has a sealed structure, but if the inside is set to a positive pressure, there is a slight gas leak, and there are individual variations and changes over time, and this is also an uncertain factor in hot water level control. It becomes.

このような背景において、本発明は、炉の空気室の実効体積およびリーク量を正確に把握し、より正確な湯面のモデル予測を可能とする技術を提供することを目的とする。   In such a background, an object of the present invention is to provide a technique capable of accurately grasping an effective volume and a leak amount of an air chamber of a furnace and enabling a more accurate model prediction of a molten metal surface.

請求項1に記載の発明は、炉内の空気室を大気圧よりも高い所定の圧力にした際における前記空気室からの気体のリーク量を計測するリーク量計測工程と、所定の流量の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から前記所定の圧力になるまでに掛かる経過時間を計測する経過時間計測工程と、前記リーク量、前記計測された前記経過時間、前記経過時間計測工程において前記空気室に供給された前記気体の流量、前記所定圧力および前記空気室の温度に基づき、気体の状態方程式を利用して前記空気室の体積を算出する空気室体積算出工程とを有することを特徴とする炉の空気室の体積を算出する方法である。   According to the first aspect of the present invention, there is provided a leak amount measuring step for measuring a leak amount of gas from the air chamber when the air chamber in the furnace is set to a predetermined pressure higher than the atmospheric pressure, and a gas having a predetermined flow rate. When the air chamber is sent to the air chamber at atmospheric pressure, the elapsed time measuring step of measuring the elapsed time until the air chamber is changed from the atmospheric pressure to the predetermined pressure, the leak amount, the measured Based on the elapsed time, the flow rate of the gas supplied to the air chamber in the elapsed time measuring step, the predetermined pressure, and the temperature of the air chamber, the volume of the air chamber is calculated using a gas equation of state. An air chamber volume calculating step, wherein the volume of the air chamber of the furnace is calculated.

請求項1に記載の発明によれば、空気室からの気体のリーク量に着目し、このリーク量から気体の状態方程式を利用して空気室の体積を算出する。この方法によれば、炉内壁への気体の浸透分や、設計図からは把握されない炉内のスペース等に存在する気体の量を含む実際に空気室に存在している気体の量(実効値)を高い精度で取得することができる。ここで、空気室は、炉内の溶湯と接し、溶湯の上部における気体で占められている空間のことをいう。気体は、溶湯の湯面に圧力を加えるために利用されるもので、一般に空気が用いられるが、窒素ガス等の不活性気体であってもよい。また、溶湯としては、アルミニウム、アルミニウム合金、鉄等の溶湯が挙げられるが、鋳造で用いる金属であれば特に限定されない。   According to the first aspect of the present invention, attention is paid to the amount of gas leakage from the air chamber, and the volume of the air chamber is calculated from the amount of leakage using the gas equation of state. According to this method, the amount of gas that actually exists in the air chamber (effective value), including the amount of gas permeating into the furnace inner wall and the amount of gas present in the space in the furnace that cannot be grasped from the design drawing. ) Can be obtained with high accuracy. Here, the air chamber refers to a space that is in contact with the molten metal in the furnace and is occupied by gas in the upper part of the molten metal. The gas is used to apply pressure to the surface of the molten metal, and generally air is used, but an inert gas such as nitrogen gas may be used. Moreover, as a molten metal, molten metals, such as aluminum, an aluminum alloy, and iron, are mentioned, However, If it is the metal used by casting, it will not specifically limit.

請求項2に記載の発明は、請求項1に記載の発明において、前記空気室体積算出工程では、下記「数1」に基づいて前記空気室の体積Vの算出が行われることを特徴とする。   The invention according to claim 2 is characterized in that, in the invention according to claim 1, in the air chamber volume calculating step, the volume V of the air chamber is calculated based on the following “Equation 1”. .

Figure 0005912547
P:空気室の圧力、GIn:空気室への気体の流入量、GLeak:空気室からの気体のリーク量
V:空気室の体積、T:空気室の温度、R:気体定数
Figure 0005912547
P: pressure of the air chamber, G In : amount of gas flowing into the air chamber, G Leak : amount of gas leakage from the air chamber V: volume of the air chamber, T: temperature of the air chamber, R: gas constant

請求項3に記載の発明は、請求項1または2に記載の発明において、前記炉を用いた鋳造を行う前と後に前記空気室体積算出工程を行うことを特徴とする。請求項3に記載の発明によれば、鋳造を行うことで空気室の体積に変化が生じても、すなわち当該体積に経時変化が生じても、その変化を取り込んだ上で正確な空気室の体積を得ることができる。時間軸上で隣接する第1の空気室体積算出工程と第2の空気室体積算出工程との間における鋳造の回数は、通常複数回であるが、その数は限定されない。   The invention according to claim 3 is the invention according to claim 1 or 2, wherein the air chamber volume calculation step is performed before and after casting using the furnace. According to the third aspect of the present invention, even if a change occurs in the volume of the air chamber due to casting, that is, even if a change with time occurs in the volume, the change of the air chamber is taken into account after taking in the change. Volume can be obtained. The number of castings between the first air chamber volume calculating step and the second air chamber volume calculating step adjacent on the time axis is usually a plurality of times, but the number is not limited.

請求項4に記載の発明は、請求項1〜3のいずれか一項に記載の発明によって算出した前記炉の前記空気室の体積に基づいて、前記空気室に流入する気体による圧力増加分と、前記空気室に流入した前記気体により前記空気室において溶湯が押し下げられる結果、前記空気室の容積が増加することによる圧力減少分と、前記空気室からのリークによって生じる圧力減少分とを算出する圧力変化分算出工程と、前記圧力変化分算出工程における演算の結果に基づき、前記炉からストークを介して送り出される溶湯の湯面の高さ位置を算出する湯面位置算出工程とを有することを特徴とする鋳造方法である。   According to a fourth aspect of the present invention, based on the volume of the air chamber of the furnace calculated by the invention according to any one of the first to third aspects, the pressure increase due to the gas flowing into the air chamber and As a result of the molten metal being pushed down in the air chamber by the gas flowing into the air chamber, a pressure decrease due to an increase in the volume of the air chamber and a pressure decrease caused by a leak from the air chamber are calculated. A pressure change calculating step, and a molten metal surface position calculating step of calculating the height position of the molten metal sent from the furnace via the stalk based on the calculation result in the pressure change calculating step. This is a characteristic casting method.

請求項4に記載の発明によれば、正確な空気室の体積を用い、更に空気室からの気体のリーク分を勘案した演算が行われるので、ストーク内の湯面の位置の算出精度をより高くできる。   According to the invention described in claim 4, since the calculation is performed using the accurate volume of the air chamber and further taking into account the gas leak from the air chamber, the calculation accuracy of the position of the hot water surface in the stalk is further improved. Can be high.

請求項5に記載の発明は、請求項4に記載の発明において、前記圧力変化分算出工程が下記「数2」に基づいて行われ、前記湯面位置算出工程が下記「数3」に基づいて行われることを特徴とする。   The invention according to claim 5 is the invention according to claim 4, wherein the pressure change calculation step is performed based on the following “Equation 2”, and the molten metal surface position calculation step is based on the following “Equation 3”. It is characterized by being performed.

Figure 0005912547
P:空気室の予測ガス圧力、P:空気室の実測ガス圧力、V:空気室の体積、R:気体定数
T:空気室の温度、G:空気室に流入するガス流量、GLeak:空気室からのガスリーク量
Figure 0005912547
P 1 : Predicted gas pressure in the air chamber, P: Actual gas pressure in the air chamber, V: Volume of the air chamber, R: Gas constant T: Temperature of the air chamber, G: Gas flow rate flowing into the air chamber, G Leak : Gas leakage from the air chamber

Figure 0005912547
P:空気室の予測ガス圧力、S:ストーク内部の水平断面積、Pi:溶湯の湯面にかかる背圧
mA:ストーク内部の溶湯の質量、g:重力加速度、h:送り出された溶湯の湯面の高さ
μ(h):溶湯の粘性係数
Figure 0005912547
P 1 : Predicted gas pressure in the air chamber, S: Horizontal cross-sectional area inside the stalk, P i : Back pressure applied to the molten metal surface
m A : Mass of molten metal inside stalk, g: Gravitational acceleration, h 1 : Height of molten metal surface sent out μ (h): Viscosity coefficient of molten metal

請求項6に記載の発明は、炉内の空気室を大気圧よりも高い所定の圧力にした際における前記空気室からの気体のリーク量を計測するリーク量計測手段と、所定の流量の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から前記所定の圧力になるまでに掛かる経過時間を計測する経過時間計測手段と、前記リーク量、前記計測された前記経過時間、前記気体の前記所定の流量、前記所定圧力および前記空気室の温度に基づき、気体の状態方程式を利用して前記空気室の体積を算出する空気室体積算出手段とを有することを特徴とする炉の空気室の体積を算出する装置である。   According to a sixth aspect of the present invention, there is provided a leak amount measuring means for measuring a leak amount of gas from the air chamber when the air chamber in the furnace is set to a predetermined pressure higher than atmospheric pressure, and a gas having a predetermined flow rate. , An elapsed time measuring means for measuring an elapsed time until the air chamber is changed from atmospheric pressure to the predetermined pressure, and the amount of leakage is measured. Air chamber volume calculating means for calculating the volume of the air chamber using a gas equation of state based on the elapsed time, the predetermined flow rate of the gas, the predetermined pressure, and the temperature of the air chamber. This is a device for calculating the volume of the air chamber of the furnace.

請求項7に記載の発明は、コンピュータに読み取らせて実行させるプログラムであって、コンピュータに読み取らせて実行させるプログラムであって、CPUに、炉内の空気圧を大気圧よりも高い所定の圧力保つようにした際に、前記空気室への空気の流入量を計測することで前記空気室からの気体のリーク量 Leak を計測するリーク量計測ステップと、前記リーク量の計測後に前記空気室を大気圧に戻し、前記リーク量G Leak より多い流量G In で前記空気室に前記圧力Pになるまで空気を流入させるように比例弁を制御する制御ステップと、所定の流量G In の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から前記所定の圧力Pになるまでに掛かる経過時間の範囲で上記「数1」の計算を行い、前記空気室の体積Vを算出する空気室体積算出ステップとを実行させることを特徴とする炉の空気室の体積を算出するためのプログラムである。 The invention according to claim 7 is a program that is read and executed by a computer, and is a program that is read and executed by a computer. The CPU causes the air pressure in the furnace to be a predetermined pressure P higher than atmospheric pressure. It said air upon way kept, and the leakage amount measuring step of measuring the leakage amount G leak of gas from said air chamber by measuring the inflow of air into the air chamber, after the measurement of the leakage quantity in A control step for controlling the proportional valve so that the chamber is returned to atmospheric pressure and air is introduced into the air chamber at a flow rate G In greater than the leak amount G Leak until the pressure P is reached, and a gas having a predetermined flow rate G In When the air chamber is sent to the air chamber at atmospheric pressure, the calculation of the above-mentioned “Equation 1” is performed within the range of the elapsed time from the atmospheric pressure to the predetermined pressure P. Volume of A program for calculating the volume of the air chamber of the furnace, characterized in that to execute the air chamber volume calculation step of calculating a.

本発明によれば、炉の空気室の実効体積およびリーク量が正確に把握され、より正確な湯面のモデル予測を可能とする技術が提供される。   ADVANTAGE OF THE INVENTION According to this invention, the effective volume and leak amount of the air chamber of a furnace are grasped | ascertained correctly, and the technique which enables the more accurate model prediction of the hot_water | molten_metal surface is provided.

発明を利用した低圧鋳造装置の一例を示す概念図である。It is a conceptual diagram which shows an example of the low pressure casting apparatus using invention. 空気室の体積算出装置のブロック図である。It is a block diagram of the volume calculation apparatus of an air chamber. 空気室の体積を求める処理の一例を示すフローチャートである。It is a flowchart which shows an example of the process which calculates | requires the volume of an air chamber. 実施形態のモデル予測制御を説明するためのタイムチャートである。It is a time chart for demonstrating the model prediction control of embodiment. 実施形態のモデル予測制御を説明するためのタイムチャートである。It is a time chart for demonstrating the model prediction control of embodiment. 図5に示す状態から1ステップ進んだ状態を示すタイムチャートである。It is a time chart which shows the state which advanced 1 step from the state shown in FIG. 図6に示す状態から1ステップ進んだ状態を示すタイムチャートである。It is a time chart which shows the state advanced 1 step from the state shown in FIG. 図7に示す状態から1ステップ進んだ状態を示すタイムチャートである。It is a time chart which shows the state which advanced 1 step from the state shown in FIG. 図8に示す状態から1ステップ進んだ状態を示すタイムチャートである。It is a time chart which shows the state advanced 1 step from the state shown in FIG.

以下、本発明の一実施形態を説明する。   Hereinafter, an embodiment of the present invention will be described.

(システムの構成)
図1には、本発明を利用した低圧鋳造装置10が示されている。低圧鋳造装置10は、炉11とそれに付随する周辺装置を備えている。まず、低圧鋳造装置10は、炉(坩堝)11を備えている。炉11の上には、下型12とこの下型12に対して上下方向に接近離間可能な上型13とが配置され、下型12および上型13によって鋳型のキャビティ14が形成されている。
(System configuration)
FIG. 1 shows a low pressure casting apparatus 10 utilizing the present invention. The low pressure casting apparatus 10 includes a furnace 11 and peripheral devices associated therewith. First, the low-pressure casting apparatus 10 includes a furnace (crucible) 11. On the furnace 11, a lower mold 12 and an upper mold 13 that can be approached and separated in the vertical direction with respect to the lower mold 12 are arranged. A mold cavity 14 is formed by the lower mold 12 and the upper mold 13. .

炉11の頂部には、軸線を上下方向に向けた中空パイプ構造のストーク15が貫通して配置され、ストーク15の上端部は下型12を貫通してキャビティ14に臨んでいる。加熱炉11の側壁の上部には、図示しない空気圧縮装置に接続された配管16が比例弁(流量調整弁)17を介して接続され、内部に保持された溶湯Mの上部の空間であり、空気室の一例である上部空間Aに空気を流入させるようになっている。また、炉11には、溶湯Mの上部空間Aの圧力を検出する圧力センサ(圧力検出手段)18が設けられている。この低圧鋳造装置10では、炉11における溶湯Mの上部空間Aの圧力(炉内圧力)を高めることにより、炉11内の溶湯Mがストーク15に送り出され、ストーク15内を溶湯が上昇することで、溶湯がキャビティ14に充填される。   At the top of the furnace 11, a stalk 15 having a hollow pipe structure with the axis line directed in the vertical direction is disposed so as to pass through, and the upper end portion of the stalk 15 passes through the lower mold 12 and faces the cavity 14. A pipe 16 connected to an air compressor (not shown) is connected to the upper part of the side wall of the heating furnace 11 via a proportional valve (flow rate adjusting valve) 17 and is an upper space of the molten metal M held inside, Air is allowed to flow into the upper space A, which is an example of an air chamber. The furnace 11 is provided with a pressure sensor (pressure detection means) 18 for detecting the pressure in the upper space A of the molten metal M. In this low pressure casting apparatus 10, by increasing the pressure (in-furnace pressure) of the upper space A of the molten metal M in the furnace 11, the molten metal M in the furnace 11 is sent to the stalk 15, and the molten metal rises in the stalk 15. Then, the molten metal is filled in the cavity 14.

低圧鋳造装置10は、制御部20およびモデル予測コントローラ21を備えている。制御部20は、低圧鋳造装置10の動作を制御するとともに、モデル予測コントローラ21に、溶湯Mの湯面位置を指示する湯面位置指示信号を出力する。モデル予測コントローラ21は、圧力センサ18から入力される圧力情報に基づき、現在から数ステップ(例えば3ステップ)先までの湯面位置と炉内圧力を予測するともに、予測された湯面位置と指示された湯面位置とを比較し、予測された湯面位置が指示された湯面位置に近似するように比例弁17を制御する。なお、比例弁17は、通過する気体の流量を検出する流量検出センサ19を備えており、その検出結果はモデル予測コントローラ21に入力される。   The low pressure casting apparatus 10 includes a control unit 20 and a model prediction controller 21. The control unit 20 controls the operation of the low-pressure casting apparatus 10 and outputs a molten steel surface position instruction signal that instructs the model prediction controller 21 of the molten metal surface position of the molten metal M. Based on the pressure information input from the pressure sensor 18, the model prediction controller 21 predicts the molten metal surface position and the furnace pressure from the present to several steps (for example, three steps) ahead, and the predicted molten metal surface position and instruction. The proportional valve 17 is controlled so that the predicted hot water surface position is approximated to the instructed hot water surface position. The proportional valve 17 includes a flow rate detection sensor 19 that detects the flow rate of the passing gas, and the detection result is input to the model prediction controller 21.

図1に示すシステムは、空気室の体積算出装置30を備えている。空気室の体積算出装置30は、コンピュータとしての機能を有し、CPU、半導体メモリやハードディスク装置等のデータ記憶装置、各種インターフェースを備えている。図2には、空気室の体積算出装置30のブロック図が示されている。空気室の体積算出装置30は、ソフトウェアー的(勿論、一部に専用のハードウェアを含んでもよい)に構成された機能部として、気体リーク量計測部31、経過時間計測部32、気体流入量計測部33、空気室圧力計測部34、空気室温度計測部35、空気室体積算出部36および制御信号出力部37を有している。   The system shown in FIG. 1 includes an air chamber volume calculation device 30. The volume calculation device 30 of the air chamber has a function as a computer, and includes a CPU, a data storage device such as a semiconductor memory and a hard disk device, and various interfaces. FIG. 2 shows a block diagram of the air chamber volume calculation device 30. The air chamber volume calculation device 30 includes a gas leak amount measurement unit 31, an elapsed time measurement unit 32, a gas inflow as functional units configured in software (of course, some may include dedicated hardware). It has a quantity measuring unit 33, an air chamber pressure measuring unit 34, an air chamber temperature measuring unit 35, an air chamber volume calculating unit 36, and a control signal output unit 37.

気体リーク量計測部31は、流量検出センサ19の出力に基づき、上部空間Aからの空気のリーク量を算出する。経過時間計測部32は、配管16を用いて上部空間Aに空気を流入させる際における空気の流入時間を計測する。気体流入量計測部33は、流量検出センサ19の出力に基づき、配管16から上部空間Aに流入する空気の流入量を計測する。空気室圧力計測部34は、空気室である上部空間Aの圧力を計測する。この空気室Aの圧力の計測は、圧力センサ18の出力またはストーク15内の湯面の高さ位置に基づいて行われる。空気室温度計測部35は、上部空間Aの温度を計測する。空気室体積算出部36は、後述する方法により、設計図面から得られない上部空間Aの実効体積の算出を行う。制御信号出力部37は、後述する図3の処理を実行する際に必要な、比例弁19その他への制御信号を生成し、それを出力する。   The gas leak amount measuring unit 31 calculates the amount of air leak from the upper space A based on the output of the flow rate detection sensor 19. The elapsed time measurement unit 32 measures the inflow time of air when the air is caused to flow into the upper space A using the pipe 16. The gas inflow measuring unit 33 measures the inflow of air flowing from the pipe 16 into the upper space A based on the output of the flow rate detection sensor 19. The air chamber pressure measuring unit 34 measures the pressure in the upper space A that is an air chamber. The pressure in the air chamber A is measured based on the output of the pressure sensor 18 or the height position of the hot water surface in the stalk 15. The air chamber temperature measurement unit 35 measures the temperature of the upper space A. The air chamber volume calculation unit 36 calculates the effective volume of the upper space A that cannot be obtained from the design drawing by a method described later. The control signal output unit 37 generates a control signal for the proportional valve 19 and others necessary for executing the process of FIG. 3 described later, and outputs it.

(上部空間Aの体積の算出)
前述したように、炉11の空気室である上部空間Aの体積は、炉の設計情報からは正確な値は判らず、またその値は炉11の使用時間の経過に従って変動する。また、上部空間Aは、シール構造とされ密閉されているが、高温の状態で使用されることもあり、上部空間Aからの空気のリークが存在する。ここで、上部空間Aの圧力をP、上部空間Aへの空気の流入量をGIn、上部空間Aからの空気のリーク量をGLeak、上部空間Aの体積をV、上部空間Aの温度をT(上部空間Aに含まれる空気の温度)、気体定数をRとすると、気体の状態方程式より、下記の数4が得られる。なお、GInおよびGLeakは、移動した空気の質量(質量流量)である。
(Calculation of volume of upper space A)
As described above, the volume of the upper space A, which is the air chamber of the furnace 11, cannot be accurately determined from the design information of the furnace, and the value fluctuates with the passage of time of use of the furnace 11. The upper space A is sealed and sealed, but may be used at a high temperature, and air leaks from the upper space A. Here, P is the pressure of the upper space A, G In is the amount of air flowing into the upper space A, G Leak is the amount of air leakage from the upper space A, V is the volume of the upper space A, and the temperature of the upper space A is Is T (temperature of air contained in the upper space A), and R is a gas constant, the following equation 4 is obtained from the gas equation of state. G In and G Leak are masses (mass flow rate) of the moved air.

Figure 0005912547
P:上部空間Aの圧力、GIn:上部空間Aへの気体の流入量、GLeak:上部空間Aからの気体のリーク量、V:上部空間Aの体積、T:上部空間Aの温度、R:気体定数
Figure 0005912547
P: pressure of the upper space A, G In : amount of gas flowing into the upper space A, G Leak : amount of gas leakage from the upper space A, V: volume of the upper space A, T: temperature of the upper space A, R: Gas constant

ここで、上部空間Aの圧力は、ストーク15内における湯面の位置から算出される(圧力センサ18で計測してもよい)。GInは、気体流入量計測部33において計測される。ここでは、後述する方法でGLeakを求め、「数4」からVを算出する。 Here, the pressure in the upper space A is calculated from the position of the molten metal surface in the stalk 15 (may be measured by the pressure sensor 18). G In is measured by the gas inflow measuring unit 33. Here, G Leak is obtained by a method described later, and V is calculated from “ Equation 4”.

図3には、Vを求める手順の一例が示されている。図3に示す処理を実行するプログラムは、図1における空気室の体積算出装置0の記憶部に記憶されており、適当な記憶領域に読み出されて、空気室の体積算出装置0おいて実行される。なお、このプログラムを適当な記憶媒体に記憶させ、そこから提供する形態であってもよい。 FIG. 3 shows an example of a procedure for obtaining V. The program for executing the processing shown in FIG. 3 is stored in the storage unit of the air chamber volume calculation device 30 in FIG. 1, and is read out to an appropriate storage area so that the air chamber volume calculation device 30 And executed. The program may be stored in a suitable storage medium and provided from there.

図3の処理が開始されると(ステップS301)、まず上部空間Aからの気体(この場合は空気)のリーク量(GLeak)の計測が行われる。以下、GLeakを得る手順を説明する。まず、炉11の内部に所定の溶湯を入れた状態で、上部空間Aが存在する状態とする。ここで、炉11に入れる溶湯の量は、実際に鋳造を行う状態に近い値とする。 When the processing of FIG. 3 is started (step S301), first, a leak amount (G Leak ) of gas (in this case, air) from the upper space A is measured. The procedure for obtaining G Leak will be described below. First, it is assumed that the upper space A exists in a state where a predetermined molten metal is put in the furnace 11. Here, the amount of the molten metal put into the furnace 11 is set to a value close to a state in which casting is actually performed.

次に、空気室の体積算出部20から比例弁7に制御信号を送り、図示しない空気圧縮装置に接続された配管16から、上部空間Aに空気を流入させる。この際、上部空間Aの圧力を特定の正圧P(例えば、大気圧+11kP)の値に保つように、比例弁17を操作する。ここでは、予め求めておいたストーク15内における湯面の位置と上部空間Aの圧力との関係に基づき、ストーク15内における湯面の位置を一定に保つように比例弁17の操作を行い、上部空間Aの圧力を特定の正圧の状態とする(ステップS302)。   Next, a control signal is sent from the volume calculation unit 20 of the air chamber to the proportional valve 7 so that air flows into the upper space A from the pipe 16 connected to an air compressor (not shown). At this time, the proportional valve 17 is operated so as to keep the pressure of the upper space A at a specific positive pressure P (for example, atmospheric pressure +11 kP). Here, the proportional valve 17 is operated so as to keep the position of the hot water surface in the stalk 15 constant based on the relationship between the position of the hot water surface in the stalk 15 and the pressure in the upper space A, which has been obtained in advance. The pressure in the upper space A is set to a specific positive pressure state (step S302).

そして、この際における上部空間Aへの空気の流入量G1を気体流入量計測部33において計測する(ステップS303)。仮に、上部空間Aにリークが無ければ、圧力一定の条件で上部空間Aへの気体の流入はない。したがって、圧力一定の状態で、上部空間Aに流入する空気の流量G1は、上部空間Aで発生しているリーク量GLeakとなる。こうして、GLeakの値が得られる(ステップS304)。 Then, the inflow amount G1 of air into the upper space A at this time is measured by the gas inflow amount measuring unit 33 (step S303). If there is no leak in the upper space A, no gas flows into the upper space A under a constant pressure condition. Accordingly, the flow rate G1 of the air flowing into the upper space A in a state where the pressure is constant becomes the leak amount G Leak generated in the upper space A. In this way, the value of G Leak is obtained (step S304).

GLeakを得たら、一旦、上部空間Aの圧力を大気圧に戻す。そして、比例弁17を制御し、図示しない空気圧縮装置に接続された配管16から、上部空間Aに空気を一定流量GInで流入させる(ステップS305)。GInの値は、気体流入量計測部33で計測し、またGIn>GLeakの条件となるようにGInの流量を選択する。そして、上部空間Aの圧力が先のリーク量GLeakを求めた場合に利用した特定の正圧P(例えば、空気室大気圧+11kP)になるまでの経過時間Δtを、経過時間計測部32において計測する(ステップS306)。ここで、上部空間Aの圧力Pは、ストーク15内の湯面の高さから求められる。この処理は、空気室圧力計測部34において行われる。この際の圧力の値Pが「数4」におけるPの値となる。そして、上記のΔtの期間において、「数4」におけるGInおよびGLeakの積分をΔtの範囲で実行し、その積分値を算出する。また、図示しない温度センサにより、上部空間A内の空気の温度Tを測定する(ステップS307)。この処理は、空気室温度計測部35において行われる。 Once G Leak is obtained, the pressure in the upper space A is once returned to atmospheric pressure. Then, the proportional valve 17 is controlled to allow air to flow into the upper space A at a constant flow rate G In from the pipe 16 connected to an air compressor (not shown) (step S305). The value of G In is measured by the gas inflow measuring unit 33, and the flow rate of G In is selected so that the condition of G In > G Leak is satisfied . Then, an elapsed time Δt until the specific positive pressure P (for example, air chamber atmospheric pressure +11 kP) used when the pressure in the upper space A is obtained in the previous leak amount G Leak is calculated in the elapsed time measuring unit 32. Measurement is performed (step S306). Here, the pressure P in the upper space A is obtained from the height of the hot water surface in the stalk 15. This process is performed in the air chamber pressure measurement unit 34. The pressure value P at this time is the value of P in “Equation 4”. Then, during the period of Δt, the integration of G In and G Leak in “ Equation 4” is executed in the range of Δt, and the integration value is calculated. Further, the temperature T of the air in the upper space A is measured by a temperature sensor (not shown) (step S307). This process is performed in the air chamber temperature measurement unit 35.

こうして、「数4」におけるV以外の値が得られる。そして、「数4」に基づくVの算出を空気室体積算出部36において行い(ステップS308)、処理を終了する(ステップS309)。ステップS308において得られるVの値は、設計図面から得られない内部寸法の誤差、および炉11の内壁に浸透する空気の量を含んだ実効値となる。   Thus, values other than V in “Equation 4” are obtained. Then, the calculation of V based on “Equation 4” is performed in the air chamber volume calculation unit 36 (step S308), and the process ends (step S309). The value of V obtained in step S308 is an effective value including an error in internal dimensions that cannot be obtained from the design drawing and the amount of air that permeates the inner wall of the furnace 11.

(高精度加圧制御)
以下、図4の手順により求めた上部空間Aの実効体積Vを用いて、指示した湯面の位置に対して遅れやオーバーシュートを生じることなく、金型内での溶湯の挙動を適正に制御する方法の一例を説明する。
(High precision pressurization control)
In the following, the effective volume V of the upper space A obtained by the procedure of FIG. 4 is used to appropriately control the behavior of the molten metal in the mold without causing a delay or overshoot with respect to the position of the indicated molten metal surface. An example of the method of performing will be described.

この例において、モデル予測コントローラ21は、圧力センサ18から入力される圧力情報に基づき、現在から3ステップ先までの湯面位置と炉内圧力を予測するともに、予測された湯面位置と指示された湯面位置とを比較し、予測された湯面位置が指示された湯面位置に近似するように比例弁17を制御する。   In this example, the model prediction controller 21 predicts the molten metal surface position and the furnace pressure from the present three steps ahead based on the pressure information input from the pressure sensor 18, and is also instructed as the predicted molten metal surface position. The proportional valve 17 is controlled so that the predicted hot water surface position approximates the instructed hot water surface position.

まず、モデル予測における演算の概要を説明する。ここでは、下記「数4」および「数5」により現在から所定時間経過した時点におけるストーク15内における溶湯の湯面の高さを予測する   First, an outline of calculation in model prediction will be described. Here, the height of the molten metal surface in the stalk 15 at the time when a predetermined time has elapsed from the present time is predicted by the following “Equation 4” and “Equation 5”.

Figure 0005912547
P:上部空間Aの予測ガス圧力、S:ストーク内部の水平断面積、Pi:溶湯の湯面にかかる背圧
A:ストーク内部の溶湯の質量、g:重力加速度、P:予測ガス圧力
:送り出された溶湯の湯面の高さ、μ(h):溶湯の粘性係数
Figure 0005912547
P 1 : Predicted gas pressure in the upper space A, S: Horizontal cross-sectional area inside the stalk, P i : Back pressure applied to the molten metal surface m A : Mass of the molten metal inside the stalk, g: Gravitational acceleration, P 1 : Predicted Gas pressure h 1 : Height of molten metal surface sent out, μ (h): Viscosity coefficient of molten metal

Figure 0005912547
P:上部空間Aの予測ガス圧力、P:上部空間Aの実測ガス圧力、V:上部空間Aの体積
R:気体定数、T:上部空間A内の温度、G:上部空間Aに流入するガス流量
GLeak:上部空間Aからのガスリーク量
Figure 0005912547
P 1 : Predicted gas pressure in the upper space A, P: Actual gas pressure in the upper space A, V: Volume of the upper space A R: Gas constant, T: Temperature in the upper space A, G: Inflow into the upper space A Gas flow rate
G Leak : Gas leakage from the upper space A

ここで、Gは、流量調整弁の開度に比例するものであるが、電磁式の流量調整弁の種類によっては、操作電圧が下のしきい値を超えるまで流路が開かず、流路が開いてからは操作電圧に対して流量が線形に変化するものがある。また、このような流量調整弁では、流量調整弁を全開にした状態から絞る場合にも、操作電圧が上のしきい値を下回るまで流路が絞られず、流路が絞られ始めてからは操作電圧に対して流量が線形に変化する。したがって、操作電圧から流量を求める場合には、上記のような作動ヒステリシスを加味する必要がある。   Here, G is proportional to the opening degree of the flow rate adjusting valve. However, depending on the type of the electromagnetic flow rate adjusting valve, the channel is not opened until the operating voltage exceeds the lower threshold value. After opening, there is a thing whose flow rate changes linearly with respect to the operating voltage. Also, with such a flow rate adjustment valve, even when the flow rate adjustment valve is throttled from the fully open state, the flow path is not throttled until the operating voltage falls below the upper threshold, and the flow control valve is operated after the flow path starts to be throttled. The flow rate varies linearly with voltage. Therefore, when the flow rate is obtained from the operation voltage, it is necessary to consider the operation hysteresis as described above.

ここで、「数6」の右辺の第1項は、炉11に流入するガスによる圧力増加分であり、第2項は、炉11に流入したガスにより上部空間Aにおいて溶湯が押し下げられる結果、上部空間Aの容積が増加することによる圧力減少分である。また、第3項は、上部空間Aからのリークによって生じる圧力減少分である。   Here, the first term on the right side of “Equation 6” is a pressure increase due to the gas flowing into the furnace 11, and the second term is a result of the molten metal being pushed down in the upper space A by the gas flowing into the furnace 11, This is a pressure decrease due to an increase in the volume of the upper space A. The third term is a pressure decrease caused by a leak from the upper space A.

現在の時刻t1からt秒後のガス圧力Pは、「数6」をt1から(t1+t)まで積分することで求めることができる。この場合、ガス温度Tは一定とすることができ、ガス流量Gは、時刻(t1+t)のものとして予め設定されている。また、時刻(t1+t)のときのガス体積Vは、時刻t1のときのガス体積(既に求められている)にAh1を加えることで求めることができる。   The gas pressure P after t seconds from the current time t1 can be obtained by integrating “Equation 6” from t1 to (t1 + t). In this case, the gas temperature T can be kept constant, and the gas flow rate G is preset as that at time (t1 + t). The gas volume V at time (t1 + t) can be obtained by adding Ah1 to the gas volume at time t1 (which has already been obtained).

「数5」において背圧Piは大気圧とすることができる。また、ストーク内溶湯の質量mAは、時刻t1のときのものであり、過去の計算で既に求められている。さらに、粘性係数μ(h)は溶湯の材質および温度から既知である。そして、それらの値と「数6」で求めたガス圧力Pを「数5」に代入して積分を2回行うと、送り出された溶湯の時刻(t1+t)における湯面の高さhが求められる。この場合において、時刻(t1+t)における算出された湯面の高さhが指定した高さよりも高い場合には、流量調整弁の操作電圧を下げ、上部空間Aの圧力を下げることでストーク15内における湯面の位置が下がるように制御を行う。一方、時刻(t1+t)における湯面の高さhが指定した高さよりも低い場合には、流量調整弁の操作電圧を上げ、上部空間Aの圧力を挙げることでストーク15内における湯面の位置が上がるように制御を行う。このような制御を行うことで、時刻(t1+t)における湯面の実際の高さhを指定した高さに近付けることができる。なお、操作電圧の下げ幅は、湯面の高さhと指定した高さとの差に比例させることができる。 Back pressure P i in the "Number 5" can be atmospheric pressure. Further, the mass mA of the molten metal in the stalk is obtained at the time t1, and has already been obtained in the past calculation. Further, the viscosity coefficient μ (h) is known from the material and temperature of the molten metal. Then, when these values and the gas pressure P 1 obtained by “Equation 6” are substituted into “Equation 5” and integration is performed twice, the height h 1 of the molten metal surface at the time (t1 + t) of the delivered molten metal. Is required. In this case, when the time (t1 + t) greater than the height the height h 1 is designated by the calculated water surface in lowers the operating voltage of the flow control valve, Stoke by lowering the pressure in the upper space A 15 Control is performed so that the position of the hot water surface in the inside is lowered. On the other hand, when the height h 1 of the hot water surface at time (t1 + t) is lower than the specified height, the operating voltage of the flow control valve is increased and the pressure in the upper space A is raised to increase the temperature of the hot water surface in the stalk 15. Control so that the position is raised. Such control by performing, can be brought close to the height specified the actual height h 1 of the molten steel surface at the time (t1 + t). Incidentally, of decrease of the operation voltage may be proportional to the difference between the height and the specified height h 1 of the molten metal surface.

以上の方法により、ストーク15内での湯面の波立ち(湯面の位置の変動)を抑制することができる。以下、このモデル予測制御について図4〜図9を参照して詳細に説明する。なお、以下に示す例では、現在から3ステップ先まで予測するが、何ステップ先まで予測するかは適宜決定することができるので、以下の例に限定されるものではない。   By the method described above, it is possible to suppress the undulation of the molten metal surface in the stalk 15 (variation in the position of the molten metal surface). Hereinafter, this model prediction control will be described in detail with reference to FIGS. In the example shown below, the prediction is made up to 3 steps ahead from the present, but the number of steps ahead can be determined as appropriate, and is not limited to the following example.

図4はモデル予測制御を説明するためのタイミングチャートであり、横軸はモデル予測制御が行われる各ステップにおける時刻を示す。なお、図4では、1ステップの間隔を0.1の目盛で示している。現在時刻は時刻tで示され、この時刻tのときにコントローラ20からモデル予測コントローラ21に湯面位置を指示する情報が入力される。   FIG. 4 is a timing chart for explaining the model predictive control, and the horizontal axis indicates the time at each step in which the model predictive control is performed. In FIG. 4, the interval of one step is indicated by a scale of 0.1. The current time is indicated by time t, and information indicating the hot water surface position is input from the controller 20 to the model prediction controller 21 at this time t.

湯面の位置を支持する情報が入力されると、モデル予測コントローラ21は、図5に示すように、比例弁17に対する操作電圧を上げて比例弁17を全開となるように作動させる。比例弁17は、作動ヒステリシスを有するため、時刻(t+1)から流路を開き始め、時刻(t+1)から空気が流れ始める。また、モデル予測コントローラ21は、現在の時刻tの時点で圧力センサ18から入力された炉内圧力を前述の「数6」に代入し、時刻(t+1)、時刻(t+2)、時刻(t+3)の時点における予測炉内圧力を得るとともに、得られた予測炉内圧力を「数5」に代入して時刻(t+1)、時刻(t+2)、時刻(t+3)時点でのストーク15内における予測湯面位置を得る。なお、図4には、説明の便宜のために予測炉内圧力と予測湯面位置を時刻(t+3)よりも将来の時刻のものまで示しているが、実際には時刻(t+3)までの予測しかしていない。これについては図5以降の図面も同様である。   When the information that supports the position of the molten metal surface is input, the model prediction controller 21 increases the operating voltage for the proportional valve 17 to activate the proportional valve 17 as shown in FIG. Since the proportional valve 17 has an operation hysteresis, the flow path starts to open from time (t + 1), and air starts to flow from time (t + 1). Further, the model prediction controller 21 substitutes the in-furnace pressure input from the pressure sensor 18 at the time of the current time t into the above-mentioned “Equation 6”, and the time (t + 1), time (t + 2), time (t + 3) Is obtained, and the obtained predicted furnace pressure is substituted into “Equation 5”, and the predicted hot water in the stalk 15 at time (t + 1), time (t + 2), and time (t + 3) is obtained. Get the surface position. In FIG. 4, for convenience of explanation, the predicted furnace pressure and the predicted molten metal surface position are shown up to the future time from the time (t + 3), but actually the prediction up to the time (t + 3) is shown. But not. This also applies to the drawings subsequent to FIG.

図6は、図5に示す状態から次のステップ(1.1)に進んだ状態を示している。モデル予測コントローラ21は、現在の時刻tの時点で圧力センサ18から入力された炉内圧力と前述の「数6」および「数5」から時刻(t+1)、時刻(t+2)、時刻(t+3)の時の予測湯面位置を得る。モデル予測コントローラ21による予測では、現在の時刻tと時刻(t+1)の中間から炉内圧力が上がり始め、その結果、加熱炉11内の溶湯Mの湯面が押され、時刻(t+1)の少し手前からストーク15内に溶湯Mが押し込まれてストーク15内における湯面が上がり始める。この予想される状態においては、予測湯面位置が指示湯面位置よりも大幅に低い。したがって、モデル予測コントローラ21は、比例弁17に印加する操作電圧を維持し、ストーク15内における湯面の上昇を促す。   FIG. 6 shows a state that proceeds from the state shown in FIG. 5 to the next step (1.1). The model prediction controller 21 calculates the time (t + 1), the time (t + 2), and the time (t + 3) from the in-furnace pressure input from the pressure sensor 18 at the time of the current time t and the above-described “Equation 6” and “Equation 5”. The predicted hot water surface position at the time of is obtained. In the prediction by the model prediction controller 21, the furnace pressure starts to increase from the middle of the current time t and time (t + 1), and as a result, the surface of the molten metal M in the heating furnace 11 is pushed, and a little at time (t + 1). The molten metal M is pushed into the stalk 15 from the front, and the molten metal level in the stalk 15 starts to rise. In this expected state, the predicted hot water surface position is significantly lower than the indicated hot water surface position. Therefore, the model prediction controller 21 maintains the operation voltage applied to the proportional valve 17 and promotes the rise of the molten metal level in the stalk 15.

図7は、図6に示す状態から次のステップ(1.2)に進んだ状態を示している。この場合も、モデル予測コントローラ21は、現在の時刻tの時点で圧力センサ18から入力された炉内圧力と前述の「数6」および「数5」から時刻(t+1)、時刻(t+2)、時刻(t+3)の時の予測湯面位置を得る。この際、モデル予測コントローラ21による予測では、時刻(t+2)で予測湯面位置が指示湯面位置に近くなるので、時刻(t+1)で比例弁17に印加する操作電圧を下げる必要があるが、時刻tでは制御電圧を維持する。   FIG. 7 shows a state that proceeds from the state shown in FIG. 6 to the next step (1.2). Also in this case, the model prediction controller 21 calculates the time (t + 1), the time (t + 2), the in-furnace pressure input from the pressure sensor 18 at the time of the current time t, and the above-mentioned “Equation 6” and “Equation 5”. The predicted hot water surface position at time (t + 3) is obtained. At this time, in the prediction by the model prediction controller 21, since the predicted hot water surface position is close to the indicated hot water surface position at time (t + 2), it is necessary to lower the operation voltage applied to the proportional valve 17 at time (t + 1). At time t, the control voltage is maintained.

図8は、図7に示す状態から次のステップ(1.3)に進んだ状態を示している。この状態において、モデル予測コントローラ21は、現在の時刻tの時点で圧力センサ18から入力された炉内圧力と前述の「数4」および「数3」から時刻(t+1)、時刻(t+2)、時刻(t+3)の時の予測湯面位置を得る。この状態におけるモデル予測コントローラ21の予測では、時刻(t+1)で予測湯面位置が指示湯面位置に近くなるので、ストーク15内の湯面をさげるべく、1ステップ前のタイミングである時刻tで比例弁17に印加する操作電圧を下げる。   FIG. 8 shows a state that proceeds from the state shown in FIG. 7 to the next step (1.3). In this state, the model prediction controller 21 determines the time (t + 1), the time (t + 2), the furnace pressure inputted from the pressure sensor 18 at the time of the current time t, and the above-mentioned “Equation 4” and “Equation 3”. The predicted hot water surface position at time (t + 3) is obtained. In the prediction of the model prediction controller 21 in this state, the predicted hot water surface position is close to the indicated hot water surface position at time (t + 1), so that the hot water surface in the stalk 15 is lowered at time t, which is the timing one step before. The operating voltage applied to the proportional valve 17 is lowered.

つまり、図8に示す状況において、予測湯面位置が指示湯面位置に接近する時刻(t+1)の手前の時刻tにおいて、早めに比例弁17を絞り、上部空間Aの圧力を早めに低下させ、ストーク15内の湯面を下げようとする操作を行う。上部空間Aの圧力を変化させると、空気の圧縮性および溶湯の慣性に起因して少し遅れたタイミングでストーク15内の湯面の変化が現れる。よこで、上記のタイミングの遅れを考慮して、少し早いタイミングで圧力変化を与えることで、ストーク15内の湯面の実際の位置が、時刻(t+1)で指示した位置に近接するようにする。   In other words, in the situation shown in FIG. 8, at the time t before the time (t + 1) when the predicted hot water surface position approaches the indicated hot water surface position, the proportional valve 17 is throttled early to reduce the pressure in the upper space A early. Then, an operation for lowering the hot water level in the stalk 15 is performed. When the pressure in the upper space A is changed, a change in the hot water level in the stalk 15 appears at a timing slightly delayed due to the compressibility of the air and the inertia of the molten metal. In consideration of the above timing delay, the pressure change is applied at a slightly earlier timing so that the actual position of the hot water surface in the stalk 15 is close to the position indicated at time (t + 1). .

図9は、図8に示す状態から次のステップ(1.4)に進んだ状態を示している。この状態においても、モデル予測コントローラ21は、現在の時刻tの時点において圧力センサ18から入力された炉内圧力と前述の「数6」および「数5」から時刻(t+1)、時刻(t+2)、時刻(t+3)の時の予測湯面位置を得る。モデル予測コントローラ21の予測では、時刻t以降では、予測湯面位置が指示湯面位置に近似するが、溶湯の慣性に起因して湯面の高さ位置は一定値とならず、上下に揺れる。例えばこの場合、時刻(t+2)で予測湯面位置が指示湯面位置よりも少し低い位置となる。そこで、モデル予測コントローラ21は、時刻tで比例弁17の操作電圧を少し上げ、時刻(t+2)よりも少し早いタイミングでストーク15内の湯面の位置をあげるべく、上部空間Aの圧力を少し上げる。このようにして、予測されるストーク15内の湯面の位置に対応させて、上部空間Aの圧力を小刻みに動的に制御することで、時刻(t+3)以降における予測湯面位置を、指示湯面位置に近接させ、ストーク15内の湯面の位置を指示した位置近傍で安定させる。   FIG. 9 shows a state that proceeds from the state shown in FIG. 8 to the next step (1.4). Even in this state, the model prediction controller 21 determines the time (t + 1) and the time (t + 2) from the in-furnace pressure input from the pressure sensor 18 at the time of the current time t and the above-mentioned “Equation 6” and “Equation 5”. The predicted hot water surface position at time (t + 3) is obtained. In the prediction of the model prediction controller 21, after the time t, the predicted molten metal surface position approximates the indicated molten metal surface position, but due to the inertia of the molten metal, the molten metal surface height position does not become a constant value and fluctuates up and down. . For example, in this case, the predicted hot water surface position is slightly lower than the indicated hot water surface position at time (t + 2). Therefore, the model prediction controller 21 slightly increases the operation voltage of the proportional valve 17 at time t, and slightly increases the pressure in the upper space A to raise the position of the hot water surface in the stalk 15 at a timing slightly earlier than time (t + 2). increase. In this manner, the predicted hot water position after time (t + 3) is indicated by dynamically controlling the pressure in the upper space A in small increments in correspondence with the predicted hot water position in the stalk 15. The hot water surface is brought close to the surface, and the position of the hot water surface in the stalk 15 is stabilized in the vicinity of the designated position.

以上のようなモデル予測制御により、ストーク15内における溶湯Mの湯面の位置を安定させることができる。そして、湯面の位置を安定させることで、ストーク15内での湯面の上昇速度を速くすることができ、鋳造時間を短縮することができる。また、溶湯がキャビティ14に充填されるときも安定状態がほぼ維持されるので、PID制御のように遅れやオーバーシュートの発生を抑制することができ、溶湯の脈動の発生とそれに起因する湯境や未充填、あるいは空気の巻き込みなどの発生を未然に防止することができる。   By the model predictive control as described above, the position of the molten metal surface of the molten metal M in the stalk 15 can be stabilized. Then, by stabilizing the position of the molten metal surface, the rising speed of the molten metal surface in the stalk 15 can be increased, and the casting time can be shortened. Further, since the stable state is substantially maintained even when the molten metal is filled in the cavity 14, it is possible to suppress the occurrence of delay and overshoot as in PID control, and the occurrence of molten metal pulsation and the resulting hot water boundary. It is possible to prevent the occurrence of unfilling or air entrainment.

ところで、溶湯がストーク15内を通過した後も比例弁17の操作電圧を一定にしておくと、湯口やキャビティ14内で水平断面積が大きく変化する場合に、湯面の上昇速度が変動して湯面が乱れることが想定される。この場合、モデル予測コントローラ21がストーク15およびキャビティ14の水平断面積情報を記憶し、これを外乱として操作電圧を補正する。   By the way, if the operating voltage of the proportional valve 17 is kept constant even after the molten metal has passed through the stalk 15, the rise rate of the molten metal surface fluctuates when the horizontal cross-sectional area greatly changes in the gate or the cavity 14. It is assumed that the hot water surface is disturbed. In this case, the model prediction controller 21 stores the horizontal cross-sectional area information of the stalk 15 and the cavity 14 and corrects the operation voltage using this as a disturbance.

具体的には、ストーク15の水平断面積をAとし、湯口やキャビティ14の水平断面積をSとしたときに、操作電圧にS/Aを乗算して補正する。これにより、湯口やキャビティ14内で湯面の上昇速度は一定となり、湯面の乱れが防止される。   Specifically, when the horizontal cross-sectional area of the stalk 15 is A and the horizontal cross-sectional area of the gate or cavity 14 is S, the operation voltage is multiplied by S / A for correction. As a result, the rising speed of the hot water surface in the pouring gate and the cavity 14 becomes constant, and the hot water surface is prevented from being disturbed.

(Vの再算出処理)
上部空間Aの実効体積Vは、炉11内壁の材質の変化等に起因して、炉11の使用に従って変化する。また、リーク量もパッキンの劣化等に起因して変化する。したがって、鋳造を繰り返し行う過程において、「数6」による予測ガス圧力Pの算出結果の精度は低下する。そこで、適切なタイミングで図3に示す処理を再度行うことで、長期に渡りVの値の精度を確保することができ、「数6」の算出結果の精度を維持できる。図3の処理を再度行うタイミングは、鋳造回数をカウントし、鋳造回数の累積値が規定の値となった段階で行う方法、炉の使用時間が一定の累積時間が経過した段階で行う方法等が挙げられる。この処理が自動的に行われるようにプログラムを設定することは有効である。勿論、オペレータの操作により図3の処理を行うタイミングを決めることもできる。
(V recalculation process)
The effective volume V of the upper space A changes according to the use of the furnace 11 due to a change in the material of the inner wall of the furnace 11 or the like. Further, the leak amount also changes due to deterioration of the packing or the like. Therefore, in the process of repeating the casting, the accuracy of the predicted gas pressure P 1 of the calculation result by "6" is reduced. Therefore, by performing the process shown in FIG. 3 again at an appropriate timing, the accuracy of the value of V can be ensured over a long period of time, and the accuracy of the calculation result of “Equation 6” can be maintained. The timing of performing the process of FIG. 3 again is a method of counting the number of castings, a method of performing when the cumulative value of the number of castings reaches a specified value, a method of performing when the furnace has been used for a certain cumulative time, etc. Is mentioned. It is effective to set a program so that this process is automatically performed. Of course, the timing for performing the processing of FIG. 3 can be determined by the operation of the operator.

(優位性)
以上例示した技術によれば、炉内壁へのガスの浸透や不定形な部分の体積を考慮にいれた上部空間Aの実効体積Vを得ることができるので、「数5」および「数6」を用いて行う図5〜図9に示す制御の精度を高めることができる。仮に、Vの値として、設計値の値を用いると、そこには大きな誤差が含まれるので、「数6」により求めたPの値の誤差が大きくなり、「数5」により算出されるストーク15内の湯面の高さ位置の予測値も大きな誤差を含むことになる。そして、図5〜図9に示す制御の精度も低下し、低圧鋳造における溶湯の脈動の発生とそれに起因する湯境や未充填、あるいは空気の巻き込みなどの発生を防ぐ効果が低下する。
(Superiority)
According to the techniques exemplified above, the effective volume V of the upper space A can be obtained in consideration of the gas permeation into the furnace inner wall and the volume of the irregular portion. Therefore, “Equation 5” and “Equation 6” The accuracy of the control shown in FIGS. If, as the value and V, the use of the value of the design value, since there is included a large error, the error value of P 1 obtained by "6" increases, is calculated by the "Number 5" The predicted value of the height position of the hot water surface in the stalk 15 also includes a large error. And the precision of control shown in Drawing 5-Drawing 9 also falls, and the effect which prevents generation | occurrence | production of the pulsation of the molten metal in low-pressure casting, and the hot water boundary and unfilling resulting from it, or entrainment of air, etc. falls.

また、「数6」を用いてのPの算出において、上部空間Aにおける外部へのリークの影響が取り込まれているので、このリークの影響を取り込まずに計算を行う場合に比較して、Pの予測精度および「数5」から算出されるストーク15内における湯面の予測位置の精度が向上する。このことも、図5〜図9に示す制御の精度を向上させ、低圧鋳造における溶湯の脈動の発生とそれに起因する湯境や未充填、あるいは空気の巻き込みなどの発生を防ぐ効果を高くすることに寄与する。 In addition, in the calculation of P 1 using “Equation 6”, the influence of the leak to the outside in the upper space A is taken in. Compared to the case of performing the calculation without taking in the influence of this leak, The prediction accuracy of P 1 and the accuracy of the predicted position of the molten metal surface in the stalk 15 calculated from “Equation 5” are improved. This also improves the accuracy of the control shown in FIGS. 5 to 9 and enhances the effect of preventing the occurrence of molten metal pulsation and the occurrence of molten metal boundary, unfilled, or air entrainment in low pressure casting. Contribute to.

本発明は、鋳造技術に用いることができる。   The present invention can be used in casting techniques.

10…低圧鋳造装置、11…炉、12…下型、13…上型、14…キャビティ、15…ストーク、16…配管、A…上部空間(炉内における溶湯の上の空間)、M…溶湯。   DESCRIPTION OF SYMBOLS 10 ... Low pressure casting apparatus, 11 ... Furnace, 12 ... Lower mold, 13 ... Upper mold, 14 ... Cavity, 15 ... Stoke, 16 ... Pipe, A ... Upper space (space above the molten metal in the furnace), M ... Molten metal .

Claims (7)

炉内の空気室を大気圧よりも高い所定の圧力にした際における前記空気室からの気体のリーク量を計測するリーク量計測工程と、
所定の流量の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から
前記所定の圧力になるまでに掛かる経過時間を計測する経過時間計測工程と、
前記リーク量、前記計測された前記経過時間、前記経過時間計測工程において前記空気室に供給された前記気体の流量、前記所定圧力および前記空気室の温度に基づき、気体の状態方程式を利用して前記空気室の体積を算出する空気室体積算出工程と
を有することを特徴とする炉の空気室の体積を算出する方法。
A leak amount measuring step of measuring a leak amount of gas from the air chamber when the air chamber in the furnace is set to a predetermined pressure higher than the atmospheric pressure;
An elapsed time measuring step of measuring an elapsed time required for the air chamber to change from atmospheric pressure to the predetermined pressure when a gas having a predetermined flow rate is sent to the air chamber at atmospheric pressure;
Based on the amount of leak, the measured elapsed time, the flow rate of the gas supplied to the air chamber in the elapsed time measuring step, the predetermined pressure, and the temperature of the air chamber, using a gas equation of state An air chamber volume calculating step for calculating the volume of the air chamber. A method for calculating the volume of the air chamber of the furnace.
前記空気室体積算出工程において、下記「数1」に基づいて前記空気室の体積Vの算出が行われることを特徴とする請求項1に記載の炉の空気室の体積を算出する方法。
Figure 0005912547
P:空気室の圧力、GIn:空気室への気体の流入量、GLeak:空気室からの気体のリーク量
V:空気室の体積、T:空気室の温度、R:気体定数
The method of calculating the volume of the air chamber of the furnace according to claim 1, wherein, in the air chamber volume calculating step, the volume V of the air chamber is calculated based on the following "Equation 1".
Figure 0005912547
P: pressure of the air chamber, G In : amount of gas flowing into the air chamber, G Leak : amount of gas leakage from the air chamber V: volume of the air chamber, T: temperature of the air chamber, R: gas constant
前記炉を用いた鋳造を行う前と後に前記空気室体積算出工程を行うことを特徴とする請求項1または2に記載の炉の空気室の体積を算出する方法。   The method for calculating the volume of the air chamber of the furnace according to claim 1, wherein the air chamber volume calculating step is performed before and after casting using the furnace. 請求項1〜3のいずれか一項に記載の発明によって算出した前記炉の前記空気室の体積に基づいて、前記空気室に流入する気体による圧力増加分と、前記空気室に流入した前記気体により前記空気室において溶湯が押し下げられる結果、前記空気室の容積が増加することによる圧力減少分と、前記空気室からのリークによって生じる圧力減少分とを算出する圧力変化分算出工程と、
前記圧力変化分算出工程における演算の結果に基づき、前記炉からストークを介して送り出される溶湯の湯面の高さ位置を算出する湯面位置算出工程と
を有することを特徴とする鋳造方法。
Based on the volume of the air chamber of the furnace calculated by the invention according to any one of claims 1 to 3, the pressure increase due to the gas flowing into the air chamber and the gas flowing into the air chamber As a result of the molten metal being pushed down in the air chamber, a pressure change calculation step for calculating a pressure decrease due to an increase in the volume of the air chamber and a pressure decrease caused by a leak from the air chamber;
And a molten metal surface position calculating step of calculating the height position of the molten metal surface fed from the furnace via the stalk based on the calculation result in the pressure change component calculating step.
前記圧力変化分算出工程が下記「数2」に基づいて行われ、
前記湯面位置算出工程が下記「数3」に基づいて行われることを特徴とする請求項4に記載の鋳造方法。
Figure 0005912547
P:空気室の予測ガス圧力、P:空気室の実測ガス圧力、V:空気室の体積、R:気体定数
T:空気室の温度、G:空気室に流入するガス流量、GLeak:空気室からのガスリーク量
Figure 0005912547
P:空気室の予測ガス圧力、S:ストーク内部の水平断面積、Pi:溶湯の湯面にかかる背圧
A:ストーク内部の溶湯の質量、g:重力加速度、h:送り出された溶湯の湯面の高さ
μ(h):溶湯の粘性係数
The pressure change calculation step is performed based on the following “Equation 2”,
The casting method according to claim 4, wherein the molten metal surface position calculating step is performed based on the following “Equation 3”.
Figure 0005912547
P 1 : Predicted gas pressure in the air chamber, P: Actual gas pressure in the air chamber, V: Volume of the air chamber, R: Gas constant T: Temperature of the air chamber, G: Gas flow rate flowing into the air chamber, G Leak : Gas leakage from the air chamber
Figure 0005912547
P 1 : Predicted gas pressure in the air chamber, S: Horizontal cross-sectional area inside the stalk, P i : Back pressure applied to the molten metal surface, m A : Mass of molten metal inside the stalk, g: Gravitational acceleration, h 1 : Delivery Molten metal surface height μ (h): Melt viscosity coefficient
炉内の空気室を大気圧よりも高い所定の圧力にした際における前記空気室からの気体のリーク量を計測するリーク量計測手段と、
所定の流量の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から
前記所定の圧力になるまでに掛かる経過時間を計測する経過時間計測手段と、
前記リーク量、前記計測された前記経過時間、前記気体の前記所定の流量、前記所定圧力および前記空気室の温度に基づき、気体の状態方程式を利用して前記空気室の体積を算出する空気室体積算出手段と
を有することを特徴とする炉の空気室の体積を算出する装置。
A leak amount measuring means for measuring a leak amount of gas from the air chamber when the air chamber in the furnace is set to a predetermined pressure higher than the atmospheric pressure;
An elapsed time measuring means for measuring an elapsed time required for the air chamber to change from atmospheric pressure to the predetermined pressure when a gas having a predetermined flow rate is sent to the air chamber at atmospheric pressure;
An air chamber that calculates the volume of the air chamber using a gas equation of state based on the leak amount, the measured elapsed time, the predetermined flow rate of the gas, the predetermined pressure, and the temperature of the air chamber A device for calculating the volume of the air chamber of the furnace.
コンピュータに読み取らせて実行させるプログラムであって、
CPUに、
炉内の空気圧を大気圧よりも高い所定の圧力保つようにした際に、前記空気室への空気の流入量を計測することで前記空気室からの気体のリーク量 Leak を計測するリーク量計測ステップと、
前記リーク量の計測後に前記空気室を大気圧に戻し、前記リーク量G Leak より多い流量G In で前記空気室に前記圧力Pになるまで空気を流入させるように比例弁を制御する制御ステップと、
所定の流量G In の気体を大気圧とされた前記空気室に送った際に、前記空気室が大気圧から前記所定の圧力Pになるまでに掛かる経過時間の範囲で下記「数4」の計算を行い、前記空気室の体積Vを算出する空気室体積算出ステップと
を実行させることを特徴とする炉の空気室の体積を算出するためのプログラム。
Figure 0005912547
T:空気室の温度、R:気体定数
A program that is read and executed by a computer,
CPU
When the air pressure in the furnace is kept at a predetermined pressure P higher than the atmospheric pressure, the amount of air leakage G Leak from the air chamber is measured by measuring the amount of air flowing into the air chamber. A leak amount measuring step ;
A control step of controlling the proportional valve so that the air chamber is returned to atmospheric pressure after the leakage amount is measured and air is allowed to flow into the air chamber at the flow rate G In higher than the leakage amount G Leak until the pressure P is reached. ,
When a gas having a predetermined flow rate G In is sent to the air chamber that is set to atmospheric pressure, the following “Equation 4” is satisfied within the range of elapsed time that the air chamber takes from atmospheric pressure to the predetermined pressure P. An air chamber volume calculating step for calculating and calculating a volume V of the air chamber;
Program for calculating the volume of the air chamber of the furnace, characterized in that for the execution.
Figure 0005912547
T: temperature of air chamber, R: gas constant
JP2012003052A 2012-01-11 2012-01-11 Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber Expired - Fee Related JP5912547B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012003052A JP5912547B2 (en) 2012-01-11 2012-01-11 Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012003052A JP5912547B2 (en) 2012-01-11 2012-01-11 Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber

Publications (2)

Publication Number Publication Date
JP2013141691A JP2013141691A (en) 2013-07-22
JP5912547B2 true JP5912547B2 (en) 2016-04-27

Family

ID=49038547

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012003052A Expired - Fee Related JP5912547B2 (en) 2012-01-11 2012-01-11 Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber

Country Status (1)

Country Link
JP (1) JP5912547B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019010513A1 (en) * 2017-07-13 2019-01-17 Fill Gesellschaft M.B.H. Casting device for pressure casting

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104889364B (en) * 2015-06-03 2018-01-05 天津市航宇嘉瑞科技股份有限公司 A kind of High-Voltage Electrical Appliances accessory atmospheric air chamber product mold
CN108380846B (en) * 2018-03-21 2021-05-07 海纳川(滨州)轻量化汽车部件有限公司 Riser tube leak test detection method and device

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62248551A (en) * 1986-04-19 1987-10-29 Mazda Motor Corp Molten metal supplying control apparatus for low pressure casting
JP2615094B2 (en) * 1987-11-27 1997-05-28 マツダ株式会社 Hot water supply control device for water heater
JP3268304B2 (en) * 1994-04-04 2002-03-25 愛知機械工業株式会社 Pressure control device in low pressure casting machine
JP3038321B2 (en) * 1997-04-21 2000-05-08 小川技研株式会社 Leak amount measurement device
JP2001041843A (en) * 1999-08-02 2001-02-16 Hitachi Ltd Leak checker for sealed vessel
JP5339987B2 (en) * 2009-03-27 2013-11-13 本田技研工業株式会社 Casting machine control device

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019010513A1 (en) * 2017-07-13 2019-01-17 Fill Gesellschaft M.B.H. Casting device for pressure casting

Also Published As

Publication number Publication date
JP2013141691A (en) 2013-07-22

Similar Documents

Publication Publication Date Title
US8875960B2 (en) Tilting-type automatic molten metal pouring method, tilting control system, and storage medium having tilting control program stored therein
CN105073305B (en) Pour into a mould control method and be stored with for making computer as the storage medium of the program of cast control unit function
JP5912547B2 (en) Method for calculating volume of furnace air chamber, casting method, apparatus for calculating volume of furnace air chamber, and program for calculating volume of furnace air chamber
US10705546B2 (en) Flow rate control apparatus, flow rate control method, and program recording medium
KR102483417B1 (en) Flow rate control apparatus, flow rate control method, and program recording medium
JP6607215B2 (en) Flow state estimation method for molten steel, flow state estimation device, on-line display device for flow state of molten steel, and continuous casting method for steel
KR20100017355A (en) Automatic pouring control method, control system of servo motor of automatic pouring device and medium storing tilting control program for ladle
JP7059053B2 (en) Flow control device, flow control method, and program for flow control device
JPWO2019163676A1 (en) Flow control device and flow control method
TWI695244B (en) Fluid control system and flow measurement method
JP2010282243A (en) Simulation device for vacuum container
JP2009241139A (en) Forecasting method for molten steel temperature within tundish, and management method
JP5339987B2 (en) Casting machine control device
US9993868B2 (en) Control device for the advancing motion of a casting plunger
JP5998914B2 (en) Model predictive control method for level of molten metal in continuous casting equipment.
US5031805A (en) Processes and device for dosing free-flowing media
JP2010037770A (en) Inflow rate predicting method and dam management control system
TWI552816B (en) Active liquid level control system and control method thereof
JP6217324B2 (en) Method of controlling the height of the pouring bath of the molten metal holding furnace for low pressure casting
Jabri et al. Suppression of periodic disturbances in the continuous casting process
JP5998919B2 (en) Mold level control method in mold
JP7474541B1 (en) Flow Control Device
KR20190002610A (en) Pressure type flow control device, flow rate calculation method and flow rate control method thereof
JPH07256420A (en) Device for controlling molten metal surface level in continuous caster
JPH01170568A (en) Method for controlling molten metal surface level

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141127

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150727

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150730

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150928

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160303

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160401

R150 Certificate of patent or registration of utility model

Ref document number: 5912547

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees