Fiber metal laminates, that is made of two aluminum sheets in exterior layers and a self-reinforced polypropylene in interior layer, have favorable characteristics such as excellent fatigue resistance, high impact absorbability, possibility to reform and recycle components, reduction of manufacture processing time, and so on. For these advantages, fiber metal laminates have been widely researched and applied to aerospace and automotive industries. To manufacture commercial products with fiber metal laminates, one is required to examine the formability of fiber metal laminates using numerical analysis. In order to analyze the formability numerically, the flow stress of fiber metal laminates should be identified experimentally or analytically by constructing the flow stress model. In case of the fiber metal laminates containing thermoplastic composite, however, temperature-dependent flow stress modeling should be required analytically to avoid inefficient experiments. In this study, the uniaxial tensile tests were carried out to obtain the observed flow stresses at various temperature conditions. The representative flow stress models such as Hollomon, Ludwik, and Johnson-Cook models were constructed by regression fitting techniques. Next, the modified Hollomon and Ludwik models were proposed while the material constants based on the Hollomon and Ludwik models were fitted by linear regression fitting. Moreover, the newly modified Ludwik models were proposed to improve the accuracy to predict the flow stress while the material constants based on the Ludwik model were fitted by quadratic and exponential equations with respected to the temperature. Then, the flow stress models were validated graphically as well as numerically. In the graphical validation, the flow stresses predicted by the flow stress models were compared with the observed flow sresses obtained by the uniaxial tensile tests. In the numerical validation, average maximum absolute error and R-squared error were calculated to compare the errors of the flow stress models quantitatively. It is found that the flow stress model in which the material constants were formulated as only quadratic equations, together with the model where the material constants were formulated as a combination of exponential (yield stress) and quadratic equations (strength coefficient and work hardening exponent) with respect to the temperature, predict the flow stress more appropriately than the other flow stress models.