1 : C TRESCA SOURCE FANDEUR 10/12/14 21:20:01 6812 2 : SUBROUTINE TRESCA 3 : *--------------------------------------------------------------------- 4 : * 5 : * CALCUL DE LA CONTRAINTE DE TRESCA 6 : * 7 : * SYNTAXE : 8 : * --------- 9 : * CHAL1 = TRESCA MODL1 SIG1 (CAR1) (MOT1) 10 : * 11 : * MODL1 = Modele de calcul (type MMODEL) 12 : * SIG1 = Objet de type MCHALM de CONTRAINTES 13 : * CAR1 = Objet de type MCHAML de CARACTERISTIQUE 14 : * MOT1 = Mot qui indique pour les coques ou l'on calcule les 15 : * contraintes 16 : * CHEL1 = Objet de type MCHAML de sous type SCALAIRE (aux pts 17 : * de gauss pour les contraintes) 18 : * 19 : * Passage aux nouveaux chamelem par jm CAMPENON le 04/91 20 : * 21 : *-------------------------------------------------------------------- 22 : * 23 : IMPLICIT INTEGER(I-N) 24 : IMPLICIT REAL*8(A-H,O-Z) 25 : * 26 : -INC CCOPTIO 27 : * 28 : CHARACTER*4 MOMOT 29 : * 30 : MOMOT=' ' 31 : IMIL=0 32 : * 33 : IPMODL=0 34 : IPCHE1=0 35 : IPCHE2=0 36 : * 36 : CALL LIROBJ ('MMODEL',IPMODL,1,IRET) 37 : IF(IERR.NE.0) RETURN 38 : * 38 : CALL LIROBJ('MCHAML',IPCHE1,1,IRET) 39 : IF(IERR.NE.0) RETURN 40 : * 41 : * On lit les CARACTERISTIQUES (facultatif) 42 : * 42 : CALL LIROBJ('MCHAML',IPCHE2,0,IRET) 43 : IF(IERR.NE.0) RETURN 44 : * 45 : * Lecture d'un mot si necessaire 46 : * 46 : CALL LIRCHA(MOMOT,0,IRET) 47 : IF(IERR.NE.0) RETURN 48 : IF (IRET.NE.0) THEN 49 : IF (MOMOT.EQ.'SUPE') IMIL=1 50 : IF (MOMOT.EQ.'INFE') IMIL=-1 51 : IF (MOMOT.NE.'SUPE'.AND.MOMOT.NE.'INFE'.AND.MOMOT.NE.'MOYE') THEN 52 : MOTERR(1:4)=MOMOT 52 : CALL ERREUR(7) 53 : RETURN 54 : ENDIF 55 : ENDIF 56 : * 57 : IRET=0 57 : CALL TRESK(IPMODL,IPCHE1,IPCHE2,IMIL,IPSCAL,IRET) 58 : * 58 : IF(IRET.NE.0) CALL ECROBJ('MCHAML',IPSCAL) 59 : * 60 : RETURN 61 : END 62 : 63 :
© Cast3M 2003 - All rights reserved.
Disclaimer