損傷モデル選択に関する基本的な基準ランダム応答疲労解析
ランダムな荷重を受ける構造物の疲労寿命に関するスタディです。
ランダム応答解析から得られたパワースペクトル密度(PSD)の結果は、応力範囲に対するサイクル数の確率密度関数を生成するために使用されるモーメントmn の計算に使用されます。
PSDモーメントは、ランダム応答解析から生成されたPSD応力をもとに算出されます。
パワースペクトル密度(PSD)モーメント

以下を元にモーメントを算出します:
- fk
- 周波数値。
- Gk
- 周波数 fk での PSD 応答値。
応力範囲の発生確率の計算
各Binセクション内の初期応力範囲値と最終応力範囲値の間の応力範囲の発生確率をユーザー定義で計算します。
(ΔSi − δS/2) と(ΔSi + δS/2) の間に発生する応力範囲の確率P(ΔSi) :
確率密度関数(応力範囲に対するサイクル回数の確率密度)
上記のように算出されたPSDモーメントは、応力範囲の確率密度関数 f(mn) の生成に使用されます。この関数は、指定された損傷モデルに基づいています。損傷モデルを定義するために、DIRLIK、LALANNE、NARROW、Steinberg 3 band(THREE)のオプションが用意されています。
- DIRLIK
- DIRLIKは、確率密度関数を決定するための閉形式の解を次のように仮定します:p(S) = D1Qe−ZQ + D2ZR2e−Z22R2 + D3Ze−Z222√m0
- LALANNE
- LALANNEランダム疲労損傷モデルは、確率密度関数を次のように表現しています:p(S) = 1√m0 √1 − γ2√2πe−S28m0(1 − γ2) + Sγ4√m0(1 + erf(Sγ2√2m0(1 − γ2)))
- NARROW
- 狭帯域ランダム疲労損は、以下の確率関数を使用します:p(S) = (S4m0e−(S28m0))
- THREE
- Steinberg 3-band ランダム疲労損傷モデルは、以下の確率関数を使用します:P(S) = {0.683 at 2√m00.271 at 4√m00.043 at 6√m0
上限応力範囲係数
応力範囲の上限の算出:
応力範囲上限値 = 2 * RMS応力 * 応力範囲上限値係数
ランダム応答サブケースからRMS応力を出力します。対象となる応力範囲は、応力範囲の上限値で制限されます。上限を超える応力は、ランダム疲労損傷の計算では考慮されません。
上限応力範囲
上限応力範囲を直接指定します。
Bin数
確率が算出される応力範囲 DS = δS の幅を算出する(図 2参照)。デフォルトは100で、最初のBinは0.0からδS までです。応力範囲の幅は、DS=上限応力範囲/Bin数として計算されます。
応力範囲幅
応力範囲幅( δS )を直接定義します。
応力範囲の発生確率の計算
- DIRLIK, LALANNE, NARROW
-
(ΔSi − δS/2)
と
(ΔSi + δS/2)
の間に発生する応力範囲の確率
P(Si)
: P(Si) = pi(Si)δS
損傷モデルの選択
- 応力のPSDモーメントから対応するモーメントを算出し、これを用いて応力範囲の確率密度関数を決定します。
- DIRLIKとLALANNEモデルは、より広い応力範囲スペクトラム分布での確率を生成します。したがって、これらのモデルは、入力されるランダム信号が複数の周波数にわたる様々な応力範囲で構成されている場合に適しています。DIRLIKとLALANNEを使用した場合の確率密度関数の情報は、応力範囲分布のより広い範囲を把握できます。
- NARROWモデルは、応力範囲が特定の応力範囲分布と高い確率で密接に関連すると予想されるランダム信号に使用使用することを想定しています。したがって、入力されたランダムデータの応力範囲分布が広くなく、主に特定の応力範囲に分布が集中している場合、NARROWの使用が適しています。このモデルでは、特定の応力範囲に最も高い確率で応力範囲が存在すると予想されます。
- THREEモデルはNARROWモデルと似ていますが、ランダム信号の分布に、1*RMSとの関連に加え、2*RMS、3*RMSとの関連もある程度含まれると予想されます。したがって、入力されたランダムデータが主に1*RMSの応力範囲に集まっており、さらに少ないながらも2*RMS、3*RMSに関係してくる場合は、THREEの使用が適しています。
ピーク数とゼロクロッシング
- NARROW、THREE
- オリジナルの時間領域ランダム荷重(ここから周波数ベースのランダムPSD荷重が生成される)の1秒あたりのゼロクロッシングの数は、次のように決定されます:nzcross = √m2m0ここで、
- mn
- 対応するモーメントが計算されます。
- DIRLIK、LALANNE
- オリジナルの時間領域ランダム荷重(ここから周波数ベースのランダムPSD荷重が生成される)の1秒あたりのピーク数は、次のように決定されます:npeaks = √m4m2ここで、
- mn
- 対応するモーメントが計算されます。
サイクル数
- NARROW、THREE
- サイクルの総数を次の式から求めます:NT = nzcrossTここで、
- T
- 暴露時間。
- DIRLIK、LALANNE
- サイクルの総数を次の式から求めます:NT = npeaksTここで、
- T
- 暴露時間。
- 特定応力範囲における総サイクル数
- 応力範囲
ΔSi
を持つ総サイクル数は、次のように計算されます:Ni = P(ΔSi)NT
疲労寿命と損傷
疲労寿命(材料が破壊するまでの特定の応力範囲 Si の最大サイクル数)は、材料SN曲線に基づいて次のように計算されます:
適用されたランダム荷重の結果としてのトータルの疲労損傷は、次のように計算されます:
平均応力がゼロとならない荷重による平均応力補正を考慮するために、そのような荷重(通常は重力荷重)からなる静的サブケースを定義することができます。この静的サブケースは、ランダム疲労解析のイベント設定で参照することができます。