機械要素の応力解析はじめに各種機械や構造物を設計するとき、その構成部材が通常の使用内で破壊しないように材質や形状寸法が決定される。負荷による変形の大きさや破壊は部材内に生じる応力（単位面積当たりに作用する内力）によって評価される。つまり、垂直応力あるいはせん断応力が材料固有のある限界値に達すると破壊や永久変形が生じると考える。したがって、応力解析は設計上、非常に重要な問題である。一般に機械に用いられる部品（機械要素）の形状、作用する負荷、支持条件などは、おおよそ材料力学で扱ったものより複雑である。このような材料力学により応力や変形を評価することが困難な問題の解析に適用される手法として、実験的に解析する方法と計算機（コンピュータ）を用いた数値解析がある。実験的に解析する方法としてひずみゲージ法や光弾性法が代表的に知られている。ひずみゲージ法は、実部材表面上で測定したい箇所にひずみゲージを接着し、直接的に部材表面上のひずみを測定し、応力に換算する。２次元面の応力、ひずみ分布を測定するには大きな労力が必要である。一方、光弾性法は部材の表面の応力分布が得られ、複雑形状によって起こる応力集中の領域とその応力値を測定することができるのが特徴である。３次元光弾性法によれば、表面のみならず内部の応力分布も得られる。しかし、光弾性の性質を利用するために光弾性体である透明樹脂製の模型を用いるモデル実験であり、実部材の形状、負荷、支持条件などの再現に課題がある。計算機を用いた数値解析の中で、最も広く用いられるのが有限要素法（Finite Element Method、FEM）である。近年の計算機の進歩によって、機械工学の分野においては、設計、製造に CAD、CAM が普及し、より複雑形状な機械要素において最適な設計条件や性能予測が必要になる中、FEM は CAE（Computer Aided Engineering：計算機支援工学）として欠かすことができない重要なツールとなっている。１．目的実験手法（光弾性実験）と数値解析手法（FEM）により比較的簡単な形状の部材の二次元弾性応力解析を行い、機械要素や構造部材に生ずる応力や変形の評価方法を学ぶ。また、光弾性実験法および有限要素法の基礎、特徴を理解する。２．有限要素法の基本概念現実の物体は、すべての点で力のつり合いと変形の連続を満たしているが、有限要素法では解析対象領域を適当な形状と寸法の有限要素（本解析では四辺形）の集合により物体を近似する。図１の例で示すように節点で複数の要素がつながっている。すなわち、それら要素は節点だけで互いに結合され力を伝えると仮定する。解析対象を要素に分割したものを有限要素メッシュ（finite element mesh）とよぶ。要素内の変形様式（変位分布）を仮定すれば、 1 図１引張を受ける円孔を有する帯板の幾何形状と荷重条件、有限要素メッシュ節点に作用する力（節点力）を節点変位の関数として表すことができる。この変位法では、節点力が未知の節点変位の連立一次方程式として表され、その係数は要素の材料特性、座標値等により定められる。一般に物体内部のつり合っている節点では節点力の合計は０であり、外力が与えられた節点での値は既知であるので、連立方程式を解くことにより、節点変位が求められる。一般に、近似解の精度は要素（および節点）の数が増えるほど改善されるが、それにともない計算時間が長くなる。３．境界条件の導入多くの機械要素は図２に示すように対称形状を持ち、対称な負荷を受けることが多い。この場合、対称性より、図の斜線部分のみを解析すればよい。Ａ C O B 図２形状、負荷の対称性図２の場合の境界条件は次のように与えられる。変位境界条件：OA に沿って u（x 方向の変位）＝0、OB に沿って v（y 方向の変位）＝0 負荷境界条件：円孔上および BC に沿って外力＝0、AC 上で与えられる外力の値 2 ４．FEM 解析４．１ソフトウェア現在、多くの汎用有限要素法ソフトウェアが開発され市販されている。それぞれ解析対象、機能やグラフィックユーザーインターフェイスに特徴がある。本解析では、構造解析、応力解析において代表的な FEM ソフトである ABAQUS(Student Edition)を用いる。４．２解析問題解析対象として、図３に示すような比較的簡単な形状で負荷も単純な円板の直径方向圧縮と円孔を有する帯板の引張りの二次元弾性解析を行う。それぞれの問題に対して、その対称性から解析領域および変位と負荷の境界条件を考慮しなければならない。円柱の直径方向圧縮有孔帯板の引張り図３解析問題４．３解析手順ここでは，FEM 解析の大まかな手順を示す．解析に用いる材料特性、解析条件、実際に使用する ABAQUS の操作手順の詳細は別に示す． ① 解析対象のモデル（形状）の作成 ② 材料特性値の定義 ③ 要素特性の定義と選定 ④ モデルの要素分割 ⑤ 各要素の系全体へのアセンブリ（集合化） ⑥ 変位境界条件の設定 ⑦ 負荷境界条件の設定 ⑧ 解析実行 ⑨ 解析結果の表示 3 ５．解析結果と考察通常、解析結果はコンピュータのディスプレイ上に応力分布やひずみ分布の等高線（コンターcontour 図）などで可視化できる。図４に垂直応力σy のコンター図の例を示す。応力の分布が詳細にわかる。図４垂直応力σy のコンター図（１）解析結果の応力、ひずみの分布をコンター図として出力する。（２）応力分布について考察しなさい。・分布の状況、平均の応力値と比較してどのような状態なのかなど（３）ひずみ分布について考察しなさい。・分布の状況、応力分布と合わせて、破壊の危険部位、どのような破壊が予想されるかなど６．研究課題（１）応力分布、ひずみ分布からどのような知見が得られ、それらをどこでどのように活用するか述べなさい。（２）有限要素法の特徴を述べなさい。７．参考文献例えば、山田貴博（訳）、有限要素法 ABAQUS(Student Edition)付 4 丸善有限要素法プログラム ABAQUS による２次元弾性応力解析はじめにここでは、ABAQUS／CAE(Student Edition)を使って２次元弾性応力解析の手順を説明する。解析は、（１）円板の直径方向圧縮、（２）円孔を有する帯板の引張の二つを取り上げる。一般に汎用有限要素法ソフトは、使用する単位系は指定されていない。これにより相対的な解析が可能である。解析者が正しく統一された単位系を想定すればよい。本解析では、以下の単位を想定して解析を進める。長さ： mm、力： N 、応力：Ｎ／mm2＝MPa 準備 ABAQUS／CAE(Student Edition)が起動したら、「セッションの開始」ダイアログボックスの「モデルデータベースの作成」を選択して開始状態にする。１．円板の直径方向圧縮図１に円板の直径方向圧縮の幾何および荷重条件を示す。モデルや荷重条件の対称性を利用すると、全体の 1/4 領域だけをモデル化すればよいことがわかる。円板の直径ｄ＝50 mm 円板の厚さｔ＝ 5 mm とし、材質はエポキシ樹脂として、その材料特性値は縦弾性係数（ヤング率）E＝4000 MPa ポアソン比 ν＝0.38 とする。また、圧縮荷重 P＝1000 N とする。図１円板の直径方向圧縮の幾何 1.1 モデルの作成モデルデータベースにモデルを作成する。モデルツリーの「Model-1」と表示されている「モデル」の上でマウスを右クリックし、「名前の変更」を選択。表示されたダイアログボックスに compression（例として）と入力し、「OK」をクリック。以下の操作は compression モデルに対して行われる。 1.2 パートの作成２次元平面のジオメトリ（幾何）を作成する。（１）モデルツリーの「パート」をダブルクリック。（２）ダイアログボックスで「名前」を circle plate と入力、「モデリング空間」を 2 次元平面、「タイプ」を変形体、「ベースフィーチャ」をシェル、「近似サイズ」に 200 を入力し、「続ける」クリック。（３）「スケッチャー」ツールボックスの右上にある、「円の作成」ツールを選択し、ビューポート上の円の中心となる点（０，０）をクリックし、円の半径位置（２５，０）をクリック。 5 円が書けたら右クリック「作業のキャンセル」で終了。プロンプトエリアの「完了」をクリックして「スケッチャー」を終了。右上の 1/4 領域を残してそれ以外の領域を削除する。（４）「パートモジュール」のメニューから「形状」→「切り取り」→「押し出し」を選択。「スケッチャー」ツールボックスの右上にある、「直線の作成」ツールを選択し、一続きの閉じた直線で削除領域を囲み、右クリック「作業のキャンセル」、プロンプトエリアの「完了」をクリックして終了すると、解析領域の「パート」ができる。 1.3 材料特性の定義（１）モデルツリーの「材料特性」をダブルクリック。（２）ダイアログボックスで「名前」を epoxy と入力、メニューから「機械的」→「弾性」→ 「弾性」を選択、「ヤング率」に 4000、「ポアソン比」に 0.38 を入力。（３）「OK」をクリックで終了。 1.4 要素特性の定義と割り当て 1.3 で定義した材料特性を適用する要素特性を定義し，これを「パート：circle plate」に割り当てる．定義：（１）モデルツリーの「要素特性」をダブルクリック。（２）ダイアログボックスで「名前」を plateSection と入力し，「カテゴリ」から「ソリッド」，「タイプ」の一覧から「均質」，「続ける」をクリック．（３）「要素特性の編集」ダイアログボックスで「材料特性」として「epoxy」，「平面応力／ひずみ要素厚の値」で５を入力，「OK」をクリック．割り当て：（１）モデルツリーの「circle plate」と名づけた「パート」を展開（＋印をクリック）．（２）「要素特性の割り当て」をダブルクリックする。ビューポート上で要素特性を割り当てるパート領域をクリックすると，選択領域がハイライト表示される．（３）プロンプトエリアの「完了」をクリックして終了．「要素特性割り当ての編集」ダイアログボックスですでに定義した要素特性「plateSection」を選択し，「OK」をクリック． 1.5 モデルのメッシュ分割 1.4 で要素特性を定義し，それを割り当てた「circle plate」を選択する要素で分割する．（１）「circle plate」の下にある「メッシュ」をダブルクリック．（２）「メッシュモジュール」のメニューから「メッシュ」→「コントロール」を選択。（３）「メッシュコントロール」ダイアログボックスで「要素形状」から「四辺形」，「テクニック」から「構造」を選択．「OK」をクリックして終了．要素タイプの選択：（１）メニューから「メッシュ」→「要素タイプ」を選択．（２）「要素タイプ」ダイアログボックスで「要素ライブラリ」から「Standard」，「ジオメトリ次数」から「2 次」，「ファミリ」から「平面応力」を選択．「OK」をクリック 6 して終了．「パート：circle plate」に CPS8R 要素が選択される．要素分割：（１）メニューから「シード」→「パート」を選択，「全体シード」ダイアログボックスの「近似全体サイズ」に 2.0 を入力，「OK」をクリック．プロンプトエリアの「完了」をクリックして終了．ここで，節点が必ず配置される場所を正方形で，おおよそ節点が配置される場所を丸印で表される．（２）メニューから「メッシュ」→「パート」を選択，プロンプトエリアの「はい」をクリックするとメッシュ分割される． 1.6 モデルのアセンブリ要素分割された「パート：circle plate」をアセンブリする．（１）モデルツリーの「アセンブリ」を展開し，「インスタンス」をダブルクリック．（２）「インスタンスの作成」ダイアログボックスでパート「circle plate」を選択し，「OK」をクリック． 1.7 解析条件の設定弾性解析であるので静的解析ステップを設定する．（１）モデルツリーの「ステップ」をダブルクリック，「ステップの作成」ダイアログボックスで「プロシージャタイプ：一般」から「Static,General」を選択，「続ける」をクリック．（２）「ステップの編集」ダイアログボックスで「基本」タブの説明に Circle plate under compression 1.8 と入力．「OK」をクリック．境界条件の設定変位と荷重の境界条件を設定するために，設定箇所の定義し，必要な条件を設定する．集合の定義：（１）モデルツリーの「アセンブリ」の下の「集合」をダブルクリック，「集合の作成」ダアログボックスで「名前」を left と入力，「タイプ」はジオメトリを選択し，「続ける」をクリック．ビューポート上でモデルの左側エッジ（稜線）をクリックして選択し，プロンプトエリアの「完了」をクリック．（２）同様の手順でモデルの底辺を bottom として「集合」を定義する．（３）同様の手順でモデルの頂上の点を top として「集合」を定義する．変位境界条件の設定：（１）モデルツリーの「境界条件」をダブルクリック，「境界条件の作成」ダイアログボックスで「名前」を left と入力，「ステップ」として「Step-1」を選択，「カテゴリ」から「機械的」，「選択されたステップに対するタイプ」は「対称/反対称/完全固定」を選択，「続ける」．（２）プロンプトエリアの「集合」をクリック，「領域の選択」ダイアログボックスで left を選択し，「ビューポート内で選択をハイライトする」のチェックを有効にすれば選択されたエッジがハイライト表示される．（3）「続ける」をクリック，「境界条件の編集」ダイアログボックスで「XSYMM」を選択し，「OK」をクリックして終了すると，設定されたエッジに沿って記号が表示される． 7 （4）同様に bottom としたエッジに「YSYMM」の条件を設定する．荷重境界条件の設定：（１）モデルツリーの「荷重」をダブルクリック，「荷重の作成」ダイアログボックスで「名前」を compression と入力，「ステップ」として「Step-1」を選択，「カテゴリ」から「機械的」，「選択されたステップに対するタイプ」から「集中力」を選択，「続ける」．（２）プロンプトエリアの「集合」をクリック，「領域の選択」ダイアログボックスで top を選択すると，荷重作用点がハイライト表示される．（３）「続ける」をクリック，「荷重の編集」ダイアログボックスで CF1 に 0，CF2 に-1000 を入力し，「OK」をクリックして終了すると，指定された点にベクトルが表示される． 1.9 解析ジョブ（解析の単位）の作成と投入（実行）（１）モデルツリーの「ジョブ」をダブルクリック，「ジョブの作成」ダイアログボックスで「名前」を circle-plate と入力，一覧から「compression」を選択し，「続ける」．（２）「ジョブの編集」ダイアログボックスではそのまま「OK」をクリック．（３）モデルツリーの「ジョブ」を展開し，circle-plate と名づけた「ジョブ」の上で右クリックして「ジョブの投入」を選択して，計算を実行させる． 1.10 解析結果の表示（１）ジョブが問題なく完了したら，ステータス表示は「完了」となり，解析結果を結果表示モジュールを使って表示させることができる．（２）ジョブ名を右クリックして「結果」を選択し，応力分布のコンター図を作成する．（３）メニューから「結果」→「フィールド出力」，表示させる出力変数を選択して，所要の結果を得る． 8 ２．円孔を有する帯板の引張図２に幾何および荷重条件を示す。モデルや荷重条件の対称性を利用すると、全体の 1/4 領域だけをモデル化すればよいことがわかる。円孔の直径ｄ＝ 10 mm 帯板の厚さｔ＝ 5 mm b＝ 50 mm 帯板の幅帯板の長さｌ＝60 mm、材質はエポキシ樹脂として、その材料特性値は縦弾性係数（ヤング率）E＝4000 MPa ポアソン比 ν＝0.38 とする。負荷条件として引張応力σy＝2 MPa（P = 500N）とする。図２円孔帯板の引張の幾何 2.1 モデルの作成モデルデータベースにモデルを作成する。モデルツリーの「Model-1」と表示されている「モデル」の上でマウスを右クリックし、「名前の変更」を選択。表示されたダイアログボックスに tension （例として）と入力し、「OK」をクリック。以下の操作は tension モデルに対して行われる。 2.2 パートの作成２次元平面のジオメトリ（幾何）を作成する。（１）モデルツリーの「パート」をダブルクリック。（２）ダイアログボックスで「名前」を plate と入力、「モデリング空間」を 2 次元平面、「タイプ」を変形体、「ベースフィーチャ」をシェル、「近似サイズ」に 300 を入力し、「続ける」クリック。（３）「スケッチャー」ツールボックスの右上にある、「直線の作成：矩形(4 ライン)」ツールを選択し、ビューポート上に矩形の最初のコーナーとなる点（-25，-30）をクリックし、反対のコーナーとなる点（25，30）をクリック。（４）「スケッチャー」ツールボックスの右上にある、「円の作成」ツールを選択し、ビューポート上の円の中心となる点（０，０）をクリックし、円の半径位置（5，０）をクリック。円が書けたら右クリック「作業のキャンセル」で終了。プロンプトエリアの「完了」をクリックして「スケッチャー」を終了。右上の 1/4 領域を残してそれ以外の領域を削除する。（５）「パートモジュール」のメニューから「形状」→「切り取り」→「押し出し」を選択。「スケッチャー」ツールボックスの右上にある、「直線の作成」ツールを選択し、一続きの閉じた直線で削除領域を囲み、右クリックで「作業のキャンセル」、プロンプトエリアの「完了」をクリックして終了すると、解析領域の「パート」ができる。 9 2.3 材料特性の定義（１）モデルツリーの「材料特性」をダブルクリック。（２）ダイアログボックスで「名前」を epoxy と入力、メニューから「機械的」→「弾性」→ 「弾性」を選択、「ヤング率」に 4000、「ポアソン比」に 0.38 を入力。（３）「OK」をクリックで終了。 2.4 要素特性の定義と割り当て 2.3 で定義した材料特性を適用する要素特性を定義し，これを「パート： plate」に割り当てる．定義：（１）モデルツリーの「要素特性」をダブルクリック。（２）ダイアログボックスで「名前」を plateSection と入力し，「カテゴリ」から「ソリッド」，「タイプ」の一覧から「均質」，「続ける」をクリック．（３）「要素特性の編集」ダイアログボックスで「材料特性」として「epoxy」，「平面応力／ひずみ要素厚の値」で５を入力，「OK」をクリック．割り当て：（１）モデルツリーの「plate」と名づけた「パート」を展開（＋印をクリック）．（２）「要素特性の割り当て」をダブルクリックして，ビューポート上で要素特性を割り当てるパート領域をクリックすると，選択領域がハイライト表示される．（３）プロンプトエリアの「完了」をクリックして終了．「要素特性割り当ての編集」ダイアログボックスですでに定義した要素特性「plateSection」を選択し，「OK」をクリック． 2.5 モデルのメッシュ分割 2.4 で要素特性を定義し，それを割り当てた「plate」を適切な要素で分割する．まず、要素分割を適切に制御するために，以下の操作で「パート： plate」を対角上で 2 分割する．（１）「plate」の下にある「メッシュ」をダブルクリック．メニューから「ツール」→「パーティション」を選択．「パーティションの作成」ダイアログボックスでタイプとして「フェイス」，方法として「スケッチ」を選択．（２）「スケッチャー」ツールボックスの右上にある、「直線の作成：結合」ツールを選択し、孔の中心と板の右上を結ぶ直線を作成．（３）右クリック「作業のキャンセル」で終了。プロンプトエリアの「完了」をクリックして「スケッチャー」を終了。メッシュコントロールの割り当て（メッシュ分割の方法を指定）：（１）「メッシュモジュール」のメニューから「メッシュ」→「コントロール」を選択。（２）メッシュコントロールを割り当てるパート領域をビューポート上でクリックして選択すると，選択領域がハイライト表示される．プロンプトエリアの「完了」をクリック。（３）「メッシュコントロール」ダイアログボックスで「要素形状」から「四辺形」，「テクニック」から「構造」を選択．「OK」をクリック。（４）同様にもう一つのパート領域もメッシュコントロールを割り当てる。（５）プロンプトエリアの「完了」をクリックして終了。 10 要素タイプの選択：（１）メニューから「メッシュ」→「要素タイプ」を選択．（２）要素タイプを割り当てる領域をビューポート上でクリックして選択すると，選択領域がハイライト表示される．プロンプトエリアの「完了」をクリック。（３）「要素タイプ」ダイアログボックスで「要素ライブラリ」から「Standard」，「ジオメトリ次数」から「2 次」，「ファミリ」から「平面応力」を選択．「OK」をクリック。「パート：circle plate」に CPS8R 要素が選択される．（４）同様にもう一つの領域も要素タイプを割り当てる。（５）プロンプトエリアの「完了」をクリックして終了。要素分割：（１）メニューから「シード」→「パート」を選択，「全体シード」ダイアログボックスの「近似全体サイズ」に 2.0 を入力，「OK」をクリック．プロンプトエリアの「完了」をクリックして終了．ここで，節点が必ず配置される場所を正方形で，おおよそ節点が配置される場所を丸印で表される．（２）メニューから「メッシュ」→「パート」を選択，プロンプトエリアの「はい」をクリックするとメッシュ分割される． 2.6 モデルのアセンブリ要素分割された「パート：plate」をアセンブリする．（１）モデルツリーの「アセンブリ」を展開し，「インスタンス」をダブルクリック．（２）「インスタンスの作成」ダイアログボックスでパート「plate」を選択し，「OK」をクリック． 2.7 解析条件の設定弾性解析であるので静的解析ステップを設定する．（１）モデルツリーの「ステップ」をダブルクリック，「ステップの作成」ダイアログボックスで「プロシージャタイプ：一般」から「Static,General」を選択，「続ける」をクリック．（２）「ステップの編集」ダイアログボックスで「基本」タブの説明に plate with hole under tension と入力． 2.8 「OK」をクリック．境界条件の設定変位と荷重の境界条件を設定するために，設定箇所の定義し，必要な条件を設定する．集合とサーフェスの定義：（１）モデルツリーの「アセンブリ」の下の「集合」をダブルクリック，「集合の作成」ダアログボックスで「名前」を left と入力，「タイプ」はジオメトリを選択し，「続ける」をクリック．ビューポート上でモデルの左側エッジ（稜線）をクリックして選択し，プロンプトエリアの「完了」をクリック．（２）同様の手順でモデルの底辺を bottom として「集合」を定義する．（３）「アセンブリ」の下の「サーフェス」をダブルクリック，「サーフェスの作成」ダイアログ 11 ボックスで「名前」を pull と入力，「続ける」をクリック，モデルの上辺を選択し，プロンプトエリアの「完了」をクリック．変位境界条件の設定：（１）モデルツリーの「境界条件」をダブルクリック，「境界条件の作成」ダイアログボックスで「名前」を left と入力，「ステップ」として「Step-1」を選択，「カテゴリ」から「機械的」，「選択されたステップに対するタイプ」は「対称/反対称/完全固定」を選択，「続ける」．（２）プロンプトエリアの「集合」をクリック，「領域の選択」ダイアログボックスで left を選択し，「ビューポート内で選択をハイライトする」のチェックを有効にすれば選択されたエッジがハイライト表示される．（４）「続ける」をクリック，「境界条件の編集」ダイアログボックスで「XSYMM」を選択し，「OK」をクリックして終了すると，設定されたエッジに沿って記号が表示される．（５）同様に bottom としたエッジに「YSYMM」の条件を設定する．荷重境界条件の設定：（１）モデルツリーの「荷重」をダブルクリック，「荷重の作成」ダイアログボックスで「名前」を tension と入力，「ステップ」として「Step-1」を選択，「カテゴリ」から「機械的」，「選択されたステップに対するタイプ」から「圧力」を選択，「続ける」．（２）プロンプトエリアの「サーフェス」をクリック，「領域の選択」ダイアログボックスで pull を選択すると，サーフェスがハイライト表示される．（３）「続ける」をクリック，「荷重の編集」ダイアログボックスで大きさに -2.0 を入力し，「OK」をクリックして終了すると，指定されたサーフェスにベクトルが表示される． 2.9 解析ジョブ（解析の単位）の作成と投入（実行）（１）モデルツリーの「ジョブ」をダブルクリック，「ジョブの作成」ダイアログボックスで「名前」を plate-hole と入力，一覧から「tension」を選択し，「続ける」．（２）「ジョブの編集」ダイアログボックスではそのまま「OK」をクリック．（３）モデルツリーの「ジョブ」を展開し，plate-hole と名づけた「ジョブ」の上で右クリックして「ジョブの投入」を選択して，計算を実行させる． 2.10 解析結果の表示（１）ジョブが問題なく完了したら，ステータス表示は「完了」となり，解析結果を結果表示モジュールを使って表示させることができる．（２）ジョブ名を右クリックして「結果」を選択し，応力分布のコンター図を作成する．（３）メニューから「結果」→「フィールド出力」，表示させる出力変数を選択して，所要の結果を得る． 12 機械要素の応力解析はじめに各種機械や構造物を設計するとき、その構成部材が通常の使用内で破壊しないように材質や形状寸法が決定される。負荷による変形の大きさや破壊は部材内に生じる応力（単位面積当たりに作用する内力）によって評価される。つまり、垂直応力あるいはせん断応力が材料固有のある限界値に達すると破壊や永久変形が生じると考える。したがって、応力解析は設計上、非常に重要な問題である。一般に機械に用いられる部品（機械要素）の形状、作用する負荷、支持条件などは、おおよそ材料力学で扱ったものより複雑である。このような材料力学により応力や変形を評価することが困難な問題の解析に適用される手法として、実験的に解析する方法と計算機（コンピュータ）を用いた数値解析がある。実験的に解析する方法としてひずみゲージ法や光弾性法が代表的に知られている。ひずみゲージ法は、実部材表面上で測定したい箇所にひずみゲージを接着し、直接的に部材表面上のひずみを測定し、応力に換算する。２次元面の応力、ひずみ分布を測定するには大きな労力が必要である。一方、光弾性法は部材の表面の応力分布が得られ、複雑形状によって起こる応力集中の領域とその応力値を測定することができるのが特徴である。３次元光弾性法によれば、表面のみならず内部の応力分布も得られる。しかし、光弾性の性質を利用するために光弾性体である透明樹脂製の模型を用いるモデル実験であり、実部材の形状、負荷、支持条件などの再現に課題がある。計算機を用いた数値解析の中で、最も広く用いられるのが有限要素法（Finite Element Method、FEM）である。近年の計算機の進歩によって、機械工学の分野においては、設計、製造に CAD、CAM が普及し、より複雑形状な機械要素において最適な設計条件や性能予測が必要になる中、FEM は CAE（Computer Aided Engineering：計算機支援工学）として欠かすことができない重要なツールとなっている。１．目的実験手法（光弾性実験）と数値解析手法（FEM）により比較的簡単な形状の部材の二次元弾性応力解析を行い、機械要素や構造部材に生ずる応力や変形の評価方法を学ぶ。また、光弾性実験法および有限要素法の基礎、特徴を理解する。２．光弾性実験の概要 2.1 光弾性実験装置図１は光弾性実験装置を示す．光源部からフィルタ，絞りを通り出た光を視野レンズで平行光線にした光を偏光子（偏光レンズ）と 1/4 波長板で偏光にして，負荷を与えた光弾性試験片に垂直入射し，1/4 波長板と偏光子を通して光弾性縞を測定するものである．図２に円板の直径方向圧縮の光弾性縞を示す．明と暗から成る縞が見られる．負荷を大きくしていくと縞の数が増えていくのがわかる．試験片は光弾性材料（透明エポキシ樹脂など） 1 で作られた光弾性モデルであり，これに荷重を加えて，モデル内の応力分布の様相を測定するものである．図１光弾性実験装置 (a) 荷重 30kgf の等色線 (b) 荷重 60kgf の等色線 (c) 荷重 90kgf の等色線図２円板の圧縮の光弾性縞の例 2.2 光弾性効果無応力状態の光学的に等方な透明材料（光弾性材料）に外力を加えると，異方性すなわち複屈折性を帯びる．これを光弾性効果という．図３は主応力σ1 とσ2 の応力状態の光弾性材料に光が入射したときの複屈折性を相対的光路差 J で示したものである．すなわち，入射した光は主応力の方向に振動する二つの平面偏光となって進み，かつ両者の速度 v1 と v2 が異なるので相対的な位相の差δが生じる．無負荷の状態のモデルの屈折率を n0，負荷後における主応力σ1，σ2 の方向に振動する平面偏光の屈折率をそれぞれ n1，n2 とすれば，次の関係式が成り立つ． n1 ― n0＝c1σ1＋c2σ2 n2 ― n0＝c1σ2＋c2σ1 から n1 ― n2＝（c1-c2）（σ1－σ2）＝c（σ1－σ2） 2 図３主応力σ1、σ2 と相対的光路差Ｊの関係ここで c1，c2 および c（＝c1-c2）は材料，温度，使用光の波長により定まる光弾性定数である．屈折率は，真空中（空気中）の光の速度 v0 とその物質中の速度 v1，v2 の比で与えられるから，モデルの厚さを t とすれば相対的光路差 J は， J＝t（v0/v1－v0/v2）＝t（n1－n2）＝t c（σ1－σ2）光弾性実験では単色光を使用するので，単色光の波長λを用いて，J/λ＝Ｎを用いる．すなわち，光路差をλの何倍にあたるかで計量する．このようなＮを縞次数あるいはフリンジ数(fringe order)といい，次式のようになる．Ｎ＝(ｃ/λ)t（σ1－σ2）＝αt（σ1－σ2）この式は光弾性実験の最も基本になる式で，Ｎが計測できれば主応力差が求められることになる．ここで，α＝ｃ/λは光弾性（応力）感度と呼ばれる．３．光弾性感度αの検定等色線から主応力差を知るためには，光弾性感度αを測定する必要がある．ここでは，円板の直径に沿って圧縮する実験から求める．図４に示すように直径Ｄ＝2Ｒ、厚さ t の円板に中心Ｏを通り直角の xy 軸をとり、y 軸方向の直径に集中荷重Ｗを加えるとき、円板中心Ｏにおける主応力は弾性理論より、 σ1＝σx＝２Ｗ／πtＤ σ2＝σy＝－６Ｗ／πtＤよって σ1－σ2＝８Ｗ／πtＤしたがって、αは、 α＝Ｎ／t（σ1－σ2）＝πＤＮ／８Ｗ＝（πＤ／８）・（Ｎ／Ｗ） (mm/N) 円板の中央において、縞次数と荷重との関係を求め、図５のような感度検定曲線からαは求まる。図４円板の圧縮図５感度検定曲線の例 3 ４．円板圧縮の等色線の測定部材内に生じている応力分布を示すために、等色線図からただちに主応力差分布が得られる。本実験では、等色線図を観察し、主応力差分布の様相を知ることにより応力分布を考察する。・円板圧縮の荷重Ｗが（）kgf のときの等色線図をスケッチする。・等色線図から、x 軸上におけるフリンジ分布図を作図する。これにより x 軸上の主応力差分布がわかる。５．考察円板圧縮の x 軸上の主応力差分布から主応力分布の状態を推定し、主応力分布について考察、検討する。６．研究課題（１）ＦＥＭによる応力解析から得られた応力分布と比較する。（２）光弾性実験の特徴を述べなさい。（３）主応力について述べなさい。７．参考文献例えば、山田貴博（訳）、有限要素法 ABAQUS(Student Edition)付 4 丸善