JP2018179730A - Stress measurement device and stress measurement method - Google Patents

Stress measurement device and stress measurement method Download PDF

Info

Publication number
JP2018179730A
JP2018179730A JP2017079121A JP2017079121A JP2018179730A JP 2018179730 A JP2018179730 A JP 2018179730A JP 2017079121 A JP2017079121 A JP 2017079121A JP 2017079121 A JP2017079121 A JP 2017079121A JP 2018179730 A JP2018179730 A JP 2018179730A
Authority
JP
Japan
Prior art keywords
stress
temperature
load
temperature fluctuation
fluctuation amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2017079121A
Other languages
Japanese (ja)
Inventor
勇佐 大久保
Yusuke Okubo
勇佐 大久保
朋彦 春山
Tomohiko Haruyama
朋彦 春山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JTEKT Corp
Original Assignee
JTEKT Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by JTEKT Corp filed Critical JTEKT Corp
Priority to JP2017079121A priority Critical patent/JP2018179730A/en
Publication of JP2018179730A publication Critical patent/JP2018179730A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a stress measurement device and a stress measurement method capable of calculating more accurate temperature fluctuation quantity and searching for more accurate stress distribution in consideration of heat loss by a stress measurement device which calculates a stress distribution from a temperature fluctuation quantity at each position of a measurement subject based on temperature distribution and measuring a temperature distribution using infrared thermography, giving periodic load to a measurement subject.SOLUTION: A calculation device includes a temperature-fluctuation-quantity calculation part which calculates the temperature fluctuation quantity at each position of the surface of a measurement subject based on the temperature distribution by thermometry equipment, a heat-loss calculation part which calculates the thermal-diffusion loss information about the loss of the heat by diffusion of the heat on the surface of a measurement subject according to the cycle or frequency of the load given by the load giving equipment, an amount calculation part of the corrected temperature fluctuation quantity based on heat-diffusion loss information, and a stress calculation part which calculates the stress distribution of the surface of the measurement subject according to the load based on the amount of the corrected temperature change.SELECTED DRAWING: Figure 3

Description

本発明は、応力計測装置及び応力計測方法に関する。   The present invention relates to a stress measurement device and a stress measurement method.

近年、車両の燃費をより向上させるために、車両に搭載される種々の部品の軽量化や小型化が望まれている。例えば、高い剛性が必要とされて鉄等の金属で構成された比較的重量が有る部品(対象物)は、剛性が必要な部分の肉厚を残し、剛性をそれほど必要としない部分の肉厚を薄くすることで、軽量化に大きく貢献できる。そこで、対象物に負荷がかかった際、どの位置に、どのくらいの応力がかかっているか、を認識するために、種々の応力計測装置や応力計測方法が開示されている。   In recent years, in order to further improve the fuel consumption of a vehicle, weight reduction and miniaturization of various parts mounted on the vehicle are desired. For example, a relatively heavy part (object) requiring high rigidity and made of metal such as iron leaves a thickness of a portion requiring rigidity and a thickness of a portion not requiring much rigidity Can make a great contribution to weight reduction. Therefore, various stress measuring devices and stress measuring methods have been disclosed in order to recognize which position and how much stress is applied when an object is loaded.

例えば特許文献1には、等速ジョイント等の回転装置を計測対象物として、当該計測対象物に繰り返し荷重変動を生じさせ、計測対象物の温度分布を赤外線サーモグラフィで測定し、繰り返し荷重変動に応じた温度変動に基づいて応力分布を求める応力測定方法が開示されている。特許文献1では、計測対象物の初期温度を測定した後、繰り返し荷重変動を生じさせ、荷重変動の周期における第1の1周期分を含む第1周期と、第2の1周期分を含む第2周期を設定し、第1周期における第1平均温度、第2周期における第2平均温度を求めている。そして、初期温度、第1平均温度、第2平均温度から決定した温度特性を用いて、赤外線サーモグラフィで計測した温度分布を補正し、補正した温度分布から応力分布を求めている。   For example, in Patent Document 1, a rotating device such as a constant velocity joint is used as an object to be measured, load fluctuation is repeatedly generated on the object to be measured, temperature distribution of the object to be measured is measured by infrared thermography, There is disclosed a stress measurement method for obtaining a stress distribution based on temperature variations. In Patent Document 1, after measuring the initial temperature of the measurement object, repeated load fluctuation is caused, and a first cycle including a first one cycle in a cycle of the load fluctuation and a second cycle including a second cycle Two cycles are set, and the first average temperature in the first cycle and the second average temperature in the second cycle are obtained. Then, using the temperature characteristics determined from the initial temperature, the first average temperature, and the second average temperature, the temperature distribution measured by infrared thermography is corrected, and the stress distribution is determined from the corrected temperature distribution.

特開2012−103124号公報JP 2012-103124 A

特許文献1に記載の応力測定方法では、繰り返し荷重変動を生じさせている際(周期的な負荷を与えている際)、計測対象物には荷重変動に応じた(周期的な)温度変動が発生するが、荷重変動の開始時点から時間の経過に応じて計測対象物の全体の温度が徐々に上昇する。この温度上昇分を補正して、より正確な温度分布を求め、より正確な応力分布を求めている。なお、温度分布中の各位置の温度変動量から、各位置の応力を求める際、一般的には、下記の(式1)から求めている。なお(式1)中において、ΔTは温度変動量、kは計測対象物を構成している材料に固有の材料定数(熱弾性係数)、Tは材料の絶対温度、Δσは応力変動量である。
ΔT=−kTΔσ (式1)
In the stress measurement method described in Patent Document 1, when cyclic load fluctuation is caused (when cyclic load is applied), the measurement object has (periodical) temperature fluctuation corresponding to load fluctuation. Although it generate | occur | produces, the whole temperature of a measurement object will rise gradually according to progress of time from the start time of load fluctuation. The temperature rise is corrected to obtain a more accurate temperature distribution and a more accurate stress distribution. In addition, when calculating | requiring the stress of each position from the temperature fluctuation amount of each position in temperature distribution, it is generally calculated | required from the following (Formula 1). In Equation (1), ΔT is the amount of temperature change, k is the material constant (thermoelastic coefficient) specific to the material constituting the measurement object, T is the absolute temperature of the material, and Δσ is the amount of stress change .
ΔT = −kTΔσ (Expression 1)

上記の(式1)は、断熱状態のもとで成立する式である。しかし、実際には、計測対象物の各位置(応力がかかる各位置)を断熱状態とすることは不可能であり、周囲に熱が伝播して熱損失が発生している。従って、実際の温度変動量に対して、計測した温度変動量は熱損失の分だけ低くなっている。すなわち、実際の応力変動量に対して、求めた応力変動量は低くなってしまっている。   The above (Equation 1) is an equation that holds under the adiabatic condition. However, in actuality, it is impossible to make each position (each position to which stress is applied) of the object to be measured adiabatically, and heat propagates to the periphery to generate heat loss. Therefore, the measured temperature fluctuation is lower than the actual temperature fluctuation by the amount of heat loss. That is, the calculated stress fluctuation amount is lower than the actual stress fluctuation amount.

本発明は、このような点に鑑みて創案されたものであり、計測対象物に周期的な負荷を与えながら赤外線サーモグラフィを用いて温度分布を計測し、温度分布に基づいた計測対象物の各位置の温度変動量から応力分布を求める応力計測装置において、熱損失を考慮して、より正確な温度変動量を求め、より正確な応力分布を求めることができる、応力計測装置及び応力計測方法を提供することを目的とする。   The present invention has been made in view of such a point, and measures the temperature distribution using infrared thermography while periodically applying a load to the measurement object, and measures each of the measurement objects based on the temperature distribution. Stress measurement apparatus and stress measurement method capable of obtaining more accurate temperature variation and more accurate stress distribution in consideration of heat loss in a stress measurement apparatus for obtaining stress distribution from temperature variation of position Intended to be provided.

上記課題を解決するため、本発明の第1の発明は、応力の計測対象物に対して周期的に変動する負荷を与える負荷付与装置と、前記負荷付与装置にて前記計測対象物に前記負荷を与えながら赤外線サーモグラフィを用いて前記計測対象物の表面の温度分布を測定する温度測定装置と、前記負荷付与装置から与えた前記負荷と、前記温度測定装置による前記温度分布と、に基づいて前記計測対象物の表面の応力分布を求める演算装置と、を有する応力計測装置において、前記演算装置は、前記温度測定装置による前記温度分布に基づいて、前記計測対象物の表面の各位置の温度変動量を求める温度変動量演算部と、前記負荷付与装置から与えた前記負荷の周期または周波数に応じて、前記計測対象物の表面における熱の拡散による熱の損失に関する熱拡散損失情報を求める熱損失演算部と、前記温度変動量を前記熱拡散損失情報に基づいて補正した、補正温度変動量を求める補正温度変動量演算部と、前記補正温度変動量に基づいて、前記負荷に応じた前記計測対象物の表面の応力分布を求める応力演算部と、を有する、応力計測装置である。   In order to solve the above problems, according to a first aspect of the present invention, there is provided a load applying device for applying a load that periodically changes to a measurement object of stress, and the load application device applies the load to the measurement object. The temperature measurement device for measuring the temperature distribution of the surface of the measurement object using infrared thermography while giving the above, the load applied from the load application device, and the temperature distribution by the temperature measurement device. In a stress measuring device having a computing device for obtaining a stress distribution on the surface of the measurement object, the computing device is characterized in that the temperature fluctuation of each position of the surface of the measurement object is based on the temperature distribution by the temperature measuring device. Heat loss due to diffusion of heat on the surface of the object to be measured according to the temperature fluctuation computing unit for obtaining the amount and the period or frequency of the load applied from the load applying device; A heat loss calculation unit for obtaining heat diffusion loss information, a correction temperature change amount calculation unit for obtaining a correction temperature change amount obtained by correcting the temperature change amount based on the heat diffusion loss information, and And a stress calculation unit for obtaining a stress distribution on the surface of the measurement object according to the load.

次に、本発明の第2の発明は、上記第1の発明に係る応力計測装置であって、前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、前記演算装置は、前記熱損失演算部にて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、前記補正温度変動量演算部にて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、応力計測装置である。   A second invention of the present invention is the stress measurement device according to the first invention, wherein the temperature correction characteristic is stored in the arithmetic device, and the temperature correction characteristic is caused by the load. The thermal diffusion loss information corresponding to area-related information on the size of a stress generation area which is an area in which stress is generated and the period or frequency of the load is included, and the arithmetic device is configured to The heat loss calculation unit extracts the stress generation area, which is an area having a large temperature fluctuation compared to the surroundings in the temperature distribution, and the extracted area related information of the stress generation area and the period of the load Or, based on the frequency and the temperature correction characteristic, the generation area thermal diffusion loss information which is the thermal diffusion loss information corresponding to the stress generation area is determined, and the correction temperature variation calculation unit calculates the stress generation area Before The variation with temperature, on the basis of the said generation region thermal diffusion loss information corresponding to the stress generating areas to correct, obtaining the correction temperature variation in the stress generated region, a stress measuring apparatus.

次に、本発明の第3の発明は、上記第2の発明に係る応力計測装置であって、mを1以上の整数として、前記温度補正特性は、前記応力発生領域から広がるm個の領域である第1〜第mの応力周辺領域のそれぞれの広さに関する第1〜第mの周辺関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、前記演算部は、前記熱損失演算部にて、第1〜第mの前記応力周辺領域を抽出し、抽出した第1〜第mの前記応力周辺領域に対応する第1〜第mの前記周辺関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて、第1〜第mの前記応力周辺領域に対応する前記熱拡散損失情報である第1〜第mの前記熱拡散損失情報を求め、前記補正温度変動量演算部にて、第1〜第mの前記応力周辺領域におけるそれぞれの前記温度変動量を、第1〜第mの前記応力周辺領域のそれぞれに対応する第1〜第mの前記熱拡散損失情報に基づいて補正し、第1〜第mの前記応力周辺領域のそれぞれの前記補正温度変動量を求める、応力計測装置である。   Next, a third invention of the present invention is the stress measurement device according to the second invention, wherein m is an integer of 1 or more, and the temperature correction characteristic includes m regions extending from the stress generation region. The thermal diffusion loss information corresponding to first to mth peripheral related information regarding the width of each of the first to mth stress peripheral areas, and the period or frequency of the load; The arithmetic unit extracts the first to m-th stress peripheral regions in the heat loss arithmetic unit, and extracts the first to m-th corresponding to the extracted first to m-th stress peripheral regions. The first to m-th heat diffusion loss information corresponding to the 1st to m-th stress peripheral regions on the basis of peripheral related information, the period or frequency of the load, and the temperature correction characteristic. The thermal diffusion loss information is obtained, and the first to m-th The respective temperature fluctuation amounts in the force peripheral region are corrected on the basis of the first to m th thermal diffusion loss information corresponding to the first to m th stress peripheral regions, respectively. It is a stress measurement device which calculates | requires each said correction | amendment temperature fluctuation amount of the said stress periphery area | region.

次に、本発明の第4の発明は、上記第3の発明に係る応力計測装置であって、第1〜第mの前記応力周辺領域は、前記応力発生領域に対して同心状に広がるm個のドーナツ状の領域である、応力計測装置である。    Next, according to a fourth aspect of the present invention, in the stress measurement device according to the third aspect, the first to m-th stress peripheral regions extend concentrically with respect to the stress generation region. It is a stress measurement device which is an area of individual donuts.

次に、本発明の第5の発明は、上記第1の発明〜第4の発明のいずれか1つに係る応力計測装置であって、前記演算装置は、前記温度変動量演算部にて、前記負荷付与装置から与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、応力計測装置である。   Next, according to a fifth aspect of the present invention, there is provided the stress measurement device according to any one of the first to fourth aspects, wherein the arithmetic device comprises the temperature fluctuation amount calculation unit. It is a stress measurement device which calculates the said temperature fluctuation amount of each position on the surface of the said measurement object using lock-in processing based on the load given from the load giving device and the temperature distribution.

次に、本発明の第6の発明は、負荷付与装置を用いて、応力の計測対象物に、周期的に変動する負荷を与えながら、赤外線サーモグラフィを用いて、前記計測対象物の表面の温度分布を測定し、演算装置を用いて、与えた前記負荷と、測定した前記温度分布と、に基づいて前記計測対象物の表面の応力分布を求める応力計測方法であって、前記演算装置を用いて、前記温度分布に基づいて、前記計測対象物の表面の各位置の温度変動量を求める温度変動量演算ステップと、前記演算装置を用いて、前記負荷の周期または周波数に応じて、前記計測対象物の表面における熱の拡散による熱の損失に関する熱拡散損失情報を求める熱損失演算ステップと、前記演算装置を用いて、前記温度変動量を前記熱拡散損失情報に基づいて補正した、補正温度変動量を求める補正温度変動量演算ステップと、前記演算装置を用いて、前記補正温度変動量に基づいて、前記負荷に応じた前記計測対象物の表面の応力分布を求める応力演算ステップと、を有する、応力計測方法である。   Next, according to a sixth aspect of the present invention, the temperature of the surface of the measurement object is measured using infrared thermography while applying a load that periodically changes to the measurement object of stress using a load applying device. A stress measurement method of measuring a distribution and determining a stress distribution of a surface of the measurement object based on the applied load and the measured temperature distribution using an arithmetic device, the arithmetic device being used The temperature fluctuation amount calculating step for obtaining the temperature fluctuation amount of each position of the surface of the measurement object based on the temperature distribution, and the measurement according to the cycle or frequency of the load using the arithmetic device A heat loss calculating step for obtaining heat diffusion loss information on heat loss due to heat diffusion on the surface of the object, and the temperature change amount corrected based on the heat diffusion loss information using the arithmetic device, a corrected temperature A corrected temperature fluctuation amount calculating step for obtaining a fluctuation amount; and a stress calculating step for calculating a stress distribution of the surface of the measurement object according to the load based on the corrected temperature fluctuation amount using the arithmetic device. It is a stress measurement method.

次に、本発明の第7の発明は、上記第6の発明に係る応力計測方法であって、前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、前記熱損失演算ステップにおいて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、前記補正温度変動量演算ステップにおいて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、応力計測方法である。   A seventh invention of the present invention is the stress measurement method according to the sixth invention, wherein the temperature correction characteristic is stored in the arithmetic device, and the temperature correction characteristic is caused by the load. The heat diffusion loss information corresponding to the region related information on the width of the stress generation region which is a region in which stress is generated and the cycle or frequency of the load is included, and in the heat loss calculating step Extracting the stress generation area which is an area having a large temperature fluctuation compared to the surroundings in the temperature distribution, the area related information of the stress generation area extracted, the period or frequency of the load, and the temperature The thermal diffusion loss information as the thermal diffusion loss information corresponding to the stress generation area is determined based on the correction characteristic, and the correction temperature fluctuation amount calculating step calculates the thermal diffusion loss information in the stress generation area. The temperature variation, the stress is corrected on the basis of the generation region thermal diffusion loss information corresponding to the generating region, obtaining the correction temperature variation in the stress generated region, a stress measurement method.

次に、本発明の第8の発明は、上記第7の発明に係る応力計測方法であって、mを1以上の整数として、前記温度補正特性は、前記応力発生領域から広がるm個の領域である第1〜第mの応力周辺領域のそれぞれの広さに関する第1〜第mの周辺関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、前記熱損失演算ステップにおいて、第1〜第mの前記応力周辺領域を抽出し、抽出した第1〜第mの前記応力周辺領域に対応する第1〜第mの前記周辺関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて、第1〜第mの前記応力周辺領域に対応する前記熱拡散損失情報である第1〜第mの前記熱拡散損失情報を求め、前記補正温度変動量演算ステップにおいて、第1〜第mの前記応力周辺領域におけるそれぞれの前記温度変動量を、第1〜第mの前記応力周辺領域のそれぞれに対応する第1〜第mの前記熱拡散損失情報に基づいて補正し、第1〜第mの前記応力周辺領域のそれぞれの前記補正温度変動量を求める、応力計測方法である。   An eighth invention of the present invention is the stress measurement method according to the seventh invention, wherein m is an integer of 1 or more, and the temperature correction characteristics are m areas extending from the stress generation area. The thermal diffusion loss information corresponding to first to mth peripheral related information regarding the width of each of the first to mth stress peripheral areas, and the period or frequency of the load; In the heat loss calculating step, first to m-th peripheral related information extracted, and extracted first to m-th peripheral related information corresponding to the extracted first to m-th stress peripheral region; Based on the period or frequency of the load and the temperature correction characteristic, the first to m th thermal diffusion loss information corresponding to the 1 st to m th stress peripheral region is determined The correction temperature variation calculation step, the first to the m-th The respective temperature fluctuation amounts in the stress peripheral region are corrected based on the first to m th thermal diffusion loss information corresponding to the first to m th stress peripheral regions, respectively. In the stress measuring method, the corrected temperature fluctuation amount of each of the stress peripheral region is calculated.

次に、本発明の第9の発明は、上記第8の発明に係る応力計測方法であって、第1〜第mの前記応力周辺領域は、前記応力発生領域に対して同心状に広がるm個のドーナツ状の領域である、応力計測方法である。    Next, according to a ninth aspect of the present invention, in the stress measurement method according to the eighth aspect, the first to m-th stress peripheral regions extend concentrically with respect to the stress generation region. It is a stress measurement method which is an area of individual donuts.

次に、本発明の第10の発明は、上記第6の発明〜第9の発明のいずれか1つに係る応力計測方法であって、前記温度変動量演算ステップにおいて、前記計測対象物に与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、応力計測方法である。   Next, a tenth invention of the present invention is the stress measurement method according to any one of the sixth invention to the ninth invention, wherein the temperature fluctuation amount calculation step is performed on the measurement object. It is a stress measurement method which calculates | requires the said temperature fluctuation amount of each position of the surface of the said measurement object using a lock-in process based on the said load and the said temperature distribution.

第1の発明、及び第6の発明によれば、負荷の周期または周波数に応じて、計測対象物の表面における熱拡散損失情報を求め、当該熱拡散損失情報に基づいて温度変動量を補正して補正温度変動量を得る。そして補正温度変動量に基づいて応力分布を求める。これにより、計測対象物に周期的な負荷を与えながら赤外線サーモグラフィを用いて温度分布を計測し、温度分布に基づいた計測対象物の各位置の温度変動量から応力分布を求める応力計測装置において、熱損失を考慮して、より正確な温度変動量を求め、より正確な応力分布を求めることができる。   According to the first invention and the sixth invention, the thermal diffusion loss information on the surface of the measurement object is determined according to the cycle or frequency of the load, and the temperature fluctuation amount is corrected based on the thermal diffusion loss information. To obtain the corrected temperature fluctuation amount. Then, the stress distribution is determined based on the corrected temperature fluctuation amount. In this way, in a stress measurement device that measures the temperature distribution using infrared thermography while applying cyclical load to the measurement object, and obtains the stress distribution from the temperature fluctuation amount of each position of the measurement object based on the temperature distribution, In consideration of heat loss, more accurate temperature variation can be determined, and more accurate stress distribution can be determined.

第2の発明、及び第7の発明によれば、抽出した応力発生領域の領域関連情報と、負荷の周期または周波数と、に対応した熱拡散損失情報を有する温度補正特性を予め記憶しておくことで、応力発生領域における適切な熱拡散損失情報を得ることができる。従って、応力発生領域における補正温度変動量を適切に求めることが可能であり、応力発生領域における、より正確な応力分布を、比較的容易に求めることができる。   According to the second invention and the seventh invention, the temperature correction characteristic having the thermal diffusion loss information corresponding to the extracted region related information of the stress generation region and the period or frequency of the load is stored in advance. Thus, appropriate heat diffusion loss information in the stress generation region can be obtained. Therefore, it is possible to appropriately determine the correction temperature fluctuation amount in the stress generation region, and to obtain a more accurate stress distribution in the stress generation region relatively easily.

第3の発明、及び第8の発明によれば、応力発生領域の周囲となる応力周辺領域について、応力発生領域と同様に、適切な熱拡散損失情報を得ることができる。従って、応力周辺領域における補正温度変動量を適切に求めることが可能であり、応力周辺領域における、より正確な応力分布を、比較的容易に求めることができる。   According to the third invention and the eighth invention, it is possible to obtain appropriate heat diffusion loss information about the stress peripheral area around the stress generation area as in the stress generation area. Therefore, it is possible to appropriately determine the correction temperature fluctuation amount in the stress peripheral region, and to obtain a more accurate stress distribution in the stress peripheral region relatively easily.

第4の発明、及び第9の発明によれば、応力周辺領域の位置及び形状を適切に定め、各応力周辺領域の補正温度変動量を適切に求めることができる。   According to the fourth invention and the ninth invention, the position and the shape of the stress peripheral area can be appropriately determined, and the correction temperature fluctuation amount of each stress peripheral area can be appropriately determined.

第5の発明、及び第10の発明によれば、ロックイン処理を用いることで、温度分解能を向上させ、より分解能の高い温度変動量を得ることが可能であり、より分解能の高い応力分布を得ることができる。従って、より正確な応力分布を求めることができる。   According to the fifth invention and the tenth invention, by using the lock-in process, it is possible to improve the temperature resolution and obtain a temperature fluctuation amount with higher resolution, and a stress distribution with higher resolution. You can get it. Therefore, more accurate stress distribution can be determined.

応力計測装置の全体構成を説明する図である。It is a figure explaining the whole structure of a stress measurement device. 図1に応力計測装置における演算装置の内部構成の例を説明する図である。FIG. 1 is a view for explaining an example of the internal configuration of a computing device in a stress measurement device. 演算装置の処理手順の例を説明するフローチャートである。It is a flowchart explaining the example of the processing procedure of an arithmetic unit. 負荷付与装置から計測対象物に与える、周期的に変動する負荷の例を説明する図である。It is a figure explaining the example of the load which changes to a measurement object from a load giving device, and which changes periodically. 温度測定装置にて取得した温度分布情報の例を説明する図である。It is a figure explaining an example of temperature distribution information acquired with a temperature measurement device. 応力発生領域、及びその周囲の各画素に対応する温度変動量の例を説明する図である。It is a figure explaining the example of the temperature fluctuation amount corresponding to each stress generation area | region and its surrounding pixel. 温度補正特性の例を説明する図である。It is a figure explaining the example of a temperature correction characteristic. 計測対象物における応力発生領域を抽出する例を説明する図である。It is a figure explaining the example which extracts the stress generating area in a measurement subject. 図8におけるIX矢視方向から見た図であり、応力発生領域の例を説明する図である。It is the figure seen from the IX arrow direction in FIG. 8, and is a figure explaining the example of a stress generation area | region. 計測対象物に与えた負荷と、計測対象物の応力発生領域内の所定位置における温度変動量と補正温度変動量の例を説明する図である。It is a figure explaining the load given to the measurement object, the example of the temperature fluctuation amount in the predetermined position in the stress generating area of the measurement object, and the correction temperature fluctuation amount. 応力発生領域、及びその周囲の各画素に対応する補正温度変動量の例を説明する図である。FIG. 8 is a diagram for explaining an example of a corrected temperature fluctuation amount corresponding to a stress generation region and each pixel around the stress generation region. 応力発生領域、及びその周囲の各画素に対応する応力の例を説明する図である。It is a figure explaining the example of the stress generation area | region and the stress corresponding to each pixel of the periphery. 計測対象物の応力分布の例を説明する図である。It is a figure explaining the example of stress distribution of a measurement subject.

以下、本発明を実施するための形態について図面を用いて説明する。本実施の形態の説明では、計測対象物として等速ジョイント80を用いた例を説明するが、計測対象物は等速ジョイントに限定されず、種々のものを計測対象物とすることができる。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. In the description of the present embodiment, an example in which the constant velocity joint 80 is used as the measurement object will be described, but the measurement object is not limited to the constant velocity joint, and various objects can be used as the measurement object.

●[応力計測装置1の全体構成(図1)と、演算装置の内部構成(図2)]
図1は、本発明の応力計測装置1の全体構成と、計測対象物である等速ジョイント80の外観の例を示している。応力計測装置1は、負荷付与装置10と、温度測定装置20と、演算装置40等を有している。
● [Overall configuration of stress measuring device 1 (FIG. 1) and internal configuration of arithmetic device (FIG. 2)]
FIG. 1 shows an example of the overall configuration of a stress measuring device 1 according to the present invention and an appearance of a constant velocity joint 80 which is an object to be measured. The stress measurement device 1 includes a load application device 10, a temperature measurement device 20, an arithmetic device 40, and the like.

まず計測対象物の例である等速ジョイント80の構成について説明する。本実施の形態にて説明する等速ジョイント80は、いわゆる固定式等速ジョイントであり、アウターレース81、ボールケージ82、インナーレース83、ボール84、アウターシャフト85、インナーシャフト86を有している。   First, the configuration of a constant velocity joint 80 which is an example of a measurement object will be described. The constant velocity joint 80 described in the present embodiment is a so-called fixed constant velocity joint, and has an outer race 81, a ball cage 82, an inner race 83, a ball 84, an outer shaft 85, and an inner shaft 86. .

アウターレース81は、内部が空洞とされた半球状の形状を有しており、内壁にはアウターシャフト85の長手方向に沿った6本の案内溝が形成されている。またアウターレース81には、アウターシャフト85を挿通して固定するためのアウターシャフト孔が形成されている。本実施の形態では、このアウターレース81の表面に発生する応力の分布である応力分布を求める。また温度計測にあたり、アウターレース81の表面からの放射率をほぼ1.0にするために、つや消し黒色に塗装(黒体塗装)されていることが好ましい。   The outer race 81 has a hemispherical shape in which the inside is hollow, and six guide grooves along the longitudinal direction of the outer shaft 85 are formed on the inner wall. Further, an outer shaft hole for inserting and fixing the outer shaft 85 is formed in the outer race 81. In the present embodiment, a stress distribution which is a distribution of stress generated on the surface of the outer race 81 is determined. Moreover, in order to make the emissivity from the surface of the outer race 81 about 1.0 in temperature measurement, it is preferable that the matte black coating (black body coating) is carried out.

インナーレース83は、球体状の形状を有しており、外周面には、インナーシャフト86の長手方向に沿った6本の案内溝が形成されて、アウターレース81に収容されている。またインナーレース83には、インナーシャフト86を挿通して固定するためのインナーシャフト孔が形成されている。   The inner race 83 has a spherical shape, and six guide grooves along the longitudinal direction of the inner shaft 86 are formed on the outer peripheral surface, and the inner race 83 is accommodated in the outer race 81. Further, in the inner race 83, an inner shaft hole for inserting and fixing the inner shaft 86 is formed.

ボール84は、アウターレース81の内壁の案内溝と、インナーレース83の外周面に形成された案内溝と、の間に配置されているとともに、それぞれの案内溝に接している。またインナーレース83の外周面とアウターレース81の内壁との間には、各ボール84を保持するボールケージ82が設けられている。ボールケージ82は、球面板状の形状を有し、ボール84に対応する位置に孔部が形成されている。   The balls 84 are disposed between the guide grooves of the inner wall of the outer race 81 and the guide grooves formed on the outer peripheral surface of the inner race 83, and are in contact with the respective guide grooves. Further, between the outer peripheral surface of the inner race 83 and the inner wall of the outer race 81, a ball cage 82 for holding each ball 84 is provided. The ball cage 82 has a spherical plate shape, and holes are formed at positions corresponding to the balls 84.

アウターシャフト85は、アウターレース81に形成されたアウターシャフト孔に挿通されて固定されている。またインナーシャフト86は、インナーレース83に形成されたインナーシャフト孔に挿通されて固定されている。インナーシャフト86は、インナーシャフト軸線86J回りの回転動力が入力された場合、インナーレース83に回転動力を伝達する。インナーレース83に回転動力が伝達された場合、インナーレース83の回転動力は、ボール84を介してアウターレース81に伝達され、アウターシャフト85をアウターシャフト軸線85J回りに回転させる。   The outer shaft 85 is inserted into and fixed to an outer shaft hole formed in the outer race 81. The inner shaft 86 is inserted into and fixed to an inner shaft hole formed in the inner race 83. The inner shaft 86 transmits rotational power to the inner race 83 when rotational power around the inner shaft axis 86J is input. When the rotational power is transmitted to the inner race 83, the rotational power of the inner race 83 is transmitted to the outer race 81 via the ball 84 to rotate the outer shaft 85 around the outer shaft axis 85J.

負荷付与装置10は、アウターシャフト軸線85Jに対するインナーシャフト軸線86Jの傾斜角度θが、0度から所定角度の範囲内で変化するように、インナーシャフト86を揺動して周期的に変動する負荷(荷重)を与える装置である。この負荷は、所定周期で、負荷ゼロから最大負荷Fmaxの間を繰り返すように与えられる(図4参照)。負荷ゼロの状態では、アウターシャフト軸線85Jとインナーシャフト軸線86Jが一致し、最大負荷Fmaxでは、アウターシャフト軸線85Jに対するインナーシャフト軸線86Jの傾斜角度θが最大となる。なお、最大負荷Fmaxは、実際の車両に取り付けられた等速ジョイント80から取得した負荷や、シミュレーションにて求めた負荷等から、適切な値が設定されている。また負荷付与装置10は、例えば演算装置40から、負荷の周波数f[Hz]、負荷の最大負荷Fmax[MPa]、動作の開始が指示されると、負荷ゼロから最大負荷Fmaxまで周期的に変動する負荷を、周波数f[Hz]にて(図4参照)、等速ジョイント80に与える。そして負荷付与装置10は、時間の経過とともに変動させて与えている負荷の大きさに関する負荷信号(図4参照)を演算装置40に出力する。   The load applying device 10 swings the inner shaft 86 to change periodically such that the inclination angle θ of the inner shaft axis 86J with respect to the outer shaft axis 85J changes within a range of 0 degrees to a predetermined angle ( Load)). This load is given to repeat between zero load and maximum load Fmax at a predetermined cycle (see FIG. 4). In the zero load state, the outer shaft axis 85J and the inner shaft axis 86J coincide with each other, and at the maximum load Fmax, the inclination angle θ of the inner shaft axis 86J with respect to the outer shaft axis 85J is maximum. The maximum load Fmax is set to an appropriate value from the load acquired from the constant velocity joint 80 attached to the actual vehicle, the load obtained by simulation, and the like. Also, for example, when the load frequency f [Hz], the maximum load Fmax [MPa], and the start of operation are instructed from the arithmetic device 40, the load application device 10 periodically changes from zero load to the maximum load Fmax. Load is applied to the constant velocity joint 80 at a frequency f [Hz] (see FIG. 4). Then, the load application device 10 outputs a load signal (see FIG. 4) regarding the magnitude of the load which is fluctuated with the passage of time to the arithmetic device 40.

温度測定装置20は、例えば赤外線カメラであり、赤外線サーモグラフィを用いて、等速ジョイント80のアウターレース81(計測対象物)の表面の温度分布を測定する。例えば、温度測定装置20は、2次元状に配置された複数の画素を有し、各画素に温度情報を対応させた温度分布画像(図5参照)を次々に作成し、作成した温度分布画像を次々と演算装置40に出力する。また温度測定装置20は、例えば演算装置40から、撮影速度(200[fps]等)と、撮影時間と、動作の開始が指示されると、指示された撮影速度で、指示された撮影時間の間、次々と温度分布画像を作成し、作成した温度分布画像を演算装置40に出力する。なお「温度分布画像」は、「温度分布」に相当している。   The temperature measurement device 20 is, for example, an infrared camera, and measures the temperature distribution on the surface of the outer race 81 (object to be measured) of the constant velocity joint 80 using infrared thermography. For example, the temperature measurement device 20 has a plurality of pixels arranged in a two-dimensional manner, and creates temperature distribution images (see FIG. 5) in which the temperature information is made to correspond to each pixel one after another. Are output to the arithmetic unit 40 one after another. Further, the temperature measuring device 20 may receive, for example, the imaging speed (200 fps), the imaging time, and the instructed imaging speed at the instructed imaging speed when instructed to start the operation from the arithmetic device 40, for example. In the meantime, the temperature distribution image is created one after another, and the created temperature distribution image is output to the computing device 40. The “temperature distribution image” corresponds to the “temperature distribution”.

演算装置40は、例えばパーソナルコンピュータであり、CPU等の演算手段を有している。そして演算装置40には、負荷付与装置10から負荷信号が入力され、温度測定装置20から温度分布画像が入力される。そして演算装置40は、負荷信号と温度分布画像に基づいて、等速ジョイント80(計測対象物)の応力分布を求めて表示する。なお図1では、演算装置40が、いわゆるロックイン処理を行うロックイン機能を有している(内蔵している)例を示しているが、ロックイン処理を行うロックイン装置を、演算装置の外部に設けるようにしてもよい。本実施の形態の説明では、演算装置40内にロックイン処理の機能を有している例にて説明する。なお、ロックイン処理は既存の処理であり、外的に繰り返し荷重変動(負荷変動)を付与し、荷重信号に同期する温度変化のみを演算し、他の不定期な周囲の温度変化を排除するものである。このロックイン処理を用いることで、例えば赤外線カメラで計測した温度分布画像であっても、温度分解能を0.0001[℃]程度まで可能とすることができる。例えば、既存のロックイン機能を有するロックインアンプの場合、計測信号(この場合、温度分布画像)と参照信号(この場合、周期的に繰り返す負荷信号)の2つを入力し、それらをPSD(Phase Sensitive Detector)で乗算し、その後、LPF(Low Pass Filter)を通して、目的の信号(この場合、温度変動量の信号)の強度成分を持つ出力信号を得る。参照信号と同じ周波数を持つ信号(目的の信号)は、PSDでの乗算により、直流成分に変換され、LPFを通して取り出される。また、雑音は周波数成分を保持するため、LPFで除去される。   The arithmetic device 40 is, for example, a personal computer, and includes arithmetic means such as a CPU. Then, a load signal is input from the load application device 10 to the arithmetic device 40, and a temperature distribution image is input from the temperature measurement device 20. Then, the computing device 40 obtains and displays the stress distribution of the constant velocity joint 80 (measurement object) based on the load signal and the temperature distribution image. Although FIG. 1 shows an example in which the computing device 40 has (includes) a lock-in function for performing so-called lock-in processing, the lock-in device for performing lock-in processing It may be provided outside. In the description of the present embodiment, an example in which the function of lock-in processing is included in the arithmetic device 40 will be described. The lock-in process is an existing process, externally applies load fluctuation (load fluctuation), calculates only the temperature change synchronized with the load signal, and eliminates other irregular ambient temperature changes. It is a thing. By using this lock-in process, for example, even in the case of a temperature distribution image measured by an infrared camera, the temperature resolution can be made about 0.0001 [° C.]. For example, in the case of a lock-in amplifier having an existing lock-in function, two of a measurement signal (in this case, a temperature distribution image) and a reference signal (in this case, a periodically repeating load signal) are input Multiplication is performed by a Phase Sensitive Detector, and then an LPF (Low Pass Filter) is used to obtain an output signal having an intensity component of a target signal (in this case, a signal of a temperature fluctuation amount). A signal (the target signal) having the same frequency as the reference signal is converted to a direct current component by multiplication with the PSD and extracted through the LPF. Also, noise is removed by the LPF in order to retain frequency components.

また、演算装置40は、図2に示すように、負荷付与装置駆動指示部41A、温度測定装置駆動指示部41B、情報取得部41C、温度変動量演算部41D、熱損失演算部41E、補正温度変動量演算部41F、応力演算部41G、入力手段42A(キーボードやマウス等)、表示手段42B(液晶モニタ等)、記憶手段42C(Hard Disk Drive等)などを有している。   In addition, as shown in FIG. 2, the arithmetic device 40 is, as shown in FIG. 2, a load application device drive instruction unit 41A, a temperature measurement device drive instruction unit 41B, an information acquisition unit 41C, a temperature fluctuation amount operation unit 41D, a heat loss operation unit 41E, and a correction temperature It has a fluctuation amount calculation unit 41F, a stress calculation unit 41G, an input unit 42A (a keyboard, a mouse, etc.), a display unit 42B (a liquid crystal monitor, etc.), a storage unit 42C (Hard Disk Drive, etc.) and the like.

●[演算装置40の処理手順(図3)と、各処理の概要(図4〜図13)]
演算装置40は、作業者から動作の開始が指示されると、図3のフローチャートに示す処理を開始する。
[[Processing procedure of the computing device 40 (FIG. 3) and outline of each processing (FIG. 4 to 13)]]
Arithmetic device 40 starts the process shown in the flowchart of FIG. 3 when instructed by the operator to start the operation.

ステップS10にて演算装置40は、負荷付与装置10に動作開始情報を出力し、負荷付与装置10から等速ジョイント80に、周期的に変動する負荷を与えさせ、付与している負荷の状態を示す負荷信号を出力させ、ステップS20に進む。例えば負荷付与装置10への動作開始情報には、作業者から入力された負荷の周波数f[Hz]や最大負荷Fmax[Mpa]、及び動作を開始するコマンド等が含まれている。このステップS10を実行している演算装置40の演算手段(CPU)は、図2に示す負荷付与装置駆動指示部41Aとして機能する。   In step S10, the computing device 40 outputs operation start information to the load applying device 10, causes the constant velocity joint 80 to apply a periodically fluctuating load from the load applying device 10, and gives the state of the applied load. The load signal shown is output, and the process proceeds to step S20. For example, the operation start information to the load application device 10 includes the frequency f [Hz] of the load input from the operator, the maximum load Fmax [Mpa], a command for starting the operation, and the like. The computing means (CPU) of the computing device 40 executing this step S10 functions as a load applying device drive instructing unit 41A shown in FIG.

ステップS20にて演算装置40は、温度測定装置20に設定情報を出力し、温度測定装置20の設定を行い、ステップS30に進む。例えば温度測定装置20への設定情報には、作業者から入力された撮影速度(200[fps]等)や撮影時間等が含まれている。   In step S20, the arithmetic device 40 outputs setting information to the temperature measurement device 20, performs setting of the temperature measurement device 20, and proceeds to step S30. For example, the setting information for the temperature measurement device 20 includes a shooting speed (200 [fps] or the like), a shooting time, and the like input from the operator.

ステップS30にて演算装置40は、負荷付与装置10から図4の例に示すような負荷信号をリアルタイムに取得する。なお、図4は、演算装置40に入力される負荷信号をグラフ化した例を示しており、周波数f[Hz]にて、負荷ゼロから最大負荷Fmax[MPa]の負荷が、等速ジョイント80に与えられていることを示している。そして演算装置40は、負荷信号の位相に同期させて、温度測定装置20の動作を指示する。例えば、演算装置40は、取得した負荷信号の位相が所定の位相であることを検出した場合に、温度計測装置20に動作情報を出力する。この動作情報には、設定されている撮影速度で、設定されている時間の間、温度分布画像を作成し、作成した温度分布画像を演算装置40に出力するコマンド等が含まれている。演算装置40は、所定時間分あるいは所定周期分の温度分布画像を取得すると、ステップS40に進む。なお図5は、図4における時刻[t1](負荷=0)の場合、時刻[t2](負荷=F1)の場合、時刻[t3](負荷=Fmax)の場合、の各温度分布画像の例を示しており、濃い色の部分は温度が高いことを示している。ステップS20、S30を実行している演算装置40の演算手段(CPU)は、図2に示す温度測定装置駆動指示部41Bとして機能する。またステップS30を実行している演算装置40の演算手段(CPU)は、負荷信号や温度分布情報を取得する情報取得部41C(図2参照)として機能する。   In step S30, the arithmetic device 40 acquires a load signal as shown in the example of FIG. 4 from the load application device 10 in real time. FIG. 4 shows an example in which the load signal input to the arithmetic unit 40 is graphed, and at a frequency f [Hz], the load from zero load to the maximum load Fmax [MPa] is constant velocity joint 80 It shows that it is given to. Then, the arithmetic device 40 instructs the operation of the temperature measurement device 20 in synchronization with the phase of the load signal. For example, when it is detected that the phase of the acquired load signal is a predetermined phase, the arithmetic device 40 outputs the operation information to the temperature measurement device 20. The operation information includes a command for generating a temperature distribution image at a set imaging speed and for a set time, and outputting the generated temperature distribution image to the arithmetic device 40. When the arithmetic device 40 acquires a temperature distribution image for a predetermined time or a predetermined cycle, the processing proceeds to step S40. FIG. 5 shows temperature distribution images of time [t1] (load = 0), time [t2] (load = F1), time [t3] (load = Fmax) in FIG. An example is shown, where darker parts indicate higher temperatures. The computing means (CPU) of the computing device 40 executing steps S20 and S30 functions as a temperature measurement device drive instruction unit 41B shown in FIG. The computing means (CPU) of the computing device 40 executing step S30 functions as an information acquisition unit 41C (see FIG. 2) for acquiring load signals and temperature distribution information.

ステップS40にて演算装置40は、ロックイン処理を用いて、温度分布画像に基づいて温度変動量を求め、ステップS60に進む。温度分布画像は、2次元状に配置された複数の画素に対応させて温度情報が記憶されている。演算装置40は、ロックイン処理を用いて、各画素に対応させて、環境温度の影響や摩擦熱の影響等のノイズ成分を排除した微小な温度変動を求める。そして演算装置40は、各画素に対応させた(すなわち、各画素に相当する等速ジョイントの表面の位置に対応させた)温度変動量を求める。これにより、赤外線サーモグラフィの分解能を桁違いに向上させ、分解能が高く、より正確な温度変動を求める。図6は、図5における領域A1〜A3及びその周囲の画素と、各画素に対応させて求めた温度変動量の例のイメージを示している。2次元状に配置された各画素G(u、v)のそれぞれに、温度変動量が対応させて記憶されている。図6は、例えば領域A1内の画素G(u5、v5)に対応させて記憶されている温度変動量が、Δt[G(u5、v5)]であることを示している。ステップS40の処理は、温度分布画像に基づいて温度変動量を求める温度変動量演算ステップに相当している。そしてステップS40を実行している演算装置40の演算手段(CPU)は、図2に示す温度変動量演算部41Dとして機能する。   In step S40, the arithmetic unit 40 obtains the temperature fluctuation amount based on the temperature distribution image using lock-in processing, and proceeds to step S60. In the temperature distribution image, temperature information is stored corresponding to a plurality of pixels arranged in a two-dimensional manner. The arithmetic unit 40 uses lock-in processing to obtain minute temperature fluctuations excluding noise components, such as the influence of environmental temperature and the influence of frictional heat, corresponding to each pixel. Then, the arithmetic unit 40 obtains the temperature fluctuation amount corresponding to each pixel (that is, corresponding to the position of the surface of the constant velocity joint corresponding to each pixel). Thereby, the resolution of the infrared thermography is improved by an order of magnitude, the resolution is high, and more accurate temperature fluctuation is determined. FIG. 6 shows images of the regions A1 to A3 in FIG. 5 and surrounding pixels, and an example of the temperature variation obtained corresponding to each pixel. The amount of temperature change is stored in correspondence with each of the pixels G (u, v) arranged in a two-dimensional manner. FIG. 6 shows that, for example, the amount of temperature change stored corresponding to the pixel G (u5, v5) in the area A1 is Δt [G (u5, v5)]. The process of step S40 corresponds to a temperature fluctuation calculation step of calculating the temperature fluctuation based on the temperature distribution image. The computing means (CPU) of the computing device 40 executing step S40 functions as a temperature fluctuation amount computing unit 41D shown in FIG.

ステップS60にて演算装置40は、温度分布画像内において周囲と比較して温度変動量が大きい円状の領域である応力発生領域を抽出する。図6の例では、周囲と比較して温度変動量が大きい領域A1を応力発生領域として抽出する。図6の例は説明のためのイメージであり、領域A1が4個の画素で構成されているので円状になっていないが、実際には多数の画素で領域A1が構成されて円状になる。なお、周囲と比較して温度変動量が大きい領域を抽出する方法として、以下の(1)、(2)、(3)の例を示すが、この例に限定されるものではない。
(1)適宜設定した第1閾値以上の温度変動量を有する画素を含む領域を、応力発生領域として抽出する。
(2)図5の例に示すように、予めボール84を含む領域B1、B2、B3を設定しておき、各領域B1、B2、B3内で最大温度変動量ΔTmaxを求める。そして、例えば0.8*ΔTmax以上の温度変動量を有する画素を含む領域を、応力発生領域として抽出する。
(3)図8の例に示すように、応力発生領域は、ボール84とアウターレース81との接触面Sである。なお、図9に示すように、ボール84はアウターレース81の内壁を接触しており、この接触面をS1とすると、図8に示す接触面Sは、図9に示す接触面S1をアウターレース81の表面に投影した面となる。図8に示す接触面Sは、予め計測しておくことが可能である。従って、予め計測した接触面Sの領域を、応力発生領域として抽出する。
In step S60, the arithmetic unit 40 extracts a stress generation area which is a circular area having a large amount of temperature fluctuation in the temperature distribution image compared to the surroundings. In the example of FIG. 6, a region A1 in which the amount of temperature fluctuation is larger than that of the surroundings is extracted as a stress generation region. The example of FIG. 6 is an image for explanation, and although the area A1 is not circular because it is composed of four pixels, the area A1 is actually composed of a large number of pixels and is circular. Become. In addition, although the example of the following (1), (2), and (3) is shown as a method of extracting the area | region where temperature variation is large compared with circumference | surroundings, it is not limited to this example.
(1) A region including a pixel having a temperature variation amount equal to or more than the appropriately set first threshold value is extracted as a stress generation region.
(2) As shown in the example of FIG. 5, areas B1, B2 and B3 including the ball 84 are set in advance, and the maximum temperature fluctuation amount ΔTmax is determined in each of the areas B1, B2 and B3. Then, for example, a region including a pixel having a temperature fluctuation amount of 0.8 * ΔTmax or more is extracted as a stress generation region.
(3) As shown in the example of FIG. 8, the stress generation region is the contact surface S between the ball 84 and the outer race 81. As shown in FIG. 9, the ball 84 is in contact with the inner wall of the outer race 81. Assuming that this contact surface is S1, the contact surface S shown in FIG. 8 has the contact surface S1 shown in FIG. It is a plane projected on the surface of 81. The contact surface S shown in FIG. 8 can be measured in advance. Therefore, the area of the contact surface S measured in advance is extracted as a stress generation area.

そして演算装置40は、応力発生領域の広さに関する領域関連情報と、負荷の周波数(または周期)と、予め記憶手段に記憶されている温度補正特性と、に基づいて、応力発生領域に対応する熱拡散損失情報である発生領域熱拡散損失情報を求める。なお、上記の手順にて抽出した応力発生領域は円状であり、領域関連情報は、例えば、円状の応力発生領域の直径である。   The arithmetic unit 40 corresponds to the stress generation area based on the area related information on the width of the stress generation area, the frequency (or period) of the load, and the temperature correction characteristic stored in advance in the storage means. The generation area heat diffusion loss information which is the heat diffusion loss information is determined. In addition, the stress generation area | region extracted in said procedure is circular, and area | region relevant information is a diameter of a circular stress generation area | region, for example.

温度補正特性は、図7の例に示す特性であり、横軸は負荷の周波数、縦軸は相対温度(熱拡散損失情報であり、この場合、0〜1.0の値)とされ、計測対象物の材質に特有の特性である。「相対温度」は、理想的な断熱状態に対して、実際に発生した熱の拡散による熱の損失にて、どれくらいの割合まで温度が低下しているか、を比率で示している。例えば、領域関連情報が「直径D1」、負荷周波数が「f[Hz]」の場合、図7の例では、相対温度である(発生領域)熱拡散損失情報は、K1(D1、f)である、と求めることができる。この場合、演算装置40は、応力発生領域に対応する発生領域熱拡散損失情報は、K1(D1、f)である、と求めてステップS65に進む。   The temperature correction characteristic is the characteristic shown in the example of FIG. 7, the horizontal axis is the frequency of the load, and the vertical axis is the relative temperature (heat diffusion loss information, in this case, a value of 0 to 1.0). It is a characteristic unique to the material of the object. “Relative temperature” indicates a ratio of how much the temperature is reduced by the loss of heat due to the diffusion of heat actually generated to the ideal adiabatic state. For example, when the area related information is “diameter D1” and the load frequency is “f [Hz]”, in the example of FIG. 7, the thermal diffusion loss information which is the relative temperature (generation area) is K1 (D1, f) It can be asked that there is. In this case, the arithmetic unit 40 determines that the heat generation area heat diffusion loss information corresponding to the stress generation area is K1 (D1, f), and proceeds to step S65.

ステップS65にて演算装置40は、応力発生領域の温度変動量を、発生領域熱拡散損失情報に基づいて補正して、応力発生領域の補正温度変動量を求め、ステップS70に進む。例えば、図6の例にて、応力発生領域(領域A1)内の画素G(u5、v5)の温度変動量がΔt[G(u5、v5)]であり、発生領域熱拡散情報がK1(D1、f)である場合、補正温度変動量をΔT[G(u5、v5)]とすると、以下の(式2)が成立していることになる。この(式2)を整理した(式3)にて、補正温度変動量ΔT[G(u5、v5)]を求めることができる。演算装置40は、応力発生領域内の各画素に対して、それぞれの温度変動量を補正した補正温度変動量を求める。なお、負荷と、温度変動量と、補正温度変動量の関係は、図10に示すとおりである。負荷における+(正)方向への変動は引っ張りを示し、引っ張りの際には温度が下降する側に変動する。また、負荷における−(負)方向への変動は圧縮を示し、圧縮の際には温度が上昇する側に変動する。
ΔT[G(u5、v5)]*K1(D1、f)=Δt[G(u5、v5)] (式2)
ΔT[G(u5、v5)]=Δt[G(u5、v5)]*(1/K1(D1、f)) (式3)
In step S65, the arithmetic unit 40 corrects the temperature fluctuation amount of the stress generation region based on the heat generation loss information of the generation region, obtains the correction temperature fluctuation amount of the stress generation region, and proceeds to step S70. For example, in the example of FIG. 6, the temperature variation of the pixel G (u5, v5) in the stress generation area (area A1) is .DELTA.t [G (u5, v5)], and the generation area thermal diffusion information is K1 ( In the case of D1, f), assuming that the correction temperature fluctuation amount is ΔT [G (u5, v5)], the following (Expression 2) is satisfied. The corrected temperature fluctuation amount ΔT [G (u5, v5)] can be obtained by (Expression 3) which rearranges this (Expression 2). Arithmetic unit 40 obtains a corrected temperature change amount obtained by correcting the temperature change amount for each pixel in the stress generation region. The relationship between the load, the temperature fluctuation amount, and the correction temperature fluctuation amount is as shown in FIG. Fluctuations in the + (positive) direction in the load indicate a pull, and during pull the temperature fluctuates to the falling side. In addition, fluctuation in the negative (-) direction in the load indicates compression, and during compression, it fluctuates to the side where the temperature rises.
ΔT [G (u5, v5)] * K1 (D1, f) = Δt [G (u5, v5)] (Expression 2)
ΔT [G (u5, v5)] = Δt [G (u5, v5)] * (1 / K1 (D1, f)) (Expression 3)

なお、図7に示す温度補正特性は、例えば、以下の(式4)、(式5)、(式6)から作成するようにしてもよいし、シミュレーション等を用いて作成するようにしてもよい。なお、(式4)において、cは計測対象物の比熱、ρは計測対象物の密度、(∂u/∂t)は温度の時間偏微分、kは計測対象物の熱伝導率、rはエネルギーを与えている半径(この場合、ボールとアウターレースとの接触面の半径)、(∂2u/∂r2)は温度が半径方向に広がる加速度、(∂u/∂r)は温度が半径方向に広がる速度、qは投入エネルギーである。また(式5)及び(式6)において、qは投入エネルギー、rはエネルギーを与えている半径(この場合、ボールとアウターレースとの接触面の半径)、bはエネルギーの大きさ、aはエネルギーを与える円状面の直径、ω=2πf(fは負荷周波数)である。なお本実施の形態では、投入エネルギーqを、Asinωtの形式としたが、投入エネルギーqは、Asinωtの形式に限定されるものではなく、種々の形式の投入エネルギーを与えることが可能である。例えば、投入エネルギーqを、ガウス分布状に与えてもよい。
cρ(∂u/∂t)−k[(∂2u/∂r2)+(1/r)(∂u/∂r)]=q(r、t) (式4)
q(r、t)=[(−(b/a)r+b)sinωt (ただし、0≦ωt≦π)(式5)
q(r、t)=0 (ただし、π<ωt) (式6)
Note that the temperature correction characteristic shown in FIG. 7 may be created from, for example, the following (Equation 4), (Equation 5), (Equation 6), or may be created using simulation etc. Good. In equation (4), c is the specific heat of the object to be measured, ρ is the density of the object to be measured, (∂u / ∂t) is the time partial derivative of temperature, k is the thermal conductivity of the object to be measured, r is The radius of energy (in this case, the radius of the contact surface between the ball and the outer race), (∂ 2 u / ∂r 2 ) is the acceleration at which the temperature spreads radially, and (∂u / ∂r) is the temperature The radially spreading velocity, q, is the input energy. Further, in (Equation 5) and (Equation 6), q is the input energy, r is the radius giving the energy (in this case, the radius of the contact surface between the ball and the outer race), b is the size of the energy, a is The diameter of the circular surface giving energy, ω = 2πf (f is the loading frequency). In the present embodiment, although the input energy q is in the form of Asin ωt, the input energy q is not limited to the form of Asin ωt, and various types of input energy can be given. For example, the input energy q may be given in a Gaussian distribution.
cρ (∂u / ∂t) -k [ (∂ 2 u / ∂r 2) + (1 / r) (∂u / ∂r)] = q (r, t) ( Equation 4)
q (r, t) = [(− (b / a) r + b) sin ωt (where 0 ≦ ωt ≦ π) (Equation 5)
q (r, t) = 0 (where π <ωt) (Equation 6)

ステップS70にて演算装置40は、応力発生領域に対して同心状に広がるm個のドーナツ状の領域である第1〜第mの応力周辺領域を抽出する。なお、mは1以上の整数である。図6は、m=2に設定した例を示しており、領域A2を第1の応力周辺領域、領域A3を第2の応力周辺領域として抽出した例を示している。以下、m=2に設定した例で説明する。図6の例は説明のためのイメージであり、領域A2、領域A3ともに円状になっていないが、実際には多数の画素で領域A1、領域A2が構成されて円状になる。なお、応力発生領域(図6の例では領域A1)に対して同心状に広がる領域A2、A3の径を決定する方法として、以下の(a)、(b)の例を示すが、この例に限定されるものではない。
(a)第1閾値、第2閾値、第3閾値(第1閾値>第2閾値>第3閾値(第1閾値>・・・>第m+1閾値))を適宜設定する。上述したように、第1閾値以上の温度変動量を有する画素を含む領域を応力発生領域(領域A1)として抽出する。第1閾値未満かつ第2閾値以上の温度変動量を有する画素を含む領域を第1の応力周辺領域(領域A2)として抽出する。第2閾値未満かつ第3閾値以上の温度変動量を有する画素を含む領域を第2の応力周辺領域(領域A3)として抽出する。
(b)直径D1を有する円状の応力発生領域(領域A1)に対して、同心状の外径D1+ΔDa(ΔDaは適宜設定)かつ内径D1のドーナツ状の領域を第1の応力周辺領域(領域A2)として抽出する。同様に、外径D1+ΔDa+ΔDb(ΔDbは適宜設定)かつ内径D1+ΔDaのドーナツ状の領域を第2の応力周辺領域(領域A2)として抽出する(m個のドーナツ状の応力周辺領域を抽出する)。
In step S70, the arithmetic unit 40 extracts first to m-th stress peripheral regions which are m doughnut-shaped regions concentrically extending with respect to the stress generation region. Here, m is an integer of 1 or more. FIG. 6 shows an example in which m is set to 2 and shows an example in which the region A2 is extracted as a first stress peripheral region and the region A3 is extracted as a second stress peripheral region. Hereinafter, an example in which m is set to 2 will be described. The example of FIG. 6 is an image for explanation, and although both the area A2 and the area A3 are not circular, an area A1 and an area A2 are actually configured by a large number of pixels and become circular. In addition, although the example of the following (a) and (b) is shown as a method of determining the diameter of area | region A2 and A3 which extend concentrically with respect to a stress generation area | region (area A1 in the example of FIG. 6), this example It is not limited to
(A) The first threshold, the second threshold, and the third threshold (first threshold> second threshold> third threshold (first threshold>...> M + 1th threshold)) are appropriately set. As described above, the region including the pixel having the temperature fluctuation amount equal to or more than the first threshold is extracted as the stress generation region (region A1). An area including a pixel having a temperature fluctuation amount smaller than the first threshold and equal to or larger than the second threshold is extracted as a first stress peripheral area (area A2). A region including a pixel having a temperature fluctuation amount smaller than the second threshold and equal to or larger than the third threshold is extracted as a second stress peripheral region (region A3).
(B) A first stress peripheral region (a region having a donut shape of concentric inner diameter D1 + ΔDa (ΔDa is appropriately set) and an inner diameter D1 with respect to a circular stress generation region (region A1) having a diameter D1 Extract as A2). Similarly, a doughnut-shaped area having an outer diameter D1 + ΔDa + ΔDb (ΔDb is appropriately set) and an inner diameter D1 + ΔDa is extracted as a second stress peripheral area (area A2) (m donut-shaped stress peripheral areas are extracted).

そして演算装置40は、第1〜第2(第m)の応力周辺領域の広さに関する第1〜第2(第m)の周辺関連情報と、負荷の周波数(または周期)と、予め記憶手段に記憶されている温度補正特性と、に基づいて、第1〜第2(第m)の応力周辺領域に対応する熱拡散損失情報である第1〜第2(第m)の熱拡散損失情報を求める。なお、上記の手順にて抽出した第1〜第2(第m)の応力周辺領域はドーナツ状であり、例えば、第1の周辺関連情報は第1の応力周辺領域(領域A2)の外径であり、第2の周辺関連情報は第2の応力周辺領域(領域A3)の外径である。   The computing device 40 stores first to second (m-th) peripheral related information related to the width of the first to second (m-th) stress peripheral area, the frequency (or period) of the load, and storage means in advance. The first to second (mth) thermal diffusion loss information, which is the thermal diffusion loss information corresponding to the first to second (mth) stress peripheral region, based on the temperature correction characteristics stored in Ask for The first and second (m-th) stress peripheral regions extracted according to the above procedure are donut-shaped, and, for example, the first peripheral related information is the outer diameter of the first stress peripheral region (region A2) The second peripheral related information is the outer diameter of the second stress peripheral region (region A3).

応力発生領域も応力周辺領域も、どちらも同じ材質のアウターレース81の一部であるので、温度補正特性は、応力発生領域にて説明した図7の例に示す温度補正特性と同じである。例えば、第1の周辺関連情報が「外径D2」の場合、外径D2を「直径D2」とみなし、負荷周波数が「f[Hz]」の場合、図7の例では、相対温度である第1の熱拡散損失情報は、K2(D2、f)である、と求めることができる。また、第2の周辺関連情報が「外径D3」の場合、外径D3を「直径D3」とみなし、負荷周波数が「f[Hz]」の場合、図7の例では、相対温度である第2の熱拡散損失情報は、K3(D3、f)である、と求めることができる。この場合、演算装置40は、第1〜第2(第m)の応力周辺領域に対応する第1〜第2(第m)の熱拡散損失情報は、それぞれK2(D2、f)、K3(D3、f)である、と求めてステップS75に進む。   Since both the stress generation region and the stress peripheral region are part of the outer race 81 of the same material, the temperature correction characteristics are the same as the temperature correction characteristics shown in the example of FIG. 7 described in the stress generation region. For example, when the first peripheral related information is "outside diameter D2", the outside diameter D2 is considered as "diameter D2", and when the load frequency is "f [Hz]", in the example of FIG. The first heat diffusion loss information can be determined as K2 (D2, f). When the second peripheral related information is "outside diameter D3", the outside diameter D3 is considered as "diameter D3", and when the load frequency is "f [Hz]", the relative temperature is in the example of FIG. The second heat diffusion loss information can be determined as K3 (D3, f). In this case, the first to second (mth) thermal diffusion loss information corresponding to the first to second (mth) stress peripheral regions are K2 (D2, f), K3 ( Since it is D3, f), the process proceeds to step S75.

ステップS60、S70の処理は、負荷の周波数(または周期)に応じて、熱拡散損失情報を求める熱損失演算ステップに相当している。そしてステップS60、S70を実行している演算装置40の演算手段(CPU)は、図2に示す熱損失演算部41Eとして機能する。   The processes of steps S60 and S70 correspond to a heat loss calculation step of obtaining heat diffusion loss information according to the frequency (or period) of the load. The computing means (CPU) of the computing device 40 executing steps S60 and S70 functions as the heat loss computing unit 41E shown in FIG.

ステップS75にて演算装置40は、第1〜第2(第m)の応力周辺領域のそれぞれの温度変動量を、第1〜第2(第m)の熱拡散損失情報に基づいて補正して、第1〜第2(第m)の応力周辺領域のそれぞれの補正温度変動量を求め、ステップS80に進む。例えば、図6の例にて、第1の応力周辺領域(領域A2)内の画素G(u3、v3)の温度変動量がΔt[G(u3、v3)]であり、第1の熱拡散情報がK2(D2、f)である場合、補正温度変動量をΔT[G(u3、v3)]とすると、以下の(式7)が成立していることになる。この(式7)を整理した(式8)にて、補正温度変動量ΔT[G(u3、v3)]を求めることができる。また、例えば、図6の例にて、第2の応力周辺領域(領域A3)内の画素G(u7、v6)の温度変動量がΔt[G(u7、v6)]であり、第2の熱拡散情報がK3(D3、f)である場合、補正温度変動量をΔT[G(u7、v6)]とすると、以下の(式9)が成立していることになる。この(式9)を整理した(式10)にて、補正温度変動量ΔT[G(u7、v6)]を求めることができる。演算装置40は、第1〜第2の応力周辺領域内の各画素に対して、それぞれの温度変動量を補正した補正温度変動量を求める。図11は、図6に示す温度変動量を補正した、補正温度変動量の例を示している。
ΔT[G(u3、v3)]*K2(D2、f)=Δt[G(u3、v3)] (式7)
ΔT[G(u3、v3)]=Δt[G(u3、v3)]*(1/K2(D2、f)) (式8)
ΔT[G(u7、v6)]*K3(D3、f)=Δt[G(u7、v6)] (式9)
ΔT[G(u7、v6)]=Δt[G(u7、v6)]*(1/K3(D3、f)) (式10)
In step S75, the arithmetic unit 40 corrects each temperature fluctuation amount of the first to second (m-th) stress peripheral region based on the first to second (m-th) thermal diffusion loss information. The corrected temperature fluctuation amount of each of the first to second (m-th) stress peripheral regions is determined, and the process proceeds to step S80. For example, in the example of FIG. 6, the temperature variation of the pixel G (u3, v3) in the first stress peripheral region (region A2) is Δt [G (u3, v3)], and the first thermal diffusion is When the information is K2 (D2, f), the following (Equation 7) is satisfied, assuming that the correction temperature fluctuation amount is ΔT [G (u3, v3)]. The corrected temperature fluctuation amount ΔT [G (u3, v3)] can be obtained by (Expression 8) which rearranges (Expression 7). Further, for example, in the example of FIG. 6, the temperature variation of the pixel G (u7, v6) in the second stress peripheral region (region A3) is Δt [G (u7, v6)], and the second In the case where the thermal diffusion information is K3 (D3, f), assuming that the correction temperature fluctuation amount is ΔT [G (u7, v6)], the following (Expression 9) is satisfied. The corrected temperature fluctuation amount ΔT [G (u7, v6)] can be obtained by (Expression 10) in which this (Expression 9) is arranged. Arithmetic unit 40 obtains a corrected temperature change amount obtained by correcting the temperature change amount for each pixel in the first and second stress peripheral regions. FIG. 11 shows an example of the corrected temperature change amount obtained by correcting the temperature change amount shown in FIG.
ΔT [G (u3, v3)] * K2 (D2, f) = Δt [G (u3, v3)] (Equation 7)
ΔT [G (u3, v3)] = Δt [G (u3, v3)] * (1 / K2 (D2, f)) (Equation 8)
ΔT [G (u7, v6)] * K3 (D3, f) = Δt [G (u7, v6)] (Equation 9)
ΔT [G (u7, v6)] = Δt [G (u7, v6)] * (1 / K3 (D3, f)) (Equation 10)

ステップS65、S75の処理は、温度変動量を熱拡散損失情報にて補正して補正温度変動量を求める補正温度変動量演算ステップに相当している。そしてステップS65、S75を実行している演算装置40の演算手段(CPU)は、図2に示す補正温度変動量演算部41Fとして機能する。   The processes of steps S65 and S75 correspond to a correction temperature fluctuation amount calculation step of correcting the temperature fluctuation amount with the thermal diffusion loss information to obtain a correction temperature fluctuation amount. The computing means (CPU) of the computing device 40 executing steps S65 and S75 functions as a correction temperature fluctuation amount computing unit 41F shown in FIG.

ステップS80にて演算装置40は、補正温度変動量に基づいて応力分布を求め、ステップS90に進む。なお、応力発生領域と応力周辺領域については、上述したように各画素に対応するそれぞれの補正温度変動量を求めた。しかし、応力発生領域と応力周辺領域を除いた領域は、もともと温度変動量がほとんど無い領域であるので、特に温度変動量を補正する必要が無く、温度変動量=補正温度変動量とみなしてよい。演算装置40は、画素毎に、その画素に対応する補正温度変動量ΔTを用いて、以下の(式11)を用いてその画素に対応する応力変動量Δσを求める。なお、負荷は負荷ゼロから最大負荷Fmaxの間で与えているので、Δσはゼロからの変動量である。従って、画素に対応するΔσは、その画素に対応する最大応力となる。なお(式11)中において、ΔTは補正温度変動量、kは計測対象物を構成している材料に固有の材料定数(熱弾性係数)、Tは材料の絶対温度、Δσは応力変動量である。図12は、図11に示す補正温度変動量から求めた応力変動量(最大応力)の例を示している。
ΔT=−kTΔσ (式11)
In step S80, the arithmetic unit 40 obtains a stress distribution based on the corrected temperature fluctuation amount, and proceeds to step S90. In addition, about the stress generation | occurence | production area | region and stress periphery area | region, as above-mentioned, each corrected temperature fluctuation amount corresponding to each pixel was calculated | required. However, since the region excluding the stress generation region and the stress peripheral region is originally a region where there is almost no temperature variation, there is no need to particularly correct the temperature variation, and temperature variation may be regarded as the correction temperature variation. . The arithmetic unit 40 obtains the stress fluctuation amount Δσ corresponding to the pixel using the following (Equation 11), using the correction temperature fluctuation amount ΔT corresponding to the pixel for each pixel. Since the load is given between zero load and the maximum load Fmax, Δσ is a variation from zero. Therefore, Δσ corresponding to a pixel is the maximum stress corresponding to that pixel. In equation (11), ΔT is a correction temperature fluctuation amount, k is a material constant (thermoelastic coefficient) inherent to the material constituting the measurement object, T is an absolute temperature of the material, and Δσ is a stress fluctuation amount. is there. FIG. 12 shows an example of the stress fluctuation (maximum stress) obtained from the corrected temperature fluctuation shown in FIG.
ΔT = −kTΔσ (Equation 11)

ステップS80の処理は、補正温度変動量に基づいて応力分布を求める応力演算ステップに相当している。そしてステップS80を実行している演算装置40の演算手段(CPU)は、図2に示す応力演算部41Gとして機能する。   The process of step S80 corresponds to a stress calculation step of obtaining a stress distribution based on the corrected temperature fluctuation amount. The computing means (CPU) of the computing device 40 executing step S80 functions as a stress computing unit 41G shown in FIG.

ステップS90にて演算装置40は、画素毎に求めた応力変動量Δσ(すなわち、画素毎の最大応力)を、応力分布として表示手段に表示してステップS95に進む。なお図13に、演算装置40の表示手段に表示された応力分布の例を示す。作業者は、この応力分布を確認することで、肉厚を厚くして剛性を確保するべき個所と、それほど剛性を必要としない個所と、を適切に判断することができる。従って、剛性をそれほど必要としない個所の肉厚を薄くして軽量化することに寄与することができる。   In step S90, the arithmetic unit 40 displays the stress fluctuation amount Δσ (that is, the maximum stress for each pixel) obtained for each pixel as a stress distribution on the display means, and proceeds to step S95. FIG. 13 shows an example of the stress distribution displayed on the display means of the arithmetic device 40. By confirming this stress distribution, the operator can appropriately determine where the thickness should be increased to ensure the rigidity and where the rigidity is not so required. Therefore, it is possible to contribute to reducing the thickness and reducing the weight of a portion which does not require much rigidity.

ステップS95にて演算装置40は、終了指示の有無を判定し、終了指示が入力された場合(Yes)は処理を終了し、終了指示が入力されていない場合(No)はステップS90に戻り、応力分布の表示を維持する。   In step S95, the arithmetic device 40 determines the presence / absence of the termination instruction, and when the termination instruction is input (Yes), the processing is terminated, and when the termination instruction is not input (No), the processing returns to step S90. Maintain an indication of stress distribution.

●[本願の効果]
以上、本実施の形態にて説明した応力計測装置及び応力計測方法は、熱損失を考慮して、より正確な(補正)温度変動量から求めた、より正確な応力分布である。従って、作業者は、剛性が必要な個所と、それほど剛性を必要としない個所と、を適切かつより正確に判断することができる。また熱損失を考慮していない従来では、実際の温度変動量よりも低い温度変動量(熱損失にて温度が低くなる)を計測していたので、実際の応力変動量よりも低い応力変動量を求めていた。これに対して本実施の形態では、熱損失を考慮して、低く計測された温度変動量を実際の温度変動量となるように高い側に補正して、応力を高い側に補正するので、応力に対する剛性不足となることを防止することができる。
● [Effect of this application]
As described above, the stress measurement device and the stress measurement method described in the present embodiment are more accurate stress distributions obtained from more accurate (corrected) temperature fluctuation amounts in consideration of heat loss. Therefore, the operator can appropriately and more accurately determine where stiffness is required and where stiffness is not required. Also, in the prior art in which the heat loss was not taken into consideration, the temperature fluctuation amount lower than the actual temperature fluctuation amount was measured (the temperature is lowered by the heat loss), so the stress fluctuation amount lower than the actual stress fluctuation amount Was seeking. On the other hand, in the present embodiment, in consideration of the heat loss, the low measured temperature fluctuation amount is corrected to the high side to become the actual temperature fluctuation amount, and the stress is corrected to the high side. Insufficient rigidity against stress can be prevented.

本発明の応力計測装置の構造、構成、形状、外観等、及び本発明の応力計測方法の処理、処理手順等は、本発明の要旨を変更しない範囲で種々の変更、追加、削除が可能である。例えば温度分布画像は、図5の例に示すような画像に限定されるものではなく、表示手段に表示された応力分布は、図13の例に限定されるものではない。また処理手順は、図3に示すフローチャートに限定されるものではない。   The structure, configuration, shape, appearance, etc. of the stress measuring device of the present invention, and the processing of the stress measuring method of the present invention, processing procedure, etc. can be variously changed, added, deleted without departing from the scope of the present invention. is there. For example, the temperature distribution image is not limited to the image as shown in the example of FIG. 5, and the stress distribution displayed on the display means is not limited to the example of FIG. Also, the processing procedure is not limited to the flowchart shown in FIG.

本実施の形態の説明では、個々の画素毎に、温度変動量、補正温度変動量、応力変動量(応力)を求める例を説明したが、隣り合う複数の画素をまとめた複数の画素グループを構成し、画素グループ毎に、(平均)温度変動量、(平均)補正温度変動量、(平均)応力変動量(応力)を求めるようにしてもよい。   In the description of the present embodiment, an example in which the temperature fluctuation amount, the correction temperature fluctuation amount, and the stress fluctuation amount (stress) are obtained for each individual pixel has been described, but a plurality of pixel groups in which a plurality of adjacent pixels are put together The (average) temperature variation, the (average) corrected temperature variation, and the (average) stress variation (stress) may be determined for each pixel group.

本実施の形態では、計測対象物として等速ジョイント80を例にして説明したが、本実施の形態にて説明した応力計測装置及び応力計測方法は、等速ジョイントに限定されず、種々の機器の応力分布の計測に適用することができる。   In the present embodiment, the constant velocity joint 80 has been described as an example of the measurement object, but the stress measurement device and the stress measurement method described in the present embodiment are not limited to the constant velocity joint, and various devices It can be applied to the measurement of stress distribution of

本実施の形態では、第1〜第mの応力周辺領域の例として、応力発生領域に対して同心状に広がるm個のドーナツ状の領域としたが、同心状でなくてもよいし、ドーナツ状でなくてもよい。つまり、第1〜第mの応力周辺領域は、応力発生領域から広がるm個の領域であればよく、位置や形状を限定するものではない。   In the present embodiment, as an example of the first to m-th stress peripheral regions, m donut-shaped regions that extend concentrically with respect to the stress generation region are used. It does not have to be in the form. That is, the first to m-th stress peripheral regions may be m regions extending from the stress generation region, and the position and the shape are not limited.

1 応力計測装置
10 負荷付与装置
20 温度測定装置
40 演算装置
80 等速ジョイント(計測対象物)
81 アウターレース
82 ボールケージ
83 インナーレース
84 ボール
85 アウターシャフト
85J アウターシャフト軸線
86 インナーシャフト
86J インナーシャフト軸線
A1 領域(応力発生領域)
A2 領域(第1応力周辺領域)
A3 領域(第2応力周辺領域)
G(u、v) 画素
1 Stress Measurement Device 10 Load Application Device 20 Temperature Measurement Device 40 Arithmetic Device 80 Constant Velocity Joint (Measurement Object)
81 outer race 82 ball cage 83 inner race 84 ball 85 outer shaft 85 J outer shaft axis 86 inner shaft 86 J inner shaft axis A 1 area (stress generation area)
A2 area (area around the first stress)
A3 area (area around the second stress)
G (u, v) pixels

Claims (10)

応力の計測対象物に対して周期的に変動する負荷を与える負荷付与装置と、
前記負荷付与装置にて前記計測対象物に前記負荷を与えながら赤外線サーモグラフィを用いて前記計測対象物の表面の温度分布を測定する温度測定装置と、
前記負荷付与装置から与えた前記負荷と、前記温度測定装置による前記温度分布と、に基づいて前記計測対象物の表面の応力分布を求める演算装置と、
を有する応力計測装置において、
前記演算装置は、
前記温度測定装置による前記温度分布に基づいて、前記計測対象物の表面の各位置の温度変動量を求める温度変動量演算部と、
前記負荷付与装置から与えた前記負荷の周期または周波数に応じて、前記計測対象物の表面における熱の拡散による熱の損失に関する熱拡散損失情報を求める熱損失演算部と、
前記温度変動量を前記熱拡散損失情報に基づいて補正した、補正温度変動量を求める補正温度変動量演算部と、
前記補正温度変動量に基づいて、前記負荷に応じた前記計測対象物の表面の応力分布を求める応力演算部と、
を有する、
応力計測装置。
A load applying device that applies a periodically varying load to a measurement object of stress;
A temperature measurement device that measures the temperature distribution on the surface of the measurement object using infrared thermography while applying the load to the measurement object with the load application device;
An arithmetic device for obtaining a stress distribution on the surface of the measurement object based on the load given from the load application device and the temperature distribution according to the temperature measurement device;
In a stress measuring device having
The arithmetic device is
A temperature fluctuation amount calculation unit which calculates a temperature fluctuation amount at each position of the surface of the measurement object based on the temperature distribution by the temperature measurement device;
A heat loss calculation unit for obtaining heat diffusion loss information on heat loss due to heat diffusion on the surface of the measurement object according to the load cycle or frequency given from the load application device;
A correction temperature fluctuation amount calculation unit for calculating a correction temperature fluctuation amount, in which the temperature fluctuation amount is corrected based on the thermal diffusion loss information;
A stress calculation unit for obtaining a stress distribution on the surface of the measurement object according to the load based on the correction temperature fluctuation amount;
Have
Stress measuring device.
請求項1に記載の応力計測装置であって、
前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記演算装置は、
前記熱損失演算部にて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、
前記補正温度変動量演算部にて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、
応力計測装置。
The stress measurement device according to claim 1, wherein
A temperature correction characteristic is stored in the arithmetic unit, and the temperature correction characteristic is area related information related to the size of a stress generation area which is an area in which a stress is generated due to the load, a cycle of the load or The heat diffusion loss information corresponding to the frequency,
The arithmetic device is
The heat loss calculation unit extracts the stress generation area, which is an area having a large temperature fluctuation compared to the surroundings in the temperature distribution, and the extracted area related information of the stress generation area and the load Based on the period or frequency and the temperature correction characteristic, generation area heat diffusion loss information, which is the heat diffusion loss information corresponding to the stress generation area, is determined
The correction temperature fluctuation calculation unit corrects the temperature fluctuation amount in the stress generation area based on the heat generation area heat diffusion loss information corresponding to the stress generation area, and the correction temperature fluctuation in the stress generation area. Determine the quantity,
Stress measuring device.
請求項2に記載の応力計測装置であって、
mを1以上の整数として、
前記温度補正特性は、前記応力発生領域から広がるm個の領域である第1〜第mの応力周辺領域のそれぞれの広さに関する第1〜第mの周辺関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記演算部は、
前記熱損失演算部にて、第1〜第mの前記応力周辺領域を抽出し、抽出した第1〜第mの前記応力周辺領域に対応する第1〜第mの前記周辺関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて、第1〜第mの前記応力周辺領域に対応する前記熱拡散損失情報である第1〜第mの前記熱拡散損失情報を求め、
前記補正温度変動量演算部にて、第1〜第mの前記応力周辺領域におけるそれぞれの前記温度変動量を、第1〜第mの前記応力周辺領域のそれぞれに対応する第1〜第mの前記熱拡散損失情報に基づいて補正し、第1〜第mの前記応力周辺領域のそれぞれの前記補正温度変動量を求める、
応力計測装置。
The stress measurement device according to claim 2, wherein
Let m be an integer greater than or equal to 1,
The temperature correction characteristic includes first to mth peripheral related information regarding the width of each of the first to mth stress peripheral regions which are m regions extending from the stress generation region, and the period or frequency of the load And the heat diffusion loss information corresponding to
The arithmetic unit is
The heat loss calculation unit extracts the first to m-th stress peripheral regions, and extracts the first to m-th peripheral related information corresponding to the extracted first to m-th stress peripheral regions, and Based on the period or frequency of the load and the temperature correction characteristic, the first to m th thermal diffusion loss information corresponding to the 1 st to m th stress peripheral region is determined ,
In the correction temperature fluctuation amount calculation unit, each of the temperature fluctuation amounts in the first to m-th stress peripheral regions corresponds to each of the first to m-th stress peripheral regions. The correction is made on the basis of the thermal diffusion loss information, and the correction temperature fluctuation amount of each of the first to m-th stress peripheral regions is determined.
Stress measuring device.
請求項3に記載の応力計測装置であって、
第1〜第mの前記応力周辺領域は、前記応力発生領域に対して同心状に広がるm個のドーナツ状の領域である、
応力計測装置。
The stress measurement device according to claim 3, wherein
The first to m-th stress peripheral regions are m donut-shaped regions concentrically extending with respect to the stress generation region,
Stress measuring device.
請求項1〜4のいずれか一項に記載の応力計測装置であって、
前記演算装置は、
前記温度変動量演算部にて、前記負荷付与装置から与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、
応力計測装置。
It is a stress measurement device according to any one of claims 1 to 4,
The arithmetic device is
The temperature fluctuation amount calculation unit calculates the temperature fluctuation amount of each position on the surface of the measurement object using lock-in processing based on the load applied from the load application device and the temperature distribution. Ask,
Stress measuring device.
負荷付与装置を用いて、応力の計測対象物に、周期的に変動する負荷を与えながら、
赤外線サーモグラフィを用いて、前記計測対象物の表面の温度分布を測定し、
演算装置を用いて、与えた前記負荷と、測定した前記温度分布と、に基づいて前記計測対象物の表面の応力分布を求める応力計測方法であって、
前記演算装置を用いて、前記温度分布に基づいて、前記計測対象物の表面の各位置の温度変動量を求める温度変動量演算ステップと、
前記演算装置を用いて、前記負荷の周期または周波数に応じて、前記計測対象物の表面における熱の拡散による熱の損失に関する熱拡散損失情報を求める熱損失演算ステップと、
前記演算装置を用いて、前記温度変動量を前記熱拡散損失情報に基づいて補正した、補正温度変動量を求める補正温度変動量演算ステップと、
前記演算装置を用いて、前記補正温度変動量に基づいて、前記負荷に応じた前記計測対象物の表面の応力分布を求める応力演算ステップと、
を有する、
応力計測方法。
Using a load applying device, while applying a cyclically varying load to a measurement object of stress,
Measuring the temperature distribution of the surface of the measurement object using infrared thermography;
A stress measurement method for determining a stress distribution on a surface of the measurement object based on the applied load and the measured temperature distribution using an arithmetic device,
A temperature fluctuation amount calculating step of calculating a temperature fluctuation amount at each position of the surface of the measurement object based on the temperature distribution using the arithmetic device;
A heat loss calculating step of obtaining heat diffusion loss information on heat loss due to heat diffusion on the surface of the measurement object according to the load cycle or frequency using the calculation device;
A corrected temperature fluctuation amount calculating step of calculating a corrected temperature fluctuation amount by correcting the temperature fluctuation amount based on the thermal diffusion loss information using the arithmetic device;
A stress calculation step of obtaining a stress distribution on the surface of the measurement object according to the load based on the corrected temperature fluctuation amount using the calculation device;
Have
Stress measurement method.
請求項6に記載の応力計測方法であって、
前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記熱損失演算ステップにおいて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、
前記補正温度変動量演算ステップにおいて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、
応力計測方法。
The stress measurement method according to claim 6, wherein
A temperature correction characteristic is stored in the arithmetic unit, and the temperature correction characteristic is area related information related to the size of a stress generation area which is an area in which a stress is generated due to the load, a cycle of the load or The heat diffusion loss information corresponding to the frequency,
In the heat loss calculating step, the stress generation area, which is an area having a large temperature fluctuation compared to the surroundings in the temperature distribution, is extracted, and the area related information of the stress generation area extracted, and the period of the load Or generating region heat diffusion loss information which is the heat diffusion loss information corresponding to the stress generation region based on a frequency and the temperature correction characteristic,
In the correction temperature change amount calculation step, the temperature change amount in the stress generation area is corrected based on the heat generation area heat diffusion loss information corresponding to the stress generation area, and the correction temperature change amount in the stress generation area Ask for
Stress measurement method.
請求項7に記載の応力計測方法であって、
mを1以上の整数として、
前記温度補正特性は、前記応力発生領域から広がるm個の領域である第1〜第mの応力周辺領域のそれぞれの広さに関する第1〜第mの周辺関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記熱損失演算ステップにおいて、第1〜第mの前記応力周辺領域を抽出し、抽出した第1〜第mの前記応力周辺領域に対応する第1〜第mの前記周辺関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて、第1〜第mの前記応力周辺領域に対応する前記熱拡散損失情報である第1〜第mの前記熱拡散損失情報を求め、
前記補正温度変動量演算ステップにおいて、第1〜第mの前記応力周辺領域におけるそれぞれの前記温度変動量を、第1〜第mの前記応力周辺領域のそれぞれに対応する第1〜第mの前記熱拡散損失情報に基づいて補正し、第1〜第mの前記応力周辺領域のそれぞれの前記補正温度変動量を求める、
応力計測方法。
The stress measurement method according to claim 7, wherein
Let m be an integer greater than or equal to 1,
The temperature correction characteristic includes first to mth peripheral related information regarding the width of each of the first to mth stress peripheral regions which are m regions extending from the stress generation region, and the period or frequency of the load And the heat diffusion loss information corresponding to
In the heat loss calculating step, the first to m-th peripheral related information corresponding to the extracted first to m-th stress peripheral region is extracted and the load is extracted. The first to m th thermal diffusion loss information, which is the thermal diffusion loss information corresponding to the 1 st to m th stress peripheral region, is determined based on the period or the frequency of the and the temperature correction characteristic,
In the correction temperature fluctuation amount calculating step, each of the temperature fluctuation amounts in the first to m-th stress peripheral regions corresponds to the first to m-th corresponding to each of the first to m-th stress peripheral regions. The correction is made based on the thermal diffusion loss information, and the correction temperature fluctuation amount of each of the first to m-th stress peripheral regions is obtained.
Stress measurement method.
請求項8に記載の応力計測方法であって、
第1〜第mの前記応力周辺領域は、前記応力発生領域に対して同心状に広がるm個のドーナツ状の領域である、
応力計測方法。
The stress measurement method according to claim 8, wherein
The first to m-th stress peripheral regions are m donut-shaped regions concentrically extending with respect to the stress generation region,
Stress measurement method.
請求項6〜9のいずれか一項に記載の応力計測方法であって、
前記温度変動量演算ステップにおいて、前記計測対象物に与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、
応力計測方法。

The stress measurement method according to any one of claims 6 to 9, wherein
In the temperature variation calculation step, the temperature variation of each position of the surface of the measurement object is determined using lock-in processing based on the load applied to the measurement object and the temperature distribution. ,
Stress measurement method.

JP2017079121A 2017-04-12 2017-04-12 Stress measurement device and stress measurement method Pending JP2018179730A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017079121A JP2018179730A (en) 2017-04-12 2017-04-12 Stress measurement device and stress measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017079121A JP2018179730A (en) 2017-04-12 2017-04-12 Stress measurement device and stress measurement method

Publications (1)

Publication Number Publication Date
JP2018179730A true JP2018179730A (en) 2018-11-15

Family

ID=64275039

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017079121A Pending JP2018179730A (en) 2017-04-12 2017-04-12 Stress measurement device and stress measurement method

Country Status (1)

Country Link
JP (1) JP2018179730A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006153865A (en) * 2004-10-29 2006-06-15 Kobe Steel Ltd Stress measurement method and intensity evaluation method utilizing infrared imager
US20120062751A1 (en) * 2009-05-14 2012-03-15 Christian Homma Capture of thermal images of an object
JP2012103124A (en) * 2010-11-10 2012-05-31 Jtekt Corp Stress measuring method
JP2015001392A (en) * 2013-06-13 2015-01-05 パナソニック株式会社 Infrared stress measurement method and infrared stress measurement system
JP2015087375A (en) * 2013-09-25 2015-05-07 株式会社ジェイテクト Stress measurement system, stress measurement method and stress measurement processor
JP2016142679A (en) * 2015-02-04 2016-08-08 株式会社ジェイテクト Infrared stress measurement method and infrared stress measurement device
JP6034524B1 (en) * 2016-04-08 2016-11-30 Jfeテクノリサーチ株式会社 Resonance frequency estimation method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006153865A (en) * 2004-10-29 2006-06-15 Kobe Steel Ltd Stress measurement method and intensity evaluation method utilizing infrared imager
US20120062751A1 (en) * 2009-05-14 2012-03-15 Christian Homma Capture of thermal images of an object
JP2012103124A (en) * 2010-11-10 2012-05-31 Jtekt Corp Stress measuring method
JP2015001392A (en) * 2013-06-13 2015-01-05 パナソニック株式会社 Infrared stress measurement method and infrared stress measurement system
JP2015087375A (en) * 2013-09-25 2015-05-07 株式会社ジェイテクト Stress measurement system, stress measurement method and stress measurement processor
JP2016142679A (en) * 2015-02-04 2016-08-08 株式会社ジェイテクト Infrared stress measurement method and infrared stress measurement device
JP6034524B1 (en) * 2016-04-08 2016-11-30 Jfeテクノリサーチ株式会社 Resonance frequency estimation method

Similar Documents

Publication Publication Date Title
Sirohi et al. Measurement of helicopter rotor blade deformation using digital image correlation
Poozesh et al. Large-area photogrammetry based testing of wind turbine blades
Baqersad et al. Photogrammetry and optical methods in structural dynamics–A review
Ali et al. Digital image correlation (DIC) technique in measuring strain using opensource platform Ncorr
Wang et al. Dynamic modeling of robotic fish with a base-actuated flexible tail
Ong et al. Development of an economic wireless human motion analysis device for quantitative assessment of human body joint
Warren et al. Optical non-contacting vibration measurement of rotating turbine blades II
Peng et al. Phase-based noncontact vibration measurement of high-speed magnetically suspended rotor
Lee et al. 3-D measurement of structural vibration using digital close-range photogrammetry
Mashimo et al. Experimental verification of elliptical motion model in traveling wave ultrasonic motors
EP3561468B1 (en) Stress measurement method, stress measurement device, and stress measurement system
Poozesh et al. A multi-camera stereo DIC system for extracting operating mode shapes of large scale structures
JP2011149471A (en) Calculation method, device and program for calculating torque of four-point contact ball bearing
JP5375153B2 (en) Stress measuring method and stress measuring apparatus for constant velocity joint
Khadka et al. Strain monitoring of wind turbines using a semi-autonomous drone
Kumar et al. Experimental vibration analysis of large structures using 3D DIC technique with a novel calibration method
Chumacero-Polanco et al. Basin of attraction and limit cycle oscillation amplitude of an ankle-hip model of balance on a balance board
Ye et al. Full-field motion and deformation measurement of high speed rotation based on temporal phase-locking and 3D-DIC
Chen et al. Development of an upper limb rehabilitation system using inertial movement units and kinect device
US20190310174A1 (en) Fatigue limit stress specification system, fatigue limit stress specification device, and fatigue limit stress specification method
JP2018179730A (en) Stress measurement device and stress measurement method
Royo Sánchez et al. Development of a new calibration procedure and its experimental validation applied to a human motion capture system
Yang et al. Harmonic transfer function to measure translational and rotational velocities with continuous-scan laser Doppler vibrometry
Lawson et al. Measurement of deformation of rotating blades using digital image correlation
Zhang et al. High-accuracy three-dimensional reconstruction of vibration based on stereo vision

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200317

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210106

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210126

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20210720