Keywords: TDC, Submarine, Torpedo, World War II

Transcription

1 Self Position Keeper Modeling Abstract During World War II, all the submarines of all warring nations used analog computers to assist submarine commanders with aiming straight-running torpedoes (e.g. UK used their "fruit machine."). US submarines used a computer called the Torpedo Data Computer (TDC) that was considered the best of these devices.the Torpedo Data Computer (TDC) included a function known as Position Keeping that was a mechanical solver for a system of differential equations.this technical note compares the output from a simple Mathcad model of the Position Keeper with the output from a kinematic model for a simple test case. I solve the system of ODEs using Mathcad's standard ODE solver and a homebrew routine that may be more appropriate for people who are implementing software versions of the TDC's position keeper function. Table of Contents Update TOC Abstract Table of Contents TOC Introduction Analysis of Bearing Angle and ange Versus Time Simulation Approach Field of Battle Kinematic Model Differential Equation Model Euler Method Solution Graphical Display Standard Differential Equation Solver Conclusion eference Introduction Jump to region The model I am using here comes from following web site. See reference section for excerpt I am going to focus on the Position Keeper function of the TDC, which is made up of two coupled differential equations and two subsidiary equations. The key function of the Position Keeper is to provide estimates of the target's bearing and range. My simulation will focus on these variables, which are critical to obtaining a fire control solution, which I will not cover here. A more thorough modeling would include more complex scenarios, but this illustrates the basic approach used. Tactically, the submarine would estimate the target course, range, and speed from periscope and hydrophone readings. Position Keeper Model for the TDC 1 of 6 22-February-212

2 Self The TDC would put out a constantly updated position for the target, which the sub skipper would compare against further periscope readings. The TDC target model would be updated until it could accurately predict the target's motion. With the TDC, the submarine could maneuver and the TDC would keep track of the target's position relative to the submarine. This greatly improved the accuracy of the whole fire control operation. Analysis of Bearing Angle and ange Versus Time Simulation Approach I am going to work this problem assuming the simple case of a target and submarine (referred to as "Own Ship") both pursuing constant velocity courses. This will allow me to determine the exact distance and bearing numbers using a simple kinematic model and the differential equation model. This way, I can verify that my differential equation solution is reasonable. I could use the Mathcad differential equation solvers, but I decided to put together a simple solution using Euler's method (i.e. the simplest possible way). This would be easy for someone to code using any number of programming languages. Field of Battle Assume that we are going to place our submarine and target on a 1 meter by 1 meter grid indexed as shown in Figure 1. 1 meters Y Axis meters meters X Axis 1 meters Figure 1: Torpedo Data Computer Simulation Grid Position Keeper Model for the TDC 2 of 6 22-February-212

3 Self Kinematic Model I am going to work this problem assuming the simple case of a target and submarine both pursuing constant velocity courses. This will allow me to determine the exact distance and bearing numbers using a simple kinematic model and the differential equation model. Target Own Ship 1 x coordinate 9 T m 9 O m y coordinate 1 x coordinate y coordinate 1 m x component 1 m V T V 1 s O y component 5 s x component y component T () t T V T t O () t O V O t t () T () t O () t Target Position as a function of time. Submarine Position (i.e. Own Ship) Position as a function of time. Distance between target and submarine. T t θ( t) 18deg atan2 T () t O () t Differential Equation Model () O () t 1 Bearing of target from submarine. S O V O m s S T V T m s Speed of the Submarine (Own Ship) Speed of the Target B' ( BrA ) S O sin( Br) ' ( BrA) S O cos( Br) S T sin( A) S T cos( A) Bearing Angle Differential Equation ange Differential Equation C O atan2 V O V O1 C T atan2 V T V T1 18deg 18deg deg 225deg Own Ship Course (sub knows its course) Target Ship Course (sub estimates this from periscope reading of angle on the bow) Position Keeper Model for the TDC 3 of 6 22-February-212

4 Self Euler-Method Solution B Br A dt.1s B θ( ) 1259m It is a bit crude, but it appears to work. Algorithm time step Initialize Bearing Initialize ange A B 18deg C T Initialize Target Angle Br B C O Initialize elative Bearing for i 1 9 Main Loop B dt B' Br A i i1 i1 i1 B i 1 Predict the Bearing Br B C i i O A B 18deg C i i T dt ' Br A i i i i1 Compute the Next elative Bearing Compute the Next Target Angle Compute Next ange B dt B' Br A B Corrector for Bearing i i i i 1 Graphical Display B Br m A eturn all the values I needed to remove the units from the range vector (a Mathcad 15 limitation) Here is a graph comparing my kinematic model with the model from the differential equation solution. They are identical, which they should be. The nice thing about the differential equation solution is that it can handle changes in course by the submarine. These course changes were automatically fed into the TDC by the sub's gyrocompass. Position Keeper Model for the TDC 4 of 6 22-February-212

5 Self i Each time increment corresponds to.1 seconds. Kinematic Model vs Simulated TDC Output 25 ange (meters) Bearing Angle ( ) Time Increment (.1 seconds) Exact Target ange (Kinematic) TDC Target ange Exact Target Bearing (Kinematic) TDC Target Bearing Standard Differential Equation Solver Mathcad has excellent ODE solvers. Here is the same problem worked using one of their ODE solver routine. T1 9 Maximum time -- the solver does not like units S O S T S O S m T m The solver does not like units s s Given d du Bu ( ) d du u ( ) S O sin B( u) C O S T sin B( u) π C T = B ( ) =.785 u ( ) = S O cos B( u) C O S T cos B( u) π C T ( ) = 1259 f g B Odesolve u T1 9 z1.1 9 Position Keeper Model for the TDC 5 of 6 22-February-212

6 Self ange (m) Solution using the Mathcad Standard Solver Bearing Angle ( ) ange Bearing Time (sec) My crude solver got the same result. Conclusion This model appears to provide a reasonable example for the operation of the TDC, at least as I read it in the old manual. eference Old Navy TDC manual excerpt Position Keeper Model for the TDC 6 of 6 22-February-212