To analyze the nonflatness and residual stress in accelerated cooled plate, a numerical analysis model has been developed. Two factors, i.e. temperature and phase transformation, are considered in calculating the stress distribution that develops during cooling. The plastic strain and plate-buckling, which are often shown in accelerated cooled plate, were determined from this stress. Mean temperature in through thickness direction and temperature difference in width direction are considered in the model to simplify the calculation. The temperature and stress distribution changes caused by phase transformation are involved in terms of the effective specific heat and the effective thermal expansion coefficient. With the model, accelerated cooling of 10mm(t) $^{*}$3000mm(w) plate was simulated. The condition of accelerated cooling was .deg. C/sec from just after hot rolling to 500.deg. C. The initial temperature-difference ratio, .DELTA.Tr, in width direction is an important factor in evaluating the stress distribution. When .DELTA.Tr is 0.08, buckling occurs during cooling and 7kgf/m $m^{2}$ of residual stress develops at the edge of plate. To secure the flatness, .DELTA.Tr should be less than 0.07. Small scaled cooling test was conducted to verify the exactness of the model and the results proved the usefulness of this numerical analysis model.l.