-
Notifications
You must be signed in to change notification settings - Fork 26
/
main.f90
111 lines (108 loc) · 2.3 KB
/
main.f90
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
PROGRAM TEST_GENERAL_UMAT
use,intrinsic :: ISO_Fortran_env
INCLUDE 'ABA_PARAM.INC'
INCLUDE 'PARAM_UMAT.INC'
!C ADD COMMON BLOCKS HERE IF NEEDED ()
!C COMMON /KBLOCK/KBLOCK
PARAMETER(NTENS = 6, NSTATEV = NSDV, NPROPS = 6, NDI=3, NSHR=3)
PARAMETER(NOEL = 1, NPT = 8)
!
CHARACTER*8 CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATEV),DDSDDE(NTENS,NTENS),DDSDDT(NTENS), &
DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), &
PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
!
i=1.0d0
j=1.0d0
DO i=1,NTENS
DO j=1,NTENS
DDSDDE(i,j)=0.0D0
ENDDO
STRESS(i)=0.0D0
ENDDO
!
! DEFORMATION GRADIENT
DFGRD1(1,1)= 1.1D0
DFGRD1(1,2)= 0.0D0
DFGRD1(1,3)= 0.0D0
DFGRD1(2,1)= 0.0D0
DFGRD1(2,2)= 1.0D0/DFGRD1(1,1)
DFGRD1(2,3)= 0.0D0
DFGRD1(3,1)= 0.0D0
DFGRD1(3,2)= 0.0D0
DFGRD1(3,3)= 1.0D0/DFGRD1(1,1)
!
time(1)=0.d0
time(2)=0.d0
call UEXTERNALDB(0,0,time,0.D0,0,0)
!
! MATERIAL PROPERTIES
!
! k PENALTY PARAMETER
PROPS(1)=2.d0/100000.d0
! ISOTROPIC MATRIX
! C10=
PROPS(2)=1.00d0
! C01
PROPS(3)=1.00d0
!k1
PROPS(4)=1.00d0
!k2
PROPS(5)=0.1d0
!kappa
PROPS(6)=0.1d0
!
! !viscous parameters - maxwell
! ! v - number of dashpots
! PROPS(17)=0
! !tau1 %
! PROPS(18)=2.0d0
! !teta1
! PROPS(19)=0.835d0
! !tau2 %
! PROPS(20)=1.2d0
! !teta2
! PROPS(21)=7.0d0
! !tau3 %
! PROPS(22)=12.d0
! !teta3
! PROPS(23)=2.0d0
! !
STATEV=0.D0
!
erf=0.d0
RHO=0.D0
!
!
DFGRD1(1,1)= 1.0D0
DFGRD1(1,2)= 0.0D0
DFGRD1(1,3)= 0.0d0
DFGRD1(2,1)= 0.0d0
DFGRD1(2,2)= 1.0D0
DFGRD1(2,3)= 0.0d0
DFGRD1(3,1)= 0.0d0
DFGRD1(3,2)= 0.0d0
DFGRD1(3,3)= 1.0D0
!
!################################################################################################!
!! TENSILE MONOTONIC LOAD TEST
DFGRD1(1,1)= 1.2D0
DFGRD1(1,2)= 0.0D0
DFGRD1(1,3)= 0.0d0
DFGRD1(2,1)= 0.0d0
DFGRD1(2,2)= 1/sqrt(DFGRD1(1,1))
DFGRD1(2,3)= 0.0d0
DFGRD1(3,1)= 0.0d0
DFGRD1(3,2)= 0.0d0
DFGRD1(3,3)= 1/sqrt(DFGRD1(1,1))
!
CALL UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT, DRPLDE,DRPLDT,STRAN, &
DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,NDI,NSHR,NTENS,NSTATEV,PROPS, &
NPROPS,COORDS,DROT,PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
write(*,*) STRESS
write(*,*)
write(*,*) DDSDDE
close(150)
!################################################################################################!
!
END PROGRAM