Ceramic-polymer composite materials are used throughout materials science, but theoretical models which take into account the microstructural design to predict their mechanical properties are still not fully developed. We present an approach to use the discrete element method to model the mechanical behavior under bending load of dense composite materials made from ceramic particles which are bonded together by polymeric layers. Unlike many other modeling approaches, the internal particulate structure of the material, including the particle size distribution, packing structure, and pore structure can thus be considered. Linear-elastic bonds are created between all contacting particles to model the polymeric binder. A three-dimensional beam with a packing density of 63% is generated and its mechanical properties are tested in 3- and 4-point bending. The loading speed, loading scheme (position of the supports), and the mechanical properties of the polymeric bonds are varied and their effect on the modulus of elasticity of the material is investigated. (C) 2013 The Authors. Published by Elsevier B.V. All rights reserved.