C DFOURI SOURCE CB215821 20/11/25 13:25:20 10792 SUBROUTINE DFOURI C======================================================================= C C OPERATEUR DFOURI C C SYNTAXE I CHPT I = DFOURI I CHPT I ANGL C I CHEL I I CHEL I C C CHPT : CHPOINT DE TYPE FORCES OU DEPLACEMENTS C CHEL : CHAMP/ELEMENT DE CONTRAINTES OU DEFORMATIONS (MCHAML) C ANGL : VALEUR DE L ANGLE C C DANS LE CAS D UNE ANALYSE EN SERIE DE FOURIER , CALCULE LES VALEURS C DU CHPOINT OU DU MCHALM POUR L ANGLE DONNE C C J BROCHARD MARS 87 C PASSAGE AUX NOUVEAUX MCHAMLS P DOWLATYARI SEPT. 91 C======================================================================= IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC PPARAM -INC CCOPTIO -INC SMCHPOI C C ******** ON LIT L ANGLE C CALL LIRREE(ANGL,1,IRETOU) IF (IERR.NE.0) RETURN C CALL LIROBJ('CHPOINT',IPCHP1,0,IRT1) IF(IERR.NE.0)RETURN IF(IRT1.NE.0) THEN C C C'EST BIEN UN CHPOINT QU'IL FAUT LIRE C CALL DFOUR1(IPCHP1,ANGL,IPCHP2) IF (IERR.NE.0) RETURN * ATTRIBUTION D'UNE NATURE IDENTIQUE AU CHPO ARGUMENT MCHPOI = IPCHP1 MCHPO1 = IPCHP2 SEGACT ,MCHPOI,MCHPO1 NAT = MAX ( 1 , MCHPO1.JATTRI(/1) ) NSOUPO=MCHPO1.IPCHP(/1) SEGADJ MCHPO1 IF ( JATTRI(/1) .GE. 1 ) THEN MCHPO1.JATTRI(1) = JATTRI(1) ELSE JATTRI(1) = 0 ENDIF SEGDES ,MCHPOI,MCHPO1 CALL ECROBJ('CHPOINT',IPCHP2) ELSE C C LE CHAMP EST UN CHAMP/ELEMENT C CALL LIROBJ('MCHAML',IPCHE1,1,IRT3) IF(IERR.NE.0)RETURN CALL DFOUR2(IPCHE1,ANGL,IPCHE2) IF (IERR.EQ.0) CALL ECROBJ('MCHAML',IPCHE2) ENDIF RETURN END