Frozen fractured rock mass is often encountered during the implementation of geotechnical engineering in cold regions. The tensile strength parameters of frozen rock play an important role in the construction of rock slopes involving tensile failure. In order to study the tensile characteristics of a frozen fractured rock mass in a cold region, original rock specimens were mined and processed in the Yulong Copper Mine, and artificial, frozen fractured marble specimens were made. The effects of different ice-filled crack angles, lengths, and widths on the force–displacement curve and the tensile strength of frozen rock were studied by laboratory Brazilian splitting experiments and RFPA^3D, and the evolution law of the tensile strength of frozen rock was revealed. At the same time, wing crack initiation and cracking mode after tensile failure were analyzed by high-speed camera; the whole process of the Brazilian splitting of frozen rock was reconstructed, and the development of microcrack initiation in frozen rock was analyzed. The following conclusions were drawn from the test results: the frozen rock specimens have typical brittle-failure characteristics. The tensile strength of frozen rock decreases gradually with the increase in the width and length of ice-filled cracks, and decreases first and then increases with the increase in the angle of the ice-filled crack. The ice-filled crack incurs damage first, and then the wing cracks start from the tip of the ice-filled crack and extend continuously. The tensile strength of frozen rock is significantly affected by the angle and length of ice-filled cracks.

Multilayered geological structures are common in geotechnical engineering, where cooling shrinkage induces polygonal cracks in interlayers, compromising rock mass strength, permeability, and long-term stability. Existing thermo-mechanical studies on layered rock cracking insufficiently address the combined effects of weak interlayer geometry or interface-regulated mechanisms. To address this gap, based on meso-damage mechanics and thermodynamics, this study adopts a thermo-mechanical coupling simulation method considering rock heterogeneity, innovatively focusing on the understudied stress transfer effect at strong–weak interlayer interfaces. Systematic investigations were conducted on the initiation, propagation, and saturation of polygonal cracks in plate-like layered rocks under surface cooling, analyzing the influences of weak interlayer thickness, number, and position, and comparing surface vs. inner interlayer behaviors. Results showed that stress transfer interruption at weak–strong layer interfaces can inhibit crack propagation. Inter weak interlayers produce significantly more cracks and fragments than surface weak interlayers, with a stratified crack length distribution, and the maximum fragment area increases exponentially with weak interlayer thickness. Crack development is strongly influenced by weak interlayer thickness, with thinner layers dominated by non-penetrating cracks and thicker layers tending to develop penetrating lattice-like cracks. The inter layer stress and crack distribution exhibit fractal characteristics, with crack density decreasing layer by layer and no new cracks forming after saturation. This study clarifies the regulatory mechanism of weak interlayer features and surface cooling on crack evolution, quantifies interface effects and fractal characteristics, and provides a theoretical basis for instability prediction of layered rock structures in low-temperature geotechnical engineering.