Télécharger endom.eso

Retour à la liste

Numérotation des lignes :

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

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