単軸疲労解析
周期荷重を受ける構造の寿命(荷重サイクル数)を予測する単軸疲労解析が、S-N(応力-寿命)、E-N(ひずみ-寿命)アプローチを用いて、OptiStructで実行できます。
応力-寿命法は、構造の応力レベルがほとんど弾性範囲に収まる際の疲労寿命予測に機能します。このような繰り返し荷重条件の下では、構造は一般的に多数の荷重サイクルに耐えることができ、これは高サイクル疲労として知られています。繰り返しによるひずみが塑性範囲まで延びた場合、構造の疲労耐久性は著しく低下し、これは低サイクル疲労として特性づけられます。高サイクルと低サイクル疲労の間の一般的に受け入れられている境界点はおよそ10,000荷重サイクルぐらいになります。低サイクル疲労では、ひずみ-寿命 (E-N) 法が適用され、塑性ひずみがダメージ計算の重要な因子として考慮されます。
疲労解析が実行されるモデルのセクションは、FATDEFバルクデータエントリで特定される必要があります。適当なFATDEFバルクデータエントリが疲労サブケース定義からFATDEFサブケース情報エントリを通して参照される必要があります。
応力-寿命(S-N)アプローチ
S-N曲線
セグメント1に対して
ここで、公称応力のレンジ、は破壊する疲労サイクル、は最初の強度指数、は疲労強度指数です。
S-N アプローチは弾性的な繰り返し荷重に基づき、S-N 曲線は寿命軸上の、1000 サイクル以上の数の範囲にあるはずとの推定に基づいています。このことは大きな塑性は起こらない、ということを保証します。これは一般に高サイクル疲労と言われます。
与えられた材料のS-N曲線は、MATFATバルクデータエントリで用意します。これは、構造材料定義と材料IDによって共有される (MID) を通して参照されます。
レインフローサイクルカウント
サイクルカウントはランダム荷重履歴から離散的で単純な "等価"な一定振幅サイクルを取り出すために用いられます。"サイクルカウント"を理解する1つの方法は応力-ひずみ関係を、時刻との対比に変更したものとして捉えることです。サイクルカウントは応力-ひずみ関係のヒステリシスループの数をカウントし、その範囲/平均値または最大値および最小値の記録をつけます。
- 単純な荷重履歴:
この荷重履歴は連続的であるため、山と谷のみから成る荷重履歴に変換されます。荷重履歴においてポイント4は応力の山(最大)であることが明らかであり、再配列中に前面に移動されます(図 5)。再配列後、山と谷は便宜のため再番号付けされます。
次に、最初の3つの応力値(1、2および3)をピックし、サイクルが存在するかを見極めます。
が応力値、ポイントを表すとすると:(2) (3) であるため、ポイント2から3までのサイクルが抽出されます。これで1つのサイクルが抽出され、2つのポイントがグラフから消去されます。残りのポイントについて、同じプロセスが適用されます:(6) (7) この場合、であるため、ポイント1から4までの別のサイクルが抽出されます。これら2つのポイントも破棄された後は、ポイント5だけが残ります; したがって、レインフローカウントプロセスは完了となります。
2つのサイクル(2→3および1→4)が、この荷重履歴から抽出されました。最も高い山 / 谷を選択し、荷重履歴を再配列する主な理由の1つは、最も大きなサイクルが常に抽出される(この場合は1→4)ことを確実にすることにあります。再配列の前に荷重履歴を見て同じレインフローカウントプロセスを実行すると、1→4のサイクルが抽出されないことは明らかです。
- 複雑な荷重履歴レインフローカウントプロセスは、荷重履歴のポイント数にかかわらず同じです。しかしながら、再配列に使用される最も高い山 / 谷の位置によっては、再配列プロセスがどのように行われるかがはっきりしない場合もあります。図 7 は、より複雑な荷重履歴の再配列プロセスを示しています。後続のレインフローカウントは、上記の単純な例で挙げているプロセスの外挿に過ぎず、ここでは繰り返しません。この荷重履歴は連続的であるため、山と谷のみから成る荷重履歴に変換されます。明らかに、荷重ポイント11は荷重の最大値であり、したがって、荷重履歴はここで再配列、再番号付けされます。
荷重履歴は、最大荷重後も含めすべてのポイントが荷重履歴の最初に移動され、荷重履歴の最後から消去されるよう再配置されます。
レインフローサイクルカウントに関係するパラメータは、FATPARMバルクデータエントリで定義することができます。疲労サブケース定義からFATPARMサブケース情報エントリを通して、適切なFATPARMバルクデータエントリを参照することができます。
等価相当応力
S-N 理論は単軸応力を扱うため、それぞれの計算点、それぞれの時間ステップにおける応力成分は1つの統合された値にする必要があり、そこで、等価な公称応力が用いられ、S-N 曲線に適用されます。
さまざまな応力の組み合わせタイプが利用可能で、デフォルトは “Absolute maximum principal stress (絶対値最大の主応力)”です。脆性材料には "Absolute maximum principal stress"が推奨されますが、延性材料には"Signed von Mises stress"が推奨されます。符号付きパラメータの符号には最大絶対値主値の符号が用いられます。
応力の組み合わせに関係するパラメータは、FATPARMバルクデータエントリで定義できます。疲労サブケース定義からFATPARMサブケース情報エントリを通して、適切なFATPARMバルクデータエントリを参照することができます。
平均応力補正
一般的に S-N 曲線は完全反転の繰り返し荷重の標準の試験から得られます。しかし、本来の疲労荷重は、完全なものではなく、平均応力は部品の疲労性能上で大きな影響を与えます。疲労強度に関して引張りの軸平均応力は有害ですが、圧縮の軸平均応力は有益です。平均応力修正は0でない平均応力の影響を考慮するために用いられます。
Haigh'座標のGerber放物線とGoodman線が平均応力の影響を考慮する際に広く用いられており、これらは次のように表すことができます:
- 下記で与えられる平均応力;
- 下記で与えられる応力範囲;
- 平均応力補正後の応力範囲(応力範囲には、平均応力には)
- 最大強度
Gerber法では正と負の平均応力修正を平均応力が疲労破壊を加速するものとして同じ方法で取り扱うのに対し、Goodman法は負の平均応力を無視します。どちらの方法とも圧縮平均応力に対して保守的な結果を与えます。Goodman法は脆性材料に対して推奨されるのに対し、Gerber法は延性材料に対して推奨されます。Goodman法では、引張り平均応力が UTSより大きい場合、損傷は1.0より大きくなります。Gerber法では、平均応力がUTSより大きい場合、引張または圧縮のいずれでも、損傷は1.0より大きくなります。
平均応力の影響に関係するパラメータは、FATPARMバルクデータエントリで定義することができます。疲労サブケース定義からFATPARMサブケース情報エントリを通して、適切なFATPARMバルクデータエントリを参照することができます。
FKM:
OptiStructにおけるFKM補正には2つのオプションが使用でき、それらはFATPARMエントリでUCORRECTをFKM/FKM2に、もしくはMCORRECT(MCi)フィールドをFKMに設定することでアクティブ化されます。
- Regime 1 (R > 1.0)
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- Regime 4 (R ≥ 0.5)
- 平均応力補正後の応力振幅(耐久応力)
- 平均応力
- 応力振幅
- Regime 1 (R > 1.0)かつRegime 4 (R ≥ 0.5)
- 平均応力補正は適用されません
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- 平均応力補正後の応力振幅(耐久応力)
- 平均応力
- 応力振幅
- = MSS2
平均応力補正について4つすべてのMSSiフィールドが指定されている場合、対応する平均応力感度値は4つすべてのレジュームを制御する勾配です。FKM-Guidelinesに基づき、Haigh図は応力比()の値をベースに4つのレジームに分割されます。続いて、補正値を使って損傷および寿命計算ステージのためにS-Nカーブが選択されます。
OptiStructにおけるFKM補正には2つのオプションが使用でき、それらはFATPARMエントリでUCORRECTをFKM/FKM2に、かつ、MCORRECT(MCi)フィールドをFKMに設定することでアクティブ化されます。
- Regime 1 (R > 1.0)
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- Regime 4 (R ≥ 0.5)
- 平均応力補正後の応力振幅(耐久応力)
- 平均応力
- 応力振幅
- = MSSi
- Regime 1 (R > 1.0)かつRegime 4 (R ≥ 0.5)
- 平均応力補正は適用されません
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- 平均応力補正後の応力振幅(耐久応力)
- 平均応力
- 応力振幅
- = MSSi
- Regime 1 (R > 1.0)
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- Regime 4 (R ≥ 0.5)
損傷加算モデル
- ある応力の大きさの組み合わせと平均応力レベルにおけるS-N曲線からの材料の疲労寿命(破壊のサイクル数)。
- 荷重レベルにおける応力のサイクル数。
- 荷重サイクルにおける累積ダメージです。
線形損傷総和則は、周期疲労荷重による損傷の加算で荷重順序の効果を考慮しません。しかしながら、多くの応用でうまく機能することが実証されてきています。
ひずみ-寿命(E-N)アプローチ
ひずみ-寿命解析は、切欠きの根元のような多くの重大な場所では応力集中を伴うという事実に基づいており、そしてこれらは疲労破壊に達する前の周期荷重の間に明らかな塑性変形を伴います。このため、弾塑性ひずみの結果がひずみ-寿命解析の実行に不可欠になります。
Neuber補正
Neuber補正は弾性ひずみの結果を弾塑性の結果に補正する最も有名な手法です。
ここで、、は弾性解析で得られる局所弾性応力と局所弾性ひずみで、、は塑性ひずみが存在するときの応力とひずみです。とは、周期応力-ひずみ曲線とヒステリシスループのための式と共に式 15から計算できます。
単調増加の応力-ひずみ挙動
ここで、は現在の断面積、は現在の試験体の長さ、は試験体の初期長さ、とはそれぞれ真応力と真ひずみで、図 12は真応力-真ひずみ空間での単調増加の応力-ひずみ曲線を示します。全ての過程において、試験体がCで破壊するまで応力は増加し続けます。
繰り返し応力-ひずみ曲線
- 安定状態
- 周期硬化
- 周期軟化
- ひずみレンジによって軟化または硬化
- 周期強度係数
- ひずみ周期硬化指数
ヒステリシスループの形状
平均応力補正
実際には平均応力は避けられないにもかかわらず、実験室での疲労試験は常に完全反転で実施されます、このため、完全反転の実験で確立された疲労則は工学問題に適用される前に補正される必要があります。
Morrowの式は、低い塑性ひずみでは平均応力の影響が顕著で、高い塑性ひずみでは小さいという考察に一致しています。
SWT法では最大応力が0または負の時、損傷は0と予測することになりますが、これは現実と一致しません。
2つの方法を比較すると、SWT法は荷重で引張りが大部分の場合、保守的な寿命を予測するのに対し、Morrowのアプローチは荷重で圧縮が支配的な場合にも現実的な結果をもたらします。
損傷加算モデル
E-Nアプローチでは、S-Nアプローチと同じ損傷加算モデルを用いるものとし、Palmgren-Minerの線形総和則が用いられます。