35
36
37
38
42
43
44
45#include "implicit_f.inc"
46
47
48
49
50
51
52 INTEGER NTYP
53 INTEGER DEF_INTER(100)
54 TYPE(SUBMODEL_DATA) LSUBMODEL(NSUBMOD)
55
56
57
58#include "units_c.inc"
59
60
61
62 INTEGER
63 . IGSTI,IREM24I2,ILEV,IDEL24,INTPLY,
64 . IEDGE,IGAP0,IPEN0,INACTI,MFROT,IFQ,
65 . IDSENS,INTFRIC,IPRINT,IASSIGN
67 . eten_m,edg_angl,gapmax_s,gapmax_m,stmin,stmax,penmax,penmin,
68 . visc,viscf,bumult,
alpha,gap,startt,stopt,stfac,fric
69 CHARACTER(LEN=NCHARFIELD)::BCFLAG, BCFLAGM
70 CHARACTER(LEN=NCHARKEY)::KEY,KEY1
71 LOGICAL
72
73
74
75
76
77 ilev = 0
78 igsti = 0
79 idel24= 0
80 intply = 0
81 inacti = 0
82 idsens = 0
83 mfrot=0
84 ifq=0
85 intfric = 0
86
87 fric = zero
88 gap = zero
89 startt = zero
90 stopt=ep30
91 viscf = zero
92 visc = zero
93
94 iassign = 0
95
96
97 iprint = 0
98 is_available = .false.
99
100
101
102 CALL hm_get_intv(
'Istf',igsti,is_available,lsubmodel)
103 CALL hm_get_intv(
'Irem_i2',irem24i2,is_available,lsubmodel)
104 CALL hm_get_intv(
'TYPE24_Idel',idel24,is_available,lsubmodel)
105 CALL hm_get_intv(
'TYPE24_Iedge',iedge,is_available,lsubmodel)
106 CALL hm_get_intv(
'Itied',igap0,is_available,lsubmodel)
107 CALL hm_get_intv(
'Ipen',ipen0,is_available,lsubmodel)
108 CALL hm_get_intv(
'Irs',inacti,is_available,lsubmodel)
109
110 key1='IGSTI'
111 CALL definter(key1 ,igsti ,iassign ,iprint ,
112 . ntyp ,def_inter)
113 key1='IDEL24'
114 CALL definter(key1 ,idel24 ,iassign ,iprint ,
115 . ntyp ,def_inter)
116 key1='IREM24I2'
117 CALL definter(key1 ,irem24i2 ,iassign ,iprint ,
118 . ntyp ,def_inter)
119
120 key1='IEDGE'
121 CALL definter(key1 ,iedge ,iassign ,iprint ,
122 . ntyp ,def_inter)
123
124 key1='IGAP0'
125 CALL definter(key1 ,igap0 ,iassign ,iprint ,
126 . ntyp ,def_inter)
127
128 key1='IPEN0'
129 CALL definter(key1 ,ipen0 ,iassign ,iprint ,
130 . ntyp ,def_inter)
131
132 key1='INACTI'
133 CALL definter(key1 ,inacti ,iassign ,iprint ,
134 . ntyp ,def_inter)
135
136
137
138
139 iprint = 1
140
141 key1='IGSTI'
142 CALL definter(key1 ,igsti ,iassign ,iprint ,
143 . ntyp ,def_inter)
144
145 key1='IDEL24'
146 CALL definter(key1 ,idel24 ,iassign ,iprint ,
147 . ntyp ,def_inter)
148
149 key1='IEDGE'
150 CALL definter(key1 ,iedge ,iassign ,iprint ,
151 . ntyp ,def_inter)
152
153 key1='IGAP0'
154 CALL definter(key1 ,igap0 ,iassign ,iprint ,
155 . ntyp ,def_inter)
156
157 key1='IPEN0'
158 CALL definter(key1 ,ipen0 ,iassign ,iprint ,
159 . ntyp ,def_inter)
160
161 key1='INACTI'
162 CALL definter(key1 ,inacti ,iassign ,iprint ,
163 . ntyp ,def_inter)
164
165 key1='IREM24I2'
166 CALL definter(key1 ,irem24i2 ,iassign ,iprint ,
167 . ntyp ,def_inter)
168
169 WRITE(iout,5024) igsti,idel24,iedge,igap0,irem24i2,inacti
170
171
172 RETURN
173
174 5024 FORMAT(/
175 . ' INTERFACE TYPE 24 ' //,
176 . ' STIFFNESS FORMULATION. . . . . . . . . . ',i5/,
177 . ' DELETION FLAG ON FAILURE OF MAIN ELEMENT',i5/,
178 . ' EDGE FLAG . . . . . . . . . . . . . . . . ',i5/,
179 . ' GAP MODIFICATION FLAG FOR SECONDARY SHELL ',/,
180 . ' NODES ON THE FREE EDGES . . . . . . . . . ',i5/,
181 . ' DE-ACTIVATION FLAG IF CONFLICTING W/ TYPE2',i5/,
182 . ' DE-ACTIVATION OF INITIAL PENETRATIONS . . ',i5/,
183 . ' FRICTION FORMULATION. . . . . . . . . . . ',i5/)
184
185
subroutine definter(key, ival, flag, iprint, ityp, def_inter)
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
integer, parameter ncharkey
integer, parameter ncharfield