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 fxed parameters. Layer thickness, infll 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 identifed layer thickness as the major controlled variable for all responses. Interaction of infll percentage and extruder temperature had a signifcant infuence on elongation at break and therefore, tough fracture of printed parts. The input parameters were optimized to materialize tow criteria; the frst 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 efciency of the optimization method. There was a good agreement between empirical results and response surface method predictions which confrmed the reliability of predictive models. The optimal setting to fulfll the frst criterion could bring on a specimen with more than 1500 (N) maximum failure load and less than 9 (g) weight.