JP6874485B2

JP6874485B2 - Extruded product internal temperature distribution estimation device and extrusion product internal temperature distribution estimation method during hot extrusion, billet heating control device and billet heating control method and program

Info

Publication number: JP6874485B2
Application number: JP2017078035A
Authority: JP
Inventors: 真友岸; 純一郎山口; 淳一西森
Original assignee: Nippon Steel Corp
Current assignee: Nippon Steel Corp
Priority date: 2017-04-11
Filing date: 2017-04-11
Publication date: 2021-05-19
Anticipated expiration: 2037-04-11
Also published as: JP2018176210A

Description

本発明は、熱間押出加工時の押出製品内部温度分布推定装置および押出製品内部温度分布推定方法、ビレット加熱制御装置およびビレット加熱制御方法ならびにプログラムに関する。 The present invention relates to an extruded product internal temperature distribution estimation device and an extruded product internal temperature distribution estimation method, a billet heating control device, a billet heating control method, and a program during hot extrusion processing.

炭素鋼、ステンレス鋼、アルミニウム合金または銅合金等の金属加工品の製造方法の１つとして、熱間押出加工がしばしば用いられている。熱間押出加工においては、素材となるビレットを加熱した後、コンテナ内に挿入し、ビレットをラムでダイスに向けて押出すことで、所望の形状の押出製品が得られる。 Hot extrusion is often used as one of the methods for producing processed metal products such as carbon steel, stainless steel, aluminum alloys and copper alloys. In hot extrusion processing, a billet as a material is heated, inserted into a container, and the billet is extruded toward a die with a ram to obtain an extruded product having a desired shape.

熱間押出加工を行うに際して、押出製品温度の過度な温度上昇に伴う、横切れと呼ばれる押出製品の外内面におけるキズの発生が問題となる。 When hot extrusion is performed, there is a problem of scratches on the outer and inner surfaces of the extruded product, which is called cross-cutting, due to an excessive rise in the temperature of the extruded product.

例えば、特許文献１においては、押出の進行に伴って、ラムからビレットに対して加えられる押出荷重（ラム荷重）が一定の割合で減少するように、ビレットの設定温度を制御することによって、ダイスから押し出される押出製品の温度を一定にして、かかる押出材の品質を、高度に且つ安定して維持することのできる押出方法が開示されている。 For example, in Patent Document 1, the die is controlled by controlling the set temperature of the billet so that the extrusion load (ram load) applied from the ram to the billet decreases at a constant rate as the extrusion progresses. There is disclosed an extrusion method capable of maintaining the quality of the extruded material highly and stably by keeping the temperature of the extruded product extruded from the extruded product constant.

特開２００９−２４８１８８号公報JP-A-2009-248188

ところで、熱間押出加工における操業上の問題の１つとして、押出初期における押出速度の急変（オーバーシュート）および押出力の増大が挙げられる。押出初期にオーバーシュートが生じると、加工発熱量が過大となり、押出製品の外内面に横切れを生じさせる結果となる。また、押出力の増大により設備制約を超過した場合、加工途中詰まり（半押し）または最悪の場合は設備が破壊されるなどの大きなトラブルを引き起こすことになる。 By the way, one of the operational problems in hot extrusion processing is a sudden change in extrusion speed (overshoot) and an increase in extrusion output at the initial stage of extrusion. If overshoot occurs at the initial stage of extrusion, the amount of heat generated by processing becomes excessive, resulting in cross-cutting on the outer and inner surfaces of the extruded product. In addition, if the equipment restrictions are exceeded due to an increase in push output, major troubles such as clogging during processing (half-push) or, in the worst case, equipment destruction will occur.

横切れ発生リスクを抑制するためには、加熱温度の低温化が有効であり、一方、最大押出力を低減するには、加熱温度の高温化が有効である。すなわち、両者はトレードオフの関係にあり、最適な操業条件を見出すためには試行錯誤が必要となる。 In order to suppress the risk of cross-cutting, it is effective to lower the heating temperature, while in order to reduce the maximum push output, it is effective to raise the heating temperature. That is, there is a trade-off relationship between the two, and trial and error is required to find the optimum operating conditions.

そのため、効率よく最適な加熱条件を見出すためには、押出加工後の押出製品の温度を正確に把握できることが重要である。しかし、押出製品の内部温度を正確に測定することは困難であり、表面温度の測定結果から内部温度分布を推定することにも技術的なハードルが存在する。 Therefore, in order to efficiently find the optimum heating conditions, it is important to be able to accurately grasp the temperature of the extruded product after extrusion. However, it is difficult to accurately measure the internal temperature of an extruded product, and there is a technical hurdle in estimating the internal temperature distribution from the measurement result of the surface temperature.

そこで、押出製品の内部温度分布を推定する方法が求められている。さらに、その推定結果に基づいてビレット加熱条件を適正化することで、品質面（外内面キズ抑制）および操業面（安定操業）での改善効果が期待できる。 Therefore, there is a demand for a method for estimating the internal temperature distribution of extruded products. Furthermore, by optimizing the billet heating conditions based on the estimation result, improvement effects in terms of quality (suppression of scratches on the outer and inner surfaces) and operation (stable operation) can be expected.

特許文献１によれば、熱間加工時の製品の温度を一定に保つことができるとされているが、内部温度分布を推定する方法については一切開示がなされていない。 According to Patent Document 1, it is said that the temperature of the product during hot working can be kept constant, but the method of estimating the internal temperature distribution is not disclosed at all.

本発明は上記の問題を解決し、熱間押出加工時における押出製品内部の温度分布を高精度で推定する装置および方法を提供することを目的とする。 An object of the present invention is to solve the above problems and to provide an apparatus and a method for estimating the temperature distribution inside an extruded product with high accuracy during hot extrusion processing.

本発明は、上記の課題を解決するためになされたものであり、下記の熱間押出加工時の押出製品内部温度分布推定装置および押出製品内部温度分布推定方法、ビレット加熱制御装置およびビレット加熱制御方法ならびにプログラムを要旨とする。 The present invention has been made to solve the above problems, and is described below in an extruded product internal temperature distribution estimation device and an extruded product internal temperature distribution estimation method, a billet heating control device, and a billet heating control during hot extrusion processing. The gist is the method and program.

（１）ビレットに対して所定の押出速度および押出力の条件で熱間押出加工を施し、押出製品を製造した際の、前記押出製品の長手方向における内部の温度分布を推定する装置であって、
前記押出製品の長手方向の温度分布を複数設定する、設定部と、
前記設定部によって設定された前記複数の温度分布ごとに、前記温度分布ならびに予め定められた押出速度変化目標値および押出力変化目標値を含む情報に基づいて、複数の押出速度推定値および複数の押出力推定値の時系列データを算出する、計算部と、
前記計算部によって得られた、前記複数の押出速度推定値および前記複数の押出力推定値の時系列データを、押出速度実測値および押出力実測値の時系列データと比較し、押出速度推定値および押出力推定値の最適な時系列データに対応する１つの温度分布を抽出する、抽出部と、
前記抽出部によって抽出された前記１つの温度分布を、前記押出製品の長手方向における内部の温度分布に決定する、決定部と、を備える、
熱間押出加工時の押出製品内部温度分布推定装置。 (1) An apparatus for estimating the internal temperature distribution in the longitudinal direction of an extruded product when a billet is hot-extruded under predetermined extrusion speed and extrusion power conditions to manufacture an extruded product. ,
A setting unit that sets a plurality of temperature distributions in the longitudinal direction of the extruded product, and
For each of the plurality of temperature distributions set by the setting unit, a plurality of extrusion speed estimation values and a plurality of extrusion speed estimation values are based on the information including the temperature distribution and predetermined extrusion speed change target values and extrusion output change target values. A calculation unit that calculates time-series data of extrusion output estimates,
The time series data of the plurality of extrusion speed estimated values and the plurality of push output estimated values obtained by the calculation unit are compared with the time series data of the extrusion speed measured value and the push output measured value, and the extrusion speed estimated value is compared. And an extraction unit that extracts one temperature distribution corresponding to the optimum time series data of the push output estimate,
A determination unit is provided, which determines the temperature distribution of the one extracted by the extraction unit as the internal temperature distribution in the longitudinal direction of the extruded product.
Extruded product internal temperature distribution estimation device during hot extrusion processing.

（２）前記設定部は、前記複数の温度分布の設定を複数回行い、２回目以降の前記温度分布の設定を、前記抽出部によって抽出された前記１つの温度分布に基づいて行う、
上記（１）に記載の熱間押出加工時の押出製品内部温度分布推定装置。 (2) The setting unit sets the plurality of temperature distributions a plurality of times, and sets the temperature distribution from the second time onward based on the one temperature distribution extracted by the extraction unit.
The device for estimating the internal temperature distribution of an extruded product during hot extrusion processing according to (1) above.

（３）前記設定部は、前記２回目以降の前記温度分布の設定を、粒子群最適化手法を用いて行う、
上記（２）に記載の熱間押出加工時の押出製品内部温度分布推定装置。 (3) The setting unit sets the temperature distribution from the second time onward by using a particle swarm optimization method.
The device for estimating the internal temperature distribution of an extruded product during hot extrusion processing according to (2) above.

（４）ビレットに対して所定の押出速度および押出力の条件で熱間押出加工を施し、押出製品を製造した際の、前記押出製品の長手方向における内部の温度分布を推定する方法であって、
（ａ）前記押出製品の長手方向の温度分布を複数設定するステップと、
（ｂ）前記（ａ）のステップにおいて設定された前記複数の温度分布ごとに、前記温度分布ならびに予め定められた押出速度変化目標値および押出力変化目標値を含む情報に基づいて、複数の押出速度推定値および複数の押出力推定値の時系列データを算出するステップと、
（ｃ）前記（ｂ）のステップにおいて得られた、前記複数の押出速度推定値および前記複数の押出力推定値の時系列データを、押出速度実測値および押出力実測値の時系列データと比較し、押出速度推定値および押出力推定値の最適な時系列データに対応する１つの温度分布を抽出するステップと、
（ｄ）前記（ｃ）のステップにおいて抽出された前記１つの温度分布を、前記押出製品の長手方向における内部の温度分布に決定するステップと、を備える、
熱間押出加工時の押出製品内部温度分布推定方法。 (4) A method of estimating the internal temperature distribution in the longitudinal direction of the extruded product when the billet is hot-extruded under predetermined extrusion speed and extrusion power conditions to manufacture the extruded product. ,
(A) A step of setting a plurality of temperature distributions in the longitudinal direction of the extruded product, and
(B) For each of the plurality of temperature distributions set in the step (a), a plurality of extrusions are performed based on the temperature distribution and information including a predetermined extrusion speed change target value and push output change target value. Steps to calculate time series data for velocity estimates and multiple extrusion estimates,
(C) Compare the time series data of the plurality of extrusion speed estimates and the plurality of push output estimates obtained in the step (b) with the time series data of the extrusion speed measurement value and the push output measurement value. Then, the step of extracting one temperature distribution corresponding to the optimum time series data of the extrusion speed estimate and the push output estimate, and
(D) The step comprises determining the temperature distribution of the one extracted in the step (c) to the internal temperature distribution in the longitudinal direction of the extruded product.
A method for estimating the internal temperature distribution of an extruded product during hot extrusion.

（５）前記（ａ）のステップにおいて、前記複数の温度分布の設定を複数回行い、２回目以降の前記温度分布の設定を、前記抽出部によって抽出された前記１つの温度分布に基づいて行う、
上記（４）に記載の熱間押出加工時の押出製品内部温度分布推定方法。 (5) In the step (a), the plurality of temperature distributions are set a plurality of times, and the second and subsequent temperature distributions are set based on the one temperature distribution extracted by the extraction unit. ,
The method for estimating the internal temperature distribution of an extruded product during hot extrusion processing according to (4) above.

（６）前記（ａ）のステップにおいて、前記２回目以降の前記温度分布の設定を、粒子群最適化手法を用いて行う、
上記（５）に記載の熱間押出加工時の押出製品内部温度分布推定方法。 (6) In the step (a), the temperature distribution is set from the second time onward by using the particle swarm optimization method.
The method for estimating the internal temperature distribution of an extruded product during hot extrusion processing according to (5) above.

（７）熱間押出加工時のビレットの加熱温度を制御する装置であって、
上記（１）から（３）までのいずれかに記載される押出製品内部温度分布推定装置と、
温度制御部と、を備え、
前記温度制御部は、前記押出製品内部温度分布推定装置が推定した前記押出製品の温度分布、ならびに前記押出速度変化目標値および前記押出力変化目標値を含む情報を用いて算出される、押出速度推定値および押出力推定値の時系列データに基づき、
前記押出速度推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を下げ、
前記押出力推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を上げる、
ビレット加熱制御装置。 (7) A device that controls the heating temperature of billets during hot extrusion.
The extruded product internal temperature distribution estimation device according to any one of (1) to (3) above, and
Equipped with a temperature control unit
The temperature control unit calculates the extrusion speed using information including the temperature distribution of the extrusion product estimated by the extrusion product internal temperature distribution estimation device, the extrusion speed change target value, and the extrusion output change target value. Based on time series data of estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
Billet heating control device.

（８）熱間押出加工時のビレットの加熱温度を制御する方法であって、
上記（４）から（６）までのいずれかに記載される（ａ）〜（ｄ）のステップと、
（ｅ）ビレットの加熱温度を変更するステップと、を備え、
前記（ｅ）のステップにおいて、前記（ｄ）のステップで決定された前記押出製品の温度分布、ならびに前記押出速度変化目標値および前記押出力変化目標値を含む情報を用いて算出される、押出速度推定値および押出力推定値の時系列データに基づき、
前記押出速度推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を下げ、
前記押出力推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を上げる、
ビレット加熱制御方法。 (8) A method of controlling the heating temperature of billets during hot extrusion.
The steps (a) to (d) described in any of the above (4) to (6) and
(E) The step of changing the heating temperature of the billet is provided.
In step (e), extrusion calculated using information including the temperature distribution of the extruded product determined in step (d), the extrusion speed change target value, and the extrusion output change target value. Based on time series data of velocity estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
Billet heating control method.

（９）コンピュータによって、ビレットに対して所定の押出速度および押出力の条件で熱間押出加工を施し、押出製品を製造した際の、前記押出製品の長手方向における内部の温度分布を推定するためのプログラムであって、
前記コンピュータに、
（ａ）前記押出製品の長手方向の温度分布を複数設定するステップと、
（ｂ）前記（ａ）のステップにおいて設定された前記複数の温度分布ごとに、前記温度分布ならびに予め定められた押出速度変化目標値および押出力変化目標値を含む情報に基づいて、複数の押出速度推定値および複数の押出力推定値の時系列データを算出するステップと、
（ｃ）前記（ｂ）のステップにおいて得られた、前記複数の押出速度推定値および前記複数の押出力推定値の時系列データを、押出速度実測値および押出力実測値の時系列データと比較し、押出速度推定値および押出力推定値の最適な時系列データに対応する１つの温度分布を抽出するステップと、
（ｄ）前記（ｃ）のステップにおいて抽出された前記１つの温度分布を、前記押出製品の長手方向における内部の温度分布に決定するステップと、を実行させる、
プログラム。 (9) To estimate the internal temperature distribution in the longitudinal direction of the extruded product when the billet is hot-extruded under predetermined extrusion speed and extrusion power conditions by a computer to manufacture the extruded product. Program
On the computer
(A) A step of setting a plurality of temperature distributions in the longitudinal direction of the extruded product, and
(B) For each of the plurality of temperature distributions set in the step (a), a plurality of extrusions are performed based on the temperature distribution and information including a predetermined extrusion speed change target value and push output change target value. Steps to calculate time series data for velocity estimates and multiple extrusion estimates,
(C) Compare the time series data of the plurality of extrusion speed estimates and the plurality of push output estimates obtained in the step (b) with the time series data of the extrusion speed measurement value and the push output measurement value. Then, the step of extracting one temperature distribution corresponding to the optimum time series data of the extrusion speed estimate and the push output estimate, and
(D) The step of determining the temperature distribution inside the extruded product in the longitudinal direction of the one temperature distribution extracted in the step (c) is executed.
program.

（１０）前記（ａ）のステップにおいて、前記複数の温度分布の設定を複数回行い、２回目以降の前記温度分布の設定を、前記抽出部によって抽出された前記１つの温度分布に基づいて行う、
上記（９）に記載のプログラム。 (10) In the step (a), the plurality of temperature distributions are set a plurality of times, and the second and subsequent temperature distributions are set based on the one temperature distribution extracted by the extraction unit. ,
The program described in (9) above.

（１１）前記（ａ）のステップにおいて、前記２回目以降の前記温度分布の設定を、粒子群最適化手法を用いて行う、
上記（１０）に記載のプログラム。 (11) In the step (a), the temperature distribution is set from the second time onward by using the particle swarm optimization method.
The program according to (10) above.

（１２）コンピュータによって、熱間押出加工時のビレットの加熱温度を制御するためのプログラムであって、
前記コンピュータに、
上記（９）から（１１）までのいずれかに記載される（ａ）〜（ｄ）のステップと、
（ｅ）ビレットの加熱温度を変更するステップと、を実行させ、
前記（ｅ）のステップにおいて、前記（ｄ）のステップで決定された前記押出製品の温度分布、ならびに前記押出速度変化目標値および前記押出力変化目標値を含む情報を用いて算出される、押出速度推定値および押出力推定値の時系列データに基づき、
前記押出速度推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を下げ、
前記押出力推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を上げる、
プログラム。 (12) A program for controlling the heating temperature of billets during hot extrusion by a computer.
On the computer
The steps (a) to (d) described in any of the above (9) to (11), and
(E) Perform the steps of changing the heating temperature of the billet,
In step (e), extrusion calculated using information including the temperature distribution of the extruded product determined in step (d), the extrusion speed change target value, and the extrusion output change target value. Based on time series data of velocity estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
program.

本発明によれば、熱間押出加工時における押出製品内部の温度分布を高精度で推定することが可能となる。また、推定された押出製品内部の温度分布に基づいて、ビレットの加熱温度を適切に制御することが可能になる。 According to the present invention, it is possible to estimate the temperature distribution inside the extruded product during hot extrusion with high accuracy. In addition, the heating temperature of the billet can be appropriately controlled based on the estimated temperature distribution inside the extruded product.

図１は、本発明の一実施形態に係る押出製品内部温度分布推定装置の概略構成を示す図である。FIG. 1 is a diagram showing a schematic configuration of an extruded product internal temperature distribution estimation device according to an embodiment of the present invention. 図２は、押出製品内部温度分布推定装置を備えた押出機の概略構成を示す図である。FIG. 2 is a diagram showing a schematic configuration of an extruder provided with an extrusion product internal temperature distribution estimation device. 図３は、シミュレーションモデルの構成の一例を説明するための図である。FIG. 3 is a diagram for explaining an example of the configuration of the simulation model. 図４は、本発明の一実施形態に係るビレット加熱制御装置の概略構成を示す図である。FIG. 4 is a diagram showing a schematic configuration of a billet heating control device according to an embodiment of the present invention. 図５は、ビレット加熱制御装置を備えた押出機の概略構成を示す図である。FIG. 5 is a diagram showing a schematic configuration of an extruder provided with a billet heating control device. 図６は、本発明の一実施形態に係る押出製品内部温度分布推定方法を示すフロー図である。FIG. 6 is a flow chart showing a method for estimating the internal temperature distribution of an extruded product according to an embodiment of the present invention. 図７は、本発明の一実施形態に係るビレット加熱制御方法を示すフロー図である。FIG. 7 is a flow chart showing a billet heating control method according to an embodiment of the present invention.

本発明の一実施形態に係る熱間押出加工時の押出製品内部温度分布推定装置および押出製品内部温度分布推定方法、ビレット加熱制御装置およびビレット加熱制御方法ならびにプログラムについて、図１〜７を参照しながら説明する。 Refer to FIGS. 1 to 7 for an extrusion product internal temperature distribution estimation device and an extrusion product internal temperature distribution estimation method, a billet heating control device, a billet heating control method, and a program during hot extrusion processing according to an embodiment of the present invention. I will explain while.

図１は、本発明の一実施形態に係る押出製品内部温度分布推定装置の概略構成を示す図である。押出製品内部温度分布推定装置は、ビレットに対して所定の押出速度および押出力の条件で熱間押出加工を施し、押出製品を製造した際の、押出製品の長手方向における内部の温度分布を推定する装置である。 FIG. 1 is a diagram showing a schematic configuration of an extruded product internal temperature distribution estimation device according to an embodiment of the present invention. The extruded product internal temperature distribution estimation device estimates the internal temperature distribution in the longitudinal direction of the extruded product when the billet is hot-extruded under predetermined extrusion speed and extrusion power conditions to manufacture the extruded product. It is a device to do.

ここで、本発明において、押出製品内部の温度とは、押出製品の長手方向に垂直な断面における平均温度を意味し、温度分布とは、押出製品の先端からの距離ごとの当該平均温度の分布を意味する。 Here, in the present invention, the temperature inside the extruded product means the average temperature in the cross section perpendicular to the longitudinal direction of the extruded product, and the temperature distribution is the distribution of the average temperature for each distance from the tip of the extruded product. Means.

また、図２は、押出製品内部温度分布推定装置１０を備えた押出機２００の概略構成を示す図である。図２に示すように、押出機２００は、押出製品内部温度分布推定装置１０に加えて、メインシリンダー２１０、ピアサシリンダー２１５、メインラム２２０、ピアサラム２２５、マンドレル２３０、コンテナ２４０、ダイス２５０、ダイスホルダー２５５、メインバルブ２６０、ピアサバルブ２６５、メイン圧力検出器２７０、ピアサ圧力検出器２７５、および速度検出器２８０を備える。 Further, FIG. 2 is a diagram showing a schematic configuration of an extruder 200 provided with an extrusion product internal temperature distribution estimation device 10. As shown in FIG. 2, in addition to the extrusion product internal temperature distribution estimation device 10, the extruder 200 includes a main cylinder 210, a piercer cylinder 215, a main ram 220, a piercer ram 225, a mandrel 230, a container 240, a die 250, and a die. It includes a holder 255, a main valve 260, a piercer valve 265, a main pressure detector 270, a piercer pressure detector 275, and a speed detector 280.

熱間押出加工を行うに際しては、まず、図示しない中空状のビレットを加熱後、コンテナ２４０内に挿入し、ビレットの中空状の部分にマンドレル２３０を挿入する。そして、その状態から、メインラム２２０でダイス２５０に向けて押出加工を行うことで、押出製品が得られる。 When performing hot extrusion, first, a hollow billet (not shown) is heated, then inserted into the container 240, and the mandrel 230 is inserted into the hollow portion of the billet. Then, from that state, an extruded product can be obtained by extruding the main ram 220 toward the die 250.

メインシリンダー２１０は、図示しないアキュームレーターにより増圧された高圧流体がメインバルブ２６０を介して供給されることで、メインラム２２０を図中の左方に移動させる。 The main cylinder 210 moves the main ram 220 to the left in the drawing by supplying a high-pressure fluid boosted by an accumulator (not shown) via the main valve 260.

また、メインシリンダー２１０およびピアサシリンダー２１５には、それぞれメイン圧力検出器２７０およびピアサ圧力検出器２７５が接続されている。メインシリンダー２１０およびピアサシリンダー２１５の内部の圧力を検出することによって、熱間押出加工時のビレットに付与される押出力を測定し、その値を押出製品内部温度分布推定装置１０に入力する。 Further, a main pressure detector 270 and a piercer pressure detector 275 are connected to the main cylinder 210 and the piercer cylinder 215, respectively. By detecting the pressure inside the main cylinder 210 and the piercer cylinder 215, the push force applied to the billet during hot extrusion is measured, and the value is input to the extrusion product internal temperature distribution estimation device 10.

また、メインラム２２０には、速度検出器２８０が接続されている。メインラム２２０の移動速度を検出することによって、熱間押出加工時のビレットの押出速度を測定し、その値を押出製品内部温度分布推定装置１０に入力する。 Further, a speed detector 280 is connected to the main ram 220. By detecting the moving speed of the main ram 220, the extrusion speed of the billet during hot extrusion is measured, and the value is input to the extrusion product internal temperature distribution estimation device 10.

図１に示すように、押出製品内部温度分布推定装置１０は、設定部１と、計算部２と、抽出部３と、決定部４とを備える。 As shown in FIG. 1, the extrusion product internal temperature distribution estimation device 10 includes a setting unit 1, a calculation unit 2, an extraction unit 3, and a determination unit 4.

設定部１は、押出製品の長手方向の温度分布を複数（ｎｐ個）設定する。そして、計算部２は、設定部１によって設定された温度分布ならびに押出速度変化目標値および押出力変化目標値を含む情報に基づいて、ビレットの押出速度推定値および押出力推定値の時系列データを算出する。 The setting unit 1 sets a plurality (np) of temperature distributions in the longitudinal direction of the extruded product. Then, the calculation unit 2 performs time-series data of the extrusion speed estimation value and the extrusion speed estimation value of the billet based on the information including the temperature distribution set by the setting unit 1 and the extrusion speed change target value and the extrusion speed change target value. Is calculated.

なお、押出速度変化目標値および押出力変化目標値は、予め定められた押出開始から押出終了までの押出速度変化および押出力変化の目標値のことであり、時系列データとして表わされる。また、時系列データとは、押出開始から押出終了までの所定時間ごとの一連のデータを意味する。 The extrusion speed change target value and the extrusion output change target value are predetermined target values of the extrusion speed change and the extrusion output change from the start of extrusion to the end of extrusion, and are represented as time series data. Further, the time series data means a series of data at predetermined time intervals from the start of extrusion to the end of extrusion.

上記の押出速度推定値および押出力推定値の時系列データは、ｎｐ個の温度分布ごとに算出される。すなわち、ｎｐ回の計算が行われ、押出速度推定値の時系列データおよび押出力推定値の時系列データのそれぞれが、ｎｐ個ずつ得られることとなる。 The time-series data of the extrusion speed estimation value and the extrusion output estimation value are calculated for each of the np temperature distributions. That is, the calculation is performed np times, and np pieces of each of the time-series data of the extrusion speed estimated value and the time-series data of the push output estimated value are obtained.

押出速度推定値および押出力推定値の時系列データの計算方法については特に制限はなく、例えば、図３に構成を示すようなシミュレーションモデルにより計算することができる。また、シミュレーションには、例えば、ＭａｔｈＷｏｒｋｓ社製ＭＡＴＬＡＢ／Ｓｉｍｕｌｉｎｋを使用することができる。 There is no particular limitation on the method of calculating the time series data of the extrusion speed estimated value and the push output estimated value, and the calculation can be performed by, for example, a simulation model as shown in FIG. Further, for the simulation, for example, MATLAB / Simulink manufactured by MathWorks can be used.

抽出部３は、計算部２によって得られた、ｎｐ個の押出速度推定値およびｎｐ個の押出力推定値の時系列データを、押出速度実測値および押出力実測値の時系列データと比較する。そして、押出速度推定値および押出力推定値の最適な時系列データを選択し、それに対応する１つの温度分布を抽出する。 The extraction unit 3 compares the time-series data of the np extrusion speed estimation value and the np push output estimated value obtained by the calculation unit 2 with the time-series data of the extrusion speed measurement value and the push output measurement value. .. Then, the optimum time series data of the extrusion speed estimated value and the extrusion output estimated value are selected, and one temperature distribution corresponding thereto is extracted.

なお、押出速度実測値および押出力実測値の時系列データとしては、上述のメイン圧力検出器２７０、ピアサ圧力検出器２７５および速度検出器２８０から入力されるデータを用いることができる。 As the time series data of the extrusion speed actual measurement value and the push output actual measurement value, the data input from the above-mentioned main pressure detector 270, piercer pressure detector 275 and speed detector 280 can be used.

押出速度推定値および押出力推定値の最適な時系列データを選択する方法については、特に制限はない。例えば、実測値と推定値とで、押出速度の最大値および押出力の最大値、ならびにそれらの最大値に対応する時間から選択される１種以上の指標について比較を行い、平均２乗誤差が最小となる押出速度推定値および押出力推定値の時系列データを選択することができる。上記の指標のうち複数を採用する場合においては、重要な指標に適宜重み付けして比較を行ってもよい。 There is no particular limitation on the method of selecting the optimum time series data of the extrusion speed estimated value and the extrusion output estimated value. For example, the measured value and the estimated value are compared with each other for one or more indexes selected from the maximum value of extrusion speed and the maximum value of extrusion output, and the time corresponding to those maximum values, and the average square error is obtained. Time series data of the minimum extrusion speed estimate and push output estimate can be selected. When a plurality of the above indicators are adopted, important indicators may be appropriately weighted for comparison.

決定部４は、抽出部３によって抽出された１つの温度分布を、押出製品の長手方向における内部の温度分布に決定する。 The determination unit 4 determines one temperature distribution extracted by the extraction unit 3 as the internal temperature distribution in the longitudinal direction of the extruded product.

設定部１は、温度分布の設定を複数回（ｋ_ｍａｘ回）行うことが好ましい。この際には、抽出部３によって抽出された１つの温度分布に基づいて、２回目以降の温度分布の設定を行うことが好ましい。これにより、推定精度が向上し、より真値に近い推定値を得ることが可能になる。 It is preferable that the setting unit 1 sets the temperature distribution a plurality of times (k _{max times).} At this time, it is preferable to set the temperature distribution from the second time onward based on one temperature distribution extracted by the extraction unit 3. As a result, the estimation accuracy is improved, and it becomes possible to obtain an estimated value closer to the true value.

また、温度分布の設定をｋ_ｍａｘ回行う場合においては、２回目以降の温度分布の設定は、粒子群最適化手法（ＰａｒｔｉｃｌｅＳｗａｒｍＯｐｔｉｍｉｚａｔｉｏｎ：ＰＳＯ）を用いて行うことが好ましい。ＰＳＯとは、多目的最適化手法の１つである。 Further, when the temperature distribution is set k _max times, it is preferable to set the temperature distribution from the second time onward by using a particle swarm optimization (PSO). PSO is one of the multi-objective optimization methods.

なお、１回目の温度分布の設定は任意の方法で行えばよく、例えば、乱数を採用すればよい。 The first temperature distribution may be set by any method, for example, a random number may be adopted.

図４は、本発明の一実施形態に係るビレット加熱制御装置の概略構成を示す図である。図４に示すように、ビレット加熱制御装置１００は、押出製品内部温度分布推定装置１０と、温度制御部２０とを備える。 FIG. 4 is a diagram showing a schematic configuration of a billet heating control device according to an embodiment of the present invention. As shown in FIG. 4, the billet heating control device 100 includes an extruded product internal temperature distribution estimation device 10 and a temperature control unit 20.

また、図５は、ビレット加熱制御装置１００を備えた押出機３００の概略構成を示す図である。押出機３００の構成は、図２に示す押出機２００の構成と同様であるが、ビレットヒータ２４５をさらに備えている。ビレットヒータ２４５は、温度制御部２０からの指示に応じてビレットの加熱温度を調整する。ビレットヒータ２４５は、ビレット長手方向の加熱温度を調節するいわゆる「傾斜加熱」が可能な機能を有することが好ましい。 Further, FIG. 5 is a diagram showing a schematic configuration of an extruder 300 provided with a billet heating control device 100. The configuration of the extruder 300 is the same as that of the extruder 200 shown in FIG. 2, but further includes a billet heater 245. The billet heater 245 adjusts the heating temperature of the billet according to an instruction from the temperature control unit 20. The billet heater 245 preferably has a function capable of so-called "tilt heating" for adjusting the heating temperature in the longitudinal direction of the billet.

温度制御部２０は、押出製品内部温度分布推定装置１０が推定した押出製品の温度分布、ならびに上述の押出速度変化目標値および押出力変化目標値を含む情報を用いて算出される、押出速度推定値および押出力推定値の時系列データに基づき、ビレットの加熱温度の制御を行う。 The temperature control unit 20 estimates the extrusion speed, which is calculated by using the temperature distribution of the extrusion product estimated by the extrusion product internal temperature distribution estimation device 10 and the information including the above-mentioned extrusion speed change target value and push output change target value. The billet heating temperature is controlled based on the time-series data of the value and the extrusion output estimated value.

具体的には、押出速度推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を下げ、押出力推定値の最大値が、予め設定される値を超える場合には、ビレットの加熱温度を上げる制御を行う。 Specifically, when the maximum value of the extrusion speed estimation value exceeds the preset value, the heating temperature of the billet is lowered, and when the maximum value of the extrusion speed estimation value exceeds the preset value. Controls to raise the heating temperature of the billet.

次に、本発明の一実施形態に係る押出製品内部温度分布推定方法を、図６を用いて説明する。図６は、本発明の一実施形態に係る押出製品内部温度分布推定方法を示すフロー図である。なお、以降の説明においては、ＰＳＯを採用する場合を例に挙げるが、本発明はこれに限定されるものではない。 Next, the method for estimating the internal temperature distribution of the extruded product according to the embodiment of the present invention will be described with reference to FIG. FIG. 6 is a flow chart showing a method for estimating the internal temperature distribution of an extruded product according to an embodiment of the present invention. In the following description, the case where PSO is adopted will be given as an example, but the present invention is not limited to this.

図６に示すように、設定部１が乱数にて、ｎｐ個の温度分布ｘ_ｉ（ｉ＝１〜ｎｐ）を設定する（ステップＡ１）。温度分布ｘ_ｉを設定するに際しては、押出製品の先端からの距離に応じた温度を１つの粒子がもつ座標として定義する。 As shown in FIG. 6, the setting unit 1 sets np temperature distributions x _i (i = 1 to np) with random numbers (step A1). When setting the temperature distribution x _i , the temperature corresponding to the distance from the tip of the extruded product is defined as the coordinates of one particle.

次に、計算部２は、設定部１によって設定されたｋ回目の温度分布ｘ_ｉ ^ｋならびに押出速度変化目標値および押出力変化目標値を含む情報に基づいて、ｎｐ回のシミュレーションにより、ｎｐ個のビレットの押出速度推定値および押出力推定値の時系列データを算出する（ステップＡ２）。 Next, the calculation unit 2 performs ^{np simulations based on the information including the kth} _{temperature distribution x ik} set by the setting unit 1 and the extrusion speed change target value and the push output change target value. The time series data of the extrusion speed estimated value and the push output estimated value of the billet of No. 1 are calculated (step A2).

次に、抽出部３は、ｎｐ個の押出速度推定値およびｎｐ個の押出力推定値の時系列データと、押出速度実測値および押出力実測値の時系列データとを、評価関数ｆ（ｘ_ｉ ^ｋ）を用いて比較する（ステップＡ３）。評価関数ｆ（ｘ_ｉ ^ｋ）としては、実測値と推定値との、押出速度の最大値および押出力の最大値、ならびにそれらの最大値に対応する時間から選択される１種以上の平均２乗誤差を用いることができる。 Next, the extraction unit 3 uses the evaluation function f (x) to evaluate the time-series data of the np extrusion speed estimated value and the mp push-output estimated value and the time-series data of the extrusion speed actual measurement value and the push-output actual measurement value. _i ^k) compared with (step A3). The evaluation function f (x i _^k), between the actual measurement value and the estimated value, the maximum value of the maximum value and the extrusion force of the extrusion rate, and the average of one or more selected from the time corresponding to the maximum value thereof 2 The power error can be used.

次に、抽出部３は、ｎｐ回のシミュレーション結果のうち、評価関数ｆ（ｘ_ｉ ^ｋ）が最小となった粒子の座標から、個の最適点ｘ_ｉ ^{ｂｅｓｔ，ｋ}および群の最適点ｘ_{ｓｗａｒｍ} ^{ｂｅｓｔ，ｋ}を抽出する（ステップＡ４）。 Next, the extraction unit 3 extracts _{the optimum points x i} ^{best, k} and the optimum points x _{swarm of the} individual points x i best, k and the group from the coordinates of the particles having the smallest _{evaluation function f (x i} ^{k) in the simulation results of np times.} ^{The best and k} are extracted (step A4).

温度分布の設定がｋ_ｍａｘ回目以下である場合（ステップＡ５でＮｏの場合）、設定部は抽出された個の最適点ｘ_ｉ ^{ｂｅｓｔ，ｋ}および群の最適点ｘ_{ｓｗａｒｍ} ^{ｂｅｓｔ，ｋ}に近づくようにｎｐ個の温度分布を更新し（ステップＡ６）、再度、ｎｐ個のビレットの押出速度推定値および押出力推定値の時系列データを算出する（ステップＡ２）。 When the temperature distribution is set to the k _maxth time or less (No in step A5), the setting unit should approach the _{extracted individual optimum points x i} ^{best, k} and the group optimum points x _swarm ^{best, k.} The temperature distribution of the np pieces is updated (step A6), and the time series data of the extrusion rate estimated value and the push output estimated value of the np number billets are calculated again (step A2).

温度分布の設定がｋ_ｍａｘ回目を超える場合（ステップＡ５でＹｅｓの場合）、決定部４は、温度分布を決定する（ステップＡ７）。このとき、ｋ_ｍａｘ回目の群の最適点を温度分布として採用する。 When the temperature distribution setting _{exceeds the kmaxth} time (Yes in step A5), the determination unit 4 determines the temperature distribution (step A7). At this time, _{the optimum point of the kmax} th group is adopted as the temperature distribution.

次に、本発明の一実施形態に係るビレット加熱制御方法を、図７を用いて説明する。図７は、本発明の一実施形態に係るビレット加熱制御方法を示すフロー図である。 Next, the billet heating control method according to the embodiment of the present invention will be described with reference to FIG. 7. FIG. 7 is a flow chart showing a billet heating control method according to an embodiment of the present invention.

図７に示すように、押出製品内部温度分布推定装置１０は、推定された押出製品の温度分布ならびに押出速度変化目標値および押出力変化目標値を含む情報を用いて押出速度推定値および押出力推定値の時系列データを算出する（ステップＢ１）。 As shown in FIG. 7, the extrusion product internal temperature distribution estimation device 10 uses the estimated temperature distribution of the extruded product and the information including the extrusion rate change target value and the extrusion rate change target value to estimate the extrusion rate and the extrusion output. Calculate the time series data of the estimated value (step B1).

次に、温度制御部２０は、押出速度推定値の最大値ＯＳ_ｓｉｍおよび押出力推定値の最大値Ｐｔｐ_ｓｉｍを求める（ステップＢ２）。 Next, the temperature control unit 20 obtains the maximum value _{Ptp sim} maximum value _{OS sim} and extrusion force estimated value of the extrusion speed estimation value (step B2).

そして、温度制御部２０は、押出速度推定値の最大値ＯＳ_ｓｉｍが、予め設定される値ＯＳ_ｍａｘを超える場合（ステップＢ３でＮｏの場合）には、ビレットの加熱温度を下げる制御を行う（ステップＢ４）。また、押出力推定値の最大値Ｐｔｐ_ｓｉｍが、予め設定される値Ｐｔｐ_ｍａｘを超える場合（ステップＢ５でＮｏの場合）には、ビレットの加熱温度を上げる制御を行う（ステップＢ６）。 Then, the temperature control unit 20 controls to lower the heating temperature of the billet when _{the maximum value OS sim} of the extrusion speed estimated value exceeds the preset value OS _{max (when No in step B3).} Step B4). Further, when the maximum value Ptp _sim of the push output estimated value exceeds the preset value Ptp _max (No in step B5), control is performed to raise the heating temperature of the billet (step B6).

押出速度推定値の最大値ＯＳ_ｓｉｍが、予め設定される値ＯＳ_ｍａｘ以下で（ステップＢ３でＹｅｓの場合）、かつ押出力推定値の最大値Ｐｔｐ_ｓｉｍが、予め設定される値Ｐｔｐ_ｍａｘ以下である場合（ステップＢ５でＹｅｓの場合）には、温度制御部２０は加熱温度を維持する（ステップＢ７）。 _{The maximum value OS sim} of the extrusion speed estimation value is equal to _{or less than the preset value OS max} (in the case of Yes in step B3), and the maximum value Ptp _sim of the push output estimated value is equal to or less than the _{preset value Ptp max.} In some cases (yes in step B5), the temperature control unit 20 maintains the heating temperature (step B7).

本発明の一実施形態に係るプログラムは、コンピュータに、図６に示すステップＡ１〜Ａ７を実行させるプログラムであればよい。このプログラムをコンピュータにインストールし、実行することによって、本実施の形態における押出製品内部温度分布推定装置１０を実現することができる。この場合、コンピュータのＣＰＵ（Central Processing Unit）は、設定部１、計算部２、抽出部３および決定部４として機能し、処理を行う。 The program according to the embodiment of the present invention may be any program that causes a computer to execute steps A1 to A7 shown in FIG. By installing this program on a computer and executing it, the extrusion product internal temperature distribution estimation device 10 according to the present embodiment can be realized. In this case, the CPU (Central Processing Unit) of the computer functions as a setting unit 1, a calculation unit 2, an extraction unit 3, and a determination unit 4 to perform processing.

また、本発明の他の実施形態に係るプログラムは、コンピュータに、図７に示すステップＢ１〜Ｂ７を実行させるプログラムであればよい。このプログラムをコンピュータにインストールし、実行することによって、本実施の形態におけるビレット加熱制御装置１００を実現することができる。この場合、コンピュータのＣＰＵ（Central Processing Unit）は、押出製品内部温度分布推定装置１０および温度制御部２０として機能し、処理を行う。 Further, the program according to another embodiment of the present invention may be a program that causes a computer to execute steps B1 to B7 shown in FIG. By installing this program on a computer and executing it, the billet heating control device 100 according to the present embodiment can be realized. In this case, the CPU (Central Processing Unit) of the computer functions as the extrusion product internal temperature distribution estimation device 10 and the temperature control unit 20 to perform processing.

１設定部
２計算部
３抽出部
４決定部
１０押出製品内部温度分布推定装置
２０温度制御部
１００ビレット加熱制御装置
２００押出機
２１０メインシリンダー
２１５ピアサシリンダー
２２０メインラム
２２５ピアサラム
２３０マンドレル
２４０コンテナ
２４５ビレットヒータ
２５０ダイス
２５５ダイスホルダー
２６０メインバルブ
２６５ピアサバルブ
２７０メイン圧力検出器
２７５ピアサ圧力検出器
２８０速度検出器
３００押出機

1 Setting unit 2 Calculation unit 3 Extraction unit 4 Decision unit 10 Extruded product internal temperature distribution estimation device 20 Temperature control unit 100 Billet heating control device 200 Extruder 210 Main cylinder 215 Piercer cylinder 220 Main ram 225 Piersaram 230 Mandrel 240 Container 245 Billet Heater 250 Die 255 Die Holder 260 Main Valve 265 Piercer Valve 270 Main Pressure Detector 275 Piercer Pressure Detector 280 Speed Detector 300 Extruder

Claims

A device that estimates the internal temperature distribution in the longitudinal direction of an extruded product when a billet is hot-extruded under predetermined extrusion speed and extrusion power conditions to manufacture an extruded product.
A setting unit that sets a plurality of temperature distributions in the longitudinal direction of the extruded product, and
For each of the plurality of temperature distributions set by the setting unit, a plurality of extrusion speed estimation values and a plurality of extrusion speed estimation values are based on the information including the temperature distribution and predetermined extrusion speed change target values and extrusion output change target values. A calculation unit that calculates time-series data of extrusion output estimates,
The time series data of the plurality of extrusion speed estimates and the plurality of push output estimates obtained by the calculation unit are compared with the time series data of the extrusion speed measurement value and the push output measurement value, and the maximum extrusion speed is maximized. maximum values and extrusion force, and for to compare one or more indicators being selected from the time corresponding to the maximum value thereof, when the mean square error extrusion speed estimate and extrusion force estimate becomes the minimum An extraction unit that extracts one temperature distribution corresponding to series data,
A determination unit is provided, which determines the temperature distribution of the one extracted by the extraction unit as the internal temperature distribution in the longitudinal direction of the extruded product.
Extruded product internal temperature distribution estimation device during hot extrusion processing.

The setting unit sets the plurality of temperature distributions a plurality of times, and sets the temperature distribution from the second time onward based on the one temperature distribution extracted by the extraction unit.
The device for estimating the internal temperature distribution of an extruded product during hot extrusion according to claim 1.

The setting unit sets the temperature distribution from the second time onward by using a particle swarm optimization method.
The device for estimating the internal temperature distribution of an extruded product during hot extrusion according to claim 2.

This is a method of estimating the internal temperature distribution in the longitudinal direction of the extruded product when the billet is hot-extruded under predetermined extrusion speed and extrusion power conditions to manufacture the extruded product.
(A) A step of setting a plurality of temperature distributions in the longitudinal direction of the extruded product, and
(B) For each of the plurality of temperature distributions set in the step (a), a plurality of extrusions are performed based on the temperature distribution and information including a predetermined extrusion speed change target value and push output change target value. Steps to calculate time series data for velocity estimates and multiple extrusion estimates,
(C) Compare the time series data of the plurality of extrusion speed estimated values and the plurality of push output estimated values obtained in the step (b) with the time series data of the extrusion speed measured value and the push output measured value. Then, the maximum value of the extrusion speed and the maximum value of the push output, and one or more indexes selected from the time corresponding to those maximum values are compared, and the estimated value of the extrusion speed and the estimated value of the extrusion speed that minimizes the average square error and the mean square error are compared. A step to extract one temperature distribution corresponding to the time series data of the push output estimated value, and
(D) The step comprises determining the temperature distribution of the one extracted in the step (c) to the internal temperature distribution in the longitudinal direction of the extruded product.
A method for estimating the internal temperature distribution of an extruded product during hot extrusion.

In the step (a), the plurality of temperature distributions are set a plurality of times, and the second and subsequent temperature distributions are set based on the one temperature distribution extracted by the extraction unit.
The method for estimating the internal temperature distribution of an extruded product during hot extrusion according to claim 4.

In the step (a), the temperature distribution is set from the second time onward by using the particle swarm optimization method.
The method for estimating the internal temperature distribution of an extruded product during hot extrusion according to claim 5.

A device that controls the heating temperature of billets during hot extrusion.
The extruded product internal temperature distribution estimation device according to any one of claims 1 to 3.
Equipped with a temperature control unit
The temperature control unit calculates the extrusion speed using information including the temperature distribution of the extrusion product estimated by the extrusion product internal temperature distribution estimation device, the extrusion speed change target value, and the extrusion output change target value. Based on time series data of estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
Billet heating control device.

A method of controlling the heating temperature of billets during hot extrusion.
The steps (a) to (d) according to any one of claims 4 to 6 and
(E) The step of changing the heating temperature of the billet is provided.
In step (e), extrusion calculated using information including the temperature distribution of the extruded product determined in step (d), the extrusion speed change target value, and the extrusion output change target value. Based on time series data of velocity estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
Billet heating control method.

A program for estimating the internal temperature distribution in the longitudinal direction of the extruded product when the billet is hot-extruded under predetermined extrusion speed and extrusion power conditions by a computer to manufacture the extruded product. There,
On the computer
(A) A step of setting a plurality of temperature distributions in the longitudinal direction of the extruded product, and
(B) For each of the plurality of temperature distributions set in the step (a), a plurality of extrusions are performed based on the temperature distribution and information including a predetermined extrusion speed change target value and push output change target value. Steps to calculate time series data for velocity estimates and multiple extrusion estimates,
(C) Compare the time series data of the plurality of extrusion speed estimated values and the plurality of push output estimated values obtained in the step (b) with the time series data of the extrusion speed measured value and the push output measured value. Then, the maximum value of the extrusion speed and the maximum value of the push output, and one or more indexes selected from the time corresponding to those maximum values are compared, and the estimated value of the extrusion speed and the estimated value of the extrusion speed that minimizes the average square error and the mean square error are compared. A step to extract one temperature distribution corresponding to the time series data of the push output estimated value, and
(D) The step of determining the temperature distribution inside the extruded product in the longitudinal direction of the one temperature distribution extracted in the step (c) is executed.
program.

In the step (a), the plurality of temperature distributions are set a plurality of times, and the second and subsequent temperature distributions are set based on the one temperature distribution extracted by the extraction unit.
The program according to claim 9.

In the step (a), the temperature distribution is set from the second time onward by using the particle swarm optimization method.
The program according to claim 10.

A program for controlling the heating temperature of billets during hot extrusion by a computer.
On the computer
The steps (a) to (d) according to any one of claims 9 to 11.
(E) Perform the steps of changing the heating temperature of the billet,
In step (e), extrusion calculated using information including the temperature distribution of the extruded product determined in step (d), the extrusion speed change target value, and the extrusion output change target value. Based on time series data of velocity estimates and extrusion estimates
When the maximum value of the extrusion speed estimation value exceeds a preset value, the heating temperature of the billet is lowered.
When the maximum value of the push output estimated value exceeds a preset value, the billet heating temperature is raised.
program.