特許 5912541 対称性喪失レベルの定量的把握方法及び定量的把握装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】5912541

(24)【登録日】2016年4月8日

(45)【発行日】2016年4月27日

(54)【発明の名称】対称性喪失レベルの定量的把握方法及び定量的把握装置

(51)【国際特許分類】

   G06T 7/60 20060101AFI20160414BHJP

【ＦＩ】

   G06T7/60 180C

【請求項の数】9

【全頁数】28

(21)【出願番号】特願2012-2363(P2012-2363)

(22)【出願日】2012年1月10日

(65)【公開番号】特開2013-142976(P2013-142976A)

(43)【公開日】2013年7月22日

【審査請求日】2014年12月26日

【新規性喪失の例外の表示】特許法第３０条第１項適用 平成２３年７月１５日、（社）日本機械学会にて発行されたＭ＆Ｍ２０１１材料力学カンファレンス講演論文集のＯＳ２０２２にて発表 平成２３年７月１７日、（社）日本機械学会主催のＭ＆Ｍ２０１１材料力学カンファレンスにて発表

(73)【特許権者】

【識別番号】504136568

【氏名又は名称】国立大学法人広島大学

(74)【代理人】

【識別番号】110001427

【氏名又は名称】特許業務法人前田特許事務所

(72)【発明者】

【氏名】有尾 一郎

【審査官】 佐藤 実

(56)【参考文献】

【文献】 特開２０１１−０８６３５０（ＪＰ，Ａ）

【文献】 特開２００５−１４４４１２（ＪＰ，Ａ）

【文献】 田中玲子 外１名，システム制御理論における群論的対称性，システム／制御／情報，システム制御情報学会，２００１年 ４月１５日，第45巻 第4号，第21−29頁

【文献】 有尾一郎 外１名，ブロック対角化法の数値解析効率評価，計測自動制御学会論文集，社団法人計測自動制御学会，１９９８年 ６月３０日，第34巻 第6号，第593−600頁

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｔ ７／００ − ７／６０

(57)【特許請求の範囲】

【請求項1】

測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する工程ａと、
前記基本構造を表す行列を対角化する行列である変換フィルタを求める工程ｂと、
行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する工程ｃと、
前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する工程ｄとを備えていることを特徴とする対称性喪失レベルの定量的把握方法。

【請求項2】

請求項１に記載の対称性喪失レベルの定量的把握方法において、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とする対称性喪失レベルの定量的把握方法。

【請求項3】

請求項１に記載の対称性喪失レベルの定量的把握方法において、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とする対称性喪失レベルの定量的把握方法。

【請求項4】

測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａと、
前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂと、
行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃと、
前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する手段Ｄとを備えていることを特徴とする対称性喪失レベルの定量的把握装置。

【請求項5】

請求項４に記載の対称性喪失レベルの定量的把握装置において、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とする対称性喪失レベルの定量的把握装置。

【請求項6】

請求項４に記載の対称性喪失レベルの定量的把握装置において、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とする対称性喪失レベルの定量的把握装置。

【請求項7】

コンピュータを、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａ、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂ、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃ、及び、前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する手段Ｄとして機能させるためのプログラム。

【請求項8】

請求項７に記載のプログラムにおいて、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とするプログラム。

【請求項9】

請求項７に記載のプログラムにおいて、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とするプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、例えば周期対称性を持つ材料や構造における隠れた対称性のレベルを定量的に把握する方法及び装置に関する。

【背景技術】

【0002】

古典的な対称性の研究は、静的な現象を前提として、対象とする系に対称性が存在するか否かということに力点がおかれてきた。

【0003】

物理学では、対称性が低下することは「対称性の自発的破れ」又は単に「対称性の破れ」と呼ばれ、レイリー・ベナール対流やクエット・テイラーの実験などが知られている。

【0004】

静力学の分野における分岐階層構造の研究においては、系はその対称性の部分群に階層的に支配されていることが明らかにされ、構造解析に群の標準分解定理を取り入れた群論的分岐理論への展開が行われてきた（非特許文献１〜４参照）。

【先行技術文献】

【非特許文献】

【0005】

【非特許文献1】Ikeda,K.他、Bifurcation analysis of symmetric structures using block-diagonalization、Computer Methods in Applied Mechanics and Engineering、86(2)、pp.215-243、1991年

【非特許文献2】Ikeda,K.他、Block-diagonalization analysis of symmetric plates、Int. J. of Solids and Structures、29(22)、pp.2779-2793、1992年

【非特許文献3】Ario,I.他、Block-diagonalization method for symmetric structures with rotational displacements、J. of Structural mechanics and Earthquake Engineering、JSCE、No.489/I-27、pp.27-36、1994年

【非特許文献4】有尾一郎、群積による直方体ソリッド要素の剛性行列のブロック対角化法、土木学会応用力学論文集、Vol.1、1998年

【発明の概要】

【発明が解決しようとする課題】

【0006】

対称性の高い格子構造、例えば均一な物質の格子構造（ナノストラクチャー、マイクロストラクチャー、又は規則的に配置されたコロイド粒子等を含む）等において、格子欠陥等の周期性を持つ不完全パターンが存在する場合、そのような状態の認知や理解は単純に見えて一般に極めて複雑であり、例えば、どの程度の格子欠陥状態にあるのかを定量化することは困難である。特に、周期性を持つ不完全性を人間が認識できるのは、格子構造の自由度が非常に小さい場合に限られており、通常解析しようとする構造の自由度では、周期性を持つ不完全性を人間が認識することは不可能である。特に、複数のパターンによって不完全性が繰り返されている場合、人間の認識では対称性のレベルを理解することは困難で認知エラーを起こしやすい。

【0007】

ところが、元の完全な対称性を基準として、どの程度の不完全性が存在しているのかを定量的に把握できる方法、例えば、破壊等の進行に伴い不完全性を表す乱れパターンが変化して元の高い対称性が段階的に失われているような場合における対称性レベルを定量的に把握できる方法は皆無であった。

【0008】

前記に鑑み、本発明は、例えば周期的な配列を持つ格子構造等のように高い対称性を持つ構造において、ある格子欠陥等の乱れが生じている場合に、その乱れに起因する対称性喪失のレベルを定量的に把握できるようにすることを目的とする。

【課題を解決するための手段】

【0009】

前記の目的を達成するために、本発明に係る対称性喪失レベルの定量的把握方法は、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する工程ａと、前記基本構造を表す行列を対角化する行列である変換フィルタを求める工程ｂと、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する工程ｃとを備えている。

【0010】

また、本発明に係る対称性喪失レベルの対称性喪失レベルの定量的把握装置は、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａと、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂと、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃとを備えている。

【0011】

さらに、本発明に係るプログラムは、コンピュータを、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａ、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂ、及び、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃとして機能させる。

【0012】

すなわち、本発明は、物質内部における対称性の分岐階層構造を、破壊現象や転移現象などの構造不安定性現象を対象とする構造工学的見地から分析するために、均質場からの摂動パターンを想定して、対称性の階層的喪失の追跡を行うものである。

【0013】

また、本発明は、隠れた対称性の階層構造を確認するために、例えば、群の直積（群のテンソル積）を記述表現した上で、フーリエ級数列からなる選点直交性を利用して、場の（偏）微分支配方程式をその対称性を張る空間に変換し、既往のブロック対角化（ＢＤＭ）の状況（ブロック数、行列サイズ、非対角成分（連成項）の状態など）に基づいて、隠れた対称性を見つけ出し、対称性の階層レベルを定量的に把握するものである。

【0014】

さらに、本発明は、構造物の幾何学的対称性に基づく従来のＢＤＭとは異なり、新たに均質場や規則的に配置されている系からの摂動パターンを動的現象のある段階として捉え、周期位相構造における「隠れた対称性」の構造及び分岐の系列を見つけ出すことを主眼とするものである。すなわち本発明は、例えば、完全な直交関数列から構成される座標変換フィルターを、境界条件や力学的変形を伴う位相格子系に適用し、変換後の行列表示から、該系における階層的な対称性の喪失（対称性の喪失レベル）を得るものである．
以上に説明した本発明によると、例えば周期構造を持つような高い対称性を持つパターンの認識をコンピュターシステムによって指標化することが可能になると共に、その対称性の変化に伴うパターン崩壊の定量的な把握が各種装置やコンピュターシステムの利用によって可能になる。

【0015】

また、例えば新しい材料開発では、原子・分子の構造や粒子（ナノ粒子やマイクロ粒子を含む）が規則正しく配置された組成における配列の乱れ方や質的変化を把握することによって、物理化学的に重要で且つ本質的な現象の変化が判明する場合が多いが、本発明は、このような配列の乱れ方や質的変化を体系的に定量化することを可能にする。言い換えると、本発明は、多様な組成を持つ物質における対称性崩壊に起因する多様なパターン変化を魚群（物質）全体の動きの変化に例えると、漁網の種類（対称性の認識方法）によって、その時点での魚群の形（対称性の階層レベル）を定量的に捕獲することを可能にする。従って、本発明は、新しい材料開発のための指標化や予測等を与えることを可能にする。また、本発明の手法を用いることによって、構造解析において不変を付与し、超並列計算をも可能とする。

【0016】

さらに、本発明は、例えば材料開発分野における物質の周期微視構造の「対称性の破れ」に関する技術、方法、認識等を計算によって具現化するものであって、対称性の体系的な分類と対称性のレベルの定量化を実現するものであるが、本発明は、開発された材料の物理化学変化反応後に適用することも可能である。

【発明の効果】

【0017】

本発明によると、ブロック行列の個数から、対象物に残った対称性の数や種類が分かり、ブロック行列の大きさから対称性の大きさが分かるので、これらの情報から、例えば、対象物へ印加されている力の大きさや方向が分かる。尚、対象物として、例えば分子構造設計、材料設計、建築物などにも本発明は適用可能である。

【0018】

また、本発明は、例えば規則正しく配列された格子構造中の欠陥を正しく認識することが困難であったという課題を、幾何学的な情報に基づいて欠陥を定量的に把握することによって解決するものであり、この欠陥の度合いから対称性のレベルを指標化できるので、複雑な対称性の分類を機械的に実施できるという利点がある。

【0019】

さらに、本発明は、対称性に対する高度で体系的な群論の表現論を応用することによって、物質配列構造の隠れた対称性に基づく質的変化の予測を可能とするものであり、本発明によって、例えば、その時点でのエネルギー準位及び特異点、変形パターンの数及び種類、並びに対称性のレベル及び系統などの情報を得ることができる。すなわち、本発明には、物質配列構造の質的変化についての科学的予見を可能とし、且つ隠れた対称性の存在を体系的に把握できるという利点がある。これによって、対称性が見え隠れする事象を的確に捉えることができると共に、例えば物理変形パターンの種類などを正確に識別することができるので、例えばナノテクノロジー等の分野において、微細構造の崩壊や構造色による可視光の変化等が説明可能となる。

【図面の簡単な説明】

【0020】

【図1】図１は、実施形態に係る対称性喪失レベルの定量的把握方法のフロー図である。

【図2】図２は、実施形態に係る対称性喪失レベルの定量的把握方法を実施するための装置の概略構成を示す図である。

【図3】図３は、節点列Ｄn の対称性を示す図である。

【図4】図４は、群Ｄn ×Ｄn'の既約表現の次数を示す図である。

【図5】図５は、群積により拡張された既約表現の空間を示す図である。

【図6】図６は、位相格子構造（ｐ×ｑ＝６×６）を示す図である。

【図7】図７（ａ）は、図６に示す構造の剛性行列を示す図であり、図７（ｂ）は、図７（ａ）に示す剛性行列に変換フィルタを適用することによって得られた行列を示す図である。

【図8】図８（ａ）は、図６に示す構造に乱れが生じた構造の一例を示す図であり、図８（ｂ）は、図８（ａ）に示す構造の剛性行列を示す図であり、図８（ｃ）は、図８（ｂ）に示す剛性行列に変換フィルタを適用することによって得られた行列を示す図である。

【図9】図９（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図９（ｂ）は、図９（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図10】図１０（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１０（ｂ）は、図１０（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図11】図１１（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１１（ｂ）は、図１１（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図12】図１２（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１２（ｂ）は、図１２（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図13】図１３（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１３（ｂ）は、図１３（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図14】図１４（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１４（ｂ）は、図１４（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図15】図１５（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１５（ｂ）は、図１５（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図16】図１６（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１６（ｂ）は、図１６（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。

【図17】図１７は、図７（ｂ）、図８（ｃ）、図９（ｂ）、図１０（ｂ）、図１１（ｂ）、図１２（ｂ）、図１３（ｂ）、図１４（ｂ）、図１５（ｂ）及び図１６（ｂ）に示すブロック対角化行列の既約表現組成を示す図である。

【図18】図１８は、図７（ｂ）、図８（ｃ）、図９（ｂ）、図１０（ｂ）、図１１（ｂ）、図１２（ｂ）、図１３（ｂ）、図１４（ｂ）、図１５（ｂ）及び図１６（ｂ）に示すブロック対角化行列における階層的対称性喪失の推移を示す図である。

【図19】図１９（ａ）は、平面ｆｃｃ格子構造の面心格子構造を示す図であり、図１９（ｂ）は、図１９（ａ）に示す構造に対応する点群構成を示す図である。

【図20】図２０は、図１９（ａ）、（ｂ）に示す構造の剛性行列を示す図である。

【図21】図２１（ａ）〜（ｄ）は、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列の成分（基底成分）を調和関数により抽出した結果を示す図である。

【図22】図２２は、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列を示す図である。

【図23】図２３は、図２２に示す座標変換行列を、図２０に示す剛性行列に適用することによって得られた行列を示す図である。

【図24】図２４は、図２３に示す行列に所定の行・列操作をすることによって得られた行列を示す図である。

【図25】図２５は、図２４に示す行列を整理することによって得られたブロック対角化行列を示す図である。

【図26】図２６は、４質点結合（トラス）の構造を示す図である。

【図27】図２７は、図２６に示す構造の基本構造を示す図である。

【図28】図２８は、図２７に示す基本構造に乱れが生じた構造の一例を示す図である。

【図29】図２９は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。

【図30】図３０は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。

【図31】図３１は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。

【図32】図３２は、３質点結合（トラス）の構造を示す図である。

【図33】図３３は、図３２に示す構造の基本構造を示す図である。

【図34】図３４は、図３３に示す基本構造に乱れが生じた構造の一例を示す図である。

【図35】図３５は、図３３に示す基本構造に乱れが生じた構造の他例を示す図である。

【図36】図３６は、図３３に示す基本構造に乱れが生じた構造の他例を示す図である。

【図37】図３７（ａ）は、正方形要素（８×８）のＸＹ平面の対称性を説明するための基本構造の剛性行列を模式的に示す図であり、図３７（ｂ）は、該ＸＹ平面に存在する対称性のバリエーションを示す図であり、図３７（ｃ）は、正方形要素（８×８）に乱れが生じて連成項が出現した剛性行列を模式的に示す図である。

【図38】図３８（ａ）は、正方形要素（６×６）に所定の荷重が作用している様子を示す図であり、図３８（ｂ）は、図３８（ａ）に示す正方形要素（６×６）に生じた変形を示す図であり、図３８（ｃ）〜（ｈ）はそれぞれ、正方形要素（６×６）に存在する対称性のバリエーションを示す図である。

【図39】図３９は、図３８（ａ）に示す荷重を求めるためのベクトル場の剛性方程式を示す図である。

【図40】図４０は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる測定対象物の一例を示す斜視図である。

【図41】図４１は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。

【図42】図４２は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。

【図43】図４３は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。

【発明を実施するための形態】

【0021】

以下、本発明の一実施形態に係る対称性喪失レベルの定量的把握方法及び定量的把握装置について、図面を参照しながら説明する。

【0022】

図１は、本実施形態に係る対称性喪失レベルの定量的把握方法のフロー図である。

【0023】

まず、ステップＳ１において、測定対象物から配列情報（例えば格子点の幾何学的位置情報及び格子点同士の結合情報）を取得して行列化すると共に、該測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する。すなわち、ステップＳ１においては、物質内部における対称性の分岐階層構造を構造工学的見地から分析するために、言い換えると、対称性の階層的喪失を追跡するために、均質場からの摂動パターンを想定する。ここで、測定対象物からの配列情報の取得には、例えばＳＥＭ（走査型電子顕微鏡）写真等の撮影及び例えば格子点を目標ターゲットとしたデジタル画像解析が可能な画像認識装置を用いてもよい。

【0024】

次に、ステップＳ２において、ステップＳ１で抽出した基本構造を表す行列を対角化する行列である変換フィルタを求める。すなわち、ステップＳ２においては、隠れた対称性の階層構造を確認するために、例えば、群の直積（群のテンソル積）を記述表現した上で、フーリエ級数列からなる選点直交性を利用して、測定対象物の基本構造（元の周期構造）を対角化する変換フィルター（完全な直交関数列から構成される座標変換フィルター）を作成する。

【0025】

次に、ステップＳ３において、ステップＳ１で行列化された測定対象物の配列情報に、ステップＳ２で求めた変換フィルタを適用し、ブロック行列に分解する。すなわち、ステップＳ３においては、場（測定対象物の配列情報）の支配方程式をその対称性を張る空間に変換し、既往のブロック対角化（ＢＤＭ）の状況（ブロック数、行列サイズ、非対角成分（連成項）の状態など）に基づいて、隠れた対称性を見つけ出し、対称性の階層レベルを定量的に把握する。

【0026】

図２は、図１に示す本実施形態に係る対称性喪失レベルの定量的把握方法を実施するための装置（コンピュターシステム）の概略構成を示す図である。

【0027】

図２に示すように、例えば、演算手段としてのＣＰＵ（中央演算装置）１、メモリ等の記憶装置２、及びＳＥＭ写真等を撮影可能な画像認識装置３等がネットワーク４を介して接続されることにより、コンピュターシステムが構成されている。ここで、前述のステップＳ１〜Ｓ３の処理は、図２に示すコンピュターシステムをステップＳ１〜Ｓ３の各処理を実施する手段として機能させるプログラムによって行われる。

【0028】

以上に説明した本実施形態によると、ステップＳ３で得られたブロック行列の個数から、測定対象物に残った対称性の数が分かると共に、ステップＳ３で得られたブロック行列の大きさから、測定対象物に残った対称性の大きさが分かるので、これらの情報から、例えば、測定対象物へ印加されている力の大きさや方向を知ることができる。

【0029】

尚、本実施形態において、測定対象物の種類は特に限定されるものではないが、例えば分子構造設計、材料設計、建築物などを測定対象物とすることも可能である。

【0030】

また、本実施形態によると、例えば規則正しく配列された格子構造中の欠陥を正しく認識することが困難であったという課題を、幾何学的な情報に基づいて欠陥を定量的に把握することによって解決できると共に、この欠陥の度合いから対称性のレベルを指標化できるので、複雑な対称性の分類を機械的に実施することを可能にする。

【0031】

さらに、本実施形態によると、対称性に対する高度で体系的な群論の表現論を応用することによって、物質配列構造の隠れた対称性に基づく質的変化の予測を行うので、例えば、その時点でのエネルギー準位及び特異点、変形パターンの数及び種類、並びに対称性のレベル及び系統などの情報を得ることができる。すなわち、本実施形態には、物質配列構造の質的変化についての科学的予見を可能とし、且つ隠れた対称性の存在を体系的に把握できるという利点がある。これによって、対称性が見え隠れする事象を的確に捉えることができると共に、例えば物理変形パターンの種類などを正確に識別することができるので、例えばナノテクノロジーや材料設計などの分子・原子配列に関係する分野において、微細構造の崩壊や構造色による可視光の変化等が説明可能となる。

【0032】

以下、本実施形態で利用している各種物理数学的手法及びそれを用いた各種計算例について、詳細に説明する。

【0033】

（直交関数列と群対称性）
まず、ノルム完備な空間を満足する直交関数列を定義した上で、位相空間における該関数列の対称性を利用して、該関数列を間接的に変位関数に適合させる方法について、説明する。

【0034】

変位関数（直交関数列）上の節点列に対する有限な対称性として、例えば図３に示すような、周期軌道上の変位関数（三角関数列）が離散化された節点列の対称性、具体的には、三角関数が振舞う節点変位の対称群（節点列）Ｄn について考える。節点列Ｄn の全長をｎ等分割すると、正規長さｘ／Ｌ（０≦ｘ≦２Ｌ）は、節点並びによって２ｊ／ｎ（ｊ＝１、・・・、ｎ）と置換される。ここで、節点列Ｄn に対応する直交関数列は、

・・・（式１）
と表される。また、１次既約表現（１，３）_Ｄn に対応する関数列は、

・・・（式２）
と表される。このような節点列に対応する有限群の対称性は、変位関数の直交空間を求める際の組成となる。

【0035】

（位相格子空間の対称性）
群の構造については、群の直積分解を利用することにより、いくつかの群の組成を知ることができる。以下、群の直積分解を用いて、位相平面上にある支配方程式の「隠れた対称性」を把握する方法について説明する。

【0036】

以下、簡単のため、群積について、二面体群の直積Ｇ＝Ｄn ×Ｄn'に限定して説明する。群Ｇの既約表現全体は、

・・・（式３）
と与えられる。ここで、μは群Ｄn の既約表現を、μ’は別の群Ｄn'の既約表現をそれぞれ表す。群Ｄn ×Ｄn'の既約表現（μ，μ’）の次数を図４に示す。このＤn ×Ｄn'不変系の空間Ｘは、

・・・（式４）
のように拡張される。ここで、Ｘ^（μ）、Ｘ^（μ’）は群Ｄn 、Ｄn'の既約表現の空間を表す。２つの群Ｄn 、Ｄn'からなる群元｛ｇ∈Ｄn ｜ｇ_１，・・・，ｇ_ｎ｝、｛ｇ’∈Ｄn'｜ｇ_１’，・・・，ｇ_ｎ’｝の積は、

・・・（式５）
と表される。また、この群元の積からなる表現行列は既約表現μ，μ’の表現行列の直積

・・・（式６）
と表される。この表現行列を用いて同値変換

・・・（式７）
が可能となり、この座標変換行列Ｈの成分は基底関数φ^（μ）、ψ^（μ’）から構成される。

【0037】

以下、この有限群の直積Ｄn ×Ｄn'の組成を前述の「直交関数列と群対称性」で述べた直交関数列を用いて表す。ここで、この拡張された既約表現の空間を図５に示す。２つの群Ｄn 、Ｄn'の基底関数となるφ^（μ）、ψ^（μ’）がそれぞれ、｛φ^（μ）｜φ^{（１，ｉ）Ｄn}，・・・，φ^{（2，ｉ）Ｄn}，・・・｝、｛ψ^（μ’）｜ψ^{（１，j）Ｄn'}，・・・，ψ^{（2，j）Ｄn'}，・・・｝と表されるとき、φ^（μ）、ψ^（μ’）の直積は、

・・・（式８）
と表される。この基底関数の直積を各節点に対応させた成分からなる既約表現（μ，μ’）の座標変換行列は、

・・・（式９）
となる。

【0038】

ここで、このＨはユニタリーな行列であり、Ｈ^ＴＨ＝Ｉと単位行列となる。２変数関数ｗ（ｘ，ｙ）を未知関数とする偏微分方程式Ｆ（ｗ）と、（式５）の群元の積からなる表現行列Ｔとの間で、

・・・（式１０）
の群同変条件が成り立つとき、ＢＤＭが可能となる。このＢＤＭの座標変換によって既約表現毎にブロック分解すると、

・・・（式１１）
と変換される。この同値変換によって関数ｗ（ｘ，ｙ）の隠れた対称性を調べることができる。

【0039】

以下、３次元空間にも適用可能ではあるが、簡単のため、格子の例として、２次元の節点列の対称性を利用して、有限群の直積によって格子上の関数の対称性を表現する。ここで、完全な格子の対称性を記述するために、周期的な境界を持たせた位相空間（左右及び上下の端節点を巡回させた位相構造）を考えた上で、これをすべて固有の対称性を張る空間へ変換後に境界条件を導入するものとする。すなわち、境界条件は２次的要因による制約された対称条件として取り扱う。

【0040】

周期２πに正規化された位相領域

・・・（式１２）
に平面格子を導入する。ここで、Ｌ_ｘ、Ｌ_ｙはｘ，ｙ方向の長さを表す。平面上の格子点（ｐ，ｑ）（ｐ＝１，・・・，ｎ，ｑ＝１，・・・，ｎ’）を設定すると、格子間隔はｈ＝２π／ｎ、ｈ’＝２π／ｎ’となり、ｎ，ｎ’はそれぞれｘ，ｙ方向の分割数である。従って、離散化された格子座標は、

・・・（式１３）
と設定できる。格子座標（ｐｈ，ｑｈ’）の任意関数が二重の直交関数列
ｗ_ｐ，ｑ≡ｗ（ｐｈ，ｑｈ’）＝ｕ（２πｐ／ｎ）×ｕ（２πｑ／ｎ’）
と表されるとき、（式１）より、

・・・（式１４）
と分解できる。ここで、ａ^{（μ，μ’）}は既約表現（μ，μ’）に対する未定係数である。例えば、（式１４）の左辺が既知の関数であれば、（ｐ，ｑ）に対する支配方程式は、

と表される。ここで、ｗ_０，ｑ＝ｗ_ｎ，ｑ、ｗ_ｐ，０＝ｗ_ｐ，ｎ’である。このように全ての格子点について方程式を組み立てれば、

のような連立方程式を得る。ここで、

である。従って、未定係数ａ^{（μ，μ’）}は例えばガウスの消去法などの数値計算によって求めることが可能である。また、行列Ａは方程式Ｆ（ｗ）中の関数ｗの対称性（直交関数列）に依存する。例えば、関数ｗが一定値や二重フーリエ級数列であれば、Ａは対角行列となり、関数ｗが非対称であれば、Ａは密で連成された行列となる。

【0041】

（「隠れた対称性」の数値計算）
前述のように、均質場に「隠れた対称性」が存在する場合があると共に、均質領域が分岐に伴い周期的なパターンを形成しながら乱れていく場合がある。以下、均質な位相格子モデルにおける剛性分布パターンと位相対称性との関係において、剛性方程式のブロック対角化行列Ｋ’の直交性によって「隠れた対称性」を数値的に把握できることについて説明する。尚、以下の説明においては、均質場の離散化を伴う要素解析は場の位相構造に依存し、該位相構造から周期パターン変形が形成されるものとして、その階層的な対称性を探るものとする。具体的には、均等に分割された格子上の８隣接の部材剛性をモデル化し、均質な剛性分布を持つ構造の面外変位を自由度（１節点１自由度）として有限要素解析を行う。

【0042】

まず、前述のステップＳ１の処理により、図６に示す閉じた位相格子空間（ｐ×ｑ＝６×６）の構造から、図７（ａ）に示すような帯行列化されたＫ（剛性行列）を形成する。ここで、図７（ａ）において、＋は正値を、−は負値を、・はゼロ成分をそれぞれ表す（以下の他の図面において同様）。

【0043】

次に、前述のステップＳ２の処理により求めた変換フィルタを適用し、前述のステップＳ３の処理によって、図７（ａ）に示す剛性行列Ｋを座標変換して、図７（ｂ）に示すＫ’（パターンＰ_０）を得る。ここで、Ｋ’は完全に対角行列となり、最初の対角成分以外はそれぞれ直交し異なる正値となった。すなわち、図６に示す構造の剛性行列が線形独立な既約表現と一致し、この離散構造では最も対称位数が高い構造となった。図７（ｂ）から明らかなように、この対称性は、１次元のＤ_３６不変配下の群と等価な群を持つ。尚、図７（ｂ）に示すＫ’の最初の対角成分０は剛体変形モードに対応する。

【0044】

以上に説明したように、ＢＤＭによるブロックの数が対称性の高低を表すことを利用し、対称性の階層構造を調べることができる。また、以上のような位相格子構造となる大規模な離散系では対称性の利用価値が高くなる。従って、均質場の分岐階層構造に、パターン化された位相格子構造の群論的な変換則を用いることによって、該構造が持つ隠れた対称性及び階層的な対称性の仕組みを明らかにすることが可能となる。

【0045】

例えば、図８（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の一例（チェッカーボードパターン）を示している。図８（ａ）に示す剛性分布からは、前述のステップＳ１の処理により、図８（ｂ）に示すような剛性行列が得られる。また、前述のステップＳ２の処理により求めた変換フィルタを適用し、前述のステップＳ３の処理によって、図８（ｂ）に示す剛性行列を座標変換すると、図８（ｃ）に示すブロック対角化行列（パターンＰ）が得られる。

【0046】

また、図９（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図９（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図９（ｂ）に示すブロック対角化行列（パターンＨ）が得られる。

【0047】

また、図１０（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１０（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１０（ｂ）に示すブロック対角化行列（パターンＡ）が得られる。

【0048】

また、図１１（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１１（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１１（ｂ）に示すブロック対角化行列（パターンＢ１）が得られる。

【0049】

また、図１２（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１２（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１２（ｂ）に示すブロック対角化行列（パターンＢ２）が得られる。

【0050】

また、図１３（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１３（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１３（ｂ）に示すブロック対角化行列（パターンＢ’２）が得られる。

【0051】

また、図１４（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１４（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１４（ｂ）に示すブロック対角化行列（パターンＣ）が得られる。

【0052】

また、図１５（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１５（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１５（ｂ）に示すブロック対角化行列（パターンＤ）が得られる。

【0053】

また、図１６（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１６（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１６（ｂ）に示すブロック対角化行列（パターンＥ）が得られる。

【0054】

図１７は、前述の各ブロック対角化行列（パターンＰ_０、Ｐ、Ｈ、Ａ、Ｂ１、Ｂ２、Ｂ’２、Ｃ、Ｄ、Ｅ）の既約表現組成を示しており、図１８は、前述の各ブロック対角化行列（パターンＰ_０、Ｐ、Ｈ、Ａ、Ｂ１、Ｂ２、Ｂ’２、Ｃ、Ｄ、Ｅ）における階層的対称性喪失の推移を示している。

【0055】

（平面ｆｃｃ格子構造のブロック対角化）
図１９（ａ）は、平面ｆｃｃ格子構造の面心格子構造を示し、図１９（ｂ）は、図１９（ａ）に示す構造に対応する点群構成を示している。図１９（ａ）、（ｂ）において、数字は格子番号を示している。尚、以下の説明では、１節点１自由度のｎ，ｍ＝４の周期構造を前提とする。

【0056】

図１９（ａ）、（ｂ）に示す構造からは、前述のステップＳ１の処理により、図２０に示すような剛性行列（帯行列）が得られる。ここで、格子点の構成種類（○、●）から、該構造の幾何学的な関係を示す剛性行列は、

と表される。

【0057】

一方、群の直積から得られる、白節点格子点に対応する座標変換Ｈ_ｅは、中点に位置する黒節点とも点群では同型であるので、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列は、

と表される。

【0058】

図２１（ａ）〜（ｄ）は、前記座標変換行列の成分（基底成分）を調和関数により抽出した結果を示す。この基底成分を用いて、前述のステップＳ２の処理により、図２２に示す座標変換行列（変換フィルタ）が得られる。そして、前述のステップＳ３の処理により、つまり、図２２に示す座標変換行列を、図２０に示す剛性行列に適用することによって、図２３に示す行列が得られる。図２３に示す行列に所定の行・列操作をすると、図２４に示す行列が得られ、該行列をさらに整理すると、図２５に示すブロック対角化行列が得られる。

【0059】

（４質点結合構造のブロック対角化）
図２６は、４質点結合（トラス）の構造を示している。図２６において、数字は質点番号を、ｋは質点同士の結合状態（剛性係数）を、ｗは変位を示している。尚、以下の説明では、１質点１自由度の周期境界条件を前提とする。

【0060】

図２７は、図２６に示す構造の基本構造を示す図である。図２７に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_１が得られる。ここで、ｋは剛性係数を表す（以下の他の行列において同じ）。また、前述のステップＳ２の処理により、

に示す座標変換行列（変換フィルタ）Ｈが得られる。そして、前述のステップＳ３の処理により、つまり、変換フィルタＨを剛性行列Ｋ_１に適用することによって、

となり、完全に対角化される。

【0061】

また、図２８は、図２７に示す基本構造に乱れが生じた構造の一例を示す。図２８に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_２が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_２に適用することによって、

となり、完全に対角化される。

【0062】

また、図２９は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図２９に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_３が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_３に適用することによって、

となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。

【0063】

また、図３０は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図３０に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_４が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_４に適用することによって、

となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。

【0064】

さらに、図３１は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図３１に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_５が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_５に適用することによって、

となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。

【0065】

（３質点結合構造のブロック対角化）
図３２は、３質点結合（トラス）の構造を示している。図３２において、ａ、ｂ、ｃは質点番号を、ｋは質点同士の結合状態（剛性係数）をしている。尚、以下の説明では、１質点１自由度の周期境界条件を前提とする。

【0066】

図３３は、図３２に示す構造の基本構造を示す図である。尚、図３３において、数字は質点番号を、ｋは質点同士の結合状態（剛性係数）を、ｗは変位を示している。図３３に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_１が得られる。また、前述のステップＳ２の処理により、

に示す座標変換行列（変換フィルタ）Ｈが得られる。そして、前述のステップＳ３の処理により、つまり、変換フィルタＨを剛性行列Ｋ_１に適用することによって、

となり、完全に対角化される。

【0067】

また、図３４は、図３３に示す基本構造に乱れが生じた構造の一例を示す。図３４に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_２が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_２に適用することによって、

となり、完全に対角化される。

【0068】

また、図３５は、図３３に示す基本構造に乱れが生じた構造の一例を示す。図３５に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_３が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_３に適用することによって、

となり、完全に対角化される。

【0069】

また、図３６は、図３３に示す基本構造に乱れが生じた構造の他例を示す。図３６に示す構造からは、前述のステップＳ１の処理により、

に示す剛性行列Ｋ_４が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ_４に適用することによって、

となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。

【0070】

（非対角成分の分析）
以下、本実施形態に係る対称性喪失レベルの定量的把握方法のブロック対角化により非対角成分（連成項）が出現した場合の分析例について説明する。

【0071】

図３７（ａ）は、正方形要素（８×８）のＸＹ平面の対称性を説明するための基本構造の剛性行列を模式的に示している。また、図３７（ｂ）は、該ＸＹ平面に存在する対称性のバリエーションを示している。ここで、例えば、（１，１）の対称性は、行列対角成分（１，１）に対応し、（１，２）の対称性は、行列対角成分（３，３）に対応し、（１，３）の対称性は、行列対角成分（２，２）に対応し、（１，４）の対称性は、行列対角成分（４，４）に対応し、（２，１）＋（左側）の対称性は、行列対角成分（５，５）に対応し、（２，１）＋（右側）の対称性は、行列対角成分（７，７）に対応し、（２，１）−（左側）の対称性は、行列対角成分（６，６）に対応し、（２，１）−（右側）の対称性は、行列対角成分（８，８）に対応する。

【0072】

図３７（ｃ）は、正方形要素（８×８）に乱れが生じて連成項が出現した剛性行列を模式的に示している。図３７（ｃ）に示すように、行列成分（２，１）及び（１，２）に連成項が出現している場合、（１，１）及び（１，３）の対称性が重なって、対称性が低下していることが分かる。

【0073】

また、図３８（ａ）は、正方形要素（６×６）に所定の荷重が作用している様子を示している。また、図３８（ｃ）〜（ｈ）はそれぞれ、該正方形要素（６×６）に存在する対称性のバリエーションを示している。ここで、各バリエーションにおけるスケールは同じではない。図３８（ｃ）〜（ｈ）に示す各対称性が完全に分解できるとすると、該各対称性を重ね合わせることによって、図３８（ｂ）に示すような変形が生じる。これは、非対称荷重を作用させたときの非対称変形に相当する。すなわち、荷重の対称性が低くても（変換後の行列に連成項がある場合は要素の対称性が低下する）、各対称性の空間で完全に独立して図３８（ｂ）に示す変形が得られれば、図３８（ａ）に示す荷重、つまり、力のベクトルは、図３９に示す変換則（ベクトル場の剛性方程式）を用いて、容易に求めることができる。ここで、図３９に示すように、例えば、剛性行列の成分（４，１）や（１，４）に連成項が出現している場合において、該連成項＞０のとき、モードｃとモードｆとがカップリングして対称性の低下が生じる。

【0074】

尚、以上に説明した実施形態においては、直積を２つ用いた２次元での対称性喪失レベルの定量的把握について説明したが、同様に、直積を３つ用いて３次元での対称性喪失レベルの定量的把握を行えることは言うまでもない。

【0075】

以下、本実施形態に係る対称性喪失レベルの定量的把握方法を用いた具体的な応用解析例について、図面を参照しながら説明する。

【0076】

（応用解析例１）
図４０は、本実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる測定対象物の一例を示す斜視図である。図４０に示すように、測定対象物１０は、規則的に配列された複数の粒子１１の一部が他の粒子又は空洞１２に置換された構造を持つ。このような場合も、前述のステップＳ１の処理において、測定対象物１０の幾何学的状態の画像データ等に対して、例えば有限要素法（ＦＥＭ）等を適用することにより、測定対象物１０の配列情報を解析データとして行列化（離散モデル化）することが可能である。すなわち、測定対象物１０において、他の粒子又は空洞１２は、粒子１１とは異なる剛性（空洞１２の場合は０剛性）を持つので、前述の周期格子構造における不均一な剛性分布の場合と同様に取り扱うことができる。具体的には、他の粒子又は空洞１２をターゲット（自由度）として追跡（計測）することによって、周期対称性が破れる様相（変形パターン）を認識することができるので、測定対象物１０の対称性を指標化することができると共に、測定対象物１０におけるひずみ（又は応力）の大きさや主ひずみの方向等を把握することができる。

【0077】

（応用解析例２）
応用解析例１と同様に、規則的に配列された複数の貫通孔を有するシート状構造（以下、多孔シートという）に対しても、例えばＦＥＭを用いて配列情報を行列化（離散モデル化）することが可能である。ここで、配列情報を抽出するためのメッシュを再分割することによって、解像度を向上させてもよい。

【0078】

図４１、図４２及び図４３は、多孔シートの例を示す斜視図である。尚、図４１、図４２及び図４３において、多孔シートに荷重を作用させる前の形状を配線で示している。

【0079】

図４１に示す多孔シート２０は、格子状に配列された複数の貫通孔２１を有しており、応力２５を受けて多孔シート２０が延びるように変形している。

【0080】

図４２に示す多孔シート３０は、千鳥状に配列された複数の貫通孔３１を有しており、応力３５を受けて多孔シート３０が延びるように変形している。

【0081】

図４３に示す多孔シート４０は、千鳥状に配列された複数の貫通孔４１を有しており、応力（せん断力）４５を受けて多孔シート４０がねじれるように変形している。

【0082】

以上に述べたような多孔シートの利点として、貫通穴のサイズや配置などの幾何学的な条件を利用することによって、局所応力集中を誘発させることができると共に、応力やひずみの感度を格段に向上させることができる。

【0083】

このような多孔シートの特性を利用して、例えば応力集中の程度に応じて変色するシートであれば、「ひずみ分布」の模様を誘発させる際に、貫通孔の存在によって変色を確実に促進させることが可能となる。このように、対称性喪失レベルの測定において多孔シートを用いて応力集中現象を効果的に生起させることにより、ひずみ分布や応力分布において例えばせん断帯発生によって対称性を低下させると共にひずみや応力の局所化を促すことができる。すなわち、多孔シートを用いることによって、例えばひずみ計測の対象物におけるひずみレンジの感度を格段に向上させることが可能となる。これは、言うまでもなく、様々な適用が可能である。

【0084】

また、前述のような多孔シートを用いた場合、ひずみ分布や応力分布は、対称性が低下する様相を反映したものとなるので、多孔シートに生じる例えば縞模様や放射的な模様（例えば応力が高い箇所）のように、格子構造における変形パターンと類似する多様な対称パターンが出現する可能性がある。従って、多孔シートのように幾何学的に周期性を有する構造体を用いることにより、対象物の形状に関わりなく、対称性の指標化や構造解析を汎用的に実施することができる。

【産業上の利用可能性】

【0085】

本発明は、例えば周期対称構造を持つ構造系における隠れた対称性のレベルを定量的に把握する方法及び装置に関し、機械工学、情報工学、品質管理、材料開発分野などの他分野において好適である。

【符号の説明】

【0086】

１ ＣＰＵ
２ 記憶装置
３ 画像認識装置
４ ネットワーク
１０ 測定対象物
１１ 粒子
１２ 他の粒子又は空洞
２０、３０、４０ 多孔シート
２１、３１、４１ 貫通孔
２５、３５、４５ 応力

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】

【図25】

【図26】

【図27】

【図28】

【図29】

【図30】

【図31】

【図32】

【図33】

【図34】

【図35】

【図36】

【図37】

【図38】

【図39】

【図41】

【図42】

【図43】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図40】