The dead water resistance F1=12CdwρSU2 (ρ fluid density, U ship speed, S wetted body surface, Cdw resistance coefficient) on a ship moving at subcritical speed along the upper layer of a two-layer fluid is calculated by a strongly nonlinear method assuming potential flow in each layer. The ship dimensions correspond to those of the Polar ship Fram. The ship draught, b0, is varied in the range 0.25h0–0.9h0 (h0 the upper layer depth). The calculations show that Cdw/(b0/h0)2 depends on the Froude number only, in the range close to critical speed, Fr = U/c0 ∼ 0.875–1.125 (c0 the linear internal long wave speed), irrespective of the ship draught. The function Cdw/(b0/h0)2 attains a maximum at subcritical Froude number depending on the draught. Maximum Cdw/(b0/h0)2 becomes 0.15 for Fr = 0.76, b0/h0 = 0.9, and 0.16 for Fr = 0.74, b0/h0 = 1, where the latter extrapolated value of the dead water resistance coefficient is about 60 times higher than the frictional drag coefficient and relevant for the historical dead water observations. The nonlinear Cdw significantly exceeds linear theory (Fr < 0.85). The ship generated waves have a wave height comparable to the upper layer depth. Calculations of three-dimensional wave patterns at critical speed compare well to available laboratory experiments. Upstream solitary waves are generated in a wave tank of finite width, when the layer depths differ, causing an oscillation of the force. In a wide ocean, a very wide wave system develops at critical speed. The force approaches a constant value for increasing time.