In a nuclear industry area, high temperature treatment of materials is a factor which requires special attention. Assessment needs to be conducted on the properties of the materials used, including the strength of the materials. The measurement of material properties under thermal processes may reflect residual stresses. The use of Genetic Algorithm (GA) to determine the optimal residual stress is one way to determine the strength of a material. In residual stress modeling with several parameters, it is sometimes difficult to solve for the optimal value through analytical or numerical calculations. Here, GA is an efficient algorithm which can generate the optimal values, both minima and maxima. The purposes of this research are to obtain the optimization of variable in residual stress models using GA and to predict the center of residual stress distribution, using fuzzy neural network (FNN) while the artificial neural network (ANN) used for modeling. In this work a single-material 316/316L stainless steel bar is modeled. The minimal residual stresses of the material at high temperatures were obtained with GA and analytical calculations. At a temperature of 650 0 C, the GA optimal residual stress estimation converged at–711.3689 MPa at adistance of 0.002934 mm from center point, whereas the analytical calculation result at that temperature and position is-975.556 MPa. At a temperature of 850 0 C, the GA result was-969.868 MPa at 0.002757 mm from the center point, while with analytical result was-1061.13 MPa. The difference in residual stress between GA and analytical results at a temperatureof650 0 C is about 27%, while at 850 0 C it is 8.67%. The distribution of residual stress showed a grouping concentrated around a coordinate of (-76; 76) MPa. The residuals stress model is a degree-two polynomial with coefficients of 50.33,-76.54, and-55.2, respectively, with a standard deviation of 7.874.