Télécharger dfouri.eso

Retour à la liste

Numérotation des lignes :

  1. C DFOURI SOURCE FANDEUR 10/12/14 21:15:52 6812
  2. SUBROUTINE DFOURI
  3. C=======================================================================
  4. C
  5. C OPERATEUR DFOURI
  6. C
  7. C SYNTAXE I CHPT I = DFOURI I CHPT I ANGL
  8. C I CHEL I I CHEL I
  9. C
  10. C CHPT : CHPOINT DE TYPE FORCES OU DEPLACEMENTS
  11. C CHEL : CHAMP/ELEMENT DE CONTRAINTES OU DEFORMATIONS (MCHAML)
  12. C ANGL : VALEUR DE L ANGLE
  13. C
  14. C DANS LE CAS D UNE ANALYSE EN SERIE DE FOURIER , CALCULE LES VALEURS
  15. C DU CHPOINT OU DU MCHALM POUR L ANGLE DONNE
  16. C
  17. C J BROCHARD MARS 87
  18. C PASSAGE AUX NOUVEAUX MCHAMLS P DOWLATYARI SEPT. 91
  19. C=======================================================================
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. -INC CCOPTIO
  23. -INC SMCHPOI
  24. C
  25. C ******** ON LIT L ANGLE
  26. C
  27. CALL LIRREE(ANGL,1,IRETOU)
  28. IF (IERR.NE.0) RETURN
  29. C
  30. CALL LIROBJ('CHPOINT',IPCHP1,0,IRT1)
  31. IF(IERR.NE.0)RETURN
  32.  
  33. IF(IRT1.NE.0) THEN
  34. C
  35. C C'EST BIEN UN CHPOINT QU'IL FAUT LIRE
  36. C
  37. CALL DFOUR1(IPCHP1,ANGL,IPCHP2)
  38. IF (IERR.NE.0) RETURN
  39. * ATTRIBUTION D'UNE NATURE IDENTIQUE AU CHPO ARGUMENT
  40. MCHPOI = IPCHP1
  41. MCHPO1 = IPCHP2
  42. SEGACT ,MCHPOI,MCHPO1
  43. NAT = MAX ( 1 , MCHPO1.JATTRI(/1) )
  44. NSOUPO=MCHPO1.IPCHP(/1)
  45. SEGADJ MCHPO1
  46. IF ( JATTRI(/1) .GE. 1 ) THEN
  47. MCHPO1.JATTRI(1) = JATTRI(1)
  48. ELSE
  49. JATTRI(1) = 0
  50. ENDIF
  51. SEGDES ,MCHPOI,MCHPO1
  52. CALL ECROBJ('CHPOINT',IPCHP2)
  53. ELSE
  54. C
  55. C LE CHAMP EST UN CHAMP/ELEMENT
  56. C
  57. CALL LIROBJ('MCHAML',IPCHE1,1,IRT3)
  58. IF(IERR.NE.0)RETURN
  59. CALL DFOUR2(IPCHE1,ANGL,IPCHE2)
  60. IF (IERR.EQ.0) CALL ECROBJ('MCHAML',IPCHE2)
  61. ENDIF
  62. RETURN
  63. END
  64.  
  65.  
  66.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales