Graphene can at present be grown at large quantities only by the chemical vapor deposition method, which produces polycrystalline samples. Here, we describe a method for constructing realistic polycrystalline graphene samples for atomistic simulations, and apply it for studying their mechanical properties. We show that cracks initiate at points where grain boundaries meet and then propagate through grains predominantly in zigzag or armchair directions, in agreement with recent experimental work. Contrary to earlier theoretical predictions, we observe normally distributed intrinsic strength (∼ 50% of that of the mono-crystalline graphene) and failure strain which do not depend on the misorientation angles between the grains. Extrapolating for grain sizes above 15 nm results in a failure strain of ∼ 0.09 and a Young's modulus of ∼ 600 GPa. The decreased strength can be adequately explained with a conventional continuum model when the grain boundary meeting points are identified as Griffith cracks.Grain boundaries define the electronic and mechanical properties of polycrystalline materials. Since chemical vapor deposition (CVD) is currently the only way for producing industry-scale graphene membranes, and leads to polycrystalline samples, study of grain boundaries in graphene has become of fundamental importance during the recent years. In a two-dimensional material, such as graphene, the boundaries also have a critical contribution to the chemical reactivity. Because of this, although atomic scale imaging can in principle reveal their exact structure, the boundaries tend to be covered by adsorbates with only short segments available for direct imaging. Nevertheless, experiments [1-6] have revealed meandering serpent-like boundaries which are typically formed from pentagon-heptagon-pairs in the parts not covered by the adsorbates.Mechanical properties of graphene sheets have been a topic of intense research already for two decades in the context of carbon nanotubes (see Ref.[7] for a topical review). More recently, in 2007 [8], Liu and co-workers utilized ab initio calculations to study the elastic moduli and fracture characteristic of graphene. Young's modulus was found to be 1.05 TPa, and failure strain, depending on the pulling direction, 0.194-0.266. Intrinsic strength was estimated to be 110-121 GPa, similarly depending on the pulling direction. The role of pre-existing defects on these properties was also studied [9]. It was noticed that their effect does not depend on the exact atomic structure of the defects but rather on their size. The authors also showed that the intrinsic strength of graphene with crack-like defects can be described with a continuum model using the Griffith formula for defect sizes down to 10Å. Soon after this, Frank et al. used a tip of an atomic force microscope to obtain a Young's modulus of 0.5 TPa for suspended stacks of graphene sheets [10]. A year later, Lee and co-workers reported on several mechanical properties of graphene using a similar technique [11], establishing graphene ...