Close this search box.
Clear all

progressive damage in VUMAT subroutine

2 Posts
2 Users
0 Reactions
Posts: 2
Topic starter
New Member
Joined: 14 years ago

I am making a code of vumat for progressive damage, the Isotropic Hardening Plasticity is used.
I dont know how to deal with the material point with damage, especially when damage is equal to 1. I do it like that following:

if (damg .ge. one) then
stressNew(k,1) = zero
stressNew(k,6) = zero
e=(1-damg) *props(1)
xnu = props(2)
twomu = e / ( one + xnu )
alamda = xnu * twomu / ( one - two * xnu )
thremu = op5 * twomu
eps=(1-damg)*eps Get the damaged yield stress!
stateNew(k,5)=... Update the damg variable!
end if

Am i correct?
Thanks in advance.

Topic Tags
1 Reply
Posts: 2
New Member
Joined: 13 years ago

dear t_y2004,

im actually doing the similar subroutine as you, and what ive done is also based the same idea as you posted here.
but unfortunately the result is not all good.
ive got very nice responce for one sigle element.
but for a real sample for example tensile test sample (5000 elements), the response is bad and the deformation of some damaged elements are totally distorted.
have you got the samilar problem?


Topic Tags
1 Reply