Search
Close this search box.
Notifications
Clear all

Definition of tickness strain in VUMAT

2 Posts
2 Users
0 Reactions
744 Views
Posts: 17
Topic starter
(@thegiving21)
Active Member
Joined: 8 years ago

Dear All,

I have implemented a damage model into user subroutine VUMAT, but I dont know how to define the tickness strain in the case of plane stress. Can anyone help me int this regard ?

subroutine vumat(
1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
2 stepTime, totalTime, dt, cmname, coordMp, charLength,
3 props, density, strainInc, relSpinInc,
4 tempOld, stretchOld, defgradOld, fieldOld,
3 stressOld, stateOld, enerInternOld, enerInelasOld,
6 tempNew, stretchNew, defgradNew, fieldNew,
5 stressNew, stateNew, enerInternNew, enerInelasNew )

include vaba_param.inc
C integer ndir, nshr, nstatev, nfieldv, nprops, lanneal
C real stepTime, totalTime, dt
C

dimension props(nprops), density(nblock),
1 coordMp(nblock,*),
2 charLength(*), strainInc(nblock,ndir+nshr),
3 relSpinInc(*), tempOld(*),C(ndir+nshr,ndir+nshr),
4 stretchOld(*), defgradOld(*),epsO(nblock,ndir+nshr),
5 fieldOld(*), stressOld(nblock,ndir+nshr),eigVal(nblock,3),
6 stateOld(nblock,nstatev), enerInternOld(nblock),
7 enerInelasOld(nblock), tempNew(*), ,
8 stretchNew(*), defgradNew(*), fieldNew(*),
9 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev),
1 enerInternNew(nblock), enerInelasNew(nblock),
2 epsN(nblock,ndir+nshr), stressinc(nblock,ndir+nshr)
character*80 cmname
*USER SUBROUTINE

PARAMETER (ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,
+ SIX=6.D0)
C

E=props(1)
xnu=props(2)

C
A=ONE
C ----- plane stress
DDSDE(1,1)=E/(ONE-xnu*xnu)
DDSDE(1,2)=(E*xnu)/(ONE-xnu*xnu)
DDSDE(1,3)=ZERO
DDSDE(2,1)=(E*xnu)/(ONE-xnu*xnu)
DDSDE(2,2)=E/(ONE-xnu*xnu)
DDSDE(2,3)=ZERO
DDSDE(3,1)=ZERO
DDSDE(3,2)=ZERO
DDSDE(3,3)=(E/(ONE-xnu*xnu))*(ONE-xnu)

do i = 1, nblock

C ---- definition of the tickness strain
.
.
C--- main code
.
.
.

enddo
return
end

Topic Tags
1 Reply
Posts: 80
(@FrankMonkey)
Trusted Member
Joined: 15 years ago

Subscribe to and seek assistance from the
ABAQUS mailing list
[url] https://groups.yahoo.com/neo/groups/Abaqus/info [/url]
Search the archive of the list before posting in it.
The list does not accept attachments.

Or
[url] http://www.eng-tips.com/threadminder.cfm?pid=799 [/url]
or from the ResearchGate discussion forum.

For advice on mechanics in general, including modeling,
you can consult
[url] http://imechanica.org/ [/url]

For an intro to subroutines get the file
[url] http://imechanica.org/files/Writing [/url] User Subroutines with ABAQUS.pdf

Good luck

Frank

Topic Tags
1 Reply
Share: