Use Vumat to simulate the hot ductile fracture criterion of Freudenthal

Dear everybody，
I have just studied Abaqus and would like to use Vumat to simulate the hot ductile

fracture criterion of Freudenthal which have the equation: D=integral(0-fracture strain)

of (Smises stress)*d(strain). The fracture will be appear when I reach 1 unit.
Could you please show me how to write algorithm adding temperature to calculate (D) and

out put (D) in Vumat.

thanks you so much

I am looking forward to hearing from you

hongfengtl