1 : C MUTU SOURCE CHAT 09/10/09 21:21:26 6519 2 : SUBROUTINE MUTU 3 : * 4 : *_______________________________________________________________________ 5 : * 6 : * OPERATEUR MUTU 7 : * 8 : * ELEMENTS FINIS : 9 : * ============== 10 : * 11 : * SYNTAXE : 12 : * ----------- 13 : * 14 : * RIG1=MUTU MOD1 CAR1 GEO2 15 : * 16 : * MOD1 MODELE DE CALCUL, TYPE MMODEL 17 : * CAR1 CHAMP PAR ELEMENT DE CARACTERISTIQUES 18 : * GEOMETRIQUES ET MATERIELLES, TYPE MCHAML 19 : * GEO2 MAILLAGE INDUCTEUR SUPPORT DU COURANT, TYPE MELEME 20 : * 21 : *----------------------------------------------------------------------- 22 : * Y. STEPHAN (FEVRIER 97) 23 : *_______________________________________________________________________ 24 : IMPLICIT INTEGER(I-N) 25 : IMPLICIT REAL*8(A-H,O-Z) 26 : * 27 : -INC CCOPTIO 28 : -INC SMRIGID 29 : C 30 : C LECTURE DU MODELE 31 : C 32 : IPMODL=0 33 : IPCHEL=0 34 : IPGEOM=0 35 : IPMUTU=0 36 : C 37 : MOTERR(1:8)=' MODELE ' 37 : CALL MESLIR(-137) 37 : CALL LIROBJ('MMODEL',IPMODL,1,IRT1) 38 : IF(IERR.NE.0) RETURN 39 : C_______________________________________________________________________ 40 : C LECTURE DU 1 IER MCHAML 41 : C 42 : MOTERR(1:8)=' MCHAML ' 42 : CALL MESLIR(-137) 42 : CALL LIROBJ('MCHAML',IPCHEL,1,IRT2) 43 : IF(IERR.NE.0) RETURN 44 : C 45 : C_______________________________________________________________________ 46 : C LECTURE DU MAILLAGE 47 : C 48 : MOTERR(1:8)=' MAILLAGE ' 48 : CALL MESLIR(-137) 49 : IPGEOM=0 49 : CALL LIROBJ('MAILLAGE',IPGEOM,0,IRT3) 50 : IF(IERR.NE.0) RETURN 51 : C 52 : C CALCUL DE LA MATRICE DE MUTUELLE 52 : CALL MUTU1(IPMODL,IPCHEL,IPGEOM,IPMUTU) 53 : IF(IERR.NE.0) RETURN 54 : C 55 : C ECRITURE DU RESULTAT 55 : CALL ECROBJ('RIGIDITE',IPMUTU) 56 : C 57 : RETURN 58 : END 59 : 60 : 61 : 62 : 63 : 64 : 65 : 66 : 67 : 68 :
© Cast3M 2003 - All rights reserved.
Disclaimer