OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_prop06.F File Reference
#include "implicit_f.inc"
#include "units_c.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "scr17_c.inc"
#include "sphcom.inc"
#include "tablen_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_prop06 (geo, igeo, prop_tag, multi_fvm, igtyp, ig, idtitl, unitab, lsubmodel, rtrans, sub_id, iskn, ipart, sub_index, defaults_solid)

Function/Subroutine Documentation

◆ hm_read_prop06()

subroutine hm_read_prop06 ( dimension(npropg), intent(inout) geo,
integer, dimension(npropgi), intent(inout) igeo,
type(prop_tag_), dimension(0:maxprop), intent(inout) prop_tag,
type(multi_fvm_struct), intent(in) multi_fvm,
integer, intent(in) igtyp,
integer, intent(in) ig,
character(len=nchartitle), intent(in) idtitl,
type (unit_type_), intent(in) unitab,
type(submodel_data), dimension(nsubmod), intent(in) lsubmodel,
dimension(ntransf,nrtrans), intent(in) rtrans,
integer, intent(in) sub_id,
integer, dimension(liskn,siskwn/liskn), intent(in) iskn,
integer, dimension(lipart1,npart+nthpart), intent(in) ipart,
integer, intent(in) sub_index,
type(solid_defaults_), intent(in) defaults_solid )

Definition at line 41 of file hm_read_prop06.F.

44C-----------------------------------------------
45C ROUTINE DESCRIPTION :
46C ===================
47C READ PROPERTY TYPE06 WITH HM READER
48C-----------------------------------------------
49C DUMMY ARGUMENTS DESCRIPTION:
50C ===================
51C
52C NAME DESCRIPTION
53C
54C GEO PROPERTY ARRAY(REAL)
55C IGEO PROPERTY ARRAY(INTEGER)
56C PROP_TAG PROPERTY TAGS STRUCTURE
57C MULTI_FVM MULTI_FVM STRUCTURE
58C IGTYP PROPERTY TYPE
59C IG PROPERTY ID
60C IDTITL TITLES ARRAY
61C UNITAB UNITS ARRAY
62C LSUBMODEL SUBMODEL STRUCTURE
63C RTRANS TRANSFORMATION ARRAY(REAL)
64C SUB_ID SUBMODEL ID
65C ISKN SKEW ARRAY(INTEGER)
66C IPART PART ARRAY(INTEGER)
67C SUB_INDEX SUBMODEL INDEX
68C-----------------------------------------------
69C============================================================================
70C-----------------------------------------------
71C M o d u l e s
72C-----------------------------------------------
73 USE unitab_mod
74 USE elbuftag_mod
75 USE submodel_mod
76 USE message_mod
77 USE multi_fvm_mod
78 USE defaults_mod
79 USE names_and_titles_mod , ONLY : nchartitle
80C-----------------------------------------------
81C I m p l i c i t T y p e s
82C-----------------------------------------------
83#include "implicit_f.inc"
84C-----------------------------------------------
85C C o m m o n B l o c k s
86C-----------------------------------------------
87#include "units_c.inc"
88#include "com01_c.inc"
89#include "com04_c.inc"
90#include "param_c.inc"
91#include "scr17_c.inc"
92#include "sphcom.inc"
93#include "tablen_c.inc"
94#include "tabsiz_c.inc"
95C-----------------------------------------------
96C D u m m y A r g u m e n t s
97C-----------------------------------------------
98 INTEGER,INTENT(INOUT) :: IGEO(NPROPGI)
99 INTEGER,INTENT(IN) :: IG
100 INTEGER,INTENT(IN) :: ISKN(LISKN,SISKWN/LISKN)
101 INTEGER,INTENT(IN) :: IPART(LIPART1,NPART+NTHPART)
102 INTEGER,INTENT(IN) :: IGTYP
103 INTEGER,INTENT(IN) :: SUB_ID
104 INTEGER,INTENT(IN) :: SUB_INDEX
105C
106 my_real,INTENT(INOUT) :: geo(npropg)
107 my_real,INTENT(IN) :: rtrans(ntransf,nrtrans)
108C
109 CHARACTER(LEN=NCHARTITLE),INTENT(IN) :: IDTITL
110C
111 TYPE(PROP_TAG_) , DIMENSION(0:MAXPROP),INTENT(INOUT) :: PROP_TAG
112 TYPE(MULTI_FVM_STRUCT),INTENT(IN) :: MULTI_FVM
113 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(NSUBMOD)
114 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
115 TYPE(SOLID_DEFAULTS_), INTENT(IN) :: DEFAULTS_SOLID
116C-----------------------------------------------
117C L o c a l V a r i a b l e s
118C-----------------------------------------------
119 INTEGER IHBE,ISMSTR,IPLAS,ICPRE,ICSTR,IINT,IP,ISK,IREP,JCVT,
120 . NPT,NPTR,NPTS,NPTT,ITET4,IET,ISTR,IHBE_OLD,ITET10,
121 . ISHEAR,ISORTH,ICONTROL_D,ICONTROL
122 my_real
123 . cvis,qa,qb,qh,vx,vy,vz,angle,dtmin,px,py,pz,vn,
124 . vdefmin,vdefmax,aspmax,asptet
125 INTEGER K,NSPHDIR,ID_PARTSPH,IPARTSPH,J
126 INTEGER IHBE_DS,ISST_DS,IPLA_DS,IFRAME_DS, IMAS_DS,
127 . ITET4_D,ITET10_D,ICPRE_D
128 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
129C-----------------------------------------------
130 is_encrypted = .false.
131 is_available = .false.
132
133 isorth=1
134 ismstr=0
135 ihbe=0
136!--- defaults values
137 ihbe_ds= defaults_solid%ISOLID
138 isst_ds= defaults_solid%ISMSTR
139 icpre_d= defaults_solid%ICPRE
140 itet4_d= defaults_solid%ITETRA4
141 itet10_d= defaults_solid%ITETRA10
142 iframe_ds= defaults_solid%IFRAME
143 icontrol_d=defaults_solid%ICONTROL
144C======================================================================|
145C
146 geo(3)=ismstr
147
148 igeo( 1)=ig
149 igeo(10)=ihbe
150 igeo(11)=igtyp
151 geo(12) =igtyp+em01
152 geo(171)=ihbe
153C
154 icstr = 0
155C--------------------------------------------------
156C EXTRACT DATA (IS OPTION CRYPTED)
157C--------------------------------------------------
158 CALL hm_option_is_encrypted(is_encrypted)
159C--------------------------------------------------
160C EXTRACT DATAS (INTEGER VALUES)
161C--------------------------------------------------
162 CALL hm_get_intv('ISOLID',ihbe,is_available,lsubmodel)
163 CALL hm_get_intv('Ismstr',ismstr,is_available,lsubmodel)
164C
165C IPLAS=0 ! New HM reader <=> possible numerical change for decks translated from old formats (IPLAS is no more read since format 2018)
166c CALL HM_GET_INTV('Iplas',IPLAS,IS_AVAILABLE,LSUBMODEL)
167C
168 CALL hm_get_intv('Itetra10',itet10,is_available,lsubmodel)
169 CALL hm_get_intv('Itetra4',itet4,is_available,lsubmodel)
170 CALL hm_get_intv('Icpre',icpre,is_available,lsubmodel)
171 CALL hm_get_intv('Inpts_R',nptr,is_available,lsubmodel)
172 CALL hm_get_intv('Inpts_S',npts,is_available,lsubmodel)
173 CALL hm_get_intv('Inpts_T',nptt,is_available,lsubmodel)
174 CALL hm_get_intv('SKEW_CSID',isk,is_available,lsubmodel)
175 CALL hm_get_intv('REFPLANE',ip,is_available,lsubmodel)
176 CALL hm_get_intv('Iframe',jcvt,is_available,lsubmodel)
177 CALL hm_get_intv('ORTHTROP',irep,is_available,lsubmodel)
178c CALL HM_GET_INTV('Istrain',ISTR,IS_AVAILABLE,LSUBMODEL)
179c CALL HM_GET_INTV('IHKT',IET,IS_AVAILABLE,LSUBMODEL)
180 CALL hm_get_intv('Ndir',nsphdir,is_available,lsubmodel)
181 CALL hm_get_intv('SPHPART_ID',id_partsph,is_available,lsubmodel)
182 CALL hm_get_intv('Icontrol',icontrol,is_available,lsubmodel)
183C--------------------------------------------------
184C EXTRACT DATAS (REAL VALUES)
185C--------------------------------------------------
186 CALL hm_get_floatv('dn',cvis,is_available,lsubmodel,unitab)
187 CALL hm_get_floatv('Vx',vx,is_available,lsubmodel,unitab)
188 CALL hm_get_floatv('Vy',vy,is_available,lsubmodel,unitab)
189 CALL hm_get_floatv('Vz',vz,is_available,lsubmodel,unitab)
190 CALL hm_get_floatv('MAT_BETA',angle,is_available,lsubmodel,unitab)
191 CALL hm_get_floatv('qa',qa,is_available,lsubmodel,unitab)
192 CALL hm_get_floatv('qb',qb,is_available,lsubmodel,unitab)
193 CALL hm_get_floatv('h',qh,is_available,lsubmodel,unitab)
194 CALL hm_get_floatv('deltaT_min',dtmin,is_available,lsubmodel,unitab)
195 CALL hm_get_floatv('vdef_min',vdefmin,is_available,lsubmodel,unitab)
196 CALL hm_get_floatv('vdef_max',vdefmax,is_available,lsubmodel,unitab)
197 CALL hm_get_floatv('ASP_max',aspmax,is_available,lsubmodel,unitab)
198 CALL hm_get_floatv('COL_min',asptet,is_available,lsubmodel,unitab)
199 CALL hm_get_floatv('Px',px,is_available,lsubmodel,unitab)
200 CALL hm_get_floatv('Py',py,is_available,lsubmodel,unitab)
201 CALL hm_get_floatv('Pz',pz,is_available,lsubmodel,unitab)
202C--------------------------------------------------
203 IF(isk == 0 .AND. sub_index /= 0 ) isk = lsubmodel(sub_index)%SKEW
204C-----removed flags:
205 istr=1
206 iplas = 2
207 iet = 0
208C----------------------
209 IF(itet4 == 0) itet4 = itet4_d
210 IF(itet10 == 0) itet10 = itet10_d
211C
212C
213 IF (sub_id /= 0)
214 . CALL subrotvect(vx,vy,vz,rtrans,sub_id,lsubmodel)
215c
216 ipartsph=0
217 IF(id_partsph > 0) THEN
218 DO j=1,npart
219 IF(ipart(4,j)==id_partsph) THEN
220 ipartsph=j
221 GOTO 175
222 ENDIF
223 ENDDO
224 CALL ancmsg(msgid=1037,
225 . msgtype=msgerror,
226 . anmode=aninfo,
227 . i1=ig,
228 . c1=idtitl,
229 . i2=id_partsph)
230 CALL arret(2)
231175 CONTINUE
232 END IF
233C
234C---------------------------------------------------
235C--- Default values
236C
237 iint = 1
238C------ new solid17, IHBE will be changed to 17 in sgrtails.F
239 IF (ihbe == 18 ) iint = 2
240 cvis = zero
241C
242C ihbe
243C
244 IF (ihbe == 0) ihbe = ihbe_ds
245C
246C old formats < 5.1 (BTW these formats are no more supported)
247 IF(ihbe == 101) THEN
248 ihbe=1
249 jcvt=2
250 END IF
251 IF(ihbe == 102) THEN
252 ihbe=2
253 jcvt=2
254 END IF
255 IF(ihbe == 104) THEN
256 ihbe=24
257 jcvt=2
258 END IF
259 IF(ihbe == 112) THEN
260 ihbe=12
261 jcvt=2
262 END IF
263C------- used for elasto-platic critia parameter no necessary with IET---
264 IF(n2d>0 .AND. ihbe/=0 .AND. ihbe/=2.AND. ihbe/=17)THEN
265 ihbe_old=ihbe
266 ihbe=0
267 CALL ancmsg(msgid=321,
268 . msgtype=msgwarning,
269 . anmode=aninfo_blind_2,
270 . i1=ig,
271 . c1=idtitl,
272 . i2=ihbe_old,
273 . i3=ihbe)
274 ELSEIF (ihbe/= 1.AND.ihbe/= 2.AND.ihbe/=14.AND.ihbe/= 24
275 . .AND.ihbe/= 17.AND.ihbe/= 18) THEN
276 CALL ancmsg(msgid=533,
277 . msgtype=msgwarning,
278 . anmode=aninfo_blind_1,
279 . i1=ig,
280 . c1=idtitl,
281 . i2=ihbe)
282 ihbe=1
283 ENDIF
284C
285C jcvt
286C
287 IF (jcvt == 0) jcvt = iframe_ds
288 IF (ihbe == 14.OR.ihbe == 18) jcvt = 2
289 IF (ihbe == 24) jcvt = 2
290 IF (iframe_ds == -2.OR.jcvt<0) jcvt = -1
291C
292C smstr
293C
294 IF(ismstr == 0) ismstr=isst_ds
295 IF(ismstr == 0.AND.ihbe /= 18) ismstr=4
296 IF (isst_ds == -2) ismstr = -1
297C
298C IPLAS
299C
300c IF (IPLAS == 0)IPLAS=IPLA_DS
301c IF (IPLAS /=1 .AND. IPLAS /=2 .AND. IPLAS /=3 ) THEN
302c CALL ANCMSG(MSGID=416,
303c . MSGTYPE=MSGERROR,
304c . ANMODE=ANINFO_BLIND_1,
305c . I1=IG,
306c . C1=IDTITL,
307c . I2=IPLAS)
308c ENDIF
309c IF (IPLAS == 1)THEN
310c CALL ANCMSG(MSGID=417,
311c . MSGTYPE=MSGWARNING,
312c . ANMODE=ANINFO_BLIND_1,
313c . I1=IG,
314c . C1=IDTITL,
315c . I2=IPLAS)
316c ENDIF
317C
318C icpre/icstr
319C
320 IF (icpre == 0) icpre = icpre_d
321 IF((n2d > 0 .AND. ihbe == 17) ) THEN
322 IF(icpre/=1 .AND. icpre/=2) icpre=0
323C------ no effet for Axi Isolid17 for the moment
324 IF(n2d == 1 .AND. ihbe == 17) icpre=0
325 ELSE
326 IF (ihbe /= 14 .AND. ihbe /= 24 .AND. ihbe /= 17 .AND. ihbe /= 18) icpre = 0
327 IF (ihbe == 17 ) THEN
328 IF (icpre == 0 ) THEN
329 icpre = 1
330 ELSEIF(icpre == 3 ) THEN
331 icpre = 0
332 ENDIF
333 ENDIF
334 IF (icpre == 3 .AND. ihbe /= 18) icpre =0
335 END if!((N2D > 0 .AND. IHBE == 17) .OR.
336 icstr = 0
337 IF (icpre_d == -2) icpre = -1
338C
339C NPT
340C
341 SELECT CASE (ihbe)
342 CASE(14,16)
343 IF(nptr == 0) nptr = 2
344 IF(npts == 0) npts = 2
345 IF(nptt == 0) nptt = 2
346 npt=nptr*100+npts*10+nptt
347 IF (ihbe == 14 .AND.
348 . (nptr < 1 .OR. npts < 1 .OR. nptt < 1 .OR.
349 . nptr > 9 .OR. npts > 9 .OR. nptt > 9)) THEN
350 CALL ancmsg(msgid=563,
351 . msgtype=msgerror,
352 . anmode=aninfo_blind_1,
353 . i1=ig,
354 . c1=idtitl,
355 . i2=npt,
356 . i3=ihbe)
357 ELSEIF (ihbe == 16 .AND.
358 . (nptr < 1 .OR. npts < 1 .OR. nptt < 1 .OR.
359 . nptr > 3 .OR. npts > 9 .OR. nptt > 3)) THEN
360 CALL ancmsg(msgid=563,
361 . msgtype=msgerror,
362 . anmode=aninfo_blind_1,
363 . i1=ig,
364 . c1=idtitl,
365 . i2=npt,
366 . i3=ihbe)
367 ENDIF
368 CASE(1,2,24)
369 npt = 1
370 CASE(12,13,17,18)
371 npt = 8
372 END SELECT
373C
374 IF((n2d > 0 .AND. ihbe == 17) .OR.
375 . (n2d == 1 .AND. ihbe == 22)) THEN
376 npt = 4
377 ENDIF
378C
379C viscosity / hourglass
380C
381C igeo(31) flag for default qa qb for law 70 can be used for other law
382 IF(qa == zero .AND. qb == zero) igeo(31) = 1
383 IF (qa == zero) qa = onep1
384 IF (qb == zero) qb = fiveem2
385 IF (qh == zero) qh = em01
386C
387C hourglass
388C
389 IF (ihbe == 24) THEN
390 IF (cvis == zero) cvis = em01
391 geo(13) = cvis
392 qh = zero
393 iint = iet
394 ELSEIF (ihbe==1.OR.ihbe==2) THEN
395 geo(13) = qh
396 ELSE
397 qh = zero
398 geo(13) = zero
399 ENDIF
400C
401C orthotropy
402C
403C IP =20, 21:Pj;23:Vj+MAT_BETA;24:Vj+Pj
404 IF (ip == 23 .OR. ip == 24) THEN
405C---- error-out if Vj=0+IP=23,24
406 vn = vx*vx+vy*vy+vz*vz
407 IF (vn<em20) THEN
408 CALL ancmsg(msgid=1918,
409 . msgtype=msgerror,
410 . anmode=aninfo_blind_1,
411 . i1=ig,
412 . c1=idtitl,
413 . i2=ip)
414 ENDIF
415 END IF
416 IF (ip == 0) THEN
417 DO k=0,numskw+min(1,nspcond)*numsph+nsubmod
418 IF(isk == iskn(4,k+1)) THEN
419 ip=-(k+1)
420 isk=k+1
421 GO TO 100
422 ENDIF
423 ENDDO
424 CALL ancmsg(msgid=137,anmode=aninfo,msgtype=msgerror,
425 . c1='PROPERTY',
426 . c2='PROPERTY',
427 . i1=ig,i2=isk,c3=idtitl)
428100 CONTINUE
429 ENDIF
430C Istrain
431c IF (ISTR == 0)ISTR =ISTRA_D
432 IF (icontrol==0) icontrol=icontrol_d
433 IF (icontrol>1) icontrol=0
434C----
435 igeo(2) = ip
436 igeo(4) = npt
437 igeo(5) = ismstr
438 igeo(6) = irep
439 igeo(9) = iplas-1
440 igeo(10) = ihbe
441 igeo(12) = istr
442 igeo(13) = icpre
443 igeo(14) = icstr
444 igeo(15) = iint
445 igeo(16) = jcvt-1
446 igeo(37) = nsphdir
447 igeo(38) = ipartsph
448 igeo(97) = icontrol
449C
450 geo(1) = angle
451 geo(7) = vx
452 geo(8) = vy
453 geo(9) = vz
454 geo(14) = qa
455 geo(15) = qb
456 geo(172)= dtmin
457 geo(33) = px
458 geo(34) = py
459 geo(35) = pz
460 geo(190)= vdefmin
461 geo(191)= vdefmax
462 geo(192)= aspmax
463 geo(193)= asptet
464C----
465 IF(.NOT. is_encrypted)THEN
466 IF(igeo(31) == 1)THEN
467 WRITE(iout,1100)ig,ihbe,ismstr,iplas,npt,jcvt,itet4,itet10,
468 . icpre,icstr,cvis,qa,qb,qh,dtmin,irep,istr,icontrol
469 ELSE
470 WRITE(iout,1000)ig,ihbe,ismstr,iplas,npt,jcvt,itet4,itet10,
471 . icpre,icstr,cvis,qa,qb,qh,dtmin,irep,istr,icontrol
472 ENDIF
473 IF((vdefmin+vdefmax+aspmax+asptet)>zero) THEN
474 IF (vdefmax==zero) vdefmax=ep10
475 IF (aspmax==zero) aspmax=ep10
476 WRITE(iout,3000) vdefmin,vdefmax,aspmax,asptet
477 END IF
478 IF(ip < 0) THEN
479 WRITE(iout,1001) iskn(4,isk)
480 ELSEIF(ip == 1 .OR. ip == 2 .OR. ip == 3) THEN
481 WRITE(iout,1002) ip,angle
482 ELSEIF(ip == 11 .OR. ip == 12 .OR. ip == 13) THEN
483 WRITE(iout,1003) ip,vx,vy,vz
484 ELSEIF(ip == 20) THEN
485 WRITE(iout,2001) ip
486 ELSEIF(ip == 21) THEN
487 WRITE(iout,2002) ip,px,py,pz
488 ELSEIF(ip == 23) THEN
489 WRITE(iout,2003) ip,angle,vx,vy,vz
490 ELSEIF(ip == 24) THEN
491 WRITE(iout,2004) ip,px,py,pz,vx,vy,vz
492 ENDIF
493 IF (iet > 0) WRITE(iout,2010) iet
494 IF(nsphdir/=0)WRITE(iout,2020)nsphdir, id_partsph
495 ELSE
496 WRITE(iout,1099) ig
497 ENDIF
498
499 IF (itet4 == 1000) itet4 = 0
500 igeo(20) = itet4
501 IF (itet10 == 1000) itet10 = 0
502 igeo(50) = itet10
503C
504 prop_tag(igtyp)%G_SIG = 6
505 prop_tag(igtyp)%L_SIG = 6
506 prop_tag(igtyp)%G_EINT = 1
507 prop_tag(igtyp)%G_QVIS = 1
508 prop_tag(igtyp)%L_EINT = 1
509 prop_tag(igtyp)%G_VOL = 1
510 prop_tag(igtyp)%L_VOL = 1
511 prop_tag(igtyp)%L_QVIS = 1
512 IF (multi_fvm%IS_USED) THEN
513 prop_tag(igtyp)%G_MOM = 3
514 ENDIF
515C
516 prop_tag(igtyp)%G_FILL = 1
517 prop_tag(igtyp)%L_STRA = 6
518 IF (n2d /= 0 .AND. multi_fvm%IS_USED) THEN
519! 2D analysis
520 prop_tag(igtyp)%G_AREA = 1
521 ENDIF
522 prop_tag(igtyp)%G_GAMA = 6
523 prop_tag(igtyp)%L_SIGL = 6
524 IF (geo(16) /= zero .OR. geo(17) /= zero) THEN
525 igeo(33) = 1 ! ISVIS flag
526 ENDIF
527
528C-------------------------------
529 igeo(1) =ig
530 igeo(11)=igtyp
531 igeo(17)=isorth
532 IF(geo( 3)/=zero.AND.igeo( 5)== 0) igeo( 5)=nint(geo( 3))
533 IF(geo(39)/=zero.AND.igeo( 9)== 0) igeo( 9)=nint(geo(39))
534 IF(geo(171)/=zero.AND.igeo(10)== 0) igeo(10)=nint(geo(171))
535C
536C----
537 RETURN
538C---
539 1000 FORMAT(
540 & 5x,'ORTHOTROPIC SOLID PROPERTY SET'/,
541 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
542 & 5x,'SOLID FORMULATION FLAG. . . . . . . . .=',i10/,
543 & 5x,'SMALL STRAIN FLAG . . . . . . . . . . .=',i10/,
544 & 5x,'SOLID STRESS PLASTICITY FLAG. . . . . .=',i10/,
545 & 5x,'NUMBER OF INTEGRATION POINTS. . . . .=',i10/,
546 & 5x,'COROTATIONAL SYSTEM FLAG. . . . . . . .=',i10/,
547 & 5x,'TETRA4 FORMULATION FLAG. . . . . . . .=',i10/,
548 & 5x,'TETRA10 FORMULATION FLAG . . . . . . .=',i10/,
549 & 5x,'CONSTANT PRESSURE FLAG. . . . . . . . .=',i10/,
550 & 5x,'CONSTANT STRESS FLAG. . . . . . . . . .=',i10/,
551 & 5x,'HOURGLASS NUMERICAL DAMPING . . . . . .=',1pg20.13/,
552 & 5x,'QUADRATIC BULK VISCOSITY. . . . . . . .=',1pg20.13/,
553 & 5x,'LINEAR BULK VISCOSITY . . . . . . . . .=',1pg20.13/,
554 & 5x,'HOURGLASS VISCOSITY . . . . . . . . . .=',1pg20.13/,
555 & 5x,'BRICK MINIMUM TIME STEP................=',1pg20.13/,
556 & 5x,'LOCAL ORTHOTROPY SYSTEM FORMULATION . .=',i10/,
557 & 5x,'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
558 & 5x,'SOLID DISTORTION CONTROL FLAG . . . . .=',i10/)
559 1099 FORMAT(
560 & 5x,'ORTHOTROPIC SOLID PROPERTY SET'/,
561 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i8/,
562 & 5x,'CONFIDENTIAL DATA'//)
563 1100 FORMAT(
564 & 5x,'ORTHOTROPIC SOLID PROPERTY SET'/,
565 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
566 & 5x,'SOLID FORMULATION FLAG. . . . . . . . .=',i10/,
567 & 5x,'SMALL STRAIN FLAG . . . . . . . . . . .=',i10/,
568 & 5x,'SOLID STRESS PLASTICITY FLAG. . . . . .=',i10/,
569 & 5x,'NUMBER OF INTEGRATION POINTS. . . . .=',i10/,
570 & 5x,'COROTATIONAL SYSTEM FLAG. . . . . . . .=',i10/,
571 & 5x,'tetra4 formulation flag. . . . . . . .=',I10/,
572 & 5X,'tetra10 formulation flag . . . . . . .=',I10/,
573 & 5X,'constant pressure flag. . . . . . . . .=',i10/,
574 & 5x,'CONSTANT STRESS FLAG. . . . . . . . . .=',i10/,
575 & 5x,'HOURGLASS NUMERICAL DAMPING . . . . . .=',1pg20.13/,
576 & 5x,'DEFAULT VALUE FOR QUADRATIC BULK. . . . ',/,
577 & 5x,' VISCOSITY (QA) WILL BE USED. . . .=',1pg20.13/,
578 & 5x,'EXCEPT IN CASE LAW 70 WHERE QA = 0. ',/,
579 & 5x,'DEFAULT VALUE FOR LINEAR BULK . . . . . ',/,
580 & 5x,' VISCOSITY (QB) WILL BE USED . . . =',1pg20.13/,
581 & 5x,'EXCEPT IN CASE LAW 70 WHERE QB = 0. ',/,
582 & 5x,'HOURGLASS VISCOSITY . . . . . . . . . .=',1pg20.13/,
583 & 5x,'BRICK MINIMUM TIME STEP................=',1pg20.13/,
584 & 5x,'LOCAL ORTHOTROPY SYSTEM FORMULATION . .=',i10/,
585 & 5x,'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
586 & 5x,'SOLID DISTORTION CONTROL FLAG . . . . .=',i10/)
587 1001 FORMAT(
588 & 5x,'ORTHOTROPIC SKEW FRAME. . . . . . . . .=',i10)
589 1002 FORMAT(
590 & 5x,'ORTHOTROPIC PLANE NUMBER. . . . . . . .=',i10/,
591 & 5x,' 1=(R,S) 2=(S,T) 3=(T,R) ',/,
592 & 5x,'ORTHOTROPIC ANGLE . . . . . . . . . . .=',1pg20.13)
593 1003 FORMAT(
594 & 5x,'ORTHOTROPIC PLANE NUMBER. . . . . . . .=',i10/,
595 & 5x,' 1=(R,S) 2=(S,T) 3=(T,R) ',/,
596 & 5x,'REFERENCE VECTOR VX . . . . . . . . . .=',1pg20.13/,
597 & 5x,'REFERENCE VECTOR VY . . . . . . . . . .=',1pg20.13/,
598 & 5x,'REFERENCE VECTOR VZ . . . . . . . . . .=',1pg20.13)
599 2010 FORMAT(
600 & 5x,'HOURGLASS MODULUS FLAG. . . . . . . . .=',i10/)
601 2020 FORMAT(
602 & 5x,'NUMBER OF SPH PARTICLES PER DIRECTION .=',i10/,
603 & 5x,'CORRESPONDING PART FOR SPH PARTICLES. .=',i10/)
604 2001 FORMAT(
605 & 5x,'ORTHOTROPIC DIRECTIONS BY ELEMENT CONNECTIVITY,IP='i10)
606 2002 FORMAT(
607 & 5x,'ORTHOTROPIC DIRECTION FLAG IP. . . . . =',i10/,
608 & 5x,'REFERENCE POINT PX . . . . . . . . . . =',1pg20.13/,
609 & 5x,'REFERENCE POINT PY . . . . . . . . . . =',1pg20.13/,
610 & 5x,'REFERENCE POINT PZ . . . . . . . . . . =',1pg20.13)
611 2003 FORMAT(
612 & 5x,'orthotropic direction flag ip . . . . .=',I10/,
613 & 5X,'orthotropic angle . . . . . . . . . . .=',1PG20.13/,
614 & 5X,'reference vector vx . . . . . . . . . .=',1PG20.13/,
615 & 5X,'reference vector vy . . . . . . . . . .=',1PG20.13/,
616 & 5X,'reference vector vz . . . . . . . . . .=',1PG20.13)
617 2004 FORMAT(
618 & 5X,'orthotropic direction flag ip. . . . . =',I10/,
619 & 5X,'reference point px . . . . . . . . . . =',1PG20.13/,
620 & 5X,'reference point py . . . . . . . . . . =',1PG20.13/,
621 & 5X,'reference point pz . . . . . . . . . . =',1PG20.13/,
622 & 5X,'reference vector vx . . . . . . . . . .=',1PG20.13/,
623 & 5X,'reference vector vy . . . . . . . . . .=',1PG20.13/,
624 & 5X,'reference vector vz . . . . . . . . . .=',1PG20.13)
625 3000 FORMAT(
626 & 5X,'solid minimum volumetric strain........=',1PG20.13/,
627 & 5X,'solid maximum volumetric strain........=',1PG20.13/,
628 & 5X,'solid maximum aspect ratio.............=',1PG20.13/,
629 & 5X,'solid minimum collapse ratio...........=',1PG20.13/)
630C---
my_real
#define my_real
Definition cppsort.cpp:32
if
if(complex_arithmetic) id
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
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
min
#define min(a, b)
Definition macros.h:20
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
arret
subroutine arret(nn)
Definition arret.F:87
istr
integer function istr(str)
Definition istr.F:31
subrotvect
subroutine subrotvect(x, y, z, rtrans, sub_id, lsubmodel)
Definition subrot.F:54