With the intensifying exploitation of the resources and the waterway of the Arctic region in worldwide, the variable inclination propulsion system are gradually used in the ice-sailing ship, its design and strength study are of great engineering significance and application value. A variable inclination propulsion shafting with universal joint is introduced and its mechanical model is established in this paper. By combine torsional stress and inclination angle calculation, a decoupling algorithm based on Newmark is designed, and then the analysis of ice loads response in both frequency domain and time domain are completed. The results show that due to the complex triangle function relation among the various components of the shafting, there is an interaction between the propeller and the transmission shafting. In particular, when considering the dynamic change of inclination, nonlinear bending-torsional coupling vibration may occur and further lead to the enhancement effect of torsional vibration. Among the three working conditions defined by DNV, the bending-torsional coupling effect vary with rotation speed and inclination angle, and is mainly concentrated in high frequency bands. In the ice loads condition 2, with the largest time domain excitation amplitude and the smallest high order component, it’s vibration response is the strongest, and the change law of dangerous phase angle is similar to the bending moment. With the increase of initial inclination angle, the amplitude of low frequency vibration decreases, yet the amplitude of high frequency vibration increase. This trend is more obvious with the increase of rotation speed, and more high frequency resonance points appear in the condition 1 and condition 3.