Setting the intial value of a STATEV

[fusuym]

Hi all,

I'm using a the user subroutine HETVAL to simulate epoxy curing.

I have a STATEV(1) = DEGREE OF CURE which is updated in the subroutine.

All I am wondering is how you set this value to zero at the very start of the solution.

I feel it is a really simple thing but am a bit lost!

Thanks

Yusuf

 

[fusuym]

In fact,

is it that you can use subroutine SDVINI before HETVAL?

SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,
1 LAYER,KSPT)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION STATEV(NSTATV),COORDS(NCRDS)


user coding to define STATEV(NSTATV)


RETURN
END

 

[andraina]

Hi Yusuf,

This subroutine enables the determination of the heat produced by the thermosetting reaction, for each Gauss point, thanks to the degree of cure and degree of cure rate provided by the FIELD array.

Cf. “A FEM coupling model for properties prediction during the curing of an epoxy matrix” (N. Rabearison, C. Jochum, JC. Grandidier) – Computational Material Science45 (2009) 715-724

See you