A methods to calculate non-linear moment-curvature relations of high-strength PSC flexural members for numerical analysis has been proposed. The moment-curvature relations were calculated with assumptions of design codes and by the layer method. The results of the proposed procedures for moment-curvature relations and numerical analysis were compared with those of pre-existing tests. The absorption energy rate of the design codes was about 30% larger than that of the layer method. The ultimate load and the external work of the layer method were 90% and 85% of those of tests, respectively The ultimate load of the strength design method was 97% of that of tests, but the external work was over-estimated with 122%. The ultimate load and external work by the proposed equation of the CEB-FIP Model Code were 113% and 173% of those of tests, respectively. It show that the use of ultimate strain of 0.0035 should be over-estimated for high-strength concrete. The procedure of non-linear numerical analysis of this research could be stably simulated the behavior of concrete flexural members until the ultimate state, and calculate results of the load-deflection relation and cracking pattern were very similar with those of tests.