A three-dimensional meshfree method for modeling arbitrary crack initiation and crack growth in reinforced concrete structure is presented. This meshfree method is based on a partition of unity concept and formulated for geometrically non-linear problems. The crack kinematics are obtained by enriching the solution space in order to capture the correct crack kinematics. A cohesive zone model is used after crack initiation. The reinforcement modeled by truss or beam elements is connected by a bond model to the concrete. We applied the method to model the fracture of several reinforced concrete structures and compared the results to experimental data.