In this work a model, based on a X-ray diffraction contrast tomography data of a stainless steel wire with a diameter of 0.4mm is presented. As measured 3D grain geometry and crystallographic orientation of individual grains are directly transferred into a finite element model. Anisotropic elasticity and crystal plasticity constitutive laws are used for the bulk grain material while the grain boundaries are explicitly modeled using the cohesive zone approach. A parametric study on the effects of the grain boundary strength and other cohesive zone parameters on the macroscopic response and damaged grain boundary area of a polycrystalline aggregate is presented. Recommendations for the cohesive zone parameters values aimed at achieving low damaged grain boundary area during numerical tensile tests are given while at the same time taking into account the numerical stability of the simulations.