Télécharger la plaquette Cast3M


fr en

Télécharger prchl2.eso

Retour à la liste

Numérotation des lignes : 

prchl2

  1. C PRCHL2    SOURCE    PV090527  25/01/07    14:42:54     12115          
  2. C
  3. C-----------------------------------------------------------------------
  4. C Sous-programme de l'operateur PRESSION :
  5. C
  6. C    IPCHE2 = PRES IPMOD1 MOT1 VAL1 (MAIL1) ;
  7. C
  8. C   En entree :
  9. C   -----------
  10. C     IPMOD1 : MMODEL de CHARGEMENT PRESSION, deja reduit sur MAIL1
  11. C     MOT1   : MOT, nom de la composante de pression
  12. C     XP     : FLOTTANT, valeur de la pression
  13. C
  14. C   En sortie :
  15. C   -----------
  16. C     IPCHE1 : contient le MCHAML de pression resultat / 0 si echec
  17. C    
  18. C-----------------------------------------------------------------------
  19.       SUBROUTINE PRCHL2(IPMOD1,MOT1,XP,IPCHE1)
  20. C
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24.       CHARACTER*(*) MOT1
  25.  
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. -INC SMCHAML
  30. -INC SMMODEL
  31. -INC SMCOORD
  32. C
  33.       IPCHE2=0
  34. C
  35. C  Initialisation du MCHAML :
  36.       MMODEL=IPMOD1
  37.       SEGACT,MMODEL
  38.       N1 = MMODEL.KMODEL(/1)
  39.       N3 = 6
  40.       L1 = 11
  41.       SEGINI,MCHEL1
  42.       MCHEL1.TITCHE='CONTRAINTES'
  43.       MCHEL1.IFOCHE=IFOUR
  44. C
  45. C Boucle sur les sous-zones
  46.       DO I=1,N1
  47.          IMODEL=MMODEL.KMODEL(I)
  48.          SEGACT,IMODEL
  49. C
  50. C On verifie que MOT1 est un nom de composante de contrainte
  51.          NOMID=IMODEL.LNOMID(4)
  52.          SEGACT,NOMID
  53.          NCOMP=NOMID.LESOBL(/2)
  54.          DO J=1,NCOMP
  55.             IF (MOT1.EQ.NOMID.LESOBL(J)) THEN
  56.                GOTO 10
  57.             ENDIF
  58.          ENDDO
  59. C Pas de composante de nom MOT1 : on sort !
  60.          MOTERR(1:4)=MOT1
  61.          CALL ERREUR(197)
  62.          GOTO 999
  63. C
  64.  10      CONTINUE
  65. C
  66. C On remplit MCHEL1
  67.         MCHEL1.CONCHE(I)=IMODEL.CONMOD(1:24)
  68.         MCHEL1.IMACHE(I)=IMODEL.IMAMOD
  69.         MCHEL1.INFCHE(I,1)=0
  70.         MCHEL1.INFCHE(I,2)=0
  71.         MCHEL1.INFCHE(I,3)=NIFOUR
  72.         MCHEL1.INFCHE(I,4)=IMODEL.INFMOD(7)
  73.         MCHEL1.INFCHE(I,5)=0
  74.         MCHEL1.INFCHE(I,6)=5
  75. C
  76. C Definition du ICHAML
  77.          N2 = 1
  78.          SEGINI,MCHAM1
  79.          MCHEL1.ICHAML(I)=MCHAM1
  80.          MCHAM1.NOMCHE=MOT1
  81.          MCHAM1.TYPCHE='REAL*8'
  82. C
  83. C et du MELVAL
  84.          N1PTEL = 1
  85.          N1EL   = 1
  86.          N2PTEL = 0
  87.          N2EL   = 0
  88.          SEGINI,MELVAL
  89.          MELVAL.VELCHE(1,1)=XP
  90.          MCHAM1.IELVAL=MELVAL
  91. C
  92.          SEGDES,MELVAL,MCHAM1
  93.       ENDDO
  94. C
  95.       IPCHE1=MCHEL1
  96. C
  97.  999  SEGDES,MCHEL1
  98. C
  99.       RETURN
  100.       END
  101.  
  102.  
  103.  
  104.  

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