(57)【特許請求の範囲】
【請求項１】
測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する工程ａと、
前記基本構造を表す行列を対角化する行列である変換フィルタを求める工程ｂと、
行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する工程ｃと、
前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する工程ｄとを備えていることを特徴とする対称性喪失レベルの定量的把握方法。
【請求項２】
請求項１に記載の対称性喪失レベルの定量的把握方法において、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とする対称性喪失レベルの定量的把握方法。
【請求項３】
請求項１に記載の対称性喪失レベルの定量的把握方法において、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とする対称性喪失レベルの定量的把握方法。
【請求項４】
測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａと、
前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂと、
行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃと、
前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する手段Ｄとを備えていることを特徴とする対称性喪失レベルの定量的把握装置。
【請求項５】
請求項４に記載の対称性喪失レベルの定量的把握装置において、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とする対称性喪失レベルの定量的把握装置。
【請求項６】
請求項４に記載の対称性喪失レベルの定量的把握装置において、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とする対称性喪失レベルの定量的把握装置。
【請求項７】
コンピュータを、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａ、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂ、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃ、及び、前記ブロック行列の個数と大きさから前記測定対象物の対称性喪失レベルを判別する手段Ｄとして機能させるためのプログラム。
【請求項８】
請求項７に記載のプログラムにおいて、
前記測定対象物は、規則的に配列された複数の粒子の一部が他の粒子又は空洞に置換された構造を持つことを特徴とするプログラム。
【請求項９】
請求項７に記載のプログラムにおいて、
前記測定対象物は、複数の貫通孔が規則的に配列されたシート状構造を持つことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば周期対称性を持つ材料や構造における隠れた対称性のレベルを定量的に把握する方法及び装置に関する。
【背景技術】
【０００２】
古典的な対称性の研究は、静的な現象を前提として、対象とする系に対称性が存在するか否かということに力点がおかれてきた。
【０００３】
物理学では、対称性が低下することは「対称性の自発的破れ」又は単に「対称性の破れ」と呼ばれ、レイリー・ベナール対流やクエット・テイラーの実験などが知られている。
【０００４】
静力学の分野における分岐階層構造の研究においては、系はその対称性の部分群に階層的に支配されていることが明らかにされ、構造解析に群の標準分解定理を取り入れた群論的分岐理論への展開が行われてきた（非特許文献１〜４参照）。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Ikeda,K.他、Bifurcation analysis of symmetric structures using block-diagonalization、Computer Methods in Applied Mechanics and Engineering、86(2)、pp.215-243、1991年
【非特許文献２】Ikeda,K.他、Block-diagonalization analysis of symmetric plates、Int. J. of Solids and Structures、29(22)、pp.2779-2793、1992年
【非特許文献３】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年
【非特許文献４】有尾一郎、群積による直方体ソリッド要素の剛性行列のブロック対角化法、土木学会応用力学論文集、Vol.1、1998年
【発明の概要】
【発明が解決しようとする課題】
【０００６】
対称性の高い格子構造、例えば均一な物質の格子構造（ナノストラクチャー、マイクロストラクチャー、又は規則的に配置されたコロイド粒子等を含む）等において、格子欠陥等の周期性を持つ不完全パターンが存在する場合、そのような状態の認知や理解は単純に見えて一般に極めて複雑であり、例えば、どの程度の格子欠陥状態にあるのかを定量化することは困難である。特に、周期性を持つ不完全性を人間が認識できるのは、格子構造の自由度が非常に小さい場合に限られており、通常解析しようとする構造の自由度では、周期性を持つ不完全性を人間が認識することは不可能である。特に、複数のパターンによって不完全性が繰り返されている場合、人間の認識では対称性のレベルを理解することは困難で認知エラーを起こしやすい。
【０００７】
ところが、元の完全な対称性を基準として、どの程度の不完全性が存在しているのかを定量的に把握できる方法、例えば、破壊等の進行に伴い不完全性を表す乱れパターンが変化して元の高い対称性が段階的に失われているような場合における対称性レベルを定量的に把握できる方法は皆無であった。
【０００８】
前記に鑑み、本発明は、例えば周期的な配列を持つ格子構造等のように高い対称性を持つ構造において、ある格子欠陥等の乱れが生じている場合に、その乱れに起因する対称性喪失のレベルを定量的に把握できるようにすることを目的とする。
【課題を解決するための手段】
【０００９】
前記の目的を達成するために、本発明に係る対称性喪失レベルの定量的把握方法は、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する工程ａと、前記基本構造を表す行列を対角化する行列である変換フィルタを求める工程ｂと、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する工程ｃとを備えている。
【００１０】
また、本発明に係る対称性喪失レベルの対称性喪失レベルの定量的把握装置は、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａと、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂと、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃとを備えている。
【００１１】
さらに、本発明に係るプログラムは、コンピュータを、測定対象物から配列情報を取得して行列化すると共に、前記測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する手段Ａ、前記基本構造を表す行列を対角化する行列である変換フィルタを求める手段Ｂ、及び、行列化された前記測定対象物の前記配列情報に変換フィルタを適用し、ブロック行列に分解する手段Ｃとして機能させる。
【００１２】
すなわち、本発明は、物質内部における対称性の分岐階層構造を、破壊現象や転移現象などの構造不安定性現象を対象とする構造工学的見地から分析するために、均質場からの摂動パターンを想定して、対称性の階層的喪失の追跡を行うものである。
【００１３】
また、本発明は、隠れた対称性の階層構造を確認するために、例えば、群の直積（群のテンソル積）を記述表現した上で、フーリエ級数列からなる選点直交性を利用して、場の（偏）微分支配方程式をその対称性を張る空間に変換し、既往のブロック対角化（ＢＤＭ）の状況（ブロック数、行列サイズ、非対角成分（連成項）の状態など）に基づいて、隠れた対称性を見つけ出し、対称性の階層レベルを定量的に把握するものである。
【００１４】
さらに、本発明は、構造物の幾何学的対称性に基づく従来のＢＤＭとは異なり、新たに均質場や規則的に配置されている系からの摂動パターンを動的現象のある段階として捉え、周期位相構造における「隠れた対称性」の構造及び分岐の系列を見つけ出すことを主眼とするものである。すなわち本発明は、例えば、完全な直交関数列から構成される座標変換フィルターを、境界条件や力学的変形を伴う位相格子系に適用し、変換後の行列表示から、該系における階層的な対称性の喪失（対称性の喪失レベル）を得るものである．
以上に説明した本発明によると、例えば周期構造を持つような高い対称性を持つパターンの認識をコンピュターシステムによって指標化することが可能になると共に、その対称性の変化に伴うパターン崩壊の定量的な把握が各種装置やコンピュターシステムの利用によって可能になる。
【００１５】
また、例えば新しい材料開発では、原子・分子の構造や粒子（ナノ粒子やマイクロ粒子を含む）が規則正しく配置された組成における配列の乱れ方や質的変化を把握することによって、物理化学的に重要で且つ本質的な現象の変化が判明する場合が多いが、本発明は、このような配列の乱れ方や質的変化を体系的に定量化することを可能にする。言い換えると、本発明は、多様な組成を持つ物質における対称性崩壊に起因する多様なパターン変化を魚群（物質）全体の動きの変化に例えると、漁網の種類（対称性の認識方法）によって、その時点での魚群の形（対称性の階層レベル）を定量的に捕獲することを可能にする。従って、本発明は、新しい材料開発のための指標化や予測等を与えることを可能にする。また、本発明の手法を用いることによって、構造解析において不変を付与し、超並列計算をも可能とする。
【００１６】
さらに、本発明は、例えば材料開発分野における物質の周期微視構造の「対称性の破れ」に関する技術、方法、認識等を計算によって具現化するものであって、対称性の体系的な分類と対称性のレベルの定量化を実現するものであるが、本発明は、開発された材料の物理化学変化反応後に適用することも可能である。
【発明の効果】
【００１７】
本発明によると、ブロック行列の個数から、対象物に残った対称性の数や種類が分かり、ブロック行列の大きさから対称性の大きさが分かるので、これらの情報から、例えば、対象物へ印加されている力の大きさや方向が分かる。尚、対象物として、例えば分子構造設計、材料設計、建築物などにも本発明は適用可能である。
【００１８】
また、本発明は、例えば規則正しく配列された格子構造中の欠陥を正しく認識することが困難であったという課題を、幾何学的な情報に基づいて欠陥を定量的に把握することによって解決するものであり、この欠陥の度合いから対称性のレベルを指標化できるので、複雑な対称性の分類を機械的に実施できるという利点がある。
【００１９】
さらに、本発明は、対称性に対する高度で体系的な群論の表現論を応用することによって、物質配列構造の隠れた対称性に基づく質的変化の予測を可能とするものであり、本発明によって、例えば、その時点でのエネルギー準位及び特異点、変形パターンの数及び種類、並びに対称性のレベル及び系統などの情報を得ることができる。すなわち、本発明には、物質配列構造の質的変化についての科学的予見を可能とし、且つ隠れた対称性の存在を体系的に把握できるという利点がある。これによって、対称性が見え隠れする事象を的確に捉えることができると共に、例えば物理変形パターンの種類などを正確に識別することができるので、例えばナノテクノロジー等の分野において、微細構造の崩壊や構造色による可視光の変化等が説明可能となる。
【図面の簡単な説明】
【００２０】
【図１】図１は、実施形態に係る対称性喪失レベルの定量的把握方法のフロー図である。
【図２】図２は、実施形態に係る対称性喪失レベルの定量的把握方法を実施するための装置の概略構成を示す図である。
【図３】図３は、節点列Ｄn の対称性を示す図である。
【図４】図４は、群Ｄn ×Ｄn'の既約表現の次数を示す図である。
【図５】図５は、群積により拡張された既約表現の空間を示す図である。
【図６】図６は、位相格子構造（ｐ×ｑ＝６×６）を示す図である。
【図７】図７（ａ）は、図６に示す構造の剛性行列を示す図であり、図７（ｂ）は、図７（ａ）に示す剛性行列に変換フィルタを適用することによって得られた行列を示す図である。
【図８】図８（ａ）は、図６に示す構造に乱れが生じた構造の一例を示す図であり、図８（ｂ）は、図８（ａ）に示す構造の剛性行列を示す図であり、図８（ｃ）は、図８（ｂ）に示す剛性行列に変換フィルタを適用することによって得られた行列を示す図である。
【図９】図９（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図９（ｂ）は、図９（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１０】図１０（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１０（ｂ）は、図１０（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１１】図１１（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１１（ｂ）は、図１１（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１２】図１２（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１２（ｂ）は、図１２（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１３】図１３（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１３（ｂ）は、図１３（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１４】図１４（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１４（ｂ）は、図１４（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１５】図１５（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１５（ｂ）は、図１５（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１６】図１６（ａ）は、図６に示す構造に乱れが生じた構造の他例を示す図であり、図１６（ｂ）は、図１６（ａ）に示す構造に対して実施形態に係る対称性喪失レベルの定量的把握方法を適用することによって得られた行列を示す図である。
【図１７】図１７は、図７（ｂ）、図８（ｃ）、図９（ｂ）、図１０（ｂ）、図１１（ｂ）、図１２（ｂ）、図１３（ｂ）、図１４（ｂ）、図１５（ｂ）及び図１６（ｂ）に示すブロック対角化行列の既約表現組成を示す図である。
【図１８】図１８は、図７（ｂ）、図８（ｃ）、図９（ｂ）、図１０（ｂ）、図１１（ｂ）、図１２（ｂ）、図１３（ｂ）、図１４（ｂ）、図１５（ｂ）及び図１６（ｂ）に示すブロック対角化行列における階層的対称性喪失の推移を示す図である。
【図１９】図１９（ａ）は、平面ｆｃｃ格子構造の面心格子構造を示す図であり、図１９（ｂ）は、図１９（ａ）に示す構造に対応する点群構成を示す図である。
【図２０】図２０は、図１９（ａ）、（ｂ）に示す構造の剛性行列を示す図である。
【図２１】図２１（ａ）〜（ｄ）は、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列の成分（基底成分）を調和関数により抽出した結果を示す図である。
【図２２】図２２は、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列を示す図である。
【図２３】図２３は、図２２に示す座標変換行列を、図２０に示す剛性行列に適用することによって得られた行列を示す図である。
【図２４】図２４は、図２３に示す行列に所定の行・列操作をすることによって得られた行列を示す図である。
【図２５】図２５は、図２４に示す行列を整理することによって得られたブロック対角化行列を示す図である。
【図２６】図２６は、４質点結合（トラス）の構造を示す図である。
【図２７】図２７は、図２６に示す構造の基本構造を示す図である。
【図２８】図２８は、図２７に示す基本構造に乱れが生じた構造の一例を示す図である。
【図２９】図２９は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。
【図３０】図３０は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。
【図３１】図３１は、図２７に示す基本構造に乱れが生じた構造の他例を示す図である。
【図３２】図３２は、３質点結合（トラス）の構造を示す図である。
【図３３】図３３は、図３２に示す構造の基本構造を示す図である。
【図３４】図３４は、図３３に示す基本構造に乱れが生じた構造の一例を示す図である。
【図３５】図３５は、図３３に示す基本構造に乱れが生じた構造の他例を示す図である。
【図３６】図３６は、図３３に示す基本構造に乱れが生じた構造の他例を示す図である。
【図３７】図３７（ａ）は、正方形要素（８×８）のＸＹ平面の対称性を説明するための基本構造の剛性行列を模式的に示す図であり、図３７（ｂ）は、該ＸＹ平面に存在する対称性のバリエーションを示す図であり、図３７（ｃ）は、正方形要素（８×８）に乱れが生じて連成項が出現した剛性行列を模式的に示す図である。
【図３８】図３８（ａ）は、正方形要素（６×６）に所定の荷重が作用している様子を示す図であり、図３８（ｂ）は、図３８（ａ）に示す正方形要素（６×６）に生じた変形を示す図であり、図３８（ｃ）〜（ｈ）はそれぞれ、正方形要素（６×６）に存在する対称性のバリエーションを示す図である。
【図３９】図３９は、図３８（ａ）に示す荷重を求めるためのベクトル場の剛性方程式を示す図である。
【図４０】図４０は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる測定対象物の一例を示す斜視図である。
【図４１】図４１は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。
【図４２】図４２は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。
【図４３】図４３は、実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる多孔シートを示す斜視図である。
【発明を実施するための形態】
【００２１】
以下、本発明の一実施形態に係る対称性喪失レベルの定量的把握方法及び定量的把握装置について、図面を参照しながら説明する。
【００２２】
図１は、本実施形態に係る対称性喪失レベルの定量的把握方法のフロー図である。
【００２３】
まず、ステップＳ１において、測定対象物から配列情報（例えば格子点の幾何学的位置情報及び格子点同士の結合情報）を取得して行列化すると共に、該測定対象物に乱れが生じていない場合の基本構造を表す行列を抽出する。すなわち、ステップＳ１においては、物質内部における対称性の分岐階層構造を構造工学的見地から分析するために、言い換えると、対称性の階層的喪失を追跡するために、均質場からの摂動パターンを想定する。ここで、測定対象物からの配列情報の取得には、例えばＳＥＭ（走査型電子顕微鏡）写真等の撮影及び例えば格子点を目標ターゲットとしたデジタル画像解析が可能な画像認識装置を用いてもよい。
【００２４】
次に、ステップＳ２において、ステップＳ１で抽出した基本構造を表す行列を対角化する行列である変換フィルタを求める。すなわち、ステップＳ２においては、隠れた対称性の階層構造を確認するために、例えば、群の直積（群のテンソル積）を記述表現した上で、フーリエ級数列からなる選点直交性を利用して、測定対象物の基本構造（元の周期構造）を対角化する変換フィルター（完全な直交関数列から構成される座標変換フィルター）を作成する。
【００２５】
次に、ステップＳ３において、ステップＳ１で行列化された測定対象物の配列情報に、ステップＳ２で求めた変換フィルタを適用し、ブロック行列に分解する。すなわち、ステップＳ３においては、場（測定対象物の配列情報）の支配方程式をその対称性を張る空間に変換し、既往のブロック対角化（ＢＤＭ）の状況（ブロック数、行列サイズ、非対角成分（連成項）の状態など）に基づいて、隠れた対称性を見つけ出し、対称性の階層レベルを定量的に把握する。
【００２６】
図２は、図１に示す本実施形態に係る対称性喪失レベルの定量的把握方法を実施するための装置（コンピュターシステム）の概略構成を示す図である。
【００２７】
図２に示すように、例えば、演算手段としてのＣＰＵ（中央演算装置）１、メモリ等の記憶装置２、及びＳＥＭ写真等を撮影可能な画像認識装置３等がネットワーク４を介して接続されることにより、コンピュターシステムが構成されている。ここで、前述のステップＳ１〜Ｓ３の処理は、図２に示すコンピュターシステムをステップＳ１〜Ｓ３の各処理を実施する手段として機能させるプログラムによって行われる。
【００２８】
以上に説明した本実施形態によると、ステップＳ３で得られたブロック行列の個数から、測定対象物に残った対称性の数が分かると共に、ステップＳ３で得られたブロック行列の大きさから、測定対象物に残った対称性の大きさが分かるので、これらの情報から、例えば、測定対象物へ印加されている力の大きさや方向を知ることができる。
【００２９】
尚、本実施形態において、測定対象物の種類は特に限定されるものではないが、例えば分子構造設計、材料設計、建築物などを測定対象物とすることも可能である。
【００３０】
また、本実施形態によると、例えば規則正しく配列された格子構造中の欠陥を正しく認識することが困難であったという課題を、幾何学的な情報に基づいて欠陥を定量的に把握することによって解決できると共に、この欠陥の度合いから対称性のレベルを指標化できるので、複雑な対称性の分類を機械的に実施することを可能にする。
【００３１】
さらに、本実施形態によると、対称性に対する高度で体系的な群論の表現論を応用することによって、物質配列構造の隠れた対称性に基づく質的変化の予測を行うので、例えば、その時点でのエネルギー準位及び特異点、変形パターンの数及び種類、並びに対称性のレベル及び系統などの情報を得ることができる。すなわち、本実施形態には、物質配列構造の質的変化についての科学的予見を可能とし、且つ隠れた対称性の存在を体系的に把握できるという利点がある。これによって、対称性が見え隠れする事象を的確に捉えることができると共に、例えば物理変形パターンの種類などを正確に識別することができるので、例えばナノテクノロジーや材料設計などの分子・原子配列に関係する分野において、微細構造の崩壊や構造色による可視光の変化等が説明可能となる。
【００３２】
以下、本実施形態で利用している各種物理数学的手法及びそれを用いた各種計算例について、詳細に説明する。
【００３３】
（直交関数列と群対称性）
まず、ノルム完備な空間を満足する直交関数列を定義した上で、位相空間における該関数列の対称性を利用して、該関数列を間接的に変位関数に適合させる方法について、説明する。
【００３４】
変位関数（直交関数列）上の節点列に対する有限な対称性として、例えば図３に示すような、周期軌道上の変位関数（三角関数列）が離散化された節点列の対称性、具体的には、三角関数が振舞う節点変位の対称群（節点列）Ｄn について考える。節点列Ｄn の全長をｎ等分割すると、正規長さｘ／Ｌ（０≦ｘ≦２Ｌ）は、節点並びによって２ｊ／ｎ（ｊ＝１、・・・、ｎ）と置換される。ここで、節点列Ｄn に対応する直交関数列は、
0005912541.tif 000002
・・・（式１）
と表される。また、１次既約表現（１，３）Ｄn に対応する関数列は、
0005912541.tif 000003
・・・（式２）
と表される。このような節点列に対応する有限群の対称性は、変位関数の直交空間を求める際の組成となる。
【００３５】
（位相格子空間の対称性）
群の構造については、群の直積分解を利用することにより、いくつかの群の組成を知ることができる。以下、群の直積分解を用いて、位相平面上にある支配方程式の「隠れた対称性」を把握する方法について説明する。
【００３６】
以下、簡単のため、群積について、二面体群の直積Ｇ＝Ｄn ×Ｄn'に限定して説明する。群Ｇの既約表現全体は、
0005912541.tif 000004
・・・（式３）
と与えられる。ここで、μは群Ｄn の既約表現を、μ’は別の群Ｄn'の既約表現をそれぞれ表す。群Ｄn ×Ｄn'の既約表現（μ，μ’）の次数を図４に示す。このＤn ×Ｄn'不変系の空間Ｘは、
0005912541.tif 000005
・・・（式４）
のように拡張される。ここで、Ｘ（μ）、Ｘ（μ’）は群Ｄn 、Ｄn'の既約表現の空間を表す。２つの群Ｄn 、Ｄn'からなる群元｛ｇ∈Ｄn ｜ｇ１，・・・，ｇｎ｝、｛ｇ’∈Ｄn'｜ｇ１’，・・・，ｇｎ’｝の積は、
0005912541.tif 000006
・・・（式５）
と表される。また、この群元の積からなる表現行列は既約表現μ，μ’の表現行列の直積
0005912541.tif 000007
・・・（式６）
と表される。この表現行列を用いて同値変換
0005912541.tif 000008
・・・（式７）
が可能となり、この座標変換行列Ｈの成分は基底関数φ（μ）、ψ（μ’）から構成される。
【００３７】
以下、この有限群の直積Ｄn ×Ｄn'の組成を前述の「直交関数列と群対称性」で述べた直交関数列を用いて表す。ここで、この拡張された既約表現の空間を図５に示す。２つの群Ｄn 、Ｄn'の基底関数となるφ（μ）、ψ（μ’）がそれぞれ、｛φ（μ）｜φ（１，ｉ）Ｄn，・・・，φ（2，ｉ）Ｄn，・・・｝、｛ψ（μ’）｜ψ（１，j）Ｄn'，・・・，ψ（2，j）Ｄn'，・・・｝と表されるとき、φ（μ）、ψ（μ’）の直積は、
0005912541.tif 000009
・・・（式８）
と表される。この基底関数の直積を各節点に対応させた成分からなる既約表現（μ，μ’）の座標変換行列は、
0005912541.tif 000010
・・・（式９）
となる。
【００３８】
ここで、このＨはユニタリーな行列であり、ＨＴＨ＝Ｉと単位行列となる。２変数関数ｗ（ｘ，ｙ）を未知関数とする偏微分方程式Ｆ（ｗ）と、（式５）の群元の積からなる表現行列Ｔとの間で、
0005912541.tif 000011
・・・（式１０）
の群同変条件が成り立つとき、ＢＤＭが可能となる。このＢＤＭの座標変換によって既約表現毎にブロック分解すると、
0005912541.tif 000012
・・・（式１１）
と変換される。この同値変換によって関数ｗ（ｘ，ｙ）の隠れた対称性を調べることができる。
【００３９】
以下、３次元空間にも適用可能ではあるが、簡単のため、格子の例として、２次元の節点列の対称性を利用して、有限群の直積によって格子上の関数の対称性を表現する。ここで、完全な格子の対称性を記述するために、周期的な境界を持たせた位相空間（左右及び上下の端節点を巡回させた位相構造）を考えた上で、これをすべて固有の対称性を張る空間へ変換後に境界条件を導入するものとする。すなわち、境界条件は２次的要因による制約された対称条件として取り扱う。
【００４０】
周期２πに正規化された位相領域
0005912541.tif 000013
・・・（式１２）
に平面格子を導入する。ここで、Ｌｘ、Ｌｙはｘ，ｙ方向の長さを表す。平面上の格子点（ｐ，ｑ）（ｐ＝１，・・・，ｎ，ｑ＝１，・・・，ｎ’）を設定すると、格子間隔はｈ＝２π／ｎ、ｈ’＝２π／ｎ’となり、ｎ，ｎ’はそれぞれｘ，ｙ方向の分割数である。従って、離散化された格子座標は、
0005912541.tif 000014
・・・（式１３）
と設定できる。格子座標（ｐｈ，ｑｈ’）の任意関数が二重の直交関数列
ｗｐ，ｑ≡ｗ（ｐｈ，ｑｈ’）＝ｕ（２πｐ／ｎ）×ｕ（２πｑ／ｎ’）
と表されるとき、（式１）より、
0005912541.tif 000015
・・・（式１４）
と分解できる。ここで、ａ（μ，μ’）は既約表現（μ，μ’）に対する未定係数である。例えば、（式１４）の左辺が既知の関数であれば、（ｐ，ｑ）に対する支配方程式は、
0005912541.tif 000016
と表される。ここで、ｗ０，ｑ＝ｗｎ，ｑ、ｗｐ，０＝ｗｐ，ｎ’である。このように全ての格子点について方程式を組み立てれば、
0005912541.tif 000017
のような連立方程式を得る。ここで、
0005912541.tif 000018
である。従って、未定係数ａ（μ，μ’）は例えばガウスの消去法などの数値計算によって求めることが可能である。また、行列Ａは方程式Ｆ（ｗ）中の関数ｗの対称性（直交関数列）に依存する。例えば、関数ｗが一定値や二重フーリエ級数列であれば、Ａは対角行列となり、関数ｗが非対称であれば、Ａは密で連成された行列となる。
【００４１】
（「隠れた対称性」の数値計算）
前述のように、均質場に「隠れた対称性」が存在する場合があると共に、均質領域が分岐に伴い周期的なパターンを形成しながら乱れていく場合がある。以下、均質な位相格子モデルにおける剛性分布パターンと位相対称性との関係において、剛性方程式のブロック対角化行列Ｋ’の直交性によって「隠れた対称性」を数値的に把握できることについて説明する。尚、以下の説明においては、均質場の離散化を伴う要素解析は場の位相構造に依存し、該位相構造から周期パターン変形が形成されるものとして、その階層的な対称性を探るものとする。具体的には、均等に分割された格子上の８隣接の部材剛性をモデル化し、均質な剛性分布を持つ構造の面外変位を自由度（１節点１自由度）として有限要素解析を行う。
【００４２】
まず、前述のステップＳ１の処理により、図６に示す閉じた位相格子空間（ｐ×ｑ＝６×６）の構造から、図７（ａ）に示すような帯行列化されたＫ（剛性行列）を形成する。ここで、図７（ａ）において、＋は正値を、−は負値を、・はゼロ成分をそれぞれ表す（以下の他の図面において同様）。
【００４３】
次に、前述のステップＳ２の処理により求めた変換フィルタを適用し、前述のステップＳ３の処理によって、図７（ａ）に示す剛性行列Ｋを座標変換して、図７（ｂ）に示すＫ’（パターンＰ０）を得る。ここで、Ｋ’は完全に対角行列となり、最初の対角成分以外はそれぞれ直交し異なる正値となった。すなわち、図６に示す構造の剛性行列が線形独立な既約表現と一致し、この離散構造では最も対称位数が高い構造となった。図７（ｂ）から明らかなように、この対称性は、１次元のＤ３６不変配下の群と等価な群を持つ。尚、図７（ｂ）に示すＫ’の最初の対角成分０は剛体変形モードに対応する。
【００４４】
以上に説明したように、ＢＤＭによるブロックの数が対称性の高低を表すことを利用し、対称性の階層構造を調べることができる。また、以上のような位相格子構造となる大規模な離散系では対称性の利用価値が高くなる。従って、均質場の分岐階層構造に、パターン化された位相格子構造の群論的な変換則を用いることによって、該構造が持つ隠れた対称性及び階層的な対称性の仕組みを明らかにすることが可能となる。
【００４５】
例えば、図８（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の一例（チェッカーボードパターン）を示している。図８（ａ）に示す剛性分布からは、前述のステップＳ１の処理により、図８（ｂ）に示すような剛性行列が得られる。また、前述のステップＳ２の処理により求めた変換フィルタを適用し、前述のステップＳ３の処理によって、図８（ｂ）に示す剛性行列を座標変換すると、図８（ｃ）に示すブロック対角化行列（パターンＰ）が得られる。
【００４６】
また、図９（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図９（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図９（ｂ）に示すブロック対角化行列（パターンＨ）が得られる。
【００４７】
また、図１０（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１０（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１０（ｂ）に示すブロック対角化行列（パターンＡ）が得られる。
【００４８】
また、図１１（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１１（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１１（ｂ）に示すブロック対角化行列（パターンＢ１）が得られる。
【００４９】
また、図１２（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１２（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１２（ｂ）に示すブロック対角化行列（パターンＢ２）が得られる。
【００５０】
また、図１３（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１３（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１３（ｂ）に示すブロック対角化行列（パターンＢ’２）が得られる。
【００５１】
また、図１４（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１４（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１４（ｂ）に示すブロック対角化行列（パターンＣ）が得られる。
【００５２】
また、図１５（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１５（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１５（ｂ）に示すブロック対角化行列（パターンＤ）が得られる。
【００５３】
また、図１６（ａ）は、図６に示す基本構造に乱れ（図中■）が生じた構造の他例を示している。図１６（ａ）に示す剛性分布からは、前述のステップＳ１〜Ｓ３の処理により、図１６（ｂ）に示すブロック対角化行列（パターンＥ）が得られる。
【００５４】
図１７は、前述の各ブロック対角化行列（パターンＰ０、Ｐ、Ｈ、Ａ、Ｂ１、Ｂ２、Ｂ’２、Ｃ、Ｄ、Ｅ）の既約表現組成を示しており、図１８は、前述の各ブロック対角化行列（パターンＰ０、Ｐ、Ｈ、Ａ、Ｂ１、Ｂ２、Ｂ’２、Ｃ、Ｄ、Ｅ）における階層的対称性喪失の推移を示している。
【００５５】
（平面ｆｃｃ格子構造のブロック対角化）
図１９（ａ）は、平面ｆｃｃ格子構造の面心格子構造を示し、図１９（ｂ）は、図１９（ａ）に示す構造に対応する点群構成を示している。図１９（ａ）、（ｂ）において、数字は格子番号を示している。尚、以下の説明では、１節点１自由度のｎ，ｍ＝４の周期構造を前提とする。
【００５６】
図１９（ａ）、（ｂ）に示す構造からは、前述のステップＳ１の処理により、図２０に示すような剛性行列（帯行列）が得られる。ここで、格子点の構成種類（○、●）から、該構造の幾何学的な関係を示す剛性行列は、
0005912541.tif 000019
と表される。
【００５７】
一方、群の直積から得られる、白節点格子点に対応する座標変換Ｈｅは、中点に位置する黒節点とも点群では同型であるので、図１９（ａ）、（ｂ）に示す構造構造に対応する座標変換行列は、
0005912541.tif 000020
と表される。
【００５８】
図２１（ａ）〜（ｄ）は、前記座標変換行列の成分（基底成分）を調和関数により抽出した結果を示す。この基底成分を用いて、前述のステップＳ２の処理により、図２２に示す座標変換行列（変換フィルタ）が得られる。そして、前述のステップＳ３の処理により、つまり、図２２に示す座標変換行列を、図２０に示す剛性行列に適用することによって、図２３に示す行列が得られる。図２３に示す行列に所定の行・列操作をすると、図２４に示す行列が得られ、該行列をさらに整理すると、図２５に示すブロック対角化行列が得られる。
【００５９】
（４質点結合構造のブロック対角化）
図２６は、４質点結合（トラス）の構造を示している。図２６において、数字は質点番号を、ｋは質点同士の結合状態（剛性係数）を、ｗは変位を示している。尚、以下の説明では、１質点１自由度の周期境界条件を前提とする。
【００６０】
図２７は、図２６に示す構造の基本構造を示す図である。図２７に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000021
に示す剛性行列Ｋ１が得られる。ここで、ｋは剛性係数を表す（以下の他の行列において同じ）。また、前述のステップＳ２の処理により、
0005912541.tif 000022
に示す座標変換行列（変換フィルタ）Ｈが得られる。そして、前述のステップＳ３の処理により、つまり、変換フィルタＨを剛性行列Ｋ１に適用することによって、
0005912541.tif 000023
となり、完全に対角化される。
【００６１】
また、図２８は、図２７に示す基本構造に乱れが生じた構造の一例を示す。図２８に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000024
に示す剛性行列Ｋ２が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ２に適用することによって、
0005912541.tif 000025
となり、完全に対角化される。
【００６２】
また、図２９は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図２９に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000026
に示す剛性行列Ｋ３が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ３に適用することによって、
0005912541.tif 000027
となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。
【００６３】
また、図３０は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図３０に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000028
に示す剛性行列Ｋ４が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ４に適用することによって、
0005912541.tif 000029
となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。
【００６４】
さらに、図３１は、図２７に示す基本構造に乱れが生じた構造の他例を示す。図３１に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000030
に示す剛性行列Ｋ５が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ５に適用することによって、
0005912541.tif 000031
となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。
【００６５】
（３質点結合構造のブロック対角化）
図３２は、３質点結合（トラス）の構造を示している。図３２において、ａ、ｂ、ｃは質点番号を、ｋは質点同士の結合状態（剛性係数）をしている。尚、以下の説明では、１質点１自由度の周期境界条件を前提とする。
【００６６】
図３３は、図３２に示す構造の基本構造を示す図である。尚、図３３において、数字は質点番号を、ｋは質点同士の結合状態（剛性係数）を、ｗは変位を示している。図３３に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000032
に示す剛性行列Ｋ１が得られる。また、前述のステップＳ２の処理により、
0005912541.tif 000033
に示す座標変換行列（変換フィルタ）Ｈが得られる。そして、前述のステップＳ３の処理により、つまり、変換フィルタＨを剛性行列Ｋ１に適用することによって、
0005912541.tif 000034
となり、完全に対角化される。
【００６７】
また、図３４は、図３３に示す基本構造に乱れが生じた構造の一例を示す。図３４に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000035
に示す剛性行列Ｋ２が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ２に適用することによって、
0005912541.tif 000036
となり、完全に対角化される。
【００６８】
また、図３５は、図３３に示す基本構造に乱れが生じた構造の一例を示す。図３５に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000037
に示す剛性行列Ｋ３が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ３に適用することによって、
0005912541.tif 000038
となり、完全に対角化される。
【００６９】
また、図３６は、図３３に示す基本構造に乱れが生じた構造の他例を示す。図３６に示す構造からは、前述のステップＳ１の処理により、
0005912541.tif 000039
に示す剛性行列Ｋ４が得られる。そして、前述のステップＳ３の処理により、つまり、前述のステップＳ２の処理により得られた変換フィルタＨを剛性行列Ｋ４に適用することによって、
0005912541.tif 000040
となる。すなわち、非対角成分（連成項）が出現しており、基本構造からの対称性の低下が見られる。
【００７０】
（非対角成分の分析）
以下、本実施形態に係る対称性喪失レベルの定量的把握方法のブロック対角化により非対角成分（連成項）が出現した場合の分析例について説明する。
【００７１】
図３７（ａ）は、正方形要素（８×８）のＸＹ平面の対称性を説明するための基本構造の剛性行列を模式的に示している。また、図３７（ｂ）は、該ＸＹ平面に存在する対称性のバリエーションを示している。ここで、例えば、（１，１）の対称性は、行列対角成分（１，１）に対応し、（１，２）の対称性は、行列対角成分（３，３）に対応し、（１，３）の対称性は、行列対角成分（２，２）に対応し、（１，４）の対称性は、行列対角成分（４，４）に対応し、（２，１）＋（左側）の対称性は、行列対角成分（５，５）に対応し、（２，１）＋（右側）の対称性は、行列対角成分（７，７）に対応し、（２，１）−（左側）の対称性は、行列対角成分（６，６）に対応し、（２，１）−（右側）の対称性は、行列対角成分（８，８）に対応する。
【００７２】
図３７（ｃ）は、正方形要素（８×８）に乱れが生じて連成項が出現した剛性行列を模式的に示している。図３７（ｃ）に示すように、行列成分（２，１）及び（１，２）に連成項が出現している場合、（１，１）及び（１，３）の対称性が重なって、対称性が低下していることが分かる。
【００７３】
また、図３８（ａ）は、正方形要素（６×６）に所定の荷重が作用している様子を示している。また、図３８（ｃ）〜（ｈ）はそれぞれ、該正方形要素（６×６）に存在する対称性のバリエーションを示している。ここで、各バリエーションにおけるスケールは同じではない。図３８（ｃ）〜（ｈ）に示す各対称性が完全に分解できるとすると、該各対称性を重ね合わせることによって、図３８（ｂ）に示すような変形が生じる。これは、非対称荷重を作用させたときの非対称変形に相当する。すなわち、荷重の対称性が低くても（変換後の行列に連成項がある場合は要素の対称性が低下する）、各対称性の空間で完全に独立して図３８（ｂ）に示す変形が得られれば、図３８（ａ）に示す荷重、つまり、力のベクトルは、図３９に示す変換則（ベクトル場の剛性方程式）を用いて、容易に求めることができる。ここで、図３９に示すように、例えば、剛性行列の成分（４，１）や（１，４）に連成項が出現している場合において、該連成項＞０のとき、モードｃとモードｆとがカップリングして対称性の低下が生じる。
【００７４】
尚、以上に説明した実施形態においては、直積を２つ用いた２次元での対称性喪失レベルの定量的把握について説明したが、同様に、直積を３つ用いて３次元での対称性喪失レベルの定量的把握を行えることは言うまでもない。
【００７５】
以下、本実施形態に係る対称性喪失レベルの定量的把握方法を用いた具体的な応用解析例について、図面を参照しながら説明する。
【００７６】
（応用解析例１）
図４０は、本実施形態に係る対称性喪失レベルの定量的把握方法の適用対象となる測定対象物の一例を示す斜視図である。図４０に示すように、測定対象物１０は、規則的に配列された複数の粒子１１の一部が他の粒子又は空洞１２に置換された構造を持つ。このような場合も、前述のステップＳ１の処理において、測定対象物１０の幾何学的状態の画像データ等に対して、例えば有限要素法（ＦＥＭ）等を適用することにより、測定対象物１０の配列情報を解析データとして行列化（離散モデル化）することが可能である。すなわち、測定対象物１０において、他の粒子又は空洞１２は、粒子１１とは異なる剛性（空洞１２の場合は０剛性）を持つので、前述の周期格子構造における不均一な剛性分布の場合と同様に取り扱うことができる。具体的には、他の粒子又は空洞１２をターゲット（自由度）として追跡（計測）することによって、周期対称性が破れる様相（変形パターン）を認識することができるので、測定対象物１０の対称性を指標化することができると共に、測定対象物１０におけるひずみ（又は応力）の大きさや主ひずみの方向等を把握することができる。
【００７７】
（応用解析例２）
応用解析例１と同様に、規則的に配列された複数の貫通孔を有するシート状構造（以下、多孔シートという）に対しても、例えばＦＥＭを用いて配列情報を行列化（離散モデル化）することが可能である。ここで、配列情報を抽出するためのメッシュを再分割することによって、解像度を向上させてもよい。
【００７８】
図４１、図４２及び図４３は、多孔シートの例を示す斜視図である。尚、図４１、図４２及び図４３において、多孔シートに荷重を作用させる前の形状を配線で示している。
【００７９】
図４１に示す多孔シート２０は、格子状に配列された複数の貫通孔２１を有しており、応力２５を受けて多孔シート２０が延びるように変形している。
【００８０】
図４２に示す多孔シート３０は、千鳥状に配列された複数の貫通孔３１を有しており、応力３５を受けて多孔シート３０が延びるように変形している。
【００８１】
図４３に示す多孔シート４０は、千鳥状に配列された複数の貫通孔４１を有しており、応力（せん断力）４５を受けて多孔シート４０がねじれるように変形している。
【００８２】
以上に述べたような多孔シートの利点として、貫通穴のサイズや配置などの幾何学的な条件を利用することによって、局所応力集中を誘発させることができると共に、応力やひずみの感度を格段に向上させることができる。
【００８３】
このような多孔シートの特性を利用して、例えば応力集中の程度に応じて変色するシートであれば、「ひずみ分布」の模様を誘発させる際に、貫通孔の存在によって変色を確実に促進させることが可能となる。このように、対称性喪失レベルの測定において多孔シートを用いて応力集中現象を効果的に生起させることにより、ひずみ分布や応力分布において例えばせん断帯発生によって対称性を低下させると共にひずみや応力の局所化を促すことができる。すなわち、多孔シートを用いることによって、例えばひずみ計測の対象物におけるひずみレンジの感度を格段に向上させることが可能となる。これは、言うまでもなく、様々な適用が可能である。
【００８４】
また、前述のような多孔シートを用いた場合、ひずみ分布や応力分布は、対称性が低下する様相を反映したものとなるので、多孔シートに生じる例えば縞模様や放射的な模様（例えば応力が高い箇所）のように、格子構造における変形パターンと類似する多様な対称パターンが出現する可能性がある。従って、多孔シートのように幾何学的に周期性を有する構造体を用いることにより、対象物の形状に関わりなく、対称性の指標化や構造解析を汎用的に実施することができる。
【産業上の利用可能性】
【００８５】
本発明は、例えば周期対称構造を持つ構造系における隠れた対称性のレベルを定量的に把握する方法及び装置に関し、機械工学、情報工学、品質管理、材料開発分野などの他分野において好適である。
【符号の説明】
【００８６】
１                ＣＰＵ
２                記憶装置
３                画像認識装置
４                ネットワーク
１０                測定対象物
１１                粒子
１２                他の粒子又は空洞
２０、３０、４０    多孔シート
２１、３１、４１    貫通孔
２５、３５、４５    応力
【図１】
0005912541.tif 000041
【図２】
0005912541.tif 000042
【図３】
0005912541.tif 000043
【図４】
0005912541.tif 000044
【図５】
0005912541.tif 000045
【図６】
0005912541.tif 000046
【図７】
0005912541.tif 000047
【図１７】
0005912541.tif 000048
【図１８】
0005912541.tif 000049
【図１９】
0005912541.tif 000050
【図２０】
0005912541.tif 000051
【図２１】
0005912541.tif 000052
【図２２】
0005912541.tif 000053
【図２３】
0005912541.tif 000054
【図２４】
0005912541.tif 000055
【図２５】
0005912541.tif 000056
【図２６】
0005912541.tif 000057
【図２７】
0005912541.tif 000058
【図２８】
0005912541.tif 000059
【図２９】
0005912541.tif 000060
【図３０】
0005912541.tif 000061
【図３１】
0005912541.tif 000062
【図３２】
0005912541.tif 000063
【図３３】
0005912541.tif 000064
【図３４】
0005912541.tif 000065
【図３５】
0005912541.tif 000066
【図３６】
0005912541.tif 000067
【図３７】
0005912541.tif 000068
【図３８】
0005912541.tif 000069
【図３９】
0005912541.tif 000070
【図４１】
0005912541.tif 000071
【図４２】
0005912541.tif 000072
【図４３】
0005912541.tif 000073
【図８】
0005912541.tif 000074
【図９】
0005912541.tif 000075
【図１０】
0005912541.tif 000076
【図１１】
0005912541.tif 000077
【図１２】
0005912541.tif 000078
【図１３】
0005912541.tif 000079
【図１４】
0005912541.tif 000080
【図１５】
0005912541.tif 000081
【図１６】
0005912541.tif 000082
【図４０】
0005912541.tif 000083