Numerical modeling of failure modes in bolted composite laminates