Rigid Body

A rigid body is defined by a main node and its associated secondary nodes. Mass and inertia may be added to the initial main node location. The main node is then moved to the center of mass, taking into account the main node and all secondary node masses. Figure 1 shows an idealized rigid body.


Figure 1. Idealized Rigid Body

Rigid Body Mass

The mass of the rigid body is calculated by:(1)
m=mM+ImI
The rigid body's center of mass is defined by:(2)
xG=mMxM+mIxIm
(3)
yG=mMyM+mIyIm
(4)
zG=mMzM+mIzIm
Where,
mM
Main node mass
mI
Secondary node masses
xG , yG , zG
Coordinates of the mass center

Rigid Body Inertia

The six components of inertia of a rigid body are computed by:(5)
Ixx=JxxM+mM((yMyG)2+(zMzG)2)+i(Ixxi+mi((yiyG)2+(zizG)2)) MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG4bGaamiEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadIhacaWG4baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamyEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaGccqGHRaWkdaqadaqaaiaadQhadaWgaaWcbaGaamytaaqabaGccqGHsislcaWG6bWaaSbaaSqaaiaadEeaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaIYaaaaaGccaGLOaGaayzkaaGaey4kaSYaaabuaeaadaqadaqaaiaadMeadaqhaaWcbaGaamiEaiaadIhaaeaacaWGPbaaaOGaey4kaSIaamyBamaaCaaaleqabaGaamyAaaaakmaabmaabaWaaeWaaeaacaWG5bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamyEamaaBaaaleaacaWGhbaabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaaGOmaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaGaaGOmaaaaaaaakiaawIcacaGLPaaaaiaawIcacaGLPaaaaSqaaiaadMgaaeqaniabggHiLdaaaa@7007@
(6)
Iyy=JyyM+mM((xMxG)2+(zMzG)2)+i(Iyyi+mi((xixG)2+(zizG)2)) MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG5bGaamyEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadMhacaWG5baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaGccqGHRaWkdaqadaqaaiaadQhadaWgaaWcbaGaamytaaqabaGccqGHsislcaWG6bWaaSbaaSqaaiaadEeaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaIYaaaaaGccaGLOaGaayzkaaGaey4kaSYaaabuaeaadaqadaqaaiaadMeadaqhaaWcbaGaamyEaiaadMhaaeaacaWGPbaaaOGaey4kaSIaamyBamaaCaaaleqabaGaamyAaaaakmaabmaabaWaaeWaaeaacaWG4bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamiEamaaBaaaleaacaWGhbaabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaaGOmaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaGaaGOmaaaaaaaakiaawIcacaGLPaaaaiaawIcacaGLPaaaaSqaaiaadMgaaeqaniabggHiLdaaaa@7009@
(7)
Izz=JzzM+mM((xMxG)2+(yMyG)2)+i(Izzi+mi((xixG)2+(yiyG)2)) MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG6bGaamOEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadQhacaWG6baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaGccqGHRaWkdaqadaqaaiaadMhadaWgaaWcbaGaamytaaqabaGccqGHsislcaWG5bWaaSbaaSqaaiaadEeaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaIYaaaaaGccaGLOaGaayzkaaGaey4kaSYaaabuaeaadaqadaqaaiaadMeadaqhaaWcbaGaamOEaiaadQhaaeaacaWGPbaaaOGaey4kaSIaamyBamaaCaaaleqabaGaamyAaaaakmaabmaabaWaaeWaaeaacaWG4bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamiEamaaBaaaleaacaWGhbaabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaaGOmaaaakiabgUcaRmaabmaabaGaamyEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaGaaGOmaaaaaaaakiaawIcacaGLPaaaaiaawIcacaGLPaaaaSqaaiaadMgaaeqaniabggHiLdaaaa@700C@
(8)
Ixy=JxyM+mM((xMxG)+(yMyG))+i(Ixyimi((xixG)+(yiyG))) MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG4bGaamyEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadIhacaWG5baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamyEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaaaOGaayjkaiaawMcaaiabgUcaRmaaqafabaWaaeWaaeaacaWGjbWaa0baaSqaaiaadIhacaWG5baabaGaamyAaaaakiabgkHiTiaad2gadaahaaWcbeqaaiaadMgaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamyEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaaaaaaaaOGaayjkaiaawMcaaaGaayjkaiaawMcaaaWcbaGaamyAaaqab0GaeyyeIuoaaaa@6D1E@
(9)
Iyz=JyzM+mM((yMyG)+(zMzG))+i(Iyzimi((yiyG)+(zizG))) MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG5bGaamOEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadMhacaWG6baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamyEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaaaOGaayjkaiaawMcaaiabgUcaRmaaqafabaWaaeWaaeaacaWGjbWaa0baaSqaaiaadMhacaWG6baabaGaamyAaaaakiabgkHiTiaad2gadaahaaWcbeqaaiaadMgaaaGcdaqadaqaamaabmaabaGaamyEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadMhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaaaaaaaaOGaayjkaiaawMcaaaGaayjkaiaawMcaaaWcbaGaamyAaaqab0GaeyyeIuoaaaa@6D2C@
(10)
Ixz=JxzM+mM((xMxG)+(zMzG))+i(Ixzimi((xixG)+(zizG))) MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBaaaleaacaWG4bGaamOEaaqabaGccqGH9aqpcaWGkbWaa0baaSqaaiaadIhacaWG6baabaGaamytaaaakiabgUcaRiaad2gadaahaaWcbeqaaiaad2eaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGnbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaaaOGaayjkaiaawMcaaiabgUcaRmaaqafabaWaaeWaaeaacaWGjbWaa0baaSqaaiaadIhacaWG6baabaGaamyAaaaakiabgkHiTiaad2gadaahaaWcbeqaaiaadMgaaaGcdaqadaqaamaabmaabaGaamiEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadIhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaaaakiabgUcaRmaabmaabaGaamOEamaaBaaaleaacaWGPbaabeaakiabgkHiTiaadQhadaWgaaWcbaGaam4raaqabaaakiaawIcacaGLPaaadaahaaWcbeqaamaaCaaameqabaaaaaaaaOGaayjkaiaawMcaaaGaayjkaiaawMcaaaWcbaGaamyAaaqab0GaeyyeIuoaaaa@6D25@
Where,
Iij
Moment of rotational inertia in the ij MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyAaiaadQgaaaa@37D3@ direction
JijM
Main node added inertia

Rigid Body Force And Moment Computation

The forces and moments acting on the rigid body are calculated by:(11)
F=FM+iFi
(12)
M=MM+iMi+iSiG×Fi
Where,
FM
Force vector at the main node
Fi
Force vector at the secondary nodes
MM
Moment vector at the main node
Mi
Moment vector at the secondary nodes
G
Vector from secondary node to the center of mass

Resolving these into orthogonal components, the linear and rotational acceleration may be computed as:

Linear Acceleration(13)
γi=Fim
Rotational Acceleration(14)
I1α1=M1(I3I2)ω2ω3
(15)
I2α2=M2(I1I3)ω1ω3
(16)
I3α3=M3(I2I1)ω1ω2
Where,
Ii
Principal moments of inertia of the rigid body
α1
Rotational accelerations in the principal inertia frame (reference frame)
ωi
Rotational velocity in the principal inertia frame (reference frame)
Mi
Moments in the principal inertia frame (reference frame)

Time Integration

Time integration is performed to find velocities of the rigid body at the main node:(17)
ν(t+Δt2)=ν(tΔt2)+γ(t)Δt
(18)
ω(t+Δt2)=ω(tΔt2)+α(t)Δt

Where, v is the linear velocity vector. Rotational velocities are computed in the local reference frame.

The velocities of secondary nodes are computed by:(19)
νi=νM+SiGxω
(20)
ωi=ωM

Boundary Conditions

The boundary conditions given to secondary nodes are ignored. The rigid body has the boundary conditions given to the main node only.

A kinematic condition is applied on each secondary node, for all directions. A secondary node is not allowed to have any other kinematic conditions.

No kinematic condition is applied on the main node. However, the rotational velocities are computed in a local reference frame. This reference frame is not compatible with all options imposing rotation such as imposed velocity, rotational, rigid link.

The only exception concerns the rotational boundary conditions for which a special treatment is carried out. Connecting shell, beam or spring with rotation stiffness to the main node, is not yet allowed either.