WO2018158828A1 - シート材生産ラインの数学モデル算出装置および制御装置 - Google Patents

シート材生産ラインの数学モデル算出装置および制御装置 Download PDF

Info

Publication number
WO2018158828A1
WO2018158828A1 PCT/JP2017/007823 JP2017007823W WO2018158828A1 WO 2018158828 A1 WO2018158828 A1 WO 2018158828A1 JP 2017007823 W JP2017007823 W JP 2017007823W WO 2018158828 A1 WO2018158828 A1 WO 2018158828A1
Authority
WO
WIPO (PCT)
Prior art keywords
mathematical model
motor
roll
sheet material
production line
Prior art date
Application number
PCT/JP2017/007823
Other languages
English (en)
French (fr)
Inventor
鈴木　敦
Original Assignee
東芝三菱電機産業システム株式会社
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 東芝三菱電機産業システム株式会社 filed Critical 東芝三菱電機産業システム株式会社
Priority to KR1020197013943A priority Critical patent/KR102186020B1/ko
Priority to CN201780063842.3A priority patent/CN109843764B/zh
Priority to PCT/JP2017/007823 priority patent/WO2018158828A1/ja
Priority to JP2019502326A priority patent/JP6835200B2/ja
Publication of WO2018158828A1 publication Critical patent/WO2018158828A1/ja

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
    • B65H—HANDLING THIN OR FILAMENTARY MATERIAL, e.g. SHEETS, WEBS, CABLES
    • B65H23/00—Registering, tensioning, smoothing or guiding webs
    • B65H23/04—Registering, tensioning, smoothing or guiding webs longitudinally
    • B65H23/18—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web
    • B65H23/182—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web in unwinding mechanisms or in connection with unwinding operations
    • B65H23/185—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web in unwinding mechanisms or in connection with unwinding operations motor-controlled
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
    • B65H—HANDLING THIN OR FILAMENTARY MATERIAL, e.g. SHEETS, WEBS, CABLES
    • B65H23/00—Registering, tensioning, smoothing or guiding webs
    • B65H23/04—Registering, tensioning, smoothing or guiding webs longitudinally
    • B65H23/18—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web
    • B65H23/195—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web in winding mechanisms or in connection with winding operations
    • B65H23/198—Registering, tensioning, smoothing or guiding webs longitudinally by controlling or regulating the web-advancing mechanism, e.g. mechanism acting on the running web in winding mechanisms or in connection with winding operations motor-controlled (Controlling electrical drive motors therefor)
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
    • B65H—HANDLING THIN OR FILAMENTARY MATERIAL, e.g. SHEETS, WEBS, CABLES
    • B65H2557/00—Means for control not provided for in groups B65H2551/00 - B65H2555/00
    • B65H2557/20—Calculating means; Controlling methods
    • B65H2557/24—Calculating methods; Mathematic models
    • B65H2557/242—Calculating methods; Mathematic models involving a particular data profile or curve

Definitions

  • the mathematical model calculation device for a sheet material production line is arranged between a plurality of roll devices on the basis of a history of actual measurement values of the sheet material production line and directed from one of the plurality of roll devices to the other.
  • the torque of the second motor device that drives the other of the plurality of roll devices and the plurality of roll devices
  • a mathematical model calculation unit is provided that calculates a mathematical model using a deviation from the torque of the first motor device that drives one as an input and using a measurement value of the tensiometer as an output.
  • the mathematical model calculation device for a sheet material production line is arranged between a plurality of roll devices on the basis of a history of actual measurement values of the sheet material production line and directed from one of the plurality of roll devices to the other.
  • the plurality of roll devices at a rotational angular velocity of a second motor device that drives the other of the plurality of roll devices with respect to a sheet material production line controlled based on a measurement value of a tensiometer that measures the tension of the traveling sheet material
  • the tensiometer is inputted with a deviation between a value obtained by multiplying the other radius of the first roll and the rotation angular velocity of the first motor device that drives one of the plurality of roll devices by one radius of the plurality of roll devices.
  • a mathematical model calculation unit that calculates a mathematical model using the measured value of as an output.
  • the mathematical model calculation device receives the torque of the first motor device and the torque of the second motor device as inputs while the sheet material production line is in operation.
  • the control part which acquires the tension
  • the mathematical model is calculated based on the history of actual measurement values of the sheet material production line. For this reason, it is possible to calculate a highly accurate mathematical model for controlling the tension of the sheet material.
  • FIG. 1 is a configuration diagram of a winder slit facility to which a mathematical model calculation apparatus according to Embodiment 1 of the present invention is applied.
  • the sheet material production line in Fig. 1 is a winder slit facility.
  • the unwinder 1 is provided upstream of the winder slit facility as one of the roll devices.
  • the winder 2 is provided on the downstream side of the winder slit facility as the other of the roll devices.
  • the winder 2 includes a front drum 3, a rear drum 4, and a rider roll 5.
  • the slitter 6 is provided between the unwinder 1 and the winder 2.
  • the intermediate roll 7 is provided between the unwinder 1 and the slitter 6.
  • the control device 21 includes a mathematical model calculation device 22.
  • the mathematical model calculation device 22 includes a mathematical model calculation unit 23.
  • the mathematical model calculation unit 23 is provided so that a mathematical model of the winder slit facility can be calculated based on a history of actual measurement values of the winder slit facility.
  • the control device 21 includes a control unit 24.
  • the control unit 24 is provided so as to control the operation of the winder slit facility based on the mathematical model calculated by the mathematical model calculation unit 23.
  • the paper material 25 is disposed between the unwinder 1 and the winder 2.
  • the paper material 25 advances from the unwinder 1 toward the winder 2.
  • the paper material 25 is wound around the unwinder 1 in advance as a sheet material.
  • the winding diameter of the paper material 25 is large.
  • the paper material 25 is heavy.
  • the paper material 25 is cut into a preset width by the slitter 6 via the intermediate roll 7.
  • the paper material 25 is wound around the winder 2.
  • the paper material 25 is wound up to a preset outer diameter. The winding diameter of the paper material 25 at this time is smaller than the winding diameter of the paper material 25 in the unwinder 1.
  • FIG. 2 is a diagram for explaining the control of the winder slit facility to which the mathematical model calculation apparatus according to Embodiment 1 of the present invention is applied.
  • the input unit of the control device 21 includes an output unit of the tension meter 20, an output unit of the front drum drive unit 15, an output unit of the rear drum drive unit 16, and an output of the rider roll drive unit 17. Connected to the part.
  • the output unit of the control device 21 includes an input unit of the unwinder drive unit 14, an input unit of the front drum drive unit 15, an input unit of the rear drum drive unit 16, an input unit of the rider roll drive unit 17, and a slitter unit.
  • the input unit of the drive unit 18 and the input unit of the intermediate roll drive unit 19 are connected.
  • speed control is performed. Specifically, the control unit 24 outputs the rotational angular velocity reference value ⁇ fd ref (rad / s) to the front drum drive device 15.
  • the front drum drive device 15 drives the front drum motor 9 based on the rotational angular velocity reference value ⁇ fd ref .
  • the front drum motor 9 drives the front drum 3.
  • speed control is performed. Specifically, the control unit 24 outputs the rotational angular velocity reference value ⁇ rr ref (rad / s) to the rider roll drive device 17.
  • the rider roll drive device 17 drives the rider roll motor 11 based on the rotational angular velocity reference value ⁇ rr ref .
  • the rider roll motor 11 drives the rider roll 5.
  • speed control is performed. Specifically, the control unit 24 outputs the rotation angular velocity reference value ⁇ sl ref (rad / s) to the slitter drive device 18.
  • the slitter drive device 18 drives the slitter motor 12 based on the rotational angular velocity reference value ⁇ sl ref .
  • the slitter motor 12 drives the slitter 6.
  • speed control is performed. Specifically, the control unit 24 outputs the rotation angular velocity reference value ⁇ int ref (rad / s) to the intermediate roll drive device 19.
  • the intermediate roll drive device 19 drives the intermediate roll motor 13 based on the rotation angular velocity reference value ⁇ int ref .
  • the intermediate roll motor 13 drives the intermediate roll 7.
  • the control unit 24 receives an input of the torque reference value ⁇ fd ref (N ⁇ m) from the front drum drive device 15.
  • the control unit 24 receives an input of the torque reference value ⁇ rd ref (N ⁇ m) from the rear drum drive device 16.
  • the control unit 24 receives an input of the torque reference value ⁇ rr ref (N ⁇ m) from the rider roll drive device 17.
  • the control unit 24 receives an input of a tension response value T res (MPa) that is a measurement value of the tension meter 20.
  • tension control is performed. Specifically, the control unit 24 applies the torque reference to the unwinding drive device 14 based on the torque reference value ⁇ fd ref , the torque reference value ⁇ rd ref , the torque reference value ⁇ rr ref, and the tension response value T res.
  • the value ⁇ uw ref (N ⁇ m) is output.
  • the unwinder drive device 14 drives the unwinder motor 8 based on the torque reference value ⁇ uw ref .
  • the unwinder motor 8 drives the unwinder 1.
  • the paper material 25 receives a certain tension when being wound up by the winder 2.
  • the tension applied to the paper material 25 is constant, wrinkles of the paper material 25 and cutting of the paper material 25 are suppressed.
  • the operation of the winder slit facility is stabilized.
  • the winder slit equipment is stably operated, product quality is improved.
  • FIG. 3 is a block diagram for explaining tension control of a paper material using a mathematical model calculated by the mathematical model calculation apparatus according to Embodiment 1 of the present invention.
  • control unit 24 includes a tension predictor 26 and a PI controller 27.
  • the tension predictor 26 uses the transfer function G T (s) corresponding to the mathematical model, and the sum of the torque reference value ⁇ fd ref , the torque reference value ⁇ rd ref, and the torque reference value ⁇ rr ref and the torque reference value ⁇ uw. calculating the tension predicted value T m (MPa) from the deviation between ref.
  • PI controller 27 calculates a proportional gain K p and the integral gain K I and the Laplace operator s torque reference value from the deviation between the tension reference value T ref (MPa) and tensile predicted value T m based on the tau uw ref To do.
  • FIGS. 4 to 6 are block diagrams for explaining an outline of a mathematical model calculation method performed by the mathematical model calculation apparatus according to Embodiment 1 of the present invention.
  • the mathematical model calculation unit 23 creates a mathematical model for estimating the tension applied to the paper material 25 from experimentally obtained data. Specifically, the mathematical model calculation unit 23 calculates the torque response value ⁇ fd res (N ⁇ m) of the front drum motor 9, the torque response value ⁇ rd res (N ⁇ m) of the rear drum motor 10, and the rider roll the difference between the torque response value of the motor 11 ⁇ rr res (N ⁇ m ) torque response value of the sum and the unwinder motor 8 and ⁇ uw res (N ⁇ m) as an input, the tension meter 20 of the tension response value T res Is used as an output to calculate a differential equation, transfer function, and state equation from input to output.
  • a (z) in the expression (1) is expressed by the following expression (2).
  • G (z) is a transfer function from the input u (k) to the output y (k). Specifically, G (z) is expressed by the following equation (5).
  • the predicted value of the output y (k) at the current time k is expressed by the following equation (8) using past data up to the time (k ⁇ 1).
  • Equation (8) The second term on the right side of equation (8) is defined by the following equation (9).
  • the current output is calculated as a linear combination of past input and output.
  • the prediction error ⁇ using the one-step predicted value is defined by the following equation (12).
  • a (z), B (z), and C (z) are determined by the prediction error method using the equation (12). Specifically, A (z), B (z), and C (z) are determined so as to minimize the evaluation function composed of the prediction error ⁇ .
  • the transfer function G ′ (s) from the input u (k) to the output y (k) can be obtained by converting the equation (13).
  • the transfer function G ′ (s) is expressed by the following equation (14).
  • FIGS. 7 to 9 are diagrams for explaining an example of a mathematical model calculation method by the mathematical model calculation apparatus according to Embodiment 1 of the present invention.
  • FIG. 7 shows the unwinding motor and the sum of the torque response value ⁇ fd res of the front drum motor 9, the torque response value ⁇ rd res of the rear drum motor 10, and the torque response value ⁇ rr res of the rider roll motor 11.
  • 8 shows a deviation from the torque response value ⁇ uw res of 8. 7
  • the torque response value ⁇ fd res of the front drum motor 9 the torque response value ⁇ rd res of the rear drum motor 10
  • the torque response value ⁇ rr res of the rider roll motor 11 and the torque of the unwinder motor 8.
  • the response value ⁇ uw res is a value converted with the rated torque as 100%.
  • the upper part of FIG. 8 is a diagram showing a transient state from the start in the upper part of FIG. 7 to the target tension. Specifically, the upper part of FIG. 8 is an enlarged view from time 0 (s) to 20 (s). In the upper part of FIG. 8, a dotted line is a line obtained by linearly approximating data.
  • the lower part of FIG. 8 is a diagram showing a transitional state from the start of the lower part of FIG. 7 to the target tension. Specifically, the lower part of FIG. 8 is an enlarged view of the time from 0 (s) to 20 (s). In the lower part of FIG. 8, a dotted line is a line obtained by linearly approximating data.
  • the lower part of FIG. 9 is data obtained by removing the average value and the slope, which are low-frequency disturbances, from the lower part of FIG. 8 as preprocessing of data before being input to the mathematical model.
  • the mathematical model calculation unit 23 performs processing corresponding to the upper and lower stages of FIG. 9 offline after data collection.
  • the mathematical model calculation unit 23 calculates a transfer function using the upper and lower data in FIG.
  • the transfer function G T (s) is expressed by the following equation (15).
  • the obtained transfer function G T (s) is used for tension control of the paper material 25 from the next time onward.
  • the mathematical model calculator 23 removes from the input data when estimating the transfer function G T (s) in the previous paper 25 from the deviation data of the winding and rewinding torques calculated from the torque reference value of each motor. Remove the average value and slope online.
  • the mathematical model calculation unit 23 inputs the value to the transfer function G T (s).
  • the mathematical model calculation unit 23 calculates a value obtained by adding the average value and the slope removed from the output data when calculating the transfer function G T (s) to the output of the transfer function G T (s). Calculate as m .
  • the predicted tension value Tm is calculated so as to be along the actual tension response value Tres .
  • the tension of the paper material 25 is controlled with high accuracy.
  • the mathematical model is calculated based on the history of the actual measurement values of the sheet material production line. For this reason, it is possible to calculate a highly accurate mathematical model for controlling the tension of the sheet material.
  • the mathematical model is the sum of the torque response value ⁇ fd res of the front drum motor 9, the torque response value ⁇ rd res of the rear drum motor 10, and the torque response value ⁇ rr res of the rider roll motor 11.
  • the deviation from the torque response value ⁇ uw res of the unwinder motor 8 is used as an input, and the tension response value T res of the tensiometer 20 is calculated as an output. For this reason, a highly accurate mathematical model can be calculated for the winder slit installation.
  • the mathematical model may be calculated after weighting the value ⁇ uw res . In this case, a highly accurate mathematical model can be calculated according to the actual situation.
  • the mathematical model is calculated using a value obtained by approximating an input and an output with a linear line and subtracting the corresponding linear line from the input and the output. For this reason, a more accurate mathematical model can be calculated.
  • the mathematical model is calculated by using a value obtained by applying a low-pass filter to the measured value of the tensiometer 20 as an output. For this reason, the noise contained in the measured value in the tension meter 20 can be removed. As a result, a more accurate mathematical model can be calculated.
  • a mathematical model may be calculated by using an average value of a plurality of sampling values of the measurement value of the tension meter 20 as an output. In this case, noise included in the measurement value obtained by the tensiometer 20 can be removed. As a result, a more accurate mathematical model can be calculated.
  • the value obtained by multiplying the roll radius corresponding to the fastest rotation speed among the front drum motor 9, the rear drum motor 10, and the rider roll motor 11 by the rotation of the unwinder motor 8 is calculated.
  • a deviation from a value obtained by multiplying the speed by the radius of the unwinder 1 may be input.
  • a transfer function from the peripheral speed of the roll to the tension of the paper material 25 may be calculated as a transfer function corresponding to the mathematical model. Also in this case, a highly accurate mathematical model can be calculated.
  • FIG. 11 is a hardware configuration diagram of a control device including the mathematical model calculation device according to Embodiment 1 of the present invention.
  • Each function of the control device 21 can be realized by a processing circuit.
  • the processing circuit comprises at least one processor 28a and at least one memory 28b.
  • the processing circuit comprises at least one dedicated hardware 29.
  • each function of the control device 21 is realized by software, firmware, or a combination of software and firmware.
  • At least one of software and firmware is described as a program.
  • At least one of software and firmware is stored in at least one memory 28b.
  • the at least one processor 28a reads out and executes the program stored in the at least one memory 28b, thereby realizing each function of the control device 21.
  • the at least one processor 28a is also referred to as a CPU (Central Processing Unit), a central processing unit, a processing unit, an arithmetic unit, a microprocessor, a microcomputer, and a DSP.
  • the at least one memory 28b is a nonvolatile or volatile semiconductor memory such as a RAM, a ROM, a flash memory, an EPROM, or an EEPROM, a magnetic disk, a flexible disk, an optical disk, a compact disk, a mini disk, or a DVD.
  • a nonvolatile or volatile semiconductor memory such as a RAM, a ROM, a flash memory, an EPROM, or an EEPROM, a magnetic disk, a flexible disk, an optical disk, a compact disk, a mini disk, or a DVD.
  • the processing circuit comprises at least one dedicated hardware 29, the processing circuit is implemented, for example, as a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an ASIC, an FPGA, or a combination thereof.
  • the control device 21 is realized by a processing circuit.
  • each function of the control device 21 is collectively realized by a processing circuit.
  • a part of each function of the control device 21 may be realized by dedicated hardware 29, and the other part may be realized by software or firmware.
  • the function of the mathematical model calculation unit 23 is realized by a processing circuit as dedicated hardware 29, and at least one processor 28a is stored in at least one memory 28b for functions other than the function of the mathematical model calculation unit 23. It may be realized by reading out and executing the program.
  • the processing circuit realizes each function of the control device 21 by hardware 29, software, firmware, or a combination thereof.
  • the mathematical model calculation unit 23 converts each calculated transfer function into a state equation, and improves the prediction accuracy by using a Kalman filter.
  • the Kalman filter is an online algorithm that sequentially estimates the state of the system using a state equation indicating the dynamic characteristics of the system and measured data that is given every moment.
  • the Kalman filter uses an assumption that the influence of the disturbance applied to the system and the noise contained in the sensor follows a normal distribution.
  • the state at the current time is estimated based on the information of the time before one sampling.
  • the correct state is estimated by correcting the predicted tension value based on the actual measurement value at the current time.
  • the mathematical model calculation unit 23 converts each transfer function calculated in the first embodiment into a state equation of Expression (17), and implements a Kalman filter. Mathematical model calculating unit 23 sequentially corrects the tension predicted value T c m by using a deviation between the tension response value T res and tension predicted value T c m for obtaining at each sampling. Specifically, the mathematical model calculating unit 23 calculates the tension prediction value T c m by performing at each sampling operations (18) from equation (22).
  • the mathematical model calculation unit 23 converts the mathematical model into a state equation, and uses the Kalman filter to calculate the mathematical model based on the measurement value of the tensiometer 20 during the operation of the winder slitter facility. model by correcting the tension predicted value T c m. For this reason, the precision of tension control of paper material 25 can be raised more.
  • control device similar to the control device 21 of the first embodiment may be applied to the film production line.
  • the accuracy of tension control can be increased based on a highly accurate mathematical model.
  • FIG. FIG. 13 is a configuration diagram of rolling equipment to which the mathematical model calculating apparatus according to Embodiment 3 of the present invention is applied.
  • symbol is attached
  • the forward rolling stand 30 is provided on the upstream side of the rolling equipment as one of the roll devices.
  • the front rolling roll 31 is provided on the front rolling stand 30.
  • the forward reduction device 32 is provided above the forward rolling stand 30.
  • the rear rolling stand 33 is provided on the downstream side of the rolling equipment as the other of the roll devices.
  • the rear rolling roll 34 is provided on the rear rolling stand 33.
  • the rear reduction device 35 is provided above the rear rolling stand 33.
  • the front motor 36 is provided corresponding to the front rolling stand 30 as a first motor device.
  • the rear motor 37 is provided corresponding to the rear rolling stand 33 as a second motor device.
  • the tension meter 40 is provided between the front rolling stand 30 and the rear rolling stand 33.
  • the mathematical model calculation unit 23 receives the deviation between the torque response value of the rear motor 37 and the torque response value of the front motor 36 and calculates the mathematical model using the tension response value of the tensiometer 40 as an output.
  • the mathematical model is calculated based on the history of the actual measurement values of the rolling equipment. For this reason, a mathematical model with high accuracy can be calculated for the rolling equipment.
  • the deviation between the value obtained by multiplying the rotational speed of the rear motor 37 by the radius of the rear rolling roll 34 and the value obtained by multiplying the rotational speed of the front motor 36 by the radius of the front rolling roll 31 is input. Also good.
  • a transfer function from the peripheral speed of each roll to the tension of the rolled material 41 may be calculated as a transfer function corresponding to the mathematical model. Also in this case, a highly accurate mathematical model can be calculated.

Landscapes

  • Controlling Rewinding, Feeding, Winding, Or Abnormalities Of Webs (AREA)
  • Control Of Metal Rolling (AREA)

Abstract

精度の高い数学モデルを算出することができるシート材生産ラインの数学モデル算出装置を提供する。シート材生産ラインの数学モデル算出装置は、シート材生産ラインの実測値の履歴に基づいて、複数のロール装置の間に配置されて前記複数のロール装置の一方から他方へ向けて進行するシート材の張力を計測する張力計の計測値に基づいて制御されるシート材生産ラインに対し、前記複数のロール装置の他方を駆動する第２モータ装置のトルクと前記複数のロール装置の一方を駆動する第１モータ装置のトルクとの偏差を入力とし、前記張力計の計測値を出力として数学モデルを算出する数学モデル算出部、を備えた。

Description

シート材生産ラインの数学モデル算出装置および制御装置
　この発明は、シート材生産ラインの数学モデル算出装置および制御装置に関する。
　特許文献１は、シート材生産ラインの張力制御装置を開示する。当該張力制御装置は、誘導電動機の界磁遅れ時間を考慮したワインダ径の結果に基づいてシート材の張力を制御する。
日本特開２００９－１１３９１１号公報
　しかしながら、特許文献１に記載の張力制御装置においては、張力計に計測値に含まれるノイズの影響により制御器のゲインを上げることができない。このため、シート材の張力制御の精度を上げることができる。
　この発明は、上述の課題を解決するためになされた。この発明の目的は、シート材の張力制御に対する精度の高い数学モデルを算出することができるシート材生産ラインの数学モデル算出装置および制御装置を提供することである。
　この発明に係るシート材生産ラインの数学モデル算出装置は、シート材生産ラインの実測値の履歴に基づいて、複数のロール装置の間に配置されて前記複数のロール装置の一方から他方へ向けて進行するシート材の張力を計測する張力計の計測値に基づいて制御されるシート材生産ラインに対し、前記複数のロール装置の他方を駆動する第２モータ装置のトルクと前記複数のロール装置の一方を駆動する第１モータ装置のトルクとの偏差を入力とし、前記張力計の計測値を出力として数学モデルを算出する数学モデル算出部、を備えた。
　この発明に係るシート材生産ラインの数学モデル算出装置は、シート材生産ラインの実測値の履歴に基づいて、複数のロール装置の間に配置されて前記複数のロール装置の一方から他方へ向けて進行するシート材の張力を計測する張力計の計測値に基づいて制御されるシート材生産ラインに対し、前記複数のロール装置の他方を駆動する第２モータ装置の回転角速度に前記複数のロール装置の他方の半径を乗じた値と前記複数のロール装置の一方を駆動する第１モータ装置の回転角速度に前記複数のロール装置の一方の半径を乗じた値との偏差を入力とし、前記張力計の計測値を出力として数学モデルを算出する数学モデル算出部、を備えた。
　この発明に係るシート材生産ラインの制御装置は、前記シート材生産ラインの稼働中において、前記数学モデル算出装置が前記第１モータ装置のトルクと前記第２モータ装置のトルクとを入力として数学モデルに基づいて算出した張力予測値を取得し、当該張力予測値に基づいて前記第１モータ装置のトルクを制御する制御部、を備えた。
　これらの発明によれば、シート材生産ラインの実測値の履歴に基づいて、数学モデルが算出される。このため、シート材の張力制御に対する精度の高い数学モデルを算出することができる。
この発明の実施の形態１における数学モデル算出装置が適用されるワインダスリット設備の構成図である。 この発明の実施の形態１における数学モデル算出装置が適用されるワインダスリット設備の制御を説明するための図である。 この発明の実施の形態１における数学モデル算出装置により算出された数学モデルを用いた紙材の張力制御を説明するためのブロック図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の概要を説明するためのブロック図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の概要を説明するためのブロック図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の概要を説明するためのブロック図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の例を説明するための図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の例を説明するための図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の例を説明するための図である。 この発明の実施の形態１における数学モデル算出装置による数学モデルに基づいた張力予側値の例を説明するための図である。 この発明の実施の形態１における数学モデル算出装置を備えた制御装置のハードウェア構成図である。 この発明の実施の形態２におけるシート材生産ラインの数学モデル算出装置が算出した数学モデルによるシミュレーション結果を示す図である。 この発明の実施の形態３における数学モデル算出装置が適用される圧延設備の構成図である。
　この発明を実施するための形態について添付の図面に従って説明する。なお、各図中、同一または相当する部分には同一の符号が付される。当該部分の重複説明は適宜に簡略化ないし省略する。
実施の形態１．
　図１はこの発明の実施の形態１における数学モデル算出装置が適用されるワインダスリット設備の構成図である。
　図１のシート材生産ラインは、ワインダスリット設備である。図１に示されるように、アンワインダ１は、ロール装置の一方としてワインダスリット設備の上流側に設けられる。ワインダ２は、ロール装置の他方としてワインダスリット設備の下流側に設けられる。ワインダ２は、フロントドラム３とリアドラム４とライダーロール５とを備える。スリッタ６は、アンワインダ１とワインダ２との間に設けられる。中間ロール７は、アンワインダ１とスリッタ６との間に設けられる。
　アンワインダ用モータ８は、第１モータ装置としてアンワインダ１に対応して設けられる。フロントドラム用モータ９は、第２モータ装置の一部としてフロントドラム３に対応して設けられる。リアドラム用モータ１０は、第２モータ装置の一部としてリアドラム４に対応して設けられる。ライダーロール用モータ１１は、第２モータ装置の一部としてライダーロール５に対応して設けられる。スリッタ用モータ１２は、スリッタ６に対応して設けられる。中間ロール用モータ１３は、中間ロール７に対応して設けられる。
　アンワインダ用ドライブ装置１４は、アンワインダ１に対応して設けられる。フロントドラム用ドライブ装置１５は、フロントドラム３に対応して設けられる。リアドラム用ドライブ装置１６は、リアドラム４に対応して設けられる。ライダーロール用ドライブ装置１７は、ライダーロール５に対応して設けられる。スリッタ用ドライブ装置１８は、スリッタ６に対応して設けられる。中間ロール用ドライブ装置１９は、中間ロール７に対応して設けられる。
　張力計２０は、アンワインダ１とワインダ２との間においてスリッタ６の下流側に設けられる。
　制御装置２１は、数学モデル算出装置２２を備える。数学モデル算出装置２２は、数学モデル算出部２３を備える。数学モデル算出部２３は、ワインダスリット設備の実測値の履歴に基づいて、ワインダスリット設備の数学モデルを算出し得るように設けられる。
　制御装置２１は、制御部２４を備える。制御部２４は、数学モデル算出部２３により算出された数学モデルに基づいてワインダスリット設備の動作を制御し得るように設けられる。
　ワインダスリット設備において、紙材２５は、アンワインダ１とワインダ２との間に配置される。紙材２５は、アンワインダ１からワインダ２に向けて進行する。
　紙材２５は、シート材としてアンワインダ１に予め巻き付いている。アンワインダ１において、紙材２５の巻き付き径は大きい。紙材２５の重量は重い。紙材２５は、中間ロール７を介してスリッタ６で予め設定された幅に裁断される。紙材２５は、ワインダ２に巻き取られる。ワインダ２において、紙材２５は、予め設定された外径まで巻き取られる。この際の紙材２５の巻き取り径は、アンワインダ１における紙材２５の巻き付き径よりも小さい。
　次に、図２を用いて、ワインダスリット設備の制御を説明する。
　図２はこの発明の実施の形態１における数学モデル算出装置が適用されるワインダスリット設備の制御を説明するための図である。
　図２に示されるように、制御装置２１の入力部は、張力計２０の出力部とフロントドラム用ドライブ装置１５の出力部とリアドラム用ドライブ装置１６の出力部とライダーロール用ドライブ装置１７の出力部に接続される。制御装置２１の出力部は、アンワインダ用ドライブ装置１４の装置の入力部とフロントドラム用ドライブ装置１５の入力部とリアドラム用ドライブ装置１６の入力部とライダーロール用ドライブ装置１７の入力部とスリッタ用ドライブ装置１８の入力部と中間ロール用ドライブ装置１９の入力部とに接続される。
　フロントドラム３においては、速度制御が行われる。具体的には、制御部２４は、フロントドラム用ドライブ装置１５に対して回転角速度基準値ωｆｄ ｒｅｆ（ｒａｄ／ｓ）を出力する。フロントドラム用ドライブ装置１５は、回転角速度基準値ωｆｄ ｒｅｆに基づいてフロントドラム用モータ９を駆動する。フロントドラム用モータ９は、フロントドラム３を駆動する。
　リアドラム４においては、速度制御が行われる。具体的には、制御部２４は、リアドラム用ドライブ装置１６に対して回転角速度基準値ωｒｄ ｒｅｆ（ｒａｄ／ｓ）を出力する。リアドラム用ドライブ装置１６は、回転角速度基準値ωｒｄ ｒｅｆに基づいてリアドラム用モータ１０を駆動する。リアドラム用モータ１０は、リアドラム４を駆動する。
　ライダーロール５においては、速度制御が行われる。具体的には、制御部２４は、ライダーロール用ドライブ装置１７に対して回転角速度基準値ωｒｒ ｒｅｆ（ｒａｄ／ｓ）を出力する。ライダーロール用ドライブ装置１７は、回転角速度基準値ωｒｒ ｒｅｆに基づいてライダーロール用モータ１１を駆動する。ライダーロール用モータ１１は、ライダーロール５を駆動する。
　スリッタ６においては、速度制御が行われる。具体的には、制御部２４は、スリッタ用ドライブ装置１８に対して回転角速度基準値ωｓｌ ｒｅｆ（ｒａｄ／ｓ）を出力する。スリッタ用ドライブ装置１８は、回転角速度基準値ωｓｌ ｒｅｆに基づいてスリッタ用モータ１２を駆動する。スリッタ用モータ１２は、スリッタ６を駆動する。
　中間ロール７においては、速度制御が行われる。具体的には、制御部２４は、中間ロール用ドライブ装置１９に対して回転角速度基準値ωｉｎｔ ｒｅｆ（ｒａｄ／ｓ）を出力する。中間ロール用ドライブ装置１９は、回転角速度基準値ωｉｎｔ ｒｅｆに基づいて中間ロール用モータ１３を駆動する。中間ロール用モータ１３は、中間ロール７を駆動する。
　制御部２４は、フロントドラム用ドライブ装置１５からトルク基準値τｆｄ ｒｅｆ（Ｎ・ｍ）の入力を受け付ける。制御部２４は、リアドラム用ドライブ装置１６からトルク基準値τｒｄ ｒｅｆ（Ｎ・ｍ）の入力を受け付ける。制御部２４は、ライダーロール用ドライブ装置１７からトルク基準値τｒｒ ｒｅｆ（Ｎ・ｍ）の入力を受け付ける。制御部２４は、張力計２０の計測値である張力応答値Ｔｒｅｓ（ＭＰａ）の入力を受け付ける。
　アンワインダ１においては、張力制御が行われる。具体的には、制御部２４は、トルク基準値τｆｄ ｒｅｆとトルク基準値τｒｄ ｒｅｆとトルク基準値τｒｒ ｒｅｆと張力応答値Ｔｒｅｓとに基づいてアンワインダ用ドライブ装置１４に対してトルク基準値τｕｗ ｒｅｆ（Ｎ・ｍ）を出力する。アンワインダ用ドライブ装置１４は、トルク基準値τｕｗ ｒｅｆに基づいてアンワインダ用モータ８を駆動する。アンワインダ用モータ８は、アンワインダ１を駆動する。
　その結果、紙材２５は、ワインダ２において巻き取られる際に一定の張力を受ける。紙材２５にかかる張力が一定であると、紙材２５のしわと紙材２５の切断とが抑制される。紙材２５のしわと紙材２５の切断とが抑制されると、ワインダスリット設備の操業が安定する。ワインダスリット設備の操業が安定すると、製品品質が向上する。
　次に、図３を用いて、紙材２５の張力制御を説明する。
　図３はこの発明の実施の形態１における数学モデル算出装置により算出された数学モデルを用いた紙材の張力制御を説明するためのブロック図である。
　図３に示されるように、制御部２４は、張力予測器２６とＰＩ制御器２７とを備える。
　張力予測器２６は、数学モデルに対応した伝達関数ＧＴ（ｓ）を用いてトルク基準値τｆｄ ｒｅｆとトルク基準値τｒｄ ｒｅｆとトルク基準値τｒｒ ｒｅｆとの和とトルク基準値τｕｗ ｒｅｆとの偏差から張力予測値Ｔｍ（ＭＰａ）を算出する。
　ＰＩ制御器２７は、比例ゲインＫｐと積分ゲインＫＩとラプラス演算子ｓとに基づいて張力基準値Ｔｒｅｆ（ＭＰａ）と張力予測値Ｔｍとの偏差からトルク基準値τｕｗ ｒｅｆを算出する。
　この際、アンワインダ用ドライブ装置１４とアンワインダ用モータ８とは、トルク基準値τｕｗ ｒｅｆに基づいて電流制御を行う。アンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓ（Ｎ・ｍ）は、当該電流制御に基づいて決定する。アンワインダ１の回転角速度応答値ωｕｗ ｒｅｓ（ｒａｄ／ｓ）は、アンワインダ１の慣性モーメントＪｕｗ（ｋｇｍ２）とラプラス演算子ｓとに基づいてトルク応答値τｕｗ ｒｅｓと外乱トルクτｄｉｓ（Ｎ・ｍ）とから決定される。
　次に、図４から図６を用いて、数学モデルの算出方法の概要を説明する。
　図４から図６はこの発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の概要を説明するためのブロック図である。
　図４に示されるように、数学モデル算出部２３は、実験的に得られたデータから紙材２５にかかる張力を推定するための数学モデルを作成する。具体的には、数学モデル算出部２３は、フロントドラム用モータ９のトルク応答値τｆｄ ｒｅｓ（Ｎ・ｍ）とリアドラム用モータ１０のトルク応答値τｒｄ ｒｅｓ（Ｎ・ｍ）とライダーロール用モータ１１のトルク応答値τｒｒ ｒｅｓ（Ｎ・ｍ）との和とアンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓ（Ｎ・ｍ）との偏差を入力とし、張力計２０の張力応答値Ｔｒｅｓを出力として、入力から出力までの差分方程式、伝達関数、状態方程式を算出する。
　例えば、図５に示されるように、数学モデル算出部２３は、線形差分方程式であるＡＲＭＡＸ（Ａｕｔｏ－Ｒｅｇｒｅｓｓｉｖｅ　Ｍｏｖｉｎｇ　Ａｖｅｒａｇｅ　ｅＸｏｇｏｎｏｕｓ）モデルを算出する。ＡＲＭＡＸモデルは、次の（１）式で表される。
Figure JPOXMLDOC01-appb-M000001
　（１）式のＡ（ｚ）は、次の（２）式で表される。
Figure JPOXMLDOC01-appb-M000002
　（１）式のＢ（ｚ）は、次の（３）式で表される。
Figure JPOXMLDOC01-appb-M000003
　（１）式のＣ（ｚ）は、次の（４）式で表される。
Figure JPOXMLDOC01-appb-M000004
　ＡＲＭＡＸモデルにおいては、多項式有理関数Ｇ（ｚ）が定義される。Ｇ（ｚ）は、入力ｕ（ｋ）から出力ｙ（ｋ）までの伝達関数である。具体的には、Ｇ（ｚ）は、次の（５）式で表される。
Figure JPOXMLDOC01-appb-M000005
　ＡＲＭＡＸモデルにおいては、多項式有理関数Ｈ（ｚ）が定義される。Ｈ（ｚ）は、雑音ｗ（ｋ）から外乱項ｖ（ｋ）までの伝達関数である。具体的には、Ｈ（ｚ）は、次の（６）式で表される。
Figure JPOXMLDOC01-appb-M000006
　その結果、図５のブロック図は、図６のブロック図に変換される。この際、出力ｙ（ｋ）は、次の（７）式で表される。
Figure JPOXMLDOC01-appb-M000007
　現時刻ｋにおける出力ｙ（ｋ）の予測値は、時刻（ｋ－１）までの過去のデータを用いて次の（８）式で表される。
Figure JPOXMLDOC01-appb-M000008
　なお、（８）式の右辺の第２項は、次の（９）式で定義される。
Figure JPOXMLDOC01-appb-M000009
　（８）式が（７）式に代入されると、次の（１０）式が得られる。
Figure JPOXMLDOC01-appb-M000010
　（７）式と（１０）式とにより雑音ｗ（ｋ）が消去されると、次の（１１）式が得られる。
Figure JPOXMLDOC01-appb-M000011
　（１１）式に示されるように、現在の出力は、過去の入力と出力との線形結合として算出される。この際、１段階予測値を用いた予測誤差εは、次の（１２）式で定義される。
Figure JPOXMLDOC01-appb-M000012
　Ａ（ｚ）とＢ（ｚ）とＣ（ｚ）とは、（１２）式を用いた予測誤差法で決定される。具体的には、Ａ（ｚ）とＢ（ｚ）とＣ（ｚ）とは、予測誤差εから構成される評価関数を最小にするように決定される。
　離散時間系において、入力ｕ（ｋ）から出力ｙ（ｋ）までの伝達関数Ｇ（ｚ）は、次の（１３）式で表される。
Figure JPOXMLDOC01-appb-M000013
　連続時間系において、入力ｕ（ｋ）から出力ｙ（ｋ）までの伝達関数Ｇ´（ｓ）は、（１３）式を変換することにより得られる。伝達関数Ｇ´（ｓ）は、次の（１４）式で表される。
Figure JPOXMLDOC01-appb-M000014
　次に、図７から図９を用いて、数学モデルの算出方法の例を説明する。
　図７から図９はこの発明の実施の形態１における数学モデル算出装置による数学モデルの算出方法の例を説明するための図である。
　図７の上段は、フロントドラム用モータ９のトルク応答値τｆｄ ｒｅｓとリアドラム用モータ１０のトルク応答値τｒｄ ｒｅｓとライダーロール用モータ１１のトルク応答値τｒｒ ｒｅｓとの和とアンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓとの偏差を示す。図７の上段において、フロントドラム用モータ９のトルク応答値τｆｄ ｒｅｓとリアドラム用モータ１０のトルク応答値τｒｄ ｒｅｓとライダーロール用モータ１１のトルク応答値τｒｒ ｒｅｓとアンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓとは、定格トルクを１００％として換算された値である。
　図７の下段は、張力計２０の張力応答値Ｔｒｅｓを示す。図７の下段において、張力計２０の張力応答値Ｔｒｅｓは、ローパスフィルタが施された後の値である。
　図８の上段は、図７の上段における開始直後から目標張力に至るまでの過渡的な状態を示す図である。具体的には、図８の上段は、時刻が０（ｓ）から２０（ｓ）までの拡大図である。図８の上段において、点線は、データを直線近似した線である。
　図８の下段は、図７の下段における開始直後から目標張力に至るまでの過渡的な状態を示す図である。具体的には、図８の下段は、時刻が０（ｓ）から２０（ｓ）までの拡大図である。図８の下段において、点線は、データを直線近似した線である。
　図９の上段は、数学モデルに入力する前のデータの前処理として、図８の上段のデータから低周波外乱である平均値と傾きとを取り除いたデータである。
　図９の下段は、数学モデルに入力する前のデータの前処理として、図８の下段のデータから低周波外乱である平均値と傾きとを取り除いたデータである。
　図９の上段と下段とにおいて、データの平均値と傾きとは０である。数学モデル算出部２３は、データの収集後にオフラインで図９の上段と下段とに対応した処理を行う。数学モデル算出部２３は、図９の上段と下段とのデータを用いて伝達関数を算出する。例えば、伝達関数ＧＴ（ｓ）は、以下の（１５）式で表される。
Figure JPOXMLDOC01-appb-M000015
　得られた伝達関数ＧＴ（ｓ）は、次回以降の紙材２５の張力制御に使用される。数学モデル算出部２３は、各モータのトルク基準値から算出した巻き取りと巻き戻しのトルクの偏差のデータから前回の紙材２５において伝達関数ＧＴ（ｓ）を推定する際に入力データから取り除いた平均値と傾きとをオンラインで取り除く。数学モデル算出部２３は、当該値を伝達関数ＧＴ（ｓ）に入力する。数学モデル算出部２３は、伝達関数ＧＴ（ｓ）の出力に対して、伝達関数ＧＴ（ｓ）を算出する際に出力データから取り除いた平均値と傾きを加えた値を張力予測値Ｔｍとして算出する。
　次に、図１０を用いて、数学モデルに基づいた張力予側値Ｔｍの例を説明する。
　図１０はこの発明の実施の形態１における数学モデル算出装置による数学モデルに基づいた張力予側値の例を説明するための図である。
　図１０に示されるように、張力予測値Ｔｍは、実際の張力応答値Ｔｒｅｓに沿うように算出される。その結果、紙材２５の張力は、精度よく制御される。
　以上で説明した実施の形態１によれば、シート材生産ラインの実測値の履歴に基づいて、数学モデルが算出される。このため、シート材の張力制御に対する精度の高い数学モデルを算出することができる。
　具体的には、数学モデルは、フロントドラム用モータ９のトルク応答値τｆｄ ｒｅｓとリアドラム用モータ１０のトルク応答値τｒｄ ｒｅｓとライダーロール用モータ１１のトルク応答値τｒｒ ｒｅｓとの和とアンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓとの偏差を入力とし、張力計２０の張力応答値Ｔｒｅｓを出力として算出される。このため、ワインダスリット設備に対し、精度の高い数学モデルを算出することができる。
　なお、フロントドラム用モータ９のトルク応答値τｆｄ ｒｅｓとリアドラム用モータ１０のトルク応答値τｒｄ ｒｅｓとライダーロール用モータ１１のトルク応答値τｒｒ ｒｅｓとの和とアンワインダ用モータ８のトルク応答値τｕｗ ｒｅｓに対して重み付けを行った上で数学モデルを算出してもよい。この場合、実情に合わせて精度の高い数学モデルを算出することができる。
　また、数学モデルは、入力と出力とを１次直線で近似し、当該入力と当該出力とから対応した１次直線を差し引いた値を用いて算出される。このため、より精度の高い数学モデルを算出することができる。
　また、数学モデルは、張力計２０の計測値にローパスフィルタを施した値を出力として算出される。このため、張力計２０での計測値に含まれるノイズを除去することができる。その結果、より精度の高い数学モデルを算出することができる。
　また、張力計２０の計測値の複数のサンプリング値の平均値を出力として数学モデルを算出してもよい。この場合、張力計２０での計測値に含まれるノイズを除去するができる。その結果、より精度の高い数学モデルを算出することができる。
　また、アンワインダ用モータ８は、数学モデル算出装置２２に算出された数学モデルから出力された張力予測値Ｔｍに基づいて制御される。このため、張力計２０の計測値に含まれるノイズの影響を抑制することができる。その結果、紙材２５の張力制御の精度を高めることができる。
　なお、数学モデルを算出する際、フロントドラム用モータ９とリアドラム用モータ１０とライダーロール用モータ１１とのうちの最も速い回転速度に対応したロールの半径を乗じた値とアンワインダ用モータ８の回転速度にアンワインダ１の半径を乗じた値との偏差を入力としてもよい。この際、数学モデルに対応した伝達関数として、ロールの周速度から紙材２５の張力までの伝達関数を算出すればよい。この場合も、精度の高い数学モデルを算出することができる。
　また、フィルムの生産ラインに実施の形態１の制御装置２１と同様の制御装置をフィルムの生産ラインに適用してもよい。この場合も、精度の高い数学モデルに基づいて張力制御の精度を高めることができる。
　次に、図１１を用いて、制御装置２１の例を説明する。
　図１１はこの発明の実施の形態１における数学モデル算出装置を備えた制御装置のハードウェア構成図である。
　制御装置２１の各機能は、処理回路により実現し得る。例えば、処理回路は、少なくとも１つのプロセッサ２８ａと少なくとも１つのメモリ２８ｂとを備える。例えば、処理回路は、少なくとも１つの専用のハードウェア２９を備える。
　処理回路が少なくとも１つのプロセッサ２８ａと少なくとも１つのメモリ２８ｂとを備える場合、制御装置２１の各機能は、ソフトウェア、ファームウェア、またはソフトウェアとファームウェアとの組み合わせで実現される。ソフトウェアおよびファームウェアの少なくとも一方は、プログラムとして記述される。ソフトウェアおよびファームウェアの少なくとも一方は、少なくとも１つのメモリ２８ｂに格納される。少なくとも１つのプロセッサ２８ａは、少なくとも１つのメモリ２８ｂに記憶されたプログラムを読み出して実行することにより、制御装置２１の各機能を実現する。少なくとも１つのプロセッサ２８ａは、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、ＤＳＰともいう。例えば、少なくとも１つのメモリ２８ｂは、ＲＡＭ、ＲＯＭ、フラッシュメモリ、ＥＰＲＯＭ、ＥＥＰＲＯＭ等の、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、ＤＶＤ等である。
　処理回路が少なくとも１つの専用のハードウェア２９を備える場合、処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、またはこれらの組み合わせで実現される。例えば、制御装置２１の各機能は、それぞれ処理回路で実現される。例えば、制御装置２１の各機能は、まとめて処理回路で実現される。
　制御装置２１の各機能について、一部を専用のハードウェア２９で実現し、他部をソフトウェアまたはファームウェアで実現してもよい。例えば、数学モデル算出部２３の機能については専用のハードウェア２９としての処理回路で実現し、数学モデル算出部２３の機能以外の機能については少なくとも１つのプロセッサ２８ａが少なくとも１つのメモリ２８ｂに格納されたプログラムを読み出して実行することにより実現してもよい。
　このように、処理回路は、ハードウェア２９、ソフトウェア、ファームウェア、またはこれらの組み合わせで制御装置２１の各機能を実現する。
実施の形態２．
　図１２はこの発明の実施の形態２におけるシート材生産ラインの数学モデル算出装置が算出した数学モデルによるシミュレーション結果を示す図である。なお、実施の形態１の部分と同一又は相当部分には同一符号が付される。当該部分の説明は省略される。
　実施の形態２において、数学モデル算出部２３は、算出した各伝達関数を状態方程式に変換し、カルマンフィルタを用いることによって予測精度を改善する。カルマンフィルタは、システムの動特性を示す状態方程式と時々刻々与えられる実測データとを用いてシステムの状態を逐次的に推定するオンラインアルゴリズムである。カルマンフィルタは、システムに加わる外乱とセンサーとに含まれるノイズの影響が正規分布に従うとの仮定を用いる。
　カルマンフィルタの実装には、制御対象の動特性を状態方程式で表現する必要がある。そこで、数学モデル算出部２３は、実施の形態１で算出された二次系の伝達関数を状態方程式に変換する。伝達関数から状態方程式への変換は無限通りの方法がある。例えば、以下の（１６）式の二次系の伝達関数は、（１7）式の状態方程式で表現される。
Figure JPOXMLDOC01-appb-M000016
Figure JPOXMLDOC01-appb-M000017
　なお、（１７）式において、ｘは状態変数である。ｕは入力（トルク差）である。ｙは出力（張力）である。
　カルマンフィルタは、１サンプリング直前までの情報と現時刻で取得した情報とに基づいてシステムの最適な状態を推定する。
　ただし、状態変数ｘと出力ｙとは、ノイズを含むとされる。
　カルマンフィルタにおいては、サンプリング時刻が更新される度に、予測処理と更新処理とが行われる。
　予測処理においては、現時刻での状態は、１サンプリング前の時刻の情報に基づいて推定される。
　予測処理においては、事前状態推定は、次の（１８）式で表される。
Figure JPOXMLDOC01-appb-M000018
　予測処理においては、事前誤差共分散は、次の（１９）式で表される。
Figure JPOXMLDOC01-appb-M000019
　更新処理においては、正確な状態は、現時刻の実測値に基づいて張力予測値を修正することにより推定される。
　更新処理においては、カルマンゲインは、次の（２０）式で表される。
Figure JPOXMLDOC01-appb-M000020
　更新処理においては、状態推定は、次の（２１）式で表される。
Figure JPOXMLDOC01-appb-M000021
　更新処理においては、事後誤差共分散は、次の（２２）式で表される。
Figure JPOXMLDOC01-appb-M000022
　（１８）式から（２２）式の演算は逐次行われる。この際、張力予測値Ｔｃ ｍは、次の（２３）式で表される。
Figure JPOXMLDOC01-appb-M000023
　数学モデル算出部２３は、実施の形態１で算出された各伝達関数を（１７）式の状態方程式に変換し、カルマンフィルタを実装する。数学モデル算出部２３は、サンプリング毎に取得する張力応答値Ｔｒｅｓと張力予測値Ｔｃ ｍとの偏差を用いて逐次張力予測値Ｔｃ ｍを補正する。具体的には、数学モデル算出部２３は、（１８）式から（２２）式の演算をサンプリング毎に行うことにより張力予測値Ｔｃ ｍを算出する。
　以上で説明した実施の形態２によれば、数学モデル算出部２３は、数学モデルを状態方程式に変換し、カルマンフィルタを用いて、ワインダスリッタ設備の稼働中における張力計２０の計測値に基づいて数学モデルによる張力予測値Ｔｃ ｍを補正する。このため、紙材２５の張力制御の精度をより高めることができる。
　また、フィルムの生産ラインに実施の形態１の制御装置２１と同様の制御装置をフィルムの生産ラインに適用してもよい。この場合も、精度の高い数学モデルに基づいて張力制御の精度を高めることができる。
実施の形態３．
　図１３はこの発明の実施の形態３における数学モデル算出装置が適用される圧延設備の構成図である。なお、実施の形態１の部分と同一又は相当部分には同一符号が付される。当該部分の説明は省略される。
　図１３に示されるように、前方圧延スタンド３０は、ロール装置の一方として圧延設備の上流側に設けられる。前方圧延ロール３１は、前方圧延スタンド３０に設けられる。前方圧下装置３２は、前方圧延スタンド３０の上方に設けられる。
　後方圧延スタンド３３は、ロール装置の他方として圧延設備の下流側に設けられる。後方圧延ロール３４は、後方圧延スタンド３３に設けられる。後方圧下装置３５は、後方圧延スタンド３３の上方に設けられる。
　前方モータ３６は、第１モータ装置として前方圧延スタンド３０に対応して設けられる。後方モータ３７は、第２モータ装置として後方圧延スタンド３３に対応して設けられる。
　前方モータ用ドライブ装置３８は、前方圧延スタンド３０に対応して設けられる。後方モータ用ドライブ装置３９は、後方圧延スタンド３３に対応して設けられる。
　張力計４０は、前方圧延スタンド３０と後方圧延スタンド３３との間に設けられる。
　前方圧延スタンド３０と後方圧延スタンド３３との間において、圧延材４１の張力は、後方モータ３７のトルクと前方モータ３６のトルクとの偏差に基づいて決定される。
　数学モデル算出部２３は、後方モータ３７のトルク応答値と前方モータ３６のトルク応答値との偏差を入力とし、張力計４０の張力応答値を出力として数学モデルを算出する。
　以上で説明した実施の形態３によれば、圧延設備の実測値の履歴に基づいて、数学モデルが算出される。このため、圧延設備に対し、精度の高い数学モデルを算出することができる。
　具体的には、数学モデルは、後方モータ３７のトルク規準値と前方モータ３６のトルク規準値との偏差を入力とし、張力計４０の張力応答値を出力として算出される。このため、圧延設備に対し、精度の高い数学モデルを算出することができる。
　また、数学モデルを算出する際、後方モータ３７の回転速度に後方圧延ロール３４の半径を乗じた値と前方モータ３６の回転速度に前方圧延ロール３１の半径を乗じた値との偏差を入力としてもよい。この際、数学モデルに対応した伝達関数として、各ロールの周速度から圧延材４１の張力までの伝達関数を算出すればよい。この場合も、精度の高い数学モデルを算出することができる。
　なお、実施の形態３において、実施の形態２と同様にカルマンフィルタを用いてもよい。この場合、圧延材４１の張力制御の精度をより高めることができる。
　以上のように、この発明に係るシート材生産ラインの数学モデル算出装置および制御装置は、精度の高い数学モデルを算出するシステムに利用できる。
　１　アンワインダ、　２　ワインダ、　３　フロントドラム、　４　リアドラム、　５　ライダーロール、　６　スリッタ、　７　中間ロール、　８　アンワインダ用モータ、　９　フロントドラム用モータ、　１０　リアドラム用モータ、　１１　ライダーロール用モータ、　１２　スリッタ用モータ、　１３　中間ロール用モータ、　１４　アンワインダ用ドライブ装置、　１５　フロントドラム用ドライブ装置、　１６　リアドラム用ドライブ装置、　１７　ライダーロール用ドライブ装置、　１８　スリッタ用ドライブ装置、　１９　中間ロール用ドライブ装置、　２０　張力計、　２１　制御装置、　２２　数学モデル算出装置、　２３　数学モデル算出部、　２４　制御部、　２５　紙材、　２６　張力予測器、　２７　ＰＩ制御器、　２８ａ　プロセッサ、　２８ｂ　メモリ、　２９　ハードウェア、　３０　前方圧延スタンド、　３１　前方圧延ロール、　３２　前方圧下装置、　３３　後方圧延スタンド、　３４　後方圧延ロール、　３５　後方圧下装置、　３６　前方モータ、　３７　後方モータ、　３８　前方モータ用ドライブ装置、　３９　後方モータ用ドライブ装置、　４０　張力計、　４１　圧延材

Claims (12)

  1. 　シート材生産ラインの実測値の履歴に基づいて、複数のロール装置の間に配置されて前記複数のロール装置の一方から他方へ向けて進行するシート材の張力を計測する張力計の計測値に基づいて制御されるシート材生産ラインに対し、前記複数のロール装置の他方を駆動する第２モータ装置のトルクと前記複数のロール装置の一方を駆動する第１モータ装置のトルクとの偏差を入力とし、前記張力計の計測値を出力として数学モデルを算出する数学モデル算出部、
    を備えたシート材生産ラインの数学モデル算出装置。
  2. 　前記複数のロール装置の一方は、シート材を巻き戻すアンワインダであり、
    　前記複数のロール装置の他方は、フロントドラムとリアドラムとライダーロールとを用いてシート材を巻き取るワインダであり、
    　前記第１モータ装置は、前記アンワインダを駆動するモータを有し、
    　前記第２モータ装置は、前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータを有し、
    　前記数学モデル算出部は、前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータのトルクの和と前記アンワインダを駆動するモータのトルクとの偏差を入力として数学モデルを算出する請求項１に記載のシート材生産ラインの数学モデル算出装置。
  3. 　前記数学モデル算出部は、前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータのトルクに対して重み付けを行った上で前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータのトルクの和と前記アンワインダを駆動するモータのトルクとの偏差を入力として数学モデルを算出する請求項２に記載のシート材生産ラインの数学モデル算出装置。
  4. 　前記複数のロール装置は、複数の圧延スタンドであり、
    　前記第１モータ装置は、前記複数の圧延スタンドの一方のロールを駆動するモータを有し、
    　前記第２モータ装置は、前記複数の圧延スタンドの他方のロールを駆動するモータを有し、
    　前記数学モデル算出部は、前記複数の圧延スタンドの他方のロールを駆動するモータのトルクと前記複数の圧延スタンドの一方のロールを駆動するモータのトルクとの偏差を入力として数学モデルを算出する請求項１に記載のシート材生産ラインの数学モデル算出装置。
  5. 　シート材生産ラインの実測値の履歴に基づいて、複数のロール装置の間に配置されて前記複数のロール装置の一方から他方へ向けて進行するシート材の張力を計測する張力計の計測値に基づいて制御されるシート材生産ラインに対し、前記複数のロール装置の他方を駆動する第２モータ装置の回転角速度に前記複数のロール装置の他方の半径を乗じた値と前記複数のロール装置の一方を駆動する第１モータ装置の回転角速度に前記複数のロール装置の一方の半径を乗じた値との偏差を入力とし、前記張力計の計測値を出力として数学モデルを算出する数学モデル算出部、
    を備えたシート材生産ラインの数学モデル算出装置。
  6. 　前記複数のロール装置の一方は、シート材が巻かれているアンワインダであり、
    　前記複数のロール装置の他方は、フロントドラムとリアドラムとライダーロールとを用いてシート材を巻き取るワインダであり、
    　前記第１モータ装置は、前記アンワインダを駆動するモータを有し、
    　前記第２モータ装置は、前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータを有し、
    　前記数学モデル算出部は、前記フロントドラムと前記リアドラムと前記ライダーロールとをそれぞれ駆動する複数のモータの回転速度のうちの最も速い回転角速度に対応したロールの半径を乗じた値と前記アンワインダを駆動するモータの回転角速度に前記アンワインダの半径を乗じた値との偏差を入力として数学モデルを算出する請求項５に記載のシート材生産ラインの数学モデル算出装置。
  7. 　前記複数のロール装置は、複数の圧延スタンドであり、
    　前記第１モータ装置は、前記複数の圧延スタンドの一方のロールを駆動するモータを有し、
    　前記第２モータ装置は、前記複数の圧延スタンドの他方のロールを駆動するモータを有し、
    　前記数学モデル算出部は、前記複数の圧延スタンドの他方を駆動するモータの回転角速度に前記複数の圧延スタンドの他方のロールの半径を乗じた値と前記複数の圧延スタンドの一方を駆動するモータの回転角速度に前記複数の圧延スタンドの一方のロールの半径を乗じた値との偏差を入力として数学モデルを算出する請求項５に記載のシート材生産ラインの数学モデル算出装置。
  8. 　前記数学モデル算出部は、入力と出力とを１次直線で近似し、当該入力と当該出力とから対応した１次直線を差し引いた値を用いて数学モデルを算出する請求項１から請求項７のいずれか一項に記載のシート材生産ラインの数学モデル算出装置。
  9. 　前記数学モデル算出部は、前記張力計の計測値にローパスフィルタを施した値を出力として数学モデルを算出する請求項１から請求項８のいずれか一項に記載のシート材生産ラインの数学モデル算出装置。
  10. 　前記数学モデル算出部は、前記張力計の計測値の複数のサンプリング値の平均値を出力として数学モデルを算出する請求項１から請求項９のいずれか一項に記載のシート材生産ラインの数学モデル算出装置。
  11. 　前記シート材生産ラインの稼働中において、請求項１から請求項１０のいずれか一項に記載の数学モデル算出装置が前記第１モータ装置のトルクと前記第２モータ装置のトルクとを入力として数学モデルに基づいて算出した張力予測値を取得し、当該張力予測値に基づいて前記第１モータ装置のトルクを制御する制御部、
    を備えたシート材生産ラインの制御装置。
  12. 　前記数学モデル算出装置は、数学モデルを状態方程式に変換し、カルマンフィルタを用いて、前記シート材生産ラインの稼働中における前記張力計の計測値に基づいて数学モデルによる張力予測値を補正する請求項１１に記載のシート材生産ラインの制御装置。
PCT/JP2017/007823 2017-02-28 2017-02-28 シート材生産ラインの数学モデル算出装置および制御装置 WO2018158828A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
KR1020197013943A KR102186020B1 (ko) 2017-02-28 2017-02-28 시트재 생산 라인의 수학 모델 산출 장치 및 제어 장치
CN201780063842.3A CN109843764B (zh) 2017-02-28 2017-02-28 片材生产线的数学模型计算装置及控制装置
PCT/JP2017/007823 WO2018158828A1 (ja) 2017-02-28 2017-02-28 シート材生産ラインの数学モデル算出装置および制御装置
JP2019502326A JP6835200B2 (ja) 2017-02-28 2017-02-28 シート材生産ラインの数学モデル算出装置および制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2017/007823 WO2018158828A1 (ja) 2017-02-28 2017-02-28 シート材生産ラインの数学モデル算出装置および制御装置

Publications (1)

Publication Number Publication Date
WO2018158828A1 true WO2018158828A1 (ja) 2018-09-07

Family

ID=63371196

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/007823 WO2018158828A1 (ja) 2017-02-28 2017-02-28 シート材生産ラインの数学モデル算出装置および制御装置

Country Status (4)

Country Link
JP (1) JP6835200B2 (ja)
KR (1) KR102186020B1 (ja)
CN (1) CN109843764B (ja)
WO (1) WO2018158828A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022202686A1 (ja) * 2021-03-24 2022-09-29 ファナック株式会社 制御装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03120526A (ja) * 1989-10-04 1991-05-22 Fuji Photo Film Co Ltd マイクロフイルム検索装置
JPH09110251A (ja) * 1995-10-23 1997-04-28 Yaskawa Electric Corp ウェブ張力制御方法
JPH09276929A (ja) * 1996-04-15 1997-10-28 Sumitomo Metal Ind Ltd 熱延巻取機の制御方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1993015008A1 (en) * 1992-01-27 1993-08-05 Beloit Technologies, Inc. Closed loop control for a web winding machine
DE59306687D1 (de) 1993-08-24 1997-07-10 Beloit Technologies Inc Verfahren und Wickelvorrichtung zum Wickeln von Bahnen
KR100202695B1 (ko) 1995-10-11 1999-06-15 이종수 장력 제어장치
GB0525676D0 (en) * 2005-12-16 2006-01-25 Rue De Int Ltd Roll storage module and method for its operation
JP4556966B2 (ja) * 2007-04-27 2010-10-06 トヨタ自動車株式会社 ウェブの搬送制御方法、及び搬送制御装置
JP4967152B2 (ja) 2007-11-06 2012-07-04 株式会社安川電機 巻取り・巻戻し装置及びその制御方法
CN100572242C (zh) * 2008-01-18 2009-12-23 济南大学 用于预应力模具制造的变张力钢带缠绕装置及其控制方法
WO2010103180A1 (en) * 2009-03-10 2010-09-16 Abb Oy Determining elastic modulus for continuous material web
CN201685285U (zh) * 2010-05-07 2010-12-29 重庆再升科技发展有限公司 分切机
FI123687B (fi) * 2010-11-18 2013-09-30 Abb Oy Menetelmä ja järjestely rullainkäytön yhteydessä

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03120526A (ja) * 1989-10-04 1991-05-22 Fuji Photo Film Co Ltd マイクロフイルム検索装置
JPH09110251A (ja) * 1995-10-23 1997-04-28 Yaskawa Electric Corp ウェブ張力制御方法
JPH09276929A (ja) * 1996-04-15 1997-10-28 Sumitomo Metal Ind Ltd 熱延巻取機の制御方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022202686A1 (ja) * 2021-03-24 2022-09-29 ファナック株式会社 制御装置

Also Published As

Publication number Publication date
KR102186020B1 (ko) 2020-12-03
CN109843764A (zh) 2019-06-04
JP6835200B2 (ja) 2021-02-24
JPWO2018158828A1 (ja) 2019-11-07
CN109843764B (zh) 2020-08-18
KR20190062577A (ko) 2019-06-05

Similar Documents

Publication Publication Date Title
JP6269428B2 (ja) アルミ箔ダブリング圧延設備の制御装置
JP6380650B2 (ja) 圧延材の板厚制御装置
KR101593812B1 (ko) 압연 제어 장치, 압연 제어 방법 및 기록 매체
JP6583559B2 (ja) アンワインダの制御装置
WO2018158828A1 (ja) シート材生産ラインの数学モデル算出装置および制御装置
JP6183299B2 (ja) 圧延設備の制御システム
JP6094494B2 (ja) 圧延材の板厚制御装置
EP2822176B1 (en) Vector control apparatus of induction motor
JP4967152B2 (ja) 巻取り・巻戻し装置及びその制御方法
KR102264946B1 (ko) 압연 라인의 수학 모델 산출 장치 및 제어 장치
JP6478890B2 (ja) 蛇行制御装置、圧延システムおよび蛇行制御方法
JP6962860B2 (ja) 圧延機の板厚制御装置および該方法ならびに圧延システム
JPWO2019053826A1 (ja) 圧延ラインの数学モデル算出装置および制御装置
JP5949658B2 (ja) 通板設備の制御システム
JP6358124B2 (ja) ブロックミルのモータ制御装置
JP6708268B2 (ja) 複数ロール駆動の速度制御システム
JPH09323851A (ja) 巻取装置
JP2001058212A (ja) 帯状材の張力制御方法及び装置
JP2542698B2 (ja) 圧延機の制御装置
JP3564510B2 (ja) ブライドルロールの張力制御方法
JP6221971B2 (ja) ルーパ装置
JPS63154593A (ja) シ−ト状材料の巻き戻し機または巻き取り機の張力制御方法および張力制御装置
JPS59107716A (ja) 可逆式圧延処理ラインの尾端停止制御方法
JPH07323312A (ja) 電動機の速度制御装置及び圧延材の張力制御装置
JPH0650747A (ja) 圧延機におけるコイル径演算方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17898876

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 20197013943

Country of ref document: KR

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2019502326

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17898876

Country of ref document: EP

Kind code of ref document: A1