Télécharger la plaquette Cast3M


fr en

Télécharger vmises.eso

Retour à la liste

Numérotation des lignes : 

vmises

  1. C VMISES    SOURCE    OF166741  24/05/02    21:15:04     11928          
  2.  
  3. C=======================================================================
  4. C=                            V M I S E S                              =
  5. C=                            -----------                              =
  6. C=                                                                     =
  7. C=  OPERATEUR CAST3M "VMISES" :                                        =
  8. C=  ---------------------------                                        =
  9. C=   CHAM1  =  'VMIS'  MODL1  SIG1  ( CAR1 ) ;                         =
  10. C=                                                                     =
  11. C=  Cet operateur sert a calculer les contraintes de Von Mises.        =
  12. C=                                                                     =
  13. C=  ARGUMENTS :                                                        =
  14. C=  -----------                                                        =
  15. C=   MODL1   (MMODEL)   Modele associe a la structure etudiee          =
  16. C=   SIG1    (MCHAML)   Champ de contraintes                           =
  17. C=   CAR1    (MCHAML)   Caracteristiques geometriques     (facultatif) =
  18. C=                      Sous-type 'CARACTERISTIQUES'                   =
  19. C=                                                                     =
  20. C=  RESULTAT :                                                         =
  21. C=  ----------                                                         =
  22. C=   CHAM1   (MCHAML)   Contrainte equivalente de Von Mises            =
  23. C=======================================================================
  24.  
  25.       SUBROUTINE VMISES
  26.  
  27.       IMPLICIT INTEGER(I-N)
  28.       IMPLICIT REAL*8 (A-H,O-Z)
  29.  
  30.  
  31. -INC PPARAM
  32. -INC CCOPTIO
  33.  
  34. C  1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  35. C ==========================================
  36. C  1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  37. C =====
  38.       CALL LIROBJ('MMODEL  ',IPMODL,1,iOK)
  39.       CALL ACTOBJ('MMODEL  ',IPMODL,1)
  40.       IF (IERR.NE.0) RETURN
  41.  
  42. C =====
  43. C  1.2 - Lecture OBLIGATOIRE du champ de contraintes (IPCHA1)
  44. C =====
  45.       IPCHA1 = 0
  46.       CALL LIROBJ('MCHAML  ',IPIN,1,iOK)
  47.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  48.       IF (IERR.NE.0) RETURN
  49.       CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  50.       IF(IR   .NE. 1) CALL ERREUR(KER)
  51.       IF(IERR .NE. 0) RETURN
  52. C =====
  53. C  1.3 - Lecture FACULTATIVE des caracteristiques geometriques (IPCHA2)
  54. C =====
  55.       IPCHA2 = 0
  56.       CALL LIROBJ('MCHAML  ',IPIN,0,iOK)
  57.       IF (IERR.NE.0) RETURN
  58.       IF (iOK .EQ. 1) THEN
  59.         CALL ACTOBJ('MCHAML  ',IPIN,1)
  60.         CALL REDUAF(IPIN,IPMODL,IPCHA2,0,IR,KER)
  61.         IF(IR   .NE. 1) CALL ERREUR(KER)
  62.         IF(IERR .NE. 0) RETURN
  63.       ENDIF
  64. C =====
  65. C  1.4 - Mise en ordre des MCHAML lus
  66. C        1 -> IPCHE1 = Contraintes, 2 -> IPCHE2 = Carac. geom.
  67. C =====
  68.       IPCHE1=0
  69.       IPCHE2=0
  70.       CALL RNGCHA(IPCHA1,IPCHA2,'CONTRAINTES','CARACTERISTIQUES',
  71.      &            IPCHE1,IPCHE2)
  72.       IF (IERR.NE.0) RETURN
  73.  
  74.  
  75. C  2 - CALCUL DU MCHAML DE CONTRAINTE DE VON MISES
  76. C =================================================
  77.       isouc=0
  78.       CALL VMISPO(IPMODL,IPCHE1,IPCHE2,IPCHE3,iRet,isouc)
  79.  
  80. C  3 - ECRITURE DU MCHAML RESULTAT
  81. C =================================
  82.       IF (iRet.EQ.1) THEN
  83.         CALL ACTOBJ('MCHAML  ',IPCHE3,1)
  84.         CALL ECROBJ('MCHAML  ',IPCHE3)
  85.       ENDIF
  86.  
  87.       END
  88.  
  89.  
  90.  

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