1 : C CEREPS SOURCE CHAT 05/01/12 21:52:20 5004 2 : SUBROUTINE CEREPS(COORD,DEPLA,DEFO,IDEUX,KERRE) 3 : C----------------------------------------------------------------------- 4 : C CALCUL DES DEFORMATIONS POUR LA CERCE 5 : C 6 : C ENTREE 7 : C DEPLA DEPLACEMENTS 8 : C COORD COORDONNEES DU NOEUD 9 : C IDEUX INDICE POUR DEFORMATIONS DU 2-EME ORDRE 10 : C 11 : C SORTIE 12 : C DEFO DEFORMATION AU NOEUD 13 : C----------------------------------------------------------------------- 14 : IMPLICIT INTEGER(I-N) 15 : IMPLICIT REAL*8(A-H,O-Z) 16 : -INC CCOPTIO 17 : -INC CCREEL 18 : DIMENSION DEFO(*),COORD(3,*),DEPLA(*) 19 : * 20 : KERRE=0 21 : R = COORD(1,1) 22 : IF(ABS(R).LE.XPETIT) THEN 23 : KERRE=1 24 : RETURN 25 : ENDIF 26 : * 27 : IF(IFOUR.EQ.0) THEN 28 : EPS=DEPLA(1)/R 29 : IF(IDEUX.NE.0) THEN 30 : EPS = EPS + 0.5D0*(DEPLA(1)/R)**2 31 : ENDIF 32 : ELSE IF(IFOUR.EQ.1) THEN 33 : AN=DBLE(NIFOUR) 34 : EPS=(DEPLA(1) + AN * DEPLA (3))/R 35 : * IF(IDEUX.NE.0) THEN 36 : * CECI N'A PAS GRAND SENS CAR IL Y A DES TERMES 37 : * EN SIN**2 ET EN COS**2 38 : * ENDIF 39 : ELSE 40 : KERRE=2 41 : RETURN 42 : ENDIF 43 : DEFO(1)=EPS 44 : RETURN 45 : END 46 : 47 : 48 :
© Cast3M 2003 - All rights reserved.
Disclaimer