Kinematic Joint TYPE33 (PROP/KJOINT)

Kinematic joints are declared by /PROP/KJOINT. Joints are defined by a spring and two local coordinate axes, which belong to connected bodies.

Assume that the connected bodies are rigid to ensure the orthogonality of their local axis. However, deformable bodies may also be connected with a joint. If the axis becomes non-orthogonal during deformation, the stability of the joint cannot be ensured.

There are several kinds of kinematic joints available in Radioss, which are listed in Kinematic Joint Types.

Kinematic Joint Types

Type No. Joint Type dx dy dz θ x θ y θ z
1 Spherical x x x 0 0 0
2 Revolute x x x 0 x x
3 Cylindrical 0 x x 0 x x
4 Planar x 0 0 0 x x
5 Universal x x x x 0 0
6 Translational 0 x x x x x
7 Oldham x 0 0 x x x
8 Rigid x x x x x x
9 Free 0 0 0 0 0 0

x: denotes a blocked degree of freedom

0: denotes a free (user-defined) degrees of freedom

Joint properties are defined in a local frame computed with respect to two connected coordinate systems. They do not need to be initially coincident. If the initial position of the local coordinate axis coincides at any time, the joint local frames are defined at a mean position. Then the joint local frame will be computed with respect to these rotated axes.

There are a total of six joint degrees of freedom: δ X , δ Y , δ Z , θ X , θ Y MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH0oazce WGybGbauaacaGGSaGaeqiTdqMabmywayaafaGaaiilaiabes7aKjqa dQfagaqbaiaacYcacqaH4oqCceWGybGbauaacaGGSaGaeqiUdeNabm ywayaafaaaaa@460B@ and θ Z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCce WGAbGbauaaaaa@3900@ . They are computed in the local skew frame.


Figure 1. Kinematic Joint Definition

In each type of joint you distinguish the blocked degrees of freedom and the free degrees of freedom. The blocked degrees of freedom are characterized by a constant stiffness. Selecting a high value with respect to the free degrees of freedom stiffness is recommended. The free degrees of freedom have user-defined characteristics, which can be linear or nonlinear elastic, combined with a sub-critical viscous damping.

The translational and rotational degrees of freedom are defined as:(1)
δ=d x 2 d x 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH0oazcq GH9aqpcaWGKbGaamiEamaaBaaaleaacaaIYaaabeaakiabgkHiTiaa dsgacaWG4bWaaSbaaSqaaiaaigdaaeqaaaaa@3F9C@
Where, d x 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGKbGaam iEamaaBaaaleaacaaIXaaabeaaaaa@392C@ and d x 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGKbGaam iEamaaBaaaleaacaaIXaaabeaaaaa@392C@ are total displacement of two joint nodes in the local coordinate system.(2)
θ = θ 2 θ 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCcq GH9aqpcqaH4oqCdaWgaaWcbaGaaGOmaaqabaGccqGHsislcqaH4oqC daWgaaWcbaGaaGymaaqabaaaaa@3F4D@
Where, θ 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCda WgaaWcbaGaaGymaaqabaaaaa@38FC@ and θ 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCda WgaaWcbaGaaGymaaqabaaaaa@38FC@ are total relative rotation of two connected body axes, with respect to the local joint coordinate frame.


Figure 2. Schematic Representation of Kinematic Joints

Forces and Moments Calculation

The force in direction δ is computed as:
  • Linear spring: 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@ : Translational stiffness (Ktx, Kty, and Ktz)
  • C t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHdbWaaS baaSqaaiaadshaaeqaaaaa@3850@ : Translational viscosity (Ctx, Cty, and Ctz)
  • Nonlinear spring: F= K t f( δ )+ C t g( δ ˙ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHgbGaey ypa0JaaC4samaaBaaaleaacaWG0baabeaakiGacAgadaqadaqaaiaa hs7aaiaawIcacaGLPaaacqGHRaWkcaWHdbWaaSbaaSqaaiaadshaae qaaOGaci4zamaabmaabaGabCiTdyaacaaacaGLOaGaayzkaaaaaa@4488@
The moment in direction θ is computed as:
  • Linear spring: 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@ : Rotational stiffness (Krx, Kry, and Krz)
  • C r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHlbWaaS baaSqaaiaadshaaeqaaaaa@3858@ : Rotational viscosity (Crx, Cry, and Crz)
  • Nonlinear spring: M= K r f( θ )+ C r g( θ ˙ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHnbGaey ypa0JaaC4samaaBaaaleaacaWGYbaabeaakiGacAgadaqadaqaaiaa hI7aaiaawIcacaGLPaaacqGHRaWkcaWHdbWaaSbaaSqaaiaadkhaae qaaOGaci4zamaabmaabaGabCiUdyaacaaacaGLOaGaayzkaaaaaa@4493@
The joint length may be equal to 0. It is recommended to use a zero length spring to define a spherical joint or a universal joint. To satisfy the global balance of moments in a general case, correction terms in the rotational degrees of freedom are calculated as:(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

Joints do not have user-defined mass or inertia, so the nodal time step is always used.