1 : C IDDILI SOURCE CHAT 11/03/16 21:25:20 6902 2 : 3 : SUBROUTINE IDDILI (ITROP,ITYPE, MOCOMP,NBROBL,NBRFAC) 4 : 5 : *----------------------------------------------------------------------* 6 : * Noms de composantes de materiau pour la DIFFUSION lineaire * 7 : *----------------------------------------------------------------------* 8 : * ITROP (E) "tropie" du materiau (IMATEE du MODELE) * 9 : * ITYPE (E) =0 si MOCOMP est un LISTMOTS, =1 si NOMID * 10 : * MOCOMP (S) pointeur selon ITYPE contenant les composantes requises * 11 : * (obligatoires et facultatives) * 12 : * NBROBL (S) nombre de composantes obligatoires * 13 : * NBRFAC (S) nombre de composantes facultatives (=0 ici) * 14 : *----------------------------------------------------------------------* 15 : IMPLICIT INTEGER(I-N) 16 : IMPLICIT REAL*8(A-H,O-Z) 17 : 18 : -INC CCOPTIO 19 : -INC SMMODEL 20 : -INC SMLMOTS 21 : 22 : NBROBL = 0 23 : NBRFAC = 0 24 : MLMOTS = 0 25 : JGN = 8 26 : JGM = 0 27 : 28 : C- Diffusion isotrope 29 : IF (ITROP .EQ. 1) THEN 30 : JGM = 1 31 : SEGINI,MLMOTS 32 : C*8 MOTS(1) = 'KDIFF ' 33 : MOTS(1) = 'KD ' 34 : 35 : C- Diffusion orthotrope 36 : ELSE IF (ITROP .EQ. 2) THEN 37 : C--- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE 38 : IF (IDIM.EQ.2) THEN 39 : IF (IFOMOD.NE.1) THEN 40 : JGM = 4 41 : SEGINI,MLMOTS 42 : C*8 MOTS(1) = 'KDIFF1 ' 43 : C*8 MOTS(2) = 'KDIFF2 ' 44 : MOTS(1) = 'KD1 ' 45 : MOTS(2) = 'KD2 ' 46 : MOTS(3) = 'V1X ' 47 : MOTS(4) = 'V1Y ' 48 : C--- Elements massifs bidimensionnels FOURIER 49 : ELSE 50 : JGM = 5 51 : SEGINI,MLMOTS 52 : C*8 MOTS(1) = 'KDIFF1 ' 53 : C*8 MOTS(2) = 'KDIFF2 ' 54 : C*8 MOTS(3) = 'KDIFF3 ' 55 : MOTS(1) = 'KD1 ' 56 : MOTS(2) = 'KD2 ' 57 : MOTS(3) = 'KD3 ' 58 : MOTS(4) = 'V1X ' 59 : MOTS(5) = 'V1Y ' 60 : ENDIF 61 : C--- Elements massifs TRIDimensionnels 62 : ELSE IF (IDIM.EQ.3) THEN 63 : JGM = 9 64 : SEGINI,MLMOTS 65 : C*8 MOTS(1) = 'KDIFF1 ' 66 : C*8 MOTS(2) = 'KDIFF2 ' 67 : C*8 MOTS(3) = 'KDIFF3 ' 68 : MOTS(1) = 'KD1 ' 69 : MOTS(2) = 'KD2 ' 70 : MOTS(3) = 'KD3 ' 71 : MOTS(4) = 'V1X ' 72 : MOTS(5) = 'V1Y ' 73 : MOTS(6) = 'V1Z ' 74 : MOTS(7) = 'V2X ' 75 : MOTS(8) = 'V2Y ' 76 : MOTS(9) = 'V2Z ' 77 : ENDIF 78 : 79 : C- Diffusion anisotrope 80 : ELSE IF (ITROP.EQ.3) THEN 81 : C--- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE 82 : IF (IDIM.EQ.2) THEN 83 : IF (IFOMOD.NE.1) THEN 84 : JGM = 5 85 : SEGINI,MLMOTS 86 : C*8 MOTS(1) = 'KDIFF11 ' 87 : C*8 MOTS(2) = 'KDIFF22 ' 88 : C*8 MOTS(3) = 'KDIFF21 ' 89 : MOTS(1) = 'KD11 ' 90 : MOTS(2) = 'KD22 ' 91 : MOTS(3) = 'KD21 ' 92 : MOTS(4) = 'V1X ' 93 : MOTS(5) = 'V1Y ' 94 : C--- Elements massifs bidimensionnels FOURIER 95 : ELSE 96 : JGM = 6 97 : SEGINI,MLMOTS 98 : MOCOMP=MLMOTS 99 : C*8 MOTS(1) = 'KDIFF11 ' 100 : C*8 MOTS(2) = 'KDIFF22 ' 101 : C*8 MOTS(3) = 'KDIFF21 ' 102 : C*8 MOTS(4) = 'KDIFF33 ' 103 : MOTS(1) = 'KD11 ' 104 : MOTS(2) = 'KD22 ' 105 : MOTS(3) = 'KD21 ' 106 : MOTS(4) = 'KD33 ' 107 : MOTS(5) = 'V1X ' 108 : MOTS(6) = 'V1Y ' 109 : ENDIF 110 : C--- Elements massifs TRIDimensionnels 111 : ELSE IF (IDIM.EQ.3) THEN 112 : JGM = 12 113 : SEGINI,MLMOTS 114 : C*8 MOTS( 1) = 'KDIFF11 ' 115 : C*8 MOTS( 2) = 'KDIFF22 ' 116 : C*8 MOTS( 3) = 'KDIFF33 ' 117 : C*8 MOTS( 4) = 'KDIFF21 ' 118 : C*8 MOTS( 5) = 'KDIFF31 ' 119 : C*8 MOTS( 6) = 'KDIFF32 ' 120 : MOTS( 1) = 'KD11 ' 121 : MOTS( 2) = 'KD22 ' 122 : MOTS( 3) = 'KD33 ' 123 : MOTS( 4) = 'KD21 ' 124 : MOTS( 5) = 'KD31 ' 125 : MOTS( 6) = 'KD32 ' 126 : MOTS( 7) = 'V1X ' 127 : MOTS( 8) = 'V1Y ' 128 : MOTS( 9) = 'V1Z ' 129 : MOTS(10) = 'V2X ' 130 : MOTS(11) = 'V2Y ' 131 : MOTS(12) = 'V2Z ' 132 : ENDIF 133 : ENDIF 134 : 135 : MOCOMP = MLMOTS 136 : IF (MOCOMP.EQ.0) THEN 136 : CALL ERREUR(5) 137 : RETURN 138 : ENDIF 139 : NBROBL = JGM 140 : 141 : IF (ITYPE.NE.0) THEN 142 : SEGINI,NOMID 143 : DO i = 1, NBROBL 144 : LESOBL(i) = MOTS(i) 145 : ENDDO 146 : SEGSUP,MLMOTS 147 : C* SEGDES,NOMID 148 : MOCOMP = NOMID 149 : ENDIF 150 : 151 : RETURN 152 : END 153 : 154 : 155 :
© Cast3M 2003 - All rights reserved.
Disclaimer