(57)【特許請求の範囲】
【請求項１】
Ｘ線源と複数のＸ線検出器とを有し、被試験体の断面に関するＸ線データを収集するＸ線ＣＴ装置と、
前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理部と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像を記憶する画像記憶部と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部と、
前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算部と
を有することを特徴とする材料解析システム。
【請求項２】
前記変位量演算部が、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする請求項１記載の材料解析システム。
【請求項３】
前記変位量演算部が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする請求項１または請求項２に記載の材料解析システム。
【請求項４】
前記力学演算部が、材料の変形が微小な場合、微小変形のひずみに対応した演算式を用い、材料の変形が大きい場合、大変形のひずみに対応した演算式を用い、前記変位量に基づいて、ひずみを演算することを特徴とする請求項１から請求項３のいずれかに記載の材料解析システム。
【請求項５】
前記力学演算部が、前記変位量に基づき、材料が線形材料の場合、線形構成則による第１の演算式を用い、第１の構成マトリクスを求め、この第１の構成マトリクスに基づき、前記第１の演算式により応力を演算し、演算材料が非線形材料の場合、非線形・増分形の構成則による第２の演算式を用いて、第２の構成マトリクスを求め、この第２の構成マトリクスに基づき前記第２の演算式により応力を演算することを特徴とする請求項１から請求項４のいずれかに記載の材料解析システム。
【請求項６】
Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集過程と、
前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理過程と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶過程と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算過程と、
前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算過程と
を有することを特徴とする材料試験方法。
【請求項７】
前記変位量演算過程において、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする請求項６に記載の材料試験方法。
【請求項８】
前記変位量演算過程において、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする請求項６または請求項７に記載の材料試験方法。
【請求項９】
前記力学演算過程において、材料の変形が微小な場合、微小変形のひずみに対応した演算式を用い、材料の変形が大きい場合、大変形のひずみに対応した演算式を用い、前記変位量に基づいて、ひずみを演算することを特徴とする請求項６から請求項８のいずれかに記載の材料試験方法。
【請求項１０】
前記力学演算過程において、前記変位量に基づき、材料が線形材料の場合、線形構成則による第１の演算式を用い、第１の構成マトリクスを求め、この第１の構成マトリクスに基づき、前記第１の演算式により応力を演算し、演算材料が非線形材料の場合、非線形・増分形の構成則による第２の演算式を用いて、第２の構成マトリクスを求め、この第２の構成マトリクスに基づき前記第２の演算式により応力を演算することを特徴とする請求項６から請求項９のいずれかに記載の材料試験方法。
【請求項１１】
請求項１から請求項５のいずれかに記載の材料解析システムを動作させる材料解析プログラムであって、
Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集処理と、
前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶処理と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理と、
前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算処理と
をコンピュータに実行させることを特徴とする材料解析プログラム。
【請求項１２】
請求項１から請求項５のいずれかに記載の材料解析システムを動作させる材料解析プログラムを記録した記録媒体であって、
Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集処理と、
前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶処理と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理と、
前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算処理と
をコンピュータに実行させる材料解析プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１３】
被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加する試料固定部と、
Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線ＣＴ装置と、
前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理部と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶部と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部と、
前記変位量から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算部と
を有することを特徴とする材料試験システム。
【請求項１４】
前記試料固定部が、駆動部と、この駆動部によって回転が与えられるスプライン軸と、このスプライン軸の回転を前記各掴み具に伝える回転伝達機構とから構成され、前記上部掴み具の回転伝達機構の駆動側要素が前記スプライン軸に対してスプライン結合されていることを特徴とする請求項１３に記載の材料試験システム。
【請求項１５】
前記変位量演算部が、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする請求項１３または請求項１４に記載の材料試験システム。
【請求項１６】
前記変位量演算部が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする請求項１３〜請求項１５のいずれかに記載の材料試験システム。
【請求項１７】
被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加させる試料固定過程と、
Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集過程と、
前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理過程と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像が画像憶部に記憶される画像記憶過程と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算過程と、
前記変位量から、荷重により前記試験体の各部における応力及びひずみが演算される力学演算過程と
を有することを特徴とする材料試験方法。
【請求項１８】
前記変位量演算過程が、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする請求項１７に記載の材料試験方法。
【請求項１９】
前記変位量演算過程が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする請求項１７または請求項１８に記載の材料試験方法。
【請求項２０】
請求項１３から請求項１６のいずれかに記載の材料試験システムを動作させる材料試験プログラムであって、
被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加する試料固定処理と、
Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線データ収集処理と、
前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶処理と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部と、
前記変位量から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算処理と
をコンピュータに実行させることを特徴とする材料試験プログラム。
【請求項２１】
請求項１３から請求項１６のいずれかに記載の材料試験システムを動作させる材料試験プログラムを記録する記録媒体であって、
被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加する試料固定処理と、
Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線データ収集処理と、
前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理と、
荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶処理と、
この第１の３次元透視画像の各ピクセルを所定の演算により３次元的に移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部と、
前記変位量から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算処理と
をコンピュータに実行させる材料試験プログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、荷重が加わる前後の材料の３次元画像に基づき、この材料の各部にかかる応力及びひずみを解析する材料解析システム及び材料試験システムに係わるものである。
【０００２】
【従来の技術】
材料試験機には目的に対応して複数の種類が存在するが、一般的には、材料の試験片（被試験体）に対して、引張／圧縮／曲げなどの荷重を印加し、圧縮試験，引張試験，曲げ試験等を行えるようにした試験機が広く用いられている。
すなわち、材料試験機は、上記試験片に対して、圧縮試験，引張試験，曲げ試験等の結果に基づき、強度／延性／靱性などの力学的特性の評価を行う。
【０００３】
また、上記材料試験機は、荷重を印加する場合に、被試験体の温度を変化させたり、繰り返し試験を行うなどして、被試験体の周囲の環境を変化させて、各種の環境に対応した圧縮試験，引張試験，曲げ試験等を行う。
上述したように、材料試験機を用いて各種試験を行うことにより、人工臓器（人工骨，人工関節，人工靱帯，人工歯根，義歯など）材料，建築材料，機械材料などの各種材料の使用用途に応じた試験を行うことが可能である（例えば、非特許文献１参照）。
【０００４】
【非特許文献１】
高堂積，姫井志乃夫，有田宏志，他、Prceeding of 6th Symposium on "Microjoining and Assembly Technology in Electronics"，85-90
【０００５】
【発明が解決しようとする課題】
上述した従来の材料試験機は、被試験体に対して印加されている荷重を測定する荷重センサと、被試験体の変位を測定する変位センサとにより、応力及びひずみの定量的な測定を行っている。
すなわち、従来の材料試験は、荷重センサと変位センサとの測定値から材料の変形状態を求め、この変形状態と被試験体の形状／材質とに基づきコンピュータによる所定の演算等により、各部に発生する応力及びひずみの評価を行う。
【０００６】
しかしながら、上記材料試験機を用いた材料試験は、被試験体の材質が均一であり、かつ単純な形状なものに限られ、材質が不均一（材質の異なる材料が混在している）な場合や、被試験体が複雑な形状の場合に、印加した荷重に対応した変位−ひずみ及びひずみ−応力の関係を、この被試験体の各部（被試験体内部を含めて）においてコンピュータ等による演算から求めることができない。
また、上記材料試験機は、人体に埋め込まれる（インプラントされる）人工臓器のように、別の材料のなかに封止状態で用いられる材料について、封止された後に、封止された被試験体（例えば、人工臓器）に、周囲の別の材料から加えられる荷重を測定することができない。
【０００７】
本発明はこのような背景の下になされたもので、他の材料に埋め込まれた、材質が不均質であり、形状が複雑な被試験体の各部（内部を含めた）の荷重による変形を、周囲の構造部材の力学的特性を含めて３次元で評価し、各部におけるひずみ及び応力を定量的に測定することが可能な材料解析システムを提供する事にある。
また、本発明は、材質が不均質であり、形状が複雑な被試験体の各部（内部を含めた）において、荷重による変形を３次元で評価し、各部における応力及びひずみを定量的に測定することが可能な材料試験システムを提供する事にある。
【０００８】
【課題を解決するための手段】
本発明の材料解析システムは、Ｘ線源と複数のＸ線検出器とを有し、被試験体の断面に関するＸ線データを収集するＸ線ＣＴ装置（Ｘ線ＣＴ装置１）と、前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理部（画像処理部２）と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像を記憶する画像記憶部（画像記憶部３）と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部（変位量演算部４）と、前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算部（力学演算部５）とを有することを特徴とする。
【０００９】
本発明の材料解析システムは、前記変位量演算部が、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする。
【００１０】
本発明の材料解析システムは、前記変位量演算部が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする。
【００１１】
本発明の材料解析システムは、前記力学演算部が、前記変位量に基づき、材料の変形が微小な場合、微小変形のひずみに対応した演算式（（７）式：変位−ひずみ関係式）を用い、材料の変形が大きい場合、大変形のひずみに対応した演算式（（５）式：変位−ひずみ関係式）を用い、前記変位量に基づいて、ひずみを演算することを特徴とする。
本発明の材料解析システムは、前記力学演算部が、前記変位量に基づき、材料が線形材料（変形に線形性を示す材料）の場合、線形構成則による第１の演算式（（８）式）を用い、第１の構成マトリクス（応力とひずみとを関係付けるマトリクス、［Ｄｅ］）を求め、この第１の構成マトリクスに基づき、前記第１の演算式により応力を演算し、演算材料が非線形材料（変形に非線形性を示す材料）の場合、非線形・増分形の構成則による第２の演算式（（６）式）を用いて、第２の構成マトリクス（［Ｄｐ］）を求め、この第２の構成マトリクスに基づき前記第２の演算式により応力を演算することを特徴とする。
【００１２】
本発明の材料試験方法は、Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集過程と、前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理過程と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶過程と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算過程と、前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算過程とを有することを特徴とする。
【００１３】
本発明の材料試験方法は、前記変位量演算過程において、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする。
【００１４】
本発明の材料試験方法は、前記変位量演算過程において、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする。
【００１５】
本発明の材料試験方法は、前記力学演算過程において、前記変位量に基づき、材料の変形が微小な場合、微小変形のひずみに対応した演算式（変位−ひずみ関係式）を用い、材料の変形が大きい場合、大変形のひずみに対応した演算式（変位−ひずみ関係式）を用い、前記変位量に基づいて、ひずみを演算することを特徴とする。
本発明の材料試験方法は、前記力学演算過程において、前記変位量に基づき、材料が線形材料の場合、線形構成則による第１の演算式を用い、第１の構成マトリクスを求め、この第１の構成マトリクスに基づき、前記第１の演算式により応力を演算し、演算材料が非線形材料の場合、非線形・増分形の構成則による第２の演算式を用いて、第２の構成マトリクスを求め、この第２の構成マトリクスに基づき前記第２の演算式により応力を演算することを特徴とする。
【００１６】
本発明の材料解析プログラムは、上記材料解析システムを動作させる材料解析プログラムであって、Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集処理と、前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶処理と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理、前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算処理とをコンピュータに実行させることを特徴とする。
【００１７】
本発明のコンピュータ読み取り可能な記録媒体は、上記材料解析システムを動作させる材料解析プログラムを記録した記録媒体であって、Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集処理と、前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記被試験体の第２の３次元透視画像が画像記憶部に記憶される画像記憶処理と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理、前記変位量から、荷重による前記被試験体の各部におけるひずみ及び応力を演算する力学演算処理とをコンピュータに実行させる材料解析プログラムを記録している。
【００１８】
本発明の材料試験システムは、被試験体の上部及び下部を対向する掴み具（掴み具２１，２２）により各々把持し、この掴み具間に荷重を発生させ、この被試験体（被試験体１４）に荷重を印加する試料固定部（固定装置２０）と、Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線ＣＴ装置（Ｘ線ＣＴ装置１）と、前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理部（画像処理部２）と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶部（画像記憶部３）と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算部（変位量演算部４）と、前記係数から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算部（力学演算部５）とを有することを特徴とする。
本発明の材料試験システムは、前記試料固定部が、駆動部と、この駆動部によって回転が与えられるスプライン軸と、このスプライン軸の回転を前記各掴み具に伝える回転伝達機構とから構成され、前記上部掴み具の回転伝達機構の駆動側要素が前記スプライン軸に対してスプライン結合されていることを特徴とする。本発明の材料試験システムは、前記変位量演算部が、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする。
本発明の材料試験システムは、前記変位量演算部が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする。
本発明の材料試験方法は、被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加させる試料固定過程と、Ｘ線源と複数のＸ線検出器とを有したＸ線ＣＴ装置から、被試験体の断面に関するＸ線データを収集するデータ収集過程と、前記Ｘ線データに基づいて、被試験体の３次元透視画像が生成される画像処理過程と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像が画像憶部に記憶される画像記憶過程と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算過程と、前記変位量から、荷重により前記試験体の各部における応力及びひずみが演算される力学演算過程とを有することを特徴とする。
本発明の材料試験方法は、前記変位量演算過程において、前記所定の演算式の係数及び定数を逐次変更しつつ、各ピクセルの移動を行うことで、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることを特徴とする。
本発明の材料試験方法は、前記変位量演算部が、第２の３次元透視画像と移動後の第１の３次元透視画像とにおける各ピクセルを各々の階調度に基づき対応させて、第１の３次元透視画像と第２の３次元透視画像との形状の一致を判定することを特徴とする請求項１７または請求項１８に記載の材料試験方法。
本発明の材料試験プログラムは、上記記載の材料試験システムを動作させる材料試験プログラムであって、被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加する試料固定処理と、Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線データ収集処理と、前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶処理と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理と、前記変位量から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算処理とをコンピュータに実行させることを特徴とする。
本発明のコンピュータ読み取り可能な記録媒体は、上記記載の材料試験システムを動作させる材料試験プログラムを記録する記録媒体であって、被試験体の上部及び下部を対向する掴み具により各々把持し、この掴み具間に荷重を発生させ、この被試験体に荷重を印加する試料固定処理と、Ｘ線源と複数のＸ線検出器とを有し、前記被試験体の断面に関するＸ線データを収集するＸ線データ収集処理と、前記Ｘ線データに基づき被試験体の３次元透視画像を生成する画像処理と、荷重が印加されていない状態の前記被試験体の第１の３次元透視画像，及び所定の荷重が印加された状態の前記試験体の第２の３次元透視画像を記憶する画像記憶処理と、この第１の３次元透視画像の各ピクセルを所定の演算により移動させ、第１の３次元透視画像と第２の３次元透視画像との形状を一致させることで、各ピクセルの荷重による変位量を求める変位量演算処理と、前記変位量から、荷重により前記試験体の各部における応力及びひずみを演算する力学演算処理とをコンピュータに実行させる材料試験プログラムを記録する。
【００１９】
【発明の実施の形態】
本発明は、荷重が印加されていない状態での被試験体の３次元透視画像（以下、非荷重透視画像）と、所定の荷重が印加された状態での被試験体の３次元透視画像（以下、荷重透視画像）とを、Ｘ線ＣＴ装置からのＸ線の透過強度に基づくＸ線データから生成する。
以降、Ｘ線ＣＴ装置を用いることで得られる３次元画像により、本発明の材料解析システムの説明を行う。
しかしながら、本発明はＸ線ＣＴ装置に限らず、ＭＲＩ（磁気共鳴画像）や超音波装置による３次元透視画像を用いることも可能である。
そして、後に述べる変位量演算部４が所定の演算式を用いて非荷重透視画像の各ピクセルを３次元的（空間的）に移動させる。
【００２０】
このとき、この移動後の非荷重透視画像と上記荷重透視画像との形状が一致するまで、演算式の係数及び定数の数値を逐次変更させながらシミュレーションを行う。
すなわち、上記変位量演算部４は、移動後の非荷重透視画像と上記荷重透視画像とにおいて、同様な階調度を有するピクセルの分布形状が略一致したときに、上記荷重透視画像と各ピクセルを移動させた非荷重透視画像との形状が一致したとする。
【００２１】
そして、このとき、変位量演算部４が、各ピクセルを移動させるために用いた演算式の係数を各々抽出して、この係数に基づいて荷重透視画像における各ピクセルの変位量を求め、後に説明する力学演算部５が、この変位量に基づいて、所定の演算式から被試験体の各部にかかる応力及びひずみの定量的な数値を演算する。この結果、人間の体内などに埋め込まれる人工骨など、他の材料に埋め込まれた被試験体（内部の構造が力学的に不均質な構造体）の各部にかかるひずみ及び応力の定量的な解析を行うことが可能となる。
【００２２】
以下、図面を参照して本発明の実施形態について説明する。
＜第１の実施形態：材料試験システム＞
図１は本発明の第１の実施形態による材料解析システムの構成例を示すブロック図である。この図において、本発明の材料解析システムは、少なくとも、Ｘ線ＣＴ装置１，画像処理部２，画像記憶部３，変位量演算部４及び力学演算部５を有している。
Ｘ線ＣＴ装置１は、被試験体に対してＸ線を照射し、被試験体を透過するＸ線の強度、すなわち透過強度に基づいたＸ線データを出力する。
【００２３】
次に、図２を参照して、上記Ｘ線ＣＴ装置１の構成を簡単に説明する。図２はＸ線ＣＴ装置の構成例を示す概念図である。
Ｘ線ＣＴ装置１は、Ｘ線を放射するＸ線管１１と、このＸ線の強度を検出するＸ線検出器１３とを対向する位置に有し、被試験体１４がＸ線管１１とＸ線検出器１３との間に、図示しない固定装置により固定されている。
上記固定装置は、つかみ具などにより、図における被試験体１４の上部及び下部を固定して、引張／圧縮／ねじり等の荷重を被試験体１４に印加する。
【００２４】
この固定装置は、上記荷重を印加した状態で、被試験体１４の体軸に沿って、この被試験体をＲ方向に、３６０°の角度範囲において回転させることができる。
ここで、Ｘ線管１１が扇状に広がるＸ線１２を放射するため、固定装置は、この扇状に広がったＸ線１２が被試験体１４の横幅となるように、Ｘ線管１１及びＸ線検出器１３の位置を調整する。
そして、Ｘ線検出器１３は、例えば、被試験体１４の断面１４ｎの部分において、被試験体１４を透過した後のＸ線１２のＸ線強度を測定する。
【００２５】
また、上記固定装置は、上記体軸に沿って、Ｚ1方向またはＺ2方向に移動可能となっている。
そして、Ｘ線ＣＴ装置１は、固定装置の回転及び上下の移動を制御し、被試験体１４を回転及び移動させつつ、この被試験体１４を透過するＸ線１２のＸ線強度の測定を行う。
このとき、Ｘ線ＣＴ装置１は、被試験体１４に対し、上記荷重を印加した状態，及び印加しない状態の双方の状態において、被試験体１４を透過するＸ線１２のＸ線強度の測定を行う。
Ｘ線ＣＴ装置は、３次元透視画像が得られるものであれば、どのような構成のものでもかまわない。
【００２６】
図１に戻り、上述した構成により、Ｘ線ＣＴ装置１は、被試験体（図２における被試験体１４）の断面のＸ線データを出力する。
そして、画像処理部２は、Ｘ線ＣＴ装置１から出力される、被試験体が３６０°回転する間のＸ線データ（Ｘ線透視画像）を集積し、この集積されたＸ線データを再構成演算を行うことにより、被試験体の断面、すなわち体軸に垂直な被試験体の断層像を生成する。
【００２７】
また、画像処理部２は、体軸に沿って、同期してＺ1方向またはＺ2方向に、微小間隔に移動させて得られた上記断層像を、体軸方向に重ねて、被試験体の３次元透視画像（すなわち非荷重透視画像及び荷重透視画像）を生成する。
このとき、Ｘ線ＣＴ装置において、上述の処理により、荷重を印加していない非荷重透視画像と、所定の荷重を印加した状態の荷重透視画像とを、画像処理部２に生成させる。
さらに、画像処理部２は、得られたこれらの非荷重透視画像，荷重透視画像を画像記憶部３に記憶させるとともに、これら非荷重透視画像，荷重透視画像を、図示しない画像表示装置の表示画面に表示する。
【００２８】
以下、変位量演算部４と力学演算部５とで行われる演算の演算式及びアルゴリズムの説明を、図３を用いて行う。図３は、変位量演算部４の変位量算出の動作例を示すフローチャートである。
ステップＳ１において、画像処理部２は、画像記憶部３から上記非荷重透視画像及び上記荷重透視画像を読み出し、以降の演算に対応するように、各３次元透視画像を所定の画像サイズ（例えば、「ピクセル数／単位体積」で表される分解能）で各ピクセルの階調度を抽出し、この各ピクセルの階調度を示すピクセル値の分布で構成される３次元画像（３次元ボクセルモデル）を生成する。
ここで、非荷重透視画像から生成した３次元ボクセルモデル（以下、非荷重ボクセルモデル）において、ｘ，ｙ，及びｚの座標で示す位置にある各ピクセルの階調度をピクセル値ｐ（ｘ，ｙ，ｚ）で表し、荷重透視画像から生成した３次元ボクセルモデル（以下、荷重ボクセルモデル）において、ｘ，ｙ，及びｚの座標で示す位置にある各ピクセルの階調度をピクセル値Ｐ（ｘ，ｙ，ｚ）で表す。
【００２９】
次に、ステップＳ２において、変位量演算部４は、以下に示す（１）式，（２）式，（３）式を用いることにより、各ピクセルのｘ，ｙ及びｚ方向それぞれにおける変位量ｕ（ｘ，ｙ，ｚ），ｖ（ｘ，ｙ，ｚ），ｗ（ｘ，ｙ，ｚ）を求める。
このとき、（１），（２），（３）式で各々用いる係数ａｕｎ〜ｊｗｍ，定数ｐｕ〜ｒｗは、仮の値とする。
ここで、上記仮の値とは、すでに解析が終了した被試験体において、材質が近いと思われる材料の係数ａｕｎ〜ｊｗｍ，定数ｐｕ〜ｒｗを代用した値である。
【数１】
0003855005.tif 000002
【数２】
0003855005.tif 000003
【数３】
0003855005.tif 000004
【００３０】
この各々の式において、各大括弧のなかの第１項（中括弧の項）は変位の周期性を示すフーリエ級数であり、第２項はべき乗的な変位の量を示しており、第３項は上記定数である。
上述した各係数及び定数の仮の値とは、実験的に求められた値や、以前に同様な材料の解析によって求められた値などの、ある程度実際に近いと考えられる値が用いられる。
【００３１】
次に、ステップＳ３において、変位量演算部４は、（１），（２），（３）式から求められた変位量ｕ（ｘ，ｙ，ｚ），ｖ（ｘ，ｙ，ｚ），ｗ（ｘ，ｙ，ｚ）に基づき、各々のピクセルを移動させる。
すなわち、上記変位量ｕ，ｖ，ｗに基づき、ｐ（ｘ，ｙ，ｚ）を移動させることにより、（ｘ，ｙ，ｚ）の位置にあったピクセルは、（ｘ＋ｕ，ｙ＋ｖ，ｚ＋ｗ）の位置に移動し、すなわち、ｐ（ｘ，ｙ，ｚ）がｐ（ｘ＋ｕ，ｙ＋ｖ，ｚ＋ｗ）となる。
【００３２】
このピクセル値ｐ（ｘ＋ｕ，ｙ＋ｖ，ｚ＋ｗ）を、ピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）とする。ここで、Ｘ＝ｘ＋ｕ，Ｙ＝ｙ＋ｖ，Ｚ＝ｚ＋ｗである。
そして、変位量演算部４は、上述したピクセルの３次元座標における移動の処理を、非荷重ボクセルモデルを構成する各ピクセルに対して行う。
ここで、ピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）は、変位量演算部４により求められた仮想の荷重状態にある試料の荷重ボクセルモデル、すなわち、シミュレーション結果としての荷重ボクセルモデル（シミュレーションボクセルモデル）である。
【００３３】
次に、ステップＳ４において、変位量演算部４は、荷重ボクセルモデルにおける座標値（Ｘ，Ｙ，Ｚ）におけるピクセル値Ｐ（Ｘ，Ｙ，Ｚ）と、上記ピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）との数値の比較を、荷重ボクセルモデルにおける各座標値において行う。
すなわち、変位量演算部４は、荷重ボクセルモデルの各座標値におけるピクセル値Ｐ（Ｘ，Ｙ，Ｚ）の数値と、移動後の非荷重ボクセルモデルの各座標値におけるピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）の数値とが一致（または所定の範囲を有していて、略一致する場合も含む）するか否かの判定を行う。
【００３４】
ここで、変位量演算部４は、荷重ボクセルモデルの各座標値におけるピクセル値Ｐ（Ｘ，Ｙ，Ｚ）の数値と、移動後の非荷重ボクセルモデルの各座標値におけるピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）の数値とが一致した場合、各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗを出力して処理を終了するが、一致しない場合、処理をステップＳ５へ進める。
このとき、変位量演算部４は、下記に示す（４）式により、ボクセルモデル全体のピクセル値に対して、ピクセル値Ｐ（ｘ，ｙ，ｚ）が移動されたピクセル値Ｐ（Ｘ，Ｙ，Ｚ）と、ｐ＊（Ｘ，Ｙ，Ｚ）との差分を合計して、２５５（階調度）×Ｍにより正規化して、ｆとして出力する。
【数４】
0003855005.tif 000005
ここで、Ｍは総ボクセル数であり、Σは画像内の全ボクセルに対して総和を取ることを表している。
そして、変位量演算部４は、誤差関数ｆの誤差関数値が求められる毎に、順次前回のｆの数値と比較して、上記誤差関数値の谷（最も誤差関数値の小さい点、すなわち極小値）を検出し、この点をボクセルモデル全体の各ピクセル値の階調度の一致した点とする。
【００３５】
次に、ステップＳ５において、変位量演算部４は、各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗの変更量を、非線形計画法などの最適化手法によって、このアルゴリズに従った演算により、上記係数及び定数の数値を逐次変更する。
そして、変位量演算部４は、変更された各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗにより、新たな変位量ｕ（ｘ，ｙ，ｚ），ｖ（ｘ，ｙ，ｚ），ｗ（ｘ，ｙ，ｚ）を求めるため、処理をステップＳ２へ戻す。
【００３６】
上述したように、変位量演算部４は、荷重ボクセルモデルの各座標値におけるピクセル値Ｐ（Ｘ，Ｙ，Ｚ）の数値と、移動後の非荷重ボクセルモデルの各座標値におけるピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）の数値とが一致するまで、ステップＳ２からステップＳ５の処理を繰り返して行う。
そして、変位量演算部４は、最終的に、荷重ボクセルモデルの各座標値におけるピクセル値Ｐ（Ｘ，Ｙ，Ｚ）の数値と、移動後の非荷重ボクセルモデルの各座標値におけるピクセル値ｐ＊（Ｘ，Ｙ，Ｚ）の数値とが一致（階調度が一致）すると、このときの各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗを出力する。
【００３７】
次に、力学演算部５は、上記各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗを入力し、これらに基づき、（１），（２），（３）式により、変位量ｕ（ｘ，ｙ，ｚ），ｖ（ｘ，ｙ，ｚ），ｗ（ｘ，ｙ，ｚ）を求める。
そして力学演算部５は、上記変位量ｕ，ｖ，ｗに基づき、変位ひずみの式から、ひずみεｉｊを求める演算を行う。
【００３８】
このとき、力学演算部５は、幾何学的非線形を伴う大変形理論で対応しなければならない材料であると、以下に示す変位−ひずみの関係式である（５）式を用いてひずみεｉｊを求める。
【数５】
0003855005.tif 000006
そして、力学演算部５は、以下に示す応力−ひずみ関係を示す演算式である（６）式において、ひずみと応力との間の構成則（構成マトリクス）を与えて、単位面積当たりの力である応力σｉｊを求める演算を、荷重増分解析により行う。
【００３９】
すなわち、材料に非線形性がある場合には、応力状態に依存して変化する行列、すなわち、構成マトリクス[Ｄｐ]を用いて増分形により（６）式で与えられ、応力「０」の状態から計算を開始して、順次、応力の変化量を算出し、得られた応力を加算することにより全体の応力を求める。
【数６】
0003855005.tif 000007
ここで、（５）式は大変形理論に対応した変位−ひずみ関係式として用い、また、（６）式は非線形・増分形の演算式であり、非線形材料の応力を求めるときに用いる。
【００４０】
一方、力学演算部５は、幾何学的に線形な微小変形理論で対応する材料であると、以下に示す変位−ひずみの演算式である（７）式を用いてひずみεｉｊを求める。
【数７】
0003855005.tif 000008
そして、力学演算部５は、以下に示す応力−ひずみ関係式である（８）式において、ひずみと応力との間の構成則を与え、単位面積当たりの力である応力σｉｊを求める演算を行う。
例えば、線形材料の場合には、応力状態によらず、一定の行列、すなわち構成マトリクス[Ｄｅ]を用いた（８）式により求めることができる。
【数８】
0003855005.tif 000009
ここで（７）式は微小変形理論に対応した変位−ひずみ関係式として用い、また、（８）式は線形材料の応力を求めるときに用いる。
【００４１】
すなわち、力学演算部５は、材料（被試験体１４）の変形が微小な場合、微小変形理論に基づく、変位−ひずみ関係式に対応した（７）式を用い、材料の変形が大きい場合、大変形理論に基づく変位−ひずみ関係式に対応した（５）式を用いて、ひずみを演算する。
そして、力学演算部５は、材料（被試験体１４）が線形材料の場合、線形構成則による（８）式を用い、印加した荷重に基づく被試験体１４各部の応力を演算し、また、材料が非線形材料の場合、非線形・増分形の構成則の（６）式を用いて、印加した荷重に基づく被試験体１４各部の応力を演算する。
【００４２】
ここで、通常の材料試験では、応力あるいはひずみを簡単に算出できる単純な形状の試験片を用い、単純負荷により、荷重と変位との関係を計測し、その後に荷重から応力を換算し、変位からひずみを換算して、構成マトリクス［Ｄｐ］や［Ｄｅ］を求めている。
上述の説明では、構成マトリクスを求める流れにおいて、微小変形理論と線形材料の線形構成則とを対応させ、また大変形理論と非線形材料とを対応させて説明したが、当然、荷重の状態及び材料の材質により、微小変形理論と非線形材料との組み合わせ、または大変形理論と線形材料との組み合わせで演算する場合もある。
【００４３】
しかしながら、（６）式及び（８）式における構成マトリクス［Ｄｐ］および［Ｄｅ］は、内部が不均一な材料に対応した関係であるため、上述したように、荷重と変位とから単純に、応力−ひずみの関係を算出することはできない。
このため、以下に、材質が不均一でかつ形状が単純でない被試験体に対応する応力−ひずみ関係を求める方法の説明を、図１及び図４を用いて行う。
図４は、力学演算部５の応力−ひずみの関係を算出する動作例を説明するフローチャートである。
ステップＳ１１において、力学演算部５は、変位量演算部４から、所定の荷重に対応した被試験体１４各部の変位量ｕ（ｘ，ｙ，ｚ），ｖ（ｘ，ｙ，ｚ），ｗ（ｘ，ｙ，ｚ）を入力する。
【００４４】
次に、ステップＳ１２において、力学演算部５に対して、過去の実験例から近いと思われる材料を例として、この材料のモデルを仮に設定する。
被試験体１４が微小変形理論で扱い得る線形材料であれば、演算に用いるモデルとして、応力−ひずみ関係式を（８）式と仮定する。
一方、被試験体１４が非線形材料であれば、演算に用いるモデルとして、応力−ひずみ関係式を増分形の（６）式と仮定する。
この設定は、図示しない入力装置から、被試験体１４が線形材料または非線形材料のいずれかであるかのデータの入力を行うことにより設定する。
【００４５】
また、被試験体１４が非線形材料のとき、応力状態等による構成マトリクス［Ｄｐ］の変化を適時モデル化して、算出方法を設定する。例えば、弾塑性材料の場合には、結合流れ則に基づき設定する。
次に、ステップＳ１３において、解析対象の被試験体１４に対応する構成マトリクス［Ｄｅ］または［Ｄｐ］に含まれる各材料パラメータ（各材料のヤング率、ポアソン比等）を、未知の変数として仮に設定する。
【００４６】
次に、ステップＳ１４において、力学演算部５は、仮のパラメータが設定された構成マトリクス［Ｄｅ］または［Ｄｐ］に基づき、有限要素解析を行う。
このとき、試験対象の被試験体１４が線形材料であれば、以下に示す（９）式による演算から要素剛性マトリクス［Ｋｅ］を算出する。
そして、力学演算部５は、仮に設定された構成マトリクス［Ｄｅ］により求めた上記要素剛性マトリクス［Ｋｅ］に基づき、印加された荷重と被試験体１４の各部の変位量（以下、演算変位量）との関係を求める。
【００４７】
この印加される荷重は、ステップＳ１１において入力された所定の荷重と同等の数値が用いられる。
【数９】
0003855005.tif 000010
この（９）式において、マトリクス［Ｂ］は、変位−ひずみマトリクスであり、∫ｖｅdvは有限要素法における解析の１要素に関する体積分を表している。
次に、力学演算部５は、試験対象の被試験体１４に対して、以下に示す（１０）式による演算から、応力状態等に依存する要素剛性マトリクス［Ｋｐ］を算出する。
【００４８】
すなわち、力学演算部５は、仮に設定された構成マトリクス［Ｄｐ］により求めた上記要素剛性マトリクス［Ｋｐ］に基づき、印加された荷重と被試験体１４の各部の変位量（以下、演算変位量）との関係を、増分解析により求める。
【数１０】
0003855005.tif 000011
さらに、大変形（幾何学的非線形性）に対する考慮が必要な場合には、ひずみの定義を、以下に示す式（１１）として、接線剛性マトリクスを算出して、増分解析を行う。
【数１１】
0003855005.tif 000012
【００４９】
次に、ステップＳ１５において、力学演算部５は、有限要素法により求められた被試験体１４各部の演算変位量と、ステップＳ１１で入力された被試験体１４各部の変位量とが一致するか否かの判定を行う。
このとき、力学演算部５は、有限要素法から求めた演算変位量と、入力された変位量とが一致しない場合、処理をステップＳ１６へ進める。
【００５０】
次に、ステップＳ１６において、力学演算部５は、所定の最適化法により材料に対する構成マトリクス［Ｄｅ］または［Ｄｐ］の各パラメータの変更量を求め、新たな各パラメータを設定し、処理をステップＳ１４へ進め、再度、有限要素法により演算変位量を求める。
また、ステップＳ１５において、力学演算部５は、有限要素法から求めた演算変位量と、入力された変位量とが一致した場合、このときの構成マトリクス［Ｄｅ］または［Ｄｐ］を求める構成関係として、図示しない出力装置へ出力し、構成関係を求める処理を終了する。
【００５１】
そして、力学演算部５は、得られた構成マトリクス［Ｄｅ］または［Ｄｐ］を用いて（８）式または（６）式により、被試験体１４の各部にかかる応力の計算を行う。
すなわち、力学演算部５は、変位量演算部４の求めた変位量から、材料（被試験体１４）が線形材料の場合、線形構成則による（８）式を用い、構成関係を求め、この線形材料の構成関係に基づき、上記（８）式により被試験体１４各部の応力を演算し、材料が非線形材料の場合、非線形・増分形の構成則の（６）式を用いて、非線形材料の構成関係を求め、この非線形材料の構成関係に基づき（６）式により、被試験体１４各部の応力を演算する。
【００５２】
上記ステップＳ３において、変位量演算部４の行う各ピクセル値Ｐ＊（Ｘ，Ｙ，Ｚ）の値の求め方について、以下に説明する。
まず、変位量演算部４は、図７（ａ）に示すボクセルモデルにおいて、各画素に対応して空間格子を作成し、図７（ｂ）に示すように、各ピクセル値を対応する格子の重心に各々割り付ける。
そして、変位量演算部４は、演算した「ｕ，ｖ，ｗ」により上記重心を、図７（c）に移動させる（重心を変形写像により移動させる、（Mapping，x=D(X)））。
移動後の重心（以下、移動点（Mapped Points）：黒丸）と、変形画像を構成する空間固定格子の重心（以下、対応点（Fixed Points）：白丸）とは、図７（ｃ）で分かるように一般的には一致しない。
次に、変位量演算部４は、最終的に、各空間固定格子毎に、移動点（：変形写像）のピクセル値を、対応点に割り付けることにより、図７（ｄ）に示すように、シミュレーションボクセルモデルを生成する。
【００５３】
このため、以下の条件により、各移動点のピクセル値から、各空間固定格子の対応点のピクセル値を決定する。
基本的に、本実施形態においては、各対応点に対して最も近傍に存在する移動点のピクセル値を、この対応点のピクセル値とする、最近隣内挿法を基礎とした８近傍最近隣近似法を用いている。
すなわち、変位量演算部４は、移動点が、対応点を含む空間固定格子内へ移動したとき、移動点がこの空間固定格子内に１つの場合、移動点のピクセル値を対応点のピクセル値とし、移動点が空間固定格子内に複数ある場合、これら移動点のピクセル値の平均値を演算し、この平均値を対応点のピクセル値とする。
【００５４】
また、対応点の空間固定格子内に移動点が存在しない場合、その周囲３×３格子内の８近傍格子内の探索を行う。
このとき、周囲３×３格子内に移動点が存在しない場合、ピクセル値を「０」とし、移動点が存在する場合、その移動点のなかから対応点の位置に最も近傍にある移動点のピクセル値を、この対応点のピクセル値とする。
上述の説明は説明を簡略化するために２次元で行ったが、三次元の場合には、その空間固定格子内に移動点が無い場合、近接格子周囲３×３×３格子内の２６格子として移動点の探索を行い、８近傍格子の場合と同様の手法により対応点の格子点を決定する。
【００５５】
次に、ステップＳ５における各係数及び定数各々の数値の求め方について説明する。
すでに述べたように、ステップＳ５において、変位量演算部４は、上述したように、各係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗの変更量を、非線形計画法などの最適化手法によって、このアルゴリズに従った演算により、上記係数及び定数の数値を逐次変更する。
ここで、ボクセル値は離散値であるため、また最近隣近似に基づいて離散的な画像処理を行うため、目的関数の感度情報（各係数及び定数の変化量に対応する数値ｆの変化の度合い）を得ることができない。
【００５６】
このため、本実施形態において、変位量演算部４は、遺伝的アルゴリズム（Genetic Algorithm）等（以下説明するa)及びb)のアルゴリズム）の離散的最適化方法により誤差最小化処理を行う。
a)遺伝的アルゴリズムを用いる場合
変位量演算部４は、例えば、変数及び定数のデータ容量として、各々８ビットバイナリを割り当て、各変数及び定数毎に２５６段階の数値を有しているとする。
この２５６段階の数値は、予め変形の範囲（探索範囲として）として想定される所定の範囲（ボクセルモデルのピクセル数や寸法など）を、２５５分割したものである。
そして、遺伝的アルゴリズムのシミュレーションの条件として、各個体は、変位量演算部４において、上記各変数及び定数の数値から選択されて、１０００個の組み合わせとして構成される。
各個体の組み合わせを変更する遺伝子操作としては、一様交叉（交叉率７５％），線形ランキング選択，エリート戦略及び突然変異（突然変異率１.０％）などを用いた。
【００５７】
b)SCE-UA(Shuffled Complex Evolution Method developed at the University Arizona)法
このSCE-UA方法は、遺伝的アルゴリズム類似の進化手法を取り入れ、物理的な意味を持つ解の探索領域を設定できることと、計算の途中で、複数の個体群で独立に進化させるため、最終的に解を得る場合、計算不能や発散などが発生せず、必ず解が求まることを特徴としている。
【００５８】
すなわち、SCE-UA法は、競争進化及び集団混合の概念を組み合わせた帯域探索による最適化手法である。以下に計算手順を簡単に示す。
▲１▼ 決定変数（（１）〜（３）式における係数ａｕｎ〜ｊｗｍ及び定数ｐｕ〜ｒｗ、または（１'）〜（３'）式における変数ａｕ１３（k１,k２,k３）等）の個数ｎ，集団数ｐ，各集団内の個体数ｍ(＝２ｎ＋１)を設定し、各個体を決定変数を座標値とするｎ次元空間上の点とみなす。
例えば、係数および定数の数が合計ｎ個ある場合、ｐ個の集団に、各々「２ｎ＋１」個の個体群を振り分ける。
▲２▼ 初期世代として、各決定変数値をランダムに与え、例えば、（４）式の誤差関数の誤差関数値が低い順に個体を並べる。
【００５９】
▲３▼ ▲２▼において並べた順に、先頭の個体から、１番目の個体を第１集団に含め、２番目の個体を第２集団へ、…、ｐ番目の個体を第ｐ集団へ、ｐ＋１番目の個体を第１番目の集団へ、２ｐ(すなわち、ｐ＋ｐ)番目の個体を第ｐ集団へ、…と、全ての集団へ順次個体を振り分ける。
すなわち、この振り分け方法は、ｐ（Ｘ，Ｙ，Ｚ）と、このピクセル値に対応する各係数および定数を選択して組み合わせて演算されたｐ＊（Ｘ，Ｙ，Ｚ）とにより、誤差関数ｆを求めて、誤差関数値が小さい順にｐ個の各集団に順次振り分けていく。
例えば、係数および定数の数が合計ｎ個ある場合、ｐ個の集団に、各々「２ｎ＋１」個の個体群を振り分ける。
【００６０】
▲４▼ 各集団毎に独立させて、CCE（Conditional Class Entropy）アルゴリズムにより、個体群の進化を行わせる。
上記CCEアルゴリズムにおいては、誤差関数ｆの誤差関数値が高い個体を、より誤差関数値の低い個体とするよう、係数および定数を変化させて、各個体を進化させてゆく。
▲５▼ 最終的に、ｐ個の集団を、すなわち全集団を１つにまとめて、全集団に含まれていた全個体を、再度、誤差関数値の低い順番に並べ替えて、この全個体の内の最も優れた個体の誤差関数値が、予め設定された収束判定値の範囲内になるか、または、係数および定数の最適化処理の繰り返し回数が予め設定された繰り返し設定値となるか、のいずれかに対応するまで、▲３▼〜▲５▼を繰り返して、係数および定数の最適化処理を行う。
【００６１】
上述してきたように、本発明の材料解析システムは、３次元透視画像を取ることが可能であり、ピクセルのピクセル値が異なった値で分布している材料（構造が不均一）であれば、非荷重透視画像と荷重透視画像とを得ることにより、各不均一なひずみ及び応力を３次元的に解析することが可能である。
【００６２】
次に、図５及び図６を用いて、実際の非荷重透視画像と荷重透視画像との応力及びひずみの評価の例として、被試験体にテキスタイル材料である八つ打ちロープ（材料の形状が複雑であり、かつ材質が均一でない材料）を例にとり、３次元透視画像を用いた、荷重印加状態における被試験体内の応力及びひずみに対する評価の流れを以下に説明する。
例えば、ロープの直径は１０mmものを使用し、Ｘ線の照射条件は電圧８５ｋＶ，電流５０μＡとし、ロープ自重分の１Ｎ（ニュートン）を印加した状態（以下、非荷重状態とする）と、５０Ｎを印加した状態（以下、荷重状態）とで、Ｘ線ＣＴ装置による撮像を評価を行った。
【００６３】
図５（ｂ）は荷重を印加しない状態（上記非荷重状態）の断層像であり、図５（ａ）は非荷重状態の３次元ボクセルモデル、すなわち非荷重透視画像である。図６（ｂ）は荷重を印加した状態（上記荷重状態）の断層像であり、図６（ａ）は荷重状態の３次元ボクセルモデル、すなわち荷重透視画像である。
図５，６の双方ともに、Ｘ線の透過量が少ないほど白く表示されており、これらの図からのみでも、八つ打ちロープが荷重を上下方向に印加したことにより、対軸方向に引っ張られることで伸びてロープの繊維が詰まっていることが確認できる。
【００６４】
この例においては、図５（ｂ），図６（ｂ）の断層像を１０２４×１０２４ピクセル（ピクセルの階調度を２５６としている）とし、この断層像を０.２mm間隔（体軸方向に対して）で上下の評点間隔約１７mmの範囲で撮像を行っている。図５（ａ），図６（ａ）は、上述した評点間隔の範囲で撮像した、図５（ｂ），図６（ｂ）の複数の断層像を、上記軸方向の解像度（軸方向の撮像枚数）に対応させるため、複数のピクセルの階調度のデータの平均化処理を行い、１０２４×１０２４ピクセルから１００×１００ピクセルに、画像の解像度を低下させて変換した３次元ボクセルモデルである。
ここで、図５（ａ）が非荷重ボクセルモデルであり、図６（ａ）が荷重ボクセルモデルである。
すなわち、画像処理部２は、画像記憶部３から上記非荷重透視画像及び上記荷重透視画像を読み出し、各３次元透視画像を所定の画像サイズの１００×１００ピクセルの解像度に変更し、各ピクセルの階調度を抽出し、この各ピクセルの階調度を示すピクセル値の分布で構成される、各々の３次元ボクセルモデルを生成する。
【００６５】
そして、変位量演算部４は、図５（ａ）の非荷重ボクセルモデルにおける各ピクセルにおいて、すでにのべた（１）〜（３）式を用いたシミュレーション（図３のフローチャートの各ステップに基づいて行われる）から、変位量ｕ，ｖ，ｗを求める。
次に、力学演算部５は、八つ打ちロープが大きく変形するため、大変形理論のに基づく（５）式により、上記変位量ｕ，ｖ，ｗから、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおけるひずみを求める。
そして、力学演算部５は、八つ打ちロープが非線形材料のため、（６）式により、変位量ｕ，ｖ，ｗに基づき、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおける構成関係（応力とひずみとの関係を、すでに説明した図４のフローチャートの各ステップに従い求める。
【００６６】
次に、変位量演算部４が変位量u,v,wの算出に用いている（１），（２）及び（３）式を、下記に示す各々（１'），（２'）及び（３'）に換えて、変位量u,v,wを算出する場合について説明する。
【数１２】
0003855005.tif 000013
【数１３】
0003855005.tif 000014
【数１４】
0003855005.tif 000015
【００６７】
これら（１'），（２'）及び（３'）式は、（１），（２）,（３）式各々に対して、フーリエ級数の項に、複数の座標軸における変位の影響を含ませたものである。
（１'），（２'）及び（３'）式各々において、Ｍu，Ｍv，Ｍw及びＮu，Ｎv，Ｎwは各々整数値で問題により適宜設定する定数であり、ωは基本周波数であり、同様に問題により適宜設定する定数である。
また、αｉｊｋ，βｉｊｋ，γｉｊｋは、i,j,kが異なれば、異なる変数として定義され、ＧＡ等の最適化アルゴリズムにより、誤差関数から求められる誤差関数値が最小となるように決定される変数である。
そして、ａｕ１３(ｋ1,ｋ2,ｋ3)，ａｖ１３(ｋ1,ｋ2,ｋ3)，ａｗ１３(ｋ1,ｋ2,ｋ3)等も、ｋ1,ｋ2,ｋ3が異なれば、異なる変数として定義され、ＧＡ等の最適化アルゴリズムにより、誤差関数から求められる誤差関数値が最小となるように決定される変数である。
【００６８】
このため、（１'），（２'）及び（３'）式を用いることにより、各々の座標軸方向の変位に対して、他座標軸方向からの影響、すなわち補正量を含ませることとなり、（１），（２）,（３）式を用いた場合に比較して、変位量u,v,wを算出する際に、より精度の高いシミュレーションが行える。
例えば、材料解析の試料をプチトマトとし、荷重をかけない非荷重状態において、Ｘ線ＣＴ装置１により、１ボクセルの１辺の長さを０.４ｍｍとして、透視画像を撮像することで、非荷重透視画像が得られる。
そして、画像処理部２は、非荷重状態の試料（プチトマト）に対して、Ｘ線ＣＴ装置１により得られた非荷重透視画像から、図８（図８(a)：全体の1/2の部分、図８(b)：全体の1/8の部分）に示す初期画像の３次元ボクセルモデル（非荷重ボクセルモデル：１ボクセルの１辺の長さを０.４ｍｍ）を生成する。
【００６９】
次に、試料のプチトマトの上下を平行な板で挟み、この板に対して垂直に４００ｍＮの荷重を印加して圧縮し、Ｘ線ＣＴ装置１により、１ボクセルの１辺の長さを０.４ｍｍとして、透視画像を撮像することで、荷重透視画像が得られる。
そして、画像処理部２は、荷重状態の試料に対して、Ｘ線ＣＴ装置１により得られた荷重透視画像から、図９（図９(a)：全体の1/2の部分、図９(b)：全体の1/8の部分）に示す３次元ボクセルモデル（荷重ボクセルモデル，ボクセルサイズ：１ボクセルの１辺の長さを０.４ｍｍ）を生成する。
そして、変位量演算部４は、図３のフローチャートに従い、非荷重ボクセルモデルのピクセルを変位させて求めたｐ＊（Ｘ，Ｙ，Ｚ）と、荷重ボクセルモデル（図９）から得られるｐ（Ｘ，Ｙ，Ｚ）との比較を行い、（４）式の誤差関数が極小となる点を探索し、探索された極小点における各係数及び定数を最適値として出力する。
【００７０】
このとき得られるシミュレーション結果のシミュレーションボクセルモデル（ボクセルサイズ：１ボクセルの１辺の長さを０.４ｍｍ）を図１０（図１０(a)：全体の1/2の部分、図１０(b)：全体の1/8の部分）に示す。
また、図９の荷重ボクセルモデルと、図１０のシミュレーションボクセルモデルとにおける、各位置におけるピクセル値の差分を示す差分ボクセルモデルを図１１（全体の1/2の部分）に示す。
図１１において、差分の多い部分の階調度を高くしてあり、結果として階調度の高い部分がほぼないことから、図１１の差分値のグラフから（１'），（２'）及び（３'）式及び図３のフローチャートに基づいて、各演算により求められたシミュレーションボクセルモデルが、実際の荷重ボクセルモデルに良く一致していることが分かる。
特に、図９及び図１０とにおいて、楕円で囲んで示しているプチトマトの「へた」と「芯」の部分が非常に一致したボクセルモデルとなっていることが確認できる。
【００７１】
次に、力学演算部５は、最適値として得られたαｉｊｋ，βｉｊｋ，γｉｊｋと、ａｕ１３(ｋ1,ｋ2,ｋ3)，ａｖ１３(ｋ1,ｋ2,ｋ3)，ａｗ１３(ｋ1,ｋ2,ｋ3)等を入力して、変位量ｕ，ｖ，ｗを各々求める。
そして、力学演算部５は、求めた変位量ｕ，ｖ，ｗにより、（５）あるいは（７）式を用いて、ひずみεｉｊを算出する。
以降は、すでに述べた（１）〜（３）式を用いて行った、変位−ひずみマトリクス，要素剛性マトリクス［Ｋｐ］，構成マトリクス［Ｄｐ］，接線剛性マトリクス算出の演算と同様であるため、説明を省略する。
【００７２】
上述したように、本発明の材料解析システムを用いることにより、従来、材質が均質でなく、形状が複雑であり、材料試験において全体としての挙動が判るのみで、局所的な状態が不明である被試験体に対して、荷重を印加することにより被試験体内部の各部位に発生する応力及びひずみの分布を定量的に求めることができるため、被試験体を構造部材とする構造体の構造設計を、使用する用途に応じて精度良く行うことが可能となる。
【００７３】
また、被試験体に荷重をかける装置を具備する構成のＸ線ＣＴ装置１を用いて、本発明の構成を説明してきたが、本発明によれば、例えば、人体に埋め込む人工臓器（人工骨など）により、周囲の人体組織（構造部材と考えられる）にかかる荷重の評価を行うことが可能である。
この場合、まず、画像処理部２は、通常の医療用等に用いられるＸ線ＣＴ装置の撮像する断層像から、変形させない状態において、人工臓器を埋め込む部位の筋肉組織等の非荷重ボクセルモデル、すなわち非荷重透視画像をすでに説明した方法により作成し、画像記憶部３に格納する。
【００７４】
次に、画像処理部２は、医療用のＸ線ＣＴ装置の撮像する断層像から、所定の荷重を印加して変形させた状態において、人工臓器を埋め込む部位の筋肉組織等の荷重ボクセルモデル、すなわち荷重透視画像を、すでに説明した方法により作成し、画像記憶部３へ格納する。
このとき、被試験体である筋肉組織等の上記非荷重ボクセルモデルと上記荷重ボクセルモデルとは、同様な解像度で撮像した断層像を用い、同一の縮尺であり、かつ、ピクセル構成が同一の数、例えば双方とも１００×１００ピクセルとする。
【００７５】
そして、荷重透視画像については、所定の荷重を印加される筋肉組織等の部位を、体形を変化させながら異なった荷重がこの筋肉組織等にかかる各々の状態毎に荷重透視画像を作成する。
そして、画像処理部２は、非荷重ボクセルモデルと上記複数の荷重ボクセルモデルとの画像データを画像記憶部３へ記憶させる。
次に、変位量演算部４は、すでに述べた（１），（２），（３）式、あるいは（１'），（２'），（３'）式により、変位量ｕ，ｖ，ｗを、各々の体形の荷重ボクセルモデル毎に求める。
【００７６】
次に、力学演算部５は、荷重を印加された筋肉組織等を被試験体とする場合、変形としては大変形であるため、（５）式により変位量ｕ，ｖ，ｗに基づき、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおけるひずみを、図３のフローチャートの各ステップに従い求める。
そして、力学演算部５は、筋肉組織等が非線形材料のため、（６）式により、変位量ｕ，ｖ，ｗに基づき、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおける構成関係を、すでに説明した図４のフローチャートの各ステップに従い求める。
【００７７】
ここで、力学演算部５は、有限要素法により、被試験者の体形毎に筋肉組織等にかかるひずみ及び応力を求める。
すなわち、力学演算部５は、被試験者の複数の体形各々において、人体に埋め込まれる人工臓器により、周囲の筋肉組織等にかかる荷重を想定し、この筋肉組織等（血管や神経等も含め）に発生する応力及びひずみを演算する。
上述した演算により、各患者毎における生体組織の力学的な構造の特性、すなわち応力及びひずみの関係を得て、埋め込む人工臓器と、被試験体（上述の場合には筋肉組織等）とを構造部材とする構造体（ここでは人体）の構造設計を精度良く行うことができる。
【００７８】
そして、埋め込み周囲の筋肉組織等の応力及びひずみの関係に基づき、埋め込む人工臓器の構造設計を行い、この筋肉組織に対応する大きさや応力及びひずみ関係を有する人工臓器を作成する。
そして、この作成された人工臓器を人体に埋め込む場合に、筋肉組織にした上述の試験と同様な試験をこの人工臓器に対して行う。
【００７９】
すなわち、画像処理部２は、通常の医療用等に用いられるＸ線ＣＴ装置の撮像する断層像から、埋め込む前の人工骨の非荷重ボクセルモデル、すなわち非荷重透視画像をすでに説明した方法により作成し、画像記憶部３に格納する。
次に、画像処理部２は、医療用のＸ線ＣＴ装置の撮像する断層像から、人体に埋め込んだ後の人工骨の荷重ボクセルモデル、すなわち荷重透視画像を、すでに説明した方法により作成し、画像記憶部３へ格納する。
【００８０】
このとき、被試験体である人工骨の上記非荷重ボクセルモデルと上記荷重ボクセルモデルとは、同様な解像度で撮像した断層像を用い、同一の縮尺であり、かつ、ピクセル構成が同一の数、例えば双方とも１００×１００ピクセルとする。
そして、荷重透視画像については、人工骨を埋め込んだ人体の部位を、体形を変化させながら異なった荷重がこの人工骨にかかる各々の状態毎に荷重透視画像を作成する。
【００８１】
そして、画像処理部２は、非荷重ボクセルモデルと上記複数の荷重ボクセルモデルとの画像データを画像記憶部３へ記憶させる。
次に、変位量演算部４は、すでに述べた（１），（２），（３）式により、変位量ｕ，ｖ，ｗを、各々の体形の荷重ボクセルモデル毎に求める。
次に、力学演算部５は、人体に埋め込まれた人工骨を被試験体とする場合、変形としては微小な変形であるため、（６）式により変位量ｕ，ｖ，ｗに基づき、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおけるひずみを、図３のフローチャートの各ステップに従い求める。
そして、力学演算部５は、人工骨が線形材料のため、（７）式により、変位量ｕ，ｖ，ｗに基づき、非荷重ボクセルモデルに対応した荷重ボクセルモデルにおける応力を、すでに説明した図４のフローチャートの各ステップに従い求める。
【００８２】
ここで、力学演算部５は、非荷重ボクセルモデルの各ピクセルに対する、荷重ボクセルモデル毎の各ピクセルの変位量に基づき、被試験者の体形毎に人工骨にかかるひずみ及び応力を求める。
すなわち、力学演算部５は、被試験者の複数の体形各々において、人体の他の生体組織から人工骨にかかる荷重により、この人工骨各部（内部も含め）に発生する応力及びひずみを演算する。
そして、すでに演算されている筋肉組織等の応力及びひずみの関係に対応しているか否かの判定、すなわち、この筋肉組織等に不都合な荷重がかかっているか否かの判定を行い、人工臓器の構造設計の妥当性を検討する。
これにより、発明の材料解析システムを用いることにより、各患者の人体の個体差に合わせた人工臓器の構造設計を、他の構造部材（例えば、筋肉組織等）に対する適合性を向上させ、高い精度で設計することが可能となる。
【００８３】
上述した本発明の材料解析システムを用いることにより、従来評価できなかった人体に埋め込んだ人工臓器の評価において、人工臓器を埋め込む前に、人工臓器を埋め込むことにより周囲の生体組織にかかる応力およびひずみと、周囲の生体組織から印加される荷重によりこの人工臓器の各部にかかる応力およびひずみとを、３次元において定量的に求めることにより、構造部材としての生体組織や人工臓器の強度などに対する構造設計及び材料設計を、人体の個体差及び使用環境に対応して行うことができ、埋め込む人体の各個体差に対して適合性の高い人工臓器の製造を可能とすることができる。
【００８４】
また本発明の材料解析システムを用いることにより、健常人の臓器、例えば骨などに、外部から荷重を印加したり、体型を変えて荷重のかかり具合を調整して、３次元ボクセルモデルで、各状態の荷重に基づいてこの骨の各部に発生する応力及びひずみを、３次元において定量的に求め、この評価結果を人工骨の構造設計及び材料設計にフィードバックすることにより、より実際の骨に近い人工骨の開発・製造を行うことが可能である。
【００８５】
＜第２の実施形態：材料試験システム＞
本発明の第２の実施形態の材料試験システムは、第１の実施形態の材料解析システムの構成が用いられており、異なる部分は固定装置２０が第１の実施形態に含まれていない点である。
本発明は、荷重が印加されていない状態での被試験体の３次元透視画像（以下、被荷重透視画像）と、所定の荷重が印加された状態での被試験体の３次元透視画像（以下、荷重透視画像）とを、Ｘ線ＣＴ装置からのＸ線の透過強度に基づくＸ線データから生成する。
しかしながら、本発明はＸ線ＣＴ装置に限らず、ＭＲＩ（磁気共鳴画像）や超音波装置による３次元透視画像を用いることも可能である。
そして、後に述べる変位量演算部４が所定の演算式を用いて非荷重透視画像の各ピクセルを３次元的（空間的）に移動させる。
【００８６】
このとき、第２の実施形態の材料試験システムは、第１の実施形態の材料解析システムと同様に、上記移動後の非荷重透視画像と上記荷重透視画像との形状が一致するまで、演算式の係数及び定数を逐次変化させながらシミュレーションを行い、他の材料に埋め込まれた被試験体（内部の構造が力学的に不均質な構造体）の各部にかかる応力及びひずみの定量的な解析を行うことが可能となる。
【００８７】
以下、図面を参照して、第２の実施形態について説明する。図１２は本発明の第２の実施形態による材料試験システムの構成例を示すブロック図である。この図において、本発明の材料試験システムは、少なくとも、Ｘ線ＣＴ装置１，画像処理部２，画像記憶部３，変位量演算部４，力学演算部５及び固定装置２０を有している。
上記Ｘ線ＣＴ装置１，画像処理部２，画像記憶部３，変位量演算部４，力学演算部５各々は、第１の実施形態の図１における各構成と同様のため（同様の処理を行うため）、同一符号を付して説明を省略する。
【００８８】
次に、図２を参照して、上記Ｘ線ＣＴ装置１の構成を簡単に説明する。図２はＸ線ＣＴ装置の構成例を示す概念図である。
Ｘ線ＣＴ装置１は、Ｘ線を放射するＸ線管１１と、このＸ線の強度を検出するＸ線検出器１３とを対向する位置に有し、被試験体１４がＸ線管１１とＸ線検出器１３との間に、図示しない固定装置２０により固定されている。
上記固定装置２０は、つかみ具などにより、図における被試験体１４の上部及び下部を固定して、引張／圧縮／ねじり等の荷重を被試験体１４に印加する。
【００８９】
この固定装置２０は、上記荷重を印加した状態で、被試験体１４の体軸に沿って、この被試験体をＲ方向に、３６０°の角度範囲において回転させることができる。
ここで、Ｘ線管１１が扇状に広がるＸ線１２を放射するため、固定装置２０は、この扇状に広がったＸ線１２が被試験体１４の横幅となるように、Ｘ線管１１及びＸ線検出器１３の位置を調整する。
そして、Ｘ線検出器１３は、例えば、被試験体１４の断面１４ｎの部分において、被試験体１４を透過した後のＸ線１２のＸ線強度を測定する。
【００９０】
また、上記固定装置２０は、上記体軸に沿って、同期してＺ1方向またはＺ2方向に移動可能となっている。
そして、Ｘ線ＣＴ装置１は、この固定装置２０を、回転及び移動させつつ、被試験体１４を透過するＸ線１２のＸ線強度の測定を行  う。
このとき、Ｘ線ＣＴ装置１は、被試験体１４に対して、上記荷重を印加した状態，及び印加しない状態の双方の状態において、被試験体１４を透過するＸ線１２のＸ線強度の測定を行う。
【００９１】
上述した固定装置２０を図１３を用いて説明する。図１３は固定装置２０の構成例を示す概念図である。
被試験体１４は、掴み具２１により上部が把持され、掴み具２２により下部が把持されている。
アーム２３は掴み具２１を有し、またアーム２４は掴み具２２を有しており、固定装置２０は、掴み具２１及び掴み具２２を、これら掴み具の対向方向に沿う荷重軸に対して平行に、上下に駆動させる負荷機構を内部に有している。
また、掴み具２１及び２２は、固定装置２０内の駆動部により、同期して駆動され、掴んでいる試料を図２のＲ方向へ、被試験体１４の上記荷重軸に平行な中心軸を中心として回転させられる。
【００９２】
また、固定装置２０には、上記他の駆動部により回転が与えられるスプライン軸２５が設けられている。
そして、アーム２３及びアーム２４の掴み具２１，掴み具２２には、上記スプライン軸２５の回転を伝達される回転伝達機構が設けられている。
これにより、掴み具２１，掴み具２２は、試料を掴んだ状態において、同期してＲ方向への回転を行うこととなる。
ここで、上部にある掴み具２１の回転伝達機構の駆動側要素がスプライン軸２５に対してスプライン結合されている。
こうすることで、掴み具２１及び２２は、固定装置２０内に設けられた負荷機構により試料に荷重をかけつつ、同期した状態において、Ｒ方向へ回転させることができる。
【００９３】
したがって、掴み具２１が上記負荷機構により上部方向に所定の力で引っ張られるとき、掴み具２２が同様に負荷機構により下部方向に同一の力で引っ張られ、被試験体１４に対して引張荷重が印加されることとなる。
また、掴み具２１が負荷機構により下部方向に所定の力で押し出されるとき、掴み具２２が同様に負荷機構により上部方向に押し出され、被試験体１４に対して圧縮荷重が印加されることとなる。
【００９４】
この被試験体１４に印加される荷重は、掴み具２１に埋め込まれたロードセルにより測定される。
固定装置２０の駆動制御部は、この上記ロードセルにより測定された荷重検出結果を用いて、被試験体１４に印加される荷重が予め設定された数値となるように、上記負荷機能を制御する。
【００９５】
加えて、掴み具２１及び２２は、被試験体１４に荷重をかけた状態で、固定装置２０内の駆動部が同期させつつ、微小角度毎に角度を変化させ、この被試験体１４を３６０°回転することができる。
これにより、Ｘ線管１１から放射されるＸ線１２が、荷重をかけられた状態の被試験体１４を透過して、Ｘ線検出装置１３で検出され、Ｘ線検出装置１３から、図１２の画像処理部２において層画像を生成するために用いられるＸ線データとして出力される。
【００９６】
図１２に戻り、上述した構成により、Ｘ線ＣＴ装置１は、被試験体（図２における被試験体１４）の断面のＸ線データを出力する。
そして、画像処理部２は、Ｘ線ＣＴ装置１から出力される、被試験体が３６０°回転する間のＸ線データ（Ｘ線透視画像）を集積し、この集積されたＸ線データを再構成演算を行うことにより、被試験体の断面、すなわち体軸に垂直な被試験体の断層像を生成する。
【００９７】
また、画像処理部２は、体軸に沿って、同期してＺ1方向またはＺ2方向に、微小間隔に移動させて得られた上記断層像を、体軸方向に重ねて、被試験体の３次元透視画像（３次元ボクセルモデル）を生成する。
このとき、Ｘ線ＣＴ装置において、上述の処理により、荷重を印加していない非荷重透視画像の非荷重ボクセルモデルと、所定の荷重を印加した状態における荷重透視画像の荷重ボクセルモデルとを、画像処理部２に生成させる。
【００９８】
さらに、画像記憶部２は、得られたこれらの非荷重ボクセルモデル，荷重ボクセルモデルを画像記憶部３に記憶させるとともに、これら非荷重ボクセルモデル，荷重ボクセルモデルを、図示しない画像表示装置の表示画面に表示する。
以下、変位量演算部４と力学演算部５とで行われる変位量，応力及びひずみの算出を行う演算の演算式及びアルゴリズム（図３及び図４のフローチャートの処理）については、第１の実施形態と同様のため（同一符号の構成が同様の処理を行うため）省略する。
【００９９】
上述したように、本発明の材料試験システムを用いることにより、従来、評価が行えなかった、荷重をかけた状態における被試験体に対して、荷重をかけた状態で被試験体１４を３６０°回転させることが可能な固定装置２０により、荷重をかけた状態での被試験体１４の３次元ボクセルモデル、すなわち荷重透視画像を生成することが可能となるため、内部の構造が力学的に不均質であり、複雑な形状の被試験体１４の、荷重による３次元変形の形状を検出することが可能である。
【０１００】
また、本発明の材料試験システムを用いることにより、材質が均一でなく、形状が複雑であり、従来、破壊試験を用いることでしか強度／延性／靱性の評価が行えない被試験体の試験を行う場合にも、被試験体の破壊を行うことなく、荷重を印加することにより被試験体内部の各部位に発生する応力及びひずみの分布を定量的に、３次元において求めることができるため、被試験体の構造設計及び材料設計を使用する用途に応じて精度良く行うことが可能となる。
【０１０１】
このように、本発明の材料試験システムを用いることにより、従来、材質が均質でなく、形状が複雑であり、材料試験において全体としての挙動が判るのみで、局所的な状態が不明である被試験体に対して、荷重を印加することにより被試験体内部の各部位に発生する応力及びひずみの分布を定量的に求めることができるため、被試験体を構造部材とする構造体の構造設計を、使用する用途に応じて精度良く行うことが可能となる。
【０１０２】
なお、上記した本発明の実施形態においては、材料解析システム及び材料試験システムにおいて実行される手順をコンピュータ読取り可能な記録媒体に記録し、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより本発明の材料解析システムが実現されるものとする。ここでいうコンピュータシステムとは、ＯＳや周辺機器等のハードウアを含むものである。
【０１０３】
更に、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のシステムやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【０１０４】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【０１０５】
以上、本発明の第１及び第２の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計変更等があっても本発明に含まれる。
【０１０６】
【発明の効果】
本発明の材料解析システムによれば、他の材料に埋め込まれた（例えば、人体に埋め込まれた人工臓器）のように、材質が不均一であり、形状が複雑な被試験体の各部（内部を含めた）の変形を３次元ボクセルモデルによりで評価することができ、各部における応力及びひずみを定量的に測定することが可能となり、使用する環境、すなわち荷重のかかり具合に対応した構造設計及び材料設計を行うことができる。
【０１０７】
本発明の材料試験システムによれば、材質が不均一であり、形状が複雑な被試験体の各部（内部を含めた）の変形を３次元ボクセルモデルによりで評価することができ、各部における応力及びひずみを定量的に測定することが可能となり、使用する環境、すなわち荷重のかかり具合に対応した構造設計及び材料設計を行うことができる。
【図面の簡単な説明】
【図１】  本発明の第１の実施形態による材料解析システムの構成例を示すブロック図である。
【図２】  図１のＸ線ＣＴ装置１の構成を示す概念図である。
【図３】  図１における変位量演算部４の変位量算出の動作例を示すフローチャートである。
【図４】  図１における力学演算部５の応力算出の動作例を示すフローチャートである。
【図５】  本発明におけるＸ線ＣＴ装置１により撮像された層画像から生成された非荷重状態の八つ打ちロープの透視画像である。
【図６】  本発明におけるＸ線ＣＴ装置１により撮像された層画像から生成された荷重状態の八つ打ちロープの透視画像である。
【図７】  変位量演算部４における各ピクセル値Ｐ＊（Ｘ，Ｙ，Ｚ）の値の算出方法の説明を行う概念図である。
【図８】  Ｘ線ＣＴ装置１により撮像された断層画像から生成された試料（プチトマト）の非荷重状態での３次元ボクセルモデル（非荷重ボクセルモデル）の図である。
【図９】  Ｘ線ＣＴ装置１により撮像された断層画像から生成された試料（プチトマト）の荷重状態での３次元ボクセルモデル（荷重ボクセルモデル）の図である。
【図１０】  非荷重ボクセルモデルから、シミュレーションによって得られた荷重ボクセルモデル（１ボクセルの１辺の長さ０.４ｍｍ）の図である。
【図１１】  図９の荷重ボクセルモデルと、図１０のシミュレーションボクセルモデルとにおける、各位置におけるピクセル値の差分を示す差分ボクセルモデルの図である。
【図１２】  本発明の第２の実施形態による材料試験システムの構成例を示すブロック図である。
【図１３】  図１２における被試験試料１４を固定する固定装置２０の構成例を示す概念図である。
【符号の説明】
１  Ｘ線ＣＴ装置
２  画像処理装置
３  画像記憶部
４  変位量演算部
５  力学演算部
１１  Ｘ線管
１２  Ｘ線
１３  Ｘ線検出装置
１４  被試験体
２０  固定装置
２１，２２  掴み具
２３，２４  アーム
２５  スプライン軸
【図１】
0003855005.tif 000016
【図２】
0003855005.tif 000017
【図３】
0003855005.tif 000018
【図４】
0003855005.tif 000019
【図５】
0003855005.tif 000020
【図６】
0003855005.tif 000021
【図７】
0003855005.tif 000022
【図８】
0003855005.tif 000023
【図９】
0003855005.tif 000024
【図１０】
0003855005.tif 000025
【図１１】
0003855005.tif 000026
【図１２】
0003855005.tif 000027
【図１３】
0003855005.tif 000028