We use discrete element model simulations to model the full boudinage process from initial fracturing of intact material to post-fracture flow of material into gaps between fragments and to investigate the role which the material properties of the weak and strong layers play in this process. The models are deformed in 2D plane strain under a range of confining stresses, in coaxial bulk flow. Results show that the material properties, i.e. Mohr-Coulomb or quasi-viscous for the matrix and elastic-brittle for the competent layer, lead to the development of natural looking boudin morphologies and deformation patterns in the matrix. The details of the matrix rheology only have a minor influence on the morphology of the boudins. By varying the material properties of the competent layer between fully brittle and semi-ductile we obtain a wide range of deformation patterns ranging from pinch-and-swell structures to a variety of boudin types including drawn, shear band and straight sided torn boudins. In a number of models we observe rotation of the boudin blocks despite the applied deformation being purely coaxial. These rotations are generally related to asymmetrical (rhombic) boudin shapes. Some features observed in natural boudins such as concave block faces or the formation of veins between fragments are not modeled because pore fluids are not yet included in our model.