Search
Definition of tickn...

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,
3 stressOld, stateOld, enerInternOld, enerInelasOld,
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),
5 fieldOld(*), stressOld(nblock,ndir+nshr),eigVal(nblock,3),
6 stateOld(nblock,nstatev), enerInternOld(nblock),
7 enerInelasOld(nblock), tempNew(*), ,
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
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
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