Type91 Aerial Torpedo Vectors of Motion Equations.
This drowing was for the set of aerial torpedo motion equations in the air, once taught in the class of Rear Admiral Shoji Naruse, who was the chief designer of Type 91 aerial torpedo, since late 1930s till 1945. The equations and parameters are shown as followings.The torpedo motion equation is the set of simultaneous ordinary differential equations, which is to model pitch motion of airborne aerial torpedo as follows. Clear mistake errors in equations have been corrected, which the original private book Koku Gyorai Note had in typescript symbol letters in privately printed edition service.Eq.1: falling velocity of the torpedo equation
Eq.2: horizontal vector velocity of the torpedo mass equation
Eq.3: vertical vector velocity of the torpedo equation
Eq.4: vertical vector acceleration of the torpedo mass equation
Eq.5: angular velocity equation with respect to time
Eq.6: angular velocity differential equation with respect to timedx/dt=VX⋯(Eq.1)W/g×(dVX/dt)=−Dcos⁡φ−Lsin⁡φ⋯(Eq.2)dz/dt=VZ…(Eq.3)W/g×(dVZ/dt)=Dsin⁡φ−Lcos⁡φ+W⋯(Eq.4)dθ/dt=ω⋯(Eq.5)I×(dω/dt)=57.3M−bVω⋯(Eq.6){\displaystyle {\begin{array}{lcr}dx/dt&amp;=&amp;V_{X}\cdots (Eq.1)\\W/g\times \left(dV_{X}/dt\right)&amp;=&amp;-D\cos \varphi -L\sin \varphi \cdots (Eq.2)\\dz/dt&amp;=&amp;V_{Z}\ldots (Eq.3)\\W/g\times \left(dV_{Z}/dt\right)&amp;=&amp;D\sin \varphi -L\cos \varphi +W\cdots (Eq.4)\\d\theta /dt&amp;=&amp;\omega \cdots (Eq.5)\\I\times \left(d\omega /dt\right)&amp;=&amp;57.3M-bV\omega \cdots (Eq.6)\end{array}}}where constant 57.3 of Eq.6 is the coefficient from 1 (radian) = 57.2958 °b V ω of Eq.6 is the damping moment, where b is defined as;b=−(δCmgH/δαH)×(ρSl/2)×lH{\displaystyle b=-(\delta C_{mgH}/\delta \alpha _{H})\times (\rho S{\mathit {l}}/2)\times {\mathit {l}}_{H}\,}Since the angular moment of a torpedo here is the lifting movement as follows;ωlH/V(radian)=(δCmgH/δαH)×(ρSl/2)×lH×V×ω{\displaystyle \omega {\mathit {l}}_{H}/V\,(radian)\,=(\delta C_{mgH}/\delta \alpha _{H})\times (\rho S{\mathit {l}}/2)\times {\mathit {l}}_{H}\times V\times \omega \,}V : the velocity of the torpedo
VH : the horizontal axis velocity of the torpedo
VZ : the vertical axis velocity of the torpedo
φ : the moving vector angle of the torpedo in reference to horizontal axis
θ : the posture angle of the torpedo in reference to the horizontal axis
α : the internal angle between φ and θ, which is equal to the lift angle of tail fins of the torpedo
G : the gravity center of the torpedo
W : weight of the torpedo
g : acceleration of gravity, 9.8 m/sec2
I : inertia coefficient with respect to lift moment at the gravity center of the torpedo
ω : lifting angular velocity
D : drag moment of force
L : lift moment of force
M : roll moment of force around the longitudial axis of the torpedo
ρ : air density
S : cross-section area of the torpedo
l : total length of the torpedo
lH : the length between the gravity center and the center of lift moment of tail fins of the torpedo
ωlH : the angular velocity at the center of lift force of tail fins of the torpedoThe value lH is measured in wind-tunnel test,
as drag moment of force coefficient differences between the torpedoes with and without box type attachment of wooden tail stabilizer plates;lH=(−CmgHl+CtHd)/CnH{\displaystyle {\mathit {l}}_{H}=(-C_{mgH}\,{\mathit {l}}+C_{tHd})/C_{nH}\,}where each coefficient is defined as;CnH = CXH sin α + CZH cos α 
CtH = CXH cos α + CZH sin α CX : drag moment of force coefficient D / (1/2 ρ V2 S )
CZ : lift moment of force coefficient L / (1/2 ρ V2 S)
Cmg : roll moment of force coefficient around the gravity center of the torpedo M / (1/2 ρ V2 l )
CXH : drag moment of force coefficient of tail stabilizer box of torpedo
CZH : lift moment of force coefficient of tail stabilizer box of torpedo
CmgH : roll moment of force coefficient around the gravity center of tail stabilizer box of the torpedofor further details, see Type 91 torpedo or Koku Gyorai Note
Own work by uploader, and converted my PNG file to SVG file
Shun Zero