1 : C RECANG SOURCE CHAT 05/01/13 02:46:30 5004 2 : SUBROUTINE RECANG(ITRES,IPOS,IPTR,XANGLE) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : *--------------------------------------------------------------------* 6 : * * 7 : * Recherche le déplacements généralisé de rotation * 8 : * * 9 : * Param}tres: * 10 : * * 11 : * e ITRES table résultat issue de l'opérateur Dyne * 12 : * e IPOS position de XTEMP dans le listreel des temps * 13 : * e IPTR point de repère du mode de rotation,dans la base modale* 14 : * es XANGLE valeur de la rotation * 15 : * Auteur, date de création: * 16 : * * 17 : * Samuel DURAND le 15 Octobre 1996. * 18 : * * 19 : *--------------------------------------------------------------------* 20 : -INC CCOPTIO 21 : -INC SMCHPOI 22 : -INC SMELEME 23 : LOGICAL L0,L1 24 : * Recherche du champoint des déplacements généralisés 24 : CALL ACCTAB(ITRES,'ENTIER',IPOS,X0,' ',L0,IP0, 25 : & 'TABLE',I1,X1,' ',L1,ITDEP) 25 : CALL ACCTAB(ITDEP,'MOT',I0,X0,'DEPLACEMENT',L0,IP0, 26 : & 'CHPOINT',I1,X1,' ',L1,ICHPT) 27 : MCHPOI = ICHPT 28 : IF (MCHPOI.EQ.0) THEN 29 : * le CHPOINT des déplacements modaux est nul 30 : MOTERR(1:8) = 'RCDEPL' 30 : CALL ERREUR(170) 31 : RETURN 32 : ENDIF 33 : SEGACT MCHPOI 34 : NSOU = IPCHP(/1) 35 : DO 10 ISOU = 1,NSOU 36 : MSOUPO = IPCHP(ISOU) 37 : SEGACT MSOUPO 38 : * on cherche la valeur du champ correspondant au point de 39 : * rotation 40 : MELEME = IGEOC 41 : SEGACT MELEME 42 : MPOVAL = IPOVAL 43 : SEGACT MPOVAL 44 : N2 = NUM(/2) 45 : DO 12 I = 1,N2 46 : IF (NUM(1,I).EQ.IPTR) THEN 47 : XANGLE=VPOCHA(I,1) 48 : ENDIF 49 : 12 CONTINUE 50 : SEGDES MPOVAL,MELEME,MSOUPO 51 : 10 CONTINUE 52 : SEGDES MCHPOI 53 : * 54 : END 55 : 56 : 57 :
© Cast3M 2003 - All rights reserved.
Disclaimer