Télécharger ccerac.eso

Retour à la liste

Numérotation des lignes :

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

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