Regarding the potential damage and failure issues that may occur in high-temperature transient environments during the service process of ceramic matrix composite laminates, the present work proposes a multi-layer thermomechanical induced crack initiation model for ceramic matrix composite laminates subjected to thermal gradients and transients. This model can determine the history of temperature, deformation, and stress distribution within each layer of the material, as well as the steady-state energy release rate of all possible crack locations. Then the influence of bending constraints on material stress distribution and energy release rate is investigated, finding that maximum stress and energy release rate values are significantly lower without bending constraints compared to those under bending constraints. Furthermore, the study explores the impact of material structural parameters and boundary conditions on the energy release rate of ceramic matrix composite laminates, providing direct insights for designing laboratory tests and evaluating the lifespan of these materials under in-service conditions.