说明：  abaqus二次开发vumat子程序，适用于三维梁单元的弹塑性本构模型，采用Kent-Park受压应力应变全程曲线
(abaqus secondary development vumat routine for three-dimensional beam element constitutive model, using the Kent-Park compressive stress-strain curve of the whole)