The　numerical　manifold　method　(NMM),　combined　with　the　precise　time　integration　method　(PTIM),　is　proposed　for　thermal　shock　fracture　analysis.　The　temperature　and　displacement　discontinuity　across　crack　faces　is　naturally　portrayed　attributing　to　the　cover　systems　in　the　NMM.　The　crack　tip　singularities　are　characterized　through　the　use　of　asymptotic　bases　in　the　approximations.　The　discrete　equations　for　transient　thermal　analysis　are　firstly　solved　with　the　PTIM　and　then　the　thermoelastic　study　is　performed.　With　the　interaction　integral,　the　stress　intensity　factors　are　computed.　Several　examples　are　tested　and　the　nice　consistency　between　the　present　and　existing　results　is　found.