In the present study additive manufacturing of Polylactic acid by fused deposition modeling were investigated based on statistical analysis. The honeycomb internal pattern was employed to build inside of specimens due to its remarkable capability to resist mechanical loads. Simplify 3D was utilized to slice the 3D model and to adjust fixed parameters. Layer thickness, infill percentage, and extruder temperature were considered as controlled variables, while maximum failure load (N), elongation at break (mm), part weight (g), and build time (min) were selected as output responses and analysed by response surface method. Analysis of variance results identified layer thickness as the major controlled variable for all responses. Interaction of infill percentage and extruder temperature had a significant influence on elongation at break and therefore, tough fracture of printed parts. The input parameters were optimized to materialize tow criteria; the first one was to rise maximum failure load and the second was to attain tough fracture and lessen build time and part weight at a time. Optimal solutions were examined by experimental fabrication to evaluate the efficiency of the optimization method. There was a good agreement between empirical results and response surface method predictions which confirmed the reliability of predictive models. The optimal setting to fulfill the first criterion could bring on a specimen with more than 1500 (N) maximum failure load and less than 9 (g) weight.