“…Experimental approaches have been widely used for such studies [1][2][3][4][5][6][7][8][9], while numerical simulations of increasing sophistication are enhancing our understanding of the often complex mechanical response. Recently, three-dimensional (3D) finite element (FE) models have been used to capture the non-uniform, through-thickness stress distributions [10][11][12][13][14][15][16][17], and models to predict the onset and growth of material damage have been incorporated [18][19][20][21][22][23][24][25][26]. The majority of these analyses used implicit FE solvers, which generally have difficulties in converging all the way to the ultimate failure point of CBJs [27], mainly due to the complexity of contact definition between bolts and a damaging laminate.…”