Industrial systems, such as logistics and supply chain networks, are complex systems because they comprise a big number of interconnected actors and significant nonlinear and stochastic features. This paper analyzes a distribution network design problem for a four-echelon supply chain. The problem is represented as an inventory-location model with uncertain demand and a continuous review inventory policy. The decision variables include location at the intermediate levels and product flows between echelons. The related safety and cyclic inventory levels can be computed from these decision variables. The problem is formulated as a mixed integer nonlinear programming model to find the optimal design of the distribution network. A linearization of the nonlinear model based on a piecewise linear approximation is proposed. The objective function and nonlinear constraints are reformulated as linear formulations, transforming the original nonlinear problem into a mixed integer linear programming model. The proposed approach was tested in 50 instances to compare the nonlinear and linear formulations. The results prove that the proposed linearization outperforms the nonlinear formulation achieving convergence to a better local optimum with shorter computational time. This method provides flexibility to the decision-maker allowing the analysis of scenarios in a shorter time.