Télécharger la plaquette Cast3M


fr en

Télécharger ecou20.eso

Retour à la liste

Numérotation des lignes : 

ecou20

  1. C ECOU20    SOURCE    CB215821  24/04/12    21:15:41     11897          
  2.       SUBROUTINE ECOU20(MATE,INPLAS,MELE,IPMAIL,NBPTEL,IMAT,ICAR,
  3.      1     NUMAT,NUCAR,IVASTR,IVARI,IVADEF,IVADET,
  4.      1     IVADS,IVAMAT,IVACAR,
  5.      2     IPH1,IPH2,IPH3,ITHHER,LHOOK,NSTRS,NVARI,NMATT,NCARR,
  6.      3     CMATE,PRECIS,JECHER,IPOTAB,ISTEP,NPINT,JNOID,LOGSUC,
  7.      4     N2EL,N2PTEL,NBNO,NBPGAU,LW,IVASTF,IVARIF,IVADEP,KERRE)
  8. ***********************************************************************
  9. *      ecoulement inelastique   appele par ecoul1
  10. c ppu  modif pour les materiaux unidirectionels en plastique
  11. *      MATERIAUX: - VISCOPLASTIQUES ET FLUAGE INTEGRES PAR CONSTI
  12. ***********************************************************************
  13. * entrees :
  14. *
  15. *  mate   = numero de materiau elastique
  16. *  inplas = numero de materiau inelastique
  17. *  mele   = numero  element fini
  18. *  ipmail = pointeur du maillage
  19. *  nbptel = nombre de points par element
  20. *  imat   = pointeur sur un segment mptval de materiau (utilise par calsig)
  21. *  icar   = pointeur sur un segment mptval de caracteristiques
  22. *           geometriques (utilise par calsig)
  23. *  numat  = nb de composantes du melval de imat
  24. *  nucar  = nb de composantes du melval de icar
  25. *  ivastr =pointeur sur un segment mptval de contraintes
  26. *  ivari  =pointeur sur un segment mptval de variables internes
  27. *  ivadef =pointeur sur un segment mptval de deformations
  28. *  ivads  =pointeur sur un segment mptval de contraintes  (increments)
  29. *  ivamat =pointeur sur un segment mptval de materiau
  30. *  ivacar =pointeur sur un segment mptval de cacarteristiques geometrique
  31. *  iph1  = pointeur sur un mchaml de temperatures au debut du pas
  32. *  iph2  = pointeur sur un mchaml de temperatures a la fin du pas
  33. *  iph3  = pointeur sur un mchaml de temperatures de reference
  34. *  ithher = 0 si pas de chargement thermique
  35. *         = 1 si chargement thermique mais materiau constant
  36. *         = 2 si chargement thermique et mat. dependant de la temperature
  37. *  ipch1,ipch2,ipch3,ithher ne servent que pour les materiaux
  38. *        endommageables de lemaitre quand ils dependent de la temperature
  39. *  lhook  =taille de la matrice de hooke
  40. *  nstrs  =nombre de composantes de contraintes
  41. *  nvari  =nombre de composantes de variables internes
  42. *  nmatt  =nombre de composnates de proprietes de materiau
  43. *  ncarr  =nombre de composnates de caracteristiques geometriques
  44. *  cmate  =nom du materiau
  45. *  precis =precision dans les iterations internes
  46. *  jecher =0 ou 1 pour action dans ecoule
  47. *  jnoid  =0 ou 1 pour action dans ecoule
  48. *  ipotab =pointeur sur segment table
  49. *  istep  =indicateur d'action pour calcul nonlocal
  50. *        =0 dans le cas d'un calcul local (normal)
  51. *        =1 ou 2 dans le cas d'un calcul nonlocal
  52. *        =1 pour calcul des fonctions seuil uniquement
  53. *        =2 pour calcul des variables dissipatives a partir
  54. *           des fonctions seuil moyennees prealablement par nloc
  55. *
  56. * sorties :
  57. *  ivastf  =pointeur sur un segment mptval de contraintes
  58. *  ivarif  =pointeur sur un segment mptval de variables internes
  59. *  ivadep  =pointeur sur un segment mptval de deformations inelastiques
  60. *  kerre   =indicateur d'erreur
  61. *
  62. *  p dowlatyari  fev. 1992
  63. *
  64. *  c. la borderie fev 92 restructuration et reecriture de certains
  65. *                        passages pour une meilleure lisibilite
  66. *
  67. *                 avril 92 ajout istep pour le non local
  68. *                 dec 92 modif pour poutres timoschenko
  69. *
  70. ************************************************************************
  71.       IMPLICIT INTEGER(I-N)
  72.       IMPLICIT REAL*8(A-H,O-Z)
  73. *
  74.  
  75. -INC PPARAM
  76. -INC CCOPTIO
  77. -INC SMCHAML
  78. -INC SMELEME
  79. -INC SMCOORD
  80. -INC SMMODEL
  81. -INC SMINTE
  82. -INC CCHAMP
  83. -INC CECOU
  84. c=======================================================================
  85. c  la variable kerre regit les impressions d erreurs dans plast
  86. c      toutes erreurs de ecoule gerees dans ce sous programme
  87. c   kerre=0 tout ok
  88. c        de 1 a 6  s aligner sur valeurs donnees par ecoinc
  89. c     =   7  un element tuyau a une epaisseur nulle
  90. c     =  21  on ne trouve pas d intersection avec la surface de charge
  91. c     =  22  sig0 a l exterieur de la surface de charge
  92. c
  93. c             anomalies avec la courbe de traction
  94. c     =  30 limite elastique nulle
  95. c     =  31 trop de points
  96. c     =  32 pas assez de points
  97. c     =  33 pente incorrecte
  98. c     =  34 module d'young nul
  99. c     =  35 manque l'origine
  100. c     =  36 pente a l'origine non egale a e
  101. c     =  37 manque la courbe de traction
  102. c     =  38 nu devrait etre nul
  103. c
  104. c     =  48 donnees erronnees pour drucker-prager
  105. c     =  49 matrice singuliere dans iter internes drucker-prager
  106. c     =  51 pb dans drucker prager option non disponible
  107. c     =  52 pb dans drucker prager donnees incompatibles
  108. c     =  53 pb dans drucker prager solution impossible
  109. c     =  54 les valeurs admissibles pour istep sont 0 1 ou 2
  110. c     =  55 modele non implante en non local
  111. c     =  56 probleme dans l'integration du modele mazars
  112. c     =  57  ....
  113. c     =  58  ....
  114. c     =  59  ....
  115. c     =  60 pb donnees du cam-clay
  116. c
  117. c     =  99 cas non encore disponible
  118. c=======================================================================
  119. *
  120.       SEGMENT MPTVAL
  121.         INTEGER IPOS(NS)    ,NSOF(NS)
  122.         INTEGER      IVAL(NCOSOU)
  123.         CHARACTER*16 TYVAL(NCOSOU)
  124.       ENDSEGMENT
  125. *
  126.       SEGMENT WRK0
  127.        REAL*8 XMAT(NCXMAT)
  128.       ENDSEGMENT
  129. *
  130.       SEGMENT WR00
  131.        CHARACTER*16 TYMAT(NCXMAT)
  132.        REAL*8 XMAT1(NCXMAT),XMAT2(NCXMAT)
  133.       ENDSEGMENT
  134. *
  135.       SEGMENT WRK1
  136.         REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  137.         REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  138.         REAL*8 DEFP(NSTRS),XCAR(ICARA)
  139.       ENDSEGMENT
  140. *
  141.       SEGMENT WRK2
  142.         REAL*8 TRAC(LTRAC)
  143.       ENDSEGMENT
  144. *
  145.       SEGMENT WRK22
  146.       REAL*8 XXE(3,NBNN)
  147.       ENDSEGMENT
  148. *
  149.       SEGMENT WRK3
  150.         REAL*8 WORK(LW),WORK2(LW2)
  151.       ENDSEGMENT
  152. *
  153.       SEGMENT WRK4
  154.         REAL*8 XE(3,NBBB)
  155.       ENDSEGMENT
  156. *
  157.       SEGMENT WRK5
  158.         REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
  159.       ENDSEGMENT
  160. *
  161.       SEGMENT WRK6
  162.         REAL*8 BB(NSTRS,NNVARI),R(NSTRS),XMU(NSTRS)
  163.         REAL*8 S(NNVARI),QSI(NNVARI),DDR(NSTRS),BBS(NSTRS)
  164.         REAL*8 SIGMA(NSTRS),SIGGD(NSTRS),XMULT(NSTRS),PROD(NSTRS)
  165.       ENDSEGMENT
  166. *
  167.       SEGMENT WRK7
  168.         REAL*8 F(NCOURB,2),W(NCOURB),TRUC(NCOURB)
  169.       ENDSEGMENT
  170. *
  171.       SEGMENT WRK8
  172.         REAL*8 DD(NSTRS,NSTRS),DDV(NSTRS,NSTRS),DDINV(NSTRS,NSTRS)
  173.       ENDSEGMENT
  174. *
  175.       SEGMENT WRK9
  176.         REAL*8 YOG(NYOG),YNU(NYNU),YALFA(NYALFA),YSMAX(NYSMAX)
  177.         REAL*8 YN(NYN),YM(NYM),YKK(NYKK),YALFA1(NYALF1)
  178.         REAL*8 YBETA1(NYBET1),YR(NYR),YA(NYA),YKX(NYKX),YRHO(NYRHO)
  179.         REAL*8 SIGY(NSIGY)
  180.         INTEGER NKX(NNKX)
  181.       ENDSEGMENT
  182. *
  183.       SEGMENT WR10
  184.         INTEGER IABLO1(NTABO1)
  185.         REAL*8  TABLO2(NTABO2)
  186.       ENDSEGMENT
  187. *
  188.       SEGMENT WR11
  189.         INTEGER IABLO3(NTABO3)
  190.         REAL*8  TABLO4(NTABO4)
  191.       ENDSEGMENT
  192. *
  193.       SEGMENT WTRAV
  194.         REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  195.         REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  196.         REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  197.         REAL*8 XLOC(3,3),XGLOB(3,3)
  198.         REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  199.       ENDSEGMENT
  200. *
  201.       SEGMENT WPOUT
  202.          REAL*8  X(2),Y(2),Z(2)
  203.       ENDSEGMENT
  204.       LOGICAL LOGVIS,LOGIN,LOGRE,LOGSUC
  205.       LOGICAL LUNI1,LUNI2
  206.       DIMENSION BID(6),BID2(6),CRIGI(12),CMASS(12)
  207. *
  208.       CHARACTER*72 CHARRE
  209.       CHARACTER*8  CMATE
  210. c
  211. *
  212. *   mise à disposition des temperatures tini tfin tref
  213. *   aux points de gauss
  214. *
  215.         TETA1=-1.E35
  216.          TETA2=-1.E35
  217.          TETREF=-1.E35
  218.          TREFA=-1.E35
  219.          IF (ITHHER.EQ.1.OR.ITHHER.EQ.2) THEN
  220.             MCHAM3=IPH1
  221.             MCHAM4=IPH2
  222.             MCHAM5=IPH3
  223.             SEGACT MCHAM3
  224.             SEGACT MCHAM4
  225.             SEGACT MCHAM5
  226.             MELVA3=MCHAM3.IELVAL(1)
  227.             MELVA4=MCHAM4.IELVAL(1)
  228.             MELVA5=MCHAM5.IELVAL(1)
  229.             SEGACT MELVA3
  230.             SEGACT MELVA4
  231.             SEGACT MELVA5
  232.          ENDIF
  233. c
  234. c   Initialisations de variables
  235. c---------------------------------
  236. c      - mise à zéro des variables du commun NECOU si besoin
  237. c      - modèles viscoplastiques:
  238. c                  . on récupère le pas de temps
  239. c                  . on récupère le nombre maximal de sous-pas
  240. c                  . on met IND=1
  241. c      - initialisation des dimensions des tableaux des segments
  242. c   Sorties: en plus du commun NECOU, on range les autres données
  243. c   initialisées dans les COMMON IECOU et XECOU
  244. c   Sauf pour KERRE,LW,LOGVIS,LUNI1 et LUNI2 qui sont sortis comme
  245. c   argument de DEFINI
  246. c
  247.       CALL DEFINI(MELE,NCARR,NSTRS,NMATT,CMATE,MATE,
  248.      .           ISTEP,INPLAS,NPINT,IPOTAB,IVADEF,
  249.      .           IPMAIL,IVAMAT,
  250.      .           ITHHER,NUMAT,NUCAR,LOGVIS,
  251.      .           LUNI1,LUNI2,LW,KERRE)
  252.       IF (KERRE.EQ.999) RETURN
  253. c
  254.       IF ((CMATE.EQ.'ORTHOTRO'.OR.CMATE.EQ.'ANISOTRO'.OR.
  255.      1     CMATE.EQ.'UNIDIREC').AND.(MFR.EQ.1.OR.MFR.EQ.31
  256.      1     .OR.MFR.EQ.33)) THEN
  257.         CALL RESHPT(1,NBNO,IELE,MELE,NPINT,IPTR1,IRT1)
  258.         MINTE2=IPTR1
  259.         SEGACT MINTE2
  260.         SEGINI WRK22
  261.       ENDIF
  262. c
  263.       IF (LOGVIS) SEGINI WRK8
  264.       SEGINI WRK0,WR00,WRK1,WRK2,WRK3,WRK5
  265.       IF(MFR.EQ.7.OR.MFR.EQ.13.OR.LUNI1)THEN
  266.          SEGINI WRK4
  267.       ENDIF
  268. c
  269.       SEGINI WTRAV
  270. *
  271. *
  272. *       boucle sur les elements
  273. *
  274.       DO 1000 IB=1,NBELEM
  275.  
  276. *
  277. *  Matériaux orthotropes, anisotropes et unidirectionnels
  278. *  en formulation massive:
  279. *     - on cherche  les coordonnees des noeuds de l element ib
  280. *     - calcul des axes locaux
  281. *  Cas particulier de l'ACIER_UNI
  282. *
  283.       CALL DEFROT(CMATE,MFR,NBNN,IB,MELE,LUNI1,IPTR1,
  284.      .     MELEME,WRK4,WRK22,WTRAV)
  285. *
  286. *
  287. *       boucle sur les points de gauss
  288. *
  289.          DO 1100 IGAU=1,NBPTEL
  290. *
  291. *  -recuperation de valmat et de valcar
  292. *  -on recupere les contraintes initiales
  293. *  -on recupere les variables internes
  294. *  -on recupere les deformations inelastiques initiales si besoin
  295. *  -on recupere les increments de deformations totales
  296. *  -on cherche la section de l'element ib
  297. *  -prise en compte de l'epaisseur et de l'excentrement
  298. *  dans le cas des coques minces avec ou sans cisaillement
  299. *  transverse
  300. *
  301.       CALL DEFVAL(NUMAT1,NBPTEL,NDEF,
  302.      .     IMAT,IVACAR,ICAR,IVASTR,IVARI,IVADEF,IVADET,
  303.      .     IVADS,MFR,CMATE,INPLAS,IB,IGAU,IND,
  304.      .     WTRAV,WRK1,WRK5,SECT,EPAIST)
  305. *
  306. *     on recupere les constantes du materiau
  307. *     calcul des contraintes effectives en milieu poreux
  308. *
  309.       CALL DEFMAT(NMATT,NSTRS,MFR,MELE,INPLAS,
  310.      .        IVAMAT,IB,IGAU,CMATE,MATE,LUNI1,LUNI2,
  311.      .        WRK1,WRK5,WRK0,WR00,WTRAV,CMASS,CRIGI,COB,XMOB,
  312.      .        BID,BID2,KERR0)
  313.       IF (KERR0.EQ.99) THEN
  314.          KERRE=99
  315.          GOTO 1000
  316.       ELSE IF (KERR0.EQ.10) THEN
  317.          GOTO 1000
  318.       ENDIF
  319. *
  320. * >>>>>>>>>>   fin du traitement du materiau
  321. *
  322. *        on recupere les caracteristiques geometriques
  323. *
  324.       CALL DEFCAR(NCARR,IB,IGAU,MFR,MELE,IVACAR,
  325.      .           WRK1)
  326. *
  327. *
  328. *          quelques impressions si iimpi = 99
  329. *
  330. *      IF(IIMPI.EQ.99) THEN
  331. *         WRITE(IOIMP,66770) IB,IGAU
  332. *66770    format(////////2x,'element  ',i6,2x,'point  ',i3//)
  333. *         WRITE(IOIMP,66771) MATE,INPLAS
  334. *66771    format('0  mate=',i4,2x,'inplas=',i4/)
  335. *         WRITE(IOIMP,66772) (SIG0(I),I=1,NSTRS)
  336. *66772    format(2x,'  sig0 '/(6(1x,1pe12.5)))
  337. *         WRITE(IOIMP,66773) (VAR0(I),I=1,NVARI)
  338. *66773    format(2x,'  var0 '/(6(1x,1pe12.5)))
  339. *         WRITE(IOIMP,66774) (DEPST(I),I=1,NSTRS)
  340. *66774    format(2x,'  depst '/(6(1x,1pe12.5)))
  341. *         WRITE(IOIMP,66775) (XMAT(I),I=1,NMATT)
  342. *66775    format(2x,'  xmat  '/(6(1x,1pe12.5)))
  343. *         IF(IVACAR.NE.0)THEN
  344. *            WRITE(IOIMP,66776) (XCAR(I),I=1,ICARA)
  345. *66776       format(2x,'  xcar  '/(6(1x,1pe12.5)))
  346. *         ENDIF
  347. *      ENDIF
  348. *
  349. *   mise à disposition des temperatures tini tfin tref
  350. *   aux points de gauss
  351. *
  352.          IF (ITHHER.EQ.1.OR.ITHHER.EQ.2) THEN
  353.             IBMN=MIN(IB,MELVA3.VELCHE(/2))
  354.             IGMN=MIN(IGAU,MELVA3.VELCHE(/1))
  355.             TETA1=MELVA3.VELCHE(IGMN,IBMN)
  356.             IBMN=MIN(IB,MELVA4.VELCHE(/2))
  357.             IGMN=MIN(IGAU,MELVA4.VELCHE(/1))
  358.             TETA2=MELVA4.VELCHE(IGMN,IBMN)
  359.             IBMN=MIN(IB,MELVA5.VELCHE(/2))
  360.             IGMN=MIN(IGAU,MELVA5.VELCHE(/1))
  361.             TETREF=MELVA5.VELCHE(IGMN,IBMN)
  362.          ENDIF
  363. *
  364. *
  365. *---------------------------------------------------------------------
  366. *
  367. *                  ecoulement selon les modeles
  368. *
  369. *---------------------------------------------------------------------
  370. *
  371. *     modeles de viscoplasticite integres par consti
  372. *
  373.             IF ( INPLAS .EQ. 17 .OR.
  374.      2                (INPLAS .GE. 19 .AND. INPLAS .LE. 25) .OR.
  375.      4                INPLAS .EQ. 61 .OR.
  376.      4                INPLAS .EQ. 63 .OR.
  377.      1                INPLAS .EQ. 53 .OR. INPLAS .EQ. 102 .OR.
  378.      8                INPLAS .EQ. 44 .OR. INPLAS .EQ. 76  .OR.
  379.      9                INPLAS .EQ. 45 .OR. INPLAS .EQ. 77  .OR.
  380.      9                INPLAS .EQ. 84 .OR. INPLAS .EQ. 85  .OR.
  381.      9                INPLAS .EQ. 86 .OR. INPLAS .EQ. 70  ) THEN
  382. *
  383.                   IF (INPLAS.EQ.44.AND.VAR0(NVARI).EQ.0.0) THEN
  384.                         VAR0(NVARI)=XMAT(20)
  385.                   ENDIF
  386.                   IF (INPLAS.EQ.45.AND.VAR0(NVARI).EQ.0.0) THEN
  387.                         VAR0(NVARI-2)=XMAT(20)
  388.                         VAR0(NVARI-1)=XMAT(21)
  389.                         VAR0(NVARI)=XMAT(27)
  390.                   ENDIF
  391. *
  392.                   SEGINI WRK7
  393.                   SEGINI WRK9
  394.              CALL CONSTI(WRK0,WR00,WRK1,WRK5,WRK7,WRK8,WRK9,WTRAV,
  395.      1 INPLAS,MFR1,DT,NSTRSS,NVARI,NMATT,PRECIS,MSOUPA,JECHER,DTT,
  396.      2 NSSINC,INV,KERRE,ICARA,IFOURB,NYOG,NYNU,NYALFA,NYSMAX,NYN,
  397.      3 NYM,NYKK,NYALF1,NYBET1,NYR,NYA,NYKX,NNKX,NYRHO,NSIGY,TETA1,
  398.      5 TETA2,TREFA,TLIFE,ITHHER,NCOURB,CMATE,N2EL,N2PTEL,IB,IGAU,
  399.      6 EPAIST,NBPGAU,MELE,NPINT,NBGMAT,NELMAT,SECT,LHOOK,CRIGI,KERREU1)
  400. c
  401. c
  402. c          write(6,*) istep
  403.           IF (MFR1.EQ.17.AND.INPLAS.EQ.19) THEN
  404.           IF (KERREU1.NE.0.AND.NSSINC.EQ.1) THEN
  405.           CALL ERREUR(KERREU1)
  406.           ENDIF
  407.           ENDIF
  408.                   SEGSUP WRK7
  409.                   SEGSUP WRK9
  410.                   DTOPTI  = MIN(DTOPTI,DTT)
  411.                   NINCMA = MAX(NINCMA,NSSINC)
  412.                   NCOMP = NCOMP + 1
  413.                   TSOM = TSOM + DTT
  414.                   NSOM = NSOM + NSSINC
  415.                   NINV = NINV + INV
  416.                   TCAR = TCAR + DTT* DTT
  417.                   IF(KERRE.NE.0.AND.KERRE.NE.99) THEN
  418.                      KERR1=1
  419.                   END IF
  420. c
  421.             ELSE
  422.                KERRE = 99
  423.             ENDIF
  424. *
  425. *   Erreurs
  426. *      - problèmes de convergence
  427. *
  428.       CALL DEFER1(JNOID,KERR1,KERRE,LOGSUC)
  429. *
  430. *      - autres problèmes
  431. *
  432.  1990 CALL DEFER2(INPLAS,MFR,MELE,IB,IGAU,
  433.      .           KERR1,KERRE)
  434.  1998 IF (KERRE.NE.0) THEN
  435.          IF (LOGVIS) SEGSUP WRK8
  436.          SEGSUP WRK0,WR00,WRK1,WRK2,WRK3,WRK5,WTRAV
  437.          IF(MFR.EQ.7.OR.MFR.EQ.13.OR.LUNI1) THEN
  438.             SEGSUP WRK4
  439.          ENDIF
  440.          IF ((CMATE.EQ.'ORTHOTRO'.OR.CMATE.EQ.'ANISOTRO'.OR.
  441.      1     CMATE.EQ.'UNIDIREC').AND.(MFR.EQ.1.OR.MFR.EQ.31
  442.      1     .OR.MFR.EQ.33)) THEN
  443.            SEGDES MINTE2
  444.            SEGSUP WRK22
  445.          ENDIF
  446.          IF (ITHHER.EQ.1.OR.ITHHER.EQ.2) THEN
  447.             SEGDES MELVA3
  448.             SEGDES MELVA4
  449.             SEGDES MELVA5
  450.             SEGDES MCHAM3
  451.             SEGDES MCHAM4
  452.             SEGDES MCHAM5
  453.          ENDIF
  454.          RETURN
  455.       ENDIF
  456. c
  457. c
  458. c     remplissage du segment contenant les contraintes a la fin
  459. *     ( rearrangement pour milieu poreux ),
  460. c     les variables internes finales
  461. c     et les increments de deformations plastiques
  462. c
  463.       CALL DEFSIG(MFR,NDEF,
  464.      .     INPLAS,IND,WRK1,WRK5,WTRAV,
  465.      .     IVASTF,IVARIF,IVADEP,COB,XMOB,IB,IGAU,
  466.      .     CMATE,MATE,MELE,KERRER)
  467.       IF (KERRER.NE.0) GOTO 1000
  468. c
  469. c    fin de la boucle sur les points de gauss
  470. c
  471. 1100     continue
  472. c
  473. c    special poutres et tuyaux  sauf timoschenko
  474. c
  475.       CALL DEFPOU(MFR,MELE,MELEME,IB,WRK4,IVASTF)
  476. c
  477. c     fin de la boucle sur les elements
  478. c
  479. 1000  continue
  480. c
  481. *  FIN: modèles visqueux, on stocke le pas de temps
  482. *  optimal en indice 'dtopti'
  483. *
  484.       CALL DEFFIN(INPLAS,TSOM,NSOM,NCOMP,NINV,NINCMA,
  485.      .     TCAR,DTOPTI,IPOTAB,KERRE)
  486.       IF (LOGVIS) SEGSUP WRK8
  487. *
  488. *
  489.       SEGSUP WRK0,WR00,WRK1,WRK2,WRK3,WRK5,WTRAV
  490.       IF(MFR.EQ.7.OR.MFR.EQ.13.OR.LUNI1) THEN
  491.          SEGSUP WRK4
  492.       END IF
  493.       IF ((CMATE.EQ.'ORTHOTRO'.OR.CMATE.EQ.'ANISOTRO'.OR.
  494.      1     CMATE.EQ.'UNIDIREC').AND.(MFR.EQ.1.OR.MFR.EQ.31
  495.      1     .OR.MFR.EQ.33)) THEN
  496.         SEGDES MINTE2
  497.         SEGSUP WRK22
  498.       ENDIF
  499. *
  500.       IF (ITHHER.EQ.1.OR.ITHHER.EQ.2) THEN
  501.             SEGDES MELVA3
  502.             SEGDES MELVA4
  503.             SEGDES MELVA5
  504.             SEGDES MCHAM3
  505.             SEGDES MCHAM4
  506.             SEGDES MCHAM5
  507.       ENDIF
  508. *
  509.       RETURN
  510.       END
  511.  
  512.  
  513.  
  514.  
  515.  
  516.  
  517.  
  518.  
  519.  
  520.  
  521.  
  522.  
  523.  
  524.  
  525.  
  526.  
  527.  
  528.  
  529.  
  530.  
  531.  
  532.  
  533.  
  534.  
  535.  
  536.  
  537.  
  538.  
  539.  
  540.  
  541.  
  542.  
  543.  
  544.  

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