軸方向力が作用する破壊抵抗曲げモーメントの算出について

求められた解答ではないのですが，一言。
破壊時は普通コンクリートの圧縮破壊で決まるので，圧縮側のひずみをコンクリートの終局ひずみとします。
さて，この状態で中立軸の位置を圧縮縁からxとすれば，部材内部のひずみ分布は直線で求まりますね。断面での各位置でのひずみがわかれば，鉄筋やコンクリートの応力度も求まり，これに断面積をかければ圧縮力や引張力がxの式として出てきます。
そうしておいて，全圧縮力＝全引張力＋軸方向圧縮力とすれば，xが求まるはずです。
xが求まれば，あとは圧縮力や引張力の値が確定しますから，中立軸からの距離を乗じて全部を足せば，曲げ破壊モーメントが求まるということになります。
と書くのは簡単なのですが，実際やってみると確か3次方程式（だったかな）になってしまいますので，求めるのは結構厄介です。こういう事情で具体的な計算例は見たことがありません。軸方向圧縮力がなければ2次方程式ですむはずですが。
ま，力業で解くか，EXCELなどの表計算で求めてみるしかないと思います。