C MUTU SOURCE CB215821 18/09/21 21:16:16 9930 SUBROUTINE MUTU * *_______________________________________________________________________ * * OPERATEUR MUTU * * ELEMENTS FINIS : * ============== * * SYNTAXE : * ----------- * * RIG1=MUTU MOD1 CAR1 GEO2 * * MOD1 MODELE DE CALCUL, TYPE MMODEL * CAR1 CHAMP PAR ELEMENT DE CARACTERISTIQUES * GEOMETRIQUES ET MATERIELLES, TYPE MCHAML * GEO2 MAILLAGE INDUCTEUR SUPPORT DU COURANT, TYPE MELEME * *----------------------------------------------------------------------- * Y. STEPHAN (FEVRIER 97) *_______________________________________________________________________ IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) * -INC PPARAM -INC CCOPTIO -INC SMRIGID C C LECTURE DU MODELE C IPMODL=0 IPCHEL=0 IPGEOM=0 IPMUTU=0 C MOTERR(1:8)=' MODELE ' CALL MESLIR(-137) CALL LIROBJ('MMODEL',IPMODL,1,IRT1) IF(IERR.NE.0) RETURN CALL ACTOBJ('MMODEL ',IPMODL,1) C_______________________________________________________________________ C LECTURE DU 1 IER MCHAML C MOTERR(1:8)=' MCHAML ' CALL MESLIR(-137) CALL LIROBJ('MCHAML',IPIN,1,IRT2) IF(IERR.NE.0) RETURN CALL ACTOBJ('MCHAML ',IPIN,1) CALL REDUAF(IPIN,IPMODL,IPCHEL,0,IR,KER) IF(IR .NE. 1) CALL ERREUR(KER) IF(IERR .NE. 0) RETURN C C_______________________________________________________________________ C LECTURE DU MAILLAGE C MOTERR(1:8)=' MAILLAGE ' CALL MESLIR(-137) IPGEOM=0 CALL LIROBJ('MAILLAGE',IPGEOM,0,IRT3) IF(IERR.NE.0) RETURN C C CALCUL DE LA MATRICE DE MUTUELLE CALL MUTU1(IPMODL,IPCHEL,IPGEOM,IPMUTU) IF(IERR.NE.0) RETURN C C ECRITURE DU RESULTAT CALL ECROBJ('RIGIDITE',IPMUTU) C RETURN END