運動学的ジョイント TYPE33(/PROP/KJOINT)

運動学的ジョイントは/PROP/KJOINTで宣言されます。ジョイントは1つのスプリングと結合された物体に属する2つの局所座標軸によって定義されます。

ここで局所軸の直交性を保証するため、結合されている物体は剛体であることを仮定します。しかしながら、変形可能な物体もジョイントに結合することができます。変形の間に軸の直交性が崩れた場合、ジョイントの安定性は保証できません。

Radiossではいくつもの種類のジョイントが使用可能です。運動学的ジョイントタイプにそれらをリストします。

運動学的ジョイントタイプ

タイプ番号 ジョイントタイプ dx dy dz θ x θ y θ z
1 x x x 0 0 0
2 回転 x x x 0 x x
3 円筒 0 x x 0 x x
4 平面 x 0 0 0 x x
5 ユニバーサル x x x x 0 0
6 並進 0 x x x x x
7 オールダム x 0 0 x x x
8 剛性 x x x x x x
9 フリー 0 0 0 0 0 0

x: 固定自由度を示します。

0: フリー自由度(ユーザー定義)を示します。

ジョイントのプロパティは2つの結合された座標系に関して計算された局所座標系で定義されます。これらは最初の時点で一致している必要はありません。局所座標軸の初期位置が常に一致する場合は、ジョイントの局所フレームは平均位置に定義されます。次にその局所座標系はこれらの回転軸に関して計算されます。

6つのジョイント自由度があります: δ X , δ Y , δ Z , θ X , θ Y MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH0oazce WGybGbauaacaGGSaGaeqiTdqMabmywayaafaGaaiilaiabes7aKjqa dQfagaqbaiaacYcacqaH4oqCceWGybGbauaacaGGSaGaeqiUdeNabm ywayaafaaaaa@460B@ および θ Z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCce WGAbGbauaaaaa@3900@ の4つのパラメータのみを指定して使用できます。これらは局所スキューフレームで計算されます。


図 1. 運動学的ジョイントの定義

各ジョイントタイプでは、固定自由度とフリー自由度を区別します。固定自由度には剛性が一定しているという特徴があります。フリー自由度の剛性には大きい値を選択することを推奨します。フリー自由度では、ユーザー定義の特徴(線形または非線形弾性の場合がある)に臨界前粘性減衰を組み合わせます。

並進および回転自由度は次のように定義されます:(1)
δ=d x 2 d x 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH0oazcq GH9aqpcaWGKbGaamiEamaaBaaaleaacaaIYaaabeaakiabgkHiTiaa dsgacaWG4bWaaSbaaSqaaiaaigdaaeqaaaaa@3F9C@
ここで、 d x 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGKbGaam iEamaaBaaaleaacaaIXaaabeaaaaa@392C@ d x 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGKbGaam iEamaaBaaaleaacaaIXaaabeaaaaa@392C@ は局所座標系における2つのジョイント節点の全変位です。(2)
θ = θ 2 θ 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCcq GH9aqpcqaH4oqCdaWgaaWcbaGaaGOmaaqabaGccqGHsislcqaH4oqC daWgaaWcbaGaaGymaaqabaaaaa@3F4D@
ここで θ 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCda WgaaWcbaGaaGymaaqabaaaaa@38FC@ θ 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCda WgaaWcbaGaaGymaaqabaaaaa@38FC@ は局所ジョイント座標フレームに関する2つの結合された物体の全相対回転です。


図 2. 運動学的ジョイントの概要図

力とモーメントの計算

方向 δ の力は次のように計算されます:
  • 線形スプリング: F= K t δ+ C t δ ˙ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHgbGaey ypa0JaaC4samaaBaaaleaacaWG0baabeaakiaahs7acqGHRaWkcaWH dbWaaSbaaSqaaiaadshaaeqaaOGabCiTdyaacaaaaa@3F9D@
  • K t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHlbWaaS baaSqaaiaadshaaeqaaaaa@3858@ :並進剛性(KtxKtyおよびKtz
  • C t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHdbWaaS baaSqaaiaadshaaeqaaaaa@3850@ :並進粘度(CtxCtyおよびCtz
  • 非線形スプリング: F= K t f( δ )+ C t g( δ ˙ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHgbGaey ypa0JaaC4samaaBaaaleaacaWG0baabeaakiGacAgadaqadaqaaiaa hs7aaiaawIcacaGLPaaacqGHRaWkcaWHdbWaaSbaaSqaaiaadshaae qaaOGaci4zamaabmaabaGabCiTdyaacaaacaGLOaGaayzkaaaaaa@4488@
方向 θ のモーメントは次のように計算されます:
  • 線形スプリング: M= K r θ+ C r θ ˙ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHnbGaey ypa0JaaC4samaaBaaaleaacaWGYbaabeaakiaahI7acqGHRaWkcaWH dbWaaSbaaSqaaiaadkhaaeqaaOGabCiUdyaacaaaaa@3FA8@
  • K r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHlbWaaS baaSqaaiaadshaaeqaaaaa@3858@ :回転剛性(KrxKryKrz
  • C r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHlbWaaS baaSqaaiaadshaaeqaaaaa@3858@ :回転粘度(CrxCryCrz
  • 非線形スプリング: M= K r f( θ )+ C r g( θ ˙ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHnbGaey ypa0JaaC4samaaBaaaleaacaWGYbaabeaakiGacAgadaqadaqaaiaa hI7aaiaawIcacaGLPaaacqGHRaWkcaWHdbWaaSbaaSqaaiaadkhaae qaaOGaci4zamaabmaabaGabCiUdyaacaaacaGLOaGaayzkaaaaaa@4493@
ジョイントの長さは0にすることもできます。球面ジョイントやユニバーサルジョイントの定義には長さ0のスプリングを用いることを推奨します。一般的なケースでのモーメントの全体での釣り合いを満足するため、回転自由度の補正項が次のように計算されます:(3)
M θ x = M θ x + L y × F z L z × F y
(4)
M θ y = M θ y + L z × F x L x × F z
(5)
M θ z = M θ z + L x × F y L y × F x

ジョイントはユーザー定義の質量や慣性をもつことができません、このため節点時間ステップが常に用いられます。