User-defined, ANSYS

Hi, I am trying to use USER300. I want to test it with a simple analysis. But couldn't do it. Seeing your post thought that you could help me. I am attaching code and error messages here. Thanks.

 

Error:

*** ERROR ***                           CP =       1.047   TIME= 19:44:49

 The inversion subroutine was supplied with a singular matrix.           

 *** ERROR ***                           CP =       1.047   TIME= 19:44:49

 Negative element Jacobian value 0 at element 25.  This is due to wrong  

 element order or bad mesh. 

!------------------------------------------------ 

 

 

/Prep7  

/UIS,MSGPOP,3

!-----------------

!ELEMENT TYPES

!-----------------

ET,1,USER300

USRELEM,8,3,BRICK,0,20,20,1,8,3,0

USRDOF,DEFINE,UX,UY,UZ,

!-----------------------

! Area meshing element

ET,2,200

KEYOPT,2,1,6

KEYOPT,2,2,0

!--------------------

ET,3,185

!-------------------------

!MATERIAL PROPERTIES

!--------------------

!MAIN Beam Material Properties

MP,DENS,1,2700 

MP,EX,1,69.5e9   

MP,PRXY,1,0.3

! Defining Beam geometry

B_L=0.071

B_W=0.01

B_T=0.00068

! Key point creation for areas to drag

! Beam

K,1,0,0,0,   

K,2,0,B_W,0,

K,3,0,B_W,B_T

K,4,0,0,B_T

K,1000,B_L,0,0

LSTR,1,1000

A,1,2,3,4

! Give even number for divisions for this line

LSEL,S,loc,Y,B_W/2

LESIZE,all, , ,6, , , , ,1 

alls

LSEL,S,loc,Z,B_T/2

LESIZE,all, , ,4, , , , ,1 

alls

LSEL,S,LINE,,1

LESIZE,all, , ,20, , , , ,1 

alls

! Meshing areas 1 & 2 with MESH200

TYPE,2

AMAP,1,1,2,3,4

! Extrudeing al beam using USER 300 and area 1

TYPE,3

MAT,1

VDRAG,       1, , , , , ,1

ESEL,S,TYPE,,3

EMODIF,ALL,TYPE,1,

ALLS

EMODIF,ALL,MAT,1,

ACLEAR,1

ALLS

NSEL,S,LOC,X,0

D,ALL,UX,0

D,ALL,UY,0

D,ALL,UZ,0  

ALLS

EPLOT

/VIEW,1,,-1 

/REP,FAST 

NPLOT

TIP=NODE(B_L,B_W/2,B_T/2)

F,TIP,FZ,-1

/SOLU

ANTYPE,0

SOLCONTROL,OFF

SOLVE