-
Notifications
You must be signed in to change notification settings - Fork 1
/
2-Inputs_Outputs.txt
99 lines (94 loc) · 5.98 KB
/
2-Inputs_Outputs.txt
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
It contains 35 input properties (nprops=35), being the last one always ZERO by default. It is not used, but it must be specified equal to ZERO in the inputs.
E1=props(1) ! LONGITUDINAL ELASTIC MODULUS
E2=props(2) ! TRANSVERSAL ELASTIC MODULUS
E3=props(3) ! OUT-PLANE ELASTIC MODULUS
G12=props(4) ! SHEAR ELASTIC MODULUS
G13=props(5) ! SHEAR ELASTIC MODULUS
G23=props(6) ! SHEAR ELASTIC MODULUS
v12=props(7) ! POISSON COEFICIENT
v13=props(8) ! POISSON COEFICIENT
v23=props(9) ! POISSON COEFICIENT
Xt=props(10) ! LONGITUDINAL TENSILE RESISTENCE
Xc=props(11) ! LONGITUDINAL COMPRESSIVE RESISTENCE
Yt=props(12) ! TRANSVERSAL TENSILE RESISTENCE
Yc=props(13) ! TRANSVERSAL COMPRESSIVE RESISTENCE
Zt=props(14) ! OUT-PLANE TENSILE RESISTENCE
Zc=props(15) ! OUT-PLANE COMPRESSIVE RESISTENCE
Sl=props(16) ! LONGITUDINAL SHEAR RESISTENCE
St=props(17) ! TRANSVERSAL SHEAR RESISTENCE
Si=props(18) ! OUT-PLANE SHEAR RESISTENCE
Gft=props(19) ! FRACTURE ENERGY FOR FIBRE TENSION
Gfc=props(20) ! FRACTURE ENERGY FOR FIBRE COMPRESSION
Gmt=props(21) ! FRACTURE ENERGY FOR MATRIX TENSION
Gmc=props(22) ! FRACTURE ENERGY FOR MATRIX COMPRESSION
Git=props(23) ! FRACTURE ENERGY FOR INTERLAMINAR TENSION
Gic=props(24) ! FRACTURE ENERGY FOR INTERLAMINAR COMPRESSION
etaft=props(25) ! VISCOUS REGULARIZATION COEFICIENT
dmax=props(26) ! MAXIMUM ALLOWED DAMAGE FOR CONVERGENCY (insert 1.0)
alpha_f=props(27) ! INTERACTION COEFICIENT FOR FIBRE TO SIMULATE 2D BEHAVIOUR (insert 1.0)
alpha_i=props(28) ! INTERACTION COEFICIENT FOR INTERLAMINAR TO SIMULATE 2D BEHAVIOUR (insert 1.0)
presft=props(29) ! FIBRE TENSILE RESIDUAL STRESS (10% then insert 0.1)
presfc=props(30) ! FIBRE COMPRESSIVE RESIDUAL STRESS
presmt=props(31) ! MATRIX TENSILE RESIDUAL STRESS
presmc=props(32) ! MATRIX COMPRESSIVE RESIDUAL STRESS
presit=props(33) ! INTERLAMINAR TENSILE RESIDUAL STRESS
presic=props(34) ! INTERLAMINAR COMPRESSIVE RESIDUAL STRESS
props(35) ! NOT USED
It contains 51 output state variables (nstatv=51), being the last 9 irrelevant at the moment, and should be disregarded.
dfto=statev(1) ! FIBRE TENSION DAMAGE FROM PREVIOUS INCREMENT
dfco=statev(2) ! FIBRE COMPRESSION DAMAGE FROM PREVIOUS INCREMENT
dmto=statev(3) ! MATRIX TENSION DAMAGE FROM PREVIOUS INCREMENT
dmco=statev(4) ! MATRIX COMPRESSION DAMAGE FROM PREVIOUS INCREMENT
dito=statev(5) ! INTERLAMINAR TENSION DAMAGE FROM PREVIOUS INCREMENT
dico=statev(6) ! INTERLAMINAR COMPRESSION DAMAGE FROM PREVIOUS INCREMENT
ds12o=statev(7) ! SHEAR 12 DAMAGE FROM PREVIOUS INCREMENT
ds13o=statev(8) ! SHEAR 13 DAMAGE FROM PREVIOUS INCREMENT
ds23o=statev(9) ! SHEAR 23 DAMAGE FROM PREVIOUS INCREMENT
Fafto=statev(10) ! FIBRE TENSION CRITERIA FROM PREVIOUS INCREMENT
Fafco=statev(11) ! FIBRE COMPRESSION CRITERIA FROM PREVIOUS INCREMENT
Famto=statev(12) ! MATRIX TENSION CRITERIA FROM PREVIOUS INCREMENT
Famco=statev(13) ! MATRIX COMPRESSION CRITERIA FROM PREVIOUS INCREMENT
Faito=statev(14) ! INTERLAMINAR TENSION CRITERIA FROM PREVIOUS INCREMENT
Faico=statev(15) ! INTERLAMINAR COMPRESSION CRITERIA FROM PREVIOUS INCREMENT
dvfto=statev(16) ! VISCOUS FIBRE TENSION DAMAGE FROM PREVIOUS INCREMENT
dvfco=statev(17) ! VISCOUS FIBRE COMRESSION DAMAGE FROM PREVIOUS INCREMENT
dvmto=statev(18) ! VISCOUS MATRIX TENSION DAMAGE FROM PREVIOUS INCREMENT
dvmco=statev(19) ! VISCOUS MATRIX COMPRESSION DAMAGE FROM PREVIOUS INCREMENT
dvito=statev(20) ! VISCOUS INTERLAMINAR TENSION DAMAGE FROM PREVIOUS INCREMENT
dvico=statev(21) ! VISCOUS INTERLAMINAR COMPRESSION DAMAGE FROM PREVIOUS INCREMENT
dvs12o=statev(22) ! VISCOUS SHEAR 12 DAMAGE FROM PREVIOUS INCREMENT
dvs13o=statev(23) ! VISCOUS SHEAR 13 DAMAGE FROM PREVIOUS INCREMENT
dvs23o=statev(24) ! VISCOUS SHEAR 23 DAMAGE FROM PREVIOUS INCREMENT
seqft0=statev(25) ! EQUIVALENTE FIBRE TENSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
seqfc0=statev(26) ! EQUIVALENTE FIBRE COMPRESSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
seqmt0=statev(27) ! EQUIVALENTE MATRIX TENSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
seqmc0=statev(28) ! EQUIVALENTE MATRIX COMPRESSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
seqit0=statev(29) ! EQUIVALENTE INTERLAMINAR TENSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
seqic0=statev(30) ! EQUIVALENTE INTERLAMINAR COMPRESSION AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqft0=statev(31) ! EQUIVALENTE FIBRE TENSION DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqfc0=statev(32) ! EQUIVALENTE FIBRE COMPRESSION DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqmt0=statev(33) ! EQUIVALENTE MATRIX TENSION DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqmc0=statev(34) ! EQUIVALENTE MATRIX COMPRESSON DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqit0=statev(35) ! EQUIVALENTE INTERLAMINAR TENSION DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
ueqic0=statev(36) ! EQUIVALENTE INTERLAMINAR COMPRESSON DISPLACEMENT AT THE ONSET OF DAMAGE FROM PREVIOUS INCREMENT
if (ndi==3) then !!!! FOR 3D ANALYSIS
sijeo(1)=statev(37) ! EFFECTIVE STRESS 11 FROM PREVIOUS INCREMENT
sijeo(2)=statev(38) ! EFFECTIVE STRESS 22 FROM PREVIOUS INCREMENT
sijeo(3)=statev(39) ! EFFECTIVE STRESS 33 FROM PREVIOUS INCREMENT
sijeo(4)=statev(40) ! EFFECTIVE STRESS 12 FROM PREVIOUS INCREMENT
sijeo(5)=statev(41) ! EFFECTIVE STRESS 13 FROM PREVIOUS INCREMENT
sijeo(6)=statev(42) ! EFFECTIVE STRESS 23 FROM PREVIOUS INCREMENT
else !!!! FOR 2D ANALYSIS
sijeo(1)=statev(37) ! EFFECTIVE STRESS 11 FROM PREVIOUS INCREMENT
sijeo(2)=statev(38) ! EFFECTIVE STRESS 22 FROM PREVIOUS INCREMENT
sijeo(3)=statev(39) ! EFFECTIVE STRESS 12 FROM PREVIOUS INCREMENT
end if
statev(43) ! NOT USED
statev(44) ! NOT USED
statev(45) ! NOT USED
statev(46) ! NOT USED
statev(47) ! NOT USED
statev(48) ! NOT USED
statev(49) ! NOT USED
statev(50) ! NOT USED
statev(51) ! NOT USED