損傷モデル選択に関する基本的な基準ランダム応答疲労解析

ランダムな荷重を受ける構造物の疲労寿命に関するスタディです。

ランダム応答解析から得られたパワースペクトル密度(PSD)の結果は、応力範囲に対するサイクル数の確率密度関数を生成するために使用されるモーメントmn の計算に使用されます。

PSDモーメントは、ランダム応答解析から生成されたPSD応力をもとに算出されます。

パワースペクトル密度(PSD)モーメント

PSDモーメントは、ランダム応答解析で生成されたPSD応力をもとに、次のように計算されます:


Figure 1. PSDモーメントの計算

以下を元にモーメントを算出します:

mn=Nk=1fnkGkδf

ここで、
fk
周波数値。
Gk
周波数 fk での PSD 応答値。

応力範囲の発生確率の計算

各Binセクション内の初期応力範囲値と最終応力範囲値の間の応力範囲の発生確率をユーザー定義で計算します。

(ΔSiδS/2) (ΔSi+δS/2) の間に発生する応力範囲の確率P(ΔSi)

P(ΔSi)=pδSi

確率密度関数(応力範囲に対するサイクル回数の確率密度)

上記のように算出されたPSDモーメントは、応力範囲の確率密度関数 f(mn) の生成に使用されます。この関数は、指定された損傷モデルに基づいています。損傷モデルを定義するために、DIRLIK、LALANNE、NARROW、Steinberg 3 band(THREE)のオプションが用意されています。

DIRLIK
DIRLIKは、確率密度関数を決定するための閉形式の解を次のように仮定します:
p(S)=D1QeZQ+D2ZR2eZ22R2+D3ZeZ222m0
ここで、
D1=2(xmγ2)1+γ2
D2=1γD1+D211R
D3=1D1D2
Z=S2m0
Q=1.25(γD3D2R)D1
R=γxmD211γD1+D21
xm=m1m0m2m4
γ=m2m0m4
不規則性係数。
S
応力範囲。
LALANNE
LALANNEランダム疲労損傷モデルは、確率密度関数を次のように表現しています:
p(S)=1m01γ22πeS28m0(1γ2)+Sγ4m0(1+erf(Sγ22m0(1γ2)))
ここで、
γ=m2m0m4
不規則性係数。
S
応力範囲。
NARROW
狭帯域ランダム疲労損は、以下の確率関数を使用します:
p(S)=(S4m0e(S28m0))
ここで、
S
応力範囲。
デフォルトでは、 狭帯域には、ゼロクロッシング数 npeaks=m4/m2 はピーク数 nzcross=m2/m0 の代わりに使用されます。これは、 m4 を含む数値計算が不安定な数値挙動につながることがあるためです。理想的な 狭帯域であれば、ゼロクロッシングの数とピークの数はほぼ同じになるはずです。
THREE
Steinberg 3-band ランダム疲労損傷モデルは、以下の確率関数を使用します:
P(S)={0.683at2m00.271at4m00.043at6m0
ここで、
S
応力範囲。
他の損傷モデルとは異なり、THREEバンドでは、これらの値は確率(確率密度ではない)です。これは、他の損傷モデルの下域 p(S) に対して、上域 P(S) を使用していることからも明らかです。
THREE損傷モデルの場合、これらの確率は、 P(S) を、全時間履歴におけるゼロクロッシングの総数とかけることで、サイクル数を計算するために直接使用されます。他のダメージモデル(THREEを除く)では、まず確率密度値に DS (Binサイズ)をかけることで確率を求めます。


Figure 2. 確率密度関数. 応力範囲に対するサイクル回数の確率密度
確率密度関数は、ランダム応答疲労解析の設定で定義された以下のパラメータに基づいて調整することができます。

上限応力範囲係数

応力範囲の上限の算出:

応力範囲上限値 = 2 * RMS応力 * 応力範囲上限値係数

ランダム応答サブケースからRMS応力を出力します。対象となる応力範囲は、応力範囲の上限値で制限されます。上限を超える応力は、ランダム疲労損傷の計算では考慮されません。

上限応力範囲

上限応力範囲を直接指定します。

Bin数

確率が算出される応力範囲 DS=δS の幅を算出する(図 2参照)。デフォルトは100で、最初のBinは0.0からδS までです。応力範囲の幅は、DS=上限応力範囲/Bin数として計算されます。

応力範囲幅

応力範囲幅( δS )を直接定義します。

応力範囲の発生確率の計算

各Binセクション内の初期応力範囲値と最終応力範囲値の間の応力範囲の発生確率は、損傷モデルに基づいて計算されます。
DIRLIK, LALANNE, NARROW
(ΔSiδS/2) (ΔSi+δS/2) の間に発生する応力範囲の確率 P(Si)
P(Si)=pi(Si)δS
THREE
式 6をご参照ください。
THREE損傷モデルの場合、Binは3つのみです。各応力範囲(2*RMS、4*RMS、6*RMS)におけるサイクル数は、対応する確率にゼロクロッシングの総数を直接乗じて算出します(ゼロクロッシングの算出については下記をご参照ください)。

損傷モデルの選択

損傷モデル選択に関する基本的な基準:
  1. 応力のPSDモーメントから対応するモーメントを算出し、これを用いて応力範囲の確率密度関数を決定します。
  2. DIRLIKとLALANNEモデルは、より広い応力範囲スペクトラム分布での確率を生成します。したがって、これらのモデルは、入力されるランダム信号が複数の周波数にわたる様々な応力範囲で構成されている場合に適しています。DIRLIKとLALANNEを使用した場合の確率密度関数の情報は、応力範囲分布のより広い範囲を把握できます。
  3. NARROWモデルは、応力範囲が特定の応力範囲分布と高い確率で密接に関連すると予想されるランダム信号に使用使用することを想定しています。したがって、入力されたランダムデータの応力範囲分布が広くなく、主に特定の応力範囲に分布が集中している場合、NARROWの使用が適しています。このモデルでは、特定の応力範囲に最も高い確率で応力範囲が存在すると予想されます。
  4. 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曲線に基づいて次のように計算されます:

Nf(Si)=(SiSf)1b

適用されたランダム荷重の結果としてのトータルの疲労損傷は、次のように計算されます:

D=Ni=1NiNf(Si)

平均応力がゼロとならない荷重による平均応力補正を考慮するために、そのような荷重(通常は重力荷重)からなる静的サブケースを定義することができます。この静的サブケースは、ランダム疲労解析のイベント設定で参照することができます。