-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathmatutil.for
executable file
·114 lines (114 loc) · 4.56 KB
/
matutil.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
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
112
113
114
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C11111111122222222223333333333444444444455555555556666666666777777777777
C23456789012345678901234567890123456789012345678901234567890123456789012
C C
C---C O N S T I T U T I V E M A T E R I A L S S U B R O U T I N E S--C
C C
C23456789012345678901234567890123456789012345678901234567890123456789012
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C11111111122222222223333333333444444444455555555556666666666777777777777
C23456789012345678901234567890123456789012345678901234567890123456789012
C C
C SUBROUTINE UMAT C
C UMAT.for Classical Isotropic Elasticity C
C UNIVERSIDAD EAFIT C
C 2011 C
C CANNOT BE USED FOR PLANE STRESS C
C NTENS: LENGTH OF STRESS VECTOR C
C NDI: NUMBER OF NORMAL STRESS COMPONENTS C
C C
C PROPS(1) - VS C
C PROPS(2) - VP C
C PROPS(3) - R0 C
C C
C23456789012345678901234567890123456789012345678901234567890123456789012
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SUBROUTINE UMAT(DDSDDE,NDI,NTENS,PROPS,NPROPS)
C
IMPLICIT REAL*8(A-H,O-Z)
C
DIMENSION DDSDDE(NTENS, NTENS),PROPS(NPROPS)
C
PARAMETER (ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, SIX=6.0D0,
1 ENUMAX=0.4999D0)
C
C **********************************************************************
C I S O T R O P I C E L A S T I C I T Y
C P L A N E S T R A I N A N A L Y S I S
C **********************************************************************
C
C Wave propagation velocities.
C
VP=PROPS(1)
VS=PROPS(2)
RO=PROPS(3)
C
VPS=VP*VP
VSS=VS*VS
C
C Elasticity tensor.
C
DO K1=1, NDI
DO K2=1, NDI
DDSDDE(K2, K1)=(VPS-2.0*VSS)*RO
END DO
DDSDDE(K1, K1)=VPS*RO
END DO
DO K1=NDI+1, NTENS
DDSDDE(K1, K1)=VSS*RO
END DO
C
RETURN
C
END
C
C23456789012345678901234567890123456789012345678901234567890123456789012
C C
C SUBROUTINE UMATSH C
C Constitutive Matrix for SH wave C
C UNIVERSIDAD EAFIT C
C 2014 C
C C
C NTENS: LENGTH OF STRESS VECTOR C
C NDI: NUMBER OF NORMAL STRESS COMPONENTS C
C C
C PROPS(1) - VS C
C PROPS(2) - VP C
C PROPS(3) - R0 C
C C
C23456789012345678901234567890123456789012345678901234567890123456789012
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SUBROUTINE UMATSH(DDSDDE,NDI,NTENS,PROPS,NPROPS)
C
IMPLICIT REAL*8(A-H,O-Z)
C
DIMENSION DDSDDE(NTENS, NTENS),PROPS(NPROPS)
C
PARAMETER (ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, SIX=6.0D0,
1 ENUMAX=0.4999D0)
C
C **********************************************************************
C S H W A V E A N A L Y S I S
C **********************************************************************
C
C Wave propagation velocities.
C
VP=PROPS(1)
VS=PROPS(2)
RO=PROPS(3)
C
UMIU=VS*VS*RO
C
C Elasticity tensor for SH.
C
DO K1=1, NDI
DDSDDE(K1, K1)=UMIU
END DO
C
RETURN
C
END
C