/RWALL

ブロックフォーマットのキーワード 無限平面、無限円筒、球形、および平行四辺形の剛壁を定義します。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/RWALL/type/rwall_ID/unit_ID
rwall_title
node_ID Slide grnd_ID1 grnd_ID2            
Dsearch fric ϕ ffac ifq  
node_ID = 0の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
XM YM ZM        
node_ID0の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Mass VX0 VY0 VZ0    
typeがPlane、Cyl、Paralの場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
XM1 YM1 ZM1        
typeがParalの場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
XM2 YM2 ZM2        

定義

フィールド 内容 SI単位の例
type 剛壁タイプのキーワード

(下の表をご参照ください)

 
rwall_ID 剛壁の識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
rwall_title 剛壁のタイトル

(文字、最大100文字)

 
node_ID 節点識別子(移動剛壁)

(整数)

 
Slide スライディングフラグ。
= 0
スライディング
= 1
タイド
= 2
摩擦を伴うスライディング

(整数)

 
grnd_ID1 剛壁に追加するセカンダリ節点を定義する節点グループ

(整数)

 
grnd_ID2 剛壁から削除するセカンダリ節点を定義する節点グループ

(整数)

 
Dsearch セカンダリ検索の距離

(実数)

[ m ]
fric 摩擦

(実数)

 
ϕ 球の直径

(実数)

[ m ]
ffac フィルタリング係数

デフォルト値はifqフラグに依存します (実数)

 
ifq フィルタリングフラグ 5
= 0(デフォルト)
フィルターを使用しません。
= 1
f f a c = α MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaadA gacaWGHbGaam4yaiabg2da9iabeg7aHbaa@3C40@ 単純なフィルター値は0~1です。
= 2
ffac =N、フィルターを適用する時間ステップ。
= 3
ffac=freqフィルターを適用するカットオフ周波数

デフォルト = 0(整数)

 
XM MのX座標

(実数)

 
YM MのY座標

(実数)

 
ZM MのZ座標

(実数)

 
Mass 剛壁の質量 8

質量が入力されていない場合、剛壁は一定の強制速度を有するようになります。

(実数)

[ kg ]
VX0 X方向の初速度

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
VY0 Y方向の初速度

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
VZ0 Z方向の初速度

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
XM1 M1のX座標

(実数)

 
YM1 M1のY座標

(実数)

 
ZM1 M1のZ座標

(実数)

 
XM2 M2のX座標

(実数)

 
YM2 M2のY座標

(実数)

 
ZM2 M2のZ座標

(実数)

 

剛壁タイプ

タイプ
内容
PLANE
¯ 平面
CYL
¯ 右記の直径の円柱; ϕ
SPHER
右記の直径の球; ϕ
PARAL
平行四辺形

サーフェス入力タイプ

タイプ
内容
PLANE
MM1は法線方向を定義します。

rwall_lagmul_plane2
CYL
MM1は円柱の軸を定義します。
rwall_cyl
SPHER
Mは球の中心です。
rwall_spher
PARAL
M M 1 M M 2 は平行四辺形を定義
rwall_paral

コメント

  1. 剛壁を定義するために最初に入力するのは、ポイントMの座標です。移動剛壁の場合は、node_IDを入力します。
  2. 次に入力するのはポイントM1の座標で、ポイントM2も入力できます(移動壁の場合、M1およびM2のモーションはnode_IDと同じになります)。
  3. 剛壁のセカンダリ節点は、節点のグループや初期段階で剛壁からDsearchほど下にある節点として定義できます。
  4. 摩擦フィルタリングオプションは、摩擦を伴ってスライドする剛壁に対してのみ使用できます(Slide=2)。
  5. フィルタリングフラグifg

    ifq0の場合、接触部内の各セカンダリ節点における接線摩擦力は、次の単純なルールを使用してフィルタリングされます:

    (1)
    F T f = α F T ( t ) + ( 1 α ) F T f ( t d t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaGccqGH9aqpcqaHXoqycaWHgbWaaSba aSqaaiaadsfaaeqaaOGaaiikaiaadshacaGGPaGaey4kaSYaaeWaae aacaaIXaGaeyOeI0IaeqySdegacaGLOaGaayzkaaGaaCOramaaBaaa leaacaWGubGaamOzaaqabaGccaGGOaGaamiDaiabgkHiTiaadsgaca WG0bGaaiykaaaa@4D30@
    ここで、
    F T f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaaaaa@38B5@
    フィルタリングされた接線力。
    F T ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubaabeaakiaacIcacaWG0bGaaiykaaaa@3A26@
    フィルター前の時間tにおける計算された接線力。
    F T f ( t d t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaGccaGGOaGaamiDaiabgkHiTiaadsga caWG0bGaaiykaaaa@3DE0@
    前の時間ステップでフィルタリングされた接線力
    t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaads haaaa@37D8@
    現在のシミュレーション時間
    d t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaads haaaa@37D8@
    現在のシミュレーション時間ステップ
    α MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3796@
    フィルタリング係数

    フラグifqは、フィルタリング手法である 係数を定義します。

    ifq =1の場合、フィルタリング係数はユーザーが直接入力します: α = f f a c

    ifq =2の場合、はユーザー定義周波数(タイムステップ数を単位として定義された周波数)の3dBフィルタリングレベルに対応します:(2)
    α = 2 π N

    ここで、 1 f r e q = T = N d t N = f f a c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOtaiabg2 da9iaadAgacaWGMbGaamyyaiaadogaaaa@3B73@

    ifq =3の場合、はユーザー定義周波数の3dBフィルタリングレベルに対応します:(3)
    α = 2 π d t f r e q MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8WjY=viVeYth9vqqj=hEeuD0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqySdeMaeyypa0JaaGOmaiabec8aWjaadsgacaWG0bGaeyyXICTa amOzaiaadkhacaWGLbGaamyCaaaa@46D6@
    ここで、
    d t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaads haaaa@37D8@
    時間ステップ
    f r e q MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaadk hacaWGLbGaamyCaaaa@39B8@
    ffac
  6. 球(SPHER)と平行四辺形(PARAL)は無限ではありません。平行四辺形では、以下を使用して法線を定義します:(4)
    M M 1 × M M 2

    と球の直径は、 ϕ を使用して定義します。

  7. 剛壁セカンダリ節点の節点板厚は考慮されません。
  8. MASS=0または空白の移動剛壁の場合、剛壁は一定の強制速度を有し、初期速度はありません。