Télécharger la plaquette Cast3M


fr en

Télécharger ubiqui.eso

Retour à la liste

Numérotation des lignes : 

ubiqui

  1. C UBIQUI    SOURCE    BP208322  17/10/03    21:16:56     9580           
  2.       SUBROUTINE UBIQUI(DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL,
  3.      1 IB,IGAU,EPAIST,MELE,NPINT,SECT,LHOOK,TXR,XLOC,
  4.      2 XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,SIG0,DEPST,VAR0,XMAT,
  5.      3           NBPGAU,NCOMAT,XCAR,DSIGT,SIGF,VARF,DEFP,KERRE
  6.      4, ecou,necou,iecou)
  7. C
  8. C----------------------------------------------------------------------
  9. C             PLASTICITE MODELE UBIQUITOUS
  10. C
  11. C ENTREES
  12. C     SIG0(NSTRS) = CONTRAINTES INITIALES
  13. C     NSTRS       = NOMBRE DE CONTRAINTES
  14. C     DEPST(NSTRS)= INCREMENT DE DEFORMATIONS TOTALES
  15. C     VAR0(NVARI) =  VARIABLES INTERNES   DEBUT
  16. C         (  1  ) = EPSE
  17. C         (  2  ) = EPN1
  18. C         (  3  ) = GAP1
  19. C         (  4  ) = LAM1
  20. C         (  5  ) = EPN2
  21. C         (  6  ) = GAP2
  22. C         (  7  ) = LAM2
  23. C     XMAT(NCOMAT)=  COMPOSANTES DE MATERIAU
  24. C     NCOMAT      = NOMBRE DE COMPOSANTES DE MATERIAU
  25. C     XCAR(ICARA) =  CARACTERISTIQUES
  26. *     CMATE = NOM DU MATERIAU
  27. *     VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU
  28. *     VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES
  29. *     N2EL  = NBRE D ELEMENTS DANS SEGMENT DE HOOKE
  30. *     N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE
  31. *     MFR   = NUMERO DE LA FORMULATION
  32. *     IFOU  = OPTION DE CALCUL
  33. *     IB    = NUMERO DE L ELEMENT COURANT
  34. *     IGAU  = NUMERO DU POINT COURANT
  35. *     EPAIST= EPAISSEUR
  36. *     NBPGAU= NBRE DE POINTS DE GAUSS
  37. *     MELE  = NUMERO DE L ELEMENT FINI
  38. *     NPINT = NBRE DE POINTS D INTEGRATION
  39. *     NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES
  40. *     NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES
  41. *     SECT  = SECTION
  42. *     LHOOK = TAILLE DE LA MATRICE DE HOOKE
  43. *     TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES
  44. *     UTILISES POUR LE CALCUL DE LA MATRICE DE HOOKE
  45. *
  46. C  SORTIES
  47. C     SIGF(NSTRS)   = CONTRAINTES FINALES
  48. C     VARF(NVARI)   = VARIABLES INTERNES FINALES
  49. C     DEFP(NSTRS)   = DEFORMATIONS PLASTIQUES
  50. C     KERRE   = 0    TOUT OK
  51. C        1   SI DLAMBDA NEGATIF
  52. C        2   NOMBRE MAX D ITERATIONS INTERNES DEPASSE
  53. C       21   ON NE TROUVE PAS L INTERSECTION AVEC LA SURFACE DE CHARGE
  54. C       22   SIG0 A L EXTERIEUR DE LA SURFACE DE CHARGE
  55. C
  56. C-----------------------------------------------------------------------
  57. C      VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU  ET NECOU
  58. C
  59. C  IFOUR     INDICE DU TYPE DE PROBLEME
  60. C            -2  CONTRAINTES PLANES
  61. C            -1  DEFORMATIONS PLANES
  62. C             0  AXISYMETRIQUE
  63. C             1  SERIE DE FOURIER
  64. C             2  TRIDIMENSIONNEL
  65. C  ITYP      TYPE DE FORMULATION MECANIQUE
  66. C       ITYP=1   CAS DES ELEMENTS MASSIFS
  67. C       ITYP=2   CAS DES COQUES
  68. C       ITYP=3   CAS DES MEMBRANES
  69. C       ITYP=4   CAS DES CABLES ET DES BARRES
  70. C       ITYP=5   CAS QUELCONQUE
  71. C       ITYP=6   CAS DES CONTRAINTES PLANES
  72. C       ITYP=7   CAS DES COQUES A NU=0. OU CONTRAINTES PLANES
  73. C       ITYP=8   CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES
  74. C       ITYP=9   CAS DES COQUES EPAISSES
  75. C       ITYP=10  CAS DES JOINTS
  76. C       ITYP=11  CAS DES POUTRES
  77. C       ITYP=12  CAS DES TUYAUX
  78. C       ITYP=13  CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE
  79. C
  80. C-----------------------------------------------------------------------
  81. C
  82.       IMPLICIT INTEGER(I-N)
  83.       IMPLICIT REAL*8(A-H,O-Z)
  84.  
  85. -INC PPARAM
  86. -INC CCOPTIO
  87.       SEGMENT IECOU
  88. *      COMMON/IECOU/NYOG,NYNU,NYALFA,NYSMAX,NYN,NYM,NYKK,
  89.        INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  90. C      INTEGER NYOG, NYNU, NYALFA,NYSMAX,NYN, NYM,  NYKK,
  91.      1  icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  92. C    .  NYALF1,NYBET1,NYR,  NYA,   NYRHO,NSIGY,  NNKX,  NYKX,  IND,
  93.      2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  94. C    .  NSOM,  NINV, NINCMA,NCOMP, JELEM, LEGAUS,INAT,  NCXMAT,
  95.      3 icow25,icow26,icow27,icow28,icow29,icow30,icow31,
  96. C    .  LTRAC, MFR,  IELE,  NHRM,   NBNN, NBELEM,ICARA,
  97.      4  icow32,icow33,NSTRS,MFR,NBGMAT,NELMAT,icow38,
  98. C    .  LW2,   NDEF,  NSTRSS,MFR1,  NBGMAT,NELMAT,MSOUPA,
  99.      5  icow39,icow40,icow41,icow42,icow43,icow44
  100. C    .  NUMAT1,LENDO, NBBB,  NNVARI,KERR1, MELEME
  101.         INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  102.      .          icow51,icow52,icow53,icow54,icow55,icow56
  103.      .          icow57,icow58
  104.       ENDSEGMENT
  105.         SEGMENT ECOU
  106. ***     COMMON/ECOU/TEST,ALFAH,
  107. C       REAL*8 ecow00,ecow0,
  108.        REAL*8 TEST, ALFAH,
  109.      1 ecow1,ecow2,ecow3(6),ecow4(9),ecow5(6),
  110. C     1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6),
  111.      2 ecow6(12),ecow7(6),ecow8(6),ecow9(6),ecow10(6),ecow11(6),
  112.      2 ecow12(6),
  113. C     2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),W1(6),W2(6),
  114.      1 ecow13(6),ecow14(6),ecow15(12),ecow16(3),
  115. C     1 DALPHA(6),EPSPLA(6),E(12),XINV(3),
  116.      2 ecow17(6),ecow18(6),ecow19,ecow20
  117. C     2 SIPLAD(6),DSIGP0(6),TET,TETI
  118.       ENDSEGMENT
  119. C      COMMON/ECOU/TEST,ALFAH,ZOZO(112)
  120.       SEGMENT NECOU
  121. *      COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO,
  122. C       INTEGER ncow1,  ncow2,ncow3,ncow4,  ncow5,
  123.       INTEGER NCOURB,IPLAST,IT,   IMAPLA,ISOTRO,
  124.      1 ITYP, IFOURB, ncow8,
  125. C    . ITYP, IFOURB,IFLUAG,
  126.      2 ncow9,ncow10,ncow11,ncow12,ncow13,
  127. C    . ICINE,ITHER, IFLUPL,ICYCL, IBI,
  128.      3 ncow14,ncow15,ncow16,
  129. C    . JFLUAG,KFLUAG,LFLUAG,
  130.      4 ncow17,ncow18,ncow19,ncow20,ncow21
  131. C    . IRELAX,JNTRIN,MFLUAG,JSOUFL,JGRDEF
  132.       ENDSEGMENT
  133. C      COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO,ITYP,IFOURB,IZOZO(14)
  134. C
  135.       DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*)
  136.       DIMENSION SIGF(*),VARF(*),DEFP(*)
  137.       DIMENSION SIGMA(6),DSIGMA(6),STOT(6),SIGFIN(6)
  138.       DIMENSION XLAMBD(6),DEFPLA(6),DSIGT(*)
  139.       DIMENSION SPHER(6),AUXIL(6),DSIGZE(6)
  140.       DIMENSION VALMAT(*),VALCAR(*)
  141.       DIMENSION TXR(IDIM,*),CRIGI(12)
  142.       DIMENSION DDAUX(LHOOK,*),DDHOMU(LHOOK,*)
  143.       DIMENSION XLOC(3,3),XGLOB(3,3)
  144.       DIMENSION D1HOOK(LHOOK,*),ROTHOO(LHOOK,*)
  145.       CHARACTER*8 CMATE
  146. C
  147. C TEST DE CONSISTANCE DES DONNEES
  148. C
  149.       DO 100 I=1,NCOMAT
  150.        IF(XMAT(I).LT.0.D0) KERRE=538
  151.   100 CONTINUE
  152.       IF(KERRE.EQ.538) RETURN
  153. C
  154.       TRA1=0.D0
  155.       PHI1=0.D0
  156.       PSI1=0.D0
  157.       COHE1=0.D0
  158.       TPHI1=0.D0
  159.       TPSI1=0.D0
  160.       HACHE1=0.D0
  161.       TRA2=0.D0
  162.       PHI2=0.D0
  163.       PSI2=0.D0
  164.       COHE2=0.D0
  165.       TPHI2=0.D0
  166.       TPSI2=0.D0
  167.       HACHE2=0.D0
  168.       CO21=0.D0
  169.       CO22=0.D0
  170.       SI21=0.D0
  171.       SI22=0.D0
  172. C
  173.       YOUN=XMAT(1)
  174.       XNU =XMAT(2)
  175.       NCRI= INT(XMAT(3))
  176.       ANG1=XMAT(4)
  177.       TRA1=XMAT(5)
  178.       PHI1=XMAT(6)
  179.       PSI1=XMAT(7)
  180.       RHO =XMAT(8)
  181.       ALPH=XMAT(9)
  182.       ANG2=XMAT(10)
  183.       TRA2=XMAT(11)
  184.       PHI2=XMAT(12)
  185.       PSI2=XMAT(13)
  186. C
  187.       IF(IIMPI.EQ.28) WRITE(IOIMP,3000) (XMAT(I),I=1,13)
  188. C
  189. C  PETIT TEST SUR NU POUR CERTAINS CAS
  190. C
  191.       IF(MFR.EQ.3.AND.IFOURB.EQ.-2.AND.XNU.NE.0.D0) THEN
  192.         KERRE=38
  193.         RETURN
  194.       ENDIF
  195. C
  196. C
  197. C INITIALISATION
  198. C
  199.       NCAS=2
  200.       UNIT=0.0174532925199432957692D0
  201. C
  202. C
  203. C CALCUL ET TEST DES PARAMETRES DU MATERIAU UBIQUITOUS
  204. C
  205. C       TEST D'ERREUR SUR LE NOMBRE DE CRITERE
  206. C-----------------------------------------------
  207. C -----------------------------------------------------------
  208. C ON SE LIMITE AUX DEFORMATIONS PLANES , CONTRAINTES PLANES
  209. C ET AXISYMETRIQUE
  210. C------------------------------------------------------------
  211.       IF(IFOURB.LE.0)  GOTO 101
  212. C
  213.       KERRE=539
  214.       RETURN
  215. C-------------------------------------------------------------
  216.   101 IF(IFOURB.EQ.-2) NCAS=1
  217.       IF(NCRI.EQ.1.OR.NCRI.EQ.2) GO TO 102
  218. C
  219.       KERRE=540
  220.       RETURN
  221. C----------------------------------------------------------------
  222. C      TEST SUR L'EGALITE DE ANG1 ET ANG2 (CAS DE DEUX CRITERES)
  223. C----------------------------------------------------------------
  224.   102 IF (NCRI.EQ.1) GO TO 104
  225.       ANG=MIN (ABS(ANG1-ANG2),ABS(ABS(ANG2-ANG1)-180.D0))
  226.       IF (ANG.GT.1D-3) GO TO 103
  227. C
  228.       KERRE=541
  229.       RETURN
  230. C------------------------------------
  231. C        TEST SUR LA VALEUR DE ANG2
  232. C------------------------------------
  233.   103 IF(ANG2.GE.0.D0.AND.ANG2.LE.180.D0) GO TO 104
  234. C
  235.       KERRE=542
  236.       RETURN
  237. C
  238. C------ CAS D'UN SEUL CRITERE
  239. C
  240. C------------------------------------
  241. C        TEST SUR LA VALEUR DE ANG1
  242. C------------------------------------
  243.   104 IF(ANG1.GE.0.D0.AND.ANG1.LE.180.D0) GO TO 105
  244. C
  245.       KERRE=543
  246.       RETURN
  247. C===============================================
  248. C   CALCUL DES PARAMETRES GEOMETRIQUES
  249. C===============================================
  250.   105 IF (NCRI.EQ.1) GO TO 106
  251. C
  252. C      CAS DE DEUX CRITERES
  253. C
  254.       ANG2=ANG2*UNIT
  255.       CO22=COS(ANG2)
  256.       SI22=SIN(ANG2)
  257.       ANG1=ANG1*UNIT
  258.       CO11=COS(ANG1)
  259.       SI11=SIN(ANG1)
  260.       CO21=COS(ANG2-ANG1)
  261.       SI21=SIN(ANG2-ANG1)
  262.       GO TO 300
  263. C
  264. C      CAS D'UN SEUL CRITERE
  265. C
  266.   106 ANG1=ANG1*UNIT
  267.       CO11=COS(ANG1)
  268.       SI11=SIN(ANG1)
  269. C
  270.   300 CONTINUE
  271. C***************************************************************
  272.       IF (NCRI.EQ.1) GO TO 303
  273.       TRA2=XMAT(11)
  274.       PHI2=XMAT(12)
  275.       PSI2=XMAT(13)
  276. C---------------------------------------------------------------
  277. C        TEST SUR PHI2 (PHI2 COMPRIS ENTRE 0 ET 90 STRICTEMENT)
  278. C---------------------------------------------------------------
  279.       IF(PHI2.GE.0.D0.AND.PHI2.LT.90.D0) GO TO 301
  280.       KERRE=544
  281.       RETURN
  282. C---------------------------------------------------------------
  283. C        TEST SUR PSI2 (PSI2 COMPRIS ENTRE 0 ET 90 STRICTEMENT)
  284. C---------------------------------------------------------------
  285.   301 IF(PSI2.GE.0.D0.AND.PSI2.LT.90.D0) GO TO 302
  286.       KERRE=545
  287.       RETURN
  288. C--------------------------------------------
  289. C        TEST SUR TRA2 (TRA2 POSITIF OU NUL)
  290. C--------------------------------------------
  291.   302 IF(TRA2.GE.0.D0) GO TO 303
  292.       KERRE=546
  293.       RETURN
  294. C
  295. C---------------------------------------------------------------
  296. C        TEST SUR PHI1 (PHI1 COMPRIS ENTRE 0 ET 90 STRICTEMENT)
  297. C---------------------------------------------------------------
  298.   303 TRA1=XMAT(5)
  299.       PHI1=XMAT(6)
  300.       PSI1=XMAT(7)
  301.  
  302.       IF(PHI1.GE.0.D0.AND.PHI1.LT.90.D0) GO TO 304
  303.       KERRE=547
  304.       RETURN
  305. C---------------------------------------------------------------
  306. C        TEST SUR PSI1 (PSI1 COMPRIS ENTRE 0 ET 90 STRICTEMENT)
  307. C---------------------------------------------------------------
  308.   304 IF(PSI1.GE.0.D0.AND.PSI1.LT.90.D0) GO TO 305
  309.       KERRE=548
  310.       RETURN
  311. C-----------------------------------------------------
  312. C        TEST SUR TRA1 (TRA1 DOIT ETRE POSTIF OU NUL)
  313. C-----------------------------------------------------
  314.   305 IF(TRA1.GE.0.D0) GO TO 306
  315.       KERRE=549
  316.       MOTERR(1:4)='TRA1'
  317.       RETURN
  318. C===========================================================
  319. C   CALCUL DES PARAMETRES PHISIQUES
  320. C===========================================================
  321.   306 IF(NCRI.EQ.1) GO TO 307
  322. C---------------------------
  323. C   CALCUL DES PARAMETRES  2
  324. C---------------------------
  325.       TPSI2=TAN(PSI2*UNIT)
  326.       TPHI2=TAN(PHI2*UNIT)
  327.       COHE2=TPHI2*TRA2
  328. C------------------------------------------------------------
  329. C        CALCUL DE HACHE2(CAS=1 CONT PLANE;CAS=2 DEFO PLANE)
  330. C------------------------------------------------------------
  331.       IF (NCAS.EQ.1) HACHE2=TPSI2*TPHI2*YOUN/(1.D0-XNU*
  332.      1 XNU)+YOUN/(XNU+1.D0)
  333.       IF (NCAS.EQ.2) HACHE2=TPSI2*TPHI2*YOUN*(1.D0-XNU)/(1.D0+XNU)/
  334.      1 (1.D0-2.D0*XNU)+YOUN/(XNU+1.D0)
  335. C---------------------------------------
  336. C------------------------------
  337. C      CAS D'UN SEUL CRITERE
  338. C------------------------------
  339. C---------------------------
  340. C   CALCUL DES PARAMETRES  1
  341. C---------------------------
  342.   307 TPSI1=TAN(PSI1*UNIT)
  343.       TPHI1=TAN(PHI1*UNIT)
  344.       COHE1=TPHI1*TRA1
  345. C------------------------------------------------------------
  346. C        CALCUL DE HACHE1(CAS=1 CONT PLANE;CAS=2 DEFO PLANE)
  347. C------------------------------------------------------------
  348.       IF (NCAS.EQ.1) HACHE1=TPSI1*TPHI1*YOUN/(1.D0-XNU*
  349.      1 XNU)+YOUN/(XNU+1.D0)
  350.       IF (NCAS.EQ.2) HACHE1=TPSI1*TPHI1*YOUN*(1.D0-XNU)/(1.D0+XNU)/
  351.      1 (1.D0-2.D0*XNU)+YOUN/(XNU+1.D0)
  352. C
  353. C IMPRESSION DES VARIABLES
  354. C
  355.       IF(IIMPI.EQ.28) THEN
  356.         WRITE(IOIMP,1000)
  357.         WRITE(IOIMP,1001) NCRI,NCAS
  358.         WRITE(IOIMP,1002) ANG1,TRA1,COHE1
  359.         WRITE(IOIMP,1003) PHI1,PSI1,HACHE1
  360.         WRITE(IOIMP,1004) ANG2,TRA2,COHE2
  361.         WRITE(IOIMP,1005) PHI2,PSI2,HACHE2
  362.       ENDIF
  363. C
  364. C CALCUL DE L INCREMENT DE CONTRAINTE
  365. C
  366.       ifou=ifourb
  367.       nstrbi=nstrs
  368.       nbgmab=nbgmat
  369.       nlmatb=nelmat
  370.       mfrbi=mfr
  371.       CALL CALSIG(DEPST,DDAUX,NSTRbi,CMATE,VALMAT,VALCAR,
  372.      1    N2EL,N2PTEL,MFR,IFOU,IB,IGAU,EPAIST,NBPGAU,
  373.      2       MELE,NPINT,NBGMAb,NEMATb,SECT,LHOOK,TXR,XLOC,
  374.      3        XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,IRTD)
  375.  
  376. *
  377.       mfr=mfrbi
  378.       ifourb=ifou
  379.       nstrs=nstrbi
  380.       nbgmat=nbgmab
  381.       nelmat=nlmatb
  382.       IF(IRTD.NE.1) THEN
  383.          KERRE=69
  384.          GOTO 1900
  385.       ENDIF
  386. *
  387. C
  388. C ECOULEMENT PLASTIQUE
  389. C
  390.       MCOD=1
  391.       nstrbi=nstrs
  392.       mfrbi=mfr
  393.       CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  394.      .      SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  395.      .      DSIGZE,ICENT2,MCOD,IBOU,MFRbi,NSTRbi,XCAR,CMATE,ecou,necou)
  396.       mfr=mfrbi
  397.       nstrs=nstrbi
  398.       IF(ITYP.EQ.0) THEN
  399.         KERRE=269
  400.         RETURN
  401.       ENDIF
  402. C
  403. C CAS DES CONTRAINTES PLANES :
  404. C             IFOURB = -2
  405. C AVEC DEUX TYPES DE FORMULATION MECANIQUE :
  406. C          CAS DES COQUES :             ITYP = 2 ==> ALFAH = 0
  407. C                                   OU  ITYP = 7 ==> ALFAH = 0
  408. C          CAS DES CONTRAINTES PLANES : ITYP = 6
  409. C
  410. C DANS LE CAS DES COQUES ITYP = 2  , NOUS METTONS ALFAH = 0
  411. C
  412.       IF(ITYP.EQ.2.OR.ITYP.EQ.7) ALFAH=0.D0
  413. C
  414. C INITIALISATIONS DES PARAMETRES
  415. C
  416.       EPSE     =VAR0(1)
  417.       XLAMBD(1)=VAR0(2)
  418.       XLAMBD(2)=VAR0(3)
  419.       XLAMBD(5)=VAR0(4)
  420.       XLAMBD(3)=VAR0(5)
  421.       XLAMBD(4)=VAR0(6)
  422.       XLAMBD(6)=VAR0(7)
  423. C
  424.       DO 200 IBA=1,IBOU
  425.        STOT(IBA)=SIGMA(IBA)+DSIGMA(IBA)
  426.   200 CONTINUE
  427.       IF(IIMPI.EQ.28) THEN
  428.         WRITE(IOIMP,2000)
  429.         WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  430.         WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  431.         WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  432.       ENDIF
  433. C
  434. C ECOULEMENT PLASTIQUE
  435. C
  436.       CALL ECUBI(SIGMA,DSIGMA,XLAMBD,DEFPLA,YOUN,XNU,
  437.      .           NCRI,ANG1,TRA1,TPHI1,TPSI1,HACHE1,COHE1,CO11,SI11,
  438.      .           NCAS,ANG2,TRA2,TPHI2,TPSI2,HACHE2,COHE2,CO22,SI22,
  439.      .           CO21,SI21,
  440.      .           SIGFIN,IDAM,KERRE)
  441. C
  442.       IF(IIMPI.EQ.28) THEN
  443.         WRITE(IOIMP,2000)
  444.         WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  445.         WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  446.         WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  447.         WRITE(IOIMP,2004) (SIGFIN(I),I=1,IBOU)
  448.         WRITE(IOIMP,2005) (DEFPLA(I),I=1,IBOU)
  449.         WRITE(IOIMP,2006) (XLAMBD(I),I=1,6)
  450.       ENDIF
  451. C
  452.       IF(KERRE.NE.0) RETURN
  453. C
  454. C PAS DE PLASTICITE COMPORTEMENT ELASTIQUE
  455. C
  456.       IF(IDAM.EQ.0) THEN
  457.         DO 400 IBA=1,IBOU
  458.          SIGFIN(IBA)=STOT(IBA)
  459.          DEFPLA(IBA)=0.D0
  460.   400   CONTINUE
  461.       ENDIF
  462. C
  463. C NOUS AVONS ENDOMMAGE LE MATERIAU
  464. C
  465. C MISE A JOUR DES VARIABLES INTERNES
  466. C
  467.       VARF(1)=EPSE+(SQRT(XLAMBD(1)*XLAMBD(1)+XLAMBD(2)*XLAMBD(2)+
  468.      .                   XLAMBD(3)*XLAMBD(3)+XLAMBD(4)*XLAMBD(4)))
  469.       VARF(2)=XLAMBD(1)
  470.       VARF(3)=XLAMBD(2)
  471.       VARF(4)=XLAMBD(5)
  472.       VARF(5)=XLAMBD(3)
  473.       VARF(6)=XLAMBD(4)
  474.       VARF(7)=XLAMBD(6)
  475. C
  476.       MCOD=2
  477.       nstrbi=nstrs
  478.       mfrbi=mfr
  479.       CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  480.      .        SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  481.      .    DSIGZE,ICENT2,MCOD,IBOU,MFRbi,NSTRbi,XCAR,CMATE,ecou,necou)
  482.       mfr=mfrbi
  483.       nstrs=nstrbi
  484. C
  485. C LES FORMATS D IMPRESSION
  486. C
  487.  1000 FORMAT(1X,'CONSTANTES DU MATERIAU UBIQUITOUS')
  488.  1001 FORMAT(1X,'NCRI =',I4,3X,'NCAS =',I4)
  489.  1002 FORMAT(1X,'ANG1 =',E12.5,2X,'TRA1 =',E12.5,2X,'COHE1 =',E12.5)
  490.  1003 FORMAT(1X,'PHI1 =',E12.5,2X,'PSI1 =',E12.5,2X,'HACH1 =',E12.5)
  491.  1004 FORMAT(1X,'ANG2 =',E12.5,2X,'TRA2 =',E12.5,2X,'COHE2 =',E12.5)
  492.  1005 FORMAT(1X,'PHI1 =',E12.5,2X,'PSI1 =',E12.5,2X,'HACH2 =',E12.5)
  493.  2000 FORMAT(1X,'RESULTATS DE L ECOULEMENT PLASTIQUE')
  494.  2001 FORMAT(1X,'SIGMA  =',6(1X,1PE12.5))
  495.  2002 FORMAT(1X,'DSIGMA =',6(1X,1PE12.5))
  496.  2003 FORMAT(1X,'STOT   =',6(1X,1PE12.5))
  497.  2004 FORMAT(1X,'SIGFIN =',6(1X,1PE12.5))
  498.  2005 FORMAT(1X,'DEFPLA =',6(1X,1PE12.5))
  499.  2006 FORMAT(1X,'XLAMBD =',6(1X,1PE12.5))
  500.  3000 FORMAT(1X,'XMAT   =',/,3(6(1X,1PE12.5),/))
  501. C
  502.  1900 RETURN
  503.       END
  504.  
  505.  
  506.  
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  

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