Télécharger la plaquette Cast3M


fr en

Télécharger tresca.eso

Retour à la liste

Numérotation des lignes : 

tresca

  1. C TRESCA    SOURCE    CB215821  23/07/12    21:15:12     11704          
  2.       SUBROUTINE TRESCA
  3. *---------------------------------------------------------------------
  4. *
  5. *              CALCUL DE LA CONTRAINTE DE TRESCA
  6. *
  7. *     SYNTAXE :
  8. *     ---------
  9. *     CHAL1 = TRESCA MODL1 SIG1 (CAR1) (MOT1)
  10. *
  11. *       MODL1 =  Modele de calcul (type  MMODEL)
  12. *       SIG1  =  Objet de type MCHALM  de CONTRAINTES
  13. *       CAR1  =  Objet de type MCHAML  de CARACTERISTIQUE
  14. *       MOT1  =  Mot qui indique pour les coques ou l'on calcule les
  15. *                contraintes
  16. *       CHEL1 =  Objet de type MCHAML de sous type SCALAIRE (aux pts
  17. *                de gauss pour les contraintes)
  18. *
  19. * Passage aux nouveaux chamelem par jm CAMPENON le 04/91
  20. *
  21. *--------------------------------------------------------------------
  22. *
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25. *
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. *
  29.       CHARACTER*4 MOMOT
  30. *
  31.       MOMOT='    '
  32.       IMIL=0
  33. *
  34.       IPMODL=0
  35.       IPCHE1=0
  36.       IPCHE2=0
  37. *
  38.       CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
  39.       CALL ACTOBJ ('MMODEL  ',IPMODL,1)
  40.       IF(IERR.NE.0) RETURN
  41. *
  42.       CALL LIROBJ('MCHAML  ',IPCHE1,1,IRET)
  43.       CALL ACTOBJ('MCHAML  ',IPCHE1,1)
  44.       IF(IERR.NE.0) RETURN
  45. *
  46. *     On lit les CARACTERISTIQUES (facultatif)
  47. *
  48.       CALL LIROBJ('MCHAML  ',IPCHE2,0,IRET)
  49.       IF(IRET .EQ. 1) CALL ACTOBJ('MCHAML  ',IPCHE2,1)
  50.       IF(IERR.NE.0) RETURN
  51. *
  52. *     Lecture d'un mot si necessaire
  53. *
  54.       CALL LIRCHA(MOMOT,0,IRET)
  55.       IF(IERR.NE.0) RETURN
  56.       IF (IRET.NE.0) THEN
  57.        IF (MOMOT.EQ.'SUPE') IMIL=1
  58.        IF (MOMOT.EQ.'INFE') IMIL=-1
  59.        IF (MOMOT.NE.'SUPE'.AND.MOMOT.NE.'INFE'.AND.MOMOT.NE.'MOYE') THEN
  60.           MOTERR=MOMOT
  61.           CALL ERREUR(7)
  62.           RETURN
  63.        ENDIF
  64.       ENDIF
  65. *
  66.       IRET=0
  67.       CALL TRESK(IPMODL,IPCHE1,IPCHE2,IMIL,IPSCAL,IRET)
  68. *
  69.       IF(IRET.NE.0) THEN
  70.         CALL ACTOBJ('MCHAML  ',IPSCAL,1)
  71.         CALL ECROBJ('MCHAML  ',IPSCAL)
  72.       ENDIF
  73.       END
  74.  
  75.  
  76.  
  77.  

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