UELMAT
1,2
*
ELEMENT, TYPE=U1, ELSET=SOLID
1, 1,2,3,4
...
*
UEL PROPERTY, ELSET=SOLID, MATERIAL=MAT
...
*
MATERIAL, NAME=MAT
*
ELASTIC
7.00E+010, 0.33
The user element defi ned above is a 4-node, fully inte grated plane strain element, similar to the Abaqus
CPE4 element.
The next excerpt shows th
e listing of the user subroutine. I nside the subroutine, a loop over the
integration points is perfo
rmed. For each integration point the utility routine MATERIAL_LIB_MECH
is called, which returns stres
s and Jacobian at the integration point. These quantities are used to compute
the right-hand-side vector and th
e element Jacobian.
c***********************************************************
subroutine uelmat(rhs,amatrx,svars,energy,ndofel,nrhs,
1 nsvars,props,nprops,coords,mcrd,nnode,u,du,
2 v,a,jtype,time,dtime,kstep,kinc,jelem,params,
3 ndload,jdltyp,adlmag,predef,npredf,lflags,mlvarx,
4 ddlmag,mdload,pnewdt,jprops,njpro,period,
5 materiallib)
c
include 'aba_param.inc'
C
dimension rhs(mlvarx,*), amatrx(ndofel, ndofel), props(*),
1 svars(*), energy(*), coords(mcrd, nnode), u(ndofel),
2 du(mlvarx,*), v(ndofel), a(ndofel), time(2), params(*),
3 jdltyp(mdload,*), adlmag(mdload,*), ddlmag(mdload,*),
4 predef(2, npredf, nnode), lflags(*), jprops(*)
parameter (zero=0.d0, dmone=-1.0d0, one=1.d0, four=4.0d0,
1 fourth=0.25d0,gaussCoord=0.577350269d0)
parameter (ndim=2, ndof=2, nshr=1,nnodemax=4,
1 ntens=4, ninpt=4, nsvint=4)
c
c ndim ... number of spatial dimensions
c ndof ... number of degrees of freedom per node
c nshr ... number of shear stress component
c ntens ... total number of stress tensor components
c (=ndi+nshr)
1.1.24–12
Abaqus ID:
Printed on: