Esope source files

   1 : C VMISES    SOURCE    CHAT      06/08/24    22:14:28     5529           
   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 : -INC CCOPTIO
  31 : 
  32 : C  1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  33 : C ==========================================
  34 : C  1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  35 : C =====
  35 :       CALL LIROBJ('MMODEL',IPMODL,1,iOK)
  36 :       IF (IERR.NE.0) RETURN
  37 : C =====
  38 : C  1.2 - Lecture OBLIGATOIRE du champ de contraintes (IPCHA1)
  39 : C =====
  39 :       CALL LIROBJ('MCHAML',IPCHA1,1,iOK)
  40 :       IF (IERR.NE.0) RETURN
  41 : C =====
  42 : C  1.3 - Lecture FACULTATIVE des caracteristiques geometriques (IPCHA2)
  43 : C =====
  43 :       CALL LIROBJ('MCHAML',IPCHA2,0,iOK)
  44 :       IF (IERR.NE.0) RETURN
  45 : C =====
  46 : C  1.4 - Mise en ordre des MCHAML lus
  47 : C        1 -> IPCHE1 = Contraintes, 2 -> IPCHE2 = Carac. geom.
  48 : C =====
  49 :       IPCHE1=0
  50 :       IPCHE2=0
  50 :       CALL RNGCHA(IPCHA1,IPCHA2,'CONTRAINTES','CARACTERISTIQUES',
  51 :      .            IPCHE1,IPCHE2)
  52 :       IF (IERR.NE.0) RETURN
  53 : 
  54 : C  2 - CALCUL DU MCHAML DE CONTRAINTE DE VON MISES
  55 : C =================================================
  56 :       iConv=0
  56 :       CALL VMISPO(IPMODL,IPCHE1,IPCHE2,iConv,IPCHE3,iRet)
  57 : 
  58 : C  3 - ECRITURE DU MCHAML RESULTAT
  59 : C =================================
  59 :       IF (iRet.EQ.1) CALL ECROBJ('MCHAML',IPCHE3)
  60 : 
  61 :       RETURN
  62 :       END
  63 : 
  64 :  
  65 :  
  66 :