OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_prop25.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| hm_read_prop25 ../starter/source/properties/spring/hm_read_prop25.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_properties ../starter/source/properties/hm_read_properties.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.f
31!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.f
32!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
33!||--- uses -----------------------------------------------------
34!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
35!|| message_mod ../starter/share/message_module/message_mod.F
36!|| submodel_mod ../starter/share/modules1/submodel_mod.F
37!||====================================================================
38 SUBROUTINE hm_read_prop25(GEO, IGEO, UNITAB, ISKN,IGTYP,
39 . IG,PROP_TAG,IDTITL,LSUBMODEL,SUB_ID)
40C-----------------------------------------------
41 USE unitab_mod
42 USE elbuftag_mod
43 USE message_mod
44 USE submodel_mod
45 USE names_and_titles_mod , ONLY : nchartitle
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "units_c.inc"
54#include "param_c.inc"
55#include "com04_c.inc"
56#include "sphcom.inc"
57#include "tablen_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
62 INTEGER IGEO(NPROPGI),ISKN(LISKN,*),IGTYP,IG,SUB_ID
63C REAL
64 my_real
65 . geo(npropg)
66 CHARACTER(LEN=NCHARTITLE) :: IDTITL
67 TYPE(prop_tag_) , DIMENSION(0:MAXPROP) :: PROP_TAG
68 TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 INTEGER J, IFUNC, IFUNC2,IFUNC3,IECROU, IFV, ISK,
73 . isens,ifl,ifail,ileng,ifail2,k
74C REAL
75 my_real
76 . a, b, d, e, f, xm, xin, xk, xc, dn, dx, fwv,lscale, fscale,
77 . vt0, vr0, cc(6), cn(6), xa(6), xb(6),gf3,pun,a_unit,d_unit,f_unit,
78 . e_unit,lscale_unit,gf3_unit,crit_scale(4),vr0_unit,vt0_unit
79 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
80C=======================================================================
81C
82 pun = em01
83 fwv = zero
84 ifail2 = 0
85C
86 is_encrypted = .false.
87 is_available = .false.
88C
89 igeo(1)=ig
90 igeo(11)=igtyp
91 geo(12) =igtyp+pun
92C
93C--------------------------------------------------
94C EXTRACT DATA (IS OPTION CRYPTED)
95C--------------------------------------------------
96 CALL hm_option_is_encrypted(is_encrypted)
97C--------------------------------------------------
98C EXTRACT DATAS (INTEGER VALUES)
99C--------------------------------------------------
100 CALL hm_get_intv('SKEW_CSID',isk,is_available,lsubmodel)
101 IF(isk == 0 .AND. sub_id /= 0 ) isk = lsubmodel(sub_id)%SKEW
102 CALL hm_get_intv('ISENSOR',isens,is_available,lsubmodel)
103 CALL hm_get_intv('ISFLAG',ifl,is_available,lsubmodel)
104 CALL hm_get_intv('Ifail',ifail,is_available,lsubmodel)
105 CALL hm_get_intv('Ileng',ileng,is_available,lsubmodel)
106 CALL hm_get_intv('Ifail2',ifail2,is_available,lsubmodel)
107C--------------------------------------------------
108C EXTRACT DATAS (REAL VALUES)
109C--------------------------------------------------
110 CALL hm_get_floatv('MASS',xm,is_available,lsubmodel,unitab)
111 CALL hm_get_floatv('INERTIA',xin,is_available,lsubmodel,unitab)
112C----
113 IF (ifl == 1) isens=-isens
114C----
115 DO k=0,numskw+min(1,nspcond)*numsph+nsubmod
116 IF(isk == iskn(4,k+1)) THEN
117 isk=k+1
118 GO TO 100
119 ENDIF
120 ENDDO
121 CALL ancmsg(msgid=137,anmode=aninfo,msgtype=msgerror,
122 . c1='PROPERTY',
123 . c2='PROPERTY',
124 . i1=igeo(1),i2=isk,c3=idtitl)
125100 CONTINUE
126C
127 igeo(3)=isens
128 geo(79)=ifail
129 geo(80)=ifl
130 geo(93)=ileng
131 geo(95)=ifail2
132 geo(1) =xm
133 geo(2) =isk+em01
134 geo(8) =6
135 geo(9) =xin
136 igeo(2)=isk
137C----
138 IF(is_encrypted)THEN
139 WRITE(iout,1000)ig
140 1000 FORMAT(
141 & 5x,'SPRING PROPERTY SET'/,
142 & 5x,'-------------------'/,
143 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
144 & 5x,'CONFIDENTIAL DATA'//)
145 ELSE
146 WRITE(iout,1801)ig,xm,xin,iskn(4,isk),abs(isens),ifl,ileng
147 ENDIF
148C
149!-------------------------------------------------------
150! Translations
151!-------------------------------------------------------
152C
153!-----------------
154 ! Tension
155!-----------------
156C--------------------------------------------------
157C EXTRACT DATAS (INTEGER VALUES)
158C--------------------------------------------------
159 CALL hm_get_intv('FUN_A1',ifunc,is_available,lsubmodel)
160 CALL hm_get_intv('HFLAG1',iecrou,is_available,lsubmodel)
161 CALL hm_get_intv('FUN_B1',ifv,is_available,lsubmodel)
162 CALL hm_get_intv('FUN_C1',ifunc2,is_available,lsubmodel)
163 CALL hm_get_intv('FUN_D1',ifunc3,is_available,lsubmodel)
164C--------------------------------------------------
165C EXTRACT DATAS (REAL VALUES)
166C--------------------------------------------------
167 CALL hm_get_floatv('STIFF1',xk,is_available,lsubmodel,unitab)
168 CALL hm_get_floatv('DAMP1',xc,is_available,lsubmodel,unitab)
169 CALL hm_get_floatv('acoeft1',A,IS_AVAILABLE,LSUBMODEL,UNITAB)
170 CALL HM_GET_FLOATV('bcoeft1',B,IS_AVAILABLE,LSUBMODEL,UNITAB)
171 CALL HM_GET_FLOATV('dcoeft1',D,IS_AVAILABLE,LSUBMODEL,UNITAB)
172 CALL HM_GET_FLOATV('min_rup1',DN,IS_AVAILABLE,LSUBMODEL,UNITAB)
173 CALL HM_GET_FLOATV('max_rup1',DX,IS_AVAILABLE,LSUBMODEL,UNITAB)
174 CALL HM_GET_FLOATV('prop_tens_f',F,IS_AVAILABLE,LSUBMODEL,UNITAB)
175 CALL HM_GET_FLOATV('prop_tens_e',E,IS_AVAILABLE,LSUBMODEL,UNITAB)
176 CALL HM_GET_FLOATV('scale1',LSCALE,IS_AVAILABLE,LSUBMODEL,UNITAB)
177 CALL HM_GET_FLOATV('f_emc_a',GF3,IS_AVAILABLE,LSUBMODEL,UNITAB)
178C
179 !units for default values
180 CALL HM_GET_FLOATV_DIM('acoeft1',A_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
181 CALL HM_GET_FLOATV_DIM('dcoeft1',D_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
182 CALL HM_GET_FLOATV_DIM('prop_tens_f',F_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
183 CALL HM_GET_FLOATV_DIM('prop_tens_e',E_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
184 CALL HM_GET_FLOATV_DIM('scale1',LSCALE_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
185 CALL HM_GET_FLOATV_DIM('f_emc_a',GF3_UNIT,IS_AVAILABLE,LSUBMODEL,UNITAB)
186 CALL HM_GET_FLOATV_DIM('min_rup1',CRIT_SCALE(1),IS_AVAILABLE,LSUBMODEL,UNITAB)
187C
188C----
189.AND..OR. IF(IECROU==4(IFUNC==0IFUNC2==0))THEN
190 CALL ANCMSG(MSGID=231,
191 . MSGTYPE=MSGERROR,
192 . ANMODE=ANINFO_BLIND_1,
193 . I1=IG,
194 . C1=IDTITL)
195 ENDIF
196.AND. IF(IECROU==4 ISK==ZERO)THEN
197 CALL ANCMSG(MSGID=230,
198 . MSGTYPE=MSGERROR,
199 . ANMODE=ANINFO_BLIND_1,
200 . I1=IG,
201 . C1=IDTITL)
202 ENDIF
203.AND..OR. IF(IECROU==5(IFUNC==0IFUNC2==0))THEN
204 CALL ANCMSG(MSGID=231,
205 . MSGTYPE=MSGERROR,
206 . ANMODE=ANINFO_BLIND_1,
207 . I1=IG,
208 . C1=IDTITL)
209 ENDIF
210.AND..OR. IF(IECROU==6(IFUNC==0IFUNC2==0))THEN
211 CALL ANCMSG(MSGID=1057,
212 . MSGTYPE=MSGERROR,
213 . ANMODE=ANINFO_BLIND_1,
214 . I1=IG,
215 . C1=IDTITL)
216 ENDIF
217.AND. IF(IECROU==7IFUNC==0)THEN
218 CALL ANCMSG(MSGID=1058,
219 . MSGTYPE=MSGERROR,
220 . ANMODE=ANINFO_BLIND_1,
221 . I1=IG,
222 . C1=IDTITL)
223
224.AND. ELSEIF(IECROU==7IFUNC2==0)THEN
225 CALL ANCMSG(MSGID=1059,
226 . MSGTYPE=MSGWARNING,
227 . ANMODE=ANINFO_BLIND_1,
228 . I1=IG,
229 . C1=IDTITL,
230 . I2=IECROU)
231 IECROU = 2
232 ENDIF
233.AND..AND.! IF(IFUNC /= 0 IECROU >= 1 XK == ZERO)THEN
234! CALL ANCMSG(MSGID=230,
235! . MSGTYPE=MSGERROR,
236! . ANMODE=ANINFO_BLIND_1,
237! . I1=IG,
238! . C1=IDTITL)
239! ENDIF
240.AND..AND. IF (IFUNC == 0 A /= ZERO A /= ONE) THEN
241 CALL ANCMSG(MSGID=663,
242 . MSGTYPE=MSGWARNING,
243 . ANMODE=ANINFO_BLIND_1,
244 . I1=IG,
245 . C1=IDTITL)
246 ENDIF
247C----
248 IF (A == ZERO) A = ONE * A_UNIT
249 IF (D == ZERO) D = ONE * D_UNIT
250 IF (E == ZERO) E = ONE * E_UNIT
251 IF (F == ZERO) F = ONE * F_UNIT
252 IF (LSCALE == ZERO) LSCALE = ONE * LSCALE_UNIT
253 IF (GF3 == ZERO) GF3 = ONE * GF3_UNIT
254 IF (IFUNC == 0) THEN
255 A = ONE
256 B = ZERO
257 E = ZERO
258 ENDIF
259C
260.OR. IF ((IFAIL2 == 0)(IFAIL2 ==1)) THEN
261 DN = DN * LSCALE / CRIT_SCALE(1)
262 DX = DX * LSCALE / CRIT_SCALE(1)
263 ENDIF
264 IF (DN == ZERO) DN=-EP30* CRIT_SCALE(1)
265 IF (DX == ZERO) DX= EP30* CRIT_SCALE(1)
266.OR. IF ((IFAIL2 == 0)(IFAIL2 ==1)) CRIT_SCALE(1) = LSCALE / CRIT_SCALE(1)
267C----
268 GEO(3) =XK / A
269 GEO(4) =XC
270 GEO(41)=A
271 GEO(42)=B
272 GEO(43)=D
273 GEO(40)=E
274 GEO(132)= GF3
275 GEO(44)=ONE/F
276 GEO(39)=ONE/LSCALE
277 GEO(65)=DN
278 GEO(66)=DX
279C-- If H=6 - additional internal variables must be stored in UVAR - - GEO(25) = NUVAR
280 IF (IECROU == 6) THEN
281 GEO(25) = 6
282 ENDIF
283C
284 IGEO(101) =IECROU
285 IGEO(102) =IFUNC
286 IGEO(103) =IFUNC2
287 IGEO(104) =IFV
288 IGEO(119) = IFUNC3
289C----
290.NOT. IF(IS_ENCRYPTED)THEN
291 IF (IECROU /= 5) THEN
292 IF (IFAIL2 == 3) THEN
293 WRITE(IOUT,1813)'tension',XK,XC,IFUNC,LSCALE,IFUNC2,F,IECROU,
294 . A,B,D,E,GF3,IFV,IFUNC3,DX
295 ELSEIF (IFAIL2 == 2) THEN
296 WRITE(IOUT,1812)'tension',XK,XC,IFUNC,LSCALE,IFUNC2,F,IECROU,
297 . A,B,D,E,GF3,IFV,IFUNC3,DN,DX
298 ELSE
299 WRITE(IOUT,1810)'tension',XK,XC,IFUNC,LSCALE,IFUNC2,F,IECROU,
300 . A,B,D,E,GF3,IFV,IFUNC3,DN,DX
301 ENDIF
302 ELSE
303 IF (IFAIL2 == 3) THEN
304 WRITE(IOUT,1823)'tension',XK,XC,IFUNC,LSCALE,IFUNC2,F,IECROU,
305 . A,B,D,E,GF3,IFV,IFUNC3,DX
306 ELSEIF (IFAIL2 == 2) THEN
307 WRITE(IOUT,1822)'tension',XK,XC,IFUNC,LSCALE,IFUNC2,F,IECROU,
308 . A,B,D,E,GF3,IFV,IFUNC3,DN,DX
309 ELSE
310 WRITE(IOUT,1820)'tension',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
311 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
312 ENDIF
313 ENDIF
314 ENDIF
315C
316!-----------------
317 ! Cisaillement
318!-----------------
319C--------------------------------------------------
320C EXTRACT DATAS (INTEGER VALUES)
321C--------------------------------------------------
322 CALL hm_get_intv('FUN_A2',ifunc,is_available,lsubmodel)
323 CALL hm_get_intv('HFLAG2',iecrou,is_available,lsubmodel)
324 CALL hm_get_intv('FUN_B2',ifv,is_available,lsubmodel)
325 CALL hm_get_intv('FUN_C2',ifunc2,is_available,lsubmodel)
326 CALL hm_get_intv('FUN_D2',ifunc3,is_available,lsubmodel)
327C--------------------------------------------------
328C EXTRACT DATAS (REAL VALUES)
329C--------------------------------------------------
330 CALL hm_get_floatv('STIFF2',xk,is_available,lsubmodel,unitab)
331 CALL hm_get_floatv('DAMP2',xc,is_available,lsubmodel,unitab)
332 CALL hm_get_floatv('Acoeft2',a,is_available,lsubmodel,unitab)
333 CALL hm_get_floatv('Bcoeft2',b,is_available,lsubmodel,unitab)
334 CALL hm_get_floatv('Dcoeft2',d,is_available,lsubmodel,unitab)
335 CALL hm_get_floatv('MIN_RUP2',dn,is_available,lsubmodel,unitab)
336 CALL hm_get_floatv('MAX_RUP2',dx,is_available,lsubmodel,unitab)
337 CALL hm_get_floatv('Prop_Shear_F',f,is_available,lsubmodel,unitab)
338 CALL hm_get_floatv('Prop_Shear_E',e,is_available,lsubmodel,unitab)
339 CALL hm_get_floatv('scale2',lscale,is_available,lsubmodel,unitab)
340 CALL hm_get_floatv('F_EMC_n',gf3,is_available,lsubmodel,unitab)
341C!units for default values
342 CALL hm_get_floatv_dim('MIN_RUP2',crit_scale(2),is_available,lsubmodel,unitab)
343C----
344 IF(iecrou/=0 .AND. iecrou/=1 .AND. iecrou/=5.AND.
345 . iecrou/=6.AND. iecrou/=7) THEN
346 CALL ancmsg(msgid=905,
347 . msgtype=msgerror,
348 . anmode=aninfo_blind_1,
349 . i1=ig,
350 . c1=idtitl)
351 ENDIF
352! IF(IECROU>=1 .AND. IFUNC/=0 .AND. XK == ZERO)THEN
353! CALL ANCMSG(MSGID=230,
354! . MSGTYPE=MSGERROR,
355! . ANMODE=ANINFO_BLIND_1,
356! . I1=IG,
357! . C1=IDTITL)
358! ENDIF
359 IF(iecrou==5.AND.(ifunc==0.OR.ifunc2==0))THEN
360 CALL ancmsg(msgid=231,
361 . msgtype=msgerror,
362 . anmode=aninfo_blind_1,
363 . i1=ig,
364 . c1=idtitl)
365 ENDIF
366 IF(iecrou==6.AND.(ifunc==0.OR.ifunc2==0))THEN
367 CALL ancmsg(msgid=1057,
368 . msgtype=msgerror,
369 . anmode=aninfo_blind_1,
370 . i1=ig,
371 . c1=idtitl)
372 ENDIF
373 IF(iecrou==7.AND.ifunc==0)THEN
374 CALL ancmsg(msgid=1058,
375 . msgtype=msgerror,
376 . anmode=aninfo_blind_1,
377 . i1=ig,
378 . c1=idtitl)
379
380 ELSEIF(iecrou==7.AND.ifunc2==0)THEN
381 CALL ancmsg(msgid=1059,
382 . msgtype=msgwarning,
383 . anmode=aninfo_blind_1,
384 . i1=ig,
385 . c1=idtitl,
386 . i2=iecrou)
387 iecrou = 2
388 ENDIF
389 IF (ifunc == 0 .AND. a /= zero .AND. a /= one) THEN
390 CALL ancmsg(msgid=663,
391 . msgtype=msgwarning,
392 . anmode=aninfo_blind_1,
393 . i1=ig,
394 . c1=idtitl)
395 ENDIF
396C----
397 IF (a == zero) a = one * a_unit
398 IF (d == zero) d = one * d_unit
399 IF (e == zero) e = one * e_unit
400 IF (f == zero) f = one * f_unit
401 IF (lscale == zero) lscale = one * lscale_unit
402 IF (gf3 == zero) gf3 = one * gf3_unit
403 IF (ifunc == 0) THEN
404 a = one
405 b = zero
406 e = zero
407 ENDIF
408C
409 IF ((ifail2 == 0).OR.(ifail2 ==1)) THEN
410 dn = dn * lscale / crit_scale(2)
411 dx = dx * lscale / crit_scale(2)
412 ENDIF
413 IF (dn == zero) dn=-ep30* crit_scale(2)
414 IF (dx == zero) dx= ep30* crit_scale(2)
415 IF ((ifail2 == 0).OR.(ifail2 ==1)) crit_scale(1) = lscale / crit_scale(2)
416C----
417 geo(10) =xk / a
418 geo(11) =xc
419 geo(45) =a
420 geo(46) =b
421 geo(47) =d
422 geo(180)=e
423 geo(133)= gf3
424 geo(48) =one/f
425 geo(174)=one/lscale
426 geo(67) =dn
427 geo(68) =dx
428C-- If H=6 - additional internal variables must be stored in UVAR - - GEO(25) = NUVAR
429 IF (iecrou == 6) THEN
430 geo(25) = 6
431 ENDIF
432C
433 igeo(105) =iecrou
434 igeo(106) =ifunc
435 igeo(107) =ifunc2
436 igeo(108) =ifv
437 igeo(120) = ifunc3
438C----
439 IF(.NOT.is_encrypted)THEN
440 IF (iecrou/=5) THEN
441 IF (ifail2==3) THEN
442 WRITE(iout,1813)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
443 . a,b,d,e,gf3,ifv,ifunc3,dx
444 ELSEIF (ifail2==2) THEN
445 WRITE(iout,1812)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
446 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
447 ELSE
448 WRITE(iout,1810)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
449 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
450 ENDIF
451 ELSE
452 IF (ifail2==3) THEN
453 WRITE(iout,1823)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
454 . a,b,d,e,gf3,ifv,ifunc3,dx
455 ELSEIF (ifail2==2) THEN
456 WRITE(iout,1822)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
457 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
458 ELSE
459 WRITE(iout,1820)'Y SHEAR',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
460 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
461 ENDIF
462 ENDIF
463 ENDIF
464C
465!-------------------------------------------------------
466! Rotations
467!-------------------------------------------------------
468C
469!-----------------
470 ! Torsion
471!-----------------
472C--------------------------------------------------
473C EXTRACT DATAS (INTEGER VALUES)
474C--------------------------------------------------
475 CALL hm_get_intv('FUN_A3',ifunc,is_available,lsubmodel)
476 CALL hm_get_intv('HFLAG3',iecrou,is_available,lsubmodel)
477 CALL hm_get_intv('FUN_B3',ifv,is_available,lsubmodel)
478 CALL hm_get_intv('FUN_C3',ifunc2,is_available,lsubmodel)
479 CALL hm_get_intv('FUN_D3',ifunc3,is_available,lsubmodel)
480C--------------------------------------------------
481C EXTRACT DATAS (REAL VALUES)
482C--------------------------------------------------
483 CALL hm_get_floatv('STIFF3',xk,is_available,lsubmodel,unitab)
484 CALL hm_get_floatv('DAMP3',xc,is_available,lsubmodel,unitab)
485 CALL hm_get_floatv('Acoeft3',a,is_available,lsubmodel,unitab)
486 CALL hm_get_floatv('Bcoeft3',b,is_available,lsubmodel,unitab)
487 CALL hm_get_floatv('Dcoeft3',d,is_available,lsubmodel,unitab)
488 CALL hm_get_floatv('MIN_RUP3',dn,is_available,lsubmodel,unitab)
489 CALL hm_get_floatv('MAX_RUP3',dx,is_available,lsubmodel,unitab)
490 CALL hm_get_floatv('Prop_Tor_F',f,is_available,lsubmodel,unitab)
491 CALL hm_get_floatv('Prop_Tor_E',e,is_available,lsubmodel,unitab)
492 CALL hm_get_floatv('scale3',lscale,is_available,lsubmodel,unitab)
493 CALL hm_get_floatv('F_EMC_b0',gf3,is_available,lsubmodel,unitab)
494 !units for default values
495 CALL hm_get_floatv_dim('Acoeft3',a_unit,is_available,lsubmodel,unitab)
496 CALL hm_get_floatv_dim('Dcoeft3',d_unit,is_available,lsubmodel,unitab)
497 CALL hm_get_floatv_dim('Prop_Tor_F',f_unit,is_available,lsubmodel,unitab)
498 CALL hm_get_floatv_dim('Prop_Tor_E',e_unit,is_available,lsubmodel,unitab)
499 CALL hm_get_floatv_dim('scale3',lscale_unit,is_available,lsubmodel,unitab)
500 CALL hm_get_floatv_dim('F_EMC_b0',gf3_unit,is_available,lsubmodel,unitab)
501 CALL hm_get_floatv_dim('MIN_RUP3',crit_scale(3),is_available,lsubmodel,unitab)
502C
503C----
504! IF(ifunc/=0.AND.iecrou>=1.AND.xk==zero)THEN
505! CALL ANCMSG(MSGID=230,
506! . MSGTYPE=MSGERROR,
507! . ANMODE=ANINFO_BLIND_1,
508! . I1=IG,
509! . C1=IDTITL)
510! ENDIF
511 IF(iecrou==4.AND.(ifunc==0.OR.ifunc2==0))THEN
512 CALL ancmsg(msgid=231,
513 . msgtype=msgerror,
514 . anmode=aninfo_blind_1,
515 . i1=ig,
516 . c1=idtitl)
517 ENDIF
518 IF(iecrou==4.AND.isk==zero)THEN
519 CALL ancmsg(msgid=230,
520 . msgtype=msgerror,
521 . anmode=aninfo_blind_1,
522 . i1=ig,
523 . c1=idtitl)
524 ENDIF
525 IF(iecrou==5.AND.(ifunc==0.OR.ifunc2==0))THEN
526 CALL ancmsg(msgid=231,
527 . msgtype=msgerror,
528 . anmode=aninfo_blind_1,
529 . i1=ig,
530 . c1=idtitl)
531 ENDIF
532 IF(iecrou==6.AND.(ifunc==0.OR.ifunc2==0))THEN
533 CALL ancmsg(msgid=1057,
534 . msgtype=msgerror,
535 . anmode=aninfo_blind_1,
536 . i1=ig,
537 . c1=idtitl)
538 ENDIF
539 IF(iecrou==7.AND.ifunc==0)THEN
540 CALL ancmsg(msgid=1058,
541 . msgtype=msgerror,
542 . anmode=aninfo_blind_1,
543 . i1=ig,
544 . c1=idtitl)
545
546 ELSEIF(iecrou==7.AND.ifunc2==0)THEN
547 CALL ancmsg(msgid=1059,
548 . msgtype=msgwarning,
549 . anmode=aninfo_blind_1,
550 . i1=ig,
551 . c1=idtitl,
552 . i2=iecrou)
553 iecrou = 2
554 ENDIF
555 IF (ifunc == 0 .AND. a /= zero .AND. a /= one) THEN
556 CALL ancmsg(msgid=663,
557 . msgtype=msgwarning,
558 . anmode=aninfo_blind_1,
559 . i1=ig,
560 . c1=idtitl)
561 ENDIF
562C----
563 IF (a == zero) a = one * a_unit
564 IF (d == zero) d = one * d_unit
565 IF (e == zero) e = one * e_unit
566 IF (f == zero) f = one * f_unit
567 IF (lscale == zero) lscale = one * lscale_unit
568 IF (gf3 == zero) gf3 = one * gf3_unit
569 IF (ifunc == 0) THEN
570 a = one
571 b = zero
572 e = zero
573 ENDIF
574C
575 IF ((ifail2 == 0).OR.(ifail2 ==1)) THEN
576 dn = dn * lscale / crit_scale(3)
577 dx = dx * lscale / crit_scale(3)
578 ENDIF
579 IF (dn == zero) dn=-ep30* crit_scale(3)
580 IF (dx == zero) dx= ep30* crit_scale(3)
581 IF ((ifail2 == 0).OR.(ifail2 ==1)) crit_scale(3) = lscale / crit_scale(3)
582C----
583 geo(19) =xk / a
584 geo(20) =xc
585 geo(53) =a
586 geo(54) =b
587 geo(55) =d
588 geo(182)=e
589 geo(135)= gf3
590 geo(56) =one/f
591 geo(176)=one/lscale
592 geo(71) =dn
593 geo(72) =dx
594C-- If H=6 - additional internal variables must be stored in UVAR - - GEO(25) = NUVAR
595 IF (iecrou == 6) THEN
596 geo(25) = 6
597 ENDIF
598C
599 igeo(109) =iecrou
600 igeo(110) =ifunc
601 igeo(111) =ifunc2
602 igeo(112) =ifv
603 igeo(121) = ifunc3
604C----
605 IF(.NOT.is_encrypted)THEN
606 IF (iecrou/=5) THEN
607 IF (ifail2==3) THEN
608 WRITE(iout,1833)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
609 . a,b,d,e,gf3,ifv,ifunc3,dx
610 ELSEIF (ifail2==2) THEN
611 WRITE(iout,1832)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
612 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
613 ELSE
614 WRITE(iout,1830)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
615 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
616 ENDIF
617 ELSE
618 IF (ifail2==3) THEN
619 WRITE(iout,1843)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
620 . a,b,d,e,gf3,ifv,ifunc3,dx
621 ELSEIF (ifail2==2) THEN
622 WRITE(iout,1842)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
623 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
624 ELSE
625 WRITE(iout,1840)'TORSION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
626 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
627 ENDIF
628 ENDIF
629 ENDIF
630C
631!-----------------
632 ! Rotation Y
633!-----------------
634C--------------------------------------------------
635C EXTRACT DATAS (INTEGER VALUES)
636C--------------------------------------------------
637 CALL hm_get_intv('FUN_A4',ifunc,is_available,lsubmodel)
638 CALL hm_get_intv('HFLAG4',iecrou,is_available,lsubmodel)
639 CALL hm_get_intv('FUN_B4',ifv,is_available,lsubmodel)
640 CALL hm_get_intv('FUN_C4',ifunc2,is_available,lsubmodel)
641 CALL hm_get_intv('FUN_D4',ifunc3,is_available,lsubmodel)
642C--------------------------------------------------
643C EXTRACT DATAS (REAL VALUES)
644C--------------------------------------------------
645 CALL hm_get_floatv('STIFF4',xk,is_available,lsubmodel,unitab)
646 CALL hm_get_floatv('DAMP4',xc,is_available,lsubmodel,unitab)
647 CALL hm_get_floatv('Acoeft4',a,is_available,lsubmodel,unitab)
648 CALL hm_get_floatv('Bcoeft4',b,is_available,lsubmodel,unitab)
649 CALL hm_get_floatv('Dcoeft4',d,is_available,lsubmodel,unitab)
650 CALL hm_get_floatv('MIN_RUP4',dn,is_available,lsubmodel,unitab)
651 CALL hm_get_floatv('MAX_RUP4',dx,is_available,lsubmodel,unitab)
652 CALL hm_get_floatv('Prop_Bend_F',f,is_available,lsubmodel,unitab)
653 CALL hm_get_floatv('Prop_Bend_E',e,is_available,lsubmodel,unitab)
654 CALL hm_get_floatv('scale4',lscale,is_available,lsubmodel,unitab)
655 CALL hm_get_floatv('F_EMC_c',gf3,is_available,lsubmodel,unitab)
656 !units for default values
657 CALL hm_get_floatv_dim('MIN_RUP4',crit_scale(4),is_available,lsubmodel,unitab)
658C----
659 IF(iecrou/=0 .AND. iecrou/=1 .AND. iecrou/=5) THEN
660 CALL ancmsg(msgid=513,
661 . msgtype=msgerror,
662 . anmode=aninfo_blind_1,
663 . i1=ig,
664 . c1=idtitl,i2=iecrou)
665 ENDIF
666! IF(IECROU>=1 .AND. IFUNC/=0 .AND. XK == ZERO)THEN
667! CALL ANCMSG(MSGID=230,
668! . MSGTYPE=MSGERROR,
669! . ANMODE=ANINFO_BLIND_1,
670! . I1=IG,
671! . C1=IDTITL)
672! ENDIF
673 IF(iecrou==5.AND.(ifunc==0.OR.ifunc2==0))THEN
674 CALL ancmsg(msgid=231,
675 . msgtype=msgerror,
676 . anmode=aninfo_blind_1,
677 . i1=ig,
678 . c1=idtitl)
679 ENDIF
680 IF(iecrou==6.AND.(ifunc==0.OR.ifunc2==0))THEN
681 CALL ancmsg(msgid=1057,
682 . msgtype=msgerror,
683 . anmode=aninfo_blind_1,
684 . i1=ig,
685 . c1=idtitl)
686 ENDIF
687 IF(iecrou==7.AND.ifunc==0)THEN
688 CALL ancmsg(msgid=1058,
689 . msgtype=msgerror,
690 . anmode=aninfo_blind_1,
691 . i1=ig,
692 . c1=idtitl)
693
694 ELSEIF(iecrou==7.AND.ifunc2==0)THEN
695 CALL ancmsg(msgid=1059,
696 . msgtype=msgwarning,
697 . anmode=aninfo_blind_1,
698 . i1=ig,
699 . c1=idtitl,
700 . i2=iecrou)
701 iecrou = 2
702 ENDIF
703 IF (ifunc == 0 .AND. a /= zero .AND. a /= one) THEN
704 CALL ancmsg(msgid=663,
705 . msgtype=msgwarning,
706 . anmode=aninfo_blind_1,
707 . i1=ig,
708 . c1=idtitl)
709 ENDIF
710C----
711 IF (a == zero) a = one * a_unit
712 IF (d == zero) d = one * d_unit
713 IF (e == zero) e = one * e_unit
714 IF (f == zero) f = one * f_unit
715 IF (lscale == zero) lscale = one * lscale_unit
716 IF (gf3 == zero) gf3 = one * gf3_unit
717 IF (ifunc == 0) THEN
718 a = one
719 b = zero
720 e = zero
721 ENDIF
722C
723 IF ((ifail2 == 0).OR.(ifail2 ==1)) THEN
724 dn = dn * lscale / crit_scale(4)
725 dx = dx * lscale / crit_scale(4)
726 ENDIF
727 IF (dn == zero) dn=-ep30* crit_scale(4)
728 IF (dx == zero) dx= ep30* crit_scale(4)
729 IF ((ifail2 == 0).OR.(ifail2 ==1)) crit_scale(4) = lscale / crit_scale(4)
730C----
731 geo(23) =xk / a
732 geo(24) =xc
733 geo(57) =a
734 geo(58) =b
735 geo(59) =d
736 geo(183)=e
737 geo(136)= gf3
738 geo(60) =one/f
739 geo(177)=one/lscale
740 geo(73)=dn
741 geo(74)=dx
742C-- If H=6 - additional internal variables must be stored in UVAR - - GEO(25) = NUVAR
743 IF (iecrou == 6) THEN
744 geo(25) = 6
745 ENDIF
746C
747 igeo(113) =iecrou
748 igeo(114) =ifunc
749 igeo(115) =ifunc2
750 igeo(116) =ifv
751 igeo(122) =ifunc3
752C----
753 IF(.NOT.is_encrypted)THEN
754 IF (iecrou/=5) THEN
755 IF (ifail2==3) THEN
756 WRITE(iout,1833)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
757 . a,b,d,e,gf3,ifv,ifunc3,dx
758 ELSEIF (ifail2==2) THEN
759 WRITE(iout,1832)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
760 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
761 ELSE
762 WRITE(iout,1830)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
763 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
764 ENDIF
765 ELSE
766 IF (ifail2==3) THEN
767 WRITE(iout,1843)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
768 . a,b,d,e,gf3,ifv,ifunc3,dx
769 ELSEIF (ifail2==2) THEN
770 WRITE(iout,1842)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
771 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
772 ELSE
773 WRITE(iout,1840)'FLEXION',xk,xc,ifunc,lscale,ifunc2,f,iecrou,
774 . a,b,d,e,gf3,ifv,ifunc3,dn,dx
775 ENDIF
776 ENDIF
777 ENDIF
778C
779!-------------------------------------------------------
780c--- New uni/multiaxial failure criteria
781!-------------------------------------------------------
782C
783 cc(1:4) = zero
784C
785C--------------------------------------------------
786C EXTRACT DATAS (REAL VALUES)
787C--------------------------------------------------
788 CALL hm_get_floatv('TRANS_VEL0',vt0,is_available,lsubmodel,unitab)
789 CALL hm_get_floatv('ROT_VEL0',vr0,is_available,lsubmodel,unitab)
790 CALL hm_get_floatv('C1',cc(1),is_available,lsubmodel,unitab)
791 CALL hm_get_floatv('REL_VEL_EXP1',cn(1),is_available,lsubmodel,unitab)
792 CALL hm_get_floatv('ALPHA1',xa(1),is_available,lsubmodel,unitab)
793 CALL hm_get_floatv('BETA1',xb(1),is_available,lsubmodel,unitab)
794 CALL hm_get_floatv('C2',cc(2),is_available,lsubmodel,unitab)
795 CALL hm_get_floatv('REL_VEL_EXP2',cn(2),is_available,lsubmodel,unitab)
796 CALL hm_get_floatv('ALPHA2',xa(2),is_available,lsubmodel,unitab)
797 CALL hm_get_floatv('BETA2',xb(2),is_available,lsubmodel,unitab)
798 CALL hm_get_floatv('C3',cc(3),is_available,lsubmodel,unitab)
799 CALL hm_get_floatv('REL_VEL_EXP3',cn(3),is_available,lsubmodel,unitab)
800 CALL hm_get_floatv('ALPHA3',xa(3),is_available,lsubmodel,unitab)
801 CALL hm_get_floatv('BETA3',xb(3),is_available,lsubmodel,unitab)
802 CALL hm_get_floatv('C4',cc(4),is_available,lsubmodel,unitab)
803 CALL hm_get_floatv('REL_VEL_EXP4',cn(4),is_available,lsubmodel,unitab)
804 CALL hm_get_floatv('ALPHA4',xa(4),is_available,lsubmodel,unitab)
805 CALL hm_get_floatv('BETA4',xb(4),is_available,lsubmodel,unitab)
806C
807 !units for default values
808 CALL hm_get_floatv_dim('TRANS_VEL0',vt0_unit,is_available, lsubmodel, unitab)
809 CALL hm_get_floatv_dim('ROT_VEL0',vr0_unit,is_available, lsubmodel, unitab)
810C-----
811C
812 IF (ifail2 > 0) THEN
813C
814 DO j = 1,4
815 IF (cn(j) == zero) cn(j) = one
816 IF (xa(j) == zero) xa(j) = one
817 IF (xb(j) == zero) xb(j) = two
818 ENDDO
819C----
820 IF ((ifail2 == 0).OR.(ifail2 == 1)) THEN
821 DO j = 1,4
822 cc(j) = cc(j) * crit_scale(j)
823 ENDDO
824 ENDIF
825C
826 ELSE
827C
828 DO j = 1,4
829 cc(j) = zero
830 cn(j) = one
831 xa(j) = one
832 xb(j) = two
833 ENDDO
834 ENDIF
835C
836 IF (vt0 == zero) vt0 = one * vt0_unit
837 IF (vr0 == zero) vr0 = one * vr0_unit
838C-----
839 IF (.NOT.is_encrypted)THEN
840 WRITE(iout,1849) 'RUPTURE',ifail,ifail2,vt0,vr0
841 IF (ifail == 0) THEN
842 WRITE(iout,1850) (cc(j),j=1,4),(cn(j),j=1,4)
843 ELSE
844 WRITE(iout,1851) (xa(j),j=1,4),(xb(j),j=1,4)
845 ENDIF
846 ENDIF
847C-----
848 geo(101) = vt0
849 geo(102) = vr0
850 geo(103) = cc(1)
851 geo(104) = cc(2)
852 geo(105) = cc(3)
853 geo(106) = cc(4)
854 geo(109) = cn(1)
855 geo(110) = cn(2)
856 geo(111) = cn(3)
857 geo(112) = cn(4)
858 geo(115) = xa(1)
859 geo(116) = xa(2)
860 geo(117) = xa(3)
861 geo(118) = xa(4)
862 geo(121) = xb(1)
863 geo(122) = xb(2)
864 geo(123) = xb(3)
865 geo(124) = xb(4)
866C
867C-----------------------------
868C PROPERTY BUFFER
869C-----------------------------
870C
871 prop_tag(igtyp)%G_FOR = 3
872 prop_tag(igtyp)%G_MOM = 3
873 prop_tag(igtyp)%G_LENGTH = 3
874 prop_tag(igtyp)%G_TOTDEPL = 3
875 prop_tag(igtyp)%G_TOTROT = 3
876 prop_tag(igtyp)%G_FOREP = 3
877 prop_tag(igtyp)%G_MOMEP = 3
878 prop_tag(igtyp)%G_DEP_IN_TENS = 3
879 prop_tag(igtyp)%G_DEP_IN_COMP = 3
880 prop_tag(igtyp)%G_ROT_IN_TENS = 3
881 prop_tag(igtyp)%G_ROT_IN_COMP = 3
882 prop_tag(igtyp)%G_POSX = 5
883 prop_tag(igtyp)%G_POSY = 5
884 prop_tag(igtyp)%G_POSZ = 5
885 prop_tag(igtyp)%G_POSXX = 5
886 prop_tag(igtyp)%G_POSYY = 5
887 prop_tag(igtyp)%G_POSZZ = 5
888 prop_tag(igtyp)%G_YIELD = 4
889 prop_tag(igtyp)%G_LENGTH_ERR = 3
890 prop_tag(igtyp)%G_SKEW = 3
891 prop_tag(igtyp)%G_SKEW_ERR = 3
892 prop_tag(igtyp)%G_E6 = 6
893 prop_tag(igtyp)%G_NUVAR = max(prop_tag(igtyp)%G_NUVAR,nint(geo(25))) ! additional internal variables for h=6
894 prop_tag(igtyp)%G_DEFINI = 6
895 prop_tag(igtyp)%G_FORINI = 6
896 prop_tag(igtyp)%G_RUPTCRIT = 1
897C
898 RETURN
899C-----------
900 1801 FORMAT(
901 & 5x,'SPRING PROPERTY SET (BEAM TYPE)'/,
902 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
903 & 5x,'SPRING MASS . . . . . . . . . . . . . .=',1pg20.13/,
904 & 5x,'SPRING INERTIA. . . . . . . . . . . . .=',1pg20.13/,
905 & 5x,'SKEW FRAME NUMBER (0:GLOBAL). . . . . .=',i10/,
906 & 5x,'SENSOR NUMBER (0:NOT USED). . . . . . .=',i10/,
907 & 5x,'SENSOR FLAG (0:ACTIV 1:DEACT 2:BOTH). .=',i10/,
908 & 5x,'UNIT LENGTH FLAG. . . . . . . . . . . .=',i10/,
909 & 5x,'IF=1 UNIT LENGTH MASS,STIFFNESS AND INPUT',/,
910 & 5x,' CURVE ARE STRAIN DEPENDING',/)
911C----
912 1810 FORMAT(
913 & 5x,a,/,
914 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
915 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
916 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
917 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
918 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
919 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
920 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
921 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
922 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
923 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
924 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
925 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
926 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
927 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
928 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
929 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
930 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
931 & 5x,'FUNCTION IDENTIFIER FOR ',/,
932 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
933 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
934 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
935 & 5x,'NEGATIVE FAILURE DISPLACEMENT . . . . .=',1pg20.13/,
936 & 5x,'POSITIVE FAILURE DISPLACEMENT . . . . .=',1pg20.13/)
937 1812 FORMAT(
938 & 5x,a,/,
939 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
940 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
941 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
942 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
943 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
944 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
945 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
946 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
947 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
948 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
949 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
950 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
951 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
952 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
953 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
954 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
955 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
956 & 5x,'FUNCTION IDENTIFIER FOR ',/,
957 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
958 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
959 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
960 & 5x,'NEGATIVE FAILURE FORCE. . . . . . . . .=',1pg20.13/,
961 & 5x,'POSITIVE FAILURE FORCE. . . . . . . . .=',1pg20.13/)
962 1813 FORMAT(
963 & 5x,a,/,
964 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
965 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
966 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
967 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
968 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
969 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
970 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
971 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
972 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
973 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
974 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
975 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
976 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
977 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
978 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
979 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
980 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
981 & 5x,'FUNCTION IDENTIFIER FOR ',/,
982 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
983 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
984 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
985 & 5x,'FAILURE ENERGY. . . . . . . . . . . . .=',1pg20.13/)
986 1820 FORMAT(
987 & 5x,a,/,
988 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
989 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
990 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
991 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
992 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
993 & 5x,'PERMANENT DISPL./MAX. DISPL. CURVE(H=5)=',i10/,
994 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
995 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
996 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
997 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
998 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
999 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1000 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1001 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1002 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1003 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1004 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1005 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1006 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1007 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1008 & 5x,'NEGATIVE FAILURE DISPLACEMENT . . . . .=',1pg20.13/,
1009 & 5x,'POSITIVE FAILURE DISPLACEMENT . . . . .=',1pg20.13/)
1010 1822 FORMAT(
1011 & 5x,a,/,
1012 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1013 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1014 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1015 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
1016 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1017 & 5x,'PERMANENT DISPL./MAX. DISPL. CURVE(H=5)=',i10/,
1018 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1019 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
1020 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
1021 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
1022 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
1023 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1024 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1025 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1026 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1027 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1028 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1029 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1030 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1031 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1032 & 5x,'NEGATIVE FAILURE FORCE. . . . . . . . .=',1pg20.13/,
1033 & 5x,'POSITIVE FAILURE FORCE. . . . . . . . .=',1pg20.13/)
1034 1823 FORMAT(
1035 & 5x,a,/,
1036 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1037 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1038 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1039 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
1040 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1041 & 5x,'PERMANENT DISPL./MAX. DISPL. CURVE(H=5)=',i10/,
1042 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1043 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
1044 & 5x,'0:ELASTIC 1:ISOTROPIC 2:UNCOUPLED',/,
1045 & 5x,'4:KINEMATIC 5:UNCOUPLED NL (UN/RE)LOADING',/,
1046 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
1047 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1048 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1049 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1050 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1051 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1052 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1053 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1054 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1055 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1056 & 5x,'FAILURE ENERGY. . . . . . . . . . . . .=',1pg20.13/)
1057C----
1058 1830 FORMAT(
1059 & 5x,a,/,
1060 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1061 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1062 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1063 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
1064 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1065 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
1066 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
1067 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1068 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
1069 & 5x,'0:ELASTIC 1:ISOTROPIC',/,
1070 & 5x,'5:UNCOUPLED NL (UN/RE)LOADING',/,
1071 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
1072 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1073 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1074 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1075 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1076 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1077 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1078 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1079 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1080 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1081 & 5x,'NEGATIVE FAILURE ROTATION . . . . . . .=',1pg20.13/,
1082 & 5x,'POSITIVE FAILURE ROTATION . . . . . . .=',1pg20.13/)
1083 1832 FORMAT(
1084 & 5x,a,/,
1085 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1086 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1087 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1088 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
1089 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1090 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
1091 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
1092 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1093 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
1094 & 5x,'0:ELASTIC 1:ISOTROPIC',/,
1095 & 5x,'5:UNCOUPLED NL (UN/RE)LOADING',/,
1096 & 5x,'6:ELASTO PLASTIC WITH HARDENING 7: ELASTIC HYSTERESIS',/,
1097 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1098 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1099 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1100 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1101 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1102 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1103 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1104 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1105 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1106 & 5x,'NEGATIVE FAILURE MOMENT . . . . . . . .=',1pg20.13/,
1107 & 5x,'POSITIVE FAILURE MOMENT . . . . . . . .=',1pg20.13/)
1108 1833 FORMAT(
1109 & 5x,a,/,
1110 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1111 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1112 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1113 & 5x,'FORCE-DISPLACEMENT CURVE. . . . . . . .=',i10/,
1114 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1115 & 5x,'FUNCTION IDENTIFIER FOR UNLOADING ',/,
1116 & 5x,'FORCE-DISPLACEMENT CURVE (H=4,5,7). . .=',i10/,
1117 & 5x,'ABSCISSA SCALE FACTOR ON CURVE . . . . =',1pg20.13/,
1118 & 5x,'HARDENING FLAG H. . . . . . . . . . . .=',i10/,
1119 & 5x,'0:ELASTIC 1:ISOTROPIC',/,
1120 & 5x,'5:UNCOUPLED NL (UN/RE)LOADING',/,
1121 & 5x,'6:elasto plastic with hardening 7: elastic hysteresis',/,
1122 & 5X,'dynamic amplification factor a. . . . .=',1PG20.13/,
1123 & 5X,'dynamic amplification factor b. . . . .=',1pg20.13/,
1124 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1125 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1126 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1127 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1128 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1129 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1130 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1131 & 5x,'FAILURE ENERGY. . . . . . . . . . . . .=',1pg20.13/)
1132 1840 FORMAT(
1133 & 5x,a,/,
1134 & 5x,'SPRING STIFFNESS. . . . . . . . . . . .=',1pg20.13/,
1135 & 5x,'SPRING DAMPING. . . . . . . . . . . . .=',1pg20.13/,
1136 & 5x,'FUNCTION IDENTIFIER FOR LOADING ',/,
1137 & 5x,'force-displacement curve. . . . . . . .=',I10/,
1138 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1139 & 5X,'permanent rot./max. rot. curve(h=5). .=',I10/,
1140 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1141 & 5X,'hardening flag h. . . . . . . . . . . .=',I10/,
1142 & 5X,'0:elastic 1:isotropic',/,
1143 & 5X,'5:uncoupled nl(un/re)loading',/,
1144 & 5X,'6:elasto plastic with hardening 7: elastic hysteresis',/,
1145 & 5X,'dynamic amplification factor a. . . . .=',1PG20.13/,
1146 & 5X,'dynamic amplification factor b. . . . .=',1PG20.13/,
1147 & 5X,'dynamic amplification factor d. . . . .=',1PG20.13/,
1148 & 5X,'dynamic amplification factor e. . . . .=',1PG20.13/,
1149 & 5X,'dynamic amplification factor gf3 . . .=',1PG20.13/,
1150 & 5X,'FUNCTION identifier for ',/,
1151 & 5X,'force-velocity curve. . . . . . . . . .=',I10/,
1152 & 5X,'function identifier for the additional ',/,
1153 & 5X,'force-velocity curve. . . . . . . . . .=',I10/,
1154 & 5X,'negative failure rotation . . . . . . .=',1PG20.13/,
1155 & 5X,'positive failure rotation . . . . . . .=',1PG20.13/)
1156 1842 FORMAT(
1157 & 5X,A,/,
1158 & 5X,'spring stiffness. . . . . . . . . . . .=',1PG20.13/,
1159 & 5X,'spring damping. . . . . . . . . . . . .=',1PG20.13/,
1160 & 5X,'function identifier for loading ',/,
1161 & 5X,'force-displacement curve. . . . . . . .=',I10/,
1162 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1163 & 5X,'permanent rot./max. rot. curve (H=5). .=',I10/,
1164 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1165 & 5X,'hardening flag h. . . . . . . . . . . .=',I10/,
1166 & 5X,'0:elastic 1:isotropic',/,
1167 & 5X,'5:uncoupled nl (UN/RE)loading',/,
1168 & 5X,'6:elasto plastic with hardening 7: elastic hysteresis',/,
1169 & 5X,'dynamic amplification factor a. . . . .=',1PG20.13/,
1170 & 5X,'dynamic amplification factor b. . . . .=',1PG20.13/,
1171 & 5X,'dynamic amplification factor d. . . . .=',1PG20.13/,
1172 & 5X,'dynamic amplification factor e. . . . .=',1PG20.13/,
1173 & 5X,'dynamic amplification factor gf3 . . .=',1PG20.13/,
1174 & 5X,'function identifier for ',/,
1175 & 5X,'force-velocity curve. . . . . . . . . .=',I10/,
1176 & 5X,'function identifier for the additional ',/,
1177 & 5X,'force-velocity curve. . . . . . . . . .=',I10/,
1178 & 5X,'negative failure moment . . . . . . . .=',1PG20.13/,
1179 & 5X,'positive failure moment . . . . . . . .=',1PG20.13/)
1180 1843 FORMAT(
1181 & 5X,A,/,
1182 & 5X,'spring stiffness. . . . . . . . . . . .=',1PG20.13/,
1183 & 5X,'spring damping. . . . . . . . . . . . .=',1PG20.13/,
1184 & 5X,'function identifier for loading ',/,
1185 & 5X,'force-displacement curve. . . . . . . .=',I10/,
1186 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1187 & 5X,'permanent rot./max. rot. curve (H=5). .=',I10/,
1188 & 5X,'abscissa scale factor on curve . . . . =',1PG20.13/,
1189 & 5X,'hardening flag h. . . . . . . . . . . .=',I10/,
1190 & 5X,'0:elastic 1:isotropic',/,
1191 & 5X,'5:uncoupled nl (UN/RE)loading',/,
1192 & 5X,'6:elasto plastic with hardening 7: elastic hysteresis',/,
1193 & 5x,'DYNAMIC AMPLIFICATION FACTOR A. . . . .=',1pg20.13/,
1194 & 5x,'DYNAMIC AMPLIFICATION FACTOR B. . . . .=',1pg20.13/,
1195 & 5x,'DYNAMIC AMPLIFICATION FACTOR D. . . . .=',1pg20.13/,
1196 & 5x,'DYNAMIC AMPLIFICATION FACTOR E. . . . .=',1pg20.13/,
1197 & 5x,'DYNAMIC AMPLIFICATION FACTOR GF3 . . .=',1pg20.13/,
1198 & 5x,'FUNCTION IDENTIFIER FOR ',/,
1199 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1200 & 5x,'FUNCTION IDENTIFIER FOR THE ADDITIONAL ',/,
1201 & 5x,'FORCE-VELOCITY CURVE. . . . . . . . . .=',i10/,
1202 & 5x,'FAILURE ENERGY. . . . . . . . . . . . .=',1pg20.13/)
1203C----
1204 1849 FORMAT(
1205 & 5x,a,/,
1206 & 5x,'FAILURE FLAG (0:UNCOUPLED 1:COUPLED). .=',i10/,
1207 & 5x,'FAILURE CRITERION (DISPL/FORCE/ENERGY).=',i10/,
1208 & 5x,' 1:DISPLACEMENT 2:FORCE 3:ENERGY ' ,/,
1209 & 5x,'TRANSL. REF. DEFORMATION VELOCITY . . .=',1pg20.13/,
1210 & 5x,'ROT. REF. DEFORMATION VELOCITY. . . . .=',1pg20.13)
1211 1850 FORMAT(
1212 & 5x,'C1 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1213 & 5x,'C2 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1214 & 5x,'C3 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1215 & 5x,'C4 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1216 & 5x,'N1 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1217 & 5x,'N2 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1218 & 5x,'N3 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1219 & 5x,'N4 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/)
1220 1851 FORMAT(
1221 & 5x,'A1 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1222 & 5x,'A2 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1223 & 5x,'A3 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1224 & 5x,'A4 COEFFICIENT. . . . . . . . . . . . .=',1pg20.13/,
1225 & 5x,'B1 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1226 & 5x,'B2 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1227 & 5x,'B3 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/,
1228 & 5x,'B4 EXPONENT . . . . . . . . . . . . . .=',1pg20.13/)
1229 RETURN
1230 END
my_real
#define my_real
Definition cppsort.cpp:32
damping
subroutine damping(nodft, nodlt, v, vr, a, ar, damp, ms, in, igrnod, dim, itask, weight, tagslv_rby, wfext)
Definition damping.F:882
the
end diagonal values have been computed in the(sparse) matrix id.SOL
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_read_prop25
subroutine hm_read_prop25(geo, igeo, unitab, iskn, igtyp, ig, prop_tag, idtitl, lsubmodel, sub_id)
Definition hm_read_prop25.F:40
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
elbuftag_mod
Definition elbuftag_mod.F:227
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
submodel_mod
Definition submodel_mod.F:712
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:720
unitab_mod
Definition unitab_mod.F:513
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
nl
character *2 function nl()
Definition message.F:2354
starter
program starter
Definition starter.F:39
elbuftag_mod::prop_tag_
Definition elbuftag_mod.F:374
submodel_mod::submodel_data
Definition submodel_mod.F:723