Télécharger ccerac.eso

Retour à la liste

Numérotation des lignes :

ccerac
  1. C CCERAC SOURCE PV 17/12/08 21:15:14 9660
  2. SUBROUTINE CCERAC(wrk52,wrk53,wrk54,NVARI,
  3. 1 NSSINC,INV,IFOURB,IB,IGAU,NBPGAU,iecou,xecou)
  4. C CERACA SOURCE AM 00/12/13 21:15:32 4045
  5. c SUBROUTINE CERACA(WRK0,WR00,WRK1,WRK5,WRK7,WRK8,WRK9,WTRAV,
  6. c 1 INPLAS,MFR,DT,NSTRS,NVARI,NCOMAT,PRECIS,MSOUPA,JECHER,DTT,
  7. c 2 NSSINC,INV,KERRE,ICARA,IFOURB,NYOG,NYNU,NYALFA,NYSMAX,
  8. c 3 NYN,NYM,NYKK,NYALF1,NYBET1,NYR,NYA,NYKX,NNKX,NYRHO,NSIGY,T0,TF,
  9. c 5 TREF,TLIFE,ITHHER,NCOURB,CMATE,N2EL,N2PTEL,IB,IGAU,EPAIST,
  10. c 7 NBPGAU,MELE,NPINT,NBGMAT,NELMAT,SECT,LHOOK,CRIGI,
  11. c 8 KERREU1)
  12. C
  13. C---------------------------------------------------------------------
  14. C Objet: Calculer au cours d'un pas de temps DT, l'evolution des
  15. C variables internes a l'aide d'un schema Runge-Kutta 1.2
  16. C ---------------------------------------------------------------------
  17. C MFR1 <- MFR, XCARB <- XCAR, NSTRS1 <- NSTRS,
  18. C
  19. C---------------------------------------------------------------------
  20. C Entree: INPLAS type de materiau
  21. C MFR1 indice de la formulation mecanique(seulement massif ou coque
  22. C pour les materiaux endommageables)
  23. C DEPST(NSTRS1) increment des deformations totales
  24. C SIG0(NSTRS1) contraintes initiales
  25. C EPIN0(NSTRS1) deformations viscoplastiques initiales
  26. C VAR0(NVARI) variables internes initiales
  27. C NVARI nombre de variables internes
  28. C SIGY(NSIGY) courbe de la limite elastique en fonction de T°C
  29. C XMAT(NCOMAT) materiau
  30. C XCARB(ICARA) caracteristiques geometriques
  31. C YSMAX(NYSMAX) intervient ds. le test de convergence des iter.
  32. C PRECIS precision relative sur SIGMA
  33. C MSOUPA nombre maximal de sous pas autorises
  34. C JECHER = 0 avancer
  35. C = 1 rechercher sortie avec DTT
  36. C IFOURB = -2 EN CONTR.PLANES
  37. C -1 EN DEFORM. PLANES
  38. C 0 EN AXISYMETRIE
  39. C 1 EN SERIE DE FOURIER
  40. C 2 EN TRIDIM
  41. C CMATE = NOM DU MATERIAU
  42. C VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU
  43. C VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES
  44. C N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE
  45. C N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE
  46. C IB = NUMERO DE L ELEMENT COURANT
  47. C IGAU = NUMERO DU POINT COURANT
  48. C EPAIST= EPAISSEUR
  49. C NBPGAU= NBRE DE POINTS DE GAUSS
  50. C MELE = NUMERO DE L ELEMENT FINI
  51. C NPINT = NBRE DE POINTS D INTEGRATION
  52. C SECT = SECTION
  53. C LHOOK = TAILLE DE LA MATRICE DE HOOKE
  54. C DT pas de temps
  55. C ITHHER = 0 pas de chargement thermique et materiau constant
  56. C = 1 chargement thermique et materiau constant
  57. C = 2 chargement thermique et materiau(T)
  58. C-----------------------------------------------------------------------
  59. C
  60. C-----------------------------------------------------------------------
  61. C Sortie: SIGF(NSTRS) contraintes finales
  62. C EPINF(NSTRS) deformations viscoplastiques finales
  63. C VARF(NVARI) variables internes finales
  64. C DTT sous-increment de temps optimal (si JECHER=1)
  65. C NSSINC nombre de sous-increments si JECHER=0
  66. C INV = 1 si inversion
  67. C 0 sinon
  68. C KERRE = 0 si tout OK
  69. C <> 0 si entrees incoherentes
  70. C-----------------------------------------------------------------------
  71. C
  72. IMPLICIT INTEGER(I-N)
  73. IMPLICIT REAL*8(A-H,O-Z)
  74.  
  75. -INC PPARAM
  76. -INC CCOPTIO
  77. -INC DECHE
  78.  
  79. SEGMENT IECOU
  80. * COMMON/IECOU/NYOG,NYNU,NYALFA,NYSMAX,NYN,NYM,NYKK,
  81. INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  82. C INTEGER NYOG, NYNU, NYALFA,NYSMAX,NYN, NYM, NYKK,
  83. 1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  84. C . NYALF1,NYBET1,NYR, NYA, NYRHO,NSIGY, NNKX, NYKX, IND,
  85. 2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  86. C . NSOM, NINV, NINCMA,NCOMP, JELEM, LEGAUS,INAT, NCXMAT,
  87. 3 icow25,icow26,icow27,icow28,icow29,icow30,ICARA,
  88. C . LTRAC, MFR, IELE, NHRM, NBNN, NBELEM,ICARA,
  89. 4 icow32,icow33,NSTRS1,MFR1,NBGMAT,NELMAT,MSOUPA,
  90. C . LW2, NDEF, NSTRSS,MFR1, NBGMAT,NELMAT,MSOUPA,
  91. 5 icow39,icow40,icow41,icow42,icow43,icow44
  92. C . NUMAT1,LENDO, NBBB, NNVARI,KERR1, MELEME
  93. INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  94. . icow51,icow52,icow53,icow54,icow55,icow56
  95. . icow57,icow58
  96. ENDSEGMENT
  97.  
  98. SEGMENT XECOU
  99. * COMMON/XECOU/DTOPTI,TSOM,TCAR,DTT,DT,TREFA,TEMP00
  100. REAL*8 xcow1, xcow2,xcow3,DTT,DT,xcow6, xcow7
  101. C REAL*8 DTOPTI,TSOM, TCAR, DTT, DT, TREFA,TEMP00
  102. ENDSEGMENT
  103.  
  104.  
  105.  
  106. *
  107. *
  108. DIMENSION SIG(8),EPSV(8),VAR(100)
  109. DIMENSION DSPT(8),XX(8),SIG1(8),EPSV1(8)
  110. DIMENSION VAR1(100),EVP1(8),VARP1(100)
  111. DIMENSION EVP2(8),VARP2(100)
  112. DIMENSION CRIGI(12)
  113. C**** debut ajout Eloi
  114. logical dtlibr
  115. C**** fin ajout Eloi
  116. NCOMAT = NMATT
  117. C**** debut ajout Eloi
  118. dtlibr=.TRUE.
  119. C**** fin ajout Eloi
  120. C ==================================================
  121. C RECHERCHE DU NUMERO DE LA VARIABLE INTERNE EGALE A 1
  122. C SI ON A ENDOMMAGEMENT GENERALISE
  123. C
  124. IBID2 = NVARI-1
  125. C*************************************
  126. C On cherche le numero de la propriété du matériau correspondant à ENDG
  127. C
  128. IF (MFR1.EQ.1.AND.IFOMOD.EQ.2) THEN
  129. IBID1 = 20
  130. ELSE
  131. IBID1 = 15
  132. ENDIF
  133. C Test pour savoir si on a dépassé la limite de déformation de fluage
  134. C
  135. IF (VAR0(IBID2).EQ.1.D0) THEN
  136. DO 10 IJ=1,NVARI
  137. VARF(IJ) = VAR0(IJ)
  138. 10 CONTINUE
  139. DO 20 K=1,NSTRS1
  140. EPINF(K) = EPIN0(K)
  141. SIGF(K) = 0.D0
  142. 20 CONTINUE
  143. GOTO 1000
  144. ENDIF
  145. DO 30 I=1,NSTRS1
  146. TOTO = ABS(EPIN0(I))
  147. IF (TOTO.GE.XMAT(IBID1)) THEN
  148. DO 40 K=1,NSTRS1
  149. EPINF(K) = EPIN0(K)
  150. SIGF(K) = 0.D0
  151. 40 CONTINUE
  152. DO 50 IJ=1,(NVARI-1)
  153. VARF(IJ) = VAR0(IJ)
  154. 50 CONTINUE
  155. VARF(IBID2) = 1.D0
  156. GOTO 1000
  157. ENDIF
  158. 30 CONTINUE
  159. C*************************************
  160.  
  161.  
  162. C =========================================================
  163. KERRE = 0
  164. IF (MFR1.NE.1.AND.MFR1.NE.3.OR.IFOURB.EQ.1) THEN
  165. KERRE = 99
  166. RETURN
  167. ENDIF
  168.  
  169. IF (MFR1.EQ.3) THEN
  170. THICK = XCARB(1)
  171. ALFA = XCARB(2)
  172. ENDIF
  173. BORNE = 2.0
  174. RMAX = 1.3
  175. RMIN = 0.7
  176. DIV = 7.0
  177. FAC = 3.0
  178. C
  179. C CALCUL DES INCREMENTS DE CONTRAINTES ELASTIQUES
  180. C
  181. CALL CALSIG(DEPST,DDAUX,NSTRS1,CMATE,VALMAT,VALCAR,N2EL,N2PTEL,
  182. 1 MFR1,IFOURB,IB,IGAU,EPAIST,NBPGAU,MELE,NPINT,NBGMAT,
  183. 2 NELMAT,SECT,LHOOK,TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,
  184. 3 CRIGI,DSIGT,IRTD)
  185. C PRINT *,'DSIGT =',DSIGT
  186. *
  187. IF (IRTD.NE.1) THEN
  188. KERRE=69
  189. GOTO 1000
  190. ENDIF
  191. *
  192. IF (MFR1.EQ.3) THEN
  193. DO 60 I=1,NSTRS1/2
  194. SIG0( I) = SIG0( I)/THICK
  195. SIG0(NSTRS1/2+I) = SIG0(NSTRS1/2+I)*6.0D0/THICK/THICK
  196. DSIGT( I) = DSIGT( I)/THICK
  197. DSIGT(NSTRS1/2+I)= DSIGT(NSTRS1/2+I)*6.0D0/THICK/THICK
  198. 60 CONTINUE
  199. IF (IFOURB.EQ.-2) THEN
  200. SIG0 (2) = 0.D0
  201. SIG0 (4) = 0.D0
  202. DSIGT(2) = 0.D0
  203. DSIGT(4) = 0.D0
  204. ENDIF
  205. ENDIF
  206. C
  207. C------------------------------------------
  208. C CONTROLE DE LA COHERENCE DES ENTREES
  209. C Supprimée dans ce cas
  210. C------------------------------------------
  211. IF (DT.LE.0.D0) KERRE = 414
  212. MOTERR(1:8) = ' CONST '
  213. C IF (INCOMAT.NE. 15)KERRE = 146
  214. C IF (IFOURB.EQ.-2.AND.NCOMAT.NE. 9)KERRE = 146
  215. IF (IFOURB.NE.2.AND.IFOURB.NE.-2.AND.IFOURB.NE.0
  216. & .AND.IFOURB.NE.-1) THEN
  217. KERRE = 194
  218. MOTERR(1:8) = 'FLUAGE'
  219. ENDIF
  220. IF (MFR1.EQ.1.AND.IFOMOD.EQ.2) THEN
  221. IBID = 15
  222. ELSE
  223. IBID = 10
  224. ENDIF
  225. XMAX=XMAT(IBID)
  226. *
  227. * TEST SUR XMAX MILL 8/3/91
  228. *
  229. IF (XMAX.EQ.0.D0) THEN
  230. XMAX=XMAT(1)*1.D-3
  231. ENDIF
  232. C
  233. C-----------------------------
  234. IF (KERRE.NE.0) THEN
  235. GOTO 1100
  236. ENDIF
  237. C
  238. C===========================================================
  239. C A PARTIR DE MAINTENANT, LES DEFORMATIONS
  240. C DE CISAILLEMENT NE SONT PLUS
  241. C DEFINIES PAR DES GAMA.
  242. C ON DIVISE DONC LES TERMES DE CISAILLEMENT PAR 2.
  243. C SEULES LES FORMULATIONS SUIVANTES SONT ACCEPTEES PAR CERACARO:
  244. C MFR1=1 (MASSIF)
  245. C MFR1=3 (COQUES MINCES)
  246. C
  247. C
  248. C Cas de la formulation massive
  249. C Les termes de cisaillement apparaissent
  250. C au delà de la troisieme composante
  251. C
  252. IF (MFR1.EQ.1) THEN
  253. DO 70 I=1,NSTRS1
  254. A=1.D0
  255. IF (I.GT.3) A=2.D0
  256. EPIN0(I)=EPIN0(I)/A
  257. 70 CONTINUE
  258. C
  259. C Cas des coques minces
  260. C Les termes de cisaillement apparaissent
  261. C pour la troisieme et la sixieme composante
  262. C uniquement dans les cas de calculs
  263. C tridimensionnels ou d'analyse de Fourier
  264. C
  265. ELSE IF (MFR1.EQ.3) THEN
  266. IF (IFOURB.EQ.2) THEN
  267. DO 80 I=1,NSTRS1
  268. A=1.D0
  269. IF (I.EQ.3) A=2.D0
  270. IF (I.EQ.6) A=2.D0
  271. EPIN0(I)=EPIN0(I)/A
  272. 80 CONTINUE
  273. ENDIF
  274. ENDIF
  275. C
  276. C===========================================================
  277. C
  278. C ----------------
  279. C INITIALISATION
  280. C ----------------
  281. C**** debut ajout Eloi
  282. ITERO = 0
  283. 1200 CONTINUE
  284.  
  285. itero = 1 + itero
  286. if (itero.ne.1) THEN
  287. c write(6,*) 'itero ib igau', itero,ib,igau
  288. dtlibr = .true.
  289. precis = precis * 7.
  290. c write(6,*) ' precision modifiée ', precis
  291. if (itero.gt.3) then
  292. **** kerre = 460
  293. kerre = 268
  294. return
  295. endif
  296. endif
  297. C**** fin ajout Eloi
  298. DTLEFT = DT
  299. TAU = DTLEFT
  300. ASIG = SQRT(PROCON(SIG0,SIG0,NSTRS1))
  301. ERRABS = PRECIS*ASIG
  302. IF (XMAX.GT.ASIG) ERRABS = PRECIS*XMAX
  303. DO 90 I=1,NSTRS1
  304. SIG(I) = SIG0(I)
  305. EPSV(I) = EPIN0(I)
  306. DSPT(I) = DSIGT(I)/DT
  307. 90 CONTINUE
  308. C
  309. DO 100 I=1,NVARI
  310. VAR(I)=VAR0(I)
  311. 100 CONTINUE
  312. C
  313. C ---------------------------------------------------------------------
  314. INV = 0
  315. NSSINC = 0
  316. C ---------------------------------------------------------------------
  317. C DEBUT DES ITERATIONS EN SSINCREMENTS /FIN SI DTLEFT = 0
  318. C ---------------------------------------------------------------------
  319. 1400 CONTINUE
  320. NSSINC = NSSINC + 1
  321. IF (NSSINC.GT.MSOUPA) THEN
  322. C**** debut ajout Eloi
  323. DTLIBR=.FALSE.
  324. GOTO 1200
  325. C**** fin ajout Eloi
  326. C*** KERRE = 460
  327. C*** GOTO 1100
  328. ENDIF
  329. C
  330. C---------------------------------------------------------------------
  331. C START OF CALCULATIONS
  332. C_____________________________________________________________________
  333.  
  334. IF(MFR1.NE.3) THEN
  335. CALL INCR11(TAU,SIG,EPSV,VAR,EVP1,VARP1,XMAT,NSTRS1,NVARI,MFR1)
  336. ELSE IF(MFR1.EQ.3) THEN
  337. CALL INCR33(TAU,SIG,EPSV,VAR,XMAT,EVP1,VARP1,ALFA,NSTRS1,NVARI)
  338. ENDIF
  339. C
  340. NITERA = 0
  341.  
  342. C --------------------------------------------------------------------
  343. C DEBUT DES ITERATIONS SUR TAU OPTIMAL /FIN SI RA PETIT
  344. C --------------------------------------------------------------------
  345.  
  346. 1300 CONTINUE
  347. NITERA = NITERA + 1
  348. IF (MFR1.NE.3) THEN
  349.  
  350. CALL ADVA11(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT,EVP1,VARP1,
  351. & XMAT,NSTRS1,NVARI,MFR1)
  352. CALL INCR11(TAU,SIG1,EPSV1,VAR1,EVP2,VARP2,XMAT,NSTRS1,NVARI,
  353. & MFR1)
  354. DO 110 I=1,NSTRS1
  355. EVP2(I) = 0.5D0*(EVP1(I)+EVP2(I))
  356. 110 CONTINUE
  357. CCC Eloi : on peut effectuer la boucle suivante sur
  358. CCC les variables internes propres a Ottosen
  359. CCC qui n'ont pas ete modifiees
  360. DO 120 I=1,NVARI
  361. VARP2(I)= 0.5*(VARP1(I) + VARP2(I))
  362. 120 CONTINUE
  363. CALL ADVA11(TAU,SIG,EPSV,VAR,SIGF,EPINF,VARF,DSPT,EVP2,VARP2,
  364. & XMAT,NSTRS1,NVARI,MFR1)
  365. ELSE
  366. C----------------------------------------------------------------------
  367. C CALCULATIONS FOR GENERALISED STRESS/STRAIN FORMULATIONS
  368. C----------------------------------------------------------------------
  369. CALL ADVA22(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT,EVP1,VARP1,
  370. & XMAT,NSTRS1,NVARI)
  371. CALL INCR33(TAU,SIG1,EPSV1,VAR1,XMAT,EVP2,VARP2,ALFA,NSTRS1,
  372. & NVARI)
  373. DO 140 I=1,NSTRS1,1
  374. EVP2(I) = 0.5D0*(EVP1(I)+EVP2(I))
  375. 140 CONTINUE
  376. CCC Eloi : on peut effectuer la boucle suivante sur
  377. CCC les variables internes propres a Ottosen
  378. CCC qui n'ont pas ete modifiees
  379. DO 150 I=1,NVARI
  380. VARP2(I) = 0.5D0*(VARP1(I)+VARP2(I))
  381. 150 CONTINUE
  382. CALL ADVA22(TAU,SIG,EPSV,VAR,SIGF,EPINF,VARF,DSPT,EVP2,VARP2,
  383. & XMAT,NSTRS1,NVARI)
  384. ENDIF
  385.  
  386. C ---------------------------------------------------------------------
  387. C CALCUL DU RAPPORT : ERREUR CALCULEE / ERREUR ADMISE
  388. C ---------------------------------------------------------------------
  389. DO 170 I=1,NSTRS1
  390. XX(I) = SIGF(I)-SIG1(I)
  391. 170 CONTINUE
  392. RA = SQRT(PROCON(XX,XX,NSTRS1))/(ERRABS)
  393. SQRA = SQRT(RA)
  394.  
  395. C ---------------------------------------------------------------------
  396. C TEST DE FIN D'ITERATIONS / MISE A JOUR DE TAU /OPTION JECHER
  397. C DIV =7 BORNE = 2
  398. C SI SQRA>7 TAU = TAU/7 ET NOUVEL ESSAI
  399. C SI 2<RA<7*7 ON VISE RA = 1 ET NOUVEL ESSAI
  400. C ------------------------------------------------------------------
  401. C**** debut ajout Eloi
  402. IF (dtlibr) Then
  403. C**** fin ajout Eloi
  404. IF (RA.GT.DIV*DIV ) THEN
  405. TAU = TAU/DIV
  406. GOTO 1300
  407. ELSEIF (RA.GT.BORNE) THEN
  408. TAU = TAU/SQRA
  409. GOTO 1300
  410. ENDIF
  411. C**** debut ajout Eloi
  412. endif
  413. C**** fin ajout Eloi
  414.  
  415. C ---------------------------------------------------------------------
  416. C ici ra < borne cas JECHER :
  417. C ---------------------------------------------------------------------
  418. IF (JECHER.EQ.1) THEN
  419. DTT = TAU
  420. NSSINC = NITERA
  421. IF ((NSSINC.EQ.1).AND.(RA.EQ.0.0)) GOTO 1100
  422. IF (NITERA.GE.8) GOTO 1100
  423. IF (FAC*SQRA.LT.1.0) THEN
  424. TAU = TAU*FAC
  425. GOTO 1300
  426. ELSEIF ((SQRA.LT.RMIN).OR.(SQRA.GT.RMAX)) THEN
  427. TAU = TAU/SQRA
  428. GOTO 1300
  429. ENDIF
  430. C ---------------------------------------------------------------------
  431. C ici rmin < sqra < rmax et nitera < 8
  432. C pas de mise @ jour des variables
  433. C ---------------------------------------------------------------------
  434. GOTO 1100
  435. ENDIF
  436. C ----------------------------------------------------------------------
  437. C FIN D'ITERATIONS / MISE A JOUR DES VARIABLES
  438. C ici RA < BORNE
  439. C fin des boucles sur tau optimal
  440. C on avance en temps
  441. C mise @ jour de SIG etc...
  442. C -------------------------------------------------------------------
  443. C --
  444. C -- ajout de ivtest = 1 par chloe
  445. C --
  446. IVTEST = 1
  447. INV = INV + IVTEST
  448. DO 180 I=1,NSTRS1
  449. SIG(I) = SIGF(I)
  450. EPSV(I) = EPINF(I)
  451. 180 CONTINUE
  452. DO 190 I=1,NVARI
  453. VAR(I) = VARF(I)
  454. 190 CONTINUE
  455. C*************************************
  456. C Test pour savoir si on a dépassé la limite de déformation totale
  457. C
  458. DO 200 I =1,NSTRS1
  459. TOTO = ABS(EPINF(I))
  460. IF(TOTO.GE.XMAT(IBID1)) THEN
  461. VARF(IBID2) = 1.D0
  462. DO 210 KI=1,NSTRS1,1
  463. SIGF(KI) = 0.D0
  464. 210 CONTINUE
  465. GOTO 1000
  466. ELSE
  467. VARF(IBID2) = 0.D0
  468. ENDIF
  469. 200 CONTINUE
  470. C*************************************
  471. C
  472. C --------------------------------------------------------------------
  473. C TEST DE FIN SS INCREMENTS / MISE A JOUR DE TAU
  474. C si SQRA<1/3 TAU = TAU*3
  475. C si 1/3<SQRA<RMIN on vise RA = 1
  476. C si RMIN<SQRA<RMAX TAU inchangé
  477. C si SQRA>RMAX on vise RA = 1
  478. C fin des boucles en ss increments si tau = dtleft
  479. C --------------------------------------------------------------------
  480. C
  481. IF (TAU.LT.DTLEFT) THEN
  482. DTLEFT = DTLEFT - TAU
  483. IF (FAC*SQRA.LT.1.D0) THEN
  484. TAU=TAU*FAC
  485. ELSEIF ( (SQRA.LT.RMIN).OR.(SQRA.GT.RMAX) ) THEN
  486. TAU=TAU/SQRA
  487. ENDIF
  488. IF (TAU.GT.DTLEFT) TAU = DTLEFT
  489. GOTO 1400
  490. ENDIF
  491. C
  492. IF (ABS(TAU-DTLEFT).GT.(TAU/1000.)) THEN
  493. WRITE ( IOIMP,* ) ' PROBLEME TAU > DTLEFT '
  494. KERRE = 223
  495. ENDIF
  496. C-----------------------------------------------------------------------
  497. 1100 CONTINUE
  498. IF (MFR1.EQ.3) THEN
  499. DO 220 I=1,NSTRS1/2
  500. SIGF( I) =SIGF( I)*THICK
  501. SIGF(NSTRS1/2+I) =SIGF(NSTRS1/2+ I)*THICK*THICK/6.0
  502. * DSIGT( I)=DSIGT( I)*THICK
  503. * DSIGT(NSTRS1/2+I)=DSIGT(NSTRS1/2+I)*THICK*THICK/6.0
  504. 220 CONTINUE
  505. ENDIF
  506.  
  507. C
  508. C===========================================================
  509. C RETOUR A LA DEFINITION NORMALE DES DEFORMATIONS
  510. C A SAVOIR: LES DEFORMATIONS DE CISAILLEMENT SONT
  511. C DEFINIES PAR DES GAMA.
  512. C SEULES LES FORMULATIONS SUIVANTES SONT ACCEPTEES PAR Ceracaro:
  513. C MFR1=1 (MASSIF)
  514. C MFR1=3 (COQUES MINCES)
  515. C
  516. C Cas de la formulation massive
  517. C Les termes de cisaillement apparaissent
  518. C au delà de la troisieme composante
  519. C
  520. IF (MFR1.EQ.1) THEN
  521. DO 230 I=1,NSTRS1
  522. A=1.D0
  523. IF (I.GT.3) A=2.D0
  524. EPIN0(I)=EPIN0(I)*A
  525. EPINF(I)=EPINF(I)*A
  526. 230 CONTINUE
  527. C
  528. C Cas des coques minces
  529. C Les termes de cisaillement apparaissent
  530. C pour la troisieme et la sixieme composante
  531. C uniquement dans les cas de calculs
  532. C tridimensionnels ou d'analyse de Fourier
  533. C
  534. ELSE IF (MFR1.EQ.3) THEN
  535. IF (IFOURB.EQ.2) THEN
  536. DO 240 I=1,NSTRS1
  537. A=1.D0
  538. IF (I.EQ.3) A=2.D0
  539. IF (I.EQ.6) A=2.D0
  540. EPIN0(I)=EPIN0(I)*A
  541. EPINF(I)=EPINF(I)*A
  542. 240 CONTINUE
  543. ENDIF
  544. ENDIF
  545. C
  546. C===========================================================
  547. C
  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.  

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