
CAEソフト【 Femtet 】-ムラタソフトウェア株式会社
亀裂の入った板を引っ張た時の応力を解析します。
計算の結果得られたひずみエネルギーからエネルギー解放率と応力拡大係数を計算します。
これらの値は破壊強度を評価する時に使用します。
同様の解析を異なる方法で解析した例題46 J積分を用いた亀裂先端のエネルギー解放率解析も
参考にしてください。
表に記載されていない条件は初期設定の条件を使用します。
項目 |
条件 |
解析空間 |
2次元 |
モデル単位 |
mm |
項目 |
条件 |
ソルバー |
応力解析[Galileo] |
解析の種類 |
静解析 |
解析オプション |
選択なし |
下図のように中央に亀裂が入った板を引っ張ったときの応力を計算します。
強く引っ張った時に、先端から破壊されますが、この時の破壊強度を評価
するのに応力拡大係数Kが使用されます。このモデルの場合のKは、次式の
理論解が知られています。
この応力拡大係数を、Femtetを使って計算してみます。
出展)
鷲津久一郎他、有限要素法ハンドブックⅡ応用編、培風館, p.329
モデルは左右上下対称なので、右上の1/4領域だけを解析します。
対称面には、「対称/不連続タブ」で「対称面」を設定します。
ただし、亀裂面は開放状態であるため、対称面を設定する必要はありません。
BODY No./BODY Type |
BODY名 |
材料名 |
0/Sheet |
Plate |
Fe |
ボディ属性で、シートボディの厚みを1mにします
ボディ属性名 |
タブ |
設定 |
Plate |
厚み/幅 |
シートボディの厚み1000 |
材料定数で、ヤング率とポアソン比を設定します
材料名 |
ヤング率 |
ポアソン比 |
PlateMat※ |
70.56×109[Pa] |
0.33 |
※材料定数は、上記記載の出典で使用しているヤング率とポアソン比を設定しています。
境界条件名/Topology |
タブ |
境界条件の種類 |
条件 |
Sym_X/Edge |
対称/不連続 |
対称 |
対称面 |
Sym_Z/Edge |
機械 |
対称 |
対称面 |
P/Edge |
機械 |
圧力 |
-98×106[Pa] |
変形図に応力の最大主応を重ねた図を示します。
コンター図の表示範囲を0~400MPaに設定しています。
対称面境界を設定しているので、全体モデルが表示されています。
亀裂の先端に強い応力が発生していることがわかります。
亀裂先端の応力は、理論的には無限大ですが、実際には大きい値になっており、その値は
メッシュ分割によって敏感に変化します。したがって、ここの応力の値で破壊強度を評価する
ことはできません。破壊強度の評価には、エネルギー解放率gや応力拡大係数Kなどが使用
され、その値の大小で破壊強度を評価します。
エネルギー解放率gは、ひずみエネルギーUと亀裂の長さaの変化率で次式で定義されます。
(1)
ここではaの変化を0.1mmとし、a=4.05とa=3.95の場合のひずみエネルギーを計算します。
上の1/4モデルで、a=4.05[mm]のとき、U=12.1805977[J]、a=3.95のとき、
12.0728236であるので、全体モデルでは
(2)
分子に4をかけたのは、全体モデルのエネルギーは1/4モデルの4倍になるからです。
また、分母に2をかけたのは、全体モデルで亀裂の長さは2倍変化するからです。
2次元の場合は、エネルギー解放率gと応力拡大係数Kとは次の関係があります。
平面応力
(3)
平面ひずみ
(4)
今回の場合、平面応力なので式(3)を使用し、式(2)の結果から応力拡大係数は
(5)
この値は、冒頭の理論解と2%程度の差で一致しています。
まずはFemtetを試してみたい
試用版・無償版はこちらもっとFemtetについて詳しく知りたい
イベント・セミナー情報はこちら