Télécharger la plaquette Cast3M


fr en

Télécharger c2dp.eso

Retour à la liste

Numérotation des lignes : 

c2dp

  1. C C2DP      SOURCE    CHAT      05/01/12    21:44:59     5004
  2.       SUBROUTINE C2DP(SIG,DSIG,YOUN,ANU,RT1,RT2,RDP,ADP,HDP,
  3.      1 XLAM1,XLAM2,XLAM3,IDAM,ANG,EPPLDP,ALPHA,KERRE)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9.       DIMENSION SIG(3),DSIG(3),SGG(3),DSGG(3),SFG(3),DSFG(3),
  10.      1 A(3),B(3),IDAM(3),EPC(3),SSI(3),SI(3),EPC1(3),
  11.      1 EPDP(3),EC(3),EPPLDP(3),SG(3),AT(4),BT(2),TENS(3)
  12.       ITER=0
  13.       Y=YOUN/(1.D0-ANU*ANU)
  14.       SIREF=1.E-6*YOUN
  15.       DO 10 ITYP=1,3
  16.       EPDP(ITYP)=0.D0
  17.    10 EPPLDP(ITYP)=0.D0
  18.       DL3=0.D0
  19.       XT=0.D0
  20.       IF(IIMPI.EQ.9) WRITE(IOIMP,9999)
  21. 9999  FORMAT(1X,'C2DP COUPLAGE COMP 2 DRUCKER',/)
  22.  9998 FORMAT(1X,'C2DP ',I4,'ITERATIONS INTERNES',/)
  23. C
  24. C------------------------------------------------
  25. C     COUPLAGE COMPRESSION 2 DRUCKER
  26. C------------------------------------------------
  27. C----------------------------------------------
  28. C     DEFINITION D UNE CONTRAINTE DE REFERENCE
  29. C     SI RMAX DU DRUCKER INFERIEURE A CETTE VALEUR
  30. C     ALORS SIGMA=0
  31. C----------------------------------------------
  32. C
  33.       RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  34.       IF(SIREF.GT.RMAX) THEN
  35.       DO 20 ITYP=1,3
  36.   20  SIG(ITYP)=SIG(ITYP)+DSIG(ITYP)
  37.       RDP=0.D0
  38.       HDP=0.D0
  39.       ADP=0.D0
  40.       CALL CPHOIN(SIG,EPPLDP,YOUN,ANU)
  41.       CALL NORME(SIG,DL3)
  42.       DL3=DL3/YOUN*10.D0
  43.       XLAM3=XLAM3+DL3
  44.       DO 30 ITYP=1,3
  45.       IDAM(ITYP)=0
  46.       DSIG(ITYP)=0.D0
  47.   30  SIG(ITYP)=0.D0
  48.       IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  49.       RETURN
  50.       ENDIF
  51. C-----------------------------------------------
  52. C     ON SEPLACE DANS LE REPERE DE FISSURATION
  53. C-----------------------------------------------
  54. C
  55.       CALL CHREP(SIG,SFG,ANG)
  56.       CALL CHREP(DSIG,DSFG,ANG)
  57. C
  58. C-----------------------------------------------------
  59. C     ESTIMATION DU PAS D'INCREMENT DE CONTRAINTES
  60. C-----------------------------------------------------
  61. C
  62.       CALL SIJ(SFG,SI,SIETOI)
  63.       CALL EPSDP(SI,SIETOI,ADP,EPC)
  64.       CALL CPHOMO(EPC,SSI,YOUN,ANU,ALPHA)
  65.       CALL SCAL2(EPC,SSI,VAL)
  66.       VAL=VAL-HDP
  67.       IF(VAL.LT.0.D0) THEN
  68.       KERRE=459
  69.       RETURN
  70.       ENDIF
  71.       CALL SCAL(SSI,DSFG,VAL)
  72.       CALL NORME(DSFG,VA1)
  73.       CALL NORME(SSI,VA2)
  74.       IF(VA1.EQ.0.D0) THEN
  75.       IDAM(1)=0
  76.       IDAM(2)=0
  77.       IDAM(3)=0
  78.       RETURN
  79.       ENDIF
  80.       RMIN=MIN((RDP/1.73),(RDP/(1.D0+2.D0*ADP)))
  81.       X=VAL/VA1/VA2
  82.       IF(X.GE.1.D0) THEN
  83.       X=1.D0
  84.       GOTO 41
  85.       ENDIF
  86.       X=1.D0/SQRT(1.0001D0-X*X)*RMIN/VA1/8.D0
  87.       IF(X.GE.1.D0) X=1.D0
  88. C
  89. C----------------------------------------------------------
  90. C     ON ECOULE ET ON REGARDE LA VARIATION DE NORMALE
  91. C----------------------------------------------------------
  92. C
  93. C
  94. C---------------------------------------
  95. C     RESOLUTION DU SYSTEME EN DL3
  96. C---------------------------------------
  97.    41 X=X/2.D0
  98.       EC(1)=0.D0
  99.       EC(2)=1.D0
  100.       EC(3)=0.D0
  101.       CALL CPHOOB(EC,SG,YOUN,ANU)
  102.   100 X=X*2.D0
  103.       CALL SIJ(SFG,SI,SIETOI)
  104.       CALL EPSDP(SI,SIETOI,ADP,EPC)
  105.       CALL CPHOMO(EPC,SSI,YOUN,ANU,ALPHA)
  106.       CALL SCAL2(EPC,SSI,VAL)
  107.       VAL=VAL-HDP
  108.       IF(VAL.LT.0.D0) THEN
  109.       KERRE=459
  110.       RETURN
  111.       ENDIF
  112.       IF((XT+X).GT.1.D0) X=1.D0-XT
  113.   110 RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  114. C
  115. C --------------------------------------------------
  116. C     CAS OU LE RAYON MAXIMAL DU DRUCKER-PRAGER
  117. C     INFERIEUR A LA CONTRAINTE DE REFERENCE
  118. C --------------------------------------------------
  119. C
  120.       IF(SIREF.LT.RMAX) GOTO 131
  121.       DO 120 ITYP=1,3
  122.   120 SFG(ITYP)=SFG(ITYP)+(1.D0-XT)*DSFG(ITYP)
  123.       RDP=0.D0
  124.       ADP=0.D0
  125.       HDP=0.D0
  126.       CALL NORME(SFG,DL3)
  127.       DL3=DL3/YOUN*10.D0
  128.       XLAM3=XLAM3+DL3
  129.       CALL CPHOIN(SFG,EPPLDP,YOUN,ANU)
  130.       DO 130 ITYP=1,3
  131.       IDAM(ITYP)=0
  132.       EPDP(ITYP)=EPDP(ITYP)+EPPLDP(ITYP)
  133.       DSIG(ITYP)=0.D0
  134.   130 SIG(ITYP)=0.D0
  135.       CALL CHREP(EPDP,EPPLDP,-ANG)
  136.       IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  137.       RETURN
  138.   131 ITER=ITER+1
  139.       IF(ITER.GT.200) THEN
  140.       KERRE=460
  141.       RETURN
  142.       ENDIF
  143. C
  144. C---------------------------------
  145. C     ON REALISE LE COUPLAGE
  146. C---------------------------------
  147. C
  148.       E=-(SFG(2)+X*DSFG(2))/SG(2)
  149.       F=SSI(2)/SG(2)
  150.       DO 140 ITYP=1,3
  151.       DSGG(ITYP)=X*DSFG(ITYP)
  152.       A(ITYP)=SFG(ITYP)+DSGG(ITYP)+E*SG(ITYP)
  153.   140 B(ITYP)=-SSI(ITYP)+F*SG(ITYP)
  154.       CALL DLAMDP(A,B,DL3,RDP,ADP,HDP,ITEST)
  155.       IF(ITEST.EQ.1) THEN
  156.       X=X/2.D0
  157.       GOTO 131
  158.       ENDIF
  159.       DL2=-E-F*DL3
  160.       DO 145 ITYP=1,3
  161.   145 SGG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  162. C
  163. C----------------------------------------------
  164. C     SI UN DES DLAMDA < 0  FAUX COUPLAGE
  165. C----------------------------------------------
  166. C
  167.       IF(DL2.GT.0.D0) IDAM(2)=0
  168.       IF(DL3.LT.0.D0) IDAM(3)=0
  169.       IF(IDAM(2).EQ.0.OR.IDAM(3).EQ.0) THEN
  170.       DL2=0.D0
  171.       DL3=0.D0
  172.       X=0.D0
  173.       GOTO 3000
  174.       ENDIF
  175. C
  176. C-----------------------------------------------------
  177. C     ON VERIFIE QUE LA NORMALE SUIVANT LE DRUCKER
  178. C     PRAGER TOURNE PEU
  179. C-----------------------------------------------------
  180. C
  181.       CALL SIJ(SGG,SI,SIETOI)
  182.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  183.       CALL SCAL(EPC,EPC1,VAL)
  184.       CALL NORME(EPC,VA1)
  185.       CALL NORME(EPC1,VA2)
  186.       CO=VAL/VA1/VA2
  187.       IF(CO.LT.0.99) THEN
  188.       X=X/2.D0
  189.       GOTO 131
  190.       ENDIF
  191.       IF((DL2+XLAM2).GE.0.D0) GOTO 200
  192. C
  193. C------------------------------------------------
  194. C     CAS OU LA FISSURE EST COMPLETEMENT FERMEE
  195. C------------------------------------------------
  196. C
  197.       IDAM(2)=0
  198.       DL2=-XLAM2
  199.       IF(DSFG(2).EQ.0.D0) GOTO 160
  200.       E=(DL2*SG(2)-SFG(2))/DSFG(2)
  201.       F=SSI(2)/DSFG(2)
  202.       DO 150 ITYP=1,3
  203.       A(ITYP)=SFG(ITYP)-DL2*SG(ITYP)+E*DSFG(ITYP)
  204.   150 B(ITYP)=F*DSFG(ITYP)-SSI(ITYP)
  205.       CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  206.       DO 151 ITYP=1,3
  207.       SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  208.   151 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  209.       CALL SIJ(SGG,SI,SIETOI)
  210.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  211.       CALL SCAL(EPC,EPC1,VAL)
  212.       CALL NORME(EPC,VA1)
  213.       CALL NORME(EPC1,VA2)
  214.       CO1=VAL/VA1/VA2
  215.       CALL SIJ(DSGG,SI,SIETOI)
  216.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  217.       CALL SCAL(EPC,EPC1,VAL)
  218.       CALL NORME(EPC,VA1)
  219.       CALL NORME(EPC1,VA2)
  220.       CO2=VAL/VA1/VA2
  221.       XIN1=E+F*DL31
  222.       XIN2=E+F*DL32
  223.       IF(XIN1.GT.-1.D-10.AND.XIN1.LE.X.AND.CO1.GT.0.9) THEN
  224.       DL3=DL31
  225.       X=XIN1
  226.       GOTO 200
  227.       ENDIF
  228.       IF(XIN2.GT.-1.D-10.AND.XIN2.LE.X.AND.CO2.GT.0.9) THEN
  229.       DL3=DL32
  230.       X=XIN2
  231.       GOTO 200
  232.       ENDIF
  233.       IF(IIMPI.EQ.9) WRITE(IOIMP,10101)
  234. 10101 FORMAT(1X,'ERREUR DANS C2DP FERM ',/)
  235.       DL3=DL31
  236.       X=XIN1
  237.       GOTO 200
  238.   160 DL3=(SFG(2)-DL2*SG(2))/SSI(2)
  239.       RD=RDP-HDP*DL3
  240.       IF(RD.LT.0.D0) THEN
  241.       RDP=0.D0
  242.       GOTO 110
  243.       ENDIF
  244.       DO 170 ITYP=1,3
  245.       A(ITYP)=SFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  246.   170 B(ITYP)=DSFG(ITYP)
  247.       CALL XDP(A,B,RD,ADP,X,ITEST)
  248.   200 CONTINUE
  249.       DO 205 ITYP=1,3
  250.       DSGG(ITYP)=X*DSFG(ITYP)
  251.   205 SGG(ITYP)=SFG(ITYP)+DSGG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  252. C
  253. C-----------------------------------------------
  254. C     ON REGARDE SI ON ENDOMMAGE PAS LES AUTRES
  255. C     CRITERES
  256. C-----------------------------------------------
  257. C-------------------------------------
  258. C     ON ENDOMMAGE LA TRACTION 1
  259. C-------------------------------------
  260. C
  261.       CALL CTRAF(SGG(1),RT1,VCTR)
  262.       IF(VCTR.LE.0.D0) GOTO 1000
  263.       IDAM(1)=1
  264.       IDAM(2)=-1
  265.       IDAM(3)=1
  266.       AT(1)=-DSFG(1)
  267.       AT(4)=Y
  268.       AT(2)=Y*ANU
  269.       AT(3)=-DSFG(2)
  270.       BT(1)=SFG(1)-RT1
  271.       BT(2)=SFG(2)
  272.       CALL SYLIN2(AT,BT,X1,DL21)
  273.       BT(1)=-SSI(1)
  274.       BT(2)=-SSI(2)
  275.       CALL SYLIN2(AT,BT,X2,DL22)
  276.       DO 230 ITYP=1,3
  277.       A(ITYP)=SFG(ITYP)-DL21*SG(ITYP)+X1*DSFG(ITYP)
  278.   230 B(ITYP)=-DL22*SG(ITYP)+X2*DSFG(ITYP)-SSI(ITYP)
  279.       CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  280.       DO 231 ITYP=1,3
  281.       SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  282.   231 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  283.       CALL SIJ(SGG,SI,SIETOI)
  284.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  285.       CALL SCAL(EPC,EPC1,VAL)
  286.       CALL NORME(EPC,VA1)
  287.       CALL NORME(EPC1,VA2)
  288.       CO1=VAL/VA1/VA2
  289.       CALL SIJ(DSGG,SI,SIETOI)
  290.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  291.       CALL SCAL(EPC,EPC1,VAL)
  292.       CALL NORME(EPC,VA1)
  293.       CALL NORME(EPC1,VA2)
  294.       CO2=VAL/VA1/VA2
  295.       DLIN1=DL21+DL31*DL22
  296.       DLIN2=DL21+DL32*DL22
  297.       XIN1=X1+DL31*X2
  298.       XIN2=X1+DL32*X2
  299.       IF(XIN1.GT.-1.D-10.AND.DLIN1.LT.1.E-10.AND.XIN1.LE.X
  300.      1 .AND.CO1.GT.0.9) THEN
  301.       DL3=DL31
  302.       X=XIN1
  303.       DL2=DLIN1
  304.       GOTO 3000
  305.       ENDIF
  306.       IF(XIN2.GT.-1.D-10.AND.DLIN2.LT.1.E-10.AND.XIN2.LE.X
  307.      1 .AND.CO2.GT.0.9) THEN
  308.       DL3=DL32
  309.       X=XIN2
  310.       DL2=DLIN2
  311.       GOTO 3000
  312.       ENDIF
  313.       IF(IIMPI.EQ.9) WRITE(IOIMP,20202)
  314. 20202 FORMAT(1X,'ERREUR DANS C2DP TRAC ',/)
  315.       DL3=DL31
  316.       DL2=DLIN1
  317.       X=XIN1
  318.       GOTO 3000
  319. C
  320. C-------------------------------------
  321. C     ON ENDOMMAGE LA COMPRESSION 1
  322. C-------------------------------------
  323. C
  324.  1000 CALL CCOAF(SGG(1),XLAM1,VCCO1)
  325.       IF(VCCO1.LE.0.D0) GOTO 2000
  326.       IDAM(1)=-1
  327.       IDAM(3)=1
  328.       IDAM(2)=-1
  329.       AT(1)=-DSFG(1)
  330.       AT(4)=Y
  331.       AT(2)=Y*ANU
  332.       AT(3)=-DSFG(2)
  333.       BT(1)=SFG(1)
  334.       BT(2)=SFG(2)
  335.       CALL SYLIN2(AT,BT,X1,DL21)
  336.       BT(1)=-SSI(1)
  337.       BT(2)=-SSI(2)
  338.       CALL SYLIN2(AT,BT,X2,DL22)
  339.       DO 1010 ITYP=1,3
  340.       A(ITYP)=SFG(ITYP)-DL21*SG(ITYP)+X1*DSFG(ITYP)
  341.  1010 B(ITYP)=-DL22*SG(ITYP)+X2*DSFG(ITYP)-SSI(ITYP)
  342.       CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  343.       DO 1011 ITYP=1,3
  344.       SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  345.  1011 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  346.       CALL SIJ(SGG,SI,SIETOI)
  347.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  348.       CALL SCAL(EPC,EPC1,VAL)
  349.       CALL NORME(EPC,VA1)
  350.       CALL NORME(EPC1,VA2)
  351.       CO1=VAL/VA1/VA2
  352.       CALL SIJ(DSGG,SI,SIETOI)
  353.       CALL EPSDP(SI,SIETOI,ADP,EPC1)
  354.       CALL SCAL(EPC,EPC1,VAL)
  355.       CALL NORME(EPC,VA1)
  356.       CALL NORME(EPC1,VA2)
  357.       CO2=VAL/VA1/VA2
  358. *     DLIN1=DL11+DL31*DL12   CORRECTION PV ???
  359. *     DLIN2=DL11+DL32*DL12   SUITE ERREUR APPOLO
  360.       DLIN1=DL21+DL31*DL22
  361.       DLIN2=DL21+DL32*DL22
  362.       XIN1=X1+DL31*X2
  363.       XIN2=X1+DL32*X2
  364.       IF(XIN1.GT.-1.D-10.AND.DLIN1.LT.1.E-10.AND.XIN1.LE.X
  365.      1 .AND.CO1.GT.0.9) THEN
  366.       DL3=DL31
  367.       X=XIN1
  368.       DL2=DLIN1
  369.       GOTO 3000
  370.       ENDIF
  371.       IF(XIN2.GT.-1.D-10.AND.DLIN2.LT.1.E-10.AND.XIN2.LE.X
  372.      1 .AND.CO2.GT.0.9) THEN
  373.       DL3=DL32
  374.       X=XIN2
  375.       DL2=DLIN2
  376.       GOTO 3000
  377.       ENDIF
  378.       IF(IIMPI.EQ.9) WRITE(IOIMP,30303)
  379. 30303 FORMAT(1X,'ERREUR DANS C2DP COMP ',/)
  380.       DL3=DL31
  381.       DL2=DLIN1
  382.       X=XIN1
  383.       GOTO 3000
  384.  2000 CONTINUE
  385. C
  386. C------------------------------------------------
  387. C     CAS OU AUCUN AUTRE CRITERE EST ENDOMMAGE
  388. C------------------------------------------------
  389. C-----------------------------------
  390. C     REMISE A JOUR DES VARIABLES
  391. C     PUIS NOUVEAU PAS SI TOUT
  392. C     L INCREMENT N EST PAS ECOULE
  393. C-----------------------------------
  394. C
  395.       XT=XT+X
  396.       RDP=RDP-DL3*HDP
  397.       RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  398.       IF(SIREF.GT.RMAX) THEN
  399.       RDP=0.D0
  400.       GOTO 110
  401.       ENDIF
  402.       DO 2010 ITYP=1,3
  403.       SFG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  404.  2010 EPDP(ITYP)=DL3*EPC(ITYP)+EPDP(ITYP)
  405.       XLAM2=XLAM2+DL2
  406.       XLAM3=XLAM3+DL3
  407.       IF(XLAM2.LT.1.E-10) XLAM2=0.D0
  408.       IF(XLAM2.LT.1.E-10) IDAM(2)=0
  409.       C=1.D0-1.D-10
  410.       IF(XT.LT.C) GOTO 2020
  411.       DO 2015 ITYP=1,3
  412.       IDAM(ITYP)=0
  413.  2015 DSIG(ITYP)=0.D0
  414.       CALL CHREP(SFG,SIG,-ANG)
  415.       CALL CHREP(EPDP,EPPLDP,-ANG)
  416.       IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  417.       RETURN
  418.  2020 IF(IDAM(2).EQ.-1) GOTO 100
  419.       XLAM2=0.D0
  420.       GAMDP=10.D0
  421.       GAMTR1=10.D0
  422.       GAMTR2=10.D0
  423.       GAMCO1=10.D0
  424.       DO 2030 ITYP=1,3
  425.       DSFG(ITYP)=DSFG(ITYP)*(1.D0-XT)
  426.  2030 TENS(ITYP)=SFG(ITYP)+DSFG(ITYP)
  427.       CALL CTRAF(TENS(1),RT1,VCTR1)
  428.       CALL CTRAF(TENS(2),RT2,VCTR2)
  429.       CALL CCOAF(TENS(1),XLAM1,VCCO1)
  430.       CALL CDP(TENS,ADP,RDP,VCDP)
  431.       IF(VCTR1.GT.0.D0) CALL GAMTAF(SFG(1),DSFG(1),RT1,GAMTR1)
  432.       IF(VCTR2.GT.0.D0) CALL GAMTAF(SFG(2),DSFG(2),RT2,GAMTR2)
  433.       IF(VCCO1.GT.0.D0) CALL GAMCAF(SFG(1),DSFG(1),GAMCO1)
  434.       IF(VCDP.GT.0.D0) CALL GDP(SFG,DSFG,RDP,ADP,GAMDP)
  435.       IDAM(1)=0
  436.       IDAM(2)=0
  437.       IDAM(3)=0
  438.       GAM=MIN(GAMTR1,GAMTR2,GAMCO1,GAMDP)
  439.       IF(GAM.GE.1.D0) THEN
  440.       DO 2040 ITYP=1,3
  441.       SFG(ITYP)=SFG(ITYP)+DSFG(ITYP)
  442.  2040 DSIG(ITYP)=0.D0
  443.       CALL CHREP(SFG,SIG,-ANG)
  444.       RETURN
  445.       ENDIF
  446.       IF(ABS(GAM-GAMTR1).LE.1.E-10) IDAM(1)=1
  447.       IF(ABS(GAM-GAMTR2).LE.1.E-10) IDAM(2)=1
  448.       IF(ABS(GAM-GAMCO1).LE.1.E-10) IDAM(1)=-1
  449.       IF(ABS(GAM-GAMDP).LE.1.E-10) IDAM(3)=1
  450.       DO 2050 ITYP=1,3
  451.       SFG(ITYP)=SFG(ITYP)-GAM*DSFG(ITYP)
  452.  2050 DSFG(ITYP)=(1.D0-GAM)*DSFG(ITYP)
  453.       CALL CHREP(SFG,SIG,-ANG)
  454.       CALL CHREP(DSFG,DSIG,-ANG)
  455.       CALL CHREP(EPDP,EPPLDP,-ANG)
  456.       RETURN
  457.  3000 RDP=RDP-DL3*HDP
  458.       RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  459.       IF(SIREF.GT.RMAX) THEN
  460.       RDP=0.D0
  461.       GOTO 110
  462.       ENDIF
  463.       XT=XT+X
  464.       DO 3010 ITYP=1,3
  465.       SFG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  466.       DSFG(ITYP)=(1.D0-XT)*DSFG(ITYP)
  467.  3010 EPDP(ITYP)=DL3*EPC(ITYP)+EPDP(ITYP)
  468.       CALL CHREP(SFG,SIG,-ANG)
  469.       CALL CHREP(DSFG,DSIG,-ANG)
  470.       CALL CHREP(EPDP,EPPLDP,-ANG)
  471.       XLAM2=XLAM2+DL2
  472.       XLAM3=XLAM3+DL3
  473.       IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  474.       RETURN
  475.       END
  476.  
  477.  

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