41 . RTRANS,LSUBMODEL,IDTITL ,IG ,SUB_ID ,ISKN ,
55#include "implicit_f.inc"
65#include "tablen_c.inc"
70 TYPE (UNIT_TYPE_),
INTENT(IN) ::UNITAB
72 . igtyp , igeo(*),iskn(liskn,*), ig
74 my_real GEO(*),RTRANS(NTRANSF,*)
76 CHARACTER(LEN=NCHARTITLE)::IDTITL
78 TYPE(
prop_tag_) ,
DIMENSION(0:MAXPROP) :: PROP_TAG
79 TYPE(SHELL_DEFAULTS_),
INTENT(IN) :: DEFAULTS_SHELL
83 INTEGER I, ISMSTR, NIP, J,
84 . ISHEAR, IP, ISTRAIN,
85 . ihbe,iplast,ithk,ibid,idf,ihbeoutp,k,n,
87 . pid1,ipid1, ihgflu, ihbe_old,nstack,igmat
88 INTEGER FLAG_FMT,FLAG_FMT_TMP
89 INTEGER ISH3N,ISROT,SUB_ID,IRP,IDSK,ISK,IUN,IPOS
90 INTEGER IHBE_D,IPLA_D,ISTR_D,ITHK_D,ISHEA_D,ISST_D,
91 . ish3n_d, istra_d,npts_d,idril_d
94 . angl,pun,cvis,rbid,vx,vy,vz,
95 . pthk, an, phi,zshift
96 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
97 CHARACTER(LEN=NCHARTITLE) :: TITR
117 is_encrypted = .false.
118 is_available = .false.
120 ihbe_d = defaults_shell%ishell
121 ish3n_d= defaults_shell%ish3n
122 isst_d = defaults_shell%ismstr
123 ipla_d = defaults_shell%iplas
124 ithk_d = defaults_shell%ithick
125 idril_d= defaults_shell%idrill
144 CALL hm_get_intv(
'Ishell',ihbe,is_available,lsubmodel)
145 CALL hm_get_intv(
'Ismstr',ismstr,is_available,lsubmodel)
146 CALL hm_get_intv(
'ISH3N',ish3n,is_available,lsubmodel)
147 CALL hm_get_intv('idrill
',ISROT,IS_AVAILABLE,LSUBMODEL)
148 CALL HM_GET_INTV('nip
',NIP,IS_AVAILABLE,LSUBMODEL)
150 CALL HM_GET_INTV('ithick
',ITHK,IS_AVAILABLE,LSUBMODEL)
151 CALL HM_GET_INTV('iplas',iplast,is_available,lsubmodel)
152 CALL hm_get_intv(
'SKEW_CSID',idsk,is_available,lsubmodel)
158 CALL hm_get_floatv(
'P_Thick_Fail',pthk,is_available,lsubmodel,unitab)
159 CALL hm_get_floatv(
'Hm',geo(13),is_available,lsubmodel,unitab)
160 CALL hm_get_floatv(
'Hf',geo(14),is_available,lsubmodel,unitab)
161 CALL hm_get_floatv(
'Hr',geo(15),is_available,lsubmodel,unitab)
162 CALL hm_get_floatv(
'Dm',geo(16),is_available,lsubmodel,unitab)
163 CALL hm_get_floatv(
'Dn',geo(17),is_available,lsubmodel,unitab)
164 CALL hm_get_floatv(
'THICK',geo(1),is_available,lsubmodel,unitab)
165 CALL hm_get_floatv(
'AREA_SHEAR',geo(38),is_available,lsubmodel,unitab)
169 CALL hm_get_floatv(
'MAT_BETA',phi,is_available,lsubmodel,unitab)
174 .
CALL subrotvect(vx,vy,vz,rtrans,sub_id,lsubmodel)
179 IF (pthk == zero) pthk = one-em06
180 pthk =
min(pthk, one)
181 pthk =
max(pthk,-one)
184 IF(ihbe==0)ihbe=ihbe_d
186 IF (ihbe == 4 .AND. ish3n==0 .AND. ish3n_d == 1)
THEN
188 . msgtype=msgwarning,
189 . anmode=aninfo_blind_1,
193 IF (ihbe==22.OR.ihbe==23)
THEN
195 . msgtype=msgwarning,
196 . anmode=aninfo_blind_1,
201 IF(ish3n==0) ish3n = ish3n_d
203 IF (geo(16) == zero) igeo(31) = 1
206 IF (cvis==zero) cvis=one
207 IF (geo(17)==zero) geo(17)=zep015
208 IF (geo(16)==zero)
THEN
214 IF(ismstr==0)ismstr=isst_d
215 IF (isst_d == -2) ismstr = -1
217 IF(geo(13)==zero)geo(13)=em01
218 IF(geo(14)==zero)geo(14)=em01
219 IF(geo(15)==zero)geo(15)=em02
221 IF(geo(13)==zero)geo(13)=em02
222 IF(geo(14)==zero)geo(14)
223 IF(geo(15)==zero)geo(15)=em02
225 IF(isrot==0)isrot=idril_d
226 IF(isrot==2) isrot = 0
229 IF (ismstr==10.AND.isrot>0.AND.idrot==0) idrot = 1
232 IF(n2d>0.AND.ihbe/=0.AND.ihbe/=2)
THEN
236 . msgtype=msgwarning,
237 . anmode=aninfo_blind_2,
261 IF (ihbe>11.AND.ihbe<29)
THEN
269 IF(ismstr==0)ismstr=2
278 IF (geo(38) == zero)geo(38)=five_over_6
279 IF (nip == -1)nip=npts_d
280 IF (nip == 0) nip = 1
281 IF (nip == 1) geo(38)= zero
282 an=sqrt(vx*vx+vy*vy+vz*vz)
314 IF(idsk == iskn(4,j+1))
THEN
319 CALL fretitl2(titr,igeo(npropgi-ltitr+1),ltitr)
330 IF ((irp==22.OR.irp==25).AND.isk==0)
THEN
342 IF(ithk == 0)ithk=ithk_d
343 IF(ithk_d==-2)ithk=-1
344 IF(ishear == 0)ishear=ishea_d
345 IF(iplast == 0)iplast=ipla_d
346 IF(ipla_d==-2) iplast=-1
356 CALL fretitl2(titr,igeo(npropgi-ltitr+1),ltitr)
361 & 5x,
'ORTHOTROPIC SHELL PROPERTY SET'/,
362 & 5x,
'------------------------------'/,
363 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
364 & 5x,
'CONFIDENTIAL DATA'//)
366 IF (ihbe>11.AND.ihbe<29)
THEN
367 WRITE(iout,1011)ig,nip,istrain,geo(1),iss,ihbe,
369 . geo(16),geo(13),geo(38),geo(42),ishear,ithk,iplast,
370 . geo(7),geo(8),geo(9),geo(10),idsk,igeo(14),ipos
372 WRITE(iout,1010)ig,nip,istrain,geo(1),iss,ihbe,
374 . geo(13),geo(14),geo(15),geo(16),
375 . geo(38),geo(42),ishear,ithk,iplast,
376 . geo(7),geo(8),geo(9),geo(10),idsk,igeo(14),ipos
381 CALL fretitl2(titr,igeo(npropgi-ltitr+1),ltitr)
384 . anmode=aninfo_blind_1,
402 ELSEIF(ihbe>=3.AND.ihbe<100.AND.ihbe/=4)
THEN
410 ELSEIF(ishear==1)
THEN
412 ELSEIF(ishear==2)
THEN
420 ELSEIF (ipos==4)
THEN
428 prop_tag(igtyp)%G_SIG = 0
429 prop_tag(igtyp)%G_FOR = 5
430 prop_tag(igtyp)%G_MOM = 3
431 prop_tag(igtyp)%G_THK = 1
432 prop_tag(igtyp)%G_EINT= 2
433 prop_tag(igtyp)%G_EINS= 0
434 prop_tag(igtyp)%G_AREA= 1
435 prop_tag(igtyp)%L_SIG = 5
437 prop_tag(igtyp)%L_THK = 0
438 prop_tag(igtyp)%L_EINT= 2
439 prop_tag(igtyp)%L_EINS= 0
440 prop_tag(igtyp)%G_VOL = 0
441 prop_tag(igtyp)%L_VOL = 0
442 prop_tag(igtyp)%LY_DMG = 2
444 prop_tag(igtyp)%LY_GAMA = 6
445 prop_tag(igtyp)%LY_DIRA = 2
447 prop_tag(igtyp)%LY_PLAPT = 1
448 prop_tag(igtyp)%LY_SIGPT = 5
449 prop_tag(igtyp)%G_FORPG = 5
450 prop_tag(igtyp)%G_MOMPG = 3
451 prop_tag(igtyp)%G_STRPG = 8
460 IF(geo( 3)/=zero.AND.igeo( 5)== 0)igeo( 5)=nint(geo( 3))
462 IF(geo(39)/=zero.AND.igeo( 9)== 0)igeo( 9)=nint(geo(39
463 IF(geo(171)/=zero.AND.igeo(10)== 0)
464 . igeo(10)=nint(geo(171))
470 & 5x,
'ORTHOTROPIC SHELL PROPERTY SET'/,
471 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
472 & 5x,
'NUMBER OF INTEGRATION POINTS. . . . . .=',i10/,
473 & 5x,
'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
474 & 5x,
'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
475 & 5x,
'SMALL STRAIN . . . . . . . . . . . . .=',i10/,
476 & 5x,
'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
477 & 5x,
'3NODE SHELL FORMULATION FLAG. . . . . .=',i10
478 & 5x,
'SHELL HOURGLASS MEMBRANE DAMPING. . . .=',1pg20.13/,
479 & 5x,
'SHELL HOURGLASS FLEXURAL DAMPING. . . .=',1pg20.13/,
480 & 5x,
'SHELL HOURGLASS ROTATIONAL DAMPING. . .=',1pg20.13/,
481 & 5x,
'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
482 & 5x,
'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
483 & 5x,
'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
484 & 5x,
' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
485 & 5x,
' < 0.0 : FRACTION OF FAILED INTG. POINTS',/,
486 & 5x,
'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
487 & 5x,
'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
488 & 5x,
'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
489 & 5x,
'X COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
490 & 5x,
'Y COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
491 & 5x,
'Z COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
492 & 5x,
'ANGLE (DIR 1,PROJ(DIR 1 / SHELL). . . .=',1pg20.13/,
493 & 5x,
'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10/,
494 & 5x,
'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/,
495 & 5x,
'SHELL OFFSET POSITION FLAG . . . . . . =',i10/)
497 & 5x,
'ORTHOTROPIC SHELL PROPERTY SET'/,
498 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
499 & 5x,
'NUMBER OF INTEGRATION POINTS. . . . . .=',i10/,
500 & 5x,
'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
501 & 5x,
'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
502 & 5x,
'SMALL STRAIN . . . . . . . . . . . . .=',i10/,
503 & 5x,
'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
504 & 5x,
'3NODE SHELL FORMULATION FLAG. . . . . .=',i10/,
505 & 5x,
'DRILLING D.O.F. FLAG . . . . . . . . .=',i10/,
506 & 5x,
'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
507 & 5x,
'SHELL NUMERICAL DAMPING . . . . . . . .=',1pg20.13/,
508 & 5x,
'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
509 & 5x,
'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
510 & 5x,
' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
511 & 5x,
' < 0.0 : FRACTION OF FAILED INTG. POINTS',/,
512 & 5x,
'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
513 & 5x,
'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
514 & 5x,
'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
515 & 5x,
'X COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
516 & 5x,
'Y COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
517 & 5x,
'Z COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
518 & 5x,
'ANGLE (DIR 1,PROJ(DIR 1 / SHELL). . . .=',1pg20
519 & 5x,
'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10/,
520 & 5x,
'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/,
521 & 5x,
'SHELL OFFSET POSITION FLAG . . . . . . =',i10/)
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)