Télécharger la plaquette Cast3M


fr en

Télécharger sste.eso

Retour à la liste

Numérotation des lignes : 

sste

  1. C SSTE      SOURCE    OF166741  24/10/07    21:15:49     12016          
  2.  
  3. *---------------------------------------------------------------------
  4. *                    ECOULEMENT  INELASTIQUE
  5. *---------------------------------------------------------------------
  6. *
  7. *    SYNTAXE POUR LES LOIS PLASTIQUES:
  8. *
  9. *    SIGF VARF DEPIN RI1 = 'SUBSTEP'  MODL SIG0 VAR0 DEPST CARAC
  10. *                                    (PRECIS) (NMAXSSTEPS) (NITMAX)
  11. *                                    (NNUMER) (DELTAX);
  12. * IN=
  13. *  MMODEL     |    MODL    OBJET MODELE
  14. *  MCHAML     |    SIG0    CONTRAINTES AU DEBUT DU PAS
  15. *             |    VAR0    VARIABLES INTERNES AU DEBUT DU PAS
  16. *             |    DEPST   INCREMENT DE DEFORMATIONS TOTALES
  17. *             |    CARAC   CONSTANTES DU MATERIAU
  18. *  FLOTTANT   |    PRECIS  PRECISION POUR ITERATIONS INTERNES
  19. *  INTEGER    |    NMAXSSTEPS Maximum number of substeps
  20. *  INTEGER    |    NITMAX  Maximum number of iterations
  21. *  INTEGER    |    NNUMER  Numerical differentiation type (MRSLade)
  22. *  FLOTTANT   |    DELTAX  Numerical differentiation parameter (MRSLade)
  23. *
  24. * OUT=
  25. *  MCHAML     |    SIGF    CONTRAINTES A LA FIN DU PAS
  26. *             |    VARF    VARIABLES INTERNES A LA FIN DU PAS
  27. *             |    DEPIN   INCREMENT DE DEFORMATIONS INELASTIQUES
  28. *             |    RI1     CONSISTENT TANGENT MATRIX
  29. *
  30. *---------------------------------------------------------------------
  31.  
  32.       SUBROUTINE SSTE
  33.  
  34.       IMPLICIT INTEGER(I-N)
  35.       IMPLICIT REAL*8(A-H,O-Z)
  36.  
  37. -INC PPARAM
  38. -INC CCOPTIO
  39.  
  40. -INC SMCHAML
  41.  
  42. * comprovar que es def plana o axisimetrico
  43.       IF (IFOUR.gt.0 .or. IFOUR.lt.-1) then
  44.         write(ioimp,*) ' Formulacion no disponible'
  45.         call erreur(21)
  46.         return
  47.       endif
  48. *
  49. * Lectura de datos
  50. *
  51. * modelo
  52.       CALL LIROBJ('MMODEL  ',IPMODL,1,IRT)
  53.       CALL ACTOBJ('MMODEL  ',IPMODL,1)
  54.       IF (IERR.NE.0) RETURN
  55.  
  56. * tensiones iniciales
  57.       CALL LIROBJ('MCHAML  ',IPIN,1,IRT)
  58.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  59.       IF (IERR.NE.0) RETURN
  60.       CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  61.       IF(IR   .NE. 1) CALL ERREUR(KER)
  62.       IF(IERR .NE. 0) RETURN
  63.  
  64.       MCHELM=IPCHE1
  65. c! Usa L1, N1, N3
  66.       IF(MCHELM.TITCHE.NE.'CONTRAINTES')THEN
  67.         MOTERR(1:32)='CONTRAINTES'
  68.         CALL ERREUR(565)
  69.         RETURN
  70.       ENDIF
  71.  
  72. * variables internas iniciales
  73.       CALL LIROBJ('MCHAML  ',IPIN,1,IRT)
  74.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  75.       IF (IERR.NE.0) RETURN
  76.       CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  77.       IF(IR   .NE. 1) CALL ERREUR(KER)
  78.       IF(IERR .NE. 0) RETURN
  79.  
  80.       MCHELM=IPCHE2
  81.       IF(TITCHE.NE.'VARIABLES INTERNES')THEN
  82.         MOTERR(1:32)='VARIABLES INTERNES'
  83.         CALL ERREUR(565)
  84.         RETURN
  85.       ENDIF
  86.  
  87. * incremento de deformacion
  88.       CALL LIROBJ('MCHAML  ',IPIN,1,IRT)
  89.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  90.       IF (IERR.NE.0) RETURN
  91.       CALL REDUAF(IPIN,IPMODL,IPCHE4,0,IR,KER)
  92.       IF(IR   .NE. 1) CALL ERREUR(KER)
  93.       IF(IERR .NE. 0) RETURN
  94.  
  95.       MCHELM=IPCHE4
  96.       IF(TITCHE.NE.'DEFORMATIONS')THEN
  97.         MOTERR(1:32)='DEFORMATIONS'
  98.         CALL ERREUR(565)
  99.         RETURN
  100.       ENDIF
  101.  
  102. * caracteristicas materiales y geometricas
  103.       CALL LIROBJ('MCHAML',IPIN,1,IRT)
  104.       CALL ACTOBJ('MCHAML',IPIN,1)
  105.       IF (IERR.NE.0) RETURN
  106.       CALL REDUAF(IPIN,IPMODL,IPCAR,0,IR,KER)
  107.       IF(IR   .NE. 1) CALL ERREUR(KER)
  108.       IF(IERR .NE. 0) RETURN
  109.  
  110.       MCHELM=IPCAR
  111.       IF(TITCHE.NE.'CARACTERISTIQUES')THEN
  112.         MOTERR(1:32)='CARACTERISTIQUES'
  113.         CALL ERREUR(565)
  114.         RETURN
  115.       ENDIF
  116.  
  117. * precision para iteraciones internas
  118.       CALL LIRREE(PRECIS,0,IRT)
  119.       IF(IRT.EQ.0) PRECIS=1.D-10
  120.  
  121. * num max de substeps
  122.       CALL LIRENT(NMAXSSTEPS,0,IRT)
  123.       IF (IRT.EQ.0) NMAXSSTEPS=1
  124.  
  125. * num max de iteraciones internas
  126.       CALL LIRENT(NITMAX,0,IRT)
  127.       IF (IRT.EQ.0) NITMAX=25
  128.  
  129. * tipo de derivacion numerica
  130.       CALL LIRENT(NNUMER,0,IRT)
  131.       IF (IRT.EQ.0) NNUMER=0
  132.  
  133. * delta de x relativo para las derivadas numericas
  134.       CALL LIRREE(DELTAX,0,IRT)
  135.       IF (IRT.EQ.0) DELTAX=0.D0
  136.  
  137.       CALL SSTE1(IPMODL,IPCHE1,IPCHE2,IPCHE4,IPCAR,
  138.      .              PRECIS,NITMAX,NMAXSSTEPS,NNUMER,DELTAX,
  139.      .              IPCHE7,IPCHE8,IPCHE9,IPRIGI)
  140.  
  141.       IF(IERR.EQ.0)THEN
  142.         CALL ACTOBJ('MCHAML  ',IPCHE9,1)
  143.         CALL ACTOBJ('MCHAML  ',IPCHE8,1)
  144.         CALL ACTOBJ('MCHAML  ',IPCHE7,1)
  145.         CALL ACTOBJ('RIGIDITE',IPRIGI,1)
  146.  
  147.         CALL ECROBJ('MCHAML  ',IPCHE9)
  148.         CALL ECROBJ('MCHAML  ',IPCHE8)
  149.         CALL ECROBJ('MCHAML  ',IPCHE7)
  150.         CALL ECROBJ('RIGIDITE',IPRIGI)
  151.       ENDIF
  152.  
  153.       END
  154.  
  155.  
  156.  

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