1 : C ERRARE SOURCE FANDEUR 12/04/11 21:15:18 7350 2 : SUBROUTINE ERRARE(IPMODL,IPCHE1,IPCHE2,XERR,IPCHRR) 3 : *_______________________________________________________________________ 4 : * 5 : * CALCUL D'ERREUR SUR LES CONTRAINTES 6 : * 7 : * METHODE ZIENKIEWICZ - VERSION BARZIC ET RICHARD 8 : * 9 : *_______________________________________________________________________ 10 : * 11 : * ENTREES : 12 : * --------- 13 : * IPMODL POINTEUR SUR UN MMODEL 14 : * IPCHE1 POINTEUR SUR UN MCHAML DE CONTRAINTES 15 : * IPCHE2 POINTEUR SUR UN MCHAML DE MATERIAU 16 : * 17 : * 18 : * SORTIES : 19 : * --------- 20 : * XERR ERREUR GENERALE SCALAIRE 21 : * IPCHRR MCHAML D'ERREUR PAR ELEMENT 22 : * = 0 SI ERREUR 23 : * 24 : * PASSAGE AUX NOUVEAUX CHAMELEM PAR JM CAMPENON LE 02/91 25 : * 26 : *------------------------------------------------------------------ 27 : * 28 : IMPLICIT INTEGER(I-N) 29 : IMPLICIT REAL*8 (X) 30 : * 31 : -INC CCOPTIO 32 : * 33 : IPCHRR=0 33 : CALL SIGMOY(IPMODL,IPCHE1,IPOSGM) 34 : IF (IPOSGM.EQ.0) GOTO 310 35 : * 35 : CALL ADCHEL(IPCHE1,IPOSGM,IPODSI,-1) 36 : IF (IPODSI.EQ.0) GOTO 300 37 : * 38 : * CHAMP DE VARIABLES INTERNES : IPCH3 = 0 ici (a lire dans ERRE ?) 39 : * 40 : IPCHE3 = 0 41 : * 41 : CALL ELAS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPOEPS,IRET) 42 : IF (IRET.NE.1) GOTO 290 43 : * 43 : CALL ELAS1(IPMODL,IPODSI,IPCHE2,IPCHE3,IPODEP,IRET) 44 : IF (IRET.NE.1) GOTO 280 45 : * 45 : CALL ENERCA(IPMODL,IPCHE1,IPOEPS,IPOENE) 46 : IF (IPOENE.EQ.0) GOTO 270 47 : * 47 : CALL ENERCA(IPMODL,IPODSI,IPODEP,IPODEN) 48 : IF (IPODEN.EQ.0) GOTO 260 49 : * 49 : CALL INTGCA(IPMODL,IPOENE,0,1,IRET,XRET,IPORET) 50 : IF (IRET.EQ.0) GOTO 250 51 : * 51 : CALL INTGCA(IPMODL,IPODEN,0,1,IRET,XDRE,IPODRE) 52 : IF (IRET.EQ.0) GOTO 240 53 : * 54 : XERR=(XDRE/(XRET+XDRE))**(.5) 54 : CALL ADCHEL(IPORET,IPODRE,IPOSOM,1) 54 : CALL PUCHEL(IPOSOM,-1,IPOFAC,IRET) 55 : IF(IRET.EQ.0) GOTO 220 56 : * 57 : IPLMOT = -1 58 : IPLREE = -1 58 : CALL MUCHSC(IPMODL,IPOFAC,IPODRE,IPOCAR, 59 : & IPLMOT,IPLMOT,IPLMOT,IPLREE,1) 60 : IF (IPOCAR.EQ.0) GOTO 210 61 : * 61 : CALL PUCHEL(IPOCAR,1,IPCHRR,IRET) 62 : IF (IRET.NE.0) GOTO 200 63 : * 64 : IPCHRR=0 64 : 200 CALL DTCHAM(IPOCAR) 64 : 210 CALL DTCHAM(IPOFAC) 64 : 220 CALL DTCHAM(IPOSOM) 64 : 230 CALL DTCHAM(IPODRE) 64 : 240 CALL DTCHAM(IPORET) 64 : 250 CALL DTCHAM(IPODEN) 64 : 260 CALL DTCHAM(IPOENE) 64 : 270 CALL DTCHAM(IPODEP) 64 : 280 CALL DTCHAM(IPOEPS) 64 : 290 CALL DTCHAM(IPODSI) 64 : 300 CALL DTCHAM(IPOSGM) 65 : 310 RETURN 66 : END 67 : 68 :
© Cast3M 2003 - All rights reserved.
Disclaimer