idmatr
C IDMATR SOURCE OF166741 24/05/06 21:15:19 11082 *--------------------------------------------------------------------* * Noms de composantes de materiaux * *--------------------------------------------------------------------* * * * ENTREES: * * MFR Numero de formulation * * IPMODL objet modele elementaire ( segment actif ) * * * * SORTIES: * * IPNOMC pointeur sur les listes de noms de composantes * * obligatoires et facultatives * * NBROBL nombre de composantes obligatoires * * NBRFAC nombre de composantes facultatives * * * *--------------------------------------------------------------------* IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC PPARAM -INC CCOPTIO C==DEB= FORMULATION HHO == Include specifique ========================== -INC CCHHOPA C==FIN= FORMULATION HHO ================================================ -INC SMLMOTS -INC SMMODEL POINTEUR NOMID1.NOMID ******Modif SELLIER **************************************************** * augmentation nombre de caracteristiques PARAMETER (ITA=500) CHARACTER*16 MOMODL(500) integer I,J,K,L c nombre de variables et de parametres par variables c include './nombre_helmholtz.h' -INC HNBRHEL CHARACTER*1 motnl1,motnlv1,motnld1,motnla1 CHARACTER*2 motnl2 CHARACTER*3 motnl3 CHARACTER*4 motnl4 ******fin modif SELLIER ************************************************ CHARACTER*(LOCOMP) TABOBL(ITA),TABFAC(ITA) IMODEL=IPMODL * Le segment existe-t-il deja? IF (lnomid(6).NE.0) THEN nomid = lnomid(6) SEGACT,nomid nbrobl = lesobl(/2) nbrfac = lesfac(/2) IF (nbrobl+nbrfac.EQ.0) GOTO 765 SEGINI,nomid1=nomid ipnomc=nomid1 RETURN ENDIF 765 CONTINUE JGOBL=0 JGFAC=0 ipnomc=0 IRET = 1 * NMAT=MATMOD(/2) NFOR=FORMOD(/2) MELE=NEFMOD C== FORMULATION HHO == MECANIQUE + idem qu'en MASSIF ==> MFR = 1 ======= C== Les composantes specifiques seront ajoutees dans HHOIDC via INOMID. IF (MFR.EQ.HHO_MFR_ELEMENT) MFR = 1 C== FORMULATION HHO ==================================================== *-------------------------------------------------------------------- * CAS DE LA FORMULATION THERMOHYDRIQUE *-------------------------------------------------------------------- IF (ithehy.ne.0) then if ( matmod(1).eq.'SCHREFLER') then JGOBL = 20 TABOBL(1)='KGG' TABOBL(2)='KGC' TABOBL(3)='KGT' TABOBL(4)='KCG' TABOBL(5)='KCC' TABOBL(6)='KCT' TABOBL(7)='KTG' TABOBL(8)='KTC' TABOBL(9)='KTT' TABOBL(10)='CGG' TABOBL(11)='CGC' TABOBL(12)='CGT' TABOBL(13)='CCG' TABOBL(14)='CCC' TABOBL(15)='CCT' TABOBL(16)='CTG' TABOBL(17)='CTC' TABOBL(18)='CTT' TABOBL(19)='KTGG' TABOBL(20)='KTCG' else IRET = 0 endif GO TO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION LIQUIDE ET ELEMENT DE RACCORD LITU *-------------------------------------------------------------------- IF (ILIQU.NE.0) THEN IF (MFR.NE.41.AND.MFR.NE.43) THEN JGOBL = 6 TABOBL(1)='RHO ' TABOBL(2)='CSON' TABOBL(3)='RORF' TABOBL(4)='CREF' TABOBL(5)='LCAR' TABOBL(6)='G ' ELSEIF (MFR.EQ.41) THEN JGOBL = 5 TABOBL(1)='RHO ' TABOBL(2)='CSON' TABOBL(3)='RORF' TABOBL(4)='CREF' TABOBL(5)='LCAR' c* ELSEIF (MFR.EQ.43) THEN ELSE JGOBL = 3 TABOBL(1)='RHO ' TABOBL(2)='LCAR' TABOBL(3)='RORF' ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION THERMIQUE * CAS DE LA FORMULATION DARCY *-------------------------------------------------------------------- IF (ITHER.NE.0 .OR. IDARCY.NE.0) THEN IF (ITHER.NE.0) THEN C Cas particuliers de THERMIQUE CONVECTION : IF (ICONV.NE.0) THEN C IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN CC Coques thermiques C JGOBL =2 C TABOBL(1)='HINF' C TABOBL(2)='HSUP' C JGFAC =2 C TABFAC(1)='TCINF' C TABFAC(2)='TCSUP' C else C Elements massifs JGOBL = 1 TABOBL(1)='H' JGFAC=1 TABFAC(1)='TC' C endif GOTO 9999 ENDIF C Cas particuliers de THERMIQUE RAYONNEMENT : IF (IRAYE.NE.0) THEN C Cas particuliers de THERMIQUE RAYONNEMENT : IF (ICAVE.NE.0) THEN C RAYONNEMENT en CAVITE : IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN JGOBL=2 TABOBL(1)='EINF' TABOBL(2)='ESUP' else JGOBL=1 TABOBL(1)='EMIS' endif JGFAC=3 TABFAC(1)='CABS' TABFAC(2)='TABS' TABFAC(3)='H' ELSEIF (IFACA.NE.0) THEN C RAYONNEMENT FACE_A_FACE : JGOBL=1 TABOBL(1)='EMIS' JGFAC = 1 TABFAC(1)='H' ELSEIF (IINFI.NE.0) THEN C RAYONNEMENT a l'INFINI : IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN JGOBL=3 TABOBL(1)='EINF' TABOBL(2)='ESUP' TABOBL(3)='T_IN' ELSE JGOBL=2 TABOBL(1)='EMIS' TABOBL(2)='T_IN' ENDIF JGFAC = 2 TABFAC(1)='E_IN' TABFAC(2)='H' ELSE IRET = 0 ENDIF GOTO 9999 ENDIF C Cas particuliers de THERMIQUE SOURCE : ISRCE = 0 IF (ISRCE.NE.0) THEN IF (INATUU.EQ.1.AND.IMATEE.EQ.1) THEN JGOBL = 1 TABOBL(1)='QVOL' JGFAC = 2 TABFAC(1)='QINF' TABFAC(2)='QSUP' GOTO 9999 ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.1) THEN JGOBL = 3 TABOBL(1)='QTOT' TABOBL(2)='ORIG' TABOBL(3)='RGAU' GOTO 9999 ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.2) THEN JGOBL = 5 TABOBL(1)='QTOT' TABOBL(2)='ORIG' TABOBL(3)='RGAU' TABOBL(4)='DIRE' TABOBL(5)='ZGAU' GOTO 9999 ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.3) THEN JGOBL = 7 TABOBL(1)='QTOT' TABOBL(2)='ORIG' TABOBL(3)='RGAU' TABOBL(4)='DIRE' TABOBL(5)='ZGAU' TABOBL(6)='DIRL' TABOBL(7)='LGAU' GOTO 9999 ELSE c write(6,*) 'INATUU, IMATEE =',INATUU, IMATEE WRITE(6,*) ' Dans IDMATR : numero IMATEE non prevu' RETURN ENDIF ENDIF ENDIF C* Cas THERMIQUE et DARCY IF (IORTH.EQ.0.AND.IANIS.EQ.0)THEN JGOBL = 1 IF(MFR .EQ. 75)THEN C Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES C ==================== TABOBL(1)='KT' ELSE TABOBL(1)='K ' ENDIF ELSEIF (IORTH.NE.0) THEN * ELEMENTS COQUES IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN JGOBL=5 TABOBL(1)='K1 ' TABOBL(2)='K2 ' TABOBL(3)='K3 ' TABOBL(4)='V1X ' TABOBL(5)='V1Y ' ELSEIF (MFR.EQ.1.OR.MFR.EQ.45) THEN * ELEMENTS MASSIFS et HYBRIDES IF (IDIM.EQ.2) THEN IF (IFOMOD.EQ.1) THEN * ELEMENT MASSIF DE FOURIER JGOBL=5 TABOBL(1)='K1 ' TABOBL(2)='K2 ' TABOBL(3)='V1X ' TABOBL(4)='V1Y ' TABOBL(5)='K3 ' ELSE JGOBL=4 TABOBL(1)='K1 ' TABOBL(2)='K2 ' TABOBL(3)='V1X ' TABOBL(4)='V1Y ' ENDIF ELSEIF (IDIM.EQ.3) THEN JGOBL=9 TABOBL(1)='K1 ' TABOBL(2)='K2 ' TABOBL(3)='K3 ' TABOBL(4)='V1X ' TABOBL(5)='V1Y ' TABOBL(6)='V1Z ' TABOBL(7)='V2X ' TABOBL(8)='V2Y ' TABOBL(9)='V2Z ' ENDIF ENDIF ELSEIF (IANIS.NE.0) THEN * ELEMENTS MASSIFS IF (MFR.EQ.1.OR.MFR.EQ.45) THEN IF (IDIM.EQ.2) THEN IF (IFOMOD.EQ.1) THEN * ELEMENT MASSIF DE FOURIER JGOBL=6 TABOBL(1)='K11 ' TABOBL(2)='K22 ' TABOBL(3)='K21 ' TABOBL(4)='V1X ' TABOBL(5)='V1Y ' TABOBL(6)='K33 ' ELSE JGOBL=5 TABOBL(1)='K11 ' TABOBL(2)='K22 ' TABOBL(3)='K21 ' TABOBL(4)='V1X ' TABOBL(5)='V1Y ' ENDIF ELSEIF (IDIM.EQ.3) THEN JGOBL=12 TABOBL(1)='K11 ' TABOBL(2)='K22 ' TABOBL(3)='K33 ' TABOBL(4)='K21 ' TABOBL(5)='K31 ' TABOBL(6)='K32 ' TABOBL(7)='V1X ' TABOBL(8)='V1Y ' TABOBL(9)='V1Z ' TABOBL(10)='V2X ' TABOBL(11)='V2Y ' TABOBL(12)='V2Z ' ENDIF ENDIF ENDIF C* C* Cas THERMIQUE CONDUCTION, THERMIQUE PHASE ou THERMIQUE ADVECTION : IF (ITHER.NE.0) THEN IF (IPHA.NE.0) THEN JGOBL0 = JGOBL JGOBL= JGOBL0+4 TABOBL(JGOBL0+1)='RHO ' TABOBL(JGOBL0+2)='C ' TABOBL(JGOBL0+3)='QLAT' TABOBL(JGOBL0+4)='TPHA' JGFAC=1 TABFAC(1)='H ' GOTO 9999 ENDIF if (iadve .ne. 0) then C Cas des Tuyaux 1D (MFR=79) C =========================== IF (MFR .EQ. 79) THEN JGFAC = 3 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITE' C Cas Massif (MFR=1) C =========================== ELSEIF (MFR .EQ. 1) THEN C Bidimensionnel PLAN DPGE C =========================== IF (IFOUR.EQ.-3) THEN JGFAC = 5 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITY' TABFAC(5)='VITZ' C Bidimensionnel PLAN (CP/DP) C =========================== ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN JGFAC = 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITY' C Axisymetrie C =========================== ELSEIF (IFOUR .EQ. 0) THEN JGFAC = 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITR' TABFAC(4)='VITZ' C Fourier C =========================== ELSEIF (IFOUR .EQ. 1) THEN RETURN C JGFAC = 5C C TABFAC(1)='RHO ' C TABFAC(2)='C ' C TABFAC(3)='VITR' C TABFAC(4)='VITZ' C TABFAC(5)='VITT' C Tridimensionnel C =========================== ELSEIF (IFOUR .EQ. 2) THEN JGFAC = 5 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITY' TABFAC(5)='VITZ' C Unidimensionnel (1D) C =========================== ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN IF (IFOUR.LE.6) THEN JGFAC = 3 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN JGFAC = 4 C Verifier l'utilite des composantes au dela de 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITY' ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN JGFAC = 4 C Verifier l'utilite des composantes au dela de 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITZ' ELSEIF (IFOUR.EQ.11) THEN JGFAC = 5 C Verifier l'utilite des composantes au dela de 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITX' TABFAC(4)='VITY' TABFAC(5)='VITZ' ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN JGFAC = 3 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITR' ELSEIF (IFOUR.EQ.14) THEN JGFAC = 4 C Verifier l'utilite des composantes au dela de 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='VITR' TABFAC(4)='VITZ' ENDIF ENDIF ENDIF GOTO 9999 endif C composantes facultatives THERMIQUE CONDUCTION seule IF(MFR .EQ. 75)THEN C Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES C ==================== JGFAC = 4 TABFAC(1)='M' TABFAC(2)='C' TABFAC(3)='TINI' ELSE JGFAC = 4 TABFAC(1)='RHO ' TABFAC(2)='C ' TABFAC(3)='H ' TABFAC(4)='TINI' ENDIF ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION MECANIQUE *-------------------------------------------------------------------- IF (IMECA.NE.0) THEN * IF (IRET.EQ.0) GOTO 9999 IF (NMAT.GE.2) THEN if ((IPLAC.EQ.9.OR.IPLAC.EQ.10).AND.NMAT.GT.2) then INMAT = 3 goto 19 endif IF (IPLAC.NE.0) THEN INMAT=3 IF (IPLAC.EQ.1) THEN INMAT=3 GOTO 10 ENDIF IF (IPLAC.EQ.2) THEN * * MATERIAU ELASTIQUE ORTHOTROPE * IF (MFR.EQ.75) THEN * JOINT UNIDIMENSIONNEL JOI1 * IF(IFOUR.EQ.2)THEN JGOBL=12 TABOBL(1)='V1X ' TABOBL(2)='V1Y ' TABOBL(3)='V1Z ' TABOBL(4)='V2X ' TABOBL(5)='V2Y ' TABOBL(6)='V2Z ' TABOBL(7)='KN ' TABOBL(8)='KS1 ' TABOBL(9)='KS2 ' TABOBL(10)='QN ' TABOBL(11)='QS1 ' TABOBL(12)='QS2 ' * JGFAC=10 TABFAC(1)='MASS' TABFAC(2)='JX ' TABFAC(3)='JY ' TABFAC(4)='JZ ' TABFAC(5)='ALPN' TABFAC(6)='ALP1' TABFAC(7)='ALP2' TABFAC(8)='ALQN' TABFAC(9)='ALQ1' TABFAC(10)='ALQ2' * ELSEIF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1)THEN JGOBL=5 TABOBL(1)='V1X ' TABOBL(2)='V1Y ' TABOBL(3)='KN ' TABOBL(4)='KS ' TABOBL(5)='QS ' * JGFAC=6 TABFAC(1)='MASS' TABFAC(2)='JZ' TABFAC(4)='ALPN' TABFAC(5)='ALPS' TABFAC(6)='ALQS' ENDIF * ELSEIF (MFR.EQ.3) THEN * COQUES MINCES * JGOBL=6 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='NU12' TABOBL(4)='G12 ' TABOBL(5)='V1X ' TABOBL(6)='V1Y ' * IF(IFOUR.EQ.-2) THEN JGFAC=4 ELSE JGFAC=3 ENDIF TABFAC(1)='ALP1' TABFAC(2)='ALP2' TABFAC(3)='RHO ' IF(IFOUR.EQ.-2) TABFAC(4)='DIM3' ELSEIF (MFR.EQ.9.OR.MFR.EQ.5) THEN * COQUES AVEC CISAILLEMENT TRANSVERSE * JGOBL=8 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='NU12' TABOBL(4)='G12 ' TABOBL(5)='G23 ' TABOBL(6)='G13 ' TABOBL(7)='V1X ' TABOBL(8)='V1Y ' * JGFAC=3 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' ELSEIF (MFR.EQ.1.OR.MFR.EQ.31) THEN * ELEMENTS MASSIFS * IF(IDIM.EQ.3)THEN * ELEMENTS 3D JGOBL=15 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='G23 ' TABOBL(9)='G13 ' TABOBL(10)='V1X ' TABOBL(11)='V1Y ' TABOBL(12)='V1Z ' TABOBL(13)='V2X ' TABOBL(14)='V2Y ' TABOBL(15)='V2Z ' * JGFAC=4 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' ELSEIF (IDIM.EQ.2) THEN IF(IFOUR.EQ.-2)THEN * CONTRAINTE PLANE JGOBL=9 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='NU12' TABOBL(4)='G12 ' TABOBL(5)='V1X ' TABOBL(6)='V1Y ' TABOBL(7)='YG3 ' TABOBL(8)='NU23' TABOBL(9)='NU13' * JGFAC=4 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='DIM3' ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN * DEFORMATION PLANE ,AXISYMETRIE JGOBL=9 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='V1X ' TABOBL(9)='V1Y ' * JGFAC=4 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' ELSEIF (IFOUR.EQ.1) THEN * AXISYMETRIE DE FOURIER JGOBL=15 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='G23 ' TABOBL(9)='G13 ' TABOBL(10)='V1X ' TABOBL(11)='V1Y ' TABOBL(12)='V1Z ' TABOBL(13)='V2X ' TABOBL(14)='V2Y ' TABOBL(15)='V2Z ' * JGFAC=4 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' ENDIF ELSEIF (IDIM.EQ.1) THEN C= Dans le cas UNID SPHErique, on doit avoir YG2=YG3 et NU12=NU13 C= et dans le cas thermomecanique ALP2=ALP3. JGOBL=6 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' IF (IFOUR.EQ.6) THEN JGFAC=2 ELSEIF (IFOUR.EQ.5 .OR. IFOUR.EQ.10) THEN JGFAC=3 TABFAC(3)='ALP3' ELSEIF (IFOUR.EQ.4.OR.IFOUR.EQ.8.OR. & IFOUR.EQ.13) THEN JGFAC=3 TABFAC(3)='ALP2' ELSE JGFAC=4 TABFAC(3)='ALP2' TABFAC(4)='ALP3' ENDIF TABFAC(1)='RHO ' TABFAC(2)='ALP1' ENDIF ELSEIF (MFR.EQ.35) THEN * ELEMENTS JOINTS IF (IFOUR.EQ.2) THEN JGOBL=5 TABOBL(1)='KS1 ' TABOBL(2)='KS2 ' TABOBL(3)='KN ' TABOBL(4)='V1X ' TABOBL(5)='V1Y ' * JGFAC=2 TABFAC(1)='RHO ' TABFAC(2)='ALPN' ENDIF ENDIF ELSEIF (IPLAC.EQ.3)THEN * MATERIAU ANISOTROPE ELASTIQUE * IF(MFR.EQ.75)THEN * JOINT UNIDIMESIONNEL JOI1 * IF (IFOUR.EQ.2) THEN JGOBL=27 TABOBL(1)='V1X ' TABOBL(2)='V1Y ' TABOBL(3)='V1Z ' TABOBL(4)='V2X ' TABOBL(5)='V2Y ' TABOBL(6)='V2Z ' TABOBL(7)='D11 ' TABOBL(8)='D22 ' TABOBL(9)='D33 ' TABOBL(10)='D44 ' TABOBL(11)='D55 ' TABOBL(12)='D66 ' TABOBL(13)='D21 ' TABOBL(14)='D31 ' TABOBL(15)='D32 ' TABOBL(16)='D41 ' TABOBL(17)='D42 ' TABOBL(18)='D43 ' TABOBL(19)='D51 ' TABOBL(20)='D52 ' TABOBL(21)='D53 ' TABOBL(22)='D54 ' TABOBL(23)='D61 ' TABOBL(24)='D62 ' TABOBL(25)='D63 ' TABOBL(26)='D64 ' TABOBL(27)='D65 ' * JGFAC=10 * MASS: masse totale de l'élément joint TABFAC(1)='MASS' TABFAC(2)='JX ' TABFAC(3)='JY ' TABFAC(4)='JZ ' TABFAC(5)='ALP1' TABFAC(6)='ALP2' TABFAC(7)='ALP3' TABFAC(8)='ALQ1' TABFAC(9)='ALQ2' TABFAC(10)='ALQ3' * ELSEIF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN JGOBL=8 TABOBL(1)='V1X ' TABOBL(2)='V1Y ' TABOBL(3)='D11 ' TABOBL(4)='D22 ' TABOBL(5)='D33 ' TABOBL(6)='D21 ' TABOBL(7)='D31 ' TABOBL(8)='D32 ' * JGFAC=5 * MASS: masse totale de l'élément joint TABFAC(1)='MASS' TABFAC(2)='JZ ' TABFAC(3)='ALP1' TABFAC(4)='ALP2' TABFAC(5)='ALQ3' ENDIF ENDIF * * ELEMENTS MASSIFS * IF(MFR.EQ.1.OR.MFR.EQ.31)THEN IF(IDIM.EQ.3)THEN * ELEMENTS 3D JGOBL=27 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='D51 ' TABOBL(12)='D52 ' TABOBL(13)='D53 ' TABOBL(14)='D54 ' TABOBL(15)='D55 ' TABOBL(16)='D61 ' TABOBL(17)='D62 ' TABOBL(18)='D63 ' TABOBL(19)='D64 ' TABOBL(20)='D65 ' TABOBL(21)='D66 ' TABOBL(22)='V1X ' TABOBL(23)='V1Y ' TABOBL(24)='V1Z ' TABOBL(25)='V2X ' TABOBL(26)='V2Y ' TABOBL(27)='V2Z ' * JGFAC=7 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' TABFAC(5)='AL12' TABFAC(6)='AL13' TABFAC(7)='AL23' ELSEIF (IDIM.EQ.2) THEN IF (IFOUR.EQ.-2) THEN * CONTRAINTE PLANE JGOBL=12 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D41 ' TABOBL(5)='D42 ' TABOBL(6)='D44 ' TABOBL(7)='V1X ' TABOBL(8)='V1Y ' TABOBL(9)='D31 ' TABOBL(10)='D32 ' TABOBL(11)='D33 ' TABOBL(12)='D43 ' * JGFAC=5 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='DIM3' ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN * DEFORMATION PLANE ,AXISYMETRIE JGOBL=12 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='V1X ' TABOBL(12)='V1Y ' * JGFAC=5 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='ALP3' ELSEIF (IFOUR.EQ.1) THEN * AXISYMETRIE DE FOURIER JGOBL=15 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='D55 ' TABOBL(12)='D65 ' TABOBL(13)='D66 ' TABOBL(14)='V1X ' TABOBL(15)='V1Y ' * JGFAC=5 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='ALP3' ENDIF ENDIF ENDIF ELSEIF (IPLAC.EQ.4) THEN IF (MFR.EQ.33) THEN * MILIEU POREUX ISOTROPE * JGOBL=12 TABOBL(1) ='YOUN' TABOBL(2) ='NU ' TABOBL(3) ='RHO ' TABOBL(4) ='ALPH' TABOBL(5) ='COB ' TABOBL(6) ='MOB ' TABOBL(7) ='ALPM' TABOBL(8) ='PERM' TABOBL(9) ='VISC' TABOBL(10)='KF ' TABOBL(11)='RHOF' TABOBL(12)='ALPF' * IF (IFOUR.EQ.-2) THEN JGFAC=1 TABFAC(1)='DIM3' ENDIF ENDIF ELSEIF (IPLAC.EQ.5) THEN IF (MFR.EQ.37) THEN * MILIEU HOMOGENEISE * JGOBL=16 TABOBL( 1)='B11 ' TABOBL( 2)='B22 ' TABOBL( 3)='B12 ' TABOBL( 4)='ROF ' TABOBL( 5)='ROS ' TABOBL( 6)='YOUN' TABOBL( 7)='CSON' TABOBL( 8)='RORF' TABOBL( 9)='CREF' TABOBL(10)='LCAR' TABOBL(11)='E111' TABOBL(12)='E112' TABOBL(13)='E121' TABOBL(14)='E122' TABOBL(15)='E221' TABOBL(16)='E222' ENDIF ELSEIF (IPLAC.EQ.6) THEN C C MATERIAU ELASTIQUE UNIDIRECTIONNEL C IF ((MFR.EQ.1.OR.MFR.EQ.31).AND.IDIM.EQ.3) THEN JGOBL=7 TABOBL(1)='YOUN' TABOBL(2)='V1X ' TABOBL(3)='V1Y ' TABOBL(4)='V1Z ' TABOBL(5)='V2X ' TABOBL(6)='V2Y ' TABOBL(7)='V2Z ' ELSE JGOBL=3 TABOBL(1)='YOUN' TABOBL(2)='V1X ' TABOBL(3)='V1Y ' ENDIF IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND. . IFOUR.EQ.-2) THEN JGFAC=3 TABFAC(3)='DIM3' ELSE JGFAC=2 ENDIF TABFAC(1)='RHO ' TABFAC(2)='ALPH' ELSEIF (IPLAC.EQ.7) THEN C MODELE ET MATERIAU DE LA SECTION DU MODELE A FIBRE C JGOBL=2 TABOBL(1)='MODS' TABOBL(2)='MATS' JGFAC=1 TABFAC(1)='MANO' C ---- IPLAC 8: Ajoute par Jiang, 22/08/1995 et fleuret 28/05/96 ELSEIF (IPLAC.EQ.8) THEN C MODELE ET MATERIAU POUR LE CABLE PRECONTRAINT C JGOBL=1 TABOBL(1)='YOUN' JGFAC=8 TABFAC(1)='FF' TABFAC(2)='PHIF' TABFAC(3)='GANC' TABFAC(4)='RMU0' TABFAC(5)='FPRG' TABFAC(6)='RH10' TABFAC(7)='ALPH' TABFAC(8)='RHO ' ELSEIF (IPLAC.EQ.9) THEN C MODAL C JGOBL=3 TABOBL(1)='FREQ' TABOBL(2)='MASS' TABOBL(3)='DEFO' JGFAC = 9 TABFAC(1) = 'AMOR' TABFAC(2) = 'CGRA' TABFAC(3) = 'MADE' TABFAC(4) = 'RICR' TABFAC(5) = 'MAIB' TABFAC(6) = 'MACR' TABFAC(7) = 'AMCR' TABFAC(8) = 'ALP0' TABFAC(9) = 'ECRO' ELSEIF (IPLAC.EQ.10) THEN C STATIQUE C JGOBL=3 TABOBL(1)='DEFO' TABOBL(2)='RIDE' TABOBL(3)='MADE' JGFAC = 8 TABFAC(1) = 'AMOR' TABFAC(2) = 'RICR' TABFAC(3) = 'MAIA' TABFAC(4) = 'MAIB' TABFAC(5) = 'MACR' TABFAC(6) = 'AMCR' TABFAC(7) = 'BET0' TABFAC(8) = 'ECRO' ELSEIF (IPLAC.EQ.11) THEN C ZONE_COHESIVE C JGOBL=2 TABOBL(1)='KS' TABOBL(2)='KN' ENDIF INMAT=3 GOTO 20 ELSE INMAT=2 GOTO 10 ENDIF ENDIF INMAT=0 10 CONTINUE * * MATERIAU ELASTIQUE ISOTROPE * IF (MFR.EQ.35) THEN * VALABLE EN 2D COMME EN 3D JGOBL=2 TABOBL(1)='KS ' TABOBL(2)='KN ' JGFAC=2 TABFAC(1)='RHO ' TABFAC(2)='ALPN' * ELSEIF (MFR.EQ.78) THEN * VALABLE EN 2D COMME EN 3D JGOBL=2 TABOBL(1)='KS ' TABOBL(2)='KN ' * * JOINT CISAILLEMENT (2D) * ELSEIF (MFR.EQ.53) THEN * VALABLE EN 2D "COMME EN 3D" JGOBL=1 TABOBL(1)='KS ' * JGFAC=2 TABFAC(1)='RHO ' TABFAC(2)='ALPN' * * TOUS LES AUTRES CAS * ELSE JGOBL=2 TABOBL(1)='YOUN' TABOBL(2)='NU ' * IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.63).AND. . IFOUR.EQ.-2) THEN JGFAC=4 TABFAC(3)='DIM3' TABFAC(4)='VISQ' ELSE JGFAC=3 TABFAC(3)='VISQ' ENDIF TABFAC(1)='RHO ' TABFAC(2)='ALPH' * AM CAS FORMULATION NON-LOCALE INLOC = 0 MN3 = INFMOD(/1) IF (MN3.GE.13) INLOC=-1*INFMOD(13) IF (INLOC.GT.0) THEN JGOBL0 = JGFAC * moyenne IF(INLOC.EQ.1) THEN JGFAC=JGOBL0+1 TABFAC(JGOBL0+1)='LCAR' * stress-based ELSEIF(INLOC.EQ.2) THEN JGFAC=JGOBL0+2 TABFAC(JGOBL0+1)='LCAR' TABFAC(JGOBL0+2)='SBFT' * helmholtz ELSE IF(INLOC.EQ.3) THEN * sellier 03/07/20 * JGFAC=JGOBL0+6 * TABFAC(JGOBL0+1)='LCAR' * TABFAC(JGOBL0+2)='LCF1' * TABFAC(JGOBL0+3)='LCF2' * TABFAC(JGOBL0+4)='LCF3' * TABFAC(JGOBL0+5)='LCF4' * TABFAC(JGOBL0+6)='LCF5' JGFAC=JGOBL0+NB_PARA_HELM c variables de Helmholtz a diffusions ORThotropes c modif liste para heilmholtz sellier 10//11/22 c print*,'Dans IDMATR ', NB_HELM, 'formulationss Helmholtz' do I=1,NB_HELM c numero de la variable de Helmholtz write (motnl1,'(I1)') I c Capacité (infinie si BH=1) motnl3='CAP' motnl4=motnl3//motnl1 TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+1)=motnl4 c print*,motnl4 c Blocage de Dirichlet si 1 sinon 0 motnl3='BLO' motnl4=motnl3//motnl1 TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+2)=motnl4 c print*,motnl4 c Valeur imposee a la zone de blocage de Dirichlet motnl3='DEP' motnl4=motnl3//motnl1 TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+3)=motnl4 c print*,motnl4 c Valeur initiale de la variable de Helmholtz motnl3='INI' motnl4=motnl3//motnl1 TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+4)=motnl4 c print*,motnl4 c indicateur logique de lineratite pour ne pas passer par istep=3 si 1 motnl3='LIN' motnl4=motnl3//motnl1 TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+5)=motnl4 c print*,motnl4 do J=1,3 c coefficients de diffusion write (motnld1,'(I1)') J c diffusion motnl2='DH' c 1er indice: numero de la variable d helmholtz motnl3=motnl2//motnl1 c 2eme indice : direction motnl4=motnl3//motnld1 TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+ # 5+(J-1)*4+1)=motnl4 c print*,motnl4 do K=1,3 write (motnla1,'(I1)') K c direction principales de diffusion motnlv1='V' c indice numero de la variable motnl2=motnlv1//motnl1 c indice direction principale motnl3=motnl2//motnld1 c projection sur axe indice de base fixe motnl4=motnl3//motnla1 TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+ # 5+(J-1)*4+1+K)=motnl4 c print*,motnl4 end do end do end do Cc on rajoute 2 inutiles pour compenser un decalage constatte dans castem20 C TABFAC(JGOBL0+NB_PARA_HELM+1)='FAC1' C TABFAC(JGOBL0+NB_PARA_HELM+2)='FAC2' * fin modif sellier * cas non prevu ELSE IRET = 0 ENDIF ENDIF * ENDIF * IF (INMAT.EQ.0) GOTO 9999 * 19 CONTINUE 20 CONTINUE DO jm = 1,matmod(/2) IF (matmod(jm).eq.'IMPEDANCE') THEN imate = imatee INMAT = INMAT + 1 JGOBL = 0 JGFAC = 0 IF(CMATEE.EQ.'IMPELAST') THEN *IMPE_ELAS JGOBL0 = JGOBL JGOBL= JGOBL0+1 TABOBL(JGOBL0 + 1) ='RAID' JGOBL0 = JGFAC JGFAC = JGOBL0+4 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'ZNU' TABFAC(JGOBL0+3) = 'MASS' TABFAC(JGOBL0+4) = 'ALPH' GOTO 22 ELSEIF (CMATEE.EQ.'IMPVOIGT'.or.CMATEE.EQ.'IMPREUSS') THEN *IMPE_VOIGT ou IMPE_REUSS JGOBL0 = JGOBL JGOBL= JGOBL0+2 TABOBL(JGOBL0 + 1) ='RAID' TABOBL(JGOBL0 + 2) ='VISC' JGOBL0 = JGFAC JGFAC = JGOBL0+2 TABFAC(JGOBL0+1) = 'MASS' TABFAC(JGOBL0+2) = 'AMOR' GOTO 22 ELSEIF (CMATEE.EQ.'IMPCOMPL') then *IMPE_COMPLEXE JGOBL0 = JGOBL JGOBL= JGOBL0+1 TABOBL(JGOBL0 + 1) ='RAID' JGOBL0 = JGFAC JGFAC = JGFAC+4 TABFAC(JGOBL0+1) = 'MOCO' TABFAC(JGOBL0+2) = 'VISC' TABFAC(JGOBL0+3) = 'MASS' TABFAC(JGOBL0+4) = 'AMOR' GOTO 22 ELSE ENDIF 22 CONTINUE * if (mele.eq.45) then JGOBL0 = JGFAC JGFAC = JGOBL0+3 IF(CMATEE.EQ.'IMPELAST') THEN TABFAC(JGOBL0+1) = TABFAC(JGOBL0-2) TABFAC(JGOBL0+2) = TABFAC(JGOBL0-1) TABFAC(JGOBL0+3) = TABFAC(JGOBL0) TABFAC(JGOBL0-2) = 'CPLE' TABFAC(JGOBL0-1) = 'INER' TABFAC(JGOBL0) = 'AROT' ELSE TABFAC(JGOBL0+1) = 'CPLE' TABFAC(JGOBL0+2) = 'INER' TABFAC(JGOBL0+3) = 'AROT' ENDIF * endif IF (NMAT.GE.INMAT) inmat = inmat+1 * GOTO 9999 ENDIF ENDDO * IF (NMAT.GE.INMAT) THEN IF (IPLAC.EQ.1) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * ISOTROPE IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF ELSEIF (IPLAC.EQ.2) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * NORTON IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF ELSEIF (IPLAC.EQ.3) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * ONERA IPLAC=2 * ELSE * INMAT=INMAT+1 ENDIF C C Cas particulier des modeles GATT_MONERIE et UO2 C 'RHO ' et 'ALPH' sont obligatoires IF (IPLAC.EQ.18.OR.IPLAC.EQ.19) THEN JGOBL0 = JGOBL JGOBL=JGOBL0+2 TABOBL(JGOBL0+1)='RHO ' TABOBL(JGOBL0+2)='ALPH' JGOBL0=JGFAC JGFAC=JGOBL0-2 IF (JGFAC.GT.0) THEN DO 200 I=1,JGFAC TABFAC(I)=TABFAC(I+2) 200 CONTINUE ENDIF C === C Modeles SYMONDS & COWPER SYCO1 et SYCO2 C === ELSEIF (IPLAC.EQ.28) THEN JGOBL0=JGOBL JGOBL=JGOBL0+1 TABOBL(JGOBL0+1)='ECRO' ELSEIF (IPLAC.EQ.29) THEN JGOBL0=JGOBL JGOBL=JGOBL0+1 TABOBL(JGOBL0+1)='ECRO' ENDIF C ELSEIF (IPLAC.EQ.4) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * MAZARS IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF ELSEIF (IPLAC.EQ.5) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * TRIAXIAL P/Y IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF ELSEIF (IPLAC.EQ.6) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * ISOTROPE IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF ELSEIF (IPLAC.EQ.7) THEN C Pas de composantes supplementaires pour une loi 'VISCO_EXTERNE' IRET = 1 ELSEIF(IPLAC.EQ.0.AND.IMATEE.EQ.31) THEN GOTO 9999 ELSE GOTO 9999 ENDIF C IRET code retour de IDPLAS,IDFLUA,IDVISC,IDENDO,IDPLEN,IDELNL IF (IRET.EQ.0) GOTO 9999 ENDIF * Parametres CRIP et FUSION : do jma=1,matmod(/2) if(matmod(jma).eq.'CRIP') then JGOBL = JGOBL + 1 TABOBL(JGOBL) = 'LIMP' elseif(matmod(jma).eq.'FUSION') then JGOBL = JGOBL + 1 TABOBL(JGOBL) = 'TFUS' endif enddo * GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION POREUX *-------------------------------------------------------------------- IF (IPORE.NE.0) THEN * IF (IRET.EQ.0) GOTO 9999 * * D'ABORD : CAS NON ISOTROPE * IF (NMAT.GE.2) THEN IF (IPLAC.NE.0) THEN INMAT=3 IF (IPLAC.EQ.1) THEN INMAT=3 GOTO 30 ENDIF IF (IPLAC.EQ.2) THEN * * MATERIAU ELASTIQUE ORTHOTROPE * IF (MFR.EQ.33) THEN * * ELEMENTS MASSIFS * IF(IDIM.EQ.3)THEN * ELEMENTS 3D JGOBL=19 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='G23 ' TABOBL(9)='G13 ' TABOBL(10)='V1X ' TABOBL(11)='V1Y ' TABOBL(12)='V1Z ' TABOBL(13)='V2X ' TABOBL(14)='V2Y ' TABOBL(15)='V2Z ' TABOBL(16)='COB1' TABOBL(17)='COB2' TABOBL(18)='COB3' TABOBL(19)='MOB ' * JGFAC=12 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' TABFAC(5)='ALPM' TABFAC(6)='PER1' TABFAC(7)='PER2' TABFAC(8)='PER3' TABFAC(9)='VISC' TABFAC(10)='KF ' TABFAC(11)='RHOF' TABFAC(12)='ALPF' ELSEIF (IDIM.EQ.2) THEN IF(IFOUR.EQ.-2)THEN * CONTRAINTE PLANE JGOBL=12 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='NU12' TABOBL(4)='G12 ' TABOBL(5)='V1X ' TABOBL(6)='V1Y ' TABOBL(7)='COB1' TABOBL(8)='COB2' TABOBL(9)='MOB ' TABOBL(10)='YG3 ' TABOBL(11)='NU23' TABOBL(12)='NU13' * JGFAC=11 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALPM' TABFAC(5)='PER1' TABFAC(6)='PER2' TABFAC(7)='VISC' TABFAC(8)='KF ' TABFAC(9)='RHOF' TABFAC(10)='ALPF' TABFAC(11)='DIM3' ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN * DEFORMATION PLANE ,AXISYMETRIE JGOBL=13 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='V1X ' TABOBL(9)='V1Y ' TABOBL(10)='COB1' TABOBL(11)='COB2' TABOBL(12)='COB3' TABOBL(13)='MOB ' * JGFAC=11 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' TABFAC(5)='ALPM' TABFAC(6)='PER1' TABFAC(7)='PER2' TABFAC(8)='VISC' TABFAC(9)='KF ' TABFAC(10)='RHOF' TABFAC(11)='ALPF' * ELSEIF (IFOUR.EQ.1) THEN * AXISYMETRIE DE FOURIER JGOBL=15 TABOBL(1)='YG1 ' TABOBL(2)='YG2 ' TABOBL(3)='YG3 ' TABOBL(4)='NU12' TABOBL(5)='NU23' TABOBL(6)='NU13' TABOBL(7)='G12 ' TABOBL(8)='G23 ' TABOBL(9)='G13 ' TABOBL(10)='V1X ' TABOBL(11)='V1Y ' TABOBL(12)='COB1' TABOBL(13)='COB2' TABOBL(14)='COB3' TABOBL(15)='MOB ' * JGFAC=12 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' TABFAC(5)='ALPM' TABFAC(6)='PER1' TABFAC(7)='PER2' TABFAC(8)='PER3' TABFAC(9)='VISC' TABFAC(10)='KF ' TABFAC(11)='RHOF' TABFAC(12)='ALPF' ENDIF ENDIF ELSE * * CAS NON PREVU * IRET = 0 GOTO 9999 ENDIF ELSEIF (IPLAC.EQ.3)THEN * * MATERIAU ANISOTROPE ELASTIQUE * IF(MFR.EQ.33)THEN * * ELEMENTS MASSIFS * IF(IDIM.EQ.3)THEN * ELEMENTS 3D JGOBL=34 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='D51 ' TABOBL(12)='D52 ' TABOBL(13)='D53 ' TABOBL(14)='D54 ' TABOBL(15)='D55 ' TABOBL(16)='D61 ' TABOBL(17)='D62 ' TABOBL(18)='D63 ' TABOBL(19)='D64 ' TABOBL(20)='D65 ' TABOBL(21)='D66 ' TABOBL(22)='V1X ' TABOBL(23)='V1Y ' TABOBL(24)='V1Z ' TABOBL(25)='V2X ' TABOBL(26)='V2Y ' TABOBL(27)='V2Z ' TABOBL(28)='COB1' TABOBL(29)='COB2' TABOBL(30)='COB3' TABOBL(31)='CO12' TABOBL(32)='CO13' TABOBL(33)='CO23' TABOBL(34)='MOB ' * JGFAC=18 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='ALP3' TABFAC(5)='AL12' TABFAC(6)='AL13' TABFAC(7)='AL23' TABFAC(8)='ALPM' TABFAC(9 )='PER1' TABFAC(10)='PER2' TABFAC(11)='PER3' TABFAC(12)='PE12' TABFAC(13)='PE13' TABFAC(14)='PE23' TABFAC(15)='VISC' TABFAC(16)='KF ' TABFAC(17)='RHOF' TABFAC(18)='ALPF' ELSEIF (IDIM.EQ.2) THEN IF (IFOUR.EQ.-2) THEN * CONTRAINTE PLANE JGOBL=16 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D41 ' TABOBL(5)='D42 ' TABOBL(6)='D44 ' TABOBL(7)='V1X ' TABOBL(8)='V1Y ' TABOBL(9 )='COB1' TABOBL(10)='COB2' TABOBL(11)='CO12' TABOBL(12)='MOB ' TABOBL(13)='D31 ' TABOBL(14)='D32 ' TABOBL(15)='D33 ' TABOBL(16)='D43 ' * JGFAC=13 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='ALPM' TABFAC(6)='PER1' TABFAC(7)='PER2' TABFAC(8)='PE12' TABFAC(9)='VISC' TABFAC(10)='KF ' TABFAC(11)='RHOF' TABFAC(12)='ALPF' TABFAC(13)='DIM3' * ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN * DEFORMATION PLANE ,AXISYMETRIE JGOBL=17 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='V1X ' TABOBL(12)='V1Y ' TABOBL(13)='COB1' TABOBL(14)='COB2' TABOBL(15)='CO12' TABOBL(16)='COB3' TABOBL(17)='MOB ' * JGFAC=13 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='ALP3' TABFAC(6)='ALPM' TABFAC(7)='PER1' TABFAC(8)='PER2' TABFAC(9)='PE12' TABFAC(10)='VISC' TABFAC(11)='KF ' TABFAC(12)='RHOF' TABFAC(13)='ALPF' * ELSEIF (IFOUR.EQ.1) THEN * AXISYMETRIE DE FOURIER JGOBL=20 TABOBL(1)='D11 ' TABOBL(2)='D21 ' TABOBL(3)='D22 ' TABOBL(4)='D31 ' TABOBL(5)='D32 ' TABOBL(6)='D33 ' TABOBL(7)='D41 ' TABOBL(8)='D42 ' TABOBL(9)='D43 ' TABOBL(10)='D44 ' TABOBL(11)='D55 ' TABOBL(12)='D65 ' TABOBL(13)='D66 ' TABOBL(14)='V1X ' TABOBL(15)='V1Y ' TABOBL(16)='COB1' TABOBL(17)='COB2' TABOBL(18)='CO12' TABOBL(19)='COB3' TABOBL(20)='MOB ' * JGFAC=14 TABFAC(1)='RHO ' TABFAC(2)='ALP1' TABFAC(3)='ALP2' TABFAC(4)='AL12' TABFAC(5)='ALP3' TABFAC(6)='ALPM' TABFAC(7)='PER1' TABFAC(8)='PER2' TABFAC(9)='PE12' TABFAC(10)='PER3' TABFAC(11)='VISC' TABFAC(12)='KF ' TABFAC(13)='RHOF' TABFAC(14)='ALPF' ENDIF ENDIF ELSE * * CAS NON PREVU * IRET = 0 GOTO 9999 ENDIF ELSEIF (IPLAC.EQ.6) THEN C C MATERIAU ELASTIQUE UNIDIRECTIONNEL C IF(MFR.EQ.33)THEN IF (IDIM.EQ.3) THEN JGOBL=9 TABOBL(1)='YOUN' TABOBL(2)='V1X ' TABOBL(3)='V1Y ' TABOBL(4)='V1Z ' TABOBL(5)='V2X ' TABOBL(6)='V2Y ' TABOBL(7)='V2Z ' TABOBL(8)='COB ' TABOBL(9)='MOB ' ELSE JGOBL=5 TABOBL(1)='YOUN' TABOBL(2)='V1X ' TABOBL(3)='V1Y ' TABOBL(4)='COB ' TABOBL(5)='MOB ' ENDIF * IF(IFOUR.EQ.-2) THEN JGFAC=9 TABFAC(9)='DIM3' ELSE JGFAC=8 ENDIF TABFAC(1)='RHO ' TABFAC(2)='ALPH' TABFAC(3)='ALPM' TABFAC(4)='PERM' TABFAC(5)='VISC' TABFAC(6)='KF ' TABFAC(7)='RHOF' TABFAC(8)='ALPF' ELSE * * CAS NON PREVU * IRET = 0 GO TO 9999 ENDIF * ENDIF INMAT=3 GOTO 40 ELSE INMAT=2 GOTO 30 ENDIF ENDIF INMAT=0 30 CONTINUE * * CAS MATERIAU POREUX ELASTIQUE ISOTROPE * IF(MELE.GE.79.AND.MELE.LE.83)THEN * JGOBL=4 TABOBL(1) ='YOUN' TABOBL(2) ='NU ' TABOBL(3) ='COB ' TABOBL(4) ='MOB ' * IF(IFOUR.EQ.-2) THEN JGFAC=9 TABFAC(9)='DIM3' ELSE JGFAC=8 ENDIF TABFAC(1)='RHOF' TABFAC(2)='ALPF' TABFAC(3)='ALPM' TABFAC(4)='PERM' TABFAC(5)='VISC' TABFAC(6)='KF ' TABFAC(7)='RHO ' TABFAC(8)='ALPH' * ELSEIF(MELE.GE.108.AND.MELE.LE.110)THEN * * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D) * JGOBL=4 TABOBL(1)='KS ' TABOBL(2)='KN ' TABOBL(3)='COB ' TABOBL(4)='MOB ' * JGFAC=4 TABFAC(1)='PERT' TABFAC(2)='PERH' TABFAC(3)='PERB' TABFAC(4)='VISC' * ELSEIF(MELE.GE.173.AND.MELE.LE.177)THEN * JGOBL=10 TABOBL(1) ='YOUN' TABOBL(2) ='NU ' TABOBL(3) ='COP1' TABOBL(4) ='COP2' TABOBL(5) ='CPP1' TABOBL(6) ='CPP2' TABOBL(7) ='KK11' TABOBL(8) ='KK12' TABOBL(9) ='KK21' TABOBL(10)='KK22' * IF(IFOUR.EQ.-2) THEN JGFAC=9 TABFAC(9)='DIM3' ELSE JGFAC=8 ENDIF TABFAC(1)='RHOF' TABFAC(2)='ALPF' TABFAC(3)='ALPM' TABFAC(4)='PK11' TABFAC(5)='PK12' TABFAC(6)='PK21' TABFAC(7)='PK22' TABFAC(8)='ALPH' * ELSEIF(MELE.GE.185.AND.MELE.LE.187)THEN * * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D) * JGOBL=10 TABOBL(1)='KS ' TABOBL(2)='KN ' TABOBL(3)='COP1' TABOBL(4)='COP2' TABOBL(5)='CPP1' TABOBL(6)='CPP2' TABOBL(7)='KK11' TABOBL(8)='KK12' TABOBL(9)='KK21' TABOBL(10)='KK22' * IF(IFOUR.EQ.-2) THEN JGFAC=18 TABFAC(18)='DIM3' ELSE JGFAC=17 ENDIF TABFAC(1)='RHOF' TABFAC(2)='ALPF' TABFAC(3)='ALPM' TABFAC(4)='PT11' TABFAC(5)='PH11' TABFAC(6)='PB11' TABFAC(7)='PT12' TABFAC(8)='PH12' TABFAC(9)='PB12' TABFAC(10)='PT21' TABFAC(11)='PH21' TABFAC(12)='PB21' TABFAC(13)='PT22' TABFAC(14)='PH22' TABFAC(15)='PB22' TABFAC(16)='RHO ' TABFAC(17)='ALPH' ELSEIF(MELE.GE.178.AND.MELE.LE.182)THEN * JGOBL=17 TABOBL(1)='YOUN' TABOBL(2)='NU ' TABOBL(3)='COP1' TABOBL(4)='COP2' TABOBL(5)='COP3' TABOBL(6)='CPP1' TABOBL(7)='CPP2' TABOBL(8)='CPP3' TABOBL(9)='KK11' TABOBL(10)='KK12' TABOBL(11)='KK13' TABOBL(12)='KK21' TABOBL(13)='KK22' TABOBL(14)='KK23' TABOBL(15)='KK31' TABOBL(16)='KK32' TABOBL(17)='KK33' * IF(IFOUR.EQ.-2) THEN JGFAC=15 TABFAC(15)='DIM3' ELSE JGFAC=14 ENDIF TABFAC(1)='RHOF' TABFAC(2)='ALPF' TABFAC(3)='ALPM' TABFAC(4)='PK11' TABFAC(5)='PK12' TABFAC(6)='PK13' TABFAC(7)='PK21' TABFAC(8)='PK22' TABFAC(9)='PK23' TABFAC(10)='PK31' TABFAC(11)='PK32' TABFAC(12)='PK33' TABFAC(13)='RHO ' TABFAC(14)='ALPH' ELSEIF(MELE.GE.188.AND.MELE.LE.190)THEN * * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D) * JGOBL=17 TABOBL(1)='KS ' TABOBL(2)='KN ' TABOBL(3)='COP1' TABOBL(4)='COP2' TABOBL(5)='COP3' TABOBL(6)='CPP1' TABOBL(7)='CPP2' TABOBL(8)='CPP3' TABOBL(9)='KK11' TABOBL(10)='KK12' TABOBL(11)='KK13' TABOBL(12)='KK21' TABOBL(13)='KK22' TABOBL(14)='KK23' TABOBL(15)='KK31' TABOBL(16)='KK32' TABOBL(17)='KK33' * IF(IFOUR.EQ.-2) THEN JGFAC=33 TABFAC(33)='DIM3' ELSE JGFAC=32 ENDIF TABFAC(1)='RHOF' TABFAC(2)='ALPF' TABFAC(3)='ALPM' TABFAC(4)='PT11' TABFAC(5)='PH11' TABFAC(6)='PB11' TABFAC(7)='PT12' TABFAC(8)='PH12' TABFAC(9)='PB12' TABFAC(10)='PT13' TABFAC(11)='PH13' TABFAC(12)='PB13' TABFAC(13)='PT21' TABFAC(14)='PH21' TABFAC(15)='PB21' TABFAC(16)='PT22' TABFAC(17)='PH22' TABFAC(18)='PB22' TABFAC(19)='PT23' TABFAC(20)='PH23' TABFAC(21)='PB23' TABFAC(22)='PT31' TABFAC(23)='PH31' TABFAC(24)='PB31' TABFAC(25)='PT32' TABFAC(26)='PH32' TABFAC(27)='PB32' TABFAC(28)='PT33' TABFAC(29)='PH33' TABFAC(30)='PB33' TABFAC(31)='RHO ' TABFAC(32)='ALPH' ENDIF * IF (INMAT.EQ.0) THEN ** IRET = 0 GOTO 9999 ENDIF * 40 CONTINUE IF (NMAT.GE.INMAT) THEN IF (IPLAC.EQ.1) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * ISOTROPE IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF IF (IRET.EQ.0) GOTO 9999 ELSEIF (IPLAC.EQ.2) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * NORTON IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF IF (IRET.EQ.0) GOTO 9999 ELSEIF (IPLAC.EQ.3) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * ONERA IPLAC=2 * ELSE * INMAT=INMAT+1 ENDIF IF (IRET.EQ.0) GOTO 9999 ELSEIF (IPLAC.EQ.4) THEN INMAT=INMAT+1 IF (IPLAC.EQ.0) THEN * MAZARS IPLAC=1 * ELSE * INMAT=INMAT+1 ENDIF IF (IRET.EQ.0) GOTO 9999 ELSE * GOTO 9999 ENDIF * * En cas de creation de materiaux combinant plusieurs materiaux * deja existant Ex ELASTIQUE ISOTROPE PLASTIQUE PARFAIT FLUAGE N * GOTO 40 * * GOTO 9999 ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION CONTACT *-------------------------------------------------------------------- IF (ICONT.NE.0) THEN iplla =0 do iou = 1,NMAT iplla = max(iplac,iplla) enddo IF (iplla.EQ.0) THEN IRET = 0 GOTO 9999 ENDIF IF ((IPLLA.EQ.1) .OR. (IPLLA.EQ.2) .OR. (IPLLA.EQ.3)) THEN * CONTACT SIMPLE * -------------- INMAT=1 JGOBL=0 JGFAC=2 TABFAC(1)='JEU' TABFAC(2)='ADHE' ELSEIF (IPLLA .EQ. 4) THEN * FROTTEMENT DE COULOMB * --------------------- INMAT=1 JGOBL=1 TABOBL(1)='MU ' JGFAC=3 TABFAC(1)='COHE' TABFAC(2)='ADHE' TABFAC(3)='JEU' ELSEIF (IPLLA .EQ. 5) THEN * FROTTEMENTS DE CABLES * --------------------- INMAT=1 JGOBL=2 TABOBL(1)='FF ' TABOBL(2)='PHIF' JGFAC=0 ELSE IRET=0 RETURN ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION CONTRAINTE *-------------------------------------------------------------------- IF (ICNTR.NE.0) THEN CALL MODCLI(MOMODL,NMOD) iplla =0 do iou = 1,NMAT iplla = max(iplac,iplla) enddo IF (iplla.EQ.0) THEN IRET = 0 GOTO 9999 ENDIF IF (IPLLA .EQ. 1) THEN * ROTATION * --------------------- INMAT=1 JGOBL=1 TABOBL(1)='ANGL' ELSEIF (IPLLA .EQ. 2) THEN * DEPLACEMENT * --------------------- INMAT=1 JGOBL=1 TABOBL(1)='AMPL' ENDIF ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION MAGNETODYNAMIQUE *-------------------------------------------------------------------- IF (ICONV.NE.0) THEN * * FORMULATION EN COQUES IF(MFR.EQ.3) THEN IF(ISOT.NE.0) THEN JGOBL=3 TABOBL(1)='ETA' TABOBL(2)='PERM' TABOBL(3)='EPAI' ELSE IF(IORTH.NE.0) THEN JGOBL=4 TABOBL(1)='ETA1' TABOBL(2)='ETA2' TABOBL(3)='PERM' TABOBL(4)='EPAI' ENDIF ENDIF ENDIF * GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION FISSURE *-------------------------------------------------------------------- IF (ICONV.NE.0) THEN * * si POISEU_BLASIUS ou POISEU_COLEBROOK ou par defaut JGOBL=1 TABOBL(1)='RUGO' * si FROTTEMENT1 ou FROTTEMENT2 IF(IFT1.NE.0.OR.IFT2.NE.0) THEN JGOBL=7 TABOBL(2)='REC' TABOBL(3)='FK' TABOBL(4)='FA' TABOBL(5)='FB' TABOBL(6)='FC' TABOBL(7)='FD' ELSE * si FROTTEMENT3 ou FROTTEMENT4 IF(IFT1.NE.0.OR.IFT2.NE.0) THEN JGOBL=2 TABOBL(2)='FK' ENDIF ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION NAVIER_STOKES *-------------------------------------------------------------------- IF (INAST.GT.0) THEN iplac = imatee if (iplac.eq.0) then iret = 0 goto 9999 endif * JGOBL = 0 JGFAC = 0 if (iplac.eq.4) then *NLIN JGOBL0 = JGOBL JGOBL = JGOBL0+2 TABOBL(JGOBL0 + 1) ='RHO ' TABOBL(JGOBL0 + 2) ='MU ' JGOBL0 = JGFAC JGFAC = JGOBL0+1 TABFAC(JGOBL0 + 1) ='FREQ' endif GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION MELANGE *-------------------------------------------------------------------- IF (ICONV.NE.0) THEN if (iplac.eq.0) then goto 9999 endif INMAT=1 * * a priori elements massifs ou coques : pas de verif * IF (IPLAC.NE.3.AND.IPLAC.NE.4) THEN JGOBL = 0 JGFAC = 0 IF (IRET.EQ.0) GOTO 9999 ELSE ENDIF GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION LIAISON *-------------------------------------------------------------------- IF (ICONV.NE.0) THEN iplac = imatee if (iplac.eq.0) then iret = 0 goto 9999 endif * JGOBL = 0 JGFAC = 0 * 'SORT' facultatif dans tous les cas IF (iplac.EQ.1) THEN *PO_PL_FL JGOBL0 = JGOBL JGOBL = JGOBL0+7 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='INER' TABOBL(JGOBL0 + 3) ='CONV' TABOBL(JGOBL0 + 4) ='VISC' TABOBL(JGOBL0 + 5) ='PELO' TABOBL(JGOBL0 + 6) ='PRAP' TABOBL(JGOBL0 + 7) ='JFLU' ELSEIF (iplac.eq.2) then *PO_PL_FR JGOBL0 = JGOBL JGOBL = JGOBL0+7 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='GLIS' TABOBL(JGOBL0 + 5) ='ADHE' TABOBL(JGOBL0 + 6) ='RTAN' TABOBL(JGOBL0 + 7) ='ATAN' JGOBL0 = JGFAC JGFAC = JGOBL0 +2 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'LOIC' ELSEIF (iplac.eq.3) then *PO_PL JGOBL0 = JGOBL JGOBL= JGOBL0+3 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='JEU' JGOBL0 = JGFAC JGFAC = JGOBL0+4 TABFAC(JGOBL0+1) = 'LOIC' TABFAC(JGOBL0+2) = 'PERM' TABFAC(JGOBL0 + 3) ='SPLA' TABFAC(JGOBL0 + 4) ='AMOR' ELSEIF (iplac.eq.4) then *PO_PO_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='POIB' TABOBL(JGOBL0 + 5) ='ADHE' TABOBL(JGOBL0 + 6) ='RTAN' TABOBL(JGOBL0 + 7) ='ATAN' TABOBL(JGOBL0 + 8) ='GLIS' JGOBL0 = JGFAC JGFAC = JGOBL0 + 3 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'LOIC' TABFAC(JGOBL0+3) = 'MODE' ELSEIF (iplac.eq.5) then *PO_PO_DP JGOBL0 = JGOBL JGOBL= JGOBL0 + 6 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='ECRO' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='POIB' TABOBL(JGOBL0 + 5) ='PERM' TABOBL(JGOBL0 + 6) ='LOIC' JGOBL0 = JGFAC JGFAC = JGOBL0 + 1 TABFAC(JGOBL0+1) = 'AMOR' ELSEIF (iplac.eq.6) then *PO_PO_RP JGOBL0 = JGOBL JGOBL= JGOBL0 + 6 TABOBL(JGOBL0 + 1) ='AXRO' TABOBL(JGOBL0 + 2) ='ECRO' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='POIB' TABOBL(JGOBL0 + 5) ='PERM' TABOBL(JGOBL0 + 6) ='LOIC' JGOBL0 = JGFAC JGFAC = JGOBL0 + 2 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'ELAS' ELSEIF (iplac.eq.7) then *PO_PO JGOBL0 = JGOBL JGOBL= JGOBL0 + 5 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='POIB' TABOBL(JGOBL0 + 5) ='PERM' JGOBL0 = JGFAC JGFAC = JGOBL0 + 2 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'LOIC' ELSEIF (iplac.eq.8) then *PO_CE_MO JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='PCER' TABOBL(JGOBL0 + 4) ='RAYO' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' JGOBL0 = JGFAC JGFAC= JGOBL0 + 2 TABFAC(JGOBL0+1) = 'CINT' TABFAC(JGOBL0+2) = 'AMOR' ELSEIF (iplac.eq.9) then *PO_CE_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='EXCE' TABOBL(JGOBL0 + 4) ='RAYO' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' JGOBL0 = JGFAC JGFAC = JGOBL0 + 2 TABFAC(JGOBL0+1) = 'CINT' TABFAC(JGOBL0+2) = 'AMOR' ELSEIF (iplac.eq.10) then *PO_CE JGOBL0 = JGOBL JGOBL= JGOBL0 + 4 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='EXCE' TABOBL(JGOBL0 + 4) ='RAYO' JGOBL0 = JGFAC JGFAC = JGOBL0 + 1 TABFAC(JGOBL0+1) = 'AMOR' ELSEIF (iplac.eq.11) then *CE_PL_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='JEU' TABOBL(JGOBL0 + 4) ='RAYS' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' JGOBL0 = JGFAC JGFAC = JGOBL0+1 TABFAC(JGOBL0+1) = 'AMOR' ELSEIF (iplac.eq.12) then *CE_CE_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 9 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='EXCE' TABOBL(JGOBL0 + 4) ='RAYS' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' TABOBL(JGOBL0 + 9) ='RAYB' JGOBL0 = JGFAC JGFAC = JGOBL0 + 2 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'CINT' ELSEIF (iplac.eq.13.or.iplac.eq.14) then *PR_PR_IN ou PR_PR_EX JGOBL0 = JGOBL JGOBL= JGOBL0 + 5 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='RAID' TABOBL(JGOBL0 + 3) ='PFIX' TABOBL(JGOBL0 + 4) ='PMOB' TABOBL(JGOBL0 + 5) ='ERAI' ELSEIF (iplac.eq.15) then *LI_LI_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='LIMA' TABOBL(JGOBL0 + 3) ='LIES' TABOBL(JGOBL0 + 4) ='RAID' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' JGOBL0 = JGFAC JGFAC = JGOBL0 + 4 TABFAC(JGOBL0+1) = 'JEU' TABFAC(JGOBL0+2) = 'AMOR' TABFAC(JGOBL0+3) = 'RECH' TABFAC(JGOBL0+4) = 'SYME' ELSEIF (iplac.eq.16) then *LI_CE_FR JGOBL0 = JGOBL JGOBL= JGOBL0 + 8 TABOBL(JGOBL0 + 1) ='NORM' TABOBL(JGOBL0 + 2) ='LIMA' TABOBL(JGOBL0 + 3) ='LIES' TABOBL(JGOBL0 + 4) ='RAID' TABOBL(JGOBL0 + 5) ='GLIS' TABOBL(JGOBL0 + 6) ='ADHE' TABOBL(JGOBL0 + 7) ='RTAN' TABOBL(JGOBL0 + 8) ='ATAN' JGOBL0 = JGFAC JGFAC = JGOBL0 + 5 TABFAC(JGOBL0+1) = 'AMOR' TABFAC(JGOBL0+2) = 'RECH' TABFAC(JGOBL0+3) = 'RAYO' TABFAC(JGOBL0+4) = 'ACTN' TABFAC(JGOBL0+5) = 'INVE' ELSEIF (iplac.eq.17) then *PA_FL_RO JGOBL0 = JGOBL JGOBL= JGOBL0+10 TABOBL(JGOBL0 + 1) ='LONG' TABOBL(JGOBL0 + 2) ='RAYO' TABOBL(JGOBL0 + 3) ='VISC' TABOBL(JGOBL0 + 4) ='RHOF' TABOBL(JGOBL0 + 5) ='PADM' TABOBL(JGOBL0 + 6) ='VROT' TABOBL(JGOBL0 + 7) ='EPSI' TABOBL(JGOBL0 + 8) ='PHII' TABOBL(JGOBL0 + 9) ='AFFI' TABOBL(JGOBL0 + 10)='TLOB' JGOBL0 = JGFAC JGFAC = JGOBL0 + 1 TABFAC(JGOBL0+1) = 'AMOR' ELSEIF (iplac.eq.23) then *NEWMARK MODAL JGOBL0 = JGOBL JGOBL= JGOBL0+3 TABOBL(JGOBL0 + 1) ='JEU' TABOBL(JGOBL0 + 2) ='MASS' TABOBL(JGOBL0 + 3) ='FREQ' JGOBL0 = JGFAC JGFAC = JGOBL0 + 3 TABFAC(JGOBL0+1) = 'EXCE' TABFAC(JGOBL0+2) = 'FROT' TABFAC(JGOBL0+3) = 'MOFR' ENDIF * 'SORT' facultatif dans tous les cas JGOBL0 = JGFAC JGFAC= JGOBL0+1 TABFAC(JGOBL0+1) = 'SORT' GOTO 9999 ENDIF *-------------------------------------------------------------------- * CAS DE LA FORMULATION ELECTROSTATIQUE *-------------------------------------------------------------------- IF (IELEC.NE.0) THEN C -- Permittivite isotrope IF (IMATEE.EQ.1) THEN C* IF (CMATEE.EQ.'ISOTROPE') THEN C* IF (MATMOD(1).EQ.'ISOTROPE ') THEN JGOBL = 1 TABOBL(1)='PEL ' C -- Permittivite orthotrope ELSEIF (IMATEE.EQ.2) THEN C* ELSEIF (CMATEE.EQ.'ORTHOTRO') THEN C* ELSEIF (MATMOD(1).EQ.'ORTHOTROPE ') THEN C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IDIM.EQ.2) THEN IF (IFOMOD.NE.1) THEN JGOBL = 4 TABOBL(1) = 'PE1 ' TABOBL(2) = 'PE2 ' TABOBL(3) = 'V1X ' TABOBL(4) = 'V1Y ' C ---- Elements massifs bidimensionnels FOURIER ELSE JGOBL = 5 TABOBL(1) = 'PE1 ' TABOBL(2) = 'PE2 ' TABOBL(3) = 'PE3 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' ENDIF C ---- Elements massifs TRIDimensionnels ELSEIF (IDIM.EQ.3) THEN JGOBL = 9 TABOBL(1) = 'PE1 ' TABOBL(2) = 'PE2 ' TABOBL(3) = 'PE3 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' TABOBL(6) = 'V1Z ' TABOBL(7) = 'V2X ' TABOBL(8) = 'V2Y ' TABOBL(9) = 'V2Z ' ENDIF C -- Permittivite anisotrope ELSEIF (IMATEE.EQ.3) THEN C* ELSEIF (CMATEE.EQ.'ANISOTRO') THEN C* ELSEIF (MATMOD(1).EQ.'ANISOTROPE ') THEN C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IDIM.EQ.2) THEN IF (IFOMOD.NE.1) THEN JGOBL = 5 TABOBL(1) = 'PE11 ' TABOBL(2) = 'PE22 ' TABOBL(3) = 'PE21 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' C ---- Elements massifs bidimensionnels FOURIER ELSE JGOBL = 6 TABOBL(1) = 'PE11 ' TABOBL(2) = 'PE22 ' TABOBL(3) = 'PE21 ' TABOBL(4) = 'PE33 ' TABOBL(5) = 'V1X ' TABOBL(6) = 'V1Y ' ENDIF C ---- Elements massifs TRIDimensionnels ELSEIF (IDIM.EQ.3) THEN JGOBL = 12 TABOBL( 1) = 'PE11 ' TABOBL( 2) = 'PE22 ' TABOBL( 3) = 'PE33 ' TABOBL( 4) = 'PE21 ' TABOBL( 5) = 'PE31 ' TABOBL( 6) = 'PE32 ' TABOBL( 7) = 'V1X ' TABOBL( 8) = 'V1Y ' TABOBL( 9) = 'V1Z ' TABOBL(10) = 'V2X ' TABOBL(11) = 'V2Y ' TABOBL(12) = 'V2Z ' ENDIF ELSE IRET = 0 ENDIF GOTO 9999 ENDIF *----------------------------------------------------------------------- * CAS DE LA FORMULATION 'DIFFUSION' *----------------------------------------------------------------------- IF (IDIFF.NE.0) THEN JGFAC = 1 TABFAC(1)='CDIF' C- Diffusion isotrope IF (IMATEE .EQ. 1) THEN JGOBL = 1 C*8 TABOBL(1) = 'KDIFF ' TABOBL(1) = 'KD ' C- Diffusion orthotrope ELSEIF (IMATEE .EQ. 2) THEN C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IDIM.EQ.2) THEN IF (IFOMOD.NE.1) THEN JGOBL = 4 TABOBL(1) = 'KD1 ' TABOBL(2) = 'KD2 ' TABOBL(3) = 'V1X ' TABOBL(4) = 'V1Y ' C --- Elements massifs bidimensionnels FOURIER ELSE JGOBL = 5 TABOBL(1) = 'KD1 ' TABOBL(2) = 'KD2 ' TABOBL(3) = 'KD3 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' ENDIF C --- Elements massifs TRIDimensionnels ELSEIF (IDIM.EQ.3) THEN JGOBL = 9 TABOBL(1) = 'KD1 ' TABOBL(2) = 'KD2 ' TABOBL(3) = 'KD3 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' TABOBL(6) = 'V1Z ' TABOBL(7) = 'V2X ' TABOBL(8) = 'V2Y ' TABOBL(9) = 'V2Z ' ENDIF C - Diffusion anisotrope ELSEIF (IMATEE.EQ.3) THEN C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IDIM.EQ.2) THEN IF (IFOMOD.NE.1) THEN JGOBL = 5 TABOBL(1) = 'KD11 ' TABOBL(2) = 'KD22 ' TABOBL(3) = 'KD21 ' TABOBL(4) = 'V1X ' TABOBL(5) = 'V1Y ' C --- Elements massifs bidimensionnels FOURIER ELSE JGOBL = 6 TABOBL(1) = 'KD11 ' TABOBL(2) = 'KD22 ' TABOBL(3) = 'KD21 ' TABOBL(4) = 'KD33 ' TABOBL(5) = 'V1X ' TABOBL(6) = 'V1Y ' ENDIF C --- Elements massifs TRIDimensionnels ELSEIF (IDIM.EQ.3) THEN JGOBL = 12 TABOBL( 1) = 'KD11 ' TABOBL( 2) = 'KD22 ' TABOBL( 3) = 'KD33 ' TABOBL( 4) = 'KD21 ' TABOBL( 5) = 'KD31 ' TABOBL( 6) = 'KD32 ' TABOBL( 7) = 'V1X ' TABOBL( 8) = 'V1Y ' TABOBL( 9) = 'V1Z ' TABOBL(10) = 'V2X ' TABOBL(11) = 'V2Y ' TABOBL(12) = 'V2Z ' ENDIF ENDIF JGOBL0 = JGOBL JGFAC0 = JGFAC * Loi FICK : * ------------ IF (INATUU.EQ.0) THEN * Rien a ajouter * Loi SORET : * ------------- ELSEIF (INATUU.EQ.1) THEN C -- Diffusion isotrope IF (IMATEE.EQ.1) THEN JGOBL = JGOBL+1 C*8 TABOBL(JGOBL0+1) = 'KSORET ' TABOBL(JGOBL0+1) = 'KS ' C -- Diffusion orthotrope ELSEIF (IMATEE.EQ.2) THEN C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IFOMOD.LE.0) THEN JGOBL = JGOBL + 2 C*8 TABOBL(JGOBL0+1) = 'KSORET1 ' C*8 TABOBL(JGOBL0+2) = 'KSORET2 ' TABOBL(JGOBL0+1) = 'KS1 ' TABOBL(JGOBL0+2) = 'KS2 ' C ---- Elements massifs bidimensionnels FOURIER C ---- Elements massifs TRIDimensionnels ELSE JGOBL = JGOBL + 3 C*8 TABOBL(JGOBL0+1) = 'KSORET1 ' C*8 TABOBL(JGOBL0+2) = 'KSORET2 ' C*8 TABOBL(JGOBL0+3) = 'KSORET3 ' TABOBL(JGOBL0+1) = 'KS1 ' TABOBL(JGOBL0+2) = 'KS2 ' TABOBL(JGOBL0+3) = 'KS3 ' ENDIF C -- Diffusion anisotrope ELSEIF (IMATEE.EQ.3) THEN C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE IF (IFOMOD.LE.0) THEN JGOBL = JGOBL + 3 C*8 TABOBL(JGOBL0+1) = 'KSORET11' C*8 TABOBL(JGOBL0+2) = 'KSORET22' C*8 TABOBL(JGOBL0+3) = 'KSORET21' TABOBL(JGOBL0+1) = 'KS11 ' TABOBL(JGOBL0+2) = 'KS22 ' TABOBL(JGOBL0+3) = 'KS21 ' C ---- Elements massifs bidimensionnels FOURIER ELSEIF (IFOMOD.EQ.1) THEN JGOBL = JGOBL + 4 C*8 TABOBL(JGOBL0+1) = 'KSORET11' C*8 TABOBL(JGOBL0+2) = 'KSORET22' C*8 TABOBL(JGOBL0+3) = 'KSORET21' C*8 TABOBL(JGOBL0+4) = 'KSORET33' TABOBL(JGOBL0+1) = 'KS11 ' TABOBL(JGOBL0+2) = 'KS22 ' TABOBL(JGOBL0+3) = 'KS21 ' TABOBL(JGOBL0+4) = 'KS33 ' C ---- Elements massifs TRIDimensionnels ELSEIF (IFOMOD.EQ.2) THEN JGOBL = JGOBL + 6 C*8 TABOBL(JGOBL0+1) = 'KSORET11' C*8 TABOBL(JGOBL0+2) = 'KSORET22' C*8 TABOBL(JGOBL0+3) = 'KSORET33' C*8 TABOBL(JGOBL0+4) = 'KSORET21' C*8 TABOBL(JGOBL0+5) = 'KSORET31' C*8 TABOBL(JGOBL0+6) = 'KSORET32' TABOBL(JGOBL0+1) = 'KS11 ' TABOBL(JGOBL0+2) = 'KS22 ' TABOBL(JGOBL0+3) = 'KS33 ' TABOBL(JGOBL0+4) = 'KS21 ' TABOBL(JGOBL0+5) = 'KS31 ' TABOBL(JGOBL0+6) = 'KS32 ' ENDIF ELSE C -- Erreur ne devant pas survenir RETURN ENDIF * Loi UTILISATEUR : * ------------------- * ELSEIF (INATUU.EQ.2) THEN * Par defaut pas de composantes autres que celles de l'utilisateur ELSEIF (INATUU.EQ.5) THEN * diffusion avec dissipation affine JGOBL = JGOBL+2 TABOBL(JGOBL0+1)='DAF0' TABOBL(JG