-
Notifications
You must be signed in to change notification settings - Fork 4
/
USDFLD_cartesian_asym.for
54 lines (42 loc) · 2.02 KB
/
USDFLD_cartesian_asym.for
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
C ABAQUS Subroutine for functionally graded material (FGM) in a cartesian coordinate system 2D
C
C
C More Infos can be found at
C https://classes.engineering.wustl.edu/2009/spring/mase5513/abaqus/docs/v6.6/books/sub/default.htm?startat=ch01s01asb40.html#sub-xsl-usdfld
C
C by Dr. Ronald Wagner - University Braunschweig, Germany, [email protected]
C---------------------------------------------------------------------------------------
C Start of Base code, DO NOT change
C---------------------------------------------------------------------------------------
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),T(3,3),TIME(2)
DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C---------------------------------------------------------------------------------------
C End of Base code
C---------------------------------------------------------------------------------------
C---------------------------------------------------------------------------------------
C Start of USER code
C---------------------------------------------------------------------------------------
c read x, y, z coordinate
x = COORD(1)
y = COORD(2)
z = COORD(3)
c calculate Young module
E = y
c define depend state variable (E(z))
FIELD(1) = E
c save E value for state dependend variable
C STATEV(1) = 19164.61 + (81303.16 - 19164.61)/(1 + (E/1.733635)**4.936655)
STATEV(1) = -19750*E+ 80000
C---------------------------------------------------------------------------------------
C End of USER code
C---------------------------------------------------------------------------------------
RETURN
END