abaqus材料用户子程序vumat用于梁单元时，
应力如何更新，子程序用B31梁单元做。