OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_prop17.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_prop17 ../starter/source/properties/shell/hm_read_prop17.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!|| fretitl2 ../starter/source/starter/freform.F
30!|| hm_get_float_array_2indexes ../starter/source/devtools/hm_reader/hm_get_float_array_2indexes.F
31!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
32!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
33!|| hm_get_int_array_2indexes ../starter/source/devtools/hm_reader/hm_get_int_array_2indexes.F
34!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
35!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
36!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
37!|| subrotvect ../starter/source/model/submodel/subrot.F
38!||--- uses -----------------------------------------------------
39!|| defaults_mod ../starter/source/modules/defaults_mod.F90
40!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
41!|| message_mod ../starter/share/message_module/message_mod.F
42!|| stack_mod ../starter/share/modules1/stack_mod.F
43!|| submodel_mod ../starter/share/modules1/submodel_mod.F
44!||====================================================================
45 SUBROUTINE hm_read_prop17(GEO ,IGEO ,PM ,IPM ,ISKN ,
46 . UNITAB ,RTRANS ,LSUBMODEL,SUB_ID ,IDTITL ,
47 . PROP_ID ,PROP_TAG , STACK_INFO,DEFAULTS_SHELL)
48C-----------------------------------------------
49C M o d u l e s
50C-----------------------------------------------
51 USE unitab_mod
52 USE elbuftag_mod
53 USE submodel_mod
54 USE message_mod
55 USE stack_mod
56 USE names_and_titles_mod , ONLY : nchartitle
57 USE defaults_mod
58C-----------------------------------------------
59C I m p l i c i t T y p e s
60C-----------------------------------------------
61#include "implicit_f.inc"
62C-----------------------------------------------
63C C o m m o n B l o c k s
64C-----------------------------------------------
65#include "units_c.inc"
66#include "com01_c.inc"
67#include "com04_c.inc"
68#include "param_c.inc"
69#include "scr16_c.inc"
70#include "scr17_c.inc"
71#include "sphcom.inc"
72#include "tablen_c.inc"
73C-----------------------------------------------
74C D u m m y A r g u m e n t s
75C-----------------------------------------------
76 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
77 INTEGER ::
78 . igeo(npropgi),ipm(npropmi,*),sub_id,iskn(liskn,*),
79 . prop_id
80 my_real
81 . geo(npropg), pm(npropm,*),rtrans(ntransf,*)
82 TYPE(submodel_data) LSUBMODEL(*)
83 CHARACTER(LEN=NCHARTITLE)::IDTITL
84 TYPE(prop_tag_) , DIMENSION(0:MAXPROP) :: PROP_TAG
85 TYPE(stack_info_) , TARGET :: STACK_INFO
86 TYPE(shell_defaults_), INTENT(IN) :: DEFAULTS_SHELL
87C-----------------------------------------------
88C L o c a l V a r i a b l e s
89C-----------------------------------------------
90 CHARACTER LAW_ID*4
91 INTEGER :: IGMAT, IHBE , ISMSTR, ISROT, ISHXFEM ,ISTRAIN, ITHK, IPLAST, IDSK,
92 . iorth, ipos , lamin , nsub , nisub , nn, nply, idsub , i, ply_id,
93 . imid_pi,inter ,ipid1 , ipid2 , irep, ishear ,ihbeoutp ,ipid0,
94 . isk , j, iss, n1, imid, k, nc, kk, ii, npt_sub, m1, imat,
95 . n,ipid,ishell,ish3n,id,id1,iun,nply_sub,igtyp,fail_shell,is,irp
96 my_real ::
97 . an, vx, vy, vz, thickt, dt, tmin, tmax,pthk,cvis,pun,zshift, thk ,
98 . ashear, ang, pos, hm, hf, hr, dn, dm,visc_int
99 INTEGER IHBE_D,IPLA_D,ISTR_D,ITHK_D,ISHEA_D,ISST_D,
100 . ISH3N_D, ISTRA_D,NPTS_D,IDRIL_D
101 INTEGER, DIMENSION(:),ALLOCATABLE :: IDMAT_INTP,IDMAT_SUB
102
103 CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
104 LOGICAL :: IS_AVAILABLE, IS_ENCRYPTED
105C-----------------------------------------------
106C D a t a
107C-----------------------------------------------
108 DATA iun/1/
109 DATA pun/0.1/
110C======================================================================|
111C----------------------------------------------------------------
112C COMPOSITE LAYERED SHELL
113C LAYERS WITH : -VARIABLE THICKNESS
114C -VARIABLE MATERIAL (BUT LAW 25 OR 27 ONLY)
115C----------------------------------------------------------------
116C======================================================================|
117 is_available = .false.
118 is_encrypted = .false.
119c
120 igtyp = 17
121 igmat = 1
122 cvis = one
123 istrain = 1
124 irp = 0
125 idsk = 0
126!
127 ihbe_d = defaults_shell%ishell
128 ish3n_d= defaults_shell%ish3n
129 isst_d = defaults_shell%ismstr
130 ipla_d = defaults_shell%iplas
131 ithk_d = defaults_shell%ithick
132 idril_d= defaults_shell%idrill
133 ishea_d = 0
134 npts_d = 0
135 istra_d = 1
136C--------------------------------------------------
137C EXTRACT DATA (IS OPTION CRYPTED)
138C--------------------------------------------------
139 CALL hm_option_is_encrypted(is_encrypted)
140C======================================================================|
141 ! Card 1
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)
146 CALL hm_get_intv('plyxfem',ishxfem,is_available,lsubmodel)
147 CALL hm_get_floatv('Z0',zshift,is_available,lsubmodel,unitab)
148 CALL hm_get_floatv('Vinterply',visc_int,is_available,lsubmodel,unitab) ! global interply viscosity for plyxfem formulation
149 ! Card 2
150
151 CALL hm_get_floatv('Hm',hm,is_available,lsubmodel,unitab)
152 CALL hm_get_floatv('Hf',hf,is_available,lsubmodel,unitab)
153 CALL hm_get_floatv('Hr',hr,is_available,lsubmodel,unitab)
154 CALL hm_get_floatv('Dm',dm,is_available,lsubmodel,unitab)
155 CALL hm_get_floatv('Dn',dn,is_available,lsubmodel,unitab)
156 ! card3
157!! CALL HM_GET_INTV ('ISTRAIN' ,ISTRAIN ,IS_AVAILABLE, LSUBMODEL) ! removed
158 CALL hm_get_floatv('P_Thick_Fail', pthk ,is_available, lsubmodel, unitab)
159 CALL hm_get_floatv('THICK' ,thk ,is_available, lsubmodel, unitab)
160 CALL hm_get_floatv('AREA_SHEAR',ashear ,is_available, lsubmodel, unitab)
161 CALL hm_get_intv ('ITHICK' ,ithk ,is_available, lsubmodel)
162 CALL hm_get_intv ('IPLAS' ,iplast ,is_available, lsubmodel)
163!! CALL HM_GET_FLOATV('Fexp' ,FAIL_SHELL ,IS_AVAILABLE, LSUBMODEL, UNITAB) not existing
164 ! card4
165 CALL hm_get_floatv('V_X' ,vx ,is_available, lsubmodel, unitab)
166 CALL hm_get_floatv('V_Y' ,vy ,is_available, lsubmodel, unitab)
167 CALL hm_get_floatv('V_Z' ,vz ,is_available, lsubmodel, unitab)
168 CALL hm_get_intv('SKEW_CSID' ,idsk ,is_available, lsubmodel)
169 CALL hm_get_intv('Iorth' ,iorth ,is_available, lsubmodel)
170 CALL hm_get_intv('Ipos' ,ipos ,is_available, lsubmodel)
171 CALL hm_get_intv('Ip',irp,is_available,lsubmodel)
172c--------------------------------------------
173c Read ply input cards from laminate.cfg & sub_laminate.cfg
174c either using list of plies or list of sub_stacks with interfaces
175c Fill up STACK_INFO data base
176c--------------------------------------------
177 CALL hm_get_intv('laminateconfig' ,lamin, is_available, lsubmodel)
178c
179 nsub = 0
180 nisub = 0
181 IF (lamin > 0) THEN
182 CALL hm_get_intv('sublaminateidlistmax' ,nsub, is_available, lsubmodel)
183 CALL hm_get_intv('interfacepairsize' ,nisub, is_available, lsubmodel)
184 nply = 0
185 DO is = 1,nsub
186 CALL hm_get_int_array_index('plyidlistmax',nply_sub,is,is_available,lsubmodel)
187 CALL hm_get_int_array_index('dummy',IDSUB,IS,IS_AVAILABLE,LSUBMODEL) !
188 !!
189 STACK_INFO%SUB(2*(IS - 1) + 1) = IDSUB
190 STACK_INFO%SUB(2*(IS - 1) + 2) = NPLY_SUB
191 DO I = 1,NPLY_SUB
192 CALL HM_GET_INT_ARRAY_2INDEXES ('plyidlist',PLY_ID,IS,I,IS_AVAILABLE,LSUBMODEL)
193 CALL HM_GET_FLOAT_ARRAY_2INDEXES('prop_phi' ,ANG ,IS,I,IS_AVAILABLE,LSUBMODEL,UNITAB)
194 CALL HM_GET_FLOAT_ARRAY_2INDEXES('prop_zi' ,POS ,IS,I,IS_AVAILABLE,LSUBMODEL,UNITAB)
195 CALL HM_GET_INT_ARRAY_2INDEXES ('prop_mi' ,IMID_PI,IS,I,IS_AVAILABLE,LSUBMODEL)
196 ! stockage
197 STACK_INFO%PID(NPLY + I) = PLY_ID
198 STACK_INFO%ANG(NPLY + I) = ANG
199 STACK_INFO%POS(NPLY + I) = POS
200 STACK_INFO%MID_IP(NPLY + I) = IMID_PI
201 END DO
202 NPLY = NPLY + NPLY_SUB
203 END DO
204c
205 IF (NISUB > 0) THEN
206 DO INTER = 1,NISUB
207 CALL HM_GET_INT_ARRAY_2INDEXES ('interfacepairplyids' ,IPID1 ,1 ,INTER,IS_AVAILABLE,LSUBMODEL)
208 CALL HM_GET_INT_ARRAY_2INDEXES ('interfacepairplyids' ,IPID2 ,2 ,INTER,IS_AVAILABLE,LSUBMODEL)
209 CALL HM_GET_INT_ARRAY_INDEX ('interfacepair_prop_mi' ,IMID_PI,INTER,IS_AVAILABLE,LSUBMODEL)
210 !!
211 STACK_INFO%ISUB(3*(INTER - 1) + 1) = IPID1
212 STACK_INFO%ISUB(3*(INTER - 1) + 2) = IPID2
213 STACK_INFO%ISUB(3*(INTER - 1) + 3) = IMID_PI
214 END DO
215 END IF
216 ELSE ! property defined by a list of plies
217 CALL HM_GET_INTV('plyidlistmax' ,NPLY ,IS_AVAILABLE ,LSUBMODEL)
218 DO I=1,NPLY
219 CALL HM_GET_INT_ARRAY_INDEX ('plyidlist' ,PLY_ID,I,IS_AVAILABLE,LSUBMODEL)
220 CALL HM_GET_FLOAT_ARRAY_INDEX('prop_phi',ANG,I,IS_AVAILABLE,LSUBMODEL,UNITAB)
221 CALL HM_GET_FLOAT_ARRAY_INDEX('prop_zi' ,POS,I,IS_AVAILABLE,LSUBMODEL,UNITAB)
222 CALL HM_GET_INT_ARRAY_INDEX ('prop_mi' ,IMID_PI,I,IS_AVAILABLE,LSUBMODEL) !
223 !!
224 STACK_INFO%PID(I) = PLY_ID
225 STACK_INFO%ANG(I) = ANG
226 STACK_INFO%POS(I) = POS
227 STACK_INFO%MID_IP(I) = IMID_PI
228 END DO
229 END IF
230c
231 IREP = IORTH
232C IPOS = 0
233C IDSK = 0
234C----------------------
235C Default value
236C----------------------
237 IF (PTHK == ZERO) PTHK = ONE-EM06
238 PTHK = MIN(PTHK, ONE)
239 PTHK = MAX(PTHK,-ONE)
240C ISHEAR NEVER USED
241 ISHEAR = 0
242C CVIS
243 CVIS = ZERO
244
245C Double stockage temporaire - supprimer GEO(12,I)=IGTYP apres tests
246 IGEO( 1)= PROP_ID
247C
248C recommended formulation Iplyxfem=2 (old formulation iplyxfem=1)
249 IF(ISHXFEM > 0) IPLYXFEM = 2
250.AND. IF(ISHXFEM > 0 ISHXFEM /= 2 ) THEN
251 CALL ANCMSG(MSGID=1607,
252 . MSGTYPE=MSGWARNING,
253 . ANMODE=ANINFO_BLIND_2,
254 . I1=PROP_ID,
255 . C1=IDTITL)
256 ENDIF
257 ISHXFEM = MIN(1,ISHXFEM)
258C
259.AND. IF(IHBE /= 12 ISHXFEM >0) THEN
260 ISHXFEM = 0
261 IPLYXFEM = 0
262 CALL ANCMSG(MSGID=726,
263 . MSGTYPE=MSGWARNING,
264 . ANMODE=ANINFO_BLIND_2,
265 . I1=PROP_ID,
266 . C1=IDTITL)
267 ENDIF
268 IF(IHBE==0)IHBE=IHBE_D
269 IHBEOUTP=IHBE
270.AND..AND. IF (IHBE == 4 ISH3N==0 ISH3N_D == 1) THEN
271 CALL ANCMSG(MSGID=680,
272 . MSGTYPE=MSGWARNING,
273 . ANMODE=ANINFO_BLIND_1,
274 . I1=PROP_ID,
275 . C1=IDTITL)
276 ENDIF
277.OR. IF (IHBE==22IHBE==23) THEN
278 CALL ANCMSG(MSGID=539,
279 . MSGTYPE=MSGWARNING,
280 . ANMODE=ANINFO_BLIND_1,
281 . I1=PROP_ID,
282 . C1=IDTITL)
283 IHBE=24
284 ENDIF
285 IF(ISH3N==0) ISH3N = ISH3N_D
286 IGEO(18) = ISH3N
287 IF (DM == ZERO) IGEO(31) = 1
288
289 IF (IHBE==24) THEN
290 IF (CVIS==ZERO) CVIS = ONE
291 IF (DN == ZERO) DN = ZEP015
292 IF (DM == ZERO) THEN
293C-------------remove in cgrtails
294C IF (IGTYP==1.OR.IGTYP==9) GEO(16,I)=ZEP015
295 END IF
296 ENDIF
297C---
298 IF(ISMSTR==0)ISMSTR=ISST_D
299 IF(IHBE==3)THEN
300 IF(HM == ZERO)HM = EM01
301 IF(HF == ZERO)HF = EM01
302 IF(HR == ZERO)HR = EM02
303 ELSE
304 IF(HM == ZERO)HM = EM02
305 IF(HF == ZERO)HF = EM02
306 IF(HR == ZERO)HR = EM02
307 ENDIF
308
309 IF(ASHEAR == ZERO) ASHEAR=FIVE_OVER_6
310C
311.AND. IF (IHBE>11IHBE<29) THEN
312C---------- geo (13, i) is used for DN Porte; Geo (17, I), Cvis changes -----
313 HM = DN
314 DN = CVIS
315C IGEO(20,I)=ISROT
316 ENDIF
317C
318C Double stockage temporaire - supprimer GEO(12,I)=IGTYP apres tests
319 IGEO(10)=IHBE
320 GEO(171)=IHBE
321C
322 IF(ISROT==0)ISROT=IDRIL_D
323 IF(ISROT==2) ISROT = 0
324 IGEO(20)=ISROT
325C-------to have DR----
326.AND..AND. IF (ISMSTR == 10 ISROT > 0 IDROT == 0 ) IDROT = 1
327 IF(ISMSTR == 0)ISMSTR=2
328.AND. IF(ISMSTR == 3. AND.IHBE /= 0 IHBE /= 2) THEN
329 ISMSTR = 2
330 CALL ANCMSG(MSGID=319,
331 . MSGTYPE=MSGWARNING,
332 . ANMODE=ANINFO_BLIND_2,
333 . I1=PROP_ID,
334 . C1=IDTITL)
335 ENDIF
336 GEO(3) =ISMSTR
337 IGEO(5)=ISMSTR
338C
339 IF (SUB_ID /= 0)
340 . CALL SUBROTVECT(VX,VY,VZ,RTRANS,SUB_ID,LSUBMODEL)
341C---
342 GEO(1) = THK
343 IGEO(32) = 0 ! ISLV not used (may be needed for plyxfem)
344 IF(GEO(38) == ZERO) GEO(38) = FIVE_OVER_6
345 !! IF(NLY == 1) GEO(38)= ZERO
346 FAIL_SHELL = 0
347 AN=SQRT(VX*VX+VY*VY+VZ*VZ)
348 IF(AN < EM10)THEN
349 VX=ONE
350 VY=ZERO
351 VZ=ZERO
352 IF (IRP==23) THEN
353 CALL ANCMSG(MSGID=1922,
354 . MSGTYPE=MSGERROR,
355 . ANMODE=ANINFO,
356 . C1='property',
357 . I1=PROP_ID,
358 . C2='property',
359 . C3=TITR,
360 . I2=IRP)
361 END IF
362 ELSE
363 VX=VX/AN
364 VY=VY/AN
365 VZ=VZ/AN
366 ENDIF
367 ISK = 0
368 IF (IDSK/=0) THEN
369 DO J=0,NUMSKW+MIN(IUN,NSPCOND)*NUMSPH+NSUBMOD
370 IF(IDSK == ISKN(4,J+1)) THEN
371 ISK=J+1
372 GO TO 10
373 ENDIF
374 END DO
375 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
376 CALL ANCMSG(MSGID=184,
377 . MSGTYPE=MSGERROR,
378 . ANMODE=ANINFO,
379 . C1='property',
380 . I1=PROP_ID,
381 . C2='property',
382 . C3=TITR,
383 . I2=IDSK)
384 10 CONTINUE
385 ENDIF
386.OR..AND. IF ((IRP==22IRP==25)ISK==0) THEN
387 CALL ANCMSG(MSGID=1923,
388 . MSGTYPE=MSGERROR,
389 . ANMODE=ANINFO,
390 . C1='property',
391 . I1=PROP_ID,
392 . C2='property',
393 . C3=TITR,
394 . I2=IRP)
395 END IF
396c check duplicated py IDs
397 IPID0 = STACK_INFO%PID(1)
398 DO K=2,NPLY
399 IF (STACK_INFO%PID(K) == IPID0) THEN
400 CALL ANCMSG(MSGID=1584,MSGTYPE=MSGERROR,ANMODE=ANINFO_BLIND_2,
401 . I1=PROP_ID,
402 . I2=IPID0)
403 ENDIF
404 ENDDO
405c
406 IF(ITHK == 0)ITHK=ITHK_D
407 IF(ISHEAR == 0)ISHEAR=ISHEA_D
408 IF(IPLAST == 0)IPLAST=IPLA_D
409!
410 IHBE = IGEO(10)
411 ISS = IGEO(5)
412C======================================================================|
413 IF (DM == ZERO) IGEO(31) = 1
414 IGEO( 1) = PROP_ID
415 IGEO(2) = ISK
416 IGEO(5) = ISMSTR
417 IGEO(6) = IORTH ! IREP
418 IGEO(11) = IGTYP
419 IGEO(18) = ISH3N
420 IGEO(19) = ISHXFEM
421 IGEO(20) = ISROT
422 IGEO(14) = IRP
423!! IGEO(47) = IINT
424 !! is not necessary it done in translation : before 2017 igmat=-1
425 !! the new global mat is by default since 2017 version IGMAT IS USED GENERALLY > 0
426 IGEO(98) = IGMAT
427c
428 GEO(3) = ISMSTR
429 GEO(7) = VX
430 GEO(8) = VY
431 GEO(9) = VZ
432 GEO(11) = ISTRAIN
433 GEO(12) = IGTYP
434 GEO(13) = HM
435 GEO(14) = HF
436 GEO(15) = HR
437 GEO(16) = DM
438 GEO(17) = DN
439 GEO(20) = VISC_INT ! for ply-xfem formulation
440 GEO(35) = ITHK
441 GEO(37) = ISHEAR
442 GEO(38) = ASHEAR
443 GEO(39) = IPLAST
444 GEO(42) = PTHK
445 GEO(43) = ONE
446 GEO(171)= IHBE ! double stockage
447 GEO(199)= ZSHIFT
448C
449 IHBE=NINT(GEO(171))
450 IF(IHBE==0)THEN
451 GEO(171)=0
452 ELSEIF(IHBE==1)THEN
453 GEO(171)= 1
454 ELSEIF(IHBE==2)THEN
455 GEO(171)=0
456.AND..AND. ELSEIF(IHBE>=3IHBE<100IHBE/=4)THEN
457 GEO(171)=IHBE-1
458 ENDIF
459
460 ISTRAIN=NINT(GEO(11))
461 IF(ISTRAIN==0)THEN
462 GEO(11)=0
463 ELSEIF(ISTRAIN==1)THEN
464 GEO(11)=1
465 ELSEIF(ISTRAIN==2)THEN
466 GEO(11)=0
467 ENDIF
468
469 ISHEAR = GEO(37)
470 IF(ISHEAR==0)THEN
471 GEO(37)=0
472 ELSEIF(ISHEAR==1)THEN
473 GEO(37)=1
474 ELSEIF(ISHEAR==2)THEN
475 GEO(37)=0
476 ENDIF
477!! N1 it not used only with old input format
478!! IGEO(IPPID + N1) = IPOS
479 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
480C-----------------------------------------------------------------
481 IF(IS_ENCRYPTED)THEN
482 WRITE(IOUT,1000)PROP_ID
483 1000 FORMAT(
484 & 5X,'composite layered shell property set'/,
485 & 5X,'------------------------------------'/,
486 & 5X,'property set number . . . . . . . . . .=',I10/,
487 & 5X,'confidential data'//)
488 ELSE
489 IF(ISK==0)THEN
490.AND. IF (IHBE > 11IHBE < 29) THEN
491 WRITE(IOUT,2112)PROP_ID,ISTRAIN,GEO(1),ISS,IHBE,
492 . ISH3N,IGEO(20),
493 . GEO(16),GEO(13),GEO(38),PTHK,FAIL_SHELL,ISHEAR,ITHK,
494 . IPLAST,IREP,GEO(7),GEO(8),GEO(9),IGEO(14)
495 ELSE
496 WRITE(IOUT,2110)PROP_ID,ISTRAIN,GEO(1),ISS,IHBE, ! nly of less
497 . ISH3N,
498 . GEO(13),GEO(14),GEO(15),GEO(16),
499 . GEO(38),PTHK,FAIL_SHELL,ISHEAR,ITHK,IPLAST,IREP,
500 . GEO(7),GEO(8),GEO(9),IGEO(14)
501 ENDIF
502 ELSE
503.AND. IF (IHBE > 11IHBE < 29) THEN
504 WRITE(IOUT,2113)PROP_ID,ISTRAIN,GEO(1),ISS,IHBE,
505 . ISH3N,IGEO(20),
506 . GEO(16),GEO(13),GEO(38),PTHK,FAIL_SHELL,ISHEAR,ITHK,
507 . IPLAST,IREP,IDSK,IGEO(14)
508 ELSE
509 WRITE(IOUT,2111)PROP_ID,ISTRAIN,GEO(1),ISS,IHBE,
510 . ISH3N,
511 . GEO(13),GEO(14),GEO(15),GEO(16),GEO(38),PTHK,
512 . FAIL_SHELL,ISHEAR,ITHK,IPLAST,IREP,IDSK,IGEO(14)
513 ENDIF
514 ENDIF
515 IF(ISHXFEM > 0) WRITE(IOUT, 2114)
516 ENDIF
517
518C---
519!!! IF (NLY>NLYMAX) THEN
520!! CALL ANCMSG(MSGID=28,
521!! . MSGTYPE=MSGERROR,
522!! . ANMODE=ANINFO_BLIND_1,
523!! . I1=IG,
524!! . C1=TITR,
525!! . I2=NLYMAX)
526!! ENDIF
527C
528 NC = NPLY
529 ALLOCATE( IDMAT_INTP(NPLY),IDMAT_SUB(NISUB))
530 IDMAT_INTP = 0
531 IDMAT_SUB = 0
532 N1 = NPLY
533 GEO(6)=N1 + EM01
534 IGEO(4) = N1
535C
536 IGEO(43) = NSUB ! number of substack
537 IGEO(44) = NISUB ! number of interface
538 DO 250 K = 1, NPLY
539 IMID = STACK_INFO%MID_IP(K)
540 IDMAT_INTP(K) = IMID
541.AND. IF( IMID == 0 ISHXFEM > 0) THEN
542 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
543 CALL ANCMSG(MSGID=2019,
544 . MSGTYPE=MSGERROR,
545 . ANMODE=ANINFO_BLIND_2,
546 . I1=PROP_ID,
547 . C1=TITR,
548 . I2=IMID)
549 ENDIF
550 DO J=1,NUMMAT
551 IF(IPM(1,J) == IMID) THEN
552 STACK_INFO%MID_IP(K) = J
553 GO TO 250
554 ENDIF
555 ENDDO
556 IF(ISHXFEM > 0) THEN
557 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
558 CALL ANCMSG(MSGID=2019,
559 . MSGTYPE=MSGERROR,
560 . ANMODE=ANINFO_BLIND_2,
561 . I1=PROP_ID,
562 . C1=TITR,
563 . I2=IMID)
564 ENDIF
565 STACK_INFO%MID_IP(K) = 0
566 250 CONTINUE
567C isub stack
568 IF(IGEO(44) > 0) THEN
569 DO 300 K = 1, NISUB
570 IMID = STACK_INFO%ISUB (3*(K-1) + 3)
571 IDMAT_SUB(K) = IMID
572.AND. IF( IMID == 0 ISHXFEM > 0) THEN
573 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
574 CALL ANCMSG(MSGID=31,
575 . MSGTYPE=MSGERROR,
576 . ANMODE=ANINFO_BLIND_2,
577 . I1=PROP_ID,
578 . C1=TITR,
579 . I2=IMID)
580 ENDIF
581 DO J=1,NUMMAT
582 IF(IPM(1,J) == IMID) THEN
583 STACK_INFO%ISUB (3*(K-1) + 3) = J
584 GO TO 300
585 ENDIF
586 ENDDO
587 IF(ISHXFEM > 0) THEN
588 CALL FRETITL2(TITR,IGEO(NPROPGI-LTITR+1),LTITR)
589 CALL ANCMSG(MSGID=31,
590 . MSGTYPE=MSGERROR,
591 . ANMODE=ANINFO_BLIND_2,
592 . I1=PROP_ID,
593 . C1=TITR,
594 . I2=IMID)
595 ENDIF
596 STACK_INFO%ISUB (3*(K-1) + 3) = 0
597 300 CONTINUE
598 ENDIF
599 IGEO(99) = IPOS
600 IF(NSUB > 0) THEN
601 KK = 0
602 DO II = 1,NSUB
603 IDSUB = STACK_INFO%SUB ( 2*(II - 1) + 1)
604 NPT_SUB = STACK_INFO%SUB ( 2*(II - 1) + 2)
605C
606.EQV. IF(IS_ENCRYPTED .FALSE. )WRITE(IOUT,3000)IDSUB
607 DO K=1,NPT_SUB
608 M1= KK + K
609.EQV. IF(IS_ENCRYPTED .FALSE. ) THEN
610 WRITE(IOUT,2121)K,STACK_INFO%PID(M1),STACK_INFO%ANG(M1),
611 . STACK_INFO%POS(M1)
612 IF(ISHXFEM > 0) WRITE(IOUT,2122)IDMAT_INTP(M1),GEO(20)
613 ENDIF
614 STACK_INFO%ANG(M1)=STACK_INFO%ANG(M1)*PI/HUNDRED80
615 ENDDO
616 KK = KK + NPT_SUB
617 ENDDO
618 NISUB= IGEO(44)
619 DO K=1,NISUB
620 IPID1 =STACK_INFO%ISUB (3*(K - 1) + 1)
621 IPID2 =STACK_INFO%ISUB (3*(K - 1) + 2)
622 IMAT =STACK_INFO%ISUB (3*(K - 1) + 3)
623.EQV. IF( IS_ENCRYPTED .FALSE. ) THEN
624 WRITE(IOUT,3100)K
625 IF(ISHXFEM > 0) THEN
626 WRITE(IOUT,3200)IPID1,IPID2,IMAT
627 ELSE
628 WRITE(IOUT,3300)IPID1,IPID2
629 ENDIF
630 ENDIF
631 ENDDO
632 ELSE
633 DO N=1,N1
634.EQV. IF(IS_ENCRYPTED .FALSE. )THEN
635 WRITE(IOUT,2121)N,STACK_INFO%PID(N),STACK_INFO%ANG(N),
636 . STACK_INFO%POS(N)
637.AND. IF(N < N1 ISHXFEM > 0) WRITE(IOUT,2122)IDMAT_INTP(N),GEO(20)
638 ENDIF
639 STACK_INFO%ANG(N)=STACK_INFO%ANG(N)*PI/HUNDRED80
640 ENDDO
641 ENDIF
642C
643 DEALLOCATE (IDMAT_INTP,IDMAT_SUB)
644C----
645 RETURN
646 2110 FORMAT(
647 & 5X,'composite stack shell property set ',
648 & 'with variable thicknesses and materials'/,
649 & 5X,'property set number . . . . . . . . . .=',I10/,
650!! & 5X,'number of layers. . . . . . . . . . . .=',I10/,
651 & 5X,'post processing strain flag . . . . . .=',I10/,
652 & 5X,'shell thickness . . . . . . . . . . . .=',1PG20.13/,
653 & 5X,'small strain flag . . . . . . . . . . .=',I10/,
654 & 5X,'shell formulation flag. . . . . . . . .=',I10/,
655 & 5X,'3node shell formulation flag. . . . . .=',I10/,
656 & 5X,'shell hourglass membrane damping. . . .=',1PG20.13/,
657 & 5X,'shell hourglass flexural damping. . . .=',1PG20.13/,
658 & 5X,'shell hourglass rotational damping. . .=',1PG20.13/,
659 & 5X,'shell membrane damping. . . . . . . . .=',1PG20.13/,
660 & 5X,'shear area reduction factor . . . . . .=',1PG20.13/,
661 & 5X,'element deletion PARAMETER. . . . . . .=',1PG20.13/,
662 & 5X,' > 0.0 : fraction of failed thickness ',/,
663 & 5X,' < 0.0 : fraction of failed layers/plys ',/,
664 & 5X,'shell failure flag. . . . . . . . . . .=',I10/,
665 & 5X,'shear formulation flag. . . . . . . . .=',I10/,
666 & 5X,'thickness variation flag. . . . . . . .=',I10/,
667 & 5X,'plasticity formulation flag . . . . . .=',I10/,
668 & 5X,'local ortothropy system flag. . . . . .=',I10/,
669 & 5X,'x component of dir 1 of orthotropy. . .=',1PG20.13/,
670 & 5X,'y component of dir 1 of orthotropy. . .=',1PG20.13/,
671 & 5X,'z component of dir 1 of orthotropy. . .=',1PG20.13/,
672 & 5X,'reference direction flag in shell plane=',I10/)
673 2111 FORMAT(
674 & 5X,'composite stack shell property set ',
675 & 'with variable thicknesses and materials'/,
676 & 5X,'property set number . . . . . . . . . .=',I10/,
677!! & 5X,'number of layers. . . . . . . . . . . .=',I10/,
678 & 5X,'post processing strain flag . . . . . .=',I10/,
679 & 5X,'shell thickness . . . . . . . . . . . .=',1PG20.13/,
680 & 5X,'small strain flag . . . . . . . . . . .=',I10/,
681 & 5X,'shell formulation flag. . . . . . . . .=',I10/,
682 & 5X,'3node shell formulation flag. . . . . .=',I10/,
683 & 5X,'shell hourglass membrane damping. . . .=',1PG20.13/,
684 & 5X,'shell hourglass flexural damping. . . .=',1PG20.13/,
685 & 5X,'shell hourglass rotational damping. . .=',1PG20.13/,
686 & 5X,'shell membrane damping. . . . . . . . .=',1PG20.13/,
687 & 5X,'shear area reduction factor . . . . . .=',1PG20.13/,
688 & 5X,'element deletion PARAMETER. . . . . . .=',1PG20.13/,
689 & 5X,' > 0.0 : fraction of failed thickness ',/,
690 & 5x,' < 0.0 : FRACTION OF FAILED LAYERS/PLYS ',/,
691 & 5x,'SHELL FAILURE FLAG. . . . . . . . . . .=',i10/,
692 & 5x,'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
693 & 5x,'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
694 & 5x,'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
695 & 5x,'LOCAL ORTOTHROPY SYSTEM FLAG. . . . . .=',i10/,
696 & 5x,'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10/,
697 & 5x,'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/)
698 2112 FORMAT(
699 & 5x,'COMPOSITE STACK SHELL PROPERTY SET ',
700 & 'WITH VARIABLE THICKNESSES AND MATERIALS'/,
701 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
702!! & 5X,'NUMBER OF LAYERS. . . . . . . . . . . .=',I10/,
703 & 5x,'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
704 & 5x,'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
705 & 5x,'SMALL STRAIN FLAG . . . . . . . . . . .=',i10/,
706 & 5x,'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
707 & 5x,'3NODE SHELL FORMULATION FLAG. . . . . .=',i10/,
708 & 5x,'DRILLING D.O.F. FLAG . . . . . . . . .=',i10/,
709 & 5x,'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
710 & 5x,'SHELL NUMERICAL DAMPING . . . . . . . .=',1pg20.13/,
711 & 5x,'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
712 & 5x,'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
713 & 5x,' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
714 & 5x,' < 0.0 : FRACTION OF FAILED LAYERS/PLYS ',/,
715 & 5x,'SHELL FAILURE FLAG. . . . . . . . . . .=',i10/,
716 & 5x,'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
717 & 5x,'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
718 & 5x,'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
719 & 5x,'LOCAL ORTOTHROPY SYSTEM FLAG. . . . . .=',i10/,
720 & 5x,'X COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
721 & 5x,'Y COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
722 & 5x,'Z COMPONENT OF DIR 1 OF ORTHOTROPY. . .=',1pg20.13/,
723 & 5x,'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/)
724 2113 FORMAT(
725 & 5x,'COMPOSITE STACK SHELL PROPERTY SET ',
726 & 'WITH VARIABLE THICKNESSES AND MATERIALS'/,
727 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
728!! & 5X,'NUMBER OF LAYERS. . . . . . . . . . . .=',I10/,
729 & 5x,'POST PROCESSING STRAIN FLAG . . . . . .=',i10/,
730 & 5x,'SHELL THICKNESS . . . . . . . . . . . .=',1pg20.13/,
731 & 5x,'SMALL STRAIN FLAG . . . . . . . . . . .=',i10/,
732 & 5x,'SHELL FORMULATION FLAG. . . . . . . . .=',i10/,
733 & 5x,'3NODE SHELL FORMULATION FLAG. . . . . .=',i10/,
734 & 5x,'DRILLING D.O.F. FLAG . . . . . . . . .=',i10/,
735 & 5x,'SHELL MEMBRANE DAMPING. . . . . . . . .=',1pg20.13/,
736 & 5x,'SHELL NUMERICAL DAMPING . . . . . . . .=',1pg20.13/,
737 & 5x,'SHEAR AREA REDUCTION FACTOR . . . . . .=',1pg20.13/,
738 & 5x,'ELEMENT DELETION PARAMETER. . . . . . .=',1pg20.13/,
739 & 5x,' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
740 & 5x,' < 0.0 : FRACTION OF FAILED LAYERS/PLYS ',/,
741 & 5x,'SHELL FAILURE FLAG. . . . . . . . . . .=',i10/,
742 & 5x,'SHEAR FORMULATION FLAG. . . . . . . . .=',i10/,
743 & 5x,'THICKNESS VARIATION FLAG. . . . . . . .=',i10/,
744 & 5x,'PLASTICITY FORMULATION FLAG . . . . . .=',i10/,
745 & 5x,'LOCAL ORTOTHROPY SYSTEM FLAG. . . . . .=',i10/,
746 & 5x,'SKEW OF THE FIRST ORTHOTROPY DIRECTION.=',i10/,
747 & 5x,'REFERENCE DIRECTION FLAG IN SHELL PLANE=',i10/)
748 2114 FORMAT(
749 & 5x,'PLY XFEM SHELL COMPOSITE FORMULATION IS ACTIVATED '/)
750 2121 FORMAT(
751 & 5x,' PLY ',i3/,
752 & 5x,' PLY PID NUMBER . . . . . . . . .=',i10/
753 & 5x,' ANGLE (DIR 1,PROJ(DIR 1 / SHELL).=',1pg20.13/,
754 & 5x,' POSITION. . . . . . . . . . . . .=',1pg20.13/)
755 2122 FORMAT(
756 & 5x,' INTER-PLY MID NUMBER . . . . . .=',i10/
757 & 5x,' INTER-PLY SHELL DAMPING . . . . . .=',1pg20.13/)
758 3000 FORMAT(
759 & 5x,' COMPOSITE SUBSTACK SHELL ID . . . . . . =',i10/ )
760 3100 FORMAT(
761 & 5x,' INTERFACE NUMBER BETWEEN-SUBSTACK . . .:',i10/ )
762 3200 FORMAT(
763 & 5x,' INTER-PLY_1 PID NUMBER . . . . . =',i10/,
764 & 5x,' INTER-PLY_2 PID NUMBER . . . . . .=',i10/,
765 & 5x,' INTER-PLY MID NUMBER . . . . . . =',i10/)
766 3300 FORMAT(
767 & 5x,' INTER-PLY_1 PID NUMBER . . . . . =',i10/,
768 & 5x,' INTER-PLY_2 PID NUMBER . . . . . .=',i10/)
769 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
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_int_array_index
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
Definition hm_get_int_array_index.F:164
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_read_prop17
subroutine hm_read_prop17(geo, igeo, pm, ipm, iskn, unitab, rtrans, lsubmodel, sub_id, idtitl, prop_id, prop_tag, stack_info, defaults_shell)
Definition hm_read_prop17.F:48
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
id
initmumps id
Definition multiplerhs_example.m:4
elbuftag_mod
Definition elbuftag_mod.F:230
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
stack_mod
Definition stack_mod.F:75
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
elbuftag_mod::prop_tag_
Definition elbuftag_mod.F:377
stack_mod::stack_info_
Definition stack_mod.F:91
submodel_mod::submodel_data
Definition submodel_mod.F:728