A servo-controlled tension-internal pressure testing machine with an optical 3D deformation analysis system (ARAMIS , GOM) was used to measure the multiaxial plastic deformation behavior of a high-strength steel sheet with a tensile strength of 590 MPa for a strain range from initial yield to fracture. Tubular specimens were fabricated by roller bending and laser welding the as-received flat sheet materials. Many linear stress paths in the first quadrant of the stress space were applied to the tubular specimens to measure the forming limit curve (FLC), forming limit stress curve (FLSC), and forming limit plastic work per unit volume (FLPW) of the as-received sheet material in addition to the contours of plastic work and the directions of the plastic strain rates. Differential hardening behavior was observed; the shapes of the work contours constructed in the principal stress space changed with an increase in plastic work. The observed differential hardening behavior was approximated by changing the material parameters and the exponent of the Yld2000-2d yield function as functions of the reference plastic strain. Marciniak-Kuczyński-type forming limit analyses were performed using both the differential hardening model and isotropic hardening models based on the Yld2000-2d yield function. It was found that the material model that is capable of reproducing both the work contours and the directions of the plastic strain rates measured for a strain range close to the fracture limit can give a more effective constitutive model for accurately predicting the FLC, FLSC, and FLPW.T. Hakoyama JSPS Research Fellow (Doctoral Course Students),