JP2003014702A - Apparatus and method for measurement of young's modulus of paper - Google Patents
Apparatus and method for measurement of young's modulus of paperInfo
- Publication number
- JP2003014702A JP2003014702A JP2001203093A JP2001203093A JP2003014702A JP 2003014702 A JP2003014702 A JP 2003014702A JP 2001203093 A JP2001203093 A JP 2001203093A JP 2001203093 A JP2001203093 A JP 2001203093A JP 2003014702 A JP2003014702 A JP 2003014702A
- Authority
- JP
- Japan
- Prior art keywords
- paper
- frequency
- modulus
- young
- maximum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
【０００１】[0001]
【発明の属する技術分野】本発明は、シート物のヤング
率を、特に紙のヤング率を求めるデータ処理装置とデー
タ処理方法、およびそのデータ処理方法をコンピュータ
において読み出し可能なプログラムとして格納した記録
媒体に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a data processing apparatus and a data processing method for determining the Young's modulus of a sheet, particularly the Young's modulus of paper, and a recording medium storing the data processing method as a computer-readable program. Regarding
【０００２】[0002]
【従来の技術】従来、複写機、プリンター、原稿送り装
置や印刷機などの機器の設計に際して、それらの中を搬
送される用紙などシート物の挙動を、計算機シミュレー
ションによって解析し、把握する試みが広く、行われて
いる。その際に、紙などのシート物の搬送挙動を解析
し、予め搬送路を最適化しておくことが重要である。実
際の試験においては、搬送されるべき紙などのシート物
が折れ曲がるなどして、その力学特性が変化することが
あり、様々な条件が発生するため、実際の試験で製品の
信頼性を確認するためには、多くの回数が必要になり、
開発期間の長期化やコスト増大につながる。そのため、
ますます計算機シミュレーションのニーズは高まってい
る。2. Description of the Related Art Conventionally, when designing equipment such as copiers, printers, document feeders and printers, it has been attempted to analyze and grasp the behavior of sheets such as paper conveyed in them by computer simulation. Widely practiced. At that time, it is important to analyze the transport behavior of a sheet such as paper and optimize the transport path in advance. In an actual test, the mechanical properties of a sheet such as paper to be conveyed may be bent and the mechanical properties of the product may change, and various conditions occur, so confirm the reliability of the product in an actual test. In order to do that,
This leads to a longer development period and higher costs. for that reason,
The need for computer simulation is increasing more and more.
【０００３】この種の機器においてシート物を搬送する
シミュレーションのためには、ヤング率のデータが必要
となってくる。しかし、ヤング率の測定や算出は、必ず
しも容易なものではない。[0003] Young's modulus data is required for a simulation of transporting a sheet in this type of equipment. However, measuring and calculating Young's modulus is not always easy.
【０００４】紙のヤング率は、紙種、湿度により変化す
るため、そのデータを得るには、多くの計測データを必
要とする。しかも、シミュレーション計算を行う前に、
予め獲得して準備しておかなくてはならない。従って、
ヤング率を求めるための多くのデータ解析を短時間に必
要とする。このような問題を解決すべく従来技術では、
簡便なヤング率の計測装置やその計測方法が試みられて
きた。Since the Young's modulus of paper changes depending on the paper type and humidity, a large amount of measurement data is required to obtain the data. Moreover, before performing the simulation calculation,
You must acquire and prepare in advance. Therefore,
It requires a lot of data analysis in a short time to obtain the Young's modulus. In the prior art to solve such a problem,
A simple Young's modulus measuring device and its measuring method have been tried.
【０００５】例えば、特開平06-003237号公報（公開日1
994年01月11日）は「ヤング率測定装置及びヤング率測
定方法」と題して、「紙葉類のヤング率測定に適したヤ
ング率測定装置及びヤング率測定方法に関し、簡単な構
成でより正確なヤング率を測定できるヤング率測定装置
及びヤング率測定方法」を開示している。しかしなが
ら、同公報には、シミュレーション計算を用いてヤング
率データを大量に短時間に求める方法については、開示
していない。For example, Japanese Patent Laid-Open No. 06-003237 (published date 1
(January 11, 994) titled "Young's modulus measuring device and Young's modulus measuring method", "Young modulus measuring device and Young's modulus measuring method suitable for measuring Young's modulus of paper sheets, A Young's modulus measuring device and Young's modulus measuring method capable of accurately measuring Young's modulus are disclosed. However, this publication does not disclose a method of obtaining a large amount of Young's modulus data in a short time by using a simulation calculation.
【０００６】一般に、ヤング率を求める方法としては、
振動リード法と静的たわみ法があるが、振動リード法の
方が精度も高い。また従来は大型であった装置もノート
パソコンの普及により、環境試験室内に持ち込んで振動
リード法により測定することが可能になっている。Generally, as a method for obtaining the Young's modulus,
There are vibration lead method and static deflection method, but the vibration lead method is more accurate. Also, with the spread of notebook computers, devices that were large in the past can now be brought into environmental test rooms and measured by the vibration lead method.
【０００７】また振動リード法は実際の測定と後のデー
タ処理を分離して行うことが可能である。また、データ
処理もＶｉｓｕａｌ Ｂａｓｉｃによるプログラムによ
り、半自動化されてきた。しかし、まだ完全に自動化は
されてはいない状態である。本発明者は、８種類の紙種
に対して温湿度がわずか９条件であるデータ処理でさ
え、半自動化されたプログラムでは１ヶ月以上費やし
た。このことは、紙の計算機シミュレーションにとって
は非常に重大な問題である。In the vibration lead method, actual measurement and subsequent data processing can be performed separately. Data processing has also been semi-automated by a program by Visual Basic. However, it has not been fully automated yet. The inventor spent more than a month on semi-automated programs, even for data processing where the temperature and humidity were only 9 conditions for 8 paper types. This is a very serious problem for paper computer simulation.
【０００８】[0008]
【発明が解決しようとする課題】本発明は、上記事情に
鑑みてなされたものであって、その目的はヤング率を求
めるためのデータ処理を自動処理するために必要な極大
値の判定方法、およびその方法を用いた装置、そしてそ
の方法がプログラムとして格納された記録媒体を提供す
るものである。SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and an object thereof is a method of determining a maximum value necessary for automatically processing data processing for obtaining Young's modulus, And an apparatus using the method, and a recording medium in which the method is stored as a program.
【０００９】[0009]
【課題を解決するための手段】上記課題は、以下に述べ
る本発明に係る以下の演算処理手段を有する装置、演算
処理工程を有する方法、および演算処理工程をプログラ
ムとして格納した記録媒体によって達成される。The above object can be achieved by an apparatus having the following arithmetic processing means according to the present invention described below, a method having arithmetic processing steps, and a recording medium storing the arithmetic processing steps as a program. It
【００１０】即ち請求項１に記載のデータ処理装置は、
紙の共振周波数を求めて、それを用いて、振動リード法
により紙のヤング率を求める方法において、周波数を横
軸にとり、該横軸において該周波数領域を第１の一定長
さに分割した複数の第１閉区間を設け、該各第１閉区間
における伝達関数の平均値を縦軸にとったグラフを作
り、次に、該周波数領域を該第１の一定長さよりも長い
他の長さを有する第２閉区間に分け、各該第２区間にお
ける前記グラフの最大値を求めて、隣接する該最大値同
士を結ぶ線分の傾きが、正から、負又は０に変わる点を
該グラフ全体の極大点の一つであると判定することを特
徴とする。That is, the data processing apparatus according to claim 1 is
In a method for obtaining the Young's modulus of a paper by the vibration lead method using the resonance frequency of paper and using it, the frequency is plotted on the horizontal axis, and the frequency domain is divided into a first constant length on the horizontal axis. A first closed section is provided, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis is created, and then the frequency domain is set to another length longer than the first fixed length. Is divided into second closed sections, the maximum value of the graph in each of the second sections is obtained, and the point at which the slope of the line segment connecting the adjacent maximum values changes from positive to negative or 0 in the graph. The feature is that it is determined to be one of the maximum points of the whole.
【００１１】このような構成の装置によって、自動的に
極大値を求める事ができるので、ヤング率を用いるシミ
ュレーション計算においてスピードアップを行い、省力
化を達成することができる。Since the maximum value can be automatically obtained by the apparatus having such a configuration, the simulation calculation using the Young's modulus can be speeded up and labor saving can be achieved.
【００１２】請求項２に記載のデータ処理装置は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、周波数を横軸に
とり、該横軸において該周波数領域を第１の一定長さに
分割した複数の第１閉区間を設け、該各第１閉区間にお
ける伝達関数の平均値を縦軸にとったグラフを作り、次
に、該周波数領域を該第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各該第２閉区間における
前記グラフの最大値と最小値とを求めて、隣接する該最
大値同士を結ぶ線分の傾きが正から、負または０に変わ
る点をそのグラフ全体の極大点の候補の１つであると判
定し、かつ、隣接する該最小値同士を結ぶ線分の傾き
が、負から、正または０に変わる点をそのグラフ全体の
極小点の一つであると判定し、更に、ある一つの該極大
について、該極大値を与える周波数の左右で最も近い位
置にある極小値を求め、該極大値と該左右それぞれの極
小値との差Ａ，Ｂを求め、また、該極大値と該左右それ
ぞれの極小値との平均値Ｃ、Ｄを求め、Ａ／Ｃ、Ｂ／Ｄ
のいずれか大なる方がある一定値よりも小さい時、該極
大を共振による極大とは判定しないことを特徴とする。According to a second aspect of the present invention, in the data processing apparatus, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by the vibration lead method using the resonance frequency. A plurality of first closed sections, which are obtained by dividing the frequency domain into a first constant length, are provided, and a graph in which the average value of the transfer function in each of the first closed sections is plotted on the vertical axis is created. Is divided into second closed sections having other lengths longer than the first fixed length, the maximum value and the minimum value of the graph in each of the second closed sections are obtained, and the adjacent maximum values are The point at which the slope of the connecting line segment changes from positive to negative or 0 is determined as one of the candidates for the maximum points of the entire graph, and the slope of the line segment connecting the adjacent minimum values is negative. From one to positive or zero is one of the minimum points in the whole graph Judgment, further, for a certain one of the maxima, obtain the minimum value at the closest positions on the left and right of the frequency giving the maximum value, and obtain the difference A, B between the maximum value and the respective minimum values of the left and right, Further, the average values C and D of the maximum value and the left and right minimum values are calculated to obtain A / C and B / D.
Is larger than a certain value, the maximum is not judged to be a maximum due to resonance.
【００１３】このような構成の装置によって、ねじれモ
ードの極大を除くことができるので、より正確なヤング
率を求める事が自動的にでき、シミュレーション計算に
おいてスピードアップ、省力化を達成することができ
る。Since the maximum of the twist mode can be removed by the device having such a configuration, it is possible to automatically obtain a more accurate Young's modulus, and it is possible to speed up the simulation calculation and save labor. .
【００１４】請求項３に記載のデータ処理装置は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、複数の紙片の長
さから、該紙のヤング率を求める場合、第１の長さにつ
いて該紙のヤング率を求めた後、さらに、第２の長さの
該紙のヤング率を求める時、既知のヤング率から共振周
波数を逆算して得られた周波数を中心に含むある一定周
波数の範囲内において最大値が存在する場合、該最大値
を極大値と判定することを特徴とする。According to a third aspect of the present invention, there is provided a data processing apparatus, wherein a resonance frequency of a paper is obtained and a Young's modulus of the paper is obtained by a vibration lead method using the resonance frequency of the paper. When the Young's modulus of the paper is obtained for the first length and then the Young's modulus of the paper for the second length is obtained, the resonance frequency is calculated backward from the known Young's modulus. When a maximum value exists within a certain constant frequency range including the frequency obtained as a center, the maximum value is determined to be the maximum value.
【００１５】このような構成の装置によって、決められ
た範囲内の最大値を求めるだけなので、ヤング率を使っ
たシミュレーション計算において、さらにスピードアッ
プすることができる。With the apparatus having such a configuration, only the maximum value within the determined range is obtained, so that the simulation calculation using the Young's modulus can be further speeded up.
【００１６】請求項４に記載のデータ処理装置は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、複数の紙片の長
さそれぞれに対して、一つまたは複数の共振周波数から
求められる近似直線の傾きから算出される一つのヤング
率の値に対して、該近似直線からある一定以上離れた共
振周波数を除いて、ヤング率を算出することを特徴とす
る。According to a fourth aspect of the present invention, there is provided a data processing device, wherein a resonance frequency of paper is obtained, and the Young's modulus of the paper is obtained by using the resonance frequency of the paper. Calculating the Young's modulus by excluding the resonance frequency that is a certain distance or more from the approximate straight line, with respect to one Young's modulus value calculated from the slope of the approximate straight line obtained from one or more resonance frequencies Is characterized by.
【００１７】このような構成の装置では、近似直線から
はずれた共振周波数を除くので、さらに正確なヤング率
を求めることができる。In the apparatus having such a configuration, since the resonance frequency deviating from the approximate straight line is excluded, a more accurate Young's modulus can be obtained.
【００１８】請求項５に記載のデータ処理方法は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、周波数を横軸に
とり、該横軸において該周波数領域を第１の一定長さに
分割した複数の第１閉区間を設け、該各第１閉区間にお
ける伝達関数の平均値を縦軸にとったグラフを作り、次
に、該周波数領域を該第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各該第２区間における前
記グラフの最大値を求めて、隣接する該最大値同士を結
ぶ線分の傾きが、正から、負又は０に変わる点を該グラ
フ全体の極大点の一つであると判定することを特徴とす
る。According to a fifth aspect of the data processing method of the present invention, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by the vibration lead method using the resonance frequency. A plurality of first closed sections, which are obtained by dividing the frequency domain into a first constant length, are provided, and a graph in which the average value of the transfer function in each of the first closed sections is plotted on the vertical axis is created. Is divided into second closed sections having other lengths longer than the first fixed length, the maximum value of the graph in each of the second sections is determined, and the slope of a line segment connecting the adjacent maximum values to each other. Is determined to be one of the maximal points of the entire graph, which is changed from positive to negative or zero.
【００１９】このような方法によって、自動的に極大値
を求める事ができるので、ヤング率を用いるシミュレー
ション計算においてスピードアップを行い、省力化を達
成することができる。Since the maximum value can be automatically obtained by such a method, it is possible to speed up the simulation calculation using the Young's modulus and achieve labor saving.
【００２０】請求項６に記載のデータ処理方法は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、周波数を横軸に
とり、該横軸において該周波数領域を第１の一定長さに
分割した複数の第１閉区間を設け、該各第１閉区間にお
ける伝達関数の平均値を縦軸にとったグラフを作り、次
に、該周波数領域を該第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各該第２閉区間における
前記グラフの最大値と最小値とを求めて、隣接する該最
大値同士を結ぶ線分の傾きが正から、負または０に変わ
る点をそのグラフ全体の極大点の候補の１つであると判
定し、かつ、隣接する該最小値同士を結ぶ線分の傾き
が、負から、正または０に変わる点をそのグラフ全体の
極小点の一つであると判定し、更に、ある一つの該極大
について、該極大値を与える周波数の左右で最も近い位
置にある極小値を求め、該極大値と該左右それぞれの極
小値との差Ａ，Ｂを求め、また、該極大値と該左右それ
ぞれの極小値との平均値Ｃ、Ｄを求め、Ａ／Ｃ、Ｂ／Ｄ
のいずれか大なる方がある一定値よりも小さい時、該極
大を共振による極大とは判定しないことを特徴とする。According to a sixth aspect of the present invention, in the data processing method, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by the vibration lead method using the resonance frequency. A plurality of first closed sections, which are obtained by dividing the frequency domain into a first constant length, are provided, and a graph in which the average value of the transfer function in each of the first closed sections is plotted on the vertical axis is created. Is divided into second closed sections having other lengths longer than the first fixed length, the maximum value and the minimum value of the graph in each of the second closed sections are obtained, and the adjacent maximum values are The point at which the slope of the connecting line segment changes from positive to negative or 0 is determined as one of the candidates for the maximum points of the entire graph, and the slope of the line segment connecting the adjacent minimum values is negative. From one to positive or zero is one of the minimum points in the whole graph Judgment, further, for a certain one of the maxima, obtain the minimum value at the closest positions on the left and right of the frequency giving the maximum value, and obtain the difference A, B between the maximum value and the respective minimum values of the left and right, Further, the average values C and D of the maximum value and the left and right minimum values are calculated to obtain A / C and B / D.
Is larger than a certain value, the maximum is not judged to be a maximum due to resonance.
【００２１】このような方法によって、ねじれモードの
極大を除くことができるので、より正確なヤング率を求
める事が自動的にでき、シミュレーション計算において
スピードアップ、省力化を達成することができる。Since the maximum of the twist mode can be removed by such a method, a more accurate Young's modulus can be automatically obtained, and the simulation calculation can be speeded up and labor can be saved.
【００２２】請求項７に記載のデータ処理方法は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、複数の紙片の長
さから、該紙のヤング率を求める場合、第１の長さにつ
いて該紙のヤング率を求めた後、さらに、第２の長さの
該紙のヤング率を求める時、既知のヤング率から共振周
波数を逆算して得られた周波数を中心に含むある一定周
波数の範囲内において最大値が存在する場合、該最大値
を極大値と判定することを特徴とする。According to a seventh aspect of the present invention, in the data processing method, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by using the resonance frequency of the paper. When the Young's modulus of the paper is obtained for the first length and then the Young's modulus of the paper for the second length is obtained, the resonance frequency is calculated backward from the known Young's modulus. When a maximum value exists within a certain constant frequency range including the frequency obtained as a center, the maximum value is determined to be the maximum value.
【００２３】このような方法によって、決められた範囲
内の最大値を求めるだけなので、ヤング率を使ったシミ
ュレーション計算において、さらにスピードアップする
ことができる。By such a method, only the maximum value within the determined range is obtained, so that the speed can be further increased in the simulation calculation using the Young's modulus.
【００２４】請求項８に記載のデータ処理方法は、紙の
共振周波数を求めて、それを用いて、振動リード法によ
り紙のヤング率を求める方法において、複数の紙片の長
さそれぞれに対して、一つまたは複数の共振周波数から
求められる近似直線の傾きから算出される一つのヤング
率の値に対して、該近似直線からある一定以上離れた共
振周波数を除いて、ヤング率を算出することを特徴とす
る。A data processing method according to claim 8 is a method for obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of the paper and using the resonance frequency of the paper, for each length of a plurality of paper pieces. Calculating the Young's modulus by excluding the resonance frequency that is a certain distance or more from the approximate straight line, with respect to one Young's modulus value calculated from the slope of the approximate straight line obtained from one or more resonance frequencies Is characterized by.
【００２５】このような方法では、近似直線からはずれ
た共振周波数を除くので、さらに正確なヤング率を求め
ることができる。In such a method, since the resonance frequency deviating from the approximate straight line is excluded, a more accurate Young's modulus can be obtained.
【００２６】請求項９に記載のコンピュータが読み込み
可能なプログラムを格納した記録媒体は、紙の共振周波
数を求めて、それを用いて、振動リード法により紙のヤ
ング率を求める方法において、周波数を横軸にとり、該
横軸において該周波数領域を第１の一定長さに分割した
複数の第１閉区間を設け、該各第１閉区間における伝達
関数の平均値を縦軸にとったグラフを作り、次に、該周
波数領域を該第１の一定長さよりも長い他の長さを有す
る第２閉区間に分け、各該第２区間における前記グラフ
の最大値を求めて、隣接する該最大値同士を結ぶ線分の
傾きが、正から、負又は０に変わる点を該グラフ全体の
極大点の一つであると判定することを特徴とする。According to a ninth aspect of the present invention, there is provided a recording medium having a computer-readable program stored therein, wherein a resonance frequency of a paper is obtained, and the resonance frequency of the paper is obtained. A graph in which the horizontal axis is provided with a plurality of first closed sections obtained by dividing the frequency domain into first constant lengths on the horizontal axis, and the average value of the transfer function in each of the first closed sections is taken on the vertical axis is shown. Then, the frequency domain is divided into second closed sections having other lengths longer than the first fixed length, the maximum value of the graph in each of the second sections is calculated, and the maximum value of the adjacent maximum It is characterized in that the point at which the inclination of the line segment connecting the values changes from positive to negative or zero is determined as one of the maximum points of the entire graph.
【００２７】このような構成の記録媒体によって、自動
的に極大値を求める事ができるので、ヤング率を用いる
シミュレーション計算においてスピードアップを行い、
省力化を達成することができる。Since the maximum value can be automatically obtained by the recording medium having such a structure, the speed is increased in the simulation calculation using the Young's modulus.
Labor saving can be achieved.
【００２８】請求項１０に記載のコンピュータが読み込
み可能なプログラムを格納した記録媒体は、紙の共振周
波数を求めて、それを用いて、振動リード法により紙の
ヤング率を求める方法において、周波数を横軸にとり、
該横軸において該周波数領域を第１の一定長さに分割し
た複数の第１閉区間を設け、該各第１閉区間における伝
達関数の平均値を縦軸にとったグラフを作り、次に、該
周波数領域を該第１の一定長さよりも長い他の長さを有
する第２閉区間に分け、各該第２閉区間における前記グ
ラフの最大値と最小値とを求めて、隣接する該最大値同
士を結ぶ線分の傾きが正から、負または０に変わる点を
そのグラフ全体の極大点の候補の１つであると判定し、
かつ、隣接する該最小値同士を結ぶ線分の傾きが、負か
ら、正または０に変わる点をそのグラフ全体の極小点の
一つであると判定し、更に、ある一つの該極大につい
て、該極大値を与える周波数の左右で最も近い位置にあ
る極小値を求め、該極大値と該左右それぞれの極小値と
の差Ａ，Ｂを求め、また、該極大値と該左右それぞれの
極小値との平均値Ｃ、Ｄを求め、Ａ／Ｃ、Ｂ／Ｄのいず
れか大なる方がある一定値よりも小さい時、該極大を共
振による極大とは判定しないことを特徴とする。According to a tenth aspect of the present invention, there is provided a recording medium having a computer-readable program stored therein, wherein a resonance frequency of a paper is obtained, and the resonance frequency of the paper is obtained. On the horizontal axis,
A plurality of first closed sections obtained by dividing the frequency domain into a first constant length are provided on the horizontal axis, and a graph having the vertical axis of the average value of the transfer function in each of the first closed sections is created. , The frequency domain is divided into second closed sections having other lengths longer than the first fixed length, the maximum value and the minimum value of the graph in each of the second closed sections are obtained, and the adjacent The point where the slope of the line segment connecting the maximum values changes from positive to negative or 0 is determined to be one of the candidates for the maximum point of the entire graph,
And, the slope of the line segment connecting the adjacent minimum values is determined to be one of the minimum points of the entire graph from the point where the slope changes from negative to positive or 0, and further, for one certain maximum, The local minimum values at the closest positions on the left and right of the frequency giving the local maximum value are obtained, the differences A and B between the local maximum value and the respective left and right local minimum values are obtained, and the local maximum value and the respective left and right local minimum values are obtained. It is characterized in that the average values C and D are calculated, and when either A / C or B / D, whichever is larger, is smaller than a certain value, the maximum is not judged to be the maximum due to resonance.
【００２９】このようなプログラムを有する記録媒体に
よって、ねじれモードの極大を除くことができるので、
より正確なヤング率を求める事が自動的にでき、シミュ
レーション計算においてスピードアップ、省力化を達成
することができる。Since the maximum of the twist mode can be eliminated by the recording medium having such a program,
It is possible to automatically obtain a more accurate Young's modulus, which can speed up and save labor in simulation calculation.
【００３０】請求項１１に記載のコンピュータが読み込
み可能なプログラムを格納した記録媒体は、紙の共振周
波数を求めて、それを用いて、振動リード法により紙の
ヤング率を求める方法において、複数の紙片の長さか
ら、該紙のヤング率を求める場合、第１の長さについて
該紙のヤング率を求めた後、さらに、第２の長さの該紙
のヤング率を求める時、既知のヤング率から共振周波数
を逆算して得られた周波数を中心に含むある一定周波数
の範囲内において最大値が存在する場合、該最大値を極
大値と判定することを特徴とする。A recording medium storing a computer-readable program according to claim 11 is a method of obtaining a Young's modulus of a paper by a vibration read method by obtaining a resonance frequency of the paper and using the resonance frequency of the paper. When the Young's modulus of the paper is obtained from the length of the piece of paper, when the Young's modulus of the paper is obtained for the first length and then the Young's modulus of the paper of the second length is obtained, a known value is obtained. When the maximum value exists within a certain frequency range including the frequency obtained by back-calculating the resonance frequency from the Young's modulus, the maximum value is determined to be the maximum value.
【００３１】このような記録媒体に格納されたプログラ
ムは、決められた範囲内の最大値を求めるだけなので、
ヤング率を使ったシミュレーション計算において、さら
にスピードアップすることができる。Since the program stored in such a recording medium only finds the maximum value within a predetermined range,
It is possible to further speed up the simulation calculation using Young's modulus.
【００３２】請求項１２に記載のコンピュータが読み込
み可能なプログラムを格納した記録媒体は、紙の共振周
波数を求めて、それを用いて、振動リード法により紙の
ヤング率を求める方法において、複数の紙片の長さそれ
ぞれに対して、一つまたは複数の共振周波数から求めら
れる近似直線の傾きから算出される一つのヤング率の値
に対して、該近似直線からある一定以上離れた共振周波
数を除いて、ヤング率を算出することを特徴とする。According to a twelfth aspect of the present invention, there is provided a recording medium having a computer-readable program stored therein, wherein a plurality of resonance frequencies of the paper are obtained, and a Young's modulus of the paper is obtained by a vibration read method using the resonance frequency. For each length of paper, for one Young's modulus value calculated from the slope of the approximate straight line obtained from one or more resonant frequencies, exclude the resonant frequency that is a certain distance or more from the approximate straight line. Then, the Young's modulus is calculated.
【００３３】このような記録媒体の有するプログラムで
は、近似直線からはずれた共振周波数を除くので、さら
に正確なヤング率を求めることができる。In the program included in such a recording medium, since the resonance frequency deviated from the approximate straight line is excluded, more accurate Young's modulus can be obtained.
【００３４】請求項１３に記載のデータ処理方法は、紙
の共振周波数を求めて、それを用いて、振動リード法に
より紙のヤング率を求める方法において、周波数を横軸
にとり、該横軸において該周波数領域を第１の一定長さ
に分割した複数の第１閉区間を設け、該各第１閉区間に
おける伝達関数の平均値を縦軸にとったグラフを作り、
次に、該周波数領域を該第１の一定長さよりも長い他の
長さを有する第２閉区間に分け、各該第２区間における
前記グラフの最大値を求めて、隣接する該最大値同士を
結ぶ線分の傾きが、正から、負又は０に変わる点を該グ
ラフ全体の極大点の候補であると判定する方法におい
て、第１次共振による１次極大値の候補をもとにして第
２次共振による２次極大値、３次共振による３次極大値
など、該１次以外の共振による極大値を予想し、該１次
以外の共振による極大値を与えると予想される周波数を
含むある範囲において極大点の候補が存在する場合、該
極大点による極大値の候補を、該２次、３次などの極大
値であると判定することを特徴とする。According to a thirteenth aspect of the present invention, in the data processing method, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by the vibration lead method using the resonance frequency. A plurality of first closed sections provided by dividing the frequency domain into first constant lengths are provided, and a graph in which the average value of the transfer function in each of the first closed sections is plotted on the vertical axis,
Next, the frequency domain is divided into second closed sections having other lengths longer than the first constant length, the maximum value of the graph in each of the second sections is obtained, and the adjacent maximum values are In the method of determining the point at which the slope of the line segment connecting the points from positive to negative or 0 is a candidate for the maximum point of the entire graph, based on the candidates for the primary maximum value due to the first resonance. A frequency that is expected to give a maximum value due to a resonance other than the primary resonance, such as a secondary maximum value due to the second resonance and a third maximum value due to the third resonance, is predicted. When there is a maximum point candidate in a certain range that includes the maximum point, the maximum value candidate based on the maximum point is determined to be the maximum value such as the secondary or tertiary order.
【００３５】このような方法によって、自動的に極大を
求めることができるので、シミュレーション計算におい
てスピードアップを行い、省力化を達成することができ
る。Since the maximum value can be automatically obtained by such a method, the simulation calculation can be speeded up and labor saving can be achieved.
【００３６】請求項１４に記載のデータ処理方法は、紙
の共振周波数を求めて、それを用いて、振動リード法に
より紙のヤング率を求める方法において、周波数を横軸
にとり、該横軸において周波数領域を第１の一定長さに
分割した複数の第１閉区間を設け、該各第１閉区間にお
ける伝達関数の平均値を縦軸にとったグラフを作り、次
に、該周波数領域を該第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各該第２閉区間にける前
記グラフの最大値を求め、隣接する該最大値同士を結ぶ
線分の傾きが正から、負又は０に変わる点をそのグラフ
全体の極大点の候補と判定する方法において、該伝達関
数をもとにしたグラフに基づいてクロス位相１階微分し
たクロス位相１回微分のグラフを作り、該極大の候補を
与える周波数を含む一定の範囲内において、該クロス位
相１階微分値が一定量以上変化している時、該極大点の
候補を極大と判定することを特徴とする。According to a fourteenth aspect of the present invention, in the data processing method, the resonance frequency of the paper is obtained, and the Young's modulus of the paper is obtained by the vibration lead method using the resonance frequency. The frequency domain is divided into a plurality of first constant lengths, a plurality of first closed sections are provided, and a graph in which the vertical axis is the average value of the transfer function in each of the first closed sections is created. Divide into second closed sections having other lengths longer than the first fixed length, find the maximum value of the graph in each of the second closed sections, and tilt the line segment connecting the adjacent maximum values. In a method of determining a point where is changed from positive to negative or to 0 as a candidate for the maximum point of the entire graph, a graph of cross phase one-time differentiation obtained by performing cross phase first derivative based on a graph based on the transfer function. And include the frequencies that give the maximum candidate. In a certain range, when the cross phase the first derivative values change more than a certain amount, and judging the maximum candidates of ultra large point.
【００３７】このような方法によって、自動的に極大を
求めることができ、さらにクロス位相の１階微分を参照
しているので、シミュレーション計算において、正確
に、しかもスピードアップを行い、省力化を達成するこ
とができる。By such a method, the maximum value can be automatically obtained, and since the first derivative of the cross phase is referred to, the simulation calculation can be accurately and speeded up, and the labor saving can be achieved. can do.
【００３８】[0038]
【発明の実施の形態】以下、図面を参照しながら、本発
明の実施の形態について説明する。DETAILED DESCRIPTION OF THE INVENTION Embodiments of the present invention will be described below with reference to the drawings.
【００３９】第１の実施の形態（請求項５）
本発明の第１の実施の形態を図面を見ながら説明する。
まず、図１に、振動リード法により得られる伝達関数の
グラフを、横軸を周波数にとって示す。従来は、このグ
ラフの極大付近の左右の点に印を付けることにより、半
自動的に、極大を求めることができた。この作業を完全
に自動化するのが本発明である。First Embodiment (Claim 5) A first embodiment of the present invention will be described with reference to the drawings.
First, FIG. 1 shows a graph of a transfer function obtained by the vibration lead method, with the horizontal axis representing frequency. Conventionally, by marking the left and right points near the maximum of this graph, the maximum could be calculated semi-automatically. The present invention fully automates this work.
【００４０】図２に本発明の第一の発明を説明するため
のフローチャートを示す。まず、図１に見られる細かい
ノイズの影響を取り除く必要がある。そのために、連続
するある一定長さを有する横軸（周波数）の閉区間を設
定し、その範囲における伝達関数の平均値を求めて、改
めて、グラフにする。この工程をノイズ除去グラフ作成
工程と呼ぶ。FIG. 2 shows a flowchart for explaining the first invention of the present invention. First, it is necessary to remove the effect of fine noise seen in FIG. For that purpose, a continuous closed section on the horizontal axis (frequency) having a certain fixed length is set, the average value of the transfer function in that range is obtained, and the graph is made again. This process is called a noise removal graph creation process.
【００４１】次に、そのグラフを該一定長さよりも長い
他の長さを有する横軸上の区間を設定し、全体を分け
て、一つ一つの区間における該グラフの最大値を求め
る。次に、隣接する区間における最大値を結ぶ線分の傾
きを計算する。周波数の低いところから、高いところへ
順に見ていった時、傾きの符号が正から、負又は０に変
化する点を極大点と判定する。以上の操作によって、極
大点を自動的に発見することが可能となる。Next, a section on the horizontal axis having another length longer than the fixed length is set in the graph, and the whole is divided to obtain the maximum value of the graph in each section. Next, the slope of the line segment connecting the maximum values in the adjacent sections is calculated. A point where the sign of the slope changes from positive to negative or 0 when the frequency is viewed from a low frequency to a high frequency is determined to be a maximum point. With the above operation, the maximum point can be automatically found.
【００４２】第２の実施の形態（請求項６）
次に図３に、本発明の第２の実施の形態を説明するため
のフローチャートを示す。Second Embodiment (Claim 6) Next, FIG. 3 shows a flow chart for explaining the second embodiment of the present invention.
【００４３】図１では、左から２番目の極大点（Ｂ）
は、ねじれのモードの共振なので、共振周波数を求める
ときには、排除しなければならない。図１のような、ね
じれモードの共振による極大点は、山と谷の差が他の極
大よりも小さいのでこれを利用する。まず、細かいノイ
ズをとり除くために連続するある一定の長さを有する区
間に細分した周波数領域での伝達関数の平均値を縦軸に
とり、周波数を横軸に取ったグラフを作成する（第１の
実施の形態の時に説明したのと同じノイズ除去グラフ作
成工程）。In FIG. 1, the second maximum point (B) from the left.
Is a torsional mode resonance and must be eliminated when determining the resonance frequency. The maximum point due to the resonance of the torsion mode as shown in FIG. 1 is used because the difference between the peak and the valley is smaller than the other maximums. First, a graph is created in which the vertical axis represents the average value of the transfer function in the frequency domain subdivided into continuous sections having a certain length to remove fine noise, and the horizontal axis represents the frequency (first The same noise removal graph creation process as described in the embodiment).
【００４４】次にこうして作成されたグラフに対して、
前記一定の長さより長い他の長さを有する横軸上の閉区
間に全体を分け、一つ一つの閉区間における伝達関数の
最大値と最小値を求める。そして、隣接する最大値を結
ぶ線分の傾きが、周波数の低いところから高いところへ
見た時、正から、負または０に変わる点を、そのグラフ
全体の極大の一つであると仮定する。同様に、隣接する
最小値を結ぶ線分の傾きが周波数の低いところから高い
ところへ見た時、負から、正または０に変わる点をその
グラフの極小の一つであると仮定する。Next, for the graph thus created,
The whole is divided into closed sections on the horizontal axis having another length longer than the predetermined length, and the maximum value and the minimum value of the transfer function in each closed section are obtained. Then, it is assumed that the point where the slope of the line segment connecting the adjacent maximum values changes from positive to negative or zero when viewed from a place where the frequency is low to a place where it is one of the maxima of the entire graph. . Similarly, when the slope of a line segment connecting adjacent minimum values is viewed from a place where the frequency is low to a place where the frequency is high, the point where the slope changes from negative to positive or 0 is assumed to be one of the minimums in the graph.
【００４５】このようにして求めた極大値と極小値の差
が小さい極大値を除けば良い。まず、極大点を周波数の
低い順に見ていく。一つの極大点について左右の極小点
を選ぶ。次に、極大値と左右の極小値の差Ａ，Ｂを求め
る。さらに、極大値と左右の極小値との平均Ｃ，Ｄを求
める。Ａ／ＣとＢ／Ｄの最大値（ＭＡＸ）を求めＥとお
く。Ｅがある一定値よりも小の場合は極大とは見なさな
いことにする。最後の極大まで上記の手順を行う。これ
により、山と谷の差が小さい極大点は取り除くことがで
きる。It is only necessary to remove the maximum value having a small difference between the maximum value and the minimum value thus obtained. First, let's look at the maximum points in ascending order of frequency. Select left and right minimum points for one maximum point. Next, the differences A and B between the maximum value and the left and right minimum values are obtained. Further, the averages C and D of the maximum value and the left and right minimum values are obtained. The maximum value (MAX) of A / C and B / D is calculated and set as E. If E is smaller than a certain value, it is not considered as the maximum. Repeat the above procedure until the final maximum. As a result, it is possible to remove the maximum point where the difference between the peak and the valley is small.
【００４６】第３の実施の形態（請求項７）
次に図４に本発明の第３の実施の形態を説明するフロー
チャートを示す。まず複数の紙片の長さからヤング率を
出す場合を考える。Ｅをヤング率とすると次式が成り立
つ。Third Embodiment (Claim 7) Next, FIG. 4 shows a flow chart for explaining the third embodiment of the present invention. First, consider the case where the Young's modulus is obtained from the lengths of a plurality of paper pieces. When E is Young's modulus, the following equation holds.
【００４７】[0047]
【数１】 [Equation 1]
【００４８】ここでνは共振周波数、Lは試料の長さ、
ρは密度、ｄは紙の厚さ、ｍは振動モードに関する定数
でｍ０＝１．８７５、ｍ１＝４．６９４、ｍ２＝７．８
５５、ｍ３＝１０．９９５などである。ｍ０が１次、ｍ
１が２次、ｍ２が３次、ｍ３が４次など共振周波数に対
応する。Where ν is the resonance frequency, L is the sample length,
ρ is the density, d is the thickness of the paper, m is a constant relating to the vibration mode, m 0 = 1.875, m 1 = 4.694, m 2 = 7.8.
55, m 3 = 1.995, etc. m 0 is primary, m
1 corresponds to the second order, m 2 corresponds to the third order, m 3 corresponds to the fourth order, etc.
【００４９】ここで一つの紙片の長さからヤング率Eを
求めると、他のLの場合、Eが同じであるとすると、共振
周波数νを逆算することができる。よって、逆算された
周波数を中心にある範囲内の周波数における伝達関数の
最大値を極大値にとると、容易に極大値を求めることが
できる。Here, when the Young's modulus E is obtained from the length of one piece of paper, the resonance frequency ν can be calculated backward when E is the same for other L. Therefore, if the maximum value of the transfer function at a frequency within a range centered on the back-calculated frequency is set to the maximum value, the maximum value can be easily obtained.
【００５０】尚、実際には、Eを最終的に求めるには、
（１）式から、Eが一定であるとするとνL２／ｍ２は一
定となる。そこで縦軸にν、横軸にｍ２／L２をとる
と、ほぼ直線上に並ぶことになる。その直線の傾きをｔ
ａｎθとすると、（１）式から、Actually, in order to finally obtain E,
From Eq. (1), if E is constant, then ν L 2 / m 2 is constant. Therefore, when ν is taken on the vertical axis and m 2 / L 2 is taken on the horizontal axis, they are arranged on a substantially straight line. The slope of that straight line is t
If an θ, then from equation (1),
【００５１】[0051]
【数２】 によってヤング率を求めることができる。[Equation 2] The Young's modulus can be calculated by
【００５２】第４の実施の形態（請求項８）
次に図５に、本発明の第４の実施の形態を説明するフロ
ーチャートを示す。まず、第３の実施の形態の説明の後
半で述べた方法によって、近似直線の傾きからヤング率
を算出する。Fourth Embodiment (Claim 8) Next, FIG. 5 shows a flow chart for explaining the fourth embodiment of the present invention. First, the Young's modulus is calculated from the slope of the approximate straight line by the method described in the latter half of the description of the third embodiment.
【００５３】近似直線の例を図６に示す。図６に示すよ
うに、右から２番目の点は近似直線から大きくずれてい
る。このような点を入れて近似直線の傾きを計算すると
全体の精度が落ちてしまうため、このような点は除いて
近似直線の傾きを計算する必要がある。そこで、図５に
示すフローチャートのように、近似直線からある一定量
以上、離れた共振周波数を排除して、改めて近似直線を
引きなおし、グラフの傾きからヤング率を算出すれば、
良好な精度を維持しながらヤング率を求めることができ
る。An example of the approximate straight line is shown in FIG. As shown in FIG. 6, the second point from the right deviates greatly from the approximate straight line. If the slope of the approximate straight line is calculated by including such points, the accuracy of the whole is degraded. Therefore, it is necessary to calculate the slope of the approximate straight line excluding such points. Therefore, as shown in the flowchart of FIG. 5, if the resonance frequency that is separated from the approximate straight line by a certain amount or more is eliminated, the approximate straight line is drawn again, and the Young's modulus is calculated from the slope of the graph,
It is possible to obtain the Young's modulus while maintaining good accuracy.
【００５４】第５乃至８の実施の形態（それぞれ請求項
１乃至４に対応）
図７に、上述した第１乃至第４の実施の形態（方法）を
用いる処理装置の構成図の一例を示す。構成図は請求項
１乃至４に対応する。Fifth to Eighth Embodiments (Corresponding to Claims 1 to 4 respectively) FIG. 7 shows an example of a block diagram of a processing apparatus using the above-described first to fourth embodiments (methods). . The block diagram corresponds to claims 1 to 4.
【００５５】図７において、入力部は、情報を入力する
ためのものであり、データやコマンドキー入力するため
のキーボード、タブレット、マウスなどである。ＣＰＵ
は各種処理のための演算、論理判定などを行い、バスに
接続された各構成要素を制御する。複数のＣＰＵを一つ
の装置において保有してもよい。また、Ｆｏｒｔｒａｎ
やＣによる開発環境では、並列処理専用コマンドを備え
たもの、自動的に並列処理を最適化して行えるような実
行ファイル作成機能を有するものもある。In FIG. 7, the input unit is for inputting information, and is a keyboard, tablet, mouse or the like for inputting data or command keys. CPU
Performs arithmetic operations for various kinds of processing, logical judgments, etc., and controls each component connected to the bus. A plurality of CPUs may be held in one device. Also, Fortran
Some of the development environments based on C and C have a command dedicated to parallel processing and some have an execution file creation function that can automatically optimize parallel processing.
【００５６】プログラムメモリは本発明で前述した処理
手順を含むプログラムを格納するメモリである。プログ
ラムメモリはＲＯＭであってもよいし、外部記憶装置な
どからプログラムがロードされるＲＡＭであってもよ
い。データメモリは各種処理で生じたデータを格納す
る。外部記憶装置は、計算機によって得られたデータが
保存されている。上記の各構成要素は、常法により構成
可能である。The program memory is a memory for storing a program including the processing procedure described above in the present invention. The program memory may be a ROM or a RAM loaded with a program from an external storage device or the like. The data memory stores data generated by various processes. The external storage device stores data obtained by the computer. Each of the above components can be constructed by a conventional method.
【００５７】第９乃至１２の実施の形態（請求項９乃至
１２に対応）
第９乃至１２の実施の形態は、請求項５乃至８に記載の
方法を実行するプログラムが格納された記録媒体であ
る。そのようなデータ処理方法をプログラム化して格納
する記録媒体は、常法によって作成、あるいは複製可能
である。Ninth to twelfth embodiments (corresponding to claims 9 to 12) The ninth to twelfth embodiments are recording media in which a program for executing the method according to claims 5 to 8 is stored. is there. A recording medium for storing such a data processing method by programming can be created or copied by a conventional method.
【００５８】第１３の実施の形態（請求項１３）
まず、図８において振動リード法によって得られる横軸
を周波数、縦軸を伝達関数をとったグラフを示す。従来
は、このグラフの極大付近の左右の点に印を付けること
によって、半自動的に、極大を求めることができた。Thirteenth Embodiment (Claim 13) First, FIG. 8 shows a graph in which the horizontal axis represents the frequency and the vertical axis represents the transfer function obtained by the vibration lead method. Conventionally, it was possible to find the maximum semi-automatically by marking the left and right points near the maximum of this graph.
【００５９】図９において、図８に対応する横軸に周波
数、縦軸にクロス位相の１階微分をとったグラフを示
す。従来、図８と図９とを見比べながら、図８の極大と
図９で大きくクロス位相の１階微分が変化するところを
みて、図８の極大を求めていた。In FIG. 9, a graph corresponding to FIG. 8 is shown with the horizontal axis representing frequency and the vertical axis representing the first derivative of the cross phase. Conventionally, by comparing FIG. 8 and FIG. 9, the maximum of FIG. 8 has been obtained by observing that the first derivative of the cross phase greatly changes between the maximum of FIG. 8 and FIG. 9.
【００６０】図１０に第１３の実施の形態を説明するた
めのフローチャートを示す（請求項１３に対応）。FIG. 10 shows a flow chart for explaining the thirteenth embodiment (corresponding to claim 13).
【００６１】まず第１次の共振周波数を求める。図８に
見られるような細かいノイズの影響を取り除くために、
既に詳述したノイズ除去グラフ作成工程をおこなう。そ
して、既に詳述したと同様に、一定長さの区間に分割し
てその区間内において前記グラフの最大値を求める。First, the first-order resonance frequency is obtained. In order to remove the effect of fine noise as seen in FIG.
The noise removal graph creation process described in detail above is performed. Then, in the same manner as described in detail above, the segment is divided into sections of constant length, and the maximum value of the graph is obtained within the section.
【００６２】次に、隣接する最大値を結ぶ線分の傾きを
計算する。周波数の低いところから、高いところへ順に
見ていった時、傾きの符号が正から、負または０に変化
する点を極大点の候補とみなす。極大点の候補のうち、
横軸の周波数が最小のものを第１次の共振周波数とす
る。Next, the slope of the line segment connecting the adjacent maximum values is calculated. When the frequency is viewed from a low frequency area to a high frequency area, the point where the sign of the slope changes from positive to negative or 0 is regarded as a maximum point candidate. Of the candidates for maximum points,
The one with the minimum frequency on the horizontal axis is the primary resonance frequency.
【００６３】次に第１次の共振周波数からヤング率を計
算する。Ｅをヤング率とすると次式が成り立つ。Next, Young's modulus is calculated from the first-order resonance frequency. When E is Young's modulus, the following equation holds.
【００６４】[0064]
【数３】 [Equation 3]
【００６５】ここでνは共振周波数、Lは試料の長さ、
ρは密度、ｄは紙の厚さ、ｍは振動モードに関する定数
でｍ０＝１．８７５、ｍ１＝４．６９４、ｍ２＝７．８
５５、ｍ３＝１０．９９５などである。ｍ０が１次、ｍ
１が２次、ｍ２が３次、ｍ３が４次など共振周波数に対
応する。Where ν is the resonance frequency, L is the sample length,
ρ is the density, d is the thickness of the paper, m is a constant relating to the vibration mode, m 0 = 1.875, m 1 = 4.694, m 2 = 7.8.
55, m 3 = 1.995, etc. m 0 is primary, m
1 corresponds to the second order, m 2 corresponds to the third order, m 3 corresponds to the fourth order, etc.
【００６６】ここでｍ＝ｍ０としてνに第１次の共振周
波数を入れると、Ｅを求めることができる。（ただしこ
こでは仮に求めるだけで最終的なＥの出し方は後で述べ
る。）Ｅを算出したら、今度はｍ＝ｍ１、ｍ＝ｍ２、…
…として、逆に第２次の共振周波数、第３次の共振周波
数を求めることができる。If m = m 0 and the first-order resonance frequency is entered in ν, E can be obtained. (However, here, how to obtain the final E by merely obtaining it will be described later.) After calculating E, this time, m = m 1 , m = m 2 , ...
On the contrary, the second-order resonance frequency and the third-order resonance frequency can be obtained on the contrary.
【００６７】次に逆算された周波数を中心にある範囲内
の周波数における伝達関数の最大値を極大値とする。Next, the maximum value of the transfer function at a frequency within a range centered on the back-calculated frequency is set as the maximum value.
【００６８】これにより、伝達関数の極大を自動的に求
めることができた。尚、実際には、Ｅを最終的に求める
には、（１）式から、Ｅが一定であるとするとνＬ２／
ｍ２ は一定となる。そこで縦軸にν、横軸にｍ２／Ｌ
２をとると、ほぼ直線上に並ぶことになる。その直線の
傾きをｔａｎθとすると、（１）式から、As a result, the maximum of the transfer function could be automatically obtained. Actually, in order to finally obtain E, from Eq. (1), if E is constant, νL 2 /
m 2 is constant. Therefore, the vertical axis is ν and the horizontal axis is m 2 / L
If you take 2 , the lines will be almost straight. If the inclination of the straight line is tan θ, from equation (1),
【００６９】[0069]
【数４】 でヤング率を求めることができる。[Equation 4] The Young's modulus can be calculated with.
【００７０】第１４の実施の形態（請求項１４）
図１１に本発明の第１４の実施の形態を説明するフロー
チャートを示す。Fourteenth Embodiment (Claim 14) FIG. 11 shows a flowchart for explaining the fourteenth embodiment of the present invention.
【００７１】まず、縦軸に伝達関数、横軸に周波数のグ
ラフを作成する。それから、極大の候補を求める。これ
は、以下のようにする。図８に見られるような細かいノ
イズの影響を取り除く。既に詳述したノイズ除去グラフ
作成工程を行う。次に、既に詳述したと同様に各区間に
おける最大値を求め、隣接する最大値を結ぶ線分の傾き
を計算する。周波数の低いところから、高いところへ順
に見ていった時、傾きの符号が正から、負または０に変
化する点を極大の候補とみなす。First, a graph of transfer function is plotted on the vertical axis and frequency is plotted on the horizontal axis. Then seek the maximal candidate. This is done as follows. Eliminate the effects of fine noise as seen in FIG. The noise removal graph creation step already described in detail is performed. Next, as described above in detail, the maximum value in each section is obtained, and the slope of the line segment connecting the adjacent maximum values is calculated. A point where the sign of the slope changes from positive to negative or zero when viewed from a place where the frequency is low to a place where the frequency is high is regarded as a maximum candidate.
【００７２】次に縦軸にクロス位相の１階微分、横軸に
周波数をとったグラフを作成する。ここで、極大点の候
補に対応する周波数の前後のある一定の範囲内におい
て、クロス位相の１階微分がある一定範囲以上変化して
いる時、その極大の候補を極大と判定する。なぜなら、
クロス位相１回微分の性質上、その値が大であれば極大
値である可能性が大であるからである。本実施の形態で
は、クロス位相１回微分は、ＰＣ ＦＦＴアナライザー
ＯＲ２５（東陽テクニカ製高速フーリエ変換器）を使用
して、通常のウィンドウズ（登録商標）準拠のＰＣで解
析した。Next, a graph is created with the vertical axis representing the first derivative of the cross phase and the horizontal axis representing the frequency. Here, within a certain range before and after the frequency corresponding to the maximum point candidate, when the first derivative of the cross phase changes by a certain range or more, the maximum candidate is determined to be the maximum. Because
This is because, due to the property of the first cross-phase differentiation, if the value is large, there is a high possibility that it will be the maximum value. In the present embodiment, the cross-phase one-time differential is analyzed by a normal PC compatible with Windows (registered trademark) using a PC FFT analyzer OR25 (Fast Fourier Transform manufactured by Toyo Technica).
【００７３】この判別によって、自動的に極大点を求め
ることができる。尚、図９に示すような、ねじれのモー
ド（図８において左から２番目（Ｚ）の極大点の候補）
でも、クロス位相の１階微分は多少変化するので、それ
以上の変化があるときのみを極大とするように注意す
る。By this discrimination, the maximum point can be automatically obtained. A twist mode as shown in FIG. 9 (the second (Z) maximum point candidate from the left in FIG. 8)
However, the first derivative of the cross phase changes a little, so be careful to maximize only when there is a further change.
【００７４】また、第１３の実施の形態と第１４の実施
の形態を組み合わせて使用しても良い。具体的には、第
１の共振周波数を求め、ヤング率を算出し、第２、第３
……の共振周波数を逆算し、その周波数の前後のある範
囲内でクロス位相の１階微分がある一定の値以上変化し
ているとき、その周波数を共振周波数とすれば良い。こ
れによって、さらに、正確に極大値を求めることができ
る。Further, the thirteenth embodiment and the fourteenth embodiment may be combined and used. Specifically, the first resonance frequency is obtained, the Young's modulus is calculated, and the second and third
The resonance frequency of ... is calculated backward, and when the first derivative of the cross phase changes by a certain value or more within a certain range before and after that frequency, that frequency may be used as the resonance frequency. With this, the maximum value can be obtained more accurately.
【００７５】[0075]
【発明の効果】請求項５に係る発明によって、自動的に
極大値を求める事ができるので、ヤング率を用いるシミ
ュレーション計算においてスピードアップを行い、省力
化を達成することができる方法を提供できる。また、そ
の方法を用いた装置（請求項１に係る発明）、その方法
をプログラムとして有する記録媒体（請求項９に係る発
明）においては、同様の効果を有する装置と記録媒体と
を、それぞれ提供できる。According to the invention of claim 5, since the maximum value can be automatically obtained, it is possible to provide a method capable of speeding up the simulation calculation using Young's modulus and achieving labor saving. Further, in an apparatus using the method (invention of claim 1) and a recording medium having the method as a program (invention of claim 9), an apparatus and a recording medium having similar effects are provided. it can.
【００７６】請求項６に係る発明においては、ねじれモ
ードの極大値を除くことができるので、より正確なヤン
グ率を自動的に求める事ができ、シミュレーション計算
においてスピードアップ、省力化を達成できる方法を提
供できる。また、その方法を用いた装置（請求項２に係
る発明）、プログラムとして格納された記録媒体（請求
項１０に係る発明）においても同様の効果を有する装置
と記録媒体とをそれぞれ提供できる。In the invention according to claim 6, since the maximum value of the twisting mode can be eliminated, a more accurate Young's modulus can be automatically obtained, and a speedup and labor saving can be achieved in the simulation calculation. Can be provided. Further, an apparatus using the method (an invention according to claim 2) and a recording medium stored as a program (an invention according to claim 10) can also provide an apparatus and a recording medium having similar effects.
【００７７】請求項７に係る発明においては、決められ
た範囲内において最大値を求めるだけなので、ヤング率
を使ったシミュレーション計算において、さらにスピー
ドアップできる方法を提供できる。また、その方法を用
いた装置（請求項３に係る発明）、プログラムとして格
納した記録媒体（請求項１１に係る発明）においては、
同様の効果を有する装置と記録媒体とをそれぞれ提供で
きる。In the invention according to claim 7, since the maximum value is simply obtained within the determined range, it is possible to provide a method capable of further speeding up the simulation calculation using the Young's modulus. Further, in an apparatus using the method (the invention according to claim 3) and a recording medium stored as a program (the invention according to claim 11),
It is possible to provide an apparatus and a recording medium that have similar effects.
【００７８】請求項８に係る発明においては、近似直線
からはずれた共振周波数を除くので、さらに正確なヤン
グ率を求める方法を提供できる。また、その方法を用い
た装置（請求項４に係る発明）、そのプログラムを格納
した記録媒体（請求項１２に係る発明）においても、同
様の効果を有する装置と記録媒体とをそれぞれ提供でき
る。In the invention according to claim 8, since the resonance frequency deviated from the approximate straight line is excluded, it is possible to provide a more accurate method of obtaining the Young's modulus. Also, in an apparatus using the method (invention of claim 4) and a recording medium storing the program (invention of claim 12), an apparatus and a recording medium having the same effect can be provided.
【００７９】請求項１３に係る発明においては、自動的
に極大を求めることができるので、シミュレーション計
算において、更にスピードアップを行い、省力化を達成
する方法を提供できる。According to the thirteenth aspect of the present invention, since the maximum value can be automatically obtained, it is possible to provide a method for further speeding up the simulation calculation and achieving labor saving.
【００８０】請求項１４に係る発明においては、自動的
に極大を求めることができ、追加的にクロス位相の１階
微分を参照しているので、シミュレーション計算におい
て、さらに正確に、しかもスピードアップを行い、省力
化を達成する方法を提供できる。According to the fourteenth aspect of the present invention, the maximum value can be automatically obtained and the first derivative of the cross phase is additionally referred, so that the simulation calculation can be performed more accurately and speed up. It is possible to provide a method for achieving labor saving.
【図１】振動リード法によって得られる伝達関数のグラ
フである。FIG. 1 is a graph of a transfer function obtained by a vibration reed method.
【図２】第１の実施の形態を説明するためのフローチャ
ートである。FIG. 2 is a flowchart for explaining the first embodiment.
【図３】第２の実施の形態を説明するためのフローチャ
ートである。FIG. 3 is a flowchart for explaining a second embodiment.
【図４】第３の実施の形態を説明するためのフローチャ
ートである。FIG. 4 is a flowchart for explaining a third embodiment.
【図５】第４の実施の形態を説明するためのフローチャ
ートである。FIG. 5 is a flowchart for explaining a fourth embodiment.
【図６】第４の実施の形態における共鳴周波数を近似す
る直線のグラフである。FIG. 6 is a graph of a straight line approximating the resonance frequency in the fourth embodiment.
【図７】第５乃至８の実施形態における処理装置の一例
を示す構成図である。FIG. 7 is a configuration diagram illustrating an example of a processing device according to fifth to eighth embodiments.
【図８】振動リード法によって得られる伝達関数のグラ
フである。FIG. 8 is a graph of a transfer function obtained by the vibration reed method.
【図９】伝達関数のクロス位相１回微分を示したグラフ
である。FIG. 9 is a graph showing a cross phase once derivative of a transfer function.
【図１０】第１３の実施の形態を説明するためのフロー
チャートである。FIG. 10 is a flowchart for explaining a thirteenth embodiment.
【図１１】第１４の実施の形態を説明するための、フロ
ーチャートである。FIG. 11 is a flowchart for explaining the fourteenth embodiment.
Claims (14)
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２区間における
前記グラフにおけるの最大値を求めて、隣接する前記最
大値同士を結ぶ線分の傾きが、正から、負又は０に変わ
る点を該グラフ全体の極大点の一つであると判定するこ
とを特徴とするデータ処理装置。1. A method of obtaining a Young's modulus of a paper by a vibration lead method using a resonance frequency of a paper and using the resonance frequency, wherein a frequency is plotted on a horizontal axis, and the frequency region has a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Dividing into a second closed section having a length, obtaining a maximum value of in the graph in each of the second sections, and a slope of a line segment connecting the adjacent maximum values changes from positive to negative or zero. Is determined as one of the maximum points of the entire graph.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２閉区間におけ
る前記グラフにおける最大値と最小値とを求めて、隣接
する前記最大値同士を結ぶ線分の傾きが正から、負また
は０に変わる点をそのグラフ全体の極大点の候補の１つ
であると判定し、かつ、 隣接する前記最小値同士を結ぶ線分の傾きが、負から、
正または０に変わる点をそのグラフ全体の極小点の一つ
であると判定し、 更に、ある一つの前記極大点について、前記極大値を与
える周波数の左右で最も近い位置にある極小値を求め、
前記極大値と該左右それぞれの極小値との差Ａ，Ｂを求
め、また、前記極大値と該左右それぞれの極小値との平
均値Ｃ、Ｄを求め、 Ａ／Ｃ、Ｂ／Ｄのいずれか大なる方がある一定値よりも
小さい時、前記極大点を共振による極大点とは判定しな
いことを特徴とするデータ処理装置。2. A method for obtaining a Young's modulus of a paper by a vibration lead method using a resonance frequency of a paper, the frequency being plotted on the horizontal axis, and the frequency region having a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Dividing into a second closed section having a length, the maximum value and the minimum value in the graph in each of the second closed sections are obtained, and the slope of the line segment connecting the adjacent maximum values is changed from positive to negative or 0. Is determined as one of the candidates for the maximum point of the entire graph, and the slope of the line segment connecting the minimum values adjacent to each other is negative,
The point that changes to positive or 0 is determined to be one of the minimum points of the entire graph, and for one certain maximum point, the minimum value at the closest position on the left and right of the frequency giving the maximum value is obtained. ,
The difference A, B between the maximum value and the respective left and right minimum values is obtained, and the average value C, D of the maximum value and the left and right respective minimum values is calculated, and either A / C or B / D is obtained. A data processing device, wherein the maximum point is not determined as a maximum point due to resonance when the larger one is smaller than a certain value.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さから、該紙のヤング率を求める
場合、 第１の長さについて該紙のヤング率を求めた後、さら
に、第２の長さの該紙のヤング率を求める時、既知のヤ
ング率から共振周波数を逆算して得られた周波数を中心
値として含むある一定周波数の範囲内において最大値が
存在する場合、該最大値を極大値と判定することを特徴
とするデータ処理装置。3. A method of obtaining a Young's modulus of a paper by a vibration lead method using the resonance frequency of a paper, wherein the Young's modulus of the paper is obtained from the lengths of a plurality of paper pieces. After obtaining the Young's modulus of the paper for the length of, and further obtaining the Young's modulus of the paper of the second length, the frequency obtained by back-calculating the resonance frequency from the known Young's modulus is taken as the central value. A data processing device, characterized in that, when a maximum value exists within a certain constant frequency range including the maximum value, the maximum value is determined as a maximum value.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さそれぞれに対して、一つまたは
複数の共振周波数から求められる近似直線の傾きから算
出される一つのヤング率の値に対して、該近似直線から
ある一定以上離れた共振周波数を除いて、ヤング率を算
出することを特徴とするデータ処理装置。4. A method of obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency for each length of a plurality of paper pieces from one or a plurality of resonance frequencies. A data processing device, wherein a Young's modulus is calculated with respect to one Young's modulus value calculated from the obtained slope of the approximate straight line, by removing a resonance frequency that is apart from the approximate straight line by a certain amount or more.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２区間における
前記グラフの最大値を求めて、隣接する前記最大値同士
を結ぶ線分の傾きが、正から、負又は０に変わる点を該
グラフ全体の極大点の一つであると判定することを特徴
とするデータ処理方法。5. A method for obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency, wherein a frequency is plotted on a horizontal axis, and the frequency region has a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Divide into second closed sections having a length, find the maximum value of the graph in each of the second sections, and determine the point at which the slope of the line segment connecting the adjacent maximum values changes from positive to negative or zero. A data processing method, characterized in that it is judged to be one of the maximum points of the entire graph.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２閉区間におけ
る前記グラフの最大値と最小値とを求めて、隣接する前
記最大値同士を結ぶ線分の傾きが正から、負または０に
変わる点をそのグラフ全体の極大点の候補の１つである
と判定し、かつ、 隣接する前記最小値同士を結ぶ線分の傾きが、負から、
正または０に変わる点をそのグラフ全体の極小点の一つ
であると判定し、 更に、ある一つの前記極大について、前記極大値を与え
る周波数の左右で最も近い位置にある極小値を求め、前
記極大値と該左右それぞれの極小値との差Ａ，Ｂを求
め、また、前記極大値と該左右それぞれの極小値との平
均値Ｃ、Ｄを求め、 Ａ／Ｃ、Ｂ／Ｄのいずれか大なる方がある一定値よりも
小さい時、前記極大を共振による極大とは判定しないこ
とを特徴とするデータ処理方法。6. A method of obtaining a resonance frequency of a paper and using the resonance frequency of the paper to obtain a Young's modulus of the paper by a vibration lead method, wherein a frequency is plotted on a horizontal axis and the frequency region is a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Dividing into a second closed section having a length, obtaining the maximum value and the minimum value of the graph in each of the second closed sections, the slope of the line segment connecting the adjacent maximum values is from positive to negative or 0. Is determined as one of the candidates for the maximum point of the entire graph, and the slope of the line segment connecting the minimum values adjacent to each other is negative,
The point that changes to positive or 0 is determined to be one of the minimum points of the entire graph, and further, for one certain maximum value, the minimum value at the closest position on the left and right of the frequency giving the maximum value is obtained, The difference A, B between the maximum value and the respective left and right minimum values is obtained, and the average value C, D of the maximum value and the left and right respective minimum values is calculated, and either A / C or B / D is obtained. When the larger one is smaller than a certain value, the maximum is not judged to be the maximum due to resonance.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さから、該紙のヤング率を求める
場合、 第１の長さについて該紙のヤング率を求めた後、さら
に、第２の長さの該紙のヤング率を求める時、既知のヤ
ング率から共振周波数を逆算して得られた周波数を中心
に含むある一定周波数の範囲内において最大値が存在す
る場合、該最大値を極大値と判定することを特徴とする
データ処理方法。7. A method of obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of the paper and using the resonance frequency, when the Young's modulus of the paper is obtained from the lengths of a plurality of pieces of paper. After obtaining the Young's modulus of the paper with respect to the length of, and further determining the Young's modulus of the paper of the second length, the resonance frequency is calculated backward from the known Young's modulus, and the frequency is included in the center. A data processing method characterized in that, when a maximum value exists within a certain frequency range, the maximum value is determined as a maximum value.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さそれぞれに対して、一つまたは
複数の共振周波数から求められる近似直線の傾きから算
出される一つのヤング率の値に対して、該近似直線から
ある一定以上離れた共振周波数を除いて、ヤング率を算
出することを特徴とするデータ処理方法。8. A method for obtaining Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of paper and using the resonance frequency for each length of a plurality of paper pieces from one or a plurality of resonance frequencies. A data processing method, wherein Young's modulus is calculated for one value of Young's modulus calculated from the obtained slope of the approximate straight line by removing the resonance frequency which is apart from the approximate straight line by a certain amount or more.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２区間における
前記グラフの最大値を求めて、隣接する前記最大値同士
を結ぶ線分の傾きが、正から、負又は０に変わる点を該
グラフ全体の極大点の一つであると判定することを特徴
とするコンピュータが読み込み可能なプログラムを格納
した記録媒体。9. A method of obtaining a Young's modulus of a paper by a vibration lead method using a resonance frequency of a paper and using the resonance frequency, wherein a frequency is plotted on a horizontal axis, and the frequency region has a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Divide into second closed sections having a length, find the maximum value of the graph in each of the second sections, and determine the point at which the slope of the line segment connecting the adjacent maximum values changes from positive to negative or zero. A recording medium storing a computer-readable program, which is characterized as being one of the maximum points of the entire graph.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２閉区間におけ
る前記グラフの最大値と最小値とを求めて、隣接する前
記最大値同士を結ぶ線分の傾きが正から、負または０に
変わる点をそのグラフ全体の極大点の候補の１つである
と判定し、かつ、 隣接する前記最小値同士を結ぶ線分の傾きが、負から、
正または０に変わる点をそのグラフ全体の極小点の一つ
であると判定し、 更に、ある一つの前記極大について、前記極大値を与え
る周波数の左右で最も近い位置にある極小値を求め、前
記極大値と該左右それぞれの極小値との差Ａ，Ｂを求
め、また、前記極大値と該左右それぞれの極小値との平
均値Ｃ、Ｄを求め、 Ａ／Ｃ、Ｂ／Ｄのいずれか大なる方がある一定値よりも
小さい時、前記極大を共振による極大とは判定しないこ
とを特徴とするコンピュータが読み込み可能なプログラ
ムを格納した記録媒体。10. A method for obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency, wherein the frequency is plotted on the horizontal axis, and the frequency region has a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. It is divided into second closed sections having a length, the maximum value and the minimum value of the graph in each of the second closed sections are obtained, and the slope of the line segment connecting the adjacent maximum values is changed from positive to negative or 0. Is determined as one of the candidates for the maximum point of the entire graph, and the slope of the line segment connecting the minimum values adjacent to each other is negative,
The point that changes to positive or 0 is determined to be one of the minimum points of the entire graph, and further, for one certain maximum value, the minimum value at the closest position on the left and right of the frequency giving the maximum value is obtained, The difference A, B between the maximum value and the left and right minimum values is calculated, and the average value C, D of the maximum value and the left and right minimum values is calculated, and either A / C or B / D is calculated. A recording medium storing a computer-readable program, wherein the maximum is not determined to be a maximum due to resonance when the larger one is smaller than a certain value.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さから、該紙のヤング率を求める
場合、 第１の長さについて該紙のヤング率を求めた後、さら
に、第２の長さの該紙のヤング率を求める時、既知のヤ
ング率から共振周波数を逆算して得られた周波数を中心
に含むある一定周波数の範囲内において最大値が存在す
る場合、該最大値を極大値と判定することを特徴とする
コンピュータが読み込み可能なプログラムを格納した記
録媒体。11. A method of obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of the paper and using the resonance frequency, when the Young's modulus of the paper is obtained from the lengths of a plurality of paper pieces. After obtaining the Young's modulus of the paper with respect to the length of, and further determining the Young's modulus of the paper of the second length, the resonance frequency is calculated backward from the known Young's modulus, and the frequency is included in the center. A recording medium storing a computer-readable program, wherein a maximum value is determined to be a maximum value when the maximum value exists within a certain frequency range.
て、振動リード法により紙のヤング率を求める方法にお
いて、複数の紙片の長さそれぞれに対して、一つまたは
複数の共振周波数から求められる近似直線の傾きから算
出される一つのヤング率の値に対して、該近似直線から
ある一定以上離れた共振周波数を除いて、ヤング率を算
出することを特徴とするコンピュータが読み込み可能な
プログラムを格納した記録媒体。12. A method of obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency from one or a plurality of resonance frequencies for each length of a plurality of paper pieces. A computer-readable method characterized in that the Young's modulus is calculated by excluding a resonance frequency that is apart from the approximate straight line by a certain value with respect to one Young's modulus value calculated from the obtained slope of the approximate straight line. A recording medium that stores a program.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、前記横軸において該周波数領域を
第１の一定長さに分割した複数の第１閉区間を設け、前
記各第１閉区間における伝達関数の平均値を縦軸にとっ
たグラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２区間における
前記グラフの最大値を求めて、隣接する前記最大値同士
を結ぶ線分の傾きが、正から、負又は０に変わる点を該
グラフ全体の極大点の候補値であると判定する方法にお
いて、 第１次共振による１次極大値の候補をもとにして第２次
共振による２次極大値、３次共振による３次極大値な
ど、該１次以外の共振による極大値を予想し、前記１次
以外の共振による極大値を与えると予想された周波数を
含むある範囲において極大点の候補が存在する場合、該
極大点による極大値の候補を、該２次、３次などの極大
値であると判定することを特徴とするデータ処理方法。13. A method of obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency, wherein a frequency is plotted on a horizontal axis, and the frequency region has a first constant length on the horizontal axis. A plurality of first closed sections divided into two, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is longer than the first constant length. Divide into second closed sections having a length, find the maximum value of the graph in each of the second sections, and determine the point at which the slope of the line segment connecting the adjacent maximum values changes from positive to negative or zero. In the method of determining the candidate value of the maximum point of the whole graph, the secondary maximum value of the secondary resonance and the third maximum value of the third resonance are determined based on the candidate of the primary maximum value of the first resonance. Predict the maximum value due to resonance other than the first order, such as , If there is a candidate for a maximum point in a certain range including a frequency expected to give a maximum value due to resonance other than the first order, the candidate for the maximum point due to the maximum point is converted to the second, third, etc. A data processing method characterized by determining that the value is a value.
て、振動リード法により紙のヤング率を求める方法にお
いて、 周波数を横軸にとり、該横軸において周波数領域を第１
の一定長さに分割した複数の第１閉区間を設け、前記各
第１閉区間における伝達関数の平均値を縦軸にとったグ
ラフを作り、 前記周波数領域を前記第１の一定長さよりも長い他の長
さを有する第２閉区間に分け、各前記第２閉区間にける
前記グラフの最大値を求め、隣接する前記最大値同士を
結ぶ線分の傾きが正から、負又は０に変わる点をそのグ
ラフ全体の極大点の候補と判定する方法において、 前記伝達関数をもとにしたグラフに基づいてクロス位相
１階微分したクロス位相１回微分のグラフを作り、前記
極大の候補を与える周波数を含む一定の範囲内におい
て、該クロス位相１階微分値が一定量以上変化している
時、前記極大点の候補を極大値と判定することを特徴と
するデータ処理方法。14. A method for obtaining a Young's modulus of a paper by a vibration lead method by obtaining a resonance frequency of a paper and using the resonance frequency, wherein a frequency is plotted on a horizontal axis and a frequency domain is set to a first frequency region on the horizontal axis.
A plurality of first closed sections divided into constant lengths, and a graph in which the average value of the transfer function in each of the first closed sections is taken on the vertical axis, and the frequency domain is set to be longer than the first constant length. Divide into second closed sections having other long lengths, find the maximum value of the graph in each of the second closed sections, and change the slope of the line segment connecting the adjacent maximum values from positive to negative or zero. In a method of determining a changing point as a candidate for a maximum point of the entire graph, a graph of a cross phase one-time derivative obtained by differentiating the cross phase by the first order is created based on the graph based on the transfer function, and the candidate for the maximum is obtained. A data processing method, wherein the candidate of the maximum point is determined to be a maximum value when the cross phase first-order differential value changes by a certain amount or more within a certain range including a given frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001203093A JP2003014702A (en) | 2001-07-04 | 2001-07-04 | Apparatus and method for measurement of young's modulus of paper |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001203093A JP2003014702A (en) | 2001-07-04 | 2001-07-04 | Apparatus and method for measurement of young's modulus of paper |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2003014702A true JP2003014702A (en) | 2003-01-15 |
Family
ID=19039787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001203093A Pending JP2003014702A (en) | 2001-07-04 | 2001-07-04 | Apparatus and method for measurement of young's modulus of paper |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2003014702A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100343661C (en) * | 2005-06-21 | 2007-10-17 | 南京大学 | Method for measuring complex Young's modulus in fluid dynamic spectrum |
| JP2010223925A (en) * | 2009-03-25 | 2010-10-07 | Toshiba Corp | Elastic modulus measuring device of paper sheet and paper sheet processor |
| CN107655756A (en) * | 2017-09-18 | 2018-02-02 | 天津商业大学 | A kind of word bomb Modulus Test Method based on speckle printing |
-
2001
- 2001-07-04 JP JP2001203093A patent/JP2003014702A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100343661C (en) * | 2005-06-21 | 2007-10-17 | 南京大学 | Method for measuring complex Young's modulus in fluid dynamic spectrum |
| JP2010223925A (en) * | 2009-03-25 | 2010-10-07 | Toshiba Corp | Elastic modulus measuring device of paper sheet and paper sheet processor |
| CN107655756A (en) * | 2017-09-18 | 2018-02-02 | 天津商业大学 | A kind of word bomb Modulus Test Method based on speckle printing |
| CN107655756B (en) * | 2017-09-18 | 2020-01-10 | 天津商业大学 | Paper elastic modulus testing method based on speckle printing |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103679190B (en) | Sorter, classification method and electronic equipment | |
| Martin | Computational improvements to estimating kriging metamodel parameters | |
| CN116432091B (en) | Equipment fault diagnosis method based on small sample, construction method and device of model | |
| Anaissi et al. | Regularized tensor learning with adaptive one-class support vector machines | |
| CN115359846A (en) | Batch correction method and device for group data, storage medium and electronic equipment | |
| CN113469151B (en) | Method, device, equipment and medium for screening sensor in industrial manufacturing | |
| WO2019235603A1 (en) | Relationship analysis device, relationship analysis method, and recording medium | |
| CN113646714A (en) | Processing parameter setting method and device for production equipment and computer readable medium | |
| CN110598305B (en) | Sensitivity analysis method for comparing scanning simulation increment of circuit | |
| JP2003014702A (en) | Apparatus and method for measurement of young's modulus of paper | |
| Kleppe | Modified Cholesky Riemann manifold Hamiltonian Monte Carlo: exploiting sparsity for fast sampling of high-dimensional targets | |
| CN112347776B (en) | Medical data processing method and device, storage medium and electronic equipment | |
| JP4994199B2 (en) | Machine learning apparatus and machine learning method | |
| WO2022091639A1 (en) | Abnormality diagnosing model construction method, abnormality diagnosing method, abnormality diagnosing model construction device, and abnormality diagnosing device | |
| Su et al. | An online outlier detection method based on wavelet technique and robust RBF network | |
| US20130191071A1 (en) | System and method for automatic modal parameter extraction in structural dynamics analysis | |
| CN115222114A (en) | Automobile part assembly index value prediction method, terminal device and storage medium | |
| CN112505337B (en) | Data processing method for assisting analysis of sample | |
| Volz et al. | Surrogate-assisted partial order-based evolutionary optimisation | |
| CN113962499A (en) | Material quality detection and analysis method, device, equipment and storage medium | |
| CN113642209A (en) | Structure implantation fault response data acquisition and evaluation method based on digital twinning | |
| JP2940470B2 (en) | Structure analysis method and structure analysis device | |
| JP2004227279A (en) | Failure cause diagnosing method using mahalanobis distance and program | |
| CN112783772B (en) | Software reliability analysis method, system, medium, equipment, terminal and application | |
| CN113742849B (en) | Variable sensitivity analysis method and device for solid-liquid-like aircraft overall design |