Equação de Movimento do Navio: Manobrabilidade e Referenciais

A Equação de Movimento do Navio

  • O objetivo principal de construir a equação de movimento no estudo da manobrabilidade do navio é conhecer as causas dos movimentos.

  • A fundação dessas equações reside na Segunda Lei de Newton, onde as forças externas que atuam sobre o corpo são relacionadas à aceleração resultante.

Sistemas Referenciais de Posição

  • O navio se move frequentemente em seis graus de liberdade (66 DOF): três translações em eixos ortogonais e três rotações em torno desses eixos.

  • As forças externas dependem significativamente da geometria do casco. Em um referencial inercial, as funções que descrevem essa geometria variariam constantemente com o tempo devido ao movimento do casco.

  • Simetria do Casco: Praticamente todos os navios e veículos hidrodinâmicos possuem um plano de simetria (plano da linha central). Boreste e bombordo são espelhados geometricamente. Essa simetria é observada em navios, submarinos, foguetes, aviões, peixes e pássaros.

  • Assimetria Dinâmica: Pode ocorrer devido ao sentido de rotação de um único propulsor, mas é um efeito controlado.

  • Escolha do Sistema de Eixos: Tira proveito da simetria. Dois eixos ficam no plano de simetria e o terceiro é perpendicular a ele. Em corpos de revolução (torpedos, foguetes), existe um segundo plano de simetria perpendicular ao primeiro.

Definição dos Eixos no Referencial do Navio
  • Eixo x: Eixo longitudinal no plano de simetria, positivo para frente (vante). Geralmente paralelo à quilha ou linha d'água. Vetor unitário: i\mathbf{i}.

  • Eixo y: Eixo transversal, perpendicular ao plano de simetria, positivo para boreste. Vetor unitário: j\mathbf{j}.

  • Eixo z: Eixo vertical, no plano de simetria, positivo para baixo (em direção à quilha). Vetor unitário: k\mathbf{k}.

Referenciais Inerciais e o Sistema NED
  • Referencial Inercial: Referencial onde as leis de Newton são válidas. Quando o navio rotaciona, os eixos fixos a ele deixam de ser inerciais, complicando o cálculo do lado direito da equação de movimento (variação da quantidade de movimento).

  • Sistema North-East-Down (NED): Sistema inercial geográfico utilizado para análise de navegação, trajetórias, ângulo de rumo e forças ambientais (vento, corrente, ondas).

    • Eixo N (ou x): Norte.

    • Eixo E (ou y): Leste (East).

    • Eixo D (ou z): Down (para baixo).

O Vetor Posição e Velocidade do Navio

Vetor Posição
  • O vetor posição do navio no sistema inercial (RbR_b) define a localização do navio no tempo.

  • Para um elemento infinitesimal de massa dmdm, o vetor posição R0R_0 no sistema NED é dado pela soma:   R0=Rb+RGR_0 = R_b + R_G   Onde RbR_b acompanha o centro de gravidade no referencial inercial e RGR_G acompanha a posição de dmdm em relação ao centro de gravidade.

Variação dos Vetores Unitários com a Rotação
  • Quando o navio rotaciona, os vetores unitários i,j,k\mathbf{i}, \mathbf{j}, \mathbf{k} mudam de sentido, logo suas derivadas no tempo não são zero:

    • Rotação θ\theta (Pitch) no eixo y: didt=kdθdt\frac{d\mathbf{i}}{dt} = -\mathbf{k} \frac{d\theta}{dt}; djdt=0\frac{d\mathbf{j}}{dt} = 0; dkdt=idθdt\frac{d\mathbf{k}}{dt} = \mathbf{i} \frac{d\theta}{dt}

    • Rotação ψ\psi (Yaw) no eixo z: didt=jdψdt\frac{d\mathbf{i}}{dt} = \mathbf{j} \frac{d\psi}{dt}; djdt=idψdt\frac{d\mathbf{j}}{dt} = -\mathbf{i} \frac{d\psi}{dt}; dkdt=0\frac{d\mathbf{k}}{dt} = 0

    • Rotação ϕ\phi (Roll) no eixo x: didt=0\frac{d\mathbf{i}}{dt} = 0; djdt=kdϕdt\frac{d\mathbf{j}}{dt} = \mathbf{k} \frac{d\phi}{dt}; dkdt=jdϕdt\frac{d\mathbf{k}}{dt} = -\mathbf{j} \frac{d\phi}{dt}

Vetor Velocidade Angular
  • Definido por Ω=pi+qj+rk\vec{\Omega} = p\mathbf{i} + q\mathbf{j} + r\mathbf{k}, onde:

    • p=dϕdtp = \frac{d\phi}{dt} (velocidade angular de roll).

    • q=dθdtq = \frac{d\theta}{dt} (velocidade angular de pitch).

    • r=dψdtr = \frac{d\psi}{dt} (velocidade angular de yaw).

  • Relação cinemática fundamental: didt=Ω×i\frac{d\mathbf{i}}{dt} = \vec{\Omega} \times \mathbf{i}, djdt=Ω×j\frac{d\mathbf{j}}{dt} = \vec{\Omega} \times \mathbf{j}, dkdt=Ω×k\frac{d\mathbf{k}}{dt} = \vec{\Omega} \times \mathbf{k}.

Vetor Velocidade Linear
  • O vetor velocidade do navio VGV_G (ou vGv_G) medido no sistema solidário ao navio é:   VG=ui+vj+wkV_G = u\mathbf{i} + v\mathbf{j} + w\mathbf{k}

    • uu: Velocidade longitudinal (vante).

    • vv: Velocidade transversal (boreste).

    • ww: Velocidade vertical (baixo).

Equação de Movimento no Referencial Inercial

  • O navio é tratado como um corpo rígido (tamanho e forma não variam).

  • Quantidade de movimento (PP): P=dmvNEDP = \sum dm \cdot v_{NED}.

  • Em translação pura, o corpo é tratado como uma grande partícula de massa Δ\Delta (deslocamento) concentrada no CG.

  • Segunda Lei de Newton: F=dPdt=ΔdVNEDdtF = \frac{dP}{dt} = \Delta \frac{dV_{NED}}{dt}.

  • Em navios mercantes, a taxa de variação de massa devido ao combustível (dΔdt\frac{d\Delta}{dt}) é negligenciada durante a manobra (dΔdt=0\frac{d\Delta}{dt} = 0).

Equações de Movimento Translacionais (Origem no CG)

Considerando a massa Δ\Delta constante e substituindo as derivadas dos vetores unitários, as forças nas direções xx, yy e zz (referencial do navio) são:

  • Força em x (Surge):   X=Δ(u˙+wqvr)X = \Delta (\dot{u} + wq - vr)

  • Força em y (Sway):   Y=Δ(v˙+urwp)Y = \Delta (\dot{v} + ur - wp)

  • Força em z (Heave):   Z=Δ(w˙+vpuq)Z = \Delta (\dot{w} + vp - uq)

Nota: Os termos u˙,v˙,w˙\dot{u}, \dot{v}, \dot{w} são acelerações translacionais, enquanto (wqvr),(urwp),(vpuq)(wq - vr), (ur - wp), (vp - uq) são acelerações centrípetas.

Equações de Movimento para Momentos (Origem no CG)

  • A quantidade de movimento angular (HGH_G) é dada pelo tensor de inércia.

  • Como os eixos são os eixos principais de inércia, os produtos de inércia são nulos.

  • Momento Angular: HG=Ixpi+Iyqj+IzrkH_G = I_x p\mathbf{i} + I_y q\mathbf{j} + I_z r\mathbf{k}.

  • Variação Temporal dos Momentos (Equações de Euler):

    • Momento em x (Roll - K):     K=Ixp˙+(IzIy)qrK = I_x \dot{p} + (I_z - I_y)qr

    • Momento em y (Pitch - M):     M=Iyq˙+(IxIz)prM = I_y \dot{q} + (I_x - I_z)pr

    • Momento em z (Yaw - N):     N=Izr˙+(IyIx)pqN = I_z \dot{r} + (I_y - I_x)pq

Os termos (IzIy)qr(I_z - I_y)qr, etc., representam efeitos giroscópicos (precessão e nutação).

Equações com Origem fora do Centro de Gravidade

Frequentemente, a origem OO é fixada no centro geométrico do navio para facilitar o cálculo das forças hidrodinâmicas, mesmo que o centro de gravidade GG esteja em outra posição (xG,yG,zGx_G, y_G, z_G).

Forças Translacionais em Origem Não Inercial (OGO \neq G)

Novos termos aparecem devido às forças de Euler e centrífugas:

  • X:   X=Δ[u˙vr+wqxG(q2+r2)+yG(pqr˙)+zG(pr+q˙)]X = \Delta [\dot{u} - vr + wq - x_G(q^2 + r^2) + y_G(pq - \dot{r}) + z_G(pr + \dot{q})]

  • Y:   Y=Δ[v˙wp+uryG(r2+p2)+zG(qrp˙)+xG(qp+r˙)]Y = \Delta [\dot{v} - wp + ur - y_G(r^2 + p^2) + z_G(qr - \dot{p}) + x_G(qp + \dot{r})]

  • Z:   Z=Δ[w˙uq+vpzG(p2+q2)+xG(rpq˙)+yG(rq+p˙)]Z = \Delta [\dot{w} - uq + vp - z_G(p^2 + q^2) + x_G(rp - \dot{q}) + y_G(rq + \dot{p})]

Momentos em Origem Não Inercial (OGO \neq G)

Os momentos K,M,NK, M, N em relação à origem OO incluem termos da aceleração de Coriolis e reações inerciais:

  • K (Roll):   K=Ixp˙+Ixy(prq˙)Ixz(r˙+pq)+(IzIy)qr+Iyz(r2q2)+Δ[yG(w˙+vpuq)zG(v˙+urwp)]K = I_x \dot{p} + I_{xy}(pr - \dot{q}) - I_{xz}(\dot{r} + pq) + (I_z - I_y)qr + I_{yz}(r^2 - q^2) + \Delta [y_G(\dot{w} + vp - uq) - z_G(\dot{v} + ur - wp)]

  • M (Pitch):   M=Iyq˙+Iyz(qpr˙)Ixy(p˙+qr)+(IxIz)pr+Ixz(p2r2)+Δ[zG(u˙+wqvr)xG(w˙+vpuq)]M = I_y \dot{q} + I_{yz}(qp - \dot{r}) - I_{xy}(\dot{p} + qr) + (I_x - I_z)pr + I_{xz}(p^2 - r^2) + \Delta [z_G(\dot{u} + wq - vr) - x_G(\dot{w} + vp - uq)]

  • N (Yaw):   N=Izr˙+Ixz(qrp˙)Iyz(q˙+rp)+(IyIx)pq+Ixy(q2p2)+Δ[xG(v˙+urwp)yG(u˙+wqvr)]N = I_z \dot{r} + I_{xz}(qr - \dot{p}) - I_{yz}(\dot{q} + rp) + (I_y - I_x)pq + I_{xy}(q^2 - p^2) + \Delta [x_G(\dot{v} + ur - wp) - y_G(\dot{u} + wq - vr)]

Simplificações Comuns: Plano Horizontal

Referencial fora do CG (3 DOF em Surge, Sway, Yaw)
  • X: X=Δ[u˙vrxGr2+yG(pqr˙)+zG(pr+q˙)]X = \Delta [\dot{u} - vr - x_G r^2 + y_G(pq - \dot{r}) + z_G(pr + \dot{q})] (Simplificado dependendo dos planos de interesse).

Referencial no CG (3 DOF em Surge, Sway, Yaw)

As equações reduzem-se drasticamente ao analisar apenas o movimento no plano horizontal (p=q=w=0p=q=w=0):

  • Surge: X=Δ(u˙vr)X = \Delta (\dot{u} - vr)

  • Sway: Y=Δ(v˙+ur)Y = \Delta (\dot{v} + ur)

  • Yaw: N=Izr˙N = I_z \dot{r}

Resumo de Definições Técnicas (Convenção SNAME)

  • Surge (Avanço): Força XX, aceleração u˙\dot{u}.

  • Sway (Deriva): Força YY, aceleração v˙\dot{v}.

  • Heave (Arfagem): Força ZZ, aceleração w˙\dot{w}.

  • Roll (Balanço): Momento KK, aceleração angular p˙\dot{p}.

  • Pitch (Arfagem): Momento MM, aceleração angular q˙\dot{q}.

  • Yaw (Guinada): Momento NN, aceleração angular r˙\dot{r}.