JP2018179730A - Stress measurement device and stress measurement method - Google Patents
Stress measurement device and stress measurement method Download PDFInfo
- 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
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
Abstract
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
特許文献1に記載の応力測定方法では、繰り返し荷重変動を生じさせている際(周期的な負荷を与えている際)、計測対象物には荷重変動に応じた(周期的な)温度変動が発生するが、荷重変動の開始時点から時間の経過に応じて計測対象物の全体の温度が徐々に上昇する。この温度上昇分を補正して、より正確な温度分布を求め、より正確な応力分布を求めている。なお、温度分布中の各位置の温度変動量から、各位置の応力を求める際、一般的には、下記の(式1)から求めている。なお(式1)中において、ΔTは温度変動量、kは計測対象物を構成している材料に固有の材料定数(熱弾性係数)、Tは材料の絶対温度、Δσは応力変動量である。
ΔT=−kTΔσ (式1)
In the stress measurement method described in
Δ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.
以下、本発明を実施するための形態について図面を用いて説明する。本実施の形態の説明では、計測対象物として等速ジョイント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
まず計測対象物の例である等速ジョイント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
アウターレース81は、内部が空洞とされた半球状の形状を有しており、内壁にはアウターシャフト85の長手方向に沿った6本の案内溝が形成されている。またアウターレース81には、アウターシャフト85を挿通して固定するためのアウターシャフト孔が形成されている。本実施の形態では、このアウターレース81の表面に発生する応力の分布である応力分布を求める。また温度計測にあたり、アウターレース81の表面からの放射率をほぼ1.0にするために、つや消し黒色に塗装(黒体塗装)されていることが好ましい。
The
インナーレース83は、球体状の形状を有しており、外周面には、インナーシャフト86の長手方向に沿った6本の案内溝が形成されて、アウターレース81に収容されている。またインナーレース83には、インナーシャフト86を挿通して固定するためのインナーシャフト孔が形成されている。
The
ボール84は、アウターレース81の内壁の案内溝と、インナーレース83の外周面に形成された案内溝と、の間に配置されているとともに、それぞれの案内溝に接している。またインナーレース83の外周面とアウターレース81の内壁との間には、各ボール84を保持するボールケージ82が設けられている。ボールケージ82は、球面板状の形状を有し、ボール84に対応する位置に孔部が形成されている。
The
アウターシャフト85は、アウターレース81に形成されたアウターシャフト孔に挿通されて固定されている。またインナーシャフト86は、インナーレース83に形成されたインナーシャフト孔に挿通されて固定されている。インナーシャフト86は、インナーシャフト軸線86J回りの回転動力が入力された場合、インナーレース83に回転動力を伝達する。インナーレース83に回転動力が伝達された場合、インナーレース83の回転動力は、ボール84を介してアウターレース81に伝達され、アウターシャフト85をアウターシャフト軸線85J回りに回転させる。
The
負荷付与装置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
温度測定装置20は、例えば赤外線カメラであり、赤外線サーモグラフィを用いて、等速ジョイント80のアウターレース81(計測対象物)の表面の温度分布を測定する。例えば、温度測定装置20は、2次元状に配置された複数の画素を有し、各画素に温度情報を対応させた温度分布画像(図5参照)を次々に作成し、作成した温度分布画像を次々と演算装置40に出力する。また温度測定装置20は、例えば演算装置40から、撮影速度(200[fps]等)と、撮影時間と、動作の開始が指示されると、指示された撮影速度で、指示された撮影時間の間、次々と温度分布画像を作成し、作成した温度分布画像を演算装置40に出力する。なお「温度分布画像」は、「温度分布」に相当している。
The
演算装置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
また、演算装置40は、図2に示すように、負荷付与装置駆動指示部41A、温度測定装置駆動指示部41B、情報取得部41C、温度変動量演算部41D、熱損失演算部41E、補正温度変動量演算部41F、応力演算部41G、入力手段42A(キーボードやマウス等)、表示手段42B(液晶モニタ等)、記憶手段42C(Hard Disk Drive等)などを有している。
In addition, as shown in FIG. 2, the
●[演算装置40の処理手順(図3)と、各処理の概要(図4〜図13)]
演算装置40は、作業者から動作の開始が指示されると、図3のフローチャートに示す処理を開始する。
[[Processing procedure of the computing device 40 (FIG. 3) and outline of each processing (FIG. 4 to 13)]]
ステップS10にて演算装置40は、負荷付与装置10に動作開始情報を出力し、負荷付与装置10から等速ジョイント80に、周期的に変動する負荷を与えさせ、付与している負荷の状態を示す負荷信号を出力させ、ステップS20に進む。例えば負荷付与装置10への動作開始情報には、作業者から入力された負荷の周波数f[Hz]や最大負荷Fmax[Mpa]、及び動作を開始するコマンド等が含まれている。このステップS10を実行している演算装置40の演算手段(CPU)は、図2に示す負荷付与装置駆動指示部41Aとして機能する。
In step S10, the
ステップS20にて演算装置40は、温度測定装置20に設定情報を出力し、温度測定装置20の設定を行い、ステップS30に進む。例えば温度測定装置20への設定情報には、作業者から入力された撮影速度(200[fps]等)や撮影時間等が含まれている。
In step S20, the
ステップ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
ステップ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
ステップ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
(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
(3) As shown in the example of FIG. 8, the stress generation region is the contact surface S between the
そして演算装置40は、応力発生領域の広さに関する領域関連情報と、負荷の周波数(または周期)と、予め記憶手段に記憶されている温度補正特性と、に基づいて、応力発生領域に対応する熱拡散損失情報である発生領域熱拡散損失情報を求める。なお、上記の手順にて抽出した応力発生領域は円状であり、領域関連情報は、例えば、円状の応力発生領域の直径である。
The
温度補正特性は、図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
ステップ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
Δ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
(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
応力発生領域も応力周辺領域も、どちらも同じ材質のアウターレース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
ステップ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
ステップ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
Δ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
ステップS80にて演算装置40は、補正温度変動量に基づいて応力分布を求め、ステップS90に進む。なお、応力発生領域と応力周辺領域については、上述したように各画素に対応するそれぞれの補正温度変動量を求めた。しかし、応力発生領域と応力周辺領域を除いた領域は、もともと温度変動量がほとんど無い領域であるので、特に温度変動量を補正する必要が無く、温度変動量=補正温度変動量とみなしてよい。演算装置40は、画素毎に、その画素に対応する補正温度変動量ΔTを用いて、以下の(式11)を用いてその画素に対応する応力変動量Δσを求める。なお、負荷は負荷ゼロから最大負荷Fmaxの間で与えているので、Δσはゼロからの変動量である。従って、画素に対応するΔσは、その画素に対応する最大応力となる。なお(式11)中において、ΔTは補正温度変動量、kは計測対象物を構成している材料に固有の材料定数(熱弾性係数)、Tは材料の絶対温度、Δσは応力変動量である。図12は、図11に示す補正温度変動量から求めた応力変動量(最大応力)の例を示している。
ΔT=−kTΔσ (式11)
In step S80, the
Δ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
ステップS90にて演算装置40は、画素毎に求めた応力変動量Δσ(すなわち、画素毎の最大応力)を、応力分布として表示手段に表示してステップS95に進む。なお図13に、演算装置40の表示手段に表示された応力分布の例を示す。作業者は、この応力分布を確認することで、肉厚を厚くして剛性を確保するべき個所と、それほど剛性を必要としない個所と、を適切に判断することができる。従って、剛性をそれほど必要としない個所の肉厚を薄くして軽量化することに寄与することができる。
In step S90, the
ステップS95にて演算装置40は、終了指示の有無を判定し、終了指示が入力された場合(Yes)は処理を終了し、終了指示が入力されていない場合(No)はステップS90に戻り、応力分布の表示を維持する。
In step S95, the
●[本願の効果]
以上、本実施の形態にて説明した応力計測装置及び応力計測方法は、熱損失を考慮して、より正確な(補正)温度変動量から求めた、より正確な応力分布である。従って、作業者は、剛性が必要な個所と、それほど剛性を必要としない個所と、を適切かつより正確に判断することができる。また熱損失を考慮していない従来では、実際の温度変動量よりも低い温度変動量(熱損失にて温度が低くなる)を計測していたので、実際の応力変動量よりも低い応力変動量を求めていた。これに対して本実施の形態では、熱損失を考慮して、低く計測された温度変動量を実際の温度変動量となるように高い側に補正して、応力を高い側に補正するので、応力に対する剛性不足となることを防止することができる。
● [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
81
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.
前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記演算装置は、
前記熱損失演算部にて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、
前記補正温度変動量演算部にて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、
応力計測装置。 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.
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.
第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.
前記演算装置は、
前記温度変動量演算部にて、前記負荷付与装置から与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、
応力計測装置。 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.
前記演算装置には、温度補正特性が記憶されており、前記温度補正特性は、前記負荷による応力が発生している領域である応力発生領域の広さに関する領域関連情報と、前記負荷の周期または周波数と、に対応させた前記熱拡散損失情報を有しており、
前記熱損失演算ステップにおいて、前記温度分布内において周囲と比較して温度変動量が大きい領域である前記応力発生領域を抽出し、抽出した前記応力発生領域の前記領域関連情報と、前記負荷の周期または周波数と、前記温度補正特性と、に基づいて前記応力発生領域に対応する前記熱拡散損失情報である発生領域熱拡散損失情報を求め、
前記補正温度変動量演算ステップにおいて、前記応力発生領域における前記温度変動量を、前記応力発生領域に対応する前記発生領域熱拡散損失情報に基づいて補正し、前記応力発生領域における前記補正温度変動量を求める、
応力計測方法。 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.
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.
第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.
前記温度変動量演算ステップにおいて、前記計測対象物に与えた前記負荷と、前記温度分布と、に基づいてロックイン処理を用いて、前記計測対象物の表面の各位置の前記温度変動量を求める、
応力計測方法。
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.
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)
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 |
-
2017
- 2017-04-12 JP JP2017079121A patent/JP2018179730A/en active Pending
Patent Citations (7)
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 |