I am writing UEXPAN subroutine in ABAQUS 2017 to calculate thermal strain increment in my domain for powder bed electron beam additive manufacturing. This is how my subroutine looks:
subroutine uexpan(expan,dexpandt,temp,time,dtime,predef,dpred,
1 statev,cmname,nstatv,noel)
c...
I am simulating a 2-dimensional additive manufacturing process where there is a moving Gaussian heat source on top of a powder layer. The heating step is followed by a cooling step after which a new powder layer is added and the same heat source moves on top of this new layer. The image shown...