ビーム要素(/PROP/BEAM, /PROP/INT_BEAM)

Radiossでは2種類のビーム要素が利用可能で、1次元構造とフレームに用いられます。これは軸力、せん断力、曲げ、およびねじりモーメントを伝達します(反対にトラスは軸力のみ支持します)。

古典的ビーム(/PROP/BEAM

デフォルトの定式化はTimoshenko定式化に基づいているため、横せん断ひずみが考慮されます。この定式化は横せん断エネルギーを無視する通常のEuler-Bernoulli定式化に縮退させることができます。

節点1と2が局所x-軸の定義に用いられます。局所y-軸はx-軸に直角で、t=0における節点1、2、3の平面内に置かれます。次に、平均のx-軸周りの回転を考慮して、それぞれのサイクルで位置は修正されます。局所 z-軸は右手則を用いて得られます。


図 1. ビーム要素と局所座標系
Radiossではビーム形状はその断面面積と3つの断面積慣性モーメントで定義されます。局所 Y軸と Z軸に関する面積慣性モーメントは曲げのためで、次の式を用いて計算できます:(1)
I y = A z 2 d y d z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG5baabeaakiabg2da9maapifabaGaamOEamaaCaaaleqa baGaaGOmaaaakiaadsgacaWG5bGaamizaiaadQhaaSqaaiaadgeaae qaniabgUIiYlabgUIiYdaaaa@439C@
(2)
I z = A y 2 d y d z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG6baabeaakiabg2da9maapifabaGaamyEamaaCaaaleqa baGaaGOmaaaakiaadsgacaWG5bGaamizaiaadQhaaSqaaiaadgeaae qaniabgUIiYlabgUIiYdaaaa@439C@

局所X軸に関する面積慣性モーメントはねじりのためです。これは単純に、 I y MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG6baabeaaaaa@37F0@ I z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG6baabeaaaaa@37F0@ の和として得ることができます。ねじりモデルは中実断面で反りが無視できる場合のみ有効です。

ビーム要素の最小時間ステップは次の式を用いて決められます:(3)
Δ t = a L c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeiLdiaads hacqGH9aqpdaWcaaqaaiaadggacaWGmbaabaGaam4yaaaaaaa@3BBF@
ここで、
c
音速
E / ρ
a = 1 2 min ( min ( 4 , 1 + b 12 ) F 1 , b 3 F 2 )
F 1 = 1 + 2 d 2 2 d
F 2 = min ( F 1 , 1 + 2 d s 2 2 d s )
b = A L 2 max ( I y , I z )
d = max ( d m , d f )
d s = d max ( 1 , 12 b 1 + 12 E 5 6 G b ( 1 I s h e a r ) )
ビーム断面を定義するユーザー入力パラメータは3つの面積慣性モーメントと断面積です。安定性と精度のために、次の制限を尊重することが推奨されます:(4)
L > A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitaiabg6 da+maakaaabaGaamyqaaWcbeaaaaa@38B1@
  • 0.01 A 2 < I y < 100 A 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGimaiaac6 cacaaIWaGaaGymaiaadgeadaahaaWcbeqaaiaaikdaaaGccqGH8aap caWGjbWaaSbaaSqaaiaadMhaaeqaaOGaeyipaWJaaGymaiaaicdaca aIWaGaamyqamaaCaaaleqabaGaaGOmaaaaaaa@4279@
  • 0.01 A 2 < I z < 100 A 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGimaiaac6 cacaaIWaGaaGymaiaadgeadaahaaWcbeqaaiaaikdaaaGccqGH8aap caWGjbWaaSbaaSqaaiaadMhaaeqaaOGaeyipaWJaaGymaiaaicdaca aIWaGaamyqamaaCaaaleqabaGaaGOmaaaaaaa@4279@
  • 0.1 ( I y + I z ) < I x < 10 ( I y + I z ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGimaiaac6 cacaaIXaWaaeWaaeaacaWGjbWaaSbaaSqaaiaadMhaaeqaaOGaey4k aSIaamysamaaBaaaleaacaWG6baabeaaaOGaayjkaiaawMcaaiabgY da8iaadMeadaWgaaWcbaGaamiEaaqabaGccqGH8aapcaaIXaGaaGim amaabmaabaGaamysamaaBaaaleaacaWG5baabeaakiabgUcaRiaadM eadaWgaaWcbaGaamOEaaqabaaakiaawIcacaGLPaaaaaa@4A7C@
ビーム要素には材料則1と2のみ利用可能です。材料則2では、内力の関数にグローバル塑性モデルが使用されます。主要な仮定は、ビーム断面が中実で矩形であることです。断面と断面のモーメントの間の最適な関係は:
  • 12 I y I z = A 4 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGymaiaaik dacaWGjbWaaSbaaSqaaiaadMhaaeqaaOGaamysamaaBaaaleaacaWG 6baabeaakiabg2da9iaadgeadaahaaWcbeqaaiaaisdaaaaaaa@3E2A@
  • I x = I y + I z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG4baabeaakiabg2da9iaadMeadaWgaaWcbaGaamyEaaqa baGccqGHRaWkcaWGjbWaaSbaaSqaaiaadQhaaeqaaaaa@3DDB@

このモデルは円形または楕円断面に対しても良い結果を与えます。薄肉断面に対しては、グローバル塑性モデルは誤った結果を与えるかもしれません。フレーム構造の1つのラインに対し1つのビーム要素を用いることは推奨しません。質量は節点に集中されます。そのため、正しい質量の分布を得るためには、細かいメッシュが必要です。これは特に動的効果が重要な場合には顕著になります。

加えて、 Radiossビーム要素では、モーメントがビームの長さ方向に変化しません。モーメントは一定であるという仮定の下でビームの中心で評価されます。応力も同様です。

その結果として、固定された単純梁では、モーメントがビームのつけ根ではなく中央で計算されるために、ビームに若干高い力を生みます。
注: ビーム要素の出力は局所座標系で表現されます。いくつかの結果において、それが平均X回転の考慮によって更新されるために、混乱を生じるかもしれません。例えば、ビームの1節点が完全に固定され、速度Vが他節点に与えられている場合、ビームはVの速度で回転しますが、局所座標系は速度V/2で回転します。これは、特にせん断力や曲げモーメントに間違った解釈を引き起こすかもしれません。

新しいビーム (/PROP/INT_BEAM

要素の断面は100 までの積分点を用いて定義されます(図 2)。要素の断面特性、面積慣性モーメントと面積はRadiossによって次のように計算されます:(5)
A = A i = ( d y i d z i ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaiabg2 da9maaqaeabaGaamyqamaaBaaaleaacaWGPbaabeaaaeqabeqdcqGH ris5aOGaeyypa0ZaaabqaeaadaqadaqaaiaadsgacaWG5bWaaSbaaS qaaiaadMgaaeqaaOGaamizaiaadQhadaWgaaWcbaGaamyAaaqabaaa kiaawIcacaGLPaaaaSqabeqaniabggHiLdaaaa@4670@
(6)
I z = A i ( y i 2 + 1 12 d y i 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG6baabeaakiabg2da9maaqaeabaGaamyqamaaBaaaleaa caWGPbaabeaakmaabmaabaGaamyEamaaBaaaleaacaWGPbaabeaakm aaCaaaleqabaGaaGOmaaaakiabgUcaRmaalaaabaGaaGymaaqaaiaa igdacaaIYaaaaiaadsgacaWG5bWaaSbaaSqaaiaadMgaaeqaaOWaaW baaSqabeaacaaIYaaaaaGccaGLOaGaayzkaaaaleqabeqdcqGHris5 aaaa@48BE@
(7)
I y = A i ( z i 2 + 1 12 d z i 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaWG5baabeaakiabg2da9maaqaeabaGaamyqamaaBaaaleaa caWGPbaabeaakmaabmaabaGaamOEamaaBaaaleaacaWGPbaabeaakm aaCaaaleqabaGaaGOmaaaakiabgUcaRmaalaaabaGaaGymaaqaaiaa igdacaaIYaaaaiaadsgacaWG6bWaaSbaaSqaaiaadMgaaeqaaOWaaW baaSqabeaacaaIYaaaaaGccaGLOaGaayzkaaaaleqabeqdcqGHris5 aaaa@48BF@
ビームモデルは横せん断ひずみを考慮してねじりの反りはないTimoshenko理論に基づいています。深いビーム(短いビーム)の場合に用いることができます。断面積に複数の積分点を用いることで、それぞれの積分点で von Mises基準での弾塑性モデルを得ることを可能にし、古典的なビーム要素と異なり、断面は部分的に塑性化することができます(TYPE3)。材料則1と2に加えて、材料則36も用いることができます。しかしながら、要素は長さ方向には1積分点のみ持つため、フレーム構造の1つの線に1要素を用いることは、深さ方向だけでなく長さ方向の塑性の進展を考慮するために、推奨されません。


図 2. 積分ビームの断面定義