JP2003232688A - 二次元応力場計測システム及び二次元応力場計測プログラム - Google Patents
二次元応力場計測システム及び二次元応力場計測プログラムInfo
- Publication number
- JP2003232688A JP2003232688A JP2002032171A JP2002032171A JP2003232688A JP 2003232688 A JP2003232688 A JP 2003232688A JP 2002032171 A JP2002032171 A JP 2002032171A JP 2002032171 A JP2002032171 A JP 2002032171A JP 2003232688 A JP2003232688 A JP 2003232688A
- Authority
- JP
- Japan
- Prior art keywords
- stress
- distribution
- traction
- component
- test body
- 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.)
- Granted
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
(57)【要約】
【課題】 材料試験体の一様でない変形や応力の状態を
計測し、破壊直前又は破壊時の材料試験体の特性をより
正しく計測する。 【解決手段】 二次元応力場計測システムは、例えば、
載荷系と、撮影系と、解析系とを備える。まず、撮影系
は、載荷系により荷重を受けた試験体の表面の変位分布
を測定し取得する（Ｓ１０１）。解析系の演算部は、こ
の変位分布から試験体の表面のひずみ分布を計算して、
ひずみ分布ファイルに記憶し、トラクション分布をトラ
クション分布ファイルに記憶する（Ｓ１０２）。演算部
は、ひずみ分布ファイルと、試験体の縁のトラクション
分布ファイルとを、メモリからそれぞれ読み出す（Ｓ１
０３１、Ｓ１０３２）。つぎに、演算部は、読み出され
たデータを基に、応力成分に関する境界値問題を設定す
る（Ｓ１０３３）。演算部は、設定された境界値問題
を、いわゆる有限要素法と同様、弱形式を用いた定式化
と数値解法の適用によって解き、応力分布の計算を行い
（Ｓ１０３４）、この計算結果を応力分布ファイルに書
き込む（Ｓ１０３５）。
計測し、破壊直前又は破壊時の材料試験体の特性をより
正しく計測する。 【解決手段】 二次元応力場計測システムは、例えば、
載荷系と、撮影系と、解析系とを備える。まず、撮影系
は、載荷系により荷重を受けた試験体の表面の変位分布
を測定し取得する（Ｓ１０１）。解析系の演算部は、こ
の変位分布から試験体の表面のひずみ分布を計算して、
ひずみ分布ファイルに記憶し、トラクション分布をトラ
クション分布ファイルに記憶する（Ｓ１０２）。演算部
は、ひずみ分布ファイルと、試験体の縁のトラクション
分布ファイルとを、メモリからそれぞれ読み出す（Ｓ１
０３１、Ｓ１０３２）。つぎに、演算部は、読み出され
たデータを基に、応力成分に関する境界値問題を設定す
る（Ｓ１０３３）。演算部は、設定された境界値問題
を、いわゆる有限要素法と同様、弱形式を用いた定式化
と数値解法の適用によって解き、応力分布の計算を行い
（Ｓ１０３４）、この計算結果を応力分布ファイルに書
き込む（Ｓ１０３５）。
Description
【０００１】
【発明の属する技術分野】本発明は、二次元応力場計測
システム、二次元応力場計測方法及びプログラムに係
り、特に、画像解析と逆解析理論を利用した二次元応力
場計測システム、二次元応力場計測方法及びプログラム
に関する。
システム、二次元応力場計測方法及びプログラムに係
り、特に、画像解析と逆解析理論を利用した二次元応力
場計測システム、二次元応力場計測方法及びプログラム
に関する。
【０００２】
【従来の技術】材料試験装置を用いたいわゆる材料試験
は、例えば、材料試験体に荷重をかけ、そのときの材料
試験体の伸びや縮みを計測するものである。この材料試
験では、材料試験体が一様に変形し、さらに、一様に応
力を受けていると仮定し、材料試験体の平均的な応力を
計測している。
は、例えば、材料試験体に荷重をかけ、そのときの材料
試験体の伸びや縮みを計測するものである。この材料試
験では、材料試験体が一様に変形し、さらに、一様に応
力を受けていると仮定し、材料試験体の平均的な応力を
計測している。
【０００３】
【発明が解決しようとする課題】しかしながら、材料試
験体が破壊に至る直前の過程では、実際には、変形や応
力は一様とはならない場合が想定される。これは、破壊
前には亀裂等の損傷が材料試験体に発生し、そこに変形
や応力が集中するためである。各種材料試験体の破壊の
状態を正しく計測することは、各種材料を限界により近
い状態まで効率的に使うことにつながる。したがって、
破壊直前（又は、破壊時）の、特に損傷を受ける部分の
一様でない変形や応力の状態を計測することが必要とな
っている。
験体が破壊に至る直前の過程では、実際には、変形や応
力は一様とはならない場合が想定される。これは、破壊
前には亀裂等の損傷が材料試験体に発生し、そこに変形
や応力が集中するためである。各種材料試験体の破壊の
状態を正しく計測することは、各種材料を限界により近
い状態まで効率的に使うことにつながる。したがって、
破壊直前（又は、破壊時）の、特に損傷を受ける部分の
一様でない変形や応力の状態を計測することが必要とな
っている。
【０００４】本発明は、以上の点に鑑み、材料試験体の
一様でない変形や応力の状態を計測し、破壊直前又は破
壊時の材料試験体の特性をより正しく計測することを目
的とする。また、本発明は、材料試験体の表面の変位の
分布を高精度・高分解能の画像解析によって計測し、新
しく確立された逆解析理論を用いて応力分布を推定・計
測することを目的とする。
一様でない変形や応力の状態を計測し、破壊直前又は破
壊時の材料試験体の特性をより正しく計測することを目
的とする。また、本発明は、材料試験体の表面の変位の
分布を高精度・高分解能の画像解析によって計測し、新
しく確立された逆解析理論を用いて応力分布を推定・計
測することを目的とする。
【０００５】
【課題を解決するための手段】本発明の第１の解決手段
によると、複数の標点を含む試験体に荷重を加え、試験
体の縁又は境界の各位置において加えられる力を示すト
ラクションを測定する載荷系と、前記載荷系により加え
られた荷重による前記試験体の各標点の変位を計測する
撮影系と、前記撮影系により計測された標点の各位置に
対応する少なくとも２成分のひずみ分布と、前記載荷系
により計測された前記試験体の縁又は境界の各位置に対
応する圧力を示す少なくとも２成分のトラクション分布
とに基づき、前記試験体に関する所定の境界値問題を解
くことにより、標点の各位置に対応する３成分の応力分
布を計算する解析系とを備えた二次元応力場計測システ
ムが提供される。
によると、複数の標点を含む試験体に荷重を加え、試験
体の縁又は境界の各位置において加えられる力を示すト
ラクションを測定する載荷系と、前記載荷系により加え
られた荷重による前記試験体の各標点の変位を計測する
撮影系と、前記撮影系により計測された標点の各位置に
対応する少なくとも２成分のひずみ分布と、前記載荷系
により計測された前記試験体の縁又は境界の各位置に対
応する圧力を示す少なくとも２成分のトラクション分布
とに基づき、前記試験体に関する所定の境界値問題を解
くことにより、標点の各位置に対応する３成分の応力分
布を計算する解析系とを備えた二次元応力場計測システ
ムが提供される。
【０００６】本発明の第２の解決手段によると、複数の
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備え、前記解析系は、
さらに、前記撮影系により計測された、標点の各位置に
対応する少なくとも２成分のひずみ量を、ひずみ分布と
して記憶する第１ファイルと、前記載荷系により計測さ
れた、前記試験体の縁又は境界の各位置に対応する圧力
を示す少なくとも２成分のトラクション量を、トラクシ
ョン分布として記憶する第２ファイルと、前記第１ファ
イル及び前記第２ファイルを読み出して、ひずみ分布及
びトラクション分布に基づき設定された前記試験体に関
する境界値問題を解くことにより、標点の各位置に対応
する少なくとも３成分の応力を計算する処理部と、前記
処理部により計算された応力を、応力分布として記憶す
る第３ファイルとを有する二次元応力場計測システムが
提供される。
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備え、前記解析系は、
さらに、前記撮影系により計測された、標点の各位置に
対応する少なくとも２成分のひずみ量を、ひずみ分布と
して記憶する第１ファイルと、前記載荷系により計測さ
れた、前記試験体の縁又は境界の各位置に対応する圧力
を示す少なくとも２成分のトラクション量を、トラクシ
ョン分布として記憶する第２ファイルと、前記第１ファ
イル及び前記第２ファイルを読み出して、ひずみ分布及
びトラクション分布に基づき設定された前記試験体に関
する境界値問題を解くことにより、標点の各位置に対応
する少なくとも３成分の応力を計算する処理部と、前記
処理部により計算された応力を、応力分布として記憶す
る第３ファイルとを有する二次元応力場計測システムが
提供される。
【０００７】本発明の第３の解決手段によると、複数の
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備えた二次元応力場計
測システムに用いられる二次元応力場計測方法であっ
て、前記撮影系により計測された、標点の各位置に対応
する少なくとも３成分のひずみ量を、ひずみ分布として
第１ファイルに記憶する機能と、前記載荷系により計測
された、前記試験体の縁又は境界の各位置に対応する圧
力を示す少なくとも２成分のトラクション量を、トラク
ション分布として第２ファイルに記憶する機能と、前記
第１ファイル及び前記第２ファイルを読み出して、ひず
み分布及びトラクション分布に基づき設定された前記試
験体に関する境界値問題を解くことにより、標点の各位
置に対応する３成分の応力を計算し、応力分布として第
３ファイルに記憶する機能とを含む二次元応力場計測方
法が提供される。
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備えた二次元応力場計
測システムに用いられる二次元応力場計測方法であっ
て、前記撮影系により計測された、標点の各位置に対応
する少なくとも３成分のひずみ量を、ひずみ分布として
第１ファイルに記憶する機能と、前記載荷系により計測
された、前記試験体の縁又は境界の各位置に対応する圧
力を示す少なくとも２成分のトラクション量を、トラク
ション分布として第２ファイルに記憶する機能と、前記
第１ファイル及び前記第２ファイルを読み出して、ひず
み分布及びトラクション分布に基づき設定された前記試
験体に関する境界値問題を解くことにより、標点の各位
置に対応する３成分の応力を計算し、応力分布として第
３ファイルに記憶する機能とを含む二次元応力場計測方
法が提供される。
【０００８】本発明の第４の解決手段によると、複数の
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備えた二次元応力場計
測システムに用いられる二次元応力場計測プログラムで
あって、標点の各位置に対応する少なくとも３成分のひ
ずみ量を、ひずみ分布として第１ファイルに記憶する手
順と、前記撮影系により計測された、前記試験体の縁又
は境界の各位置に対応する圧力を示す少なくとも２成分
のトラクション量を、トラクション分布として第２ファ
イルに記憶する手順と、前記載荷系により計測された、
前記第１ファイル及び前記第２ファイルを読み出して、
ひずみ分布及びトラクション分布に基づき設定された前
記試験体に関する境界値問題を解くことにより、標点の
各位置に対応する３成分の応力を計算し、応力分布とし
て第３ファイルに記憶する手順とをコンピュータに実行
させるための二次元応力場計測プログラムが提供され
る。
標点を含む試験体に荷重を加え、試験体の縁又は境界の
各位置において加えられる力を示すトラクションを測定
する載荷系と、前記載荷系により加えられた荷重による
前記試験体の各標点の変位を計測する撮影系と、前記載
荷系及び前記撮影系で計測された値に基づいて、標点毎
の応力分布を計測する解析系とを備えた二次元応力場計
測システムに用いられる二次元応力場計測プログラムで
あって、標点の各位置に対応する少なくとも３成分のひ
ずみ量を、ひずみ分布として第１ファイルに記憶する手
順と、前記撮影系により計測された、前記試験体の縁又
は境界の各位置に対応する圧力を示す少なくとも２成分
のトラクション量を、トラクション分布として第２ファ
イルに記憶する手順と、前記載荷系により計測された、
前記第１ファイル及び前記第２ファイルを読み出して、
ひずみ分布及びトラクション分布に基づき設定された前
記試験体に関する境界値問題を解くことにより、標点の
各位置に対応する３成分の応力を計算し、応力分布とし
て第３ファイルに記憶する手順とをコンピュータに実行
させるための二次元応力場計測プログラムが提供され
る。
【０００９】
【発明の実施の形態】以下、図面を用いて本発明の実施
の形態を詳細に説明する。図１は、二次元応力場計測シ
ステム１００の概略構成図である。二次元応力場計測シ
ステム１００は、例えば、載荷系１と、撮影系２と、解
析系３とを備える。載荷系１は、例えば、一般的な載荷
装置１１と圧力計１０とを含む。この載荷装置１１は、
例えば、この載荷装置１１に設置された材料試験体（試
験体）１２に荷重を加える。圧力計１０は、例えば、試
験体１２に加えられた圧力値（例えば、トラクション分
布と呼ばれる試験体１２と載荷装置１１荷重部との境界
に働く単位面積毎の力の分布）を、信号ケーブルを介し
て接続された解析系３に出力する。圧力値は、試験体１
２の上端、下端、測端又は複数端に加えられた値を適宜
用いる。なお、圧力計１０としてトランスデューサーを
用いてもよい。撮影系２は、例えば、ディジタルビデオ
カメラ２１を含む。このディジタルビデオカメラ２１
は、例えば、載荷系１に含まれる試験体１２の表面の変
位分布を計測する。なお、撮影系２では、ディジタルビ
デオカメラ２１に限られず、ビデオカメラ、レーザー変
位計等の適宜の撮影手段又は測定手段を用いることがで
きる。また、解析系３は、例えば、ワークステーション
３１を含む。このワークステーション３１は、例えば、
撮影系２のディジタルビデオカメラ２１より得られた変
位分布からひずみ分布を計算し、さらに、載荷系１の圧
力計１０よりトラクション分布を得て、逆解析によって
応力分布を推定・計測するコンピュータである。
の形態を詳細に説明する。図１は、二次元応力場計測シ
ステム１００の概略構成図である。二次元応力場計測シ
ステム１００は、例えば、載荷系１と、撮影系２と、解
析系３とを備える。載荷系１は、例えば、一般的な載荷
装置１１と圧力計１０とを含む。この載荷装置１１は、
例えば、この載荷装置１１に設置された材料試験体（試
験体）１２に荷重を加える。圧力計１０は、例えば、試
験体１２に加えられた圧力値（例えば、トラクション分
布と呼ばれる試験体１２と載荷装置１１荷重部との境界
に働く単位面積毎の力の分布）を、信号ケーブルを介し
て接続された解析系３に出力する。圧力値は、試験体１
２の上端、下端、測端又は複数端に加えられた値を適宜
用いる。なお、圧力計１０としてトランスデューサーを
用いてもよい。撮影系２は、例えば、ディジタルビデオ
カメラ２１を含む。このディジタルビデオカメラ２１
は、例えば、載荷系１に含まれる試験体１２の表面の変
位分布を計測する。なお、撮影系２では、ディジタルビ
デオカメラ２１に限られず、ビデオカメラ、レーザー変
位計等の適宜の撮影手段又は測定手段を用いることがで
きる。また、解析系３は、例えば、ワークステーション
３１を含む。このワークステーション３１は、例えば、
撮影系２のディジタルビデオカメラ２１より得られた変
位分布からひずみ分布を計算し、さらに、載荷系１の圧
力計１０よりトラクション分布を得て、逆解析によって
応力分布を推定・計測するコンピュータである。
【００１０】図２は、試験体１２の拡大図である。試験
体１２の表面には、例えば、複数の標点１３が施されて
いる。なお、この標点１３は、試験体１２の表面の適宜
の位置に施されてもよい。撮影系２は、ディジタルビデ
オカメラ２１により材料実験中の標点１３の動きを撮影
することで、試験体１２の表面の変位分布を計測するこ
とになる。
体１２の表面には、例えば、複数の標点１３が施されて
いる。なお、この標点１３は、試験体１２の表面の適宜
の位置に施されてもよい。撮影系２は、ディジタルビデ
オカメラ２１により材料実験中の標点１３の動きを撮影
することで、試験体１２の表面の変位分布を計測するこ
とになる。
【００１１】図３は、ワークステーション３１の構成図
である。ワークステーション３１は、例えば、インター
フェース（Ｉ／Ｆ）部３５と、演算部（ＣＰＵ）３２
と、メモリ３８とを含む。Ｉ／Ｆ部３５は、例えば、各
種データ（ここでは、載荷系１からのトラクション分布
に関するデータ、撮影系２からのひずみ分布に関するデ
ータ）を入力する。演算部３２は、入力部３５から入力
された各種データに基づいて、計測・演算を行う。メモ
リ３８は、例えば、ひずみ分布ファイル４０、トラクシ
ョン分布ファイル５０、応力分布ファイル６０を含む。
また、演算部３２は、演算した値等を、メモリ３８に含
まれる上述の所定のファイルに書き込み及び／又は読み
出す。つぎに、解析系３での処理について主に説明す
る。解析系３では、撮影系２で撮影された載荷系１より
荷重を受けながら変形する試験体１２の画像の解析を行
う。
である。ワークステーション３１は、例えば、インター
フェース（Ｉ／Ｆ）部３５と、演算部（ＣＰＵ）３２
と、メモリ３８とを含む。Ｉ／Ｆ部３５は、例えば、各
種データ（ここでは、載荷系１からのトラクション分布
に関するデータ、撮影系２からのひずみ分布に関するデ
ータ）を入力する。演算部３２は、入力部３５から入力
された各種データに基づいて、計測・演算を行う。メモ
リ３８は、例えば、ひずみ分布ファイル４０、トラクシ
ョン分布ファイル５０、応力分布ファイル６０を含む。
また、演算部３２は、演算した値等を、メモリ３８に含
まれる上述の所定のファイルに書き込み及び／又は読み
出す。つぎに、解析系３での処理について主に説明す
る。解析系３では、撮影系２で撮影された載荷系１より
荷重を受けながら変形する試験体１２の画像の解析を行
う。
【００１２】図４は、二次元応力場計測システム１００
のフローチャートである。まず、撮影系２のディジタル
ビデオカメラ２１は、載荷系１より荷重を受けた試験体
１２の表面の標点の変位分布を測定し取得する（Ｓ１０
１）。解析系３の演算部３２は、この変位分布をＩ／Ｆ
部３５を介して入力し、その変位分布から試験体１２の
表面のひずみ分布を計算するための、高精度・高分解能
の画像解析を行い、ひずみ分布を計算し、ひずみ分布フ
ァイル４０に記憶する（Ｓ１０２）。さらに、演算部３
２は、圧力計１０で計測したトラクション分布をトラク
ション分布ファイル５０に記憶する（Ｓ１０２）。解析
系３の演算部３２は、ひずみ分布ファイル４０からひず
み分布を読み出し、トラクション分布ファイル５０から
トラクション分布を読み出し、これらの各分布から応力
分布を推定・計測するための、新しい逆解析理論に基づ
く数値解析を行い、応力分布を推定・計測する（Ｓ１０
３）。
のフローチャートである。まず、撮影系２のディジタル
ビデオカメラ２１は、載荷系１より荷重を受けた試験体
１２の表面の標点の変位分布を測定し取得する（Ｓ１０
１）。解析系３の演算部３２は、この変位分布をＩ／Ｆ
部３５を介して入力し、その変位分布から試験体１２の
表面のひずみ分布を計算するための、高精度・高分解能
の画像解析を行い、ひずみ分布を計算し、ひずみ分布フ
ァイル４０に記憶する（Ｓ１０２）。さらに、演算部３
２は、圧力計１０で計測したトラクション分布をトラク
ション分布ファイル５０に記憶する（Ｓ１０２）。解析
系３の演算部３２は、ひずみ分布ファイル４０からひず
み分布を読み出し、トラクション分布ファイル５０から
トラクション分布を読み出し、これらの各分布から応力
分布を推定・計測するための、新しい逆解析理論に基づ
く数値解析を行い、応力分布を推定・計測する（Ｓ１０
３）。
【００１３】ステップＳ１０３においては、通常、撮影
系２により撮影された画像から計測されるものは変位の
みであるため、解析系３は、載荷系１より荷重を受けた
材料試験体１２の変形の様子を知ることはできるが、そ
れだけでは、材料試験としては不十分である。なぜな
ら、通常変形のみならず、応力も同時に測らなければ、
材料の特性を知ることができないためである。そこで、
本実施の形態の二次元応力場計測システム１００は、高
い空間分解能や高い精度で変形を計測することに加え、
新しく構築された逆解析理論を応用することで、応力分
布の推定・計測を可能とするものである。このため、二
次元応力場計測システム１００では、分解能や精度を極
限まで上げた高度な定量的画像解析が必要となり、解析
データ取得に見合う分解能や精度での画像解析が重要と
なる。
系２により撮影された画像から計測されるものは変位の
みであるため、解析系３は、載荷系１より荷重を受けた
材料試験体１２の変形の様子を知ることはできるが、そ
れだけでは、材料試験としては不十分である。なぜな
ら、通常変形のみならず、応力も同時に測らなければ、
材料の特性を知ることができないためである。そこで、
本実施の形態の二次元応力場計測システム１００は、高
い空間分解能や高い精度で変形を計測することに加え、
新しく構築された逆解析理論を応用することで、応力分
布の推定・計測を可能とするものである。このため、二
次元応力場計測システム１００では、分解能や精度を極
限まで上げた高度な定量的画像解析が必要となり、解析
データ取得に見合う分解能や精度での画像解析が重要と
なる。
【００１４】以下に、ステップＳ１０３での逆解析理論
とその応用について説明する。この解析系３による数値
解析は特に重要な解析であって、例えば、上述の画像解
析によって得られるデータ（後述のひずみ分布ファイル
４０、トラクション分布ファイル５０等）を利用して設
定される境界値問題を解くことで、応力分布を表す関数
を求めるものである。応力分布の境界値問題を設定する
ところに、新しく確立された逆解析理論（応力逆解析理
論）が利用される。
とその応用について説明する。この解析系３による数値
解析は特に重要な解析であって、例えば、上述の画像解
析によって得られるデータ（後述のひずみ分布ファイル
４０、トラクション分布ファイル５０等）を利用して設
定される境界値問題を解くことで、応力分布を表す関数
を求めるものである。応力分布の境界値問題を設定する
ところに、新しく確立された逆解析理論（応力逆解析理
論）が利用される。
【００１５】図５は、解析系３の応力逆解析理論に基づ
く数値解析のフローチャートである。解析系３の応力逆
解析理論による数値解析では、まず、演算部３２は、画
像解析より計測されたひずみ分布に基づくひずみ分布フ
ァイル４０を、メモリ３８から読み出し（Ｓ１０３
１）、さらに、試験体１２の縁のトラクション計測の結
果又は仮定に基づくトラクション分布ファイル５０を、
メモリ３８から読み出す（Ｓ１０３２）。
く数値解析のフローチャートである。解析系３の応力逆
解析理論による数値解析では、まず、演算部３２は、画
像解析より計測されたひずみ分布に基づくひずみ分布フ
ァイル４０を、メモリ３８から読み出し（Ｓ１０３
１）、さらに、試験体１２の縁のトラクション計測の結
果又は仮定に基づくトラクション分布ファイル５０を、
メモリ３８から読み出す（Ｓ１０３２）。
【００１６】つぎに、演算部３２は、ステップＳ１０３
１、Ｓ１０３２から読み出されたデータを基に、応力成
分に関する境界値問題を設定する（Ｓ１０３３）。演算
部３２は、ステップＳ１０３３で設定された境界値問題
を、いわゆる有限要素法と同様、弱形式を用いた定式化
と数値解法の適用によって解き、応力分布の計算を行い
（Ｓ１０３４）、この計算結果を応力分布ファイル６０
に書き込む（Ｓ１０３５）。
１、Ｓ１０３２から読み出されたデータを基に、応力成
分に関する境界値問題を設定する（Ｓ１０３３）。演算
部３２は、ステップＳ１０３３で設定された境界値問題
を、いわゆる有限要素法と同様、弱形式を用いた定式化
と数値解法の適用によって解き、応力分布の計算を行い
（Ｓ１０３４）、この計算結果を応力分布ファイル６０
に書き込む（Ｓ１０３５）。
【００１７】以下、上述のフローチャートの主なステッ
プの処理を詳細に説明する。 （境界値問題の設定：ステップＳ１０３３について）図
６は、ひずみ分布ファイル４０及び要素構成ファイル４
９の説明図である。図６（ａ）に示すように、ひずみ分
布ファイル４０は、例えば、試験体１２の標点１３の位
置に関するひずみが記憶されているファイルである。こ
のひずみ分布ファイル４０は、例えば、ポイントＩＤ４
２、位置（ｘ、ｙ）４３、ひずみ量（εｘｘ、εｙｙ、
εｘｙ）４４を含む。
プの処理を詳細に説明する。 （境界値問題の設定：ステップＳ１０３３について）図
６は、ひずみ分布ファイル４０及び要素構成ファイル４
９の説明図である。図６（ａ）に示すように、ひずみ分
布ファイル４０は、例えば、試験体１２の標点１３の位
置に関するひずみが記憶されているファイルである。こ
のひずみ分布ファイル４０は、例えば、ポイントＩＤ４
２、位置（ｘ、ｙ）４３、ひずみ量（εｘｘ、εｙｙ、
εｘｙ）４４を含む。
【００１８】また、ひずみ分布ファイル４０の代わりに
要素構成を示すファイルを用いてもよい。要素構成ファ
イル４９は、例えば、ひずみ分布ファイル４０と同様
に、試験体１２の標点１３の位置に関するひずみの種類
（ここでは、「要素構成」という概念）が記憶されてい
るファイルである（図６（ｂ））。この要素構成ファイ
ル４９は、例えば、要素番号４６、要素タイプ（例え
ば、１５節点三角形）４７、要素構成ＩＤ（例えば、１
〜１５）４８を含む。なお、要素構成ファイル４９は、
例えば、上述のステップＳ１０２で作成することになる
が、ユーザの選択により、適宜作成しないようにするこ
ともできる。
要素構成を示すファイルを用いてもよい。要素構成ファ
イル４９は、例えば、ひずみ分布ファイル４０と同様
に、試験体１２の標点１３の位置に関するひずみの種類
（ここでは、「要素構成」という概念）が記憶されてい
るファイルである（図６（ｂ））。この要素構成ファイ
ル４９は、例えば、要素番号４６、要素タイプ（例え
ば、１５節点三角形）４７、要素構成ＩＤ（例えば、１
〜１５）４８を含む。なお、要素構成ファイル４９は、
例えば、上述のステップＳ１０２で作成することになる
が、ユーザの選択により、適宜作成しないようにするこ
ともできる。
【００１９】図７は、トラクション分布ファイル５０の
説明図である。トラクション分布ファイル５０は、例え
ば、ポイントＩＤ５１、位置（ｘ、ｙ）５２、トラクシ
ョン量（ｔｘ、ｔｙ）５３を含む。トラクションとは、
例えば、上述のように、「境界に働く単位面積毎の力」
とされる。境界は、試験体１２の縁を指し、例えば、試
験体１２が四角形であれば四つの辺が境界となる。トラ
クションの計測では、例えば、縁に何も接触していない
場合、トラクションを０とし、一方、縁に何か接触して
いる場合、その接触している物体に圧力計１０（例え
ば、トランスデューサー）を当てて圧力を計測し、トラ
クションを近似的に評価する。トラクション分布ファイ
ル５０は、例えば、材料試験体１２の縁（境界）でのト
ラクションを計測又は仮定することで得られるファイル
である。なお、このトラクション分布ファイル５０に含
まれるトラクションは、応力成分の境界条件となる。
説明図である。トラクション分布ファイル５０は、例え
ば、ポイントＩＤ５１、位置（ｘ、ｙ）５２、トラクシ
ョン量（ｔｘ、ｔｙ）５３を含む。トラクションとは、
例えば、上述のように、「境界に働く単位面積毎の力」
とされる。境界は、試験体１２の縁を指し、例えば、試
験体１２が四角形であれば四つの辺が境界となる。トラ
クションの計測では、例えば、縁に何も接触していない
場合、トラクションを０とし、一方、縁に何か接触して
いる場合、その接触している物体に圧力計１０（例え
ば、トランスデューサー）を当てて圧力を計測し、トラ
クションを近似的に評価する。トラクション分布ファイ
ル５０は、例えば、材料試験体１２の縁（境界）でのト
ラクションを計測又は仮定することで得られるファイル
である。なお、このトラクション分布ファイル５０に含
まれるトラクションは、応力成分の境界条件となる。
【００２０】以下、演算部３２での計算過程を詳細に説
明する。載荷系１で荷重を加えられ、平面応力状態や平
面ひずみ状態にある材料試験体１２では、考慮すべき材
料内の応力成分は３成分、即ち、ｘｘ成分、ｙｙ成分、
ｘｙ成分のみである。この３成分は、以下の２つの釣り
合い式（１）（２）を満たしている。 δσｘｘ／δｘ＋δσｘｙ／δｙ＝０ （１） δσｘｙ／δｘ＋δσｙｙ／δｙ＝０ （２）
明する。載荷系１で荷重を加えられ、平面応力状態や平
面ひずみ状態にある材料試験体１２では、考慮すべき材
料内の応力成分は３成分、即ち、ｘｘ成分、ｙｙ成分、
ｘｙ成分のみである。この３成分は、以下の２つの釣り
合い式（１）（２）を満たしている。 δσｘｘ／δｘ＋δσｘｙ／δｙ＝０ （１） δσｘｙ／δｘ＋δσｙｙ／δｙ＝０ （２）
【００２１】ここで、２つの釣り合い式（１）（２）
は、ｘとｙを座標系として、応力成分ｘｘ成分、ｙｙ成
分、ｘｙ成分を、それぞれσｘｘ、σｙｙ、σｘｙとし
たものである。なお、σｘｘとσｙｙは、ｘ方向とｙ方
向の直応力であり、σｘｙは、せん断応力である。
は、ｘとｙを座標系として、応力成分ｘｘ成分、ｙｙ成
分、ｘｙ成分を、それぞれσｘｘ、σｙｙ、σｘｙとし
たものである。なお、σｘｘとσｙｙは、ｘ方向とｙ方
向の直応力であり、σｘｙは、せん断応力である。
【００２２】また、ひずみとこの応力３成分の間に一つ
の関係式（次式（３））が分かっていれば、ひずみを計
測し、これを解析データとすると、３つの応力成分に対
し、つりあい式とひずみの解析データの３つの式が成り
立つことになる。すなわち、十分な数の条件式が応力成
分に与えられることになる。 σｘｘ＋σｙｙ＝κ（εｘｘ￣＋εｙｙ￣） （３）
の関係式（次式（３））が分かっていれば、ひずみを計
測し、これを解析データとすると、３つの応力成分に対
し、つりあい式とひずみの解析データの３つの式が成り
立つことになる。すなわち、十分な数の条件式が応力成
分に与えられることになる。 σｘｘ＋σｙｙ＝κ（εｘｘ￣＋εｙｙ￣） （３）
【００２３】ここで、εｘｘ￣とεｙｙ￣は、ひずみの
ｘｘとｙｙ成分、κは、適宜仮定された定数（例えば、
弾性定数）である。なお、明細書中、「￣」は、直接又
は間接的に計測可能な既知の値を示すものであり、実際
には、各式中文字のそれぞれ真上に付されたものと同一
である。
ｘｘとｙｙ成分、κは、適宜仮定された定数（例えば、
弾性定数）である。なお、明細書中、「￣」は、直接又
は間接的に計測可能な既知の値を示すものであり、実際
には、各式中文字のそれぞれ真上に付されたものと同一
である。
【００２４】以上より、演算部３２は、ステップＳ１０
１、Ｓ１０２において、試験体１２のひずみ分布と縁の
トラクション分布を計測することで得られたひずみ分布
ファイル４０及びトラクション分布ファイル５０を、メ
モリ３８から読み出し（Ｓ１０３１、Ｓ１０３２）、さ
らに、「３つの条件式と境界条件」に基づいて、応力３
成分に対する境界値問題を設定する（Ｓ１０３３）。
１、Ｓ１０２において、試験体１２のひずみ分布と縁の
トラクション分布を計測することで得られたひずみ分布
ファイル４０及びトラクション分布ファイル５０を、メ
モリ３８から読み出し（Ｓ１０３１、Ｓ１０３２）、さ
らに、「３つの条件式と境界条件」に基づいて、応力３
成分に対する境界値問題を設定する（Ｓ１０３３）。
【００２５】ここで、「３つの条件式と境界条件」によ
る境界値問題の具体例について説明する。この境界値問
題は、偏微分方程式を支配方程式としており、以下、数
理的に表現する。 まず、いわゆるエイリの応力関数Ａ
は、次式のように、釣り合い状態にある応力成分を与え
る。
る境界値問題の具体例について説明する。この境界値問
題は、偏微分方程式を支配方程式としており、以下、数
理的に表現する。 まず、いわゆるエイリの応力関数Ａ
は、次式のように、釣り合い状態にある応力成分を与え
る。
【００２６】
【数６】
【００２７】つぎに、エイリの応力関数Ａに対する境界
値問題は、以下の支配方程式及び境界条件により表現さ
れる。支配方程式は、
値問題は、以下の支配方程式及び境界条件により表現さ
れる。支配方程式は、
【００２８】
【数７】
となる。境界条件は、
【００２９】
【数８】
となる。ここで、ｎｘとｎｙを、境界の単位法線ベクト
ルのｘ方向成分とｙ方向成分、κは、仮定された弾性定
数である。
ルのｘ方向成分とｙ方向成分、κは、仮定された弾性定
数である。
【００３０】また、ε￣とｒ￣は、計測されたひずみε
ｘｘ、εｙｙとトラクションで決定される次の関数であ
る。
ｘｘ、εｙｙとトラクションで決定される次の関数であ
る。
【００３１】
【数９】
【００３２】ここで、εｘｘ￣、εｙｙ￣は、直ひずみ
のｘ方向成分とｙ方向成分であり、ｔｘ￣、ｔｙ￣は、
トラクションのｘ方向成分とｙ方向成分である。なお、
ひずみと応力の関係式が線形であれば、応力３成分に対
する境界値問題は線形となる。また、ひずみと応力の関
係式が線形でない場合でも、ひずみと応力の増分の間に
は線形関係が成立すれば、応力増分の３成分に関して
は、同様に線形の境界値問題を設定することができる。
のｘ方向成分とｙ方向成分であり、ｔｘ￣、ｔｙ￣は、
トラクションのｘ方向成分とｙ方向成分である。なお、
ひずみと応力の関係式が線形であれば、応力３成分に対
する境界値問題は線形となる。また、ひずみと応力の関
係式が線形でない場合でも、ひずみと応力の増分の間に
は線形関係が成立すれば、応力増分の３成分に関して
は、同様に線形の境界値問題を設定することができる。
【００３３】（応力分布の計算：ステップＳ１０３４に
ついて）演算部３２により、ステップＳ１０３３で設定
された応力ないし応力増分に対する線形境界値問題は、
その線形性により、数値計算によって解を求めることが
できる。なお、計測データとなるひずみ分布やトラクシ
ョン分布をこの境界値問題に取り入れたり、また、線形
問題を正確に計算するには、相応の工夫が必要となる。
ついて）演算部３２により、ステップＳ１０３３で設定
された応力ないし応力増分に対する線形境界値問題は、
その線形性により、数値計算によって解を求めることが
できる。なお、計測データとなるひずみ分布やトラクシ
ョン分布をこの境界値問題に取り入れたり、また、線形
問題を正確に計算するには、相応の工夫が必要となる。
【００３４】演算部３２は、応力逆解析理論による数値
解析コードにおいて、例えば、堅牢性や適用範囲を考慮
し、いわゆる有限要素法と同様に、弱形式を用いた定式
化と、値積分や形状関数に高度な数値テクニックを適用
した解析手法で上述の境界値問題を解く。即ち、上述の
境界値問題に対し、支配方程式の弱形式の主要部（領域
内の積分）は、次のように積分表示することができる。
解析コードにおいて、例えば、堅牢性や適用範囲を考慮
し、いわゆる有限要素法と同様に、弱形式を用いた定式
化と、値積分や形状関数に高度な数値テクニックを適用
した解析手法で上述の境界値問題を解く。即ち、上述の
境界値問題に対し、支配方程式の弱形式の主要部（領域
内の積分）は、次のように積分表示することができる。
【００３５】
【数１０】
【００３６】ここで、ψは、重み関数であり、ε
ｘｘ￣、εｙｙ￣は、計測された変位のｘ方向成分とｙ
方向成分である。なお、この弱形式によれば、計測され
た変位を直接データとして用いることができる。その
後、演算部３２は、境界値問題を解くことで得られる応
力分布の計算結果を、応力分布ファイル６０に書き込
む。
ｘｘ￣、εｙｙ￣は、計測された変位のｘ方向成分とｙ
方向成分である。なお、この弱形式によれば、計測され
た変位を直接データとして用いることができる。その
後、演算部３２は、境界値問題を解くことで得られる応
力分布の計算結果を、応力分布ファイル６０に書き込
む。
【００３７】図８は、応力分布ファイル６０の説明図で
ある。応力分布ファイル６０は、例えば、境界値問題を
解くことで得られる応力分布の計算結果を示しており、
ポイントＩＤ６１、位置（ｘ、ｙ）６２、応力
（σｘ ｘ、σｙｙ、σｘｙ）６３を含む。
ある。応力分布ファイル６０は、例えば、境界値問題を
解くことで得られる応力分布の計算結果を示しており、
ポイントＩＤ６１、位置（ｘ、ｙ）６２、応力
（σｘ ｘ、σｙｙ、σｘｙ）６３を含む。
【００３８】
【発明の効果】本発明によると、以上説明した通り、材
料試験体の一様でない変形や応力の状態を計測し、破壊
直前又は破壊時の材料試験体の特性をより正しく計測す
ることができる。また、本発明によると、材料試験体の
表面の変位の分布を高精度・高分解能の画像解析によっ
て計測し、新しく確立された逆解析理論を用いて応力分
布を推定・計測することができる。また、本発明による
と、材料を限界に近い状態まで使うことができ材料利用
の効率化を図ることができる。
料試験体の一様でない変形や応力の状態を計測し、破壊
直前又は破壊時の材料試験体の特性をより正しく計測す
ることができる。また、本発明によると、材料試験体の
表面の変位の分布を高精度・高分解能の画像解析によっ
て計測し、新しく確立された逆解析理論を用いて応力分
布を推定・計測することができる。また、本発明による
と、材料を限界に近い状態まで使うことができ材料利用
の効率化を図ることができる。
【図１】二次元応力場計測システム１００の概略構成
図。
図。
【図２】試験体１２の拡大図。
【図３】ワークステーション３１の構成図。
【図４】二次元応力場計測システム１００のフローチャ
ート。
ート。
【図５】解析系３の応力逆解析理論に基づく数値解析の
フローチャート。
フローチャート。
【図６】ひずみ分布ファイル４０の説明図。
【図７】トラクション分布ファイル５０の説明図。
【図８】応力分布ファイル６０の説明図。
１ 載荷系
２ 撮影系
３ 解析系
１０ 圧力計
１１ 載荷装置
１２ 試験体
２１ ディジタルビデオカメラ
３１ ワークステーション
４０ ひずみ分布ファイル
５０ トラクション分布ファイル
６０ 応力分布ファイル
１００ 二次元応力場計測システム
Claims (10)
- 【請求項１】複数の標点を含む試験体に荷重を加え、試
験体の縁又は境界の各位置において加えられる力を示す
トラクションを測定する載荷系と、 前記載荷系により加えられた荷重による前記試験体の各
標点の変位を計測する撮影系と、 前記撮影系により計測された標点の各位置に対応する少
なくとも２成分のひずみ分布と、前記載荷系により計測
された前記試験体の縁又は境界の各位置に対応する圧力
を示す少なくとも２成分のトラクション分布とに基づ
き、前記試験体に関する所定の境界値問題を解くことに
より、標点の各位置に対応する３成分の応力分布を計算
する解析系とを備えた二次元応力場計測システム。 - 【請求項２】複数の標点を含む試験体に荷重を加え、試
験体の縁又は境界の各位置において加えられる力を示す
トラクションを測定する載荷系と、 前記載荷系により加えられた荷重による前記試験体の各
標点の変位を計測する撮影系と、 前記載荷系及び前記撮影系で計測された値に基づいて、
標点毎の応力分布を計測する解析系とを備え、 前記解析系は、さらに、 前記撮影系により計測された、標点の各位置に対応する
少なくとも２成分のひずみ量を、ひずみ分布として記憶
する第１ファイルと、 前記載荷系により計測された、前記試験体の縁又は境界
の各位置に対応する圧力を示す少なくとも２成分のトラ
クション量を、トラクション分布として記憶する第２フ
ァイルと、 前記第１ファイル及び前記第２ファイルを読み出して、
ひずみ分布及びトラクション分布に基づき設定された前
記試験体に関する境界値問題を解くことにより、標点の
各位置に対応する少なくとも３成分の応力を計算する処
理部と、 前記処理部により計算された応力を、応力分布として記
憶する第３ファイルとを有する二次元応力場計測システ
ム。 - 【請求項３】前記撮影系により計測されるひずみ量は、
試験体のｘｘ成分及びｙｙ成分の変位を含み、 前記載荷系により計測されるトラクション量は、試験体
のｘ成分及びｙ成分のトラクションを含むことを特徴と
する請求項１又は２に記載の二次元応力場計測システ
ム。 - 【請求項４】前記ひずみ分布は、試験体のｘ方向の直応
力、y方向の直応力、せん断応力を含むことを特徴とす
る請求項１乃至３のいずれかに記載の二次元応力場計測
システム。 - 【請求項５】前記解析系は、前記試験体において、次の
つり合い式及び関係式に基づき、前記トラクション量を
境界条件として境界値問題を計算するようにした請求項
１乃至４のいずれかに記載の二次元応力場計測システ
ム。 δσｘｘ／δｘ＋δσｘｙ／δｙ＝０ δσｘｙ／δｘ＋δσｙｙ／δｙ＝０ σｘｘ＋σｙｙ＝κ（εｘｘ￣＋εｙｙ￣） （ここで、σｘｘ、σｙｙ、σｘｙは、応力成分のｘｘ
成分（直応力）、ｙｙ成分（直応力）、ｘｙ成分（せん
断応力）、また、εｘｘ￣、εｙｙ￣は、ひずみのｘｘ
成分、ｙｙ成分、また、κは、仮定された弾性定数） - 【請求項６】前記解析系は、 前記境界値問題を、釣り合い状態にある応力成分を与え
る次式と、 【数１】 （ここで、Ａは、エイリの応力関数、σｘｘ、σｙｙ、
σｘｙは、応力成分のｘｘ成分（直応力）、ｙｙ成分
（直応力）、ｘｙ成分（せん断応力）） 支配方程式である次式と、 【数２】 （ここで、κは仮定された弾性定数、ε￣は、ひずみの
ｘ成分εｘｘ￣及びｙ成分εｙｙ￣に基づいて得られる
値） 境界条件である次式と 【数３】 （ここで、ｎｘとｎｙは境界の単位法線ベクトルのｘ方
向成分とｙ方向成分、ｒ￣はトラクションのｘ成分ｔｘ
￣とｙ成分ｔｙ￣に基づいて得られる値） により、境界値問題を計算するようにした請求項１乃至
５のいずれかに記載の二次元応力場計測システム。 - 【請求項７】値ε￣と値ｒ￣は、それぞれ計測されたひ
ずみ量及びトラクション量に基づいて、次式により決定
されるようにした請求項６に記載の二次元応力場計測シ
ステム。 【数４】 （ここで、ｎｘとｎｙは境界の単位法線ベクトルのｘ方
向成分とｙ方向成分、ε ｘｘ￣とεｙｙ￣はひずみのｘ
成分とｙ成分、ｔｘ￣とｔｙ￣はトラクションのｘ成分
とｙ成分、ｌは領域） - 【請求項８】前記解析系は、設定された境界値問題を、
前記支配方程式の弱形式を次式のように積分表示するこ
とにより、前記標点の変位を直接データとして用いるよ
うにした請求項６又は７に記載の二次元応力場計測シス
テム。 【数５】 （ここで、ψは、重み関数、εｘｘ￣、εｙｙ￣は、計
測された変位のｘ方向成分とｙ方向成分） - 【請求項９】複数の標点を含む試験体に荷重を加え、試
験体の縁又は境界の各位置において加えられる力を示す
トラクションを測定する載荷系と、 前記載荷系により加えられた荷重による前記試験体の各
標点の変位を計測する撮影系と、 前記載荷系及び前記撮影系で計測された値に基づいて、
標点毎の応力分布を計測する解析系とを備えた二次元応
力場計測システムに用いられる二次元応力場計測方法で
あって、 前記撮影系により計測された、標点の各位置に対応する
少なくとも３成分のひずみ量を、ひずみ分布として第１
ファイルに記憶する機能と、 前記載荷系により計測された、前記試験体の縁又は境界
の各位置に対応する圧力を示す少なくとも２成分のトラ
クション量を、トラクション分布として第２ファイルに
記憶する機能と、 前記第１ファイル及び前記第２ファイルを読み出して、
ひずみ分布及びトラクション分布に基づき設定された前
記試験体に関する境界値問題を解くことにより、標点の
各位置に対応する３成分の応力を計算し、応力分布とし
て第３ファイルに記憶する機能とを含む二次元応力場計
測方法。 - 【請求項１０】複数の標点を含む試験体に荷重を加え、
試験体の縁又は境界の各位置において加えられる力を示
すトラクションを測定する載荷系と、 前記載荷系により加えられた荷重による前記試験体の各
標点の変位を計測する撮影系と、 前記載荷系及び前記撮影系で計測された値に基づいて、
標点毎の応力分布を計測する解析系と を備えた二次元応力場計測システムに用いられる二次元
応力場計測プログラムであって、 標点の各位置に対応する少なくとも３成分のひずみ量
を、ひずみ分布として第１ファイルに記憶する手順と、 前記撮影系により計測された、前記試験体の縁又は境界
の各位置に対応する圧力を示す少なくとも２成分のトラ
クション量を、トラクション分布として第２ファイルに
記憶する手順と、 前記載荷系により計測された、前記第１ファイル及び前
記第２ファイルを読み出して、ひずみ分布及びトラクシ
ョン分布に基づき設定された前記試験体に関する境界値
問題を解くことにより、標点の各位置に対応する３成分
の応力を計算し、応力分布として第３ファイルに記憶す
る手順とをコンピュータに実行させるための二次元応力
場計測プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002032171A JP3626460B2 (ja) | 2002-02-08 | 2002-02-08 | 二次元応力場計測システム及び二次元応力場計測プログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002032171A JP3626460B2 (ja) | 2002-02-08 | 2002-02-08 | 二次元応力場計測システム及び二次元応力場計測プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003232688A true JP2003232688A (ja) | 2003-08-22 |
| JP3626460B2 JP3626460B2 (ja) | 2005-03-09 |
Family
ID=27775367
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002032171A Expired - Fee Related JP3626460B2 (ja) | 2002-02-08 | 2002-02-08 | 二次元応力場計測システム及び二次元応力場計測プログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3626460B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007510913A (ja) * | 2003-11-05 | 2007-04-26 | イノヴェイティブ サイエンティフィック ソリューションズ，インコーポレイテッド | 表面接触力を判定するための方法 |
| JP2007298343A (ja) * | 2006-04-28 | 2007-11-15 | Tokyo Institute Of Technology | 溶接変形および残留応力の評価方法 |
| JP2008224605A (ja) * | 2007-03-15 | 2008-09-25 | Sumitomo Chemical Co Ltd | 変形挙動の測定方法及び装置 |
| JP2009236754A (ja) * | 2008-03-27 | 2009-10-15 | Fukuoka Prefecture | ひずみ計測方法、ひずみ計測システム |
| JP2010085105A (ja) * | 2008-09-29 | 2010-04-15 | Toshiba Corp | 電動弁の状態監視診断装置及び方法 |
| CN107726991A (zh) * | 2017-10-09 | 2018-02-23 | 北京航空航天大学 | 一种基于弱形式的数字图像相关应变场计算方案 |
| WO2020059304A1 (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価方法 |
| WO2020059303A1 (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価システム |
| CN112986011A (zh) * | 2021-02-09 | 2021-06-18 | 北京工业大学 | 一种平面二维应力波通过节理前后的传播特性测试*** |
| CN114838852A (zh) * | 2022-05-10 | 2022-08-02 | 重庆科技学院 | 一种确定地质应力场方向的实验装置及其实验方法 |
-
2002
- 2002-02-08 JP JP2002032171A patent/JP3626460B2/ja not_active Expired - Fee Related
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007510913A (ja) * | 2003-11-05 | 2007-04-26 | イノヴェイティブ サイエンティフィック ソリューションズ，インコーポレイテッド | 表面接触力を判定するための方法 |
| JP2007298343A (ja) * | 2006-04-28 | 2007-11-15 | Tokyo Institute Of Technology | 溶接変形および残留応力の評価方法 |
| JP2008224605A (ja) * | 2007-03-15 | 2008-09-25 | Sumitomo Chemical Co Ltd | 変形挙動の測定方法及び装置 |
| JP2009236754A (ja) * | 2008-03-27 | 2009-10-15 | Fukuoka Prefecture | ひずみ計測方法、ひずみ計測システム |
| JP2010085105A (ja) * | 2008-09-29 | 2010-04-15 | Toshiba Corp | 電動弁の状態監視診断装置及び方法 |
| CN107726991A (zh) * | 2017-10-09 | 2018-02-23 | 北京航空航天大学 | 一种基于弱形式的数字图像相关应变场计算方案 |
| CN107726991B (zh) * | 2017-10-09 | 2020-05-05 | 北京航空航天大学 | 一种基于弱形式的数字图像相关应变场计算方案 |
| WO2020059303A1 (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価システム |
| JP2020046265A (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価方法 |
| JP2020046264A (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価システム |
| WO2020059304A1 (ja) * | 2018-09-18 | 2020-03-26 | 横浜ゴム株式会社 | ホースの耐疲労性評価方法 |
| US11125666B1 (en) | 2018-09-18 | 2021-09-21 | The Yokohama Rubber Co., Ltd. | Hose fatigue resistance evaluation method |
| US11243153B2 (en) | 2018-09-18 | 2022-02-08 | The Yokohama Rubber Co., Ltd. | Hose fatigue resistance evaluation system |
| JP7099216B2 (ja) | 2018-09-18 | 2022-07-12 | 横浜ゴム株式会社 | ホースの耐疲労性評価方法 |
| JP7099215B2 (ja) | 2018-09-18 | 2022-07-12 | 横浜ゴム株式会社 | ホースの耐疲労性評価システム |
| CN112986011A (zh) * | 2021-02-09 | 2021-06-18 | 北京工业大学 | 一种平面二维应力波通过节理前后的传播特性测试*** |
| CN114838852A (zh) * | 2022-05-10 | 2022-08-02 | 重庆科技学院 | 一种确定地质应力场方向的实验装置及其实验方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3626460B2 (ja) | 2005-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nicolella et al. | Machine vision photogrammetry: a technique for measurement of microstructural strain in cortical bone | |
| CN110532591B (zh) | 基于dic-efg联合仿真分析裂纹尖端应变场的方法 | |
| JP6240206B2 (ja) | 変位場およびひずみ場の測定方法、および、材料試験機 | |
| CN101369199B (zh) | 将触控面板上的坐标校正到真实显示坐标的方法 | |
| WO2011027838A1 (ja) | ひずみ計測方法、ひずみ計測装置およびプログラム | |
| US20110178728A1 (en) | Hardness test method, hardness tester, and computer-readable storage medium storing program | |
| JP5343219B2 (ja) | ひずみ計測方法、ひずみ計測システム | |
| US8204699B2 (en) | Analyzing apparatus, analyzing method, and computer-readable recording medium storing an analyzing program | |
| EP3714240B1 (en) | Weighting method and storage medium thereof | |
| JP3626460B2 (ja) | 二次元応力場計測システム及び二次元応力場計測プログラム | |
| JP3897477B2 (ja) | 応力−ひずみ関係シミュレート方法およびスプリングバック量予測方法 | |
| Neggers et al. | Improving full-field identification using progressive model enrichments | |
| Neggers et al. | Simultaneous full-field multi-experiment identification | |
| Ilg et al. | Application of a full-field calibration concept for parameter identification of HS-steel with LS-OPT® | |
| JP6464815B2 (ja) | 歪み測定方法及び装置、並びにプログラム及び記録媒体 | |
| Lecompte et al. | Parameter identification for anisotropic plasticity model using digital image correlation: Comparison between uni-axial and bi-axial tensile testing | |
| CN110795820B (zh) | 工程结构裂纹问题求解方法以及装置 | |
| KR20180019396A (ko) | 변형률비의 변동을 고려한 상사법칙을 이용한 철근콘크리트 구조물의 축소모형 유사동적실험방법 | |
| Syed-Muhammad et al. | Optimization of a mechanical test on composite plates with the virtual fields method | |
| Riegel et al. | Analysis of the stress state in QFN package during four bending experiment utilizing piezoresistive stress sensor | |
| Hagara et al. | Stress Analysis Performed in the Near Surrounding of Small Hole by a Digital Image Correlation Method | |
| WO2020255231A1 (ja) | 変位測定装置、変位測定方法、コンピュータ読み取り可能な記録媒体 | |
| JP3799419B2 (ja) | 非接触力学特性測定システム | |
| JP2017009432A (ja) | 情報処理装置、測定値の補正方法、及びプログラム | |
| Joshi et al. | The use of digital processing in studying stretch-forming sheet metal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A711 | Notification of change in applicant |
Effective date: 20031031 Free format text: JAPANESE INTERMEDIATE CODE: A712 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20031210 |
|
| A977 | Report on retrieval |
Effective date: 20040519 Free format text: JAPANESE INTERMEDIATE CODE: A971007 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Effective date: 20041109 Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Effective date: 20041202 Free format text: JAPANESE INTERMEDIATE CODE: A61 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |