Elastoplastic analysis by return-mapping method with rigorous loading criterion for extended subloading surface model