42 . RTRANS,LSUBMODEL,IDTITL ,IG ,SUB_ID ,ISKN ,
56#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
84 . ISHEAR, ISTRAIN,IHBE,IPLAST,ITHK,IHBEOUTP,,
86 INTEGER ISH3N,ISROT,SUB_ID,NLY,IRP,IDSK,ISK,IUN
87 INTEGER NLYMAX,N1,M1, IPANG,IPTHK,IPPOS,IPWEIGHT,IPOS
88 INTEGER IHBE_D,IPLA_D,ITHK_D,ISHEA_D,ISST_D,
89 . ish3n_d, istra_d,npts_d,idril_d
92 . pthk, an, phi_i,thk10,zshift
93 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
94 CHARACTER(LEN=NCHARTITLE) :: TITR
105 ihbe_d = defaults_shell%ishell
106 ish3n_d= defaults_shell%ish3n
107 isst_d = defaults_shell%ismstr
108 ipla_d = defaults_shell%iplas
109 ithk_d = defaults_shell%ithick
110 idril_d= defaults_shell%idrill
125 is_encrypted = .false.
126 is_available = .false.
142 CALL hm_get_intv(
'Ishell',ihbe,is_available,lsubmodel)
143 CALL hm_get_intv(
'Ismstr',ismstr,is_available,lsubmodel)
144 CALL hm_get_intv(
'ISH3N',ish3n,is_available,lsubmodel)
145 CALL hm_get_intv(
'Idrill',isrot,is_available,lsubmodel)
148 CALL hm_get_intv(
'ITHICK',ithk,is_available,lsubmodel)
149 CALL hm_get_intv(
'IPLAS',iplast,is_available,lsubmodel)
150 CALL hm_get_intv(
'SKEW_CSID',idsk,is_available,lsubmodel)
151 CALL hm_get_intv(
'Ipos',ipos,is_available,lsubmodel)
156 CALL hm_get_floatv(
'P_Thick_Fail',pthk,is_available,lsubmodel,unitab)
157 CALL hm_get_floatv(
'Hm',geo(13),is_available,lsubmodel,unitab)
158 CALL hm_get_floatv(
'Hf',geo(14),is_available,lsubmodel,unitab)
159 CALL hm_get_floatv(
'Hr',geo(15),is_available,lsubmodel,unitab)
160 CALL hm_get_floatv(
'Dm',geo(16),is_available,lsubmodel,unitab)
161 CALL hm_get_floatv(
'Dn',geo(17),is_available,lsubmodel,unitab)
162 CALL hm_get_floatv(
'THICK',geo(1),is_available,lsubmodel,unitab)
163 CALL hm_get_floatv(
'AREA_SHEAR',geo(38),is_available,lsubmodel,unitab)
171 .
CALL subrotvect(vx,vy,vz,rtrans,sub_id,lsubmodel)
180 IF (pthk == zero) pthk = one-em06
181 pthk =
min(pthk, one)
182 pthk =
max(pthk,-one)
187 IF(ihbe==0)ihbe=ihbe_d
189 IF (ihbe == 4 .AND. ish3n==0 .AND. ish3n_d == 1)
THEN
191 . msgtype=msgwarning,
192 . anmode=aninfo_blind_1,
196 IF (ihbe==22.OR.ihbe==23)
THEN
198 . msgtype=msgwarning,
199 . anmode=aninfo_blind_1,
209 ELSEIF(ihbe>=3.AND.ihbe<100.AND.ihbe/=4)
THEN
212 IF(geo(171)/=zero.AND.igeo(10)== 0)igeo(10)=nint(geo(171))
213 IF(ish3n==0) ish3n = ish3n_d
216 IF(ismstr==0)ismstr=isst_d
217 IF (isst_d == -2) ismstr = -1
218 IF(ismstr==0)ismstr=2
219 IF(ismstr==3.AND.ihbe/=0.AND.ihbe/=2)
THEN
222 . msgtype=msgwarning,
223 . anmode=aninfo_blind_2,
230 IF(isrot==0)isrot=idril_d
231 IF(isrot>=2) isrot = 0
234 IF(ithk == 0)ithk=ithk_d
235 IF(ithk_d==-2)ithk=-1
236 IF(ishear == 0)ishear=ishea_d
237 IF(iplast == 0)iplast=ipla_d
238 IF(ipla_d==-2) iplast=-1
243 IF(ishear==2) geo(37)=0
245 IF(geo(39)/=zero.AND.igeo(9)== 0)igeo(9)=iplast
247 IF (ismstr==10.AND.isrot>0.AND.idrot==0) idrot = 1
250 IF(geo(13)==zero)geo(13)=em01
251 IF(geo(14)==zero)geo(14)=em01
252 IF(geo(15)==zero)geo(15)=em02
254 IF(geo(13)==zero)geo(13)=em02
255 IF(geo(14)==zero)geo(14)=em02
256 IF(geo(15)==zero)geo(15)=em02
258 IF(ismstr==3) geo(5)=ep06
260 IF (geo(16) == zero) igeo(31) = 1
262 IF (ihbe==12) geo(13)=geo(17)
264 IF (geo(17)==zero) geo(17)=zep015
274 IF(geo(38) == zero) geo(38)=five_over_6
275 IF(nly == -1)nly=npts_d
277 IF(nly == 1) geo(38)= zero
278 an=sqrt(vx*vx+vy*vy+vz*vz)
281 geo(ipweight+i) = one
311 IF(idsk == iskn(4,j+1))
THEN
316 CALL fretitl2(titr,igeo(npropgi-ltitr+1),ltitr)
327 IF ((irp==22.OR.irp==25).AND.isk==0)
THEN
344 & 5x,
'ORTHOTROPIC LAYERED SHELL PROPERTY SET'/,
345 & 5x,
'--------------------------------------'/,
346 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
347 & 5x,
'CONFIDENTIAL DATA'//)
349 IF (ihbe>11.AND.ihbe<29)
THEN
350 WRITE(iout,2011)ig,nly,istrain,geo(1),iss,ihbe,
352 . geo(16),geo(13),geo(38),geo(42),ishear,ithk,iplast,
353 . geo(7),geo(8),geo(9),idsk,igeo(14),ipos
355 WRITE(iout,2010)ig,nly,istrain,geo(1),iss,ihbe,
357 . geo(13),geo(14),geo(15),geo(16),
358 . geo(38),geo(42),ishear,ithk,iplast,
359 . geo(7),geo(8),geo(9),idsk,igeo(14),ipos
365 . anmode=aninfo_blind_1,
372 . anmode=aninfo_blind_1,
381 ELSEIF (ipos==4)
THEN
396 geo(ippos+k) = -half*(one+thk10)+k*thk10 + zshift
397 IF(.NOT.is_encrypted)
WRITE(iout,2020) k,geo(m1)
398 geo(m1)=geo(m1)*pi/hundred80
406 prop_tag(igtyp)%G_SIG = 0
407 prop_tag(igtyp)%G_FOR = 5
408 prop_tag(igtyp)%G_MOM = 3
409 prop_tag(igtyp)%G_THK = 1
410 prop_tag(igtyp)%G_EINT= 2
411 prop_tag(igtyp)%G_EINS= 0
412 prop_tag(igtyp)%L_SIG = 5
414 prop_tag(igtyp)%L_THK = 0
415 prop_tag(igtyp)%L_EINT= 2
416 prop_tag(igtyp)%L_EINS= 0
417 prop_tag(igtyp)%G_VOL = 0
418 prop_tag(igtyp)%L_VOL = 0
419 prop_tag(igtyp)%LY_DMG = 2
421 prop_tag(igtyp)%LY_GAMA = 6
422 prop_tag(igtyp)%LY_DIRA = 2
424 prop_tag(igtyp)%LY_PLAPT = 1
425 prop_tag(igtyp)%LY_SIGPT = 5
426 prop_tag(igtyp)%G_FORPG = 5
427 prop_tag(igtyp)%G_MOMPG = 3
428 prop_tag(igtyp)%G_STRPG = 8
436 & 5x,
'ORTHOTROPIC LAYERED SHELL PROPERTY SET'/,
437 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
438 & 5x,
'NUMBER OF LAYERS. . . . . . . . . . . .=',i10/,
439 & 5x,
'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
440 & 5x,
'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
441 & 5x,
'SMALL STRAIN. . . . . . . . . . . . . .=',i10/,
442 & 5x,
'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
443 & 5x,
'3NODE SHELL FORMULATION FLAG. . . . . .=',i10/,
444 & 5x,
'SHELL HOURGLASS MEMBRANE DAMPING. . . .=',1pg20.13/,
445 & 5x,
'SHELL HOURGLASS FLEXURAL DAMPING. . . .=',1pg20.13/,
446 & 5x,
'SHELL HOURGLASS ROTATIONAL DAMPING. . .=',1pg20.13/,
447 & 5x,
'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
448 & 5x,
'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
449 & 5x,
'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
450 & 5x,
' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
451 & 5x,
' < 0.0 : FRACTION OF FAILED LAYERS ',/,
452 & 5x,
'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
453 & 5x,
'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
454 & 5x,
'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
455 & 5x,
'X COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
456 & 5x,
'Y COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
457 & 5x,
'Z COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
458 & 5x,
'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10/,
459 & 5x,
'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/,
460 & 5x,
'SHELL OFFSET POSITION FLAG . . . . . . =',i10/)
462 & 5x,
'ORTHOTROPIC LAYERED SHELL PROPERTY SET'/,
463 & 5x,
'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
464 & 5x,
'NUMBER OF LAYERS. . . . . . . . . . . .=',i10/,
465 & 5x,
'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
466 & 5x,
'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
467 & 5x,
'SMALL STRAIN. . . . . . . . . . . . . .=',i10/,
468 & 5x,
'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
469 & 5x,
'3NODE SHELL FORMULATION FLAG. . . . . .=',i10/,
470 & 5x,
'DRILLING D.O.F. FLAG . . . . . . . . .=',i10/,
471 & 5x,
'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
472 & 5x,
'SHELL NUMERICAL DAMPING . . . . . . . .=',1pg20.13/,
473 & 5x,
'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
474 & 5x,
'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
475 & 5x,
' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
476 & 5x,
' < 0.0 : FRACTION OF FAILED LAYERS ',/,
477 & 5x,
'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
478 & 5x,
'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
479 & 5x,
'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
480 & 5x,
'X COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
481 & 5x,
'Y COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
482 & 5x,
'Z COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
483 & 5x,
'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10
484 & 5x,
'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/,
485 & 5x,
'SHELL OFFSET POSITION FLAG . . . . . . =',i10/)
487 & 5x,
'LAYER,ANGLE (DIR 1,PROJ(DIR 1 / SHELL).=',i10,e12.5)
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)