Télécharger cendom.eso

Retour à la liste

Numérotation des lignes :

  1. C CENDOM SOURCE BP208322 17/03/01 21:15:26 9325
  2. SUBROUTINE CENDOM(wrk52,wrk53,wrk54,WRK6,WRK7,WRK8,WRK9,
  3. 1 NVARI, T0,T1,TREF,IB,IGAU,ifourb,ncourb,iecou)
  4.  
  5. C ENDOM SOURCE NICOLAS 99/02/10 21:16:58 3476
  6. * SUBROUTINE ENDOM(WRK0,WR00,WRK1,WRK6,WRK7,WRK8,WRK9,WTRAV,
  7. * 1 NSTRS,NXMATT,ICARA,INPLAS,NVARI,PRECIS,MFR,IFOURB,KERRE,
  8. * 2 NNVARI,NYOG,NYNU,NYRHO,NYALFA,NNKX,NYKX,NCOURB,NYN,NYM,
  9. * 3 NYKK,T0,T1,TREF,ITHHER,CMATE,N2EL,N2PTEL,IB,IGAU,
  10. * 4 EPAIST,MELE,NPINT,NBGMAT,NELMAT,SECT,LHOOK,CRIGI)
  11. C
  12. c NSTRSS <= NSTRS , MFR1 <- MFR , XCARB <- XCAR
  13. C ==================================================================
  14. C CE SOUS-PROGRAMME EST APPELE DANS "PLAST".
  15. C C'EST UN INTEGRATEUR GENERAL DES LOIS DE COMPORTEMENT DES
  16. C MATERIAUX ENDOMMAGEABLES EN REGIME ANISOTHERME.
  17. C
  18. C ENTREES:
  19. C -------
  20. C WRK0,WR00,WRK1,WRK6,WRK7,WRK8,WRK9,WTRAV = SEGMENTS DE TRAVAIL
  21. C NSTRSS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  22. C SIG0(NSTRSS) = CONTR. AU DEBUT DU PAS D'INTEGRATION
  23. C DEPST(NSTRSS) = INCREMENT DES DEFORM. ELASTIQUES
  24. C NVARI = NBR. DE VARIABLES INTERNES
  25. C VAR0(NVARI) = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  26. C CE TABLEAU CONTIENT DANS L'ORDRE:
  27. C - LA DEFORM. PLASTIQUE CUMULEE
  28. C - LES VARIABLES INTERNES PILOTANT LES EQUATIONS DU
  29. C MODELE
  30. C - LES DEFORM. PLASTIQUES
  31. C MFR1 = INDICE DE LA FORMULATION MECANIQUE; SEULEMENT
  32. C MASSIF OU COQUE POUR LES MATERIAUX ENDOMMAGEABLES
  33. * CMATE = NOM DU MATERIAU
  34. * N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE
  35. * N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE
  36. * MFR1 = NUMERO DE LA FORMULATION
  37. * IFOURB= OPTION DE CALCUL
  38. * IB = NUMERO DE L ELEMENT COURANT
  39. * IGAU = NUMERO DU POINT COURANT
  40. * EPAIST= EPAISSEUR
  41. * NBPGAU= NBRE DE POINTS DE GAUSS
  42. * MELE = NUMERO DE L ELEMENT FINI
  43. * NPINT = NBRE DE POINTS D INTEGRATION
  44. * NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES
  45. * NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES
  46. * SECT = SECTION
  47. * LHOOK = TAILLE DE LA MATRICE DE HOOKE
  48. C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  49. C XCARB(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  50. C INPLAS = INDICE DU MODELE DECRIVANT LE MATERIAU ENDOMMAGEABLE
  51. C NYKX = DIMENSION DE YKX
  52. C YKX(NYKX) = TABLEAU CONTENANT LES COURBES DE TRACTION(EPS,SIGMA)
  53. C A DES TEMPERATURES DONNEES (Ti);IL EST DE LA FORME
  54. C SUIVANTE :
  55. C (T1,X1,F(X1,T1),X2,F(X2,T1),..,T2,Y1,F(Y1,T2),Y2,..)
  56. C NNKX = NBRE DE COURBES DE TRACTION
  57. C NKX = TABLEAU CONTENANT LE NBRE DE PTS POUR CHAQUE COURBE
  58. C DE TRACTION (EX:SI NKX(1)=4,LA COURBE DE TRACTION
  59. C A LA TEMPERATURE T1 SERA DEFINIE PAR 4 POINTS)
  60. C PRECIS = PRECISION DES ITERATIONS INTERNES
  61. C TO = TEMPERATURE INITIALE
  62. C T1 = TEMPERATURE FINALE
  63. C Tref = TEMPERATURE DE REFERENCE
  64. C NYOG = 2*NBRE DE POINTS DE LA COURBE T-->YOG(T)
  65. C YOG(NYOG) = TABLEAU CONTENANT LES VALEURS CONNUES DU MODULE
  66. C D'YOUNG E A DES TEMPERATURES PRECISES T
  67. C NYNU = 2*NBRE DE POINTS DE LA COURBE T-->YNU(T)
  68. C YNU(NYNU) = TABLEAU CONTENANT LES VALEURS CONNUES DU COEFF. DE
  69. C POISSON YNU A DES TEMPERATURES PRECISES T
  70. C NYALFA = 2*NBRE DE POINTS DE LA COURBE T-->YALFA(T)
  71. C YALFA(NYALFA) = TABLEAU CONTENANT LES VALEURS CONNUES DU COEFF. DE
  72. C DILAT. THERMIQUE YALFA A DES TEMPERATURES PRECISES T
  73. C NYN = 2*NBRE DE POINTS DE LA COURBE T-->YN(T)
  74. C YN(NYN) = TABLEAU CONTENANT LES VALEURS CONNUES DU SEUIL
  75. C D'ENDOMT YN A DES TEMPERATURES PRECISES T
  76. C NYM = 2*NBRE DE POINTS DE LA COURBE T-->YM(T)
  77. C YM(NYM) = TABLEAU CONTENANT LES VALEURS CONNUES DE L'ENDOMT
  78. C CRITIQUE YM A DES TEMPERATURES PRECISES T
  79. C NYKK = 2*NBRE DE POINTS DE LA COURBE T-->YKK(T)
  80. C YKK(NYKK) = TABLEAU CONTENANT LES VALEURS CONNUES DE LA DEFORM.
  81. C A RUPTURE YKK A DES TEMPERATURES PRECISES T
  82. C
  83. C SORTIES:
  84. C -------
  85. C NCOURB = DIMENSION DE TRUC
  86. C TRUC = TABLEAU DE DIMENSION LA PLUS GRANDE PARMI CELLES
  87. C DEFINIES DANS NKX;IL CONTIENDRA LA COURBE DE
  88. C TRACTION A LA TEMPERATURE T,OBTENUE GRACE AU SOUS-
  89. C PROGRAMME TRACTT.
  90. C SIGF(NSTRSS)= CONTR. A LA FIN DU PAS D'INTEGRATION
  91. C VARF(NVARI)= VARIABLES INTERNES A LA FIN DU PAS D'INTEGRATION
  92. C DEFP(NSTRSS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS
  93. C D'INTEGRATION
  94. C KERRE = INDICE QUI REGIT LES ERREURS
  95. C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (2IEME VAR.
  96. C INT.) EST EN DEHORS DE LA COURBE DE TRACTION, DS.
  97. C LE CAS DE L'ECROUISSAGE ET DE L'ENDOMM. ISOTROPES.
  98. C CECI PEUT SE PRODUIRE SUITE A L'APPEL A "CRIDAM"
  99. C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE
  100. C POUR LE MODELE CONSIDERE OU S'IL Y A INCOMPATIBILITE
  101. C ENTRE MFR1 ET IFOURB
  102. C ==================================================================
  103. C
  104. IMPLICIT INTEGER(I-N)
  105. IMPLICIT REAL*8 (A-H,O-Z)
  106. -INC CCOPTIO
  107. -INC DECHE
  108. SEGMENT IECOU
  109. * COMMON/IECOU/NYOG,NYNU,NYALFA,icow1,NYN,NYM,NYKK,
  110. INTEGER NYOG,NYNU,NYALFA,NYSMAX,NYN,NYM,NYKK,
  111. 1 icow8,icow9,icow10,icow11,NYRHO,icow13,NNKX,NYKX,icow16,
  112. * . NYALF1,NYBET1,NYR, NYA, NYRHO,NSIGY, NNKX, NYKX, IND,
  113. 2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  114. * . NSOM,NINV,NINCMA,NCOMP,JELEM,LEGAUS,INAT,NCXMAT,
  115. 3 icow25,icow26,icow27,icow28,icow29,icow30,icara,
  116. * . LTRAC,MFR,IELE,NHRM,NBNN,NBELEM,ICARA,
  117. 4 icow32,icow33, NSTRSS,MFR1,NBGMAT,NELMAT,icow38,
  118. * . LW2, NDEF, NSTRSS,MFR1,NBGMAT,NELMAT,MSOUPA,
  119. 5 icow39,icow40,icow41,NNVARI,icow43,icow44
  120. * . NUMAT1,LENDO,NBBB,NNVARI,KERR1,MELEME
  121. INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  122. . icow51,icow52,icow53,icow54,icow55,icow56
  123. . icow57,icow58
  124. ENDSEGMENT
  125.  
  126. *
  127. SEGMENT WRK6
  128. REAL*8 BB(NSTRS,NNVARI),R(NSTRS),XMU(NSTRS)
  129. REAL*8 S(NNVARI),QSI(NNVARI),DDR(NSTRS),BBS(NSTRS)
  130. REAL*8 SIGMA(NSTRS),SIGGD(NSTRS)
  131. ENDSEGMENT
  132. *
  133. SEGMENT WRK7
  134. REAL*8 F(NCOURB,2),W(NCOURB),TRUC(NCOURB)
  135. ENDSEGMENT
  136. *
  137. SEGMENT WRK8
  138. REAL*8 DD(NSTRS,NSTRS),DDV(NSTRS,NSTRS),DDINV(NSTRS,NSTRS)
  139. REAL*8 DDINVp(NSTRS1,NSTRS1)
  140. ENDSEGMENT
  141. *
  142. SEGMENT WRK9
  143. REAL*8 YOG(NYOG),YNU(NYNU),YALFA(NYALFA),YSMAX(NYSMAX)
  144. REAL*8 YN(NYN),YM(NYM),YKK(NYKK),YALFA1(NYALF1)
  145. REAL*8 YBETA1(NYBET1),YR(NYR),YA(NYA),YKX(NYKX),YRHO(NYRHO)
  146. REAL*8 SIGY(NSIGY)
  147. INTEGER NKX(NNKX)
  148. ENDSEGMENT
  149.  
  150. *
  151. DIMENSION DEPS(6)
  152. *
  153. C
  154. C ------------------------------------------------
  155. C DEFINITION DU NBR. MAXIMUM D'ITERATIONS INTERNES
  156. C ------------------------------------------------
  157. DATA ITMAX/30/
  158. c
  159. NXMATT = NMATT
  160. C
  161. ZEERO=0.D0
  162. PRECI=PRECIS
  163. PRECIG=-2.D0*PRECIS
  164. PRECG=ABS(PRECIG)
  165. IALPH=0
  166. T=T1
  167. YUNGV=0.D0
  168. XNUV=0.D0
  169. C
  170. C ================================================================
  171. C MAJ. POUR LE CAS ELASTIQUE ET TEST POUR SAVOIR SI ON A PLASTIFIE
  172. C ================================================================
  173. CALL ZDANUL(DEFP,NSTRSS)
  174. ALPH=1.D0
  175. ALPH0=0.D0
  176. ALPH1=1.D0
  177. X=VAR0(2)
  178.  
  179. CALL INITTE(NYOG,YOG,NYNU,YNU,NYRHO,YRHO,NYALFA,YALFA,NYN,
  180. 1 YN,NYM,YM,NYKK,YKK,NYKX,YKX,NNKX,NKX,X,T0,
  181. 2 XMAT1,NXMATT,TRUC,NCOURB)
  182. CALL INITTE(NYOG,YOG,NYNU,YNU,NYRHO,YRHO,NYALFA,YALFA,NYN,
  183. 1 YN,NYM,YM,NYKK,YKK,NYKX,YKX,NNKX,NKX,X,T1,
  184. 2 XMAT2,NXMATT,TRUC,NCOURB)
  185. ALFA0=XMAT1(4)
  186. ALFA1=XMAT2(4)
  187. CTEPS=ALFA0*(T0-TREF)-ALFA1*(T1-TREF)
  188. C
  189. C==================================================
  190. C CALCUL DE L INCREMENT DE CONTRAINTES
  191. C DEPST EST L INCREMENT DE DEFORMATIONS ELASTIQUES
  192. C==================================================
  193. C
  194. CALL ELAST1(1,IFOUR,VAR0,NNVARI,XMAT2,NXMATT,YUNGV,XNUV,
  195. 1 XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  196. C
  197. CALL ZDANUL(DSIGT,NSTRSS)
  198. DO 13 I=1,NSTRSS
  199. DO 13 J=1,NSTRSS
  200. DSIGT(I)=DSIGT(I)+DD(I,J)*DEPST(J)
  201. 13 CONTINUE
  202. C
  203. C===============================================
  204. C CALCUL DE L INCREMENT DE DEFORMATIONS TOTALES
  205. C===============================================
  206. C
  207. CALL ELAST1(2,IFOUR,VAR0,NNVARI,XMAT1,NXMATT,YUNGV,XNUV,
  208. 1 XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  209. CALL ELAST1(2,IFOUR,VAR0,NNVARI,XMAT2,NXMATT,YUNGV,XNUV,
  210. 1 XCARB,ICARA,MFR1,NSTRSS,DDINV,DDV,KERRE,INDIC,ITHHER)
  211. CALL ZDANUL(DEPS,NSTRSS)
  212. AA=1.D0
  213. DO 45 I=1,NSTRSS
  214. DEPS(I)=DEPST(I)
  215. IF (I.GT.3) AA=0.D0
  216. DO 46 J=1,NSTRSS
  217. DEPS(I)=DEPS(I)+(DDINV(I,J)*SIG0(J))
  218. DEPS(I)=DEPS(I)-(DD(I,J)*SIG0(J))
  219. 46 CONTINUE
  220. DEPS(I)=DEPS(I)-(AA*CTEPS)
  221. 45 CONTINUE
  222. C
  223. C ------
  224. C
  225. CALL INITTE(NYOG,YOG,NYNU,YNU,NYRHO,YRHO,NYALFA,YALFA,NYN,
  226. 1 YN,NYM,YM,NYKK,YKK,NYKX,YKX,NNKX,NKX,X,T,
  227. 2 XMAT,NXMATT,TRUC,NCOURB)
  228.  
  229. 5 DO 10 I=1,NSTRSS
  230. SIGF(I)=SIG0(I)+DSIGT(I)
  231. SIGMA(I)=SIG0(I)
  232. 10 CONTINUE
  233. DO 20 I=1,NVARI
  234. VARF(I)=VAR0(I)
  235. 20 CONTINUE
  236. C
  237. IF (ITHHER.EQ.2.AND.IALPH.GT.0) GOTO 32
  238. 30 CALL CRIDAM(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,INPLAS,
  239. 1NYKX,YKX,NNKX,NKX,NNVARI,MFR1,XF,KERRE,SIGGD,YY,T,TRUC,
  240. 2 NNFUS)
  241. IF (XF.LT.PRECI) GOTO 1000
  242. CALL CRIDOM(SIG0,NSTRSS,VAR0,XMAT,NXMATT,XCARB,ICARA,INPLAS,
  243. 1NNVARI,MFR1,XG,KERRE,SIGGD,Y)
  244. C
  245. C............ Y DESIGNE LE SEUIL D'ENDOMMAGEMENT .................
  246. IF (Y.EQ.ZEERO) THEN
  247. GOTO 32
  248. ENDIF
  249. TTEST=XG/Y
  250. IF (TTEST.GT.PRECIG) INDIC=1
  251. C
  252. C
  253. C ==============
  254. C ON A PLASTIFIE
  255. C ==============
  256. C ------------------------------------------------------------------
  257. C DEFINITION DE DD, BB, R, S, XMU, ET QSI:
  258. C . . .
  259. C SIG = DD*EPSELAS + BB*Q AVEC:
  260. C - SIG = CONTR.
  261. C - EPSELAS = DEFORM. ELASTIQUES
  262. C - DD = MATR. DE HOOK ENDOMMAGEE
  263. C - Q = VARIABLES INTERNES PILOTANT LES LOIS DU MODELE
  264. C R = LOI D'ECOULEMENT PLASTIQUE
  265. C S = LOIS D'EVOLUTION DES VARIABLES Q
  266. C XMU = DERIVEE DU CRITERE PAR RAPPORT AUX CONTR.
  267. C QSI = DERIVEE DU CRITERE PAR RAPPORT AUX VARIABLES Q
  268. C ------------------------------------------------------------------
  269. 32 ICO=0
  270. C
  271. C
  272. C
  273. C ------------------------------------------
  274. C DEBUT DE LA BOUCLE DES ITERATIONS INTERNES
  275. C --------------------------------------------------------------------
  276. C A CHAQUE ITERATION, CORRESPONDANT A SIGF ET VARF CALCULES AU RANG I,
  277. C "ENDOM1" CALCULE DD, BB, R, S, XMU ET QSI EN FONCTION DU MODELE DU
  278. C MATERIAU ENDOMMAGEABLE.
  279. C --------------------------------------------------------------------
  280. C
  281. 35 CALL TRACTT(YKX,NKX,TRUC,NYKX,NNKX,NNFUS,T)
  282. NFUST=2*NNFUS
  283. CALL DERTRA(NYOG,YOG,T,TYUNG,YUNGV,TINF,TSUP)
  284. CALL DERTRA(NYNU,YNU,T,TNU,XNUV,TINF,TSUP)
  285. CALL ENDOM1(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,
  286. 1NNFUS,TRUC,NNVARI,MFR1,IFOURB,DD,DDINV,
  287. 2BB,R,S,XMU,QSI,KERRE,SIGGD,W,F,INDIC,ITHHER,
  288. 3YUNGV,XNUV)
  289. IF (KERRE.NE.0) GOTO 1000
  290. C
  291. CALL ZDANUL(DDR,NSTRSS)
  292. DO 50 J=1,NSTRSS
  293. DO 40 K=1,NSTRSS
  294. DDR(J)=DDR(J)+DD(J,K)*R(K)
  295. 40 CONTINUE
  296. 50 CONTINUE
  297. C
  298. CALL ZDANUL(BBS,NSTRSS)
  299. DO 70 J=1,NSTRSS
  300. DO 60 K=1,NNVARI
  301. BBS(J)=BBS(J)+BB(J,K)*S(K)
  302. 60 CONTINUE
  303. 70 CONTINUE
  304. C
  305. DLAMA1=0.D0
  306. DLAMA2=0.D0
  307. DO 80 J=1,NSTRSS
  308. DLAMA1=DLAMA1+XMU(J)*(DDR(J)-BBS(J))
  309. 80 CONTINUE
  310. DO 90 J=1,NNVARI
  311. DLAMA2=DLAMA2+QSI(J)*S(J)
  312. 90 CONTINUE
  313. C
  314. C -------------------------------------------------------------
  315. C CALCUL DE LA VARIATION DU MULT. PLAST. PAR LE CRIT. DE PLAST.
  316. C -------------------------------------------------------------
  317. DT=T-T0
  318. DLAMDA=XF/(DLAMA1-DLAMA2)
  319. C
  320. C ------------------------------------------------------------------
  321. C CALCUL DES CONTR., DES VAR. INT. Q ET DE LA DEFORM. PLAST. CUMULEE
  322. C AU RANG I+1
  323. C ------------------------------------------------------------------
  324. DO 100 J=1,NSTRSS
  325. SIGF(J)=SIGF(J)-DLAMDA*(DDR(J)-BBS(J))
  326. 100 CONTINUE
  327. VARF(2)=VARF(2)+DLAMDA*S(1)
  328. C
  329. IF (INDIC.EQ.1) THEN
  330. S2=S(2)
  331. IF (IFOURB.NE.-2) THEN
  332. XA=1.D0+VARF(3)
  333. CALL ENDOM1(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,
  334. 1NNFUS,TRUC,NNVARI,MFR1,IFOURB,DD,DDINV,
  335. 2BB,R,S,XMU,QSI,KERRE,SIGGD,W,F,INDIC,ITHHER,
  336. 3YUNGV,XNUV)
  337.  
  338. IF (KERRE.NE.0) GOTO 1000
  339. XB=VARF(3)+DLAMDA*S(2)*(1.D0-VARF(3))
  340. XDELTA=(XA*XA)-4.D0*XB
  341. IF (XDELTA.LT.0.D0) THEN
  342. VARF(3)=VARF(3)+DLAMDA*S2
  343. ELSE
  344. XDELTA=SQRT(XDELTA)
  345. VARF(3)=(XA-XDELTA)/2.D0
  346. ENDIF
  347. ELSE
  348. C
  349. C CAS DES CONTR. PLANES - MAJ DE D
  350. C
  351. XA=VARF(3)+0.5D0*DLAMDA*S(2)
  352. CALL ENDOM1(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,
  353. 1NNFUS,TRUC,NNVARI,MFR1,IFOURB,DD,DDINV,
  354. 2BB,R,S,XMU,QSI,KERRE,SIGGD,W,F,INDIC,ITHHER,
  355. 3YUNGV,XNUV)
  356. XB=0.5D0*DLAMDA*S(2)*(1.D0-VARF(3))
  357. XDELTA=(1.D0-XA)*(1.D0-XA)-4.D0*XB
  358. IF (XDELTA.LT.0.D0) THEN
  359. VARF(3)=VARF(3)+DLAMDA*S2
  360. ELSE
  361. XDELTA=SQRT(XDELTA)
  362. VARF(3)=(1.D0+XA-XDELTA)/2.D0
  363. ENDIF
  364. ENDIF
  365. ENDIF
  366. C
  367. VARF(1)=VARF(1)+DLAMDA/(1.D0-VARF(3))
  368. C
  369. C CONTR. PLANES - MAJ DES DEFORM. PLAST.
  370. C
  371. IF (IFOURB.EQ.-2) THEN
  372. CALL ENDOM1(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,
  373. 1 NNFUS,TRUC,NNVARI,MFR1,IFOURB,DD,DDINV,
  374. 2 BB,R,S,XMU,QSI,KERRE,SIGGD,W,F,INDIC,ITHHER,
  375. 3 YUNGV,XNUV)
  376. DO 110 J=1,NSTRSS
  377. VARF(1+NNVARI+J)=VARF(1+NNVARI+J)+DLAMDA*R(J)
  378. 110 CONTINUE
  379. ENDIF
  380. C
  381. C
  382. C
  383. C -------------------------------------------------------------
  384. C SI LE CRIT. PLAST. EST ATTEINT, ON ARRETE LES ITERATIONS
  385. C SI LE NBR. MAXIMUM D'ITERATIONS EST DEPASSE, IDEM MAIS ERREUR
  386. C -------------------------------------------------------------
  387. 115 CALL CRIDAM(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,INPLAS,
  388. 1NYKX,YKX,NNKX,NKX,NNVARI,MFR1,XF,KERRE,SIGGD,YY,T,TRUC,
  389. 2 NNFUS)
  390. TEST=ABS(XF/YY)
  391. IF (KERRE.NE.0) GOTO 1000
  392. IF (TEST.LT.PRECI.AND.INDIC.EQ.1) GOTO 120
  393. IF (TEST.LT.PRECI.AND.INDIC.EQ.0) THEN
  394. CALL CRIDOM(SIGF,NSTRSS,VARF,XMAT,NXMATT,XCARB,ICARA,INPLAS,
  395. 1NNVARI,MFR1,XG,KERRE,SIGGD,Y)
  396. TTEST=XG/Y
  397. ATTEST=ABS(TTEST)
  398. IF (ATTEST.LT.PRECG) GOTO 120
  399. IF (TTEST.LE.PRECIG) THEN
  400. IF (ALPH.EQ.1.D0) GOTO 120
  401. A=1.D0
  402. ELSE
  403. A=-1.D0
  404. ENDIF
  405. ALPH=ALPH1+0.5D0*A*ABS(ALPH1-ALPH0)
  406. ALPH0=ALPH1
  407. ALPH1=ALPH
  408. IALPH=IALPH+1
  409. T=(1-ALPH)*T0+ALPH*T1
  410. X=VARF(2)
  411. CALL INITTE(NYOG,YOG,NYNU,YNU,NYRHO,YRHO,NYALFA,YALFA,
  412. 1 NYN,YN,NYM,YM,NYKK,YKK,NYKX,YKX,NNKX,NKX,X,T,
  413. 2 XMAT,NXMATT,TRUC,NCOURB)
  414. ALFFA=XMAT(4)
  415. C
  416. CALL ELAST1(1,IFOURB,VAR0,NNVARI,XMAT,NXMATT,YUNGV,
  417. 1 XNUV,XCARB,ICARA,MFR1,NSTRSS,DDINV,DDV,KERRE,INDIC,ITHHER)
  418. CALL ELAST1(2,IFOURB,VAR0,NNVARI,XMAT1,NXMATT,YUNGV,
  419. 1 XNUV,XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  420. CALL MULMAT(DDV,DDINV,DD,NSTRSS,NSTRSS,NSTRSS)
  421. C
  422. CTEPS=ALFA0*(T0-TREF)-ALFFA*(T-TREF)
  423. CALL ZDANUL(DSIGT,NSTRSS)
  424. C
  425. DO 8 I=1,NSTRSS
  426. AA=1.D0
  427. DO 9 J=1,NSTRSS
  428. IF (J.GT.3) AA=0.D0
  429. DSIGT(I)=DSIGT(I)+DDINV(I,J)*((ALPH*DEPS(J))+(AA*CTEPS))
  430. DSIGT(I)=DSIGT(I)+DDV(I,J)*SIG0(J)
  431. 9 CONTINUE
  432. DSIGT(I)=DSIGT(I)-SIG0(I)
  433. 8 CONTINUE
  434. IF (IALPH.GT.50) GOTO 1000
  435. GOTO 5
  436. ENDIF
  437. ICO=ICO+1
  438. IF (ICO.GT.ITMAX) THEN
  439. KERRE=2
  440. GOTO 1000
  441. ENDIF
  442. GOTO 35
  443. C ----------------------------------------
  444. C FIN DE LA BOUCLE DES ITERATIONS INTERNES
  445. C ----------------------------------------
  446. C
  447. C
  448. C ====================================================
  449. C FIN DES ITERATIONS INTERNES - LE CRITERE EST VERIFIE
  450. C MAJ DES VARIABLES
  451. C ====================================================
  452. C
  453. C ----------------------------------------------------------------
  454. C CALCUL DES DEFORM. ELAST. A LA FIN DU PAS - ELLES SERONT DS. DDR
  455. C ----------------------------------------------------------------
  456. 120 CALL ELAST1(2,IFOURB,VARF,NNVARI,XMAT,NXMATT,YUNGV,XNUV,
  457. 1XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  458. CALL ZDANUL(DDR,NSTRSS)
  459. DO 140 J=1,NSTRSS
  460. DO 130 K=1,NSTRSS
  461. DDR(J)=DDR(J)+DD(J,K)*SIGF(K)
  462. 130 CONTINUE
  463. 140 CONTINUE
  464. C
  465. C -----------------------------------------
  466. C CALCUL DES DEFORM. PLAST. A LA FIN DU PAS
  467. C ----------------------------------------------------------------
  468. C ON CALCULE D'ABORD, DS. BBS, LES DEFORM. TOTALES A LA FIN DU PAS
  469. C ----------------------------------------------------------------
  470. CALL ELAST1(2,IFOURB,VAR0,NNVARI,XMAT1,NXMATT,YUNGV,XNUV,
  471. 1XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  472. CALL ZDANUL(BBS,NSTRSS)
  473. DO 170 J=1,NSTRSS
  474. DO 160 K=1,NSTRSS
  475. BBS(J)=BBS(J)+DD(J,K)*SIGMA(K)
  476. 160 CONTINUE
  477. 170 CONTINUE
  478. AA=1.D0
  479. DO 180 J=1,NSTRSS
  480. IF (J.GT.3) AA=0.D0
  481. BBS(J)=BBS(J)+VAR0(1+NNVARI+J)+(ALPH*DEPS(J))+(AA*CTEPS)
  482. 180 CONTINUE
  483. C
  484. C
  485. DO 190 J=1,NSTRSS
  486. IF (IFOURB.EQ.-2) GOTO 195
  487. VARF(1+NNVARI+J)=BBS(J)-DDR(J)
  488. 195 DEFP(J)=VARF(1+NNVARI+J)-VAR0(1+NNVARI+J)
  489. 190 CONTINUE
  490. IF (ALPH.EQ.1.D0) GOTO 900
  491. IF (INDIC.EQ.0) THEN
  492. ALPH=1.D0-ALPH
  493. DO 520 J=1,NVARI
  494. VAR0(J)=VARF(J)
  495. 520 CONTINUE
  496. DO 525 J=1,NSTRSS
  497. SIGMA(J)=SIGF(J)
  498. 525 CONTINUE
  499. C
  500. T0=T
  501. T=T1
  502. DO 530 J=1,NXMATT
  503. XMAT1(J)=XMAT(J)
  504. XMAT(J)=XMAT2(J)
  505. 530 CONTINUE
  506. ALFA0=XMAT1(4)
  507. ALFFA=XMAT(4)
  508. CALL ELAST1(1,IFOURB,VAR0,NNVARI,XMAT,NXMATT,YUNGV,
  509. 1 XNUV,XCARB,ICARA,MFR1,NSTRSS,DD,DDV,KERRE,INDIC,ITHHER)
  510. CALL ELAST1(2,IFOURB,VAR0,NNVARI,XMAT1,NXMATT,YUNGV,
  511. 1 XNUV,XCARB,ICARA,MFR1,NSTRSS,DDINV,DDV,KERRE,INDIC,
  512. 2 ITHHER)
  513. CALL MULMAT(DDV,DD,DDINV,NSTRSS,NSTRSS,NSTRSS)
  514. C
  515. CTEPS=ALFA0*(T0-TREF)-ALFFA*(T-TREF)
  516. CALL ZDANUL(DSIGT,NSTRSS)
  517. C
  518. DO 515 I=1,NSTRSS
  519. AA=1.D0
  520. DO 510 J=1,NSTRSS
  521. IF (J.GT.3) AA=0.D0
  522. DSIGT(I)=DSIGT(I)+DD(I,J)*((ALPH*DEPS(J))+(AA*CTEPS))
  523. DSIGT(I)=DSIGT(I)+DDV(I,J)*SIGMA(J)
  524. 510 CONTINUE
  525. DSIGT(I)=DSIGT(I)-SIGMA(I)
  526. SIGF(I)=SIGF(I)+DSIGT(I)
  527. 515 CONTINUE
  528. GOTO 30
  529. ENDIF
  530. C
  531. C ==================================================================
  532. C SI LE PT. DE GAUSS EST ROMPU, ON ANNULE PROGRESSIVEMENT SES CONTR.
  533. C ==================================================================
  534. 900 DCT=XMAT(7)
  535. DCC=1.1D0*DCT
  536. XD=VARF(3)
  537. IF (XD.GE.DCC) THEN
  538. VARF(3)=1.D0
  539. CALL ZDANUL(SIGF,NSTRSS)
  540. ENDIF
  541. IF (XD.GT.DCT.AND.XD.LT.DCC) THEN
  542. VARF(3)=1.D0
  543. BETA=11.D0-10.D0*(XD/DCT)
  544. DO 950 J=1,NSTRSS
  545. 950 SIGF(J)=SIGF(J)*BETA
  546.  
  547. ENDIF
  548. 1000 RETURN
  549. END
  550. *******************************************************************
  551.  
  552.  
  553.  
  554.  
  555.  
  556.  
  557.  
  558.  
  559.  
  560.  
  561.  
  562.  
  563.  
  564.  
  565.  
  566.  
  567.  
  568.  
  569.  
  570.  

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