In this paper, a finite element model is developed for 3000 series clad aluminum alloy brazing sheet to account for temperature and strain rate dependency, as well as plastic anisotropy. The current work considers a novel implementation of the Barlat YLD2000 yield surface in conjunction with the Bergstrom hardening model to accurately model aluminum alloy sheet during warm forming. The Barlat YLD2000 yield criterion is used to capture the anisotropy while the Bergstrom hardening rule predicts the temperature and strain rate dependency. The results are compared with those obtained from experiments. The measured stress–strain curves of the AA3003 aluminum alloy sheet at elevated temperatures and different strain rates are used to fit the Bergstrom parameters and measured R-values and directional yield stresses are used to fit the yield function parameters. Isothermal uniaxial tensile tests and nonisothermal deep drawing experiments are performed and the predicted response using the new constitutive model is compared with measured data. In simulations of tensile tests, the material behavior is predicted accurately by the numerical models. Also, the nonisothermal deep drawing simulations are able to predict the load–displacement response and strain distributions accurately.