Thank you dear @Yogi Firmansyah for your replay, it is the only solution but i have problem how can i calculate tagent mudulus for the case of bilinear kinematic hardening behavior (is it Et= ((Su-Sy)/(eu-ei)) ? ) if i have this input http://www.matweb.com/search/datasheet_print.aspx?matguid=abc4415b0f8b490387e3c922237098da&src=xls